From d94a8171a80353ef4820f1443fe32674c072ea1a Mon Sep 17 00:00:00 2001 From: f7daf5521663c3c5d45bc52e37749a77 Date: Wed, 12 May 2021 07:22:57 +0000 Subject: [PATCH] =?UTF-8?q?exo=20r=C3=A9ponses=202009=20et=202002?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- module3/exo2/exercice.ipynb | 2251 ++++++++++++++++++++++++++++++++++- 1 file changed, 2248 insertions(+), 3 deletions(-) diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb index 0bbbe37..0c5a3de 100644 --- a/module3/exo2/exercice.ipynb +++ b/module3/exo2/exercice.ipynb @@ -1,5 +1,2251 @@ { - "cells": [], + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + " %matplotlib inline\n", + "import matplotlib.pyplot as plt\n", + "import pandas as pd\n", + "import isoweek" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [], + "source": [ + "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-7.csv\"" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020211774750290965917410FRFrance
120211674780289166697410FRFrance
2202115711215762714803171222FRFrance
3202114711197799414400171222FRFrance
420211379714628913139151020FRFrance
5202112711520841514625171222FRFrance
620211179386667812094141018FRFrance
720211079056645211660141018FRFrance
8202109710988793814038171222FRFrance
9202108711281836114201171321FRFrance
102021077135611031516807211626FRFrance
11202106713401981016992201525FRFrance
12202105712210898815432181323FRFrance
13202104712026882615226181323FRFrance
142021037891363751145113917FRFrance
152021027779554301016012816FRFrance
16202101710525775013300161220FRFrance
17202053711978840615550181323FRFrance
18202052712012828515739181224FRFrance
19202051710564757413554161121FRFrance
20202050770634744938211715FRFrance
2120204975026314569078511FRFrance
22202048766834312905410614FRFrance
2320204774999296370358511FRFrance
242020467375219635541639FRFrance
252020457369620165376639FRFrance
2620204474391237564077410FRFrance
2720204374376250562477410FRFrance
282020427400019796021639FRFrance
292020417396120995823639FRFrance
.................................
15571991267176081130423912312042FRFrance
15581991257161691070021638281838FRFrance
15591991247161711007122271281739FRFrance
1560199123711947767116223211329FRFrance
1561199122715452995320951271737FRFrance
1562199121714903897520831261636FRFrance
15631991207190531274225364342345FRFrance
15641991197167391124622232291939FRFrance
15651991187213851388228888382551FRFrance
1566199117713462887718047241632FRFrance
15671991167148571006819646261834FRFrance
1568199115713975978118169251832FRFrance
1569199114712265768416846221430FRFrance
157019911379567604113093171123FRFrance
1571199112710864733114397191325FRFrance
15721991117155741118419964271935FRFrance
15731991107166431137221914292038FRFrance
1574199109713741878018702241533FRFrance
1575199108713289881317765231531FRFrance
1576199107712337807716597221529FRFrance
1577199106710877701314741191226FRFrance
1578199105710442654414340181125FRFrance
15791991047791345631126314820FRFrance
15801991037153871048420290271836FRFrance
15811991027162771104621508292038FRFrance
15821991017155651027120859271836FRFrance
15831990527193751329525455342345FRFrance
15841990517190801380724353342543FRFrance
1585199050711079666015498201228FRFrance
15861990497114302610205FRFrance
\n", + "

1587 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202117 7 4750 2909 6591 7 4 \n", + "1 202116 7 4780 2891 6669 7 4 \n", + "2 202115 7 11215 7627 14803 17 12 \n", + "3 202114 7 11197 7994 14400 17 12 \n", + "4 202113 7 9714 6289 13139 15 10 \n", + "5 202112 7 11520 8415 14625 17 12 \n", + "6 202111 7 9386 6678 12094 14 10 \n", + "7 202110 7 9056 6452 11660 14 10 \n", + "8 202109 7 10988 7938 14038 17 12 \n", + "9 202108 7 11281 8361 14201 17 13 \n", + "10 202107 7 13561 10315 16807 21 16 \n", + "11 202106 7 13401 9810 16992 20 15 \n", + "12 202105 7 12210 8988 15432 18 13 \n", + "13 202104 7 12026 8826 15226 18 13 \n", + "14 202103 7 8913 6375 11451 13 9 \n", + "15 202102 7 7795 5430 10160 12 8 \n", + "16 202101 7 10525 7750 13300 16 12 \n", + "17 202053 7 11978 8406 15550 18 13 \n", + "18 202052 7 12012 8285 15739 18 12 \n", + "19 202051 7 10564 7574 13554 16 11 \n", + "20 202050 7 7063 4744 9382 11 7 \n", + "21 202049 7 5026 3145 6907 8 5 \n", + "22 202048 7 6683 4312 9054 10 6 \n", + "23 202047 7 4999 2963 7035 8 5 \n", + "24 202046 7 3752 1963 5541 6 3 \n", + "25 202045 7 3696 2016 5376 6 3 \n", + "26 202044 7 4391 2375 6407 7 4 \n", + "27 202043 7 4376 2505 6247 7 4 \n", + "28 202042 7 4000 1979 6021 6 3 \n", + "29 202041 7 3961 2099 5823 6 3 \n", + "... ... ... ... ... ... ... ... \n", + "1557 199126 7 17608 11304 23912 31 20 \n", + "1558 199125 7 16169 10700 21638 28 18 \n", + "1559 199124 7 16171 10071 22271 28 17 \n", + "1560 199123 7 11947 7671 16223 21 13 \n", + "1561 199122 7 15452 9953 20951 27 17 \n", + "1562 199121 7 14903 8975 20831 26 16 \n", + "1563 199120 7 19053 12742 25364 34 23 \n", + "1564 199119 7 16739 11246 22232 29 19 \n", + "1565 199118 7 21385 13882 28888 38 25 \n", + "1566 199117 7 13462 8877 18047 24 16 \n", + "1567 199116 7 14857 10068 19646 26 18 \n", + "1568 199115 7 13975 9781 18169 25 18 \n", + "1569 199114 7 12265 7684 16846 22 14 \n", + "1570 199113 7 9567 6041 13093 17 11 \n", + "1571 199112 7 10864 7331 14397 19 13 \n", + "1572 199111 7 15574 11184 19964 27 19 \n", + "1573 199110 7 16643 11372 21914 29 20 \n", + "1574 199109 7 13741 8780 18702 24 15 \n", + "1575 199108 7 13289 8813 17765 23 15 \n", + "1576 199107 7 12337 8077 16597 22 15 \n", + "1577 199106 7 10877 7013 14741 19 12 \n", + "1578 199105 7 10442 6544 14340 18 11 \n", + "1579 199104 7 7913 4563 11263 14 8 \n", + "1580 199103 7 15387 10484 20290 27 18 \n", + "1581 199102 7 16277 11046 21508 29 20 \n", + "1582 199101 7 15565 10271 20859 27 18 \n", + "1583 199052 7 19375 13295 25455 34 23 \n", + "1584 199051 7 19080 13807 24353 34 25 \n", + "1585 199050 7 11079 6660 15498 20 12 \n", + "1586 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 10 FR France \n", + "1 10 FR France \n", + "2 22 FR France \n", + "3 22 FR France \n", + "4 20 FR France \n", + "5 22 FR France \n", + "6 18 FR France \n", + "7 18 FR France \n", + "8 22 FR France \n", + "9 21 FR France \n", + "10 26 FR France \n", + "11 25 FR France \n", + "12 23 FR France \n", + "13 23 FR France \n", + "14 17 FR France \n", + "15 16 FR France \n", + "16 20 FR France \n", + "17 23 FR France \n", + "18 24 FR France \n", + "19 21 FR France \n", + "20 15 FR France \n", + "21 11 FR France \n", + "22 14 FR France \n", + "23 11 FR France \n", + "24 9 FR France \n", + "25 9 FR France \n", + "26 10 FR France \n", + "27 10 FR France \n", + "28 9 FR France \n", + "29 9 FR France \n", + "... ... ... ... \n", + "1557 42 FR France \n", + "1558 38 FR France \n", + "1559 39 FR France \n", + "1560 29 FR France \n", + "1561 37 FR France \n", + "1562 36 FR France \n", + "1563 45 FR France \n", + "1564 39 FR France \n", + "1565 51 FR France \n", + "1566 32 FR France \n", + "1567 34 FR France \n", + "1568 32 FR France \n", + "1569 30 FR France \n", + "1570 23 FR France \n", + "1571 25 FR France \n", + "1572 35 FR France \n", + "1573 38 FR France \n", + "1574 33 FR France \n", + "1575 31 FR France \n", + "1576 29 FR France \n", + "1577 26 FR France \n", + "1578 25 FR France \n", + "1579 20 FR France \n", + "1580 36 FR France \n", + "1581 38 FR France \n", + "1582 36 FR France \n", + "1583 45 FR France \n", + "1584 43 FR France \n", + "1585 28 FR France \n", + "1586 5 FR France \n", + "\n", + "[1587 rows x 10 columns]" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "raw_data = pd.read_csv(data_url, skiprows=1)\n", + "raw_data" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020211774750290965917410FRFrance
120211674780289166697410FRFrance
2202115711215762714803171222FRFrance
3202114711197799414400171222FRFrance
420211379714628913139151020FRFrance
5202112711520841514625171222FRFrance
620211179386667812094141018FRFrance
720211079056645211660141018FRFrance
8202109710988793814038171222FRFrance
9202108711281836114201171321FRFrance
102021077135611031516807211626FRFrance
11202106713401981016992201525FRFrance
12202105712210898815432181323FRFrance
13202104712026882615226181323FRFrance
142021037891363751145113917FRFrance
152021027779554301016012816FRFrance
16202101710525775013300161220FRFrance
17202053711978840615550181323FRFrance
18202052712012828515739181224FRFrance
19202051710564757413554161121FRFrance
20202050770634744938211715FRFrance
2120204975026314569078511FRFrance
22202048766834312905410614FRFrance
2320204774999296370358511FRFrance
242020467375219635541639FRFrance
252020457369620165376639FRFrance
2620204474391237564077410FRFrance
2720204374376250562477410FRFrance
282020427400019796021639FRFrance
292020417396120995823639FRFrance
.................................
15571991267176081130423912312042FRFrance
15581991257161691070021638281838FRFrance
15591991247161711007122271281739FRFrance
1560199123711947767116223211329FRFrance
1561199122715452995320951271737FRFrance
1562199121714903897520831261636FRFrance
15631991207190531274225364342345FRFrance
15641991197167391124622232291939FRFrance
15651991187213851388228888382551FRFrance
1566199117713462887718047241632FRFrance
15671991167148571006819646261834FRFrance
1568199115713975978118169251832FRFrance
1569199114712265768416846221430FRFrance
157019911379567604113093171123FRFrance
1571199112710864733114397191325FRFrance
15721991117155741118419964271935FRFrance
15731991107166431137221914292038FRFrance
1574199109713741878018702241533FRFrance
1575199108713289881317765231531FRFrance
1576199107712337807716597221529FRFrance
1577199106710877701314741191226FRFrance
1578199105710442654414340181125FRFrance
15791991047791345631126314820FRFrance
15801991037153871048420290271836FRFrance
15811991027162771104621508292038FRFrance
15821991017155651027120859271836FRFrance
15831990527193751329525455342345FRFrance
15841990517190801380724353342543FRFrance
1585199050711079666015498201228FRFrance
15861990497114302610205FRFrance
\n", + "

1587 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202117 7 4750 2909 6591 7 4 \n", + "1 202116 7 4780 2891 6669 7 4 \n", + "2 202115 7 11215 7627 14803 17 12 \n", + "3 202114 7 11197 7994 14400 17 12 \n", + "4 202113 7 9714 6289 13139 15 10 \n", + "5 202112 7 11520 8415 14625 17 12 \n", + "6 202111 7 9386 6678 12094 14 10 \n", + "7 202110 7 9056 6452 11660 14 10 \n", + "8 202109 7 10988 7938 14038 17 12 \n", + "9 202108 7 11281 8361 14201 17 13 \n", + "10 202107 7 13561 10315 16807 21 16 \n", + "11 202106 7 13401 9810 16992 20 15 \n", + "12 202105 7 12210 8988 15432 18 13 \n", + "13 202104 7 12026 8826 15226 18 13 \n", + "14 202103 7 8913 6375 11451 13 9 \n", + "15 202102 7 7795 5430 10160 12 8 \n", + "16 202101 7 10525 7750 13300 16 12 \n", + "17 202053 7 11978 8406 15550 18 13 \n", + "18 202052 7 12012 8285 15739 18 12 \n", + "19 202051 7 10564 7574 13554 16 11 \n", + "20 202050 7 7063 4744 9382 11 7 \n", + "21 202049 7 5026 3145 6907 8 5 \n", + "22 202048 7 6683 4312 9054 10 6 \n", + "23 202047 7 4999 2963 7035 8 5 \n", + "24 202046 7 3752 1963 5541 6 3 \n", + "25 202045 7 3696 2016 5376 6 3 \n", + "26 202044 7 4391 2375 6407 7 4 \n", + "27 202043 7 4376 2505 6247 7 4 \n", + "28 202042 7 4000 1979 6021 6 3 \n", + "29 202041 7 3961 2099 5823 6 3 \n", + "... ... ... ... ... ... ... ... \n", + "1557 199126 7 17608 11304 23912 31 20 \n", + "1558 199125 7 16169 10700 21638 28 18 \n", + "1559 199124 7 16171 10071 22271 28 17 \n", + "1560 199123 7 11947 7671 16223 21 13 \n", + "1561 199122 7 15452 9953 20951 27 17 \n", + "1562 199121 7 14903 8975 20831 26 16 \n", + "1563 199120 7 19053 12742 25364 34 23 \n", + "1564 199119 7 16739 11246 22232 29 19 \n", + "1565 199118 7 21385 13882 28888 38 25 \n", + "1566 199117 7 13462 8877 18047 24 16 \n", + "1567 199116 7 14857 10068 19646 26 18 \n", + "1568 199115 7 13975 9781 18169 25 18 \n", + "1569 199114 7 12265 7684 16846 22 14 \n", + "1570 199113 7 9567 6041 13093 17 11 \n", + "1571 199112 7 10864 7331 14397 19 13 \n", + "1572 199111 7 15574 11184 19964 27 19 \n", + "1573 199110 7 16643 11372 21914 29 20 \n", + "1574 199109 7 13741 8780 18702 24 15 \n", + "1575 199108 7 13289 8813 17765 23 15 \n", + "1576 199107 7 12337 8077 16597 22 15 \n", + "1577 199106 7 10877 7013 14741 19 12 \n", + "1578 199105 7 10442 6544 14340 18 11 \n", + "1579 199104 7 7913 4563 11263 14 8 \n", + "1580 199103 7 15387 10484 20290 27 18 \n", + "1581 199102 7 16277 11046 21508 29 20 \n", + "1582 199101 7 15565 10271 20859 27 18 \n", + "1583 199052 7 19375 13295 25455 34 23 \n", + "1584 199051 7 19080 13807 24353 34 25 \n", + "1585 199050 7 11079 6660 15498 20 12 \n", + "1586 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 10 FR France \n", + "1 10 FR France \n", + "2 22 FR France \n", + "3 22 FR France \n", + "4 20 FR France \n", + "5 22 FR France \n", + "6 18 FR France \n", + "7 18 FR France \n", + "8 22 FR France \n", + "9 21 FR France \n", + "10 26 FR France \n", + "11 25 FR France \n", + "12 23 FR France \n", + "13 23 FR France \n", + "14 17 FR France \n", + "15 16 FR France \n", + "16 20 FR France \n", + "17 23 FR France \n", + "18 24 FR France \n", + "19 21 FR France \n", + "20 15 FR France \n", + "21 11 FR France \n", + "22 14 FR France \n", + "23 11 FR France \n", + "24 9 FR France \n", + "25 9 FR France \n", + "26 10 FR France \n", + "27 10 FR France \n", + "28 9 FR France \n", + "29 9 FR France \n", + "... ... ... ... \n", + "1557 42 FR France \n", + "1558 38 FR France \n", + "1559 39 FR France \n", + "1560 29 FR France \n", + "1561 37 FR France \n", + "1562 36 FR France \n", + "1563 45 FR France \n", + "1564 39 FR France \n", + "1565 51 FR France \n", + "1566 32 FR France \n", + "1567 34 FR France \n", + "1568 32 FR France \n", + "1569 30 FR France \n", + "1570 23 FR France \n", + "1571 25 FR France \n", + "1572 35 FR France \n", + "1573 38 FR France \n", + "1574 33 FR France \n", + "1575 31 FR France \n", + "1576 29 FR France \n", + "1577 26 FR France \n", + "1578 25 FR France \n", + "1579 20 FR France \n", + "1580 36 FR France \n", + "1581 38 FR France \n", + "1582 36 FR France \n", + "1583 45 FR France \n", + "1584 43 FR France \n", + "1585 28 FR France \n", + "1586 5 FR France \n", + "\n", + "[1587 rows x 10 columns]" + ] + }, + "execution_count": 6, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "raw_data[raw_data.isnull().any(axis=1)]\n", + "data = raw_data.dropna().copy()\n", + "data" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Pas de données manquantes" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [], + "source": [ + "def convert_week(year_and_week_int):\n", + " year_and_week_str = str(year_and_week_int)\n", + " year = int(year_and_week_str[:4])\n", + " week = int(year_and_week_str[4:])\n", + " w = isoweek.Week(year, week)\n", + " return pd.Period(w.day(0), 'W')\n", + "\n", + "data['period'] = [convert_week(yw) for yw in data['week']]" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [], + "source": [ + "sorted_data = data.set_index('period').sort_index()" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [], + "source": [ + " periods = sorted_data.index\n", + "for p1, p2 in zip(periods[:-1], periods[1:]):\n", + " delta = p2.to_timestamp() - p1.end_time\n", + " if delta > pd.Timedelta('1s'):\n", + " print(p1, p2)" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm4HEW5/vvNnC072QNZCJBASEACCWFTZJUAKoig4V4BFW8U8XddUbyKXhdQriJeVFAQZVEERAUuAgJhh7AkEJYkhJyQQEJC9j05J+ecqd8f3TVdXf3V0n3mbJl+n2eemamurqrurq6vvp2EEMiRI0eOHDlUFLp6ADly5MiRo/shJw45cuTIkSOBnDjkyJEjR44EcuKQI0eOHDkSyIlDjhw5cuRIICcOOXLkyJEjgZw45MiRI0eOBHLikCNHjhw5EsiJQ44cOXLkSKCmqweQFUOGDBFjx47t6mHkyJEjR4/C3Llz1wkhhrrq9VjiMHbsWMyZM6erh5EjR44cPQpE9LZPvVyslCNHjhw5EsiJQ44cOXLkSCAnDjly5MiRI4GcOOTIkSNHjgRy4pAjR44cORLIiUOOHDly5EggJw45cuTIkSOBnDjk2C1RKgncOWc5WtpKXT2UHDl6JHLikGO3xL2vrMQ373oV1z2+pKuHkiNHj0ROHHLslti0YxcAYP225i4eSY4cPRM5ccixW0N09QBy5OihyIlDjt0SRNTVQ8iRo0fDSRyIqIGIXiCiV4hoPhH9ICz/byJ6l4jmhZ/TlHO+TUSNRLSIiE5RyqcQ0WvhsWsofIOJqJ6I7gjLnyeisZW/1Bw5cuTI4QsfzqEZwAlCiEMATAYwnYiODI9dLYSYHH7uBwAimghgBoBJAKYDuJaIimH96wDMBDA+/EwPyy8EsFEIMQ7A1QCubP+l5ciRI0eOrHASBxFgW/i3NvzYRLlnALhdCNEshFgKoBHANCLaE0B/IcRsIYQAcAuAM5Vzbg5/3wXgRMrlAjly5MjRZfDSORBRkYjmAVgD4GEhxPPhoS8R0atE9AciGhiWjQSwXDl9RVg2Mvytl8fOEUK0AtgMYDAzjplENIeI5qxdu9brAnPkyJEjR3p4EQchRJsQYjKAUQi4gIMQiIj2QyBqWgXgqrA6t+MXlnLbOfo4rhdCTBVCTB061JnIKEeO1BBC4L5XV6I1d57LUeVIZa0khNgE4HEA04UQq0OiUQJwA4BpYbUVAEYrp40CsDIsH8WUx84hohoAAwBsSHUlOXJUAA+8/h6+dNvLufNcjqqHj7XSUCLaI/zdC8BJAN4IdQgSHwPwevj7XgAzQgukfRAonl8QQqwCsJWIjgz1CecDuEc554Lw99kAHg31EjlydCrWbw+c597b0tTFI8mRo2vhk0N6TwA3hxZHBQB3CiHuI6JbiWgyAvHPMgCfBwAhxHwiuhPAAgCtAC4WQrSFbV0E4CYAvQA8EH4A4EYAtxJRIwKOYUYFri1HDmTdYuQ7kxzVDidxEEK8CuBQpvw8yzmXA7icKZ8D4CCmvAnAOa6x5MjR0chN5HLkCJB7SOfYLaEaQm9rbsX0Xz6J19/d7H1eLtTMUe3IiUOO3R4vLtuAN97bip/9a5GzLpV5h5w65Khu5MQhRw4GOeeQo9qRE4ccORTkfvk5cgTIiUOOHF2E595ajx/dt6Crh5EjB4ucOOTY/ZFBRNQZYqUZ1z+HG59e2vEd5ciRATlxyJFDQS5VypEjQE4ccuRgIHJrpRxVjpw45NitoS7yPsrmXCGdI0eAnDjk2C2Rr/E5crQPOXHIkYNB7ueQo9qRE4cegrVbm7Fg5ZauHkaPRBr9AeU8R44cAHLi0GNw4lWP47RrnurqYfQctFN5kDMOOaodOXHoIdjS1NrVQ6gO5IxDjhwAcuKQo4qQr/s5cvgjJw45qgZpREW5Qnr3wrzlm/D4ojVdPYwehZw45MihIArY3XnUIc+I2/E478bn8ek/vtjVw+hRyIlDjt0a6rrrI1aicrafDhkOi5w2dDy25jq71HASByJqIKIXiOgVIppPRD8IywcR0cNEtDj8Hqic820iaiSiRUR0ilI+hYheC49dQ+GbSET1RHRHWP48EY2t/KXmqCZk1S9Ug15i7KX/xBX3L+zqYeTo5vDhHJoBnCCEOATAZADTiehIAJcCmCWEGA9gVvgfRDQRwAwAkwBMB3AtERXDtq4DMBPA+PAzPSy/EMBGIcQ4AFcDuLIC11ZRrNnahL+/tKKrh5FjN0RXMA7XP/lWF/SaoyfBSRxEgG3h39rwIwCcAeDmsPxmAGeGv88AcLsQolkIsRRAI4BpRLQngP5CiNkiELLeop0j27oLwIlE3SvKzWdvehFfu/MVrN/W3NVDydEJ6MwFO9c55OiO8NI5EFGRiOYBWAPgYSHE8wCGCyFWAUD4PSysPhLAcuX0FWHZyPC3Xh47RwjRCmAzgMFZLqijsHpLQBTaSvmLvDuje21JcuToOngRByFEmxBiMoBRCLiAgyzVuddLWMpt58QbJppJRHOIaM7atWtdw86RA0A2hW9n7uZzLiVHd0QqayUhxCYAjyPQFawORUUIv6UR8QoAo5XTRgFYGZaPYspj5xBRDYABADYw/V8vhJgqhJg6dOjQNEPP0Y3R3NqGNVuaunoYALqGc6j0ej1r4Wps2L6rU/rKsfvCx1ppKBHtEf7uBeAkAG8AuBfABWG1CwDcE/6+F8CM0AJpHwSK5xdC0dNWIjoy1Cecr50j2zobwKOim25xHpz/Hv7xcq6YriQu+tNLmHbFLO/667Y144xfP41Vm3d22Jg6dTdfwd62Nbfiwpvn4DN/fIE9Xuqer1WObggfzmFPAI8R0asAXkSgc7gPwE8BnExEiwGcHP6HEGI+gDsBLADwIICLhRBtYVsXAfg9AiX1EgAPhOU3AhhMRI0AvobQ8qk74nv3zMdX73ilYu3NW74JrW2lirXXE/HoG+k8V/86ZwVeWbEZNz27LNV5PjYOPTUq66Ydu/Dais1lsdGStdvZernKLIcvalwVhBCvAjiUKV8P4ETDOZcDuJwpnwMgoa8QQjQBOMdjvLsVXn93M878zTP44nH74ZvTJ3T1cHZLPPD6ezhhQmAr0U2Z0YqIes757WwsXrMN838QuBW1lvgNR7WnPxVCeG0ScuQe0l2KtaFZ7Pw8TwOAjlm8TbJ3F7opHTFi8Zptsf+tbfwF9LTrqjSq/frTICcOXQi5f8nlwAEqeRu4zaGXWKmHbyrlLWw1yI+qfa5V+/WnQU4cuhA5extH6he3A9/zzjUvrWRb9sb0wy1tJXzzrlewYuOOyg2iGyMnDf7IiUMXoiBjvOUzFkD1vriV1APYWmpqacNyjQg83bgOd85Zgf/6x+sVG0N3Rs45+MOpkM7RcZCWMfmEDVDJ+6BaHVXT7RUWw7eZt87Fk2/GnUclp1GsEia2muZCe5FzDl0IyTnkxCGA723wkcYVw5m979A+qcYgRX2d6iFdwa7aLI3phAEApFFToUpEnPm75o+cOHQh5EKU254HqKzsPfhuqCnaK2qIkv10HirZV9rYX3KxrBb9V04b/JEThy4EZViJuqutfiVQyV0d11Kq5a9Tk/1UrrO091DSkkJ10Iacc0iBnDh4oiPeHcnKuxSS6uKxO3MZaS/NVj/rGtDTN9BpOQc5t6pHrNTVI+g5yImDJzpiTlFZ5+DoW6i/d9/ZXVnOIWjL5CncXlx29+uY+uNHYmWvrdiMRs0ZzQddK1YKvgvVshLsvq9PxVEtU6JbIjJltc9YddHcHed2OW1zynXctteVt+zN1dvw4rJEgF8nXNzcrc+9jXVa4qeP/PppnPSLJ9L35fFQV29pwq8fXczOlV/NWlz+3dzaljhuQ7XpHHKxkj9y4tCl8FNIl2KcQwcOp4sQqV5SikQsx9RF4HcpUmJ218B7/++2l/Hzh97EglXJUCvXPBoRh8/c9KJXe6RZylWPWGk3fIE6CDlx6EJk4xx2v8nNWW29uXorbp29LHObLYbYQhyef2s9tja1AACeeDOIENupa4hHX1ubW4OqTF2VoC3fkC6MeUQcUp3WY5HrHPyREwdP6O9OJWT/BU9TVlEtnINycdN/+SQuu2d+5jbbPHUN25pb8cnrn8MX/jQXpZLAnXM6P1eHgMCP71uAhQxXUK5TFv8wBzMs7EIEHtPyNhWrhHPYHTdXHYWcOGREJRbpsqzdMWF3F1Z4ydpt+LUiApHgFPM2gumzjPlyDvLevrJ8M7aE3ENnY8P2Xfj900vxbzc856zLib2y7vonXPZg1ekcdpNXqVOQEwdP6HOqEnOsbMrqaEz1eu3Jk/vc65/Dzx96E5t38ItwJS2xTCGrE32GO+fm1jZsVMbVmffZpys5Hm4Nb4++4JK7Xg3byNxEj8I7G3bgxqeXdvUwegTy2EoZESxklXmjnGIlRULSk9niphbekibYDVf2ynzFSnLn3NImUlv6VArCY/cu7w5HCCoxC6tFIX3+jS9gZ0sbPnrIXhjar76rh9OtkXMOGVGJhYwyKKR7skLNOPSMMaZs963F80apXJnKbXQmEX51xWYA9kW+1EGcQ7mNKlkJdoYblLT+INUI55QgotFE9BgRLSSi+UT05bD8v4noXSKaF35OU875NhE1EtEiIjpFKZ9CRK+Fx66hcKtERPVEdEdY/jwRja38pbYPSYV0+9uUbbjailkr9WS5kgTxf7lLy3q9XF5ubg2NE96uubdfu9Odk3znrmBRY8lABTb9ldI5vLe5Cau3NFWkrY6ELUBhjgA++4VWAF8XQhwI4EgAFxPRxPDY1UKIyeHnfgAIj80AMAnAdADXEpGMfnYdgJkAxoef6WH5hQA2CiHGAbgawJXtv7SORWVj8Ls4B7Xu7gfd5l5F1nfYVyGttq9mT+uKtcO2Pr+7aaexTmXEShVoBMCRP5mFI66YVZnGOhClnHNwwkkchBCrhBAvhb+3AlgIYKTllDMA3C6EaBZCLAXQCGAaEe0JoL8QYrYItoO3ADhTOefm8PddAE6kSm1lKoSEQrqCc8ttyrp7KKRNlE1a4HDXlnU37ys2UOt1vajBJ40pY61UgZW9s0xZG9dsxfE/fzxzbu9KYXexAOxIpJI0huKeQwE8HxZ9iYheJaI/ENHAsGwkgOXKaSvCspHhb708do4QohXAZgCD04ytJ0LOT9dEja1Zu+GcjnQvyWNZL9eUQ1nFW2u34eEFq8v/VeLQXW9zhdwckm10EnG49rElWLpuOx57Y02n9GdC128Euj+8iQMR9QXwNwBfEUJsQSAi2g/AZACrAFwlqzKnm0x75BOyHVPHMJOI5hDRnLVrk4lLOhMV5RwcE7W7eUi/9M5GPKIsqt4wrD+y2Fes5LOO+egqTrjqCXz/3sjRrq0bi5WiOoy1kuHEu+b6O/R1lrXSrlAXVNPFqedy2uCGF3EgoloEhOHPQoi/A4AQYrUQok0IUQJwA4BpYfUVAEYrp48CsDIsH8WUx84hohoAAwAkoqUJIa4XQkwVQkwdOnSo3xV2ECqxSMs2XPLxUjcTK5117bP43C1z0p/oEp+xZdkuOMt9WrO1axWpPsslR/RMUqVv/NWt6JYodpK1krQIm7NsI+ZkCIhYKewWhh0dDB9rJQJwI4CFQohfKOV7KtU+BkBmKL8XwIzQAmkfBIrnF4QQqwBsJaIjwzbPB3CPcs4F4e+zATwqutnT60hrpRbGssbUV6eG/BECX/zzXNw5Z7m7cgrom9QotlLlFNJZiMpX71AX086fftk37xUwZe0kzkGGUL/1ubdx9m9nd3h/jWu24vDLH0mU59ZKbvg4wR0D4DwArxHRvLDsvwCcS0STEbxFywB8HgCEEPOJ6E4ACxBYOl0shJDeRRcBuAlALwAPhB8gID63ElEjAo5hRvsuKz3WbWvGm+9txdHjhnjVr+TU2uUgDl1lyrpjVxvuf+09PDR/NT4xdbT7BAdMI+diK5XPyUocdtN3n7usSlgadZbOYVeKgIiVwB+fWYa1W6PQ6kTB3Mh1Dm44iYMQ4mnwW5P7LedcDuBypnwOgIOY8iYA57jG0pG44A8vYP7KLWi8/FTUMDx20lqpEmKlAC2tduLQVYpS2ZePcjdVu3pzVoV0RrFSprO6Fj7hwrl7VBEnuE5SAXD+JzbMX7kZB47on9kiS783gS9+Thx8UCV+kW68vX4HAGDzTr/ga53LOai/O29Sp+nrxWUbsHTddr/KWrORQpobg/cQ4l2w5/kvMN1VIc3NvEps+jtLrOQSoap4cdkGnH7N0+2KhWQSYVZ6w7M7IicOIfo1BEzURkNQOB2V0Tn4KaRFFykd0lzjOb+djeN//rijPb5B+cJyXIKNQ7ONrztYdaWFn0I623kudBbnkGZNfndj4Pj3+srNmfvjOAfAPzBjNSMnDiHqaoJboQZfm/7LJ3HYjx4GwLyAnTi3uspDuqP0G/rCXfaQZjaVmUfQzqF3xdLhI/fnxlUJfUFn6RzS9CKHVEkRkGwpd4JzI4/KGkJOFnXOvPHeVmP9ypiy+qGrTFkryXk/vmgNtofxgUzXwForpcwrvbsjq9+HC90xKmux4BfS3gbTZeXhM9zIOYcQaSdgexfpuW9vxPINO7zqdpUTnOSaKoFP/9Gc29i2MHHX66W49RtWjwN7PypCHNx1XnpnI15/N7uIB0g3Vjkv2sM5mMRKuSmrGznnEMIWIXX1liY0axZF7Z1aH7/uWe+6XZEm1FektHpLEwb0qk3XtvY/rYe0Vx/tvFFd4Wbj5/mdLKtMyG53G2ddG8zZZT89PXM/PoS9PCaL/4t/f3HIlnJrJTdyziFEWazELPtHXDErYcXUmYtHnHPoHPhe3hFXzMLMW+embJtvvJKB93riq5+VOFQmtlIFGqkghBB4e31g/WZbx+9++V2cd+PzxuP6dUmikDMObuTEIYRvboVy/Yz9lEoCn/hdOs/QmEK6k2Z1ml6efDNdnCtdl2MN2W1px3Rs5aadeKZxfaoxpem3o+AnLuPESu1f2TtyWpVKAu9tDkOTeA71ltlv4ycPvBGcbxncV+6Yh6cWrzO+FyauKucc3MiJQ4iIc/BD1pdp445deGFpupgyXaGQ7khrjn//vb7Tk+KDZN0swzj6p49i3bbmRLlpDd13SJ/0nXQiGtdExLQn7nh//VgjjvzJLCzfsMNKG+548R2MvfSf2LRjFxas3FIu95mLuti3DEOHuc7BjZw4hEhr4pZVMbxjV/o8xV2Rz8Gnn0pxMdGinc7PIS26+3pgIl4/um+h/TzLMd/7Z6q3c1ebt2OoCU8tDjjLlWHCIhNufvZtAMCKjTvRtyFSh/ps8idc9mCMoEiYuLHcWsmNnDiEkC+H9wKScW5tbWpNfY7qVNpZ1ko+/VRqseU8pMs5HirThd8AFHSJh7RHnbTjypIRT8VJv3gCh/zgoXSdalAXaF8JWJ/6iDjYCFytEvr7ltnLkn3nnENm5MQhRDRXfDmHbGhqTc85dIVYyaefSomerMl+LF2klbSbFgqvZCJdCOdYLDfCd76Z+njXsdtXMfftjd51XWNQM9PZ9AM1hYK1ntmnJvXwqg45cQjBOcHZcPfL72bqJ4uYpLtaK1XqBaOyzkEoZeE4LFdcqXvBKXS7xpTVpDyNWMetzS1oavHfYPjWrcTlbm1yi598lO5EcUJu4zbU+cFxA6b5YxIr/eHppax4qhqRE4cQcq74viP3zFvprsQgy0sY93PoPmKlSiutY8SBzErq7oRKPg/TGrhiY7R7/7cbnsdHfvW0d5vNLX4u5q7n/dD89zzacB/3FSv5coXq7U+jRzBxIz+8bwFOu+Yp73Zc2NLUgtVbujaJVFZUPXF4aP57OOQHD2Fniz20Q6WQpfnuyjlUTOcQsQlMH9k76VVb9Os/cw8Vni/MQN5Zv6McMVhi8ZptrtPKSMNl2MD5sry7aSee8DFjznCDY5yDpQH19rPrveH5dJbO4aSrnsARV8zqlL4qjaonDlfcvxCbd7ZgV2gK19E78yxWEnE/hwoOxgKfbiqmcyi3lyxrTywhX+V9e9wEOvpx+KQutY3BaOKpt5HhQj70iydwwR9e8BtI2Ic356BUtJ6j9MmLlXhw76H67v/sX29UZC1YszVpUt1TUPXEQZfzdrQYo72cQ2fwDqWSwK8fbXTXY16e2UvSO59xIbttSmpf+J6bJqRDso9sA9zMhIbnRnHlg29kal+iI506t3uaZXek87U6Z9JsvFy5Q37z2JKYOK8akRMH7X9Hm4pm2W37+jm0lQR+8dAijL30n/jNY/bF/Y33thgXtqcb1+G3Tyxxjot7wVSHrbRgneDa8Tz0y0uzSHW0p/zVj7yZKNM3Ku+s34EXl7ktgOwba98Laf+8V/tqzmCVp8KbO1Q5Bwc3oILjMvIw3nFUPXEwRubqKGRoX81zYNsc3f/aKlwT7vhveOotY72nFq/F9F8+hb+8sJw9vosRRTSVdTIqoWqHzIc5JW3gPde77Puyt0uslOF5Lly1hXUI04dx8W0vZRuUAt/NdCWmvXovLv7zy2yd9nBpbJ/K7+XMTt9oysrcGJ24dLd4U50NJ3EgotFE9BgRLSSi+UT05bB8EBE9TESLw++ByjnfJqJGIlpERKco5VOI6LXw2DUUbpWIqJ6I7gjLnyeisZW/VMP1af99X5KsE6f9CmlzC2oKRtvwZDrPhav8TfYOv/wRAMDd8yIT3oqZsjIKac68NS3aMzzbfd7W3OpVz4RT//cpPLRgdaI8mdIyddMJeHv8V+BZqm08sjC4vs07WsppcNPcq7jjnEUhrXS6cNUW7Njl52TKcRn6vap2RsKHc2gF8HUhxIEAjgRwMRFNBHApgFlCiPEAZoX/ER6bAWASgOkAriUiaTZyHYCZAMaHn+lh+YUANgohxgG4GsCVFbi2TPCWU2d8cbMsdr4K6QL5vVAucKdKz+4VG6LdGXctWXplCUEFPKR99QHcvbKlkfzJ/VE4i45cQBp8ra2si6dfXx0lTj3khw/h5Xc2lf9bfRaUwT7x5proHEv7+qi3N7dZj0twc1cnGNUenM9JHIQQq4QQL4W/twJYCGAkgDMA3BxWuxnAmeHvMwDcLoRoFkIsBdAIYBoR7QmgvxBitghmwS3aObKtuwCcSO1Z3VJA76ajdQ6ml9WeK9lP51BUYvJ31M2Lmw4yxKHCIpr2LL6+7zY3ZFsC+p0O89CsSmpd5JI2QGN7xtKddskEwnNvRdduJyjx/7va/KyzuLmrp6m1zYFqQCqdQyjuORTA8wCGCyFWAQEBATAsrDYSgCrMXhGWjQx/6+Wxc4QQrQA2AxjM9D+TiOYQ0Zy1a9OFifZFV/k5WBXNnmKlGHHoIOoQd8hLHs8iU+Z0DpEpa9e8oK2WRaZfvT1H1jueGf50dMQz61SdQ4U6uf7JuDFEmtui68tM04d7vLqSupVLal5F8CYORNQXwN8AfEUIYRNWm0LV2ELYeIW3EUJcL4SYKoSYOnToUNeQvaB33NEWC6b2bf1mESt1rAFhgEpxDl0deI8b8y6LWEnlNrvTjptDZ3IOPn35CATuzhh9AIjr3YB075suRrKJFqsBXsSBiGoREIY/CyH+HhavDkVFCL+lkHAFgNHK6aMArAzLRzHlsXOIqAbAAADt56k9oM9Vb4W0Y/Ede+k/8V//eC15wGQ9YbXK8RtV0ScRsAd8Y9lUTiFtdpG2WitViHRw12vjHFxjyJq2syMkqa5n9PHDRqGupvOMFrNcYbHgPz7O0o6Dj9lrLlZyIJT93whgoRDiF8qhewFcEP6+AMA9SvmM0AJpHwSK5xdC0dNWIjoybPN87RzZ1tkAHhUdKE/4+HXPlv0AEot8BXu97fl3EmWmBc220PlGZS0qTzPLOiOEcBIiVyybbArpsD3GWqkzQpRzhH7MoN5e51bImjeBSuUbcD7P8P5W4j77tOAVbkNDTYpNj+4RrhKBb02fUP7Ncg5aWVsFxUoLVm7xJlzdBT4k+RgA5wE4gYjmhZ/TAPwUwMlEtBjAyeF/CCHmA7gTwAIADwK4WAghNXgXAfg9AiX1EgAPhOU3AhhMRI0AvobQ8qkjIITA3Lc34mf/WgSA4xw6WKxkmB925zalni1KadLYJ4GX3tmI6x7nHdxmXP8c9vn2/eaBIL4A3Pj00sTxTAsj4w1d1kNY3qc/PrMsQ2d+6GPRK6jXyD0NIsLYS/+JS/76Sqo+1VuXJrR7Gmse4/ntmPa+uVCydlFT9J9U+gKsLvhq7gcfP4dKipVOu+Yp/PC++RVrrzNg16wBEEI8DfP8O9FwzuUALmfK5wA4iClvAnCOayyVwJadcTtoPTCZa4J/cP+h1t2PEALzlm8yHzeeZ+5T3cG4PKQlTIv0Wdc+azz/eR/rGGUAzzSuSxzOpJAOv/kc0tlf0MF96tCrrugMg0AEHLv/0FgubFuv6jVyO/M5y4L7+Ne5K/Czcw7xHq/6zHZmyBjIwalDE0G/7VkGo5hJHbOxqi36i5X0BV59PioHwnpId7C1km1d6I6oOg9p3QJBdWgC7IvvvkP74FNH7g0AOGT0gHJ5qSTKMuo7XlyOj1kWYBObb3uJ1WxedsV1hWTwngu8Hh20vUjrIQ0Eup17X4kUmPr9/e2nppR/27gafXNqE8e4OIcv3z7P3JEnsqST5SAdHm0gULuswuSZTjqUsYs0YqWkaEghDkU1MZD7XNP7tH5bM4644pFUTqQA8Pq7PStPRNURB/1x68/ftsAWiXDcAYGV1Ij+DeXyb9z1CsZ9J5CQLVptjy1k2ozY+lUnuO39Uic8t8DrC0AmvYTjeJY2OQWuLSqrjv/8SxSqQa2/3dNbFkgqPfVul67bjgMvexBvr98eu7OV1IzFOIcUobZt9/x799hFGcUCBZxDBcRKmxy5prNuXtIo+HVxkfpOxMRKHtZKpnf18UVrsXpLM2540hyixoSuMs3OAqdYaXeD/mwSxMJyboGoPFHVifP3l6KQEu5nb1JIm6Gyt7b2Va6ItcCpAJtcqXhGKuRY//Tc2+hVW8ShYwaykVq9xqf8bvJMdAMwu1Ot27+/tAI7W9pwz7yV8XtbSeKgkJ1KcQ4u1BSDXtuS6VV/AAAgAElEQVRzGXJacfGiVLjChxvjIKWYU3axUrQB4HQO2zUpgskoQDazxIMr01ESSS61u6IKOQfz5An+m88tFAhyDcm6CzKtz8Ly3sStJsz9XvlAFN6Zm3/e1hLtmLyeFqAsXly2ETNvnYvzbny+XNbejZbvprOovbG+RKlxrZ1TfHbJOnz8umcT9vcc1LF2lmVLTaEAImof5xDeK9eYs0ZqTTO0hFhJJQ5Fu85h2fr4Ym8KnyE3iK8s34R129Lla+hJjnVVRxz0mcapQE0oUGSLblzkU5iCxnuNH7j28UZ8+o9BIhUfzmFbcytWbo4Sw3A28z4LVDgYyyH79WUhmrrYYPmGHVYltW8gNsBff1LU2kxwmDFLqqju316y5xK/5K+vYu7bG71SRcbFVZ0jfigWJOeQvb+v3/kKfnL/QvzOIWZJw8mpSDOnkmIl1VrJzjno/Zj6VZ9/WiLek+I1VZ9YSfuvTxKf2EUFMr+8rkdvTHiuFf/Pg4vKv9sUhbSp/X+/4TlHz/4T2fYyZjFXXM6EkxBClF8y3VSxQJQ58F6WV48AJKRKSkPPv7UeN89eVq6rVvVNRZoW+nxwWcllRW0xuNem5zpyj1541yEuuu/VVV59fSOlaW8ZKR6qzhGof13WSr76R3WupJ1vPcmxruo4h+QE0I5bzpWLWYHIOHHcMnm+vNFi+ePDObyyYrO9Y/injLR6azvO5YgmZ/Kq9qErg2PWQCnfJb2+2ta/5ifDZAPBNZ08cYQ2PhG2J/DJ658rR6XVUVHioIbl0O70aQeP0GuzUB29fFAsFKy8Ve+6jiF+aWCaAtxcs0VWjXEOTKOutUFC5TLbUvpCpK3flag+4gDz5AEcnEM4JwLiwNdxscAmjuNcy84/Zq3kuVpykpf123d5ndsuzsFz7pvsz4FAt1OpwHu+OofT37dn7L+8548sXJOo+5ziD9KvoXLMd2yo2mXbQlyoorOPTxnJ1jHdR5eZaHfIjpZmI6bX3aFYfdU4nOD0c00iIFWslFaHkHMO3Rj6hPK1bQYi2TiRySY/+4O3ySJjnINne0IEbf78X4uwMSQK6zyTndtCN1RK5xDnHJJipS3hTj29WCn9M+CWR3nPt2jmmUSBIlJi36F9U/dnwpadLeV7r1+FryPYsH4NuOWz0xLlm5ic1SP6N+A/jt03VEgn79u9r6xks6t1JvYZ0seip0tCV6ttbYquW7VWYmMraf99xEppdQg9SedQfcRB+5/gHDzaKBAZdh7u800TzrYDVYPApdnJPb5oDX79WCO+d+983P/aKvzKkVc66sO7iwR8Yw2p11Gb0DnY20vTf9a0lHJHyMmm4/1V7mV/a912/OLhN8N248fSiHe4AIwtzA73F588BAN61bIe0kII/OdfXu5wq6n1DmsfQrrIqvp7qRL3Woe1kn4TzMRB5RzMz/+QHzyUKMutlbox9JfZZhetY8naQC9QUyR2Unz59pex1rE7NzXfv6GWLV+3rTnel+daVChE1knNLW344p9fiu14bbAuiE6dSvrFMqlzsDsrpVnuswbBk3F19PmwRnu+ld4IPvD6qrDdeMN96vxiPXH/AbDPrcwJw26d1ZG4KiSGJthCe3Bj1N9L1UKqxmGtpHOdpnWcHJxDc2sbdu5qw2bGKTDnHLoxXJPeZpooWfv6miJrlnffq6vwMJMb2Kf//r144jD1x49g9pL10fnW1iMUFdv1tAukSay0accuLFmb3vHHJRvWh+fiHNK8XllowzHjBpcXGV1Mccvst1ONxmU+/IUP7hf7L3f9equ2QIA6OI9ibpSRmJSSC6Ny4z92KK/HsOGUScNTn6OP4YzJexlFXoBfZFUVqn7FRyHt4hoB3hLvxKuewIHfe5Ctn+sceggOvCz5AJ+0mAtK4tBQW8js0GPaWR86Zg/jOWpIDt8dXTxlKL9EJhe6AKYF7fRrni4njjfBl3NQPYD1cwoWqx0XKrHjrSkUyqI81wLhetf1wI46Zh67b+y/vHZ9QezbTuLA1wu+Oc5Bva5+DTW47T+O8O5/SN86DO5b71XXdntPmDCMHZvt3KRpOn8ya8rqaKt8rlJ+0Z9fwltr45aGtkCPPSmBUNURB3VOpIlfA0Qyy/qaApozOvSYpsYAA+eQPN9vctV4Jkj5ywvJnBMm1tdl7w74L85fuHWu8ZzY4sbpMCztqonpgWxipVpFbNjm2Pm7rrfJIbPX1QNl4qDV613vr3PgDJBsuiBOdKMS7JpCAf3q/eZn2GIKIm2uWCwE4Wp8HUcBXod41qEj8cdPH46DRsaDZSba04ou/TuTrAvJzcwJVz3BD9Dj3O6M6iMOKXaiI/foxZbX1xQzcw7tjcVzw1NLveoVCuTV1a9mLU6UtYf19T11ztsby7/1Z6ISIZb9t7R719wVWkl66lAsUHmH57oe18suFbqmajpXV6bpYX25IUnjT8ERRG7eR7od+2JeW6RURDZN+G5bvwUio2UgwD8bte5rKzZj1eYmtJYEjp8QpLj/20VHo65YMFgr+Y250gYb3RXVRxxSPJxE3KXwu6G2UPFQAL7j4sRe76xPyj1rCtl1Dr4pMjlwL5jr0mwGHGx77fDD8EFNscBa97D9tbMv0t7AYplzCFq++bPT8PBXj02VypMLL8LdF9U0W0eMcyiSt6gK4MVUJriJg3mTw80Dddwf+fXTABAL6T5l74EYN6wvb+nkOeb2KJVzzqEbI82jSXhPh/+LBco8QYwTvR3LzLE/eyxRNnoQz/X4jKc9nEOWuZ/W6Y679zIAmn4kk1hJeb6uVAK+pqzGOD3a/0LYoaRNfetrMH54P9TXFLFH76Ro56212/DGe/Hgf3qcKMCkkDbXUG9xEJyPHT6L9oYAlyiQJDRpOAePdguGep6Dbg9x6EG0oQqJg+PpHLhn/6iuHsE1/M9Zd/jCuBAqxZWwnU+TPUuHK4CaDVnGbjvDd4f36greTNdrTVNWvpoCoaCIlQoO6tDeHbK+y9d1DqrYaeregxLnc/JubpevK03VerxCOiqoTc05ZH8/VMhcE8a1WCnvE/qB+CzcReI3d74jzhwjCu3bBHY2nCsIEf2BiNYQ0etK2X8T0btaTml57NtE1EhEi4joFKV8ChG9Fh67hsK3gojqieiOsPx5Ihpb2UuMw/ZoBvSqRa/a6JboL4xU8gaOORn7d9OGipi7BQ55QTsPvP5eqvG4fDVc/aaFVUzElNlELMmERh6LWnjOLZ+dhse+cVxsIXSdb3rZRw2Mc24+3raAIlYSciMS780H+pC3NLXgvBtfMNbjdvpqCPm6moKTg9Lb9SWaLjFLDcOl72otYcnabeVzv33qBLzwnZOCcXt0XCjwsdF8xpzwO0qJHmTJ6sU53ARgOlN+tRBicvi5HwCIaCKAGQAmhedcS0RSk3YdgJkAxocf2eaFADYKIcYBuBrAlRmvxQu2CVCjKHFn3jIn4fC0XxgqoaCYd7h2Kvpknb+SD5Cn1rO1ef5Re1v7U9vzmew+O5kTJgzzVohmEyvZGkwWcd7Ccoctq4/coxfmfe/kVOM4dv+hGD2oNwqKEpQT0agwqSZ0p0b1MvopZqm6Qlp21w6pR2KXv6OZN56QfXM7fXXx7FVXk06shPbrYoBgk1RTLCS8ir9/73yceNUTWBuKEhtqi2gI56fP4mvkHJgbrEfCvezu1xN1TOdy2K10DkKIJwF4ZJ4HAJwB4HYhRLMQYimARgDTiGhPAP2FELNFcBdvAXCmcs7N4e+7AJxIXtu9rDA/nIKixH2IcWZrCLkK1YLiD08vtfb23FvxW7dgFZ8cRp0zJj+D+poCenmGUajkHEyTwzfL5Le9WN7taUMc0rcOe/SuSy1WAuJRd506h/BbV+LrlsTqZahXRAT8+XORD0HZCU4kh+Z7Z/W3x+Wrwe301fvep67ox4GV2zNbPx0yakDsv1otMYaSQE2BYjnUAWD2kiDKr/RALlD0nHzESgWDzpA784I/xDmu9dv44JW+eogeRBvapXP4EhG9GoqdBoZlIwEsV+qsCMtGhr/18tg5QohWAJsBDG7HuKywPZwiYxnxndMOxA8+OglAJMcvKPXeYTwkVeiLW2tbCR8YPyThGavWsmWgSsOuz35rvbOej2FS0cCGc8iy47WGCGeOsTb7hvOzbDPIEnVXh7wvP7xvQaw8mTxIsL+JgGPGDcGeA4Kc5JHzoiROxJ5ng8452AIpAvxOXz2ld10xlc4BCDhSbrwXHD1W68c8ttaSQG2xkCC8iTNCq6aAyCXbSwZ29J9XOrY08XmyfUVNPSmHdFbicB2A/QBMBrAKwFVhuSmqiy3ai2ckGICIZhLRHCKas3ZttsQntkdTZGbNoD51qA9l3HIHrXIODbX2W5h4UUXABh8wIh7NU+3WNNEKlG6HctvzSQc3HT7Z4Xx9JkwvZ9Yw5kB68US5qXYwn2oyJ2f/YQU94Y2+0xaG31K0o4ayAKLFOctl6LYIRvNp1cBC1zkoo2yoLaaPZyXs5rM+aCuJmENieWwi/i3XflMofZ3zLRaI5aaeXeLeTJnyeuzyNP/e3XQOCQghVgsh2oQQJQA3AJAxglcAGK1UHQVgZVg+iimPnUNENQAGwCDGEkJcL4SYKoSYOnTo0CxDt9tVF5KLQe+6IlrCJ1pbk1RWNzhk8XqWMyFEaKJnfkmMnEOKHbxvPR9iU1NgXGgZmDgbs2+HCI+b22QVhxaHLn3HkSUqq7rIuO6jPK6HydDFUWo7apN6vbqiJlZSxu+7ruhxunwWpGRu9eh3bbGQzlopnC5ct2mIndQ56JFh5Vgvv39h0GaZwPIiNN1yr2DQObhCwwDxEOCxsXqGxdjtOYdQhyDxMQBSS3MvgBmhBdI+CBTPLwghVgHYSkRHhvqE8wHco5xzQfj7bACPig68gzYFrAxWp77oDXXFcmgL6TFdUIKBucJUJDkHUfb8NI3LtJtPJ1byq+fDORTJjygFVlzJeiYCVN4BWp5J2pmQtFZKdz6AmELa1b+J/dWfe1znoIqV4tZJcj5Fu3q+DRuG9WvAHTOPLP93OV7Knb7E5p0t+NWjkee8NCn1BYFCgwgzEY8GYW6nrVTC8g078MZ7W/HyO4pHfXiOjDIc5xzcYqU0YlIdpmyKvvnZexLn4IzmRUR/AXAcgCFEtALA9wEcR0STETzaZQA+DwBCiPlEdCeABQBaAVwshJAr7UUILJ96AXgg/ADAjQBuJaJGBBzDjEpcmAlWnUMhsNo4/ZqnymW9aos4bv+hEELg9IMDmqjGomlzeNLqL0hJRJ6fpnGZdQ5pHHD86vnISm1eqioKhnqmMZeEAErA6+9usbTKcA4eg2mPSQMRGZPu6JCLjMlfQa8HBOM/d9ponHXYqEQU1poE5xAhTS6ACaq/jodCWsX373kdd8+LvIpTEwcLo6k3YxtZa0mUHfz+OncFDh0TqDb1yyGFOHCXmhArOTY7s77+QXz/nvl4mklvy+XKAPyJQ0/iHJzEQQhxLlN8o6X+5QAuZ8rnADiIKW8CcI5rHJWCkzgIxMJS9w4tNc6YHIUtVieha23VF9+SECBKihNs56jjy5JpzQYvp6GC56SmdJxDSQCNa3nrLQlWcch3nRncrr/cr1NfEnwnIss6rJX696rF4WMjpzZ5vhSBlDkSZXBpQrao88snzLo6+o1a1rgCZQufwd25JEdlvr+D+9SVLal2NPOy/qA/RazEzDVdtFtTJLS0mvsd0b+hfI5+2SbTZt2iyoSeQxqq0UPaZsrK7Dw4+35VfOJarHVZpBD8y6a+JCb5ZRYrmkogIErukBwFg4H77S8uTxbCb4xp2XB9x+3lA6f9Vw0OXN3L56YvSvoOM3apIqkLkafLXW7kBBfVc4X/VgNFqud9XomAy0GKgST0XXBNISVxCDlIXiEd/79oddJzGwCuP28KTpk0oizSVTkZfd7Uh0YhJp2cLvrt31BrtDoKxq8quzWOsJ2cw27l57C7wS1WioNTOFMqziE+aUqhQvpdLea72syXb3+ZbatAbrPEcnuWanuFZpO+iHIMJI/1UfwuCMmXs1QSWLqO37kKAcxVorNymKUpCXe1lth7EDmPxRfVLC4zqv7dV+eg10sQf+UJC4ikzilsQL/XajUXcVAJUmqPZuW/vtCl1zkE18MbDsT/L1wViRRVZfKHJo0AEbEZEvV7XV8TzEGTWEkn1AN61bJZ2qLxR+tAwoO93cQhWTZn2Qajc6yKppa2ThVLVR1xsCEQK8VvPud0FttZOh6WvqOUCuk1W5ti5WozeiC18vg8FcOyHxPSLpi6YxYAfPh9eybqcfbj1z2xxDrG7/yD9ziV+PvL8cx8+3/3AXz/3vmOEacDJ1byfb6SUOlWMvo9jmV6Fck+64qRgyWQJHJAPO4Xh5oYcUgvBpLYpXGunAGF3p/E+8cNAWw6B8u4uHzVcu4N7lNXLtOJjjQnN4mV9Pp9G2qwY5d5oVVNsl26JAmbWOlvFx0djYXp8+zfzsbp1zxtPB8ICMOEyx7ELx9JhtjvKFQdcbCasjIPvp6J46O+E67FWp80pVLSrnxQnzqvMBZpxEqV3GDILGTqtZ4zdXR5TLbxvbbCvCPyScPoC5PJahZdhHodphH2a6iJHU/jyyGQ3EHfcmFgDT4wXAQ5zuF/zn4fpo0dxIYPAdxBAnUM7hv0pRsctGiLtClkt76Lbqgt4KbPHB6MWaSfgxxxkFCbWr0lHtZGin5VnZyNqJdDrcQINl9fv6UqkVJhC3OvEtGsU75xTSB+kznGOwPVRxwsi7CaA0GCS8+oik9efoePBiqhi5Wkn4OqdNZ3biYEoYZ9FV8W3UrKpy5fvm2KUpBTzKkclU9fwl+/6kR5x51B5zBvefwZxpzgDLdRvvBlIqLV46zU1GM6MRs3rB/qioWEg5c6/t51NTh0zB7GOaA+ExfnMPvbJ2DPAYGOQoqBJHQRSYF40qsTh161RdQUC9aoxSYRI9evHBtgX+yl6Fcl6rZNVMSdRVCbV8VT+n3s3ytYD/R1weYERwT88pOTw3Flow4yrW7vOv90se1F9REHG+dQiE/qjxyyF8sGFwpRO3McMnM96Jk0ZT15YpSA3VcCkEqsZFl40zqGyfGp3qHcos+1ahMjVFI5Z1tU06LgwTkUpT+C4RqS12bnHACE4hgRq60vTjbusZBC5yAJQ9RvBE7n4MM5lPU8CK2VmHHarONMPgSA3TBAEgdVJ2frhyM4KherbnI4f5VJe/XHM5eeECu3iZWEiIJ2Zp3ykqsy6Tw6AtVHHCzHihTnHIqG58ApXk14aEE8XHZJCBQKQTwdAGUi4WptzKDexvAAHOw6B782ovoMgWSJZvK+rN9mDv/tcw9PZ3QbHKJFNZQVKxFH0yLuBMePsVbzR0iMRyt36RwAlMUxQX2zmMM0JtUROFWgPKVfILnQmayVdJ2Dyq2Z7ovNoELujjnYpooUs6m6IpNHuhwfEH/n1PoqJ895ug/sXZfI+e7Knij7zO58F9ybziQOncejdBPY2FPdWsnEmtucfHQ8sjCe8L4kohf3zR+fipoCYdoVs9jJ/8dPH46Je/XH5p0tGNavHh+/7ll/ayXLsbTTi7sNkZI66omzFtGj0gKRtYjPpegvoQmmtnxepnOnjYn9lzoh11wBzOK7pJhJfsc5AxUxJ0IDBxTsavkxHTJqD6Utw8AZ6GIgNYc3EHIkljkgsX1XwFnKEODcvUntIexxHRHn4Bf2JPJKV8ZVih+XY9f1OCauT3JbXHiNgb3rsLW5JRyX9VKMKHMO7WGJUyLnHBQUNGsl0+7LFpLYBWnKCoRJVMpmgskG+/eqxfD+Ddh/eD/s0bsuXXRUS720ETa53TfXBsFvZyTltT5meaZopjoS/YbD8wk3XqexiKo5qalLaV1UEsAdLyYDHCZMesP/MvHSjUyod1LEIv8IrbQ4axnTffj+RyYpbaXkHEIsXp20lCPwxEYnDrqTHi9WSu6wtzS1YMFKm5e8/dlLi8KC4qypipVMt0IlXibOQT9VCP7eSm5LT5R16Jg9MGZwb+V9yco5BPdNd+jrSFQfcbApqhB/dKZJpSvwAOCD+ycDAZ528IhEWUmIBPX3Vkh76BwuOm4/HDNusH2HonQvzSNt2d+4+yDXBd1ayWfqywnuFRBO8L+T9YKDMt6+HFVaCx4guraS4He/QNy891t/ey05nsT4EI4vuM+cCEXOv1IpCreuj96mc7BlyHNBjm8r44nc0lYyiJX4/qRYiRtmL0ahet6NL+C0a57CuGF9E8ci/YBx6GioiULpS92BSoOMYiWVc9B0DmZxoWCZGck56EH+Dhjerzy2oB/zddggOYe0G7v2oOqIg41yk0YdTOsKF6T05s9OSySAH9K3PlFWKonEzsM0GZMiBXLmX2gNX2TTTuu7px8Ym9xSlvml214ytumza5Tj9eEGiuUXxV1XrWOrLxeDxaHJn7x3PpxDYneuWCKZuiwWHNeQECuJ8DzzKyc5Un2hio2tvLBlXGXYfpMKfRWDwjAWOqyWaODHeOH794n9P2hk/3IAPZspq7Tu49qsCRfkokI4bWbSHCdsEitxoddtYiV9Zx/F3or/T4N5yzfhm397FYC/mLUSqDriwD2bIX3rAQQ+DTHOwWQ7b9jBtzHOQ2q1NVuasIWJB++bkJ2LcbRd2+lt3tlS3l2qE+mQ0ZE8+qzDoujpzaEowOUxqoNVSFt2XCokwUwbPoOrLYMhmlrKosBTX2RTuy72XneelO3Yzgv2JiJ2zYnwDe3cgfL92udfv4Za9nmbHp9N7KpyN/sN7YM+CifR1NKGKXsPxG1KZjy5OEsxkctEVc4pq7WSg3MAgAkjAo56WL/6WLkQ/Nw3mzQH3wWmT19c+cAb5d86Z9KRqD7iwJRd/rGDMO97J6O+pqjpHPg2TDt9LmCeOum+FVL/F5bGlbQ2NlYFJ1a6/sm3Yv/71teWxV7qwtg/dNoiInzuA9HuTXIOKk6YMCwxPm4sOnysuGYeuy9mHB4ogH2uWW3v6cXJKJkyro6pXxPnYJNxq05Sstn/Om1CrI7kAEwGAudMHRX7HwXWsxCHcB7o8m8VadJhmvDTsw5m+w0Qb1cu1NwcMBKHsB3XCGuL8c1YU0sbJu7ZH0eHlnwqJCfg0qWVLBxGfHzxEeqcxiXTDwAAjB7UO1ZeMoiVIuIVb0f+i5I4pX9ui9dEeiBXFOhKovqIA/Ns6moKQb5hTVxkVEiD3xnpDm+6mMpkx21aLhK6CSLo5tT6InHJKQeUxV7qMTWgm7qwc5E+//sjk/DCd06M9ZsYG7PoqpzD9uZWvPRO0gdkxuGjFZ1DNL5LTjkgUTcYb/DduGYrPnPTi0yf0XWpkAu8iXP45SNvRnUNohtV53DqQXGT2vL9ZFtP7vDk8Gz5P4hkHgS1MF5HEqWsxOFHZx6EGZp1FmC+Dili47hHe/gJ4O119hS69TWFWBtNLaWE3kT26sM5qAYb6mKfVucABMH5Dtyzf6zOph27MH/lFvY9NjlNlsVKhrH4QD3n7nkr8feXVpgrVxBVSBwYmWX5BYg/CBPnYLI11zkHU34Dbn1Q60lx0Pu0ZOxqv+9tbsJ7m5vw68caY3V6hfl+S0LEdrXqIllbLODxbxyHzx6zDxvMjShIGFP+z10DJY+pytJJ3/8Xzrr2WebMaEFftj5aPPoYQkLIl2vzTj5ks+z/u3fH4y3JZ2ci8Or90Be+MsGBIhZgEsao49NxosZ9RWlli+F3chLIzUmcc4j3KzmPlqw7SGa8qhhIt6KS18nRWNNCLa/jR/9cwFcIUVssxNrY1VZiw9XE+7SJixSdg9UJLknYudup1/p/fwkCYnJ5Hkycgzw9mlPtlwd+7c5X2t2GD6rPz4Epky+Abu9t2tGbbM3lvJAxcP703NvsZNZfeF1GO6RvHd4/bkhiYVM9pI/8ySzD6KRJpO71GW9r7JA+GNCrFq0lgTZNSc4lZNdRYHbOwaIgjKkU1XoAcLcWVI+Da4MsX7p1Fmc79jw13o02K3yUhxEnljz26aPHlpWkUR8B5OL3209NSZwnNyc2hbTsV+q3Dh87EC8us3vpuyCf5PINO3D/a3GnzYg4cHJ2/v7IzclKzV9CR22xkBBrJjiHFCoj1UPaGiGAUeqz76nCCW/asQtPMWLN6Pz4t4T8K7llm9I9DVrbSok5VmlUIeeQLJOsvs45GFltsivwjtt/KMYO6WPUJbjyIge5c3klsM+GUZqUqpwMt8DLPvRdFheyQYc8p19Zl6HIzC1jFOBFQaa76ZLR1taYjAasp8VEdkbOoaSKC+Lj0DO4qeB2v1G+7OBbz/McjJkSCmn9MorhgtBSKmHt1ma8scqeLMkH8rlxOjN5n1idg7E9slp6SdQUk2a5OnGwZdRLjNUgVjJB7ZvjNFTO35VoyZTfRf7vVx88762MQYoL3L2vZNBKE6qPODBTuqjISGLEwcI2lwTw+rvxiKOnTApCYcTiyyv9DewdRHT85OGjE+2pw2ptE6xsmshvUkjxk0msFOsXoWxdaVfvmpucMhjfRw/ZC0CUhN4Uxz8+Pn5h/dOFR+CJS46LlblMNo/aN1JeqjqOT0wdzVUvw2bFFNM5CDmOeB0iMpvuWpS3kfWKWWdjc8SsDQfX2iZwxBWPlP0SnrzkeOP1xMbBlNk4JZWrTrRleDZF033RwJlc12m74etCDksmMrI1q4qVdCs+vV7QWFTGEUbV+ombLt/78MTyb1O4Ffm3b7iJykIcOHQCbXATByL6AxGtIaLXlbJBRPQwES0Ovwcqx75NRI1EtIiITlHKpxDRa+Gxayh8QkRUT0R3hOXPE9HYyl6iBpZzMCvdOEgx0J1z4hnO/nfGobj74mMwQPo2aOKn/r1qMLhPHc4/aqzWXnxYraUSa9XC5ZtQce+XjilfR0kIjXNgiINilaOKAKS+4VvTJ+DKj/9T75gAACAASURBVB/MLgyjB/XGU988Ht+aPgEfPWQv3PyZaaGHtHF4YZ9ms773jx+CvQf3iZW5OKX3K5Yt5/3++fJvNa0rB/V+6MSsoOgT5BF9rPsP62tMLsNBD+vA06ZgcVPv4R4ahyFFCa1tcQ5jzOC4VU1aCPALuo2Imq69QIS2kjDO1fGhs1uBEc/Wa8m19hnSB2cdFj1LG+cgY2ItXr0VH/6VOT9CWTmsPHdO3KNKErhN2XEHRI6v0v/IJFYqFgh96oqxyMbtQbcgDgBuAjBdK7sUwCwhxHgAs8L/IKKJAGYAmBSecy0Ryad9HYCZAMaHH9nmhQA2CiHGAbgawJVZL8YHdp2D/oLwTyCYNALbtYirDbVFTFb8CUgLd9nUUkLvei7taHzRb20TBmsgs3dsTYHwvjC2TqHgFhUFZcG3gCj7X5w5ea/y8YuO2w+fPHyMkWSOHtQbhQLhmnMPxVH7DS6z4TYCIRDt3GJ32nCOS6xUryh203iP2sS1kdmhsuPX5sI3p08whguxWfbIe8MSawp6km3+6IxJCUW43MjolnG+4BkdxkoK5nFKGHUOoXhHHpUc9b5DA8J/1xeOxr++ciwrnq1nHow6Bqu1UqjrWLDKHoqDs1ZiM7kpYiWXdZhLrAQEIrNsOgdmPlVAse2CkzgIIZ4EoEdPOwPAzeHvmwGcqZTfLoRoFkIsBdAIYBoR7QmgvxBitgjeklu0c2RbdwE4kbitaoVgewG4nQwHudO3sa6yPfUhrtnaxOek1jiHlrYS6+xCZJ6k6gu0cNXWRFJ5E7EB4tds6tcL4c7N9SKl8fJ1PQ9Vvm8a58NfPTZRpt4PfRjq+OTzGzGgARNG9CvXqaspGK3RJDfzhQ/uV44qq3MO3Fh1hTT3Gkg90dUeGcHUnNI2yPnHXYuNCBlDh1PA8QkR1bvls9Nwx8yjAAADetfigBH9AjGQ1jwXAiRGhC3zQerkXMtHxDlE4MxT5SYQ8CAOBkW4OpLaYiFBhNR2fYNqApV1gjQhq85huBBiFQCE39JubyQAVdayIiwbGf7Wy2PnCCFaAWwGMDjjuJzgI0WGL6PmKWpjm0tCWBN8AHGrpr/OWY5nGtfjnQ1J229dEd5aEqzzlk2spL4PXEIV8041vrvhLCB8d+SFcJVx7fY5ncOpShwqlYDK692wfVeinef/68TYQmCKozR+eL9Eme1lV4mmvJS6YgEPfkUjMpS81kU/no73jw+Iw6WnTsCXjh8nryR2Pdw93byzBas2N1n1EnLc//fKSuP4JU6ZlIztxUHOP+6xtRryFPzmsUZs3tmCUQN74YEvfyB2rFjmHKJrPnb/oRiqeRsTc/+4zYkqvnOFom8TvJNavJ4UpwZtNbW0YdOOYH5Jj3vZnoSLOJRjOmnjU4MhBsQhflwlFkvWbjOMN1nWGbmkK62Q5p6LsJTbzkk2TjSTiOYQ0Zy1a9dmGiB3T+WD162LrMShJJyLoCouenxRMF7O6kHf6bSWBIomayVDl+pO+Osn788cT56jWuVIcEQpBeMQ+FcY7sthY/bA6IG9I4W3Uk1NPvMxRsb8H7fMSbQ3vH9D7L+8nlED3TvmLQa/iaCd4Hv7rlbc9kIQcZXPaYHETJXGCBIRAY5/cwt/c2sJT7y51qqXWM5sLkzw5/ikODD53FTOQYqFAOBn/1oEAPjw+/ZK5LaWm6eyvN7I7SbL2HmqZEC0JQQybZ4+c8zY2H9dH33KL5/EhTfPSdSN6RxYhXVS3KV3rxLEupo457B6SxO+ese82Pg5cKXdmXNYHYqKEH7LpAUrAKhmIqMArAzLRzHlsXOIqAbAACTFWAAAIcT1QoipQoipQ4cmo6D6gLunKquvHjdZBgU7o2giHDSST/wuPZV//ehiZ+5XtafWthJqGWulgkWsdIyimNXj+gDRjkxdlDjOwWbV5ILc4ZnG+PcvHoO6moKic2i/WCnef/Bt2u2qOFhxMEx4V4fju+qhRdZotU0tpUTug+SYIqU/YFdIB6IqOwH5RGjp5iMySkPUhYGoq1zXaQftmRi3Kmort0eBJ79szSx+SoZb4bnU4H37x8srrP493Obpps8cjs99YF+ttQCy67cVZ8y6mJgykiRwa0HMT8LkBKegtkgx4vC9e14vh3AHUsYB68bE4V4AF4S/LwBwj1I+I7RA2geB4vmFUPS0lYiODPUJ52vnyLbOBvCo6ECeiWs62oHGrU9MMsBigdBaKrnjpISL5c8ferM8aWWQv1g1ZVwBR8IHaLOF7P7BRyfF6un4ykn74/yj9sa/HTEmUU9t8VNHJkMr+KqApJjApSvVF8xEO8pv3dPbp10fZa2aA1hvXbZjC0Yocd+rdqIvryVSWIblzD09ccIwjBvWt3y93G0f1q8BJ08czvpJJPr2FEdwClogiLHVvyHqp1DeFNmfh27Kapqzqti1fK4hLAsgEg56XHuBQ6e9PVg2Jyrnp4aD4TYcqmWVSSGtQtc56Ndu8zjX0S0U0kT0FwCzARxARCuI6EIAPwVwMhEtBnBy+B9CiPkA7gSwAMCDAC4WQkiTnosA/B6BknoJgAfC8hsBDCaiRgBfQ2j51FFQb+k+Q/rgxe+chL3CXZi0FpHgbJ+BYMK1eYmVktATy8iKAoEX5mX3BBbDLtmrDlUcxL0QA3rV4odnHFQO3xC0F3yr1zFuGLMT5LtM1gvFEy5fjDQRKoUAq6dRRRwSMl3kum1J/USyXYvsOvzWRURZoFtmRTqHZN3IBDT6z8EUvsXUtwruLF0RLvHN6fF4V1GodXu/bQJ4472t5YXQTBySYiCTbszkVPfI1yI9kBQrqdZipqRUQLDgqzG2AKCPYk2oBpLkrmHkHr3wvzMmAwCufXxJWC85RonaYgG7FCKjX7spqB4bYrwTOAdn+AwhxLmGQydyhUKIywFczpTPAXAQU94E4BzXOCoG5aYWCxSTCeqKYZN4pEYSh/BZ2nQTiTKTTF8AVz64CH8JZdzsDqoQLbwTRvTDG+9F3rHqdfALT7LMN1Kkr0JackBua6U0YiW+TqTojbD34D6xeE32dqPfCWulkC674vz4INqVxxcZcw5upY6he7loHbnvIDYNa7k9b44vEJ9IPYKEDFstIfcr6vPlnuGTbwb6NWm2aWLkuM0O935EznLJvtTNDCdWshli3DNvJX6pWX2pIcRVp1jTRvGMySPx5dsjvYGNaNcVC2hRdCZ6VVMfXaWQrsLYStFN1ech54zGQeocXJEiuYfKy/SDl1PdOXCKYVWspM6NS045wGm1YzNl9ZCOWf9H/Uqrl2SDt3x2mvN8Drb4PTp8Xpd31u/AmMG9Y+0mnODCtjnLrRvOn4qdTLBCExI6h1K8XEVRiuUsBASQClq3g6Bv6B3Ziy12UNBvtJmQxhvjGU5Th5FzYI6xmxj47ZQ5nRz7voVXvHNX0iihd4xzUPQmnlt1K+dQQ+X8KUHdeGUfXZlPP5VC9YXPUG5qIoaQ5oxm4xwAoNllysqUcQnCJceijieNSV8yUB5PgJJlybY4+C7mckeripUuPn4/vHTZyThWSaPqS5SCsfHlHAH02U29tY43F1Rh07GcPHF4OWSID5I6B7M+QYqVbHoJWS6EcOpW5r7tH5DP51nIudtWEjjr0FGoLRIOGjnAcZaNwCcJOidC4cRPfHtUJlzlMTMrnDzewkwuNXwHKdRhXpitzoaWthKufvhN4/HaYiHWp969kXNQfv/ozED40i10DrsbrMRBm6xGE7xwYZJss2nisiIkg35MiPhiwCukozGpPeoEx9fowde6J0lEze0JxO9bXbGIQX3q2PZ8Ep8Y7y0zCJ+dlzQldoVhsPWdBjohLOsTWLFhYCJt00vINoUwLyYSm3YkFercJekm3CbITUibCLjcvTyd7NJYK3GLnrxe1xhlPdU51aZz0N/vr5+8f+wdVP2eFq92byrueHF5Ofc3AAzW5n2tLlbSzjfqHFSpgPzZHXQOuxtii2ohueipE9D08tWUiYO/eMHUZ9BvMAnVecyKlULnotVbmtC4ZlusPFYvhawZgNOZzxuhzDeejzdZrUyUPHjjkjCI55jCJo/nIUNEq885IfdmrLiyQufObKasRSK0iYhzsCmkS0Jgi4c1lQ8KYb8+9YBAxNJa4kO8cDj+AIPZOblFY0BkBecaodTJXXLXq+Uym2m2vpnYL4z5pNaTt2XskKQBhA7VB+PTR4/FxZperE6zVtI3H7qDHAfJWeVipQ5ALPpoQpbuJ1aS2biay5wD35dJdGCqp46Hi94od0Y//L8FWrm7Dw6yniveiy7eMIk7pBVNLI+EpT3pGGiDzdJFx85dbuIgkxvZXq6su7M/fubwRFlZIR3+ty38UiH9yMLV8XEkxhfUW894jcf79psH9bXFhHOZzu0BCudQCowOTClYVZx60Ah88bik8QBg0Bsx91xypPpi+uUTx8f+c7ndbQEn9Z26HhFWlSTYOM3TD94T+w3tE/ORGD+8b8IjXPdz0Nv0ye6nxkPraFQfcVB+u0z93JxDSByMAfoYpaPhhdJ1Dq9p4cCBSKykT6qEzsHzqcrTbF6naj0XJOelvnQ2zsEHtkCDOnwUxXJ3FldIx5EltNclpxyA4w8YlijXQzXoaSNVBDGJIqshc5raoB0Z4fOcKaOM9XxQX1NAc0sbjtx3ULns7i8ekxyfIlYKOAf3RBvSt94Y1kRyQAc79Bbcov/LT07GV7VIAEXO+omV4wZf+vvN5ZLwia0kibpq3TZmUDJKrh4+Q+eafDj4iBN1Vm03qo84KDf1/Voic1UBBZgtFMo6hzY758CHBzAQB8Qnss6SynO58BT6C1BxzkHPXGdpT0Bg9ZbIq5gPOxEv02Wz+w6N2HtTtNB6Js2mD+fAwSxW8n8DTbe8zISETclFhpsHxUJcvGN6jkHAuuC+fOWk8fjZOYekGpOO+jBa6P6KN7RM4hQbX1msFGYi86DytjFIA4aByvO36eQeUzhN06bDtXEK+g2gL/icnq8kknW/8SGdKAXHVeLwgfFJUVptTSFGAPT5ZZq/XNrY3JS1QxDc1OvPm4ITDxweOxLI/iNw1gxAknMwwV+sRKFCOirT4/jLepwzkI+1km18ruvYsCMuvjDtBCmUIashJzhrFn18N1wwNfb/M0ePxcQ9++P3T72F1Vub2CWac1Dj8mGbYDNNTOOkJ2HKBaJHvpViBVOAQ3VcpsW3QNHc5HQv5TExh7hLagjFSuqullsoy34Owl/nYJuLnOkzL1ZKcgSmTUfSlDXZnjxX5xwG99EDA1JMrEQELP3J6cx1BP02MBGXVdQWIrFSW0kkRMc7HMThfz7+PqM3e0egajmH0YN6JxXSBE3nYPZzAFSxEg9uwTDZcQNxz05+8ZB28Np4tBfFV1EYKaTti+p8RcS1z5A++McXjza2JyBixOaD+yd3ULo0QuccCgXCUfsNhgzpzMl7udDOrlSOQFL2Hy8N+88gVjJyDpqMWC5IXDIn3YnL9BwLRGgNFxkToQb4+cftOOtrglzO6nw3mVIDAWFtKwn2GgBgrJJ4yMo5UNxfyFgPyTnAiuUKfmIljnMYP6wv9h/eN1lPESuZCLHMI2G6HxI1xUJZCX7ZPa9j/sp43gnO7yLoGzh32hh84vDRipjS2lVFUH3EIfw2iSID8U7w3yQCjBKuBK35hNGWMFlPBDoHpQ9u50Yy2qW2mCU4h+j39z8yEb/5t8PY8cl6Ls5hu7Kj+drJ+2PSXryMWAiB597a4JSdcv4lpvFx12uCj8z2in8uhJ7K1NS8/vLaYFoW9J2eXNS5NLDFQtzxkpsDQEBc5SJjE+340ri6mgKaW0ox6x2OOMjxtDLWSkfvN7jstX79+REnaCO0gY5KWPU/sg0BYGDv2lhZsp6nWCksUpXDXzx+vwQ3Uluksj6uTQgjIS6H03HGFIvG97e5KxLHt1s4B/k4ymLKXCFdeci5wyqLixSb9C7OodymoS8febtsT598XFRWKWtOvgDmfieM6F9OOGMai3wBfnxmIroJAOBcJee1bTF6cVngdOVyvuL8S0z1AgV38piPZQeHnS1teHXFZqtY6b0tTQD4mE4mmHUOwYHv3h3EzJKiG5ZzKFBsgeYISNAXlQmhjUvkjqiZCiXkpkMVo9o86iNrpWh8t/3HkfjGKQckzrXRJ5/MgUBkyho/l2/PJwOiLFLrNjBiypF79MKKjUHU3ZKFc5Am5s4wNIVIsc7dX9OcLomob3nLu0Vspd0NkuJyz7m2UAgUbcUCAGENvKc1yoJlfZmOX10RiG3UPATmqKzJiaG3qU5i0+4zOC/4lsRhhJYfQWLK2HKKcKsYQ4JNucj0K2EkDmEc/5/9643EMVUUdf5Re+OW2W87xyVxxm+eif3X3+nWDH4fNu4n1naphGKBjBuHVscCHfQVLSRWERhzbOrYQYmyQrg5MW2G9PHI/OSm8akbCNt88eUMpU4uXsa0V2BiK1nCZ6j3mtMX9O9Vi+27WrGrtYQbnlpq2cSAfS+T9aLxqc9Nch4mJ84g0izFxt4Tk/10e0ScQxI1xeDhcTLJWL1wwvUJo4D+74xD2XpprJWAeGwbPhMc7wxkM2W1y8+DYzILlomQqFyMTQEq4fJU7tcQV7ZbrXKEwCML18TKH/7qsTGnJM5sMA30O6pvCg5gMsnpMOsc4gda28z+Afq9NSukIyJiI/6+mhOZXtPlhCXHt3j1NryyfJNxLsc4B5fOoZTUoSXqhd9qPVPypQ3b4/k3OCMFjnPgovxKzvVvLwUiINN6XCQZTTeo8OBXPsDWk6bK8nf5fMVEWEdTSxu2NrWWOc1Ih9XxqD7iEH6znEMon4nkwybOIajXWhJ4/7ghscQxKmzxjNixOSc/r3Ow9WsTA8lDV9wf7Mw5OTMQ3335KLtdnspcukh+fHyIcj3tZ6VTjuucz10XHeU8x+YYGG9beN1nwKJzIPW3+dp9c0gXQ89iF8ckx3fxbS8BALYxjppAXBxm1TlQMtzK1L0HJuqV83Qo4+Na3bB9V4LQ7DXAfA/U52xz+vO5L6VS5Pxp9myPxEpbmtQQH6GZOkMlZawm6TCq+810JKqPOJRvavIByoVU2tAfZ3D7Lwfeay3ZZb4p1yzXwmuy3tGJmLoDtbWpT2LTohU7x4c4hFZDv/3UFPa4r+ktp2Dk0F7SoCvYh/WLxGvFAiU4nTRjSIisSiXjop/kHMw6B3V8Jlz24YnGY3q/bSWzGNU0vmbDJiCVzkEzZe1Tn5R0b9oZcLaqwpa77vVMHg8uK6K8fypR4g0ESJ5guoSgHvmGPQkkE/p9a20LdArc/V8d6r+kg6dsObdW6kBwz0+aj25tasWHJg7HfyvZ1VSoE9NuLZI8ZrOFdhETuVg2ayab+k5Xl2ca29Oevk1EUW7PMsjPHxukY2xuacPeg3tj+kH+Ce7Z8VFShszX8+rGiPOP2jv2/4zJUdRVb+9wz3otbcK46Ot92UxZXXUAfmFk2wvbcOmK9L5M1lwxnYPDWokzzdaxalOTcyx6XyP6N7DhTGS/AJz6Hdmeaw5InYFPwEQAuPRvr8XKW0silidbhdz36Wl+c7FSB8Cmc1AtSA4eOcC4k1YnUlprkZ+cdbCx/ltrtxuPyb5KQmDttrhcVY++yckz2fFpL64eW4aDLWLCkfsNBhDISX3aMo0jKvfjHHy4mTT9E1FZn2Rb3L543H7ROYY6A/vUoa5YKKeHbW5tMyYRSoiVLM6GEmmd4DjINpxhVLTxGA02lPfI9mikX4LrGW9n7P9NybAkzjh0LzacCcDrHEw6PsDt96JbK9k4YQB4ujGZN6OmUGB1nJHpc1zn4PNetBfVRxzK1kqcWElRvFp20erEtIli9Ely8MgBbA5pX0inobe1bGejB8Xlqr46Av1IezkHeb1NrSUvEVV5HCl1DonzvXvyh+QibffvP5XAb9YcEJOGo3+vQFzS3FpiQ39wfZkipfoSf+8wKmEbL79jz1ngY4wAxBda232RxN/FOaghuG1jUa/XluK1bK2kcEqmDHRBfTukWM6WyMnUBxB4PheIN4DRPeplE7lYqQNgtVbyFBepx1Zs2mmupy22vRzu9S7oQxrevx5PXnI8Tpk0QqunEAePxVzCJO5QYV+Mgu+mljbWg9l3HGp7XjqHCiukgegZW8M/kLoImtsqSM0rgOaWEmtTH9SL/zcplL3Fhr6cg2dFfXp8+uixzvas1kpI6hw4cDnBbYs5ADQYCLA6JhenJK/DNQOlCa0k5rbNDgBsVKLpHjZmD3zi8NFl0ZQOaUEW5Z4Pvrs950BEy4joNSKaR0RzwrJBRPQwES0Ovwcq9b9NRI1EtIiITlHKp4TtNBLRNdQRb3uIMnFgeqiJscN+nMMrlgxRqpjq0DF74JpzeZNXX+gLvRDAmMG9E4ujumlPp5B233a73XrIOaQkDqYWuWQwnC+G/qjSiLRMKJaJg7mOesx256RsHQjFSoaFS30ePzpjklE/5auQNvle6PDlCNR6/RpqjDo5X2ulAiFMj2tf6HSdEGAOXChh5xwCuKL46mFyjPXK1lRmBzcgmi+qOO7SUw8sn8NxitJjXnIOZVPWHsI5HC+EmCyEkD7zlwKYJYQYD2BW+B9ENBHADACTAEwHcC0RySd4HYCZAMaHn+kVGBeLyFbJPrlsnIPvTksNafyt6RMwYgDvZOYLfWH28cxOs7PsXef2ifSRcevhi12w+Qio68aI/g2Y/e0TEvXUBegnZx2MJ795vHffJvjszGL31rkIAmu2NOGpxeuwcBWvyI3vuM3tzX17Q/l3JcJn+OpsfN+PuKmtpV8itLQJLF5jz7L2mWP2Yc+1ldnmn6zmiuLrTRzCruRCbrNW0jFtn0HlvtoY0/kWLUxKlrhfWdERYqUzANwc/r4ZwJlK+e1CiGYhxFIAjQCmEdGeAPoLIWaLgL+8RTmn4hAW1i+2qFp2n77EoTalf4ALukOPOYWm30usL0B9GTNCHb6cCBfy2QTTLo8ofo2fPmYsu2iq+QAOHbNHu4kwEIkythhs+YPxqZymuS3pzCdTSJoCBKr31rYIfGhiJEZMo/MywYNhDNrz1WXFxG0VIF5MPd66KPptj5Aacbg2yI2QM1ZY2drLbq1ku15pDqtD6kXkubKJbi9WQrAheoiI5hLRzLBsuBBiFQCE39JkYCSA5cq5K8KykeFvvTwBIppJRHOIaM7ate4sYqYBmxBT9FmepI9sHohHVk1D8c86jL18/O6Jt2L/uZwPQHzs1qid2iGbnLbcnqcMfu/B7rSKEkbrHYp7xZoIjq/1TkfCJsIJiJxbN+JLbCaN7F/+bfWQ9rwVvhuXmgybHavOQTv4588dYayXMPN1cA7D+pkNP2S1Jk/rrN8/9Za1XiRWkgs5f9F2s3LeCW5s+B6dO21MWC8o7wlipWOEEIcBOBXAxUR0rKUud2eEpTxZKMT1QoipQoipQ4ca8tK6YNE5+O64vV8mT5M+NTDef544Hr/4xGS2nrqD+e7pB7LsNhC/Njv7Hz/mo+rxFVOZlK7seRZ7fiFEOTzGjMPHGOu5xnfn53kv51s+O817nDa4FNJCuAlXjKhb6qqGDbaNSlprJQndg708PpU4uMw7yT0G/VEdoyXfiteNV7ZFWwVgtQos7749nf42hqbiPzIEpix7UjOhMVTIa+D0ZjUF3gluQO9a9GuowRmTgw2j5LKfWZI0h6002kUchBArw+81AP4BYBqA1aGoCOG3DIyzAsBo5fRRAFaG5aOY8g6BzZQ1Jiv1tFayIRaTyHLOTxXfh2lMYDSJyz8WTU6bCEh9kVyBzyQ+OXW0sZ4Km643FrbDV1ZhgYyyKSDwsUNHtsspTMp2VZwzZRSOZfJNZIHtagmBGMA1bdR7a1t7Vee2uhr/xfdDE4cb6kUVT5wwDLO+/kG2nirCspl6q21axW3tMELmnrO6m7YZV+iLuW8fpxqcOvWMii6F9CGjmQRYBoV0S1spZmAxbZ9BOHfaGJx+MB9puZLITByIqA8R9ZO/AXwIwOsA7gVwQVjtAgD3hL/vBTCDiOqJaB8EiucXQtHTViI6MrRSOl85p+KwmbL6xiTKxjn4ialsi6oqW+diwXDj89U5jNcSnZjbtu1Uo98+lk8uyB23LVhdMCb1t3+/aSyqXPDhHC7680uONvzmi8qV2TgHtb3fnTcllmdBRbMiex81sBf6G8KFxIiDk3Og2Ddfx9pEDK4w9UGZ36bI97n7OiXKfqVpLJeoC4ieB6fgLlJSrLR8ww786bl3sF4xfa0tFvCTsw5OJbbNivaE7B4O4B/hBdcAuE0I8SARvQjgTiK6EMA7AM4BACHEfCK6E8ACAK0ALhZCyFl5EYCbAPQC8ED46RDYAu+pRbZFxlshrSyQvmEsbApGtY3BFrY5bi3i93L6iiBseon4Dr79C6/0c2hpE8YXDtCJete47tgXaU/rpxjHZ+krA+G1nREPI2Hu2DcUN+Cn70jj2e4KU6+3ZyNevg6avhkWZbFc9F1ERM2FrR7TTXqfYTypOxOZiYMQ4i0AiczmQoj1AE40nHM5gMuZ8jkAeIFehWFL9uPrXOTvLOenkFbr2Wz0Y8TBwjkUMnAOvu+pzX5cvUZfzsEk6gAiVrutVLK2V/BcVHVUUndtVwzHc5ObFOu+1krqfHHFQ/JpTxVR2u6zuqi6fEnkLtrWr2qJ9ompo4z1OLBJs5Qi2/vryzno9YwxscK+ZDRiY9Rdy3zjiMOarc2G2p2D6vOQtiT7KXiKJ9RjtolWE+MczGNSu7ItMupCP6ivhThk8HPw5YZ8PE8B/x28KReGbK9UFiv5ibPScA7tkXnrcEXnVRdCkxJcvX++TpiuHAwStnmqGkRY559yrN7T27/W0u/SdVEYmP85O7HPtGIbE1Kj4MnZ+DpJ6lZ0pmcsF/Xbnn/HWs9m8ME5wb3xnn+a2o5A9REHT52D/AGN6wAAEXZJREFUDTEdgWUS1nqaslJsx+0ntulrcViL+2b57bh9r92Xc/AVfdgihxZDa6UWS5hrvd80pqwV5RwcRFh97/v34mX6vroidXf/PkMuER024tCvoRYnTBgWjsFcTzWwaPDcfddb5vP67dl3xuOGJXVkRc95oC/6ptDyej3TM3EFzJRwEXydc2hxeFx3NKowTWgIjnPwVAiq78+PP2aOsuobuliFbVcTE1OlMFE1IeYf4DkBfXUnvhZdNsiQ3cKhkI5dRwp5/B6GRToLbIsqIR4GhMsPDsQXNKvDlHIvuPwHHFyiFMnZ1Frus0qgfUUztRZrKpeHsg2ck1tcIW0+Vx372VNGGUPL6xsm0zu3g4kay8FmYSiD96lobSuhd10Rz3wrGRWgM1B1nIPcwvE6h+h3b8uOVl2kP3rIXsZ6MeuODBZOOtIsfD5QiYjLrO93502xXisQX9xsCmRfFAjlJDS29rJyDjM/uJ/x2P/O4H1NTHByDsr/PvX83Fq3XQ0w56dz8IVLlLIjXKh7WxYwdW6qKW05jB3cO+zX/B5JD+VLTjnA2pYvfHU2Krdqmy+umE8S25v9iNzw/hYjEoZzaC0JHDCiHwZa9IsdiaojDjZrJXVCDbDsKjOZsnpbOFkU0p4Ln6/3pNrcE4wFhYpTJo1wBg70NQX2hXrPbDtaXx2LDttOTjod+cIVTkI1U+xrUEhL0U5wTra+THDFuloX5gixeRabOB4Oe4URZW0Kbhn47qQDzUYJaeDrpDekT3SNtvfStgao+IJlk6HCllGwppAMMtnSVkp1zyuN6iMOVp1D9HvSXv2ZGgF8F764tZLXKRUxoTXlAdBR6SBeWRTc9vb8OJEs4qzffuowZ51/P2IMTjvYL5udS+ylPhKT3kYN0W27iiyE1yUGkgpekz4ESGd66tPv8SExtBGkNFCJjFXsqhx7utG8KZq4V3+vkDIHjxpgtR6UsBHoIuMh3domKuJMmhXVp3MoB95L3nRV1GRV5Gbxc/BciG31fBdc36BcanOf+8A+XufY28sm3jG3F/32N2VtP4cmcblFnySx54AGrNrc5BR7pQ2FU4n5p8JFHOS6ZOPQZDuuKKUqbCKwb54yAf/xgX29xCbD+tU7TTtVb3ffd2XT9hbr8WPHD8VDC1Y725HP5KQD+exzgJ04FBTucseuVsxesh6tJYFeXaSMBqqROITf7b3lc797kjPxfCzwnudDHtC7/eKsQb2Dl40LGaFCLkAj+jfgyH0He7VtQ42n058vfKOAZuGAKsU1yVZcGfdKQqB3XRH/No2PD6XjmHHtfx4qetfaX3W5MLl0RbUFQjL1ThKbdwaLri06b7FA3pkR1RwYe3lE3fXdnHApSFW4dHEScj82aS+z9ZjN0k81Zf3u3a/j7y+9i7qaAt5viTfV0ag+4mAJvOfKSKXC5qEskcVayQafkNpAkLd43vdONoZB0MdUKVO52gzE0Ab1cdh35hkar9CGTBJYq0K6EIYBKZWsdv8qfHJrpIHNZBiIuE3XXBjUtw7bN+zE9z480Vrvh2cchGsfa8REi3g2DSRXf+MFU425oVX4iupPtjhhAsD8lZsBuJ06pc5mnyHmsBYq93bXF47CYWPKedACsVJourp0XWAau6u1VBHdXVZUH3EIvzlrJV9ZvS9UNrISYpZigfCRQ/ZCLw856B693ay6HFGlVA9xYtj+9tQ0jjZxR5bEgZV65eSGwuohjSgMiEts44v9h/fF2VP8vYpdYiVpKeNaBOWG4tj97TvaKXsPxI2fPtx7fC7IR1xbLHhtPHzft5+e9T7r8dVbgkX/OA+CBADDmYirEuoz0K+jSJFCWjVeeHihW6TVUag+4mDRSHtykN5QF63eBvPFtPhVO1ONqpCTsVK7E3V3P37Y/2/vXGOtqK4A/K17L6C8US9CeRcRi1Qe91algtY2KlRTrO/WCEqjtbGpmjQWax9pbBNoqvEZ7W2FaB9qGtuIrdXSpohWW5EKKiKKioqS+kLkERFk9cfs4c49j5k5586ZB6wvOblz9t1nn71mzZk1e++11+rX7faCc9tJjxziuilG4ecECBuliXQuNibh4gvwtysrR06tF/9nETufdMLODFH43xb3Wo3bv7j7RObPOiJWvYExp4VLz3NwE1wwGVQaeRuqsd95K/lUunYamV0pagH0hnMnc/fFxzbs+yvhy5vEFBB09aVPIhtb0DjEXZAO46ErZuw9jkouHxd/I1fY3LpIpzHK0vskDP9aiBuULnXj4E/fxe1fHU4jlfjqFM+leWxrvKjFYd5eQUqvg5bmzgcIP0ZT1ux3xiHMlbWWNYe4zJ91RBf/9WqcPmUY08YmuwgZhR+zLansaUnf+HYGfiRhXi9xn3aHDuh0Fa3F4yYMPzZ/eH6NzuO4cX3SZk+NI4e0ibPwD+GbUiu2G3HtX3f2JNb9NH5K+7i500tHQD2bO73AurNzPEn2v2mlkGQ/e5K5X3Th0hPGxt4kkzb+02xSN4SkF8+CN/CowHaxCNj+pIxDx5x2Xn9vR3hI8cCjSNQ5mjttVNUc0/Xw6FUnxtJvS7PAruiNbllOc0D0+bv+nElVM7YFObd9BCs2vB9Zr6lJ6NUUf0o4bliR0oedXj2a9o5mgy67tiCdImEjB39oPShk3nBfYk+MxdRaqGdhOIxgWtSw4X/c6KpBd8idMUNdR9H/gB5MHBYe/C74+17x2mYurJLeFeAns5ONXD/CpViNok/PFrZ+tDvRYIRJ4ofaiFojaGluYsCB0TfohWcd1ZCZgrgjw1KD3auluUvSJZ+4ayKNIJ9j3AYSFj7DH1p/LiRV577EhKH9mTttFLd+PXq3cC0cGDOccxTBp/vQ0Oixw403M3uyN+1Q6YfYKN7Z1rkzYPP2OLsE0sePJRa1FnPh50cDcEhCu5rj4u+baI25LyIOST/MQHzjUPrVvVqa+PiTPewueWixkUOKzJo4hMMP7VtxQ4r/JJH2YltWNDVJ4k+qd8xt5/BDu++pBF2Nw6AQ19ymJmHDglNjtXncYYdw/6q3GD4o3hN1Enz4Uecu3DOm1pbUJi065rTxq+WvhvrpA8ybPoZ507u/m75WrjtnMrctW0//A/N9y4q7EF4aVbZnSxO7PlHOuO3xLuXVUrumQb7PdAMYdXCfqvlXDx/i3dROPjKZQGD7I19KKIgadJ1WGj8kGYNzdttwxg3uy5TABqRGE7xdnDm1toB+aXHY4H4sPCvc5z9LZk4cUjW0dpG471vT2Lx9V9nO8HfcOsMzG7fsLTtxfCtto9K7TkvJjXEQkZnAjUAz8GtVXZB2H8a29uWFa2dWjBVvJMu9lxwbuXh38YxP7w0NHdfFMgoRSdUwQNcRUCOmMozi0Daq8pR16W9h9Y9ODg2lkwa5MA4i0gzcCpwEbARWiMgSVX0+7b6YYUiHY2LEcgoGUtsXqDVHhFEcHr3qRF57b0d0xSqUPvxkbRggJ8YBOBpYr6qvAIjIPcBsIHXjYOSLpVceHxqwrAj4GwJ3x8z3bBSPEQf1ju0ZVomZE4fQsfyVBHvUffJiHIYBbwTebwSOyagvRo4Yl9DidpZ895Tx9Gpp4rRJQ7PuipFTpo4cxIYFp/Lutp1dYitlSV6MQ8UNy2WVRC4BLgEYOTJe6GPDyJr+B/TgmlPDo5gaBhA7hHka5GWfw0ZgROD9cOCt0kqq2qGq7ara3tq6b81HG4Zh5Im8GIcVwDgRGSMiPYHzgCUZ98kwDGO/JRfTSqq6W0S+DTyM58q6SFXXZNwtwzCM/ZZcGAcAVX0QeDDrfhiGYRj5mVYyDMMwcoQZB8MwDKMMMw6GYRhGGWYcDMMwjDKkEQkv0kBEtgLrKvxrJPB6jCYGAFsia2VXz+TIVz2TI1/1TI76641X1ejQA6payBfwVJXyd2J+viPn9UyOfNUzOfJVz+Sos161e2fpa1+cVvogZr0Hcl7P5MhXPZMjX/VMju7Vi6TI00pPqWpZmqRq5UXD5MgXJke+MDka/51FHjl01FheNEyOfGFy5AuTo8HfWdiRg2EYhtE4ijxyMAzDMBpE7o2DiCwSkbdF5LlA2SQReUJEnhWRB0SkvyvvKSKLXflqEflC4DNtrny9iNwkKSfzTVCOZSKyTkRWudfglOUYISL/FJG1IrJGRC535QeJyFIRecn9HRT4zNXuvK8TkVMC5ZnpJGE5MtNJrXKIyMGu/jYRuaWkrcLoI0KOIunjJBFZ6c77ShH5YqCtTO9Zke5MWb+A44GpwHOBshXACe54HnCtO74MWOyOBwMrgSb3/klgGl5iob8CswoqxzKgPUN9DAWmuuN+wIvABODnwHxXPh9Y6I4nAKuBXsAY4GWgOWudJCxHZjqpQ44+wHTgUuCWkraKpI8wOYqkjynAp9zxRODNPOhDtQCurKq6HHi/pHg8sNwdLwXOdMcTgH+4z72N5ybWLiJDgf6q+oR6Z/0u4PRG9z1IEnKk0M1IVHWTqv7XHW8F1uKleZ0N3Omq3Unn+Z0N3KOqO1X1VWA9cHTWOklKjrT6W41a5VDV7ar6GPBRsJ2i6aOaHFlThxxPq6qf2GwNcICI9MpaH1CAaaUqPAd8xR2fTWcWudXAbBFpEZExQJv73zC8bHM+G11Z1tQqh89iN1z+YepDzQAiMhrvyec/wKGqugm8HwjeiAcq5wcfRo500k05fDLXSUw5qlE0fURRRH2cCTytqjvJgT6KahzmAZeJyEq8odvHrnwR3kl8CrgBeBzYTcwc1RlQqxwA56vqZ4EZ7nVBqj12iEhf4D7gClX9MKxqhTINKU+VBOSAHOikBjmqNlGhLM/6CKNw+hCRI4GFwDf9ogrVUtVHIY2Dqr6gqierahtwN978L6q6W1WvVNXJqjobGAi8hHejHR5oomKO6rSpQw5U9U33dyvwezKY2hCRHngX/u9U9Y+u+H9uKOxPUbztyqvlB89cJwnJkblOapSjGkXTR1WKpg8RGQ78CZijqi+74sz1UUjj4HsfiEgT8APgdve+t4j0cccnAbtV9Xk3jNsqIse6IeYc4P5set9JrXK4aaZDXHkP4DS8qak0+yzAHcBaVb0+8K8lwFx3PJfO87sEOM/No44BxgFPZq2TpOTIWid1yFGRAuqjWjuF0oeIDAT+Alytqv/yK2etD78TuX7hPVFvAnbhWdNvAJfjeQG8CCygczPfaLxIrWuBvwOjAu20410kLwO3+J8pkhx4HhorgWfwFq9uxHnMpCjHdLzh7TPAKvf6MnAw3iL6S+7vQYHPXOPO+zoCHhdZ6iQpObLWSZ1ybMBzjtjmrsUJBdVHmRxF0wfeQ+H2QN1VwOCs9aGqtkPaMAzDKKeQ00qGYRhGYzHjYBiGYZRhxsEwDMMow4yDYRiGUYYZB8MwDKMMMw6G0QBE5FIRmVND/dESiNhrGFnTknUHDGNfQ0RaVPX2rPthGN3BjINhVMAFTXsIL2jaFLyNinOAzwDXA32Bd4ELVXWTiCzDi4F1HLBERPoB21T1FyIyGW/3e2+8DU3zVHWziLThxdHaATyWnnSGEY1NKxlGdcYDHap6FPAhXp6Nm4Gz1IuHtQj4WaD+QFU9QVWvK2nnLuB7rp1ngR+78sXAd1R1WiOFMIx6sJGDYVTnDe2Md/Nb4Pt4CVmWuijQzXghUXzuLW1ARAbgGY1HXNGdwB8qlP8GmJW8CIZRH2YcDKM6pbFltgJrQp70t9fQtlRo3zByg00rGUZ1RoqIbwi+BvwbaPXLRKSHi8NfFVXdAmwWkRmu6ALgEVX9ANgiItNd+fnJd98w6sdGDoZRnbXAXBH5JV40zZuBh4Gb3LRQC14ypjUR7cwFbheR3sArwEWu/CJgkYjscO0aRm6wqKyGUQHnrfRnVZ2YcVcMIxNsWskwDMMow0YOhmEYRhk2cjAMwzDKMONgGIZhlGHGwTAMwyjDjINhGIZRhhkHwzAMowwzDoZhGEYZ/wdypNwFR9m0QwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 12, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXl8JGd95/95qqqrb90ajaTR3OOxPePx2B4fmMvGBBsScNYB1iQEhyU/Z40Tkl9OyLGQvOIQICG7kIXFCSyGJIADIZjEF76wgfEx43s8h+fSSKP7aqmvqq6qZ/+oeqqru6sPSX1Udz/v12te0pS6WtWl7uf7fL4noZSCw+FwOBwnQqMvgMPhcDjegxsHDofD4RTAjQOHw+FwCuDGgcPhcDgFcOPA4XA4nAK4ceBwOBxOAdw4cDgcDqcAbhw4HA6HUwA3DhwOh8MpQGr0BayVvr4+unXr1kZfBofD4TQVhw8fnqOU9pd7XNMah61bt+LQoUONvgwOh8NpKggho5U8jruVOBwOh1MANw4cDofDKYAbBw6Hw+EUwI0Dh8PhcArgxoHD4XA4BXDjwOFwOJwCuHHgcDgcTgHcOHA4ZXh5fAkvjS01+jI4nLrCjQOHU4a/uv8o7rr/aKMvg8OpK01bIc3h1IuUqkOntNGXweHUFW4cOJwyKJqBjG40+jI4nLrCjQOHUwZVM6Bo3Dhw2gsec+BwyqBoBtIZvdGXweHUFW4cOJwyKJqOFDcOnDaDGwcOpwxKxkAqo4PyoDSnjeDGgcMpg6IZoBRQeVCa00Zw48DhlMAwqG0U0io3Dpz2gRsHDqcETrXA4w6cdqKscSCEjBBCHieEHCWEHCGE/LZ1vIcQ8iNCyOvW127HOZ8ghJwkhBwnhNzoOH4FIeQV62dfIIQQ67ifEPId6/gzhJCt1X+pHM7qUTLcOHDak0qUgwbg9yilFwG4BsCdhJCLAXwcwKOU0l0AHrX+D+tntwLYA+AmAF8ihIjWc30ZwO0Adln/brKOfwTAIqV0J4C/A/CZKrw2DmfdKFrWIPB0Vk47UdY4UEonKaXPW9+vADgKYBjAzQDusR52D4BftL6/GcC3KaUKpfQMgJMAriKEDALooJQepGbaxzfyzmHP9V0ANzBVweE0EmfxG1cOnHZiVTEHy91zGYBnAAxQSicB04AA2GA9bBjAmOO0cevYsPV9/vGccyilGoAYgN7VXBuHUwtylIPKjQOnfajYOBBCIgC+B+B3KKXLpR7qcoyWOF7qnPxruJ0QcogQcmh2drbcJXM468apHNIaNw6c9qEi40AI8cE0DP9MKf036/C05SqC9XXGOj4OYMRx+iYAE9bxTS7Hc84hhEgAOgEs5F8HpfRuSukBSumB/v7+Si6dw1kXOW4lnsrKaSMqyVYiAL4K4Cil9POOH90H4Dbr+9sA/MBx/FYrA2kbzMDzs5braYUQco31nB/KO4c913sBPEZ5OSrHA/BsJU67UklX1jcC+FUArxBCXrSO/TGAvwZwLyHkIwDOAXgfAFBKjxBC7gXwGsxMpzsppexTdQeArwMIAnjA+geYxuebhJCTMBXDret8XRxOVXDGHLhx4LQTZY0DpfQncI8JAMANRc65C8BdLscPAdjrcjwNy7hwOF7C6VZSuHHgtBG8QprDKUFuzKEy46AbtOLHcjhehRsHDqcETrVQqVvpn54exXV/83itLonDqQvcOHA4JVhLEdzp2TimlxVovIsrp4nhxoHDKQEzDrIkIJ2pbLFfSmUA8BbfnOaGGwcOpwQsW6kr6Ku4t1KMGQeH6hidT2A5nan+BXI4NYIbhxbkv3/zML78xKlGX0ZLwOocOoK+ioPMS8lc46BoOt79xZ/g7x87WZuL5HBqQCV1Dpwm49DoAnReQ1gVVN2ALAoIy2LF7TOYcmAuqadPL2A5rWFiKVWz6+Rwqg1XDi1IQtHtBYqzPpSMAb8kIOATK1YOsbyYwyOvTecc53CaAW4cWgzdoEhldCzzhagqKJoOv09AUBYrijlQSrPKIWOAUopHj5rGYTGp1vRaOZxqwo1Di5FUNQDuu9SfnZzDwVPz9b6kpkbRDPglEQFJrCiVNa5o0A3TpafqBl6bXMZELI2AT8BightsTvPAjUOLkVDMBYwFRRmKpuO3vvUCPvfQsUZcVtNiGgemHMqnpjrvu6oZeOa02Vz4hosGuFuJ01Rw49BixBVTOaQyek4q5f2vTGI+oWIxyReo1aBkdMgs5lCBcnAaAFUz7L/Hzv4I4oqW8zfhcLwMNw4tBnMrAbkL1TcPjgLgfu/VwpRDwCdUNAkuxzjoOhRNhygQ9EX9AIClFL//nOaAG4cWg+1UgexCdXxqBc+fW0JfREYslbF94pzyKJoOvyQi6KsslTVfOaQzBgKSgK6gz/w5V26cJoEbhxYjqWQXMLZQnZlLAACu270BlPKUytWgaIaZreQTkdEpMmVaYjhjDopmWNlOIrpDMgBwtx6naeDGocVIONxKLJ11xWrbsKUnBIC7llYDq3MIyiIAlE1ndRpexakcQqZy4Pee0yzwCukWI+FQDsy/zVxNI5ZxWOILVMWoupnK6veZxiGV0REN+AoeNxdXMLmUzokpqJphKQ/RNg783nOaBa4cWoyEM+aQZMqBGYcgAPB8+1VgxhxMtxKQO1PayRcffR3/9e6DmFtRIYvmx8qMOZjnM7dSfooxh+NVuHFoMXID0pp9LOAT0B8JAOCujdWgZLIxB6D4TIfRhSSSqo6Dp+bQb2UmqbqpHAI+ESFZhCwKPObAaRq4cWgxkqqGoE9EWBZt//dKOoOI34fuMPd7V8ozp+fx6vmYXSEdlM2PSrH+Sqyp3kQsnTUODuVACEFXyMfdSpymgcccWoy4oiPsN3epWeOgoSMgIeKXIAmE714BvDIewz89PYpP33IJBIEU/PxTP3wN3SGf7VYKSMUD0pRSnF/MdlztDcsQBWLHHFgaq2kc+L3nNAdcObQYSVVD2C+hMyTnGIdIQLJ2rzLfvQJ46uQsvnNorGha73IqgzNzCbsILuQ391FOt132sRoSqg5mYzpDPsiiAEXToVjKAQC6QjJXbZymgRuHFiOhaAjLEjqDEmKObKVowFzcukM+HpAGkNHMQkBn6q+TuKJhMpYGpYDfJ2Ko04zXnHeZycCOXbujDwDQGfRBlgTbrRSw4hXdXDlwmghuHFqMhOVW6gz68mIOzDi01u41rmj43e+8iJmV9KrOY8VsSZcYAqU0J+vLLwnoj/rhlwScm08WPJ7FG95z6RAAoCsom8bBCkjbyiEo8/YZnKaBG4cWI8HcSg7jEE9riPhNv3d3OLt7pZTi4SNTFc9G9iKvTSzj3144jx++NLmq85hxSLi4iRTNgOZoMcICyiM9IYwtuhiHmGkcrtvdj4+8aRvesWfAcivlKoeusA+LyQwon9LHaQK4cWgx4rZbyaEcctxKMhYs5fDc2UXc/s3DePToTMOud70oVr+jn56cW9V5agnlkB9XkK2d/0h3EGMLWbfSQ0em8Hc/OoHzSynIkoC+iB9/9gsX46LBDvgtt5JTOXSHZKiaUVF3Vw6n0XDj0GIkHW6ldMbcuTpjDiwgTSnFkydmARTvtfTo0Wnc/o1Ddbv2tcCK0p45PV+275GTUsoh/5jfylTa3BPC2ELS3vnf99IEvvDY63h+dBFDnYGcrCdZKlQOnaz5Hu9txWkCuHFoMRJK1q0EwA6qOgPSGZ0ioep46vVZ+xw3nju7iIdfm/Z0F1emABKqjpfHlyo+r1RAmlWUM9jOf6QnhBVFy8kCo9S8T0NdwYJzkqoGg2bPD/jMr8WqrDkcL8GNQwtBKTVjDrKELqtdw+i82ZHVjjlYx8/MJvDy+RgA9/RMILu79vKAGsXRRvunJysfgZpVDoUuHmYsh60F3+/LGgcAtmuJNTQEUGAcZEmwjQxTDkyBqKtQOBxOo+DGoYVIZwwYFAj7JQx1mamXJ6ZXADiUQ9g0Dve9dB4sLlpMOWjWIqZUMMegUbBdeF/Ej0ePTlcc7M3GHApfOzOWe4c7AGQX9ZFu0zicWzCD0itpza5tGHYxDqwrLjMurOcSVw6cZoAbhxaCuUjCfhHDXeZCdmzKNA4RyzgMdJitHf7hqTPoCEjoCctIFGkJoermQqt4WjmY1/ahN2zBS+MxPHRkqqLzSikHZhwuGe4EkN35s8aFLGNpJZ3Bdbs3oCvkw75NnTnPIYsO5WAZF2YkvGxsORwGb5/RQjAFEJYlbIj64RMJjk2axqHDMg6XDHfiHz90AIdGF3Hhxij+5yMncpTD2EIS5xaSeOPOPnsB9fJOly20t127Ffe/Mom//M+juG73BntBL0bGMnyllMPN+4fRGZKxf6QLABAN+NAd8tnKIZ7WsK0vjH/80IGCFhyyJGA5nascmALxsrHlcBhcObQQbFEL+yUIAsFQVxAnZ+IAsjEHQgjefvEAPv7OC/GLlw0j7JdyjMNXf3IGv/WtFwBkd9de3ukywxWWRXz8nRdifDGFJ46XT83NOALZ+bD70RXy4Vev2QLRsfCPWBlLumEG9aMBybU3kyyJtgFiRoEFpr18PzkcBjcOLcLPTs7h9WnTEIT95mK0qTto+9ZZzCGfsCzlZOwkVc3eTWeNg3d3uopmQBQIJFHAroEoALPXUTlYkD3pEm+JW66msFx4z4a7gpiMpRG3XEZug3+AbHwByGYp+Xm2EqeJ4G6lFuGOf37e3pGGrVYZziBppJhx8IuYi2dbOiiaYe94VY3FHLy701U03V6IQ5Yryc1VlA8zfHG3mENaQ1gWXRVBT1jGYkK1XUZRv/t9ZYYAcCoH7lbiNA9cObQAmm4glsogbe1IWR+lTVZ2DeC+CwZQ4FZKZ3ToBoVu0CaJORj2QszmPCcrqEAuFXNIKFpRY9oTNntTsVqHYorMVTlI2QlxHI7X4cahiVhJZ3D93zxRkJHDsmIka6cbshZJphwifinHb+4k4pdy6hzYrjajG83hVspk21OYPZCAdJHsKyelYg5xq5DQje6QDINmO7EWcyuxazK/N/8eMo85cJqItjYO6YyOH1stJJqBH7w4gTNziYJKYLaL/dgNu/D777jANgqburPGoRghWcrpL8RUguleygak44qGw6OL1XsxVcIcxmMuvoQQhHyia7+kfOw6B9eYg1bUXdQbMetEWHfWospBKq4cvGxsORxGWeNACPkaIWSGEPKq49inCCHnCSEvWv/e5fjZJwghJwkhxwkhNzqOX0EIecX62RcIIcQ67ieEfMc6/gwhZGt1X2JxvvqTM7jta89iLq7U61euGUop/uWZcwCAhbx5DMz/fdFgB37zbbtg3VoMW8ah2AIGABG/iISq2cVjaWtXm9GNnDqHe58bw/u/ctBzfYFU3cjZpQdlsUK3UvHGe4kyygEARhfMyvNK3Ep+H485cJqPSpTD1wHc5HL87yil+61/9wMAIeRiALcC2GOd8yVCCEs4/zKA2wHssv6x5/wIgEVK6U4AfwfgM2t8LavmkaPTAIrPBfYSr5yP4bXJZQDAQiLXmLEFm/VTYmzsCEAUSFH/OWDGHCjNLpJMOaiagYyWjTkspTLQDYpxl5bVjUTJGDnB36AsVuZWKtFbqZRbqSfM2pKY96HYvc1RDtb3PpGAEEDhXVk5TUBZ40ApfRLAQoXPdzOAb1NKFUrpGQAnAVxFCBkE0EEpPUjNLeo3APyi45x7rO+/C+AGpipqyXxcwYtjpntG83BjOca/HhpHwCdgz1AHFhK5A2NY6mZHMHehkkQBGzsCRf3iAOzxl2yRVBzKwelWYjMfxhcLJ6E1ErMldrbgLVihW8lWDkUqpIu5lfKNQ0exVFapUDkQQuC3urVyOF5nPTGH3ySEvGy5nbqtY8MAxhyPGbeODVvf5x/POYdSqgGIAeh1+4WEkNsJIYcIIYdmZ9cXK3ji+KzdW2g1rZ4bgWFQPHRkCtfv3oCtveEC41BMOQDAR6/fgQ9cOVL0uSNWTQRrI5F2KgdHQJpl9Zz3nHHQ89xKUkVuJRZzUHWjIHuoXEAaMAPSPpHk/G4nToMVyAtOc+PAaQbWahy+DGAHgP0AJgH8rXXcbcdPSxwvdU7hQUrvppQeoJQe6O/vX90V5/HYsWwVrdeNwwtjS5hZUXDjno3oCcuFysGKObjtYn/l6i145yWDRZ+bpbiydFa2cKl6tt5B0QykVPO42wzlRuIcpgOYtQ6VZiuxuIAznZWNCC3mLgrKIoI+EbpBEQ34UEzkMuXACvScx7lx4DQDazIOlNJpSqlOKTUA/AOAq6wfjQNwblM3AZiwjm9yOZ5zDiFEAtCJyt1Ya+bp0/PoszJPNN3bbqWHj0zBJxJcf+EGdIdl2//PiKUykARip7CuBpbJFFdy3UqqZti7aza0BoA3Yw5Ot5IsIpmppAiOojNkGlNnOisrAiyV4cVcS6UC/cw4BPKUhelW4jEHjvdZk3GwYgiM/wKAZTLdB+BWKwNpG8zA87OU0kkAK4SQa6x4wocA/MBxzm3W9+8F8Bit8ZBdSikWkyqGrSIxzfDuTo5SigePTOENO/rQGfShNyyDUmApmVUPsVQGncHiu9hSsJgD2z2zgHRGpzkxB9ut5DnloBcEpMvFHFiRX5flhnOmszIFVYlxKPUYpkryGwDymAOnWSjbPoMQ8i0A1wHoI4SMA/gkgOsIIfthun/OAvgNAKCUHiGE3AvgNQAagDsppeyTegfMzKcggAesfwDwVQDfJISchKkYbq3GCytFQtVhUKDH2jmyNhFeZDauYHQ+idvesBVAdh7DQkJFb8Rsv72cyqDDJd5QCSzmEFd0GAbN+uLzspXY3GPvxRyMnLTRStxKzOix+IFTObDYS7GYA5D9G5RSDs7CvNzjoqcrzjkcRlnjQCn9gMvhr5Z4/F0A7nI5fgjAXpfjaQDvK3cd1YRN8OoJm4url5VD2vL1s2Bzr8M4MGLrMA5sEUwoWs6EskxBzMFcNBeTmZJ1APXG2T4DqKzOgRkH5lZyKocVxXxvlFIFvbZxKH7PbbdSvnLwcbcSpzloywpp1m6itwliDmzB9lmLTY+LcVhOa/a8htXiNA7pTK7vXXW4lVIZ3Z565iXXkpLRC2MOZZWD+fdmbiU35VDKODDFUUnMQXaNOXh3M8LhMNrUOJi7Q/Yh9/JMX5ZmKYvmysyMw7zTOFgxh7XAOpkmFD1n0UplcvstJVUdm60Zyl5yLRVmK0lQNSMnYJ8PUw5dLCDtUA5xSzmwtudu9ITN84rVOABZd1JhzIGnsnKag7Y0DqxojH3Ivawc2ELGdqDMoC1WyThIooCAT0BC1XJ84c5W1krGzFbaucGcl+CVjCXdoNAMmqcczPuUKuFaYga3y7qXr56P4VP3HbG72wIo6aZj7shKlEN+zEGWBN6VldMUeMNxXGeWmyjmwIyDT8y6KaJ+yVYOlNJ1xRyAbGfWtMMXnszp1KojZSkHWRQw7hG3Eltkc2MO2eyrYq4hO+Zg3bOv/fQMDAr8ytWbMRUzW5MMdASK/l62qVh7thKPOXC8T1sqBxZzYB/yjIeVA1sAfY6MnJ5IthAuqerQDLpm5QCYcYekkqscnK4WJWMgmdER9ovoj/oxu+yNRoVskc2pkLYWYxbId4P9vdk9Yx6oiVgaU7EUogGpwphD+YA0z1biNCttbhxM5eDlCmk1z60EZAfOAKWroyslJEuIK3rOjtYZpF1OZ0CpuQvui8iY9UgXW+a7d7qVQvbAn+KFcOzvHfCZ1c5slz+xlMJELI2hzmDRcwFzjrQoELsluhuls5W8+37jcBhtahzMimK2O9Q8bBzYLteZy98TkjFvjfYs1VepUiJ+0cpWclcOS0nzd4RkEX0Rv/27Gw3bgee37AbcW3Ez7AwwkeBNu/rwhzfthigQTCylMBVLY2NncZcSAAx1BfGzj78Nb97VV/QxfpG16W7dCumJpRQ+8+AxGE3QuNJr1LjOtyq0pXFYTmcQDUj2gtsMbqWiysEKrq/XrZRQNVflQAiwlDJ/V9BnGgevzL+w3Uo+N7dS4QI8FUvj4SNTdnGfLAr4hw8dwK+/eTs2dgRwfimFyVgKQ12ljQNgxiRKVaQXVQ4tlK30yNFpfPmJU55KbfY6lFJ88dHXceVdj9hZk16lLY3DSlpDR9AHyUoP9bJbKT8gDQC9ET+ml9P44++/gmNT5oyH/HbdqyEsm3OknYuW3UZClmxFEZRF9FrxDi/sFku6lVyMwxcfex13/PPz9nk+h8Ed7Azg7FwCc3EVGztKu5UqoXjMwcxWaoadYzny+3FxyvPFx07ib390AnNxFZOxdKMvpyRtma20ktYQDUj2guvleQ7ZgHR2l/rBazZjZiWNf3t+3P75epRDNCBhJZ1bBGcbh4CEFet7phw0w8yQYm0kGgVblJyqKhtzKFywDp1dhG5dO5BrcIe6gnjQms09WIFyKIcoEAx1BrC5N5RznF2rOcFu9Y0SvQR7j6R5gL0iDIPiS0+cxHBXEOeXUjnp6F6kTZVDBlG/z15wvawc3ALSm7pD+Pz79+Pe33gDQlbq5noC0h1BH2KpTK5yUAsb0IVkya4qn0803rXkFnMIFHErxZIZHJ9eAZBtWug0uENdQdvQDpaJOVTK439wHT549ZacY600Rzo7A4Qrh0qYWVGQzhh4ywVmrGrJYyN382lL47CcMpUDIQSiQJqjCE4s/FPt29SFb3zkKnz0uh12te9a6Az6oGgGlq03qygQe0Kac65BUBbQbzX7m11p/K4n61ZyKofcLrOMw+eyXeAXrQC7854OO9RCtYyDXxIhCLlxCTYVrhXSWZlbqVTBISfL6Lw5d/zSTV0AzA2Ll2lTt1LGzlGXBOJt5eASkHZy+eZuXL652/VnlcIK6GZXTDUQDUj2B9+pHAI+0b5vXghKZ+sc3FJZcxesQ2cX7e9Z9lW+W4mxsUwq63rIKofmX1Djae5WWg1stOylI6ZxWEw2foNVirZUDmZA2lz0ZFHwdLaSW0C62rB4xbRlHCJ+yQ7oOltEhGTJ7kg67wnjUFgh7ZcEEFLoVjp0dhGitYu33UpSoXEoVwC3XlrKraQy49D8hq4ejC4kIAkEuzZE4BMJdyt5DcOgiKtaVjmIxNPtM1TLcEnC6gf5VArr6DqznIZPJAj6RNeYQ9AnojskQyDAnAdqHdzcSoSY1+/MVlI1Ay+OL+EKS2EtusUcLLVQrgBuvTCV00puJW4cKmN0Ponh7iAkUUBXSM4Z2OVF2s44xFUNlGYXREkUPO9WkiVhTVPeKoUph9m4Ar8kwicKYJmWEX82lhGUTR96T9jf0IC0YVD87NQclEyhWwkwjZjTrXRmLgFVM3DN9h4A7jGHjqCEsCyWLYBbL63kVrKzlVpABdWD0fmk3dm4K+iz3Ztepe2MAwu6MndJM7iV3ILR1cQ2DssK/JKQE9/ICUhbwdS+iNzQgPTTZ+bxy//wDH52ah5ArlsJMI2Y063EspQu22Iqh2y2Uq7ieM/+Ybz94oGaXjszDq3QmZVlKylcOVTE6HwCW6zU5u6Q7PmYQ9sFpFlfpQ6nW8njysHp/qgFzDisKBqGg8EcYxS13EqSQGyj0R9tbJU0y/J4eTwGoLDQLJQ38OfE1ApEgeCS4U4AWeWQH8f59C2X1OyaGcyQtULMgbuVKmcpqWI5rWFrbxiAOYVwbMEbre+L0XbKgRmHnGwlDxfBZXSjaKZStXC2+y6mHIKONhC9YbmhbiXWWpy1bchXVvlupePTK9jWF7a7qbJmhbU2um7YMYcmNw6UUl4EtwrOWplK3K3kYVg/E+ZW8omC3WvHi6i6UdNMJcC8BywFVJaEnEWTBaQDctY49EX8mGugW8m5GPld4jH5bqUT0yvYPRCFKJjBakpNw1DLOE4xWiXmoGiG3VmAK4fysBqHLZZy6A7Lds8yr9KGxoEph6xx8Hr7jForByDrWvL7RFflEHIah6gfqYxeUGhWL5yLkdu9CcmS3bI7peo4t5DEBQPmFDs2M7vWBrcYrZKt5Ozam25yQ1drYqkMvvT4KUQDkh1z6Az6kLYmLHqVtjMOy7ZyyMYcvJytVI+ANJA1DgFJyFk4WczB6VbqtqqxFxrUG8ZZkevWn8iZynpyJg5Kgd0bIwDM9uRAA41Di8QcEo4xsqkSg5XagefPLeIzDx7D/a9MuiYafOxbL+D0XBxf/pUr7PYuzMXpZddS2xkH9qa2lYPg7VTWjE7rspCxAL1TOQgkOx8hKDurkM17lyoxM6GW5LuV8gnJon1tLFOJKQd27Y0yDszQN7tbaUXJLmrtrhz+8anT+PITp/DRf34e//b8eM7PdIPixydm8eE3bsObHPM/WLsbL2cstZ1xuOO6HTj+lzfZFtwnebu3Ur3cSiwo7ZcEewHziYK9M3cqh7C//ECdWuJMncxPYwWs+RSW2+P1mRXIomD7elkMRW5AMBrIXm+zp7I6lUO7p7KeX0rjmu098IkEo3kZSKyYdEPUn3OcGQeuHDyG0xUhCYKns5XMgHTtFzLbreRQDqZxML93xhyCPnOBTTQw5pCdl1DoVmKprJRSLCZU9IRlu3UGM2y+OhhcN7LKodmNg/m3FwXS9tlKk0spbO4JYaAjgMm8wUes/1R+S5auoOlWink4KN12dQ75+Dxe55DRjZr2+mGwXlN+R8zBJxJ7pxtwUw5K49xK3SEfKHV3K4X9EjSDQtEMJBTdvl72M6BxbiVJFCAJpOndSqzGoScsezqoWmtUzcBsXMFgZxBDnUFM5A3wYUY0nPcZ7g4ztxJXDp5F8ljM4ZnT83j06LT9f1UzXBfAatPpdCtJpd1KdlvsBi0KqYyOoE/EZZu7XEd6hh3T4FYULce4hhsccwCsOdJNvttmi15vWG7rmMP0chqUAkNdAQx2BTAZy1UOK46hWU6YcnBzK33yB6/ij7//So2uuHK4cpAEz8QcFE3Hx779AtIZA8/9ydshS6bhqsdCljUOokM5uLuVssqhcW6lgE/E/7r1MriVKrBdWkLRkFC0nF1buMExB8AM+je7W4kph76IH1PL3h53WUvYqM/BziA2dgYwFUvDMKg9xyPh0voeAAI+cxPm1nzv0OiiJ2KEnNPmAAAgAElEQVRSba8cfAJBpgFdWdMZvcC18O8vnMf0soJYKoOfnpoDUP86h4AvaxBkKfu9swguZMccGuRW0gz4fSICPtE15mAbB7XQODQ6lRUw4w7N7lZiAeneSHu7lZhSGOoKYKgziIxOMefoHlAs5kAIQXfIh1mXNjSLCdU+Tqk51nZsIVn3ueNtbxzM3kr1Vw6/fs8h/Mn3X7X/rxsUX3nyNC7cGEXUL+H+lycBNCCVVRLtADirIv6DG3fj3fuG7MeytNZUAwPSgRIGk6mchKIjnu9WanDMATAX1LGFVPkHepiEqiHgExCSpbYOSE8sZZUDmyA4uZRVUm5DsxhXbevFA69MYTpPeS0mM1hKZpDRDfzV/Udx6Z8/jDd/9nE8+fpcrV6GK21vHHwNatl9Zi6Bl8aW7P+/ej6G07MJ3P6W7Xj7xQN4+LVpZHSjLu0zALMRGGCmWjpTWQHgzut3Yq/VtA4wFYUsCo1TDhk9p+4in0iBW8mhephxaFC2EgBcu6MXh0cXG1ZhXg2Y0Q34hLZOZZ2MpdARkBD2S/bAKGfcoZRx+P13XADdoPjbh4/bx1Kqbhd5zsdVHB5dRJ81t51NaqwX3Dg0qGX3cjqD0fkkdCuNlk2F2tIbxrsuGUQslcGzZxYaEpD2SbnGwY2gLDY25uDiTmI450ib2UqFbqVGxhzeckE/VN3A06fnG3YN64W56wI+sa1nSE8spW2jwJTDhEM5FMtWAszP+m3XbsG/Hh7HhJUC6yyKm4srmIqlsXujWcBZ701s2xsHSah/KqtuUKykNai6Yb8p4o6eTxdab4aJpZQVkK79QjYQDSDoEzHUlW3ZXaptRzivLXY9SWcMBFyK3xhsl7aUzEDVDURkb2UrXbm1BwGfgCdP1NdNUE3iaQ1hWUJAEqEZ1NPp4LVkMpayjUJPWIZfEnKUw4qimUq7yAbvrRdsAKWw23c7jcNULI3pFQWbe8wCznoHqblxaIByYJ1hAdO9BABxqx1B2J+dYRxXtPplK4V8ePqPb8A7Lh7IprJKxY1SyDFnut6wbKVihCx1MGPJ8Fzl0HjjEPCJuHpbL548Mduwa1gvTrcS0L7T4CaWUhi0lAMhBIOdATuDCTCVQ6k6pV7LZcT6lC0msmvDsall6Aa1m/Vx41BnZLH+2UrLqaw75qzVynfFkdXAFrPllIaMTuuSrQSYriVCSEHMwY2QLDasQjpVxjgwdcACfV4LSAOma+n0XMJWjs1GQjVjOezv0I4ZSylVx2IygyHHaNnBzmCOcYinKzMOc5ZxWHAoBzbMaos1A0LlbqX6IlnzkvU6ttCIpbK7g9OzpnFgqYERv2SnkC66jLOsB74KjUPjeisZJY1DwCdAIO7KgQWn62Vwi7Frg9kldnyxSY2DFcsJtrFxYJuPjZ1B+9hgV24LjXhezCsf1p113kpdXXR0On71vGkcRnq4cmgIkuXPr2ewhxkHQrLKIa5kEJJFuwdQNCDZUrMeAWknbOEsHXOQGpJtoxsUql465kAIQViWMGN9eN3aZzQyIA2Yo1YBNHTc6npgbiXWXqUd01mZ2u90TFLsj/gx71jg40rGbnvvhk8U0BXyYT5uuZWsDeFwV7YVx3BX0JxY6TXlQAj5GiFkhhDyquNYDyHkR4SQ162v3Y6ffYIQcpIQcpwQcqPj+BWEkFesn32BWGO4CCF+Qsh3rOPPEEK2VvclloYtgPW88WymxI7+iCPmkF+sJXlAORRfQM1spfrvFlnxWCnlAJhGgCkHL7qVmHGod3piNdANigWroWE7u5WYWzXsSKuOBiQommG/T/N7e7nhHLu7mFDRGfRho+WqCvhM4yFLgieVw9cB3JR37OMAHqWU7gLwqPV/EEIuBnArgD3WOV8ihLA782UAtwPYZf1jz/kRAIuU0p0A/g7AZ9b6YtaCZO3U61kIx5TD/pEujC+moGoG4oqes8OIBCR7N1HvhUyuIJXVVA71XxDYnIZgGeMQ8ov2wpvjVrLiEVKDjUN3SIZAmlM5zMUV6AbFYFfQNg7jiyn8jx+82vSV36uBKeeQ4/3FWt8zVRFXNEQCvsKTHfRG/PZnfSGZQXfIh/6IuXkY7AyCENKQeqyynxBK6ZMAFvIO3wzgHuv7ewD8ouP4tymlCqX0DICTAK4ihAwC6KCUHqRmDfg38s5hz/VdADcwVVEP2CJRz6D0smUcLh3pgm5QjC0mEU9ncppzRfxZt1K9/eN2QLpUFbK/MQFplhVTyq0EmPePjX91KgdRIPidt+/CO/durN1FVoAoEPRG/E2pHFgQfagzYFeq//ClCXzj4CiOTa408tLqCosT5isHIPsZN91vlSgH87O+lFTRHZbRFzVjERs7TAUhS0LTBKQHKKWTAGB93WAdHwYw5njcuHVs2Po+/3jOOZRSDUAMQO8ar2vV+OyYQ32VgygQOyg5FUubbiVHPn404LMzF+pR5+CkkpiDc2ZCPWHui3JupdxGgbk+3995+wW4dKSr+he3Svqa1DiwbJyNnQH773BsahlAbrZNq+OmHKL+POVQJlsJMDOWWEB6IaGiJySjjykHq+uwLApQtebureS2itESx0udU/jkhNxOCDlECDk0O1udHHHmOqlnEU8slUFn0IfesJWpkFARV/Qc5RD1S7aPsd4BaWdvpWKEZMkODtcTZhzcGu45cRracj7fRtEf9XvKrbSYUO1irFIw4zDUmXUrnZ1P2s/RLjC3qlM5ON1KukGRypTOVgKA3rAfi8kMNN3AYsJUDiwmxQrsmkk5TFuuIlhfZ6zj4wBGHI/bBGDCOr7J5XjOOYQQCUAnCt1YAABK6d2U0gOU0gP9/f1rvPRcbLdSHZXDclpDR0BCj2UcFhNqQVaD01B4M+bQmIE/zDiU6q0EOAPPpKwhaRR9EdlTyuGzDx3Dh7/+XNnHTS6l7EApc++xVPCFNjQOoRzFb7mV0pmSfZWcsN5JC0kVC0kV3SGfrRxYmqxPJMh4MCDtxn0AbrO+vw3ADxzHb7UykLbBDDw/a7meVggh11jxhA/lncOe670AHqN19FX4WEC6jjEHphy6QjIIsZRDujBbyb7GehuHSuocHG2x6wlLmSzVlRXIqoV6TNFbK6ZyUOvumivG+GKqIiUzGUvbgdL8xIB2Mg4JRYNPJDkxwaxyyBSd5ZBPT9g0BOcXU+aUw7CMzVZtwzZr9nkjlEPZTw4h5FsArgPQRwgZB/BJAH8N4F5CyEcAnAPwPgCglB4hhNwL4DUAGoA7KaVsa3kHzMynIIAHrH8A8FUA3ySEnISpGG6tyiurELYAZuroz1tOZdAR9EEUCLqCPiwkFCurwV051D0g7ZjnUIyQ3ba7McqhbCqrtZsrJ+kbSX/ED1U3sJzS7K64jWQurtoLWimc/YT8eX+HxbaKOeg5qgFwBqS1rHIIlI85AMDJmTgAoCck46LBDtz/sTfjokGzz5pPrH8qa9lPDqX0A0V+dEORx98F4C6X44cA7HU5noZlXBqBXQRX52yl4W5TLvaEZUzFFGR0mrPDiDZQOQRlEdGAhA2W39MNtvjWu223rRzKprKa1+d15QAAs3HFE8ZhPm6+DxVNL+mKm4ylce2OPgCFWWPtphzCee7NiCyBEFM5xEt0ZHXSl2ccuqyq6YuHOuzHyGLzxBxahmxAur7ZSqyqsics49yCWQgXLaYc6mwc/JKIp/7wetxy+aaij2HKYXQ+gU/dd6RuuxrWHrpcnQP70HpdOQDeKIQzrMI2oHQcSdMNzKwo9uxuWRTsUa0bOwJtZRySqp6TqQQAgkAQ8UtYTmvZTssVBKQB4LmzZqiVGQsnXi2Ca2lYEVy9CkwopVhO5xsHM9PDmWET8Wd3ko3oA9QVku1WHm4wOf3tZ8fw9Z+dxYnp+uS3Z91K5WIO3ncr9XmohcZyOmPXhcRdXEuGQXHT/3wSn3nwGHSD2hW8hBAEJLPty97hzrYyDgm1UDkA5lTFZUfModx7sNNyMT9/bglbekPY75JmLTegCM67n5w6wQq96nXjUxkdGZ3aYzl7wrLtKskvgrOvscF9gNxgbbGfP7cIINsSpNbYqaxl22ewgLQ3M5UAbykHp4FySzI4PRfHsakVHJsyNwFDjmZzAZ+AvqiMgQ6//X5oB5JKYcwBMD0AK2kNKxUGpAWBoDskYy6u4M7rd7pW7zdCOXDjINTXrcTadTuVA8MpP6MNDEhXAlM5ivWGZUU/tUapsELaDki7fHi9QmfQB0kgrkPm681cPLvjdwtKP39uKef/rDgLMOM/W3rC6AnLWEqq0A1aUnW2CglVswPzTjoCvlVlKwHAQIcfIVnEf7ls2PXnjWif4d1PTp1gAel6pbKyvkpZ45AN+hZLZa13zKES8usM6mUcUqoOgZS/J83gVhIEgr6IH3MeUA7zOcahMObw4tgSogEJv7BvEPceGrdHYwLAey4dws4NEcQVDQY1Ey66w4V+81bDLVsJMDd2k7G0HXOo5D34mV/ah4BPKJp8wpVDA2B/DLUOykHVDFu+dwTNW98TzsYWiqWyNrqDqBuhPOOwnKqfWyngE1Gu/Ra7Pi9nKwHAhg4/ppbT5R9YY1hXUMBdObxwbgn7R7rwyXfvwQeu2my7RQHgE++6CADwgxfPW8+ltoVxMOdou8Qcgj6cmFnB1HIaHQGpIuW/d7iz5M99olCXNcqJ91adOsP8+fVon/H+rxzEh772LAB35RAtphw86FbyiblzcaulHF4eX8Lv/+tLBcOX0hkd9700UXYKHCPSBMoBADZ1B3HeAwN/nG6l/IB0UtVwfGoZl410IeATsW+Te18qNrimXWodSimH5ZSGs/MJbOuPVOV3+SUBap073npv1akzUh1TWc8tJLGxI4ArtnRjW59Z+djr2GE5FzK/JDh6HHnzzxSSRWyImr7SagWknzwxi+8eHseoNQSJ8eCrU/jYt17AE8dny1ZHA1nlFS1TgNRoRrpDGF9MwajjJEI35uOKnbmXrxxeHo/BoMD+zaWbFbL4WTtkLFFKi2YrRQMS4oqG07MJbLc+5+vFJ5K6z7r35qpTR1j7jHoUmMTTGt596RC+d8e1iFqynMlvQnJdNYQQe/frxWwlwPwQ7NvUaQfgqgGrY8hPjWXN3s4vpSpSDoOdQXzuvfvw7n1DVbmuWrGpJwRVNzC90ljX0lxcsQsz8wsbWQbS/pHugvOcOHuFtTrpjAFKUVDnAJgBad2gmIylsbW3OsahmRrvtQz16sqqaDpU3SjYyfZYUjzilwr86JGAZBUZedM4fO69l+IT77rITt2rBinV/Dscn4rnHJ92+OUrMQ4A8L4DI56oPC4F66EzttBY19J8XMVgZwCSQAqUw6NHZ3DRYEdOZp0b9jzkNjAOblPgGFFHPGZbf7WUgwDdoHWddd/2xiGbrVTbm84yF/IDpEFZRNAnugZOI36fZ1UDAFyzvRc7+iPoCPqq5lZKZcz7lK8cZlacxqF13rYj1m69klbZtWQ+oaIv4kfYL+UYh5nlNJ4/t1jRcCT2Xm4H5cCqyN1iDizZBEDV3EpyneuxAG4cso33auzPK9W+tycsux6P+ivLdGg01VUO5ofueL5xWFbsXk+VKodmYLg7CEKAscXGGoe5uIK+iB8Rv4S4I5X14demQSlw457KJuf1hOW2GPhjKweXbCWncthaLeNgrVNKHdNZvb/y1JiscajtTWeLp1uHxt6I7Ho8EpA8G4x2Eg34qpbKymIOZ+YSOfOIp1fSeMOOXgx3BXPSKJsdvyRiY0fAbqHSCBRNx0paQ19Etib8ZQ39Q0emsK0vjAsGKsu66Y/67TGirYw9Ba5IthJg3otqpVI3Qjl4O5WjDogCASG1jzkw5eDWhOuDV29xPacjIDXFLrmjisqBDVDRDYrTswlcNNgBSimmlxUMdATwfz98JQIeHd6zVka6QxhvYMyBZRf1Wm4l9l5NZ3QcPDWP//ambRXHvfYMdeC+FydgGBRCC1dJ2/Oj3eocrM3LtiqpBiCrHOpZCNf2xgEwW2hkahxzSJTo7f7+K0cKjgHAR6/fmROI9SrRgA8raQ2U0nUHz9MZHX0RGXNxFcenVnDRYAdiqQxUzcCGqB8XDESrdNXeYVNPEAdPzTfs97PqaObeZO/V+YQKzaCr8ptfuqkL//zMOZydT2B7lXL8vUgp5dBhfca3VSlTCaifh8OJ930WdUASSd2Uw2qKsi4YiOLNu6ozDrWWdAQlqLpRFX9oKqPjosEOSALBi2NmP58Zq73EQEdhH5tWYKQ7hKnldI4brZ6worXesOlWYrviBctorKbaed+IWen78nisylfpLWzl4GIcOkM+dAQk+15UA+ZWqqdy4MYBrKlVbZXDSoW93ZsRFoCrRsZSUtXREfThxr0b8a1nz+HcfNJWT6WGDzUzm3tCoBQNq5RmbqWukKUcrF3xgsNoVMrO/ggCPgEvjS+Vf3ATYysHF7eSXxLx1B+9Dbdeublqvy/b5ocbh7piVh/WRzmUGxnYjDAZXY24Q1rVEfKJ+LOfvxiSQPCpHx7B9HJrKwf2umrZuvvTDxzF3U+ecv0ZSz3tCcs5qawLVr+lcvUNTiRRwN6hztZXDmpx5QBkZzRUCz9XDo1BEoSat8+IpzUIpPwEs2aEBeCqkbGUyugIyiI2dgbwsRt24bFjM3jktWkAZpO6VoR1uGWZWrXg8WMzeOTojOvPFpIZEGIuaCF/1q3EYhG94dXd932bunBkIlaXfmWNIqloIKR+NTf1Srl3wo0DTJ/q4XOL9iCZWhBXNNcq6FYgWkXlkFR124DeetVmyJKAh16bQjQguQb/WgH2emv5/lM0o2jPo8WEii5rpxuRzfiRaj1eFEhOUVcl7NvUiXTGwMnZePkHNykJVUdYrt/nmcccGsQf3bQbJ2fi+PyPTtTsd6yktZzimFaiI1idmINhUCiaYe+kO4M+3LhnIyhtXZcSkO2plVRrZxzSGb24cUhmW2yzhImkqpnHQ/KqF0A2AMc5I6LVSKpaQdv6WsI6JfBspTpz3e4N+JWrN+MfnjqNyVhtgoIJSzm0ItVSDmkrW8fpevuly83JWK0ajAbq41ZKZwwsWlPa8mFGAMhW8McVDfNxdVXBaAYzMG5zIVqFhKLXtR08Uw68QroB3HDRBlAKO/hZbeJFBoO0AkwRrbczK9s5O6fMvXlXPzZ1B7G9Sg3MvIhtHGqoHBRNB6XAkktri4VExjYOLPsmoZhKYzXBaEY9lFCjWUln6vp59vPeSo2DLXDxKvjNf/2e5/DQkamcYyuKhkiLupXCsgiBZOdjrxW2ODqVgygQ/PA334Q/edfF63puL8Neb62MA6UU6Yy5qLh1TF1MqPZEwrBDOSwkVPRE1qEc1NZVDqbhrJ+a9TWgQpobB4usnF7f7jed0fHI0Rk8fiw3MySezrRkjQNgzp6IVmGmAwvI5s+n7g7LBcdaCZ9oDnZK1sit5MyNz48DUEqx4Ig5RBwxh/nE+txKSZdZ1K3CQnJt92at2AFprhzqD/tQLK9TOSwlzQUyv5FavIVjDoBZJb3ee5d0UQ7tQtAn1kw5MNUAmDve07NxnJwxu96mMjpUzbDnijCXUCyVQSyVdTetBvb3azXl8OkHjuIv/+M1AGb1+FruzVppRPuM1l2tVgkLqq7XrRRLFTEOaa0lC+AYUf/6lUOqiHJoB4Jy7YyD4lAkCwkFf/S9Mzg2tYLvf/Rau7FjvnJg1dq9a3AriQJB0Ce2XED6x8dnoRsUv5/RkVD1Nd2btcJTWRuIM0tjPbCA32QsbVt53aBIqHpLKwezYdv6FjfbOLShcgjJUs3cSs4Ml/mEihPTcaykNXzknkMYnTc3MWwX3BfxQxIInj2zAGB11dFOwn6xYNxoszO7omBmRbF7Ua313qwFmbfPaBySKCDoE9e9+2XKQTcoJpfMnkBMXnt92P16CPlz5wCshZRLtlK7EFinW0k3aNE6Bmdx3cmZOGKpDG7cM4DR+ST+6elRAMgJSB/Y2o0nTsxax9e2AIZkCckWUg4Z3cB8QkUslcGE9bluhFuJK4cGEQlI61cOjhYSzLWUWENH1mYjLK//3rHFMeRr3ftUjJAs2iNS18L3Do/jLZ993NVAO5XDobOLAIBbr9yMjR0BPGolTjgXuut3b7DrIVbbOoMR9kstpRzm4tkU9+NTZrymnm4lUSAQhdr3gHPCjYODaoy7XHYxDsXmR7cSYb+47rx25lYKyO33tlxvQPrMfAJxRXOdwuZUDlNWh9vt/WG8eVefvRN1KoTrL9xgf98dXlv6dVhev5L0EjOO+qdjU8sA6utWAkzXElcODSLqX79xWEpmIBCz3J3NBV5p4Y6sjFAVlUM7xhyC8vqMK+usylweTli2UnfIXOhlUcCm7hDefIE5K0QgyBm9umtDBMNdQeucNbqV8mZRNzszjo65xyZN5dBTR7cSwLpH88Z7DaEabqVYKoPOoA/DXcEC5dCqdQ6AqYqSqg5K1/7mbeeAdNAnrqvx3qKdCFGoHNgQoSFrwd/SG4IoELxpZx8IMQ2Ac6QnIQTv3LsRIz3BNc8wD8tiS8UcZlayRvfo1DIEq4ttPZElkbfPaBRRv68qqaydQR9GekIYY8ahHZSDX4RuNc5bK6mMDlkUIK1xQWpmQutWDqY7s5RyYMaBtSLpCcvYO9TpOuntD2+6EP/xm29e8/WEZKml2mcwt5IkEKyktQKDWg/8ksBjDo0iEpDWna20lMqgMyRjc0/IVg7sOYsNBmkF2GtbT257StXr1h/fawRlcV2N99jUtomlFB4/PoP9f/GwnTlnKwerW6pztvNf3LwHn3x3YWsSWRLQGVr7zjjiF1uqCG5mRUFvWLYbQNY73gCYbiUec2gQEb9kxwfWClMOW3pDWEpmsJBQcWo2AVkUWrrtdLbV89oXuJSqt+zMhnKsNyDtrK956sQclpIZnJwx5ykUKIe+bBPDyzZ312ROecgvtVT7jNmVNPqjfvRbn+FGGAeZK4fGEbViDuvxm8eS5uCUfZu6AAAvji3ipbElXDTUYVc5tiJhqzZhPTEbNgWuHQnJIjSDrmlnSCnFotW2ZSKWwtFJM5uGuTVZLGPPUCcEAvu9WUvCsmgPDWoFZlYUbOgINFg58GylhhENSKAU68rPZsph36ZOiALB4dFFvHo+hn3DnVW8Uu/hHBKzVpKqbrdzaDfY616La2k5rUE3KHwiweRSGketVEvm1mRxoANbu/HC/3gHdm+MVumqi8MUYKuks84sK9gQ9WOgo3HGQZaE5qmQJoScJYS8Qgh5kRByyDrWQwj5ESHkdetrt+PxnyCEnCSEHCeE3Og4foX1PCcJIV8gDZqlGfGvr223YVDEUhl0hXwIyRIu3BjFv78wgYSq45JNrW4cmHJYu2FNZ/S6TtfyEmwxXYtriaWxXjAQRSqjFzR/ZMpBFoW6Zdiw90MrFMIZBsVc3DQOG6KmW6meHVkZzagcrqeU7qeUHrD+/3EAj1JKdwF41Po/CCEXA7gVwB4ANwH4EiGErQRfBnA7gF3Wv5uqcF2rhmUTrbVtd1zVYNBsittlm7tw3ipK2tfixsHeKa7XrdSmyiFoFf6tRTmwNNaLBzvsYwGfkKMcZEmoa3ZNNd4PXmEhqUIzaI5ycMvwqjX+ZlIORbgZwD3W9/cA+EXH8W9TShVK6RkAJwFcRQgZBNBBKT1ITWf/Nxzn1JX1jruMWTs2NlP58s2maAr6ROx0ZIi0ItVoXNjObqWgb+1uGGYc9gxljcNbdvXnxBz8dY53ReyBP82vHFga60BHwFYOjYo5NFNAmgJ4mBBymBByu3VsgFI6CQDWV1aLPwxgzHHuuHVs2Po+/3gBhJDbCSGHCCGHZmdn13nphbAitdUYh4WEapfTs9TBrjzjsGeoo+Vz96sxGrKd3UosEL+WQrgFq8bh4iFTnY70BHHxUAemltNIZ3QoWv2Nrv1+aAHlMG0VwG3o8OPioQ5s6wvjkgbEEJutfcYbKaWXA3gngDsJIW8p8Vg3TUtLHC88SOndlNIDlNID/f3VT7/LupUqf0P//WMn8cF/fBZAdtAPcytt6Q1hS28Ib9zZV+Ur9R7VGA2ZUtvXrbQe48rSWC8YiEASCC7a2IHNPSFQCpxfSkHJGHVXDuEqtcD3AuOWAhvpDmGgI4DHf/+6nFqReuGTBHs4Uz1YV1I5pXTC+jpDCPk+gKsATBNCBimlk5bLiM3LHAcw4jh9E4AJ6/gml+N1Zy1zpCdjKczFFehWMBoAuqyeK4QQPPQ7b1lzC4Jmwi8JEAWyriK4pKq1bSrreuZILyRUiAJBZ9CHO67bgSu2dNtunXMLSaQbqRxawK10biEJvySgP1q/mdFudAYljC2ksOeTD+LP37MXv3z15pr+vjUbB0JIGIBAKV2xvn8HgL8AcB+A2wD8tfX1B9Yp9wH4F0LI5wEMwQw8P0sp1QkhK4SQawA8A+BDAL641utaD9lRoZUHpNnA9uVUBksp83tnRki7+NAJIQjJ4poH/lBKkcq0ccxBXnsq62Iyg+6QD4QQ/N47dgMAZqzuq2MLSaQbqBxaoUr63EISm3tCaFASpc0f3nQhrtzag2NTK7jYEV+qFetRDgMAvm/dMAnAv1BKHySEPAfgXkLIRwCcA/A+AKCUHiGE3AvgNQAagDsppeyTcAeArwMIAnjA+ld31hJUZQNWlqyZuwDQtY62A82MOQ1ubYtBXNGQ0ak9dKbdYDvttaay5ndP7Y/64ZcEnJtPNiTmYNe9tECV9NhCCiM9oUZfBjoCPty8fxg31+n3rdk4UEpPA7jU5fg8gBuKnHMXgLtcjh8CsHet11ItRMHc/a7GrcSMQyyVwVIyA78ktO3udz3N41hLZJYN0m4wt9Ja7t9istA4EEIw3B3ERCyFdMaoe88q9nqaSTkkFA3fODiKyVgK77pkENds7wWlFGMLSVy1rafRl1d32hU6jSkAABaISURBVLORTQlWM/BHN6idRriUVDEfVxtSHOMVIn5pzYsBSxfc0GC/bqNYn1tJxTZHvyRGX8SPubgKRdPr3l5aFAiCvvUPgKonn//RCXz1J2cgEODE9Aq+ffsbsJTMYEXRPKEc6g03DnmYzfcqizksJVWwNkyxVAbzCQW9kfZc3ACz8GmtbqUZR7pgOyKLAgSyOrfSA69M4n89+jrGFlJ22rST/ojfmj1AGtLtNuwXmyZbaSWdwXeeG8PN+4fQE5bxL8+cg6oZdiHhSHewwVdYf1o/jWaVDHQEXHviu+Ec6B5LmR1YG1Ec4xXC/rUHpGctt1J/m7qVzID+6mYgPHJ0BsemVpDK6BjsLFy8eiMy5uOqVQRXf1dnSJaaps7hu4fHEVc0fPiN23Dl1h4omoFXJ2L2NMfNvVw5tD3b+8O478UJUErLZifMO4zDUjKD+bja8pXQpQivx620osAvCeho4YFI5Qj4VjfT4fRcHG/Y3os/+fmLsLWIWymWyoAQNEg5SJ6vkD47l8BfP3AMT70+i8s3d2H/SJc9IvW5MwvQLdfASHf7GQeuHPLY0R/BclrDXFzFZx88hk8/cLToY/OVg+lWal/lYLqVzCKdSip9KaX434+fxJm5BGaWzX75jU4XbCQhWUSqQuNKKcXp2QS294exd7jTzrRz0me5OM1Eiforh7AsrqvupR788KUJPHhkCjfu2YjPvncfADPTa1tfGM+dXcTYQhK9YdnOvmon2u8Vl2GHtfM/NRvHd54bw4qi4c7rd+YMYGcw5SBLAqZiaaQzBnrC7ekzB6zpX4qGO/7pMJ45s4BbLh/Gx995YdEBPjMrCj730HHEUhmzX36bBqMZwVUoh4WEilgqU7JS17lR8TdAOXQGfZiMVeaibRTzCRXRgITP/9f9OccPbOnGw69NY0PU35bBaIArhwLYfN0njs9iPqFC1Qw89OqU62MX4qZx2NITwqlZc+pWO2crhWQJqYyOH5+YRUdAwjcOjuLp0/NFH88awx2dXLaMQ3vGGxjBVaQCn5pNAMi+X93ocyRHBBqgHHrCco669iKzcSXnPjGu2d6LWCqDs/MJvHPvxgZcWePhxiGPoc4gAj4B33/B7AUYkkXc95J7N4+FhIKOgITeiIwzc+aHtZ3dSqyHv2ZQ/NYNuwAAy6nibgUW7Ds6uYKZ5XTbZioxhruCeG1iuaKMpdPWZqRUjKvfseg1Qjn0RvxYSKjrmqxYa+bjCvpcPrM37x/C1z98JQ796c/hN966owFX1ni4cchDEAi290Uwvawg6BNx27Vb8dOTc3aqpZP5hIreiB9dQdmettXe2Uqm+yjil/DWC8zGiCslWpGMLZizLubiCpbTWtu7lW67divmEyruPTRW9rGn5xKQJcGeC+1GXzT7XmyEcugNy1B1Y91z2WuJWZtU+L6TRAHX7d5Q9/oQL8GNgwtMqu/b1Ikb92yEQYHnR5cKHsdSV53tMtzeaO1C2IotvHFnr12xW2phYG4lRru7la7c2o0rtnTj7idPl+3bf3o2jm29YYglBviEZMmuVG6McjDfA8z96kXm4u2dRFIKbhxcYEHpyzZ3Y5NV/DK9XKgcmHFw7i562viNxpTDdbs3IOATIAmkZLX52GIyx2fe3+ZuJUIIbn/LdpxfSuHgqeKxGgB2plI5mHpoVMwBAOYTSt1/dyVouoHFZMY15sDhxsGVHRuYcehCT0iGLAquWRfzCbNdRqelHPySgHCbtpwGgCu2dOOWy4bxrr2DIIRYrUhKu5X2DXdiY4epGNrdrQQAl410AQDOzidcfz69nMZf3X8UowvJyoyDtfA1RDlYKnreo8phwWp94xZz4PBUVlduuHAD/v+3X4C3XtAPQSDY0OHHVCyV8xhKKRbzlENvWG7rPP2esJyTEhgN+Io2MczoBiZjKYz0DGMplcHUcrrt3UpAtptqvsuN8eUnTuEbB8/ibRcO4LY3bC37fGyBbkjMIcKUgzeNw9yKeV3t3PKmFNw4uBD2S/jtt++y/z/YGcBUnltpOaVBM6gZcwiaH4J2dim5EfEXb2I4uZSGQc3KU58o4PnRxbYO5jMIIdjUHcT4Ysr157NxBVv7wvjH2w5U9Hz9zK3UgE7B7O/p1XRW5u7ibiV3uHGogI2dQbwynhuQfn1mBQCwuSdk+9rbORjtRqkOtyyNdVNPEO/ZP4T3XrGpZHC1ndjUHbLvTz6LCRU9ocqNaCPdSgGfiIhfwlzcmzEH5u7iAWl3eMyhAgY7A5iMpXPytV8cM43F/s1dOW4lTpZowFc0W2nMMZc34BNLpmS2GyM9QTvNN5+FhIruVbzPmHFohFsJqF8hHKUUPz4xC8OovKaCGa0+vqlzhRuHChjoCEDRDCwls8HVF8aWMNwVxIZowDYO3C2SS6mA9LGpFYgCwWAnjzPkM9IdQiyVcR1Xu5hcnXLY3m+muzZqd1wv4/DC2BJu+9qz+PGJ2YrPmYur8IkEHUHuQHGDG4cKYAuYM+7w4rkl7LcyS3rCMnwiwSDf/eZQzK304KuTuOfgWdy0ZyMkkb8F82G9fPKD0mYSRGZVyuFNO/vws4+/rWHKrC8iY64O2UqjVnZXsSwvN+bjCnrD7d3ssRT8k1kBGy3jcGYugU8/cBSHzi7g/FLKNg5hv4Tvf/SN+OWrNjfyMj1HNCAhrmg57rjZFQW//e0XsX+kC3/zvoIpsxxk20PnB6UTqg5VN1Y1Z5sQgoGOxqkzUznUPuYwbrnhzhcJ5LvBC+BKw41DBbA8/K/+5Ay+8uPT+LX/+xwAM97A2DvcaY965JhE/D7oBs3pNPrQkSkomoFP33IJv19FYIWX+cph0XLP5M+L9jL16q90fillf11OZ/CBu5/GSStphEEpxecfPo7Xp83j8wmVZyqVgBuHCuiP+iEQ4PDoIjqs3bAoEOwd6mz0pXmaqDW4x+laeujIFLb1hbF7INqoy/I8XSEfIn6pQDkw330zxbZ6wzIyOsVyhXPZ1wozDhNLKbw0toSDp+fx05O5Veavz8TxhcdO4nvPnwdg9VXiyqEoPBJTAT5RQH/Uj+llBb927Vb4fSLGF5N851sGp3EY6DBnbh88NY//7y3buZ+3BKzWoUA5WBW9q4k5NBq7EC6u1LSJHXMnnV9K4cS02bF2Iq9w9RmrffzYQhKGQYu26+aYcONQIRs7g5heVnDzZcN27yVOabLGwcy6eeToDDSD4qY97dkffzWM9IRwdi43uMqMw2qylRoNG361kFCxvb82v8MwKMaXUpAEgrm4ilfPxwCYhZZOnj6zAAAYXUhgIpaCqhnY0oazoSuFu5UqZP+mTrxpZx83DKsgak3PY26ln7w+iw1RP/Zt4u64cmzvC2N0Pgndkbe/kDCNbDMpB9Yvq5YT4eYSClTNwN5h8331pJXOOulQDpRSWzmMzidxcqb8PIx2hyuHCvnzm/d6emiJF2HKIW4Vwp1bSGJHf4S7lCpgR38Eqm5gfDGJLb1mg73FhApRIOgINM/Hlu3MR1eRYrpamEvpqm09eHFsye7lNOFQDqdmE5iLq9i1IYLXZ+I4PLoIINtkk1MIVw6rgC9qqyOrHMwd7/mllJ2JwykN67h6eja7qC4kVXSHmqu5Y0iWMNDhx9l593Yg1YAFo6/a2mMf6wz6ML2ctpXXM2dM1fD+AyMAgEePzqAr5ONdDUrAjQOnZkT82YC0oumYXlYwzI1DRWy33B1sNjlg9VVaRY2DV9jSG66pcmBZXVds6QZrz/WmnX3QDGq3yHjm9AL6o368dbcZ+Hhtchk7uYotCTcOnJrhNA4sOLipmwcAK4FNGDztCEovJFR0NVEwmrG1N1QT5bCSzuBvHjqOp0/PoyMgoTss2zVJbEztxFLKjDecmcfV23qwuSf7/tvJXUolaR7nJafpEAVit+1mu7th3mKkYrb3hXF6No5nTs9jdCGJxaSK7X3Nt6Bt7Qtj9tA44opmbxiqwfdfOI+/f/wkAOCiwQ4AwHB3EKmMbgenJ2NpdIeSmF5WcPX2XgR8IgY6zLR0bhxKw40Dp6aYxiGD80tWi27uVqqY7f0RPHF8Fn/4vZdxbiEJWRRwxZae8id6jK1WQH10PoE9VSwcfeToDLb2hnDrVZvt33HL5ZswvZzGUJepICaWUnbM65pt5r3b0hPG9LLCg9Fl4MaBU1NYf6XxxRQEku1TxSnP9v4wvnt4HHNxQCCAoq2ur5JXyGYsJatmHOKKhqdPzeO2a7fgv791h338A1Z/M0opgj4RE0tpLCXNcb5MKWzuDeHZsws8jbUMPObAqSmsM+v5xRQGO4Pw8S6sFcNcSNv6wvjdn7sAQHP1VWKwVNzVdEzNh1KK3/3Oi/jKj08BMGtmVN3ADRcNuD6eEILBrgAmYyk8fXoeV2/vsYPPb9jeiws3RrmLswxcOXBqSiTgw9yKWaTEP4yrY89QBwQC3HHdDrx73xBG55N2oLWZiPgl9Ef9GJ0zXYvPnV3A2bkE3ndgBDMraYwtpHDFlu6Sz/Gfr0zi3144D4EAb9jRix++PInOoA8HSpw31BnEQ0emYFDgt3Zl79svXbEJv3TFpuq8uBaGGwdOTbl6Ww8+99BxBHwC3rV3sNGX01SM9IRw8BM32C23P9fELc639oZwxsq8+swDx/DyeAzvvGQQd/3nUfzHy5P4zu3X4MBW93hKOqPjrx84ht0DUSwmVbz/KweRzhj48Bu3lpwHctnmLhydXMbHbthl1zdwKodrfE5N+bVrt6I3LCOdMXiNwxpo5CyGanLppi68OLaEsYUknj+3CFU38OjRaTx6dAa6QfFb33rBbkmez/eeH8f4Ygp/9gsX49O3XIKIX8Kfv2cP/uznLy75O3/vHbtx+M9+Drddu5XPJ18D3DhwakrYL+HO63cC4JlK7czP7xuEqhv45H1HYFAzwP7ZB48jrmj43Z+7ALMrCv7m4eNIZ3T86lefwf+xYgsA8O8vnMcFAxG8cWcvbrhoAIf+1FzwBb7g1xTuVuLUnF+5ZjMUzcCNvBtr27J/pAvDXUE8dmwGHQEJb929AT98aQIRv4TfeOt2TC2nce+hMRgUeOr1OTz1+hxUzcAtlw/jubOL+IMbd/Nq5jrjGeVACLmJEHKcEHKSEPLxRl8Pp3r4JRF3XLejKat7OdWBEIJf2GfGnN58Qb/dtv1tF26AXxJx5/U7QUDwrWfP4cY9A7jl8mF8/kcn8BvfPAwAeM+lQw279nbFE8aBECIC+N8A3gngYgAfIISUdihyOJym4t3WAv9zFw3gut39uGJLNz54zRYAZuX8B6/ZgrAs4k9//mJ87r2X4pbLhnFkYhmXb+7CSA9vu1JviBfaUBNC3gDgU5TSG63/fwIAKKWfLnbOgQMH6KFDh+p0hRwOpxqcno1jW1/Y1UVkGBSxVMaeV6EbFP/nx6dwzfbesqmunMohhBymlB4o9zivxByGAYw5/j8O4OoGXQuHw6kR20tUJQsCyRlkJArETmbg1B9PuJUAuEWaCiQNIeR2QsghQsih2dnZOlwWh8PhtCdeMQ7jAJxVKpsATOQ/iFJ6N6X0AKX0QH9/81WKcjgcTrPgFePwHIBdhJBthBAZwK0A7mvwNXE4HE7b4omYA6VUI4T8JoCHAIgAvkYpPdLgy+JwOJy2xRPGAQAopfcDuL/R18HhcDgc77iVOBwOh+MhuHHgcDgcTgHcOHA4HA6nAE9USK8FQsgKgONrOLUTQKyKl1Lt5wOAPgBzVXieZnittXjeat0/wPv30Mv3DvD+622n+7cZwDkAWyil5WsBKKVN+Q/AoTWed3eVr6Oqz7ee19aMr7VG11mV+9cM99DL965JXm/b3D8As6t5fDu6lX7o8eerJs3yWvk99M7zVRuvv952un9Lq3lwM7uVDtEKmkc1I6382uoBv39rh9+79eHl+7faa2tm5XB3oy+ghrTya6sH/P6tHX7v1oeX79+qrq1plQOHw+FwakczKwcOh8Ph1AhuHOoAIWSEEPI4IeQoIeQIIeS3reM9hJAfEUJet752W8d7rcfHCSF/n/dcHyCEvEIIeZkQ8iAh5P+1d3chUpVxHMe/P1opfMvyJZQK6aYyCc0gLSMovLCbAoOKaDe7yQqiuzSCuvHCpSTUC4s0tCIsLLKiooQkKwtM8SXBNIQUSSLzlaLo38XzDA07O7vu7tk5M7O/Dxzm7DPnPDzPn9nzP+fMmeeZVEafGqng+N2fY7dfUncZ/WmkQcRugaSd+TO2U9KdVXXNyeWHJK3SCJjUueD4LZf0i6SzZfVnQIp87MpL3UfIpgI35fVxwEHSdKjdwNJcvhRYkdfHAPOBJcCaqno6gBPApPx3N2kGvdL72CLxm0h6znty/nsDcFfZ/Wuy2M0GpuX1mcCxqrq+B+aR5l/5BFhYdv9aLH5zc31ny+7XhSy+cmiAiDgeET/k9TPAAdLsd/eQDlDk13vzNuciYjvwZ4+qlJcx+axtPL3Me9FuCozfNcDBiKjMFPUFsGiYm1+qQcRuV0RUPlP7gUskXSxpKjA+Ir6NdKTbWNmnnRUVv/zejog43sj2D4WTQ4NJmk46u/gOuKLyYcmvU/raNyL+Bh4H9pKSwgxg3TA2t+kMJX7AIeA6SdMldZD+oa/qZ5+2MYjYLQJ2RcRfpAPi0ar3juayEWOI8Ws5Tg4NJGkssBl4OiJOD2L/UaTkMBuYBuwBlhXayCY21PhFxElS/DYBXwFHgH+KbGOzGmjsJN0ArAAeqxT1stmIedSxgPi1HCeHBskH9s3AWxHxXi7+NV+uk19P9FPNLICIOJwv7d8Bbh2mJjeVguJHRHwYEbdExDzS2Fw/DVebm8VAYyfpSuB9oDMiDufio6Tpeyt6ncq3HRUUv5bj5NAA+fuBdcCBiFhZ9dYWoCuvdwEf9FPVMWCGpMqgWQtI90DbWoHxQ9KU/HoZ8ATwWrGtbS4DjZ2kCcDHwLKI+Lqycb51ckbS3FxnJxcQ71ZXVPxaUtnfiI+EhfTkTJBuA+3Oy92kp2e2ks5etwKXV+1zBPgdOEs6a5uRy5eQEsIe0rgrE8vuX4vF723gx7w8UHbfmi12wHPAuaptdwNT8ns3A/uAw8Aa8o9o23kpOH7d+bP4b359oez+9bX4F9JmZlbDt5XMzKyGk4OZmdVwcjAzsxpODmZmVsPJwczMajg5mA0DSUskdQ5g++mS9g1nm8wGoqPsBpi1G0kdEbG27HaYDYWTg1kv8iBrn5IGWZtNGqq5E7geWAmMBX4DHomI45K+BL4BbgO2SBpHGpr5RUmzgLXAaNIPyB6NiJOS5gDrgfPA9sb1zqx/vq1kVt+1wKsRcSNwGngSWA3cFxGVA/vyqu0nRMQdEfFSj3o2As/kevYCz+fy14GnIo3zZNZUfOVgVt8v8f/4OG8Cz5ImcPk8T4J2EVA9Pv+mnhVIupSUNLblog3Au72UvwEsLL4LZoPj5GBWX8+xZc4A+/s40z83gLrVS/1mTcO3lczqu1pSJRE8COwAJlfKJI3K4/bXFRGngJOSbs9FDwPbIuIP4JSk+bn8oeKbbzZ4vnIwq+8A0CXpFdLom6uBz4BV+bZQB/AyaTrIvnQBayWNBn4GFufyxcB6SedzvWZNw6OymvUiP630UUTMLLkpZqXwbSUzM6vhKwczM6vhKwczM6vh5GBmZjWcHMzMrIaTg5mZ1XByMDOzGk4OZmZW4z/0VAEapTpymwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'][-200:].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "metadata": {}, + "outputs": [], + "source": [ + "first_sept_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n", + " for y in range(1985,\n", + " sorted_data.index[-1].year)]" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": {}, + "outputs": [], + "source": [ + " year = []\n", + "yearly_incidence = []\n", + "for week1, week2 in zip(first_sept_week[:-1],\n", + " first_sept_week[1:]):\n", + " one_year = sorted_data['inc'][week1:week2-1]\n", + " yearly_incidence.append(one_year.sum())\n", + " year.append(week2.year)\n", + "yearly_incidence = pd.Series(data=yearly_incidence, index=year)" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 17, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGDpJREFUeJzt3XGQndV53/HvI60sYVuyEZYckDAiYw21IGNsdgSp29TYDYI4Y5GJ7VFNQG3o0DqkdTKZYlHokLHVGdPJxImGmsIU29gYMCFmTGsrtizHE5pSiVVsF4OsSAlEyMjWepaC8IwEQk//uGfNlbzaPbt79+57734/Mzt777nvee/ZI+3+3vec8743MhNJkmrMm+0GSJJ6h6EhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKnawGw3oNPe9KY35apVq2a7GZLUU3bt2vWTzFw20XZ9FxqrVq1iaGhotpshST0lIv6hZjuHpyRJ1QwNSVI1Q0OSVM3QkCRVMzQkSdUMDUnTcuiFI3zojkc5dPjIbDdFXWBoSJqWLdv38tjTI2z55t7Zboq6oO+u09CJDr1whN+97zvc9uF3sHzxotlujvrIeTdv5eix4z97fs+O/dyzYz8LB+axZ/MVs9gyzSTPNPqcR4GaKY/ccCnvv/AsFi1o/RlZtGAe6y88i0c+dunPbesQVv/wTKNPeRSombZ8ySIWLxzg6LHjLByYx9Fjx1m8cGDMM9r2g5fNv/FLs9BadYqh0aceueFSNn9tN9944kccefk4ixbMY935v8BN73vbbDdNfeQnLx7lqovP4cNr38K9O/czfNKZhAcv/cfQ6FOTOQqUpuqOqwd/9njzlRf83OsevPQfQ6OPTXQUKM00D176j6HRxyY6CpS6wYOX/hKZOdtt6KjBwcH01uiSNDkRsSszByfarmrJbUT8fkQ8ERHfj4j7ImJRRCyNiG0Rsbd8P71t+xsjYl9E7ImIdW3lF0XE4+W1LRERpXxhRHyplO+IiFVtdTaW99gbERsn0wmSpM6aMDQiYgXw74HBzLwAmA9sADYB2zNzNbC9PCci1pTXzwcuBz4dEfPL7m4HrgNWl6/LS/m1wHOZ+VbgU8CtZV9LgVuAi4G1wC3t4SRJ6q7ai/sGgNMiYgB4LfAssB64u7x+N3BlebweuD8zj2bmU8A+YG1EnAksycxHszUm9vmT6ozu60HgveUsZB2wLTNHMvM5YBuvBo0kqcsmDI3M/CHwR8B+4CDwfGZ+A3hzZh4s2xwElpcqK4Bn2nZxoJStKI9PLj+hTmYeA54HzhhnXyeIiOsiYigihoaHhyf6kSRJU1QzPHU6rTOBc4GzgNdFxG+NV2WMshynfKp1Xi3IvDMzBzNzcNmyCT8XXZKmba7eGqVmeOqfA09l5nBmvgx8GfjHwI/LkBPl+6Gy/QHg7Lb6K2kNZx0oj08uP6FOGQJ7AzAyzr6kvjRX/xD1orl6X7ea0NgPXBIRry3zDO8FdgMPA6OrmTYCXymPHwY2lBVR59Ka8N5ZhrAOR8QlZT/XnFRndF8fAL5V5j2+DlwWEaeXM57LSpnUl+bqH6Ject7NW1m16avcs2M/ma1bo6za9FXOu3nrbDetKya8uC8zd0TEg8DfAMeA7wB3Aq8HHoiIa2kFywfL9k9ExAPAk2X76zPzlbK7jwCfA04DtpYvgLuAL0TEPlpnGBvKvkYi4hPAY2W7j2fmyLR+YqmBvEdT75jrt0bx4j6pAQ69cOSUf4i85Ubz3PTQ49y7cz+vmT+Pl145zlVr39Lzd++tvbjP24hIDeA9mnrLXL41iqGhvtW0Ty2cqD1z+Q9Rr5nL93UzNNS3mvbBPxO1Zy7/IVLvcE5DfefkSeVRszWp3LT2TEbTztY0czp6w0Kpl0zms6vnYnsmwyXAOpnDU+o7TZtUblp7argEWKfimYb60uik8kO/8y6uuvgchl88ansmoZfPjjSzPNNQX2rapHLT2jORXjw7UncYGpLG5BJgjcXVU5IkV09JkjrP0JDmKG/DrqkwNKQ5ymswNBVOhEtzjNdgaDo805DmGK/B0HQYGtIc08/XYDhPM/MMDWkO6rUr1Gs5TzPzvE5DUs/r5TsJN4XXaUiaMyYzT+MQ1vQYGmocf6k1WZOZp3EIa3pccqvGqfnEPT8cSCeb6F5ZLjXuDOc01BiTGZe++aHH+eLO/Vy19i2N+ChXNd+hF46w+Wu7+cYTP+LIy8dZtGAe687/BW5639s88KB+TsMzDTXGIzdcespf6lEeLWqq+nmpcTc5p6HGqPml9sI0TUe/LjXuJs80VKVbcwgTjUt7tKjp6LUPw2oiQ0NVaianO6Hml9oPB5JmjxPhGtdkJqdd0ST1Li/uU0dMZg7B9e9S/zM0NO7FdDVzCOfdvJVVm77KPTv2k9la0bRq01c57+at3fwx1MO8oLN3GBqa8AxhohUnrmjSdHmW2juc05jDOnmTt5seepx7d+7nNfPn8dIrx73obpb1yvySNxpsDuc0NKFOniG4/r17aoZyeuXI3bPU3uOS2zmsk9c8uP69e8Zb/txrV8xP5v9gr5w99TtDY47zmofeURMINbdiaZra/4PdulZI43NOQ3NazdFrU45wa2+412/zS857dIdzGlKFmrH/pswP1A7l9Nv8kvMezeLwlOakmqGeJs4P1Azl9Nv8kvcba5aqM42IeGNEPBgRP4iI3RHxyxGxNCK2RcTe8v30tu1vjIh9EbEnIta1lV8UEY+X17ZERJTyhRHxpVK+IyJWtdXZWN5jb0Rs7NyPrrms5ui1iUe4d1w9yOYrL2DNWUvYfOUFJwREP+v22ZMXG55a7fDUnwJ/kZn/CHg7sBvYBGzPzNXA9vKciFgDbADOBy4HPh0R88t+bgeuA1aXr8tL+bXAc5n5VuBTwK1lX0uBW4CLgbXALe3hJE1VzdGrR7jN0e2wbMqQZBNNODwVEUuAXwH+JUBmvgS8FBHrgXeXze4Gvg18DFgP3J+ZR4GnImIfsDYingaWZOajZb+fB64EtpY6f1j29SBwWzkLWQdsy8yRUmcbraC5bxo/swTUDfW4umxuaeKQZNPUzGn8IjAMfDYi3g7sAj4KvDkzDwJk5sGIWF62XwH8n7b6B0rZy+XxyeWjdZ4p+zoWEc8DZ7SXj1FHmpaasf9+mx/Q+HpxyXK31QxPDQDvBG7PzHcAP6UMRZ1CjFGW45RPtc6rbxhxXUQMRcTQ8PDwOE2TpFNzSHJiNaFxADiQmTvK8wdphciPI+JMgPL9UNv2Z7fVXwk8W8pXjlF+Qp2IGADeAIyMs68TZOadmTmYmYPLli2r+JEkaWxNXLLcpIn5CUMjM38EPBMR55Wi9wJPAg8Do6uZNgJfKY8fBjaUFVHn0prw3lmGsg5HxCVlvuKak+qM7usDwLeyddXh14HLIuL0MgF+WSmTpBnRxBVqTZqYr71O498BX4yI1wB/D/wrWoHzQERcC+wHPgiQmU9ExAO0guUYcH1mvlL28xHgc8BptCbARz9w4S7gC2XSfITW6isycyQiPgE8Vrb7+OikuCT1uyZOzHsbEUlqqNpbx3SCtxFRIzVpbFZquiZOzBsa6qomjc1KvaBpE/MOT6krvFOp1GwOT80BvTTU08T7OEmaPEOjh/XSUE8Tx2almdZLB3a1vDV6D2riMrwa3sdJc00/ftqgcxo9qJvL8CRN3mTm8JryyZDOafQxh3qkZpvMHF4vDTODw1M9y6EeqblqDux6dZjZ0OhR3rJbaraJDux69TbshoYkzYCJDux6dZjZ0JCkWdKLw8yunpI6oCkrYKSpcvWU1EW9tgJGmiqHp6Rp6NUVMNJUeaYhTYP31NJcY2hI09CrK2CkqXJ4SpqmXlwBI02Vq6ckSa6ekiR1nqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqpVh0ZEzI+I70TE/yzPl0bEtojYW76f3rbtjRGxLyL2RMS6tvKLIuLx8tqWiIhSvjAivlTKd0TEqrY6G8t77I2IjZ34oSVJUzOZM42PArvbnm8CtmfmamB7eU5ErAE2AOcDlwOfjoj5pc7twHXA6vJ1eSm/FnguM98KfAq4texrKXALcDGwFrilPZwkSd1VFRoRsRJ4H/Df24rXA3eXx3cDV7aV35+ZRzPzKWAfsDYizgSWZOajmZnA50+qM7qvB4H3lrOQdcC2zBzJzOeAbbwaNJKkLqs90/gT4AbgeFvZmzPzIED5vryUrwCeadvuQClbUR6fXH5Cncw8BjwPnDHOviRJs2DC0IiIXwcOZeauyn3GGGU5TvlU67S38bqIGIqIoeHh4cpmSpImq+ZM413A+yPiaeB+4D0RcQ/w4zLkRPl+qGx/ADi7rf5K4NlSvnKM8hPqRMQA8AZgZJx9nSAz78zMwcwcXLZsWcWPJEmaiglDIzNvzMyVmbmK1gT3tzLzt4CHgdHVTBuBr5THDwMbyoqoc2lNeO8sQ1iHI+KSMl9xzUl1Rvf1gfIeCXwduCwiTi8T4JeVMknSLBiYRt1PAg9ExLXAfuCDAJn5REQ8ADwJHAOuz8xXSp2PAJ8DTgO2li+Au4AvRMQ+WmcYG8q+RiLiE8BjZbuPZ+bINNosSZqGaB3Q94/BwcEcGhqa7WZIUk+JiF2ZOTjRdl4RLkmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkaktQHDr1whA/d8SiHDh+Z0fcxNCSpD2zZvpfHnh5hyzf3zuj7DMzo3iVJM+q8m7dy9Njxnz2/Z8d+7tmxn4UD89iz+YqOv59nGpLUwx654VLef+FZLFrQ+nO+aME81l94Fo987NIZeb8JQyMizo6Iv4yI3RHxRER8tJQvjYhtEbG3fD+9rc6NEbEvIvZExLq28osi4vHy2paIiFK+MCK+VMp3RMSqtjoby3vsjYiNnfzhJanXLV+yiMULBzh67DgLB+Zx9NhxFi8cYPniRTPyfjVnGseAP8jMtwGXANdHxBpgE7A9M1cD28tzymsbgPOBy4FPR8T8sq/bgeuA1eXr8lJ+LfBcZr4V+BRwa9nXUuAW4GJgLXBLezhJkuAnLx7lqovP4aHfeRdXXXwOwy8enbH3mnBOIzMPAgfL48MRsRtYAawH3l02uxv4NvCxUn5/Zh4FnoqIfcDaiHgaWJKZjwJExOeBK4Gtpc4fln09CNxWzkLWAdsyc6TU2UYraO6bzg8tSf3kjqsHf/Z485UXzOh7TWpOowwbvQPYAby5BMposCwvm60AnmmrdqCUrSiPTy4/oU5mHgOeB84YZ1+SpFlQHRoR8Xrgz4Hfy8wXxtt0jLIcp3yqddrbdl1EDEXE0PDw8DhNkyRNR1VoRMQCWoHxxcz8cin+cUScWV4/EzhUyg8AZ7dVXwk8W8pXjlF+Qp2IGADeAIyMs68TZOadmTmYmYPLli2r+ZEkSVNQs3oqgLuA3Zn5x20vPQyMrmbaCHylrXxDWRF1Lq0J751lCOtwRFxS9nnNSXVG9/UB4FuZmcDXgcsi4vQyAX5ZKZMkzYKai/veBVwNPB4R3y1l/xH4JPBARFwL7Ac+CJCZT0TEA8CTtFZeXZ+Zr5R6HwE+B5xGawJ8aym/C/hCmTQfobX6iswciYhPAI+V7T4+OikuSeq+aB3Q94/BwcEcGhqa7WZIUk+JiF2ZOTjRdl4RLkmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqoZGpKkaoaGJKmaoSFJqmZoSJKqGRqSpGqGhiSpmqEhSapmaEiSqhkakqRqhoYkqZqhIUmqZmhIkqr1RGhExOURsSci9kXEptlujyTNVY0PjYiYD/xX4ApgDfAvImLNTLzXoReO8KE7HuXQ4SOzvk3NPiSp2xofGsBaYF9m/n1mvgTcD6yfiTfasn0vjz09wpZv7p31bWr2IUndFpk5220YV0R8ALg8M/91eX41cHFm/u5Y2w8ODubQ0NCk3uO8m7dy9NjxnytfODCPPZuv6Oo2NfuQpE6LiF2ZOTjRdr1wphFjlJ2QdBFxXUQMRcTQ8PDwpN/gkRsu5f0XnsWiBa3uWLRgHusvPItHPnZp17ep2YckzZZeCI0DwNltz1cCz7ZvkJl3ZuZgZg4uW7Zs0m+wfMkiFi8c4Oix4ywcmMfRY8dZvHCA5YsXdX2bmn1I0mwZmO0GVHgMWB0R5wI/BDYAH+70m/zkxaNcdfE5fHjtW7h3536Gx5iA7tY2NfuQpNnQ+DkNgIj4NeBPgPnAZzLzP59q26nMaUjSXFc7p9ELZxpk5teAr812OyRpruuFOQ1JUkMYGpKkaoaGJKmaoSFJqmZoSJKq9cSS28mIiGHgH07x8puAn3SxOdPVa+0F29wttnnm9Vp7YXptPiczJ7w6uu9CYzwRMVSzDrkpeq29YJu7xTbPvF5rL3SnzQ5PSZKqGRqSpGpzLTTunO0GTFKvtRdsc7fY5pnXa+2FLrR5Ts1pSJKmZ66daUiSpqGnQyMiPhMRhyLi+21lb4+IRyPi8Yj4HxGxpJQviIi7S/nuiLixrc63I2JPRHy3fC1vSJtfExGfLeXfi4h3t9W5qJTvi4gtETHWh1U1rc1d6eeIODsi/rL8Oz8RER8t5UsjYltE7C3fT2+rc2Ppyz0Rsa6tvCv93OE2N7KfI+KMsv2LEXHbSfua8X7ucHub2se/GhG7Sl/uioj3tO2rM32cmT37BfwK8E7g+21ljwH/rDz+beAT5fGHgfvL49cCTwOryvNvA4MNbPP1wGfL4+XALmBeeb4T+GVan2y4FbiiB9rclX4GzgTeWR4vBv4WWAP8F2BTKd8E3FoerwG+BywEzgX+DpjfzX7ucJub2s+vA/4J8G+B207a14z3c4fb29Q+fgdwVnl8AfDDTvdxT59pZOZfASMnFZ8H/FV5vA34zdHNgddFxABwGvAS8EI32tlukm1eA2wv9Q4B/w8YjIgzgSWZ+Wi2/jd8HriyyW2eqbaNJTMPZubflMeHgd3ACmA9cHfZ7G5e7bP1tA4ojmbmU8A+YG03+7lTbZ6JtnWqzZn508z8X8AJnyrWrX7uVHu7aQpt/k5mjn6y6RPAoohY2Mk+7unQOIXvA+8vjz/Iqx8V+yDwU+AgsB/4o8xs/0P42XKa+Z9mcqjnFE7V5u8B6yNiIFqfXHhReW0FrY/BHXWglHXTZNs8qqv9HBGraB197QDenJkHofXLSOtMCFp990xbtdH+nJV+nmabRzWxn0+l6/08zfaOanof/ybwncw8Sgf7uB9D47eB6yNiF63TuZdK+VrgFeAsWqfzfxARv1heuyozfwn4p+Xr6u42+ZRt/gytf9whWp9c+L+BY7ROL0/W7WVwk20zdLmfI+L1wJ8Dv5eZ451Vnqo/u97PHWgzNLefT7mLMcpmrJ870F5oeB9HxPnArcC/GS0aY7Mp9XHfhUZm/iAzL8vMi4D7aI31QmtO4y8y8+UybPLXlGGTzPxh+X4YuJfun+aP2ebMPJaZv5+ZF2bmeuCNwF5af5RXtu1iJfDsyfttWJu72s8RsYDWL9kXM/PLpfjH5TR9dEjkUCk/wIlnQ6P92dV+7lCbm9zPp9K1fu5QexvdxxGxEngIuCYzR//+dayP+y40RlcxRMQ84Gbgv5WX9gPviZbXAZcAPyjDKG8qdRYAv05r6GXW2xwRry1tJSJ+FTiWmU+W09HDEXFJOS2+BvhKk9vczX4ufXIXsDsz/7jtpYeBjeXxRl7ts4eBDWXs91xgNbCzm/3cqTY3vJ/H1K1+7lR7m9zHEfFG4KvAjZn516Mbd7SPpzJ73pQvWke4B4GXaSXptcBHaa0w+Fvgk7x6AePrgT+jNTn0JPAf8tUVEruA/1te+1PKKpQGtHkVsIfW5Nc3ad2FcnQ/g7T+o/4dcNtonaa2uZv9TGvFS5b3+m75+jXgDFqT9HvL96VtdW4qfbmHtlUl3ernTrW5B/r5aVqLKl4s/5fWdKufO9XeJvcxrQO4n7Zt+11geSf72CvCJUnV+m54SpI0cwwNSVI1Q0OSVM3QkCRVMzQkSdUMDUlSNUNDklTN0JAkVfv/ZquYoDCrbPcAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.plot(style='*')" + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "1986 0\n", + "1987 0\n", + "1988 0\n", + "1989 0\n", + "1990 0\n", + "2020 221186\n", + "2002 516689\n", + "2018 542312\n", + "2017 551041\n", + "1991 553090\n", + "1996 564901\n", + "2019 584066\n", + "2015 604382\n", + "2000 617597\n", + "2001 619041\n", + "2012 624573\n", + "2005 628464\n", + "2006 632833\n", + "2011 642368\n", + "1993 643387\n", + "1995 652478\n", + "1994 661409\n", + "1998 677775\n", + "1997 683434\n", + "2014 685769\n", + "2013 698332\n", + "2007 717352\n", + "2008 749478\n", + "1999 756456\n", + "2003 758363\n", + "2004 777388\n", + "2016 782114\n", + "2010 829911\n", + "1992 832939\n", + "2009 842373\n", + "dtype: int64" + ] + }, + "execution_count": 18, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "yearly_incidence.sort_values()" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 19, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEICAYAAABPgw/pAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAEchJREFUeJzt3X2QXXV9x/H3l0TbwGJEg4sN6lLHUh2jImvxYYq7QVsRH1rHVigqqQ9xxqoZm840dXSoVcdMlY4Og+2kVaGFcavBFgzWwqiLz0gC1gWjQiFiggYpEl2KQOTbP85J3QTYmz3nZs/d375fMzu55+4593z3m7Of/d1zz0NkJpKkhe+wrguQJPWHgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqxNL5XNmKFStyZGSk0bJ33XUXRxxxRH8LKoj96c0ezc7+zK7L/mzbtu32zDy613zzGugjIyNs3bq10bKTk5OMjY31t6CC2J/e7NHs7M/suuxPRPzgYOZzl4skFcJAl6RCGOiSVAgDXZIKYaBLUiF6BnpEfCwibouI62Y896iIuCIibqj/PerQlilJ6uVgRujnAy864LkNwOcz80nA5+tpSVKHegZ6Zn4JuOOAp18OXFA/vgD4gz7XJUmaoziYe4pGxAiwJTOfWk/fmZmPnPH9n2bmg+52iYi1wFqA4eHhEycmJhoVOj09zdDQUKNlFwP705s92t/Urj37TQ8vg913z8+6V61cPj8r6qMut5/x8fFtmTnaa75DfqZoZm4CNgGMjo5m0zOtPIttdvanN3u0vzUbLttvev2qvZwzNT8nj+84c2xe1tNPC2H7aXqUy+6IeCxA/e9t/StJktRE00C/FDirfnwWcEl/ypEkNXUwhy1+Avg6cHxE7IyI1wMbgRdGxA3AC+tpSVKHeu4wy8wzHuJbp/S5FklSC54pKkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSpEq0CPiLdHxPURcV1EfCIifr1fhUmS5qZxoEfESuBtwGhmPhVYApzer8IkSXPTdpfLUmBZRCwFDgdubV+SJKmJxoGembuADwK3AD8C9mTm5f0qTJI0N5GZzRaMOAq4GHgVcCfwKWBzZl54wHxrgbUAw8PDJ05MTDRa3/T0NENDQ42WXQzsT2/2aH9Tu/bsNz28DHbfPT/rXrVy+fysqI+63H7Gx8e3ZeZor/mWtljHC4CbM/MnABHxaeC5wH6BnpmbgE0Ao6OjOTY21mhlk5OTNF12MbA/vdmj/a3ZcNl+0+tX7eWcqTaRcPB2nDk2L+vpp4Ww/bTZh34L8OyIODwiAjgF2N6fsiRJc9VmH/pVwGbgGmCqfq1NfapLkjRHrd5fZebZwNl9qkWS1IJnikpSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklSI+bnflCTNMHLA7e/my46Np3Wy3vniCF2SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIK0SrQI+KREbE5Ir4bEdsj4jn9KkySNDdtb0H3YeBzmfnKiHg4cHgfapIkNdA40CPiEcDJwBqAzLwXuLc/ZUmS5ioys9mCEc8ANgHfAZ4ObAPWZeZdB8y3FlgLMDw8fOLExESj9U1PTzM0NNRo2cXA/vQ2qD2a2rWn6xIAGF4Gu+/uuorB1bY/q1Yub7zs+Pj4tswc7TVfm0AfBb4BPC8zr4qIDwM/y8x3PdQyo6OjuXXr1kbrm5ycZGxsrNGyi4H96W1QezSy4bKuSwBg/aq9nDPVdi9sudr2Z8fG0xovGxEHFehtPhTdCezMzKvq6c3AM1u8niSphcaBnpk/Bn4YEcfXT51CtftFktSBtu+v3gpcVB/hchPwp+1LkiQ10SrQM/NbQM/9OpKkQ88zRSWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSrEgrmB4NSuPazp6N6Lbe4FKEnzxRG6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQrQM9IpZExLURsaUfBUmSmunHCH0dsL0PryNJaqFVoEfEscBpwD/1pxxJUlORmc0XjtgMvB84EviLzHzJg8yzFlgLMDw8fOLExESjdd12xx5239241FZWrVzezYrnYHp6mqGhoa7LGGiD2qOpXXu6LgGA4WV09ju2ELTtT5scGR8f35aZo73ma3yT6Ih4CXBbZm6LiLGHmi8zNwGbAEZHR3Ns7CFnndW5F13COVPd3NN6x5ljnax3LiYnJ2na28ViUHvU1c3PD7R+1d7OfscWgrb9mY8cabPL5XnAyyJiBzABrI6IC/tSlSRpzhoHemb+VWYem5kjwOnAFzLz1X2rTJI0Jx6HLkmF6MsOs8ycBCb78VqSpGYcoUtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEaB3pEPC4ivhgR2yPi+ohY18/CJElzs7TFsnuB9Zl5TUQcCWyLiCsy8zt9qk2SNAeNR+iZ+aPMvKZ+/HNgO7CyX4VJkuYmMrP9i0SMAF8CnpqZPzvge2uBtQDDw8MnTkxMNFrHbXfsYffd7epsatXK5d2seA6mp6cZGhrquoyBNqg9mtq1p+sSABheRme/YwtB2/60yZHx8fFtmTnaa77WgR4RQ8CVwPsy89OzzTs6Oppbt25ttJ5zL7qEc6ba7CFqbsfG0zpZ71xMTk4yNjbWdRkDbVB7NLLhsq5LAGD9qr2d/Y4tBG370yZHIuKgAr3VUS4R8TDgYuCiXmEuSTq02hzlEsBHge2Z+Xf9K0mS1ESbEfrzgNcAqyPiW/XXi/tUlyRpjhrvEMrMrwDRx1okSS14pqgkFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmF8AaCA2wu95pcv2ova/p0b8qFcA9VSQ/kCF2SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVIhWgR4RL4qI70XEjRGxoV9FSZLmrnGgR8QS4DzgVOApwBkR8ZR+FSZJmps2I/TfAW7MzJsy815gAnh5f8qSJM1VZGazBSNeCbwoM99QT78GOCkz33LAfGuBtfXk8cD3Gta6Ari94bKLgf3pzR7Nzv7Mrsv+PCEzj+41U5ubRMeDPPeAvw6ZuQnY1GI91coitmbmaNvXKZX96c0ezc7+zG4h9KfNLpedwONmTB8L3NquHElSU20C/WrgSRFxXEQ8HDgduLQ/ZUmS5qrxLpfM3BsRbwH+E1gCfCwzr+9bZQ/UerdN4exPb/ZodvZndgPfn8YfikqSBotnikpSIQx0SSrEggn0iIiZ/0qS9jfQgR4RSyJiXURsBv4MIN3p3zcRsToijuu6jkFlf3qzR7Ob7/60ObFoPrwA+H3gQ8CG+voxF2WmZ7O1UF9zZwK4E7g/Is7OzCs7Lmtg2J/e7NHsuurPQI/QgdcCl2fm5cA7gccCL+u2pIUnIo6NiEfMeOpVwMWZeTLVRndGRDynm+q6N5f+LNZdfm5DsxuU/gx6oH8dGKkfXwvcAJxQj9TVQ0Q8OSI+C3wF+JuI2HfxtF8Ah9ePP0l1fYpnL7awatKfxbbLz21odoPWn0EP9JuAoYh4dGbeXU9DdblePYiIOGLG5DOAnZk5AnwB+GD9/B3APRFxZGbeAXwfGOZXfzyLNUt/vsgD+zNU9+cGFkl/wG2ol0Huz6AH+hRwD3BKPf1Tqqb8pLOKBlBEHBUR50fE1cDGiDi6Hgk8DfhqPbK8FLgzIk6j2riOBFbVL3ED1ZXk7u2i/kNtlv6s4lf9uYQH9udp9Ut8n4L7A25DvSyU/gx6oO8CvgG8rZ7+CVWg/7yzigbTycBe4MVUV8F8B/AIqv/fY2bsJrgA+BPgm1Q9PBUgM78OrAZ+Nr9lz5uH6s8S7M8+bkOzWxj9ycyB/wIuAi4DbgPeCiztuqaO+rAEeBNwJdU15lfUz38SeFv9+DhgY/39Z1Ht21tSf2+I6o/iEmAlsB14C/Bx4CPA4V3/jPbHHtmf5l+DPkLfZw3VX8QnZOa5mbm343q68hKqo3zeDTwH+Nv6+SuA59aPfwh8GTg1M6+mGk2MA2TmNHAV8KzM3AW8huodz4+Bd2bm/87Tz3Go2J/e7NHsFnR/Bv04dAAy8z7gv7quYz7sO5IiIp5F9dbty8BlmXkP8FvATZn5hYi4GfhARPwesA34w4hYkZm3R8QNwF0R8XjgXODVEfEYqmvW/w/V20Eycyuwdd5/yBbsz8GzRw+u5G1ooYzQF4UZG9rJwMeoDn16AfD+epb7ge9HxLLMvJnq7dzTqPbV3Up17CvAL6ne8h0GXEx12c8zgROBTZl5/zz9SH1V/9wZEWPYnwcVEY+u/30ucD72aD8R8fB6G3o+1W6QovqzIEboJYuIw6lOoDoZ+GJEnA88E/iHzDwvIo4CroyIE6juEnUC1Vu4HVRv7V5GdfTPvwPvj4gLgScAy4Fb63c3X4uIly7gX8IjqU7O2A28jqoH9of9BgF/BKwD9kbEm6hu4v4RewQR8TDgjVQfUO6MiHOoDjf8+9L64wi9QxFxDLAFGAP+hepDlldQ7bvbC5CZPwUuoTrS5wrgMcCT65f4EvB84N7M3AJ8FNgMnEf1aft9+9bV9YbW0jLg14AnRsQK4IlUo6NF3586zJcDfwx8KDPHMvN7VCPFffMs6h5RXQdqNdUlRJLqd2w11Si7rP4cyk9c/er5ifoy4KQZ02uoPoQ5C/jmjOd/g2okANXGOQEcVS//GeDxM+Zd0fXPdQj6dBbwAeBdwOuBlwJX25///5neDLynfrzvpjWvAK6yR0n98722fvx6qiPlXlXi75gj9G79AvjmjNOBv00V8BdQjUaPAcjMW4HrI+KkzDwPuBH4V6q3h1/JzFv2vWAWdOGyGX05DPhv4HpgdWZ+BjguIoZh8fZnhtuB342IM4FtEXEB1ejzt+t3NIu9R/8BnBURnwL+Gng6cB3VPZGPhnL64z70DuW+IdWvgmsd8On68WeBtwN/GRGPAn4A/Kj+3ruB44EbM/MX81fx/NrXH6qTOU6nOsZ3XUS8k+owsLXAexZrf2a4lmqX1ElUZ1WfSr2bAHhdRHwQeCSLtEeZ+ZGIuJfqsMMTqI4LfyXVeS1vjIiNFNIfR+gDIDMzIo4FjqEaTUA1kiAitgCfpzpx4ZZ6/vsy87qFtKE1FRFDVCPQfwQ+R3XEwUnAGcBREfEZFnF/ardQhffSrPYHb6mf20J1gahLsUdPASbr0fXHqT6D+QTVLpVi+uMIfXCcQHXtmpsj4g1Ub/XeQbWv77uZeU2XxXVoL1Uo3Ud1CeVfAu/IzCngzyPi1cD1mXlthzV2KjPvqUeZ762fOozqg/WNmXntYu9RVFdn3Ul1JNk/A3dRneH55sy8uaT+7PsARR2LiK8Cv0l1qNStwLsz89udFjWA6hM5XgFMZOaPu65nkETE+6i2oROAy4Gz6xH7ohcRT6Q6Vvxeqh79G/DerM7sLIaBPgDq42TPBm4GLszqjDXNUI+y7k832FlFxPHADxbaroL5UA8GngR8LavLcRfHQJekQvihqCQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5Jhfg/y3YIomGrCoMAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.hist(xrot=20)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], "metadata": { "kernelspec": { "display_name": "Python 3", @@ -16,10 +2262,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.3" + "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 } - -- 2.18.1