From 6ddbc61529e63ca2de172898b2beb7c0b6486254 Mon Sep 17 00:00:00 2001 From: 7c8b98dd05a86d272dc44ef63cc416c0 <7c8b98dd05a86d272dc44ef63cc416c0@app-learninglab.inria.fr> Date: Tue, 22 Jun 2021 10:16:19 +0000 Subject: [PATCH] test_exo1 --- module3/exo1/analyse-syndrome-grippal.ipynb | 2210 ++++++++++++++++++- 1 file changed, 2173 insertions(+), 37 deletions(-) diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5..1d0b5a5 100644 --- a/module3/exo1/analyse-syndrome-grippal.ipynb +++ b/module3/exo1/analyse-syndrome-grippal.ipynb @@ -9,7 +9,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -28,10 +28,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 2, + "metadata": {}, "outputs": [], "source": [ "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" @@ -61,9 +59,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
0202123375215069.09973.0117.015.0FRFrance
1202122378795495.010263.0128.016.0FRFrance
2202121378275403.010251.0128.016.0FRFrance
32021203102787540.013016.01612.020.0FRFrance
4202119395396860.012218.01410.018.0FRFrance
52021183121359165.015105.01814.022.0FRFrance
62021173120588891.015225.01813.023.0FRFrance
720211631650512735.020275.02519.031.0FRFrance
820211531930615398.023214.02923.035.0FRFrance
920211432107317099.025047.03226.038.0FRFrance
1020211332641322094.030732.04033.047.0FRFrance
1120211233065825919.035397.04639.053.0FRFrance
1220211132498820718.029258.03832.044.0FRFrance
1320211031953915951.023127.03025.035.0FRFrance
1420210931757213926.021218.02721.033.0FRFrance
1520210832088216907.024857.03226.038.0FRFrance
1620210732239318303.026483.03428.040.0FRFrance
1720210632318319134.027232.03529.041.0FRFrance
1820210532242618445.026407.03428.040.0FRFrance
1920210432580421491.030117.03932.046.0FRFrance
2020210332181017894.025726.03327.039.0FRFrance
2120210231732013906.020734.02621.031.0FRFrance
2220210132179917778.025820.03327.039.0FRFrance
2320205332122016498.025942.03225.039.0FRFrance
2420205231642812285.020571.02519.031.0FRFrance
2520205132161917370.025868.03327.039.0FRFrance
2620205031684513220.020470.02620.032.0FRFrance
272020493129399923.015955.02015.025.0FRFrance
2820204831380410641.016967.02116.026.0FRFrance
2920204731908515285.022885.02923.035.0FRFrance
.................................
188119852132609619621.032571.04735.059.0FRFrance
188219852032789620885.034907.05138.064.0FRFrance
188319851934315432821.053487.07859.097.0FRFrance
188419851834055529935.051175.07455.093.0FRFrance
188519851733405324366.043740.06244.080.0FRFrance
188619851635036236451.064273.09166.0116.0FRFrance
188719851536388145538.082224.011683.0149.0FRFrance
18881985143134545114400.0154690.0244207.0281.0FRFrance
18891985133197206176080.0218332.0357319.0395.0FRFrance
18901985123245240223304.0267176.0445405.0485.0FRFrance
18911985113276205252399.0300011.0501458.0544.0FRFrance
18921985103353231326279.0380183.0640591.0689.0FRFrance
18931985093369895341109.0398681.0670618.0722.0FRFrance
18941985083389886359529.0420243.0707652.0762.0FRFrance
18951985073471852432599.0511105.0855784.0926.0FRFrance
18961985063565825518011.0613639.01026939.01113.0FRFrance
18971985053637302592795.0681809.011551074.01236.0FRFrance
18981985043424937390794.0459080.0770708.0832.0FRFrance
18991985033213901174689.0253113.0388317.0459.0FRFrance
190019850239758680949.0114223.0177147.0207.0FRFrance
190119850138548965918.0105060.0155120.0190.0FRFrance
190219845238483060602.0109058.0154110.0198.0FRFrance
1903198451310172680242.0123210.0185146.0224.0FRFrance
19041984503123680101401.0145959.0225184.0266.0FRFrance
1905198449310107381684.0120462.0184149.0219.0FRFrance
190619844837862060634.096606.0143110.0176.0FRFrance
190719844737202954274.089784.013199.0163.0FRFrance
190819844638733067686.0106974.0159123.0195.0FRFrance
19091984453135223101414.0169032.0246184.0308.0FRFrance
191019844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1911 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202123 3 7521 5069.0 9973.0 11 7.0 \n", + "1 202122 3 7879 5495.0 10263.0 12 8.0 \n", + "2 202121 3 7827 5403.0 10251.0 12 8.0 \n", + "3 202120 3 10278 7540.0 13016.0 16 12.0 \n", + "4 202119 3 9539 6860.0 12218.0 14 10.0 \n", + "5 202118 3 12135 9165.0 15105.0 18 14.0 \n", + "6 202117 3 12058 8891.0 15225.0 18 13.0 \n", + "7 202116 3 16505 12735.0 20275.0 25 19.0 \n", + "8 202115 3 19306 15398.0 23214.0 29 23.0 \n", + "9 202114 3 21073 17099.0 25047.0 32 26.0 \n", + "10 202113 3 26413 22094.0 30732.0 40 33.0 \n", + "11 202112 3 30658 25919.0 35397.0 46 39.0 \n", + "12 202111 3 24988 20718.0 29258.0 38 32.0 \n", + "13 202110 3 19539 15951.0 23127.0 30 25.0 \n", + "14 202109 3 17572 13926.0 21218.0 27 21.0 \n", + "15 202108 3 20882 16907.0 24857.0 32 26.0 \n", + "16 202107 3 22393 18303.0 26483.0 34 28.0 \n", + "17 202106 3 23183 19134.0 27232.0 35 29.0 \n", + "18 202105 3 22426 18445.0 26407.0 34 28.0 \n", + "19 202104 3 25804 21491.0 30117.0 39 32.0 \n", + "20 202103 3 21810 17894.0 25726.0 33 27.0 \n", + "21 202102 3 17320 13906.0 20734.0 26 21.0 \n", + "22 202101 3 21799 17778.0 25820.0 33 27.0 \n", + "23 202053 3 21220 16498.0 25942.0 32 25.0 \n", + "24 202052 3 16428 12285.0 20571.0 25 19.0 \n", + "25 202051 3 21619 17370.0 25868.0 33 27.0 \n", + "26 202050 3 16845 13220.0 20470.0 26 20.0 \n", + "27 202049 3 12939 9923.0 15955.0 20 15.0 \n", + "28 202048 3 13804 10641.0 16967.0 21 16.0 \n", + "29 202047 3 19085 15285.0 22885.0 29 23.0 \n", + "... ... ... ... ... ... ... ... \n", + "1881 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1882 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1883 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1884 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1885 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1886 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1887 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1888 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1889 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1890 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1891 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1892 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1893 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1894 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1895 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1896 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1897 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1898 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1899 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1900 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1901 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1902 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1903 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1904 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1905 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1906 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1907 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1908 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1909 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1910 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 15.0 FR France \n", + "1 16.0 FR France \n", + "2 16.0 FR France \n", + "3 20.0 FR France \n", + "4 18.0 FR France \n", + "5 22.0 FR France \n", + "6 23.0 FR France \n", + "7 31.0 FR France \n", + "8 35.0 FR France \n", + "9 38.0 FR France \n", + "10 47.0 FR France \n", + "11 53.0 FR France \n", + "12 44.0 FR France \n", + "13 35.0 FR France \n", + "14 33.0 FR France \n", + "15 38.0 FR France \n", + "16 40.0 FR France \n", + "17 41.0 FR France \n", + "18 40.0 FR France \n", + "19 46.0 FR France \n", + "20 39.0 FR France \n", + "21 31.0 FR France \n", + "22 39.0 FR France \n", + "23 39.0 FR France \n", + "24 31.0 FR France \n", + "25 39.0 FR France \n", + "26 32.0 FR France \n", + "27 25.0 FR France \n", + "28 26.0 FR France \n", + "29 35.0 FR France \n", + "... ... ... ... \n", + "1881 59.0 FR France \n", + "1882 64.0 FR France \n", + "1883 97.0 FR France \n", + "1884 93.0 FR France \n", + "1885 80.0 FR France \n", + "1886 116.0 FR France \n", + "1887 149.0 FR France \n", + "1888 281.0 FR France \n", + "1889 395.0 FR France \n", + "1890 485.0 FR France \n", + "1891 544.0 FR France \n", + "1892 689.0 FR France \n", + "1893 722.0 FR France \n", + "1894 762.0 FR France \n", + "1895 926.0 FR France \n", + "1896 1113.0 FR France \n", + "1897 1236.0 FR France \n", + "1898 832.0 FR France \n", + "1899 459.0 FR France \n", + "1900 207.0 FR France \n", + "1901 190.0 FR France \n", + "1902 198.0 FR France \n", + "1903 224.0 FR France \n", + "1904 266.0 FR France \n", + "1905 219.0 FR France \n", + "1906 176.0 FR France \n", + "1907 163.0 FR France \n", + "1908 195.0 FR France \n", + "1909 308.0 FR France \n", + "1910 213.0 FR France \n", + "\n", + "[1911 rows x 10 columns]" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data = pd.read_csv(data_url, skiprows=1)\n", "raw_data" @@ -78,9 +1043,73 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
167419891930NaNNaN0NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1674 198919 3 0 NaN NaN 0 NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1674 FR France " + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -94,9 +1123,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 5, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
0202123375215069.09973.0117.015.0FRFrance
1202122378795495.010263.0128.016.0FRFrance
2202121378275403.010251.0128.016.0FRFrance
32021203102787540.013016.01612.020.0FRFrance
4202119395396860.012218.01410.018.0FRFrance
52021183121359165.015105.01814.022.0FRFrance
62021173120588891.015225.01813.023.0FRFrance
720211631650512735.020275.02519.031.0FRFrance
820211531930615398.023214.02923.035.0FRFrance
920211432107317099.025047.03226.038.0FRFrance
1020211332641322094.030732.04033.047.0FRFrance
1120211233065825919.035397.04639.053.0FRFrance
1220211132498820718.029258.03832.044.0FRFrance
1320211031953915951.023127.03025.035.0FRFrance
1420210931757213926.021218.02721.033.0FRFrance
1520210832088216907.024857.03226.038.0FRFrance
1620210732239318303.026483.03428.040.0FRFrance
1720210632318319134.027232.03529.041.0FRFrance
1820210532242618445.026407.03428.040.0FRFrance
1920210432580421491.030117.03932.046.0FRFrance
2020210332181017894.025726.03327.039.0FRFrance
2120210231732013906.020734.02621.031.0FRFrance
2220210132179917778.025820.03327.039.0FRFrance
2320205332122016498.025942.03225.039.0FRFrance
2420205231642812285.020571.02519.031.0FRFrance
2520205132161917370.025868.03327.039.0FRFrance
2620205031684513220.020470.02620.032.0FRFrance
272020493129399923.015955.02015.025.0FRFrance
2820204831380410641.016967.02116.026.0FRFrance
2920204731908515285.022885.02923.035.0FRFrance
.................................
188119852132609619621.032571.04735.059.0FRFrance
188219852032789620885.034907.05138.064.0FRFrance
188319851934315432821.053487.07859.097.0FRFrance
188419851834055529935.051175.07455.093.0FRFrance
188519851733405324366.043740.06244.080.0FRFrance
188619851635036236451.064273.09166.0116.0FRFrance
188719851536388145538.082224.011683.0149.0FRFrance
18881985143134545114400.0154690.0244207.0281.0FRFrance
18891985133197206176080.0218332.0357319.0395.0FRFrance
18901985123245240223304.0267176.0445405.0485.0FRFrance
18911985113276205252399.0300011.0501458.0544.0FRFrance
18921985103353231326279.0380183.0640591.0689.0FRFrance
18931985093369895341109.0398681.0670618.0722.0FRFrance
18941985083389886359529.0420243.0707652.0762.0FRFrance
18951985073471852432599.0511105.0855784.0926.0FRFrance
18961985063565825518011.0613639.01026939.01113.0FRFrance
18971985053637302592795.0681809.011551074.01236.0FRFrance
18981985043424937390794.0459080.0770708.0832.0FRFrance
18991985033213901174689.0253113.0388317.0459.0FRFrance
190019850239758680949.0114223.0177147.0207.0FRFrance
190119850138548965918.0105060.0155120.0190.0FRFrance
190219845238483060602.0109058.0154110.0198.0FRFrance
1903198451310172680242.0123210.0185146.0224.0FRFrance
19041984503123680101401.0145959.0225184.0266.0FRFrance
1905198449310107381684.0120462.0184149.0219.0FRFrance
190619844837862060634.096606.0143110.0176.0FRFrance
190719844737202954274.089784.013199.0163.0FRFrance
190819844638733067686.0106974.0159123.0195.0FRFrance
19091984453135223101414.0169032.0246184.0308.0FRFrance
191019844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1910 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202123 3 7521 5069.0 9973.0 11 7.0 \n", + "1 202122 3 7879 5495.0 10263.0 12 8.0 \n", + "2 202121 3 7827 5403.0 10251.0 12 8.0 \n", + "3 202120 3 10278 7540.0 13016.0 16 12.0 \n", + "4 202119 3 9539 6860.0 12218.0 14 10.0 \n", + "5 202118 3 12135 9165.0 15105.0 18 14.0 \n", + "6 202117 3 12058 8891.0 15225.0 18 13.0 \n", + "7 202116 3 16505 12735.0 20275.0 25 19.0 \n", + "8 202115 3 19306 15398.0 23214.0 29 23.0 \n", + "9 202114 3 21073 17099.0 25047.0 32 26.0 \n", + "10 202113 3 26413 22094.0 30732.0 40 33.0 \n", + "11 202112 3 30658 25919.0 35397.0 46 39.0 \n", + "12 202111 3 24988 20718.0 29258.0 38 32.0 \n", + "13 202110 3 19539 15951.0 23127.0 30 25.0 \n", + "14 202109 3 17572 13926.0 21218.0 27 21.0 \n", + "15 202108 3 20882 16907.0 24857.0 32 26.0 \n", + "16 202107 3 22393 18303.0 26483.0 34 28.0 \n", + "17 202106 3 23183 19134.0 27232.0 35 29.0 \n", + "18 202105 3 22426 18445.0 26407.0 34 28.0 \n", + "19 202104 3 25804 21491.0 30117.0 39 32.0 \n", + "20 202103 3 21810 17894.0 25726.0 33 27.0 \n", + "21 202102 3 17320 13906.0 20734.0 26 21.0 \n", + "22 202101 3 21799 17778.0 25820.0 33 27.0 \n", + "23 202053 3 21220 16498.0 25942.0 32 25.0 \n", + "24 202052 3 16428 12285.0 20571.0 25 19.0 \n", + "25 202051 3 21619 17370.0 25868.0 33 27.0 \n", + "26 202050 3 16845 13220.0 20470.0 26 20.0 \n", + "27 202049 3 12939 9923.0 15955.0 20 15.0 \n", + "28 202048 3 13804 10641.0 16967.0 21 16.0 \n", + "29 202047 3 19085 15285.0 22885.0 29 23.0 \n", + "... ... ... ... ... ... ... ... \n", + "1881 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1882 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1883 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1884 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1885 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1886 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1887 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1888 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1889 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1890 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1891 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1892 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1893 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1894 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1895 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1896 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1897 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1898 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1899 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1900 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1901 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1902 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1903 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1904 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1905 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1906 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1907 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1908 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1909 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1910 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 15.0 FR France \n", + "1 16.0 FR France \n", + "2 16.0 FR France \n", + "3 20.0 FR France \n", + "4 18.0 FR France \n", + "5 22.0 FR France \n", + "6 23.0 FR France \n", + "7 31.0 FR France \n", + "8 35.0 FR France \n", + "9 38.0 FR France \n", + "10 47.0 FR France \n", + "11 53.0 FR France \n", + "12 44.0 FR France \n", + "13 35.0 FR France \n", + "14 33.0 FR France \n", + "15 38.0 FR France \n", + "16 40.0 FR France \n", + "17 41.0 FR France \n", + "18 40.0 FR France \n", + "19 46.0 FR France \n", + "20 39.0 FR France \n", + "21 31.0 FR France \n", + "22 39.0 FR France \n", + "23 39.0 FR France \n", + "24 31.0 FR France \n", + "25 39.0 FR France \n", + "26 32.0 FR France \n", + "27 25.0 FR France \n", + "28 26.0 FR France \n", + "29 35.0 FR France \n", + "... ... ... ... \n", + "1881 59.0 FR France \n", + "1882 64.0 FR France \n", + "1883 97.0 FR France \n", + "1884 93.0 FR France \n", + "1885 80.0 FR France \n", + "1886 116.0 FR France \n", + "1887 149.0 FR France \n", + "1888 281.0 FR France \n", + "1889 395.0 FR France \n", + "1890 485.0 FR France \n", + "1891 544.0 FR France \n", + "1892 689.0 FR France \n", + "1893 722.0 FR France \n", + "1894 762.0 FR France \n", + "1895 926.0 FR France \n", + "1896 1113.0 FR France \n", + "1897 1236.0 FR France \n", + "1898 832.0 FR France \n", + "1899 459.0 FR France \n", + "1900 207.0 FR France \n", + "1901 190.0 FR France \n", + "1902 198.0 FR France \n", + "1903 224.0 FR France \n", + "1904 266.0 FR France \n", + "1905 219.0 FR France \n", + "1906 176.0 FR France \n", + "1907 163.0 FR France \n", + "1908 195.0 FR France \n", + "1909 308.0 FR France \n", + "1910 213.0 FR France \n", + "\n", + "[1910 rows x 10 columns]" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -122,7 +2118,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 6, "metadata": {}, "outputs": [], "source": [ @@ -152,10 +2148,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 7, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -179,9 +2173,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n" + ] + } + ], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", @@ -199,9 +2201,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXe4HcV5/7/vaberdwmQAFEEmKbIEIptOrZjnATnB4lBcXBwbFxwEifgEhJsHJw4LrgRbGQDxsYY26GDRTVgIRCIJgmhhgpqV7oqV7eeMr8/dmbP7jkzO3vO2dPufT/Po+ceze7Mzs7uzjtvmRkSQoBhGIZhKiVW7wowDMMwIwMWKAzDMEwksEBhGIZhIoEFCsMwDBMJLFAYhmGYSGCBwjAMw0QCCxSGYRgmEligMAzDMJHAAoVhGIaJhES9K1BLJk2aJGbPnl3vajAMwzQVL7300i4hxGTbeaNKoMyePRvLli2rdzUYhmGaCiLaGOY8NnkxDMMwkcAChWEYhokEFigMwzBMJLBAYRiGYSKBBQrDMAwTCVaBQkSLiGgnEb3hSZtARIuJaI38O95z7FoiWktEq4nofE/6yUT0ujx2ExGRTG8hol/J9KVENNuTZ6G8xhoiWuhJnyPPXSPzpipvCoZhGKYSwmgoPwNwQUHaNQAeF0LMBfC4/D+IaB6ASwAcI/P8kIjiMs+PAFwJYK78p8q8AsAeIcThAL4N4BuyrAkArgPwbgALAFznEVzfAPBtef09sgyGYRimjlgFihDiDwB6CpIvAnCb/H0bgA970u8SQgwJITYAWAtgARFNBzBGCLFEOHsO316QR5V1D4CzpfZyPoDFQogeIcQeAIsBXCCPnSXPLbz+iKanbxgPvb6t3tVgGIbRUq4PZaoQYhsAyL9TZPpMAJs9522RaTPl78J0Xx4hRAbAPgATA8qaCGCvPLewrBHNx297EZ+682XsPjBU76owDMMUEbVTnjRpIiC9nDxBZRVXiOhKIlpGRMu6u7tNpzUFb+/uBxBwswzDMHWkXIGyQ5qxIP/ulOlbABzkOW8WgK0yfZYm3ZeHiBIAxsIxsZnK2gVgnDy3sKwihBC3CCHmCyHmT55sXYqmoclkcwCAREwnUxmGYepLuQLlPgAq6mohgHs96ZfIyK05cJzvL0izWC8RnSJ9IJcX5FFlXQzgCelneRTAeUQ0XjrjzwPwqDz2pDy38PojmmzO0U0EqygMwzQg1sUhieiXAN4LYBIRbYETeXUjgLuJ6AoAmwB8BACEECuI6G4AKwFkAFwlhMjKoj4JJ2KsDcDD8h8A3ArgDiJaC0czuUSW1UNEXwXwojzveiGECg74VwB3EdHXACyXZYx40lKg5FiiMAzTgFgFihDiUsOhsw3n3wDgBk36MgDHatIHIQWS5tgiAIs06evhhBKPKlwNpc71YBiG0cEz5ZuILGsoDMM0MCxQmhGWJwzDNCAsUJqQHAsUhmEaEBYoTQibvBiGaURYoDQhLE4YhmlEWKA0ITm2eTEM04CwQGkiiCfIMwzTwLBAaSJiUqKwD4VhmEaEBUoToRQUtngxDNOIsEBpQgRrKAzDNCAsUJoI5UNhDYVhmEaEBUoTQa7RiyUKwzCNBwuUJoI1FIZhGhkWKE1EXqCwRGEYpvFggdJEKJNXI8uTbE5g9jUP4kdPrat3VRiGqTEsUJqIZtBQ0nKb4m8/9lada8IwTK1hgdJEuC75xpUnDMOMYligNBFEjW/yYhhm9MICpYnIz5RvAonSBFVkGCZaWKA0E03gQ2EYZvTCAoVhGIaJBBYoTQjrJwzDNCIsUJqIZojyauS6MQxTXVigNBHN1FeLpqotwzBRwAKlKWnczpoFCcOMXligMJHCJi+GGb2wQGlCmqHTboY6MgwTLSxQmohm2A2lkevGMEx1YYHSRDRTZ91MdWUYJhpYoDQhJnOSEALv/vpj+NWLm2pboYI6MAwzOqlIoBDR54loBRG9QUS/JKJWIppARIuJaI38O95z/rVEtJaIVhPR+Z70k4nodXnsJpKrIBJRCxH9SqYvJaLZnjwL5TXWENHCSu5jpCAEsGP/EP71N6/Xrw51uzLDMPWmbIFCRDMBfBbAfCHEsQDiAC4BcA2Ax4UQcwE8Lv8PIponjx8D4AIAPySiuCzuRwCuBDBX/rtApl8BYI8Q4nAA3wbwDVnWBADXAXg3gAUArvMKrpFOM2gBzVBHhmGipVKTVwJAGxElALQD2ArgIgC3yeO3Afiw/H0RgLuEEENCiA0A1gJYQETTAYwRQiwRTi90e0EeVdY9AM6W2sv5ABYLIXqEEHsALEZeCI1YbE75RujCWY4wzOilbIEihHgHwDcBbAKwDcA+IcTvAUwVQmyT52wDMEVmmQlgs6eILTJtpvxdmO7LI4TIANgHYGJAWaOahtAKGqAKDMPUh0pMXuPhaBBzAMwA0EFEHw3KokkTAenl5ims55VEtIyIlnV3dwdUr3kwOuVrW41AGqkuDMPUhkpMXucA2CCE6BZCpAH8FsCfAtghzViQf3fK87cAOMiTfxYcE9kW+bsw3ZdHmtXGAugJKKsIIcQtQoj5Qoj5kydPLvNWGwNbJ90YCkoDVIJhmLpQiUDZBOAUImqXfo2zAawCcB8AFXW1EMC98vd9AC6RkVtz4DjfX5BmsV4iOkWWc3lBHlXWxQCekH6WRwGcR0TjpaZ0nkwbFZg6be7MGYapJ4lyMwohlhLRPQBeBpABsBzALQA6AdxNRFfAETofkeevIKK7AayU518lhMjK4j4J4GcA2gA8LP8BwK0A7iCitXA0k0tkWT1E9FUAL8rzrhdC9JR7L82Ca+czzkOpVU3MNEIdGIapD2ULFAAQQlwHJ3zXyxAcbUV3/g0AbtCkLwNwrCZ9EFIgaY4tArCoxCozVYblCcOMXnimfBPSDJ02ayoMM/pggdKEmJdeqW099HVogEowDFMXWKCMIBrBKV//GjAMUy9YoDQhxigv7s0ZhqkjLFBGEI0gT6IQaoPpLC67dSlWb++tvDCGYWoGC5QmJGj5+pHA8k178cyaXfjKvW/UuyoMw5QAC5QmpKEXh4ywFrr1dRiGaVxYoIwgGkJBaYQ6MAxTF1igNCFG0xZ35gzD1BEWKEykRCHTGiH8mWGY0mGB0oSYfSj174ijNLsRO1EYpqlggdKMNPBM+UgYKffBMKMMFigjiEboh6ON8mIVhWGaCRYoTYh5pnwjiBSGYUYrLFBGEI0gTqKQaY1wHwzDlA4LlCakoVcbjrAsdsozTHPBAqUJMU9DaQCJwjDMqIUFChMp7MdhmNELC5QmxNhlN0BfHokPRZbBJi+GaS5YoDQRJHtYkxbQAPIkEpTpjsOGGaa5YIEygmBrE8Mw9YQFShPSyEuvMAwzemGB0kTYHN6NoKE0Qh0YhqkPLFCakEZevT4KLYmFEsM0JyxQmghyw55G9tIr6i44yothmgsWKEyk1Eqm7e0fRt9QpjYXYxgmFCxQSuCNd/Zh8cod9a7GqFl6JYgTrl+M9/z3kzW6GsMwYUjUuwLNxAe/9ywA4O0bP1Dnmoxswprudh0YrnJNGIYpBdZQmhBj2HAjaCgRVoLYicIwTQULlCaEF4dkGKYRYYEygmgIDaXeFWAYpm5UJFCIaBwR3UNEbxLRKiI6lYgmENFiIloj/473nH8tEa0lotVEdL4n/WQiel0eu4mkrYOIWojoVzJ9KRHN9uRZKK+xhogWVnIfzYZxx8Ya10MHb7DFMKOXSjWU7wJ4RAhxFIDjAawCcA2Ax4UQcwE8Lv8PIpoH4BIAxwC4AMAPiSguy/kRgCsBzJX/LpDpVwDYI4Q4HMC3AXxDljUBwHUA3g1gAYDrvIJrpGOO8hohXbFabbi+tWAYpkTKFihENAbAmQBuBQAhxLAQYi+AiwDcJk+7DcCH5e+LANwlhBgSQmwAsBbAAiKaDmCMEGKJcHrE2wvyqLLuAXC21F7OB7BYCNEjhNgDYDHyQmjU0hjiJEqnfGRFMQxTAyrRUA4F0A3gp0S0nIh+QkQdAKYKIbYBgPw7RZ4/E8BmT/4tMm2m/F2Y7ssjhMgA2AdgYkBZo4LGEBx6RoqSxDBM6VQiUBIATgLwIyHEiQD6IM1bBnTjTRGQXm4e/0WJriSiZUS0rLu7O6B6zQ935gzD1JNKBMoWAFuEEEvl/++BI2B2SDMW5N+dnvMP8uSfBWCrTJ+lSfflIaIEgLEAegLKKkIIcYsQYr4QYv7kyZPLuM3Gw+wrqb9EqX8NGIapF2ULFCHEdgCbiehImXQ2gJUA7gOgoq4WArhX/r4PwCUycmsOHOf7C9Is1ktEp0j/yOUFeVRZFwN4QvpZHgVwHhGNl87482TaqGakaCg8n4ZhmpNKl175DIA7iSgFYD2Aj8ERUncT0RUANgH4CAAIIVYQ0d1whE4GwFVCiKws55MAfgagDcDD8h/gOPzvIKK1cDSTS2RZPUT0VQAvyvOuF0L0VHgvoRFCNOQs7kbohiPdU77yohiGqSEVCRQhxCsA5msOnW04/wYAN2jSlwE4VpM+CCmQNMcWAVhUSn2jIieAeB17u8ZeHJKXXmGY0QrPlC+DbK4Bem4NbCoanfT0DWP2NQ/i3lfeqXdVmFEOC5QyqLdAMc6UbwB5EqXJiwnH+u4DAIDb/vh2fSvCjHpYoJRBts49XiN3uFHWjQ1ewPPrd2P3gaHAc9gyyDQKLFDKoN4aiolGFjRM6QghcMktz+PSHz8f7vwq14dhbLBAKYNcvU1eDbx8fSPUYaSQke/ZWzsOWM6sv4qSzQksXrlj5Kwnx5QFC5QyqLfJy0SDVqtkGuU2lm/ag+vvX1m362eypbVEPZ//j59Zj7+/fRkeeWN7/SrB1B0WKGVQdw2lrlcPJopOLScLqbdv4C9/9Ecsem4D0tlcXa6fzoW7br3bCQA2dPcBAPYOpOtcE6aesEApg3ppKKrjMJkVRoyG4t5HfXvKloSzu8JAOms5szqUrKFUqR5hGJZCNxVv7i5l30Aas695ED9/fmO9q9KUNPfTrxP1csrbBMbI8V8E30eUdvqevmH0DupH1a1J5/MYHK6XQAmpoVS5HmEYzkiBkmjuLkVF1P3kmfV1rklz0txPv06EtERUDePSkA0gT2oxDyVKeX7SVxfjjP96Unus3hrKcKmmtjq+AEMjRKC0pZxn3lenQUSz09xPv05k6iRRXFt5AwgOE1FoSbYSoo4k2tuv11BalIaSrs/zDmvyaoQlapTfK94AdYmC/qFMvavQlLBAKYNcI6gCGhqzVqXjLg5p6JtqdZ/xmFOBejnlMyWqYiPl+dcT1eQNOtWs4WGBUgZ16l9czEuv1P8riMTkZfWhVH6NMCSkQCm1Y/fSO5jGknW7y8wdUkMps/QoaYQ6RIGK4Bw5/sjawgKlDBp2pny9KxARqnlNnVStPvZEzPk8KtFQPv2L5bj0x89jT99wyXlLFZyVCNqtewfw9FuV72g6Ut5BpjxYoJRBvU1ejb18fQRlWG6kZhqK3KMgnSlfoKze3gsAGMyU7uQNe5tRuC0+cNMzWLjohbLzjxDXifttN8K31IywQCmDIA3l5qfXYcOuvqpe33z1+n8FtTC7RXUJW11dH0oFGqkqoxytttSBSyWa2x5DYEKp1Nvs+smfv4QbHix/dQMWJJXBAqUMTBMb9w+mcePDb+LSW8It5hc1I+VjsDvlo7lRWzhwUpm8KtBQ1D2UExgY9nk2xnNvDBXl4Te248fPbCg7f72tD80OC5QyMC29ot7FA1UOOTQvDll/IjF51cgp3zcULFCkPKkoTNzVUMqodNgsjWSmaYAqVIT6tJv9PuoFC5QyMJov1MtYpy+7ITqUmkxsjOZG+4eDBT/JUXclMRgxKt/kFVYTa4DH7lkWSH/8rR29dQu/Lo1GaM3mhQVKGZhGmyq92kFgIz2kMR/lpTejRHX3acvEQddcVYEAi1nWXwsivMkruveh3IVP80+qOP/q7b0479t/wPefWFt2vWpFgwZwNg0sUMrANNpU6dW2w5qjvBrha4hgprxlteFaOeXzAqX8ayiTVzU7qigfe7nvbpCGsmrbfgDAum7bvi71x73/RviUmhAWKCWgRps2gVKvd3GkfAPW+4joRsN28pUIamXyKscPE96HUtr5YcoqFZM2CeTXJFNro9WCcp9ZQ4zJmhgWKCWgOgfTKC7rOkerrKGY0hvgY4ikDpYyojL5hR2NV6JxKg2l1KXogRJ8KFGavCosS5dbfTe1NNWWKxjdAIcRMzyrLSxQSiDvYNUfz7kmr1rVyE8jfAS1iPKKqn1tnacadVfSx1bihylZQyn5CrqyKjN5aY+pHzV8Pcu9j0YYlDUzLFBKgMKavKr9VjZy3HAEWPd9iah9S+2wK6GcMsJ2ilEOJCq916Aq1/L1ZIFSH1iglIDN5KUWEaz2cuJGk1dVrxqOaBaHdKj2asO2ukYR5eW5WtVyuH7kCOpZsVNeU+t6LMtSblM00pyeZoQFSgnYnPLuXui1qlABjfARRNmpGcOGa2Tyyl8vinsqPU89ZsqLCqeK6Opi27q6GpS7gGsDfEJNDQuUErA65eVLHKu2hjLC3/pambzCO+XLv4Y7ObKsQqofNFBIOTP6Ae8kUI2GonxR5VerZMptE156pTJYoJRAWB9KtVUUU4c6cpzy1b8GYBcUZBlARHktHaE1lNKLNlL2vYZ452vZV5c7CFDfVv2/pOaEBUoJxCwbLtV7n5QRM7gSwYK5VhMbFVE81rJmyoc8L0q7f9k+lID8ef9K7eB5KPWhYoFCRHEiWk5ED8j/TyCixUS0Rv4d7zn3WiJaS0Srieh8T/rJRPS6PHYTyaEhEbUQ0a9k+lIimu3Js1BeYw0RLaz0PsKg9ss2mS+UuaDaPpTR4pQ3H4/K5BV8XD3HKExsZRm8Qnvlyyi80muWkb+WPpTy56E4fxtj1YnmIwoN5XMAVnn+fw2Ax4UQcwE8Lv8PIpoH4BIAxwC4AMAPiUhNnf0RgCsBzJX/LpDpVwDYI4Q4HMC3AXxDljUBwHUA3g1gAYDrvIKrWigTiMnOnKuRD8VEI3wEUXT2dh9KxZcAEGIeirv0fH2ip8JrUMpMU3k9y9Wy8+bBiqsQCexDqQ8VCRQimgXgAwB+4km+CMBt8vdtAD7sSb9LCDEkhNgAYC2ABUQ0HcAYIcQS4XxBtxfkUWXdA+Bsqb2cD2CxEKJHCLEHwGLkhVDViFk6mHzYcHXrMcKnoVg70qg++po45StYD6zUsOEoqHimvCZ/PZbHqnQeSrVD/0cqlWoo3wHwLwC8wYZThRDbAED+nSLTZwLY7Dlvi0ybKX8XpvvyCCEyAPYBmBhQVlWxLUWuBE21XkXXBGM6wXPgydU7q1QLC5HY8Z2/xj3lI/OhBB8P8guUSjlllLo0TCTmxgp98rr89Vhwsdz7cJ3yrKmURdkChYg+CGCnEOKlsFk0aSIgvdw8/osSXUlEy4hoWXd3d6iKmnDnoRjeNdeHUqXRTSm+hY/99MWq1MFGlFFe1R4l2k1eESy9Iv+W1UHVIcqr0sASndktL09q6UMpU0OJuB7l8tjKHXh2za56V6NkKtFQTgPwISJ6G8BdAM4iop8D2CHNWJB/1VB5C4CDPPlnAdgq02dp0n15iCgBYCyAnoCyihBC3CKEmC+EmD958uTy7lRCNqd8lTUUxUgfPdnuLzofStjzyr+gylnWFsBhz4sw1LXSmfK6Nq3H7PNKF4esNx+/fRk+euvSelejZMoWKEKIa4UQs4QQs+E4258QQnwUwH0AVNTVQgD3yt/3AbhERm7NgeN8f0GaxXqJ6BTpH7m8II8q62J5DQHgUQDnEdF46Yw/T6ZVFbUlrHGDrRrNQzHRCN9CLeoQ1Udv3Q/FvV4E1yonT1gNJVIfSnn5gkxe+aVhyiu7HMoNpOAtgCujGvNQbgRwLhGtAXCu/D+EECsA3A1gJYBHAFwlhFCben8SjmN/LYB1AB6W6bcCmEhEawH8I2TEmBCiB8BXAbwo/10v06qKzYdSKw3FRBQf7N7+YTz0+rbKC6oAa5RXja6jqESAVeKHCWsiijLUNaiMLXv6jceCJoGq+6ilyatSHwpTHokoChFCPAXgKfl7N4CzDefdAOAGTfoyAMdq0gcBfMRQ1iIAi8qtczlYnfLyZVQTIKMmaBQIRNPRfu6uV/D0W9149l/fh1nj20vOH0nYsKWMWi+9Uu2O2pwn5HkRdtQm7fv3K7bjyjtewqK/nY+zjppadDzIV1SPUGJebbg+8Ez5ErAvvSLPq1F9Comi49u2bwAA0DeUtZxpqkPFVbBNlK/Z0iuu/6OSC1YwPyNslkj3QzH4el7ZvBcAsHLr/sD8ujrUx4dSrsmr9nUdSbBAKYOhjP6rU9u8Vn/5ev3bHsU3kIw7r0TatItYDXBdUcalV2qrodQrbDjsfdZix8b8MzE8lIA5WvWZh1JePhYklcECpRTky9Y/nNEertXy9UaTVwQfgxIow2UKlGjChmsT5WWPJotutFqT5esj1A6L06U51zJY0mWPsh3DUq6QbZQor2aFBUoJqFfNZA5yTV51m2Rb+ceQUhqKQQuz1iASf4PleMVXcAhr8qqbDyXsnvIRtIhrzrVEMMYNPUZ++friY/XxoZSXj8VJZbBAKQHVKQymTQJFdcJ12rExCg0l4dS9XA0lSqo9U942GlWdUjSrDVcvj3rtKnL1qLIsJi+rhqJ1ykeoQoWkfKc8i5RKYIFSAupVyxg8l/XXUCrHFslmIxKTl80UVaPVhlU96hY2HDrKq3LyqwIERzCafCj5XRmLj9VlHkrZJq+IKzLKYIFSAuodNQ3elbmgSlHDRfUoSo+gbCVQ6vld2fqCcmad669j01CUQKn8WtWN8lI+iiiCBwx1kWXHg33yWmGfXw25unjvv/x5KBFVZpTCAqUE1MeSNfRo+cUh6xTlFcHHQPmeoTyiMA9Zj9cmyss1JTV4lFcUzRGzhcRb5lgFLb2S11Cq21t7i+fl6+sDC5QyMO3YaEqvFZE4ZyssK4o62MwrkflQLJpOXkOpgZSsIEsUGkDQnvBOujzP6kPRpKHy+oXBW/dmX8urWWGBUgJ5k5fho3M32KpNPcKml0J+AczKyyoX28TGqMgLruB6RGPyamwfiu2a+bBh/fEggZRfGqbi6gXifU713o57tMICpQTUB2HcU94ysg7LwHAW+wbSJeeLxodSWVlRDuZNRZXSOfcPZ/D3ty/D5p7idaisvpoonPIVbLAV9ilEMqq2rgKhfCilv9s186F4rsDzUOoDC5QysH10lXL2/zyF4//j9yXni8ZGHWz6sNchgipYnMylXOPxVTuxeOUO3PjIm0XHbBNRXcFWJw0l7OsUoTwpO2w4aLvk+vhQKi+DKR0WKCWgPgijhhLRFsBb9w2Gqkc1cDWUOn5Zdqd8eNQyOS2J4lfd5heIREORf6u6OGSEM9FtWqHt3dblD3Pv2ZzAud96Grf8YZ31XOO1I3HKl315BixQSkK9a9YNtpp4HkrQfIIwRKigBEzgDH+VoYwzCbUlES86Fn5iY/l35b4zZTnlw2XKmwgrb32bVmi6Qt60Z/ahBDGUyWLNzgP4+kPFmmRY/E756pq8Fi56Af+l0XpHOyxQSsDmQ3FHtAEO7dXbe7HrwFAk9QibXgrKuVrPgZptxdeSNJS0WUNRnad5Rr7SUEq4oKWs0vKEO6/czaRKuaYaLNnn7ujS7BpU1JuYlf0thMz39Fvd+OFTpWtTuZzAJbcswVOrd9pPbkJYoJSAbR5KJsRHd/53/oD3/vdTkdcNiGaEqnalrOfSFTanfCnXUKsmJzUz8myrGkcxYbCSXR/DZonC12MLHgg9yVNTibx2Y84chf+xlhpKufQNZ/D8+h58+hfLq3qdesECpQSsGkou3Ed3YEi/WnHoepjSI9RQbPfwyua9+sipyqtgdeKWcp/qWcVjOh+K0lAMPpSc/28lVHf5+pKLNmJ2yqvBVPAz0U9stGvuYe917c5ePLd2V2AdTPUIQ7V9KEE7WyqaeT2xSHZsHC2ox2yL8qr2KKe681CcvzYzyod/8BwA4O0bP1D5RQuwLl9fQlk5V6BojlkKiiZsWK2RVXYRVqLcFMpURiYbrK3lQ4PNPhTTSsZAeA3lnG/9AYD+vRMRaCjV3qa4Vua/esEaSgm4GkrWIFAi6ID81yutnCiuGmYEFViHSFQU35+KrqE+Tt38Cbd9K5jY+NDr23DNb14LUY/q+VCixPTOqQ7ftlupVkNRM+UDbiiKSYnebOXPQykrW2iEuzJ0dc1/9YIFSkkEC4xsSJNXWIzmBeNaXhH4UCwT3OxE50MxFVVK5/zC27sB6EODbfNQgkbdik/d+TLuenGztR7VjPKKUiM2lZS2vNvKr6iriquhBDSC9x527A8Omzfh01DKNFNW29yUH3Saz2nmyZUsUErA5kOJ2uRVaCKwmU8i0VBUWXV8p4WlI/fZygO+TMfe7giUhGbNkLAmr2jmd5ReSNhOsRY+FCUwjMeVVqkNGxa+c2zXjUJDqaaGbRM6T7y5A69s3qs9FiZarpkFCvtQSsDmQ8mHDUfzQhR2KLVw1imhFWTvDiKKKtp3UvR0PkIgZtAx9g3kgx90q+SGX8srAq2rLA0lHLmAzrzkaxqKSGdt2rnUUALKDKqf95sqfy8erw+lrCJCPWvbIrB/97NlAIL9PEGXYZPXKMGdKW/YECU/+ojmeqZO3Vh8BNetZEOoqMh3QKYT8j+DPj6vVhLXCJS8C8U2U95cVxtue5ZRSOgorxAPftnbPXhrR6/1vHLNue5x7dIrwf4XJ1/+d7mrdnurXq5wDZMtHbCbqWk3V0UYP2s9F2atFNZQSsCmoagPwSQISu1UCq9jyx3N3uJqteFyR4mVk3fi2q8R1El5hYjOKW+7xyhmyheWVQphs+TneZi5+OYlACqPyrNFOOqO5tvRXG4Uc0iiCRu2Z0wH2O627h2wlO/8DbpKM5u8WEMpgbDzUGybFIXF2OFFMD/DRGWr45YYgZUTOOubT+H+V7dqyzAJSJ+9PeCCOq3EX47zN4qJjXbhVD2JEsVWxQpbGbYoMP3SK3ZTsPeMWE0qAAAgAElEQVQ5mqIobUQtlEwEayjB6oVqg0CTFwuU0YFNdVffgc1sEJbCF8s2Eo3iNQwT5RWVj2gwk8X6XX34wj2vao8bNRTvSDSkhhLU0ZlwfSghTBCmQYYo+FsK4aO8nL9R2N5tfZmpL80EdJRhfFFeQRXmPnSCzZtSzXkoQQIl6BgQ7t6iXEqn1rBAKQHX5GVxTEa1vH3hi2UbKUeioYRYvj5oBFWK2c3kK7Hep+d3kL3dW4yuzraaljKx0RaoEXRPS9fvxqMrthfXL7SG4vyNYsfQsJFvpnRtR+9GeQW8U55+OMwIXbtmmCex3G8hTBMGaVDDFoES5l2KYj5NvWCBUgqeDlA3ivA6LnUvQqmqrFFDMY3cI/GhOH+DXvygD6qUW7QFHRg1Ma9pI8iMYulgbKaiUpzyGYMaE2Z0/v9ueR6fuOOl4rz2y/rKLtdU5Lum0Zwa3FbuTHpt/VQZ5uv6w4btKqGpvXXllYI3n6ktgoTGcKbyemVF8HvbyLBAKQHbyNg3ytIcr7pTPgoNRYUNB3wXQR9zKR+ya082lBHG5BU0Ks9ZBI/tOqWE45qapJIVi0NrKPKvzdwShLvKtKUt7MLX3M5BGrr3WBjBqP2+fD4UaxFa/JFi+nOC2tmmoYR5RN53tdkc9CxQSsBm57U5i0s1SZhMQcYJfyWVrieMhhJVnLxbTtF9ur+0+XzzUEKavLQmEv3lPfnDm7zSJg3FvVbpbRY2j6pnJQLFdk2bYMzkzMfDtEHYQIvC6/nL0JdXCmHWAws0eUWgoZTaFo0EC5QS8D5am8DQ9S+VaihhTAdhGM7k8Iulm7T1UX7soLqG9VvYyBoEZL4Dsl8jbCcVNHI2L3gYXA/TtfTXsJdRSNgswlPPSoW9KbfbFpYIR919uqtwh9RQynVce+eHlR32bhmEAMFaiK3uoe4thJbUqJQtUIjoICJ6kohWEdEKIvqcTJ9ARIuJaI38O96T51oiWktEq4nofE/6yUT0ujx2E0m7CxG1ENGvZPpSIprtybNQXmMNES0s9z5KQYj8ZLmsZpTifYl1AqdiH4qtiwlZ/s1Pr8MXf/c6frv8naJjYZavD9QKStCT7EuhBx8HgoWb3xZdfJ7NJ1WKU95k8lLpYcxmReeEfJ7e+lWqpZRj0gKC9wIK44vymbxCdLq6c7zzQ8qfh5L/bXqX0wFaiK39wzzSKFYNqBeVaCgZAP8khDgawCkAriKieQCuAfC4EGIugMfl/yGPXQLgGAAXAPghEal9WX8E4EoAc+W/C2T6FQD2CCEOB/BtAN+QZU0AcB2AdwNYAOA6r+CqFgICSbkOus6PYHsRKo/yUvUw1S8cu+WOkb2DaeM5QcIv6KPxNoutEzW3hwg87k0NGon6zRfFx/P715Q+6i7EFmAQ5tEXTpgLraF4ftts+NayjMLV+WsqPmiVCJUn6J3ya/f2O9e9G95v0vRM13UfwOxrHsRLG3u0x/1Oef21lxvW6TLVy3e8xIGF6T7e3tWHO5dutJZVa8oWKEKIbUKIl+XvXgCrAMwEcBGA2+RptwH4sPx9EYC7hBBDQogNANYCWEBE0wGMEUIsEU5L3l6QR5V1D4CzpfZyPoDFQogeIcQeAIuRF0JVIyeARFw5rYM1kKAosLCUHOUVsnh1mm4+X5ilxoM1lDzWEFR39B4uPV9usCaou77e5JX/G6TBhNNQTNpUeC2ncJAS9nn6NBSLDd+E7bnb7qPSiY3eY1FoKKb7eHaNsznX/y3fqj3uf3/1Zdz4sHkveVvdw/QBvn7EcPrFNy/Bl373hnEZqHoRiQ9FmqJOBLAUwFQhxDbAEToApsjTZgLwrvO9RabNlL8L0315hBAZAPsATAwoS1e3K4loGREt6+7uLu8GJULkNRStSctm8qo4yivcR29DnaZb0l0dC6pr0GiylE2O8j4UfXook1eAg9RmgrRFBZVk8rLUNcyjKdJQSh0haMooxGZmNOW2rZJrepZAuHb0vlOhOl3Nffp8KIYi1Pdr0rJtWq21XpZMpUYMmoTwvoFhAMBQmQOIalGxQCGiTgC/AXC1EGJ/0KmaNBGQXm4ef6IQtwgh5gsh5k+ePDmgenaEyO9NruvIsgEdLWAPYy06v2i14eLr+I5bS1TnOWfqlhxRZYS1dxdGN4UZ4eXL0X8MtiVsvFcJOwFTd5r349aZMEtxyltXoA6joWTNbRlEKT4UY0eqyrLch3lSr/k+bZtzFZYbTkMpvg/fQMdQT/X9mgRvpcu3eJ+hdi5aKKe8vQ5qFQjbYpS1piKBQkRJOMLkTiHEb2XyDmnGgvy7U6ZvAXCQJ/ssAFtl+ixNui8PESUAjAXQE1BWVcl5NZQywoa933KYj6awDNtHHXpAqzQU7bEQo8mstwML0FAsg6eswWEddn4IED7iTC/gvXUxP89Qo0qrDyVMJ1nYltYsRefZfCi2984owi3C1RQC7uQJfp6AX+Mo14eSjkRD8fw2NOW586YCANqS8aJjfqFWnDfcKgDB/QgAJGLOfQyOFA1F+jJuBbBKCPEtz6H7ACyUvxcCuNeTfomM3JoDx/n+gjSL9RLRKbLMywvyqLIuBvCE9LM8CuA8IhovnfHnybSqkhPwOOWDRx82p3w5y3nYRpHhNRSJdgVe52+QKck7OiwaVXs78ZBO+cKz1KXDCM5Au3xIkxagf56lmKtM/XgpYcOFnVx4DSX/2zYPwmhztwgMmx8k0IcSwilfuoZSfE4mhA9FvfKmK9j8bgBwYNDZZ8dq9g54p4KwrfAANK6GUsny9acBuAzA60T0ikz7IoAbAdxNRFcA2ATgIwAghFhBRHcDWAknQuwqIYRqjU8C+BmANgAPy3+AI7DuIKK1cDSTS2RZPUT0VQAvyvOuF0LowzYiQr2gSmUOUu0Lf+vSwnw0hddwTV5lqOK6cnQaSqn27sIRsTeXTaCYOltr9JUIbuf8ecF18WkoOhNmCeYqewi0tYgiIV7Ofig2k4pJg7H55/Lmv2CBojsaph3934591K27z8FMvnM1XUvnN/QiQphTDww5AkXXVjZfUKFZXL/xm/d3sGCMYjJrlJQtUIQQz8K8HffZhjw3ALhBk74MwLGa9EFIgaQ5tgjAorD1rRT1kJWqafOh2KLASp285Z+lrz8//EJyznm6b0tdMmh5laBlMvwaSnAtrFspW+z5hXUpOs/Qfrq0Il+QENZNpUx10l0jzLMpN8rLphX6fEUGzdOmjeWXT9EfD2qr5ZucMNsgOeHviM3n6c5XDAx7BYo+XykLrJrKUJ247rhXA3Sep98s5n0+2/YPYua4tqIywrzf6tNNZyobXEYNz5QPiauhJMw+FH+4X+kCJ6i8sGaeMKgPO6aN8nLKfnTFDmN+vw/F//WHWVxPYds3JpQzPOAapZi8TKsSAJX5UGwjey/F81BCaiiWDsjnczIJFMs18wIn+Jnpju+S854CVzWIQkNJ2zUU6+oIlkGI99r6uTBeDcR8fQDoH8oUn1BUB+0prqZV6byjqGGBEhL1jJNSRTVNbEy5TntdGaVqKPnf3rMrdcorgkxe3b1DRvusL8orYDKe7R5tEwptkVOAbXHI4GsFdbR+f5fxEtrzvagOOkwZRT4UmcdipbHeZ1BUXv5awfW0mUKVyTDQ8R7SjBoqYEWnoXgFiqEMW1i196g1AAHFQsdv1tb1Efnf/cOG76uEfqLRTF4sUEKiPiTllDd9uKkADcbbaZWuoXhHPsGdFwBM7moxlhs08g0TQeXtlIpeaItW4MWuodgFZ3infGkaY6kBFGbh6PwN55Q3C+cg/D6U4LpZNRRLJ2rqj03rsvmuIcJpBmG+DV1nPeTZLdFUhG0l4zAhu0GRXDYfirfMvmG9hmIzYQL5wWAUWxZECe8pHxL1XBMB81ByORHotC91FdGc4cW1aShHTevCpE6zQFFl6Uxe/s5HP/rxOrAL26EUU5FZAymui+kawasNB5sfvBFRhcIzzGxl3/kWv1YoH4pJQ7Hks7VHKfNUKl3Ly2atygkgrrmhkic2as7xmn9MZShBZLpCmM68cCDi3RnU+wy1k2k9eQdMGkoI7bhRnfKsoYREjbxSAfNQsiJ4norvRSxxz4cwkTwq1RrJIorLLDwGmM0DgVFeEWgo+Sgvfb6wpkNv1XSdg3eWcZCGEmZlgqrMQwm7BbBFm/KbKIPnX5juNS8wgp9ZYZ2LfVPBWqn3WkHozhnO5NCSiIEo4D4s351X8zE9sqA1w3ymu4DIQcBs8go3uZJ9KE2Nek9cDcUgMAKXZvF9NCEcjyanvEVDiZHF9OAuvlh8LBeijkEdlPeqNi2snFnXhRepxCk/XBSRU1yHoHrcviS/OJ8tIs0oHAPaUh3TaZKmumo1lBIm1BpNXtlwGkrhO1V4T0YzZ4il531mMU1nPZTJIZWIIUZkfC9s312Y5x40T8T6LDxJT63WLwUVRqCoV6LRTF4sUEKwuacfj69yop5sa3m1SB+KdnHIrP1F8WKK9rCZBGxOXNUh66Jpwtjb/RMbC0bVXsFnc8rbfCghnPLB81CCtYxhz7yFcpzy976SX/7fZoYMY48vrIM6FtfMVTDV1bZmmd3kZbiGCBiE+DpR/wnz/u0RAHBNsKHawWiuCj5nKJNDSyKOOFFAyK9SxfTHM1n7c88EtHfQMcDfVr95eUvRceec/G9r2HCDaSjsQwnB+7/7DHpliJ8rUAzzUAJNXiWq9T4NxeJ49Z5DoEAnsDpkDTE1fdihw4bNdQgq327yyv8O9qHo66UI8qGEGSWGsdmrZJuz26mDvy1VO9sGCNZVrn0Cpfi4MLxnvrook5bmRoJMoOqQO9AyPlP7IMF2jtfkZZvRbyKUhhIwMMxY3glv2umHT9KXX8I3xCavJqTXEy/ubrClHQkiMMorTCSL6ZySTV5BAkWNNjXn+CfBlW7y8mKdKW/RUGyhuEHneI/FyGzyajE8r4yh7b14J5VZJzaaNmvyRczp62DrVOwmr+BnmhP637oybMEmpqVfWpLm7wIoR0Mpvs5w1jF5xWNk1oTcSYn2awT5k5QvtXC9r7BRXodMbPeFOXsJ45tjp/wIQflQtFuQ5nJ5gWLxoYT5sEyjZGsUjHV5CVmmRUMJ45QvXhzSW5ZlNGiJ4rItAwKEM3kl4jHttYbSObSnnJnMhR2UzdENFCxGGLxElvG4L2KuSEMJ11lYw6M9aVv3DRYd913X5iPRHPZpKAaB0pqIy+KDB1LJuFkYBM1/cq6dRSoufSiGpkvLMkzvdhhTp2OJ0Edz2qYGqPM7WxLmeSihTF7maNN6wgKlRNTSK4UdVDqbw2A6P+ItZ2mWwnTTxEbbaNjqlJeHdCPBcE55szM7zFpIbl7TEuJu5xV8n0HnOOc5f1PxmNZEOZzNoT3lWH1NGgpRuNGsWTgG30soDUWbM082J9wxhD7QIv97536NQLH4DWyz2L1tZ9qfI6yGkorHrO8FoBdcw65TPuCZyQYyvdv+FYtNdc3lV8woOMe2Xp9K6mhJGCcO24JJgPyY0aTl1AsWKCVi2rHx+0+sBQD8cd1uAHrbZhiB4n3RTbZUszCSdYyF86HYZo+HChsu+LDDmE/cvBaTmil/2AAFdX+pRMwYZprXUPQdQzIWM2oX6RIik8IMIIqd8jm3nkFkc8HBIt669Q4WT6bzrdKrEV82DcR7D6ZnqjSUIJMakdQmDQ3u14zNJq9YgMlLCTyT0DIN6Lz1zAnzBOd0SB9KRypunIcSKspL/g3aPbIesEApEdOLtLmn3/d/3QtfuobiFSj2vKpOtnko+YX+issJ40PJ+Mw0/jJKMXkZR6IRRXmpzjUVj+lt7pkc2lukhlI0QTNvgjGZabwjTFsItMnWnQk0eeWPvRq4j3kOLfGA6EKLQPFqSbrmzJYiUAzCL5kwrzABOO9RnAiJmDnk1xatNiQtBHEi43uh6md8HhZTp/edAoothGF9KJ2tSfQbZsr7BIrhPryrFO8fTGvPqQcsUEpEhXAWvo9q6Wx1vFyTl1crKFVDyQlnlEewmEmEuZww62QFjRR9Ji+LC0DlLZR/tnko4X0ozt9UImaMCmpP6n0o6vkl4jHjqHpPf/5Dti0jE6YDM5m8AOCdvQP6SsDxa5g0Z28dAGAoUzwqDppX4dQjX3edSSuMQBmUo/Eg/0g8RojHzMLA24Y6TcivoWiLcOtnXFbIt+Oivp5AfhuLwrp682t3AZXnd7bEMZjWt5V/Qq72FN91ew4M60+qAyxQSkQtDln4YaiXo6vVGfHqXnifOhzC8ejTULznBH2URNYwU/WimybBKaFo6gS9JonCDqQUDcUVKIXl26K8POXu6Td/TN49bLQj2qzZ5OXVUEz3MWNsq+da+jrkNRTTM8vXa6BgxOrVEE3mEQByyR+zycvXEWs6fH8nWpzf2/Hp3mv1PrUmY77j3vf3hIPHOfUL0MwTUqAY5z950o0+lLj0oRiukzd56d9tr08iaJvhlEHjspnM1C10pBIYzuYMUXd2k1c6KzC2LQmgsfaVZ4ESAu+8sjHyIRa+sMr8oQSK7qMYCpj3kE/X22DDqME54VGFA/pyVa5pEpxS580TG4NMXvaPQaE62cKZ4O4yHobsKn1MawKbeswj95z78ceLOjIhBIYzObRJgWJyyidiMWNTHjal090G1uZsNnVgXkHT0+c3XQTN7/Dinf+keze8na+u8/H7UIrxCjNdR65e2fZUQju355/OPQJHTO1yyjc0ZibnbDYVDzB5eQWf1ocinfLOxMZggWIS8APDWfcb1mkQyjSab++C+8jmJzcHaSgd0tSqc6qHiTDM5HLolGXYfGy1hAVKCLwdntsBeR70P//6VdcZP2Oss2GOyWmoMK3l5Td55dNts6EB5+WLkZzYGCBR8stk6M1yrUnzNseAv1NKB2ooxio4eQ0mL5XPGOUl723KmFbsGzDbj1X7peJU1IGo/3ek9AMA1TaJAB/KcCaHCR0pY11zOeG2x3AIJ3BP35DvmD/6ytyY2ZwIDFf3vos6k5dvbSrNQ+tPO5rTmNaEtvNS+duS8QKB4vxOysgrVVcdOeFoKIlAk5fZPAjkl16hgKVX1DdoivIaSGcxvt15poMBbaUGEoXbAaRz5rlNQP45drTE3esVn5P/bRyoZIWrXQ9nGyfSiwVKCLwd6/xDJgDwf3j3vJRfQuGbHzkegHkEpTC98D19eROOKUwyaNe8UCYv+TFqzR854X4spo+ufyiDMXIUly544b3/C2vyKiSsD2V8exL7AwSK0gpak8UaiupYTBqKGyEW4EMZ8kSJ2WaomzWUfPruPr/5zr9Zk7kth7M5tCbN9RjO5M13umd+YCh4p0M1X2Jce0qrKam2akvFkcnld7rMa3nkmlEDnfIxQixGobR3rclLTlRNxs1msyHZgZuODwxnMV4OEoY0nb26N/XeFL7Dmaxwn0WQP6uzxbF0DA5rtCCD79R3nVw+oGTI4IupByxQSuCD75ru7jOie1mScXLXLNKNoL7z2Br3t2nksXF3X/4cg4M8KPpJmbyC+nI3HFUbiZbvnEwfXd9wFp0tCa1vwntd25Lt+dGivjM39aHq+Ni2VLBAkQW0JuNF9Xzo9W0AgOfW7tLWIe+UN5tPgsKOAf8ztgnPZJx8gwmnzODoK2/ZSqvUnafaubMloTV5/cMdL7m/dZrtoCtQktrOS927Goi4jm/VhjFyIw/N+8Y4AiURowDt3eKUl2t5tSbjxjkeKp/JhDgwnMX4dtnZB9yr+kYKt+BNZ3OBc27U7SsNRWl/XsJNqhXolGUMNdBseRYoJTBvxhjESS29Unw8nc3PoC0ckW7Y1ef7v0mg/OPdr2rPCePQz+UEYhVqKJmsQEtSP/pS9A1l0NGSQDIeK7rPbMhO0FsPIQonzwn3fnSo2x/Xngw0eanyW5PFUV5qgp9KLZz7kHfKxwIcvFmMk+YRndPcNrMbyHdQU7pasa1gFnvYJW7S2Zw7z0P3XimzZFdrUitQtnsmO+peLaWhjG1LBs6vUqP2fCSV8zcRj+W/G8NtqLDheEy/qgFQYPLS3MdQJotUIoa2VNw44a9Q2BUykM5igjJ5acpQzyFvbip8/wVaEvZBhjK16t4b21yWdDaHbM7jlGcNpbl4zxGTAQCfOPMwyInyxk7GFCG1+4DfPh5mDR7dIoyphLmDy4r8Zj9BfbmyDevt4XYfSt9wFu0tCSRiVLR8RCkbJflDLItHZebZzk76+Pakb521ovNyORA5TvnCezlkYgcA4D8+dIysS4GG4vpQYsbw0XXdfe7abgc09bBNxPPey+FTOtHTN4ydvYO+PLMntgMIjuRJZ4VHQzH7UDpa9D6Q9x83DYCjwei0yv50XqBkPSYthStQkmrEnPWlJ+Nk9aFkcwLxOCEeCwqpD+GUj8fQFqChDBUIu8L8mZxwBwm6MlT7Kad6kckrF/wsXIES4JT3PmtdUwy4z8OpZyMtEMkCJQRf/4vj8PQX3ot4jNzoDtNDJCKk4rEiJ2yf7Hi/8/9OAODXRBSvyMlrV5w+B8k4+ZyCquMZ25bUdl6AjPIikqsNmztzNSoyTb5Uo12T3b9/KIOOVBxdrUm8vduvedkmhnkxjcS8M+X1q9s6gqKrNYnhTM7YWaezAslYDIkYFXUgqrNwneqFZjfVGRpmXT/x5k4AwONv7kR7Ko4+zTMJ2uvEey+As1ggAPQN+Z/59LFtINLb8xVB0WpA/l3taklonfLtqQRmjG2VIdLF5atwZuWsLixDXbO9UEOR72w8FnODBoL22EnEYojH9KsaAP4NqQoFbCabcxdnbUuaNRT13HUa44D7TkiTl1YLypsPgeJBmVdbNM1FI/I45TUail+gmE2QeQ2FnfJNxcxxbe6ItiURkyPz4g7k5o+eBMCxuxd2xqrDUSMTHR/+wXMAgOljW9GSiPtUWdUhTe5swd7+tLajdUxezvWDtANlGzZF7KgRlslMc2Aog/ZUAodO7sCBgpnXtolh/nM9Ha5hjxVdGcPZHJLxmNuB9Q/pP6hMNodE3LHLF37cg64ZSD9vSHVqLUl9B+f90DtbEloh79UqbXMrxrQ6nYP3vcrkBBJxQksiZtFQPCYvTYMNe+5VV87AcBZtqThipI8OVB351DGOf7Bwtr3Z5JXXUJQZyHQfWfXuxsiogavOt1OzDtawp61bA5Y1Uc9JN1ga8AQfAOFMXjqnvOtD0T0L+ayUNqfVUDxpOoEykPYLFNZQmhgiQmdrwu1IVccejxEuOHY6AMfuXviiqQl4B09o15br/YhOOXQiWpMxn4aihMDR08dgOJvDWzsOFJWhZhsnNRqSQgjhfihae7gnSsU0mtx1YAiTu1LoSCWKPoiwoa5AgYbiyed1Mv7rb17DEV96uOgayRi5o8Q+wxIWmZwMRY0XRw6pUd6YtqR2gKBG4Z0tCa1ASXo2Rt/ZO4S7XtxcJOSVOW6CITpK5QXygs07Cu8byqA9FXcGFwECZWA46w5UdJ2xEtydhrDfgbQjUMiwzP/q7b0AgGkyJL7Qb1XoF1Dvquq04zFyO1mTvX8wnUVLIm6MRAPyz3lCR6rIYa7y5E1exWUIIVxhWBidCOQ76o6WOFLxmLaM4QINRWfyagnwZw2lHT+P+sZ0Kw6HNXmNk8EDPA+lyUnFY3hWRgd1S9+I9+XpbElgf8EoTi2PMHtSOy4/9RC3AwGcTuDQLz7k/r81GfNpKEIIfPTWpQCA0w6fCABY3+0XKJt7+vHHdbsRIzLODAeAi37wnBuequuk0jmvQNHb43cdGMaUrlat8zNoN8dCfOadnHLQC6SzOUyR0XS/fmlLsfaQzSERj7mdqM7cpOqajMeQiBUHDyih2paMo6MlUaRpqc6koyWBbE4UddTxWPGno+YiKVZt2w/A6QBNZpwv/u51AMDsSY4G7F0TbmfvEKZ0tUoNRT/ifuLNHegdyng2dis+Z9gjHHXPvH84g/ZkQttOAHDXi5sB5DWU+1/d6s+fLtRgHIGjBFkiRu7cDNN97B9MY0xbItBcddsf3wYATO5qKTpHaR5tqbixjMF0Dmrped19qkFFWzLhDOh0PpSsX6AUDty8Gr55V8mY+/3r3l2vtqszZw0Umrw8z/StHb24+el11t1SqwULlDLY2TuEdd192LF/EJ/6+csAgI+dNts9fvCEdqzYus+X55cvbAIAtCTimNTZgt7BjDuy8M5jUee0eDSUP6zZ5R5THU/hgnB/9v1n3fWedBqS4rUt+XrpzvHOMdEJBNVJThsrBcpwYWefz2NbtO7N7fvd3965C0IAU2TnpPCO/tOyU1AftcmnpLYTiGvmNgxmsq4255is/B+u6vi61EjUt55VFgsXvQDAEfB/d9oc534LRu6f/sVyAMCMca3I5gT2FIQFezusmePafGUMZbLYN5DGlK4WdLYmsH9Af4+f/5Xji9s7MAwivVb57/evBOBoKNmcKOpMB9I5tKbiaE/FjXt0APmtG74nV9ZWKJOj0mBU+PNAWnXyCXfUbhIWvYMZjGlNBkZoKa18YkeqqLN/c5ujRc2d0omOlgR6B9NFnaoSdOPaU8iJYm1uXbfjD5w1vg2tybhW+L29u1+e41gaCqPNnJny5jlBQxknrFgNhgoHMoC/P9At5llk8vLU4Qu/fhU3PvwmVm7bX5SvFrBAqYCd+4ewbOMeAMBZR01x008+ZDzW7jzgPujlm/b4NjZSE6f2SjPYq1v8K8m2FGgoXtPRYZM7ARS/aHvlQoX7B9OOQNEukeF/wbt7/ZFnw5kc+oazmNTZghjpR5OPrdwBADj76CnaaJpMLj8xMmiOyP2vbsWWPfllU9ToUAm5qV2tvvO9Wkomm0Mi5vGhaDrBZ9Z04zcvb0FXa1I70W1gOIdWOWruaCl2qg+5a7MlZb3y+TfuzmsRV59zBK44QwoUgwD908OcrV5XFXzkXtORMl+o66hnM5HdxKkAABqXSURBVLmrBZM7W1xNuBBVRlsyjs5Uoui98D5Dk819cDiLtqQMtw0QKIdNcQYzap6GQpmi3jVrrO8+lZbe1ZrAzPFtiBG0plrAeVfGtCXRlkwYHNX5tLZU8Xu3UWp2h03uxIxxrUhnBXYVrDyg6qPqXzjLfasckM2Z1CHnshR/Q3ulsFRBFIWDssF0NnBuUt60F0NrMhYYpQjoBcqjb2wH4GhJ8Rj52mbNTqd9Te9LtWGBUgbKD+IN8VQRMIDzsuVE/gX98x/+0Zdfxbn3SIFy59JNvuOtyThak3kzhxoBPvCZ091ReaFJTdEnzR86H8r2gs2VtuwZ8H2Yam+FcR0pTOxswe6CVUyfWdONm55Yi5nj2hyTlzQteLWH4WwOk7qUwNR3sOu7D+Czdy33pamVe9VEsUINxTujeDiTQzJB+VGe5qO87FZHg1i9oxcJqbF569k/nEGbtPkXOtVXbduPlzc5A4VOtSKAR0B7NY1EjNAhO5A+Q3DAvBljAAAbCiLivO1T2Nnf/6oz8XJ8RwqTu1qwqze4gzj98EnoavXfRy4nsL47f021snKhH6M/7QRZdKQSRf4o1WFefc5cTOlqxSET23GmDKPP34fTHtPHtmJce9J9lp//1SsAHC2vsyWBqWNa8cBrfnOZYv9gBl2tCbSlYloN5Qu/fg2AYwnQmbTU4GVMW9L9FvcVvH9q1H6Q1C4KBxk79g+isyWBjhbH5KUTbHe/5Jj/2jTzUFTYsRqE6PZ12TeQdp91Z0uySGCo9v7Eew5FeyrualVebluyEYBjSm1JxHwaivqe9w+kMZjO4tJbnsdLctBbC1iglMH3Lj0RQD7MF8hrHUBe4Gzq6fd1Yl/78LEA8qGqj6/aWbSPCuBEkrUkYu6Hf2Ao7V4jLp3RuhcNcJx4qTihdzBd5KzzjqwVasLl5p5+LHpug1O/9hQmd7ZgV8EoR3XSyrR2YCjjmHI8H27fUAbj21NoT8WNkw7/5idL3egtFRmXN5NI88mYNl8e72hTfZSdFh+KYlxbEpmccDvb+1/disdW7cCkTuc5dBQIlAu/+wweeG2bewzwj0S99xuP5QXb1x5c6dbn3G89DcARVgtmT0BLIoa3Cya3qkCNSxcc7Dq01XW27XPa+ISDxmH62Fas39Xn0yj39A3jfd98Kl/n46ajqzXpey8+f/cruPC7zwAAfvDXJ7mdYKEAHhh2VkeY1JVygwQUqm1VFFpHKlEkOHfsH0J7Ko7OlgSmjWnFWzt6kcnmXIE5Va7KvG3fINZ392kmwzrPZkxrEu2phDaC8j7pt+kdzKA1WaxJ7RtIo6vVGbEr/0ThoOsuaXY+efZ4AH6BsmP/IN7c1uv67sa1pbB3wD+gGs7ksLknb1ZWafl2dOqkNCBdAIJXoBQOAID8YGXW+HZ0tSYCzcaTu1qQ8giU5Zv2uE7836/YgaO+8giWrN+Nz/5yubGMqGGBUgbHzRyLRIx8tuRpY/ImmoOlOryxp9/tVBeeegg+esohAPIC5b8fXY0z/utJdzKiIhWPoSOVwPb9g77IFNWBdrUmfNqD9+O6dMFBmNzVgqFMDkd8OR8d9cKGHlz64+cBAEuuPQs/+9ifAMibTLZ69tuY3NWCSV0tRSYxxTlHO+a9lVudEd8za7oBOKOjp1Z348BQBhM7U0UakULNCJ86pgXHznTMJOpDukeOAGeO9wsUb2fc05/G+PaU23E8tmqH79w7lrzt/r7yzEPd5XB+JSOxPvPL5dh1YBgTpUDpbEkYhZJa3sL74e/1LJnfP5z1bLrmtONnf7ncNT1ccfocxGKEo6ePcc2jCiWArjzzUMTksiNKoKSzOUzqbMHUMa044SCnA3z31x9z8z78xnZ3MPCRk2cBcN4L9a5ksjnc+0peGzh33lQ39H29py1PuP732HXAEQizxrfjnT0DvkHQNx5Z7bSDbGudeXDH/kFMHdMKIsJJh4zHWzt63XlXQF4YKQpXBFAd3pi2pGtq8ppn75ZBAaqt1Dmqnv/79Dr87I9vezpq569XuG7Y1Ve0gOu6XXnz26U/fh5L1u92NeOJnSnfN7Z/MI3rH1jh/r+zJYEY+c2Wu+WgZ3xHCql4zNcGCq9AGdOa8L1LQF6LmtSRwti2ZJGWv1223TUXHoXWZBxjWpPokedc+9vX3fMelEsLASi6RjVpaoFCRBcQ0WoiWktE19TqurEY4eNnHOr+/+t/fpxPKCj7/1f+7w1XKzjt8EnucWV/Vcyd0ul20oATmnzW0VOwqacfq7b1uiMtJVC27RvEfa9uxfu++RT29g9jh+y4v/yBo/HVi4511xvz8t+P5rcKnT62ze1kFz3raCXqGjPHteFds8ZicmeL64AEHKd4RyqOkw4eh+//taNV/MsFRwLIh4v+492OiWN9dx/eNXMc/vBWt29k/+TqnZh9zYMAgFMOnYDnrz3brcetz27A3v5hfPP3b7nnX3DMNBwjzUWqLkIIvLp5L9qScUyUeR96fbsbBLG5px9fudf58D/wrun44vuPxiTZHl97cJU70gWAiR1OeqGG4g0JnioHCoue24DF0n/k1VD+ZLazWOh3L3EmrJ72jSfwuJz0CACHTnY68ZMPGY/lm/bi4h/9EQPDWcy+5kG88Y7TeRwiNVonmMLpJJ98s9uNmpomR/g5kdcOvZ36DOnQ7/QIlELtMJWI4fApjv9NRQhu3zfodljTxrRi1vg2DGVyrv197c4DbjCJeqd0GsQza3a5PqDZE9uxpz/tCoGvXnSMe94P/8Z5b7xa+fruA27n154qnp+RyebwL795zT3/iKldmNjhhGGrQIX/lKZa5ZNTQSV/+9MX8Z3H3pK/X5BtTHjvkY7J7sUNPZ56+LXHCR0p32KdNz+1Dj9/3mmLy045BPEYYWJnftCVyebwV/+7BABw4kHj0F4geIcyWXz6Fy9jy54BNwpsxrg230AOAN54x3mPj5kxFgdPaMemAguGsoosmOO8d3OndOLljXsghMCR07qgo284a1w5IGqaVqAQURzADwBcCGAegEuJaF6trn/NhUe5v6cUdODe7TkvkpMVD/YIERWWq3hzey8OndyJufKDBxybOAC8/6ZncNPja5CKx1yhpbShDbv6cML1i90wzqOmjUEiHsPFJx/klrPs7R5s7ul3O46rz5kLAO7+FL9fuQPPrtnlCqV7PnkqWpNxzBjnLA1/51LHXrt25wH0DWfxkfkHufVXE8B+8ux6rNy6H2/K+Qp3f+JUnDtvKvYPZvDc2l1IZ53Z7H9cm49W+7cPHgMicstavaPX9/EcPb0LN192Mh74zOnoak3gNRm4cOSXHwGQH8ktPNXR+v76x05YtfJ9AI6pC4Br2gKAz931ivv7nHlTATgj+z39w8jlBH778hafA1619c+f34S/v30Z3t7V54747vv0ae4zUR944dSb2VIrUO/Iso178IMn/VFS6n1pS8Xx62Wb8cY7+9DTN+xqst4ByGk3PoHr7n3DN8r/s+Od+U8zxrXh9Xf2Ydu+AVfoTRvTilsuO9nXHj9/fiNyOeHO9geASxYchFlSK3zg1W3oHUy7C2cCwGkysGBMWxIbdvVBCIGtewfwv0+vw76BNFZIbVXd7w0PrQIAN+IJcOZXAcAjK7ajp28YS9btdvOptlMCWNn913rC43/7qT8FkNde39k74Bt9q0i5Lo9G9J3H1uDB17a5A7v/u+o0jGtPYVx7Eiu27sfAcNbtxAHgH891BkmTu1rQ0zeMwXQW2ZzAD59aB8CxTlwvheRkj0B5Zs0u7JIazeFTOpGIxXDH8xtd896SdbtdM6qKlpvYmfItCLpj/yC++fu3cPCEdhw8sR1zJnVgw64+V1vbN5DGP/zcWcjz6GnOQOs9R07GO3sHcOuzG3DvK1tx5NQu/OlhTjs/8y/vw41/cRwA4Pzv/MFoJo8S87TtxmcBgLVCiPUAQER3AbgIwMpaVeBzZ8/Fdx9fU6RxAMC3/up43/Iqc2S4r+LnV7zbnVsCOC/qp886HHvlJkvKPKHwOv+e+sJ78e6vP+6OQv9nsTMKmzbW6bTiMcKdH383/uYnS3HxzUuQjBMmdrTg/cdNw9XnHAHAGbFedsohuOP5jW49UomYq11dcfoc/PS5t/Gl372BL/3uDffaXoes6vCeX9+D99/k2Or/7PgZWDBngjsZ7m9/+mJR29y6cL7rqAaAz5x1OL73xFq3s7/3qtNwzAzHFEZEOPuoKfi/V7Zi275Btx1++reOye6T7z0cty3ZiH0DaVf7UVwmhc3kzmKN7cl/fq/7TA6d1IHBdM43FwgAjp05BgcVTER9r8dv8a5Z49zfXpOn4ruXnIDjD3LOmebZ3fH7HoGy+msXuL8PmdiO5Zv24oPfexYAcM7RjsCb1NmCq8+Z665WrZyyx88ai3s/fXo+v6zrqf/5hJv2w4+ehJMOdkxmCWmaW9fdhw/94FlX4z1qWhe6WpOYOc7Jf/0DK3H9A/nPaOkXz3aXTpl/yHjc/+pWzLnW31a3LpwPIC9YFX9+0kz394SOFI6bORa3L9mI2+U9KI6a1oULjp3mRkhdLsOyvaj7UIJPvXMAcPmph+D6i46V7ZXy5bvqF05o/wXHTMO86c57t7c/jfte3erTWH98+Xy3/soPetRXHsGVZzrWiHOOnoKbLj3RXTk5lYjh8Td3+t67fz7vCBARzj9mKu5cuglX3vESckLgqdXd7jmfeI9TXmeLE8Bw4vW/x+fOnosbH3E0rXfLOsyZ1ImhTA5fe3AV/mT2eF84sfKHqfp+7UFHgJ962ER85YPzkMk5Ky//1fyDcM1vX8fG3f3YsX/IJ2yrQdNqKABmAtjs+f8WmVYzPn/uEXj5K+di7tRiVfMvTpqFI6Y6GsdFJ8zwjdQA4PS5k1wzyUdPORgXHjsNY1qTPk3mxS+d4/6+68pT3N+tyTheve483HHFAjdt1vg2nxBSHx/ghKJu3z/ommcU//GhY/CZsw53/z93Sqc7Wh7XnsIXzj+y6L7UKBBwOgg10VJxsbTnz53S6XamXk48eBzOlh2l4vJTncidDbv6MGNsa5HqftEJzmN9QZoo/usv3+W2+bSxrfi8FJJeNvzn+3GUHMVN7mrB358xxw1rjcfIJ+AXzJlYlH/Zl8/BA585A8l4DMu+fE7R8UKzIhG5HR0AHDG10603AFx47HSfKQ0AHv+n9/jei0+ceajvuLeOV59zhOsrUZx1lL8dT587CYWoDlRxyqHOO/DGO/vx/PoeHH/QODz02TMA5LWDQiZ5BPKFciFJL39/xhycMdcZaIxrT+HRq8/ECQeNw/PXnu36lxRK8Hj58xNn4pGrz8SEjhRak3G3w1UcP2ssXvv389z/q9B5L9deeLT7OxGPFQ3gAMfkpoTBxQVtecqhE3DuvHx7egdOt/xhPQDne29P5cfgn35f/ttReT59lmMBuEoee+LNna4wOWhCG5Zce5Y7h+VEuS3ynv40/v3+lRhM5/AP7zkMX/6gY2g57xinPoue24BP3vkyHn9zJ86YOwlvfjU/CDlyapfbz0wd04JPve8wZ2UC+V7FYoS1N1yI/7vqNG2bRA3Z9qxoVIjoIwDOF0J8XP7/MgALhBCfKTjvSgBXAsDBBx988saNG4vKqhb9wxnc89IWfOC46a6934sQAgPprO8lLZVNu/uxu28IJ3oEiGLJut1Yv+sAFq/cgZnj2vBvfzavSLABwNqdvXj6rV34ixNn+qLVgPys9OFMDsl4fl8LL4PpLNZ39+GoaV0+cx/g7O/y3cfWIJMT+NT7DsPcKV1FQQiA40C9fclGfOj4GUVaAeDY3X+3/B1kcgKfOetwX0eVywms2XkAK7ftw7Z9g/j46Ye6I2ov6WwOD72+DefNm+aO8Lz30JJwVkBYMGeCtp36hjK48LvP4Jyjp+Lqc+cWOZsBxwx5/f0rcN2fHeNOQvWyY/8gnnhzJ2aMa3NXsfayde8ANvX0Y82OXnz0lEO07Z3LCbz4dg8WzJlQdHwwncWmnn48tmoHPnDc9CJNF3ACIC787jM4MJTBY//4Hp/2NJjOYuPufnztwZV4Z88AvnHxu4oGIt29Q7jnpS2YP3s8MlmBkw4Zp20vE/sH09i+bxBvbu/Fhu4+fOasw33vTS4n0NM/jEdXbMeMsW1475GTi+7znb0D6BvKYGA4i6Onjyl63u/sHcDmnn68s2cA2ZzAB941vWgdvf7hDHYfcCaEdrUkMba9+Hm+s3cAz63dhT19w7jyzEO1z+PVzXuxsacfFxwzzVePXyzdhLd29GLOpA6cddQUTOhIFdXhzqUbsbc/jXkzxqA1Eceph/kHN69t2Ytn1+7C9n2DOGxyJ/7y5FmuZulFbX9cLYjoJSFE8Wig8LwmFiinAvh3IcT58v/XAoAQ4j9NeebPny+WLVtWoxoyDMOMDMIKlGY2eb0IYC4RzSGiFIBLANxX5zoxDMOMWprWKS+EyBDRpwE8CiAOYJEQYoUlG8MwDFMlmlagAIAQ4iEAD1lPZBiGYapOM5u8GIZhmAaCBQrDMAwTCSxQGIZhmEhggcIwDMNEAgsUhmEYJhKadmJjORDRAICg0OKxAPYFHD8YwKaA42HKqPR4mHNs9YziGlHcR6X1rMV9NMIz57aM7niYc7gti48fKYTQL2fsRQgxav4B6LYcv6WS/CHLqOh4FPWM6BpR3Eelz6MW91H3Z85tWfP74LYsOA5gme16QohRZ/Laazl+f4X5w5RR6fEw51R6n2HOieI+Kq1nLe6jEZ45t2V0x8Ocw20Z7ngRo83ktUyEWI+mWvlrBdczOpqhjkBz1LMZ6gg0Rz1rXcew1xttGsotdc5fK7ie0dEMdQSao57NUEegOepZ6zqGut6o0lAYhmGY6jHaNBSGYRimSox6gUJEi4hoJxG94Uk7noiWENHrRHQ/EY2R6Ukiuk2mr1J7sMhjTxHRaiJ6Rf6bUqc6pojopzL9VSJ6ryfPyTJ9LRHdRLrdghqjntVsy4OI6En5/FYQ0edk+gQiWkxEa+Tf8Z4818o2W01E53vSq9aeEdezKu1Zah2JaKI8/wARfb+grIZpS0s9G6UtzyWil2SbvUREZ3nKqup3HkiYULCR/A/AmQBOAvCGJ+1FAO+Rv/8OwFfl778GcJf83Q7gbQCz5f+fAjC/Aep4FYCfyt9TALwEICb//wKAUwEQgIcBXNig9axmW04HcJL83QXgLQDzAPwXgGtk+jUAviF/zwPwKoAWAHMArAMQr3Z7RlzPqrRnGXXsAHA6gH8A8P2CshqpLYPq2ShteSKAGfL3sQDeqUVb2v6Neg1FCPEHAD0FyUcC+IP8vRjAX6rTAXQQUQJAG4BhAPsbrI7zADwu8+2EE144n4imAxgjhFginLfudgAfbrR6RlkfQx23CSFelr97AawCMBPARQBuk6fdhnzbXARnEDEkhNgAYC2ABdVuz6jqGVV9oqijEKJPCPEsgEFvOY3WlqZ6VpMy6rhcCLFVpq8A0EpELbX4zoMY9QLFwBsAPiR/fwTAQfL3PQD6AGyDM0v1m0IIbwf6U6kGf6UGaqapjq8CuIiIEkQ0B8DJ8thMAFs8+bfItGpTaj0VVW9LIpoNZ6S3FMBUIcQ2wPm44WhNgNNGmz3ZVLvVrD0rrKeiqu0Zso4mGq0tbTRaW/4lgOVCiCHU7zsHwALFxN8BuIqIXoKjfg7L9AUAsgBmwDEr/BMRHSqP/Y0Q4jgAZ8h/l9WpjovgvETLAHwHwB8BZOCov4XUIsSv1HoCNWhLIuoE8BsAVwshgrRMU7vVpD0jqCdQ5fYsoY7GIjRp9WzLIBqqLYnoGADfAPAJlaQ5rWahvCxQNAgh3hRCnCeEOBnAL+HYowHHh/KIECItzTTPQZpphBDvyL+9AH6B6psbtHUUQmSEEJ8XQpwghLgIwDgAa+B03rM8RcwCsLWw3AaoZ9XbkoiScD7aO4UQv5XJO6S5QJlgdsr0LfBrTqrdqt6eEdWzqu1ZYh1NNFpbGmmktiSiWQB+B+ByIYTqo+rynStYoGhQkRtEFAPwZQA3y0ObAJxFDh0ATgHwpjTbTJJ5kgA+CMfUU/M6ElG7rBuI6FwAGSHESqku9xLRKVJNvxzAvdWsYzn1rHZbynu/FcAqIcS3PIfuA7BQ/l6IfNvcB+ASaZ+eA2AugBeq3Z5R1bOa7VlGHbU0YFuaymmYtiSicQAeBHCtEOI5dXK9vnOXqL38zfYPzqh5G4A0HOl+BYDPwYmyeAvAjchPAO0E8Gs4TrCVAL4g8lEhLwF4TR77LmSETR3qOBvAajhOvccAHOIpZz6cD2AdgO+rPI1Uzxq05elwTACvAXhF/ns/gIlwggTWyL8TPHm+JNtsNTwRM9Vsz6jqWc32LLOOb8MJ3Dgg35F5DdqWRfVspLaEMzjr85z7CoAptfjOg/7xTHmGYRgmEtjkxTAMw0QCCxSGYRgmEligMAzDMJHAAoVhGIaJBBYoDMMwTCSwQGGYBoGI/oGILi/h/NnkWdmZYepNot4VYBjGmTQnhLjZfibDNC4sUBgmIuSifo/AWdTvRDiTOS8HcDSAb8GZGLsLwN8KIbYR0VNw1jA7DcB9RNQF4IAQ4ptEdAKcVQXa4UxQ+zshxB4iOhnOOmj9AJ6t3d0xjB02eTFMtBwJ4BYhxLvgbG1wFYDvAbhYOOuZLQJwg+f8cUKI9wgh/qegnNsB/Kss53UA18n0nwL4rBDi1GreBMOUA2soDBMtm0V+baWfA/ginA2QFsuVzuNwlqdR/KqwACIaC0fQPC2TbgPwa036HQAujP4WGKY8WKAwTLQUrmXUC2BFgEbRV0LZpCmfYRoGNnkxTLQcTERKeFwK4HkAk1UaESXlHhZGhBD7AOwhojNk0mUAnhZC7AWwj4hOl+l/E331GaZ8WENhmGhZBWAhEf0vnBVivwfgUQA3SZNVAs6GYiss5SwEcDMRtQNYD+BjMv1jABYRUb8sl2EaBl5tmGEiQkZ5PSCEOLbOVWGYusAmL4ZhGCYSWENhGIZhIoE1FIZhGCYSWKAwDMMwkcAChWEYhokEFigMwzBMJLBAYRiGYSKBBQrDMAwTCf8fRI1f+Jq/9D4AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'].plot()" ] @@ -215,9 +2240,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xt8nGWZ+P/PNTOZnJMmadJzaaHl0HIoUMvRY+XgEVxhLbpSlRXli7v6W3+rsu6+cPXLKq7CLrqyi1I5qAsIKl0R2VJEDkLbFGoPlNLQ0mPapk2a5jjH6/vHc086SSeZSTLJJJPr/XrNK5M7z/PknqfpXHPdR1FVjDHGmEz4cl0BY4wx44cFDWOMMRmzoGGMMSZjFjSMMcZkzIKGMcaYjFnQMMYYkzELGsYYYzJmQcMYY0zGLGgYY4zJWCDXFci2yZMn65w5c3JdDWOMGVfWr19/WFVr0x2Xd0Fjzpw51NfX57oaxhgzrojIrkyOs+YpY4wxGbOgYYwxJmMZBw0R8YvIqyLyW/d9tYisEpHt7mtV0rG3iEiDiGwTkSuSys8XkU3uZ3eJiLjyQhF52JWvEZE5Secsd79ju4gsz8aLNsYYMzSDyTS+CGxN+v5rwGpVnQ+sdt8jIguAZcBC4ErgRyLid+fcDdwIzHePK135DUCLqs4D7gRud9eqBm4FLgCWALcmBydjjDGjK6OgISIzgQ8AP0kqvgq43z2/H7g6qfwhVQ2p6k6gAVgiItOAClV9Sb1NPB7oc07iWo8CS10WcgWwSlWbVbUFWMXxQGOMMWaUZZpp/BvwFSCeVDZFVRsB3Nc6Vz4D2JN03F5XNsM971ve6xxVjQKtQM0A1+pFRG4UkXoRqW9qasrwJRljjBmstEFDRD4IHFLV9RleU1KU6QDlQz3neIHqPaq6WFUX19amHWZsjDFmiDLJNC4BPiwibwEPAe8RkZ8BB12TE+7rIXf8XmBW0vkzgf2ufGaK8l7niEgAqASaB7iWMRnpjsR4pH4Ptq2xMdmRNmio6i2qOlNV5+B1cD+jqn8FrAQSo5mWA4+75yuBZW5E1Fy8Du+1rgmrTUQudP0V1/c5J3Gta9zvUOAp4HIRqXId4Je7MmMysnrrIb7y6Ea2NrbluirG5IXhzAj/DvCIiNwA7AauBVDVLSLyCPAaEAVuVtWYO+cm4D6gGHjSPQDuBR4UkQa8DGOZu1aziHwLWOeO+6aqNg+jzmaCOdIRAqAzHM1xTYzJD4MKGqr6LPCse34EWNrPcbcBt6UorwfOTFHejQs6KX62AlgxmHoak9DcEQagOxJPc6QxJhM2I9zkteNBI5bmSGNMJixomLzWEzSiFjSMyQYLGiavtXRa85Qx2WRBw+S15o4IYM1TxmSLBQ2T11qsT8OYrLKgYfKWqtLsmqdCUWueMiYbLGiYvNURjhF2wcIyDWOyw4KGyVuJpimwoGFMtljQMHmruVfQsOYpY7LBgobJW4n+DLBMw5hssaBh8laiecrvE7qtI9yYrLCgYfJWonlqakWRZRrGZIkFDZO3mjvCBHzC5PJCCxrGZIkFDZO3WjrDVJUGKS7wEbKOcGOywoKGyVvNHWGqS4IUFfhtwUJjssSChslbLR0RqkoLKAz4rHnKmCxJGzREpEhE1orIn0Vki4j8syv/hojsE5EN7vH+pHNuEZEGEdkmIlcklZ8vIpvcz+5y277itoZ92JWvEZE5SecsF5Ht7rEcYzLU0hmmKpFpWPOUMVmRyc59IeA9qtouIgXACyKS2Kb1TlX9XvLBIrIAb7vWhcB04GkROdVt+Xo3cCPwMvA74Eq8LV9vAFpUdZ6ILANuBz4mItXArcBiQIH1IrJSVVuG97LNRNAdjVFc4KfA7yNkzVPGZEXaTEM97e7bAvfQAU65CnhIVUOquhNoAJaIyDSgQlVfUlUFHgCuTjrnfvf8UWCpy0KuAFaparMLFKvwAo0xaYWjcYIBH0UFPss0jMmSjPo0RMQvIhuAQ3hv4mvcj74gIhtFZIWIVLmyGcCepNP3urIZ7nnf8l7nqGoUaAVqBrhW3/rdKCL1IlLf1NSUyUsyE8DxoOG3Pg1jsiSjoKGqMVVdBMzEyxrOxGtqOgVYBDQC33eHS6pLDFA+1HOS63ePqi5W1cW1tbUDvhYzcYSjcQr8PgoL/ISicbwE1xgzHIMaPaWqR4FngStV9aALJnHgx8ASd9heYFbSaTOB/a58ZoryXueISACoBJoHuJYxaYVjx5unwPbUMCYbMhk9VSsik9zzYuC9wOuujyLhI8Bm93wlsMyNiJoLzAfWqmoj0CYiF7r+iuuBx5POSYyMugZ4xvV7PAVcLiJVrvnrcldmzIBUlUhMCfp9FAX8gC1aaEw2ZDJ6ahpwv4j48YLMI6r6WxF5UEQW4TUXvQV8DkBVt4jII8BrQBS42Y2cArgJuA8oxhs1lRiFdS/woIg04GUYy9y1mkXkW8A6d9w3VbV5GK/XTBDhmJdVJPo0wJZHNyYb0gYNVd0InJui/JMDnHMbcFuK8nrgzBTl3cC1/VxrBbAiXT2NSZbYsa8wqXnKMg1jhs9mhJu8lAgaBf6kTMPmahgzbBY0TF7q3TyVyDSsecqY4bKgYfJSJOoNr7WOcGOyy4KGyUvhmBcgggFvngZY0DAmGyxomLwU6tWnYc1TxmSLBQ2Tl3qPnvIyDVu00Jjhs6Bh8lIiaPSep2FBw5jhsqBh8lIk5jrCAz4KA9Y8ZUy2WNAweamnI9xvmYYx2WRBw+SlXpP7LNMwJmssaJi8FErq0wj4fQR8YjPCjckCCxomLyWPngJsIyZjssSChslLyR3hAEUFPttPw5gssKBh8lI4erwjHKAwYJmGMdlgQcPkpcSChQXJmYZ1hBszbBY0TF7qmdzntz4NY7Ipk+1ei0RkrYj8WUS2iMg/u/JqEVklItvd16qkc24RkQYR2SYiVySVny8im9zP7nLbvuK2hn3Yla8RkTlJ5yx3v2O7iCzHmAwcH3IrgNchnsg+jDFDl0mmEQLeo6rnAIuAK0XkQuBrwGpVnQ+sdt8jIgvwtmtdCFwJ/MhtFQtwN3Aj3r7h893PAW4AWlR1HnAncLu7VjVwK3ABsAS4NTk4GdOfcEwJBny4zyUE/L6eQGKMGbq0QUM97e7bAvdQ4Crgfld+P3C1e34V8JCqhlR1J9AALBGRaUCFqr6kqgo80OecxLUeBZa6LOQKYJWqNqtqC7CK44HGmH6Fo/GepinwmqkilmkYM2wZ9WmIiF9ENgCH8N7E1wBTVLURwH2tc4fPAPYknb7Xlc1wz/uW9zpHVaNAK1AzwLXygqpy1+rtbD/Yluuq5J1wLNYz3Ba8ZqpoXHNYI2PyQ0ZBQ1VjqroImImXNZw5wOGS6hIDlA/1nOO/UORGEakXkfqmpqYBqja2dEVi3LHqDX67sTHXVck7fTONAmueMiYrBjV6SlWPAs/iNREddE1OuK+H3GF7gVlJp80E9rvymSnKe50jIgGgEmge4Fp963WPqi5W1cW1tbWDeUk51RmOua/RHNck/4Sj8d6ZRsCap4zJhkxGT9WKyCT3vBh4L/A6sBJIjGZaDjzunq8ElrkRUXPxOrzXuiasNhG50PVXXN/nnMS1rgGecf0eTwGXi0iV6wC/3JXlhc6QFzQ6wjYUNNsiriM8ocAnPbPEjTFDF8jgmGnA/W4ElA94RFV/KyIvAY+IyA3AbuBaAFXdIiKPAK8BUeBmVU28K94E3AcUA0+6B8C9wIMi0oCXYSxz12oWkW8B69xx31TV5uG84LGkM+JlGJ0hyzSyLRSNU9CnecoyDWOGL23QUNWNwLkpyo8AS/s55zbgthTl9cAJ/SGq2o0LOil+tgJYka6e41Gieao9ZJlGtoVjqZqnLNMwZrhsRngOdVmfxogJR2MU2pBbY7LOgkYOJTIN69PIvr4d4QGfWNAwJgssaORQIsOwPo3sO6Ej3EZPGZMVFjRyqCfTsKCRdeFovGfdKUh0hCveoDxjzFBZ0Mgha54aOV5HuL/n+6ALIDYr3JjhsaCRQ12J5inrCM+6VDPCAWuiMmaYLGjkUCLTiMTUlrjIslDfjvBE0IhapmHMcFjQyKHOpGYp69fIrkgs3tMkBcebp2xPDWOGx4JGDiU3S3VYE1VWnbD2lMs0onELGsYMhwWNHErONDqtMzyrTpgRbs1TxmSFBY0c6rLmqRERiyuxuBL0Hx89VeACiDVPGTM8FjRyqDMcw+ea3Tts/amsSQwq6LvKLdjoKWOGy4JGDnWGo1SXFgLWp5FNiWyi7+Q+sKBhzHBZ0MihznCMyWVB99yCRrYkMo3CPsuIALbSrTHDZEEjhzrDMWrLvUzDlkfPnkSm0XePcLBMw5jhsqCRQ12RGLVlXtCwRQuzJ1WfRtCap4zJiky2e50lIn8Qka0iskVEvujKvyEi+0Rkg3u8P+mcW0SkQUS2icgVSeXni8gm97O73LavuK1hH3bla0RkTtI5y0Vku3ssJ490hqPUuOYpW38qe3qCRtLoqYAFDWOyIpPtXqPAl1X1FREpB9aLyCr3sztV9XvJB4vIArztWhcC04GnReRUt+Xr3cCNwMvA74Ar8bZ8vQFoUdV5IrIMuB34mIhUA7cCiwF1v3ulqrYM72XnXiyudEfilAQDlAT9lmlkUSRlR7ibEW7zNIwZlrSZhqo2quor7nkbsBWYMcApVwEPqWpIVXcCDcASEZkGVKjqS+qtT/0AcHXSOfe7548CS10WcgWwSlWbXaBYhRdoxr2uiJdZlBb6KQkGbPRUFoUGaJ6yGeHGDM+g+jRcs9G5wBpX9AUR2SgiK0SkypXNAPYknbbXlc1wz/uW9zpHVaNAK1AzwLXGvcRoqeJggLJCv83TyKKU8zSsecqYrMg4aIhIGfAY8CVVPYbX1HQKsAhoBL6fODTF6TpA+VDPSa7bjSJSLyL1TU1NA76OsSIxG7ykwMs0bMht9iRGTxX2WuXWjZ6y5iljhiWjoCEiBXgB4+eq+isAVT2oqjFVjQM/Bpa4w/cCs5JOnwnsd+UzU5T3OkdEAkAl0DzAtXpR1XtUdbGqLq6trc3kJeVcYq2pkqCfUss0siqRaRT4T2yesmVEjBmeTEZPCXAvsFVV70gqn5Z02EeAze75SmCZGxE1F5gPrFXVRqBNRC5017weeDzpnMTIqGuAZ1y/x1PA5SJS5Zq/Lndl414iaBQHrU8j2yIp52lY85Qx2ZDJ6KlLgE8Cm0Rkgyv7B+A6EVmE11z0FvA5AFXdIiKPAK/hjby62Y2cArgJuA8oxhs19aQrvxd4UEQa8DKMZe5azSLyLWCdO+6bqto8tJc6tiSao0qCAcoKA+xt6cxxjfLH8SG3J84Ij9qMcGOGJW3QUNUXSN238LsBzrkNuC1FeT1wZorybuDafq61AliRrp7jTXLzVEnQb0ujZ1HqjnDbhMmYbLAZ4TnS1atPI2BLo2dRKOrd26KCpKXRfdY8ZUw2WNDIkeOZhje5ryMcw+vGMcPVHTlx9JTPJ/h9YkHDmGGyoJEjx+dp+KkoLiAWV1tKJEsSmUZhwN+rvMAv1qdhzDBZ0MiR5D6N6lJv/amWjnAuq5Q3QtE4Pum9jAh4I6isT8OY4bGgkSOd4RgFfqHA76PGBY0jFjSyojsSozDgx62H2SPo91nzlDHDZEEjR7ojsZ6O2ioXNJo7QrmsUt4IReMUFpz4px3wi80IN2aYLGjkSCga72lz78k02i3TyIZQJE5Rn/4M8JqnLNMwZngsaORIKBrrGd1T3ZNpWNDIhu5oLGWmEfT7iMQt0zBmOCxo5EhyE0pZYYCg30dzpwWNbAhF4r2G2yYU+H1EopZpGDMcFjRyJJzUPCUiVJcGabbmqazwsrgUzVMBm6dhzHBZ0MgRr0/j+O2vLg1a81SWhKJxilJ1hPtsyK0xw2VBI0dCkVivtZFqyoI25DZLEkNu+7Iht8YMnwWNHOmbaVSVWKaRLX3vbUJBwGaEGzNcFjRyJHnILXjNUzYjPDu85ikbcmvMSLCgkSPhPsNCa0qDtIWiPesmmaHzmqf669OwTMOY4bCgkSMndISXJdafiuSqSnmjvxnhQRs9NWhHO8Mc67a/SXNcJtu9zhKRP4jIVhHZIiJfdOXVIrJKRLa7r1VJ59wiIg0isk1ErkgqP19ENrmf3eW2fcVtDfuwK18jInOSzlnufsd2EVlOnugbNI6vP2VLiQxXqJ+OcGueGrzP/2w9X/nlxlxXw4whmWQaUeDLqnoGcCFws4gsAL4GrFbV+cBq9z3uZ8uAhcCVwI9EJPE/+G7gRrx9w+e7nwPcALSo6jzgTuB2d61q4FbgAmAJcGtycBrP+r6xVZcWApZpZEN3P5lGgd9nHeGDoKps2tvKtoNtua6KGUPSBg1VbVTVV9zzNmArMAO4CrjfHXY/cLV7fhXwkKqGVHUn0AAsEZFpQIWqvqTebkMP9Dknca1HgaUuC7kCWKWqzaraAqzieKAZ11LN0wDLNIZLVXtNnExmS6MPTmNrNx3hGHtbOonZ8ivGGVSfhms2OhdYA0xR1UbwAgtQ5w6bAexJOm2vK5vhnvct73WOqkaBVqBmgGuNa6rab9CwYbfDE4qeuGtfQoHf+jQGY/uhdgAiMeXAse4c18aMFRkHDREpAx4DvqSqxwY6NEWZDlA+1HOS63ajiNSLSH1TU9MAVRsbIq6JpDBpWOik4gJ8YivdDlciaPQ75NbWnsrY9qRmqd1HOnNYEzOWZBQ0RKQAL2D8XFV/5YoPuiYn3NdDrnwvMCvp9JnAflc+M0V5r3NEJABUAs0DXKsXVb1HVRer6uLa2tpMXlJOJYbVBv2997CePqmY3c32n3M4QpHEVq/9LFhozSwZazjUjt/nfW7bY3+Xxslk9JQA9wJbVfWOpB+tBBKjmZYDjyeVL3MjoubidXivdU1YbSJyobvm9X3OSVzrGuAZ1+/xFHC5iFS5DvDLXdm41tOE0qez9pTaMnYcbs9FlfLGQM1TQdc85f1pmXQaDrWzaNYk/D6xDzOmRyCDYy4BPglsEpENruwfgO8Aj4jIDcBu4FoAVd0iIo8Ar+GNvLpZVRMz1m4C7gOKgSfdA7yg9KCINOBlGMvctZpF5FvAOnfcN1W1eYivdczo743t5NpS1r3VjKqesFWpyUwii+uveUoVYnEl4Lf7OxBVZfuhdj549jQOtXWzp8WChvGkDRqq+gKp+xYAlvZzzm3AbSnK64EzU5R344JOip+tAFakq+d4crwJpfcb28m1ZXSGYxw41s20yuJcVG3c6470n2kEXHNgJKakGFxlkjS1h2jtijCvroy3jnRYpmF62IzwHEgM++z7xnZKbSkAbx7qGPU65YtEplGYMtPwPvvYsNv0GtzIqXl1ZcyuLrE+DdPDgkYOhCL992kA1q8xDKEBMo3EUvRRCxppJUbxTa0oYlZ1CYfbw3SEojmulRkLLGjkQKJPI+jv/Wm4rryQssIAbx6yoDFU3Wn6NOD4kGfTv86wFyBKCgPMri4BsH4NA1jQyInjTSi9b7+IcHJtKTsOW/PUUA2UaRwPGpZppNMR8v5GS4P+nv61xlab4GcsaOTEQG9sJ08utUxjGNLNCAfr08hET6YRDFBVUgB4K94aY0EjB453hJ/YhDKvroz9rd28srtltKuVF9INuQXLNDLRGY5R4BeCAV/PEje2mKYBCxo50dM8leLT8LWLZ3FSTQl/9ZM1rN817qekjLqBhtwmgoatdJteZzhGSdAbkV9R5C1x02KZhsGCRk4kmqeCKd7YplQU8cvPXURJ0M+KF94a5ZqNfzbkNjs6QlFKg9499PmESSVBCxoGsKCREwO1uwPUVRSxaFYVb9g+BoOWCMhFKZcR8crCtmhhWp3hGCWFx+f+TiopsOYpA1jQyImBPg0nnDa1jJ2HO2zP8EHqjsbw+6Rn9neyxGi1kAWNtDrCUUqCSZuEWaZhHAsaORBOk2kAnDqlnGhc2WnDbwclFIn3e18TAw+6IxaI0+kMxXoFjUklQdvrxQAWNHIiFI3jEwj4+l807/SpFQBsO2BNVIPRd3OrZIkRVRY00usIRykNHm+eqi4tsEzDABY0ciIUjRMM+AZcyXbu5FICPrGgMUjdkVjK4bYAxUELGpnq26dRVRKkpTNiy8obCxq5EIrEUs7RSBYM+Di5ttQ6wwdpwEzDlSeG5Zr+JY+eAqgqDRKOxukMW8Cd6Cxo5MBAb2zJTp1SzjYLGoMSivYfkK15KnNdSfM0wOsIB5urYSxo5EQ4Gj9h3alUTp9azp7mLltddBBC0ThF/dzbRNDosqAxIFX1+jQKkzvCvaVEbNityWS71xUickhENieVfUNE9onIBvd4f9LPbhGRBhHZJiJXJJWfLyKb3M/uclu+4raFfdiVrxGROUnnLBeR7e6R2A523PMyjfS7AM2rKweO721g0useoOnP7xOCfp81T6URisaJK70zjVLLNIwnk0zjPuDKFOV3quoi9/gdgIgswNuqdaE750cikvgffDdwI96e4fOTrnkD0KKq84A7gdvdtaqBW4ELgCXArW6f8HEvFI31TDQbyLw6b3+NN5ssaGQqlCaLKyzwWfNUGonMtu+QW7CgYTIIGqr6HN6+3Zm4CnhIVUOquhNoAJaIyDSgQlVfUm/4xQPA1Unn3O+ePwosdVnIFcAqVW1W1RZgFamD17iT7o0t4aSaEgI+sUxjELx5Gv1ncUUFfpswmUais7vX5L6eRQstaEx0w+nT+IKIbHTNV4kMYAawJ+mYva5shnvet7zXOaoaBVqBmgGuNe4NNAEtWYHfx5zJpRY0BqE7GhswIBcX+OmyEUAD6nDLopcmDbmtLC5ABJo7rU9johtq0LgbOAVYBDQC33flqSYe6ADlQz2nFxG5UUTqRaS+qalpoHqPCQON8OlrXm0ZDdY8lbF0AbmowPo00klswJScafh9QmVxgWUaZmhBQ1UPqmpMVePAj/H6HMDLBmYlHToT2O/KZ6Yo73WOiASASrzmsP6ulao+96jqYlVdXFtbO5SXNKoyHXILXr/GriOdtshehtIF5KICf8+WsCa1zhSZBiQm+FnQmOiGFDRcH0XCR4DEyKqVwDI3ImouXof3WlVtBNpE5ELXX3E98HjSOYmRUdcAz7h+j6eAy0WkyjV/Xe7Kxr2wmxGeiXl1ZcTiyq4jtgZVJrx9IAYIGgG/dYSnkSrTAKgqsaVEDATSHSAi/w28C5gsInvxRjS9S0QW4TUXvQV8DkBVt4jII8BrQBS4WVUT/0NvwhuJVQw86R4A9wIPikgDXoaxzF2rWUS+Baxzx31TVfNiV6JMh9zC8RFUWw+0EYkpC6ZXjGTVxjVVpSuSJmgE/Rzrsnb5gfRkGsHebw/VpYXsbenMRZXMGJI2aKjqdSmK7x3g+NuA21KU1wNnpijvBq7t51orgBXp6jjehNJ01iY7ubYUgK8+upGuSIxffv4i3janeiSrN26FonFUj68xlUpRwMchyzQGlGr0FEBteZANe2wb4onOZoTnQKajp8CbYDV3cik+8dajemJj4wjXbvzqebMbYJ+SogJrnkonkWmU9OnTqC0rpLkjTCxuixZOZBY0ciAUy7x5CmDFp97Gqr97J+8+rZYnNzcSt/+0KfW82QX7T6Bt9FR6iT6N4j7Bd3J5IXGFIx2hXFTLjBEWNEaZqnprT2WYaYC3TPr0ScW8/6xpHDwW4lVrIkgpMf9ioOapYhs9lVZnOEpxgR9/n/1eassKATjcZp3hE5kFjVGW2Go009FTyd5zeh1Bv48nNh7IdrXyQqJ5qu8n5GRFNrkvrY5wrNdihQm15V7QaGq3TGMis6AxykIZbPXan/KiAs6dPck6I/vRXwdussICv+swtya+/nSGoimb+Cb3ZBoWNCYyCxqjLOQ6YQsH+DQ8kJNqStjT0pXNKuWNrojXpzHg6Ck3ai1kkyX71dHPXBfLNAxY0Bh1iU/DZSnS/0zMri6hqS1kTSwpHM80+u8IL7aNmNLqDEdPmA0O3gzx4gK/ZRoTnAWNUZZYDK64IO0UmZRmVZcA2CSrFLoyaJ6yjZjS6wj1P0GytrzQMo0JzoLGKEt8Gk7V0ZiJRNDYY0HjBIlAkEnzlA277V9XONbvYILa8kIOW9CY0CxojLLjG9wMMdOo8oLG7iMWNPrKpCO8KGDNU+kc645QUVyQ8meTy4I0WfPUhGZBY5QNN9OYXBakuMBvneEpJO5t0UCr3AYtaKRztDPCpH6Chpdp2DyNicyCxihLZBp9F4PLlIgwq7qY3c2WafTV5Sal+XyptmLxJAKK9WmkForG6IrEmFTSX6bhLSUSiVnz3kRlQWOUZdKEks6sqhL2WNA4Qbpl0SFpyK31aaTU6lYArnR7gveVGHbb7DZjOtZtKwZPNBY0Rtnx5qmhZRrgdYbvbemyCWp9dIVjPaOj+lNkQ24H1Oq2c63sr3nKTfBraguxo6mdxd96mofW7h61+pncs6AxyjrDUUSGNiM8YVZ1Ce2hKC22X3MvmWQaPfM0bP2plI66TGOgPg2AHYc7eHjdHsKxON/73zd6ml1N/hv6x10zJB2hGKXBAN4GhkMz2w273d3cSXVp6maEiagzzQZMkDRPI2zNU6kcdR9E+uvTWDi9ktnVJfzHMw0c6Qgxv66M7Yfa+afHN3NKbRnXLp5JXXnRaFbZjLK0H3dFZIWIHBKRzUll1SKySkS2u69VST+7RUQaRGSbiFyRVH6+iGxyP7vLbfuK2xr2YVe+RkTmJJ2z3P2O7SKS2BJ2XOsMR4fVnwFw6hS3m1/jsWxUKW90h2MDztGA5Hkalmmk0tqTaaT+MBIM+PjKlaex7WAbh9vD/MP7z+B9Z07lV6/s41+f2sbKDftHs7omBzJpI7kPuLJP2deA1ao6H1jtvkdEFuBt17rQnfMjEUn8L74buBFv3/D5Sde8AWhR1XnAncDt7lrVeFvLXgAsAW5NDk7jlbeC6PASvNnVJVQUBdi4tzVLtcoPnZHUC+0lK7LmqQEddXuA99enAfCBs6Zx7uxJzJhUzDtOreXOjy3imS+/k6DfZ7PFJ4C0QUNVn8PbuzvZVcD97vn9wNVJ5Q+pakhVdwIXRT0+AAAgAElEQVQNwBIRmQZUqOpL6vXePtDnnMS1HgWWuizkCmCVqjaraguwihOD17jjrSA6vExDRDh75iQ27TuapVrlh84MMo1EX5LNCE+ttSuCCJQX9R98RYSffuptPHrTRfh9QlGBn5Nry6gpC9Jsczjy3lB7Y6eoaiOA+1rnymcAe5KO2+vKZrjnfct7naOqUaAVqBngWuNaZzg25Dkayc6aWcm2A23WzJKkKxwbcKtX8N7wvN377L6lcrQzQmVxwYBzXQAmlQSZVlncq6y6NMiRDgsa+S7bo6dS/aXpAOVDPaf3LxW5UUTqRaS+qakpo4rmSmc4SskQZ4MnO3tGJZGYsu1AWxZqlR8yGT0Ftk/4QFq7+p8Nnk5NWSFHrHkq7w01aBx0TU64r4dc+V5gVtJxM4H9rnxmivJe54hIAKjEaw7r71onUNV7VHWxqi6ura0d4ksaHf3tVTBYZ82sBGDTPuvXSOgKxyjOIIsrCljQ6M/RrsiA/RkDmWyZxoQw1KCxEkiMZloOPJ5UvsyNiJqL1+G91jVhtYnIha6/4vo+5ySudQ3wjOv3eAq4XESqXAf45a5sXOtvV7TBmjGpmKqSAjZZZzgA0ViccCyeUUAuDvqtT6MfrZ3hfmeDp1NdGuSI9WnkvbTvXiLy38C7gMkishdvRNN3gEdE5AZgN3AtgKpuEZFHgNeAKHCzqiY+0t2ENxKrGHjSPQDuBR4UkQa8DGOZu1aziHwLWOeO+6aq9u2QH3c6wjFKs5BpiAgLp1fy+gEbdgveHA0YeH/whMKAz9ae6sfRrggn1ZQO6dyaskK6IjE3rNymgOWrtP+yqnpdPz9a2s/xtwG3pSivB85MUd6NCzopfrYCWJGujuOJ16eRnf9Q0ycV8ew269OA4xswpRs9BdanMZDWrki/E/vSqXETTY+0hymptqCRr2wZkVEUjsaJxDQrmQbA1MpimtpDtuIome3al1BU4LMFC1OIx5XWYfRp1JS5oGH9GnnNgsYo6spgD+vBmFpRhCq2KQ6DWz24vKjAVmdNoa07iurAE/sGUlOWWAHX/h7zmQWNUZTYH3yoGzD1NbXS+0964Fh3Vq43nnVF3N7rGQTkmtJgz9Le5rijXd49mTTEjvBE85Rt0pTfLGiMos5w5m9smZha4U2uOtBqQWMwmUZVaZCWzrAtLd9Ha5oVbtPpaZ6yoJHXLGiMoo6Q20sja30a3mqiFjSOB41MRk9VlwSJxJR2W867l8QKt5VD7AgvCQYoKvBZ81Ses6AxihLNU9nq06gqKSAY8FnzFIPrCK9yzSgtHdavkSzdXhqZqCkttEwjz1nQGEWdiUwjS30aIsLUiiLLNEhunkofkKtLvTfF5k57c0t2yH34SGy0NBQ1ZUEOW39RXrOgMYqynWmA10RlQQPa3GioTAJylevotWaU3va2dFFWGBjy6ClIDDKw+5rPLGiMoq5wdjMN8IbdWvMU7D/aRXlhgPKi9G94NaWJoaHWPJVsb0sXMyYVD2tXSW/RQss08pkFjVHUkeV5GgDTKr2gMdFHAu072sWMquL0BwJVrnmqxZpRehnMPexPTZm3/tRE/3vMZxY0RlFnKNE8lb1MY0pFEeFonJbOif2pOfEpORNlhQEK/GJ9Gn3sbelk5jCDRl15EeFYvGcklsk/FjRGUUc4RtDvo8CfvdueGHbb2NqVtWuOR/sH8SlZRKgqCVqmkeRYd4S27mjGgbc/da4T/ZCtUpC3LGiMomxtwJRsmgsae1smbtBo645wbJBveLbLXG/73N/PzKqSYV1nSoX393jQ+tnylgWNUdTaFRlw7+WhOH1qBQGfsGHPxN0vfN9R7w1v+iCChmUavSU+dAy3T8MyjfxnQWMUNR7tPmFf5eEqDvpZOKOS9W+1ZPW648m+IbzhVZcGrU8jyb6WToDh92lUJIKGZRr5yoLGKNp3NPPO2sFYfFIVf957lFB0Yu4Rkcg0Zg6yecoyjeP2tnRRVODrWXRwqEqCAcoLAxw6ZplGvrKgMUpiceXgse6ePohsetucKkLROJv3Tcxd/Pa1dBH0+5hclvlM5qrSIEe7IsTiNjQUvMA7fZhzNBLqKgot08hjwwoaIvKWiGwSkQ0iUu/KqkVklYhsd1+rko6/RUQaRGSbiFyRVH6+u06DiNzl9hHH7TX+sCtfIyJzhlPfXGpqCxGN66Da3TN1/knVAKzfNe53wx0S7w2vCJ8v8ze86pICVOGoNVEBXqYx3E7whLryIg5appG3spFpvFtVF6nqYvf914DVqjofWO2+R0QW4O3/vRC4EviRiCSGEt0N3AjMd48rXfkNQIuqzgPuBG7PQn1zItGEMhLNU7XlhcypKaF+gvZrJD4lD0bPooUWNFBVdjcPf45GwhTLNPLaSDRPXQXc757fD1ydVP6QqoZUdSfQACwRkWlAhaq+pN400gf6nJO41qPAUslG/pwDiXkU0yZlv3kKYNGsSWza1zoi1x7r9g1iYl9CdWli/SmbhHbgWDetXRFOn1qelevVVXiZhs0Kz0/DDRoK/K+IrBeRG13ZFFVtBHBf61z5DGBP0rl7XdkM97xvea9zVDUKtAI1fSshIjeKSL2I1Dc1NQ3zJY2M/UMYFjoY86eU09jaTccE2yPiUFs3h9pCnFJXNqjzEiu52rpdsLXR6ws7Y1pFVq5XV15IOBrnWNfE+lucKIYbNC5R1fOA9wE3i8g7Bjg2VYagA5QPdE7vAtV7VHWxqi6ura1NV+ec2H+0m/LCABUZLKg3FKfUlgKwo6ljRK4/Vq3b6TXJXTC3elDnzakpxSfQcLBtJKo1rmxt9O7BaVnMNMCG3earYQUNVd3vvh4Cfg0sAQ66Jifc10Pu8L3ArKTTZwL7XfnMFOW9zhGRAFAJjMve3v1DaHcfjFNqvU/abza1j9jvGIvW7DxCSdDPmTMqB3VeUYGfk2pK2X5oYt2vVLY2HmNmVXHWPtBMcVnceOgMV1UeW7+Xm3/+Crf8aiNH2sd+nXNtyEFDREpFpDzxHLgc2AysBJa7w5YDj7vnK4FlbkTUXLwO77WuCatNRC50/RXX9zknca1rgGd0nDaU7m/tGrH+DICTakrx+2TiBY0dzZx/UtWQ1vOaX1fGG5Zp8PqBNk6fmp2mKRhfmcaz25r48i//zPpdLTy2fh+X3/kcOybY/6HBGk6mMQV4QUT+DKwFnlDV3wPfAS4Tke3AZe57VHUL8AjwGvB74GZVTcxGuwn4CV7n+JvAk678XqBGRBqAv8ONxBqP9h/tHtFMIxjwcVJ1CQ0T5JNzdyTGkfYQ2w62ceHJJ3RzZWT+lDLeOtI5YSdFgncfdzS1s2BadpqmwBs9BdA4DjYH+/Wr+5hUUsBzX3k3v775Yo50hHlqy8FcV2tMG/JCSKq6AzgnRfkRYGk/59wG3JaivB44M0V5N3DtUOs4VnRHYjR3hEdkuG2yk2vLJkSm0RmOcsl3numZiDbY/oyEU6eUE4srOw93ZPWT9niy/WA7cc1eJzh4s8KnVhTx5hj/ANMRirLqtYN85LwZBAM+Fk6vZFplkWWfadiM8FHw1hGvc3qkg8a8ujLeOtxJNBYf0d+Ta2t3NtPSGaG00M/0yiLOnjlpSNeZX+d9un7j4Nh+cxtJW/Z7w7RPz2LQ8K5XztYDY/vN9+mtB+mKxLh60YyeslOnlLNtjNc71yxojIJ1btLdebOr0hw5PKfUlhKOxfN+mfQXGw4T9Pv43y+9kz/dspRgYGh/xifXeiOotk+QT5aqyk+e30HDoeOv96ktB5heWcRJ1dmZDZ5w+tQKGg61ERnDH2B+9co+plcWsfik4/8vT5taTkNTe95/8BoOCxqjYO3OZqZWFDGremQzjcRchXwfEfRiwxHOO2kSxcPcAbGowM+cmtIJ0xzxWuMx/u8TW/ncg+t7mkyf336YD50zfVBLsGTijGnlRGI6ZptLdx3p4LntTVy7eFav137qlHLC0Ti7mjtzWLuxzYLGCFNV1u48wtvmVmdlMbiBnDG1gsKAjz+9eXhEf08uHWkP8VrjMS6dNzkr1zttajmb9rZOiNnLj2/Y70bYdfAvv9vKE5saicaVDy+anvXflegjeb1xbAbkX6zZjU+E65bM7lV+2hTXZDnKTVTPvH6Qrz66cUxnZgkWNEbY7uZODh4LsWSInbWDURz0c+m8yTy99WDevgn+6c0jAFySpaDx7tPq2N/anfcrBMfiysoN+3nXqbV86uI5PPDSLv555Rbm1ZWxIMv9GQBzJ5cS9PvYemDs3dfuSIxH6vdw2RlTerZLTphXV4YIbBvF7LMjFOUrj27i4fo93LnqjVH7vUNlQWOErdnpzUUc6gifwXrvginsae7K287dp7cepLK4gLMGOZmvP5ctmILfJzy5uTEr1xur1u5s5sCxbq46dwb/9MEFfO/aczippoTPXDJ3RDLgAr+PeXVlYyLTONweoq37+Bpjv9vUSEtnhE9edNIJxxYH/ZxUXTLsJsudhzvYk2ET14+f38Hh9hAXnlzN3X98k5d3HBnW7x5pFjRG2NqdzVSVFDCvdnBrIw3V0tO9pb6e3pp/Y81buyL8fvMBPnTONAJDmMyXSlVpkAtPrub3mw/kbXbW1h3htt+9RnlhgMvO8ILkNefPZPWX38XHL5id/gJDdPq08p51rXIlGotz9X+8yGcfqO/5933w5V2cXFvKxaeknt9z6pRyNu87NuS/hz3NnVz1wxf40A9foOFQG6pKvJ99Ww61dXPPczt4/1lTWfGptzGzqphbH98ypjviLWiMIFXlTw2HuWBuTdY7GvtTV1HEOTMr+d8tB0bl942mJzY2EorGufb8WekPHoT3nTmNHYc78jI7C0VjfPaBel5vbOOuj5877MEDg7H4pGoOtYV46c3cfXJ+astB9rZ08fKOZv74RhOb97Xy6u6j/NUFJ/WbYV2+cCq7mzv57cbU2efBY909fQ/haJzN+1p5seEwqkooGuPmX7yCAgGfcO1/vsSib67inG/+Lzc+UM/z25t6BaN/e3o74Wicr1xxOiXBAP/wvjPYdrCNh9btSfm7xwILGiPorSOd7G/t5pL52Wl/z9RVi2bw572tebMpUzgaZ/2uZn6xdhenTinj7JnZaZpKuGLhVIIBHz95fkdWr5trqso//nozL+9o5vt/eQ7vPq0u/UlZ9BfnzaC2vJB/X527dvr7/rSTWdXFzKou5rYntvLP/7OFogIfHz1/Zr/nfOTcGZw+tZzvPvX6CasFvHW4g3f+6x/40kMb6AxHef9dz/PBH7zAJ36yhsc37OeOVW+wcW8r37v2HO779BIWTK/gyoVT+cBZ09iw5yifvHctn7lvHd2RGA2H2nh43R7+6sKTmDPZW3D0yjOnsmRuNbf//nUeWben3wwllyxojKAXG7xRTJf0kwaPlGVLZlFdGuSHzzSM6u8dCd2RGB//8ct89O6X2LzvGB972+yst8HXlhey/KKTeOyVvXkzsUtV+cEzDfxy/V7+dul8rkqawDZaigr83PTOU3h5RzOrc9BcunlfK+veamH5RXP4+ytOZ/uhdrY2tvH/vfdUKov7X5zR7xO+/oEz2NPcxYMv7eopV1W+/ptNdEfiPLGpkc8+UE/DoXZu+8iZnDd7Ev/4m838+LkdXLdkFlcsnMqZMyr5+V9fyO3XnM13Pno2z3/13Xz9/Wfwh21NfOa+dXzy3rWUFPj5m/fM6/kdIsL3rjmHU6eU85XHNvLZB+rpjoytZW4saIygP715mGmVRcx1nyJGS0kwwGcumcMftjX1zPgdj1SVv390I/W7WvjGhxbw27+5lE9dPGdEftf/edc8SgsDfPvJreO+byMeV/7xN5u5Y9UbfOTcGXxp6fyc1eXjF8xmdnUJN9xfz5cf+fOottU/vG4PhQEf1y6exYfPmc7ary9l462X87l3npL23LfPr+Wdp9Zy1+rtPVsCP7RuDy82HOGfPriAWdXFvNhwhI9fMJtPXHAS3//LRcTcds5f/8CClNcsDPj57DtO5tYPLeBPbx6hvCjAz/76Amr67G0/u6aERz9/Ed/40AKe2XaI6+9dO6bWR7OgMQjtoSjXr1jL13+9ic5w/xvMbN7Xyq9f3ctLbx7hknmTR3x+RiqfvGgOBX5h5Z/3pz94jHqkfg//8+f9/P0Vp/GpS+Zy5oxK/CPUN1RVGuRL7z2VZ7c18ZPnd47I7xgtT24+wM/X7ObGd5zM9689Z9T601IpKvDzxN9eyo3vOJnHXtnLT14YnXvbHYnx+IZ9XHnm1J6soq58cPvI3/L+02kPRfnuU9tYv6uZW1du4dJ5k/n0xXO4/aNn894z6vjqlacD3hDjx266mIc/dxFlhQMv6ffpS+by+y+9nd/+zds5Z1bqJXBEhE9dMpc7/3IRa99q5s5V2zOu90gb8oKFE000Fufmn7/CC9ubULwllRdMr+DD50znQ+ccnxzVEYry6fvW0dTmrcv/9lHuz0ioLC7g3NlV/KlhbA/f609jaxf/94mtXDC3mpsy+GSYDZ+5ZA71bzXz7Se3cu7sSSyeMzrDpLMpFlfuWLWN+XVlfPXK03MaMBLKiwq45X2n89bhDu5c9QZXLJw64tn301sPcqw7OqxBE6dPrWDZktn8Ys1ufrFmN9Mri/j3ZYvw+YSLT5nMxaf0/r+9YHrm810yXSDz6nNn8PKOI9zz3JsEfEJdRSEfe9ssCgOjN6ChLwsaGXhldwu3Pr6FTfta+fZfnMXs6hLueW4Hrx84xqrXDvJa4zH+/vLT8PmEe57bQVNbiB9cdy6FAR/vOX10Ox+TXTpvMnc+/QYtHWGq3J7YY93WxmN8+qfrOHCsm6ICH7d/9OxRe+MTEf712nNYv6uFf1+9nQdvuGBUfu9whaIxVL1P9b+s38ObTR3c/YnzRiwrGwoR4VtXn8lld/yRT/z4ZX70V+ezqJ9P2cPRHorym1f38bOXdzG9soiLhtmf+M0PL+SyBVN4eccRrjpnxglNSaPhHz+4gA17jvLDP3h9lL/d2Mh1S2ZxoDXEx5fMprJkZHYD7Y+M9/bbvhYvXqz19fVZu96R9hCX3P4Mk4qD/MMHzuDDSVlFJBbn1pVb+MWa3Vy5cCrvO2sqX3tsE+85o47/+Ph5WavDUK3f1cJH7/4TP/rEebz/rGk5rUs4GicUjVEY8HO0K0zDwXbebGqnwk3UO7m2jHA0zlX/8SJNbSE+dfFJXDxv8ogv8pjKj55t4Lu/38aTX3x7VpcMH654XFm19SDPb2/ifWdOY35dGU9tOcCdT28nEotz2RlT+M2GfSyaNYnHbro4J82i6Wza28pNP1/PoWMhfrx8Me88Nf32zOFonB/+oYHX9h8jFo9z6pRybrh0bs9mTwmxuPKpn67l+e2HKQ36ufXDC/nLxdkdnp0r8bgSjsV5assB/v6XGwm7vqGTakr44tL5FPh91JUXMqu6ZMj79ojIelVdnPY4CxoDu2v1du5Y9QZP/907mFd34kY1qsqKF9/itideI64wp6aEB2+4gFlZXjV0KKKxOIu+uYoPL5rOv3zkrJzVY82OI9z8i1c5PMBWmufNnkRcYcOeo/z4+sVctmDKKNawt9bOCBd9ZzWXL5jCvy07N2f1SNZwqI0vPrSBLfuP4fcJsaShmEvmVlMa9POHbU18+JzpfPsvzqI0Tbt6LrV0hPnET9aw43A7P/3UkrTZwLd/t5X/em4Hp00px+cTGg61cdrUch753EWs2dnMwdZuQtE4r+5u4Tcb9vMvHzmL65bMGpNBMxv2NHfS2hWhM+zNCUk0hQOcPbOSlV+4dEjXzTRojN2/rCQiciXw74Af+Imqfmc0fm8oGuOBl3bxrtNqUwYMVzduuHQuZ8+spCsc49J5k8dEOzJAwO/jwpOreWbrIfa/u4uasiDhaJzyNHtB72nu5Gcv76I7EuOiUyazYFoFM6uKB/262kNR7n62gf/64w5m15Tw2bfPJRyNU1lSwEk1pZw2pZy27girth7kqc0HiMaVv7vs1JwGDIDKkgI+ccFsfvz8ThRvZNX8urKe16+qo/aG1B2Jce8LO/nBM9spCQa482PncPmCqTy5+QDHuiIsmj2Jc2dNQkQ4eKybuvLCMf9mWVUa5MEblvCxe17m+hVr+Nr7ziAY8FFVUsB7z5hCUYGfprYQq7cepLG1m/96bgefuGA2t7kPPqu3HuSG++u58F9Wc6y794CUT18yZ0RnuY8Fs6pLSORPz/7/76KxtYu4epMOhZH/tx/zmYaI+IE38LaO3QusA65T1ddSHT/UTONYd4Rv/+51usJROsIxusIxDreHeP1AGw/esIS3z0+fRo9Ff3rzMJ+9vx6fT1D1Uv1vXb2QSEx5assBCgM+dh3ppLkjzN8unc/elk7ufWEnPhEK/D663Bjx2dUlfOE98ygvDNAWSj1yrDQY4PRp5RQGfDz3xmHuWLWNw+1hPnLuDP75qoVUpAlWY0kkFudHf3iTHzyznWhcqS0v5ItL57PurWae3HSAypICplQUUldexJSKQmrLi5g5qZilZ9RRWVzAzsMdNLWFaGhqp+FQO2+fX8vcyaU96wpFY3Eaj3Uzt6aU06dVeEtNqFJU4Gd+XTlrdzbzq1f38uy2Jpo7wly2YAq3XX3mCU0y49nRzjA3/+IVXkwarFFU4KO2vJDGo91EXTa1cHoFj37+4l6z2e957k0e37CfG99xMm+bU00w4KMw4Ev7gcj0L2+ap0TkIuAbqnqF+/4WAFX9dqrjhxo0WjsjLL3jj5QW+iku8FNaGKAk6OeU2jJu/dCCMf/pbSA7D3fwnSe3Ul0aZNeRzp6VYufVlVHg9zGlopDuSIyXd3gzyK9bMpsvLp1PdWmQTftaef3AMX728u5BryP0tjlVfP0DC0akw3O07DvaxctvHuHna3bxyu6jBAM+PnreTOJx5VBbNwePhTjUFuJIRwhVb+mI5GAL3v7t4eiJ8xP6NjMl+ATiCpNKCnjnqbV87G2zThipky8isTgb9hxlWmURu5s7eWbrIQ63h5hSWcRfnDuTqRVFlBUFxlSnfr7Kp6BxDXClqv61+/6TwAWq+oWkY24EbgSYPXv2+bt27Up5LeP9J/3pizs5eXIZS8+o6wmGqsoTmxqpLSvkgpNPbGOOx5VX9xylqMBHRVEBqWLo0c4IWxuPoQozqoq5+JSacR1sk6kqL715hJlVJcyuObG/KhKLs/1gOyv/vJ/uSIxzZlUyrbKYWdUl1JUXsnrrQY50hLl03mSKC/yICDWlQXYc7uCtwx34/YJfhNauCFv2H+PUKWV84OxpOR1aaSaWfAoa1wJX9AkaS1T1b1Idn+2OcGOMmQgyDRrjYUb4XiB53NxMYPxOczbGmHFsPASNdcB8EZkrIkFgGbAyx3UyxpgJacwPuVXVqIh8AXgKb8jtClXdkuNqGWPMhDTmgwaAqv4O+F2u62GMMRPdeGieMsYYM0ZY0DDGGJMxCxrGGGMyZkHDGGNMxsb85L7BEpE2YFuf4kogm/ueZvt6AJOBw1m6VjbrNxKvdSSum837B2P7b8bu3di6Xj7cv8lAqaqmX2RPVfPqAdSnKLsny78jq9frr95joX4j8VpH6N8ka/dvhOo3Zv9NJtK9s/s3/NcwUZqn/meMXy/bslm/kXqtE+keZvt6du/G1vWybUy/3nxsnqrXDNZPGWvGa73HCrt/Q2f3bnjy4f4N5jXkY6ZxT64rMETjtd5jhd2/obN7Nzz5cP8yfg15l2kYY4wZOfmYaRhjjBkhFjRGiIjMEpE/iMhWEdkiIl905dUiskpEtruvVa68xh3fLiI/7HOt60Rkk4hsFJHfi0h+buOWJMv372Pu3m0Rke/m4vWMpiHcu8tEZL37G1svIu9Jutb5rrxBRO6SfNlVawBZvn+3icgeEWnP1evJumwO7bJHryFs04Dz3PNyvH3OFwDfBb7myr8G3O6elwKXAp8Hfph0nQBwCJjsvv8u3va3OX+N4+T+1QC7gVr3/f3A0ly/vjF2784FprvnZwL7kq61FrgIEOBJ4H25fn3j7P5d6K7XnuvXla2HZRojRFUbVfUV97wN2ArMAK7Ce+PCfb3aHdOhqi8A3X0uJe5R6j7lVTABNqHK4v07GXhDVZvc908DHx3h6ufUEO7dq6qa+JvaAhSJSKGITAMqVPUl9d4BH0ick8+ydf/cz15W1cbRrP9Is6AxCkRkDt6nkTXAlMQfkftaN9C5qhoBbgI24QWLBcC9I1jdMWc49w9oAE4XkTkiEsD7jz4rzTl5Ywj37qPAq6oawnuj3Jv0s72ubMIY5v3LSxY0RpiIlAGPAV9S1WNDOL8AL2icC0wHNgK3ZLWSY9hw75+qtuDdv4eB54G3gGg26zhWDfbeichC4Hbgc4miFIdNmOGWWbh/ecmCxghyb/iPAT9X1V+54oMu7cd9PZTmMosAVPVN10TwCHDxCFV5TMnS/UNV/0dVL1DVi/DWJds+UnUeKwZ770RkJvBr4HpVfdMV7wVmJl12JhOgaRSydv/ykgWNEeL6H+4FtqrqHUk/Wgksd8+XA4+nudQ+YIGIJBYSuwyvjTWvZfH+ISJ17msV8H+An2S3tmPLYO+diEwCngBuUdUXEwe7Jpg2EbnQXfN6Mrjf41227l/eynVPfL4+8EbyKF5z0gb3eD/eaJ7VeJ92VwPVSee8BTQD7Xif8ha48s/jBYqNeOvI1OT69Y2z+/ffwGvusSzXr22s3TvgH4GOpGM3AHXuZ4uBzcCbwA9xE4Lz+ZHl+/dd97cYd1+/kevXN9yHzQg3xhiTMWueMsYYkzELGsYYYzJmQcMYY0zGLGgYY4zJmAUNY4wxGbOgYcwoE5HPi8j1gzh+johsHsk6GZOpQK4rYMxEIiIBVf3PXNfDmKGyoGHMILlF7H6Pt4jduXhLZ18PnAHcAZQBh4FPqWqjiDwL/Am4BOAxjMQAAAFfSURBVFgpIuV4S2V/T0QWAf8JlOBNoPuMqraIyPnACqATeGH0Xp0xA7PmKWOG5jTgHlU9GzgG3Az8ALhGVRNv+LclHT9JVd+pqt/vc50HgK+662wCbnXlPwX+Vr31sowZMyzTMGZo9ujxdYZ+BvwD3gY8q9zmdn4geR+Fh/teQEQq8YLJH13R/cAvU5Q/CLwv+y/BmMGzoGHM0PRdf6cN2DJAZtAxiGtLiusbMyZY85QxQzNbRBIB4jrgZaA2USYiBW5/hX6paivQIiJvd0WfBP6oqkeBVhG51JV/IvvVN2ZoLNMwZmi2AstF5L/wVj39AfAUcJdrXgoA/4a3/edAlgP/KSIlwA7g067808AKEel01zVmTLBVbo0ZJDd66reqemaOq2LMqLPmKWOMMRmzTMMYY0zGLNMwxhiTMQsaxhhjMmZBwxhjTMYsaBhjjMmYBQ1jjDEZs6BhjDEmY/8PCnaTPjcjZ9gAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -252,10 +2300,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 11, + "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", @@ -274,7 +2320,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +2344,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 13, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAD8CAYAAABQFVIjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X+QldWd5/H3B0FxElEg6PBDhY3GCppZHLrQLfdHxB0gPyrgjJmwOkrVWEXijy1nNrWiE7fMKFMVU5M4y7ohMWNG1FF0TCzdRIa0UWucWQSaoFE0TjMrgygjWI2KUwVj43f/eL53fLhpbt9ubve9t/vzqrrVT5/nnHNPPzT9vefH8xxFBGZmZkdrTLMbYGZmI4MDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNYQDipmZNcTYZjdgOH3sYx+LmTNnNrsZZmZtZcuWLW9FxJT+8o2qgDJz5ky6urqa3Qwzs7Yi6R/ryechLzMzawgHFDMzawgHFDMza4i6AoqkHZJekPScpK5M+7qk1zPtOUmfLeW/UdJ2Sa9IWlhKn5v1bJe0SpIy/ThJD2b6RkkzS2WWSerO17JS+qzM251ljz36y2FmZoM1kB7KhRExJyI6Smm3Z9qciHgcQNJsYClwNrAI+I6kYzL/amA5cGa+FmX6lcC+iDgDuB24LeuaBNwMnAfMA26WNDHL3JbvfyawL+swM7MmGYohr8XA2og4GBGvAtuBeZKmAhMiYkMUu3rdAywplVmTxw8DF2XvZSHQGRE9EbEP6AQW5bn5mZcsW6mrre159wC/+70N7Nl/oNlNMTMbkHoDSgA/lbRF0vJS+rWSfiHpB6Wew3TgtVKeXZk2PY+r0w8rExG9wDvA5Bp1TQbezrzVdbW1VT/rZvOOHlY90d3sppiZDUi996FcEBFvSDoZ6JT0S4rhq1spgs2twLeA3wfUR/mokc4gytSq6zAZAJcDnHbaaX1laQln3bSOg70f/Ov3923cyX0bd3Lc2DG8svIzTWyZmVl96uqhRMQb+XUP8AgwLyLejIhDEfEB8H2KOQ4oegunlorPAN7I9Bl9pB9WRtJY4ESgp0ZdbwEnZd7quqrbfmdEdEREx5Qp/d7o2TTPXH8hX5gzjfHjin+S8ePGsHjONJ5ZcWGTW2ZmVp9+A4qkj0g6oXIMLABezDmRiouBF/P4MWBprtyaRTH5vikidgP7JZ2fcyBXAI+WylRWcF0CPJnzLOuBBZIm5pDaAmB9nnsq85JlK3W1pZMnjOeE48ZysPcDjhs7hoO9H3DCcWM5+YTxzW6amVld6hnyOgV4JFf4jgXuj4i/lnSvpDkUQ007gC8DRMQ2SQ8BLwG9wDURcSjrugq4GzgeWJcvgLuAeyVtp+iZLM26eiTdCmzOfLdERE8erwDWSloJbM062tpb7x3ksvNO59J5p3H/pp3s9cS8mbURFR/2R4eOjo7ws7zMzAZG0paqW0b65DvlzcysIRxQzMysIRxQzMysIRxQzMysIRxQzMysIRxQ2pCf92VmrcgBpQ35eV9m1opG1Z7y7c7P+zKzVuYeShvx877MrJU5oLQRP+/LzFqZh7zajJ/3ZWatys/yqsOedw9w7QNbuePSc90bMLNRx8/yaiCvqjIz65+HvGrwqiozs/q5h1KDV1WZmdXPAaUGr6oyM6tfXQFF0g5JL0h6TlJXpk2S1CmpO79OLOW/UdJ2Sa9IWlhKn5v1bJe0KrcCJrcLfjDTN0qaWSqzLN+jW9KyUvqszNudZY89+svxqyqrqh65+gIuO+909r53cCjexsys7dW1ykvSDqAjIt4qpX0T6ImIb0i6AZgYESskzQYeAOYB04AngE9ExCFJm4DrgGeBx4FVEbFO0tXAb0TEVyQtBS6OiC9JmgR0AR0UWw1vAeZGxL7cZvhHEbFW0neB5yNida2fwzs2mpkN3HCs8loMrMnjNcCSUvraiDgYEa8C24F5kqYCEyJiQxRR7J6qMpW6HgYuyt7LQqAzInoiYh/QCSzKc/Mzb/X7m5lZE9QbUAL4qaQtkpZn2ikRsRsgv56c6dOB10pld2Xa9DyuTj+sTET0Au8Ak2vUNRl4O/NW12VmZk1Q77LhCyLiDUknA52Sflkjr/pIixrpgylTq67DG1MEwOUAp512Wl9ZzMysAerqoUTEG/l1D/AIxfzImzmMRX7dk9l3AaeWis8A3sj0GX2kH1ZG0ljgRKCnRl1vASdl3uq6qtt+Z0R0RETHlClT6vlxzcxsEPoNKJI+IumEyjGwAHgReAyorLpaBjyax48BS3Pl1izgTGBTDovtl3R+zoFcUVWmUtclwJM5z7IeWCBpYq4iWwCsz3NPZd7q9zczsyaoZ8jrFOCRXOE7Frg/Iv5a0mbgIUlXAjuBLwJExLZcgfUS0AtcExGHsq6rgLuB44F1+QK4C7hX0naKnsnSrKtH0q3A5sx3S0T05PEKYK2klcDWrMPMzJrED4c0M7Oa/HBIMzMbVg4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEA4oZmbWEHUHFEnHSNoq6cf5/dclvS7puXx9tpT3RknbJb0iaWEpfa6kF/Lcqtxbntx//sFM3yhpZqnMMknd+VpWSp+Vebuz7LFHdynMzOxoDKSHch3wclXa7RExJ1+PA0iaTbEn/NnAIuA7ko7J/KuB5cCZ+VqU6VcC+yLiDOB24LasaxJwM3AeMA+4WdLELHNbvv+ZwL6sw8zMmqSugCJpBvA54M/ryL4YWBsRByPiVWA7ME/SVGBCRGyIYiP7e4AlpTJr8vhh4KLsvSwEOiOiJyL2AZ3Aojw3P/OSZSt1mZlZE9TbQ/kz4Hrgg6r0ayX9QtIPSj2H6cBrpTy7Mm16HlenH1YmInqBd4DJNeqaDLydeavrMjOzJug3oEj6PLAnIrZUnVoNfByYA+wGvlUp0kc1USN9MGVq1XUYScsldUnq2rt3b19ZzMysAerpoVwAfEHSDmAtMF/SfRHxZkQciogPgO9TzHFA0Vs4tVR+BvBGps/oI/2wMpLGAicCPTXqegs4KfNW13WYiLgzIjoiomPKlCl1/LhmZjYY/QaUiLgxImZExEyKyfYnI+L3ck6k4mLgxTx+DFiaK7dmUUy+b4qI3cB+SefnHMgVwKOlMpUVXJfkewSwHlggaWIOqS0A1ue5pzIvWbZSl5mZNcHY/rMc0TclzaEYatoBfBkgIrZJegh4CegFromIQ1nmKuBu4HhgXb4A7gLulbSdomeyNOvqkXQrsDnz3RIRPXm8AlgraSWwNeswM7MmUfFhf3To6OiIrq6uZjfDzKytSNoSER395fOd8mZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGZm1hAOKGY2Kux59wC/+70N7Nl/oNlNGbEcUMxsVFj1s2427+hh1RPdzW7KiHU0z/IyM2t5Z920joO9H27ldN/Gndy3cSfHjR3DKys/08SWjTzuoZjZiPbM9RfyhTnTGD+u+HM3ftwYFs+ZxjMrLmxyy0YeBxQza5rhmNc4ecJ4TjhuLAd7P+C4sWM42PsBJxw3lpNPGD9k7zlaOaCYWdMM17zGW+8d5LLzTueRqy/gsvNOZ+97B4f0/UYrP77ezIZd9bxGhec1WpMfX29mLcvzGiOTA4qZDTvPa4xMdQcUScdI2irpx/n9JEmdkrrz68RS3hslbZf0iqSFpfS5kl7Ic6tyb3ly//kHM32jpJmlMsvyPbolLSulz8q83Vn22KO7FGY2nDyvMfLUPYci6b8BHcCEiPi8pG8CPRHxDUk3ABMjYoWk2cADwDxgGvAE8ImIOCRpE3Ad8CzwOLAqItZJuhr4jYj4iqSlwMUR8SVJk4CufN8AtgBzI2Jf7lv/o4hYK+m7wPMRsbrWz+A5FDOzgWvoHIqkGcDngD8vJS8G1uTxGmBJKX1tRByMiFeB7cA8SVMpgtGGKKLYPVVlKnU9DFyUvZeFQGdE9ETEPqATWJTn5mfe6vc3M7MmqHfI68+A64HysoxTImI3QH49OdOnA6+V8u3KtOl5XJ1+WJmI6AXeASbXqGsy8Hbmra7rMJKWS+qS1LV37946f1wzMxuofgOKpM8DeyJiS511qo+0qJE+mDK16jo8MeLOiOiIiI4pU6b0lcXMzBqgnh7KBcAXJO0A1gLzJd0HvJnDWOTXPZl/F3BqqfwM4I1Mn9FH+mFlJI0FTgR6atT1FnBS5q2uy6wmP3XWbGj0G1Ai4saImBERM4GlwJMR8XvAY0Bl1dUy4NE8fgxYmiu3ZgFnAptyWGy/pPNzDuSKqjKVui7J9whgPbBA0sRcRbYAWJ/nnsq81e9vVpOfOms2NI7macPfAB6SdCWwE/giQERsyxVYLwG9wDURcSjLXAXcDRwPrMsXwF3AvZK2U/RMlmZdPZJuBTZnvlsioiePVwBrJa0EtmYdZkfkp86aDS0/esVGjT3vHmDl4y/z023/xIH3P2D8uDEsPPvX+drnPukb6sxq8KNXzKoM5O7sdpxnacc228jigGKjSr13Z7fjPEs7ttlGFg95mZW041Nw27HN1l485GU2CO34FNx2bLONTA4oZiXt+BTcdmyzjUxHs2zYbESqzLNcOu807t+0k71tMMndjm22kcdzKGZmVpPnUMzMbFg5oJiZWUM4oJiZDYBvID0yBxQzswHwDaRH5lVe1lb2vHuAax/Yyh2XnutlsTas/HDR/rmHYm3Fnw6tWXwDaf/cQ7G2MNo/Hbpn1ny+gbR/7qFYWxjtnw7dM2sN9T5cdLRyD8VaRq1P4aP10+Fo75m1mu9d/uG9fSuXnNPElrSmfnsoksZL2iTpeUnbJP1xpn9d0uuSnsvXZ0tlbpS0XdIrkhaW0udKeiHPrcqtgMntgh/M9I2SZpbKLJPUna9lpfRZmbc7yx7bmEtiA9WoZZT9fQofjZ8OR3vPzNpLPT2Ug8D8iHhP0jjgbyVVtu69PSL+tJxZ0myKLXzPBqYBT0j6RG4DvBpYDjwLPA4sotgG+EpgX0ScIWkpcBvwJUmTgJuBDiCALZIei4h9mef2iFgr6btZx+rBXwobrHIgWHnxpwZcvt5P4aPx0+Fo7ZlZe+q3hxKF9/Lbcfmq9QCwxcDaiDgYEa8C24F5kqYCEyJiQxQPELsHWFIqsyaPHwYuyt7LQqAzInoyiHQCi/Lc/MxLlq3UZcPkrJvWMfOGn3Dfxp1EFIFg5g0/4ayb1vVfuMSfwmsbjT0za091zaFIOgbYApwB/O+I2CjpM8C1kq4AuoCv5h/96RQ9kIpdmfZ+Hlenk19fA4iIXknvAJPL6VVlJgNvR0RvH3XZMHnm+guPuEf7QPhTeG2jsWdm7amuVV4RcSgi5gAzKHob51AML30cmAPsBr6V2dVXFTXSB1OmVl2HkbRcUpekrr179/aVxQapkYHAn8JHHj+iZPQZ0CqviHhb0tPAovLciaTvAz/Ob3cBp5aKzQDeyPQZfaSXy+ySNBY4EejJ9E9XlXkaeAs4SdLY7KWU66pu853AnVA8vn4gP6/1r1H7cPhT+MhztHNr1n763Q9F0hTg/QwmxwM/pZgQ3xIRuzPPHwLnRcRSSWcD9wPzKCblfwacGRGHJG0G/iuwkWJS/n9FxOOSrgE+FRFfyUn5346I381J+S3Ab2Zzfg7MjYgeSX8F/LA0Kf+LiPhOrZ/F+6GYDT3vcT/y1LsfSj09lKnAmpxHGQM8FBE/lnSvpDkUQ007gC8DRMQ2SQ8BLwG9wDW5wgvgKuBu4HiK1V2V2du7gHslbafomSzNunok3Qpszny3RERPHq8A1kpaCWzNOsysyRo1t2btp9+AEhG/AM7tI/3yGmX+BPiTPtK7gF8Zz4iIA8AXj1DXD4Af9JH+/yh6QWbWQrzIYvTynfJm1nDe43508p7yZmZWk/eUNzOzYeWAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYmZmDeGAYjaEvCeIjSYOKGZDqLwniNlI54dDmg2B6j1B7tu4k/s27vSeIDaiuYdiNgSeuf5CvjBnGuPHFf/Fxo8bw+I503hmxYVNbpnZ0HFAMRsC3hPERiMHFBtRWmkSvLInyCNXX8Bl553O3vcONrtJZkOq34AiabykTZKel7RN0h9n+iRJnZK68+vEUpkbJW2X9IqkhaX0uZJeyHOrJCnTj5P0YKZvlDSzVGZZvke3pGWl9FmZtzvLHtuYS2LtrJUmwb93eQcrl5zD7GkTWLnkHL53eb/bSZi1tX432Mo/+h+JiPckjQP+FrgO+G2gJyK+IekGYGJErJA0G3iAYnveacATwCci4pCkTVn2WeBxYFVErJN0NfAbEfEVSUuBiyPiS5ImAV1AB8Xe9VuAuRGxL/et/1FErJX0XeD5iFhd62fxBlsjV/UkeIUnwc2OXsM22IrCe/ntuHwFsBhYk+lrgCV5vBhYGxEHI+JVYDswT9JUYEJEbIgiit1TVaZS18PARRnIFgKdEdETEfuATmBRnpufeavf30YhT4IPr1YaWrTWUdcciqRjJD0H7KH4A78ROCUidgPk15Mz+3TgtVLxXZk2PY+r0w8rExG9wDvA5Bp1TQbezrzVdVW3fbmkLklde/furefHtTbkSfDh1UpDi9Y66roPJSIOAXMknQQ8IumcGtnVVxU10gdTplZdhydG3AncCcWQV195bGSoTIJfOu807t+0k73+9Nxwvr/GahnQjY0R8bakp4FFwJuSpkbE7hzO2pPZdgGnlorNAN7I9Bl9pJfL7JI0FjgR6Mn0T1eVeRp4CzhJ0tjspZTrslGqPOm9ckmtzzw2WM9cfyErH3+Zn277Jw68/wHjx41h4dm/ztc+98lmN81aQD2rvKZkzwRJxwP/Gfgl8BhQWXW1DHg0jx8DlubKrVnAmcCmHBbbL+n8nAO5oqpMpa5LgCdznmU9sEDSxFxFtgBYn+eeyrzV729mQ8RDi1ZLPT2UqcAaScdQBKCHIuLHkjYAD0m6EtgJfBEgIrblCqyXgF7gmhwyA7gKuBs4HliXL4C7gHslbafomSzNunok3Qpszny3RERPHq8A1kpaCWzNOsxsiHlo0Y6k32XDI8lQLhve8+4Brn1gK3dceq4/rZk1UDv+32rHNtfSsGXDVh+vejEbGu34f6sd29wI7qEcJd9QZzY02vH/Vju2uR7uoQwT31BnNjTa8f9WO7a5kRxQjtJAVr347mJrB63ye9qOK8rasc2N5IDSAPU+VXa0jqtae2ml39N2fGJzO7a5UTyHMgxG6riqjSz+PbUj8RxKCxnt46rWHvx72npaZfixXg4ow6Cdx1Xb7RfaBq+df09HqlYafqzHgJ7lZYPXrncXl3+hV178qWY3x4ZYu/6ejjTt+hBOz6FYnzyebtY8e949cMSHcDajx+g5FDsqHk83a552HX70kJf1qV1/oc1GinYcfnRAsSNqx19os5GiHff38RzKKDbSnohqZkPDcyjWr3Zbkmhmrc1DXqNQuy5JbDXu4Zkdrp4tgE+V9JSklyVtk3Rdpn9d0uuSnsvXZ0tlbpS0XdIrkhaW0udKeiHPrcqtgMntgh/M9I2SZpbKLJPUna9lpfRZmbc7yx7bmEsy8nkFV2O4h2d2uHp6KL3AVyPi55JOALZI6sxzt0fEn5YzS5pNsYXv2cA04AlJn8htgFcDy4FngceBRRTbAF8J7IuIMyQtBW4DviRpEnAz0AFEvvdjEbEv89weEWslfTfrWD34SzF6eAXX0XEPz6xv/fZQImJ3RPw8j/cDLwPTaxRZDKyNiIMR8SqwHZgnaSowISI2RLES4B5gSanMmjx+GLgoey8Lgc6I6Mkg0gksynPzMy9ZtlKX1WE0PxH1aLmHZ9a3Ac2h5FDUucBG4ALgWklXAF0UvZh9FMHm2VKxXZn2fh5Xp5NfXwOIiF5J7wCTy+lVZSYDb0dEbx91WR3acUliq3APz6xvda/ykvRR4IfAH0TEuxTDSx8H5gC7gW9VsvZRPGqkD6ZMrbqq271cUpekrr179/aVxWzA3MMz+1V19VAkjaMIJn8ZET8CiIg3S+e/D/w4v90FnFoqPgN4I9Nn9JFeLrNL0ljgRKAn0z9dVeZp4C3gJEljs5dSruswEXEncCcU96HU8/Oa9cc9PLNfVc8qLwF3AS9HxLdL6VNL2S4GXszjx4CluXJrFnAmsCkidgP7JZ2fdV4BPFoqU1nBdQnwZM6zrAcWSJooaSKwAFif557KvGTZSl1mZtYE9fRQLgAuB16Q9Fym/RHwXyTNoRhq2gF8GSAitkl6CHiJYoXYNbnCC+Aq4G7geIrVXesy/S7gXknbKXomS7OuHkm3Apsz3y0R0ZPHK4C1klYCW7MOMzNrEj96xczMavKjV8zMbFg5oJiZNcFI3F7bAcVsFBmJf8Ta1Uh8dI8fDmk2ipT/iK28+FPNbs6oNJIf3eNJebNRoPqPWMVI+CPWbhq5X/xwPfHak/Jm9q/8/LHW0chH97TasJmHvMyabDg+Zfr5Y63laLfXbtVhMwcUsyYbrnmNo/0jZo1ztI/ueeb6C484bNZMDihmTTLcnzL9/LGRo1V7nJ5DMWsSz2s0zmhcDt2KT7x2D8WsSVr1U2Y7Go3LoVuxx+mAYtZEntc4Oq06OT1a+T4UM2tbjbynw47M96GY2YjnYcPW4iEvM2trHjZsHR7yMjMbwRpx46yHvMzMbFgfz1LPnvKnSnpK0suStkm6LtMnSeqU1J1fJ5bK3Chpu6RXJC0spc+V9EKeW5V7y5P7zz+Y6RslzSyVWZbv0S1pWSl9VubtzrLHNuaSmJm1v7NuWsfMG37CfRt3ElGsgJt5w08466Z1/RcepHp6KL3AVyPik8D5wDWSZgM3AD+LiDOBn+X35LmlwNnAIuA7ko7JulYDy4Ez87Uo068E9kXEGcDtwG1Z1yTgZuA8YB5wcylw3Qbcnu+/L+swMzOac+NsvwElInZHxM/zeD/wMjAdWAysyWxrgCV5vBhYGxEHI+JVYDswT9JUYEJEbIhi4uaeqjKVuh4GLsrey0KgMyJ6ImIf0AksynPzM2/1+5uZjXrNWAE3oFVeORR1LrAROCUidkMRdCSdnNmmA8+Wiu3KtPfzuDq9Uua1rKtX0jvA5HJ6VZnJwNsR0dtHXWZmxvCvgKs7oEj6KPBD4A8i4t2c/ugzax9pUSN9MGVq1XV4Y6TlFMNsnHbaaX1lMTMbkYb78Sx1rfKSNI4imPxlRPwok9/MYSzy655M3wWcWio+A3gj02f0kX5YGUljgROBnhp1vQWclHmr6zpMRNwZER0R0TFlypR6flwzMxuEelZ5CbgLeDkivl069RhQWXW1DHi0lL40V27Noph835TDY/slnZ91XlFVplLXJcCTOc+yHlggaWJOxi8A1ue5pzJv9fubmVkT1DPkdQFwOfCCpOcy7Y+AbwAPSboS2Al8ESAitkl6CHiJYoXYNRFxKMtdBdwNHA+syxcUAeteSdspeiZLs64eSbcCmzPfLRHRk8crgLWSVgJbsw4zM2sS3ylvZmY1+U55aymjcQMks9HGAcWGxXA+/sHMmsNPG7Yh5Q2QzEYP91BsSHnfdLPRwwHFhpQ3QDIbPTzkZUPOGyCZjQ5eNmxmZjV52bCZmQ0rBxQzM2sIB5QRyjcSmtlwc0AZoXwjoZkNN6/yGmF8I6GZNYt7KCOMbyQ0GzwPFR8dB5QRxjcSmg2eh4qPjoe8RiDfSGg2MB4qbgzf2Ghmo96edw+w8vGX+em2f+LA+x8wftwYFp7963ztc590754G3tgo6QeS9kh6sZT2dUmvS3ouX58tnbtR0nZJr0haWEqfK+mFPLcqtwEmtwp+MNM3SppZKrNMUne+lpXSZ2Xe7ix7bD0XxcysLx4qbox65lDuBhb1kX57RMzJ1+MAkmZTbN97dpb5jqRjMv9qYDnFHvNnluq8EtgXEWcAtwO3ZV2TgJuB84B5wM25rzyZ5/aIOBPYl3WYmQ1aZaj4kasv4LLzTmfveweb3aS20+8cSkT8TbnX0I/FwNqIOAi8mnvEz5O0A5gQERsAJN0DLKHYU34x8PUs/zBwR/ZeFgKdlT3kJXUCiyStBeYDl2aZNVl+dZ1tNDP7Fd+7/MMRnZVLzmliS9rX0azyulbSL3JIrNJzmA68VsqzK9Om53F1+mFlIqIXeAeYXKOuycDbmbe6LjMza5LBBpTVwMeBOcBu4FuZrj7yRo30wZSpVdevkLRcUpekrr179x4pm5mZHaVBBZSIeDMiDkXEB8D3KeY4oOgtnFrKOgN4I9Nn9JF+WBlJY4ETgZ4adb0FnJR5q+vqq613RkRHRHRMmTJloD+qmZnVaVABRdLU0rcXA5UVYI8BS3Pl1iyKyfdNEbEb2C/p/JwfuQJ4tFSmsoLrEuDJKNYyrwcWSJqYQ2oLgPV57qnMS5at1GVmZk3S76S8pAeATwMfk7SLYuXVpyXNoRhq2gF8GSAitkl6CHgJ6AWuiYhDWdVVFCvGjqeYjF+X6XcB9+YEfg/FKjEiokfSrcDmzHdLZYIeWAGslbQS2Jp1mJlZE/nGRjMzq6neGxtHVUCRtBf4xz5OfYxibqaduM1Dr93aC27zcGm3Nh9te0+PiH4noUdVQDkSSV31RN9W4jYPvXZrL7jNw6Xd2jxc7fXThs3MrCEcUMzMrCEcUAp3NrsBg+A2D712ay+4zcOl3do8LO31HIqZmTWEeyhmZtYQIzKgHGEPl38raUPuyfJ/JE3I9HGS1mT6y5JuLJV5Ovd1qez7cnKLtPlYSX+R6c9L+nSpTJ/7zrR4m4flOks6VdJT+e+8TdJ1mT5JUmfur9NZetjpgPf3afE2t+R1ljQ5878n6Y6qulryOvfT5iG/zoNo729J2pLXcouk+aW6GneNI2LEvYD/CPwm8GIpbTPwn/L494Fb8/hSikfuA/waxZ3/M/P7p4GOFmzzNcBf5PHJwBZgTH6ruEO0AAADpklEQVS/Cfh3FA/RXAd8pg3aPCzXGZgK/GYenwD8PTAb+CZwQ6bfANyWx7OB54HjgFnAPwDHDOd1bnCbW/U6fwT498BXgDuq6mrV61yrzUN+nQfR3nOBaXl8DvD6UFzjEdlDiYi/oXiMS9lZwN/kcSfwO5XswEdUPGzyeOBfgHeHo51lA2zzbOBnWW4P8DbQoeIZaxMiYkMUvymVfWdats1D1ba+RMTuiPh5Hu8HXqbY+mAxxb465NfKNfvX/X0i4lWgsr/PsF3nRrV5KNrWqDZHxD9HxN8CB8r1tPJ1PlKbh8sg2rs1IioP0d0GjFfxzMWGXuMRGVCO4EXgC3n8RT58kvHDwD9TPIZ/J/Cn8eEzwwD+Irut/2Moh4+O4Ehtfh5YLGmsiodwzs1ztfadGS4DbXPFsF5nFZvGnQtsBE6J4gGm5NfKEMVg9vcZMkfZ5opWvM5H0srXuT/Ddp0H0d7fAbZGsRFiQ6/xaAoovw9cI2kLRRfxXzJ9HnAImEYxRPBVSf8mz10WEZ8C/kO+Lh/eJh+xzT+g+IfvAv4M+L8UD+Mc0F4xQ2SgbYZhvs6SPgr8EPiDiKjVG23InjyN0IA2Q+te5yNW0Udaq1znWobtOg+0vZLOpthC/cuVpD6yDfoaj5qAEhG/jIgFETEXeIBibBmKOZS/joj3cyjm78ihmIh4Pb/uB+5n+IcO+mxzRPRGxB9GxJyIWAycBHRTe9+ZVm3zsF5nSeMo/gP+ZUT8KJPfzK5/ZZhlT6YPZn+fVm1zK1/nI2nl63xEw3WdB9peSTOAR4ArIqLy96+h13jUBJTKSgtJY4CbgO/mqZ3AfBU+ApwP/DKHZj6WZcYBn+fDfV+a2mZJv5ZtRdJvAb0R8VLU3nemJds8nNc5r8ldwMsR8e3SqfKePOX9dQazv09LtrnFr3OfWvw6H6meYbnOA22vpJOAnwA3RsTfVTI3/BoPdja/lV8Un4x3A+9TROArgesoVkL8PfANPryp86PAX1FMVL0E/Pf4cBXHFuAXee5/kqtlWqDNM4FXKCbinqB4Emilng6KX+B/AO6olGnVNg/ndaZYlRP5Xs/l67PAZIoFA935dVKpzNfyWr5CafXLcF3nRrW5Da7zDooFHu/l79LsNrjOv9Lm4brOA20vxYe7fy7lfQ44udHX2HfKm5lZQ4yaIS8zMxtaDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQDihmZtYQ/x/rFzgxQKYDrwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.plot(style='*')" ] @@ -314,9 +2383,55 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "2014 1600941\n", + "1991 1659249\n", + "1995 1840410\n", + "2020 2053781\n", + "2012 2175217\n", + "2003 2234584\n", + "2019 2254386\n", + "2006 2307352\n", + "2017 2321583\n", + "2001 2529279\n", + "1992 2574578\n", + "1993 2703886\n", + "2018 2705325\n", + "1988 2765617\n", + "2007 2780164\n", + "1987 2855570\n", + "2016 2856393\n", + "2011 2857040\n", + "2008 2973918\n", + "1998 3034904\n", + "2002 3125418\n", + "2009 3444020\n", + "1994 3514763\n", + "1996 3539413\n", + "2004 3567744\n", + "1997 3620066\n", + "2015 3654892\n", + "2000 3826372\n", + "2005 3835025\n", + "1999 3908112\n", + "2010 4111392\n", + "2013 4182691\n", + "1986 5115251\n", + "1990 5235827\n", + "1989 5466192\n", + "dtype: int64" + ] + }, + "execution_count": 14, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,9 +2446,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 15, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEKCAYAAAAGvn7fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGbhJREFUeJzt3X2UJXV95/H3h5kBhmkYjAONDkr7QAjIqDgXXWQ13WhcdNCcGE5QQcVIGo0PRCdnM8v6sLrLOj5MsphgkklUiAodw8PZyBjUE2hQNEgPqA2OEBdmlSEMAjLSMAuMfPePX7XctP1wq27dvsXPz+ucPn3vrbpVn/rdut9b9auqexURmJlZPvbqdwAzM6uXC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLzNJeTnzVqlUxNDQ067AHH3yQFStW9HL2lTU5GzQ7n7NV1+R8zlZd2Xxbt269JyIO6mqmEdGzv7Vr18ZcrrrqqjmH9VuTs0U0O5+zVdfkfM5WXdl8wER0WXvdFWNmlhkXdjOzzLiwm5llxoXdzCwzLuxmZpkpVdglvUfSzZJuknSRpH17FczMzKrpuLBLWg28G2hFxNHAEuB1vQpmZmbVlO2KWQosl7QU2A+4s/5IZmbWDUWJ3zyVdBZwDrAb+GpEnDrLOKPAKMDg4ODasbGxWac1NTXFwMBAlcw91+RsUH++yR27apvW4HLYubvz8desXlnbvBfyq/a61snZqiubb2RkZGtEtLqZZ8eFXdKTgEuAU4D7gX8ALo6Iz8/1nFarFRMTE7MOGx8fZ3h4uGzeRdHkbFB/vqENW2qb1vo1e9g02fk3VWzfuK62eS/kV+11rZOzVVc2n6SuC3uZrpiXA7dHxE8i4lHgUuDF3czczMzqV6aw/wj4D5L2kyTgZcC23sQyM7OqOi7sEXEdcDFwAzBZPHdzj3KZmVlFpb62NyI+CHywR1nMzKwGvvLUzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMdF3ZJR0j6TtvfzyT9US/DmZlZeR3/NF5E3AI8H0DSEmAHcFmPcpmZWUVVu2JeBvyfiPi/dYYxM7PuKSLKP0n6DHBDRPzFLMNGgVGAwcHBtWNjY7NOY2pqioGBgdLzXgxNzgb155vcsau2aQ0uh527Ox9/zeqVtc17Ie3tVucylzHf8jZ5vXO26srmGxkZ2RoRrW7mWbqwS9obuBN4TkTsnG/cVqsVExMTsw4bHx9neHi41LwXS5OzQf35hjZsqW1a69fsYdNkxz18bN+4rrZ5L6S93epc5jLmW94mr3fOVl3ZfJK6LuxVumJeSdpan7eom5lZf1Qp7K8HLqo7iJmZ1aNUYZe0H/BbwKW9iWNmZt3qvDMUiIiHgCf3KIuZmdXAV56amWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZabsT+MdKOliST+QtE3Scb0KZmZm1ZT6aTzgXOCKiDhZ0t7Afj3IZGZmXei4sEs6AHgpcDpARDwCPNKbWGZmVpUiorMRpecDm4HvA88DtgJnRcSDM8YbBUYBBgcH146Njc06vampKQYGBqon75HJHbsYXA47dy/+vNesXtnReHW33eSOXbVNq2zbdbrMdWhvtzqXuYz5lrep7wlwtm6UzTcyMrI1IlrdzLNMYW8B/wIcHxHXSToX+FlEvH+u57RarZiYmJh12Pj4OMPDw+UT99jQhi2sX7OHTZNle6m6t33juo7Gq7vthjZsqW1aZduu02WuQ3u71bnMZcy3vE19T4CzdaNsPkldF/YyB0/vAO6IiOuK+xcDL+hm5mZmVr+OC3tE3AX8WNIRxUMvI3XLmJlZg5Ttb3gX8IXijJjbgLfUH8nMzLpRqrBHxHeArvp+zMyst3zlqZlZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmSv2CkqTtwAPAz4E93f6StpmZ1a/sb54CjETEPbUnMTOzWrgrxswsM4qIzkeWbgd+CgTw1xGxeZZxRoFRgMHBwbVjY2OzTmtqaoqBgYEqmXtqcscuBpfDzt2LP+81q1d2NF7dbTe5Y1dt0yrbdp0ucx3a263OZS5jvuVt6nsCnK0bZfONjIxs7babu2xhf2pE3CnpYOBrwLsi4pq5xm+1WjExMTHrsPHxcYaHh0vG7b2hDVtYv2YPmyar9FJ1Z/vGdR2NV3fbDW3YUtu0yrZdp8tch/Z2q3OZy5hveZv6ngBn60bZfJK6LuylumIi4s7i/93AZcALu5m5mZnVr+PCLmmFpP2nbwOvAG7qVTAzM6umTH/DIHCZpOnnXRgRV/QklZmZVdZxYY+I24Dn9TCLmZnVwKc7mpllxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWWmdGGXtETSjZIu70UgMzPrTpUt9rOAbXUHMTOzepQq7JIOBdYBf9ubOGZm1i1FROcjSxcDHwH2B/44Ik6aZZxRYBRgcHBw7djY2KzTmpqaYmBgoErmnprcsYvB5bBzd7+TzK3J+cpmW7N6Ze/CzNC+zk3u2LVo82033/I29T0BztaNsvlGRka2RkSrm3ku7XRESScBd0fEVknDc40XEZuBzQCtViuGh2cfdXx8nLmG9dPpG7awfs0eNk123DSLrsn5ymbbfupw78LM0L7Onb5hy6LNt918y9vU9wQ4Wzf6ka9MV8zxwGskbQfGgBMkfb4nqczMrLKOC3tE/JeIODQihoDXAVdGxGk9S2ZmZpX4PHYzs8xU6qiNiHFgvNYkZmZWC2+xm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMdFzYJe0r6duSvivpZkkf6mUwMzOrpswvKD0MnBARU5KWAd+Q9E8R8S89ymZmZhV0XNgjIoCp4u6y4i96EcrMzKpTqtcdjiwtAbYCzwbOi4g/mWWcUWAUYHBwcO3Y2Nis05qammJgYGDOeU3u2NVxrroNLoedu/s2+wU1OV/ZbGtWr+xdmBna17l+rV/zLe9C74l+eqJma8LrXLbtRkZGtkZEq5v5lyrsv3iSdCBwGfCuiLhprvFarVZMTEzMOmx8fJzh4eE55zG0YUvpXHVZv2YPmyYr/c73omhyvrLZtm9c18M0/177Otev9Wu+5V3oPdFPT9RsTXidy7adpK4Le6WzYiLifmAcOLGbmZuZWf3KnBVzULGljqTlwMuBH/QqmJmZVVNmf/4pwAVFP/tewBcj4vLexDIzs6rKnBXzPeCYHmYxM7Ma+MpTM7PMuLCbmWXGhd3MLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8yU+c3Tp0m6StI2STdLOquXwczMrJoyv3m6B1gfETdI2h/YKulrEfH9HmUzM7MKOt5ij4h/i4gbitsPANuA1b0KZmZm1Sgiyj9JGgKuAY6OiJ/NGDYKjAIMDg6uHRsbm3UaU1NTDAwMzDmPyR27Sueqy+By2Lm7b7NfUJPzOVt1vc63ZvXKys9d6P3aT/Nl61cdaW/rsm03MjKyNSJa3cy/dGGXNABcDZwTEZfON26r1YqJiYlZh42PjzM8PDznc4c2bCmVq07r1+xh02SZXqrF1eR8zlZdr/Nt37iu8nMXer/203zZ+lVH2tu6bNtJ6rqwlzorRtIy4BLgCwsVdTMz648yZ8UI+DSwLSL+tHeRzMysG2W22I8H3gicIOk7xd+repTLzMwq6rhDLyK+AaiHWczMrAa+8tTMLDMu7GZmmXFhNzPLjAu7mVlmXNjNzDLjwm5mlhkXdjOzzLiwm5llxoXdzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsM2V+8/Qzku6WdFMvA5mZWXfKbLGfD5zYoxxmZlaTjgt7RFwD3NfDLGZmVgNFROcjS0PA5RFx9DzjjAKjAIODg2vHxsZmHW9qaoqBgYE55zW5Y1fHueo2uBx27u7b7BfU5HzOVl2v861ZvbLycxd6v/bTfNn6VUfa27ps242MjGyNiFY386+9sLdrtVoxMTEx67Dx8XGGh4fnfO7Qhi0d56rb+jV72DS5tG/zX0iT8zlbdb3Ot33jusrPXej92k/zZetXHWlv67JtJ6nrwu6zYszMMuPCbmaWmTKnO14EfAs4QtIdkt7au1hmZlZVxx16EfH6XgYxM7N6uCvGzCwzLuxmZplxYTczy4wLu5lZZlzYzcwy48JuZpYZF3Yzs8y4sJuZZcaF3cwsMy7sZmaZcWE3M8uMC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDOlCrukEyXdIumHkjb0KpSZmVVX5jdPlwDnAa8EjgJeL+moXgUzM7NqymyxvxD4YUTcFhGPAGPAb/cmlpmZVaWI6GxE6WTgxIg4o7j/RuBFEfHOGeONAqPF3SOAW+aY5CrgniqhF0GTs0Gz8zlbdU3O52zVlc13WEQc1M0Ml5YYV7M89kufChGxGdi84MSkiYholZj/omlyNmh2Pmerrsn5nK26fuQr0xVzB/C0tvuHAnfWG8fMzLpVprBfDxwu6RmS9gZeB/xjb2KZmVlVHXfFRMQeSe8EvgIsAT4TETd3Me8Fu2v6qMnZoNn5nK26JudztuoWPV/HB0/NzOyJwVeempllxoXdzCwzLuxmZpl5QhZ2Saslre53jtlIeqak90g6od9ZZmpyNmh2Pmerrsn5mpwNqud7QhV2SUOSrgauAD4u6SX9ztRO0n8Evkb6Lp23SXp7nyP9QpOzQbPzOVt1Tc7X5GzQZb6IaPQfsG/b7dcCnyhuvxn4B2BNcV99yHYC8Izp+QMfAE4r7r8I+BIw3I98Tc7W9HzOlme+JmerO18jt9glHSDpryTdCnxC0mHFoN8BflTcHgN+CJwx/bRFzHeUpO8B/w34rKQTIrX2UcAhABFxHfBN4C2Lma/J2Zqez9nyzNfkbL3K18jCDpwI7EtasEeAD0haTtoteTVARDwMXAy8pLj/WK/CSDpU0gFtD50CXBIRLyV9wLxB0uHAhdP5CpcBR0vap1f5mpyt6fmcLc98Tc62WPn6VtiVLJX0Vklfl3SWpGcVg58NPBIRe4A/A34KnAZ8FXiKpF8rxrsV+LGk43qU8UhJXwa+AXxY0vTXFP8/YL/i9heBu4B1pE/UJ7ftYdxH+nbL5/0qZWt6PmfLM1+Tsy12vr4V9mJX4zeBNwEfA/YB/qYYfBdwd/HJ9GPSwjyL1ADf5/GvBV4G3Fs8XgtJK9ruPh+4IyKGgCuBTxSP3wc8LGn/iLgP+FfgqUWObwLvLcbbG/g5sD33bE3P52x55mtytn7mW7TCLuk4SR+VdHpxX8CRwBUR8aWI+BhwmKQXAztIn2BHFk/fBgwUj/0F8CpJryZ9KAwC3+0y25MknS/pemCjpIOKfGuAayUpIv4RuF/SOtKewv7FcIr7BwOPkfYwDpb0N8BFwJ6IuDvHbE3P52zVNTlfk7M1Jd+iFHZJzwH+EngA+D1J7y3mvRp4oFhogPOBN5AK9R7gxcXjN5COGD8UEdcAG4DTgeOB/x4Rj7VNo4qXFvN7FemgxNnAAaQvOzuk2LsAuKDI9+1iWV4JEBHfKqaxNCK2AWcCNwP/MyLeQneanK3p+Zwtz3xNztaMfHOdLlP1j7RlfQZpt2Np8difAmcVt1vAJ4GTgZcDX2l77tNIuyqQCvmNpF9hOgb438BT2sYtfTpS0bBnAleTunNWFY9/EXh3cfsZwMZi+LGk/rAlbcv2k2I6q0l7Eu8EPgt8CljRRbs1NlvT8zmbX1e33b//q3WLXdLzSQc4fxv4IPC+YtAO0m+mQvrkuRb4XeCfgUMkPVfSskj96TskvSQiriR93eVHgUuBiyLi36bnFUXLlHQS8BrgQ8BxpL59SGfbTO8d/Bj4OvDKiLie9Ik7UsxzCrgOODYidgBvJHUF3QW8LyIeLBuobU/j1U3LNoPbrprGtRu47brJ9kRouzI/jfdLJL0QOBz4akT8hLQ1fmtEnC7pBcA5klrAOPCfJO0XEQ9J+i7we6RzNC8E/gD4pKTdwCRwezGLvwIujIhdJTIpIkLSsaTdnK8DWyKdHvnrwG0RcaWk20lXr74C2Ar8jqRVEXGPpH8FHpT0dODPgdMkHUz61ah7SbtORMQEMFGh3VqkvZoHgI8DdwPP7Hc2t121bE+EdnPb5dd28ym1xa5kmaQ3SbqR1LF/IDBdeH8ObC+2vm8g7VocBzzE46fwADxK2gU5hLRVfhOpf/1q4J6IuAPSVnnFov5S4DOko8ovBz5SjPIYcKuk5RFxe5HvuaQX607S+aTTy7GE1D6XFBlPBdYCm6PiOa6SVkr6bDHN24FzI+JuSXuRPsn7mW1J0Xa/SdoVbEzbFevdgKTzaVjbFfMMScM0c53bR9KKhrbdAQ1vuwFJ+0q6gIa13YI66a8BVgAvLm4fWAT75CzjnUW6DHZ1cf9kUn/6YaSvALi6eHxfUjfMqrbnHgPs3UmeGfPcD3gbj2/5LwP+CHhHMfxJwPeK6Z9C6u8aKoadVCzLquL2JLCS1L//5fY8wF5dZLuIdMXYAKlr6cy2caaPQ7wT+B+Lla3tdT2DtLKtJx3gaUrbTWe7tFivDmpY2+0PbCH9khjAe5rQbjPyfRn46+L+x4C39bvtSO+JN5Pe/5c0re3a8l0J/H3xWGPWu07/Ftxil3Q2cBuwRdJgRNxP6he6s+gbf40ev0DoW6QDoNMXGl1LOoj6UERcAPxU0udIB0VvAX7RhxQRN0bEIwvlmZHtEOByYBj4HOkAxWtJewl7iun+lHTg9d2kvq+Defw0ymtI59I/EhGXA58mXc16HumI9aNt+Up9qs7I9nfA24tstwJHSNpYbEX9vtIFV1eQ9mB6nq3It4L05jqBdP3AK0jHPY4lbSn1s+3as20mnS3wWtI1DL/R77YrLCdde/EsSatI6/ySYpp9abdZ8u1NWteeSuriOFrSR/rVdpKWkY6xnQx8PCJ+txh0TNs0+9Z2M/J9LCKmt7gngaP62XaldfAJNkzavfhb4D3FY8eSitYdRfALgU3FsHOAD7c9/3rgmOL2PqRTgI6t41OJtPK+qO3+6aQtkzcD3257/KnAncXtd5Au231S8fwvAU9vG3dVj7K9iXSk+9eBvy/+Xg/8L9K5/IuWrW16B7bd/s+kN9Op/W67WbL9MemUsWc2qO3eTOprfT/wVtKBtOv73W6z5HsfaY9nVRPajrQHduqMx04BrmtC282R7+lFhr6vdx0vRwcLOn1qzinAeHF7GWlramVx/zDS1vqxpF3Ai0lbWv9E+qTapyfhUx+X4Be/3foCHu/uuZd0zuj0uF+jKLSk3aevFuP8ySJlOwb4xvSK2zbeMtLB5ROK++f0OtuMnAeQjm/sBD5c3L8XGOxX282S7a5ivisouvn61XZtr+dbSN1srwW+UDx2T7/bbY58Y8Vj7acL92W9I3VR3ApsKub/gaJ+3Acc3IB1rj3fVaQv5jq03+td6eUoscBPJl0o9Jzi/tIZw88HTp5egUhdD2fSo6I+x8p8AY+fL/854KPF7V8j7XE8ve2FOZq2rwRepGzvaH+suH1I0XbPXexsbRn+kHS+7WZSv/Y3izec+tl2M7KdRzqt7NlNaDvSV0YvIfWhXk3aMr4JeH+/17lZ8v0z6QyzFzSk7b5C2gN7Gmkr+CzShmFT1rn2fF8gXfp/eBPartO/6aLTEUmfAn4WERuK+3uRzrt8B/Ac4JQo2U9eF0mHkvq03hURtyp9odhokWs18J2o56qybrK9PSJuKx47htQtta7I9of9yNZO6TqEM0lvsiNJK+uh9LHt2rIdTXqz/TnpLKuT6FPbSRogdXPsQ2qn3yBdeHI2aUv5cPrYbrPkO5x0fOK3SMe8XkZqv76sdypOey5uP4/0Pr2WdEl939e5GfmOJl3pfi7pm2b7tt6VUfY89s3AucVBhiNJK/HxpBfl7H4V9cIxFOfASzqD1P9/NqkL6QeRTr/sd7YfFdluJ60ce0hb8Tf2MVu7e0kHAd8XEX8n6TTg5obku5/UT3wT6XVdRv/abg/p7IlHSVvqPyet/5PAexvQbnPle1jSa0gFv2/r3XTRLNxPOu70/oi4sAFtNzPfA6SN123Af6W/613Hym6xv450oPRh0jeOXRkRt/QoWymSriUdXNtOOof0QxHxvb6GKszIdhewoUHttpK0BfcG0vffbwbOi4hH533iIpgl26cjYlN/U/2y4sKT6b7su/qdZ6Yi38nAZyOdddLvPPuQfnPhjaQ96r8EPhXpa7r7bpZ8myPiz/qbqpyOC7uk55LO57yYdLCotq/K7VaxB/FB0pbw5yNdtdYITc4GIGkpqfvlYVK+Jr2ujc0G6aIu4LEos3W0iJqcT9KZpNNqP9e01xWan28hpbbYzcys+Zr603hmZlaRC7uZWWZc2M3MMuPCbmaWGRd2M7PMuLCbmWXGhd3MLDP/H+KofDj+oV4qAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.hist(xrot=20)" ] @@ -341,9 +2479,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": true - }, + "metadata": {}, "outputs": [], "source": [] } @@ -364,7 +2500,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4, -- 2.18.1