diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5b58a3ae26346460dd39e62a39c55243d7..82fabd085054b1e8d4871500d165137d16efefdd 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,11 +59,980 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], + "execution_count": 3, + "metadata": { + "scrolled": true + }, + "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
020214334348633073.053899.06650.082.0FRFrance
120214232920423995.034413.04436.052.0FRFrance
220214132504320586.029500.03831.045.0FRFrance
320214032628621842.030730.04033.047.0FRFrance
420213932215518014.026296.03428.040.0FRFrance
520213831561412310.018918.02419.029.0FRFrance
620213731367310404.016942.02116.026.0FRFrance
72021363102897505.013073.01612.020.0FRFrance
82021353126099282.015936.01914.024.0FRFrance
92021343130159485.016545.02015.025.0FRFrance
102021333103927042.013742.01611.021.0FRFrance
1120213231558611009.020163.02417.031.0FRFrance
1220213131885513664.024046.02921.037.0FRFrance
132021303139919695.018287.02114.028.0FRFrance
142021293136269618.017634.02115.027.0FRFrance
15202128386365430.011842.0138.018.0FRFrance
162021273106936838.014548.01610.022.0FRFrance
17202126370864109.010063.0116.016.0FRFrance
18202125379425540.010344.0128.016.0FRFrance
19202124348553011.06699.074.010.0FRFrance
20202123367104455.08965.0107.013.0FRFrance
21202122378795495.010263.0128.016.0FRFrance
22202121378275403.010251.0128.016.0FRFrance
232021203102787540.013016.01612.020.0FRFrance
24202119395396860.012218.01410.018.0FRFrance
252021183121359165.015105.01814.022.0FRFrance
262021173120588891.015225.01813.023.0FRFrance
2720211631650512735.020275.02519.031.0FRFrance
2820211531930615398.023214.02923.035.0FRFrance
2920211432107317099.025047.03226.038.0FRFrance
.................................
190119852132609619621.032571.04735.059.0FRFrance
190219852032789620885.034907.05138.064.0FRFrance
190319851934315432821.053487.07859.097.0FRFrance
190419851834055529935.051175.07455.093.0FRFrance
190519851733405324366.043740.06244.080.0FRFrance
190619851635036236451.064273.09166.0116.0FRFrance
190719851536388145538.082224.011683.0149.0FRFrance
19081985143134545114400.0154690.0244207.0281.0FRFrance
19091985133197206176080.0218332.0357319.0395.0FRFrance
19101985123245240223304.0267176.0445405.0485.0FRFrance
19111985113276205252399.0300011.0501458.0544.0FRFrance
19121985103353231326279.0380183.0640591.0689.0FRFrance
19131985093369895341109.0398681.0670618.0722.0FRFrance
19141985083389886359529.0420243.0707652.0762.0FRFrance
19151985073471852432599.0511105.0855784.0926.0FRFrance
19161985063565825518011.0613639.01026939.01113.0FRFrance
19171985053637302592795.0681809.011551074.01236.0FRFrance
19181985043424937390794.0459080.0770708.0832.0FRFrance
19191985033213901174689.0253113.0388317.0459.0FRFrance
192019850239758680949.0114223.0177147.0207.0FRFrance
192119850138548965918.0105060.0155120.0190.0FRFrance
192219845238483060602.0109058.0154110.0198.0FRFrance
1923198451310172680242.0123210.0185146.0224.0FRFrance
19241984503123680101401.0145959.0225184.0266.0FRFrance
1925198449310107381684.0120462.0184149.0219.0FRFrance
192619844837862060634.096606.0143110.0176.0FRFrance
192719844737202954274.089784.013199.0163.0FRFrance
192819844638733067686.0106974.0159123.0195.0FRFrance
19291984453135223101414.0169032.0246184.0308.0FRFrance
193019844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1931 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202143 3 43486 33073.0 53899.0 66 50.0 \n", + "1 202142 3 29204 23995.0 34413.0 44 36.0 \n", + "2 202141 3 25043 20586.0 29500.0 38 31.0 \n", + "3 202140 3 26286 21842.0 30730.0 40 33.0 \n", + "4 202139 3 22155 18014.0 26296.0 34 28.0 \n", + "5 202138 3 15614 12310.0 18918.0 24 19.0 \n", + "6 202137 3 13673 10404.0 16942.0 21 16.0 \n", + "7 202136 3 10289 7505.0 13073.0 16 12.0 \n", + "8 202135 3 12609 9282.0 15936.0 19 14.0 \n", + "9 202134 3 13015 9485.0 16545.0 20 15.0 \n", + "10 202133 3 10392 7042.0 13742.0 16 11.0 \n", + "11 202132 3 15586 11009.0 20163.0 24 17.0 \n", + "12 202131 3 18855 13664.0 24046.0 29 21.0 \n", + "13 202130 3 13991 9695.0 18287.0 21 14.0 \n", + "14 202129 3 13626 9618.0 17634.0 21 15.0 \n", + "15 202128 3 8636 5430.0 11842.0 13 8.0 \n", + "16 202127 3 10693 6838.0 14548.0 16 10.0 \n", + "17 202126 3 7086 4109.0 10063.0 11 6.0 \n", + "18 202125 3 7942 5540.0 10344.0 12 8.0 \n", + "19 202124 3 4855 3011.0 6699.0 7 4.0 \n", + "20 202123 3 6710 4455.0 8965.0 10 7.0 \n", + "21 202122 3 7879 5495.0 10263.0 12 8.0 \n", + "22 202121 3 7827 5403.0 10251.0 12 8.0 \n", + "23 202120 3 10278 7540.0 13016.0 16 12.0 \n", + "24 202119 3 9539 6860.0 12218.0 14 10.0 \n", + "25 202118 3 12135 9165.0 15105.0 18 14.0 \n", + "26 202117 3 12058 8891.0 15225.0 18 13.0 \n", + "27 202116 3 16505 12735.0 20275.0 25 19.0 \n", + "28 202115 3 19306 15398.0 23214.0 29 23.0 \n", + "29 202114 3 21073 17099.0 25047.0 32 26.0 \n", + "... ... ... ... ... ... ... ... \n", + "1901 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1902 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1903 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1904 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1905 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1906 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1907 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1908 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1909 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1910 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1911 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1912 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1913 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1914 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1915 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1916 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1917 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1918 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1919 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1920 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1921 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1922 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1923 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1924 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1925 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1926 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1927 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1928 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1929 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1930 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 82.0 FR France \n", + "1 52.0 FR France \n", + "2 45.0 FR France \n", + "3 47.0 FR France \n", + "4 40.0 FR France \n", + "5 29.0 FR France \n", + "6 26.0 FR France \n", + "7 20.0 FR France \n", + "8 24.0 FR France \n", + "9 25.0 FR France \n", + "10 21.0 FR France \n", + "11 31.0 FR France \n", + "12 37.0 FR France \n", + "13 28.0 FR France \n", + "14 27.0 FR France \n", + "15 18.0 FR France \n", + "16 22.0 FR France \n", + "17 16.0 FR France \n", + "18 16.0 FR France \n", + "19 10.0 FR France \n", + "20 13.0 FR France \n", + "21 16.0 FR France \n", + "22 16.0 FR France \n", + "23 20.0 FR France \n", + "24 18.0 FR France \n", + "25 22.0 FR France \n", + "26 23.0 FR France \n", + "27 31.0 FR France \n", + "28 35.0 FR France \n", + "29 38.0 FR France \n", + "... ... ... ... \n", + "1901 59.0 FR France \n", + "1902 64.0 FR France \n", + "1903 97.0 FR France \n", + "1904 93.0 FR France \n", + "1905 80.0 FR France \n", + "1906 116.0 FR France \n", + "1907 149.0 FR France \n", + "1908 281.0 FR France \n", + "1909 395.0 FR France \n", + "1910 485.0 FR France \n", + "1911 544.0 FR France \n", + "1912 689.0 FR France \n", + "1913 722.0 FR France \n", + "1914 762.0 FR France \n", + "1915 926.0 FR France \n", + "1916 1113.0 FR France \n", + "1917 1236.0 FR France \n", + "1918 832.0 FR France \n", + "1919 459.0 FR France \n", + "1920 207.0 FR France \n", + "1921 190.0 FR France \n", + "1922 198.0 FR France \n", + "1923 224.0 FR France \n", + "1924 266.0 FR France \n", + "1925 219.0 FR France \n", + "1926 176.0 FR France \n", + "1927 163.0 FR France \n", + "1928 195.0 FR France \n", + "1929 308.0 FR France \n", + "1930 213.0 FR France \n", + "\n", + "[1931 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 = pd.read_csv(data_url, encoding = 'iso-8859-1', skiprows=1)\n", "raw_data" ] }, @@ -78,11 +1045,75 @@ }, { "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
169419891930NaNNaN0NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1694 198919 3 0 NaN NaN 0 NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1694 FR France " + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "raw_data[raw_data.isnull().any(axis=1)]" + "raw_data[raw_data.isnull().any(axis=1)] " ] }, { @@ -94,9 +1125,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
020214334348633073.053899.06650.082.0FRFrance
120214232920423995.034413.04436.052.0FRFrance
220214132504320586.029500.03831.045.0FRFrance
320214032628621842.030730.04033.047.0FRFrance
420213932215518014.026296.03428.040.0FRFrance
520213831561412310.018918.02419.029.0FRFrance
620213731367310404.016942.02116.026.0FRFrance
72021363102897505.013073.01612.020.0FRFrance
82021353126099282.015936.01914.024.0FRFrance
92021343130159485.016545.02015.025.0FRFrance
102021333103927042.013742.01611.021.0FRFrance
1120213231558611009.020163.02417.031.0FRFrance
1220213131885513664.024046.02921.037.0FRFrance
132021303139919695.018287.02114.028.0FRFrance
142021293136269618.017634.02115.027.0FRFrance
15202128386365430.011842.0138.018.0FRFrance
162021273106936838.014548.01610.022.0FRFrance
17202126370864109.010063.0116.016.0FRFrance
18202125379425540.010344.0128.016.0FRFrance
19202124348553011.06699.074.010.0FRFrance
20202123367104455.08965.0107.013.0FRFrance
21202122378795495.010263.0128.016.0FRFrance
22202121378275403.010251.0128.016.0FRFrance
232021203102787540.013016.01612.020.0FRFrance
24202119395396860.012218.01410.018.0FRFrance
252021183121359165.015105.01814.022.0FRFrance
262021173120588891.015225.01813.023.0FRFrance
2720211631650512735.020275.02519.031.0FRFrance
2820211531930615398.023214.02923.035.0FRFrance
2920211432107317099.025047.03226.038.0FRFrance
.................................
190119852132609619621.032571.04735.059.0FRFrance
190219852032789620885.034907.05138.064.0FRFrance
190319851934315432821.053487.07859.097.0FRFrance
190419851834055529935.051175.07455.093.0FRFrance
190519851733405324366.043740.06244.080.0FRFrance
190619851635036236451.064273.09166.0116.0FRFrance
190719851536388145538.082224.011683.0149.0FRFrance
19081985143134545114400.0154690.0244207.0281.0FRFrance
19091985133197206176080.0218332.0357319.0395.0FRFrance
19101985123245240223304.0267176.0445405.0485.0FRFrance
19111985113276205252399.0300011.0501458.0544.0FRFrance
19121985103353231326279.0380183.0640591.0689.0FRFrance
19131985093369895341109.0398681.0670618.0722.0FRFrance
19141985083389886359529.0420243.0707652.0762.0FRFrance
19151985073471852432599.0511105.0855784.0926.0FRFrance
19161985063565825518011.0613639.01026939.01113.0FRFrance
19171985053637302592795.0681809.011551074.01236.0FRFrance
19181985043424937390794.0459080.0770708.0832.0FRFrance
19191985033213901174689.0253113.0388317.0459.0FRFrance
192019850239758680949.0114223.0177147.0207.0FRFrance
192119850138548965918.0105060.0155120.0190.0FRFrance
192219845238483060602.0109058.0154110.0198.0FRFrance
1923198451310172680242.0123210.0185146.0224.0FRFrance
19241984503123680101401.0145959.0225184.0266.0FRFrance
1925198449310107381684.0120462.0184149.0219.0FRFrance
192619844837862060634.096606.0143110.0176.0FRFrance
192719844737202954274.089784.013199.0163.0FRFrance
192819844638733067686.0106974.0159123.0195.0FRFrance
19291984453135223101414.0169032.0246184.0308.0FRFrance
193019844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1930 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202143 3 43486 33073.0 53899.0 66 50.0 \n", + "1 202142 3 29204 23995.0 34413.0 44 36.0 \n", + "2 202141 3 25043 20586.0 29500.0 38 31.0 \n", + "3 202140 3 26286 21842.0 30730.0 40 33.0 \n", + "4 202139 3 22155 18014.0 26296.0 34 28.0 \n", + "5 202138 3 15614 12310.0 18918.0 24 19.0 \n", + "6 202137 3 13673 10404.0 16942.0 21 16.0 \n", + "7 202136 3 10289 7505.0 13073.0 16 12.0 \n", + "8 202135 3 12609 9282.0 15936.0 19 14.0 \n", + "9 202134 3 13015 9485.0 16545.0 20 15.0 \n", + "10 202133 3 10392 7042.0 13742.0 16 11.0 \n", + "11 202132 3 15586 11009.0 20163.0 24 17.0 \n", + "12 202131 3 18855 13664.0 24046.0 29 21.0 \n", + "13 202130 3 13991 9695.0 18287.0 21 14.0 \n", + "14 202129 3 13626 9618.0 17634.0 21 15.0 \n", + "15 202128 3 8636 5430.0 11842.0 13 8.0 \n", + "16 202127 3 10693 6838.0 14548.0 16 10.0 \n", + "17 202126 3 7086 4109.0 10063.0 11 6.0 \n", + "18 202125 3 7942 5540.0 10344.0 12 8.0 \n", + "19 202124 3 4855 3011.0 6699.0 7 4.0 \n", + "20 202123 3 6710 4455.0 8965.0 10 7.0 \n", + "21 202122 3 7879 5495.0 10263.0 12 8.0 \n", + "22 202121 3 7827 5403.0 10251.0 12 8.0 \n", + "23 202120 3 10278 7540.0 13016.0 16 12.0 \n", + "24 202119 3 9539 6860.0 12218.0 14 10.0 \n", + "25 202118 3 12135 9165.0 15105.0 18 14.0 \n", + "26 202117 3 12058 8891.0 15225.0 18 13.0 \n", + "27 202116 3 16505 12735.0 20275.0 25 19.0 \n", + "28 202115 3 19306 15398.0 23214.0 29 23.0 \n", + "29 202114 3 21073 17099.0 25047.0 32 26.0 \n", + "... ... ... ... ... ... ... ... \n", + "1901 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1902 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1903 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1904 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1905 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1906 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1907 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1908 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1909 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1910 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1911 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1912 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1913 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1914 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1915 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1916 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1917 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1918 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1919 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1920 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1921 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1922 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1923 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1924 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1925 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1926 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1927 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1928 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1929 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1930 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 82.0 FR France \n", + "1 52.0 FR France \n", + "2 45.0 FR France \n", + "3 47.0 FR France \n", + "4 40.0 FR France \n", + "5 29.0 FR France \n", + "6 26.0 FR France \n", + "7 20.0 FR France \n", + "8 24.0 FR France \n", + "9 25.0 FR France \n", + "10 21.0 FR France \n", + "11 31.0 FR France \n", + "12 37.0 FR France \n", + "13 28.0 FR France \n", + "14 27.0 FR France \n", + "15 18.0 FR France \n", + "16 22.0 FR France \n", + "17 16.0 FR France \n", + "18 16.0 FR France \n", + "19 10.0 FR France \n", + "20 13.0 FR France \n", + "21 16.0 FR France \n", + "22 16.0 FR France \n", + "23 20.0 FR France \n", + "24 18.0 FR France \n", + "25 22.0 FR France \n", + "26 23.0 FR France \n", + "27 31.0 FR France \n", + "28 35.0 FR France \n", + "29 38.0 FR France \n", + "... ... ... ... \n", + "1901 59.0 FR France \n", + "1902 64.0 FR France \n", + "1903 97.0 FR France \n", + "1904 93.0 FR France \n", + "1905 80.0 FR France \n", + "1906 116.0 FR France \n", + "1907 149.0 FR France \n", + "1908 281.0 FR France \n", + "1909 395.0 FR France \n", + "1910 485.0 FR France \n", + "1911 544.0 FR France \n", + "1912 689.0 FR France \n", + "1913 722.0 FR France \n", + "1914 762.0 FR France \n", + "1915 926.0 FR France \n", + "1916 1113.0 FR France \n", + "1917 1236.0 FR France \n", + "1918 832.0 FR France \n", + "1919 459.0 FR France \n", + "1920 207.0 FR France \n", + "1921 190.0 FR France \n", + "1922 198.0 FR France \n", + "1923 224.0 FR France \n", + "1924 266.0 FR France \n", + "1925 219.0 FR France \n", + "1926 176.0 FR France \n", + "1927 163.0 FR France \n", + "1928 195.0 FR France \n", + "1929 308.0 FR France \n", + "1930 213.0 FR France \n", + "\n", + "[1930 rows x 10 columns]" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -122,11 +2120,11 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 7, "metadata": {}, "outputs": [], "source": [ - "def convert_week(year_and_week_int):\n", + " def convert_week(year_and_week_int):\n", " year_and_week_str = str(year_and_week_int)\n", " year = int(year_and_week_str[:4])\n", " week = int(year_and_week_str[4:])\n", @@ -152,13 +2150,11 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 8, + "metadata": {}, "outputs": [], "source": [ - "sorted_data = data.set_index('period').sort_index()" + " sorted_data = data.set_index('period').sort_index()" ] }, { @@ -179,15 +2175,23 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "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", " delta = p2.to_timestamp() - p1.end_time\n", " if delta > pd.Timedelta('1s'):\n", - " print(p1, p2)" + " print(p1, p2) " ] }, { @@ -199,9 +2203,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": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztvXn8HEWd//96z/G580k+ST4JIQdJuMMNMYAiyCHgsUZXUNSVrOKiLq7uV3/rgq6yK6LieuCxHij3siICCsplQJArHAkEyQUJ5D4/ySfHJ597Zur3R1f1VPdUdffM9Hxm5pP38/H4PKY/fVRVd1fXu95HVZEQAgzDMAxTLolqF4BhGIYZHbBAYRiGYWKBBQrDMAwTCyxQGIZhmFhggcIwDMPEAgsUhmEYJhZYoDAMwzCxwAKFYRiGiQUWKAzDMEwspKpdgJFk4sSJYubMmdUuBsMwTF2xZMmSnUKIzrDzDiiBMnPmTCxevLjaxWAYhqkriGh9lPPY5MUwDMPEAgsUhmEYJhZYoDAMwzCxwAKFYRiGiQUWKAzDMEwshAoUIrqJiHYQ0TJt33giWkhEq+Vvh3bsKiJaQ0SvEdEF2v5TiOhVeezHRERyfyMR/Vbuf56IZmrXLJB5rCaiBdr+WfLc1fLahvIfBcMwDFMOUTSUWwBc6Nt3JYDHhBCHA3hM/g8imgPgEgDHyGt+RkRJec3PAVwO4HD5p9K8DMBuIcRhAH4I4DqZ1ngAVwM4FcA8AFdrgus6AD+U+e+WaTAMwzBVJFSgCCGeBNDt2z0fwK1y+1YA79f23ymEGBRCrAWwBsA8IpoCoF0IsUg4aw7f5rtGpXU3gHOl9nIBgIVCiG4hxG4ACwFcKI+dI8/15z+qWbK+Gyu27Kt2MRiGYYyUOrBxshBiKwAIIbYS0SS5fyqA57TzNsl9w3Lbv19ds1GmlSGivQAm6Pt910wAsEcIkTGkNar54M8XAQDWfec9VS4JwzBMIXE75cmwTwTsL+WaoLQKC0R0OREtJqLFXV1dttMYhmGYMilVoGyXZizI3x1y/yYA07XzpgHYIvdPM+z3XENEKQBj4ZjYbGntBDBOnutPqwAhxA1CiLlCiLmdnaFT0TAMwzAlUqpAuR+AirpaAOA+bf8lMnJrFhzn+wvSPNZDRKdJH8ilvmtUWhcB+Iv0szwC4Hwi6pDO+PMBPCKPPS7P9efPMAzDVIlQHwoR/QbAOwBMJKJNcCKvvgPgLiK6DMAGABcDgBBiORHdBWAFgAyAK4QQWZnUZ+FEjDUDeEj+AcCNAG4nojVwNJNLZFrdRHQNgBfled8QQqjggH8HcCcRfRPAyzINhmEYpoqQ0+E/MJg7d66o59mGZ175AAB2yjMMM7IQ0RIhxNyw83ikPMMwDBMLLFAYhmGYWGCBwjAMw8QCCxSGYRgmFligMAzDMLHAAoVhGIaJBRYoDMMwTCywQGEYhmFigQUKwzAMEwssUBiGYZhYYIHCMAzDxAILFIZhGCYWWKAwDMMwscAChWEYhokFFihMrJzz/Sfw8Rufr3YxGIapAqELbDFMMbzZ1Ys3u3qrXQyGYaoAaygMwzBMLLBAYRiGYWKBBQrDMAwTCyxQGIZhmFhggVInCCGqXQSGYZhAWKAwDMMwscACpU5gBYVhmFqHBQrDMAwTCyxQ6gRWUBiGqXVYoNQJ7JRnGKbWYYHCMAzDxAILlDqB9ROGYWodFigMwzBMLLBAqROiuFAefHUrTvivP2Mwk618gRiGYXyUJVCI6P8R0XIiWkZEvyGiJiIaT0QLiWi1/O3Qzr+KiNYQ0WtEdIG2/xQielUe+zERkdzfSES/lfufJ6KZ2jULZB6riWhBOfcxWrjmTyuwt38Yu/YPVbsoDMMcgJQsUIhoKoDPA5grhDgWQBLAJQCuBPCYEOJwAI/J/0FEc+TxYwBcCOBnRJSUyf0cwOUADpd/F8r9lwHYLYQ4DMAPAVwn0xoP4GoApwKYB+BqXXCNRkQRXhT2tzAMUw3KNXmlADQTUQpAC4AtAOYDuFUevxXA++X2fAB3CiEGhRBrAawBMI+IpgBoF0IsEk5s7G2+a1RadwM4V2ovFwBYKIToFkLsBrAQeSE0Koli8qLKF4NhGMZKyQJFCLEZwPcAbACwFcBeIcSfAUwWQmyV52wFMEleMhXARi2JTXLfVLnt3++5RgiRAbAXwISAtBjwmBWGYapDOSavDjgaxCwABwNoJaJ/CLrEsE8E7C/1Gn85LyeixUS0uKurK6B49Y90PTEMw1SFckxe5wFYK4ToEkIMA7gXwFsBbJdmLMjfHfL8TQCma9dPg2Mi2yS3/fs910iz2lgA3QFpFSCEuEEIMVcIMbezs7PEW2UYhmHCKEegbABwGhG1SL/GuQBWArgfgIq6WgDgPrl9P4BLZOTWLDjO9xekWayHiE6T6Vzqu0aldRGAv0g/yyMAzieiDqkpnS/3jVrYisUwTK2TKvVCIcTzRHQ3gJcAZAC8DOAGAG0A7iKiy+AInYvl+cuJ6C4AK+T5Vwgh1ICJzwK4BUAzgIfkHwDcCOB2IloDRzO5RKbVTUTXAHhRnvcNIUR3qffCMAzDlE/JAgUAhBBXwwnf1RmEo62Yzr8WwLWG/YsBHGvYPwApkAzHbgJwU5FFrluKChtmbYZhmCrAI+XrBBYSDMPUOixQRhEc5MUwTDVhgVInsILCMEytwwKFYRiGiQUWKHVCMaPf69nfMpzN4RM3v4BXNu6pdlEYhikSFiijiNHgQ3mzqxePv9aF/+93r1S7KAzDFAkLlDqhjpWOkhgNwpFhDjRYoNQJxZixihmzwjAMExcsUEYh9exDYRimfmGBUi9EWg+Fop5as7B2xTD1CwuUUYTyO+RGgYpCvFwYw9QdLFDqBJ7Li2GYWocFyqikfiUKC0OGqV9YoNQJRUV5jYJGmcOGGab+YIFSJxQjI+pZnowGYcgwByosUEYRqlM/GpzyDMPUHyxQ6oQoc3mRtBOxPGEYphqwQBmFsEBhGKYasECpE4rzodSvRKnnsjPMgQ4LlFFIPWsoquzEYV4MU3ewQKkToggJ1QTXs0BhGKZ+YYFSJxQ1Un4UmI1YP2GY+oMFyiiENRSGYaoBC5R6oaj1UBiGYUYeFiijCWknKmb9+VqFffIMU3+wQKkTooiI/Ej5Spak+nT3DqF/KFvtYjAM44MFShG8snEPHlu5vdrFiMDolignX7MQ7/3JU9UuBsMwPlLVLkA9Mf9/ngEArPvOe0Y87wNltuH8OJTg897o6q18YRiGKQrWUOqE4sKG6xd1n7xiI8PUHyxQRiH1rKEwDFO/sECpEyKNlJd2Ip6+nmGYalCWQCGicUR0NxGtIqKVRHQ6EY0nooVEtFr+dmjnX0VEa4joNSK6QNt/ChG9Ko/9mGTLSESNRPRbuf95IpqpXbNA5rGaiBaUcx+jDZYnDMNUg3I1lB8BeFgIcRSAEwCsBHAlgMeEEIcDeEz+DyKaA+ASAMcAuBDAz4goKdP5OYDLARwu/y6U+y8DsFsIcRiAHwK4TqY1HsDVAE4FMA/A1brgGo0UEzZcz1OvsDBkmPqlZIFCRO0AzgRwIwAIIYaEEHsAzAdwqzztVgDvl9vzAdwphBgUQqwFsAbAPCKaAqBdCLFIOCPybvNdo9K6G8C5Unu5AMBCIUS3EGI3gIXICyGmjhtlVXQe2Mgw9Uc5GspsAF0Abiail4no10TUCmCyEGIrAMjfSfL8qQA2atdvkvumym3/fs81QogMgL0AJgSkNWopZvR7HcuTUTHKn2EOVMoRKCkAJwP4uRDiJAC9kOYtC6Y+pwjYX+o13kyJLieixUS0uKurK6B4tU00p7zzW89O+fotOcMw5QiUTQA2CSGel//fDUfAbJdmLMjfHdr507XrpwHYIvdPM+z3XENEKQBjAXQHpFWAEOIGIcRcIcTczs7OEm6z/qhjeeJqKGzxYpj6o2SBIoTYBmAjER0pd50LYAWA+wGoqKsFAO6T2/cDuERGbs2C43x/QZrFeojoNOkfudR3jUrrIgB/kX6WRwCcT0Qd0hl/vtzHoL57+fUsDBnmQKfcqVf+BcAdRNQA4E0An4AjpO4iossAbABwMQAIIZYT0V1whE4GwBVCCDXD32cB3AKgGcBD8g9wHP63E9EaOJrJJTKtbiK6BsCL8rxvCCG6y7yXUUM9+yHckrNXnmHqjrIEihBiKYC5hkPnWs6/FsC1hv2LARxr2D8AKZAMx24CcFMx5a1noi0B7DTC9StOWENhmHqGR8qXQK1qADQK1kOp57JXg1OuWYh/+c3L1S4GwwBggVIS2SosOFLU5JB13CbXcdGrwq7eIfzxFWM8CsOMOCxQSiBbhRZbzzKsF1/XAkVNX289Xsc3VwQbu/vwt017ql0MhikKXg+lBHK5apcgmHpucsM0sdG+GqXi7d99HEB11t5hmFJhDaUEqqKh6Nsh2ddzLz5sga16HrQ5GhkYzuLx13aEn8gcELBAKYFq+FCKocaLF0iYvBgpgXLrs+tw39LNI5JXPfOf9y/HJ25+ESu27Kt2UZgagE1eJZCrhlNea0jDc69fiRJm8hopBeXq+5cDAOafOKqniCubVdt6AAD9w9mQM5kDAdZQSqAaJi+dA9kpzyav2mIo4zgUG1P13ZR8+JeL8JZrH612Meqe+q4FVaIqGkqEc9SKjfXc5IaVPU55smVPP4azNR5hUeMMyefXUOcC5fm13ejqGax2Meqe+q4FVaLqYcNFnFtvhGlfcWkou3uH8Nbv/AXX/GlFLOlVilr31ykNJcFT5TBggVIS1f7IbW2q+qTr2SwUKixjyqe7bwgA8NTqnTGlWBkyNR6jnv8W6rfOMfHBAqUEqvONHxgLbKnCk6XHG5eszGSdhFKJ2u5Z10vfoF7KyVQWFiglUHWnfGgkVP1+3Uq7sjbzMd2a8p2kk5X7BJZu3IMdPQNlpVFtbTiM/KJu1S0HUxuwQCmBqszldYB8sKGDNmOSKK5AqaAz+f3/8wzedf1TZaVR6c5Ld+8Q/rx8W9npxPVemPqGBUoJVNtHET5SfmTKUQnCih6XLB+WJq90hU1eu3qHyrpeVNi8+slbXsTlty/B3r7hkq53NZTadvUwIwQLlBKwaSi5nMD3HnkNO/fHH34YLWxYlqOOJUr4GJt4NZRUsrZ9KJXWUDZ09wEAhsuUCNXUUNbv6sWxVz+CdTt7q1YGxoEFSgnYBMqzb+zCTx9fg6/c+2rseUZaYMtdDyX27EcMVXRbFGpct9Y/5Izsbmmo7ckiKm1eTZRZZ9xF3apY5/7w8hbsH8zgnpc2lZ1WPfsfawEWKCVg+8hVb7LS01CE+xnql5Ey56l31NyQDDyvGoNYPflXvIFTg2FLy6cWOjFNaacZG8yUb3fj4ILyYIFSAjYzRLm9vSCKW2CrOl9FPPmGmLxiEpeuQEmHCJQq91hHLP+SNRQHUzmFECMyaaSa9mUgho5ctd93vcMCpQRsvVal/le6Utoa1VowPyhKFS7q0ZIlcDjucSjpEB9KtXuslTZ5xTXA3VTKO57fgHf/+Ck8s6ayg0ebZKdgX39pgQU6LFDKgwVKCdg+8opqKMK8bTy3SkavYsoYmobNhxLTreWfUZhAqbKGUuHoKXX3pd6lGoBqek5Pre4CAOwpMYIsKpPaGwEA62WAQTmwPCkPFiglYI28qXKUVS0NMiu1CKGDNkdYWFa7gal0XYorMtB0uQrNrvRMxDZtthSq3YGod1iglICt15io4Gy/9TA5pJ5tqSavkXLKR02n2gP2Kh02XK6Z1NVwDAmUq/1ERQmBOB5VLXTG6hkWKCVg+8iDPq5yidKwVXtyyOIWAbOkUebxyPmoKV5COrfVbmBGKsqsbA3FsC8fAVbZe4jzEbGGUh4sUErA9pGr3ZVuA6wfaA2th1K6DyX4wrga2GgelNIbmLga0UrXpbJHurvXG0VKiYkWR5wCq9IzE4x2WKCUgHWkvKzYlZjNI4rJq5IaUhQ8Jq8yxVqlm6KogqnUBiYuQRAk0OJ4z+VqtVHMWpU3ecWXT7VNnPUOC5QSsJm8MrJmx+kkLIWq+VDijPIq8XhU3PDkUJNXdTWUoGTiEFpBUVrlXj9SJq84RVa1TZz1DguUIkgEqvfa/grLkxCLV13bgXMhvo24epD5afIrEzY8Eg1TnA112eU1XD9S3ar8Gl/lP496/nZqARYoRaCiuGwaStbVUCrMCPXii0Vv7CulocTVUKt8KuWUj6thChKg8WgoKq3yTF6msozUtCxxps8CpTxYoBSBK1BC5vKqxPLakSaHlL/1/FHkneW2kfLxaiih5anhZxmHthbbOBRDWdyQ5LJSDscNG44hrRp+3XVB2QKFiJJE9DIR/Un+P56IFhLRavnboZ17FRGtIaLXiOgCbf8pRPSqPPZjkoZZImokot/K/c8T0UztmgUyj9VEtKDc+4h2r86v7ePLjpAPpVYdh97Agcr4HuK686i9+6prKAHJxNozLzH4IGgw7YhpKOo3hnzquTNWC8ShoXwBwErt/ysBPCaEOBzAY/J/ENEcAJcAOAbAhQB+RkRqZr6fA7gcwOHy70K5/zIAu4UQhwH4IYDrZFrjAVwN4FQA8wBcrQuuSpHXUMzHXYFSCQ0lyjiUMh2scVKyyatC6foJXWrYd16xjMQriCOPcuefC7p+pHx6wtVQ4vChlJ3EAU1ZAoWIpgF4D4Bfa7vnA7hVbt8K4P3a/juFEINCiLUA1gCYR0RTALQLIRYJp2bc5rtGpXU3gHOl9nIBgIVCiG4hxG4AC5EXQhUj1Ck/QiavkJlfakJtL7kIoRfG1fOPlk7pTvnKv4Q48yh5pHxA3PCITZYaq6ZWAx9PHVOuhnI9gC8D0Pvsk4UQWwFA/k6S+6cC2Kidt0numyq3/fs91wghMgD2ApgQkFZFieqUT1RComjYqrzK9tsPrcJgprJrspjwCr0Se/YhAiO+8R3OL4W8q0ppWnEQRx6qk1TuFC9GoTFCGkqcU6/UQmesnilZoBDRewHsEEIsiXqJYZ8I2F/qNd5MiS4nosVEtLirqytSQW2otsfmlM9UsHdTbMqPryrvXkvBE+VVahoh0Vdxm7ziOs9PXCOug8ehxOGUj0eLCLrcZiKOi9ES5fXX17vw+KodVcs/DsrRUN4G4H1EtA7AnQDOIaL/BbBdmrEgf9UT2gRgunb9NABb5P5phv2ea4goBWAsgO6AtAoQQtwghJgrhJjb2dlZ2p1KEongKK+c60OpsIZinUusdtZHL/W7DBtwGN84lHjPK7yu8g1TnFmUqlEGDmyUv/WkoVRToCy46QV84pYXq5Z/HJQsUIQQVwkhpgkhZsJxtv9FCPEPAO4HoKKuFgC4T27fD+ASGbk1C47z/QVpFushotOkf+RS3zUqrYtkHgLAIwDOJ6IO6Yw/X+6rKJHDhiuQdy2Hryo8RSzZVBQS5RXXY3AboTATW6mmu3gIeh7xTr1SXjqmy11hU2G/RJypswulPFIVSPM7AO4iossAbABwMQAIIZYT0V0AVgDIALhCCKEM/Z8FcAuAZgAPyT8AuBHA7US0Bo5mcolMq5uIrgGgxPk3hBDdFbgXD4kQm3CuglFeOtY6r+VbahkeXrYNb5nZgQltjaUlICk9bLi841GJOv9T6ZpWnYQNh5hxAWDznn5MHdccdHng9PVVmyy1ymkdiMQiUIQQTwB4Qm7vAnCu5bxrAVxr2L8YwLGG/QOQAslw7CYAN5Va5lII01Dyc3nFj6fzHxLlVSp7+4fxmf9dghOmjcV9nzuj6OujlLGYNEzE1VBHNZNUal2XOIjFhxKS1iPLt+HTty/Bzf/4Fpx91CTjOYD5fikmh38YcU4OyRpKefBI+SKIPPVKBVSUSCPlY8r2za7ekq6LYz0UuGbDyqp5+aUGKhNVFtvkkCUeKzofS2KvbNwDAFixdZ/xeODARjftOgobZg2lLFigFEEq6Xwig8PmsJWog+XKxWZOKvdbUCa9oRjCckoPGw5Lt6RkDekIa3562Wt5csiRiPIKC5JwR8MHvLkgc1oc5LXN8vNhgVIeLFCKQH08vYMZ43HVDlcmyksYNy1nlIT67odLFCgx+OTDTVAxzzZsyk/fV7pTPi4fSpATpfz0VU21NfoiRGPMD140HCP7sTgZCU2NiQYLlCJQla13yDxoMCsnRKq0U96K9jGUVATXDFRi9jF8jGFL88bulDckqAuRcsOfK0mcsw3b7lP4zrMR6JSvdJRXjFKANZTyYIFSBKquDQzbBIrzW5mwYW3bdk6ZfbU4J52sVEM8Ek55vQxVX2Ar8FgcTvlgk1deQ7FcHySQRmwur/jSYqd8ebBAKQHbiHjlrK9WnSz3wyr7w/QIvQr5UEpK1ZCOGxkUrKFUU1sLI87Gz5aWuo+w6YQCfSg8sPGAgQVKEajeWtiKjZWokx7/RAQfSil+nFiLXaHoqLinXjG9ylh8KLEFDwQdi8Mp7/xaZ38Ic8qr8wxuN3c9lEr7UGJMn8ehlAcLlCJQVS1jWTwi4wqU+Cul1+RVmUpfbrnjmMsrSi5xEGzy0n0opUZ5Vb5hiiML1fGw3WdoXVPX2w+NXJQXT19fdViglIC9Nyc8vyNN+QKh3PzN26WmUcrxqOQHw1XG5BVfHbCnMxK+gzCTV9DAyHpcQZSnry8PFihFoL4Lqw8lZzejlJ93eORRudnGur5GyT4U6QS2NGDxrSnvSpTAPEptYOL29ZjzKF+TCmv0w6LutMJYd1W6kY7Xh1J+GgcyLFCKQH3A1skhXYFSXq3c3TtU0nomcTrVy728chpK9IQ37+nHp25djL6hwnFDymppelciBg1lJGzxetlKzS5sVUW1NyzKy/wc1bHSyhYV9qHUDixQikDVtUoLlJOuWYhP3brYm7dl23ZOKaHL5Zu8yvehhF1XTLrfeWgVHl25HQtXbC84lre7G/LwNNRVdsoH5qELvhI1lDCBokxeibAoL9O1sgNW8SgvexlKTYspDRYoRZB3ylsEihsFVn5eT63eaS+H1eZVpg+lBnp6YdPXFJPsoBwv1JgqrOb5gY32MgClN1IjPbCx7OnnLdeHvY8gk1kUn+Lq7T04+ZqFWLJ+d9SiWvOJg0oFvBwosEApgjANJReThhKUd+A55eZR9sBIbbtMk5d9pHz0hAczjmRvTCWt6ZgbQn27PF9QuQRnX76Gkr8+JIcSlkqO4kN5fm03unuHcO9Lm6znjCRBgvnWZ9fhwuufHLnC1CEsUIoimg8l6NtevK67NP+IxwFrOUfbX8r0L/VgPi6miEOuQDFpKEEmr/J9KHFoqaF5xOBDyV8f7JS3WbyCwo7zMzrb802ETE4Zhfz4r/IrcFA5rr5/OVZt6ykp3S/etRR3PL++1GLVDSxQiqBcH8qWPf246BeLcNU9r1akfOVSvg/FvF1cGsEXFpOumuQylbSbvEKjvKo8DiXoecQxANOUlmm/fXJIeZ7hmCpT0DgU9Woy2dLLH2c/KIpQKkVw3fvSZnz198tKKVJdwQKlCCL7UCwVTn1Yz7xRin8kSvnKbVSiXb91bz8WrytcIDMOM09olFcReWQDetdBg+HiGNg4EnjHy5Rrrgx2yttwnfrGKQdUGkECxWmCwhz3f16+DVv29BuPBWmbUfBopBE0y0oP1KxnWKAUgRu1Yql1Uceh9FlmK3byiFKO0q8tJV0/537/r7joF4uC06rUXF5FJOuaIAPSMb3KXJENjInYNJSgYx4Npbx87D4UWedDosACNZRAgeL8hjXSl9++BB/42TOBZSh9up/8dpT3NlyGNjXaYYFSBKoalTpSPkrPJmw8QBAjJVCsAjEWk5fza3MBFdNQuwLe8NyDNJSoDcwX7nwZf1lVGJLsT6NSeMpWoXfvCgzLCUFzogUdUygfyn1LtwSUzUlg+77BwOMlayjadhTBrHxzTCEsUIqgXB9KNqg75z8nqBw284S2XZJTPs4orxLTCBMYUdMVQmD5ln0yTXs+YWHDQQ3MfUu34JO3LDYeK0dDiTIrQlz5hU1f7zrWQyaPNAkcEXKtfk4QYR0xdTQOf1cUE2ccK5qOVligFIGqbKFTr1jqW5QpKKKYs+znlCkQYuxVV2oJ4KgS5XeL82GoQaHBZmGjZRcS/WSjnEfpVTzsKcXqQwm5T5uVx7U2BTzHIIEQpdxhnaxy59Ar1nRoWtF0xZZ9eHjZtpLyH02kql2AeiLM5BVZQwnAbvKqvA2l2ByyOYGk5vGOxQIT2lBHS3nj7j53O3h6Fbs5zNk2px/W8JQj3KNrYdHLU2qeeYERLHCCnnFQ2SKZgUMUgiB/WBSKFcymc97946cAAOu+856SylDLwR/FwBpKMYSZvFz133x5lNDIKA2D9eMv245eXAJBa8+X7EPx/Zaari7ogrQQU3pRpjQJawjLaeBFRMkcxVSzZkcPXjRE5BWkFbKmvN0kJqzFjGKKitKA25aLyKch84tFQyndilBOGfT6VM8zHrNAKYIwDUVVBJsmEiW80Zb2SHRgis3CL1BEDF6UvBPYfDyKOQoAkqQLFLsWUupI+VBfT8RyGq8t4TxbG3TeD57ExSEReUF55s1WluMBk2xGMUVFaTtDNRSEa0JRrgeifWdB9xMUwRmcZrT0ax0WKEUQ5kNRPSlbDyOKel/OwCr9w7ANRAtOt7jz/RpXFD+Pzod+sQg3Pr3WmwaChW7UUezJpG6KswsNUxLenr85/bCPvpw1VaKaDvX7CuvFhxE29UqohhKgBQbV+yjfRNi9ucEyJTvl9e3yrAi2tiE8zdLrSy3BAqUIwjUUeTz04yu+xxal719TJq8I17+wrhvX/GmFrwzBZdH3BjVGKd3kZSimCGoItfNLNXmVp6FEO19PNsycGj4DQYjJK0RzNkd52Z9x1HIB4YJCRPiugii2MS9XQIZdxxrKAYLbEwoZKW932ofnEbbQUWD5wpOP9fph3316hF7J5oeQ4xF7k2oEtu08tyE0XBvJKR/qKA5PY2A4i+sffb1gbreomp6ebpBwB4oTgJ79hry8ZbALjbDvJUq5gPBnHWXOsCCK96F4z9GFrel+ony7cUbsVRMWKEUQeYGtEJNYEGHrUvi3vee/728fAAAgAElEQVRoB0Zgcshh3wAv73oo5fcWzehmngCBot1/8DgUu7CxHQci9JoNefm59dl1uP7R1bjlmXWBaVnz0NINGxsRZoqxva+wGbSjDWwsXiPXCXfKh+cThHfcT/EmL/3Zm9qGYoVmPU/twgKlCFRds1VwVRFsH28UM7ftnCgNdPkaSnEp+J9DsT4USyECr9cfbeCkg9qEkGYfSoDJK4IPJTzKKzwNVU/29A979nt9KNEa47DpQMI0GNvthIXC5zU9+zMuN8orethwqQJFL0/4+f4yewRKSJBHlDTrWJ6wQCkGXf03NVJRx6EE1ZcoGoo1hTIrYrFCYCjjvaBY04GxDO5v+HMIakBSUcOGQ8ehBPfMbURx9DannXVa+n2RQVGfnV72MIER7mMx7w+r00EaitvQB2Qdj1Pe3jmIQrHmJn+ZdU09a3jOxQ7erOcxKSxQiiGkdxw2XXeUHlTYFBhAuL07jEVv7MIKOS2J5/oi63GBhlJk+KWJvBPYcjyqyUsTKKZeY1A+XmFgTj/cJxHeSKWlXa4w/FpPJyiP/Lbf/Ohn2PJAw5YAzvsFzem67yBACwz0oejPyTq9S4j2FSL0wojyvnX82ej3ZxJ+xQ6WPCBNXkQ0nYgeJ6KVRLSciL4g948nooVEtFr+dmjXXEVEa4joNSK6QNt/ChG9Ko/9mOSqPUTUSES/lfufJ6KZ2jULZB6riWhBqfdRDGGNmaoI5SzAFck/UiYf+dVz7sheTx5FqjgFDWERGkpY1FAUTS3QKU9Rw4btwiYoj6KivALPNF0bUUPRTgv1oZQ4Q26+zgZ3koK0wKD3pPfo7eH4wWXPyHuPZxxK8Z0+PTjFdK/F+lDqWJ6UpaFkAHxJCHE0gNMAXEFEcwBcCeAxIcThAB6T/0MeuwTAMQAuBPAzIlJrs/4cwOUADpd/F8r9lwHYLYQ4DMAPAVwn0xoP4GoApwKYB+BqXXBVCiGABmmbD3K+WXu1JajTbt6WbW/58kdK8MkX75T3NVLF9LKs06GHlEXfHZSHfiR40F3htVEGT4ZPYqk1Ura2Xq12WHCteTuoDGX7UGyNeTa4968ae5NgDjOXAT4NpUThHVZGAPjor57D34dNfw/7cxgYzpsl/fmECcVIvtMIJq/V23tw14sbwxOrIiULFCHEViHES3K7B8BKAFMBzAdwqzztVgDvl9vzAdwphBgUQqwFsAbAPCKaAqBdCLFIOE/yNt81Kq27AZwrtZcLACwUQnQLIXYDWIi8EKoYOSGQkmYKU4MYHjZcfO9HEWUG2nI7NsULFLupJuxWw5YAsK+/EU1oeeziAWatIIe9fp6fUKd8hLEsNqHvDW6IJjTDBIZN4OSd6mbCTF6RorwCiuY1F9m0IPv1QF5DCKq/z76xCy9t2GM8FsUpr0/86H+fupkryBQehH6dre6/60dP4cv3/C00rWoSiw9FmqJOAvA8gMlCiK2AI3QATJKnTQWgi9dNct9Uue3f77lGCJEBsBfAhIC0TGW7nIgWE9Hirq6u0m5QkhN5Z6/J+eaavErsaak8TEQx9ei7qYT564uO8gocKR/S4IaaUMLvMVBDCen5Bs16EI9TPjyNfFn8OwJP167TNZSwsGHz8fy0JcF1NlRDCRAoQZq5/vxN35RzTpg5T5m8SutSFfu+/VXGIwxMlosI5Yoi1NSzqmUfS9kChYjaANwD4F+FEIWeXu1Uwz4RsL/Ua7w7hbhBCDFXCDG3s7MzoHjBqI+3IeU8MrNqKzy/BcdVNEpAi2E3seS37WaxcLU5iGIvKbTbh/fsFaHLKEe4x6APVS9aYGNnHIimnxeevrmc4WmYz/a9x6CrtINhiz7ZfCiqrbY9yrzACBY45qhHWI/l8w9/n6FO+QgmryCihHh7Jhv1vVC9Lge1C1HLEHa+bn6rNcoSKESUhiNM7hBC3Ct3b5dmLMjfHXL/JgDTtcunAdgi908z7PdcQ0QpAGMBdAekVTHU+07JEdimyuv2IMrQUOy+hfBeVLl++2I/SH8jFSVU1j1uFbqFaXmvy29HNnmZzFoyHVNDG8WeXUyUV7FpRH0N3rDhsEbXpqGoPG1ldK4Lq9Omo7kIPWpvY2wuY1hAgTuHXon1P4r2Pzic087xHgubKdjzXYR0pILKoBiVAkX6Mm4EsFII8QPt0P0AFsjtBQDu0/ZfIiO3ZsFxvr8gzWI9RHSaTPNS3zUqrYsA/EX6WR4BcD4RdUhn/PlyX8VQL1n5UIJMJUKYK06UKC9ruGwUDUXbXYrJq9gPMjDKK8wpHxJabRWaEdJwyhKsLQVpKFEEYzFmLNuZ6vn5kxLWf3xl0B7/UCa4kbH7J+wCAchrGWG+pKBlloOqQpi5CAg3GekCp1zN3FbWnsGMdo73JP07MEZ/aucPWjTJqJOeAsBADS9BXM4CW28D8HEArxLRUrnvKwC+A+AuIroMwAYAFwOAEGI5Ed0FYAWcCLErhBDqK/gsgFsANAN4SP4BjsC6nYjWwNFMLpFpdRPRNQBelOd9QwgRvuhDGah37EZ5BfhQAKcSJXyWuWLj0W37KzXXT7mTQ4oiTF5hPd4wp33QOf78g+byMo4byIXfRzHroYT5J/yHow9s1NIKucSmoYSNZlcaSrgZ13BtSNpO+uHvM8wp743A8k67E4Uo39b+AbtACdVQtH37BobR3JAsOCcbUfMG8j6jWqRkgSKEeBr2QJVzLddcC+Baw/7FAI417B+AFEiGYzcBuClqeculUEMpfKn+jyPtqzdRpraO0pZUqj4Vr6F4L4hiOlCE9UZtl0eN8tKPBUVyhc29FDbgz0YU4WqrD17tJkhoBjdkOjaTWNjzVmUMc8qbjr/Z1Rt4rZ4/YH+fYU75fl9Ib7LIoPko1V4vg7+YoT4U7R7f6NqPye1NgeeEdezCAjCqCY+Uj4h6x+mk3YcS1hBFcc7Zw2XN+XjPye8vZRyKXuaegeGAMx38H3oxWlTJJi/9OQQ5eyOavEwNQMYjjCzpF6Oh2BrKrOrd+wRzCdPXh5uFbKYWWcYQTcz4DHMiH3bsO97VM6ilYS9XlIGNYXXJL1CKJYpgzgR826FRXtq+vkGzadLbdgSX1z/lUS3BAiUieQ3FHuWVFQIqGCSoYgVVh3JMXsKy7TknYiN85T2vWs9T+COLPD3rkDpv63Sq9sUe8aOnYc8kTAC7Jq8Rccqbz3Ed3vZguZBZFaJpa0CAhuLWyeDjRp9gwHPSIwCD6lyUqVfCnPID2lxoYeax0qP67OXMhAgUfVefxaEepcOo3KLlLqZWSVigRMTVUKTEMDZEOU2DMbzzKLNf2Bc60tKJ4pS3pB/syM5vb9rTbz1PEWTCC9VQbAIjJHAh6lxeYQJYXRs0hY5zrTn98EWftG1LYz3slsE+QDQwD207TGOyNUJ6IImJoHEoHrNiQN5RTZOlaii6ozs8xLjwOUSZakf/3gtMXiFOeT3N/qFMwXEgWpSX+qbZ5DUKUC850OQlRN5pH2TyCqjztm9P3x0UWnz4pLbALKI2wgmLRNIbLv+EhEVNvRLSwNmd8uayBJ1nelxKuzKVI9JAt9BFnyIIJSVQSlxK2fO8Q01ewZ0QW0cmHwpfeCzIxKs3suU65cPaz6EiBIq5wc9vh4VPm/IIMof599nWnC/GXMwmr1GAqjNpy8BGIYTjiE/Z5/qK4pS3OZOj2HmFyKvFNoLKoNfjhCUh/fqCFRsjmA4UtkYifP2NaI1o2AeqTDJGH0pAb1QRNhbA+yyCe71B0XJRe/7hTvngNXxslwf5tLwmL3t+QUWLIlDCTDyD2Zw74DhM0Jt8SUEO9/w55u8SCNeyPD4Um0CJoBWrb5JNXqMAVYmUyctf+dW/DQEaTLFhw548tM2gCSRVpbP2tALsbnrezf4QNUPewdOuh2koweYNq8lL2x+obYV8oHkNJayhNOfxqdsWW/N28gxPQzUM/vuIOrVGURpKic87aKR81iN4vcf1XnRQvdcbx1JMXkIIDGVyaErZvzsdky8pUlSfR/B5j3l9KKYw9Pz271/ebEzfo3nbTF6yj8cmrzpnY3cf/rLKGfDvTg5paQTSKfNx2z4/tsiduEbKB/Vu9Os7WhuM5+jravh7e8UMzgr7cKOMQ4lq8jJqKBm7hpKL0MCEoV8V1uv1NxBhIc/588zbOmGNUPg4FPv70OunR6vLCXd5hLHN6RAfSmFefoKc8kpAqLEd4YLVpKGE19vgKC99csjCa/Xz1+zYb0y/GG2zlk1e5QxsPGC48Pon0StV1ZRl+nq/jyU4yivoA9Mrlp6+fo69rGSZEl0R1YcyY3yzuXzax104fb05rWLKoZK0R7tp50a8F3+jLIRwTV5hpsmIykIB0Ry9Uqj5nmPYNCpuuhEET5IIGSGsaYbdX9C0Jvqz06O69E5HUzoR2NHxNsbFaygq3yapUYeF65uEUxQBHqTF6M/WpKHoQm7K2MIxKP58be+EQAAEm7zqnV7N7mkzealGKGi9lDCHs5OO9oFZbNT2Xlg+bNkmUYIFinaedeZXXaAETQ4Z3FKVPg5F/3ijCRT/aaoRSicJ2ZwoaETimJUgLCgA0CPNAjSUgDzCwlUBIOFGJZamoeQnjwxuiHXHuF53GlPJYA0lQgchyH+h8lUmWlNbqws4k6aWCTDd6WWwrW4ZZA7Tyzx1XHOk6WWsz4tNXqOPlCWKS1WCxpQ9ykt3gFrt6llzxfKYvMpwykf1oYSNW3DOCZjLK6QhDhtvY/umoglWJ52kRfj7GyGb+dKfXzFEeRYqH/+zjtpgeBohSx6qA7R174A5jdAFtAJ8TTaBou1vTCVCorxyoeMrcgGCU+XbGKCh2L4p0z5bvcvmctYhAbqwNmoocldbYwr91iiv/LatbciHDZdYKUcAFihFotZD8TfqamR5Y8qp2GGqtT1E0iw4oph6BDSnvKVvG+xDCVbd/dcXhLtq26GLIoVGHYWbP8Lm8koSIUEBAkXa3YOc4lFmNwgrp+1ZqOcXqKEEZJ+LUE5lAt2+zyxQVO/dlE8uJ9x6Fxa9pJu89Aa2MR0mUITbCYviEPc3poNyUszmdEB0ZYCZFvBqMEERebaAm7CpV1SZWhuTnlH9OlGivJTgHZWzDR+o2JzyX7rrFQDAlr3OgEBTg+kRKJaKO2w7J0JvNJvL98rDzCwm9EP+kGA9D4V/PZRinNnW9TlCzIL63jCTF5GjUfrvWZW7pcFxIdoESjJBxo87yiSa+jnWadlVlFfW30Dlz+8dNA+E85czbLLN/ZZ08svnGq7V0jStt6I/N/24/l6aUsnQsGHVUEcxs/rrnMpX+VBMj2FQm4nZOAefZ7ZieznTSfO3FaYBqbrQ2phCJieMbUMxEXtfv2954PFqwgKlSGzroaza1gMgvxCPUaBE6F3bBlB5HN4BzstkiM0rKGJGz89mc9d7eEEaSlibG6ah2K6PrKHkBBJESCWo4F5UI9SiIoN896HeUypBRsGoPwPb4w7rtQL2cSj6c/3S714xZwA9EITs43rkOfsGCgWKEEIrW3QNxJ+//7jH5JVOhER5CTSkzKZHUz7+ZzXoN18a3pc+kt6koYQNTFTnpC3m7uGQwAJ1fluj04ExaSlRQpebtFD+WtVSWKAUSdrVULz71Qse0+RUmrB49yghknpnKsqAvlxOuE7YoJ6WjShO+bh8KJ6Bbx6tLLicHh9KmMkrIQWK77zhrFegFJicZCEakgmjNqI3nspc48fznCzrV9imf4kyAFbPI51M2DsZcv9gWCNmKKJNAzFdH+SUj2rysmrtWXs5okR52cxx+fTDTV6O4LOYvHQNxRRlJne1KoFi8KNE8aHog417DB2EWoAFSpHYnPKDbq9XmlFCKm6UiBbvOJQ8Qau+hWkowxF9KDaTl974FoyUj7g8L+DTdAy+gChRXrv77DMi6yYv/7N2e7XyXRU45WUeyaTZ5KU3zrbn7REoIc/S38BEFigiL/hsU/2otEwahid6yqCh6L3gIIHSkEz4nPL57c62xkA/VFbkG+oonSx/J8Zv8jLlpa+2aHq2usYQNFo/70PxH9PqvXGOP2nykh0Ys0DRO5vGImA4m3PTGAxZUK1asEApkvamNAB7o656W6ZGZNBiZ9YxNa5ANKe86pUD9nDTKKG2DcmE1Skf1PP2aijWbJxrLcI13OTl/LY2JLGhu8+avhDIm7xy5kaoJW12yudyjjBKktnkpRrnxpS5IQe8AtWqocjGZ3ffkNfnEjHKS9UPm4aiCxG9UXXLFTI9it7wmQSSeldN6YTR5PWTj5yEsS3pwM6F7uy2+1DsYb95gWJu7P1lN5laleBsaUhiwPCcAPgmfvVrKDk3WCfYKW83eUUJVc9khZuGScDXAixQIqB3QtuaCnu1//Dr593t9mZH4JgaEZvjUmfYco5nAJitsfU45cN7e/5z1LHGdMIamjho6Yk66enbwRLF0+s0+I3C1oWZ1N6EfQFrtmRzzpgcx4fiTavAh2IwOaUSBCIyCmbVOI9rSdun4fdoYMHCeTgrPE7zqBqKOi+dMjvl9cbTqKGEOKNVw9fakAx0yrc0pIwmr3QyAaLgzkU2J9CYDjZ56WX3jxIPCwHXz9HLpqPm1+poabD6JjK5HJpkfSnweeUEUkknojBoyiUlDEzzeUWJLMzkcm6dNb3PWoAFSgT0OnLS9HEAvJX/6TU73e3L3z4bgLkR8QgUy8ezp2/IeM5wNrzCZbWxFzaC5k5S/zenk9Zesoo6GtOYwlCBUz68jIohiw8ldKCd3D+2OY19/XaBksnlkEomjFFeqvG2hg0Lx6GfILNg1KPEImkoFuGsa7Hdvfn3Hrb+Rz7dHJIJQjphdnyrfInMPpQBzWxiuk/V8I1tThsFinoXzT6Bo8qSSpCj5YWZvAIGAwMhJq8QfxgQHuXVr92nbc33bE64ocmmcUPpRAKpRGFdA/LPSTnlTUIrinafyQnXpG7SOGsBFihFMP/EgzFeznFl+kgOmdCSV0l9lW5gOIs/r9ju/m9rNJZt2edue6aKzwYLIyGEa+YBok294m8k1LGmdNLaS1aNTHtz2jCXV347rJNtW0NCN3mZGjkRVaBkHS3D5JT/d7l4mNI8/eY9pekliIwaiPqYm9NJ40h7/T6AoIi2nDuzgUegRJxaY1iaixKWaDT1ftsaU8Ye7TV/WuFum16X29C2NBgbWlWHm9NegaI0TqfXbi6bQnd226Phwk1eyjJgMll5ymb47gaGs2hMJaTJy6ahCFcLMkXlpZLkhJkbBYrzG+SUD4vyymRzECIvlFhDGQXMmdKeH31teOnrd/W5UWD+xvaBv231/G+qNPsHM+4klP48wmLd1a5kyBvV0/ELlKzmG7A1gkpDGdeSLhCKUUIfFbbwY/3WTO2LgCMI2pvT2BuooeQ/cv+72LnfWZ5W9fb8jUzGFSi2qCGnQRjX4jRiYWGgVt9AVmDSGGduJ30ke1QNZSiTQzrpaAFmDcW57zGNKQxmcgWC76nVec3adJ/9w/JdR9BQBg0+sVRCCbvgObLcgY02M7DeCbIIlLHN9ndhCxhQ9A1l0dKQRFPaLlCGszk0ppJIkNnklUwknLoW4ENpa0xay5ixDBdQDKj7lHWOfSh1zFlHdAIAPvX22a4GYKv8KgrMX+l65Upt0+Wki6ae0kafk9nvZ1DCzLgcq/YRA3antidyx/BhAI6GYjM/KA1lXEsafcPe0MUo8fymvD3zl4WkMZwVSCcSGNucCgyddM0QBpPXP751JgBg7iEdBXkCeQ2FyBzl9crGvQDgaqumQYNhYzgA53kdcdAYAPlxTP5rD+1sNV6r0m1IOX4KUxYq37amFIQo1AA+euoMAMCE1gZjfdHf9VC2UCDpJtIhTWCpuq/8CkDQlCbhGoruV/Q3pEqQKQ3F1PsPG4fSP5xFczqJpnTC6pQfzjrlTCcThd9NVgr2hFmw+30oYWU0CnfNLOecz1Fedcu3/v44/PXf3uHYq6UG4jdpAcCJ08e5x/0VV32c//yOwwAA1z28quD6/37kNQDAl955BABv4z+cEWhKJdCQTGC/YRlRVQkTrg/F/HEGhYLqDYTN7q8E4+QxTdi8u98jOL2r9Bkv1841a1xhY3Uy2RxSSUJ7k6OhBPV8U0nH5FUw9Uo2hwmtDe6sByYfSpIIiYS5Z/0NaSqaIAVK72CYhmK3y49rTqO1IekZEa9MRqfOGh/4HIczzvxSyQSZfT3y/aqGzG+2ShChIZnA2OZ0YCOmNDH/9apj0+w6q4V7X4DjQ1HCImgga35go/mcvoA669dQTBqGvs+k/fUPZdHckERjOunxK/nzSScTBSHSgFcbDpr6JSjKS/eJmB6VugcVZcoaSh0zdVwzDpng9BTdSA1Dr/Sb7z8W6YQKgfS+8L7BDIjycyot1PwpALBy6z7X3HXstLEACiOq0qkExrWksae30NSjh/w650cQKJZ1OBrTCasdv3cwg2SCcNSUMcgJb8SKbpooZmCjd1R5fr8pCTViub05jUxOWOdGGs46ZohUkgoas4HhLJrSSTcM1LQeSVI6lIMirsa1SA3FoCkph7kqi7mMTrhpc0PK8xz1MFOTM12/XgmUoCivNkuoaf9QFk1pR8Mx+lBk3p3SLOfXCPUOCKCtginvN5VMuHPbBTm78055832qcgJ5049C3dO4AJOXN4LO4JQfdgRKUyppdXYPZhxtMG0wBytt2GbyUs+lQ9YXo0DRBJnZ/Bgs3GsFFihFkk4m0JhKuJVUffyzJrbi2Klj3SWA/Y3I7r5hjG1O4y0zxwMAZvtMGW929brb42XF0yuZchzO7mzFsi17C8qljwlw8jdXuP4gDSWbN3nZ7Pg7e4YwobXBFay6gMqGCAMd23QV+nP76u9fxayrHvBeJ80Lqqe2r98yR1XOOc+koQwO59CUTriRQf75spzjSWkus3+4qqG+88UNBcd6BjPokB+/aTDpUCaHHfsGkUwQWhqS6NO0zv2anyqo4VC2/0SYD6VJaSjehqx/KIuWhpQTHm14Ya9JM9zBcg0Pv89Kn/QQyNeFjKahKP+IzUQzMJx1y2frdfcNZTChtdGThyKKhqILQpNw7xvKoCWdkiYvuw+lIelYKIZ9ocuqA5KyOOUHXe3CHjasv2fTuyg0ebFAGTWkkwk89+YuAMCrm53Gfe1ORyA0p5NIJqhgjER37xDGtzTgbYdNxAnTx+HgsfkFrF7b1oMr/u8l938Vl6/suUvW78Zdizdh+75BzJky1iN8FI+tdLQb1SP0V3oAWLpxD771YN7UVqi6O1OJNyTtTvlt+wYwub3J7ZXq9mD9Yw2bgr1vUL/OayNXvol7X95snIgvlUigvdn5OG1jUfJRXokC4ag0lLZG5+P0m6wGMo7wTicTgavjqSixO54vFCgLV2x3NRjTu7jhyTcwlM2hbziLKWObPP6zrp5BtDWm0GGJrgKc+vTnFduxaluPnMTS0JBlQjQU2TNPJcjY0Kr7mtzuCJSHXt3qu95pqCf5NBilnSd1gWLo+WeyOfQOZTGuJY2GZMJobsrmBJ57s9t9z/4Gv3cog3SSAs1JPQOZ/BT5hnrZP5RFU4hTfkhqKA1GDUUglUwgYdFo9elh2hpTxgk/dS3K9M79Gor+Ll/asBu3L1pnLDcArNvZi9ufW4/dWiRhpWCBUgL7BzN4ZdNe7O0fxmf/dwkA4GvvnQPA+YgOHteE5Vr4rxACD7y61Z0eorOtwRMmesH1T3rSz5sJnEr0wZ8/6x7raEmjfzjrqVC9gxn862+XOtcqDcXQK37itR2e/wsCBwazaEknkUrabME5vL69B5Pbm9zevW3airC5hlZtyz8fvwN70phGz7n+8OmUR0MxCxRHo3Puxa9l9EuBonrWJg2lMZ1Ag+HaJ1/vcrc/ePI0Y97PrtmJnoEM1uzYj6Z0wo0q07npmXVu+Se0NXie146eQXSOaURrYwq9QxljI6iPfbIJhF8/tRYAXMHpb6jUc2hpSFrX6QDygvP7C1/37FeCWK1CqOq06oG3SL8EYNYcVCPa3pRGUzphLMP6XU7nST0f/znrdvbikAmtsgNARo11/+Cwa27yN/hCCLy5sxfTOpodDcXQmO/YN4D+4SwmtzeZnfJKG05aBPtwPnqyrTFlNJH+bvEmd9v07ajvTAUf6N//pTe+gK/dtxw7esxLFCzduAdf+8My7O5jgVLTvLpprxvueeqs8e7+U2Z0uOYCAPjho6sBACu2Oo1oR0tD4Mtt8mkop8hopHOPmuSO1Nd7NNu0tS5Ub9Q0Ut8vI/wRLXv6hjCupQENyYSxl7Ri6z5s3TuA846e5ApHXaCohq8hmQgM6V2xZR9eXLfb/V81QEI4U3tPavcuk6r3XIc1Hwpg1lB+//ImLF6/G21NKeM4FEdDSbjPyh+lNZjJSWFU2Bv98WOr3e2O1gZ85qxDXR+AokdL7/ip4zzC070P7fmmffl0SYHSOaYRQnjHqCh0IWvr9SqfXKcU0H4toV+azFoaUh6Tm58jZSTawb7la9U1x0mf30pZv9XzbGtMYcb4FgDA69sL11JXDeeYphSaLUJtvdTcrjj7UACFdXZ7zyCmjG0CEeGgsU3YKpeP8Oejevb+97l/MIOegQxmTmhBk1xd0n+OmjNuythmo/auzNE2n9tAJuuYxJIJtDWljFGBal9DMmGs0w/KYQcmk5e6dmePuU1R9VG1HZWEBUoJqIaoR3vxHdJMAwAzJrRiy95+DGVyyOWE2wipnt74VkdDsUUouRqKbKwPam/CIRNacMOlczFG9sz1Xo6+eJKqcFHCkvWBbWt27Me9L29Ge3Ma49sasGu/t3y7e4fwvp8+AwA4fto41+Q14JnvyTl/QluDVaDs6BnAtx5c6dmnZgdwBgmiQEPRG5HhjPKhSJOXoUf6/37rTPm+flcvUoaonL6hLJrTKddMon/gj6/agS17+90er8lcpfQlFbsAABkJSURBVNPakMRQNufJQ5+rYObEFqzdWTjnmOrlThnbJAVKPuT2hbXd2NkziM425zns6CnUcO55yenRHnXQGIxp8oZQCyE8GsGRB7XJPH0+FBku29yQLLDrK63x8+cchiljmzF1XDNOP3Si55zdfcNIJwmHdjrp7+kbwlAmh6vvd9braGtK4ZiD2wHAM75KoepIe3MazWnz4lOfuPlFAMA75xyEhlSi4JydUvgCju9xj2/CUCEEVmzdh4PHNiNBhVFe6tlOGtPkdpL82tSS9U7np6UxKc2gJsGcklFehR2xvf3Dbn1tbUx5OhwKIuAj86YXvEvF75Y471tpgyZ/097+YezoGcAHfvaM51tX7ZTS6isJC5QSuOHSUwAAr2zKO8eVAxYAZoxvgRDA5j396NLMHX/+1zOdc1sd2/gf/7bV2LNUGorqhewbGMa4lgYkE+QKM70Xs2NfPg/VE+s2aEDKfPC2wyYAAFbvyGtRf/eTpwEAa3fuR2dbI4ayOc8aGjc/s9bdntrR7H58i6QvCchHvk0d12wVKN95aJVrrvnvi453yiqj1vo1AarT3Zu/v739TnBDkIaieH37foxrTruNjBACv37qTaza1oOJbQ1OGGgq4b6DwUwWn7jlRbzZ1Zv3ofh7oz47vxp39JdVTtTevS9twuW3O2bQmz/xFsyc2Iqd+wc9nQ8AmC0b4S9feJQnn21S453a0ew2Hs9o5i0A+OJvl7qDEu/41KkFvd4bn16Lo772MADgAydNzXdCfL4iJ3oqiYltjejyCS0VHq6ec2tjsqCubt834DbEHS1pvL59v0cba0zlI+mUANTRNZQmi0BRTGhtQFPK6zTP5YSrzTnppAue8/pdfVi/qw/HTxvraJyeFUdz7rOd1N7odsb8Qukrv3dmVmhtSCGdpII6oYIjGlIJo69ob3/GTXuMQZvsH8pCCGD6+BajQNE7dp1tjUglyO0c6H6t2xatw7xrH8PLG/bgZ0+scfd/92FnOIJtqYU4YYFSAvNkpNYv/vqGu0+NugbgqvkbuvvcGXG//ffH4fDJjulARXF9/jcv46p7ncqqKhwAd0SuMmXtH8y4PRz1q1d6XUM5Qubx8yfewLNvOB+LEAKfunUxXtqwBx87dQbu+NRpmDquGeccNdm9Tn3MbY1p9wPVG5n8+BanrEp9vvHpvKD5yeNOJZ7c3mRdcnbDrnzP6e9OOBgAXPPf/a9sAQAc5DOt6D383dIspyKDFr2xy3Pu9Y/m7fzXf/hETBzTiG37BvDC2m6s3dmLbz7gaEfK8d/WmG+Mded8Vjg+AP/H7ze5KHPkvS9txsPLtuGLd+UXxDp99gQcPcXpob+8YY+7f82OHqzcug//cNoMTGxrRIMW2qzs4JedMcvt3X/7oXwgxcBwFve+vNn9f0Jbo9uQqobnJ3/JNyZfvvBIHDLBqY9vduXNTu/60VNYsXUfWhqSmNbRjF29Qx6z1w/+/Lr7fACnfvf6zGI79uUb85NndOD17T2eNTv86A1j/1AWH/nVcwCcnrNp2pOZV+Yj/KaOa/aYxYYyOcz+yoMYyubcABdTY/yO7z0BADhpRgcmtDZ4Alr+sHSLu/rhpDFNmNDm1AndxPgH7VmPa0ljrG+GBiEEdvcNYUxTCq2GZwTkO0EA0N6cKjB3q47dhNYGjG1Oe+bzA4BNux0z3tfeOwcpae5VZrjP3pEP5nlo2TZ3u7vX0RZ/Kdsoovxy0JWkrgUKEV1IRK8R0RoiunKk8k0lE7j4lLxD9leXzvUcnyk/4K/ft8xtQHUfy8yJ+ZDhZ2WD+LmzD3P3JROEs47odKdr6RnIuB/2eFnp/+HG53HPkk0Yzuawbd8AWhqSWPIf57kCBXCccYBjmnh0pdODfvvhjtlicnsjHl253W1EJkrzyvUfPtE1tej2aHUfi646BwBcM8eH5k537mPNTlcNP27aWKzf1edpwFZt24c5X38Yi9fvxvTxzVhz7btcLedHj63GcDaHr/5+GQDHST7/xIPdxlhpVkIIrNrWg9aGpGsWfGjZNjfC7oW13bhe+qs+d/ZheP9JU937+tAvF3m0qQlyv+5/0IXHk6934aD2JqzesR83P7PWzcPf6M2bNR4nzxiH7fsG8BkZoKFoSiddgXPpTS/grsUbccfz63HeD5wgjCmyIUwnE65PRQnIye1NrvYDALc/t16+k0JBPaYpheGscDXatHbdlLHN6GxrRHtTynMPyt8xub0R0zqccmyWDdf6Xb245dl1APL+l9bGQrPY02t2YqKsjzMmtGBDdx9ulsEGd3zqVPe8f3r7LADe9Wvu1jSW9qY0mn3C+5Hl+cbxPcdPQSJBGN/a6AY4vKHVLfVttDWmsHrHfsy88gH0DAzjVc2CMOfgdpx1RCcWr+t2972wNl8fpoxtcuvELqkR7+gZcINdAGD2xFZMbGv0+CqufWAldu4fwgnTx0ktzvuMvvHHFXjy9S7XfzllbDO27On3CFf1nR49pR3TxrcULMvwyibnuOrIHtbZhpekGU7VLz/PvrELtzy71u2MXP/hE43nxU3dChQiSgL4HwDvAjAHwEeIaM5I5f/fF5/gbiuNRKGcyut39blLuE7tyIcJHz0l3+irD0SNS1HmrtNmT8COnkHMvPIBrNmx3+09Te/I5/Wl372Cf7z5Bdy9ZBMmt+c/iHs+ezoAYOueATz5epenYX/nnIMAOH4QAPj07Uuwb2AYu/uG8LmzD8MZh0/ElHFOWZX2BACL1+/GWUd0uo0g4HyEj63ajgdf3YqnpOng8jNn47yjHc1n4YrtGBjOYmA4i6Ub9rgN0lffPcfTWA5lcp4ghpNmdOBHl5yEBz9/BsY2p/E32TDMuupBAMDqHc79KOF45T1/AwC8qDUWSmirBg+AK7CAvNmvrTGFbtnQ/dNti6GzRj63//rjCpz9vScwnM1hYDiL6eObsfTr73TPe8us8W74uJ/2prQ7WvzLd//NUwY1pU9zQxK9Q1n84eXNeFM2+v5xSl/7wzJ89FfPYYUWPfjdDzomw6nyfT2yfBuEEOgbyqA5ncR/vOdoAE7PtH84i9sWrUcuJzxa3aWnz3QFysPLtmFP3xCeX5t/jm87bKJ7H+t29kIIgU27+/BdOdODKu+sia3oG8q6pi0lyAHg7KMmAQAWrtiGN7scs9ia7fn33TmmEbMntmHF1n2uBvLEa/louh/JxnDquGZs3tOPXE54NIcTZzh1eYzmI3izqxd/99On5f07186c2IrdfcPu93CXjKy6+JRpaG1MuXVl8x5HaF/30Gtueo996Sykkgl0jnGEWi4n0DMwjF9LDf3sIydhfGsjXt2819Uyt+7tx03SVKwiNie0NWBgOOdqut29Q/j6fcvR3pTCcVPHYvbEVmzs7nM7Z3v7hvG5/3sZAHCE9IWdecRErNrWg1ufXYcl63fjrYdOcNug5646F58+czZ6BjKeIQLHTR2LkaDybv/KMQ/AGiHEmwBARHcCmA9gReBVMfKhudNw1+JN6GgtdHZ9+qzZ+OVf33T/Vz1qwKn43/774zwN9jEHj8WzV57jOrs/NHe6x9Sh9rc2pvD7f34rPvAzJ5T4mTVO49CqmdxOOWQ8JrY14Pbn1uP259bj9NlO4/noF89yR29/4dzDccuz6/DU6p04/j//DMDpZQJO4/C+Ew7G/a9s8ZgdPvyW6Z57VD3Cf9bU7q+8+2jXofvth1Z57gEAvnvR8Th/Tt7Udu0HjsVXf78M/yidr7dfNg/zpDZHRDjnqEn4/cubXS0FAP7noycDAP7rfcfgnO//Fc+v7faUEwDOP8bJQ2/YAKAhlcBzV53rmrxaGpJ48vUuz/WpBOFfzjkcJ80Y5z5fADj8qw8BcOYCU2NMALkqoebrPXH6OPzkIye5/5scqLd9ch6OlR/5yTOcXqbqDc/ubHXry9P/fjbOuO5xAE6vU2m0L371PFd7OOZgJ50v3LkUX7jTSeM/3nM0PiWXUgDyQRqzv/Ig5p94sLt/+vgWV6P5/sLXPaHBT335bFeLPOWQDjy0bJsr1BU//JDT2L/98E5335jGFI6Y3Ob+f/KMDrQ2JN1ZnnXOnzMZzQ1JvOu4g3D7c+tx9Ncf9hz/+5Onup2PiW0NeHTldsz+Sr4MS7/+TvddHHVQvqM2/3+ecbcf/PzbAcAVnOd8/6+ePFTnUJnOvvaHZchkc65wfOXr57uTMk4Z14xMTnjKMLm9EVPGNuHdxx2Ee17ahPN/+CTecUQn/rB0i3vOVe92hLsyfR399YfxH+852jWbz5s1HkSEWRNbkRPA1fcvxxmHTcTdSza6aag6MW+W8z2r4Ie3HjoBnz7rUBAc68ln33Eofvlkvu1555zJmDXRPidcnNSthgJgKoCN2v+b5L4R47oPHo+Xv/ZOd2CXzlXvOtrdVj1FnY/Mm4EPzXXMZj/96Ek4aGwTDh7X7EaLdbQ2uA0nAHzhvCPc7ZNmdGDVNRfiY3JyPwC48l1HedJXjRTgOM4bUgnXlq7Sv/szp3uu0c1lX77wyIIyK/OW4tNnHer5/8PyeDJB+MTbZhZcr9LQ/THzT5yK8a0N2Ll/EO1NKZwg15tRvE/6WVQAxA8+dIKrfczubMOC0w/xnD+xrRFrv/1uN6LlpBnj3OcMAB+Q+bn5n1RYZdZ86934wnmH48wjOnHjgrkFx+dI34bivccf7Pn/G/OPwXRNa1XBB4qzj+zEmUfkG+CztG3AG5QwraMFP/zwCZ7Q5PamlCtMABgbC78pRJmdAOC+pVvQ3pTCm996NwCnQTQxWSuH/x4B4N8uONJ9X7MmtuL//ulUzJs5Hs995VyPvb4pncRtl80ruP5nHzsZN0hz8emzJ7h+McXnzj4M130w/+z893TqrPEewX7G4d4oNAC4ccFc13RqMg8pYQM4fkLV0/+vPzr90k+fOdsVJgBw3tGTCibs/Ou/nQ0iwrlSM9/TN+wKk6Z0Ag98/gycfaSjpZ0wLV+/vynNZR89dYZ7n+dIbe43L2zAFf/3Eh5/rQtnHtGJ1755oXvdSTPGuQEbc6a04+Onz0Rarv0DOFMCPfrFM3HqrPF45evn41eXzh0R/wkAUNjKerUKEV0M4AIhxKfk/x8HME8I8S++8y4HcDkAzJgx45T169ePWBm7egbx4Ktb8eG3THd7ejoq5t10LCpPvt6F46aO9YQtA46d/JZn16GlIYk7ntuAb8w/BqdKTUVHCIEbn16Lkw/p8AghwImCcSYedCZMTCcL+x9CCCzbvA9HTRlTcHzfwDC+9cBK7B/M4JK3zMCps8cb0xjMZHHzM+tw3tGTcNikMQXHN3b34Z6XNkEI4F/OOcxjLstkc3ijqxcvbdiNvqEsPvm2mcaPZ0/fEBa9sQvvnDPZc73KP5MV+NumvTj90MJnBABPre7CP922GD//2Ck464hOj1AEnOd9+6L16BkYxhfPLxTGQgg8uXon1u3sxfwTD/Y0hIATeLG2qxdLN+7GOUdPds1YfpZv2YtpHS2eIA51D3v7hnHvy5vxlpkdOOWQ8QXXDgxn8aFfLsLfNu3FI/96pju+BHC0qA3dffjmAyvwRtd+XPfB4/FWX5jwjn0DuPulTZg3czyyOYGTD+kwvk8bb3btRzqZwN1LNuGIyWPwnuOnFJRvV+8Q7npxIz5w0lSPrxFwvpfVO3qQlCtp6h0gxeOrdqClIYk3unpx5EFtBc8hmxPY0N2HSTJYQ/kCdfYNDGPFln14bOV2XHbG7IIgEcDxMz2zZhfOO3qSZ9zU46t24InXdqBzTCP+7oSDMaYp7enAAI5pcenGPTj90AkgOKZbvc4u27wXT67uwtY9Azi0sxUXzZ3u+okUzpit/EzNlYaIlgghCntX/vPqWKCcDuA/hRAXyP+vAgAhxLdt18ydO1csXrzYdphhGIYxEFWg1LPJ60UAhxPRLCJqAHAJgPurXCaGYZgDlrp1ygshMkT0OQCPAEgCuEkIsbzKxWIYhjlgqVuBAgBCiAcBPBh6IsMwDFNx6tnkxTAMw9QQLFAYhmGYWGCBwjAMw8QCCxSGYRgmFligMAzDMLFQtwMbS4GI+gEEhRaPBWCe5c9hBoDCBcSLS6PSx4HwcnIZo51T6TLGkUa5ZYyjDFzGeMpQy+3PkUKIwqkJ/AghDpg/AF0hx28o5/qIaVT0eEz3OerLGMf7HqEyVLSMI1QfuIwjUMZKlgHA4rC8hRAHnMlrT8jxP5Z5fZQ0Kn0cKP8+D4QyRjmn0mWMI41yyxhHGbiM8ZShHtqfQA40k9diEWE+mkpdP1LUQzm5jPHAZYwHLmM8eR9oGsoNVb5+pKiHcnIZ44HLGA9cxhjyPqA0FIZhGKZyHGgaCsMwDFMhDniBQkQ3EdEOIlqm7TuBiBYR0atE9Eciapf700R0q9y/Uq3BIo89QUSvEdFS+TepSmVsIKKb5f5XiOgd2jWnyP1riOjHFOMybjGWsZLPcToRPS7f3XIi+oLcP56IFhLRavnboV1zlXxerxHRBdr+ijzLmMtYkWdZbBmJaII8fz8R/dSXVk08x5Ay1spzfCcRLZHPawkRnaOlVbFvuyiihIKN5j8AZwI4GcAybd+LAM6S258EcI3c/iiAO+V2C4B1AGbK/58AMLcGyngFgJvl9iQASwAk5P8vADgdAAF4CMC7arCMlXyOUwCcLLfHAHgdwBwA3wVwpdx/JYDr5PYcAK8AaAQwC8AbAJKVfJYxl7Eiz7KEMrYCOAPAZwD81JdWrTzHoDLWynM8CcDBcvtYAJsr/RyL/TvgNRQhxJMAun27jwTwpNxeCOCD6nQArUSUAtAMYAjAvhor4xwAj8nrdsAJNZxLRFMAtAshFgmnBt4G4P21VMa4yhJQxq1CiJfkdg+AlQCmApgP4FZ52q3IP5f5cDoQg0KItQDWAJhXyWcZVxnjKEtcZRRC9AohngYwoKdTS8/RVsZKUkIZXxZCbJH7lwNoIqLGSn/bxXDACxQLywC8T25fDGC63L4bQC+ArXBGrH5PCKE3ojdLlfhrI6By2sr4CoD5RJQiolkATpHHpgLYpF2/Se6rpTIqKv4ciWgmnB7f8wAmCyG2As5HDkdrApzns1G7TD2zEXmWZZZRUdFnGbGMNmrpOYZRa8/xgwBeFkIMojrfthEWKGY+CeAKIloCRxUdkvvnAcgCOBiOeeFLRDRbHvuYEOI4AG+Xfx+vUhlvglOhFgO4HsCzADJwVGE/lQ7xK7aMwAg8RyJqA3APgH8VQgRpmLZnVvFnGUMZgQo/yyLKaE3CsK9azzGImnqORHQMgOsAfFrtMpxWlfBdFigGhBCrhBDnCyFOAfAbOHZpwPGhPCyEGJammmcgTTVCiM3ytwfA/6HyZgdjGYUQGSHE/xNCnCiEmA9gHIDVcBrwaVoS0wBs8adb5TJW/DkSURrOx3uHEOJeuXu7NBsoM8wOuX8TvJqTemYVfZYxlbGiz7LIMtqopedopZaeIxFNA/B7AJcKIVS7NOLftg0WKAZUFAcRJQD8B4BfyEMbAJxDDq0ATgOwSppuJspr0gDeC8fcM+JlJKIWWTYQ0TsBZIQQK6Tq3ENEp0mV/VIA99VSGSv9HOV93whgpRDiB9qh+wEskNsLkH8u9wO4RNqpZwE4HMALlXyWcZWxks+yhDIaqbHnaEunZp4jEY0D8ACAq4QQz6iTq/FtW4nby19vf3B6zlsBDMOR9JcB+AKciIvXAXwH+QGgbQB+B8chtgLAv4l8hMgSAH+Tx34EGWlThTLOBPAaHAffowAO0dKZC+djeAPAT9U1tVLGEXiOZ8AxBfwNwFL5924AE+AECayWv+O1a74qn9dr0CJnKvUs4ypjJZ9liWVcBydoY7+sH3Nq8DkWlLGWniOcTlmvdu5SAJMq/W0X88cj5RmGYZhYYJMXwzAMEwssUBiGYZhYYIHCMAzDxAILFIZhGCYWWKAwDMMwscAChWFqBCL6DBFdWsT5M0mb3Zlhqk2q2gVgGMYZQCeE+EX4mQxTu7BAYZiYkBP8PQxngr+T4AzovBTA0QB+AGdg7E4A/yiE2EpET8CZx+xtAO4nojEA9gshvkdEJ8KZWaAFzmC1TwohdhPRKXDmQusD8PTI3R3DhMMmL4aJlyMB3CCEOB7O0gZXAPgJgIuEM6fZTQCu1c4fJ4Q4SwjxfV86twH4d5nOqwCulvtvBvB5IcTplbwJhikF1lAYJl42ivw8S/8L4CtwFkNaKGc9T8KZokbxW38CRDQWjqD5q9x1K4DfGfbfDuBd8d8Cw5QGCxSGiRf/XEY9AJYHaBS9RaRNhvQZpmZgkxfDxMsMIlLC4yMAngPQqfYRUVquZ2FFCLEXwG4iervc9XEAfxVC7AGwl4jOkPs/Fn/xGaZ0WENhmHhZCWABEf0SzmyxPwHwCIAfS5NVCs6iYstD0lkA4BdE1ALgTQCfkPs/AeAmIuqT6TJMzcCzDTNMTMgorz8JIY6tclEYpiqwyYthGIaJBdZQGIZhmFhgDYVhGIaJBRYoDMMwTCywQGEYhmFigQUKwzAMEwssUBiGYZhYYIHCMAzDxML/D5NVj/nkEJwRAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'].plot()" ] @@ -215,11 +2242,34 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEKCAYAAADaa8itAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl83GW58P/Plclk39M0XdJSurC0LIWWUkCUtUV9FDyC1gWqh2MPiAs/fTyPHM/vwYNyjqhHOBwPKEqx4AIIR6msVhYRhZYCpfuSli7pmjZL0+wzcz1/fO9JJ+kkmSQzmcnker9eeWVyz/f7zT3fpnPNfV/3IqqKMcYYEw8Zya6AMcaY9GFBxRhjTNxYUDHGGBM3FlSMMcbEjQUVY4wxcWNBxRhjTNxYUDHGGBM3FlSMMcbEjQUVY4wxcZOZ7ArE25gxY3TKlCnJroYxxowob7311mFVrRjqddIuqEyZMoXVq1cnuxrGGDOiiMiueFzHur+MMcbEjQUVY4wxcRNzUBERn4i8IyJPu5/LRGSFiGxz30sjjr1NRKpFZIuILIwonyMi69xz94qIuPJsEXnMla8UkSkR5yx2v2ObiCyOx4s2xhiTGANpqXwV2BTx8zeBF1V1BvCi+xkRmQksAmYBVwH3iYjPnXM/sASY4b6ucuU3AvWqOh24G7jLXasMuB04H5gH3B4ZvIwxxqSWmIKKiFQBHwZ+HlF8NbDMPV4GXBNR/qiqtqvqe0A1ME9ExgNFqvq6epu4PNzjnPC1ngAud62YhcAKVa1T1XpgBccDkTHGmBQTa0vlHuCfgFBEWaWq7gdw38e68onAnojjalzZRPe4Z3m3c1Q1ADQC5X1cqxsRWSIiq0VkdW1tbYwvyRhjTLz1G1RE5H8Bh1T1rRivKVHKtI/ywZ5zvED1AVWdq6pzKyqGPMzaGGPMIMXSUrkI+KiI7AQeBS4TkV8CB12XFu77IXd8DTAp4vwqYJ8rr4pS3u0cEckEioG6Pq5lRpi/bKvlvcPNya6GMSbB+g0qqnqbqlap6hS8BPxLqvpZYDkQHo21GHjKPV4OLHIjuk7GS8ivcl1kTSIy3+VLbuhxTvha17rfocALwAIRKXUJ+gWuzIwwX310DT/7y45kV8MYk2BDmVH/PeBxEbkR2A1cB6CqG0TkcWAjEABuUdWgO+dm4BdALvCc+wJ4EHhERKrxWiiL3LXqROQ7wJvuuDtUtW4IdTZJEAiGqG/poKU9kOyqGGMSbEBBRVVfAV5xj48Al/dy3J3AnVHKVwNnRClvwwWlKM8tBZYOpJ4mtTS0dqIKbZ2h/g82xoxoNqPeJFx9cwcAbYFgP0caY0Y6Cyom4erCQaXTgoox6c6Cikm4+pZwULHuL2PSnQUVk3B1zZ2AtVSMGQ0sqJiEq2tuB6A9YC0VY9KdBRWTcNZSMWb0sKBiEu54TsWCijHpzoKKSbjjo7+s+8uYdGdBxSRcXcQ8FW/1HWNMurKgYhIuHFRUoSNorRVj0pkFFZNw9S0d+DK8XQysC8yY9GZBxSRUW2eQlo4g44pyAGi3ZL0xac2CikmocNfXhBIvqFhLxZj0ZkHFJFQ4qIwvzgVsUUlj0p0FFZNQ4Tkq47taKhZUjElnFlRMQoVbKhNLXEvFur+MSWv9BhURyRGRVSLyrohsEJF/deXfFpG9IrLGfX0o4pzbRKRaRLaIyMKI8jkiss49d6/bVhi39fBjrnyliEyJOGexiGxzX4sxI0pDi7dESzhRby0VY9JbLDs/tgOXqeoxEfEDr4lIeBvgu1X1h5EHi8hMvO2AZwETgD+JyCluS+H7gSXAG8CzwFV4WwrfCNSr6nQRWQTcBXxSRMqA24G5gAJvichyVa0f2ss2wyUcRErzswBbVNKYdNdvS0U9x9yPfvfV17Toq4FHVbVdVd8DqoF5IjIeKFLV19WbVv0wcE3EOcvc4yeAy10rZiGwQlXrXCBZgReIzAjR4YJIUY4fsJaKMekuppyKiPhEZA1wCO9NfqV76ksislZElopIqSubCOyJOL3GlU10j3uWdztHVQNAI1Dex7V61m+JiKwWkdW1tbWxvCQzTDqCITIE8rJ8gAUVY9JdTEFFVYOqOhuowmt1nIHXlTUNmA3sB/7DHS7RLtFH+WDPiazfA6o6V1XnVlRU9PlazPDqCITw+zLI9nt/am3W/WVMWhvQ6C9VbQBeAa5S1YMu2ISAnwHz3GE1wKSI06qAfa68Kkp5t3NEJBMoBur6uJYZIdoDIbIyM8jxey0Vm1FvTHqLZfRXhYiUuMe5wBXAZpcjCfsYsN49Xg4sciO6TgZmAKtUdT/QJCLzXb7kBuCpiHPCI7uuBV5yeZcXgAUiUuq61xa4MjNCdAZDZGdmkJNp3V/GjAaxjP4aDywTER9eEHpcVZ8WkUdEZDZed9RO4B8BVHWDiDwObAQCwC1u5BfAzcAvgFy8UV/hUWQPAo+ISDVeC2WRu1adiHwHeNMdd4eq1g3h9Zph1hEIkeXLwO8TMsTmqRiT7voNKqq6FjgnSvn1fZxzJ3BnlPLVwBlRytuA63q51lJgaX/1NKmpI+h1f4kIOX6ftVSMSXM2o94kVDhRD3hBxdb+MiatWVAxCdXhEvUAOZkZ1v1lTJqzoGISKtz9BVj3lzGjgAUVk1DhRD1Att9nLRVj0pwFFZNQ3VsqGbRbTsWYtGZBxSRUZEslJ9O6v4xJdxZUTEJ1S9T7LVFvTLqzoGISqtMS9caMKhZUTEJ16/6yeSrGpD0LKiahOoIh/Nb9ZcyoYUHFJFR75JBiS9Qbk/YsqJiE6gh4qxSD1/3Vbi0VY9KaBRWTUJ095ql0BEOEQn3tRm2MGcksqJiECQRDhJRuiXrwusSMMenJgopJmI6gFzz8EQtKgm3UZUw6s6BiEqbDtUh6tlRsWLEx6SuW7YRzRGSViLwrIhtE5F9deZmIrBCRbe57acQ5t4lItYhsEZGFEeVzRGSde+5et60wbuvhx1z5ShGZEnHOYvc7tonIYsyI0RVUMsMLSmZ0KzfGpJ9YWirtwGWqejYwG7hKROYD3wReVNUZwIvuZ0RkJt52wLOAq4D73FbEAPcDS/D2rZ/hnge4EahX1enA3cBd7lplwO3A+cA84PbI4GVSW7j7KxxUMjO8751BCyrGpKt+g4p6jrkf/e5LgauBZa58GXCNe3w18Kiqtqvqe0A1ME9ExgNFqvq6qirwcI9zwtd6ArjctWIWAitUtU5V64EVHA9EJsX17P4K7wDZEbDRX8akq5hyKiLiE5E1wCG8N/mVQKWq7gdw38e6wycCeyJOr3FlE93jnuXdzlHVANAIlPdxLTMC9GypZGUKAIGQtVSMSVcxBRVVDarqbKAKr9VxRh+HS7RL9FE+2HOO/0KRJSKyWkRW19bW9lE1M5x6a6lY95cx6WtAo79UtQF4Ba8L6qDr0sJ9P+QOqwEmRZxWBexz5VVRyrudIyKZQDFQ18e1etbrAVWdq6pzKyoqBvKSTAL1TNSHcyrW/WVM+opl9FeFiJS4x7nAFcBmYDkQHo21GHjKPV4OLHIjuk7GS8ivcl1kTSIy3+VLbuhxTvha1wIvubzLC8ACESl1CfoFrsyMAL11f1lLxZj0lRnDMeOBZW4EVwbwuKo+LSKvA4+LyI3AbuA6AFXdICKPAxuBAHCLqoYnJtwM/ALIBZ5zXwAPAo+ISDVeC2WRu1adiHwHeNMdd4eq1g3lBZvhE26p+K37y5hRo9+goqprgXOilB8BLu/lnDuBO6OUrwZOyMeoahsuKEV5bimwtL96mtQTDirhBSWPBxXr/jImXdmMepMwPbu/rKViTPqzoGIS5sTRX5ZTMSbdWVAxCdNpLRVjRh0LKiZhekvUd1hOxZi0ZUHFJEx7j3kq4W6wgLVUjElbFlRMwoQT9eHRX5mWUzEm7VlQMQnT+zwV6/4yJl1ZUDEJ0xkM4csQfBleCyU8+sv2UzEmfVlQMQnTEQh15VEARAS/T2yVYmPSmAUVkzAdgVBXkj7M78uw7i9j0pgFFZMwHcETg0pmhlj3lzFpzIKKSZj2Ht1f4A0vttFfxqQvCyomYTqD2kv3lwUVY9KVBRWTMB2B4AktFb8vg4DlVIxJWxZUTMJES9Rn+qRrUqQxJv1YUDEJEy1Rn2XdX8aktVi2E54kIi+LyCYR2SAiX3Xl3xaRvSKyxn19KOKc20SkWkS2iMjCiPI5IrLOPXev21YYt/XwY658pYhMiThnsYhsc1+LMSNGz3kqYEOKjUl3sWwnHAC+rqpvi0gh8JaIrHDP3a2qP4w8WERm4m0HPAuYAPxJRE5xWwrfDywB3gCeBa7C21L4RqBeVaeLyCLgLuCTIlIG3A7MBdT97uWqWj+0l22GQ0dQyc3qGVTEWirGpLF+Wyqqul9V33aPm4BNwMQ+TrkaeFRV21X1PaAamCci44EiVX1dVRV4GLgm4pxl7vETwOWuFbMQWKGqdS6QrMALRGYEiNZSybTuL2PS2oByKq5b6hxgpSv6koisFZGlIlLqyiYCeyJOq3FlE93jnuXdzlHVANAIlPdxLTMCdASCXSsUh2VZ95cxaS3moCIiBcCTwK2qehSvK2saMBvYD/xH+NAop2sf5YM9J7JuS0RktYisrq2t7fN1mOETLVFv3V/GpLeYgoqI+PECyq9U9X8AVPWgqgZVNQT8DJjnDq8BJkWcXgXsc+VVUcq7nSMimUAxUNfHtbpR1QdUda6qzq2oqIjlJZlh0BEIda1MHOb3ZdgyLcaksVhGfwnwILBJVX8UUT4+4rCPAevd4+XAIjei62RgBrBKVfcDTSIy313zBuCpiHPCI7uuBV5yeZcXgAUiUuq61xa4MjMCRJ1Rb8u0GJPWYhn9dRFwPbBORNa4sn8GPiUis/G6o3YC/wigqhtE5HFgI97IsVvcyC+Am4FfALl4o76ec+UPAo+ISDVeC2WRu1adiHwHeNMdd4eq1g3upZrh5iXqfd3K/BlCIGQ5FWPSVb9BRVVfI3pu49k+zrkTuDNK+WrgjCjlbcB1vVxrKbC0v3qa1NPr0vfW/WVM2rIZ9SYhQiGlIxgix39i91eHjf4yJm1ZUHHqmjv43EOreHHTwWRXJS20u9ZIdmb37i9bpsWY9GZBxcnKzOCVLbVUHzqW7KqkhfaAl0brOU8lM0MIWFAxJm1ZUHHys3xkZWZQ19yR7KqkhXBLJcffI1GfaZMfjUlnFlQcEaE8P4sjFlTioq0zekvF78ugIxjCGzFujEk3FlQilOVnWUslTrpyKv6ey7R4AwltWLEx6cmCSoQya6nETXun6/7qkaj3uwUmLVlvTHqyoBKhLD+LegsqcdGVqPefuEoxYHkVY9KUBZUI1v0VP22dvQ0p9rq/rKViTHqyoBKhPD+LY+2Brk/ZZvB6G1Js3V/GpDcLKhHK8rMBrLUSB70OKQ4HlYB1fxmTjiyoRCjLzwLgyDELKkPV6+RH1/3VYS0VY9KSBZUI5QVeUKlvsaAyVF05lROGFHs/B0IWVIxJRxZUIpTmeUHFur+Grt1Nfux1SLF1fxmTliyoRCi37q+46W3yo991h1n3lzHpyYJKhOJcP74MsZZKHPQ2pNhvQ4oHTFWpqW9JdjWMiUks2wlPEpGXRWSTiGwQka+68jIRWSEi29z30ohzbhORahHZIiILI8rniMg699y9blth3NbDj7nylSIyJeKcxe53bBORxSRQRoZQmue3WfVx0B4I4vcJvowT96gHCyoD8cy6/XzgB6+wv7E12VUxpl+xtFQCwNdV9XRgPnCLiMwEvgm8qKozgBfdz7jnFgGzgKuA+0Qk/HH1fmAJ3r71M9zzADcC9ao6HbgbuMtdqwy4HTgfmAfcHhm8EsFm1cdHeyB0QisFjgeVgM2oj9namkaCIWX7oeZkV8WYfvUbVFR1v6q+7R43AZuAicDVwDJ32DLgGvf4auBRVW1X1feAamCeiIwHilT1dfWWqH24xznhaz0BXO5aMQuBFapap6r1wAqOB6KEsFn18dEeCJ4wnBiOd39ZTiV24T1+dtdZF5hJfQPKqbhuqXOAlUClqu4HL/AAY91hE4E9EafVuLKJ7nHP8m7nqGoAaATK+7hWwpTnZ3OkuT2Rv2JUaOsMRQ0qWdb9NWDbDjUBFlTMyBBzUBGRAuBJ4FZVPdrXoVHKtI/ywZ4TWbclIrJaRFbX1tb2UbX+leZbTiUe2gOhE2bTQ+SCkhZUYtHSEaCm3sul7LGgYkaAmIKKiPjxAsqvVPV/XPFB16WF+37IldcAkyJOrwL2ufKqKOXdzhGRTKAYqOvjWt2o6gOqOldV51ZUVMTyknpVVZpHQ0snR9s6h3Sd0a69M0hWH91ftkpxbHbUNqMKvgxhj40AMyNALKO/BHgQ2KSqP4p4ajkQHo21GHgqonyRG9F1Ml5CfpXrImsSkfnumjf0OCd8rWuBl1ze5QVggYiUugT9AleWMNMqCgDvP7MZvLZAiOwoLRXr/hqYcD5lzkml1v1lRoRYWioXAdcDl4nIGvf1IeB7wJUisg240v2Mqm4AHgc2As8Dt6hqeNnfm4Gf4yXvtwPPufIHgXIRqQa+hhtJpqp1wHeAN93XHa4sYaZW5AOw3f1nNoPT3hkkJ2pLJTyj3oJKLLYdaiIzQ3j/jDE0tHTS2GotaJPaMvs7QFVfI3puA+DyXs65E7gzSvlq4Iwo5W3Adb1caymwtL96xsvksjwyM4Qdhy2oDEV7IERRrv+E8vCMeuv+is22g8c4qTyvqwW9p66F4onFSa6VMb2zGfU9+H0ZTC7PszkBQ+TNUznxzyszw4YUD0R17TGmjy1gUlkegM2sNynPgkoUU8cUWEtliNo7e5unYpMfB+LIsQ7GFeUwudwLKpZXManOgkoU08bms/NwC8GQvfENVm9Din0Z3tItlqiPTUtHgLzsTIpy/JTk+S2omJRnQSWKaWMK6AiGrKthCHqbUQ/esGILKv3rCIToDCr5WV5wHl+cy4HGtiTXypi+WVCJomsEWK11gQ1We2f0tb8A/BkZllOJQUtHAIC8LG88TWmen/oWG/1lUpsFlSjCI22eW3cAb7qMGai2QJAcfy8tlcwMa6nEoKXDG4mfn+0F59L8LNuV1KQ8CypRlOZn8Q/vO5nfvlXDv/x+fbKrM+IEQ0pnUHtvqfjEEvUxiNpSsSWETIqzoNKLb334dD59/mR+tXK3/UceoI5edn0M8/us+ysWze3dWypleVk0tnbaABKT0iyo9EJEWDCzEoCtB5uSXJuRpa1rf/rof15ZvoyuwGN619yjpVKSl0VI4ajNqjcpzIJKH04dVwjAFgsqA3J8f/ro3V9ZmRldx5jetbiWSp4b/VWWnwVgeRWT0iyo9GFcUQ6FOZlsOWBBZSDaA96bYW9DinP8vq7WjOndiS0Vb9kbCyomlVlQ6YOIcNq4Quv+GqCulkovifocf4YFlRj0HP0VbqnUNVv3l0ldFlT6cUplIVsONNnQ4gHoyqn0kqjP9fto67Tur/40t/cc/WXdXyb1WVDpx6njCjnaFuDgUdtiOFb9t1Ss+ysW4ZZKOKdSGs6p2GhEk8IsqPTjlEovWb/5QF87KJtI7Z19DynO8ftoC1hQ6U9Lh7d7ZngRzvwsH1m+DJtVb1KaBZV+nOqCiiXrY3d8SHHvLZXWDuv+6k9LR6Br3S/wcnwlNgHSpLhYthNeKiKHRGR9RNm3RWRvj50gw8/dJiLVIrJFRBZGlM8RkXXuuXvdlsK4bYcfc+UrRWRKxDmLRWSb+wpvNzysSvOzGFOQzTbbCTJm7f1MfszxZ9Bu3V/9am4PduVTwspsqRaT4mJpqfwCuCpK+d2qOtt9PQsgIjOBRcAsd859IhL+qHU/sARvz/oZEde8EahX1enA3cBd7lplwO3A+cA84Ha3T/2wmz423xaXHICYhhRb91e/WjoCXfmUsJI8vwUVk9L6DSqq+ioQ677wVwOPqmq7qr6Htxf9PBEZDxSp6uvqDaN6GLgm4pxl7vETwOWuFbMQWKGqdapaD6wgenBLuOljC6g+dMxGgMUo3FKJtp8KeN1inUElYEu19Km5I0hedrSWiuVUTOoaSk7lSyKy1nWPhVsQE4E9EcfUuLKJ7nHP8m7nqGoAaATK+7jWsJteUUBTW4DaJhsBFotwTqW3lkpullfeZrPq+9TS3j2nAt5SLZZTMalssEHlfmAaMBvYD/yHK5cox2of5YM9pxsRWSIiq0VkdW1tbV/1HpTpY71kfbXlVWISy5BiwIYV96O5I0pOJc/LqYRsUUmTogYVVFT1oKoGVTUE/Awv5wFea2JSxKFVwD5XXhWlvNs5IpIJFON1t/V2rWj1eUBV56rq3IqKisG8pD5NH+vtr1JteZWYhINFVm85lUwLKrFo6Qh0zaYPK8nzE1JoagskqVbG9G1QQcXlSMI+BoRHhi0HFrkRXSfjJeRXqep+oElE5rt8yQ3AUxHnhEd2XQu85PIuLwALRKTUda8tcGXDrrIom4LsTGupxKi1w9tK2JcRrbF5fFSYzarvW2+jvwDqLFlvUlRmfweIyG+AS4AxIlKDNyLrEhGZjdcdtRP4RwBV3SAijwMbgQBwi6qGP47ejDeSLBd4zn0BPAg8IiLVeC2URe5adSLyHeBNd9wdqhrrgIG4EhGmuWT9jtpjjCnMpijHn4yqjAgtHcETRi1FyrXur5j0nKcCx4PKkWPtnDwmPxnVMqZP/QYVVf1UlOIH+zj+TuDOKOWrgTOilLcB1/VyraXA0v7qOBymVxTwh7X7uOJHf+ajZ0/gnkXnJLtKKau188RP2JEsp9K/UEjdfeweVMYUZANw+JgNGjGpyWbUx+j08YV0BEKU5mWxYuNBe0PsQ2tHkNw+WirHg4p1f/WmLRBElROGFI8t9IJK7THr/jKpqd+WivFcf8FJnDO5hOb2IDcsXcWrW2tZMGtcsquVkqJN2ouU05VTscDcm66thKN0f4lgw9tNyrKWSoyyM33MOamMC6aVU5Ln57n1B5JdpZTV0hHsyptE05VTsVn1vWrpsUFXWKYvg7K8LOv+MinLgsoA+X0ZLJhZyZ82HuxajsR019oZW/dXa4fdv950tVSyT7yPFYXZ1lIxKcuCyiBcPKOCpvYA2w81J7sqKam/0V9dQ4ptRn2vemupgJest5aKSVUWVAbhpPI8APbUtyS5JqmptSNIrr//0V+2UnHvmjuspWJGJgsqgzC5zAWVOgsq0fSXqLd5Kv1rae+rpeLlVGyBU5OKLKgMQnGun8LsTAsqveiv+8vv82bbt1pQ6VVzj62EI1UUZtPWGeJYuy3VYlKPBZVBEBEmleWxp7412VVJOcGQ0h4I9ZmoB8jJzLB5Kn1odTmVaPexojA8AdLmqpjUY0FlkCaV5bLbWionCLc++mqpgNuoy1oqvTrqFoyMthxQeFa95VVMKrKgMkiTSvPYU9di/do9tHR9wu57Xq0XVKyl0puGlg5y/BlRNzo73lKxoGJSjwWVQZpcnkd7IGSfFnto6/ACRV4fkx/Bm1VvLZXeNbR0UpKbFfU5a6mYVGZBZZAmldqw4mhaOsOjlqz7aygaWzspyYu+EnZpXha+DOlqqbQHgnYvTcqwoDJIk7qGFVuyPlKLG7WUE0tQsRUJetXQ2klRbvSg4ssQyvOzuloqNz3yFp/46eu2G6RJCRZUBqmqNBfAkvU9hJde6a/7K9dyKn1qbOmkpJegAl5eZcfhZnYdaeblLbWsrWnkqXf3DmMNjYnOgsog5fh9VBZlW1DpoaVrfkV/ifoMW/urDw2tHb12fwF8+KzxrHqvjn/+3ToyBKZW5PPDF7bywKvbeWXLoWGsqTHd9RtURGSpiBwSkfURZWUiskJEtrnvpRHP3SYi1SKyRUQWRpTPEZF17rl73bbCuK2HH3PlK0VkSsQ5i93v2CYi4S2HU8YplYVs2n802dVIKS19zK+IlG3dX33ycirRE/UAf3/RyYwvzuGv1Ue49NSx3PHRM9jb0Mq/PbuZby/fMIw1Naa7WFoqvwCu6lH2TeBFVZ0BvOh+RkRm4m0HPMudc5+IhN9d7geW4O1bPyPimjcC9ao6HbgbuMtdqwxv6+LzgXnA7ZHBKxWcObGYLQeaLEkaobWPmeCRcjJ9tFv3V1RtnUHaOkMU99H9leP38fUFpwLwqXmTed+MMbz6jUtZdN4kmxRpkqrfoKKqr+LtHR/pamCZe7wMuCai/FFVbVfV94BqYJ6IjAeKVPV19SZ2PNzjnPC1ngAud62YhcAKVa1T1XpgBScGt6Q6q6qYQEjZfKAp2VVJGS0xBpXcLBtS3JvG1k6APoMKwMfPncjzt17M5aePBbxh7pPK8jjWHrB7a5JmsDmVSlXdD+C+j3XlE4E9EcfVuLKJ7nHP8m7nqGoAaATK+7hWyjizqgSAdTUNSa5J6gjPqO9/mRafrf3Vi4YWL6j0lVMBb7mg08YV4XqSAW9nSIC6ZmutmOSId6JeopRpH+WDPaf7LxVZIiKrRWR1bW1tTBWNhwnFOZTnZ7G2pnHYfmeqa+kI4MsQsnx9/2mF56nYigQnCrdUepv82JdyF1SOWBeYSZLBBpWDrksL9z083KQGmBRxXBWwz5VXRSnvdo6IZALFeN1tvV3rBKr6gKrOVdW5FRUVg3xJAycinFlVzLq9FlTCWjqC5Pl93T49R5PjzyCk0Bm0oNJTQ4sXEPrr/oqm3M22P9Jss+1Ncgw2qCwHwqOxFgNPRZQvciO6TsZLyK9yXWRNIjLf5Utu6HFO+FrXAi+5vMsLwAIRKXUJ+gWuLKWcObGYbYeO2fBYp7Wj762Ew3Jsn/peNbTG1v0VjbVUTLL1PZkAEJHfAJcAY0SkBm9E1veAx0XkRmA3cB2Aqm4QkceBjUAAuEVVw+8aN+ONJMsFnnNfAA8Cj4hINV4LZZG7Vp2IfAd40x13h6r2HDCQdLMmFBEMKdWHjnFmVXGyq5N0/e1PH9YVVDqCUVfiHc0aXU6leDBBpcAFFWupmCTpN6io6qd6eeryXo6/E7gzSvlq4Iwo5W24oBTluaXA0v7qmEwTSryZ9fsbWy2o4HV/5fYzmx4igooNKz5BY2snvgyhMLvf/54nKMg1xdSjAAAelElEQVTOJMuXwRFL1JsksRn1QzSuOAeAg0fbklyT1NDaz66PYTl+70/Pur9O1NDaQVFOZr95qWhEhPKCLOv+MkljQWWIxuRnk5kh7G+0oALh/en7/4Rd6Lq8mto6E12lEaehpe/Z9P0pL8iyIcUmaSyoDFFGhlBZlMMBa6kArvsrhpaKJZR719jaOaiRX2Fl+dkcsQ28TJJYUImDyqJsDlhLBfAS9bF0f5W6oFLfYkGlp772UonFmPwsW6rFJI0FlTgYX5xrLRWnJcacSlleeOa3dX/11NAy1JaKdX+Z5LGgEgeVRTkcaGyz2eG4eSr+/nMquVk+cvwZ1lKJoqGlo8+9VPpTXpBNa2ewa8VoY4aTBZU4GF+cQ0tHkKb20f2fWFVdor7/lgp4rRX7RN1da0eQo20BKgqzB30Ny1eZZLKgEgeVbljxaM+rtHQECSnkxzi/otS6aU6wt8HbnrqqNG/Q1zg+AdLurRl+FlTiYLwFFQD2uTfECSU5MR1vff8nqqn3dhKd6LarHoyu9b9sBJhJAgsqcTCuyIIKQI0LKhNLYntDLMvPspxKD3sHeA+jCXd/HbagYpLAgkocjC3yPhmO9hFge+vdG2KMn7JLLadygpr6VjLd3KfBCudjDh21oGKGnwWVOMjO9FGenzXqZ9Xva/DeEMcWxt791dQWoDNo63+F7a1vZXxJDr6MgS/REpbj91Gc6+dQkwUVM/wsqMTJ+JKcrv7w0Wpvw8DeELsmQFprpcvehlaqSgafpA+rLMq29ehMUlhQiZMzJhTz7p4GQqHRO1dlb30rE4pjzwV0TYC0vEqXmvqWISXpw8YW5lhLxSSFBZU4mXNSKUfbAlTXHkt2VZJmb0PrgN4QS/O9CX6WV/G0B4IcamqnKh5BpSibWgsqJgksqMTJ3CllALy5M+X2ERsWncEQB4+2UTWAUUvl+V5Cud6WagFgf0MbqkMb+RXmtVRslQcz/CyoxMmU8jzGFGTx1s76ZFclKQ40thHSgc2v6GqpWPcXEDGcOA4tlcqibDqDSn2LBWwzvIYUVERkp4isE5E1IrLalZWJyAoR2ea+l0Ycf5uIVIvIFhFZGFE+x12nWkTudfvY4/a6f8yVrxSRKUOpbyKJCHNOKuXNXaOzpXJ8fkXsSebScE7FlhMBjk98nDSE2fRh4RF4lqw3wy0eLZVLVXW2qs51P38TeFFVZwAvup8RkZl4+8/PAq4C7hOR8CJR9wNLgBnu6ypXfiNQr6rTgbuBu+JQ34Q5b0oZe+paR+V/5PAclVhn0wP4fRkU5mTaBEhn15EWfBnStZvoUFS6uVOWrDfDLRHdX1cDy9zjZcA1EeWPqmq7qr4HVAPzRGQ8UKSqr6vXAfxwj3PC13oCuFwGs8fqMJk9qQSAdTWNSa7J8NvbtUTLwLpubKmW47YcaGJaRT5+39D/W1pLxSTLUP96FfijiLwlIktcWaWq7gdw38e68onAnohza1zZRPe4Z3m3c1Q1ADQC5T0rISJLRGS1iKyura0d4ksavOljCwDYcXj0jQDbsK+RqtJccvyxrVAcVlGQPepXIgjbtP8op48visu1wqs82AgwM9yGGlQuUtVzgQ8Ct4jI+/s4NloLQ/so7+uc7gWqD6jqXFWdW1FR0V+dE6YkL4sxBVlsP9SctDokg6qy6r06zj/5hHjfr2kVBVQfGn1BuKfGlk72NbZx2rj4BJUcv4+inEwOWcA2w2xIQUVV97nvh4DfAfOAg65LC/f9kDu8BpgUcXoVsM+VV0Up73aOiGQCxUBKZ8KnVhSMurkq2w4do76lk/Onlg343BmVBdQ1d4z6FXU3HTgKwOnjC+N2zcqiHA6OkPW/th5s4l//sIHPPbSKV7cmr7dhJFu/tzElPqANOqiISL6IFIYfAwuA9cByYLE7bDHwlHu8HFjkRnSdjJeQX+W6yJpEZL7Ll9zQ45zwta4FXtIUH3g/faz3yTvFqxlXK3ccAWD+IFoqp1R6b6JbDyb/P0Mybd4fDirxaamA1wV2qCn1WyqhkPL5h97k1yt3s3HfUW5Yuor/frk62dUacf7t2U3c+tg7ya7GkFoqlcBrIvIusAp4RlWfB74HXCki24Ar3c+o6gbgcWAj8Dxwi6oG3bVuBn6Ol7zfDjznyh8EykWkGvgabiRZKptWUUBja+eoSD4HQ0p7IMgb79UxvjiHSWUDn18xo9LLQ2071BTv6o0om/Y3UZafxdgh7PjYU2VhzojYjmH1rnr2NrRy18fP4tV/upQLppbz65W7k12tESUUUtbWNHYNFkqm2Lboi0JVdwBnRyk/Alzeyzl3AndGKV8NnBGlvA24brB1TIZpFfkAbK9t7tosKV1977lNPPLGLgRh4axKBjMwb1xRDoXZmWw9OLqDyuYDRzl9fOGg7mFvplbk8z/v7OVYe4CCGHfjTIan1uwl1+/jypmV5Ph9vG/GGH7wwhaa2jopzPEnu3ojwvbaYxxrDzB7Umn/ByeYzaiPs/AIsO2jIK/y/IYDFOX4CamycNa4QV1DRJhRWTCqu786AiG2HGyKW5I+LHy9LQdSN2B3BEI8s24/V86s7NqGOtwlui0F8gMjxTt7GgCYPak4yTWxoBJ3E4pzyfFnpETCLJF2H2lhT10rX7psOlu++0E+eOb4QV/rlMpCth1sGjV5qHd21/P4m8dH1/9lWy1tnSEunDbwnFRfTnNJ/81uEEAq+vPWWhpaOvno2RO6yk4N59lSOBimmnf3NFCYncnUMQXJrooFlXjLyBCmjilI+09Zf91+GIALp40Z8rVmVBZS39LJ4VGyXMsdT2/kn55cy8ubvYGRy9/dR3Gun4tnxHc4/MSSXApzMtm0P3WDyi/f2EVlUTYfOPX4a68qzSXX72PLKO8SHYg1exo4a1IxGUPY3C1eLKgkwOzJJby1s472QLD/g0eo16oPM64opyuHNBSnjfM+ma7b2zDka6W6XUeaeWd3A5kZwj89uZZdR5r544aDfOjM8WRlxve/o4hw+rgiNu9PzTfnXUea+fPWWj41b3K3VQQyMoRTKguGPc/W2NrJl3/zDu/uGVl/h22dQTYfaEqJJD1YUEmIK0+vpLkjyBs7UnpKzaCFQsrfqg9z0fQxcUksz51SSkF2Ji+sPxiH2qW2p9Z4U7B+8tk5NLZ0cvl//JnWziBXz57Qz5mDc9r4QjYfSM2uxV+v3I0vQ1h03uQTnjulspAtB4a3tX/fy9X84d19fPFXb9PYOnJWd35x0yGCIU2JJD1YUEmIC6aVk+v38aeN6fkm+fbueupbOrl4xtC7vgCyM31cdtpY/rjxAIE03q9eVfn9mr3MO7mMK2ZW8sxX3selp43lvCmlzJsy8ImjsThtXBHH2gPUuAU/kyUYUvbUHd9uu60zyGOr97BgZmXUBTRPHVfI4WPtQ5oU2x4Isuq9upgCak19Cw/9bSfnTSnlwNE2/u9T6wf9e4dDbVM7dz2/mRUbD/Ivv1/HzPFFvP+U+Px/HCoLKgmQ4/dx8Ywx/GnTwZT8hDhUv11dQ16WNwQ0Xj54xjjqWzpZ9V56tu4AHnztPXbUNvPxc72l7WZUFvKzG+by25suTFhfeDhZn+y8yk/+vJ1LfvhK16CBZ9bup6Glk+vnnxT1+PAIsPX7BldvVeUbv13LJ376Ove+6E2kDPax1feP/rgVgP9cdA5funQ6T63Zx+vbjwzqdydaMKR85TfvcP8r2/nCw6tpD4T48afPITtzYOvuJYoFlQS5YmYl+xvbWLc3vVYsbukI8PTafXz4zPFdQ0Dj4ZJTx5Lr9/Hs+v1xu2YqefKtGr77zCY+dOY4rp0zqf8T4uS0cYXkZ/lY/u6+/g9OkI5AiF/8bSfBkPLDF7w371+u3MXUinwu6GXE25yTSinPz+K+l6ujfjAL7zQatq+hlZc3H+qadPzLlbtZ/u4+plbkc/eftnLJD15mxree5SP/9Rr3v7KdY+2BrnPX723kd2v28vcXncyEklxuvmQaE0ty+c7TG/sMRMly38vVvL7jCN+5ehY/uPYsln7uPKZWJH/UV5gFlQRZMLOSvCwfD/11Z7KrEjfvHW7mZ6++R3NHkOvmxveNMTfLx4JZlfzu7b0cTrN1wN7cWcc3/2ctF04r5+5PzsY3jCN08rIyWXzhFJ5Ztz9pE0yfW7+f2qZ23jfda73/8+/W8c7uBj57/km95uTyszO59YoZrHyvjhc3Her2nKryxV+9zSU/eIWa+hbue6WaC7/3Ep//xZvc9MhbbD3YxHee3sglp1bw3Fcv5rPzJ3NSeT43vu9ksjIzuOv5zVx810u8tcvrGvv35zZRkuvn5kumAV5Pwzc/eBob9x/la4+vSalFOQ8fa+fHL1fz4TPH89n5J3Hd3EnMnxrfoehDZUElQUrysvj0vMksf3cfu4+09H9Cint01W4u/eEr3P2nrUytyOe8KfFPCn7l8hm0BUL8+KX0Wfdp68EmbnrkLSaW5HL/Z+YkpYviCxdPJc/v4wcvbKEzCTmrh/66k6lj8vnJ9XMYV5TDo6t2c96UUq6dW9XneYvmTWZqRT7/9tymbvV+bv0BVmw8SGtnkK899i53r9jKFaeP5RsLT2XVzjquvf9v5Gf5+MG1Z5Od6eO715zJsr+fx7c+PJMnb76Qp265iJK8LD7/0Jt84eG3+Gv1Eb582QyKc4/P3v9fZ43nS5dO57l1B1hwz6u8vTsx24Q3tnRy/YMreXptbC3JZX/bSUcwxP935SlxXX0hniyoJNAX3j8Vnwj3/3l7sqsyJK9tO8y//H49F88Yw29vuoAnb7owIX/Q0yoK+MTcSfxq5S52pMGKBG/sOMLH7/sbGRnCzxefR3FecpYcKc3P4qYPTGPFxoMsuPtV1g9jl+zmA0dZs6eBz8w/iYLsTFZ87f2s+/ZCfnvThRT1swSL35fBbR88nR21zTzqJosePtbO7cs3cMbEIr5y+QxW7ayjIDuTuz5+Fl+8ZBpXnF7J0bYA373mTCp6WUft7EklPPz388jN8vHqtlr+6apT+dyFU7odIyL874Wn8tytF1Oc6+czP1uZkHzfXS9s5i/bDvPVR9fw4qa+B/Y0twd4+PVdXHl6ZdfKHanIgkoCVRbl8JGzJ/D02n0p2Tcbi4aWDm597B2mVuRz32fO5bwpZZTmZyXs9916xQzysjL54q/epqUj0P8JKaqtM8hXH32HiqJsfn/LRUl/E/jSZdN56HPn0doR5CuPvkNb5/DMoXpidQ1+n/Cxc7zBCYU5/gHl4q44fSzzTi7jnhVb2X2kha/85h2Otnby/Y+fzRcvmcZVs8bxg2vPprwgGxHhnkWz+fUXzufDZ/W9wsOksjye/vLF/Pkbl/DFS6b3OlBiWkUBv73pAsYV5/DVR9+J21DjrQebePj1nfx65W4+ff5kZk0o4ku/fof9jdFH6a3f28gnfvo6ja2d3OS66VKVBZUEu+TUCpraAiMyYR8KKXc8vZGGlk7u+eQ5w7K4X2VRDvd+6hy2HGzi///9hoT/vkT55Ru7OHi0nTuvOZOJA9xiORFEhEtPG8v3rz2LHbXN/NdL2xL+OzuDIX6/Zi+Xn1ZJ2SA/iIgI//Lh02ls7eT9P3iZv20/wnevOYOZE4rI8fv4yfVzuCJiFGJBdmbMqzxUFGYzvrj/f5uxhTnc/cnZHGpq58u/eYf/frma6iGsqv236sMsvOdV/u9TG5hakc+3PnQ6//3pcwmq8m/PbubuFVu58N9f5GP3/ZXfrt7DO7vrue4nr3OoqZ37P3Mu505OjfkovUndpUvTRHg9p79WH06ZGa/9CYWUL//mHZ5bv5+Qwpcvm87MCfFd7LAvHzilgiUXT+WBv+zgS5dN5+QxQ5+1n2iqyrH2AIU5fhpbOrn/le1cNL2819FNyfL+Uyq4dk4V//3ydnwZGXz18hkJGTjw2rbDPLt+P4ePdXDtnL5zJ/05q6qEFV/7ACs2HiDX74v7IJFYzJ5Uwv9ecCrff2Ezr26t5b6Xq/nWh2fS0NrBeVPKOK+feUb7Glr57jMbuXr2RO56fjOTSvN4cPFcJpfnkZ3pIz87k5veP5V7XT7xfdPHUNfcwTeeWEtWZgbjinJ48uYLe+3SSyWSbvMo5s6dq6tXr052Nbr54H/+hdI8P7/+wvxkV4WjbZ3kZPpoCwTZ19DKuppGQqpMLMlj/tQyMn0ZLPvbTm5fvoFPzp3EuSeV8HfnVnVbRmM41Da1c9H3XuIT51Xx3WvOHNbf3Z/dR1r41apdVBRkc805E9l+6Bj3/Gkbr+84whWnV7L1YBP7Glp54uYLU/KDRFtnkG/9bj1Pvl3DR8+eEPOItD9vreV3b9dQe6ydCcW5fPDMcVx22olzlZ5Zu59bfv02GeKtDffQ588b9r+fRGnrDHKkuYMlD69mg5tDkyHw5ctmMLUin1MqC7s2WusMhvhr9WHGF+dy62Nrus0Veujz53HpqWO7Xbu1I8iSR1Zz8YwxfOHiqYQU7n1xGy9sOMBPPjuHKQn+cCUib6nq3CFfx4JK4t35zEaW/W0X796+gNys5ExQamrr5OuPv8sf+5jlX1mUzZkTi3mt+jDzp5bz0OfOS+oIk//zxFp+v2Yvr/2fy1LiE1pHIMQP/7iFn/9lByLSLU9WmufnQ2eO56k1+yjMyeTHnz6XOSelbjeFqnL/n7fz/ee38Im5Vfz7353VZ2DZerCJj/zXaxTmZDKxNI89dS00tHTwyI3nM7Ekl1U762gPhGjtCHDPn7Zx2rhCfvUP85P2955obZ1B1u9tpKo0j3/9wwaeW3+g67mzqor5xsJT+fXK3V3lGQI/vX4u6/d6H+K+vuDUZFW9V6MqqIjIVcB/Aj7g56r6vd6OTcWg8sqWQ3zuoTf5ztWzuP6CKdQ1d1Ca5+/zDbs9EOSpNft4Y/sRzp5UwrmTS5k2Np+8rIH1WKoqr2yp5Y6nN7K7roV/eN/JFGRnkuP3MbYom1kTisnxZ7Bx31EeX72HfQ1tVBRm8/1rz6Ky6MTlM4ZT9aEmrrrnL1QW5XD7R2ZywbTyrryOqg5rwFu54wjf/sNGNu0/yqfmTeLWK07h4NE2/rLtMNPHFnDBtHKKcvwcbevEJxLXiaGJ9KM/buHel6q57LSx/N25E6ltaufy0yqZXJ5HIBjiL9sOs/NIM4+u2sPhY+08f+v7qSjMprk9wMfu+yt761tpC4S6Bdiq0lyevPnCpP/9DBdVZeeRFgKuZfLgX99jT52XcP/6ladQlOtnQkluXFegSIRRE1RExAdsxduauAZ4E/iUqm6MdnwqBpX2QJDP/nwlb+6sZ2JJLnsbWrnstLHc9IFpLHt9J83tAdo6g2zcd5QLp43hw2eN53vPbWZvQytFOZkcbfNGQfl9wmfOP4kPnFJBXXMHwSj/dhkinDwmv+v33L1iK69VH2ZKeR7//ndnpVwff3/e2lXH1x9/l51HWsgQ+LtzqzhzYjH3vriNjkCIiqJsKgtzGFuUTWVRDmMLs5k7pYyzq4qpbWpnT30L+xvbWFvTSGleFh8+czwb9zdyoLGNbL+PA41tZGYIc04qJS870705KlPHFNAWCPLYm3v406aDrN97lPHFOfzrR2exYJAbkqWqX76xi9uXb+gWGCqLsukMatcM9ezMDH7y2TlcetrxLpv3Djdz8y/f4sJpY/jM/MkUZmeSlZlBQXYmmWnS3TUYbZ1BHnztPcYUZPHJKItlpqrRFFQuAL6tqgvdz7cBqOq/Rzs+FYMKeP2r//1yNWtrGplSns+y171lK0ry/Ewuy0NEmDomn2fX7ac9EGJaRT63f2QWF88Yw566Vjbsa+SVLbU88XbNgIYnF+f6+erlM/js/JPivrT6cGnr9BYGfHnLIX75xi46g8q8k8uYOb6IQ01tHDrazkH3vT3gTZIrzMmkqe34kOQsXwYdUSb+hRs70f4bZAgoMPekUhbOGsdnzj8pbbtzdtQeo6UjSFGOn+fW72dHbTOdoRALZ41j3pQycvy+tH3txjOagsq1wFWq+g/u5+uB81X1SxHHLAGWAEyePHnOrl27klLXgVi54wird9Xz2fkndZvJu6P2GK/vOMLHz60ix3/if+Ka+hYOHm2jPD+bTN+J3T+dQWXrwSYOH2unIDuTD5xSQUle4uaVDLddR5rZ29DKBVPLT+j+UvU+Wa/YeJC3dtVz2vgiplbkM7YwmxljC9ld18wrW2o5q6qEaRX5tAdClBdk0dYRYu3eBjqDITJEUIUtB5to6wzy8XOrmFSWl6RXa8zwGU1B5TpgYY+gMk9Vvxzt+FRtqRhjTCqLV1AZCf0hNUDkwPQqIHlLrhpjjOnVSAgqbwIzRORkEckCFgHLk1wnY4wxUaT8uEdVDYjIl4AX8IYUL1XVkbt+hzHGpLGUDyoAqvos8Gyy62GMMaZvI6H7yxhjzAhhQcUYY0zcWFAxxhgTNxZUjDHGxE3KT34cKBFpArYM4RJjgMNxqEoxEM+dueJ9vURcM173LizV7+Foun9279L7emOAfFWtGPKVVDWtvoDVyTw/4joPxPl1xfV6CapjXO7dSLmHo+n+2b1L++vF7f5Z91fi/CHFr5eoa8ZTqt/D0XT/7N6l9/XiJh27v1brENavGer5o5ndu6Gx+zd4du+GJp73Lx1bKg8k+fzRzO7d0Nj9Gzy7d0MTt/uXdi0VY4wxyZOOLRVjjDFJkvZBRUQmicjLIrJJRDaIyFddeZmIrBCRbe57qSsvd8cfE5Ef97jWp0RknYisFZHnRWRMMl7TcInzvfuku28bROT7yXg9w20Q9+9KEXnL/Y29JSKXRVxrjiuvFpF7pecOZWkmzvfuThHZIyLHkvV6hlu87p+I5InIMyKy2V3ne/3+8ngOS0vFL2A8cK57XIi33/1M4PvAN135N4G73ON84H3ATcCPI66TCRwCxrifv4+3zXHSX+MIuHflwG6gwv28DLg82a8vBe/fOcAE9/gMYG/EtVYBFwACPAd8MNmvbwTdu/nueseS/bpG2v0D8oBL3eMs4C/9/e0l/cUn4WY/BVyJN0FyfMQ/wJYex32uxxujH6gFTnL/sX8CLEn26xkh9+484E8RP18P3Jfs15Oq98+VC3AEyHbHbI547lPAT5P9ekbCvetRPmqCSiLun3vuP4Ev9PW70r77K5KITMGLyCuBSlXdD+C+j+3rXFXtBG4G1uHtPDkTeDCB1U0pQ7l3QDVwmohMEZFM4Bq67+aZ9gZx/z4OvKOq7cBEvB1Qw2pc2agwxHs36sXr/olICfAR4MW+ft+oCSoiUgA8CdyqqkcHcb4fL6icA0wA1gK3xbWSKWqo905V6/Hu3WN4zeedQCCedUxlA71/IjILuAv4x3BRlMNGxbDNONy7US1e9899GPwNcK+q7ujrGqMiqLiA8CTwK1X9H1d8UETGu+fH4+VL+jIbQFW3q9cOfBy4MEFVThlxuneo6h9U9XxVvQCvCb4tUXVOJQO9fyJSBfwOuEFVt7viGqAq4rJVeK3ltBanezdqxfn+PQBsU9V7+vu9aR9U3CiZB4FNqvqjiKeWA4vd48V4fY592QvMFJHwgmtXApviWddUE8d7h4iMdd9LgS8CP49vbVPPQO+f6154BrhNVf8aPth1UzSJyHx3zRuI4Z6PZPG6d6NVPO+fiHwXbwHLW2P65clOIA1Dgup9eF0Fa4E17utDeCOSXsT7xPwiUBZxzk6gDjiG9ylxpiu/CS+QrMVbe6c82a9vBN273wAb3deiZL+2VLx/wL8AzRHHrgHGuufmAuuB7cCPcROX0/Urzvfu++5vMeS+fzvZr2+k3D+8VrG6971w+T/09bttRr0xxpi4SfvuL2OMMcPHgooxxpi4saBijDEmbiyoGGOMiRsLKsYYY+LGgooxw0xEbhKRGwZw/BQRWZ/IOhkTL5nJroAxo4mIZKrqT5JdD2MSxYKKMQPkFuh7Hm+BvnPwlhW/ATgd+BFQABwGPqeq+0XkFeBvwEXAchEpxFsx94ciMhtvxes8vImNf6+q9SIyB1gKtACvDd+rM2ZorPvLmME5FXhAVc8CjgK3AP8FXKuq4YBwZ8TxJar6AVX9jx7XeRj4P+4664DbXflDwFfUWyvNmBHDWirGDM4ePb5G0i+Bf8bb3GiF25TRB+yPOP6xnhcQkWK8YPNnV7QM+G2U8keAD8b/JRgTfxZUjBmcnusbNQEb+mhZNA/g2hLl+saMCNb9ZczgTBaRcAD5FPAGUBEuExG/25uiV6raCNSLyMWu6Hrgz6raADSKyPtc+WfiX31jEsNaKsYMziZgsYj8FG/F1/8CXgDudd1XmcA9wIZ+rrMY+ImI5AE7gM+78s8DS0WkxV3XmBHBVik2ZoDc6K+nVfWMJFfFmJRj3V/GGGPixloqxhhj4sZaKsYYY+LGgooxxpi4saBijDEmbiyoGGOMiRsLKsYYY+LGgooxxpi4+X9FfYAQnZLJSwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ - "sorted_data['inc'][-200:].plot()" + " sorted_data['inc'][-200:].plot()" ] }, { @@ -252,13 +2302,11 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 12, + "metadata": {}, "outputs": [], "source": [ - "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", + " first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", " for y in range(1985,\n", " sorted_data.index[-1].year)]" ] @@ -274,7 +2322,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, "outputs": [], "source": [ @@ -298,11 +2346,34 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 17, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 17, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "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='*')" + " yearly_incidence.plot(style='*')" ] }, { @@ -314,11 +2385,57 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 18, "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": 18, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "yearly_incidence.sort_values()" + " yearly_incidence.sort_values()" ] }, { @@ -331,9 +2448,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 19, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 19, + "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)" ] @@ -364,7 +2504,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4,