diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5b58a3ae26346460dd39e62a39c55243d7..1e86b3e8d02da089f30c3e43f8d90b23e6c99d08 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": 8, "metadata": {}, "outputs": [], "source": [ @@ -28,10 +28,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 9, + "metadata": {}, "outputs": [], "source": [ "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" @@ -61,11 +59,978 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \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
020211931404710193.017901.02115.027.0FRFrance
12021183126509582.015718.01914.024.0FRFrance
22021173120588891.015225.01813.023.0FRFrance
320211631650512735.020275.02519.031.0FRFrance
420211531930615398.023214.02923.035.0FRFrance
520211432107317099.025047.03226.038.0FRFrance
620211332641322094.030732.04033.047.0FRFrance
720211233065825919.035397.04639.053.0FRFrance
820211132498820718.029258.03832.044.0FRFrance
920211031953915951.023127.03025.035.0FRFrance
1020210931757213926.021218.02721.033.0FRFrance
1120210832088216907.024857.03226.038.0FRFrance
1220210732239318303.026483.03428.040.0FRFrance
1320210632318319134.027232.03529.041.0FRFrance
1420210532242618445.026407.03428.040.0FRFrance
1520210432580421491.030117.03932.046.0FRFrance
1620210332181017894.025726.03327.039.0FRFrance
1720210231732013906.020734.02621.031.0FRFrance
1820210132179917778.025820.03327.039.0FRFrance
1920205332122016498.025942.03225.039.0FRFrance
2020205231642812285.020571.02519.031.0FRFrance
2120205132161917370.025868.03327.039.0FRFrance
2220205031684513220.020470.02620.032.0FRFrance
232020493129399923.015955.02015.025.0FRFrance
2420204831380410641.016967.02116.026.0FRFrance
2520204731908515285.022885.02923.035.0FRFrance
2620204632480120503.029099.03831.045.0FRFrance
2720204534251636857.048175.06556.074.0FRFrance
2820204434456738521.050613.06859.077.0FRFrance
2920204334373737523.049951.06657.075.0FRFrance
.................................
187719852132609619621.032571.04735.059.0FRFrance
187819852032789620885.034907.05138.064.0FRFrance
187919851934315432821.053487.07859.097.0FRFrance
188019851834055529935.051175.07455.093.0FRFrance
188119851733405324366.043740.06244.080.0FRFrance
188219851635036236451.064273.09166.0116.0FRFrance
188319851536388145538.082224.011683.0149.0FRFrance
18841985143134545114400.0154690.0244207.0281.0FRFrance
18851985133197206176080.0218332.0357319.0395.0FRFrance
18861985123245240223304.0267176.0445405.0485.0FRFrance
18871985113276205252399.0300011.0501458.0544.0FRFrance
18881985103353231326279.0380183.0640591.0689.0FRFrance
18891985093369895341109.0398681.0670618.0722.0FRFrance
18901985083389886359529.0420243.0707652.0762.0FRFrance
18911985073471852432599.0511105.0855784.0926.0FRFrance
18921985063565825518011.0613639.01026939.01113.0FRFrance
18931985053637302592795.0681809.011551074.01236.0FRFrance
18941985043424937390794.0459080.0770708.0832.0FRFrance
18951985033213901174689.0253113.0388317.0459.0FRFrance
189619850239758680949.0114223.0177147.0207.0FRFrance
189719850138548965918.0105060.0155120.0190.0FRFrance
189819845238483060602.0109058.0154110.0198.0FRFrance
1899198451310172680242.0123210.0185146.0224.0FRFrance
19001984503123680101401.0145959.0225184.0266.0FRFrance
1901198449310107381684.0120462.0184149.0219.0FRFrance
190219844837862060634.096606.0143110.0176.0FRFrance
190319844737202954274.089784.013199.0163.0FRFrance
190419844638733067686.0106974.0159123.0195.0FRFrance
19051984453135223101414.0169032.0246184.0308.0FRFrance
190619844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1907 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202119 3 14047 10193.0 17901.0 21 15.0 \n", + "1 202118 3 12650 9582.0 15718.0 19 14.0 \n", + "2 202117 3 12058 8891.0 15225.0 18 13.0 \n", + "3 202116 3 16505 12735.0 20275.0 25 19.0 \n", + "4 202115 3 19306 15398.0 23214.0 29 23.0 \n", + "5 202114 3 21073 17099.0 25047.0 32 26.0 \n", + "6 202113 3 26413 22094.0 30732.0 40 33.0 \n", + "7 202112 3 30658 25919.0 35397.0 46 39.0 \n", + "8 202111 3 24988 20718.0 29258.0 38 32.0 \n", + "9 202110 3 19539 15951.0 23127.0 30 25.0 \n", + "10 202109 3 17572 13926.0 21218.0 27 21.0 \n", + "11 202108 3 20882 16907.0 24857.0 32 26.0 \n", + "12 202107 3 22393 18303.0 26483.0 34 28.0 \n", + "13 202106 3 23183 19134.0 27232.0 35 29.0 \n", + "14 202105 3 22426 18445.0 26407.0 34 28.0 \n", + "15 202104 3 25804 21491.0 30117.0 39 32.0 \n", + "16 202103 3 21810 17894.0 25726.0 33 27.0 \n", + "17 202102 3 17320 13906.0 20734.0 26 21.0 \n", + "18 202101 3 21799 17778.0 25820.0 33 27.0 \n", + "19 202053 3 21220 16498.0 25942.0 32 25.0 \n", + "20 202052 3 16428 12285.0 20571.0 25 19.0 \n", + "21 202051 3 21619 17370.0 25868.0 33 27.0 \n", + "22 202050 3 16845 13220.0 20470.0 26 20.0 \n", + "23 202049 3 12939 9923.0 15955.0 20 15.0 \n", + "24 202048 3 13804 10641.0 16967.0 21 16.0 \n", + "25 202047 3 19085 15285.0 22885.0 29 23.0 \n", + "26 202046 3 24801 20503.0 29099.0 38 31.0 \n", + "27 202045 3 42516 36857.0 48175.0 65 56.0 \n", + "28 202044 3 44567 38521.0 50613.0 68 59.0 \n", + "29 202043 3 43737 37523.0 49951.0 66 57.0 \n", + "... ... ... ... ... ... ... ... \n", + "1877 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1878 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1879 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1880 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1881 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1882 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1883 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1884 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1885 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1886 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1887 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1888 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1889 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1890 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1891 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1892 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1893 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1894 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1895 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1896 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1897 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1898 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1899 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1900 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1901 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1902 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1903 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1904 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1905 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1906 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 27.0 FR France \n", + "1 24.0 FR France \n", + "2 23.0 FR France \n", + "3 31.0 FR France \n", + "4 35.0 FR France \n", + "5 38.0 FR France \n", + "6 47.0 FR France \n", + "7 53.0 FR France \n", + "8 44.0 FR France \n", + "9 35.0 FR France \n", + "10 33.0 FR France \n", + "11 38.0 FR France \n", + "12 40.0 FR France \n", + "13 41.0 FR France \n", + "14 40.0 FR France \n", + "15 46.0 FR France \n", + "16 39.0 FR France \n", + "17 31.0 FR France \n", + "18 39.0 FR France \n", + "19 39.0 FR France \n", + "20 31.0 FR France \n", + "21 39.0 FR France \n", + "22 32.0 FR France \n", + "23 25.0 FR France \n", + "24 26.0 FR France \n", + "25 35.0 FR France \n", + "26 45.0 FR France \n", + "27 74.0 FR France \n", + "28 77.0 FR France \n", + "29 75.0 FR France \n", + "... ... ... ... \n", + "1877 59.0 FR France \n", + "1878 64.0 FR France \n", + "1879 97.0 FR France \n", + "1880 93.0 FR France \n", + "1881 80.0 FR France \n", + "1882 116.0 FR France \n", + "1883 149.0 FR France \n", + "1884 281.0 FR France \n", + "1885 395.0 FR France \n", + "1886 485.0 FR France \n", + "1887 544.0 FR France \n", + "1888 689.0 FR France \n", + "1889 722.0 FR France \n", + "1890 762.0 FR France \n", + "1891 926.0 FR France \n", + "1892 1113.0 FR France \n", + "1893 1236.0 FR France \n", + "1894 832.0 FR France \n", + "1895 459.0 FR France \n", + "1896 207.0 FR France \n", + "1897 190.0 FR France \n", + "1898 198.0 FR France \n", + "1899 224.0 FR France \n", + "1900 266.0 FR France \n", + "1901 219.0 FR France \n", + "1902 176.0 FR France \n", + "1903 163.0 FR France \n", + "1904 195.0 FR France \n", + "1905 308.0 FR France \n", + "1906 213.0 FR France \n", + "\n", + "[1907 rows x 10 columns]" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "raw_data = pd.read_csv(data_url, skiprows=1)\n", + "raw_data = pd.read_csv(\"incidence-PAY-3.csv\", skiprows=1)\n", "raw_data" ] }, @@ -78,9 +1043,73 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
167019891930NaNNaN0NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1670 198919 3 0 NaN NaN 0 NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1670 FR France " + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -94,9 +1123,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "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
020211931404710193.017901.02115.027.0FRFrance
12021183126509582.015718.01914.024.0FRFrance
22021173120588891.015225.01813.023.0FRFrance
320211631650512735.020275.02519.031.0FRFrance
420211531930615398.023214.02923.035.0FRFrance
520211432107317099.025047.03226.038.0FRFrance
620211332641322094.030732.04033.047.0FRFrance
720211233065825919.035397.04639.053.0FRFrance
820211132498820718.029258.03832.044.0FRFrance
920211031953915951.023127.03025.035.0FRFrance
1020210931757213926.021218.02721.033.0FRFrance
1120210832088216907.024857.03226.038.0FRFrance
1220210732239318303.026483.03428.040.0FRFrance
1320210632318319134.027232.03529.041.0FRFrance
1420210532242618445.026407.03428.040.0FRFrance
1520210432580421491.030117.03932.046.0FRFrance
1620210332181017894.025726.03327.039.0FRFrance
1720210231732013906.020734.02621.031.0FRFrance
1820210132179917778.025820.03327.039.0FRFrance
1920205332122016498.025942.03225.039.0FRFrance
2020205231642812285.020571.02519.031.0FRFrance
2120205132161917370.025868.03327.039.0FRFrance
2220205031684513220.020470.02620.032.0FRFrance
232020493129399923.015955.02015.025.0FRFrance
2420204831380410641.016967.02116.026.0FRFrance
2520204731908515285.022885.02923.035.0FRFrance
2620204632480120503.029099.03831.045.0FRFrance
2720204534251636857.048175.06556.074.0FRFrance
2820204434456738521.050613.06859.077.0FRFrance
2920204334373737523.049951.06657.075.0FRFrance
.................................
187719852132609619621.032571.04735.059.0FRFrance
187819852032789620885.034907.05138.064.0FRFrance
187919851934315432821.053487.07859.097.0FRFrance
188019851834055529935.051175.07455.093.0FRFrance
188119851733405324366.043740.06244.080.0FRFrance
188219851635036236451.064273.09166.0116.0FRFrance
188319851536388145538.082224.011683.0149.0FRFrance
18841985143134545114400.0154690.0244207.0281.0FRFrance
18851985133197206176080.0218332.0357319.0395.0FRFrance
18861985123245240223304.0267176.0445405.0485.0FRFrance
18871985113276205252399.0300011.0501458.0544.0FRFrance
18881985103353231326279.0380183.0640591.0689.0FRFrance
18891985093369895341109.0398681.0670618.0722.0FRFrance
18901985083389886359529.0420243.0707652.0762.0FRFrance
18911985073471852432599.0511105.0855784.0926.0FRFrance
18921985063565825518011.0613639.01026939.01113.0FRFrance
18931985053637302592795.0681809.011551074.01236.0FRFrance
18941985043424937390794.0459080.0770708.0832.0FRFrance
18951985033213901174689.0253113.0388317.0459.0FRFrance
189619850239758680949.0114223.0177147.0207.0FRFrance
189719850138548965918.0105060.0155120.0190.0FRFrance
189819845238483060602.0109058.0154110.0198.0FRFrance
1899198451310172680242.0123210.0185146.0224.0FRFrance
19001984503123680101401.0145959.0225184.0266.0FRFrance
1901198449310107381684.0120462.0184149.0219.0FRFrance
190219844837862060634.096606.0143110.0176.0FRFrance
190319844737202954274.089784.013199.0163.0FRFrance
190419844638733067686.0106974.0159123.0195.0FRFrance
19051984453135223101414.0169032.0246184.0308.0FRFrance
190619844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1906 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202119 3 14047 10193.0 17901.0 21 15.0 \n", + "1 202118 3 12650 9582.0 15718.0 19 14.0 \n", + "2 202117 3 12058 8891.0 15225.0 18 13.0 \n", + "3 202116 3 16505 12735.0 20275.0 25 19.0 \n", + "4 202115 3 19306 15398.0 23214.0 29 23.0 \n", + "5 202114 3 21073 17099.0 25047.0 32 26.0 \n", + "6 202113 3 26413 22094.0 30732.0 40 33.0 \n", + "7 202112 3 30658 25919.0 35397.0 46 39.0 \n", + "8 202111 3 24988 20718.0 29258.0 38 32.0 \n", + "9 202110 3 19539 15951.0 23127.0 30 25.0 \n", + "10 202109 3 17572 13926.0 21218.0 27 21.0 \n", + "11 202108 3 20882 16907.0 24857.0 32 26.0 \n", + "12 202107 3 22393 18303.0 26483.0 34 28.0 \n", + "13 202106 3 23183 19134.0 27232.0 35 29.0 \n", + "14 202105 3 22426 18445.0 26407.0 34 28.0 \n", + "15 202104 3 25804 21491.0 30117.0 39 32.0 \n", + "16 202103 3 21810 17894.0 25726.0 33 27.0 \n", + "17 202102 3 17320 13906.0 20734.0 26 21.0 \n", + "18 202101 3 21799 17778.0 25820.0 33 27.0 \n", + "19 202053 3 21220 16498.0 25942.0 32 25.0 \n", + "20 202052 3 16428 12285.0 20571.0 25 19.0 \n", + "21 202051 3 21619 17370.0 25868.0 33 27.0 \n", + "22 202050 3 16845 13220.0 20470.0 26 20.0 \n", + "23 202049 3 12939 9923.0 15955.0 20 15.0 \n", + "24 202048 3 13804 10641.0 16967.0 21 16.0 \n", + "25 202047 3 19085 15285.0 22885.0 29 23.0 \n", + "26 202046 3 24801 20503.0 29099.0 38 31.0 \n", + "27 202045 3 42516 36857.0 48175.0 65 56.0 \n", + "28 202044 3 44567 38521.0 50613.0 68 59.0 \n", + "29 202043 3 43737 37523.0 49951.0 66 57.0 \n", + "... ... ... ... ... ... ... ... \n", + "1877 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1878 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1879 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1880 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1881 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1882 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1883 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1884 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1885 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1886 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1887 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1888 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1889 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1890 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1891 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1892 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1893 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1894 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1895 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1896 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1897 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1898 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1899 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1900 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1901 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1902 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1903 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1904 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1905 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1906 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 27.0 FR France \n", + "1 24.0 FR France \n", + "2 23.0 FR France \n", + "3 31.0 FR France \n", + "4 35.0 FR France \n", + "5 38.0 FR France \n", + "6 47.0 FR France \n", + "7 53.0 FR France \n", + "8 44.0 FR France \n", + "9 35.0 FR France \n", + "10 33.0 FR France \n", + "11 38.0 FR France \n", + "12 40.0 FR France \n", + "13 41.0 FR France \n", + "14 40.0 FR France \n", + "15 46.0 FR France \n", + "16 39.0 FR France \n", + "17 31.0 FR France \n", + "18 39.0 FR France \n", + "19 39.0 FR France \n", + "20 31.0 FR France \n", + "21 39.0 FR France \n", + "22 32.0 FR France \n", + "23 25.0 FR France \n", + "24 26.0 FR France \n", + "25 35.0 FR France \n", + "26 45.0 FR France \n", + "27 74.0 FR France \n", + "28 77.0 FR France \n", + "29 75.0 FR France \n", + "... ... ... ... \n", + "1877 59.0 FR France \n", + "1878 64.0 FR France \n", + "1879 97.0 FR France \n", + "1880 93.0 FR France \n", + "1881 80.0 FR France \n", + "1882 116.0 FR France \n", + "1883 149.0 FR France \n", + "1884 281.0 FR France \n", + "1885 395.0 FR France \n", + "1886 485.0 FR France \n", + "1887 544.0 FR France \n", + "1888 689.0 FR France \n", + "1889 722.0 FR France \n", + "1890 762.0 FR France \n", + "1891 926.0 FR France \n", + "1892 1113.0 FR France \n", + "1893 1236.0 FR France \n", + "1894 832.0 FR France \n", + "1895 459.0 FR France \n", + "1896 207.0 FR France \n", + "1897 190.0 FR France \n", + "1898 198.0 FR France \n", + "1899 224.0 FR France \n", + "1900 266.0 FR France \n", + "1901 219.0 FR France \n", + "1902 176.0 FR France \n", + "1903 163.0 FR France \n", + "1904 195.0 FR France \n", + "1905 308.0 FR France \n", + "1906 213.0 FR France \n", + "\n", + "[1906 rows x 10 columns]" + ] + }, + "execution_count": 12, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -122,7 +2118,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, "outputs": [], "source": [ @@ -152,10 +2148,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 14, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -179,9 +2173,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n" + ] + } + ], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", @@ -199,9 +2201,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 16, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXm8HUWZ93/PWe65+83CTUIWEnYICbKEAKOigCziEh113oBKRhlRh3EYdeYd0HFwRBxwfHVARAclbCLIIAqKgJFF9iWQQEhCyCUJ2XNvcrPcJXc559T7R1f1qe7T1d3nnD7bzfP9fJLTt7qrurq7up56lqomIQQYhmEYplRi1a4AwzAMMzZggcIwDMNEAgsUhmEYJhJYoDAMwzCRwAKFYRiGiQQWKAzDMEwksEBhGIZhIoEFCsMwDBMJLFAYhmGYSEhUuwKV5KCDDhKzZs2qdjUYhmHqildeeWWnEKIz6LgDSqDMmjULS5curXY1GIZh6goieifMcWzyYhiGYSKBBQrDMAwTCSxQGIZhmEhggcIwDMNEAgsUhmEYJhICBQoRLSaibiJ6Q0ubQERLiGit/B2v7buSiLqIaA0Rnaeln0xEK+S+G4iIZHqKiH4t018kollankXyHGuJaJGWfqg8dq3M21D6rWAYhmFKIYyGchuA811pVwB4TAhxJIDH5N8gotkAFgI4Tua5iYjiMs9PAVwK4Ej5T5V5CYDdQogjAPwIwHWyrAkArgJwKoD5AK7SBNd1AH4kz79blsEwDMNUkUCBIoR4CkCvK3kBgNvl9u0APqal3yOEGBZCrAfQBWA+ER0MoF0I8bywvjl8hyuPKus+AGdL7eU8AEuEEL1CiN0AlgA4X+47Sx7rPv+Ypn84jd8t21LtajAMw3hS7MTGyUKIbQAghNhGRJNk+jQAL2jHbZZpo3Lbna7ybJJlpYloL4CJerorz0QAe4QQaY+yxjTfuH8FHnxtKw7vbMXc6R3Vrg7DMIyDqJ3y5JEmfNKLyeNXVn6FiC4loqVEtLSnp8d0WF2wafcgAGAkk6lyTRiGYfIpVqDskGYsyN9umb4ZwAztuOkAtsr06R7pjjxElADQAcvEZiprJ4Bx8lh3WXkIIW4WQswTQszr7AxciqamyWQtuZmIcXAewzC1R7E904MAVNTVIgAPaOkLZeTWobCc7y9J81gfEZ0mfSAXu/Kosj4J4HHpZ3kUwLlENF46488F8Kjc94Q81n3+MU06YwmUeMxLSWMYhqkugT4UIrobwPsBHEREm2FFXl0L4F4iugTARgCfAgAhxEoiuhfAKgBpAJcJIZR95suwIsaaADws/wHALQDuJKIuWJrJQllWLxFdDeBledx3hBAqOOBfAdxDRN8FsEyWMeZJZ7MAgEScBQrDMLVHoEARQlxo2HW24fhrAFzjkb4UwByP9CFIgeSxbzGAxR7p62CFEh9QpKXJK04sUBiGqT3YGF9HKB8KwzBMLcICpY5QAoXFCsMwtQgLlDpCWbqygkUKwzC1BwuUOoLkFBzpm2cYhqkpWKDUEUpDEWz0YhimBmGBUkeo2C62eDEMU4uwQKkj5Ir/LFAYhqlJWKDUEeyUZximlmGBUkcokxcLFIZhahEWKHUIixOGYWoRFih1RM6HwiKFYZjagwVKHcFRXgzD1DIsUOqImNRQeEkvhmFqERYodQRHeTEMU8uwQKlDal2enHz1EvzDr16tdjUYhqkwLFDqiHpxyu8aGMEfXt9W7WowDFNhWKDUEbl5KFWtBsMwjCcsUOoIXhySYZhahgVKHZFzyle3HgzDMF6wQKkj7O+h1LgPhWGYAxMWKAzDMEwksECpR1hBYRimBmGBUkewU55hmFqGBQrDMAwTCSxQ6hD2yTMMU4uwQKkjWJAwDFPLsECpQ1iwMAxTi7BAqSNyTnmGYZjagwUKwzAMEwksUOoQv9WGP3TD07jx8bUVrA3DMIxFSQKFiL5KRCuJ6A0iupuIGoloAhEtIaK18ne8dvyVRNRFRGuI6Dwt/WQiWiH33UBynXYiShHRr2X6i0Q0S8uzSJ5jLREtKuU6xhIrt+7DD/70VrWrwTDMAUjRAoWIpgH4RwDzhBBzAMQBLARwBYDHhBBHAnhM/g0imi33HwfgfAA3EVFcFvdTAJcCOFL+O1+mXwJgtxDiCAA/AnCdLGsCgKsAnApgPoCrdME11mEfCsMwtUipJq8EgCYiSgBoBrAVwAIAt8v9twP4mNxeAOAeIcSwEGI9gC4A84noYADtQojnhWXLucOVR5V1H4CzpfZyHoAlQoheIcRuAEuQE0JjHo7yYhimFilaoAghtgD4AYCNALYB2CuE+BOAyUKIbfKYbQAmySzTAGzSitgs06bJbXe6I48QIg1gL4CJPmUxDMMwVaIUk9d4WBrEoQCmAmghos/4ZfFIEz7pxeZx1/NSIlpKREt7enp8qldPsIrCMEztUYrJ6wMA1gsheoQQowDuB/BXAHZIMxbkb7c8fjOAGVr+6bBMZJvltjvdkUea1ToA9PqUlYcQ4mYhxDwhxLzOzs4iL7W2YJMXwzC1SCkCZSOA04ioWfo1zgawGsCDAFTU1SIAD8jtBwEslJFbh8Jyvr8kzWJ9RHSaLOdiVx5V1icBPC79LI8COJeIxktN6VyZxjAMw1SJRLEZhRAvEtF9AF4FkAawDMDNAFoB3EtEl8ASOp+Sx68konsBrJLHXyaEyMjivgzgNgBNAB6W/wDgFgB3ElEXLM1koSyrl4iuBvCyPO47QojeYq+l3mAFhWGYWqRogQIAQoirYIXv6gzD0la8jr8GwDUe6UsBzPFIH4IUSB77FgNYXGCVGYZhmDLBM+XrEPahMAxTi7BAqUP4i40Mw9QiLFAYhmGYSGCBUoewyYthmFqEBUodoWZzHgjy5Mr7X8fDK7ZVuxoMwxQACxSmJrn7pU348l2vVrsaDMMUAAuUOsTveygMwzDVggVKHcFihGGYWoYFCsMwDBMJLFDqCNspz6oKwzA1CAsUhmEYJhJYoNQhPFOeYZhahAVKHcImL4ZhahEWKAzDMEwksECpQ1hDYRimFmGBwjAMw0QCC5Q6xKSg8Ax6hmGqCQuUOsQkOFieMAxTTVigjCFYnjAMU01YoNQhLDgYhqlFWKCMIcaKD2WsXAfDHGiwQKlHDP3tWOmGWZ4wTH3CAqUOMS29MlY64jFyGQxzwMECZQwxVtb4YpMXw9QnLFDqEFN/O1b64TFyGQxzwMECpZ4g64soY73DDSMYB0fS2D0wUv7KMAwTGhYoBfDGlr1YsmpHtathZOxoKMEXcsH1T+PEq5dUoDYMw4QlUe0K1BMf/vEzAIAN136oqvUwmrzGiO4SRjBu2DVY/oowDFMQrKHUE2NFBQngALlMhhlzsECpQ8Z+2PAYuRCGOcBggVJPKKc8T2xkGKYGKUmgENE4IrqPiN4kotVEdDoRTSCiJUS0Vv6O146/koi6iGgNEZ2npZ9MRCvkvhuIrJ6TiFJE9GuZ/iIRzdLyLJLnWEtEi0q5jrHCWJm/MTaugmEOPErVUK4H8IgQ4hgA7wKwGsAVAB4TQhwJ4DH5N4hoNoCFAI4DcD6Am4goLsv5KYBLARwp/50v0y8BsFsIcQSAHwG4TpY1AcBVAE4FMB/AVbrgGusYv4dS0VqUj7EiGBnmQKNogUJE7QDOAHALAAghRoQQewAsAHC7POx2AB+T2wsA3COEGBZCrAfQBWA+ER0MoF0I8bywepI7XHlUWfcBOFtqL+cBWCKE6BVC7AawBDkhNPYZ499DGSOXwTAHHKVoKIcB6AFwKxEtI6JfEFELgMlCiG0AIH8nyeOnAdik5d8s06bJbXe6I48QIg1gL4CJPmUd2NRATxyFdiGyEVSEYZiKU4pASQA4CcBPhRAnAhiANG8ZII804ZNebB7nSYkuJaKlRLS0p6fHp3r1g9nkVX2JEoWWVAvXwTBM4ZQiUDYD2CyEeFH+fR8sAbNDmrEgf7u142do+acD2CrTp3ukO/IQUQJAB4Ben7LyEELcLISYJ4SY19nZWcRlMpVmrJjuGOZAo2iBIoTYDmATER0tk84GsArAgwBU1NUiAA/I7QcBLJSRW4fCcr6/JM1ifUR0mvSPXOzKo8r6JIDHpZ/lUQDnEtF46Yw/V6YdENTy4pBRVKEGLoNhmCIodemVrwC4i4gaAKwD8DlYQupeIroEwEYAnwIAIcRKIroXltBJA7hMCJGR5XwZwG0AmgA8LP8BlsP/TiLqgqWZLJRl9RLR1QBelsd9RwjRW+K1hEYIARnZXBVMfopa6Igj8aHUgmRkGKZgShIoQojlAOZ57DrbcPw1AK7xSF8KYI5H+hCkQPLYtxjA4kLqGxWZrEAiXj2BYmKsdMRj4yoY5sCDZ8oXQabKHXctz0OJxORVCxdSZ3z0xmfw1zc9W+1qMAc4vNpwEWSytdnjjZWOeKxoWpXk9c17q10FhmENpRiqLVBqefn6aMKGGcXGXYN4YwsLC6Y+YA2lCLJVnnhn7HBroCeOQqixgpLjjP96AkD1v8HDMGFgDaUIqu1DMVGbtSqcWtC0mMJ5Y8tebN7NHz47kGGBUgTVN3nVbocbicmrRi6vb2gUl9+zjL9dH5IP//gZvOe6J6pdDaaKsEApgmyVerygQOVa6YhLpVYu41cvbsQDy7fipie7ql0VhqkLWKAUQbU1FBNjxVRUKxpYKmG9HkOjvFrlgcS3fvcGZl3xULWrUZewQCmCaguUml56pc5MXvtHMujpG/bc15i0PtcznM547mfGJne+8A6A2hnY1BMsUIqgWiavoLPWQvOvtyivT/7sOZxyzZ8996WSrKGEpdqDrHIwnObnXigsUIogXW0NxdBpj5URVSVNdyu37jPuSyUsDWWkBjqWbI132Olqx9KXgf7hdLWrUHewQCmCar3cyinPJq/KEI9Zd7zaAwigelpxWGq8egUhH/uY1LrKDQuUIqjVeShjhVq5uwnZs2RKHH0/uaa7ZC2H21zlqXUhXouwQCmCao9cjItD1kD7j2ZxyBq4EESjoby6cTf+9taXce3Db5ZUl3JblIQQeGD5FgyNFheAUCOPLBLUpylYQSkcFihFUG1zcS2/vJF8DyWCekRJOlN8jfYOjgIA1u3sL6kO5R4tv7S+F5ffsxzffWhVUfnHSsi6Tq37rWoRFihF4Gd++PlT69DV3VfB2uQYKy91IX1nObUZVXQpGmnMNpuVVs9ym7z6hiwH9Pa9Q0Xlr5VBzorNezH3249iZ793KHgYgnyVjBkWKEVg6hxGM1lc88fV+PhNz5X1/OYor7KeNhSVNnmVcxCp7nMpEUzKwVuqhlHu0XJM9gTFVrMGmh4A4JZn1qFvKI2n3uopuSz2oRQOC5QiMDU0JWiqFW5YC82/0svX17qGEqeINJQyCxSS4/Jiz1Irfi81GTWKuUMsUAqHBUoRmF5u1QDL3Q7NYcNj4wUoyORVvmrY2k8p54jKwVspc3652lDvwAi27d1flrJ1csvlFL+6AdlaZRQ1OrDg76EUgcn8UKvRXxUlEg0lfCHllKEiggGCihQr1WRV9tGy7ESL1lAC9p909RIA5f+ui9JQ9g2NllzWWBmgVRLWUApA2cNNYaQc/RUNhWkoZTR5RVBGXL5hpQqE8pu8LIr2odRI25t1UAsA60uXpcIaSuGwQCmAmLKHG96eSi0/YR45Vf8NqPRaXhXRUEq4plybKa0u5dZQlGmu6LNUv+kBABqkBC/JTAllpqyRi6ojWKAUgOocjCavsr/01u+YX3qlRnqnKK6FKJrZ9tXWfoOomWcWYVksUAqHBUoBUMAaP+WfzVze8muFQq6znC99Lc36L/tgRf4WW99aaZu5wJgSKhQwcGPMsEApAFtDqbbJq8D0ShJNJ1yeYwsliqi9qOzw5Td5lZa/FtoegEgrwhpK4bBAKYDcKqTe+6ttlqiF9h/N0isFRHmVfDafsoXzt8hSoqhK2SOO7HkoRTvla6DxQRsERFJWBIUcYLBAKYAgp3ylVoQ1+lBqZ5xYEtVYesWrnCg7pVKrWfao4TGioURRD3UrWEMpHBYoBWBPeKrSPBTbKT/Wl16p8PkA73snIhjt1sIzKYR6H5REObm4VrSueoIFSgHEAxb6q/rExhpo/5V+kaO6Zq/RaM7kVfxJohrllvvR2qPyIs22tdD2gGjrwSavwmGBUgCxgHWZKiVQxrzJq5Bjy9hhRyEMKrUcT8mUbPKqjQuMoj0EWSIYMyULFCKKE9EyIvqD/HsCES0horXyd7x27JVE1EVEa4joPC39ZCJaIffdQDJ4n4hSRPRrmf4iEc3S8iyS51hLRItKvY6Q1wrA7Cupts21FjqtaCY2Vl5D8TZ5RVFwBGWgAj4Ue3HIYr3yEVamBKJYf81dFhOeKDSUywGs1v6+AsBjQogjATwm/wYRzQawEMBxAM4HcBMRxWWenwK4FMCR8t/5Mv0SALuFEEcA+BGA62RZEwBcBeBUAPMBXKULrnIR9K3pcmsoIsKXpZapxuKQniavSMqNoJAKUu/L1+fWX6vswIaxKEmgENF0AB8C8AsteQGA2+X27QA+pqXfI4QYFkKsB9AFYD4RHQygXQjxvLCe4B2uPKqs+wCcLbWX8wAsEUL0CiF2A1iCnBAqG8HzUMosUOyNGm7oUfhQCjm2hHsRlDcKc1VUpqBym5RU+cWepVaaZCQrRNtLr5RenwONUjWU/wbwfwHorrzJQohtACB/J8n0aQA2acdtlmnT5LY73ZFHCJEGsBfARJ+yykqQhlJLJq8Hlm+pTh307QrMui7ljuuP0fOcdudUilO+6KzOqpS7aQnXb8HZa6P3jbIW1X6f65GiBQoRfRhAtxDilbBZPNKET3qxeZwnJbqUiJYS0dKentK+4kYhnfKlxvQbCQhj1V/qy+9ZXqZKhKcSk+RKm8UuPLftsiPonqKcaBeGYoV4tkThWSt9r4jALpz7Hkr1LiqdyeJHS95CXwTL8FeSUjSUdwP4KBFtAHAPgLOI6JcAdkgzFuRvtzx+M4AZWv7pALbK9Oke6Y48RJQA0AGg16esPIQQNwsh5gkh5nV2dhZ3pZJYwFLktkAp6SzB1MvikEWbTwo6NpqQXu8oL7mvlPtaYQ2leB9Iaea9Gmh6AKJ9B6r5Pj20Yhuuf2wtvv/ImupVogiKFihCiCuFENOFELNgOdsfF0J8BsCDAFTU1SIAD8jtBwEslJFbh8Jyvr8kzWJ9RHSa9I9c7MqjyvqkPIcA8CiAc4lovHTGnyvTykoubNh7f05DKY9ICWrftfBS19Py9Q7h5zcPpfhTVHyUW+z5xoq/IKcRRqddVoNh+Qnj/SV8ebIalGMeyrUAziGitQDOkX9DCLESwL0AVgF4BMBlQgh1t74My7HfBeBtAA/L9FsATCSiLgBfg4wYE0L0ArgawMvy33dkWlkJcsqrcOKyayjGmfKlvwAj6Szue2VzVaNkKvXFRj2vV4eajeB5RhbWHPKeFCsYwqwKsHtgBAPDad/8Yc9TLqIoPbf0SgSFlUi5+5KoieQTwEKIJwE8Kbd3ATjbcNw1AK7xSF8KYI5H+hCATxnKWgxgcbF1LgaleKQNX0tSE6Fi5dJQAhp4FO3/xsfX4obHu9CUjONDxx9ccP5ITF4FOeWjMXl5FRNF51cvGkqYVQFOvHoJprQ34oVv5L/eYU+bFUC8jL1klPebnfKFwzPlCyD4i40yvUwvTJCdO4r2v3vQcgLuGhguuazinfLlPwfg9qHkFxTFCDUboqMOQ8V8KAHHbd83VNwJJBWbqxXBaao5D6VWouYKhQVKAagGNpz2tmtmK+WUL2JPWBJy+DiSLm5RJ+HYroDJq6gzWOh9m5/JqzTqTUMp7rzhNZRyC5TSyyfieSjFwgKlAFT7MtmRbR9KuTSUCjRw9U3u0SI/gq6/0JXRUIq/Kc66mjWUUp5nJTolve7VcsqH9/GU94ZE9bkAqyyWKIXCAqUQZPsaHPbWUHJhw+XVUcpp8mpIKIFSva+FFXIZpZm8/M8ZTWBCyUUElhMUXBCu/HAmrzB18KPcAjaK+10LTnkRwWCmGrBAKQDVvkyhfOWe2JgLYzVEeUVwjqTUUIo2eQnv7cLKqEyUVzZAm4pyteFy4jAzlqqhFB2ZF45y+1CiCBvO3QrWUAqFBUoBqAZmWrOrUhMbTUTR/u1vvkQxOi/ah1L+cwBugVIup3w0nZLfdTpn/Bd/htz/ReSusbBhNnlVBxYoBaCal+k7CaoBli1sOG/DtT8Sh6Qqq+Sial5DCQpxrpcOxWnyKk1DKccl68+z/FFeVvmlnKbUj41FSbnN51HDAqUA1Hth0lDKHjYc8LZH8aqW+l2Mys9DKZ4gk1eUo9xSy/LLH7QmWSHlFztp1m9vUDRdlISZTxOWehlQ1BIsUIrAuNpwpb7YaEqPwiGphGGx2gX0Trp0oRR8bCkmL33bw+QVwfOsdJ9U8jwUQ/4gzcI/aKB0gReWrK2h1N+zGwuwQCkA9dKlDbqweunKZfKy62Fo6VFMhspFuEThQyl/vpI0FK2T9DZ5lVB4hGUAQRpA6R12UD2Dv/Xj5+PJbZffKZ9/zoKpgdWG61WWsUApANW+TC9FOqIor6HRDPYMjhSeMYJWqIRhFHNIKuNDCT72O79fhcff3OGR17+cKMxVFYnyisCkJAKuNSiMvNwmubCoQVU1B0RRwmHDY5ggH0oUiwkCwKd+9jxO+M4S4/nL+U7mvgURQWFFm18KODbEwYufXY/P37Y0Lz3YhyIKrk9+IaVk1orxudCgaLVw5eeXpWNav87OH6JsoPyO7ijekSi19AMNFihFYP7AlvVb6vL1K7bsLSpfJE75gBWVC6lDJT7WFGiICd0Re+3PP65QolpO3bfD1rdL9KGY2m6QycvXh6LV0O9e3vnCO3j3tY8jXcKk2lyUVxRO+ZKLKJp6lWUsUApANdagTwCX7YONAXMFomiEsRIrX4ll7wtZ3mXEp3NyzpQvj8mrEv2C0C6xaB+KLMN0701+Q7sOYX0oPvX7/sNvYsue/egdKMLc6zpXNE55cxmrtu7D7H9/BNv3lrZY5liDBUoBqOZl9KFkgn0oXd196C5yxdYgdT5Kp3w1lx0J1jr0Y/2PHhoxd4TO6KP8/VFqKMUQVnA6NYAizxWwP8jk5UdYk9ywXJ3BbxAQRM50V3QROS3dp5A7X9iAwZEMHn+z23iMiWe7duJvfvZ8SZpYrRLJ91AOFIJ8KGFml3/gh08BADZc+6EQ5xMFmc8i0VBiah5KcQjDdkFlBGQMMlXpmFaGtsrRz5lfUBTmk1ImDIbtFIPCn8MQNIAoyeSlayg+fWg8RkCm+GV/AF2rLLfJq/j5Wl/99XJ09w1j18AIJrc3eh6TM0EWXHxVYQ2lANRDDpqHEpX90/0SC/vX+wSR+FDkb1DH1NXdh9Xb9uXXIaCTDkPQS6rvDaqnv8lLE0w++0tyJEekofg93aBVk8Ody39/0Gg6inkoCTmY8VvpWgiBe5duwuCI94rfufP47g5FGOFczO2mEGHJ7EM5AMhpKN4vlxIAUUWHmEZqRpNXFOcNGTb8gR8+hQ9e/7TvMeXSUAoJTVamGi/fUNi1vKKaPFkoYbNGMRM9SIgPB2gNfm0+7DyUmC1QzOd6/u1d+L/3vY6r/7DKtx6lvINhlh8qRXOwV6MIKYTrCRYoBWD7UEyfAA5w2hdKnoYS5KyO4JylL90d3hwVXIJpf/iClfCPe0iUIMGU65zM5W/Zsx9/8z/PY6/80mX+OUrRUPzrZ+8LGUXlR9DzDvJr+I+2w7UJpeH/5a0e4zH98ltEPX3eXxQN88zC4ndNpSwooZqiX95cH1JfNi8WKAWg2pfJV5KJ2ORV8NIf0SkoRRdWiMPcXEZ4e32Q8F651TLLea1eEGTysoMgfK7jpie68NL6Xjz42hbP/aVpKOEyRzHPIzBaLkBD8fOxhI3yamu0XLpLVuVPQg1LLnCl3D4U1wkLIIzTv4QYiKrCAqUg/DWQTMQmL/fLV8k2FsklFKuhFGDyCgqEuPye5QC8NRSnUz4/b5jRLgWMNnPzUArHKZzDHVe8huKfL0ig+HWOYVcbPnnWBADA3GkdxmMowCQbRWReqRMbh0YzuPmpt42mO9Vm/O4Fm7wOAAKjvLLBHVAh5L2kAaOvKMKGg0KTA/MbtgsrI8i0F66D0okHaCiei0OGiBgKYw+PAl9TUQFRb8byA/ZHpaH4TjSVBwavG+aHfAcjiMj1q6tqd/FYfhd60xNd+N4f38QDy7d65lXast91lnvNs3LBAqUAgnwoGXtEWx4NxV2PvPQITquKqOZqrQVpKGEFStzLh+LfEYdZaDDnwPU+qLR5KCGPC3G+Xf3DeOSN7UWfLGgtL18NJeQAQO3L+EiDIN+FylqaU16tFmE+JidQ8ve90zvoW4cwGop69zlseAwT9MXGbNQmL0PYsIlozFTFm2jcdSjX0iuODrQkDUUvM7+cMPNQgjq4Up5J2AmB+j0w1fWyX72KL/3yFezsNzmz/esSZFoMq6H4nScX1GI+JqiDDVqGvxD8nrvf/egbSss6eB+T01D8Jt2GqWHtwQKlAGwNJWC14ajUVVN7M4YNR3DOUjUU5/dQgo//4p1L8YNH13jWwYRet7DmEU8fSta/rvZot4TnWYmwYR1T2+veZwmSPUVGo5UicMIIPH2fn4aiMGuEwfUJW3YYDcVL+Kk2aWqbqiWG0dbqDRYoBRA0DyWnoUTjVMtzytvag8EUFkVkiz35ouSiQhXx6ModuPGJLme+QqK8Ql5zKU55//kIQU5iVUbhN9RhkvM5Tu+kTZMCW1JWBJUKu80vI3xdvDCZgfPPE6zJ+A0SAjUU20xZgkCxywru8P0+zGYUCvIa2IdygKOPXLwam965ea4NVWAjKdjkVVDp/mWU+ilZa7tYLSf8AWE7Mq+wYWeH7dExFNA5lXLvewdG8JMnuvLaR9gy9eqZfB1J6UMyzXgvRCv0wldDCRnlVUiUpNnEWLrJK4zp2m+gkAkQKKothonyqjfTFwuUAtCfrVdj0NM89xfYOtwNOmjEHIlTvkAFxf1COQVKsZUI2q3d54g0FK93O4ynmYkjAAAgAElEQVTpI9Apr+Ym+dTtqgdX4r8eXYPn1+1ynd97O+8cDg3FW2Co6zd1YoEmrwArlF6uu6wgTTB3nDQV+S694l+PUgdE+jn8nruqo997brrXqiX6XWfGR2DVMixQCkF7tl7qaqBAKVVDqUDbyqn74Y5334dovkMR0LnpgqAUH0qA0zvMMh4UMJM5TO2GR60FLPuGXP6NkLdSP8wkUIJCVYMeW6CG4vCTmPP6vQNhVpoImjwcxdIrYYSS3xyloPBnsk1ewZ9WqDfTFwuUAtAfrVdjcwiUgP1hcJ8jWHsovfEVusKuXzhpFB/p8txfhFM+cC0vjzy2OymUhuK9P8w3cpoa4gCA/aPOlZEdAQ4hF4ccMYx6lUAxPa/AuT8Bt9lPQ9H/9I2cCjEPJWiiqG2mLGEeShjfmZ95LkhDCTMPRffH1hMsUAogqCMLGokVbPJyvRRBdtWwxQshcM9LG7F/xLy0e9iaup3Azs4jZCHuc4c0a1jnCHeShMcEtCB/TyGjXZOzOEyMg6qb2wTiOK1PAQ4fimECorp8k5klSHhGpaH4jvpDzCEJ+oSI8hFFsYaan/abDtCUAB+Tl/Kh+Ji8otC0qkHRAoWIZhDRE0S0mohWEtHlMn0CES0horXyd7yW50oi6iKiNUR0npZ+MhGtkPtuIHnHiShFRL+W6S8S0SwtzyJ5jrVEtKjY6ygEgdwS216NQRcyXo2xYKe8e6QXon5heHJND664fwWue+TN/DLsjqU4DaWYWex5dShgtBz2HI3J/KYe/E35/OMKJsRo17R8S9jz6rcgyORl1FACThUYNqxrKK4r8RM2OsoE5OdbUMeY2qca4JQysg+z5I5tnvMciFi/Jg0kFirKy3meeqEUDSUN4OtCiGMBnAbgMiKaDeAKAI8JIY4E8Jj8G3LfQgDHATgfwE1EFJdl/RTApQCOlP/Ol+mXANgthDgCwI8AXCfLmgDgKgCnApgP4CpdcJULIYCknBob6EMJMHmFES5mH4p33rBtb5+01XtNcit0Ypj75Xc6YIsUKIEair+m6EUqEc9LC3LKR7FybZhIMdN6nCEVFMf9CHLKjxouJmdK8t9vIuMjnAuN8grjZzGhBI7fcZfc9rLvpxfCDCSe7dplPEbNozG942FmyofxJ/3+ta14ffMe4/5qULRAEUJsE0K8Krf7AKwGMA3AAgC3y8NuB/Axub0AwD1CiGEhxHoAXQDmE9HBANqFEM8Lqwe6w5VHlXUfgLOl9nIegCVCiF4hxG4AS5ATQmUjK4Qdfum39hPg3Zj0xhGmIyw4yqtAH4r3CrzwPYebPA3FYd4oqDqeZXjvz22HXZTQi7BreYUpy1gHtbZUEfnDZtFNo0E+FGPYsDyZV5vQ95vQ23OeQAn5zXvb9+BzTNCniNMhNJTH3uz2/DicIjeoCn4AXoeoOpgnNgbPlA/jx/nK3cvw0RufDaxjJYnEhyJNUScCeBHAZCHENsASOgAmycOmAdikZdss06bJbXe6I48QIg1gL4CJPmV51e1SIlpKREt7eszfWQhDyRpKSFuyV3lAcPRJeB+K9evVd4QZnemdeL7JSzuuyE44yIzlMHn5RuL4Hxe0Cm6QBuOui18d/JzEOZOX+3mbO2nTcWaTl/Vr6pCDfBNBz8RvNnzYdh/UEQflB3LXH8VH0UzV0D8r7RfNaZrxH/QsrHMXPxCpJiULFCJqBfAbAP8khDCLfe9AF+GTXmweZ6IQNwsh5gkh5nV2dvpULxgBgUTc7EMpJGw4lIaSZ/JSDdVUv3D4RgyFKCXtECjuOua2ixYoAdnC+mmCNZDcttfzCPvpWj8KyZdXhwIHCIDZKa8cwaOmVR4CtTn/OqR9BIrT5OV3Dv+O2H0ev/2ldMRBnblaq8t0TG4JJu/y1aDU7yuYKu8BFTZMRElYwuQuIcT9MnmHNGNB/nbL9M0AZmjZpwPYKtOne6Q78hBRAkAHgF6fsspKVtNQgn0kHvmz3seayF96RdXDpKGEa3zqMO/Z4/LXr15+GorwvwdhCPIvOTSP0ALFf79XJ+bUUIp8sUNofLYJxC2cHds+gwBdoJh8KEqgGDqxYIGia0sewtnhlDfv89VQVNiw38g9ILpKmfRK6Ydz75n3/n6HQPGog2xLJsHYkLD6kMER72VwrDoooRZU29qilCgvAnALgNVCiB9qux4EsEhuLwLwgJa+UEZuHQrL+f6SNIv1EdFpssyLXXlUWZ8E8Lj0szwK4FwiGi+d8efKtLKSFQINSqB4NBb9PQgyeYUSKAaTV6mjFpXdS83LvbDmc+ij3FJMXsUu+R72g036I/L7Zjzg3YkV8p0RUz0KMV3kC+fALHllDwc45Y0rZQf55xzCNX+/w4fiqkLYYJQwS54EaSgjmfD32wsRIDgBl4biZYlQs+gN+VNSoAyNmsP2o/5YX6VIlJD33QA+C2AFES2Xad8AcC2Ae4noEgAbAXwKAIQQK4noXgCrYEWIXSaEUHf0ywBuA9AE4GH5D7AE1p1E1AVLM1koy+oloqsBvCyP+44QoreEawlECAEhYJu8vH0oWc9t7/3BDcVkOihlJVVAe1E8J/upY8z5dXNfKSYv45cvC3DKl2LyCppX5Jys51sl87drAka7gD5z2vt5B53fIcQDIotMi0cGham672Xc1Xick0SdZQStcacIWqXXqqfwPIdCrTpQtEIZov32DY/6HhO06rgS7vuGzBpKmOjAWqRogSKEeAbmCcBnG/JcA+Aaj/SlAOZ4pA9BCiSPfYsBLA5b31JRz9U0CQ1w2ky9BoqZAk1e+RMbVbp/5xWEOsxr2ZAwQkt/4d1RQ4VEeZnOEWzy0jQUX6e8PjL23x80UdVoZgwIs84J6OCHk3cvDdv559A1Nu9jwkZ5mc7jNjMmXVHY+vvgvpX6kvl+zyuKsGG12kCYjjibFYjFzILRVA2/awWCTXcqz9vd/ea6BZi8SvmkQjnhmfIhUY8vKdVVTydvAU75MFqGafl64wJ/BS7+5B3lZe18ck2PsRPUzz/i6qD0qhW74GBQxJFDQ/GNlNG3vZ6Xfs78kzo7juK0rTBzWUzaQ6E+Mf18blS/aZqHEmRiCdKWnPNQnAf8/V2vGvd51cFPoATNUB8qQKB4BSiE8Zs5J2p6aCiyAZvecZV/wMeHYrcbw70o7TPJ5YMFSkjUA27wM3lpPpageQ1hll3Pj/LKL8drfxBK8Hj6ULQyXt6w2zO/Hj+f70j2F6qmcnT8OierjiE1lACTlUNDCfChBIUNB4bV+ml88tx+PhS/jjiMT0ldnslpH2TqDIp40wWy360oeWJjwH22NZQQASFezzxMqHbQihhB16Hun9/SR0EC3m8OSzVhgRIS9WCVycukgagIDq/9ekMsSkPRzuNFeCeu9es9sTFXiMlpmHb4UNx2ufzzGOtheCecUUH5+4Xh2LxyAgRPIWHepg49aDHNMBPUwiyK6Ieey2hGlOkmk1eunqZRtf859P1+mrKvsBHB98FPQxFC2KG4YTQUT4ESSkPxF55BviBV7qDfWnoBAxWTL6zalOKUP6CwfShqHoqhA1Iz6YM7sOARRv7SK3LkYxpFyt/ZB7eH8tF4mbz0bEo4unHMQ8lzJOfX14TRJOAacbuXng+72nCQ+UI313mZgsKsQRW0smyQVqmf21dDMeZ2CWBTJ2afw7+TC7pOID+KC3C2Z7/H7jcAsKOjiowES2etwJkYhRuweZu8gp+535wbfb/pWtX1+QmUoGX6TQODasMaSkjUg82FDZsEijR5BTh5w7SH/CgvGMu26milB30mNetznN5ZK+HoxjEPJe12JAe/kIqgTynrdXXs15L8Og5nCGj+/hGt7hmPBxJmhreak2B0nobSULwXRQzrEwsTjp7TggxmxoDoKUeEYpCG4idQ/HwoAcIZ8I+gUs+zMRmHEMEDGi8NJUhYuM/tPkYIEahxqlvpNw8lKNycfSh1jm3y8g0bzgmU4JnyYTQU59/qZTc75S1iRP5mBzu/v8nLPELzmYcSsrO36uGdHhSuG9YpHyQQ9JnK3s9Tq6tH/p39w/jd8q3G/Va6Ob9CaQ3uNuF4zr6jfi1PoMkroJ6GZxLU0foJnA8ff3Bg/fRz+PpZfMxi6nk2yRC0oD7Xy5+ktyeTQPIL1ggyo+p5dg+OOgY1XseY7pffd4iqCQuUEGzqHcSfV+8AAF8NJCOEPWmpmJn0bkxO+aCOOkhDUS9jkMnL1GgdGoq7jo6y/Otp7vz8jwkfNuxfF4eG4ilQ/M04m3fv980vcxrP787rNkeFncCqCyKjyUsmG01eQUEFQY5oQ8DJXS++gz+8vi3XyYdwyvsNtvzMYrqGAgS3Py+hlA7xjjqFhjl/mHlWyzd5rxacG4gE16GWYB9KCC64/mn0DVvqqd/ikFmHySu/nII1lDx1OnceT2Qywd/soEZY3lFeIu84N3qn5DZ5BS3R4TjW1UmpOQHBa3D5d27e5efv130oQRMbveqhC1zjaNSeMGispl2O2y6u32d/R7cuYL2PyQkt/xFxGBOL17U4Q8lzvoFbn90AQJ8f4l0/3VRUrIaSEyg5K4F7voyOlx8i6Jm7z+1u42EESjYr0JSMY/9oBvv2j3oeY0d5hXDKCyHstdqqDWsoIVDCBNA+sGUYpSlHtpfAKH3pFX812A4HJvK1vvtpKHrRQfZ2z2N0rSBAZpqi3oKW6nDW0Xyl6rh4jDyF20g66+h83Dg1nPzydQ0naOKj3wOxO3vXScLO1HeY5gLMLEa/laqm4TxBfi29rropsbkhbjzO6/x+x+j18BIGSpCpTyqbIsEUXtpaOiCCy6qfdUxDPJZv8tLK9LvWye0pAOa5KMECvrBVNyoFC5QCSfh8DyWTyQmUYiY+5h1foIai9geavOTL6DdTHjCbR5w+FG+h5y7LC5O9OWipDr3YME7eRIw8hfBwOmObYrw6qCBTkq7hlOJDsZdtd2soIecahLHbq/Ste4YM+/3DbYNG3vp+XdA2aSoCkV9YspVuddLmNu4Xkuv2oQSZnT0HfT7ah7sOyTjltc/REB19JivQ2mgZh0yRXrYPxXQfMvp1sECpW0zzUNKZLPqG03YUmFdnHObF13G/VEGOOpUaI+8RuV1XXx+K1lCN6z7pQsfslA+2YXubnIIErzNazmc0K3KdlMnm3pCIIREjgwkzN8Pc61Kc9fSvg78z2jvKK6yGos7hFy6rjtmxz1ugKO3IdDuD7nnWIFB0DSVG3oJdL9PPB6nXw8uZHcaH4vfpBXdaUJRXMhHLe0edHb3ZvNjSYAkU0+RGldX03PX3zm8Z/ErDAqVATBrKz59eDwB4fp31aVDPCBLDKFzHMevZdYj609R5qfLjMX+Tl98sXP39ML0QaR+BEmQm0hlNa9dqMBV4CUY9xd/kJQVKIma0uTckYkjEyTDqziLhs/KBfu1BWqMIMep2m7z08v2W6bA7uHh+B5ern/Xbb1iQMKcdBY+I3cvtAC4NJaMLlJybNk75I3pFRntWQPDI3KsOeQLFS4sZ1Vd5CPKheNdVXauXycthivIRzi0pKVAMk4fV/QijMf7PX972PlEVYIFSIKaw4K179jv+9hQo+ug/wDQBmKO8guahxANsXkot9xJMjkmDJpNXxnxMIR+lcpgHDJFbQVFe/lFD1m9DIuZ5LSOZLFKJOBIxb4GTycJ3KR19ZGgaIDgngRoEtMHkpdf58nuWw4RtLkrEAieL9o+kPe9Z0KdzHU53jxGxfn/0/S0pp8nLONkvE06gqOvwer/yTF5eAkULGPCO8go2YzoEuOuQtGNgZPZBNiXjiJFZQ8lNYg5+B296kgVK3RI3OOXdL1mQySvMom+mNa2CnOXxGPk7ge3QS4+OQUsynyeX7h4p6qcN1lB0k5f3aN/bh2LW4nSyDg3Fe0TbEI8hHvPWULJCaBppfvn6aDdohjoQ7JPyM3n5oXdwQT4UIfw1jCC/AeBtYnFoMA6Tl6ahxCjQyZxKmKMo9fN4fSjMHeUVFEjhJZSc0VOeVbDr4CXA/da508uNxQjNDQmzhqI9L886aOdpafAJZaswLFAKJGkYsQ7J70yrxuzVWEcCJtK50002dLPZwPp1L8ltOoeXhpIVwo5kM3eAeidp9qEEhQ2PGsxcDrOZx8Xqxe4eGDGWr55RMu6toQw7fChenWTW9pl5XctwCKe8s4PxNyEOjjrNUWEnr2Xs6zR32HpZ3gJFOeW9z6Gv6+aloahRt7t8dd8+feohlg/FZK4Na/KS9fRqm+q8TT4+FMdkVo8ydI3Bb0BFZAV7uNuFqheRj4YjBGJkmebMTnl1Ln/BOq45iSMmtXoeUw1YoIRA75/bZXSG+8VQL1xbYxKAd2MdCWFz130JplBN87pRyuTlv/ZTrvPwGpXnXmpTB+iIlPExeQWNsEcNo7mgGe6q2I6mJN7pHTCWr7I2xGNynad8rbIhYWko3qsN55af8boS9TGncc1Jn/kf2vkC7ufugVHP9CBUO/HTUHQhoGtWCtURmgYB+kja6zoyQtgOePcKBBNaGnDNx+ci5tfJamY7/e/86xB2HbyeJwA0NpgFyohBK1bsl0K9vTGBIY/7pK4pESMpIF0airyPKYPfDrCuLU6E5oa4cQHWoLBhNUBoaUiwU77e0FflVTHueqd35f2v49GV1kz6aeOaABSvoegvq8nnEmTWiMdiob47YfrqpHJqmuqod0jRmby8BanXtaprm9yewt5B74lhet6koZOyfCjmKK9MVviGiasXuTWV8FkIMNgEotJ3DQw70k3fLsnLrzuJDX2LMyoovxNTgwdTsxkcydiDKZOGotrNiEsLUEI55mPyykV5mf0fgH/I+kgYH0ranB8A9o9Y+ye0NBg7+0zWWrA0mcgfiORMd3GfgZ81V6wpGTeu5xW0Bpx65i2puHH5lmrAAiUEemczf9YEAM5O7+6XNtnbN150IgCDSh5iItzuwZwJx2l/D+EsVBpKwFP1s0OPZgUaE+bQZyAXcTSuOemhoeS2A53yBpNXUKSNShvX3IC+YW8nM5B7uZuS3nb5DTst7SZhGNnra7N5zrTXOjCzEzm3bTJhqXoOjWYdHUyYFan1chuTcZ+1n3Lf6vHqgN7use6F6YkNjmQwvqXBmD+TzWko7oGTMhtaUV4BGoq9EoV/AIN1Tc5jlAlS1cPTShCgoQzabbvBOPJX15SMx/IGVKqNNSb9/VnxmKVJ7TdoQapYs2C10ptZQ6lfPnz8wehss2a4mh70lPZGAN6dx38+/Ka9bRqtbtw1aG87NJQQM3DV4UFhw6oxepouMrmRpqlDG5ArB4xrSub7UBwfKAoSKJo2ZjDveYYNy7RxTUkI4VzJwFm+ernzNa50Jovdg6N4eu1Os4YiBJIx/yivhngMibifecN/RGzVJdfZ7+ofCTzeTU6g5Iex6nVVk+ncHdDWPfuxVy4BYso/NJrBuKakzO+h4WSFrb07l7TJ2loe+flQXCYvkyx1fHLAVZgSFh2ynkMe9dTrroetK5RWEkpDicfy6qD+9tNQssLK35yMY79JQ8kGmLw07ZgFSp0ye2q753pTOvEYgSjf/+AOKzaNJP/ujqW5Y7QiQmkoshF6fThLR5XlNdJMZwVSsgM2aygZJONWlIrfgoaBJi9D5xA0X0eljGu2Oo6g9ZBsE4hWV90nYEV5uYMLrLWlkgnpQ/G4luF0BqlEDM0Ncdv2bqoD4KehCHugsm1vbuJhWB+KLjjNaz9l0ZryFih92twUU9saHEmjo9msoaj1qdz7lb8BsDTnoG+E+C1dBDjfK/eASAkLJVC8QnIdUV6eGkrOL2YSKCOZLJJxQkM8lncv0vazMA8yssJae6upIW6M8hr18XMCOWHZ0ZTEiIfgrBYsUELwvqM6AQBfPONwe46H6WUnIiRjMYy4OtpeVzRSmAge/eVR9vRUIjg01Fq7ylyucjZ6CxTrZYkbIp8AS0NpSSWQjFNeRxpGk1IYTV662cyjCurejZcdXJ9hsp5uCgKcHYh6kb9y1hGIxyhPMKrqJHw0lHtf3oQ+eS/6h/3DPwE/H0oWx0xpAwCs2b7Pka44enKbZ14g9xxTPvNQdIGSH+KuR2V5n2NwJKehmMKOGxIxxMj9nZmc2dBvpnzaJVBMbUd/Tu7rUH+3hxQoXs9DCZTxzQ1Gp/yo1EyT8fx2o94ZP+GunPJNybhxHoryU5peIZWvoznJGkq98b2/nou//Mv7EY+Rrb77CQSroTn3KzPRf/+fEwAA3/ztG3n5lE3/4ydOs1RubeShOpeOpiQGhjNGUxBR8PdQ1Mhr2OMalDqfMEQ+WdeSQUtDAq2NCWzYOejYpwuhQkxeJqe85/IZsl6q4zA5NtMuk5djAUP5ws6c2OI5Uz43t8OskQ7Il7o1FbefrxvHhEAfDWXGhGYAcAgm1Vm954iDfNf1Gs2oQUDMaFIacZi8nJ2YEq5zp3UYBcr+kQzGS43QK0pMmXFSiXieyUvN3bImmPpHugXNQ9GftZdAScQot6yJx+hfT/N6h4dGc1rnUNr7PRvJWNGBDQkvk1fuOvymBsRjUkMxCZSATxmrd7ijKekZ8VYtWKCEYNq4Jsyc2ALAso02xGN2Z6Jzw4WWQz4Rz39xlCO7rdH8xYD3/+BJAEBnWwqpRMzx4qqG29mWwkgm6xm/bsW3W8LA78NTqjF6aigZy29g2YdNAiWN5oY4Du9sRd+Q09zkXAvJWAXHNQFm05DfkifK5NVv6Mxzo8X8uUHqHjQmY0jE8juGsFFHgGXHNi1pMprJfSPHz+TVkoqDCC6nvHXOllTcsxPPnSOLZDxmNCllswLprECbMnm5ylJtqbkh7jkQGUlnkc4KHNSaAhHynjlgtZtEjNCQiOWbvKSGkkrEjKPpvLW8DPd7/0jGvg63BqHCwJUvx0ug6G3Fq8MfHMmguSFuf/XRtMRLMu7tlNfbjUlDGUlnkUrG7CXsvVAh6cb7MCrNzrKetfKNeRYoRdDamOtA9JHBR+SX6ayG5nzAao7BVBlW7MdZx0xCYzKOobTeAVrbxx7cDgBYsWVvXr5M1oqkSSbyz68QQtgdiGm9sXjM3+TV0z+Mg1pTaEnlx+qbnOte6C+Bfq6RTNZeuPI3r2zGrCsewq7+XEitMv8pk5dpcpgq3yuMVL3IjYk4WlOJPA1DjeKVmcjvhd0/msX2fUPo9lh4sX84jQkyOsrbxJKW9zyG5mQcA5qGojq/jiZ/s4bqBOOGVZWVdqM0FHcnqEbJLamE5yBARR42pxLoaEpij4fPKisHM43JmMP3oAQNYHWyputQeWyznEH4Do5mMKHVup9up7sKA1fP22v0r5tHvbSl/aPWCtRKsHmZvUaVhuLhQ9EDJLwElhACw+ksUnFL8JnarrpPZo3RCu9PycGS6X5VGhYoRZBKxPBs104AwB5tHoT6yE1bYwL7XKM45UOZOq4JX3zfYXmzbE+6eomjfEtDyTW2D//4GQDAu4+YCABY1+Oc0Ne9bwhPrulGLAYkY/kmN8Xl9yzHqm2Wnd5LQxmVUTle9mHFjn1DmNyeQlPSMm84fD26uSNg1KSPxt3hoJOkk/r2598BALzTO6gdKzWUprAaSn6QgeoomhriaEnF88pQ+1tsgWJ+Yf+8ypqD9IEf/sWRvn8kg/U7B2yB4lWGMn02xAkHj2vCpt256+zuG8bElga0pBKekVWAJfjuenEjdvaPIEbkvcy+fM5tRpNXTnvOZEVeGf8k1xEbGs1gz+Ao7pDPREcJtbbGpKPTHs3kTF6WhuJ9HftknskySnLIo6PtHRjBnsFRW2C4j+kfTqMpGbc1FC+net9Q2jZjes9DyaCpIVfGsEcZaoUFryiv3MRG7ygv1fGnknE7QssrWk0dN2wwuynB5xUKfv+rm/HM2p15eSoBC5Qi2LZ3COt2DmDrnv34+7teBQB85rRD7P0zJjRj1dZ9jjz3LrXmqrQ3JtDZmkI6K+xw16fe6nE47a2RR05DWbujz943S5re3ALrC3csxZvb+zA0mpVLjXh3gA++ttXe9hIog8NWQ03EYp5hw7sHRrB1z35M7mjMjQS1l05/ibxMIzpvbs9dl3vxwUltjY5jG7TJNW4fisl/4RYK+jl0k1drKunQDIBcp6s6YfdL/+8P5HxgS752BgDg5JnjHcd8+8GVAHKmv55+58RFAPjtsi32MRNaGhz3rKdvCJ1tKbRJDcqrg+rq7re3k/GYp6/lT3LSbYvB5KUm801ssYS42wyjVtD2Y3AkjeZUAu2NCezUrnP/aMaeF5JKxow+A3Xdk6RA8TIF/c9T1iKIqt24j3lzWx+OmNxmazn7PMyQfUOjaGtMGjXwt3v6MX18MxoTSijlH7Nh1wAmtzUimTA75VNJ7wAJpXmkEjG7nu72u0VGhI5vTmI0Izy1uv0jaTQ1xO2ITNVe05ksvnbva/jMLS/m5akELFBKYFf/iP2ynX3sZDt93szxWL9zwO60Vm7di7XyxSci21Sj1qFym6/cGor+4sye2o54jPI66zWa0EkmzP4Pne0eJprewRFMaGlAg8He/fib3cgK4JxjJ9tLXOj1y024ittzG7zYumc/lsiRPeA0W41mhP1FO4Vj/kDW6UPxMhts3DWIK+9fYR0nBY8uFFSdU4k4Wj00FH0WvDsvAHuU/pWzjsD08c04anKr7W9RvLHVeq4HtTaguSGON7f1wURWWHNR9OfW3TeMzrYUOttSyIr8mfQAsG9/rt4tqbinL+fr//saAKCz1bqneSYveS8mtnqbEJWZVUU7Avm2fStQI44jJrVi9bbcYKpvKG0vR3TIhBas7e73FIwq9FvN4/J6pmpQcd0n5jrqrXhn1wAO72xBYzKO8c3JvFD9XH0SxqCTLXv2Y9bEZlur9ZrLsmdgFFPHNaEhHs+bHGwPYhosbc+tXQzZ7S5mmyDdUYpPr+0BAHsiqXvwmEfHi/IAABmVSURBVMkK/G75VntxUyA3ONzZb17brhKwQCmCww6ytITuvlyHrIQEAMycaEXsbJJmmr/52fOO/MoEsksKlP96dI1jf2My7vChqAb3y0tORSoRR1tjwhgqC8B2FprWOtLRJ1L+5pXN6OkbxrjmBkxsbXBMsgOAN7fvszunudM7NNOD5v+Q55jY2uBpawcsQfrDJW8507QVAkbSWbQ3JR1ayX6Pc7Q0WKHLXiav7z60yt5WgkcXCsrc1twQl2G/aft+7eofxsMrtgPIaSgmn5Sa89PckMj7Zok675T2Rsyc2IINu8zrjrU1Jh3RgUOjGby+eS/27R+156j09OULlIdWbHWU0TeUdjx33exz4iHjAHhpKFa9J7YogeK8jrnT2jG5PYU50zrwzQuOzTtmaDSD/aMZjGtuwMEdTdg3lEY6k8VTb/Vg/c4BWyjPndaBvqE0NvY6IwOBnDahTJ1uYdE3NIofP94FADh55gRZb2cU5MBIBuOarGsY3+Ld/lZu3Ysp7Y2eQSeDI2n0DaUxqb3RDuRwa1S7B0bQJ4NSkgnKE85KELY15mvFQE5wtjcl7eACk8n2ovmHyGt37n91424AlqVE+VDUAEgfJAoh8GzXTnz2lhdDLzRaKixQiuD6hVY01/JNe+w09TICwCEyBPQd2VmriLArPngMgJxA+csa64Vz49ZQlDai8gUJlAZpI3Y31C0eI7a13dJ8MJKxhcX45iQ6W1N5HdjHf/KcVse4vWjmam3uxMBwGi0NcYxvbjBqKDc92YX7XtkMAPj+J44H4JynMzSaQWPSGX6qj86V36qtMYEWD4c6kDPvALmRnsq3fNMe3PrsBgDAxNYUWqXvQL2U//Tr5fjRny2Bp0aRuglR76SVkNgzOIKn1+60n9W3H1yJZ7ss7fXfPzIbhx3U4vmsGxIxEAELT5lhDQRkHVTbmTGhGVM6rEAOt138W797A798YSMA4I//+F60NSaQzgp7lPzkmm4c861H7OOPn24JFHe72D+aQSJGmNJhaQfdrufeP5xGu9Qy1H3VNYjufdbxk9sb7TK6evpx8eKXAOTeDTUB72v35n/bpW9oVPpgvL9kqC+cac+E157DXrujtvK7fTmANcB7u2cA75oxDok45YW4Py3v7+T2Rnsw4BZKysT9Tu+g5ZR3Ddz2j6RBlKujW8tX9exoStpty/08lJahIkvdE3fV351tqTwN5Z/lOwwAX/7lq/j0L17E02t3+mrHUcICpQjmTGtHKhGzR0wAcHBHzuZvC5TeQXskt/CUGfjS+w4HkBMM1z+2Fmf+4El7VKZoTMbR0hDHzv5hpDNZe/SmXrbWVNJhpxZC2Hb6c2ZPxsGyA5r77T/Zx6zd0YczZVjyb//+r/BnafdXL53eYTenEuhsS+XZ/I+f3gEA9nwE1ZE/sHyLXY9fPLMeAyMZTGxpwPa93p+b1WeD/80pM9DWmLAFygvrdmHXwIg9qlXonfGewRF0NCWRiMcQJ8JvpHBSvLBul+2bmNjSYJt67npxI4ZGM/jYT57F65v3IhkntDcm7HOpF1sXUErz1F96XZt6YLmlIajovWUb9+AXT6/Dbc9tsI9pa0zi2IPbsGHXgENI/2nldoyks/jX84+xJoomcmGoSvu9+PRZmDPVMjnpS/cIIXDnCznn+Oyp7bZpqW941FE3wBLcyXgM08c3YZ12L3+45C385Im30ZSM23NhNmuBAau27sMfV2y3Q3/VB7P0+7FD1nVye8r2I63Z3me31zOPmQQg10G+s8upoSzbuBs/f3o9WlIJz5Df0UwWZ/zXEwCsZeHdx+wdHMXJ3/0zgJzJrL0xkWcWfu/3rTJOnDEOE1oaHIEtz3btwhfvfMW+DuVP0qMLn+vaaZu4GxMxdMilf3RfjWq7SvC6tVZdoCgh7f4Ew6vv7Mb45iQmtHgLtTe2WAO4P3/1fbYfUbVJ3af2yMrt9vbyzXtQCepaoBDR+US0hoi6iOiKCp4Xl55xmP33tz8y237hgJzAuPoPq+yX591HHGTv14UPYAmg0w6bYP+dSsRw1rGTsbN/BC9v2G13+uoFXb1tH55euxNH/dvD6O4bQv9wGiPpLC4783D85KKTcIg0uenc9eJGe3vutA67o7zhsbUAnLb4D86Zgs62FHoHRhxmsn1DaUztaMST/3wmAGDBidMAACfOsDqR/9U69uOnj8Ob2/scWlFXdz+O//aj+MPr2zCpLYU13z0fgCXUbntuA4QQWHjzCwCAt7v78c0LjsXUDsv8oHdCT77VY9fr6CltGBjJ4OdPrQMARxkA8Mq3zsFBUqD8efUO21EOWOYqotxEOOV/mC07cCAXdfTQ69vw86fWQQjhGC3/5KKTAAA3ffokEAE3PtGF7z60Ou/+nzxzAoQATrnmz1i6oRff/O0KXCo7sOOnWYK6QYsaeqZLjZZTjrZ167PrAeTMpTrKhKKepa71fUiGtB8xqRXrenKdjnr+ne0pe6Xszb25Z3bBDU8DALqkJqs+mDWoBTGskU7yye2N9mBqXc+A3W6VkDnzmEmIxwgnzhjnqPeFP7eeV0/fsO0MH9QE1nd+nzNf3v65+WhpsMJ6lXB+9u2c5nawvIa2xgSWbdyDWVc8hKHRDN7Q/JTHTe3A+4+ahJc39NoRber6AGtwoEKTdbPvRb/IObq/du5RdrtS9Xhjy17c9eJGnDBjnC149WCPZ9buxOdvs5ZW6mhK2oMQt6/n7Z5+zJzYYgv4jS4BvHzTbhw1uRUdzUkc3tkKIuCVd3b7Tm781u/eMA7woqRuBQoRxQH8BMAHAcwGcCERza7U+b9+7tH2topM0epmb3/0RivcVzUOAI4OAgCWvrMbh3W22kIlEY/h3Ydb4cEX/vwFXP0H64VSI5r3HmkJp5F0FvOvecweoR81uQ0NiRjmzRxvN/an1/agq7vf7lw+OGcKEvGYLfTW7RzAwyu22aPM+750Og5qTWGq1HK+90erc+wfTmP1tn346AnT0CE1FNUR379sC55ZuxMvrusFAHxnwXE4Z7YVpLBk5XYMpzMYGs3gtU177NHcNy44Ns+JvVpTy888ZhK+cMZheO7Ks3HSIePxmhxh/fsDb+CdXYP26PSfz7OewzWynpt355v11IgWAO55Obcy9AVzrU5WCerewRGs3LrXNiMBuWCAF9f34po/rsYP/rQGe+Ro8NIzDrM76nHNDTh6chteWt/rOLe6D5O0IIMv3LHUIeCPm2oJlMZkHJt6rWAF1YlMc81b+o/fr8LCm593BHIov8a08dax979qtYc9gyOIEbDo9Jn2iLmlIYGVW/ehfzjtMCt94qTpaEzG0dmWwqsbd+PN7c4oxcV/ewqA3PeANuyygk6WbdyNf/udFfE2ua3RdmZfLwWVmwXvmoon3+rB/pEMnl7bg76hUYc2GosRDj2oBa9IPwEAhyZ2xlGdICJMG9eErXuGkMkKh4Y6Sw6m2lJJO23rnv122D0AHDKxGYce1IzhdNYOnPi2FFqnHTYBh05sQVsqgVQiZg+Ific1XsAaRBzc0ZTn27pQDmSmj2/ClHbrWdyphVj/3R0v29tTOhoxoaUBRECvNvXgtmfXY9veIVwwdwo6W1NoTSUc2vmfVm7HE2t6MH28dZ2dbSkcN7UdT67pwX/Ia/ji+6zB7skzx+Olb55tt0G/1TOiwjxtu/aZD6BLCLEOAIjoHgALAKzyzRUh/3Le0fivR9fgUOmk17nxohPxD79aZjv+3Gsx/frS0/B/tJH08dLhqezXE1ut6B7dRKIWprxl0Sn46I3P2OGT6mVQo2kiwt1fOBXn/OgpfPYWy449Z1o75s0cj5s+fZJ9zNULjsO3HliJL0u7MABbu/noCVPxs6fexm3PbXCYb/RIHzW/YPW2fXaY4rRxTbj49Fl2R/zt36+y66f47sfm4CPvmmr/ffvn52PR4pdw0S+s+/G9j8/FRafmwrDPOmYSvvvQarz72sftF/z6hdYSNnPl6B4AZl3xkOM8V39sDkxc8/E5ttPzsE7ri3d/fdNzeccdM6Xd8fdPnngbP3nCCl/91MnTHfvcWsOF82fgmx+yxjidmllzt9aBPPHP77cF9Mkzx+PulzbiC3KB0BMPGWcPPp674iz81bWPAwBeWNeLF6TwfupfzrSf2WwZjXXTk2/b3xk//7gp+I8FufvwttRO5lz1KP71/GPs9AUnWM+joymJJ9b04Ik1PXak3T+efSTee6T13E+Qjv2v3L3Mca0zJzbb13HqoRPwohSsd/3dqY7jPjB7Mu5ftgXH/vsjcHPVR6x79aG5B+PGJ7rynqd65oA1SHhoxTY8tGKbnbbuexfY78jRU3Lv21n/Lzc/6OHL3wsAOESa3z5647OOc9xz6en29nA6i9ue24D2piTukkLtT189A0fJd1lpGErDUvz9+4+wfUmLn12PtsYEHntzh+3bOmHGOFvTE8LSEpdt3I2jJ7fhF89YGuhZx0wCkSVcb3tuA06aOR4tDXFbqz3/uCn2+ebPmojFz663fboffddUfPUDRyGViIGI8INPvguDo2nbFF5O6lZDATANwCbt780yrWJcduYRePVb59hhlTofPn4qTpIv3zmzJztGyQBw6mETccuieQCsBvDXJ01HSyrhEE6Pf/199rY6FrAcuY/80xm470u5xp+ME+Zonav7s6BvbNmHebMmOLSnz5w2E1cvOM7+u6MpafsbGpNx/NuHjs27rlNmOedaXHbm4Y6/F54yA4DltNeFgs5nTptpCyMA+KvDJ+Lwzhbbaa40MMW5s62XRwmTC+fPwIITpsnrjuV1WgDw2lXn4rOnzbT//t7H5+KS9xxq//2Rd02174XXgOAnF52EDdd+CBNaGvDm1efj/Ud35h2ja50AcNn7nffiWx+ebY++2xuT+LTrfvzDmUc4zv1hqe0o9FUVpo5rws8+cxIOam1wHKObN1tS+eND99yYf/1gTohc94jlk3n6/55pj3j1AcMO6Ww/vDNXx1QijhNcJisAePKf329v/+oLp+HC+YfgGxcc4zD1ApZWeOH8GXn5X/rG2fjcu63n83fvPTRv/w0Xnmg/cwB4j6uNfPF9hzk+fX3GUc79AHDr355iv6vHa++K4peXONuRar83PLYWuwZGcPLM8bYwASxt6NzZkx15/vdLp2PGhGYk4zH7eV//2Frb7/GVs47A3V84Le/cT6/diV88sx4HtaZw2+dOwRGTrPOcd5xV/j/evQyX3G4NNH766ZPwN6fk7uHHT8zdl8vOPBzHTe1AYzJut++O5mRFhAkAUK0sKlYoRPQpAOcJIf5O/v1ZAPOFEF9xHXcpgEsB4JBDDjn5nXfyZ/mWi/0jGdz36macd9zkvIl6gGXvtyZ+Fa8obt87hA27BnDaYRPz9q3cuhcvre/Fso17kBUC3/vrubbZTGfHviH8dtkWLDhhal7DS2eySEjbfpzI83v12azAii17MWdah0NQAJYT8j//uBqDIxl8at50nHroRHtFWXcZdzy/AcfPGIeTDhmft39n/zAeeWM71u7ow+UfOMo22Sm6uvuwbe8QXli3C5e857C8/Yon1nTjyEmtdgeqUAsLrty6D5PbU3lmTMXXfr0c63cN4MaLTsozRwFW6Om//fYNnHvcFJw/Z0re/qHRDP7yVg929g/bZiadvftH8c6uASzbuAcfnDvFs90AVgj35LZGO4JNv479oxnc/+pmHDKh2R7puvnBo2tw4xNd+P4njnd0TkCuPdz67Hp87MRp+Po5RzuemRACD762FVkhMG1cM2ZMaCq4w1qxeS86mpK45Zl1+MxpM3GkS4MfGE5jw64BPLa6G4tOn2VrP4rhdAbregYQI0IqEcMsj0HBklU70NGURFd3P46e0pYnXDNZga179mN8SwO27tmPIye15t2rodEMurr78dCKbfjgnCl2pJzOSDqLP7y+FUdPabPNlwCwbe9+LH5mPRLxGM48ehImt6dwyIRmxzm6+4Zww2NrcdYxkzA4ksEpsybYlgZVx9+8uhlrd/ShrTGJ2Qe34wMuIQZY7ykR5b1/UUFErwgh5gUeV8cC5XQA3xZCnCf/vhIAhBD/acozb948sXTpUtNuhmEYxoOwAqWeTV4vAziSiA4logYACwE8WOU6MQzDHLDUrVNeCJEmon8A8CiAOIDFQoiVAdkYhmGYMlG3AgUAhBB/BPDHateDYRiGqW+TF8MwDFNDsEBhGIZhIoEFCsMwDBMJLFAYhmGYSGCBwjAMw0RC3U5sLAYi2g/AL7S4A8Ben/2HANjosz9MGaXuD3NMUD2jOEcU11FqPStxHfzMo9sP8DMPuz/MMZV85kcLIdp8jrMQQhww/wD0BOy/uZT8IcsoaX8U9YzoHFFcR6nPoxLXwc+cn/kB/8wBLA0qSwhxwJm8gr4y8/sS84cpo9T9YY4p9TrDHBPFdZRaz0pcBz/z6PYD/MzD7g9zTC08cwcHmslrqQixHk258lcKrmd01EMdAa5nlNRDHYHK1jPsuQ40DeXmKuevFFzP6KiHOgJczyiphzoCla1nqHMdUBoKwzAMUz4ONA2FYRiGKRMHvEAhosVE1E1Eb2hp7yKi54loBRH9nojaZXqSiG6X6avVN1jkvieJaA0RLZf/JlWxng1EdKtMf42I3q/lOVmmdxHRDeT19aXq17Fs95KIZhDRE/L5rSSiy2X6BCJaQkRr5e94Lc+V8n6tIaLztPRy3sso61kz95OIJsrj+4noRldZZbmfEdexlu7lOUT0irxnrxDRWVpZZWubvoQJBRvL/wCcAeAkAG9oaS8DeJ/c/jyAq+X2RQDukdvNADYAmCX/fhLAvBqp52UAbpXbkwC8AiAm/34JwOkACMDDAD5Yg3Us270EcDCAk+R2G4C3AMwG8H0AV8j0KwBcJ7dnA3gNQArAoQDeBhCvwL2Msp61dD9bALwHwJcA3Ogqqyz3M+I61tK9PBHAVLk9B8CWct/LoH8HvIYihHgKQK8r+WgAT8ntJQA+oQ4H0EJECQBNAEYA7KvBes4G8JjM1w0rvHAeER0MoF0I8bywWt0dAD5WS3WMqi4+ddwmhHhVbvcBWA1gGoAFAG6Xh92O3H1ZAGsQMSyEWA+gC8D8CtzLSOoZVX2iqqcQYkAI8QyAIb2cct7PqOpYboqo5zIhxFaZvhJAIxGlyt02/TjgBYqBNwB8VG5/CoD66PZ9AAYAbIM1Q/UHQgi9A71VqsHfqpCKaarnawAWEFGCiA4FcLLcNw3AZi3/ZplWS3VUlP1eEtEsWKO8FwFMFkJsA6wXG5bWBFj3Z5OWTd2zit3LEuupqJX7aaIi97PEOipq8V5+AsAyIcQwqvOeA2CBYuLzAC4joldgqZ4jMn0+gAyAqbDMCl8nosPkvk8LIeYCeK/899kq1nMxrEa0FMB/A3gOQBqW+uum3GF+hdYRqMC9JKJWAL8B8E9CCD8t03TPKnIvI6gnUFv301iER1qk9zOCOgI1eC+J6DgA1wH4okryOKwi4bwsUDwQQrwphDhXCHEygLth2aMBy4fyiBBiVJppnoU00wghtsjfPgC/QmXMDZ71FEKkhRBfFUKcIIRYAGAcgLWwOvDpWhHTAWx1l1vlOpb9XhJREtYLe5cQ4n6ZvEOaCpT5pVumb4ZTc1L3rOz3MqJ61tr9NFHW+xlRHWvuXhLRdAC/BXCxEEL1UxV/zxUsUDxQkRtEFAPwbwB+JndtBHAWWbQAOA3Am9Jsc5DMkwTwYVimnqrUk4iaZf1AROcASAshVkl1uY+ITpOq+sUAHqilOpb7XsrrvgXAaiHED7VdDwJYJLcXIXdfHgSwUNqmDwVwJICXyn0vo6pnDd5PT8p5P6OqY63dSyIaB+AhAFcKIZ5VB1fjPbeJ2stfb/9gjZq3ARiFJdkvAXA5rAiLtwBci9wE0FYA/wvLAbYKwL+IXFTIKwBel/uuh4ywqVI9ZwFYA8up92cAM7Vy5sF6Cd4GcKPKUyt1LPe9hBW9I2T5y+W/CwBMhBUksFb+TtDyfFPerzXQomXKfC8jqWeN3s8NsII3+mU7mV3O+xlVHWvtXsIaoA1oxy4HMKncbdPvH8+UZxiGYSKBTV4MwzBMJLBAYRiGYSKBBQrDMAwTCSxQGIZhmEhggcIwDMNEAgsUhqkRiOhLRHRxAcfPIm1lZ4apNolqV4BhGGvSnBDiZ8FHMkztwgKFYSJCLuj3CKwF/U6ENZnzYgDHAvghrImxOwH8rRBiGxE9CWsNs3cDeJCI2gD0CyF+QEQnwFpVoBnW5LTPCyF2E9HJsNZBGwTwTOWujmGCYZMXw0TL0QBuFkIcD+vTBpcB+DGATwprPbPFAK7Rjh8nhHifEOL/ucq5A8C/ynJWALhKpt8K4B+FEKeX8yIYphhYQ2GYaNkkcusq/RLAN2B9/GiJXOk8Dmt5GsWv3QUQUQcsQfMXmXQ7gP/1SL8TwAejvwSGKQ4WKAwTLe61jPoArPTRKAYKKJs8ymeYmoFNXgwTLYcQkRIeFwJ4AUCnSiOipPx+hREhxF4Au4novTLpswD+IoTYA2AvEb1Hpn86+uozTPGwhsIw0bIawCIi+h9Yq8P+GMCjAG6QJqsErA+KrQwoZxGAnxFRM4B1AD4n0z8HYDERDcpyGaZm4NWGGSYiZJTXH4QQc6pcFYapCmzyYhiGYSKBNRSGYRgmElhDYRiGYSKBBQrDMAwTCSxQGIZhmEhggcIwDMNEAgsUhmEYJhJYoDAMwzCR8P8BC8nVk4MZPHsAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'].plot()" ] @@ -215,9 +2240,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 17, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 17, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xt8XHWd+P/Xey65X5qkSW9pKfRGaYFiSwFBRdAW8QKu8LOuC1Xxh7K4Xr7uruJecHXxK7oKsrviolQKugqCSlUUSwERubWFQim900LTW5Lm0lzn+v7+cc6kk3QymUkmyWTm/Xw85pHJZ845c+Y0Pe/5fN6fi6gqxhhjTCo8430CxhhjJg4LGsYYY1JmQcMYY0zKLGgYY4xJmQUNY4wxKbOgYYwxJmUWNIwxxqTMgoYxxpiUWdAwxhiTMt94n0CmTZ48WWfPnj3ep2GMMRPK5s2bm1W1dqjtci5ozJ49m02bNo33aRhjzIQiIm+ksp01TxljjEmZBQ1jjDEpSzloiIhXRF4Skd+6v1eLyHoR2e3+rIrb9iYR2SMiO0VkZVz5UhHZ6r52h4iIW14oIve75c+LyOy4fVa777FbRFZn4kMbY4wZnnRqGp8Dtsf9/mVgg6rOAza4vyMiZwCrgEXAZcD3RcTr7nMncD0wz31c5pZfB7Sq6lzgNuBW91jVwM3AecBy4Ob44GSMMWZspRQ0RKQeeC/wo7jiK4C17vO1wJVx5T9X1YCq7gP2AMtFZBpQoarPqrOIx70D9okd60HgUrcWshJYr6otqtoKrOdEoDHGGDPGUq1p3A78IxCNK5uiqocB3J91bvkM4EDcdg1u2Qz3+cDyfvuoahhoB2qSHKsfEbleRDaJyKampqYUP5Ixxph0DRk0ROR9QKOqbk7xmJKgTJOUD3efEwWqd6nqMlVdVls7ZDdjY4wxw5RKTeNC4AMish/4OXCJiPwEOOo2OeH+bHS3bwBmxu1fDxxyy+sTlPfbR0R8QCXQkuRYxoxINKo8sPEAgXBkvE/FmAllyKChqjepar2qzsZJcD+uqn8DrANivZlWAw+7z9cBq9weUafiJLxfcJuwOkTkfDdfce2AfWLHusp9DwUeBVaISJWbAF/hlhkzIq8cbOcfH3qFp3Y1j/epGDOhjGRE+DeBB0TkOuBN4GoAVd0mIg8ArwFh4EZVjX2duwG4BygGfu8+AO4G7hORPTg1jFXusVpE5OvARne7r6lqywjO2RgAjnUGAOgOhsf5TIyZWNIKGqr6JPCk+/wYcOkg290C3JKgfBOwOEF5L27QSfDaGmBNOudpzFCOdQUBCISiQ2xpjIlnI8JNXmp1g0av5TSMSYsFDZOXWrrdoBGyoGFMOixomLzUV9Ow5ilj0mJBw+Sllq4QYDUNY9JlQcPkpZYup/eU1TSMSY8FDZOXWrvdmoYlwo1JiwUNk5dauiwRbsxwWNAweSccidLe49Q0bJyGMemxoGHyTqxpCqymYUy6LGiYvNPqjtEAy2kYky4LGibvxPIZXo9Y7ylj0mRBw+Sd2MC+qRVF1jxlTJosaJi8E5uscPokCxrGpMuChsk7fTWNymJrnjImTRY0TN5p6Q5SXuijvMhnK/cZkyYLGibvtHYFqSotoMjntZqGMWkaMmiISJGIvCAiL4vINhH5N7f8qyJyUES2uI/L4/a5SUT2iMhOEVkZV75URLa6r93hLvuKuzTs/W758yIyO26f1SKy232sxpgRau0OUVXip8jvsZyGMWlKZeW+AHCJqnaKiB94WkRiy7Tepqr/Eb+xiJyBs1zrImA68JiIzHeXfL0TuB54DngEuAxnydfrgFZVnSsiq4BbgQ+LSDVwM7AMUGCziKxT1daRfWyTz3pDEYr8Xor8XsJRJRyJ4vNapduYVAz5P0Udne6vfvehSXa5Avi5qgZUdR+wB1guItOAClV9VlUVuBe4Mm6fte7zB4FL3VrISmC9qra4gWI9TqAxZtiCkSgFPg9FfufPvzdsTVTGpCqlr1ci4hWRLUAjzk38efelz4jIKyKyRkSq3LIZwIG43Rvcshnu84Hl/fZR1TDQDtQkOdbA87teRDaJyKampqZUPpLJY8FwlEKfhyK/F7CpRIxJR0pBQ1UjqroEqMepNSzGaWqaAywBDgPfcTeXRIdIUj7cfeLP7y5VXaaqy2pra5N+FmOC4Sh+r4dCn1vTsKBhTMrSashV1TbgSeAyVT3qBpMo8ENgubtZAzAzbrd64JBbXp+gvN8+IuIDKoGWJMcyZthONE/FahrWPGVMqlLpPVUrIpPc58XAu4Adbo4i5oPAq+7zdcAqt0fUqcA84AVVPQx0iMj5br7iWuDhuH1iPaOuAh538x6PAitEpMpt/lrhlhkzbKFwlAKvh0KfNU8Zk65Uek9NA9aKiBcnyDygqr8VkftEZAlOc9F+4FMAqrpNRB4AXgPCwI1uzymAG4B7gGKcXlOxXlh3A/eJyB6cGsYq91gtIvJ1YKO73ddUtWUEn9eYkxLhNsDPmNQNGTRU9RXgnATl1yTZ5xbglgTlm4DFCcp7gasHOdYaYM1Q52lMqgJha54yZrisc7rJO0G3ecp6TxmTPgsaJq+o6snjNKymYUzKLGiYvBKJKqo4NQ1LhBuTNgsaJq8EI06tol9OwxLhxqTMgobJK0F3yhC/15qnjBkOCxomr8SCRoFNI2LMsFjQMHklEBc0YtOIBCxoGJMyCxomr4TcnEahz4OIUOjz2Cy3xqTBgobJK32JcHf9jCK/15qnjEmDBQ2TV+IT4QBFfg8BS4QbkzILGiavxCfCwa1pWJdbY1JmQcPklZOChs+ap4xJhwUNk1fiB/eB0zxl4zSMSZ0FDZNX+moabk6j0BLhxqTFgobJKyfXNLzW5daYNFjQMHllYE2jyOexwX3GpCGV5V6LROQFEXlZRLaJyL+55dUisl5Edrs/q+L2uUlE9ojIThFZGVe+VES2uq/d4S77irs07P1u+fMiMjtun9Xue+wWkdUYMwIDE+GFfm9fmTFmaKnUNALAJap6NrAEuExEzge+DGxQ1XnABvd3ROQMnOVaFwGXAd93l4oFuBO4Hmfd8Hnu6wDXAa2qOhe4DbjVPVY1cDNwHrAcuDk+OBmTrtCA5im/R/qarIwxQxsyaKij0/3V7z4UuAJY65avBa50n18B/FxVA6q6D9gDLBeRaUCFqj6rqgrcO2Cf2LEeBC51ayErgfWq2qKqrcB6TgQaY9IWGDC4z+/19AUSY8zQUsppiIhXRLYAjTg38eeBKap6GMD9WeduPgM4ELd7g1s2w30+sLzfPqoaBtqBmiTHynn3Prufjftbxvs0ck4wbu4pAL9PCEd0PE/JmAklpaChqhFVXQLU49QaFifZXBIdIkn5cPc58YYi14vIJhHZ1NTUlOTUJo7vrt/FAxsPDL2hScvARLjf67HmKWPSkFbvKVVtA57EaSI66jY54f5sdDdrAGbG7VYPHHLL6xOU99tHRHxAJdCS5FgDz+suVV2mqstqa2vT+UhZqzsYoTtovXoyLRiO4vMIHo/zfcSap4xJTyq9p2pFZJL7vBh4F7ADWAfEejOtBh52n68DVrk9ok7FSXi/4DZhdYjI+W6+4toB+8SOdRXwuJv3eBRYISJVbgJ8hVuW08KRKMFwlK5geLxPJeeEItG+JDiA3yuErHnKmJT5UthmGrDW7QHlAR5Q1d+KyLPAAyJyHfAmcDWAqm4TkQeA14AwcKOqxr4y3wDcAxQDv3cfAHcD94nIHpwaxir3WC0i8nVgo7vd11Q15xv6u91xA90Bq2lkWjAc7UuCg1PTiESVaFT7ah/GmMENGTRU9RXgnATlx4BLB9nnFuCWBOWbgJPyIaraixt0Ery2Blgz1Hnmkh63WaozYDWNTAueVNNwnoeiUQo93sF2M8a4bER4ForlMrqteSrjAuFoXxIcTiTErYnKmNRY0MhCsWDRZYnwjAuGo33dbQF8XqdJKmSjwo1JiQWNLNRX07DmqYw7ORHu6Ss3xgzNgkYWigWNrmCEaNSaTTJpYCK8r3nKrrMxKbGgkYV64nIZPTYDa0adlAj3WfOUMemwoJGF4gf12ViNzAoOSIT7PNY8ZUw6LGhkofigYWM1MisYTpzTsKlEjEmNBY0sFN/V1sZqZFYwov1zGm7zlE1aaExqLGhkoX41Det2m1HBcKRfl1vrPWVMeixoZKEey2mMmsFGhFvzlDGpsaCRhSynMXoGJsL9scF91jxlTEosaGSh7mCE2Nx5XZbTyKjBEuHW5daY1FjQyELdwTDVpYWANU9lWmhAIjz2PBy1oGFMKixoZKHuYITJZQV9z03mDN7l1pqnjEmFBY0s1BOMMKnEj88j1uU2g1Q1QSLcRoQbkw4LGlmoOxSmpMBHaaHPJi3MoFgPKetya8zwpbLc60wReUJEtovINhH5nFv+VRE5KCJb3MflcfvcJCJ7RGSniKyMK18qIlvd1+5wl33FXRr2frf8eRGZHbfPahHZ7T5Wkwe6gxGKC7yUFnhtevQMCrq1iYIEOQ2bsNCY1KSy3GsY+KKqvigi5cBmEVnvvnabqv5H/MYicgbOcq2LgOnAYyIy313y9U7geuA54BHgMpwlX68DWlV1roisAm4FPiwi1cDNwDJA3fdep6qtI/vY2a0nGKHE76Wk0GcLMWVQrFttrEkK4ma5teYpY1IyZE1DVQ+r6ovu8w5gOzAjyS5XAD9X1YCq7gP2AMtFZBpQoarPqqoC9wJXxu2z1n3+IHCpWwtZCaxX1RY3UKzHCTQ5rSsQprTQR2mBl04bp5ExfTUN34llXfsWYbLmKWNSklZOw202Ogd43i36jIi8IiJrRKTKLZsBHIjbrcEtm+E+H1jebx9VDQPtQE2SY+W0npDbPGU5jYw6ETQsp2HMcKUcNESkDHgI+LyqHsdpapoDLAEOA9+JbZpgd01SPtx94s/tehHZJCKbmpqakn6ObBeKRAlF1GmeKvBZTiODghHnWibqPWVdbo1JTUpBQ0T8OAHjp6r6SwBVPaqqEVWNAj8ElrubNwAz43avBw655fUJyvvtIyI+oBJoSXKsflT1LlVdpqrLamtrU/lIWSs2LsOpaXgtp5FBgb5E+InvIiKC3yuEraZhTEpS6T0lwN3AdlX9blz5tLjNPgi86j5fB6xye0SdCswDXlDVw0CHiJzvHvNa4OG4fWI9o64CHnfzHo8CK0Skym3+WuGW5azYZIUlBT6npmHNUxkTS4TH1zTAaaKy5iljUpNK76kLgWuArSKyxS37CvAREVmC01y0H/gUgKpuE5EHgNdwel7d6PacArgBuAcoxuk19Xu3/G7gPhHZg1PDWOUeq0VEvg5sdLf7mqq2DO+jTgyxaUNKCryUFXrpskR4xpzocuvtV+7ziE1YaEyKhgwaqvo0iXMLjyTZ5xbglgTlm4DFCcp7gasHOdYaYM1Q55krTtQ0nJxGTyhCJKp4PYn+CUw6EiXCY7/b1OjGpMZGhGeZ7rjmqdJC5xtxT8hqG5kQCDvXsch/cvOU5TSMSY0FjSwTS3wXuzUNsOnRM6U3FJtGpH/zlJPTsOYpY1JhQSPLxDdPVRT7AWjvCY3nKeWMWE2jcEDzlM8r1jxlTIosaGSZrrigUVPqTI/e0hUcz1PKGbEut0X+/jWNAq/HphExJkUWNLJMT1zzVLUFjYzqDSWuafi9HsI2YaExKbGgkWVi7e7F/hNB45gFjYyI1TQKT0qEi43TMCZFqYzTMGPoRLu7ty9h29JpQSMTAkkS4UFrnjImJVbTyDKBcBQR59tvgc9DeZGPlq7AeJ9WTgiEI/i9ctKYFxsRbkzqLGhkmUA4SqHPg7s+FTWlBbR0W++pTOgNRU+qZUCsecpyGsakwoJGlgmG+9/YqksLrKaRIYFw5KQkOFhNw5h0WNDIMgNvbNWlhRyznEZGBMLRk7rbAvh9FjSMSZUFjSwTCEX7zY1UU1pgXW4zJNb0N5DfJiw0JmUWNLLMwBtbVWkBrd1BnJnizUj0hiInTVYI1jxlTDosaGQZp3nqRBNKTWkBoYjSYfNPjVjy5ikLysakwoJGlgmEo/0Gn/WNCre8xogFQokT4QVW0zAmZRY0sszA5qnqMhsVnim94SiFiWoaNiI8bZ2BsOXa8lQqy73OFJEnRGS7iGwTkc+55dUisl5Edrs/q+L2uUlE9ojIThFZGVe+VES2uq/d4S77irs07P1u+fMiMjtun9Xue+wWkdXkuEA4SsGA5imw+acyIRCKUJSgpuGzmkba/vXhV1m95oXxPg0zDlKpaYSBL6rqQuB84EYROQP4MrBBVecBG9zfcV9bBSwCLgO+LyKxu+CdwPU464bPc18HuA5oVdW5wG3Are6xqoGbgfOA5cDN8cEpFw1sQok1T7Va0Bix4KA1DSenYZ0NUvdKQzs7j3YQtYke886QQUNVD6vqi+7zDmA7MAO4AljrbrYWuNJ9fgXwc1UNqOo+YA+wXESmARWq+qw6/zvvHbBP7FgPApe6tZCVwHpVbVHVVmA9JwJNTgpG+jdP1ZQWAtY8lQmDdbkt8Dqj722m29SEIlH2N3cRDEdp7LCBp/kmrZyG22x0DvA8MEVVD4MTWIA6d7MZwIG43Rrcshnu84Hl/fZR1TDQDtQkOVbOCgyY6qK4wEuR32OjwjOgd5BEuN/rlFkTVWreONbVF2DfbOke57MxYy3loCEiZcBDwOdV9XiyTROUaZLy4e4Tf27Xi8gmEdnU1NSU5NSy38DeU+DUNmxU+MgN1uXWFwsaYatppGL30c6+5xY08k9KQUNE/DgB46eq+ku3+Kjb5IT7s9EtbwBmxu1eDxxyy+sTlPfbR0R8QCXQkuRY/ajqXaq6TFWX1dbWpvKRslYgHKHA2/+fZUZVsf3nzIDB5p6KNU/Zkq+p2d3oBA0ROGB/l3knld5TAtwNbFfV78a9tA6I9WZaDTwcV77K7RF1Kk7C+wW3CatDRM53j3ntgH1ix7oKeNzNezwKrBCRKjcBvsIty1mJahpzast4vblrnM4oN0SiSiiig8xy61zvcNSCRir2NHYyY1Ix0yuLLWjkoVQWYboQuAbYKiJb3LKvAN8EHhCR64A3gasBVHWbiDwAvIbT8+pGVY24+90A3AMUA793H+AEpftEZA9ODWOVe6wWEfk6sNHd7muq2jLMz5r1VPWkWW4B5tSW0tIVpLUrSJXbm8qkJ7a4VZE/SU7DmqdSsruxk3lTyugNRawGnIeGDBqq+jSJcwsAlw6yzy3ALQnKNwGLE5T34gadBK+tAdYMdZ65INY8MrAJ5bTaUgBeb+5kaWn1mJ9XLjixal+CoOGWWfPU0CJRZW9TJxfNraGtO8Sfdk3sHKJJn40IzyJ9a1j7Tm6eAtjbaE1Uw3ViffAEzVPuSn7We2poDa3dBMNR5taVMau6hMaOAL2hyNA7mpxhQSOL9H0bHnBjq68qocDrYW9zZ6LdTApiNzbrcjsyzW4vvikVRcyqKQGcQGLyhwWNLBJrdy8c0HvK6xFmTy6xmsYIxGoag81yC9hMtynoDjqzLZcW+phZ7QQNy2vkl1QS4WaMBPuaUE6O5adNLmNXY8dYn1LO6AvICWsa1jyVqq6Acx1LCrx9sxUcbu8dz1MyY8xqGllksJwGOMnwN491241tmE5c28G73Nq1HVpPyK1pFPiYVOIHoK07NJ6nZMaYBY0skuzGNreujHBUWf/a0bE+rZwQy2kk7XJrQWNIfTWNQi9Ffi8lBV6bTDPPWNDIIoEkydoVi6ZyVn0ln/3ZS/zh1SNjfWoT3okut4nX0wDLaaSiL6dR4LRsV5UU0NJtQSOfWNDIIrFxAonWsS4r9PHTT57H7Mml/OBPe8f61Ca8QJJ8UYHVNFIWq2kUux0Kqkr91jyVZyxoZJFk34YByov8XDinht1HO2zthzT1NU8lyWnEOiKYwXUHw5QUePG4Y1uqSgpsgbA8Y0EjiyT7Nhwzf2o5XcEIDa09Y3VaOSHZtY2VBSxoDKkrGKGk4ETgrSopoM2ap/KKBY0skqxbaMyCKeUA7Dpq3W/TkezaxmofNrJ5aN2BMCUFJ3rqV5X4raaRZyxoZJFkvadi5k91gsZOCxppSXZtYwP+eixoDOmkmkZpAcd7w4QtH5Q3LGhkkVibeqJEeExFkZ/plUXsOmJBIx3JphGJlfWG7MY3lJ5ghNLC+JqGM+tyW48lw/OFBY0skkrzFDi1jZ1HbR6qdATCUQq8nr4EbjyPRyj0efq6PJvBdbmJ8JjYVP02ViN/WNDIIsmm7463YEo5exs7rYtoGpy11we/rkV+r+U0UtAdiPSN0QCodmsardbtNm9Y0MgigXAUr0f61qwezIKp5QQjUfbban4p6w1HEk6LHlPs91pOIwVdwTAlhSeuY2wqEUuG549UlntdIyKNIvJqXNlXReSgiGxxH5fHvXaTiOwRkZ0isjKufKmIbHVfu8Nd8hV3Wdj73fLnRWR23D6rRWS3+4gtB5uzBlvDeqB5dU4yPLZWsxna0DUNj+U0UtAdHFDTcJunrNtt/kilpnEPcFmC8ttUdYn7eARARM7AWap1kbvP90Uk9rXkTuB6nDXD58Ud8zqgVVXnArcBt7rHqgZuBs4DlgM3u+uE56xgOJo0CR4zp85ZyW+vBY2UBcKRpONfrHkqNV2B8EnjNACbSiSPDHmHUtWncNbtTsUVwM9VNaCq+4A9wHIRmQZUqOqz6gxlvhe4Mm6fte7zB4FL3VrISmC9qraoaiuwnsTBK2cEwsm/DceUFPiYMamYPU0WNFIVCEcTjgaPKfR76bXBfUlFokogHO03TqO4wEuR32NTieSRkeQ0PiMir7jNV7EawAzgQNw2DW7ZDPf5wPJ++6hqGGgHapIcK2c5QWPwG1u8OXVl7LGaRsp6Q8lrGsV+D71Bq2kkc2IBpv5/ozaVSH4ZbtC4E5gDLAEOA99xy0/uzwiapHy4+/QjIteLyCYR2dTUNHEXuk81pwEwt7aMvU2dRKM2B1UqhqrFFfm99IYtaCTTHYwtwNR/7baqkgLrcptHhhU0VPWoqkZUNQr8ECfnAE5tYGbcpvXAIbe8PkF5v31ExAdU4jSHDXasROdzl6ouU9VltbW1w/lIWSEQiib9Nhxvbl0ZvaEoB9tsDqpUBEKRpLW4Ip/lNIbSFRikplHqp9VyGnljWEHDzVHEfBCI9axaB6xye0SdipPwfkFVDwMdInK+m6+4Fng4bp9Yz6irgMfdvMejwAoRqXKbv1a4ZTkrGIn2TdM9lLl1ZQCW10hR94DpLway3lNDS1rTsJxG3hhyjXAR+RlwMTBZRBpwejRdLCJLcJqL9gOfAlDVbSLyAPAaEAZuVNXY17cbcHpiFQO/dx8AdwP3icgenBrGKvdYLSLydWCju93XVDXVhPyE5HQLTS2n0Rc0jnYyubSQxTMqcHsxmwS6gxGKkwSN4gKraQylr6Yx4DrWlBbQ3BkYj1My42DIoKGqH0lQfHeS7W8BbklQvglYnKC8F7h6kGOtAdYMdY65IhCO9E3LMJTq0gKqSwv47vpd9IQifG/VEq5YktP9BEakN5S8plHos8F9Q4nVNAYG38llhXT0hukNRfomfzS5y0aEZ5FUu9zGLJxWTjgapbzQx29fOTyKZzbxOc1Tg39HKvJ7+6ZxMYl19fWe6n8da8sLAThmyfC8MGRNw4ydYBpdbgG+fdXZBMJR7n12Pz99/k06A2HKCu2fdKBoVOkJRfqWKE2kyO8hGIkSiSreBJMamvicxsk1DYCmjgAzJhWP+XmZsWU1jSySbk1j+qRiTp1cyuVnTiMYjvL4jsZRPLuJK9aVNlnzVCygBKzb7aC6+3IaiWsazR2W18gHFjSySCAcSWkakYGWzqqirryQR6yJKqHB2uLj9S3EZAP8BtUVq2kM6HIbCxpNlgzPCxY0skg6vafieTzCRXMns+VA2yic1cQXCwRDNU8BNpVIEt3BMD6PnNQtvKbM6bxhNY38YEEji/QOMaleMrNqSjja0WvNKwkMNr4gXqymYd1uB9cVcHqgDezaXejzUlnst5pGnrCgkSWC4SihiA47kT2rugRVONhqI8QHis2ZlHxwnwWNoXQHwyf1nIqZXGZjNfKFBY0skUoTSjIzq0sAeLOlO2PnlCt60shpWNAYXFeSAZK15YU0WfNUXrCgkSW6BplBNFWz3KBxwGoaJxmsq2i8IrcDgk0lMrieJFOx1JYX0dxp4zTygQWNLHGiCWV4zVO1ZYUU+DwcsJrGSbpDKQQNq2kM6XhPiIoif8LXJpcVWE0jT1jQyBJdAedmNdyahscj1FcVW9BIoMcNyMVJAnKs2cVqGoNr6wn1rQk+UG15IZ2BsHVZzgMWNLJE1whrGuA0UVlO42SxG1lJsi63bldnm39qcG3dISqLE8+NFhsVbsnw3GdBI0t0x2oaIwgaM6tKrKaRQKx5KnkiPJbTsKCRiKrS3hNMWtOAEwP8OnpDOCscmFxjQSNL9NU0htk8BU5N43hvmHZb26CfnmAEEZJO0VJoOY2kekIRQhGlsniQoBE3/1R7d4i3fvNxblu/ayxP0YwRCxpZIpUePkOZWe1MFneg1Wob8bqDEUr8Jw9Ki3di7inLaSTS5n4RmTRI0Khzaxr7mrv49ZaDdPSG+Z+nXudIe++YnaMZGxY0skQqo5aHYmM1EnMWYEp+Xf1ewSM299Rg+oJGkuapc2dX8aM/7+Mnz73B7JoSoqr806+28t9P7GF/c9dYnq4ZRUMGDRFZIyKNIvJqXFm1iKwXkd3uz6q4124SkT0islNEVsaVLxWRre5rd7jLvuIuDXu/W/68iMyO22e1+x67RSS2JGxOis0gOpKaxpzaMrweYfvh45k6rZzQEwwPeV1FhCK/rd43mLYeZwzGYIlwEeGmyxfS3Blgd2Mnn3zbaVx7wWw27Gjk24/uZO2z+8fuZM2oSqWmcQ9w2YCyLwMbVHUesMH9HRE5A2e51kXuPt8Xkdj/1juB63HWDZ8Xd8zrgFZVnQvcBtzqHqsaZ2nZ84DlwM3xwSnXdAWdGW79Ka4RnkiR38u8ujJeaWjP4JlNfEOtDx5T5Pf2TaNu+ovlyQbLaQDTd3ZyAAAgAElEQVS8ZVYV7z1zGqUFXj6wZDpffs/pbPjiO5hVXWID/3LIkHcoVX0KZ+3ueFcAa93na4Er48p/rqoBVd0H7AGWi8g0oEJVn1WnS8W9A/aJHetB4FK3FrISWK+qLaraCqzn5OCVM7qD4ZPWXh6Os+or2Xqw3XquxOkJJV8fPKbY77VxGoNo70nePBXz7avP4pHPvY2KIj9+r4c5tWXUlhfS0mVdcXPFcL/WTlHVwwDuzzq3fAZwIG67Brdshvt8YHm/fVQ1DLQDNUmOlZOcGURHvuremfWTaOkKcrDNphOJSbWmUej32DiNQbSlGDRKCnycUlPar6y6tIBjVtPIGZlOhCfqnqJJyoe7T/83FbleRDaJyKampqaUTjTbODOIZqCmMaMSgK3WRNWnOxih2D90QC7yeQlY0EiorTtEgdczrAk1nRlwLWjkiuEGjaNukxPuz9g6ow3AzLjt6oFDbnl9gvJ++4iID6jEaQ4b7FgnUdW7VHWZqi6rra0d5kcaX10p9PBJxenTyvF7ha0HLWjEpJIIB2eAnzVPJdbeE6Si2J+02/JgakoLae0OEo1ak2kuGG7QWAfEejOtBh6OK1/l9og6FSfh/YLbhNUhIue7+YprB+wTO9ZVwONu3uNRYIWIVLkJ8BVuWU7qDmQmp1Ho8zJ/SrkFjTipNk8VF1jvqcG0J5l3aijVpQVEotqXFzET25BfbUXkZ8DFwGQRacDp0fRN4AERuQ54E7gaQFW3icgDwGtAGLhRVWP/C2/A6YlVDPzefQDcDdwnIntwahir3GO1iMjXgY3udl9T1YEJ+ZzRHYwwqSRxd8Z0LZpeweM7JmYz3WjoCUX6ZrFNpsjn7RuPYPpr6w4NOrBvKLHlYI91BagqzczfuBk/QwYNVf3IIC9dOsj2twC3JCjfBCxOUN6LG3QSvLYGWDPUOeaCTOU0AKZPKqa5M0AwHKUgydQZ+SLZOhDxbJzG4Nq6Q0yfVDSsfWtKndHixzqDzK0bYmOT9eyOkiW6gpnpPQUwrdL5z93YYVM4BMNRwlFNufeU5TQSa+8JUTHimoYlw3OBBY0skamcBsCUCido2Lw/8Uu9Dh2QK4r8dPRa81Qibd1BJg0yGnwoFjRyiwWNLBCNKt2hCCWFmappOBMXHjluQaM7lPr0LNWlBRzvDROKWG0jXigSpSsYGXYivMrN1R2ztTZyggWNLNAbjqA6snmn4k21mkafdGYPjiVpW7vtG3G8VEeDD8bv9TCpxG8D/HKEBY0s0LfUa4aCRkWxj2K/14IGcc1TKfSeqna/Ebd2WRNVvLYU5p0aSnVpAS3WPJUTLGhkge4MLPUaT0SYWlnEYWueSmvK+apS56ZoN7f+2vtmuB1+0JhcWmhLweYICxpZIHZjy1SXW3CaqI5aTaMvIKcyYWG1NU8ldPS4c7OPLek6HNWlBZYIzxEWNLJApmsagFPTsKBBR69zbctS6GQQa56ym1t/B1udyS/rq0qGfYyaMmueyhUWNLJAX04jkzWNyiIaO3rzfr6fQ+5sv9NSGJjWlwi3m1s/Da3dlBf6RtQ8VVPmzD8Vtp5pE54FjSwwKjWNiiJCEc37b80H23ooL/JRUTT0Dc/v9VBe5LNvxAMcbOthRlXxiI4xuawAVWixpr8Jz4JGFojVNDLV5RacmgbA0TxPhh9s7WHGpNRveNWlBZbTGKChtYf6EQaNOjcf0njckuETnQWNLDBaNQ0g7/MaB9vSu+FVlVjb+0DpBt5Easudv8emDgsaE50FjSzQNQq9p2Jt+A2t3Rk75kR0sC39moYFjRPae0J0BMIjSoIDTKlwahr5XvPNBRY0skB7TwifRyjyZS5o1JYVUldeyJYDbRk75kRzvDdER2+Y6WkEjaqSAkuEx4l96RhpTiPWXbfRahoTngWNLHC4rYeplUV4POmvijYYEWHZ7Co27W/N2DEnmlhX0XRueNWlfkvWxjnR3XZkQaPQ56WqxG8zL+cACxpZ4FBbb1rfhlO19JRqDrb1cLi9J+PHngj6gkZazVOF9IaifdOP5LuGYVzDwdSVF/UNFDQTlwWNLHCovYfplcNb4CaZc2dXAeRtbeNg2/BqGmBdQ2MOtvVQ5Pf0jZYfibqKQmueygEjChoisl9EtorIFhHZ5JZVi8h6Ednt/qyK2/4mEdkjIjtFZGVc+VL3OHtE5A53HXHctcbvd8ufF5HZIznfbBSJKkfaR6emsXBaBcV+L5v25+wquUkdauuhwOdhcmnq01/EpvFusRlZASenUV9VgvtfckTqyototET4hJeJmsY7VXWJqi5zf/8ysEFV5wEb3N8RkTNw1v9eBFwGfF9EYpnfO4HrgXnu4zK3/DqgVVXnArcBt2bgfLNKU0eAcFRHJWj4vR7OmTWJTW/kZ02joc2pwaWTK4p9o7aahuONY90jzmfETKkopKkjkPezFEx0o9E8dQWw1n2+FrgyrvznqhpQ1X3AHmC5iEwDKlT1WVVV4N4B+8SO9SBwqWTiK08WOeTmG4a7/vJQlsycxI4jHQTD+Td9w8HW9Ecy21QiJ4QiUfY2dbJganlGjldXXkg4qhaQJ7iRBg0F/igim0XkerdsiqoeBnB/xpaSnwEciNu3wS2b4T4fWN5vH1UNA+1AzcCTEJHrRWSTiGxqamoa4UcaW7G5kUajpgEwb0oZkajyZkvXqBw/W4UiUXYd7WBObVla+00uc5qybNVD2NvUSSiinDGtIiPHiy1DbKPCJ7aRBo0LVfUtwHuAG0Xk7Um2TVRD0CTlyfbpX6B6l6ouU9VltbW1Q51zVhntoBG7ae5pzK+g8erBdrqDEc479aTvGElVFvupLS9k99HOUTqziWPH4Q4ATp+amaBRVxEbq2EBeSIb0bwVqnrI/dkoIr8ClgNHRWSaqh52m54a3c0bgJlxu9cDh9zy+gTl8fs0iIgPqARyKqt7qK2X8sLUJtQbjtPcoLG3Kb9ugs/vc/5Mlp9anfa+86eUsaexI9OnNOFsP3ycAq+H02pLM3K8uvKJVdN4cmcjv3n5MKFIlH9YuYCZ1SMbFZ8rhl3TEJFSESmPPQdWAK8C64DV7margYfd5+uAVW6PqFNxEt4vuE1YHSJyvpuvuHbAPrFjXQU87uY9csahtp6Upu0errJCH9Mqi9jbmGdB4/VjzKktHdbCQfPqytnd2Jn3CdvtRzqYW1eG35uZ1OeJUeHZX9PYfbSDj/14Ixt2HGXD9qNcdvtTedsLcaCR/DVMAZ4WkZeBF4DfqeofgG8C7xaR3cC73d9R1W3AA8BrwB+AG1U1NoLqBuBHOMnxvcDv3fK7gRoR2QP8H9yeWLnkUHvPqDVNxcypLcubmkYgHCEYjrJpfyvnnZZe01TMvClldAcjfeM88tX2w8dZmKF8BkCR38ukEv+EmETz11sO4hH44xfezqNfeDsiwsNbDg29Yx4YdvOUqr4OnJ2g/Bhw6SD73ALckqB8E7A4QXkvcPVwz3EiONzWy1n1k0b1PebUlvLQiwdR1Yz0t89Wqsrl3/szrd3OJHvnDaNpCmD+FKe30O7GjrxtkmjuDNDUEWDhtMz0nIo5bXIpe7K81quqPLzlEBfOndzXpLZgajk7j1qTJdiI8HHV3hPiWFcwI1M0JDO3rozOQDjnR+Pubepkb1MXpYVeyot8vHXO5GEdZ36dc6PclcfJ8NcOHQfIaE0D4PRpFew40kE2tzK/+GYrDa09XLlkRl/Z/Cnl7Dqa3ec9VixojKPNbzhtpG+ZVTXEliNzogdVbt8En97dDMD/fvJ8tn515bDyGQCVJX7qygvZlUffLH+x6QCb4waBPrrtCEV+D0tmZrYWvHBqOe09oazu0vzLFw9S6POwcvHUvrIFU8po6w7ZeiBY0BhXz+9rwe8Vzpk1ys1TdU7Q2J3jN8G/7D3GrOqSjDQpzZ9Snjfdbps7A3z5l1u5/t5NNHcGCEWiPLL1MO9aOIXSwswtDAYnai7bDx/P6HEzpTMQ5tcvHeR9Z02nLO6zz3cHOFoTlQWNcbVxXwtn1U+iyJ+5dTQSqSsvZHplEc++fmxU32c8hSNRntt7jAvnDq9JaqAFU53miHyY7fa3Lx8iElXae0Lc9MutPLmzidbuEFfENc9kSuzmu/1wdt58f/XSQbqCEa654JR+5QvcPNfOI2N73i++2cpnf/YSx3tDY/q+yVjQGCc9wQivNLQPaxxBukSESxdO4aldzfSGcvMm+MrBdjoCYS7KUNB454I6AuEof9o1sWYYGI5fbznE6VPL+dJlp7P+taPc8JPNVBb7ecf8zA+UrSjyU19VzI4xvvmmQlX5ybNvcOaMSs6ur+z3Wk1ZIZPLCsa0yTISVW56aCvrXj7Ev/761TF736FY0BgnL73ZSjiqYxI0AN51xhR6QhGe3ZubtY3HXjuKR+CCOcPrZjvQeadVM6nEzx9ePZyR42WrN451seVAG1eeM4NPvu1U7rpmKQumlrP6glMo8I3O7eH0qRXsyILmqfbuUL+lfV/Y18LOox1cc/4pCXsZzp9Szs4RNlkeae9NuZn4oRcb2Hm0gwtOq+HXWw7x8JaDI3rvTLGgMU6e39eCR2DpKaObBI85/7RqSgu8rN9+dEzebyxFosovXzzIxQvqMrLuAzgzBL974RQ2bG8kEM7N2lkwHOVfH96G1yN84OzpiAgrFk3ld599G/9nxYJRe9+F08p5vblr3Gu9n1i7kavufIZQxJnM877n3qCiyMf7z56ecPv5U8rZdaRj2Ofd3h3iQ3c+w/v+82k27m9BVQcdQNoTjPDdP+7i7JmTuO+65SyZOYmv/3Y7nYHwsN47kyxojJNn9jazeEblqE0fMlChz8vb59fy2GtHCUdya8bbP+9u4sjxXq5eWj/0xmm4/MxpdATC/GVPc0aPmw2iUeXvf/Eyf9rVxDc+uHjUB5jGW3pKFZGo8sjW8avFvfRmK5vfaOX15i5+samBxo5e/vDqEa5eNpPigsQ5xpWLptITirD2mf0JXz/WGejLgUWiys4jHTyxs5FIVFFVvviLLTR29FJbXsjHf7yRc2/ZwKKbH+XaNS/wyNbD/QLImr/s48jxXr7yntPxeT189QOLaO4M8P0n9mT8WqTLgsY46AqEeenNtowlbVN1xZIZNHYEeOTVI2P6vqMlGlVePtDGPc/sp6rEz6ULp2T0+G+dW0NNaQE/fGpfzvXP/96G3ax7+RD/eNkCPnzurDF977fPq2XhtAr+8/E94/YF5p5n9lNW6OPs+kq+t2EXNz20lXBU+eh5g1+LC+bU8M4FtfzXE3tOmjq/tSvIytv/zEd/9BzhSJS//uFzrLz9KT7+44384E97+clzb/DY9kb+6fKF/Oz/P5+lp1Rx0dwarlpaz77mTv72py/ywe//hbbuIMc6A9z55F7etXBK36wGS2ZO4q/OmcGP/ryPu57aO65LHWS2P51JyQv7WwhHlQuHOfhsuFacMYW5dWV8/4k9vO/MaWktTpRtolHlxv99kd+7AfC6i07NeBt8oc/LZy6Zy7/95jWe2t08Konh8fDApgN8b8NuPvSWem54x5wxf3+PR/jcpfP49E8287MX3uSaC2aP6fs3Hu/ld68c5poLTuHyM6fx4f95lr/sbea6i07tm+BzMDddvpDLbn+KOx7fzc3vX9RXfssj22nuDNDcGeBjP97I8/ta+OK757P1YDu3P7YLn8fDO+bXsvqtsxER1n5ied++kajyq5cO8pVfbmX1mhfoCkboCUX48nv6NxH+03sX0tYT4huP7OCRrUf48cfO7Vv/ZSxZTWMc/GV3MwU+D8tmj00+I8bjEf724jnsONLBhh2NQ++QxW59dAe/f/UIn71kLr++8UK+dNnpo/I+f33eLGZWF/PN3+/oa/ueyP7r8d3844OvcOHcGm754OJxm1Zm5aIpLD2lin95eBufuGfjmLbVP/TiQcJR5doLZnPu7Gpe+Kd38epXV/Iv7ztjyH3nTynnw+fO4r5n32B/s7PcwJM7G3lwcwM3XDyHc2ZN4uk9zVy8oJbPXDKXb37oLCqLC/B7hVs/dFbC6+31CFctred7q5aw9WA7PcEIP/7Yucyt6z+FS01ZIWs+di7//ddv4bXDx/n//udZ2sZhQSsLGhkQikT57M9e4jP/+2LSFd/2NXfxi00HeHJXE8tOqRr18RmJfODs6dSWF7Lu5Yk7+dpTu5r4nz+9zkfPm8UX3j2fJTMnjVpPn0Kfl3+6fCHbDx/nlt9tH5X3GCsvvtnKf/xxF1csmc6PP7Z8XP7+YkSEn37yPL78ntN5Ymcj3/rDjjF5X1XlF5sPcO7sKk6d7Ez5PrmsEF8aM/l+4d3zKPB5+PffbWfnkQ4+f/8WFkwp53OXzuMbHzyTixfU8u9XOgG5urSAh264gIdueCtTK5PPZv2eM6fxxy+8gz9+4e28PUmt9r1nTeOej53LvuYu/vXhbX2f64mdjfz6pdHvYWXNUyOkqvzzr15l3cuH8HmE515v4ez6Si5eUMvfxHXdi0SVG36yua9/+l+9JfMDp1Lh83p429zJ/GlXE9GoTrgmqo5eZwDanNpS/uV9Z4zJN+XLFk/jExeeypq/7OPc2dW896xpo/6eo+G7f9xFTWkB3/jgmaMWZNNR5Pfy6XfM4ejxXn78l/28/+zpnDt7dLugv3Sgjdebuvj024ffLFdXXsSN75zLtx/dyWPbj1Je6OMH1yylyO9l4bQK7vn48n7bn1KT+nokc+tSW2nyrXMn87lL5/Gd9btQYMfh4+xu7GTR9AquWDJ9VP9fWNAYgd1HO7h53Tae2XuMv7tkLisXTeX2x3azr7mTf3m4kS0H2vnGXy2m0Ofloc0N7DjSwVfffwbTJhXztnljm8+I99a5k/nlSwfZcaSDM6ZndkK60XKkvZe//uFzvN7chQg8+OkLxvSb8k2Xn86zrx/ju+t38p7FUydEsA1FooQjSnGBlz/vbuLpPc3883sXZnxqkJH6+xULWP/aUT65dhO3r1rCOxfUDb1TmoLhKA9vOcgvNjdQ7Pdy+QgD/99ePIelp1Tx9O5mLpw7ua/WMpZuuHgOT+1uYsP2oyyeXsl3rj6b9589ugEDLGgMWyAc4aM/ep5AOMrXrljUNyDoR6uXoarcsWEPtz22iwMt3Xz8wtl8+487OWfWpL5E2Hi6cK7TI+OZvc3jHjQiUaUzEKakwMvxnhCvN3ex80gHJQVeFkwtZ9H0SlSVLz30Cofbe/n8u+bxlllVLD1lbAZFxvi9Hj79jtP43M+38OSuRi45PbM9tUZCVXl27zEeefUwF82t5dzZVTy1u4n/eHQXzZ0B3nfWdH77yiFmVZfwN+efMvQBx1hpoY+ffvI8PnXfZj5xz0a+9aGzuHrZzCH3U1V+9Od9bH6jle5QhDm1pVxz/iknJbNVla/8aisPbm6g0Ofh7y6Z229eqeEQEc4/rYbzh7lmSyb4vB4e+NQFqDKmX2Ik17oSLlu2TDdt2jSsfXtDkZS/vT64uYG//8XL3PuJ5YO2P657+RB//4uXCYaj1JU7SazFMyoTbjvWLv3Ok8ysLjmpKj2Wdhw5zqfv28z+Y92DbrNoegXlRT6ee72Fr12xiGvHuKdNvFAkytu/9QQzq0u4//rzxz34Axxu7+GLD7zMM3uP4fUIkbi+/gunVTCntpTfvnKYi+ZO5vZVS5hcNryZf8dCTzDC9fdt4i97mrl91Tl8YJBBdjFrn9nPzeu2cerkUkoLvew+2snkskJ+83cXsetoB/ubuwhGorze1MU9z+zn7y6ZyxfeNX9C1BLHg4hsVtVlQ243EYKGiFwGfA/wAj9S1W8Otu1wg0ZHb4i3f+sJ3rmgjg++ZQbzp5RTV16Y8Magqlx+x9NEolEe/fzbk948Xjt0nIbWbt55el3Gls3MhJsffpUHNjXw8Gcu5LTJpXQFIlSWJB9o2NjRy0+fe5PGjgAXzKlh8fQKZlWXpJVEBKeWtvaZ/dy2fjdlRT4+ceGphCJRKop81FeVsHB6BYFQhD/vbmbdy4foDUVYekoVX33/onH/D//jv+zj337zGm+bN5l/WLmAhdMq+v5dx3KRq3Akys83HuDbj+4kFInypctO56ql9fxpVxMHW3tYPKOS5adW4/UIjcd7mVxWOO7XLhXdwTAfW7ORjW+08NlL5jGtsgif18OKRVOoKPLT3hPisdeO0twZ4Dvrd3HhnBrWfOxcRIRXGtq46gfPUuTzcLy3f2+s9yyeyn//9VsmxDUYLzkTNETEC+zCWTq2AdgIfERVX0u0/XCDRnNngDs27OaXLx7s6/5XWuBl9uRSZk8upcDroSsQpicU4XhPiJcb2vnmX53JquVjOzAqU3Ye6eCjP3qOzkCYYr+Xtp4Qf79iATOrS3hocwN+r3C4vZcDLd1cd9Fp+H3C7Y/tJhSJUlbgo8O9RrXlhdx48RymTSqmvSfxTJyFPg8Lp1VQWujjxTda+dajOzjQ0sPFC2q59UNnMaVi9NZIzzRV5SfPv8k3frednlCE8iIfn3nnXA619fCzjQcoL/RRV1FEXXkhUyoKqSsvYtqkIt65oI6pFUXsP9bF0eMB9h/rYvvh4yw9pYolMyfx3OvH+gZsHWzrZVplEYtnVOIRiKri93qYP6Wc7YeP84vNDTyxo5HD7b0sP7Wab/7VmUOOL5hIekMR/uHBV/hNXA+/Aq+H2vJCjnUF6A0516m+qphf/e2F/dZNeXjLQf77iT1ce8Fs3nl6HQVeDwU+DxVFvqyoGWazXAoaFwBfVdWV7u83Aajq/020/Uiap8AZrb35jVb2H+vi9aYu9jV3sf9YF1FVSgt8FBd4KS3wMbWyiH+/cvG4dlscqcaOXv7vIztQVXpDUf6wzRkod0pNCaUFPqpK/RT6vDzujum4bNFUvvSe05lVXcJrh46z44hzA3thX0ta73v61HL+6b0Ledu8iTtYrqkjwHOvH+OhFxt4cmcTHoErz5lBkd9L4/Fejh4P0NjRS1NHgKiCCJQW+PqNRyjweRKO7B3YzBTjBBDny8zb5tXywbfMYMUZU3LyZqiqvPhmGzWlBbT1hPjDq0do7OilosjPB8+ZweyaUkoKvVlVe5/ociloXAVcpqqfdH+/BjhPVT8Tt831wPUAs2bNWvrGG2+My7lOZKrK/77wJiUFXj5w9gy8cdX4P+1qIhSO8q4zTk7+qiqvHnRmLJ1U4ifR/aszEGb74eMEw1Fqywt5x/y6fsef6Da/0UJ5kb9vbfF4kajyxrEu1r18iObOAGfVT6J+UjHTJxUzs7qEp/c0s7+5iwvn1lBR7Ad1xg0cbOth19EOPCJ4PEJPMMyrB48ztbKIK8+ZMeJErjED5VLQuBpYOSBoLFfVv0u0/UhrGsYYk49SDRoToW7XAMT3v6sHJu5wZmOMmcAmQtDYCMwTkVNFpABYBawb53Myxpi8lPUNo6oaFpHPAI/idLldo6rbxvm0jDEmL2V90ABQ1UeAR8b7PIwxJt9NhOYpY4wxWcKChjHGmJRZ0DDGGJMyCxrGGGNSlvWD+9IlIh3AzhQ3rwTaM/j2mT4ewGSgOUPHmgifN9PHzOT1g+y/hpk8nl27kZlI128yUKqqQ8/to6o59QA2pbHtXRl+74weL93PkyOfN9PnmLHrNxGuYSaPZ9cuf65fOuea781Tv8ny42XaRPi8dg2z63iZlO2fNZuvHWTJ583F5qlNmsL8KRNFrn2esWbXb/js2o3MRLp+6ZxrLtY07hrvE8iwXPs8Y82u3/DZtRuZiXT9Uj7XnKtpGGOMGT25WNMwxhgzSixojDERmSkiT4jIdhHZJiKfc8urRWS9iOx2f1a55TXu9p0i8l8DjvUREdkqIq+IyB9EZPJ4fKaxlOHr92H32m0TkW+Nx+cZS8O4du8Wkc3u39hmEbkk7lhL3fI9InKH5OLygQNk+PrdIiIHRKRzvD7PsGWyC5c9UurmNg14i/u8HGf98zOAbwFfdsu/DNzqPi8FLgI+DfxX3HF8QCMw2f39WzjL4o77Z5wg168GeBOodX9fC1w63p8vy67dOcB09/li4GDcsV4ALgAE+D3wnvH+fBPs+p3vHq9zvD9Xug+raYwxVT2sqi+6zzuA7cAM4AqcGxfuzyvdbbpU9Wmgd8ChxH2Uut/yKsiDxakyeP1OA3apapP7+2PAh0b59MfVMK7dS6oa+5vaBhSJSKGITAMqVPVZde6A98b2yWWZun7ua8+p6uGxPP9MsaAxjkRkNs63keeBKbE/IvdnXbJ9VTUE3ABsxQkWZwB3j+LpZp2RXD9gD3C6iMwWER/Of/SZQ+yTM4Zx7T4EvKSqAZwbZUPcaw1uWd4Y4fWb0CxojBMRKQMeAj6vqseHsb8fJ2icA0wHXgFuyuhJZrGRXj9VbcW5fvcDfwb2A+FMnmO2Svfaicgi4FbgU7GiBJvlTTfMDFy/Cc2Cxjhwb/gPAT9V1V+6xUfdaj/uz8YhDrMEQFX3uk0EDwBvHaVTzioZun6o6m9U9TxVvQBnvrLdo3XO2SLdayci9cCvgGtVda9b3ADUxx22njxoGoWMXb8JzYLGGHPzD3cD21X1u3EvrQNWu89XAw8PcaiDwBkiEptg7N04baw5LYPXDxGpc39WAX8L/CizZ5td0r12IjIJ+B1wk6r+Jbax2wTTISLnu8e8lhSu90SXqes34Y13Jj7fHjg9eRSnOWmL+7gcpzfPBpxvuxuA6rh99gMtQCfOt7wz3PJP4wSKV3DmkakZ7883wa7fz4DX3Meq8f5s2XbtgH8GuuK23QLUua8tA14F9gL/hTtQOJcfGb5+33L/FqPuz6+O9+dL9WEjwo0xxqTMmqeMMcakzIKGMcaYlFnQMMYYkzILGsYYY1JmQcMYY0zKLGgYM8ZE5NMicm0a288WkVdH85yMSZVvvE/AmHwiIj5V/cF4n4cxw2VBw5g0uZPV/UAu480AAAF8SURBVAFnsrpzcKbIvhZYCHwXKAOagY+p6mEReRJ4BrgQWCci5ThTYv+HiCwBfgCU4AyU+4SqtorIUmAN0A08PXafzpjkrHnKmOFZANylqmcBx4Ebgf8ErlLV2A3/lrjtJ6nqO1T1OwOOcy/wJfc4W4Gb3fIfA59VZ14sY7KG1TSMGZ4DemI+oZ8AX8FZaGe9u4idF4hfL+H+gQcQkUqcYPInt2gt8IsE5fcB78n8RzAmfRY0jBmegfPvdADbktQMutI4tiQ4vjFZwZqnjBmeWSISCxAfAZ4DamNlIuJ311EYlKq2A60i8ja36BrgT6raBrSLyEVu+Uczf/rGDI/VNIwZnu3AahH5H5zZTf8TeBS4w21e8gG34yzzmcxq4AciUgK8DnzcLf84sEZEut3jGpMVbJZbY9Lk9p76raouHudTMWbMWfOUMcaYlFlNwxhjTMqspmGMMSZlFjSMMcakzIKGMcaYlFnQMMYYkzILGsYYY1JmQcMYY0zK/h+14JK0aIFPzAAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -252,10 +2300,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 18, + "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", @@ -274,7 +2320,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 19, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +2344,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 20, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 20, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAD8CAYAAABQFVIjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X+QldWd5/H3B0FxElEg6PBDhY3GCppZHLrQLfdHxB0gPyrgjJmwOkrVWEXijy1nNrWiE7fMKFMVU5M4y7ohMWNG1FF0TCzdRIa0UWucWQSaoFE0TjMrgygjWI2KUwVj43f/eL53fLhpbt9ubve9t/vzqrrVT5/nnHNPPzT9vefH8xxFBGZmZkdrTLMbYGZmI4MDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNcTYZjdgOH3sYx+LmTNnNrsZZmZtZcuWLW9FxJT+8o2qgDJz5ky6urqa3Qwzs7Yi6R/ryechLzMzawgHFDMzawgHFDMza4i6AoqkHZJekPScpK5M+7qk1zPtOUmfLeW/UdJ2Sa9IWlhKn5v1bJe0SpIy/ThJD2b6RkkzS2WWSerO17JS+qzM251ljz36y2FmZoM1kB7KhRExJyI6Smm3Z9qciHgcQNJsYClwNrAI+I6kYzL/amA5cGa+FmX6lcC+iDgDuB24LeuaBNwMnAfMA26WNDHL3JbvfyawL+swM7MmGYohr8XA2og4GBGvAtuBeZKmAhMiYkMUu3rdAywplVmTxw8DF2XvZSHQGRE9EbEP6AQW5bn5mZcsW6mrre159wC/+70N7Nl/oNlNMTMbkHoDSgA/lbRF0vJS+rWSfiHpB6Wew3TgtVKeXZk2PY+r0w8rExG9wDvA5Bp1TQbezrzVdbW1VT/rZvOOHlY90d3sppiZDUi996FcEBFvSDoZ6JT0S4rhq1spgs2twLeA3wfUR/mokc4gytSq6zAZAJcDnHbaaX1laQln3bSOg70f/Ov3923cyX0bd3Lc2DG8svIzTWyZmVl96uqhRMQb+XUP8AgwLyLejIhDEfEB8H2KOQ4oegunlorPAN7I9Bl9pB9WRtJY4ESgp0ZdbwEnZd7quqrbfmdEdEREx5Qp/d7o2TTPXH8hX5gzjfHjin+S8ePGsHjONJ5ZcWGTW2ZmVp9+A4qkj0g6oXIMLABezDmRiouBF/P4MWBprtyaRTH5vikidgP7JZ2fcyBXAI+WylRWcF0CPJnzLOuBBZIm5pDaAmB9nnsq85JlK3W1pZMnjOeE48ZysPcDjhs7hoO9H3DCcWM5+YTxzW6amVld6hnyOgV4JFf4jgXuj4i/lnSvpDkUQ007gC8DRMQ2SQ8BLwG9wDURcSjrugq4GzgeWJcvgLuAeyVtp+iZLM26eiTdCmzOfLdERE8erwDWSloJbM062tpb7x3ksvNO59J5p3H/pp3s9cS8mbURFR/2R4eOjo7ws7zMzAZG0paqW0b65DvlzcysIRxQzMysIRxQzMysIRxQzMysIRxQzMysIRxQ2pCf92VmrcgBpQ35eV9m1opG1Z7y7c7P+zKzVuYeShvx877MrJU5oLQRP+/LzFqZh7zajJ/3ZWatys/yqsOedw9w7QNbuePSc90bMLNRx8/yaiCvqjIz65+HvGrwqiozs/q5h1KDV1WZmdXPAaUGr6oyM6tfXQFF0g5JL0h6TlJXpk2S1CmpO79OLOW/UdJ2Sa9IWlhKn5v1bJe0KrcCJrcLfjDTN0qaWSqzLN+jW9KyUvqszNudZY89+svxqyqrqh65+gIuO+909r53cCjexsys7dW1ykvSDqAjIt4qpX0T6ImIb0i6AZgYESskzQYeAOYB04AngE9ExCFJm4DrgGeBx4FVEbFO0tXAb0TEVyQtBS6OiC9JmgR0AR0UWw1vAeZGxL7cZvhHEbFW0neB5yNida2fwzs2mpkN3HCs8loMrMnjNcCSUvraiDgYEa8C24F5kqYCEyJiQxRR7J6qMpW6HgYuyt7LQqAzInoiYh/QCSzKc/Mzb/X7m5lZE9QbUAL4qaQtkpZn2ikRsRsgv56c6dOB10pld2Xa9DyuTj+sTET0Au8Ak2vUNRl4O/NW12VmZk1Q77LhCyLiDUknA52Sflkjr/pIixrpgylTq67DG1MEwOUAp512Wl9ZzMysAerqoUTEG/l1D/AIxfzImzmMRX7dk9l3AaeWis8A3sj0GX2kH1ZG0ljgRKCnRl1vASdl3uq6qtt+Z0R0RETHlClT6vlxzcxsEPoNKJI+IumEyjGwAHgReAyorLpaBjyax48BS3Pl1izgTGBTDovtl3R+zoFcUVWmUtclwJM5z7IeWCBpYq4iWwCsz3NPZd7q9zczsyaoZ8jrFOCRXOE7Frg/Iv5a0mbgIUlXAjuBLwJExLZcgfUS0AtcExGHsq6rgLuB44F1+QK4C7hX0naKnsnSrKtH0q3A5sx3S0T05PEKYK2klcDWrMPMzJrED4c0M7Oa/HBIMzMbVg4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEHUHFEnHSNoq6cf5/dclvS7puXx9tpT3RknbJb0iaWEpfa6kF/Lcqtxbntx//sFM3yhpZqnMMknd+VpWSp+Vebuz7LFHdynMzOxoDKSHch3wclXa7RExJ1+PA0iaTbEn/NnAIuA7ko7J/KuB5cCZ+VqU6VcC+yLiDOB24LasaxJwM3AeMA+4WdLELHNbvv+ZwL6sw8zMmqSugCJpBvA54M/ryL4YWBsRByPiVWA7ME/SVGBCRGyIYiP7e4AlpTJr8vhh4KLsvSwEOiOiJyL2AZ3Aojw3P/OSZSt1mZlZE9TbQ/kz4Hrgg6r0ayX9QtIPSj2H6cBrpTy7Mm16HlenH1YmInqBd4DJNeqaDLydeavrMjOzJug3oEj6PLAnIrZUnVoNfByYA+wGvlUp0kc1USN9MGVq1XUYScsldUnq2rt3b19ZzMysAerpoVwAfEHSDmAtMF/SfRHxZkQciogPgO9TzHFA0Vs4tVR+BvBGps/oI/2wMpLGAicCPTXqegs4KfNW13WYiLgzIjoiomPKlCl1/LhmZjYY/QaUiLgxImZExEyKyfYnI+L3ck6k4mLgxTx+DFiaK7dmUUy+b4qI3cB+SefnHMgVwKOlMpUVXJfkewSwHlggaWIOqS0A1ue5pzIvWbZSl5mZNcHY/rMc0TclzaEYatoBfBkgIrZJegh4CegFromIQ1nmKuBu4HhgXb4A7gLulbSdomeyNOvqkXQrsDnz3RIRPXm8AlgraSWwNeswM7MmUfFhf3To6OiIrq6uZjfDzKytSNoSER395fOd8mZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGY2Kux59wC/+70N7Nl/oNlNGbEcUMxsVFj1s2427+hh1RPdzW7KiHU0z/IyM2t5Z920joO9H27ldN/Gndy3cSfHjR3DKys/08SWjTzuoZjZiPbM9RfyhTnTGD+u+HM3ftwYFs+ZxjMrLmxyy0YeBxQza5rhmNc4ecJ4TjhuLAd7P+C4sWM42PsBJxw3lpNPGD9k7zlaOaCYWdMM17zGW+8d5LLzTueRqy/gsvNOZ+97B4f0/UYrP77ezIZd9bxGhec1WpMfX29mLcvzGiOTA4qZDTvPa4xMdQcUScdI2irpx/n9JEmdkrrz68RS3hslbZf0iqSFpfS5kl7Ic6tyb3ly//kHM32jpJmlMsvyPbolLSulz8q83Vn22KO7FGY2nDyvMfLUPYci6b8BHcCEiPi8pG8CPRHxDUk3ABMjYoWk2cADwDxgGvAE8ImIOCRpE3Ad8CzwOLAqItZJuhr4jYj4iqSlwMUR8SVJk4CufN8AtgBzI2Jf7lv/o4hYK+m7wPMRsbrWz+A5FDOzgWvoHIqkGcDngD8vJS8G1uTxGmBJKX1tRByMiFeB7cA8SVMpgtGGKKLYPVVlKnU9DFyUvZeFQGdE9ETEPqATWJTn5mfe6vc3M7MmqHfI68+A64HysoxTImI3QH49OdOnA6+V8u3KtOl5XJ1+WJmI6AXeASbXqGsy8Hbmra7rMJKWS+qS1LV37946f1wzMxuofgOKpM8DeyJiS511qo+0qJE+mDK16jo8MeLOiOiIiI4pU6b0lcXMzBqgnh7KBcAXJO0A1gLzJd0HvJnDWOTXPZl/F3BqqfwM4I1Mn9FH+mFlJI0FTgR6atT1FnBS5q2uy6wmP3XWbGj0G1Ai4saImBERM4GlwJMR8XvAY0Bl1dUy4NE8fgxYmiu3ZgFnAptyWGy/pPNzDuSKqjKVui7J9whgPbBA0sRcRbYAWJ/nnsq81e9vVpOfOms2NI7macPfAB6SdCWwE/giQERsyxVYLwG9wDURcSjLXAXcDRwPrMsXwF3AvZK2U/RMlmZdPZJuBTZnvlsioiePVwBrJa0EtmYdZkfkp86aDS0/esVGjT3vHmDl4y/z023/xIH3P2D8uDEsPPvX+drnPukb6sxq8KNXzKoM5O7sdpxnacc228jigGKjSr13Z7fjPEs7ttlGFg95mZW041Nw27HN1l485GU2CO34FNx2bLONTA4oZiXt+BTcdmyzjUxHs2zYbESqzLNcOu807t+0k71tMMndjm22kcdzKGZmVpPnUMzMbFg5oJiZWUM4oJiZDYBvID0yBxQzswHwDaRH5lVe1lb2vHuAax/Yyh2XnutlsTas/HDR/rmHYm3Fnw6tWXwDaf/cQ7G2MNo/Hbpn1ny+gbR/7qFYWxjtnw7dM2sN9T5cdLRyD8VaRq1P4aP10+Fo75m1mu9d/uG9fSuXnNPElrSmfnsoksZL2iTpeUnbJP1xpn9d0uuSnsvXZ0tlbpS0XdIrkhaW0udKeiHPrcqtgMntgh/M9I2SZpbKLJPUna9lpfRZmbc7yx7bmEtiA9WoZZT9fQofjZ8OR3vPzNpLPT2Ug8D8iHhP0jjgbyVVtu69PSL+tJxZ0myKLXzPBqYBT0j6RG4DvBpYDjwLPA4sotgG+EpgX0ScIWkpcBvwJUmTgJuBDiCALZIei4h9mef2iFgr6btZx+rBXwobrHIgWHnxpwZcvt5P4aPx0+Fo7ZlZe+q3hxKF9/Lbcfmq9QCwxcDaiDgYEa8C24F5kqYCEyJiQxQPELsHWFIqsyaPHwYuyt7LQqAzInoyiHQCi/Lc/MxLlq3UZcPkrJvWMfOGn3Dfxp1EFIFg5g0/4ayb1vVfuMSfwmsbjT0za091zaFIOgbYApwB/O+I2CjpM8C1kq4AuoCv5h/96RQ9kIpdmfZ+Hlenk19fA4iIXknvAJPL6VVlJgNvR0RvH3XZMHnm+guPuEf7QPhTeG2jsWdm7amuVV4RcSgi5gAzKHob51AML30cmAPsBr6V2dVXFTXSB1OmVl2HkbRcUpekrr179/aVxQapkYHAn8JHHj+iZPQZ0CqviHhb0tPAovLciaTvAz/Ob3cBp5aKzQDeyPQZfaSXy+ySNBY4EejJ9E9XlXkaeAs4SdLY7KWU66pu853AnVA8vn4gP6/1r1H7cPhT+MhztHNr1n763Q9F0hTg/QwmxwM/pZgQ3xIRuzPPHwLnRcRSSWcD9wPzKCblfwacGRGHJG0G/iuwkWJS/n9FxOOSrgE+FRFfyUn5346I381J+S3Ab2Zzfg7MjYgeSX8F/LA0Kf+LiPhOrZ/F+6GYDT3vcT/y1LsfSj09lKnAmpxHGQM8FBE/lnSvpDkUQ007gC8DRMQ2SQ8BLwG9wDW5wgvgKuBu4HiK1V2V2du7gHslbafomSzNunok3Qpszny3RERPHq8A1kpaCWzNOsysyRo1t2btp9+AEhG/AM7tI/3yGmX+BPiTPtK7gF8Zz4iIA8AXj1DXD4Af9JH+/yh6QWbWQrzIYvTynfJm1nDe43508p7yZmZWk/eUNzOzYeWAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYjaEvCeIjSYOKGZDqLwniNlI54dDmg2B6j1B7tu4k/s27vSeIDaiuYdiNgSeuf5CvjBnGuPHFf/Fxo8bw+I503hmxYVNbpnZ0HFAMRsC3hPERiMHFBtRWmkSvLInyCNXX8Bl553O3vcONrtJZkOq34AiabykTZKel7RN0h9n+iRJnZK68+vEUpkbJW2X9IqkhaX0uZJeyHOrJCnTj5P0YKZvlDSzVGZZvke3pGWl9FmZtzvLHtuYS2LtrJUmwb93eQcrl5zD7GkTWLnkHL53eb/bSZi1tX432Mo/+h+JiPckjQP+FrgO+G2gJyK+IekGYGJErJA0G3iAYnveacATwCci4pCkTVn2WeBxYFVErJN0NfAbEfEVSUuBiyPiS5ImAV1AB8Xe9VuAuRGxL/et/1FErJX0XeD5iFhd62fxBlsjV/UkeIUnwc2OXsM22IrCe/ntuHwFsBhYk+lrgCV5vBhYGxEHI+JVYDswT9JUYEJEbIgiit1TVaZS18PARRnIFgKdEdETEfuATmBRnpufeavf30YhT4IPr1YaWrTWUdcciqRjJD0H7KH4A78ROCUidgPk15Mz+3TgtVLxXZk2PY+r0w8rExG9wDvA5Bp1TQbezrzVdVW3fbmkLklde/furefHtTbkSfDh1UpDi9Y66roPJSIOAXMknQQ8IumcGtnVVxU10gdTplZdhydG3AncCcWQV195bGSoTIJfOu807t+0k73+9Nxwvr/GahnQjY0R8bakp4FFwJuSpkbE7hzO2pPZdgGnlorNAN7I9Bl9pJfL7JI0FjgR6Mn0T1eVeRp4CzhJ0tjspZTrslGqPOm9ckmtzzw2WM9cfyErH3+Zn277Jw68/wHjx41h4dm/ztc+98lmN81aQD2rvKZkzwRJxwP/Gfgl8BhQWXW1DHg0jx8DlubKrVnAmcCmHBbbL+n8nAO5oqpMpa5LgCdznmU9sEDSxFxFtgBYn+eeyrzV729mQ8RDi1ZLPT2UqcAaScdQBKCHIuLHkjYAD0m6EtgJfBEgIrblCqyXgF7gmhwyA7gKuBs4HliXL4C7gHslbafomSzNunok3Qpszny3RERPHq8A1kpaCWzNOsxsiHlo0Y6k32XDI8lQLhve8+4Brn1gK3dceq4/rZk1UDv+32rHNtfSsGXDVh+vejEbGu34f6sd29wI7qEcJd9QZzY02vH/Vju2uR7uoQwT31BnNjTa8f9WO7a5kRxQjtJAVr347mJrB63ye9qOK8rasc2N5IDSAPU+VXa0jqtae2ml39N2fGJzO7a5UTyHMgxG6riqjSz+PbUj8RxKCxnt46rWHvx72npaZfixXg4ow6Cdx1Xb7RfaBq+df09HqlYafqzHgJ7lZYPXrncXl3+hV178qWY3x4ZYu/6ejjTt+hBOz6FYnzyebtY8e949cMSHcDajx+g5FDsqHk83a552HX70kJf1qV1/oc1GinYcfnRAsSNqx19os5GiHff38RzKKDbSnohqZkPDcyjWr3Zbkmhmrc1DXqNQuy5JbDXu4Zkdrp4tgE+V9JSklyVtk3Rdpn9d0uuSnsvXZ0tlbpS0XdIrkhaW0udKeiHPrcqtgMntgh/M9I2SZpbKLJPUna9lpfRZmbc7yx7bmEsy8nkFV2O4h2d2uHp6KL3AVyPi55JOALZI6sxzt0fEn5YzS5pNsYXv2cA04AlJn8htgFcDy4FngceBRRTbAF8J7IuIMyQtBW4DviRpEnAz0AFEvvdjEbEv89weEWslfTfrWD34SzF6eAXX0XEPz6xv/fZQImJ3RPw8j/cDLwPTaxRZDKyNiIMR8SqwHZgnaSowISI2RLES4B5gSanMmjx+GLgoey8Lgc6I6Mkg0gksynPzMy9ZtlKX1WE0PxH1aLmHZ9a3Ac2h5FDUucBG4ALgWklXAF0UvZh9FMHm2VKxXZn2fh5Xp5NfXwOIiF5J7wCTy+lVZSYDb0dEbx91WR3acUliq3APz6xvda/ykvRR4IfAH0TEuxTDSx8H5gC7gW9VsvZRPGqkD6ZMrbqq271cUpekrr179/aVxWzA3MMz+1V19VAkjaMIJn8ZET8CiIg3S+e/D/w4v90FnFoqPgN4I9Nn9JFeLrNL0ljgRKAn0z9dVeZp4C3gJEljs5dSruswEXEncCcU96HU8/Oa9cc9PLNfVc8qLwF3AS9HxLdL6VNL2S4GXszjx4CluXJrFnAmsCkidgP7JZ2fdV4BPFoqU1nBdQnwZM6zrAcWSJooaSKwAFif557KvGTZSl1mZtYE9fRQLgAuB16Q9Fym/RHwXyTNoRhq2gF8GSAitkl6CHiJYoXYNbnCC+Aq4G7geIrVXesy/S7gXknbKXomS7OuHkm3Apsz3y0R0ZPHK4C1klYCW7MOMzNrEj96xczMavKjV8zMbFg5oJiZNcFI3F7bAcVsFBmJf8Ta1Uh8dI8fDmk2ipT/iK28+FPNbs6oNJIf3eNJebNRoPqPWMVI+CPWbhq5X/xwPfHak/Jm9q/8/LHW0chH97TasJmHvMyabDg+Zfr5Y63laLfXbtVhMwcUsyYbrnmNo/0jZo1ztI/ueeb6C484bNZMDihmTTLcnzL9/LGRo1V7nJ5DMWsSz2s0zmhcDt2KT7x2D8WsSVr1U2Y7Go3LoVuxx+mAYtZEntc4Oq06OT1a+T4UM2tbjbynw47M96GY2YjnYcPW4iEvM2trHjZsHR7yMjMbwRpx46yHvMzMbFgfz1LPnvKnSnpK0suStkm6LtMnSeqU1J1fJ5bK3Chpu6RXJC0spc+V9EKeW5V7y5P7zz+Y6RslzSyVWZbv0S1pWSl9VubtzrLHNuaSmJm1v7NuWsfMG37CfRt3ElGsgJt5w08466Z1/RcepHp6KL3AVyPik8D5wDWSZgM3AD+LiDOBn+X35LmlwNnAIuA7ko7JulYDy4Ez87Uo068E9kXEGcDtwG1Z1yTgZuA8YB5wcylw3Qbcnu+/L+swMzOac+NsvwElInZHxM/zeD/wMjAdWAysyWxrgCV5vBhYGxEHI+JVYDswT9JUYEJEbIhi4uaeqjKVuh4GLsrey0KgMyJ6ImIf0AksynPzM2/1+5uZjXrNWAE3oFVeORR1LrAROCUidkMRdCSdnNmmA8+Wiu3KtPfzuDq9Uua1rKtX0jvA5HJ6VZnJwNsR0dtHXWZmxvCvgKs7oEj6KPBD4A8i4t2c/ugzax9pUSN9MGVq1XV4Y6TlFMNsnHbaaX1lMTMbkYb78Sx1rfKSNI4imPxlRPwok9/MYSzy655M3wWcWio+A3gj02f0kX5YGUljgROBnhp1vQWclHmr6zpMRNwZER0R0TFlypR6flwzMxuEelZ5CbgLeDkivl069RhQWXW1DHi0lL40V27Noph835TDY/slnZ91XlFVplLXJcCTOc+yHlggaWJOxi8A1ue5pzJv9fubmVkT1DPkdQFwOfCCpOcy7Y+AbwAPSboS2Al8ESAitkl6CHiJYoXYNRFxKMtdBdwNHA+syxcUAeteSdspeiZLs64eSbcCmzPfLRHRk8crgLWSVgJbsw4zM2sS3ylvZmY1+U55aymjcQMks9HGAcWGxXA+/sHMmsNPG7Yh5Q2QzEYP91BsSHnfdLPRwwHFhpQ3QDIbPTzkZUPOGyCZjQ5eNmxmZjV52bCZmQ0rBxQzM2sIB5QRyjcSmtlwc0AZoXwjoZkNN6/yGmF8I6GZNYt7KCOMbyQ0GzwPFR8dB5QRxjcSmg2eh4qPjoe8RiDfSGg2MB4qbgzf2Ghmo96edw+w8vGX+em2f+LA+x8wftwYFp7963ztc590754G3tgo6QeS9kh6sZT2dUmvS3ouX58tnbtR0nZJr0haWEqfK+mFPLcqtwEmtwp+MNM3SppZKrNMUne+lpXSZ2Xe7ix7bD0XxcysLx4qbox65lDuBhb1kX57RMzJ1+MAkmZTbN97dpb5jqRjMv9qYDnFHvNnluq8EtgXEWcAtwO3ZV2TgJuB84B5wM25rzyZ5/aIOBPYl3WYmQ1aZaj4kasv4LLzTmfveweb3aS20+8cSkT8TbnX0I/FwNqIOAi8mnvEz5O0A5gQERsAJN0DLKHYU34x8PUs/zBwR/ZeFgKdlT3kJXUCiyStBeYDl2aZNVl+dZ1tNDP7Fd+7/MMRnZVLzmliS9rX0azyulbSL3JIrNJzmA68VsqzK9Om53F1+mFlIqIXeAeYXKOuycDbmbe6LjMza5LBBpTVwMeBOcBu4FuZrj7yRo30wZSpVdevkLRcUpekrr179x4pm5mZHaVBBZSIeDMiDkXEB8D3KeY4oOgtnFrKOgN4I9Nn9JF+WBlJY4ETgZ4adb0FnJR5q+vqq613RkRHRHRMmTJloD+qmZnVaVABRdLU0rcXA5UVYI8BS3Pl1iyKyfdNEbEb2C/p/JwfuQJ4tFSmsoLrEuDJKNYyrwcWSJqYQ2oLgPV57qnMS5at1GVmZk3S76S8pAeATwMfk7SLYuXVpyXNoRhq2gF8GSAitkl6CHgJ6AWuiYhDWdVVFCvGjqeYjF+X6XcB9+YEfg/FKjEiokfSrcDmzHdLZYIeWAGslbQS2Jp1mJlZE/nGRjMzq6neGxtHVUCRtBf4xz5OfYxibqaduM1Dr93aC27zcGm3Nh9te0+PiH4noUdVQDkSSV31RN9W4jYPvXZrL7jNw6Xd2jxc7fXThs3MrCEcUMzMrCEcUAp3NrsBg+A2D712ay+4zcOl3do8LO31HIqZmTWEeyhmZtYQIzKgHGEPl38raUPuyfJ/JE3I9HGS1mT6y5JuLJV5Ovd1qez7cnKLtPlYSX+R6c9L+nSpTJ/7zrR4m4flOks6VdJT+e+8TdJ1mT5JUmfur9NZetjpgPf3afE2t+R1ljQ5878n6Y6qulryOvfT5iG/zoNo729J2pLXcouk+aW6GneNI2LEvYD/CPwm8GIpbTPwn/L494Fb8/hSikfuA/waxZ3/M/P7p4GOFmzzNcBf5PHJwBZgTH6ruEO0AAADpklEQVS/Cfh3FA/RXAd8pg3aPCzXGZgK/GYenwD8PTAb+CZwQ6bfANyWx7OB54HjgFnAPwDHDOd1bnCbW/U6fwT498BXgDuq6mrV61yrzUN+nQfR3nOBaXl8DvD6UFzjEdlDiYi/oXiMS9lZwN/kcSfwO5XswEdUPGzyeOBfgHeHo51lA2zzbOBnWW4P8DbQoeIZaxMiYkMUvymVfWdats1D1ba+RMTuiPh5Hu8HXqbY+mAxxb465NfKNfvX/X0i4lWgsr/PsF3nRrV5KNrWqDZHxD9HxN8CB8r1tPJ1PlKbh8sg2rs1IioP0d0GjFfxzMWGXuMRGVCO4EXgC3n8RT58kvHDwD9TPIZ/J/Cn8eEzwwD+Irut/2Moh4+O4Ehtfh5YLGmsiodwzs1ztfadGS4DbXPFsF5nFZvGnQtsBE6J4gGm5NfKEMVg9vcZMkfZ5opWvM5H0srXuT/Ddp0H0d7fAbZGsRFiQ6/xaAoovw9cI2kLRRfxXzJ9HnAImEYxRPBVSf8mz10WEZ8C/kO+Lh/eJh+xzT+g+IfvAv4M+L8UD+Mc0F4xQ2SgbYZhvs6SPgr8EPiDiKjVG23InjyN0IA2Q+te5yNW0Udaq1znWobtOg+0vZLOpthC/cuVpD6yDfoaj5qAEhG/jIgFETEXeIBibBmKOZS/joj3cyjm78ihmIh4Pb/uB+5n+IcO+mxzRPRGxB9GxJyIWAycBHRTe9+ZVm3zsF5nSeMo/gP+ZUT8KJPfzK5/ZZhlT6YPZn+fVm1zK1/nI2nl63xEw3WdB9peSTOAR4ArIqLy96+h13jUBJTKSgtJY4CbgO/mqZ3AfBU+ApwP/DKHZj6WZcYBn+fDfV+a2mZJv5ZtRdJvAb0R8VLU3nemJds8nNc5r8ldwMsR8e3SqfKePOX9dQazv09LtrnFr3OfWvw6H6meYbnOA22vpJOAnwA3RsTfVTI3/BoPdja/lV8Un4x3A+9TROArgesoVkL8PfANPryp86PAX1FMVL0E/Pf4cBXHFuAXee5/kqtlWqDNM4FXKCbinqB4Emilng6KX+B/AO6olGnVNg/ndaZYlRP5Xs/l67PAZIoFA935dVKpzNfyWr5CafXLcF3nRrW5Da7zDooFHu/l79LsNrjOv9Lm4brOA20vxYe7fy7lfQ44udHX2HfKm5lZQ4yaIS8zMxtaDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQ/x/rFzgxQKYDrwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.plot(style='*')" ] @@ -314,9 +2383,55 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 21, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "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": 21, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,9 +2446,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 22, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 22, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEKCAYAAAAGvn7fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGbhJREFUeJzt3X2UJXV95/H3h5kBhmkYjAONDkr7QAjIqDgXXWQ13WhcdNCcGE5QQcVIGo0PRCdnM8v6sLrLOj5MsphgkklUiAodw8PZyBjUE2hQNEgPqA2OEBdmlSEMAjLSMAuMfPePX7XctP1wq27dvsXPz+ucPn3vrbpVn/rdut9b9auqexURmJlZPvbqdwAzM6uXC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLzNJeTnzVqlUxNDQ067AHH3yQFStW9HL2lTU5GzQ7n7NV1+R8zlZd2Xxbt269JyIO6mqmEdGzv7Vr18ZcrrrqqjmH9VuTs0U0O5+zVdfkfM5WXdl8wER0WXvdFWNmlhkXdjOzzLiwm5llxoXdzCwzLuxmZpkpVdglvUfSzZJuknSRpH17FczMzKrpuLBLWg28G2hFxNHAEuB1vQpmZmbVlO2KWQosl7QU2A+4s/5IZmbWDUWJ3zyVdBZwDrAb+GpEnDrLOKPAKMDg4ODasbGxWac1NTXFwMBAlcw91+RsUH++yR27apvW4HLYubvz8desXlnbvBfyq/a61snZqiubb2RkZGtEtLqZZ8eFXdKTgEuAU4D7gX8ALo6Iz8/1nFarFRMTE7MOGx8fZ3h4uGzeRdHkbFB/vqENW2qb1vo1e9g02fk3VWzfuK62eS/kV+11rZOzVVc2n6SuC3uZrpiXA7dHxE8i4lHgUuDF3czczMzqV6aw/wj4D5L2kyTgZcC23sQyM7OqOi7sEXEdcDFwAzBZPHdzj3KZmVlFpb62NyI+CHywR1nMzKwGvvLUzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMdF3ZJR0j6TtvfzyT9US/DmZlZeR3/NF5E3AI8H0DSEmAHcFmPcpmZWUVVu2JeBvyfiPi/dYYxM7PuKSLKP0n6DHBDRPzFLMNGgVGAwcHBtWNjY7NOY2pqioGBgdLzXgxNzgb155vcsau2aQ0uh527Ox9/zeqVtc17Ie3tVucylzHf8jZ5vXO26srmGxkZ2RoRrW7mWbqwS9obuBN4TkTsnG/cVqsVExMTsw4bHx9neHi41LwXS5OzQf35hjZsqW1a69fsYdNkxz18bN+4rrZ5L6S93epc5jLmW94mr3fOVl3ZfJK6LuxVumJeSdpan7eom5lZf1Qp7K8HLqo7iJmZ1aNUYZe0H/BbwKW9iWNmZt3qvDMUiIiHgCf3KIuZmdXAV56amWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZabsT+MdKOliST+QtE3Scb0KZmZm1ZT6aTzgXOCKiDhZ0t7Afj3IZGZmXei4sEs6AHgpcDpARDwCPNKbWGZmVpUiorMRpecDm4HvA88DtgJnRcSDM8YbBUYBBgcH146Njc06vampKQYGBqon75HJHbsYXA47dy/+vNesXtnReHW33eSOXbVNq2zbdbrMdWhvtzqXuYz5lrep7wlwtm6UzTcyMrI1IlrdzLNMYW8B/wIcHxHXSToX+FlEvH+u57RarZiYmJh12Pj4OMPDw+UT99jQhi2sX7OHTZNle6m6t33juo7Gq7vthjZsqW1aZduu02WuQ3u71bnMZcy3vE19T4CzdaNsPkldF/YyB0/vAO6IiOuK+xcDL+hm5mZmVr+OC3tE3AX8WNIRxUMvI3XLmJlZg5Ttb3gX8IXijJjbgLfUH8nMzLpRqrBHxHeArvp+zMyst3zlqZlZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmSv2CkqTtwAPAz4E93f6StpmZ1a/sb54CjETEPbUnMTOzWrgrxswsM4qIzkeWbgd+CgTw1xGxeZZxRoFRgMHBwbVjY2OzTmtqaoqBgYEqmXtqcscuBpfDzt2LP+81q1d2NF7dbTe5Y1dt0yrbdp0ucx3a263OZS5jvuVt6nsCnK0bZfONjIxs7babu2xhf2pE3CnpYOBrwLsi4pq5xm+1WjExMTHrsPHxcYaHh0vG7b2hDVtYv2YPmyar9FJ1Z/vGdR2NV3fbDW3YUtu0yrZdp8tch/Z2q3OZy5hveZv6ngBn60bZfJK6LuylumIi4s7i/93AZcALu5m5mZnVr+PCLmmFpP2nbwOvAG7qVTAzM6umTH/DIHCZpOnnXRgRV/QklZmZVdZxYY+I24Dn9TCLmZnVwKc7mpllxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWWmdGGXtETSjZIu70UgMzPrTpUt9rOAbXUHMTOzepQq7JIOBdYBf9ubOGZm1i1FROcjSxcDHwH2B/44Ik6aZZxRYBRgcHBw7djY2KzTmpqaYmBgoErmnprcsYvB5bBzd7+TzK3J+cpmW7N6Ze/CzNC+zk3u2LVo82033/I29T0BztaNsvlGRka2RkSrm3ku7XRESScBd0fEVknDc40XEZuBzQCtViuGh2cfdXx8nLmG9dPpG7awfs0eNk123DSLrsn5ymbbfupw78LM0L7Onb5hy6LNt918y9vU9wQ4Wzf6ka9MV8zxwGskbQfGgBMkfb4nqczMrLKOC3tE/JeIODQihoDXAVdGxGk9S2ZmZpX4PHYzs8xU6qiNiHFgvNYkZmZWC2+xm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMdFzYJe0r6duSvivpZkkf6mUwMzOrpswvKD0MnBARU5KWAd+Q9E8R8S89ymZmZhV0XNgjIoCp4u6y4i96EcrMzKpTqtcdjiwtAbYCzwbOi4g/mWWcUWAUYHBwcO3Y2Nis05qammJgYGDOeU3u2NVxrroNLoedu/s2+wU1OV/ZbGtWr+xdmBna17l+rV/zLe9C74l+eqJma8LrXLbtRkZGtkZEq5v5lyrsv3iSdCBwGfCuiLhprvFarVZMTEzMOmx8fJzh4eE55zG0YUvpXHVZv2YPmyYr/c73omhyvrLZtm9c18M0/177Otev9Wu+5V3oPdFPT9RsTXidy7adpK4Le6WzYiLifmAcOLGbmZuZWf3KnBVzULGljqTlwMuBH/QqmJmZVVNmf/4pwAVFP/tewBcj4vLexDIzs6rKnBXzPeCYHmYxM7Ma+MpTM7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8yU+c3Tp0m6StI2STdLOquXwczMrJoyv3m6B1gfETdI2h/YKulrEfH9HmUzM7MKOt5ij4h/i4gbitsPANuA1b0KZmZm1Sgiyj9JGgKuAY6OiJ/NGDYKjAIMDg6uHRsbm3UaU1NTDAwMzDmPyR27Sueqy+By2Lm7b7NfUJPzOVt1vc63ZvXKys9d6P3aT/Nl61cdaW/rsm03MjKyNSJa3cy/dGGXNABcDZwTEZfON26r1YqJiYlZh42PjzM8PDznc4c2bCmVq07r1+xh02SZXqrF1eR8zlZdr/Nt37iu8nMXer/203zZ+lVH2tu6bNtJ6rqwlzorRtIy4BLgCwsVdTMz648yZ8UI+DSwLSL+tHeRzMysG2W22I8H3gicIOk7xd+repTLzMwq6rhDLyK+AaiHWczMrAa+8tTMLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsM2V+8/Qzku6WdFMvA5mZWXfKbLGfD5zYoxxmZlaTjgt7RFwD3NfDLGZmVgNFROcjS0PA5RFx9DzjjAKjAIODg2vHxsZmHW9qaoqBgYE55zW5Y1fHueo2uBx27u7b7BfU5HzOVl2v861ZvbLycxd6v/bTfNn6VUfa27ps242MjGyNiFY386+9sLdrtVoxMTEx67Dx8XGGh4fnfO7Qhi0d56rb+jV72DS5tG/zX0iT8zlbdb3Ot33jusrPXej92k/zZetXHWlv67JtJ6nrwu6zYszMMuPCbmaWmTKnO14EfAs4QtIdkt7au1hmZlZVxx16EfH6XgYxM7N6uCvGzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDOlCrukEyXdIumHkjb0KpSZmVVX5jdPlwDnAa8EjgJeL+moXgUzM7NqymyxvxD4YUTcFhGPAGPAb/cmlpmZVaWI6GxE6WTgxIg4o7j/RuBFEfHOGeONAqPF3SOAW+aY5CrgniqhF0GTs0Gz8zlbdU3O52zVlc13WEQc1M0Ml5YYV7M89kufChGxGdi84MSkiYholZj/omlyNmh2Pmerrsn5nK26fuQr0xVzB/C0tvuHAnfWG8fMzLpVprBfDxwu6RmS9gZeB/xjb2KZmVlVHXfFRMQeSe8EvgIsAT4TETd3Me8Fu2v6qMnZoNn5nK26JudztuoWPV/HB0/NzOyJwVeempllxoXdzCwzLuxmZpl5QhZ2Saslre53jtlIeqak90g6od9ZZmpyNmh2Pmerrsn5mpwNqud7QhV2SUOSrgauAD4u6SX9ztRO0n8Evkb6Lp23SXp7nyP9QpOzQbPzOVt1Tc7X5GzQZb6IaPQfsG/b7dcCnyhuvxn4B2BNcV99yHYC8Izp+QMfAE4r7r8I+BIw3I98Tc7W9HzOlme+JmerO18jt9glHSDpryTdCnxC0mHFoN8BflTcHgN+CJwx/bRFzHeUpO8B/w34rKQTIrX2UcAhABFxHfBN4C2Lma/J2Zqez9nyzNfkbL3K18jCDpwI7EtasEeAD0haTtoteTVARDwMXAy8pLj/WK/CSDpU0gFtD50CXBIRLyV9wLxB0uHAhdP5CpcBR0vap1f5mpyt6fmcLc98Tc62WPn6VtiVLJX0Vklfl3SWpGcVg58NPBIRe4A/A34KnAZ8FXiKpF8rxrsV+LGk43qU8UhJXwa+AXxY0vTXFP8/YL/i9heBu4B1pE/UJ7ftYdxH+nbL5/0qZWt6PmfLM1+Tsy12vr4V9mJX4zeBNwEfA/YB/qYYfBdwd/HJ9GPSwjyL1ADf5/GvBV4G3Fs8XgtJK9ruPh+4IyKGgCuBTxSP3wc8LGn/iLgP+FfgqUWObwLvLcbbG/g5sD33bE3P52x55mtytn7mW7TCLuk4SR+VdHpxX8CRwBUR8aWI+BhwmKQXAztIn2BHFk/fBgwUj/0F8CpJryZ9KAwC3+0y25MknS/pemCjpIOKfGuAayUpIv4RuF/SOtKewv7FcIr7BwOPkfYwDpb0N8BFwJ6IuDvHbE3P52zVNTlfk7M1Jd+iFHZJzwH+EngA+D1J7y3mvRp4oFhogPOBN5AK9R7gxcXjN5COGD8UEdcAG4DTgeOB/x4Rj7VNo4qXFvN7FemgxNnAAaQvOzuk2LsAuKDI9+1iWV4JEBHfKqaxNCK2AWcCNwP/MyLeQneanK3p+Zwtz3xNztaMfHOdLlP1j7RlfQZpt2Np8difAmcVt1vAJ4GTgZcDX2l77tNIuyqQCvmNpF9hOgb438BT2sYtfTpS0bBnAleTunNWFY9/EXh3cfsZwMZi+LGk/rAlbcv2k2I6q0l7Eu8EPgt8CljRRbs1NlvT8zmbX1e33b//q3WLXdLzSQc4fxv4IPC+YtAO0m+mQvrkuRb4XeCfgUMkPVfSskj96TskvSQiriR93eVHgUuBiyLi36bnFUXLlHQS8BrgQ8BxpL59SGfbTO8d/Bj4OvDKiLie9Ik7UsxzCrgOODYidgBvJHUF3QW8LyIeLBuobU/j1U3LNoPbrprGtRu47brJ9kRouzI/jfdLJL0QOBz4akT8hLQ1fmtEnC7pBcA5klrAOPCfJO0XEQ9J+i7we6RzNC8E/gD4pKTdwCRwezGLvwIujIhdJTIpIkLSsaTdnK8DWyKdHvnrwG0RcaWk20lXr74C2Ar8jqRVEXGPpH8FHpT0dODPgdMkHUz61ah7SbtORMQEMFGh3VqkvZoHgI8DdwPP7Hc2t121bE+EdnPb5dd28ym1xa5kmaQ3SbqR1LF/IDBdeH8ObC+2vm8g7VocBzzE46fwADxK2gU5hLRVfhOpf/1q4J6IuAPSVnnFov5S4DOko8ovBz5SjPIYcKuk5RFxe5HvuaQX607S+aTTy7GE1D6XFBlPBdYCm6PiOa6SVkr6bDHN24FzI+JuSXuRPsn7mW1J0Xa/SdoVbEzbFevdgKTzaVjbFfMMScM0c53bR9KKhrbdAQ1vuwFJ+0q6gIa13YI66a8BVgAvLm4fWAT75CzjnUW6DHZ1cf9kUn/6YaSvALi6eHxfUjfMqrbnHgPs3UmeGfPcD3gbj2/5LwP+CHhHMfxJwPeK6Z9C6u8aKoadVCzLquL2JLCS1L//5fY8wF5dZLuIdMXYAKlr6cy2caaPQ7wT+B+Lla3tdT2DtLKtJx3gaUrbTWe7tFivDmpY2+0PbCH9khjAe5rQbjPyfRn46+L+x4C39bvtSO+JN5Pe/5c0re3a8l0J/H3xWGPWu07/Ftxil3Q2cBuwRdJgRNxP6he6s+gbf40ev0DoW6QDoNMXGl1LOoj6UERcAPxU0udIB0VvAX7RhxQRN0bEIwvlmZHtEOByYBj4HOkAxWtJewl7iun+lHTg9d2kvq+Defw0ymtI59I/EhGXA58mXc16HumI9aNt+Up9qs7I9nfA24tstwJHSNpYbEX9vtIFV1eQ9mB6nq3It4L05jqBdP3AK0jHPY4lbSn1s+3as20mnS3wWtI1DL/R77YrLCdde/EsSatI6/ySYpp9abdZ8u1NWteeSuriOFrSR/rVdpKWkY6xnQx8PCJ+txh0TNs0+9Z2M/J9LCKmt7gngaP62XaldfAJNkzavfhb4D3FY8eSitYdRfALgU3FsHOAD7c9/3rgmOL2PqRTgI6t41OJtPK+qO3+6aQtkzcD3257/KnAncXtd5Au231S8fwvAU9vG3dVj7K9iXSk+9eBvy/+Xg/8L9K5/IuWrW16B7bd/s+kN9Op/W67WbL9MemUsWc2qO3eTOprfT/wVtKBtOv73W6z5HsfaY9nVRPajrQHduqMx04BrmtC282R7+lFhr6vdx0vRwcLOn1qzinAeHF7GWlramVx/zDS1vqxpF3Ai0lbWv9E+qTapyfhUx+X4Be/3foCHu/uuZd0zuj0uF+jKLSk3aevFuP8ySJlOwb4xvSK2zbeMtLB5ROK++f0OtuMnAeQjm/sBD5c3L8XGOxX282S7a5ivisouvn61XZtr+dbSN1srwW+UDx2T7/bbY58Y8Vj7acL92W9I3VR3ApsKub/gaJ+3Acc3IB1rj3fVaQv5jq03+td6eUoscBPJl0o9Jzi/tIZw88HTp5egUhdD2fSo6I+x8p8AY+fL/854KPF7V8j7XE8ve2FOZq2rwRepGzvaH+suH1I0XbPXexsbRn+kHS+7WZSv/Y3izec+tl2M7KdRzqt7NlNaDvSV0YvIfWhXk3aMr4JeH+/17lZ8v0z6QyzFzSk7b5C2gN7Gmkr+CzShmFT1rn2fF8gXfp/eBPartO/6aLTEUmfAn4WERuK+3uRzrt8B/Ac4JQo2U9eF0mHkvq03hURtyp9odhokWs18J2o56qybrK9PSJuKx47htQtta7I9of9yNZO6TqEM0lvsiNJK+uh9LHt2rIdTXqz/TnpLKuT6FPbSRogdXPsQ2qn3yBdeHI2aUv5cPrYbrPkO5x0fOK3SMe8XkZqv76sdypOey5uP4/0Pr2WdEl939e5GfmOJl3pfi7pm2b7tt6VUfY89s3AucVBhiNJK/HxpBfl7H4V9cIxFOfASzqD1P9/NqkL6QeRTr/sd7YfFdluJ60ce0hb8Tf2MVu7e0kHAd8XEX8n6TTg5obku5/UT3wT6XVdRv/abg/p7IlHSVvqPyet/5PAexvQbnPle1jSa0gFv2/r3XTRLNxPOu70/oi4sAFtNzPfA6SN123Af6W/613Hym6xv450oPRh0jeOXRkRt/QoWymSriUdXNtOOof0QxHxvb6GKszIdhewoUHttpK0BfcG0vffbwbOi4hH533iIpgl26cjYlN/U/2y4sKT6b7su/qdZ6Yi38nAZyOdddLvPPuQfnPhjaQ96r8EPhXpa7r7bpZ8myPiz/qbqpyOC7uk55LO57yYdLCotq/K7VaxB/FB0pbw5yNdtdYITc4GIGkpqfvlYVK+Jr2ujc0G6aIu4LEos3W0iJqcT9KZpNNqP9e01xWan28hpbbYzcys+Zr603hmZlaRC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDP/H+KofDj+oV4qAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.hist(xrot=20)" ] @@ -341,9 +2479,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": true - }, + "metadata": {}, "outputs": [], "source": [] } @@ -364,7 +2500,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4,