diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5b58a3ae26346460dd39e62a39c55243d7..6f45eb70456b53e831c1fd1a91cb895ddfaeaac7 100644 --- a/module3/exo1/analyse-syndrome-grippal.ipynb +++ b/module3/exo1/analyse-syndrome-grippal.ipynb @@ -9,7 +9,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -28,19 +28,19 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 17, + "metadata": {}, "outputs": [], "source": [ - "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" + "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"\n", + "data_local = \"C:\\\\Users\\\\lefevrej\\\\Documents\\\\training_repro\\\\local-file.csv\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ + "Les données sont disponibles soit à partir d'un [URL](http://www.sentiweb.fr/datasets/incidence-PAY-3.csv), soit à partir d'un fichier enregistré en local.\n", "Voici l'explication des colonnes données [sur le site d'origine](https://ns.sentiweb.fr/incidence/csv-schema-v1.json):\n", "\n", "| Nom de colonne | Libellé de colonne |\n", @@ -61,11 +61,981 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
02022333133379122.017552.02014.026.0FRFrance
120223232225716158.028356.03425.043.0FRFrance
220223132182816268.027388.03325.041.0FRFrance
320223031966314779.024547.03023.037.0FRFrance
420222932426818906.029630.03729.045.0FRFrance
520222832484519214.030476.03729.045.0FRFrance
620222734074533994.047496.06151.071.0FRFrance
720222633401028521.039499.05143.059.0FRFrance
820222532337719042.027712.03528.042.0FRFrance
920222432632821829.030827.04033.047.0FRFrance
1020222332343018950.027910.03528.042.0FRFrance
1120222231895115099.022803.02923.035.0FRFrance
1220222131363210251.017013.02116.026.0FRFrance
1320222031978715756.023818.03024.036.0FRFrance
1420221931788414079.021689.02721.033.0FRFrance
1520221833035325089.035617.04638.054.0FRFrance
1620221733600630373.041639.05446.062.0FRFrance
1720221634994942836.057062.07564.086.0FRFrance
18202215310080690824.0110788.0152137.0167.0FRFrance
192022143155441143891.0166991.0234217.0251.0FRFrance
202022133191914179558.0204270.0289270.0308.0FRFrance
212022123166224155035.0177413.0251234.0268.0FRFrance
222022113122849113306.0132392.0185171.0199.0FRFrance
2320221038790479741.096067.0133121.0145.0FRFrance
2420220935018243958.056406.07667.085.0FRFrance
2520220833096325942.035984.04739.055.0FRFrance
2620220733488229446.040318.05345.061.0FRFrance
2720220634662340398.052848.07061.079.0FRFrance
2820220536297056043.069897.09585.0105.0FRFrance
2920220437220964804.079614.010998.0120.0FRFrance
.................................
194319852132609619621.032571.04735.059.0FRFrance
194419852032789620885.034907.05138.064.0FRFrance
194519851934315432821.053487.07859.097.0FRFrance
194619851834055529935.051175.07455.093.0FRFrance
194719851733405324366.043740.06244.080.0FRFrance
194819851635036236451.064273.09166.0116.0FRFrance
194919851536388145538.082224.011683.0149.0FRFrance
19501985143134545114400.0154690.0244207.0281.0FRFrance
19511985133197206176080.0218332.0357319.0395.0FRFrance
19521985123245240223304.0267176.0445405.0485.0FRFrance
19531985113276205252399.0300011.0501458.0544.0FRFrance
19541985103353231326279.0380183.0640591.0689.0FRFrance
19551985093369895341109.0398681.0670618.0722.0FRFrance
19561985083389886359529.0420243.0707652.0762.0FRFrance
19571985073471852432599.0511105.0855784.0926.0FRFrance
19581985063565825518011.0613639.01026939.01113.0FRFrance
19591985053637302592795.0681809.011551074.01236.0FRFrance
19601985043424937390794.0459080.0770708.0832.0FRFrance
19611985033213901174689.0253113.0388317.0459.0FRFrance
196219850239758680949.0114223.0177147.0207.0FRFrance
196319850138548965918.0105060.0155120.0190.0FRFrance
196419845238483060602.0109058.0154110.0198.0FRFrance
1965198451310172680242.0123210.0185146.0224.0FRFrance
19661984503123680101401.0145959.0225184.0266.0FRFrance
1967198449310107381684.0120462.0184149.0219.0FRFrance
196819844837862060634.096606.0143110.0176.0FRFrance
196919844737202954274.089784.013199.0163.0FRFrance
197019844638733067686.0106974.0159123.0195.0FRFrance
19711984453135223101414.0169032.0246184.0308.0FRFrance
197219844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1973 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202233 3 13337 9122.0 17552.0 20 14.0 \n", + "1 202232 3 22257 16158.0 28356.0 34 25.0 \n", + "2 202231 3 21828 16268.0 27388.0 33 25.0 \n", + "3 202230 3 19663 14779.0 24547.0 30 23.0 \n", + "4 202229 3 24268 18906.0 29630.0 37 29.0 \n", + "5 202228 3 24845 19214.0 30476.0 37 29.0 \n", + "6 202227 3 40745 33994.0 47496.0 61 51.0 \n", + "7 202226 3 34010 28521.0 39499.0 51 43.0 \n", + "8 202225 3 23377 19042.0 27712.0 35 28.0 \n", + "9 202224 3 26328 21829.0 30827.0 40 33.0 \n", + "10 202223 3 23430 18950.0 27910.0 35 28.0 \n", + "11 202222 3 18951 15099.0 22803.0 29 23.0 \n", + "12 202221 3 13632 10251.0 17013.0 21 16.0 \n", + "13 202220 3 19787 15756.0 23818.0 30 24.0 \n", + "14 202219 3 17884 14079.0 21689.0 27 21.0 \n", + "15 202218 3 30353 25089.0 35617.0 46 38.0 \n", + "16 202217 3 36006 30373.0 41639.0 54 46.0 \n", + "17 202216 3 49949 42836.0 57062.0 75 64.0 \n", + "18 202215 3 100806 90824.0 110788.0 152 137.0 \n", + "19 202214 3 155441 143891.0 166991.0 234 217.0 \n", + "20 202213 3 191914 179558.0 204270.0 289 270.0 \n", + "21 202212 3 166224 155035.0 177413.0 251 234.0 \n", + "22 202211 3 122849 113306.0 132392.0 185 171.0 \n", + "23 202210 3 87904 79741.0 96067.0 133 121.0 \n", + "24 202209 3 50182 43958.0 56406.0 76 67.0 \n", + "25 202208 3 30963 25942.0 35984.0 47 39.0 \n", + "26 202207 3 34882 29446.0 40318.0 53 45.0 \n", + "27 202206 3 46623 40398.0 52848.0 70 61.0 \n", + "28 202205 3 62970 56043.0 69897.0 95 85.0 \n", + "29 202204 3 72209 64804.0 79614.0 109 98.0 \n", + "... ... ... ... ... ... ... ... \n", + "1943 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1944 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1945 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1946 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1947 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1948 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1949 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1950 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1951 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1952 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1953 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1954 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1955 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1956 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1957 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1958 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1959 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1960 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1961 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1962 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1963 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1964 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1965 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1966 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1967 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1968 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1969 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1970 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1971 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1972 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 26.0 FR France \n", + "1 43.0 FR France \n", + "2 41.0 FR France \n", + "3 37.0 FR France \n", + "4 45.0 FR France \n", + "5 45.0 FR France \n", + "6 71.0 FR France \n", + "7 59.0 FR France \n", + "8 42.0 FR France \n", + "9 47.0 FR France \n", + "10 42.0 FR France \n", + "11 35.0 FR France \n", + "12 26.0 FR France \n", + "13 36.0 FR France \n", + "14 33.0 FR France \n", + "15 54.0 FR France \n", + "16 62.0 FR France \n", + "17 86.0 FR France \n", + "18 167.0 FR France \n", + "19 251.0 FR France \n", + "20 308.0 FR France \n", + "21 268.0 FR France \n", + "22 199.0 FR France \n", + "23 145.0 FR France \n", + "24 85.0 FR France \n", + "25 55.0 FR France \n", + "26 61.0 FR France \n", + "27 79.0 FR France \n", + "28 105.0 FR France \n", + "29 120.0 FR France \n", + "... ... ... ... \n", + "1943 59.0 FR France \n", + "1944 64.0 FR France \n", + "1945 97.0 FR France \n", + "1946 93.0 FR France \n", + "1947 80.0 FR France \n", + "1948 116.0 FR France \n", + "1949 149.0 FR France \n", + "1950 281.0 FR France \n", + "1951 395.0 FR France \n", + "1952 485.0 FR France \n", + "1953 544.0 FR France \n", + "1954 689.0 FR France \n", + "1955 722.0 FR France \n", + "1956 762.0 FR France \n", + "1957 926.0 FR France \n", + "1958 1113.0 FR France \n", + "1959 1236.0 FR France \n", + "1960 832.0 FR France \n", + "1961 459.0 FR France \n", + "1962 207.0 FR France \n", + "1963 190.0 FR France \n", + "1964 198.0 FR France \n", + "1965 224.0 FR France \n", + "1966 266.0 FR France \n", + "1967 219.0 FR France \n", + "1968 176.0 FR France \n", + "1969 163.0 FR France \n", + "1970 195.0 FR France \n", + "1971 308.0 FR France \n", + "1972 213.0 FR France \n", + "\n", + "[1973 rows x 10 columns]" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "raw_data = pd.read_csv(data_url, skiprows=1)\n", + "try:\n", + " raw_data = pd.read_csv(data_local, encoding = 'iso-8859-1', skiprows=1)\n", + "except ValueError:\n", + " raw_data = pd.read_csv(data_url, encoding = 'iso-8859-1', skiprows=1)\n", "raw_data" ] }, @@ -78,9 +1048,73 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
173619891930NaNNaN0NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1736 198919 3 0 NaN NaN 0 NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1736 FR France " + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -94,9 +1128,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 5, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
02022333133379122.017552.02014.026.0FRFrance
120223232225716158.028356.03425.043.0FRFrance
220223132182816268.027388.03325.041.0FRFrance
320223031966314779.024547.03023.037.0FRFrance
420222932426818906.029630.03729.045.0FRFrance
520222832484519214.030476.03729.045.0FRFrance
620222734074533994.047496.06151.071.0FRFrance
720222633401028521.039499.05143.059.0FRFrance
820222532337719042.027712.03528.042.0FRFrance
920222432632821829.030827.04033.047.0FRFrance
1020222332343018950.027910.03528.042.0FRFrance
1120222231895115099.022803.02923.035.0FRFrance
1220222131363210251.017013.02116.026.0FRFrance
1320222031978715756.023818.03024.036.0FRFrance
1420221931788414079.021689.02721.033.0FRFrance
1520221833035325089.035617.04638.054.0FRFrance
1620221733600630373.041639.05446.062.0FRFrance
1720221634994942836.057062.07564.086.0FRFrance
18202215310080690824.0110788.0152137.0167.0FRFrance
192022143155441143891.0166991.0234217.0251.0FRFrance
202022133191914179558.0204270.0289270.0308.0FRFrance
212022123166224155035.0177413.0251234.0268.0FRFrance
222022113122849113306.0132392.0185171.0199.0FRFrance
2320221038790479741.096067.0133121.0145.0FRFrance
2420220935018243958.056406.07667.085.0FRFrance
2520220833096325942.035984.04739.055.0FRFrance
2620220733488229446.040318.05345.061.0FRFrance
2720220634662340398.052848.07061.079.0FRFrance
2820220536297056043.069897.09585.0105.0FRFrance
2920220437220964804.079614.010998.0120.0FRFrance
.................................
194319852132609619621.032571.04735.059.0FRFrance
194419852032789620885.034907.05138.064.0FRFrance
194519851934315432821.053487.07859.097.0FRFrance
194619851834055529935.051175.07455.093.0FRFrance
194719851733405324366.043740.06244.080.0FRFrance
194819851635036236451.064273.09166.0116.0FRFrance
194919851536388145538.082224.011683.0149.0FRFrance
19501985143134545114400.0154690.0244207.0281.0FRFrance
19511985133197206176080.0218332.0357319.0395.0FRFrance
19521985123245240223304.0267176.0445405.0485.0FRFrance
19531985113276205252399.0300011.0501458.0544.0FRFrance
19541985103353231326279.0380183.0640591.0689.0FRFrance
19551985093369895341109.0398681.0670618.0722.0FRFrance
19561985083389886359529.0420243.0707652.0762.0FRFrance
19571985073471852432599.0511105.0855784.0926.0FRFrance
19581985063565825518011.0613639.01026939.01113.0FRFrance
19591985053637302592795.0681809.011551074.01236.0FRFrance
19601985043424937390794.0459080.0770708.0832.0FRFrance
19611985033213901174689.0253113.0388317.0459.0FRFrance
196219850239758680949.0114223.0177147.0207.0FRFrance
196319850138548965918.0105060.0155120.0190.0FRFrance
196419845238483060602.0109058.0154110.0198.0FRFrance
1965198451310172680242.0123210.0185146.0224.0FRFrance
19661984503123680101401.0145959.0225184.0266.0FRFrance
1967198449310107381684.0120462.0184149.0219.0FRFrance
196819844837862060634.096606.0143110.0176.0FRFrance
196919844737202954274.089784.013199.0163.0FRFrance
197019844638733067686.0106974.0159123.0195.0FRFrance
19711984453135223101414.0169032.0246184.0308.0FRFrance
197219844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1972 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202233 3 13337 9122.0 17552.0 20 14.0 \n", + "1 202232 3 22257 16158.0 28356.0 34 25.0 \n", + "2 202231 3 21828 16268.0 27388.0 33 25.0 \n", + "3 202230 3 19663 14779.0 24547.0 30 23.0 \n", + "4 202229 3 24268 18906.0 29630.0 37 29.0 \n", + "5 202228 3 24845 19214.0 30476.0 37 29.0 \n", + "6 202227 3 40745 33994.0 47496.0 61 51.0 \n", + "7 202226 3 34010 28521.0 39499.0 51 43.0 \n", + "8 202225 3 23377 19042.0 27712.0 35 28.0 \n", + "9 202224 3 26328 21829.0 30827.0 40 33.0 \n", + "10 202223 3 23430 18950.0 27910.0 35 28.0 \n", + "11 202222 3 18951 15099.0 22803.0 29 23.0 \n", + "12 202221 3 13632 10251.0 17013.0 21 16.0 \n", + "13 202220 3 19787 15756.0 23818.0 30 24.0 \n", + "14 202219 3 17884 14079.0 21689.0 27 21.0 \n", + "15 202218 3 30353 25089.0 35617.0 46 38.0 \n", + "16 202217 3 36006 30373.0 41639.0 54 46.0 \n", + "17 202216 3 49949 42836.0 57062.0 75 64.0 \n", + "18 202215 3 100806 90824.0 110788.0 152 137.0 \n", + "19 202214 3 155441 143891.0 166991.0 234 217.0 \n", + "20 202213 3 191914 179558.0 204270.0 289 270.0 \n", + "21 202212 3 166224 155035.0 177413.0 251 234.0 \n", + "22 202211 3 122849 113306.0 132392.0 185 171.0 \n", + "23 202210 3 87904 79741.0 96067.0 133 121.0 \n", + "24 202209 3 50182 43958.0 56406.0 76 67.0 \n", + "25 202208 3 30963 25942.0 35984.0 47 39.0 \n", + "26 202207 3 34882 29446.0 40318.0 53 45.0 \n", + "27 202206 3 46623 40398.0 52848.0 70 61.0 \n", + "28 202205 3 62970 56043.0 69897.0 95 85.0 \n", + "29 202204 3 72209 64804.0 79614.0 109 98.0 \n", + "... ... ... ... ... ... ... ... \n", + "1943 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1944 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1945 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1946 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1947 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1948 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1949 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1950 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1951 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1952 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1953 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1954 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1955 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1956 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1957 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1958 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1959 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1960 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1961 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1962 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1963 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1964 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1965 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1966 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1967 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1968 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1969 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1970 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1971 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1972 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 26.0 FR France \n", + "1 43.0 FR France \n", + "2 41.0 FR France \n", + "3 37.0 FR France \n", + "4 45.0 FR France \n", + "5 45.0 FR France \n", + "6 71.0 FR France \n", + "7 59.0 FR France \n", + "8 42.0 FR France \n", + "9 47.0 FR France \n", + "10 42.0 FR France \n", + "11 35.0 FR France \n", + "12 26.0 FR France \n", + "13 36.0 FR France \n", + "14 33.0 FR France \n", + "15 54.0 FR France \n", + "16 62.0 FR France \n", + "17 86.0 FR France \n", + "18 167.0 FR France \n", + "19 251.0 FR France \n", + "20 308.0 FR France \n", + "21 268.0 FR France \n", + "22 199.0 FR France \n", + "23 145.0 FR France \n", + "24 85.0 FR France \n", + "25 55.0 FR France \n", + "26 61.0 FR France \n", + "27 79.0 FR France \n", + "28 105.0 FR France \n", + "29 120.0 FR France \n", + "... ... ... ... \n", + "1943 59.0 FR France \n", + "1944 64.0 FR France \n", + "1945 97.0 FR France \n", + "1946 93.0 FR France \n", + "1947 80.0 FR France \n", + "1948 116.0 FR France \n", + "1949 149.0 FR France \n", + "1950 281.0 FR France \n", + "1951 395.0 FR France \n", + "1952 485.0 FR France \n", + "1953 544.0 FR France \n", + "1954 689.0 FR France \n", + "1955 722.0 FR France \n", + "1956 762.0 FR France \n", + "1957 926.0 FR France \n", + "1958 1113.0 FR France \n", + "1959 1236.0 FR France \n", + "1960 832.0 FR France \n", + "1961 459.0 FR France \n", + "1962 207.0 FR France \n", + "1963 190.0 FR France \n", + "1964 198.0 FR France \n", + "1965 224.0 FR France \n", + "1966 266.0 FR France \n", + "1967 219.0 FR France \n", + "1968 176.0 FR France \n", + "1969 163.0 FR France \n", + "1970 195.0 FR France \n", + "1971 308.0 FR France \n", + "1972 213.0 FR France \n", + "\n", + "[1972 rows x 10 columns]" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -122,7 +2123,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 6, "metadata": {}, "outputs": [], "source": [ @@ -152,10 +2153,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 7, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -179,9 +2178,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n" + ] + } + ], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", @@ -199,9 +2206,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXm4HFWZ/79vb3fLvcnNSlYSMLLLFgO4IghBcYAR0Yyj5IfM4DjMjDo6CjrKIIMDjssMKDIoSMAFEFBQZIlhiSIEErYkZCV7CNlulpu79nJ+f9Q51aeq61RVV3Xf7r55P89zn9tdVWep01XnPe9yziEhBBiGYRgmLolaV4BhGIYZHrBAYRiGYSoCCxSGYRimIrBAYRiGYSoCCxSGYRimIrBAYRiGYSoCCxSGYRimIrBAYRiGYSoCCxSGYRimIqRqXYGhZOzYsWL69Om1rgbDMExDsXTp0t1CiHFB1x1SAmX69OlYsmRJravBMAzTUBDRpjDXscmLYRiGqQgsUBiGYZiKwAKFYRiGqQgsUBiGYZiKwAKFYRiGqQiBAoWI7iCinUS0XDs2mogWENFa+b9TO3c1Ea0jotVENEc7fioRLZPnbiIiksebiOheeXwxEU3X0syTZawlonna8Rny2rUybSZ+UzAMwzBxCKOh3AngPNexqwAsFELMBLBQfgcRHQtgLoDjZJpbiCgp0/wYwBUAZso/leflAPYKId4G4AcAbpR5jQZwDYDTAMwGcI0muG4E8ANZ/l6ZB8MwDFNDAgWKEGIRgC7X4QsBzJef5wO4SDt+jxBiQAixAcA6ALOJaCKADiHEc8Lac/guVxqV1/0AzpbayxwAC4QQXUKIvQAWADhPnjtLXusuf9iyfX8fFq7cUetqMAzDGInqQ5kghNgOAPL/eHl8MoAt2nVb5bHJ8rP7uCONECIHYD+AMT55jQGwT17rzmvY8lc3P4vL5/OkTIZh6pdKO+XJ45jwOR4ljV9epRUiuoKIlhDRkl27dpkuq3t2HxyodRUYhmF8iSpQdkgzFuT/nfL4VgBTteumAHhTHp/icdyRhohSAEbCMrGZ8toNYJS81p1XCUKI24QQs4QQs8aNC1yKhmEYholIVIHyMAAVdTUPwEPa8bkycmsGLOf7C9Is1k1Ep0sfyKWuNCqvjwF4UvpZHgdwLhF1Smf8uQAel+eekte6y2cYhmFqRODikET0KwBnAhhLRFthRV7dAOA+IrocwGYAlwCAEGIFEd0H4HUAOQBXCiHyMqvPwYoYawHwqPwDgNsB3E1E62BpJnNlXl1EdB2AF+V13xJCqOCArwK4h4j+E8DLMo9DAiEEZMQ1wzBMXUHWgP/QYNasWaJRVxueftUjAID13/4wEgkWKAzDDB1EtFQIMSvoOp4p32AUDqEBAMMwjQULlAaDxQnDMPUKC5QGgzUUhmHqFRYoDQbLE4Zh6hUWKA0GCxSGYeoVFigNhmAvCsMwdQoLlAajwPKEYZg6hQVKg8FOeYZh6hUWKA0GyxOGYeoVFigNxqG0sgHDMI0FC5QGg+UJwzD1CguUBoN9KAzD1CssUBoMFicMw9QrLFAaDNZQGIapV1igNBj1Kk9++/I2TL/qEezrHax1VRiGqREsUBqMehUodzy7AQCwcU9vjWvCMEytYIHSYNSryYu3/GIYhgVKg1Gf4qR+68UwzNDBAqXBKPBiXgzD1CksUBqMOrV4scmLYRgWKExl4aVhGObQhQVKg1G3+6EQ6ygMc6jDAoVhGIapCCxQGox6tyjVefUYhqkiLFAajHrtsNngxTAMCxSmotS7BsUwTPVggdJg1GsUFfvkGYZhgcIwDMNUBBYoDUZ96icMwzAsUBoOk8VrycYuvP3rj2LPwYGhrRDDMIwklkAhoi8S0QoiWk5EvyKiZiIaTUQLiGit/N+pXX81Ea0jotVENEc7fioRLZPnbiKyLPJE1ERE98rji4loupZmnixjLRHNi3Mfw4Fbn1mPwXwBSzftrXFNWIdimEOVyAKFiCYD+BcAs4QQxwNIApgL4CoAC4UQMwEslN9BRMfK88cBOA/ALUSUlNn9GMAVAGbKv/Pk8csB7BVCvA3ADwDcKPMaDeAaAKcBmA3gGl1wDW/8O+xadefsk2cYJq7JKwWghYhSAFoBvAngQgDz5fn5AC6Sny8EcI8QYkAIsQHAOgCziWgigA4hxHPCCmG6y5VG5XU/gLOl9jIHwAIhRJcQYi+ABSgKoUMSjrJiGKbWRBYoQohtAL4LYDOA7QD2CyGeADBBCLFdXrMdwHiZZDKALVoWW+WxyfKz+7gjjRAiB2A/gDE+eQ176jRq2Kbe68cwTPWIY/LqhKVBzAAwCUAbEX3KL4nHMeFzPGoadz2vIKIlRLRk165dPtVrDIL661p16MQqEsMc8sQxeX0QwAYhxC4hRBbAgwDeBWCHNGNB/t8pr98KYKqWfgosE9lW+dl93JFGmtVGAujyyasEIcRtQohZQohZ48aNi3ir9U+9dOesoDDMoUscgbIZwOlE1Cr9GmcDWAngYQAq6moegIfk54cBzJWRWzNgOd9fkGaxbiI6XeZzqSuNyutjAJ6UfpbHAZxLRJ1SUzpXHhv21KtJqV4EGsMwtSMVNaEQYjER3Q/gJQA5AC8DuA3ACAD3EdHlsITOJfL6FUR0H4DX5fVXCiHyMrvPAbgTQAuAR+UfANwO4G4iWgdLM5kr8+oiousAvCiv+5YQoivqvQwv6lTiMAwz7IksUABACHENrPBdnQFY2orX9dcDuN7j+BIAx3sc74cUSB7n7gBwR5lVbnhMG2zViwujXjUohmGqD8+UbzCCOmzu0BmGqRUsUIYJVCdejHrRlBiGGXpYoDQY9a6B1Hv9GIapHixQmIrAmgnDMCxQGgyTU75eiLoBWKEg8Nm7l2Dx+j0VrhHDMEMFC5QGI9ApPzTVKCGuD2dfXxaPr9iBz/58aYVqxDDMUMMCZZigTE6N7sNgyxnDNC4sUIYJ9eLDaHB5xjBMDFigMJWhTgQawzC1gwVKg1HvJq2o9YvqzGcYpn5ggdJg1HuUV1x4GXyGaVxYoAwzhrvAYRimfmGB0mCYLEMqbJctRwzD1AoWKMOFGluK4hbPcpBhGh8WKA1GvXe8UU1uSrNiDwrDNC4sUBqMeo2GqpQvnX3yDNO4sEAZZtSnuGEY5lCABUqDYRIYdTOwjzoPhUUhwzQ8LFCGGbUyicXe4Muudt2IRoZhyoQFSoNhDBuuE+dDVHHG+gnDND4sUBqO+ux64652XKexBgzDlAELFKYixBUI7ENhmMaHBUqDUe8j+djzUOrDcscwTARYoAwTat0PxzZ5Va4qDMPUCBYoDUa9d7yRnfL1rnoxDBMIC5QGI3BP+Rr3y1EFQ5ilV3Z292MwV4iUP8Mw1YcFShks3dSFRWt21boantTa92CbvCKmD5JD+YLA7OsX4l/veyViCQzDVJtUrSvQSFz84+cAABtvOL9mdQjSAGoeLRVzprxJMOYL1vnHlr8VrQCGYaoOayjDhJo75dV+LDGjvIznAwQOwzC1hwVKg1HvrutqRXkVfSwsURimXmGB0mDU2ukeRPSZ8gGmPHWa5QnD1C2xBAoRjSKi+4loFRGtJKIziGg0ES0gorXyf6d2/dVEtI6IVhPRHO34qUS0TJ67ieTCVETURET3yuOLiWi6lmaeLGMtEc2Lcx9MfGI75VU+BolRkBKF5QnD1C9xNZT/BfCYEOJoACcCWAngKgALhRAzASyU30FExwKYC+A4AOcBuIWIkjKfHwO4AsBM+XeePH45gL1CiLcB+AGAG2VeowFcA+A0ALMBXKMLruFMkI+i3jUYE0H1LjTqjTHMIURkgUJEHQDeB+B2ABBCDAoh9gG4EMB8edl8ABfJzxcCuEcIMSCE2ABgHYDZRDQRQIcQ4jlh2T3ucqVRed0P4GypvcwBsEAI0SWE2AtgAYpC6JBErTZc6343+gRF/3QFXpqFYeqeOBrKEQB2AfgZEb1MRD8lojYAE4QQ2wFA/h8vr58MYIuWfqs8Nll+dh93pBFC5ADsBzDGJ6/hj2n5+qGthZG481BMAkPYJq96uVOGYdzEESgpAKcA+LEQ4mQAPZDmLQNePYHwOR41jbNQoiuIaAkRLdm1qz4nJZZDvRt+oiootgYSdJ7lCcPULXEEylYAW4UQi+X3+2EJmB3SjAX5f6d2/VQt/RQAb8rjUzyOO9IQUQrASABdPnmVIIS4TQgxSwgxa9y4cRFukymPiPNQAk1e7JRnmHonskARQrwFYAsRHSUPnQ3gdQAPA1BRV/MAPCQ/PwxgrozcmgHL+f6CNIt1E9Hp0j9yqSuNyutjAJ6UfpbHAZxLRJ3SGX+uPDbsCZ4AWFuqtcGWLVBYRWGYuiXu0iv/DOAXRJQBsB7AZbCE1H1EdDmAzQAuAQAhxAoiug+W0MkBuFIIkZf5fA7AnQBaADwq/wDL4X83Ea2DpZnMlXl1EdF1AF6U131LCNEV816GBTXbU14FBURMH3bRSxYnDFO/xBIoQohXAMzyOHW24frrAVzvcXwJgOM9jvdDCiSPc3cAuKOc+g4HjKahmPuRVIroM+X9NZACSxSGqXt4pnyDYeqw466lVSmqtpZXrW15dcJ9S7Zg+lWPoLs/W+uqMEwJLFCGGbXqeJXiUK3yeWKjxU8WrQcAvLmvv8Y1YZhSWKA0GIGLKA5JLSpfPvtQgMXr92DPwQHfazgmgalnWKAMMxp1IB86bHgY96ifuO15XHLrc6GurbVpk2G8YIESgVruf24qW/WztTYNxd0C2MRwn9hYkDe4fneP73W1XingzX19eG3rvprWgalfWKBEQO0eWI/Ub838KWogAeerWIfl2/bjmoeW12TAkC0Uyrq+VuOGd93wJC744bO1KZype1igRCBfSw3FcNzuaGs2DyVe8cEbbFXf5PWp2xdj/nObsLd36COocvlGHQowTBEWKBEoczBZWep9pnyVwoaHQilsSlmvw0AuH3Bl5WGBMjQIITD9qkfw34+vqnVVhiUsUCJQSw0liFpXLXr5/iavSt1X32AeO7u9Q26b09b2PP3ZoR8xhDV5xdUED3VUu/3oqTdqW5FhCguUCNTShxK8wVaNTF52+dHSB69RVpn7mvuT5zH7+oWe51hDGf7UOmhluMMCJQKFOnTKF6O8aluP6myvVTQzxvWgvLrFHKGUkQJlMFcDDSVfplO+5sbNxoRbrbqwQIlATZ3ydepDKe4YWaWlV4bgzpIJ63XI1UAqh9V6h/M8nKGANZTqwgIlArXUUIJnlNd6La+I6QJ2ZByK20onrLLjmJ9e2rwXuwNmu3tRboncL0aD2626sECJQD065Ws94c2mSmHDQ0HCFijRTV4fveUvuODmP5edLuxAIO6v3DuYw0OvbIuZS+NSh6/usCLufiiHJLV1ygecb9AXph5WG05JgZKN+fu+ub/8hRuH6me79uHXce+SLZjS2YJTDx89RKXWD+x7qi6soUTAFOG5q3sA33tidU1MYnWz9ErkeSghZ8pXURFLJa3XIVsDp3zYny1u2PCb+/sAAN39uWgZ1JgdB/pxwjWPY9VbByKlr3XQynCHBUoETCavqx98DTc/uQ7Pb9hTtbKDTCM1c8qr8mOavEzyYijuKykLz9Vg5mq5vq+ogjvuzpq1ZuHKnegeyOHOZzdGSl/rAddwhwVKBEwmLzUhrppzCurd5FW95etFqOvioKK8sjWYE1JuiVFH2sUleqKlrzXNaTVXKJrQr/X7MdxhgRIB0yjHNkcMYV3c1NpGHHcLYPP56lNLs2H5Jq+oGopMX4UW3dsziO3SpFYtmlJqNYNok09rHQU53GGnfARMGooyJ1SzQzJuAVzjJTnidlT14JSv9q6TfpTbbnE1lGrc48nXLQAAbLzh/MpnLlEayoGIWyCzPKkurKFEwCRQEvbbOnR1KRJvYmGliO1DMXrdh+6+aqHlhdZQYv7OjT4xsr05DQDY3NUbKT37UKoLC5QImB7KxBBoKKaOtfZLr8Rz9gZ1kEO5wVZUn3wcYV5u0ri/c6P2q3EHTA162w0DC5QIBGko1ezU68E0VI0K1EOwQVwfSpzfPaxWFLeONVWiK4Bq46jPA2so1YUFSgTMD+VQaCh+Jdf+hYlc+hAvKeOVX9GcVLk8w6ct7/rIAiWmU7/WxDVHNuhtNwwsUCJgWpljKBzjpqxFwPl6R3UUQ+VB8fuNaqGhhCW+U72x56GImBoKC5TqwgIlAiaTl3pYa+H3tF+UWmsoUU0RSkhXeYMtOz+PY3H9UHFGz0OtoTQqcZ+DWmvwwx0WKBEwP5T+o+xKUA9b5foRefn6KuVbTn5x/RNxqhhaGNmBH9HLAuJ3zLUymcUVCCxOqgsLlAgEaSjVjfHyF2a1mtgYd1LnUC8p4/UTxg3JjdPZDZWGYoe2x2zRWg1c4hZbj5vjDSdYoETAtJZX0epUu3kMtX5fqraWV8VNXuYMI5u8YtSxbAEReWJjvMADu/gaayjV2siNiQcLlAiYRjnqYa9F2HBcZ2Vc4oajDvWOjX7l1URDCXldpaL5Yo/0a9UxxxWEbPSqKrEFChEliehlIvq9/D6aiBYQ0Vr5v1O79moiWkdEq4lojnb8VCJaJs/dRHI6LxE1EdG98vhiIpqupZkny1hLRPPi3kc5BJm8arIWFOKN3CpWj8jlB5i8huK2YjvlozNUExsr5ZSvVcdsayiR01euLkwpldBQPg9gpfb9KgALhRAzASyU30FExwKYC+A4AOcBuIWIkjLNjwFcAWCm/DtPHr8cwF4hxNsA/ADAjTKv0QCuAXAagNkArtEFV7UxCQxhn69e2UYPyhD4b6pJsIZSWfyEfmSnfKxV78OVWakFLOObvOKlr1W5tR5wDXdiCRQimgLgfAA/1Q5fCGC+/DwfwEXa8XuEEANCiA0A1gGYTUQTAXQIIZ4T1q99lyuNyut+AGdL7WUOgAVCiC4hxF4AC1AUQlUnaIfYajr+ApcoqdEQLO4cnKBkxQ60MkNsr3rG39Nl6JzytV5tuFb9ctGHEjV9BSvDlBBXQ/kfAF8BoHexE4QQ2wFA/h8vj08GsEW7bqs8Nll+dh93pBFC5ADsBzDGJ68hweiUt30otTB5Wfz0zxvQHXEl1srUI56z1Lh4YcWd8qXEXS063tIr4Sj6UKKVUymnfK3mc8QvlSVKNYksUIjoIwB2CiGWhk3icUz4HI+axlko0RVEtISIluzatStURYMI0gJqMQrS3++7nts05OXHXrYkcD+USjvlS/OL21k3wtIriBneraiZT972oQy90K8Uz6zZhadW7ax1NapCHA3l3QAuIKKNAO4BcBYR/RzADmnGgvyvWm4rgKla+ikA3pTHp3gcd6QhohSAkQC6fPIqQQhxmxBilhBi1rhx46LdqQuTU75QUw2lxs74Cs1rCA4bjtNpF9P6dSw10VBClkkxJzYWzXpxTV410lAaVLPSmXfHC7jszhdrXY2qEFmgCCGuFkJMEUJMh+Vsf1II8SkADwNQUVfzADwkPz8MYK6M3JoBy/n+gjSLdRPR6dI/cqkrjcrrY7IMAeBxAOcSUad0xp8rjw0JZpOX+l9NH4rpRNWKLIvoYcPVj/Jy5FGFsOFYPpSQ11VKIMSlViP9uKsN14E8GdZUYx7KDQDOIaK1AM6R3yGEWAHgPgCvA3gMwJVCCLWP5+dgOfbXAXgDwKPy+O0AxhDROgD/ChkxJoToAnAdgBfl37fksSHBZPIaismFpk5LPxolNHQwV8Cvl2yJ7uytkG3eVPdKjCz1PPw6/1p0VkO3lpf/wyGEwLZ9IbbxrYJADkN8TZglSjWpyBbAQoinATwtP+8BcLbhuusBXO9xfAmA4z2O9wO4xJDXHQDuiFrnKBBZL755pnwNTV4xy7z1mTfw/QVrkEklcOFJ0eMbqrYFsP0pepSXXoS/ySti/rGc8uUljulCMaa/8y8bce3vXscf/uW9OHZShzEfv2dciOotQmlrKBHTszypLjxTvgzsHRlrqaGEyJsidLp7ewcBALsPDpadFqhE2HD1TV4ODaUK81CGdi2vaOUEhQ0/v34PAGBzV49vPn7FmwZclSC+76dCFWE8YYFSBpmk1VwDOe+JKOpZrclaXo7P5Zev7i0bNMmmSgQ3WQUCRrUsvHKL+/vFqWHYIiu2Y6MheXGrZf9BSTUmhoYh7hJDtQ5eGe6wQCmD1ow1sb93MO953p6HUtWJjeUdD0sqaXUgWYOwDEvcRftM2lXlNRTz+ag/X5zfvXyTVzwfiqmqxd8hqHzzuUIVxySViiZkqgMLlDJQL2PPYM7z/NA45YOPRzF5pZWGErHy1Z4pX4km1W/Nq0OOOzF1KJzy6reN+4wF7ukTwnlfft7xKQqroTdLMsGwQCkL62HsN2goikZ0yqcrZPKK7iz1T6naNI6z1xnl5XG+oK6Llv9QhA0r4pu8vNMHzQdS1MyHEjc9y5OqwgKlDNTDGLQfSnVtyP5lA9E63aCAg7DE1VCMK6/EtJ0DcCze6Gfyir58faRk5ZUZc0Vkld40bhAhBbdvlFcVTV5x1/Kq9fyd4Q4LlDJQj6J5+fp4NvhYxCwzYUf/RMOehxI5bDggyitSrk6C5qGo3y26yav6GkrciY1Fk5n/wCQRaPIyn6umhhJ/PxSmmrBAKQP1Eufy/uaCqmooVcq36AOpzSsXOA+lwhMbvYS+iDn6jTWQKDNtVE0y6He2s42hoVTVh6J+o6jp2StfVViglEGghqL+V/OZNUXnVEjUxK17tZzylSDIKR9/x82hi/KK75Q31EOZvALSu5vPsU5aNaMca5ye8YcFShkE+VDUBbUYBcVe8sQ2WUUsH/FGjsEaSsSMHXloJi9PH4r6XwsfSrllxXPKB6Uv1+Slf6/m4x/fz8UipZqwQCkD2+QVoKHEeaEGcnl09ZhnqxvX8tIORwmEihv261mRcpIFhKuWE+V1/9KtuHnhWo8y/InbWQ1J2HDMJU3IdsqbTF7h2tn9HDrNiVXUUGqkQTPhYIFSBrbJy+BDKc5Dif7UfvbupTjlugVlp4tr8lIdedy1uOpBQ/nyr1/F9xasKTke1OnFnUcUa+mVISoraBHPsNF07jbSv5uEVSUo7ocSNX3l6sKUwgKlHAJMXupFiuNAfnq1/yZgYWbKRxnFxt3+1qseZaWLeT4MTh+K13nh+F8ucdoubJmVmtgYpOkG1cf9jOv5+SX9u/kv4uoHXwtXSa9yI6e0YJNXdWGBUgZBTnl1vBIDNOMClPGz9iSuKaWooUSVKOV1YFHQ29Qrt9hLr8TRUIbIh1JM759vuRqKfr1f2PAfV+7Er17YYjwfRKEQrn4mwiTL5gs4/dsL8Ydl26MVcgjDAqUMgnwo+YAR7oH+LF7bui9UWdkyF0SK293Gnd9gO+XrWENxOo5Lcyxu3lSLUWxIDaVCvi6jpiv/B5utauRDiZk+TN329gzirQP9+OZDK8rO/5an1+Ha35WfbrjAAqUMii+bd2evRk+mEdqVv3gJF/zwWQzk/JduAYBsgJ8m7PGwFH0o8YiaXrWdeQvgCmgoAVFecdfyGgoNxV5tOKYabKxrSJOXn4ZSzSjHuEK/vHTll/Gdx1bjZ89uLDvdcIEFShnYYcMBGorpmV25/QAAYE+IPUfKX/U3rlPe+h/Ukfxp7S7s6y2tf1hnrgmVLGjplTg4781PQ4maf7R03rWpXln++YYTqu7TQZNGy+GlzXuNmvxQ7IdS1Zn+wxwWKGWgzDqBPhTD+daMtUHmwQHv1Yp1Bg2LLYULGy7fIRLGKd8zkMOnb38Bl89fYrym+js2RifIKR936ZxYS68MkQ8lSMNRR4PawF1+JaO8PnrLX3DBD5/1PBe3rw9TNbUSBsuV8mGBEgGTD0W9pKaHNmgOgH6Nca5LgO07MiFMXqpOa97qrnj5gRb7mNqDldZ/FB3GnPKdx1bh7uc3eZ6Lp6GESyxC1DFceab8o2koeoZhhN3+3mzgNZ7lxnzSnJNbvfMyDeaYYFiglIEdUmnoOXKFcC+jr0CBfxnmulV/OOXX2diHYo7ujSYvhCtg98EB4zmHhuKRT5jO9Jan38A3frvccLb6GkrcuTKKoMUhg7SMUg1FGM95sWL7/sBrvMu1/ke9ff22TLdoaygRyziUYYFSBuoBM2ootg8lhkCRPWrQemF+RAoBDhEuarpvRzYRig5D2I70nO8/YzwXvGNjuDLM+UdLB4Rvt+LiiPFMi6bfOaofSb88TDtE9rXFfsCCBZ/fnkADuTxuW/QGBmPubDpcYYFSDvL5izoPRfXzYZx+5ToGK2dyMudkT9z0yaFay5aENcXs9TGlhN8CuDr34J82pMlL/o+7CZixvJBtEFdDCTM48S23AkI/SKAkE6Ujs58sWo9v/2EVfr002lya4b4fCwuUMgjtlK+myctkqohrAgkxYSyMSS96lFeAzV7+j7Vve0BnEj/Kawg6C1uL8i6ruz+Lh19905w8oD/2M3n5La7pECghfqNaLSMfNKgAiiH7SQ9Vf+vePiufiPWv5rI09UCq1hVoJILChoNWqw0yZ1nXOPMqqYOpbsYcw1Ec+fpoKD7RLzEHjnZ6U4RaJZzyYZcTqeelV4Jmsl/z0Ao8+PI2zBjThhOmjCw5H/Q7+7Wz0wdlyBjhtKfIGoqvlhxM0KACAHI+Gkp3f06eizYWd2+hYFoMtVFhDaUM1LNgehlycsJjoEMzxMtU7kimUvH5/hpKwXiNOhSmGl/7zbKS9ZyCkqn7izNHQAjvz4pabi8b2ikv/5s2edsjV6o2BScE3aOf2U9/JuOGDUcdqccdOIWZLzPoI1CUOcw0ubmc8oejtsICpQzsTs0YNgx53ju9ejz9RmeBW7RW6RkMsxNe0YfiY/IK8cr/cvHmkvWcArWHkNf5UfUtgCOlkmlDJlZ1yxk6tLamJACfuU72PfrXw0tw68+9n8krjGCN2pnG3WJAx5SHEtZeAqXY/vFNXsNQnrBAKQc/+zIQvJaXcF3nR/kaSlmXR8pHvUTeJq+4o3vrv8kCEDcCy53Wex6KMJ4Ll39wwi1dvbjz2Q0lx8MWqYowRSKlk9YrbRI4QVFiRcFP7HB9AAAgAElEQVRdei7vIzSE4TpH2VqmpvoFEde0GkZDUW3rIU/s9zKuQHR/Hi6wD6UM1O9vXBwy5MNm2k8FgK3GmF7KMCOzKHbZMOGofkEHlXo1AtfyilFQ0KS2otkvoobiMKl528cvu/NFrNt5EBecNBmj2zK+9fEuw7puMOd9vXIkm/rrINNkcWDgoaFoz637bCHEyLsSnWmlBi56Xm6yBT8Nxfof9I6bfn/9dxmOAoU1lAh4RghpD5hJGIQJG/Za9dfx2ZAu/gzigAKgaSge54Zq0b44PhRfpzLihw2HiSDqlaao3kGnSSq0hiL/m0b4qhMLdrr7n/fqMPW2d/sBHc7uAA0eMPuAgihqUFF/I+/POmodPS/Huyo36DkMJ1R9s2hIWKCERDhehtKX2e9lcxPOKa+XHaZ++ufyn9QQ8sR2RPqZvKKvg+VfvjpfMR9KFUxeejJTPZvTlo+jd9C14nTIMlW+JpNXgpzXmdKbowjN5x0+FEM637IrMDqPGzgRRktSwjrloaHYVogAgRhkcvQrv5FhgRIS/bf3HL3pGkrAzEZfp7zHel9hRr5OgWLMHk+u2oFNe3pKjodxdoYZVVbLoR13vxUg2HGs+oC45hgrD+9rUknrB+7POgVKuWt5GU1eCf8dHe2BQ8Bz5Gny8onyCvP8OzSUiFI77PJGRhwainceg/I5T3gIFL+gBZ0wbVCruTjVJLJAIaKpRPQUEa0kohVE9Hl5fDQRLSCitfJ/p5bmaiJaR0SriWiOdvxUIlomz91EUm8noiYiulceX0xE07U082QZa4loXtT7CIv+0wdFwJQ7yzjomjCPnWMLVp8Un7lzCc787tOl6QM0BMDfblzUIPxqaabYkfmfj4ND6Hqejzn6DTECV1F87g41bJlBTnnVCRp/qwBNzzbpBJi83A0YxpSkj+qjdqZqUFNNs6Q9D8XDoefXPo48TEsnhWinRiaOhpID8CUhxDEATgdwJREdC+AqAAuFEDMBLJTfIc/NBXAcgPMA3EJESZnXjwFcAWCm/DtPHr8cwF4hxNsA/ADAjTKv0QCuAXAagNkArtEFVzXQR2xBL1vgw+Yz0lcdjilE05QyrIZiOh+mM/UbVQYu6RFA0dRi6ugi5OnKK2g291CEDdurSbuegbAlhjV5mX6H8OHZpeccAqFEoPhrf0CFNBR5337JZ1//R3z51696ntOTmdri969ZW/96OeXjBt6EGXj+5Y3deGrVTt/865XIAkUIsV0I8ZL83A1gJYDJAC4EMF9eNh/ARfLzhQDuEUIMCCE2AFgHYDYRTQTQIYR4TlhP4l2uNCqv+wGcLbWXOQAWCCG6hBB7ASxAUQhVBfUcEHm/DE6nvHceoZzyHlFeYV5WYfgcljAail+oZ3EOTtSRo7MebqIEHZR2euZzehnVtM8r3H64css0CRQ7ysuQX2DggY9QdYQNu34Px4S9AN8EEP05yRaC34Wd3QO4f+lWz3NhnOJLN+0F4C1QciEFStACsoBZS/vkTxbjsjtf9M2/XqmID0Waok4GsBjABCHEdsASOgDGy8smA9Bns22VxybLz+7jjjRCiByA/QDG+OTlVbcriGgJES3ZtWtXtBtE8UFIJxMQovRhcIw8KuCUH8rNmgDdtm5O7KdZ+Tlzw1DcSybYVOBXxxO15Ub8FjD06hCCNJSg36Qcc0bWdUFYgVnUULyvV1Fe5tB2WV6QwAnwE7pPh5mwp49HwggUr/bWBXE0rVWvj38G3gIl3GoYpsEXm7wCIKIRAB4A8AUhxAG/Sz2OCZ/jUdM4DwpxmxBilhBi1rhx43yqF460slFHcEoGvezO/IqfgwICrIuKH6N0iOXMlPfOO7gMP8qJ3vGrR0dLuiRPhSNSz+OFDzL7hZl7YCpbYT8DLg0jbOcS5ENRnaBXJKJex6DJt4FRXm5hrfuPDDejt3nQJnWmOuTywW3sh98Cl6XXmssPeocHst7tX46vNQ4bdvfgO4+tqsiKAuUQS6AQURqWMPmFEOJBeXiHNGNB/lfGwK0ApmrJpwB4Ux6f4nHckYaIUgBGAujyyatq2BpKymoyd+fi8KEEdOh+D5JtFjM8eMa8Q4xwfX0gtooSMb38H/UlUfcVbpZ1sGAESjuEoBFykIYSFNmjnw1qhhINI2S7hfWhBO7ZE5B/0FpeJW0bQpiGETpOP0vpPWZjdsjOuUje6SePajHmH3ZF8T5XFF+x/OCBZyX4u/kv4pan38CWrr6qleFFnCgvAnA7gJVCiO9rpx4GME9+ngfgIe34XBm5NQOW8/0FaRbrJqLTZZ6XutKovD4G4EnpZ3kcwLlE1Cmd8efKY1VD/faphEGglGHy8nXKe0xM06827gUeQpX31zDUyMvsJ/F9AezO2HyJH4EmL/1aP03LJ9JK/+bpBwsQ+EGrhZSznpW7ncM2m7rOtMFTQj4/JoGjHj2j094ndNo5D8VnQBVDQwnS9HXNK5JpN4QPRYV0+2lpQRpKGIFSTeVhoEYbgMVZeuXdAD4NYBkRvSKPfQ3ADQDuI6LLAWwGcAkACCFWENF9AF6HFSF2pRBCtfrnANwJoAXAo/IPsATW3US0DpZmMlfm1UVE1wFQnqtvCSG6YtxLIOpBTCe9zVaO0W+AuaNcDcUZYWaon+GzThgNw2877TCaQWQNJaRt37o2nIbiF4nkN5coyKFtIox93Ov3dacNU4ZRk1MDA8OgJcispzr9IKe8W7iGMSWF0WL0dvF63vT5O3E1FFP6brmagacPRz0jAQJFLXMfpfxKoLIe6tXxIwsUIcSfYV566WxDmusBXO9xfAmA4z2O90MKJI9zdwC4I2x946IeBLX4XolTPkT0hn1tiGG8qWMMs8aX0Qfg51RXHZWvhmI+p3KOOz/AuPWx3gZhTHcedXH6UDw0lIIouU4neLmN4A5T4TZ5lbuWl+mnUII5a7ggyGSTtzUUr3OaduA6t7+vuFOmOcrLv/3dab2e1z6HQPHMwhfHwMtUT3uJevOgI0hDWb/rIN7/9lKfrUMDq6pA8R84VAueKR8S9QOlDBpKmAclzPL1CscLr3eSYaJjDDqKqZPRy3tqtTkSzt8H49/RBRGk4QTNA7LP6YLVVReHYPaoaHG1aEMdy+jBTC9ycR5KPJOXeZM363jWMJM+KPhBtYt3lFdpPop/+uXLxnPF9MG/oS5EvDUUPcrL3OGbCDJ5FQpC86WVns/5tA9QXK6lZGkdu0x94Be+ruViPyeN5JQ/lLA1lBA+FD+zERDQMdlLrxQP6QLC7JQvrasb35G99vm1rfsipy9HQ/ESEqa2CaOlucv3Cxv2MgkFj97DayimzkAJlJKw4ZDNFrQ4oW3yMq0l5aOBAP5t4OeU1zvQMM9fmKVJvK5xmLw8zpt8R8W6+WuRfkv063UyDa5UlF2fQaCEMYsqok7+1MsxRftVCxYoIbF9KKngsOFAp7zPea895cM43PVrjhzXVna5+stlcujpnXDJfhg+ozpzmaXlm9I7Fh/0eUf8bNRBHZq9gZipQ6yAD0URWUMR6n+AhmIwbxajvPyFpp9D2q98IFxQiHmTOv/Bk/5sBq2p54XzXfKvo7eGInzLUfmbNJRywoajrsgMFH/fOEIpCixQQqJ+FzvKy/Vjq4cjk0wErjTq65SnUoEVxiEtAEzpbPG5A38fit6/ZJLej4Wfw1R9K0dN95o5HWbpFb9IND+Tht4Zla6lpZs6TJ2dsdiS8gLX8orsQ7H+G0f4so6mkWmQycuvw3TOlDcTx+Sl/y5ez+tggEAJ6oSDfqNcQB1VnYJMjn1Zk1M+uA0UgzG0i0LAc1ItWKCERGgCAzBrKOkk+UQJWf/L9aHoVxtHyULYIaMmfDti7bPXDGErvSZQDB1iOSYvr9FgmJny/hqKeRStm0P8TJamWzCF6nqVbZ7YKOsS0dkUtJaUGuGbnd7W/yANx9vc4++/cNfR73iYsGGv53UgV/CdaxO0E2RQ2G6QBpELbF/reDgfim9VY5mrVN5BJsBKwwIlJLaGklSLNzp/KPUgpVOJwBcqjHO3bJMXgkMEfcNtQ73sxXt2d4j2yDng1hwmDQ8zodHkpY/sfDszrSzXZYM+GopJI9T59B2LjeValdQ+BnYWZoHmR2DggG3yCpopb6iXr8lLL8d57uOzinOTzQ5//bcP3i/Eq00Gc3m0ZlLG80GDtWzeXygGDSz8wqqFEHaah17xnmcdRost1jWOduE/8KgWLFBCsKWrF0/K1T9T9p7dzh9KdYiZZLBA8fuRvfZDCdOZClGc1GZ6Tv0eUP2MaWSkp8+6RuthVxs2OV2LHaX/yNlKF9bk5czLaS5xLX0SYun5TXt6tetLrwkz+rSdpa706rvXpk7Oevo/Q6odTaYf209kyLs4MPDvbPX2+8Oy7bhvSXE5PtMz6vztPS8J1GKyeYGWTNJ4PkigBGsg/gLHT0MJ03ebBlRexNEuihrK0AoU3lM+BOf9zyL0SBU2bdhvQj1gmRAaSqi1vPTRriMPc5qgOUxhJwSGM0c4ryk65cO/0F726nDzUMz5+5mdwmooYax2eSGQcLW4nizIj+AW2uq4ydzormfQEjVmP565fmEjnABnWPJ9S7Y4rjPde84gkExleGoo+QJa5K6XXnkEmYlyDoHiX35JlKAmcIMWz1R1Sbn8kU5N2Leq8QSKoU7VhjWUEPRo9lBl8nI/y6E0lBB+Bq8HNpzJS9gSxRTB4+/M1q4zjGr0l9HtTyjW21iEdd6goRSjl0zpip/9BWNpnRS6k9Pt8HXu9RFC4AdoKOZOtSD/l3Y+AAL9YKqeQeHVppGpLSiCOlOvCCjtnga0tmzNJB3XGevm0Eg9L/HdMyVfEMgXhF2eV3+bczxTQZ2+hwbis+eLc2JmaeHqNx87ogmAs98o1kn/XFp+0OTbsASFj1cLFihlkjI55ZUPJZnwiRIK1lC8nNuhFocUwSavMGt5AX6zrM0j/LBLrwRpKOa1vMJ1+H42+MF8AZmUt8mynJEj4N3Gjs7CkE6V6x59hu08ggYl6vxb+/t9yw+McApwyusDiua0S6AYbqV8DcV5jSqzaPLy0lC059hDavkJDHf5fmHnXtVX17c1WfVzb/PsziPIBxQUBOKHbVodYpMXC5QySRlMXpu7LPt6UzphHB3anYFvx+68FnCZUnxMQkEmrzBreQHhNBTz8ugBdTBsA2ubckKYvEKtmuxx3WCugKZkAskE+Ud5hZgV4tXhOvw3QWZPg1M+qGx9nojXCFeVu7PbIFDk7xYU4eQ9x6P4We/s3BpKqDkmYcyq+dLfTy/PKw9dUHsJHEekmk+HnkpQiUDJFvzzVtm1yaABr8mNfmvNAc57jhM2XFzslQVKXaPmobgftq//ZjkAoDmVNK/0GsIpr844TV7Bo/OCELb93dSp+26QpZ0KM3p0P+xhfShZg5ZTXEfLO13Y+P18QdhaSImGkrM0lGSCjBqWVRffWzDWQT9kepGL+2k4C1GDkH7DPhqe9fQoQp0/OJDzX9wwwBzk1Znpz6RToDhdscY1ssKYhX2EzkDe6qBb0ilHfjp6vT01lAAfimqXdDJR8hzo74/Xq6TqqzQUrxWHg1ZT0J+LL93nvY1xGIrBH2zyqmvSSW8NRdHalDR23EEROkDxgTNusOXT4QY5dH238NVNXqE0FOc1amQdNJtcz1sXXEHOZocfKcDk1WRH4pWaTDKpBFIJKhGaQbb30vr4d1ZmLc8q193Gen2Wb9tvLDdfEJ6RgMXzsPP3WvFA1curdnmDwCjWXRc4xc5S11CSCb95WFKgpBKB+7W4y9Pr5Keh6BtbeWnRQbP9VftY88mc550TcUvzVvm1NUkNxUugaMm8nnX9udmwu6fkfFhULuyUr3NUp20yabRlUsgWCt4OwYBOE9AEiiHqyM+UYmsohrzDaEaAWfA4F+4zaSjGIqx0eYOGEqDhBK3DpddDbYLmvi6bLyCdNGgo8rtfh6jjqaFox4JW+3V3dvraXm/u894USc3mt1e89jJ56Y5zD21H3XfQ6NhLoOQNgw49kCBBPmHDsuwmn0hI54DDW6AUo7z8NZQgDcbPh5JJJUsFiqNu5rTK5NUfweQVdcKrG9vkxT6U+qb4Mnufb2tKQgjDw57zjvDRsRfvM9j0/XwMRad8sIZRmj7M6Fofwbo6ZGHuqHRMs9V1k5dn9Iv2eX/foDH/ghD2agYlI9y8rqG4NQTru5ft3AvP0aUuUEzroclr3KNXXVj3G9Kq7O3VGgJ8CAP50g5NCQ2vW9TNbZ4CJe99XgnHvz55MhJkbr8wofU5x6DFPSCwvrf4aCh6vbxMzw4NxsupL49lPFa8yBm0a4XtQwlp8vIzO6YShDFtmZLzYVF1YZNXndPRLGfpGl4aNXryCgu1R+E+HbvXTFz9cr9JY4FzGCI6s4vpzS+jSlHO7F+TvdyreCGEbW7029a0IIS9gKdbCxjMFZBJJpBMlJpcdPu+qZnam4q+Aq/bDDOXR3VKe3uzzuMB5iZA62yU2dWjEnpab6EgSuqq0J3IXj4U9dNlUomSddGSCcL3P36iJVACnPJ+AiXrY1YqjfLyv3+vMnoHc/bSLV4mQXsJpVSi5D681p7TUW2qfEreAkX77BUUkFd5JGNFeXmZzocCFihlMrLVGjXoD8Nx33zM/twkBYp7dDRoGJnr5AvFBQqds4r9I1NU2iCTl99EqZzmzDZFcA34jf5sYWksQubtbTYJWvqkUAA6WzNIEHCgP1ty3s6nUNQiS5y6ug/FsLhnKklGLWswX0Bna9ozb/e9mSJ0VKe0t2fQ87heFzd2Zyfvz73fC+DSUDw6JNVhe5WgFjRMJsizM1PPXks6WdJxp5MEIvI1GYZxymddgkpH+W1afSY26r4dL19gXzaPTvkOe4X1qjK9gmv0CDBv/5V1bESTOcoraOCkymxrSsXaxldl7beqRDVggVImx0xsB+B8MPQJTJNGNgMofZj1mcWml2lvb7GTcdjjDR2vjhBAMmhSnI+Gki8INEuBYnLK9w7mbS2hZMdBVe8ADcVpwy41eZnyELAEZntzGgf6/ARKwaglbunqhZCanGmDtJRhHpEQwjVL26NDCWE2VOn039pdV5NAU8KgOe09Fwpwtm/p5FNhm7W8HoW+QevcqJa0r1O+JZ0siaZS0Y9EPmHDIZzyfpragFtDCViN2Mvc05ctYJQcFHj5mFTdmz3C/1Vn35RKeLa9raH4zEMJMnmp+2/NWG0cdZMtDhtuAC4+ZQo6muUI1eOH/vAJhxXX+nKNbh58OXitI90Zqyf3WyVXP26bvAJGiKZzzT6dJWCZC0a2WPcfdbVhh1PeMDvde9KYNc+moyWFA4b9ulWeTSnvKK+Ne3rx6tb9SCVLo7zs0b9hhJ3NW8tu+NnvHSZBDw1FCGF3Unt7s84FOR2/sfe9qdF7c8pnHkau6ENyj3D1JVK8OqreQatdR7akvcOGZZqWTBKDuWJnmS8UbDNcMmHW8NTvnUmZJ//qz7p7xV63ySvIh+IlcPoHNQ0lZ9YgmtOlGoodVJBOei6tr37+EX5RXgFzlXQNBYg2F+W3L2+znzOO8qpjjpnY7hvllc0X7fzuB+Ha371ufzb9yKZtVLOGjlenIAQSQWHDPhEf+UIBTXLka4o06RnIo0MKFJM5IOj5NfpQdNuyp4Zg7RXT0Zx27F9ekn+hYAtG04ZgQRqKV4eofs9RsjPy6iwcE/88O2TIPNLIFwR29wyUlA+YTZOq7dT9ebVTNl/ACOnnc4/wX95c3InT6zFS9zSy1aChaOU7nN8FYU/4TRAF7iaZSSbMgR96cELWW6DYYcMe5fiaZQH0ZnP2b+iloag0rZkkcgXhaSloMmhYJT6UQa/8/U1eSoNUA7dy/ShCCHzh3lfs76yh1CHvf/s4AMBn3j3DFiheD/OC13fY9m33C6PMFJNGNnuqwgDw1oHi7GbHjN8QGkpBCCQD1vIylavqq0a+ppdd11BKopRsgeL/AJtmGwft1JeTo+CO5jS6fX0o3pqW6mj+9rRp3lFeDh9Kab5v7be0x1Hy/g8OlGpJTg3F+x4A4MhxIwAAq7Z3a+cExo6QHZ3HyBkoPhO2yctrhJ4v2CNkd2d0yrRO+7PXM6Js/iNb0t4Oa6WhuMxB+XxRO04QGTUsPcorkoYiz6mJjZ7zUHx8MIDVySs/mNf74NYQsh5aZ5Oh/sXllwhNqYTnoEP/TbzyUHXqiChQ3L8bayh1yLc/egKe+bczkUiQPRJTD57+Atz6qVO15e2dP+y00a0477jD8Ob+fry4cW/JKHtf7yAGcwWcMHkkJo9q8Xwx2jJJ361Fk9KOberT+x1mitIRepBTvmcwj8mjrF0hN+5xTrpSHUyQyTeb8xaOjglnHsXn8tYouL05hQN93iYvZVJSfg59bocajR4xboRnlJcSaOmEd2fx5V+/BgBYs9MSAt4CRbtPjzZUQmaq3FmzR8sjlxd223qNnIFih+pnmszmC3Zn6BZMKv2MsW2eZkXVAY5qSXsKtYI0q2ZSiZKVm5UPJZMk4/Ojz0MxjZz159tk8lIaipevbyCrO+VL69GfzaO9OY1kgjxNXuq9U3NJvDTqplTSs/7qfIIILZmkp8AaCClQ1MCtXMe8OxCA56HUIZNHteDwMdY+7Wr01zNg/XD6KGTOcRPs5e3d8zR6BvKOGcW7ugcc50/61gIAQHtzCk2usEz1Yoxtbypx5ioKAprvwKShmM0BuYJAKmlFQJnS9wzkML69GZ2taRx0+TH8Nh5yl+P92f9FU5MSO1rSxigvlV2TPYIv5qk6j+Z0AukkGU12TWnvzk4NAFRH0+MhUPQ8sx55qHzV6FPvMHOFAppSSSQTZOxEdKewdb/eZtd2g4aiOpvWTNJbQ5HP8rj2JvQM5DzCZgWSRMikko7VhnOaD6UpnTTWvzhpsDQk112HlnQSfYPONlb3b9KSAeDggNamHn6+vqz1HjalEp7L3Kg2szWUXOlv2pz2jlJTaZtSCVl/L4FSPOYZuq00lOZoGsqC13c4vnOUV52jHjQ1Qt15wBIMnzp9GoioaPJy/ZD7+7IYKVVt67u3YPjYqVOQSSUcoxv14B99WDu2dPVhf29ph1ooFJ3RJhu8nqf7pc8VCkglCKmkt0DpHcyhdzCPse0ZtGZSxol5wfNQvDUUfST4td8sw5Ff+4OrfkIzeXlrKMUOt3QEq+69KZVESzqJ3gFn/fWRqZ99XL3o31+wpuSagwM5uzP3GqXvk4OBdunj6NU6zJ6BPFpkR2cyeSkBpJ5Brz55MFewJ9a5f2P1m7Vmkp5pV0oT3Pj2ZhQEcNDVoVt+OqA5lXBoAmoeCgDf+vdn80gmCM0ps9DpG8wjQZafyaShjJYT/tyDGgA4OFB8N7xWm84XrA26mtNJz3qqZ0a1oZfpuSmVNJjb5DOWls+Yl4aiCTFPLXHQqaGU65R3LwoatBRSpWGBUiZKy3hp014AwE//tB5AcbKdcojqZpmBXB4HB3IY05bB7/7pPQCArp7ig//p24tbyzank45RXqEg8A8/XwoAOP2IMQBKzU2L1+9B90DOHpl7qblPrtqB2xatt7+XjNClHTydSHgKJCU4J7Q3oyWTLDHLKFtzkIrtHJUXrx3MFdAm2/bR5W+VvLDZfAGpRAIdLSkcHMh5dtg5bQQJODsUJZSb0wmMaEqVmKxUBzmiOWXNB3KbxIRTu1i/q3SdpT++vgOj2ryDFoBi0MXEkS3IJBPYqkX17T44gHHtTSXaqc7XHlwGoNjZGH0ohtFtXzaPTDIhO9PSMn73qrVtbbP8HRav73KmH7S23x3RnHII9Vy+YJuC3ZMedbr7c2hvTlnPj0Ho/PCpdSgIK5LL3SEr7aOzLYNkgjx9afrv6v4NVGfdnE6i2aihOIW211IuTQYNZUDTUFqbkp5arH7Mq53UKgleTvl7X9yM59fvKUmj0+0qs2fAu52rBQuUMiE51+OJ13dgf18W97xohWJ+5byjAABTR7cCAFa9dcBO8/qb1ufOtgw6ZYejJrbt78viT2t329c2pRJWpyJfpidef8s+N12a3dwmn0/c9rxMa3UEXqOa37zs3OPa3dn0DlqmgGTSe9LWyu3WPUzoaEZzutThqNJ4jRp11u08qKUp1mEwX8C49ibHte7lYNJJax4K4O3DUFqIMkvmCqUaSnMqibamFHpco2/1crd7OGO3dPXaA4ZPnja1JE8AeF4K9S1dfUgmCHsOlmqgy+Sij21NSXS0FDWtQkHYAqWtKWXUwFa9ZWkQdui2a4j78ua9DsE84NGhNqcts2G3T6TcBiksP3/Py47jPYM5tGaS6GhOY1/voC10+7IFtEhTYJNLu9bp7s+iozltNAfptGZKr9m0pwdj2jIY0ZRCe3PpoMAqI2ffv1vTVIstTh7Vgua0t49jU1cvWtJJWwvStdycFuXlJ1AyKTlo8fgdH3ipOH1A16YU/S4NReXZN5jHVx9Y5hh8uukbzOP/nlnvOLZ8237M+cEi7DzgvZ1BpWGBEoNfvbDZ/qwidyaPasGEjib75R/MFfDXt/wFgPWSqBh45Qvpcs2YzqScI0h94b3D5KRJ/UHVBYPdkXpoCZ2auc2dDgD29Q2iszWDTDLh6RReuGon2ptTOPXwzpIOQZ9f0T2Q81lvTOD2P2+wv+vaSjZfwPj2Zsf1zsAEqaFIDdDd6RYKAle7RvC6FmMLlLQUKG4NxW0719pw8YbiSP1dR47FdRcdX1IHPb+jD2u3f38TmSTZ9vn9fVlk8wLjRjRh7Igm7D444Jt2amdrSZkA7OdMDXoGXB1mf9Yyq4308EPpwvv8dxwGAHjfzHGOa3oH8mjLpHDk+BHoGcxjm9SwDvZnbUE8dXQr1u086Bl6fUDTUPqy+ZJr1G80cWQzWtMph0kQAHZ2D9jvgKnD7u4vhgW7NZTtctOxw8e0SrOyl1kyKxDHoFUAABobSURBVM26pSavfttkmPIOWbbNqt5aMFBccoeo9BkGgJ8v3gSgVENRpiy/PeLX7z5YcmzJpr1YvaMbv39tuzFdJWGBEgO949R3rZs2uhVb5IZbi9bsso9PH9OG1kwSmVQCXVKgfPeJ1Y48U4mEw+yhNIHHvvBeW2DoD+IurfNpa0ohmSh1OAPA1r3O9a9+qQnDx5a/hS1dfRjZksbotgz2uITco8u24/6lWzFz/Ajb/tyXLTVdqVGd10z2XL6AW55+w3Fsn+YLyuYKGNfh1FD0EeRgXiCdStgmJ3eU3OauXtsh6dWhKOHVkkliRFPS8bJv2tODv6yztERlstSdsXo9kgmyR8BeJg0AmD62DZv2mJceb04lkU4VTYt/XLnDLnvsiKaSgA13HaaPtTRV/TnQTUhnHGmZRt2aal82j5Z00p7Lo3fo6vf86nlH45RpnUglCEeOb3Ok39s7iJEtaUzsaLa/v7Z1H17avM+ebHjcpJHYfXDQ8VwquvuzaG9OoTlt+XDcJp+LfvQsAOCyd0+3hI5LQ9nVPWBrse3N6ZIJrrl8AWt2dGOyjKJzD6zUKH1ce5OnD0UIgd+8vA3d/TnbF6oPvJQJrr05BSFK56IpM3dHc9ooUJpSCVx8yhSM8NBECwWBTXusfkO1p6rjjgPO9rxvyRZ8+dfO/VIefGkbAOD8d0wsKbdFCwiqJixQIvC/c08CYN5mdWpnq92Bb5KCBQBOntYJIkJnaxort3dj8fo9eGz5W460yQQ5TF7qoexszdgOYd1OqquyTSkrgqnLIxJM7Sj59gmWJrVmhzWCFqLoo+loTmFce+kI+XO/eAlAsVNe/VY3Xtmyz+6QlIN7ohw9ek08fPaNPfjvxy3h+Tezp6E1k7TNfkII9OcKmODSUHTBdqAvi5Etaduh7R5h6y/vqJY0UgmyBdYbuw7ifxeuBQCMacugrSmF/mzBHhC8/7+fxoMvWy/jCFtDKXYk+v2kEmSPXtXM80de247L5y8BANw+bxamj7F+fz0P1VbtzSmcdcx4pJMJO2poi3xWTpg8EhNHNmPVW92O+3t1yz4c/Y3ienFK21Qd0sbdPTjq34vn/+odE0GEksADy+RlaSjZvHBoiGq0P6I5BSLrHt1O8V3dAxjX0YQJUqCs3XEQn/v5S472UsJQn8irsHwoaTus221yUlodwVy+2q+9vSlVYjJavKELvYN5e76NbhI80J/FK1v2IZUgjG7NYGRL6QRZZRLb15u1VxvQhXKftpIAUCoQ1W82sjWNER4muWy+gIFcAdNGt3oGl+gRnCrARgk0XXh8/Nbn8JX7X8P9S7faZezvzdra/1fmHIVvfuRYnHvsBDvNtr3mBVUrCQuUCHzg6PEAgLuf3+R5fsroVry5vw+DuQK27rU68tf+41z7fGdrBovW7MInbnsen5w9DQAwfYxlxkiQZRrb0zOIbL5gv+jtzSk78qRLm2Gtj1yOGNeG/mwBv1y8GY8tt1TcQkHg/Jv+hHU7D+JzZx6JJ774fgDWXATAaXZKJhLGETIAvHfmWACW6QEovlAvbuyS99Am61QqaPU8r7vwOHS2ZmzB99TqncgXBMaMcC7XvXlPURjv7R3EqJa0LVSVX8qqRx6Xz3/R/j5reidGt2XwyLLt2Nndj0tufQ6vbLFmiY8ZkSmGfg+WjiCVSVLvDPSFHIkIRx/WAQD43Wtv4qnVO3HlL1+yz59x5BjMHN+OXEFgtWb2+u0rlsD6tzlHoSmVtASKbL/+bB6ZVAJHH9aOE6aMBAD81c1/ttN+9YHX7M+3z5tla1GqjkrDAYBPnmZFGx4+utWxQdM//+plPPH6DrRkkpjSWTqX6MbHLGGvTFetmVSJQNrZPYDx7U04blIHMskE1uzsxjTpM5z7Tus5njneGrC4NbS7n9uIVW912yYvwBn2+3/PFLXX898x0XLKa8/mWd97Gm8d6LfXymtvLh3h/+1PLf/C+2aORTJBWK/d/zd+uxy/feVNdLSkkZBLw+t+rp3d/fivR1fZ3ztarHbQhU5XTxYJgi3U3Ca5A31ZEAEjMil7zTldi1G/x+g2a2DkjvR8Xfop539mNjqltq/epc3awPSFjUUT7LKt+/HKln34zuPFuh82shmfec8MhzD84VPrjPODKklDCxQiOo+IVhPROiK6aqjK7WhO49iJHfb3Z/7tTMf56WNaIYTlY9nS1YejJrTbHSHg7KwKQqAtk7RXMU4kCGcdPR77+7JYvL4L3f3Wctst6aQ9afJHT72Bqx54DX2Dedu2+ujn34szjxpv56sc/W8d6McK2fm+922WQGhvSuFnz26EEMIxEv74O6dg7IgMtu3rs0fUA7k8iCxh8pXzjgYA/Os5bwdQNPn83V3W6PzYSR1IELBobdHMl80XcP5Nf7JHWC9/4xykkgnsPjiAB1/aBiEEPnOnlf6t/f342KlTcIQUdqrDy+YL2NdraSiHS8H7n4+stNvx679ZbgvWBz73LoxqtaKANu3pxZwfLHL4qUY0pWw/iar/MdpvOXW01dnetmg9npHmyn2ukez0sW348AmHIZsTuOxnLzrOtWYsPxMAfOTmP+OJFW/hmw8txxfvte7/+MmWwNDnwjy9eifGtzfZS8sAwKY9vXaYsb6CwoyxbZrp06qXLqw/duoUAJZPT+9QVQTXhPZmO3BEdVbrdx20ncXKpNTalHQI3J6BHA7KeUiJBGHK6BZs3tOL52TU0XnHW36Xj54y2S5fIYTANx5aYbVlb9bWUJTA6M/mHZ35pFEtllNeCpz9fVk7qk6ZuUa2pLHizQOYftUj6BnI2SZLADhhykicMm0Ulm7cax976BXr/tWzMGZEBnu0gdm1v3vdNpl+75ITbaGxW7btht09uPWZN3DspA7b7KpHUN3zwmbc9OQ6a1+iBGHSyGbkZLCF4uXNVn2OOqwDU0e3OoQEALy21QraOHnaKEwa2YzRbRm8tnW/7wKRT6/ZiYt+9Cx+sdgyYc85boIdnLPTZSZzb5lQDRpWoBBREsCPAHwIwLEA/oaIjh2q8n975bvtz1Okk1Rx4tRRAIBrHl6BP67cUTIZUY3iAOAXizdjxrg2nCrV9HEjmuzw4E/dvhg/fGodiMh2tCrueXELjvnmY3hgqdURvH2CtQrygi++D4A1+lu0ZpcjvFXZ1o+ZZHWg31+wxnZU/uTSWZg4sgWHjbQ6VOXvWL5tP4SwRqDKrqxMW998eIXtmAWspWlOmDIKz72xB3t7BtGfzWOHJtA+MWuqPfJSr8gLDof3GHz3khOx8EvvR1smiWXbrJdp5tcfBWCNIts1wfx72UmqFxWAPWJWwkJ/iWZPHw0isjvkrp5B/GntLjuCTb+3e17cgnl3vIAdB/rtjv3pL59pXzfr8NGOjl5nvOYLuuLupbjruaIme/Rh1u/UnE7iqdW7sHj9HnT1ZG3NSK1mDViTXX/01DocJX/b8e1NmNzZghFNKWRSCfzs2Y0oFIRt2vurEyfh+EmWwGptSmHl9gPoGcjZ/jwA+PQZh9tCY9Ga3dhzcMAWCgAwe8ZoANagSQn09bsO4twfLHLc44T2ZjzqMtcClgb3wWPG48lVO5HNF7Dg9R0O/93ugwO2sHnuDatcZX4FgMe/YD2/bRnLZFQoCDyxoliOElyd2uZTG3b34JNSOzn6sHa0N6cxbXQbNuzuQVfPoMPX9f2PnwjA0jL6swXs78tiMFfAI9Jp/aHjD8PFp06xBYoyu6qQ7bZMChPk73vvEqsTF0LgKnleMUamV1r4/r4svvrAMmRSCcw6vBMzxrZh455eW4PZvr/PNgl3NKdBRDhp6igs3tCF6x9ZCcDSPgHgg8dMwAtfOxvtzamSqK7/+/Qs+/NP51mfP/u+I6y69Jg3pqsUqeBL6pbZANYJIdYDABHdA+BCAKXG2yqQSSVw5Lg2vLGrp2RjqxljnM7Mo7URMAB895IT8dEf/8V2wB0/aSSu/vDRmDt7qj16HN2WsR8A3fn/yjfPsWfVA8CrclSj6jBzQjuOHNeGB1/ahgdf2mYLtz995QO2ULr2guPwof/9E25+ch1ufnIdgOLI/BPvnIpbn3kD//34avsBB6zOXqFMU4+8tt1+EX9++WloySTtsk++rlhHADjvuMPw9Y8cY3//1d+fhot//JzdEXz/4yfiQydYzkQiwplHj8eDL23Dsq3F/dWVhnTjxSfgqw8sw1UPLnO8yEeMbbM7yyPGjbDNXADwj2ceiS+da4V2K3Pf+TcVzUqAFaE3bbTztzvt2wsBWL+HcoYDsM1GimMmduD/PnUqgGL4tpunvnymvXDgSVNH4YUNXXbI96dOtzqLw8e04a9OnGRrFLrf6b8+eoKd12CugO37+3GEnABKBNz8Nyfb51fLsPXjrnkcl79nRvF+Zoy2I5Tu/MtG3PmXjfa5Z686yx40HDOxA796YTOmX/WI4x5UPd915BhbEKmoN8WZR43HH1futAcCOtdecByOn9yB9qYU/v23y/Hvv13uOH+UFLhTR7diMFew7w8Anvji++yBkxLMgKUJKm6THeqo1jT29AziFO05/OuTJ+Ojp1ganFr54sRrn8Dcd1qh4J+YNRXfuug4AMVovxseXYUbNO3p2guPs027P3rqDQzmCnh8RdHk+Pfvtdpa+VnO+58/4UvnvB33Sn/bO6d3IpEgzBjbhsFcAf/5yEq8Y8pI3CR9fKdMG2XnNWt6J55ctRM/lb6Rv3nnNHzzI8faAUB//94jHBNs773idEdbThrVgo03nI+BXB5fPOftjsChatGwGgqAyQC2aN+3ymNDxmNfeJ/DN6JIJAjPX322/f3Gi09wnB8zogkPX2lNcBzRlMK/zTkK6WTCflkAq/NR3KM9KKNaM9jwXx/GtRccZx/7xzOPdOSvLwL46pZ9aMsk7XWiAKuz+Omlsxxp1Mh+RFMKX5lzlONcOkmOEeG7jhxbcs+zpltl/r93TS85B1gvom72O2lqJ44+rN0WlmdoAgsALjhxEgBgrZy3cudl77RHjR+fNdU2Kym+9uGjsfBL77e/f2XOUY66XHzqFE3ojoCb3//ze/DsVWchk0rgha+dXXL+jCOc9XPX9xvnH4NpY4qa6tUfOtpx/or3HWELMgC4RJqmFFM1LfemuSfZ/irFhSdN8v2um5gA4OoPF4W3ctYu+49zkUom7CVe3IzX5gFdesbhJed/+Xen2RriP589E1/44Exc9u7p+PTpzmv/9rRp+MBR40rSr73+Q3Zgyuc/ONNx7qSpo7D82jmO72709+M9M0ufwQc+d4b9G3zwmAmOc6Na07jx4nfY308/YrT9Wc0l+6ez3uYYDHzxg2935PGjT56Cow/rQHM6iYtk+//kTxts09WNF5+Aqz9ktbv+jH1vwRps3duHi06ahB98/CRH/e54dgO+cO8rWL+7BxefMgW//od32ek+fHwxWuvv3zsDx0/ucAiFfzl7Jr724aNx5QeOxMYbzsdprmdU0ZRKDokwAQCKuoFLrSGiSwDMEUL8nfz+aQCzhRD/7LruCgBXAMC0adNO3bTJ25FeDV7c2IV9vVmcc+wEz/O9gzl7xBqFlzfvxei2jD3aUgzmCvjZsxuQTiZwz4ubcf1fn4B3Th9dkv7gQA53P7cJ5x1/mKOzA+Ts56S1jEYmmSgxuQHAnoMD6B3MY0pni+O8EAK3LVqP17btx4lTRuLSM6Z7PtCFgsBdz23ESdM6PTuQTXt6cN+SLRjRlMZn33eEY3n+7v4sNnf1YtGa3Zgxts02hbh5Zcs+pJOE46QpSG+jdJKwZNNevH1Cuz2i1NnVPYCzvvs0/vEDb8Nl7y69h8FcAfct2YKNu3vw9fOPKWmjXL6AJ1ftxI7uAXzslCkloZu7ugewcU8PXtjQhcvfM8OzjfoG81i9o7ukfQoFgYGclf++vkF8fNZUW7tQWL6L5fj585vxk0tnOZ7DfEFg054efPsPK7HizQP43sdPLBkorNtp+VbOOXYCEtIEE5Z8QWDxhj0Y1ZLBLxZvwpfOPcoOK1d129MziIUrd2BMWxPOPmZ8SftZk2CtTcGU70nnqVU70ZROYMPuHhw/aaStjet1WL/roGXeyuUxcWRLSR7b9vVh0Zpd6M/mcdm7Z5SczxcEHnxpK6aNbnV02Fu6enH7nzdgRFMK5xw7AWPbmxyDNsCa7PrY8rfwgaPHI18o4ANHOe/xlS37sGjNLuw+OICjDmvHJadOtRdo1csXQtj+01pBREuFELMCr2tggXIGgP8QQsyR368GACHEf5nSzJo1SyxZsmSIasgwDDM8CCtQGtnk9SKAmUQ0g4gyAOYCeLjGdWIYhjlkaVinvBAiR0T/BOBxAEkAdwghVtS4WgzDMIcsDStQAEAI8QcAfwi8kGEYhqk6jWzyYhiGYeoIFigMwzBMRWCBwjAMw1QEFigMwzBMRWCBwjAMw1SEhp3YGAUi6gPgF1o8EsB+n/PTAGz2OR+UPuh8JfKodh25fo1dv0rU4VBvw+FeP69rjhJCtJsuthFCHDJ/AHYFnL+tyul9z1cij2rXkevX2PVrhDpy/WpbP69rACwJSiOEOORMXvsCzv+uyumDzlcij2rXketX2/Nx61eJOhzqbTjc6xf2mhIONZPXEhFiPZpqpR8K6r2OXL941Hv9gPqvI9evfMLW6VDTUG6rcfqhoN7ryPWLR73XD6j/OnL9yidUnQ4pDYVhGIapHoeahsIwDMNUiUNeoBDRHUS0k4iWa8dOJKLniGgZEf2OiDrk8TQRzZfHV6o9WOS5p4loNRG9Iv/G16B+GSL6mTz+KhGdqaU5VR5fR0Q3kdeOWbWtX7XabyoRPSV/rxVE9Hl5fDQRLSCitfJ/p5bmatlOq4lojna84m1Y4frVRRsS0Rh5/UEi+qErr5q3YUD9Kt6GEep3DhEtle20lIjO0vKqyntcMcKEgg3nPwDvA3AKgOXasRcBvF9+/gyA6+TnTwK4R35uBbARwHT5/WkAs2pcvysB/Ex+Hg9gKYCE/P4CgDMAEIBHAXyozupXrfabCOAU+bkdwBoAxwL4DoCr5PGrANwoPx8L4FUATQBmAHgDQLJabVjh+tVLG7YBeA+AfwDwQ1de9dCGfvWreBtGqN/JACbJz8cD2FbN9qvk3yGvoQghFgHoch0+CsAi+XkBgIvV5QDaiCgFoAXAIIADdVS/YwEslOl2wgo/nEVEEwF0CCGeE9ZTeReAi+qlfpWoh0/9tgshXpKfuwGsBDAZwIUA5svL5qPYHhfCGjQMCCE2AFgHYHa12rBS9Ytbj0rWUQjRI4T4M4B+PZ96aUNT/apFhPq9LIR4Ux5fAaCZiJqq+R5XikNeoBhYDuAC+fkSAFPl5/sB9ADYDmsm63eFEHpn+jOpJn+jyqqoqX6vAriQiFJENAPAqfLcZABbtfRb5bF6qZ+iqu1HRNNhjf4WA5gghNgOWC88LI0JsNpli5ZMtVXV2zBm/RT10IYm6qUNg6haG0ao38UAXhZCDGDo3+OyYYHizWcAXElES2GpqIPy+GwAeQCTYJkbvkRER8hzfyuEOAHAe+Xfp2tQvztgPWRLAPwPgL8AyMFSj91UM7yv3PoBVW4/IhoB4AEAXxBC+GmVpraqahtWoH5A/bShMQuPY7VoQz+q1obl1o+IjgNwI4DPqkMel9VVmC4LFA+EEKuEEOcKIU4F8CtYdmrA8qE8JoTISpPNs5AmGyHENvm/G8AvUUUzhKl+QoicEOKLQoiThBAXAhgFYC2sTnyKlsUUAG+6861h/arafkSUhvUi/0II8aA8vEOaEJQpZqc8vhVOrUm1VdXasEL1q6c2NFEvbWikWm1Ybv2IaAqA3wC4VAih+p8hfY+jwALFAxXZQUQJAP8O4FZ5ajOAs8iiDcDpAFZJE85YmSYN4COwzD5DWj8iapX1AhGdAyAnhHhdqtPdRHS6VOEvBfBQvdSvmu0n7/d2ACuFEN/XTj0MYJ78PA/F9ngYwFxps54BYCaAF6rVhpWqX521oSd11IamfKrShuXWj4hGAXgEwNVCiGfVxUP9Hkei0l7+RvuDNYLeDiALawRwOYDPw4rEWAPgBhQngI4A8GtYjrLXAfybPN4GK2LpNXnufyEjb4a4ftMBrIbl9PsjgMO1fGbBejneAPBDlaYe6lfl9nsPLLPAawBekX8fBjAGVoDAWvl/tJbm67KdVkOLoqlGG1aqfnXYhhthBWsclM/FsXXWhiX1q1Yblls/WIOwHu3aVwCMr+Z7XKk/ninPMAzDVAQ2eTEMwzAVgQUKwzAMUxFYoDAMwzAVgQUKwzAMUxFYoDAMwzAVgQUKw9QJRPQPRHRpGddPJ22VZ4apNalaV4BhGGtSnRDi1uArGaZ+YYHCMBVCLvz3GKyF/06GNbHzUgDHAPg+rImxuwH8PyHEdiJ6GtZ6Zu8G8DARtQM4KIT4LhGdBGuFgVZYk9g+I4TYS0SnwloTrRfAn4fu7hgmGDZ5MUxlOQrAbUKId8Da2uBKADcD+Jiw1ja7A8D12vWjhBDvF0J8z5XPXQC+KvNZBuAaefxnAP5FCHFGNW+CYaLAGgrDVJYtorj+0s8BfA3WJkkL5EroSVhL1SjudWdARCNhCZpn5KH5AH7tcfxuAB+q/C0wTDRYoDBMZXGvZdQNYIWPRtFTRt7kkT/D1A1s8mKYyjKNiJTw+BsAzwMYp44RUVruc2FECLEfwF4ieq889GkAzwgh9gHYT0Tvkcf/tvLVZ5josIbCMJVlJYB5RPR/sFaRvRnA4wBukiarFKzNxVYE5DMPwK1E1ApgPYDL5PHLANxBRL0yX4apG3i1YYapEDLK6/dCiONrXBWGqQls8mIYhmEqAmsoDMMwTEVgDYVhGIapCCxQGIZhmIrAAoVhGIapCCxQGIZhmIrAAoVhGIapCCxQGIZhmIrw/wHhuQzplAiyiQAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'].plot()" ] @@ -215,9 +2245,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xt8XHWZ+PHPMzO535qk6S1taekN2gIFagFBQQptUaC4wlpUqCuKIrr603UVd/3hovxWvKHoooJUC4sKgkrlaikKFNpCCxR6pek9TZumuTXXmczM8/tjzqSTdDIzSTPJzPR5v17zmsl3zjk5c5rOc77P9yaqijHGGJMI13CfgDHGmPRhQcMYY0zCLGgYY4xJmAUNY4wxCbOgYYwxJmEWNIwxxiTMgoYxxpiEWdAwxhiTMAsaxhhjEuYZ7hMYbCNHjtRJkyYN92kYY0xa2bBhwxFVrYi3XcYFjUmTJrF+/frhPg1jjEkrIrI3ke0sPWWMMSZhFjSMMcYkLOGgISJuEXlTRJ50fi4TkZUissN5Lo3Y9jYRqRKR7SKyMKL8XBF5x3nvHhERpzxHRB5xyteJyKSIfZY6v2OHiCwdjA9tjDFmYPpT0/gSsDXi528Aq1R1GrDK+RkRmQksAWYBi4B7RcTt7PML4GZgmvNY5JTfBDSq6lTgbuAu51hlwO3AecA84PbI4GSMMWZoJRQ0RGQ88CHg1xHFi4HlzuvlwDUR5X9QVa+q7gaqgHkiMhYoVtU1GlrE48Fe+4SP9Rgw36mFLARWqmqDqjYCKzkWaIwxxgyxRGsaPwH+HQhGlI1W1YMAzvMop7wS2B+xXbVTVum87l3eYx9V9QPNQHmMY/UgIjeLyHoRWV9XV5fgRzLGGNNfcYOGiFwJHFbVDQkeU6KUaYzyge5zrED1PlWdq6pzKyridjM2xhgzQInUNC4ErhaRPcAfgEtF5H+BWiflhPN82Nm+GpgQsf94oMYpHx+lvMc+IuIBSoCGGMdKK6rK4xuqafP6h/tUjDHmhMQNGqp6m6qOV9VJhBq4X1DVTwArgHBvpqXAE87rFcASp0fUZEIN3q85KawWETnfaa+4sdc+4WNd6/wOBZ4DFohIqdMAvsApSys769r46h838rcth4b7VIwx5oScyIjw7wGPishNwD7gOgBV3SwijwJbAD9wq6oGnH1uAX4L5AHPOA+AB4CHRKSKUA1jiXOsBhH5DvC6s90dqtpwAuc8LOpbvQC0+wJxtjTGmNTWr6Chqv8A/uG8rgfm97HdncCdUcrXA7OjlHfiBJ0o7y0DlvXnPFNNQ5sPgM6uYJwtjTEmtdmI8CHQ0B4OGlbTMMakNwsaQ6DRqWl4LWgYY9KcBY0h0NDWBUCn39JTxpj0ZkFjCDRaesoYkyEsaAyB+jYLGsaYzGBBYwg0Wu8pY0yGsKAxBBqspmGMyRAWNIZAd5uGNYQbY9KcBY0k6+wKdI8Et5qGMSbdWdBIsnBqCmychjEm/VnQSLJw0HC7xBrCjTFpz4JGkoXbM8YU59Lpt5qGMSa9WdBIsnBNY9yIXGvTMMakPQsaSRYOGmNL8iw9ZYxJexY0kqyxzYdLYHRxjtU0jDFpz4JGkjW0+xiRn01etgevP0hoQUJjjElPcYOGiOSKyGsislFENovIfznl3xaRAyLylvP4YMQ+t4lIlYhsF5GFEeXnisg7znv3OMu+4iwN+4hTvk5EJkXss1REdjiPpaSZxvYuRuRnkZsVutReG+BnjEljiazc5wUuVdVWEckCVotIeJnWu1X1h5Ebi8hMQsu1zgLGAc+LyHRnyddfADcDa4GngUWElny9CWhU1akisgS4C/ioiJQBtwNzAQU2iMgKVW08sY89dLxdAfKy3OR63EBogF9ulnuYz8oYYwYmbk1DQ1qdH7OcR6wcy2LgD6rqVdXdQBUwT0TGAsWqukZDOZoHgWsi9lnuvH4MmO/UQhYCK1W1wQkUKwkFmrTh9QfJ9ri6A4U1hhtj0llCbRoi4haRt4DDhL7E1zlvfUFE3haRZSJS6pRVAvsjdq92yiqd173Le+yjqn6gGSiPcaze53eziKwXkfV1dXWJfKQh4/MHyXa7utNT1hhujElnCQUNVQ2o6hxgPKFaw2xCqaYpwBzgIPAjZ3OJdogY5QPdJ/L87lPVuao6t6KiIuZnGWq+QKimkRNOT9kAP2NMGutX7ylVbQL+ASxS1VonmASB+4F5zmbVwISI3cYDNU75+CjlPfYREQ9QAjTEOFbaOL6mYekpY0z6SqT3VIWIjHBe5wGXAducNoqwDwObnNcrgCVOj6jJwDTgNVU9CLSIyPlOe8WNwBMR+4R7Rl0LvOC0ezwHLBCRUif9tcApSxtdgd5tGlbTMMakr0R6T40FlouIm1CQeVRVnxSRh0RkDqF00R7gswCqullEHgW2AH7gVqfnFMAtwG+BPEK9psK9sB4AHhKRKkI1jCXOsRpE5DvA6852d6hqwwl83iHn624ItzYNY0z6ixs0VPVt4Owo5TfE2OdO4M4o5euB2VHKO4Hr+jjWMmBZvPNMVeH0VHebhqWnjDFpzEaEJ5kvECQrIj3ltYZwY0was6CRZF7rcmuMySAWNJKsKxAkxwb3GWMyhAWNJPMdNyLcahrGmPRlQSOJ/IEgQYUst4tcj43TMMakPwsaSeQLhAJEtseFx+3C4xIbEW6MSWsWNJLI50yDnu0OXebcLLelp4wxac2CRhJF1jQAcrNclp4yxqQ1CxpJ1F3TcIJGjseN12oaxpg0ZkEjiY5PT7ls5T5jTFqzoJFEx6enrE3DGJPeLGgkUdSGcOs9ZYxJYxY0kqjLGsKNMRnGgkYSeXs1hOd6LD1ljElvFjSSKJyeyrJxGsaYDGFBI4nCQSMn3OXW0lPGmDSXyHKvuSLymohsFJHNIvJfTnmZiKwUkR3Oc2nEPreJSJWIbBeRhRHl54rIO8579zjLvuIsDfuIU75ORCZF7LPU+R07RGQpaaR376kcj7u7zBhj0lEiNQ0vcKmqngXMARaJyPnAN4BVqjoNWOX8jIjMJLRc6yxgEXCvs1QswC+AmwmtGz7NeR/gJqBRVacCdwN3OccqA24HzgPmAbdHBqdU190Q7qSnstzSXWaMMekobtDQkFbnxyznocBiYLlTvhy4xnm9GPiDqnpVdTdQBcwTkbFAsaquUVUFHuy1T/hYjwHznVrIQmClqjaoaiOwkmOBJuV1t2l4wkHDRZcN7jPGpLGE2jRExC0ibwGHCX2JrwNGq+pBAOd5lLN5JbA/Yvdqp6zSed27vMc+quoHmoHyGMdKC73HaWS5XXQFdThPyRhjTkhCQUNVA6o6BxhPqNYwO8bmEu0QMcoHus+xXyhys4isF5H1dXV1MU5taPXucpvtpKdCFS1jjEk//eo9papNwD8IpYhqnZQTzvNhZ7NqYELEbuOBGqd8fJTyHvuIiAcoARpiHKv3ed2nqnNVdW5FRUV/PlJShRu9w72nPG4XqhCw2oYxJk0l0nuqQkRGOK/zgMuAbcAKINybaSnwhPN6BbDE6RE1mVCD92tOCqtFRM532itu7LVP+FjXAi847R7PAQtEpNRpAF/glKWFLn8oOGRFpKcAugIWNIwx6cmTwDZjgeVODygX8KiqPikia4BHReQmYB9wHYCqbhaRR4EtgB+4VVXDI9puAX4L5AHPOA+AB4CHRKSKUA1jiXOsBhH5DvC6s90dqtpwIh94KPkCAdwuwe0KZdmy3KHnrmCQPNyxdjXGmJQUN2io6tvA2VHK64H5fexzJ3BnlPL1wHHtIaraiRN0ory3DFgW7zxTkc8f7G4Eh2NtG9aDyhiTrmxEeBL5/MHuQAGWnjLGpD8LGknkC/QMGh4nTWUD/Iwx6cqCRhL5/Bo1PWVTiRhj0pUFjSTqXdMIp6f8lp4yxqQpCxpJ5PMHetQ0jrVpWE3DGJOeLGgkUe+GcI/T5dbSU8aYdGVBI4l6p6fCtQ7rcmuMSVcWNJKoy6/dA/ogok3DphExxqQpCxpJ5A0EyfYcG/mdZekpY0yas6CRRL1HhGdZesoYk+YsaCSRzx/onuEWbES4MSb9WdBIouPHadiIcGNMerOgkUR9NYRb0DDGpCsLGkl0XJdbj6WnjDHpzYJGEoUawo/1nrIJC40x6c6CRhIdNzW6x9JTxpj0lshyrxNE5O8islVENovIl5zyb4vIARF5y3l8MGKf20SkSkS2i8jCiPJzReQd5717nGVfcZaGfcQpXycikyL2WSoiO5zHUtKEqvY9ItzSU8aYNJXIcq9+4Kuq+oaIFAEbRGSl897dqvrDyI1FZCah5VpnAeOA50VkurPk6y+Am4G1wNPAIkJLvt4ENKrqVBFZAtwFfFREyoDbgbmAOr97hao2ntjHTr5wYMi2hnBjTAaJW9NQ1YOq+obzugXYClTG2GUx8AdV9arqbqAKmCciY4FiVV2jqgo8CFwTsc9y5/VjwHynFrIQWKmqDU6gWEko0KS88KjvyJqG2yWIWNAwxqSvfrVpOGmjs4F1TtEXRORtEVkmIqVOWSWwP2K3aqes0nndu7zHPqrqB5qB8hjHSnk+Z9R35IhwCNU2bBoRY0y6SjhoiEgh8DjwZVU9SijVNAWYAxwEfhTeNMruGqN8oPtEntvNIrJeRNbX1dXF/BxDpTtoRMw9BaEg0uW3Ng1jTHpKKGiISBahgPGwqv4JQFVrVTWgqkHgfmCes3k1MCFi9/FAjVM+Pkp5j31ExAOUAA0xjtWDqt6nqnNVdW5FRUUiHynpwkEjcnBf+Gd/0Goaxpj0lEjvKQEeALaq6o8jysdGbPZhYJPzegWwxOkRNRmYBrymqgeBFhE53znmjcATEfuEe0ZdC7zgtHs8BywQkVIn/bXAKUt50do0IJSesjYNY0y6SqT31IXADcA7IvKWU/ZN4HoRmUMoXbQH+CyAqm4WkUeBLYR6Xt3q9JwCuAX4LZBHqNfUM075A8BDIlJFqIaxxDlWg4h8B3jd2e4OVW0Y2EcdWuGaRk6UoOGz9JQxJk3FDRqquprobQtPx9jnTuDOKOXrgdlRyjuB6/o41jJgWbzzTDV91zTEahrGmLRlI8KTxNsVqlzl9moIz3K7rE3DGJO2LGgkSWc4PZVl6SljTOawoJEk4ZpGznE1DUtPGWPSlwWNJPE6NY3cKDUNCxrGmHRlQSNJOvusabjw24SFxpg0ZUEjSbx9dbn12DQixpj0ZUEjSbqDRlbvaUSsTcMYk74saCSJ1x9OT/W8xB6XtWkYY9KXBY0k6ezqOz1lizAZY9KVBY0k8foDZHtcOIsTdrMut8aYdGZBI0m8XUFyPcdf3mzrcmuMSWMWNJLE6w8e1wgO4HGLpaeMMWnLgkaSeLsCx7VngDO4z281DWNMerKgkSRef5DcKDWNbLeLLpuw0BiTpixoJInXH6OmYekpY0yasqCRJJ1dwT6DRiCoBIIWOBLV7vNzpNU73KdhjCGx5V4niMjfRWSriGwWkS855WUislJEdjjPpRH73CYiVSKyXUQWRpSfKyLvOO/d4yz7irM07CNO+ToRmRSxz1Lnd+wQkaWkCa8/EDU95XHWDLceVIm765ltXPfLNcN9GsYYEqtp+IGvqurpwPnArSIyE/gGsEpVpwGrnJ9x3lsCzAIWAfeKSPjb8xfAzYTWDZ/mvA9wE9CoqlOBu4G7nGOVAbcD5wHzgNsjg1Mq8/qj1zSy3aEyCxqJ21jdzO4jbd2TQBpjhk/coKGqB1X1Ded1C7AVqAQWA8udzZYD1zivFwN/UFWvqu4GqoB5IjIWKFbVNaqqwIO99gkf6zFgvlMLWQisVNUGVW0EVnIs0KQ0b1fwuBluITS4D7CZbhOkquw83ApAdWP7MJ+NMaZfbRpO2uhsYB0wWlUPQiiwAKOczSqB/RG7VTtllc7r3uU99lFVP9AMlMc4Vsrr9AeOW7UPQtOIgNU0ElV71EuL1w/AvgYLGsYMt4SDhogUAo8DX1bVo7E2jVKmMcoHuk/kud0sIutFZH1dXV2MUxs6oRHhUWoartAlt+nRE7PjcEv36331FjSMGW4JBQ0RySIUMB5W1T85xbVOygnn+bBTXg1MiNh9PFDjlI+PUt5jHxHxACVAQ4xj9aCq96nqXFWdW1FRkchHSjpvnzWNcEO4pacSsaM2lJpyu4T9jR3DfDbGmER6TwnwALBVVX8c8dYKINybaSnwRET5EqdH1GRCDd6vOSmsFhE53znmjb32CR/rWuAFp93jOWCBiJQ6DeALnLKUF6vLLYDfahoJqaprpSQviykVBZaeMiYFJFLTuBC4AbhURN5yHh8EvgdcLiI7gMudn1HVzcCjwBbgWeBWVQ13e7kF+DWhxvGdwDNO+QNAuYhUAV/B6Ymlqg3Ad4DXnccdTllKU9U+u9yGg4alpxJTVdvKtFGFTCwrYL8FDZPiDjR1cNNvX+doZ9dwn0rSeOJtoKqrid62ADC/j33uBO6MUr4emB2lvBO4ro9jLQOWxTvPVOIPKkE9fi0NiOxya+mpROw43MKi2WPIzXLz6s4jqOpx080bkype2XGEVdsOs7XmKOedWj7cp5MUcYOG6b9j64Pb4L4TUd/qpbG9iykVhbhdQrsvQEObj/LCnOE+NWOiOtAUanerb/MN85kkjwWNJAgPQovaEB6uadhMt3GF/+ONKcklz0n17Wtot6BhUlbNSRA0bO6pJAjXNKJ2uQ0HDZt7Kq42Z3xGQbaHiWX5gI3VMKktXNNoaLWgYfrBG6OmkW01jYS1+0LXMT/bzdgReQAcau4czlMyJqZjNY3MnWDTgkYSHGvTOP7yWptG4sJBoyDHQ0G2myy30Nieub1STHoLBpUa56bG0lOmX461aViX2xPR7gulp/Kz3YgIpfnZNLVn7n9Gk96OtHnxOTeMlp4y/RKrppHdPbjP2jTiafMeq2kAlOZn05DBd3AmvdU0hWoZOR6XpadM/8TqcntsGhGracQTWdMAKC3IosnSUyZFhdszZo4rzuibGwsaSRBOT+XG6HJr6an4wjWN/OyImoalp0yKOuDMjXZGZQkNbT6CGdpD0oJGEsSqaYRTVt4uCxrxtPv85HhcuF2h2llpgbVpmNR1oKmDwhwPk8oLCCo0dWRmrdiCRhJ0d7mN0qYRno/KVqGLr83n727PACjNz6KxvStj7+BMeqtp6qByRB7lhdkANGRou4YFjSTormn0kZ5yu4QOCxpxtXsD3e0ZEEpPBYJKS6d/GM/KmOhqmjsYOyKX8oLQjAX1GdqDyoJGEhxr0zg+PQWQl+Wm09JTcbX7AhRkR9Y0QndwjZaiMimovtVHRWFOd00jU8dqWNBIglhdbiHUQN7pt5pGPG0+P/k5xwJvWYFT7begYVJQU3sXI/KzKC+woGH6yesPInJsTEZvOR63tWkk4LiahvOf0RrDTarp7ArQ0RVgRH52999ppg7ws6CRBN6uADkeV5/rPuRlW9BIRJvX36tNIwuAhrbM7JVi0tdRp6dUSV4WWW4XJXlZGTvAL5HlXpeJyGER2RRR9m0ROdBrJb/we7eJSJWIbBeRhRHl54rIO8579zhLvuIsC/uIU75ORCZF7LNURHY4j/BysCnP6w9G7W4blpvlsjaNBLT7Aj17T1lNw6SocPfaEc6NTXlB9kmdnvotsChK+d2qOsd5PA0gIjOBJcAsZ597RST87fkL4GZCa4ZPizjmTUCjqk4F7gbuco5VBtwOnAfMA2531glPeV5/oM/2DAhNmW41jfjafX7yImoaRTkePC7J6NG2Jj2FZyoYkRe6scnkMUVxg4aqvgQkui73YuAPqupV1d2E1gKfJyJjgWJVXaOqCjwIXBOxz3Ln9WPAfKcWshBYqaoNqtoIrCR68Eo53q5gnz2nINSryrrcxtfmDVAQETREhBH52TbTrUk54QARrmkU5Xoytmv4ibRpfEFE3nbSV+EaQCWwP2Kbaqes0nndu7zHPqrqB5qB8hjHSnmd8Woa1uU2rmBQ6egKdE8hElaan0Wj1TRMimmKaNMAKMrNsqDRyy+AKcAc4CDwI6c8Wsuvxigf6D49iMjNIrJeRNbX1dXFOu8h4e0KRh3YF5ab5eoeNW6iC9fECnJ61thKC2z+KZN6mtt7tmkU5nho6czMGvGAgoaq1qpqQFWDwP2E2hwgVBuYELHpeKDGKR8fpbzHPiLiAUoIpcP6Ola087lPVeeq6tyKioqBfKRBFappxE5PWZtGbG3dM9weX9PI1FyxSV9NHT7cLqHQ6bhRbOmpnpw2irAPA+GeVSuAJU6PqMmEGrxfU9WDQIuInO+0V9wIPBGxT7hn1LXAC067x3PAAhEpddJfC5yylNfu6zn9RW+5WS5r04ij3Ru9plFWkG1dbk3Kae7oYkReVnc3+6JcD15/sHtRpkziibeBiPweuAQYKSLVhHo0XSIicwili/YAnwVQ1c0i8iiwBfADt6pq+NvxFkI9sfKAZ5wHwAPAQyJSRaiGscQ5VoOIfAd43dnuDlVNtEF+WHX4AlQU5vT5vk0jEl9fNY2ygmwa20PTTrtc0cfBGDPUmtq7KHFSUxBq0wBo6eyiPMZ3QTqKGzRU9fooxQ/E2P5O4M4o5euB2VHKO4Hr+jjWMmBZvHNMNaEG3DjpKX8AVe1zAODJLrw+eO/rOLIwh0BQaWz3Zdx/RpO+wjWNsHCaqqXTn3F/pzYiPAnafQHysvuOx7lZblRtIaZY2rzRaxoVRaH/gEcydIoGk56a2ru6e05BKD0FZGS7hgWNJOiI06YR7o7b6bOg0ZcOX/Q2jZHOXVtdS2ZO0WDSU1OHjxHOLMwQkZ7yZl77mwWNQaaqoZHMMQb3hUc520y3fWsLB40+axoWNEzqsJqGGTCvP0hQ6TH9RW+5Hlu9L5727obwntcxHDSspmFShT8QpKXT3z1GA6C4uyHcgoaJo6OPBtxIx5Z8tfRUX9q6u9z2rGkU5XjI9rios5qGSRFHncDQoyG8u6Zh6SkTR3tXIkEjdNltrEbf2n1+XHL8QlYiQkVhDkespmFSxLF5pyLbNCw9ZRLU4aRVYvWeysuy9FQ8ockKPVG7JI8syrGahkkZ3fNORaSnstwucrNctHotaJg4OpweUfkxGsJzLGjE1Xta9EgVhTnWpmFSRve8UxHpKQhPWmjpKRNHXw24kcLpKWvT6FusqVgqinKs95RJGY1Oeqqkd9DI8XS3d2QSCxqDLNymkZtQQ7jVNPpytLOL4l7/CcMqCrNpaPMRCEad9NiYIRVeFKz3yO9MXVPDgsYgS6T3lLVpxNe733ukiqIcgkrGrsFs0ktdq5dst4vi3F49/XKzaLX0lImne86krNjTiIAFjViaO7p69EaJZKPCTSqpb/VRXph9XKcNq2mYhBzrPZVIl1tr0+hL7wngIvWef6rV67dUlRk29a1eyguPv8EJLcRkQcPE0dfsrJFsRHhswaDS1O6LmZ6CUE2jwxfgkh/8ne8+tWUoT9GYbkdafd2130jWe8okJBw0Ys095XIJ2R6XzT3Vh1afn6DSY1qGSOGgsftIK0+/c5AjrT4eWrOX3UfahvI0jQGcmkZBtKDhoc0XyLhacNz1NEz/dHQFyM1yxV0gKNfjwmvpqajC/d77qmnkZ3v4wIwKHlyzl1PK86kckUdju49v/WUTF0wp5/KZo5k+umgoT9mcpFSVI20+RkZJT4VHhbd6/X3+LaejuDUNEVkmIodFZFNEWZmIrBSRHc5zacR7t4lIlYhsF5GFEeXnisg7znv3OMu+4iwN+4hTvk5EJkXss9T5HTtEJLwkbEoLTYsePxbnZrm7e1qZnpqdEbZ9NYQDfOOK02nz+tl04Cg3XHAKN7//VFZXHeEHz23nvpd2DdWpmpNci9ePzx+Mmp4qjli9L5Mkkp76LbCoV9k3gFWqOg1Y5fyMiMwktFzrLGefe0UknKf5BXAzoXXDp0Uc8yagUVWnAncDdznHKiO0tOx5wDzg9sjglKrafYGYqamw8Op95nhNcWoaADPGFHHduRPIdrv4p3Mq+eKl01j11Ys5bUyRDfwzQ6a+NTxGI0pDeIbOPxU3aKjqS4TW7o60GFjuvF4OXBNR/gdV9arqbqAKmCciY4FiVV2jqgo82Guf8LEeA+Y7tZCFwEpVbVDVRmAlxwevlNPR1ff0F5FC64Rb0IimqSM8AVzsKv1/LZ7Fs19+H6OKcnG7hCkVhYwuzu0ebGVMstU7NyjRlnTN1EkLB9oQPlpVDwI4z6Oc8kpgf8R21U5ZpfO6d3mPfVTVDzQD5TGOldJiTX8RKTfLZdOI9KGpj7l8esvNcnNqRWGPsvLC7O67P2OSLVyrjd6mcfKmp/ojWuuvxigf6D49f6nIzSKyXkTW19XVJXSiyZJoeiony21To/ch3KbR1zQisZQXZFPf5iVUoTUmucJjhaJ3uT3WEJ5JBho0ap2UE87zYae8GpgQsd14oMYpHx+lvMc+IuIBSgilw/o61nFU9T5VnauqcysqKgb4kQZHvPXBw3Kz3HgtaETV3NFFbpare+R8f5QX5tDZFezu+mxMMoVrtWUFUWoazgJimTZp4UCDxgog3JtpKfBERPkSp0fUZEIN3q85KawWETnfaa+4sdc+4WNdC7zgtHs8BywQkVKnAXyBU5bS2n3+hHpP5Vl6qk9N7T5G5PXdcyqWcuc/r6WozFA40uplRH4WWe7jv0ozNT0V99tNRH4PXAKMFJFqQj2avgc8KiI3AfuA6wBUdbOIPApsAfzAraoavuW7hVBPrDzgGecB8ADwkIhUEaphLHGO1SAi3wFed7a7Q1V7N8innA5fIKGGcOs91bem9q64jeB9CfdiqW/zMrE8fzBPy5jj1Ld5u29UesvNcuFxScY1hMcNGqp6fR9vze9j+zuBO6OUrwdmRynvxAk6Ud5bBiyLd46ppL0rwfSUx8Zp9KW5o+8ZbuMJj8y1moYZCkdafVF7TkFoaeKiXA+tGRY0bBqRQdaR8DgNl3W57cOJBI1wbtm63ZqhcKTVS0UfQQMyc/4pCxqDKBBUvP5gQumpotwsWr1+ghk2L81gGIz01BFba8MMgfpWX9RG8LBMnOnWgsYgCnehTSQ9VVaQTVBDK9Q47SbtAAAgAElEQVSZnpo6fDGnEIklP9tDfrabBktPmSTzB4I0d3TFDBqZuKaGBY1B1N69lkb83lPhP7R6S6P00NkVoLMreEITvJUVZNt1NUnX6AxCjTaFSFhRblbG3Rha0BhE3Uu9JtCmUeoEjUb7cuvhaEf8eafiKS/MsfmnTNKF281KY9SKi3M9NrjP9C2RBZjCyvKtwTaapu4ZbgceNEYWZNt1NUkX/hvrq8sthCYttPSU6VP3AkwJBI3SgtCXYmO7fblFqj3aCRCzR0o8ZQU2/5RJvnDQKIuZngrVNDJpWhsLGoMo3LWuMCd+m0Z4PEFDW2blO09UdWMHAOPLBj4wr7wwh4Y2X0b9RzVD509vVPPp5a9T09QRc7sGp4deWYz0VFFuFoGgZtS0NhY0BlFNU+gueeyIvLjb5mW7yc1yWU2jlwONHbhdwuiigdc0RhZm4wsEacmwXLIZGn95q4bntx7mqp+tZvuhlj63C9/wlcbpPQWZNT26BY1BdKCpvV9feGX5lnvvrbqxnTHFuXiizOWTqO6xGi3WGG76r6q2hXmTy+jsCvCbV3b3uV1Dm5fiXE/UeafCwvNPtXozJ6NgQWMQHWjs6NcXXqk12B7nQFMH40vj19RiGVWUC8BhCxqmn1q9fmqaO7l4egWXzRzNc5sP0RWIPrFoQ3vsMRqQmTPdWtAYRAeaOqjsxxdemQWN4xxo7N81jGaUU9OzoGH6a+fhVgCmVBRyxeyxNLZ3sW5X9HlSG9q88YOGpadMLAcaOxifQHtGWFlBtrVpROgKBDl0tJPxpSc2O+2oYqem4fTEMiZRO5ygMW10IZfMqCA/283Tmw5G3TY0hUjsVHQmTo9uQWOQhL/w+nOXXGptGj0cau4kqPQr8EZTnOshx+Oymobpt6rDrWS5hVPK8snNcnPpaaP42+ZDUXviNbb7KCuIPZ6oe/U+q2mY3mqPhr7wKvtZ02jp9PeZMz3Z7G9sBzjhNg0RYVRxTveYD2MSVXW4hckjC7rbJc87tZwjrT5qmnv+LakqDW3xaxqFlp4yfTngjC/oV00jPJWIpaiAgV3DvowuyuXwUatpmMTUNHWw6UAzOw63Mm1UUXf5zLGh11trjvbYvsXrpyugcWsahdkeRCw9ZaI44AwEGtePmka5rf3QQ3VjByIwtuTEg8ao4hxqW6ymYeLbV9/O4v95hSt/tpq99e1MHVXY/d6MMcWIwJaDPYNGeM64eDUNl0sozPZY76kwEdkjIu+IyFsist4pKxORlSKyw3kujdj+NhGpEpHtIrIwovxc5zhVInKPs444zlrjjzjl60Rk0omcbzJ13yX3I2iU2vxTPRxo6mB0US7ZnhO/lxlVlEud1TRMHK1eP594YB1dgSCfff+pVBTlcNG0kd3vF+Z4mFRewNZeQaM+gXmnwjJtevTBqGl8QFXnqOpc5+dvAKtUdRqwyvkZEZlJaP3vWcAi4F4RCU/S9AvgZmCa81jklN8ENKrqVOBu4K5BON+kONDUwcjCbHITmOE2rKx7ptvMqbqeiL31bSfcnhE2qjiHFq+/e7p6Y6JZvaOOfQ3t/OSjc7jtg6fz+n9cxnsmlfXY5vSxRcfVNMLrtcQaDR4WWnAtc/6PJyM9tRhY7rxeDlwTUf4HVfWq6m6gCpgnImOBYlVdo6EuCg/22id8rMeA+eFaSKo50NTRr1oGHJu0sMHaNFBVth1qYcaYovgbJ2B0eICf1TZMDOt2N5Cb5eK9U0b2uc3MscXsrW/v0S4R/j+bSE0j02a6PdGgocDfRGSDiNzslI1W1YMAzvMop7wS2B+xb7VTVum87l3eYx9V9QPNQHnvkxCRm0VkvYisr6urO8GP1H+qypaao0ypKIy/cYTS/GzcLqG22XLvB5o6aOn0c/rY4kE53qjiUK7ZelCZWF7b3cA5E0tjpkTDf5OR81DVNIXa3yoSmDKoojAn7uSH6eREg8aFqnoOcAVwq4i8P8a20WoIGqM81j49C1TvU9W5qjq3oqIi3jkPup11rdS3+Tjv1LL4G0fIcrs4pTyfHYf7nhTtZLHtYOganD52cGoaNpWIiedoZxdbDx49Lh3V28xxoaCxOaIH1c66UCo1kXT0mRNK2FPfnjELrp1Q0FDVGuf5MPBnYB5Q66SccJ4PO5tXAxMidh8P1Djl46OU99hHRDxACRB9TP8wWutMM3De5OMqQXFNH1XUPQr1ZLbtUOg/5Iwxg1PTGJ1mNY0X363j3/64kc8/vIFddfb3MBQ27G0kqHDe5NhBY0xxLmOKc3l9z7GvnqrDrQlnFs6eEOoL9FZ108BPNoUMOGiISIGIFIVfAwuATcAKYKmz2VLgCef1CmCJ0yNqMqEG79ecFFaLiJzvtFfc2Guf8LGuBV7QFFwkYd3uBkYX53BKef+nv5g+upC99e14/Zkz3/5AbD3YwsSy/ITWIklESV4W2R4XdWlQ09jf0M4nf/Maz2+t5eUdR7jipy/z0rtDn2bNZF5/gKffOdijY8RruxvwuISzJ5bG2DM0WPS9U8pZs7OeYFAJBpVdda1MTTBonDm+BJfAm/tO8qABjAZWi8hG4DXgKVV9FvgecLmI7AAud35GVTcDjwJbgGeBW1U1/E15C/BrQo3jO4FnnPIHgHIRqQK+gtMTK1X4/EG6AkHW7arnvMnlDKSNfuroIgJBZVddWxLOMH1sPXSU0wapERycUeFFORxKg5rGio01qMJfv3ARz3/lYgpzPPz5zQPDfVoZ5SfP7+DzD7/BxT/4B49tqGZ/QzuPvr6fc04pTWilzQumlFPf5uPdwy0caOrA6w8yZVRiQaMgx8OMMcW8ua9xQOeuqnzlkbd4bEN1/I2HwIBv61R1F3BWlPJ6YH4f+9wJ3BmlfD0wO0p5J3DdQM8x2a771Rr2N7TTMID2jLDpo0N/eO/WtgxaI3C66fAF2HOkjavOHDeoxz21opAdtamf6lnxVg1zTyllgrNa4ezKkpiL/5j+2XOkjQde3s3F0ys42tnFv/1xI/nZbjwu4b//6YyEjvHeqaHeVa9W1TN5ZAFAj0GA8Zw9cQR/3VhDMKi4XP27uVy7q4E/vXmA2pZOrj13fPwdksxGhA/QweYONu5vIj/bTX62m/dPG1gD/OSRBbhdkhZfbsmy7dBRgjp4jeBhp48poupwa0rP7bXt0FG217aw+OzK7rIZY4qoqmvFn8LnnU7++5mteNzC9689k8c/916+deVMKopy+OUN5ybcLlE5Io9J5fm8uvMIO502p0TTUwBnTxhBS6efXUf6n1G4/+VdQKghPhWy8xY0BuiVqnoA7r9xLlvuWNR9l9hfOR43k8rzebf25LmzXLGxhld3Hun++bnNtXhcwrwBdCSI5bSxRfgCQXYP4D/qUPnzGwfwuIQPnTG2u2z66CJ8/iB7G9qH8cwyw/ZDLTy3uZbPvO9URhfn4nIJN100mRe/9oGYYzOiuWDKSNbuauD1PQ2UFWQnNLAvLJS+hj+u3x9/4whVh1t4YdthJpTl0dTeddzEicPBgsYAvVp1hPKCbGaMPvG74+mjT54eVC2dXXztjxv57EMbqGnqIBhU/rqxhoumjYy7oE1/hdN9vaeASBWdXQEeXb+fy2eO7vHZw39T71qK6oT96sWd5Ge7+eR7J53wsT5x/kR8gSDPba7tVy0DYGJ5Ph8+u5LfvLqnX2M2/vTGAdwu4fYrZwGw+UBzv35vMljQGABVZXXVEd47dWS/85PRzBhTxN76NprbM2eqgb48t7kWrz9Ihy/A1x7byNrd9Rxo6mDxnMFtzwA4dWQhWW5hW4p++T79zkEa27u44fxTepRPHVWICGwf4trnlpqj3Pq7N6hvTf0eZ4k40NTBio01LHnPxH7VCvoya1wJd1wd+vKeMqqg3/t/5fLpoPCT599NeJ9/bK9j7imlvHdqqKayuWb4b4AsaAxA1eFWDrd4uWjq4KRTLpkxiqDCyq21g3K8VPbEWweYUJbHHYtn80pVPTc+8Bq5WS4unzlm0H9XtsfFlIpCtqVoTeOhtXs5taKAC6b0/DvKy3ZzStnQpixVlf/8yzs89fZBvv742ymRO0/UW/ubuPOpLd1jfcL+tKGagCo3vW/yoP2uJfMm8v1rz+Smi07t977jS0O1jafePphQe1Xt0U62HDzKJTNGkZ/tYUpFoQWNdBX+cr9ogI3fvZ01voRxJbk828eykpnicEsnr1QdYfFZlXzsvIks/9Q8ZlWW8PHzThm08Rm9nT62OCVqGi2dXRyJuIPfdKCZN/c18YnzTonaVXv66KIT7kF1pNXLlgS/ZJ7ddIg39jVxwanlPL/1MA+v23dCvzvZNu5v4scr32XV1lpueGAd97+8m0U/eZkH1+zp3uZvW2o5e8KIfs8JF88/z53Qr55TkS6cNpI2X+C4CRAjbdzfxD2rdrByS+h75pIZoe+ZWeOK2VJj6am0o6o8tr6aeZPKBu2PUURYOHsML+04Qqs3cyY2ixQIKrc/sZmgwjVOT6GLp1fwxK0X8q0rZybt9542poiDzZ3DPoXDF373Jlf9bDWdXaGhSQ+v20tulouP9NGFcsaYIvb0miSvP9q8fj76qzVc/fPVvPhuHapKIBi99tAVCHLXs9uYPrqQB2+ax4VTy/n+s9tSdsr+uhYvn35wPfes2sFNy9dTmOPhqX+9iHNPKeX+l3cRDCoHmjp450AzC2YNfg32RIRHn7+2O/rEFqt3HGHJfWv58cp3+a+/bmZMcW73+KVZ44qpae7k5gfX9+hIMtQsaPTTG/sa2XWkjWvnDm5/6Q+eMRafP8iqDExRhVIfm3hm0yH+80OnD/gubSDmTgqN9v3r2zVxtkyeHbUtvPhuHQebO1n+6h6aO7r4y5s1XDOnkpK86Cu/XT5zNIGgct9Lu6K+39Tu6w4owaA6vWxq6QoEUVW++ed32H2kjcrSPD730AbO/+9VzPy/z/LxX6/lL28e6BFAfrduH3vq2/nGFaeR5XZx+1WzaPMF+pV7Hyqqylf/uJHmji4eumke37pyJr//zPnMGlfCDeefwv6GDtbtbmDl5kMALJg5epjPuKfRxblMKs/vnnoo0s66Vj794OucUp7P1xedRldAufT0Ud010avOGsfVZ43jrf1NfO6hDdQe7WTrwaP8+uVd/PC57bQN0Q1ncnICGUhV2XTgKL9+eTf52e4eXSQHw7kTS6kckccDq3dz9VnjBjS6PFU9sHo3v39tH5+/ZAqffl//c8En4pyJpbxnUin3/n0n/zx3Qr/WOxksv311D9keF2dWlnDvP3aysbqJjq4An+jVAB7pzPEjuOqscdz/8i4+dt7EHqsZtnr9fOie1RTlevjrFy/isw9t4IVtoSnebrlkCtNHF/LEWzV89fLpfPQ9E7jtT++Qn+OhvCCbl96t48uPvMUvX9zJgzfNIy/LzU9X7eD8U8v4wIzQhNTTRxfx8fMm8vC6fVQU5vCZ9586LNctmn+8W8dL79bx7atm8r5pFbwvIkW8cNYYinI8LHtlNwebO5g6qpBT+9nLaSjMm1zGc5trOdjcQZvXz9RRRfgDQb766EZys9w8+Kl5jCrO5X3TRvaYmmhsSR73XH82e460seinL3H9/WvZW9/efQMwpiQ35t/UYLGaRgJUla899jZX/Xw1z2w6xFVnjqNgkHPwLpfw5cum8XZ1M0+/c2hQjz2cnnr7IP/v6a0smjWGf1swY8h/v4jw5cumc+hoJ8te2T3kDbzNHV386Y0DLD5rHN++ehZtXj/Pbz3MkvdMYHZlScx9/33hDIJB+NHfet7x//C57Rxo6mDboRb+5Tev88K2w3z+kilcM2ccv3pxJ9/6y2bmnlLK5z8wlVHFuTzwyffws+vP5ttXz+L5r1zMzz92Nnvr27nxgdf46K/W0tDm45sfPL3Hjcq/LZzBgpmj+dHKd/nwva9yOElL575T3cxHf7Um4W7Rv/zHTsaW5PKx847/cszLdnPlWeNYuaWWrQdb+PRFg9cAPpjOm1xOc0cXF3//Hyz++Ss0tPl4YPVu3trfxHcWz2ZUcWiG5tmVJRTlHl8TnTSygK8tPI1ddW1cdeZYXvuP+UypKODJIapNW00jAf/z9yoe21DNZ943mSvOGMvMJE338U/njOf+l3fxw79tZ/7po1Lm7m6glq3ezXee2sLZE0bwo38+a1C6Jw/Ee6eU8/7pFXz/2e28UnWE//nYOYzIH9wxIX1ZsbGGjq4AN14widmVJaz75nyK87LIcse/X5tQls8nL5zE/S/v4lMXTmbmuGJe39PA8jV7uPGCU9h9pI2Xdxxh3qQy/m3BDNp8fl7f00hDm48fXncW7ijX2+USrjxzHEW5WXx6+euU5mdz3w3ncub4ET22K87N4hefOJcXttXyhd+9ybW/WMNjn7ug+wttMHQFgnztsY1sO9TCjcte44+fvYBJI/vuyvrmvkbW7W7gPz90ep/rX3xp/jQmlOVx9VnjGF86sAG3yXbh1JHkZbk5e+II1uyq565ntvHk2zVcdvoorjorsa7nn7pwEgtmju4eVHzlmeO454Ud1B7tZPQg/htFI+nUtS4Rc+fO1fXr1w/a8TbsbeDaX65h8VnjuPujc5KeNnp5Rx03PPAaHzlnPD+87sy0TVNtP9TCop++xGWnj+Zn15897AGwKxDk4bV7+e5TW7n6rHH8+KNzhuT3Lv75arz+IM986X0D+rdsbu/i/T/4O7PGFfOda2Zz/X1ryct28+QXL6KuxcsdT27hW1fO7J4O40BTB83tXd1rQMSy+0gbIwuzo97NRnpzXyPX37+WC04tZ9kn33PCf5Oqyut7Gnny7RoeXLOXry86jfte2smYkjz++oUL8UQJqM+8c5BvPbEZfzDI6q9fmrTedkPF6w+Q7Xbxhd+/yVNvHyTLLfzt/1zcPa9Vf1UdbuGyH7/E7VfN5F8uHFgNS0Q2RCzb3SdLT8XQ2RXg3x97m3EleXz3w2cMyRf4+6ZV8K/zp/H4G9U82s8pB1LJ3SvfpTDbw/c/cuawBwwILXj1yQsn8/lLpvCnNw/w9+2H4+90gt6tbWFjdTPXzZ0w4L+dkvwsvnL5dF7dWc/8H73I0c4ufvmJcynKzeLUikJ++y/zesyfVDkiL6GAAaF5z+IFDICzJ5by9UWn8fftdfzHXzZx30s7T6g32s9fqOKff7WGB9fs5aqzxnHLJVP47jVnsPXgUe5/eTefeXA98+58no//ei0b9jbw0Nq93PLwG4wpyeHhT5+X9gEDQtMHiQj/euk0PC7hUxdNHnDAAJg6qojTxhTx5NvJ77ZvNY0oGtp8XH/f2u4Rucs/NY+Lpw/dioDBoHLdr9ZwsKmDF//9AwmlMoabPxDEFwiSn+3hjX2N/NO9r/Lly6bx5cumD/ep9eD1B/jQPas51NzJD649kysGuUMDhK7Fk28f5PE3qlmzs56135zPyML4y4LGsmFvAy9ur+PcSWVD+rcYFgwqNz+0geed3n2TRxbw5cumsb+hnavPqmRinLVk1uys53/X7eUj51Tyuf99g0tnjOI/PnQ640vzEBFUlevvX8vaXaE1Lq44Yywb9jRwuMVLQJVLZ4zilzecmxb/F/rrYHMHY4pzT/imdNXWWtwu4ZIZo+JvHEWiNQ0LGlHc+rs3+NvmQ3zu4inMmTCC+acPfbe9VVtruWn5en66ZA6L51TG32EIrd/TwIqNNbxnUhkXTh3JqzuP8MPntlPT1MmVZ43luU2HKMjx8PxXL6Y4gTvZoVbT1MEtD7/Bxv1N/NfVs1iawLxEqspDa/eydlc9LZ1+Th1ZwMfOO4UZUdYA+faKzd09pj514WS+ccVpSfgUQ09V8fqDbDrQzM0Pbegex1GSl8W/L5pBcW4W75lUxpiSUE69ub2LN/Y1Upzn4ZO/eZ2WzlCX0MIcD6u+evFxufcdtS1840/v8KX503j/9AqaO7r4z79s4mhHqHaVyLoXZuAsaPRTc0donv0OX4DVVUf42sIZ3PqBqUk4w8QEg8pld79IXpabFV+4KGqj5lCra/Hy9cff5oVth3G7pEdf/6mjCjmjsoQn3jrAORNL+fnHzun+8khFXn+AL/7uTf62pZbv/dMZLJk3Meb2j2+o5qt/3MjEsnxK8rKoOtxKQY6HJ794EQea2qk63IrPH2R/Ywf3vbSLT753Et+6cmZK/LslQ32rl70N7RTnerj14Te7a+UugctOH81nLz6Vrz/+DlXORJwjC3P49dK5/H7dPi6aNjLhBl8zdDIqaIjIIuCngBv4tap+r69tBxw02rv46H1rcLuEWeOK+X8fPiNqg9xQ+uP6/XztsbeZN6mMr19xGmdUlnT3GlHVIWskDwaVx9+o5q5nt3G0089XL5/Ox88/hTU769lV18rsyhLmTS4jyx1aXrWsIDstviy9/gCfeXADL71bx+cunsKpFQWgsGDWaEbkZ9Pq9fP8llqOtHq5e+W7zKos4fefOd9Z/6SFa/7nFdwu4Whnz0FVF0+v4NdL52ZkKiUanz/I3vq2UIP/poMsf3UvrV4/eVluvnvNbGpbOrlk+qiE21rM8MiYoCEibuBdQkvHVgOvA9er6pZo2w9276nhpKo8tqGaO/66hRavn4JsN5+7eAqtXj/L1+whP9vDqKIcRhXnMqooh9HFOYwpyeOS6RVUjshjX0M7h452sq++nc01zcyqLOH8yeWs210fms5ChJqmDioKczhrQglul4tAMIjH5WL66CJ21rXy6Pr9vLDtMNWNHZw1YQTf/8iZUVMy6crnD/J/n9jEH14/1unA4xJGF+fS2O6j3Rea9mNUUQ5/vvXCHlPH/G3zIX7w3HaWzJvIwlmjyfG4yfa4KM71pG2vt8FQe7STX724i0WzxzBv8sBWtDRDL5OCxgXAt1V1ofPzbQCq+t/Rts+koBHW0OZj7a56nnjrAM9tDjVEXn3WOIpyPdQe9VLX0hl6bvV2p4yKcjy0REwrkO1x4fMfP7Nm7zRTmEsgqJCb5eKiqRUsnjOOD50xdtjGWiSTqrKxupmiXA8dvgDPbDrIoWYv+dlurjm7kqkVheRlu/scG2BMJkg0aKRD37VKILLvaTVwXuQGInIzcDPAxImxc9PpqKwgmw+eMZYPnjGWt/Y3ke12Ra3qB4JKdWM7f91YQ01zJ2dWljChLJ+xJblMKi9g7e56th9q4YIp5ZQX5BBUZWRhDoeOdrL90FEEweUSOrsCbDrQTHlBNh8+Z3yf8yNlChFhzoRjg9vijdQ25mSWDjWN64CFqvpp5+cbgHmq+sVo22diTcMYY5Itkwb3VQMTIn4eDwzflKXGGHMSS4eg8TowTUQmi0g2sARYMcznZIwxJ6WUb9NQVb+IfAF4jlCX22WqunmYT8sYY05KKR80AFT1aeDp4T4PY4w52aVDesoYY0yKsKBhjDEmYRY0jDHGJMyChjHGmISl/OC+/hKRFmB7gpuXAM2D+OtHAkcG8XiDfX6pfDy7difGrt+JGczrl+qfNdrxRgIFqhp/sRZVzagHsL4f2943XL97mM4vZY9n186uX6ZcvzT4rMcdrz+f/2RPT/11uE8gjsE+v1Q/3mBK9c+aytcOUv/zpvL1S/XPekLHy8T01HpNYP6UTPvd6c6u3Ymx63diTvbr15/Pn4k1jftO0t+d7uzanRi7fifmZL9+CX/+jKtpGGOMSZ5MrGkYY4xJEgsaMYjIBBH5u4hsFZHNIvIlp7xMRFaKyA7nudQpL3e2bxWRn/c61kdF5G3nON8fjs8zlAZw7S4XkQ0i8o7zfGnEsc51yqtE5B45CdZSHeTrd6eI7BeR1uH6PENtsK6fiOSLyFMiss05zveG83OlhMHsypVpD2AscI7zuojQWuUzge8D33DKvwHc5bwuAC4CPgf8POI45cA+oML5eTkwf7g/X4pdu7OBcc7r2cCBiGO9BlwACPAMcMVwf740u37nO8drHe7PlW7XD8gHPuC8zgZePhn+/mJe2+E+gXR6AE8AlxMaPDjWKRsLbO+13Sd7BY33AM9H/HwDcO9wf55UvHZOuQD1QI6zzbaI964HfjXcnyddrl+v8pMmaCTj+jnv/RT4zHB/nuF8WHoqQSIyidDdyDpgtKoeBHCeR8XZvQo4TUQmiYgHuIaeqxFmtAFcu48Ab6qql9Aa8dUR71U7ZSeNE7x+J73Bun4iMgK4CliVzPNNdWmxnsZwE5FC4HHgy6p6tL8pdVVtFJFbgEeAIPAqcOqgn2gK6u+1E5FZwF3AgnBRlM1Omi5/g3D9TmqDdf2cm73fA/eo6q4knW5asJpGHCKSReiP7mFV/ZNTXCsiY533xwKH4x1HVf+qquep6gWEqsg7knXOqaK/105ExgN/Bm5U1Z1OcTWhdeHDTpo14gfp+p20Bvn63QfsUNWfJP/MU5sFjRicXjoPAFtV9ccRb60AljqvlxLKl8Y71ijnuRT4PPDrwT3b1NLfa+dU/Z8CblPVV8IbOymEFhE53znmjSRwvdPdYF2/k9VgXj8R+S6hSf6+nOzzTgvD3aiSyg9CPaEUeBt4y3l8kFBvqFWEagurgLKIffYADUArobvkmU7574EtzmPJcH+2VLt2wH8CbRHbvgWMct6bC2wCdgI/xxmUmsmPQb5+33f+FoPO87eH+/Oly/UjVLNVYGtE+aeH+/MN58NGhBtjjEmYpaeMMcYkzIKGMcaYhFnQMMYYkzALGsYYYxJmQcMYY0zCLGgYM8RE5HMicmM/tp8kIpuSeU7GJMqmETFmCImIR1V/OdznYcxAWdAwpp+cCfCeJTQB3tmEpt2+ETgd+DFQCBwBPqmqB0XkH4TmG7sQWCEiRYRmnP2hiMwBfkloCu6dwKc0NFfZucAyoB1YPXSfzpjYLD1lzMDMAO5T1TOBo8CtwM+Aa1U1/IV/Z8T2I1T1YlX9Ua/jPAh83TnOO8DtTvlvgH/V0FxlxqQMq9hLw0cAAAEBSURBVGkYMzD79dgcRf8LfJPQ4j0rnZlU3cDBiO0f6X0AESkhFExedIqWA3+MUv4QcMXgfwRj+s+ChjED03v+nRZgc4yaQVs/ji1Rjm9MSrD0lDEDM1FEwgHiemAtUBEuE5EsZ22GPqlqM9AoIu9zim4AXlTVJqBZRC5yyj8++KdvzMBYTcOYgdkKLBWRXxGaMfVnwHPAPU56yQP8BNgc5zhLgV+KSD6wC/gXp/xfgGUi0u4c15iUYLPcGtNPTu+pJ1V19jCfijFDztJTxhhjEmY1DWOMMQmzmoYxxpiEWdAwxhiTMAsaxhhjEmZBwxhjTMIsaBhjjEmYBQ1jjDEJ+/+8SSgnFKcDNAAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -252,10 +2305,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 11, + "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", @@ -274,7 +2325,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +2349,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 13, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAD8CAYAAABQFVIjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAG3xJREFUeJzt3X+QldWd5/H3BxvBMaiAYPghwlYcN5rM6tCFTjm7K7oC2aTEzGiG1VWqxiqMmipnnVqBlS13lanVqa1xlrJitEYT1CgyTijZUUIQteLMEqAZNYiGabI6yEhs3EbBVNGx5bt/3NN6uTbdt7ufe+/z3Pt5Vd26T597znkOh+77vefH81xFBGZmZiM1qtENMDOz5uCAYmZmmXBAMTOzTDigmJlZJhxQzMwsEw4oZmaWCQcUMzPLhAOKmZllwgHFzMwy0dboBtTT6aefHjNnzmx0M8zMCmXHjh3vR8SkwfK1VECZOXMmHR0djW6GmVmhSPqnavJ5ysvMzDLhgGJmZplwQDEzs0w4oJiZWSYcUMzMLBMOKDnSdegI33pwC12HjzS6KWZmQ+aAkiOrNney/e1uVj3f2eimmJkNWUtdh5JX56zYQE/v0U9/fnzrXh7fupcxbaPYvfJrDWyZmVn1PELJgZdvn8sV509l7OjSf8fY0aNYeP5UXl46t8EtMzOrngNKDkw+ZSzjxrTR03uUMW2j6Ok9yrgxbUweN7bRTTMzq5qnvHLi/Y96uPbCs7hmzgye2LaXA16YN7OCUUQ0ug11097eHr6Xl5nZ0EjaERHtg+XzlJeZmWXCAcXMzDLhgGJmZplwQCkYX01vZnnlgFIwvprezPLK24YLwlfTm1neeYRSEL6a3szyzgGlIHw1vZnlnae8CsRX05tZnvlK+Sp0HTrCd558hfuvucAjAjNrOb5SPkPeWWVmNjhPeQ3AO6vMzKrnEcoAvLPKzKx6VQUUSW9L2inpVUkdKW2CpE2SOtPz+LL8yyXtkbRb0vyy9Nmpnj2SVklSSh8j6amUvlXSzLIyi9M5OiUtLkuflfJ2prInjrw7juWdVWZm1RvKCGVuRJxftjCzDNgcEWcDm9PPSDoXWAScBywAvivphFTmAWAJcHZ6LEjpNwAHI+JLwH3AvamuCcCdwIXAHODOssB1L3BfOv/BVEfm+nZWrbv5Yq698CwOfNRTi9OYmRVeVbu8JL0NtEfE+2Vpu4FLImK/pCnASxFxjqTlABHxP1K+jcB/A94GXoyIf5nS/0Mqf2NfnojYIqkN+BUwiVJguiQibkxlHgReAtYAB4AvRkSvpN9L5T8dDfXH34diZjZ0We/yCuAnknZIWpLSzoiI/QDpeXJKnwa8U1Z2X0qblo4r048pExG9wIfAxAHqmgh8kPJW1mVmZg1Q7S6viyPiXUmTgU2SfjFAXvWTFgOkD6fMQHUd25hSAFwCMGPGjP6ymJlZBqoaoUTEu+m5C1hHaT3jvTTVRXruStn3AWeWFZ8OvJvSp/eTfkyZNOV1KtA9QF3vA6elvJV1Vbb9oYhoj4j2SZMmVfPPNTOzYRg0oEg6WdK4vmNgHvA6sB7o23W1GHgmHa8HFqWdW7MoLb5vS9NihyVdlHZ3XV9Rpq+uq4AXorS4sxGYJ2l8WoyfB2xMr72Y8lae38zMGqCaKa8zgHVph28b8ERE/FjSdmCtpBuAvcDVABGxS9Ja4A2gF7glIj5Jdd0E/AA4CdiQHgAPA49J2kNpZLIo1dUt6W5ge8p3V0R0p+OlwBpJK4FXUh1mZtYgvpeXmZkNyPfyMjOzunJAMTOzTDigmJlZJhxQzMwsEw4oZmaWCQcUMzPLhAOKmZllwgHFzMwy4YBiZmaZcEAxM7NMOKCYmVkmHFDMzCwTDihmZpYJBxQzM8uEA4qZmWXCAcXMzDLhgGJmZplwQDEzs0w4oJiZWSYcUMzMLBMOKGZmlgkHFDMzy4QDipmZZcIBxczMMuGAYmZmmXBAMTOzTDigmJlZJhxQzKxpdB06wrce3ELX4SONbkpLckAxs6axanMn29/uZtXznY1uSktqa3QDzMxG6pwVG+jpPfrpz49v3cvjW/cypm0Uu1d+rYEtay0eoZhZzdV6Kurl2+dyxflTGTu69JY2dvQoFp4/lZeXzq3J+ax/DihmVnO1noqafMpYxo1po6f3KGPaRtHTe5RxY9qYPG5sTc5n/fOUl5nVTD2not7/qIdrLzyLa+bM4Ilteznghfm6U0Q0ug11097eHh0dHY1uhlnL6Dp0hJXPvclPdv2KIx8fZezoUcw/74vc8fUve/RQIJJ2RET7YPk85WVmNeOpqNZSdUCRdIKkVyT9bfp5gqRNkjrT8/iyvMsl7ZG0W9L8svTZknam11ZJUkofI+mplL5V0syyMovTOTolLS5Ln5XydqayJ46sK8ysFvqmotbdfDHXXngWBz7qaXSTrEaqnvKSdBvQDpwSEd+Q9OdAd0TcI2kZMD4ilko6F3gSmANMBZ4HfjsiPpG0DbgV+BnwHLAqIjZIuhn4nYj4tqRFwDcj4o8kTQA60nkD2AHMjoiDktYCP4qINZK+B7wWEQ8M9G/wlJeZ2dBlOuUlaTrwdeCvypIXAqvT8WrgyrL0NRHRExFvAXuAOZKmUApGW6IUxR6tKNNX19PAZWn0Mh/YFBHdEXEQ2AQsSK9dmvJWnt/MzBqg2imvvwRuB46WpZ0REfsB0vPklD4NeKcs376UNi0dV6YfUyYieoEPgYkD1DUR+CDlrazrGJKWSOqQ1HHgwIEq/7lmZjZUgwYUSd8AuiJiR5V1qp+0GCB9OGUGquvYxIiHIqI9ItonTZrUXxZrAb7Hk1ntVTNCuRi4QtLbwBrgUkmPA++laSzSc1fKvw84s6z8dODdlD69n/RjykhqA04Fugeo633gtJS3si6zz/E9nsxqb9CAEhHLI2J6RMwEFgEvRMR/BNYDfbuuFgPPpOP1wKK0c2sWcDawLU2LHZZ0UVoDub6iTF9dV6VzBLARmCdpfNpFNg/YmF57MeWtPL/Zp85ZsYGZy57l8a17iShdWDdz2bOcs2JDo5tm1nRGch3KPcDlkjqBy9PPRMQuYC3wBvBj4JaI+CSVuYnSwv4e4JdA31/1w8BESXuA24Blqa5u4G5ge3rcldIAlgK3pTITUx1mx2jWezx5Cs/yaEi3XomIl4CX0vH/Ay47Tr4/A/6sn/QO4Cv9pB8Brj5OXY8Aj/ST/n8pbU02O66hXFjXdegI33nyFe6/5oLcX3hXPoW38ptfbXRzzADfy8taQLX3eCrCm7Rv02555nt5WcurfJPuk8c3ad8byxrB9/Iyq1KR1ll8byzLM095Wcsr2pu0b9NueeWAYkax3qQfvO6zmYeVV35uj4tZw3gNxczMBuQ1FDOzYfA1PsPngGJmVsa36Rk+r6FY7hXpgkMrLl/jM3IeoVju+ROj1UORto/nlUcollut/InRo7L6K9r28TzyCMUaaqAF0Fb+xOhRWWP0bR9fd/PFXHvhWRz4qKfRTSoUj1Bs2LL4FD3Q/bNa8RNjK4/K8sDX+IyMA4oN20hupljtG2eRLjjMwsu3zz3uvbrM8s4BxYYsi0/R1b5xttonxlYclVnz8BqKDVkWaxt+4zw+z+NbUXmEYkOWVTBotemsahVtVOYdadbHAcWGJYtgULQ3TutfEb6YzOrDN4c0s2Ep0heT2cj45pBmVlOtfJ2Q9c8BxcyGxRsrrJLXUMxs2Lyxwsp5DcXMzAbkNRQzM6srBxQzM8uEA4qZmWXCAcXMzDLhgGJmZplwQDEzs0w4oJhlaKBvoDRrdg4oZhnyV/daK/OV8mYZ8Ff3mnmEYpYJ3yjRzAHFmkAe1i18o0SzKgKKpLGStkl6TdIuSf89pU+QtElSZ3oeX1ZmuaQ9knZLml+WPlvSzvTaKklK6WMkPZXSt0qaWVZmcTpHp6TFZemzUt7OVPbEbLrEiiYv6xb+6l5rdYPeHDK96Z8cER9JGg38HXAr8AdAd0TcI2kZMD4ilko6F3gSmANMBZ4HfjsiPpG0LZX9GfAcsCoiNki6GfidiPi2pEXANyPijyRNADqAdiCAHcDsiDgoaS3wo4hYI+l7wGsR8cBA/xbfHLK5+AuezOojs5tDRslH6cfR6RHAQmB1Sl8NXJmOFwJrIqInIt4C9gBzJE0BTomILVGKYo9WlOmr62ngshTI5gObIqI7Ig4Cm4AF6bVLU97K81uL8LpF7eVhOtGKo6o1FEknSHoV6KL0Br8VOCMi9gOk58kp+zTgnbLi+1LatHRcmX5MmYjoBT4EJg5Q10Tgg5S3si5rEV63qL28TCdaMVS1bTgiPgHOl3QasE7SVwbIrv6qGCB9OGUGquvYxkhLgCUAM2bM6C+LFZi/4Kk2vA3ahmNI16FExAeSXgIWAO9JmhIR+9N0VlfKtg84s6zYdODdlD69n/TyMvsktQGnAt0p/ZKKMi8B7wOnSWpLo5Tyuirb/BDwEJTWUIby77X8e/C6z6Z1V1450OccG4qXb5/Lyufe5Ce7fsWRj48ydvQo5p/3Re74+pcb3TTLsWp2eU1KIxMknQT8O+AXwHqgb9fVYuCZdLweWJR2bs0Czga2pWmxw5IuSmsg11eU6avrKuCFtM6yEZgnaXzaRTYP2JheezHlrTy/mY2QpxNtOKpZQ5kCvCjp58B2SmsofwvcA1wuqRO4PP1MROwC1gJvAD8GbklTZgA3AX9FaaH+l8CGlP4wMFHSHuA2YFmqqxu4O513O3BXSgNYCtyWykxMdTSMFy+t2dRzG3TR/n6K1t568XfKZ2TFup38cNterp0zg5Xf/GpNzmHWrIr291O09o5UtduGHVBGyNdCmA1f0f5+itberGR2HYoNzNdCmA1f0f5+itbeenNAGaFqFy8952p51Ojfy6It/hetvfXmgJKBahYvfYGY5VEefi+Ldg+0orW3nryGUmOtOudq+ebfSxsKr6HkRFHnXBs9FWK1VdTfy2bSjH9jDig1VtQ51zxMhVjtFPX3spk049+YvwK4Dop0vynfw6l1FOn3spk089+Y11DsGF2Hjhz3Hk7+9Go2ckX8G/Maig2Lp0LMaquZ/8YcUFrQYIuB3hZpVlvN+jfmKa8W1Gr3ITKzkal2ysuL8i2kmRcD66Hr0BG+8+Qr3H/NBU0xPWGWNU95tRBfezAyzbjN0yxLHqG0kGZeDKwlj+zMquMRSotp1sXAWvLIzqw6HqG0GH8H+9B5ZGdWHQcUsyr4qnKzwXnbsJmZDchXypuZWV05oJiZWSYcUMzMLBMOKGZNqhm/wKlIWrH/HVDMmpSv7G+sVux/7/IyazL+vvjGyrL/83L/OO/yMsupWk+F+Mr+xsqy/4s2yvGFjWZ1Vv4mUYuvD/CV/Y2VRf8X9f5xDihmdVLPNwlf2d9YI+3/l2+fe9yvCc4zr6GY1UkRv0s8T/KynlAvd6zbyRPb9nLiCaP4zSdHG/qFeF5DMcsZT0WNTNHWE0aqiHcG9wjFrI5ufKyDSePGHjMVUn4HaPs871prvGpHKA4oZpZrnipsPE95mVlT8FRhcTigmFnuFXE9oV7ydIsXT3mZmRXYinU7+eG2vTXdBZbZlJekMyW9KOlNSbsk3ZrSJ0jaJKkzPY8vK7Nc0h5JuyXNL0ufLWlnem2VJKX0MZKeSulbJc0sK7M4naNT0uKy9Fkpb2cqe2K1nWNmVnTnrNjAzGXP8vjWvUSUrmuauexZzlmxoWFtqmbKqxf404j4MnARcIukc4FlwOaIOBvYnH4mvbYIOA9YAHxX0gmprgeAJcDZ6bEgpd8AHIyILwH3AfemuiYAdwIXAnOAO8sC173Afen8B1MdZmYtIY+32Bk0oETE/oj4h3R8GHgTmAYsBFanbKuBK9PxQmBNRPRExFvAHmCOpCnAKRGxJUrzbI9WlOmr62ngsjR6mQ9siojuiDgIbAIWpNcuTXkrz29m1vTyuFlhSLdeSVNRFwBbgTMiYj+Ugo6kySnbNOBnZcX2pbSP03Flel+Zd1JdvZI+BCaWp1eUmQh8EBG9/dRlZtYS8naLnaoDiqQvAH8D/ElEHErLH/1m7SctBkgfTpmB6jq2MdISStNszJgxo78sZmaFVH5R7Morv9LAlpRUtW1Y0mhKweSHEfGjlPxemsYiPXel9H3AmWXFpwPvpvTp/aQfU0ZSG3Aq0D1AXe8Dp6W8lXUdIyIeioj2iGifNGlSNf9cMzMbhmp2eQl4GHgzIv6i7KX1QN+uq8XAM2Xpi9LOrVmUFt+3pemxw5IuSnVeX1Gmr66rgBfSOstGYJ6k8Wkxfh6wMb32YspbeX4zM2uAaqa8LgauA3ZKejWl/RfgHmCtpBuAvcDVABGxS9Ja4A1KO8RuiYhPUrmbgB8AJwEb0gNKAesxSXsojUwWpbq6Jd0NbE/57oqI7nS8FFgjaSXwSqrDzMwaxBc2Ws202u3GzZqV7+VlDddqtxs3a3X+xkbLXFG/vtTMRsYjFMtcHq/gNbPac0CxzOXxCl4zqz1PeVlN5O0KXjOrPe/yMjOzAXmXVwvL0xfumFnrcEBpQt6ua2aN4DWUJuLtumbWSB6hNBFv1zUbnKeEa8cBpYl4u67Z4DwlXDue8moy3q5r1j9PCdeetw2bWUvoOnSElc+9yU92/YojHx9l7OhRzD/vi9zx9S97FD8Ibxs2MyvjKeHa85SXmbUMTwnXlqe8zMxsQJ7yMjOzunJAMTOzTDigmJlZJhxQzMwsEw4oZmaWCQcUMzPLhAOKmZllwgHFzMwy4YBiZmaZcEAxM7NMOKCYmVkmHFDMzCwTDihmZpYJBxQzM8uEA4qZmWXCAcXMzDLhgGJmZplwQDEzs0w4oJiZWSYGDSiSHpHUJen1srQJkjZJ6kzP48teWy5pj6TdkuaXpc+WtDO9tkqSUvoYSU+l9K2SZpaVWZzO0SlpcVn6rJS3M5U9ceRdYWZmI1HNCOUHwIKKtGXA5og4G9icfkbSucAi4LxU5ruSTkhlHgCWAGenR1+dNwAHI+JLwH3AvamuCcCdwIXAHODOssB1L3BfOv/BVIeZmTXQoAElIn4KdFckLwRWp+PVwJVl6Wsioici3gL2AHMkTQFOiYgtERHAoxVl+up6GrgsjV7mA5siojsiDgKbgAXptUtT3srzm5lZgwx3DeWMiNgPkJ4np/RpwDtl+faltGnpuDL9mDIR0Qt8CEwcoK6JwAcpb2VdnyNpiaQOSR0HDhwY4j/TzKzYug4d4VsPbqHr8JGanyvrRXn1kxYDpA+nzEB1ff6FiIcioj0i2idNmnS8bGZmTWnV5k62v93Nquc7a36utmGWe0/SlIjYn6azulL6PuDMsnzTgXdT+vR+0svL7JPUBpxKaYptH3BJRZmXgPeB0yS1pVFKeV1mZgacs2IDPb1HP/358a17eXzrXsa0jWL3yq/V5JzDHaGsB/p2XS0GnilLX5R2bs2itPi+LU2LHZZ0UVoDub6iTF9dVwEvpHWWjcA8SePTYvw8YGN67cWUt/L8ZmYGvHz7XK44fypjR5fe5seOHsXC86fy8tK5NTvnoCMUSU9SGimcLmkfpZ1X9wBrJd0A7AWuBoiIXZLWAm8AvcAtEfFJquomSjvGTgI2pAfAw8BjkvZQGpksSnV1S7ob2J7y3RURfZsDlgJrJK0EXkl1mJlZMvmUsYwb00ZP71HGtI2ip/co48a0MXnc2JqdU6UP/K2hvb09Ojo6Gt0MM7O6uPGxDiaNG8s1c2bwxLa9HDh8hAevax9yPZJ2RMSgBR1QzMxsQNUGFN96xczMMuGAYmZmmXBAMTOzTDigmJlZJhxQzMwsEw4oZmaWiZbaNizpAPBP/bx0OqVbuhRFkdpbpLZCsdpbpLZCsdpbpLZC7dt7VkQMejPElgooxyOpo5o91nlRpPYWqa1QrPYWqa1QrPYWqa2Qn/Z6ysvMzDLhgGJmZplwQCl5qNENGKIitbdIbYVitbdIbYVitbdIbYWctNdrKGZmlgmPUMzMLBNNG1AkPSKpS9LrZWn/StIWSTsl/W9Jp6T00ZJWp/Q3JS0vK/OSpN2SXk2PyQ1u64mSvp/SX5N0SVmZ2Sl9j6RV6cvMMpdhe+vRt2dKejH9v+6SdGtKnyBpk6TO9Dy+rMzy1Ie7Jc0vS69p/2bc1tz1raSJKf9Hku6vqCtXfTtIW/PYt5dL2pH6cIekS8vqqsv7AgAR0ZQP4N8Avwu8Xpa2Hfi36fiPgbvT8TXAmnT8W8DbwMz080tAe47aegvw/XQ8GdgBjEo/bwN+DxClLzD7Ws7bW4++nQL8bjoeB/wjcC7w58CylL4MuDcdnwu8BowBZgG/BE6oR/9m3NY89u3JwO8D3wbur6grb307UFvz2LcXAFPT8VeAf65X35Y/mnaEEhE/pfQNkOXOAX6ajjcBf9iXHThZpe+0Pwn4DXCoHu2EIbf1XGBzKtcFfAC0S5oCnBIRW6L0W/QocGVe21uLdvUnIvZHxD+k48PAm8A0YCGwOmVbzWd9tZDSh4ueiHgL2APMqUf/ZtXWLNuUZXsj4tcR8XfAkfJ68ti3x2trvQyjva9ExLspfRcwVqWvYq/b+wI08ZTXcbwOXJGOrwbOTMdPA78G9lP6SuP/GZ993TDA99PQ9r/WdLhYXVtfAxZKapM0C5idXpsG7Csrvy+l1ctQ29unbn0raSalT3JbgTMiYj+U/ngpjZ6g1GfvlBXr68e69u8I29onb317PHns28HkuW//EHglInqoc9+2WkD5Y+AWSTsoDSN/k9LnAJ8AUylNHfyppH+RXrs2Ir4K/Ov0uK7BbX2E0i9FB/CXwP8BeikNZyvVcwvfUNsLdexbSV8A/gb4k4gYaPR5vH6sW/9m0FbIZ98et4p+0hrdtwPJbd9KOg+4F7ixL6mfbDV7X2ipgBIRv4iIeRExG3iS0pwzlNZQfhwRH6dpmb8nTctExD+n58PAE9RpSuF4bY2I3oj4TxFxfkQsBE4DOim9aU8vq2I68G5lvTlqb936VtJoSn+UP4yIH6Xk99J0QN+US1dK38exI6i+fqxL/2bU1rz27fHksW+PK699K2k6sA64PiL63tvq+r7QUgGlbzeGpFHACuB76aW9wKUqORm4CPhFmqY5PZUZDXyD0tROw9oq6bdSG5F0OdAbEW+k4e9hSRelIfj1wDP1aOtw2luvvk198TDwZkT8RdlL64HF6Xgxn/XVemBRmn+eBZwNbKtH/2bV1hz3bb9y2rfHqyeXfSvpNOBZYHlE/H1f5rq/L2S9yp+XB6VPyfuBjylF6RuAWyntlvhH4B4+u7DzC8BfU1rMegP4z/HZTo8dwM/Ta/+LtIumgW2dCeymtEj3PKW7gPbV007pl/uXwP19ZfLY3jr27e9TGuL/HHg1Pf49MJHSZoHO9DyhrMwdqQ93U7Yjptb9m1Vbc963b1Pa0PFR+t05N8d9+7m25rVvKX2I+3VZ3leByfXo2/KHr5Q3M7NMtNSUl5mZ1Y4DipmZZcIBxczMMuGAYmZmmXBAMTOzTDigmJlZJhxQzMwsEw4oZmaWif8P5AS727RchlQAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.plot(style='*')" ] @@ -314,9 +2388,56 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "2021 938731\n", + "2014 1600941\n", + "1991 1659249\n", + "1995 1840410\n", + "2020 2053781\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", + "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", + "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": 14, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,9 +2452,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 15, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEKCAYAAAAyx7/DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFu1JREFUeJzt3X2UZHV95/H31+FpoGEgQRodiGPUuCIjII0GXbUbXYNAko3hiIpGPJrW3YhEJ2cz8RhZTdxgzGSP2WjMaIweH+hjEHYTxpC4ug1iDNKDDwOixIUxMIQhgIz2MCsOfPeP3+10t/bD7emqrv51v1/n9Jl6uHXrW9+p+tS9v3tv3chMJEn1eEyvC5AkLYzBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSarMQd2Y6bHHHpsbNmzoyLz27t3LEUcc0ZF51c5eTLIX09mPSbX2Yvv27fdl5mPbTNuV4N6wYQNjY2Mdmdfo6CiDg4MdmVft7MUkezGd/ZhUay8i4rttp3WoRJIqY3BLUmUMbkmqjMEtSZUxuCWpMq2COyLeEhG3RMTNEXF5RBzW7cIkSTObN7gjYj3wZmAgM08G1gAv73ZhkqSZtR0qOQhYGxEHAYcDd3evJEnSXKLNOScj4hLg3cA+4O8z88IZphkGhgH6+/tPHxkZ6UiB4+Pj9PX1dWRetVtJvdixa8+iHt+/FnbvW/jjNq5ft6jnXa5W0ntjsWrtxdDQ0PbMHGgz7bzBHRHHAJ8BLgAeBP4KuCIzPzHbYwYGBtIjJztvJfViw+Zti3r8po372bJj4Qf+7rzs3EU973K1kt4bi1VrLyKidXC3GSp5EXBHZv5rZv4IuBJ4zmIKlCQduDbB/c/Az0fE4RERwAuBW7tbliRpNvMGd2beAFwB3ATsaB6ztct1SZJm0WqQMDMvBS7tci2SpBY8clKSKmNwS1JlDG5JqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqM29wR8RTI+JrU/6+HxG/uRTFSZJ+0rynLsvMbwOnAkTEGmAXcFWX65IkzWKhQyUvBP5vZn63G8VIkuYXmdl+4oiPADdl5p/OcN8wMAzQ399/+sjISEcKHB8fp6+vryPzqt1K6sWOXXsW9fj+tbB738Ift3H9ukU973K1kt4bi1VrL4aGhrZn5kCbaVsHd0QcAtwNPD0zd8817cDAQI6NjbWa73xGR0cZHBzsyLxqt5J6sWHztkU9ftPG/WzZMe9I30/Yedm5i3re5WolvTcWq9ZeRETr4F7IUMlLKEvbc4a2JKm7FhLcrwAu71YhkqR2WgV3RBwO/Afgyu6WI0maT6tBwsx8CPjpLtciSWrBIyclqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMm1PXXZ0RFwREd+KiFsj4sxuFyZJmlmrU5cB7wOuyczzI+IQ4PAu1iRJmsO8wR0RRwHPBy4CyMyHgYe7W5YkaTaRmXNPEHEqsBX4JnAKsB24JDP3/th0w8AwQH9//+kjIyMdKXB8fJy+vr6OzKt2ne7Fjl17Ojavpda/Fnbv63UV7W1cv66r8/dzMqnWXgwNDW3PzIE207YJ7gHgH4HnZuYNEfE+4PuZ+buzPWZgYCDHxsYWUvOsRkdHGRwc7Mi8atfpXmzYvK1j81pqmzbuZ8uOtiN9vbfzsnO7On8/J5Nq7UVEtA7uNhsn7wLuyswbmutXAM880OIkSYszb3Bn5j3AnRHx1OamF1KGTSRJPdB2XfNi4JPNHiW3A6/tXkmSpLm0Cu7M/BrQauxFktRdHjkpSZUxuCWpMga3JFXG4JakyhjcklQZg1uSKmNwS1JlDG5JqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlWl1BpyI2An8AHgE2N/2TMSSpM5re85JgKHMvK9rlUiSWnGoRJIqE5k5/0QRdwDfAxL488zcOsM0w8AwQH9//+kjIyMdKXB8fJy+vr6OzKt2ne7Fjl17Ojavpda/Fnbv63UV7W1cv66r8/dzMqnWXgwNDW1vOwzdNrgfn5l3R8RxwOeAizPzutmmHxgYyLGxsdYFz2V0dJTBwcGOzKt2ne7Fhs3bOjavpbZp43627FjISF9v7bzs3K7O38/JpFp7ERGtg7vVUElm3t38ey9wFfCsAy9PkrQY8wZ3RBwREUdOXAZeDNzc7cIkSTNrs67ZD1wVERPTfyozr+lqVZKkWc0b3Jl5O3DKEtQiSWrB3QElqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMq2DOyLWRMRXI+LqbhYkSZrbQpa4LwFu7VYhkqR2WgV3RJwAnAt8uLvlSJLmE5k5/0QRVwB/ABwJ/FZmnjfDNMPAMEB/f//pIyMjHSlwfHycvr6+jsxrudqxa0+r6frXwu59XS6mErX1YuP6dV2d/2r4nLRVay+Ghoa2Z+ZAm2nnPct7RJwH3JuZ2yNicLbpMnMrsBVgYGAgBwdnnXRBRkdH6dS8lquLNm9rNd2mjfvZsmPe/7JVobZe7LxwsKvzXw2fk7ZWQy/aDJU8F/iliNgJjABnRcQnulqVJGlW8wZ3Zv5OZp6QmRuAlwNfyMxXdb0ySdKM3I9bkiqzoEHCzBwFRrtSiSSpFZe4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMga3JFXG4JakyhjcklQZg1uSKmNwS1JlDG5JqozBLUmVmTe4I+KwiPhKRHw9Im6JiHcuRWGSpJm1OQPOD4GzMnM8Ig4Gro+Iv83Mf+xybZKkGcwb3JmZwHhz9eDmL7tZlCRpdlFyeZ6JItYA24EnA+/PzN+eYZphYBigv7//9JGRkY4UOD4+Tl9fX0fmNZcdu/Z0/TkWq38t7N7X6yqWB3sx3Vz92Lh+3dIW02Pj4+PcseeRnjz3Yno9NDS0PTMH2kzbKrj/beKIo4GrgIsz8+bZphsYGMixsbHW853L6Ogog4ODHZnXXDZs3tb151isTRv3s2XHgs7vvGLZi+nm6sfOy85d4mp6a3R0lIuu2duT515MryOidXAvaK+SzHyQcpb3sw+gLklSB7TZq+SxzZI2EbEWeBHwrW4XJkmaWZt1zccBH2vGuR8DfDozr+5uWZKk2bTZq+QbwGlLUIskqQWPnJSkyhjcklQZg1uSKmNwS1JlDG5JqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKtDnn5IkR8X8i4taIuCUiLlmKwiRJM2tzzsn9wKbMvCkijgS2R8TnMvObXa5NkjSDeZe4M/NfMvOm5vIPgFuB9d0uTJI0s8jM9hNHbACuA07OzO//2H3DwDBAf3//6SMjIx0pcHx8nL6+vo7May47du3p+nMsVv9a2L2v11UsD/ZiuuXYj43r1/XkecfHx7ljzyM9ee7FvOahoaHtmTnQZtrWwR0RfcC1wLsz88q5ph0YGMixsbFW853P6Ogog4ODHZnXXDZs3tb151isTRv3s2VHm9Gtlc9eTLcc+7HzsnN78ryjo6NcdM3enjz3Yl5zRLQO7lZ7lUTEwcBngE/OF9qSpO5qs1dJAH8B3JqZf9z9kiRJc2mzxP1c4NXAWRHxtebvnC7XJUmaxbyDYpl5PRBLUIskqQWPnJSkyhjcklQZg1uSKmNwS1JlDG5JqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKtDnn5Eci4t6IuHkpCpIkza3NEvdHgbO7XIckqaV5gzszrwMeWIJaJEktRGbOP1HEBuDqzDx5jmmGgWGA/v7+00dGRg6ooB279ky73r8Wdu87oFmtOPZikr2Ybjn2Y+P6dT153vHxce7Y80hPnnsxr3loaGh7Zg60mXbes7y3lZlbga0AAwMDOTg4eEDzuWjztmnXN23cz5YdHSuzavZikr2Ybjn2Y+eFgz153tHRUbZcv7cnz71Ur9m9SiSpMga3JFWmze6AlwNfBp4aEXdFxOu6X5YkaTbzDopl5iuWohBJUjsOlUhSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMga3JFXG4JakyhjcklQZg1uSKmNwS1JlWgV3RJwdEd+OiO9ExOZuFyVJml2bc06uAd4PvAQ4CXhFRJzU7cIkSTNrs8T9LOA7mXl7Zj4MjAC/3N2yJEmzicyce4KI84GzM/P1zfVXA8/OzDf92HTDwHBz9anAtztU47HAfR2aV+3sxSR7MZ39mFRrL56QmY9tM+G8Z3kHYobbfiLtM3MrsLXNky5ERIxl5kCn51sjezHJXkxnPyathl60GSq5CzhxyvUTgLu7U44kaT5tgvtG4CkR8cSIOAR4OfDX3S1LkjSbeYdKMnN/RLwJ+DtgDfCRzLyl65VN6vjwS8XsxSR7MZ39mLTiezHvxklJ0vLikZOSVBmDW5IqY3BLUmVWdHBHxPqIWN/rOpaDiPjZiHhLRJzV61p6zV5MshfT1dKPFRncEbEhIq4FrgHeGxHP63VNvRQR/x74HOW3Zt4YEf+pxyX1jL2YZC+mq6kfKya4I+KwKVefCdyYmRspuzG+OSI2NtPNdCToihIRZ0XEE5vLAZwFXJqZvw5sAc6JiMEp969Y9mKSvZiu5n5UHdwRcVREfDAibgP+KCKe0Nz1K8A/N5dHgO8Ar5942BKXuWQi4qSI+AbwX4G/jIizsuzveRJwPEBm3gD8A/DaiYf1otZusxeT7MV0K6EfVQc3cDZwGKXhDwPviIi1lNWdXwTIzB8CVwDPa64/2ptSOy8iToiIo6bcdAHwmcx8PuUL65UR8RTgUzT9aFwFnBwRh66UftiLSfZiupXYj2Uf3FEcFBGvi4gvRsQlEfGk5u4nAw9n5n7gvwPfA14F/D3wuIj4qWa624A7I+LMJX8BXRART4uIzwLXA++KiImf2f1/wOHN5U8D9wDnUpYcfnrKGskDlF9vPGXpqu4OezHJXky3kvux7IO7WYV5AfBrwB8ChwIfau6+B7i3+Ua8k9LkJ1H+Y77J5M/MHgzc39xepYg4YsrVU4G7MnMD8AXgj5rbHwB+GBFHZuYDwD8Bj6e87n8A3tpMdwjwCLCz+5V3nr2YZC+mWy39WHbBHRFnRsR7IuKi5noATwOuycy/ycw/BJ4QEc8BdlG+OZ/WPPxWoK+57U8pGxd+kRL6/cDXl/TFLFJEHBMRH42IG4HLIuKxTT82Al+KiMjMvwYejIhzKWsWRzb301w/DniUskZyXER8CLgc2J+Z9y71azpQ9mKSvZhuNfZjWQV3RDwd+DPgB8DLIuKtlBrXAz+YsmX3o8ArKUG8H3hOc/tNlC3DD2XmdcBm4CLgucDvZeajy23r8DyeT3l951A2jrwNOIryY1/H5+QPzXyM0o+vUHr3EoDM/HIzj4My81bgDcAtwH/LzNdSF3sxyV5Mt/r6kZk9+aMsGb+esjpzUHPbHwOXNJcHgD8BzgdeBPzdlMeeSFkFghLUX6Wcdec04H8Bj5sybfTqNS6gF2sob5ZrKcM7xza3fxp4c3P5icBlzf1nUMbt1kzp5b8281lPWfN4E/CXwAeAI3r9Gu2FvbAfnfvryRJ3RJxC2YD4y8ClwNubu3ZRznEJ5RvvS8CvAp8Hjo+IZ0TEwVnGs3dFxPMy8wuUn3F8D3AlcHlm/svEc2XzP7bMnQf8EvBO4EzKWD6UvWMm1ibuBL4IvCQzb6QsWQwBZOY4cANwRmbuAl5NGRq6B3h7Zu5dotfRCfZikr2Yzn402py6bNEi4jRgb2be1tz0LOC2zLwoIp4JvDsiBoBR4Bci4vDMfCgivg68jLJv5aeAXwf+JCL2ATuAO5r5fRD4VGbuWYrXcyCacbaMiDMoq2tfBLZl2V3x54DbM/MLEXEH5WjPFwPbgV+JiGMz876I+Cdgb0T8DPA/gFdFxHGUsxLdT1kFJDPHgLElf5Et2YtJ9mI6+9FOV5e4I+LJEXE9ZUz6HRHxa81djwI7m6XnmyirLGcCDzG5aw7AjyirNsdTlqpvbuZ1LXBfZt4FZam6ktB+PvARytbrFwF/0EzyKHBbRKzNzDso/XgGZRzubsp+p1C2cK+h/L99htKTC4HTga25zPY1nUlErGl68QLKKupq7sXapheDrPL3BZQD6uxHS50cdwGOAM6ccv084H3N5WdTvt2eALyGsrqzvrnvfMp49sR91za3H0YZJjl2yjxPAw7p1djSAnpxOPBGJtcUDgZ+E/iN5v5jgG80r+cCyrjchil920o5W/V5lLWLdZTx/M9Off3AY3r9Wlu+L15P+RBtomw4WpW9aOo8EthGOZsUwFtWcS8Obz7zn6ccFLOq+9H2r2NL3BHxNuB24LMR0d/c/AuUfavJcgjpV4CLKftUnkjZ5xrKWPaplL1BPgZ8LyI+Ttno+G3g38aeMvOrmflwp+ruhog4HrgaGAQ+TtlQ8lLKWsV+gMz8HmVD6pspY3THMblb43WUfdcfzsyrgb+gHP35fsqW8R9NPFcu86WHZr/az1M+TB8CXkzZbnEGZQlq1fRiirWU4xGeFBHHUj4Ha2B19SIiDqZsyzofeG9m/mpz12kT06ymfixIB785BymrLR8GNjW3vZUyPjUxzUnAnc3ldwPvmnLfjcBpzeVDKbv2nNHrb7YD7MVa4NlTrl9E2ZDyGuArU25/PHB3c/k3KIffHtM8/m+An5ky7bFLUXuX+nH0lMv/hfIhvHA19qKp/zXAe4HfBV5HOcz6xlXaiyuBC3/stguAG1ZjP1r3rYP/ARO73FzA5FDH0cCDwGFTpruR8o16NPBXlFWdv6V8Qx7a64Z0qBcx8ddcf+aUntxP2bd0YtrPTYQ88PuUvW3uB36716+jwz05irJ9Yjfwrub6/UD/aunFlPfDaynDaC8FPtncdt9q6sWU13ke5QCYLZSdE95BGTJ9ADhutfWj7V/Hhkoy85Hm4v8GjoyIUzLzQcq49humTHoTcGRz38WU4ZD/CQxn2XJcvWxMuekSypIFlLG3twBE+S2V7wITuy++k7KWsj4z37NE5S6JzPw+Zajs5ykbm19GGQZ7QxQrvhdT3hPnUIaNPg+cEBFvp2yUH4ZV9764mrJ32P2U3fOeDvxHynvjjavlvbFQXTnLe0R8gDJe/VvNnhT/mRLgx1AOrDlnStCvaBFxAmXs7eLMvC3KD2QNU96g64Gv5XI9OqtLIuJUypf5lynjlSdTdtVa8b2IiD7KMMmhlNf+7ygHibyNsiT+FFZJLyZM7P7bXD6F8vn4EuWQ9FXz3liIbgX3qZS9RF5IeXM+RDn8fB/wwcy8peNPukxF+a2UFwC/QxnbvIuy2ncB8K0su0OuKhFxIuXL7BWZeX9EvAq4JTO/2uPSui7KCT/+nLLh7HLKrmtvy8wXN/evml7MJMov830YuCAzH1jt/ZhNt4L75ZTd4B4Cfo+yxXhlbdVtKSK+BPws5RfG7gbemZnf6GlRPRAR6yhf5K+kbKTeCrw/M3805wNXuOYgkZcCI5l5T6/r6YWIOJTy2/oTQyV/Bnwgy881awYdD+6IeAbl8PMrKBteqv0p1cVqdne6lDKG94mVMoZ/ICLiIMrwyA8pvVi17wsoByIBj2Y3lpwqFBFvoOwe+vHV/t5ooytL3JKk7llWP+sqSZqfwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5Iq8/8BnVLOiAvrP5IAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.hist(xrot=20)" ] @@ -341,9 +2485,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": true - }, + "metadata": {}, "outputs": [], "source": [] } @@ -364,7 +2506,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4,