From 5570afb98811d885b341448a7f65da6da726f860 Mon Sep 17 00:00:00 2001 From: 8573a96fbb50caf6a42cf0daddd945a7 <8573a96fbb50caf6a42cf0daddd945a7@app-learninglab.inria.fr> Date: Wed, 9 Jul 2025 10:09:38 +0000 Subject: [PATCH] =?UTF-8?q?Lecture=20des=20donn=C3=A9es=20avec=20une=20cop?= =?UTF-8?q?ie=20locale=20des=20donn=C3=A9es.?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- module3/exo1/analyse-syndrome-grippal.ipynb | 3206 ++++++++++++++++++- 1 file changed, 3169 insertions(+), 37 deletions(-) diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5..540f7a7 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
020252632751321273.033753.04132.050.0FRFrance
120252532373718844.028630.03528.042.0FRFrance
220252432315418577.027731.03528.042.0FRFrance
320252332439119307.029475.03628.044.0FRFrance
420252231875514333.023177.02821.035.0FRFrance
520252132376018671.028849.03527.043.0FRFrance
620252032026515814.024716.03023.037.0FRFrance
720251931626412394.020134.02418.030.0FRFrance
820251831811513975.022255.02721.033.0FRFrance
920251732215017291.027009.03326.040.0FRFrance
1020251632856422550.034578.04334.052.0FRFrance
1120251533572129592.041850.05344.062.0FRFrance
1220251433757931232.043926.05647.065.0FRFrance
1320251333967333686.045660.05950.068.0FRFrance
1420251235254345627.059459.07868.088.0FRFrance
1520251135946952154.066784.08978.0100.0FRFrance
1620251036033453048.067620.09079.0101.0FRFrance
1720250938453174994.094068.0126112.0140.0FRFrance
182025083136020124824.0147216.0203186.0220.0FRFrance
192025073208952195988.0221916.0312293.0331.0FRFrance
202025063273519258159.0288879.0408385.0431.0FRFrance
212025053334395318416.0350374.0499475.0523.0FRFrance
222025043350043332885.0367201.0522496.0548.0FRFrance
232025033252772238917.0266627.0377356.0398.0FRFrance
242025023257247242991.0271503.0384363.0405.0FRFrance
252025013231549214627.0248471.0345320.0370.0FRFrance
262024523201726185870.0217582.0302278.0326.0FRFrance
272024513201697187843.0215551.0302281.0323.0FRFrance
282024503136694126369.0147019.0205190.0220.0FRFrance
29202449310848799037.0117937.0163149.0177.0FRFrance
.................................
209219852132609619621.032571.04735.059.0FRFrance
209319852032789620885.034907.05138.064.0FRFrance
209419851934315432821.053487.07859.097.0FRFrance
209519851834055529935.051175.07455.093.0FRFrance
209619851733405324366.043740.06244.080.0FRFrance
209719851635036236451.064273.09166.0116.0FRFrance
209819851536388145538.082224.011683.0149.0FRFrance
20991985143134545114400.0154690.0244207.0281.0FRFrance
21001985133197206176080.0218332.0357319.0395.0FRFrance
21011985123245240223304.0267176.0445405.0485.0FRFrance
21021985113276205252399.0300011.0501458.0544.0FRFrance
21031985103353231326279.0380183.0640591.0689.0FRFrance
21041985093369895341109.0398681.0670618.0722.0FRFrance
21051985083389886359529.0420243.0707652.0762.0FRFrance
21061985073471852432599.0511105.0855784.0926.0FRFrance
21071985063565825518011.0613639.01026939.01113.0FRFrance
21081985053637302592795.0681809.011551074.01236.0FRFrance
21091985043424937390794.0459080.0770708.0832.0FRFrance
21101985033213901174689.0253113.0388317.0459.0FRFrance
211119850239758680949.0114223.0177147.0207.0FRFrance
211219850138548965918.0105060.0155120.0190.0FRFrance
211319845238483060602.0109058.0154110.0198.0FRFrance
2114198451310172680242.0123210.0185146.0224.0FRFrance
21151984503123680101401.0145959.0225184.0266.0FRFrance
2116198449310107381684.0120462.0184149.0219.0FRFrance
211719844837862060634.096606.0143110.0176.0FRFrance
211819844737202954274.089784.013199.0163.0FRFrance
211919844638733067686.0106974.0159123.0195.0FRFrance
21201984453135223101414.0169032.0246184.0308.0FRFrance
212119844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2122 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202526 3 27513 21273.0 33753.0 41 32.0 \n", + "1 202525 3 23737 18844.0 28630.0 35 28.0 \n", + "2 202524 3 23154 18577.0 27731.0 35 28.0 \n", + "3 202523 3 24391 19307.0 29475.0 36 28.0 \n", + "4 202522 3 18755 14333.0 23177.0 28 21.0 \n", + "5 202521 3 23760 18671.0 28849.0 35 27.0 \n", + "6 202520 3 20265 15814.0 24716.0 30 23.0 \n", + "7 202519 3 16264 12394.0 20134.0 24 18.0 \n", + "8 202518 3 18115 13975.0 22255.0 27 21.0 \n", + "9 202517 3 22150 17291.0 27009.0 33 26.0 \n", + "10 202516 3 28564 22550.0 34578.0 43 34.0 \n", + "11 202515 3 35721 29592.0 41850.0 53 44.0 \n", + "12 202514 3 37579 31232.0 43926.0 56 47.0 \n", + "13 202513 3 39673 33686.0 45660.0 59 50.0 \n", + "14 202512 3 52543 45627.0 59459.0 78 68.0 \n", + "15 202511 3 59469 52154.0 66784.0 89 78.0 \n", + "16 202510 3 60334 53048.0 67620.0 90 79.0 \n", + "17 202509 3 84531 74994.0 94068.0 126 112.0 \n", + "18 202508 3 136020 124824.0 147216.0 203 186.0 \n", + "19 202507 3 208952 195988.0 221916.0 312 293.0 \n", + "20 202506 3 273519 258159.0 288879.0 408 385.0 \n", + "21 202505 3 334395 318416.0 350374.0 499 475.0 \n", + "22 202504 3 350043 332885.0 367201.0 522 496.0 \n", + "23 202503 3 252772 238917.0 266627.0 377 356.0 \n", + "24 202502 3 257247 242991.0 271503.0 384 363.0 \n", + "25 202501 3 231549 214627.0 248471.0 345 320.0 \n", + "26 202452 3 201726 185870.0 217582.0 302 278.0 \n", + "27 202451 3 201697 187843.0 215551.0 302 281.0 \n", + "28 202450 3 136694 126369.0 147019.0 205 190.0 \n", + "29 202449 3 108487 99037.0 117937.0 163 149.0 \n", + "... ... ... ... ... ... ... ... \n", + "2092 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "2093 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "2094 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "2095 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "2096 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "2097 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2098 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2099 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2100 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2101 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2102 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2103 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2104 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2105 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2106 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2107 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2108 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2109 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2110 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2111 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2112 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2113 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2114 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2115 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2116 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2117 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2118 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2119 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2120 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2121 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 50.0 FR France \n", + "1 42.0 FR France \n", + "2 42.0 FR France \n", + "3 44.0 FR France \n", + "4 35.0 FR France \n", + "5 43.0 FR France \n", + "6 37.0 FR France \n", + "7 30.0 FR France \n", + "8 33.0 FR France \n", + "9 40.0 FR France \n", + "10 52.0 FR France \n", + "11 62.0 FR France \n", + "12 65.0 FR France \n", + "13 68.0 FR France \n", + "14 88.0 FR France \n", + "15 100.0 FR France \n", + "16 101.0 FR France \n", + "17 140.0 FR France \n", + "18 220.0 FR France \n", + "19 331.0 FR France \n", + "20 431.0 FR France \n", + "21 523.0 FR France \n", + "22 548.0 FR France \n", + "23 398.0 FR France \n", + "24 405.0 FR France \n", + "25 370.0 FR France \n", + "26 326.0 FR France \n", + "27 323.0 FR France \n", + "28 220.0 FR France \n", + "29 177.0 FR France \n", + "... ... ... ... \n", + "2092 59.0 FR France \n", + "2093 64.0 FR France \n", + "2094 97.0 FR France \n", + "2095 93.0 FR France \n", + "2096 80.0 FR France \n", + "2097 116.0 FR France \n", + "2098 149.0 FR France \n", + "2099 281.0 FR France \n", + "2100 395.0 FR France \n", + "2101 485.0 FR France \n", + "2102 544.0 FR France \n", + "2103 689.0 FR France \n", + "2104 722.0 FR France \n", + "2105 762.0 FR France \n", + "2106 926.0 FR France \n", + "2107 1113.0 FR France \n", + "2108 1236.0 FR France \n", + "2109 832.0 FR France \n", + "2110 459.0 FR France \n", + "2111 207.0 FR France \n", + "2112 190.0 FR France \n", + "2113 198.0 FR France \n", + "2114 224.0 FR France \n", + "2115 266.0 FR France \n", + "2116 219.0 FR France \n", + "2117 176.0 FR France \n", + "2118 163.0 FR France \n", + "2119 195.0 FR France \n", + "2120 308.0 FR France \n", + "2121 213.0 FR France \n", + "\n", + "[2122 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
18851989193-NaNNaN-NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1885 198919 3 - NaN NaN - NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1885 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
020252632751321273.033753.04132.050.0FRFrance
120252532373718844.028630.03528.042.0FRFrance
220252432315418577.027731.03528.042.0FRFrance
320252332439119307.029475.03628.044.0FRFrance
420252231875514333.023177.02821.035.0FRFrance
520252132376018671.028849.03527.043.0FRFrance
620252032026515814.024716.03023.037.0FRFrance
720251931626412394.020134.02418.030.0FRFrance
820251831811513975.022255.02721.033.0FRFrance
920251732215017291.027009.03326.040.0FRFrance
1020251632856422550.034578.04334.052.0FRFrance
1120251533572129592.041850.05344.062.0FRFrance
1220251433757931232.043926.05647.065.0FRFrance
1320251333967333686.045660.05950.068.0FRFrance
1420251235254345627.059459.07868.088.0FRFrance
1520251135946952154.066784.08978.0100.0FRFrance
1620251036033453048.067620.09079.0101.0FRFrance
1720250938453174994.094068.0126112.0140.0FRFrance
182025083136020124824.0147216.0203186.0220.0FRFrance
192025073208952195988.0221916.0312293.0331.0FRFrance
202025063273519258159.0288879.0408385.0431.0FRFrance
212025053334395318416.0350374.0499475.0523.0FRFrance
222025043350043332885.0367201.0522496.0548.0FRFrance
232025033252772238917.0266627.0377356.0398.0FRFrance
242025023257247242991.0271503.0384363.0405.0FRFrance
252025013231549214627.0248471.0345320.0370.0FRFrance
262024523201726185870.0217582.0302278.0326.0FRFrance
272024513201697187843.0215551.0302281.0323.0FRFrance
282024503136694126369.0147019.0205190.0220.0FRFrance
29202449310848799037.0117937.0163149.0177.0FRFrance
.................................
209219852132609619621.032571.04735.059.0FRFrance
209319852032789620885.034907.05138.064.0FRFrance
209419851934315432821.053487.07859.097.0FRFrance
209519851834055529935.051175.07455.093.0FRFrance
209619851733405324366.043740.06244.080.0FRFrance
209719851635036236451.064273.09166.0116.0FRFrance
209819851536388145538.082224.011683.0149.0FRFrance
20991985143134545114400.0154690.0244207.0281.0FRFrance
21001985133197206176080.0218332.0357319.0395.0FRFrance
21011985123245240223304.0267176.0445405.0485.0FRFrance
21021985113276205252399.0300011.0501458.0544.0FRFrance
21031985103353231326279.0380183.0640591.0689.0FRFrance
21041985093369895341109.0398681.0670618.0722.0FRFrance
21051985083389886359529.0420243.0707652.0762.0FRFrance
21061985073471852432599.0511105.0855784.0926.0FRFrance
21071985063565825518011.0613639.01026939.01113.0FRFrance
21081985053637302592795.0681809.011551074.01236.0FRFrance
21091985043424937390794.0459080.0770708.0832.0FRFrance
21101985033213901174689.0253113.0388317.0459.0FRFrance
211119850239758680949.0114223.0177147.0207.0FRFrance
211219850138548965918.0105060.0155120.0190.0FRFrance
211319845238483060602.0109058.0154110.0198.0FRFrance
2114198451310172680242.0123210.0185146.0224.0FRFrance
21151984503123680101401.0145959.0225184.0266.0FRFrance
2116198449310107381684.0120462.0184149.0219.0FRFrance
211719844837862060634.096606.0143110.0176.0FRFrance
211819844737202954274.089784.013199.0163.0FRFrance
211919844638733067686.0106974.0159123.0195.0FRFrance
21201984453135223101414.0169032.0246184.0308.0FRFrance
212119844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2121 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202526 3 27513 21273.0 33753.0 41 32.0 \n", + "1 202525 3 23737 18844.0 28630.0 35 28.0 \n", + "2 202524 3 23154 18577.0 27731.0 35 28.0 \n", + "3 202523 3 24391 19307.0 29475.0 36 28.0 \n", + "4 202522 3 18755 14333.0 23177.0 28 21.0 \n", + "5 202521 3 23760 18671.0 28849.0 35 27.0 \n", + "6 202520 3 20265 15814.0 24716.0 30 23.0 \n", + "7 202519 3 16264 12394.0 20134.0 24 18.0 \n", + "8 202518 3 18115 13975.0 22255.0 27 21.0 \n", + "9 202517 3 22150 17291.0 27009.0 33 26.0 \n", + "10 202516 3 28564 22550.0 34578.0 43 34.0 \n", + "11 202515 3 35721 29592.0 41850.0 53 44.0 \n", + "12 202514 3 37579 31232.0 43926.0 56 47.0 \n", + "13 202513 3 39673 33686.0 45660.0 59 50.0 \n", + "14 202512 3 52543 45627.0 59459.0 78 68.0 \n", + "15 202511 3 59469 52154.0 66784.0 89 78.0 \n", + "16 202510 3 60334 53048.0 67620.0 90 79.0 \n", + "17 202509 3 84531 74994.0 94068.0 126 112.0 \n", + "18 202508 3 136020 124824.0 147216.0 203 186.0 \n", + "19 202507 3 208952 195988.0 221916.0 312 293.0 \n", + "20 202506 3 273519 258159.0 288879.0 408 385.0 \n", + "21 202505 3 334395 318416.0 350374.0 499 475.0 \n", + "22 202504 3 350043 332885.0 367201.0 522 496.0 \n", + "23 202503 3 252772 238917.0 266627.0 377 356.0 \n", + "24 202502 3 257247 242991.0 271503.0 384 363.0 \n", + "25 202501 3 231549 214627.0 248471.0 345 320.0 \n", + "26 202452 3 201726 185870.0 217582.0 302 278.0 \n", + "27 202451 3 201697 187843.0 215551.0 302 281.0 \n", + "28 202450 3 136694 126369.0 147019.0 205 190.0 \n", + "29 202449 3 108487 99037.0 117937.0 163 149.0 \n", + "... ... ... ... ... ... ... ... \n", + "2092 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "2093 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "2094 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "2095 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "2096 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "2097 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2098 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2099 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2100 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2101 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2102 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2103 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2104 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2105 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2106 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2107 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2108 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2109 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2110 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2111 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2112 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2113 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2114 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2115 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2116 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2117 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2118 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2119 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2120 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2121 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 50.0 FR France \n", + "1 42.0 FR France \n", + "2 42.0 FR France \n", + "3 44.0 FR France \n", + "4 35.0 FR France \n", + "5 43.0 FR France \n", + "6 37.0 FR France \n", + "7 30.0 FR France \n", + "8 33.0 FR France \n", + "9 40.0 FR France \n", + "10 52.0 FR France \n", + "11 62.0 FR France \n", + "12 65.0 FR France \n", + "13 68.0 FR France \n", + "14 88.0 FR France \n", + "15 100.0 FR France \n", + "16 101.0 FR France \n", + "17 140.0 FR France \n", + "18 220.0 FR France \n", + "19 331.0 FR France \n", + "20 431.0 FR France \n", + "21 523.0 FR France \n", + "22 548.0 FR France \n", + "23 398.0 FR France \n", + "24 405.0 FR France \n", + "25 370.0 FR France \n", + "26 326.0 FR France \n", + "27 323.0 FR France \n", + "28 220.0 FR France \n", + "29 177.0 FR France \n", + "... ... ... ... \n", + "2092 59.0 FR France \n", + "2093 64.0 FR France \n", + "2094 97.0 FR France \n", + "2095 93.0 FR France \n", + "2096 80.0 FR France \n", + "2097 116.0 FR France \n", + "2098 149.0 FR France \n", + "2099 281.0 FR France \n", + "2100 395.0 FR France \n", + "2101 485.0 FR France \n", + "2102 544.0 FR France \n", + "2103 689.0 FR France \n", + "2104 722.0 FR France \n", + "2105 762.0 FR France \n", + "2106 926.0 FR France \n", + "2107 1113.0 FR France \n", + "2108 1236.0 FR France \n", + "2109 832.0 FR France \n", + "2110 459.0 FR France \n", + "2111 207.0 FR France \n", + "2112 190.0 FR France \n", + "2113 198.0 FR France \n", + "2114 224.0 FR France \n", + "2115 266.0 FR France \n", + "2116 219.0 FR France \n", + "2117 176.0 FR France \n", + "2118 163.0 FR France \n", + "2119 195.0 FR France \n", + "2120 308.0 FR France \n", + "2121 213.0 FR France \n", + "\n", + "[2121 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,15 +2148,1004 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 7, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" ] }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "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", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
period
1984-10-29/1984-11-0419844436842220056.0116788.012537.0213.0FRFrance
1984-11-05/1984-11-111984453135223101414.0169032.0246184.0308.0FRFrance
1984-11-12/1984-11-1819844638733067686.0106974.0159123.0195.0FRFrance
1984-11-19/1984-11-2519844737202954274.089784.013199.0163.0FRFrance
1984-11-26/1984-12-0219844837862060634.096606.0143110.0176.0FRFrance
1984-12-03/1984-12-09198449310107381684.0120462.0184149.0219.0FRFrance
1984-12-10/1984-12-161984503123680101401.0145959.0225184.0266.0FRFrance
1984-12-17/1984-12-23198451310172680242.0123210.0185146.0224.0FRFrance
1984-12-24/1984-12-3019845238483060602.0109058.0154110.0198.0FRFrance
1984-12-31/1985-01-0619850138548965918.0105060.0155120.0190.0FRFrance
1985-01-07/1985-01-1319850239758680949.0114223.0177147.0207.0FRFrance
1985-01-14/1985-01-201985033213901174689.0253113.0388317.0459.0FRFrance
1985-01-21/1985-01-271985043424937390794.0459080.0770708.0832.0FRFrance
1985-01-28/1985-02-031985053637302592795.0681809.011551074.01236.0FRFrance
1985-02-04/1985-02-101985063565825518011.0613639.01026939.01113.0FRFrance
1985-02-11/1985-02-171985073471852432599.0511105.0855784.0926.0FRFrance
1985-02-18/1985-02-241985083389886359529.0420243.0707652.0762.0FRFrance
1985-02-25/1985-03-031985093369895341109.0398681.0670618.0722.0FRFrance
1985-03-04/1985-03-101985103353231326279.0380183.0640591.0689.0FRFrance
1985-03-11/1985-03-171985113276205252399.0300011.0501458.0544.0FRFrance
1985-03-18/1985-03-241985123245240223304.0267176.0445405.0485.0FRFrance
1985-03-25/1985-03-311985133197206176080.0218332.0357319.0395.0FRFrance
1985-04-01/1985-04-071985143134545114400.0154690.0244207.0281.0FRFrance
1985-04-08/1985-04-1419851536388145538.082224.011683.0149.0FRFrance
1985-04-15/1985-04-2119851635036236451.064273.09166.0116.0FRFrance
1985-04-22/1985-04-2819851733405324366.043740.06244.080.0FRFrance
1985-04-29/1985-05-0519851834055529935.051175.07455.093.0FRFrance
1985-05-06/1985-05-1219851934315432821.053487.07859.097.0FRFrance
1985-05-13/1985-05-1919852032789620885.034907.05138.064.0FRFrance
1985-05-20/1985-05-2619852132609619621.032571.04735.059.0FRFrance
.................................
2024-12-02/2024-12-08202449310848799037.0117937.0163149.0177.0FRFrance
2024-12-09/2024-12-152024503136694126369.0147019.0205190.0220.0FRFrance
2024-12-16/2024-12-222024513201697187843.0215551.0302281.0323.0FRFrance
2024-12-23/2024-12-292024523201726185870.0217582.0302278.0326.0FRFrance
2024-12-30/2025-01-052025013231549214627.0248471.0345320.0370.0FRFrance
2025-01-06/2025-01-122025023257247242991.0271503.0384363.0405.0FRFrance
2025-01-13/2025-01-192025033252772238917.0266627.0377356.0398.0FRFrance
2025-01-20/2025-01-262025043350043332885.0367201.0522496.0548.0FRFrance
2025-01-27/2025-02-022025053334395318416.0350374.0499475.0523.0FRFrance
2025-02-03/2025-02-092025063273519258159.0288879.0408385.0431.0FRFrance
2025-02-10/2025-02-162025073208952195988.0221916.0312293.0331.0FRFrance
2025-02-17/2025-02-232025083136020124824.0147216.0203186.0220.0FRFrance
2025-02-24/2025-03-0220250938453174994.094068.0126112.0140.0FRFrance
2025-03-03/2025-03-0920251036033453048.067620.09079.0101.0FRFrance
2025-03-10/2025-03-1620251135946952154.066784.08978.0100.0FRFrance
2025-03-17/2025-03-2320251235254345627.059459.07868.088.0FRFrance
2025-03-24/2025-03-3020251333967333686.045660.05950.068.0FRFrance
2025-03-31/2025-04-0620251433757931232.043926.05647.065.0FRFrance
2025-04-07/2025-04-1320251533572129592.041850.05344.062.0FRFrance
2025-04-14/2025-04-2020251632856422550.034578.04334.052.0FRFrance
2025-04-21/2025-04-2720251732215017291.027009.03326.040.0FRFrance
2025-04-28/2025-05-0420251831811513975.022255.02721.033.0FRFrance
2025-05-05/2025-05-1120251931626412394.020134.02418.030.0FRFrance
2025-05-12/2025-05-1820252032026515814.024716.03023.037.0FRFrance
2025-05-19/2025-05-2520252132376018671.028849.03527.043.0FRFrance
2025-05-26/2025-06-0120252231875514333.023177.02821.035.0FRFrance
2025-06-02/2025-06-0820252332439119307.029475.03628.044.0FRFrance
2025-06-09/2025-06-1520252432315418577.027731.03528.042.0FRFrance
2025-06-16/2025-06-2220252532373718844.028630.03528.042.0FRFrance
2025-06-23/2025-06-2920252632751321273.033753.04132.050.0FRFrance
\n", + "

