diff --git a/module3/exo1/analyse-syndrome-grippal.ipynb b/module3/exo1/analyse-syndrome-grippal.ipynb index 59d72b5b58a3ae26346460dd39e62a39c55243d7..df83241ea87a6fff5eb8f248120aad5cf44a2455 100644 --- a/module3/exo1/analyse-syndrome-grippal.ipynb +++ b/module3/exo1/analyse-syndrome-grippal.ipynb @@ -9,7 +9,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -23,18 +23,10 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "Les données de l'incidence du syndrome grippal sont disponibles du site Web du [Réseau Sentinelles](http://www.sentiweb.fr/). Nous les récupérons sous forme d'un fichier en format CSV dont chaque ligne correspond à une semaine de la période demandée. Nous téléchargeons toujours le jeu de données complet, qui commence en 1984 et se termine avec une semaine récente." - ] - }, - { - "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, - "outputs": [], - "source": [ - "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-3.csv\"" + "Les données de l'incidence du syndrome grippal sont disponibles du site Web du [Réseau Sentinelles](http://www.sentiweb.fr/). Nous les récupérons sous forme d'un fichier en format CSV dont chaque ligne correspond à une semaine de la période demandée. Nous avons téléchargé le jeu de données complet (de 1984 à une semaine récente) le 2025/08/27 en utilisant l'URL suivante : \n", + "https://www.sentiweb.fr/datasets/all/inc-3-PAY.csv.\n", + "\n", + "Pour réaliser l'analyse, nous commençons donc par charger la version locale de ce fichier." ] }, { @@ -61,11 +53,981 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 2, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020253432355517312.029798.03526.044.0FRFrance
120253331693712133.021741.02518.032.0FRFrance
220253231990014303.025497.03022.038.0FRFrance
320253131847012625.024315.02819.037.0FRFrance
420253031916614283.024049.02922.036.0FRFrance
520252931867313815.023531.02821.035.0FRFrance
620252832328518131.028439.03527.043.0FRFrance
720252732145317129.025777.03226.038.0FRFrance
820252632194517422.026468.03326.040.0FRFrance
920252532332318546.028100.03528.042.0FRFrance
1020252432315418577.027731.03528.042.0FRFrance
1120252332439119307.029475.03628.044.0FRFrance
1220252231875514333.023177.02821.035.0FRFrance
1320252132376018671.028849.03527.043.0FRFrance
1420252032026515814.024716.03023.037.0FRFrance
1520251931626412394.020134.02418.030.0FRFrance
1620251831811513975.022255.02721.033.0FRFrance
1720251732215017291.027009.03326.040.0FRFrance
1820251632856422550.034578.04334.052.0FRFrance
1920251533572129592.041850.05344.062.0FRFrance
2020251433757931232.043926.05647.065.0FRFrance
2120251333967333686.045660.05950.068.0FRFrance
2220251235254345627.059459.07868.088.0FRFrance
2320251135946952154.066784.08978.0100.0FRFrance
2420251036033453048.067620.09079.0101.0FRFrance
2520250938453174994.094068.0126112.0140.0FRFrance
262025083136020124824.0147216.0203186.0220.0FRFrance
272025073208952195988.0221916.0312293.0331.0FRFrance
282025063273519258159.0288879.0408385.0431.0FRFrance
292025053334395318416.0350374.0499475.0523.0FRFrance
.................................
210019852132609619621.032571.04735.059.0FRFrance
210119852032789620885.034907.05138.064.0FRFrance
210219851934315432821.053487.07859.097.0FRFrance
210319851834055529935.051175.07455.093.0FRFrance
210419851733405324366.043740.06244.080.0FRFrance
210519851635036236451.064273.09166.0116.0FRFrance
210619851536388145538.082224.011683.0149.0FRFrance
21071985143134545114400.0154690.0244207.0281.0FRFrance
21081985133197206176080.0218332.0357319.0395.0FRFrance
21091985123245240223304.0267176.0445405.0485.0FRFrance
21101985113276205252399.0300011.0501458.0544.0FRFrance
21111985103353231326279.0380183.0640591.0689.0FRFrance
21121985093369895341109.0398681.0670618.0722.0FRFrance
21131985083389886359529.0420243.0707652.0762.0FRFrance
21141985073471852432599.0511105.0855784.0926.0FRFrance
21151985063565825518011.0613639.01026939.01113.0FRFrance
21161985053637302592795.0681809.011551074.01236.0FRFrance
21171985043424937390794.0459080.0770708.0832.0FRFrance
21181985033213901174689.0253113.0388317.0459.0FRFrance
211919850239758680949.0114223.0177147.0207.0FRFrance
212019850138548965918.0105060.0155120.0190.0FRFrance
212119845238483060602.0109058.0154110.0198.0FRFrance
2122198451310172680242.0123210.0185146.0224.0FRFrance
21231984503123680101401.0145959.0225184.0266.0FRFrance
2124198449310107381684.0120462.0184149.0219.0FRFrance
212519844837862060634.096606.0143110.0176.0FRFrance
212619844737202954274.089784.013199.0163.0FRFrance
212719844638733067686.0106974.0159123.0195.0FRFrance
21281984453135223101414.0169032.0246184.0308.0FRFrance
212919844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2130 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202534 3 23555 17312.0 29798.0 35 26.0 \n", + "1 202533 3 16937 12133.0 21741.0 25 18.0 \n", + "2 202532 3 19900 14303.0 25497.0 30 22.0 \n", + "3 202531 3 18470 12625.0 24315.0 28 19.0 \n", + "4 202530 3 19166 14283.0 24049.0 29 22.0 \n", + "5 202529 3 18673 13815.0 23531.0 28 21.0 \n", + "6 202528 3 23285 18131.0 28439.0 35 27.0 \n", + "7 202527 3 21453 17129.0 25777.0 32 26.0 \n", + "8 202526 3 21945 17422.0 26468.0 33 26.0 \n", + "9 202525 3 23323 18546.0 28100.0 35 28.0 \n", + "10 202524 3 23154 18577.0 27731.0 35 28.0 \n", + "11 202523 3 24391 19307.0 29475.0 36 28.0 \n", + "12 202522 3 18755 14333.0 23177.0 28 21.0 \n", + "13 202521 3 23760 18671.0 28849.0 35 27.0 \n", + "14 202520 3 20265 15814.0 24716.0 30 23.0 \n", + "15 202519 3 16264 12394.0 20134.0 24 18.0 \n", + "16 202518 3 18115 13975.0 22255.0 27 21.0 \n", + "17 202517 3 22150 17291.0 27009.0 33 26.0 \n", + "18 202516 3 28564 22550.0 34578.0 43 34.0 \n", + "19 202515 3 35721 29592.0 41850.0 53 44.0 \n", + "20 202514 3 37579 31232.0 43926.0 56 47.0 \n", + "21 202513 3 39673 33686.0 45660.0 59 50.0 \n", + "22 202512 3 52543 45627.0 59459.0 78 68.0 \n", + "23 202511 3 59469 52154.0 66784.0 89 78.0 \n", + "24 202510 3 60334 53048.0 67620.0 90 79.0 \n", + "25 202509 3 84531 74994.0 94068.0 126 112.0 \n", + "26 202508 3 136020 124824.0 147216.0 203 186.0 \n", + "27 202507 3 208952 195988.0 221916.0 312 293.0 \n", + "28 202506 3 273519 258159.0 288879.0 408 385.0 \n", + "29 202505 3 334395 318416.0 350374.0 499 475.0 \n", + "... ... ... ... ... ... ... ... \n", + "2100 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "2101 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "2102 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "2103 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "2104 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "2105 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2106 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2107 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2108 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2109 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2110 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2111 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2112 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2113 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2114 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2115 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2116 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2117 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2118 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2119 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2120 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2121 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2122 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2123 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2124 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2125 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2126 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2127 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2128 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2129 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 44.0 FR France \n", + "1 32.0 FR France \n", + "2 38.0 FR France \n", + "3 37.0 FR France \n", + "4 36.0 FR France \n", + "5 35.0 FR France \n", + "6 43.0 FR France \n", + "7 38.0 FR France \n", + "8 40.0 FR France \n", + "9 42.0 FR France \n", + "10 42.0 FR France \n", + "11 44.0 FR France \n", + "12 35.0 FR France \n", + "13 43.0 FR France \n", + "14 37.0 FR France \n", + "15 30.0 FR France \n", + "16 33.0 FR France \n", + "17 40.0 FR France \n", + "18 52.0 FR France \n", + "19 62.0 FR France \n", + "20 65.0 FR France \n", + "21 68.0 FR France \n", + "22 88.0 FR France \n", + "23 100.0 FR France \n", + "24 101.0 FR France \n", + "25 140.0 FR France \n", + "26 220.0 FR France \n", + "27 331.0 FR France \n", + "28 431.0 FR France \n", + "29 523.0 FR France \n", + "... ... ... ... \n", + "2100 59.0 FR France \n", + "2101 64.0 FR France \n", + "2102 97.0 FR France \n", + "2103 93.0 FR France \n", + "2104 80.0 FR France \n", + "2105 116.0 FR France \n", + "2106 149.0 FR France \n", + "2107 281.0 FR France \n", + "2108 395.0 FR France \n", + "2109 485.0 FR France \n", + "2110 544.0 FR France \n", + "2111 689.0 FR France \n", + "2112 722.0 FR France \n", + "2113 762.0 FR France \n", + "2114 926.0 FR France \n", + "2115 1113.0 FR France \n", + "2116 1236.0 FR France \n", + "2117 832.0 FR France \n", + "2118 459.0 FR France \n", + "2119 207.0 FR France \n", + "2120 190.0 FR France \n", + "2121 198.0 FR France \n", + "2122 224.0 FR France \n", + "2123 266.0 FR France \n", + "2124 219.0 FR France \n", + "2125 176.0 FR France \n", + "2126 163.0 FR France \n", + "2127 195.0 FR France \n", + "2128 308.0 FR France \n", + "2129 213.0 FR France \n", + "\n", + "[2130 rows x 10 columns]" + ] + }, + "execution_count": 2, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "raw_data = pd.read_csv(data_url, skiprows=1)\n", + "raw_data = pd.read_csv(\"inc-3-PAY.csv\", skiprows=1)\n", + "# Pour telecharger directement les donnees il faut indiquer l'URL donnée précédemment (ou une version actualisée si besoin)\n", + "# data_url = https://www.sentiweb.fr/datasets/all/inc-3-PAY.csv\n", + "# raw_data = pd.read_csv(data_url, skiprows=1)\n", "raw_data" ] }, @@ -78,9 +1040,73 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
18931989193-NaNNaN-NaNNaNFRFrance
\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low inc100_up \\\n", + "1893 198919 3 - NaN NaN - NaN NaN \n", + "\n", + " geo_insee geo_name \n", + "1893 FR France " + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] @@ -94,9 +1120,976 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020253432355517312.029798.03526.044.0FRFrance
120253331693712133.021741.02518.032.0FRFrance
220253231990014303.025497.03022.038.0FRFrance
320253131847012625.024315.02819.037.0FRFrance
420253031916614283.024049.02922.036.0FRFrance
520252931867313815.023531.02821.035.0FRFrance
620252832328518131.028439.03527.043.0FRFrance
720252732145317129.025777.03226.038.0FRFrance
820252632194517422.026468.03326.040.0FRFrance
920252532332318546.028100.03528.042.0FRFrance
1020252432315418577.027731.03528.042.0FRFrance
1120252332439119307.029475.03628.044.0FRFrance
1220252231875514333.023177.02821.035.0FRFrance
1320252132376018671.028849.03527.043.0FRFrance
1420252032026515814.024716.03023.037.0FRFrance
1520251931626412394.020134.02418.030.0FRFrance
1620251831811513975.022255.02721.033.0FRFrance
1720251732215017291.027009.03326.040.0FRFrance
1820251632856422550.034578.04334.052.0FRFrance
1920251533572129592.041850.05344.062.0FRFrance
2020251433757931232.043926.05647.065.0FRFrance
2120251333967333686.045660.05950.068.0FRFrance
2220251235254345627.059459.07868.088.0FRFrance
2320251135946952154.066784.08978.0100.0FRFrance
2420251036033453048.067620.09079.0101.0FRFrance
2520250938453174994.094068.0126112.0140.0FRFrance
262025083136020124824.0147216.0203186.0220.0FRFrance
272025073208952195988.0221916.0312293.0331.0FRFrance
282025063273519258159.0288879.0408385.0431.0FRFrance
292025053334395318416.0350374.0499475.0523.0FRFrance
.................................
210019852132609619621.032571.04735.059.0FRFrance
210119852032789620885.034907.05138.064.0FRFrance
210219851934315432821.053487.07859.097.0FRFrance
210319851834055529935.051175.07455.093.0FRFrance
210419851733405324366.043740.06244.080.0FRFrance
210519851635036236451.064273.09166.0116.0FRFrance
210619851536388145538.082224.011683.0149.0FRFrance
21071985143134545114400.0154690.0244207.0281.0FRFrance
21081985133197206176080.0218332.0357319.0395.0FRFrance
21091985123245240223304.0267176.0445405.0485.0FRFrance
21101985113276205252399.0300011.0501458.0544.0FRFrance
21111985103353231326279.0380183.0640591.0689.0FRFrance
21121985093369895341109.0398681.0670618.0722.0FRFrance
21131985083389886359529.0420243.0707652.0762.0FRFrance
21141985073471852432599.0511105.0855784.0926.0FRFrance
21151985063565825518011.0613639.01026939.01113.0FRFrance
21161985053637302592795.0681809.011551074.01236.0FRFrance
21171985043424937390794.0459080.0770708.0832.0FRFrance
21181985033213901174689.0253113.0388317.0459.0FRFrance
211919850239758680949.0114223.0177147.0207.0FRFrance
212019850138548965918.0105060.0155120.0190.0FRFrance
212119845238483060602.0109058.0154110.0198.0FRFrance
2122198451310172680242.0123210.0185146.0224.0FRFrance
21231984503123680101401.0145959.0225184.0266.0FRFrance
2124198449310107381684.0120462.0184149.0219.0FRFrance
212519844837862060634.096606.0143110.0176.0FRFrance
212619844737202954274.089784.013199.0163.0FRFrance
212719844638733067686.0106974.0159123.0195.0FRFrance
21281984453135223101414.0169032.0246184.0308.0FRFrance
212919844436842220056.0116788.012537.0213.0FRFrance
\n", + "

