From a5ee52a4a2c2d0f4a9a3c75ce2f7fe260afa5be1 Mon Sep 17 00:00:00 2001 From: 8c82b2701c130ec4af8160ecac35027a <8c82b2701c130ec4af8160ecac35027a@app-learninglab.inria.fr> Date: Wed, 4 Jun 2025 16:39:15 +0000 Subject: [PATCH] corrected in conversion to numerical --- module3/exo1/analyse-syndrome-grippal.ipynb | 2217 ++++++++++++++++++- 1 file changed, 2184 insertions(+), 33 deletions(-) diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5..8f0699d 100644 --- a/module3/exo1/analyse-syndrome-grippal.ipynb +++ b/module3/exo1/analyse-syndrome-grippal.ipynb @@ -9,7 +9,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -28,10 +28,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 2, + "metadata": {}, "outputs": [], "source": [ "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" @@ -61,9 +59,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020252232083115052.026610.03122.040.0FRFrance
120252132406518883.029247.03628.044.0FRFrance
220252032026515814.024716.03023.037.0FRFrance
320251931626412394.020134.02418.030.0FRFrance
420251831811513975.022255.02721.033.0FRFrance
520251732215017291.027009.03326.040.0FRFrance
620251632856422550.034578.04334.052.0FRFrance
720251533572129592.041850.05344.062.0FRFrance
820251433757931232.043926.05647.065.0FRFrance
920251333967333686.045660.05950.068.0FRFrance
1020251235254345627.059459.07868.088.0FRFrance
1120251135946952154.066784.08978.0100.0FRFrance
1220251036033453048.067620.09079.0101.0FRFrance
1320250938453174994.094068.0126112.0140.0FRFrance
142025083136020124824.0147216.0203186.0220.0FRFrance
152025073208952195988.0221916.0312293.0331.0FRFrance
162025063273519258159.0288879.0408385.0431.0FRFrance
172025053334395318416.0350374.0499475.0523.0FRFrance
182025043350043332885.0367201.0522496.0548.0FRFrance
192025033252772238917.0266627.0377356.0398.0FRFrance
202025023257247242991.0271503.0384363.0405.0FRFrance
212025013231549214627.0248471.0345320.0370.0FRFrance
222024523201726185870.0217582.0302278.0326.0FRFrance
232024513201697187843.0215551.0302281.0323.0FRFrance
242024503136694126369.0147019.0205190.0220.0FRFrance
25202449310848799037.0117937.0163149.0177.0FRFrance
2620244838738178687.096075.0131118.0144.0FRFrance
2720244737628667626.084946.0114101.0127.0FRFrance
2820244635639949006.063792.08574.096.0FRFrance
2920244534734740843.053851.07161.081.0FRFrance
.................................
208819852132609619621.032571.04735.059.0FRFrance
208919852032789620885.034907.05138.064.0FRFrance
209019851934315432821.053487.07859.097.0FRFrance
209119851834055529935.051175.07455.093.0FRFrance
209219851733405324366.043740.06244.080.0FRFrance
209319851635036236451.064273.09166.0116.0FRFrance
209419851536388145538.082224.011683.0149.0FRFrance
20951985143134545114400.0154690.0244207.0281.0FRFrance
20961985133197206176080.0218332.0357319.0395.0FRFrance
20971985123245240223304.0267176.0445405.0485.0FRFrance
20981985113276205252399.0300011.0501458.0544.0FRFrance
20991985103353231326279.0380183.0640591.0689.0FRFrance
21001985093369895341109.0398681.0670618.0722.0FRFrance
21011985083389886359529.0420243.0707652.0762.0FRFrance
21021985073471852432599.0511105.0855784.0926.0FRFrance
21031985063565825518011.0613639.01026939.01113.0FRFrance
21041985053637302592795.0681809.011551074.01236.0FRFrance
21051985043424937390794.0459080.0770708.0832.0FRFrance
21061985033213901174689.0253113.0388317.0459.0FRFrance
210719850239758680949.0114223.0177147.0207.0FRFrance
210819850138548965918.0105060.0155120.0190.0FRFrance
210919845238483060602.0109058.0154110.0198.0FRFrance
2110198451310172680242.0123210.0185146.0224.0FRFrance
21111984503123680101401.0145959.0225184.0266.0FRFrance
2112198449310107381684.0120462.0184149.0219.0FRFrance
211319844837862060634.096606.0143110.0176.0FRFrance
211419844737202954274.089784.013199.0163.0FRFrance
211519844638733067686.0106974.0159123.0195.0FRFrance
21161984453135223101414.0169032.0246184.0308.0FRFrance
211719844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2118 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202522 3 20831 15052.0 26610.0 31 22.0 \n", + "1 202521 3 24065 18883.0 29247.0 36 28.0 \n", + "2 202520 3 20265 15814.0 24716.0 30 23.0 \n", + "3 202519 3 16264 12394.0 20134.0 24 18.0 \n", + "4 202518 3 18115 13975.0 22255.0 27 21.0 \n", + "5 202517 3 22150 17291.0 27009.0 33 26.0 \n", + "6 202516 3 28564 22550.0 34578.0 43 34.0 \n", + "7 202515 3 35721 29592.0 41850.0 53 44.0 \n", + "8 202514 3 37579 31232.0 43926.0 56 47.0 \n", + "9 202513 3 39673 33686.0 45660.0 59 50.0 \n", + "10 202512 3 52543 45627.0 59459.0 78 68.0 \n", + "11 202511 3 59469 52154.0 66784.0 89 78.0 \n", + "12 202510 3 60334 53048.0 67620.0 90 79.0 \n", + "13 202509 3 84531 74994.0 94068.0 126 112.0 \n", + "14 202508 3 136020 124824.0 147216.0 203 186.0 \n", + "15 202507 3 208952 195988.0 221916.0 312 293.0 \n", + "16 202506 3 273519 258159.0 288879.0 408 385.0 \n", + "17 202505 3 334395 318416.0 350374.0 499 475.0 \n", + "18 202504 3 350043 332885.0 367201.0 522 496.0 \n", + "19 202503 3 252772 238917.0 266627.0 377 356.0 \n", + "20 202502 3 257247 242991.0 271503.0 384 363.0 \n", + "21 202501 3 231549 214627.0 248471.0 345 320.0 \n", + "22 202452 3 201726 185870.0 217582.0 302 278.0 \n", + "23 202451 3 201697 187843.0 215551.0 302 281.0 \n", + "24 202450 3 136694 126369.0 147019.0 205 190.0 \n", + "25 202449 3 108487 99037.0 117937.0 163 149.0 \n", + "26 202448 3 87381 78687.0 96075.0 131 118.0 \n", + "27 202447 3 76286 67626.0 84946.0 114 101.0 \n", + "28 202446 3 56399 49006.0 63792.0 85 74.0 \n", + "29 202445 3 47347 40843.0 53851.0 71 61.0 \n", + "... ... ... ... ... ... ... ... \n", + "2088 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "2089 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "2090 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "2091 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "2092 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "2093 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2094 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2095 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2096 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2097 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2098 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2099 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2100 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2101 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2102 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2103 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2104 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2105 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2106 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2107 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2108 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2109 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2110 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2111 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2112 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2113 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2114 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2115 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2116 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2117 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 40.0 FR France \n", + "1 44.0 FR France \n", + "2 37.0 FR France \n", + "3 30.0 FR France \n", + "4 33.0 FR France \n", + "5 40.0 FR France \n", + "6 52.0 FR France \n", + "7 62.0 FR France \n", + "8 65.0 FR France \n", + "9 68.0 FR France \n", + "10 88.0 FR France \n", + "11 100.0 FR France \n", + "12 101.0 FR France \n", + "13 140.0 FR France \n", + "14 220.0 FR France \n", + "15 331.0 FR France \n", + "16 431.0 FR France \n", + "17 523.0 FR France \n", + "18 548.0 FR France \n", + "19 398.0 FR France \n", + "20 405.0 FR France \n", + "21 370.0 FR France \n", + "22 326.0 FR France \n", + "23 323.0 FR France \n", + "24 220.0 FR France \n", + "25 177.0 FR France \n", + "26 144.0 FR France \n", + "27 127.0 FR France \n", + "28 96.0 FR France \n", + "29 81.0 FR France \n", + "... ... ... ... \n", + "2088 59.0 FR France \n", + "2089 64.0 FR France \n", + "2090 97.0 FR France \n", + "2091 93.0 FR France \n", + "2092 80.0 FR France \n", + "2093 116.0 FR France \n", + "2094 149.0 FR France \n", + "2095 281.0 FR France \n", + "2096 395.0 FR France \n", + "2097 485.0 FR France \n", + "2098 544.0 FR France \n", + "2099 689.0 FR France \n", + "2100 722.0 FR France \n", + "2101 762.0 FR France \n", + "2102 926.0 FR France \n", + "2103 1113.0 FR France \n", + "2104 1236.0 FR France \n", + "2105 832.0 FR France \n", + "2106 459.0 FR France \n", + "2107 207.0 FR France \n", + "2108 190.0 FR France \n", + "2109 198.0 FR France \n", + "2110 224.0 FR France \n", + "2111 266.0 FR France \n", + "2112 219.0 FR France \n", + "2113 176.0 FR France \n", + "2114 163.0 FR France \n", + "2115 195.0 FR France \n", + "2116 308.0 FR France \n", + "2117 213.0 FR France \n", + "\n", + "[2118 rows x 10 columns]" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data = pd.read_csv(data_url, skiprows=1)\n", "raw_data" @@ -78,9 +1043,73 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
18811989193-NaNNaN-NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1881 198919 3 - NaN NaN - NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1881 FR France " + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -94,9 +1123,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 5, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020252232083115052.026610.03122.040.0FRFrance
120252132406518883.029247.03628.044.0FRFrance
220252032026515814.024716.03023.037.0FRFrance
320251931626412394.020134.02418.030.0FRFrance
420251831811513975.022255.02721.033.0FRFrance
520251732215017291.027009.03326.040.0FRFrance
620251632856422550.034578.04334.052.0FRFrance
720251533572129592.041850.05344.062.0FRFrance
820251433757931232.043926.05647.065.0FRFrance
920251333967333686.045660.05950.068.0FRFrance
1020251235254345627.059459.07868.088.0FRFrance
1120251135946952154.066784.08978.0100.0FRFrance
1220251036033453048.067620.09079.0101.0FRFrance
1320250938453174994.094068.0126112.0140.0FRFrance
142025083136020124824.0147216.0203186.0220.0FRFrance
152025073208952195988.0221916.0312293.0331.0FRFrance
162025063273519258159.0288879.0408385.0431.0FRFrance
172025053334395318416.0350374.0499475.0523.0FRFrance
182025043350043332885.0367201.0522496.0548.0FRFrance
192025033252772238917.0266627.0377356.0398.0FRFrance
202025023257247242991.0271503.0384363.0405.0FRFrance
212025013231549214627.0248471.0345320.0370.0FRFrance
222024523201726185870.0217582.0302278.0326.0FRFrance
232024513201697187843.0215551.0302281.0323.0FRFrance
242024503136694126369.0147019.0205190.0220.0FRFrance
25202449310848799037.0117937.0163149.0177.0FRFrance
2620244838738178687.096075.0131118.0144.0FRFrance
2720244737628667626.084946.0114101.0127.0FRFrance
2820244635639949006.063792.08574.096.0FRFrance
2920244534734740843.053851.07161.081.0FRFrance
.................................
208819852132609619621.032571.04735.059.0FRFrance
208919852032789620885.034907.05138.064.0FRFrance
209019851934315432821.053487.07859.097.0FRFrance
209119851834055529935.051175.07455.093.0FRFrance
209219851733405324366.043740.06244.080.0FRFrance
209319851635036236451.064273.09166.0116.0FRFrance
209419851536388145538.082224.011683.0149.0FRFrance
20951985143134545114400.0154690.0244207.0281.0FRFrance
20961985133197206176080.0218332.0357319.0395.0FRFrance
20971985123245240223304.0267176.0445405.0485.0FRFrance
20981985113276205252399.0300011.0501458.0544.0FRFrance
20991985103353231326279.0380183.0640591.0689.0FRFrance
21001985093369895341109.0398681.0670618.0722.0FRFrance
21011985083389886359529.0420243.0707652.0762.0FRFrance
21021985073471852432599.0511105.0855784.0926.0FRFrance
21031985063565825518011.0613639.01026939.01113.0FRFrance
21041985053637302592795.0681809.011551074.01236.0FRFrance
21051985043424937390794.0459080.0770708.0832.0FRFrance
21061985033213901174689.0253113.0388317.0459.0FRFrance
210719850239758680949.0114223.0177147.0207.0FRFrance
210819850138548965918.0105060.0155120.0190.0FRFrance
210919845238483060602.0109058.0154110.0198.0FRFrance
2110198451310172680242.0123210.0185146.0224.0FRFrance
21111984503123680101401.0145959.0225184.0266.0FRFrance
2112198449310107381684.0120462.0184149.0219.0FRFrance
211319844837862060634.096606.0143110.0176.0FRFrance
211419844737202954274.089784.013199.0163.0FRFrance
211519844638733067686.0106974.0159123.0195.0FRFrance
21161984453135223101414.0169032.0246184.0308.0FRFrance
211719844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2117 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202522 3 20831 15052.0 26610.0 31 22.0 \n", + "1 202521 3 24065 18883.0 29247.0 36 28.0 \n", + "2 202520 3 20265 15814.0 24716.0 30 23.0 \n", + "3 202519 3 16264 12394.0 20134.0 24 18.0 \n", + "4 202518 3 18115 13975.0 22255.0 27 21.0 \n", + "5 202517 3 22150 17291.0 27009.0 33 26.0 \n", + "6 202516 3 28564 22550.0 34578.0 43 34.0 \n", + "7 202515 3 35721 29592.0 41850.0 53 44.0 \n", + "8 202514 3 37579 31232.0 43926.0 56 47.0 \n", + "9 202513 3 39673 33686.0 45660.0 59 50.0 \n", + "10 202512 3 52543 45627.0 59459.0 78 68.0 \n", + "11 202511 3 59469 52154.0 66784.0 89 78.0 \n", + "12 202510 3 60334 53048.0 67620.0 90 79.0 \n", + "13 202509 3 84531 74994.0 94068.0 126 112.0 \n", + "14 202508 3 136020 124824.0 147216.0 203 186.0 \n", + "15 202507 3 208952 195988.0 221916.0 312 293.0 \n", + "16 202506 3 273519 258159.0 288879.0 408 385.0 \n", + "17 202505 3 334395 318416.0 350374.0 499 475.0 \n", + "18 202504 3 350043 332885.0 367201.0 522 496.0 \n", + "19 202503 3 252772 238917.0 266627.0 377 356.0 \n", + "20 202502 3 257247 242991.0 271503.0 384 363.0 \n", + "21 202501 3 231549 214627.0 248471.0 345 320.0 \n", + "22 202452 3 201726 185870.0 217582.0 302 278.0 \n", + "23 202451 3 201697 187843.0 215551.0 302 281.0 \n", + "24 202450 3 136694 126369.0 147019.0 205 190.0 \n", + "25 202449 3 108487 99037.0 117937.0 163 149.0 \n", + "26 202448 3 87381 78687.0 96075.0 131 118.0 \n", + "27 202447 3 76286 67626.0 84946.0 114 101.0 \n", + "28 202446 3 56399 49006.0 63792.0 85 74.0 \n", + "29 202445 3 47347 40843.0 53851.0 71 61.0 \n", + "... ... ... ... ... ... ... ... \n", + "2088 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "2089 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "2090 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "2091 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "2092 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "2093 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2094 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2095 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2096 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2097 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2098 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2099 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2100 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2101 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2102 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2103 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2104 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2105 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2106 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2107 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2108 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2109 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2110 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2111 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2112 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2113 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2114 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2115 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2116 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2117 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 40.0 FR France \n", + "1 44.0 FR France \n", + "2 37.0 FR France \n", + "3 30.0 FR France \n", + "4 33.0 FR France \n", + "5 40.0 FR France \n", + "6 52.0 FR France \n", + "7 62.0 FR France \n", + "8 65.0 FR France \n", + "9 68.0 FR France \n", + "10 88.0 FR France \n", + "11 100.0 FR France \n", + "12 101.0 FR France \n", + "13 140.0 FR France \n", + "14 220.0 FR France \n", + "15 331.0 FR France \n", + "16 431.0 FR France \n", + "17 523.0 FR France \n", + "18 548.0 FR France \n", + "19 398.0 FR France \n", + "20 405.0 FR France \n", + "21 370.0 FR France \n", + "22 326.0 FR France \n", + "23 323.0 FR France \n", + "24 220.0 FR France \n", + "25 177.0 FR France \n", + "26 144.0 FR France \n", + "27 127.0 FR France \n", + "28 96.0 FR France \n", + "29 81.0 FR France \n", + "... ... ... ... \n", + "2088 59.0 FR France \n", + "2089 64.0 FR France \n", + "2090 97.0 FR France \n", + "2091 93.0 FR France \n", + "2092 80.0 FR France \n", + "2093 116.0 FR France \n", + "2094 149.0 FR France \n", + "2095 281.0 FR France \n", + "2096 395.0 FR France \n", + "2097 485.0 FR France \n", + "2098 544.0 FR France \n", + "2099 689.0 FR France \n", + "2100 722.0 FR France \n", + "2101 762.0 FR France \n", + "2102 926.0 FR France \n", + "2103 1113.0 FR France \n", + "2104 1236.0 FR France \n", + "2105 832.0 FR France \n", + "2106 459.0 FR France \n", + "2107 207.0 FR France \n", + "2108 190.0 FR France \n", + "2109 198.0 FR France \n", + "2110 224.0 FR France \n", + "2111 266.0 FR France \n", + "2112 219.0 FR France \n", + "2113 176.0 FR France \n", + "2114 163.0 FR France \n", + "2115 195.0 FR France \n", + "2116 308.0 FR France \n", + "2117 213.0 FR France \n", + "\n", + "[2117 rows x 10 columns]" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -122,7 +2118,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 6, "metadata": {}, "outputs": [], "source": [ @@ -152,10 +2148,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 7, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -179,9 +2173,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n" + ] + } + ], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", @@ -199,9 +2201,41 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, "outputs": [], + "source": [ + "sorted_data['inc'] = pd.to_numeric(sorted_data['inc'], errors='coerce')" + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 13, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXm4HUWZ/7/vWe65W272hGyQIGFfJUZQQNlRcXBGmAkzSlQUdVBxGRX8qcygzKAzCuKCwwgIuLCJAiJiCIsCgZBASEggJGTfl5vl7vcs9fujq/pUd1cv55w+y715P89zn3NudddyqrvrrXepahJCgGEYhmEqJVHvBjAMwzDDAxYoDMMwTCywQGEYhmFigQUKwzAMEwssUBiGYZhYYIHCMAzDxAILFIZhGCYWWKAwDMMwscAChWEYhomFVL0bUEvGjRsnpk+fXu9mMAzDDCkWL168SwgxPuy8A0qgTJ8+HYsWLap3MxiGYYYURLQ+ynls8mIYhmFigQUKwzAMEwssUBiGYZhYYIHCMAzDxAILFIZhGCYWQgUKEd1ORDuI6DUtbQwRzSOiVfJztHbsGiJaTUQrieh8Lf1kIlomj91MRCTTM0R0r0x/kYima3nmyjpWEdFcLX2GPHeVzNtUeVcwDMMwlRBFQ/klgAtcaVcDmC+EmAlgvvwfRHQ0gDkAjpF5fkZESZnnFgBXAJgp/1SZlwPYI4Q4DMCNAL4nyxoD4FoA7wQwG8C1muD6HoAbZf17ZBkMwzBMHQkVKEKIvwLodCVfBOBO+f1OAB/S0u8RQgwIIdYCWA1gNhFNAtAhhFggrHcO3+XKo8p6AMDZUns5H8A8IUSnEGIPgHkALpDHzpLnuusfdixevwcrtuyvdzMYhmFCKXdh40QhxFYAEEJsJaIJMn0KgBe08zbJtKz87k5XeTbKsnJEtA/AWD3dlWcsgL1CiJyhrGHHh295HgCw7oYP1LklDMMwwcTtlCdDmghILydPUFneBhFdQUSLiGjRzp07/U5jGIZhKqRcgbJdmrEgP3fI9E0ApmnnTQWwRaZPNaQ78hBRCsBIWCY2v7J2ARglz3WX5UEIcasQYpYQYtb48aFb0TAMwzBlUq5AeRiAirqaC+AhLX2OjNyaAcv5vlCax7qI6BTpA7nMlUeVdTGAJ6Wf5XEA5xHRaOmMPw/A4/LYU/Jcd/0MwzBMnQj1oRDRbwG8F8A4ItoEK/LqBgD3EdHlADYAuAQAhBDLieg+ACsA5ABcKYTIy6I+CytirAXAY/IPAG4DcDcRrYalmcyRZXUS0XcAvCTPu04IoYIDvg7gHiL6LoBXZBkMwzBMHSFrwn9gMGvWLDHUdhuefvWjANgpzzBM/SCixUKIWWHn8Up5hmEYJhZYoDAMwzCxwAKFYRiGiQUWKAzDMEwssEBhGIZhYoEFCsMwDBMLLFAYhmGYWGCBwjAMw8QCCxSGYRgmFligMAzDMLHAAoVhGIaJBRYoDMMwTCywQGEYhmFigQUKwzAMEwssUJiSuW/RRky/+lF09Wfr3RSGYRoIFihMyfzfX9cAALbs7a9zSxiGaSRYoDAMwzCxwAKFYRiGiQUWKEzJENW7BQzDNCIsUJiyERD1bgLDMA0ECxSmZAisojAM44UFClM2ghUUhmE0WKAwDMMwscAChWEYhokFFihMyXCUF8MwJligMGXDPhSGYXRYoDAMwzCxwAKFKRteh8IwjA4LFIZhGCYWWKAMYXZ09eNt3/gTXtmwp95NYRiGqUygENGXiGg5Eb1GRL8lomYiGkNE84holfwcrZ1/DRGtJqKVRHS+ln4yES2Tx24msuKIiChDRPfK9BeJaLqWZ66sYxURza3kdzQqIsTr/fzq3cgXBO54bl1tGiQhDvNiGMZA2QKFiKYA+AKAWUKIYwEkAcwBcDWA+UKImQDmy/9BREfL48cAuADAz4goKYu7BcAVAGbKvwtk+uUA9gghDgNwI4DvybLGALgWwDsBzAZwrS64DhTqPa5zlBfDMDqVmrxSAFqIKAWgFcAWABcBuFMevxPAh+T3iwDcI4QYEEKsBbAawGwimgSgQwixQFhT8rtceVRZDwA4W2ov5wOYJ4ToFELsATAPRSE0bGjUAZv1E4ZhTJQtUIQQmwH8D4ANALYC2CeE+AuAiUKIrfKcrQAmyCxTAGzUitgk06bI7+50Rx4hRA7APgBjA8piGIZh6kQlJq/RsDSIGQAmA2gjoo8EZTGkiYD0cvO423kFES0iokU7d+4MaF7j0aAKSt1NbQzDNCaVmLzOAbBWCLFTCJEF8CCAdwHYLs1YkJ875PmbAEzT8k+FZSLbJL+70x15pFltJIDOgLI8CCFuFULMEkLMGj9+fJk/tbFpVMHDMMyBRSUCZQOAU4ioVfo1zgbwOoCHAaioq7kAHpLfHwYwR0ZuzYDlfF8ozWJdRHSKLOcyVx5V1sUAnpR+lscBnEdEo6WmdJ5MG1aERXkxDMM0EqlyMwohXiSiBwC8DCAH4BUAtwJoB3AfEV0OS+hcIs9fTkT3AVghz79SCJGXxX0WwC8BtAB4TP4BwG0A7iai1bA0kzmyrE4i+g6Al+R51wkhOsv9LUOdegkelncMw+iULVAAQAhxLazwXZ0BWNqK6fzrAVxvSF8E4FhDej+kQDIcux3A7SU2eUgRNl7Xaz0I+1AYhjHBK+WZsuG9vBiG0WGB0sBENSnVeljnd8ozDGOCBcoQhod1hmEaCRYoDQyblBiGGUqwQGHKppwor+6BHC699QWs390Tf4MYhqkrLFAamEYNy60kymveim1YsGY3fjjvzfgaxDBMQ8ACZThQJ8FTSbXs/2GY4QcLlCGM0hRq7WthYcAwjAkWKA1MmMlrKIbvNqoZj2GYymGBwtQFfusjwww/WKA0MI0eNsybVzIMo8MCZRhQ83GdtQuGYQywQGlgQn0o9X6nfDl5WKlhmGELCxSmZOKQY6zjMMzwgwVKA8OTeYZhhhIsUJiawkKSYYYvLFAamLAoKrLPq35bTFRUL9u8GGbYwQJlCFO3lfIsDBiGMcACpYFpfPNQ47eQYZjawQKFKZlKTG28GJJhhi8sUBqY4Tz2DsV9yBiGCYYFyjCgbk75+lTLMEyDwgKlkQkdses7yy/L5BV/MxiGaRBYoDAlo3YKrsQfwpFiDDP8YIHSwDT8bsP1bgDDMA0FC5RhQM03G1b1VqHi7oEc9vVm4y+YYZiqwwKlBJ58Yzte3bi3ZvVF3W14SEWDhbT15O/MwwnX/aU2bWEYJlZS9W7AUOITv1wEAFh3wwfq3BKLershKjHJ+bV9IFcou0yGYeoLaygNTMMrHg3fQIZhagkLFKZkinuIlU6jBxowDFM+LFAamEbfpqSS5nHYMMMMPyoSKEQ0iogeIKI3iOh1IjqViMYQ0TwiWiU/R2vnX0NEq4loJRGdr6WfTETL5LGbSS50IKIMEd0r018koulanrmyjlVENLeS3zH0qfFuw9IDwtoGwzA6lWooPwLwZyHEkQBOAPA6gKsBzBdCzAQwX/4PIjoawBwAxwC4AMDPiCgpy7kFwBUAZsq/C2T65QD2CCEOA3AjgO/JssYAuBbAOwHMBnCtLriGC40+XDe4AsUwTI0pW6AQUQeAMwDcBgBCiEEhxF4AFwG4U552J4APye8XAbhHCDEghFgLYDWA2UQ0CUCHEGKBsGw8d7nyqLIeAHC21F7OBzBPCNEphNgDYB6KQuiAYSiO5yyEGGb4UomGciiAnQDuIKJXiOgXRNQGYKIQYisAyM8J8vwpADZq+TfJtCnyuzvdkUcIkQOwD8DYgLKGFWGDrzo+FDeH5N2GGWb4UYlASQF4O4BbhBAnAeiBNG/5YBpBREB6uXmclRJdQUSLiGjRzp07A5o3FKnvdL/RgwYYhqktlQiUTQA2CSFelP8/AEvAbJdmLMjPHdr507T8UwFskelTDemOPESUAjASQGdAWR6EELcKIWYJIWaNHz++jJ9ZP6I6vWs+rFcUNswwzHClbIEihNgGYCMRHSGTzgawAsDDAFTU1VwAD8nvDwOYIyO3ZsByvi+UZrEuIjpF+kcuc+VRZV0M4EnpZ3kcwHlENFo648+TaQcUdVcQOGyYYRiNSrde+TyAXxNRE4A1AD4OS0jdR0SXA9gA4BIAEEIsJ6L7YAmdHIArhRB5Wc5nAfwSQAuAx+QfYDn87yai1bA0kzmyrE4i+g6Al+R51wkhOiv8LY1HmA+lNq3wYG8OyfoGwzAaFQkUIcQSALMMh872Of96ANcb0hcBONaQ3g8pkAzHbgdweyntHW7UXUMpg6HYZoZhosEr5RuYsLG33hoCr5SPh2seXIYjvvlY+IkM0+DwbsPDgHpFW5X3CmBWUdz8duGGejeBYWKBNZQGJuo6lFpTyeaQWikxtKTx2d09gBfW7K53MximJrBAaWDCZvPC9VlrytGMDjQfypxbX8CcW1/gNTvMAQELlCFMvQap4uaQFZRxYCgoWLWjGwCQLzSeQHlpXSd2dw/UuxnMMIIFSgPDk9p4+cXf1uBPy7bWpe5cgwkUIQQu+fkCXPp/L9S7Kcwwgp3yw4C67eVVllO+fnz30dcB1OcVzoP5AprTyfATa4QScG9u765zS5jhBGsoDUxo2HDdJ73xv1N+uJLNFerdBAcDsj3JxNC6EnNuXYBZ351X72YwPrBAGcLUKwTXjvKqu0Dzks0XsHVfX72b4aHRTF6DUqA0JYfWEPDCmk7s6h6sdzMYH4bW3XSAEeZ0r/eA3lhDpMU3HlyGU//rSfQM5OrdFAeDDaehWLsepZNDS0NhGhsWKEOYuq9DKWu74eo2+vHl2wBYmkojUai39HeRy1vtGWomL6axYYHSwIQubHR91oo4xsZqhQ0r01KjDZSNFjbcYPKNGSawQBkiBJm/6rb1SgMavZRAibtl+/uzFa14bzB50pDXjhn6sEAZwtR79XUjhg3npKlLxGzx+tdfvYw5t76AfX3ZsvJXy+T11s5uvLSu9Dc32K+Pjrk9zIENr0MZIgjhNRPVezBoxHfKF+yBMt7eeWNbFwDlzE6XnL9aAuXsHzwDoPy1NWz6YuKENZSYEELgx/NXYdOe3hjL1L4bT4itqrKot4YURNxNUy6ZQomaj8oXxYdSy/6s55W76p5X8O2HXqtjC5hqwQIlJtbv7sUP5r2JT921uGZ11nsdSiMTd88oJ3+pmgbJzoqSrZZ+FiW86jEpeGjJFty1YH3N62WqDwuUmFCPZe9gfOsfdIHRSNpAJU2p1c+Iu78SUjCUGq2lZG8UQXSgaCjM8IUFShmYHnz7Pes1nWXWrq6466+2lhP3bL/ctTelCKJ6RIKZquwbzGPtrp6at6UUGmmCxRRhgVIGpsFBDRxxOl/DfCjCcF4taeTQ07jbVvb1Vb6XSCavGmooAYuY5t6+EGf+z9M1a0s5NFoYNmPBAqUM8iYNpQ77W9V7klZW2HDNbF7xFqd8KKZrH4TtzI9k8iq5WRXgX9nCMsKQy2pBBT+40XYeYCxYoJSBKdJHaS1xDph6SaZi660hVGTyiq8ZRuLumWKUV6k+FIqcrx4aSlCN1Rb+lRTPAqUxYYFSBqabWc1cq3WbBwmPWguWeN4pX13iDxu2fnQ2X56GEkWzqWV/Rqmr2tvFVCIUWJ40JixQysA0OBRsDSW+eqLuNly/F2w17lMdt5ClEkxXznylhA3Xvj+DrmG1/RSVlN/At94BDQuUMjCZL5SQqdagYDZ51YeKwobja0YgjeK0LSlsuIYbJEcxeVVbwFVSPpu8GhMWKGVgMgWotDjDYUMfGRG/VlQKlVRLVY4bjnszTdsXMlyivCJcvWq3h30oww8WKGVgNnlZn9Xao8pE3R+pujfAn6DxphLtpeS8JZxfUx9KBHNp9U1elWgoMTaEiQ0WKGVgeg6UkIlVQ4n40NQr2qucehthYlmJ76fcvFHy1cWHEnANG9nk1ciTmQMZFihlEGTySlTJlGP0odTJKd/I75RXNIqGErB+0FB2HRY2Bp1TZZ9OZdehgW++AxgWKGVgEijVucHDorzqvA6lRnnKIeh6lHOtikK0er+gpotiEe5/K3URZ8ltYKf8sKNigUJESSJ6hYj+KP8fQ0TziGiV/BytnXsNEa0mopVEdL6WfjIRLZPHbibpsSWiDBHdK9NfJKLpWp65so5VRDS30t9RCsZ1KFVwyuuYTBNDMsqrRgNB8IK96pRbacZG22Wh+iav+uRlqkccGspVAF7X/r8awHwhxEwA8+X/IKKjAcwBcAyACwD8jIiSMs8tAK4AMFP+XSDTLwewRwhxGIAbAXxPljUGwLUA3glgNoBrdcFVbUw3c6EaUV7C/N2dNhQFS7UJjPKqoMdKXSlfigCtjw/Fn+pHedXel8VUl4oEChFNBfABAL/Qki8CcKf8fieAD2np9wghBoQQawGsBjCbiCYB6BBCLBDWXXKXK48q6wEAZ0vt5XwA84QQnUKIPQDmoSiEqo7Rh6Kc8jWM8rKpW9hw+U75au82HDxQVqfc4HzVdcqXPcAG+Zoa2IfC4qQxqVRDuQnA1wDot95EIcRWAJCfE2T6FAAbtfM2ybQp8rs73ZFHCJEDsA/A2ICyakKtTF7C57s7beG6Tmzf3x9fxREpa3PIGg0FwU752tnuSzm7MlNi/Oc3sobCPpTGpGyBQkQXAtghhIj6ikLTUCsC0svN46yU6AoiWkREi3bu3BmpoSb0mz/IKV9L/URv03V/XFGzeivZy6tmmw0HmbwqmHmXrQhUeQAvXdAJx2fc7YnCUPKhLF7fiUde3VLbSocglWgo7wbwd0S0DsA9AM4iol8B2C7NWJCfO+T5mwBM0/JPBbBFpk81pDvyEFEKwEgAnQFleRBC3CqEmCWEmDV+/PjyfimcA4I5bLjsoiPV2ZA24zLaVJwpVHmlfOCxSkxL1Tu/spXjpZ0fKWy4yrdcRQK0xhLlw7cswOd/+0pN6xyKlC1QhBDXCCGmCiGmw3K2PymE+AiAhwGoqKu5AB6S3x8GMEdGbs2A5XxfKM1iXUR0ivSPXObKo8q6WNYhADwO4DwiGi2d8efJtKrhMD+ZnPL2wsYarpSvl++khqaZatRTj/UPjaahKALDhht4t2GmManGOpQbAJxLRKsAnCv/hxBiOYD7AKwA8GcAVwoh8jLPZ2E59lcDeAvAYzL9NgBjiWg1gC9DRowJIToBfAfAS/LvOplWNRwmL+PWK/GbvBzvlA85t9R6CwWB+xdtRK4C1aosk1fEXJVqZNUy5ZRrWopWdqmt0erxybt5b5+xL4Xr09yeavtQys/LwqgxScVRiBDiaQBPy++7AZztc971AK43pC8CcKwhvR/AJT5l3Q7g9nLbXAlBUV61dKJUYrr5/Sub8dUHlmJH1wCuPPOw8uovxykfMY8QlQU4uOvRB9WKoqmqmC9uJ/Vrm/fhwh8/i+suOgaXnTq95Lp4Ly+mVHilfET0+zcwyivOOoX5e1BaVLr6swBQUXRYJQNgmLCodLzwCpTKCi++16REDaUUH0pJJTsx3ZNrdvUAABau9SrvtoYS844CpcBbrww/WKBERL9/B3NeM1E1fChhg2Alj1QqaV36Ut9ACFQa5VUbk5d7wCk4NJQKyi3TQljtzSFNv0kE3JONEDZcl3U3TFVhgRIR3bzUO5j3HM/b29fXDv2ZKlWQNUmBUpEPpYxnWg18Ya2Ne7jQB9xamrxKOb+ShYRGP0mEvg70oVR5YWNlK+VjbAgTGyxQIqLfwL2DOc/xqmy94nDKx/sEpZJWQ3NlTNfjiPIK66dKZ8cek5fWf3WJpqpi2VZeU51WYsLY19axai0AjcJQWofCRIMFShkMZL1TN+WUr9b29SYqETJFk1ejRnmVUXhAPWH+qDDUVS15Vl2KDyVmn4L90rdhaPJiH0o4c25dgDufX1fTOlmgRMSxsDHAKV+tOuN2yiepzFfaOuqvrXZTST3676zHOppICwkr2bQyIDQ48rYSnjLLbk4kKjGpsUAJ54U1nbj24eU1rZMFSkT0h91kJqrHwkadUmutUzNDtpWPZ9AHTE55/2OllVva+fVch2L3Z5nXuurrUGq4YwFTG1igRES/gU3bPlQjbNhRv7FNlT9VZZl/Knljo/C33Ts0sgp9Ru7cziivSpzy5eatbpSXSUMuaihlmryqrKIMhR0XmNJggVIGJg2l6i/YivkJUo7aWg+uwvWpE5dZypTfKazKpxp7ZhXPjdmnIJNMTvmhv7CRJUojwgIlIvrta5q5FU1eMdYZ8sxU9kyV39BqzSyFz/cya3LVK4zfSy613CivKD6UmPs16J6sdtRZFHhh4/CDBUpE9IHErKFYn+XuopsvCOzo8l+1bjR5lVWTq4wamx2C3mXuDELwLzxfEPj8b1/Ba5v3RW6b04cSqamRyg09v4Rz4x5gG93kxVuvDD9YoETEoaGYQjQr1FB+9MSbmH39fGzbVxQqYSYl58LG0upTZpBKnssKXCjG3+YweQWUsX53Dx55dQs+95uXI7etUh8KVWgirLZGYF4pb32We09We9CuzIzLEqURYYESEf3ezxm2K6nUKf/USuvlX357a5lN5OU/VMW9qcouokwNJVreoOOqr5PmFXsAvLPruPwzpftQomeIO4Rb3R9mk1f5QQLLt+zDSdf9Bcu3+GuIUeCFjcMPFihR0W7gwHUoZU4H1eCYj2ngi075s/VK3ikfeixIoMgTUwn/29eT3WHyGn4+lGANpTwnil8/Pb1yJ/b0ZvHYsm0ltNBQfgVSwS/v4vWdOO7ax7G3d7DssoPgPcSCYYESEX3wzBtWZNkPn88Nt7GzFxt29/qWbwsUnwfFNHjHcW/Xy4diwmny8j9PaYiJAA0lyIdSy1DpUk6v7Hqa7g9/rTkOE5xpYlUK1dBQfvLkanQN5PDyhj0llTf39oX482vhAnKoyJN6CT4WKBFxrJQ3rPBVgsDvRj/9+0/hjP9+yrd8k0AJuyUqskCr9SAVlFFexc76DYfkcf8iCraGEiBQXL8srnUo5ftQyjcxRSrfaBK1MCookZzy5nR1r5ajYdQ72s6PZ97cic/8anHoeZUK0VpRL5MgC5QyML4PRQjfY1FIhWgoPisb7a8mv04QhYCBXbF5bx8Wr/fO9CrSaoKORXTK5yL4UNwFxLb1Sqnnl5Ahbqe8GvDL3V/Orz1q256gjUVfWLMbb2zb70l3LBCuwnUoxzdYKz9XLalXO2N5Y+OBgH55TIN3IURDCUMNjvpmjWEDrJ5W+iaP4RrK6d97EgUBrLvhA+YSyrhpgzQjp1nKv2wldN0aivOtjO569XrqoKFU2YcSpAGZxEmU3+F3TiJs8gNgzq0vAPDeO3qOal6H0gR59HOHiDypm0BhDSUiYa+QVeN5uaq4mkWWexuUO3MONi1FKyO2eoXxqwc1kLl9KI51LAG7Dfv9rpueeBO3PbvWeEyt5aimKaGizSEN84miycsrUqKYbvx+qxLkv3lxQ9TmFdsUm+nRnG7vCl1SWcNPQ6lXM1lDiYhDQwlwyld8w/kMquELAUurxjZ5lZYNQIVvbLQ/TT4UbbAJGLnVgka3hhJk1nIeM5d90xOrAACXnzbDt+5yn9QouSp6wZZxXU9QXeVrKGn56oPBMl59UGlwRLEcc+ZyrHul7BQ+VMKVWUNpcKI65cu9juU8CI4XcJVYcXHFevk3XiUaStixoBn0dx99HYDXh+IYrFx5oi6aDKOaA0rcTvmgl75FGUT9TpHyBONHZKI2zyauF535XcSgIu99aQM27+0rKY+boaKhsFO+wQkLG67UKW+qJ8iE4z1eYj0x3HDlbQ7pL3j1vosy4AVrKMJ1TPteyfqHsn0o5Q/g0cr3pimnuckpH8mH4tMgNaEKionwoxKtOoq5TD2H7qbt7h7A13+3DFfctciTp5RrKqr8WuS4YA2l0QnRUAqVaiiqmlIchPr3kjWU6PW5y65oYhlgatPToph/kq6FjUECNshhXwrVfE79ruFDSzajZ8D72mlH3oCtbEzjflHY+JfpNyip9Eo11FIHvSh+MNtS4ErfstfagWJvb9abp4R2DJWw4cFcfSQfC5SI6LdR0MLGcm84c7ijv08gqH1RKEZbhefM+oQkx/1sRTV5KTxRXgEmwEqEr041Z36mkl/fuh9X3bMEX/vd0sC8pgE26JUKxWMB29f4mpXMg3YUnCavUvMW8dVQ7GhL5/Gt+yxT14hmr9u4FK1jqJi8oizSrAYsUCLiHOy8x5XWUumrgEtSvyPM2MLyRqnOFIQAhA8of3hlM87+wdPGxWzmIIPSTF7JZIAPxZU9Lh9KNQl6E+hbO7oD85qEZNA6EVVuoIbia/IqX0OJGhpuzus/YXC3zd32fjljb0p5h7xGjfJ68OVNuObB4ImEH9PGtMbcmmiwQCmDIA0lzD4fetwnUslsIgp/wPwI8mW4cWsodp6QzF+8dwne2tnj4yz35o0yA9VJUglRXgXzeVGpeLfhKILbYEttTicBAL2D+eDyDWlBZtgor1vwN3kF1RpMJabHKP4Xle4uW/WtyZ9UilWhlgrKl+97Fb9duLGsvOxDaXCcTnnvxQrbesU+L+RC+0UqGQVGBfeMGmCjmLzcA50tjELyqWdX13CCNCNnJF0ZJi+tme7chQoGMmc55vSV27pw36LyHn6FacGsarZJ2DjP89dQglbRB0UX+v3WSnwozvfSlDcJcpej4xccoxb+mnZXaFQNpRLq5UPhdSgRCRvs1I0cJjDyBQE56XRA2nFTncY2+bQvChGVDABe00kp5jLAqR1E3RwyksmrBB+KXz2l4pf1/Jv+CgD4x1nTzPmi+KoMmq+wB8jS2xW0NirK9jV+/RQUWBFKJQIlQl4/H8qg2lDU8HNLWf8zVNahlL5zRjywhhIR/T4yDXYFH9ut5zyfB8FkUill/6nSH85oWgbgvTltYRSSWz27kTUU7Xskk1fAOhT3ZYhifw9C1RTWrkoe5KxhVhk1HN10tOjr8J8ABUX+hoXmliOY/cLiI+WNIFD8tCfVt+WGUNvnDhGJwlFeDY7DYWy0SUfTUIIcpUBpQkSUcK43r/oSfq7fxpNRn62CwRRlXlcTrqHo53gA4CJuAAAgAElEQVQFikNn821rJVaLMGHkF95bjiYIFPsurK9NA13QPRll48gwp3w5W05WYvKKcr66V/1MXhULlCFo8qrlVvYsUCLiNHkFhA2X7ZRXG+5pdYaM9k6nvfncvb2DePDlTd68th8kwkPq/r0yS9jDRfautPH5UPTkIIESFOVVzXUoAxXMDE1bmRRn3OVrKEEhxUFSwdeH4rPWIwr6vVCNsOGimc9drxSCw8jklcsXcP+ijcYxRb+XatnmsgUKEU0joqeI6HUiWk5EV8n0MUQ0j4hWyc/RWp5riGg1Ea0kovO19JOJaJk8djPJkYiIMkR0r0x/kYima3nmyjpWEdHccn9HOQQ65UOunt9AaYwiCtNQ9O8+537x3iX48n2v4q2dzrDTUvwgniiviBFitl/IoDmYuiGaQPEXGkFdV+1NCRV+Jq+gKlVwgXEX64i+uSA/iUkYFcOGS/ehhLYloJP031iJhuInBPx8KEFrcoaqhvLL59fhqw8sxb2GYBBdQ6l0KUMpVKKh5AB8RQhxFIBTAFxJREcDuBrAfCHETADz5f+Qx+YAOAbABQB+RkTKPX0LgCsAzJR/F8j0ywHsEUIcBuBGAN+TZY0BcC2AdwKYDeBaXXBVA4f91rSXlz1rDy6nlIcxzEwTRfZs22etEB7IOhtdyuaQ7oFO1RsmPNXDawo0MD2Yjki6kBmoqYwgG7ujryoSKOVNGIJQg7opkivqRMV0IYt+Pe+xilbKh7Qp6B533gsGIZgvYPrVj+IXf1vjORbFh7JKrtdxlx1lTU4UjMK5IPDTp1Zjf793FX416ewZdHzqODWUISBQhBBbhRAvy+9dAF4HMAXARQDulKfdCeBD8vtFAO4RQgwIIdYCWA1gNhFNAtAhhFggrKt1lyuPKusBAGdL7eV8APOEEJ1CiD0A5qEohKqCflGMuw1H9KH4zayKTt9imsOBGWr+Mh9Xye6ZWSmbQ7qjj1SOqGNn3iEkzTNId3l+/eQcVNz5zcLY/X85z5dwffrhu6tAYNnW0UGjhiLPKV2eBDrP7VcAlLFSPigcWS/bRDbEFKMGwv9+fKXnWCk+Q8/1V89nwOLRKJjaPP+NHfjvx1fiu39cEbmcarN8S/HlZmF+2ziJxYciTVEnAXgRwEQhxFbAEjoAJsjTpgDQdbNNMm2K/O5Od+QRQuQA7AMwNqAsU9uuIKJFRLRo586d5f1AuKK8DNcnsskrbIbr89CYo6LCHzB1jkeglHCPeTUU/4FKh2y/kHegN5u8StNQvBtA6iaRgGPlCBR7YA/O7LerQBCqPSYNRdVXjskryIdSXIdSvsnL77cGCZScYXJhymsaBMN8hqaJi7tes0Dxba6sV7svTf4KaV7q6g/eb61c/MaUoDVE+qGhYvICABBRO4DfAfiiEML7zk/tVEOaCEgvN48zUYhbhRCzhBCzxo8fH9C8YJxrJPzDO0MFhs8M1vah6CYBQ/k6UUxetobi6rJSwob9FtVFNR05HnQfG7e7LWERRqYy9H/dA1KQqSwKagALy+sXERfUV8IeoIOEQnC9psNBee2wYc9Ew18ou9N9/RgBbdX7J6jNYZqEqWk9gznf4+qZDfJ/+hEWmWa/H6hK43bYmBLkIwOGkEAhojQsYfJrIcSDMnm7NGNBfu6Q6ZsA6Ku+pgLYItOnGtIdeYgoBWAkgM6AsqqG02HsPa5vcxH4+lqfY8U3Appn4GGaj6/JS5XvMXnJciPca1nXSZHzGnwoQVFL0Zzy5u/e/AXfY+UIFHtnAZ+sxV0ByhFWFuYoL+enb/sCNRTDgONTYJBJsViu/AzxsZhwRvx5z/MzGQLhUV56yLa7DbaGYtLWSvCLmU5VfqhyfRX5gsCVv3kZyzbtC63fRJBgjpI/TiqJ8iIAtwF4XQjxQ+3QwwDmyu9zATykpc+RkVszYDnfF0qzWBcRnSLLvMyVR5V1MYAnpZ/lcQDnEdFo6Yw/T6ZVjbDBPR9xRhB2cf1u3jATkW+IZ4hvJcq0yrP1ij3IRdRQHKYs+Wk0aejnhc/03fU7/Vz+2ks5RBWiJu1Vz288Jg+aFjZWLWzYLtecrtftLbc42zcumgz0oQTfs0F5wzSU7v6c8Vy9XNP1CffHhPVJZa+H3rSnF48u3Yp//c3i0Pqd6dan+1XYgHPSW0unfCVbr7wbwEcBLCOiJTLtGwBuAHAfEV0OYAOASwBACLGciO4DsAJWhNiVQgi1491nAfwSQAuAx+QfYAmsu4loNSzNZI4sq5OIvgPgJXnedUKIzgp+Syjq4hEFO+X1c83HfDQUwywnbP8wHd+jPoO/X7y+CW/YsCojOJ+9Ul4fRAIGOccMNIKG4u7KIB9K2EAZOmDbJkKf6wer/UEz7KBy3W1URL5OAUIjKKTYnTGoD03tLAjAtelzoIlmIFvc5NLcrgAflH66SUMZ1Mt2l6sEirdYd+SZ268UNrErjuflDdwqMCIsFNqNbc4OCYWupVO+bIEihHgW/suizvbJcz2A6w3piwAca0jvhxRIhmO3A7g9ansrRQ0kqQQFzvgAw+yoBHum34AZ5nPwm2apVHe96vQo6rB3YWO0WbPCNMMLNXn5zsqcD7/zmN5m1zGHH8dbbthPUcf9TV4ECOHfn37JAQISiO5DCTJ5+YW6Wvnc6XqZ5rr0QTlXKCCZcG5OFyQT+nP+g77eZhNB5k7AqUl77o2AoBm973Z0DWBiR7PvcVM/2wKhSuO2f8SjVaFpx2iT37IW8Er5iKiLmkokjLMo/SFzPxSl2DP9TF5mE5H23ac8OyLLvNgdK7YGxVFYeKK85Gf4SnmZ3/D7jZFHEQRvkNkjKBonyJkPRI+iCns4S93LK+w3BwmysDDaYuSh6Zi3DAD4/SubjW1ztDlEODsjFZ1l9Gcr8aEEX8OgCLKiD8UkUIrf+wyvCXBMFgOuf7nrm8ImDWGTq0bSUFigRMTWUJJkfIj0G819A0QZKIuRIuaHImwhoK+vxKdNej2vbtxrzKvwGySj3qcmc6DZ7KSX7fN7As7R//MIFIcp0UtUk6LfWSbznjN/8PUBzANH1Bm76aygVfZ+prRt8s2Gfvnc6UH7hAHe9g/kgk1elWgoQWbnqNFjfVmvQBEhWpvtjzI1OgLq+QrbTsabLusP+U1Dwil/oKGuT1MyEbh9PeCdxej/hm9vr9XpKCN4AL7w+MnG8vxMW3resFm1n4M76joU08zROBt3+IzMZQZtvaIP5kEmr7AwSxNR196UOhsMMuGF1Rc24VD9EeQ4Dw5nNqc7NG6DANWPu/tD3xLEaLYKsJeZ3vzpzBtBQzENvlqaSaCE9XPUdyH5kbWvk/l4WESeaRITppFXCxYoEVEXJZUkn11h/WcEpdgz/RY2+g2wah+opmTwpfTO5ov/N5te0KLh+4KtMjSUoNfHOrSPCE5592/SBWOQhmI2txmr87TN34diffq+LtknX5hZM2gwCAtnLTrlvcfsfb58Sw9wygdo43rZVt1BAqU0DSXU3xQw2Qhah6In9YeZvIxmuoKss7yBW90zfgsV/Sah2/db2yoZN6vV/VwlBopUAguUiKhrmkokQm3wXkdnuPpZHLD048E3skDwalnrHPPMLOze19vs9wrgcnwoQVuvOAbXKPb7AIHiFvphM7Yghy6gmbxCNczyNRTj+qaA4pwDbNDM2X8C5D6UTBSHBF+nfIg5Jch+PxiysDHIhxLqBzNEE9rt8NnW3p0WpqEYBVnAJCkKYZGBpj6et2I7Hl221cofZjFhDaXxUBclnaTATfz0c+1jITZnPY/vOhQfh60KcfS7ZezBP8gMF3JDumfdRYHiU2lAWYECJUKYtNOH4jzm2BDP83uDB4WvPbDUt1xnu43Nss17pe7lFbZZYtR1GaazoqxDcV+HQ8e3Gct31BsyQQoyiTk0lJC8bpwTk9KOq99itC5ov3PhOu/qA+c9Z5iMBJjToqDGE1O0FmAOfHhxzW5Pfp0gs2M1YYESgQ27e/HUG9aC/1QyEeoQdN9YYVExgOZg89F0zIJI2DHwYZtDepzy0G+44BvSfx1KiIaiysp7f3+Y2Sna9vX+g5VXQyl+N/XVX1ZsD6w7VCtT2lipUV56dGAFPhTTaUEh2qbNEvf3Z/H5374SWrcjai9ECLqPh20OGeRDCdNQgnwd9oaWIULsf5/x7nIc5uC2fVVluuXV81WKySsX8HwCVpvV2FBLpzy/Uz4C59z4jD1YpRLkO6NPymNB9nu/B8ZkggibsQvhP6vxtC/A5BXm1CtnfyqgqD2ZNDTzOpRwDSUobFh/sNx2Zf3/sOcraLAKezT9ZoO+DtcSzEfeY8HnqesWpIHqA/yO/QOR6nZOkILb7K47my9ALtkp2YcSNjAG+W7sKK8QU+tFJ3qDW0J9VYWC77EoqN289Sc57FkIW4qQzQs0p5PoHcyzyavR0Ge+qST5zMqKDnK/m9l0TGFSm8OckEAEH4qPaUO/YU2DoJ7mNXmp2V5w3baGYhiAzFFeRcJCJfV2KIJ9KOFlBx03+7iKFMOGS3v1apitO3DheIjJK8i8qG9Br9qcSTmHAz9lKx9y74RFeakAEvP7UKIJFLOm4e/wt4VrwEvMAKBnwBA2HCK47Sgtv4aHUE6UVtDzCVjh2S0y2Iad8g1MMpHwnZWpB8Vr8ip+931IDQ5UYTiuY2kowahc7nodGkrIoODnF4g68zE9HKbJpl6cvw9F13acx4Ls80ERQG6CJsJhQtS0iBPwN4eEmUP1Pu52va8+55h8BAuNoHqzBod1Okm+AjFsduzQIj1O+QKaUgkkKNwP4j0WrGU6Jw3mNgeZ6JqSCcc6mWJZ+v1rqlfNNvxaHoztQ9Fmh373UTGt+GNNz+dgroDWjCVQWENpYKytV7yz/c6eQaSS3nUXQHhUjJ7ut17CqKpD2Ns++N0ztg/FbfJy1B3i1PPdHNJcp7sOp01dlWESkMG/112nabCy2xzkxwp5wIIGSb+8poi2KM7Q0EWCWtrubqdJKswUo35zmNNW9Zve3qakOZrR3Saz2ar43U9DSRAZhWyQDyq8Xn8NJRshbLg1kzSulA81LcoTBkr0nxXbZtJ6/EPgAdeaK0O9A7kCWtMpR/tqAQuUElF7eemD3x+WWNtV7Om1XgHqXigY9iAA5hmU0+RlfhBMO43q+L3HQ//fNMMJMltEfdujKXKt8s0h9XY5+zlrz/S8D2HY4OtoR8DA7rvTgeFlYvr3vb3m18OGmXH0pEANxfTCKdlm87b4mkDJeQfbppQ5+ARwD2bB9457spLNF5CWAsVUvqmtprrMEXHQjpvbZDZnWmltTSnHXmPuvH71qgH91Y177bUhpWDaZXogIMAEcO6JZjo+kM2jpUlqKCxQGhelhej3lduZ6X7ICgGDczFdbb9QTAtbOV4oWIEA7nN1fDUU7V/jDEgbCNwPeVFDCb5RTXVHXocSYaW8u91qYGxJJyt6wZbZxGB9mmaTOn5+HNMrbQGnmSvMSe227+szU+N74+V9OJAN1lBUm/X7NpNKRhLqYdFW7us4mLNMXkTmvIOGwdXUZlPLHBqKq+3qtxnNu7IdrU1Jx15jdl0Bv8dd5jf/8Jq58QGYFjbq/WCyIPRnCzhuykgcPrHdKNQH8wW0SoHCW680MGrhV9BCRvfMOcrmcmaTl7mMYr0i8J3ggH+Ir/6faeuVIJNXsf7Aqo2x/1H3VPJ/H0rxu1uzUr+jJe0dDMMGI2c7vGnqofbvC/lb897fGoTpXTGmcgHvYJsLHWCtVOOs29GPfhpK8L3q/m5ql1eLFEgnyTJ5laih6MI8zFHt90wK4dUybJNXU3kmL70PTj5ktG/7/TCFDTtC4A03Rt9gHs3phHGzWiEEBnKFolOeBUrjoiK5ggZc90DnnG34aSjembsjusTHHFJch2Jur582ERaW6BgUylyHYmsoht9kyqqn7e8zm4iCQ1Kt/5sNGkopK4eDFhj6BSjYg3dWN0WE29TDTF76vTOYdw52YaYY5WA2ayhauTnlQymmJX1e0wA4hZExQjDA/j+QK6AplbSc8ibfQYCG4uxb/2tk0n6CnNyqHa1NqQhOeUObtX4b2ZL2bb8ffj4QxT/d+oLneH8uj+Z0Eukkee7JXEFACKA1k/Jtc7VggVIiHc3WRdKvkXsG6B6AHQ62MB+KjwnEb7uXZIgPRQ3/pUZ5OaKAPOYy4SnDhO13MPz+sJXym/b0eY5b+azPtGFPNfUQNqe9G3j6aX6KY6d0eNpoqtckJIQQdlv2aoIwioYS9AZKqz5/DUW/r0y/SW0j0m/YTiTIKX/XJ2Zbfijfe9UyW/m1ORugpWXzBTRJDaVUH4r+O0ymMVVX2hBQEOQTVOe2+TjlQ9eh5P2vURRUe/Q1Ze5y3MLX0lCS9to3xSOvbsG7bngSANDKGkrjM6q1CUDxYXvw5U346VNvOc7JugadwVzwLBQwb2EdZloQmskr7JZx22GDNBDA5RT0bA4JT1tNqKMOU0Sgc9T6JAK6+oOd2OlkwujwbUpaZgCvU774PUwTMJu8rERTX+nn7+0dtL9Hif8PCw/VZ+wDrkEmKCxZCGH7A9z5AKfPRQkA1d50MoEkBYUNAxkZIm/6jbrQNQqUCD4U0zxJFyim36SeyUwq4Z21O1497BYo1mdrUwr9hnIdlgIfoa/MVaW+D8fKY+gHlzba1e8MyFAmrVQy4ajzp0+txs4uy6fLTvkhwKSR1tvc1INy36KN9rHvffg4AN6HbDCChrJPzmwd0VchNuNCAaEaioptd9+0ecdD730IerWZmueGF/5tctQtP00rmINMXh3Nac8DVGyrdZLJrLWzawCD+QKSCa/2EvY+FH2ACtqewzRg6APoHi2aK8rMsDiIe8PR3fW5Z6161Je7KvV7mlIJDOQKHuGg94c7yiulNAi/wIiCQDpAQ3FoPzn35MqK8rKukbeCQduf4L2vlYC0nOcGTULmbUknPZOgICGnrndbJol8QYREaXqqRa5QsLWBIA3LDz06UeEWmN2DzudB+VDcmnqTtji1Q5rf2CnfoHzq9Bn2BVMPaErbnXXGuHYA3hn9fz76uv3dHFNewK7uQc9xRxSPadasmbz8xvaEz8wpbC+gXu0G9jj9AtqkY8/acl6B6reuBgA6WlLoGjBrKKotmZRXC3lg8SYAcjcDV5sL2izSb2asFqYGRQKFLQLdo2koJkHtqddepZ409qc+qXAPVp+6c5H93f2TlOlmlBxUTNpNc9r6vXaUl2xvMkG+GgRg9UU66fUl2m3W7ie3OViFDTenk0bfji7c3EJQmfBGtaSNGoq6Npl0wnO/O31Vznzqd7bIdRvuHYfDgmpyeWFrA2WZvOzt74tp7t93z8INjv/75Ur4ZCLhuCd1H84oFiiNzcwJI2wTk7pIuobQljHbLFdu77K/my7uL55da393rLMIiRiyTF7yu4/Rq7gLrneADRoU1IBkzfbcJgKlZYSYvORhx4aABWcZznKtzygaSiaV8DUpJcisoRS3xvHmGcgVMLLVegCDfA4mB6pe1879A3a/RNNQigIyzEntHqx6NC3Sff3twVf+JvfgnS8IOwpICSpbQ0koH4e5/YWCsCdWfhMkux2DznoHZNhwczppNC8FbZ+jrktHS9o3FDpB1qJMt98vmxe+76xRp6rn1339g8xl6vyWJsufUarJSwiB/fJeD9Lsfvzkavv7eTc+g729WcspnyDHPdKhCxR57cNe6hcnLFAicPrMcQCAS2ZNtRcSqoukC5TWJrUy1f+mMg2Ca3f22N/12Z3ui/EbgMNMXiqfKeQ0k0p66lEok1dHS8r3AQsaL3VH9aDhgTS/59xKHN3ahK7+XOBeT5lU0jOYjWvPYNqYFuNuBoWC5W+yNiZ0mSRzBXT2DNoz9h7XAsKgd8MARVPL5JHN6BrIYZtc3KZf63HtGe8P1sprTieND342X7AFofsannf0xGIbXVn77dm85fNzRy9ZGopzVp3TJkluZ69OruC/zRDg1Kp6XaYa5efKpBJGwR0ULtuXzaMplUBLU9IYjaX8M+lkwhMtltfb7GP+U89vv0sIhvnX+rMFNKfMEVdh3PbsWvzy+XUAXGu/AjSdN7d3A7DumfbmlMP0mdbGAyVQ2CnfYNzw4ePxzFffCyKyL5hpoZSa8es3lZBmqY+9azoA4Po/Fc1fgHWj3yv9MCdMG+UMMdbKMW1aly+Eb72SNwzqKl1tBmjaMK9XPuwTRjRj455ex7HiIkz/G9VpUvPar4OExajWNPIFYVxkpps13A9KR3MKJ0wdZdnnDb83mTCvf7j3JcucsLHTiizrcUX6hG3Ep45PHdNq5R9wzjjHtWdsoeAmq2soxpBUgTYZ/ume/Xa0pIs2cz/zkK11ef0CSkNxr0NJJRK+byZV56V9BmfA+dZDt9BQe3k1p81+kKCw6/5By8yTkX4hN4PSnJZOeu+NnB6Z5jZ5aT4UwGum0x3kRg1FhvA2JRMlm7x+q5my9HtW/b7vX3y843z9uWlpSmJMWxM6e4pmVn3sUeav3gGztl8NWKBEYMqoFhwy1nrxkIrtVjMvNSh/5JSD7YdMvzEG8wXkC8K+Wd3oK6itB0VzhmsmiH2GdRkFEWx6AIozHZMPRbXdNHD0yd931KQRHmFW3MbCWKWnPpMZw+zctBLHtFmzalOkl9JimlNJj4+iP6ti87029JwUKMkEeTQy9+DkfgCdfi1/H0pHs/UAd8v+Ur+n1WdGDRTNTU2phFFry+aKK57dg1VfNm8f89VQlMnLoKFkXAJF11CakgmjWcmqS5uMmO4dTSi4/RHZnGVqbU6by+/RfXeuvu7PFtCcTiCTSmLAIIyU9mNpCq7rny9q5G5BpfpVDcBuQae30zQRsu67BJpSiZKd8m9p1glTiPgZM8fjmMkddt36xKAgBMa2NaF7IIeBXB5/XLoFb2rm9WmjW9HRnMJbO7tLalMlsEApkTb5AKtB9rnVuwAA377wGHsWqt+wShCM9zF5/GX5NgDAidNGeWZeA7k8mpIJTB/Xhje27ffkFaK40NJk2xdC2A+0yQSQCoi22d0ziKZkAmPbM16TV4CWoXCY7nThItvRn8vj+H9/HF+9/1XPeSo0e7/Bj+LQUFwDzkCugIw0ifRlvUIjk7Jmt+6B2f0zfva0MwxcNymYzIPKEd/RIicb8nw16RjZkvaduaq+bc+kfBzc1qy7KZXwbD44kM3b0UXua7Flr2V2Gy370uRoVsJIDVJKQKcShEw66bvZ4UCuYJuHTD6Dvmze9u25fShZXUMxCNnuAV1D8WpdLemkJYxM/pecpTm5Q2mV+TXjo6EoYasmBO61KI43gfqZvNREpgynvEJ/hp943XrhW1MqgQuPn2zXs1+bZPUO5DFaTr52dQ/ic795Bat2FIVHR0sah4xtw3bX1lDVhAVKiagHad1ua2ahfA1NqYRtmtAdympTwLHtGVxxxqFoSibsh//RpVuxZpdVTjZfcAx2u7oH8L/PrMFgvoAjJo7Aht1Os9P63T14ce1uezbpNmkBwEU/fc5+KP1m7ClDiC1gOZfHj8jY+2I5tIy8KtNfoOiDsL4ORw3I+/uy2N+fw/0yMgsoznZHt6Y9ZbjrzqS8Zg0lNNqakh7b/aAUNpmUd78mtxlj2eZ9jv/VtW5KmgMBrvzNywCAMXLwViaz3TJyb/KoZuMACABX3bMEADCiOeVr8konCRmDOaU/W8AIOQi6Z8bqzYtqwHnqjZ3OvIN5jGuX7ZX9rK6nfU/5tLl7IIexMq9Jy7jl6bdQEFZ/6YJMCIHewTwyqSSaU2aTl+6/ck90+qQGmkklfU1eTamE5ZQ3BLTYWpXBhwZY1wCAJ1hAtw6YNM3+rPWbTJpxKaj7eX9/Fk/KN8Q2pRJol9aN7oGcoy3JBGGsvL7XPuTdQyyTSmBkS9p3TVc1YIFSIu1SaFx1zxL8/hVrMFRaS3M6iXHtTVi5rah2Llq3B4A1Uxzd2oTBfMEe0NRABFiRSfqD8qV7l9jHOlrSntn6eTf+VUaukGXvNtzISzcVB0a3wOkbzKE5nbQEimGQXLWjG+NHZGxHtcm27XZe66zR1Gyn+UtIf4E3j2pHkMlLPdCtTd4ZvTI9tGZSnrYN5PJSoDjNirl8Af/12BsAgEtnH2yn6zP+S36+AAAwcWQG+/qynr5eI80WY+Qgq4TZLrnd/CFj25ArCOMAquhoSRujmh5duhUb9/RZ5hSDyWtMWxMS5F34plCD6I1PvOlI783mMX6EpTWrSY9qX7PyU/iYlQZzBTvIwK356LQ0JW3TKWDNovuyeUwb3YK2TArdhjY7BIrrvtywuxeTR7X4tm3Tnl6MH5HxPA+qX/3MwwOa1gR4NRS1UBBwalB2/lzBNnnpguzVjXvx/Fu7POcr3FqlaufWvcUdi1uk492quyhQ2jMpfOLdMzCmzboOL2/Y6ymfiNCeSfneG9WABUqJqN2GAeBL91rmmo+ceoidduK0UbZA6RvM4xu/XwbAsmUre/bevqITTZEgtQjNumGPnzoSgCWs2jNJwwBZsD/9tAydfa469/VlMaoljbRhoNrTM4hlm/fhjJnj7IdMn9WrBzJo5vPR2xba33Wnb74g/COeCk6Tl/tB2N09gK//zurPkS1pxww2ly/YkWttTUmP30ftcus2l2zdV3x4P3rKIfjq+UcAMK/EPnbySPQO5h15dJS2oOre3TOIdJIwXfrfdvd4r7tiYkez5xq/uLbTbnvaqKFYPpT2TMp377P3HjEBAHDouDZHeu9gHm2ZFKaPbbXt7kojbM+kMHV0Czbt6fPMyFUblTDy+Bvk+Uq71QWOGpgndDTjoJEZbO8a8IRKdw/kbHOZ+57euq8PU0e3IONj8trbm8XEjozVVwYHt9/uuwPZoqlU/w2KHVJbbzcIwa7+LNbu6rGDAfR2XfTT5/DP//eip52Kl+RkU5EvCBuzwUsAABo9SURBVLzj+idw14J1dloyQWjPSI29P2c74H/zqXdiZGvannx1+txbM8a3Ye2uHuwJuPfihAVKiRwsI3l0lO3VOt6GDZ29EELg0WVbi+ljW+2FRnt6sp4HiYjkzMu6IZWD/5Vvn4f2TBp92bxRC+kdyBlVbbcT/0/LttnfX9u8Dy+t24ORLWmMbm1yLMYrFAQ+fMvzAIAjJ3WgOeWMzRdC2LOw/T6hvS+s2W1/P2Rsq32+aqMytbhR9mdlOnI/vHctWG9/H9WaRkEUBwc1cLU0JdDalEJfNm8fW7y+E29u77JNXvrsVu+nVJLsQce0p9MMOShv6Oz1HNPbrQbdW55+C9m8sH9vZ7fzodbt4e2ZFAZyheIWPAVhDwIkJxtus5YyAXW0eNftvGP6aMyeMQYzxrVhbFsTTnnbWPtYfzaPwVwBHc1pTOhoxt6+LLL5Ar7/ZytApDmdwFGTOjCYL+CtHT2OcpXQUQOZ23z46bsXAwDmnnqIx5e1u8cSKOPaMxjZYkXy9WrXoj+bx6Y9fZg0sgWAc2HoQC6P/f05jG/PGM1lg7kCVu3oBhFJ02TBUS4AtDebF/p19ecwIpOyo97ck5Gd3QOYMCKDNsPE7hd/WwsAeGLFdjQZggEA/61P/rh0iydtZ9cAfv2iFfn1geMmASgGVuzpHcRmucfdlFFWH00d3WJbERSHjG3FLz/+DgDAKYeORa4gsLpGjnkWKCXSlknh706Y7Jhl6w/zIWNb0ZfNY2fXALbtK25w2NGcthfO/XHpFmzZ59z8MJmwBIoaNHoGcmhJJ6VvRt7ohkGuL5tHOplwCAUA9o1n4sIfPwsA2Lq/H2Pbmmxbv0pTfp3pY9tsu/5amaYG7jFtTdaAYGjTHG131FGtTcVdb6XAGD+i2dgu1Y+TR1nH97s0ID1STjmb98nB8Dt/XGGnq/NUWz98ywKs292LplTCM7vd77JJp6Qgn7fCcorqQQPKJLbRIFDampI4+yhLG+hx+W8mdli/Z+G6Tjvttc37cPy//wUAbC0DKAqjL9zzCr4ozZ5/+9qZaEknPX09IJ3BI5rTnglEz0De3sjU/SZCpSmMb89gwogM1u7qsc12gDW5Udfsy/ctgY4SKB3NaaST5DF5Pb3S8tX0ys0L9XqVCXBce5Nj1q1QfoMTDx4FwOmje2GN1XfjRmQwojmNnsG8Y/BWs/pHl25FKkkO4avaoKK43IK5qz+LDjm5ArwWhJ1dmobiEihqHdg/yyhP1Sbd5/naln24+ndLbb8WYAm1p1buwOET2z0CQXHjP50IAJgshcfyLfuxaU8fWtJJW6A3p5M456iJjnxnHjHB1kyV4PnvP5vfxxM3LFDK4D2Hj7cfDsC5U+3BYy0NZkNnLzbvtQb1N75zAYDiIPizp9+y7fKKJBGam5Lo7s+hbzCP7oGcbTtVg40a/PQH6R9OmoJt+/vxp2Xb8IdXrDdHbtjdi/ff/DcAwKNfOA1Xnvk2Tz7A8qOMaWuyZ44AsEkbLGeMa7N/56NLLW1rmfTLTJe/070rsD4b+8VlszCmNW2r42qQ1rW8plQxSGFvXxatTUnb5PWSNgDv6OrHf/7pDft/5RT+49IteOTVLbhvkeXPGt3aZAdOuMN/uwfyDi1Q1anI5gs4cao1mD26bCt+/eJ6O2hg2pgWTB7VgpZ0Eq9vLUbcPfmG9Zu+dO7haJYRSL2DeQgZWvup02fgqEnW/aGEHgB8S3Oi/uqT77SvcfdADoWCwB+XFrXbqaNbMW1MC9bvLg76F/3kWWze24fmdAJTRrU41gotWteJFVv32+bK1nTKEaSwU17T8SMyOG7KSOzsGvCYL6eNtq6RPjEpFAT+RZpwRrelPWtJnpDXFwDmvONgtKSLixff2LbfNhFP7Gi2hX63tsXOv/7a8imedpi1kHi9HJQ37+3D3NstE+q49gxGt0nTsfT9bOzsxQ3SDzayJY0RzU4T4F5XpKU7YKOrP4cRzSm0NFnXTzcPFQoCK7d1YYIUZG5hs7GzFyMyKVx9wZFo0hZrnvHfT9nn/N1PnsM9L23EI69usc2WP5y3Ehs7+9Ddn8PSa8/HLy6b5Sh3XHuT7fOZMqoFY9qa8KP5b2Ljnl5MHd3i2OvsxGmjHHnfNr5o3pw62hIoC9d1Yus+/0lmXAxpgUJEFxDRSiJaTURX16re0w8f5/hfhfUBwCFysPzfv67B5r39OH7qSPvBnqYNpMoO/4+zpgKwzBpnzByPXEHgr6t2onsgbw8yaoA9/ftP4bXN+2yN4qvnH4EvnXu4XebTK60Zngo5BIAjD+qwtalnZYizima58Z9OxNj2jGPQUA/xbXNnoaUpiYtPngYAOE76dNS7GY6WsfH6oP/n17bi0G/8CQDw+bMOwzlHT0QykcDyLfvRn83j3+Rsf/yIDA4Z24q2pqS9Sh0AFry1G72DeXvW9/jy7dgm++mcHzxj17PmP99vz9C+/dByLNlYdEiObmvy1ehe37ofTakkFq7rxM3zV2EwV3DM7I+YOALHTR2Jkw4ehTW7uvH/fl8c9A8d145kgnDc1JG4c8F6/Mcjy3H3gnX4xC8X2f0MWK+R7RnIYcu+fgxI57W+YZ8SQNs1P8zUUS12hGD3QM6zESBgaYvq2uTyBbwqBfvEjmZMGdVs9xMAXCwnK+o+aWlyajfrpLZp7Sxg3ZPKaX/HxyxTybvkoH6qZip7ZeMe2w80ti3jESifvKu4t9jBY1tlvdZvufWZNfaxtkzKvgfP+eFfsWhdJ17UzKQXnTgZ6STZ11Wf7U8a2WxPzJRWftU9r9j+ll9/8p0Y25bBXhk8USgIfOwOSxi9Y7r18qtezaR10xNvYsGa3bY2NLYtg86e4j2xakc3BvMFHDK2DdPHtmLdLqd2umLrfpx0yGgQEaaObrGvkd9775Zvsa7b7xZbkz+1OHXqmBbHeUorUaQShP5sAfNWbLfNsgp1f336PYfiK+cebj+zAOyxp1YMWYFCREkAPwXwPgBHA7iUiI6uRd0TRjTjK9pArjNVzuzmrdiOv765E7u0CBElIHRmTR8DwHrT20lS1f/03YvxyKtb7BvnsAnt9vkX/vhZ3PCYtdp+5oR2EBGe/rf32seff2uXPUCfMHUkkgnCGYePBwB8/I6XsGVvH3oGcvjC2TNxzOSRmNiRQc9gHr+TM/Glm/eiOZ2w86jB+eb5qxxrYa5+31EY29aE51bvQs9ADv3ZvCM09VNnHAoAkBYk3K/tynzSwaPwzFfPxE1zTgIAvL61C3e/sB4rtJm/UtX/5y+Wqq5HuSUSZAsUwOlbOeKgEWiTGsre3kHcvWCdo7//+qbVxh/OexOX3/mSLVDuueIUe9Z34rRR9qp5vc1A0Rl9x3Pr8K2HltvHj5w0AoBlXvr1ixvw+GuWz2qay+f2iV8uwr/d/yoOP8g6f0xbE8a2Z3CQ3MX6sWXbsE/OvFMJwn2fPhWA3L8qV8CWvX22MAGAf3rHNIxtz2B/fw5/W+UMDf7mhUcBsCYQyq+3fMs+fPk+S7BP6MjYmshzq60BXWnFgGXDf3b1LvQN5vHEiu14Q4teHN2WxqHj2vDSuk4ZDly8PkootTal0DNg+bKec0U7qWsEAL9ZuMGeqBx50Ai0NqUwaWQLXt9qTUR+8tQqAFbE2hEHjbC10+37+9GfzdsRTp8+41AcO2UkxrY3QQhr5+eX1nXamowypd3+nOX3yOULuOkJq2xlyhrdlrYFVaEg8PXfLQUAnHPURLxtfDs2y+cHsDYjXb5lv60FHD5xBHb3DGL7/n7Mnj7Gvod1fvPiBhz1rT/b2/PcNOdE+bs7HOe5fUS3fOTt9vcTXBrJJSdPw+fOPAxXnT0Tnz97ph1coHjkc6fhz1883fZNVRPvCDd0mA1gtRBiDQAQ0T0ALgKwIjBXTFx55mH4wbw38b5jD3KkN6USGN2atrcxP0IOHIrb5s7C5douse879iAcP3UkZk4YgWSC0NGcsgdPFX1z2IR2fPMDR+G7ctfiPyyxnHlqcJs+rg2TRzbjD0u24A9LtmDyyGbr/yvfDQB42/iiQFIv35kqb/aPv2sGbnpiFb5y/6v4itQgmmTEClBcd7OjawAX3GSZ0b5/8fFoz6QwfkQGj722DY+9VnT4A8D9nznVDlT4/odPwOPL/2IPvj+45ASccqg161Wz34/cVoyE+eYHrEHw/y6bhfff/Dc8sHiTvYswAFxysqXRjW1zRoqNa8/gma++F22ZlL2rwd//7HnHOR9/93QM5gq20/Nvq3bhb6usge6dM8bY57W4ZnUffvtUfO7MwwAUHe8691xxiq0FKv/EddK8pV4J+8jnTsMHf2L5rtTv+eAJk/HjSy2hqu6TH81fhR/Ntwa5uz4xG7Nlu5RJ6l03PImzjpxg1z2xo9muW4+su0pOGABrkP7bql2Ycc2fHO2e2NFsO3wVs7RX2L77sHF4dNlWHPXtPzvOmT19DMa2ZfD3J03B1Q8u85R7pmzf+BEZzFuxHW/7RvH42v96PwDgbdok6cGXN9vf/0+afjZ09mJDZy+O/Fax7levPQ/pZMK+nz9620J86ETLOnDp7IPxbzJCTw3k77j+CTvvZ9/7NntAff6t3fjmH5bh4SVFp/jXZN58wfLl/MPPnsP4ERks2bgXqQTh8IntWLPTqvc/HlmOWdPH4GsPWMLmPXLypa71O/9zPgDggmMOwrc/eDS6+nO4+OSpmH71o451V/dccYpDkDzx5fdgf38W//Cz5/G+Yyc5+vTkQ8bgk6fNwKub9uIb7z/KcaylKWn/dhPKulALhqyGAmAKgI3a/5tkWk1IJAhL//083CwHBJ0nv/Je+/t1Fx3rOHb2URPx8XdPBwD8/CNvx4jmNI48qMM28zzx5ffY5z76hdPt7588/VA8pWkigHNWc+Sk4vct+/px3NSRDjvrc1efhYkdxUFYPdAjW9P4wlmHOcq96pyZ9vdkgvD+45xC82w5YFx5pjMfYGlF75heHJxHtqZxzlHFAVCPNmrPpByD4/FTR+KTp1uazdGTOxybHwLAzz9ysr230YQRGXzytBn2sb87YbJtNprhCpEFgBXXnY9rP3gMvvuhY3GF1J509L76h7dPdRy77qJjbGf9J7Q61e9VAhIAvnr+kY7jarA/bupIfNml1er9YtJej9Kuqd4m5bxeff37AACnHDoGbg7V7Oj/9I6DPcdf/ta5AKxNNud/5T04YepILLjmLEc/XDp7msM/CAAXHj8J933mVOu+ON456AFWAIHi5IOd71f//oePt8s3hY7//CNvtzW6S2dPcxz74jkzbfPNhBHFvGpy9bmzDrMnQe8+bByOmFicyE0YkcHXL7Cuy1VnW/f2r17YgP39VojyHR97B847xrrH1fuOXt6wF48v345x7U14+HOngYhsgXHfok22MPnhP56A82XeE6eNcrTt3KMn4vxjDsLFchJ09+WzAVibzf7pC6c77hvAmji+/eDRWHDNWfjC2TPh5psXHo37P/OumpuxSoHCtiBvVIjoEgDnCyE+Kf//KIDZQojPu867AsAVAHDwwQefvH79ek9Z1eDlDXuwY38/LjjW+9AVCgKD+UJZN8be3kG8uslaI6I//P3ZPG57di329Axi8YY9uHnOSR5zCwD86oX1OHziCHvmq8jJF1Nl88Jh81f0DeaxbncPjjxohKfeax9ajl3dA7jijEMxa/oYzw7IQgjc89JGHDah3SFsACvcc96K7Xjsta344jmHO8x7fYN5rN7RjSff2IETpo20I1fc/fHCmt047+iD7J2gVbnpJOG51bvxjhmj7X2cdD5+x0L0ZfO45V9OtleU67/r+kdfx8mHjMaHTvLOU15csxuL1u/BJbOmYoIWtSaEwNpdPVi6aR8yqQTed5z3+ufyBSzZuNc2d+pt7uwZxH2LNuLcoyc6BIri2w+9hrsWrMevP/lOvPuwoi+vs2cQz7y5A//xyAp84LhJ+PYHj3b85tU7uvDgy5tx2sxxGNeeweETR3jKNpHNF7B4/R70ZfN4Yc1ufPW8I2zhClhRaUs37cPyLfvw0VMPcdRZKAis3N6FghDoaE577scVW/ZjV/cANnT2YtqYVnumr+jqz6KrPycXQ7Z67ssd+/vx1Mod6BvM42Pvdgp6AHhrZzdeWLMbHzxhsq0x7+oewE+fWo1UgvCB4ydj8shmTOgoXr/+bB7//vByvOuwcUgQcPrM8Y53jCxevwfzVmxH32AOs2eMxQdcQvXN7V1Ytb0bh45v8zwrqk8SCR8HSwNDRIuFELNCzxvCAuVUAP8uhDhf/n8NAAgh/ssvz6xZs8SiRYv8DjMMwzAGogqUoWzyegnATCKaQURNAOYAeLjObWIYhjlgGbJOeSFEjog+B+BxAEkAtwshlodkYxiGYarEkBUoACCE+BOAP4WeyDAMw1SdoWzyYhiGYRoIFigMwzBMLLBAYRiGYWKBBQrDMAwTCyxQGIZhmFgYsgsby4GIugBsA7Av4LSRAccPBrChzLxBxyrNG9SuatbLfRVPXu6r6Me5r6Ifj7OvjhBChG+vIIQ4YP4ALAJwa8g5vscB7KwgbyX1huX1bVeV6+W+4r7ivjoA+grAoqCy1N+BaPJ6pILjewOOheWtpN6wvEHtqma93Ffx5OW+in6c+yr68Wr2lZEDzeS1SETYj6Za+atFI7arEdsENGa7GrFNQGO2qxHbBDRmu+JsU9SyDjQN5dY6568WjdiuRmwT0JjtasQ2AY3ZrkZsE9CY7YqzTZHKOqA0FIZhGKZ6HGgaCsMwDFMlDniBQkS3E9EOInpNSzuBiBYQ0TIieoSIOmR6mojulOmvq3ewyGNPE9FKIloi/7xvg6pOm5qI6A6Z/ioRvVfLc7JMX01EN5P7bT/1a1ecfTWNiJ6S12M5EV0l08cQ0TwiWiU/R2t5rpF9spKIztfSY+mvmNtUt74iorHy/G4i+omrrLr0VUib6tlX5xLRYtkni4norAboq6A2xdZXDqKEgg3nPwBnAHg7gNe0tJcAvEd+/wSA78jv/wzgHvm9FcA6ANPl/08DmFWHNl0J4A75fQKAxQAS8v+FAE4FQAAeA/C+BmlXnH01CcDb5fcRAN4EcDSA7wO4WqZfDeB78vvRAF4FkAEwA8BbAJJx9lfMbapnX7UBOA3AZwD8xFVWvfoqqE317KuTAEyW348FsLkB+iqoTbH1laONcRc4FP8ATIdzkNyPon9pGoAV8vulsMLpUgDGygs6phoXqIQ2/RTAR7Tz5gOYLW++N7T0SwH8b73bVY2+crXvIQDnAlgJYJJMmwRgpfx+DYBrtPMflw97VfqrkjbVu6+08z4GbfCuZ1/5talR+kqmE4DdsCYIde8rd5uq2VcHvMnLh9cA/J38fgmsgRIAHgDQA2ArrBWo/yOE6NTy3SHVx2+Vq9aW0aZXAVxERCkimgHgZHlsCoBNWv5NMi1uSm2XIva+IqLpsGZlLwKYKITYCgDyU6n0UwBs1LKpfqlKf1XYJkW9+sqPevZVGI3QVx8G8IoQYgCN01d6mxSx9xULFDOfAHAlES2GpVoOyvTZAPIAJsMyTXyFiA6Vx/5FCHEcgNPl30dr1KbbYd2kiwDcBOB5ADlYMxI31QjpK7VdQBX6iojaAfwOwBeFEPuDTjWkiYD0erYJqG9f+RZhSKtVXwVR974iomMAfA/Ap1WS4bSa9pWhTUCVxisWKAaEEG8IIc4TQpwM4LewbNqA5UP5sxAiK4TYAeA5ALNkns3yswvAb2AJn6q3SQiRE0J8SQhxohDiIgCjAKyCNZhP1YqYCmBLnG0qs12x9xURpWE9YL8WQjwok7cT0SR5fBKAHTJ9E5yakuqXWPsrpjbVu6/8qGdf+VLvviKiqQB+D+AyIYQaM+raVz5tqtp4xQLFgIp4IKIEgG8C+Lk8tAHAWWTRBuAUAG9Is844mScN4EJYpqCqt4mIWmVbQETnAsgJIVZI1beLiE6R6uxlsGyusVJqu+LuK/nbbgPwuhDih9qhhwHMld/novjbHwYwh4gy0hQ3E8DCOPsrrjY1QF8ZqXNf+ZVT174iolEAHoXlC3tOnVzPvvJrU1XHq7idMkPtD9aseiuALKzZxOUAroLlcH8TwA0oOp3bAdwPYDmAFQC+KtPbYEUxLZXHfgQZpVODNk2H5ZR7HcATAA7Rypklb5S3APxE5alnu6rQV6fBMiEsBbBE/r0fVtDEfFha0XzI4AmZ5//JPlkJLeImrv6Kq00N0lfrAHQC6JbX/OgG6CtPm+rdV7AmUz3auUsATKhnX/m1Ke6+0v94pTzDMAwTC2zyYhiGYWKBBQrDMAwTCyxQGIZhmFhggcIwDMPEAgsUhmEYJhZYoDBMg0BEnyGiy0o4fzppOz8zTL1J1bsBDMNYi82EED8PP5NhGhcWKAwTE3LDvj/D2rDvJFiLPS8DcBSAH8JaGLsLwMeEEFuJ6GlYe5y9G8DDRDQCQLcQ4n+I6ERYuw60wloQ9wkhxB4iOhnWPmm9AJ6t3a9jmHDY5MUw8XIEgFuFEMfD2tr/SgA/BnCxsPY7ux3A9dr5o4QQ7xFC/MBVzl0Avi7LWQbgWpl+B4AvCCFOreaPYJhyYA2FYeJloyjum/QrAN+A9XKjeXKH8CSs7WsU97oLIKKRsATNMzLpTgD3G9LvBvC++H8Cw5QHCxSGiRf3XkZdAJYHaBQ9JZRNhvIZpmFgkxfDxMvBRKSEx6UAXgAwXqURUVq+n8IXIcQ+AHuI6HSZ9FEAzwgh9gLYR0SnyfR/ib/5DFM+rKEwTLy8DmAuEf0vrN1ffwzrlb43S5NVCtYLx5aHlDMXwM+JqBXAGgAfl+kfB3A7EfXKchmmYeDdhhkmJmSU1x+FEMfWuSkMUxfY5MUwDMPEAmsoDMMwTCywhsIwDMPEAgsUhmEYJhZYoDAMwzCxwAKFYRiGiQUWKAzDMEwssEBh/v9GwSgYBaOAKgAAi8TzUCE/aTAAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'].plot()" ] @@ -215,9 +2249,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 14, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXmYXGWZ6H9vbb13p/fsZAWyAIGEhFVFloA6AypoGBUcuYM6OOPoeOfKXGf0yjAjMyrKjOIgIOggi7iQURADyLBnIwkhISF7Oumkl/S+1P7dP8453dXd1d3V3VVdS7+/56mnqr8636nvnK4673l3McagKIqiKIngSvcCFEVRlOxBhYaiKIqSMCo0FEVRlIRRoaEoiqIkjAoNRVEUJWFUaCiKoigJo0JDURRFSRgVGoqiKErCqNBQFEVREsaT7gUkm6qqKjNv3rx0L0NRFCWr2Lp1a7Mxpnq07XJOaMybN48tW7akexmKoihZhYgcSWQ7NU8piqIoCaNCQ1EURUkYFRqKoihKwqjQUBRFURJGhYaiKIqSMKMKDRHJF5FNIrJDRHaJyP+zx78hIsdFZLv9+EDMnNtFZL+I7BWRtTHjK0Vkp/3ePSIi9nieiDxuj28UkXkxc24WkX324+ZkHryiKIoyNhIJuQ0A7zfGdImIF3hFRJ6x37vbGPPt2I1FZCmwDlgGzASeE5HTjTER4F7gVuAN4GngauAZ4Bag1RizSETWAXcBHxeRCuDrwCrAAFtFZL0xpnVih60oiqKMh1E1DWPRZf/ptR8j9Yi9FnjMGBMwxhwC9gOrRWQGUGqMed1YPWZ/ClwXM+dh+/WTwOW2FrIW2GCMabEFxQYsQaMoipLRhCNRHt98lHAkmu6lJJWEfBoi4haR7UAj1kV8o/3WF0TkLRF5UETK7bFZQF3M9GP22Cz79eDxAXOMMWGgHagcYV+D13eriGwRkS1NTU2JHJKiKEpK2XKklf/zy538cW9uXZMSEhrGmIgxZgUwG0trWI5laloIrABOAN+xN5d4uxhhfLxzYtd3nzFmlTFmVXX1qFnwiqIoKac3GAHg3YbONK8kuYwpesoY0wa8CFxtjGmwhUkU+DGw2t7sGDAnZtpsoN4enx1nfMAcEfEAZUDLCPtSFEXJaAJhyyw15YSGiFSLyDT7dQFwBbDH9lE4fBh42369HlhnR0TNBxYDm4wxJ4BOEbnA9lfcBDwVM8eJjLoeeMH2ezwLXCUi5bb56yp7TFEUJaMJ2r6MvSdzS2gkEj01A3hYRNxYQuYJY8xvReRnIrICy1x0GPgsgDFml4g8AewGwsBtduQUwOeBh4ACrKgpJwrrAeBnIrIfS8NYZ++rRUTuADbb233TGNMygeNVFEWZFIK2pnGwqZtwJIrHnRtpcaMKDWPMW8C5ccY/NcKcO4E744xvAZbHGfcDNwyzrweBB0dbp6IoSibhCI1gJMqRlh4WVheneUXJITdEn6IoSoYRDEf6Xr+bQyYqFRqKoigpIBiTn/FuQ9cIW2YXKjQURVFSgGOemlGWn1MRVCo0FEVRUoAjNBbXllDX2pPm1SQPFRqKoigpIBCJ4nO7KM5z0xOMjD4hS1ChoSiKkgKC4Sg+j4sCr6cvOzwXUKGhKIqSAvqEhs9Fb0iFhqIoijICwbBlnir0eegJhtO9nKShQkNRFCUFBCOOecqNPxQlGh2po0T2oEJDURQlBTjmqUKfGyBnTFQqNBRFUVJAv3nKEhq5EkGlQkNRFCUF9JmnfFaJv1yJoFKhoSiKkgIC4X6fBqh5SlEURRmBYDhKnifWPJUbEVQqNBRFUVKA49MocBzhap5SFEVRhsPxaagjXFEURRmVwSG3PerTUBRFUYaj3zzlRE+pT0NRFEUZhtiMcJhC5ikRyReRTSKyQ0R2icj/s8crRGSDiOyzn8tj5twuIvtFZK+IrI0ZXykiO+337hERscfzRORxe3yjiMyLmXOz/Rn7ROTmZB68oihKqpjKGeEB4P3GmHOAFcDVInIB8FXgeWPMYuB5+29EZCmwDlgGXA38UETc9r7uBW4FFtuPq+3xW4BWY8wi4G7gLntfFcDXgTXAauDrscJJURQlU3E0jTyPC5EpFD1lLJwGt177YYBrgYft8YeB6+zX1wKPGWMCxphDwH5gtYjMAEqNMa8bYwzw00FznH09CVxuayFrgQ3GmBZjTCuwgX5BoyiKkpEYY6w8DbcLEaHQmzuNmBLyaYiIW0S2A41YF/GNQK0x5gSA/Vxjbz4LqIuZfswem2W/Hjw+YI4xJgy0A5Uj7EtRFCVjCUWsirY+j3WJLfB5ppbQMMZEjDErgNlYWsPyETaXeLsYYXy8c/o/UORWEdkiIluamppGWJqiKErqCUas/uCO0Cj0uadm9JQxpg14EctE1GCbnLCfG+3NjgFzYqbNBurt8dlxxgfMEREPUAa0jLCvweu6zxizyhizqrq6eiyHpCiKknSCYVtouG1NYyqZp0SkWkSm2a8LgCuAPcB6wIlmuhl4yn69HlhnR0TNx3J4b7JNWJ0icoHtr7hp0BxnX9cDL9h+j2eBq0Sk3HaAX2WPKYqiZCx9QsNjxQAV+Nw5Ez3lSWCbGcDDdgSUC3jCGPNbEXkdeEJEbgGOAjcAGGN2icgTwG4gDNxmjHHO1ueBh4AC4Bn7AfAA8DMR2Y+lYayz99UiIncAm+3tvmmMaZnIASuKoqSafqERa56aIkLDGPMWcG6c8VPA5cPMuRO4M874FmCIP8QY48cWOnHeexB4cLR1KoqiZArBiCUgYoVGW08onUtKGpoRriiKkmQCg30aPk/OmKdUaCiKoiQZxzyV52gaXrf201AURVHiM9inUeCbQtFTiqIoytgYnKdR4HPjV/OUoiiKEo/BeRqFXjehiCFkC5NsRoWGoihKkolnnoLcKI+uQkNRFCXJDC0j4jRiUqGhKIqiDGJwyG1/n/Dsj6BSoaEoipJkBofc5udQ9z4VGoqiKEkmXhkRICciqFRoKIqiJJl4pdFBNQ1FURQlDkNKo6vQUBRFUYYjGI7iEvC4B0VPhdQRriiKogwiGIn2maZAzVOKoijKCATD0T7TFPRHT2mehqIoijKEQDja17UPVNNQFEVRRiAYjvblaAB43S68bsmJnhoqNBRFUZLMYJ8GQIE3N1q+qtBQFEVJIpGoIRiODPBpgBVBlQtlREbtEa4oiqIkxrajrXzi/o2UF/qoKPINeK8wRxoxjappiMgcEfmjiLwjIrtE5Iv2+DdE5LiIbLcfH4iZc7uI7BeRvSKyNmZ8pYjstN+7R0TEHs8Tkcft8Y0iMi9mzs0iss9+3JzMg1cURUkmu0900BOMcLytd4h5Kj9HzFOJaBph4G+NMW+KSAmwVUQ22O/dbYz5duzGIrIUWAcsA2YCz4nI6caYCHAvcCvwBvA0cDXwDHAL0GqMWSQi64C7gI+LSAXwdWAVYOzPXm+MaZ3YYSuKoiSfho5A3+uh5qkpomkYY04YY960X3cC7wCzRphyLfCYMSZgjDkE7AdWi8gMoNQY87oxxgA/Ba6LmfOw/fpJ4HJbC1kLbDDGtNiCYgOWoFEURck4Gjv8VBX7uHJpLctmlg54r8DnzonoqTH5NGyz0bnARuBi4AsichOwBUsbacUSKG/ETDtmj4Xs14PHsZ/rAIwxYRFpBypjx+PMURRFyShOdviZXpbPj29aNeS9Qp+bxhhNJFtJOHpKRIqBXwJ/Y4zpwDI1LQRWACeA7zibxpluRhgf75zYtd0qIltEZEtTU9OIx6HkFsYY/vwnm/jvHfXpXoqi0NARoLYkP+57hT4PPVOl9pSIeLEExiPGmF8BGGMajDERY0wU+DGw2t78GDAnZvpsoN4enx1nfMAcEfEAZUDLCPsagDHmPmPMKmPMqurq6kQOSckRGjoC/HFvE9uOtqV7KYpCQ4ef2rL4QqPAlxuO8ESipwR4AHjHGPPdmPEZMZt9GHjbfr0eWGdHRM0HFgObjDEngE4RucDe503AUzFznMio64EXbL/Hs8BVIlIuIuXAVfaYogCw+0Q7AMFI9v8YlewmEI7Q0h0cVtMo8OaGIzwRn8bFwKeAnSKy3R77e+BGEVmBZS46DHwWwBizS0SeAHZjRV7dZkdOAXweeAgowIqaesYefwD4mYjsx9Iw1tn7ahGRO4DN9nbfNMa0jO9QlVxkd30HAIFQNM0rUaY6TZ2Wv2J6WV7c9wttR7gxBjvbICsZVWgYY14hvm/h6RHm3AncGWd8C7A8zrgfuGGYfT0IPDjaOpWpye4TltBwOqUpSrpo6PADUFM6vHnKGKuYoVP1NhvRMiJKVuNoGk6nNEVJF06OxvRhhEahNzcq3arQULKWrkCYIy09gHX3pijp5GS7pWnUDic07O592V5/SoWGkrXsPdmBsQOwVdNQ0k1Dpx+f20V5oTfu+06f8GyPoFKhoWQtjmlqXmWhCg0l7TS0+6kpzRvWyV2g5ilFSS+NnQFcAnMriwioI1xJMw0dgWH9GZA73ftUaChZS28wQr7XTZ7HRSAHavoo2U1Dp6VpDEefeSrLs8JVaChZiz8cocDrxudxacitknaaOwNUFw8vNBxHeG8wu7+rKjSUrKU3aMW757ld6tNQ0kogHKHDH6ZqRKHhmKdU01CUtOAPR8j3usjzujTkVkkrp7qCAFSVJGKeym5TqgoNJWvxByMU+Nz4VNNQ0kxzl5XYVzmoxWssGj2lKGmmNxQh32P7NFRoKGkkIU1DhYaipBd/yNI08jxuAuHs/iEq2U2TrWmM5Ah3uYR8r4te9WkoSnroDUXJszWNqIGwRlApacIxT43kCAcrgkp9GoqSJhxNw+exvsYadpsa/KEIv3rzGCE9v8PS3BmkyOfuc3YPRy701BhTj3BFyST8oQgFXhd5ttAIhKIUDu+HVMZBJGr44mPbeHZXAzUl+VyyuCrdS8pImrsCI/ozHApzoHufahpK1tIbsjLCVdNIHd9/7l2e3dUAwIn23jSvJnNp7gqMGDnlUODLfk1DhYaStViahhVyC1rpNhU8904j582dBli1vpT4NHcFRvVngGWeUk1DUdJANGrwh6Lked3k2aGMmuCXfLoCYeZWFFKa7+nrTKcM5VRXMGHzVI/WnlKUyccRELGahobdJp+uQJjifA+1pfkqNIYhHInS0hNMSNMo9HnUPKUo6cBvhy3GOsLVPJV8ugJhivO8ttBQ81Q8WnqCGAPVxaP7NPK8LgKh7P6ejio0RGSOiPxRRN4RkV0i8kV7vEJENojIPvu5PGbO7SKyX0T2isjamPGVIrLTfu8esbuViEieiDxuj28UkXkxc262P2OfiNyczINXshcn1t0pjQ4qNJJNIBwhGI5SnOempjSPRtU04tLcaWeDJ+jT8E+BPI0w8LfGmCXABcBtIrIU+CrwvDFmMfC8/Tf2e+uAZcDVwA9FxAlevhe4FVhsP662x28BWo0xi4C7gbvsfVUAXwfWAKuBr8cKJ2Xq4giN2DwN9Wkkl+6AdY6L8yzzVGNngGjUpHlVmUUgHOFAUxcAlQkIjfypIDSMMSeMMW/arzuBd4BZwLXAw/ZmDwPX2a+vBR4zxgSMMYeA/cBqEZkBlBpjXjfGGOCng+Y4+3oSuNzWQtYCG4wxLcaYVmAD/YJGmcI4PzwnIxxU00g23QHLYVuc76W2JI9w1NDSE0zzqjKL6+99nb96dBvAiF37HPK9LnpDEYzJXuE7puQ+22x0LrARqDXGnABLsIhIjb3ZLOCNmGnH7LGQ/XrwuDOnzt5XWETagcrY8ThzlCmMP0bTyPNYiqzmaSSXTr8tNPLcFOdZ57ihw5+QGWaqcORUNxctrOTmi+Yxt7Jw1O0LvG6iBkIRg88Tv5d4ppOwI1xEioFfAn9jjOkYadM4Y2aE8fHOiV3brSKyRUS2NDU1jbA0JVfwh2KipzwaPZUKuhxNI89LjX0X3ajO8AH4Q1HOmTONtcumJ7R9vh0e7s/i72pCQkNEvFgC4xFjzK/s4Qbb5IT93GiPHwPmxEyfDdTb47PjjA+YIyIeoAxoGWFfAzDG3GeMWWWMWVVdXZ3IISlZjpMgle91qXkqRfSbpyyfBqBhtzGEI1GCkWhfyfNEcHKKstmvkUj0lAAPAO8YY74b89Z6wIlmuhl4KmZ8nR0RNR/L4b3JNmV1isgF9j5vGjTH2df1wAu23+NZ4CoRKbcd4FfZY8oUp88RrtFTKaMz0G+eckp+a9htP853sHCUIoWxOALGn8V9whPxaVwMfArYKSLb7bG/B74FPCEitwBHgRsAjDG7ROQJYDdW5NVtxhhHrH4eeAgoAJ6xH2AJpZ+JyH4sDWOdva8WEbkD2Gxv901jTMs4j1XJIfwxIbcaPZUauvz95imfx0VlkY+GTtU0HPq13cSFRr7X+q5ms3lqVKFhjHmF+L4FgMuHmXMncGec8S3A8jjjfmyhE+e9B4EHR1unMrUYIDTcKjRSQax5CqCmNF9zNWIYj6aR75kC5ilFyUR6B0RPqXkqFTjmqUL7Trqq2EdTl4bcOjjlQMbi03D6bWRz0UIVGjlKKBLlLx/ZytvH29O9lJTgRE/le1yICD63S0Nuk0yXP0xxngeXyzI0lOR76PKH0ryqzKGvKsFYNI0+81T2fldVaOQox1p7eXrnSV7c2zj6xllIbyiC1y14bNOUz5P9NX0yje6AJTQcinyevixxpV9bKBxL9JSap5RMxbE9N+VoDwS/3YDJIc/jIhjJ3h9iJtIVCFOU13+Oi/I8fX4OpV9ojNbiNRZnWxUaSsbRZDe6z9XGOYOFhs/jUp9GkukMhCnO9/b9XZznoTsYzuoSGMmkZzyO8KmQp6FkJ07mbq5qGr3ByAAHpM/j0uipJNMdCFMSa57K8xA1/bb8qY5/PCG3dtCGP4tNqSo0cpTc1zSifU5FsM1TKjSSSpd/oHnKCb3tUhMVAD1BO7rMl3gJv77oqSwWvCo0cpRYTSMXzQm9oaGahgqN5OI0YHJwihaqM9yiN6b+WaJonoaSsTTambu9oUhO3hn2DvZpuNU8lWy6AmFK8gdGTwHqDLfptTWNWI13NFwuKzxczVNKxhHry8hFE1VAHeEpxRgzJHrKCb/NxZuQ8eBou3YD0oTJ97pU01Ayj6bOAPPs+v656AwfbJ7K87gJaHJf0vCHokSiZoB5qihPNY1YeoKRMUVOOWR79z4VGjlIKBLlVHeQZTPLgNzUNAY7wq3kvuz9IWYaXYPqTkG/0FBNw2KwiTRRVGgoGccpuz7Q0pmlQA5rGr5B5inVNJJGV0xZdIfiPk0jey94ycQ/6DuYKAVet/o0lMzCcYIvrinG53b1/Z1L+INxMsLVp5E0YsuiOxT1RU+ppgETMU+5NORWySwczaK2NJ/qkryc1DT84aFCQ6Onkke/pjE0ekrNUxa9wfGZp/LUPKVkGo4Po7okj6ocFBrhSJRQxPTFvIMVcquaRvKIJzRcLqHQ51ZNw6Y3ND5No8Dr1iq3SmbhJPZVFedRk4NCo9M2ncTmEOR53So0kkhrj+UXm1boHTBelOdRTcNmcCmbRMn3uvpKkGQjKjRykKYuP+WFVovOmpK8nIue6rB7OpQV9F/QtJ9GcmnttoRGRZFvwHixCo0+eoLjc4Tne91Z3e5VhUYO0t4bprzQ+rGXFXhp7w3lVCmR9l5LaJTGCg2Pi0jUEFbBkRRaeoL4PK4h5peiPDVPOfhD49M0CtSnoWQanf5QX3x9UZ6HSNTklJO4o9e6aJXGmqeclq8qNJJCa3eQikLfkGxnbcTUz0SS+3K63auIPCgijSLydszYN0TkuIhstx8fiHnvdhHZLyJ7RWRtzPhKEdlpv3eP2N9GEckTkcft8Y0iMi9mzs0iss9+3Jysg851uvz9NYOKczCLt888FWNvd8wEPVn8Y8wkWrpDlA8yTYGapxyMMUOqEiRKnteV847wh4Cr44zfbYxZYT+eBhCRpcA6YJk954ci4pzVe4FbgcX2w9nnLUCrMWYRcDdwl72vCuDrwBpgNfB1ESkf8xFOQTr9YUrs+PqiHEzI6nDMUzENgqbZ5rg224GrTIzWniAVRd4h48X5ViOmqY6TnFcwhrLoDgV20EY0mp0m41GFhjHmJaAlwf1dCzxmjAkYYw4B+4HVIjIDKDXGvG4s4/pPgeti5jxsv34SuNzWQtYCG4wxLcaYVmAD8YWXMohY85ST0ZtLd4fxfBoVttBo6Q6lZU25Rmt3sM8vFou2fLVwkvMKxlDh1qGve1+WOsMn4tP4goi8ZZuvHA1gFlAXs80xe2yW/Xrw+IA5xpgw0A5UjrCvIYjIrSKyRUS2NDU1TeCQcoPOmJLWfZpGDt0ddvhDuF1CUYw92YnyaenOrUixVBGNGupaevj5xqP86H8ODAmUaOkJDomcAjVPOTgNmMYVPZXl3fvGrltZ3AvcARj7+TvAZ4B4NYLNCOOMc87AQWPuA+4DWLVqVXbqfEkiGrVKWjttOnOxyFxHb5jSfM8AJ22/0FBNYzSe3XWSLz2+fYD/58PnzqK2NB+wkifbe0PxNQ2fB38oSjgSxeOeunE0TvTTuMxTvuxuxDSu/7oxpsEYEzHGRIEfY/kcwNIG5sRsOhuot8dnxxkfMEdEPEAZljlsuH0pI9AdDGMMlNj2/lx0hLf3hgaYpqA/Ca1VfRqj8vqBUxgDd354Of/n6jOB/iKXgB2iPTRHA2LqT03xgANH4I4vuS+7W76OS2jYPgqHDwNOZNV6YJ0dETUfy+G9yRhzAugUkQtsf8VNwFMxc5zIqOuBF2y/x7PAVSJSbpu/rrLHlBFwNIoh5qkcEhod/tCAxD6wfohFPjct3So0RqO+rZc5FQV8Ys1pnDt3GjBQ2Dqvh4uegtz6Po0HJ2R2PCG3eVne8nVU3UpEHgXeB1SJyDGsiKb3icgKLHPRYeCzAMaYXSLyBLAbCAO3GWOcM/N5rEisAuAZ+wHwAPAzEdmPpWGss/fVIiJ3AJvt7b5pjEnUIT9lcUps9DnC+4rMZecXNB4dvaEBkVMOFcU+FRoJcKLdz4yyAgAq+8x6/efNMfFVDOMIBxUaPfYFfzwFC/vNUznq0zDG3Bhn+IERtr8TuDPO+BZgeZxxP3DDMPt6EHhwtDUq/fTXZbIuqoU5WM66vTfE9LL8IeMVhSo0EuFEey/LZ1m9VsrjCg1H04gTcpuDPrLx4J+ApuE4wrO1adjU9WTlKJ124pvz4/a6Xfg8rpwSGh3+8BDzFFgXQPVpjIw/FKG5K9inaUwr8CIyUGg45zC+TyP38n7Gg/o0lJzB0TRiS2zkWpjksOapQt8Ah64ylJPtVkOuGbam5nG7KCvwDhC2fZpGHPOU4ytzsvKnKs4Ff1yl0bPcPKVCI8fod4QP7LiWK5qGPxQhEI4OiZ4C685YNY2RqW/vBWDmtIK+sYoiH6diNY3uIIU+d1x7vaN9TPXz7CSYxvZQT5QC1TSUTKLPPBXzZS7yeXImRNK5wy2N82MtL/LRE4xkbVTKZHCibaCmAZaG1hrr0+iJnw0O/aHNLVNco2vqDFCc56FwHHkajnaSrTdyKjRyjE5/GBEGZEsX51Dph74Kt8NoGqB3wSNxwtY0HJ8GWMK2ZZCmEc+fAVa4aEmeh5Ypfo6bugJUl+SNa262J9yq0MgxOv1hivMGZkvnUr2gPk0jniPcvjtWv8bw1Lf7qSjyDSh/UTlIaLT0xK9w61BeNFAzmYo0dY5faOR5XHjdkrW/SRUaOUanPzzESZxLjvD2OBVuHSqLVdMYjRNtvQNMU9AfdWaMIRSJcqipi5lxQppjt2/pmdqO8OYJCA0Ryeq2uSo0coyuQGhA72xwHOG5Yed3yqLHDbktHJpzoAwkNrHPobLIRyhi6AyE2Xy4hQ5/mPedUTPsPioKvVO+MGRTZ4Dq4vEJDbD8jCo0lIzAMU/FklvmKcenMdQB2efTUKExLMfbepk5bZCmUdh/3jbsbsDncXHp4qph91FRlEfrFC4M2RuM0BkIj1vTgOz2M6rQyDE6Y7r2ORTneexChtlfADheAyaHsjiJako/Ld1BOv1h5pQXDhivsM16p7qDPPdOA5csqupz1sajosg7pc9xc5elZU1EaGSz9q9CI8foCoQpHnRBLfR5iJrsjQuPpb03hM/jiptD4HYJ5YW+KR/ZMxw7j7cDsMwuIeLg1Jh64+Ap6lp6uWJJ7Yj7KS/y0RuKZHWf64nQ2DlxoVGc76VTNQ0lE+j0D/Vp5FL3vtFsyaX5nr6wXGUgO4+1AbB8VtmAcces9+Arh/G4hCuWDu/PgJgih1NUODfZQqNmQuap7E24VaGRY3TEMU/lUr2gk+3+uMUKHUoLvFO+xMVw7DzezvyqoiGmPUdoNHcFuGHVHGpKhj+/MNAHMhVpSoZ5yqc+DSUDCIQjBMPRvq59DrlUzrqhw8/00hGERr63z++hDOTt4x1DtAywMpR9Hhcel3DbZQtH3Y8jZE5NVaHRGcAlUFk0EZ+GRk8pGUCXf2jdKcidxjnGGE52+PvaksajJN/TV7RR6edUV4Djbb2cHUdoiAhLZpRy04XzmD3ISR6PqR6l1tQZoKIoD7crXkfqxCjJtzSNbAxOGW+PcCUD6WvANJymEczui2lnIExPMML0spF8GmqeiofjBI+naQD8+vMXIQleAyvi9OCYSkwkG9yhKK8/OGU89avSiWoaOUTLMH0Q+h3h2e3TaLDLequmMXbeHiZyysHlkgGlZ0aiNN+LawqHNjd1+pMiNCA7g1NUaOQQjrlgsNDIFZ/GyQ5LaIzo0yjw0hOMEIpkZ6+CVHG0pYeakry4+S1jxTXFQ5snmg0O/Tdy2RicokIjhziV60LD1jRGjJ6yI8dU2xhIY2eAmtKJXehiqZiiRQsjUUNT18TPZZEve3+TKjRyiJbhhIYve1XhWBo6EjFPWXfSnerXGEBjR2DUUNqxUD6ocdNUobHTTyhimF1eMPrGI1CcxTc3owoNEXlQRBpF5O2YsQoR2SAi++zn8pj3bheR/SKyV0TWxoyvFJGd9nv3iG1AFZE8EXncHt8oIvNi5txsf8Y+Ebk5WQedq7R2B/F5XENaULpdQoE3e5OJHBo6ApQVeONmgzuPS2fKAAAgAElEQVQ4JdM1wW8gjZ2BCSWjDWZw46apwrFWqx/J4FIsYyWbIxoT0TQeAq4eNPZV4HljzGLgeftvRGQpsA5YZs/5oYg4v/B7gVuBxfbD2ectQKsxZhFwN3CXva8K4OvAGmA18PVY4aQM5VR3kMoiX1yHZlGeO+sd4SdHydGA/h7Wqmn0E45EOdWdZKFRPDVb6x5r7QGYsKaRzRGNowoNY8xLQMug4WuBh+3XDwPXxYw/ZowJGGMOAfuB1SIyAyg1xrxurMDknw6a4+zrSeByWwtZC2wwxrQYY1qBDQwVXkoMI3Vcy4VKtw0dfmpH8GdAfyFDDbvt51R3EGOgehSBOxYqCn209oSIRrMvz2Ai1LUM7bE+HoqnYPRUrTHmBID97BSrmQXUxWx3zB6bZb8ePD5gjjEmDLQDlSPsawgicquIbBGRLU1NTeM8pOzn1EhCI4vLFjicbPczfRQHpFMyXc1T/TR2WGUvapOoaZQX+YhEzZQTzsdarSi0kUykidAXcpuLPo0xEi/Q24wwPt45AweNuc8Ys8oYs6q6ujqhheYirT3DC41s794XjkRp7gokYJ6Kr2lEooZ/f34fq/7pOfac7EjZOjORxk4rgKAmiZpG5RRN8DvW2jth0xRAodeNSO76NOLRYJucsJ8b7fFjwJyY7WYD9fb47DjjA+aIiAcowzKHDbcvZRhauoJ9xeQGU5TnpieLS1k3dwWJmtEvfCV5HkT6mzU5/N9f7+Q7G96luSvA45vrhpmdmzQmoSrrYMqntNCYmBMcrFwXq3tf9v0mxys01gNONNPNwFMx4+vsiKj5WA7vTbYJq1NELrD9FTcNmuPs63rgBdvv8SxwlYiU2w7wq+wxJQ6BsNVNrDJHfRpO45vRLnwul1Ds8wwpWvjyvmauXjadtctq+e8dJwhPoeQ/J1S5aoIJabFUTMHWupGoob4tOZoGOI2Ysu83mUjI7aPA68AZInJMRG4BvgVcKSL7gCvtvzHG7AKeAHYDvwduM8Y4ovTzwP1YzvEDwDP2+ANApYjsB76MHYlljGkB7gA2249v2mNKHNp6rItkeY6ap5weBlUJ3C2XFngHxL+HIlFOtPdy+vQSrlsxi+auAK8fPJWytWYajZ0BKop8+DzJs0Y73f6yLYLqhy/u59ofvMpnHtrcZ7ZLlJMdfsJRw5yKiWsakL2VbketlGWMuXGYty4fZvs7gTvjjG8BlscZ9wM3DLOvB4EHR1ujAqe6rB9vrmoajtBIpHxDSb5ngE+jvq2XqIE55QVcdmYNJXkefrOtnksXTw3/l5XYlzwtA/o1jWxL8PvNtuO0dAd561gbj7xxlC9deXrCc4+1JCfc1iFbb+Q0IzxHaB2mWKFDUZ6H7mAka0Mkx9L4pjTfOyBP46j9Y59TUUi+1817zqhm46Gpo2kko8DeYAp8bvK9rqxL8GvtCXHFklouWVTFE1vqiIzh9+Ak9iXDpwGW0MjGGzkVGjnCcHWnHJwCaT1Z2ie8qTNASZ4noVDH0oKBLV+d2HrHrLB8ZhnHWntp75ka4aJWNnjyIqccKgp9tHRnzzk0xtDWE2RaoY8/Wz2XE+1+/ufdxtEn2hxs7sLtEmZOS865zFbzlAqNHKHFvhMfSdOA7AzxA0vTSPRuuWRQT4261h68bukL11020yoPvutEe/IXmmFEo4amJBcrdMi2rPDuYIRQxFBe6OWKpbVUFefxyBtHE56/5XAry2aWkueZWI6GQ3GeJzczwpXsoKUnhAhMGy7kNsuLFjZ1BhJygoNV6TbWEV7X0sPMaQV9ndYcobG7PvfzNRo7A4SjhhmjZNKPh/LC7Cpa6JjSyot8eN0uPrFmLs/vaWTvyc5R5wbCEbbVtXH+vIqkraeswEtLV/acPwcVGjlCS3eAaQXeYVtQZrum0TyGbmlW9FSor5VmXWvvgAJzlcV5TC/NZ9cUEBrb61qB4Tv2TYRsK4/uaEVOLtOnL5pHoc/NvS/uH3XuzmPtBMNRVs9PntCYUZZPdzCSdXXSVGjkCK3doWHDbcGKCYcs1jS6Em98U5JvtdLstpMZj7X0MKdiYMTL0pml7KrPffPU1iOt+DyuPu0qmWSf0LDD0gutqgHlRT4+ecFprN9Rz9FTPSPO3XTYivZPpqYxw65f5fSJyRZUaOQIrT3DZ4NDbCnm7HOE+0MROv3hhDWNMrs8+qmuAN2BMKe6g0MiXpbNLOVAUzf+LA0MSJQ3j7Zx1qyypNnhY6ko9NEZCBMIZ8c5bLM1jVgT7k0XnkbUwAt7Gkacu+lQC4triof1GY4Hx2R4QoWGkg5ae0J9d1DxcMxTPVnoeBtLjgbAOXOmAfDagVP9/Q8qhgqNSNSw52QnHf4Q1/3gVf7jhX05lSkeDEfZebyd8+ZOS8n+Hc22LUui0OI1KZtdXsjMsny2HGkddl40ath6pJXzk2iagv62xdmmaYya3KdkB+09wRFNENlcinksORoAZ9SWMKeigOd2N1Bgh+jOHSI0LBv/rvp22ntDbK9rY3tdGxt2N/Cdj61gUU1xEo8gPeyqt+zw581NTRua2KKFI3VTzBRa7WARRxN1OO+0ct4cQWg0dPrp9IdZMiO5Jj7nnKmmoaSF1p4Q0wpG1zSy0RHe7JQQSVDTEBGuWFLLK/ub+eGL+1lQVcTyQQJ1dnkBpfkedtV39EVR/ctHzuJoSw8fvOdlNmZ4mZH9jV0cau4ecZs3j7YB1kUxFWRb0cK2niCl+UODRVadVk59u5/6tt648xx/x2lJKh/i4PO4qCrO42RH/M/NVFRo5AD+UITeUGRER3ih13GEZ4f9OZaxahoAVy6pJRCO8m5DF391+SI87oFfdRGxneEdvHOig1nTCrhx9Vye/dJ7KMn38NM3jiT1GJKJMYZbHt7Mn/7HK+xvHD5cdFd9O7WleSnTAhwzj1NMMtNp7QnF9UmsPM0yO20dRttwKgqcVplcoQGWX0M1DWXSabcruk4bwadhlWLOzqqajk+jsjhxJ+T58ysozfewoKqIPzl7Ztxtls0sY8+JDt4+3t5neqgpyeea5TN44Z3GjPX/7Gvs4sipHroCYf78oc19///BHGvtHWKWSyZzKwrxuV1ZE7psZYMP/Y2cOaOEAq97WKFR19KDSyberS8e08vys86noUIjgznY1MXnfraVp7YfJxge3kHrxJ9PKxj5opqtRQubOgOUF3rxuhP/unrdLn70qZX84BPnDdEyHJbNLCUQjnKwuZulM0r6xj949gx6QxFe2JN4iYnJZMNuK9Ln7o+toK6ll+d2x4/8Od7ay6wUXOgc8r1uVsyZxsZD2VF8uqU7foSh1+3inDllvHl0eE1j5rSCMX3/EkU1DSWpPL3zBL/fdZIvPrad//3kjmG3axsUfz4c2VpV82S7n+llY7/4XbSwakTnpeMMBytvw+H8eRVUl+Txu7dOjPkzJ4Pn32ng7Nll/Ok5M5lW6OWNOP6XSNRwssPPrCRVZB2O8+eX8/bx9qy4GWnrCQ2rjZ9eW8LhYXxER1t6UqaxTS/Lp703lLFabTxUaGQwexu6mDWtgFsumc/6HfXD2q+d+POyUYRGYZY2fTnelpo75oXVReTZPSZihYvbJaxdVsuLe5syripwc1eAbXVtXH5mLS6XsGZ+BW/Eqdjb0OEnEjXMmpY68xTA6vmVRKKGbbbTPZNp7Qn2lXQfzOzyAjr84bimvqMtqTPzZWOuhgqNDGbvyQ7OnF7CbZctIt/j5od/PBB3u/5M11HMUz5PVib3WUIj+c5cj9vFmdNLKM7zDCgzAlYl3N5QhOPDRNSki9cOnMIYuOxMqxfIBQsqqWvp5VjrwIxmZ92p1jRWnlaOS2BThpea94ci9ASHDxZxkj+Ptw78f/cEwzR3BZLWeGkw00uzLytchUaGEgxHOdjUzRnTS6go8vFna+by1I76uF+utgSFRjaapzr8ITr94ZQ4IQE+ecFp3HLJfFyDwjAX2nka+5u6UvK542V3fQc+t4szp1ua0YULKwF44+BAv4Jz8UuFsI2lOM/D8lllGe/XcH4jw5mnnMZKg4VvKiOngL4y609tP55xNyjDoUIjQznY3EU4ajhjuuWg/cBZ04lEDbvjlPNu6wni87jI94787yzKwlLMJ9osIZkqoXHDqjlxu7ctrLaExoHGzBIau+rbWVxb3Ne69fSaEsoLvfxh18kB2ezOBShV5y2WNfMr2Ha0jd5g5mqxg4sVDsbRNI4N0jScHI1UmadmTSvgwgWVPLHlGB+65+Ws8G2o0MhQnHLNjtCYV1kEwKHmoYXVrLpTXkTiV7h1sKKnMveHHY/jbdbxptrMMpiKIh8VRT4OZJCmYYxhd30HS2P8Ly6X8OFzZ/OH3Q186N9f6QtPPtbaS0WRj0Jf6os+XLK4mmAk2lfULxMZTWiUF3op9LmHCo2W1AoNj9vFo7dewM9uWU1rT4ind55MyeckkwkJDRE5LCI7RWS7iGyxxypEZIOI7LOfy2O2v11E9ovIXhFZGzO+0t7PfhG5R+yrn4jkicjj9vhGEZk3kfVmC/5QhL0nO/G4hAVV1h1vRZGP0nxP3AiPtp7QqKYpsLr3ZZsj/LitaaQydHQ4FlYXcaBx5KzryaSpM8Cp7uCASC+Af/jQEr6/bgV7Tnby9E4r4itVwQPxWD2vAp/HxcvvNk3K59W19Iy5Rtgpu2/FcLk+IsKc8sIh5qndJzqoKPINKT2SbC5ZVMX8qiKe2FKX0s9JBsnQNC4zxqwwxqyy//4q8LwxZjHwvP03IrIUWAcsA64GfigiTunNe4FbgcX242p7/Bag1RizCLgbuCsJ681oHtl4hLO/8Qd+s+04C6qL+swQIsL8qqK4pSPaekIJfalL8730hiJZVdm1vq0Xr1sSLlaYTBbVFGeUT2PXCSuJbumgMGIR4U/PmUltaV5frsHx1p5JExoFPjer51Xw8r7mlH/W0VM9vP87L/L4GC+uTomQkZpRzS4vGKBpRKOG/9nbxHsWV42qxU8UEeGGVbPZdKhl1PIw6SYV5qlrgYft1w8D18WMP2aMCRhjDgH7gdUiMgMoNca8bqyuOT8dNMfZ15PA5ZLq/14aiUYN9798CATq2/2cXlsy4P15wwmN3pHLojs4ESB1LSP3Dsgk6tt6mV6WP8RRPRksrC6mpTuYMbWVnBpZS+IUphQRVp5WztYjrRhjLE1jEk16ly6uYm9DJw0dqY0CenJrHaGIYUfd2EJ869t6Kcn3UJI//M2VJTT6fxs7jrVxqjvIZWfWjHu9Y+H682YjAv+9o35SPm+8TFRoGOAPIrJVRG61x2qNMScA7GfnjM8CYm8Pjtljs+zXg8cHzDHGhIF2oHKCa85YXj3QzKHmbu766Fn8+43n8uVBDtr5VUXUt/cO0RRae0KUF42uacyrsvwih0dpOJNJpDqreSScCKpM8WvsPtHBnIoCSoe58J03t5xjrb3sOdmJPxSd1PN26WIrBDiV2kYkanhyq3Wp2JNAi9ZYjrf5Rz0fs8sLB+Rq/HFPIy6B955ePb4Fj5Ga0nwWVhezfYwCcbKZqNC42BhzHnANcJuIvGeEbePdKpoRxkeaM3DHIreKyBYR2dLUNDl21VTws9ePUFHk4wNnzeBPzpnJguqB5bnnVxVhTL9zDiznaFtPkLJRSogAzLed6cNlvqaL/Y1dvP87L/KN9buGaEH1bb2TEgEUj0X2+d+fIRFUO4+1DzFNxXKuXQL99l/tBGDVvNRUt43HmdNLqCrO4+V9yf39hSNRQrb/4rUDzdS3+5lbUcjek51j8mucaB/9ezQ47Pb5PY2sPK18QNOmVHP27DLeOtbe16o4E5mQ0DDG1NvPjcCvgdVAg21ywn52CvgcA+bETJ8N1Nvjs+OMD5gjIh6gDBgSomGMuc8Ys8oYs6q6enLuCpJNb9CqdfSRc2cN22VtfpUTQdV/0e8JRghFzKglRMDKGC8v9HLoVGYJjYdeO0RdSw+PbDzCn/zHK31mmHAkapXCSJPQmDWtgJJ8D28dS39b2JPtfo629IzYbnT5rFJ8bhfb69q4amktZ89OTfOleLhcwiWLKnllX3NSs+i/8osdXPbtFznc3M23n93LtEIvn3vvQgLh6Jg0ZuvmY+ScFSfstq6ll6bOALvqO3jfGZNjmnI4Z/Y0mrsCGZ0hPm6hISJFIlLivAauAt4G1gM325vdDDxlv14PrLMjouZjObw32SasThG5wPZX3DRojrOv64EXTCaL4Anwdn074ahhzYLhrW995qUYoTFaKGG8fWSSptETDPObbfX8yTkz2fCl91LodfOJ+9/gWGsP9W1+omZycg3i4XIJq+dVxK3tNNlstDOu18wf/vuR53GzfFYpLoG/u/qMyVpaH5curuZUd5DdJ5JT9fboqR6e2lHPsdZervreS7x1vJ27Pno258yxaobtOZnY5/QEw7T2hEb9Hi2oLsLjEnYeb2PrEeve1EmenCzOnm0d21vHMtdENRFNoxZ4RUR2AJuA3xljfg98C7hSRPYBV9p/Y4zZBTwB7AZ+D9xmjHGM858H7sdyjh8AnrHHHwAqRWQ/8GXsSKxcZLtdu2fFnOHvDkvzvVQV+wZoGo6TdqSy6LHMr8wsofHbHSfoCoS5cfVc5lUV8dO+ePUTbD1q/XDPmcQ75sFcsKCSQ83daS/zsPFQC8V5HpbMKBlxu7+54nT++cNnsahm5O1SwaWLqwB4ZX9y/BoPv34YtwjfueEcvC7h769Zwtpl01lUU4zbJew5kZhfo95JEB2l6GVRnoezZ5fx+oFTbDrUSp7HxfKYopaTwZIZpXhcwo4M0G6HY9yZP8aYg8A5ccZPAZcPM+dO4M4441uA5XHG/cAN411jNrG9ro1Z0wpGbTS0sLqYd2Lu5Bwn7YLqooQ+57TKIn617Tj+UIR8b3wz2GTx1rE27n7uXRbVFLPK7i63qKaExTXFvLL/FDNK8ykr8HLm9Mm/ADpcYGt+Gw+d4toVs0bZOnVsPHiKVfPKhy3z7vCeSXLaxqOmNJ8zp5fw8r4mPvfehRPaV1cgzBOb6/jAWTP46MrZ/OmKmX2lyfM8bhZWFyWsadSPITv+ggWV3PfSQTr8YVbMmdYX8j5Z5HvdnDG9JGc1DSWJbDvayoq5o99Rnz+vgrfrO/qS9Pae7MLrFk6rTExozKuy7LZH0xx2u6u+nevvfR2XCHd/bMWAOPiLF1Wx6dApXtnfzJr5FWkJt3VYOrOUkjxPWk1UzV0BDjR1j2iayhTee0Y1Gw+2DNs6NVEefu0wnYEwt1wyH2BIL4szp5fy9vGOPif5SPQLjdHrcF24sJJw1LC/sWtE/1EqOXv2NHbUtfeVMMk0VGhkAI0dfurb/Zw7gmnKYfX8CiJR05fEta+hkwVVxQk3iInnTE8Hf9jVQDga5Te3XcxZsweaAC5eVIU/FOV4W2/fnX66cLuE1fMrhhQEnExeO2D7Mxak5yI2Fj51wWkA/Pjlg+PeR4c/xH0vHeTyM2s4Z5jfxAfOms7JDj/fWL9r1Eij+nY/LiGhtrcrTyvH67ZuUiYz+iyWj58/BxH44D0v92X4ZxIqNDKAN21/xrkJaBrnnVaO2yVstC9i7zZ2sri2eJRZ/cRzpo+FZ3ed5L+S0D978+EWlswojWuOW7OgAretXaRbaIB193mouTttVUh/saWOmWX5afXtJMrs8kKuO3cWj246Oq7e4cYYfvDCftp7Q3ELSTpcvXwGn3vvQh7ZeJT1oyTD1bf1Uluan9CNVaHPwzmzpyFi/dbSwYo503j6ry9lQU0xf/nImwkJxslEhcYk09jh56J/eZ5rf/Aqd/x2N1/7zU6+8osdlOR7BnSSGw6nFPWmQy10B8LUtfRyRm3iNv/SfC9zKgrYNI5S1uFIlH986m3u+v2eCYVVBsNR3jzaOqz6X5rv5ZzZZWn3Zzg4foKXJqm2UiyHm7t5eV8z61bP7ROkmc7n32eFxP7z794Z08UuEjV89Zc7+c+XDvKR82axfNbIv4e/W3sGs6YV8MwoRf7q23pHLB8ymFsumc+tly4YNolyMphTUcgvPnshn1gzl4deOzxsK9p0oEJjkvn2H/bS1BUAY3hk4xGe3HqMSxdX8avPX5SwY3rN/Aq217Wxy85nWDwGoQFw1dLpvLyvecy9NZ7f00hDR4BOf3hCuR5v17fjD0VZPX94c8s/fGgp3/3YOWn1Zzgsrilmeml+WoTGo5uO4nYJHz9/zugbZwgLq4v54uWL+dW24zzwyqGE5724t5HHt9Tx2fcu4N+uHxJjMwSXS7hoYSWvHzw17E2MPxThnRMdfVWiE+Gas2Zw+weWJLx9qvB5XNz+gSUUeN08ufV4upfThwqNSWRXfTu/2HqMT180j6e+cAl77riGPXdcw72fXDmmC/8FCyoIRqJ877l3ATh9DOYpgLXLphOMRHlxb+PoG8fw841HKbAF21hr/8Sy2dZyRnI0nju3nMuX1I77M5KJiHDp4ipe3d885uqqE6HDH+LxLXVcsaQmIXt8JvHX71/M2mW1/NPv3uG2n785pHpsPDYfbsXrFr50xekJa1UXLaqkvTc0bG7Ib986QWtPiI+unB33/UynOM/D1cun89u36jOm0KgKjUnk+8/to6zAyxcuWzyh/bz39BrWzK/gtQOn8HlcCUdOOaw8rZzKIh/P7mpIeE5dSw8v7WvilkvmU+hzTyhLevPhFuZXFY0aXpxJvOf0ajr84UmNn7//pYO09YQm/H1JBy6X8P115/KlK07nud0NXPbtF/nHp96mwz+0B7fDm0dbWTazbEyh4BcusHJD4kW3GWN4+LXDLK4p5qJJTtJLJh89bzad/jAbdif+e00lKjQmibqWHp57p4FPrJlLWYKJeMPhtn+Q5YVezqgtGbOt2+0Srlxayx/3NBIIJ3b38uimowjwZ2vmctassnEXVesJhnntwKlJz7SdKJcsqsIl8L9/sYOHXj2Ucsdkc1eA+185xAfOmj4kuixbyPe6+eIVi/njV97HDavm8F9vHOGa770cV0sNRaK8dayN8+aOzfk8vSyfBVVFfRFmDrvrO/j6+l3sPN7Opy+el/LS5qnkwoWVzJpWwGObj6Z7KYAKjUnjv944gojwSTskcaJML8vnF5+7kO98bHTbbzwuO7OGrkCYHXWj3zmHIlGe2HKM959Zw8xpBZwzZxq76zsIhsduqvnDrgZ6ghGuPWfmeJadNsqLfNxz47kU53v4xn/vTnkl0v964wi9oQhfvnLyy4Ekm5nTCvjnD5/Fk5+/CGMMf/P49iHfnXdOdOAPRVk5joilCxdW8sbBU30msA5/iI/95+s8tqmOq5bW8pFzs9M05eB2CTeunsOr+09lRMVlFRpJJt4dqD8U4fEtdaxdVsuMUUoZjIVFNSVDem4kygXzKxGxKoeOxobdDTR3BfizNXMBq6xHMBIdkJmeKL/adpxZ0wrSljg1ET509kwe+V9ryPe6+kp0pwJjDE9tr2fN/AoW1YzNX5XJnDe3nDs/chaHmrv56euHB7y39YgVHXTeaWMPK/70RfPwuIRP3r+Rps4AT2yuoysQ5pefv4j7blpFgS+9lQ+SwcfOn4PXLTzyxlBt41BzN3/x0y3saxhbufjxokIjiby6v5mLvvXCkPLQ//NuE209IW5cPTdNKxtKWaGXZTNLef3AyJnOwXCUH764n1nTCnjv6VbFz/PnleNzu/jab96mrSfxBkWNHX5e2dfEdefOzIioqPFQku/lmuUzWL8jdY7JncfbOdTczXVpLFuSKi47o4b3nVHN95/fR3uP5d+IRg0v72tmZln+uG6qFteW8JM/P5+GjgA3/vgNfvLqYc6fV561Zr141JTkc/XyGfxia92AqMdwJMqXn9jOht0N3PLwFlonoWGYCo0k0dId5EuPb+dEu5+/enTbgL4Qv3/7JNMKvRmRqBbLRQur2Ha0bcSL37ee2cPbxzv42geX9PlOakrz+dGnzmNvQyefeWhzwvb99TvqiRr4cJabC65faTkm/5Aix+RvttXjc7u4ZvmMlOw/3XzlqjPo9IdZ/1Y9Hf4Qn35oMy/saeRPVozfZLnytAp+8ufnc6Ktl+NtvXzm4vlJXHFm8BeXzqcrEOauZ/b0jd338kG2HW3js+9dwMl2P7f9/M2klqaPhwqNCRKKRHl881E+cf9G2npC3PuJ84hGDV/4+ZuEI1GC4SjPvdPAlUtqEy71MVlcuKCSYCTaZxqIxRjDf7ywjwdfPcSnL5rHNWcNvIC9/8xavvmny3jzaBsvJdit7dfbjnP27LKsN7lcuKCSeZWF3Pm73UNqLP1xbyPbJpCI1RuMsH5HPe87o3rCAROZyrKZpZw5vYQnt9Tx3T+8yyv7mrjjuuV89eozJ7TfCxZU8vO/uIC/ev8irlyaGeHayeTs2dP4zMXz+dkbR3h1fzPvnOjg7g3v8sGzZvDVq8/kWx89i2tXpF6LH3eVW8VSq//uybf49bbjLKgu4u6Pr+Cas2YQjhr+6tFtPPTaYRbXltDpD3P18unpXu4Qzp9vlev47M+2Ul7k5Zefu4ia0nwiUcPfPrGd32yv59oVM7n9A/F/zB85bzZ3P/cu9710YNSWmO82dLKrvoN//NDSVBzKpOJyCT/61EpuuPd1PvnARj5z8XymFXr5w64G1u+opzTfw3Nffi8148it+OGL+2nuCvAX71mQgpVnBiLCDavmcMdvd7PzeDs3rp7bV7NqopwzZ9qw9apyga9cdQYv7Gnkz3+ymeqSPMoKfNxx3XJEhI+cNzkafGbd+mYZ//rsXn697Th/e+XpPP/l9/LBs6278Q+dPYMrltTwb8/u5e+e3EGRz83Fi6rSvNqhFOd5uP2aM7l6+XSOt/bywCtWKOk31u/iN9vr+dsrT+d7H18xbCdBn8fFZy6ez6v7T/HqKD0UfvXmcdwu4U8nYILIJM6cXsp/3nr4QkQAAAniSURBVLSSYDjK137zNl/4+Tae3nmCWy6ZTyAc5e9/vTNhs11jh5+/fWIHX/3lW/znSwe5bsXMrAwUGAvXrZiJxyUU+Tx8eYQaU8pACnxufvG5C3nP6dXUt/fyrY+cRUXR5LWjBdU0xs2G3Q386H8O8Gdr5vKF9y8aEAcuItxx3XK+8osdFOd5WLtsetp7VwzH/7rUuqMNhqP81xtH6PCHeXTTUT77ngX81eWjJ5XduGYuP3vjCJ96YCM3XTiPmy+a11dJ18EfivCrN4/xnsVVVBVnT0LfaFy0sIqX/+4yDp/qIRCO2O1hvUwvzefOp9/h/pcPjaoxbD7cwl8+8iad/hBel4t8u3RErlNZnMf//eASppfmU5lD34nJoKo4jx/ftJJT3cG0/J4kk6onJoNVq1aZLVu2jGvuvoZOFtUUj5oIdLCpi4/c+xqzphXwq7+8aNg78Wxi78lO1n7vJQA+c/F8vvbBJQnbRjv8If7l6T08saWOSNTwT9ctH5CP8sMX9/Ovv9/LY7dekHHBAKkgGjV84dE3eXrnSe64bjk3nj9nSPOkaNRw38sH+bdn9zK3opAffXIli2qK8YciFOXpvZwy+YjIVmPMqlG3U6FhcbCpi6u//zJXLqnlXz56VtwKl9Go4b82HuFfnt6Dz+Pi1395EQuqs9upG8sP/rif0nwPn7zgtHFl0DZ2+PniY9t5u76dF7/yPiqL8zjVFeB9//YiaxZUcP/N56dg1ZmJPxTh0z/ZxBsHW5hRls+nLjyNq5bW4hKh0x/mnuf38fyeRj541gy+9dGzKEljRVVFARUaY54XjRr+86WDfPsPe5k5LZ/vfmwFz+1uYMM7DZxRW8LcykJ2HmvntQOneM/p1dz10bOSmqiXK+xv7GTt917mmuXTuXJpLd/5w7vUt/XyzBcvHXM13mwnEjX8cU8jP3ntEK/uH5gP43ULX/vgUm66cHwCWlGSjQqNcbL1SCt//ei2voY7Fy2s5HhbL/VtveR73Pz9B5ew7vw5+kMfgX9++h3ue8nq3Da3opB/u/5s1kwBs9RI7LOjxwyG4jwvp9cWj7nQpKKkkpwSGiJyNfB9wA3cb4z51nDbTlRoALT3hLjnhX1csqiKy860sqCNMUQNWdMIJ50YY3i3oYveUIQzp5dkbBCAoij9JCo0Mt7jJiJu4AfAlcAxYLOIrDfG7E7VZ5YVevmHQfkEIoJb5UVCiAhnZEDHPUVRkk825GmsBvYbYw4aY4LAY8C1aV6ToijKlCQbhMYsoC7m72P2mKIoijLJZIPQiGcUGuCIEZFbRWSLiGxpapr8Ps6KoihThWwQGseAOTF/zwbqYzcwxtxnjFlljFlVXT1yDSRFURRl/GSD0NgMLBaR+SLiA9YB69O8JkVRlClJxkdPGWPCIvIF4FmskNsHjTG70rwsRVGUKUnGCw0AY8zTwNPpXoeiKMpUJxvMU4qiKEqGkBUZ4WNBRDqBvSNsUga0J/Ejk72/KiCxVniJkcnHq+cus/YHeg4nSraevyqgyBgzeiSRMSanHsCWUd6/L8mfl+z9jbj+DFhf0van5y6z9qfncOqev7Gseyqap/47w/eXbDL5ePXcZdb+UkGmH3Omn8OMO95cNE9tMQkU3cpUsn396UTP3cTRczgxsvX8jWXduahp3JfuBUyQbF9/OtFzN3H0HE6MbD1/Ca875zQNRVEUJXXkoqahKIqipAgVGilGROaIyB9F5B0R2SUiX7THK0Rkg4jss5/L7fErRWSriOy0n99vjxeKyO9EZI+9n2EbUeUKyTp39nu/F5Ed9n5+ZPdpyXmSeQ5j9rleRN6e7GNJB0n+Dr4oIntFZLv9qEnXcU2IZIZz6SNuKNsM4Dz7dQnwLrAU+Ffgq/b4V4G77NfnAjPt18uB4/brQuAy+7UPeBm4Jt3Hlw3nzv671H4W4JfAunQfX7adQ3vsI8DPgbfTfWzZdv6AF4FV6T6mCZ+TdC9gqj2Ap7C6EO4FZthjM4C9cbYV4BSQF+e97wN/ke7jybZzB3ixwg4/nu7jybZzCBQDr9gXzSkhNJJ8/nJCaKh5ahIRkXlYdyIbgVpjzAkA+zmeqvpRYJsxJjBoP9OAPwGeT+V6M4lknDsReRZoBDqBJ1O85IwjCefwDuA7QE/KF5uBJOn3+xPbNPUPIpKVDaRVaEwSIlKMZRb5G2NMRwLbLwPuAj47aNwDPArcY4w5mIq1ZhrJOnfGmLVYd4V5wBBbfS4z0XMoIiuARcaYX6d0oRlKkr6DnzDGnAVcaj8+lYq1phoVGpOAiHixvnCPGGN+ZQ83iMgM+/0ZWHfAzvazgV8DNxljDgza3X3APmPM91K/8vST5HOHMcaP1Y9lyvSZT9I5vBBYKSKHsUxUp4vIi5NzBOklWd9BY8xx+7kTyy+0enKOILmo0Egxtgr6APCOMea7MW+tB262X9+MZSt1TE+/A243xrw6aF//hFVw7G9Sve5MIFnnTkSKY37gHuADwJ7UH0H6SdY5NMbca4yZaYyZB1wCvGuMeV/qjyC9JPE76BGRKvu1F/gQkJURaJrcl2JE5BKsSKedQNQe/nssu+gTwFzgKHCDMaZFRL4G3A7si9nNVVgRU3VYFzvHRvofxpj7U34QaSKJ506A32KZpdzAC8CXjDHhyTiOdJKsc2iMib2Tngf81hizPOUHkGaS+B3sBl7CCsRwA88BXzbGRCbjOJKJCg1FURQlYdQ8pSiKoiSMCg1FURQlYVRoKIqiKAmjQkNRFEVJGBUaiqIoSsKo0FCUSUZEPiciN41h+3lTpaqskvl40r0ARZlKiIjHGPOjdK9DUcaLCg1FGSN2ctvvsRK8zsUql30TsAT4LlY12Gbg08aYE3a5jdeAi4H1IlICdBljvm3XdPoRVun7A8BnjDGtIrISeBCrOOArk3d0ijIyap5SlPFxBnCfMeZs+P/t3atKRFEUh/FvoRaLFqtP4AsYxGy3iMgwWbAaTUZBtGhRvBTxCcQ0zS6DzWIwekPrMsweGAYvZw4jDPj94uKwOLucP2uHdXgF1oF9YDkzux/87Z7npzNzMTN3+vqcApulzy2wVerHwEZmzv/lIaRBOWlI9Tz07BY6p7NaYg64Lhuvx4DHnucv+htExBSdMGmV0glw+UX9DFga/hGkwRkaUj39+3fegPYPk8H7AL3ji/7SSPB6SqpnNiK6AbEC3AAz3VpETJR/KnwrM1+Ap4hYKKU1oJWZz8BLWZYHsDr815fqcdKQ6rkDGhFxSGej6T5wBeyV66VxYBdo/9KnARxExCRwDzRLvQkcRcRH6SuNBLfcSgP6T6vBpX5eT0mSKnPSkCRV5qQhSarM0JAkVWZoSJIqMzQkSZUZGpKkygwNSVJln0RdJYYhM1/tAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -252,10 +2309,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 15, + "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", @@ -274,7 +2329,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +2353,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 17, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 17, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAD8CAYAAABU4IIeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAG35JREFUeJzt3X+Q1HV+5/HnCweHjYsKCC4/1OFqORN07zROoVveD3HDj81tqXvn7rEa5SpWYVy3ypypU6hw5Z6Si1xdrTnirdFa3bjrD/Q2S0likCBCnZcQYCg1ikoGbwkSWWfIoGCq4HbkfX/0Z7Snne7p6f729LdnXo+qru7+9Ofz6c98GPo9n1/fVkRgZmaWlQnNboCZmY0tDixmZpYpBxYzM8uUA4uZmWXKgcXMzDLlwGJmZplyYDEzs0w5sJiZWaYcWMzMLFNtzW7AaDrnnHOio6Oj2c0wM2spe/bsORIR06vNP64CS0dHB11dXc1uhplZS5H0dyPJ76kwMzPLlAOLmZllyoHFzMwy5cBiZmaZcmAxM7NMObDkQM+xE3zz4R30HD/R7KaYmdXNgSUH1m3tZveBPta92N3sppiZ1W1cnWPJmwtXb+Jk/6lPnj+x8yBP7DxIe9sE9q35ahNbZmZWO49YmujluxZyzSWzmDSx8M8waeIErr1kFi/fvbDJLTMzq50DSxPNOHMSk9vbONl/iva2CZzsP8Xk9jZmTJ7U7KaZmdXMU2FNduSjk9x4+QXcsOB8ntp1kF4v4JtZi1NENLsNo6azszN8rTAzs5GRtCciOqvN76kwMzPLlAOLmZllyoHFzMwy5cDSAnwy38xaiQNLC/DJfDNrJd5unGM+mW9mrcgjlhzzyXwza0UOLDnmk/lm1oocWKrQzMXzgZP5G759JTdefgG9H50c9TaYmY2ET95XYfWG13ly10FuXHA+a77+pQa0zMwsv0Z68t6L9xV48dzMbOQ8FVaBF8/NzEauqsAi6YCk1yW9KqkrpU2VtEVSd7qfUpR/laT9kvZJWlKUflmqZ7+kdZKU0tslPZPSd0rqKCqzPL1Ht6TlRelzU97uVPb0+rtjMC+em5mN3EhGLAsj4pKiebaVwNaImAdsTc+RNB9YBlwELAW+L+m0VOYhYAUwL92WpvRbgKMR8UXgAWBtqmsqcA9wObAAuKcogK0FHkjvfzTVkTkvnpuZjUxVi/eSDgCdEXGkKG0fcFVEHJY0E9geERdKWgUQEb+f8m0GvgscALZFxC+n9G+l8rcO5ImIHZLagJ8D0ykEqKsi4tZU5mFgO7Ae6AW+EBH9kr6cyn8yOhqKL5tvZjZyjbpsfgB/IWmPpBUp7dyIOAyQ7mek9NnAu0VlD6W02elxafqgMhHRD3wITKtQ1zTgg5S3tC4zM2uianeFXRkR70maAWyR9HaFvBoiLSqk11KmUl2DG1MIhCsAzj///KGymJlZhqoasUTEe+m+B9hAYb3j/TQFRrrvSdkPAecVFZ8DvJfS5wyRPqhMmgo7C+irUNcR4OyUt7Su0rY/EhGdEdE5ffr0an5cMzOrw7CBRdIZkiYPPAYWA28AG4GBXVrLgefS443AsrTTay6FRfpdabrsuKQr0m6wm0vKDNR1PfBSFBZ/NgOLJU1Ji/aLgc3ptW0pb+n7m5lZE1UzFXYusCHtDG4DnoqIFyTtBp6VdAtwEPgGQETslfQs8CbQD9weER+num4D/hj4HLAp3QAeBX4saT+FkcqyVFefpPuA3SnfvRHRlx7fDayXtAZ4JdVhZmZN5ku6mJlZRY3aFWZmZlYVBxYzM8uUA4uZmWXKgcXMzDLlwGJmZplyYDEzs0w5sJiZWaYcWMzMLFMOLGZmlikHFjMzy5QDi5mZZcqBxczMMuXAYmZmmXJgMTOzTDmwmJlZphxYzMwsUw4sZmaWKQcWMzPLlAOLmZllyoHFzMwy5cBiZmaZcmAxM7NMObCYmVmmHFjMzCxTDixmZpYpBxYzM8uUA4uZNVXPsRN88+Ed9Bw/0eymWEYcWMysqdZt7Wb3gT7Wvdjd7KZYRtqa3QAzG58uXL2Jk/2nPnn+xM6DPLHzIO1tE9i35qtNbJnVyyMWM6uoUVNVL9+1kGsumcWkiYWPoUkTJ3DtJbN4+e6Fmb6PjT4HFjOrqFFTVTPOnMTk9jZO9p+ivW0CJ/tPMbm9jRmTJ2X6Pjb6PBVmZkMajamqIx+d5MbLL+CGBefz1K6D9HoBf0xQRDS7DaOms7Mzurq6mt0Ms5bQc+wEa/78Lf5i78858YtTTJo4gSUXfYHf/Te/4lHFOCNpT0R0Vpu/6qkwSadJekXSn6XnUyVtkdSd7qcU5V0lab+kfZKWFKVfJun19No6SUrp7ZKeSek7JXUUlVme3qNb0vKi9Lkpb3cqe3q1P4uZDc9TVVarkayx3AG8VfR8JbA1IuYBW9NzJM0HlgEXAUuB70s6LZV5CFgBzEu3pSn9FuBoRHwReABYm+qaCtwDXA4sAO4pCmBrgQfS+x9NdZhZhgamqjZ8+0puvPwCej862ewmWQuoaipM0hzgceD3gDsj4muS9gFXRcRhSTOB7RFxoaRVABHx+6nsZuC7wAFgW0T8ckr/Vip/60CeiNghqQ34OTCdQoC6KiJuTWUeBrYD64Fe4AsR0S/py6n8J6OjoXgqzMxs5Bo1FfYHwF3AqaK0cyPiMEC6n5HSZwPvFuU7lNJmp8el6YPKREQ/8CEwrUJd04APUt7SuszMrImGDSySvgb0RMSeKuvUEGlRIb2WMpXqGtwYaYWkLkldvb29Q2WxFudLgpjlSzUjliuBayQdoDAFdbWkJ4D30xQY6b4n5T8EnFdUfg7wXkqfM0T6oDJpKuwsoK9CXUeAs1Pe0roGiYhHIqIzIjqnT59exY9rrcaXBDHLl2EDS0Ssiog5EdFBYc3jpYj4DWAjMLBLaznwXHq8EViWdnrNpbBIvytNlx2XdEXaDXZzSZmBuq5P7xHAZmCxpClp0X4xsDm9ti3lLX1/GycuXL2JjpXP88TOg0QUzll0rHyeC1dvanbTzMa1ek7e3w8sktQNLErPiYi9wLPAm8ALwO0R8XEqcxvwA2A/8A4w8AnwKDBN0n7gTtIOs4joA+4DdqfbvSkN4G7gzlRmWqrDxhFfEsQsn0Z08j4itlPYlUVE/APwlTL5fo/CDrLS9C7g4iHSTwDfKFPXY8BjQ6T/XwpbkG2cqvacRc+xE3zn6Vd48IZLc3kGI+/tMxspXyvMWlo15yzyvgaT9/aZjZQv6WJjVum1rgbk5bLseW+f2YCGXdLFrNXkfQ0m7+0zq5UDi41Zeb/WVd7bZ1YrXzbfxrS8X5Y97+0zq4XXWMzMxrAsdh16jcXMxhVf0qeyZuw69FSYmbW04g/ONV//UrObkxuj8Q2g5XgqzJrOBwStFt6uXVmW3wDqqTBrOT4gaLXwdu3Kmrnr0FNh1jTNHKqPJo/IGsPbtYfXrF2HHrFYw5VbXB0vf3F6RNY4/urkyh6+qZM1113M/Flnsua6i3n4pqpns+riEYsNq96/uMstro71vzjHy4ismYo/KNdc95nr21qTOLDYsGrddVPNB+tYPiD48l0Lyy6emo1lDixWVr1/cVfzwTqW/+Ic6yMys3K8xmJl1bsG4g/W1lgD8AFDy5pHLFZWFoFhLE91VaMVRmQ+YGhZ8wFJq+jWH3cxffKkQYFhtHaWWGP5gKFVa6QHJB1YzMapLE9m29jmk/dmVhWvgVmjeI3FbBwb72tg1hieCjMzs4o8FWZmZk3lwGJmZplyYDEzs0w5sJiZWaYcWMzMLFMOLGZmlikHFrNh+CKNZiPjwGI2DH8DpNnI+OS9WRn+Bkiz2njEYlZGvd9HYzZeObBY7jVrjcMXaTSrzbCBRdIkSbskvSZpr6T/ktKnStoiqTvdTykqs0rSfkn7JC0pSr9M0uvptXWSlNLbJT2T0ndK6igqszy9R7ek5UXpc1Pe7lT29Gy6xPKmmWscrfANkGZ5M+xFKNOH/xkR8ZGkicD/Ae4A/i3QFxH3S1oJTImIuyXNB54GFgCzgBeBfxoRH0valcr+NfDnwLqI2CTp28A/i4jfkrQM+HpE/HtJU4EuoBMIYA9wWUQclfQs8NOIWC/pj4DXIuKhSj+LL0LZWvxFVGb5kPlFKKPgo/R0YroFcC3weEp/HLguPb4WWB8RJyPiZ8B+YIGkmcCZEbEjCtHsRyVlBur6CfCVFNCWAFsioi8ijgJbgKXptatT3tL3tzHCaxzZ8HZpG21VrbFIOk3Sq0APhQ/6ncC5EXEYIN3PSNlnA+8WFT+U0manx6Xpg8pERD/wITCtQl3TgA9S3tK6Stu+QlKXpK7e3t5qflzLCa9xZMPbpW20VbXdOCI+Bi6RdDawQdLFFbJrqCoqpNdSplJdgxMjHgEegcJU2FB5LL/8RVS183Zpa5YRnWOJiA8kbQeWAu9LmhkRh9M0V0/Kdgg4r6jYHOC9lD5niPTiMocktQFnAX0p/aqSMtuBI8DZktrSqKW4LhtDHr7p02ndNddV+nvGSr1818Ky32lv1kjV7AqbnkYqSPoc8GvA28BGYGCX1nLgufR4I7As7fSaC8wDdqXpsuOSrkhrJDeXlBmo63rgpbQOsxlYLGlK2nW2GNicXtuW8pa+v5kxelOJeV/DyXv76pXHn6+aNZaZwDZJfwPsprDG8mfA/cAiSd3AovSciNgLPAu8CbwA3J6m0gBuA35AYUH/HWBTSn8UmCZpP3AnsDLV1Qfcl953N3BvSgO4G7gzlZmW6miKPP7DmsHobJfO+xpO3ttXrzz+fP7O+wys3vA6T+46yI0LzmfN17+Uef1meZT37eB5b1+9RvPnG+l2YweWOoz1X1yzSnqOnSi7hpOHnXt5b1+9RvPny/wci5VX7TkLT5VZI/mSN0PLe/vqleefz4GlDtX+w+ZxDtTGDl/ypry8t69eef35PBVWp1t/3MX0yZMGnbMY2CLrqTJrJP9+2WjxGksFo32tsFaZ4+05doLvPP0KD95waa7aZZW1yu9Xq/P/D6+x5Eqe50CLeaquNbXK71er8/+PkfM3SDZYni9J4kt+tL48/361Ov//qJ2nwsYxT6WYlef/H5/yVJhVzVMpZuX5/0ftHFjGgUrnHPK6XdEsD/z/ozaeChsHfMkZM6vHSKfCvHg/hnnxsTreTmqWLU+FjWH+at/qeDupWbY8YhnDvPhYmUd0Zo3hEcsY58XH8jyiM2sMj1jGOH+1b3ke0Zk1hgOLjWs+uW6WPW83NjOzinzy3szMmsqBxczMMuXAYpZz/mrrxnL/Zs+BxSznfICzsdy/2fPivVlO+auHG8v9Wz0v3puNET7A2VhZ9a+n0j7LgcWsTo36YPEBzsbKqn89lfZZPiBpVqfiD5asv5bABzgbq57+9bXmyvMai1mNPEdfnbH6tQTj6auLvcZiNkq8BlKdsTpV5KnK8jwVZlYjf7BUNh6mijxVOTQHFrM6+IOlvJfvWlh2qmis8NXDh+bAYlYHf7CU5xHd+OU1FjNrGH/RXP1a8ZyMd4WZmeXY6g2v8+Sug9y44PzMt7NXK/NdYZLOk7RN0luS9kq6I6VPlbRFUne6n1JUZpWk/ZL2SVpSlH6ZpNfTa+skKaW3S3ompe+U1FFUZnl6j25Jy4vS56a83ans6dX+0GZmo6XWEceFqzfRsfJ5nth5kIjC5oeOlc9z4epNDWppdqqZCusHficifgW4Arhd0nxgJbA1IuYBW9Nz0mvLgIuApcD3JZ2W6noIWAHMS7elKf0W4GhEfBF4AFib6poK3ANcDiwA7ikKYGuBB9L7H011mJnlSq3brVt5O/uwi/cRcRg4nB4fl/QWMBu4FrgqZXsc2A7cndLXR8RJ4GeS9gMLJB0AzoyIHQCSfgRcB2xKZb6b6voJ8GAazSwBtkREXyqzBVgqaT1wNXBD0ft/l0LgMjNrunq3W7fy5ocRLd6nKapLgZ3AuSnoDASfGSnbbODdomKHUtrs9Lg0fVCZiOgHPgSmVahrGvBByltaV2mbV0jqktTV29s7kh/XzKxmWYw4WnXzQ9XbjSV9HvgT4Lcj4lhaHhky6xBpUSG9ljKV6hqcGPEI8AgUFu+HymNmlrUsRhytup29qhGLpIkUgsqTEfHTlPy+pJnp9ZlAT0o/BJxXVHwO8F5KnzNE+qAyktqAs4C+CnUdAc5OeUvrMjPLhVYdcdRr2BFLWut4FHgrIr5X9NJGYDlwf7p/rij9KUnfA2ZRWKTfFREfSzou6QoKU2k3A39YUtcO4HrgpYgISZuB/1q0YL8YWJVe25byri95fzOzXGjVEUe9qhmxXAncBFwt6dV0+3UKAWWRpG5gUXpOROwFngXeBF4Abo+Ij1NdtwE/APYD71BYuIdC4JqWFvrvJO0wS4v29wG70+3egYV8ChsF7kxlpqU6rAla8QCXmTWOD0ha3fJwgMvMGmekByR9rTCr2Xi4eq2ZjZyvFWY1a+UDXGbWOA4sVrNWPsBlZo3jqTCri7+PxMxKefHezMwq8nfem5lZUzmwjAE+R2JmeeLAMgbUelluM7NG8OJ9C/M5EhsPeo6d4DtPv8KDN1zqHYctwiOWFuZzJDYeeETeejxiaWE+R2JjmUfkrcsjlhY3Xi/LbWOfR+StyyOWFjdeL8ttY59H5K3LgcXMcstXdmhNPnlvZmYV+eS9mZk1lQOLmZllyoHFzMwy5cBiZmaZcmAxM7NMObCYmVmmHFjMzCxTDixmZpYpBxYzM8uUA4uZmWXKgcXMzDLlwGJmZplyYDEzs0w5sJiZWaYcWMzMLFMOLGZmlikHFjMzy5QDi5mZZWrYwCLpMUk9kt4oSpsqaYuk7nQ/pei1VZL2S9onaUlR+mWSXk+vrZOklN4u6ZmUvlNSR1GZ5ek9uiUtL0qfm/J2p7Kn198VZmaWhWpGLH8MLC1JWwlsjYh5wNb0HEnzgWXARanM9yWdlso8BKwA5qXbQJ23AEcj4ovAA8DaVNdU4B7gcmABcE9RAFsLPJDe/2iqw8zMcmDYwBIR/xvoK0m+Fng8PX4cuK4ofX1EnIyInwH7gQWSZgJnRsSOiAjgRyVlBur6CfCVNJpZAmyJiL6IOApsAZam165OeUvf38zMmqzWNZZzI+IwQLqfkdJnA+8W5TuU0manx6Xpg8pERD/wITCtQl3TgA9S3tK6PkPSCkldkrp6e3tH+GOamdlIZb14ryHSokJ6LWUq1fXZFyIeiYjOiOicPn16uWxmZpaRWgPL+2l6i3Tfk9IPAecV5ZsDvJfS5wyRPqiMpDbgLApTb+XqOgKcnfKW1mVmZk1Wa2DZCAzs0loOPFeUvizt9JpLYZF+V5ouOy7pirRGcnNJmYG6rgdeSuswm4HFkqakRfvFwOb02raUt/T9zcysydqGyyDpaeAq4BxJhyjs1LofeFbSLcBB4BsAEbFX0rPAm0A/cHtEfJyquo3CDrPPAZvSDeBR4MeS9lMYqSxLdfVJug/YnfLdGxEDmwjuBtZLWgO8kuowM7McUGEAMD50dnZGV1dXs5thZtZSJO2JiM5q8/vkvZmZZcqBxczMMuXAYmbWRD3HTvDNh3fQc/xEs5uSGQcWM7MmWre1m90H+lj3Ynezm5KZYXeFmZlZ9i5cvYmT/ac+ef7EzoM8sfMg7W0T2Lfmq01sWf08YjEza4KX71rINZfMYtLEwsfwpIkTuPaSWbx898Imt6x+DixmZk0w48xJTG5v42T/KdrbJnCy/xST29uYMXlSs5tWN0+FmZk1yZGPTnLj5Rdww4LzeWrXQXrHyAK+D0iamVlFPiBpZmZN5cBiZmaZcmAxM7NMObCYmVmmHFjMzCxTDixmZpapcbXdWFIv8HdlXj6Hwtce55XbVx+3rz5uX31avX0XRMT0aisbV4GlEkldI9mnPdrcvvq4ffVx++oz3trnqTAzM8uUA4uZmWXKgeVTjzS7AcNw++rj9tXH7avPuGqf11jMzCxTHrGYmVmmxmxgkfSYpB5JbxSl/XNJOyS9LulPJZ2Z0idKejylvyVpVVGZ7ZL2SXo13WY0oX2nS/phSn9N0lVFZS5L6fslrZOkLNqXcRsz70NJ50nalv699kq6I6VPlbRFUne6n1JUZlXqp32SlhSlZ96HGbev6f0naVrK/5GkB0vqanr/DdO+PPTfIkl7Uj/tkXR1UV156L9K7Rt5/0XEmLwB/wr4VeCNorTdwL9Oj38TuC89vgFYnx7/EnAA6EjPtwOdTW7f7cAP0+MZwB5gQnq+C/gyIGAT8NUctjHzPgRmAr+aHk8G/haYD/w3YGVKXwmsTY/nA68B7cBc4B3gtEb1Ycbty0P/nQH8C+C3gAdL6spD/1VqXx7671JgVnp8MfD3Oeu/Su0bcf9l1tF5vAEdDP5QPMan60rnAW+mx98C/pTCF59NS/8IUxv1S1lD+/4n8BtF+bYCC9Ivz9tF6d8CHs5TGxvdh0Xv9xywCNgHzExpM4F96fEqYFVR/s3pP3PD+7Ce9uWl/4ry/QeKPrjz0n/l2pe3/kvpAv6Bwh8Rueq/0vbV2n9jdiqsjDeAa9Ljb1D4YAT4CfCPwGHgIPDfI6KvqNwP0xDwP2cxTK2hfa8B10pqkzQXuCy9Nhs4VFT+UEprpJG2cUDD+lBSB4W/uHYC50bEYYB0PzBsnw28W1RsoK8a3od1tm9As/uvnLz033Dy1H//DnglIk6Sz/4rbt+AEfXfeAssvwncLmkPheHh/0vpC4CPgVkUpiF+R9I/Sa/dGBFfAv5lut3UhPY9RuEXrgv4A+CvgH4Kf1mUavQ2v5G2ERrYh5I+D/wJ8NsRcaxS1iHSokJ6JjJoH+Sj/8pWMURaM/qvktz0n6SLgLXArQNJQ2RrWv8N0T6oof/GVWCJiLcjYnFEXAY8TWEeGwprLC9ExC8iogf4S6Azlfn7dH8ceIpCEBrV9kVEf0T8x4i4JCKuBc4Guil8kM8pqmIO8F6j2ldjGxvWh5ImUvhP82RE/DQlvy9pZnp9JtCT0g8xeAQ10FcN68OM2peX/isnL/1XVl76T9IcYANwc0QMfPbkpv/KtK+m/htXgWVgN4OkCcBq4I/SSweBq1VwBnAF8Haa1jknlZkIfI3CVNCotk/SL6V2IWkR0B8Rb6ah7HFJV6Th6c0U5lIbZqRtbFQfpp/3UeCtiPhe0UsbgeXp8XI+7Y+NwDJJ7Wmqbh6wq1F9mFX7ctR/Q8pR/5WrJxf9J+ls4HkK62h/OZA5L/1Xrn0191/Wi0R5uVH4a/ow8AsKfxXcAtxBYWH+b4H7+XQR+vPA/wL2Am8C/ymln0Fhd9PfpNf+B2mnzii3r4PCottbwIsUrjQ6UE9n+od+B3hwoExe2tioPqSwAyhSva+m269T2HyxlcJoaStpE0Yq87upn/ZRtPOmEX2YVfty1n8HgD7go/T7MD9n/feZ9uWl/yj8EfaPRXlfBWbkpf/Kta/W/vPJezMzy9S4mgozM7PGc2AxM7NMObCYmVmmHFjMzCxTDixmZpYpBxYzM8uUA4uZmWXKgcXMzDL1/wERgRmeflUJDAAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.plot(style='*')" ] @@ -314,9 +2392,59 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 18, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "2021 743449\n", + "2014 1600941\n", + "1991 1659249\n", + "1995 1840410\n", + "2020 2010315\n", + "2022 2060304\n", + "2012 2175217\n", + "2003 2234584\n", + "2019 2254386\n", + "2006 2307352\n", + "2017 2321583\n", + "2001 2529279\n", + "1992 2574578\n", + "1993 2703886\n", + "2018 2705325\n", + "1988 2765617\n", + "2007 2780164\n", + "1987 2855570\n", + "2016 2856393\n", + "2011 2857040\n", + "2023 2873501\n", + "2008 2973918\n", + "1998 3034904\n", + "2002 3125418\n", + "2009 3444020\n", + "1994 3514763\n", + "1996 3539413\n", + "2004 3567744\n", + "1997 3620066\n", + "2015 3654892\n", + "2024 3670417\n", + "2000 3826372\n", + "2005 3835025\n", + "1999 3908112\n", + "2010 4111392\n", + "2013 4182691\n", + "1986 5115251\n", + "1990 5235827\n", + "1989 5466192\n", + "dtype: int64" + ] + }, + "execution_count": 18, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,9 +2459,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 19, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 19, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEKCAYAAAACS67iAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFJlJREFUeJzt3XmUZGV5x/Hv4wzCQLOZkUYHZNxiREZAGhWJ2IM74M4RFRA86miiQHA8CXJUooaIC+aoccm4Hxf6KJi4YFwCaXALMoPLiCh4YBQGAVl1gIgjT/5474RyYOjuquq+XW99P+fMmapbt+q+z1vVv7r3vffWjcxEkjT47tN2AyRJ/WGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiqxcC4Xtnjx4ly6dOlcLvLP3HrrrWy33XatLb9tw14/2AfWP5j1r1mz5vrMvP9U881poC9dupTVq1fP5SL/zOTkJOPj460tv23DXj/YB9Y/mPVHxK+mM59DLpJUCQNdkiphoEtSJQx0SaqEgS5JlZgy0CPi4xFxXUT8tGPa/SLiWxFxWfP/zrPbTEnSVKazhv5J4BmbTTsJOCczHw6c09yXJLVoykDPzPOBGzeb/BzgU83tTwHP7XO7JEkzFNO5pmhELAW+mpl7NfdvzsydOh6/KTPvcdglIlYAKwBGR0f3m5iY6EOzu7NhwwZGRkZaW37bBqH+tetvmdXXH10E195+9+nLluw4q8udLwbhMzCbBrX+5cuXr8nMsanmm/UzRTNzFbAKYGxsLNs8S2tQzxLrl0Go/9iTzp7V11+5bCOnr737x37dkeOzutz5YhA+A7Op9vq7Pcrl2oh4AEDz/3X9a5IkqRvdBvqXgWOa28cAX+pPcyRJ3ZrOYYtnAN8HHhERV0XEy4HTgKdGxGXAU5v7kqQWTTmGnpkv3sJDT+5zWyRJPfBMUUmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1Ilegr0iDgxIi6OiJ9GxBkRsU2/GiZJmpmuAz0ilgDHA2OZuRewAHhRvxomSZqZXodcFgKLImIhsC1wde9NkiR1IzKz+ydHnACcCtwOfDMzj7yHeVYAKwBGR0f3m5iY6Hp5vdqwYQMjIyOtLb9tg1D/2vW3zOrrjy6Ca2+/+/RlS3ac1eXOF4PwGZhNg1r/8uXL12Tm2FTzdR3oEbEzcBZwBHAz8AXgzMz8zJaeMzY2lqtXr+5qef0wOTnJ+Ph4a8tv2yDUv/Sks2f19Vcu28jpaxfebfq60w6d1eXOF4PwGZhNg1p/REwr0HsZcnkKcEVm/jYz/wh8EXhCD68nSepBL4H+a+DxEbFtRATwZOCS/jRLkjRTXQd6Zl4AnAlcBKxtXmtVn9olSZqhuw8mzkBmngKc0qe2SJJ64JmiklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVaKn30NXnWb7up6SZodr6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklSJngI9InaKiDMj4ucRcUlEHNCvhkmSZqbXS9C9F/h6Zh4eEfcFtu1DmyRJXeg60CNiB+Ag4FiAzLwDuKM/zZIkzVRkZndPjNgHWAX8DNgbWAOckJm3bjbfCmAFwOjo6H4TExM9NbgXGzZsYGRkpLXlt2269a9df8sctKYdo4vg2tvvPn3Zkh3nvjGNuezvzvrbrLktg5oBy5cvX5OZY1PN10ugjwH/AxyYmRdExHuB32Xmm7b0nLGxsVy9enVXy+uHyclJxsfHW1t+26Zb/9KTzp79xrRk5bKNnL727hum6047tIXWFHPZ3531t1lzWwY1AyJiWoHey07Rq4CrMvOC5v6ZwGN6eD1JUg+6DvTMvAa4MiIe0Ux6MmX4RZLUgl6PcjkO+GxzhMvlwMt6b5IkqRs9BXpm/giYclxHkjT7PFNUkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekSvR6gQupCjVfR1XDwzV0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1Ileg70iFgQET+MiK/2o0GSpO70Yw39BOCSPryOJKkHPQV6ROwGHAp8tD/NkSR1KzKz+ydHnAm8HdgeeH1mHnYP86wAVgCMjo7uNzEx0fXyerVhwwZGRkZaW37bplv/2vW3zEFr2jG6CK69ve1WtKez/mVLdmy3MS0Y1AxYvnz5mswcm2q+ri8SHRGHAddl5pqIGN/SfJm5ClgFMDY2luPjW5x11k1OTtLm8ts23fqPrfiCySuXbeT0tcN7bfTO+tcdOd5uY1pQewb0MuRyIPDsiFgHTAAHR8Rn+tIqSdKMdR3omfmGzNwtM5cCLwLOzcyj+tYySdKMeBy6JFWiL4OJmTkJTPbjtSRJ3XENXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoM78UVB8DSPl/bc+WyjVVfL1Saytr1t7T2N7DutENnfRmuoUtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1Ilug70iNg9Iv47Ii6JiIsj4oR+NkySNDO9XIJuI7AyMy+KiO2BNRHxrcz8WZ/aJkmaga7X0DPzN5l5UXP798AlwJJ+NUySNDORmb2/SMRS4Hxgr8z83WaPrQBWAIyOju43MTHR8/K6tWHDBkZGRmb8vLXrb5mF1sy90UVw7e1tt6Jdw94H86X+ZUt2bGW51914S2v191Lz8uXL12Tm2FTz9RzoETECnAecmplfvLd5x8bGcvXq1T0trxeTk5OMj4/P+HlLW7pKeL+tXLaR09f2Mso2+Ia9D+ZL/etOO7SV5b7/s19qrf5eao6IaQV6T0e5RMRWwFnAZ6cKc0nS7OrlKJcAPgZckpnv6V+TJEnd6GUN/UDgaODgiPhR8++QPrVLkjRDXQ8mZeZ3gOhjWyRJPfBMUUmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUifavRTVN/bgM3MplGzm2ksvJSdLmXEOXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekSvQU6BHxjIj4RUT8MiJO6lejJEkz13WgR8QC4APAM4E9gRdHxJ79apgkaWZ6WUN/LPDLzLw8M+8AJoDn9KdZkqSZiszs7okRhwPPyMxXNPePBh6Xma/dbL4VwIrm7iOAX3Tf3J4tBq5vcfltG/b6wT6w/sGsf4/MvP9UMy3sYQFxD9Pu9u2QmauAVT0sp28iYnVmjrXdjrYMe/1gH1h/3fX3MuRyFbB7x/3dgKt7a44kqVu9BPqFwMMj4sERcV/gRcCX+9MsSdJMdT3kkpkbI+K1wDeABcDHM/PivrVsdsyLoZ8WDXv9YB9Yf8W63ikqSZpfPFNUkiphoEtSJQx0SaqEgT5NEbEkIpa03Y62RMRDIuLEiDi47ba0YdjrB/tgEOo30KcQEUsj4jzg68C7IuKJbbdprkXEXwPfovxmz6sj4m9abtKcGvb6wT4YlPoN9HsQEdt03H0McGFmLqMconl8RCxr5runs2UHXkQcHBEPbm4HcDBwSma+EjgdOCQixjser8qw1w/2waDWb6A3ImKHiPhwRFwKvDsi9mgeeh7w6+b2BPBL4BWbnjbHzZxVEbFnRPwE+EfgExFxcJbjWvcEdgXIzAuA7wEv2/S0Nto6G4a9frAPBr1+A/0uzwC2obxxdwBvjohFlM2sZwFk5h+AM4EnNvfvbKep/RERu0XEDh2TjgDOysyDKF9eL4mIhwOfo+mDxr8De0XE1oPcB8NeP9gHtdU/VIEexcKIeHlEfDsiToiIhzYPPwy4IzM3Av8C3AQcBXwTeEBE3K+Z71Lgyog4YM4L6JOIeGREfA34DvDWiNj0s8f/C2zb3P48cA1wKGVt5C86tlpupPxq5t5z1+r+Gfb6wT6otf6hCvRm0+lJwEuBdwJbAx9pHr4GuK75xr2S8mY9lPIG/4y7fgJ4K+CGZvrAiIjtOu7uA1yVmUuBc4F3N9NvBP4QEdtn5o3AZcADKbV+D3hdM999gT8B62a/5f0x7PWDfTAM9Vcd6BFxQES8IyKObe4H8Ejg65n5lcx8J7BHRDwBWE/5Zn5k8/RLgJFm2r9SdoI8i/JlMAr8eE6L6UJE7BwRn4yIC4HTIuL+TR8sA74bEZGZXwZujohDKVsf2zeP09zfBbiTstWyS0R8BDgD2JiZ1811TTMx7PWDfTBs9Vcb6BHxKOBDwO+BF0bE6yj1LgF+37Fn+pPASygBvRF4QjP9Isqe7dsy83zgJOBY4EDgbZl553zau70FB1FqOoSy4+ZkYAfKj6ntmnf9kM+nKH3wA0p/PRMgM7/fvMbCzLwEeBVwMfDPmfky5r9hrx/sg+GqPzMH/h9lTfoVlM2ohc209wAnNLfHgPcBhwNPAb7R8dzdKZteUAL8h5QrK+0LfAl4QMe80XatW6h/AeWDdh5laGhxM/3zwPHN7QcDpzWP708ZO1zQ0X+/bV5nCWXr5LXAJ4APAtu1XaP12wfWP/W/gV9Dj4i9KTsunwOcAryxeWg95bqnUL5Rvwu8ADgH2DUiHh0RW2UZL18fEU/MzHMpP6/5DuCLwBmZ+ZtNy8rmnZ+HDgOeDbwFOICyfwDKETqbtjiuBL4NPDMzL6SsrSwHyMwNwAXA/pm5HjiaMqx0DfDGzLx1juro1rDXD/bBsNcP9HYJulZExL7ArZl5aTPpscClmXlsRDwGODUixoBJ4OkRsW1m3hYRPwZeSDmW9HPAK4H3RcTtwFrgiub1Pgx8LjNvmbuqptaM9WVE7E/ZNPw2cHaWQyn/Erg8M8+NiCsoZ7Q+DVgDPC8iFmfm9RFxGXBrRDwIeD9wVETsQrna1A2UzU0yczWwes6LvBfDXj/YB8Ne/3QMzBp6RDwsIr5DGfN+c0S8tHnoTmBds7Z9EWVT6QDgNu465Ajgj5RNql0pa+E/bV7rPOD6zLwKylr4PA7zg4CPU/a4PwV4ezPLncClEbEoM6+g9MGjKWOBV1OOrYWyV34B5X0/i9IPRwL7AatyHh1P2ykiFjT1P4myCTxU9QM0tWWUsxOH8TOwwzDXP21tj/ls6R+wHXBAx/3DgPc2tx9H+fbcAziGspm1pHnscMp4+abHzmumb0MZblnc8Zr7Avdtu9Yt1L8t8Gru2prYCvg74DXN4zsDP2lqOIIyNri0o69WUa5wfhhlC2RHyj6Cr3XWDNyn7Vrv5f1/BeWPbiVlR9bQ1N/Rvu2BsylXBAM4cVj6oPkbOKb5uz1r2Orv5t+8XEOPiJOBy4GvRcRoM/nplGPDyXLq7Q+A4yjHkO5OOWYcylj5PpSjUz4F3BQRn6bs7PwF8P9jYZn5w8y8Y/YrmpmI2BX4KjAOfJqyE+f5lC2PjQCZeRNlp+3xlHHCXbjrkMvzKcfb35GZXwU+RjnD9QOUvfl/3LSsnIdrJM3xwudQ/vg+AjyNsv9jf8qaWNX1b2YR5XyJh0bEYsrnfAHU3QcRsRVl39fhwLsy8wXNQ/tumqfm+rvW9jfKFr6ZxymbSx8FVjbTXkcZL9s0z57Alc3tU4G3djx2IbBvc3tryiFL+7dd1wzqXwQ8ruP+sZSdPMcAP+iY/kDg6ub2ayinKu/cPP8rwIM65l08F23vYx/s1HH77yl/tEcOS/0d7T4GeBfwJuDllNPPLxyGPqAcmHDkZtOOAC4Yhvq76rO2G7CFN3LToURHcNeQyU7AzcA2HfNdSPnG3gn4AmUT6z8p38Bbt11HD/XHpn/N/cd09MMNlONnN837rU3hD/wT5YifG4B/aLuOPvTDDpT9HNcCb23u3wCM1l5/x3v/MsrQ2/OBzzbTrh+SPjiMcmLP6ZSDHN5MGUq9Edil9vq7+Tcvh1wy80/Nzf8Cto+IvTPzZsq4+as6Zr0I2L557DjKsMp/ACuy7PkeSNnomHQCZW0FyvjfiQBRfl/mV8CmQyvfQtmSWZKZ75ij5s6azPwdZWjt8ZSd2S+kDJu9Kopq6+94/w+hDDudA+wWEW+k7OxfAXV/BrIMlVxBCeejgUcBz6V8Bl5d+2egG/HnuTH/RMQHKePhr2+O8vhbSrDvTDlh6JCOL4DqRMRulPG/4zLz0ig/JraC8uFeAvwo5+MZa30WEftQvsy/Txkn3YtyqFm19UfECGW4ZWtKzX9FORnmZMqa+8Opvw+2zczbmtt7Uz7736Wcml/9Z2CmBiHQ96EctfJkyof6Nspp+LcDH87Mi1ts3qyL8vsxTwLeQBlPvYqyiXkE8PMsh2pWLyJ2p3yxvTgzb4iIo4CLM/OHLTdt1kS50Mq/UXbgnUE55O7kzHxa83j1fdApyi8dfhQ4IjNvHLb6p2MQAv1FlEP3bgPeRtnjXc9e6SlExHeBh1B+1e1q4C2Z+ZNWGzVHImJHyhf5Syg7wVcBH8jMP97rEyvVnAzzfGAiM69puz1zISK2plyrYNOQy4eAD2b5mWttZl4HekQ8mnIa/pmUHUID9ZO1vWoO3TqFMo74mUHeL9CNiFhIGWb5A6X+oXr/N4mIBcCdOZ//WGdRRLyKcrjqp4f1MzBd8zrQJUnTNy+PcpEkzZyBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkirxfxo0Erf7i0deAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.hist(xrot=20)" ] @@ -364,7 +2515,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4, -- 2.18.1