2121 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 \\\n", + "period \n", + "1984-10-29/1984-11-04 198444 3 68422 20056.0 116788.0 125 \n", + "1984-11-05/1984-11-11 198445 3 135223 101414.0 169032.0 246 \n", + "1984-11-12/1984-11-18 198446 3 87330 67686.0 106974.0 159 \n", + "1984-11-19/1984-11-25 198447 3 72029 54274.0 89784.0 131 \n", + "1984-11-26/1984-12-02 198448 3 78620 60634.0 96606.0 143 \n", + "1984-12-03/1984-12-09 198449 3 101073 81684.0 120462.0 184 \n", + "1984-12-10/1984-12-16 198450 3 123680 101401.0 145959.0 225 \n", + "1984-12-17/1984-12-23 198451 3 101726 80242.0 123210.0 185 \n", + "1984-12-24/1984-12-30 198452 3 84830 60602.0 109058.0 154 \n", + "1984-12-31/1985-01-06 198501 3 85489 65918.0 105060.0 155 \n", + "1985-01-07/1985-01-13 198502 3 97586 80949.0 114223.0 177 \n", + "1985-01-14/1985-01-20 198503 3 213901 174689.0 253113.0 388 \n", + "1985-01-21/1985-01-27 198504 3 424937 390794.0 459080.0 770 \n", + "1985-01-28/1985-02-03 198505 3 637302 592795.0 681809.0 1155 \n", + "1985-02-04/1985-02-10 198506 3 565825 518011.0 613639.0 1026 \n", + "1985-02-11/1985-02-17 198507 3 471852 432599.0 511105.0 855 \n", + "1985-02-18/1985-02-24 198508 3 389886 359529.0 420243.0 707 \n", + "1985-02-25/1985-03-03 198509 3 369895 341109.0 398681.0 670 \n", + "1985-03-04/1985-03-10 198510 3 353231 326279.0 380183.0 640 \n", + "1985-03-11/1985-03-17 198511 3 276205 252399.0 300011.0 501 \n", + "1985-03-18/1985-03-24 198512 3 245240 223304.0 267176.0 445 \n", + "1985-03-25/1985-03-31 198513 3 197206 176080.0 218332.0 357 \n", + "1985-04-01/1985-04-07 198514 3 134545 114400.0 154690.0 244 \n", + "1985-04-08/1985-04-14 198515 3 63881 45538.0 82224.0 116 \n", + "1985-04-15/1985-04-21 198516 3 50362 36451.0 64273.0 91 \n", + "1985-04-22/1985-04-28 198517 3 34053 24366.0 43740.0 62 \n", + "1985-04-29/1985-05-05 198518 3 40555 29935.0 51175.0 74 \n", + "1985-05-06/1985-05-12 198519 3 43154 32821.0 53487.0 78 \n", + "1985-05-13/1985-05-19 198520 3 27896 20885.0 34907.0 51 \n", + "1985-05-20/1985-05-26 198521 3 26096 19621.0 32571.0 47 \n", + "... ... ... ... ... ... ... \n", + "2024-12-02/2024-12-08 202449 3 108487 99037.0 117937.0 163 \n", + "2024-12-09/2024-12-15 202450 3 136694 126369.0 147019.0 205 \n", + "2024-12-16/2024-12-22 202451 3 201697 187843.0 215551.0 302 \n", + "2024-12-23/2024-12-29 202452 3 201726 185870.0 217582.0 302 \n", + "2024-12-30/2025-01-05 202501 3 231549 214627.0 248471.0 345 \n", + "2025-01-06/2025-01-12 202502 3 257247 242991.0 271503.0 384 \n", + "2025-01-13/2025-01-19 202503 3 252772 238917.0 266627.0 377 \n", + "2025-01-20/2025-01-26 202504 3 350043 332885.0 367201.0 522 \n", + "2025-01-27/2025-02-02 202505 3 334395 318416.0 350374.0 499 \n", + "2025-02-03/2025-02-09 202506 3 273519 258159.0 288879.0 408 \n", + "2025-02-10/2025-02-16 202507 3 208952 195988.0 221916.0 312 \n", + "2025-02-17/2025-02-23 202508 3 136020 124824.0 147216.0 203 \n", + "2025-02-24/2025-03-02 202509 3 84531 74994.0 94068.0 126 \n", + "2025-03-03/2025-03-09 202510 3 60334 53048.0 67620.0 90 \n", + "2025-03-10/2025-03-16 202511 3 59469 52154.0 66784.0 89 \n", + "2025-03-17/2025-03-23 202512 3 52543 45627.0 59459.0 78 \n", + "2025-03-24/2025-03-30 202513 3 39673 33686.0 45660.0 59 \n", + "2025-03-31/2025-04-06 202514 3 37579 31232.0 43926.0 56 \n", + "2025-04-07/2025-04-13 202515 3 35721 29592.0 41850.0 53 \n", + "2025-04-14/2025-04-20 202516 3 28564 22550.0 34578.0 43 \n", + "2025-04-21/2025-04-27 202517 3 22150 17291.0 27009.0 33 \n", + "2025-04-28/2025-05-04 202518 3 18115 13975.0 22255.0 27 \n", + "2025-05-05/2025-05-11 202519 3 16264 12394.0 20134.0 24 \n", + "2025-05-12/2025-05-18 202520 3 20265 15814.0 24716.0 30 \n", + "2025-05-19/2025-05-25 202521 3 23760 18671.0 28849.0 35 \n", + "2025-05-26/2025-06-01 202522 3 18755 14333.0 23177.0 28 \n", + "2025-06-02/2025-06-08 202523 3 24391 19307.0 29475.0 36 \n", + "2025-06-09/2025-06-15 202524 3 23154 18577.0 27731.0 35 \n", + "2025-06-16/2025-06-22 202525 3 23737 18844.0 28630.0 35 \n", + "2025-06-23/2025-06-29 202526 3 27513 21273.0 33753.0 41 \n", + "\n", + " inc100_low inc100_up geo_insee geo_name \n", + "period \n", + "1984-10-29/1984-11-04 37.0 213.0 FR France \n", + "1984-11-05/1984-11-11 184.0 308.0 FR France \n", + "1984-11-12/1984-11-18 123.0 195.0 FR France \n", + "1984-11-19/1984-11-25 99.0 163.0 FR France \n", + "1984-11-26/1984-12-02 110.0 176.0 FR France \n", + "1984-12-03/1984-12-09 149.0 219.0 FR France \n", + "1984-12-10/1984-12-16 184.0 266.0 FR France \n", + "1984-12-17/1984-12-23 146.0 224.0 FR France \n", + "1984-12-24/1984-12-30 110.0 198.0 FR France \n", + "1984-12-31/1985-01-06 120.0 190.0 FR France \n", + "1985-01-07/1985-01-13 147.0 207.0 FR France \n", + "1985-01-14/1985-01-20 317.0 459.0 FR France \n", + "1985-01-21/1985-01-27 708.0 832.0 FR France \n", + "1985-01-28/1985-02-03 1074.0 1236.0 FR France \n", + "1985-02-04/1985-02-10 939.0 1113.0 FR France \n", + "1985-02-11/1985-02-17 784.0 926.0 FR France \n", + "1985-02-18/1985-02-24 652.0 762.0 FR France \n", + "1985-02-25/1985-03-03 618.0 722.0 FR France \n", + "1985-03-04/1985-03-10 591.0 689.0 FR France \n", + "1985-03-11/1985-03-17 458.0 544.0 FR France \n", + "1985-03-18/1985-03-24 405.0 485.0 FR France \n", + "1985-03-25/1985-03-31 319.0 395.0 FR France \n", + "1985-04-01/1985-04-07 207.0 281.0 FR France \n", + "1985-04-08/1985-04-14 83.0 149.0 FR France \n", + "1985-04-15/1985-04-21 66.0 116.0 FR France \n", + "1985-04-22/1985-04-28 44.0 80.0 FR France \n", + "1985-04-29/1985-05-05 55.0 93.0 FR France \n", + "1985-05-06/1985-05-12 59.0 97.0 FR France \n", + "1985-05-13/1985-05-19 38.0 64.0 FR France \n", + "1985-05-20/1985-05-26 35.0 59.0 FR France \n", + "... ... ... ... ... \n", + "2024-12-02/2024-12-08 149.0 177.0 FR France \n", + "2024-12-09/2024-12-15 190.0 220.0 FR France \n", + "2024-12-16/2024-12-22 281.0 323.0 FR France \n", + "2024-12-23/2024-12-29 278.0 326.0 FR France \n", + "2024-12-30/2025-01-05 320.0 370.0 FR France \n", + "2025-01-06/2025-01-12 363.0 405.0 FR France \n", + "2025-01-13/2025-01-19 356.0 398.0 FR France \n", + "2025-01-20/2025-01-26 496.0 548.0 FR France \n", + "2025-01-27/2025-02-02 475.0 523.0 FR France \n", + "2025-02-03/2025-02-09 385.0 431.0 FR France \n", + "2025-02-10/2025-02-16 293.0 331.0 FR France \n", + "2025-02-17/2025-02-23 186.0 220.0 FR France \n", + "2025-02-24/2025-03-02 112.0 140.0 FR France \n", + "2025-03-03/2025-03-09 79.0 101.0 FR France \n", + "2025-03-10/2025-03-16 78.0 100.0 FR France \n", + "2025-03-17/2025-03-23 68.0 88.0 FR France \n", + "2025-03-24/2025-03-30 50.0 68.0 FR France \n", + "2025-03-31/2025-04-06 47.0 65.0 FR France \n", + "2025-04-07/2025-04-13 44.0 62.0 FR France \n", + "2025-04-14/2025-04-20 34.0 52.0 FR France \n", + "2025-04-21/2025-04-27 26.0 40.0 FR France \n", + "2025-04-28/2025-05-04 21.0 33.0 FR France \n", + "2025-05-05/2025-05-11 18.0 30.0 FR France \n", + "2025-05-12/2025-05-18 23.0 37.0 FR France \n", + "2025-05-19/2025-05-25 27.0 43.0 FR France \n", + "2025-05-26/2025-06-01 21.0 35.0 FR France \n", + "2025-06-02/2025-06-08 28.0 44.0 FR France \n", + "2025-06-09/2025-06-15 28.0 42.0 FR France \n", + "2025-06-16/2025-06-22 28.0 42.0 FR France \n", + "2025-06-23/2025-06-29 32.0 50.0 FR France \n", + "\n", + "[2121 rows x 10 columns]" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "sorted_data" + ] + }, { "cell_type": "markdown", "metadata": {}, @@ -179,9 +3164,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,10 +3192,34 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXncHEWd/z/fmXnuJ8+T+05IgHAEBIFsAFFAEYiii6vg4ipE5beoy3pfsOuKK7qCuqB4o4CgKKcrCAKG+xASEhIICYSE3PeTPE+ePPczR/3+6Kqe6u7qY2Z6nplJvu/XCzJPdVd1dXd1fet7VBUJIcAwDMMwpZKodAUYhmGYAwMWKAzDMEwssEBhGIZhYoEFCsMwDBMLLFAYhmGYWGCBwjAMw8QCCxSGYRgmFligMAzDMLHAAoVhGIaJhVSlKzCSjB8/XsyaNavS1WAYhqkpli1btkcIMSHsvINKoMyaNQtLly6tdDUYhmFqCiLaFOU8NnkxDMMwscAChWEYhokFFigMwzBMLLBAYRiGYWKBBQrDMAwTC6EChYhuJqLdRPSqljaWiBYR0Vr57xjt2JVEtI6I1hDRuVr6SUS0Uh67gYhIpjcQ0Z0yfTERzdLyLJTXWEtEC7X02fLctTJvfemPgmEYhimFKBrKbwEscKVdAeAxIcQcAI/Jv0FEcwFcBOAYmefnRJSUeX4B4DIAc+R/qsxLAXQJIQ4HcD2Aa2VZYwFcBeBkAPMBXKUJrmsBXC+v3yXLYBiGYSpIqEARQjwNoNOVfD6AW+XvWwF8QEu/QwgxJITYAGAdgPlENAVAmxDieWHtOXybK48q6x4AZ0nt5VwAi4QQnUKILgCLACyQx94lz3Vf/4DjhfV7sW53T6WrwTAME0qxExsnCSF2AIAQYgcRTZTp0wC8oJ23Vaal5W93usqzRZaVIaJuAOP0dFeecQD2CSEyhrIOOC660XqcG685r8I1YRiGCSZupzwZ0kRAejF5gsryVojoMiJaSkRLOzo6/E5jGIZhSqRYgbJLmrEg/90t07cCmKGdNx3Adpk+3ZDuyENEKQDtsExsfmXtATBanusuy4MQ4kYhxDwhxLwJE0KXomEYhmGKpFiBcj8AFXW1EMB9WvpFMnJrNizn+xJpHusholOkD+QSVx5V1gUAHpd+lkcAnENEY6Qz/hwAj8hjT8hz3ddnGIZhKkSoD4WI/gjgTADjiWgrrMirawDcRUSXAtgM4EIAEEKsIqK7AKwGkAFwuRAiK4v6DKyIsSYAD8n/AOAmAL8jonWwNJOLZFmdRHQ1gBfled8WQqjggK8DuIOIvgNguSyDYRiGqSBkDfgPDubNmydqbbXhWVc8CICd8gzDVA4iWiaEmBd2Hs+UZxiGYWKBBQrDMAwTCyxQGIZhmFhggVLFHEz+LYZhah8WKAzDMEwssECpYlhBYRimlmCBUsWwPGEYppZggcIwDMPEAguUKoad8gzD1BIsUBiGYZhYYIFSxbB+wjBMLcEChWEYhokFFihVDLtQGIapJVigMAzDMLHAAqWKEVXqRfneQ6/Zy+ozDMMoWKAwBfOrp9ZXugoMw1QhLFCqGPahMAxTS7BAYRiGYWKBBQrDMAwTCyxQGIZhmFhggVLFsA+FYZhaggUKwzAMEwssUKqYap2HwjAMY4IFCsMwDBMLLFCqGPahMAxTS7BAYRiGYWKBBUoVwwoKwzC1BAsUhmEYJhZYoFQxvKc8wzC1BAuUGmZ9Ry9mX/kg1nf0VuT6LPAYhtEpSaAQ0ReJaBURvUpEfySiRiIaS0SLiGit/HeMdv6VRLSOiNYQ0bla+klEtFIeu4GISKY3ENGdMn0xEc3S8iyU11hLRAtLuY9qJay7/vOK7RACuG/F9hGpjxuWJwzD6BQtUIhoGoDPAZgnhDgWQBLARQCuAPCYEGIOgMfk3yCiufL4MQAWAPg5ESVlcb8AcBmAOfK/BTL9UgBdQojDAVwP4FpZ1lgAVwE4GcB8AFfpgutAIazDppGphi8sTxiG0SnV5JUC0EREKQDNALYDOB/ArfL4rQA+IH+fD+AOIcSQEGIDgHUA5hPRFABtQojnhWVDuc2VR5V1D4CzpPZyLoBFQohOIUQXgEXICyFmhGCTF8MwOkULFCHENgA/BLAZwA4A3UKIvwGYJITYIc/ZAWCizDINwBatiK0ybZr87U535BFCZAB0AxgXUNaBRZX311VePYZhRphSTF5jYGkQswFMBdBCRB8LymJIEwHpxeZx1/MyIlpKREs7OjoCqscUCisoDMPolGLyejeADUKIDiFEGsCfALwNwC5pxoL8d7c8fyuAGVr+6bBMZFvlb3e6I480q7UD6Awoy4MQ4kYhxDwhxLwJEyYUeauVIerikNyvMwxTDZQiUDYDOIWImqVf4ywArwG4H4CKuloI4D75+34AF8nIrdmwnO9LpFmsh4hOkeVc4sqjyroAwOPSz/IIgHOIaIzUlM6RaQcVVGGvPK+GzDCMTqrYjEKIxUR0D4CXAGQALAdwI4BWAHcR0aWwhM6F8vxVRHQXgNXy/MuFEFlZ3GcA/BZAE4CH5H8AcBOA3xHROliayUWyrE4iuhrAi/K8bwshOou9l2olskmpQrYnNnkxDKNTtEABACHEVbDCd3WGYGkrpvO/C+C7hvSlAI41pA9CCiTDsZsB3FxglQ8oqOKBwwzDMHl4pnwVU+0KAGsoDMPosEA5AKhUv84+FIZhdFigVDFhEwcr7ZRnGIbRYYHCFA2bvBiG0WGBUsVUe39dTP027unDR258AX1DmdjrwzBMZWGBwhRNMWt5Xfvw63h+/V48uYZXLWCYAw0WKFVM1P66UqanUi7L/h+GOfBggVLDqD65UtFW7ENhGEaHBUoVEyYoanGUz0KIYQ5cWKAwxVOCcKhBWcgwTAgsUKqZKh/N88RGhmF0WKAcAFTMKV/EdVkIMcyBCwuUKias66UKO1E4yothGB0WKEzR8J7yDMPosECpYri/ZhimlmCBUsVU+xbAxVyXhSTDHLiwQDkAqCWnfB52ojDMgQYLlComrMOutGObI7YYhtFhgcIUT1FhwwzDHKiwQKliqr3z5bBhhmF0WKAcALDpiWGYaoAFShUTugVwhR3bRc2UZ9nHMAcsLFCYomHNiGEYHRYoVUy1j+ZLqR+7UBjmwIMFClM0VS7vGIYZYVig1DCU37KxIhS3lldwnq6+YfQPZ4qrEMMwFYUFSgE8/vouvLxlX6WrYVNheVISfisln3D1Irznx8+McG0YhomDVKUrUEt88rdLAQAbrzlvRK53IPtQgti0t788BTMMU1ZYQ2FGlGoXkgzDFA8LlCqm2sNyOcqLYRgdFigHAJXa6KraBR7DMCNLSQKFiEYT0T1E9DoRvUZEpxLRWCJaRERr5b9jtPOvJKJ1RLSGiM7V0k8iopXy2A0kPbZE1EBEd8r0xUQ0S8uzUF5jLREtLOU+qpWqX22YF4dkGEajVA3lxwAeFkIcBeB4AK8BuALAY0KIOQAek3+DiOYCuAjAMQAWAPg5ESVlOb8AcBmAOfK/BTL9UgBdQojDAVwP4FpZ1lgAVwE4GcB8AFfpgosZGYrbYMvKVWlhyDBM/BQtUIioDcDpAG4CACHEsBBiH4DzAdwqT7sVwAfk7/MB3CGEGBJCbACwDsB8IpoCoE0I8bywepvbXHlUWfcAOEtqL+cCWCSE6BRCdAFYhLwQOmA4kEfzLFAY5sCjFA3lUAAdAG4houVE9BsiagEwSQixAwDkvxPl+dMAbNHyb5Vp0+Rvd7ojjxAiA6AbwLiAsg5KKrdj44Es8hiGKZRSBEoKwIkAfiGEOAFAH6R5ywfTmFQEpBebx3lRosuIaCkRLe3o6AioXvURdbXhmtpTPvZaMAxTLZQiULYC2CqEWCz/vgeWgNklzViQ/+7Wzp+h5Z8OYLtMn25Id+QhohSAdgCdAWV5EELcKISYJ4SYN2HChCJus3qptNmotLBhtnkxzIFG0QJFCLETwBYiOlImnQVgNYD7Aaioq4UA7pO/7wdwkYzcmg3L+b5EmsV6iOgU6R+5xJVHlXUBgMeln+URAOcQ0RjpjD9Hph1QVP9ovvpryDDMyFHq0iufBXA7EdUDWA/gE7CE1F1EdCmAzQAuBAAhxCoiuguW0MkAuFwIkZXlfAbAbwE0AXhI/gdYDv/fEdE6WJrJRbKsTiK6GsCL8rxvCyE6S7wXpkB4gy2GYXRKEihCiBUA5hkOneVz/ncBfNeQvhTAsYb0QUiBZDh2M4CbC6lvrXFAd75s8WKYAw6eKV/VRJMoFYvyqsxlDzjectUj+MrdL1e6GgxTMixQmKI5oDWoEaRnKIN7lm0NP5FhqhwWKFVMtXfYxazlVeW3FDuv79yPNzt6K10NhhkReD+UA4BKLdLIqw2Hs+BH1mZhI7WHDsNUEtZQqpiw/tpv18ORoto1KMafdDaHx17bVelqMAcYLFAOAGqpY6/YUvtC4Bt/XonV2/dX5PrVxnWL3sClty7F39ftqXRVmAMIFihVTOjy9SNTDV9qaT+UXfuH8PsXNuPjtyypdFWqgjd3W36d/YPpCteEOZBggcIUTSnKxkiLolTSEr+ZXO0IwXKSzuYAAHXJ2ukC9vQOYdYVD+JPL3FEXLVSO63pIKSWNIBqoXcog86+YU96KmEJlHQmN9JVqkrSWatt1ZJA2dLZDwD47d83VrYijC+105qYqqMkd0iZZOXp338CJ169yHs5eb10jgUKAAzXoIbSXG8FpQ4MZ0POZCpF7bSmg5CwDrtS+osKLitFgyqX9mXSTqzrWWSyrPUBQE6a/pKJSnviolOfsrqrfhYoVQsLlBqmchFT5SmjnPeTk2WzD6X24Y3dqhcWKFVM1O+mlgRLUJ5y3kapZW/a24c3dvXEU5kqoNJ76RSDGhSwOKleWKDUMJUeqJXkQjFpKCWUF3q9Eks/4wdP4pzrn46pNvHx1Bsd2L5voOj8tTTar6GqHrSwQKlionaCFdsCuIQv3JQzV1YVpXxFV5KFNy/BeTc8U3C+Sm8fXQy1JPwOVligxIQQAj9+dC12dBc/Wiz4mhXuDorbU94/F8uT4ujqL2JyYoVMXh09Q3jLtx7Bq9u6C857IL/DAwUWKDGxZlcPrn/0DVx++0uxlRka5VWhL8yO8iplYqMhczkFZFm1nxpmpB/Ls+s60DOYwW+eWV9wXtuHwq+yamGBEhNqesNIhjSq72qkP7BYorzKVK7v9bgTclApn3x9MgkAGCpigim/w+qHBUoRmEbXcYzaC6/HyF3LpwbF5zSGDZdQlbDrla/omsakFW7bN4DuYsxoEWiQc0mKESisZVY/LFCKwDSVISElSiUafS3thxIYNlzG+2CHrhM7bNjwWE675nGc95PCHf1RaKizupyeIhal5FdY/bBAKYKsQaLkZ4/HR/hM+Rp0ytuZDD4UNnmNGORj9FKCd2tXeYJLWhqs5VO2FVG+eoeVbveMPyxQisCkhVDAsXJhf2AV+r5Kc8p708r57FigmHE/FtNgKdbr2UKhiLwsSKoeFihFYOr47O8wxjZfrR9Q3l9UzJ7yAWHDxVaoxOvWKqWY8dQ7dLfl8q9MU/wFeNWc6ocFShGYRnEqbSTbfKX9AnEvNnwwmrxKeYdxdLDuy5db8Ko6F+d/47DhaocFShGYVkBXI704zTaR1/KK7YoRr1eDi0NWug/yu7dShEI5NJRyd9al+EEq/Q6ZcFigFEHW8NUpDSXO+P6wD6gWfSiVWhyy0iGnfpePe/maqPgtvVJ+gVJKqPnIWwGYwmCBUgRGk5ds7DSCy7hW+sOKez+UsvpQKv2wfChNQ4mhAq4yyi14SzN5xVsXJn5YoBSBaZSVK4eGEvIFBYXgjghVFuUV/Lyq099UilAu5Xn5mryKLjEapd1vjBVhygILlCIIMnmN5JoWlYpcimPOTTmc8pUyp5VCpevlvn65NZSSwoYr/bCYUFigFEGgySvG60T9fGrKhxJ4rLQbCeoMK90V+XWGcc/lKbiMMpQZeL0STF6lmMuYkaFkgUJESSJaTkQPyL/HEtEiIlor/x2jnXslEa0jojVEdK6WfhIRrZTHbiDpiCCiBiK6U6YvJqJZWp6F8hpriWhhqfdRCKYor0o08ko55eOJ8jKFeZVWZpBJpCJL4kS4ZuVMXualgsqtBcTtd2Oqizg0lM8DeE37+woAjwkh5gB4TP4NIpoL4CIAxwBYAODnRJSUeX4B4DIAc+R/C2T6pQC6hBCHA7gewLWyrLEArgJwMoD5AK7SBVe5MX3ISmtJxOiUD196Rf1bmQ8t7uuWaiMP1FAqKPCBAB9KKRpK8Vl9r1/u51TxIASmrJQkUIhoOoDzAPxGSz4fwK3y960APqCl3yGEGBJCbACwDsB8IpoCoE0I8bywhke3ufKosu4BcJbUXs4FsEgI0SmE6AKwCHkhVHaMPhQ7ymukagH7C6slk1dQL1jexSHjKid6QVHOLEXLKGkeik8Z5fehlHK/MVaEKQulaig/AvA1ALoRaJIQYgcAyH8nyvRpALZo522VadPkb3e6I48QIgOgG8C4gLLKhv4h5AzDrHyU18h5UdTRu5dtRfdAeZYbj3L9ovKaLF7l1FBiElaF1FFvM77zUEqpSwl5/cooe5RXCZGJ+ffLkqVaKVqgENH7AOwWQiyLmsWQJgLSi83jvCjRZUS0lIiWdnR0RKqoCb1DCJzYWKGdi256dsOIXauUtbwUpg6+1NFxkDklrtFtIXWMcmZJJq/CtxSx8du/p+waCorXqislRoQQ+Onja7G3d6hCNagdStFQTgPwj0S0EcAdAN5FRL8HsEuasSD/3S3P3wpghpZ/OoDtMn26Id2Rh4hSANoBdAaU5UEIcaMQYp4QYt6ECROKu1M4G7MpyqscH2KoD6XSIadF5Smfn2MkfCiF+ACcPhS/KK/KOKn9VsceqSivYqjUagdLNnTih397A1+/d2VFrl9LFC1QhBBXCiGmCyFmwXK2Py6E+BiA+wGoqKuFAO6Tv+8HcJGM3JoNy/m+RJrFeojoFOkfucSVR5V1gbyGAPAIgHOIaIx0xp8j08pGmPkiK0eLIztTvkLO+FiivEovo5AyYzN5FVBOlHMPtrBhO/S3mMwVGkCls9aF+4czlalADVGOeSjXADibiNYCOFv+DSHEKgB3AVgN4GEAlwsh1Absn4Hl2F8H4E0AD8n0mwCMI6J1AL4EGTEmhOgEcDWAF+V/35ZpZSNMQ6nEPBT94y/0uoPpLO5dtrW0MNGYO8NyzpSPa5Z1YT6U8HylVCvseQ1nctjdM2g8pgY+teSUV3UbaUWlUmbsWiQVRyFCiCcBPCl/7wVwls953wXwXUP6UgDHGtIHAVzoU9bNAG4uts6lYPKh5CrsQymU6xa9gRufXo/RzXU46+hJRZVR2pwCQ1rJJq+A68XUC8Xd4ZYU5RVy/At3LsdfV+7Ehu+911dz9oQNF12baOQnJxZ+pZqKZjxI4ZnyEdEblSnKqxxOeecI13vNUtr5Hulg7OwbLrqMohyrgWap0qhmDSWO8grN+9eVOwGY7z2/pbxLQyn7gllC+39h2BpKjLUphFoZLFYSFigR0T+84Yw3vCZv8ipPq4s7zLYuYb36TBEdiF+EUCEYF9gsa5TXyGsokXwoIzBz3LhltVocsoRIsWKo9ETO4q7LKkpUWKBERP8Q+oeznuPlMHmFrjasNfRCr1uXsjKks8X3KKV8ZuUweY2EhlJ0lJefD2UEnPKFrHHGy9f7U67B4oEEC5Qi6E97BUo5nPI6xm+phA8sJTUUFcFSUF3KNsqsfg2lkCpGOXUkBIr5PL+1vIqvT6S6lKKR2U551hiqFRYoEdHb8KBBoNi25zjX8oqtJC91SauemVI0lJglSznX8opPQynA5KWHmvvNQ6mQyUsrJPq5MVCSUz7muhzI3P/ydpx7/dMjLnxjifI6GNA/3kCnfLmuL4Sn9FKaSiJhlWWKWItcpxKub9yxsVSTl6Ms4YhsqowPJUp5xdclalVMYe75PW1cGkrx1YlEHGHDTDhfunMFMjmBoUwOjXXJ8AwxwRpKRPS2bHJkK8tRuaK8zMeL/8DisAcXZwdXZgvDsVJNXtp7cb+i2KK8Cjk3kg+l/B1sUJSX+1g1m5M4bDg6TVKImPy95YQFSkT0NmX6kMuxBbDf9e20GBp6MWXkhWa8X1rpTnn9t7OwuEa3BZUT4dSR8aEEOOXLJHj9KCX0l8OGo9NYbwmUvqGRnd3PAqUIMgZHdn75+jh9KJoN3jiiL55EHAs8xuycLz1sWPdZhF+vKAqSJ/71sdNHoHc0aih+e8qX2ylfg6P9Gqwy6uQHPlyCj7QYWKBExLF8veGrKIsPpYCWXKgJqxS5F4cgiXteDeASKGVyNhcbNux7TgndVXSTV4CGUmSZxWI/vyIuU6mlV2oRe0fOsk9UdcICJSL6azH5UEqdh5LJ5uzZ6+brx+/ELrWMYrLGUee9vUP4xC1L0OWa5e9wyrudzRV2yldyT3njxEaYZ6eWX0Mp3mzFgiQ6clZAUROXS7ruiF6thtEbs3n5+tLKv/qB1Zj3nUfRM5jfKCusyFImNqqtikupdmnCyPQMoxX462c24Ik1HfjDks2u+gRoKDFp/oU55cPPLudaXvY1DPeeN3nFV58olFL6CPeNNY0aMJj6qnLCAiUqIQKllPBbAHhk1S4AQM+g2YkWt4kovwVs8WWUNkktWpqJrOwhkwmnFNVfS9lMXgV8oMLnt985hRJV6wreJ6bwGvzv39bg/J89V3A+/XrFLQ5ZGYlSi+HKCR8fWbnheSgR0TtPo4aSK82+qzpHveyRaAvFCAUiq25FhQ0XeUxH+RlTHoGiO8FdppyIZYdRUJBXFB9KCRWLKtsC1/IqQvD+5PF10S5sIA4TX5Bgufz2lzCutR7fPt+zeHnRjLQfIg6UD4VNXlVKmMlLaSh+H+SmvX3Y2tXvW74SKCM2orD3wyi+iJJG14a0qPfuq6Foph13UXGNbovdYKsc81CivoECJsqXfRBTygZbUZ79gyt34LbnNxVU7rNr9+DCX/7d1zyUXwSjduKG7XlGIyxQWEOJiP5agvZD8Xt/Z/zgSQDAxmvOMx5XnaM+ogj7gEqb2KiuUTylLb1SfJCBev5uDSUoTLcSi0NGmodSdE3iccqXMl8nmxMeoR5GKXvKl6tv/MKdK7Cndwh7+4YwcVSj53jWtj7UjqaiZB/7UGqAoKVXim10tobi0wDinodCEZwo63b34rUd+yPVpVB87jJSXvWskwln8622iY3C94/iygssPwBjk/LZgqCQ2vitVJ3LCfxxyWbj8VhMXsUXYSRsO4ZaEiQKpU2xQKlS9EZlXnpFCpQiy1ejbX31X0cHGXPYsD1CDTjn3dc9hff8+JniL2IioNJR276aWJp0tV69cy7b0itl9KEU2nFFNxH6j0ZKCa/2mzR370tbceWfVuJXT73pvaz9nRT+QsplDg4bW9WyyavUYKFCYYESEYfJK8ApX2yjtyciFTQCzp9baFOP49uIew5FoSYvt4bieC0BPhTj7pdC4NO/W4a/r9sTeO1COlynCS5c8yz0eUY9P2gzM+9aXtGvv3ZXjzG9e8AKfd9r2A00Fo9RzH1kPoTeXPBId8pxkGCnfHUT2Slf4nwHR5SXz/XjoJSwYb+VagvBr1OPwuL1nQAKjPIS+nneMgfSWTy8aic+eeuLgdeOe6Z80HIxcZRvXcM/rze8Ovr1n1u315hOAQEfJUVDlktD8Yl4y1/WOlA7+om+IycLlKrEETZsaNhK+y9aQ0Hh+UsyecUQp16JsOF1u3uwbd8AgPwS/Kb6BM1DMQmuvJM6+PqFRXn518dEoe+ilKVX8gstFu5rmjG2CQAwtqXeeDyo4y2leyulb9za1Y+7lm4xHguLiBppP0QcVMqHwlFeUdE1FMPikKXOQ7HL8en4TMWWcimq+Ex5L1E6M2VOAYAkuQWK//PSvyvTNxbVrFGIBhpF43JqTuX5+M0CxXt999/uPWUUarfPYupbmpm0+Mz//KsXsG3fAN5/3FQ01Tv3BwnrfMuwd17ZsX0orKFUJw4fiklDCZmHEkY+zK+AOmmXGsnGHkuUV5ESMqHdqPuec67O0HlMGH+708KqUKwG6ZcraEFLAOjuT+OvK3eElh+1Hvk04fjXnR5UflhIaj5qqvjnbKKUOSxKq9UHJApVXz9/Qy0uSrlaRmeyU75KCfOhKA2l2BdoMkE5TSbBPXChlw0LlYxC3GbwKOWlNEe8u3MPivLSCw+27Y/sBxhmFvvy3Svwb7e/hI17+gx5o9XVvPacuZMUhnPchI1+R3JPoELpHfIXKL4aSg2avBSsoVQpYUuv5J3ypb3AYvMXms1vYlshhOX90aNv4NLfmp3cxZq89MCuQD9JgG/A1BFH/fCK1lB88jnD0b3qaUePtQJ1V783YqqUpVeEfcz/3LDBUdgzC3zHRTS7sPYW5R2atJBEFZq8+ocz+PAvn8cbPpF0UWGBUgMYBUrITPkw7NVBfewkRv3EcXzkGk4+yiuYHz26Fo+9vtt4zBzlFX7tIA0laDXGMB+Kn0/B77woRHknenmmjdua6y03p2kr1+iLQ/qnBUXD+RVvd8B+GkxAlJfwuW4Uwm7Xb6KljunbDdO4KhE2vHh9J5Zs7MR3H3ytqPx1SZ7YWNWEmrzkCaW+QIcfIMSkFeXj98NeAqOwbO5CYsVhbgmxzwPBIa/u3HH5UPw68b+t2okX1u91netfH9ORtEFDSSXVhFfDrPPAmur1MAlv8wAo7DkBmonIIABD61Jwjjx5M525lCgCxbyUv5qzYc5fiZny6vsscGUbG9UPcZRXlaK/FnPYsPpAQ8wAPpEzCr0jLSSEtdBmY48USzF5FTPKDMjjNrckDNb4oA4vyMHtUF4M/UbUpXP8vs/LfrcMgHOttihPR7+cSUNJRBjth2GMavNRyRxC3ad8ozYdkVIc3GFZTM/Pc45BaISF7Nub54WWHh+qmsXOzlfvjp3yVYre0ZhNXv7HdHw/UoNjMGjmt1Unvz+iEyXbiO02GKIFAs4RZiGj67AtnIM0lLBZ9n5EiZjS78E0wjalsu6AAAAgAElEQVRta1BoXUz3q3wJ7iNh7dzKIwdPYVFeMS8XFDZYM2l4Ucqwo7x8BJJKHsmuWdWzWA3FLodNXtVJVJNXWKP3U6tNo6RCOq9C203Op0MxkXZ9aGrUXEpbNXY2WprvaDFQaOhlufL5rECQP+4tI38d/3KDiKahaCYvo4ai6mB6XtEwdSp+SwVFCSRQ7d8vzDZoFYZybgEcRUMxWcXCnPKVMHnFtX4YL71SA5g+7qhO+bBBlG/YsLEDDj4OADu6B/CQYR6DOjvKt+IWgsXM6vdc39Rxa5fx+xCCOjynz8l5zOmUN7w/mWb6fh3aTQEfaCTtT/tt0lAoQHiXZPLyMT2FBS/o6cWYU9R7LWbkHSaMijZ5RZyHMrLzGsOvuWxTJ9bsDI4CqxmnPBHNIKIniOg1IlpFRJ+X6WOJaBERrZX/jtHyXElE64hoDRGdq6WfREQr5bEbSH5FRNRARHfK9MVENEvLs1BeYy0RLSz2PqISZuuO6pT3/QjtUZL5mmF18jv3w796Hp+5/SVPR1hItI175Bw0cS2MIPu7nuLXcQeuKBzw7KI65c3XNP8ORx8chN+P2YfiX7+Sll7xGQDp7dc36knlDbHfmo6W0sGFRnlFMXmZnPIhPqFKTENRVUkEaCgf+sXzOPdHT3vSo7zDclGKhpIB8GUhxNEATgFwORHNBXAFgMeEEHMAPCb/hjx2EYBjACwA8HMiUmsg/ALAZQDmyP8WyPRLAXQJIQ4HcD2Aa2VZYwFcBeBkAPMBXKULrnIQtpaXbfIKEyghZgLfpVdMI1SHich8vS2d1gxhd50LcY5msm4NpXiTV8IWRt5jUez3xTrlndcxlBtwM1FWDY56HUPhNqYOMWjV2Ogair9GFrR8fbEmr6D6qcGJX9Y7lmzGrCsexIAhTDostNs0+dON6dsNi1oL65SXberEI6t2hl67EEqZ+6JrujXjlBdC7BBCvCR/9wB4DcA0AOcDuFWediuAD8jf5wO4QwgxJITYAGAdgPlENAVAmxDieWG14NtceVRZ9wA4S2ov5wJYJIToFEJ0AViEvBAqC3qbCpqHEvYCwwSOv8krmLCOzl3nQpqZp+OIYWFJo11ff8YRRotep7z+uzBhE/TeomiCxny+f5iT0xmDQFE7eRo3q4pWGeP9quLcGkqABpjPGzx4CqqX2r7Zr+387Elrv/rdPYPeckNa7VfufjnwOGCuc9j+6/Zqwz69+4d+8Tw+JaP84iIfNly4RNmvLS9Tk055aYo6AcBiAJOEEDsAS+gAmChPmwZAX+5zq0ybJn+70x15hBAZAN0AxgWUZarbZUS0lIiWdnR0FHeDCB895yJGeflPBpPl+HR8xo+0gI7OOwlQmNMNuG37eW0qNGtAfUypWmfmY70IGkE73lGBPpQga4kIyRslnx8Ok1fALG6zmTVaPcw+P3PHHvacgPyz9WvL+W/Ae1xpKH7Ppk5OXDUFKISZaYcNAtmN6RmHTWzMa0Yj1znbVSlCQ9nS1W//rjmnPBG1ArgXwBeEEN79YrVTDWkiIL3YPM5EIW4UQswTQsybMGFCQPWCierULdaUkz+u5Qsxszk0mJDG7r6u7ViN0ODcnVl+bkTUEbJX6zJGLUXQUBxh1T735D7Pfb1C/RFhE0wj5fM7RztgcsorE6GpYyhlx8b8vBt3fcLbqr1uXYiPxawZ+Q8IgPxETpPzXAR8YwBw7LT2wPr61VnNKh/Oes1sQGmaeLGoey1GQ9GbUU1pKERUB0uY3C6E+JNM3iXNWJD/qrU3tgKYoWWfDmC7TJ9uSHfkIaIUgHYAnQFllRFtJGkYPQWZvKI4yUx7MjhGxobBV5R5Dn757bkEPvmC1pgqdC8V59waf81IT/F3ypvPt/72f0dhZitVR9PnG2XUbiLKqfop5rBh/w428pL7htNUWpBT3q94W0OJ4Odyo9+HKbvahTOd8W8fYQJFCQh3ff3qVp+yrjkwbNZwCgmxL4aHX92B7zyw2nismMgy/b3UjIYifRk3AXhNCHGdduh+AAvl74UA7tPSL5KRW7NhOd+XSLNYDxGdIsu8xJVHlXUBgMeln+URAOcQ0RjpjD9HppWNUJNJwEeWjTDqM20BHDaqFj6/TfjNN4iyqZAnyssuM+SiqizH/Zvr406LMicgyIfi8RmFCN9ADaUAwe3MZ/7td12zhqKWXjFpdFHrHKChuFpOJJNXVrV187VVepBT3q/8vEbmLTxsrlffUEYed9cnuINVAmUw7aehOP+Nm0///iX85tkNjjRV51LCq4GRd8qXsvTKaQAuBrCSiFbItP8AcA2Au4joUgCbAVwIAEKIVUR0F4DVsCLELhdCqDf4GQC/BdAE4CH5H2AJrN8R0TpYmslFsqxOIroagFrK9ttCiM4S7iWUsM7ObYrRdxN0zu4OfsFRTD1h9YtSbthcAv3Dc4/2TcIv8NoGM4fpdqKYvBwaSgE+lGzIO4gaNlzYjo2FqSimd6yaUdbUwQa4DHxXXHAdd992lNWG8xqKz4hemaYM96932sGajEmAyn998vTJyDB3uUFmUgBoSFnBpr1SILkJW0OsHKj79wsECDJlOTSyEdZQihYoQohn4a+RneWT57sAvmtIXwrgWEP6IKRAMhy7GcDNUetbKuod1SUpVKC416GKoqEofE1eIT6HcJOXq/OFv0blTneHswbNhA4rK2jNMz0l2jwUt5DUr+Osc5jZI3BdQf09GM7L+QwmHO/HpxsMWx0gSEOJ2iEHLjXjOrazOx9dFRY27DePMGiC4WA6/wBNxQctg5LXFHwGQfIlCmHV3RS9ZWrvKfm+Nu41hx2HXbccBJlgAWA4oMHWpMnrYEM1plQiETgPxf0bMPsQ3JjMSGF7sBTSVDx1ln/6xc8HaSh22HDExuoWtkCE4AVfDcV/5O2Yae+qs6PTNxQdrKGECCNR/Accpr3mw4YL1aqC6+y3ltd1i97Q8gWX7WsuDXjHQ5moGoppdeV8uSZhl/EZjJnan4467G/yCh585a9ZWuctDO3IT0MJWlnZ8e5ZoFQn6h2lImgo7m8hV8CIwa8cs4koXFCZytXP7x/Oon/Yq+o76uxqvIWu5WVy9JoXO/Svr+kcj19Iv6ZbqDv+NnQqATcTnNM/6CKKBumoc6AJy+RTCDo/REMJ0BRNZZjS/dalC9q9dCgd7JRXhE3kNPrBfO5Zr2fQOnymPWf0csPau0mLLAQ9v/rm/HwoQSHSejuqmYmNBxt5k1fCp1Hmf/vNSgeCnPLec/Uzw0xEx08fba64yu+O8grpvPUPOu06XkqUV7DJy18o2+kBDvJsQMcRpqHk1/LyfsFhDm7ns9JMOhF0yNBwZlm2+x34na8IGuAAwVpEUPlCCC3k3OfaASP6oawuUAI0lJAghCCtyzqeT9fvP0igmGbn62WFaSBR9mOJmj/vlPfTUPzrEqW/KRcsUCKSN3lF8KG4jjudZOby88uZmDsvv4ahljdPhLxJt5BzjoyD78etoYiIH5jp2kFLvuQijKycHYXznGF9hOd+B6E+BePlvNc0CSPtuvpvpw/FjAi4H+t6UhswdFZBAQph95uP8vLH9AqimG+D5qnoI2vTvjTqOwi9X9N1fQSOM1TZ/92HmbzCNZTSBIouRJXA8JuG0jOYNh+Asx2xQKlS1GupSyYCP1DA2zFE6Sjz5XivCfh8vCJ6nLo3jDb/2zTa0T9C9/GgSC0T+ocSZGoJE3LufO5T9KVL3OsyhQl1v87Eqpe/VgS4npVDQwknLKpK3aNxYqOPeQ0In09iz62I6Icx1dHPfBu0W2A6qoYSZlY2aSg+Iclh88DUsxgIFSjBbzTIUR4Fp2nO+u3nQzn7eu+ikHbeEnx6pcICJSK2hpKk0FnLHqd8FBXUZEbyUdvzh0Xgjn5+9XP/HaqhuC6ujkQ1eZmEQNiGUX5lB50TtCiecz8Ub9lBazEF+W0A/w7LT9v0yxsULBDmlHcfDY1qsx3rzvSjJo8y1s0uK8IsbFW26TvRNZRCnfJhQtLvWYZFeal6vLGr1/ieomooUZZ+CcIRCGM75b3nFbIqBjvlq5DNe/vxxOvWhP9IJq+Aziw0ysvnXOPSKwKRVZQg1dekqjtNXq68QtUv2rUzhmdjDjLQ8vjYiIPMPPoHHWR2LNRPGSZQ/CLiopgbwiJybIES4pQPNnl5r2ta2eG25zfidW1/DVO+MAc3oGmhoRqKMbt1HaPWHE1Iuo+HmYD0tD29w57j6nBY5xxlP5bA/IZ2ZLrmkEtweZcgygujEqtUMLynfATeff1TdmeVSiR8P/wEWY3P/e0Hds4SlZr1GXX6RUUFLQfvVwfr/OgfmZ9tOPIWtIYRkylv4Ts2Oo/p9SzUh6JIGsJqComg8xMufkVE1RTDNBR31rBQWXVcf2a/eWaD4xzjaD3Ewa1fz6ShpLM5EPmH/gZtdqUHXZhNj9rz9tmszXdgJunqH8aEUQ2O41Ei4oA4fCje9mt6xm7zbDqXQ0Mimf9btpXGVNJ38mm5YA0lAvrIN5Uk3w+0Lmk9zqCQ1TBTjp9d3K/zpogqitfklf9t+nhN6rddF58y/TCNLMMmavr7UPR6uDQUh3bgmtgYMmJXJE1RXo7rh2ko4SN4Z72Cy1ZJ5jBaffDhL0BNKwqow3rbbq5POs4LXUbI14znn38ok0ODXOokUEMxCpT8b7Pp0WxOcwjXEMFsmi0f3SlfmjpgChs2PWN9cqh1rvMcpcE01SfZKV/tpBJkDsPUBIrXKR88CgXyH5Bfh2me2CjyuyeG1NujoegLKYbYq71RXtE+MFNZWcN95uuk5YngQ3EPCJ0mL3cdzGW4McX9O7cA9h7XOzKnIA4fHYaZQ9Xz2rXfuz+I30Q+PZ/7GtZ18r/1UXWTS6AEmcrcvx35AkbX6WzOXuokOGzY1CaDzWV+AidMQ8kKgfqk/3peKnu5w4ZNmq7JGuIOHnALX/UdNLNAqX6SCfKO9nMC+wcz9iqn7oYVZVRnmgymtwU/VT3qUvKeSYC6hhKwejJgiPLyKdMPkzM1yF/gd9xzTpBT3tWZB5mHdBIGiVJOH4pp4Uwddb2OniFv3iCB4tOhuvPpWp1bQwlyULvLMV3bL2xYaShBTydsbpR5wUuzwNGFk9/Gbs0N1r2bBYpw/OtHqVFeevtN2xqK9zz3fBm38FWrEbBAqQFSSe/SK3cvs/b66uq3YsPdH3CUyJi8hqKrKP4dhnVu9NVI3W3dIaxCTV5+81AiXtsgUIs2eQUsrxLVhxIkfMN8KMZnpc8f8Am39uuMos6P6TOYYhw+mgCTl8eJq11nWFsKpaku5TrPW9+wiCkg38bd7Sabs0xtDXVmTV4v02Q+ihpoAJiXMQH8BmYCzXVKoJgW4VRt1ntNnQt/+XzwCSGY3plJUxvMBGsoajWCprrkiDvlWaAUSCpBHofirv3O0aNbQ/FTv3VsU5D+URiOu8tVHWBYu/Hm9xcY7vO9Gkq0EZuprCB7tJ4UxSnv1kKGMzk0yY4hcHJpQLVNPhTHswjoBIN+3/LcRuP1okZ5mWz7QXuLBAmUjM+79fhQjE75cCFpO5Rdh5UpplGavPwGSVb9w8KGowscv6AJ/XhzgyVMg0xepvcT5wrE+vtUz8pU30GpoXzoxOnGc4azOdQlyRr8slO+ulEduEOldr1Q98g5Slx4Pqojn+YYsfuYH0wjahNex6x/fa365AKPq+v74VwLTL9//7xB/pH8OXodvR+S8gMELb0S1Akoc4xOkD/JXQ+nlpT/vXqHeTPTqNFYJnNKkAlIPzbk6iT9IvjcAiVszk34xEaXwJfXsjWUgDZg0lDCQo79worTEUxeLfLeTZMbg0xecZqU9O9kOIKGMqoxZTxnKG35qZJknuJQTligFIha6tphh/aMAM0fEhAUGeNttCJk9JrN5ZfoDlNRgtYXC5uVXIzJK+0XcRPR5LXfZ2mJIPNTOhugoQSYS/TnnEwaNBT9WYWYYvTFD7MBnbqxXob3oDpD06Q5xzIvrmN6PYLmLejlKmF6+hHWVtlGLSCCn0u1F/ezcmsopjZg+w5MI3OfZ2uqj/9MeW99czmB5nqloXhPsP1+IUJscluj94QCcAwCpNB4Yk0Hnlu3x3Ge2lmyVWpVbuE7nM2iIZWw/L0jq6CwQCmUtqY6AM4G64kLd9v2AybcKdRH6Oc8NpqIhNlEYyJo6ZUwgeLvlPe/nt6ZGMOGTTP/tUpt7eo3lquKSiW8KxYMZ3JolKNfjx/LR1ADToHvF8GncO8NY10rn7ZvYFhL1zsIv7k8eh29x1UHmxPekahufnOvi6W3ySGXzV0XCu7BzpjmOnz1nCOtvwOeRZ1P+DygTcrzCZrIayjevKo+pues+w5MAjaTy0drOXxtDoFimiAq0BLglDcN9kxll770irmsPy7Z7DhP1bFVaiiqfkOZLBb86Gn8/oXNqJcCJUqkYZywQCmQ0S6BctfSLfjt3zc6zvGMzCKsX6TyRLHH62XlfSjBKkrQiN30kekdoEfttjUU/2s6o51y9vn50Z7Z9KDoGTTvnqc6MdOqz+ms8DV5BWkoYcuBBAUCuMvu6ksb0/0eVZjDWBfm7g7LMdHP9f71jtE96s45BgvOwU5dMmEvNBpk4qlPJnxNoUrQuYW6ulbeh+LNr96Fab6IHt3kFpKqvkrLSvuYbP1Mxw2pJIiCo7xM71DVk8g5cIyKI3hAexe6hul+f8osp0xe6rnu60/bKx3UpxJIJIid8tXO5PYmAPkP684Xt9jHvv+h4wB4R1dB8yMU+2SEmP4N+pmN8mWJ0FWG7bJcLcvPMavQP17vxEb/EZtdpiHaKVTjkv+2Nab8BYrqlOoSnnpt6eq3J3p6nfLe+iuGQjRIpz8p2Jexr9+sofiP5q3yrPlN/iYgwDsq39mdDwZxZ1WdTn0yEayhZJwdbypBgaHotkBJmRdJ1e/JXV/b5CXNkmY/iVkYAdZ7agqIxsrkBBrlgMKpIQd/fzlhhYs31SWDnfImgS/Lbq5LFqWh+PmkHJqjqz+xNRRp8lL59Ofd2pCSy0SNrIbCS68UwL++YzbqXbN8dZPTIeOaAXhHsd+6f5X92/SCe4cydgfgmAQZEnaqm7zCgk0C58YYPl59syF33iCbssK5FLfXLm5cekWmtTXV+S7PrT6ehlTS07lv2tvvOU8RtB+K+hD9tnd2OHhN5kHtXlXouJUe7AwG8s/GbxXrIIFy70tb7d/u56naU3tzne/M6oZUwtGpZ3MCySQFbqCWC9AQ3eW7r6s6ybYmq9sxaRnqHk2TBAfTWYxursNAd9aYN5sVGNWY8OQP81PlhECSLEEX5JQPChlvqk9ib98w9K2Ho+A3h0l/10+s6XCUO5g2O+V1IdTakEKCyHcAWy5YQymAORNH2fM+VMPUNYSWBmXTdL7F7do+3aYX/CNt29WsT+fltwxGImLYsEegyNEoYB4Nql0cm+qS3u10Axzrpuvl/QDBI3aV1NZYh/1+GooSKAYNRZE0jMwsc453zxkgr6G0N9X72NDzv8M0lN09g8Z05Xtzk5Z2f6vOhuPas/fzwwDe9680zNFNdZ7OV91/U33SYabJ5ARSiYTdxoPWFqtP+QsUNWp3P0vVSY5qrJPHvfczbBh82PeUzqJdPkfTs8jkhG1Oc0TbBeyTA1jPI0GExlTCWCfdj+Um37nXQYjCo772D+QHILom5b6/VdutKMGHVu7AD/9m9RduTU+/51GNdb5abzlhgRKBd8wZDwC4cN5022ehOkQ9bFc5hIPW9DE5yTZ3aiNrx4cQYtsXIrJT3huJlrc3mzpJ1SGNakw5nKG5XH4dqMAoL12gZLwCKChiZkxLna/JK6+hODs0dQ8X/cMMY+ecFVZnadXbefH1Hb3WdZvr0Dec9XyEQQEK+vEp7Y1Yu6vXmO+YqW3G+0lnrDkD1sKiZg1FCf4gk4qfkBzTXO+wx+v1akwlHTsoZnOWTy6RMAtePW+9YYKvQnXgA+ms41mrZ9cmR9buCXpCCE1D8ZatNBTr/ryCP53N2fNJ0gYh4memy+WsfUca680aiqqTSatWz7mt0RxxFUQ6m8P8/3ksX08fDQUAtu0bAAD8QXPQj2qwnoWao6TnGdWYYqd8tXLNh47DU189E0Rkd0rq5esfq33MYHP/+NtmAQC+eV/e/AVYjepvq3cBAOZMbHV0Gvpoqm/IbNu1nfKGduz8mL0O3QY5wjFqKOks6pKEMc312KIJvCgLXbrLNI06g5bOGN1cj94hH5OXba5JGqOoDp3QYrQd53J5jcx96UtvXeqon7tTCQqh1tNmjGl2TEBU9Rvf2uCrXaSzOaSkhuInUFoanM5XRVtjyo5qcqsoA8NZJBOEloakp+N2aCjZnP0uMjlLeNUl/AdGUTQUNWrP5oQzqMClobhDqZ2+LLeZVWAwnctrKIZAg0xO2PNJnIt0Wr8bfMx0liCVAtYkUAwatvte1T0V4kfpdQ2a3BrKv5w80/67q8/yzan1AgFgbGs9AKCzz/Kl6e2jLmkNDPp9tjUuFyxQIjBtdBMOGdcCAHZ4YZ8yCckG/PG3zULKXssr3/BU56TyudH9K/WphONDyWhmqe4BbwebywnbpxMWqTXs6hyyAmi08xrMC8NZNNUlcci4Zocwc+6KZ7wlAGaTl9+OevnzpIbS7K+hZHPW8udWh6ZHwlh1bEglkUqQN3Q7J+xwVb9RpArDVO9WETbJUz2/tqaUY4kUld5cn/TdfGk4K6SGYvbfpLP5TtJdxmA6P5HT5JRvqkuisS5p0FCsfxtSCYeZRmkoqk0Z576IcIGiC+QBgy9O+VDc5iX92bl9Vaotj26yOlG3kLSd43I+ybCj/eXNpMb1xbI51KcSaKwzm7zUcwg2eaUc50Zhb59zhQ33PJT6ZAL3//tpAPJaiLKCAMDY5npZzjA27+3HPcu22ceSCcIhY5uxo3vQNl+PBCxQCkQ12H7Zyb6wfi8A4D/PO9oePeidjhIEk3wmPf3l5e0AgHOPmYSGVMLxIQxlsmhrqsOktgas0TY+UuREPu7e1EnqHUmghmIYVe3tG8bYlno0ucwAYftz2+cZnPJ6HfYPZDDrigfx2+c2OMpOJghtjXXoHcwYtZiM1DSSCXJcQ3U4DakEmutTngX0hjM5tIWMIpX9/fk39zrSdYFqes5buyxzRFtjHfrTeZNZ33AG9ckEmur8BcpQOovGuiQSfhpKJq+huDcQG87m7Nnt7si1V7d1y3k5XjOOLugAYFBb5iOVINsUOmwwK6nn3Cy1GxOD6azth9GvrZ67MtW4fSz6c3aHDav32d5s1lDUe1EDN9MinQ0p82KJw5kc6pNJT1tXDNkCxaCh2FqXWYsM4v+Wb3P8req8s3sQPYMZNKQSmDvFMpUqgZLSnLZN9Uk01SXR1TeMd1/3lGO+SoIIs8Zbg+Dd+70Li5YLFigFohrshr19APINuS6ZsD98PUJJhQNPamvExaccYjc8ALh76Rbb+TyYzjk0lM17+3H74s3o7BvGEZNGYZO8nmLd7h6s2r5fG3V7G/Lx3/6b/dsdI5/Jaj4Uw0fW0TOICaMaPKGUYcuQKHTTj+001PJu3GPdzy+eejNfnuzQWhtTyOSEb2hoMkGeiCz10TfWJdHckPRoGcOZLEY1mTsyhXKof/6OFY70FVu6AFhmBJPJ6zoZVDGmpR5C5EfPe3uHMb613nqvPjb/Py3fhq1dA0iSeVZzOmcWKEOZ4JHx4g2dGM7msH3fALZ2DTgErOo0x7daG0kpzSCTtYI8lIZiMtOpc8e1NiCdFZ4Ouncogz29w3bHp3fQaqSs/CDu9+A0FzqvbUet+Tjl80u2ezt2pb3Up8x+n+FMDg11CTSmzGHDynFueh7qfDVYKUSgTBzlHGSq7/Br974CAOgZyiAlByTqubstFWNb6rG3b9gzSEoQ2c/Kb9WJcsACpUBaZIP93B+X495lVtimenGtDSm0NabwhuaYXbyhE4DlHB3bUo+ewYzdEX/1nlfs84isEZTqKD560wv2sfYmrwno3dc9DQC2huIX8aRwN/SBdDYfx+6eo5LNYX1HHyaOakSjS6CE+XUUb+zKa1S2o1XTInpsFT5vClQmPmWP7jH4UYbSWWudooQzyitv8kqgpT7lsR0PZ3O241TvGDbsyQvqGy+ZByD/PhW/ljsZTmprdIQF6/cGAOOkTVs9lz29QxjX2oDm+qRxcUdlFwcsE4V7/tKq7d0YTOdn/+sOdCUgVKektw9ds1Ptb9mmLk9epTWrQY8yk0URKBOkMHJ3wHcvteZlqQ5OF2RbOy1Nbs4ka9969zPRBwHu9qxCwmfLUbdbQKuVFaa2W/fkiJDUwqTdARdCWJpefTLhq83t6bVG+MOZnMHs6PSh5M27OfzuhU0eTVmnq9+53bDqF55+owOApakAlim21yVQrv/n4wFYAmV9h3OwCVjtSQ02/MzH5YAFSoGktLWevnz3ywCAj2rOsxNmjsGaXVaIX99QBlc/sBqA5RdQIzNTSGyCSI5krUZ12mFWZNkxU9vQ2pAydkiAeUQGeD909zW7B9L2Vqfu0c2aXT3Y3TOEd8wZj4a6hK3WA04NxW+uCAD8919W27/dJi81MgbyZibA+qBSyYTd8bs/hNd37setz29CY11COt4NJq+6BJrrk57l3ofSeZOX7nhVHy8AHDFpFN555ATMHNtsvKeTDhmDdbucpkd99De+xbovNRLf2zuMca31GN/agL193r3K92mjzZaGpG1GVXz0N4tleVa63pmpdzJRvkM9/FTX7E6ePRYA0FSf/9RV/VTkmRL+fUMZtDSkUJdMYFJbgx39pqPeybgW6ctwtTP1dSw89RAATg2lo3cI7U11aG+qw5jmOux0bRrm0FBcgxy1wdjhE1sBeE1eSihOH2tNPE47tLkcEmQJFM9S71q7MfmbOvuGkRPAoVKQudvVX/hSMoEAABpfSURBVFfuAKDPrbHyP/76bvzXn1/Fjx59AyaEEJ4tl3/y+Dp88OfPOc4BgFEN+Ym+Xf3DeP/xU/FPJ1grDY9tqceKLfs85Y9tqccs6fddvrnLc7xcsEApEOWc19HnGBwyrhmb9vZDCOGwkc4Y22wLlH39w54PMUGW7Vp1Gi0NKdSnErj/39+OloaUcT8MwBIGCfJ+gDu6nR+rvjzMc+v2YNPefkwYZY2e9ZHyYDqL8254FgBw1JQ2NKYs+78a2emdmt9ckYdf3Wn/bqrLR2MpwaVG8gDs0bB13HJQtzaYBcq//2E5AGu7gGSCHEJUdZKNdUm0NDg1lL+t2om9fcOej95Ec0PK14k5a1wLduwfdIyO9Y5ctYPeoQyGMlms3NaNnd2DGNtSj06DQHlEe04tDSnHCH04k7NHo2fIxRr1Z69GvmpQoL8LVc5//+Mx+NLZRwBwChlVruqc9w2k8crWfXh9Z4/tVzlsQiv+8soOT52V9jVODgrcI/pvyYHEyYeOk9fNH1cmQMDSAvcPOJ/zG9JPOLW90WPyUhuMTW5vRDJBHqf889KXqQYNen6leZn8VKqNjWpIoak+4TGVKjOo0ozcA7tHVlkRmuo9KPOu+v78zE0dvUOesnqHMnhp8z5bI/3m+48BYPmN9vWnkcsJbN83gGmjm+w8c+Q71PnUGYfiU2ccisntjZg9vgXPrN2D59btCdSW4oIFSoG0NqTwvuOm2CNDAOjXGsbMsc3oGcxgX3/aHlWpaBsVoXLfiu12XLkiYZu8rA+hbyiD9qY6GfqZcsyP0NX2/uEM6pIJdLrUZ7/FFYH8yDebE7YNVvHs2vzKprPHtdiNfo+MSFEd9ejmOkdnqhBC4NO/XwbA0srqkmRrBOreJmgaytTReTtyz2AabY11tvnAHVapd05tjXV2x7indwjf++vrAKyRWXN93ofS1TeMy363DNmcQEt9CkRODcV9D/sH0nizow+7ewYxMJzFGT94AgBw4szRmDm2GUIA2/flhbXqyC/6hxm2f61/OIsd8pwZY5sxub0R+/rTjndyy3Mb8L/S93LnZac4tNDBdBb/8N1H7fDmD8+bAcBpPlLPYoJBQ1GDj5aGVD6IRMu7p9d634dPbEVdkrChow+f+p31zpTQ+vube5HNCc9Okb1Dacus2OBdAkU3tU2Wpie9zh3SBAg4zTiKe1/aivamOswc1+wYIPUNZbB8SxfqUwmMakhhVGPKI4x+8aTli1NC3SF801k01afkMjROQaW07FGNdRjTXI/9A2mH5qvu3yRQ9PtV37Ya5Nwhl2SqTybw8Ks7cdo1j2P19vwWBkukKfIfj58KN4PpHL5yzhH2Nae2N2HV9m7s6hlEOiswfUxeoJx33BRP/q+de5S9zfL0MU1YvKETH/3NYmwJ6BPiggVKEZx+xATs1j60uVPb7d9Kg9nU2Y9t+waQTBCWf/NsAHln5I8fW4uP37LEUWYyQWiuT2Jf/zCGMzn0DmUwSo7UW+XH2zucV3sV7z56EoYyOfxh8WY7OunVbd24+Car/MX/cRYuPGm6PerX6R/OYlxrg20jBpwfTHtzHVZu7QaQj3x6Vi6lfej4FnT1D3s6BV043fOZtzkE1j3S59So7buha1LdA2m0NdXZtt/Xd+Y/wJc2d9nRVOfMnYRxrfXY0T2IV7d14yt3v4yV26x6jmm2BIoyH+kjzmTS0gL1TmWfS6AcNsEa8T3+2m58+vfLbNv9hfNm4NAJ1rt9Xdvb5KePrwMAXHzqIXbn3TeUsZ/px045BHOlaent1z5h59NNgicfOs4SKFI4LdvUZQvL+bPHYnJ7I1IJwkYZmLG3dwjv+4mlRc4a14JkghydhZpJ3dqQtMOKda2ro2cILfVJjG6ux5GTR+HNjl572aAL5013PI+d2vt5s6MXv35mA0Y1prQ1tfICY/aVf7V/K7Ol0mC+df8qLNnQaQ/EWupTjgHDqu3deGNXL7oH0hjX0oAtXf32wOnah1/HX1fuRGMqASLC2OZ6xwBKmZ0A4BBprtQ1tu7+NFobkmiuT3pG6eq8UY0pjGmuR044hbOKrjxUtot9mg9Ntd1vnHe0rWkPprPY3TOI12QbufX5Tfj075dh274B/OwJq63kcsLWti87/VC8fvUCW1tU6FGhJ8wcja7+tB0ROk0TKMdNHw03+mRrXTArs105qWmBQkQLiGgNEa0joitG6rpnShOEYsGxk+3f6sP8zTPrsX3fAE6YMdp2POvmsi3SQfk+OcJIEOHtc8ajbziLFzd22vZsIO8kPu5bf8Mbu3rQITurqz9wrD1hEgCWy2gkfcHKSW2NmNjWgN6hDN7s6HWYib58zhEY11Lv0JZUB/rol063zwHyIz/lE3rHnAnICaez9zfPrMe87zwKAPjpv5yAwya0Yv9gBg+8sgN7eofwh8VWWONH589EQyqB+mTCFhIA8MzaPSCCPQL7zoOv2R3WB3/+d/u8Gy+Zh7HShv++nzxrL0sBWEK7uT5vPtJH5ks2dCIngEdW7bRDLN1RM1e85ygQAVf8aSWe0vwrR09pwzFT25FKED5z+0u4b8U2fPDnz+HR1yyTx+ETWx0aylL5XMa31jv2ydjbaw7h1M2auoA/7y1TUJdMYPqYJmySE0yf08Kap49twoTWBrvj375vwO54VEAA4NQUXtrcZYeMzxjTjC1d/XhhvTViftdRkwBYWhPg1HQ/KzvBPb3DdptWAmPx+nydFp56SP64vK4yuT4gzWijGlNYsrETp37Pmil+2W2WhjSprQGnHjYOu/YP2T6W37+wCUC+8x/TUm8vwjmUyeLfbn8JgBWFN7GtEa0NKeyVWtir27rx4ModOGZqu6NdAMCWzn584GeWz6Ktqc42xerCauW2bqQShFMPs0x4G7VoS6VxnDBzjN1mN+7tN5o3gXxbe0nzaUxptwJf1u12+qscZi0ZxPA/Ugufrh3ThceNF5+EP/7rKY5yVOjwu4+ehFSy/N19zQoUIkoC+BmA9wCYC+AjRDR3JK49sa0Rn3vX4cZjyqH7wCs78ML6TscIWHWCOm+Tzvfjprdj3iFjAFgmqSfWdNg2fzVqBoBzrn8aP5Gj4iMnjQIR4d7PvA0AsK1rAEs35q/5wROnAQBOnm19DO/98TN2uO53PnAspo9pxsRRDVjf0Webul7Zug8zxzbj8IlWI1Zx///z4GtYtb3brscnT5sNAHhqTQcG01kMZ3K496W8z+i9x1qCUjlHH1mV9xccNrEVa77zHnzx7CPQ2TeMnd2D+NfbrNnqyzfvs01egGWCM4Vy6s9SmSUOHd+ChlQSrdKJmcnm8Lk/LrfP6x5IYziTw8a9/bjyTyvxwvq99kf+7NffCcDywcwy+MkmjGpAfSofGv75O1bgpc15Z2hDKmlHAP7q6fV4U3YQh01oxVStAzjpO4/izhfz8wU+cdosAFbHsr17EJv39tt1OumQMfgn+Q6HMzk8+MoO5HICL2tO2KMmt2Fcaz1WbuvGmx29DrPKCTNG25M1N3f2o3cogwde2Y4VW/bZnd70MU3GKKG3TG9HMkF46o0OvNnRi9d37re1HSBvAnpxoyWIlmoDi6vef4z9LHqHMo4VmFW7Uc9xR/cglm/usgc1v7/0ZMyQ39CKLfscS/384AJrNe+xLfW2AP310+vtsp/66jvt42r2uDLvThjVgNHNdVjf0YfdUlDd9vxGO++MMc0YIycKKp/i1q5+PPDKDpxxxAQcMrYZjXUJ+70OprP4f7LNHjGpFVPaGzGqMYXV2/ejUwozfYoAYGn3tzy3ARfIvecvO/1Q2wT4Hm1QCji1kLcfPt5xbJZL0/jZv5yI3106H+ccM9kWfIpvvm8ufvKRE/DrS07CSFDLqw3PB7BOCLEeAIjoDgDnA1gdmCsmvvDuI3DD4+tsDUPRWJd0zLg9Uo4uFL/82Em2jwGwOv23zhiNoyaPstdQUigtZt6ssfi3Mw/Dz6Wd+EE5ylP285OkILp98WbcvngzprQ34pipbfjfC63QQrUD31Amh7Ovt8KNVYP97FlzcMeLW/Cxmxbb11XlAfnZuGt399rO+jsvO8UWNDc/twE3P+eMVnn4C++w7+WWT/wDLvjl8/jP/3sVAPDny0+z1fl3zBmPax8GTvlefj2juz51KgDgq+ceiR88ssb+aBXf++BbAHhDe8+ZOwk/++iJAICZUks8/D8fcpzzxXcfgc9qAuaiG63Q7Ka6JKaPyUd2uf0G1//z8faI0bRiwepvnwsgH7328pZ9eHnLPrQ31cn2kMQFJ023TX5fv3clAOCHFx6PC06yTEyqEz39B3mz2N2fOtV+jmqB0UP/46+YISOZLpGRVKlkAqu278NZ//uUnXfRF0+3IubkqPTnT75ptx8AuPp8y+E7f/Y4Oyz6awuOtI8316dw/PR23PHiFtsfoPjP9x6NGWObMamtAd9/eA2+//Aa+9iCYyYjkSDbvPvtB1bj21Krvfydh+Gr5x4FwOlM/iepfR45aRQOn9hqr6qrNA8AOGryKFyo+ZLe7OjDxTctticWP/mVM23BPZDO4s8rtuPPK7bb+S87/VCs2dWD2xdvxlnXPYWjJ7dhiRSGzfVJTG5vtDXDC375PD531hzc8NhaAMCJh4xBIkGYMaYZv3l2A2aOa8YKbTChBkAnzhyDPy7ZbGu/v/3EfDy0cgc+c+ZheOy13fjava/Yps4JoxrwH+892i7jug+/FZe/sxfXL3oDizd0OgYhyQThr597Bz53x3J8+ewjHMuvAGY/iqKpPon3G/w05aJmNRQA0wDoLX2rTBsREgnCK986B9f/81s9x5752rvs3998v1NpWnDsZPyz/DBu+fg/oLEuiblT2+yO45mvvdM+9/eXnmz//tqCo/C3L55u/93eVIcZ2ihmUlve0b2jexAnzBztWEZ7kZYXAA6XWs+00U12p6a47PRD7d8T2xptM57irTMtu62KINK5+JRDcNTk/EKIJ8wcY+dvb6rDUZPzAvaoyaMcIcQfPXkm5ssw10+dfqjts1A8+ZUz8ZH5Voj2qYeNw4c1e/8HT5xmf2gnzvTalTdecx7ef/xUvH71ArtDVugCFAA+f9Ycx98qRBPIx/8rvnHe0bbvRM1cVqhOFbBG1+6R5tlzJ9m/57nqAMAxwPjKOflnvaVzAAuOmYz//kdLKJzmGpUCcHRIF59yiOPYlPZGXHzqLLsO33r/XHzoxOn4tzOdWvevLp7nKffBz70d/yrbx2fOOMxx7IKTpuOGj5wAwFpscZxLI/+c9lxPnOm937989u0gIk97SyUIt31yvv23uv4za/cgnRV4x5zxjlH7t+VzUfz6knmYOroJ7zxyIgArsksJk5Nnj8XzV5wFwDnyV8Lk42+bhUvfbmlVR8i2+837VuFPy7dhxtgmezABOLWMWeOacdz0dnzjfXMxrrUBHzxxGo6b3o5RjSlc9f65ju8csNrOsdPa8ZN/OQFPffVMj9CYO7UNj37pDLznLf7CoxqgoF33qhkiuhDAuUKI/yf/vhjAfCHEZ13nXQbgMgCYOXPmSZs2bRqR+i3d2InOvmGcc8xkz7GcXDZDn9QXlY6eIazd3WObyhS9Qxnc/OwG7OgexKa9fbjhIyc4OmvAikq56dkNOHn2OLxlervjWCabs9aTEsLTmAFrZL6ndwiHjm9xCKquvmF8475XMZTO4esLjsRhE1o9mlYmm8PNz23Au46aaJvSFD2Dafzf8m14eUs3/ut9R2N0c72j7HUdvVi0ehcuPGm6bUvW2by3H5s7+/H2Oc7noUJ7X9zQhdMOH+fZo2IwncVp1zyOD544DV86+0iHOUfd7xfuWI4vnn2Ex/GZkQt67ugexEdPnul4j8OZHN7s6MUzaztwxhETceRkb51391hLa+imTFWnNTt78MSa3fjEabM9WljvUAZn/uBJDKWzePpr78QYrcPetX8QP3p0Le5bsQ3Xffh4LDjW2fE8tHIHtnT1Y+6UdrxlerunbD+27RtA/1AGdy/bitMOH2+HMANWe9rbN4z7VmzHMVPbcMqhTsHWPZDGls5+pLM5HD2lzfGchBB4ZNVOpBIJ7O4ZwplHTnAIwVxOYGvXgBXRNZj2hOtncwJv7OrBn5dvw/lvnWYHPujl/9/ybZg6uslRr0Wrd+HR1btw9JRReNdRkzCpvcGOiAKAv7+5Bw+t3Il3HjUBzfUpR96B4SzufHEzVm3fj8MntuK846Y4NNtsTuCBV7ZjSnsTDpvQYpuzdHI54fk+agEiWiaE8I4w3OfVsEA5FcC3hBDnyr+vBAAhxPf88sybN08sXbrU7zDDMAxjIKpAqWWT14sA5hDRbCKqB3ARgPsrXCeGYZiDlpp1ygshMkT07wAeAZAEcLMQYlVINoZhGKZM1KxAAQAhxF8B/DX0RIZhGKbs1LLJi2EYhqkiWKAwDMMwscAChWEYhokFFigMwzBMLLBAYRiGYWKhZic2FgMR9QDYCaA74LT2gOMzAWz2ORaWN+hYqXmD6lXO6/KziidvueoUdpzfXzx5D4ZndaQQwrvsgxshxEHzH4ClAG4MOcf3OICOEvKWct2wvL71KvN1+VmV+VmVUm4pdeb3x8/K9ffSoPPVfwejyesvJRz3bt4cPW8p1w3LG1Svcl6Xn1U8ectVp7Dj/P7iyXuwPStfDjaT11IRYT2acuUvF9VYr2qsE1Cd9eI6Raca61WNdQLirVfUsg42DeXGCucvF9VYr2qsE1Cd9eI6Raca61WNdQLirVeksg4qDYVhGIYpHwebhsIwDMOUiYNeoBDRzUS0m4he1dKOJ6LniWglEf2FiNpkeh0R3SrTX1N7sMhjTxLRGiJaIf+bOEJ1qieiW2T6y0R0ppbnJJm+johuIPcuU5WpU2zPSZY3g4iekO9jFRF9XqaPJaJFRLRW/jtGy3OlfCZriOhcLT2W5xVznWJ5XoXWiYjGyfN7ieinrrLibFdx1qtSz+psIlomn8kyInqXVlYln1VQvWL9Dm2ihIIdyP8BOB3AiQBe1dJeBHCG/P1JAFfL3/8C4A75uxnARgCz5N9PAphXgTpdDuAW+XsigGUAEvLvJQBOBUAAHgLwniqoU2zPSZY3BcCJ8vcoAG8AmAvg+wCukOlXALhW/p4L4GUADQBmA3gTQDLO5xVznWJ5XkXUqQXA2wF8GsBPXWXF2a7irFelntUJAKbK38cC2FYlzyqoXrE8K08d4y6wFv8DMAvOjnI/8v6lGQBWy98fgRVOlwIwTr7QseV4QQXU6WcAPqad9xiA+bLxva6lfwTArypZp3I8J0Md7wNwNoA1AKbItCkA1sjfVwK4Ujv/EfnBx/68Sq1TOZ9XWJ208z4OreMu53MqpV7V8KxkOgHYC2twUBXPyl2vcj6rg97k5cOrAP5R/r4QVmcJAPcA6AOwA9YM1B8KITq1fLdI9fG/SlFtC6zTywDOJ6IUEc0GcJI8Ng3AVi3/VplWyTopyvKciGgWrFHZYgCThBA7AED+q1T6aQC2aNnUcynL8yqxTopYn1fEOvlRtnZVYr0UlX5WHwKwXAgxhOp6Vnq9FLF/hyxQzHwSwOVEtAyWajks0+cDyAKYCss08WUiOlQe+6gQ4i0A3iH/u3iE6nQzrIa6FMCPAPwdQAbWiMRN3CF9hdYJKNNzIqJWAPcC+IIQYn/QqYY0EZBeyToBMT+vAurkW4QhreR2FUO9gAo/KyI6BsC1AD6lkgynjfizMtQLKNN3yALFgBDidSHEOUKIkwD8EZZNG7B8KA8LIdJCiN0AngMwT+bZJv/tAfAHWMKn7HUSQmSEEF8UQrxVCHE+gNEA1sLq0KdrRUwHsL3CdSrLcyKiOlgf2O1CiD/J5F1ENEUenwJgt0zfCqe2pJ5LrM8rpjrF+rwKrJMfsbermOpV0WdFRNMB/B+AS4QQqr+o+LPyqVfZ+isWKAZUxAMRJQB8A8Av5aHNAN5FFi0ATgHwujTtjJd56gC8D5Y5qOx1IqJmWRcQ0dkAMkKI1VL17SGiU6Q6ewksm2vF6lSO5yTv7SYArwkhrtMO3Q9gofy9EPl7vx/ARUTUIM1xcwAsifN5xVWnOJ9XEXUyEne7iqtelXxWRDQawIOw/GDPqZMr/az86lXW/ipup0yt/QdrZL0DQBrWiOJSAJ+H5XB/A8A1yDueWwHcDWAVgNUAvirTW2BFMr0ij/0YMkpnBOo0C5ZT7jUAjwI4RCtnnmwobwL4qcpTqTrF/ZxkmW+HZUZ4BcAK+d97YQVNPAZLM3oMMnhC5vlP+UzWQIu6iet5xVWnOJ9XkXXaCKATQK9853PL0K5iqVclnxWswVSfdu4KABMr/az86hXns3L/xzPlGYZhmFhgkxfDMAwTCyxQGIZhmFhggcIwDMPEAgsUhmEYJhZYoDAMwzCxwAKFYaoEIvo0EV1SwPmzSFv9mWEqTarSFWAYxppsJoT4ZfiZDFO9sEBhmJiQC/Y9DGvBvhNgTfi8BMDRAK6DNTF2D4CPCyF2ENGTsNY5Ow3A/UQ0CkCvEOKHRPRWWCsPNMOaFPdJIUQXEZ0Ea620fgDPjtzdMUw4bPJimHg5EsCNQojjYC3vfzmAnwC4QFhrnt0M4Lva+aOFEGcIIf7XVc5tAL4uy1kJ4CqZfguAzwkhTi3nTTBMMbCGwjDxskXk1036PYD/gLW50SK5QngS1hI2ijvdBRBROyxB85RMuhXA3Yb03wF4T/y3wDDFwQKFYeLFvZZRD4BVARpFXwFlk6F8hqka2OTFMPEyk4iU8PgIgBcATFBpRFQn96fwRQjRDaCLiN4hky4G8JQQYh+AbiJ6u0z/aPzVZ5jiYQ2FYeLlNQALiehXsFZ//QmsLX1vkCarFKxNx1aFlLMQwC+JqBnAegCfkOmfAHAzEfXLchmmauDVhhkmJmSU1wNCiGMrXBWGqQhs8mIYhmFigTUUhmEYJhZYQ2EYhmFigQUKwzAMEwssUBiGYZhYYIHCMAzDxAILFIZhGCYWWKAwDMMwsfD/AcGA/IkB7OeIAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ + "sorted_data['inc'] = sorted_data['inc'].astype(int)\n", "sorted_data['inc'].plot()" ] }, @@ -215,9 +3232,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 12, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXl4XOV59/+5Z9UuS7Isy5Z3G4zN4g1jthRCgoEkP0gDjZO3gba0TmnapE3eNCFv+5KS8mvSJiGhKUlIoBCSAi4hhSyEGAg7eMHYeMX7IluWZGtfZkYz87x/nHNGI2kkjaSRZtH9ua65ZvzMOUfnHJ8533Mvz32LMQZFURRFSQZXundAURRFyR5UNBRFUZSkUdFQFEVRkkZFQ1EURUkaFQ1FURQlaVQ0FEVRlKRR0VAURVGSRkVDURRFSRoVDUVRFCVpPOnegVQzdepUM3fu3HTvhqIoSlbx9ttvnzHGVA63XM6Jxty5c9m6dWu6d0NRFCWrEJFjySyn7ilFURQlaVQ0FEVRlKRR0VAURVGSRkVDURRFSRoVDUVRFCVphhUNEckTkc0iskNEdovIP9njXxWRkyKy3X7dELfOnSJyUETeE5G1ceMrRWSn/d19IiL2uF9EnrDHN4nI3Lh1bhORA/brtlQevKIoijIykkm5DQLvN8Z0iIgXeE1EnrW/u9cY8834hUVkCbAOWArMAJ4XkXOMMRHg+8B64C3gN8B1wLPA7UCzMWahiKwDvgF8XETKgbuAVYAB3haRZ4wxzWM7bEVRFGU0DGtpGIsO+59e+zVUj9gbgceNMUFjzBHgILBaRKqBEmPMm8bqMfsT4Ka4dR6xPz8JXGNbIWuBjcaYJlsoNmIJjaIoSkbQEQzz1LZaJkvr7KRiGiLiFpHtQAPWTXyT/dVfi8i7IvKQiJTZYzOBE3Gr19pjM+3P/cf7rGOMCQOtQMUQ2+q/f+tFZKuIbG1sbEzmkBRFUVLC73af5vMbdrC3rj3duzIhJCUaxpiIMWYZUINlNZyP5WpaACwD6oBv2YtLok0MMT7adeL37wFjzCpjzKrKymFnwSuKoqSMrlAEgP31KhoDMMa0AC8B1xlj6m0xiQI/Albbi9UCs+JWqwFO2eM1Ccb7rCMiHqAUaBpiW4qiKBlBKBwFVDRiiEiliEyxP+cDHwD22TEKh48Cu+zPzwDr7IyoecAiYLMxpg5oF5E1drziVuDpuHWczKibgRftuMdzwLUiUma7v661xxRFUTKCUGRyiUYy2VPVwCMi4sYSmQ3GmF+JyKMisgzLXXQU+DSAMWa3iGwA9gBh4DN25hTAHcDDQD5W1pSThfUg8KiIHMSyMNbZ22oSka8BW+zl7jbGNI3heBVFUVJKr6XRMcySucGwomGMeRdYnmD8U0Oscw9wT4LxrcD5CcYDwC2DbOsh4KHh9lNRFCUdOKJxormLrlCYAl/OFQ/vg84IVxRFGQOOe8oYONiQ+9aGioaiKMoYcCwNmBwuKhUNRVGUMRAMRykr8OLzuCZFMDy3nW+KoijjTCgcJd/rprzQx4mmrnTvzrijloaiKMoYCEWi+DwuivweOkOR4VfIclQ0FEVRxkAoHMHncZHndRNQ0VAURVGGIhS2LI0Cn5uunnC6d2fcUdFQFEUZAz0Rg8/tosDnidWhymVUNBRFUcaAY2nk+9x0q2goiqIoQxGMRPF53JZ7SkVDURRFGYpQOIrP7SLf66a7R0VDURRFGYJQOILfdk+FwlEi0dzu4KeioSiKMgaceRoFPjcAXaHczqBS0VAURRkDMfeUXd0214PhKhqKoihjIDZPw+tYGioaiqIoyiDET+4Dcj4YrqKhKIoyBpyYRp5PLQ1FURRlCKJR0zsj3HZPaUxDURRFSYjTtc9yT1mB8EmfPSUieSKyWUR2iMhuEfkne7xcRDaKyAH7vSxunTtF5KCIvCcia+PGV4rITvu7+0RE7HG/iDxhj28Skblx69xm/40DInJbKg9eURRlLDii4czTAI1pAASB9xtjLgKWAdeJyBrgy8ALxphFwAv2vxGRJcA6YClwHXC/iLjtbX0fWA8ssl/X2eO3A83GmIXAvcA37G2VA3cBlwCrgbvixUlRFCWdOK1e+87TmOSiYSycxrde+2WAG4FH7PFHgJvszzcCjxtjgsaYI8BBYLWIVAMlxpg3jTEG+Em/dZxtPQlcY1sha4GNxpgmY0wzsJFeoVEURUkrjmh43XHZU5NdNABExC0i24EGrJv4JqDKGFMHYL9PsxefCZyIW73WHptpf+4/3mcdY0wYaAUqhtiWoihK2olZGm6rCROoewoAY0zEGLMMqMGyGs4fYnFJtIkhxke7Tu8fFFkvIltFZGtjY+MQu6YoipI64gPhfo8Ll2ggvA/GmBbgJSwXUb3tcsJ+b7AXqwVmxa1WA5yyx2sSjPdZR0Q8QCnQNMS2+u/XA8aYVcaYVZWVlSM5JEVRlFETH9MQkUnRiCmZ7KlKEZlif84HPgDsA54BnGym24Cn7c/PAOvsjKh5WAHvzbYLq11E1tjxilv7reNs62bgRTvu8RxwrYiU2QHwa+0xRVGUtBOMEw1gUjRi8iSxTDXwiJ0B5QI2GGN+JSJvAhtE5HbgOHALgDFmt4hsAPYAYeAzxhjnLN4BPAzkA8/aL4AHgUdF5CCWhbHO3laTiHwN2GIvd7cxpmksB6woipIqHEvD77ZEo8CX+z01hhUNY8y7wPIE42eBawZZ5x7gngTjW4EB8RBjTABbdBJ89xDw0HD7qSiKMtHExzQA8r25371PZ4QriqKMktAkdE+paCiKooyS/qJh9QnX7ClFURQlAaGIZVX43I57SrOnFEVRlEFIZGnkeiBcRUNRFGWUJBQNtTQURVGURARjKbdWCZE8r4qGoiiKMgj9U24LfG66eiJYc5NzExUNRVGUUZLIPRWJmpiY5CIqGoqiKKMkFI7idglul1VbNd/u3pfLLioVDUVRlFESCkdj6bZAb0+NHM6gUtFQFEUZJT2RaMw1BUyK7n0qGoqiKKMk1E80Yo2YVDQURVGU/gQHcU+ppaEoiqIMIBSO4k/onsrd+lMqGoqiKKMkFO7rnsr3avaUoiiKMgj9YxqaPaUoiqIMSv+U23yNaSiKoiiDMcA95dPsKUVRFGUQBrinvGppKIqiKP14ZX8j133nFRrbg33cUx63C5/bRVfPJM6eEpFZIvJ7EdkrIrtF5HP2+FdF5KSIbLdfN8Stc6eIHBSR90Rkbdz4ShHZaX93n4iIPe4XkSfs8U0iMjdundtE5ID9ui2VB68oijIath5rZt/pdupaA30sDbBcVIEctjQ8SSwTBr5gjNkmIsXA2yKy0f7uXmPMN+MXFpElwDpgKTADeF5EzjHGRIDvA+uBt4DfANcBzwK3A83GmIUisg74BvBxESkH7gJWAcb+288YY5rHdtiKoiijp741EPvcXzSsPuG5KxrDWhrGmDpjzDb7czuwF5g5xCo3Ao8bY4LGmCPAQWC1iFQDJcaYN41VbP4nwE1x6zxif34SuMa2QtYCG40xTbZQbMQSGkVRlLRR3x5gfmUhM6fkM6M0v893+V6rp0aukoylEcN2Gy0HNgGXA38tIrcCW7GskWYsQXkrbrVae6zH/tx/HPv9BIAxJiwirUBF/HiCdRRFUdLC6dYA86cW8r1PrugzIxws95RmTwEiUgT8HPhbY0wblqtpAbAMqAO+5SyaYHUzxPho14nft/UislVEtjY2Ng55HEp20hOJ8rnH32H3qdZ074qi0NAepKokjzyvGzs0G8NyT03iQDiAiHixBONnxpinAIwx9caYiDEmCvwIWG0vXgvMilu9Bjhlj9ckGO+zjoh4gFKgaYht9cEY84AxZpUxZlVlZWUyh6RkGfvq2nl6+yneOtyU7l1RJjnBcISmzhBVJXkJv8/3eSa3pWHHFh4E9hpjvh03Xh232EeBXfbnZ4B1dkbUPGARsNkYUwe0i8gae5u3Ak/HreNkRt0MvGjHPZ4DrhWRMhEpA661x5RJxt66NqC3vaaipIuGtiAAVSX+hN8XeN05XUYkmZjG5cCngJ0ist0e+wrwCRFZhuUuOgp8GsAYs1tENgB7sDKvPmNnTgHcATwM5GNlTT1rjz8IPCoiB7EsjHX2tppE5GvAFnu5u40x+qg5Cdlji0YwnLs/RiU7aGi3MqcGszRyPXtqWNEwxrxG4tjCb4ZY5x7gngTjW4HzE4wHgFsG2dZDwEPD7aeS2+w5pZaGkhmcbnUsjcSikaeBcEVJL9GoiVkaKhpKuqlvG8bS8Oa2paGioWQ8tc3ddAStbJSgioaSZurbAvjcLsoKvAm/L/BZMY1odECiZ06goqFkPI6VAWppKOmnvi3AtBL/gFRbh3yf5fXP1QccFQ0l49lT14ZLoLLYTyiSmz9EJXs43RYY1DUFud/yVUVDyXhOtXQzvSSPkjyPWhpK2mloCzJ9CNHIz/Hy6CoaSsbT3RMh3+fG73Fryq2SdurbAlQWJ56jAXGNmHJ0rsaIak8pSjoIhCLked143a6c9RMr2UFXKExnKDKkaBTkeMtXFQ0l4wmEI+R73bhcou4pJa2c7QgBUFmUhKWRo6Kh7ikl4+m2LQ2/Ry0NJb00dlgT+6YW+wZdpsDOnurO0e59KhpKxtPdE42JhloaSjo5026LxhCWhrqnFCXNBO1AeCQa1ZRbJa2c7bTcU0OJRq5nT6loKBlPd0+EPI+LcFRjGkp6cSyN8sLB3VMa01CUNBOwLQ2f26Upt0paOdMRpDjPQ55tTSQi191TKhpKxtPdYwXCfRrTGFfCkWismrCSmDMdoSEzpwDyPLk9T0NFQ8lojDEE7EC4isb4EY0aPr9hBzfc9yqnWrrTvTsZy5mO4JDxDACXS8j3uunWMiKKMvE4Kbb5mnI7rnznhQM8s8PqpFzbrKIxGGc6gkOm2zrkciMmFQ0lo3GCiXleFz6Pi3DU5GzJ6XTys7eOsWhaEdDbL0IZyJmOEBWFQ1saAHne3G3EpKKhZDQBO/Cdb7unAE27HQfaA2FWzS0HVDQGIxSO0trdM6x7CtTSUJS00WtpWAULIXf7FKSLYDhCKBJl5pQ8/B4XDXZaqdKXJmeORrLuKQ2EK8rE42Sg5MVZGpp2m1o6g9b5LPJ7qCrJU0tjEM50DD8b3CHP6yYwWUVDRGaJyO9FZK+I7BaRz9nj5SKyUUQO2O9lcevcKSIHReQ9EVkbN75SRHba390ndusrEfGLyBP2+CYRmRu3zm323zggIrel8uCVzCfQY1kVeV4XfrftnlJLI6V0BKwsn6I8L1UlfhWNQYjVnSoa3tLI87oJTlbRAMLAF4wx5wFrgM+IyBLgy8ALxphFwAv2v7G/WwcsBa4D7hcRZybM94H1wCL7dZ09fjvQbIxZCNwLfMPeVjlwF3AJsBq4K16clNzHeVrL97rxe1U0xgOn/3qR3820kjwa2tQ9lQinwm1yloYr9sCTawwrGsaYOmPMNvtzO7AXmAncCDxiL/YIcJP9+UbgcWNM0BhzBDgIrBaRaqDEGPOmMcYAP+m3jrOtJ4FrbCtkLbDRGNNkjGkGNtIrNMokICYa9oxw0JhGqukVDS9Vxeqe6o8xho//8E3uenoXkJxo5HvdOTu5b0S1p2y30XJgE1BljKkDS1hEZJq92EzgrbjVau2xHvtz/3FnnRP2tsIi0gpUxI8nWEeZBCSKaailkVo6HdHI81BV4qczFKEjGKbIr6XpwLoGNx1pYsXsKXzowhkUJnFeJnVMw0FEioCfA39rjBmq1oAkGDNDjI92nfh9Wy8iW0Vka2Nj4xC7pmQbTvZUflz2lKbcppb2OPdUld37Wq2NXpxr8KblM7n9inlJrTPpRUNEvFiC8TNjzFP2cL3tcsJ+b7DHa4FZcavXAKfs8ZoE433WEREPUAo0DbGtPhhjHjDGrDLGrKqsrEzmkJQsIWBbFX57ch9AMEd9xekiFgj3e5lmtzFV0eilKy7tO1ks0cjN6zSZ7CkBHgT2GmO+HffVM4CTzXQb8HTc+Do7I2oeVsB7s+3KaheRNfY2b+23jrOtm4EX7bjHc8C1IlJmB8CvtceUSUIglGhyX24+waWLePfUNNvS0GB4L47F4FSvTYY8r4tQJEokB6sXJOO0vBz4FLBTRLbbY18Bvg5sEJHbgePALQDGmN0isgHYg5V59RljjPMrvwN4GMgHnrVfYInSoyJyEMvCWGdvq0lEvgZssZe72xjTNMpjVbKQQHxMQ1NuxwXHPVXgdVNVopZGf7riHlySxbFKguFIrP1rrjDs0RhjXiNxbAHgmkHWuQe4J8H4VuD8BOMBbNFJ8N1DwEPD7aeSm3T3RPC4BK/bFUu51eyp1NIRsILeLpdQ5PdQ4HPrrPA4untGLhrOst2h3BMNnRGuZDTdPZHYD1BTbseHzmCYQr91jkWEacV+FY04YskYI3RPQW9MLpdQ0chyeiJR7t24n/ZAT7p3ZVwI9ETx26Lh15TbcaF/em1xnjdnr6fR0N0zGtGwls3FDCoVjSxn27FmvvvCAV56LzdTja1Wr9ZlGku5VdFIKe3BMEV53ti/i/yeWHBc6bU0CrzJu5mca1VFQ8k4nHo4jTnqTgj0RGLtM3sLFqpopJLOYJgif+9TdKHfQ0cw9252o8WpVpvnS/526VglKhpKxuGkRuaqD7q7JxL7AeqM8PHBCYQ7FPndamnE4aR9jySgnWdfq7k4V0NFI8txxCJXLY3uUCTmH3a7BI9LdJ5GirFiGr3uqUJ1T/UhNrnPk/ztUmMaSsbiiEVDe27m1QfC0T4zcX0el1oaKaajn3uqyO+JFTFULGvX53bhcY/cPZWLRQtVNLIcRyxy1dIIhCLke3svU5/HpTGNFGKMsUQjr9f1Uuj3EAxHCWuNLwC6Q+ERZU4BsTicuqeUjKMxx91TgXCkj6XhV0sjpQTDVqmL/u4p6O3oN9mJnyuULLF5GmppKJmGIxZNXSF6cvDJsDvU9wer7qnU0h7orXDr4HzuCKmLCqyYxkjqTgHkafaUkon0RKKc7QxRWezHmN7G97lEd09fS8PnVvdUKumIK1bo0GtpqGiAnfY9UktD52komYjT6H7pjBIgNyuTBnv6B8LdKhopxBGGQt9A0dBguMVoLA2vW3CJxjSUDMNxTTmi0diRWxlUkaghFIn2cU/5PS5twpRCYu6pvPh5GmppxBM/VyhZRIT8HG3EpKKRxTiWxfkzSvv8O1foLYveN3sqFM69H2K6cIShOD4Q7lPRiCd+rtBIyMvRPuEqGlmMU0LkvGrb0sixDKpEheL8mnKbUhwXVGG/eRrWd7l3wxsN3T0jd09B7nbvU9HIYhzLYsaUfKYUeHOulEh3bCauptyOF739weNjGtb5VkvDon8GX7LkeV0EctAqVtHIYhraA5QVePF5XFQW+XPO0ugMDfS3a8ptamnrtkqgl+QPnKehgXCL7tDIYxpgWRpBdU8pmURje5BpxVZP52kl/pwrJdLWbd20SuLKdmvKbWpp6gxR6HMPmEDpcYlaGjajmdwHGtNQMpDmrhBlhdYNtazAR3NXbjXOaY09BfdaGn6PWy2NFNLcGaKs0NdnTES0aKFNKBwlHDWjjGm4NKahZBbtgTDF9lN4LjbOiblO4i0NTblNKU1dIcr7iQY4RQtz7yl5pHTHMvhGLhqTNuVWRB4SkQYR2RU39lUROSki2+3XDXHf3SkiB0XkPRFZGze+UkR22t/dJyJij/tF5Al7fJOIzI1b5zYROWC/bkvVQecK7YEwxbb/ORefDNvslqOlcf52v8eVkz/EdNHcGaKsYKBoFGpPDWB0/cEd/JPYPfUwcF2C8XuNMcvs128ARGQJsA5Yaq9zv4g4Z/v7wHpgkf1ytnk70GyMWQjcC3zD3lY5cBdwCbAauEtEykZ8hDlMRzBMcV6caIQiRKMmzXuVOpyYRnFcIHxKgZeuUIRgDmalpIPBLI1CLY8O9Foao3JPedwEJ6N7yhjzCtCU5PZuBB43xgSNMUeAg8BqEakGSowxbxpjDPAT4Ka4dR6xPz8JXGNbIWuBjcaYJmNMM7CRxOI1KXFKWve6p6yLuiuHnmxau3so9Ln79DFw/O8tORa/SRfNnT0JLQ3tqWHRZWfwjSYQnu/LTat4LDGNvxaRd233lWMBzAROxC1Ta4/NtD/3H++zjjEmDLQCFUNsS8F6AopETSwdNReLzLUFevqkggKU2ze4XCzOOF4YY9h3uo09p9r6jAfDETqCYcoLvQPWKfTlnrtzNARiE0yTb/XqkOfJzZjGyM+ExfeBrwHGfv8W8GeAJFjWDDHOKNfpg4isx3J9MXv27KH2O2dwagY5rpuiuNz6qrTtVWpp6+7pE88AYq4UFY3heePQGb7/0iH217dT3xakOM/Du3ddix1OjFlr/bOnIDdjZKOhO2S5l8aScmuMiZ3zXGBUloYxpt4YEzHGRIEfYcUcwLIGZsUtWgOcssdrEoz3WUdEPEApljtssG0l2p8HjDGrjDGrKisrR3NIWUdvHwTb0sjBekFtgZ4+mVOgojESfrbpONuONXPZgqlcs3ga7YFwn8Cscw7LE7qn3Oqeotc9NZqYRr7PTdRATyR34owwStGwYxQOHwWczKpngHV2RtQ8rID3ZmNMHdAuImvseMWtwNNx6ziZUTcDL9pxj+eAa0WkzHZ/XWuPKUB7oG86aoHTOCeHfuit3eE+czSg96m4uUtFYzjqWrq5aNYU7v34MtYunQ7A2Y7e89Zsi8aglkbIekqezIwl5dbvsbv35VjSxrDuKRF5DLgKmCoitVgZTVeJyDIsd9FR4NMAxpjdIrIB2AOEgc8YY5wzdgdWJlY+8Kz9AngQeFREDmJZGOvsbTWJyNeALfZydxtjkg3I5zz9S1oX5WCLzrbuHs6rLu4zNsV2V6mlMTx1rQEuWzAV6Cu2s8oLACtzChg0eyoSNQTD0VHdMHMFJ+V2tAULwepz399izmaGFQ1jzCcSDD84xPL3APckGN8KnJ9gPADcMsi2HgIeGm4fJyOORVGc64Hwfj82j9vFlAKvisYwhCNR6tsCzJhilZlJ5NaLWRoJ3FPOddURDE9u0XAC4aOc3Ae514hJZ4RnKY57Kn5GOPQW+ct2IlFDeyA8IHsKLB+8isbQNLQHiRqoLs0HEotGU6d1DU0pSJw9Bbn1EDIausYwuS9maeSYe0pFI0sZEAjPMUujI+AUKxxoDJcV+jSmMQx1rd0AVDuWRoJU5eauECV5HrzugbcB53pyrrPJSqAngkhvfGIkOM3DHBdXrqCikaX0F40CrxMIz40LNFEJEYeyAl/sKVlJzKkWq+LxDNvSKMn34HZJH7G1Cl4OdE0BlNnWx2SfRNkeCFPo84wqZbbXPZUbv0kHFY0spSMYptDnxu2yLmaXSyj05U69oNYEfR4cygu9NHXmVu+QVNPf0hARW2zj3VOJ604BVBTZlskkt+jOdASpLPaPat0C+4GuSy0NJRNoD/TE4hkOuTQhK1GFW4fyQj/NnT2TPh10KE61BCjye/qcP0ts+1oaiTKnoDc43jzJY0eN7UEqi0YnGkU5mAYPKhpZS0cw3KejHeRWvaCh3FPlhV5CkSidOfYEl0rqWrupLs3rM1Ze6KM5zq03WN0psM67iKY2N47B0sjVDogqGlmK1Uujr2gU5FA561jXvvwEgXB9Ch6WutYAM6bk9xkrL/Rx1nbrBXoiNHYEmVqcWDQ8bhel+Zra3Ng+dtHIld+kg4pGltIWCMeC4A5WkbncePp2LI3EMQ0tJTIcp1p652g4xHd3fPPwWULhKJfOrxh0G+WFvkkd0wj0RGgPhEcvGj61NJQMoiPBxLdcck+1dvcgAkUJqouWqWgMSTAc4UxHMDZHw6Gi0EdLV4hI1PD8nnoKfW4uXTCEaBT4JrU119huWWWjjWm4XUK+N3esfwcVjSwlkXvKqheUGxdoW7clii7XwFTHChWNIYml2/ZzT5UV+ogaaOkK8fzeet53TiV+z+CT1soKJ/ckysYOWzRGaWmA08wqN6x/BxWNLKUjmMA95c8l99TAYoUOamkMze5TrQAsnt63bpfj1nvlQCP1bUE+cN7QRfQn+8z7mKUxBtEozsudjEYHFY0sJByJ0hWKDEi5LcqhQPiZjmDCkt0AxX4PLumdy6H0ZWdtKz63i3OqEovGw68fxe0Srl48bcjtlBdZM+8na2pzKkSjMAdLzKtoZCGONdE/5bbQ74l19Mt26tsCVJXkJfxORCjO88bqbyl92XmylcXVxfj6lb5wss521LZy47IZg87RcCgv8NETMTl300uWxvYgIr3u0NFQ6MudOKODikYW0hYrVjhwngbkRtHC060BppcmFg2wUnHbJnldpEQYY9h1spXzZ5YO+M4RCZfA37x/0bDbmuxuwMaOIBWFvj496kdKUQ5NuHVQ0chCYq1eE8Q0IPvzwrtDEdoC4UEtDYBivzc2a1zp5XhTF22BMBcMIhoel3DTspnMm1o47LYme8JBY3uQqaPMnHLIpSoNDqPtEa6kkd5eGgPLiED2i8bpNiv7Z/oQolGS75n0FVgTsfOkFQRPJBp5XjdPfHoN504vSWpbk71L4lgm9jlo9pSSETg/4v59EAp9uVHp9nSrLRpDuafyvDE3ndLLYEFwh5Vzygdk3Q1Gbzn1yXmeUyEauZSc4qCikYUM1ts5VyyNhnZLNKpKBv/BWoHw7D7O8eBQYwfzKwsHBMFHQ1mh01p38lUUNsaMqe6UQ5HfS3dPhHAkd7r3qWhkIWdt0eif1VGUIwXSHEtjqJhGSb5HYxoJqG8LDnneRkKR34PP7ZqUlkZbd5hQODrq2eAOhXal21wqrqmikYU0d4Yo8LkH9G4ujNXvz3LRaAtQ6HMPiNnEU5znpSMUJpoD6cWppKE9wLQxPh07iAhlhd5JWUqktqULgJqy/GGWHJqiHLH+4xlWNETkIRFpEJFdcWPlIrJRRA7Y72Vx390pIgdF5D0RWRs3vlJEdtrf3Sd2KywR8YvIE/b4JhGZG7fObfbfOCAit6XqoLOdwZrnFPpzI6ZR3xagaoh4BlhtYI2B9hz6MY6VSNRwpiPEtCHceiOlrMDwK3syAAAgAElEQVQXs2wnE7XNVhOrmrKCMW0nV1zG8SRjaTwMXNdv7MvAC8aYRcAL9r8RkSXAOmCpvc79IuI8Dn8fWA8ssl/ONm8Hmo0xC4F7gW/Y2yoH7gIuAVYDd8WL02SmqSsU66wWT6481ZxuDQyZOQW9zZnURdVLU6dVjHBacWrcU2D34JiE2VO9opEaSyPbXcbxDCsaxphXgKZ+wzcCj9ifHwFuiht/3BgTNMYcAQ4Cq0WkGigxxrxprJoEP+m3jrOtJ4FrbCtkLbDRGNNkjGkGNjJQvCYlg1ka+V43Lsl+0ahvCw4vGnZdKg2G95JMAsFIsRo3TUbR6KLI70nYBGwk9Foa2W39xzPamEaVMaYOwH53itjMBE7ELVdrj820P/cf77OOMSYMtAIVQ2xr0tPUGUpY2kBEsr5sQTRqaGhPxj1lWxqadhujIVYrKbWWxmTsqXGiqZuasnxsL/qo6XUZ5851mupAeKIzbIYYH+06ff+oyHoR2SoiWxsbG5Pa0WymqTM0IN3WIdtnoDZ1heiJGKqGCeY6QfJElsabh85y9y/3EAznztNdMjS2WaKRqkA4WDGNlq6enEoZTYba5q4xu6bAqlwA2R9njGe0olFvu5yw3xvs8VpgVtxyNcApe7wmwXifdUTEA5RiucMG29YAjDEPGGNWGWNWVVZWjvKQsoNAT4SuUGTQYnOFfndWm8LJTOyDXvdU/5jGD14+xCd//BYPvX6EF/Y2JFo1Z3HcU2OdWxCPc521TKLYkTGGk83dYw6CQ1zKbRY/yPVntKLxDOBkM90GPB03vs7OiJqHFfDebLuw2kVkjR2vuLXfOs62bgZetOMezwHXikiZHQC/1h6b1Dh1gAYXjex2T/U2vhlaNHotjd6bmTGG+144wBULpzKt2M/T20+O345mIA3tQUrzvQNSsceCc51NprhGW3eY9mA4JZZGYQ4GwoetJyAijwFXAVNFpBYro+nrwAYRuR04DtwCYIzZLSIbgD1AGPiMMcZ57L0DKxMrH3jWfgE8CDwqIgexLIx19raaRORrwBZ7ubuNMf0D8pOOYUXDl93uqTNJtth0KvzGV7o92xmiKxTh/YunsWhaMT996xit3T1jDmZmCw1twZS6pqD3OjvbGWL4uriZw8v7GzlQ3051aT43XDB9RLGJE83OHI2xWxp+jwuPS7L6N9mfYUXDGPOJQb66ZpDl7wHuSTC+FTg/wXgAW3QSfPcQ8NBw+ziZSMbSqLUv+mzEsTSmFg/dw8DrdpHvdfdxT51oso57VlkBK2aX8dDrR3hu12n+6OJZg20mp2hoD6R0jgb09uDINkvjCxt2cMa+lh5fv4Y18wfvhd6fVKXbgp2ckuXWf390RniWMZxoFPnddGVxyYLG9iBFfg8FvuGL6vWvdHvC/rHPKi/gwppSZpXns3Fv/bjta6ZR3xZM6RwN6L3OsimDKho1NHUG+ZPL5lKc5+GxzcdHtL7z0DUrBZYGWHM1VDSUtBETjUFaoWZ79tSZjhBTE0xcTET/SreOpeGkSi6fVcZuu1R4rmOMobE99e4pp2hhNlkabYEeosZ6ePjo8pk8u+v0iPb/RJM1R2OwHvUjpTDHKt2qaGQZTZ0hXMKgfvpsf6ppbA8knf1TnNfX0qht7qKi0BcLPi6ZUcKp1kBW3fBGS2t3D6FINKWZUwB+j5sivyerSok0d1kPEmUFXtZdPJtQOMrPt9UOs1Yv20+0sKS6ZMxzNBys7n3Za/33R0Ujy2jqsmaDu1yJL+hCv4dgOJq1efUj6WFQkt/f0uimprzXpbB0htVsaPepttTuZAZS78zRSFGF23iyrWihU/akrMDHkhklrJg9hYffOJrUb6IzGGbXqTYunpe6ikXFOdb7RUUjy2jqCA0az4DsL1twpiOUdDnq4ry+LV9PNHcxKy54uXSG1b1uT13uu6h2n7KO8ZyqopRvu7zAR1NX9tz0Wvo1KbvjqoXUNnfzy3cTTvPqwzvHW4hEDavnJR84H46qEj/1djfKXEBFI8to6hp8NjhYgXCAjiwsjx4MR2jt7km6L3NJnHsqEjWcaulmVpylUV7oo7o0b1JYGtuON1Pk97BoWuKOfWMh2+pPNdv9P5yHq2sWT+PcqmLu//2hYUvpbz7ahEtgxewpKduf6aX5NLQH6clS678/KhpZRnNnaNAgOGR3KeYzHdaNKfmYhmX2G2M43RagJ2IGZLwsnVEyOUTjWAvLZk3BPYjbciyUFfpiCRjZQG87ZOt34nIJf3nVfA40dLDtePOQ62450sSSGSVD9nIZKdWleRjTWxss21HRyDJauntiGS2JKPRl7wzUxljBveREY05FAT0Rw/76jt45GuV9c+uXzCjlcGMH3aEIgZ4IP3j5EKdaulO742mmMxhm3+m2lD4dx1NekH2i4XYJJXm92U/vW2SVF3r72OCiEQpH2Xa8mYvnlqd0f5ySOKdbc+O6U9HIIowxtHSFKM3PUUvDFo1k3VPvX2wVV35+bz3HzybOrV86o4Sogb2n23hhbwNff3Yfa7/zCj9/uxarWk32s6O2haiBFXPGp91MWaGP7p4I3Vky/6e5q4cp+d4+2U8VRX7mTS1k6xCicaixg2A4yvLZqT2P1bZo1LXmRlxDRSOL6AxF6IkYygqGsDSyuEBab92p5ESjqiSPi2pK+d2eev777RNML8kbMIs3PoNqb10bbpeweHoxX/jvHXz60bczvjVuMsL2zvEWAJbPGh/RqMiyCX4tXaFYEDyeFbPL2HasedBzety2VudWpGZSn0N1iXVNnlbRUCaa/lkhiSjK4uwpxz2VqCvhYHzgvCp2nGhhy9Fm7rhqAR5330t65pR8SvO97DnVyt66NhZUFvL4+ku58/rF/G5PPRu2nBhky+nHGMNH73+Dz2/YPmQA953jLSyoLKR0iOtiLJRlWdHC5s6ehBmGK+eUcbYzxLGzicvsOC7O2eWpFY2SfA8FPrdaGsrE02KnPU5JJhCe4U/QiTjTYVVp9XuSr9L6waVVgJXW+PEENaZEhCXVJew51caeujaWVJfgdgmf/oMFLJ5ezK931qVs/1PNvtPtbD/RwlPbTvLN37036HLHmzpZUJn6VFuHmVOsJ+XDZzrH7W+kkuauUMLfyErbfTdYXON4UxcleZ4hf1+jQUSYXpqnloaSOlq6Qjy/p55Az9DWQUtsputQKbfZHQgf6Yzmc6uK+chFM/jHDy8ZtCT40hkl7Klro641wHnVJbHxD11QzZajzRn7Y37Brpt1/fnTuf+lQxxJcNN2ej/MTEFxvcFYPL2YQp+bLUeyo8h0c1cooQt30bQiivM8g8Y1jjd1MTvFrimH6tI86jQQrqSKH716mD//yVYu//qLvLJ/8M6DzUm4p/weF+4sLcVc3xYYce0kEeHfP7GcD184Y9Blls4soSdiuXeWzOgVjRsurAbgNxlqbWzc28BFs6bwpesWA/D6wTMDlmnt7qEzFIlZA+OBx+1ixZwyNmeBaBhjaO7qSfhg5XIJS2eU8N7pxCnYx5u6Uu6acphekp+xDycjRUUjA9hzqo2ZU/Ipyffy1Wd2ExnEf51MTMPqE56d3ftOtnSPy83PmRkO9LE0FlQWsXh6Mb/ddTrlf3OsNLQH2HGihQ8snsacigKml+Tx1uGzA5ZzyniPp2gAXDKvnPfq22PXYKbS3RMhFI4OOgF2VlkBJxOkXEejhtqmvpNDU0l1aR717cFBf9vZhIpGBrC/voOL55bxxbXncvhMJ8/uSvzkG4tpDJFyC9lZtDAUjtLQHmTGONz85k8txO9xMa3YPyCdd838Cnafas249NuX37MszmvOq0JEuHRBBW8dbhqwn86ck/F0TwGxuQtbjg49OS7dxBcrTERNWQH1bcEB/eNPtwUIRaLjZ2mU5hGJmliPj2xGRSPNtAV6ONnSzTnTi7lu6XQWVBbyvRcPJryJNXf1UOhz4/MM/d9WkIXl0evbAhgzPk/MHreLVXPLEk7aWlBZSGcoEiv4lym8W9tKsd/D4ulWWZA188s50xHkUGNHn+Wcp+bxtjQumjUFn9vFlqOZ7aJyMrwGC2Y7KdmnWvq6io6PU+aUgzNX42iWJBMMhYpGmjlQ3w5YAV2XS/jUmjnsO92e0IRuGSQrpD/Z2CnMOd7xsDQAfnTrKr71RxcNGF8wzco6OtjQMeC7dLKnro3zqkti1YwvnT8VgDcP9XVRnWzuJs/rGrKIZSrI87pZNmtKQhdZJhFf4TYRjmj072453qIxd2ohAJ/88Sbu/uWecfkbE4WKRprZd9oWDfuJ8hz7/eiZgbnkLd09Q8YzHIqysOnLyebxdbMU+DwJs6sW2qmq/Z/g00k0athb19YnaD+rPJ9Z5fl85/kDbNjaO7fkZEs3M6bkp6z3w1BctrCCnSdbMzquMax7yhYFJxbkcKKpC5eM30PLgsoifvPZK7l2SRX/+caRrG7JrKKRZvafbqfQ5465F+ZPtW5iR84ONGOb7V4aw1Hoy76mL45v3jHjJ4rKYj/Ffk9GWRrHmrroCkX6iIaI8MCnVjGnooC/f/Jd3q21ZoGfGqfkgURcuagSY+D1g5lrbQznnppekofHJQNu2ocbO5kxJR+ve/xuiUtmlPB/PnQeAD9/++S4/Z3xZkxnSESOishOEdkuIlvtsXIR2SgiB+z3srjl7xSRgyLynoisjRtfaW/noIjcJ/Zjk4j4ReQJe3yTiMwdy/5mEq8dOMP/+cVOth1v4ZzpxbEnxaoSP/led0LfZ2tXT1KzfrMxEH6qtZupRb5B51qMFyLCgmlFGWVp7LGr8i6Jy/QCK/PrB3+8EiCW/nqypXtA6ZTx4qKaUorzPLx6YPC08FTR0Bbgj374ZsK5KUNR1xrA65ZY6ZP+uF3CjCn5fSyNSNTwxqEzrE5xocJE1JQVcPmCqfz32yeGLdOeqaRCVq82xiwzxqyy//1l4AVjzCLgBfvfiMgSYB2wFLgOuF9EnDvE94H1wCL7dZ09fjvQbIxZCNwLfCMF+5t2jDH886/38LNNx9l5sjUW7ATrJjanoiDhj2WwSUv9qSrNo74tkFXd+2qbJ+6JuT8LKjNMNOpa8biERQkaKk0ryWPmlHy2HW8m0BPhTEeIGaUTc948bheXLajg1QNnxj3b7IktJ9h8pImX3msY0Xp1rd1Ul+YP2tkSrLhGvGhsP9FMc1cPV9kFMMebW1bVUNvcnfHxocEYD1vsRuAR+/MjwE1x448bY4LGmCPAQWC1iFQDJcaYN411Jf6k3zrOtp4ErnGskGzm7WPN7Dvdzp9ePpeL55Zx/fnVfb6fX1k4wNKIRg2t3YknLfVnXkUh4agZkCGSyZyyffPpYMG0QurbghnTknPPqTYWTisatJzKyjllbDvW0ps5NUGWBlguqpMt3SO2AEaCMYYn7Z7e++raR7SudR0N7eK0RKPXPfXivgbcLuEP7PLp480Hl1ThEngrCyZLJmKsomGA34nI2yKy3h6rMsbUAdjvjnzPBOKrw9XaYzPtz/3H+6xjjAkDrcCAPowisl5EtorI1sbG8Tedx8qjbx2j2O/hi2vP5b//8jLed07fi3VuRSHHm7r6WArtgTBRA6X5w1saTqZGorhIOmnqDPHDlw8NEERjLIFLl2g4wfDDjek/X8YYdp/qGwTvz4rZUzjdFuAp+8Y6z/7/ngicvhSvHhg4Oz1VbD7SxLGzXfjcLvYNMnt7ME61BIa1vPrP1XhxXyMrZ5eNW8HH/hT4rA6LTlwq2xiraFxujFkBXA98RkTeN8SyiSwEM8T4UOv0HTDmAWPMKmPMqsrKiXlaGC1tgR6e3Xmaj62socDnSbjM3KmWpRBvQg+XSthnfbt+TqblhN//+4P8y7P7uOqbL3HX07tiPt2Wrh66eyJpE41FVZZ7cF9d+jv81TZ309Ae5KKawRsqOX0z7n/pEKvmlLFs1vg0X0rE7IoCZpcXpDyu8cgbR3n49SMYY3j4jaMU+T18bOVM3qtvT3oWdTgS5XTb8A8f8XM1TrcG2FvXxtUT5JpyuLCmlJ21mTepNBnGJBrGmFP2ewPwC2A1UG+7nLDfHadkLRBfhrQGOGWP1yQY77OOiHiAUiA7bTqbHSdaCEWifOC8qkGXmZ/AUkim7pRDZbGfQp97XF0IIyXQE+Hn22q56txKPrVmDo+8eYx/fHqXVXBvgiaoDcbcigIqCn0ZUVtpk70Pl8wfPCh7XnUJeV4XxsCXrl88Iem28Vy5aCpvHjqbsp7XDW0B/vnXe/jqL/fwJ/+5hWd3neYvrpzP8tllBHqiHEvSYm6wy3QMJxpz7IeqA/XtbDpixRWuXDR1bAcxQi6cNYWznaGE87EynVGLhogUikix8xm4FtgFPAPcZi92G/C0/fkZYJ2dETUPK+C92XZhtYvIGjtecWu/dZxt3Qy8aLJRmuPYfrwFEbhwVumgyzjupXhLoaV7+LLoDlYwvTDpH9tE8Nzu0zR39XD7FfO4+8alrH/ffH626TjbjjfHsoUWVE6cmyUeEWHN/AreOnw27U9+m4+cZUqBl3OmFQ+6jNft4przqvjwhdUpb02aDFcuqqQzFIk1fxorP910nHDUsHpeOS/vb+SGC6bzN+9fyHnTLRedM5dpOJwqssPFNM6fWYrf4+Ktw01sOdpEUdzM+4niwpnW7//d2tYJ/bupILF/JDmqgF/YTzke4L+MMb8VkS3ABhG5HTgO3AJgjNktIhuAPUAY+IwxxplMcAfwMJAPPGu/AB4EHhWRg1gWxrox7G9GsP1ECwsqiygZonF9RaGP4jxPn4wep0FMsvMY5k4tYO8Ig4jjRVcozEOvH2VWeT6XL5iKiPCZqxbyo1cP8+qBMxw/20VFoY+F08avJ8RwrJlfzq931nGiqXvcymMnw6YjTVw8t3zI7B+A//jkirQJ3KULKnC7hFcPNLJ63thEKxiO8F+bjnH1udP4j0+u4NlddVx/fjUuO3vMJZbb8IYLqofd1kk78WM4S8PvcbNyThlvHT5L1BiWz54yoHnXeLO4uhivW3i3tjWpY8skRi0axpjDwIC6DMaYs8A1g6xzD3BPgvGtwPkJxgPYopMLGGPYfqJlWP+piHBRzRTePtb7JLe/vp1ivyd50ago5He76wlHohP+g4inrrWb2x7azIGGDv7t5otiN8PSAi8Xzizl9YNnONnczZr5FRPuZolnzXwrv+Ktw2fTJhqnWwMcO9vFp9bMSWr5dJ2v0nwvF9WU8sLeBj7/wXPGtB9PbTvJmY4Qf3r5XPJ9bv5wRa+nOs/rZt7UQvYmaWmMZILomvkV3Pv8foC03LT9HjfnVZew40T2BcN1RvgEUtvczdnOUFKBy9Xzytl3uo1WuyzC/voOFlUVJf0DTRRMTwePvnmMw42d/OTPVnPzypo+3122cCpbjzVzqjXAmiF8+BPBwmlFVBT60po77/jXL5k3IEEw4/jDFTXsqWvjzTGcr1A4yvdePMhFs6ZwxcLEMYWLZk3htQNn2H1qeDfOqZZuSvI8FA9hxTusmV+BMWAMrJo7Pr3Vh+PSBRW8efgs//g/u4ZtwJZJqGhMIO/YTxXJioYxsOWoVQ77QH0751Ql73d10jCPjjKuUd8WSEkZ5y1Hmzh/ZilXJsiBv2LhVBwPi/Okny6cuMabaYxr/PrdOioKfZxXPbH+9dFw88oaKov93P/7Q6PexoatJzjZ0j2ktfLl6xZTmu/lLx7ZOuz1OJK07YtmlZLntRqWTWT2WTxf+OC5/MWV83j0rWN87PtvZFTiylCoaIwT3aEI3/7dezyz4xRHz3Sy/UQLj7xxlDyvK1accCiW2aWoNx9t4kxHiOaunlhqaDI4orFnFGmk0ajhEw+8xec37BjxuvEEeiLsONE6qN975ZwyfB5X2uMZDpcvnEpdayAts8NPtwZ4YV8DN6+qSas7MVnyvG7+4sp5vHbwDJtGYW1s3FPP15/dx8o5ZbxviMylaSV5/OjWVZxqDfDzt2sHXQ5GVofL73FzxcKprJpTNmjq+3jj87j4Px9awo9vXcXJlm5u+cGbA/p8ZCKZf3VmKT94+RD3vXiQzz72Dld98yVu+o/X2X+6nX/40JKkiqI5pag3HWnqUz49WaYW+Vk6o4Tn99SPeN/fOHSWw2c6eedY85jq47xb20ooEh00wyfP6+ZjK2q4eVVNWuMZDu87x7p5vbx//CauDcYTW04QiRo+uXr2hP/t0fK/LpnDrPJ8Pvv4OzS0JV994PWDZ/iLn2xlfmUh//6J5cP+319QU8rCaUVDusLCkSgnmruoHiZzKp77PrGcB//k4qSXHy8+sKSK765bzpmOIC/sHVnZlHSgojEOnG4N8MNXDvGhC6r5+R2X8a1bLuI7H1/Gq1+6mj9OMsgJlotq18nWWO7+OQlqEQ3F2qXT2Xa8ZUQ/aID/2nwMgPZgmMNjMJmdhj2r5gzuM/6XP7yAO68/b9R/I5XUlBUwv7JwyD7t40FPJMrjW45z5aKpzKlIT9rxaCj0e3jgU6to6w5z60Ob2XUyufTRF/Y24Pe4eGL9pUm7ky6dX8HmI02Dzg15fm8D7YHwoLGRRBT4PBT502Nl9OeKhVOZXpI3rDWVCahojAP3btxPNApfvn4xK+eU8bGVNdy0fGZScyzi+fBF1bgEvvf7g5Tme6ks9g+/Uhxrl04H4HcjsDYa2gP8bnc9f2CXNhlLdsfmI00smlY0aL/mTOR9iyrZdOTshAYmn9hygrrWAH92+bwJ+5up4rzqEu7/4xWc6Qjyke+9xuc3bB92ftC2481cVDOFfF/yFY0vW1BBVygy6LyGR944yswp+UNOms1k3C7hpuUzeWl/I43tmdVFsj8qGinmTEeQX7xzko9fPGvMTeoXTy/hS9ctJhI1nDOCzCmHc6qKmFtRwHO7Tye9zn9vrSUcNfzjh8+j0Odmxyjr4wR6Irx9rJmLx5jHP9H8wTmVBHqiPPT6kQkpYNgdinDfCwe4eG4ZV52b2SVwBuPqc6fxwheu4vbL5/Hrd+v44L2v8OBrRxK6NgM9EXafamX5nJEFny+xEyXePNTXdWiMVdb8zcNn+eM1c7IiHjQYH1sxk0jU8D/vjK7XxusHz0xIVYPsPcMZyuObjxOKRLntsrkp2d6fXT6PT62Zwy0rZw2/cD9EhGuXTuetw2fpCg3fXyMaNTy+5Thr5pezcFoxF9SUsmOUM1Z/v6+BjmCY68+fPqr108Wa+RUsnl7Mv/72PW747quEwuNbXv5nm47R0B7ki2snvhxIKinN9/IPH17CK39/Ne9bVMnXfrWH79jzIOLZdbKVnohh5eyRpbmWF/o4r7qE1w72ikagJ8KH//01PvmjTZTme1l38ch/I5nEoqpiVs4p4782Hx9xLPFsR5DPPb6d/xtX0228UNFIIeFIlJ++ZfmmU5UN5HIJX7vpfP5olD+IyxZU0BMxvH2sedhlXz14hhNN3XzyEivuctGsKew91TaqjI6n3jnJtGI/ly2Y2Jo+YyXf5+Y3n72Sb95yEbXN3by4b3wDk0++XcuK2VPGPLM6U6gqyeNHt67k/7toBj945XCskoGDcx2uGCLONRjXLZ3OW4ebeOi1IwA8s+MUu0+18cW157Lx796XVW7QwfjUmjkcOdPJ64cGJmMYY2Kuq5auEN95fj/7TrcR6Ilw51M7aevu4dt/tGzYagJjRUVjFNS1dvODlw8NuJm+vL+R022BpGf0TgQXzy3H4xLeODR8WuSjbx6lvNDH2qWWX3hZzRRCkSi/3FE3or/Z3BnipfcauHHZDNzjfAGPBy6XcNOyGUwr9vPkOAYm951uY9/pdm5aPnP4hbMIEeHOGxbjEvj6s/v6fLfteDOzywuYWjSy+BzAZ65ewHVLp3P3r/bw2Obj/OfrRzm3qpi/umoB00omtk3weHH9BdMpL/Tx07eODfjuG799jzX/8gLbjjdz/0uH+M7zB7juO69ywVef43d76vnfa88ZsqR+qsiM1IEsoicS5a9+to13jrdw7GwX//KHF8S+++2u0xTnebjq3IktszwUhX4PF82awpvDiMZT22p5fm8DX/jgObHmP1cvnsbFc8v4+yd34Pe4+MhFM5L6m7/aWUdPxGT1zdDjdvHRFTP58atHaGwPjjgJIRme2X4Kt0uyrvZQMlSX5rP+yvnc9+JBvnimkzkVBWzYeoKX3mvkwxcmdx31x+N28d1PLOPTj77NnU/tBODrf3hBVrv1+uP3uPn4xbP44cuH2HGihYvsiYcb99Tzg5etiZR3/3IPB+rbuXZJFctmT6GtO8wl88onLCamlsYI6AiG+f9/s5d3jrdw2YIKHtt8nCe2HAcs19TGvfV84LwqfJ7MOq2Xzq9g58nWQfuG76xt5R/+Zxer55Vzx1ULYuN5XjcP/+lqls2awj8+vSupuAjAL7bVcm5V8YAe19nGzStqiEQNP3x54KznaNSMaeZ4NGp4evsprlg4dVRP3dnAJy+Zg0vg59tqefiNo3zp5ztZMbuML11/7qi36fe4+eGnVnL9+dOpKcvP6geTwbjjqgVMK87j7598l1A4SkN7gC8+uYOlM0r4yg2L2X6ihc5QhL/74Dn81VUL+fL1i7l68bQJE8/MurtlME9tq+Xif36e/3z9KOsunsWjt1/CZQsq+Nqv9nKqpZvNR5po6eqJpblmEpctqCASNdzz6z0xkXP45Y5T3PLDN5iS7+W+dcsHZJ8U+j185YbzaOnqYcOWEwzHsbOdbDvewk3LZ2b9E+CiqmI+sXo2P37tCA+/fiTWSfH5PfVc8i8v8JVf7Bz1tp/cVmvNAl5VM/zCWcr00jyuWFTJhq0n+PbG/Vy5aCo//fNLmFY8NleS3+Pm+3+8kpf+91XkeZNP280WSvK83PPR83mvvp2/27Cdrzy1i65QhO+uW8Ztl81lbkUBH1xSxXlpeihT91QSvLivni8++S6r5pHJFd8AAApNSURBVJTx99edy4rZZYgI3/jYhVx77yv8/ZPv4ve4yPO6YvMbMokVc8qoKPTx2OYTPLb5BIuqilkxu4xnd9bx2cff4eI55dz/xysGfeJdNbeclXPKePD1I8OmNf7inZOIwI3LRueCyDTuvnEpp1q6+eov9/CtjfvxuITmrh5K8708tvkEH1xSxfsXJzc3wBjDtuPN5Hs9/Otv97Fi9hRuOD/3XFPx3LKyhr957B3cLuH/fnhJSmNc2ZxeOxzXnFfFF9eey7c37icSNXz5+sUstHus/OqzV+JJY6xQ0t10JtWsWrXKbN26NWXbO9TYwUf+/TXmVxby+PpLB8wg/fGrh/nnX+8FrCJu37xlQLX4jCAYjtAVjHD1t15i1ZwybrtsLrc/vJXzZ5bwsz9fM+xEq9/tPs36R9/mA+dN459vuoDpCcpPR6OG93/rJapL83ls/ZrxOpQJJxiOsHFPPa8fPItL4JyqYj62soaP3f8GLd0hfv3ZK4d1MXUGw3zp5+/yq3etpAIR+OVfX8H5MwdvxpULBHoivO9ff8+HLqzmro8sTffuZB3vHG/mjUNn+fT75o+7SIrI28aYVcMup6IxOMFwhI/+xxvUtXbzm89dSXWChvXGGPbWtTO12EdlkT/jXTLfff4A9z6/H69bWFBZxOPr1yQ1U90Yw4OvHeHfnnuPqDFcu3Q6X/3I0j4B4qe3n+Rzj2/nu+uWceOy3PM192fXyVZu/sEbnFtVzGPr1wxa+O5AfTt3/Gwbhxs7+NsPnENZoY+SPM+kOEdgCYfP7Rr3VFBlbKhojJHa5i6++N/v8ubhs/z41lV8YEl2lifoT0tXiKu/+RILpxXx49supjR/+N4D8Rw728nDbxzlZ5uOs3bpdP79E8sBS2Df/82XKc338qu/uWLS3CCe31PP+ke3MqeikD+9fC4fW1FDYZw1+ot3avnKU7so9Lu5b91yLhtBbSRFmUhUNEZINGrY39DO1CI/L+5t4O5f7cEYw10fWTrqiXWZSmt3D8V+z5hu7N/euJ/7XjjAk395KefPLOWffrmbxzaf4Ke3X8IVQ5S6zkVe3FfPd58/wI7aVorzPKywZzs3d4V4t9YqDf+9TyzPmbkESm6iojFCmjpDrPjaxti/18wv599uvmjM9aNyla5QmPd/82U6g2F8HhdnO0P8+RXz+IcPL0n3rqUFYwzvnGjhJ28c5ciZTgxQnOdhzbwK7rhqQU4HbZXcQEVjhHSHIry4r4EzHUGmFHj5yIUzJo2LZbS8c7yZn751nEg0ys0rZ006C0NRcolkRSMrUm5F5Drgu4Ab+LEx5uup/hv5PjcfujC30x9TzfLZZSwfYeE5RVGym4y3mUXEDfwHcD2wBPiEiExOH4iiKEqayXjRAFYDB40xh40xIeBx4MY075OiKMqkJBtEYyYQX7+i1h6LISLrRWSriGxtbJzYVp2KoiiTiWwQjUTR6D7Re2PMA8aYVcaYVZWVmVfGQ1EUJVfIBtGoBeInStQAp9K0L4qiKJOabBCNLcAiEZknIj5gHfBMmvdJURRlUpLxKbfGmLCI/DXwHFbK7UPGmN1p3i1FUZRJScaLBoAx5jfAb9K9H4qiKJOdnJsRLiLtwHtAKdCawk2nentTgYHd40dPKvdPz11ubw/0HI6VXDt/U4FCY8zwmUTGmJx6AVvt9wdSvN1Ub29rpu6fnrvc3p6eQz1/YzmebAiEj5ZfZvj2Uk0q90/PXW5vbzzI9GPO9HOYNcebi+6prSaJolvpJlv2MxPRczd29ByOjVw7fyM5nly0NB5I9w4kSbbsZyai527s6DkcG7l2/pI+npyzNBRFUZTxIxctDUVRFGWcUNFIESIyS0R+LyJ7RWS3iHzOHi8XkY0icsB+L7PHPygib4vITvv9/fZ4gYj8WkT22dtJee+QTCNV587+7rcissPezg/s0vo5TyrPYdw2nxGRXRN9LOkgxdfgSyLynohst1/T0nVc40Iq07wm8wuoBlbYn4uB/Vj9P/4V+LI9/mXgG/bn5cAM+/P5wEn7cwFwtf3ZB7wKXJ/u48uGc2f/u8R+F+DnwLp0H1+2nUN77A+B/wJ2pfvYsu38AS8Bq9J9TON2rtK9A7n6Ap4GPog10bDaHqsG3kuwrABnAX+C774L/EW6jyfbzh3gxUo7/Hi6jyfbziFQBLxm3zQnhWik+PzltGioe2ocEJG5WE8im4AqY0wdgP2eyFT9GPCOMSbYbztTgI8AL4zn/mYSqTh3IvIc0AC0A0+O8y5nHCk4h18DvgV0jfvOZiAp+v3+p+2a+kcRSdTeIWtR0UgxIlKE5Rb5W2NMWxLLLwW+AXy637gHeAy4zxhzeDz2NdNI1bkzxqzFeir0AwN89bnMWM+hiCwDFhpjfjGuO5qhpOga/F/GmAuAK+3Xp8ZjX9OFikYKEREv1gX3M2PMU/ZwvYhU299XYz0BO8vXAL8AbjXGHOq3uQeAA8aY74z/nqefFJ87jDEBrBL6k6Y1cIrO4aXAShE5iuWiOkdEXpqYI0gvqboGjTEn7fd2rLjQ6ok5golBRSNF2Cbog8BeY8y34756BrjN/nwblq/UcT39GrjTGPN6v239M1bBsb8d7/3OBFJ17kSkKO4H7gFuAPaN/xGkn1SdQ2PM940xM4wxc4ErgP3GmKvG/wjSSwqvQY+ITLU/e4EPAzmVgaaT+1KEiFyBlem0E4jaw1/B8otuAGYDx4FbjDFNIvIPwJ3AgbjNXIuVMXUC62bn+Ei/Z4z58bgfRJpI4bkT4FdYbik38CLwd8aY8EQcRzpJ1Tk0xsQ/Sc8FfmWMOX/cDyDNpPAa7ARewUrEcAPPA583xkQm4jgmAhUNRVEUJWnUPaUoiqIkjYqGoiiKkjQqGoqiKErSqGgoiqIoSaOioSiKoiSNioaiTDAi8pcicusIlp87WarNKpmPJ907oCiTCRHxGGN+8P/au3+QqsIwjuPfh3IJopZWVwlEkCYJaW53ESnTNWh1bGrLiFysQfHPIm1NRZMQ0S7SVEvDXQK7la1Pw3kvXA5m5xzMBL+f8eHwcN7lPjwvnN/93+8hdeXQkFoqH729pvrwa5IqRvsucB14QpUS+xW4l5m9EsPxHrgJvIqIy8DPzHxcsp5WqSLxPwGLmXkQETeANarQwHendzrpeF5PSd2MAS8ycwL4DtwHVoCZzBz84D8aev5qZt7KzOVan01gqfTZAx6W+jrwIDOn/uUhpLbcNKRuvgxlDm1TRU6MA29LEvYFoDf0/E69QURcoRomu6W0Abw8or4F3D75I0jtOTSkbur5Oz+A/WM2g8MWveOI/tKZ4PWU1M1oRAwGxCzwAbg2qEXESPmvhT/KzD5wEBHTpXQH2M3Mb0C/hOgBzJ3860vduGlI3XwE5iPiOVXS6QrwBnhWrpcuAk+B/b/0mQdWI+IS8BlYKPUFYC0ifpW+0plgyq3U0nmKDJfqvJ6SJDXmpiFJasxNQ5LUmENDktSYQ0OS1JhDQ5LUmENDktSYQ0OS1NhvvapaZxWDSQMAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -252,10 +3292,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 13, + "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", @@ -274,7 +3312,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +3336,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 15, + "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 +3375,59 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "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": 16, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,9 +3442,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": "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)" ] @@ -341,9 +3475,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": true - }, + "metadata": {}, "outputs": [], "source": [] } @@ -364,7 +3496,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4, -- 2.18.1