diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/chickenpox.ipynb similarity index 100% rename from module3/exo2/exercice.ipynb rename to module3/exo2/chickenpox.ipynb diff --git a/module3/exo3/exercice.ipynb b/module3/exo3/exercice.ipynb index 0bbbe371b01e359e381e43239412d77bf53fb1fb..489fb5d634ff300e4ed9847633e6e4f835c29338 100644 --- a/module3/exo3/exercice.ipynb +++ b/module3/exo3/exercice.ipynb @@ -1,5 +1,2390 @@ { - "cells": [], + "cells": [ + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [], + "source": [ + "%matplotlib inline\n", + "import matplotlib.pyplot as plt\n", + "import pandas as pd\n", + "import isoweek" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + " data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "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
0202011310170493652.0109756.0154142.0166.0FRFrance
1202010310497796650.0113304.0159146.0172.0FRFrance
22020093110696102066.0119326.0168155.0181.0FRFrance
32020083143753133984.0153522.0218203.0233.0FRFrance
42020073183610172812.0194408.0279263.0295.0FRFrance
52020063206669195481.0217857.0314297.0331.0FRFrance
62020053187957177445.0198469.0285269.0301.0FRFrance
72020043122331113492.0131170.0186173.0199.0FRFrance
820200337841371330.085496.0119108.0130.0FRFrance
920200235361447654.059574.08172.090.0FRFrance
1020200133685031608.042092.05648.064.0FRFrance
1120195232813523220.033050.04336.050.0FRFrance
1220195132978625042.034530.04538.052.0FRFrance
1320195033422329156.039290.05244.060.0FRFrance
1420194932566221414.029910.03933.045.0FRFrance
1520194832236718055.026679.03427.041.0FRFrance
1620194731866914759.022579.02822.034.0FRFrance
1720194631603012567.019493.02419.029.0FRFrance
182019453101387160.013116.01510.020.0FRFrance
19201944378225010.010634.0128.016.0FRFrance
20201943394876448.012526.0149.019.0FRFrance
21201942377475243.010251.0128.016.0FRFrance
22201941371224720.09524.0117.015.0FRFrance
23201940385055784.011226.0139.017.0FRFrance
24201939370914462.09720.0117.015.0FRFrance
25201938348972891.06903.074.010.0FRFrance
26201937331721367.04977.052.08.0FRFrance
2720193632295728.03862.031.05.0FRFrance
28201935310102.02018.020.04.0FRFrance
2920193431672279.03065.031.05.0FRFrance
.................................
181619852132609619621.032571.04735.059.0FRFrance
181719852032789620885.034907.05138.064.0FRFrance
181819851934315432821.053487.07859.097.0FRFrance
181919851834055529935.051175.07455.093.0FRFrance
182019851733405324366.043740.06244.080.0FRFrance
182119851635036236451.064273.09166.0116.0FRFrance
182219851536388145538.082224.011683.0149.0FRFrance
18231985143134545114400.0154690.0244207.0281.0FRFrance
18241985133197206176080.0218332.0357319.0395.0FRFrance
18251985123245240223304.0267176.0445405.0485.0FRFrance
18261985113276205252399.0300011.0501458.0544.0FRFrance
18271985103353231326279.0380183.0640591.0689.0FRFrance
18281985093369895341109.0398681.0670618.0722.0FRFrance
18291985083389886359529.0420243.0707652.0762.0FRFrance
18301985073471852432599.0511105.0855784.0926.0FRFrance
18311985063565825518011.0613639.01026939.01113.0FRFrance
18321985053637302592795.0681809.011551074.01236.0FRFrance
18331985043424937390794.0459080.0770708.0832.0FRFrance
18341985033213901174689.0253113.0388317.0459.0FRFrance
183519850239758680949.0114223.0177147.0207.0FRFrance
183619850138548965918.0105060.0155120.0190.0FRFrance
183719845238483060602.0109058.0154110.0198.0FRFrance
1838198451310172680242.0123210.0185146.0224.0FRFrance
18391984503123680101401.0145959.0225184.0266.0FRFrance
1840198449310107381684.0120462.0184149.0219.0FRFrance
184119844837862060634.096606.0143110.0176.0FRFrance
184219844737202954274.089784.013199.0163.0FRFrance
184319844638733067686.0106974.0159123.0195.0FRFrance
18441984453135223101414.0169032.0246184.0308.0FRFrance
184519844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1846 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202011 3 101704 93652.0 109756.0 154 142.0 \n", + "1 202010 3 104977 96650.0 113304.0 159 146.0 \n", + "2 202009 3 110696 102066.0 119326.0 168 155.0 \n", + "3 202008 3 143753 133984.0 153522.0 218 203.0 \n", + "4 202007 3 183610 172812.0 194408.0 279 263.0 \n", + "5 202006 3 206669 195481.0 217857.0 314 297.0 \n", + "6 202005 3 187957 177445.0 198469.0 285 269.0 \n", + "7 202004 3 122331 113492.0 131170.0 186 173.0 \n", + "8 202003 3 78413 71330.0 85496.0 119 108.0 \n", + "9 202002 3 53614 47654.0 59574.0 81 72.0 \n", + "10 202001 3 36850 31608.0 42092.0 56 48.0 \n", + "11 201952 3 28135 23220.0 33050.0 43 36.0 \n", + "12 201951 3 29786 25042.0 34530.0 45 38.0 \n", + "13 201950 3 34223 29156.0 39290.0 52 44.0 \n", + "14 201949 3 25662 21414.0 29910.0 39 33.0 \n", + "15 201948 3 22367 18055.0 26679.0 34 27.0 \n", + "16 201947 3 18669 14759.0 22579.0 28 22.0 \n", + "17 201946 3 16030 12567.0 19493.0 24 19.0 \n", + "18 201945 3 10138 7160.0 13116.0 15 10.0 \n", + "19 201944 3 7822 5010.0 10634.0 12 8.0 \n", + "20 201943 3 9487 6448.0 12526.0 14 9.0 \n", + "21 201942 3 7747 5243.0 10251.0 12 8.0 \n", + "22 201941 3 7122 4720.0 9524.0 11 7.0 \n", + "23 201940 3 8505 5784.0 11226.0 13 9.0 \n", + "24 201939 3 7091 4462.0 9720.0 11 7.0 \n", + "25 201938 3 4897 2891.0 6903.0 7 4.0 \n", + "26 201937 3 3172 1367.0 4977.0 5 2.0 \n", + "27 201936 3 2295 728.0 3862.0 3 1.0 \n", + "28 201935 3 1010 2.0 2018.0 2 0.0 \n", + "29 201934 3 1672 279.0 3065.0 3 1.0 \n", + "... ... ... ... ... ... ... ... \n", + "1816 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1817 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1818 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1819 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1820 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1821 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1822 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1823 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1824 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1825 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1826 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1827 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1828 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1829 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1830 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1831 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1832 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1833 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1834 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1835 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1836 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1837 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1838 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1839 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1840 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1841 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1842 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1843 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1844 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1845 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 166.0 FR France \n", + "1 172.0 FR France \n", + "2 181.0 FR France \n", + "3 233.0 FR France \n", + "4 295.0 FR France \n", + "5 331.0 FR France \n", + "6 301.0 FR France \n", + "7 199.0 FR France \n", + "8 130.0 FR France \n", + "9 90.0 FR France \n", + "10 64.0 FR France \n", + "11 50.0 FR France \n", + "12 52.0 FR France \n", + "13 60.0 FR France \n", + "14 45.0 FR France \n", + "15 41.0 FR France \n", + "16 34.0 FR France \n", + "17 29.0 FR France \n", + "18 20.0 FR France \n", + "19 16.0 FR France \n", + "20 19.0 FR France \n", + "21 16.0 FR France \n", + "22 15.0 FR France \n", + "23 17.0 FR France \n", + "24 15.0 FR France \n", + "25 10.0 FR France \n", + "26 8.0 FR France \n", + "27 5.0 FR France \n", + "28 4.0 FR France \n", + "29 5.0 FR France \n", + "... ... ... ... \n", + "1816 59.0 FR France \n", + "1817 64.0 FR France \n", + "1818 97.0 FR France \n", + "1819 93.0 FR France \n", + "1820 80.0 FR France \n", + "1821 116.0 FR France \n", + "1822 149.0 FR France \n", + "1823 281.0 FR France \n", + "1824 395.0 FR France \n", + "1825 485.0 FR France \n", + "1826 544.0 FR France \n", + "1827 689.0 FR France \n", + "1828 722.0 FR France \n", + "1829 762.0 FR France \n", + "1830 926.0 FR France \n", + "1831 1113.0 FR France \n", + "1832 1236.0 FR France \n", + "1833 832.0 FR France \n", + "1834 459.0 FR France \n", + "1835 207.0 FR France \n", + "1836 190.0 FR France \n", + "1837 198.0 FR France \n", + "1838 224.0 FR France \n", + "1839 266.0 FR France \n", + "1840 219.0 FR France \n", + "1841 176.0 FR France \n", + "1842 163.0 FR France \n", + "1843 195.0 FR France \n", + "1844 308.0 FR France \n", + "1845 213.0 FR France \n", + "\n", + "[1846 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" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "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
160919891930NaNNaN0NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1609 198919 3 0 NaN NaN 0 NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1609 FR France " + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "raw_data[raw_data.isnull().any(axis=1)]\n" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": {}, + "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
0202011310170493652.0109756.0154142.0166.0FRFrance
1202010310497796650.0113304.0159146.0172.0FRFrance
22020093110696102066.0119326.0168155.0181.0FRFrance
32020083143753133984.0153522.0218203.0233.0FRFrance
42020073183610172812.0194408.0279263.0295.0FRFrance
52020063206669195481.0217857.0314297.0331.0FRFrance
62020053187957177445.0198469.0285269.0301.0FRFrance
72020043122331113492.0131170.0186173.0199.0FRFrance
820200337841371330.085496.0119108.0130.0FRFrance
920200235361447654.059574.08172.090.0FRFrance
1020200133685031608.042092.05648.064.0FRFrance
1120195232813523220.033050.04336.050.0FRFrance
1220195132978625042.034530.04538.052.0FRFrance
1320195033422329156.039290.05244.060.0FRFrance
1420194932566221414.029910.03933.045.0FRFrance
1520194832236718055.026679.03427.041.0FRFrance
1620194731866914759.022579.02822.034.0FRFrance
1720194631603012567.019493.02419.029.0FRFrance
182019453101387160.013116.01510.020.0FRFrance
19201944378225010.010634.0128.016.0FRFrance
20201943394876448.012526.0149.019.0FRFrance
21201942377475243.010251.0128.016.0FRFrance
22201941371224720.09524.0117.015.0FRFrance
23201940385055784.011226.0139.017.0FRFrance
24201939370914462.09720.0117.015.0FRFrance
25201938348972891.06903.074.010.0FRFrance
26201937331721367.04977.052.08.0FRFrance
2720193632295728.03862.031.05.0FRFrance
28201935310102.02018.020.04.0FRFrance
2920193431672279.03065.031.05.0FRFrance
.................................
181619852132609619621.032571.04735.059.0FRFrance
181719852032789620885.034907.05138.064.0FRFrance
181819851934315432821.053487.07859.097.0FRFrance
181919851834055529935.051175.07455.093.0FRFrance
182019851733405324366.043740.06244.080.0FRFrance
182119851635036236451.064273.09166.0116.0FRFrance
182219851536388145538.082224.011683.0149.0FRFrance
18231985143134545114400.0154690.0244207.0281.0FRFrance
18241985133197206176080.0218332.0357319.0395.0FRFrance
18251985123245240223304.0267176.0445405.0485.0FRFrance
18261985113276205252399.0300011.0501458.0544.0FRFrance
18271985103353231326279.0380183.0640591.0689.0FRFrance
18281985093369895341109.0398681.0670618.0722.0FRFrance
18291985083389886359529.0420243.0707652.0762.0FRFrance
18301985073471852432599.0511105.0855784.0926.0FRFrance
18311985063565825518011.0613639.01026939.01113.0FRFrance
18321985053637302592795.0681809.011551074.01236.0FRFrance
18331985043424937390794.0459080.0770708.0832.0FRFrance
18341985033213901174689.0253113.0388317.0459.0FRFrance
183519850239758680949.0114223.0177147.0207.0FRFrance
183619850138548965918.0105060.0155120.0190.0FRFrance
183719845238483060602.0109058.0154110.0198.0FRFrance
1838198451310172680242.0123210.0185146.0224.0FRFrance
18391984503123680101401.0145959.0225184.0266.0FRFrance
1840198449310107381684.0120462.0184149.0219.0FRFrance
184119844837862060634.096606.0143110.0176.0FRFrance
184219844737202954274.089784.013199.0163.0FRFrance
184319844638733067686.0106974.0159123.0195.0FRFrance
18441984453135223101414.0169032.0246184.0308.0FRFrance
184519844436842220056.0116788.012537.0213.0FRFrance
\n", + "

1845 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202011 3 101704 93652.0 109756.0 154 142.0 \n", + "1 202010 3 104977 96650.0 113304.0 159 146.0 \n", + "2 202009 3 110696 102066.0 119326.0 168 155.0 \n", + "3 202008 3 143753 133984.0 153522.0 218 203.0 \n", + "4 202007 3 183610 172812.0 194408.0 279 263.0 \n", + "5 202006 3 206669 195481.0 217857.0 314 297.0 \n", + "6 202005 3 187957 177445.0 198469.0 285 269.0 \n", + "7 202004 3 122331 113492.0 131170.0 186 173.0 \n", + "8 202003 3 78413 71330.0 85496.0 119 108.0 \n", + "9 202002 3 53614 47654.0 59574.0 81 72.0 \n", + "10 202001 3 36850 31608.0 42092.0 56 48.0 \n", + "11 201952 3 28135 23220.0 33050.0 43 36.0 \n", + "12 201951 3 29786 25042.0 34530.0 45 38.0 \n", + "13 201950 3 34223 29156.0 39290.0 52 44.0 \n", + "14 201949 3 25662 21414.0 29910.0 39 33.0 \n", + "15 201948 3 22367 18055.0 26679.0 34 27.0 \n", + "16 201947 3 18669 14759.0 22579.0 28 22.0 \n", + "17 201946 3 16030 12567.0 19493.0 24 19.0 \n", + "18 201945 3 10138 7160.0 13116.0 15 10.0 \n", + "19 201944 3 7822 5010.0 10634.0 12 8.0 \n", + "20 201943 3 9487 6448.0 12526.0 14 9.0 \n", + "21 201942 3 7747 5243.0 10251.0 12 8.0 \n", + "22 201941 3 7122 4720.0 9524.0 11 7.0 \n", + "23 201940 3 8505 5784.0 11226.0 13 9.0 \n", + "24 201939 3 7091 4462.0 9720.0 11 7.0 \n", + "25 201938 3 4897 2891.0 6903.0 7 4.0 \n", + "26 201937 3 3172 1367.0 4977.0 5 2.0 \n", + "27 201936 3 2295 728.0 3862.0 3 1.0 \n", + "28 201935 3 1010 2.0 2018.0 2 0.0 \n", + "29 201934 3 1672 279.0 3065.0 3 1.0 \n", + "... ... ... ... ... ... ... ... \n", + "1816 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "1817 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "1818 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "1819 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "1820 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "1821 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "1822 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "1823 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "1824 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "1825 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "1826 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "1827 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "1828 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "1829 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "1830 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "1831 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "1832 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "1833 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "1834 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "1835 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "1836 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "1837 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "1838 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "1839 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "1840 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "1841 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "1842 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "1843 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "1844 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "1845 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 166.0 FR France \n", + "1 172.0 FR France \n", + "2 181.0 FR France \n", + "3 233.0 FR France \n", + "4 295.0 FR France \n", + "5 331.0 FR France \n", + "6 301.0 FR France \n", + "7 199.0 FR France \n", + "8 130.0 FR France \n", + "9 90.0 FR France \n", + "10 64.0 FR France \n", + "11 50.0 FR France \n", + "12 52.0 FR France \n", + "13 60.0 FR France \n", + "14 45.0 FR France \n", + "15 41.0 FR France \n", + "16 34.0 FR France \n", + "17 29.0 FR France \n", + "18 20.0 FR France \n", + "19 16.0 FR France \n", + "20 19.0 FR France \n", + "21 16.0 FR France \n", + "22 15.0 FR France \n", + "23 17.0 FR France \n", + "24 15.0 FR France \n", + "25 10.0 FR France \n", + "26 8.0 FR France \n", + "27 5.0 FR France \n", + "28 4.0 FR France \n", + "29 5.0 FR France \n", + "... ... ... ... \n", + "1816 59.0 FR France \n", + "1817 64.0 FR France \n", + "1818 97.0 FR France \n", + "1819 93.0 FR France \n", + "1820 80.0 FR France \n", + "1821 116.0 FR France \n", + "1822 149.0 FR France \n", + "1823 281.0 FR France \n", + "1824 395.0 FR France \n", + "1825 485.0 FR France \n", + "1826 544.0 FR France \n", + "1827 689.0 FR France \n", + "1828 722.0 FR France \n", + "1829 762.0 FR France \n", + "1830 926.0 FR France \n", + "1831 1113.0 FR France \n", + "1832 1236.0 FR France \n", + "1833 832.0 FR France \n", + "1834 459.0 FR France \n", + "1835 207.0 FR France \n", + "1836 190.0 FR France \n", + "1837 198.0 FR France \n", + "1838 224.0 FR France \n", + "1839 266.0 FR France \n", + "1840 219.0 FR France \n", + "1841 176.0 FR France \n", + "1842 163.0 FR France \n", + "1843 195.0 FR France \n", + "1844 308.0 FR France \n", + "1845 213.0 FR France \n", + "\n", + "[1845 rows x 10 columns]" + ] + }, + "execution_count": 6, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + " data = raw_data.dropna().copy()\n", + "data" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [], + "source": [ + "def convert_week(year_and_week_int):\n", + " year_and_week_str = str(year_and_week_int)\n", + " year = int(year_and_week_str[:4])\n", + " week = int(year_and_week_str[4:])\n", + " w = isoweek.Week(year, week)\n", + " return pd.Period(w.day(0), 'W')\n", + "\n", + "data['period'] = [convert_week(yw) for yw in data['week']]" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [], + "source": [ + "sorted_data = data.set_index('period').sort_index()" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n" + ] + } + ], + "source": [ + "periods = sorted_data.index\n", + "for p1, p2 in zip(periods[:-1], periods[1:]):\n", + " delta = p2.to_timestamp() - p1.end_time\n", + " if delta > pd.Timedelta('1s'):\n", + " print(p1, p2)" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 16, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm4HUWZ9+89y91v7s1OyEIChCVsAjGACCj7jAvMiCM4Co7MoH586jjjAjMo84E44DqiI8ooq44IiIIiSwiEnZCENWQhCUnInpvc5OYuucs5p74/uqpP9znV3dWn+2y57+957nP6VtfydnV1vfUu9RYJIcBgMBgMRlQkqk0Ag8FgMPYPMENhMBgMRixghsJgMBiMWMAMhcFgMBixgBkKg8FgMGIBMxQGg8FgxAJmKAwGg8GIBcxQGAwGgxELmKEwGAwGIxakqk1AJTFhwgQxc+bMapPBYDAYdYWlS5fuFEJMDMo3qhjKzJkzsWTJkmqTwWAwGHUFItpgko9VXgwGg8GIBcxQGAwGgxELmKEwGAwGIxYwQ2EwGAxGLGCGwmAwGIxYEMhQiOg2ItpBRMscaeOIaD4RrZa/Yx33riaiNUS0iojOc6SfSERvyns3ExHJ9EYi+p1MX0REMx1lLpNtrCaiyxzps2Te1bJsQ/SuYDAYDEYUmEgodwA4vyDtKgALhBCzASyQ/4OI5gC4GMBRsszPiCgpy9wC4AoAs+WfqvNyALuFEIcC+BGAm2Rd4wBcC+AkAPMAXOtgXDcB+JFsf7esg8FgMBhVRCBDEUI8A6C7IPkCAHfK6zsBXOhIv0cIMSSEWAdgDYB5RDQFwBghxIvCOnP4roIyqq77AZwlpZfzAMwXQnQLIXYDmA/gfHnvTJm3sP39HtmcwL2LNyKTzVWbFAaDwXChVBvKZCHEVgCQv5Nk+lQAGx35Nsm0qfK6MN1VRgiRAdADYLxPXeMB7JF5C+va7/G/L7+Lr//+Ddz1otE+IwaDwagY4jbKkyZN+KSXUsavrmKCiK4goiVEtKSrq8srW92gq3cIANCzb6TKlDAYDIYbpTKU7VKNBfm7Q6ZvAjDdkW8agC0yfZom3VWGiFIAOmCp2Lzq2gmgU+YtrKsIQohbhRBzhRBzJ04MDEVT88jmLFVXOqnjqwwGg1E9lMpQHgKgvK4uA/CgI/1i6bk1C5bx/WWpFuslopOlDeTSgjKqrosAPCntLI8BOJeIxkpj/LkAHpP3npJ5C9vf75HJWcJYMsEe3wwGo7YQGBySiH4L4AMAJhDRJlieVzcCuJeILgfwLoCPA4AQ4i0iuhfAcgAZAFcKIbKyqi/A8hhrBvCI/AOAXwG4m4jWwJJMLpZ1dRPR9QAWy3zXCSGUc8A3ANxDRN8G8KqsY1Qgk7UYSirBEgqDwagtBDIUIcQlHrfO8sh/A4AbNOlLABytSR+EZEiae7cBuE2T/g4sV+JRh6yUUFKs8mIwGDUG1pvUGTLShsISCoPBqDUwQ6kzqO0nMtAAg8Fg1AyYoTAYDAYjFjBDqTMoTZfnxhsGg8GoEpih1BlsTZdglsJgMGoLzFDqDCQDBTA7YTAYtQZmKHUGJaHkcsxSGAxGbYEZSp2BfbsYDEatghlKnUG5C7N8wmAwag3MUOoUbJNnMBi1BmYodQZit2EGg1GjYIZSZ7C9vFhEYTAYNQZmKHUGjrjCYDBqFcxQ6gy8r5HBYNQqmKHUGex9KHXAUZ5Yvh0zr3oY23oGq00Kg8GoAJih1BnqyW34ty+/CwB4c3NPlSlhMBiVADOUOkM9qbzY3sNgjC4wQ6k32G7DdcBRGAzGqAIzlDpD3m24yoSEALs4MxijA8xQ6gz1pUaqK2IZDEZEMEOpM/AUzWAwahXMUBhlByu8GIzRAWYojLKhvtRzDAYjKpihMMoGtsUzGKMLzFAYDAaDEQuYodQp6sEVl1VeDMboAjMURtlRB7yPwWDEAGYodQqqg+V/7VPIYDDiBDMUBoPBYMQCZij7Ie5dvBHzbniihuwstUIHg8EoJyIxFCL6ChG9RUTLiOi3RNREROOIaD4RrZa/Yx35ryaiNUS0iojOc6SfSERvyns3k9TnEFEjEf1Opi8iopmOMpfJNlYT0WVRnqMe4ccsvv77N7Cjd6jqtos60MoxGIwYUTJDIaKpAL4EYK4Q4mgASQAXA7gKwAIhxGwAC+T/IKI58v5RAM4H8DMiSsrqbgFwBYDZ8u98mX45gN1CiEMB/AjATbKucQCuBXASgHkArnUyLoYFlgsYDEYlEVXllQLQTEQpAC0AtgC4AMCd8v6dAC6U1xcAuEcIMSSEWAdgDYB5RDQFwBghxIvCWnbfVVBG1XU/gLOk9HIegPlCiG4hxG4A85FnQqMCfkZ5datWVF41QgaDwSgzSmYoQojNAL4P4F0AWwH0CCEeBzBZCLFV5tkKYJIsMhXARkcVm2TaVHldmO4qI4TIAOgBMN6nrlEDP2ahWE2u2iov9vNiMEYVoqi8xsKSIGYBOBBAKxF9yq+IJk34pJdappDOK4hoCREt6erq8iFv/wMfwsVgMCqJKCqvswGsE0J0CSFGADwA4H0Atks1FuTvDpl/E4DpjvLTYKnINsnrwnRXGalW6wDQ7VNXEYQQtwoh5goh5k6cOLHER61P1IqqqUbIYDAYZUYUhvIugJOJqEXaNc4CsALAQwCU19VlAB6U1w8BuFh6bs2CZXx/WarFeonoZFnPpQVlVF0XAXhS2lkeA3AuEY2VktK5Mm2/Rz15TtUTrQwGIzpSpRYUQiwiovsBvAIgA+BVALcCaANwLxFdDovpfFzmf4uI7gWwXOa/UgiRldV9AcAdAJoBPCL/AOBXAO4mojWwJJOLZV3dRHQ9gMUy33VCiO5Sn6WeYCJ1EBEgBHK1IqIwGIxRgZIZCgAIIa6F5b7rxBAsaUWX/wYAN2jSlwA4WpM+CMmQNPduA3BbSJJHFWqFn9QKHQwGo7zgnfJ1hjBqpGrP46zyYjBGF5ih7MeolX0oDAZjdIAZyn6IWtmHosDuywzG6AAzlP0ZvLGRwWBUEMxQ6hR+2iw79Eq1OYpEVM3bt/+8HPcu2RickcFgVBWRvLwY1YPfHG1JBqL6Kq+YBJRfPrcOAPB3c6cH5GQwGNUESyj7Mdgoz2AwKglmKHUKE15RK+ykVuhgMBjlBTOUOoWJfaTaAgqb5BmM0QVmKPsjauw8FAaDMTrADKVOUVcqL2ZsDMaoADOUOoXJFF3tedzvVEkGg7H/gRlKnSEMk+BowwwGo5JghlKvMDgCmNkJg8GoJJih1BlEwa9v3ipLKKzwYjBGF5ih1BnC8Ihqa7xYQmIwRheYodQpjGJ58YzOYDAqCGYodQaTDY0qym+1g0OyyovBGF1ghlKnqIed8gwGY3SBGUq9QTKJ+trYWG0KGAxGJcAMpc5gMjcrG0q196HwvkYGY3SBGUqdoh52yjMYjNEFZih1BrW3xIxZ1AZHqbZzAIPBqAyYoeyHUJqmap/YyBovBmN0gRlKnUFJJvXk5VUrdDAYjPKCGUqdQRRd+OWttlG+cjLK5j37KtYWg8HQgxlKCKzZ0Ys/vb6l2mQYY7RIBg++thmn3vgkXli7s9qkMBijGsxQQuDsHz6DL/721arSkFd5eUNJBtV2G1YoNxmvbNgNAFi1rbe8DTEYDF8wQ6lTmEQSrjY/qZTCSz0mOwEwGNUFM5Q6g1ksr9EJPiGSwagumKHUGUSY0Cu1ofEqu2tArTwngzHaEYmhEFEnEd1PRCuJaAURnUJE44hoPhGtlr9jHfmvJqI1RLSKiM5zpJ9IRG/KezeTXGoSUSMR/U6mLyKimY4yl8k2VhPRZVGeY79DjYReqZSopKQ2FlAYjOoiqoTyYwCPCiGOAHAcgBUArgKwQAgxG8AC+T+IaA6AiwEcBeB8AD8joqSs5xYAVwCYLf/Ol+mXA9gthDgUwI8A3CTrGgfgWgAnAZgH4Fon4yo3qn0SImAYeqXsVNQG1OtgfsJgVBclMxQiGgPgdAC/AgAhxLAQYg+ACwDcKbPdCeBCeX0BgHuEEENCiHUA1gCYR0RTAIwRQrworJn6roIyqq77AZwlpZfzAMwXQnQLIXYDmI88Eyo7stXegg5TlVf16QTKT4ddO4soDEZVEUVCORhAF4DbiehVIvolEbUCmCyE2AoA8neSzD8VwEZH+U0ybaq8Lkx3lRFCZAD0ABjvU1dFkK3iRG0yOddO6JXKTvDMThiM6iIKQ0kBOAHALUKI4wH0Q6q3PKD73oVPeqll3I0SXUFES4hoSVdXlw955qimhCLs3/pRekWhoh7coxkMhoUoDGUTgE1CiEXy//thMZjtUo0F+bvDkX+6o/w0AFtk+jRNuqsMEaUAdADo9qmrCEKIW4UQc4UQcydOnFjCYxajflRe5afDCBHoCBNVmTVeDEZ1UTJDEUJsA7CRiA6XSWcBWA7gIQDK6+oyAA/K64cAXCw9t2bBMr6/LNVivUR0srSPXFpQRtV1EYAnpZ3lMQDnEtFYaYw/V6ZVBLlcpVoqRhgmUW1+oib4KDHFwpz7UmkVG4PBcCMVsfwXAfyGiBoAvAPgH2AxqXuJ6HIA7wL4OAAIId4ionthMZ0MgCuFEFlZzxcA3AGgGcAj8g+wDP53E9EaWJLJxbKubiK6HsBime86IUR3xGcxRjVtKCawQ6/UgCQFRJOUTFyfbYbC/ITBqCoiMRQhxGsA5mpuneWR/wYAN2jSlwA4WpM+CMmQNPduA3BbGHrjQnVtKOqALYOJttzEGCKaDcWkfqnyitAOg8GIDt4pXwKqylDCqLyq7uVlIQodYdRlLKG4sbF7ADOvehhPvx2PMwqDEQRmKCWgFlRe/tGGVZ7q0wlEtKGEcD5gG4obr7xrRWG+b8nGgJwMRjxghlICqmmbsN2GfUiIQzKIAzZjKzMd9usYJfxkYDiDp1btCMzHwTIZlQYzlBBIyO+zFtyGTVBthqIQhQwjo/wos6Fc/cCb+IfbF2NtV1+1SWEwXGCGEgIJueKr7k55+WtypnyNqLyicDajoraX1+hgKe909QMA+gYzRvmrPQqEEHhi+XZkslX0t2dUBMxQQiAhRZTqSigGoVfsExvLTUsAHVJmiOTlFaq98kIIgW/+cRmWb9lb5pb8kbeRBeQrOyVmeHLlDvzjXUtwy8K11SaFUWYwQwmBWlJ51UNwyLyLc4Q6asg9uqtvCHe/tAGX3vZyhVrUI28jM3zyKg/XLXv2AQC27R2sLiGMsoMZSgjYKq8acBuun0he0VBLz9CYtE5bGBrJBuQsM8hM8qsVDeBQxlJ1NaT2j+lm5ba9mHnVw3jpnV3VJqXmsH+84QohYauSamma80GtqLyi2FAM1O6q/jgm0L6hDHb1DenbkR06mKkuQ1GPabqwqbYtTTGUxlQyIGd9YPE6KyjHn9/Qhg8c1WCGEgJUAyovkyOA8+Hra4PxRbOhVHZj41k/WIgTv/2E9p567SPZ6vares6RjD+3rZV9OfubhNLSYAUY6R+qsqRag9g/3nCFkKwBo3x+gq2fsO7RbCjx0WGC7Xv10glQfZuUgiIjYyqhVJlstW8rlagNBhcVrY2WpNU/ZOZlN5rADCUEasGGEgbVptLUG8kPYcqWe+Ksdn8qKMYW5L5eKzaU/Q3JhDVt1ooGoJbADCUEbC+vWtiH4qfysneo18aAj2RDqZFnAOKbQPYNZ/Hc6p0ll1du4bXUN36oVMQERvXBDCUE8mHhq0wIgj7O2tiHEgdC7UMp+0aUeKq5+oE38KlfLcL6nf0llc/b8szyR53In1i+HV293qrAIOyvghIzyGIwQwmBmpBQypg7bsSxMq0ltUJcDHqt3Om+d3CkpPKmqld7v0qEcTCUyeIf71qCT/1yUXDmAFTb2ywu7K8MMg4wQwmBpP0hey8NX1y7C394dVPZaAgVeqVGvt9IE0mNPAMQ34QYNeKCWtgEMds4GLoqu25XadKUi5AawK9f2oC/+vGz1SZjv0XUExtHFchmKN55LvmflwAAf3P8tLLSYmJDqRWVV7TzUGoHcTFoU4bgBbJVmhUMvhLDs9fCAueaPy6Lra4aeJyaA0soIVAT+1BCDOPqqxhiiOUVonD5w+TH00DSYGHih7DjMI7+jzKW8qq32sFIhECVNSRw1RyYoYRALe2UNwq9Un0yAVTuxMayuw3HJaEkoo0j03EYx8QXRzy2IDp27B2MZPQvBQPDpW9KrJXvqhbBDCUElKrCdENZWWDQdK3slI/j5MgwXV0v37mt8irVhiK/2kp4eYWJHVcqIfO+swDvvUEfnSBuqN36ewaGS67D3lrMnKUIzFBCwF5Z1viJjbWGSkUbLvcHHpvKSxnly2xDiUMzo9qI0rdxHGMQF2aObwEAvFOiyzaQ74taeJ5aAzOUEKilnfJ+q36KaPStJYSyoZSPDKv+2Izy0cYRRZRwwsA80I83amljY1PaCpsShUHWwOdfs2CGEgI1sQ8lRNvV/oBDn9sRFRGbCaIzbgmlZC+v0CeHRohUINVqkWwopReNHbY6OMLmZFtCYcZSBGYoIWAbQw2WKGWfRP3chmtsp3y5g0PGZTMK6q+4ujNhh04ptbz1GzQO46A3Tk/B6nsd5hFlUVg7T1F7YIYSAmplaGKUL9dkbqKCqDWVVyS3VZMNnDG0AwT3V1zdGXWxkWegZu1Ei1RQelmFWlJ5IYY4aLXyXdUimKGEQFL2lsmAKtegC7Nir/ZSKp7QK+Z5o3Z5MEOpjYmEDN2G4yA3jmemGtq4YcqM/VArkn8tghlKCIQxppZ7FeP3oZtOOJVCFFVHKJtRRA4a1FRcvRl1gjWVPNQwjUJ3nJNnbYxGC1G+Dfby8gYzlBCgMAylTBGJw6h3amXAVyr0SrkllNph0Kb5otNbS3aPOBBHWKIaGQY1CWYoIRAmBlMtTD61QANQydArZTbKx9ydUesLKm5LKJFW4yUXLWtdUcE2lPKAGUoIJEIZ5ctlQzE3tFZ73NvtRyIkjMorGipllI/qlWZaLA77RxzPHEfEhLgQh0dgHIx6f0VkhkJESSJ6lYj+LP8fR0TziWi1/B3ryHs1Ea0holVEdJ4j/UQielPeu5mkbomIGonodzJ9ERHNdJS5TLaxmogui/ocJkiGcBuuBS+vWhnwldLhR17xB9ok4u3PqGMk6P3GQW4cz6zc2HUDoVpjNI59KIxixCGhfBnACsf/VwFYIISYDWCB/B9ENAfAxQCOAnA+gJ8RUVKWuQXAFQBmy7/zZfrlAHYLIQ4F8CMAN8m6xgG4FsBJAOYBuNbJuMqFMCfl1cKgqzYFcQgolVR5VfqdldqeaSk7bEpJrYRryw/kzU+qJkVHM8rHSMh+hkgMhYimAfgQgF86ki8AcKe8vhPAhY70e4QQQ0KIdQDWAJhHRFMAjBFCvCisL+yugjKqrvsBnCWll/MAzBdCdAshdgOYjzwTKhtshmIwosoWnsVA3K5kaA4/xBH6PMxBYtFVXmbtREWlzquJRUKJgUg/n7ZK2yMo4qZSgG0ofogqofwXgK8DcK7ZJwshtgKA/J0k06cC2OjIt0mmTZXXhemuMkKIDIAeAON96ior1DiqpsorDP7jT8sxlCk9THdcqJyEUno7QOW9vEq3oZi6DZvb2yoB3SKoWqTxTvnyoGSGQkQfBrBDCLHUtIgmTfikl1rG3SjRFUS0hIiWdHV1GRHqBdWAiVG+XOoT+3wKnzzk6J6nVkZ75miIQeViUDiuyACBRvlItYdvLyriqD1Oo7wOFZdQYmi3liSUhat2YOGqHdUmw0YUCeVUAB8lovUA7gFwJhH9GsB2qcaC/FVPuwnAdEf5aQC2yPRpmnRXGSJKAegA0O1TVxGEELcKIeYKIeZOnDixtCe1K7N+zNyGozXlSYKJVb5GYKuiokgoFXzQShvloyKob+LYgBfnM+uqqnSXxqFurAXtg8Jnbl+Mz9y+uNpk2CiZoQghrhZCTBNCzIRlbH9SCPEpAA8BUF5XlwF4UF4/BOBi6bk1C5bx/WWpFuslopOlfeTSgjKqrotkGwLAYwDOJaKx0hh/rkwrK9QHbGIfqWZE4hqKdAEg6k558zzlD70Srf6w7UWlI5bQK9Gr8K2rWkw62t6c2lIl1hLKsQ/lRgDnENFqAOfI/yGEeAvAvQCWA3gUwJVCCKXg/wIsw/4aAGsBPCLTfwVgPBGtAfAvkB5jQohuANcDWCz/rpNpZYUaQGY75cuk8jIwdDv5SanMZXAkiz+8uimeTXGVsqFEnP6cr0yr74/bhhIxmoJx6JUyq3cy2Ry29Qx63verouISCsxd/73AjMQbqTgqEUIsBLBQXu8CcJZHvhsA3KBJXwLgaE36IICPe9R1G4DbSqW5FKhxZPKRlXvQ+dUfRzC+Gx9ZiTteWI8JbY04bXZEVWEEVPJMeeckkxNAsqAb41sjRIu1ZtonsYReMajiuj8vx10vbsDr3zoXHS1pDR3edVVLQomm8mKO4gXeKR8CaqVX1Z3yBpNEHBqvnX1DAIDu/ihnb1faKB+hoYK2dKv6uN9pqdWZuknHwQBNpJvH39oOAOgfznjU4V224vaIGBw48kE3mbEUghlKCKjhMzQSrKsolw0lrL2gVObSIGP1j2SjqwYiqc3CtBNZ5eWWUPzux4Go9QUVjyX0ikEe9RzqJMriOoTr13WvajaUKGWZkXiBGUoIqHHUP6RfibnzlnfQ+U6eMYgoDSlraAxnoodNrtQmsjiN8vrJL1r9ClElqkoa5cMEQvXStPotgnin/P4FZighoMaRl2hfqIMvJw3lRtqWUEpnKDHY5Cv68bqN8rr7NSahBPRsHBsbTcrmfBhGcNnKzs7xHLDFXl5eYIYSBnIE7RvW7z53qrnKf8CW9z23l1dp4kocEko87rxhJJSoKiSHhKJlKJGq920vVDlTo3yFJZRSTrys1p6OSBKK+mWGUgRmKCGgxo+XUd7pTly2A7YMDbJRofThcdiCyr0PpZS8OrgkFA3NtRdtOKh+b9tFXG0Aecncy53ej3FW2h4RR5w79vLyBjOUEAjah+JiKGUbdJUZzEquiWOXe7Q6zBF1gs4FSChxT34lSyghi1Vb5eW3CKq0hGLvQ4mhT9jLqxjMUEJADSBPCaVGVF5xtGxHZY20KzE6HWFWkvF6eWkkFIfUGQdzKdkob1x/HAzdXOXlKaFErL8ciOQ2XEuxV2oMzFBCIC+h6PVZ2axzQiovDf6ZoreTP6Qrel1ldxu22ym5maLy+tV0PO83jgCFQHC/xmNDCc6jGInX8/jZ0io9N+el5vKOydEKZighoMZgLUgopsO6VA/ihM1QYjBellyD2aQYRxBEwL3CDpr84ni/JVdhWC6ODXgm719l8WQoPvtQKr3azwn3b2l1MEvxAjOUEFDDyMiGEvFDGRjOoGffiCcN5VZ5JSgOXXNlVC52DoOG9gwM47N3LLYjATjhtqEU1yViWjDEFW4/cGNjHKFXDPKohZTXWPEdq06psBITdQDzM0EUN+n9HcxQQsAOveKxe9xtlI/W1ge+txDH/b/HPWnwQxwfZhxGeYWy22FCrDp//dIGPLlyB25/fl3RvaB9KFkXwzGgKwCl21DMClZKXWlqQ9E6OiDePg1CLoD5GYH3oXiCGUoJqISX147e4hW0E/6GzjxKjhMZg1E+DtuGGT/xVqkUYlCGzWlMJYvrcTIMTdm4VV7RNzYG1J+LQUI0UjnK9gLcvPThbJzX3o0t3dCNE66fjzU7+oIJ8kHeUSG6hMKqr2IwQwmBQBtKBdyG85O0d/1xNB3HjuI49swYbazLudvzgzoSuTFVPPSDJje3yiu4LS8o19Vyuw3HMQLDPKfX3is/Nu3sZ789T4+/tR3d/cORTydULcRxYmNQHTc9uhKfuf3lktupRzBDCYH8AVseXl4GKhEhBF5Yu7Os+mLnSr1UCSURo5tX2VfIIaIaD2WUhKJjKP7vL8itOCyiug0HkZDf2Fg6wozToD7RfTZB0Qnse/LX5CwiXxpCqEeDaAmq45aFa7FwVWlHcL/T1Ye/+vGz2DNQerTvaoAZSgiEkVC8Bv4jy7bhk/+zCP/78ruRaPBVecXAqxIxhIT3VyCFrcMnj3D/+kHFJksHMRStR5KjzRgiIZQ7lpddfYRmwrx/LwnDTyXmTPLrDxW5weToCF8E2HtMEIfaLAg/fWoNVmzdiydW1M558SZghhICQaskE5WX8tx6Y2NPNFr8VnNxqLyUgBIwG/XsG8HTb/uvwqJJKMGF8/NmcF71jpIa0c3tcVRcNi4JJWoEAdOJLO9OHZ3WMO151aG761Yzetet3lcQI1izoxfPr9npeT+OUyzzDLLkKgKR97KsLzsNM5QQCDpgK2gfAwC0NlqHZHpFLA6kIWR+KnEnir1TPqDBP7yyCZ+5/WX0aUL6x+M2bJAnRDsqeLLu7I7g0Cv6vGEhYpjUCunR3pe/cdjBTOAVmNrPiO2yofgQaiqhnP3DZ/D3v1zkeV8xtzjc4cs52cexD6waYIYSAnFIKIV1BbZZUI/JJj63t1JpA9L0zI59IzkIYZ1BX0RHRBoAQwklxASt7F96huKoM8CA7Ncv33tsJe56cb3n/ah6fNNipsZj37ZCFPVqx96TFSD1+b2/lHxfd2jcvcNAqS3j2IdSTgkl77hRvjbKgVjOlB81kC83mxMQQhSFhs8arrbCQHe2ORDg5RVDu6aeSGqC9tqbY9VROh1mRnlFS1BdAs+t2QUgWEIJcnH165f/fmotAODSU2Z6UVLUXhiYes/FoZoJQ6PXZl4lVQTbULzrVscp7B0sTbLPtxFdQqmEDSUhl/p+LdRiTDGWUELA+fp0DMNkY2NYBVQpjMk5zktVedkid0A+dUSw9iAuw4nPD0YMxTa0+nOUx97art0hX1hP4bVCXLG8IksohhNZLPaCEHm9nkeNYd1tJ2kmKq8pHU0hKPJGPDYUszp0z7Wrbwh3vrDeh45gG0ot2ldYQgkB58vPClHUeW4biv5l20KN4Vio1qAxPTdCPfOwhqFUKny9IjHo7JZtPfscZTQMwxVNWENLTEb5qCtcu1SQWjUGlVeYsl79nw8e6V+/X3+oOhIl79Q+HbYJAAAgAElEQVR1txdNDRhOyskJgWTBwu6ff/canl29E6ceOgGHTmorKmOics4W9F2ph+nFCZZQQiBIQjE5AjjsOC5sx0wF5BJRIiGoOaXO8DsquOw2FPkbJM0lk/nhrhNmgveh6POGRdgVrlf5QJWX/PUT3Dbs6vf1igrz6ryex0/lZSr1xSFtWeWD2wpC2J3yunxvb+8FAPRrnFmAvIbA74MPWgBVAyyhhIDzpem8TZxpXqu1sJNIYT0mq/44Blfe0O2fL5P1tqHYZcssoeRVXv65Uw67iX5yc7brP/lF6eOoenxTBp3zmcgVPvj9hcgJYP2NH9LXEYMNxd8ob9ZWLDG4XPXEIWEa5tcw9F5pC+r1sAkl7NBH3nDZbIVAIurqMQawhBICzg85q5lAXeHrAzzBwk4KfrQU34sOe3Ub8NUoJqpVeRmupBU+/JNncU/Bhs8wK9Igl1KnIT5IpVUXEkpAcRPJLWjlH87LS5/uJ6E4R4cvnSG/Gy+Y9p1vHfLX1L6pe261AMt4iI926COfNipzQmw4MEMJAec7G9EMBCeTMdmrYoKSVF4xaLxMx6eibyTjo/IyrGzZ5r246oE3Q9NhrzqDGIpDx6yTICsVyyv6xkazfOoZ/NSRCrp9REC4xYm3O73VfhCT9nsu9b5Mn917k6WiKcqCIJyUo8unnseLDjI4PsKlYo8hckMcYIYSAi4XR80LDCOhmKJI5WW0woq+WslvAAuSUKyOGNGpvELE2PKmwyCPULQEqLySQSovh4SiKW+6CS8I+WCW3nXsHRzBTxas1rYj7F9/GlT1undTiCGPBUGYla/X86jm9Y4QZqvsvPrOjBav9xOLyitEMFJnficUfUFj1u9u3LHl4gAzlBKhE1WdH4enx0tE/auJ/rayNhS5Ctb0Rxj1QhSVi6kKKemyoRTfD7KRZF1G0AgTkoE94LuPrsQP5r+NR5ZtLboXNvSKTh1ZCK99RGEe09vLS0oomnvmNhR3XUHwmqjjUXlFl1AUvBifbUPxK1tgQ6kFMEMJASGEPSnpBoLLKB9koDREsVE+GM48pboSmsaByviovHxssUXwPpwpuLTK47e5EnAb5fUnMsL3flz7UFRRv4lmJGPd27tPF9LG/RvUjonKyytPuGjD+nT1XrR9DmefBk+epsEhyyqhGCxgCrcYeMGboah6vOlwjVdWedUfBPKTkm5gm6hEwkZLLWRAJhN9HDt486t+/3x5t2FdRvOP16tfTLrL1I3TuYdBx9hd+4g05ePah2Kyj6G5wToAbJ82pI0aA/6wbSg+9i0FT4YSWNLRXuAkrqnfJaEE122qagySUOI5D8Ws/VIkFJNjol3RzVlCqT8IAaTlXgathOKYVIPifYU1ejtpCCpvOrR29Q3hwdc2+9YRRKcdesVH5WUilXl9OEbMUWYJYtbO2zr+V6l9KCZ7KtLS3pPx8Z4zbcfEhuLN0P3LmjDZjD3m/Zm4X1uqDtMFma7fgDzjjCRhGrw/J4P260Kv5yEDt+H9yoZCRNOJ6CkiWkFEbxHRl2X6OCKaT0Sr5e9YR5mriWgNEa0iovMc6ScS0Zvy3s0ke5OIGonodzJ9ERHNdJS5TLaxmoguK/U5wkAgr/LSqVeyBi84+j6UYJg28YVfv4Iv3/MatvUMetYRNKGryWpYq/Iyl1Aied3AbPUaFFrFvVr2l0Aj6eAN+sVPVWnK7NX94Wwu8D3q3p9JGy4mHbCI0t12O7J4t2PXEVHl5cfcTGEiYSqVpR8t1j0Pt2EDCcVZtFbiekWRUDIA/lUIcSSAkwFcSURzAFwFYIEQYjaABfJ/yHsXAzgKwPkAfkZE6mDvWwBcAWC2/Dtfpl8OYLcQ4lAAPwJwk6xrHIBrAZwEYB6Aa52Mq1xwSihBHivlUnmZnMJn+rFs3mOFItGpO0wNj7bbsNbLC7KOYFo8bSgmAorwr0MhSMIIWvHFLaGY9IsuS1ijPOAtpShdfakSigmT9YvlZerl5fSKMnl+r+cxkVC+et/rOOk7T3jeN1GxbgkI86PgZd6yz0PxY0YGNr1sTuCnT67G3sERz3riRMkMRQixVQjxirzuBbACwFQAFwC4U2a7E8CF8voCAPcIIYaEEOsArAEwj4imABgjhHhRWCPlroIyqq77AZwlpZfzAMwXQnQLIXYDmI88EyobcgJoSHrbUJxpnh9oyA1aXm7DpnotP5O8+jB1i2FTw69f6BVTKQcINqL61aPyBDGUoA8waE9EbPtQQkhuejrkr6HbMOBtI1ETl7dR3p8+l92wBJWXKZN22yT8aQK8x0KeoXhXcv/STdi+1zuIqMm+p+/8ZUU+v6/k5b+x0c+fwqTvn1y5A99//G1c/6fl3hXFiFhsKFIVdTyARQAmCyG2AhbTATBJZpsKYKOj2CaZNlVeF6a7ygghMgB6AIz3qUtH2xVEtISIlnR1lXa+s4IQwj46VjcQXCtcT/Ff1WXWZpENRdNWEZ1mVdv5dOoVO5xJAKFKV62bkEwner92TI5VFgH37XyBEoj/gsBEZ22ierDXBD5ZfRcCgS1IWlwSigdDSfgzlOA+1bfnosNWeRXfN4nQ7awD8N5d7oTXgk6pquPw0vN7f0MjeRpNGaWuDb9ndW9s9Pp+rPLqpNhyIzJDIaI2AL8H8M9CiL1+WTVpwie91DLuRCFuFULMFULMnThxog95wcgJYau8dDYUl1Hec4L09snXtlkwnowmesPK/SQIk48GCPLyshBlRWniqmurvEKoZ3QfoLMfdEZdtwTj8X4NVgpGtiWfI5jttgOact722ouiogd4uVxH7VPAKaH4l/djXiYu+U7oFnxCiHwYmCg2O4P3N31ci31dmE8YPLMq47Xh1MpTnL8QyYS3E1E5EImhEFEaFjP5jRDiAZm8XaqxIH93yPRNAKY7ik8DsEWmT9Oku8oQUQpAB4Bun7rKipzIuw1row2bSCgG4rIuv4L616+8847fNhQ/Y7axGslA5RXFKO+sNsgTLJhWZ13F951pwRKKVxsGEopw/+rgd2KfKPj1gpNeL6N7MkBCCVJXuic1fZ68Ud5fQvEPX5+nz8QOGaSSNhmTXn2SP/XRu+ykMY2ebZk4Mqjve8DnqHAT6S7lo6IvB6J4eRGAXwFYIYT4oePWQwAuk9eXAXjQkX6x9NyaBcv4/rJUi/US0cmyzksLyqi6LgLwpLSzPAbgXCIaK43x58q0skI4JZQSbShqjJa6Qcue6P1UXsZGW30bVjve95wwUXlFYigGaiaVGmyU968raLXsUu+UaMS28pjnDaLDF458XgxFLTi8JMygPZEmKkk/V3nTzaIuZmDCUDTP4xyjURiKicNKr8MIXkhuxoA5qqYHhov3IeXrDe77/L65yux8jBK+/lQAnwbwJhG9JtP+DcCNAO4lossBvAvg4wAghHiLiO4FsByWh9iVQgjVW18AcAeAZgCPyD/AYlh3E9EaWJLJxbKubiK6HsBime86IUR3hGcxghD5/QFBJzYGha8vNVJpngn40GlUs789xpROP5VXXpoKpsUkskCQeiCcDUVXT/5a96GbrAjNJBTvCbYwj58HXrD0EMyMg4zycexDMWco3m0F2be82nRin2NyNvn8vJisKuvXNX2OsPSFtJjsV1PPq4uUUJjHoqV0WuNEyQxFCPEcvG2HZ3mUuQHADZr0JQCO1qQPQjIkzb3bANxmSm8ccNpQ/A7YakgmvPXJ2XAMxdMo76fyMhw8Kp/uA83f8+cGviqvEDvlgyYjK49H2VxxXn1dDtp0dqMAG4qJesaEeZqoAk0YdRDcqhV9nmTACjbYbTg4b8a2Gxbf37InvwfKb0ybTMLuNjUMxRF1wMj1OEAN6Bs6x0e9ZuQNKsus7erzbMNkAaueoVIMhXfKh0BOAA2pYJVXOkmeH3BoCaWQoZiovAxlFD/bg2r3iRU7jHYE61Ve6re0CaCwbJCaKYzKK8gGVrKEEpNR3t534euO7d+O83ZQEEJvlVcYKcgrj/f97z22KrB8IR21IqH4Lgh8dso773m5DSv6vU50dNIBeC9k1DNEPUfGFMxQQiBQQhECREAqmQgcKKY6zZKM8kJ/XZSvgCYnnElLN+z2ps9HQoHPROJVj1964CbLUOqZ4vtBxlITe47JQsEkFpTqT730aLgYMaBXnYrsrfLyb9dEj6/GephNksV1ONoxCiVT/DxOCcXkPQU5KpgywMK23Pf05VVfmNpQAqVDllBqD04bim7AZnICqQQhmaDATV4G34SVz8Mo77dKE8hLUn7N+EooDvoHR7yZn1oB+Z2HYqTy8mjCJES33ScBnRqksnLe100mzonMJLRO0CFPfoxBPYv2JMyCX08EMEggL6F49V3O1We6+47mvJisHW3Yj9gwNpTgxZjueZULbipBkXbbBzHZwrJ+Ki+vhadK7h/O+KhXDRiKLaFUBsxQQiAnBFIBNpQEERLko/JSDMVQQikyystf3xWWAN4zvTOwblWDbjA6B7FiTjr4SSh+qo5CeOrwAyY0dztBK2DntfeEAwRLKF5NmajFhEG/2MEQI6i8TAzZYYzyuhqc781rzOfPAIogoRgwcxHwvGqzYXM6aaSa9LKhhHUfL2IoJie7Ohw3vPaiuBZbHvXY75UllNpDTgikfcLXZ20JxVslZZ/rYCiiFI5pE4OggHCcp+CTT97S0eIk33nSYSH8Q69IWkOogQrhrNY73lew1FbYhi7rUCavXgi2oejbcrrnBkkx/ioeqfLykfyCYEJvwlZ5+Y9XrzpcAQoD1IBBVJsEhwTMJAedWmw4a73f5oZkoDs0YGZXCpIMCukCCplwsPThZUdx9702S8X2nygwQwkBAUdwSA+jfCJBSCUSgWfKm6iBnPltGoQ+vRD2xjifPH4GfvMdzNaoHs4U5/GTgAoRtLr1q8fEruQsn0yQ9pmczEAbCcFA+vjoT5+zr4MCXppIKCNaxibrCZiiTd5hfqd8sISo63/n5BjEUHQ0zJ7UFlgecI8Dz139LqZT/Dzq/Talk4YqL32fmGyQdDGdguc22fXvfN71u/q1eYxsKMrLi43ytYNNuwewYMV2y4bi4+WVE1Z4+0QieKD52kB8BorJJC2E/w75oro0342TPL8T/9RK0Pc8FIOxbKLy8lrBmtiVnHR46dCHXSovjQrPYCLZ69h/ECSh+NtQpPecRt1RyjG4XsyWAlVe+WsduUardQ+G8s0/LsPqHX1oa0z5lnfWoatHl65VeWUcKq8IRvmswZjM5HL5PWs+Ki/vfSj564de0wcBcYf+19czYmi/igtRNjaOGvzVj59Fr5wo0gFHAKcShCTpV8BAfhD4DWg/byOTMCMCjtMJfUUURbf/5OmnnotN5RVBQrFtKAHtqD5LJUjL5IYzObQ0JDEwnPVwC3e6gppIXf4Sil8VGR9GnTGcJEz2Kdg2MINJWstQXFKQPx2Fu/XvfmkDAKAxlUDfUJDXVA6pBCGTE97Sf4BEZjOUhqRROHcvlZdTmvNjog3JBEay2aJ+M1Lf5QSmdjZj8559mDSmSZvH5Jjh/B6gyoAlFAP0OladRkZ5Hy8vtaL3Ywh+ISJUMV+GIoSRhOKnyzf1qlH3RiKqvEwkFK8+tc+UD2Ao6gNMJsjTKN8ij97VHqCWy4eyMJG6vLrNxA7mt7FxJIRDh7KlBa3qvY4Jdh1KppmWnP3k6eUl2/ALdKjo9UImK2znEM9ncdCi67dhh4RismL3WkiZqLz86B0xsKFkcwKtjUk0JBOu+cedJ3/t9Tx5CaUyLIUZSkgoA7VuICmjfCpBnjppEwnFz4c974rrTaOA8whRf0kGCDbK+0USVh+X1r1VSQ4Gg7kUXXQ+XeY19PJKJxPa9iwJxRLa9Ub5nO/7L8ofJFGZGOU1zCNMlAW/A+GA/ETjPBDKTYdz8iy+b2KnCWIoKt3vsXJCoDGIobjsLMVtOSUUI5VXwF4yP5qzuTxDKdXLK0GE5oYk9nkEiDTaA8Ruw7WNtE846OfW7ET3wDDSyUSgX7/fatpPx6rGuJ/UIAQcXl6e2XwnfHcYEu9nsT3FfM6HMVkcmXkZ6cuaRxu27nu9n+FsXkLR2Smywnlip29TvvTYkqHPgl3RVyj5CSGMV53OSS3I6UF3DDQQvBo3CSNiLzoy+k16yrvON/RKTqAxpd6NVzv5h9QtcIZD2lBMJBTvfSg5e6wUZnGpzDzVohZDaWlIujZkOmG0D0VtKjVchEQFM5SQSHkY2gBga88gBkesgeSlkzY5j8HEc8bTQC3rTZrovCR0k42Jysu5gtOqvGSSic+/d6jw4FWYUw3o7yZtqQLTKfIw2mbRLBmKjsFlXZNE6VKXatpflSgZSkEeZx9s2q2XKpx5G3xUtFa69eulVnG+lyCjvG4Sz/ncby3o66BQNPmwR8GSgy66skprTCcC1XMWXR5eXi6VtA+9Hn3vnBv8niWZIDSnk5675U32Ram+DVI3xgVmKCFh69ALB4ljkKWT5KmTNvFIcnu0uO8FxfJS6Sron6+EYtse9JsS82dllPbxqbsm6iFPjxoDo7xTreevNrGcFdIebt3DmRyaUkkkSD8BWyokc5WXtxpP2p38VIkq9EpBHifdC1bugB/cEpW/tOTFUIJsJEGTuKKXyJrUnHW0Nbl9gvyN8gYqrwDmNpTJIpWw3r/X9+N8Bk+GYiQZ+Km8DGwowjpN01J5eUgojqLeKi+rrR29QxWxozBDCYmUx3kogw6xNJ1M+KxugtUzLoYS0iifK2Qonq346/JzIr/C8rIHOenU6ZvzXl4+RKjyGgkHcDPUoBW/ld9nlSuN1EkPG9dwJoeGVAKppJ7hWDYytQ/Js5lAWnzjn0l4ec+F2aiWzeWQTnl7JTrThzzUUW5pWXPfwXB0q2BVf4s0hDvpb29Ku/KaSyhek6djPGrGk3q/iYR3JAvX5lYPhm9yrkom622/UguJdFIvKQPWt5MkxKDystK7+4exo3dImydOMEMJCS8JxRnvyk/lpVZOvgzFR8fqWo1r6lDfvzor3Bc2c9LcEvmQK57qu4APOIxRftilWtE/Y5ANBfA/Ozsn9dKpZEIrHQxlcmhMJZBKkN6G4rPq1LYX4Ebry1AkfYWTiRdz17cDe1EQtLPfSyUyEiChOBdSOhuJWhA1S2cHpwQwsc061fAz75tp0RiwyAry8nLbUIppUe83mfBWWTr7wUsdNTCcRVM6iBbnWNHT2ZTytuVkc8oon/JUeZlsbFRj7JMnzajIrnlmKCGRSpClEil4gc6VjZ/KS31QxhJKoVHeuRrXqSBkWspWeXm3kz8CWG9Qt1eEnhKKlU4U/TyUEQ81gEm8IiGAyfLI1Q0eu4pVvgSRJ8NQK9hkgjxsKMLRr7r6rURFS1CAUH/vOYu+7r5hbVkT5HL+0bEtGvwZSpC9oF96ILU0JLUqL+XKq5wdCvPMPWgsrvzgoRYtftKlgcrLGa3By224IZWwYu0Zqbz0eQZHshjX0iCvvWwgOU8bivq/MZ3wZygJQnM64anyMvE2y2QFJo9pxHf+5hhM7WzWZ4oRzFBCIi0nHOf8+fyanXj/TU/l8/iovNSA9TPI+onUfi7FzvxJClZ5KegmqZGsQHPae08GkA+219aQ0htkhfvXlwYPV0p3cEhvNd/4VmsS9wv3ba36LMcKr53UDakE0kn9h54V/m64alIeKyebOFRevUMZIzWMvo6cr5oolxP2ZtxsTmgXDs5JXueCria7zua0limpca4YSqEEkE4mbBqHPFQ7gGIoIby8dEb5bA6NqaQM3uoloTj7WtMf2RxGsgJjWxVD8VIVCtvBo/BbV4yqMZX0iUsmkCRCS0MKAyN6+5ZzCHouGHI5W01bCTBDCYkPHj4JCXJvjHvsrW329T++f5bllhowYHPCe4J0rkoL63FKPvpgfVaaicpL7VXRrixzOVus9/LHV6vTjpa0PoihbUMJp/IKK6HkBOzwHQMePvuANUGkkgmkE3qGv3nPPjQkrQWDdqd81mmUL65fTZZq8vQ0ygtzlRdg6b8VCsv49e1IVqBJLgp0+VSft0vjuI4hOCU+3eMoBt7Z0qAfR7bKq1hCGckKpJJkSx66RYldj0NC8fq2nH3jZZS3FgzkvUgKkFDU846TDMXLvpHN5Rdkhf2iGFWTj4SinGIso3zpDivOMVsJMEMJgQ8fOwUHdjYXBRdMOFx0TztsItLJhPbjent7L9bvGrD/9xpMa3f2e+YZyeXsPSZ+EXH9VDMKSR+GkskFSygqCmpnSzpWLy+vfThec2dOWLuKAaBvyHuVOyK9tHTBITd2W+/lnsUbPTem5gIlFKvtVsnctBtGc/m9O0Eqr7EtltF6l2OBUUi33yQ8ks3Z71Cn4skzlLSkv9hFednmvXnaNXWoCXVsq15CsY3yiqE4bBtq42WjLaH4SWw5hw0l2N7j5TbckEygOZ3EYCarXdC5GIrukK4ChuInoShmPpwt/Iat/5vS3hJKNme5uDenvTc2ur8NbyarxmwlwAwlBI44oB0AimJ1pRzSQKNcAekm2Kt+/4brf6/B9M0/LrOvi92T/VedauIwk1CsXy1DyQbbUPrl5D22pQGZnCiiR/WRicrLxVA89uF4MiYBHNBh6Ye37vHem5HJWuJ/KllsI+kbcobX0atEMjlh70PSTkZyQuxotibofRpVhcszLkBCmSxjOG12PFOh+sRrQgOs99rsJ6HI966ku0JPr8LVt15CySCZILQ2pDwkXSmhpFOSXrdbbjJBICI0pBK+eyVyOf/jt1V9umuFoUwOjekEmhqSnueMOJmabkGgJGCl1tQxQWvzac5e5BTaU9X31JhK+B5rnUxYGxsHRvTMz+28o60G+0bye6sqAWYoBjjjsIkAgM+fcQgAa7J2Tm6pghWAlw1lijSKnTZ7AgC9hFI4cArtCVnHykdf3vrNb2z0WQX7nAqoXGTTSfI0luYlFOvjKlzRqfrNJBS9VGK0Csvl0N6UQntTCrv6h7V5FD2pJEmVZKFe2/r/excd63n8gNPIreuSNTv6AAATpPeSbm+HeoaEhyODTWtOYM6UMQCAFVvzUoLqp/Fyhew3CQ9nc77eSKp9tR+kcHIsVB/q+r9/KIuWdBKN6aTnOAKcEopz4ZBXxzQmE56uy1benG1D8Y6sG8xQlIQCQGvsdtGnqUOpvMb7qLwyUgpVkmphv6jvIkhCScrQK57Mzxkd2+Pb2DecteeLSoAZigG+87fH4OmvfcBmHIUMwymhnDBjrHfoFQEcMrEVRx3YAQC468UNRVkeXWbZY7523uEyumqxCN4pV8D9GvWO+oCDNjZmc8Ie6HqVlxW3KpVIeEooalV/gPRqWrfT7WGl+sjIhuJSNTgZiptmHZSeuK0x5XkYkao3rWwkBe9HTS5Kpel1xLOfyusf7lgMIK8C7dPQoiaQpnQSI1nvnf2ZbA5jWxvQnE66nkn1wUkHjwMQLKHYiw9NM6rP2xv1NhTVJ6ccPB6AfiztG86ipTGJxlRCS0umkKEUnIqpDMaNab2aGLDGT04g0IaiJvuGZEIrOQxncmhMJ21aBjT0OpmMbiGVV/F5q7wKJb9CCUV9x03ppKf6Tnl5tfgwP6e05zWOBkey9vNWAsxQDDC1sxkHjW+1/29vSrl09Wp1cMKMTtvop1ut9Q1l0NqYwnK54nzwtc2u+7mcwBd+8woAYNaEVsm4nCoS63piuzWB7+ov3qhkuyQGfHzOD0HvnSMDXWrUQwpqonvP9LEAil1cS1V5uaPGOj+a4nJCKE+lBFobU7ajgA6Wyou0G8oGHYcvWTYwvX0qyHUVAMY0WxOJjrmpyaXFJ8QLYE1mqSShtTGJfuckJ/ujtcHbkA7k1S5+6lFVV94o75601AStVtpalddIFi0NKYxpSuslMgcDLaR3RL4PwPJ4CjrqNqjvB4by9o1BjbSjbChNPpN0ryOsvc7138QoP+x4xzqX+ryE4u28oxYDzT7Mz/m+vITdgeGsLZFVAsxQSkBbY8o18B54ZRMA4K7LTwLgrfLaMzCMsS0N+NAxBwAAZk9ud93vHshPyMdO67AYk2NQKyagJJw3N/cUtaFWzo1pNWH5e2gBepVXJiuQTCS06iGFrr4hpBKEKZ2Wrr/w47IllLBG+Zx70lET3qsb92DmVQ/j1Xd3O+5bdTekLIbia5TPWnppndSlJpemdAKtDUmtt9jgSDa/6vQxqKs8//6HZUX3lNSiJiRd3w4MZzCcsSbb5oYkBhyMqbC8lyF7KJNDTuSZhU4loibwNg8vrzxD8fZaW7yuG6kEobMljb6hjOfOftUnzoWMUkECFrPwYiiqTJuHJKWgxvT4tgYtsxjKZNGYStgRpfUMJd/XusleGcjzRnm9ezIA2wV9qIih5KyYcklvG4rahNls06obj/4en4C0oTBDqW20NaawcFWXLWZu3ztkpwPWR5wT7tUOAOzsG8b4tgZ84r0zMGNcC5xm8+7+Ycz99hP2/2ql7PxA1f1jp3WgvSmFtVJnryCEsE93U6s5L4Yy74YF9rWXyiudVKH49YN1a88gJrU32s9dtKtbfiyWTtmfqThX4c4PeTiTs43T9y3ZCACYv3y7fV89XypBaGtM+qq81L6HlMYtWK32mtNJtDWltOqqoUzOnnz97B/5/QeiSCWipFO1b0YXYeAa6ZSRTiZwYEcz3u3OewbukGNt+rgWANCuxAHgl8++AwAYIz249DYUK8328iqYHBXtqo5C5re1Zx+27R3E6h19tq1m8bpuVx7FlNRmz72OSAaZXA5JqfJqSCU896GoUzAntjeCyFvN9z/PvGPn26eZ6PuGMtaqX0komnoUQ+lo1nsuFkoofiqvhmTCOmSr4B2P5ATSCW/3dFVvYyphqyN10l//UMb+zr36ZJCN8rWPF9/ZBcCa2G57bh0Atx3loPHWx+50uezZN4LNe/bZxrwJbQ2u/QVf+ha5ulIAAByOSURBVO2rrjYa5epGDWpnSJH2phTGtjQUhRl5fVMPfjD/bQDwdFkEinX7OobSP2QNxEK1m4IQAss29+CAjib7Ay3cVJjJWq6P2Zzw9NdXWOkwPGcKVF6TpIpvqwyx7tQJq7zpZAKtDSn0eQQ5BKyPqylteXl52VCa0km0NRbXI4TFHFobUlqD+uL1+Yn0pFnj7Qm2qyB+0ufuXgoAtlTX1VccNv6BVzbLNi0Dv/NY4e29Vv6DJ1oq2L0eoWa+/7g1DvKeesXv8Fdy7CoppjBciW2Abmtw/a/g7KOVW3sBAF+597WCOqw8B0qHFKfTxMBw1o443JhK2GrHQvQMWM/Y0ZxGi0/03S1yfIxpShdNsN39w9i+dwiHTW6zFwW6vusdHEFzOmmpozR9tmZHH9JJwvSxLZ7M7R1pSzygowkNqURRv1r7ocjeUKrDkFR5KVp1C5zfLdloS2tewT33scqrftA7mMF1f14OALj8tFl2+tyZlsFUrUZzOYHj/t/jAPIuoONaG1wf18ptva66m9JJK8y6HNTOwX/UgR0Y05wqGkS7HSozNUno9MCF4Um6NZ5RuweGMb61QbpzFn80G3YN4O3tfThnzgH2Cqjw4xrJ5WyPpz0D3jG2Hl22FWu78jQ53W2HMjm0NqbsiQeAy2tFGTiVUV734QHAvUs24qV3utHelJbBH/Wr8eZ0Eu1NKfQW1JNxGIbTyUSRmvD259fZ14dOasNPLjnB97mVp99bW/Zq7wNq34tb7blj7xDaG1O2TS8o4F9jyvJq6hsqpuNPr1vS7KR2i7l5eXmpd1joBJLf1JjG354wDQDwvkMmuPIoiXFKRzMmtTdiuXzeTDaHgeGsLR3NGN+KVdv2aiVZdVxvR3PaM7aVen/5vRvuPMphZPbkNhzYWeyOrdA7mEF7U0rG4yv+djbv2YcpHc1obkiiKZXUMpTd8nua2tkswzAV2OtGLGbhJaGoxUtjKmF/x14MA1Aq+OL763f2o3+YJZSahwpm54xooGL7AJZLYXtjCu/KiXvJhrzO/2TpMTOutQHdDqP6zj73xGBt+MpHGlUf1c2XHI8DOprQ3pguOhd7u+OQJOUJphPb1+8ccP3/3Jqdrv9vWbgWA8NZjG1twPjWhiLa9gwM4wPfXwgAOH5GZ95rxqm2yubsFTbgZnZOrN7ei6/e596f093vMIxmrQ1tyqsGyIdkB/LSVTrpb5T/+v1WGxt29SOdoCKpTKncmhuKJZSB4Qz+IKWGpnRSq8YoRH6jpZ6eIw4YA6JizzgnxjSl0ZByqz3veGE9eocydmBFpQJzwjkpHz3VUo/u3Zdx3XdOhEdPtdyTvby8JtgSivtZ1LP94lMn4tRDx4MImDbWHS9KvcvOljQO7Gy2x8EvpHpKTZjHTevA9r1D2NlXPE7UYmpMcxotDfpJ/NLbXgYAfPHM2dbu8oI8qo7Olgb7W9Ux+zc292BKR5Onp+a2nkFbfadrBwD+8uZWAJYzg86eqgzlyQRpbSgDw1lkcgIdzWlbnVwoMavx+/kzDsGYplTRXAAAn7nd6hOdN2i5wAylBHzlnMOQIGDVtrwNI+lQeRERZoxvsXfFr9uZz/fpkw8CAIxrbUR3/zDmL9/uGS6kszltq7XUhDBBTqztTcWrkm178wxFTcA66WO9ZHRqnwOQZ1jbegZx06MrLRpbGjBpTGPRKvg3i961r6ePa0GT3B9wz8v5dDWI1QTT47FSv+C/n7cnpt/+08mS5nx7gyOWcdK5cN3lYHBqUuhsSVsMxceGAgBru/rR2dKAXocBecGK7fjjq5vR3piSKq809o1kbXXEj59Yja/LTamN6QTSqeJJotA4qyaCny1cI+kcxgclE25rTOHIKWMwtbMZ73QVM5Sk9ES75KQZrgnJySgaUglMbG/Eo46wPwCwcttezPnWY/b/iqH0OiSUHzz+No745qMAgL+bOw2dzdZYKZTK1AJB2Wuc4wvI2wjbm9IgsjY3FjLQbXsHQWTZNQ4Y04SV23ohhMD3HlsFAJjS0ST7x6rrmj++WdQfSuU3piktD5wqfscvS9tNTu7TKpzo9+yzvoOOZktCbWlIFk3C23oGsWLrXhw/YyxSieLNyWu7+rCxe8C26XU2p7FbM64flza+1saUJeEXMZSMPC9ezxz37MuPaSXBFdKqvuvp45rR3pQustcCedvYKYeML7pXLjBDKQEdzWnMPWgcfv70WjvtsAKPrZnjW23V0gbJWJZcc7YdP2tCWwNGsgL/dNcSXKPxBgIsKebd7gHkcsIeMGOk5NHelC6SHNQxrqkE2XacXzzzjivW2C+eXovvPbYKB3Y04S9fPg3fvvBoAPkJ33m2+LjWBkxqbyqyAzj1vgd2NNm78p1hZV58x5J6Dp5gqWa2ehwxqyatgye04j3TOwHk9eyL3tmFzXv2oa0x5VJPOEPTqBXv2JYGaUgVdhgVwJqEb3h4uf3/f33iPZjY3gghgPuWWN55l9+5BKt39Nm2gkLpIuWIhdSUsoy6hR9w4cRw0DjruZ9dvRPrdvbjc3cvtaWRz51+MADg6AM7sGR9t4tRPPjaZmRzAt84/wi0yRWukhzUouKaDx0JwLLPrNi61/V+r/vTcnsy/eaH5wCAnHDyk/BPn1pjX3/tvCMwqb0RLQ1Jl5PHL599x1bnWi7shI3d+Xcwks3h879+RdafsvttoGA1vL1nEBPaGpFOJvDeWePQ1TvkmoSVenj25DYAxRLbpt0D+Op9rwOw3LGbG4ptKN97bKV9/Xdzp6M5bUU+Vp58T63aga/8zqrjAMkMxmgm4ZP/03JUOfGgsRjf1oA1O/rsd7NvOIuzfvA0tvQM2nWMb2twLW4GR7L48ROr7f9b0kl0NKexp0A637BrQEoflkt4oZTitBm1NVo2u0Jp6q0tlofnhLZGdLToGdvgSBbvO2Q8zpkzueheuVDXDIWIzieiVUS0hoiuqmTbf3/yDPt6amczTpe76RVmTrAklKUburGhewAzx7fY6h8gb1QFLPffhlTCnmgUzj3qAGzavQ+vbtxjTwjq402Q5TU286qH8daWHgghsG3vII6cMgYv//vZtl4cAF5wqLT+8xHr41NGUrWn5V7pQaUkogltDTjlkPGY2N6I3sGMSye9fmc/UgnCon87y2aQsya04nDJVIUQ9mRz2OR2dDSnsWBl3jMLAO5fugkzr3oYAHDmEZPw5Fc/YOt6fyQdCz5x60sALIPvzz91AlobkpgxrgXrHKt6pZMf19qA0w+z9PenffcpWyXw0Otb8D/PWvaNr513OC48fqr9zP/2hzexcFX+1EP1ftoLDKFOo2ZX3xAO6GjCM6t34udPr7WZbaGDREdLGl8881AkE4QPfn8hFjm8nw6XIXzeO2sctvQM4rTvPoVNuwdwxV1L8OV7LKO2Yq5OldcrGy3VqVohKyhDP+CO+6Xeh6Xy0kuIE9sbkUgQDp3UhrVdeYby7YdXALBsU83pJA7sbMam3XlG/e0/55m0koZbG1PoK5Aenluz01a/zppgLXIWSaeW46Z32u/io8cdCAA4eEKbq/wNkg7AkuxaChjKjr2D+O+n8gu7GeNb7IWBUp/9y+/yjgJqT82Y5hTuXbIJc789HwDw2sY9dp6jp3bgzCMmYfWOPnui3rwn/+wTJM0T2hpdKrp7Xn4XP3rCGrvHz+hEIkGY2NboWpDdt2Qjlm/di5MPHm8b3J17TF55dzf++uZnAQAdzQ1IJghTOpqL7D2vb+qxaZ0xrgXrCxjx3sERrN7RZ29KrRTqlqEQURLAfwP4KwBzAFxCRHMq1f4F75lqX+uMXodMtD6Mj93yIh5+Y2vRhHP01A77evWOPgxncjhkUpvL+KwMtx+75QX8q1yljZeT3vEzxtr5PnTzc7j1mXewcFUXpnY22y6Nj3z5NACWQe/Z1V0udcSNHzsWAHDyLGvA/eTJNVi8vhs7pFrjD//nVLQ3pe0zFK5xxBdbtK4bpx820TWxHXFAO7btHcRPn1ztGvwfOnYKzjpyEl5etxvbegYxKOMSzV+eX1X/x0eOcvXNSFZgu0O9cu5Rk3H+0VPw1nXn48wjJmFNVx96B0ewdEM3bviLNeGMa22w9+cAwBLpdfXCml122kxpyFaTGAB85vbF9vV5R1n7g5SqYFffEDbs6rc9pgBLdbV0w2509w/jxkdW4r03PAEhBHb3D2P6uGa8/O9n2XlPPni81otHqZCU99qm3fvwX0+stlUlADDnQEsd2ZxOYnAkh3uXbMQbG/OTCAB87oz8AuSiW17AH17dhAnt1rtvSidwlKxj2tgWvL6pB1v27LMXBpPHNOL6C/L93pRK4tnVOzE4knWpei46cRqILK+mtV39NjP4tVR7njZ7gq3e62hO2yq8zXv24TeLNmDznn1YLSWfWZJZqM27q7fnHVGICB857kAsWLkdgyNZvLh2F3b1DbkYZCqZwCET27Bsc4+t2nxjU34v1pP/eoZ83mbZrwMYymS1q3f1jnf2DWNgOIML//t5q40EYdaEVpvWxVKC/NnCPNO6UH77k9obsWXPPghhucWrxdrkMY2493OnWHnG5BlKLifwNWnLmzdrnB1/7neLN9p1f+P+vD3xELnonNDe6HLg2bxnH25esBqzJ7VhamczDpnYhh29Q7bElc0JnPHdpwAAx0zLfxOVQCo4S81iHoA1Qoh3AICI7gFwAYDlvqVixBP/cgbO/uHT+Lu504runX/0AfiXe1+3/1cThMKk9iZcd8FR+NaDb9lpx03rxHPfONNWWUwe04SJ7e4Vjvp4L5k3HSu27sXdL1nhW9RgVuoaADhyyhicMKMTD7y6GQ+8utmOSfbdjx2LQydZH0xHSxrnHTUZj721HR//+YsArFWpkmDOP/oAHPXsGPz+lU34vdzACQCffX/eqw2wxP+efSP4/uNv2xPww196P5rSScyZMgYPvLLZVik4cffl8zBDqucA4PsfPw5fve91XPI/lnRy498eg4vn5aXBc+dMxh0vrMfJ31lgS0fHTO2wGe3nzjgYv3j6HXzyl4swqd1t/zl7ziQA+ZMCnfjjlafaUoFaDHz0p88X5bvgPVMxlMm5wubM/fYT2NU/jEtPOcglGao+duLWT5+II6XtStkPAEtiU3j538+yN9/NnWktHL7umGhmSTXiN847Ar942jJuL9mw23b+uGTedPzn3x5r51dqx/fd+KSddu1HjsJfHzPF/v8VuVn0iG8+aqtBAeAy6YAypaMJz63ZiU/c+hIa5HkxR04Zg7vlZl4AOH76WNz2/Dpb8lT45aVzAQAzxrW40p/66gdc/587ZzL+9PoW277jhHKE+chxB+LulzbgqGsfc91vTidxsHxv08Za7VwkxzNgTf6L/i3P7I84oB1LZX85bU5//tL7AeSlKaf0BwCvf+tcdMgo0JM7mjAwnMWsq/+CL501G0OZHE6aNQ6/+sx7bceRKR3N2D0wUtQn7505zvaqu/7PyzEwlMFjy7fZzPfAjiZMkgu2BAHPvN2FM3+wEJ89dRa+9aC1uHuftI0oxnP1A2/inDmTsXnPPpuJnnDQWFQSdSuhAJgKYKPj/00yrWI4dFIbXv3mOfin0w4uutfSkML8r5xu///Dv3tPUZ5LT5mJT0nV2XcvOhaHH9COsa0N9mQOuD+6my853r4mIlx/4dF45msfdE1cn3RMvgBwkkPkffrtLgDAiTPdg+zmS47HJ+ZOdzxXu+1k0JRO4lsfLhb8LjrBzUQ/ffLMojxHHmBNnB+R6oxCXPze6ThttltV+KFjpuDAjiZ7pfv+2W431BMOGoupnc3oH86ibyiDca0N+NMX32/f/9q5h9tRoRUz+dCxU7DuP//aDi44tbMZ/3DqTFsCBPIqJiD/gTpx7UfmYP2NH8LRUztw3QVH4zrH6l6tHufNGucqM6ndzbjOPnISzpVSEGDp6o+Z6l5BXnfBUS6m9P5D3c/f2ZI/hz2RICy95mwc56AdAD7ueJcAcGpBHYD7eQH32FLS6IJ/PQNHyHfopFu5TBeeAPjJk9xjDwCuOP1gnC11+MkE4fVrz8X7DhmP3/zjSUWqu48cdyAuPeWgojoe/8rp+I+PWv393pljbfugwoeOmYJXvnmO/X/hfQC493On2AsQq55xRXnu+/wp9vPOHF88Bv7lnMNsZgIAHzk2P65vXmDZTr4ubV8KHz52ir0AUHjtW+egKZ1ER0safy2jZvxg/tv2vrUrTj8YjzrmDkXTO139uOaPyyAAfPXcw/C1848AkP9G/vzGVnz5ntfw3UdXYWpnM5Zec7a9KbVSoKAdzLUKIvo4gPOEEP8o//80gHlCiC8W5LsCwBUAMGPGjBM3bCgOyFhOvLh2FwaGMzjrSL1hLCeDNEaJCJrJ5vDEih0484hJ9kY2hWxO4Pbn12FCWyN+/vRa/NNpB+NjJxZLVIDlIfTi2l248D1TXW66qo2UXJk6zwUpxMbuATSlk5jQ1uD6gAHLlvPIsm3ICYEvnjkbk6T+Xvc8d7ywHvNmjcOx0zqL7vcNZfDi2l14atUOXPze6UV5evaNYFvPIF5YuxNjmtKez5vNCTz85lacecQk1yQAwA59smp7L9qbUvaqtxCPLtuG/3joLfzwE8fhlIPHFz2zEAJvbOrBT55cg+/8zdH2qrOQ3hfW7ETvUAZ/e/zUoujVA8MZbNg1gDc39+DYaR32BKPrl3e6+rR9NjiSxeY9+7BgxXacfthEzzp+/MRq/OiJt/GdvzmmiEF09Q5h8fpu/PTJNZjS0YRv/83RmNLhZiqbdg/gode34H2HTMDgSBbHz+i0Gbkplm3uwbjWBvzy2XX4yHFTXOpdwIpq0DMwgvuWbsLpsydq1Tpbe/Zh7z7Lk2/WhFbbdqKQzQk8umwbxrc1YMOufhw5ZUxRv2VzAjv7htCUSmJ77yBmjm8t+r4Aa8y/9M4udPUN4fOnH6Id09v3DuL5NTsxb9Y411hav7Mfd7+0Aa0NSZx+2ERMaGvEQeNbXOMolxP43uOrcMjENnQ2pzGurQEnFPTJ6u29eHb1Tqzt6sMhE9twzpzJtmo1DhDRUiHE3MB8dcxQTgHwH0KI8+T/VwOAEOI/vcrMnTtXLFmypEIUMhgMxv4BU4ZSzyqvxQBmE9EsImoAcDGAh6pME4PBYIxa1K1RXgiRIaL/C+AxAEkAtwkh3gooxmAwGIwyoW4ZCgAIIf4C4C/VpoPBYDAY9a3yYjAYDEYNgRkKg8FgMGIBMxQGg8FgxAJmKAwGg8GIBcxQGAwGgxEL6nZjYykgol4Aq3yydADo8bkPADMAvOtz36SOoDxx1AHUD61BdMbVzmjrV4BpLSUP01qMw4UQ7QF5YEfKHA1/AJYE3L/VoI6uGOrwzRNHHfVEaxCdtURrPfUr01rW5xlttPrOneqPVV5u/Mkgz56A+yZ1BOWJow6gfmgNojOudkZbvwJMayl5mNYSMdpUXkuEQTyactdRKdQLrfVCp0I90cu0lgejjVbTOkabhHJrjdRRKdQLrfVCp0I90cu0lgejjVajOkaVhMJgMBiM8mG0SSgMBoPBKBNGPUMhotuIaAcRLXOkHUdELxLRm0T0JyIaI9PTRHSnTF+hzmCR9xYS0Soiek3+TaoyrQ1EdLtMf52IPuAoc6JMX0NEN1PhqVC1RWsl+nU6ET0l3+lbRPRlmT6OiOYT0Wr5O9ZR5mrZf6uI6DxHeln7NmZay9q3YWklovEyfx8R/bSgrprq1wBaa61fzyGipbL/lhLRmY664u1XE1ew/fkPwOkATgCwzJG2GMAZ8vqzAK6X158EcI+8bgGwHsBM+f9CAHNriNYrAdwurycBWAogIf9/GcApAAjAIwD+qoZprUS/TgFwgrxuB/A2gDkAvgvgKpl+FYCb5PUcAK8DaAQwC8BaAMlK9G3MtJa1b0ugtRXA+wF8HsBPC+qqtX71o7XW+vV4AAfK66MBbC5Xv456CUUI8QyA7oLkwwE8I6/nA/iYyg6glYhSAJoBDAPYWwk6gdC0zgGwQJbbAct1cC4RTQEwRgjxorBG1F0ALqxFWuOmyQtCiK1CiFfkdS+AFQCmArgAwJ0y253I99MFsBYWQ0KIdQDWAJhXib6Ni9Y4aYqLViFEvxDiOQCDznpqsV+9aK0ESqD1VSHEFpn+FoAmImosR7+OeobigWUAPiqvPw5gury+H0A/gK2wdp5+XwjhnDRvlyLuN8uhRgpJ6+sALiCiFBHNAnCivDcVwCZH+U0yrRZpVahYvxLRTFgrukUAJgshtgLWRwxLegKs/troKKb6sKJ9G5FWhYr0rSGtXqjFfg1CrfbrxwC8KoQYQhn6lRmKHp8FcCURLYUlUg7L9HkAsgAOhKU++FciOlje+3shxDEATpN/n64yrbfBGiBLAPwXgBcAZGCJtoWolKtfWFqBCvYrEbUB+D2AfxZC+EmeXn1Ysb6NgVagQn0bglbPKjRp1e5XP9RkvxLRUQBuAvA5laTJFqlfmaFoIIRYKYQ4VwhxIoDfwtI7A5YN5VEhxIhUzTwPqZoRQmyWv70A/heVUytoaRVCZIQQXxFCvEcIcQGATgCrYU3c0xxVTAOwpbDeGqG1Yv1KRGlYH+dvhBAPyOTtUi2g1C47ZPomuCUo1YcV6duYaK1I34ak1Qu12K+eqMV+JaJpAP4A4FIhhJrPYu9XZigaKK8MIkoAuAbAz+WtdwGcSRZaAZwMYKVU1UyQZdIAPgxLvVM1WomoRdIIIjoHQEYIsVyKwr1EdLIUxS8F8GAt0lqpfpX98CsAK4QQP3TcegjAZfL6MuT76SEAF0s99CwAswG8XIm+jYvWSvRtCbRqUaP96lVPzfUrEXUCeBjA1UKI51XmsvRrFIv+/vAHa6W8FcAILI59OYAvw/KceBvAjchvAG0DcB8sw9ZyAF8TeY+PpQDekPd+DOlJU0VaZ8KKrLwCwBMADnLUMxfWIF8L4KeqTK3RWsF+fT8sUf8NAK/Jv78GMB6Ws8Bq+TvOUebfZf+tgsMzptx9GxetlejbEmldD8uZo0+Omzk13K9FtNZiv8JavPU78r4GYFI5+pV3yjMYDAYjFrDKi8FgMBixgBkKg8FgMGIBMxQGg8FgxAJmKAwGg8GIBcxQGAwGgxELmKEwGDUCIvo8EV0aIv9MckRzZjCqjVS1CWAwGNaGOCHEz4NzMhi1C2YoDEZMkIH6HoUVqO94WBs4LwVwJIAfwtoYuxPAZ4QQW4loIay4ZacCeIiI2gH0CSG+T0TvgRVJoAXWprPPCiF2E9GJsGKfDQB4rnJPx2AEg1VeDEa8OBzArUKIY2EdbXAlgJ8AuEhYMcxuA3CDI3+nEOIMIcQPCuq5C8A3ZD1vArhWpt8O4EtCiFPK+RAMRin4/+3dMS6EURSG4fcLGo1KawWWwAIsQUTENiyDhkYiNKLViGo6GxCdxgIQob2K/04kk0FMTjKK9ylPcXL/6su5+XOuE4pU66l97Uu6AA4YHjW67VvMFxhW0oxdTjZIssIQNKNeOgOuptTPga36T5BmY6BItSZ3Gb0B9z9MFO9/6J0p/aV/wysvqdZaknF4bAN3wOq4lmSpv0vxrdbaK/CcZLOXdoFRa+0FeE2y0es79ceXZueEItV6APaSnDBsfT0CboDDfmW1yPCI2P0vffaA4yTLwCOw3+v7wGmSj95X+jfcNiwV6X95XbfW1ud8FGkuvPKSJJVwQpEklXBCkSSVMFAkSSUMFElSCQNFklTCQJEklTBQJEklPgFtv8n7ZLIz8wAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + " sorted_data['inc'].plot()\n", + " " + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": {}, + "outputs": [ + { + "ename": "SyntaxError", + "evalue": "invalid syntax (, line 1)", + "output_type": "error", + "traceback": [ + "\u001b[0;36m File \u001b[0;32m\"\"\u001b[0;36m, line \u001b[0;32m1\u001b[0m\n\u001b[0;31m \u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n" + ] + } + ], + "source": [ + " " + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [ + { + "ename": "SyntaxError", + "evalue": "invalid syntax (, line 1)", + "output_type": "error", + "traceback": [ + "\u001b[0;36m File \u001b[0;32m\"\"\u001b[0;36m, line \u001b[0;32m1\u001b[0m\n\u001b[0;31m \u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n" + ] + } + ], + "source": [ + " sorted_data['inc'][-200:].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "metadata": {}, + "outputs": [ + { + "ename": "SyntaxError", + "evalue": "invalid syntax (, line 1)", + "output_type": "error", + "traceback": [ + "\u001b[0;36m File \u001b[0;32m\"\"\u001b[0;36m, line \u001b[0;32m1\u001b[0m\n\u001b[0;31m \u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n" + ] + } + ], + "source": [ + "" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 19, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztvXl8nOV57/29ZkajXbJWb7Kx8RIwEExwDIRsDQkm6WmgDbTmkwa3pYc0IT3J25yeQE/fl4SU05CTJi1JQ0ICDdA0QFZICqEOBAhLDIaw2QYsgxfZRpKtXaPZr/eP537kkTSSRtJoGen6fj7zmdE9z33PPY9m5vdcy33doqoYhmEYRi4EZnsChmEYRuFgomEYhmHkjImGYRiGkTMmGoZhGEbOmGgYhmEYOWOiYRiGYeSMiYZhGIaRMyYahmEYRs6YaBiGYRg5E5rtCeSb+vp6XbVq1WxPwzAMo6B49tlnj6lqw3jHzTvRWLVqFTt37pztaRiGYRQUInIgl+PMPWUYhmHkjImGYRiGkTM5i4aIBEXkdyLyC/d3rYhsF5G97r4m49hrRaRZRF4VkS0Z7WeLyEvuuZtERFx7sYjc7dp3iMiqjD7b3GvsFZFt+XjThmEYxuSYiKXxaWBPxt/XAA+p6jrgIfc3IrIB2AqcBlwEfFNEgq7PzcBVwDp3u8i1Xwl0qupa4GvAjW6sWuA64BxgM3BdpjgZhmEYM0tOoiEiTcDvA9/NaL4YuN09vh24JKP9LlWNqeobQDOwWUSWAlWq+pR6m3jcMayPP9aPgAucFbIF2K6qHaraCWznhNAYhmEYM0yulsY/A/8LSGe0LVbVowDuvtG1LwcOZRzX4tqWu8fD24f0UdUk0A3UjTHWEETkKhHZKSI729vbc3xLhmEYxkQZVzRE5L8Bbar6bI5jSpY2HaN9sn1ONKjeoqqbVHVTQ8O4acaGYRjGJMnF0jgf+LCI7AfuAt4nIv8OtDqXE+6+zR3fAqzI6N8EHHHtTVnah/QRkRBQDXSMMZZhGAXIb/a288ax/tmehjEFxhUNVb1WVZtUdRVegPthVf1T4D7Az2baBtzrHt8HbHUZUavxAt5POxdWr4ic6+IVVwzr4491qXsNBR4ELhSRGhcAv9C1GYZRgHzmrue55bF9sz0NYwpMZUX4l4B7RORK4CBwGYCq7hKRe4DdQBK4WlVTrs8ngO8BpcAD7gZwK3CniDTjWRhb3VgdIvJF4Bl33PWq2jGFORuGMUuk0kpHJE4knhr/YGPOMiHRUNVHgEfc4+PABaMcdwNwQ5b2ncDpWdqjONHJ8txtwG0TmadhGHOPrkgcVYgmTDQKGVsRbhjGjNAZiQMQTaTHOdKYy5hoGIYxI3T0JwCzNAodEw3DMGaEjn5naSTN0ihkTDQMw5gRfNGImaVR0JhoGIYxI5yIaZhoFDImGoZhzAiD7ikLhBc0JhqGYcwInYMxDbM0ChkTDcMwZoTj/eaemg+YaBiGMSNkrtPwqgQZhYiJhmEYM4If0wCIWdptwWKiYRjGjNDZHycU8HY7iFkwvGAx0TAMY9qJJlL0x1Msrirx/rZgeMFiomEYxrTjxzOWLXKiYcHwgsVEwzCMacePZyytLgVsrUYhY6JhGMa00+mKFS5b5IuGWRqFiomGYRjTzvH+GGDuqfnAuKIhIiUi8rSIvCAiu0TkC6798yJyWESed7cPZfS5VkSaReRVEdmS0X62iLzknrvJbfuK2xr2bte+Q0RWZfTZJiJ73W0bhmEUHF0Rz9JYMhgIN/dUoZLLzn0x4H2q2iciRcDjIuJv0/o1Vf1K5sEisgFvu9bTgGXAr0Rkvdvy9WbgKuC3wP3ARXhbvl4JdKrqWhHZCtwI/ImI1ALXAZsABZ4VkftUtXNqb9swjJkk5rKlasrDgFkahcy4loZ69Lk/i9xtrOWcFwN3qWpMVd8AmoHNIrIUqFLVp9RbDnoHcElGn9vd4x8BFzgrZAuwXVU7nFBsxxMawzAKiLizLKpKigATjUImp5iGiARF5HmgDe9HfId76lMi8qKI3CYiNa5tOXAoo3uLa1vuHg9vH9JHVZNAN1A3xljD53eViOwUkZ3t7e25vCXDMGaQeDKNCJQXBwFb3FfI5CQaqppS1Y1AE57VcDqeq2kNsBE4CvyTO1yyDTFG+2T7ZM7vFlXdpKqbGhoaxnwvhmHMPLFUmnAwQEmRJxq2uK9wmVD2lKp2AY8AF6lqqxOTNPAdYLM7rAVYkdGtCTji2puytA/pIyIhoBroGGMswzAKiHhymGiYe6pgySV7qkFEFrnHpcD7gVdcjMLnD4GX3eP7gK0uI2o1sA54WlWPAr0icq6LV1wB3JvRx8+MuhR42MU9HgQuFJEa5/660LUZhlFAxJNpwqEAJSHvJ8cW9xUuuWRPLQVuF5Egnsjco6q/EJE7RWQjnrtoP/BxAFXdJSL3ALuBJHC1y5wC+ATwPaAUL2vKz8K6FbhTRJrxLIytbqwOEfki8Iw77npV7ZjC+zUMYxZIpDzRCAUDhAJilkYBM65oqOqLwFlZ2j82Rp8bgBuytO8ETs/SHgUuG2Ws24DbxpunYRhzF9/SACgpCpqlUcDYinDDMKaduAuEA5QUBSwQXsCYaBiGMe3Ek2mKnGgUh4LmnipgTDQMw5h2YkPcUwFbp1HAmGgYhjHt+IFw8GMaZmkUKiYahmFMO/FkmuJM0bCYRsFiomEYxrQTT52IaZQUBSx7qoAx0TAMY9rxV4SDBcILHRMNwzCmnfiwQLiJRuFiomEYxrSTSOkJ0QjZ4r5CxkTDMIxpJzPltrgoOLgpk1F4mGgYhjHtxJOpoSvCzdIoWEw0DMOYduK2TmPeYKJhGMa0k5k9VRIKkkwryZRZG4WIiYZhGNNKKq2klSHZUwDRpIlGIWKiYRjGtBJ34pDpngLbva9QMdEwDGNa8UUjc0U4mGgUKrls91oiIk+LyAsisktEvuDaa0Vku4jsdfc1GX2uFZFmEXlVRLZktJ8tIi+5525y277itoa927XvEJFVGX22udfYKyLbMAyjoIilPHEYTLkNeZZGIqWzNidj8uRiacSA96nqmcBG4CIRORe4BnhIVdcBD7m/EZENeNu1ngZcBHzTbRULcDNwFd6+4evc8wBXAp2quhb4GnCjG6sWuA44B9gMXJcpToZhzH18S6PYWRq+xZGwQHhBMq5oqEef+7PI3RS4GLjdtd8OXOIeXwzcpaoxVX0DaAY2i8hSoEpVn1JVBe4Y1scf60fABc4K2QJsV9UOVe0EtnNCaAzDKAB8i8K3NEJBAU6IiVFY5BTTEJGgiDwPtOH9iO8AFqvqUQB33+gOXw4cyuje4tqWu8fD24f0UdUk0A3UjTGWYeSd3+xt55bH9s32NOYdwwPhYbM0CpqcRENVU6q6EWjCsxpOH+NwyTbEGO2T7XPiBUWuEpGdIrKzvb19jKkZxuj89LnDfPvR12d7GvOO4YFw/z6ZtphGITKh7ClV7QIewXMRtTqXE+6+zR3WAqzI6NYEHHHtTVnah/QRkRBQDXSMMdbwed2iqptUdVNDQ8NE3pJhDBKJp+iPJ2d7GvOO+LBAeJFzTyXMPVWQ5JI91SAii9zjUuD9wCvAfYCfzbQNuNc9vg/Y6jKiVuMFvJ92LqxeETnXxSuuGNbHH+tS4GEX93gQuFBEalwA/ELXZhh5J5JIEU2kbaVynon57qmgH9Pw7uN2nguSUA7HLAVudxlQAeAeVf2FiDwF3CMiVwIHgcsAVHWXiNwD7AaSwNWq6idkfwL4HlAKPOBuALcCd4pIM56FsdWN1SEiXwSeccddr6odU3nDhjEaA87KiCRSVAVtCVO+GD2mYe6pQmRc0VDVF4GzsrQfBy4Ypc8NwA1Z2ncCI+IhqhrFiU6W524DbhtvnoYxVfpj3rVNJJaiqqRolmczf/DFwd8jvCjk3FNmaRQkdjllGI4Bt0LZ4hr5ZbRAuIlGYWKiYRiOiBOL/piJRj4ZHgg391RhY6JhGI5I3FkaMauJlE+GxzT8xX1maRQmJhqGAagqA040IuaeyivxYdlT5p4qbEw0DAMv/dNfbNZn7qm8Eh9WRqTI3FMFjYmGYcCglQEn3FRGfhhuaVgZkcLGRMMwGCoUFgjPL8NjGrYivLAx0TAMhoqGWRr5JZ5KEQwIwYAnFv69WRqFiYmGYTDUPWWWRn6JJ9ODLikAESEcDAzGOozCwkTDMBi6oM8W9+WXREoHXVI+RUGxGl8FionGPOeRV9u44T93z/Y05jxDAuG2TiOvxJJpwqHgkLaiUMDcUwWKicY8JpVWvvDz3XznN2/Q0hmZ7enMafw4RjAgZmnkmXgyPVh3yicUMPdUoWKiMY/51Z5W3jjWD8Bjrx2b5dnMbfwFfbXlYVsRnmfiqfRg5pRPOChmaRQoJhrzmFsee52mmlKWVZfw6Gtt43dYwPjFChsqis3SyDPxZGpIIBw895TFNAoTE415ykN7Wnn2QCf//V0n8563NPJE83G7shsD37poqCy2mEaeSaR0sBy6T1EwYCvCCxQTjXlINJHi8z/fxdrGCi7fvJL3rG+gL5bkuQOdsz21OctAPIkI1JWHrYxInhmecgueaNjOfYVJLtu9rhCRX4vIHhHZJSKfdu2fF5HDIvK8u30oo8+1ItIsIq+KyJaM9rNF5CX33E1u21fc1rB3u/YdIrIqo882EdnrbtswxuXfntjPoY4Brr/4NMKhAO9YWwfAswdNNEYjEk9RWhSkvDhkBQvzTDw5MqZRZDGNgiWX7V6TwGdV9TkRqQSeFZHt7rmvqepXMg8WkQ1427WeBiwDfiUi692WrzcDVwG/Be4HLsLb8vVKoFNV14rIVuBG4E9EpBa4DtgEqHvt+1TVfv3G4IVDXaxpKOcda+oBqCopIhwM0Bu1H8PRiCRSlIWDlBUH6bcV4XkllkpTHR66E6LnnjLRKETGtTRU9aiqPuce9wJ7gOVjdLkYuEtVY6r6BtAMbBaRpUCVqj6lqgrcAVyS0ed29/hHwAXOCtkCbFfVDicU2/GExhiDtt4oi6tKhrSVFQdtpfMYDMRTlIaDVIRDxJNp+0HLI9ndU2IxjQJlQjEN5zY6C9jhmj4lIi+KyG0iUuPalgOHMrq1uLbl7vHw9iF9VDUJdAN1Y4xljEFbb2yEaJSHQ5ZKOgb9sSTl4RBlxZ7xbcHw/JFIpQlnDYSbMBciOYuGiFQAPwY+o6o9eK6mNcBG4CjwT/6hWbrrGO2T7ZM5t6tEZKeI7Gxvbx/zfcx3VJW2nhiNlcVD2svN0hiTgYRnaZSHvZXLlnabP0YLhJtoFCY5iYaIFOEJxvdV9ScAqtqqqilVTQPfATa7w1uAFRndm4Ajrr0pS/uQPiISAqqBjjHGGoKq3qKqm1R1U0NDQy5vad7SPZAgnkrTMEI0QvZDOAaRuBfTKPctDTtXeWPUQHjS3FOFSC7ZUwLcCuxR1a9mtC/NOOwPgZfd4/uArS4jajWwDnhaVY8CvSJyrhvzCuDejD5+ZtSlwMMu7vEgcKGI1Dj314WuzRiFtt4YAI1Z3VP2QzgaXvZUiPJiz9LoM/dU3si2IrwoGCCRNkujEMkle+p84GPASyLyvGv7O+ByEdmI5y7aD3wcQFV3icg9wG68zKurXeYUwCeA7wGleFlTD7j2W4E7RaQZz8LY6sbqEJEvAs+4465X1Y7JvdWFQVuPE41hlkZZOMixvthsTKkgiMSTXvZU2I9pmMDmC889NbRgYdjcUwXLuKKhqo+TPbZw/xh9bgBuyNK+Ezg9S3sUuGyUsW4DbhtvnoZHW28UGCkaFeaeGhPfPVXh3FOWdps/4sl09hXh5p4qSGxF+DxjNPeUl3JrP4SjMRBPURYOUeYC4RbTyA/ptBJPpSkZVho9ZIv7ChYTjXlGW09syBWzT3mxxTRGQ1UH3VN+INxKieQHv1RISdGw/TSsjEjBYqIxz2jtjY5wTYEXCI8l01ZZNAuxZJq0Qmk4SFWJt3K5eyAxy7OaH0Rd9eDh+2mEQwGStrivIDHRmGe098RGuKaAwSto89WPxN+AqSwcpDQcpLQoSEdffJZnNT+IJb2LlOIiqz01XzDRmGe0jWppuEVr5nYZgR+/8OMZteVhOvpNNPJBLOHcU8NjGoEAybSSTpu1UWiYaMwz2npjNFaOtDTKbNHaqEQTQ/3udRVhjpto5IVo0rmnika6pwBbq1GAmGjMI/piSSLxFI1VIy2NClu0Niox/4ctZJZGvvEtjeLQ8EC4l4JrRQsLDxONeUS7S7dtqBgpGrZobXTiw/zuJhr5wxfkkhExDe9vS8woPEw05hHH3YrvuorwiOds0droDAZrncuktizM8X5bPZ8PTpzbkSm3gKXdFiAmGvMI3w9fn9XSsED4aIwQjYow0USaARPYKTNayq25pwoXE415xHGXJlpbPtLSOJFya6IxnFhiaEyjzp0/szamzugpty4QnjRLo9Aw0ZhHdLgfuTFFwyyNEYywNMo9S83iGlNnMKYxinsqadlTBYeJxjzieH+c8nBwRMkGgLIi3z1lLpfhxIf53WsHLQ0TjanipzOPZmnErWhhwWGiMY/o6I9TlyWeARAICGVh270vG8NdKL57ylaFT53hrj8ff/tXWxVeeJhozCM6+uNZXVM+ZeGQZU9l4cQ6De/rUOOLhlkaU8YX5OEpt6GAi2mYaBQcJhrziGN98cGr5GxU2D7hWfF/2PxVylUlIYqCYu6pPDB4brPsEQ6WcjsRkqk0L7Z0zfY0ctrudYWI/FpE9ojILhH5tGuvFZHtIrLX3ddk9LlWRJpF5FUR2ZLRfraIvOSeu8lt+4rbGvZu175DRFZl9NnmXmOviGzDGJWO/ti4loaVERmJv2rZ/2ETEWrKwnSaaEyZaCJFKCCEhomG756ySre5c+/zR/jwN57gUEdkVueRi6WRBD6rqqcC5wJXi8gG4BrgIVVdBzzk/sY9txU4DbgI+KaI+A7Nm4Gr8PYNX+eeB7gS6FTVtcDXgBvdWLXAdcA5wGbgukxxMk6gqmPGNADKbSOmrMRTI3/Yasut/lQ+iCXTWRMzBlNuzdLImd1HewDY1943q/MYVzRU9aiqPuce9wJ7gOXAxcDt7rDbgUvc44uBu1Q1pqpvAM3AZhFZClSp6lOqqsAdw/r4Y/0IuMBZIVuA7araoaqdwHZOCI2RQW8sSSKlY7qnym3L16zEEukRi8/qKsKDKczG5IklUyPOLVhMYzLsbfPEohAsjUGc2+gsYAewWFWPgicsQKM7bDlwKKNbi2tb7h4Pbx/SR1WTQDdQN8ZYxjA6xljY51Mett37shFLpikedjVcW15sgfA8EM0iyHDCPRU391TO7HOicbBQRENEKoAfA59R1Z6xDs3SpmO0T7ZP5tyuEpGdIrKzvb19jKnNX/zVy7VZ6k75mHsqO7FkakSgtq48zPG+OJ5RbEyW8dxTVrAwN/piSQ53DQAFIhoiUoQnGN9X1Z+45lbncsLdt7n2FmBFRvcm4Ihrb8rSPqSPiISAaqBjjLGGoKq3qOomVd3U0NCQy1uad/glROrLR49peCm3ZmkMx7M0hn4VmmpK6Y0lzdqYIrFEajArLROLaUwM38oIBwMc7BiY1bnkkj0lwK3AHlX9asZT9wF+NtM24N6M9q0uI2o1XsD7aefC6hWRc92YVwzr4491KfCwi3s8CFwoIjUuAH6hazOG4f+4jW9pJO3qeRjx5EgXyprGCgBeP9Y/G1OaN2Rz/UFmyq19FnOh2YnGOSfXcqgjMqvf4VwsjfOBjwHvE5Hn3e1DwJeAD4jIXuAD7m9UdRdwD7Ab+CVwtar6PpFPAN/FC47vAx5w7bcCdSLSDPwNLhNLVTuALwLPuNv1rs0Yhp/pM1YgvCwcIq1WWXQ4sWR6xIrlNfVONGY5U6XQiSayB8IHq9xawcKc2NvWR1FQeNe6evpiSboiiVmbS2i8A1T1cbLHFgAuGKXPDcANWdp3AqdnaY8Cl40y1m3AbePNc6HT0R+nbJS6Uz6+3z6WzO4yWKhky/BZXlNKOBRgX7tZGlMhlkxTVVo0ot3cUxOjua2X1fXlrHYXMwc7IoOVC2Ya++WYJxzvG3thH5yorRSzq7shxBIjYxrBgLC6rtwsjSkSy+L6g8wqt2b15kJzWx9rGytYWVsGzG4w3ERjntATTbKobOQVXSb+l9dEYyixZHpE9hTAmsZyszSmSGwc91TcPos50RlJ0FhZQlNNKWCiYeSB/lhycB/w0fD99vZFHUo8S0wD4OT6Cg52ROx8TYHRUm5FhFBAzD2VI5F4krJwkPLiEPUV4Vld4GeiMU/ojycH9wEfjROWhq3VyCSWTI1wT4FnaaTSysEOszYmy2grwsFzUZlojE88mSaR0sEtm5dUl9DWO3vVCkw05gn9sdTg7nyj4Qe//QJ9hsdofvc1DV7Q0VxUk8cr0ZI9OaMoKJbJlwN+kVHfk1BTFp7V9UMmGvOEvliSiuLRM6fghHvKYhpDyZZyC3ByQwUBgR/ubDEX1SSJJlMj9tLwCYfM0sgFfw+ccvf9rikL0xUx0TCmSCSXmEaRuaeyMdqq5YriEH/3oVP51Z5Wrv6P52xR5ARJpZVESsewNEw0ciESG25pFM2qpTHuOg1j7pNOK/3x8d1TvgvGrpqHEk9ld08B/OW7TqYzEudff72P9r4YjZUlMzy7wiWezL4/uE/I3FM5ERluaZSH6YkmSabSI/YpmQnM0pgHRNw+zOaemjjjXQ0DnL+2HoDX3rQ1GxMhmhi6je5wioIBu4DJgf4sMQ2AroHZWRVuojEP8Mud5xwIN/fUIONdDQO8ZXElAK+8OVZxZ2M4J/YHzy7IJaGgfRZzIOIqU5f7ouEW8c7WzpImGvOAQdEYd52GZU8Nx//RGu1qGKCuopj6ijCvtfbO1LTmBeOd25KiAFH7LI7LoKXhPAm1ztLonKX6UyYa8wB/j4ycYxoWfBzEvxoeyz0FsH5xJa+2mntqIox3bkvDQQYSZmmMx2BMw10U+pUfZisYbqIxD+gbdE+NE9NwbgKzNE7gu6fGK+D4liWV7G3tJW21knLGj2mMlnJbEgoOHmOMju9JKHWL+/wac7OVdmuiMQ/wP1S2Inzi5OKeAi+uEYmnaOmc3Q1wConxLI2SIhONXPAtDX9FuB8I7zDRMCbL8OyK0QgFBBHLnsrE96mPJxrrl3jB8FctrpEzvkU7WpJBscU0ciIS99YR+ZWBS8NBSooCFgg3Jo8f0xjP0hARikMBE40MBq+Gx9iHBGCd28nvVcugyplB99RoMQ2zNHIiEk9SHh56DmvKwnM3EC4it4lIm4i8nNH2eRE5PGwnP/+5a0WkWUReFZEtGe1ni8hL7rmb3JavuG1h73btO0RkVUafbSKy19387WCNYfTnGNMAz1VgufEnGEy5HcfSqCwpYll1CXvbLBieK7Fx0pnNPZUb/bHUCC9CTVl4Tlsa3wMuytL+NVXd6G73A4jIBmArcJrr800R8X/JbgauwtszfF3GmFcCnaq6FvgacKMbqxa4DjgH2Axc5/YJN4aRq3sKcJaGfVF9/HORy06GaxorBvdqNsYnp5TbZNrKs4xDJJ4ccUFYU15E51yNaajqY0Cu+3JfDNylqjFVfQNvL/DNIrIUqFLVp9T7hNwBXJLR53b3+EfABc4K2QJsV9UOVe0EtpNdvBY8/bEkpUVBgoHRduU9QXFRwLKnMojlaGkArG2sYF97n2VQ5Ugui/v8FfnG6PTHR7E05qp7agw+JSIvOveVbwEsBw5lHNPi2pa7x8Pbh/RR1STQDdSNMZYxjL4cyqL7hIMW08gk13Ua4IlGNJHmcJdlUOXCeGVE/BTSqFm+YxKJjbQ0asvDc9fSGIWbgTXARuAo8E+uPdulro7RPtk+QxCRq0Rkp4jsbG9vH2ve85L+HMqi+xRb6YYhxMb5Yctkrdtfo9n2Dc+J8QTZTz6wuMbY9MdTlBYNvShcVBameyBBchYW6k5KNFS1VVVTqpoGvoMXcwDPGliRcWgTcMS1N2VpH9JHREJANZ47bLSxss3nFlXdpKqbGhoaJvOWChrP55mbpVFcZJZGJv7q+LFqT/msdRlU+yyukRO+GIwWLyqxsjY5kS2mUVtWhCp0z0LRwkmJhotR+Pwh4GdW3QdsdRlRq/EC3k+r6lGgV0TOdfGKK4B7M/r4mVGXAg+7uMeDwIUiUuPcXxe6NmMYfbHkuHWnfCzldiiDawlycE/VVRRTU1ZkwfAcGYh7W72OFmvzYx1WSmRsItliGn7RwllwUY37SyMiPwDeC9SLSAteRtN7RWQjnrtoP/BxAFXdJSL3ALuBJHC1qvqfiE/gZWKVAg+4G8CtwJ0i0oxnYWx1Y3WIyBeBZ9xx16tqrgH5BUV/LEV9RTinY4tDwVkrqTwXmUggHDxro7mtjyeaj7G6vpxli0qnc3oFjfdjN7oYl5p7KicisZHrNOrKiwE43hdnbePMzmdc0VDVy7M03zrG8TcAN2Rp3wmcnqU9Clw2yli3AbeNN8eFTn8syUl1ZTkdGw4FBv34RkbKbY6b2axtrOCuZw7x0e/u4IJTGrn1z94+ndMraLJdIWdSMigaZvmORjqtRBIpyoa5n+srvYvE9r7YjM/JVoTPA/rjyXFXg/sUhwJW5TaDWDJNOBggkEO6MsDGFYtQhTUN5fxm7zF6oma1jcZAIjmYIZUNv5ChuadGJ5pMocoIS6OhwrM0jvWaaBiTINuK0dEoDgUt8JhBPDn6Vq/ZuOzsFfz22gv48qVnEk+leXhP2zTOrrAZzz1VYu6pcfFLBA23NGrKwgQDYpaGMXFU1VkaOabcWvbUEGLJVE6ZUz6BgLCkuoSzVixiSVUJ9790dBpnV9gMxFODcYtsmGiMTyTub7A29DwGAkJdeZhjvTMfCDfRKHAicWe+TsA9Zes0ThBLpHOOZ2QSCAgXnb6ER15rH6z9ZQxlIDGepWEpt+MxaGlk8SQ0VBabpWFMHL/uVM4rwi3ldgixZHrcCrejcd6aOuLJNPtssV9Wcg2EW0xjdCKDdeVGfkbrK4ppt5iGMVFObPWa+4rwuBWJGySaSE0oppHJylovY+1Qh5V8i5IAAAAgAElEQVQVycZAPDVq3SmwlNtcGNzqNcv3u6GymGNmaRgTZbAs+gQW94HtE+4Tiedet2s4K5xoHOyI5HNK84ZIPJljINw+i6MRGaOCtS8aM30BaKJR4Pgpn5UlRTkdf2LLV/uigufeG+uHbSwqikPUloc51GmikY3xsqeCAaEoKFawcAwGPQlZRKO+ophESme8lIiJRoHTM+BdiVSV5lp7yvsSW/DRIxJL5WylZWNFTSmHzNIYQSqtxJLpMddpgFcefSBuojEafswy23lsqPTWasx0XMNEo8DpdZZGVa6WRtC3NOyLCs7SyDEelI0VtWUmGlnwg9vjWXElYau6PBY9zorIdlHolw4y0TAmRE/UWRq5ioZLc7QtXz0i8SlaGrVlHO4aIGUbMw0hMniFPPa5LSkKWExjDLoHEpSFg1kLajb6lsYMB8NNNAoc39KoKJlYINxiGh7jBWvHY2VtGYmU8mZPNI+zKnx8l1PZOOnM5p4am65IgurS7BeEDRUlgFkaxgTpGfDqTuWy1SucKAFuouH53aOJdM4lWLKxosZlUB03F1UmObunioIWCB+DroHRRaOqNEQ4GOBY38yuCjfRKHB6owmqcrQyIMPSsNz4EyUaphTT8EqjWwbVUPz1BeMFwkuLgrZOYwy6IwkWlWUXDRGhviJsloYxMXqiiZzTbeHELmpmaZz4YZuKpbFsUSkBwYLhw/BdTmPVngIvxmYxjdHpHkiwqHT0vXIWV5fwyps9M7pWw0SjwOkZSOacbgsn3FMWCM9YGDkFS6MoGGDZolIOmHtqCLkKcolZGmPSNRAf1T0F8MebVrDrSA+/msFqy+OKhojcJiJtIvJyRlutiGwXkb3uvibjuWtFpFlEXhWRLRntZ4vIS+65m9y2r7itYe927TtEZFVGn23uNfaKiL8lrJFBb2xiloafPWWWRn4sDYBTllSx60h3PqY0b4iMsb4gExONsekawz0FcNnZTZxcX87/ffCVGcvgy8XS+B5w0bC2a4CHVHUd8JD7GxHZgLdd62muzzdFxP/U3Axchbdv+LqMMa8EOlV1LfA14EY3Vi3e1rLnAJuB6zLFyfDoGUhOLqZhwcdBS2Mq2VMAb22q5vVj/YOZbEZG9tS4MQ1zT41GNJEilkxTPYZohIIB/ueWt/Baax//8J+7Z8RNNa5oqOpjeHt3Z3IxcLt7fDtwSUb7XaoaU9U3gGZgs4gsBapU9Sn13tUdw/r4Y/0IuMBZIVuA7araoaqdwHZGiteCp3eCMQ3LnjpBJMcMn/E4o6kaVdh1pCcf05oXRHIUDcueGh2/PMhYMQ2AD56+hL84fzX/9sR+bnqoedrnNVm7fLGqHgVQ1aMi4m9tvhz4bcZxLa4t4R4Pb/f7HHJjJUWkG6jLbM/Sx8DbgKknOtGYhi3u84kMVgiemnvqjOXVALx8uJtzT66b8rzmA37KbS7uKVunkZ2uiCcaY8U0wMui+vvfP5XugQS7jnSTSmvOKfiTYWrflpFkm6mO0T7ZPkNfVOQqPNcXK1euHH+W84SBRIpUWieZPWVf1P4x9iqYCPUVxSyrLuHFFotr+ETiSYIBGXeDq5KiIDFXqt+FOQ1HV8RbfzFWTMMnEBBu/MgZANMqGDD57KlW53LC3fuh+xZgRcZxTcAR196UpX1IHxEJAdV47rDRxhqBqt6iqptUdVNDQ8Mk31LhMViscELuKdstzScywbLyY3FGUzUvHTbR8InEU5QVBccVghJLzBiVroHcLA2fUDBAaBK7UE6Uyb7CfYCfzbQNuDejfavLiFqNF/B+2rmyekXkXBevuGJYH3+sS4GHXdzjQeBCEalxAfALXZvh8MuiT8Q9FQoGCAbEvqRAv+93n0LKrc9bmxbxxrH+GS9TPVeJJlKU5GDBlbgYm7moRtLt3FO5WBozybi/NiLyA+C9QL2ItOBlNH0JuEdErgQOApcBqOouEbkH2A0kgatV1f80fAIvE6sUeMDdAG4F7hSRZjwLY6sbq0NEvgg84467XlWHB+QXNL0T3EvDpywcpM/2tc7ZhZILpy6tBKC5rY+zT7Ikv/H20vAZ3IjJ3KUjGAyEl40dCJ9pxhUNVb18lKcuGOX4G4AbsrTvBE7P0h7FiU6W524DbhtvjguVE+6piblXGipmZ0P6uUZ/zPthy4cvfUmVV06k1QoXAp5ojLcaHKA07Am2pd2OpGsgTjAglE8x5pZvbEV4ATPRXft86iuKOTYLG9LPNQamWBY9k6XVXsXRo90mGuCd25wsDXNPjUpXJMGi0qI5lyBgolHADO6lMYGYBkB9ZXhWNqSfa0x1A6ZMFpUVURwKmKXh8ErOj/+59IO8PbYwcgRdA4kxF/bNFiYaBcxEd+3zqSsv5nj/zJZTnotMdQOmTESEJdUlZmk4IvHUuGs0AGrKPX99h30eR9Az4Fkacw0TjQKmZyBJOBgYTKPNlfqKYroiCRKphe1H7o9NbQOm4SypKqHVRAPw1hDlcm5rTTRGZawNmGYTE40CxiuLHpqwz7O+0vuiHp/hzVvmGrlm+OTKkuoSjvYM5G28QibXc+unk3aaaIygayA+5zKnwESjoOmNJqmaxJVIfYW3t/BCj2t4MY38FUVYUl1Ca3dsRvc2mKsMxFOUFo1/botDQSqKQ3RETDSGY5aGkXd6BjxLY6L4orHQ024jsVRe0xmXVJUQT6UXvKtFVYnEk4PptONRWx42S2MYfbEkvdEki6tKZnsqIzDRKGA6I3FqJmG+NviWxgJPu801wydX/LTbNxd4BlVfLElac0/QqCkP0xGx7KlMDnd6bs6mmtJZnslITDQKmON9ceoqJi4afkxjpjekn0t4V8OpKe3aNxz/qvDNBR4MP9zl/eAtW5TbD15tWREd/Qv7AmY4h7u8nSCXm2gY+UJVOd4fG3Q1TYSycIiycHBBxzTiqTTJtObZ0vC+4Avd0vCvknP9waspD9PZb5ZGJi1maRj5JhJPEU2kB1MWJ0p9RfGCFg1/L418Zk/VV4QJiFkavqXRlKOlUVceXvBxoOEc7hwgHApQXz7xi8LpxkSjQPHTZesmLRoLe1W4v5dGvhb3gVdBuLGyxESjc4BwMJCzFVxTHmYgkbJSIhm0dA6wfFEpgWneG2MymGgUKMedD3gy7im/37HehXt1F8ljWfRMli0q4cDxSF7HLDRaugZYtqgk5x+8WpfM0Wlpt4O0dA3MSdcUmGgULL6lMVn3VN0Cd0/5KZ75zoM/c8UiXjzctaBX2x/uHJhQANdKiYzkcGeE5Tm692YaE40Cxbc0JpM9BdBQEaYjEieVXpgL0Y50TyzDJ1fevqqWaCLNriM9eR23kDjcNTChHzwrJTKUaCLFsb64WRpGfvELDtZNMlDWUFmMKhxfoNbGkS4v7rCsOr9fzE1uA6ad+xfmfmHRRIr23hjLF5Xl3McXDXNPebRMMPtsppmSaIjIfhF5SUSeF5Gdrq1WRLaLyF53X5Nx/LUi0iwir4rIloz2s904zSJyk9sSFrdt7N2ufYeIrJrKfOcTx/vilIeDOVUSzUajW1PQtkAX+B3uGqC2PDzp8zcajVUlrKwtY+f+zryOWyj4VX4n8oPnxzTM0vAYzD6ryV14Z5J8WBq/p6obVXWT+/sa4CFVXQc85P5GRDbgbeV6GnAR8E0R8b+xNwNX4e0pvs49D3Al0Kmqa4GvATfmYb7zguN9MWon6ZqCEwvRFur+D0ddsHY62HRSDTsPdCzIGlRH3A/eRNxTVaVFBMSKFvq0dLqFfQsopnExcLt7fDtwSUb7XaoaU9U3gGZgs4gsBapU9Sn1vmV3DOvjj/Uj4AKZa9tYzRLH++OTdk0BLK7y+rb2LExL40hXdHAxXr7ZtKqWY31x9i/ALKrJlL8IBoRFZWErWujYf6yfcCgwJ+tOwdRFQ4H/EpFnReQq17ZYVY8CuPtG174cOJTRt8W1LXePh7cP6aOqSaAbqJvinOcFx/vi1E/B0qivKEZk4VoaRyYYrJ0Ib22qBmDP0YUXDG/pjCDiVfydCLW2wG+QV97sZf3iCoJzcI0GwFRXNp2vqkdEpBHYLiKvjHFstjOgY7SP1WfowJ5gXQWwcuXKsWc8TzjeH+P05VWT7l8UDFBXXkxb78ITjZ5ogt5YctrcU2saKgDY19Y3LePPZV5o6WZtQwVFwYldjzZWFnO4a+F9FrOx52gPv/eWxvEPnCWmZGmo6hF33wb8FNgMtDqXE+6+zR3eAqzI6N4EHHHtTVnah/QRkRBQDYxIS1HVW1R1k6puamhomMpbKghUlY7+OHWTXNjns7iqeEG6p45MsKDeRCkNB1m+qJTm9oUlGslUmp37Ozjn5NoJ913XWMG+tr4FGQfKpL03xrG+OKcsnfwF4XQzadEQkXIRqfQfAxcCLwP3AdvcYduAe93j+4CtLiNqNV7A+2nnwuoVkXNdvOKKYX38sS4FHtaF/qkCeqJJEimddAkRn8VVJQvSPXXUT7edxkDjmsYK9i0w0dh1pIf+eIpzVk/cg7xucSV9sSRHFngJllfe9Fyapy6tnOWZjM5U3FOLgZ+6uHQI+A9V/aWIPAPcIyJXAgeBywBUdZeI3APsBpLA1arqF5v5BPA9oBR4wN0AbgXuFJFmPAtj6xTmO2/w11ZMdmGfz+KqYl5s6c7HlAqKwdLd0xQIB1jTUM4zb3SQTuucrB80Hex44zjApCyN9Yu9H8nXWnvnbNbQTODHwU5dMnctjUmLhqq+DpyZpf04cMEofW4AbsjSvhM4PUt7FCc6xgmODZYQmZp7qrGyhOP9MRKp9IR90IXMka4BQgGhoXL6KoiuaahgIJHizZ7otFo0c4HugQStPVF2vN7ByfXlNFZOPFa0frEXB9rb2jun/fnTzStHe1lSVTJYWmUukr8Sn8aMsf9YPwAn1U5t8c/iqhJUvb3Cpyv9dC5ypGuAJdUl05qdsrbRBcPb++a9aFz/8938+LkWRGDr21eM3yELi8rCNFQW81rrwnLpDWf30R5OmcOuKbAyIgXJa629FIcCrJiyaCzMtRrN7X2snOK5Gw8/g6p5nmdQqSqPvtbGKUsqeWvTIv7wrKbxO43C+sUV7G3tzePsCov+WJJ97X2cOoeD4GCiUZDsbetjTcPU87gX4qrw3miC3Ud62LRq4n73iVBfEaaqJDTvg+GvtvZyrC/Ole9czb1Xn8/m1ZM/r+saK9nb1kd6gRTRvP+lo3z+vl2DGWO/2tNKIqVz3j1nolGA7G3tHfQBT4VGZ2m0LSDRePZAJ2mFc6bw45YLIsKaxgpee3N+i8bje48BcP7a+imPtX5xJZF4arAC8XxGVfnq9tf43pP7+fffHgDg3uePsKy6ZLDo5VzFRKPA6I0mONIdZd3iqfs968qLCQZkQbmnnn6jg1BAOGvloml/rXNW1/HcwU56ovN3/+snmo9xckN5XuI2G5Z5bplnD8z/Yo+7jvTQ3NbHorIibrh/Dw/taeWx19r5g43L5ny2nYlGgbHX+cjX50E0ggFhSVUJhzoXTo2kp9/o4IymasryuM3raHxgQyPJtPLoq+3T/lqzQTyZZscbHbwzD1YGwFuXV9NQWcwvX34zL+PNZe59/jBFQeHuq86jurSIK2/fSTKtfPjMZbM9tXEx0Sgw/EBhPtxTAKcsqVwwNZKiiRQvtHRNye8+ETauqKGuPMyv9rTOyOvNNL/Z204knuI96/NThSEQEC46bQmPvNpOxO3hPh9JptLc98IR3rO+kbcsqeShz76XT1+wjivOO4kNczwIDiYaBcfe1j5KigJ5q7V/6tIq9rX3E02kxj+4wPndwS4SKWXzNAfBfYIB4X2nNPLrV9rm5favP9zZQn1FmHfnSTQAPnj6EgYSqXlrnakq//unL9PaExtMT64oDvH/fGA91198OoVQxNtEo4DY29rLY3vb85I55bNhWRWptM771FCAJ/cdIxiQGbM0AN6/YTE90SRP7js+Y685E3T0x3nolVYu2bg8rwtDN6+upaasiP986WjexpxLfOW/XuXunYf46/et5f0bFs/2dCaFiUaB8F+73uTCf36MQx0D/MX5q/M2rp8TvnsB7Gn9ePMx3tpUTWVJ0Yy95nvWN1BfUcwtj+2bsdecTlJp5f6XjvLFX+wmkVIu2zS5xXyjEQoG+MOzmnjg5Tfn3YXMz353mH/99T4u37yCv/nA+tmezqQx0SgAeqIJ/v5nL3PKkiqeuOZ9fOTsyS+gGs5JtWWUhYPsnudxjZ5oghdbuvMWtM2VkqIg//1dq3mi+Ti/O1j4WUFfemAPn/z+c/z0d4d559p63rIk/6uXP/l7aygJBfjKg6/mfezZ4le7W/lfP36Rc1bXFowbajRMNOYwqspT+47z2Xte4FhfjC/90RnU5rkmTSAgnLKksiBFQ1WJJXOLxex4vYNUWnnHmpkVDYCPnnsS1aVF3PxIYVsb/7HjIN/5zRt87NyT2PWFLdx55eZpeZ36imKuevcafrnrTZ7cd2xaXmMmuWfnIa66cyenLqnk5j89u+DrvBX27OcJvdEEP3muhS//8hXufubg4F7J//bEfi7/zm/59SttfPqC9Zy5YnrWFmxYVsWeoz0FtZeBqvL3P3uZs67fzo+ebRl17m09Uf7fn73Md3/zOiVFAd520vSvzxhORXGIj7ytiUdeay/YhIP9x/q5/he7eNe6eq77gw2UF4em9Wr5L9+1mpMbyvkfP/gdRwt4sV9nf5zrf76bc1bX8YOrzs37Rd9sYKIxy7zZHeXSm5/ib+55gW89uo/P/fgl3nnjwzzZfIyvP7yX89fW8fx1F/Lp96+btjmc2bSI3miSP/72UwWxsEpV+frDzXx/x0FqysL8zx++wNu+uJ0//vZTvOGKOfr8n/v3cOdvD7j1BA0Uh4KzMud3ra8nnkyzc//cP7+ZRBMpdh3p5nM/fpGiQID/e+mZhGbgSrm8OMQtHzubgXiKz9z1/LS/3nTxrcf20R9P8oWLT5uRtUEzwfx4FwVKTzTBZd9+ko6+OLdu28S71zew52gPn/z+c3z01h2owjUXnUpF8fT+m/7obU10DyS49fE3+IvvPcMv/vqdUy6GOF0cON7P3/7wRZ7e38ElG5fxlcvO5O6dh9h9pIcHXn6TS/71Cb71p2dz3po6Xmzp4mfPH+GT713Dn52/isrimQuAD2fzqlqKgsLjzcd457qZd5FNhmgixcXfeIJX3dqgGz9yxoT3/p4Kaxsr+dstb+HzP9/Nswc6OPukmct6myyReJIXDnVz9kk1HDjez+1P7ufiM5flZTHuXEEKySWRC5s2bdKdO3fO9jRy4nM/epEfPnuIez5+3pACentbe/nIzU/y3rc0ctPlZ83YfA4c7+e/ff1xasrCpNJKSVGAyzev5M/PXz1nNrm/4ran+d2BTq750ClsffvKIfM6eDzCX9z+DPuP9XPlu1azfVcr3QMJHvnb985oxtRo/PG3nmIgkeLnf/3OWZ1HJJ7k+UNdlBQFiSZSHOmK8vLhblJpZdmiUv7gzKU01ZTxD7/YzXcff4PrLz6Nt62s4fTl1bMy1/P+8WHOPbmWb39s04y//mik08r2Pa38+28PUFce5uKzlvN6ez+3PLaP1p4YK2pLOdYbp7w4yE8+cT4r6+bmRVgmIvKsqo57kgvC0hCRi4B/AYLAd1X1S9P1Wm8c6+fx5mMkU2nWNlZw+rJq0qpUlhQRDnlm+a4j3Ty+9xh/dv6qSbk7VJV7nz/C3TsP8VfvWTOi4uq6xZX85nPvozw8s66Uk+rK+eofb+RzP36RTSfV0NEf5x/+cw/lxSEu37wyr6+VTKXZdaSHipIQTTWlOZ3H11p7eey1dj77gfV89JyTRjy/sq6Mn3zyHVz9/ef49qOvs6ahnK/+ycY5IRjgFfX754deo7M/Pmub7Lze3sdVdz47Ip21LBwkHArQFUnw5QdfYXFlCa29UT527klccd6qWZmrN68QHzv3JP71kWZePtw9K8I1nLaeKH9zzws83nyM5YtKee5AJz97/ggAZzZV85n3r+f7Ow6woqaMr/3JxsFq0vOFOW9piEgQeA34ANACPANcrqq7sx0/WUujO5Lgq9tf5fs7DpLMUpo5ILC0upSGymJeaOlCFT6wYTHf/OjbKAoGON4X43cHuzjWF2PLaUuG/CjEk2meP9TFusYKOiNxvvDz3Tz6Wjtvbarmno+fR0nR7PjZx0NV+fA3nqA3muChz743q7WRTKV5en8HS6pKWF1fPmZwNJ1Wdh7o5Ml9x7jnmUOD+0HXlBXxR29rYuOKRayuL2d1fTnlWVxy1/z4RX76u8M8de0FYwYUk6k0+9r7Wb+4Yk6lNj57oIOP3PwUHzv3JD7z/nW8fKSHs1YuoipD1NJp5XDXAI1VxRSHghw8HuHWx1/nyX3HOXPFIt62soZTl1ayccWiCb23gXiK7/7mdW5+dB8lRUE+/+HTqCwOUVwUoLHS+98FA8LhrgF+uPMQR7ui1FWE+dT71s66L769N8YH/+UxeqJJPv7uk3nXugbOPqmGYEBo6YyQTCmNVcXTPs9oIsV/7DjIvzy0l1gyxd///ga2vn0F/fEUzx/qYm1jBcuqS+bUZ24i5GppFIJonAd8XlW3uL+vBVDVf8x2/GRF43hfjAu++ii/f8ZS/uo9aygNB3n5cDf72vsJBYTjfTEOdQ5wuHOAs1Yuoq4izP+5/xXOWF7N21fV8oOnDzLgMmMqikN8eOMyTq4vZ29rH7/a08rx/jihgCACJaEgn37/Ora9Y9WcT7974KWjfOL7z3HT5Wfx4TOXcaRrgCf3HefZAx30x1I8d7CTlk4vu+Xk+nL+dstbWN1QztGuKEe7o6RUWVRaxPrFldz4y1d4+JU2RLzS5JdvXkkqrWzf3cp/7W4llSHWjZXF1JSFiSSSRGIpIvEUA4kUl29ewT/+0Vtn63RMCVXluvt2ccdTBwbbTq4v50/PPYk7ntpPXyxJLJmmN5pkdX05F25YzL89sR9FefuqWnYd6aF7wKuYe/ZJNbxnfQMd/XHSWb7DARGv+mx1KUe6B7j5kX0c7Y6y5bTF/H9/cFrB7cPd1hvlf//0Zbbv9up4bVhaxWnLqvjRcy34b7+yOERjVTGLq0pYXFXC206q4d3r6mlu6+NQR4SjPVFePNRNbXmYD56xhD1He3izO0ZxUYDW7iiBgHD2STWUF4dIpxVV5eSGChKpND94+iC/2XuMWDLtMshOG9ydcb4wn0TjUuAiVf1L9/fHgHNU9VPZjp9KTKM3mpiQK+Mnz7XwjV8383p7Px88fQlXvnM1xaEg33p0H4++1k5fLEl1aRHvWFPHh85YysuHu0mklL9678mT2kd5NkillQ987VFeb++ntjxMh0sHri4tYlFZEcsXlXL55pV0DyS446n9Y27XWRQUrvngqVz6tiaqy4ae52gixf7j/bzR3s/rx/p5vb2f3miC8uIQZeEgZeEgFcVFXH7OioI5d6Pxq92tvHi4m9X1ZfzDL/ZwvD/OmU3VbFhWTTAAq+sr+N6Tb3CoY4APnbGE6/7gNBZXlQxaIY++1s7XH95La0+MyuIQoeDIK9tESumLnSj699amav73h07lnJPrZvKt5p2O/jiPvNrGVx58ldbeGFecdxKnLaumrTdKW0+Mtt4orT0xjnYNDFqyPqGAcOrSKlo6I3RGEgQDwuLKYmLJNA2VxQwkUhw4nr3ic0NlMb9/xlK2nLaEc0+uLVhrYizmk2hcBmwZJhqbVfWvM465CrgKYOXKlWcfOHAg61jTgarSGUmMcJeoKh39cWrLwwX/ATvcNcB9zx9hX3sfpyyp5B1r6jllSeWIuv/JVJrtu1tJqbK0upSl1SWEgkJbT4wXW7rZuGLR4J4Jhseb3VFebe3l3evqh3xOIvEkzW19vLUp+7qSVFpJpNKjujZVlTd7orT3xqgoDrGqrnzO79MwEaKJFAPx1KixIVVl15EenjvYySlLqji5oZyasjDBgBBNeO6kU5dWUV069OKlsz9OIpUmEBBU4dU3e4kmUrx7fcNgTHO+Mp9EY0bcU4ZhGAuZXEWjEKTzGWCdiKwWkTCwFbhvludkGIaxIJnzKbeqmhSRTwEP4qXc3qaqu2Z5WoZhGAuSOS8aAKp6P3D/bM/DMAxjoVMI7inDMAxjjmCiYRiGYeSMiYZhGIaRMyYahmEYRs6YaBiGYRg5M+cX900UEekFhm8uXA105/Fl5vp49UC+9smc6+813+P55OscFsL7ncufP5j753A+nL96oFxVG8YdTVXn1Q3YmaXtljy/xlwfb8Q5mENzm9Pj5fscFsL7ncufv0I4h/Ph/E3kNReKe+rnC2y8fDLX3+tcPndQGO/XzuHcGi/f5HV+89E9tVNzqJ8yn7FzMHXsHE4eO3dTYzbO30Recz5aGrfM9gTmAHYOpo6dw8lj525qzMb5y/k1552lYRiGYUwf89HSMAzDMKYJE40CQERWiMivRWSPiOwSkU+79loR2S4ie919jWuvc8f3icg3MsapFJHnM27HROSfZ+t9zST5OofuuctF5CUReVFEfiki9bPxnmaKPJ+7P3HnbZeIfHk23s9MM4nz9wERedZ9xp4VkfdljHW2a28WkZtkNnZ4y2dql92m5wYsBd7mHlcCrwEbgC8D17j2a4Ab3eNy4J3AXwHfGGPcZ4F3z/b7K6RziFcZug2od39/GW+TsFl/jwVw7uqAg0CD+/t24ILZfn9z8PydBSxzj08HDmeM9TRwHiDAA8AHZ/r9mKVRAKjqUVV9zj3uBfYAy4GL8b54uPtL3DH9qvo4EM0yHAAisg5oBH4zjVOfM+TxHIq7lburvCrgyPS/g9kjj+fuZOA1VW13f/8K+Mg0T3/WmcT5+52q+p+pXUCJiBSLyFKgSlWfUk9B7vD7zCQmGgWGiKzCuxLZASxW1aPgfTDxRCBXLgfudsmWMtoAAAOVSURBVB++BcVUzqGqJoBPAC/hicUG4NZpnO6cYoqfv2bgFBFZJSIhvB+8FdM327nHJM7fR4DfqWoMT2haMp5rcW0ziolGASEiFcCPgc+oas8Uh9sK/GDqsyospnoORaQITzTOApYBLwLX5nWSc5SpnjtV7cQ7d3fjWbj7gWQ+5ziXmej5E5HTgBuBj/tNWQ6b8Ys+E40Cwf1Y/Rj4vqr+xDW3OpMVd9+W41hnAiFVfXZaJjtHydM53AigqvuclXYP8I5pmvKcIV+fP1X9uaqeo6rn4dWI2ztdc55LTPT8iUgT8FPgClXd55pbgKaMYZuYBdeoiUYB4HzntwJ7VPWrGU/dB2xzj7cB9+Y45OUsMCsjj+fwMLBBRPzCbh/A81HPW/L5+RORRndfA3wS+G5+Zzv3mOj5E5FFwH8C16rqE/7BzoXVKyLnujGvIPfvfP6Y7cwCu41/w8tEUTxXyPPu9iG8bJSH8K7WHgJqM/rsBzqAPrwrlA0Zz70OnDLb76tQzyFeVtAeN9bPgbrZfn8FdO5+AOx2t62z/d7m4vkD/h7ozzj2eaDRPbcJeBnYB3wDt0B7Jm+2ItwwDMPIGXNPGYZhGDljomEYhmHkjImGYRiGkTMmGoZhGEbOmGgYhmEYOWOiYRgzjIj8lYhcMYHjV4nIy9M5J8PIldBsT8AwFhIiElLVb832PAxjsphoGMYEcUXnfolXdO4svFLXVwCnAl8FKoBjwJ+p6lEReQR4EjgfuE9EKoE+Vf2KiGwEvgWU4S3Y+gtV7RSRs4HbgAjw+My9O8MYG3NPGcbkeAtwi6q+FegBrga+Dlyqqv4P/g0Zxy9S1feo6j8NG+cO4HNunJeA61z7vwH/Q70aTYYxZzBLwzAmxyE9URfo34G/w9swZ7vbTC0IHM04/u7hA4hINZ6YPOqabgd+mKX9TuCD+X8LhjFxTDQMY3IMr7/TC+wawzLon8DYkmV8w5gTmHvKMCbHShHxBeJy4LdAg98mIkVuP4RRUdVuoFNE3uWaPgY8qqpdQLeIvNO1fzT/0zeMyWGWhmFMjj3ANhH5Nl6V0q8DDwI3OfdSCPhnvO06x2Ib8C0RKcOrPvznrv3PgdtEJOLGNYw5gVW5NYwJ4rKnfqGqp8/yVAxjxjH3lGEYhpEzZmkYhmEYOWOWhmEYhpEzJhqGYRhGzphoGIZhGDljomEYhmHkjImGYRiGkTMmGoZhGEbO/P+25kNxUd5IsAAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'][-200:].plot()\n" + ] + }, + { + "cell_type": "code", + "execution_count": 20, + "metadata": {}, + "outputs": [], + "source": [ + "first_sept_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n", + " for y in range(1985,\n", + " sorted_data.index[-1].year)]" + ] + }, + { + "cell_type": "code", + "execution_count": 22, + "metadata": {}, + "outputs": [], + "source": [ + "year = []\n", + "yearly_incidence = []\n", + "for week1, week2 in zip(first_sept_week[:-1],\n", + " first_sept_week[1:]):\n", + " one_year = sorted_data['inc'][week1:week2-1]\n", + " assert abs(len(one_year)-52) < 2\n", + " yearly_incidence.append(one_year.sum())\n", + " year.append(week2.year)\n", + "yearly_incidence = pd.Series(data=yearly_incidence, index=year)" + ] + }, + { + "cell_type": "code", + "execution_count": 24, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 24, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAD8CAYAAAC7IukgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X+w1fV95/HnCy+CTUCBgFUQoZV1oukWKwPu2N2N2AL5MUGnJmU1yozM4lptzU52VRJ2cJTulE4TWtdqJTURNQZdW0c3SsmNmonbJSAEG0VioZUikQr2omJnIF547x/f9ynnnl4O516+3HPOva/HzJn7PZ/v5/M5n/vlct/38+P7/SgiMDMzK9OwZjfAzMwGHwcXMzMrnYOLmZmVzsHFzMxK5+BiZmalc3AxM7PSObiYmVnpHFzMzKx0Di5mZla6jmY3YCB97GMfiylTpjS7GWZmbWXz5s3vRMT4vpQZUsFlypQpbNq0qdnNMDNrK5L+oa9lPCxmZmalc3AxM7PSObiYmVnpGgouknZKekXSy5I2Zdodkn6WaS9L+nRV/iWSdkh6XdLcqvSLs54dku6WpEwfIemxTN8gaUpVmYWStudrYVX61My7PcueeuKXw8zMytCXnstlETE9ImZUpa3MtOkR8SyApAuABcCFwDzgXkmnZP77gMXAtHzNy/RFwP6IOA9YCazIusYCy4BZwExgmaQxWWZFfv40YH/WYWZmLeBkDIvNB9ZExKGIeAPYAcyUdBYwOiLWR7FD2UPAFVVlVufxE8Dl2auZC3RGRFdE7Ac6gXl5bnbmJctW6mp7e98/yBfuX8/eAweb3RQzs35pNLgE8D1JmyUtrkq/WdJPJH2zqkcxEXizKs/uTJuYx7XpPcpERDfwHjCuTl3jgHczb21dPUhaLGmTpE379u1r8Nttrruf285LO7u4+/vbm90UM7N+afQ+l0sj4i1JE4BOST+lGOK6iyLw3AV8DbgeUC/lo046/ShTr66eiRGrgFUAM2bMaOk9nc9fupZD3Uf+5f0jG3bxyIZdjOgYxuvLP9XElpmZ9U1DPZeIeCu/7gWeBGZGxNsRcTgijgDfoJgTgaIXcU5V8UnAW5k+qZf0HmUkdQCnA1116noHOCPz1tbVtl689TI+N/1sRg4v/llGDh/G/Oln8+JtlzW5ZWZmfXPc4CLpI5JGVY6BOcCrOYdScSXwah4/DSzIFWBTKSbuN0bEHuCApEtyzuQ64KmqMpWVYFcBz+e8zDpgjqQxOew2B1iX517IvGTZSl1ta8LokYwa0cGh7iOM6BjGoe4jjBrRwYRRI5vdNDOzPmlkWOxM4MlcNdwBPBoRfyXpYUnTKYajdgI3AETEVkmPA68B3cBNEXE467oReBA4DVibL4AHgIcl7aDosSzIurok3QW8lPnujIiuPL4NWCNpObAl62h773xwiGtmncvVMyfz6MZd7POkvpm1IRWdgKFhxowZ4WeLmZn1jaTNNbehHJfv0Dczs9I5uJiZWekcXMzMrHQOLmZmVjoHFzMzK52Di5mZlc7BpU354ZZm1socXNqUH25pZq2s0QdXWovwwy3NrB2459Jm/HBLM2sHDi5txg+3NLN24GGxNuSHW5pZq/ODK83MrC4/uNLMzFqCg0sDfE+JmVnfOLg0wPeUmJn1TUPBRdJOSa9IelnSpkwbK6lT0vb8OqYq/xJJOyS9LmluVfrFWc8OSXfndsfklsiPZfoGSVOqyizMz9guaWFV+tTMuz3Lnnril6On85euZcrtz/DIhl1EFPeUTLn9Gc5fuvb4hc3MhrC+9Fwui4jpVZM6twPPRcQ04Ll8j6QLKLYpvhCYB9wr6ZQscx+wGJiWr3mZvgjYHxHnASuBFVnXWGAZMAuYCSyrCmIrgJX5+fuzjlL5nhIzs/45kWGx+cDqPF4NXFGVviYiDkXEG8AOYKaks4DREbE+iiVqD9WUqdT1BHB59mrmAp0R0RUR+4FOYF6em515az+/NL6nxMysfxq9zyWA70kK4P6IWAWcGRF7ACJij6QJmXci8KOqsrsz7cM8rk2vlHkz6+qW9B4wrjq9psw44N2I6O6lrlL5nhIzs75rNLhcGhFvZQDplPTTOnnVS1rUSe9PmXp19WyMtJhiKI7Jkyf3lqWu+689urR7+RWf6HN5M7OhqKFhsYh4K7/uBZ6kmP94O4e6yK97M/tu4Jyq4pOAtzJ9Ui/pPcpI6gBOB7rq1PUOcEbmra2rtu2rImJGRMwYP358I9+umZmdoOMGF0kfkTSqcgzMAV4FngYqq7cWAk/l8dPAglwBNpVi4n5jDqEdkHRJzplcV1OmUtdVwPM5L7MOmCNpTE7kzwHW5bkXMm/t55uZWZM1Mix2JvBkrhruAB6NiL+S9BLwuKRFwC7g8wARsVXS48BrQDdwU0QczrpuBB4ETgPW5gvgAeBhSTsoeiwLsq4uSXcBL2W+OyOiK49vA9ZIWg5syTrMzKwF+NliZmZWl58tZmZmLcHBxczMSufgYmZmpXNwMTOz0jm4mJlZ6RxczMysdA4uZmZWOgcXMzMrnYOLmZmVzsHFzMxK5+BiZmalc3AxM7PSObiYmVnpHFzMzKx0Di5mZlY6BxczMyudg4uZmZWu4eAi6RRJWyR9N9/fIelnkl7O16er8i6RtEPS65LmVqVfLOmVPHe3cu9kSSMkPZbpGyRNqSqzUNL2fC2sSp+aebdn2VNP7FKYmVlZ+tJzuQXYVpO2MiKm5+tZAEkXAAuAC4F5wL2STsn89wGLgWn5mpfpi4D9EXEesBJYkXWNBZYBs4CZwDJJY7LMivz8acD+rMPMzFpAQ8FF0iTgM8CfN5B9PrAmIg5FxBvADmCmpLOA0RGxPiICeAi4oqrM6jx+Arg8ezVzgc6I6IqI/UAnMC/Pzc68ZNlKXWZm1mSN9lz+GLgVOFKTfrOkn0j6ZlWPYiLwZlWe3Zk2MY9r03uUiYhu4D1gXJ26xgHvZt7aunqQtFjSJkmb9u3b1+C3a2ZmJ+K4wUXSZ4G9EbG55tR9wC8D04E9wNcqRXqpJuqk96dMvbp6JkasiogZETFj/PjxvWUxM7OSNdJzuRT4nKSdwBpgtqRHIuLtiDgcEUeAb1DMiUDRizinqvwk4K1Mn9RLeo8ykjqA04GuOnW9A5yReWvrMjOzJjtucImIJRExKSKmUEzUPx8RX8w5lIorgVfz+GlgQa4Am0oxcb8xIvYAByRdknMm1wFPVZWprAS7Kj8jgHXAHEljcthtDrAuz72QecmylbrMzKzJOo6f5Zj+UNJ0iuGoncANABGxVdLjwGtAN3BTRBzOMjcCDwKnAWvzBfAA8LCkHRQ9lgVZV5eku4CXMt+dEdGVx7cBayQtB7ZkHWZm1gJUdAKGhhkzZsSmTZua3Qwzs7YiaXNEzOhLGd+hb2ZmpXNwMTOz0jm4mJlZ6RxczMysdA4uZmZWOgcXMzMrnYOLmZmVzsHFzMxK5+BiZkPe3vcP8oX717P3wMFmN2XQcHAxsyHv7ue289LOLu7+/vZmN2XQOJFni5mZtbXzl67lUPfRbaoe2bCLRzbsYkTHMF5f/qkmtqz9uediZkPWi7dexuemn83I4cWvwpHDhzF/+tm8eNtlTW5Z+3NwMbMha8LokYwa0cGh7iOM6BjGoe4jjBrRwYRRI5vdtLbnYTEzG9Le+eAQ18w6l6tnTubRjbvY50n9UviR+2ZmVpcfuW9mZi2h4eAi6RRJWyR9N9+PldQpaXt+HVOVd4mkHZJelzS3Kv1iSa/kubtzu2NyS+THMn2DpClVZRbmZ2yXtLAqfWrm3Z5lTz2xS2FmZmXpS8/lFmBb1fvbgeciYhrwXL5H0gUU2xRfCMwD7pV0Spa5D1gMTMvXvExfBOyPiPOAlcCKrGsssAyYBcwEllUFsRXAyvz8/VmHmQ0Cvqmx/TUUXCRNAj4D/HlV8nxgdR6vBq6oSl8TEYci4g1gBzBT0lnA6IhYH8VEz0M1ZSp1PQFcnr2auUBnRHRFxH6gE5iX52Zn3trPN7M255sa21+jq8X+GLgVGFWVdmZE7AGIiD2SJmT6ROBHVfl2Z9qHeVybXinzZtbVLek9YFx1ek2ZccC7EdHdS11m1qZ8U+Pgcdyei6TPAnsjYnODdaqXtKiT3p8y9erq2RhpsaRNkjbt27evtyxm1iJ8U+Pg0ciw2KXA5yTtBNYAsyU9ArydQ13k172ZfzdwTlX5ScBbmT6pl/QeZSR1AKcDXXXqegc4I/PW1tVDRKyKiBkRMWP8+PENfLs21Hh8v3X4psbB47jBJSKWRMSkiJhCMVH/fER8EXgaqKzeWgg8lcdPAwtyBdhUion7jTmEdkDSJTlncl1NmUpdV+VnBLAOmCNpTE7kzwHW5bkXMm/t55v1icf3W0vlpsYnf+dSrpl1Lvs+ONTsJlk/9OkmSkmfBP5bRHxW0jjgcWAysAv4fER0Zb6vAtcD3cCXImJtps8AHgROA9YCvxsRIWkk8DBwEUWPZUFE/H2WuR74Sjbh9yPiW5n+SxQ9qbHAFuCLEVH3p9A3UVq12vH9Co/vm/XUn5sofYe+DVl73z/I8me38b2t/8jBD48wcvgw5l74i3z1Mx/3MIxZFd+hb9YHHt83O3kcXGxIG+zj+16sYM3iYTGzBux9/yA3f2cL91x9UVv1bJY++Qrf3riLa2ZOZvmVv9Ls5lib6s+wmB+5b9aA6hVl7fBL2jcjWrO552JWR7uuKPNiBSuTJ/TNStaud4x7sYI1m4fFzOpo51/S3mHRmsnBxew42vWX9P3XHh3FWH7FJ5rYEhuKPOdiZmZ1ec7FzMxagoOLmVmJfONqwcHF2pL/A1ur8lO2C57Qt7bUbjc12uDnG1d78oS+tZV2vanRBr/BfOOqJ/Rt0GvXmxpPBg8NtpZ2vifqZHBwsbbi/8BHtevY/mAOioP9Kdt94WExaznHewLxDQ9vYvyokT1uaqy+YXCwa/ehQT+puf2clJ0ocwviHwIjKBYAPBERyyTdAfxnYF9m/UpEPJtllgCLgMPA70XEuky/mKPbHD8L3JLbHI8AHgIuBv4J+O2I2JllFgJL8zOWR8TqTJ/K0W2OfwxcGxE/r/e9OLicXGU9lt6/fOpr17H9dg+KQ9nJmnM5BMyOiF8FpgPzJF2S51ZGxPR8VQLLBcAC4EJgHnCvpFMy/33AYmBavuZl+iJgf0ScB6wEVmRdY4FlwCxgJrBM0pgssyI/fxqwP+uwJjrRYZrzl65lyu3P8MiGXUQUq22m3P4M5y9dW3JL21u7Dg16vmxoOW5wicIH+XZ4vup1d+YDayLiUES8AewAZko6CxgdEeuj6C49BFxRVWZ1Hj8BXC5JwFygMyK6ImI/0EkR3ATMzrxk2UpdNsDKCgr+5dO4dhzbb9egaP3T0H0u2fPYDJwH/GlEbJD0KeBmSdcBm4AvZwCYCPyoqvjuTPswj2vTya9vAkREt6T3gHHV6TVlxgHvRkR3L3XZAHvx1suOOUzTF/7l07h2fShluz4E1PquoeASEYeB6ZLOAJ6U9AmKIa67KHoxdwFfA64H1FsVddLpR5l6dfUgaTHFUByTJ0/uLYudoDKDgn/5DG7tGhSt7/p0h35EvCvpB8C8iPijSrqkbwDfzbe7gXOqik0C3sr0Sb2kV5fZLakDOB3oyvRP1pT5AfAOcIakjuy9VNdV2+ZVwCooJvT78v1a48oKCv7lYzY4HDe4SBoPfJiB5TTgN4AVks6KiD2Z7Urg1Tx+GnhU0teBsykm7jdGxGFJB3IxwAbgOuB/VZVZCKwHrgKez1Vk64D/WTWJPwdYkudeyLxrsuxTJ3Ad7AQ5KJhZtUZ6LmcBq3PeZRjweER8V9LDkqZTDEftBG4AiIitkh4HXgO6gZtyWA3gRo4uRV6bL4AHgIcl7aDosSzIurok3QW8lPnujIiuPL4NWCNpObAl6zAzsxbgmyjNzKwuP1vMzMxagoOLmZmVzsHFzMxK5+BiZmalc3AxM7PSObiYmVnpHFzMzKx0Di5mZlY6BxczMyudg4uZmZXOwcXMzErn4GI2wPa+f5Av3L+evd6rxgYxBxezAXb3c9t5aWcXd39/e7ObYnbS9GmzMDPrv/OXruVQ95F/ef/Ihl08smEXIzqG8fryTzWxZWblc8/FbIC8eOtlfG762YwcXvy3Gzl8GPOnn82Lt13W5JaZlc/BxQa1VprfmDB6JKNGdHCo+wgjOoZxqPsIo0Z0MGHUyGY3zax0Di42qLXa/MY7Hxzimlnn8uTvXMo1s85l3weHmt0ks5PiuDtRShoJ/BAYQTFH80RELJM0FngMmEKxzfEXImJ/llkCLAIOA78XEesy/WKObnP8LHBLRISkEcBDwMXAPwG/HRE7s8xCYGk2Z3lErM70qcAaYCzwY+DaiPh5ve/FO1EOHbXzGxWe3zDru5O1E+UhYHZE/CowHZgn6RLgduC5iJgGPJfvkXQBsAC4EJgH3CvplKzrPmAxMC1f8zJ9EbA/Is4DVgIrsq6xwDJgFjATWCZpTJZZAazMz9+fdZgBnt9ollYahrTmOm5wicIH+XZ4vgKYD6zO9NXAFXk8H1gTEYci4g1gBzBT0lnA6IhYH0V36aGaMpW6ngAulyRgLtAZEV3ZK+qkCG4CZmfe2s838/xGk7TaMKQ1T0NLkbPnsRk4D/jTiNgg6cyI2AMQEXskTcjsE4EfVRXfnWkf5nFteqXMm1lXt6T3gHHV6TVlxgHvRkR3L3XVtn0xRW+JyZMnN/Lt2iBRmd+4euZkHt24i33+a/qk8TJrq9VQcImIw8B0SWcAT0r6RJ3s6q2KOun9KVOvrp6JEauAVVDMufSWxwan+689OkS8/Ip6P7J2ol689TKWP7uN7239Rw5+eISRw4cx98Jf5Kuf+Xizm2ZN0qfVYhHxLvADirmSt3Ooi/y6N7PtBs6pKjYJeCvTJ/WS3qOMpA7gdKCrTl3vAGdk3tq6zGyAeRjSah03uEganz0WJJ0G/AbwU+BpYGFmWwg8lcdPAwskjcgVXdOAjTmEdkDSJTlncl1NmUpdVwHP57zMOmCOpDE5kT8HWJfnXsi8tZ9vZk3gZdZWrZGlyP+WYsL8FIpg9HhE3ClpHPA4MBnYBXw+IrqyzFeB64Fu4EsRsTbTZ3B0KfJa4HdzKfJI4GHgIooey4KI+Psscz3wlWzO70fEtzL9lzi6FHkL8MWIqPvTfDKXIu99/yA3f2cL91x9kf9aM7NBpT9LkY8bXAaTkxlclj75Ct/euItrZk5m+ZW/clI+w8x6atc/6tqt3SfrPher4/yla5ly+zM8smEXEcUqmSm3P8P5S9c2u2lmg167Ln1u13b3hXsuJ2jv+wePuUqmHf4iMWtH7foEhnZtt3suTeBVMjZUtNLd9+36BIZ2bXd/OLiUwKtkbChopaGcdv2jrl3b3R/eLKwEjd6s126TeGbQunfft+sTGNq13X3lOZcB5BVl1o48r2j9mXNxz2UAtOpffmaNGEpDOVYez7kMgMEwiddKk7k28Dyv2Jpa+f+ley4DYDD85Vc9meshvaHHDwFtTa38/9LBZYC06ySeh/TMWk87/L/0hL7V5clcs9Yz0P8vfROllW4wDOmZDTbt8P/SwWWIa2RC0JO5Zq2n1f9felhsiPO9N2Z2PL7PxRrWDhOCZta+PCw2RA2Ge29aTSvfc2A20BrZ5vgcSS9I2iZpq6RbMv0OST+T9HK+Pl1VZomkHZJelzS3Kv1iSa/kubtzu2NyS+THMn2DpClVZRZK2p6vhVXpUzPv9ix7ajmXZGhohwnBdtNKD3Y0a7ZGhsW6gS9HxI8ljQI2S+rMcysj4o+qM0u6AFgAXAicDXxf0r+JiMPAfcBi4EfAs8A8iu2OFwH7I+I8SQuAFcBvSxoLLANmAJGf/XRE7M88KyNijaQ/yzru6/+lGHra9d6bVuMhRrN/7bjBJSL2AHvy+ICkbcDEOkXmA2tyP/s3JO0AZkraCYyOiPUAkh4CrqAILvOBO7L8E8A92auZC3RGRFeW6QTmSVoDzAauzjKrs7yDSx/4rutyvHjrZce858BsqOrTnEsOV10EbMikmyX9RNI3JY3JtInAm1XFdmfaxDyuTe9RJiK6gfeAcXXqGge8m3lr6zIbUB5iNPvXGg4ukj4K/AXwpYh4n6KX8MvAdIqezdcqWXspHnXS+1OmXl217V4saZOkTfv27esti9kJa/V7DswGWkNLkSUNpwgs346IvwSIiLerzn8D+G6+3Q2cU1V8EvBWpk/qJb26zG5JHcDpQFemf7KmzA+Ad4AzJHVk76W6rh4iYhWwCor7XBr5fs36ykOMZj01slpMwAPAtoj4elX6WVXZrgRezeOngQW5AmwqMA3YmHM3ByRdknVeBzxVVaayEuwq4Pko7u5cB8yRNCaH3eYA6/LcC5mXLFupy8zMmqyRnsulwLXAK5JezrSvAP9J0nSK4aidwA0AEbFV0uPAaxQrzW7KlWIANwIPAqdRTOSvzfQHgIdz8r+LYrUZEdEl6S7gpcx3Z2VyH7gNWCNpObAl6zAzsxbgx7+YmVldfiqymZm1BAcXM7MW1O6PE3JwMbO2/0U2GLX744T8VGQza+m92IeawfI4IU/omw1htb/IKtrtF9lgUvYWxnvfP8jN39nCPVdf1O+nRnhC38z6xFsvtJ6yHyfUrOE1D4uZtaAy/tpshJ+L1prKeGJ5s4fXHFzMWtBAzoF464XWU8bjhJr9tG4HF7MW0oy/Ngf7c9EGqhfYaprdK/Wci1kL8RxI+dp9Se+JaObTut1zMWshzf5rczBp9pxDK2hmr9Q9F7MW471hyuFeYHO552LWYgb7HMhAcS+wuRxczGzQ8kq45vEd+mZmVpfv0Dczs5bg4GJmZqU7bnCRdI6kFyRtk7RV0i2ZPlZSp6Tt+XVMVZklknZIel3S3Kr0iyW9kufulqRMHyHpsUzfIGlKVZmF+RnbJS2sSp+aebdn2VPLuSRmZnaiGum5dANfjoiPA5cAN0m6ALgdeC4ipgHP5Xvy3ALgQmAecK+kU7Ku+4DFwLR8zcv0RcD+iDgPWAmsyLrGAsuAWcBMYFlVEFsBrMzP3591mJlZCzhucImIPRHx4zw+AGwDJgLzgdWZbTVwRR7PB9ZExKGIeAPYAcyUdBYwOiLWR7GK4KGaMpW6ngAuz17NXKAzIroiYj/QCczLc7Mzb+3nm5lZk/VpziWHqy4CNgBnRsQeKAIQMCGzTQTerCq2O9Mm5nFteo8yEdENvAeMq1PXOODdzFtbl5mZNVnDwUXSR4G/AL4UEe/Xy9pLWtRJ70+ZenX1bIy0WNImSZv27dvXWxYzMytZQ8FF0nCKwPLtiPjLTH47h7rIr3szfTdwTlXxScBbmT6pl/QeZSR1AKcDXXXqegc4I/PW1tVDRKyKiBkRMWP8+PGNfLtmZnaCGlktJuABYFtEfL3q1NNAZfXWQuCpqvQFuQJsKsXE/cYcOjsg6ZKs87qaMpW6rgKez3mZdcAcSWNyIn8OsC7PvZB5az/fzMyarJHHv1wKXAu8IunlTPsK8AfA45IWAbuAzwNExFZJjwOvUaw0uykiDme5G4EHgdOAtfmCIng9LGkHRY9lQdbVJeku4KXMd2dEdOXxbcAaScuBLVmHmZm1AD/+xczM6vLjX6zl7X3/IF+4fz17/QBBs0HNwcUG1FDeFdBsKPEj921AeFdAs6HFPRcbEN4V0GxocXCxAeFdAc2GFg+L2YDxroBmQ4eXIpuZWV1eimxmZi3BwcXMzErn4DKI+YZFM2sWB5dBzDcsmlmzeLXYIOQbFs2s2dxzGYR8w6JZ+TzM3DcOLoOQb1g0K5+HmfvGw2KDlG9YNCuHh5n7xzdRmpnVsff9gyx/dhvf2/qPHPzwCCOHD2Puhb/IVz/z8SEzGnBSbqKU9E1JeyW9WpV2h6SfSXo5X5+uOrdE0g5Jr0uaW5V+saRX8tzdudUxuR3yY5m+QdKUqjILJW3P18Kq9KmZd3uWPbUv37SZWaM8zNw/jcy5PAjM6yV9ZURMz9ezAJIuoNii+MIsc6+kUzL/fcBiYFq+KnUuAvZHxHnASmBF1jUWWAbMAmYCyySNyTIr8vOnAfuzDjOzk6IyzPzk71zKNbPOZd8Hh5rdpJZ33DmXiPhhdW/iOOYDayLiEPCGpB3ATEk7gdERsR5A0kPAFcDaLHNHln8CuCd7NXOBzojoyjKdwDxJa4DZwNVZZnWWv6/BNpqZ9cn91x4dEVp+xSea2JL2cSKrxW6W9JMcNqv0KCYCb1bl2Z1pE/O4Nr1HmYjoBt4DxtWpaxzwbuatrcvMzFpAf4PLfcAvA9OBPcDXMl295I066f0pU6+uf0XSYkmbJG3at2/fsbKZmVmJ+hVcIuLtiDgcEUeAb1DMiUDRizinKusk4K1Mn9RLeo8ykjqA04GuOnW9A5yReWvr6q2tqyJiRkTMGD9+fF+/VTMz64d+BRdJZ1W9vRKorCR7GliQK8CmUkzcb4yIPcABSZfkfMp1wFNVZSorwa4Cno9iffQ6YI6kMTnsNgdYl+deyLxk2UpdZmbWAo47oS/pO8AngY9J2k2xguuTkqZTDEftBG4AiIitkh4HXgO6gZsi4nBWdSPFyrPTKCby12b6A8DDOfnfRbHajIjoknQX8FLmu7MyuQ/cBqyRtBzYknWYmVmL8E2UZmZWV39uohxSwUXSPuAfejn1MYq5nHbjdg8st3tgtWu7oX3bfqx2nxsRfZq0HlLB5VgkbeprVG4FbvfAcrsHVru2G9q37WW2209FNjOz0jm4mJlZ6RxcCqua3YB+crsHlts9sNq13dC+bS+t3Z5zMTOz0rnnYmZmpRuUweUYe9D8qqT1uafM/5E0OtOHS1qd6dskLakq84Pcl6ayb82EFmr3qZK+lel/I+mTVWV63TunDdo90Nf7HEkv5L/7Vkm3ZPpYSZ25X1Bn1YNZ+7xfURu0e8CueV/bLWlc5v9A0j01dQ30z3iZbW/la/6bkjbntd0saXZVXX275hEx6F7AfwB+DXi1Ku0l4D/m8fXAXXl8NcU2AQC/QPHEgSn5/gfAjBZt903At/J4ArAZGJZNCNQeAAADr0lEQVTvNwL/juIhn2uBT7VJuwf6ep8F/FoejwL+FrgA+EPg9ky/HViRxxcAfwOMAKYCfwecMtDXvOR2D9g170e7PwL8OvBfgHtq6hron/Ey297K1/wi4Ow8/gTws/5e80HZc4mIH1I8Sqba+cAP87gT+K1KduAjKh6EeRrwc+D9gWhnrT62+wLguSy3F3gXmKHiuW+jI2J9FD8Rlb1zWrrdJ7N9xxIReyLix3l8ANhGsX3DfIp9gsivlev3L/sVRcQbQGW/ogG95mW1+2S171j62u6I+OeI+L/Awep6mvQzXkrbB1o/2r0lIioPAt4KjFTxrMg+X/NBGVyO4VXgc3n8eY4+cfkJ4J8ptg7YBfxRHH2GGcC3suv6P0521/sYjtXuvwHmS+pQ8ZDQi/Ncvb1zBlJf213RlOutYkO8i4ANwJlRPGyV/FoZtujPfkUn1Qm2u2LAr3mD7T6Wpv6Mn2DbK9rhmv8WsCWKzR/7fM2HUnC5HrhJ0maK7uHPM30mcBg4m2LI4MuSfinPXRMRvwL8+3xdO7BNBo7d7m9S/ANvAv4Y+H8UDwvt0343J1Ff2w1Nut6SPgr8BfCliKjXay1lj6GylNBuaMI170O7j1lFL2kD8jNeQtuhDa65pAsptpO/oZLUS7a613zIBJeI+GlEzImIi4HvUIw7QzHn8lcR8WEO0/w1OUwTET/LrweAR2nOUEKv7Y6I7oj4rxExPSLmA2cA26m/d04rt7sp11vScIr/dN+OiL/M5LdzGKAyBLM30/uzX1Ert3vAr3kf230sTfkZL6ntLX/NJU0CngSui4jK78k+X/MhE1wqKzIkDQOWAn+Wp3YBs1X4CHAJ8NMctvlYlhkOfJaj+9Y0vd2SfiHbi6TfBLoj4rWov3dOy7a7Gdc7r88DwLaI+HrVqeo9hqr3C+rPfkUt2+6Bvub9aHevmvEzXlbbW/2aSzoDeAZYEhF/Xcncr2teb7a/XV8UfynvAT6kiLiLgFsoVkr8LfAHHL2B9KPA/6aYvHoN+O9xdLXHZuAnee5PyBU2LdLuKcDrFBN036d4ammlnhkUP7B/B9xTKdPK7W7S9f51iq79T4CX8/VpYBzFooPt+XVsVZmv5nV9narVMgN5zctq90Bf8362eyfFYpEP8mfrgib9jJfS9la/5hR/CP5zVd6XgQn9uea+Q9/MzEo3ZIbFzMxs4Di4mJlZ6RxczMysdA4uZmZWOgcXMzMrnYOLmZmVzsHFzMxK5+BiZmal+//q8gNev7GS5AAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.plot(style='*')" + ] + }, + { + "cell_type": "code", + "execution_count": 23, + "metadata": {}, + "outputs": [ + { + "ename": "SyntaxError", + "evalue": "invalid syntax (, line 1)", + "output_type": "error", + "traceback": [ + "\u001b[0;36m File \u001b[0;32m\"\"\u001b[0;36m, line \u001b[0;32m1\u001b[0m\n\u001b[0;31m \u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n" + ] + } + ], + "source": [ + "" + ] + }, + { + "cell_type": "code", + "execution_count": 25, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "2014 1601698\n", + "1991 1663610\n", + "1995 1828304\n", + "2012 2183912\n", + "2003 2234514\n", + "2019 2254363\n", + "2006 2297262\n", + "2017 2322818\n", + "2001 2540826\n", + "1992 2590314\n", + "1993 2699482\n", + "2018 2701716\n", + "1988 2759663\n", + "2007 2786458\n", + "2011 2852504\n", + "2016 2859019\n", + "1987 2867464\n", + "2008 2984311\n", + "1998 3047298\n", + "2002 3115484\n", + "1994 3514133\n", + "1996 3540251\n", + "2009 3558474\n", + "2004 3572810\n", + "1997 3624129\n", + "2015 3647492\n", + "2000 3808190\n", + "2005 3831409\n", + "1999 3914003\n", + "2010 3992174\n", + "2013 4176872\n", + "1986 5050543\n", + "1990 5214494\n", + "1989 5461328\n", + "dtype: int64" + ] + }, + "execution_count": 25, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "yearly_incidence.sort_values()" + ] + }, + { + "cell_type": "code", + "execution_count": 26, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 26, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEKCAYAAADpfBXhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGfNJREFUeJzt3X2UHFWZx/HvQxIgZCCggUYDMr4gghkE0ugCq84AumBQj8gRFFRc2UEFRY1nN8uCrO6yRjS6qPgSlRcVGJWXs0oU9QgDgohMABkwgi5klSBB3gIDWWDg2T9utfQO3aSru2rqcvl9zpkz/VJd9avb1U9X3brdbe6OiIikZaOqA4iISPFU3EVEEqTiLiKSIBV3EZEEqbiLiCRIxV1EJEEq7iIiCVJxFxFJkIq7iEiCZpY583nz5nl/f3/L+x566CHmzJlT5uK7FnM2iDtfzNkg7nzK1r2Y8+XNtnLlyrvdfeueF+zupf0tXLjQ27n00kvb3le1mLO5x50v5mzucedTtu7FnC9vNmDMC6i/6pYREUmQiruISIJU3EVEEqTiLiKSIBV3EZEE5SruZvYRM7vJzG40s3PNbNOygomISPc6Lu5mNh/4EFB39wXADOCwsoKJiEj38nbLzARmm9lMYDPgjuIjiYhIr8xz/IaqmR0HnAysB37q7oe3mGYYGAao1WoLR0ZGWs5rYmKCvr6+bjKXLuZsUHy+8TXrCptXbTasXd/59APz5xa27E402q7Idc6r3TrHvN3FnA3izpc329DQ0Ep3r/e63I6Lu5ltBZwPHArcD3wfOM/dv9PuMfV63cfGxlreNzo6yuDgYN680yLmbFB8vv4lKwqb1+KBSZaNd/6tFquXLips2Z1otF2R65xXu3WOebuLORvEnS9vNjMrpLjn6ZbZH7jN3f/i7o8BFwB79xpARESKl6e4/xH4GzPbzMwM2A9YVU4sERHpRcfF3d2vBs4DrgXGs8cuLymXiIj0INdX/rr7ScBJJWUREZGC6BOqIiIJUnEXEUmQiruISIJU3EVEEqTiLiKSIBV3EZEEqbiLiCRIxV1EJEEq7iIiCVJxFxFJkIq7iEiCVNxFRBKk4i4ikiAVdxGRBKm4i4gkSMVdRCRBKu4iIgnquLib2U5mdn3T3wNm9uEyw4mISHc6/pk9d78Z2A3AzGYAa4ALS8olIiI96LZbZj/gv939f4oMIyIixTB3z/8gs9OBa939Sy3uGwaGAWq12sKRkZGW85iYmKCvry/3sqdDzNmg+Hzja9YVNq/abFi7vvPpB+bPLWzZnWi0XZHrnFe7dY55u4s5G8SdL2+2oaGhle5e73W5uYu7mW0M3AG83N3XPt209Xrdx8bGWt43OjrK4OBgrmVPl5izQfH5+pesKGxeiwcmWTbecW8fq5cuKmzZnWi0XZHrnFe7dY55u4s5G8SdL282MyukuHfTLXMgYa/9aQu7iIhUp5vi/nbg3KKDiIhIcXIVdzPbDHgdcEE5cUREpAidd44C7v4w8NySsoiISEH0CVURkQSpuIuIJEjFXUQkQSruIiIJUnEXEUmQiruISIJU3EVEEqTiLiKSIBV3EZEEqbiLiCRIxV1EJEEq7iIiCVJxFxFJkIq7iEiCVNxFRBKk4i4ikiAVdxGRBOX9mb0tzew8M/udma0ys73KCiYiIt3L9TN7wKnAxe5+iJltDGxWQiYREelRx8XdzLYAXgMcCeDujwKPlhNLRER6Ye7e2YRmuwHLgd8CrwBWAse5+0NTphsGhgFqtdrCkZGRlvObmJigr6+v++QlijkbFJ9vfM26wuZVmw1r13c+/cD8uYUtuxONtitynfNqt84xb3cxZ4O48+XNNjQ0tNLd670uN09xrwO/AvZx96vN7FTgAXc/sd1j6vW6j42NtbxvdHSUwcHB/ImnQczZoPh8/UtWFDavxQOTLBvvvLdv9dJFhS27E422K3Kd82q3zjFvdzFng7jz5c1mZoUU9zwnVG8Hbnf3q7Pr5wF79BpARESK13Fxd/c7gT+Z2U7ZTfsRumhERCQyeUfLfBA4OxspcyvwnuIjiYhIr3IVd3e/Hui5L0hERMqlT6iKiCRIxV1EJEEq7iIiCVJxFxFJkIq7iEiCVNxFRBKk4i4ikiAVdxGRBKm4i4gkSMVdRCRBKu4iIglScRcRSZCKu4hIglTcRUQSpOIuIpIgFXcRkQSpuIuIJCjXLzGZ2WrgQeBxYLKIX+gWEZHi5f0NVYAhd7+78CQiIlIYdcuIiCTI3L3zic1uA+4DHPiauy9vMc0wMAxQq9UWjoyMtJzXxMQEfX193WQu1fiaddRmw9r107/sgflzO5qu6LYbX7OusHnlbbtO17kojbYrcp3zarfOsb4mIO5sEHe+vNmGhoZWFtHlnbe4P9/d7zCzbYCfAR9098vbTV+v131sbKzlfaOjowwODuaMW77+JStYPDDJsvFueqx6s3rpoo6mK7rt+pesKGxeeduu03UuSqPtilznvNqtc6yvCYg7G8SdL282MyukuOfqlnH3O7L/dwEXAq/sNYCIiBSv4+JuZnPMbPPGZeD1wI1lBRMRke7l6XuoAReaWeNx57j7xaWkEhGRnnRc3N39VuAVJWYREZGCaCikiEiCVNxFRBKk4i4ikiAVdxGRBKm4i4gkSMVdRCRBKu4iIglScRcRSZCKu4hIglTcRUQSpOIuIpIgFXcRkQSpuIuIJEjFXUQkQSruIiIJUnEXEUmQiruISIJyF3czm2Fm15nZRWUEEhGR3nWz534csKroICIiUpxcxd3MtgMWAd8oJ46IiBTB3L3zic3OAz4FbA58zN0PajHNMDAMUKvVFo6MjLSc18TEBH19fd1kLtX4mnXUZsPa9VUnaS/mfHmzDcyfW16YFhrb3fiaddO63Gbt1jnW1wTEnQ3izpc329DQ0Ep3r/e63JmdTmhmBwF3uftKMxtsN527LweWA9TrdR8cbD3p6Ogo7e6r0pFLVrB4YJJl4x03zbSLOV/ebKsPHywvTAuN7e7IJSumdbnN2q1zrK8JiDsbxJ2vqmx5umX2Ad5kZquBEWBfM/tOKalERKQnHRd3d/9nd9/O3fuBw4BL3P2I0pKJiEjXNM5dRCRBXXXcuvsoMFpoEhERKYz23EVEEqTiLiKSIBV3EZEEqbiLiCRIxV1EJEEq7iIiCVJxFxFJkIq7iEiCVNxFRBKk4i4ikiAVdxGRBKm4i4gkSMVdRCRBKu4iIglScRcRSZCKu4hIgjou7ma2qZn92sx+Y2Y3mdknygwmIiLdy/NLTI8A+7r7hJnNAq4wsx+7+69KyiYiIl3quLi7uwMT2dVZ2Z+XEUpERHpjoWZ3OLHZDGAl8BLgNHf/pxbTDAPDALVabeHIyEjLeU1MTNDX19d2WeNr1nWcq2i12bB2fWWL36CY8+XNNjB/bnlhWmhsd1VuX+3WeUOviSrFnA3izpc329DQ0Ep3r/e63FzF/a8PMtsSuBD4oLvf2G66er3uY2NjLe8bHR1lcHCw7TL6l6zInasoiwcmWTbe1W+HT4uY8+XNtnrpohLTPFVju6ty+2q3zht6TVQp5mwQd7682cyskOLe1WgZd78fGAUO6DWAiIgUL89oma2zPXbMbDawP/C7soKJiEj38hzbPw84K+t33wj4nrtfVE4sERHpRZ7RMjcAu5eYRURECqJPqIqIJEjFXUQkQSruIiIJUnEXEUmQiruISIJU3EVEEqTiLiKSIBV3EZEEqbiLiCRIxV1EJEEq7iIiCVJxFxFJkIq7iEiCVNxFRBKk4i4ikiAVdxGRBKm4i4gkKM9vqG5vZpea2Sozu8nMjiszmIiIdC/Pb6hOAovd/Voz2xxYaWY/c/fflpRNRES61PGeu7v/2d2vzS4/CKwC5pcVTEREumfunv9BZv3A5cACd39gyn3DwDBArVZbODIy0nIeExMT9PX1tV3G+Jp1uXMVpTYb1q6vbPEbFHO+vNkG5s8tL0wLje2uyu2rnbKf117aekOv16q1y1fV89zc1nnbbmhoaKW713vNkLu4m1kfcBlwsrtf8HTT1ut1Hxsba3nf6Ogog4ODbR/bv2RFrlxFWjwwybLxPD1W0yvmfHmzrV66qMQ0T9XY7qrcvtop+3ntpa039HqtWrt8VT3PzW2dt+3MrJDinmu0jJnNAs4Hzt5QYRcRkerkGS1jwDeBVe7+ufIiiYhIr/Lsue8DvBPY18yuz/7eUFIuERHpQccdfO5+BWAlZhERkYLoE6oiIglScRcRSZCKu4hIglTcRUQSpOIuIpIgFXcRkQSpuIuIJEjFXUQkQSruIiIJUnEXEUmQiruISIJU3EVEEqTiLiKSIBV3EZEEqbiLiCRIxV1EJEEq7iIiCcrzG6qnm9ldZnZjmYFERKR3efbczwQOKCmHiIgUqOPi7u6XA/eWmEVERApi7t75xGb9wEXuvuBpphkGhgFqtdrCkZGRltNNTEzQ19fXdlnja9Z1nKtotdmwdn1li9+gmPPFnA3izld2toH5c7t+7IZer09nOl7LsT2vzW2dt+2GhoZWunu91wyFF/dm9Xrdx8bGWt43OjrK4OBg28f2L1nRca6iLR6YZNn4zMqWvyEx54s5G8Sdr+xsq5cu6vqxG3q9Pp3peC3H9rw2t3XetjOzQoq7RsuIiCRIxV1EJEF5hkKeC1wF7GRmt5vZe8uLJSIivei4k8rd315mEBERKY66ZUREEqTiLiKSIBV3EZEEqbiLiCRIxV1EJEEq7iIiCVJxFxFJkIq7iEiCVNxFRBKk4i4ikiAVdxGRBKm4i4gkSMVdRCRBKu4iIglScRcRSZCKu4hIglTcRUQSlKu4m9kBZnazmf3BzJaUFUpERHqT5zdUZwCnAQcCuwBvN7NdygomIiLdy7Pn/krgD+5+q7s/CowAby4nloiI9MLcvbMJzQ4BDnD3o7Lr7wRe5e7HTpluGBjOru4E3NxmlvOAu7sJPQ1izgZx54s5G8SdT9m6F3O+vNl2cPete13ozBzTWovbnvLO4O7LgeUbnJnZmLvXcyx/2sScDeLOF3M2iDufsnUv5nxVZcvTLXM7sH3T9e2AO4qNIyIiRchT3K8BdjSzF5rZxsBhwA/KiSUiIr3ouFvG3SfN7FjgJ8AM4HR3v6mHZW+w66ZCMWeDuPPFnA3izqds3Ys5XyXZOj6hKiIizxz6hKqISIJU3EVEEqTiLiKSoGdkcTez+WY2v+ocrZjZi8zsI2a2b9VZWok5n7J1L+Z8MWeDuPP1ku0ZVdzNrN/MLgMuBj5jZq+uOlMzM/tb4GeE7955n5m9v+JI/0/M+ZStezHnizkbxJ2v52zuHvUfsGnT5YOBz2aX3w18HxjIrlsF2fYFXthYPvBx4Ijs+quAHwKDyqdsz5Z8MWeLPV/R2aLcczezLczsq2Z2C/BZM9shu+stwB+zyyPAH4CjGg+bxny7mNkNwL8CZ5jZvh5afBdgWwB3vxr4JfAe5VO21PPFnC32fGVli7K4AwcAmxJW7lHg42Y2m3CI8kYAd38EOA94dXb9ibLCmNl2ZrZF002HAue7+2sIbzLvMLMdgXMa+TIXAgvMbJNnaz5lSzNfzNlizzdd2Sor7hbMNLP3mtkvzOw4M3txdvdLgEfdfRL4PHAfcATwU+B5ZvacbLpbgD+Z2V4lZdzZzH4EXAF80swaX3H8v8Bm2eXvAXcCiwjvrM9tOtK4l/CtmK94tuVTtjTzxZwt9nzTna2y4p4ddrwWeBdwCrAJ8PXs7juBu7J3qD8RVujFhEb4LU9+pfAs4J7s9kKY2Zymq7sBt7t7P3AJ8Nns9nuBR8xsc3e/F/g98Pwsxy+Bj2bTbQw8Dqx+NuRTtjTzxZwt9nxVZpu24m5me5nZp83syOy6ATsDF7v7D939FGAHM9sbWEN4J9s5e/gqoC+77UvAG8zsjYQ3hhrwmx6zbWVmZ5rZNcBSM9s6yzcAXGlm5u4/AO43s0WEI4bNs/vJrm8DPEE40tjGzL4OnAtMuvtdqeZTtu7FnC/mbLHniyXbtBR3M3s58BXgQeBtZvbRbNnzgQezFQc4E3gHoVhPAntnt19LOJP8sLtfDiwBjgT2Af7N3Z9omkc3XpMt7w2EExXHA1sQviBt2+woA+CsLN+vs3U5EMDdr8rmMdPdVwFHAzcB/+Hu76F3MedTtjTzxZwt9nxxZGs3jKbbP8Ie9lGEQ5CZ2W2fA47LLteBLwCHAPsDP2l67PaEwxYIxfw6wq857Q78F/C8pmlzD1PKGvdo4DJC18687PbvAR/KLr8QWJrdvyehf2xG07r9JZvPfMIRxbHAGcCXgTk9tl20+ZRNz6va7pmTzb3goZBmthvhpOebgZOAE7K71hB+gxXCO9CVwFuBnwPbmtmuZjbLQ//6GjN7tbtfQviqzE8DFwDnuvufG8vyrHVyOgh4E/AJYC9CXz+EUTiNo4Q/Ab8ADnT3awjvvEPZMieAq4E93X0N8E5Ct9CdwAnu/lAXmWg66nhjjPky0bXdM6TdQG3XC7Vdl/L8zN5TmNkrgR2Bn7r7Xwh75be4+5FmtgdwspnVgVHg78xsM3d/2Mx+A7yNMIbzHOAfgC+Y2XpgHLgtW8RXgXPcfV2OTObubmZ7Eg55fgGs8DB08qXAre5+iZndRviU6+uBlcBbzGyeu99tZr8HHjKzFwBfBI4ws20Ivz51D+EwCncfA8a6bLs64QjnQeAzwF3Ai6rM90xouxjbTW2ntqu6nrSSa8/dgllm9i4zu47Q2b8l0Ci+jwOrs73wawmHGXsBD/Pk8B6AxwiHI9sS9s5vJPS3Xwbc7e63Q9g777KwvwY4nXC2eX/gU9kkTwC3mNlsd78ty7cr4Qm7gzDetLEeMwjtc36W8XBgIbDcexj/amZzzeyMbL63Aae6+11mthHhXb2SfGY2I2u71xIOC6Npu2y76zOzM4ms3bJ8s7O2GySy7c7MNjGzORG33RYRt12fmW1qZmcRYdttUCd9N8AcYO/s8pZZuC+0mO44wkdm52fXDyH0r+9A+LqAy7LbNyV0ycxreuzuwMad5JmyzM2A9/HkEcAs4MPAMdn9WwE3ZPM/lND/1Z/dd1C2LvOyy+PAXEJ//4+a8wAb5c02Jd+5hE+X9RG6mo5umqZxbuJY4N+nK1/2vB5F2OAWE076RNF2TdkuyLarrWNpt6bHbg6sIPwqGcBHYmi7pmw/Ar6WXT8FeF8MbUd4TbybUAPOj6ntmrJdAnw3uy2q7a7Tvw3uuZvZ8cCtwAozq7n7/YR+ojuyvvI32ZMfIrqKcFK08WGkKwknVh9297OA+8zs24QTpTcDf+1Tcvfr3P3RDeWZkm1b4CJgEPg24aTFwYSjhclsvvcRTsZ+iNAXtg1PDrG8nDDW/lF3vwj4JuFTr6cRzmQ/1pSvm3f+5nzfAt6f5bsF2MnMlmZ7VH9v4YNZFxOOZkrPZ2H87c8JG97XgdcTzoPsSdhbqqztpmRbThhFcDDhMw4vq7LdpphN+HzGi81sHmG7n5HNt7LtrinbxoTt7PmEro4FZvapKtvOzGYRzrsdAnzG3d+a3bV703yr2u6as53i7o0973Fgl6rbLrcO3skGCYca3wA+kt22J6Fo3Z6FPwdYlt13MvDJpsdfA+yeXd6EMDxozyLemQgb8Kuarh9J2EN5N/DrptufD9yRXT6G8BHfrbLH/xB4QdO084rI1ibfuwhnwV8KfDf7ezvwn4Tx/tOdb8umy/9IeDEdHknbNWf7GGE42YtiaLemeb6b0P96IvBewgm2a6puuynZTiAc+cyLpe0IR2OHT7ntUODqqtuuTbYXZMuvvO1yrUsHK9sYtnMoMJpdnkXYq5qbXd+BsNe+J+Fw8DzCHtePCe9Ym5QSPvR5Gfz1t2D34Mmun3sIY0ob0/6MrNASDqV+mk3zT6U17lPz7Q5c0dh4m6abRTjpvG92/eTpyJctawvC+Y61wCez6/cAtSrbbkq2O7PlziHr8quy3Zqez/cQutwOBs7Obru7yrZrk20ku615KHGV29xBhKPXZVmGj2c15F5gmyq3uynZLiV8mdd2MWx3udclx0o/l/Bhopdn12dOuf9M4JDGRkTogjiakgr7lGU3NuizeHI8/beBT2eXn0M48nhB05OzgKavE57GfMc035Zd3jZrv10ryvcBwnjc5YR+7l9mLziLoO0a2U4jDDl7SUTt9n1CN8xcwmCAEwiDA06serubku3nhJFne0TUdj8hHI1tT9gjPo6wg1j5djcl29mErwnYMZa26/SvUXQ6YmZfBh5w9yXZ9Y0I4zKPAV4OHOo5+82LYmbbEfq4Pujut1j4ErLhLNd84Hov5pNxveZ7v7vfmt22O6GbalGW7wNV5cvy7EZ4Q76K0I+4gLDXUmnbZdkWEF5wXySMvjqICtvNzPoI3R6bENrqZYQPqBxP2GvekYrarkW2HQnnK15HOAe2H6H9KtvmLBsWnV1+BeG1eiXhI/iVbndTsi0gfCL+VMI31Fa63eWRd5z7cuDU7MTDzoSNeB/Ck3J8VYU9szvZGHkzO4pwPuB4QnfS7zwMzaxSI98fs3y3ETaSScLe/HVVhsvcQzgpeIK7f8vMjgBuiiTb/YR+4xsJz+ssqm23ScLIiscIe+yPE14D48BHK267dtkeMbM3EYp+pdtco3hm7iechzrR3c+perubku1Bwg7sKuBfqH6761jePffDCCdPHyF8U9kl7n5zSdlyMbMrCSfcVhPGmH7C3W+oNFSTKfnuBJbE0HZmNpewJ/cOwvfnLwdOc/fHnvaB06BFtm+6+7JqU7WWfUCl0b99Z9V5mmXZDgHO8DASpXJmtgnhdxveSTi6/grwZQ9f812pFtmWu/vnq02VX8fF3cx2JYz3PI9w8qiwr9ntVXYkcRJhb/g7Hj7ZFo2Y85nZTEJXzCOEbDE9r9FmazCzGcATnmcvaZrEnA3AzI4mDLv9dmzPbczZOpVrz11ERJ4ZYv2ZPRER6YGKu4hIglTcRUQSpOIuIpIgFXcRkQSpuIuIJEjFXUQkQf8Hi7KSI4WeoTEAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.hist(xrot=20)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], "metadata": { "kernelspec": { "display_name": "Python 3", @@ -16,10 +2401,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.3" + "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 } -