2129 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202534 3 23555 17312.0 29798.0 35 26.0 \n", + "1 202533 3 16937 12133.0 21741.0 25 18.0 \n", + "2 202532 3 19900 14303.0 25497.0 30 22.0 \n", + "3 202531 3 18470 12625.0 24315.0 28 19.0 \n", + "4 202530 3 19166 14283.0 24049.0 29 22.0 \n", + "5 202529 3 18673 13815.0 23531.0 28 21.0 \n", + "6 202528 3 23285 18131.0 28439.0 35 27.0 \n", + "7 202527 3 21453 17129.0 25777.0 32 26.0 \n", + "8 202526 3 21945 17422.0 26468.0 33 26.0 \n", + "9 202525 3 23323 18546.0 28100.0 35 28.0 \n", + "10 202524 3 23154 18577.0 27731.0 35 28.0 \n", + "11 202523 3 24391 19307.0 29475.0 36 28.0 \n", + "12 202522 3 18755 14333.0 23177.0 28 21.0 \n", + "13 202521 3 23760 18671.0 28849.0 35 27.0 \n", + "14 202520 3 20265 15814.0 24716.0 30 23.0 \n", + "15 202519 3 16264 12394.0 20134.0 24 18.0 \n", + "16 202518 3 18115 13975.0 22255.0 27 21.0 \n", + "17 202517 3 22150 17291.0 27009.0 33 26.0 \n", + "18 202516 3 28564 22550.0 34578.0 43 34.0 \n", + "19 202515 3 35721 29592.0 41850.0 53 44.0 \n", + "20 202514 3 37579 31232.0 43926.0 56 47.0 \n", + "21 202513 3 39673 33686.0 45660.0 59 50.0 \n", + "22 202512 3 52543 45627.0 59459.0 78 68.0 \n", + "23 202511 3 59469 52154.0 66784.0 89 78.0 \n", + "24 202510 3 60334 53048.0 67620.0 90 79.0 \n", + "25 202509 3 84531 74994.0 94068.0 126 112.0 \n", + "26 202508 3 136020 124824.0 147216.0 203 186.0 \n", + "27 202507 3 208952 195988.0 221916.0 312 293.0 \n", + "28 202506 3 273519 258159.0 288879.0 408 385.0 \n", + "29 202505 3 334395 318416.0 350374.0 499 475.0 \n", + "... ... ... ... ... ... ... ... \n", + "2100 198521 3 26096 19621.0 32571.0 47 35.0 \n", + "2101 198520 3 27896 20885.0 34907.0 51 38.0 \n", + "2102 198519 3 43154 32821.0 53487.0 78 59.0 \n", + "2103 198518 3 40555 29935.0 51175.0 74 55.0 \n", + "2104 198517 3 34053 24366.0 43740.0 62 44.0 \n", + "2105 198516 3 50362 36451.0 64273.0 91 66.0 \n", + "2106 198515 3 63881 45538.0 82224.0 116 83.0 \n", + "2107 198514 3 134545 114400.0 154690.0 244 207.0 \n", + "2108 198513 3 197206 176080.0 218332.0 357 319.0 \n", + "2109 198512 3 245240 223304.0 267176.0 445 405.0 \n", + "2110 198511 3 276205 252399.0 300011.0 501 458.0 \n", + "2111 198510 3 353231 326279.0 380183.0 640 591.0 \n", + "2112 198509 3 369895 341109.0 398681.0 670 618.0 \n", + "2113 198508 3 389886 359529.0 420243.0 707 652.0 \n", + "2114 198507 3 471852 432599.0 511105.0 855 784.0 \n", + "2115 198506 3 565825 518011.0 613639.0 1026 939.0 \n", + "2116 198505 3 637302 592795.0 681809.0 1155 1074.0 \n", + "2117 198504 3 424937 390794.0 459080.0 770 708.0 \n", + "2118 198503 3 213901 174689.0 253113.0 388 317.0 \n", + "2119 198502 3 97586 80949.0 114223.0 177 147.0 \n", + "2120 198501 3 85489 65918.0 105060.0 155 120.0 \n", + "2121 198452 3 84830 60602.0 109058.0 154 110.0 \n", + "2122 198451 3 101726 80242.0 123210.0 185 146.0 \n", + "2123 198450 3 123680 101401.0 145959.0 225 184.0 \n", + "2124 198449 3 101073 81684.0 120462.0 184 149.0 \n", + "2125 198448 3 78620 60634.0 96606.0 143 110.0 \n", + "2126 198447 3 72029 54274.0 89784.0 131 99.0 \n", + "2127 198446 3 87330 67686.0 106974.0 159 123.0 \n", + "2128 198445 3 135223 101414.0 169032.0 246 184.0 \n", + "2129 198444 3 68422 20056.0 116788.0 125 37.0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 44.0 FR France \n", + "1 32.0 FR France \n", + "2 38.0 FR France \n", + "3 37.0 FR France \n", + "4 36.0 FR France \n", + "5 35.0 FR France \n", + "6 43.0 FR France \n", + "7 38.0 FR France \n", + "8 40.0 FR France \n", + "9 42.0 FR France \n", + "10 42.0 FR France \n", + "11 44.0 FR France \n", + "12 35.0 FR France \n", + "13 43.0 FR France \n", + "14 37.0 FR France \n", + "15 30.0 FR France \n", + "16 33.0 FR France \n", + "17 40.0 FR France \n", + "18 52.0 FR France \n", + "19 62.0 FR France \n", + "20 65.0 FR France \n", + "21 68.0 FR France \n", + "22 88.0 FR France \n", + "23 100.0 FR France \n", + "24 101.0 FR France \n", + "25 140.0 FR France \n", + "26 220.0 FR France \n", + "27 331.0 FR France \n", + "28 431.0 FR France \n", + "29 523.0 FR France \n", + "... ... ... ... \n", + "2100 59.0 FR France \n", + "2101 64.0 FR France \n", + "2102 97.0 FR France \n", + "2103 93.0 FR France \n", + "2104 80.0 FR France \n", + "2105 116.0 FR France \n", + "2106 149.0 FR France \n", + "2107 281.0 FR France \n", + "2108 395.0 FR France \n", + "2109 485.0 FR France \n", + "2110 544.0 FR France \n", + "2111 689.0 FR France \n", + "2112 722.0 FR France \n", + "2113 762.0 FR France \n", + "2114 926.0 FR France \n", + "2115 1113.0 FR France \n", + "2116 1236.0 FR France \n", + "2117 832.0 FR France \n", + "2118 459.0 FR France \n", + "2119 207.0 FR France \n", + "2120 190.0 FR France \n", + "2121 198.0 FR France \n", + "2122 224.0 FR France \n", + "2123 266.0 FR France \n", + "2124 219.0 FR France \n", + "2125 176.0 FR France \n", + "2126 163.0 FR France \n", + "2127 195.0 FR France \n", + "2128 308.0 FR France \n", + "2129 213.0 FR France \n", + "\n", + "[2129 rows x 10 columns]" + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "data = raw_data.dropna().copy()\n", "data" @@ -122,7 +2115,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 5, "metadata": {}, "outputs": [], "source": [ @@ -152,10 +2145,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 6, + "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" @@ -179,9 +2170,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 7, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1989-05-01/1989-05-07 1989-05-15/1989-05-21\n" + ] + } + ], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", @@ -199,10 +2198,34 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 8, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXecHMWZ93/PzGzWKksgJIEEEjkYEDJgYxsTjf2C42t8DrKNjQN39vn83h2cAz4wB9zZxodtfMaHSLYJxvjgTBQimSyJpBxQzmGVVqvdnVDvH13VXd1TXd09YWdm9Xw/H312VN0Vuru6nnpCVZMQAgzDMAxTLqlaN4BhGIYZHLBAYRiGYSoCCxSGYRimIrBAYRiGYSoCCxSGYRimIrBAYRiGYSoCCxSGYRimIrBAYRiGYSoCCxSGYRimImRq3YCBZPTo0WLSpEm1bgbDMExDMW/evO1CiDFR5x1QAmXSpEmYO3durZvBMAzTUBDRmjjnscmLYRiGqQgsUBiGYZiKwAKFYRiGqQgsUBiGYZiKwAKFYRiGqQiRAoWIZhLRViJaoKWNJKJZRLRc/h2hHbuKiFYQ0VIiukBLP5WI5stjNxMRyfQWIrpPpr9KRJO0PDNkHcuJaIaWPlmeu1zmbS7/VjAMwzDlEEdDuQPAhYG0KwHMFkJMBTBb/h9EdCyASwEcJ/PcQkRpmefXAC4HMFX+U2VeBmCnEGIKgJsA3CjLGgngagDvBjAdwNWa4LoRwE2y/p2yDIZhGKaGRAoUIcTzALoCyZcAuFP+vhPAR7X0e4UQfUKIVQBWAJhOROMADBVCvCycbw7fFcijynoAwDlSe7kAwCwhRJcQYieAWQAulMc+KM8N1j/oeHrJFmzctb/WzWAYhomkVB/KQUKITQAg/46V6eMBrNPOWy/TxsvfwXRfHiFEDsBuAKMsZY0CsEueGyxr0PHlO+bi4l++UOtmMAzDRFJppzwZ0oQlvZQ8trKKG0R0ORHNJaK527ZtCzutLnGUOWB7d3+NW8IwDBNNqQJlizRjQf7dKtPXA5ionTcBwEaZPsGQ7stDRBkAw+CY2MLK2g5guDw3WFYRQohbhRDThBDTxoyJ3IqmrsgVQuUkwzBM3VGqQHkYgIq6mgHgIS39Uhm5NRmO8/01aRbbS0SnSx/IFwJ5VFmfBPC09LM8AeB8IhohnfHnA3hCHntGnhusf1CRy7NAYRimcYjcHJKI7gHwAQCjiWg9nMirGwDcT0SXAVgL4FMAIIRYSET3A1gEIAfgCiFEXhb1DTgRY20AHpP/AOA2AHcT0Qo4msmlsqwuIroWwBx53jVCCBUc8M8A7iWiHwN4Q5Yx6OjPF2rdBIZhmNhEChQhxGdCDp0Tcv51AK4zpM8FcLwhvRdSIBmOzQQw05C+Ek4o8aAmJwVKJmVyGzEMw9QXvFK+jlE+lDQLFIZhGgAWKHWMEigpYoHCMEz9wwKFYRiGqQgsUOoYtQ5FhC+zYRiGqRtYoDQAguUJwzANAAuUOkYJEpYnDMM0AixQGgGWKAzDNAAsUBoA9qEwDNMIsECpY1yTF8sThmEaABYoDQDLE4ZhGgEWKHUMm7oYhmkkWKA0AIJtXgzDNAAsUOoYDhtmGKaRYIHSALCCwjBMI8ACpY6pVznywZ88i8/f9mqtm8EwTJ0R+T0Uhgmycvs+rNy+r9bNYBimzmANpY5hZzzDMI0ECxSGYRimIrBAqWNYP2EYppFggcIwDMNUBBYodQy7UBiGaSRYoDAMwzAVgQVKXcMqCsMwjQMLFIZhGKYisECpY9iHwjBMI8EChWEYhqkILFDqGFZQGIZpJFigMAzDMBWBBUodwz4UhmEaCRYoDcyzS7fimB88ju6+XK2bwjAMU55AIaLvENFCIlpARPcQUSsRjSSiWUS0XP4doZ1/FRGtIKKlRHSBln4qEc2Xx24mIpLpLUR0n0x/lYgmaXlmyDqWE9GMcq6jXon6pvzPZi3D/mweK7Z2D1CLGIZhwilZoBDReADfAjBNCHE8gDSASwFcCWC2EGIqgNny/yCiY+Xx4wBcCOAWIkrL4n4N4HIAU+W/C2X6ZQB2CiGmALgJwI2yrJEArgbwbgDTAVytC64DBap1AxiGYTTKNXllALQRUQZAO4CNAC4BcKc8fieAj8rflwC4VwjRJ4RYBWAFgOlENA7AUCHEy8L5AMhdgTyqrAcAnCO1lwsAzBJCdAkhdgKYBU8IDRrYh8IwTCNRskARQmwA8BMAawFsArBbCPEkgIOEEJvkOZsAjJVZxgNYpxWxXqaNl7+D6b48QogcgN0ARlnKYhiGYWpEOSavEXA0iMkADgHQQUSfs2UxpAlLeql5gu28nIjmEtHcbdu2WZpXf7CGwjBMI1GOyetcAKuEENuEEFkADwI4E8AWacaC/LtVnr8ewEQt/wQ4JrL18ncw3ZdHmtWGAeiylFWEEOJWIcQ0IcS0MWPGlHip9Q1/KphhmHqgHIGyFsDpRNQu/RrnAFgM4GEAKupqBoCH5O+HAVwqI7cmw3G+vybNYnuJ6HRZzhcCeVRZnwTwtPSzPAHgfCIaITWl82XaoCIqygvEbnmGYeqHTKkZhRCvEtEDAF4HkAPwBoBbAQwBcD8RXQZH6HxKnr+QiO4HsEief4UQIi+L+waAOwC0AXhM/gOA2wDcTUQr4Ggml8qyuojoWgBz5HnXCCG6Sr2WRof1E4Zh6oGSBQoACCGuhhO+q9MHR1sxnX8dgOsM6XMBHG9I74UUSIZjMwHMTNjkhiLKksX6CcMw9QSvlGcYhmEqAguUQQD75BmGqQdYoDQw7JNnGKaeYIFSx7DmwTBMI8EChRlQ7p+zDv/26OJaN4NhmCrAAqWOiVyHop3ZKPzTn97Grc+vrHUzGIapAixQ6hgOG2YYppFggTIIYF8LwzD1AAuUOiZKThCHeTEMU0ewQGEYhmEqAguUOibuLsJs8WIYph5ggdLAsMGLYZh6ggVKHcOaB8MwjQQLFIZhGKYisECpYzgcmGGYRoIFyiCABQ/DMPUAC5S6xi4p1DIU/qY8wzD1AAuUBoY4zothmDqCBUodw4oHwzCNBAuUQQDLHYZh6gEWKHVMpKBgixfDMHUECxSGYRimIrBAqWPYh8IwTCPBAoVhGIapCCxQ6pjYuw2zJsMwTB3AAqWBUT75+N+er3827tpf6yYwDFMiLFASMHvxFry9fteA1Rf9xUZ53iCRJ6+s3IEzb3ga//PGhlo3hWGYEmCBkoDL7pyLi3/5Yq2bUcRgEShLNu0BALyxdmeNW8IwTCmwQKlj4gqKwWLyGhxXwTAHLixQBgGDRUNREPGKTYZpRFig1DFRmofaHHKQyROGYRqUsgQKEQ0nogeIaAkRLSaiM4hoJBHNIqLl8u8I7fyriGgFES0logu09FOJaL48djPJKSoRtRDRfTL9VSKapOWZIetYTkQzyrmORmewbF8/SC6DYQ5YytVQ/hPA40KIowGcBGAxgCsBzBZCTAUwW/4fRHQsgEsBHAfgQgC3EFFalvNrAJcDmCr/XSjTLwOwUwgxBcBNAG6UZY0EcDWAdwOYDuBqXXANGmL7UBiGYWpPyQKFiIYCeB+A2wBACNEvhNgF4BIAd8rT7gTwUfn7EgD3CiH6hBCrAKwAMJ2IxgEYKoR4WThT7bsCeVRZDwA4R2ovFwCYJYToEkLsBDALnhA68GCJwjBMHVCOhnI4gG0AbieiN4jov4moA8BBQohNACD/jpXnjwewTsu/XqaNl7+D6b48QogcgN0ARlnKGlTEXocySCTK4LgKhjlwKUegZACcAuDXQoiTAeyDNG+FYArdEZb0UvP4KyW6nIjmEtHcbdu2WZrXeCiBUijUth2VhoO8GKYxKUegrAewXgjxqvz/A3AEzBZpxoL8u1U7f6KWfwKAjTJ9giHdl4eIMgCGAeiylFWEEOJWIcQ0IcS0MWPGlHCZtSP+OhSGYZjaU7JAEUJsBrCOiI6SSecAWATgYQAq6moGgIfk74cBXCojtybDcb6/Js1ie4nodOkf+UIgjyrrkwCeln6WJwCcT0QjpDP+fJl2QDJ4orwGx3UwzIFKpsz8fwfg90TUDGAlgC/BEVL3E9FlANYC+BQACCEWEtH9cIRODsAVQoi8LOcbAO4A0AbgMfkPcBz+dxPRCjiayaWyrC4iuhbAHHneNUKIrjKvpe6I6xsZbMMw8acoGaYhKUugCCHeBDDNcOickPOvA3CdIX0ugOMN6b2QAslwbCaAmUnaO9hwFzYONonCMExDwivl65j4goIlSqNyx4urMOnKR9CbzUefzDB1DguUQQBrKI3LL59ZAQDY05utcUsYpnxYoNQxceXEYJMngyVsOJsv4MmFmyOCDQbJxTIMWKDUNVFRT4PtA1uDjV/MXo7L756H55YNrvVPDBMGC5RBQGGQSJRBchku63c6nzPe0d0fffIAX/urK3dg574Y7WKYBLBAqWMOVJNXNXh7/S786OGFA7vWxd0ax3JKDSxeuXwBn771FXzx9tcGvnJmUMMCZRAwWBYEVnNw/fRvXsEdL61GT//ARVN5Yd319Xz6cs5ePUu37K1xS5jBBguUeqa+xqGqU81xN5N2BvdcfuBuKsXRUBB9TqXplwKlKd1Yr/+kKx/BtX9ZVOtmMBYaq0cxPtSncutsAlw25SgqG3ftRzZfvFtmsxw8+w3Hqk6dPR91D1oyjff63/bCqlo3gbHQeD3qACL+1it1NmLViF09/Tjzhqdxzf8Wz2JdDWUAt2b2tI/o5zOQk4K+bGNqKEz9wz1qEDBYNJRyBWOXjFr66/LiMN1Myunq2dzA3SzvGwzhOpf7CYIBfIhZKVRTg2XBD1M3sECpY2JvXz9IBIqi1HFOmXKaDaacJqmhDKTJS3gSJfrcqrYkUNcg6y9M/cACpYGphUO3nlHah0mgZJQPJVeaQHlmyVbkEgojSvCACoUBFSkDWFdlqLdIOcYMC5Q6JuodqoW5pJqUexn9eSck2OQbUGkmh30Uzy3bhi/dMQe/eHpFonypGJ9oLmfH6HVdPXhpxfbE+RqxuzRimw9Eyv0eClMPDLKXjUq0ean1Fc1GgeKUWYpA6drXBwBYs2NfonxKWNiUj3ImBWf9+zMAgNU3fDhx3kZjsEyaBjusoVSIQkHgp08uxba9fRUrM/5KeX7ZAM+cZTR5SXWhFJOXcl4nXcIi4wBiza5rYfCqhRnpvJ89h/vmrE2cj3t4Y8ACpULMWd2FXzy9Av/0wFsDXjdP3hyycsQ3aSiZMtahKI0p+Sw5fr6BnIGrqmrRbZZv7cY//2l+4nysoTQGLFAqRF52+Epu7RF3BjlYXrVyr8OmoXgmr+S1uL6QhINaopXyAylQGrDHsDxpDFigVAjXuTqgdToMtpet1NURyj9icspTAm0hSFppKAmVm1SCBzSgQV4SU7P29maxdkdPleor/SIHWx8frLBAKQHTi1HqLNZaT+zzGu9tq8aMPC9H5XQqXCSVUi25PpSEGkosp/zAb5/jmbyKK73kly/iff/xTFXrLSlvA/bxAxEWKCVgGiBqMTB4tv2Bq7OaVOre2TScUgSZElBJ8yaZZNTCh2Ji5fZkkWyJ6i0j72Dp44MdFiglkDf07lQV1oTELqoB7QHV+CpunO+OlDIwpUrMm0TgD6hAkXdqoLtNOdfICxsbAxYoJWB6MeI4YCtNI6+Ur0ab3UHHIpBKGdTcsOGEEiXJGpNamLwGmnLqZQ2lMWCBUgLmAaIaJi97Ye6A1YBvm2nGWa6dPM5mjCXVUOJkIUkgQC0G+YGusqzn23hd/ICEBUoJmGaqaoCshWo+2N41m0CwEsv5PXB3K4mZrRbrLAa6SnbKD35YoJSAKXzUFTIV3BK8HncbrtSAXI0mV2sBYalPVOWz1VmJ/diSPhM2eTHVggVKCZjCR1XawH5hovS1FbXG1ORyL8Mm093BfQA/2JjkcsoZMJPet1rN9tkpP/hhgVICZpNX5euJKrIW30caiPe61OsqxBDqtRiWbPfMu9bSW1b6QB2erxoDOIcND35YoJSA6WVTQsaypq5qDKjJy1dvGTPOKgzt8cKG69O0VM6AmTSvu7DRkq8qE6Q66y9M5WGBUgJWk9cA+lDi2OgrjT4olOVkreKAZXsEpQxqlXymYZQTqZf0+YvAXxNJdwWIw0Ca9ZjaULZAIaI0Eb1BRH+R/x9JRLOIaLn8O0I79yoiWkFES4noAi39VCKaL4/dTPINJqIWIrpPpr9KRJO0PDNkHcuJaEa515EEk8lLDQi1+Er3QL5r9fxee88g/CmUMqiVOrOOM6uOsz1LZD2JNZTqBC9EV1xG1nrueIxLJTSUbwNYrP3/SgCzhRBTAcyW/wcRHQvgUgDHAbgQwC1ElJZ5fg3gcgBT5b8LZfplAHYKIaYAuAnAjbKskQCuBvBuANMBXK0Lrmqgv4S2KK9KTmajBiR3MWWtonaqVG6pt9C6Z5b8W8q9Kvf+2gZwb0HswGso1nOqIk9q4SdiBpKyBAoRTQDwYQD/rSVfAuBO+ftOAB/V0u8VQvQJIVYBWAFgOhGNAzBUCPGycN68uwJ5VFkPADhHai8XAJglhOgSQuwEMAueEKoKen82dW7PITxwOko5O+iWil5Vve0ea/ehlH6vSh0Ik1RVXkhtqe0Lz1eNPlWWFla5ZjBVpFwN5ecA/gmAPmc/SAixCQDk37EyfTyAddp562XaePk7mO7LI4TIAdgNYJSlrAHB6ENRd6CSGkrMt+g/nliaeFuQUtEH13JqNA3S5UYWVcuHUm6ocZwqyxnAS3bKW86pRn8q5/k24m4QByIlCxQi+giArUKIeXGzGNKEJb3UPP5KiS4norlENHfbtm2xGmpCL9zUuWuzDsVjwYbdA1KPX0OpTh1lhw3bBEoJ5VZzKPMCK0ovI/lAHceHUlpbBrrMaiOEwK+ffQdb9vTWuikNQTkaynsAXExEqwHcC+CDRPQ7AFukGQvy71Z5/noAE7X8EwBslOkTDOm+PESUATAMQJelrCKEELcKIaYJIaaNGTOmtCuF/6U1aSiFqvhQ7NRiHYpOOTbxagijOD6UUma65Zp/4vksaqChWMOGqxHW3Xg+lKVb9uLGx5fg7/7wRk3qbzRKFihCiKuEEBOEEJPgONufFkJ8DsDDAFTU1QwAD8nfDwO4VEZuTYbjfH9NmsX2EtHp0j/yhUAeVdYnZR0CwBMAzieiEdIZf75Mqxp+DaX4eD5GhFGlaeSFjaZiynd+q1/FN6ac7evjDMDmfDGivCrwHZ1q+FCqYkJtwCivnPxkdHdfrjYNaDAyVSjzBgD3E9FlANYC+BQACCEWEtH9ABYByAG4QgihPsD+DQB3AGgD8Jj8BwC3AbibiFbA0UwulWV1EdG1AObI864RQnRV4Vpcopzy+RjmluR1xn+Lkta7afd+LNm0F2cfPTb65BDKW4diiXwqUShXay+vco1e1fehmPNu2LUfhwxrLVpHIwJ/zWWW3JyqlMlRXo1BRQSKEOJZAM/K3zsAnBNy3nUArjOkzwVwvCG9F1IgGY7NBDCz1DYnRVfXretQamyGisvHb3kJm3b3YvUNH06Ur55XLMdxypdCqQNhkmyVXoeydPNeXPDz5/H9Dx+Dr5x1eOT5xWXWl8mrfnsdo8Mr5WOiv1/2zSEHTqKU885v2u04GZP6FCrllDeavEovzslfJU2gbFOc5cq89TGV1VBW73A+5fvqqmLFPd7CxpKbU5Uya705ZKNMFGsNC5QS6MsWO1Gq4pQfgHeoP58sJlb4flfWKa/SSo/ykvljnJOs3FJ9FEnqKKmK0LyutmY4XxT9KKYaW6/U27qleq63UWGBUgL7s8UOuoFaB6Kjd/akmpHaxDKbVKBUaC8vY9llfus8TthwaQsby6PaJiaTlhlHOFt9KFVZh1JG3so1oyRYQ4kHC5SY6C9DT3++6HjefYEruDlkxGtUjoaQSTuPvj9XuoZSFhYNpeQiRbRAKm3rFVluwquPdX6JnxfWMQlJlZIy9Ecvai281ups3ll6XnbKNwYsUGKiDw69BpOXzcRQLcp5x5qkipLU5OWrv/TqrSvlS97qRP61Ta5L0QRKvc9xVqRXYsdo0/XatLU497caA3i9rVsarNzy7Apcdsec6BOrQDXChgclvrBh00r5GvhQynnHUlKgJDXVVXMvL2E5FgdvEAwvoKR1KGUKuDhUK6TWaAqNUVd1fCil52UNJT7//vjSmtXNGkpM9O6cq8OtV5JS8u67wvizIsSZ0dtQj8Xm8C9lYCp1Ly+RQEKWJ5zDfSi2DmmrsRpRVeV93KyCDUlSb829N40FC5SYxN96ZeA+sFWJlyyxhqJvDllhJ2ulnPIV96GU1pxYOSuzUt6UFh3lZauz2gvlq/n1y0riBjc0zFTRoxah1ixQYqI/mrzB76CSBrbbld5hlMkr8bc0fKNCydUbO7urYZS5XbzRP2OpN4o4gsrenvh1lIItr80pb6Pauw0nvpc10hQa2dRmsqRUGxYoMfEvbCw+XhWbc9Rx7YSkilEldrmtNDFcIBH5wwd+NTCUdL1lOuVtVOI5mExycZzytkG6Kk553Q+ZsPxa9dN6ej+SUoulDCxQSsDklFdppT7CvlweO/f1J8pTie6S2PTg+11ZY5A30Nl56M0N+OmTxY5H9VjMJiD1t3QNJSlJTHiV1lA8U42hXTGqqkrYsPY7+Q7JtRnZa71CvxxYoNQz2rOxOeVLHRi+etc8nHztLH+VVezMyhSSVLOq1MJGm+M86rq/fe+b+MXTK4rSXdOUNSQ5OQl864E6S6ikBGz30ubTs/tQKt/39DIbRUOpxaBcKaphNYmCBUpM9EHK+AngQhkmFQDPL0v+8a9yBI67nXvCCKYKuVDsg365JiarySt54SW3x/0bXoC3rX452l64z8jqlLeUGTaQPr9sG06+5klsLeGDU/6Q88S5E9dXCdztfBrPJ1+Tr1yyQImJ/gLkDE4U9QJW6jO2QAwfSlk1OZTjlK+0BlV+2HC4FuJpP6WXm5Qk2aq1l5dxGUqMhoW1557X1mJnTxZvr0/+hdBqayg//ssifOXOuUmbFVFv42oo7JSvY/RHY9ttOKwDrtjajc27o2d1STpBeX29tHDVcr4p7xvILH6OcjUUcwSZfD4lvGSl3uaB8qEY16HIv8Yor6If8coEvP5Zijkl4vHHzhvGf7+wCk8t3pKo3LfW7cLFv3wB+w3bKQH28Ot6hzWUOkZ/wWxO+TAT0rk/ew6nXz87sh6f9lPF/qDGmXLsrOX4FIxaRIzoIxs2DcULSU5OyZpYjGxqfUN55kNDmm0gjFFZmMkrLTtOuYNVcg2lOi/DtX9ZhLfX78aCjWaNK19msE0tYR9KHaM/mmo45RVJ9tYqp6Y4e0i9tqqr2FZeRqWRJrwyNRRblJco4/mU70OJU0eFNRSZZNJQ3HNK2KImnXbKs2nSTyzcjE279xvK1LTbpL67Ko2Nni/RXEEDW7yMpvlqwwIlJlEx9EoOlNsB9d1/I3cbroBT3lbG//3Ny7jo5hf8dfp+J/W/aAOKze5fIp5T32TyCq83Clv0WJz22Iga0OLVU5xmcybHEaphbVcaii366Wt3z8PHb3nJ2s7kvrvqDI5RGmJDm7xYQ2kMTJLfi/IaWFOAIvnCRrVS3n7e9u4+3/+FX6IkIkoYuTP6Mp3gNhNQ2P196Z3t+Nrdc+37YiVtT4L8lXbK2/qRt5FpeKcJa4/6js5jCzaZ88mMmwz+Qr3IxAIl0dkJcCdW5sONHDZcC6c87zYck8iw4QqZvPQOHFVUOVW5PpRy9vJKWGdU2KhKs9l+9/RmQ48VLBIlSkO57I652J/No6c/j44W/2tR6m1O8nzK8qFY+qNJaMSZdYc9A3VvZi0yO79tz84f5WWp3EAwurBSe+a5m6SGPAG3nQ0YN8xO+XpGezY5g+ddvUjlPsKBUlOjfChxtIRyIsRsx21uJFtYqG0vsCgNxTUBmtrlmtLC22XC86GUbmJKUo+O0qJN46C6vzb/Sth9GtvZCgA4/9iDjcdtE5RyQs59/peEt6prXz/ueHGVsc6U7cGjNmajSlELpzxrKDHRH41pwCt3YaNXjrlOc5vK8aHYw4bDrsMfqZXUDq7nDa/TNrN6Y+1OS/nyGRj3tvL/DWITsKWHMcfxoZS2SaeOTWM2iQw10KQsk+6wtnsTp7AZvU2glKGhBOpIJfBq/MP9b+LZpdtwxhGjcdTBnb5j3sJSc95G9qGwU76O8TkUjVFe6rwyTV4J8uunlrq9dtgAYNLCgMrZsm2+Cts9iDOrNg12UR/fIks4bOl7eVmrDNRhTp+1aAu27e0zH7SUr+znpvuVl8/Wei9DtET13MK0SKuGopdfjoaSKCewcts+AM5+eUE8BcVu8mpEPYWd8nWM3uFMzi41EJXqxDNtwVFLH0pYellbkEeerzSM8BMzlmm1KzIM2VVa2ECpSjU925JvcxIfiuHc/lwBX71rLj7z21eseU23y/YFUfdTC5Y5SJhQL7gCxXwjbVv5lBPlhTLydvflfH91UlGaOjvlE8ECJSbRYcPlCRS37CQr5csyeckyQooI64xRZiudj93yIu56ebV2vl0YqcHIqqFYBIptNhk3aML0/MrebdhyTsoV7MUjsap3xdbuWPX48iqBYjg/zsfgQk1e7jMy54vrlE96S8vJm5WNtr2bYc+4Fiavrn39+PgtL2LDruK1PElgp3wdoz8a28LGUrUG1WH1F7KaHxXywoZDBo4Y9tco894ba3fhhw8t1M6PKM91ypemodg+hOVqKCFF2zQ2WziyjTh9IS2vJ2sKRY/ZmUyn5SxCI44PJdyHZtci4zrly9lDLipvsF/arAeuLzGkrFqEDf/5jQ14fe0u/Pb5lWWVw9vXNwjG2WQZexz5y/F+V9Pk5dZXwkyz1Pr10+2DfnjBaZvJy+IniVqcSJYFe1GCc09vFr+YvbworyuIrI5q568xcjDmgGBzypvqVuXafChhdatyw3xsdoHiHUs61iXRUIJtsO1i4XanUAEau4kVQ4jo5xOGfu0sUOoY/WUw7jYcMXOLwhTtE1WSfjxp30uV7ENJVo8/rz2zqtLmccjSAAAgAElEQVT2IsT5vodx65XAOWHYNJQwfvrEUvx01jI8Ot+/2C9WuLA8x6yhRGaXZRSjrsO26NHWZ8JNQPJvmFPeJjz130k1lBht846b/296b13LQIQAHchlKOrySqkzq0VL8F5edYz+bIwmrzI1FFPYauRLV0Z/8cKGzYWE+lAi/CA2fAOKcaW8Mk1YyrDUaZuVR/lQ1MtbilNe7b8WXHTpaSjhedXAnDVcdOxB13CaJ1CKD3prVJJ/fCuqn9smVH6zVehpIXnjazfBa7btYuFqphECdCDHZtVOm0kyjD5t6ybWUBoEY5SXG/2S3GSik2SlfCUINXmF+FAqZQc3ZnU1jHhmkyJbuXr5LXVHrUMpxSnfkkkDAHqzfqEQ5+6oc0wCpRyTl02gqLR0QKD4B237jD1swhHX8Z1YQ/H1Hfs7FjzsvZvF5apBO9QnVINZvmpKKSYvn4bCAqV+8WkolpffGL0U47maw4btGeOYVBZt3INXVu4ors9tW7IXKa7pwfiCCuPPovJsL0LOInBtK9pt2ovpPH+56oc5jwoUCPrW4oxFqj3G/eHiKiiG89R9MpmmwkKK9fqiTKGlDMB6mcl9KObf5rLNEw1z2+y7J7ufARhIk5fqaJY6/7p8G9Z19RSl9zeqhkJEE4noGSJaTEQLiejbMn0kEc0iouXy7wgtz1VEtIKIlhLRBVr6qUQ0Xx67maQeSkQtRHSfTH+ViCZpeWbIOpYT0YxSryMu0etQnL+mhxjnwZo2a4zKFTnjB3DRzX/Fpbca1jFEfHo2bJ2Bju2yjANzxH5ownJMoc/AwgYOm7CKGuTNAkXmDXXoq/OKcoa2xyvb+VtOlFfSveW8KC8y5tHbFVZuqM8h5mQg+eaQ9olWnLJNfTpqPVYtQm9V820ayudvew3v/49nitIbVqAAyAH4rhDiGACnA7iCiI4FcCWA2UKIqQBmy/9DHrsUwHEALgRwCxGlZVm/BnA5gKny34Uy/TIAO4UQUwDcBOBGWdZIAFcDeDeA6QCu1gVXNdD7qFFoqBfN8nLHIcksTvh+J+s8URpK1Iwt+DuI6T74zGUW57ftRdBfmOBptu3rPZNXhIZikKOeoDPnCfNFxNJQXKd8+DqU6DKKsTrlIxzQtrqj1lvZNJScZTIAAPPX78akKx/B4k17DG0z/1Z07eu3HgfsJq9wp7y5LMW8NTvx+ILN9pMSYv04moapbQ3rlBdCbBJCvC5/7wWwGMB4AJcAuFOedieAj8rflwC4VwjRJ4RYBWAFgOlENA7AUCHEy8K5k3cF8qiyHgBwjtReLgAwSwjRJYTYCWAWPCFUFfQXwGbyitq6I2rGk8jkVUaH8bYaMR8P2wdIT7VqKMb9tOzCMqnJK0xDsZmtQjfDDClTT4s2QZr/b8umjpUTNmwq3zPBmjRm528q8PbrTQgbjKI+AVyquRJwPswFmHcyjprI3Kqt2Qh7TqZ3T1kGoqK8wvjEr1/C1383z3pOUlSVpZjZNmqfDWg0DcVFmqJOBvAqgIOEEJsAR+gAGCtPGw9gnZZtvUwbL38H0315hBA5ALsBjLKUZWrb5UQ0l4jmbtu2rbQLhH+wsDrlLQ7SsOO+chJ+yc5tX8K+E6WhhO9CHH0OEKap2fPaBnXFJe86xNgWPZ/NsR51n2yaVZRDvzhf9ENR7erP2TU6G7a1JmaHvdPJgvu/xdlevmCZODnp4e3UhaapXRn1NciIr5ZGzcxDNRRDnbbovrB2Vpsop7xtUrpT09QaUqAQ0RAAfwLw90KIYl1VO9WQJizppebxJwpxqxBimhBi2pgxYyzNsxO1DsW1WVv8K4DlIRu+8R7Vl8vpLqYgAJ3wfYCizSJAdLSUMa8Iz6toyXhdNliGzWTm+RRCi5ZlhAu6aI0xJJ/lSakjZWkohjTPBBt+LNguf98L0UBcjSrsuEVDyevPv/h4U9p5ttmIvmO6n2M6W4zn+uo3lJuRdfYbNo4Ma2e1UdcXNlGxfSa8oRc2ElETHGHyeyHEgzJ5izRjQf7dKtPXA5ioZZ8AYKNMn2BI9+UhogyAYQC6LGVVDb2Pml9+568xAkVEP2TjOpQIkVHO5Cnqi41xFjba6o8y/ZnlSfgg6JVhLk/Pb9YgVb32m2b8NIE7AIcQ8ubHej7ynHJ8KEnDht2+GqhS/857uAYSrgXa0oFox7mKljNpKH7NuLjsoa1N1rIB8zU1Sa2oN2cepKutofz5jfX46ZNLfWmeycvcsawCJcZYU03KifIiALcBWCyE+Jl26GEAM+TvGQAe0tIvlZFbk+E431+TZrG9RHS6LPMLgTyqrE8CeFr6WZ4AcD4RjZDO+PNlWtVQz6Y5nTLPgC0vWjKTV7TZQVEJDSWs08X5lkJSk1fU7KngDnTh5drMMrZIu9hOeYvmFLmHVOCJeJpNdH1Zg8mrnLBhmw8lTqh46OaPFkFlSwf8gsLULtuGpb6JhOHG6DsJhzXB1C+Uxru/P0RDqfKg/J373sIvnl7hS4sKVc6GCD89L1Abp3w5H9h6D4DPA5hPRG/KtH8BcAOA+4noMgBrAXwKAIQQC4nofgCL4ESIXSGEUE/xGwDuANAG4DH5D3AE1t1EtAKOZnKpLKuLiK4FMEeed40QoquMa4lEvShNaTKGeMaN8grroKYB3q8NGD57msA8VlyffXPIKOe1c054+bYorrDybRqGqYzihY3hQj1Sy1DnleB/UdpekckrxkNRZ5QTHWg6ywsSCT8Wdv9MxxTq2ZSysFFf+Jl0nI76bMK+fk+gJOnTKmmfYWt7J4/zdyB3G1atDPvGkWn8UeRjaJnVpGSBIoR4AeH3+ZyQPNcBuM6QPhfA8Yb0XkiBZDg2E8DMuO0tF9XxmjJ2DcX48a0Ydk3zOhR/PuW0rATSdByqicTZvr4cDcVUvKtFxNRQigfw4nqCdYcVbXO8Rznlg3UkwaY5xRYoiU1e5nvh93HYB+VSdhvuzXpagKl420fifBMJgwiN8s8A9mCa1Tv2GfPUxinv1Bm29Uq/RUPR3zv+Hkodox5yJpWyflPe9gLr50XVA/hfImP0ke93ss6jtt2Ys9qs2OmLwHyzw4gFZl7+cC0uLK/KYv+mhv7bPMO2fQAt0odiaVdYXlezCti244xFtgG6nLBhWxCCq6EE+sxtL3iht1G+tVLChns1x7d1ex1jml1g+CYrCYSdStqftTvlqzGZCMPzoZiP23wo+rU3lA/lQKU5TdbNIQvCbkoI1VAitl4xLrgT5t9xUCavx0IWZekzPl3FjnKOuseM7Y2aRdpnv06+8DLUf82faA4X+GHle+UqYRSSJ8QM5OWzDJ4WIRr3mRpn7DYNJUTYrOvyPugUtS4jbK8322Sgz2fyKj7PvvtxcRt89VpMxW4+y0QjyocSNRExBVWUimtmC3PK2zSUiMCHasMCJSbq4TRnUsaXyeZM1/8fNWsIezGitjJJiu27IsF2hH37wjbwm01e5vIVtsHVOydcy7HZ9217renYfD9R4ahFAiXG4/E+WFV8LPbmkJa8RlNajHsR/sVGu4Zi6xP9Pqe8rW57WpKtV/Tm2Exe+7P2KK+oR1GuQDH167B7aavLtwCbNZT6RT2bpnQKWYvJCygeCPzquLl8c9iwuXxfPio+Nw7BnWaD6PWVoqFErUOxRR9Zd6zVnY6B01T0i30dSpSGUpwmIvKq+oL+KHW6rUZXq6qwUz6OycvWsig/RClhw/5tc5JpZJG7LISYlaMWU6psvZEmrygNpbzB2zeRVGkhdXb19BvTgcB9YoFSv6hBpSkkbNgXXWEzeYUKhuKtUHyDd8jMuVrRJz4NRZ9ZxvShRJkljAO3Ohbb5OU/T83cbIsEw4rukVFCRs0p8DeIp6EEt6+PfqFte2PFHQ9M91oJN6PDPkTYNGmLRqM2Syxlq5L+iNXstrxRa5j827qYJ3c2v16YyUtYBLOv/jI1lJyhnab3oD9XwJdun1OU7uU1lzlQsECJiauhZFJW1Tn4G0i2ejXsm/JhJiK1PUPSfb2iZlz6bDtsD6akGkqkCU8mleqUVwIl+G4LIayO9QUbdrszTJvmFNYsZQIN9guvTnM+vexyNBST3LKWG2IOO/3wkV6RYdqY0iJDtbXwZmYjoshsA2BUBFocM5fZT+Wkbd7TG7Flj/1Z9Fn8GnEwCRRTe3dZtBOANZSGYO2OHjyzxFnw35wmCGE3a1lNXmEaijoeMosPm9HH3UAu+EJEdTW9zWE227gzSrdMi9DV0+xfbAyfqfa7Ji9/Afq9MzX5rfW7jG0M5on0oQQaHscHEndXYKtj3+aUN9xLk9lqy55efO/PC7y6Q6rT/S+24AcT+vYmtlBn0/VEvUOhPhT9vTStH9OOr9zWbThefJ6JcrWBvC8IxjwxAoo/4lZUToQArTblLGw8YDj3pufcwUrtN5QrFJBOpd1zfIIg8Mx9jrKIXXx9/TJiAPYWO0YbVwoC0JexJAmfzYX6UCwCJcLJHNdebjsWPK8/RFOIGoyiBHeUpqEEWLDenGWADJYdJYCD65CiIuZc05RFcOsThRVbu43nWNskBFLBDSYtA6tPQ7EIOlPVUXvc6ZOIQojZJ2rD0t37s4bjZm0uSLlOed0vawuoCAtvVqhJTXOIJaXasIYSA92Z6AqUgGDIC4FmeSzYceNoKGqA8Ju8zGXoxPWhBPOr/7Y3pw1n+1/EsCgv2ztmMy+Y2qOXZ3XKW7QNb2YXvNZ4g2/wt1ZTUTk6noYSfO4xPlJm8Uno1YWZ04LnuXVbTDWqLn2Abwv0gzgmrahthoL05wruVidmv0+0DwwIEyjmdyxqfYp+f7oNq+VtXwHVsYXyxsFvEQjvE2HBA4o+eY/TRGzyagQ8DaVYaKgZpG1AC910Uf0NmYlFrV+J6vDBF1j9L0xjymsjhy/Kq4zdhqMEqyo7rlM+OOj1h0R5mSJodKI2LVRJ4RqKeQDwHOPmfIDmII+4X1v29PqORUUvqbpNt1I346g61GQIcELKo7SxsHqtJq+8LlBMZYcPpFE+lDCBE6Wh5AvC3ZTSZE6yaZC+9pU5ePu/BR8erRiloSiBkkmZ18tVGxYoCWnOFAuNfEGgL1dAc8Yzh+n4ZnWhGopXliIqskUI/ZsJ9s5T1DllgWHah09DKSFs2BT1EmlaUgOdzV+glxEoJExDiVyhb/F/OW2Np6EEzR5l+1C0+rbu7Qst21SLddcArVzTPWtOp2KtNrftSGAimy+gpcnRhEzPQU1cbHvlBdtvapdv0pGPeraelm6a/YdFxAXpixjoozBpKKbr1NvYZNiKyRmH0kilWENpCDyTlzd43P7iKgDArp6sPBauoYQ9ZHfr9pCBIswWHrWeJCy/N/Myt8k36OjbsIScEyRqIDN/dyR8hmoqw6dZFIRvgWHYCukoh67pktThqA0Rw3wo63fuL8oTbE/USvng5oW+uizXZBro9Pvb75qZvLRMisIXLloEOmAX3rrJy1R6fz5f1D6vzfaQ4zAt0+Sb0CkIgfZmx5VsEihRa5AUn771FevxKHQh2ieDF0yCVW9jNi9C77HtGVYTFigJMZm8ggNGcKaqv3hRmy6KkHEiLIY+lTLvdOukhQsy3XRlalO4hmIXCoqo73uYboM6HNcpbxoYTSaVqIWlUeuE3K1VQhyvYSvllfb33LLwL4Xa1hzo7d7X5x/s8nn7vVTlmWz7vvtmMBMSRZv3gr/jHM/m7T4UtYW/qe/47619shL27MO+ddPeEq6hxPmkQiXQ29kXEq0IeGa5j53sfKA22Ky+XB4tmRRSKYqlIVcaFigJaXI/U1rcARRF/pU4GoooPjfKXyGEfQsV28sdKaz0dSihYcOhVRt9M1EhjUmd8vp9VoOQci6H1WUqOW5Iaqi/SUV5BU1eEaun9fUxUf4INXtX+LRGm4YSIVBMi0HTKQqP8kogUILvQX+ugFZp8jI94v4QkyUQvZtw2KQncv1KAeiQGopp+xXbjgOVFDK6EFVC3tTflA9lSEumKJ/K2yyd8ixQGgA9bFjRF/h8aPAh+xxuYT4UUySRbl4IsRunLFuv5C0agf5/01YyvsE6xAxn06hNvpm4Jq+CiF65DPhncMo80N5ULFD8uw8U16uLZNM7qO6j6T4B3osfnOFGOUXDZtIK0yBjOt9UixIkNr8A4GkFenkpiilQLMEAwXMB5xm5Ji9DXnWNUaYe2ye4nd9ae7RnZtzjTQi0NcXxoYQLbQA4aGhL0fEownx3SrA++MYGLNy425dHtXFIa6aoDYB0yjelkU6RdT1XtWCBkpDhbc6nRk0qqiLY4fXBIDTKS81Uw5yLIfbflMWHYpt5+wbmiBfU/6U98zlBTIOCfx+ucA0FANZ0hX2fwvtOhH6f1T1WM2B9IIly6Oq30LYjbbiG4qTvCqxjiJohRpkATddnOmbKqwYek4ZSMAxeamB68Jtngih8MMoL4WrpRiEYct9VXZ6GUpw3jqlHb7NOrmBuV9DPFkQIgXSK0JJJ+bbX947LvKaJhsF0mISsT9hpE1TtWp9YuMWXJ6ihqHd3294+fOA/nsFzy7ahJZ2SAmXgJQoLlISM6GgG4A2Yd7+8Go+8vcl3TnB23h+hdjumDzVT1Msxz7K9fJ7JyzS250JeLKee8GOAfyDyCwf7tXhlFptk/Ka/4jx6eXt7zV/QyxeEqyWaTDfugKULr4J98NWv1RZBFHT2K9T9U0EZwfQwokxt+sAZFAz+MNnwGX1QewYCpjQ1iMtn3JRKIZ0K948VCtp6qwjTVHCi4luHYnj+WUOAQPB6gPBAg5aMYTIRYaIryHeotSmNXoNWbNt6JWfQkJPg295Inzzozz1wrUqwuiYv2YZV2/dh9Y4eAM6ixnSKQj/jXE1YoCRk3LA2AF7n/N0ra91j//HJEwEUd66sT0MpLrOnP2+0pftezlCTl0VDsfhB/AKluFE92idV/cIBxt9BsgbTh4jQFASATqnK7+0tXrWsyvK0EF1LdF48FQIaNqjYTC3Bc01ptq3xd/X0J/rAkRpA0yEOVJ/JK/D8erTBL3hJ2XwB2bxwP1ddtC7K6EMRblusJi8h3E0kjQLFoP3odbUZno97vsV30JvNu8KoN8Qv1NpUvOg4ytycLwgQAW1N6cTrUFQ9bU1p60evwgjbL0/vj8H+1pvNo7UpVbTAWn8WQ1oyrKE0Ape/73C0yE6rZgYpzSk+foQjbIIvxHf/+Jb72/SQ9UVroWGxYWHDqfDNIU1Oa7cdIYJLoQ9YPodnoP4wTB/osu3IDDgD49BWx6QYpqHkCp4dXh+UbnthFQBghdyPKdQpH6KhqMVtYdF0puvS2+RcE7C3L0wQmwSGk9aSSVkHLKDYpHLlg297ZaN40AGAYW2ONh3UUnIF4a6ZCq5DyaSVQClqjnueGsxsq931dujX0NnSJNtkECiWlfK9uTyGt8u8ERqKPrhHaSjK5NXaZDZ52XYyUGNAe3Ma/blC4g1asyH+Hf1Zq36tcARK2u2v6vnp7/eQ1gw75RuBKWOGFA082gJjo8od7GSmicz3/2eBdtysVYStFXEFiqG9ejv6c/4zcnkvr3FFbn/enfFlQ9ph66++umX+6LBhgaFtEQIlX3CFuj7gbt7tCGUl3k2287Dopb5cwa3XNEv1mdaM24J4adv29mrpds1GPV9HoBQdtjrl317vOWuDl6SuwRuAi7VT5YgODuLpFIEoPIIpr5m8osLN9VXdap2Q0kBNDnCbyWt/fx7DXQFp9qG42k/IzD/M5JUix+RlCgRRbTJOROR97QxxkEexR/O55Sza6IZdzrKE219chbteXoNMKuXuymFaVNvZmmGnfD1z1tTRAIBPTZuATMo/mOkdLV5IcfFTVoOhczzpyxDe7j6LI9M/MJhMXnlXW+gNMQlZv33h88GYBIpZExjRbhco+YJAa6Y4kuuDR48FAHzlrMOLjikNI5Myz7zvnbMWXfv60dqU8pn6TG0N01AOHtoKAFi2pduXrrDNyFsyaePzVcfTKSoSKH/z7kO19vnzqcFaBZAEZ94F4ZmHgutQMimS9vdwgdJiMXnpA5s+QCtB3NkaLrjtJq8ChrWHazfZfAEdUqD4zFzKN5Q2TyZUpGRrU9poSjNNhhRK81PXlGSDyJ7+HD740+fc//v6SkDYbpIC5VfPrAAAbO/uc30o3fI90ftHJ5u86psbPnEinvvHD4CItJmB87D0l8YVNtqDVC/3F844DADw7Xvf9JXdm81j5fZ9aG1KYVhbk3EgBoAe06Irn8mruN19IWU5bRTeTD9kRnjwMGeQXL+zx5dPYVPxTeY2vzAy5MkL1/7b3Wf2oWQLBWO71bWq2aJf0/MGS9POv+oyOpoz2GcQKDmLL0rVNXGkY+5UGwwWCgJCACNlEIcpCkiVq66n6JPGKhS6OV30/FozaXcCE7wmpRnE0VCU2U33oWRSFBrRVhCaycswXumah66hqOtvb3bMNSYNRT3DMKf8sDYlUAyaRM7zz+j3Sr2LYUJbRUq2NqWMbVLtNvXX3oCGks2Fvw9Buvb5v2sSdMp/8wNHFJ2rxhcAGDXE6Vc79vW5eRRDWjNIpQj7QkLvqwkLlBiMH96Gw0Z1AIBr8lKdvlN28ivOPsJ9wXWHtPJFKEdykMvvnuceb86kfC9LvuAJjD2GrbXzmi3c9LL4tmkIDGgF4c30jT6UbA6jh7RgZEezb9sP/6aExkuSZWrCTL5oYYvP3PPyBTSlU+hszYRrKHnPVm4K3TYt+FLntTSlQwdKAGhvSaOnz2yfd9sYEtmktLkeea9U/1BBAsZBMLC637SyXM2egxpmby7vCoVik1fQhxJ89l5/LPKhpFJozqRDPxil9zmTZqtPsILbhACOptAa4gBXwtgktHWBYtRu8sJdoJg1mLxaQrZzVwsBnTbZBEr4+6WefV8+/gAettlnoSCQzTv3+NbPnwoA7iRHaZUAMKLdebY7e/qxcONuPKmFF3e2NuGwke1Ytd0cel9NWKAkJGjyemud83GmfzjvKGQMJiT1jYVDR7Yby3tebssx44xJaMmkfC9yX66Aw0a1oylNWLxpb1FeIbxdYk22fd8MMaih5M0zfUVXdz9GtDejrSmN/f3ajC8iBFqRDcy4gm146M2NmHTlI9ga8Dlk0oQhLRlXlQ+S06N5AmYCIu/lNl370NaMdb3Auq79ePCNDUXCTt/WPDjYZfMFbNi13/XBqFmhMp2pQdBUrxoYO9SaAsO6jaZ0Cs3pVNEA39ufd/egCrZ3kzShKvNhcIfavOZvUIOirqE0Z1LGqKVsvuDbTNGo2crnAMDXb9T1N2fSoQ5wNXEx+vSyjgANTrr0trW3KF9G8WSiOWPe8FJt6toWIlDUfTd19SINJUGc7n1z1vn+r96rldu73faefOgIAJ5Za9QQb/HkqA7n947ufnz45hfwyHxv6UJHSwaTRndge3ffgDvmWaAkRGkMmwMzjHSK0CH3BNJn1zvlJzsPH9OBi086xCdYfv3sO+7v3lxevixOJ335nR14eslWbNi5HxNHtmNtYKHfG2t3orsv5wUCBDqzEAKXahvWBTu73xZeKMq7rbsPYzpb0BJ4+U2mJBP6IOw6W7VBat6anQD8H3bK5gUyqRQ6W5uwxyZQXA1FMy3KNQ5qcNb3vlKD2dC2psjvSQDAG+u8Lzjm8gUs2bzXfe7B+3jHi6sBOGHOTWlyB8Xt3c5zV2Hmphn/vz26GIBngw8OWrm84+dqyaSKBFJvLo+OljSIioXV16TWO1KaRWYvLl4cN1oOTvsCWkE6RWhJp4yRVOpcldd0L59YuMW9Dl2QKQHbkkmhJVM8eBcKwtXmwyZHrU3OvQia8ADnHqhdEvxmY0+gmPxCznbv6VCtaY8MX+/PF4oEkroG14eSYHHj5NFDfP9X79WX7nC+F18oCM9PIvuyrv0NbXNMw+u6ehAkkyIMa2uCEAidmFULFigJUWatH/zPAtw3x1mDohyyozpa0JJJYdkWT5t4ccV2AI6KOrKj2RUwAHDj40vc3wRCSybtviyf+a0jDPpyBQxrayoyAX3slpcAwDM/BGaUwa/PBSO19mfzmnnI/6Ks2dGDbF5gTGcLWjNp3+Ciz0pNHyRSzN/gRSF5oY3y5dZC43RTYL7gfJuiszUT6kPpy+a1tQx+DaW1Ke0Kdd257gqU1qaigf2nTy51f39m+kQA/gnBb/+6ym0bUPxN7z+9vh6Ac7/bmzPuoLij27FtHzbKmUCY7tXLK3fIdhWb6YQQuO2FVdjbl3M0hqBAyRbQ3pzBkJZMqPBVk41fPL3Cl76/P4+xnY5Q2CkXY6rB39MCwk1Syn4fjIoKCgldoKyTG6hOHNlu1EB131VwQeTWPb3ozRYwblib844E2rajuw/7s3kcMlyG7fuCWTyzojHwQW6m2NqUMn5rZJv22YCgf01NtIYEnl8uX8Ddr6yxTl6C/UgJ0XVdzn3a3u0Eiej+RPVO3/7F00BEGNHejDfX+7dmAZxAnU75bu8JWc9VLVigJEQ3Bfzzn+YDAD59mjMQpVOEkyYMx9LNjkDZvT+LXz3jaCEj2psxvN0RDCYbcYrgMzWcNHG4U/a0ic4LGDJ4q84cNFEEd0DW15V09+dQEMAYOagEB6vnlztmuPdNHS2dleYorzA/R9e+fry2qsv9fzZg8hotByQguBNzAZk0obPF7EO586XV2Nfv2dKDPpSWTMo1A+kOSS8aJ4NcQfjuvxpsrzj7CHzpPZPldXkv4YZd/hng8sCncr0INqCjOe3O4ndIR+pRB3U696TbP4DojJcDof6MtnV7A5nJBKXCuoe2Fk82FO8/cgwAYIJcHwU4gqqnP4cR7c0YN6zVvZ59fTnpr0lh4sg2rNq+r8iUprQ+V0MJ9BsV3vo+Wa++8ny7HJjHdrZg7NCWIh+CrlEG/VTKGjBhRJs0C/sHarXljVoHlg2YjQGgrTljFCiq34ZMTSoAABpVSURBVJh8KN19OfT05zFJTgqCOz4/Nn8zAG9CoOp6ctEW/OB/FuCXAUGuKBQEfvP8Sl/avz2yGF+8/TVfGpFn/hVCYGdPP778nsk4W0Y0jupodk3uOh0tGUwe4/h83zAcryYsUBJi8oWoAQ4ADh3VjjVyC4S7X17tpo/tbHHDOPf05tyXT6H2E1LawIThbTh0ZDuu//gJTvSRJlB8qq9Ut4Pmp6BA+X/a4krV0SeNdjqdrjWt2bEPP3xoIQDgiDFDil40fWAzzX5y+QJ+/tQyX5r34SQn70hNoOizzf6c3Sl/9cMLfdeqC8KefsfG7moorumpD08uckw+rvPUGHbqRT71GO7v0Qd3IpMirA2YGNQgNXFEG9pbPA3llmede3zkwY5AUdE4Cl2QHhrQYoQQWLXNM3E2pc0mr9amNDpbM0XP4cwjRuHkQ4djytghGNPZgvdOGe0e6+5zJhND2zIYO7QVu3r60ZvN41fPvIOCcAaxow8eit37s0Uf9QpqKMEB+BwZBvupUycA8Gso26WAHN3ZYtS4lelz3LDWooFf5R3T2YLWpmKT18KNewBovgx9LYx8Hp0txQKlP1dAX66AIS0Z6SvM+4So0k4my/ckqDUrv8VYaaFQ/XujfLfDdntYt7PYTLWvP49nl3qfOfiGjPIa3t6EnT1Z7OrJoqc/7wpNAJhy0JCicv714uPw3imjccqhI5BJEZZs2mNsQ7VggZKQztYmfOj4gzFOhtQC/kH2sJHt2Lq3D/v789iyxzN7pFKE4TIy449z12F9YGAicgSKKqu7L4fh7U1IpQgdARPBdm32qvwgOwMqtAr1VeGsit37s7hVzo4Oly/KDm32rOzvgLMLwPwNuzFX+jsAb6Ae1tZkjDx7btk23PXyGgDAR04cB8AbwPtyBTSlyfdRMBUSWSgIdPflMLQ1gyGt4RoZAHT19IPIi3x7eskW/OXtjRje3lykoXzt7nl48PUNAJxBVG+Pzspt3a7/5RdPLwcA/Gneetwi/Vx3fGk6xo9oK7JZT5QTjO99+FhHQ5FmETVgKg3l+eXb3Tzb9vbhyO8/5v7/oE6nL6lJww2PLXE/2PSHr74bbU1pn5ADHJNXa1MaQ1uLn8O+fs+c2d7sz6t8O6OHtGBsZwtWbd+Huat3+vKryc51jyz2pbsCpSPchwI4odLN6ZSv3h37+h0fV7MjBPcGnu/DbznPaPrkkT7TXzZfwKsru9w2O/41//V+6543ADjmuqa0f82OEmpD2zJFWp4a8DtbMxjR3oyc7IMKT6AMkdcfcr3t/v39HpjnmEHTqRQenb8JZ14/22cGf/kdx9R50QkHG8u79pLjcJAUUuOGtWLhxt3uBFFpswDwkRPGFeWdceYkEDlriQ4a2opbnn0Hp//bbPz8qWU+8121YIFSAu87cowbSQMAR8lZKODNNtft7MHWvb0Y0pLBU//wfgDeuoDrH1uCHzy0ADppIrQ3p7G9uw9CCOzry7mDwpCWNLo0gaHPHE8/fBQA4MePLHY7zOMLNuHHjyxGSyaF5//xbBx9cCfePXkkAL8wGj+8HSnyx8QfLa9leGCBoZpR3/2KIywmj+4wfolwp7ZB4nfOOxIA0LWvD919Ofzh1bVobUr7/CtqtrbXnTk3oaM5g659/b6B5WezPK3nO+ceieFtTXj4rY3Y25vFl++Yi4Jc86EikJTg06/N01DyvmsCgKcWb3U1zXVd+/HSiu2+LXMOHtaKI8YMwZLN3sCwo7sPf12+HeceMxZjOlscH4ocdMZ2tuD0w0e6/p5H3t7kDlYPvbnBtfP/+KPHu4JMTRp0c8i7J4/CxJFtPs1oxszXsHjTHrQ2pYuE3JvrduGtdbtc/0lQGLmawpAWnDB+GNbv3O8ujlWmyCljnAFUX3/UnytgxkzHJHPIcGew0wWK/nxOnDDMt67jL29vxK3Pr0RLJgUiQkezf4JUKAjcP9cZhI8YMwTdfTm3nffPXefejzGdThi7/kyVjxKA1G6bfJrBrp4s2pvT6GjOFPl8VN/ubG1yN33duc/Lu1QKgcOl+Uj3e2yVZrh/ueho14/Zm81j0+79bh+Z+eIqfPP3r2Pj7l78Vl5DNl/AlQ86pvL/d/5RWHLthW40nkIJEwA45dAReGfbPry80rlO3Xx5ymEjfPnUeKFQE4PNe3rx86eWh34SopI0tEAhoguJaCkRrSCiKweq3rOPGhv6f2USu+vl1diwaz9OPWyEuxBMN5epFdXTJzkDfYqA904dg3Vd+/HOtn3o7su5A02uINCbLWDSlY9gw6797st255en48LjvVmOskt/5z5nIOzLFUBEGDu0Fa+u6sK+vpz7ErU2pTBt0giM7Gj2md929jjRSs/949kA4C6w2p/NY0d3H96Rpph3Hz4SCzfu9jm/L7tjjmtae/OH57kvxl0vr8HMF1Zh9/4s9vbmcNsXT3MHfiWU1JYlw9qaXBv9V++aCwBYuHE3bp7taA3f//AxOPaQoejLFbB6Rw/Ov+l5t/7h7U1umK3SUNT9BeCG9v7Xs+9g6ea9PrPLjZ84AekU4R+kEPyb/34VQU6aMBwrtnbj4l++gBVbu3Hqj58CAOzZ75TT0eJoKEIIzFm9070OxfUyqkt/sYe3N2nRPMVaWTpFOGxUB7r29WNvbxZCCPcLkGM7W3DI8FZs3dvnCseP/upFAJ7Zpb05jf1Zr1y1NmH0kBZXe73xCScw4Vd/cwoAZ0cIAJgy1jOpPLlos/tbOb9VWHA2X3CfD+AM0G3N3lYmf/sHR4NQwQNDWjPYn81j0pWPYE9vFv+p5T3jiFEQApgvnc0vaJpda1MaI9qbfbs6f1Z7TqccNhyjOppdjXtdVw/ueGk1jjq4E+2afwtwJgMX3fxX2d4MRnY4fUOfuM1f7/gflC9qtbauQ5nZTj50hDvIr+nq8QkkHTXRekUGYgBOBGBrU9o3CQO8+wsAU6VZ698edQJ4dIEyRutfP/nUSbjva6f7ylEWCMXYEr7ZkpSGFShElAbwKwAfAnAsgM8Q0bEDUffBw1px2XsnG49Nkgsgf/fKWizYsMc3eJj8L+cdexAAp+OcNsmZcZz7s+ewZPNed+Zy0oTh7vnvueFp/P5VJ7pMaRPfu+gYAMDsxVuxZPMed33ITZ8+yXfe393zhrvNyz1fPR2tTWmM6WzFn9/YgK17eyGEwPwNu3HRCePc2boy7V398ELXWX/0wZ049dARyOYFnl+2Db3ZPAoFgdlLtgIALjjuIAxvb3YHyjfW7vKZIc4+aiwWXXMhTpowDIs27kGhIHDuzxzBkCsIfEZuK6Jsym+t8zSaCSOce6hm3bqmqAR7R0sau/f3Y11XD+6b68X7K/PgnS+vwQU/f971a5x/7EH49GlOnUce5GmbCrWlhwpieHv9bpz7M2/bjM+efqisN4OFG/e437AYJWe9agb7+1fXYsbM13xa4hmHj3J3JHjw9Q2+gIH/+pwzwKtnsXDjHt9itU+fNhGjOlqQKwg8vnCzz/5//cdPAOAM7uu69kMIgQUbduOfHnjbvZaJ8l4qx65a1U9EOPOIUXhu2Tbs6c3i+WXbfOHdI9ubMXFkG+aucUxRepue/X8fkPcsg+5+cwCKPpN+e91uV6Dc8PET3Da9uW4XhBB4bIEjyJRPYdSQZmzd2+v0Ny0c+sFvnomWTNqnwfxI+tyaUimM6WzBnt4cFkjt+KE3N7p9aNywNm+hoMy7ZU8v7p+7HmceMQoTRrShszXjTqZ6s3k3vPfIgzoxblgrOlszWLRxj1t3W2Ah81OLt+A3z72Dz9/maHn/cN6RrvY6LaBp6ELjQ8f7zVq6vzaVIlz/8RPw0BXvwSdPnYDjDhnmO/f+r5+BmV+chsXXXIjH//6s0MXVlSQTfUrdMh3ACiHESgAgonsBXAJg0UBUftWHjsZtL6zCJ06Z4EtXqrPi2EOGur8z6RRu+vRJrgYBAF85azLOnDKqqDMA3idJPzVtAp5dthWPyqgStRhSDQBfOWsyrnt0MW56ahluemoZOlszuOiEcfjYyU7bvv7+I3Dr8yvx9JKteFoO+sq597dnT8EVf3gd06+b7db7romeAFPX88C89a5t+KG/fQ827nIG8q//7nVfmw8f04EbP3Gi+/9LT5uIe+esw2//uhKjh7Rg7vfPdY+9d+po/OqZd3DSvz7ppn3s5PG+jj/pykd85X/gKGe2eJgW/AA42pqaSe7an8U9r63DPa95wqSzNYPzjzvYNTcAcPdS+r/TJrppRwYcnUcd1Im7LpsOAG6orc4Hjx6LS97lfN9bhY9//Xfz5PU57Xn9B+fh+KufAOD4mJ6T9Tz5nff7ynp84WZM+Z7jW7n+4yfgQjmYKJ/Tpbe+gnOP8bThqWOHYLF0un7z995z+Nr7D3ejBI86uBPPLduGyVc96qtrTGdLkYlE10jOPGIUXnpnB0780ZO+cy447iCkUoSLTzoEv3rmnaLnowI9mjMpPPL2Jt+3glbf8GEAwBFaPZ+7zdMwPnryeFeb/8/Zy32ayz9feDQAJ6Iqmxc4/F8ede/3n75xJk6RiwCb0im8sGK7r13/eslxKAiBnzy5DB/5xQt475TReEEzlR0/fqhrUvzSHXNwxdlHuNGZJx86HESESaM6cPcrazD1oCGu9gR4A/wph47APa+txT2vOZO933/13Xj07U345tlT8MjbG/GDhxbi+seWuPf5W+dMdcu448vTsX5nD37wPwuwZkePT2i0NqVx7+Wn43t/no9rLzkeFPhcxWemH4owRg9pwQePdiasRx88NPS8StKwGgqA8QD05abrZdqAkEmn8NbV5+PGT5xQdOyv/3S2+/s75x7pO/axkye4g8IDXz8DROQTJv/7t+91f19z8XEAnBnjLZ89FQ9+80z32FlTR7svX7CT7e3N4eRDPaEwsqMZ93zVU4eb0oTR0rF60QkHu1EsXhu926gGacWQlgxaMmlMGtVeNPgCwL9/4kQ3+ABwtvwHHPObPhgCwPnHOuY65aC9/2tnuMJEvw+KVddf5B7/0cXH4ZJ3HeIeO/OIUe7vL5x+mC/f9R8/AfN/dAFGdjTjL39XXK4+wKktdhQ//D/HuqY7Jcx0Zn7xNPf3/znpEN8xNdMc0pLBLZ89xXfs5In+WanSYhQnTvD6xJlTvGt7arEzIZj/o/NBRDhNM+kpdFOHCmnXefOH5wEA2prTeOzbZ+GkicMx9/vn+vrRNz8wxV3LoPjWB6fgN5+fBgDuZEW/1reuPt/9v9pQVfHot85yf5va/Nq/nINW+enaIH/8+hnub/0eb97Ti5EdzThVm+F/+9ypvrzXXnIcjhk3FMcdMswtWwmTY8YNxes/OA9E5C5ABeAKk0tPm4grzp4CwPOT/vChhfjjvPU4ZFgrFv7rBW4e3fQ8eXQHThw/DN//yLEY2dGMT592KI4+uBPD25tw7SXH+e4F4PSPow8eitu/NB1Pfud9Re/z6YePwuzvfgBnTvHf03qEku7hXy8Q0acAXCCE+Ir8/+cBTBdC/F3gvMsBXA4Ahx566Klr1qwZkPa99M529OUKRf4WwAk1zRUKruM0CWt27MOunqw7A1Vs2dOLe15biwUbdqOztQnXf/yEIhW3uy+Hu19eg4+cOM6NTgLUFyOdrd/z2uZ/Ouu6etDSlMKYIS2+Dj9vzU78/KllmDy6A989/yjf7ErP+9iCTfjc6Ye5UViK9Tt78JvnVmL8iDZ87X2H+8p+Z1s3lm3ei0Wb9uCKs6cYVfanl2zBseOGuWYjdT39+YJ0OFORwASA/31rI757/1u4/+tn+DQywHESPzBvPRZs3I1/vfi4ohd8XVcPHnpzA9535BicOKH4OSzauAcrt+/Dl98zqSgvAMxZ3YVTDh3hGzxz+YITuPDaWpw2aaRx0P3l08vxkyeX4ZbPnoKLtAifXT39+Ovy7fjRwwtxzjFjcc0lx/vu1Yqte/Hg6xvw3qmjMWZIC6YazHomerN5zN+wGxt37cf6nfvxzQ8c4bue3T1ZvLpqB3b1ZPGJUyf4rqcvl8fyLd3IFwQmjGjzbRsCOM70/lzB9TMeM86bQQshsGVPH5ozKezq6cfhY/wTFyEEVm3fh4ff2ojTJo3EewwD7eMLNqMpTTjnmIPctJdWbMdDb27E4WM6cMFxB+PgYa2++zRvTRf+9PoGnHP0WLQ2pX3l9mbzuG/OOizcuBuHjxmCD5/gf4fyBYG/vL0R44a1YcrYIa71QKdQEL7vJzUSRDRPCDEt8rwGFihnAPiREOIC+f+rAEAIcX1YnmnTpom5c+cOUAsZhmEGB3EFSiObvOYAmEpEk4moGcClAB6ucZsYhmEOWBrWKS+EyBHR3wJ4AkAawEwhxMIaN4thGOaApWEFCgAIIR4F8GjkiQzDMEzVaWSTF8MwDFNHsEBhGIZhKgILFIZhGKYisEBhGIZhKgILFIZhGKYiNOzCxlIgor0ANgMo/m6mxzDL8UMBrC0xr+1YuXlt7apWvfXYpqjj9fj86rFNUe2qx+dXTrlRx/n5AUcJIaK3WBBCHDD/AMwFcGvEOaHHAWwrI2859UblDW1XteqtxzY14vOrxzY14vMrp1x+ftHHAcy15VH/DkST1/+WcTzqA822vOXUG5XX1q5q1VuPbYo6Xo/Prx7bBDTe8yun3Kjj/PxicqCZvOaKGPvRVCt/tajHdnGb4lGPbQLqs13cpvhUul1xyzvQNJRba5y/WtRju7hN8ajHNgH12S5uU3wq3a5Y5R1QGgrDMAxTPQ40DYVhGIapEge8QCGimUS0lYgWaGknEdHLRDSfiP6XiIbK9CYiulOmL1bfYJHHniWipUT0pvxX/GWt6rSpmYhul+lvEdEHtDynyvQVRHQzmb72NPBtquR9mkhEz8hnsZCIvi3TRxLRLCJaLv+O0PJcJe/HUiK6QEuvyL2qcJtqdq+IaJQ8v5uIfhkoqyb3KqJNFblXJbTpPCKaJ+/HPCL6YKXvUxXaVbF+VUScULDB/A/A+wCcAmCBljYHwPvl7y8DuFb+/hsA98rf7QBWA5gk//8sgGk1aNMVAG6Xv8cCmAcgJf//GoAz4HyM8TEAH6qDNlXyPo0DcIr83QlgGYBjAfw7gCtl+pUAbpS/jwXwFoAWAJMBvAMgXcl7VeE21fJedQB4L4CvA/hloKxa3Stbmypyr0po08kADpG/jwewodL3qQrtqli/KmpnNQpttH8AJsE/UO6B51+aCGCR/P0ZOCF1GQCj5EMdWY2HlKBNvwLwOe282QCmyw64REv/DIDf1LJN1bhPgfY9BOA8AEsBjJNp4wAslb+vAnCVdv4T8oWv+L0qt021vlfaeV+ENnjX8l6Ftama9ypum2Q6AdgBZ3JQtftUTruq3a8OeJNXCAsAXCx/fwrOYAkADwDYB2ATnFWoPxFCdGn5bpcq5A/KUW8TtuktAJcQUYaIJgM4VR4bD2C9ln+9TKtlmxQVv09ENAnOrOxVAAcJITYBgPyrVPrxANZp2dQ9qcq9KrNNilrdqzBqea+iqOi9KqFNnwDwhhCiD1V8/8psl6IqYxULFDNfBnAFEc2Do172y/TpAPIADoFjnvguER0uj31WCHECgLPkv88PUJtmwumscwH8HMBLAHJwZiVBKh3Sl7RNQBXuExENAfAnAH8vhNhjO9WQJizptWwTUNt7FVqEIW2g7pWNit6rpG0iouMA3AjgayrJcFrZ718F2gVUcaxigWJACLFECHG+EOJUAPfAsWsDjg/lcSFEVgixFcCLAKbJPBvk370A/gBH+FS9TUKInBDiO0KIdwkhLgEwHMByOAP6BK2ICQA21rhNFb9PRNQE5wX7vRDiQZm8hYjGyePjAGyV6evh15TUPanovapQm2p9r8Ko5b0KpZL3KmmbiGgCgD8D+IIQQo0VFX//KtSuqo5VLFAMqKgHIkoB+D6A/5KH1gL4IDl0ADgdwBJp2hkt8zQB+Agcc1DV20RE7bItIKLzAOSEEIuk+ruXiE6XKu0X4Nhda9amSt8neV23AVgshPiZduhhADPk7xnwrvthAJcSUYs0xU0F8Fol71Wl2lQH98pIje9VWDkVu1dJ20REwwE8AscP9qI6udLvX6XaVfWxqhqOmUb6B2dmvQlAFs6s4jIA34bjcF8G4AZ4juchAP4IYCGARQD+UaZ3wIlkelse+0/ISJ0BaNMkOI65xQCeAnCYVs402VneAfBLladWbarCfXovHDPC2wDelP8ughMwMRuOVjQbMnBC5vmevB9LoUXdVOpeVapNdXKvVgPoAtAtn/mxdXCvitpUyXuVtE1wJlL7tHPfBDC2Cu9fRdpV6X4V/Mcr5RmGYZiKwCYvhmEYpiKwQGEYhmEqAgsUhmEYpiKwQGEYhmEqAgsUhmEYpiKwQGGYOoGIvk5EX0hw/iTSdn9mmFqTqXUDGIZxFpwJIf4r+kyGqV9YoDBMhZCb9j0OZ9O+k+Es+PwCgGMA/AzOwtjtAL4ohNhERM/C2efsPQAeJqJOAN1CiJ8Q0bvg7DzQDmdh3JeFEDuJ6FQ4e6X1AHhh4K6OYaJhkxfDVJajANwqhDgRzvb+VwD4BYBPCmfPs5kArtPOHy6EeL8Q4qeBcu4C8M+ynPkArpbptwP4lhDijGpeBMOUAmsoDFNZ1glv76TfAfgXOB84miV3CU/D2cJGcV+wACIaBkfQPCeT7gTwR0P63QA+VPlLYJjSYIHCMJUluJfRXgALLRrFvgRlk6F8hqkb2OTFMJXlUCJSwuMzAF4BMEalEVGT/EZFKEKI3QB2EtFZMunzAJ4TQuwCsJuI3ivTP1v55jNM6bCGwjCVZTGAGUT0Gzg7wP4Czmd9b5Ymqwycj44tjChnBoD/IqJ2ACsBfEmmfwnATCLqkeUyTN3Auw0zTIWQUV5/EUIcX+OmMExNYJMXwzAMUxFYQ2EYhmEqAmsoDMMwTEVggcIwDMNUBBYoDMMwTEVggcIwDMNUBBYoDMMwTEVggcIwDMNUhP8PK/1OHMapunEAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ + "sorted_data['inc'] = sorted_data['inc'].astype('int')\n", "sorted_data['inc'].plot()" ] }, @@ -215,9 +2238,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXl4XOV59/+5Z9Uua/Mm2ZbxAthADDbGhoQlDpjQtEAhjdMkOCktJCVpk+bXtyHt+5I3KWlo0yy8TUhIIJClAZoNsrAYE3bwglm875u8aLH2Zfbn98c5ZzQajaSRNJJGo/tzXXNp9Mx5js45mpnvuddHjDEoiqIoSjq4JvoAFEVRlMmDioaiKIqSNioaiqIoStqoaCiKoihpo6KhKIqipI2KhqIoipI2KhqKoihK2qhoKIqiKGmjoqEoiqKkjWeiDyDTVFZWmtra2ok+DEVRlEnFG2+80WSMqRpqu5wTjdraWrZu3TrRh6EoijKpEJGj6Wyn7ilFURQlbVQ0FEVRlLRR0VAURVHSRkVDURRFSRsVDUVRFCVthhQNEckTkc0i8raI7BSR/2uPf0lETojIW/bjuoQ5d4rIARHZKyJrE8aXi8h2+7V7RUTscb+IPGqPbxKR2oQ560Vkv/1Yn8mTVxRFUYZHOim3QeC9xphOEfECL4vIk/Zr3zTGfD1xYxFZAqwDlgKzgWdFZLExJgrcB9wGvA78AbgWeBK4FWgxxiwUkXXAPcCHRKQcuAtYARjgDRF5whjTMrrTVhRFUUbCkJaGsei0f/Xaj8HWiL0eeMQYEzTGHAYOACtFZBZQYox5zVhrzP4YuCFhzsP2818Aa2wrZC2wwRjTbAvFBiyhURRFGRMON3Xx0v7GiT6MrCWtmIaIuEXkLaAB60t8k/3Sp0XkHRF5UETK7LFq4HjC9Dp7rNp+njzeZ44xJgK0ARWD7Cv5+G4Tka0isrWxUf/ZiqKMnPtfPMSnfroN695WSSYt0TDGRI0xy4AaLKvhPCxX0wJgGXAK+E97c0m1i0HGRzon8fjuN8asMMasqKoasgpeURRlQHpCETqDEU609kz0oWQlw8qeMsa0As8D1xpj6m0xiQE/AFbam9UBcxKm1QAn7fGaFON95oiIBygFmgfZl6IoypgQisYA2F/fOcSWU5N0sqeqRGSa/TwfeB+wx45RONwI7LCfPwGsszOi5gOLgM3GmFNAh4issuMVtwCPJ8xxMqNuBp6z4x5PA9eISJnt/rrGHlMURRkTQhFLNPbWd0zwkWQn6WRPzQIeFhE3lsg8Zoz5nYj8RESWYbmLjgC3AxhjdorIY8AuIALcYWdOAXwKeAjIx8qacrKwHgB+IiIHsCyMdfa+mkXkK8AWe7svG2OaR3G+iqIogxK0RWOfikZKhhQNY8w7wIUpxj82yJy7gbtTjG8FzksxHgA+OMC+HgQeHOo4FUVRMoFjaah7KjVaEa4oipJAPKbR0EEsphlUyahoKIqiJOBYGoFwjOMt3RN8NNmHioaiKEoCoUiMWaV5AOw9rXGNZFQ0FEVREghFYyyaUQzA8Rat1UhGRUNRFCWBUCRGeYEXsAr9lL6oaCiKoiQQisQo9Htwu4SecHToCVMMFQ1FUZQEQpEYPo+LAq+b7pCKRjIqGoqiKAkEo5Zo5Pvc9Kho9ENFQ1EUxcYYQygSw+92UeBTSyMVKhqKoig24ahVzGdZGh4VjRSoaCiKotg41eA+j4t8r4uABsL7oaKhKIpi41SD+9wuCnweujXlth8qGoqiKDZx0fC4ydeYRkpUNBRFUWx6RcMKhGudRn9UNBRFUWxCUUskHNFQS6M/KhqKoig2wYSYRr7XQ0BFox8qGoqiKDaOe8rvcZHvc9EdjmKtPK04qGgoiqLY9I1peIjGTDwNV7FQ0VAURbHpW6fhBtBWIkkMKRoikicim0XkbRHZKSL/1x4vF5ENIrLf/lmWMOdOETkgIntFZG3C+HIR2W6/dq+IiD3uF5FH7fFNIlKbMGe9/Tf2i8j6TJ68oihKIn3rNCzR0GB4X9KxNILAe40x7wKWAdeKyCrgC8BGY8wiYKP9OyKyBFgHLAWuBb4rIm57X/cBtwGL7Me19vitQIsxZiHwTeAee1/lwF3AJcBK4K5EcVIURckkie6pfBWNlAwpGsai0/7Vaz8McD3wsD3+MHCD/fx64BFjTNAYcxg4AKwUkVlAiTHmNWNFln6cNMfZ1y+ANbYVshbYYIxpNsa0ABvoFRpFUZSM4rinvHZFOKCtRJJIK6YhIm4ReQtowPoS3wTMMMacArB/Trc3rwaOJ0yvs8eq7efJ433mGGMiQBtQMci+FEVRMk4wMXvKq5ZGKtISDWNM1BizDKjBshrOG2RzSbWLQcZHOqf3D4rcJiJbRWRrY2PjIIemKIoyMKndU9p/KpFhZU8ZY1qB57FcRPW2ywn7Z4O9WR0wJ2FaDXDSHq9JMd5njoh4gFKgeZB9JR/X/caYFcaYFVVVVcM5JUVRlDipAuGaPdWXdLKnqkRkmv08H3gfsAd4AnCymdYDj9vPnwDW2RlR87EC3pttF1aHiKyy4xW3JM1x9nUz8Jwd93gauEZEyuwA+DX2mKIoSsYJR/v2ngJ1TyXjSWObWcDDdgaUC3jMGPM7EXkNeExEbgWOAR8EMMbsFJHHgF1ABLjDGONc9U8BDwH5wJP2A+AB4CcicgDLwlhn76tZRL4CbLG3+7Ixpnk0J6woijIQqdxT2rSwL0OKhjHmHeDCFONngDUDzLkbuDvF+FagXzzEGBPAFp0Urz0IPDjUcSqKooyWUDSGCHhcosV9A6AV4YqiKDahSAyf24WIxFNu1T3VFxUNRVEUm2Akhs9jfS26XYLP46I7rNlTiahoKIqi2ISiMfye3q/FAp9b3VNJqGgoiqLYOO4phwKvLsSUjIqGoiiKTSjBPQWQr0u+9kNFQ1EUxSalaKil0QcVDUVRFJtQtK9oFHg92kYkCRUNRVEUm+SYhloa/VHRUBRFsUl2TxX4NBCejIqGoiiKTTAaw+dxx3/XQHh/VDQURVFs+qXcqnuqHyoaiqIoNqFItE9xX77WafRDRUNRFMUmOXsq3+ehJxwlFuu39tuURUVDURTFJpV7CiAQUWvDQUVDURTFJlX2FGin20RUNBRFUWz6VYTrmhr9UNFQFEWx6R/T0NX7klHRUBRFAWIxQzhqUsY01D3Vi4qGoigKlpUBJLmnnNX7tP+Uw5BrhCuKouQ63/njAV450ATQbxEm0JhGIkNaGiIyR0T+KCK7RWSniPy9Pf4lETkhIm/Zj+sS5twpIgdEZK+IrE0YXy4i2+3X7hURscf9IvKoPb5JRGoT5qwXkf32Y30mT15RFAXglQNNvHrwDEDK7CmNafSSjqURAT5vjNkmIsXAGyKywX7tm8aYryduLCJLgHXAUmA28KyILDbGRIH7gNuA14E/ANcCTwK3Ai3GmIUisg64B/iQiJQDdwErAGP/7SeMMS2jO21FUZReTrcH4s+Tu9yCxjQSGdLSMMacMsZss593ALuB6kGmXA88YowJGmMOAweAlSIyCygxxrxmjDHAj4EbEuY8bD//BbDGtkLWAhuMMc22UGzAEhpFUZSM0dAe5JL55fg9LmZNy4+Pa8ptf4YVCLfdRhcCm+yhT4vIOyLyoIiU2WPVwPGEaXX2WLX9PHm8zxxjTARoAyoG2ZeiKEpG6AxG6AxGuOqc6bx91zVcsbgq/lqBzwmEq2g4pC0aIlIE/BL4rDGmHcvVtABYBpwC/tPZNMV0M8j4SOckHtttIrJVRLY2NjYOeh5KdlDX0s1tP95KV1CzUpSJpd52Tc0sySPP6+7zWp7XhQj0aPZUnLREQ0S8WILxM2PMrwCMMfXGmKgxJgb8AFhpb14HzEmYXgOctMdrUoz3mSMiHqAUaB5kX30wxtxvjFlhjFlRVVWV/LKShTyzs55ndtVzsLFzog9FmeLUt1miMb3E3+81EdFOt0mkkz0lwAPAbmPMNxLGZyVsdiOww37+BLDOzoiaDywCNhtjTgEdIrLK3uctwOMJc5zMqJuB5+y4x9PANSJSZru/rrHHlEnOrlPtgNW2QVEmkvqOXksjFQW6EFMf0smeugz4GLBdRN6yx74IfFhElmG5i44AtwMYY3aKyGPALqzMqzvszCmATwEPAflYWVNP2uMPAD8RkQNYFsY6e1/NIvIVYIu93ZeNMc0jO1Ulm9h10hKNoIqGMsHUtwcBmDGAaOg64X0ZUjSMMS+TOrbwh0Hm3A3cnWJ8K3BeivEA8MEB9vUg8OBQx6lMHkKRGAcaOuPPFWUiOd0WoNjvodCf+utQ3VN90TYiyrhzsLEz3rJBLQ1lomnoCKSMZzjk+zx0q3sqjoqGMu44rimAoC5uo0wwp9sCA7qmAAq8bs2eSkBFQxl3nCA4qHtKmXjq24MDBsFBA+HJqGgo486uk+3MKrU+pI6bSlEmgljM2O6pgUUj36cxjURUNJRx52RbD4tmFANqaSgTS0t3iHDUMHOwmIZXs6cSUdFQxp2eUJSyAi+ggXBlYnHSbQezNArU0uiDioYy7vSEo5TmW6KhloYykTR1WqJRVTx49pRaGr2oaCjjTjAco9DvwSUqGsrE4ohGZdHAolHgcxOKxoho/A1Q0VDGmWjMEIrGyPO48XlcmnKrTCi9ouEbcBtdiKkvKhrKuBKwP3j5Phd+j1stDWVCaeoM4fe4KBqgGhx6F2JSF5WFioYyrjh3a/ley9LQlFtlImnqDFJZ5MdeeTolzkJMGgy3UNFQxhXnbs3vdeNzuzR7SplQmjpDg7qmoNc9paJhoaKhjCtODCPf68bvVdFQJpamjuCgQXCwsqcAesLaSgRUNJRxpidkiUSebWloTGPsaOsJc/RM10QfRlbjuKcGIx4ID+l7FVQ0lHEmkGhpeFQ0xoq2njA33/cqH/nhpok+lKwlFjM0d4WoGMI91RvTUEsD0luESVEyhhPTyPe5NOV2jIjFDJ/66Rvsb+jEJRCJxvC49f4wmbaeMJGYSd/S0JRbQC0NZZxxUm79Hrem3I4RR8508erBMyyaXkTMwJmu0EQfUlYSr9EYpBocelNuNRBuoaKhjCvxlFufptyOFR0By42yorYcgPr2wEQeTtbS6IhG4RDZU17LIaOiYaGioYwrjqWhgfCxoytoicaCqkKgtymf0pcznZYFlq6loQsxWahoKONKIGyJhKbcjh0dtmicFRcNtTRSkU7fKQCvW3C7JP7eneoMKRoiMkdE/igiu0Vkp4j8vT1eLiIbRGS//bMsYc6dInJARPaKyNqE8eUist1+7V6xyzBFxC8ij9rjm0SkNmHOevtv7BeR9Zk8eWX86YlbGi61NMaITts9Na+iEJdAg4pGSpo6g7hdwjS74/JAiAh5HlfcSp7qpGNpRIDPG2POBVYBd4jIEuALwEZjzCJgo/079mvrgKXAtcB3RcRt7+s+4DZgkf241h6/FWgxxiwEvgncY++rHLgLuARYCdyVKE7K5CPunrIbFqpoZJ4u241Smu+lssiv7qkBONMZorzQh8s1cAsRhzyvLvnqMKRoGGNOGWO22c87gN1ANXA98LC92cPADfbz64FHjDFBY8xh4ACwUkRmASXGmNeMMQb4cdIcZ1+/ANbYVshaYIMxptkY0wJsoFdolElITziK3+PC5RI75VZFI9M4gfAiv4cZJXnUd6ilkciJ1h6u+vrz/OrNE0O6phzyvG51T9kMq07DdhtdCGwCZhhjToElLCIy3d6sGng9YVqdPRa2nyePO3OO2/uKiEgbUJE4nmKOMgkJhKLk2cVSmnI7NnQGI3jdgt/jYkaJnxOtKhqJ7K/v4HBTF9cuncnNy2vSmpPndcULU6c6aQfCRaQI+CXwWWNM+2Cbphgzg4yPdE7isd0mIltFZGtjY+Mgh6ZMNIFwLF5h66TcWoankim6ghEK/R5EhOkleRrTSMJxkf7dmkW8b8mMtObked0E1T0FpCkaIuLFEoyfGWN+ZQ/X2y4n7J8N9ngdMCdheg1w0h6vSTHeZ46IeIBSoHmQffXBGHO/MWaFMWZFVVVVOqekTBA94Sh5Xutt5/dYP9VFlVk6A5H4+hAzivM40xVSiy6B7lBvrVC6aEyjl3SypwR4ANhtjPlGwktPAE4203rg8YTxdXZG1HysgPdm25XVISKr7H3ekjTH2dfNwHN23ONp4BoRKbMD4NfYY8okJRBOdE9Zbz8t8MssHcEE0SixfPZOIZvSm8FXMAzRyNeYRpx0YhqXAR8DtovIW/bYF4GvAY+JyK3AMeCDAMaYnSLyGLALK/PqDmOMI9GfAh4C8oEn7QdYovQTETmAZWGss/fVLCJfAbbY233ZGNM8wnNVsoCecDR+h+dzREPvgjNKVx/RyAOsWo3qafkTeVhZg9P/zLl5SYc8r4uWbm3HAmmIhjHmZVLHFgDWDDDnbuDuFONbgfNSjAewRSfFaw8CDw51nMrkIBi21gcH8LlVNMaCzmCEcrs1RpVd7dygabdxHNEYjqXhV/dUHK0IV8aVVJaGxjQyS2fACoQDTHfcU5p2G6cnHMXjErzD6Pyb73UTVPcUoKKRdZxo7eH7LxzM2YyivoFwSzzU0sgsncEIxbZolORZ1c7tAe2b5NAdisYz+NIlz6sV4Q4qGlnGb948wb89uSdnq3gTA+Ea0xgbOhNiGn6PC49L4k0MFes9OJzMKbA6GKh7ykJFI8to7Aj2+ZlrpBINXYgpc0Rjhu5QNO6eEhEK/R4VjQS6QyMQDa+bQDiasx6A4aCikWU02L7nhhz1QScW9/nV0sg4Tt+p4rzeHJciv4fOoAqzQ094+O6pfJ+bmIFwVEVDRSPLcLJcctXSSPzAxi0NrdPIGJ0JfaccCv1utTQS6BmBpeHc4KiLSkUj63CKsBpyUDTC0RjRmIkHwjXlNvN02uJQ2Ec0PHELRLG++IeTbgu9NR3aSkRFI6swxuS0pdET7ltUpW1EMo8jGkX93FMqGg49I8iecrbXqnAVjayiMxiJf7HmYkwjEEoWDU25zTSOe6o40dLwaSA8EatWaFgNvuPvWXVPqWhkFYnWRS5aGolLvYKm3I4FA7qnNBAex7I0hvfV57hUtVZDRSOrcOIY04v9OdlgzrlL618Rrh/ETBF3T/kT3VNudU8l0B2KUDBCS0NFQ0Ujq3BEY+nsEhragzmXEx5f6jWpNbpaGpkjdfaUFdPItffTSAmEY8NqVgjqnkpERSOLaIyLRinBSIyOHLs7TA6Eq3sq8wzknorGjCYcAJFojFA0NoLsKcc9pddQRSOLaOgI4HO7WDi9yPo9x1qJJIuGxyWI6HoamaQrGMHvccUFGXqtDnVRJbhIR2hpqCtVRSOraOwIUlXsZ3qxP/57LhFM+sCKWOtY6x1w5khcgMnBsTo0gyrhxmWYlka+xjTiqGhkEXHRsNtZ51rabbKlAVaBn7qnMkdiW3QHtTR6ia+lMdKYRkhFQ0Uji2hot0SjqshabS3XLA2n/1HinbDP41ZLI4O09YSZVuDtM1YUtzT0Cy85gy9d4jENfa+qaGQTjZ1Bphf7Kcn34PO4ck402nvCQN9mepZ7Sr/MMkVLd4iyAl+fsUK/9QWp7imrwy2MQDQ86p5yUNHIEmIxQ0t3iPJCHyJCeYGP5q7cWpO4vSeM3+Pq457ye9Q9lUmau0LxpV4d1D3Vi9OVYLiBcJdL8HlcmnKLikbW0B2OYkzvXXih3x2/K8oV2gNhSvL7uk58KhoZpaUrlaWhgXCH7hGsD+6Q53Hpkq+kIRoi8qCINIjIjoSxL4nICRF5y35cl/DanSJyQET2isjahPHlIrLdfu1eERF73C8ij9rjm0SkNmHOehHZbz/WZ+qks5GOgOO6sb5Uc7HJXHtPhNIk0fB7XOonzhCBcJSuUJTywr7XuFAtjTgjTbmF3oWYpjrpWBoPAdemGP+mMWaZ/fgDgIgsAdYBS+053xUR579zH3AbsMh+OPu8FWgxxiwEvgncY++rHLgLuARYCdwlImXDPsNJQnIlby6uttYeCFOS1zezp7TAR1t3brnhJorWbuvGoyzJPVXoc2Ia+oU30kC4M0fdU2mIhjHmRaA5zf1dDzxijAkaYw4DB4CVIjILKDHGvGasXgY/Bm5ImPOw/fwXwBrbClkLbDDGNBtjWoANpBavnKA90HfFtcIctDTaevq7p8oLvDSraGQEJwZWnuSe8rhd5HlduqYGvSmzI7I0PGppwOhiGp8WkXds95VjAVQDxxO2qbPHqu3nyeN95hhjIkAbUDHIvnKSXveUJRpFObhwTntPmJK8vqJRVuijpSs8QUc0OQlHY2w90szJ1p4+4622+CZbGpCb7s6R4FgKw21YCFbarbYRgeFfOYv7gK8Axv75n8BfAZJiWzPIOCOc0wcRuQ3L9cXcuXMHO+6sxflAOzENa4nO3LqraQ/0j2lUFProDEYIRqLx9TWU1Pz09aM88dZJdp9upyMQYc0503ng4xfHX3cstuTsKchNd+dIcALhfs/w75fzvOqeghFaGsaYemNM1BgTA36AFXMAyxqYk7BpDXDSHq9JMd5njoh4gFIsd9hA+0p1PPcbY1YYY1ZUVVWN5JQmnI4cd08ZYyxLI7/vfYpzV6zWxtD8v+f2c6K1hz9912zOmVnM6fa+HQNabPdUcvYU6EJMDoFwlDyvC5cr1T3p4OR53brcKyMUDTtG4XAj4GRWPQGsszOi5mMFvDcbY04BHSKyyo5X3AI8njDHyYy6GXjOjns8DVwjImW2++saeywn6RcI93kIRWKEc6SZX3coSiRm+rmnHP97rtWkZJpwNEZDR5Cbltfw1RvP57zq0n7XrNkW3uSKcLDeV86NyVRmJGtpOKh7ymLIqyciPweuBCpFpA4ro+lKEVmG5S46AtwOYIzZKSKPAbuACHCHMcaR5k9hZWLlA0/aD4AHgJ+IyAEsC2Odva9mEfkKsMXe7svGmHQD8pOOjkAYEUssoG9u/bQUd46TjXY7ZpPsnopbGhoMH5T69gDGwOxSq8VMeaFV/GmMwc5ep6U7REmeB6+7/71god+dkwt7DZeeUGxEQXBQ95TDkKJhjPlwiuEHBtn+buDuFONbgfNSjAeADw6wrweBB4c6xlygPRChyOeJm81FduuHzlwRjR7rLjc5e6qiUC2NdDjVZrmiZk3LByzRCEZidIei8RuMVNXgDoV+D0fOdI/PwWYxgXB0ROm2YGVcafaUVoRnDZ3BSJ+eTM4XQa5UhbfZfadSZU+BWhpD4WRKxS2NFG69lu5QyswpsGJluRQjGyndocioLA0VDRWNrKEjEI5nTkHuVfE6zQqTA+HTbMvjTKeKxmCcbO1vaUBfsW3uCvWr0XAo9HnicbOpTM8oLA2/V7sXgIpG1tAZjFCUYGkU5Vi/oIFiGh63i9J8r1oaQ3CqrYfiPE/8feFYFGcSLY2u0ICuzLJCHz3h6JS/U+4MRuIV8sMl3+smFIkRjU3ttdZVNLKEjkCSe8qXY6IxgHsKeoO6ysCcbA0wuzQ//nvc0ki4bs3doX59p/ptP8XFuakjRJW9MuZwcT6T3TlWdDtcVDSyhI5A32U6e9tZ58adYVtP3zqURMoLfVP+y2woTrX1MGtaXvz38qQEgp5QlEA4NmBMo0xTm4nFDE2dwZGLhi5mBahoZA2WpZEY08ithXPaA2GK/B48KdJBywp88RoDJTWn2gLMSrA0SvI8eFwSFwFHdAeKaZRrESWtPWEiMUNV0UhFw8lonLrXEFQ0soaOpA6wuRgIT+5w61Be6O3jZlH6EghHae4KUZ1gaYiI1bfLFosTdnZV5QBfiI7b6kzX1K3VcFbCrCrOG2LL1OSa9T9SVDSygFAkRjAS6+Oe8ntcuF2SU5ZGco2GQ1lCoZrSn3iNRoKlAZZV4WSdPb+3AbdLuLi2POU+ygstMZnK4twrGqN1T+XGZ3KkqGhkAb3NCntFQ0Qo9Llz5g3alqLDrUN5gY9QNEZXjtSkZJpTthWRGNMAKCvszTp7dlcDK2vLKU3RQgSsrDURaO6euq6Vxk5LfEcqGrpsroWKRhYQ7zuV9KVqtbPOjS/S9p7IgJZGqkwgpZe6Fks0qqf1tTQqCv00d4U4dqabvfUdvG/JjAH34XYJ0/K9NE9h91RD++gsjVxLgx8pKhpZQHvSWhoOhX5PzqT3We6pgWIa/WsOlF52nmyj0OdmTllBn/GyQi/NXSE27K4H4H3nTh90P1N97ZLGjiD5XveI6zTUPWWhopEFJLdFd8iV9ujGWKmOA2X2OF1ZnVYjSl/eOdHG0tml/dp5lxf4aO0J8+s361g8o4h5FYWD7qdiitfDNNrptk6Dx+HiWBodOfCZHA0qGllAPKbh7++eyoW7mvZAhEA4xszS1FkrTqqxs3qh0kskGmP3qXbOqy7t91p5oQ9jYMeJdv7qsvlD7qusYGrXwzR2jLxGA6zW6C5RS0NFIwtIXurVIVdW76u3FwuaUZJaNJwAudMJV+nlQGMngXCM82tK+r3mFPJVT8vnzy+q6fd6MuWFvintAmzsCI64RgPs5BS/Jyc+k6NBRSMLcCyNohx1T522U0YHsjScWEe7Whr92F7XBsD5KSyNmbYI33HVQnxpLF9aXuijZQqnNjeOohrcoThHPpOjYaRrhCsZZKCYRpHfQ1cOBMKdZUlnDFBUle9143aJuqdSsOOEFQSfX1nU77WLa8v52V9fwuqzKtLaV3mhj0jM0BGMDJj+nKsEI1Fau8OjFg1da10tjaygpStEvteN39M3q6MgR9Z1rrctjeklqT+wIkJJnkfdUylwguDuFGtau1zCZQsr017v2uk/NRVTm50iyEyIxlS3NFQ0sgCrO2n/zKIiv5tw1BCMTG4f6un2AGUFXvIGWfymOM+r7qkUHGzo5OyZxRnZ11RObY5Xg48ipgFO7ZSKhjLBDLRMZ6501axvDwwYBHcoyffE3XSKRSAcpT0QYcYAFtpwmcpFlA2jbCHiYCWnTO33qYpGFtDSlXqZzlwpJjrdHhgwCO5QkueNr7mhWDh3x9OHENx0SW6nPpU40WKtj15dlj/EloOj2VNpiIaIPCgiDSKyI2GsXEQ2iMgbickVAAAgAElEQVR++2dZwmt3isgBEdkrImsTxpeLyHb7tXvFrrAREb+IPGqPbxKR2oQ56+2/sV9E1mfqpLON5u4QFSndU7ZoTPJg+Om24IBBcIfiPLU0kmnosGNBo7w7diibwqJR19JDnteV8nM2HDR7Kj1L4yHg2qSxLwAbjTGLgI3274jIEmAdsNSe810RcRzZ9wG3AYvsh7PPW4EWY8xC4JvAPfa+yoG7gEuAlcBdieKUSzR3huJBykRywdIIR2Oc6QoyIx1LQ2MafXB6JU0fYSvvZAp9bnxuF81TsMCvrqWHmrKCEVeDOzjZU1M1bRnSEA1jzItAc9Lw9cDD9vOHgRsSxh8xxgSNMYeBA8BKEZkFlBhjXjPW1f5x0hxnX78A1thWyFpggzGm2RjTAmygv3hNegLhKF2hKBVFqQPhMLn79zd2BDGmt6ZgIEry1T2VTEPcPZUZS0NE4rUaU4261m5qRumaAks0IjFDMBLLwFFNTkYa05hhjDkFYP90OqVVA8cTtquzx6rt58njfeYYYyJAG1AxyL5yCqetQ65aGk6NxszSwb/4ivM8dIWiRKJT98OYTH17AI9LBuzZNRKstUumnjgfb+7JiGhoe/TMB8JT2X5mkPGRzun7R0VuE5GtIrK1sbExrQPNFhz/csrsKd/kf4M6NRpDZk/ZxWbJ52qM4dEtx/jhS4fG5gCzmIaOIJVF/rTrMNKhvHDqtUdvD4Rp6wlTk9QleCTkwo3caBmpaNTbLifsnw32eB0wJ2G7GuCkPV6TYrzPHBHxAKVY7rCB9tUPY8z9xpgVxpgVVVVVIzyliWEw0ciF/v2nh+g75eBUwycGw6Mxw6d+uo1/+uV27v7D7ng7kqlCQ0cwY64pB6tp4dSyNE7Y65GopZEZRioaTwBONtN64PGE8XV2RtR8rID3ZtuF1SEiq+x4xS1Jc5x93Qw8Z8c9ngauEZEyOwB+jT2WUwwmGgV2TGMyi0ZjRxB3Gi4WZ4GmxPbou0+189TO09x0UQ3GwO/eSXnPkLM0tAcyljnlMBXbozuLWCWvRzISinKkdmo0pJNy+3PgNeBsEakTkVuBrwFXi8h+4Gr7d4wxO4HHgF3AU8Adxhjn6n4K+CFWcPwg8KQ9/gBQISIHgH/AzsQyxjQDXwG22I8v22M5xWCi4fe48bplUgfCmzqDVBT6hnSxxDvdJmRQHW+2cus/cVktF9SU8vhbU0s0rFbemcmccigr9NHWEyY8iWJHjovywZcPs+XI8L8C6uwajcwEwif/jdxoGbJhoTHmwwO8tGaA7e8G7k4xvhU4L8V4APjgAPt6EHhwqGOczLR0hXCJtYZzKiZ7g7R01zBI5Z46bn/Y55QV8Gfvms2//n43Bxs7WVDVv3lfrmGlKocyVg3u4NycZKJ533hx5Ew3//TL7YDV3HLTP68ZVsPFupYe8r3ulDdmw0UXYtKK8AnnTFeIaQW+lA3pwAqGT+bivqbOUFpfTo5oJqbd1rX0UJznobTAywcumA3ARntp01ynqTOzNRoO8aaFk6hW44x9Lf5uzSJ6wlGeGKbFebzZSrcdbY0GaCAcVDQmnJYBmhU6TPbV+xrtDKChcCyN9kRLo7k77oeeWZrH7NI8dpxoH5sDzTJ6C/syH9OAyVUV7gTu33fudM6dVcLPNx8b1vzjLZlJtwUVDVDRmHDOdIYGDRJP5tX7YjFrbfB0LI242Z8Y02jpYU5574d9yewSdp2aGqLhrHaY8eypSSkavbVMf7lyDjtPtvNOXWtac7uCEfbVd7Bkdv+VD0eCZk+paEw4Ld0hygoH9s9O5v79bT1hIjGTVjtqj9tFkb93TQ1jDHUt3X0yXpbMLuVQYyc9ockposOhvmNs3FOTsWlhqy0a0wq8XH9hNYU+Nw+8fDituW8eayUaM1xcW56RY3G7hEKfe0qv/aKiMcFYbdEH/lKdzO6pxs7htaO2mhaG43MD4RhzyntFY+nsEmIGdp/OfWtj18k2SvI8GXdPTSuwblAmUyuR5q4wHpdQ5PdQkuflI6vm8du3T3KkqWvIuZuPNOMSWD4vc23rZpTkUd8xtWqGElHRmEBiMUNLd5jyISyNSSsa9t1yOjEN6Nu08HiznVuf4J5aarsYdp7MfdHYdrSVZXPLMloNDlYad7HfM6maFrZ2W0sHOIHsv373fDxuF9974eCQc7ccbubcWSUUZ3B525mleVOu0DQRFY0JpCMQIRozKftOOUzmlcKaRmBpOGZ/XUK6rUP1tHxK873sOtmW4SPNLtoDYfY1dHDR3Gljsv+ySVbg19Idoqyg90t/ekkeN11Uw6/ePDHoqpahSIw3j7dkzDXloKKhTBitPQM3K3Qo8LnpCkUnZSvmxmGulja3ooDdp9uJRGPxwr7EfkEiwpJZJeyyLY3jzd088PLhnItxvH28FWPgorljsxLA5BONMNOSPiNXLK4iFIkNmk23/UQbgXCMS+ZnVjRmleZR3x4gGpt8n8lMoKIxgTiphNMKBndPRSdpK+bGjiA+t4uSvCFrSAF437kzaO0O88bRFo6e6aayyE++r++64ktnl7DndAeRaIwHXj7MV363i+vufYk3jraMxSlMCNuOtiICy8bI0igv8E6qOo3WJEsD4KJ51rXZNsj/3bFIM30dZ5bmE4mZeP3IVENFYwJpiWeFDO6egsmZ4tdop9umW1R1+eIqfG4Xv9p2gqd2nubi2v532kurSwhGYhxs7GL3qXZqyvIJRWJ88Huv8o0N+7LeIkvn+LYda2Hx9OJhVT0Ph/JCP82dk0c0mrvC/azx6cV5zC0vGPRm4VhzN36Pa8hVI4fLLLv55qkp6qJS0ZhA2mxLI/kuKhGnmKh7EtZqNHYEqRxG9k+R38OqBRU8uvU4HYEId1y1sN82S2eXArDjRBu7TrVz+eIqnvrse/jABbO5d+P+rA6S17cHWPVvG3n41SMDbmOM4e26Vi4cIysD7Pbok8TSMMbEA+HJrJhXxhvHWgYU4mPN3cwpL8h4MoGz3r2KhjLupGdpOKv3TT5Lo6kzRFWKFQkH4+pzrfW8rl4yg/OqS/u9flZlIX6Pi2d319MRiLDEzoz50p8txe0S/rD9VEaOfSx4Zudp6tuDfOm3O3l2V+p2KO2BCK3d4THtrzV7Wj6BcIxTbT1j9jcyRWcwQiRmUt5YXTSvjMaOYDzTLpljzT3MKx99Z9tkZtmicXoSXL+xQEVjDDjU2MlrB88M6YpotS2NgZoVQkLbgknYfyrdZoWJvP/8Waw+q4J/XHt2ytc9bhfnzCzmWbsH1bmzrDTc8kIfly6o4PfbT2Wti+rZ3Q3MLS9g6ewSvvCr7SmP01n7oTpDbS9S4WQTbT6c/U2jW+Nxv/43H07txdaj/c/DGGO1oRkD0Sgv9OFzuzjVrpaGkiG+8MvtfPgHr7P2Wy/Gs4BS0dodoiTPM2CzQugVjclmaVhdWoff2ruyyM/Pb1vF4hnFA26zZHYp4ahBBM6Z2bvdBy6YxdEz3VnpouoMRnjt4BmuWTKD9atraeoMsq++s992J1pt0Zg2dqJx7qwSivyeSSEaTpZXqgzDxTOK8bld7D3d0e+1lu4wncEIc8dANERkSqfdqmhkGGMMu0+1s2JeGceau/n2xv0DbtvSHU7pq01ksq7ed7otgDFQMwZffk6RX21FYVxUAa5ZMhNPlrqoXt7fSCgaY825M1h1VgUArx8602+7E3Z9yuwxFA23S1g+r2xEa1OMN44LN1UBrNslVJflU9fa3010zL5ZGwvRACuuoTENJSOcbAvQEYxww4XVrLt4Lr9588SA1kZrT//882Qma1dN5455LL78HNFYMqtvE7qyQh9nzyzOSkvjj3saKcnzsKK2jDnlBVRPy+e1g/1F42RbAL/HReUwY0HDZeX8cvbVd2Z9O5HB3FNgLazkrMyXyNEzVouRuRVjIxqz1NJQMsVeuy/S2TOLuf2KsxCB77+Yut1Ba3eIaYPEMwAKfU4gfHJlT52Mi0Zm0x0BzplZEv8CTmZBVREHG/u7fSaat+tauWheGV639ZFbvaCCTYfPEEsqEDvR0kP1tMys/TAYK+2Ct2y3NhI73Kaipiw/bp0l4tyoZWKJ11Q47qlAeHJ9LjOBikaG2Xva+sJaPKOYWaX5XLN0Jht3N6TcNrk9Qiomq6VxcgwtjXyfm5f+13u5ZXVtv9cWTi/iRGtPVlWJByNRDjR09rGMVp1VQUt3mL31ff3xda09Y+qacrigphSfx8Xrh7JcNLpCyCArW9aUFdDUGer3/z7W3E1Vcf/i0EyxoLKIUDTG8q9s4Ndv1o3J38hWVDQyzN7T7cwqzYu/yc+eUcyptkDKL7HWFO0RkvG6Xfg8rkknGidae6gs8pHnHZsPbWmBN2UCwYKqIoyBQ03ZY23sr+8kEjN91nRYvaACl8DnHn2rz92+Y2mMNX6Pm5W15bxyoGnM/9ZoaOkOU5KX+n8Nvet+n2jta20ca+4es3gGwM3La/jxX61kbkUh//nMvn4WYy6jopFh9tZ3cnZCRk9tZSEAR5v7tnGORGN0BCKDthBxmIxNC0+0BsbljjmZBdOt632gIXtEw1k4yilMBCs76r6PLqcjEOETP9pCKBIjEI7S1Bkc03TbRN6zqJK99R3xBZ+ykeYhrHFHNI4nxDWMMRxu6hpT0XC5hMsXV/HJK86irqUnZVJDrjIq0RCRIyKyXUTeEpGt9li5iGwQkf32z7KE7e8UkQMisldE1iaML7f3c0BE7hXboSsifhF51B7fJCK1ozneseShVw5z78b9HGzo5OyEdNH5FdaXWHLv/7Yepxp86IBnod896UTjZGsPs0vHXzRqKwpxCRxsHHqthfFi18l2CnzufoVma5fO5IvXnUtnMMLuU+3xbJzxsDQA3rOoCoCX9o+9tfH6oTN87IFNg3alTcXptkC8AjsVTkPLxGD4/oZO6tuDGe9um4q1S2dSkufh0a3Hx/xvZQuZsDSuMsYsM8assH//ArDRGLMI2Gj/jogsAdYBS4Frge+KiOO7uA+4DVhkP661x28FWowxC4FvAvdk4HgzTlNnkK/+YQ/f2LCPUDSWZGlYb+pDSaKRTrNCh5kleZxqzd67wWSMMZabZZzumBPJ87qZU16QVcHwXafaOXdWScp2FvHGe8da4oV942WhnTOzmMoiHy/vbxzzv/Xd5w/y0v4m9qeoTRmMU0PEeKqK/PjcrngrfYDn9lgxxCvPrhrZwQ6DPK+bP1s2myd3nO6zVHEuMxbuqeuBh+3nDwM3JIw/YowJGmMOAweAlSIyCygxxrxmrBLZHyfNcfb1C2CNjHVayQh4bOtxQtEY/7j2bJbNmcalCyrjrxXneaks8vezNFrTaCHiUFtRyJEz2XPnPBSt3WF6wtEJcU+BnUGVJe4pYwy7T7b3Sw92mFWaz6zSPLYda4375WvGSWxdLuHdCyt5+UDTmPrkT7X18JItTLuHscZ7JBrjdHtgUMvL5dRqJFgaz+1p4JyZxeP2/lu7dCahSIy3j+f2Oi8OoxUNAzwjIm+IyG322AxjzCkA++d0e7waSLTh6uyxavt58nifOcaYCNAGVCQfhIjcJiJbRWRrY+PY3zUlEo0Z/nvTMVafVcEdVy3kN3dc1s+cnl9ZwJGmvoG61jSaFTrUVhbS0BHMumD4zpNt/Pi1I/1y/XurmjOfbpsOC6cXcaipKyvWOzje3ENHMNInCJ7MRXPL2Ha0hdcOniHP6xrUHZNp3rOoiqbO0JguofurbScwBjwuYU+K6u2BqO8IEjOWsA5GYq1Gm91a/73nTB90Tia5oNqyFt+uax23vzmRjFY0LjPGXAS8H7hDRC4fZNtUFoIZZHywOX0HjLnfGLPCGLOiqmrsTdJEXjt4hrqWHj66at6A29RWFKZwTw29AFPifCDrrI0v/noH/+fxnaz6t438T4JPdyzTbdNh4fQiQpEYh9NYQ3qseeOYlRl1QU3/5osOF80r40RrD4+/fZL1l9bGaznGg3cvsqziTMY1ojHDl57YyYv7GukORXhkyzFWzi9nyeySlC0/BiLdWp/EWo2XDjQSjZlxFY3SAi+1FQW8o6IxNMaYk/bPBuDXwEqg3nY5Yf90ihTqgDkJ02uAk/Z4TYrxPnNExAOUAlmVWL758BlcMrj/dH5VIU2dwT4+z3izwrQsDSsukmytTCQ7T7bx9vFWbn33fFbUlvG/fvkOv3zDMhjHo3/SYKywG9llQ2+lTYeaKcnzcM7MwSwN60612O/hb6/o3w5+LJlRksfZM4p5OYOi8dyeBh569Qi3/WQrH//RFk609PCZ9y7knJnF7BmGRXMyzffRvIpCmjpDNHeFeP3QGQp9bpbNGbvW8qm4oGYa79Spe2pQRKRQRIqd58A1wA7gCWC9vdl64HH7+RPAOjsjaj5WwHuz7cLqEJFVdrzilqQ5zr5uBp4zWdbC9K26NhbPKO7TAykZJ4Pq6JneL/3WnhAel1A8yDyHbLQ0Htl8HJ/HxWfeu5AH1l/MJfPL+effbCcQjrLrZDvTCryUD9FXa6yYX1nI9GJ/VqRBbj7czMr55YM2pVw6u5Tqafl8/pqz07qJyDTvWVTJ5iPNGSuI/NErh5lZkseMkjw2H27mX/5kCe9ZVMU5M0to6gzFlwEeipN28sesIUTDWaxr06EzbD3SwkXzyvCMo7UGliV5qi1AQ8fkSVgZKaO5sjOAl0XkbWAz8HtjzFPA14CrRWQ/cLX9O8aYncBjwC7gKeAOY4zzLv0U8EOs4PhB4El7/AGgQkQOAP+AnYmVLRhjePv40AvmzK+yvvQTM3qONfcwoyQvrXYRhX4PVcX9g+kTRV1LN7958wR/cv4sphVYBXy3X76AQDjGG0dbeP3wGS6ZXz7mrTAGQkRYdVYFrx8auj39WNLQHuBQUxeXzO8XhuuDz+Pi5X+6ivWX1o7PgSXxHnu97c0ZaCmy53Q7rx48w/pLa3nktlV89yMX8YnLagE4Z1ZxfJt0ONnaQ0meJ960cyAuqJlGvtfNM7vq2VvfMS6ptqmOAWD7FLA2RiwaxphDxph32Y+lxpi77fEzxpg1xphF9s/mhDl3G2MWGGPONsY8mTC+1Rhznv3apx1rwhgTMMZ80Biz0Biz0hhzaDQnm2kON3XR1hMe0hReWFVEgc/dZ2nK/fUdLJ6R/kI787Mkg2rToTP8yb0vA/A37zkrPr5yfjkel/DY1uMcb+6Jd3KdKFadVUFDR3BC4xqbbPeY0+dpMCYyKXBlbTk+j4vndqdeGGo43P/iIfK8Lj68cg6zSvO57vxZ8XNzXHR7TqUX1ziZZksVr9vFitoynnj7JMaQsifZWHNedQkugbeP535cQyvCR8Fb9htk2ZzB36Qet4vl88riPvZINMahxi4Wzxx4zYhkaisLOJwFMY2vP7OXIr+H337m3X0yggr9Hi6cO40n3rbCURMvGtYX9UT2Vtp8uJlCnzvelTdbyfe5ue68mfzPG3XxVPCRcKixk9+8eYKPrZqXMpW8vNDHnPJ8fr7lWLy4dTBOtKbfUmXVWRVEYwaPS8Y9ngFQ4PNw4dwyvvfiIX76+tGsXQgsE6hojIK3jrdS6HOzcPrQFsMl88vZc7qD1u4QR850E4rGWDx9OKJhBdNHWhl+sLFz1B05A+Eobx9v47rzZ8bboyRy6YJKjLEKFs8eZBGl8cCJa7w2QXGNcDTGUztPs3pB5bj710fCp65cSHcoykODrF8+FPdu3I/f4+b2KxYMuM3Xb34Xx5u7+fR/bxuyNuRUW/qtaFYvsG5SllaXUuAbOk44Ftz/seWsOquCf/nNDj798zdzttgv+9/NWcLRM138x9N7eHFfI8fOdPPcnnqe2VnPBTXTBg1yOqy0/dpbjrSw3+5sOtjqdMmcVWkJ03CKoxxOtvZw7bde5Acvjs67905dG6FobECf8WULrfTNS+aXp6x+Hk9ErMK1V8a4cG0gNu6up7EjyIdXzhl64yzg7JnFXL1kBj965Uh8tbx0icUM927cz2/eOsktl86jsmjgJX4vOauCL153Li/tb+LNQVw5ncEIbT3htEXj/OpSKot8XL6ocuiNx4iKIj8Pffxi/te1Z/PUjtN89pG3JuxYxhIVjTQwxnDnr7bznT8e5JYHN3P5f/yRv3poKx638Nn3LUprH04r6s2Hz7C3vgMR0rJQHFYvqMDjEp7dNXy/8yNbjhOOGrYmxFRGgtONdSDRWDZnGu+qKeX6ZdUpXx9vLl9cRXNXiB0nxz84+bNNx5hdmseVZ49fvcBo+dz7FtMTjnLHz7YRjsbSnvdffzzANzbs48YLq/nc+xYPub3z/hgsu81ZRCnd9Vi8bhcbPncFn3lvep/HscLlEv72yoX87ZUL+OPehpxcqElFIw027m7g1YNn+OJ15/DQJy7mPz/4Lr730eU89/kruSRN332e18odf2l/E/vqO5hTVjCsXv+l+V5WL6jg6Z2nh+UvjURjPLrlGADv1LWOyte6+XAzi6YXDbhErc/j4vFPv5vrzp814r+RScaicC0dDjd18dL+JtatnJuWFZotLJldwr/deD6vHTrDZ/77TRrS7H777O56Lq4t4xt/8a60WuGXF/o4Z2Yxrx4c+P/y620ncLtkyMyzRMoKffg82fGVdtNFNcQM/PrNExN9KBknO65wFhOLGf7tyd2cVVXIJy6bz5VnT+em5TVce97MYb9Bb7qomj2nO3h6Z/2wMqcc1i6dyZEz3ewfRl+l5/Y0UN8e5IrFVbR0hzne3H9pzHSIxgzbjrZwcRqZQNlCZZGf86pLeGHf+LaW+faz+8jzulg3SVxTidy0vIY7338Oz+1p4KqvP883ntk7aNC6J2TV5VxcO7wU69ULKth6pCVl19uuYIRHtx7n/efNHNeWKpmktrKQFfPK+OW2upwLiqtoDMHz+xo42NjF369ZNOr2Dn+xYg5rl84gGjMsGkGg+JolMxCBp3ecTnvOf28+xvRiP5+72nIbvDXCVgfbT7TREYywcgJy4EfD5Yuq2Ha0hT/uaRiWy2Wk7DndzuNvn+QTl81nevHk/MK7/YoFPP25y7l8cRX3PneAq7/xwoDC+05dK5GYYfm84aW5XrqgkmAkxpvH+r4fw9EYD716hI5AhI9PUN1Kpvjzi2o40NDJ2zlWu6GiMQQPv3qU6cX+jLhcRIR/v+ldvO/cGaxdOnPY86eX5LFszjSe25t6+dhk6lq6eWFfIx+6eA5LZ5fg97h4Z4R55I+/dQKf28VVk8hHD3DjhdUU+j184qEtfO7RsQ9MfmvDfor8Hj55+cAZRJOB+ZWF3PfR5Tzx6cuYVuBl/YObUwrHG8esONmFc4cnGivnl+MS+qwceKChk1Vf3ch/PL2X5fPKhi1E2cafvmsWBT43P3v96EQfSkZR0RiEw01dvLCvkY9cMi9jTeRKC7z8cP2KEeeSX7agknfq2tJKvX10i9VE8EMXz8HrdrF0dsmIOnFGojF++/ZJ3nvO9AlpczEaFs0oZvM/r+Fjq+bx5I7TafvpR0JzV4hnd9fzlyvnTrrrNBAX1EzjiU+/m/mVhXz5tzv7WWvbjrZyVmXhsFvGlOZ7uXRBJfe/eChev3T/iwfpDkX57kcu4sd/tXJCCx4zQXGelxsurOaJt0+mrH8JR2PxDtG7T7Xz3ecPpN1iJRV1Ld0ZawUzGCoag/DoluN4XJJVaZOrF1hFTFuGaMYXCEd5dMtxrlxcFV/d7F1zrKZqw+k0CvDSgSaaOkPceFF2ZEUNF7/HzccvqyUaM/zmrbELTP5++ykiMcMNF07O6zQQeV43/3zduRxs7OKnCXfNxhi2HWsZtpXh8O11y6guy+fWh7awcXc9v3nrJDctr+a682cN2sttMvHRS+YRjMT4xRt1fcYj0Rgf+eEmrv7mC7QHwnzhV9v596f2ctnXnuPCLz/DZV97jodfPUJ7IMyOE23c8J1X+M4fDwz6t7746x3c8J1XxjyGoqIxAMYYntpxitULKphekj2+6eXzyvC5XUMWrX35d7to6AjyN5f3tvr4+KW1TCvwsu7+14YlHL958wSl+d5xWQltrFhQVcRFc6fxizfGLjD5xFsnWDyjiHOGUek/WVhz7nRWn1XB9184RDRmCEaifPl3u2juCrFy/shEo6LIz09vvYSyQh+3PryVUCTGxy+dn+Ejn1iWzC7h4toyfvDSoT4JBV9/Zh+bDzfT1Bnijp9t4+3jrdxx1QLWXzqPD1wwm+pp+dz1xE4u+NIz/Nl/vcw7da1869l9A7bFOd7czUv7G7n2vJljbqHlhpxngMaOIP/86+383ZpFnFddyr76To6c6e7zpZsN5HndXDh32qDpio9uOcZ/bzrGJ69Y0GcVwXkVhTx622puuu9VvvK7Xfz0ry8Z8u91BiM8vfM0N11Ug9+TfopwNnLz8jl88dfbeWZXfb+YUjRmRpUeW9fSzZYjLfzj2rMnvVslFSLCx1bP429/to2XDzTx6JZj/GH7aW5ZPY8bL6wZegcDMHtaPo/dvppbHtzEwulFw6pdmiz87w8s4YbvvMJXf7+be26+gFcPNvG9Fw7y4ZVzaQ+E+f07p6gs8vOZ9y6KpywbY9hypIUtdvfh65fN5obvvMKdv3qHS+ZXUJxnNTH93guHqJ6Wx1lVRQiWK3qsUdGw8XlcvHm8lX/65Ts8fsdlPLXjNCJw9ZIZE31o/Vi9oIJvb9zPNzfsY0VtGe9ZZFkAsZjhnqf38P0XDnHZwgo+f03/QqvaykL++j1ncc9Te9hxoo3zqgdeHAisTK1AOMaNOeByufHCan6++RiffeQtHvrExaycX07MwA9eOsS3n93Pl69fygdXjOxD9+1n9+NxCdcvm53ho84e1pw7nWkFXv71d7vY39DJZ9+3iM+mUcw3FDNL83j6s5eTBQstjgkX1EzjtssX8L0XDlJdls//vHGc+ZWF/J8PLApSOtsAAApISURBVOFUWw8bdtVz++Vn9alxERFWzi/v0+zy79Ys4t+e3MOmw804xvLs0jyrS8TuBtacM31cFj6TXMshXrFihdm6deuI5j614xSf/Ok2blk9j5f2N1FZ5ON/Pnlpho9w9Ow40caf/dfLxAyU5Hl49c41FPk93P37XfzgpcN85JK53PWnSwesI2nrCXPZ155jzbnT+fa6Cwf9Wx/94SaONnfx4j9elRN30A0dAf78u69S19JDZZGfnlCErlCU0nwv0Zjhmc9dnvYHLxCO8uaxVoKRKB//0RZuu/wsvnjduWN8BhPLXY/v4OHXjlI9LZ+Nn78irWI+xXqvfP6xt/n99lOIwGO3r453VmjuClFW4B3y82WM4UBDJzVlBbR0hzjc1MXK+eV8/4WDfP2ZfTz8Vyu5YvHIXcgi8oYxZsWQ26lo9OX2n2zl6Z1Wq457bjqfD108N1OHllG6QxF2n+rgpvte5QvvP4dozPAfT+9l/ep5fOnPlg75BvzqH3bzg5cOcdvlZ/HZNYtTVqefbguw+msb+cxVC/mHa84eq1MZd5q7Qjy14zRvHG2hJN/Dytpyls4uZe23XmRFbRk/+vjFQzYZPHamm9t/+ka8F1hVsZ/nPn8FxXm5kTU1ELtOtnP9d17mWx+6kD+5IDsq/ycLxhh+v/0UwXCMm5aP3KWXisaOIFXFA/f8SgcVjRESCEc51NhFdVk+pfnZ/wXwsQc2selQM6FojOuXzeYbf7EsLd98dyjCl3+7i0e2HKc038u6i+fwj2vP7vNl+S+/2c4jm4/z3OevZG5FwVieRlbw6JZj/NMvt/PhlXP56o3nDSi8z+6q5x8es2o+/uVPltDSHWJFbRnL502uwseR0h2KTFgnWWXsSFc09D+fRJ7X3WediGznM+9dxCsHXuNjqywLI91gboHPw9duuoCbltfwo1cO8/0XD1FTls/HVtcCViv1n28+zkcumTslBAPgQxfP5eiZbr77/EEONXbyicvmc/WSGfFrGo7G+Poze/n+C4c4r7qE7/7l8ilzbRJRwZjaqKWRA7R2h1IuepMuxhj+8geb2HO6nef/v6sIRKL87c+2sfd0B8//45WDtrrONYwxPPjKER58+XB8EaDFM4qIxAx1LT0cburio6vm8i9/skT9+UpOoe4pZVjsOtnOB/7fS5QXWsHhqDHcc9MFWdPmfLyJRGNs2FXPo1uP09wVQkSYlu/l5uU1/Om7cjdDSpm6qGgow+bXb9bx0r4mvG4Xn7xyAfNTrM6nKEpuklMxDRG5Fvg24AZ+aIz52gQfUk5y44U1oyrUUhQl98n6NiIi4ga+A7wfWAJ8WESWTOxRKYqiTE2yXjSAlcABY8whY0wIeAS4foKPSVEUZUoyGUSjGjie8HudPRZHRG4Tka0isrWxcXxXaVMURZlKTAbRSFV40Cd6b4y53xizwhizoqpq8nZiVRRFyXYmg2jUAYld5GqAkxN0LIqiKFOaySAaW4BFIjJfRHzAOuCJCT4mRVGUKUnWp9waYyIi8mngaayU2weNMTsn+LAURVGmJFkvGgDGmD8Af5jo41AURZnq5FxFuIh0AHvT2LQUaMvQn60EBl5Kb/hk8tgmw/70+o0evYajI5PXL9vPNdX+KoFCY8zQmUTGmJx6AFvT3O7+8f6bw9hfxo5tkuxPr59ew4neX8au3yQ41377G875T4ZA+Fjx24k+gEHI9LFl+/4yTbafb7ZfP8j+c87ma5jt5zqq/eWie2qrSaPp1mT/m7mEXr/Ro9dwdEz16zec889FS+P+KfI3cwm9fqNHr+HomOrXL+3zzzlLQ1EURRk7ctHSUBRFUcYIFY0UiMgcEfmjiOwWkZ0i8vf2eLmIbBCR/fbPMnv8ahF5Q0S22z/fa48XiMjvRWSPvZ8psQ5Ipq6f/dpTIvK2vZ/v2a3yc55MXsOEfT4hIjvG+1wmggy/B58Xkb0i8pb9mD5R55UVZDKVK1cewCzgIvt5MbAPay2Pfwe+YI9/AbjHfn4hMNt+fh5wwn5eAFxlP/cBLwHvn+jzmyzXz/69xP4pwC+BdRN9fpPtGtpjfw78N7Bjos9tsl0/4HlgxUSfU7Y8JvwAJsMDeBy4GqtocJY9NgvYm2JbAc4A/hSvfRv4m4k+n8l4/QAvVqrghyb6fCbbNQSKgJftL80pIRoZvn4qGgkPdU8NgYjUYt2FbAJmGGNOAdg/U5mpNwFvGmOCSfuZBvwpsHEsjzfbyMT1E5GngQagA/jFGB9y1pGBa/gV4D+B7jE/2CwkQ5/hH9muqf8tIqmWa5gyqGgMgogUYblEPmuMaU9j+6XAPcDtSeMe4OfAvcaYQ2NxrNlIpq6fMWYt1l2hH+jnq89lRnsNRWQZsNAY8+sxPdAsJUPvwY8YY84H3mM/PjYWxzpZUNEYABHxYr3ZfmaM+ZU9XC8is+zXZ2Hd/Trb1wC/Bm4xxhxM2t39wH5jzLfG/sizgwxfP4wxAayW+FNmqd8MXcPVwHIROYLlolosIs+PzxlMLJl6DxpjTtg/O7DiQivH5wyyExWNFNjm5wPAbmPMNxJeegJYbz9fj+UndVxPvwfuNMa8krSvf8VqEPbZsT7ubCFT109EihI+4B7gOmDP2J/BxJOpa2iMuc8YM9sYUwu8G9hnjLly7M9gYsnge9AjIpX2cy/wAWBKZKANhBb3pUBE3o2V6bQdiNnDX8TyiT4GzAWOAR80xjSLyL8AdwL7E3ZzDVbG1HGsLzrHP/pfxpgfjvlJTCAZvH4C/A7LLeUGngM+Z4yJjMd5TCSZuobGmMQ76Vrgd8aY88b8BCaYDL4Hu4AXsRIx3MCzwD8YY6LjcR7ZiIqGoiiKkjbqnlIURVHSRkVDURRFSRsVDUVRFCVtVDQURVGUtFHRUBRFUdJGRUNRxhkR+aSI3DKM7WunSndaJfvxTPQBKMpUQkQ8xpjvTfRxKMpIUdFQlGFiF8n9/+3dMUhVcRTH8e8hXYLIpdVVhAjCSUSa3V1C0mwVXBud3Cyit7waFKsl2poKJ0GiPcRJlwYXQS1tPQ33/+DxULvvoiL0/YyHy+H9l3s4f3i/+4Xqj2IPqWK3Z4FR4CVVquwB8DQz90tsxzdgAvgcEXeAk8xcKdlQbaoY/V3gWWYeRsQYsEoVMrh1faeTLub1lNTMCPA2Mx8Av4AFoAVMZ2bnhb/c9fxQZj7KzBc9fd4Bz0ufH8BSqa8Bi5k5fpWHkPrlpiE187Mro+gDVUTFfWCjJGffAva7nv/Y2yAi7lINk81SWgc+nVF/D0xd/hGk/jk0pGZ683d+A9sXbAanffSOM/pLN4LXU1IzwxHRGRCPge/AvU4tIgbLtxnOlZnHwGFETJbSE2AzM4+A4xK6BzBz+T9fasZNQ2pmB5iLiDdUyagt4CvwulwvDQCvgO1/9JkD2hFxG9gD5kt9HliNiD+lr3QjmHIr9el/ihiXenk9JUmqzU1DklSbm4YkqTaHhiSpNoeGJKk2h4YkqTaHhiSpNoeGJKm2v2Txl8IJN51AAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_data['inc'][-200:].plot()" ] @@ -252,10 +2298,8 @@ }, { "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, + "execution_count": 10, + "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 8, 1), 'W')\n", @@ -274,7 +2318,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "metadata": {}, "outputs": [], "source": [ @@ -298,9 +2342,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 12, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAD8CAYAAABU4IIeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAG35JREFUeJzt3X+Q1HV+5/HnCweHjYsKCC4/1OFqORN07zROoVveD3HDj81tqXvn7rEa5SpWYVy3ypypU6hw5Z6Si1xdrTnirdFa3bjrD/Q2S0likCBCnZcQYCg1ikoGbwkSWWfIoGCq4HbkfX/0Z7Snne7p6f729LdnXo+qru7+9Ofz6c98GPo9n1/fVkRgZmaWlQnNboCZmY0tDixmZpYpBxYzM8uUA4uZmWXKgcXMzDLlwGJmZplyYDEzs0w5sJiZWaYcWMzMLFNtzW7AaDrnnHOio6Oj2c0wM2spe/bsORIR06vNP64CS0dHB11dXc1uhplZS5H0dyPJ76kwMzPLlAOLmZllyoHFzMwy5cBiZmaZcmAxM7NMObDkQM+xE3zz4R30HD/R7KaYmdXNgSUH1m3tZveBPta92N3sppiZ1W1cnWPJmwtXb+Jk/6lPnj+x8yBP7DxIe9sE9q35ahNbZmZWO49YmujluxZyzSWzmDSx8M8waeIErr1kFi/fvbDJLTMzq50DSxPNOHMSk9vbONl/iva2CZzsP8Xk9jZmTJ7U7KaZmdXMU2FNduSjk9x4+QXcsOB8ntp1kF4v4JtZi1NENLsNo6azszN8rTAzs5GRtCciOqvN76kwMzPLlAOLmZllyoHFzMwy5cDSAnwy38xaiQNLC/DJfDNrJd5unGM+mW9mrcgjlhzzyXwza0UOLDnmk/lm1oocWKrQzMXzgZP5G759JTdefgG9H50c9TaYmY2ET95XYfWG13ly10FuXHA+a77+pQa0zMwsv0Z68t6L9xV48dzMbOQ8FVaBF8/NzEauqsAi6YCk1yW9KqkrpU2VtEVSd7qfUpR/laT9kvZJWlKUflmqZ7+kdZKU0tslPZPSd0rqKCqzPL1Ht6TlRelzU97uVPb0+rtjMC+em5mN3EhGLAsj4pKiebaVwNaImAdsTc+RNB9YBlwELAW+L+m0VOYhYAUwL92WpvRbgKMR8UXgAWBtqmsqcA9wObAAuKcogK0FHkjvfzTVkTkvnpuZjUxVi/eSDgCdEXGkKG0fcFVEHJY0E9geERdKWgUQEb+f8m0GvgscALZFxC+n9G+l8rcO5ImIHZLagJ8D0ykEqKsi4tZU5mFgO7Ae6AW+EBH9kr6cyn8yOhqKL5tvZjZyjbpsfgB/IWmPpBUp7dyIOAyQ7mek9NnAu0VlD6W02elxafqgMhHRD3wITKtQ1zTgg5S3tC4zM2uianeFXRkR70maAWyR9HaFvBoiLSqk11KmUl2DG1MIhCsAzj///KGymJlZhqoasUTEe+m+B9hAYb3j/TQFRrrvSdkPAecVFZ8DvJfS5wyRPqhMmgo7C+irUNcR4OyUt7Su0rY/EhGdEdE5ffr0an5cMzOrw7CBRdIZkiYPPAYWA28AG4GBXVrLgefS443AsrTTay6FRfpdabrsuKQr0m6wm0vKDNR1PfBSFBZ/NgOLJU1Ji/aLgc3ptW0pb+n7m5lZE1UzFXYusCHtDG4DnoqIFyTtBp6VdAtwEPgGQETslfQs8CbQD9weER+num4D/hj4HLAp3QAeBX4saT+FkcqyVFefpPuA3SnfvRHRlx7fDayXtAZ4JdVhZmZN5ku6mJlZRY3aFWZmZlYVBxYzM8uUA4uZmWXKgcXMzDLlwGJmZplyYDEzs0w5sJiZWaYcWMzMLFMOLGZmlikHFjMzy5QDi5mZZcqBxczMMuXAYmZmmXJgMTOzTDmwmJlZphxYzMwsUw4sZmaWKQcWMzPLlAOLmZllyoHFzMwy5cBiZmaZcmAxM7NMObCYmVmmHFjMzCxTDixmZpYpBxYzM8uUA4uZNVXPsRN88+Ed9Bw/0eymWEYcWMysqdZt7Wb3gT7Wvdjd7KZYRtqa3QAzG58uXL2Jk/2nPnn+xM6DPLHzIO1tE9i35qtNbJnVyyMWM6uoUVNVL9+1kGsumcWkiYWPoUkTJ3DtJbN4+e6Fmb6PjT4HFjOrqFFTVTPOnMTk9jZO9p+ivW0CJ/tPMbm9jRmTJ2X6Pjb6PBVmZkMajamqIx+d5MbLL+CGBefz1K6D9HoBf0xQRDS7DaOms7Mzurq6mt0Ms5bQc+wEa/78Lf5i78858YtTTJo4gSUXfYHf/Te/4lHFOCNpT0R0Vpu/6qkwSadJekXSn6XnUyVtkdSd7qcU5V0lab+kfZKWFKVfJun19No6SUrp7ZKeSek7JXUUlVme3qNb0vKi9Lkpb3cqe3q1P4uZDc9TVVarkayx3AG8VfR8JbA1IuYBW9NzJM0HlgEXAUuB70s6LZV5CFgBzEu3pSn9FuBoRHwReABYm+qaCtwDXA4sAO4pCmBrgQfS+x9NdZhZhgamqjZ8+0puvPwCej862ewmWQuoaipM0hzgceD3gDsj4muS9gFXRcRhSTOB7RFxoaRVABHx+6nsZuC7wAFgW0T8ckr/Vip/60CeiNghqQ34OTCdQoC6KiJuTWUeBrYD64Fe4AsR0S/py6n8J6OjoXgqzMxs5Bo1FfYHwF3AqaK0cyPiMEC6n5HSZwPvFuU7lNJmp8el6YPKREQ/8CEwrUJd04APUt7SuszMrImGDSySvgb0RMSeKuvUEGlRIb2WMpXqGtwYaYWkLkldvb29Q2WxFudLgpjlSzUjliuBayQdoDAFdbWkJ4D30xQY6b4n5T8EnFdUfg7wXkqfM0T6oDJpKuwsoK9CXUeAs1Pe0roGiYhHIqIzIjqnT59exY9rrcaXBDHLl2EDS0Ssiog5EdFBYc3jpYj4DWAjMLBLaznwXHq8EViWdnrNpbBIvytNlx2XdEXaDXZzSZmBuq5P7xHAZmCxpClp0X4xsDm9ti3lLX1/GycuXL2JjpXP88TOg0QUzll0rHyeC1dvanbTzMa1ek7e3w8sktQNLErPiYi9wLPAm8ALwO0R8XEqcxvwA2A/8A4w8AnwKDBN0n7gTtIOs4joA+4DdqfbvSkN4G7gzlRmWqrDxhFfEsQsn0Z08j4itlPYlUVE/APwlTL5fo/CDrLS9C7g4iHSTwDfKFPXY8BjQ6T/XwpbkG2cqvacRc+xE3zn6Vd48IZLc3kGI+/tMxspXyvMWlo15yzyvgaT9/aZjZQv6WJjVum1rgbk5bLseW+f2YCGXdLFrNXkfQ0m7+0zq5UDi41Zeb/WVd7bZ1YrXzbfxrS8X5Y97+0zq4XXWMzMxrAsdh16jcXMxhVf0qeyZuw69FSYmbW04g/ONV//UrObkxuj8Q2g5XgqzJrOBwStFt6uXVmW3wDqqTBrOT4gaLXwdu3Kmrnr0FNh1jTNHKqPJo/IGsPbtYfXrF2HHrFYw5VbXB0vf3F6RNY4/urkyh6+qZM1113M/Flnsua6i3n4pqpns+riEYsNq96/uMstro71vzjHy4ismYo/KNdc95nr21qTOLDYsGrddVPNB+tYPiD48l0Lyy6emo1lDixWVr1/cVfzwTqW/+Ic6yMys3K8xmJl1bsG4g/W1lgD8AFDy5pHLFZWFoFhLE91VaMVRmQ+YGhZ8wFJq+jWH3cxffKkQYFhtHaWWGP5gKFVa6QHJB1YzMapLE9m29jmk/dmVhWvgVmjeI3FbBwb72tg1hieCjMzs4o8FWZmZk3lwGJmZplyYDEzs0w5sJiZWaYcWMzMLFMOLGZmlikHFrNh+CKNZiPjwGI2DH8DpNnI+OS9WRn+Bkiz2njEYlZGvd9HYzZeObBY7jVrjcMXaTSrzbCBRdIkSbskvSZpr6T/ktKnStoiqTvdTykqs0rSfkn7JC0pSr9M0uvptXWSlNLbJT2T0ndK6igqszy9R7ek5UXpc1Pe7lT29Gy6xPKmmWscrfANkGZ5M+xFKNOH/xkR8ZGkicD/Ae4A/i3QFxH3S1oJTImIuyXNB54GFgCzgBeBfxoRH0valcr+NfDnwLqI2CTp28A/i4jfkrQM+HpE/HtJU4EuoBMIYA9wWUQclfQs8NOIWC/pj4DXIuKhSj+LL0LZWvxFVGb5kPlFKKPgo/R0YroFcC3weEp/HLguPb4WWB8RJyPiZ8B+YIGkmcCZEbEjCtHsRyVlBur6CfCVFNCWAFsioi8ijgJbgKXptatT3tL3tzHCaxzZ8HZpG21VrbFIOk3Sq0APhQ/6ncC5EXEYIN3PSNlnA+8WFT+U0manx6Xpg8pERD/wITCtQl3TgA9S3tK6Stu+QlKXpK7e3t5qflzLCa9xZMPbpW20VbXdOCI+Bi6RdDawQdLFFbJrqCoqpNdSplJdgxMjHgEegcJU2FB5LL/8RVS183Zpa5YRnWOJiA8kbQeWAu9LmhkRh9M0V0/Kdgg4r6jYHOC9lD5niPTiMocktQFnAX0p/aqSMtuBI8DZktrSqKW4LhtDHr7p02ndNddV+nvGSr1818Ky32lv1kjV7AqbnkYqSPoc8GvA28BGYGCX1nLgufR4I7As7fSaC8wDdqXpsuOSrkhrJDeXlBmo63rgpbQOsxlYLGlK2nW2GNicXtuW8pa+v5kxelOJeV/DyXv76pXHn6+aNZaZwDZJfwPsprDG8mfA/cAiSd3AovSciNgLPAu8CbwA3J6m0gBuA35AYUH/HWBTSn8UmCZpP3AnsDLV1Qfcl953N3BvSgO4G7gzlZmW6miKPP7DmsHobJfO+xpO3ttXrzz+fP7O+wys3vA6T+46yI0LzmfN17+Uef1meZT37eB5b1+9RvPnG+l2YweWOoz1X1yzSnqOnSi7hpOHnXt5b1+9RvPny/wci5VX7TkLT5VZI/mSN0PLe/vqleefz4GlDtX+w+ZxDtTGDl/ypry8t69eef35PBVWp1t/3MX0yZMGnbMY2CLrqTJrJP9+2WjxGksFo32tsFaZ4+05doLvPP0KD95waa7aZZW1yu9Xq/P/D6+x5Eqe50CLeaquNbXK71er8/+PkfM3SDZYni9J4kt+tL48/361Ov//qJ2nwsYxT6WYlef/H5/yVJhVzVMpZuX5/0ftHFjGgUrnHPK6XdEsD/z/ozaeChsHfMkZM6vHSKfCvHg/hnnxsTreTmqWLU+FjWH+at/qeDupWbY8YhnDvPhYmUd0Zo3hEcsY58XH8jyiM2sMj1jGOH+1b3ke0Zk1hgOLjWs+uW6WPW83NjOzinzy3szMmsqBxczMMuXAYpZz/mrrxnL/Zs+BxSznfICzsdy/2fPivVlO+auHG8v9Wz0v3puNET7A2VhZ9a+n0j7LgcWsTo36YPEBzsbKqn89lfZZPiBpVqfiD5asv5bABzgbq57+9bXmyvMai1mNPEdfnbH6tQTj6auLvcZiNkq8BlKdsTpV5KnK8jwVZlYjf7BUNh6mijxVOTQHFrM6+IOlvJfvWlh2qmis8NXDh+bAYlYHf7CU5xHd+OU1FjNrGH/RXP1a8ZyMd4WZmeXY6g2v8+Sug9y44PzMt7NXK/NdYZLOk7RN0luS9kq6I6VPlbRFUne6n1JUZpWk/ZL2SVpSlH6ZpNfTa+skKaW3S3ompe+U1FFUZnl6j25Jy4vS56a83ans6dX+0GZmo6XWEceFqzfRsfJ5nth5kIjC5oeOlc9z4epNDWppdqqZCusHficifgW4Arhd0nxgJbA1IuYBW9Nz0mvLgIuApcD3JZ2W6noIWAHMS7elKf0W4GhEfBF4AFib6poK3ANcDiwA7ikKYGuBB9L7H011mJnlSq3brVt5O/uwi/cRcRg4nB4fl/QWMBu4FrgqZXsc2A7cndLXR8RJ4GeS9gMLJB0AzoyIHQCSfgRcB2xKZb6b6voJ8GAazSwBtkREXyqzBVgqaT1wNXBD0ft/l0LgMjNrunq3W7fy5ocRLd6nKapLgZ3AuSnoDASfGSnbbODdomKHUtrs9Lg0fVCZiOgHPgSmVahrGvBByltaV2mbV0jqktTV29s7kh/XzKxmWYw4WnXzQ9XbjSV9HvgT4Lcj4lhaHhky6xBpUSG9ljKV6hqcGPEI8AgUFu+HymNmlrUsRhytup29qhGLpIkUgsqTEfHTlPy+pJnp9ZlAT0o/BJxXVHwO8F5KnzNE+qAyktqAs4C+CnUdAc5OeUvrMjPLhVYdcdRr2BFLWut4FHgrIr5X9NJGYDlwf7p/rij9KUnfA2ZRWKTfFREfSzou6QoKU2k3A39YUtcO4HrgpYgISZuB/1q0YL8YWJVe25byri95fzOzXGjVEUe9qhmxXAncBFwt6dV0+3UKAWWRpG5gUXpOROwFngXeBF4Abo+Ij1NdtwE/APYD71BYuIdC4JqWFvrvJO0wS4v29wG70+3egYV8ChsF7kxlpqU6rAla8QCXmTWOD0ha3fJwgMvMGmekByR9rTCr2Xi4eq2ZjZyvFWY1a+UDXGbWOA4sVrNWPsBlZo3jqTCri7+PxMxKefHezMwq8nfem5lZUzmwjAE+R2JmeeLAMgbUelluM7NG8OJ9C/M5EhsPeo6d4DtPv8KDN1zqHYctwiOWFuZzJDYeeETeejxiaWE+R2JjmUfkrcsjlhY3Xi/LbWOfR+StyyOWFjdeL8ttY59H5K3LgcXMcstXdmhNPnlvZmYV+eS9mZk1lQOLmZllyoHFzMwy5cBiZmaZcmAxM7NMObCYmVmmHFjMzCxTDixmZpYpBxYzM8uUA4uZmWXKgcXMzDLlwGJmZplyYDEzs0w5sJiZWaYcWMzMLFMOLGZmlikHFjMzy5QDi5mZZWrYwCLpMUk9kt4oSpsqaYuk7nQ/pei1VZL2S9onaUlR+mWSXk+vrZOklN4u6ZmUvlNSR1GZ5ek9uiUtL0qfm/J2p7Kn198VZmaWhWpGLH8MLC1JWwlsjYh5wNb0HEnzgWXARanM9yWdlso8BKwA5qXbQJ23AEcj4ovAA8DaVNdU4B7gcmABcE9RAFsLPJDe/2iqw8zMcmDYwBIR/xvoK0m+Fng8PX4cuK4ofX1EnIyInwH7gQWSZgJnRsSOiAjgRyVlBur6CfCVNJpZAmyJiL6IOApsAZam165OeUvf38zMmqzWNZZzI+IwQLqfkdJnA+8W5TuU0manx6Xpg8pERD/wITCtQl3TgA9S3tK6PkPSCkldkrp6e3tH+GOamdlIZb14ryHSokJ6LWUq1fXZFyIeiYjOiOicPn16uWxmZpaRWgPL+2l6i3Tfk9IPAecV5ZsDvJfS5wyRPqiMpDbgLApTb+XqOgKcnfKW1mVmZk1Wa2DZCAzs0loOPFeUvizt9JpLYZF+V5ouOy7pirRGcnNJmYG6rgdeSuswm4HFkqakRfvFwOb02raUt/T9zcysydqGyyDpaeAq4BxJhyjs1LofeFbSLcBB4BsAEbFX0rPAm0A/cHtEfJyquo3CDrPPAZvSDeBR4MeS9lMYqSxLdfVJug/YnfLdGxEDmwjuBtZLWgO8kuowM7McUGEAMD50dnZGV1dXs5thZtZSJO2JiM5q8/vkvZmZZcqBxczMMuXAYmbWRD3HTvDNh3fQc/xEs5uSGQcWM7MmWre1m90H+lj3Ynezm5KZYXeFmZlZ9i5cvYmT/ac+ef7EzoM8sfMg7W0T2Lfmq01sWf08YjEza4KX71rINZfMYtLEwsfwpIkTuPaSWbx898Imt6x+DixmZk0w48xJTG5v42T/KdrbJnCy/xST29uYMXlSs5tWN0+FmZk1yZGPTnLj5Rdww4LzeWrXQXrHyAK+D0iamVlFPiBpZmZN5cBiZmaZcmAxM7NMObCYmVmmHFjMzCxTDixmZpapcbXdWFIv8HdlXj6Hwtce55XbVx+3rz5uX31avX0XRMT0aisbV4GlEkldI9mnPdrcvvq4ffVx++oz3trnqTAzM8uUA4uZmWXKgeVTjzS7AcNw++rj9tXH7avPuGqf11jMzCxTHrGYmVmmxmxgkfSYpB5JbxSl/XNJOyS9LulPJZ2Z0idKejylvyVpVVGZ7ZL2SXo13WY0oX2nS/phSn9N0lVFZS5L6fslrZOkLNqXcRsz70NJ50nalv699kq6I6VPlbRFUne6n1JUZlXqp32SlhSlZ96HGbev6f0naVrK/5GkB0vqanr/DdO+PPTfIkl7Uj/tkXR1UV156L9K7Rt5/0XEmLwB/wr4VeCNorTdwL9Oj38TuC89vgFYnx7/EnAA6EjPtwOdTW7f7cAP0+MZwB5gQnq+C/gyIGAT8NUctjHzPgRmAr+aHk8G/haYD/w3YGVKXwmsTY/nA68B7cBc4B3gtEb1Ycbty0P/nQH8C+C3gAdL6spD/1VqXx7671JgVnp8MfD3Oeu/Su0bcf9l1tF5vAEdDP5QPMan60rnAW+mx98C/pTCF59NS/8IUxv1S1lD+/4n8BtF+bYCC9Ivz9tF6d8CHs5TGxvdh0Xv9xywCNgHzExpM4F96fEqYFVR/s3pP3PD+7Ce9uWl/4ry/QeKPrjz0n/l2pe3/kvpAv6Bwh8Rueq/0vbV2n9jdiqsjDeAa9Ljb1D4YAT4CfCPwGHgIPDfI6KvqNwP0xDwP2cxTK2hfa8B10pqkzQXuCy9Nhs4VFT+UEprpJG2cUDD+lBSB4W/uHYC50bEYYB0PzBsnw28W1RsoK8a3od1tm9As/uvnLz033Dy1H//DnglIk6Sz/4rbt+AEfXfeAssvwncLmkPheHh/0vpC4CPgVkUpiF+R9I/Sa/dGBFfAv5lut3UhPY9RuEXrgv4A+CvgH4Kf1mUavQ2v5G2ERrYh5I+D/wJ8NsRcaxS1iHSokJ6JjJoH+Sj/8pWMURaM/qvktz0n6SLgLXArQNJQ2RrWv8N0T6oof/GVWCJiLcjYnFEXAY8TWEeGwprLC9ExC8iogf4S6Azlfn7dH8ceIpCEBrV9kVEf0T8x4i4JCKuBc4Guil8kM8pqmIO8F6j2ldjGxvWh5ImUvhP82RE/DQlvy9pZnp9JtCT0g8xeAQ10FcN68OM2peX/isnL/1XVl76T9IcYANwc0QMfPbkpv/KtK+m/htXgWVgN4OkCcBq4I/SSweBq1VwBnAF8Haa1jknlZkIfI3CVNCotk/SL6V2IWkR0B8Rb6ah7HFJV6Th6c0U5lIbZqRtbFQfpp/3UeCtiPhe0UsbgeXp8XI+7Y+NwDJJ7Wmqbh6wq1F9mFX7ctR/Q8pR/5WrJxf9J+ls4HkK62h/OZA5L/1Xrn0191/Wi0R5uVH4a/ow8AsKfxXcAtxBYWH+b4H7+XQR+vPA/wL2Am8C/ymln0Fhd9PfpNf+B2mnzii3r4PCottbwIsUrjQ6UE9n+od+B3hwoExe2tioPqSwAyhSva+m269T2HyxlcJoaStpE0Yq87upn/ZRtPOmEX2YVfty1n8HgD7go/T7MD9n/feZ9uWl/yj8EfaPRXlfBWbkpf/Kta/W/vPJezMzy9S4mgozM7PGc2AxM7NMObCYmVmmHFjMzCxTDixmZpYpBxYzM8uUA4uZmWXKgcXMzDL1/wERgRmeflUJDAAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.plot(style='*')" ] @@ -314,9 +2381,59 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "2021 743449\n", + "2014 1600941\n", + "1991 1659249\n", + "1995 1840410\n", + "2020 2010315\n", + "2022 2060304\n", + "2012 2175217\n", + "2003 2234584\n", + "2019 2254386\n", + "2006 2307352\n", + "2017 2321583\n", + "2001 2529279\n", + "1992 2574578\n", + "1993 2703886\n", + "2018 2705325\n", + "1988 2765617\n", + "2007 2780164\n", + "1987 2855570\n", + "2016 2856393\n", + "2011 2857040\n", + "2023 2873501\n", + "2008 2973918\n", + "1998 3034904\n", + "2002 3125418\n", + "2009 3444020\n", + "1994 3514763\n", + "1996 3539413\n", + "2004 3567744\n", + "1997 3620066\n", + "2015 3654892\n", + "2024 3670417\n", + "2000 3826372\n", + "2005 3835025\n", + "1999 3908112\n", + "2010 4111392\n", + "2013 4182691\n", + "1986 5115251\n", + "1990 5235827\n", + "1989 5466192\n", + "dtype: int64" + ] + }, + "execution_count": 13, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "yearly_incidence.sort_values()" ] @@ -331,21 +2448,35 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 14, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEKCAYAAAACS67iAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFJlJREFUeJzt3XmUZGV5x/Hv4wzCQLOZkUYHZNxiREZAGhWJ2IM74M4RFRA86miiQHA8CXJUooaIC+aoccm4Hxf6KJi4YFwCaXALMoPLiCh4YBQGAVl1gIgjT/5474RyYOjuquq+XW99P+fMmapbt+q+z1vVv7r3vffWjcxEkjT47tN2AyRJ/WGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiqxcC4Xtnjx4ly6dOlcLvLP3HrrrWy33XatLb9tw14/2AfWP5j1r1mz5vrMvP9U881poC9dupTVq1fP5SL/zOTkJOPj460tv23DXj/YB9Y/mPVHxK+mM59DLpJUCQNdkiphoEtSJQx0SaqEgS5JlZgy0CPi4xFxXUT8tGPa/SLiWxFxWfP/zrPbTEnSVKazhv5J4BmbTTsJOCczHw6c09yXJLVoykDPzPOBGzeb/BzgU83tTwHP7XO7JEkzFNO5pmhELAW+mpl7NfdvzsydOh6/KTPvcdglIlYAKwBGR0f3m5iY6EOzu7NhwwZGRkZaW37bBqH+tetvmdXXH10E195+9+nLluw4q8udLwbhMzCbBrX+5cuXr8nMsanmm/UzRTNzFbAKYGxsLNs8S2tQzxLrl0Go/9iTzp7V11+5bCOnr737x37dkeOzutz5YhA+A7Op9vq7Pcrl2oh4AEDz/3X9a5IkqRvdBvqXgWOa28cAX+pPcyRJ3ZrOYYtnAN8HHhERV0XEy4HTgKdGxGXAU5v7kqQWTTmGnpkv3sJDT+5zWyRJPfBMUUmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1Ilegr0iDgxIi6OiJ9GxBkRsU2/GiZJmpmuAz0ilgDHA2OZuRewAHhRvxomSZqZXodcFgKLImIhsC1wde9NkiR1IzKz+ydHnACcCtwOfDMzj7yHeVYAKwBGR0f3m5iY6Hp5vdqwYQMjIyOtLb9tg1D/2vW3zOrrjy6Ca2+/+/RlS3ac1eXOF4PwGZhNg1r/8uXL12Tm2FTzdR3oEbEzcBZwBHAz8AXgzMz8zJaeMzY2lqtXr+5qef0wOTnJ+Ph4a8tv2yDUv/Sks2f19Vcu28jpaxfebfq60w6d1eXOF4PwGZhNg1p/REwr0HsZcnkKcEVm/jYz/wh8EXhCD68nSepBL4H+a+DxEbFtRATwZOCS/jRLkjRTXQd6Zl4AnAlcBKxtXmtVn9olSZqhuw8mzkBmngKc0qe2SJJ64JmiklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVaKn30NXnWb7up6SZodr6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklSJngI9InaKiDMj4ucRcUlEHNCvhkmSZqbXS9C9F/h6Zh4eEfcFtu1DmyRJXeg60CNiB+Ag4FiAzLwDuKM/zZIkzVRkZndPjNgHWAX8DNgbWAOckJm3bjbfCmAFwOjo6H4TExM9NbgXGzZsYGRkpLXlt2269a9df8sctKYdo4vg2tvvPn3Zkh3nvjGNuezvzvrbrLktg5oBy5cvX5OZY1PN10ugjwH/AxyYmRdExHuB32Xmm7b0nLGxsVy9enVXy+uHyclJxsfHW1t+26Zb/9KTzp79xrRk5bKNnL727hum6047tIXWFHPZ3531t1lzWwY1AyJiWoHey07Rq4CrMvOC5v6ZwGN6eD1JUg+6DvTMvAa4MiIe0Ux6MmX4RZLUgl6PcjkO+GxzhMvlwMt6b5IkqRs9BXpm/giYclxHkjT7PFNUkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekSvR6gQupCjVfR1XDwzV0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1Ileg70iFgQET+MiK/2o0GSpO70Yw39BOCSPryOJKkHPQV6ROwGHAp8tD/NkSR1KzKz+ydHnAm8HdgeeH1mHnYP86wAVgCMjo7uNzEx0fXyerVhwwZGRkZaW37bplv/2vW3zEFr2jG6CK69ve1WtKez/mVLdmy3MS0Y1AxYvnz5mswcm2q+ri8SHRGHAddl5pqIGN/SfJm5ClgFMDY2luPjW5x11k1OTtLm8ts23fqPrfiCySuXbeT0tcN7bfTO+tcdOd5uY1pQewb0MuRyIPDsiFgHTAAHR8Rn+tIqSdKMdR3omfmGzNwtM5cCLwLOzcyj+tYySdKMeBy6JFWiL4OJmTkJTPbjtSRJ3XENXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoM78UVB8DSPl/bc+WyjVVfL1Saytr1t7T2N7DutENnfRmuoUtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1Ilug70iNg9Iv47Ii6JiIsj4oR+NkySNDO9XIJuI7AyMy+KiO2BNRHxrcz8WZ/aJkmaga7X0DPzN5l5UXP798AlwJJ+NUySNDORmb2/SMRS4Hxgr8z83WaPrQBWAIyOju43MTHR8/K6tWHDBkZGRmb8vLXrb5mF1sy90UVw7e1tt6Jdw94H86X+ZUt2bGW51914S2v191Lz8uXL12Tm2FTz9RzoETECnAecmplfvLd5x8bGcvXq1T0trxeTk5OMj4/P+HlLW7pKeL+tXLaR09f2Mso2+Ia9D+ZL/etOO7SV5b7/s19qrf5eao6IaQV6T0e5RMRWwFnAZ6cKc0nS7OrlKJcAPgZckpnv6V+TJEnd6GUN/UDgaODgiPhR8++QPrVLkjRDXQ8mZeZ3gOhjWyRJPfBMUUmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUifavRTVN/bgM3MplGzm2ksvJSdLmXEOXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekSvQU6BHxjIj4RUT8MiJO6lejJEkz13WgR8QC4APAM4E9gRdHxJ79apgkaWZ6WUN/LPDLzLw8M+8AJoDn9KdZkqSZiszs7okRhwPPyMxXNPePBh6Xma/dbL4VwIrm7iOAX3Tf3J4tBq5vcfltG/b6wT6w/sGsf4/MvP9UMy3sYQFxD9Pu9u2QmauAVT0sp28iYnVmjrXdjrYMe/1gH1h/3fX3MuRyFbB7x/3dgKt7a44kqVu9BPqFwMMj4sERcV/gRcCX+9MsSdJMdT3kkpkbI+K1wDeABcDHM/PivrVsdsyLoZ8WDXv9YB9Yf8W63ikqSZpfPFNUkiphoEtSJQx0SaqEgT5NEbEkIpa03Y62RMRDIuLEiDi47ba0YdjrB/tgEOo30KcQEUsj4jzg68C7IuKJbbdprkXEXwPfovxmz6sj4m9abtKcGvb6wT4YlPoN9HsQEdt03H0McGFmLqMconl8RCxr5runs2UHXkQcHBEPbm4HcDBwSma+EjgdOCQixjser8qw1w/2waDWb6A3ImKHiPhwRFwKvDsi9mgeeh7w6+b2BPBL4BWbnjbHzZxVEbFnRPwE+EfgExFxcJbjWvcEdgXIzAuA7wEv2/S0Nto6G4a9frAPBr1+A/0uzwC2obxxdwBvjohFlM2sZwFk5h+AM4EnNvfvbKep/RERu0XEDh2TjgDOysyDKF9eL4mIhwOfo+mDxr8De0XE1oPcB8NeP9gHtdU/VIEexcKIeHlEfDsiToiIhzYPPwy4IzM3Av8C3AQcBXwTeEBE3K+Z71Lgyog4YM4L6JOIeGREfA34DvDWiNj0s8f/C2zb3P48cA1wKGVt5C86tlpupPxq5t5z1+r+Gfb6wT6otf6hCvRm0+lJwEuBdwJbAx9pHr4GuK75xr2S8mY9lPIG/4y7fgJ4K+CGZvrAiIjtOu7uA1yVmUuBc4F3N9NvBP4QEdtn5o3AZcADKbV+D3hdM999gT8B62a/5f0x7PWDfTAM9Vcd6BFxQES8IyKObe4H8Ejg65n5lcx8J7BHRDwBWE/5Zn5k8/RLgJFm2r9SdoI8i/JlMAr8eE6L6UJE7BwRn4yIC4HTIuL+TR8sA74bEZGZXwZujohDKVsf2zeP09zfBbiTstWyS0R8BDgD2JiZ1811TTMx7PWDfTBs9Vcb6BHxKOBDwO+BF0bE6yj1LgF+37Fn+pPASygBvRF4QjP9Isqe7dsy83zgJOBY4EDgbZl553zau70FB1FqOoSy4+ZkYAfKj6ntmnf9kM+nKH3wA0p/PRMgM7/fvMbCzLwEeBVwMfDPmfky5r9hrx/sg+GqPzMH/h9lTfoVlM2ohc209wAnNLfHgPcBhwNPAb7R8dzdKZteUAL8h5QrK+0LfAl4QMe80XatW6h/AeWDdh5laGhxM/3zwPHN7QcDpzWP708ZO1zQ0X+/bV5nCWXr5LXAJ4APAtu1XaP12wfWP/W/gV9Dj4i9KTsunwOcAryxeWg95bqnUL5Rvwu8ADgH2DUiHh0RW2UZL18fEU/MzHMpP6/5DuCLwBmZ+ZtNy8rmnZ+HDgOeDbwFOICyfwDKETqbtjiuBL4NPDMzL6SsrSwHyMwNwAXA/pm5HjiaMqx0DfDGzLx1juro1rDXD/bBsNcP9HYJulZExL7ArZl5aTPpscClmXlsRDwGODUixoBJ4OkRsW1m3hYRPwZeSDmW9HPAK4H3RcTtwFrgiub1Pgx8LjNvmbuqptaM9WVE7E/ZNPw2cHaWQyn/Erg8M8+NiCsoZ7Q+DVgDPC8iFmfm9RFxGXBrRDwIeD9wVETsQrna1A2UzU0yczWwes6LvBfDXj/YB8Ne/3QMzBp6RDwsIr5DGfN+c0S8tHnoTmBds7Z9EWVT6QDgNu465Ajgj5RNql0pa+E/bV7rPOD6zLwKylr4PA7zg4CPU/a4PwV4ezPLncClEbEoM6+g9MGjKWOBV1OOrYWyV34B5X0/i9IPRwL7AatyHh1P2ykiFjT1P4myCTxU9QM0tWWUsxOH8TOwwzDXP21tj/ls6R+wHXBAx/3DgPc2tx9H+fbcAziGspm1pHnscMp4+abHzmumb0MZblnc8Zr7Avdtu9Yt1L8t8Gru2prYCvg74DXN4zsDP2lqOIIyNri0o69WUa5wfhhlC2RHyj6Cr3XWDNyn7Vrv5f1/BeWPbiVlR9bQ1N/Rvu2BsylXBAM4cVj6oPkbOKb5uz1r2Orv5t+8XEOPiJOBy4GvRcRoM/nplGPDyXLq7Q+A4yjHkO5OOWYcylj5PpSjUz4F3BQRn6bs7PwF8P9jYZn5w8y8Y/YrmpmI2BX4KjAOfJqyE+f5lC2PjQCZeRNlp+3xlHHCXbjrkMvzKcfb35GZXwU+RjnD9QOUvfl/3LSsnIdrJM3xwudQ/vg+AjyNsv9jf8qaWNX1b2YR5XyJh0bEYsrnfAHU3QcRsRVl39fhwLsy8wXNQ/tumqfm+rvW9jfKFr6ZxymbSx8FVjbTXkcZL9s0z57Alc3tU4G3djx2IbBvc3tryiFL+7dd1wzqXwQ8ruP+sZSdPMcAP+iY/kDg6ub2ayinKu/cPP8rwIM65l08F23vYx/s1HH77yl/tEcOS/0d7T4GeBfwJuDllNPPLxyGPqAcmHDkZtOOAC4Yhvq76rO2G7CFN3LToURHcNeQyU7AzcA2HfNdSPnG3gn4AmUT6z8p38Bbt11HD/XHpn/N/cd09MMNlONnN837rU3hD/wT5YifG4B/aLuOPvTDDpT9HNcCb23u3wCM1l5/x3v/MsrQ2/OBzzbTrh+SPjiMcmLP6ZSDHN5MGUq9Edil9vq7+Tcvh1wy80/Nzf8Cto+IvTPzZsq4+as6Zr0I2L557DjKsMp/ACuy7PkeSNnomHQCZW0FyvjfiQBRfl/mV8CmQyvfQtmSWZKZ75ij5s6azPwdZWjt8ZSd2S+kDJu9Kopq6+94/w+hDDudA+wWEW+k7OxfAXV/BrIMlVxBCeejgUcBz6V8Bl5d+2egG/HnuTH/RMQHKePhr2+O8vhbSrDvTDlh6JCOL4DqRMRulPG/4zLz0ig/JraC8uFeAvwo5+MZa30WEftQvsy/Txkn3YtyqFm19UfECGW4ZWtKzX9FORnmZMqa+8Opvw+2zczbmtt7Uz7736Wcml/9Z2CmBiHQ96EctfJkyof6Nspp+LcDH87Mi1ts3qyL8vsxTwLeQBlPvYqyiXkE8PMsh2pWLyJ2p3yxvTgzb4iIo4CLM/OHLTdt1kS50Mq/UXbgnUE55O7kzHxa83j1fdApyi8dfhQ4IjNvHLb6p2MQAv1FlEP3bgPeRtnjXc9e6SlExHeBh1B+1e1q4C2Z+ZNWGzVHImJHyhf5Syg7wVcBH8jMP97rEyvVnAzzfGAiM69puz1zISK2plyrYNOQy4eAD2b5mWttZl4HekQ8mnIa/pmUHUID9ZO1vWoO3TqFMo74mUHeL9CNiFhIGWb5A6X+oXr/N4mIBcCdOZ//WGdRRLyKcrjqp4f1MzBd8zrQJUnTNy+PcpEkzZyBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkirxfxo0Erf7i0deAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "yearly_incidence.hist(xrot=20)" ] - }, - { - "cell_type": "code", - "execution_count": null, - "metadata": { - "collapsed": true - }, - "outputs": [], - "source": [] } ], "metadata": { @@ -364,7 +2495,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.1" + "version": "3.6.4" } }, "nbformat": 4,