diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5b58a3ae26346460dd39e62a39c55243d7..3620489168ed14771bf2cb969e74064e80ef9c9d 100644 --- a/module3/exo1/analyse-syndrome-grippal.ipynb +++ b/module3/exo1/analyse-syndrome-grippal.ipynb @@ -9,7 +9,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -28,10 +28,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 2, + "metadata": {}, "outputs": [], "source": [ "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" @@ -61,9 +59,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
02021183126509582.015718.01914.024.0FRFrance
12021173120588891.015225.01813.023.0FRFrance
220211631650512735.020275.02519.031.0FRFrance
320211531930615398.023214.02923.035.0FRFrance
420211432107317099.025047.03226.038.0FRFrance
520211332641322094.030732.04033.047.0FRFrance
620211233065825919.035397.04639.053.0FRFrance
720211132498820718.029258.03832.044.0FRFrance
820211031953915951.023127.03025.035.0FRFrance
920210931757213926.021218.02721.033.0FRFrance
1020210832088216907.024857.03226.038.0FRFrance
1120210732239318303.026483.03428.040.0FRFrance
1220210632318319134.027232.03529.041.0FRFrance
1320210532242618445.026407.03428.040.0FRFrance
1420210432580421491.030117.03932.046.0FRFrance
1520210332181017894.025726.03327.039.0FRFrance
1620210231732013906.020734.02621.031.0FRFrance
1720210132179917778.025820.03327.039.0FRFrance
1820205332122016498.025942.03225.039.0FRFrance
1920205231642812285.020571.02519.031.0FRFrance
2020205132161917370.025868.03327.039.0FRFrance
2120205031684513220.020470.02620.032.0FRFrance
222020493129399923.015955.02015.025.0FRFrance
2320204831380410641.016967.02116.026.0FRFrance
2420204731908515285.022885.02923.035.0FRFrance
2520204632480120503.029099.03831.045.0FRFrance
2620204534251636857.048175.06556.074.0FRFrance
2720204434456738521.050613.06859.077.0FRFrance
2820204334373737523.049951.06657.075.0FRFrance
2920204233514529812.040478.05345.061.0FRFrance
.................................
187619852132609619621.032571.04735.059.0FRFrance
187719852032789620885.034907.05138.064.0FRFrance
187819851934315432821.053487.07859.097.0FRFrance
187919851834055529935.051175.07455.093.0FRFrance
188019851733405324366.043740.06244.080.0FRFrance
188119851635036236451.064273.09166.0116.0FRFrance
188219851536388145538.082224.011683.0149.0FRFrance
18831985143134545114400.0154690.0244207.0281.0FRFrance
18841985133197206176080.0218332.0357319.0395.0FRFrance
18851985123245240223304.0267176.0445405.0485.0FRFrance
18861985113276205252399.0300011.0501458.0544.0FRFrance
18871985103353231326279.0380183.0640591.0689.0FRFrance
18881985093369895341109.0398681.0670618.0722.0FRFrance
18891985083389886359529.0420243.0707652.0762.0FRFrance
18901985073471852432599.0511105.0855784.0926.0FRFrance
18911985063565825518011.0613639.01026939.01113.0FRFrance
18921985053637302592795.0681809.011551074.01236.0FRFrance
18931985043424937390794.0459080.0770708.0832.0FRFrance
18941985033213901174689.0253113.0388317.0459.0FRFrance
189519850239758680949.0114223.0177147.0207.0FRFrance
189619850138548965918.0105060.0155120.0190.0FRFrance
189719845238483060602.0109058.0154110.0198.0FRFrance
1898198451310172680242.0123210.0185146.0224.0FRFrance
18991984503123680101401.0145959.0225184.0266.0FRFrance
1900198449310107381684.0120462.0184149.0219.0FRFrance
190119844837862060634.096606.0143110.0176.0FRFrance
190219844737202954274.089784.013199.0163.0FRFrance
190319844638733067686.0106974.0159123.0195.0FRFrance
19041984453135223101414.0169032.0246184.0308.0FRFrance
190519844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1906 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202118 3 12650 9582.0 15718.0 19 14.0 \n", + "1 202117 3 12058 8891.0 15225.0 18 13.0 \n", + "2 202116 3 16505 12735.0 20275.0 25 19.0 \n", + "3 202115 3 19306 15398.0 23214.0 29 23.0 \n", + "4 202114 3 21073 17099.0 25047.0 32 26.0 \n", + "5 202113 3 26413 22094.0 30732.0 40 33.0 \n", + "6 202112 3 30658 25919.0 35397.0 46 39.0 \n", + "7 202111 3 24988 20718.0 29258.0 38 32.0 \n", + "8 202110 3 19539 15951.0 23127.0 30 25.0 \n", + "9 202109 3 17572 13926.0 21218.0 27 21.0 \n", + "10 202108 3 20882 16907.0 24857.0 32 26.0 \n", + "11 202107 3 22393 18303.0 26483.0 34 28.0 \n", + "12 202106 3 23183 19134.0 27232.0 35 29.0 \n", + "13 202105 3 22426 18445.0 26407.0 34 28.0 \n", + "14 202104 3 25804 21491.0 30117.0 39 32.0 \n", + "15 202103 3 21810 17894.0 25726.0 33 27.0 \n", + "16 202102 3 17320 13906.0 20734.0 26 21.0 \n", + "17 202101 3 21799 17778.0 25820.0 33 27.0 \n", + "18 202053 3 21220 16498.0 25942.0 32 25.0 \n", + "19 202052 3 16428 12285.0 20571.0 25 19.0 \n", + "20 202051 3 21619 17370.0 25868.0 33 27.0 \n", + "21 202050 3 16845 13220.0 20470.0 26 20.0 \n", + "22 202049 3 12939 9923.0 15955.0 20 15.0 \n", + "23 202048 3 13804 10641.0 16967.0 21 16.0 \n", + "24 202047 3 19085 15285.0 22885.0 29 23.0 \n", + "25 202046 3 24801 20503.0 29099.0 38 31.0 \n", + "26 202045 3 42516 36857.0 48175.0 65 56.0 \n", + "27 202044 3 44567 38521.0 50613.0 68 59.0 \n", + "28 202043 3 43737 37523.0 49951.0 66 57.0 \n", + "29 202042 3 35145 29812.0 40478.0 53 45.0 \n", + "... ... ... ... ... ... ... ... \n", + "1876 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1877 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1878 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1879 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1880 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1881 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1882 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1883 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1884 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1885 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1886 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1887 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1888 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1889 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1890 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1891 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1892 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1893 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1894 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1895 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1896 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1897 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1898 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1899 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1900 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1901 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1902 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1903 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1904 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1905 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 24.0 FR France \n", + "1 23.0 FR France \n", + "2 31.0 FR France \n", + "3 35.0 FR France \n", + "4 38.0 FR France \n", + "5 47.0 FR France \n", + "6 53.0 FR France \n", + "7 44.0 FR France \n", + "8 35.0 FR France \n", + "9 33.0 FR France \n", + "10 38.0 FR France \n", + "11 40.0 FR France \n", + "12 41.0 FR France \n", + "13 40.0 FR France \n", + "14 46.0 FR France \n", + "15 39.0 FR France \n", + "16 31.0 FR France \n", + "17 39.0 FR France \n", + "18 39.0 FR France \n", + "19 31.0 FR France \n", + "20 39.0 FR France \n", + "21 32.0 FR France \n", + "22 25.0 FR France \n", + "23 26.0 FR France \n", + "24 35.0 FR France \n", + "25 45.0 FR France \n", + "26 74.0 FR France \n", + "27 77.0 FR France \n", + "28 75.0 FR France \n", + "29 61.0 FR France \n", + "... ... ... ... \n", + "1876 59.0 FR France \n", + "1877 64.0 FR France \n", + "1878 97.0 FR France \n", + "1879 93.0 FR France \n", + "1880 80.0 FR France \n", + "1881 116.0 FR France \n", + "1882 149.0 FR France \n", + "1883 281.0 FR France \n", + "1884 395.0 FR France \n", + "1885 485.0 FR France \n", + "1886 544.0 FR France \n", + "1887 689.0 FR France \n", + "1888 722.0 FR France \n", + "1889 762.0 FR France \n", + "1890 926.0 FR France \n", + "1891 1113.0 FR France \n", + "1892 1236.0 FR France \n", + "1893 832.0 FR France \n", + "1894 459.0 FR France \n", + "1895 207.0 FR France \n", + "1896 190.0 FR France \n", + "1897 198.0 FR France \n", + "1898 224.0 FR France \n", + "1899 266.0 FR France \n", + "1900 219.0 FR France \n", + "1901 176.0 FR France \n", + "1902 163.0 FR France \n", + "1903 195.0 FR France \n", + "1904 308.0 FR France \n", + "1905 213.0 FR France \n", + "\n", + "[1906 rows x 10 columns]" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data = pd.read_csv(data_url, skiprows=1)\n", "raw_data" @@ -78,9 +1043,73 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
166919891930NaNNaN0NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1669 198919 3 0 NaN NaN 0 NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1669 FR France " + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -94,9 +1123,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 5, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
02021183126509582.015718.01914.024.0FRFrance
12021173120588891.015225.01813.023.0FRFrance
220211631650512735.020275.02519.031.0FRFrance
320211531930615398.023214.02923.035.0FRFrance
420211432107317099.025047.03226.038.0FRFrance
520211332641322094.030732.04033.047.0FRFrance
620211233065825919.035397.04639.053.0FRFrance
720211132498820718.029258.03832.044.0FRFrance
820211031953915951.023127.03025.035.0FRFrance
920210931757213926.021218.02721.033.0FRFrance
1020210832088216907.024857.03226.038.0FRFrance
1120210732239318303.026483.03428.040.0FRFrance
1220210632318319134.027232.03529.041.0FRFrance
1320210532242618445.026407.03428.040.0FRFrance
1420210432580421491.030117.03932.046.0FRFrance
1520210332181017894.025726.03327.039.0FRFrance
1620210231732013906.020734.02621.031.0FRFrance
1720210132179917778.025820.03327.039.0FRFrance
1820205332122016498.025942.03225.039.0FRFrance
1920205231642812285.020571.02519.031.0FRFrance
2020205132161917370.025868.03327.039.0FRFrance
2120205031684513220.020470.02620.032.0FRFrance
222020493129399923.015955.02015.025.0FRFrance
2320204831380410641.016967.02116.026.0FRFrance
2420204731908515285.022885.02923.035.0FRFrance
2520204632480120503.029099.03831.045.0FRFrance
2620204534251636857.048175.06556.074.0FRFrance
2720204434456738521.050613.06859.077.0FRFrance
2820204334373737523.049951.06657.075.0FRFrance
2920204233514529812.040478.05345.061.0FRFrance
.................................
187619852132609619621.032571.04735.059.0FRFrance
187719852032789620885.034907.05138.064.0FRFrance
187819851934315432821.053487.07859.097.0FRFrance
187919851834055529935.051175.07455.093.0FRFrance
188019851733405324366.043740.06244.080.0FRFrance
188119851635036236451.064273.09166.0116.0FRFrance
188219851536388145538.082224.011683.0149.0FRFrance
18831985143134545114400.0154690.0244207.0281.0FRFrance
18841985133197206176080.0218332.0357319.0395.0FRFrance
18851985123245240223304.0267176.0445405.0485.0FRFrance
18861985113276205252399.0300011.0501458.0544.0FRFrance
18871985103353231326279.0380183.0640591.0689.0FRFrance
18881985093369895341109.0398681.0670618.0722.0FRFrance
18891985083389886359529.0420243.0707652.0762.0FRFrance
18901985073471852432599.0511105.0855784.0926.0FRFrance
18911985063565825518011.0613639.01026939.01113.0FRFrance
18921985053637302592795.0681809.011551074.01236.0FRFrance
18931985043424937390794.0459080.0770708.0832.0FRFrance
18941985033213901174689.0253113.0388317.0459.0FRFrance
189519850239758680949.0114223.0177147.0207.0FRFrance
189619850138548965918.0105060.0155120.0190.0FRFrance
189719845238483060602.0109058.0154110.0198.0FRFrance
1898198451310172680242.0123210.0185146.0224.0FRFrance
18991984503123680101401.0145959.0225184.0266.0FRFrance
1900198449310107381684.0120462.0184149.0219.0FRFrance
190119844837862060634.096606.0143110.0176.0FRFrance
190219844737202954274.089784.013199.0163.0FRFrance
190319844638733067686.0106974.0159123.0195.0FRFrance
19041984453135223101414.0169032.0246184.0308.0FRFrance
190519844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1905 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202118 3 12650 9582.0 15718.0 19 14.0 \n", + "1 202117 3 12058 8891.0 15225.0 18 13.0 \n", + "2 202116 3 16505 12735.0 20275.0 25 19.0 \n", + "3 202115 3 19306 15398.0 23214.0 29 23.0 \n", + "4 202114 3 21073 17099.0 25047.0 32 26.0 \n", + "5 202113 3 26413 22094.0 30732.0 40 33.0 \n", + "6 202112 3 30658 25919.0 35397.0 46 39.0 \n", + "7 202111 3 24988 20718.0 29258.0 38 32.0 \n", + "8 202110 3 19539 15951.0 23127.0 30 25.0 \n", + "9 202109 3 17572 13926.0 21218.0 27 21.0 \n", + "10 202108 3 20882 16907.0 24857.0 32 26.0 \n", + "11 202107 3 22393 18303.0 26483.0 34 28.0 \n", + "12 202106 3 23183 19134.0 27232.0 35 29.0 \n", + "13 202105 3 22426 18445.0 26407.0 34 28.0 \n", + "14 202104 3 25804 21491.0 30117.0 39 32.0 \n", + "15 202103 3 21810 17894.0 25726.0 33 27.0 \n", + "16 202102 3 17320 13906.0 20734.0 26 21.0 \n", + "17 202101 3 21799 17778.0 25820.0 33 27.0 \n", + "18 202053 3 21220 16498.0 25942.0 32 25.0 \n", + "19 202052 3 16428 12285.0 20571.0 25 19.0 \n", + "20 202051 3 21619 17370.0 25868.0 33 27.0 \n", + "21 202050 3 16845 13220.0 20470.0 26 20.0 \n", + "22 202049 3 12939 9923.0 15955.0 20 15.0 \n", + "23 202048 3 13804 10641.0 16967.0 21 16.0 \n", + "24 202047 3 19085 15285.0 22885.0 29 23.0 \n", + "25 202046 3 24801 20503.0 29099.0 38 31.0 \n", + "26 202045 3 42516 36857.0 48175.0 65 56.0 \n", + "27 202044 3 44567 38521.0 50613.0 68 59.0 \n", + "28 202043 3 43737 37523.0 49951.0 66 57.0 \n", + "29 202042 3 35145 29812.0 40478.0 53 45.0 \n", + "... ... ... ... ... ... ... ... \n", + "1876 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1877 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1878 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1879 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1880 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1881 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1882 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1883 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1884 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1885 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1886 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1887 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1888 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1889 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1890 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1891 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1892 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1893 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1894 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1895 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1896 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1897 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1898 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1899 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1900 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1901 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1902 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1903 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1904 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1905 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 24.0 FR France \n", + "1 23.0 FR France \n", + "2 31.0 FR France \n", + "3 35.0 FR France \n", + "4 38.0 FR France \n", + "5 47.0 FR France \n", + "6 53.0 FR France \n", + "7 44.0 FR France \n", + "8 35.0 FR France \n", + "9 33.0 FR France \n", + "10 38.0 FR France \n", + "11 40.0 FR France \n", + "12 41.0 FR France \n", + "13 40.0 FR France \n", + "14 46.0 FR France \n", + "15 39.0 FR France \n", + "16 31.0 FR France \n", + "17 39.0 FR France \n", + "18 39.0 FR France \n", + "19 31.0 FR France \n", + "20 39.0 FR France \n", + "21 32.0 FR France \n", + "22 25.0 FR France \n", + "23 26.0 FR France \n", + "24 35.0 FR France \n", + "25 45.0 FR France \n", + "26 74.0 FR France \n", + "27 77.0 FR France \n", + "28 75.0 FR France \n", + "29 61.0 FR France \n", + "... ... ... ... \n", + "1876 59.0 FR France \n", + "1877 64.0 FR France \n", + "1878 97.0 FR France \n", + "1879 93.0 FR France \n", + "1880 80.0 FR France \n", + "1881 116.0 FR France \n", + "1882 149.0 FR France \n", + "1883 281.0 FR France \n", + "1884 395.0 FR France \n", + "1885 485.0 FR France \n", + "1886 544.0 FR France \n", + "1887 689.0 FR France \n", + "1888 722.0 FR France \n", + "1889 762.0 FR France \n", + "1890 926.0 FR France \n", + "1891 1113.0 FR France \n", + "1892 1236.0 FR France \n", + "1893 832.0 FR France \n", + "1894 459.0 FR France \n", + "1895 207.0 FR France \n", + "1896 190.0 FR France \n", + "1897 198.0 FR France \n", + "1898 224.0 FR France \n", + "1899 266.0 FR France \n", + "1900 219.0 FR France \n", + "1901 176.0 FR France \n", + "1902 163.0 FR France \n", + "1903 195.0 FR France \n", + "1904 308.0 FR France \n", + "1905 213.0 FR France \n", + "\n", + "[1905 rows x 10 columns]" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -122,7 +2118,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 6, "metadata": {}, "outputs": [], "source": [ @@ -152,10 +2148,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 7, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -179,9 +2173,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n" + ] + } + ], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", @@ -199,9 +2201,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXm4HUWZ/7/v2e6a3Gw3ISSBJBDBsAmERUREkE1QHJQRdSSjzOAIo46zKPgbBxVQcBwYcARHkdUFEGEA2YxhkR1CCIQkhIQkZF9vkntz93NO/f7oqj7V51R19zmnz3bzfp7nPqdvL1XV3VX11rtUNQkhwDAMwzDlEqt1ARiGYZiRAQsUhmEYJhJYoDAMwzCRwAKFYRiGiQQWKAzDMEwksEBhGIZhIoEFCsMwDBMJLFAYhmGYSGCBwjAMw0RCotYFqCYTJkwQ06dPr3UxGIZhGorXXnttuxCiM+i8vUqgTJ8+HQsWLKh1MRiGYRoKInovzHls8mIYhmEigQUKwzAMEwksUBiGYZhIYIHCMAzDRAILFIZhGCYSAgUKEd1KRFuJ6C1t3zgimkdEK+TvWO3Y5US0koiWE9EZ2v6jiWixPHYjEZHc30RE98j9LxPRdO2auTKPFUQ0V9s/Q567Ql6bKv9RMAzDMOUQRkO5HcCZefsuAzBfCDELwHz5P4hoNoALABwir7mJiOLympsBXAxglvxTaV4EYKcQ4kAA1wO4VqY1DsAVAI4DcCyAKzTBdS2A62X+O2UaDMMwTA0JFChCiL8A6MrbfS6AO+T2HQA+pe2/WwgxKIRYDWAlgGOJaDKA0UKIF4XzzeE7865Rad0H4FSpvZwBYJ4QoksIsRPAPABnymOnyHPz8x/R9A2lcf/C9eDPNjMMU4+UOrFxkhBiEwAIITYR0US5fwqAl7Tz1st9w3I7f7+6Zp1MK01EuwGM1/fnXTMewC4hRNqQ1ojmew8twb0L1mO/ca2YM31crYvDMAzjIWqnPBn2CZ/9pVzjl1ZhgYguJqIFRLRg27ZtttMagvU7+wEAA8PZGpeEYRimkFIFyhZpxoL83Sr3rwcwTTtvKoCNcv9Uw37PNUSUANABx8RmS2s7gDHy3Py0ChBC/EIIMUcIMaezM3ApmromnXXkZiJukqkMwzC1pVSB8hAAFXU1F8CD2v4LZOTWDDjO91ekeayHiI6XPpAL865RaX0GwJPSz/IEgNOJaKx0xp8O4Al57Cl5bn7+I5qMEigxFigMw9QfgT4UIvodgJMBTCCi9XAir64BcC8RXQRgLYDzAUAIsYSI7gWwFEAawKVCiIxM6qtwIsZaADwm/wDgVwDuIqKVcDSTC2RaXUR0JYBX5Xk/EEKo4IBvA7ibiK4C8LpMY8STzjimrjgLFIZh6pBAgSKE+Jzl0KmW868GcLVh/wIAhxr2D0AKJMOxWwHcati/Ck4o8V6Fa/KK8XxUhmHqD+6ZGghl8mJ5wjBMPcJdUwOhBArDMEw9wgKlgSDpOuF5jQzD1CMsUBoIklNwsixRGIapQ1igNBCsoTAMU8+wQGlAWENhGKYeYYHSQMgV/+3rzDAMw9QQFigNRMw1ebFIYRim/mCB0kAoHwpHDzMMU4+wQGlAWEFhGKYeYYHSQHDYMMMw9QwLlAaCw4YZhqlnWKA0EG6UF0sUhmHqEBYoDYRatJ6d8gzD1CMsUBoI1+RV5zNRPvzjJ/Hl218NPpFhmBFF4PdQmPqhUTSUdV39WNfVX+tiMAxTZVhDaSCUD4WjvBiGqUdYoDQQ7od/WZ4wDFOHsEBpIHIz5VmiMAxTf7BAaSByExtrXBCGYRgDLFAYhmGYSGCB0oDwxEaGYeoRFiiNhDsPhWEYpv5ggdJIsCRhGKaOYYHSgLDFi2GYeoQFCsMwDBMJLFAaElZRGIapP1igNBL8PRSGYeoYFiiNBAsShmHqGBYoDYifXPnUz57H9fPeqVpZGIZhFGUJFCL6JhEtIaK3iOh3RNRMROOIaB4RrZC/Y7XzLyeilUS0nIjO0PYfTUSL5bEbSS6rS0RNRHSP3P8yEU3Xrpkr81hBRHPLuY+RxKJ1u3DD/BW1LgbDMHshJQsUIpoC4OsA5gghDgUQB3ABgMsAzBdCzAIwX/4PIpotjx8C4EwANxFRXCZ3M4CLAcySf2fK/RcB2CmEOBDA9QCulWmNA3AFgOMAHAvgCl1wjXTYh8IwTD1SrskrAaCFiBIAWgFsBHAugDvk8TsAfEpunwvgbiHEoBBiNYCVAI4loskARgshXhTOmiJ35l2j0roPwKlSezkDwDwhRJcQYieAecgJoZFLg3yxkWGYvZOSBYoQYgOAnwBYC2ATgN1CiD8BmCSE2CTP2QRgorxkCoB1WhLr5b4pcjt/v+caIUQawG4A433SYhiGYWpEOSavsXA0iBkA9gXQRkR/43eJYZ/w2V/qNfnlvJiIFhDRgm3btvkUr3FgkxfDMPVIOSavjwFYLYTYJoQYBnA/gBMAbJFmLMjfrfL89QCmaddPhWMiWy+38/d7rpFmtQ4AXT5pFSCE+IUQYo4QYk5nZ2eJt1pfsDxhGKYeKUegrAVwPBG1Sr/GqQCWAXgIgIq6mgvgQbn9EIALZOTWDDjO91ekWayHiI6X6VyYd41K6zMAnpR+licAnE5EY6WmdLrcxzAMw9SIRKkXCiFeJqL7ACwEkAbwOoBfAGgHcC8RXQRH6Jwvz19CRPcCWCrPv1QIkZHJfRXA7QBaADwm/wDgVwDuIqKVcDSTC2RaXUR0JYBX5Xk/EEJ0lXovjQZ/D4VhmHqkZIECAEKIK+CE7+oMwtFWTOdfDeBqw/4FAA417B+AFEiGY7cCuLXIIjMMwzAVgmfKMwzDMJHAAqUBYYsXwzD1CAsUhmEYJhJYoDQgPFOeYZh6hAVKA7I3mLy+//AS3L9wffCJDMPUDWVFeTFMpbjt+TUAgPOOmup/IsMwdQNrKA3I3qChMAzTeLBAYRiGYSKBBUoDwgoKwzD1CAuUBkItscxLrzAMU4+wQGkgWIwwDFPPsEBpQFiwMAxTj7BAaSDcr4qxRGEYpg5hgcIwDMNEAguUBoSXXmEYph5hgdJAsBhhGKaeYYHSgHDUMMMw9QgLlAbCnYdS01IwDMOYYYHCMAzDRAILlAaETV4Mw9QjLFAYhmGYSGCB0oBw2DDDMPUIC5QGhE1eDMPUIyxQRhC8CjHDMLWEBUoDYhMbLE8YhqklLFAaEYvk2JvkycBwBjv2DNa6GAzDaLBAKYK3NuzGn5ZsrnUxrOxNJq/P//IlHH3Vn2tdDIZhNBK1LkAjcc5PnwMArLnm7JqWw2ryqmopasvCtbtqXQSGYfJgDYVhGIaJBBYoDYjNsrUXWbwYhqlDWKA0IDZfCU94ZBimlpQlUIhoDBHdR0RvE9EyIvogEY0jonlEtEL+jtXOv5yIVhLRciI6Q9t/NBEtlsduJCKS+5uI6B65/2Uimq5dM1fmsYKI5pZzHyOFkaKh7E3BBQwzkihXQ7kBwONCiIMBHAFgGYDLAMwXQswCMF/+DyKaDeACAIcAOBPATUQUl+ncDOBiALPk35ly/0UAdgohDgRwPYBrZVrjAFwB4DgAxwK4QhdcI52R3t2yPGGYxqRkgUJEowGcBOBXACCEGBJC7AJwLoA75Gl3APiU3D4XwN1CiEEhxGoAKwEcS0STAYwWQrwonKHpnXnXqLTuA3Cq1F7OADBPCNElhNgJYB5yQmjEQvKDKCPdhzJCboNh9jrK0VBmAtgG4DYiep2IbiGiNgCThBCbAED+TpTnTwGwTrt+vdw3RW7n7/dcI4RIA9gNYLxPWns1I8WHwiYvhmlMyhEoCQBHAbhZCHEkgF5I85YFMuwTPvtLvcabKdHFRLSAiBZs27bNp3iNw0hfemWE3AbD7HWUI1DWA1gvhHhZ/n8fHAGzRZqxIH+3audP066fCmCj3D/VsN9zDRElAHQA6PJJqwAhxC+EEHOEEHM6OztLuM36YaQIjCD2lvtkmJFGyQJFCLEZwDoiOkjuOhXAUgAPAVBRV3MBPCi3HwJwgYzcmgHH+f6KNIv1ENHx0j9yYd41Kq3PAHhS+lmeAHA6EY2VzvjT5b69AnvY8MhgpJjuGGZvo9ylV74G4DdElAKwCsCX4Aipe4noIgBrAZwPAEKIJUR0LxyhkwZwqRAiI9P5KoDbAbQAeEz+AY7D/y4iWglHM7lAptVFRFcCeFWe9wMhRFeZ9xKabFYgFjNZ3SoLBWQ5UnwPI+Q2GGavoyyBIoRYBGCO4dCplvOvBnC1Yf8CAIca9g9ACiTDsVsB3FpMeaMiIwRiRjdObRkp/TALFIZpTHimfAlksrXt8UZ+2PAIuZEqcvGdC/Cha56sdTGYvRxebbgEsjXquQOzHSH98EgRjNXkT0u31LoIDMMaSinUXEOxSI6RMrIfGXcRDVu6B/DaeztrXQyGCQULlBLIZmuT714zU76IGxkpgQg2PnbdM/j0zS/UuhgMEwoWKCWQrpVECWCkdK3F3McIlyfoGUjXugihWbm1Byu37ql1MZgawgKlBDI17sVGeB9alJCopD9rYDiDb96zCJt291csj5HEx677Cz523TO1LgZTQ1iglECtFRS7yWtkiJqiTF4VLMdTb2/FA69vwPceWlLBXMKRrbHfjmHCwAKlBGqtodioz1IVTzGPt5KvoinpNI+B4dqbOIdrPYrZi7jrxTWYftkj6BkYrnVRGg4WKCVQ69GiNcprhEiUonwoZYrRgeEMtnQPGI81J5zP9QymM8bj1WQ4M0JebgNw10vvAQA27GJTZ7GwQCmBWocN26iHsOEozG7FRXmVl9c//Po1HPfD+cZj9aShDKVrX4a9hdaUMz2vd7D2A4lGgwVKCdTa5GXNvvbyJBKqeRtPL7d/0iAVdzSUeujM6zWycCTSmnLee99Q40TY1QssUEqg1iYvG/VQqihkbb34UBT1oJGOFHNmIxCXC7/WwWtvOFiglECtNRQbdVqsoinGdFfJsGGVdrnawQsrt6N3sLzRbj0Itb2NkRI1WU1YoJRArRu3/XsotW8AUZSgKA0lgvyC0k6X8b437x7A5295Gf9876KyylKNOvf4W5uxs3eo4vnUOxT0nQjGCguUEqjZ0ityyfx6HjhF45Svbn62dNS+dBkRVgPDjmP37c09JafhlKWsywPZvmcQ//Dr1/CVu16rbEYNRB03s7qFBUoJ+JlAbnl2FVZuLa/zsBGkgdSzoCmGYjStqG7ZpACofeVoB8oeX66GUWkzqwo8WNvVV9F8Ks2W7gEcdsUTWLapu+Q0WD8pHRYoJWCz22eyAlc9sgyf+lllF/Or5yCvqpu8Irpps6ajfCilZ6KsJ+UGclTa5BWTBa0Hs2k5PPX2VvQMpnHb86trXZS9EhYoJZCxKChKc9lTpgO2VOrBiRhJlFfFTi4uGeFqKKXbOF0NpcwHU+lv8AStZN0oNCfVZNQI7NIN/ixqAQuUErCNFqvlW+Hl67VzI2r1pg5bveYocrANQsJSrY+6VSqXnoFhrN1ReXNaU0JNRi19UqIrXFmiFA0LlBKwmrwqPYpE/ZsloihbcasNl52dNU8l2Mp5ra5QKrNuVNrkpfwGlarC5/7seZz0n09VJnENpaF095dvJRgpA7RqwgKlCKT1wtq4ax9OXNPsI6MWUV7GtKNIQ5avXA2j4tqv64muzPNcta23Iunm09bkLJsSRXABT/0pHhYoRaAclzZNpFoCxWryqgPNJRofSi2ivEwmL6WhlJ5Lzg9T3z4UuJpUZbOpNCKCd6Zka7XMjCMJFihFoASKLWKn0gJFdbTWKK8RUv9rE+Vl2hlduvUeNjxCqk6k91EPQS6NBguUIqAAkxePaLyU2iCLuSoqraxC8iSn5ZSbTkiBVOozj6qctSaK+1Az5dnkVTwsUIrA1VAsjbac+QphcLO1Lr1Se/SilSpfi+oUI7ppX5NXGelWcuJlOeflo26/4UflERafB4jFwwKlCHJOefPxWq9CXA+dga4xlFqa4jSUaDBHeZWfieqUyp19HbZzK7UTHDkaivMbRVNgDaV4WKAUgdJQbEuvVN6H4v21Ha8XSjZ5FXFdZKNIk0CJItmIOuqwt1nq86jWWKTSg56cn7F8p3w9DNAaDRYoReAuo1GjeShBkTj1UP89Jq8I0ojyXN90DKWNxORV5XdSupmxvOvDUulRf7QaSh00qAaDBUoR5Bb6Mx+v9TyUetBR9BKU3LlV6Fw/jK8uAr9CVB1c2BF3yRqKTL/SnWjFtfgIy88fySweFihFUD9hwxanfO3liYdSzQ71snx9FJ1rtecGleuUrzSVFlhRJB9kiWDslC1QiChORK8T0R/l/+OIaB4RrZC/Y7VzLyeilUS0nIjO0PYfTUSL5bEbScbtEVETEd0j979MRNO1a+bKPFYQ0dxy7yPkvQKo/cRGG17toDZl0fMtXUMJf2EFXSg5M1AZ6UY2yq2wDyUbxc0Wk0+F04/G/xVBInsZUWgo3wCwTPv/MgDzhRCzAMyX/4OIZgO4AMAhAM4EcBMRxeU1NwO4GMAs+Xem3H8RgJ1CiAMBXA/gWpnWOABXADgOwLEArtAFV6UIWnqlWqOvMD6UmlvfyqAWDdn07qIoRrVHuaJEARZBQFsoKl0vo3zcrKEUT1kChYimAjgbwC3a7nMB3CG37wDwKW3/3UKIQSHEagArARxLRJMBjBZCvCic4e2dedeotO4DcKrUXs4AME8I0SWE2AlgHnJCqGLUeh6KIkwuNdNQPGUoMY0aOOVNDzW39EpEeZRB2CKUHuUlPL+VotJafDTvjCc2lkq5Gsp/A/gWAH1cNEkIsQkA5O9EuX8KgHXaeevlvilyO3+/5xohRBrAbgDjfdKqKEFf36v0PJSg1HVT0f0LN1S0LNYyeKK8yjS/hMmvkjPlXStQOU758sNY9bKEza/U9Cvdh1Y+bDg6WEMpnpIFChGdA2CrECLsR6hNc7uEz/5Sr/FmSnQxES0gogXbtm0LVVAbQUuvVCuCJYzJ61t/eLOiZQlDNdpjVI/c3HlE4JSvcp9UanauyatkrTLchdWL8io9n9zHxmorUH713Gqsa7BPMpejoXwIwCeJaA2AuwGcQkS/BrBFmrEgf7fK89cDmKZdPxXARrl/qmG/5xoiSgDoANDlk1YBQohfCCHmCCHmdHZ2lnankiCTV9VWG67nKC9h3CwuiZpEeRXuiyLkNyqzWaXDhsvVpMJmW/FFLiP1oUSXVrHs7hvGlX9cis/f8lLtClECJQsUIcTlQoipQojpcJztTwoh/gbAQwBU1NVcAA/K7YcAXCAjt2bAcb6/Is1iPUR0vPSPXJh3jUrrMzIPAeAJAKcT0VjpjD9d7qsoQUuvqMZC5a6zYaEYk1etiOQDWzVYvr5iUV7V1lCq4LcyXl+lfIIYKRMbVV/SM1Cbz4mXSqICaV4D4F4iugjAWgDnA4AQYgkR3QtgKYA0gEuFEOo7nV8FcDuAFgCPyT8A+BWAu4hoJRzN5AKZVhcRXQngVXneD4QQXRW4Fw/u91ACll6pkDzJEcLkVSrpTBYPvL4Bnz5qKmKx8u6k9KVXKnOuHyb/VzTrcEUzyKi0DyWMJrW7fxgA0NGSLDhWNyavCIYYue+hlJ1U2VS8L4mYSASKEOJpAE/L7R0ATrWcdzWAqw37FwA41LB/AFIgGY7dCuDWUstcChSgoahGGauQilKNAdOtz6/GDx99G1kh8Nlj9iv6+kiWXqng2cUQxQg1yhFzMfkVSxht7Ijv/wkAsOaaswuvD5lP5eehRJdWrX0ojQjPlC+CoC82pjOVNXkpbNU8ivq/q88ZhW7rGSw7rUo7eMvJI0w6UaRdyYmXJkqNNIzAlx2KSi9nEtVinEBtTV6NGmHGAqUI1CseTGeMx3MmkgppKOrXUtmiUPeTcadKDGVK7Jis/4SnmD6xnDvWn6Pf91DKodrLwpe7OkGlnfKNtfRK+WmVCguUvQDVAfUOmh1lrimsUhpKFSpZKuFUiWGbXS8Az9IrVTB6ldPwgsxzqkMpR+OMToOqbJRXtaLQKh3llfMFRTcYqAWNujAlC5QiUNWrb9CsoajvpFTakRZmHkqppKSGMpwuv0bX+0x5vcOo2OKQUYU1hzyvkk55P8JeV/GJjVFoKLIF11JJyEWMNpZbngVKMcgK1j8cYPKqcNiw1YcSQR7JuFP4oVI1lAjKU8x15QkU/zyjGeWWnURV8quSC8Ua0BIV0Xy22bm6ll9gzVYrYjRiWKAUgapetjW7VGOplA8liCg6wKDlZYLLUH55itJQIlgWxcnTpKGUnHQu3ai66AprAOWu5VU3PpQI02IfSvGwQCkC1dhsIxe1v8zpGz75e38LjkeQh1Kxo0irZA2lSlFeXuFXeDzKsOFqeeUrGTbse31YH0qNlycKA8F/RYxqUOtPYZQKC5QiCNJQ1P5a2T2jjHApN1qovDSqg0dDMR4vP4/ofCgVdsqr3wo//EZI341447DhomGBUgTqHVsXh6y4Oi88v6YzyiVnritfGlRlteGInPLG76FE4pQvO4miKNkpX6b0DJtt5aO8nN8oPttcW5OX89tgPnkWKMWgOsi0JabPdaTVyuQVoYYSSdhiFbzy5YzkPE55k8krgh6l2svX10or9ARj+BSiat9DKSMN4QoUNnkVCwuUIgjSUNIRRWYMDGewq2+o6OuiqILK/1PyBDfLdqlpRHHuT55Yjj++WbgYddDniqMxeZWfRjGU/QngEgn76edG+B5KlLPtS6VRBUolFoccsai2YPOhZCPyoXz+ly9h4dpdBWsmVaNzogi/VledeSjBJ//PUysBAOccvq9nv36PlZ4pXy7ho6hKzaDE6wyX+yVVte+hlBOsIX/rY+mVxrJ5sYZSAkE+lHJNXgvX7irpumiXnShRQ4nAh1Kt5euD7jEaB280hE2nXKd8qYgA4Zw7Zk/jT0s246gr52Fnb/HauakcpVIPn35uUAWFBUoxqNFP0BcbK7bySl45Co5H0AJc7apkW3x1HdnRzZT3P14qUYSxlpJfsQTda6A/KbRAsR+766X30NU7hHe29Pjn5UMUPhTXKe9zz129Qzjsiifw+tqdZeRkp1FNXixQikC94kCB4qOirNzagy3dA6XlH9BYoqiCZcZ4eSi1Ey2uI/c/169hBo2qo2jTUfUL4dfyKjV9/+PDAVEaYcPF/d6tWj+unEiwSLRKNd/MJ62XV+1Az2AaP3/m3aLTX729F2f+9198NbFKr7pRKVigFEGQDyWMhvKx6/6C4344P1R+xY5SomhMsZhKKwqTV2kUJU4CTh7yWZMseB5KdbWtKKiUEC+mLvol5ZdOQla+ckbnUSwOGSbKq5yO/udPv4u3N/fgiSWbrefUctmXcmCBUgRqFGar8FEvVW5b8dc+Uz4Ck1dIp/yqbXvw1obdhjJo21WIYQ0qp9+aZN6w4cKE6mnF2lr7UIYDPmcQ1nfmVzy17E86IK9H3twUyfd6bBQjlEp53GGWx6/0fJ1KwQKlCHIairmTUppLVJ3IoGV0bW2wkTiRwwnFU/7rGZzz0+f806qCHyao0af9BIrWok2NO5q1vKpL6cEU/tf5aXr5+fp2lD4HlUDx+3TCnsE0Lv3tQsy99RXfcpTz3HMait9ZpS9RRCFC83n5+r0A14diGUGpDioqh1p+OkF9RZR+j0p1TMWUIdS5Acdt5snCfEw+lOCC7Owdwl///EVs2t1vySMarTXs3I5KhWoHrT6d8ZTBp6P0Oaa+iPp/izbY85Ftb8Mu8/OOxu8VPDAsx+TlrpcXQkNpMBcKC5RicCc2WmqCEgBRaas2wWQ1eUXikFQbZV6PynVuOkG25uWb7RFDwVFewfnf99p6vLKmC7c8u9p4vBqWi6D5NGEoxxcFeOuq33PzK9+oZmda3KOL7b6FIKJsA373UU5H7wa+lCh46xkWKEXhr4FkQoxsiiE/nciWQvchJ08i0DQimG0fRJCt+UKLaQQIdsqHifYJWkwzuo7BT9MK15n7EVTOYgSK3wv0U3QmjmoCAJxy8ET79QH+jSjmkIRZHNIveSEEbnt+NXb3DRuPK03MV/CyU37kEzbKK6pOxGbysqUehRCIclJXNcKG/Uxawflo24Z0inLOWp59FIsVOtfbj0WioQQcL05DKW3knQlhMrb5L6MktzikTzmk6S1usH298O4OfP/hpbjqkaXGa93ljUI8Cw4bHsEE+VByAiWa/Gph8kKIxlRkUsVfF6HJyz8ffw0lzKqzQfbwaiyB7p0DYs6vdzCNB318E4FO+SAfStb+LL0alE8nKoIFipIntjNUfShncBVGM1WBA/F4YY+v5pnZvuxKYTSUxlRQeC2vYlAVzbqWl6qIEdWGQpOXP9HIk/I0FK8PpWSREvrMqDQUsw8lWOMMGkBGJU/8TSy5bdsA/kePLcOvX1qLfUY347iZ433TMFHMPJWCehtSg3J9lD7vNMjEGcXjDqOhKAFr0lC65ITF0S1J47VhFqNgH8pegKuh2FYblppLVDHk9oZl01wi0CpCjMp9r9dHyyHO/+Y9i/CDh72mgcBoNn3EG8EkuPw088shhP3ZUoD5oip+Ly0Lm4DdKe35Wy3zN4JKGfSc9Tqf3xnq//spOmFMXq51wHJKJD6UEGm4Gorh86wDUjNJxc3dq5rrFcrk1WBxXixQiiDnQzG3ioxWESsxKS6ookdo8YoktTCP4IHXN+DW51dbyhCcbjnCO3CmvG7GsWQTtFRNdEuv2I/p92Grm+0pxxixZzAdmIb5uF6WwnM9kxHzDuv/hvGv+PlJ3DZmOR7kZwyDujaUD8UgUIbksZjFARLmq6isoewF6LZVU6PSOyDTKKtYIZM/mgu8OoI6WKyG4j9XprQCBWooPvnbMLXtoBDnMM7uIB9KmIl2fUNp3PDnFb4T+vzQ07alkUyoWejFrb6QO+4vXP0mNnqO+flHlA/FpyyZoDXFIpAoWa2d2xj2MXmpAAZbnQnzzSEWKHsBQR2Z3hBMpodiR6vWKC+rhhKdVhRW+OV3YPpVJfthAu5DL1tYgWJq+EGfAA4z+zto1nOYZ/Dzp9/F9X9+B3e/stZ6TtjOZzhtPk/df9AWvz81AAAgAElEQVSyQfY8ctsmrVCv7/ll9fpQfPKQVclPaKjqZjcxqnzK0FyzwWkoH0rMpKFIgWLTtMI45UscW9QcFijFECAw9IZgqozFzqC3mbwqicj7DSJfoASZkcIQ9JiyAe/BhKnhB02dCBOdlJukZs43zDtT97DLMm8hCD0LWzSWuv8o5qkYha/PxEaPDyVElJffWl5BYcNhAimCcNtACJNXwlCv3FWTLQ87qM4AufJz2PAIJlBDCTB5FVvJbVFe5YyGAwmh7uvkLxoYNqLHvwgBGop2/2Gd8kEaitGEGWASA+C2eNvhMKVrScYB2MNMffOHt+w2k1duMl2pGoq/ycujofhEefnODg8xjyswbFh4f0vBNW37yC71nE0dvtJQrAIl4F04ebPJa8QjAkZaels2VcZiNZT884MiWML230II3PvqOvQaHLRhRmc6hTZ5vbMPV55iKcUpbxpJFucX8M/HNpB050X4XJ6IBy/b7pe7nrZdoMjy2ExFAY/R88wDBlP5aenvKExIsJ/WmQloA3p9LLVTzvkRg01epnNU+W33EWZiY4PKk9IFChFNI6KniGgZES0hom/I/eOIaB4RrZC/Y7VrLieilUS0nIjO0PYfTUSL5bEbSYpwImoionvk/peJaLp2zVyZxwoimlvqfRSDAJCUE5lMkxt1k1eQnTlMZfeL57eVLwyvrO7Ct/7wJr730JLCNES4vBT5JpZoNJTw54YV0k1SC9Dxmrz8NZRS7yVM8cr8SKanbLYZ7crkZbPNB91ekHD1EyhBqzorVPvwFToBoxRdYy79nQVr6e4UAUNxlFnOdh+xIhaHbDTK0VDSAP5FCPF+AMcDuJSIZgO4DMB8IcQsAPPl/5DHLgBwCIAzAdxERKqV3wzgYgCz5N+Zcv9FAHYKIQ4EcD2Aa2Va4wBcAeA4AMcCuEIXXJVCCCApR5NGH0rAKM4TBRaiwhREeanO3lq+cJWwd8jRTLbtKZyTkFsdN1xa+SavaJYBCTB5BTxnEy2pwqru6egMHUMUa2SFucym3Xg1qHDay5DF/xAv2+SlbQdo34XzUPRrfYSF20n7CRTn11ZH0h4/pjmNK/+4FAf9+2PWPMIEptz+whrrOUGCMcz3UMK05edXbsdTy7cGnldNShYoQohNQoiFcrsHwDIAUwCcC+AOedodAD4lt88FcLcQYlAIsRrASgDHEtFkAKOFEC8K5ynemXeNSus+AKdK7eUMAPOEEF1CiJ0A5iEnhCpGVgjXdBIkMIyjuCKjkwqjvAJMXoEpqnScX1OcfG7tqXBp5Zu8PP6NEjvhIFNZtsjnCADNCX8NJSjKq/QvWIa/zs/34IcnyivIh1Lkcj6mPAI1FJ9jYZzy5azllQ6hofzqudXWbw0B4SLF/NbtCxKMublLwen78YVbXsaXbns18LxqEokPRZqijgTwMoBJQohNgCN0AKilQ6cAWKddtl7umyK38/d7rhFCpAHsBjDeJy1T2S4mogVEtGDbtm2l3aBE11DMPpTcPmPYcIDTvuB8m1O+RDt4/nmmkXGYxqRTCZNXMctrlDcPxf99eAWOrTDhtaliygYUYQLTThy2mbwCRsVBzzxosJT2OR7WF6XegZ8PxXXKW07RNeZgv5B/OwpTtcwmr+D7cK71uc8SB2O1pmyBQkTtAP4A4J+EEN1+pxr2CZ/9pV7j3SnEL4QQc4QQczo7O32KF4yAyAkUow9FN6EEaCihTF7mhmm/Nlwt9Dur2JF4foinR6CU2CqCrgsKjtBpTjrvy5RkUPhxmI7QragWoVCMUM0/Ney70Itu01BUZJHteZU7D8VPmws7sTGK1YbTAaH7OjYtpZjQY7PJSzrsLfehnp+flrRXRnkRURKOMPmNEOJ+uXuLNGNB/ioj33oA07TLpwLYKPdPNez3XENECQAdALp80qooWaHNODZU7CCTlt752lYstqUH5DqcYlchLjzPOZEMvWAYp7xe2QsnNuqdcLjy5BP09T9P5xbwHHNhpP6doHllA/O5nrIGCT/fo/7n6v/7vQ/9mdt8KK7maZvYGHAfxYUNe4+F/fhWGIES1MmHMXkp+ofMYdpBbSBoqf4gTUs1mb4h8zI4tnQbgXKivAjArwAsE0Jcpx16CMBcuT0XwIPa/gtk5NYMOM73V6RZrIeIjpdpXph3jUrrMwCelH6WJwCcTkRjpTP+dLmvomSFQDJmD/EMshWHneDlnp/vQzGkUwoqWdOoOszienpDqYTJyyMwTI0yYLTsScsnbNe7Bpb/+7JJhiCBEuYZ2Dqe8D6U3LYtyks9h2HLCD9I+AcJX9+JjVqW/g73EBqKFBi2M/T6GHRPNm0uSEPR10Mzye90gA9FDZL6h4LXLGs0ylm+/kMAvghgMREtkvu+A+AaAPcS0UUA1gI4HwCEEEuI6F4AS+FEiF0qhFBDhK8CuB1AC4DH5B/gCKy7iGglHM3kAplWFxFdCUB5pH4ghOgq414CEUIU5UMxmryK9KHYnPLWihqYovdMk5XGdcqHdBj6mrxKFSgBI9qwJhT9XFNZ9F2mNa7CjKyD7jHMI3A7yfyOWNcKfN5HmBn9qr7almYJ1LQC3qvf0ithvzefG9n7mIICHuiANjk0yGQ4HCDIwwgUPw3F+mVXub97wL4ygp+fs54pWaAIIZ6D/X5PtVxzNYCrDfsXADjUsH8AUiAZjt0K4Naw5S0X9YITcbXIXpkaSiinvLkMpc4lyD/PT0PxK57e4H1NXiVObAwyKQRpFt5zS08njA8laM2lYpZe8Rcafnno5TGfGLSSb9BSH14hbjietR/XTUt+I++c78F6Su55W5LRBUpQE7MvlOnfBvTrTO93OGAeitr97rY91rIVOwm6XuCZ8iFRrzfpM6s5SGBkQqr+pvT0MtgdxOEqoTrL9K0F1UBee2+ntTPUy16JpVeCOnJ9j18HJYI6wYCOOIw93u2ILWOrcJFC2YL8gOIHCICP012ZvCw+lqLmoRSpoeifwvX1oQiVll2iBDnl+z0CJUBDsT4L59dW//V79TN9BwnvvkH7UjvFBXPUj/BhgRIS9YLVTHnz4pACqYRd4OgVrBSTlxvlVbZT3vk1ayi5bdukqbRHoNgXhwzspAJMAvnlMaXr55T3zt4uT0MJ46A1UYzJK78D8n7a1359GBOgek22UbnrVwuTh8nM4/OsVm7NjcT9o7zUcib2TtL18VkGT7o2FGQOtXX4Km3b5UH1M/ehPf/rfZ3yWf8y2MpTa1ighCQnUGKe/3UyWYEmPw1Gq79hVsm1Lb0SFMIaNl3zxMZcKr2WEZQ+ki7sBLXtgALZRtNB0XJhRuRAngYS4PMyLenhN7ciPw1bBxdmcqTKx0/b88OjsQWavPzvwxT5B+T7tQzP0mMG8h7Tk/Sfh1JYnnxs/iaFHoobbPLy11BsZU1n/N9pzocSoKFYoszClMFTHhYojUfOhyKXXjH5UIRAMuEjcIoYvQMlOOXl8dmTR2PmhDZrum7eARqK0rby8fhQ0vllzG0HjZysnV+AZuExq4UMMTWdpkdEmcwfxTjlbcfDLN+itIbCVQfM27Yy5G+bzrEtb++GklvyyATch/74wgqNwmO5g7ZOMqjdDKWz7iTOoEgpW5RXkA8lKLhm2GLCVOjC3RaVlwkog47tndYCFighUXUzpRaHtJi83MUjA9Y78vvmg+l8vQxBHXXQNxRU52k8TWuEtm9i6/kXVubwQtMWFRc0Ug3rlNezNwkmvTEH5RMU8RP0wSe/NFwNxTKACCLMCgyuXd/S+ajdtroT1InqwiD/aFKrR2FmytvyAPQABjNDmSya5UKgwSavgDys78tfC/JblkU/DtjnwhQzudK2OkItYIESEvViEzG1OKQ5zNTPaV9s2HDRTnm5O0bkO6JVjd9k3tCLpSLa8kl7BKPfPBSfQiDPtGXZDvKh+HZQAekMZvxHxEGmjaF0Frc9v8bJK8SIOsiEkv8s9SrmG27reXbmc9xRcYlO+aAABT9H9XlHOqsiERUhUAKCC2wMDmfd78sEO+UtwjUTXiAYn4XrE7M969z2up195nMChJInPzZ5NRbruvrw52VbAMDfpKU55f3i04FwE5cKVxsWBemYCNJQVAU0z0MJHv37RXnplwSNsG0j66BVdvVdfppeUGfu1VAKO5cgk5nuVA0Tym17HO7aT3n3EhTR5J6nd/ZWM6Lza5t7ERzlpedhSt/8Lv/yzjbc+ORKAMDY1lTAsipaGpb3GqQ9eDSUgCZmqzsqD9vjD/KtBc1D0Z/VI4s3mc8R3rT8sJnNakE5Exv3Gj5+w7PokaaNZMw+DyUrhGsmCg4rDq4E+ULHrWSWOqYcwwT/zlyVPSjKy9bg9JFd/szrYpZ892g6lg7JJHjDOuWF54Nn/gLFFrXnXm/IJ0yIqp6s3eTllCO/s/eG4toJM78pN7ExaHa4OY+gTtSzeoKWx0Nv5FZEihH51omspT7o5IIgzMcyWYGWVNxzrg2bwC7GZGVKIuh7KBkhMLo5ge6BNMa3pXzzsFXvMF/prAWsoYSgR5sZq8xAQRpKcFRRcL6FS68I4373uNxNASYvPw1Fr6h2e7tWmfOd8vp5xWgo+jpnxZi8Qjjl4zEyNkzV8TUnYyEW+yy8Xo/SCQqUcLbN5cxpKFnjfudi87X5eQeZiuxzIwrLa8vDOEnU41fLPZfWVO6zAfGY//sKmt+h5x30Pp1z/OufbR6KekZB7ysVj/kKV7+JjWOlINlj+Gqqk7e/UNPLbruPWsACpUji7lpehcfSWU1DsUSB5c4NoaEU6ZTPCRT/dFXHZfahaBU1wGnpbNt9KEEN2vY8vLOuDRqKpSy29BMxMna0g+kM4jFCKh6zfDAtQEMZCtZQ9L12H4o5Ksgm0PMJNQ9FnrNx14D5uDsqN+cR6DewaCitqZwRJE7ka1rLCv8BWX4++ah8lQ8laEBjtDJkRWDIrrIuJOPmgUqgDyXrfFepKRGzOuWDvg3jt1pFLWGBUiRuFJcobPw9A2m3QRid9h5bd3Bepc5DiZG5ortlVRpKoMkrWEMpXBzSv+PxlEPXSizaW9CSKUEdFACrwBhKZ5GKx5CIx4zvK50R2ve/C9P3znnw1wycbXM51b3nmw/9Zp97r3d+Y+SjocjdW7oHzMuFBCwnEhShqOc76BEoOQ2FyCzY3XSz2jwuW/CAxdcGAINSM8pFeVmzkvmVNvdI3X8qETPej3qWVn+WEIjHCK2puMds6j3H+Q3SkgD/ZfCrDQuUIlFRXvmV5VfPrQYAvPDuDgDmZcSL1VAKo7wCVGm5Px6jUIs7mpYL8TjlbSYBXd32MXkF3aIujGzagHliY3AZnfOc31QiZkxnKJNFKhFzNBiLzyvhM5E1KOzYuS63bXvnymRRqKHk/t/VZ19IUOWdjMcCl6dPZ4WxA/JGtNnzAMzzHnQBMGQRKDbTo56H34BMld+0refbHDrKy19w2i73PG8f4evnB4oRoSUZt05uVO/LqqlpZf/67143F7QGsEApkqRlHsrGXf2e/03Oz6DZxoC3s8xvtzmnvH3kA+S+H25DaR4moeOJoAqhbvuZvIoxOdi0FbMPRd8OFpxKoOSPaIfSWTRJgWJb7DPlI1AG07oPxVaG4Ml6roaSl4j+/48ee9ucAfLuM4QPxLTKbdCHqWwmLTd97Rpd4OirMcTI31EeFCUJeNtQ/vMazDd5BTnlDS8tzAe60h6B4j0mhAjhQ3EESrOvhiL8y6CVfUNe31NLWKAUSTxmNnnlayRB6rRtZO1nq3YddQF28qCRYC4s0jwqz50XPEr0W2042ClqHuEHLXfuEVohnPK5UW+hQEklYojHzRpKJiu0IIzC9MOYvDwdrcU0oZ5DfvnCrtGk8kjGY3bBppukhgM0FNP1Wl0wChSLwNH3x2L+PpS0R0Mxn6f79fK143wfiikvT3SUIQ9dawnyJ6UShRphuMACuCavAZsPJavKa1vROFzdqDYsUIrEXcsr74UOypFGk2wQJpPXUIgOyK9Cqv+CInliMX8NxY2TN47+HYchYI8e0c0b/qsN+xbDM5L1fvHPLFz0Miq6fExB6rSkZbmcQdfkZXfKKxNn0Ez7MCGm1ol08px8B23YeSjqvSfjZB1s6GU1mayCfAcDw/7XZ7LC7chN4djHzxwnw4Z9BgBZ4bYfu7k1l/Zgxvu8XIGSsguUYY9WXHgfusYQNFvfFOWlh+T7rVoQI/iavILC71XZx7Qm0d5UP7M/WKCEQO+fRzU7Ly+/oxuQ5o9RzUkA5spq60B1bAJFfeALsPsm1Olx8l9YUDUqs4aijegtHWBUqw3bnfJaGqZJdDLdjpYk1u7otaavO+WdcnsTU075eIysQRTKxGm6E6WhdLQk7Z2PZyFN/3fe1TtkvdYP3aZvG2zo78mkoXid8oXX6x2f6fpMVrj+Eo9Akene+eXjnCgvHxmpayilhMuq9tWUVGkUXu9pg4bn2y8nq45uTni+raLjRnklCjVbFVjRlDAHejjlEojFCC2phNXkFRRhqO69LZWoq4mNLFBCoNuBVRik3hlffv9iPLp4MwBgytgWAObRaBgn7rBl1Btm2Ra1X4U223C/v2HpRJVT09YB6r6DQpNXjqBlMjwTJD3mLy0NH5PXxFFN6B5IB86dSFpGvYE+FKGZvAz3MugOIhJlmbzUe+gfznhm34ddUsMzYrbVK12gpAs7sSCnfP9wGqPlYGrQoqGkEjEQmQdOyTg5o/YA7bop4V/3PHUmXThAAPx9KINaB276HLL6LO+4tpRVoKiytSTjhWbKTO4+shZzlfKhtCTtYcN6sqb7UHWmrSmOoUy2br6JwgIlBHqlOXbGWADel/y7V9a62zd94SgAZnNRGIHS1TdkPEdPL9ApH/BWbculA04ja5YdsM1Eo0arY1qThWmUaPKyaTZ+AkUt5dFraZS5hm+2y6/r6kNWOJNVzT4U/XMFhvLL99maigcKeVP+7n7tGe7Yk3v/YQWKek/Nybi1bihtTC+3zhpN0zMFa/QNZdzJeDYfiprTM+gxeWURjxGISPr2gs1IKj3jeXo0WcYsUJSmZBokBGkoSqCPaU1hwLYSsCxbczJeGEiRVe/Cf02/uIzysjrlPdaJwuOq7GqAWy+hwyxQiuCcwyejs70ZgNn/AACTRzvHTZ3x1Y8uc7dtncW6rtxicXrnoI+m/OLbAfg2XCBXGW0dQyoR843I6ZWze8e0JP0/ARzklLcsfRKkjbkmr1bHvLi73+xH0TtaJw9vWVdt78XiDbsRt/lQslkkfXwoblRRKhFKoNgEtG7q0c1eYSc25u7THMYKOD699mZz55POZPH25h73f9Ot9A9lMKbFed62KK9EjJBKxAp8KCqQJUZmwe3kKQVKgA9Ff4b55VD3NabFEXwmDUM315neh+rgx7U5gxXTOapsLcl4gZaU1jQUwNzOswKIxZx6Y/OhBE0kVeVS/pN6WcKeBUoRzN53NJQ1ydapx2LOSCy/QWze7Z2hbLv+b2971XiOXnHtoaGyDAFhw0o4mToG1QEk4jGrU159eKujNeU7dyJogG2LegtaekXtUh3cbotjXjVKZQKxrR5sm4eSyQokEz5RXsOqA0uG6hj8nPITRzUBADbtzoWAlqKh2HwUQ+mM2/nkC5S+vI7Xtm7ZaB+Bks46foGmRJ6GkhHu+nexGFkHYq6GEuRD8XmeypSnymka/Q95zKyG+9S0b8AslFS+ral4QUeu6lhz0n4f2WxuHordrObf3lWQRIe8V5NfqxawQAnBR97XCQD4ykkHWMOGdRIxKqjsO3oHPf8HLZcQj5E3IsW1MZsn6QG5ypuI+a/lpSqfyRaezmSRiMWQjJF1hNw3lEZTIobmRKzgM6Zeh3qAhuKxtZs1MJNmoMo1LmA9pGG3cReOFlWn9/nj9pPP2qytuRNZDeV46I0NAJxG3WspgzeAwT7qfv/k0QDg0RR0Aaj8Fyb0CX1W31xG5EazeQIhP3TV5pQf1ZxAPEZGH4xaTiQVjxWYdtXk0Bj5hLxnvALF6kPx+SiayrfDT6AYAgZ01MBgXKvScuwBDO3NhQ7xnMnLrqFkhDNoa0nF0D+c8dV+Ae8ipwoliJSWzhpKA/HD8w7DM/92sjNylx2M30dtUvFYwQtWFfWGCz4AAPjug0sKrlPmrnMOn4xxbSmj87tDjoaNzj5ZeWMB81BUQ/OzhScsy5UAQO9QGm1NCYxqTmDNDu/3HPQRZJiv6+n5Bm27eciORI1Eey3f5k7nmbz0uRSqQR7Y2W7XUEQuyst0L+9uc/wObU0Jq1CzaZmefLICY1uTaEnGPYJJCdkTD5zgGwo+pJlg/KK8ciYvb0er6uZhUzoAmIV4/1AGLclEgQaSK6sz6m5KxvO0gKwbhm5qFwpVbhU2bFuNW9em8u9jKOMVKKbRvy5kTPVbHVf+IlMaQ+ksSIb95gu1nMnLZ00/+axapanU9Ex0jcP0TlU5cxqK/XPC1YQFSgimjGnB/uOdT+qmEjGkEjGjI1gJi0S80OSlOopRPiPND//4KQDAhPYmp+EOF47GOkc1IZMV6O4v7MD0RulnLhlwBYoh2ifrdKImLSt3Lxm0puI4oLMdewa85dBHfUWZvCzbpjRUuZRZIkg7UOYHvfGrkWdzMi7Ne2bhqmzhfpMMRzXbBcpQJpubW+Gz9Eo8FkNryjsvQdWhtqa4r0lDlb0paY/yGkpnMcpi8lKdk3JmmzWUNFpTcYxqTqDHMNNeTQJ1NJRMwX6nfHGr8zithds615nvtX8o7d5HvvbgaiiyXpgiqPT6avShyGvGSg3FpI0NZrJIxp1+ID+N3H3YF6gcSmfRnIy5A52BocJyDGj5+i1MqgQKaygNzKimBPYMOo1KH8198oh9ATiRQfkVTa3FNLmjJTD902ZPQnMy7qlUSgjMlqaRxRt2F1ynh176fS9bdVq2tYwcDcUcSgsA23oGMaG9CW1NCQxlsh4hEvRlPx1PqHDeWlBqQP7Eks2YftkjWK992U41WtXobQJFpW8KI1XPszkZQ3tTvGCA4Dhkc2YivyXCs1mBgeEsVm7tKTi2ZzDtmuZMk10HhjPYvmcQiRihJU+gKM1rbGsKA2mzVgrkNIzWlE+Ul6ah5Gum6vo2ea/57y2TFdjZN4zWVBxjWlLGIAgVudScjHk6+uFMzmzoDJLMI2klMG1mOUXvYAbj2s3ag7pGmQdNJi992RnjPBSlobg+FIPJK+0sYpmUWrznOy55PhRj6HI6i1Qi7tbLvuHC+qsPIPy+xeMXKFELWKCUQFMihudXOotA6o1LLQff3pwoWC9ph4zemdzRjEtOPqDAnnzcD//sST9fQznrhmcBAB86cAIAYNX2PZ70t+8ZxPy3t0qznF27+PYf3sSidbsAmCvhsPShJGIxY5w+AGzuHsA+o5vdBqGHV+rXBE3MU859wDt6H8pkMXGUEy13+wtrAABvaQJUde6qMenp6KTzfCj6M1HCuiUZR1sqUaBpqZFpmytQ7A32zfVO2T523V88+4czWby+dpcrUEw2+xvnr3COZQX27WjxRPlt7R5EjJxBiBBmoSaEcNNIxOz+teFM1tWO80fdqmNWx/O125895Xxx8e3NPVi+pQdPLNlSoAn1D2XQkopjVHPSo8GosGHAqde2jk+1o0kyStIkDAbTGSzd1O2+1/xzegfTrimqKREzpqE0yWScjPW7byiDVDzmCl+jySuTQSoRc0PKvfNuvBqKyVIwOJxBUyKXh2lApGu8Ju10IN/kpT3XZ97Zhv97fUPBNdWABUoJbNw9gNXbe7FxVz8u/e1CAI5zV7HfuFYs2djtueb3C9YBAEY3JzGhvQlZkRstvfjuDmzpzjntm5Nxj4ayZntujsD0CY7prTtvlHjJbxZiycZuGZkUs3bm9y5Y726b1OQ+2TEkLXMzugeGsbarD/t0NKNZmkh004J+jck0ovPOFs0BnfVqKBNHN3nOVUIByAkK14di0VAG8oSCXjZV5uZkHO3NiYI0VCNWnWy+QPmvPy13t3/+xaMBwHWsK1RHrBYO3dbjDcwAgJuefheAY1Ia15byDES29gxgQnsTOlqcMpgWddyuzVtJWrTK197biYHhLNrUnIVhs4Yyvq3J87/i+ZXbAXifc0/e8+odSqMtlcDoloQn9LlvKOOa0poS9nkXPVKgT/QRKK+t2Qkgtxhifme/bHMPZk5oQyIeQ3tTwmgWVvmMNUQoAsC72/Zg6tiWnDnK0Jmv2d6HiaOb3Tkz3gm60kqgorwsGkpTIuaa7noGCgNbNuzqd9PPPw7k3lGHQUOZe+sr+Kd7FoVeCy5KWKCUQVfvkKupnPb+Se7+o/cbi/d29LkRUEs3drvRO7EYuSNW1fDeyjNf5WsoutP54H1GoSkRK6hkK7fmNJakdH4GzZ7NX+pD7RvXlkIqT0NSPPvOdgyls/jY+yflNBR99rG7JETcOj8EcEJ9H9BGUXonNpzJYnxbCjY/tBoFtqTiSMXN/qzNuwfwj791lvVW5guTD6Up6XQ+e4a8M+7ViK/doqH8VH4n/bNzpmFcWwrHTh/ndvyK9Tudji+dEdi3o9kTwZVPVjgz+vUybu0ZROeoJnRKbW37nkKBtLs/9w7bmpylPPI1oU/f/AIAp+7FqHAgoTrv8dKUlB+5N136D7/4wf3dffnn9A1l0NqUwIGd7Vizo88V2HsG0q5Q3n98KzZ3DxjDvJWwnCTDp/sNgRabu53Q+8vPOhhAoUB5b0cvDuhsBwBMHtNcsAI44AxyYuQMFEw+rY27+jF9QhuaE4V1W7GzbwhTxjS7EWneeiUHMSmztqfOaUrEXQ0l3//29mY5GKVcmfNRn8tobTJrnbbrKg0LlBI4oNNpYFt7cnNLlIMYAPaXWsS6LqdCf+6XL3muH5snUPQJj0ChhqKEx21fOgbNSces0G0YtShS0gkaZvXaFZqW8KZZEnYAABmwSURBVPAbG7F+Zz/GtCYxvq2pINR55dY9rkZ2+LQOV6D0D3uFAQBMGNWEXRaB0jMwjOv//I5n3y5thYDhTBbNybg7+gK8I0XV2afiMbQ1xY0ayk1Pr3S3la9FFwqqQ2xNJdDWlIAQOaHWMzCMh+W30FVnaPJ/OGWR/oumwoX+1MBh5sR2TJ/QhtXb7euOdbY3eUJus1mBp5dvw8Zd/eiUnaxJw3l2xXZ3W60jp3dQeod41H5j0ZQodIyrzlt93zz/PpS2+Mkj9sWNnzsSQKGZcWffEMa0JDF5TIv7/+L1u/Hiqh1okZ3rYVM7IASwcN3Ogvvozjd5GRzV/3zvGwCAj82eZDynuz/ttsOxrSlj/VuysRuTO1qkn7PwnW7pHsSk0U2uDyRfUxpMZ7BkY7cTWm/QUJQgddf8yxNaw5kseocy6GhJuufkm1tVv/Cvp78PQKGG0jeUdgdr7mK08p3qQmRX3zDWbO/FeTc9b6w7lYAFSgnccIHTqBat3eXuU+YCwDF5AcBaaQ9XL//fzjgIQC7G/S/vbMOqbV5fCFCooagKpxr8aEOkjT66TlpUZX3SnEKNmoczWXxNfqhnTEsKnaOaCirh392Rm3Q5ujnpRu8oHwLgmEVi5HQMNg3ljhfWuL6R754zG6l4zPUxAU4jbsobrevCTXU+o1uSaGsqNFcB8KzAOnVsrpMDgGWbunGzNDWNb0u556p0rvrjMlfIq1GkPurXn7UybfYNZvDm+t3uIOP6ee/gF39ZBQC45cI5mDGhDau27SnQGqeNc8r2tVMORCqR832pJXimjm3F5A6nk316+TbPtdf9aTm+//BSAMDPPn+U20Gp9754/W4c/N3H3fNPmz0JLal4wYhYdYLTZL3d2u197z0DzjpeRIQ2ab7SNZSB4Qx29Q1jn45m7CMFwvItPfja75zBx2Cevf9L2uRdPQ/AEQSpRKzAUa0PhjpHOcJX7+yFENjdP4zRUqiOzvPlqPt8evk2HDGtw3Go52lqSzd2Y/ueQUwc1exGiuULpVuedTSDx5dsdkPKdXNTX16UWL7JTNXdjpaEW+/yzYdq+R2lbeWbOvX3o7QkNUi44c8r3GP/+aflOPknT2Ph2l2477X1qAYsUErgkH1HozkZw41P5kbB+45pdrf3lw3zvR29bmM9/+ipuPSjBwKAG6Vy45Mrccp/PYMpY7yRX03JOFpTCWzfM4jhTBY9MqJMVcD25oRnzScg51c4+aBOd5R41JXz3OOrtu3Bidc6Ycn3XHw8nv3WRwHkOgb9i4CpRAwT2gsFyoETRwHIfRNG2Yfvlf4hwDEFZYUz4t7Sbf5+uW73v+jEGRjXlsJOKVBee68L67r6C+Z26H6krt4hxGOE0c0JpOIx3P/6Bk9H/eb6Xa5vAoA7wr/n1XXoHUzjrBuexYL3nFHy+PacQDGF/qplPPRRpN6RrZCmxgMmOlrp8yu3466X3sMN83MNu3NUE2bvOxrdA2l37goAvLFuF9Z19eMrJ83E2LaUJzpQPbtLP3qA29H/6rnVntGwXv8+ftg+bnST6pwffWuTe/wrJ80E4Jid9EHMb19ei+9JoTRrktOB6RF1G3b14/YX1rjarlo7Sn9Wqp5MHNWEo/Z31rp7Z3OPW25lKlMDqXzWbO/FD/7olKG9OeHMIB/yCosDvvOo+//o5qSMJlPRilnMuPxR9A9nsI8UvqNbEgUDqr+/cwEAJ1JyXFsKq7b3uvVmXVcfPn6jE/gyaXSzW9YdmpnxrQ273UHC2NakKyD1tqO0i0kdZl+QGmR1tCZdH6CunQPAwrVO3ZwhLR35X+t8a6MzgHvgkhNcAarazy3SFAYAj7yZe/+mwWQlaGiBQkRnEtFyIlpJRJdVMV9c/OGZ7v//cc5sdzYwkDN/XfXIMldLOXHWBPe4GsUppoxtwbHTx7n/NyViOPX9E7Gzbxgvr+pyG4YyaazYsgcvrtqB6Zc9gi3dAxgYzqBnMI2vfGQm/veLR2OGtHnr/GFhboRy5H5j3TJeP8/p+HRt4ryjpqBzVBN6hzKe0f+O3kFMGdOChd89DUDO9PBhGXn2+Fub3XMPm9qB93b04V2t81rX1Yc5V83D7S+swejmBN6+8kwAjm/g3gXrkc0KfPrmFwE4jffaTx+G8W2pggmUz6/cjkxWgIgwe1/HEX7N47kvGn7yf553t9dcc7bb8F94d0fB51KVyQvIdZITRuU6PvXdij8t3YKfPbUS6UwWO7UG/i+nOWaJ73/yUKQSMTzy5mZ89//eQj5z9nfe78euewZPvr0FN85fgXN/5pRz1iRHUOtmmJdXdQHIOakV33vImRCbP7omIrQ3OfepRrT6SPazx0wD4EzkXKUJte88sNjd7mxvQkdL0vX9AMCHrnkSQG7krYRvn2byWiWF/aTRzehoSWJsaxJrdvS55rhzDnfC6adPaMO0cS04eJ9RnrJffn+uDOrDU7pf7Lbn17jbygw0Xhvw6KZEZQIa1ZzEtp5BTL/sEWzY1Y9Nu/vxnAwuOGLaGJx8UCfe29GHrYY09h3TjEQ8hjGtSc/A7ZyfPue2kzu/fFzOFLnHEf7rd/bh6keXYXJHs6sV64J3xZYenH69EwnY0ZLEqKYERjUlsHGXd+C1cusejGlNYvr4NrQ3JfBe3icaFq3dhaZEDIdO6cCE9hTGtiaxUFpLlNaUz50vvleVFYkbVqAQURzAzwCcBWA2gM8R0exq5f/Ppx/kbudHJJG2ltYnfvocgJw5AciN8BWvrO7CzM42fFgKnaZEDCccMB4A8De/ehn/IWfVK+eyCh0GgON+ON/VEGZNHIWmRByHTe1wK/RTy7dixZYet2F85H2dSCVibsewuXsADy7agK1yRPybvzsOk0Y3uxrX9x928h4YzmDx+t04bfYkt4NWk7v++OYmPLV8K15812mw/3r6+3CaFDbzlm7BwHAGA8MZvLVht6udfPusg91IGvW43lifMyGeceg++Owx++G1756G42eOx+L1uyGEwPXz3sEbmontm7JD/99nnJGjyeGrv4/5b291t0+SS+ooU1FX7xBWb+/Fz57KaTfdA8PoH85g8Ybd+M8nluPrd7/ujgbPO3IK/vEUR+tMJWI4ar8xeHn1Dk/eh0iBp9bqAoAv374A183L+ZDmyFF9S9IxRz3w+npX88nvfH/z8lp85uYX3A4EAC45+QAAuU8nqNWv1cj3vCOnuBNzO1qS2NoziO17Bj0dzNmHTQYRYerYFizZuBtvrMulDwA//szhAJyRP+B0wOlMFks3dmPura8AyPk/xrc3eVbg1vnE4fvi7c092Lx7AAvWdGFr90DBKH5mZxsWvrfTLd+Di3LBG5ec7DzvKWNasH5XP7JZgfsX5o4fM8MR3KM0k+dL7+7AR3/ytPv/h2d1utGSC2Tk2NWPOCbOlmQcx81w2l5ne5Pr2F+2KRe1+e0zD8ZhUzsKfFvXPu5E/iXi5A4adRPUtY8vdzW9/ca1gogwvj3lCY558d0deOHdHTjr0H0Qi5FjKtWE3bJN3bjludVokmHLRIST3teJF1ftwC3PrsJwRuD8o6c685qScbzynVPx9x+eAQB4Pe+dVoL6+dRX8RwLYKUQYhUAENHdAM4FsLRaBfi3Mw7Cfz6xHDMntBccu+kLR+GS3yx0I2ryO4bf/8MHcf7PX3T/P2xqB/7jE7OxpXsQRIQxrSlMGdPi+V600oJu/pujcN5NL7iTG5XAmaQJtl9fdBxO/snTrr36A9PG4PCpHbj9S8cAcDrZn5x/BP7192/gG3cvcq9T/p8zDtkHB+8zCvcuWO8JNVbrmrllihGWb+lx8xndnMA/njLLtXlf89jbuCbve+iXnXUw/nrONPf/uy/+ID598wu4UHZM/372+/F3mgZ46sETMW/pFhz7w/lu473y3EMAADMn5LSx6Zc94slH+axMXH7Wwfh7mcdMGWTxtwbb/gkHTMC4tlyjf3TxZvfbN587bj+PsBpMZz1mlk8csS9+8EmnnKO1AAOdBy45we3cjtp/DADgm/c4zud9Rje7JqYl3z8Dh1zxBABgwXs78df/69Sd+y85AUft5wgkZWp9cNFGPLjICSo4aNIoXPfZD7j5Ke1lzlV/xvdl2QDgYmkSm9DehGfe2YZzf/a8G+57wTHT3Pel6sfVjy4rCCZ5nzSZnXjgBDfq8EfnHeY55/RD9sFNT7+L4380v+BZfFUKxo8fNhn/74G3MOPyRz3Hv37qLHcJmvHtKTy3cjtmaqaw1797mhvwcpDW3v7l92+427/9++M8z+rS3y7Epb/N5fHm9053fZCbdw9gxdY9+PZ9ublbd3z5WJwkB35KoHz7D4vx7T/ktKwbLjgS+0qz86J1u3Dj/BWYv2yLOxBqS8Vd8/FwRuChNzZi1fY9+Phhk/FjKZROPNBpZzM72/Dgoo24d8E6TGhP4aI7HLPdmYfu4+Z37IxxeHDRRlwlheKJsybgqr86FMlYDLEY4V9OPwifO3Y/zOws7KeipmE1FABTAKzT/l8v91WNSz96IBZ+9zTX7KLz8cMm45jpTkM/9eCJbsegOGb6ONwmO/ezD5uM84+ehtZUwrWbAsCfvnmSu/2/cq4D4GgGD3/tRPzhqyd40jx86hh3e//xrZ5ji9btwjHTx3k6wE8fNQXXaA2+JRl3G0JzMo5/P7tQ4dO1IwDuCF3xuWOd+TjxGLkjo3y+ctJMt9ECwJHTxuDQKaPdzvjkg7xC65SDJwLIjQTPPnwyvvjB6QAcwXj/Jd7nAACv/L9TXZ8VAFz310e4I3kA+Ksjp7id08RRzZ4oPQC46lOHYs01Z2PauFa8/J1Tcf7RUwvyOCCvgepmUAD4wScPcTu4eIzwjVNneY7/9ZypOFIKAwA4+X0TPceVPwBwQoJ/+3fHuR264gPaO4/FyI1OUhw+tcPz/+ePy4X+XiHNZw9ccgKOmOakc+r7c2VQZi79PonI+Cze+v4Zbt264hOzccnJB+ArH5np1ge3vNPG4DIZ9qvz538+Cd8+09l//tHTCo7/xzmz8U/a8ztGMxEDTqeqnjUAnJBXTwHgW2cehBMOcPbvbzAL//CvDvPUy2/Jct6zYB2Wb+lBSzKOj7yv073PpkQcF53orePXf/YIV8B/7xNO+7lO06rPO2oKXrjsVPd89T7f2tCNHz++HC3JOP7zM4fjLCkwPvUBp0v71n1v4su3L4AQwDc/9j5cc97hbhqnvX+Sq+WffdhknHXoZDQl4m79bk7GqyJMAIDq5UtfxUJE5wM4Qwjxd/L/LwI4VgjxtbzzLgZwMQDst99+R7/33ntVK2P/UAb3LVyPM2ZPKrCFK5w1kkpXFLd2D2Dltj344MzxHmEBOPHsz63YjqUbu9E/nME15x3uRq/obN8ziN8vWI9PHDEZU8d6O6x0JussFJnJuh9JykcIgcUbduOQfTsKju8ZTOOHjy5D72Aa535gX5x4YKcbmZKfxm9eXotZE9tx3MzxBce7eocwb+lmLN6wG18/ZVbB83x32x5s6R7Asyu240snTLc+7xdWbsc+Hc0FDWw4k0WMCCu29qC9KVHwHBTff3gJnl+5Hb+8cI6xU+ofyuDqR5fiyGlj8WlDx5vOZPHMO0448KeOnOL6xRS9g2ms3t6LN9bvwocP7MR+483lWLl1D0Y3JwrucziTxWA6i4ff2IhRzQmcdehk4zu768U1+O6DS/CtMw9yzUiKrd0DmLdsC37+zLs48cAJ+PezZ7t+JsWfl27B5u4BHDqlAx0tSc9AKAxvb+7G6OYkfvnsKnziiH3dTljRP5TBpt39eHDRRnzhuP0K7jOTFVixtQeJGCErgPdN8loAAOddJxMxrNiyB9PGteDDszoL0tjWM4iWZBxbegZwQGd7wbNKZ7JYtb0Xj7+1GYdN7cBHD/IKfcAJTnliyWaMa0t56m7PwDB++ZdV6B3K4EMHjsfMCe2YNq7Vk8fAcAbXPPY2jpk+DvGYE/hy4ERv3Xz8rU1YurEbw1mBWRPb8VdHTilo69mskIuZVkZHIKLXhBBzAs9rYIHyQQDfE0KcIf+/HACEED+yXTNnzhyxYMGCKpWQYRhmZBBWoDSyyetVALOIaAYRpQBcAOChGpeJYRhmr6VhnfJCiDQR/SOAJwDEAdwqhCj8yAjDMAxTFRpWoACAEOJRAI8GnsgwDMNUnEY2eTEMwzB1BAsUhmEYJhJYoDAMwzCRwAKFYRiGiQQWKAzDMEwkNOzExlIgon4AfqHFHQB2+xzfD4B51bvwaZR7PMw5QeWMIo8o7qPcclbjPqrxzqNIox7eeZhzGuGdA41Rzmr2RwcJIQqXI8hHCLHX/AHYFnD8F+VcHzKNso5HUc6I8ojiPsp9H9W4j4q/82qUsxplqEY5q3gfdV/OavZHABYE5SWE2OtMXkHrNz9c5vVh0ij3eJhzyr3PMOdEcR/llrMa91GNdx5FGvXwzsOc0wjvHGiMctZLf+Syt5m8FogQ69FU6vpqweWMjkYoI8DljJpGKGc1yxg2r71NQ/lFja+vFlzO6GiEMgJczqhphHJWs4yh8tqrNBSGYRimcuxtGgrDMAxTIfZ6gUJEtxLRViJ6S9t3BBG9SESLiehhIhot9yeJ6A65f5n6Bos89jQRLSeiRfKv8Es81Stniohuk/vfIKKTtWuOlvtXEtGNlP+lnvooY8WeJRFNI6Kn5PtbQkTfkPvHEdE8Ilohf8dq11wun9dyIjpD21/JZxllOevmeRLReHn+HiL6n7y06uZ5BpSzIs+zhDKeRkSvyWf2GhGdoqVVsWfpS5hQsJH8B+AkAEcBeEvb9yqAj8jtLwO4Um5/HsDdcrsVwBoA0+X/TwOYUyflvBTAbXJ7IoDXAMTk/68A+CAAAvAYgLPqsIwVe5YAJgM4Sm6PAvAOgNkAfgzgMrn/MgDXyu3ZAN4A0ARgBoB3AcSr8CyjLGc9Pc82ACcC+AcA/5OXVj09T79yVuR5llDGIwHsK7cPBbChGs/S72+v11CEEH8B0JW3+yAAf5Hb8wB8Wp0OoI2IEgBaAAwB6K7Dcs4GMF9etxVOeOEcIpoMYLQQ4kXh1Lo7AXyqnsoYVVl8yrhJCLFQbvcAWAZgCoBzAdwhT7sDuedyLpxBxKAQYjWAlQCOrcKzjKScUZUnqnIKIXqFEM8BGNDTqbfnaStnJSmhjK8LITbK/UsANBNRU6WfpR97vUCx8BaAT8rt8wFMk9v3AegFsAnODNWfCCH0DvQ2qQJ/t0oqpq2cbwA4l4gSRDQDwNHy2BQA67Xr18t99VRGRcWfJRFNhzPKexnAJCHEJsBp2HC0JsB5Puu0y9Qzq9qzLLOcinp5njbq7XkGUdHnWUIZPw3gdSHEIGrTzgGwQLHxZQCXEtFrcFTPIbn/WAAZAPvCMSv8CxHNlMe+IIQ4DMCH5d8Xa1jOW+FUogUA/hvACwDScNTffCod5ldsGYEqPEsiagfwBwD/JITw0zJtz6wqzzKCcgL19TytSRj21fJ5+lHR51lsGYnoEADXAviK2mU4rSrhvCxQDAgh3hZCnC6EOBrA7+DYowHHh/K4EGJYmmmehzTTCCE2yN8eAL9FdcwNxnIKIdJCiG8KIT4ghDgXwBgAK+B04FO1JKYC2Jifbo3LWPFnSURJOA32N0KI++XuLdJUoMwvW+X+9fBqTuqZVfxZRlTOenueNurteVqp5PMstoxENBXAAwAuFEKofqrq7VzBAsWAitogohiAfwfwc3loLYBTyKENwPEA3pZmmwnymiSAc+CYempSTiJqleUDEZ0GIC2EWCrV5R4iOl6q6RcCeLCeyljpZynv+1cAlgkhrtMOPQRgrtyei9xzeQjABdI2PQPALACvVPpZRlXOOnyeRurwedrSqdjzLLaMRDQGwCMALhdCPK9OrkU7d4nay99of3BGzZsADMOR7BcB+AacCIt3AFyD3ATQdgC/h+MAWwrg30QuIuQ1AG/KYzdARtjUqJzTASyH49T7M4D9tXTmwGkA7wL4H3VNvZSx0s8STuSOkOkvkn8fBzAeTpDACvk7Trvm/8nntRxatEyFn2Uk5azT57kGTvDGHllPZtfp8ywoZyWfZ7FlhDNA69XOXQRgYqWfpd8fz5RnGIZhIoFNXgzDMEwksEBhGIZhIoEFCsMwDBMJLFAYhmGYSGCBwjAMw0QCCxSGqROI6B+I6MIizp9O2srODFNrErUuAMMwzoQ5IcTPg89kmPqFBQrDRIRc0O9xOAv6HQlnMueFAN4P4Do4E2O3A/hbIcQmInoazhpmHwLwEBGNArBHCPETIvoAnFUFWuFMTvuyEGInER0NZx20PgDPVe/uGCYYNnkxTLQcBOAXQojD4Xza4FIAPwXwGeGsZ3YrgKu188cIIT4ihPivvHTuBPBtmc5iAFfI/bcB+LoQ4oOVvAmGKQXWUBgmWtaJ3LpKvwbwHTgfP5onVzmPw1meRnFPfgJE1AFH0Dwjd90B4PeG/XcBOCv6W2CY0mCBwjDRkr+WUQ+AJT4aRW8RaZMhfYapG9jkxTDRsh8RKeHxOQAvAehU+4goKb9fYUUIsRvATiL6sNz1RQDPCCF2AdhNRCfK/V+IvvgMUzqsoTBMtCwDMJeI/hfO6rA/BfAEgBulySoB54NiSwLSmQvg50TUCmAVgC/J/V8CcCsR9cl0GaZu4NWGGSYiZJTXH4UQh9a4KAxTE9jkxTAMw0QCaygMwzBMJLCGwjAMw0QCCxSGYRgmEligMAzDMJHAAoVhGIaJBBYoDMMwTCSwQGEYhmEi4f8Dpbdw2ooDd/AAAAAASUVORK5CYII=\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": 10, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztvXt8nGWZ//++ZibnU5M2SdOm0EJbaDkVqKUKKoi0oKuggtbdlbriVhFXd7+6u7Bfd3Fx+SmecD0vQrVUvwKLKKAgVk4CQksqh1JKaXpOT0mbQ5NMMpOZuX5/PM+kk2QyM0kmyWTmer9e85qZa577nnueps/nuQ73fYuqYhiGYRip4JnsARiGYRhTBxMNwzAMI2VMNAzDMIyUMdEwDMMwUsZEwzAMw0gZEw3DMAwjZUw0DMMwjJQx0TAMwzBSxkTDMAzDSBnfZA8g3cyYMUPnzp072cMwDMOYUmzevPmoqlYnOy7rRGPu3Lk0NDRM9jAMwzCmFCKyN5XjLDxlGIZhpIyJhmEYhpEyKYuGiHhF5CUR+a37vkpENojIDve5MubYm0SkUUS2i8jKGPv5IrLF/ey7IiKuvUBE7nXtG0Vkbkyb1e537BCR1en40YZhGMboGImn8XlgW8z7G4HHVXUB8Lj7HhFZDKwCzgAuB34oIl63zY+ANcAC93G5a78OaFPV+cDtwG1uX1XAzcAFwDLg5lhxMgzDMCaWlERDROqB9wJ3xpivBNa5r9cBV8XY71HVgKruBhqBZSJSB5Sr6vPqbOJx96A20b7uBy51vZCVwAZVbVXVNmADJ4TGMAzDmGBS9TS+A/wLEImx1arqIQD3uca1zwb2xxzX5Npmu68H2we0UdUQ0AFMT9DXAERkjYg0iEhDS0tLij/JMAzDGClJRUNE/gpoVtXNKfYpcWyawD7aNicMqneo6lJVXVpdnbTM2DAMwxglqXgaFwLvF5E9wD3Au0Tk58ARN+SE+9zsHt8EzIlpXw8cdO31cewD2oiID6gAWhP0ZRhpYfvhTjbuOjbZwzCMKUNS0VDVm1S1XlXn4iS4n1DVvwUeAqLVTKuBB93XDwGr3IqoeTgJ701uCKtTRJa7+YprB7WJ9nW1+x0KPAasEJFKNwG+wrUZRlr41h+286XfvDbZwzCMKcNYZoR/DbhPRK4D9gHXAKjqVhG5D3gdCAE3qGrYbXM98DOgCHjUfQDcBawXkUYcD2OV21eriHwFeNE97hZVbR3DmA1jAMe6g/iD4eQHGoYBjFA0VPUp4Cn39THg0mGOuxW4NY69ATgzjr0XV3TifLYWWDuScRpGqrR2BwmETDQMI1Wybu0pwxgJrd1BwpEhtRWGYQyDiYaRs4TCETp6+vB54hXpGYYRD1t7yshZ2nv6AAhFlFA4kuRowzDARMPIYdq6g/2ve0MmGoaRCiYaRs5yLFY0+iwZbhipYKJh5CxtJhqGMWJMNIycpdUfKxoWnjKMVDDRMHIW8zQMY+SYaBg5S2xOwyb4GUZqmGgYOctAT8PCU4aRCiYaRs7S6u/D607ss/CUYaSGiYaRs7R1B5lZXgiYp2EYqWKiYeQsrd1BZk2LioZ5GoaRCiYaRs7S2h2krqIIgF5LhBtGSphoGDlJTzBMT1+YumkWnjKMkWCiYeQkbe7EvllRT8PCU4aREklFQ0QKRWSTiLwiIltF5D9d+5dF5ICIvOw+3hPT5iYRaRSR7SKyMsZ+vohscT/7rrvtK+7WsPe69o0iMjemzWoR2eE+VmMYaSAqGrXlBYhAwETDMFIilf00AsC7VLVLRPKAZ0Ukuk3r7ar6zdiDRWQxznatZwCzgD+KyEJ3y9cfAWuAF4BHgMtxtny9DmhT1fkisgq4DfiIiFQBNwNLAQU2i8hDqto2tp9t5DrRcFRhnpcCn8dWuTWMFEnqaahDl/s2z30k2ursSuAeVQ2o6m6gEVgmInVAuao+r6oK3A1cFdNmnfv6fuBS1wtZCWxQ1VZXKDbgCI1hjImgKxL5Pg+FeV4LTxlGiqSU0xARr4i8DDTjXMQ3uh99VkReFZG1IlLp2mYD+2OaN7m22e7rwfYBbVQ1BHQA0xP0NXh8a0SkQUQaWlpaUvlJRo4TdDddKvB5KPSZaBhGqqQkGqoaVtUlQD2O13AmTqjpVGAJcAj4lnt4vL0zNYF9tG1ix3eHqi5V1aXV1dUJf4thwAlPI8/roTDPY9VThpEiI6qeUtV24CngclU94opJBPgJsMw9rAmYE9OsHjjo2uvj2Ae0EREfUAG0JujLMMZEX9jCU4YxGlKpnqoWkWnu6yLg3cAbbo4iygeA19zXDwGr3IqoecACYJOqHgI6RWS5m6+4Fngwpk20Mupq4Ak37/EYsEJEKt3w1wrXZhhjoj+n4fVQkOe1RLhhpEgq1VN1wDoR8eKIzH2q+lsRWS8iS3DCRXuATwGo6lYRuQ94HQgBN7iVUwDXAz8DinCqpqJVWHcB60WkEcfDWOX21SoiXwFedI+7RVVbx/B7DQMYlAj3eczTMIwUSSoaqvoqcG4c+8cStLkVuDWOvQE4M469F7hmmL7WAmuTjdMwRkJgUHiqPWYXP8MwhsdmhBs5SWx4yhLhhpE6JhpGTjIkEW4LFhpGSphoGDnJAE/D5mkYRsqYaBg5STAUwSPgs/CUYYwIEw0jJwmGI+T7nD9/m6dhGKljomHkJMFQhDyv8+dfkOclEIrgTA0yDCMRJhpGThIMRyjo9zSc54BN8DOMpJhoGDlJMBQhP+pp+LyAbcRkGKlgomHkJMFQbE7DebZkuGEkx0TDyEkGiIbraQRsroZhJMVEw8hJ+sInEuGFedHwlHkahpEMEw0jJxlYchsNT5mnYRjJMNEwcpJATCL8hKdhomEYyTDRMHKSuIlwK7k1jKSYaBg5STB0Yp6GldwaRuqYaBg5Sd+gZUTARMMwUiGV7V4LRWSTiLwiIltF5D9de5WIbBCRHe5zZUybm0SkUUS2i8jKGPv5IrLF/ey77ravuFvD3uvaN4rI3Jg2q93v2CEiqzGMNBCMqZ6KehxBC08ZRlJS8TQCwLtU9RxgCXC5iCwHbgQeV9UFwOPue0RkMc52rWcAlwM/dLeKBfgRsAZn3/AF7ucA1wFtqjofuB24ze2rCrgZuABYBtwcK06GMVpiZ4RHxaMvbGtPGUYykoqGOnS5b/PchwJXAutc+zrgKvf1lcA9qhpQ1d1AI7BMROqAclV9Xp2V4e4e1Cba1/3Apa4XshLYoKqtqtoGbOCE0BjGqIlNhOd5BTixMZNhGMOTUk5DRLwi8jLQjHMR3wjUquohAPe5xj18NrA/pnmTa5vtvh5sH9BGVUNABzA9QV85weGOXr76yDbCEbsDTjex8zTyfFFPw0TDMJKRkmioalhVlwD1OF7DmQkOl3hdJLCPts2JLxRZIyINItLQ0tKSYGhTiyfeaOZ//rSLvce6J3soWUesp5Fv4SnDSJkRVU+pajvwFE6I6IgbcsJ9bnYPawLmxDSrBw669vo49gFtRMQHVACtCfoaPK47VHWpqi6trq4eyU/KaPzBkPtsVT3pRFUdT2NITsM8DcNIRirVU9UiMs19XQS8G3gDeAiIVjOtBh50Xz8ErHIroubhJLw3uSGsThFZ7uYrrh3UJtrX1cATbt7jMWCFiFS6CfAVri0n6HHFojsQmuSRZBehiKJ6wsPwegQREw3DSAVfCsfUAevcCigPcJ+q/lZEngfuE5HrgH3ANQCqulVE7gNeB0LADaoavVW+HvgZUAQ86j4A7gLWi0gjjoexyu2rVUS+ArzoHneLqraO5QdPJbpd0TBPI71ES2uj4SlwvI2giYZhJCWpaKjqq8C5cezHgEuHaXMrcGscewMwJB+iqr24ohPns7XA2mTjzEZ63PBUd9A8jXQS9ShiRSPf6yFkOQ3DSIrNCM9g/BaeGhfiexpi4SnDSAETjQzG3xcVDQtPpZPoXuDRBDiAz+sx0TCMFDDRyGB6+nMa5mmkk2juomBQeCoYsvCUYSTDRCODiYalui0Rnlb6w1NeC08Zxkgx0chgevospzEexEuE53k9hCImGoaRDBONDOZEItw8jXQybMmthacMIykmGhmM5TTGh2CcRLiFpwwjNUw0Mhh/0HIa40FgmPCUiYZhJMdEI4OJhqf8ltNIK31xE+E2uc8wUsFEI0MJR7R/PkGXiUZaiVdym+ezZUQMIxVMNDKU2DyGrT2VXuIlwvMtp2EYKWGikaFEk+AesUR4uomXCPd5LKdhGKlgopGhRL2LqpICK7lNM8F4iXCfxzZhMowUMNHIUKKiMaM0n56+sG35mkZswULDGD0mGhlKT58TkqouKwAsRJVO+j0N78C1p0w0DCM5JhoZSjQkVV0aFQ0LUaWL+GtPWXjKMFIhle1e54jIkyKyTUS2isjnXfuXReSAiLzsPt4T0+YmEWkUke0isjLGfr6IbHE/+6677Svu1rD3uvaNIjI3ps1qEdnhPlaTI0RFIupp2PpT6SMYiuDzCB6P9Nt8Xumfv2EYxvCkst1rCPiCqv5FRMqAzSKywf3sdlX9ZuzBIrIYZ7vWM4BZwB9FZKG75euPgDXAC8AjwOU4W75eB7Sp6nwRWQXcBnxERKqAm4GlgLrf/ZCqto3tZ2c+0fDUjNKoaJinkS76wpEB+Qxwl0a38JRhJCWpp6Gqh1T1L+7rTmAbMDtBkyuBe1Q1oKq7gUZgmYjUAeWq+ryqKnA3cFVMm3Xu6/uBS10vZCWwQVVbXaHYgCM0Wc8QT8NyGmkjGBoqGs4qtxaeMoxkjCin4YaNzgU2uqbPisirIrJWRCpd22xgf0yzJtc223092D6gjaqGgA5geoK+sp6e/uopS4Snm2A4MiCfAY5ohCNqVWqGkYSURUNESoFfAf+oqsdxQk2nAkuAQ8C3oofGaa4J7KNtEzu2NSLSICINLS0tCX/HVKG/5LYsH7DwVDoJxPE0fF7nT80qqAwjMSmJhojk4QjGL1T1AQBVPaKqYVWNAD8BlrmHNwFzYprXAwdde30c+4A2IuIDKoDWBH0NQFXvUNWlqrq0uro6lZ+U8XQHQ+T7PJQX5jnvLRGeNoKhoZ5G9L2JhmEkJpXqKQHuArap6rdj7HUxh30AeM19/RCwyq2ImgcsADap6iGgU0SWu31eCzwY0yZaGXU18ISb93gMWCEilW74a4Vry3p6gmGK872U5Du1CrY8evqIlwjP6/c0LDxlGIlIpXrqQuBjwBYRedm1/RvwURFZghMu2gN8CkBVt4rIfcDrOJVXN7iVUwDXAz8DinCqph517XcB60WkEcfDWOX21SoiXwFedI+7RVVbR/dTpxb+YJjiPC/FBV7nvXkaaSNuItx9HzJPwzASklQ0VPVZ4ucWHknQ5lbg1jj2BuDMOPZe4Jph+loLrE02zmyjJximKN9LntdDvs9jnkYaGS4RHv3MMIzhsRnhGYo/GKKkwNH0knyv5TTSSKAvQmGed4DNwlOGkRomGhmKPximyL2wFef7TDTSSG8oPGADJjjhaVgi3DASY6KRofjdRDhAeVEex3v7JnlE2UOgL0JBnomGYYwGE40MxR8MUexWTk0vyedYd3CSR5Q9BEIRCn0Dw1MnSm4tPGUYiTDRyFCiiXCAypJ82kw00kZvX9g8DcMYJSYaGUpvKNKf0zBPI70EQhEKBnka/TPCbaVbw0iIiUaGEug7kaytKsmnszfUvw+EMTYCCRLhVnJrGIkx0chQAqETydqqEmf9qTa/eRtjRVXp7YtQkBc/pxGynIZhJMREIwMJR5RQRPtDKNNd0TjWZaIxVqKexBBPw2cLFhpGKphoZCDRMFRseAqg1fIaYyYQGkY0LDxlGClhopGBBELOkiH5g0XDwlNjJtDniMKQGeEeK7k1jFQw0chATtwNOxe2ftHoCkzamLKF3j5HkC08ZRijw0QjA4neDUcvbNOK8xGx8FQ66BfkIWtP2Sq3hpEKJhoZSDQ8Fa2e8nqEymKbq5EO+s/tsDkNC08ZRiJMNDKQweEpcEJU5mmMnd7hchq23athpISJRgYSFY3YjYKqzNNIC8k8DZsRPjKa2vw4m2wauUIq273OEZEnRWSbiGwVkc+79ioR2SAiO9znypg2N4lIo4hsF5GVMfbzRWSL+9l33W1fcbeGvde1bxSRuTFtVrvfsUNEVpMDxLuwmaeRHoYrufV5zNMYKQfae3jnN57i968dnuyhGBNIKp5GCPiCqi4ClgM3iMhi4EbgcVVdADzuvsf9bBVwBnA58EMRicYCfgSswdk3fIH7OcB1QJuqzgduB25z+6oCbgYuAJYBN8eKU7YS78JWVWqLFqaDQH/11MDwlIiQ7/XQF7G75lTZeqCDcER543DnZA/FmECSioaqHlLVv7ivO4FtwGzgSmCde9g64Cr39ZXAPaoaUNXdQCOwTETqgHJVfV4df/buQW2ifd0PXOp6ISuBDaraqqptwAZOCE3WcqJ66sSFbXpJPm3+IBG7qI2JqCAX5g3908/zioWnRkBjSxcA+1v9kzwSYyIZUU7DDRudC2wEalX1EDjCAtS4h80G9sc0a3Jts93Xg+0D2qhqCOgApifoK6sZXD0FTngqotDeY5sxjYV+QR6UCAfweT0WnhoBjUcc0dhnopFTpCwaIlIK/Ar4R1U9nujQODZNYB9tm9ixrRGRBhFpaGlpSTC0qUF0GZHoInpwYoLfMZvgNyaGS4SDkwy3ktvUiXoaJhq5RUqiISJ5OILxC1V9wDUfcUNOuM/Nrr0JmBPTvB446Nrr49gHtBERH1ABtCboawCqeoeqLlXVpdXV1an8pIzmxAS0E/889ZVFgP0HHSvDldwC5HvFJvelSCSiNDZ34fUIzZ2B/pn2RvaTSvWUAHcB21T12zEfPQREq5lWAw/G2Fe5FVHzcBLem9wQVqeILHf7vHZQm2hfVwNPuHmPx4AVIlLpJsBXuLasJt48jVNmlAKwq6V7UsaULST0NHwWnkqVgx09+INhzj/ZqUtparObmVzBl8IxFwIfA7aIyMuu7d+ArwH3ich1wD7gGgBV3Soi9wGv41Re3aCq0duQ64GfAUXAo+4DHFFaLyKNOB7GKrevVhH5CvCie9wtqto6yt86ZYh3YassyaeqJJ9dR7sma1hZQSAUwSMnSmxjyfN6bMHCFGlsdv4OLz29hk27W9nX6md+Tdkkj8qYCJKKhqo+S/zcAsClw7S5Fbg1jr0BODOOvRdXdOJ8thZYm2yc2cTgtaeinFpdws5m8zTGQm9fmAKfF3eK0AB8HrGl0VMkKhrvOr2Grz76BvuOmaeRK9iM8AwkGI6Q7/UMubCdMqPUPI0xEghF4pbbgjMD38JTqdHY3EVVST7za0opyvOyv61nsodkTBAmGhlIoC8SN+Z+ak0JR7uCdPit7Ha0OOd2aBIcnPCUbfeaGke7gtSWFyIinFRVbAUaOYSJRgYSCIUHVE5FiSbDd5q3MWp6hzm34Ezus/BUaviDIUryHfGdU1VsE/xyCBONDCQQin83fGqNKxrNJhqjZTgvDqKJcBONVOgOhikucFKis6YVcqijd5JHZEwUJhoZSCAUGbDCbZQ5lUXkeYVdRy0ZPloCoXDcORpgojES/IETnsa04nyO9/YRtiVucgITjQwkGArHvRv2eT3MnV7C09tb6AqEJmFkUx/Hixs+PNUXsgtfKviDYYrzHU+jqjgPVeiwJW5yAhONDCTRhe0f372Q7Uc6+Zs7N9ITtFm4IyVachuPPK+Hvoh5GqngD4YoKXDOY2V0D3tbhTknMNHIQBJV+Lz37Dq+dc05vLK/nWd2TP11tiaahCW3Fp5Kme4YT6Oy2BGNdr+JRi5gopGBDFc9FeXdi2sB2GEJ8REzXJEBuAsW2tLoSekLRwiGIv05jSrzNHIKE40MJBCKDFjhdjClBT7qK4ts85tR4ISn4p/bgjxP/7pfxvD43bBotHpqWnEeAO02fygnMNHIQIKhSEJPA+C02jLeNNEYMYEE57Ywz2urtaaAP+gUYQzxNCw8lROYaGQgiUIoURbOLGNnS5eFU0ZIIEEivNDnobcvgrPAsjEc3QFHWItc0SjK85Lv89h2xDmCiUYGEhim5DaW02rLCEWUPcdszsZI6E3kabgXQQtRJeaEp+GEp0SEqmJnO2Ij+zHRyEASldxGOW2mswz1dgtRpYyqOqG/YT0Nx24hqsREPY3ighPnsbIkn9Zuy2nkAiYaGUigL/6M8FhOqS7B6xETjREQ9SCGK7mNzhSP7u5nxKenb6CnAVBZnGcltzmCiUYGEgwnz2kU+LzMm1FiFVQjIN6OiLFExaTHPI2ERD2NksGeholGTpDKdq9rRaRZRF6LsX1ZRA6IyMvu4z0xn90kIo0isl1EVsbYzxeRLe5n33W3fMXdFvZe175RRObGtFktIjvcR3Q72KwmFI4QjmjS8BTAgppSGptNNFIl0Df8Vq/gJHTBwlPJiOY0iod4GhaeygVS8TR+Blwex367qi5xH48AiMhinK1az3Db/FBEorcjPwLW4OwZviCmz+uANlWdD9wO3Ob2VQXcDFwALANudvcJz2r674aTlNwCzK8pZV+rv397WCMxJzyNZOEpO5+J6Pc0YkSjqjifdn+QiC1amPUkvTKp6p9w9u1OhSuBe1Q1oKq7gUZgmYjUAeWq+rw69Yx3A1fFtFnnvr4fuNT1QlYCG1S1VVXbgA3EF6+sInphSzS5L8r8mlIiCnuO2l4GqRAV1+FWuY0KteU0EhP1NKIlt+CsdBtRON5r3ka2M5acxmdF5FU3fBX1AGYD+2OOaXJts93Xg+0D2qhqCOgApifoK6uJXtgKhrmwxXJqtbO/RqMtJ5ISvcPsvR7FPI3U6A6Gyfd6BhRr2FIiucNoReNHwKnAEuAQ8C3XLnGO1QT20bYZgIisEZEGEWloaZnai/gFklzYYjm1uhQRE41USSbIVnKbGv5AaEC5LZxY6dbmamQ/oxINVT2iqmFVjQA/wck5gOMNzIk5tB446Nrr49gHtBERH1CBEw4brq9447lDVZeq6tLq6urR/KSMIbrdaLLqKXDCA/WVRTS2mGikQjJPIxpu6bUcUUK6g2GKBwlvpbv+VJvN1ch6RiUabo4iygeAaGXVQ8AqtyJqHk7Ce5OqHgI6RWS5m6+4Fngwpk20Mupq4Ak37/EYsEJEKt3w1wrXltWMxNMAmF9dap5GivQvtJefuOTWchqJ8QdD/YsVRokuj25lt9mPL9kBIvJL4GJghog04VQ0XSwiS3DCRXuATwGo6lYRuQ94HQgBN6hq9LbtepxKrCLgUfcBcBewXkQacTyMVW5frSLyFeBF97hbVDXVhPyUJRpCSTa5L8r8mlL+vPMYzcd76Ysos6cVjefwpjTR+RfDioaFp1KiOxDuX6wwyvRSRzSOdZloZDtJRUNVPxrHfFeC428Fbo1jbwDOjGPvBa4Zpq+1wNpkY8wmkpWFDmZ+TSmBUIR3fuMpppfm86d/vgSPJ146yOjpr/qJ/2cfTYTb5L7E9MRswBSlON9Hcb6Xo12BSRqVMVHYjPAMYyTVUwCL6soB5+65qa2HV5rax21sU53+8NRwJbc+C0+lQnfMVq+xVJcVmGjkACYaGUZwhJ7G2fXT+NX1b+X3//gO8rzCI1sOjefwpjRR0SgaJjzl8QgFPk//zHEjPv44ngbAjNICWjpNNLIdE40MY6ThKYDzT66iuqyAi+bP4JEth20/iGHoCYbxSOJzaxsxJac7MIynUWqeRi5gopFhRKunUk2Ex3LFWXUcaO/h1aaOdA8rK/AHwxTleXGXPYtLYZ7HchpJGNbTKMs3TyMHMNHIMPpzGinM0xjMOxY4c1Re3m95jXj09IWGTYJHKcrzWk4jAarq5DTihPiqSwtp8/fRF7bzl82YaGQY/RPQUliwcDC15QUU+Dzsb7W1qOLh3CEnFmMLTyXG2Q6XIfM0wPE0wMpusx0TjQwjmqwtSXJHHA8RYU5VMfvbTDTikYpoFOR56bXtXoelu39Z9Pg5DcBCVFmOiUaG4Q+GKPB58I5yrsWcyiL2tfakeVTZQW9feNjKqSiFPg+9QfM0hsMf3eo1bk7DEQ1Lhmc3JhoZhlMDP3IvI8pJVcU0tfqtgioOqXgaRfleW3sqAf4+8zRyHRONDMMfSH5hS8ScqmI6AyE6emzhuME41VOJBbnQZzmNRBzvcUSjvDBvyGfVrqfRYp5GVmOikWE4lSmj9zTqK4sB2GfJ8CH0BEMpJMI9Vj2VgHZ3QcJpxUNFozDPS1mBzzyNLMdEI8PwB8ND9ioYCSdVOaKx3/IaQ0i1esrmaQxPu+vBVhQNFQ1wvA3zNLIbE40MozswNk9jTpWzyq1VUA2lJ5hCItxKbhNy3BWNeJ4GOEuJHHU9jUhE6bTtX7MOE40MI5W74USUFeZRWZxn4alBqCr+PmdGeCIK87z9s/KNobT7+/B6hNJhijViPY0f/2knb/vaExw53juRQzTGGRONDGOsogFOMtwm+A0kGI4QjmhKOY3oscZQ2nuCVBTlDbsUS3VZAYc7eukKhPj583vp7A3x34/vmOBRGuOJiUaGEW9XtJFiojGUnv4VbpNUT+XZRkyJaPf3MW2YfAbAe8+uwx8M86n1DRzs6GVhbSn3vrifHz+9k19s3DuBIzXGi6SiISJrRaRZRF6LsVWJyAYR2eE+V8Z8dpOINIrIdhFZGWM/X0S2uJ991932FXdr2Htd+0YRmRvTZrX7HTtEJLolbFYTb1e0kXJabRl7W/10BUJpGtXUJ9mufVEK+/fUMNGIR0dPHxXD5DMA3jK3issW1/Jc4zGqSvJZ94llFOd5+dqjb/B/f/0azZ0WqprqpOJp/Ay4fJDtRuBxVV0APO6+R0QW42zXeobb5ociEv1f+iNgDc6+4Qti+rwOaFPV+cDtwG1uX1U4W8teACwDbo4Vp2wkHFF6+uKvIDoSzppdgSpsPWCr3UZJtj94lGii3JYSiU9HT9+wlVNRbrzidHwe4Zrz66mrKOLxL7yT/7rK2bTzaKetSzXVSSoaqvonnL27Y7kSWOe+XgdcFWO/R1UDqrobaASWiUgdUK6qz6szVfnuQW2ifd0PXOrXsPfuAAAgAElEQVR6ISuBDaraqqptwAaGildWEb0bjrdXwUg4c3YFAFtMNPrpD0+lkAgH8zSGI1l4CuDU6lI2/J938k+XLQSgpryQhbVlALR2m2hMdUZ7S1urqocAVPWQiNS49tnACzHHNbm2Pvf1YHu0zX63r5CIdADTY+1x2mQl/v7F4MbmaVSXFTCrotD21YjhhKeR+NxGl6TvsfWn4tLuDzKtOD/pcfNmlAx4X1XiroDbbXM4pjrpToTHK6nQBPbRthn4pSJrRKRBRBpaWlpSGmgmEl0MbqyeBsBZ9RXmacQQFeTk8zSc/xIBW39qCOGI0hkIJQ1PxWNGqS2bni2MVjSOuCEn3Odm194EzIk5rh446Nrr49gHtBERH1CBEw4brq8hqOodqrpUVZdWV1eP8idNPtFlp5Otj5QKZ9dPY/fRbluDyqUn1ZxGf3jKchqD6eztQ3X42eCJKC/Mw+cR8zSygNGKxkNAtJppNfBgjH2VWxE1DyfhvckNZXWKyHI3X3HtoDbRvq4GnnDzHo8BK0Sk0k2Ar3BtWUv/Xhrp8DTcvIYlwx1STYRbTmN42v2JZ4MnwuMRKkvyzdPIApLe0orIL4GLgRki0oRT0fQ14D4RuQ7YB1wDoKpbReQ+4HUgBNygqtH/fdfjVGIVAY+6D4C7gPUi0ojjYaxy+2oVka8AL7rH3aKqgxPyWUV3ID05DYAzZpUD8MbhTt42f8aY+5vq+PtGlgi39aeG0p5kCZFkTC/J55glwqc8Sa9OqvrRYT66dJjjbwVujWNvAM6MY+/FFZ04n60F1iYbY7aQTk+jqiSfAp+Hw7aEA+CscAup5zQsPDWUjv7FCpMnwuMxvTSfY7aY4ZTHZoRnEFFPYywLFkYREeoqCjnUYaIB0BN0RCCZF2fhqeGJLos+mpwGwPSSAiu5zQJMNDKIVGctp0pteSFHTDQAZ8e5/BS20S0rdESls9dm0w+mY6zhqVLLaWQDJhoZRHd/ye3YPQ2AmRWFHDpu+2qAUz2VihgX5/sozPPQalU+Q4gmwkfvaeTTGQhZOfMUx0Qjg/AHQ4hAgS89/ywzKwo50hGw/cJxVw9OkgSPUlWcT2u3lSoPpqOnj9ICH3ne0f19Tnf3ELcQ1dTGRCODcBYr9A277PRImVleSDAcsf+kpLYBU5TKknza/HbOBtPuT77uVCL6Z4VbiGpKY6KRQfhT2MN6JNRVFAJYBRXRc5ta2K+qJN+ENg4d7l4ao6V/Vrid2ymNiUYG0R0Mpy2fAU4iHOCwJcPxj8TTKDZPIx5HjgeoLisYdfuqEqetld1ObUw0Mgh/IN2ehrNfuHkaTjXUcFuUDqaqJJ9WC6EM4UB7D7Mri0bdfrrraZgXN7Ux0cgg/MFwWuZoRJlRmo9HzNMAONjR0x+uS0aVW+UTtD01+vEHQ7R2B6kfg2iUFfjI93o4aoI8pTHRyCCcrV7T52n4vB5qygpzXjS6AyHa/X0p3yVXugnbdgtR9XOgzSndnj1t9KIhIkwvzeeohaemNCYaGUR3inMJRkJtRWHOh6cOtI/sglfl7hfRaqLRT5MrGmPxNABqygpo7jTRmMqYaGQQTk4jfeEpgJnlBTnvaURFI9ULXmWJUyFksfcTNPWfw+Ix9VNdVkhzjt/ETHVMNDKI7mCYkjR7GnUVRRxo7yESyd0JfidCK6ld8KaX2CS0wRxo6yHf66G6dPTVUwC15eZpTHVMNDKEcETpCoQoKxx9HXw8zppdgT8Y5s3mzrT2O5U40N5DnleoSbFcNOpptJlo9NPU5qduWiGeJGt3JaOmrJDW7qAVGUxhTDQyhObOXsIRpW5aahU+qbJ0biUADXva0trvVOJAWw91FUUpX/AqozkNW0qknwPtPWPOZwDUlDvC3WLJ8CmLiUaGcNCNGc8aQ3VKPE6qKmZGaQGb9+awaLT3jKjqJ8/roazQZxP8YmhqG9k5HI5aVzQsrzF1MdHIEA62O/+JZlWkVzREhLfMreTFPVm96WFCDrT1jFiMp9tSIv309oVp6QyknBNKRE2Z40kfOW6exlRlTKIhIntEZIuIvCwiDa6tSkQ2iMgO97ky5vibRKRRRLaLyMoY+/luP40i8l13H3Hcvcbvde0bRWTuWMabyZzwNNIbngI4/+RKmtp6OJKDd3d94QhHOntHPJO50kSjn4MjrD5LRH94qjP3/hazhXR4Gpeo6hJVXeq+vxF4XFUXAI+77xGRxTj7f58BXA78UESipUI/AtYAC9zH5a79OqBNVecDtwO3pWG8GcnB9h7KCn1pT4QDvGVuFZCbeY3DHb2oQv0IPQ1neXQTDYC9rX4gPaIxvaQAj5inMZUZj/DUlcA69/U64KoY+z2qGlDV3UAjsExE6oByVX1enY0f7h7UJtrX/cClkq51wzOMgx29aQ9NRVk8q5x8r4dXD7SPS/+ZTHRS2mg8DctpOGw/7FTenT6zfMx9eT1CdVkBzeZpTFnGKhoK/EFENovIGtdWq6qHANznGtc+G9gf07bJtc12Xw+2D2ijqiGgA5g+eBAiskZEGkSkoaWlZYw/aXI42N4zLqEpcBK7c2cUs6ule1z6z2S2HuwA4NTq0hG1qy4r4GhXgFDYSkO3HTrOrIpCKka5zetgasoKzdOYwoxVNC5U1fOAK4AbROQdCY6N5yFoAnuiNgMNqneo6lJVXVpdXZ1szBmJIxrj42mAc9Hc2dw1bv1nKpt2t3Ly9GJmprhYYZRTq0vpC2t/aCaXeeNQJ6fXjd3LiGJLiUxtxiQaqnrQfW4Gfg0sA464ISfc52b38CZgTkzzeuCga6+PYx/QRkR8QAWQdWVAPcEwbf6+cReNva3+nJpUFYkom/a0sszN6YyEhbWOZ7LjSO5OigQIhMLsbOliUV1Z2vqsKZ86S4l09PTx/Sd28Hc/3cS6P+/J6ZUVooxaNESkRETKoq+BFcBrwEPAavew1cCD7uuHgFVuRdQ8nIT3JjeE1Skiy918xbWD2kT7uhp4QrNww+uDHeNXORXl1JoSwhFlX2vuhKjebO6k3d/HBacMiWgmZX6NIxpvHsk97yyWxuYuQhFNSz4jSk1ZAce6g/RNgdDffz68lW/+4U3ePNLFzQ9t5bp1L5KFl6ARMRZPoxZ4VkReATYBv1PV3wNfAy4TkR3AZe57VHUrcB/wOvB74AZVDbt9XQ/ciZMc3wk86trvAqaLSCPwf3ArsbKN/nLbcUqEw4mYfmOOhKh6gmE27nKc0gvmjdzTKM73UV9ZxJs57mlsO+T8/kVpDE9Fd5TM9BBVTzDMY68dZtVb5vDsv17Cp955Ck9ub8n5VaNHvaSqqu4CzoljPwZcOkybW4Fb49gbgDPj2HuBa0Y7xqnCoejEvnEOTwHszIFk+GNbD/Ppn2+moiiP2dOKmFM1uklpC2vLckZkh+ONQ8cp8HmYN6MkbX2eUu30teNIZ1pmmY8XG7YdoTsY5solsxERLj29lv95ehfbD3f274qZi9iM8Axg97FuvB7pvwMbD0oKfNRVFOZEMvyx1w5TlOfFK8LKM2aOup8FtaXsaunO6QqqrQePc9rMMrxjXKgwltNnOvmRNw5nthf34EsHmFle2O+pRvNcue59mmhkAA17WjlrdgX5vvH95zi1upSdLdktGqrKs41HedfpNWz+98v4j/ctHnVfC2vKCIYj7DmWGxVUzZ29fP+JHf3FEq3dQV7c08rbTp2R1u+ZVpxPXUUhbxw6ntZ+08nRrgBPv9nC+5fM6l/oclpxPrXlBWw/nN3/h5JhojHJ9PaFeWV/x6ji7iPl1OoSGpu7sroCZGdLF82dAS6cP/YL3cJa5444Vyqo7npmN9/8w5v89+NvAvDIlkOEIsr7z5mV9u9aVFee0Z7GfQ37CUWUDy+tH2BfWFtmnsZkDyDXeXl/O8FwpH+pj/FkyUnT6A6GefVAx7h/12Tx7I6jAFyUBtGYX1OK1yNZfb6iRCLKQ68cxOcRfvTUTjbuOsZDLx9kfk1pWstto5w+08kXZWIJeDii/OKFfbz1lOnMrxn420+rLWNHcyfhCbzxUlVuefh1HvhLU/KDJwATjUlm0+5WRJgQ0bjktBq8HuGPrx8Z9++aLJ5tPMZJVcWjTn7HUpTvZfkpVTy29XDWl1lu3N3KoY5evnLVmcyuLOKjP3mBTXtaufKcWYzHyj2n15UTimhGhkuf2t7MgfYePvbWk4d8tnBmGb19EfZP4KTPP25rZu1zu/nXX73KaxlwA2OiMcls2t3KabVlaVuiIRHTivNZenIlf9yWnaLRHQjx551HuWhB+mLwl59Zx66WbnZkeQHBgy8foCTfy1VLZvObz1zI6rfN5eTpxXzgvNnJG4+CRf3J8MnPa+w52j3gpmD9C3upKSvgssW1Q449zQ1Zbh9DiEpVeXFPK7194aTHhsIRvvboNuZOL2Z6SQGf++VLKbUbT0w0JpG+cITNe9tYNgH5jCiXLa7ljcOdE3qnNFE8suUQ/mCYD5ybvgvdyjNqEYFHtxxOW5+ZxsZdx3jgpQNcfmYdRfleppcWcPP7zuDpf76E+sqxe2zxmDejhHyvp38eyGSxcdcxLv7mU9zzorMs3t5j3Tz9ZgsfXXYSed6hl8fopM+x3PH/fOM+rvnx83z655vpC0cShrrua2hiZ0s3N16xiG99+Bx2He3mzmd2jfq704GJxiTy8v52evrCvHUUM5ZHy6WLnLunP2RhiOp/Nzcxb0YJS0+uTH5witSUFbL05Eoefe1Q2vrMJF4/eJxPrmvgpKpi/u97F03Y9/q8HpacNI1Hthya1LzGXc/uBuA7f3yT3r4w/2/jPjwifHTZSXGPLynw8fYFM1j/wl46eoZuB6yqNLWduCFr9wf5c+PR/jDclqYOvvLw65wyo4SntrdwyTef4rQvPco7v/Ek//nwVg64E33B8Zxv/+ObLD25kpVn1HLh/BlcfsZMfvjUzkndG8dEYxJ5rvEoIvDWUydONObNKOGcOdO4+/k9WTP/oLmzl8e2HmbT7lauPr8+7TH49y+ZzRuHO3lmx9RcQXk4WjoDfHLdi5QW+lh/3TKqSvIn9Puvv/hUmtp6Ji3Bu7/Vzx+3HeHC+dM5cjzAF+57hXte3M+KxbUJF7i86YpFdPT08cMnG4d8duczu7notid5bOthnms8yrJbH+ev79zINT9+nv2tfj53z0vMKM3nV9e/jS+/bzH1lUV8/G1zWVBTys9f2MvF33iS9S/s7e+rpTPATe9Z1P83fdN7TicUVtas3zxp+Q0TjUnkz43HOHNWBdOKJ/Y/62cuPpW9x/z8bsvUv3ve0tTBxd94ik+t30y+18MHxyEG/+Gl9cyeVsRtv38ja8qVO/x9rFnfQKs/yE+uXTopM5wvXljNOXOm8b0nGmmfhL1L1r+wFxHhm9ecw7sX1fC7LYeYVpzHZ981P2G7xbPK+dB59fz0uT0Dwrz7jvn51obtAHzlt69z0wNbmF1ZxH+vWkJXIMR7v/sMu492880Pn0NlST4fv3Ae96x5K1/6q8XcufotPP3Pl3DR/Bn8+29e4x9++RLfe2IHV5w5k/NjPOeTp5fwrQ+fw/5WP+///rP8ctO+8Tk5CTDRmCT8wRAv7W/jbfMnzsuIctmiWhbWlvKDJxundFXQoY4erlv3IpXF+ay/bhlPfPGd43LxK/B5+cKKhbx24DgPvnIg7f1PNAfbe/jQj//Mawc6+M5HzuXM2RWTMg4R4V9Wnsahjh7e8fUnefDliTu3wVCE+zc3cdmiWuoqivjB35zHS/9+GU//8yWcMSv5+fjCioV4PPDNPzgiEY4oNz7wKj6Ph29/+Bya2nrY1+rnqx88iyuXzOaLKxZyvDfEx982d9jJkrOmFfGjvz2fC+ZV8fArB7nq3Nl87UNnDznufefM4skvXsw7FlZz0wNbWOuG2CaKUa89ZQzlxT2t/MeDW/n0O0/hyiXx73hVld++eog9R7vpC2ta5hOMFI9H+MSF87jxgS1sO9TJ4lnpW4xuolBVvnDfK3QHQjzwmQs5bWb65xLEcuWS2dz9/F7+/TdbObt+2og3dcokvvzQVg6197D+ugtYPoH5tHhcOH8Gj3z+7dz0wBb+5f5XOWt2BadMwLl9cnszrd1BPvwWZ/Jegc9Lgc+bpNUJ6iqK+ORFp/D9Jxv54Hn1vLi7lT/vPMZtHzqLD55Xz47mLorzvP3n95MXncKiuvKkRS+FeV7WfWIZ+1v9LKgd/m+6oiiPn1y7lBt+8Rf+v0e28Za5VZxVPzHiL1P5TjMeS5cu1YaGhgn/3p0tXXzwh3+mOxAiFFHOPWka82aU8A/vWjBgsbc/vn6ET97tjK/A5+Hl/1hBUX7qf6zp4lBHD2/96hN86b2L+OTbT5nw7x8LqsovN+3n3369hVs/cCZ/c8HQevrx4GB7D+/73rNUluTzu89dNKKLTKbwyv52rvzBc/yfyxbyuUsXTPZw+jlyvJd3f/tpFtWVc8/fL+9fumO8+Pu7G3h5fzvP3/gufHGqpFKhKxDi0m891b8L4YeX1vP1q4es4TqudPj7WPmdP1GY5+GM2RVUFufxX1edNaq+RGSzqi5NdpyFp8ZIJKL89LndfOAHz+HzCH/4p3fwxRUL8XmEDVuPcNUPnuP5nccAp+b6q49u45QZJfzy75dzz5rlkyIY4NwpnVJdwrONRyfl+0fLnc/s4rQv/Z5/+/UWlp9SxUffEr/KZTyYNa2I2z50No3NXTz8ytTJB3UFQkQiSjiifOOx7VQW5/GJi+ZN9rAGUFteyJfeu4hNu1v59M83c7x3aGVSOnjzSCc/eLKRJ99o5oPnzh61YACUFvh45HNv5+sfOpvPXjKfW64cslD3uFNRnMc3rjmbA+09vNrUPiHfaZ7GGPnlpn3c9MAW3r5gBl9+/xkDwhb7jvn5xLoX2XO0my+9dxG7jnZz9/N7ueNj57NiDKuvpov/ePA1/rehiVduXjHuiyUmozsQwusRVJ1qqNcOHKfNH6S6rIC3nTqdssI8Xm1q5wM//DPLT6niktNq+OB59RNe8aOqXPHfzwDw6OffPi6zpUdLR08f/2/jProCfax6y0l0B0P87Lk93Nuwn0Uzyykv8vHCrla+/L7FfPzCzBINcM7tXc/u5quPvsH86lLuWbOcyhT+fXe2dPGTP+3iQHsPFUV5LJtXxd9ecPIQb2VnSxdXfv85ugIhTplRwrpPLEvLygGZQDiiY16JOFVPw0TDRVU53huivNCX8oUgElEuu/1pivK9PPzZi+K2O97bxw2/+AvPuGsiffC82XzrmnMy4mLz2NbDfGr9Zu5ds3xUu9ulg0hE+c7jO/j+EzsYrjCpON/LsnlVvOFOBHvsn95BRdH4z6Afjvs3N/HF/32Fn/7dW7jktJpJG0cUVeX+zU3c8vDrdAZCeIT+c+nzCB86r57ndh6lpTPAV646kw8vnZO4w0nmmR0tXLeugdNnlrH+ExckXC2hOxDir773LEeO97KgppSjXUEOtPfwxRUL+dvlJ/Pk9mb8wTDBUISfv7CXNn8fv/nMhZw0PTvEIp1klWiIyOXAfwNe4E5V/dpwx45WNNr9QZbcsoE8rzBvRgkfXjqHC+fPYO70kmFDSE9tb+bjP32R2z9yDh84tz7uMeCEpX776iGWzJnG3DRuZjNWOnr6OPeWP3Dlktl84+qz6Q6EKSnwJnTZVZVndhzlwZcPMnd6McvmVXHazLJRlQ2/eaST/3x4K881HuP958zitJlliMCMkgIW1ZVTXVbAvlY/92/ez5YDxynK83DTexZNyDpdiQiGIlzyzado8wf5l5Wn8Z6z66gpO1HXr6oTdlOw91g3t/5uG394/QgXzKviP963mGnF+Tz8ykFqywtYNm86s6cVEQxF6AqEJtwzGy1PvHGET63fzJzKYr6wwqmwOu/kSs6dMw2AV5s6eGlfG882HuPxN47wy79fzvJTpqOq/NO9L/PgKwcpzffRGQj191mc7+XO1UvTvtR7tpA1oiEiXuBNnK1jm4AXgY+q6uvxjh+taHQFQtyzaR/HuoNs3HWMv+w7ER+cVVHIvOoSassLCYQi9ATDdAdC7DrajQDP/uu7Jj28M1r+67evc+ezu6mrKOTw8V7mzSjhax88m9++epDdR7vxiPD6oePUVRTymYvnc+czu2jY20ZZgY+uYIjon8+lp9fw1xecRGdviOAwkwbrKgo5tboUfzDE2uf2cM+mfZQW+PjXK07nr5edlBHeV6ocbO/hX+5/tT8ntGxeFR9ZOoc7/rSLPce6qSkvoLaskJryAmrc50Uzy7lowQx6+sLsbummpTPAawc76O2LcNW5s2jtDrL9cCdFeV5a/UG6AyGWzKmkpqyAsCqRiFJbXkhlST6/fukAj712mBd2HSPP6+Fzly5gzTtOSetmSZNNNL/R2n1iDkdlcR6FeV4OdZyYEf1P717I5999IqnfEwyzZn0DpQU+PvXOU5lZXki+z0NxvpfCvKlXvDBRZJNovBX4sqqudN/fBKCqX413fLpyGo3NXWw7dJzdR7vZfbSbXUe7aT7eS1G+l+J8L8X5PorzvXxk6RyuOKtuzN83mTz48gHu39zEGbMquK9hP63dQfK8wuK6cvrCymkzy3iu8SjNnQEqivK48YrT+eB5s+nqDfHqgQ5e2tvGT/+8h87eUPIvc/F5hL9dfjKfv3RBSnHrTERVeWl/Oy/sOsadz+ymtTvISVXFrFhcS0tXgObjAY509tJyPNB/x1ta4KM7RmxFwCtCKE5szuuRuOsSRcNPC2pKedeiGj5x4bxx3fVxMjnWFWD30W7qK4v505stvNzUTldviGXzqlixuJaifC9lhZMXqswmskk0rgYuV9VPuu8/Blygqp+NOWYNsAbgpJNOOn/v3r2TMtZsYH+rn3tf3M81S+s5efqJUFpHTx+/eekAV5w5k5o4F6gOfx9vHD7O9NICCvOGel2qsK/Vz55j3RT4vCw9uTKjQnVjpcPfR8PeVi6cPyPu3ayzAu8xHt92hFnTijh9ZhkzygpYUFNKIBThkS2HmFleyLknVdIXjlBRlIfXI2w50MHxnj68HsEjwt5WP02tfq44q44lbqjGMNJBNonGNcDKQaKxTFX/Id7xkzVPwzAMYyqTTfM0moDYco964OAkjcUwDCOnmQqi8SKwQETmiUg+sAp4aJLHZBiGkZNk/NpTqhoSkc8Cj+GU3K5V1a2TPCzDMIycJONFA0BVHwEemexxGIZh5DpTITxlGIZhZAgmGoZhGEbKmGgYhmEYKWOiYRiGYaRMxk/uGyki0glsH2GzCiBdu7Sns69YZgDp2Pwi3eMbr9+bzn7Tde6iZPo5THd/dv7GRjrP33j91hlAiapWJ22hqln1ABpG0eaONH5/2voa6++aiPGN4+9N579JWs7dVDmH49Cfnb8MOX/j9VtHMkYLTzk8nKF9jQfpHt94/d5MPo+Zfg4z+dxB5v/eTD5/k/5bszE81aAprJ8y1cjW3zUR2LkbG3b+xsZUOH8jGWM2ehp3TPYAxols/V0TgZ27sWHnb2xMhfOX8hizztMwDMMwxo9s9DQMwzCMccJEY5IQkTki8qSIbBORrSLyeddeJSIbRGSH+1zp2qe7x3eJyPcH9fVREdkiIq+KyO9FJKs3QU7zufuIe962isjXJ+P3TDSjOH+Xichm929ss4i8K6av8117o4h8V6bSnr2jJM3n71YR2S8iXZP1e0ZMOsu37DGiUrc64Dz3dRnOPuiLga8DN7r2G4Hb3NclwEXAp4Hvx/TjA5qBGe77r+Nsjzvpv3EKnLvpwD6g2n2/Drh0sn9fBp6/c4FZ7uszgQMxfW0C3goI8ChwxWT/vil2/pa7/XVN9u9K9WGexiShqodU9S/u605gGzAbuBLn4oX7fJV7TLeqPgv0DupK3EeJe5dXTpZvUpXGc3cK8Kaqtrjv/wh8aJyHP+mM4vy9pKrRv6mtQKGIFIhIHVCuqs+rcwW8O9omm0nX+XM/e0FVD03k+MeKiUYGICJzce5GNgK10T8i97kmUVtV7QOuB7bgiMVi4K5xHG5GMZZzBzQCp4vIXBHx4fwnn5OkTVYxivP3IeAlVQ3gXCibYj5rcm05wxjP35TERGOSEZFS4FfAP6rq8VG0z8MRjXOBWcCrwE1pHWSGMtZzp6ptOOfuXuAZYA8QSucYM5mRnj8ROQO4DfhU1BTnsJwpx0zD+ZuSmGhMIu4F/1fAL1T1Add8xHX7cZ+bk3SzBEBVd7ohgvuAt43TkDOGNJ07VPVhVb1AVd+Ks2bZjvEacyYx0vMnIvXAr4FrVXWna24C6mO6rSfLQ6NR0nT+piQmGpOEm3+4C9imqt+O+eghYLX7ejXwYJKuDgCLRSS60NhlODHWrCWN5w4RqXGfK4HPAHemd7SZx0jPn4hMA34H3KSqz0UPdkMwnSKy3O3zWlI451OddJ2/KctkZ+Jz9YFTzaM44aSX3cd7cCp6Hse5430cqIppswdoBbpw7vIWu/ZP4wjFqzhryUyf7N83hc7dL4HX3ceqyf5tmXj+gC8B3THHvgzUuJ8tBV4DdgLfx50wnM2PNJ+/r7t/jxH3+cuT/fuSPWxGuGEYhpEyFp4yDMMwUsZEwzAMw0gZEw3DMAwjZUw0DMMwjJQx0TAMwzBSxkTDMCYYEfm0iFw7guPnishr4zkmw0gV32QPwDByCRHxqeqPJ3schjFaTDQMY4S4i9T9HmeRunNxlsa+FlgEfBsoBY4CH1fVQyLyFPBn4ELgIREpw1kK+5sisgT4MVCMM0HuE6raJiLnA2sBP/DsxP06w0iMhacMY3ScBtyhqmcDx4EbgO8BV6tq9IJ/a8zx01T1nar6rUH93A38q9vPFuBm1/5T4HPqrIllGBmDeRqGMTr264l1hH4O/BvOBjsb3M3rvEDsPgn3Du5ARCpwxORp17QO+N849vXAFXhIGu8AAADcSURBVOn/CYYxckw0DGN0DF5/pxPYmsAz6B5B3xKnf8PICCw8ZRij4yQRiQrER4EXgOqoTUTy3P0ThkVVO4A2EXm7a/oY8LSqtgMdInKRa/+b9A/fMEaHeRqGMTq2AatF5H9wVjX9HvAY8F03vOQDvoOzvWciVgM/FpFiYBfwd67974C1IuJ3+zWMjMBWuTWMEeJWT/1WVc+c5KEYxoRj4SnDMAwjZczTMAzDMFLGPA3DMAwjZUw0DMMwjJQx0TAMwzBSxkTDMAzDSBkTDcMwDCNlTDQMwzCMlPn/AXHrgACWpGl1AAAAAElFTkSuQmCC\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": 11, + "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": 12, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +2344,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 13, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAD8CAYAAABQFVIjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/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": 14, "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": 14, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,9 +2446,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 15, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "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,