From 6b3d080f40d1243813c4cf796ebab2f9ecef178f Mon Sep 17 00:00:00 2001 From: e6cfd8a10624671909b6755e01761c6e Date: Thu, 8 May 2025 12:12:18 +0000 Subject: [PATCH] exercice 3 part 1 --- module3/exo1/analyse-syndrome-grippal.ipynb | 2216 ++++++++++++++++++- 1 file changed, 2176 insertions(+), 40 deletions(-) diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index dfd0b86..a5729b0 100644 --- a/module3/exo1/analyse-syndrome-grippal.ipynb +++ b/module3/exo1/analyse-syndrome-grippal.ipynb @@ -28,20 +28,11 @@ }, { "cell_type": "code", - "execution_count": 7, + "execution_count": 8, "metadata": {}, - "outputs": [ - { - "ename": "SyntaxError", - "evalue": "(unicode error) 'unicodeescape' codec can't decode bytes in position 2-3: truncated \\UXXXXXXXX escape (, line 1)", - "output_type": "error", - "traceback": [ - "\u001b[0;36m File \u001b[0;32m\"\"\u001b[0;36m, line \u001b[0;32m1\u001b[0m\n\u001b[0;31m data_url = pd.read_csv(\"C:\\Users\\DELL\\Downloads\\inc-3-PAY.csv\")\u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m (unicode error) 'unicodeescape' codec can't decode bytes in position 2-3: truncated \\UXXXXXXXX escape\n" - ] - } - ], + "outputs": [], "source": [ - "data_url = pd.read_csv(\"C:\\Users\\DELL\\Downloads\\inc-3-PAY.csv\")" + "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" ] }, { @@ -68,9 +59,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "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
020251832065815191.026125.03123.039.0FRFrance
120251732257217604.027540.03427.041.0FRFrance
220251632856422550.034578.04334.052.0FRFrance
320251533572129592.041850.05344.062.0FRFrance
420251433757931232.043926.05647.065.0FRFrance
520251333967333686.045660.05950.068.0FRFrance
620251235254345627.059459.07868.088.0FRFrance
720251135946952154.066784.08978.0100.0FRFrance
820251036033453048.067620.09079.0101.0FRFrance
920250938453174994.094068.0126112.0140.0FRFrance
102025083136020124824.0147216.0203186.0220.0FRFrance
112025073208952195988.0221916.0312293.0331.0FRFrance
122025063273519258159.0288879.0408385.0431.0FRFrance
132025053334395318416.0350374.0499475.0523.0FRFrance
142025043350043332885.0367201.0522496.0548.0FRFrance
152025033252772238917.0266627.0377356.0398.0FRFrance
162025023257247242991.0271503.0384363.0405.0FRFrance
172025013231549214627.0248471.0345320.0370.0FRFrance
182024523201726185870.0217582.0302278.0326.0FRFrance
192024513201697187843.0215551.0302281.0323.0FRFrance
202024503136694126369.0147019.0205190.0220.0FRFrance
21202449310848799037.0117937.0163149.0177.0FRFrance
2220244838738178687.096075.0131118.0144.0FRFrance
2320244737628667626.084946.0114101.0127.0FRFrance
2420244635639949006.063792.08574.096.0FRFrance
2520244534734740843.053851.07161.081.0FRFrance
2620244433603930122.041956.05445.063.0FRFrance
2720244334657239928.053216.07060.080.0FRFrance
2820244236778560009.075561.010290.0114.0FRFrance
2920244137943571386.087484.0119107.0131.0FRFrance
.................................
208419852132609619621.032571.04735.059.0FRFrance
208519852032789620885.034907.05138.064.0FRFrance
208619851934315432821.053487.07859.097.0FRFrance
208719851834055529935.051175.07455.093.0FRFrance
208819851733405324366.043740.06244.080.0FRFrance
208919851635036236451.064273.09166.0116.0FRFrance
209019851536388145538.082224.011683.0149.0FRFrance
20911985143134545114400.0154690.0244207.0281.0FRFrance
20921985133197206176080.0218332.0357319.0395.0FRFrance
20931985123245240223304.0267176.0445405.0485.0FRFrance
20941985113276205252399.0300011.0501458.0544.0FRFrance
20951985103353231326279.0380183.0640591.0689.0FRFrance
20961985093369895341109.0398681.0670618.0722.0FRFrance
20971985083389886359529.0420243.0707652.0762.0FRFrance
20981985073471852432599.0511105.0855784.0926.0FRFrance
20991985063565825518011.0613639.01026939.01113.0FRFrance
21001985053637302592795.0681809.011551074.01236.0FRFrance
21011985043424937390794.0459080.0770708.0832.0FRFrance
21021985033213901174689.0253113.0388317.0459.0FRFrance
210319850239758680949.0114223.0177147.0207.0FRFrance
210419850138548965918.0105060.0155120.0190.0FRFrance
210519845238483060602.0109058.0154110.0198.0FRFrance
2106198451310172680242.0123210.0185146.0224.0FRFrance
21071984503123680101401.0145959.0225184.0266.0FRFrance
2108198449310107381684.0120462.0184149.0219.0FRFrance
210919844837862060634.096606.0143110.0176.0FRFrance
211019844737202954274.089784.013199.0163.0FRFrance
211119844638733067686.0106974.0159123.0195.0FRFrance
21121984453135223101414.0169032.0246184.0308.0FRFrance
211319844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2114 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202518 3 20658 15191.0 26125.0 31 23.0 \n", + "1 202517 3 22572 17604.0 27540.0 34 27.0 \n", + "2 202516 3 28564 22550.0 34578.0 43 34.0 \n", + "3 202515 3 35721 29592.0 41850.0 53 44.0 \n", + "4 202514 3 37579 31232.0 43926.0 56 47.0 \n", + "5 202513 3 39673 33686.0 45660.0 59 50.0 \n", + "6 202512 3 52543 45627.0 59459.0 78 68.0 \n", + "7 202511 3 59469 52154.0 66784.0 89 78.0 \n", + "8 202510 3 60334 53048.0 67620.0 90 79.0 \n", + "9 202509 3 84531 74994.0 94068.0 126 112.0 \n", + "10 202508 3 136020 124824.0 147216.0 203 186.0 \n", + "11 202507 3 208952 195988.0 221916.0 312 293.0 \n", + "12 202506 3 273519 258159.0 288879.0 408 385.0 \n", + "13 202505 3 334395 318416.0 350374.0 499 475.0 \n", + "14 202504 3 350043 332885.0 367201.0 522 496.0 \n", + "15 202503 3 252772 238917.0 266627.0 377 356.0 \n", + "16 202502 3 257247 242991.0 271503.0 384 363.0 \n", + "17 202501 3 231549 214627.0 248471.0 345 320.0 \n", + "18 202452 3 201726 185870.0 217582.0 302 278.0 \n", + "19 202451 3 201697 187843.0 215551.0 302 281.0 \n", + "20 202450 3 136694 126369.0 147019.0 205 190.0 \n", + "21 202449 3 108487 99037.0 117937.0 163 149.0 \n", + "22 202448 3 87381 78687.0 96075.0 131 118.0 \n", + "23 202447 3 76286 67626.0 84946.0 114 101.0 \n", + "24 202446 3 56399 49006.0 63792.0 85 74.0 \n", + "25 202445 3 47347 40843.0 53851.0 71 61.0 \n", + "26 202444 3 36039 30122.0 41956.0 54 45.0 \n", + "27 202443 3 46572 39928.0 53216.0 70 60.0 \n", + "28 202442 3 67785 60009.0 75561.0 102 90.0 \n", + "29 202441 3 79435 71386.0 87484.0 119 107.0 \n", + "... ... ... ... ... ... ... ... \n", + "2084 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "2085 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "2086 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "2087 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "2088 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "2089 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2090 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2091 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2092 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2093 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2094 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2095 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2096 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2097 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2098 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2099 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2100 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2101 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2102 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2103 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2104 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2105 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2106 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2107 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2108 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2109 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2110 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2111 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2112 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2113 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 39.0 FR France \n", + "1 41.0 FR France \n", + "2 52.0 FR France \n", + "3 62.0 FR France \n", + "4 65.0 FR France \n", + "5 68.0 FR France \n", + "6 88.0 FR France \n", + "7 100.0 FR France \n", + "8 101.0 FR France \n", + "9 140.0 FR France \n", + "10 220.0 FR France \n", + "11 331.0 FR France \n", + "12 431.0 FR France \n", + "13 523.0 FR France \n", + "14 548.0 FR France \n", + "15 398.0 FR France \n", + "16 405.0 FR France \n", + "17 370.0 FR France \n", + "18 326.0 FR France \n", + "19 323.0 FR France \n", + "20 220.0 FR France \n", + "21 177.0 FR France \n", + "22 144.0 FR France \n", + "23 127.0 FR France \n", + "24 96.0 FR France \n", + "25 81.0 FR France \n", + "26 63.0 FR France \n", + "27 80.0 FR France \n", + "28 114.0 FR France \n", + "29 131.0 FR France \n", + "... ... ... ... \n", + "2084 59.0 FR France \n", + "2085 64.0 FR France \n", + "2086 97.0 FR France \n", + "2087 93.0 FR France \n", + "2088 80.0 FR France \n", + "2089 116.0 FR France \n", + "2090 149.0 FR France \n", + "2091 281.0 FR France \n", + "2092 395.0 FR France \n", + "2093 485.0 FR France \n", + "2094 544.0 FR France \n", + "2095 689.0 FR France \n", + "2096 722.0 FR France \n", + "2097 762.0 FR France \n", + "2098 926.0 FR France \n", + "2099 1113.0 FR France \n", + "2100 1236.0 FR France \n", + "2101 832.0 FR France \n", + "2102 459.0 FR France \n", + "2103 207.0 FR France \n", + "2104 190.0 FR France \n", + "2105 198.0 FR France \n", + "2106 224.0 FR France \n", + "2107 266.0 FR France \n", + "2108 219.0 FR France \n", + "2109 176.0 FR France \n", + "2110 163.0 FR France \n", + "2111 195.0 FR France \n", + "2112 308.0 FR France \n", + "2113 213.0 FR France \n", + "\n", + "[2114 rows x 10 columns]" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data = pd.read_csv(data_url, skiprows=1)\n", "raw_data" @@ -85,9 +1043,73 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "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
18771989193-NaNNaN-NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1877 198919 3 - NaN NaN - NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1877 FR France " + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -101,9 +1123,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "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
020251832065815191.026125.03123.039.0FRFrance
120251732257217604.027540.03427.041.0FRFrance
220251632856422550.034578.04334.052.0FRFrance
320251533572129592.041850.05344.062.0FRFrance
420251433757931232.043926.05647.065.0FRFrance
520251333967333686.045660.05950.068.0FRFrance
620251235254345627.059459.07868.088.0FRFrance
720251135946952154.066784.08978.0100.0FRFrance
820251036033453048.067620.09079.0101.0FRFrance
920250938453174994.094068.0126112.0140.0FRFrance
102025083136020124824.0147216.0203186.0220.0FRFrance
112025073208952195988.0221916.0312293.0331.0FRFrance
122025063273519258159.0288879.0408385.0431.0FRFrance
132025053334395318416.0350374.0499475.0523.0FRFrance
142025043350043332885.0367201.0522496.0548.0FRFrance
152025033252772238917.0266627.0377356.0398.0FRFrance
162025023257247242991.0271503.0384363.0405.0FRFrance
172025013231549214627.0248471.0345320.0370.0FRFrance
182024523201726185870.0217582.0302278.0326.0FRFrance
192024513201697187843.0215551.0302281.0323.0FRFrance
202024503136694126369.0147019.0205190.0220.0FRFrance
21202449310848799037.0117937.0163149.0177.0FRFrance
2220244838738178687.096075.0131118.0144.0FRFrance
2320244737628667626.084946.0114101.0127.0FRFrance
2420244635639949006.063792.08574.096.0FRFrance
2520244534734740843.053851.07161.081.0FRFrance
2620244433603930122.041956.05445.063.0FRFrance
2720244334657239928.053216.07060.080.0FRFrance
2820244236778560009.075561.010290.0114.0FRFrance
2920244137943571386.087484.0119107.0131.0FRFrance
.................................
208419852132609619621.032571.04735.059.0FRFrance
208519852032789620885.034907.05138.064.0FRFrance
208619851934315432821.053487.07859.097.0FRFrance
208719851834055529935.051175.07455.093.0FRFrance
208819851733405324366.043740.06244.080.0FRFrance
208919851635036236451.064273.09166.0116.0FRFrance
209019851536388145538.082224.011683.0149.0FRFrance
20911985143134545114400.0154690.0244207.0281.0FRFrance
20921985133197206176080.0218332.0357319.0395.0FRFrance
20931985123245240223304.0267176.0445405.0485.0FRFrance
20941985113276205252399.0300011.0501458.0544.0FRFrance
20951985103353231326279.0380183.0640591.0689.0FRFrance
20961985093369895341109.0398681.0670618.0722.0FRFrance
20971985083389886359529.0420243.0707652.0762.0FRFrance
20981985073471852432599.0511105.0855784.0926.0FRFrance
20991985063565825518011.0613639.01026939.01113.0FRFrance
21001985053637302592795.0681809.011551074.01236.0FRFrance
21011985043424937390794.0459080.0770708.0832.0FRFrance
21021985033213901174689.0253113.0388317.0459.0FRFrance
210319850239758680949.0114223.0177147.0207.0FRFrance
210419850138548965918.0105060.0155120.0190.0FRFrance
210519845238483060602.0109058.0154110.0198.0FRFrance
2106198451310172680242.0123210.0185146.0224.0FRFrance
21071984503123680101401.0145959.0225184.0266.0FRFrance
2108198449310107381684.0120462.0184149.0219.0FRFrance
210919844837862060634.096606.0143110.0176.0FRFrance
211019844737202954274.089784.013199.0163.0FRFrance
211119844638733067686.0106974.0159123.0195.0FRFrance
21121984453135223101414.0169032.0246184.0308.0FRFrance
211319844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2113 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202518 3 20658 15191.0 26125.0 31 23.0 \n", + "1 202517 3 22572 17604.0 27540.0 34 27.0 \n", + "2 202516 3 28564 22550.0 34578.0 43 34.0 \n", + "3 202515 3 35721 29592.0 41850.0 53 44.0 \n", + "4 202514 3 37579 31232.0 43926.0 56 47.0 \n", + "5 202513 3 39673 33686.0 45660.0 59 50.0 \n", + "6 202512 3 52543 45627.0 59459.0 78 68.0 \n", + "7 202511 3 59469 52154.0 66784.0 89 78.0 \n", + "8 202510 3 60334 53048.0 67620.0 90 79.0 \n", + "9 202509 3 84531 74994.0 94068.0 126 112.0 \n", + "10 202508 3 136020 124824.0 147216.0 203 186.0 \n", + "11 202507 3 208952 195988.0 221916.0 312 293.0 \n", + "12 202506 3 273519 258159.0 288879.0 408 385.0 \n", + "13 202505 3 334395 318416.0 350374.0 499 475.0 \n", + "14 202504 3 350043 332885.0 367201.0 522 496.0 \n", + "15 202503 3 252772 238917.0 266627.0 377 356.0 \n", + "16 202502 3 257247 242991.0 271503.0 384 363.0 \n", + "17 202501 3 231549 214627.0 248471.0 345 320.0 \n", + "18 202452 3 201726 185870.0 217582.0 302 278.0 \n", + "19 202451 3 201697 187843.0 215551.0 302 281.0 \n", + "20 202450 3 136694 126369.0 147019.0 205 190.0 \n", + "21 202449 3 108487 99037.0 117937.0 163 149.0 \n", + "22 202448 3 87381 78687.0 96075.0 131 118.0 \n", + "23 202447 3 76286 67626.0 84946.0 114 101.0 \n", + "24 202446 3 56399 49006.0 63792.0 85 74.0 \n", + "25 202445 3 47347 40843.0 53851.0 71 61.0 \n", + "26 202444 3 36039 30122.0 41956.0 54 45.0 \n", + "27 202443 3 46572 39928.0 53216.0 70 60.0 \n", + "28 202442 3 67785 60009.0 75561.0 102 90.0 \n", + "29 202441 3 79435 71386.0 87484.0 119 107.0 \n", + "... ... ... ... ... ... ... ... \n", + "2084 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "2085 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "2086 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "2087 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "2088 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "2089 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2090 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2091 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2092 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2093 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2094 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2095 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2096 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2097 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2098 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2099 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2100 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2101 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2102 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2103 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2104 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2105 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2106 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2107 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2108 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2109 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2110 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2111 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2112 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2113 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 39.0 FR France \n", + "1 41.0 FR France \n", + "2 52.0 FR France \n", + "3 62.0 FR France \n", + "4 65.0 FR France \n", + "5 68.0 FR France \n", + "6 88.0 FR France \n", + "7 100.0 FR France \n", + "8 101.0 FR France \n", + "9 140.0 FR France \n", + "10 220.0 FR France \n", + "11 331.0 FR France \n", + "12 431.0 FR France \n", + "13 523.0 FR France \n", + "14 548.0 FR France \n", + "15 398.0 FR France \n", + "16 405.0 FR France \n", + "17 370.0 FR France \n", + "18 326.0 FR France \n", + "19 323.0 FR France \n", + "20 220.0 FR France \n", + "21 177.0 FR France \n", + "22 144.0 FR France \n", + "23 127.0 FR France \n", + "24 96.0 FR France \n", + "25 81.0 FR France \n", + "26 63.0 FR France \n", + "27 80.0 FR France \n", + "28 114.0 FR France \n", + "29 131.0 FR France \n", + "... ... ... ... \n", + "2084 59.0 FR France \n", + "2085 64.0 FR France \n", + "2086 97.0 FR France \n", + "2087 93.0 FR France \n", + "2088 80.0 FR France \n", + "2089 116.0 FR France \n", + "2090 149.0 FR France \n", + "2091 281.0 FR France \n", + "2092 395.0 FR France \n", + "2093 485.0 FR France \n", + "2094 544.0 FR France \n", + "2095 689.0 FR France \n", + "2096 722.0 FR France \n", + "2097 762.0 FR France \n", + "2098 926.0 FR France \n", + "2099 1113.0 FR France \n", + "2100 1236.0 FR France \n", + "2101 832.0 FR France \n", + "2102 459.0 FR France \n", + "2103 207.0 FR France \n", + "2104 190.0 FR France \n", + "2105 198.0 FR France \n", + "2106 224.0 FR France \n", + "2107 266.0 FR France \n", + "2108 219.0 FR France \n", + "2109 176.0 FR France \n", + "2110 163.0 FR France \n", + "2111 195.0 FR France \n", + "2112 308.0 FR France \n", + "2113 213.0 FR France \n", + "\n", + "[2113 rows x 10 columns]" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -129,7 +2118,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, "outputs": [], "source": [ @@ -159,10 +2148,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 13, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -186,9 +2173,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "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", @@ -206,10 +2201,34 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 16, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXm8HFWZ//95ertrbvaEbJAEoixhjwHZXEAWRdFRfMUZBRWNOsyMzszPGXAZRvjiwOiIoIIyArKogAwOjAgYwr6GhC1kIwGyh+Qm9ya59+Zu3X1+f9Q51aeqTi1dXd19Ozzv1+u+um9VnaVPVZ3nPMs5h4QQYBiGYZhKSdW7AgzDMMz+AQsUhmEYJhFYoDAMwzCJwAKFYRiGSQQWKAzDMEwisEBhGIZhEoEFCsMwDJMILFAYhmGYRGCBwjAMwyRCpt4VqCUTJkwQM2fOrHc1GIZhGoply5btFEJMDLvuXSVQZs6ciaVLl9a7GgzDMA0FEW2Ich2bvBiGYZhEYIHCMAzDJAILFIZhGCYRWKAwDMMwicAChWEYhkmEUIFCRDcT0Q4iel07No6IFhHRWvk5Vjt3KRGtI6I1RHSWdvx4Ilouz11HRCSPNxHRXfL4C0Q0U0tzoSxjLRFdqB2fJa9dK9PmKm8KhmEYphKiaCi/AXC269glABYLIeYAWCz/BxEdDmABgCNkmuuJKC3T3ABgIYA58k/leRGAbiHEIQCuAXC1zGscgMsAnABgPoDLNMF1NYBrZPndMg+GYRimjoQKFCHEkwC6XIfPA3Cr/H4rgE9qx+8UQgwKId4GsA7AfCKaAqBDCPGcsPYcvs2VRuV1D4DTpfZyFoBFQoguIUQ3gEUAzpbnPiyvdZe/37F88x68uml3vavBMAwTStyJjZOFENsAQAixjYgmyePTADyvXbdZHhuW393HVZpNMq88Ee0BMF4/7kozHsBuIUTekNd+x8d//jQAYP1VH6tzTRiGYYJJ2ilPhmMi4HicNEF5eStEtJCIlhLR0s7OTr/LGIZhmAqJK1C2SzMW5OcOeXwzgBnaddMBbJXHpxuOO9IQUQbAaFgmNr+8dgIYI6915+VBCHGjEGKeEGLexImhS9EwDMMwMYkrUO4HoKKuLgRwn3Z8gYzcmgXL+b5Emsd6iOhE6QO5wJVG5fUZAI9KP8vDAM4korHSGX8mgIflucfkte7yGYZhmDoR6kMhot8D+CCACUS0GVbk1VUA7iaiiwBsBHA+AAghVhDR3QBWAsgDuFgIUZBZfQNWxFgLgAflHwDcBOB2IloHSzNZIPPqIqIrALwor7tcCKGCA/4VwJ1E9P8AvCzzYBiGYeoIWQP+dwfz5s0Tjbba8MxLHgDATnmGYeoHES0TQswLu45nyjMMwzCJwAKFYRiGSQQWKAzDMEwisEBhGIZhEoEFCsMwDJMILFAYhmGYRGCBwjAMwyQCCxSGYRgmEVigMAzDMInAAoVhGIZJBBYoDMMwTCKwQGEYhmESgQUKwzAMkwgsUBiGYZhEYIHClM2L67sw85IHsG5HT72rwjDMCIIFClM2f3rV2nH5qbU761wThmFGEixQGIZhmERggcIwDMMkAgsUpmyIqN5VYBhmBMIChYmNEPWuAcMwIwkWKAzDMEwisEBhYsMKCsMwOixQGIZhmERggcIwDMMkAgsUpmw4yIthGBMsUJjYCA7zYhhGgwUKUzYEVlEYhvHCAoVhGIZJBBYoDMMwTCKwQGlg+gbzOPzfHsJja3bUuyoMwzCVCRQi+kciWkFErxPR74momYjGEdEiIlorP8dq119KROuIaA0RnaUdP56Ilstz15FcLIqImojoLnn8BSKaqaW5UJaxlogurOR3NCrrdvRi31ABP/nLGzUtl6O8GIYxEVugENE0AP8AYJ4QYi6ANIAFAC4BsFgIMQfAYvk/iOhwef4IAGcDuJ6I0jK7GwAsBDBH/p0tj18EoFsIcQiAawBcLfMaB+AyACcAmA/gMl1wvVuod8fOQV4Mw+hUavLKAGghogyAVgBbAZwH4FZ5/lYAn5TfzwNwpxBiUAjxNoB1AOYT0RQAHUKI54QVh3qbK43K6x4Ap0vt5SwAi4QQXUKIbgCLUBJCTJVhBYVhGBOxBYoQYguAHwPYCGAbgD1CiL8AmCyE2Cav2QZgkkwyDcAmLYvN8tg0+d193JFGCJEHsAfA+IC8mBoieDUvhmE0KjF5jYWlQcwCMBVAGxF9PiiJ4ZgIOB43jbueC4loKREt7ezsDKgeE5V6m9oYhhmZVGLyOgPA20KITiHEMIB7AZwEYLs0Y0F+qhCkzQBmaOmnwzKRbZbf3ccdaaRZbTSAroC8PAghbhRCzBNCzJs4cWLMnzqyYU2BYZiRQCUCZSOAE4moVfo1TgewCsD9AFTU1YUA7pPf7wewQEZuzYLlfF8izWI9RHSizOcCVxqV12cAPCr9LA8DOJOIxkpN6Ux57F0Fz1hnGGYkkYmbUAjxAhHdA+AlAHkALwO4EUA7gLuJ6CJYQud8ef0KIrobwEp5/cVCiILM7hsAfgOgBcCD8g8AbgJwOxGtg6WZLJB5dRHRFQBelNddLoToivtbmHhwlBfDMDqxBQoACCEugxW+qzMIS1sxXX8lgCsNx5cCmGs4PgApkAznbgZwc5lV3i+pdcfOe8ozDGOCZ8o3MPXu11lBYRhGhwXKfkDNNZTaFscwTIPAAoVhGIZJBBYoDMMwTCKwQGFiE8fUJoTA3/52GZ59c2fyFWIYpq6wQGHKpwInyr6hAv68/B1c9JulydWHYZgRAQuU/YB6RVtVMkO/3hFqDMMkDwuUBqZenTLP0GcYxgQLFKam8NwVhtl/YYGyHyAacA0U1nEYZv+DBQoTmwaUYwzDVBEWKA1MvXwZ7FBnGMYEC5QGphE79kY0zzEMEw0WKEzZVCLH7K04G1EaMgwTCAsUpqawgsIw+y8sUBqYhuycG7HODMNEggXKfkC9BEscf4iaXc8GL4bZ/2CB0sBUsvRJJVTi/mhIrYphmEiwQGlg6t05x1ptOPlqMAwzQmCBwpRNJfNfOGyYYfZfWKCMYMI6X3W6XqavONg1ZScKw+x3sEDZD2ikQX8j1ZVhmPJggTKCCet8662ZxCm93nVmGKZ6sEBpYOo12q9okrusM1u8GGb/gwXKCCZMXtR7rM9RXgzD6LBAaRBMDvp6RUxVpKCEVHkwX0Bnz2AFJTAMUy9YoJTBo6u349VNu2tWXlSB0Uij/jAfytduX4b3XflIjWrDMEySZOpdgUbiy79ZCgBYf9XHal62EF7fRSMJEoWSkX6rDT++prOGtWEYJklYQxnBCJ/v9rF6z5SPIdIaUQgyDBMNFigNTeOFefFMeYbZf2GBMoLR+16zU76GlTEQK8rLNnklWxeGYepPRQKFiMYQ0T1EtJqIVhHR+4loHBEtIqK18nOsdv2lRLSOiNYQ0Vna8eOJaLk8dx1JAzsRNRHRXfL4C0Q0U0tzoSxjLRFdWMnvaHRqPepnWcAwjIlKNZRrATwkhDgUwNEAVgG4BMBiIcQcAIvl/yCiwwEsAHAEgLMBXE9EaZnPDQAWApgj/86Wxy8C0C2EOATANQCulnmNA3AZgBMAzAdwmS649hd0H4XRhxJwrhbEminPFi+G2W+JLVCIqAPAaQBuAgAhxJAQYjeA8wDcKi+7FcAn5ffzANwphBgUQrwNYB2A+UQ0BUCHEOI5YQ21b3OlUXndA+B0qb2cBWCREKJLCNENYBFKQuhdQyN2zrz0CsPsv1SiocwG0AngFiJ6mYh+TURtACYLIbYBgPycJK+fBmCTln6zPDZNfncfd6QRQuQB7AEwPiCv/QqnD8V0vr6dcxzTl+1DCb2OBQ/DNBqVCJQMgOMA3CCEOBZAH6R5ywdTHyICjsdN4yyUaCERLSWipZ2d+9cch3p3ufEWh4x4Xb1/HMMwZVOJQNkMYLMQ4gX5/z2wBMx2acaC/NyhXT9DSz8dwFZ5fLrhuCMNEWUAjAbQFZCXByHEjUKIeUKIeRMnTozxM0cGgaaiGne+lW0BvP/N/mcYxiK2QBFCvANgExG9Vx46HcBKAPcDUFFXFwK4T36/H8ACGbk1C5bzfYk0i/UQ0YnSP3KBK43K6zMAHpV+locBnElEY6Uz/kx5bL/FbPKqfT0qrYCteoZIJTZ5MUzjUenSK38P4LdElAPwFoAvwRJSdxPRRQA2AjgfAIQQK4joblhCJw/gYiFEQebzDQC/AdAC4EH5B1gO/9uJaB0szWSBzKuLiK4A8KK87nIhRFeFv2XEMVL3Q6lsC+CI18UugWGYelGRQBFCvAJgnuHU6T7XXwngSsPxpQDmGo4PQAokw7mbAdxcTn33Oxqy141W6SJrKAzTcPBM+RFMmAbSiF1u9CivqldlxPDsmzsx85IHsH5nX72rwjAVwQKlQQjqYBtqYmPitWh87n1pCwBgydv7ndWWeZfBAmUEE+pDEeqzxkuvVBDlpUxZYXnsLxqKEAIPvf4OikX/H8RL2TD7CyxQRjDO5esNi0PWebxfyeKQodftJ7rMPcs24+t3LMNvX9hQ76owTNVhgdLA1GsUn8wWwGFhwxUUMoLYIbcz3rpnIPTaWgvRzd378PqWPTUtk9m/YYEygtFNWcZ5KDWsS1JE7TSr8dsWr9qOXzy2rgo5hxMkIOu1lP8pVz+Gc3/2dH0KZ/ZLeAvgBiGog20kwRLZ5FUFFeWiW60tnC/+0CGJ5+2HEhZRBOn+opUx715YQxnBhPUv9Z5NXomJJtQpHzvnkYU9CTRIQ2G3fGRe2tiNmZc8gFc37a53VRgDLFAaBOOOjXWoB6CNuqvplK/gx23q2hc/ccJQuDxhyuDJN6wFXh9Ztb3ONWFMsEAZwYR2qg3YS0XWamL+todXvINT//MxPLJyZHQ4SveIok024O2sOa05a0++/qFCyJVMPWCB0iCYd2ys01peFXiRqx02/Io0hazZ3hMrfdJEWQyzEo3v3UZL1hIo+4ZZoIxEWKCMZCJPbKx+VZzlxi8waKMbZxnx8s8XigCATMq/BHVNLYiy1Ey9orwakXTK6rLq7T9kzLBAaRCCl16pl6ZSfppq74eSlzPSM2n/R3u4EGPZfSHw6OrtKATMeDfBPpTqwPJkZMICZQQTujhkA75U0XdsjPfj8lJYBGkoQ/nyNZSHV7yDL/9mKX791FtlpSvHhxKHgeEC/vjyZh6xMyMCnofSKOwvExuVCahKYcPDyuSVDhAoMUxenXLG+6bu8iLIUvKHRunv42iaP3p4DW56+m2Mbc3hg++dVHZ6hkkS1lASQgiBXzy2Dht3JReyqndCxrW86j0PJVbxEU1eMX+aEhbZAJNXHIGinOrlJlWCM9hSFl3ouNm+11rSZe9AvvzEdeLZN3fiuCsWoWdguOy07G8a2bBASYgdPYP40cNr8MVbltSszPrNQxm5UV6Dw1aPr6KBTAzHMHmlU6rTj1evavu56j24KIefLlqLrr4hrNy6N3YeDfRz31WwQEkI9YD3DiY3UnSsNhzklG/AKK/kLnTSL8NJm4MESgwNRblkynXKl2PyikMlwr1e5DJWtzMYQ7AzIxsWKDEw7W1Rj2ieeoUNK+JFecm0YasNx6gPYDmpASDA4hVifjKjOu5y05bmmATsh1LBs1MtcbL6nb1lC8+oNLFA2W9hgRKDgqFzqEY0j2O1YfMViZVVKyKHDcedhyI7waTDrEuaRkwNJeCaJIRCkoOKlzZ24+yfPoVbnnk7uUw1mrJWtxPLh5J0ZZhEYYESg6Lh7VVCppbaQiPakSOHDVcoLIPaphhjYKw0HtO9j0KkdDHyLmc146gs32ztkfJ2lfa4nzG2FQCwubs/dh77ywZs+xssUGJg6pCUeSBuh2PC6UMZOYtD2uVXsDhk1bYAFo4PI3HukdI0yp0TKao80KjGiF35mJSvI2lUvo04IGKCYYESA5PJSwmZar0jQfnWfk/5CqK8qrzBlso/qE3iNJcSKEF7w5vr4/w0kYRfPclHgKodSFCdbJkRAAuUGJicldUweYXlVa8RXkUCLGrYcIU/LlAAxxBXcU1LpcCJ8HSxnPJV6Pxr1eFXch+YkQkLlBiYRqnVMHnpmLcAlkKsKiWWV5fIaWtQhpXeP4MqBS+Z66E+A8qkCiY22sEg5SetHwkIQTaXjUxYoMTAJDTUsSQHUCN1LS9VbGVhw9Gui5t/oFM+jvM75p2ttg8FEcKSy86yCnk68q8oLasoIxkWKDEw+VCUhlKtiWbGpVeqUlI4lWkoUX0o8QqJ4rOox6zyaplqbO0m0TyReJ4mWMnY/2CBEoOgKK9EO47I/oYEyyyDkboFcFj6eoR2Bz0XlTwzVIXev9qz79kPsv/CAiUGxigveSxVrbfF5ENR5pQaj/UqKa9kLqvOTPlS+ur4UMoVRuW0VRzNqZp9c9UFb4wC6uU3ZKLBAiUGQU75OigoNacyZ2pEk1fcRRiFCo7wniut/Fu7lo1SVCIz5RN8WqoxWdKRfwVmOnbGj2wqFihElCail4noT/L/cUS0iIjWys+x2rWXEtE6IlpDRGdpx48nouXy3HUkh69E1EREd8njLxDRTC3NhbKMtUR0YaW/oxxMYcO2U752Ckr9tgCuQdr481DkZ2CUVw0Fiv2tOg9GNfajLy0jlFyeSVGs0zPPRCMJDeWbAFZp/18CYLEQYg6AxfJ/ENHhABYAOALA2QCuJyK1JOwNABYCmCP/zpbHLwLQLYQ4BMA1AK6WeY0DcBmAEwDMB3CZLriqjXHpFduvklzH4dgPJSBsuObYprY4acsqojrECc+NeVvL+R2x5qEk4JT3CN8R7OTgJVdGNhUJFCKaDuBjAH6tHT4PwK3y+60APqkdv1MIMSiEeBvAOgDziWgKgA4hxHPCerJvc6VRed0D4HSpvZwFYJEQoksI0Q1gEUpCqOqYBUp1NZQgaj4PRX3GccpXWUcJ0trUranlPJQommslkxOTeN7c7VHtKK9K6lzLe8eUT6Uayk8B/AsAPe5pshBiGwDIT7Uv6TQAm7TrNstj0+R393FHGiFEHsAeAOMD8qoJpu00qj0Pxbxjo/XZ2TOItzp7Eyw5GnFGizWL8jLUrbQEfQVBBTGTVnuckaRfqxpmNHO5tUrE1IrYAoWIzgWwQwixLGoSwzERcDxuGmehRAuJaCkRLe3s7IxU0TCMS6/UQUPR362/+93LtS+3mmHD5WftSBe0onAcgRJ/omV1n4skHOg113CVFhmj5GIFaSuhWBS49pG12LOv/CX3301UoqGcDOATRLQewJ0APkxEdwDYLs1YkJ875PWbAczQ0k8HsFUen2447khDRBkAowF0BeTlQQhxoxBinhBi3sSJE+P9UjhHcUHL1yc5kzfch1IfKgndLIUNh1xXsYYS75wfcbWaavtQktB9/H9bdZ6wioRfnTSUx9/YgWseeQOX3f96XcpvFGILFCHEpUKI6UKImbCc7Y8KIT4P4H4AKurqQgD3ye/3A1ggI7dmwXK+L5FmsR4iOlH6Ry5wpVF5fUaWIQA8DOBMIhornfFnymNVQ3+OjUuv1EVDqc/LVZOw4bidjr3UicHkVWYdnPWJWR27bP8HI4mlTpJcF6uStcWiUEmkVr0GUUNyd8l9Q4U61aAxyFQhz6sA3E1EFwHYCOB8ABBCrCCiuwGsBJAHcLEQQt2dbwD4DYAWAA/KPwC4CcDtRLQOlmayQObVRURXAHhRXne5EKKrCr/FRn+QA01eVSoz7EWKI8jue2ULTj9sMtqbynsMooTmhqUNvS6uPIlQTpwNtuLPi4mVLDJh2wd39w0hl0mhLeAeu+tY2gIggQoGFBjHwV6McoOZupGIQBFCPA7gcfl9F4DTfa67EsCVhuNLAcw1HB+AFEiGczcDuDluncslzORViuapnYpSyTv1+pY9+Oadr+DjR0/Fzz53bHnlVjLCrJFTPqhxaupDQYQoL8QPFrCz9Ul77BWLMKG9CUu/d0ZoHUOySoyShhLnPtRXkozgiOoRAc+Uj4hTQ/GeNx2ruEzt5TG+SBW8W4NShd/UtS92HjHH7AAi+FAqHIIGpa/EhxJ3P5Skr1VE6eB29g6WVW61l/RR+VYi2FlBGZmwQImI/uwP5r12VNspn+AIppoTG7Npq6L5GPafSkwi1dZQVLogc0q8kXHM+pRxbUVrjMVP6unYSybNCjINIMo98k1bJ1HC0crRYIESA5Njrj5O+dL3csvNpKxbny93k3TA7nEqCfus1r4WUYRdvI4sJhF6okrWGEvCge5Oqp7lavWhRVugxH9+6gXvxxIMC5SI6J1nv0GglJzyjfHAKQ1lOIatrpJ3OopPAUgibNg/g3gdWTytrJzRfqzVhpPYj971CFS7z66FhsvUBxYoEdEf5IFhf5NXqpYaSgVpUyll8orvGK2qyasKS6+UWwdXhWIRpawkloSpKOTY9eOqvQBjyeRViWBnyRLGyxu7ccIPH6npZEwWKDEw7odShR0bQ30oFXVA1qcpBDpq2ljl1qCMsHIq6cjKr0d0B3csp7xdTnxq7pS3w4YbRyg0Tk1L/PzRddi+dxBL1ld1RoUDFigR0Z994zyUOrwclb3wlde3mmGfcWtX0lACorxq6EMpp6xYPpQEBjAep7xtp6s4ayOVhJ2b9iKqJY0UNtycsxZz7zdYVKoFC5SI6J23cT+UqkxsFMbv9jHdKV9myUnMVq7k1Q6rbeyJhHZ6/2tqq6GE16d0TQVmqySd8raGUh1KTvny09ZLnDSQMmXTkrUEykANZ/ezQIlImIZiH6rWBlsJP9BRRvLhaeOXG3pd+Vm7yklYQ6miD6WSMNokdqH0CxuuFiWnfAU+lERrFJ1G0lDswJs4S0PEhAVKRPQHuGZLr4S8NZUZvJJz4lYjbaVbACftQ4mvMUWvT5x2TSSq0FVstR3flTjlKxnMJEEjaSql7RpqVyYLlBiYIqMqXXpFCIEdewf8z5sTxSoLKK1nFSeH0mzl+B1CWDuF/bT/fGg1HnhtW6z0NTXzleVDKTfzZDQUd8pqz0YvOeXjp22gfr1uqIjTWkbEsUCJiGMtrypoKDc/sx7zf7gY63aUNsoKewz08+XKsSQWAKynyev6x9/Exb97qaz06lht1/IKT59E1JOpc47akXid8pU/G0FUdB/UZ51UhUYyeVV71WgTLFAiot8Tk4ZS6dIrz6zbCQBYv7PPXL7hqUhCGNR6Y6aoaSr1WQT6UGLkG38/lPB0Scz7MKWNGhLuDRsOT3P+L5/FVQ+ujpS/m7iTRCtNWwmNuJc9aygjGP2eBO6HElNHScu7r4cfl/MgDJe5hEpljvX4L3XksOEKXwJT8tIEwhpqKHa6IAEnNZQKjN1Bm76Vm7a0K6I/L67vxi+feDNq9Rwk4UMJSvtX1z+D6x9fF6dqoTTKShgA+1BGNtpNMWoo0ifhp6Fs7t7nq30AQFom9BtVGk042ktlWrAyiCR2XYyTWqUIDRsuO2eVTjg+TcTbhyOeEI1yebV8FrE1lAqisKJQSdhwlHXAXtq4G//50Jqy8n3o9Xdwwc1LfM83kjPeTS2rXo0NtvZL9A7K6EMJeeJOufoxAMD6qz5mPJ9OewVKOT6UctfkSiL0N9ZLFtWHUrHJy7/oWpoAorRVJYsllsoxDXLi5RdFQ6mM+D6j0tL3iVYIX79jmZW/EMaAkUaUJylSPhQ2eY04QuehqBVaY967TCpEQzF1kNqx4Xx5BZdeZv903X1Dtm/HUW4twoaDTEQRGjnoikps9+US7ffGf3aC1gGLKlA8vy2BsOFdvYP4v1e3Gs9F8XP5Ue2wYb8mq/cM/TiUtpauXZksUCKi35Mgp3zcjkf5UHRNw5mVN9+wOgVRGq37X/PFW5bgb379gmcxzCS0m0psXlH8RUH3oR4+lEANRd7yyjQU77G4Ji+VLK6GAwBfvW0p/v73Lxs391K/s5KtmKs16vb7zfa9aRwXSskpX0P9igVKDIKc8nE7BVs9jV2rMgVKBLPGmu09APx/U6x5KAlcF2TeCxKUYav6PvT6O/j2H141nqvmADWJEO4knfLlmJX87sWW3f2+5ytxyoeZBysVNH4CpRLhWi/YKT+C0R9Us1O+8k7BysD8j9nkVTpY7kMTZaSnIlrcGkEU7ca/3GjXBb3AbwcEN0QJF/a75ut3LMMflm32SRtuIgyqT7AJz/qsaF6G4VxkDcX1f9ie77r5Z0t3v/GaoDkQFa3lFZK23GhHN347mDagglJaiZoFyshDvydBTvm4967SB7XckVmU7lHZYPOuUWYlUUlRtZqg0fW5P3s6IH/5GSRYqigI/eoTRBIdrKmC8U1ewRqKLvj87pNtvzeVV0EUWdhaXkFbWq/d3oM/vWb269j5+yRvpKX2FSUNpXZ15yiviOj3xLj0SoUmL7scH63E+GKGnA8sJ4JtX5nhvL83vjYWdaQX2wka4XfVdC2vCPVJooOtxCnvKVcEp9cPh90n42+qYEASmC+CNZSPXPMkAODco6Z6zhFZ98hPINVje4pKqcesftZQIhI1bDiss/J7EeJEZPgJn0hpI5m8LNx28ErerahJK7VZB6WuSBMoux4REoaYvJZt6MIb0p9VTjnxTV7Bz3IkDUXlbTR5xR98haWNEj4/lPdeo+obRYg2Css37wFQ2wg1FihRCdFQ1HMcFrni98CaZuDqV4a9e+W+nJHmGtgmL+dVJZNX9XrmSjW9SqK8TC9gbIUpQjuHTZr89A3P4Uw5uvbm76+hBLaBPt/JraCEaFV6vu5nQ0G+2m2pLSryofi8Z3710ekbzHuOBdXXKq/xJMrTMuSfnfIjEP2eBG0BHDYqDFOdy7n3jqzK1VAi2B2UiPO1S4eU+fLGbhx/xSLs3jfkSRK22nCZ8zQ9+Vcij0z3KPYe92VcE0crC1oHLCiUXP+NfkuvRDJ5hTR0wfDshDn9g0hCQzG1S7iGYh1PcovvMPKFIi64eQmWxtzCt6PZ8mjU0lzHAiUi+j0pGEZB6qaF3Ty/vtn0nDp9JP4jPff3SEQYOauXxxvlJSKV+fNH12FX3xBeXN9dShuxonFfgkj7oYQJ/QANpXzTYng6dW44hkAJWmomSEAVgjR2PByVAAAgAElEQVSUkJnsDpNXaFt6j1WyurJwfbqJMh/LVG4qZOmjeoQNb9szgCff6MS37nolVvpjDxwLwCzUqwULlBgEaShho65yOsowH4lD4MQ1eQWkS4WZvELKNAvJiCavCpcNCerBw7I2dh4hneC+oTyufWStYYQcvYNzR9NFIUhgBXWC+u/w21PerwmFVs2wjtak3domr1gTG9WnudwobWiss3rWy1ipolakYmpFylcUV9uPAwuUiER1yoeP2EJMXsL8PbROMU1ekTQU15sfNkr0lKX9kKhp4o4Iy9kh0bds04AhRKP7xWPrcM0jb+DupZuc9bHT+ZepzsSZQ1GSn/E1FE+e5Tjl/XyC5H++Ig0lwGcERGtDU53CVqK2TV7hVUyMkpktXno1uKll2DALlIiEhQ2rhzSuQCnF7ZtHjuUIlyhEMcXYPhRfDSWsFEOggfA74ySKJpcyZBK0jEnUZTtMJs2w2eOqjfb0D7vKDCxKZQ6g/AU+rfzVb/KeC+pIgkxetg/FJ71+POw+mQWKudwo2E55Pw0lgtpjqpMdIu8jkCqdZxaHYsR3xY8h+TxFCVRIChYoEdFviXHplZCRk8Jf4HhnFocKkTKudRNlH3O/iY3uPMIQhu9hKaOYvDIp7+Mb3MFGK9u8xbMzfzdN2TQAYGDYrc2FF2qbvGLYgEodrPdcUEeiP4flLr2i35s4GkoSYcN+Sf2EsggRgkH1DSqvmqjfWqnJizWUEYjjgaxAQwkNWfU5H6aB+J0dGC7g90s2ejrCkqnEP8+Sycvc4YQ9p5X4UKJoKJm0twA/oeFcpiakLY3ajfPTTdZeLbr8OTu2U77MFaMd9TI55YM0lAAzpOqT/e5VIeRdAEph8EHafDW2APYzeRVChGBYRGNpA73aYf/GgEJXbt2Llzd2G88pgRJnoBKX2AKFiGYQ0WNEtIqIVhDRN+XxcUS0iIjWys+xWppLiWgdEa0horO048cT0XJ57jqSPRkRNRHRXfL4C0Q0U0tzoSxjLRFdGPd3RCXc5CU/Y5gAAPMIKVxD0VUU8zU/fngNLr13ORat3O5KHJjMqpP89GgoIWYHbz1Nx0IEawQNpSnjfXz9RrDOUNfgfI0++ZDIJ7+JqVE0MpWn21cVhaARe6BTXivKb0/5KBFP4WHD3vOq04/jJguzBPhpZfkQgZIKWaakHtNQVJ2DNJSPXvcUPnX9s8ZzyuTVKE75PIB/FkIcBuBEABcT0eEALgGwWAgxB8Bi+T/kuQUAjgBwNoDriSgt87oBwEIAc+Tf2fL4RQC6hRCHALgGwNUyr3EALgNwAoD5AC7TBVe1MU98Uy92TIFi5+PN08rXm0b4XKvTKZcP7xtyTuayrw/UUKxPv8UhwzoT02sQ1X4exSmfSRtMXnY55k4SiOeUD7Pd+81PiKShyM84tu6gEXugUz7g2QrrtPV8/Sc2+tehtJCqOe13/rgcx1z+F3PhIffhlU3m0XqYhuI3iddOE+JDWbahCw+veMfnbDxUXUy+wijYJq9GmCkvhNgmhHhJfu8BsArANADnAbhVXnYrgE/K7+cBuFMIMSiEeBvAOgDziWgKgA4hxHPCesJuc6VRed0D4HSpvZwFYJEQoksI0Q1gEUpCqCo4ZgdXwSlvl6NrKDHrF4UgU4miZLZwm3HKHWFG78wVUfrWtKET94sgctQ1TEMJiEwKD45w/x8+0LBNXjGGkkGLJUadKe9OHWryiqChBE0ULEUfmev2uxc2Yve+YeO5MGH347+8YTweNJEzrL5AeEDHp294Dl+7fZm5UjFR/UycfeyFENjRM+jIpxYk4kORpqhjAbwAYLIQYhtgCR0Ak+Rl0wDoMZWb5bFp8rv7uCONECIPYA+A8QF5meq2kIiWEtHSzs7OeD8Q4dqA/qAHdR5hq7M6tZIQDcXnu+MaecL9UEZ5xMhn1BY0Kjal17Gf7ZB3JMqoasKonG/+7r45aN6FNw/T/bU+yzb7R9FQlMmrgnkopjoHaTx5h1Aw1yda2HBw/UydmTpWDR+KH3rknnGmfNjSKxUEEsRF+ePi+OS7+kqrUzSUU56I2gH8D4BvCSH2Bl1qOCYCjsdN4zwoxI1CiHlCiHkTJ04MqF4wYbH3hYAX1JFPmAbj4ywNMsO4vzuukZ/uhzIsWgYo2W79OrqoA58opjk3Qb6oA8e1evJ15x/kHI/lQwkxtZjKAaL5UNS5sJGkqRMNuo9xw4bD5lTpj0OYw9eUh/qdZfkIJXEFu0OAGuqUlQEepoUj9XJr6N+2Tc1xlntxmCUbRUMhoiwsYfJbIcS98vB2acaC/Nwhj28GMENLPh3AVnl8uuG4Iw0RZQCMBtAVkFfVcHTuIQIlqu1aR2kQfov2GeP5Iwx//UZyUTtIwPtARk1rUtWj2nOjzPIO6vjddS5HQwkKdfVL6ffORxlJ2xqKT2emMEUwlSKmvNcHdSRB7RHWaetCJMynZBQoESfcDRl3ewzWFM6ZewAAYMa4FsfxsPezKWO5c/td210rktqewo8n3ujEd/643HFM1TOOD8Vh4msEgSJ9GTcBWCWE+Il26n4AF8rvFwK4Tzu+QEZuzYLlfF8izWI9RHSizPMCVxqV12cAPCr9LA8DOJOIxkpn/JnyWNUICxuO2mGFRXnpp8McyZEcvvIab6RIcAep18mzwZZdp/Dy3WVE3pY2IPOSY9//PrjbWf8vzsjYnb9vWpjLDUqmqjoU4jgKmiRoupNx1/IK6zz1UbqfySsoDFeZ4sLug8lkFzaY6ZUrCbuL1ethapecjBgc8BMoZfsNy+PCm5fgdy9sdBxTloE4Jq+oA9ykqWSDrZMBfAHAciJSq5d9B8BVAO4moosAbARwPgAIIVYQ0d0AVsKKELtYCKHu3jcA/AZAC4AH5R9gCazbiWgdLM1kgcyri4iuAPCivO5yIUS8JTkjou4JkblTjBpKGTZxyvGi6+XHfChUB+c1eTnLNeEb5RXioHSnd7RNxN8RZXQdJGTd7ezsQEM67jJNS47jwud4ACUBHXzxcLGIFqQdxxIxebkqbQtlP6e8Y3AVrFWZ6jCYL4TWD/ARKCEDoX1DVt7ue+zoYA3lqhD0XsPS9lZdYcw3aYpFgVTKqd3FccoXHWbJBhAoQoin4e9aPd0nzZUArjQcXwpgruH4AKRAMpy7GcDNUetbKeo5yqZSZpOXdijOGkqmeQ4Ok5ep8wyqsCsPj1M+QmLfKC/5GVX9d2pv/uVHDe2N0om627lY9NbBN3/jCrneOvrl70jn+jReE9Hha+pgg4Srfr0QwmGLDwpJf83emMlcD+fo16eyAWG4ajWB0N9rWvrejkDzSeMTQRY2Yh/dkgUAbNi1z5hvrZzyQ4UimlPWoGG4jLBh9/0Ni2qrFjxTPiLqpmTSZDZ5OUbhznNROkrb6SecnYBfnvKCsGqXBIrroQxb7BIoPcge272qa4iD0hQlFrSzZZjPyK5vgNmjtIyJa4RaxgtmrFvYWl7K0exOV8Y9irNybyGggw0SGvmi/7ktu/vlcR8NJWS0DwSH4Q7klRZhTGqso13XkAmmfhFkYRYEdaR/KNiHEhbOXqkGo/uN7DpHsHm5mypsmkO1YIESEXV/MikyrzYc8JLpl/uN6Eqzk33SxdVQlMnLdVzP+w/LNsGEHUrp8xZFdcqbOnNzkEGJIF9L0JwQu3N23wMfU6KJwP1QfFKrNHFCrO1OMkygGH0KZkFm1an03eN4j2Ci9TV56c96SNywqTMbjKyh+N8H3zr7BCnoeQW14z5fH4rzOj/irBito0eZqQFEFA3FHYnpHOCyQBlxqAc4m075POj+KnVYyDFQGu07RpUhEwL1Q188aaYxX18NRUvcO2h+idQ1XpOXf4fuIMCHEmryiqmh2BNMXS+2c2Z4DA0lRCtTy6Z4tIgI73LJ0WyoS0hHGGSK0f0bbuHg56tz5u13XB88ma9RgxFTvYZCJjYqTIuSht0Hv1n4YRqKOu+roUQ0ecWZS+SXXt1vP3mi/8YgjbyWTnkWKBGxNZQ0BT6QQLC6HbbFqKMj1b76dbDpFIEIGNVsdoeVRsjG0wCA5qz5MVAPqdcp76xzGCabe9DkQf06E34RN8WisOvsecECTJJ++euE+TmUAPOEWEeIpgtyguvlmdb6soszFBBlrol1zly7SDPly9SqCkUReXFI02g/akCFO++84RnUUadDo7xCnp1KBYrDPKyc8j4mL31la88Aqk5RXixQIqI6hoyvU95foDhHdH4CxXtevzJoHkqQRuxnAtL/bc44I4fcZXrChiMKFJMdPTBCK+ISLX5hrcOO0FC3CUD7HqqheI9Ftff7tVVgeUrIBph4gOBwdaNTPmCQo7eHXx39+qEok+b8fCh6hxvWNkHvWbkmryBtzaqL1FB856EgsFyFae5MOejtORxi8jrtR495rlU4Qrsr9OuUAwuUiKj7nE2T2eQVMCKIMqKzR8A+o0o/MwwZrnVdZayTnl9z1ixQ/Eb7YQ5qNyaBYt5r3JzG7zr3b9btz+7BbdBy7QCweFVpNeY4e3gUfNrKvj6grYI6SYeGYmi0oAmXpug6d5l+5QL+96Cs1YZd5we1exTHfBSm3fid17MyhTqrei7bYF5cMsrKEoD/TPuo5A1Oeb+w4U65VpeVzvV7A8xh1YQFSgQ27tqHR1dbE/4z6VTgFsCAVy12TgTz01C8nXRQhwBYnQhR8NIMftpEtJGzdZH7xS516CEaimG2dOnFNAtIRdCoyq9zH3aYC5x1DusEn3yjM/C8sM+Z61TSUPwHE37oo2p3uziidUw+lIAoryCnvGPErtXxvF8845vGvj6Cfd5vtWH9WYqzYoFqgzDtyZ21c2KjN53+uJjMXiYLgolKd0d0aCj20ivh6dzvqD7LvpZO+UomNr5rOOOaJ+yRRyZFZlt30fJnFIrC6wCNYPIqGB5Y/Uq/WdLW6EX4Rh+poz5bmgDwX4/Jr5NURNZQHB2QSmvqtMM1Ob1cr0DxOjRL5Wrtaup8AwYEeln+fgW194TzfCmcOEjjcg4c9H3DwsxLUZ3y7vbUnwddEL+6abdWL3N9iyF1ArQ5TK77oN6jdIrCnfLGMOlgTUHVp1wrgd5+XX1DmDqmxXi+2k75vGFQFEWgeJ47WY/mbJp9KCMNXY3NpMnY4RSKwl5gLo7JSz08TpNX8KgaAEBhPhRhLDdK5Jk67rbPlkbrIRqKIf+S/8NUV2/ZJvzMDw6TV9A9MHa+3vydlfPWUacUwOBsqygj1qARv/5vUNSTqRS/EHR3OUMFs9/A7x6ELbRo1Ud17M46q/ZpyqRCNVxT2xVCBLuff80RNhwiUPoMs+WjLr1SqQ9Ff9dsDSXCmnhu4atMi605FigjmnQqZV56RQhk5YZPQU55vxFdSZUvQ0PRrvB7N20NxX2Bo6Py6ziKxvMlM5q5TDfRnfJamoDOxs+M55gUFnAPTFk7fGCBJq9g4evVUIq+ZZbSmuvprtd2zWbuzT9EQwloD9PWw7l0yve3RgkysdvDfY/yJYESrqF4Lwhb+t5vJWM9CsocTaev5+UfTRdHCJaDKTrLpKEM5t2C2tzOrbkMO+VHMpmUN2xYCIHOnkFkUubl3qM4MU3OarcpxIOwtACi8ifr6cLIdw+IojrvfsGCR4l2enk+b+isw8wOwSYvc+ddiYbinDnuX7cwH8pwQLl+BE0y1P/fsXfAkzbIBORwzAaYAE2j6lwm5SsEo4XBy3Ld98jWUNLhnbPJV6k0eZ+k/hqKv/lPHWvL+a84HNXk5e7oy0XvO4LMZ+46uu8vaygNQiZFEMLZ6fx5ubX1Z7fcZc59c6OYl0ojK/NI2q8DTqcocPE4v9G8np3fAn/qJYw7D8XWuuKYvAKyNgkqq54lm3NQ1ItRQzH4eZzn1TmfEbmah+I2ecnr397ZZ0znrluQyavPMAE1aG8RvSru6CP995oikwpFESn6K8xc6p4fYWso2QgaiuFGmDR5x/kIJi+/KL7WnOVSNjrlVb4h8uKzv3ou+IIQ9LopoWASrO46ek1e1nn2oYxwMmk1A7h0rLPHOXJ039woUV4mE4GzgzWN1krL0oeavNwaivavaQJZsSjs3+iJ8lLXhHUIRe8L4edQd2SMaE55vwiiFsNLFLbUSJgG4xfxpgibULn6nR5jOm/dXOd02/6Q17YfFEarDxRUB2PXV7vnpt80VChGW3qljOfZKsv6vzmTDh+QBGmxftqTjyYZVueCEGhrsjQUc5SXfzsnuQKxfl+U8DXVV2konzne2krKowmyhtIYpFNWkwUtXeHuoKPE/KsHyS9fv1FVihDolffb/103eYVN1Iy7ltew/Zu8vo1wLcGctwiolxrRtWTTwUuNGLLWHbEmYWYSjs78lb+pMqe8NxorWJNQwsBUil5Xt19AL9NvvodbEzelDRUoPh1dU9bfpGbX39B2wyEaivM+679Rq7NxwABbQzGbvNRn8GCkUvR7ZmsohvujhJ5aIcN9D9nk1SAoP0lQ5+cXKgkEdUjejtZp/vIRKLI+vmHD6kUIMKWYlvTQf5N3Hoq/UNDJ26G03rRhi13u6R825qnX260Jqk6jxfAShQn1tqaMdt5brvoNYRqKp+OO8DI7Rs+eEb2uZZQXRqtf79ZQiiGC6gPvmeibb1gAA+AvgNXv8dNQnGHS/h2pcSViYS3rYtqszm/ejZ62VfpQBg1Oeb+Bmbsuk0Y1eS8og3zBe89e3bwH972yxXGdLVDkc+vX57BTfoTTIUcERZ/RD+DtoB2TuXx9KF711vlC+Gkowdvv2OGb7odKF4gBW8tadXOeV//F0lACTDS60Ny8O3hfimzamsegt6faQrc5mw40O5puwbi2XKmOAU55P41D/a7d+4Ycx/V6hG3HrJej0J8dU8dv8r0pBrSFDoM0FH2nyKmjm/FXx03DvIPGeq6z08oycz6TfK16y/vsOj/o8KF40zoGX8a9VAqyzkXDJFDrU0Vb+j3Dfhp5kIYS9Nw6JyNWOA/FR9D/4rF1juvU/WyX/ZH+m75w0wv44YOrALCGMuIZ02p1POom3f/qVlz90GrHNZ7RgmmPAxf2hCwfzcfUjxWFtrVvyDPjN+HO/V0x6Hix3RqKKj+4UJXOtDCfyZyiLksR0DNg3jlP/Q6703B0jCWTl3twG7YYoj6CN4fgSoESMglUBWa40+nXeNPqAw7nOX2wYorGsjUUQ74D2m8K0lAcE0KLArl0ytZ8gzrQXMa88jZgvvd6WU2ZtFGwh70rqrMXwtseqh3V7ot6u5oWXdQpCoT4UGCX60YNyFJU+fL1eYffq/TdPSBQqyK3N1kbg+n38Km1O61Jz2S1BQuUEczUMc0ASi/+/yzbbJ+76q+OBODtgIPWmFLslh2Rc0nq8FDHFAXPpFX6S1Aos6mT3Kc5gD1RXnb5/uVa+XpH9X5RbFa+1oGOlix6BvKBHbtaf0z/HXulmawlZ9JQtHINdR0MCDnWjw0XRGC9+ocLjg5Jbzu/kNJ8werEgWCT11De29Ft22MFhJg6/v6hgt25Bpni3OHW6RQFBnuosnIZfw1F3Xu3VjVka5Epo29Av94zoVYIDAwX7fBed3uqy1uyzh0P9fro9Xem1aO8vPUKcsqrerZk0wms5eV/X3SUYFUaiklTbc9lkE7536NqwAKlDL5yyiz7xVcPVlpbCvTgSe0AvHMR/usva+zvfoJhY5dl4nH6LkJeAhk2DPgrKCqdZwZ3yMhZ3xfCsx+K8D68JlT99bIDfRny347mLApFYTQ9qLqWRqGlPC65dzkAc6ikPknM1JaD+aI9sc08/yF4hKun2a1pKWFOdcAaaavf434+gkxeL28sLWRouhX9wwWMabVGsJ4oL1lOSzbt0VCy6ZS9wm3QqgLZtHkZIj1/92hfldXRkkVReJ9LXevwC7pQVgK3r0M9p822QNE75OABQ1FYK13k0qnAeSjGe2/77jJGU1w5+Jm8tuzud9xDjw9FG/Ao2pszSKf8NeNqwAKlDOZMbrdNAQWDQGnLKQeZ80F/cX3pxTfd3D++XHK46e9Q3qcjto/JsGEC+T7E+YL3QbPSKtMRGX0o+6RAac2lY8+UV/XXOwmHX8jH5NXRYrWjyeyVt00mckRvqHuzwRSjysqmzPMfhvJFe1/xQZMN3aE5mu5F6Tfu0MLI/ToInXxBoClrXmVhSJvF7jbxqEEIYB5Q9A8XMabF6nzdo25VTksu7TEz6RqKX+cL+Ju8lHPcqoOzLVVZKjrJI3ACAljUtR0tZiFZ0l7lwMARGl0aUPiZvFJEaMqmjCYvld70mtn1cmkLUdGtASanvOLZdbsAWAOJb9/zGgAtyku2my5Ec5mU78oe1YIFSgROnTMBAPDZeTPsF02NJNOavUnZYINCRU1Ou7U7eu3v+gOlazomzUYIIVcb9q+7eijdnZn98mXSxhdACZSO5qzvJKqoC+X5BSW4k6t8x8oRaM+AN9LLbfLStafzjpkKAJg2tsV3HkombRa+i1Zuxx6pWfQZdu3z8znYdS8ITJMLCr6xvXQ/g+aC6Pk1ZbwmPKD0+7Jp8txD1U6Aj19ouIDRPhqK+gnNmZRn29lMiko+FIMMVM931scprw9e3M+OKquj2aqXn8CxynEWroTiGClQ3EJSj2yy6uHVSvx8CsWi9S61ZNNGgaLyNj3zylelOvdyzV5zL3vY/u4YgLh+vzJvLl61wz42od2KKuuVYe962UP5ItKpZMOaw2CBEoGrPn0Unvj2B0FE9gKQqrPXw3VtVdv1Frbm0vjSyTMBAD9Z9IbjXLEo8Msn3gQAnHzIeKfdVHug9pk6Od3k5fPMKDOB14GpjzL9fSiTOpqwbe+AMfom7DkdNmhHQRMI1bkxtkDxaijDWsfgzk8IYOb4VmtFaB8NJS1XOtBRneJQoQgiYJ9hccCwPckLRYHpYy2Bos9pyRcFJspQUlNHI4S1y2RJQ3H9Xnnf2poyjmgsoNRZ+y2T0j9cwNhWc+er2r45l/Z0vOmU/pwbAgGUhpI2d866kHALDPUsqNG822yld+bugY7Ka7SPhjJkt1Xa8b+eVzbtvxZfmgjNcQTKsNK6vA7yMIbyRWcovMuHcvYRB9j/d8sIwoy2HPWUMc0gAnb1DXnKThEhTda7YPJXVQMWKBGYNqYFB41vA1Aa/ahOR0UbLXjfDHuOitsJ3T9cQHuTeaeAXz75pv09l045o6uKpRfXNC+jINX0dMq86RdQekGHPRqKNRLNpL2dL1ASYIcd0AEhnB1DUKiqTmnplqLnGGDqPK0D42QnaBIoysSlRvTOyXsFNGfTyKRTvqv+NmXSgf6ktlzGqKE4fFs+83ZUR9erCZRCUZvfYJyY6Pw9fiavtlzGI5AGtMlrfk55ZfLyNQ9l0h7ndSZFtq/QJARVMEjaILgBp1BwD4RUfqrzdXfeetu5Bbe6tuQXctZN/Y7WnHduRr4gn3efOg9JH5qloRiCBQpKoHhOaaa4jOPaKLjfa/fExqZsCs9e8mEAzrZRNGXSGNOSRVffILr7hnDLs+vtc6kUcKDst3StuZqwQCkTNfpRN1dtzPSD845ARr6Eeoe1dyAPIUqjKjePr7bST2hvQlPGGSWiXpg5k9uNS3eomfLZtFnLUFEx7joB1oObThEyqZQx1LFLjnjUvhD6ix8Uk+8oQ+brUMO1sj59/bP46LVPaXWSJq82fw3FDg012MkH85ZzuzVrdZImh3ZHS8bTEekdTO9gHjc9/ban3L6Aji5fKGL9rj605tJIp8jW7gpFgYHhgjai9t4jO1qnyWzyGtZG3e4oL3VPWrNpb8ScENjRM4jWprTR0azKac2VRuTFojU7Pp1K2e1rEiiD+QJymRQy6ZSx81QBHZkUOebCqN+TTpEtZN2dt75emXuNOVtDUQJl2C1QZFvlTE55gUyajEKwWLS0xFwmheas2SmvBLJxvo/tILfqVY7J6+EV7zj+V5rErt5BbNi1D7l0ClPHtGBsaxa98n1wL68/ti2H7r5h/O1vX8INj5cGqCkizJrQCgDY3uNdWLQasEApE+V4X7/LWuxvr7zJTZm0/ZLoHaGa6DamNYd/+PAhAEov82NrdmDJ+i4AQDplmS7Ug7tvKI8fPWxFhx01fTQ27HIuLrhtTz+eeqMTKbJGXSYzzN/9/mW8I1eo9TjlC5ZA8RtldvYMIp0iO0zaGQprPfRBviIhBHrkgz/sYxdfs70HK7ft1c5Z+SnfQO+gVysrRXl5fSiD+YJ1H6Q2qI+OB1V0UXPWO7I1COOtu/vt70P5It7a2VeKAnP97jue34CisGzcrbm03Sl27xtCUQBTR7fY9XPzuxc2AiiZS/1WybVMXmbz0KjmrOfc3Us3AQCeXrsTQ4UifvXEW560TZkU2pszdgel2iGdAnJpf62qd7CA9qYsWrIp46zyF+UzPb495+mc+4byaMqk7N874GoTvbN0R0sq4aQEtDut7UNp8vpQhgsCmVQKKfKuFq6vgOxn8trbb9XL9HtLJi9VbnSfRYtr+231fKu5bVvkc9jenLEHsW6tZnxbDrv6BvHcW7scx9NEdlvt9Vl5ImlYoJSJWqLjH+961V4OQZkHsukUDuhoxtodJW3i5Y3WDnhjW7P2yFvd3C/d8qJ9XYoITZmSyeuHf15ln+toydqCS/HX//0C9g7kUZBhnqaH+IHXttnf3R3OvuECWrJpZNJktPm+sb0HE9pzxth8JYD6hvK+Me7b9w7aL7het+FC0Y7CcaNepjEBJi/1QrcZRvQDw5aJQI1Q9c5p0DZLZD0d+7f/8Kr9/YRZ4wAAO3tLe49cu9jye6nf0+WaDX//q1sBWC96e1PG1lBUHgeNt0aJvYbfEzYp9rrF1gxpo8lL/qZJHU2eAIblW/YAAHYY9lABLGHbmkujvSljT8YcGLLyb86mbeFpEoJ9g3m0N8nO13BeRSBlUt7R/qaufswY2zFEyjMAABw/SURBVGq/R+42cZsLdTbIqLbZEywzjrtz39xtnZ862hoEubeETvuYvNQ7l5OCzi1QBoYL2NM/DCKgd8g7P0q1kYo+U+9ToShw+/MbHFFcbrp9VlZ4Re6caQuUpqxHoPzbuYcDsFZ5UBYFnVSKbNOi30ThpGGBUibqRQOAb975CgDggvcfZB87ZsYY2zw1MFzAt+6yrhnblrM7yt2G0UKKCDkt4mamtH0C1gSloXzR0fGr5dC79g35CgWdbtcDt6d/GKNbs8ilU56OSgiBJ9d24tQ5E+3O32zyMq+ACwBfvGWJ/d09UlSRKW6U9qKc8m4hmi8U8e17rM5fjbw8NmeHhlJKr4erujuih1dsBwBc8cm5+KePvMcqu7+UVrW1ihBfu91pftQ7qNZc2vbB7Oq12vw9k0cBgPGlV6jf7G5PVfa4tpznPvVLATBxVJOjvgAwsd3qVL966ixjeX1DebTmMpg9oQ0bu/ZhYLiAXll2e1PGDjB4q9O77H7fYB5tTRk0Z8yjeSWUP3HMVAwMFx2Djs7eQUzqaMIBstN/x7XHixoEpFPeZ3rbngEQAbMmWPO93MJOzf+ZPtYS4O510Joy1goAbn+jyqfJx4eiBgazJrRBCK9f6E9y4NbhivJ6et1OfP9/X8fVDzoHDTq3PLPe8f9tz23AGT95wr6fSnaNasrYg4buvmHMnzUOXz7Furfj2pqMPpIpo5sxri2HjuaMPbCtNixQykS9aDodmn/koPGt2NzVj0JR4JFV2x3plIPUPSoB1DIJaY+J4aXvf8QezZm2Ju0bzBt9KG7NYemG0lyYTV378MBr29DelMGY1qxjIh4A/NPdr2JguIhDDxhVMk1oHceQy0fkZtW2vbZQHd2S9fgyxvsIFJVvLp1Ce1PGM3pds73HzndMizeipn8oj5ZcWtNQrDq/vbMPz71pmQM6mjO+4buZFNltrQsjZRs/+ZAJSJHVfjqqqQ8ab4281X267bn1AEoCZVeAQFGdsP6b9dH6hPYmg1O+gGyaMK4159FQlPnl8ycehAvef5A9mFHs7R/GqOYMJnVYHfue/mH8Utrf25oydp2fXNvpqWuvFCgtubRHA1m7vQcvvG2ZvEy+o129g5jQ3mSfc9d7zfYetOXS6GjOeDSJnb2DGNeas7VT98Bg0UrrfVPOccezMVxAi/Rxud8NNXof1Wz9JrdQV1qe0ozc76EqV7Wleo63Se3Cb6HTfUN5WwPR67JuR689ePnFXx8HwPIbqfd0y+5+TNf2vJ8jJ1TrfPq46bh2wbHIplM49sCxeGO7//YJScICpUxacxl88pipGK8tJqg/MAeOb8VQoYh39g5gS3fpYZk0qtl2Jj70+jvY7hqZpVPWpCrVaajOsKM5Y0eIGcNoC5bJa/e+YYcqvrPPa+pQHdTnb3oBAPDa5j0Y397kuDZfKNoTLWdNaLM1Ml0g9Q8V7MUUTbbZczRH+wffO9G16m3RsyKr6jjUqKyjJYNRzRlPZ6OPDFX0imr7O57fgPW79mFca9Y206mO4UM/fhxPrd0JAGjJZhwjUF1Q7u0fRosURjc+afkcfvvCBtwl/RE/+9yxmDK6xTGhEIAdFnz5eXPRmktjn7x3qm7vOaAdo5oyeGx1af7AYL6As6550v7/5EMmOOr8x5c32/MTvnjSTLTkUhjIFx2dbP+QFdU2qjmLvqGCowNVnV5rzuok3aPqzt4hTBzVZNf9rc4+3P78BgDWirkqHP3el5yr3ALWc9TelEFzNuUZzau5EkDJP1Baf0tgZ++gZUrNpi0TkvZMF4sC9760BUdMHS0j9Uq/tatvCCu27sXEUU22GWev6/lYLNtXCStdAPcPWSbenCGQQL1XHc1ZjGnN2vOR7LbqKWkoAGzfoPpNClWuiqj8/YvWc9OSy+CRldtxytWP4tVNJU1BTXj+q2Onwc3WPQO46JRZOHL6aACWGW/djl70DuaxbU+/Y2B7zpEHeNL/69nvtd/RqWOaHfelmrBAicEH3jvRMdo8ctpo+/tB46yHbsOuPtuxu/LyswCUnM03PvkWPnad1ekq/0uaCK1Za9bynv5h9A4OSx9Hyh41u18gwDK3bdzVh2ff3IXr5Qizq28I869cDAC47cvz8YNPHAHAPFFwfFvONs0Azg7h0Ckd9nyP25+zOpsdPQPY1TdkT+Lb3O0cYel84cSDMKG9CV19QxBCYPnmPdjZO2jbuBXK7KGEk9JqFmsd8HChiK/ettT+/4ipHQCAO5dswp59w/je/74OwDIdqRHsPsMOh03ZlEND0QcDi1Zutx3oSzd047E1O/DdP75unx/TmsPBk5wRd919Q3h09Q6cOmcCJnc0o70pYwuF9qYMDpvSgaZMGj2DeSzd0I110r+2elsP1shR46eOnaZFD1p1u+P5jXYZF3/oEMwY24pCUdiDlF8/9RZ+8+x6NGVSmCY7F6U59Q8V8F9yvlMuk0Jr1jKZ6sJoZ88gJrY34SjZYb3ZWTKZHC7b1sQfX96MFVv3YmxrzvI3DDlNoRfcbJk6P3zoJI9Aed+Vj2BguIjJHc1IpQjtuYyjc1Yhr4OFIsa2ZrGxq2Ru+/ody/Dqpt2Y0N6EjuYM0ilymBD1yLzxbZaQdAbHWP6tVoNwVc/dqOaMpe0N5h3PyOpt1n06eGK7nZdiq3xf/u3cw0t70ueL6O4bsoXH75dsxFduW4rN3f34sVyGSQiBC2VbffZ9M7D6irPxvpljHfU6oKP0nhx30FjkiwK3PrseRVEy6wHAlNFeq8kkLe30sa3Y2Tto3EI6aRpaoBDR2US0hojWEdEltSr3lEMmOv5XM7SBkgP2hsffxJbd/Tj0gFH2iFlFTAHATtmJf/5Ey/9CBJwq96BYvGo7egcLtiAZ22aNfD523dN4cX2XbY5ZeNpsXPbxI2yb/T1yocqXNG3iuIPG2iOVB+VWxepF/91XTsD4tibs6R+2R7cqeu375x6OaWNacNT0MQCATx1njaKUoBrXlkM2TVgqI3oA4NVNuzHzkgcAAB85fDKu+ORcjG7JYt9QAW/t7MPHf/40ACsm/vRDJ9mmKdURKkey5U+wOgwVMfTDP6+yX+QH/uEUTJYvzAPLt+HSP75m12FcW86jocyeWPJHqXkXNzz+Jrbu7ncIlMvPm4uWXBoLT5uNXCblCJpQKB/Zp65/Bm/v7MOxVywCUDJjtuYydof11NqdmNCec6T/29++BAAODbU5m7a1UDVi1xc/mDiqCTPl6HiD7GT/3wNW0EZ7U8Z+rpRg/pf/eQ06KvpQPTdCCGzZ3Y8Jo5oweVQzsmnC8s1W279ncrutAXxBPpu6+e8f77J8WLv3DaE5a5m81ChdDzC45JxD0az2aB/KY+/AsP3MKx9Ge3MGtzyzHl+/fRmKRYEr/rQSAHDqIRPw/tnj8crG3XbeS6QZrW8oDyLC2NacHUzwVmevnfb4g8ZivGzzXVLzfvKNTixZ34W500Zb2po20Fi7vccWgh0tWYyTaXWhsXzLHrTm0pgnO/z12nbOK+Qze+yBY2ytYcOuPrtsN0rI6YOSqaNb0JxNe/xJUx1mLcsEqSI/pxlM7wBw18IT8cvPH+84ds7cA7DgfTM8UXPVoGEFChGlAfwCwDkADgfwOSI6vBZlTxzVhEvPOVSvi/1dPQRPrd2JR1btcDwkKtRV58jp1mjwhNnjMVeODP/p7lfx+yUbbe3gkIklG+n5v3zOHo0fMqkd6RRhyXdPB2Cp5E+v3YlN3SWTTHtTBicdPB4AcPmfVmLjrn3Y2TuIz82fgZMOmYADRlujuZ89akUTvSY7lk8cXRKSHc0Z3P7cBvulBoCrP30Ujpg6Gs+/3YU9/cMYGC44zGLfPuu9AEpzSq5bvNY+9+FDJ+GmL74PD/zDqQAsn8srm3bbJpe2XBoLT5tl/14AWLGlFF48fUyr7UMBgD8vL8Xynzh7vN057+kfxovruxyOZSUwr35oNU666lG74/jKKbPskfmxM8Z4/BVTpFalTEQvb9yND/34cfv8l06y6tvWlMbbO/vwxBudGMwXba1U+Uje2N6L83/5rMOm/cWTZqIlm0Y2Tbjj+Q0QQtja6Hc+aj1nypyiL7sBAOccOcUOcnhk5Q70DebxqPTdfW7+gQBK/pRtewawdXc/vnb7MgCWlpJKWWvBKbOe7tw9TQ5wnl67E8+/tcth79/RM4jZE9qQLwo7IkmZCQHLb6QGDL2DBTyxpuSLOUWa99SA6aEV7+D6x0v7fXzrjDmYMa4VfUMFvNnZ59hj5l/PttpjfFvOFsrXas/Wb79yAtqbMsilU7YVQQmMQw8YhfamDJas77JNgj99pJT2gNHNGCfvl9J+3ursxSOrtuNDh07CQePbkEkR1kltbmC4gIWyLd8zeRQO6GjGqOYMVm7da2v9bt/VK5t2447nN9hm4fOPn44Z46w+4/RDJzuu1YXGYVNG2VYBALYpTHHHRSfg5399LE6YPR5nz3WawGZPbMdVnz7Kkb5amKdvNwbzAawTQrwFAER0J4DzAKysReFfPXU2/uPB1fjI4c6HIJ0iTBvTYr98yrmpuPXL821VFwBOnTMRi/7xNOthTadw8MQ2vCk7QJXHpI5m/Nf5R+OfZXirsmtPlB3JpFHNOHH2ODy6egceXb0DB3Q0ozmbwsofnA0AGN9u2co7ewZx2o8eAwD74frE0dPwo4ffwHWL1zo6/Yman2NguIihQhGf/ZXVuS88bTYOGN2MA8e14v5Xt+LoH/zF8Rt/fcE8+3f/9fwD8f3/fR33vWKF1l78oYPx1dNmA7C0ueljW/CD/yvdsvmzxoGI8Nl5M/DDP1sjXqX1AJa2MbrVPEl02ffOwPj2Jtvp+m/3rXCcP3rGGHzp5JmOxTjVbzrvmJIde8a4Vke6o6aPxk0Xvk+2tTeg4ORDxuPTcm9vZYpQ91h1yr/9ygk45LsPArBs5y+u70YuncLqK862180aLliaw6xL/wwAOPeoKVh42sEAYAvJ3zy73mHKOf/46XaE2M3PvI2bn7FMPzPHt+I/5HYKxxxoaZlnaj4bAPj/pNCfO60DL8kooKf+5UP2+fkzxyFFsDtNnX//xOF4z+RR+PY9r+FT1z/rOHftgmMcbfHJXzxjn/vT35+CudJEPG1MC9bJdex+/BfLRHfJOYcik07ZZpwzfvKEnfZrH5iNE2dbg6M12y2T4b/fv8J+tpZ97ww7iETNvdHn35xz5BQUhMAjq3bg4z9/GqNbsnb007QxLehoLoX2n3PtU/jWGXNsgTPvoLHIykmGNzz+JqaObsYybQClhOPxB43FnS9uwp3Sf/LfF8zDX1a8g6+eNhvPvbkL37zzFXtAmE0TfnT+0XYe3/noYfjUsdNw63Prcd8rW3GwplkTEe5ceCK+ettSfP7Eg+y10BSnyPUG603DaigApgHYpP2/WR6rCakU4bV/PxPX/81xnnMPfutU+7vaI0XxgfdMxMUfsjqJaxccgwntTZgzeZRtf733Gyfb1953cen7p4+fjmcu+bC9xhJgmV8UyjQFWKaP9x7QYXdUAPCXb53m8PUcIiNDWnJp/POZ73HU8ZOaCQ+wbPg6Hz/KOv9FuT6ZmzM0IZtOkf17AUs7URARPnbUFEfau7/2fgCW2eubp89xnPvsvOlY/E8fsP//t3MPx4Gy82/Npe3osVSKbB+L4t6/PQn3XXwyjpo+Br/6gtMkAAAHTSgJEfcg4PLz5toC9nSt/oo7LjrB/v6xI52/Z4YcZWbSKdy58ERotwSHT3XeozMOcw5O5mr36wDN7/Q/L1mmzd999QTMnthur9elow8IlLlE53sfO8zWpv/w9ZPw8aOn4upPH+kQpqNbs/j1hfM8aVf84Cwcf9A4jGrO2js7Kq447whbOM+c0OpJq9+XQ6d46/X1D1jPihKCigntTfiXs0pWge9+9DAAloAFLC1Pjx68/LwjHOnv/duTMLoli/OPnwHACkJQwuSMwybjIfnOHqJFTClh8o0PHoy/OcEy/6ln4/v3rcD/vrIVB41vtX2kAPDRuaX7P2NcC46ePgbf/djhmDSqGeceNRXzZ45Day6N7370MLx2WSkdYPm7jp4xBj/81JF47tIP26bHUn6teOhbp9lm8pEIVbJ2fz0hovMBnCWE+Ir8/wsA5gsh/t513UIACwHgwAMPPH7Dhg01qd9rm3djU1e/p8MErGiWoULRHk2VQ99gHi+8vQsffM8kR2dUKAr891NvYd9gHotW7cB1C47BnMneF/belzZjfHsTTpszwWGqyxeKyKRT9mxzN8OFIt7Y3oPDXIKqWBS46qHVeKuzF3913HR85PDJ9vpmOg+9/g5SBJx5hFMdF0LguTd34fcvbrLMcAeXRlrFosDqd3rw/Fu70JpLY4E04ejkC0U8sHwbzjriAEd7FooCRSGwcuteTBvbYpz78qsn3sSvnnwLf/j6+22Hq57v75ZsxJbuflwqOy+dTV37cN8rW3DaeyY6hDlgRQWtfmcvVmzdi4Wnzna0l+KVTbsdYdmqzgPDBdz78hZMaMvh7LkHOO4RADz/1i4suPF5XHTKLHzvY4fZ5weGC1ixdS8uu/91EAg3fP44h+O2f6iAm55+SwqgHN43c6y9VFAYr27ajbamDG56+m0sPG22He0EWM/Fjp5B/P6FjfjcCQd6zCrb9vRLk2gRR08f7fg9/UMFPLm2E0P5InoH8zj3qCmOTrRYFHhrZy/GtTVhuFC0/WaKgeECVm3biwde24YvnTLLU/ZwoYh7lm3G3KmjHSai+17Zgufe3IWjZ4zByQdPwJQxzY5n9vE1O/DY6h34wHsnYlRzFu+bOc4+1zeYx91LN2Hl1r04ZFI7zj16qqPcQlHgT69txZTRLZg9sc343BWLwvhMjGSIaJkQwju6cF/XwALl/QD+XQhxlvz/UgAQQvyHX5p58+aJpUuX+p1mGIZhDEQVKI1s8noRwBwimkVEOQALANxf5zoxDMO8a2lYp7wQIk9EfwfgYQBpADcLIVaEJGMYhmGqRMMKFAAQQvwZwJ/rXQ+GYRimsU1eDMMwzAiCBQrDMAyTCCxQGIZhmERggcIwDMMkAgsUhmEYJhEadmJjHIioB8A7APYEXDY64PyBADb6nAtLG3Su0rRB9apmudVqq0rqVUlbVZo338NkyuW2Si5tUm31XiGEd+kNN0KId80fgKUAbgy5xvc8gM4K0lZSblha33pVudyqtFUl9aqkrfgecltxW5nPA1gadK36ezeavP6vgvNhGzMHpa2k3LC0QfWqZrnVaquw89Vqq0rz5nuYTLncVsmlrVZbGXm3mbyWigjr0VQrfbUYifUaiXUCuF7lMBLrBIzMeo3EOgHJ1StqPu82DeXGOqevFiOxXiOxTgDXqxxGYp2AkVmvkVgnILl6RcrnXaWhMAzDMNXj3aahMAzDMFXiXS9QiOhmItpBRK9rx44moueIaDkR/R8RdcjjWSK6VR5fpfZgkeceJ6I1RPSK/PNu7VedOuWI6BZ5/FUi+qCW5nh5fB0RXUfu3ZrqV68k22oGET0m78cKIvqmPD6OiBYR0Vr5OVZLc6lskzVEdJZ2PLH2SrheibRXuXUiovHy+l4i+rkrr7q1VUi96tVWHyGiZbJNlhHRh0dIWwXVK7H30CZKKNj+/AfgNADHAXhdO/YigA/I718GcIX8/tcA7pTfWwGsBzBT/v84gHl1qNPFAG6R3ycBWAYgJf9fAuD9AAjAgwDOGSH1SrKtpgA4Tn4fBeANAIcD+E8Al8jjlwC4Wn4/HMCrAJoAzALwJoB00u2VcL0Saa8YdWoDcAqArwP4uSuverZVUL3q1VbHApgqv88FsGWEtFVQvRJpK0f9ksysUf8AzISzk9yLkn9pBoCV8vvnYIXSZQCMlzdzXDVuThl1+gWAz2vXLQYwXz54q7XjnwPwq3rXqxpt5arffQA+AmANgCny2BQAa+T3SwFcql3/sHzZq9Jeldarmu0VViftui9C67jr3VZ+9RoJbSWPE4BdsAYHI6Kt3PWqVlu9601ePrwO4BPy+/mwOkoAuAdAH4BtsGaf/lgI0aWlu0Wqjt+vRK0ts06vAjiPiDJENAvA8fLcNACbtfSb5bGkKbdeisTbiohmwhqRvQBgshBiGwDIT6XOTwOwSUum2qVq7VVhvRSJtlfEOvlR77YKo95t9WkALwshBjGy2kqvlyLRtmKBYubLAC4momWw1MoheXw+gAKAqbDMEv9MRLPlub8RQhwJ4FT594Ua1elmWA/pUgA/BfAsgDys0YibaoT0lVsvoAptRUTtAP4HwLeEEHuDLjUcEwHH610vIOH2KqNOvlkYjtWyrYKoa1sR0REArgbwNXXIcFnN28pQL6AK7yELFANCiNVCiDOFEMcD+D0sezZg+VAeEkIMCyF2AHgGwDyZZov87AHwO1jCp+p1EkLkhRD/KIQ4RghxHoAxANbC6syna1lMB7A1yTrFrFfibUVEWVgv12+FEPfKw9uJaIo8PwXADnl8M5yakmqXxNsroXol2l5l1smPereVL/VsKyKaDuCPAC4QQqg+o+5t5VOvqvRZLFAMqGgHIkoB+B6AX8pTGwF8mCzaAJwIYLU060yQabIAzoVlCqp6nYioVdYFRPQRAHkhxEqp9vYQ0YlSlb0Alr01UcqtV9JtJX/bTQBWCSF+op26H8CF8vuFKP32+wEsIKImaYqbA2BJ0u2VVL2SbK8YdTIyAtrKL5+6tRURjQHwACw/2DPq4nq3lV+9qtZnJemQacQ/WKPqbQCGYY0mLgLwTVgO9zcAXIWS07kdwB8ArACwEsC35fE2WFFMr8lz10JG6NSgTjNhOeRWAXgEwEFaPvPkQ/ImgJ+rNPWsVxXa6hRYJoTXALwi/z4KK2hiMSytaDFk8IRM813ZJmugRdwk2V5J1SvJ9opZp/UAugD0ynt++AhpK0+96tlWsAZTfdq1rwCYVO+28qtXkm2l//FMeYZhGCYR2OTFMAzDJAILFIZhGCYRWKAwDMMwicAChWEYhkkEFigMwzBMIrBAYZgRAhF9nYguKOP6maSt/Mww9SZT7wowDGNNNBNC/DL8SoYZubBAYZiEkIv1PQRrsb5jYU32vADAYQB+Amti7E4AXxRCbCOix2GtcXYygPuJaBSAXiHEj4noGFirDrTCmhD3ZSFENxEdD2udtH0Anq7dr2OYcNjkxTDJ8l4ANwohjoK1tP/FAH4G4DPCWu/sZgBXatePEUJ8QAjxX658bgPwrzKf5QAuk8dvAfAPQoj3V/NHMEwcWENhmGTZJEprJt0B4DuwNjZaJFcHT8NavkZxlzsDIhoNS9A8IQ/dCuAPhuO3Azgn+Z/AMPFggcIwyeJey6gH+P/bu2OUBoMgDMPvB9rY23qCHEEPIiIhV9ImaTyFF7DLBYKdV9AU2k6Kf9NFRRlIivcppxi2+5hl2WHzw0Tx+YfeOdBfOhleeUm9rpLsw+MWWAOX+1qS87Gb4ltVtQXek9yM0j3wUlUfwDbJ9ajf9R9f+j8nFKnXKzBPsmL6+fWRaZ3vw7iyOmNaOLb5pc8cWCa5AN6AxagvgKckX6OvdDL8bVhqMl55PVfV7MhHkY7CKy9JUgsnFElSCycUSVILA0WS1MJAkSS1MFAkSS0MFElSCwNFktRiB9xH3+8CYy9SAAAAAElFTkSuQmCC\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()" ] }, @@ -222,9 +2241,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": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXd8XGeZ6P99pqj34irbckuxncR2HMdpQAqJaWtKcjGwxJBczLJhy92Fu8Dv7g0bNkDuLrBkIWEDySahBEKATRYSwGmkuyVO3FvcJNmSbPUy/f39cc4Zj6QZaSSNNEXP9/OZz4zfOe+r9xzPzHOeLsYYFEVRFCUZXOnegKIoipI9qNBQFEVRkkaFhqIoipI0KjQURVGUpFGhoSiKoiSNCg1FURQlaVRoKIqiKEmjQkNRFEVJGhUaiqIoStJ40r2BVFNTU2Pq6+vTvQ1FUZSsYvv27aeNMbUjHZdzQqO+vp5t27alexuKoihZhYgcS+Y4NU8piqIoSaNCQ1EURUkaFRqKoihK0qjQUBRFUZJGhYaiKIqSNCMKDREpEJEtIvKmiOwWkX+yx78qIo0issN+vDdmzpdF5JCI7BeRG2LGLxaRnfZ7d4uI2OP5IvILe3yziNTHzNkgIgftx4ZUnryiKIoyOpIJufUD1xhjekTEC7wkIk/Z733HGPOvsQeLyBJgPbAUmAU8LSLnGGPCwL3ARuA14ElgLfAUcCvQboxZJCLrgbuAj4pIFXA7sAowwHYRecIY0z6+01YURVHGwoiahrHosf/ptR/D9YhdB/zcGOM3xhwBDgGrRWQmUGaMedVYPWYfBj4YM+ch+/VjwLW2FnIDsMkY02YLik1YgkZRFCUrMcbwq+0N9PpD6d7KmEjKpyEibhHZAbRg/Yhvtt/6vIi8JSIPiEilPTYbOBEzvcEem22/Hjw+YI4xJgR0AtXDrDV4fxtFZJuIbGttbU3mlBRFUdJCQ3s/f//LN3lse8PIB2cgSQkNY0zYGLMcqMPSGpZhmZoWAsuBk8C37MMl3hLDjI91Tuz+7jPGrDLGrKqtHTELXlEUJW30B8MA7G/uTvNOxsaooqeMMR3A88BaY0yzLUwiwA+B1fZhDcCcmGl1QJM9XhdnfMAcEfEA5UDbMGspiqJkJYFQBICDuSo0RKRWRCrs14XAdcA+20fh8CFgl/36CWC9HRE1H1gMbDHGnAS6RWSN7a+4GXg8Zo4TGXUj8Kzt9/gDcL2IVNrmr+vtMUVRlKzEbwuN/ae6sX7msotkoqdmAg+JiBtLyDxqjPmtiPxYRJZjmYuOAp8FMMbsFpFHgT1ACLjNjpwC+BzwIFCIFTXlRGHdD/xYRA5haRjr7bXaRORrwFb7uDuMMW3jOF9FUZS04mgaXb4QzV1+ZpQXpHlHo2NEoWGMeQtYEWf8k8PMuRO4M874NmBZnHEfcFOCtR4AHhhpn4qiKNlAIByJvj7Q3J11QkMzwhVFUSYRR9MAS2hkGyo0FEVRJhEVGoqiKErSBMKWi3dWeQH7m3tGODrzUKGhKIoyiTiaxuLppTS296V5N6NHhYaiKMok4giNquI8+gLhEY7OPFRoKIqiTCJOnkZ5oZf+YDjrcjVUaCiKokwiTshtRZEXY8AXjIwwI7NQoaEoijKJBGI0DYC+QHZVu1WhoSiKMokEQhE8LqE438qtzja/hgoNRVGUSSQQipDncVHodQNnq95mCyo0FEVRJpFA2BIaRXm20FBNQ1EURUlEIBQhz+2i0BYaap5SFEVREuKYp4ryLJ9Gf1Ad4YqiKEoC/IPMU6ppKIqiKAmJmqe8KjQURVGUEQiEIuSrI1xRFEVJBn8oPMCnoZqGoiiKkhBL03CT77F+fvs1I1xRFEVJhJOn4XIJhV537iX3iUiBiGwRkTdFZLeI/JM9XiUim0TkoP1cGTPnyyJySET2i8gNMeMXi8hO+727RUTs8XwR+YU9vllE6mPmbLD/xkER2ZDKk1cURZlsHEc4QFGeOyfNU37gGmPMRcByYK2IrAG+BDxjjFkMPGP/GxFZAqwHlgJrgXtExG2vdS+wEVhsP9ba47cC7caYRcB3gLvstaqA24FLgdXA7bHCSVEUJdtw8jQACvPcuecINxZOT0Kv/TDAOuAhe/wh4IP263XAz40xfmPMEeAQsFpEZgJlxphXjVVA/uFBc5y1HgOutbWQG4BNxpg2Y0w7sImzgkZRFCXriBUauappICJuEdkBtGD9iG8GphtjTgLYz9Psw2cDJ2KmN9hjs+3Xg8cHzDHGhIBOoHqYtRRFUbISx6cBUJjnoS/XfBoAxpiwMWY5UIelNSwb5nCJt8Qw42Odc/YPimwUkW0isq21tXWYrSmKoqQXf6xPw+vO7egpY0wH8DyWiajZNjlhP7fYhzUAc2Km1QFN9nhdnPEBc0TEA5QDbcOsNXhf9xljVhljVtXW1o7mlBRFUSYVJ7kPctQ8JSK1IlJhvy4ErgP2AU8ATjTTBuBx+/UTwHo7Imo+lsN7i23C6haRNba/4uZBc5y1bgSetf0efwCuF5FK2wF+vT2mKIqSdRhjBpinCrLQEe5J4piZwEN2BJQLeNQY81sReRV4VERuBY4DNwEYY3aLyKPAHiAE3GaMca7K54AHgULgKfsBcD/wYxE5hKVhrLfXahORrwFb7ePuMMa0jeeEFUVR0kUoYjCGAeapbNM0RhQaxpi3gBVxxs8A1yaYcydwZ5zxbcAQf4gxxoctdOK89wDwwEj7VBRFyXSc/uCx0VM5l9ynKIqipIbBQqMwz5N15ikVGoqiKJNEIDxU0wiEI4Ts8WxAhYaiKMokEdU0YsqIAFmVq6FCQ1EUZZLwDzFPZV9PDRUaiqIok4SjacTmaUB29dRQoaEoijJJDPZpnG35mj1Z4So0FEVRJomzPg1LWBTa3fvUPKUoiqIMIV6eBpBVuRoqNBRFUSaJQNgSDkPNUyo0FEVRlEEkCrlV85SiKIoyhMEht0W2T0M1DUVRFGUIg0NuC/M0ekpRFEVJQKKQWzVPKYqiKEMY7NPI87jwuETLiCiKoihDGRxyC5aJSjUNRVEUZQjxhEaRCg1FURQlHoFwBBHwuCQ6VpTnUfOUoiiKMpRAKEKe24XIWaFR6HXTr9FTiqIoSiyHWnrwhyIDTFNgmaeyKU9jxB7hiqIoyvg41NLNdd9+gbICTzRHw6Ewz023L4c0DRGZIyLPicheEdktIn9jj39VRBpFZIf9eG/MnC+LyCER2S8iN8SMXywiO+337hZbRxORfBH5hT2+WUTqY+ZsEJGD9mNDKk9eURRlMjhyug+ALl8oGm7rYJmnckvTCAF/b4x5XURKge0issl+7zvGmH+NPVhElgDrgaXALOBpETnHGBMG7gU2Aq8BTwJrgaeAW4F2Y8wiEVkP3AV8VESqgNuBVYCx//YTxpj28Z22oijK5NHc5QNAhPjmqWAOaRrGmJPGmNft193AXmD2MFPWAT83xviNMUeAQ8BqEZkJlBljXjXGGOBh4IMxcx6yXz8GXGtrITcAm4wxbbag2IQlaBRFUbKG5i4fLoGvfmAp65YP/PkszPPknKYRxTYbrQA2A1cAnxeRm4FtWNpIO5ZAeS1mWoM9FrRfDx7Hfj4BYIwJiUgnUB07HmeOogyLLxjG63bhjglvVJR00Nzlo6Yknw2X1w95L2fzNESkBPgV8LfGmC4sU9NCYDlwEviWc2ic6WaY8bHOid3bRhHZJiLbWltbhz0PZerwiR9t5p9/tyfd21AUTnX5mVFeEPc9yzwVxjLAZD5JCQ0R8WIJjJ8aY34NYIxpNsaEjTER4IfAavvwBmBOzPQ6oMker4szPmCOiHiAcqBtmLUGYIy5zxizyhizqra2NplTUnKc/kCYN46309jen+6tKAotXT6mlcYXGoV5bow5WzY900kmekqA+4G9xphvx4zPjDnsQ8Au+/UTwHo7Imo+sBjYYow5CXSLyBp7zZuBx2PmOJFRNwLP2n6PPwDXi0iliFQC19tjijIs+5u7iZizVUUVJZ2c6vIxozw/7ntFWda9LxmfxhXAJ4GdIrLDHvsK8DERWY5lLjoKfBbAGLNbRB4F9mBFXt1mR04BfA54ECjEipp6yh6/H/ixiBzC0jDW22u1icjXgK32cXcYY9rGdqrKVGJ3UydwttaPoqQLXzBMR1+Q6Qk0jbONmEJUFedN5tbGxIhCwxjzEvF9C08OM+dO4M4449uAZXHGfcBNCdZ6AHhgpH0qSix7mrqA7FH5ldylpcsPwPQEPo2CLGv5qmVElJxkz0lLaKimoaSb5m4rR2N6WQJNI8vMUyo0lJwjHDHsO9kNqNBQ0s+pTktozEgkNPJUaChKWjl6ppf+YBi3S/CHsuOLqOQuTjZ4IqHh9An3ZUl5dBUaSs7haBnnTi9VTUNJO81dPvI9LsoK47uQzzrCVWgoSlro6A8AMLuyUENulbTTbCf2xfbQiOWseSo76k+p0FByDicKpbzQiz+oQkNJL81dPqaVxs/RgLPmqX41TylKenBsw+WFXvyqaShp5nSPn5qSxEJDHeGKkmZ8wQhul1Cc7yEQimRNTR8lNznTGxhWaBR4VGgoSlrpD4Yp8LiiHdLUr6Gki2A4QkdfcFih4XIJBV5X1vQJV6Gh5Bz9wTCFee6zQkMjqJQ0cabHCsqoKR2+PEhRnkc1DUVJF75gmAKvO9ohTYWGki5O91glRIbTNCC7Wr6q0FByjqjQsHsxa/0pJV20Jik0ivLcGj2lKOmiPxCm0Osm36uahpJeTndbQqM2CaGh5ilFSRO+YIRCr5s8txWVoprGxNDWG+ALv3yTo6d7072VjOW07dOoLhnep1GYRS1fVWgoOUd/MEy+16U+jQnEFwzzPx/aymPbG3jhoLZYTsTpHj+FXjfF+cN3oSjK89AX1OgpRUkLvqBtnoqG3GbHHVw28e1NB3jjRAciZwvyKUM50+MfMXIKLEd4tpinkuncpyhZxeDoKS0lknoOtfSwdFYZZ3oCNNtNhpShnO4ZPrHPIZvMUyo0lJyj39Y0okJDk/tSTo8vRGm+F4/LpZrGMJzu8TOnqmjE49QRrihpxBeMaHLfBNPtD1FS4GF6WX60nakylJHqTjkUasitoqQPxxHuCA2Nnko9Pf4gpfkeppcVRNuZKgMJRwxtvQFqR4icAqv+VCAUIRLJ/DppIwoNEZkjIs+JyF4R2S0if2OPV4nIJhE5aD9Xxsz5sogcEpH9InJDzPjFIrLTfu9usQvMi0i+iPzCHt8sIvUxczbYf+OgiGxI5ckruUc4YgiEBobcqqaRenr9YYptodHRF8yarnOTSVtvgIiB6iQ0jQK7T7gvCzpNJqNphIC/N8acD6wBbhORJcCXgGeMMYuBZ+x/Y7+3HlgKrAXuERG3vda9wEZgsf1Ya4/fCrQbYxYB3wHusteqAm4HLgVWA7fHCidFGYzT3lWT+yaWHp9lnnL6RKiJaiAHmrv52ebjwMjZ4ACF9mfVlwVBGyMKDWPMSWPM6/brbmAvMBtYBzxkH/YQ8EH79Trg58YYvzHmCHAIWC0iM4EyY8yrxqpV/fCgOc5ajwHX2lrIDcAmY0ybMaYd2MRZQaMoQ3AiUAaWEcn8u7dswh8KEwhHKLE1DUBNVIP4yq938p2nD+ASWDiteMTjo5pGFmhso4qess1GK4DNwHRjzEmwBIuITLMPmw28FjOtwR4L2q8HjztzTthrhUSkE6iOHY8zJ3ZfG7E0GObOnTuaU1JyDMeZWKgFCyeMXr91jQcIDY2gGkCXL8jV59byb+tXUF7oHfH4bOrel7QjXERKgF8Bf2uM6Rru0DhjZpjxsc45O2DMfcaYVcaYVbW1tcNsTcl1HPW+QKOnJowen5W5bAkNy/SiuRoD6Q+GqSzKS0pgAOR7skfTSEpoiIgXS2D81Bjza3u42TY5YT+32OMNwJyY6XVAkz1eF2d8wBwR8QDlQNswaylKXJwvXYHHhcftwiUaPZVquv1BAEoKPJQXesnzuGhRTWMA/YEwBXnukQ+0KYj6NHJAaNi+hfuBvcaYb8e89QTgRDNtAB6PGV9vR0TNx3J4b7FNWd0issZe8+ZBc5y1bgSetf0efwCuF5FK2wF+vT2mKHGJmqfsL2yex6Wd+1KMo2mU5nsQEaaX5at5ahD9gTBF3uSFRmHUp5H5n9VkfBpXAJ8EdorIDnvsK8A3gUdF5FbgOHATgDFmt4g8CuzBiry6zRjjiM/PAQ8ChcBT9gMsofRjETmEpWGst9dqE5GvAVvt4+4wxrSN8VyVKYAvxqcBltqv5qnU0mu3JXWK8E0vLVDzVAzGGPrs7pHJklOOcGPMS8T3LQBcm2DOncCdcca3AcvijPuwhU6c9x4AHhhpn4oCA6OnwNI0NHoqtXQ7Po0CW2iUFbD35HBuzqmFPxTBGMYkNHLKEa7kDn2BEOEsyDwdC86XLio03C71aaSYHv9Z8xRAZbGX9r5AOreUUTg3LoU5ap5SoTHFMMZw3bf+xI9efDvdW5kQnIq2jmMx36tCI9U4Pg3HPFWc74mG4Spnb1yKpqojXMktevwhmjp97G/uTvdWJoT+QT6NPLdLfRopptcfQuTsj2JpvodAOKLX2aZvkIk0GfKzyKehQmOK0WL3LG7tzk3H5eDoqXyPCo1U0+0PUWJHTsFZjaPXnx2d5yYaX1TTSD53ulCFhpKpODWCclVonM3TOBs9pY7w1GL10jj7g+gIjR4VGsBZTWM0Pg2vW3CJ+jSUDKTFrhHUkqNCoz8YJs/jwuWy7oLzVNNIOT3+0ICe1yWOphFQoQFDtd1kEBEKvdnRU0OFxhTD0TDaegMEczDpzRcID7jD0+S+1NNjN2ByiGoaPhUaAP228ByNpgGWD0TNU0rGEWuWOt2Te9qGLxiJRqKA5dPQHuGppcf2aTiU5Luj48rYoqfAEhqqaSgZR6xZKhd7IDj9wR1U00g9Pb6BQuOsIzzzf/Amg6hPY9RCIztucFRoTDFaun3RkuG56AzvD4YHhDpqyG3qGaxpFOdp9FQs/WMWGmqeUjKQ1m4/588oBXLTGe4bJDQ0uS/1DPZplBZo9FQsY8kIBzVPKRlKS7efJbPK7Ne5V5nUN9g85daChanEGEOPP37IrWoaFv3BMF634HWP7ue1UDUNJdPwh8J09AWZVV5IVXFezpqnYs0CGnKbWvoCYYxhQMit1+0iz+OiR0NuAesajSYb3KHA69I8DSWzcIREbWk+tSX5OWqeGho9FQhHiORogcbJxjFBxZqnwMrV0JBbC18wPOrIKVCfhpKBOEJiWlk+08ryc1PTGHSXF+0TrhFUKSEqNPIHCo3ifLeap2z6BuUKJYsKDSXjcITEtNICaktzVGgM8mnkq9BIKV39VqvX0kGaRnGehx4NuQUcE2nydaccCrwufFlgSlWhMYWIahql+VGhYXXVzQ2MMXT1Bykv9EbHHKGRDfHv2YDTN6OyKG/AeEm+RzUNm/5AmELv6H9aC73uaORVJqNCYwpxxs4AryrOo7Ykn0A4Qld/7nzR+wJhQhFDWYzQUPNUamnrtTSN6uL8AePF+R6tPWXTHwyPqsKtQ4HXjS8UzvgbORUaU4iufispy+N2Rc0L3f5gmneVOrp81rmUFQwVGv4ssBVnA+29tqZR7B0wXpLv0TwNm7FHT7kxhozPKxpRaIjIAyLSIiK7Ysa+KiKNIrLDfrw35r0vi8ghEdkvIjfEjF8sIjvt9+4Wuxi/iOSLyC/s8c0iUh8zZ4OIHLQfG1J10lOVbl8w6sDMxdIPjtY00DxlfXlV00gNbX0BvG4Z4ghX89RZxhM9BZlvSk1G03gQWBtn/DvGmOX240kAEVkCrAeW2nPuERHn6t0LbAQW2w9nzVuBdmPMIuA7wF32WlXA7cClwGrgdhGpHPUZKlF6/KGohpGLPRA6bSdtWeHQst3dGg6aEtp7A1QW5UUbMDkUa8htlL5AaIzRU3bL1wzv/zKi0DDGvAC0JbneOuDnxhi/MeYIcAhYLSIzgTJjzKvGMtg9DHwwZs5D9uvHgGttLeQGYJMxps0Y0w5sIr7wUpKk23dWaJTkYBavE9kTa56qKrYctm22WUUZH229geg1jaUk301vIKz5MNiO8LFoGrZWnOnO8PH4ND4vIm/Z5itHA5gNnIg5psEem22/Hjw+YI4xJgR0AtXDrKWMkW5fkBL7BzUXi8w5Po1Y81Sl/QPXrkIjJbT3BYZETsFZzbVPfUdDqhIkizMn6zWNBNwLLASWAyeBb9njEudYM8z4WOcMQEQ2isg2EdnW2to63L6nNN3+oZpGbpqnYjQN+wfujAqNEQmEIvz3m01864/7+ci9r3Dx1zbR0jWwPlkiTUPrT1kEwxGCYUPReMxTOeDTGIIxptkYEzbGRIAfYvkcwNIG5sQcWgc02eN1ccYHzBERD1COZQ5LtFa8/dxnjFlljFlVW1s7llOaEnT7QpRFfRrWhzqXvuSOI7wsJvGsMM9NodetmkYS/PatJv7qkTf4/nOHaO7ycaY3wKHWngHHtPcFh0ROQW7ehIyFsbR6dXAc4TlpnrJ9FA4fApzIqieA9XZE1Hwsh/cWY8xJoFtE1tj+ipuBx2PmOJFRNwLP2n6PPwDXi0ilbf663h5Txkjc6KkM/4COhi5fkOI8N55B1UWrivNo61OhMRLHzvQhAju/egM/2rAKgPbesyHZ4Yihoy8Q1d5iUU3DwjfGXhpwVmhkunlqxAwUEXkEeBdQIyINWBFN7xKR5VjmoqPAZwGMMbtF5FFgDxACbjPGOFfgc1iRWIXAU/YD4H7gxyJyCEvDWG+v1SYiXwO22sfdYYxJ1iGvDCIYjuALRii1fRr5Hhcel+TUl7yzPzjANOVQVZynjvAkONnZT21JPsX5npgAgrOlZrr6g0TMWT9RLMXa8hWI6do3FvOU7Qj3ZfiN3IhCwxjzsTjD9w9z/J3AnXHGtwHL4oz7gJsSrPUA8MBIe1RGxgmHdHwaImJl8ebQl3xwCRGHyuI8NU8lwclOHzMrCoGzZULaYjQNR1uLHz2Ve3k/Y2Gs/cEh9x3hSpbh5CnEJmVZWbyZ/QEdDZ39wQHhtg5VRV41TyVBU0c/s8oLAKtHRlmBZ4CmEc0Gj2OeOuvTyJ0KA2PB0TTG2k8DctQRrmQfTrmQ0pgf1VwrZ93lCw1I7HOoKs4fYJtXhmKMoanDx8zywuiY5QuK0TR6R9Y0cukmZCz4oprGGGpPTYE8DSWLcDSN2MiiXCsy15XQp+Glxx/Cn+Fqfzrp7A/SHwwzq6IgOmb5gs5qGm3RulNDhYZz3TunuEbn3ISNxaeh5iklo4iapwoGm6dyTGjEMU+dTfBTbSMRTR1WPsasikGaRjyfRhzzVIHXTXGee8DxU5GoNlYy9BqNRL7HhUugL8O1NRUaU4SeeOapvNxxhIcjhm5/KK6mUa2lREbkZGc/ADPLE2sa7b0BCryuhOGklcV50X4bUxWnsVnNGISGiNjNrDL7O6lCY4rQPSh6CmzzVIbf1SSLEx0WN3qqSIXGSDR1WEIjVtOwos6C0f4Obb3BuFqGQ1Vx3pTPvG/p9lNR5I1WVx4t2RDRqEJjihA/esqd8Xc1yRItIVIQzxFuC40pfhc8HE2dPjwuoabkbHOl6uI8AuFI9DNyusdPxQhCY6qHNrd2+6ktyR/5wAQU57sz3s+oQmOK0O0Lked2DQgFdO5qMr1TWDLEK1boUKVFC0fkZEc/M8oLcLvOlnxzNLT23iD+UJjtx9q5sK484RpVRZpE2drjp7Z07EKjpMCb8RFoKjSmCN2+4ADTFFhCIxQxGd8pLBm64hQrdCgv9CKiRQuHo6nTx6yYcFuA6hKn2KOf195uo8cf4t1LpidcQzPvbU1jPEIj302PL7ODCVRoTBFie2k45FJPjc44vTQcPG4X5YVe1TSGobG9n5kx4bYQo2n0BXh6TzOFXjdXLKpJuEZlcR79wXDG5xlMFMYYWrv9TBuH0LCCUzL7+qnQmCL0+EMDwm0ht1q+Rs1TRUOFBmjRwuHo8gVp7Ohn8bSSAePVxdaP3+meAE/vbeYd59QMm+k81X1HPf4Q/cHwOM1TGj2lZAjdviCl+QN/UEtyqMickx9QEcc8BdZds2oa8dnV2AnAstkD/RVOCfTn9rVwstPHdecnNk2B+o6ccNvxmacyP+FWhcYUIZ556mx59Mz+kCZDc5eP0nxP9JwGU1bg0T7hCXCExgWDhEZJvoc8t4undp2ivNDL2mUzhl1nqrfWjQqNkoIRjkyM02s9k4NTVGhMEbp9ic1TuaBpnOr0Mb088Ze1tMBLd4Y7GNPFWw2dzK4opHpQqKiIRLWN/3nl/AGJofGY8kKjJzWaRqYHp6jQmCJ0+4aW2MglR/ipLh8zyhILjbJCD12qacRlV2PnEC3Doao4n/JCLxuuqB9xnaopnkSZKvMUZPZ3cvSlGJWswxhDj38Y81QGf0CTpbnLx8KFiSN7ygq8dPVb2c1W80gFrKizo2f6uGnVnLjvf+H6c3C7JG5U2mDKC724hClbSqSl24/XLQn9askQG5xSXTLCwWlChcYUoNsfImKGhqOW5OVGOetwxNDS7WdGeeI7vNICL6GIwReMjKkVZ66yO4E/w+HaEZzfsbhcQmXR1C0l0trtp6YkH5dr7DclTnBKdwb3JVHz1BSgrSd+HwSnRWe2axpnevyEI2ZE8xScDc1VLPae6gZgyayylKw3lbskjjexD7IjDF6FxhQgWtJ6UOVNj9tFvseV9ULjVJdV1nv6cELD1rKczHHFornLR57HFa0EPF6mclZ4i61pjIds8Gmo0JgCRDWNOMXminOgp8apTktozBg2esrRNLL7XFNNS5ePaaX5KfPzTOX6U00d/QOaWI2FkiyIaBxRaIjIAyLSIiK7YsaqRGSTiBy0nytj3vuyiBwSkf0ickPM+MUistN+726xP6Uiki8iv7DHN4tIfcycDfbfOCgiG1J10lON4dp05kLL12Zb0xjePGVrGmqeGkDLOMteDGaq9tTo8gXp7A9SV1k0rnWyIQw+GU3jQWDtoLEvAc8YYxYDz9j/RkSWAOuBpface0TE8TreC2wEFtsPZ81bgXZjzCKVcMNfAAAgAElEQVTgO8Bd9lpVwO3ApcBq4PZY4aQkT9Q8FU9o5Hmy3hF+qsuH2yVD8gxiUfNUfCyhMb6741iqi/No7wsSiWRuctpE0Nhu9SOZkyKhkck3ciMKDWPMC0DboOF1wEP264eAD8aM/9wY4zfGHAEOAatFZCZQZox51Vipjg8PmuOs9Rhwra2F3ABsMsa0GWPagU0MFV5KErT1Bsj3uCiKEzVUkgVNX0biVKd1t+weJmrF6bOhWeEDaenyMa0stZpGOGKmnEZ3oq0PgLrKwhGOHJ6cME8lYLox5iSA/TzNHp8NnIg5rsEem22/Hjw+YI4xJgR0AtXDrKWMkrbeANXFeXHt1sVZUOtmJJq7fMM6wUHNU/HwBcN0+UIjXrvRMFVb6zbYmsZ4hYbbJRR6M9tknGpHeLxbPTPM+FjnDPyjIhtFZJuIbGttbU1qo1OJtt4AlQmiY0pywRE+QjY4QL7HRZ7bRVf/0HPt8Yd4+NWj0fLqU4WWrvFnMA+mcgoLjaI8d1wT8GjJ9OCUsQqNZtvkhP3cYo83ALGppXVAkz1eF2d8wBwR8QDlWOawRGsNwRhznzFmlTFmVW1t7RhPKXdp6w0k/DDniiN8uMgpsOoolRZ4htSf2neqi/fd/SL/9/Hd3P/i2xO5zYyjpdsKIEilI3yqlhI50d5HXWVhSqLQrDbMmetnHKvQeAJwopk2AI/HjK+3I6LmYzm8t9gmrG4RWWP7K24eNMdZ60bgWdvv8QfgehGptB3g19tjyigZXmhkftOX4egPhOn2hZK6Wy4r9A4Juf3WHw/Q1R/kvBml/NeOpoyuLppqWuxaSal0hDu5QFMtgqqhvX/ckVMOJQWZ7WdMJuT2EeBV4FwRaRCRW4FvAu8WkYPAu+1/Y4zZDTwK7AF+D9xmjHF+kT4H/AjLOX4YeMoevx+oFpFDwN9hR2IZY9qArwFb7ccd9pgySoYTGk79/mz9sTw9isqiZQWeIdFTx870ckl9FbdcOZ/jbX28caJjQvaZibREkyJTr2lkWymR1m4/m/Y0s/ntM4THEPnV0N7HnHH6MxysiMbMFRoj1p4yxnwswVvXJjj+TuDOOOPbgGVxxn3ATQnWegB4YKQ9Konxh8L0+ENxE/vA0jSMgb5AOGEvikwmWo46iUzcweXRjTGcaOvnqsW1rF02g//zX7t4/I1GVs6dGpHdzd1+PHa9qFRRmOem0OvOulIi//y7PTy+w7J+37FuKTdfVp/03M6+IN2+UOo0jXwPJ+2E1UxEM8JznHa7o93gEiIO2RAXPhyjKUc9uDz66Z4A/cEwcyoLKSvwcu150/jD7uYJ22um0dJl1UoaT4G9eFilRLIrqOBUp4+ls8pYOquMn752fFSa94n21ITbOmR6RKMKjRzHcUgmqi2U7S1fRyU07PLoDs6XfU6VdYe4qr6KU12+6Jq5Tku3L6VOcIfKYi9tvdl1DTv6gsyuKOTP18xjf3M3rx9P3kzZMOhzNF6y3qehZDeO0EhkgijOy/yqmsPR2u1HJH62+2BKB7V8dRKynC/7UrvS6+6mzgnYaebR2u1nWgpzNByqivNp68suTaO9L0BlUR4fuGgWxXluHtlyPOm5bzV04nEJC2qLU7KXTA+DV6GR45yx7/iqE5insiEDdThO9/ipLMrD6x75o1xW4KU/GCZgt9JsGFT6YUlUaHRN0G4zi+auidE0qoqyS9MwxlhCoziPknwP61bM5ok3m6IhySOx9Wgby2aXU5SXGp9gSb4HXzAS/ZxmGio0cpz2kTSNHPBpJOMEh7OVbh1n+Im2PmpK8qNNmcoKvMytKpoSmkZzl4/2viALalPfHq6qOD/qS8sGegNhgmFDZZFVNWDjVQsIhSPc/+KREef6gmHePNHJ6vlVKduPI8iTFVqTjQqNHKetN4AIVIwkNDLY8TYcrT3JN75xSok4JqoT7X3MqRrovFw6q2xKaBqvH2sHYOXcipSvXVXspccfwh/KDpPn4Bur+ppiPnDRLH7y2jE6Rsg3efNEB4FwhEvqUyc0ZlZYn8lTGRpBpUIjx2nrC1BR6E1YzK8kCzqFDcdouqVFK93amsbxtr4hVUmXzirj2Jm+nK9Rtf1YO3keF0tnxW/zOh6cUiLZom102P6X2FI7f/muRfQGwvzq9cZh5245YqWOXVKfujDtmXZ1g0wNu1WhkeN09AWHjcPP5pavxhhO9/ipSeCvGcws+w5u/6luQuEITR2+OJqG9SO619Y2ntvXwpHTvSncdWbw+vF2LpxdTp4n9T8B2Va00Gkd4JinAM6dUUpdZSHbjg6fT7zlaBvnTi9NqMmPBackjmoaSlro6AtSHvNlGIwTPZWNjvAefwhfMJK0pnH+zFJmlhfw9N5mTnb6CEdMXE0DLGd4e2+ATz+4lbX/9gIPvHQkZ3pE+ENhdjV2sXLexCQxOjcp2VJKxDFBDf7hXzWvkm3H2ofN2XiroZOV81Jr4ivN91Cc56apsz+l66YKFRo5jhNKmAiXSyjKy86ihaPJ0QCraOF150/nhQOneXSbVXX/ojkDv/DTygqoKclnd1MXe05a2kZ9dTF3/HYPH/vha9EugdnMrsYuAuHIhGW+O5F62VJKpD1BZ8uL66to7fZHo+wG09lndeubX5OaUFsHEWFGeYFqGkp66OgLUjGMpgGZn4GaiKjQKEk+1+C6JdPpD4b53nOHuPa8aZw/s2zIMctml7G7qZM9tonqZ5+5lP9344W82dDBXU/tS83mJ4jvPn2Qmx/YQt8w/5877PpaE+EEhxhNI0uERltfEBEoLxz4PbnYFqrbjsU3UR2383zmVqVWaIBlSlWfhpIWOkbQNMBJJso+R/jpHutHaTT9INYsqKLErrf1N9ctjnvM0lllHGrpYceJDmaUFVBdks//WDWHD1w4i017mvEFM/NahSOGh189ygsHWvnbn+9IWHjv6Oleygo8E5LYB5aZJ9/j4tiZvglZP9V09AUoKxgaLHLujFJK8z1styPNBnNWaKQmEzyWGWWqaSgpxBjD5rfPDHs3CRAIRegNhKkoHEnTyFbzlPWlStYRDpDvcfPJy+bxsdVzubAu/p320lnlhCKGp/c2c/7M0uj4+y+aRbc/xAsHMrPR144THZzpDXDV4hr+uKeZ/34zbvsZGjv6mZ2i4nrxcLuEFXMr2DqCEzlTaO8Lxq0o4HYJy+dWsO3oCEKjOvXXcmZ5AS3dPkLhzEvwU6GRhbx86Awfve811nz9GX7y2rGEx0UdfCOU2Mj0UsyJaO7243WPvkrrP6w9j298+IKE7zvOcH8oEs0SB7h8YTWVRV5+t/Pk2DY8wTy9txmPS/jex1ZSVZzHiwdPxz2usb2f2RWpKa6XiNXzq9nd1Dmk6VUm0t4bSGjCvbCunIMtPXF/vI+39VJtZ5GnmhnlhUTM2Z4nmYQKjSxkl52xvGhaCXf8dk+0L8Jg2p348xF8GiX5mV0gLRGN7f3MLC9MeZXWOZVFlNo/BEtmns1j8LpdrF02g6f3NGdk4toze5tZPb+K8iIvl86v4rW3zwyJ/DHG0NjRn7KKrIm4dH4VEUNC004mMVywyNyqIsIRE9e/cLytL2VFCgeTybkaKjSykP2nuplRVsB3PrqcUDjCDxO0Ke3oG76EiENxlgqNpo5+ZlWk3i7vcgnn2xpGrHkK4LKFNfQGwrzdmlm5Gyfa+jjQ3MO1508H4LKF1TR29HOibWDkT1d/iB5/aMI1jRVzK/C4JJr8lskMl8vk9MiIF0F1vK1vQvwZkNm5Gio0spB9p7o5d0Yp86qLWbd8Nj957XjcRCpH0xgcFTKY4ix1hDd29DO7YmK+tJfUV1JdnMe86oGRMYvsWk2HWnom5O+OlTcbrIioNQuschaXLagG4LW3zww4rqHDssPPnmBNoyjPwwV15VkhNNp6Awm1cUcjc8qfOwTt5NB5E+DPAJhVHv/vZgIqNLKMUDjC4ZYezp1h3QF/+op6+oNhXjw41Dkb1TRG8GmUZKEjPBiO0NzlY/YEaBoAf3XNYp7626uGRNQsqC1GBA63ZpbQ2N3UhdctLJ5mfS4WTSuhpiSPVwcJjUb7jnmiNQ2A1fOreLOhg/5A5t6Q+IJh+oPhhN+RmeWFiAzVNJo6+q3k0AnSNMoKPVQUefnm7/dxy4NbM8ocqkIjyzh6ppdAOMK5060fh3OmlyJC3FIXyfo0ivM99AfDY+qNnC5OdfqImLOlQVJNgdfNtNKhAqnA66ausjDjNI09TV0smlYaLQsiIly+sIYnd57k7mcORstsN3bYQmOCNQ2AKxfVEAwbNh85M/LBaSJadyqBeSrP42JmWUG0YZeDEzk1b4KEhojwX395BbdcMZ9n97Xw9J6WCfk7Y0GFRpax71Q3QFTTKPC6mVVeyNE4QqOjP0Cex0Wh1z3smiVZWOm2aRJ//AazqLaEwxnm09hzsoslgxIV/7/3nc81503j25sO8LPNVpRdY3s/BV5Xwk6OqeSS+iryPS5eOBA/iisTONukLPGNVV1l0RBNw7lJm4hwW4f6mmK+8t7zmVVewC/sCgaZwLiEhogcFZGdIrJDRLbZY1UisklEDtrPlTHHf1lEDonIfhG5IWb8YnudQyJyt4iIPZ4vIr+wxzeLSP149pvNnOzs50u/eotNe5pxiWV+cJhfUxxX0+joDVJR6MW+nAlxyqP3+LJIaNh1eSZK0xiOhbUlvN3akzGaWUu31aJ26ayBQmN6WQH3/vnFzKkqZLPtW2jq7GdWReGIn4lUUOB1s3p+VVzTaaoxxvC/frGDP+4+Nap5p7qsz9H08sRmzrrKwqhZz+HFg6eZXVHIjAlKkHRwu4QbL67jxYOt0RuldJMKTeNqY8xyY8wq+99fAp4xxiwGnrH/jYgsAdYDS4G1wD0i4twC3wtsBBbbj7X2+K1AuzFmEfAd4K4U7Dcr+cHzh/n51hM8vqOJ+ppiCmK0B0doDA6vHKnulIPT9OVkhhZIi8dk2uYHs2haCf5QJGO+xE65kyWzhpZEAascxna78N5k5GjE8o7FtRxs6Znwz9a2Y+385o1GHk+Q0JiIxg4rOmm4a1JXWcjJzn6Cdq6GPxTm5UOnuea8aZMifG9aNQdj4NevN0z430qGiTBPrQMesl8/BHwwZvznxhi/MeYIcAhYLSIzgTJjzKvG+tV7eNAcZ63HgGtlMv6XMoxef4hfvd7INedN47rzp/HB5bMHvF9fU0yXLzQkgqqjf+S6U858gCOnMy9SIxGNHT6qi/MGCM/JwtHyMsWv4RRWjFdHC2DlvEpauv00dvRPSo5GLFedUwOQMNEwVTy61TLf7Ds5ugZaTR39eN0ybPfHusoiIgZO2gJm89tt9AXCXHPetLFveBTMqSpi2eyyIUEN6WK8QsMAfxSR7SKy0R6bbow5CWA/O1d2NhBrmGuwx2bbrwePD5hjjAkBnUD14E2IyEYR2SYi21pbM7PEw3j4rx2N9PhD3Hb1In604RL++tqBNZMW2D/6R88MNFElU3cKrGQ2lxDXL5JOwhHDQ68cZfuxtiFalFUKY/K1DLDMU5A5QmN3Uxd1lYUJQ6udaraPbDnO6Z5AyquyDse500upLc2fUKHR4w/xu50n8bqFI6d7R1UbrKmjnxnlBcMmiA4Ou312XwsFXheXLRzyUzRhXFhXwVsnOjOiPP94hcYVxpiVwHuA20TkHcMcG+9/xQwzPtycgQPG3GeMWWWMWVVbWzvSnrOOX2w9wZKZZQmrkjqawuCEs/YkKtyCFSFSV1nEkTOZJTSe3tvM7U/s5iP3vsrHfvjagG56TR390Vj2yaayOI+akrxoUEI6Mcaw/Wj7kBLvsZw3o5SiPDf3PH+Y0gIPN108Z9L2JyJctbiGlw62pvQH79XDZ/in/95NIBTh0a0n6AuE+dTl9UQMHGhO/v/lZIdvxM/R4AS/Z/e1cMXCmknVci+qK6fbH8qI7+i4hIYxpsl+bgF+A6wGmm2TE/azEyvWAMR+WuuAJnu8Ls74gDki4gHKgczPFkohvf4Quxo7uW7J9IT207rKQjwuGaBpGGPo6Ask3VGsvqY44zSNR7YcZ3pZPrd/YAnbj7Vz8/1b6PYFMcbY2eDpERoAF8+rzIhQ0uNtfZzq8rFmfuIe1R63i4vqKjAG/uKdC0fM20k171hcS3tfMGW9140x3P7ELv7z5aN86j+38I2n9nLV4ho+fuk8APadTF5oWAmiw3+OZlYUkOd2cbClm8aOfo639XHl4ppxncNocYprvmUncaaTMQsNESkWkVLnNXA9sAt4AthgH7YBeNx+/QSw3o6Imo/l8N5im7C6RWSN7a+4edAcZ60bgWfNcG20cpC3GjqJGFgxzJ2k1+1iTlXRgAiq3kCYYNiMmKPhML+6iKNxnOnp4kRbH3860MpHV83h01fM53sfX8mOEx38+LVjHD3TR18gPCCCbLJZs6CahvZ+TrSl1w+0+W3rHurSBcObSt69ZDoLaou55Yr5k7GtAVyxyPqBfSFFUVQvHzrDgeYe1iyo4pXDZ5hbVcT3Pr6SeVVFFHrd7D2VnHAKRwynunwj3nx43S6Wz6ngtbfbou1fL6lPLKQngsXTSijwunjzROek/t14jKc843TgN/bdrwf4mTHm9yKyFXhURG4FjgM3ARhjdovIo8AeIATcZoxxjI+fAx4ECoGn7AfA/cCPReQQloaxfhz7zUqchjnDmR/AiqCKtbE7P2Yzk7wbr68ppjcQprXHHzepbTKJRAz3PH8YgI+ungvADUtncP7MMl46eJqKQutO2SmZkQ4ce/bmI20TlhWcDK8dOUNVcR6LRxCgt1w5n09fUT8p0T6DqS3NZ8nMMl440MptVy8a93oPvHyEmpI8HrplNc/ubWHF3MqoP+fcGaVJaxot3VbL32Q01jULq/neswd5Zm8LJfkezptROuKcVOJxu1g2qzy7NQ1jzNvGmIvsx1JjzJ32+BljzLXGmMX2c1vMnDuNMQuNMecaY56KGd9mjFlmv/d5R5swxviMMTcZYxYZY1YbY+JX5sthdpxoZ151Udx6/7FcVFfBwZYeOu0MV8eue8705O7GHb/I0TRHUAVCET7z8DYe2XKcT1w6d4Dp4IqF1Ww72s5z+1uYVpo/qQ7dwZwzrZTKIi+vHk6viWrz222srq9KShikM/DwqnNq2H6sfdzd/Paf6ubZfS38+Zp55HvcvOeCmdHifmAVmNx7qispjdkJmU6m6OVlC6qJGPjdzpNWMUb35OdFX1hXwe6mrmh2f7rQjPAMxhjDG8c7hjVNOayeX4UxRBvfHGjuxuMSFtQkJzTmVztCI71+jZcOtfLMvhb+99pz+dq6ZQPeu2JxDYFwhKf3NnPZwuq0/gi6XMKl86uHFAScTBra+2js6OfSNGpcyfKRlXWEIoYHXzk6rnX+7ekDlOZ7+NTl9XHfXz6ngo6+II9sGTmDOpkcDYcVcyvI87gIRwyrJ9k05XDV4hr8oQgf/+Frac0RUqGRwZzs9NHS7Wd5EkJjxdwK8twutkSFRg/1NcXRWkQj4TjTxxqd0eULcnAUUSuJ2HKkHa9b+PTl84cIhdX1VXjdgjGWTyHdrFlQZTlG09TW9HdvWc2grlqc+RGD50wv5fol03nwlaNjbvi1u6mTp3ad4pYr5ycM8PjIyjreeU4t//fxXSMKdOeHNxkTboHXHe0ZvipNQuPq86bx3fXL2Xuyi/fe/SKb9jSnZR8qNDIEYww/+NNh7nvhMG81dLCrsZO7fr8PgOVzK0eYbX2oL5pTzmb7i3KguTta1DAZPG4X9TXF0ezi0XL747v50D2vRLNmx8qWI2e4YHY5hXlDwxmL8z2ssK9FJgiNq86xfqz/NAllMgYTiRh+tuU4q+dXpTUgYDR8/ppFdPYH+VGC/i/Dsauxk8/+eDvlhV5uuTKxM9/jdvHvH19BbWk+970w/N9p6uinvNCbdOe9a8+fRnmhN6mbuIli3fLZ/Pavr2J2RSGfeXgbuxon3zGuQiNDeH5/K998ah9ff3Iff/a9l3n/v7/EkztPcuuV87lwdvnICwCXzq9mV1MXp3v8HG/rY3GS/gyHd55Ty6uHz4y6RWdbb4DfvXWSHn+I/ePIXfAFw+xs7OSSYcJHP7pqDlefW0v9BBaKS5YFNcXMrihMS8/wlw+f5tiZPj5x6dxJ/9tj5cK6Cj5w0SzufubgqK7Z6R4//+M/XrWSPW9ZPWJ/mLICL1efN40tR9qG7bF95HRvtENeMnz6ivm8+A9Xx72hmUzm1xTzs8+sId/j4tE0FDJUoZEBBMMR/vl3e5hfU8zzX3gXP/jzlXz/4yt5/otX84/vX5J0O9NLF1RFs6iNYVSaBsDaZTMIhCM8v390P4KPbT9BwP5yvjmO6I43jncQDBsuHUZofOTiOv7z06vT6s9wEBHeYQva8WpYo+XhV49RVZzH2mUzJvXvjpdvfvgCzpleym0/e53HtjcklfDnlO34/idWJn2Xf/nCanr8IXYmuBNv7OjnlcNneOe5yZv23C6hrCC5EPaJprzQy/VLZ/D4jqZJ77WhQiMD+OW2Bg639vLl95xHfU0xa5fN5H0Xzhx1YblL6quYV13E9547BMDiUQqNlXMrqSnJ4w+jqBQaiRge2XKCVfMqqSzy8uaJsQuNLUfaEIGL52W+Y9fhnefU0uMP8fok9sJ+q6GDTXuaoxFE2URxvocfbVjFomklfOGXb7Lu+y/z3L6WYaOdth9rp8Dr4oIkNW44a758JUF0209eO4Yxhk+umTe6E8ggbry4js7+IM/undxeGyo00owxhvtfepsL68p595Lp41qrwOvm+x9fidflIs/tGrUJx+0S3r1kOs/ta0m6fs+rb5/hyOlePrFmLhfNqRhX8tHLh05z3oyyEc0PmcTli6pxu4QfvXQkJYEAyfAvf9hPZZGXz1w1+Yl6qaCusohf/cXlfOumi2jvC/DpB7fy+Z+9kTAc9/Xj7VxYV4F3FGGuNSX5nDejdEhItDGGA83dPLLlONcvmREtEZKNXLmohull+ZPea0OFRpp56dBpDrf28qnLU5N4tWx2Of9y04V89p0LxhRLfv2SGfQGwtHQ3ZH42ebjVBR5ec+ymSyfU8GBlu4xRcc0tPex5Wgb78kyc0tZgZdPrpnH03ubec93X2R308Q6JrcebePFg6e57epFlGaIqWQsuFzCRy6u47kvvIsv3nAuf9xzik8/uHWIucoXDLO7qTNadHE0XLawmq1H2wa0Mv6rR97g+u+8QI8vxGfekZ1C18HtEj62ei7P72+d1OKZKjTSzEOvHKWmJI/3XTgzZWuuWz6bv7/+3DHNXT2/Co9Lkkpaa+n28Yfdp7hxZZ0dvWXVN9rZMPofzsd3WOXGPrRi9ghHZh5f/bOlvPwP1+B2Cb/YOrF3fb/cdoLiPDefuDR7zSqxeN0ubrt6Ed/48IXsONHBb95oHPD+zsZOgmHDxfNGLzTed8FMguEIf/nT1wmEIuw71cVv3zrJJy6dywv/++qsMoMm4s/XzCPP4+KBl49M2t9UoTEJ9PpDfOPJvUO667V0+3hmXwvrL5mbMbbp4nwPF82pSKp2/09fO04oYviYHcFzkV1U7Sebj42qq50xhl+/3sAl9ZVpLckxHmZVFHKD7ZgcTWnu0eALhnlq1yluWDYj7RE8qebDK2Zz0ZwK7vr9vgGawXbbV5SowvNwrKqv4hsfvoA/HWjlrx95g/v+9DYFXhdfvOHctBa7TCU1Jfl8eMVsfrW9gWOTVAFXhcYk8NUndvMfL7zNrQ9tHVDe+4+7mzEG/mz5rDTubiiXL6zmrYbOYUNvd5zo4J7nD/GeZTOi/SWqivP4u3efw+/eOskXfvlm0sUPdzZ2cri1lw+tqBv54AzGcUw+M0GOyef3t9DtC7FuefZpYyPhcgn/533n09Lt5792WNrG03ua+Y8/Heac6SVUD9MkaTg+eslcvvqBJfx+9yl+/UYjH1lZl3Tl52xh4zsWkOdx8f67X+LxHY0jTxgnKjQmkG5fkB/86TC/3N7A2qUzOH6mjy88evbH9Pe7TrGgpnjEYnOTzWULqglHTEK/xom2Pm776etMKy3gGx++YMB7f33tYv7qmkX85o1GthxJzi/y69cbyXO7eN8FqTPRpYMrFtUws7yAe54/RF9goF8nEIoMuIMeC795o5GakjyumMTmP5PJqnmVnDu9lMe2N/DKodP8z4e3MbO8kHs+cfG41v3UFfP5+ocuYHZFIbcOkxiYrSyoLeHJv76Kc2aU8tPXjk94oyYVGhPEK4dPs+brz/DNp/Zx+cJq/v3jK/iHtefxxz3NPPFmE+29AV59+ww3LJuRETkHsaycV0mex8W/P3uIrz6xe4B29OrhM3zgey/R7QtyzydWxr1ru+3qRVQV542YkQtWjsp/v9lkZdsmWcY9U3G7hK+tW8bek118/mdvcLKzn3DE8Mrh01zzred5z3dfHHMJjTdPdPDHPc3ctGpOWorlTQYiwk2r6njjeAdffOwt6ioL+fVfXp6SjPePXzqXl790DQtqM+sGLVXMqSriFxvX8B+fvDjpvK6xMp7S6EoCdjd1svHh7cyqKOSbH7mQlXMrEBFuuXI+v9t5kn/67z1cfe40whGTkdFCBV437102gz/uaeaN4x2UFXr5u3efwxvH27nlwa3UVRbyw5tXRSvjxpu/4bJ6vvP0AQ40d3POMPkiLx5s5UxvgA+vzG7TlMN1S6bzT+uW8Y//tYvLvvEsbpcQjhhmVxRyoqOPrz+5l69/6IKRF7J55fBpfMEw333mEDUl+fzluxZO4O7TzwdXzOabT+2jsaOfez6xMi094LMVj9s1KQ22VGikmM6+IBsf3k5pgYeHblk9wOHmdgl3feRC/ux7L/Er2/E7moSlyeTf1q8AYOPD23jw5SNcsbCav/jJdmpL8/npZy4dsefGzZfN474XDidca5cAAAmFSURBVHPz/Vu4Y91S3p2g8+CvtjdSWeTlnedkftG9ZPnkmnmsmV/Fnw600t4XoK6yiD+7aBb/9vQBfvjiEa45dxrXjZCTEwhZVQIefvVYdOzb/+OirA6zTYaaknw+tGI2rT3+jLyhUkAypVNbqli1apXZtm1bWv62MYbP/eR1nt7bzK8+d3nCxkmdfUHyva6suIt680QH677/MgCzygt4ZOMa5lUn18firYYOvvjLt9jf3M286iLuWLdsgHDY09TF+/79RT77joV86T3nTcj+MwlfMMxNP3iVQy09PLJxTcKSGCfa+vj8I2/w5okONr5jAVefO42Wbh9/dtGsjDNlKrmDiGw3xqwa8TgVGsNzusfP1367hxVzKtgwTAJeS5ePr/xmJ0/vbeEr7z2Pje/IHTPC536ynaNn+njgU6uYWT66UMVAKMJv32rie88eossX4rkvvDN6t/zJ+zezs7GTP33x6qzKAh8Prd1+PnzvyzR3+nn/hTPZcHn9gJuLTXua+ftHd2AM/MtNF7J2WXYHByjZgwqNcRCOGH66+RhHT/fx5M6TNHf7MAauO38a/3LjRZzq8vH4jiYW1BQzv7aY42f6uOO3e/AFw/zvtedxS5raak4UkYhBZHyd33ac6OCD33+Zv3jnQv7Xuxfz3acPcs/zh/nH9y/JyYiW4TjZ2c9//OltfrntBL2BMHOriijwuuj1h2ns6GfZ7DK+//GVSWt0ipIKVGiMkVA4whcfe4vfvNFIcZ6bOVVF/OtNF7H1aBvfeHIf5UVeOvuC0aquDivmVvCvN10UzVlQhvJ3j+7g16834hKIGCuv4esfuiDpRlG5RrcvyGPbG9h2tJ1wxFCU72ZhbQm3Xjk/K0yXSm6RU0JDRNYC3wXcwI+MMd9MdOxYhUZLl4+bH9jC0TO9+IIRvnD9OXz+msUDjtnV2MnfP/omC6cVc8e6ZXT1B2lot8Iq33FOLe4JDnXLdrp8QX62+Thd/UEuXVCdU85vRcl2ckZoiIgbOAC8G2gAtgIfM8bsiXf8WIVGIBThL3+6nfrqYi5dUD3uirOKoijZRLJCIxtCblcDh4wxbwOIyM+BdUBcoTFW8jwufrThklQuqSiKknNkgzF5NhBbOrTBHlMURVEmmWwQGvEcBQNsaiKyUUS2ici21tbJ79esKIoyVcgGodEAzIn5dx3QFHuAMeY+Y8wqY8yq2lp1riqKokwU2SA0tgKLRWS+iOQB64En0rwnRVGUKUnGO8KNMSER+TzwB6yQ2weMMbvTvC1FUZQpScYLDQBjzJPAk+neh6IoylQnG8xTiqIoSoagQkNRFEVJmozPCB8tItIN7B/D1HKgM4VbSfV6NcDpFK6X6eebyvX02o0fvYbjI9Ov31wseTBy+KkxJqcewLYxzrsvxftI9XpjOq8sPt+UrafXTq9hBqyX6devNdlj1Tx1lv/O8PVSTaafbyZfv0w/10y+dg6Zfs6Zfg1Tvb+OZA/MRfPUNpNE0a1sI1fPazLQazd+9BqOj0y/fqPZXy5qGvelewMTRK6e12Sg12786DUcH5l+/ZLeX85pGoqiKMrEkYuahqIoijJBqNBIEyIyR0SeE5G9IrJbRP7GHq8SkU0ictB+rrTH3y0i20Vkp/18jT1eJCK/E5F99joJuxrmCqm6dvZ7vxeRN+11fmA3/cp5UnkNY9Z8QkR2Tfa5pIMUfwafF5H9IrLDfkxL13klRSrDtvQxqhC3mcBK+3UpVnfCJcD/A75kj38JuMt+vQKYZb9eBjTar4uAq+3XecCLwHvSfX7ZcO3sf5fZzwL8Clif7vPLtmtoj30Y+BmwK93nlm3XD3geWJXuc0r63NO9AX3Y/xHwOFZL2/3ATHtsJrA/zrECnAHy47z3XeAz6T6fbLt2gBcrjPGj6T6fbLuGQAnwkv2jOSWERoqvX1YJDTVPZQAiUo91J7IZmG6MOQlgP8dTVT8CvGGM8Q9apwL4APDMRO43k0jFtRORPwAtQDfw2ARvOeNIwTX8GvAtoG/CN5uBpOj7+5+2aeofRSRe47mMQYVGmhGREiyzyN8aY7qSOH4pcBfw2UHjHuAR4G5j91PPdVJ17YwxN2DdFeYDQ2z1ucx4r6GILAcWGWN+M6EbzVBS9Bn8hDHmAuAq+/HJidhrqlChkUZExIv1gfupMebX9nCziMy035+JdQfsHF8H/Aa42RhzeNBy9wEHjTH/NvE7Tz8pvnYYY3xYzb3WTfTeM4UUXcPLgItF5CiWieocEXl+cs4gvaTqM2iMabSfu7H8Qqsn5wzGhgqNNGGroPcDe40x34556wlgg/16A5at1DE9/Q74sjHm5UFr/TNWAbO/neh9ZwKpunYiUhLzBfcA7wX2TfwZpJ9UXUNjzL3GmFnGmHrgSuCAMeZdE38G6SWFn0GPiNTYr73A+4GMjkDT5L40ISJXYkU67QQi9vBXsOyij2JVnTwO3GSMaROR/wN8GTgYs8z1WBFTJ7B+7Bwb6feMMT+a8JNIEym8dgL8Fsss5QaeBf6XMSY0GeeRTlJ1DY0xsXfS9cBvjTHLJvwE0kwKP4O9wAtYgRhu4Gng74wx4ck4j7GgQkNRFEVJGjVPKYqiKEmjQkNRFEVJGhUaiqIoStKo0FAURVGSRoWGoiiKkjQqNBRlkhGRvxCRm0dxfP1UqR6rZD6edG9AUaYSIuIxxvwg3ftQlLGiQkNRRomdxPZ7rESuFVhlsW8Gzge+jVX19TTwKWPMSbusxivAFcATIlIK9Bhj/tWu3fQDrBL3h4FbjDHtInIx8ABWEcCXJu/sFGV41DylKGPjXOA+Y8yFQBdwG/DvwI3GGOcH/86Y4yuMMe80xnxr0DoPA/9gr7MTuN0e/0/gr40xl03kSSjKaFFNQ1H+//buGKWBKIrC8H9QGxsrW1fgFlyASxCRkB2JNtoIoo24ArFKlw0EOxsLSw2i7VjMC4Qh6suQIuD/lZfhMlMd3ivO9PMy1yF0S1shsQ88lmbrDeB17vm77oIkO7RhMiqja+B+wfwGOFz9J0jLMzSkfrr9Ox/A5JeTwecSu7Ngv7QWvJ6S+tlLMguII2AM7M5mSbbKvxN+1DTNFHhLclBGJ8CoaZp3YFpK8QCOV//6Uj+eNKR+noBBkkva5tJz4AE4K9dLm8ApMPljzwC4SLINPAPDMh8CV0m+yl5pLdhyKy3pP1WAS11eT0mSqnnSkCRV86QhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqp9AyVXdR8p0aV7AAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -259,10 +2301,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 18, + "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", @@ -281,7 +2321,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 19, "metadata": {}, "outputs": [], "source": [ @@ -305,9 +2345,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 20, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 20, + "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='*')" ] @@ -321,9 +2384,59 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 21, "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": 21, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -338,9 +2451,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 22, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 22, + "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)" ] -- 2.18.1