From 195566114437af3cba6fb0db0b3f342f74120072 Mon Sep 17 00:00:00 2001 From: c84ad4823ff7a24ecd9aadc23649b02f Date: Fri, 18 Aug 2023 11:24:08 +0000 Subject: [PATCH] Update analyse-syndrome-grippal.ipynb --- module3/exo1/analyse-syndrome-grippal.ipynb | 2252 ++++++++++++++++++- 1 file changed, 2214 insertions(+), 38 deletions(-) diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5..87f3451 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,15 +28,34 @@ }, { "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\"" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Les données sont téléchargées et sauvegardées en locale pour qu’elles soient pérennes." + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + "data_file = \"syndrome-grippal.csv\"\n", + "\n", + "import os\n", + "import urllib.request\n", + "if not os.path.exists(data_file):\n", + " urllib.request.urlretrieve(data_url, data_file)" + ] + }, { "cell_type": "markdown", "metadata": {}, @@ -61,11 +80,978 @@ }, { "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", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \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
020233231678411210.022358.02517.033.0FRFrance
120233131543610820.020052.02316.030.0FRFrance
22023303132058647.017763.02013.027.0FRFrance
32023293111227113.015131.01711.023.0FRFrance
4202328391795703.012655.0149.019.0FRFrance
5202327389995763.012235.0149.019.0FRFrance
6202326390235934.012112.0149.019.0FRFrance
72023253100906739.013441.01510.020.0FRFrance
82023243113087639.014977.01711.023.0FRFrance
920232331430010661.017939.02217.027.0FRFrance
1020232231830313822.022784.02821.035.0FRFrance
1120232131646012188.020732.02519.031.0FRFrance
1220232031616211963.020361.02418.030.0FRFrance
1320231931690112577.021225.02518.032.0FRFrance
1420231831992915402.024456.03023.037.0FRFrance
1520231732700721779.032235.04133.049.0FRFrance
1620231632787522767.032983.04234.050.0FRFrance
1720231533745530993.043917.05646.066.0FRFrance
1820231434806040671.055449.07261.083.0FRFrance
1920231336485956800.072918.09886.0110.0FRFrance
2020231237275064499.081001.010997.0121.0FRFrance
2120231137463866420.082856.0112100.0124.0FRFrance
2220231037636868243.084493.0115103.0127.0FRFrance
2320230936206254778.069346.09382.0104.0FRFrance
2420230837639168065.084717.0115102.0128.0FRFrance
2520230738985180397.099305.0135121.0149.0FRFrance
2620230639736887636.0107100.0146131.0161.0FRFrance
2720230539546986268.0104670.0144130.0158.0FRFrance
2820230437490166916.082886.0113101.0125.0FRFrance
2920230336957061893.077247.010593.0117.0FRFrance
.................................
199419852132609619621.032571.04735.059.0FRFrance
199519852032789620885.034907.05138.064.0FRFrance
199619851934315432821.053487.07859.097.0FRFrance
199719851834055529935.051175.07455.093.0FRFrance
199819851733405324366.043740.06244.080.0FRFrance
199919851635036236451.064273.09166.0116.0FRFrance
200019851536388145538.082224.011683.0149.0FRFrance
20011985143134545114400.0154690.0244207.0281.0FRFrance
20021985133197206176080.0218332.0357319.0395.0FRFrance
20031985123245240223304.0267176.0445405.0485.0FRFrance
20041985113276205252399.0300011.0501458.0544.0FRFrance
20051985103353231326279.0380183.0640591.0689.0FRFrance
20061985093369895341109.0398681.0670618.0722.0FRFrance
20071985083389886359529.0420243.0707652.0762.0FRFrance
20081985073471852432599.0511105.0855784.0926.0FRFrance
20091985063565825518011.0613639.01026939.01113.0FRFrance
20101985053637302592795.0681809.011551074.01236.0FRFrance
20111985043424937390794.0459080.0770708.0832.0FRFrance
20121985033213901174689.0253113.0388317.0459.0FRFrance
201319850239758680949.0114223.0177147.0207.0FRFrance
201419850138548965918.0105060.0155120.0190.0FRFrance
201519845238483060602.0109058.0154110.0198.0FRFrance
2016198451310172680242.0123210.0185146.0224.0FRFrance
20171984503123680101401.0145959.0225184.0266.0FRFrance
2018198449310107381684.0120462.0184149.0219.0FRFrance
201919844837862060634.096606.0143110.0176.0FRFrance
202019844737202954274.089784.013199.0163.0FRFrance
202119844638733067686.0106974.0159123.0195.0FRFrance
20221984453135223101414.0169032.0246184.0308.0FRFrance
202319844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2024 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202332 3 16784 11210.0 22358.0 25 17.0 \n", + "1 202331 3 15436 10820.0 20052.0 23 16.0 \n", + "2 202330 3 13205 8647.0 17763.0 20 13.0 \n", + "3 202329 3 11122 7113.0 15131.0 17 11.0 \n", + "4 202328 3 9179 5703.0 12655.0 14 9.0 \n", + "5 202327 3 8999 5763.0 12235.0 14 9.0 \n", + "6 202326 3 9023 5934.0 12112.0 14 9.0 \n", + "7 202325 3 10090 6739.0 13441.0 15 10.0 \n", + "8 202324 3 11308 7639.0 14977.0 17 11.0 \n", + "9 202323 3 14300 10661.0 17939.0 22 17.0 \n", + "10 202322 3 18303 13822.0 22784.0 28 21.0 \n", + "11 202321 3 16460 12188.0 20732.0 25 19.0 \n", + "12 202320 3 16162 11963.0 20361.0 24 18.0 \n", + "13 202319 3 16901 12577.0 21225.0 25 18.0 \n", + "14 202318 3 19929 15402.0 24456.0 30 23.0 \n", + "15 202317 3 27007 21779.0 32235.0 41 33.0 \n", + "16 202316 3 27875 22767.0 32983.0 42 34.0 \n", + "17 202315 3 37455 30993.0 43917.0 56 46.0 \n", + "18 202314 3 48060 40671.0 55449.0 72 61.0 \n", + "19 202313 3 64859 56800.0 72918.0 98 86.0 \n", + "20 202312 3 72750 64499.0 81001.0 109 97.0 \n", + "21 202311 3 74638 66420.0 82856.0 112 100.0 \n", + "22 202310 3 76368 68243.0 84493.0 115 103.0 \n", + "23 202309 3 62062 54778.0 69346.0 93 82.0 \n", + "24 202308 3 76391 68065.0 84717.0 115 102.0 \n", + "25 202307 3 89851 80397.0 99305.0 135 121.0 \n", + "26 202306 3 97368 87636.0 107100.0 146 131.0 \n", + "27 202305 3 95469 86268.0 104670.0 144 130.0 \n", + "28 202304 3 74901 66916.0 82886.0 113 101.0 \n", + "29 202303 3 69570 61893.0 77247.0 105 93.0 \n", + "... ... ... ... ... ... ... ... \n", + "1994 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1995 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1996 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1997 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1998 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1999 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2000 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2001 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2002 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2003 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2004 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2005 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2006 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2007 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2008 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2009 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2010 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2011 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2012 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2013 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2014 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2015 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2016 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2017 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2018 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2019 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2020 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2021 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2022 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2023 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 33.0 FR France \n", + "1 30.0 FR France \n", + "2 27.0 FR France \n", + "3 23.0 FR France \n", + "4 19.0 FR France \n", + "5 19.0 FR France \n", + "6 19.0 FR France \n", + "7 20.0 FR France \n", + "8 23.0 FR France \n", + "9 27.0 FR France \n", + "10 35.0 FR France \n", + "11 31.0 FR France \n", + "12 30.0 FR France \n", + "13 32.0 FR France \n", + "14 37.0 FR France \n", + "15 49.0 FR France \n", + "16 50.0 FR France \n", + "17 66.0 FR France \n", + "18 83.0 FR France \n", + "19 110.0 FR France \n", + "20 121.0 FR France \n", + "21 124.0 FR France \n", + "22 127.0 FR France \n", + "23 104.0 FR France \n", + "24 128.0 FR France \n", + "25 149.0 FR France \n", + "26 161.0 FR France \n", + "27 158.0 FR France \n", + "28 125.0 FR France \n", + "29 117.0 FR France \n", + "... ... ... ... \n", + "1994 59.0 FR France \n", + "1995 64.0 FR France \n", + "1996 97.0 FR France \n", + "1997 93.0 FR France \n", + "1998 80.0 FR France \n", + "1999 116.0 FR France \n", + "2000 149.0 FR France \n", + "2001 281.0 FR France \n", + "2002 395.0 FR France \n", + "2003 485.0 FR France \n", + "2004 544.0 FR France \n", + "2005 689.0 FR France \n", + "2006 722.0 FR France \n", + "2007 762.0 FR France \n", + "2008 926.0 FR France \n", + "2009 1113.0 FR France \n", + "2010 1236.0 FR France \n", + "2011 832.0 FR France \n", + "2012 459.0 FR France \n", + "2013 207.0 FR France \n", + "2014 190.0 FR France \n", + "2015 198.0 FR France \n", + "2016 224.0 FR France \n", + "2017 266.0 FR France \n", + "2018 219.0 FR France \n", + "2019 176.0 FR France \n", + "2020 163.0 FR France \n", + "2021 195.0 FR France \n", + "2022 308.0 FR France \n", + "2023 213.0 FR France \n", + "\n", + "[2024 rows x 10 columns]" + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "raw_data = pd.read_csv(data_url, skiprows=1)\n", + "raw_data = pd.read_csv(data_file, encoding = 'iso-8859-1', skiprows=1)\n", "raw_data" ] }, @@ -78,9 +1064,73 @@ }, { "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", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
17871989193-NaNNaN-NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1787 198919 3 - NaN NaN - NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1787 FR France " + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -94,9 +1144,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 6, "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
020233231678411210.022358.02517.033.0FRFrance
120233131543610820.020052.02316.030.0FRFrance
22023303132058647.017763.02013.027.0FRFrance
32023293111227113.015131.01711.023.0FRFrance
4202328391795703.012655.0149.019.0FRFrance
5202327389995763.012235.0149.019.0FRFrance
6202326390235934.012112.0149.019.0FRFrance
72023253100906739.013441.01510.020.0FRFrance
82023243113087639.014977.01711.023.0FRFrance
920232331430010661.017939.02217.027.0FRFrance
1020232231830313822.022784.02821.035.0FRFrance
1120232131646012188.020732.02519.031.0FRFrance
1220232031616211963.020361.02418.030.0FRFrance
1320231931690112577.021225.02518.032.0FRFrance
1420231831992915402.024456.03023.037.0FRFrance
1520231732700721779.032235.04133.049.0FRFrance
1620231632787522767.032983.04234.050.0FRFrance
1720231533745530993.043917.05646.066.0FRFrance
1820231434806040671.055449.07261.083.0FRFrance
1920231336485956800.072918.09886.0110.0FRFrance
2020231237275064499.081001.010997.0121.0FRFrance
2120231137463866420.082856.0112100.0124.0FRFrance
2220231037636868243.084493.0115103.0127.0FRFrance
2320230936206254778.069346.09382.0104.0FRFrance
2420230837639168065.084717.0115102.0128.0FRFrance
2520230738985180397.099305.0135121.0149.0FRFrance
2620230639736887636.0107100.0146131.0161.0FRFrance
2720230539546986268.0104670.0144130.0158.0FRFrance
2820230437490166916.082886.0113101.0125.0FRFrance
2920230336957061893.077247.010593.0117.0FRFrance
.................................
199419852132609619621.032571.04735.059.0FRFrance
199519852032789620885.034907.05138.064.0FRFrance
199619851934315432821.053487.07859.097.0FRFrance
199719851834055529935.051175.07455.093.0FRFrance
199819851733405324366.043740.06244.080.0FRFrance
199919851635036236451.064273.09166.0116.0FRFrance
200019851536388145538.082224.011683.0149.0FRFrance
20011985143134545114400.0154690.0244207.0281.0FRFrance
20021985133197206176080.0218332.0357319.0395.0FRFrance
20031985123245240223304.0267176.0445405.0485.0FRFrance
20041985113276205252399.0300011.0501458.0544.0FRFrance
20051985103353231326279.0380183.0640591.0689.0FRFrance
20061985093369895341109.0398681.0670618.0722.0FRFrance
20071985083389886359529.0420243.0707652.0762.0FRFrance
20081985073471852432599.0511105.0855784.0926.0FRFrance
20091985063565825518011.0613639.01026939.01113.0FRFrance
20101985053637302592795.0681809.011551074.01236.0FRFrance
20111985043424937390794.0459080.0770708.0832.0FRFrance
20121985033213901174689.0253113.0388317.0459.0FRFrance
201319850239758680949.0114223.0177147.0207.0FRFrance
201419850138548965918.0105060.0155120.0190.0FRFrance
201519845238483060602.0109058.0154110.0198.0FRFrance
2016198451310172680242.0123210.0185146.0224.0FRFrance
20171984503123680101401.0145959.0225184.0266.0FRFrance
2018198449310107381684.0120462.0184149.0219.0FRFrance
201919844837862060634.096606.0143110.0176.0FRFrance
202019844737202954274.089784.013199.0163.0FRFrance
202119844638733067686.0106974.0159123.0195.0FRFrance
20221984453135223101414.0169032.0246184.0308.0FRFrance
202319844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2023 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202332 3 16784 11210.0 22358.0 25 17.0 \n", + "1 202331 3 15436 10820.0 20052.0 23 16.0 \n", + "2 202330 3 13205 8647.0 17763.0 20 13.0 \n", + "3 202329 3 11122 7113.0 15131.0 17 11.0 \n", + "4 202328 3 9179 5703.0 12655.0 14 9.0 \n", + "5 202327 3 8999 5763.0 12235.0 14 9.0 \n", + "6 202326 3 9023 5934.0 12112.0 14 9.0 \n", + "7 202325 3 10090 6739.0 13441.0 15 10.0 \n", + "8 202324 3 11308 7639.0 14977.0 17 11.0 \n", + "9 202323 3 14300 10661.0 17939.0 22 17.0 \n", + "10 202322 3 18303 13822.0 22784.0 28 21.0 \n", + "11 202321 3 16460 12188.0 20732.0 25 19.0 \n", + "12 202320 3 16162 11963.0 20361.0 24 18.0 \n", + "13 202319 3 16901 12577.0 21225.0 25 18.0 \n", + "14 202318 3 19929 15402.0 24456.0 30 23.0 \n", + "15 202317 3 27007 21779.0 32235.0 41 33.0 \n", + "16 202316 3 27875 22767.0 32983.0 42 34.0 \n", + "17 202315 3 37455 30993.0 43917.0 56 46.0 \n", + "18 202314 3 48060 40671.0 55449.0 72 61.0 \n", + "19 202313 3 64859 56800.0 72918.0 98 86.0 \n", + "20 202312 3 72750 64499.0 81001.0 109 97.0 \n", + "21 202311 3 74638 66420.0 82856.0 112 100.0 \n", + "22 202310 3 76368 68243.0 84493.0 115 103.0 \n", + "23 202309 3 62062 54778.0 69346.0 93 82.0 \n", + "24 202308 3 76391 68065.0 84717.0 115 102.0 \n", + "25 202307 3 89851 80397.0 99305.0 135 121.0 \n", + "26 202306 3 97368 87636.0 107100.0 146 131.0 \n", + "27 202305 3 95469 86268.0 104670.0 144 130.0 \n", + "28 202304 3 74901 66916.0 82886.0 113 101.0 \n", + "29 202303 3 69570 61893.0 77247.0 105 93.0 \n", + "... ... ... ... ... ... ... ... \n", + "1994 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1995 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1996 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1997 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1998 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1999 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2000 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2001 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2002 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2003 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2004 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2005 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2006 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2007 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2008 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2009 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2010 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2011 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2012 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2013 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2014 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2015 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2016 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2017 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2018 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2019 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2020 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2021 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2022 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2023 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 33.0 FR France \n", + "1 30.0 FR France \n", + "2 27.0 FR France \n", + "3 23.0 FR France \n", + "4 19.0 FR France \n", + "5 19.0 FR France \n", + "6 19.0 FR France \n", + "7 20.0 FR France \n", + "8 23.0 FR France \n", + "9 27.0 FR France \n", + "10 35.0 FR France \n", + "11 31.0 FR France \n", + "12 30.0 FR France \n", + "13 32.0 FR France \n", + "14 37.0 FR France \n", + "15 49.0 FR France \n", + "16 50.0 FR France \n", + "17 66.0 FR France \n", + "18 83.0 FR France \n", + "19 110.0 FR France \n", + "20 121.0 FR France \n", + "21 124.0 FR France \n", + "22 127.0 FR France \n", + "23 104.0 FR France \n", + "24 128.0 FR France \n", + "25 149.0 FR France \n", + "26 161.0 FR France \n", + "27 158.0 FR France \n", + "28 125.0 FR France \n", + "29 117.0 FR France \n", + "... ... ... ... \n", + "1994 59.0 FR France \n", + "1995 64.0 FR France \n", + "1996 97.0 FR France \n", + "1997 93.0 FR France \n", + "1998 80.0 FR France \n", + "1999 116.0 FR France \n", + "2000 149.0 FR France \n", + "2001 281.0 FR France \n", + "2002 395.0 FR France \n", + "2003 485.0 FR France \n", + "2004 544.0 FR France \n", + "2005 689.0 FR France \n", + "2006 722.0 FR France \n", + "2007 762.0 FR France \n", + "2008 926.0 FR France \n", + "2009 1113.0 FR France \n", + "2010 1236.0 FR France \n", + "2011 832.0 FR France \n", + "2012 459.0 FR France \n", + "2013 207.0 FR France \n", + "2014 190.0 FR France \n", + "2015 198.0 FR France \n", + "2016 224.0 FR France \n", + "2017 266.0 FR France \n", + "2018 219.0 FR France \n", + "2019 176.0 FR France \n", + "2020 163.0 FR France \n", + "2021 195.0 FR France \n", + "2022 308.0 FR France \n", + "2023 213.0 FR France \n", + "\n", + "[2023 rows x 10 columns]" + ] + }, + "execution_count": 6, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -122,7 +2139,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 7, "metadata": {}, "outputs": [], "source": [ @@ -152,10 +2169,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 8, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -179,9 +2194,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n" + ] + } + ], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", @@ -194,15 +2217,55 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "Un premier regard sur les données !" + "Un premier regard sur les données ! La commande suivante ne marche pas parce que les données sont des chaînes de caractères." ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "metadata": {}, "outputs": [], "source": [ + "# sorted_data['inc'].plot()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Les données doivent être converti en nombre éntiers :" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZkAAAEKCAYAAADAVygjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXm8HUWZ93/PWe6efSMkIWGJ7LJFFkVlUYKjA7woIzpKVGZ4hxd3xxHmHWUGZMR5R0EGQRlBFhdExAFFwBCQNQQCAUISyL4v9yY3y11yl3NOvX90VZ/uPlVdffrsN8/388kn51Z3V1VXd9dTz1JVJIQAwzAMw1SCRK0rwDAMw4xcWMgwDMMwFYOFDMMwDFMxWMgwDMMwFYOFDMMwDFMxWMgwDMMwFYOFDMMwDFMxWMgwDMMwFYOFDMMwDFMxUrWuQDWZOHGimDVrVq2rwTAM01C8+uqrO4UQk+Jce0AJmVmzZmHx4sW1rgbDMExDQUQb4l7L5jKGYRimYrCQYRiGYSoGCxmGYRimYrCQYRiGYSoGCxmGYRimYliFDBHdRUSdRPSWJ208Ec0nolXy/3GeY9cQ0WoieoeI5nrSTyGipfLYLUREMr2ZiH4j0xcR0SzPNfNkGauIaJ4n/VB57ip5bVPpTcEwDMOUmyiazN0Azg+kXQ1ggRBiNoAF8m8Q0TEALgVwrLzmNiJKymtuB3AFgNnyn8rzcgC7hRBHALgJwPdlXuMBXAvgNACnArjWI8y+D+AmWf5umQfDMAxTZ1iFjBDiWQDdgeQLAdwjf98D4CJP+v1CiEEhxDoAqwGcSkRTAYwWQiwUzn7P9wauUXk9COBcqeXMBTBfCNEthNgNYD6A8+Wxc+S5wfJHJP1DGTz02uZaV4NhGKZo4k7GnCKE2AYAQohtRDRZpk8D8JLnvM0ybVj+DqarazbJvDJEtBfABG964JoJAPYIITKavEYk33l4GR58dTNmTmjHKTPH2S9gGIapE8rt+CdNmghJj3NNWF6FFSK6gogWE9Hirq4u02l1zZbd+wEAg8PZGteEYRimOOIKmR3SBAb5f6dM3wxghue86QC2yvTpmnTfNUSUAjAGjnnOlNdOAGPlucG8ChBC3CGEmCOEmDNpUqyld2pONufI0GRCJ18ZhmHql7hC5hEAKtprHoCHPemXyoixQ+E4+F+WprUeIjpd+lQuC1yj8voEgKek3+YJAOcR0Tjp8D8PwBPy2NPy3GD5I5LhXA4AkEqykGEYprGw+mSI6NcAzgIwkYg2w4n4uhHAA0R0OYCNAC4BACHEMiJ6AMByABkAVwkhlI3nSjiRaq0AHpP/AOBOAPcR0Wo4GsylMq9uIroewCvyvOuEECoA4VsA7iei7wJYIvMYseQ1GZ7WxDBMY2EVMkKITxkOnWs4/wYAN2jSFwM4TpM+ACmkNMfuAnCXJn0tnLDmA4JM1hEyKTaXMQzTYPDQuAFgnwzDMI0KC5kGICccISOMMXQMwzD1CQuZBoCkAiPMkdoMwzB1CQuZBoDk1CDWZBiGaTRYyDQQLGQYhmk0WMg0AMpclmMpwzBMg8FCpoFgEcMwTKPBQqYBkFvvsCbDMEzDwUKmAVCzY1jGMAzTaLCQaQDcEGaWMgzDNBgsZBqA/DwZhmGYxoKFTAOg5snkcixmGIZpLFjINACsyTAM06iwkGkA2PHPMEyjwkKmESC1rEz9SZllW/di1tWP4s3Ne2pdFYZh6hAWMg1E/YkYYMEKZ+ftPy/bUeOaMAxTj7CQaQCUuaweJ2PyDjcMw4TBQqYByM+TqW09GIZhioWFTANQz5oMwzBMGCxkGgCi+t9PhjdUYxhGBwuZBqCe/R5Uz5VjGKbmsJBpIOpZW6hnLYthmNrBQoYpCWJVhmGYEFjINBD1rC3UcdUYhqkhLGQaiHoWMgzDMDpYyDAMwzAVg4VMA8GKDMMwjQYLGYZhGKZisJBpIOpxFWZFHVeNYZgawkJmBPDsyi4c/e3H0TMwXPWyOYKZYZgwShIyRPQ1IlpGRG8R0a+JqIWIxhPRfCJaJf8f5zn/GiJaTUTvENFcT/opRLRUHruF5OQLImomot/I9EVENMtzzTxZxioimlfKfTQKJmXhpidXYv9wFit39FS1Pl7qeaIowzC1I7aQIaJpAL4MYI4Q4jgASQCXArgawAIhxGwAC+TfIKJj5PFjAZwP4DYiSsrsbgdwBYDZ8t/5Mv1yALuFEEcAuAnA92Ve4wFcC+A0AKcCuNYrzEYqJpNULZUJqutFbxiGqTWlmstSAFqJKAWgDcBWABcCuEcevwfARfL3hQDuF0IMCiHWAVgN4FQimgpgtBBioXCcDvcGrlF5PQjgXKnlzAUwXwjRLYTYDWA+8oKJqQWsyDAMoyG2kBFCbAHwnwA2AtgGYK8Q4s8ApgghtslztgGYLC+ZBmCTJ4vNMm2a/B1M910jhMgA2AtgQkheI5z668nZJ8MwTBilmMvGwdE0DgVwMIB2IvpM2CWaNBGSHveaYD2vIKLFRLS4q6srpHoMwzBMuSnFXPYhAOuEEF1CiGEADwF4L4Ad0gQG+X+nPH8zgBme66fDMa9tlr+D6b5rpEluDIDukLwKEELcIYSYI4SYM2nSpJi3Wh/YwoQ5jJhhmHqjFCGzEcDpRNQm/STnAlgB4BEAKtprHoCH5e9HAFwqI8YOhePgf1ma1HqI6HSZz2WBa1RenwDwlPTbPAHgPCIaJzWq82TaiMTdftl4nEKPVwOWbwzD6EjFvVAIsYiIHgTwGoAMgCUA7gDQAeABIrocjiC6RJ6/jIgeALBcnn+VECIrs7sSwN0AWgE8Jv8BwJ0A7iOi1XA0mEtlXt1EdD2AV+R51wkhuuPeS6NT2+gyhmEYM7GFDAAIIa6FE0rsZRCOVqM7/wYAN2jSFwM4TpM+ACmkNMfuAnBXkVVuaOrZHFbPqxEwDFM7eMZ/A1DP/TdHlzEMEwYLmQaCZ9UzDNNosJAZQdRS46lnbYthmNrBQqaBMC4ro6LPatDT87IyDMOEwUKmAbCGMDdwR791z3588qcLsbe/+itIMwxTeVjIMDXl1qdXY9G6bjzypnYuLcMwDQ4LmQainsOES61Z4+piDMOEwUKGKQkOYWYYJgwWMiOImi4rU79KFsMwNYSFTANh7MjJcryOacQ6MwwTHRYyI4B6sFiVOlGUzW4MMzJhIdNA1OOMf2LpwDBMCCxkmBpTf4KTYZjywUKmgRjJ/otGnlDKMIwZFjJMWRjJApBhmPiwkGkgrGuX1cD0xPoHwzBhsJBpIKxrlzWgNqEEJ8cPMMzIhIUMU1NcIVPbajAMUyFYyDQQ9bh2GWsgDMOEwUKGKQv1KAAZhqk9LGQaiHrsxktVZOpxginDMOWDhUwDYHPs2zY1awTY7MYwIxMWMiMAauAFMhmGGdmwkGkg6tm0VL81YximlrCQaQDqWbiUukAma18MM7JhIdNA2DrkWgqjUoWFae2y4WwOnfsGSsucYZiawUKmCBau2YWFa3bVuhpGaqEVVNph/63fvYlT/30BBjPZyhbEMExFSNW6Ao3Ep/77JQDA+hs/WpPybTKkptsvV6j0x9/aDgAYzgo089vKMA0HazINgDIlGRfIdI833gKZbo0NGfGyMwzT2LCQGQE0cgizTYgoDYnn0TBMY8JCpoGwmaTqOQotLo0oOBmGyVOSkCGisUT0IBG9TUQriOgMIhpPRPOJaJX8f5zn/GuIaDURvUNEcz3ppxDRUnnsFpJxsUTUTES/kemLiGiW55p5soxVRDSvlPsYKYzEDlnd0ki8N4Y5EChVk/kRgMeFEEcBOAHACgBXA1gghJgNYIH8G0R0DIBLARwL4HwAtxFRUuZzO4ArAMyW/86X6ZcD2C2EOALATQC+L/MaD+BaAKcBOBXAtV5hNlKxhjDXsCOOW3ZU7YtlDMM0JrGFDBGNBvABAHcCgBBiSAixB8CFAO6Rp90D4CL5+0IA9wshBoUQ6wCsBnAqEU0FMFoIsVA4nut7A9eovB4EcK7UcuYCmC+E6BZC7AYwH3nBNGKpy+iyMjlLbJM6c6zKMExDUoomcxiALgA/J6IlRPQzImoHMEUIsQ0A5P+T5fnTAGzyXL9Zpk2Tv4PpvmuEEBkAewFMCMnrgKaWy+1XrGSZMcsYhmlMShEyKQAnA7hdCHESgD5I05gB3VBVhKTHvcZfKNEVRLSYiBZ3dXWFVK8BMPS0SgvI1WIyZoXzd81pLGQYpiEpRchsBrBZCLFI/v0gHKGzQ5rAIP/v9Jw/w3P9dABbZfp0TbrvGiJKARgDoDskrwKEEHcIIeYIIeZMmjQpxm02Eg24rExEP9NIjJxjmAOB2EJGCLEdwCYiOlImnQtgOYBHAKhor3kAHpa/HwFwqYwYOxSOg/9laVLrIaLTpb/lssA1Kq9PAHhK+m2eAHAeEY2TDv/zZNqIxuqTaeBlZWzZsLmMYRqTUhfq+BKAXxJRE4C1AD4PR3A9QESXA9gI4BIAEEIsI6IH4AiiDICrhBBqQaorAdwNoBXAY/If4AQV3EdEq+FoMJfKvLqJ6HoAr8jzrhNCdJd4Lw3PSOyHReB/hmEai5KEjBDidQBzNIfONZx/A4AbNOmLARynSR+AFFKaY3cBuKuY+jY65mVlwo83MiqYoZZBDQzDxIdn/DcQto62tn6LeGW7ER6mtctKyr3x+c8n3sGsqx+tdTUYJjYsZEYQNfHJlBhfpgSnUcgc4CHMtz69utZVYJiSYCHTQJj6WdVB13LCYqWLHonRZUs378Wm7v5aV4NhKgrv0MGURKnRZZFFx8iTMfjrW58HULv9iRimGrAm00CMxLXLFDaz2wiUMXXPzt5BLF5/wAdtMiXCQqaBMJrL3OONt2lZVCplCly4Zhf+4/G3K5J3OalFdN2Ft76AT/xkYdXLZUYWLGRGELXQZKpVZKXu7VP//RJu+8uaymReRmqxZNCWPfurXygz4mAh0wi4O19aQpgb0KYUtc6VvrVcLXrxIuBVqItj+94BzLr6UfzhDe1qU0wVYSHTCFj6l/wCmdXviMq1tpgtgKAUc9FQJodte8NH5QOZbOjxWsNCpjjWdPUCAH65aEONa8KwkIlBvc4+r0WtlHApdb6MtZwSbu6fHnwDZ3zvKQwMmwXJwHAufgFVoE5fubqlrcnZD3H/UH0PHg4EWMjEIFuvppVaRpfVcfzXY29tB6DvqJtSzicwWGVNptiBCmsyxdEqhUw/C5maw0ImBtkaffDWtctq0NGX2hRRLy+lnEzIoCCVcFovk61u2xU7TmEZUxxKs67Vt8rkYSETg1y1LSvK8W/pkkd0dFkJJSnNU6cNJKWQGcrGe6j7BoaxcM2uoq9rJE2mXs3DYdSzZn2gwUImBvU6OqpJrUpsi6gdWDmaXJdFqZrMF3+1BJ/675ewu2+o5LqEEddC+9LaXVi/sy/exSWWXUvcgWAD1n2kwcvKxKBWPpl6nPFfPU2mDHloNRlnnDUcU5N5Z/s+AMVHpxX7rOJqE5fe8RKA0pauccqu1rTb8sCaTP3AmkwMTHMqXlnfjd8u3lTl2tR2gczq+WRKv7dQTSbmwCHpho8XW5dizWXF5V9OalX0hbc+j5+/sC7WtXVqbDggYU0mBiZz2SVyCY5L5syoSLnW7ZcrUqqlTLVUf8lL/luOl5S7zEOjrCRdc1k8Tcado1SkFCi2E6ylT6ZWZb+xeS/e2LwXn3/foUVfy0KmfmBNJgb1NzvcXRKgZjUo1TxhD2oohyZjdvwPx/TJqOsr3RE38uKntUA9jwas+oiDhUwMaiVjqjHaL5ZylVkNf5MuD4q4ZI8JKWOK9tMVK5RqGeGlK7pz3wC6egarX5mIsHCpH9hcFoOazZMxfDquT6YG0q/kpqji2mVhecRtukQipk+maHNZceeXE917d+q/LwBQv3vhuJpMI6phIwzWZGJQf+YyB9ZkbHkUZlLqRNZEzHXjig9hrqVPpmZFx4ZlS/3AQiYG9RrCXItqlWukaA9qKL0cXftQzOgwhYouK/ad4MmYlUWwT6ZuYCETg3qbjOmOxuusXsVQjW0MdIJKtV3cTpxi+mSKLa2Wj7YRNZlGrPNIhYVMDGzmsmqb02zL5FeS0ufJVG/Gf1hRcQV03G0Wip+MWdz5ZaUBO+x8aD1Ta1jIxMCmycSd2GfD1BGqOSqNvFJvNcxl2hzcoIkS8y5aNSnu9JqayxpQyqhPsPFqPvJgIRMDm2mkVh1CbZaVKXF+jLq8Co5/3XMp2VwW8/riZ/yz478YGlEwjlRYyMTANuqtVH9gXOrfXVamMuWGUa57rUanEFbXUksf0SHMNVmuqDyDlwZW7kcMLGRiYDOX1UyTqcV+MqVeH3WeTFkc/2H5l9qpVdrxf2BpMqWWycKlfmAhE4NamctMuVINV5UpmyZjM5eVwyejM5eVGMLs5l2GuoRRb5MxK15miS9WflkZlja1hoVMDGxCpNodgnL818SsUaaP2JZLOdpUu6yMm39pIcxFL5BZZDm1nSdT/TKrLfSZylGykCGiJBEtIaI/yr/HE9F8Ilol/x/nOfcaIlpNRO8Q0VxP+ilEtFQeu4Xk8JKImonoNzJ9ERHN8lwzT5axiojmlXofxWDrUCrV2dflfjJlCmG231s5NJnCtHL5syrtkznQFsgsdfDSyJGWI41yaDJfAbDC8/fVABYIIWYDWCD/BhEdA+BSAMcCOB/AbUSUlNfcDuAKALPlv/Nl+uUAdgshjgBwE4Dvy7zGA7gWwGkATgVwrVeYVRq7T6a85VmXPqmh479alGeajDmX6vtkGim6rBbmslIzKEs1mDJQkpAhoukAPgrgZ57kCwHcI3/fA+AiT/r9QohBIcQ6AKsBnEpEUwGMFkIsFM6Xem/gGpXXgwDOlVrOXADzhRDdQojdAOYjL5gqji26rNwfZdTcaml/rvSkzIqtwlymOUZFX91ImkwtyixZs4ymITOVp1RN5mYA/wTA2+1OEUJsAwD5/2SZPg2Ad9vIzTJtmvwdTPddI4TIANgLYEJIXlWhVtFltmxvfnIVBovcBrhUyrZ2mTWbMpjLSipfT9xdSRvLJyNC/64EJQt9Fi51Q2whQ0QfA9AphHg16iWaNBGSHvcaf6FEVxDRYiJa3NXVFamiNuw+mbIUE4snl3dWtTx1r3GXtqlmCLN2MmaD+GTqyfFfjaqUWkStNZkXV+/En5dtr03hdUYpmsz7AFxAROsB3A/gHCL6BYAd0gQG+b/q9TYD8O5LPB3AVpk+XZPuu4aIUgDGAOgOyasAIcQdQog5Qog5kyZNinenAeothLmW9udyFW3Lpyw+mZBMSh85V9bHUtvJmIG/q1Jmlc2XZebTP1uEK+6LOv4e2cQWMkKIa4QQ04UQs+A49J8SQnwGwCMAVLTXPAAPy9+PALhURowdCsfB/7I0qfUQ0enS33JZ4BqV1ydkGQLAEwDOI6Jx0uF/nkyrCtV2/LtfjKHcRvbFRM2oPOVU0vFfrpqU8YoyPZygQKyGVlX6ZMxaixlGUYmdMW8E8AARXQ5gI4BLAEAIsYyIHgCwHEAGwFVCCOVAuBLA3QBaATwm/wHAnQDuI6LVcDSYS2Ve3UR0PYBX5HnXCSG6K3AvWuptFWbv9xTHbPXWlr3ICYF3Tx9bfNmlrl0W+N94XoVCmBVxH1ncwIFqTMaMUsT2vQOY2NGEVNI83gxmUxXTXZXNl0zlKIuQEUL8BcBf5O9dAM41nHcDgBs06YsBHKdJH4AUUppjdwG4K26dS8GmyZT7G3TnklSovI/91/MA4m2lW3JUWVSfTGnFAAjveErtOCvuk4nRa9qu2Lt/GKd/bwE+e/pMXH9RweeXL7vA8V90VYqGHf8jB57xH4N6W4W5HpbOqLSwKc/aZWHmshLzrvAzjzMyt72H+/YPAwCeejs8WKTWjv84bcuTMesHFjIxGMyET5Qp+zwZ4f+/nohq7rLmY9MOy7J2WfzybVRKk8mvSxfHJxPtuM3EWhDCXIVBjfcbimUqLGNdmNJgIROD/UPhc1GqbQ/2+WSqW3TZJJ/dJ1OGMirhk4kpBIqf8V/U6ZHKUMcTFilT6JMpvi7F4m3OOIM29Tw4AKD2sJApgoT8FvuGMqHnlfvFzmsLpuiy2mEJfCsiB8tZZbjJsM6q3n0ycbQHWxmqznZNJnhd5d844dNkyn/vTPVgIVMEzSlnqbWB4XBzmS0wIC6mbOth8cSy7ZBpOl5hUVr6qr/FajLRyK8SXVT2ThlWc5nwlWGiFo5/v0+m+OvzS/0ztYaFTBGojsQewlzmcq1fWf543Jn3cXE715I7acvxCpvLqu+TKdZcVn7ntzpqNZcVOP6r7ZNhTaaRYSFTBOrFzVQ5uqxczvVKIMojYyI4/ksnTNuIvZ+MyrtCa5epTdViOf5tx13Pf/h5NdFkhP53VDi6rH5gIVME6rXNWlSVar/fNTWXuXWIV4nI82QqPBmz1OyLn4xZZP4xtGNbm2VlnsUqv9WZ8V+iJlPOyjAlwUKmGOSbm9V88F4TWrVDmOvhg6p0HcqRf7jjv7S8i3/k1TCXhR/PSMlFFnNZ4bIyRVelaPzRZXGuV9FlZaoQExsWMjHQaTLZEkdeUTBGl/nKq65TplxzeIpwO8UvI+RY/UeXxcAWXSZf40SR0WXe97AaSyhVYo4QUz1YyBSB+rh00WNZnyYTv4wd+wY05drqVTtsS97Yr/fnYyvHxModPbjyF69iWKdmqjwq6vivbHRZPJ9M+DWuJmMZmIT5ZCo1oCp1Mqa6ph5WwzjQYSFTBOq91y0r402L22G9uHonTvv3BXhs6TZDBcLrVRNK1mCiZWA77R9/+wYee2s7lm/dF5J/Jc1llfHJlLLfje0a1ZFb58kE/y7RlBWFUgUZO/7rBxYyRaBe20xWI2RKHHkBwPJtTgf5yvrdgYLr94Mp1fHv5mO53HZcPZNkwPYTtUOM3Sm50V/FXVb8jP84JiOLJpNVQqa4EOYoTvmunkG854YnMX/5jgg1LYQd/yMHFjJFoD5anbmsHI5/1UGarjflWlNFpkyT3ko1CSpNMpX0d5jettQ1q0qr9s6YUaPF8lsJFFkhRG8zmxcvKKyiCIANu/rQ1TOI2/+y2lZNfZm+8mNcH+Gip97egZOvn4+B4epuWX6gwUKmCNRrazOXmT68t7bsxZ7+IWP+KSlkMoEeyNoBe8orZTJmSQ7WCjv+bYJbtVkqERQy3jIK88hPJq22T6ZY81r5TUZqsFSKucy0InlCPgeN0h+Jaiwrc8OjK9DdN4RN3f1F5b1o7S5cfNsLGLIslMs4sJApgrDJmFnLiBlw9m254NYXjPknE87jMH24lZ5pXUr2cR2sUa+y1U21WXD2urdeuixUU5e+rEyR51uEX5BY2rHNxJiLJmSCEWTF+GTiRp+V6vdxl5UJudad6Fpk3lc/tBSvbdyDzbuLE04HKixkYqD7cLzKR1iHsDFk1ORqMln9R22cJ+NJLyWAOazeqzt7C5zqgNcnU0LBiCKkoo3Kw53Umucm06o9479YKiBj3Pe42FWYfYLbZNqV6WGrYwwMZ3H/yxu1eXgvi7VhW4RL8s+uuLxLWU/uQKQS2y+PSLwfgk2TifvyJV1zWUDIlBjeG5WsEMYX4kM/fAZA4e6Z1ZonY9VksmrkavYfaJupXD6ZIjPwazIh2kRJ0WU2E2M0n0zYZEyT1q0iycNWx7hp/kr89Nm1GNvWhPOPO8h3zC/ILBXUEGWBTHebhmK/H7fBWMpEgTWZiPhGxFpNpjQbMmDfm6TSr3S80bLw/V8pbLmrDjP4aHIWGZM3q+hLGBjO4jM/W4S3txdqcaZyomAz4ynyo+byj+bdd7ZIp0yUOSzqnJU7eo3ZdvUOAgB6BoY1ddOXVwl00aJhxNWADlRYyETE+z5pNRnLPJlSZkaXOsqPSikdWdw6RL0ubMvrXE6gs8fpsApH3eEj4rxPRp//axt34/nVO3Htw8tC61eKTyba+TGejeV4dE0mWBfvMX0pUd53FTmnVzBLG7RF8snI8sMm8OpIFOHLqcaKCPUOC5mIeD9ybXSZd3SneWej7DGTV99NdSguvVjCOvJKY+tEwzqau19cnz8v0PbCMiK2dUZhHSFQys6Ynt8h15Y0GdNykRvCbFVkAubbCJFf5VwLLk5WxTyOYY0m8+zKLry6oVt7fv6Z2AvhSaHsk4mM91XRdcY2c1mUDtzt0AwjR/POmPl028S6MEpaiDB2qSqf8ONh7bdhV5/7O1STCSnXlL3bnDbTU7FCpsjzK9FZmSLyggSL9jnlLeayMMLattjou8Ly7eeosHdd/pfd9TKAQh8k4NFkIpTBigxrMrHQdXgZn5ApvKaoj65IymYuq6BJz3id5bhqkzAhk/DMjSnsEC1mTEt0Wb4ftGlaoYcLEIbfxvMNJ63Ytg9LNu7WHosquIveGTPCApmRNHdNftoy4gx+IrSqqmKwrlG/gygDR9ZkWJOJjPddsS+QGU+Tccsy/G00l3l+VyqE2US4nhDheosmlCRCRojQ9vNOwCzsEPW/g+ebRssUcdRavI8l2nl57VZ/wUd+9BwA/Yg76mTMYkOYozjlo2gfeVNj4TG/thTjvXQvMV/rhr0HTtnZNxiat2ovNpdFgzWZiHhHRjrHv/dl0pvTKlOvIKW80qXsIx/lW/rC3a/gR0+uCs0niKvJhBSgJrEC4eYyXeOUazJm8Z2JV8Myn6Xu37Yba3gJetw8rUv9mzUZk/Av9X33+X0qsGEbkI8qCz67weHwAot5JtU2l23c1Y+Lb3sBe/cXRuzVChYyEfFpMrr9ZGyaTBGdUMFHbR1GR846lHiaTLgm4uWptztx05MrQ/MJojSJMBNG0vMWF2iBlhGxalvTvScswRjBfKIS9fSkYYJutDIsJr7Ia5eZ/zY9Fv8KGMXXvVRNJkrnnjeV+tOHLNFmUUy4wTKqxY8WrMJrG/fgiWXbq1puGCxkYqDdGdMWfRbhhYzv18hfWHVzmavJVOZjSkT4oJPk9cmYNZnwEGZ93lEjiYydbU7glgWrsC8wF8RvxrObAm0dn61OuucTeVkZEfzbLkD8EWj6fMMj96JpeiaihDCrdyr4bG0hzco6pBFrAAAgAElEQVRcFmX5J1Hl5c1U2TYTaDVhIRMRuybj/R1uTjNhmqVs8xtUM4S50HTi/z8upntwP+iQAvzmMv8x22TM/P2E34E5xDn8+IIVO/DD+Svx3T8uj5RfEHVvcTQZW0ed12Rs0WVBzTr/26ShR1kVoDo+GTOqXsH7G86EX6xay2gqjNA+lUKVZtvttJqwkImIzQ5dDsd/ztDfudqCsW55SluF2X5OwZI3wv9/3DJNlyfdUaN5SOg1lwXNat6/w9bIMmUf1aZujE6T9d/V61992zfaDSlDPU/b6Fq7wrSlo1bPMmHpBcI0GVO7ed93m5DRYau7jSih9SrfYNMOZcOX/lfP1OSTibIie6WIuhFdNWEhExG/JhOuqej6g2LstwWT36x1K8+LHOWDKBxR600OUbEttZ+3f5vzCNNkbP4DWwizO9o1lO2OCwzXt6Sdug1k/B1X1NZS2dqEjP6dDC8vG7L9ctiOot6/zJMxPd+D5d3QmQtL3X5ZXRL2XprMZUM2TcZ9J/XPpNQN10pBFcfmsgbE+6rYQpi1m5pFMZdZbLylvK+5nMCvFm0M3QMjiiAcNnxYldJkKEK4qM/xHxZdppuPoTQZQ/ZRfTFGTczguI/e+Tjn6Wale9GNqm1Lsyi5peuP/P6cwLEII3Xf92Csuzk8vNTtzN3BQ8g7bTKX2b4D1YGbTJi+gU2VfTL1GDIdW8gQ0QwiepqIVhDRMiL6ikwfT0TziWiV/H+c55priGg1Eb1DRHM96acQ0VJ57BaSPQsRNRPRb2T6IiKa5blmnixjFRHNi3sfUfEtK2Pbfjmm4990SjmCyx5+Ywv++fdLcevT5p0KI82SNmxDEFuTkZcZN7+K4vgP02Q0ZXmx+busnYRFEzJfF+00dT82TUZ33FslfUcebcAR5ucyPZaMZdBlw3s/lQqtz7rPzp8e3DQwiM0n41+RnTWZUjSZDIBvCCGOBnA6gKuI6BgAVwNYIISYDWCB/Bvy2KUAjgVwPoDbiCgp87odwBUAZst/58v0ywHsFkIcAeAmAN+XeY0HcC2A0wCcCuBarzCrBFZNxtP5luz4N3SUxmVlIrzHe/ud6KawnTmjjBiDnVnJQgb6D12RsNi/gfDJmDazi9VcFlGTMWpCsrmCh23CL39MaTLhHZ9uVG2LrMuG3Lt/OR7zSN88Tyafbuq0w9Z9896P6dl8+IfP4O/vXawv3/JeCiEwIOfDBM+xDQht82SimMv++OZWvLFpT2g5cVDPakQIGSHENiHEa/J3D4AVAKYBuBDAPfK0ewBcJH9fCOB+IcSgEGIdgNUATiWiqQBGCyEWCudtuzdwjcrrQQDnSi1nLoD5QohuIcRuAPORF0wVQViix2yjlygRqCazmM1c5k02vVvqnLBXL8qIs0DIWISEDdfcZOoM1Hlhy8p4bqpgxr/ludjMXba5JPn71+eQMUiZqDI5r8mEX6AzY/o0Gc0dutsjaN7NsGWS/EELEUbzlndf11l7BZOpbVd19mL+8h3aY7b3sncw48nfXh8vyocVZQVq03P+4q+W4MIfm3fJjUsuxARaK8rik5FmrJMALAIwRQixDXAEEYDJ8rRpADZ5Ltss06bJ38F03zVCiAyAvQAmhOSlq9sVRLSYiBZ3dXXFu0HYo8vKsUBmvsMLdJSB/wsrZ89bnRK2gGYU+7F5186Ymowb4WMecQLRTS4F/gPh/V2YR34ypj4/23Nz28xwmmqvwmdq1hR09auEJqPeWZuPscBn4QtyMQiZKJpMSB6ZCB11GDYNu28wH4gRHMDY/EH5iD+TJqPPqxqod6mOZEzpQoaIOgD8DsBXhRBhOzvp7luEpMe9xp8oxB1CiDlCiDmTJk0KqV446l1JkH1nTNtkTXMZ4eq3MTDA89s05yHKqx5lHapgh+FqGjG/pXyEj/54FE0mTJB4/9Z3xPrr3PpZI6PCrzeNiqN2nOq0OELGFkqr6mbb/yhccOvr4yvbUPUws5P3fuIJGSXc9fg1Gf9ZPn+Spm5pGWliCqKph+iyEaPJEFEajoD5pRDiIZm8Q5rAIP/vlOmbAczwXD4dwFaZPl2T7ruGiFIAxgDoDsmrYqiXJZVMWB37cTUZ82g+eh4m1EcX9vJF2RskGN4Z1SdjNqtYylaj7VAh42374PWesnQRWBbHv01Di6sJeVPDfTLO//bossIOzx7CbK67T5MIXG2be+TN21Q3ID8g0rWR3/Ff/HtvM8OGmcu8ARG6QUZzyuk29w/r59P4rRqRqls0Qgh87TevY/F6/543qrwYC0RUjFKiywjAnQBWCCF+6Dn0CIB58vc8AA970i+VEWOHwnHwvyxNaj1EdLrM87LANSqvTwB4SvptngBwHhGNkw7/82RaxVAvejpBsRbIjGLuMUXEuELG2BFas3YJm91tGnV7OwlTh2H7mGwOYqMWZalbsOxQTUb73FQ99HnbPlabkM3vWRK8LtpDU+fZlpXRhjB7zWUhc7dsjv9g20SZA1PMhESdFpYp4nod+eeiP947EKLJWIJ4mqSQGTAJGa9wr5Am0zeUxe+XbMHnfv5K4IgclFVZgwqjlKX+3wfgswCWEtHrMu2fAdwI4AEiuhzARgCXAIAQYhkRPQBgOZzItKuEEOopXQngbgCtAB6T/wBHiN1HRKvhaDCXyry6ieh6AKqFrxNC6LexKxfymaVTCQx6XlCFbVmZ4iZjGo4bR8XRBVioJhPBvh4cUbsOVsv9mZcfsQiZCJpSmCPa+2cwZDfSfA+LBmhbI8s0mVMYfhfmX1hXHTpzml+TMb+TNs08eDQsvNk9xyLcvXXShVJ7O/A4/aXtvfJqMmHzZHTfrdq/KLgeXbBsoHKdven+or4v1SS2kBFCPA+zf+lcwzU3ALhBk74YwHGa9AFIIaU5dheAu6LWt1TUM0slEvrJlmV0/Bemq9GJ/niU91jlEbamkVmT8XQYwc4sghAAwjQF9bHoj9sCA4JlF0RCebrI4P1lQ65z0y3PTR01jVhdx3+BKhOarec0+/17y/FiEwZhbR/m/N7Vlw+DNw5MLNqAky7/17TdwHBp5jJvRKYQoiDgxauRB19pm09G5b2uq6/gWPCaSk3GHM6Ea8i13Eo9CM/4j4hrLksShNBEpFjMZVHs2MZ5MpbRSZRvUJ0SFl22ckePNt3bYRRqMrJuljrYNBlzm/jPCztHd15YpI//b7sWpy87XEjaRvFAuEklyv075ehCmMPLCDOX+YWM/9g1Dy311C9cAwzm5UW9z7o28vo74vSXvidrEaJh82TC9oYy+WREyPtoyrNYVJsV7Oop/6wncxkLmYioR6YiS4IP0basTJRRs2lFYNtcDG+qzewUFnTynYeXadO9nUBwPobNca6wrVhrPm4fmYXPTg/RZCKYfayjaIsQMM2qj9wHRBQyusAAv/DV1C2kbX37wYSoXVECHkzPzu0oNXUf9HTgcTplW4SXP3rNf9y2WoFK6x8y+WTMAtrJv3T1ZsjVZIJCRvYVrMk0HuqhpZL6iBi/ihw+ejG9ZEptNzr+I7w46aT+kbodRYzQRl/djZpMTCFjM5epckPuPWzEHtbZRVn+w752mfDVM4gbJhxI92YbVkJJ5jKLoMiHMBfmF0UAh9Ur0qoAwqzJDGZKNZd5y9HUL2TQ54suCzGX7TcIGZ8WZxFwcQl7dqZyawULmYioZ5aW62SFCRnbPBnzsvLho/2weSwT2pv0mZYBX3SZcVmZ8Dxs6zzZ5ujEnScznPHWPTDqi+LAVjOoDWWry2yhvDFdMh7zR/h5uhn/1mVlDNsPA+HmMi9RZvzbBJHuuDeaThsVaDVjeuoYEvTgnBuiyYQMGM3mMosWVQ5zmWyfoPWbNZkGRr2oSpMJvijq4SYTZPhoPb9tHWrBcftoVvlaTGfkzWXFqzK+6DLLqNSEbbVem78prJMNE+CDWe+o1H8wykq//UOFkYS6sk0DB/NkzGiduK193HI0DWQrI0zARw1BNoYwR3H8h2gy3omOOpOjraO2Le0StmyOfR1C53/Tu+ETcFpNpgzmMiVkAt9zMVaPasFCJiLqmSlzlGkpinSSrCHMpqXPTaNWW0crhGcvektnH2cmsM8nE5jlnB/Jh+dhm2hq84nE1WR8nVXQlBbBXPbdR1cYywXy9Tf5LUwdSpQ+IEqItcK2W2uYINGak0IE8AfelV85w6gBWkxGQP470NV9KJMzbpNgqx9gn2eT9U72DNyETZNR+e3YN2j/1i3rwsXFbRODJlMGi1zZYCETgU3d/ViwwlmIL23QZNSHlE4aQpwjfHSmDlX9FTZPxiY83Bn/4afp6xXiTzJrX35sJhObqTBsZBY2Yh8OMbt4byWuL9YmJE0+GV87RumoLR2T3vHvaRdd/iFtbyr7/pc34tmV+TUAo834D9dkdB3xUCaHtnTSWD+veVAvJC3HQwYm1ugyz/m6iExr0EEZhMywq8n4yWsy9TPlv5TJmAcM5930rGt/VaOrglBZ+eI0RVh2plj/RJQIK6U2m85Ql4YJo2Omjtam+zSZmI5/W6hr2OKVwTqE5R2qyQTq7u3obRNak4YJRjYhm/c7mE11JoqaNV/mEGbvO+z1j9z94npjHU3ppsFRztA2gNOJtjQl0TOY0a8V6AurzyGZSPqOWzWZECFkjS7zHN+7v3BCZs6iQZbDXGZaZkiljohlZQ4kvA4+N4S5wFyWP24b/Rg/TNd8EOjI3dGy4cUSdjOYujS4z4S381H+poJ6eUelpsmYlpfabC5T2pvuWP53qF/A91H7j0XWZCx9vml/DnWdqXqqzKDJxxdCaxBwUd4ZhS2EWeuTyZnb3lueV1C3Nfk7c5uZE7BrMlqfTDbnlqW7d98EYcugTrekTpgQskWX5QTQLuvWq1n9wx7CHC78o2BaMDWqZaGasJApErVBlikcNpmgWNszA579PQLHbS+OgKcTNLxbZhOdp3zD6ChckxG+OpqwOv51mozmPB2+XUvDNBnDMwvmoUszajKW+1dlBjtCnxZlemYhzukg2/cOFKTZlsuPOhlzyPPM25r8xg+TluLbdMyi7eie7WAmh9Z0mJAxRzw6xy2aTNR5MpqyhRBob3baQRdhZhsc2NZGi4K6/+DYhx3/I4CUcvwHXszX5S536STZl/go8qMTgeNBooyG1OisUEvyjghN83fy6YWTMZ3/be/0sKHTCTXZRHDMF5wXqIiab9GcShQ9GdNrJjIVbwrWUKgRZ7AjjGKXt5ldgLzw6+zRCBnLSsZKyOq0IK8A9grqlnRQkwnXUoAQTcYggFWZrWGaTDY8f58mY6lfgbnMIgSyubyQ0S2SaffJeL6nmB56tRp6MLosE9KmtYKFTJGkXJ+MP105Q1PJhGHZdfvoRV0XPGyLLsuJfGdjjHJSAiyoJXnPMWky3g/asGlZMSvt6j5wXZNENRd5DwXNCOojbm1KFtj+bZ3BYEhkmiK/VpRJQOs7clsnGCzTNqFRZ7bJWMpQxwczhR1lqeayKP6kME1mOJvXZGyTNeOs4pzN5YNlguXbo8vy7aAXMvnfWnOZp77f+O3rhSdEwKTJqGdVT45/FjJFkjSYyxTpZKJgzxXn/Pxv40KUhslxtgmJOSFc4Wc0vWT1eUSynYf5ZDT56BjK6s1WYavlepNsjn+TGVN9dG3pZKHJyrePfGG+g95FGi1tYxqR5keWAU3Gc/7vl2wJzdtUvhDCbaPewcLOzn9/unfSqdOgZvMt39woz7Nrb/YLmbCl/lMhIcjeOhk1mRBzmVfw6fwTfse+XgilEwkQacxlhnfVm197U4i5zCLgvO/Cn5ZuLzgeBdWmQSOuGjCUY1WBcsFCpkjSBnNZ/jhpNZkoczLimsuyOeEuP27C1WSCfntvR24QIN5OYMjgk7Fp56YPN2w063P8h/pk8s8l2GE5kUeEdKowIMP3HLSaTFZ/roe8Y9+gyWTV8eBoOX/+KsPCpLawd+/96LWRcHOfqtNQJhe6BbFfk/H7ZMIWe1X7rpgE0bBh4AM4g5Iwc5ltRYBhz0DP5I9KJBxfZoG5zKrJCLRJYbt/KPxbt5n64uJOxgyoMio9rhmuErCQKRLTiNl7XPcShTmgFRmTkJF/hkWXJcmiyRh8H1G0hTBNJrK5zJuHZqSnu9xr+tujCRX15uHOXwq0/VA2h3SSkNRsNmf1yfjMZYayLTZwlR40q3jPN+6wKBslQfr6eTsS3VbA3uM6YeC9JqjNeJ+nbsO0/3XSNAAh73LWI2Qs5mHdoGwok3P9P3FWBBjwCF3T9s6pREK2bXAAYBPuzqCmKZnQPjt/VGTBYV/+7ztiQuEJETAJKqV93/XCOqzp6o2Vd7lhIVMko1vTAEKETDKhVd+HLSq4N73YeTLRNBmD4z9kvxVdevAc9Zdt2RNvR6VbM0zv+M//3ry735h3Tgg3ICPY4QxlcmhKJpBKkHF7hpRhKSDfIo2mthd5bUCHave+oaxP28haBIRzL845acPcK2+b6gSBze/j7SDDhGDQ1DmqJYV/veBYAOEhzE2GcH+3DEPIPuB8L02phBOtaRBC+XMLr/cuXjlocM4nEwSiwmjQjEVLEkIgSYTmdMLq+A9bjaA5lTA+exumEGbvO/v8qp2x8i43LGSKZFybI2TUuzOYyWLW1Y+6x9NJ0r4AUTQZ9TEHD9tnxQuoxZdNXb36oAtDac1aisLbOQbvTVXJFuDmHXkNa+YhhEWXjW1LY2A4F9IZO2YPx1RZqMk4nVVhdJl3ZW2tT0aWpxvtBusfRUDv6R/WptvybjKsIuF9FlpNxuKX6B/KuoIgqMm4ZQc6wkzOEdq2ZYwyObsm40beGRz7TcmEVgMFgKGsR1PRCZnhrDuXRedzyuRySCXI2U69wJTpeTaaspWprTWd1AoZm4asvrO2pmSB+Tkqphn/3m/VtD10tWEhUySHTeoAkB/B9nscricdMhbpZGFnBoSvoaXY0+/sOBjc4Cw/f0Zfp5wcWYWRd1AHOpNsviMzdZTefTMKBWi4AFT4RocaW3zYZMyxUnvsMWx3m83lTWLBDm373gEMDjsdiimKKJ1MaDtLNQJuTSet2qfNJwMA3Z4dJcN2ZlSoEXqzoXzvs9B1pJ37Bt3fusezfziLce1O2wY7JPU8W9NJv5DJOhqAmpcV1i52c5lZCxzK5NCc0mugwWt0K1APDGcxtq1Je2+qTgnpqysIL8/mXCFq8skkiNDalLSay/SDAyetrSlVgiYj8/V89kKIwIAgngArNyxkiuBTpx5SMAvZ+3F/9PipSCUSWvX91qdXu79NE7ze2LzXOe55Mb15mcNY8/uOh40snfz08zWaUwmjnbdf7ofe3pQMCWHWXuoypBl120KI1b2MkUJmnyZMF8h3fLq2f+rtTvQMZrQj4pxHyISZy9qaU8ZnppJt0WUAsLt/SJtuEtDqWbWk9eYyr3Nb11l5l4DRrc81lMlhnOqIA4ED6jm3ppMFqyakEuRZXklb9aLMZcGOWgghfWlmTWbQYhnYP5TF2DYlQDWaTNa5j3QyUaBNZHICzSlz0EFOOJpzazqp3VPGZ6YMiS5ra0oazV42MhpNJpMTvufx/554J1be5YaFTBEcddCogrXLvB9IczopJ2P6X5xN3X5/gu7FXbh2lyfG3WPCsixxodLdEGZD3VU+wc5I5dmcThonY/bLexzTmi5Y6j8sBNlXfmCtKcAecaeSlB9sn8H5P5wTbodkmh+Q0hzzrtIQ5vgf25rWjlhNz8l/jjMiB4Ad+/ITJr0aTmtggqNCtVNzKqF3/Msyg9qGjmDzDnieqfO33vHf2pQsMPWmkgl3foYx6i4r3Ps2R5c5+QY76mzOCc3O+2QsmkyBCVdg/3BeyJjmASUThCaNDzWTFWhOm+uu/DnN6SQGDKY477lBVFuPbk3H0mSEEOiRAz9v9jptNsyXWS1YyERALW3+2dNnumYp9eJ7VfHVO3rkZEz/i+U1kwB6G/Karj7tcd9M+RCfjGltLYV6sQs1Gc9oWeht0P2DWSQIaG9OFXQIth36guUAHt+Tpyo6RUB1ouPlhmw9Bk0mm5NhyhqfDADMmTnOGREbJkSmE6Qdcb60dhcAYFx7E/oGC8s2aZxeMjnhmljf2d7rS1cj/YPGtGivzQuZZKhPpqMlpe1Iw1Am0HEGk5J6TZpTiYJw4VSCCr6DIE7EXzRzWVCAq/KapLnMGsKsiSjMCXjMZYWd76D01el8qNmcR0AafDJEQGs6gQGNJuPXsgoOu209uiWeueynz67FHc+uBRAI0JB5qfcKALbs3l90/uWGhUwEvnfx8Xjmm2chkSA3ikm9mN6P88qzjkA6QQUvTp/c3Ogfz3sXAOAr9y/xHRdC4Nv/8xYAxyTn7TCUVpQgaDs6wBEM7ox/Q2ev6hk0DbiaTEo/z0TVv70phY6WVIFWZppAGsRbrqux2aJw5DWqI9xn8MkMZ/MmnGCnMKY1jWMPHo2UZp8fVedUMqFtt589vw4AMKG9CX0Ws4jJ7JHJCnQ0JzG2Le3b5CqTy6GjJYXx7U1GAaEm9drMZR3N+s5q0qhmTB7V7LtXhXofxkkBXuj4l1pSU9An4wj0/F4vpqg6u09GDVj6A23r7Sx1gwOn/v5gBN8xOXdlXIgmM6yiDjWazHAuF2ouy+YcH2hrWu+TsW0d7QqZ1nQsc9m9HjOo997Vff7bhce6acEBbi1gIROBaWNbMXNCO4D8chJ90uHvfckOGtMiHf/+F0d9TOp92hV48L94aYP7u70p6V8yQ77kEzqasW8go/2oc8K8OrRC5TkUtL3nVEem5iQU5t8/mEVbcxKHTmxHb0DQ5ZfCsZnLCkNOrWtrybrlNRmT498ZNet8MgPDWbSkk/JYoVkE0K9r5mVcexOGMjnjkjXevHR1SyUSaEsnfZ2pMnE2pxK+lQX8+SstU6/JqNG8E6WkcX4PZd01toKXq/d2XJve8a/aPph3JpcP/W1OJbQmGsAZ+JgmyKrjqg5DmZx28mc6lUBLOlngLwL8A65g26t8x7ZKAapp32FXkylcocOryZiCeJrTCaPj39smusFTXpOJZy7b6lkMdTgr3DLUfTanEvjZZXMAoOB7rQUsZIqkQ3606iVfvnUfAOAnnzkFgBMOG3zp1Z4Tx0/X79fynCeevTnt/3DVi3PUQaMAAG9vL5wdnhX5j0I3MhJCuIIu2Al7O1pA/1Ht6hvE+PZmtDUVhmzaIt+C5Xjr6O289g9lcex3Hsc1D71ZcN44i7ls2DO69s9yFxiUUUptTcmCEbNq5/bmVOgeHy1yVKsWQVWod8C4ykNOYO3OPiQTMhLJU37vYBatTckCc5QXVd+2phSEKOywVPu0N6cwnBUF2s7+4azr7wkKcbVqs9FcJk9vTad8nXTWIzyaQoTMYMbpxHU+D3UccLREVVeFao/mZAIdzSntumzezjP4zrpCJkSTGco6mkyTxlw27PHJmCLbmpJJtBgc/0M+TabgMAaGneg10+CgWFQZr23cDcAxr86ZNQ4AC5mGRI0M39rqRIKp7XnVaLujOYWegYyvQ1Aq6ykzx+Ocoya7AgMAnnp7B/68fIf7d3PKCVdVnd4H/t/TAIDTDh0PAFi3M++7AYBFMmBAaSK6D/r//s9bWL7NEYbBkVPeJyM1Gc2IvLNnEJNHNWvNA1HWNQP8JpHg0hfJBGHv/gz6hrL49cubPHVzjo9pTYPI7PjPZJ0Z/6mAT8ZdgTmdRHtzqsDcmBcyyYKABi8bu502v+QnC33pqzsdH8u4tibtiPTh17dgZ+8glm3diymjW7B5T94+vrNnEBM7mp2O2qDJ/MMvXgUAjG513rlgFV2fjHwnh3zaYg6ZnHAd+0FhcNldLwMAJo5y3tt7F27wHVdCZ9KoZp8GOZzNuZGMzSlzJ9kzkMGolpRxwuKK7c77OEYKAq8pUb0rzekERrWkXCe3F2/nGXznVV5hPpmhjBO9prM8ZD3mMuMcnlQC7U0p1xTuZa8nilA3EXRACv8mOcCIu6eMQt3/1x94A4Aj/FU/ZRqYVRMWMkWiPuifv7De7WQA4ORDxgIApo9vw1A2h7UeYfC715wFEEe3pDChvck3Ke8Ldy92fxPlNYrBTM43Apsxvg1A4UvzyTteAgDX/q1zQP9q0Ub3d/CDVCOxUS3OfRXYt4ezWLezD5NHNTumi+Gs76Pwlhc2alrZmdfAVIc8LP+f1NGMnb2DBdeoujbJEa0phHkgk0VzKlngJFadd0s6ifamZIFfRdWjozlV0C7zPYL/B5ecCAA4dGK77xzVUR80pgX9Q9mCzlStzL27fxjj2pt8nfXO3kFM7GhCe3MKPYPmJXOAfARYsI4/e26dW3/AL0jUYGDKaMcnY+pspoxygg5e3bDbl943mEEqQZjY0YSewYzPB6m2RW5O6QUI4PjxOppTxgmL//J7xwe5Q2pUA541wDbucvx+Mye0o6M5rddkBjLuXJbgwGiDvP6wSc7z0mkyG7v7MXl0M9LJhC8UXOWnnOe6xWqHso523NFSOHABgB37Bl0tSvdN7Jcm3CbpC3QnJOcEfvHSBqNZGNCHqmdywvdNKjNgSzrBmkwj0pLON9nVv8ubdlRAwEkzHGGzUpq1Vnf2YoXUIogI49qbfPMlvCSknRtwOgzvx3XSDKX+5l9A7wunOhpTKK2iq8ffmSuBN1l2NsFR9Wsbd6NnIIOzj5qMlnQSORFcZiT/W7cVrUp/9M1t7t/DriYjhYx0TgdRnUc6SRjdktZ2lI++uQ1LNu5Be3MSyYBPRtnyW9IJV5Pxfoyq83HMZf7O5MbHVri/x7Slcfph4zGpQ1/POTPHI5MTWNvl1zLVKFgXKtvVO4hJo5oxsaMZu3rDnbPTxrYCKAz8eEYKsYmyXt73YcAVMs5zNWmBRxm23O4fyqKtKYnDJrVDCGDVDmdA1TOQQYcckMwY34o1nYXrYwkh0DsghUyT3mIIHNkAABq6SURBVKR01FRHm7/6r452yhvO31unfEenjG6Wmkxh3XsHM66mEvRDbpOC67CJ7UgmSKvJ7O4fxkGjW5DWmCsHMzl3telg3nntOCEHJ6JAiHX2DGD6uFY0pRIFWpgQAr9ctBH7BobdgaEqf8mm3fiX/3kL1z6yrKC+igUrdhSk/dWPnsNNT64qSB/Vkg4VWNWChUyReFc9VS+JF6VxbJahg1/81Wu+42Pb0hjM5LB/KFugJjvrIamlMLLuKOTGi4/HjPGtIPLvG+Id/avRbnCkEyxjV9+Qm282J1yTjBqle6NRuvuG8On/XgQAOPbg0a593/vRDmeFK3j37S8UAoOZLG5+cqUvLeiTMQkZdV4q6ZhNdNFlV8n2XdPVVxCOqjrlVmkuy+SE20ms7erFwjVOiLJOkwmubmsyjQDA4ZOdttuwyy9kVNMfM3W0bzO7rXv2Y0//MDJZgYkGLc4rsJXfpNedGyF82sERk50QaW9nqTSDyUrIBNru/bMn4oTpYzC+vQlnHDYB75E2fMWe/iGMbk3joNGOgNuzfwivbdyNt7f3uMu1HD11NN7YvLfgHRvMOKa6jhalyfjbdvveATwktfvp45z8vYJIDYQmtDdrfTJCCLy9bZ8rfIMz/jt7BtCUTGBMa1qrbe3sHXQ3RUsnCn0y/VILAwqFjGrHUS1p4xbMO/YNYsqoFozS1H2LNJkOZ/O+LfXNbpSRm2ERYX/0DNYUG7v7ccuCvJBRpvVxbWmOLmtUrvnIUQD0+76PaU1jVEsKm+QkKOWov+rswwHkO4wHX9uMnYERrM9cNpxzR+7j2ptAROho8tunvZP71ITFoJak0y7WS1Pen5Zuc0fbqqPa1Zfv8FQsPuCMptXH6A1j7hvM4OAxrcayHn1zG37+wnoAwOfeO8u9NyA/KlQmHSB//0D+gx7TmpaajHlUtrqzF2Na024d9vYP47o/Lgfg+MtUh6Ds/ef84BnXF9auETLBjmnfwDCWbd2HXb2DEELgs3cuco9dcMLBAID1u/zh3SrM99ZPn+SYZWQZb0t/xOGTOnDQ6Bbs7B1yOx8AeHldN074tz8DcEysSnNQQub3S7bgqG8/DgA4/9iD3FF3v+fdUOYyFcLsNdGqvEa1OO9MW1PSjZZUdPU6frjJ8tms6ezF56R5UL33SpP6Q6DjU+9tR3MKzRo/nleoqoGLT8j0DmBsWxpN0iTVG9BA13T1YeveAZwkTdReLbS7bwjLtuzDpFHNIHI04OB7+Q/3OQOrjbv6ZXRZwIQ8nDX6stRAynkuzjnBtuvsGcBkqYUFzVXebz6oyfzypY0y7zQWr+/GR295Dr9fsjlfr6EsXlyzE6cfNt59n4P87sozXP/qwWNb8daWfdrzqgkLmRicfdRkAMDzq52osDkz/aPAGePasHn3fggh0NGcwrlHTcY/nnckgHzI6Lf/5y18+dfOfBllO04myA2R7uwZdF/QUXJU1d6c8o1MOj2mLxVM8OOn1+DNzXvk8QGceN18AMCd8+bgl393GoC86cS7L7zSZHb25PP3dvipZAIbpHB5UqrsnT0D6OwZdE0yui2Ad3s6t2+d7whnFcL9a+krUqNlAL5NpHb3OdeObU0jmSC8tLbbl/ddch4LAPzgkhMwob0Jr2/ag7e378N3H12Ov7zjdILj25tcR6jOht7RnEJO+Du6oJBRIezPr96JHy1Y5UYEnnfMFIxqSeOg0S1ucAXgCLlH3tiKkw4Zi5kT2n2dmTKPzT32ILx7+hgAwPtufMq9VmmXAHDf5ae5o2o1KvaGvF/9kaNwiNSelZDbsW8Ac29+1r23qWNa3FEyACxe340lG/e4nVFbc6pAEHTuc8x5h0/qwKiWFFZ39mKWfEcuONERqkdOcd65tZ4l5YUQ+PzdjjAa396E1rR/OfxMNoeP/dfzAIBZE9rc910J/7tfWIdfvLTRXa8ub5LKd/Yf+uEzAIAzDpuAplTCFwzzpV+/hudX78REKVzHtzf5vpnXNu7GYul/ak4nMbrVryH3D2UwMJzDhI5mJBPkC0gA8gOp0XIw6U1Tv3f2DmHyqBaMakkXCPe35Tuy4BsfdIXM4HAOu3oH3Xo98sZWfOInC7Fs6z587TdvuNd+7TevY3f/MMa0prHkO+fhd1eegSBTx+S/pffMGo8te/bj1Q27sdKwZ1E1aGghQ0TnE9E7RLSaiK6uVrmHBRzA93zhVN/fk0Y146m3O7FhVz96BzM44/AJrvlFdVaAs5QMAMw7YxYAZ4R4xuETkUoQnli23e1UVAe5fd8AHn59Ky768QvoG8y4QubRL5+J9x0x0c337W3OC7V4fd6ZO2fWeFeLumfhegB5TeGLZx/h2vW9y1Aok9+ifz4XAHD5mYfKe3A6tVNvWADAGZGNak75QrGHMjl85EfP4XqpTbz1b3PdTahufnIlNnX347evOqO0Dx45CRPamzBldDMGhnPokiPdpVucCL6xbWl3ZPmTZ9YAcExTSlP5/Ptm4eOnTMcEeQ/n3/ycG70E+IWMEtyqYwbyo/2bnlzpCmg1OlUC/voLj0MqQfi3PyzHzR779/RxTj4nHTIWf3hjK/7PL1/Fup19OOE6RxNRi1R6VyNQg5NJo5px8Nh8p6CEvrdTnDG+zRUySkPwWnBmTWx3VxRQGuqTHrv95NHNmD6uFdv25jWlT8goOaW1tqWTBcJ3VWeva8acMa4NG7v78aZcW++sI51B1g/+5gSnPp4KLd+2zx09T2hvLljf6z8862n94G9OdIWMEkT/+gfnmSqBqTT0o779OJZu3uvrLE+YMRbvmTUOSzbm3/MXVjvflBpITeho8s1Lu/i2F93fX/3QbEyQPjE1sPk/v3TMr8dPH1Og4a3a0YOP3+5cP7ol7Wrw3rZ9ZZ0zEJo9pQPTx7W6Fg3Fm1v2YnRLCodNbMcM+e6s3dlXMHdOx+PLnF00U8kEmlIJ993zogZ8APAuOQj4+O0v4rybni0Iwa8WDStkiCgJ4McAPgLgGACfIqJjqlF2KpnAi1ef4/6tOjCF+njP+s+/AMjbxQHHRBJktjRVvX/2RIxpTePgsa248/l1+Lt7ncgzFR6tIl5e37QHx177BH7/mtNJqxGlEnbv7OjBS2t3uQJjVHMKY1rTmD3FKeeJZTvw6oZudO5zIpz+ce6RaG1KYnRLCrc8tdrtFF7ftBvvnz3RfXFVh3vLgtWuyQdwVkQ4c/ZEvLK+Gzt7BzEwnMWOfQNuwMPn3zfL7SgBxx7tne9z8NgWvPrtD+PGi98NAFi2ZR9eXtftLvA4uiWNL55zBADgxsfeBgB3MVEgH7Sg2gmA29GNaknhoNEtbvm7+4fw7Mou38heCb87nl2LC259AfuHsmhJJzB9XCsWSgHb2pTEMQePLrBxHy0d2EpQ/WnpdpwtnzsAXHmWYyZtlZMxH3ljK7bvHUBTypnMN3Vs/t04/XsLcMuCVa52M3tyB8bK9wEA7n/FCe9WHeg35x4p28e5t18s2gAhBF7bkO9Mjj14DCaNasb6nf14fdMen9npuguOc65vTaGzZxD7h7LYvLsfX3/A2Xd+hxSQU8e04GmpFXppa0rh8EntWPB2J3oHM3hmZRdWbOvxHE/isEkdeGd7D3bLdvOaYE+ZOc7dabN3MOMTRlPlUjvj2/LP9KElm3HeTc+6fx88thWHjG/H2q4+dPcN+Z6N2u9mQnuTa1ZWAwgA+NGlJ2LK6BZMaG9CJiewb38GW/bsd7Xf9x4+AUI4i4wqIfr9x/MCcsb4VncFa2WiHhjO4styNY+zjpyMmRPasam73/XrvLVlL361aCNam5IgIhwjgy7e2rLX1W5VX+DlT0u3+bYTUeZ6r0BRKBMtAJx95CT398wJbThBvlfVJmU/pW45FcBqIcRaACCi+wFcCGB5NQr3jkCD3PbpU9z5LQB8TtWmVAI//JsT3Jh2AHj3jLFY8I0Puo7Mvz5hKn789Br3uBpRvnHteTj6O4+76a9tdD4aFdn2QbnG2p3Pr8Odz6/DMVNHY0xrGq9/58MAnFUBLj55Gh56bQs+frszmj1RRsMBwJfPnY3vPrrCV8bcYw9yf6tR5bqdfTj/5ucAOKPBg8e2YvbkDjz21nbM+e6Tvrb4yHEHuaZCALjuwmPxnYeXuQ77J7/+QVdInHH4BCQI+Pzdr7jnv3/2RCQS5Jtb5P3gVHsB+dUYFOccNRm3/e3JSCUTrjlQBTIojjpoFM70aIEA3Pu/8MRpPuHo9YEBwJfOOQKfOGU6AP9AQvGho6fgM6fPBOCMugG4JtIPHzMFgCNArzzrcNz+F+d5/3C+EyTxuffOcjtK9fyfXLHDvfdPzpmBq852BK/Skjfs6scnf/oS1ssAhL+TmueY1jS27xvART9+wa3bN+ceieNlp3PyIeMArPM9dwD4l486kV+nHzYBC97uBAD8+/863nfOmUdMxD0LN+C4a5/wpR83bTTeNWUUPn7ydNz5/DqcdP183/EfXOJoQaqjvOahpbjmoaUAnLUC75UDJhWFBsD17QHAH790JgBHy+0ZzOBkT/5fOucI91tIJhLYtncAc2961g0QuPvz73GPK/OS0jwB4NRDx6M5lcS5R0/Gw69vxf/+xavY1Tvofm8TO5oweVSL6yP81u+WYumWvVi6ZR/6h7I4YnIHOppTmDWhDcNZgev/uBwnzhiLr/7GEd7nSHP7mLY0Zk/u8K2W/IO/OQF/Wrodf3vaIdi7fxgf+6/nXe0KAJ795tk+DebZb54NIuDDNz2DD8zOCxXA6Rce/+r78aVfLcFVZx9REMxSLRpWkwEwDcAmz9+bZVrVePv687H8urkF6YdMaMOXz50NwBnheu2kAHDxydPd7Wvv+twcdDSncPikDtdG/o0PH4mD5Ujun84/0k1vbUpi/Y0fxc2fPNHNS3UkioM9iy0u37YP754+xvdy/cfH343Pyo4PAN41JT9y+uwZM90Z2IpPn3aI+7slncSHjp7sO/7X0ul98cnTC9oBcEaUXk3vopOmYZRca2vmhDbMmpD/YFrSSXzkuKnu303JBO673PEjpZIJ3PW5OQX5r/veX7kf3V8dP9Xt9J06TXPbbtrYVgS/sce+8n48/tUPYGxbE/78tQ8U5H3Mwf7wXtWpK776oXe5bfvpUw9BkJ/Ny9dXmZgUUz3P6VvnH4XzpNBRKL+HQs27UCint+L0w5yIopfXd6OzZxCfOf0Q/F8pJIKdD+A3+Z45e2LB8Wv/+hjXtPv3HzgMX//wu/C5987yvQ8A8O2PHVPQrheccDD++KX3u9qf2hpb8b2Lj8fH5XNKJqjg3m6SZrhgPRXXX3gsjpvmCEjVYSumjmnB1z70Lvfvy85w3vV3dvRgbZez+sJZR052n9v7jpjgG8AAwK///nQA+ec9f/kOV8B89N1T8dQ/ngUAvgHIL17aiDc27cGZR0zEzz/3HgDAuUc7z/TuF9e7AuaKDxyGGy7KC+qPvjv/vh8xuQNHTx2Nqz9yFGaMb8Nx08a4mvC8M2bimW+ehUMm+E1kh0xow4zxbXjun87BLZ86qaCtjjpoNOZ//YO46KSqdo0+qNTZprWCiC4BMFcI8Xfy788COFUI8aXAeVcAuAIADjnkkFM2bNhQkFclyGRzePj1rTjxkLFaE1k2JzCczc/UL5a3tuxFa1OyIO/ewQzufmEd+oeyeGntLtz0yRN9fiDFy+u6sX5nHy448WBfHdQeKQRnuRoVZullYDiLDbv6MXtyR8G2z/ctXI/nV+/Eu6ePxeVnHqq9v77BDO5ZuB4XnHBwgV05k83huVU78T+vb8FlZ8zEKTPH+46v3NGDNZ29eGdHD/7hg4dr89+wqw+buvcXdJ5DGWfzslc37sa7poxy/TxeXtu4G5f8ZCF+/OmTcd4xUwrur28wg68/8DrmnTEL7z2isHNe09WLPy/bgbnHTnF9JYqegWGs7uzFi2t24ZPvmeH6wbz0D2Wwakevq/l422X/cBa/XbwZsya24WxPR6nYsW8AZ3xvAWZOaMcfvnSmrxPs6hnE429tw81PrsLFJ0/DP51/lO/Z7ukfwk+fXYtTDhmHUS0pzJk13md6CaN/KIPlW/ehu28Ib2zeg69/+EjftQPDWRm2vBmfOX1mgda3p38Im3fvRzYnMHtKh2tCU7yxaQ/6BjNY09WLwyZ1+PyPgPMtrdvZi3FtTcjmREH+Qgis6erDY0u34bhpY9zAHS8bd/XjudVdmHvsQe5zEULg7hedSdenHjoeJ84Yi2ljW13LAQCs7uzBPS9uwDlHTcZgJotzj57ia9clG3fj6Xe6sLtvCMdPG4OPnzK9oG0efXMbZk1sx6ET230mX0UuZ99evdIQ0atCiMJRXpRrG1jInAHgX4UQc+Xf1wCAEOJ7pmvmzJkjFi9ebDrMMAzDaChFyDSyuewVALOJ6FAiagJwKYBHalwnhmEYxkPDOv6FEBki+iKAJwAkAdwlhDCvx8AwDMNUnYYVMgAghPgTgD/Vuh4MwzCMnkY2lzEMwzB1DgsZhmEYpmKwkGEYhmEqBgsZhmEYpmKwkGEYhmEqRsNOxowDEfUA2A5gr+GUMSHHAOAQABtDjtuur/TxsPrVc91Kzb+e62Y7Xs91i3Kc37n6q1upx3V1O1IIMUp3shUhxAHzD8BiAHeEHDcek8e7LMdt11f6uLF+9Vy3UvOv57qVoe78ztVh29Vz3cpwbwV1A7A4LL+wfweiuewPMY8BgG1DBtv1lT4eVr96rlup+ddz3WzH67luUY7zOxfveD0/17JuPHOgmcsWi5jr75Tj+kpTz/XjusWjnusG1Hf9uG7x0NWtlPoeaJrMHTW+vtLUc/24bvGo57oB9V0/rls8dHWLXd8DSpNhGIZhqsuBpskwDMMwVeSAFzJEdBcRdRLRW560E4hoIREtJaI/ENFomZ4montk+gq1h4089hcieoeIXpf/CndGqmzdmojo5zL9DSI6y3PNKTJ9NRHdQmXYh7WMdatEu80goqflM1pGRF+R6eOJaD4RrZL/j/Ncc41sn3eIaK4nvaxtV+a61bztiGiCPL+XiG4N5FXTtrPUraxtF6NuHyaiV2X7vEpE53jyqnW7hdWt+HaLG5Y2Uv4B+ACAkwG85Ul7BcAH5e8vALhe/v40gPvl7zYA6wHMkn//BcCcGtbtKgA/l78nA3gVQEL+/TKAM+BsePkYgI/UUd0q0W5TAZwsf48CsBLAMQD+A8DVMv1qAN+Xv48B8AaAZgCHAlgDIFmJtitz3eqh7doBnAngHwDcGsir1m0XVreytl2Mup0E4GD5+zgAW+qo3cLqVnS7HfCajBDiWQDdgeQjATwrf88H8HF1OoB2IkoBaAUwBGBfndTtGAAL5HWdcMIQ5xDRVACjhRALhfOW3AvgonqoW6l1CKnbNiHEa/J3D4AVAKYBuBDAPfK0e5BvhwvhDB4GhRDrAKwGcGol2q5cdSulDuWsnxCiTwjxPIABbz710HamulWCGHVbIoTYKtOXAWghouY6aTdt3eKWf8ALGQNvAbhA/r4EwAz5+0EAfQC2wZkR+59CCG9H+3OpQn67VBU3Rt3eAHAhEaWI6FAAp8hj0wBs9ly/WabVQ90UFWs3IpoFZ2S2CMAUIcQ2wPnw4GhVgNMemzyXqTaqaNuVWDdFrdvORD20nY2KtF2Mun0cwBIhxCDqr928dVMU1W4sZPR8AcBVRPQqHPVySKafCiAL4GA4potvENFh8tjfCiGOB/B++e+zVa7bXXBeyMUAbgbwIoAMHJU7SKVCCoutG1DBdiOiDgC/A/BVIUSYxmlqo4q1XRnqBtRH2xmz0KRVu+3CqEjbFVs3IjoWwPcB/G+VpDmtJu2mqRsQo91YyGgQQrwthDhPCHEKgF/DsYMDjk/mcSHEsDT7vABp9hFCbJH/9wD4FSpk0jDVTQiREUJ8TQhxohDiQgBjAayC07lP92QxHcDWYL41qlvF2o2I0nA+qF8KIR6SyTukOUKZczpl+mb4NSvVRhVpuzLVrV7azkQ9tJ2RSrRdsXUjoukAfg/gMiGE6mPqot0MdYvVbixkNKiICSJKAPgXAD+RhzYCOIcc2gGcDuBtaQaaKK9JA/gYHNNR1epGRG2yTiCiDwPICCGWSzW4h4hOl6rtZQAeroe6Vard5H3eCWCFEOKHnkOPAJgnf89Dvh0eAXCptIkfCmA2gJcr0XblqlsdtZ2WOmk7Uz5lb7ti60ZEYwE8CuAaIcQL6uR6aDdT3WK3WzAS4ED7B2fEvQ3AMJxRxOUAvgInAmMlgBuRn7TaAeC3cJxhywF8U6a3w4mYelMe+xFkBFAV6zYLwDtwnHpPApjpyWeOfBnWALhVXVPrulWw3c6EY2J4E8Dr8t9fAZgAJwBhlfx/vOea/yvb5x14onnK3Xblqludtd16OEEgvfJdOKaO2q6gbpVou2LrBmcQ1uc593UAk+uh3Ux1i9tuPOOfYRiGqRhsLmMYhmEqBgsZhmEYpmKwkGEYhmEqBgsZhmEYpmKwkGEYhmEqBgsZhqkTiOgfiOiyIs6fRZ5VsBmmHknVugIMwzgT3YQQP7GfyTCNBQsZhikTcvHBx+EsPngSnEmplwE4GsAP4Uzm3Qngc0KIbUT0FzjruL0PwCNENApArxDiP4noRDgrJrTBmZT3BSHEbiI6Bc5acP0Anq/e3TFMPNhcxjDl5UgAdwgh3g1nG4irAPwXgE8IZ023uwDc4Dl/rBDig0KIHwTyuRfAt2Q+SwFcK9N/DuDLQogzKnkTDFMuWJNhmPKySeTXe/oFgH+Gs/HTfLkqehLOcjyK3wQzIKIxcITPMzLpHgC/1aTfB+Aj5b8FhikfLGQYprwE12nqAbAsRPPoKyJv0uTPMHUNm8sYprwcQkRKoHwKwEsAJqk0IkrLfTqMCCH2AthNRO+XSZ8F8IwQYg+AvUR0pkz/2/JXn2HKC2syDFNeVgCYR0Q/hbO67X8BeALALdLclYKzcdsySz7zAPyEiNoArAXweZn+eQB3EVG/zJdh6hpehZlhyoSMLvujEOK4GleFYeoGNpcxDMMwFYM1GYZhGKZisCbDMAzDVAwWMgzDMEzFYCHDMAzDVAwWMgzDMEzFYCHD/P/26lgAAAAAYJC/9Sj2lUQAG8kAsAljWwnqynzq9wAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'] = sorted_data['inc'].astype(int)\n", "sorted_data['inc'].plot()" ] }, @@ -215,9 +2278,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 12, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXmcXGWV97+ntt67093pJJ3ubGQhhCUJhLC6IIZNUXEZQUeioCgjM+IyLvO+vjoqCuMoDuOgoqCAyjLIDCibYZE1EEIIZE86+9JJ72t1V3VVPe8fdW+luru6+/ZSa5/v51Ofqjx17+3n3lTV757lOUeMMSiKoiiKE1zpnoCiKIqSPahoKIqiKI5R0VAURVEco6KhKIqiOEZFQ1EURXGMioaiKIriGBUNRVEUxTEqGoqiKIpjVDQURVEUx3jSPYGJZurUqWbu3LnpnoaiKEpW8cYbbzQZY6pG2i7nRGPu3LmsX78+3dNQFEXJKkRkv5Pt1D2lKIqiOEZFQ1EURXGMioaiKIrimBFFQ0RmichzIrJNRLaIyJes8e+KyGER2Wg9Lovb51siUiciO0Tk4rjxM0Rkk/XebSIi1nieiDxgjb8mInPj9lktIrusx+qJPHlFURRldDgJhIeArxpjNohICfCGiKyx3rvVGPPv8RuLyBLgSuBkYCbwtIgsMsaEgV8A1wGvAo8DlwBPANcCrcaYBSJyJXAL8HERqQC+A6wAjPW3HzXGtI7vtBVFUZSxMKKlYYypN8ZssF53AtuAmmF2+SBwvzEmYIzZC9QBK0WkGig1xqw10c5P9wAfitvnbuv1Q8CFlhVyMbDGGNNiCcUaokKjKIqipIFRxTQst9Fy4DVr6AYReVtE7hKRcmusBjgYt9sha6zGej1wvN8+xpgQ0A5UDnMsRVEUJQ04Fg0RKQb+BNxojOkg6mqaDywD6oGf2Jsm2N0MMz7WfeLndp2IrBeR9Y2NjcOeh6Ioik0kYnhw/UF6+8LpnkrW4Eg0RMRLVDD+YIx5GMAYc8wYEzbGRIBfAyutzQ8Bs+J2rwWOWOO1Ccb77SMiHqAMaBnmWP0wxtxhjFlhjFlRVTXigkZFURQA3jzYytcfepuXdjWleypZg5PsKQHuBLYZY34aN14dt9kVwGbr9aPAlVZG1DxgIbDOGFMPdIrI2dYxrwYeidvHzoz6KPCsFfd4CrhIRMot99dF1piiKMq42dvkB6A3pJaGU5xkT50HfArYJCIbrbF/Aa4SkWVE3UX7gM8DGGO2iMiDwFaimVdftDKnAK4HfgcUEM2aesIavxO4V0TqiFoYV1rHahGR7wOvW9t9zxjTMrZTVRRF6c+BlqhohMKDvN7KEIwoGsaYl0gcW3h8mH1uAm5KML4eOCXBeC/wsSGOdRdw10jzVBRFGS0HmrsB6AtH0jyT7EFXhCuKMmmJWRoRtTScoqKhKMqkxRYNtTSco6KhKMqkpCsQoqkrCECfxjQco6KhKMqk5KBlZQCE1NJwjIqGoiiTkv3NcaKhMQ3HqGgoijIpibc0giG1NJyioqEoyqRkf0s3ZQVe3C4hFFHRcIqKhqIok5LDrT3UlhfgcYku7hsFKhqKokxKevrCFOV58LldBDUQ7hgVDUVRJiXBUIQ8jwuPWy2N0aCioSjKpCQYjuBzu/C4XRrTGAUqGoqiTEqCoQg+jwuvS3Rx3yhQ0VAUZVISEw2PS8uIjAIVDUVRJiXBkOWe0uypUaGioSjKpCQYjuD1uPC61dIYDSoaiqJMSgKWpaGiMTpUNBRFmZT0S7nV2lOOUdFQFGXSYYyJptx6XHhdammMBhUNRVEmHaGIwRisdRoaCB8NKhqKokw67Kq2Pg2EjxoVDUVRJh39RUMX940GFQ1FUSYddoFCn8eFx6VlREaDioaiKJOOmKWhMY1Ro6KhKMqkI97S0NLoo0NFQ1GUSYdtaWhp9NGjoqEoyqQjPhCupdFHh4qGoiiTDtsd5XVb7qmQioZTVDQURZl09AuEu7SMyGhQ0UghgVCYjt6+dE9DUSY9g9xTGtNwjIpGCvnJX3fyd79cm+5pKMqkJzBwcV8kgjEqHE5Q0UghO491sqepWz+cipJm7JhGnlVGxBgIq4vKESOKhojMEpHnRGSbiGwRkS9Z4xUiskZEdlnP5XH7fEtE6kRkh4hcHDd+hohsst67TUTEGs8TkQes8ddEZG7cPqutv7FLRFZP5MmnmqPtvQRDEboCoXRPRVEmNcdjGm48bgHQuIZDnFgaIeCrxpiTgLOBL4rIEuCbwDPGmIXAM9a/sd67EjgZuAS4XUTc1rF+AVwHLLQel1jj1wKtxpgFwK3ALdaxKoDvAGcBK4HvxItTttHQGQCgpTuY5pkoyuSmX+0pV/RnUIsWOmNE0TDG1BtjNlivO4FtQA3wQeBua7O7gQ9Zrz8I3G+MCRhj9gJ1wEoRqQZKjTFrTdQ/c8+AfexjPQRcaFkhFwNrjDEtxphWYA3HhSarCITCMbFoVtFQlLQSDIUBOxBuWRoaDHfEqGIalttoOfAaMN0YUw9RYQGmWZvVAAfjdjtkjdVYrweO99vHGBMC2oHKYY6VdTRaVgZAS5eKhqKkk/gyIl63WhqjwbFoiEgx8CfgRmNMx3CbJhgzw4yPdZ/4uV0nIutFZH1jY+MwU0sfxzriREMtDUVJK/HrNLyWpdGnMQ1HOBINEfESFYw/GGMetoaPWS4nrOcGa/wQMCtu91rgiDVem2C83z4i4gHKgJZhjtUPY8wdxpgVxpgVVVVVTk4p5Rzr6I29buoODLOloijJJmi5orxuwWPFNEJqaTjCSfaUAHcC24wxP41761HAzmZaDTwSN36llRE1j2jAe53lwuoUkbOtY149YB/7WB8FnrXiHk8BF4lIuRUAv8gayzps0RBR95SipJtgKILP7UJE8HrUPTUaPA62OQ/4FLBJRDZaY/8C3Aw8KCLXAgeAjwEYY7aIyIPAVqKZV180xoSt/a4HfgcUAE9YD4iK0r0iUkfUwrjSOlaLiHwfeN3a7nvGmJYxnmtaOdYRwOsWppXkq3tKUdJMMBTBZ4mF12W5pzQQ7ogRRcMY8xKJYwsAFw6xz03ATQnG1wOnJBjvxRKdBO/dBdw10jwznYaOXqaV5DO12KfZU4qSZoLhcEw0PG7bPaWi4QQnloYyARzr7GV6aR5lBV4auzSmoSjpxHZPAbGU2z4tj+4ILSOSIo629zK9NJ+KojyNaShKmol3T9ni0afl0R2hopEiGjoCTC/Np9JyT2n9KUVJH8HwcdHwuLSMyGhQ91QK6A6E6AyEmFaah1uEQCiCPximKE8vv6Kkg/7uKc2eGg1qaaSAJiuGUVWcR0WRD4BmdVEpGUxDZy9f+++36AmGR944CwnEZ09pGZFRoaKRAnr7oncwhT4PlcWWaOgCPyWDWbu7mYfeOMTW+uGKP2Qv/VJu1dIYFSoaKSBgFUfL87ioKMoDtJSIktm09/RZz7n5OQ2GI+QNsDS0jIgzVDRSQHwZ5kp1TylZQJu/r99zrtEvpqFlREaFikYKsFtL5nlclBV6AbRXuJLR2GJhWxy5RjAUibmltIzI6FDRSAHxlkahN9qPqjuQmwFGJTdos9xSuWpp9IW1jMhYUdFIAYF+DV9c5Hlc+IPa8lXJXDp6ct/SGFxGRC0NJ6hopIDj7qmolVGc56FbRUPJYHLePRW/uE97hI8KFY0UEB/TACjMc6t7Sslo2nrsQHhuJmwE4gLh9nNQLQ1HqGikgOAA0SjyeegOqKWhZC6x7KlctTRCx1NuY2VENKbhCBWNFBAfCAcoyvPgz9GVtkr2Y4zJ6ZiGMaafe8odEw21NJygopECBsY0Cn1ujWkoGUtPXzjmqmnPweypUMRgzHG3lIjgdYsu7nOIikYKGGRpqHtKyWBs11RlkY+2nr6cq8g88PsI0VIiWhrdGSoaKSAQCuNxScwM1kC4ksnYojG7spBwxNCdY67URKLhcYlmTzlERSMFxOeEQzTlVtdpKJmKvbBvbmVR9N85lkFlu94GWRoa03CEikYKCAwQjUKfRy0NJWOxg+CzKwqB3FsVHrM03CoaY0FFIwXEp/cBFPncBMOR2IdXUTIJWyTmTo2KRkeOZVAFErmn3KIptw5R0UgBgVC4v6VhdezL1QY3SnbTFrM0ivr9O1cY0tLQmIYjVDRSQLR2vzv27+K86OsujWsoGUibvw+vW5g5JT/271wiUUzD4xJdp+EQFY0UEF+7H6IxDQC/pt0qGUh7T5CyAh9TCnzWv3NLNPo0ED4uVDRSQCAUIc8bF9OwLI1cS2VUcoP2nj6mFHrJ97rweVyxbKpcobfP7qR53Pr3ukVLoztERSMFBAZYGkVqaSgZTJu/j7ICLyJCWYE351aF2yV87Js3iJZHD0XU0nCCikYKiFoaxz+gRVYgvEtFQ8lA2vx9TCmIdpicUuDNOfeUvUbKdhNDNKahloYzVDRSwOCYRlRAtGihkol0BUIU50d/UEvyPTnXmtheI1XkO34j5/NoTMMpKhopIBAKD4hpRL+QWrRQyUT8wXDsxibf66a3L7d+TGOWRl5/S0PXaThDRSMFBEMR8twJREPdU0oG0tsXpsAb/YxGRSO3LGLbwi/w9o9pqKXhDBWNFDCwjIj9YdVSIkqmYYzBHwxR4It+XgtyVDTyva5YAVGILvTTzn3OGFE0ROQuEWkQkc1xY98VkcMistF6XBb33rdEpE5EdojIxXHjZ4jIJuu920RErPE8EXnAGn9NRObG7bNaRHZZj9UTddKpZmAZEbdLKPC6tWihknEEwxEi5niQOM/ryjn3VHcgFMtgtCnKc6vl7xAnlsbvgEsSjN9qjFlmPR4HEJElwJXAydY+t4uIbQP+ArgOWGg97GNeC7QaYxYAtwK3WMeqAL4DnAWsBL4jIuWjPsMMYGCVW7A+pBoIVzIMu7RNvjc+ppFbn1N/MExhXLotQEm+l65eFQ0njCgaxpgXgBaHx/sgcL8xJmCM2QvUAStFpBooNcasNdGOLvcAH4rb527r9UPAhZYVcjGwxhjTYoxpBdaQWLwynkAo3G8hEUTjGnpno2QaPZZA2IHw3HRPDbY0ivM8dAfDhLX+1IiMJ6Zxg4i8bbmvbAugBjgYt80ha6zGej1wvN8+xpgQ0A5UDnOsrCJkmfsDLQ0tj65kIgODxPleF72hSE517/MHwxT4BloaunbKKWMVjV8A84FlQD3wE2tcEmxrhhkf6z79EJHrRGS9iKxvbGwcbt4p53h/8AHuKZ/GNJTMw3ZP2T+q+R434YjJqYVviWIatmh05tialGQwJtEwxhwzxoSNMRHg10RjDhC1BmbFbVoLHLHGaxOM99tHRDxAGVF32FDHSjSfO4wxK4wxK6qqqsZySkkjUWtJiOaIq3tKyTQGuaes595Q7ljF8etQbEryoyvg1dIYmTGJhhWjsLkCsDOrHgWutDKi5hENeK8zxtQDnSJythWvuBp4JG4fOzPqo8CzVtzjKeAiESm33F8XWWNZxXFLo/+HtFgD4UoG0jPAPWWXv8mluIY/GI6tlbI5bmmoaIyEZ6QNROQ+4N3AVBE5RDSj6d0isoyou2gf8HkAY8wWEXkQ2AqEgC8aY+xP2/VEM7EKgCesB8CdwL0iUkfUwrjSOlaLiHwfeN3a7nvGGKcB+YxhSEvD59GChUrG4R/knop+bnuDuZN2G12HMvAmzoppqGiMyIiiYYy5KsHwncNsfxNwU4Lx9cApCcZ7gY8Ncay7gLtGmmMmEwjZZZgHxzTUFFYyDduisC2NXHRPdQfC/epOwXH3VK7V2UoGuiI8ySTqRwxQnB9N8culrBQl+7EtDXtxX74nt9xTkYihpy/cr8ItqHtqNKhoJJmhRKMk30s4YrTSrZJR9PQNTLmNPudKP3v7/IoGLe7TlFunqGgkmeAQKbfF2lNDyUB6rDTwglj2lBXTCOVGTKM7dn79LY0Crxu3SzTl1gEqGknGLoI2UDQ0L1zJRHr6wrhdgtcdXSaVl2PuKX+CXhoAIkJxnkcD4Q5Q0UgygQT9iAFKY4E3/ZAqmYM/GKbQ68aqJxpzT+WMaAyI2cRTku/RmIYDVDSSjG1pDI5paOBNyTx6+8Lkx92F51vNw3JHNKLft4ExDYi6jDvVXTwiKhpJJtA3lHsqammoe0rJJAauli6IWRq5EtPov+I9ntJ8r34fHaCikWRGsjTUh6pkEj3BcL+OdjnnnrIsiUTuqeJ8jyamOEBFI8kMFdMoVveUkoH09PWvABtLuc0V0QjagXCNaYwVFY0kM5SlUezzIKLuKSWzGGhpuF2Cz5073fvsmMbAJkygouEUFY0kM1RMw+USin0ezZ5SMopEFWCjLV9zw9IYLqZRnKfd+5ygopFkguEIIuBxDW4Ponc2SqbR2xcetPAtl1q++gMhRI6XR4mnJN9DMBzJmXNNFioaSSYYiuBzu2J57/GU5HvpCqh7Sskc/MEwBd7+Pwu51PLVXofiGuImDrRKw0ioaCSZQCgyyDVlU6yWhpJhJCrml+/NnZhGdzBMYV7i4t66dsoZKhpJJhCK4EtgCoO6p5TMoycYjmVM2eR73TlTGt0fDCWMZwCU5Fnd+/Q7OSwqGkkmEAoPaWmU6GIiJYMIhSMEw5FBP6r5XnfOVLmNBvoTWxrFWg/OESoaSSY4jHtKLQ0lkxhYFt0mamnkhnvKHwwNKlZoE3NPaUxjWFQ0kkxvXzjWZ3kgKhpKJhETjYGWhscVW6Sa7XQFholp5NmlffQ7ORwqGkkmUd67TWm+l2A4EmsJqyjpxHZBDbQ0CnzunFkR3uYPMqXAm/A9WyztniJKYlQ0ksxwomE3YtI7GyUTsIVhUEzDkzspt63dQSqKfAnfsyvfajfN4VHRSDI9w4iGpvgpmYT9Y5k/KBCeGym3oXCEjt4QUwoTWxr2gr9uFY1hUdFIMv6+0JDZGloeXckkeu0SGwMD4TninmrriX7PygsTWxoul1Dgdat7agRUNJKMPxAeFFi0UUtDySRsS2NwINxNMBQhEjHpmNaE0eYPAlA+hHsKoq45dU8Nj4pGkrHLFiRioGg0dPRqUFxJG0PGNKzPbyDL025b/balkdg9BdHqt7myJiVZqGgkkUjEWGUZhs6egqh7KhiKsOrWF/iPp3elcoqKEsMuGz64YGH0ZyLbXVQt3ZalMYR7CqDQ66Fb3VPDoqKRROzSCwO/hDbx2VNvH2qjvaeP53c2AhDOcleAkn3YFq9tAdsU5Ej3PifuqQJ1T42IikYSiXUJS9DwBaJfznyvi7rGLl7b2wLA1voOWruDXPyzF7jtGbU6lNTR0RstG16coDQ6ZL9oOHJP+dQ9NRKJb4GVCWGoxVI2HreLC0+azpObj3JSdQk+j4tgKMJP1+ykrqGLPY1dqZyuMsnp7O2j2OcZVDbcdk9le9pta3cQn8c15PcRor3D2/w9KZxV9qGWRhKxfaNDpdwCXH7aTFq6g7xc18yHls3E53Hx+9f2A9GSB4qSKrp6Q7GiffHkSp/wVn+QikJfwt42NtHsKY1pDIeKRhLxD9Na0ubdJ1bFYhvvXFTF6bOnYKxwhjZoUlJJZ29oUDwD4rKnsl40+oZc2GejKbcjo6KRRHqGyHuPJ9/r5qKTpwOwcl4F586fCkBteQHdamkoKaQz0BdbcBpPLKaR5engrd3BYTOnwKqzpaIxLCOKhojcJSINIrI5bqxCRNaIyC7ruTzuvW+JSJ2I7BCRi+PGzxCRTdZ7t4llI4pInog8YI2/JiJz4/ZZbf2NXSKyeqJOOlU4sTQA/vniE/mvT5zOtJJ8rjl/Hn/87FmcPrtc204qKWUoS8OOAfQEszym4R+67pRNkc+Dvy+MMZq9OBROLI3fAZcMGPsm8IwxZiHwjPVvRGQJcCVwsrXP7SJi/2L+ArgOWGg97GNeC7QaYxYAtwK3WMeqAL4DnAWsBL4TL07ZgN9BTAOguqyA951WDUTTcM9dMJWiPC2brqSWzt5QzFUaz/FAeHbfgTtxTxX43IQjJusXMiaTEUXDGPMC0DJg+IPA3dbru4EPxY3fb4wJGGP2AnXAShGpBkqNMWtNVMLvGbCPfayHgAstK+RiYI0xpsUY0wqsYbB4ZTQ9Di2NRJTke+hWS0NJIVFLIzfdU5GIoc0/snuqMFYePXvPNdmMNaYx3RhTD2A9T7PGa4CDcdsdssZqrNcDx/vtY4wJAe1A5TDHyhq6xyEaRT4PPX1hQmG941FSQ2dvH6XDZU9l8Q9pZ2+IiBl+YR8c/676s9yqSiYTHQhPlMtmhhkf6z79/6jIdSKyXkTWNzY2OppoKuiJlWUYvWjYqY9apllJBcFQhEAoMkT2VPRnIptdNi32avAR3VPR89dKt0MzVtE4ZrmcsJ4brPFDwKy47WqBI9Z4bYLxfvuIiAcoI+oOG+pYgzDG3GGMWWGMWVFVVTXGU5p4/MEwbpfgc4/+Mhdbq8g1GK6kArs8f6KYhs/tQiS7Yxqt/pHrTgGx/uGJMhef3HyU//fI5kkfJB+raDwK2NlMq4FH4savtDKi5hENeK+zXFidInK2Fa+4esA+9rE+CjxrxT2eAi4SkXIrAH6RNZY12BVuh1tMNBTFVr/iLg2GKyngeN2pwXfiItE+E9ksGnbdKSeBcBjcve+RjYf5hz+8wT1r97Pz2OSu1OAk5fY+YC1woogcEpFrgZuBVSKyC1hl/RtjzBbgQWAr8CTwRWOMffWvB35DNDi+G3jCGr8TqBSROuArWJlYxpgW4PvA69bje9ZY1tATHLqXxkgUqaWhpBD7c5bIPQXRuEY2rwhvs+pOTRkxEG65p/qOf++CoQhff+htFs8oBWDt7qYkzTI7GLH2lDHmqiHeunCI7W8Cbkowvh44JcF4L/CxIY51F3DXSHPMVPx9YYoSmPtOsL+8KhpKKuiw3FOJLA2AfE92t3xtt7r2lRWMvCIc+lsaTV0BAqEInzpnDrf/rY5Xdjfz6fPmJW+yGY6uCE8i/kBo2OJow2G7pzTtVkkFQ5VFt8nPcveULRqJssPiiYlGXEyjsTMAQFVxHuecUMlre1uyvovheFDRSCL+4NANmEYi5p7SmIaSAiaDaBTnefCMkJRiu6fiixbGRKMkj3PmV9Le08fW+o7kTTbDUdFIIv6+scc0SixLo1MtDSUFdI7knvJmv3tqJNcUJF6n0dgVJxonRGvDvbqnOQmzzA60n0YS6QmGmFGaN6Z9bUtD3VNKKsh1S6PDoWjkeaLpxfELGZssS6Oy2Eeex8300jy21Xcmba6ZjloaSSTqnhqbLnvcLvK9Lg2EKymhKxAi3+vCO4T7piDLs6ecWhoiQpHP02+dRmNXgLICL3me6I3cnIoiDrR0J22umY6KRhLpGUdMA6ILrVQ0lFTQ2dsXS75IRLZbGm1+Z6IBVnn0vv4xjaqS4x6DOZWF7G/2T/gcswUVjSTSHQyNXzQ0EK6kgI7e0LCZRXmTJKYBgxsxNXYGqCruLxoNnYFJ2+FPRSNJRCKG3r5IrJbNWChSS0NJEUP10rDJ9hXh7T19lI2wGtymwDtANLr6WxqzK4sAONAyOa0NFY0kYft/1T2lZAOdvYm79tlks3uqty9MIBRxbGkU5Xn6WRFNnQGmxlsaFYUAk9ZFpaKRJJx27RuOknx1TympoWuIBkw2+V4XvVla5bbDXtg3BvdUdyBEdzA8KKYBcEBFQ5lIYv3Bx7giHKJ3PN2T1G+qpBYn7qlwxNCXhf1dnJYQsSnwHu8T3hS3RsNmSqGPsgIv+ydpBpWKRpLw9zlr9TocGghXUkVXIBTr4ZKIWCOmLHRRjVY04i2N+NXg8UzmDCoVjSRhL5Ya7os4EsV5Hl0RriSdSMTQHQxRMox7Ks9u+ToJRKOiKI+Gzl5C4chxS6O4v2jMrlDRUCaY1u5o/f6KEUoxD0dxnodgKEIwS33JSnbg7wtjDMNWZLbdrIEsTLsdrWgsnVVGb1+EHcc6Y5bG1JL+3+M5lYUcbuuhLxxh3d4Wzrv52UlTWkRFI0nEOoUVOfugJsL+EmspESWZ2J+v4UTDbvmaze6pKQ5F4/TZ5QBsONDGsY4ALoHKogHuqYoiwhHDkbYeXt/XwuG2HlbftY7ndjQkOmROoaKRJFq6ox/UihEa2Q9HsfbUUFLASA2YAPI92e+ecpo9VVtewNRiH28eaOWFXY2cVF2K29W/+6adQbW/2c/BFj9TCr0snF7Mdfes57G36yf2BDIMFY0k0eoPkudxjS97ygqiawaVkkzsZIuiYZI27GrN2bgqvM3fR0meZ9AP/1CICMtnl/PMtgbePtTOFctrBm0zx1rgt7+5mwMtfk6YWsQfP3c2S2un8I/3bWBvU+5mVqloJInW7iDlhb4x9Qe3OV7pNvvu7pTswXZPDZ89lb3uqY6ePsdWhs3y2VNo7+nDJXD50pmD3p9Wkkeex8X+Zj/7m/3MriikNN/LrR9fRsTAS7saJ2r6GYeKRpJo9QcpH4drCo77mCdrjRslNdgZesMt7svLcveU0yC4zfJZ0bjGeQumMr00f9D7Lpcwu6KQ3Y1d1Lf3MNtaJV5bXkB1WT6v7W0Z/8QzFBWNJNHSHaRiHEFwOL6aXAPhSjLpdiAa+Vmecjta0Vg2awonVZdyzTC9wOdUFrJubwsRA7Ms0RARVs6rYN3eFozJzZawKhpJotXfR/k40m3h+JdY3VPKWHh8Uz3/eN+b1DV0Dbtdl4PsKTumka0pt6MVjQKfmye+9A4uWDxtyG3mVBbRbS0CtGMcACvnVdDQGcjZdRwqGkmi1R8ct2gk6lesKE753zcP8+e3jnDpf7zAi8P42J1lT2VvTKOlO0hF8fi+i4mwM6iAmHsK4Kx5FQCsy1EXlYpGEgiFI7T39I07pmFbGl1qaShj4FBrDyvmlFNR5OPuV/YNuV1Xbwi3S8jzDP1zkK3uqXDE0OIPMnWc38VE2ELh87iYFldmZH5VMRVFPl7dm5uL/VQ0kkB7Tx/GQIXD+v1Dke+N9itWS0MZC4da/SyZWcrlp83k+Z2NtFkLTgfSHQhR5HMPm+l3XDSyyz3V6g9iDFQOKAMyEdguqVnlBbji0nlFhJNnlrJ7BLdgtqKikQSrU0rHAAAgAElEQVSOrwYf391Non7FiuKE9p4+OnpD1JYX8IFlM+kLG57cfDThtl2B8LC9NADcLsHndmWde6rFKudTmQT3VM2UAlzSP55hU12Wz5H23gn/m5mAikYSaPVHV6CON6YB0bUamj2ljJbDrT0A1JYXcmpNGXMrC3n0rSMJt+0K9MXWBA1HtOVrdomGXXBwYBmQicDncXHRkhm8+8SqQe9VlxXQ1BXIybpxYy/BqgyJfXcznhIiNkU+7amhjJ5DrdHMndryAkSEVUum87tX9hGJmH6uFIhm5w2XbmtT4HUTCGWXaDR3Jc/SAPjlp85IOD5zSj7GwLGO3lg6bq6glkYSsCvcjtc9BVCY179fsaI44VCcpWE/94UNzd2D4xqdgdCw6bY2+XHNibKF5pilkRzRGIrqsgIAjrT1pPTvpgIVjSTQ4h9/WXSbIp/2CVdGz6HWHgp9bsqtZAx7VfOxjsF+9u7A8F37bPK9rqwLhDd3B3FJtNteKpk5JXq963MwrqGikQTa/H3RYoXj6A9uM7DJvaI44VCrn1nlhbGMqBllQ4tGV29o2GKFNgVeN71Z5p5q6gpSUeRzXKxwoohZGu25Z2loTCMJREuITMydTaHPrdlTimN6rYZKB1t7qC0viI3PsCyNo0NYGk46TOZloXuqpTuQlCD4SBTleSgr8FLfppZGP0Rkn4hsEpGNIrLeGqsQkTUisst6Lo/b/lsiUiciO0Tk4rjxM6zj1InIbWLdHolInog8YI2/JiJzxzPfVGFXuJ0IivM8mj2lOOZbD2/i/FueZV9Tdz/RmFrswyVwbIC7xBhDVzDkKBBenJd9rtLmrmDSguAjUV2WT30OWhoT4Z66wBizzBizwvr3N4FnjDELgWesfyMiS4ArgZOBS4DbRcT23/wCuA5YaD0uscavBVqNMQuAW4FbJmC+SaehM8DUkom5uyn0eTQQrjiity/MU1uO0twdpKcvHAuCA3jcLqpK8gZZGv7gyK1ebaaV5NFgtT/NFpq7g0lZ2OeEmVMKOKKWhiM+CNxtvb4b+FDc+P3GmIAxZi9QB6wUkWqg1Biz1kTLQt4zYB/7WA8BF9pWSCZzpK2HmimDyymPheI8N93BUM5WzFQmjtf2tuAPhvnWpYtZWlvGWSdU9Ht/Rmk+Rzv6/+g7qXBrM600n6auAH3h7AmGN3UFUp45ZRNd4KeWxkAM8FcReUNErrPGphtj6gGsZ7tMZA1wMG7fQ9ZYjfV64Hi/fYwxIaAdqBznnJNKb1+Y5u5gLBA2XgrzPBiTnYXilNTyzLZj5HtdrD53Lo/ccD6n1U7p9/600vxB7iknvTRsZpRG1x7YC+YynUAoTGdvKG2iMXNKAW3+vqyLA43EeEXjPGPM6cClwBdF5J3DbJvIQjDDjA+3T/8Di1wnIutFZH1jY3o7Zh21vpTVZRNjaRT5tHufMjLGGJ7Z1sD5C6pidaIGErU0+ovGaCyN6aVRN8/RDEsjfWN/K19+YCM7jnb2Gz9eQiQ97in7NyDXrI1xiYYx5oj13AD8D7ASOGa5nLCeG6zNDwGz4navBY5Y47UJxvvtIyIeoAwYVG/YGHOHMWaFMWZFVdXgJf2pxP6AzJwyMZZGUaynRnYFIJXUsruxm8NtPbxnmP4PM8ryae/p61cKJNYf3JFo2Gm7mWVp3PXSXv7nzcNcdtuL/PG1A7HxZK8GHwk7pnSwJbf6aoxZNESkSERK7NfARcBm4FFgtbXZauAR6/WjwJVWRtQ8ogHvdZYLq1NEzrbiFVcP2Mc+1keBZ02GO/ftFLuJEg27p4aWElGGY+ex6F320lllQ25j/+jHWwpOemkM3D/RWo90EQpHeGFXI5edOoNz51fy3Ue3sK2+AyC2+n1qmkRj0fRigEEWULYzHktjOvCSiLwFrAMeM8Y8CdwMrBKRXcAq698YY7YADwJbgSeBLxpj7Fue64HfEA2O7waesMbvBCpFpA74ClYmViZjlw2YMPeUVUhOM6iU4bDLcM+bOrjiqk2itRr2zYgTS6OyyIfHJRklGhsOtNHZG+Ly02bys48vo6zQy433byQcMTRY80zHOg2IrkKvLstne46JxpgX9xlj9gBLE4w3AxcOsc9NwE0JxtcDpyQY7wU+NtY5poMj7b1UFPmG9CuPlqJYIya1NJSh2dPUTc2UgphlmogZZdEfz/gf/ePuqZE/ry6XMC1B2m46eW5HAx6XcN7CqZTme/naRYv4xp82UdfQxa6GLnweV7/1Kqlm8YySmOWTK2gZkQmmvr1nwqwMIFbewa+BcGUYdjd2cULV0FYGJHZPHesI4HaJ4zpp08vyacigmMbfdjRyxpxySq1+IMtnR9cSbznSzrb6DhZNL8bjTt/P3OLqUuoaunKqRLqKxgRT39Y7Yem2cPwOUGMaylAYY9jd0MX8quJhtyvJ9zK12MfuxuMd5Q639TCjNN/xD+v0ksEZWOniv56rY1t9B5eeMiM2dsLUIvI8LrYc6WD70U4WzyhN4wyjlkYoYvpd82xHRWOCOdI+cQv74LilkW3ZU+FIRucr5BQNnQG6g+ERLQ2AE2eU9AvMHm7toWYU7psZZYPXeqSDP791hB8/tYMPLZvJ1efMjY173C4WV5fy4q5GGjsDLJ5Rkr5JAkuqo6K1/WjuuKhUNCaQzt4+OntDVE9Q5hRE+2lA9gTCjTH865+3cNYPn86pL0omYwfBR7I0AE6cXsrOY11ELFE/3NZD7Sg+r9NK8+gMhNJ+E/PUlqNUl+Xz448tHdRU6uSZ0XMEOKk6vZbGvKlF+NwuttXnTjBcRWMCqZ/ghX0AeR43Xrek/UvqlF+9sIffvryPjt4Qq+9aF+sgpySP3U3dgDPRWDyjhJ6+MAda/PSFI9S3j9LSyJC0253HOjl5ZineBG61JXFCkW5Lw+N2sXB6cU4Fw1U0JpB91pc3vlDcRFDoy45Ktw2dvfzbk9u57NQZPHrDebT6+/jdy/vSPa2cZ3dDF0U+d2zF9nAssn5Etx/t5Gh7LxEDNaOwNKYPU2LdCb98fjcv7Woa0742wVCEPY3dLJqeWBBOnhkVjWkleWlbDR7PkupSthzpyJn6cSoa4+Boey/t/r7YvzcebMPjktiHdqIo8rnpzgL31LPbGogYuOGChSyeUcoJU4vY19yd7mnlFIda/dzwxw08vOEQR9t7eW57A395+wiLZpTgpJbnounFiETv1A9ba4pGY2nMtdaB1DWMPrDb0h3klie3c++r+0a9bzx7m7oJRQwnDmFFLJ5RikuimUuZwKm1ZbR0BzmSAbGgiUCbMI2DT/z6VWrKC7j32rMA2HCglSUzSydsjYZNeZEvK4rErdl6jJopBZxUHf0yz64oVNGYYH7z4l7+8nY9f3m7Pja2YFoxN3/4NEf7F/o8zK4oZMfRzpiFMRrLeGZZPjNK81m/r7VfANoJL+5qxBjY0zi+z8QOa/X7UJZGgc/N1efM5fQ55QnfTzWn1ERX6W861D4qqy5TUdEYIw0dvexp6mZPUze7jnUyb2oRbx9q52Nn1I688yiZU1nI9gwPpPmDIV6qa+KqlbNjd7yzKwp5YVcjxhhHd8HK8PQEw/xpwyE+sHQmV62czd6mbvI8Lt53WvWoblROnF7C9qMdsTv10cTgRIQz5pTzxv7WUc//2e3RMnT7m/2EI2bMLVh3Hu3E7ZJhs8W++4GTx3TsZLCkuhS3S9h8uJ1L4tKDsxV1T42RDQeOf2l++8o+dh7rwh8MxxYXTSRzK4s40OInlMF9DF7c1UQgFOGiJdNjY7MrC+nti9CYBVZSNvDnt4/Q2Rvik2fN5pz5lXzirNl85IzaUVu2J1WXsrepm9f3tVBVkjfq/U+fU87htp5RVbsNRwzP72yk0OcmGI5wuHXslV93WDdpeZ6JteiTRb7XzcJpxWw63J7uqUwIKhpjZMOBNnxuF1csr+HhDYd45K3DAJyeJNEIRUxGdwFbs/UYpfkezpx3vPHPrIrcrPKZLu5fd4AF04pZOa9i5I2H4e/OnIXH7eLFXU1jcpessNw+8TdOI7HxYCtt/r6YJb67aeyL3XYe6+TEIVxTmcopNWVsPtyeE8FwFY0xsmF/K6fUlHLDexaQ53Hzq+f3UFnkY1bFxPss51RGf3wzNT4Qjhie3d7ABYun9UuBnG2JxgEVjXHT0NHLhgNtfGjZzHG7+mqmFHD12XOir8dQlykat3Oxfp9z0XhueyNul7D63LkA7B1jXKM7EOJAi5+F00dOL84kTqsto7k7GEvLz2ZUNMZAMBTh7cPtnD67nPlVxTz8D+dywtQi3rN4WlJ893bl0kwVjTf2t9LSHWRVnGsKoj9OIlEftjI+nt4WjQesWjIxPvF/uGABUwq9LB7DHbvX7eK02im8MQpL47kdDZwxu5x5U4sozfewZ4yWxvr9rRhDUtzAycT2QHzqztd4bnvDCFtnNioaY2DLkXaCoQhnWGb6/Kpinvnqu7j5I84yWEZLVUkeBV43+5oy88d3zdajeN3Cuxb1b4CV73UzozQ/5ZbG0fZefrpmZ04ViXt62zFmVRTEejSMl4oiH8//8wVc/+75Y9r/tJoyttd3OIqzHevoZcuRDi6wbqpOqCoecwbVK3VNeN3CmXOzSzROqSnj9k+eTsTAlx/cmNVldlQ0xsCD6w/hcQlnxH1wRWTM2SAjISLMqSxkfwZaGsYY1mw9xjnzp1JiVRqNZ3ZFYcpjGr99ZS+3PbOLB9YfHHnjLMDOTHvvSdMn1JItK/COuQLsSdWlBEIRR9bv33ZE76wvWBy9qTihqmjMovHy7iaWzy4ftgR8pnLZqdXc+N6FtPn72HIke4PiKhqjZPPhdu5//QBXnzOXaSUTVy5kJOZWFrE3A0RjW30Hr9QdX9Fb19DFvmb/INeUzeyKwnFbGsYYehwubjTG8Ji1huG/nq3r19o0G0gUKH1pVxPBUIRVJyW+xunArum01UEq+HPbG6kuy48Fr+dXFXO0o5fO3r4R9uxPmz/IliMdnDd/6ugnnCGctyA69xfHuSo+nahojAJjDN/7y1bKC3186b0LU/q3504t4mCLP+1m7ff/spXP3bM+1hTqr1uPAQz5gza7opBjHQH8Yyztbozhhj++yUU/e96RAGw63M6h1h6uWF7D0Y5e7lt3YMR9MoV7X93Pe37yPK1Wm1Kb1/a2kOdxsWLu+LKmJpIF04rxuoWtR4avqdTSHeSFXY0x1xTAWVb2178/tWNUf/PVPc0YA+cuqBzbpDOAqcV5nFRdOu5SKulERWMUPL+zkXV7W/jyexdSVjDYFZNM5lYW0hc248pvHy/hiOGtg210B8P8z5vRFOM1W49xWm0ZM4ZYIGavhn3zQNuY/uaTm4/y2KZ6Drb08PtX94+4/WOb6vG4hO9cvoQz55bzu1f2xSq6ZjLtPX38+Mnt7G3q5mdP7+z33vr9rSytnYLPkzlfV5/HxYJpI3el+7ll7X3GypoCWDG3gs+eP4+71+7nyc1HHf/Nl+uaKfS5WVo7ZazTzgjesXAqb+xvHfONVLrJnE9hhmOM4Sd/3UlteQEfP3N2yv++XRLhuR3py7zY3dhFdzCM2yX84dX9NHT0svFg27BukzPnVeB2CWt3Nzv6G8YY1u1tIRwxdPb28f8e3cKS6lLOnV/JL/62mzte2M2/PbmdXccGu0UCoTB/eauecxdMZUqhj0+cNZv9zX5e3ePsb6eTO1/aS0dviPMXTOX3rx2InV9PMMyWw+394meZwknVw4vGwRY/v391Px87YxYLB2Rpff2SxSyaXsx/PVfn+O+9vLuJlfMqMko8x8L5C6YSDEey1kWV3Vc/hfzl7Xo2HW7nSxcuTMuHdtH0EhbPKOGRjYdT/rdtNlrWwjXnzWX70U6uu/cNAFadPLRoFOd5OLWmjLUOf7if29HA3/1qLbc9s4tfPb+Hxs4AP/zwqXzt4hNp7g7yw8e386sX9rDq1hcG/eD85sW9HG7r4drz5wFw6SnVlOZ7uO/1zA2It/f08aMntvHL53dz2akzuO2q5RT63NxqWRtvHWojFDGxBXWZxJLqUho6AzQnWPHvD4b4x/vexOWCG1cNduX6PC4+cnotmw63O0qUONrey57G7qyOZ9icdUIFsysKueWJ7f1cruGI4b51B2KFJDMVFQ0H1Lf38O1HNnNqTRlXLK9J2zw+sGwmGw60cSBN6x7ePNhGab6HL69axAeWzmRfczen1JSOuDr33PmVUbeWg/Luv381GoO4/W91/OalPXxg6UyWzZrC6bPLefgfzuWlb1zAun+5kPedWs2//3UHL+5qBKL/Rz9/to6LlkyPpf7me918+PRantp8NGMzz2744wbueGEPl50yg3/9wClUFPn41NlzeHLzUfY1dcdqPJ2RgaJhB8Ov+d3rfO/PW2PjkYjhn+57k7cPtfEfVy4fsv3xZadWA/DE5vqE78fzyu7oXfk587M3nmGT53Fz0xWnsKepm1ue3E4gFCYUjvDP//0W33p4E9/+381jOm5dQ2fsOiUTFQ0HfO2/36IvFOG2q5antUn9B5bOBKI1iNLBWwfbWDprCoU+D7ddtZw3v72KP99w/ohpoOfMryQUMby+r2XY7Q62+HluRwOfOGs2xXkewhHDP198Yuz902eXU1teSGVxHj/+2GksnFbMjfdvpLU7yK1rdhKOGL79/iX9jvmZ8+ZSmOfmY79cm9CllU4eeuMQL+5q4ruXn8zPrlxOVUm098Onz5uLx+Xitmd38cLORhZMK2ZKoS/Nsx3M0llTWDGnnFZ/H3e9vDdmMfxh3QGe3tbAt9+/hItPHnox4qyKQk6tKeOxTSPHNV6ua6a80NuvwVI2846FVfzdilp++/I+Vnz/aU757lM8/OZhTqst49ntDaNKyTXGcO/afbzvtpf4ziNbkh7DU9EYgbqGTl6ua+bLqxbFVmani9ryQs6aV8G9a/enPJW0Jxhmx7FOls06HoQUEUfrBlbMqcDrFu5+Zd+waZb3rTuAAF+8YAG/+8xK7vjUilj9qoHYwtXW08dXHtzInzYc5pNnzx60/ZzKIh78/DkY4HP3rM+YFNzW7iA/eGwbZ84t51NWSQ+baSX5fOSMGh7ecJjX9rZwXobeXRfneXjo+nP5w2ejrQH+8nY9h9t6uPnxbZy/YCqfjgt+D8Wlp87grYNtw8adjDG8sruJc+ZXDmrtms3c/OHTuPfalbzvtGr+/qw5/PrqFdx77VmU5Hm4/W+7HR/n58/W8e1HtnD2CZX84XNnJf0aqWiMwBPWXdDl1l1+uvnyqkUc7ejld6/sS+nf3XS4nXDE9BMNpxT43PzzxSfy/M5G3nfbSxxJ4LNt7Q5y79r9rFoynZopBSydNYULFk8b9riLZ5Ty6XPn8tyORnxuF//w7gUJt1s0vYSffXwZ+5r93D6KwGsy+eXzu+no7eMHHzo14Zf86xcv5odXnMp9nzub//O+JQmOkDnMqihk+ewpPLLxMDfe/yYG+NGHT3V0Q/HxFbM4oaqI1Xet43/fPJxwncqOY53Ut/dybg7EM+JxuYR3LKzi5o+cxv99/xJWLZlOWYGXq8+dw+Ob6h2Vn39k42F+smYnVyyv4XefOTMla8dUNEbgic1HOWNOeazNZbo5+4RKLjixitufq6PNHxx5hwnirYPRIPjSMYgGwHXvnM8Dnz+H1u4g1969njte2M0H/+tlvnDvGzzw+gF+/lwdXcEQX1l14sgHi+PG9y5k0fRibnzvwph7JxHnLZjKFctr+MXzu2NtedPF0fao6F+xrGbI7nPlRdHsr3PmV2ZFttDlp81k+9FOXt/Xyo8+fOqQFuJAKovzeOgL53JSdSk3PrCRj/zilUGf6/vXHcTndsViILnO9e9eQHVpPt/809vDlsJ5dvsxvvrgW6ycV8HNH3Em0hNB5n8a08iBZj9b6zu4NMMap3z9ksV0BkL8YhQm7HjZeLCN2vICpo6j5/KZcyv4+SdPZ+exTn74+HaMMWyt7+Abf9rEnS/t5cPLa4f8ER2KknwvT934Tj7/rpFrKH3rssWICL9+cc9YT2FI/MEQT2yqH7EWU29fmG8+/DYRY/jyqkUTPo908f7TqinJ83DdO0/gg8tGlyxSUeTjoS+cw48+fCobDrRx79rj63F6+6Jrgi46eToVRZkX10kGxXkefnDFKexq6OKGP25I2Lfkld1NfOH3GzipupTfrF6R0t4i2VfAJcnUNXSyrb6T959Wzd1r9wEMG8xLBydVl3LF8hp++8o+Vp87l5kpaCG58WAby2ePf1HVuxZVce81K/G4XaycV4Exhic3H+V/3jzM1y4e24+o0zusaSX5fHh5DQ+9cYivrFpE5TgEMJ5gKMLn732DF3c1cfnSmfzs48sS1iHrCoS47p71vLK7mZuuOMXx3Xg2MK00n9f+z4Vjrgnlcbu4auVsHt9Uz33rDnD9u+fjcbt4cvNR2nv6+MTK1K+NSifvWTydb1yymFuf3smFP/kb/3ThQi5fOpOIMWyv7+Sf7n+TORWF3H3NSkoT1HxLJioacby6p5nP3b2ezkCI+18/wMt1zVy1cnBwNRP4yqpF/OWter764Fv8n/edFFt5nQwaOns53NbDZ86bOyHHO3fBcd+0iHDpqdVcmiLXw2ffMY/7Xz/I3a/s4ysXjc4VNpC6hk5+8ted7G3qZvvRTi4+eTp/fusIXrfw448u5Rd/q2NbfSfnLqjkhKnF3PLkdjYdbufWjy/liuUT3xY43UxEEcG/P3sOn7/3DZ7d3kBNeQE/eGwrJ0wt4uwTMjMZIJlc/+75XHbqDL7356386Int/OiJ7bH35lQW8ofPnpUW60tFw6KuoYur71rH7IpCPrpwKr99eR8r51XwrxnUazie2vJCvv3+k7j5ie28/z9f4u5rVg4qTT5RvHUwmv43liB4prFgWgmXnjKDnz9XR215IafPmcJTW47xxOZ6brhgoeMezt2BENfd+waNnQEWzyjhh1ecyifOms1tz+zip2t2svFgG3sau6ko8vHYpug6BJ/bxe2fPD3jLNdM4sLF06guy+eGP76JwVBVnMevV6/Iqayp0TCnsog7P30ma3c3s7+5GxGoLMrjzHkVKS9lZCO50H4wnhUrVpj169ePej9jDHe9vI8PL6+hvMjHmwdaWTi9hOK8zNbV9p4+LvnZC5xQVcQfPnu24/0OtvhZt7eF9540nbLC4T98P35qO798fg9b/vXiUfeTzkSiP/jrebnueJrn1GIfrf4+fvzR0/jQspohf6S6AiEe31TPIxsPs3Z3M3/47NmDFpz9+1M7+PlzdVz/7vl8/eIT2dfsZ39zN7XlhSyYll0d59LB2t3NPL0tWgjzmvPnjaklrTJ6ROQNY8yKEbdT0ch+bv9bHf/25A6euvGdjgLJj2+q5xsPvU1nIERxnoevXbSIT583L+G2L+1q4sYH3qRmSgGP3HD+RE89bfT2hfnfNw+T73Vzam0Z00vz+fRd61i/v5UTphbx92fP4SNn1Pa7m9t8uJ1/uu9N9jR1U5IXXRl/zfmJr9vR9t4hizgqSiaSU6IhIpcA/wG4gd8YY24eatvJKBqt3UHO/tEzvO+0an7ysaXDBoaf2XaMz96znmWzpnDjexfx25f38rcdjXz63Ll889LFMUtiw4FWbn+ujqe3NbBgWjG3f/J0Fo2hNWg2EQxFeHxTPXev3cebB9oo8Lr5wNKZBEJhXt/XyuG2HqpK8vjp3y3lvPlTJ63LRMlNckY0RMQN7ARWAYeA14GrjDFbE20/GUUD4Ad/2cpvXtrLx1fM4prz5zGnsjAmAJGIIWwMz25v4CsPbGReVRH//flzKfC5CUcMP3hsK799eR+15QV8aFkN9e29/GnDoVgdpC+8az4Fvux3S42GzYfbuWftPh7ZeISSfC8r55Vz5twKLl86c1xpx4qSqeSSaJwDfNcYc7H1728BGGN+lGj7ySoakYjh1qd38p/PRlc8i0B1aT6dvSE6AyHcLiEcMdRMKeCh688ZVETu5bombnlyO5sPtyMiXHPeXG587yKKMjymk2z6whE8LmflUhQlm3EqGtnwi1ADxNe2PgSclaa5ZCwul/DVi07kg8tmsrW+kz2NXRxo9lNa4KWswEswHOH02eW8a1FVwhXG5y2YyqM3nE84YugLR3Ii4D0ReNNYoFJRMpFsEI1Et3j9zCMRuQ64DmD27Mm1CGggC6aVsGDa2GMPbpfgdqlgKIqSmGy4jToEzIr7dy3Qrza4MeYOY8wKY8yKqqrkrFVQFEVRskM0XgcWisg8EfEBVwKPpnlOiqIok5KMd08ZY0IicgPwFNGU27uMMVvSPC1FUZRJScaLBoAx5nHg8XTPQ1EUZbKTDe4pRVEUJUNQ0VAURVEco6KhKIqiOEZFQ1EURXFMxpcRGS0i0gnsGOVuZUD7BE1hKtA0QceCiZ1bNhxPr9/4mMjrl+nnmsnXDjL/fOOPNxUoMsaMvNDNGJNTD2D9GPa5I51/P1Vzy5Lj6fXLkOuXBeeasdcuS873jrjXjs9d3VNR/pzuCQzDRM8t04830WT6+Wby9cv0c83kaweZf75jOl4uuqfWGweVGnP172c7ev3Gh16/sTOZr91ozj0XLY07Jvnfz3b0+o0PvX5jZzJfO8fnnnOWhqIoipI8ctHSUBRFUZKEisYIiMgsEXlORLaJyBYR+ZI1XiEia0Rkl/Vcbo2vEpE3RGST9fyeuGOdYY3XichtMgnawU3w9btJRA6KSFe6zifVTNT1E5FCEXlMRLZbx7k5neeVCib4s/ekiLxlHeeXVhvqyclEpnDl4gOoBk63XpcQ7Ve+BPg34JvW+DeBW6zXy4GZ1utTgMNxx1oHnEO0sdQTwKXpPr8su35nW8frSvd5Zdv1AwqBC6zXPuDFXP/8TfBnr9R6FuBPwJXpPr+0Xdd0TyDbHsAjwCqiCwirreDUHuQAAAO3SURBVLFqYEeCbQVoBvKsbbbHvXcV8Kt0n0+2XL8B45NGNJJx/az3/gP4XLrPJ9uuHeAlmqr68XSfT7oe6p4aBSIyl+jdyGvAdGNMPYD1PC3BLh8B3jTGBIj2Oj8U994ha2zSMM7rN+mZqOsnIlOAy4FnkjnfTGIirp2IPAU0AJ3AQ0mecsaiouEQESkmapbeaIzpcLD9ycAtwOftoQSbTZrUtQm4fpOaibp+IuIB7gNuM8bsScZcM42JunbGmIuJWiZ5wHsS7DopUNFwgIh4iX7o/mCMedgaPiYi1db71UTvQOzta4H/Aa42xuy2hg8R7W9uM6jXea4yQddv0jLB1+8OYJcx5mfJn3n6mejPnjGml2i76Q8me+6ZiorGCFgZTncC24wxP41761FgtfV6NVF/qW36PwZ8yxjzsr2xZQZ3isjZ1jGvtvfJZSbq+k1WJvL6icgPiBapuzHZ884EJuraiUhxnMh4gMuA7ck/gwwl3UGVTH8A5xN1I70NbLQelwGVRH3Cu6znCmv7/wt0x227EZhmvbcC2AzsBn6Otbgylx8TfP3+jajFFrGev5vu88uW60fUsjXAtrjxz6b7/LLk2k0HXreOswX4T8CT7vNL10NXhCuKoiiOUfeUoiiK4hgVDUVRFMUxKhqKoiiKY1Q0FEVRFMeoaCiKoiiOUdFQlBQjIl8QkatHsf1cEdmczDkpilM86Z6AokwmRMRjjPlluuehKGNFRUNRRolV/O5JosXvlhMtuX01cBLwU6AYaAI+bYypF5G/Aa8A5wGPikgJ0Uq9/y4iy4BfEi1dvhu4xhjTKiJnAHcBfuCl1J2dogyPuqcUZWycCNxhjDkN6AC+SHSl8EeNMfYP/k1x208xxrzLGPOTAce5B/iGdZxNwHes8d8C/2SMOSeZJ6Eoo0UtDUUZGwfN8fpEvwf+hWjjnjVWQ0Y3UB+3/QMDDyAiZUTF5Hlr6G7gvxOM3wtcOvGnoCijR0VDUcbGwPo7ncCWYSyD7lEcWxIcX1EyAnVPKcrYmC0itkBcBbwKVNljIuK1+jIMiTGmHWgVkXdYQ58CnjfGtAHtInK+Nf7JiZ++oowNtTQUZWxsA1aLyK+IVkv9T+Ap4DbLveQBfka0KupwrAZ+KSKFwB7gM9b4Z4C7RMRvHVdRMgKtcqsoo8TKnvqLMeaUNE9FUVKOuqcURVEUx6iloSiKojhGLQ1FURTFMSoaiqIoimNUNBRFURTHqGgoiqIojlHRUBRFURyjoqEoiqI45v8Dv6icVnQKuPYAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -252,10 +2338,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 13, + "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", @@ -274,7 +2358,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +2382,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": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAD8CAYAAABQFVIjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAG2JJREFUeJzt3X2Q1dWd5/H3BxubjEEFBMOD2myFZaLJro5daMp9ELMCmUmJ2TFZRjdSNVbhGFPlrFOrUGHLbGR2ZGsrzrBWHK1oxsQHdDOxZMcgQYSKM8sCzaqjaJjGDYOMxG4WFEwVPWn57h/3tN6+dje3u899/ryqbt17zz3n3MPh3vvt8/D7/RQRmJmZjdeEWjfAzMyagwOKmZll4YBiZmZZOKCYmVkWDihmZpaFA4qZmWXhgGJmZlk4oJiZWRYOKGZmlkVbrRtQTeecc050dHTUuhlmZg1l9+7dhyNi+qnytVRA6ejooKurq9bNMDNrKJL+vpx8nvIyM7MsHFDMzCwLBxQzM8vCAcXMzLJwQDEzsywcUOpEz7ETfPWB7fQcP1HrppiZjYkDSp1Yt6WbXfuPsO757lo3xcxsTFrqOJR6NH/1Rvr6T374/NEdB3h0xwHa2yawd80Xa9gyM7PR8Qilxl68YyHXXDyLSRML/xWTJk5g6cWzePHOhTVumZnZ6Dig1NiMMycxub2Nvv6TtLdNoK//JJPb25gxeVKtm2ZmNiqe8qoDh9/v44bLLuD6Befz+M4D9Hph3swakCKi1m2oms7OzvC5vMzMRkfS7ojoPFU+T3mZmVkWDihmZpaFA4qZmWXhgNJAfDS9mdUzB5QG4qPpzayeedtwA/DR9GbWCDxCaQA+mt7MGoEDSgPw0fRm1ggcUMpQD4vhA0fTP/31K7jhsgvofb+vZm0xMxuKj5Qvw+qnX+WxnQe4YcH5rPny5yrQMjOz+lXukfJelB+BF8PNzMrnKa8ReDHczKx8ZQUUSfslvSrpZUldKW2qpM2SutP9lKL8qyTtk7RX0uKi9EtTPfskrZOklN4u6cmUvkNSR1GZ5ek9uiUtL0qfm/J2p7Knj787BvNiuJlZ+UYzQlkYERcXzaOtBLZExDxgS3qOpAuBZcBFwBLgu5JOS2XuB1YA89JtSUq/CTgaEZ8G7gXWprqmAncBlwELgLuKAtda4N70/kdTHdl5MdzMrDxlLcpL2g90RsThorS9wJURcUjSTGBbRMyXtAogIv4k5dsEfAvYD2yNiN9M6b+Xyt88kCcitktqA34JTKcQmK6MiJtTmQeAbcB6oBf4VET0S/p8Kv/haGgoPn29mdno5T59fQA/lbRb0oqUdm5EHAJI9zNS+mzgraKyB1Pa7PS4NH1QmYjoB94Dpo1Q1zTg3ZS3tC4zM6uBcnd5XRERb0uaAWyW9PMR8mqItBghfSxlRqprcGMKAXAFwPnnnz9UFjMzy6CsEUpEvJ3ue4CnKaxnvJOmukj3PSn7QeC8ouJzgLdT+pwh0geVSVNeZwFHRqjrMHB2yltaV2nbH4yIzojonD59ejn/XDMzG4NTBhRJZ0iaPPAYWAS8BmwABnZdLQeeSY83AMvSzq25FBbfd6ZpseOSLk+7u24sKTNQ13XAC1FY3NkELJI0JS3GLwI2pde2pryl729mZjVQzpTXucDTaYdvG/B4RDwnaRfwlKSbgAPAVwAiYo+kp4DXgX7g1oj4INV1C/AXwCeAjekG8BDwQ0n7KIxMlqW6jki6G9iV8n07Io6kx3cC6yWtAV5KdZiZWY341CtmZjai3Lu8zMzMRuSAYmZmWTigmJlZFg4oZmaWhQOKmZll4YBiZmZZOKCYmVkWDihmZpaFA4qZmWXhgGJmZlk4oJiZWRYOKGZmloUDipmZZeGAYmZmWTigmJlZFg4oZmaWhQOKmZll4YBiZmZZOKCYmVkWDihmZpaFA4qZmWXhgGJmZlk4oJiZWRYOKGZmloUDipmZZeGAYmZmWTigmFlD6Tl2gq8+sJ2e4ydq3RQr4YBiZg1l3ZZudu0/wrrnu2vdFCvRVusGmJmVY/7qjfT1n/zw+aM7DvDojgO0t01g75ov1rBlNsAjFDPLppLTUS/esZBrLp7FpImFn61JEyew9OJZvHjnwuzvZWPjgGJm2VRyOmrGmZOY3N5GX/9J2tsm0Nd/ksntbcyYPCn7e9nYeMrLzMatWtNRh9/v44bLLuD6Befz+M4D9Hphvq4oImrdhqrp7OyMrq6uWjfDrOn0HDvBmp+8wU/3/JITvz7JpIkTWHzRp/jm73zGI4gmIGl3RHSeKl/ZU16STpP0kqS/Ss+nStosqTvdTynKu0rSPkl7JS0uSr9U0qvptXWSlNLbJT2Z0ndI6igqszy9R7ek5UXpc1Pe7lT29HL/LWaWl6ejDEa3hnIb8EbR85XAloiYB2xJz5F0IbAMuAhYAnxX0mmpzP3ACmBeui1J6TcBRyPi08C9wNpU11TgLuAyYAFwV1HgWgvcm97/aKrDzGpkYDrq6a9fwQ2XXUDv+321bpJVWVlTXpLmAI8AfwzcHhFfkrQXuDIiDkmaCWyLiPmSVgFExJ+kspuAbwH7ga0R8Zsp/fdS+ZsH8kTEdkltwC+B6RQC05URcXMq8wCwDVgP9AKfioh+SZ9P5T8cDQ3FU15mZqOXe8rrT4E7gJNFaedGxCGAdD8jpc8G3irKdzClzU6PS9MHlYmIfuA9YNoIdU0D3k15S+syM7MaOGVAkfQloCcidpdZp4ZIixHSx1JmpLoGN0ZaIalLUldvb+9QWaxJ+RQdZtVVzgjlCuAaSfspTDVdJelR4J001UW670n5DwLnFZWfA7yd0ucMkT6oTJryOgs4MkJdh4GzU97SugaJiAcjojMiOqdPn17GP9eahU/RYVZdpwwoEbEqIuZERAeFNY0XIuLfAxuAgV1Xy4Fn0uMNwLK0c2suhcX3nWla7Liky9PurhtLygzUdV16jwA2AYskTUmL8YuATem1rSlv6ftbi5u/eiMdK5/l0R0HiCgcE9Gx8lnmr95Y66aZNbXxHCl/D3C1pG7g6vSciNgDPAW8DjwH3BoRH6QytwDfA/YBbwID3/CHgGmS9gG3k3aMRcQR4G5gV7p9O6UB3AncnspMS3WY+RQdZjUyqiPlI2IbhV1WRMT/A74wTL4/prAjrDS9C/jsEOkngK8MU9fDwMNDpP9fCluJzQYp95iInmMn+MYTL3Hf9ZfU/fESjdRWa10+l5c1pXKOiWikNZZGaqu1Lp96xVpO6XmnBtTjadAbqa3WvLKfesWsWTTSGksjtdXMAcVaTiOdd6qR2mrm09dbS2qk06A3UluttXkNxczMRuQ1FDOzYfi0PJXhgGJmLcfbsCvDayhWl3wgn1VCtS5V3Ko8QrG65L8grRK8DbuyPEKxutKKf0F6NFY93oZdWR6hWE0Mtyjain9BejRWXb5UceV4hGKjluMv6uIf0TVf/tyH6a30F2QrjsbqwQNf+2j365prP3auWhsHBxQbteGCQTnK+RFtlQP5XrxjIWt+8gY/3fNLTvz6JJMmTmDxRZ/im7/zmVo3zWxMHFCsbDn+oi7nR7RV/oJspdGYtQavoVjZcqxv+Ed0sEaaz/fBgHYqHqFY2XIFg1aZ0ipHI43GxjPVaa3B5/KyUbn5h11MnzxpUDAo/lG05uNrsli55/JyQDGzEfUcOzHsulerTlW2Gp8c0syy8LqXlctrKGZ2Sl73snJ4ysvMzEbkKS8zM6sqBxQzM8vCAcXMzLJwQDEzsywcUMzMLAsHFDMzy8IBxWwcfMJEs484oJiNg6+2aPYRHylvNga+2qLZx3mEYjYGOa4NY9ZsHFCsIdV67cInTDT7uFMGFEmTJO2U9IqkPZL+c0qfKmmzpO50P6WozCpJ+yTtlbS4KP1SSa+m19ZJUkpvl/RkSt8hqaOozPL0Ht2Slhelz015u1PZ0/N0iTWCeli7aKSrLZpVwylPDpl+9M+IiPclTQT+GrgN+LfAkYi4R9JKYEpE3CnpQuAJYAEwC3ge+KcR8YGknans/wZ+AqyLiI2Svg78s4j4A0nLgC9HxL+TNBXoAjqBAHYDl0bEUUlPAT+OiPWS/hx4JSLuH+nf4pNDNj5f7Mms+rKdHDIK3k9PJ6ZbAEuBR1L6I8C16fFSYH1E9EXEL4B9wAJJM4EzI2J7FKLYD0rKDNT1I+ALKZAtBjZHxJGIOApsBpak165KeUvf35qY1y4qp9bTiNb4ylpDkXSapJeBHgo/8DuAcyPiEEC6n5GyzwbeKip+MKXNTo9L0weViYh+4D1g2gh1TQPeTXlL6ypt+wpJXZK6ent7y/nnWh3z2kXl1MM0ojW2srYNR8QHwMWSzgaelvTZEbJrqCpGSB9LmZHqGpwY8SDwIBSmvIbKY43FF3vKy1ugLZdRHYcSEe9K2gYsAd6RNDMiDqXprJ6U7SBwXlGxOcDbKX3OEOnFZQ5KagPOAo6k9CtLymwDDgNnS2pLo5TiuqzJPfC1j6Zy11w70t82Vo4X71g47DXjzUajnF1e09PIBEmfAP4N8HNgAzCw62o58Ex6vAFYlnZuzQXmATvTtNhxSZenNZAbS8oM1HUd8EJaZ9kELJI0Je0iWwRsSq9tTXlL39/MRqGa04iNtE7TSG2tF+WsocwEtkr6W2AXhTWUvwLuAa6W1A1cnZ4TEXuAp4DXgeeAW9OUGcAtwPcoLNS/CWxM6Q8B0yTtA24HVqa6jgB3p/fdBXw7pQHcCdyeykxLddSEP3jW6Kq1BbqR1mkaqa31wteUz2D106/y2M4D3LDgfNZ8+XPZ6zdrdI203buR2lot5W4bdkAZB3/wzMrTc+zEsOs09bZDr5HaWi3ZjkOx4Y3mmAhPi1mt1fIz2EjbvRuprfXGAWUcRvPB83ys1VqtP4ONdKqaRmprPfGU1zjd/MMupk+eNOiYiOJtrZ4Ws1rzZ9DGy2soQ6jFubwabT6259gJvvHES9x3/SV12T4bvUb7DDaTZvk+eQ2lTjTafGytp0Usv0b7DDaTVvs++YqNVdAIpwrx6TeaWyN8BptJq36fPOVlgKdFzHJqtu+Tp7xsVDwtYpZPq36fHFBazEjHInirpFk+rfh98pRXi/FpYsxstMqd8vKifIto1UXC8WqWbZ9m1eAprxbhS+eOTatt+zQbD49QWkSrLhKOlUd0ZqPnEUoLacVFwrHyiM5s9DxCaSG+dG75PKIzGz0HFLNh+Ohys9HxtmEzMxuRj5Q3M7OqckAxM7MsHFDMmoAvMV197vOPc0AxawI+ALP63Ocf50V5swbmy/tWXyv2uRflzVqAD8Csvpx93mzTZg4oZhVWyR8NH4BZfTn7vNmmzXxgo1mFFf9oVOKSAT4As/rG2+fNeq44r6GYVUgrzrXn0AqXDGi0SwR7DcWsxry+MTbNNg00lGadqvSUl1mFNOuPRqU06zTQcJpxqtIBxayCmvFHo1JevGPhsNNAzagZz/7tgGJWQc34o1EpHtE1Pq+hmFnd8EXg8qvmsS7e5WVm1sRWP/0qj+08wA0Lzh/ztvVsu7wknSdpq6Q3JO2RdFtKnypps6TudD+lqMwqSfsk7ZW0uCj9UkmvptfWSVJKb5f0ZErfIamjqMzy9B7dkpYXpc9NebtT2dPL7Rwzs0YwntHF/NUb6Vj5LI/uOEBEYZNDx8pnmb96YwVaWlDOlFc/8EcR8RngcuBWSRcCK4EtETEP2JKek15bBlwELAG+K+m0VNf9wApgXrotSek3AUcj4tPAvcDaVNdU4C7gMmABcFdR4FoL3Jve/2iqw8ysaYxnC3Uttq2fclE+Ig4Bh9Lj45LeAGYDS4ErU7ZHgG3AnSl9fUT0Ab+QtA9YIGk/cGZEbAeQ9APgWmBjKvOtVNePgPvS6GUxsDkijqQym4ElktYDVwHXF73/tygELDOzhpZjC3UtNjmMalE+TUVdAuwAzk3BZiDozEjZZgNvFRU7mNJmp8el6YPKREQ/8B4wbYS6pgHvpryldZW2eYWkLkldvb29o/nnmpnVRK7RRbU3OZS9bVjSJ4G/BP4wIo6l5Y8hsw6RFiOkj6XMSHUNTox4EHgQCovyQ+UxM6snuUYX1d62XtYIRdJECsHksYj4cUp+R9LM9PpMoCelHwTOKyo+B3g7pc8ZIn1QGUltwFnAkRHqOgycnfKW1mVm1vAacQv1KUcoaS3jIeCNiPhO0UsbgOXAPen+maL0xyV9B5hFYfF9Z0R8IOm4pMspTJndCPz3krq2A9cBL0RESNoE/JeihfhFwKr02taUd33J+5uZNbxGPCi2nBHKFcDXgKskvZxuv00hkFwtqRu4Oj0nIvYATwGvA88Bt0bEB6muW4DvAfuANyksyEMhYE1LC/i3k3aMpcX4u4Fd6fbtgQV6ChsAbk9lpqU6rA4020WDzKw8PrDRsstxIJWZ1Y9yD2z0ubwsm1Y7W6yZDeZzeVk2vv6HWWtzQLFsfLZYs9bmKS/Lytf/MGtdXpQ3M7MR+ZryZmZWVQ4oTcTHf5hZLTmgNJHxnOrazGy8vCjfBHz8h9lgPcdO8I0nXuK+6y/xLsMq8gilCfj4D7PBPFqvDY9QmoCP/zAr8Gi9tjxCaRKNeKprs9w8Wq8tj1CaRCOe6tosN4/Wa8sBxcyais/WUDs+Ut7MzEbkI+XNzKyqHFDMzCwLBxQzM8vCAcXMzLJwQDEzsywcUMzMLAsHFDMzy8IBxczMsnBAMTOzLBxQzMwsCwcUMzPLwgHFzMyycEAxM7MsHFDMzCwLBxQzM8vCAcXMzLJwQDEzsywcUMzMLItTBhRJD0vqkfRaUdpUSZsldaf7KUWvrZK0T9JeSYuL0i+V9Gp6bZ0kpfR2SU+m9B2SOorKLE/v0S1peVH63JS3O5U9ffxdYWZm41HOCOUvgCUlaSuBLRExD9iSniPpQmAZcFEq811Jp6Uy9wMrgHnpNlDnTcDRiPg0cC+wNtU1FbgLuAxYANxVFLjWAvem9z+a6jAzsxo6ZUCJiJ8BR0qSlwKPpMePANcWpa+PiL6I+AWwD1ggaSZwZkRsj4gAflBSZqCuHwFfSKOXxcDmiDgSEUeBzcCS9NpVKW/p+5uZWY2MdQ3l3Ig4BJDuZ6T02cBbRfkOprTZ6XFp+qAyEdEPvAdMG6GuacC7KW9pXR8jaYWkLkldvb29o/xnmplZuXIvymuItBghfSxlRqrr4y9EPBgRnRHROX369OGymZnZOI01oLyTprFI9z0p/SBwXlG+OcDbKX3OEOmDykhqA86iMMU2XF2HgbNT3tK6zMysRsYaUDYAA7uulgPPFKUvSzu35lJYfN+ZpsWOS7o8rYHcWFJmoK7rgBfSOssmYJGkKWkxfhGwKb22NeUtfX8zM6uRtlNlkPQEcCVwjqSDFHZe3QM8Jekm4ADwFYCI2CPpKeB1oB+4NSI+SFXdQmHH2CeAjekG8BDwQ0n7KIxMlqW6jki6G9iV8n07IgY2B9wJrJe0Bngp1WFmZjWkwh/8raGzszO6urpq3Qwzs4YiaXdEdJ4qn4+UNzOzLBxQzMwsCwcUM7MG1nPsBF99YDs9x0/UuikOKGZmjWzdlm527T/Cuue7a92UU+/yMjOz+jN/9Ub6+k9++PzRHQd4dMcB2tsmsHfNF2vSJo9QzMwa0It3LOSai2cxaWLhZ3zSxAksvXgWL965sGZtckAxM2tAM86cxOT2Nvr6T9LeNoG+/pNMbm9jxuRJNWuTp7zMzBrU4ff7uOGyC7h+wfk8vvMAvTVemPeBjWZmNiIf2GhmZlXlgGJmZlk4oJiZWRYOKGZmloUDipmZZeGAYmZmWbTUtmFJvcDfD/PyORQuL1zv3M68GqWd0DhtdTvzqod2XhAR00+VqaUCykgkdZWzz7rW3M68GqWd0DhtdTvzapR2gqe8zMwsEwcUMzPLwgHlIw/WugFlcjvzapR2QuO01e3Mq1Ha6TUUMzPLwyMUMzPLomkDiqSHJfVIeq0o7Z9L2i7pVUn/U9KZKX2ipEdS+huSVhWV2SZpr6SX021GDdt5uqTvp/RXJF1ZVObSlL5P0jpJytnOzG2tWJ9KOk/S1vT/uEfSbSl9qqTNkrrT/ZSiMqtSv+2VtLgovaJ9mrmtddOnkqal/O9Luq+kror1aeZ21lN/Xi1pd+q33ZKuKqqr4t/7UYmIprwB/wr4LeC1orRdwL9Oj38fuDs9vh5Ynx7/BrAf6EjPtwGdddLOW4Hvp8czgN3AhPR8J/B5QMBG4It13NaK9SkwE/it9Hgy8HfAhcB/BVam9JXA2vT4QuAVoB2YC7wJnFaNPs3c1nrq0zOAfwH8AXBfSV0V69PM7ayn/rwEmJUefxb4h2r051huTTtCiYifAUdKkucDP0uPNwO/O5AdOENSG/AJ4B+BY3XYzguBLalcD/Au0ClpJnBmRGyPwqfsB8C19djW3G0aoo2HIuL/pMfHgTeA2cBS4JGU7RE+6p+lFP6Y6IuIXwD7gAXV6NNcbc3ZphztjIhfRcRfA4Ou9lTpPs3VzkobQztfioi3U/oeYJKk9mp970ejaQPKMF4DrkmPvwKclx7/CPgVcAg4APy3iCj+4fx+Gvb+pyoNKYdr5yvAUkltkuYCl6bXZgMHi8ofTGnVMNq2Dqh4n0rqoPDX3Q7g3Ig4BIUvNIVRExT66a2iYgN9V9U+HWdbB9RLnw6nan06znYOqMf+/F3gpYjoo7bf+yG1WkD5feBWSbspDDX/MaUvAD4AZlGYSvgjSf8kvXZDRHwO+Jfp9rUatvNhCh+aLuBPgf8F9FMY7paq1va90bYVqtCnkj4J/CXwhxEx0mhzuL6rWp9maCvUV58OW8UQadn7NEM7oQ77U9JFwFrg5oGkIbLVdNtuSwWUiPh5RCyKiEuBJyjMQUNhDeW5iPh1mp75G9L0TET8Q7o/DjxOdaYYhmxnRPRHxH+IiIsjYilwNtBN4Yd7TlEVc4C3S+utk7ZWvE8lTaTwRX0sIn6ckt9JUwQDUy89Kf0gg0dOA31XlT7N1NZ669PhVLxPM7Wz7vpT0hzgaeDGiBj43arZ9344LRVQBnZqSJoArAb+PL10ALhKBWcAlwM/T9M156QyE4EvUZjiqUk7Jf1Gah+Srgb6I+L1NDw+LunyNDS/EXim0u0cS1sr3afp3/8Q8EZEfKfopQ3A8vR4OR/1zwZgWZqTngvMA3ZWo09ztbUO+3RIle7TXO2st/6UdDbwLLAqIv5mIHMtv/fDyr3KXy83Cn8tHwJ+TSGS3wTcRmFHxd8B9/DRgZ2fBP4HhQWv14H/GB/tAtkN/G167c9Iu2pq1M4OYC+FRbznKZwBdKCeTgof+jeB+wbK1FtbK92nFHbtRKr/5XT7bWAahU0C3el+alGZb6Z+20vRLplK92muttZpn+6nsIHj/fRZubDSfZqrnfXWnxT+UPtVUd6XgRnV+t6P5uYj5c3MLIuWmvIyM7PKcUAxM7MsHFDMzCwLBxQzM8vCAcXMzLJwQDEzsywcUMzMLAsHFDMzy+L/A9t1yP3IuE0bAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.plot(style='*')" ] @@ -314,9 +2421,57 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 16, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "2021 743449\n", + "2014 1600941\n", + "1991 1659249\n", + "1995 1840410\n", + "2020 2010315\n", + "2022 2060304\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": 16, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,9 +2486,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 17, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 17, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEKCAYAAAAyx7/DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAE6lJREFUeJzt3XuQZGV5x/Hvk10uC8MtWRl0QcZbjMjKbVCRiLNoeWFRo1KCgrqUuJgoEFwrWSmVqCHiBVNqvGS9lyJTCiQqGC+BDCoaZBfUFVfRglVYBAQEXSDCypM/3jOhazOz07tMT/fb/f1UTVV3n9N93uftnl+f855z+kRmIkmqx590uwGSpK1jcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqM78TL7pw4cIcGRnpxEvP6O6772bnnXfuyrJ7xaD3waDXD/ZBjfWvWbPmtsx8WDvzdiS4R0ZGWL16dSdeekYTExOMjY11Zdm9YtD7YNDrB/ugxvoj4pftzutQiSRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyHTkBR5rJyMqLO/baKxZvYtk0r7/+7KUdW640V1zjlqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMga3JFXG4JakyhjcklQZg1uSKmNwS1JlDG5JqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkirTVnBHxOkRcU1E/DgizouIHTvdMEnS1GYM7ohYBJwKjGbm/sA84LhON0ySNLV2h0rmAwsiYj6wE3BT55okSdqSyMyZZ4o4DTgLuBf4RmYeP8U8y4HlAMPDw4eMj4/PclPbs3HjRoaGhrqy7F5RQx+s3XBXx157eAHccu/U0xYv2q1jy+0lNXwGOqnG+pcsWbImM0fbmXfG4I6IPYALgGOBO4EvAudn5ueme87o6GiuXr26/RbPoomJCcbGxrqy7F5RQx+MrLy4Y6+9YvEmzlk7f8pp689e2rHl9pIaPgOdVGP9EdF2cLczVPIs4PrM/E1m3g9cCDztoTRQkrTt2gnuXwFPjYidIiKAZwLrOtssSdJ0ZgzuzLwCOB+4CljbPGdVh9slSZrG1AOBm8nMM4EzO9wWSVIbPHNSkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMga3JFXG4JakyhjcklQZg1uSKtPWFXDUnzp5pXVJneMatyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMga3JFXG4JakyhjcklSZtoI7InaPiPMj4qcRsS4iDut0wyRJU2v3YsHvB76WmcdExPbATh1skyRpC2YM7ojYFTgCWAaQmfcB93W2WZKk6URmbnmGiAOBVcBPgAOANcBpmXn3ZvMtB5YDDA8PHzI+Pt6RBs9k48aNDA0NdWXZvaLdPli74a45aM3cG14At9w79bTFi3ab28Z0yaD/H9RY/5IlS9Zk5mg787YT3KPAfwOHZ+YVEfF+4HeZ+ZbpnjM6OpqrV6/emjbPmomJCcbGxrqy7F7Rbh+MrLy4843pghWLN3HO2qk3JtefvXSOW9Mdg/5/UGP9EdF2cLezc/JG4MbMvKK5fz5w8LY2TpL00MwY3Jl5M3BDRDy+eeiZlGETSVIXtHtUySnAuc0RJdcBJ3auSZKkLWkruDPzB0BbYy+SpM7yzElJqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqTLuXLpP6Qr9e2X5zKxZvYllT66Bc2X6QuMYtSZUxuCWpMga3JFXG4JakyhjcklQZg1uSKmNwS1JlDG5JqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVabt4I6IeRFxdURc1MkGSZK2bGvWuE8D1nWqIZKk9rQV3BGxN7AU+HhnmyNJmklk5swzRZwPvBPYBXhjZh49xTzLgeUAw8PDh4yPj89yU9uzceNGhoaGurLsXtFuH6zdcNcctGbuDS+AW+7tdiu6q7UPFi/arbuN6YIac2DJkiVrMnO0nXnnzzRDRBwN3JqZayJibLr5MnMVsApgdHQ0x8amnbWjJiYm6Naye0W7fbBs5cWdb0wXrFi8iXPWzvjR7mutfbD++LHuNqYL+j0H2hkqORx4QUSsB8aBIyPicx1tlSRpWjMGd2a+KTP3zswR4Djg0sw8oeMtkyRNyeO4JakyWzUQmJkTwERHWiJJaotr3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMga3JFXG4JakyhjcklSZwb4Udo8YmeWrra9YvKlvr+AutWPthru68j+w/uylc7Ic17glqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZUxuCWpMga3JFXG4JakyswY3BGxT0T8V0Ssi4hrIuK0uWiYJGlq7VwseBOwIjOviohdgDUR8c3M/EmH2yZJmsKMa9yZ+evMvKq5/XtgHbCo0w2TJE0tMrP9mSNGgG8B+2fm7zabthxYDjA8PHzI+Pj47LVyK2zcuJGhoaGtft7aDXd1oDXdMbwAbrm3263onkGvH3qjDxYv2q1ry771jru6Uv9DqXnJkiVrMnO0nXnbDu6IGAIuA87KzAu3NO/o6GiuXr26rdedbRMTE4yNjW3180ZWXjz7jemSFYs3cc7adkbB+tOg1w+90Qfrz17atWV/8NwvdaX+h1JzRLQd3G0dVRIR2wEXAOfOFNqSpM5q56iSAD4BrMvM93W+SZKkLWlnjftw4BXAkRHxg+bvqA63S5I0jRkHgTLzO0DMQVskSW3wzElJqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5IqY3BLUmUMbkmqTM9dCvuhXm19xeJNLOujK7ZL0uZc45akyhjcklQZg1uSKmNwS1JlDG5JqozBLUmVMbglqTIGtyRVxuCWpMoY3JJUGYNbkipjcEtSZQxuSaqMwS1JlTG4JakyBrckVcbglqTKGNySVBmDW5Iq01ZwR8RzI+JnEfGLiFjZ6UZJkqY3Y3BHxDzgQ8DzgP2Al0XEfp1umCRpau2scT8Z+EVmXpeZ9wHjwAs72yxJ0nQiM7c8Q8QxwHMz86Tm/iuAp2Tm6zebbzmwvLn7eOBns9/ctiwEbuvSsnvFoPfBoNcP9kGN9e+bmQ9rZ8b5bcwTUzz2/9I+M1cBq9pZaCdFxOrMHO12O7pp0Ptg0OsH+6Df629nqORGYJ+W+3sDN3WmOZKkmbQT3FcCj4uIR0XE9sBxwJc72yxJ0nRmHCrJzE0R8Xrg68A84JOZeU3HW7btuj5c0wMGvQ8GvX6wD/q6/hl3TkqSeotnTkpSZQxuSaqMwS1JlTG4NxMRiyJiUbfb0S0R8eiIOD0ijux2W7rB+ge7fqijDwzuRkSMRMRlwNeA90TE07vdprkWEX8JfJPymzSvjYi/7nKT5pT1D3b9UE8fDHRwR8SOLXcPBq7MzMWUQx9PjYjFzXxTnT1avYg4MiIe1dwO4EjgzMx8DXAOcFREjLVM7yvWP9j1Q719MHDBHRG7RsRHI+Ja4L0RsW8z6UXAr5rb48AvgJMmnzbHzeyoiNgvIn4E/APwqYg4MstxofsBewFk5hXAd4ETJ5/WjbZ2gvUPdv1Qfx8MXHADzwV2pLxB9wFvjYgFlM2j5wNk5h+A84GnN/cf6E5TZ0dE7B0Ru7Y8dCxwQWYeQfmSenlEPA74PE0fNP4N2D8idqi5D6x/sOuH/uuDvgzuKOZHxKsj4tsRcVpEPKaZ/FjgvszcBPwz8FvgBOAbwMMj4k+b+a4FboiIw+a8gFkSEU+IiK8C3wHeHhGTP8f7P8BOze0vADcDSylrF3/WshVyB+VXHg+Yu1bPHusf7Pqhf/ugL4O72eR5BvBK4N3ADsDHmsk3A7c236A3UN6Ux1DeyJ/w4E/Tbgfc3jxejYjYueXugcCNmTkCXAq8t3n8DuAPEbFLZt4B/Bx4BKXW7wJvaObbHvgjsL7zLZ8d1j/Y9cNg9EFfBHdEHBYR74qIZc39AJ4AfC0zv5KZ7wb2jYinARso37RPaJ6+DhhqHvsXys6I51NCfxj44ZwWsw0iYo+I+HREXAmcHREPa/pgMXB5RERmfhm4MyKWUrYmdmmm09zfE3iAshWyZ0R8DDgP2JSZt851TVvD+ge7fhi8Pqg+uCPiicBHgN8DL42IN1DqWgT8vmVP8KeBl1OCeBPwtObxqyh7ku/JzG8BK4FlwOHAOzLzgV7amzyNIyg1HUXZgXIGsCvlR8H2ygd/kOYzlD74PqW/ngeQmd9rXmN+Zq4DTgauAf4pM0+k91n/YNcPg9YHmVnNH2XN+CTK5s/85rH3Aac1t0eBDwDHAM8Cvt7y3H0om0xQgvpqypV6DgK+BDy8Zd7odq3T1D+P8oG6jDKks7B5/AvAqc3tRwFnN9MPpYztzWvpv980r7OIsrXxeuBTwIeBnbtdo/Vbv30w8181a9wRcQBlB+ILgTOBNzeTNlCuiwnlG/Jy4CXAJcBeEfGkiNguy3j2hoh4emZeSvnZx3cBFwLnZeavJ5eVzTvcg44GXgC8DTiMMn4P5YiYyS2IG4BvA8/LzCspax9LADJzI3AFcGhmbgBeQRkOuhl4c2bePUd1bCvrH+z6wT4A2rt0WVdExEHA3Zl5bfPQk4FrM3NZRBwMnBURo8AE8JyI2Ckz74mIHwIvpRyL+XngNcAHIuJeYC1wffN6HwU+n5l3zV1VM2vG4jIiDqVs0n0buDjLIYp/DlyXmZdGxPWUMzyfDawBXhQRCzPztoj4OXB3RDwS+CBwQkTsSbl60e2UzUQyczWwes6L3ALrH+z6wT5oR8+tcUfEYyPiO5Qx6bdGxCubSQ8A65u156somziHAffw4KE8APdTNoX2oqxV/7h5rcuA2zLzRihr1T0c2kcAn6Ts4X4W8M5mlgeAayNiQWZeT+mDJ1HG6m6iHJsKZS/4PMr7ewGlH44HDgFWZQ8dj9oqIuY19T+Dsuk6aPUvaOofYwDffygnyA16H7Sl22M1wM7AYS33jwbe39x+CuXbcF/gVZTNo0XNtGMo49mT0y5rHt+RMkyysOU1DwK273at09S/E/BaHtw62A74W+B1zfQ9gB81NRxLGbsbaemrVZQrWh9N2aLYjTKG/9XWmoE/6XatW3j/T6L8c62g7FAamPqbtu0CXEy5uhTA6QNW/07N//AllJNiBq4Ptvavq2vcEXEGcB3w1YgYbh5+DuXYarKccvp94BTKMZj7UI65hjKWfSDlaJDPAL+NiM9Sdjr+DPi/sarMvDoz7+t8RVsnIvYCLgLGgM9Sdqa8mLIlsQkgM39L2Xl6KmUcb08ePJTxW5Tj1e/LzIuAT1DO+PwQZe/5/ZPLyh5cw2iOt72E8k/2MeDZlP0Th1LWrPq6/hYLKOcaPCYiFlI+4/Og/+uPiO0o+6aOAd6TmS9pJh00OU+/98E26fI37RhlM+fjwIrmsTdQxrMm59kPuKG5fRbw9pZpVwIHNbd3oBwKdGi3vw23ov4FwFNa7i+j7Gx5FfD9lscfAdzU3H4d5RTdPZrnfwV4ZMu8C+ei7bPYB7u33P47yj/n8YNSf9PmVwHvAd4CvJpyyvWVA1T/hcDxmz12LHDFoPTBVvdZl9+wyUN0juXBoY7dgTuBHVvmu5LyDbw78EXKptF/UL5Rd+h2Jz6E+mPyr7l/cEs/3E45/nRy3m9Ohjzwj5QjbG4H/r7bdcxCP+xK2Q9xC/D25v7twHA/19/yvp9IGS57MXBu89ht/V5/S21HU06AOYdysMFbKUOgdwB7DkIfbO1fV4dKMvOPzc3/BHaJiAMy807KuPbJLbNeBezSTDuFMhzy78DyLHuaq5SNlodOo6x9QBmfOx0gyu+n/BKYPGTxbZQtk0WZ+a45am7HZObvKENiT6XsVH4pZbjr5Cj6sv6W9/4oylDRJcDeEfFmyg735TAQ7/9FlKO9bqccnvdE4K8on4HX9vNnYFv1zFXeI+LDlPHqNzZHVfwNJcD3oJxYc1RL0PediNibMj53SmZeG+VHsZZTPsSLgB9kL57BNcsi4kDKl/b3KOOY+1MO4erL+iNiiDJMsgOl3r+gnDByBmVN/HH0cf2TJg/nbW4fQPnsX045Jb2vPwPbopeC+0DKUSLPpHyA76Gcfn4v8NHMvKaLzeu4KL+P8gzgTZQxzxspm4bHAj/Ncghk34uIfShfYC/LzNsj4gTgmsy8ustN64goF/P4V8pOtPMoh7GdkZnPbqb3df1TifLLfB8Hjs3MOwaxD2bSS8F9HOWQuHuAd1D2MPfPXuAZRMTlwKMpv0J2E/C2zPxRVxs1RyJiN8oX9sspO6NXAR/KzPu3+MQ+1Jww8mJgPDNv7nZ75kpE7ED5rfzJoZKPAB/O8vPL2kxPBHdEPIly+vn5lJ0zVf2U6kPVHBJ1JmWc73M1j9tvi4iYTxke+QOl/oF6/6GcfAQ8kL3wD9klEXEy5TDQzw7iZ2Br9ERwS5La13OnvEuStszglqTKGNySVBmDW5IqY3BLUmUMbkmqjMEtSZX5XwHKcVej17E2AAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.hist(xrot=20)" ] @@ -341,9 +2519,7 @@ { "cell_type": "code", "execution_count": null, - "metadata": { - "collapsed": true - }, + "metadata": {}, "outputs": [], "source": [] } @@ -364,7 +2540,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4, -- 2.18.1