From 53fddfa693e9f49e1dc8cead3679d177cc58d57d Mon Sep 17 00:00:00 2001 From: 1cda4be30895c9886fda013d21530393 <1cda4be30895c9886fda013d21530393@app-learninglab.inria.fr> Date: Thu, 22 Oct 2020 01:16:58 +0000 Subject: [PATCH] exercisemod32 --- module3/exo2/exercice.ipynb | 2294 ++++++++++++++++++++++++++++++++++- 1 file changed, 2291 insertions(+), 3 deletions(-) diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb index 0bbbe37..62a9f43 100644 --- a/module3/exo2/exercice.ipynb +++ b/module3/exo2/exercice.ipynb @@ -1,5 +1,2294 @@ { - "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=\"https://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
02020427343513505520528FRFrance
12020417401021305890639FRFrance
2202040720786753481315FRFrance
3202039710492371861213FRFrance
4202038722537823724315FRFrance
5202037715844052763204FRFrance
620203679191001738102FRFrance
7202035782801694102FRFrance
8202034722723714173306FRFrance
9202033712841772391204FRFrance
10202032726506894611417FRFrance
11202031713031002506204FRFrance
1220203071385752695204FRFrance
132020297841101672102FRFrance
14202028772801515102FRFrance
1520202779861491823102FRFrance
16202026769401454102FRFrance
1720202572280597001FRFrance
1820202473880959102FRFrance
19202023755811115102FRFrance
2020202272770633001FRFrance
212020217602361168102FRFrance
222020207824201628102FRFrance
2320201973100753001FRFrance
242020187849981600102FRFrance
2520201772720658001FRFrance
262020167758781438102FRFrance
27202015719186753161315FRFrance
282020147387922275531639FRFrance
29202013773265236941611814FRFrance
.................................
15291991267176081130423912312042FRFrance
15301991257161691070021638281838FRFrance
15311991247161711007122271281739FRFrance
1532199123711947767116223211329FRFrance
1533199122715452995320951271737FRFrance
1534199121714903897520831261636FRFrance
15351991207190531274225364342345FRFrance
15361991197167391124622232291939FRFrance
15371991187213851388228888382551FRFrance
1538199117713462887718047241632FRFrance
15391991167148571006819646261834FRFrance
1540199115713975978118169251832FRFrance
1541199114712265768416846221430FRFrance
154219911379567604113093171123FRFrance
1543199112710864733114397191325FRFrance
15441991117155741118419964271935FRFrance
15451991107166431137221914292038FRFrance
1546199109713741878018702241533FRFrance
1547199108713289881317765231531FRFrance
1548199107712337807716597221529FRFrance
1549199106710877701314741191226FRFrance
1550199105710442654414340181125FRFrance
15511991047791345631126314820FRFrance
15521991037153871048420290271836FRFrance
15531991027162771104621508292038FRFrance
15541991017155651027120859271836FRFrance
15551990527193751329525455342345FRFrance
15561990517190801380724353342543FRFrance
1557199050711079666015498201228FRFrance
15581990497114302610205FRFrance
\n", + "

1559 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202042 7 3435 1350 5520 5 2 \n", + "1 202041 7 4010 2130 5890 6 3 \n", + "2 202040 7 2078 675 3481 3 1 \n", + "3 202039 7 1049 237 1861 2 1 \n", + "4 202038 7 2253 782 3724 3 1 \n", + "5 202037 7 1584 405 2763 2 0 \n", + "6 202036 7 919 100 1738 1 0 \n", + "7 202035 7 828 0 1694 1 0 \n", + "8 202034 7 2272 371 4173 3 0 \n", + "9 202033 7 1284 177 2391 2 0 \n", + "10 202032 7 2650 689 4611 4 1 \n", + "11 202031 7 1303 100 2506 2 0 \n", + "12 202030 7 1385 75 2695 2 0 \n", + "13 202029 7 841 10 1672 1 0 \n", + "14 202028 7 728 0 1515 1 0 \n", + "15 202027 7 986 149 1823 1 0 \n", + "16 202026 7 694 0 1454 1 0 \n", + "17 202025 7 228 0 597 0 0 \n", + "18 202024 7 388 0 959 1 0 \n", + "19 202023 7 558 1 1115 1 0 \n", + "20 202022 7 277 0 633 0 0 \n", + "21 202021 7 602 36 1168 1 0 \n", + "22 202020 7 824 20 1628 1 0 \n", + "23 202019 7 310 0 753 0 0 \n", + "24 202018 7 849 98 1600 1 0 \n", + "25 202017 7 272 0 658 0 0 \n", + "26 202016 7 758 78 1438 1 0 \n", + "27 202015 7 1918 675 3161 3 1 \n", + "28 202014 7 3879 2227 5531 6 3 \n", + "29 202013 7 7326 5236 9416 11 8 \n", + "... ... ... ... ... ... ... ... \n", + "1529 199126 7 17608 11304 23912 31 20 \n", + "1530 199125 7 16169 10700 21638 28 18 \n", + "1531 199124 7 16171 10071 22271 28 17 \n", + "1532 199123 7 11947 7671 16223 21 13 \n", + "1533 199122 7 15452 9953 20951 27 17 \n", + "1534 199121 7 14903 8975 20831 26 16 \n", + "1535 199120 7 19053 12742 25364 34 23 \n", + "1536 199119 7 16739 11246 22232 29 19 \n", + "1537 199118 7 21385 13882 28888 38 25 \n", + "1538 199117 7 13462 8877 18047 24 16 \n", + "1539 199116 7 14857 10068 19646 26 18 \n", + "1540 199115 7 13975 9781 18169 25 18 \n", + "1541 199114 7 12265 7684 16846 22 14 \n", + "1542 199113 7 9567 6041 13093 17 11 \n", + "1543 199112 7 10864 7331 14397 19 13 \n", + "1544 199111 7 15574 11184 19964 27 19 \n", + "1545 199110 7 16643 11372 21914 29 20 \n", + "1546 199109 7 13741 8780 18702 24 15 \n", + "1547 199108 7 13289 8813 17765 23 15 \n", + "1548 199107 7 12337 8077 16597 22 15 \n", + "1549 199106 7 10877 7013 14741 19 12 \n", + "1550 199105 7 10442 6544 14340 18 11 \n", + "1551 199104 7 7913 4563 11263 14 8 \n", + "1552 199103 7 15387 10484 20290 27 18 \n", + "1553 199102 7 16277 11046 21508 29 20 \n", + "1554 199101 7 15565 10271 20859 27 18 \n", + "1555 199052 7 19375 13295 25455 34 23 \n", + "1556 199051 7 19080 13807 24353 34 25 \n", + "1557 199050 7 11079 6660 15498 20 12 \n", + "1558 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 8 FR France \n", + "1 9 FR France \n", + "2 5 FR France \n", + "3 3 FR France \n", + "4 5 FR France \n", + "5 4 FR France \n", + "6 2 FR France \n", + "7 2 FR France \n", + "8 6 FR France \n", + "9 4 FR France \n", + "10 7 FR France \n", + "11 4 FR France \n", + "12 4 FR France \n", + "13 2 FR France \n", + "14 2 FR France \n", + "15 2 FR France \n", + "16 2 FR France \n", + "17 1 FR France \n", + "18 2 FR France \n", + "19 2 FR France \n", + "20 1 FR France \n", + "21 2 FR France \n", + "22 2 FR France \n", + "23 1 FR France \n", + "24 2 FR France \n", + "25 1 FR France \n", + "26 2 FR France \n", + "27 5 FR France \n", + "28 9 FR France \n", + "29 14 FR France \n", + "... ... ... ... \n", + "1529 42 FR France \n", + "1530 38 FR France \n", + "1531 39 FR France \n", + "1532 29 FR France \n", + "1533 37 FR France \n", + "1534 36 FR France \n", + "1535 45 FR France \n", + "1536 39 FR France \n", + "1537 51 FR France \n", + "1538 32 FR France \n", + "1539 34 FR France \n", + "1540 32 FR France \n", + "1541 30 FR France \n", + "1542 23 FR France \n", + "1543 25 FR France \n", + "1544 35 FR France \n", + "1545 38 FR France \n", + "1546 33 FR France \n", + "1547 31 FR France \n", + "1548 29 FR France \n", + "1549 26 FR France \n", + "1550 25 FR France \n", + "1551 20 FR France \n", + "1552 36 FR France \n", + "1553 38 FR France \n", + "1554 36 FR France \n", + "1555 45 FR France \n", + "1556 43 FR France \n", + "1557 28 FR France \n", + "1558 5 FR France \n", + "\n", + "[1559 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": 4, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
\n", + "
" + ], + "text/plain": [ + "Empty DataFrame\n", + "Columns: [week, indicator, inc, inc_low, inc_up, inc100, inc100_low, inc100_up, geo_insee, geo_name]\n", + "Index: []" + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "raw_data[raw_data.isnull().any(axis=1)]" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \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
02020427343513505520528FRFrance
12020417401021305890639FRFrance
2202040720786753481315FRFrance
3202039710492371861213FRFrance
4202038722537823724315FRFrance
5202037715844052763204FRFrance
620203679191001738102FRFrance
7202035782801694102FRFrance
8202034722723714173306FRFrance
9202033712841772391204FRFrance
10202032726506894611417FRFrance
11202031713031002506204FRFrance
1220203071385752695204FRFrance
132020297841101672102FRFrance
14202028772801515102FRFrance
1520202779861491823102FRFrance
16202026769401454102FRFrance
1720202572280597001FRFrance
1820202473880959102FRFrance
19202023755811115102FRFrance
2020202272770633001FRFrance
212020217602361168102FRFrance
222020207824201628102FRFrance
2320201973100753001FRFrance
242020187849981600102FRFrance
2520201772720658001FRFrance
262020167758781438102FRFrance
27202015719186753161315FRFrance
282020147387922275531639FRFrance
29202013773265236941611814FRFrance
.................................
15291991267176081130423912312042FRFrance
15301991257161691070021638281838FRFrance
15311991247161711007122271281739FRFrance
1532199123711947767116223211329FRFrance
1533199122715452995320951271737FRFrance
1534199121714903897520831261636FRFrance
15351991207190531274225364342345FRFrance
15361991197167391124622232291939FRFrance
15371991187213851388228888382551FRFrance
1538199117713462887718047241632FRFrance
15391991167148571006819646261834FRFrance
1540199115713975978118169251832FRFrance
1541199114712265768416846221430FRFrance
154219911379567604113093171123FRFrance
1543199112710864733114397191325FRFrance
15441991117155741118419964271935FRFrance
15451991107166431137221914292038FRFrance
1546199109713741878018702241533FRFrance
1547199108713289881317765231531FRFrance
1548199107712337807716597221529FRFrance
1549199106710877701314741191226FRFrance
1550199105710442654414340181125FRFrance
15511991047791345631126314820FRFrance
15521991037153871048420290271836FRFrance
15531991027162771104621508292038FRFrance
15541991017155651027120859271836FRFrance
15551990527193751329525455342345FRFrance
15561990517190801380724353342543FRFrance
1557199050711079666015498201228FRFrance
15581990497114302610205FRFrance
\n", + "

1559 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202042 7 3435 1350 5520 5 2 \n", + "1 202041 7 4010 2130 5890 6 3 \n", + "2 202040 7 2078 675 3481 3 1 \n", + "3 202039 7 1049 237 1861 2 1 \n", + "4 202038 7 2253 782 3724 3 1 \n", + "5 202037 7 1584 405 2763 2 0 \n", + "6 202036 7 919 100 1738 1 0 \n", + "7 202035 7 828 0 1694 1 0 \n", + "8 202034 7 2272 371 4173 3 0 \n", + "9 202033 7 1284 177 2391 2 0 \n", + "10 202032 7 2650 689 4611 4 1 \n", + "11 202031 7 1303 100 2506 2 0 \n", + "12 202030 7 1385 75 2695 2 0 \n", + "13 202029 7 841 10 1672 1 0 \n", + "14 202028 7 728 0 1515 1 0 \n", + "15 202027 7 986 149 1823 1 0 \n", + "16 202026 7 694 0 1454 1 0 \n", + "17 202025 7 228 0 597 0 0 \n", + "18 202024 7 388 0 959 1 0 \n", + "19 202023 7 558 1 1115 1 0 \n", + "20 202022 7 277 0 633 0 0 \n", + "21 202021 7 602 36 1168 1 0 \n", + "22 202020 7 824 20 1628 1 0 \n", + "23 202019 7 310 0 753 0 0 \n", + "24 202018 7 849 98 1600 1 0 \n", + "25 202017 7 272 0 658 0 0 \n", + "26 202016 7 758 78 1438 1 0 \n", + "27 202015 7 1918 675 3161 3 1 \n", + "28 202014 7 3879 2227 5531 6 3 \n", + "29 202013 7 7326 5236 9416 11 8 \n", + "... ... ... ... ... ... ... ... \n", + "1529 199126 7 17608 11304 23912 31 20 \n", + "1530 199125 7 16169 10700 21638 28 18 \n", + "1531 199124 7 16171 10071 22271 28 17 \n", + "1532 199123 7 11947 7671 16223 21 13 \n", + "1533 199122 7 15452 9953 20951 27 17 \n", + "1534 199121 7 14903 8975 20831 26 16 \n", + "1535 199120 7 19053 12742 25364 34 23 \n", + "1536 199119 7 16739 11246 22232 29 19 \n", + "1537 199118 7 21385 13882 28888 38 25 \n", + "1538 199117 7 13462 8877 18047 24 16 \n", + "1539 199116 7 14857 10068 19646 26 18 \n", + "1540 199115 7 13975 9781 18169 25 18 \n", + "1541 199114 7 12265 7684 16846 22 14 \n", + "1542 199113 7 9567 6041 13093 17 11 \n", + "1543 199112 7 10864 7331 14397 19 13 \n", + "1544 199111 7 15574 11184 19964 27 19 \n", + "1545 199110 7 16643 11372 21914 29 20 \n", + "1546 199109 7 13741 8780 18702 24 15 \n", + "1547 199108 7 13289 8813 17765 23 15 \n", + "1548 199107 7 12337 8077 16597 22 15 \n", + "1549 199106 7 10877 7013 14741 19 12 \n", + "1550 199105 7 10442 6544 14340 18 11 \n", + "1551 199104 7 7913 4563 11263 14 8 \n", + "1552 199103 7 15387 10484 20290 27 18 \n", + "1553 199102 7 16277 11046 21508 29 20 \n", + "1554 199101 7 15565 10271 20859 27 18 \n", + "1555 199052 7 19375 13295 25455 34 23 \n", + "1556 199051 7 19080 13807 24353 34 25 \n", + "1557 199050 7 11079 6660 15498 20 12 \n", + "1558 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 8 FR France \n", + "1 9 FR France \n", + "2 5 FR France \n", + "3 3 FR France \n", + "4 5 FR France \n", + "5 4 FR France \n", + "6 2 FR France \n", + "7 2 FR France \n", + "8 6 FR France \n", + "9 4 FR France \n", + "10 7 FR France \n", + "11 4 FR France \n", + "12 4 FR France \n", + "13 2 FR France \n", + "14 2 FR France \n", + "15 2 FR France \n", + "16 2 FR France \n", + "17 1 FR France \n", + "18 2 FR France \n", + "19 2 FR France \n", + "20 1 FR France \n", + "21 2 FR France \n", + "22 2 FR France \n", + "23 1 FR France \n", + "24 2 FR France \n", + "25 1 FR France \n", + "26 2 FR France \n", + "27 5 FR France \n", + "28 9 FR France \n", + "29 14 FR France \n", + "... ... ... ... \n", + "1529 42 FR France \n", + "1530 38 FR France \n", + "1531 39 FR France \n", + "1532 29 FR France \n", + "1533 37 FR France \n", + "1534 36 FR France \n", + "1535 45 FR France \n", + "1536 39 FR France \n", + "1537 51 FR France \n", + "1538 32 FR France \n", + "1539 34 FR France \n", + "1540 32 FR France \n", + "1541 30 FR France \n", + "1542 23 FR France \n", + "1543 25 FR France \n", + "1544 35 FR France \n", + "1545 38 FR France \n", + "1546 33 FR France \n", + "1547 31 FR France \n", + "1548 29 FR France \n", + "1549 26 FR France \n", + "1550 25 FR France \n", + "1551 20 FR France \n", + "1552 36 FR France \n", + "1553 38 FR France \n", + "1554 36 FR France \n", + "1555 45 FR France \n", + "1556 43 FR France \n", + "1557 28 FR France \n", + "1558 5 FR France \n", + "\n", + "[1559 rows x 10 columns]" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "data = raw_data.dropna().copy()\n", + "data" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "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": 7, + "metadata": {}, + "outputs": [], + "source": [ + "sorted_data = data.set_index('period').sort_index()" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "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": 9, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm8HUWZ9vOec7ckZCcJIQHCEpaEJZAQwAWBgERQQYQx4wIi86EOoDNuA66oMOKgoI6Cg4LgBqKg7MoWZE8IawgkJEAgIfu+3Nybe8+p74/u6q6ufmvps9x7k9NPfjfnnOrqquru6nrr3UkIgRw5cuTIkUNFobcHkCNHjhw5+h5y4pAjR44cOVLIiUOOHDly5EghJw45cuTIkSOFnDjkyJEjR44UcuKQI0eOHDlSyIlDjhw5cuRIIScOOXLkyJEjhZw45MiRI0eOFJp6ewCVYtdddxXjxo3r7WHkyJEjxw6FZ599do0QYoSr3g5LHMaNG4c5c+b09jBy5MiRY4cCEb3lUy8XK+XIkSNHjhRy4pAjR44cOVLIiUOOHDly5EghJw45cuTIkSOFnDjkyJEjR44UcuKQI0eOHDlSyIlDjhw5cuRIIScOOXZ6/GPeCqza3NHbw8iRY4dCThxy7NTo7C7hs797Fp/41azeHkqOHDsUcuKQY6eGEMHn2+vae3cgOXLsYMiJQ46GgOjtAeTIsYMhJw45cuTIkSMFJ3EgojYimk1ELxLRPCL6blh+KRG9Q0QvhH+nKOdcQkSLiGgBEZ2slE8mornhsZ8REYXlrUT0p7B8FhGNq/2l5siRI0cOX/hwDp0AThBCHAZgEoDpRHR0eOxqIcSk8O9eACCiCQBmAJgIYDqAa4ioGNa/FsD5AMaHf9PD8vMArBdC7AfgagA/rP7ScuSIdQ45cuTIBidxEAG2hD+bwz/bK3cagFuEEJ1CiDcBLAIwlYhGAxgkhHhKCCEA/BbA6co5N4Xf/wJgmuQqcuSoBiLXNuTIURG8dA5EVCSiFwCsAvCAEELaBV5IRC8R0Q1ENDQsGwNgiXL60rBsTPhdL0+cI4ToBrARwHBmHOcT0RwimrN69WqvC8zR2Mg5hxw5KoMXcRBClIQQkwCMRcAFHIxARLQvAlHTcgA/DqtzO35hKbedo4/jOiHEFCHElBEjnImMcuSoim94dfkmzF+xqWZjyZFjR0ImayUhxAYAjwCYLoRYGRKNMoBfAZgaVlsKYA/ltLEAloXlY5nyxDlE1ARgMIB1ma4kRw4GQrIOFVCJD/z0MUz/yWO1HVCOHDsIfKyVRhDRkPB7PwAnApgf6hAkPgLg5fD7nQBmhBZIeyNQPM8WQiwHsJmIjg71CWcDuEM555zw+5kAHhYiFwjkqB75JMqRozL45JAeDeCm0OKoAOBWIcTdRPQ7IpqE4P1bDOCzACCEmEdEtwJ4BUA3gAuEEKWwrc8DuBFAPwD3hX8AcD2A3xHRIgQcw4waXFuOHLnOIUeOCuEkDkKIlwAczpR/ynLO5QAuZ8rnADiYKe8AcJZrLDlyZEZOHHLkqAi5h3SOnRq6KetFNz+P3z/9Vi+NJkeOHQc5ccixU0MXK9314jJ8828v85Vz5MgRIScOOXZq5FKlHDkqQ04ccuzUyI3ecuSoDDlxyLFTQ0SffZtIPPTqSlz9wGu9PYwcOSLkxCHHTo0dhXE476Y5+OlDC3t7GDlyRMiJQ46dGn2dY8iRo68iJw45dm7ktCFHjoqQE4ccOzVy2pAjR2XIiUOOnRo7is4hR46+hpw45NipkesccuSoDDlx2IGweM1WLFnX3tvD2KGQcw45clQGn6isOfoIjvvRIwCAxVec2rsD2YEQ+TnkRCJHjkzIOYccOzVyD+kcOSpDThxy7NTIaUOOHJUhJw45cuRoOJTLAjfPfhurNnf09lD6LHLikGOnRs455OAwf8VmXHL7XHzzr3n4dhNy4pBjp0ZuypqDQ3e5DAB4Z8O2Xh5J34WTOBBRGxHNJqIXiWgeEX03LB9GRA8Q0cLwc6hyziVEtIiIFhDRyUr5ZCKaGx77GRFRWN5KRH8Ky2cR0bjaX2qORkTOOdhxxX3zMe7ie3p7GD0OOS/y+WGGD+fQCeAEIcRhACYBmE5ERwO4GMBDQojxAB4Kf4OIJgCYAWAigOkAriGiYtjWtQDOBzA+/Jselp8HYL0QYj8AVwP4YQ2uraYolwV++9RidHSVensoOTIgf/ft+OU/X+/tIfQKhPaZIw0ncRABtoQ/m8M/AeA0ADeF5TcBOD38fhqAW4QQnUKINwEsAjCViEYDGCSEeEoE9oW/1c6Rbf0FwDTJVfQV3D13Ob59xzxc/WAec39HgjRlzReBHCqieZGzDkZ46RyIqEhELwBYBeABIcQsAKOEEMsBIPwcGVYfA2CJcvrSsGxM+F0vT5wjhOgGsBHA8EouqF7Y2tkNANjY3tXLI8mRBfmrnyNHZfAiDkKIkhBiEoCxCLiAgy3VuR2/sJTbzkk2THQ+Ec0hojmrV692DTtHjlymnINF7jnvRiZrJSHEBgCPINAVrAxFRQg/V4XVlgLYQzltLIBlYflYpjxxDhE1ARgMYB3T/3VCiClCiCkjRozIMvQcOwA2d3RhQ/v2Grfa2G///BWbMHfpRme9RhOvRArpBp8fNvhYK40goiHh934ATgQwH8CdAM4Jq50D4I7w+50AZoQWSHsjUDzPDkVPm4no6FCfcLZ2jmzrTAAPiz46W99YsxVX3De/4V6mnsAR338Ak773gHf9je1dOO0XT+DtteZghI3+mKb/5DF86OePO+s13n2SOodeHkYfhg/nMBrATCJ6CcAzCHQOdwO4AsBJRLQQwEnhbwgh5gG4FcArAP4O4AIhhDTx+TyAXyNQUr8O4L6w/HoAw4loEYAvIbR86ouY/eY6/PKfr2Pp+trYR2/u6MJrKzfXpK0dHV2lbG/qvS8vx4tLNuCaRxYZ6zTiuy+EwOML12TawDTafcqJghvOqKxCiJcAHM6UrwUwzXDO5QAuZ8rnAEjpK4QQHQDO8hjvTodP/+YZPPvW+jzSap0gF4FSuXFWg1ueWYJLbp+Ln3xskvc5ASHpUwaCdUVuyupG7iFdIco12no8+9b6mrTTyLA9ikaUKcucH0vX++f+aLS7FDvBNdqV+yMnDhWiq1SuaXv5JK0PGvG2FkIXoSzMUiPeJ6DxiGIW5MShQnR01ZY4NJDUo0fRiIuedB+t5NpnLlhl1eHsLBCxuVIOA/JMcBVie405h1JZoFhoHJlvLeBztxpRrCSDC7iu/a21W6Pvsu65v3kGAPDvx+1Xp9H1DeQ6BzdyzqFCdNacc8inaT3QiLdVEk0bN7p6cyfed+Uj0e9Gu0+5zsGNnDhUiFot5gVPEUAe8C+HL6KoZJZJVXtnwx0LklPKSYMZOXGoELXacFCkPLQ3+MiCPFyICTbxSSNuDHOFtAca7XorQE4cKkStZNmxCCCfrfVAI+ocJDdqm1N6zONGvE9AAxLFDMiJQ4Wo1aSqZJeXwx+N+PLHCmkz9PvSaPdJXm6+KTMjJw4VolZTSu7gyjl1qAvUu1ov5eO4i+9JmX/eP29FFOa9p0EM5+C69kabfXkmODdy4lAharXQFDx1DjnMIItRq/qcfnx//RI1/c/fF0TfF63agvN/9yy+dttLdelr2/YSrvzHfHR2p40UXn5nYzSWTdti4qRPr5RYqcHmX6OK0bIgJw4VouacQz5XK4btRVeJ7s9n9oxzl+QYZBiLWuMXMxfhFzNfxx9nvZ069s/XYsOFm2fHx0va4p8SK9V2iH0eDUYLK0JOHCpFjXUOjbZz88Wzb63HX59f6q5oQJZIr0vXt2P+ik0AgOUbq4+6W69Hui00a+7OcG0uzrTRpl+c7KfBLjwDcuJQIWpvrVST5nY6fPTaJ/Gff3qx4vOzRGN9zw9nYvpPHgMAvP/qRyvus5rs50IIXHrnPCy0hHGX6xnXT8HQeVnz2dSrtW/vbqiFMs8t7kZOHCpE7fwcgk+d7d+Z8feXl+P+eSuqbsdnEa40QOLmDj9lcq0X1KXrt+HGJxfj02EYC7ZPy5JmisCicw76sI/5wcO45ZklaDQ00GuXGTlxqBA1M2UN3+ZGslb63O+fw/m/e7ZH+soieqkE3Dyo94Ij2+e4BBPB9Nl8XHL73GqGtUMhjq0kcPdLy/DsW6msxA2PPPBehaiZQlq21zi0oUfRXWeiW+vWbSKjuI4w1jFZbonahgLb8aGYsl74x+cBIE+4pSHnHCpEbsrad2C7dd26sN2rveqeRTU6h1dDhbiVOMh+MvTdSGJLH+SmrG44iQMR7UFEM4noVSKaR0RfDMsvJaJ3iOiF8O8U5ZxLiGgRES0gopOV8slENDc89jMKXTmJqJWI/hSWzyKicbW/1Nqi9qas+WStByoRK3V2+xOUWuscPushbpNBGAuMgoFMCuleml9vr23H+q19L8hfns7BDR/OoRvAl4UQBwE4GsAFRDQhPHa1EGJS+HcvAITHZgCYCGA6gGuIqBjWvxbA+QDGh3/Tw/LzAKwXQuwH4GoAP6z+0uqL2r1refiMeqIShfS27f4RcK0hKqpYemyOfbfOWRrW4c7j0Vs6rWOvnIn3XTmzV/q2IfeQdsNJHIQQy4UQz4XfNwN4FcAYyymnAbhFCNEphHgTwCIAU4loNIBBQoinRLDd+i2A05Vzbgq//wXANDJtgfoMahuyO+ccssO2gEpkMWWV2FZleHSfcTnb8GmCqWS2VqpuPNVgk8Xya/GarXjflTOxanNHD45IfXvz986ETDqHUNxzOIBZYdGFRPQSEd1AREPDsjEAVJu4pWHZmPC7Xp44RwjRDWAjgOFZxtbTqH3gvXyS1gM+CulSWeC3Ty2OfmchDvV6bF60gSnjRE1A39U5XP/4m3hrbTvum1u9aXMl6KO3pU/AmzgQ0S4AbgPwH0KITQhERPsCmARgOYAfy6rM6cJSbjtHH8P5RDSHiOasXt27+Q1qNadkatB6m1zWGt2lMq55ZFGfT0LkoxP485wl+PYd86Lf2cRK9XluPowza8pqqHvpnfPQXePUtrWAFPs1FXtWUJA7wbnhRRyIqBkBYfiDEOJ2ABBCrBRClIQQZQC/AjA1rL4UwB7K6WMBLAvLxzLliXOIqAnAYAApw2MhxHVCiClCiCkjRozwu8I6oVY7jubwpajUWau3cOucpfifvy/AL3ooXpENdrm/G7rDWy0c9KpFhVIlI1F54JWVeGzhmuoGVQfI8CYvLdmIJ1/vufHl4TPc8LFWIgDXA3hVCHGVUj5aqfYRAC+H3+8EMCO0QNobgeJ5thBiOYDNRHR02ObZAO5Qzjkn/H4mgIdFgzy15mLwCLZnsJDpafxx1tu46ObnEy/Sxm1dAPr2uIHKiPjPHvYneL05S7OYsgJ903xTmhr/ac4SfPxXsxy1a4Nv/HUuvntnwCn2vTvSd+DjBPduAJ8CMJeIXgjLvg7gX4loEoL7uxjAZwFACDGPiG4F8AoCS6cLhBCST/88gBsB9ANwX/gHBMTnd0S0CAHHMKO6y8qOLZ3dmP3mWpxw4Civ+rV60SLi4OQcem8af/2vgefs1f9yWMT+S6LQ0tT7rjK2XfYOu8fwYB14sVIft+PQ0Bsc8x+YaLY50nASByHE4+Cn6r2Wcy4HcDlTPgfAwUx5B4CzXGOpJ759x8u4/bl3cP9/Hov9Rw101q/VmiMXV9cOvC+YupaEiCaMfKklcetNVCtWqpddXDVzxGtIbOA9W/W+RziyRM0FAuum4bu0YGBbc03631H3Dj2B3n+z+whWbeoEAKzY6GdSV6s51eIpVqrEJLPWUF8kyem0KpxDR1cJd7zwTp/arfdUnCMVtSA2Pgpp1lrJcl6fFCtl5ByO+9EjmHHd0zXrvy/N1b6GPLZSiIFtwa3Y1NHlVb9WkyriHBwvSV8wdVXH0BlaKalipR/9YwF+/fibGNK/Be/bv3cNBiTqfdfqZq3kU4cPruSF3p9NASox0pu3bFPN+u8Dr1WfRc45hJCLnGpS+r27XsG4i++pa7/SWskVsqEvTGKVe9ke3idVrLR8U8B1bdrmR2BN8Ca8HgthvXeGdfNz8NI5cGXmE/uiWKnvjSiHRE4cQsh1T323bnjiTWP9ntc59D51SEq2gh+2xUiHEALXP26+p3G9jAOrM/qq6CGLn0Ol9eqN3o6D0DefbN9AThxCyMXXN2pHteKEUlngrheXoclT59BbKoev/DnOwlZtfJ7n3t6A79/9irNeLQlhLZqytVGvx+Kzy+fuk3X6KsdqMe67X1qGdX0wqF4W9FXC3xeQ6xwkogQq6UMb27tSopJq59RvnngTl93zavS7r+oc/vJsHPHENIbt3WWs3drpbKuz28/zOOuV2hfv6u9bWQgUTHkS6vRcfPYonJGClZOr4VBXberAhX98HkeOG4o/f+5dFbfT2xxMThrMyDmHEHLh416uw753P35w3/xEWbVrwspNSasoF+fQF3Y4anwedThf+fOLOOYHD6MrvIZqR9rnOIde7NuGQf3S5pw22lDL+yo3M8s2VBcwL0t8zTfXbHXW2dTRhWk/fgTzlm30ajNLePZGQ04cQkRiJc/61bxmTyxag189lpS991WxkgqT2eYDr6wEkN2haWsnH62zVmvY1s5u3Pdy9aEwekOs5IOs96ke5tCVbFpKZZHaHLnw9tp2HP+jR5z1Zr2xDq+v3oqr7n/Neyw5eOTEIUTWOV7NTv4LNz+fKnOJlfrCJDbtPOXmzxZFMShPHlHFaipqZaw08Tv/wAtLNrjbUXavE3cflB6PhQTUz1qJv7rVm1XxXbbO+8IcAoAf3b8AR/33Q1i1qcP6DF9csgHjLr4HLy7Z4G1iLv1uco6geuTEIYTQPn3rV4J2JurnDiFWYhYXIeJFOuvaYzJ57U3LrCKjdOqN4ZgWzT/Meiv6nnVctcinvXZLZ2KuVpJ2Zeb8VQCA1Vs6rWKwh8J6MxesQv+WormiAkkcHl+0Bt/829zMY8sRIycOIaIQvj1AHbh8Aa4FsS9s+tQhqi+1XCBqlW2Ma6WniGMW01wAdZMr+QyD69q2WKuLeqW3c/JlD+KLt6Q5Xw4mTqWeebxUp8zfP53HUKoGOXEIIeexr3VLrT1j3cSh96mD8WUPP+UYTSP1XRO4a+0p4piVc6hfPgd3naxTokOxFqtm3Koex0a0b/DwafHV8vmOti9sonYW5MQhRDkj5/Bft9WWZXX12xcmvbpoJ8Ybvt9ZxRamBUowErae4hyKzKrcGzoHH2zrKhmV+mz9DEmMOGR9Bj6pP61EMDHf/PqutV6lo6uEH9z7Ktq3+9/nnQU5cQgRJRzvgb64F8I1p/uCzsHkdFXwFCv5W4Ixug3Pc6tFgXkjKr311YzZ5AS3SlFIf+XPL2Lid/6hnWdGZ5ViJfXxvrW2PejPFujPOaf9+iWQd91aE4ffP/0W/u/RN/DLR16vWZvbtpewdH17zdqrFxqeOMx+cx0O+OZ92NAeeHr2xCLM7k53ALGS6b2Tl1OrPMVcP7brr1a0oz4NVqxkObdeT4Vbczu7S/hjFbkI1JSulRGH+KRP/HpW2E5cNnfpRrz8jtu/gAzf7X371as1cZBWhNtrmMb3Mzc+g/f8cGbN2qsXGp44/GLmInR2l/Hi0mBS98QazCWBd/XbF7KIGk1Zw0/5Ym70DLxnumaOUPYUbeQU0jbCXTcPaaas2oUvoZCugKy5Nigf+vnj+OD/Pq704UY1eigOtdqg6HhkwSrcOmdJTdp66o21NWmn3mh44qBPzp6Iec9xDn1RIf13zYHMZX0ij3/rby+z9XzRm/qVpoycA4esj2rZhm3pQmaO/OTBhRlHkkS1c4h3gswuVqrEWMl36NVYzHHnSvHe/BWb8bW/vFRx2zsicuKg/e6JNZgTXdRK59BdKuPqB17DuIvvwU8di8kbq7ckRA06Pvf7ZxO/TUHW5NUsXLXFa4wSWTgHXt9RG5NItRlWrGTh2qwiJ89n9tW/vJgq00exalMHrnv0Da/2zOPhv/uCewY9pSj2JWzV+HLUi+vYUZETB22B6S3i4OJYfOf83S8tx08fCojCdY+alWjt27txwo//iS/fml6YTHjq9Zgdfn11TAjkLXQ58vmC93OoSdNOsGKljNZKWbjP595ej3Vb3WK4S++a592mCdVyZFnPV+/Dj+5fwNapdY4JnQhlCenS0x7kfcHIxAYncSCiPYhoJhG9SkTziOiLYfkwInqAiBaGn0OVcy4hokVEtICITlbKJxPR3PDYzyhcmYmolYj+FJbPIqJxtb9Uw/Vpv3vLWsltyupn1qeG4bDtrDu7gnpPvL7G3rGCax55HYvD4GfPLF4f9KH874LvTp/3c+iZF6mpmE0f9CITnkPW97neM655Eq8uT2c2009t4syoGNQz8F7W89Xqjy8M5tmG9u1RQq0szfn2rdfz87Uw91HPfBN9wTzdBp8Z1w3gy0KIgwAcDeACIpoA4GIADwkhxgN4KPyN8NgMABMBTAdwDRFJ3/drAZwPYHz4Nz0sPw/AeiHEfgCuBvDDGlybF/SH31tWQW5rJbWuuZ6qz6jHvH7u7fWpskpfINMOm3tprHenho+stSkdpsHW/L/9do7xWDU7Q/2W+oaP8O3SbWbKiJCYTbiNAKrvkmSWJ33vASxYuVk53z4OifnLN7srIS1WuvPFZYnftmfCcQ71DCneFywQbXASByHEciHEc+H3zQBeBTAGwGkAbgqr3QTg9PD7aQBuEUJ0CiHeBLAIwFQiGg1gkBDiKRE8od9q58i2/gJgGvluM2uNnrBWYhXS9nPUSW2rmhBZ1eEOVpG22Hsh6E1rpdam9CuRdZGXteev8FvQOOjTX8YZqgYJ7tMx0blLzqpzSNQ3iOt800987TY/ZXBXFdGNuWP15Rx2cOKgIhT3HA5gFoBRQojlQEBAAIwMq40BoNp8LQ3LxoTf9fLEOUKIbgAbAQxn+j+fiOYQ0ZzVq1dnGboFms6hB6gDN9+yWCvZXkhOn1Fv1PoFYuX4PfQisQtgxjZqMVb9liajsVr6thzrabGSuthWM0X+9sI73nUXrU4aRegLvk2vUKvYYL4o9wHzdBu8iQMR7QLgNgD/IYRIC0mVqkyZsJTbzkkWCHGdEGKKEGLKiBEjXEP2QsqUtZeIuatbVa9mm8MqcagHmeAUiLVWKvLEgRtLbaC2w1lcZZ0TtZhCFYvqrD4Z/He2LlOWWSHNiJV02OaOPN0nyY+E1KVJ6Au+jcBx1kq1ntu+Y+kL8CIORNSMgDD8QQhxe1i8MhQVIfyUfO9SAHsop48FsCwsH8uUJ84hoiYAgwGsy3oxlaBeCumfP7wQ4y6+B92MtQSvkHax+X4iAVVk1ROSOSL/hczXbLgWu/dKwes7+vZL7IuE3spQ50sn7W88P7N4LSFVMlKHmkKfO/qCb1uQe5xz2NGJQyj7vx7Aq0KIq5RDdwI4J/x+DoA7lPIZoQXS3ggUz7ND0dNmIjo6bPNs7RzZ1pkAHhZ1lCOc9csn8b+huWe9OIdfzAzMSPkkPozOwcFiGoPeaUhwDnXYfXLY0uEXlExv1dRLb1orsdeeecdc/TjqsWNNbDAybEaisoz9cQrpZB/AY6/VSjxsGIO24Otipe+dNjE+1sPWSotWbcFmzyRGvQEfzuHdAD4F4AQieiH8OwXAFQBOIqKFAE4Kf0MIMQ/ArQBeAfB3ABcIIaSn1ecB/BqBkvp1APeF5dcDGE5EiwB8CaHlUz0ghMAzi9fjxw/waQRrtQjJ3WaXZ0wW1+60kvAZpnm9YMVm/OXZQP3TrY3vT8+8jb0vuRerDGkcOWK6OUNkUB9wd8L2WG5/3l8m7cJpk8akyrLPiPiMt9ZuxbiL78E/5mVMV1oxYTcf85nbslterFSNzoG/oE2eGwv/Pl2cQ7K+upny8XOoJXfxkWuexCfDGFV9EU2uCkKIx2GeqtMM51wO4HKmfA6Ag5nyDgBnucZSC7gmo+vR3/SZqTjnhtnWOkvXt6MjlH1yTjjcLso157rLfhE1S2W3n8PJP3k0+r5FW9j/EAZ2e4cL6VAlvKOysrvW6l7KC47fN+LmbDhmn5QdRFWcwN+eXxZ+voOTJ+7mfZ56r2rmXOghVrL7SWTrT12oWVFqtua8oD8rffOjL+5Zvca7ymW0FvzMin0gY7r1RTSch7Q+OVI7GscM2XfEALZcJQJqxEV9cgKV6RxU+23bQqnSoko2n5JY7NLq3DcAAL571yvebfsuBr4KaRXjLr4Hv37MHF7iqycfaDymEtFmzgku4zKmjvXqBwMONYunro4seRhsY31S8XB3+zmky7Lumh2WrHVBinPQxUracWGpy8F03xav2YojL3+Qj5PlwJotfpZoPY2GIw76s00H3rOjX3MR7x2/Kw7fc0hUtmjVFoz/xn24S3O4AfwXBdfLWlKIjJVzcOzWXJDJY5qL/NTQuREu5akJ+rjNCmn3uRwuu+fVxG8uiJ4LHLelj+fSO+fhgj88Z2yDG6qveDEeR/y9vctf9GK7T+9s2Ib5K2yGhvH1szk1Mm711RzhvG9P7XkHvcWUWCnFOQhjXQ4mAnLz7LexenNnyunOB5UQlJ5A4xEHh/WCa34UiFAgSiwYMoTCQ6+uTNXnFNKc/NX1oiQ5BzOSkz/74ijFYfVh+f1avfulZbj9uaXJ/AMVjGhwv+bM5/hwdTc+uRj3zF1ubIO7zO6MRu3qHGmvMoObig3tcsHOfj+zLuaqaJIofR85h7Vq7VBS73fKlFWvb64LIBWY0hi2Ppw4b631N7uV8PV+72n4yQ52IqR2Ftr8dE3OQoFQ0Ca6zLDV1px+yJxYyWdcOlRdgu0lvewefzGPdTyml6CaNrXfjyzgvX7/9+FFAICXFHlsJXrAgW1NWGuIJGsC63BTA0rZVRK46cnFmPvORvzorMPc41AGkkXnULvwGemyrMRB5Zq5DRHHdQpRnQhKp8E6N6Bz1mUH57B4bTJjm4nGS0b75tlLcNnph2RyRu2rFq0NyDkkf+sT3rUIFQgh5xBXlLsLjjhwYiXfNKGfufEZ/GJmsFB2qZyDYYybOrqwZkuwB2rmAAAgAElEQVS8GFbzkvWIstBxs9UcxByxcl0fJ8pwoRa+IdxYu0plfOfOeZGVWBZkWZS99TqG8loopE/7+eP41aNv4PXV8S464ByS9bhw8T7XyumFJHQOc/+RA5Pta++Rek2cyMjXiU4lflm5xGrCjNcTjUcctMmjTwjXYyoWCESU2EFssxAHXqzEjIuZdA/PX4Ur/xGEOlZ1DqZBfuQXTzj78UU9djPX/jNbHl51oa5kPLW6BNm3EAJX/mN+Rf1ym4RKvZnfvV/aosoXrjlh91j2u6MvLt2Iy+9N6n+IKHVfPvf7tN7GpwebsYS+zqY3f2adA7emu3QYEiqjkDX0d0+HCvdFwxEH/WmniINLrESBWEmdZPLFb2F2NGx2KWZ79tjCNdZJ4mOtpO7Ugn7SdfxluibZqvvMFiZ4HQA8mtHhqVAtcfA8Sb0mW9yrv7+8wssclgMnXqzUJ+FfpuxhPKZe8zPfONFYzyWq4sVKTD1rK3qb7to+wQH7t5iJg95FSqxki63EmVA7dBgS6judlRPIOYc+grTOIduDkQppbmfnY+liw9OW3LIlTz8HFVy9WCFpR1XztUZzXaW1lSikKxkGR/zky/t5zULJtNhxxSznYBmHekxvj4scy51HBJxyCO9boWf5U88J2mFELFWwk9u7/ZTqPl1kyemtL+a6ct+lc0hzIny/CWe6jJZppT4aga/xiINjZ+GanMUCoVDgo6RyCwsbI96w++5nsVrwtVZSIcKx/fTBhVi+MbAc8bWpNuaL9hBW1cpEsWqxkuUcs8I9fX2me/HIAp4T4hZWbndoW+Q2tG+PFja9PZOZMYdrPjGZLd/MOIP+6KzDjE93c0cXfvvUW6lyX9Hltq6y17x1zZ225oI3UQXS7/cmLVyF+lg4rqASsVJXVp1DRmLSU2g84qA97pTds2MKFyhYtHwjOGYxzRtgYZe7S2li5IIQAvNXbMbVD76GC//4PJ58fQ2+d7efNZOROHisBrUiDuqCWmub+CzNmXxVjPm3mba3djKWOZY+X1u5BT8J43/pY7VtItRGs+qczpwcx8XU+/z2HfNw8+y3bd1ZsW17t9c9f2e92+Y/izhOf3Sq74XeFjfnKxEr2aQRJ171z1RZrnPoI3BxDq7nRERoLlBysQ4/fz/rLdyhxZ6vRQKRDe3bExPIdyoVCxSdt217CR//1Sw8ttAvLWg1E9b3TJcfQiEhVkrDdR+tuZ/VdhxtmmTCb61rZ8s5bNyWNql1LZb/eDmIx6QveANbzfdNveZKLK9isVIS6zKaBOvo6Cp7iQa/+hd7Up9Bbc3WjULKGlF7dh1aSG91TNyGTy8x9a3qxzhOYHt3Ge3bu7GICQmf6xz6CFIPW3sw14Smoza0NhXRqchQ5XxZvbkTX7zlhUTdLDHiTS/PpO89kGBVfXfRRUU3QmSOqc/BZI5H4PMmq6jVJt+lkK7EVr8ScGHXAeCK+3jLJa5bzkNaf956uOxCQXorJ9G/tX5OU+a5GeObpx6k1PfDhw/b3aueTdcyuF8zph00MlNCo3S4DLMYmZ3y+mbSsJCr79byjemglWdc+wQmfPsf7LlZTV97Co1HHJTZMOWyB1OBr3yiRLY1FxI7ENuOKEvScttilrAu8Vz0ipr1lMkEkJvw27v5Tp5fsgGnaSazKuRL4iP62rjNP1xxJZ6z+inqfTfqHDJwDr79AsD+o3Zx1vvCtPGJ3/G9TNZrsegc1LrVmTKbr/mwPYZgyl5Dg3qOdga2BXOutangR6wtlU6btHtgElsF52A7zucR8ZMsqBuZf/m/p7BE4ypffscctiRraJWeQgMSh/i7Szk7YfQgtrytmeccOLCmrB5j07EtEUrCD0VtpTOJGc698RllDEHrpt3MCmZXpKIpXLh819PHFvqZt2Z9fZaub7dGljW1x92jWigMOyuIrBpbwAhDeRrVjtQkVlLRVCCMHtIPgJuLjXxEPPu31esuizA6ge38+ODug9tShF2ee/05UzDzK8dh+C6t0TFWrKQV/e6pxWy/+rR535UzzYNM9ZEThx0OY4b2Y8tbmwLOwc9uO11WSZ7nDsUE78nX/fQGhQJFLwsX20aC8z8wLYiuoctgd74TfrElBeRflTwNWd+fF5dUHwpZeuJmtT7hOElOee26JrkblfVGDQoWMt2PJAsXVCnUPpqLhWiD49bRyfN9zVTNx0olAQJZCZL6qNpaigmuuDv0UgeA/UcNxN67DsCMI/fAue8eF57rJg43MRZbQPqdzsJs9lHa0HjEIcuDMNVtDT2h5W4wiwwUMJsi2sRTy5Qdu8wRoOL11WlFV1NB0Tkg7aFqA5/Bzh2SQhKHWuvYsojnTFCrJ5pTU6sqxYeMGQwgu9067+fgFlnokOuNvJdXnHEo7v3Ce7GrstvV+0uKldI3yFdJbXtPmoqUWPRtiAhc+M/Zr6VOxDlkOF8lav/78KJU/pJCgfDRIwIrLTZ8Rg+s3H1UH92AxCGbTydbWowWQJH45JAlBr7vPOScoKb9OG0it9dwPveED0ycg2txkWKlWudd9r03HJFk2/PwAP/Oh4IUklkVhlzLbAIjb84hqNjSVMCE3QNR5yeO2jOqJ+ffpo4u3PeyOVqsD8ggV1Kfe1Mh5hxcHKIvEZFwmakGoWtsOj61seRlvLEmGetJQt7nanKXu0SPNsfFvppLuvGIg+M5qMo+s/LJry1TGz7nPfvWOuMxLoYTh1QYiwxz8II/mvMV2BCLlfzqf+uOeV71fInN/fPSYdPZ9gzNEQhjQnm69CeohalhJaEndGsllSxP3H1wqu2v3Ppi0lTZQMc3WrzkffiK5iJFxMK1sKlak6xiJb26fA5258b4+9hh/VlnVSDJAcvNHscs+85jPZaUDltukb5JGjyIAxHdQESriOhlpexSInpHyyktj11CRIuIaAERnayUTyaiueGxn1E4u4iolYj+FJbPIqJxtb3EJGwPYsTAVgwb0BLXdXjRRhPPJif1sICI+4u/P/gqH84a8I//LkTc09x3NnrlenZNVNciLUVmtd4M9eTm6rbPvws3fHpKpk2ACn3enHjQSNYyy7XrLmo6B3XlVp+DnIe+qV3f88OHnXVs4pmWpiw6h/gmVvsIS+Uymgq8A+qiVVvQVQr0gBN3H4S5l74fY4b0M45P5RzkfpD3fk6WjWX0kG94cKs2jntHVkjfCGA6U361EGJS+HcvABDRBAAzAEwMz7mGiORKdi2A8wGMD/9km+cBWC+E2A/A1QB+WOG1eMH2IJoUBe4vZi7CTEN4BJ3zfnmZWfn5kuYTsL27jNdW8pNJfX1s82XKuGHmg2p7NZxzpx4yGoA7fEZT0W9HmRW+zclno76LT18yLSxz74uJgN0Gt+GEA0d5745TY9V+mwLFuVotRpxDUNPk9yG/67tT0+XaNgk+YqD+LU0RofJZ2Igk5+CjczCjuyTQVCykRDjLN27DiVf9E5fd/QoEgvs0sK05ZYSRUDMp85g08V1iPFrR0vXbUkYUn/z1LOs1Bf2ZscOKlYQQjwIwyziSOA3ALUKITiHEmwAWAZhKRKMBDBJCPCWCJ/BbAKcr59wUfv8LgGnk8xZXCNtORw2o9+P7g1DZA5hdejSZQjb0iUXmgHm/fvzNxO95FkKSZKnTA1W5Gh9UMulMN97XwqpocNyqFtXsO3cb3JZuzyhWiqFbC/li3ZakN7Hp3nHt3njukXH/moguoVBPfA9+NWmGDpW8RD7n9G8pRourz70hS73Job+EhI2AlMoCzUVCVzlpKbg2vN+zF68P9RJBeRA92TAmlXMgKVby0zlccvvcxO/VVeaA7qM+cFXpHC4kopdCsZN8wmMALFHqLA3LxoTf9fLEOUKIbgAbAVQesN4J8+QrFmKLHkkAjt1/BJ7/1kmJepG4ASJz0nUpN73ijEPsI2OajRThnn2WBfD82+szjc8EX+vb5oL0c6gtefDV3VS9ICYUlbJvv2uRC9aX//xiotxIWJlmjztgJIb2b070HxGHpBdf9FXem0pMpE2wXXFzsRAtrj46hyCXA0/eP3n0nkwpj+6yQHMxcKZjA1oiuBfyPhV0hznla1KsZCEOzPUN0JxJTU5siXMtj6Zv8g2VE4drAewLYBKA5QB+HJZzt0BYym3npEBE5xPRHCKas3p1ttwAUcOWJ1FMmH4GGNyvORUDSJW3ZnVukgv7nsP6R96j8dh4FjgaX0YxhxDAd++qTdrQYsFvqkScA3NbjtnHTvNt/GIluRlsMForKdMxtmLxa9NUz6SMNI1B9hvpHBRfFb7f4LieIa2y2Ep+52Txc5CcA0/M/cdYKotIbJmIUhwRz2CeyNut53o39Vso8O+VEIIVLQ9q88uurI5xpxQrcRBCrBRClIQQZQC/AjA1PLQUgJqJZCyAZWH5WKY8cQ4RNQEYDIMYSwhxnRBiihBiyogRIyoZupVKB3MkWaO5WIgmT1RPcfRq3+5W8qooRxPZbrDN7WJiE1q/vioKOWEo940SLRcobuFziYZsLxCrKsy42CTEMj7iEG13PJoRTyXHwzeaRawExNeqK/eNYqWy7Kd6sVI8Nsez8uaqKNI5cA9Rp0W25rrL5YgzVSPlynk1b9kmvL2uPbpuPUd0Mihh3G4sVkr2d9dLvFlwm8EgZLdByfmRyKG9kyqkUwh1CBIfASAtme4EMCO0QNobgeJ5thBiOYDNRHR0qE84G8AdyjnnhN/PBPCwqOPd8uEcukvliOqrDmtycVB3TXryEHf/QbucM4/6m3vpIvbX8/bUckdSIPOir8IWPsMlW7W1XOspkVBOEv9d9zNwDSE752AYm5A6hKT+xpTfQj5nUz+zvz7N0FManJvDk4vWpKL5RjoHbcxsmwjeK27u6IumbX5JnQMA/Pj+1+JzlFPeWtsePTfZLwdVuV8wWCutN0Si1cPSmMbuGzOpj9IGOPkjIroZwHEAdiWipQC+A+A4IpqEYG4sBvBZABBCzCOiWwG8AqAbwAVCCLl6fh6B5VM/APeFfwBwPYDfEdEiBBzDjFpcmAm2yVcsBIlEvvm3yGoXzU3BRPjdeVNxwKggWTkpC6Up1r8JcgEpFNIBxDgLFBVZQ1PU0vPSV99hG6PNmWzjti5v+/WozDIOb1ml5TxdrOQijEYxkZFzCOpfePx+eP/EUanjLcWk/kZdk9R5ZyIOsnyksqN1SY1iy5247OOMNU5s1eRBOMkmQkv+tnMOItp83PjkYlz64cBJUd8EqQppszm6OgZ+zqpje+xrx+O9/zOTHbNp7EnOgT8H6Lse0k7iIIT4V6b4ekv9ywFczpTPAXAwU94B4CzXOGoFO+cQTJB75sbsZGs4Gd87PhZjqeZ+WR+snMj8BIsbs4mVfHMt1HLOyb5dDmHNFs7BtpH6yYOvmQ/CJFbidqLmNhI5MUyLhrqj1EQnLppsOq4v2kKIUEkbYOSgVhw6dkh8XJ4ndQiMWEmN1SQvq0nTObQ2pcUfLuIuW/AXKyXHzNWT0lquSVc4FhV7Dx+Q0qtwfUuuplCw6ByUZuSc1fNqq3Nhj2H92XLbOBLEgR8GgJ1M57Ajw/YcZPwh9UFycZBUE0fXS3T0PkmfhFK0CyR0aJORExWo6Amdgwm2EAMqdNt8iVJZWPNAuAheVqswF7y4iBTnYIfp3uhGC7Iap0sA4mtt1kKRqIsSFzJeX7RSHvIAbntuaapMRSw+dCG527bNi8jPwdhKjIVMMhwgsO677CMHo4kxjNC7bm0uRP0mPaTVfuOed2ltAlE6XL+JcJnK9XH4RvPtm6ShEYmD5VEE1g7Jl6yZecFinYNwLtT6hJEvU5EotSCqv/4wK52S0cT+mmB7YY87wKDQN5wScQ6GCT889MGIFNJatccX2SPJ3jpnifX4nLeSJrnbu8t12XElxQ3Bp6/OwXT8bk2xGcnpI5ZA5yyCT5tCWuUc5HGfXfi6rfYcGt4mqhrnYKpPkLJ/YeD0zGNW34/jDhhpdCbUJ63kmHSdQ4I4KK91sUAY2NqEje1JHYNpaGaL4eQ4tlehkF66vh0PveoXCqZeaDzi4HzBkxWsnAPcL5FOAKTYnXuRXWOTYgNvsZKl2ujBfDhyE/RggwAwpH86XWXEOWh9n3PDbGv7evpGHVf+Y0Hi9/7fvA//ddtcQ203TPeGDchWjpdzG3xSSCbqGTiHdEjudDtqGAe7qFIfi99x11zUmzHVP+OIsdGmK+t4VDGPzZJab1vGxCpoug5V56V3O7CtGVu0PN+msZl0SPpraRMrXXnmofF5zPv8gZ88hvNumsMPoIfQcMTBBkLwvqrvcgsj45RPulwWTuLQpRMHRrn4oTCFoosjyG53b67IxYixQWaRU3UO0w6MlaixnLw+TnBZ4dxEG0QMybAKwWckVjJc0lmTx+pNJmBSusr6+lhvPv9oAMCAMB0oN2fOPmYczjhiTGJ8PpyDqcoRewY6j1TcMM92uOoHjR6Er518QPReZZ0R6s47skBixq+HA2kLiavu57Bd4Xq5nbwuVTDt9k33UM+z3dWt9sefA/Dvs08ctHqj4YiDc80SSSp/xhFjU1XUl3DVZrvrvJ5/mHuR+zVL2bId3O7dhlqK6fcMFXLqbo4jZpG1Uu267lkoD1+PrWS6JpvjH2DmHGJxUfL4/qMGJnJxcP0WChQ5Fcqdp49ul5s6d134Hvzx/x2daMPNORi4IQUDWoooFChM7cm3+baWTlNFgnOwXNznfvds4rfKOajj6rbs5ElSMLXM0J9uyiqj+E7cPZk50pYkSgC474vvjb73RTQecbDqHCiKGQ8EcnTdVR5I6hzO/c0zqeMqtmo7gMjPQbnzupWECVnDZ9RSIS3DhD/5ehxHiguHLBWGtVYg1xo+sZrSO37+HJMSXj9uGgO37hUUa6aIiFC6jno8i+WPilGDWqPny5myckjdG7ZOuNsHYEr2Y9vodCU4B4TjS9fTFf5RSHtK6hxUfZl+rzif1CwK6eMPGIFbP3tMcvyJcaXPkdz7TuUEtyPDaq0UThBSfnOQC7vPM128th2bO2IloMo5/OrsKQCAg8OMY7I9k++EzEzmb61kPma6NpODHWdCqI5DWOrZIH1HXPjs+/bxquf7niWUkwkxYvxK6BZaLs7Bx2wy0Z5B5yAL44UzJCJIL2hqeybnLBdUGboaN8wGXeTCLvIUf5qei81NSN00ZQkF0q855hyAePFVd/J6c1JpnijzVEiXhcCuu7RaYy5xbVUa9ben4BckZCeCU6okPByFPOWyEj+4bz7++yOHJM4pEHDShFGY//3peP7tDeHYgmPqS/HkxSdgYFsTFq7agkPGDMb1j79ZVw9pkx+DHvHT1D4XsnuFkuJUxQkHjsQSi1hBRatn/A7f6K2mWv2U0AgpJzijUlXutn05uuBTbgLYhUMZZMQZaLdANzutkDYk/DB0PcvKTfyz0+GgDWGldB3bHE0SB+bcsmCVw7FYKX5+RbKblmbhHFJe3YKvK4lRZ3cJqzXx8x5D+3sp/6VPTG+gATkH+wssUp4OaXAhBjhMn7gbgGQO4li5GDTS1lyMJ35YTV2gdx/SDwPbmnHEnkODOE/kvoZq8jibLKG40AzcMKRYST2mvxgS/VqKvaI/8YW+MxdCYO9d06lXpXXRtq4SfqRZVantSMj2rnogcPyb9UY6lJi6WN31UhCGzMU5qIvIM984kb2m+19Jm0cmOYckoTv9F0+w7ejg5uSsN5PXxT1CbsHu7C7h+bfXs4Hr1Htg2iT1b0lyDvL+qDqBlFgJ6fls5hzSXBNXV4qVtmkhdga2NuGYfYd7GZj0JlPReMTBckyfIGZzx/RO8dRDR2OgxlZeeMJ+AIBDxsYpHTn5sEYbrKaqhVAvYsN/nrQ/3rf/iEyKrlI5sENX+x43PPYK5YgD52DUxCjNO7r5+FO2qJn2vswnyUNbO+0xr3x2+bpMXwB43/5p/xApSrvpycX4+cxFqeP6oi57llwTZ9Sghpu+d+6KoB2jziHmRiVGDGw1XFUa/LMNPpcbuL40wTO3TyT9HNLHpMOaikvvnIePXPNkguPkmje9J7r+RM4dNb1uSqyk6Hhc4PQtLHEopYk2gCgPuIQ1B30vUofGIw6uSQw3e86xg2dOHov/+9TkRD09cBrA26TrSkBbDKICkVVOCwTiChuHMX3ibqkFq6OrhEe14GrqzkwPzQAkJ27k1duUDp+hOmyp0K1JbFDbsy1EMpG7nlNBh0+v0c5TcZHm5oYrNWrK7DN8flJUx+mYCOnr1LvWd56V5nNIBKHT2VgDUgTP8hwlF8SJ/D7z7r0Tv0cNasXcd4KEWKs2x8SBM2U1EYd+EXGQYws+1bmm3ylyXIMKXZQViKIZsZJ8rlqz+gbR1m9v2nU0HHHgJv1+I3cBEJqkifROXkesc0g2q08QLgjd7c+9A0DjHDQloJVzKKQnkx49cuO2rsjySt39S7t4AJi6dzIDV2d3WbOuSEZR5fI5nKu82JKQtEYOXArnYHBwK3pwQRLqNXPnfODgQITX2uyXX9sH6cB76UURiIlDm6Fv/Qz5nJuZ3ATROZSOKJrmHBCOT244KiMORUasZJqCv/zkZHYsVs4BZsKpOvxNGD0IQ/rF2Q63d5exz4gBuOYTR0R+NirknDvxoGTQwv6azkH2rb5XqcWc0Tno80zP6yIhhGAF0fJ8vR19g2h7BXLOoQfB3et/P25fPPetkzBmaD8I8PJDFZxFh4BIsZtcmkkZRkLtIyJGknOwKM44sdLl976a+N2/pRhyDsmXVp3ck/dKxnzq6CqlFJ5HK8l5ONGDKsOVL54MXaD2yxG7Gz49JSBgGt04aPSgVN2gvbiNxxemQ3E0Z3S+i0RFQqTSPkpwOgdubkiC2MaISDjEO/2gvu4LA8Tmn8nx6DqH4LevxzyHw8YONiikDVxnSIRTBM/KOYSZ4BzDbG4qoKyInzq7yxgzpB9OOWQ0Wz/W5SUb1q2VuEVan86UbibFof/tgnez4zBJG+RzSSURisaQ6xz6FLh73VQsYNiAlmiHE+fH5Z9M9AIpk0eINLsZK2fT7XCcg4RL56BPWl1sc9EJ40POIdmWzdSxo6uU2BWfeNAo/PcZcRBdTmShLlbyBZALpbzkRas2Y/nGbalzTzhwFCv6+tapB7Hjk5cxb9lGnHtj2rdEt7JxwYdLIwoT1SimrNwdlISJi4J6+J5D0K+Fz/jXHN5TLu4/EaMg1erIRV1eQyWMwx0XvoeNROtalFKiMlvdsL0FKzdb22wpJuX+3WWhcKJplKJFP1nelrJWSj9rjtDqxFhf1KVfgurDs3jN1pRntEQsztLLk88r1zn0EfAWNrE8U90FGMVKkQmhwjkIZjdiWbC4xVZWk2KG45ngeKqc/u217Vi7pTMV2K2tuRgtvOoLoROvpy+ZFuWy7uwuJ174YiFY7EYNChSbnP+C2pzsp6Up3sGv2NiBE6961JiqVCqk71VDpBt23/KaN7TzgePkQvCtv73sJzsWiQ8jVKV5IFtO15HEgTt20oRR+NmMSYky2d6oMHnUvqFYU0WsxDUvaFIPJIlLpX4OWs/hGAUWr9lqruXj5xDVDe7zV//8krXnlpBzUJvmCK6EaWfeTxPvyftt09Vxoq+U34MsV8qmXfVPAMCSdekNkCReejuJbJCwz8HeJA6N5+fA3Gy5UEtHGF3Mo4M7HhCV4MioQa347ocnYlBbc3Qs1Ubi3UpyKnLSnzl5D+gIYtQHx4+9cqZxZ0UI6qlKZV0mvdvgNowKE8F0dJUSL3ycOyL4zYVKVjkNOeFVzuGzv382dU7yWoLJryZXMsbKz7CT3Z4hAZPr5VOJcSByTI+PU9ZL/L/37pMO2R3OCOnU+LWTD0idFylxleHpvUiiJOcL54uSFSrncOYvnzLWS3EOltsoCd2aLfZQMy3FQqodPQihev/ldZvOKWi7PCcBc4iVdMORp15fG42BM7qIxVnJcvVnc5GsERJyhXQPghUrpTgHu1iJ0yUIIaL8toP7NWP6waOjsMC+YqWYcwgmC8dd6DoHfeGJ6hWC8SXESszTVuM1cRZUkfctK1aKv8t7qHIOegjk9PnBrlxly0333LWIt2RcGCNptZPoaJwD0voX2469SJQ6rlvPcARe7mQTFjZaN3EY9TI6u0t4ZjGbej0T1AVwg+X5+cRWkigQnCligYDY6c9fvzfqrZc19b7le2jTOeiQmykVuh9FzDmEnPu6mLPauC3N0ZYNnI16jYPamrGl0xxGvTdDazQecbCIeDjWkgMnKxQA9to1sAwaL9OJan2qC7VqfRHNd60epwT29Q2Q1koJsRKziKlKMfV4W8jOy/NNocv3CZ3CYmulmKV3i2zSk/+A3Qbhmk8cgbsufE+i3EUcTjhwZPR96fqYxT94DK/gjhXS/mMUAJJhLcI6lreIyLzL1sULyX4DGbi6QOkLshT1dZUFvnvXK3g2zHlx+7+/y35RFqjGFllEGraqJqs0feEPxEpIlak4acKo6LpN8bskJ1dQNj7dpTI6LWHh5cZQhW4ooJvGqty0FAV/8ug9ozJTrgv158C2JmzaZo7A2qc5ByK6gYhWEdHLStkwInqAiBaGn0OVY5cQ0SIiWkBEJyvlk4lobnjsZxS+EUTUSkR/CstnEdG42l5iEpytdbQr1mypTc/FdHxQWzPu+cJ78MOPHhrWS+omJOv5tekHaA45yXpyQS4a5Py2oHZ3XxQsqpLQuYlD8FkuiwRr/O0PTgAA/ObcI3HGEWPY3A1EwJ0XvQfPfvNE/OG8o3D6pN0TTnDuhTe9aBQIOOWQ0RGhlXC9JO8PvdEB4LK7Yx3HN06ZYD3PLVZSxijkDjNZhzNvjY4Rpbgu1fop6IM7L7hmdXijh7Ql6jQp1k6vLNsUlR+xZ9JMOQtUPZntzmQVK9lCvhw5Lhhvc7EQbmjiYzpH2NpUxA/OODTRp/4MJXGQQwSGYKIAACAASURBVCwL4MSr/om/z1thHaM+RFPaUFlNFSfKMZx6yO5RmUnspa5BQR4JG3Ho25zDjQCma2UXA3hICDEewEPhbxDRBAAzAEwMz7mGiOQqeC2A8wGMD/9km+cBWC+E2A/A1QB+WOnFeIG519EOLvzt0uvJSaK6xctnOHH3wRFXoFt+bAsXX91m28RhcJxDsWD2DTh4zKAoiJ+M6qk61LFOUspioCYXGRwSgyP2HIqr/mWSkevYpbUJw3dpxbv22xU/mXF4pqCEHBck++mvKRUrZa9N90q+oH7EIT5Hvw1NBXJanaTESnJsZUkcuAlHCbHSJR84MMW9qXm9bY6TWaCKVK0Lvvbbdh+LheTC+/Gjgt31PiMCRfwNnz4S937hvaFXOFBSroUzUNBNyfVLl2IlddO1eK09hhenGuZ0AaSIF1TOgbt8o5+DSvyaCg6dQx8mDkKIRwHowszTANwUfr8JwOlK+S1CiE4hxJsAFgGYSkSjAQwSQjwlgrf8t9o5sq2/AJgmuYp6gLvVejz86HkYnosc3KYOVVaYrqw7y0m5ZJtmgaH3K1lUs86BH5e6CL22cjPeXLNV0zmYxUruKJxMGVdPsXbxadM0+XXlapa1T50+qm/HrK9Pi74rzIC9Lag7/eD3zz9+eHT80LGDFf2MaTxBWlbpqBfLooPjHHEIioS1jiQW/3HLC858xT/52CTr8ajf8NPFqenDcekc1HnYVCDceO6R+P15UwEEu+cJuw9iLexaimlrpbRzIs85cLpB2/WkOAeH57r6PnHXb1JIq2LI5iJZDSiWb/ALfFgPVKpzGCWEWA4A4acU+I4BoCYDXhqWjQm/6+WJc4QQ3QA2AhiOOoGbKJGduCEGjg456VTOgTO5i8VPQUvTfvxPALEddlRP27fEnAOjqLSIldRFcf6KwKbcFTbY9wUyL2B8mfBsM2VpYlAsyxeNs3h56pITkmNQvktOCkBkmSXHB8CYoCdqi5KEhAj44KG7R7H7C0TRTlZd1A4dOxivXfaBsA3CjedOxTTNkzcKmMdc8sZtXVixsSOuw9xruQhu6yo5HeH2HZE2l+XgGz5jaP+WxO/3X/0ogCDG2J8/l8xrUNC43bIQOO6AkRi+SzL+k9z4qNfS3GSed6bFVxJNXSFtg65z2LitC+3h+/3HfztKqRf7Q6jm3fy6wvd/zcfjMDvNxYIxRD+AhCVfT6PWCmluMyks5bZz0o0TnU9Ec4hozurVqysaILebLRleQJe1ksoOcgHZUpxIiLaUeV7yuI1zsImV1OoXf+DA9Lk2nYPjBWK5BIuC2ySW2GNYP3zufftGfeuLmuqLoXpLy2pfvOWFVJtZ82GrcIqVCoEZ5pOvr0GpLFIOkgWimFvSoojqylT9XsuuOcK7vbuMmQtWR8SLq6PeO1Oo9Wg8nry4qnOw4Zx3jQMQ++LI/r9wwnjsNiipG9F1SzYOqyxE4lrWbUlbTKkWVdu7y6nNkmp96HMtQGzGLnHYd+/H9Y+/iYGtTXjXfrsq9XhjBo6I67oliT2VkDYtGnEolQW+8dfKc6PXEpUSh5WhqAjh56qwfCkA1Th/LIBlYflYpjxxDhE1ARiMtBgLACCEuE4IMUUIMWXEiPRi7ANuLYg8FsPfLrd2uWh3qpmqOFPPsMWl69tx0c3PR+X6jkkfh5S5ZrVWOkoJd9G/Jc3JSFGNSnRcIpG4nv24Xs/U3hen7R8RLn3RkESDQxadg89YExZIFsj7/fFfzUq0HVsaxYv+VoWT5K5fP5cLwijREoZnt9UZNzwOH27bfWaBjZM84/A4NldzsRD5aagYNqAl9S4UibRoAmbOVyBJZLduTytrZevbtpew/zfvw5zQSivqLyIOaaJ9+J5DsPiKU5m+05w7wPlZqCbncT2WYGjiQw7NxUIi1/TcdzbiD7PeTvTXW6iUONwJ4Jzw+zkA7lDKZ4QWSHsjUDzPDkVPm4no6FCfcLZ2jmzrTAAPizoa93IN7xW+ZPJBxA5gJnl4SBwM0UYlZHs3z16Cu15cFpVLr2PT2KT82ORbYLL8+PJJ+0ffuZ3m+yeMwqffNQ7fPW2iUi/4VBfpDx6ajmXjqwZSFxfuMSbi22iKyl1azd6wulmuHe6xmixddOiRYyOhi8I5yGH98p+vx+0zM02PwmnTJ0w7aCT2G7lLTByYudDSVMBFJ+zHcmA6vIl7+Mndl6s0vUWhQNBVHa3NhXTEU23OmsYqdQ5JLsjMnSZ1fmp/SZ2DCpNPisoRdHar4mKNOIDYuSOf6QBlDptMWVU0NyU5B71uL9IGL1PWmwE8BeAAIlpKROcBuALASUS0EMBJ4W8IIeYBuBXAKwD+DuACIYS8058H8GsESurXAdwXll8PYDgRLQLwJYSWT/WCumCNGtSK2d+YFokv5KSKrUAMDmZhPZMDml5PB7cbCcYGzHpjLW58cjEA3vO2SOl0hhJNTIpLvd9LPzwRuyqcC+fw9/Gpe6bOVcdpg7q4cKNUd3GphdfSQVkALy3dkCqXyd1VcFyTuV3XoqpxauEYS9GibVBGWqaGvmiYdEGlslA4FMOiFo5vyrhh7PGonucy45tDGgiyq+kinbamYmrube8u443VW6LfJiZHmu+WtDmiwyYKPXniKGs9jsjKzgWAhSs34zdPLI7HnnKTjgm/SsTk10PHDsGPzjoMAHD1g0EyJ9tc0BXSqXe7F1kHZ/gMIcS/Gg5N4wqFEJcDuJwpnwPgYKa8A8BZrnHUCuqtLxJh5MC2VB25uzBtxqS4x0UcTI9V373EL67Ax657OtWPCrloAMCQ/s1RrKGRA3XlXrpfm5+DOidN6wJZjkV1EpxD+rh6T3Wxkj2uPX/siyeOT5UNG9DC1EzCZZEWjzE5Loo+g2+tTUWeCDLj1ZW9NlPWQIlr94VQy3cfkp7HKnxzPcQbleT4f/jRQ9g2dS6guUipsb6hxWiy6fKEEAluxGYIwXEgJyv+Lpx42EgbwnGdFCrWJdZoOo+C8hKoJrfqVujMyWPxFSWfiG0DouscUrTBeGb90XAe0uqbnI7MmPw0IdI5GDKc6e3pSKUoZBbooB/GxrsQ72TVXZu+SHI7JJspq/oCmRXe7qmqLi4c7/DhSbsrdf0zwZneL1aZ6yF///nMhQDcuhadgMnujtl3OM4/dh9cccYh3vqQtM4h+OREHcVQVGQTPanlqtyagy9x0PNSS3CcR7GQdm4LItlKmT/fh21++XAOsn3OfFck3u90f6b74C92i+eDOs1s88g2PQJrpTQHIrHX8P7oLTQccVAXLFPCdteLFOscHJyD6YUuGIiD3g+7mPPKVL2ujUvg+k7u4Nlhe71A6uKitvPCt0/C/O9PTzgApqPY2sRK/KA4y1c9aRGHe8JItlmc4ADg0ddWh/0Svn7KQRg5iN+x2zgHXWbNipVCqzSbQhqId5YuhbQ/cUiOzdWmXMjHDOmH947fNdGGqUebS4auW+LmhCzhNgFJCWBaRGYisgT7Ih63Gdd7UMnHbdogdJfK+OlDC43tBQpps85hPBOxt6fQgFFZ4+9p8U5Y7iIOnmIlgHeu0duPzSPt9YCkWMk26blFk+MmfFIVRuMkt2BJJTZqzSH906KeLJnLglDO5vGr6PB4LrKOW0yWfGGfezut9+At4Ji2wk/VxJHIbBJcLguFgNg3Gk7i4Hmvs5h/qnORCBgR6rJsCuGgbQtnKmCNJKyWcdcsEuemy8ycQzqfAwA8962TkvUQz5n5K+KQJab79bcXluHBV1fyBxH4cWy3KKRdJsr1RONxDuqCatjBuxYtKe5xiZUAfveU3jGHY9MmJ0scQrHS/BWbEjFZ9Lq+lhqcWMmmc3BB1WG46I3nZhaA+eXjrrPdEqtGIsoF4Hj5OEc9HdyiwqUM1XeyZWGea0UKRDayrnHHGxab8lxI2IIDJupJhbuuRWW6V31uSmURzUEtVFkK7953V7ZcGiioAe94rir45MNbxCfIa7l3bhxTyTTlTJzDgFY9mkE8H6buHRsBmAieOsbTJ+2ectjUdQ76VHJ5vtcTjcc5KN9Tsv9w6nDiHBVykfXhHPjgcjzHkurHIFYqC4Gzr59trWvbcSX6ZsVK/IT023xKYiPg5jLsDapHXQ6JKmyBzHT4EDCXiEVfR4cNaMEvPzU5VU+/13qY9ES/hSDj38PzV0Xj4McXHHDlSuC87TlIT3994eU8rIsK59BdFpG1XKxzSA96xpF74F+npvOUBPUDgulSIMv3VJXVtzYVsMew/vjQYaOVegFmL14blRk3PgyHD6Q99gPOgdE5aJuMEw8ahXc2bEtYJu4zYpeUw2YQbDAmrroOJ+ccehAc2ykRcQ4u4hDqHGwBs/Q2E+cb2tcnp9kJTqR20voiyVvApPuMZczKONjR+YmB4qBo7ontK1YqkNshUUX7djNHt++IAYnfrvhPPiHS9cO3ff5drImtfr1lYTYxLRaCufqdO+ex58ZtBp8ybpdMZZmq5/mmyzzY+sZn8l7pSK8FRedQKotovtp0DiMGtppFZKFo5yOKs51tQ6PuuI87YAQe/NL7EmFsOJ2DUZ8GXqyUGishpTMC0nOgWAgIhuonwekPZKgPeS3qXGttsofWqDcajzgo348dz3tZu+Sz8iXgsj/p4F6ElK5DWVBVcHJ66W2q76TTnEN6LNx1sYmNjC+QG2p7TlGH1qAer0fmYhjcrznQOTBttDFRO7daOIfpB++W+O3jGe7iHEzpJFNtRX3KBbVs3CgUKLmLdFm+dZXKOOHAkXj8v05g6/nqHFpDcdg2C4FV25TX0lUqK2IlyTmkz7Fxi9J6zZW4Sbb/2MI1qbJkvfS5xidp4By4vvWMjQDw3Q9PTNSTXIDKOXzgkLRzqYzPJPUOqjhvYFuzdaNTbzSeWCl8ntd9anIiQQyQdoIzITZl9eAcuLJUYXqB3nvXAXql6NySSCdi0Yec1ZQ1SRv4t2Srx0Tl/CbMdZPjmXFkUtzwvdMOxseO3BPfv/sV4yKuOxQCdrGSPi4fayWXsl4/alzItXDmXSVhTDGqKnsBs1iIIuIgrJyYv1gpqOcjmisaOAc5anbXb2kvEO0Ia+4HWQ9AQtFrE6MmOQeDyBRu44RojOH3khDYd8QAPPTl49i+dc6BQ8Q5hGvJuq3xhqp/SxHtTPiQnkLDcQ7y0Y4d2j8VFlpOL5dYKau1UrqM5xy4seiQi1UlYiXeKib49DFlVfH1U9KB/dR+1fGpirvkeOLvXEyetuYiJu81NOWIpoLbZdp2W3orbp0DOcOF62Mzm0uq+pjAA5/LsAekcyCYiIgs7S6Xrboyf7FSwDls7fTgHAq8ziE2v/XTe8XHpM7BTxyZ2Nyx71m6zNQ2Z4j30xnpMOcERaykKOF1SM7BJTqVz7W7LPDQqysTznMBccg5hx6DfLD8ou3XRsQ5eIiVvBzHtLElCpm+t5fcYiV1zfnWBydgSL90Jjd1fOpivscwu+PNUXsPw/nHmoPkAcBba2OvWN17W+8bsO8oOaW+BFfqm1lrxcYOJ1ewYOVmLFi52VrH0wcuJT7s6hbGBV2Pl2SrBwRWLTaO19fPQRJbnx2rqkBVOYcBLU3YbVAbLgk3EH/4t6PwiV/PSoyXg9z4uCScsglVHm/lHJRWthiIXqBzSO4CVG9riaZinJynVDYv/lJZ70voykJg1pvJeKP9W4pe4r16oeE4B/mouGcqJ7eeO1aHfNG8FNIeY5Isv7oYvLF6K1uXwkVDn3Nps9z499679sdHJ48FB9WCZmBbE06aMAr7hzmwjeM17GKBeLJfds+rUZlJhqwO2SaLlrt37kXzD8YXQG3iS7e+4MX9qTiAuTe+YiVpXfOHp4Oom10WzqGgWa4YdROKn4ON4/VV/nMi07OP2cvYptR/qaasxQLh6a9Pw2mTxqTGbhuj1Dn4xLsCkmaeNgdP1arpxSVpPxVZV++WEwmNHdoPb68LssqVhZkgFwqBWMm1cVBD9ejPqLlY8PL2rxcajzhIzoFZtiVb7LItzqRz8LHwCZ/C9+95xV4RUtyQ3pHoCkd1onFhOPTxCSHQXRIYZ3HXl1YktoWGO2QTndjOU4+VFcsdFWoQwW+eepC5kRAqp/Dk62vxwf993HmOip/+a1rUoC8Apme+KbQouu25IO9Vd0kkEsaoKIYyawnTPVR1DjaxEnfsps9MTZXJaqoPD2d5BQTcaakcR1E19d/k+5zDTxfBj7glRd7HNRsHyPQxHOGIfLrVQW3NaN/ejVWbO/Dgq6ucfiq+0XIDEZXWRoGcfjj1RMMRB1s4gkiX4KDWTZETXFDv/GP3Mdb12bDJ9lzWPUBsWqkr7VJiJZU4OOS8ANDRVca2rlJKD6NCcgA2EQX3MnBKYyCwxtDHwY8xkL+/tHRjovz+/zwW+ynmgaMMoSxU+IqATBjB5OLQFfimS9HvTXe5bLzfelIns1VT3JaNaHPnczvjIDZScuNjalaKle6fFyiGbdcStWXhpWMHPD/OYXuCc+B0bMHn+q3u90pP9rPrLnzwRnnNX789SMgz952NbD3ppyKf4d0XvYevp4iV0lEO0n4PPYmGIw42sZJuOWBCtLsKdQ4XHLefsa4PM+8rD5Z9l8oipSRNK6T92pfnyV25aYcKxCkbbcSGk/eb2txtcLyYWxcNQ1hsXfzlIzqpZiP2wUNHs4ma0pwDf76+g7Xt9vWghGadgx/nwJpUm5SpRF4iU2mRc8Efn7OOUbWUsm8Cgk/VP4YTacl6W9R8Dky78t68s2FbVLb7YH4DoXMOI5hozbLNUtktNSiGc9YVOFE+A06sxEW97Uk0HnGI3mROrBSUdYQv8T4jTOakhKYCRROkaJHB+4iVshGHYEfp4hxU2a6tfX14zZa6Mtm7TW7MTeYvTOOJpxpiwsU5+OYBdsHl9GbDgbvxupi0nwM/kA4tUGN3yWKtpF2MaVeu3jen86aHF71sR138bAubj14koXPwEEmqc2js0LSYU75Tryt6OY4wcV19YVo6xLtEIrSOkVsKfBFsmyggVkhHIdcN1VWd36JVWxLHVA/03kDDEQcJXiEd3I4l69oxYfQgIysISKshczpPCZ91n3upPh3m6E23F4hY9Emjz9WkzsHNOUg0W+yyJedgu95jwlSlUk59yQcOZJ35gCR59nGOcsHnXttozF0Xpp/36VqIcbZNz3EMbEsaBwbmn6YFNfnbyGFAfc58vxIy33NcvzrOQd/ZmpJjqXoVu27JV6zEj0WHPrfPffc4nHGEyTCDrKF1JJoKhYTy3QSpkI6SQpkIcVh+94vLcc/c5YljOefQw4gV0mnISdxVEjhwt4Ho32K29FUnRxZ2/ooz0klTXLGcVBQKfJrQTduS4hx1ochCHGxjkToHq1VMgXDs/iOie2nSN+h92xYNm59DEu77aGtnSP+0ue94RXRluo+pBccwjH+ZEjj5nRp6ynZ2lY1OUqk8zEbxU/zd5eimz0XTsy4WyEuJWyTC5o543i1cuYWv58k5yGMmImNrw6fs2x+cYJyPBCR2DqYpLgNfmgwJJKRC2jcfx4KVm1LHcuLQw5BiBW4XqLLuvl7Srrr6oRlMCs6sYiVuVzdUy36mXp+VeGm/bYu5K6xBPMZYhGJrT30EduJQO84ha1KWxHM2DHKX1qZEpjSTWKlQoEBUGR7u7C6x0VuB9GJiCtWi1nPpXPSjpvpbOrvx9Buxzb2JY1q+sSPxe8M2XvGrEi2XVRrgthbkRIymIJUSLU0FB3eqBeU0clXSp8MhVipoYiWLmIrDWZPHpsyZexqNRxxsnIPKDTh2BknzPHNdn7AFWYmDGjvoAwfvhke/ejwm7TEkVY/7zrWnwjZeKXIqOV7eAlGkt7ERlARxsFqx+CWf8VFI23QOXB++Zpj+XFC8CnVYOAedEJl8T9TH5ZqzKaOFKt9+XSn7sSl8tFVVJ+fykAbcz3rN5nQEWq5ddbPX5ghjQfBLCiSDDbp0DulkTSbOjy+/8qzD0LQjcw5EtJiI5hLRC0Q0JywbRkQPENHC8HOoUv8SIlpERAuI6GSlfHLYziIi+hn5aHErhM1DWp1MvjkdXJAvLBHw3x9Ji5QA/6BoQb+EbcouUghgT8Y3QW3TtmiQLtu21JULvWs3E3AOIXGwcQ4KQbAt2pzOgbO997mNVs6BKWvy5hAV4mDpnxAvfp3dpUQU0UR7Sl/fPPUgDGjlRZxZOAd9ynonADKUq+vjiQeNxIkTRrH11HvoMjwA3NF8D2AMA7hno25MWg0cmgRpRg/GnX5Yz0usVBZRWG+XWEl9pU48aFSijd5CLTiH44UQk4QQU8LfFwN4SAgxHsBD4W8Q0QQAMwBMBDAdwDVEJJ/YtQDOBzA+/Jteg3GxiG2V0g9LnewuPYCvnkDuMKaOG4aPH5UWKQFuKxMVREkRg2lRVediFp3DAIueRXIOPglypFjJFnhMXayWrNtmrqfpHA7bYwge/6/jU/XUFKRfm34AHvry+1J1bDoH7ljRU9RY9OQkpeXVS0s3YPHadqMoRl17bO2t3bI9+u5asPR2fOedqXuVy6zFvZF4nsm2p2Kv4QMiw4d4LAxxaIrLXAHwCEhsukzjlOIiL4W0iN8V016SY0Bk3LKd0QnuNAA3hd9vAnC6Un6LEKJTCPEmgEUAphLRaACDhBBPieDt/K1yTs0hFwDu2asTwsUZ+IqC5KTNIjqyYeWmjsTu1yfPgY8TnMQubWbi0OrNOcSNqo5uOmyiJL09dVf3ial7si/vvopD3P4jB7IJamxDd3EO1kXfc3dMFIzhpiffAgA8o8XT4dqzrflHK4ukUyGt/a6ec1C4U0vflYjmXLCZbkuooh+TbkcdV7sSd8mukHaLlYoaF5RFrCTbLhZoh072IwDcT0TPEtH5YdkoIcRyAAg/ZejEMQCWKOcuDcvGhN/18hSI6HwimkNEc1avXl3xgE1Imn/a2/EmDh5exRzeb2DRn1i0NvHbJOet1JR1F4P4AohNWZ3hDZR7N84Qejzo29qMUi8pVvIRsZh8T6wKUaZMzRfxxmreGifoW23HIsYLr0XWN+32E+Iiy43afUjsrJVV5+A7J01EMSEu8uQcfHQOEtd+4ghzXd2ai2lXFWmagj/GoESwQVdAvRufXGxtTa4f0vLKJVbikvpIQtRbqJY4vFsIcQSADwC4gIiOtdTl7o6wlKcLhbhOCDFFCDFlxAg+UY8TFp2DOt9cnIOvM1VzpHPwJw5zL30/3rUfn2dXxS8/Odko501ei22xSv4eNcj8EskdjW94A8CuCPS9J0SBnH5oaGp6yiHpaJmAds2Gtr96Mh9q/KIT9mNzaPRTdpy2RT/Rn3V3DABxqASTTsZXFKPuiN1JcrTfnsTB5gwqYY/r5OkhrR3kkuOY2nFxDiMcxIEoqeswrcn6bv5vF7ybH19BLvoOa6XwHqq6J3nuDq2QFkIsCz9XAfgrgKkAVoaiIoSfq8LqSwGo29yxAJaF5WOZ8rrAZsqahXNwxfiXiMRKjvfwfz56aPTdJoo59dDRSj3zLl99WeymrPGxAS1F1iNVosWTOKj30RaryZdeyjANRIRPHMWLlAC/ax42oIXljr78/gPYdvspOhhrNFpP0UlA6GKltDEqK/HfdaiydDVMBN9msqFhBudEFdd84gi815Ax0XcD4kvosvDWOvHnNgPqvXWanmom4qYdu64H0K0E9fFIjsBEiGVxU8KiS37uoGIlIhpARAPldwDvB/AygDsBnBNWOwfAHeH3OwHMIKJWItobgeJ5dih62kxER4dWSmcr59QcNlNWdX65kmz4snu+YiVfFl8VNw03BAcD/MUS6qEjmDzBKiKFtIe1UnSOR3hvF4qhcq/LEm5Cb8/KLXn1GqB/ixriw+9abO3LnAV/fjaQpJqJg58oRl1sl4ShpM2Vg48TDxqJxVecin4tdjk8ABw0epDxWOJ+W825s1+LC/qc5gi3Sjhdm72WpkJiIbaFDOnyWLDjnC8BcWh2ZPLbzoQr6W2FdDXJfkYB+Gt4cU0A/iiE+DsRPQPgViI6D8DbAM4CACHEPCK6FcArALoBXCCEkCvw5wHcCKAfgPvCv7rAFnhPnRCu/NAbDVYmOnzFSr7EQR3j8AFmVjlhyuqpc3DRu8hD2mUymejbZsoaQ08RmqgXipW6HcHlfHeyWahDP8/4T+riYyNgulmuMWS3p4JbhWsZiZ+L/w3wzS7n8pqXcBFOX+jd6XGJgORzcL1funjPNG1NTqhcPSCO02ZOBxt83vfyiqhM3oamXlZIV0wchBBvADiMKV8LYJrhnMsBXM6UzwFwcKVjyQJbPocsC6Vv+j65OLosQ3zlv+rLOpQJ9yDha8qq1nMnWQk+nZ64ymFfwvSe8WYdixQrlZRUlK72asU5qLtr292hDH2r9/kozSSTa8N3R+0KMRKLK7yaA+AmdBIuU9FoDJZq6mXqpqo69Ht8z0vLcfXHJhnruOZsijgY6tv8dri+I87BI5y53jenkH5zzVZ86vpZuO3z7/IKUV8NGs9DOtI5pI/VyNo0gWZPsZJv1wm5uuXF9Q3vob7gLuIgD7v8/6LJTXaipz6DUx3Kx7KQmdP82rNyLBl2qKpYyXZ//Dm15O+vn8InKFLHWKtp6VKCS1x0QhxF1zc0TGuz31JiIzaq5//N5x9tbUdfvF0Z0xassKd61ZX5pnGu2tTBluuQuoYbnngzbI+/j/I5c6I3TiH9+6ffwtL123DnC3VTy8bjqHsPfQx2nYM/5+ALOSlcnIHveuXrfKe+PL4xkVwcbMlhsx337ddewsDH0mZgIVKGEPZFP6ln8evXBVWsZJsT3kpXSu4G+xns71Vi40vM/uPE/a3HZTMu4qAGErQSY+UtMnl667DNxRUb/RZeID0HTdEHJOa8td56XOV89hzW39je+vbtbLkOPQih0YgiLFbFRxQdC0SQKkcoNyu+kotq0HjE8nsyLQAAEK9JREFUQX5xcA62HA1Z0BTJ6WvSXAbdRHoMLrjEEq44MRK+zm3+pqxxfgGbxVBy9+6n63AhIVay3B9/y6ukpZs5B4LStl/TmLC7WXkc9B205BIB+c4dtV6bJ+dgI0ztDj2fCv2+HbXPMO9zOajjuvwjBycSUamwRWpW4XstnGWiqpAGktaBcj62d6WTatUaDUcc5PaPW8DUxaq/w6NS4vun21UlJYcTjDYsJ7yJQwXUyBbCAgCm7j0MR+09DN9w5Gru8rTz9XeCi605fMVKvqI0ALj3C+811lV3ujZOSA+Zbuwb5JU0fp2S2rJWXKUUlbj9IeJ2fP1UfDlaW98d4W745Im8744K/X1yXZMrtIi66Nvmjq+V4lYmIyIHzq9IxjuLiIPSpzTD9jFDrhYNRxzs1krx95EWZ7Bkg/bJcu/cwArh0dcq8+jWUYlVky9WOOSp/Vua8KfPHmOMECrh67jjbcqqcg6euoQsHum2HTcR4VsfnADAvjD4LhpE8LJ2OU5JzOP7KF2cmFzAXWKlhO7Gsugeu39sROC6fLkI2vqWsY0+fBgbICGBVDIkx+JvSjgloXqa24xH9hpu9vhXceEJ5tTBKjifGxvnIJ0eT7Ho6GqFxiMOFp2Dulid++69/drz7Ndnt+gD33g4fslx6gNXPH4J3+VbXfRcO0AJ3xwWh40d7Gxr+sG7oaWpgBlH8oETAWD6RN5rW4evKeTuStRZXzGdC1JU6iIOvkT7hAPjHb6LGMvpaOtbxonaf1Q6JpYOXbzjei9sDqNAknjYCM1Fnov+EXsOZZNH6eCdccNxMMRBvluu2E61QAMSB7OHtFrke/NN3qM6fE39XHDtkCQqMY/+8kl2haYvuDgxHHx1Dr4e1ypsi1XCdtxjDGOG9MNrl32ADRUdjdGTUyGK78/4ke5F0HOIXoislYp2kWklQSJ9z7ERh7OP2QuzvzEtkX3PhLFDkyHbTXPp/v8MI/o43ofmhIey+VqyLMqSYJ1w4EhrPf3WRVaBlCQOj762OlKI1yqQpw3VOMHtkIhDdqeRRRTzyvdORldJYHA/9+4A8LcY2tXi9Qz455GQVg0HWhY0HRdZkq9ngS9h8r3dahwrX9m27eVRU2D67OJriYIiIvvMe/y4U9fi4ouusN9+Lf5iJV+4iLZ0GrXJyokIIwf62e7rnIOJM5DzwCX2U9/PWqWTkT0ePMbOnRYLhHJJoH9LEe3bSynz91JZYEP7dpx9w+zoHF8Ouho0HnGwBt7zv+G+VgsSvjveo/a2O/8MD9OBuoY6oLUJL377/dYQ3PWCr0jL934v3xDrQnw5JxtxUAmCT67kWoIoFjH6EjpbrK0skP1mUUj74oMOGfgvPzUZtz27FCNr5LglF8f3jt8VN5471fi85bW4NiwqR+ASQUl8VM8drmHd1mCXv/eu5nhlgByjwCWnHISPT90zuhZVIV3qTl6A73pSDRqPOISfvId07fv7zocm4Lt3veIfXdExhjFD+uHwPYfg+APcu8nBHjLPesBbOevZnmqt4ZN2NahnsziJv/c050CKzqFWcuMzDh/jlWtY9tvi8EmohDjoOcx1HH/ASK856wvf0BjykCuKstoelweEw8csIV9U7DYonbVQRVOB0IkgTTHnvMqtHb4bi2rQeMTBopGuR3ZSqWRzeZDKnZDTBr1A+Ou/82GC+wp86aDvIqQGH/Nlp206gKl7D8PsMMlOz4uVUHPicJUWNsIEH3NgoD6bpFqjyTNmmdwEOuOGNWW/6IPH2P1KJFyK6ShEtzYfiorOQScQPUEcGk4hLdFT4TNkApFBDtHAiQeNwgXH74tvh2aTOzIk53DTZ6baK3re71pzDr/59JHR986e5hygiJV6QG6sQsr9Bzn0ZHVM4V4z+M4DeSku4pCFUI8e3IZBbU3eomXX/ZYcgj5n5fzoLonUJiYXK9UBvqastcIBuw3CmCH9nI5jTcWCMRHNjgbbPVbhS4zVF6MWOgc1k1xvKKTlLrAnlIoq/n979x5jRXUHcPz7Y1nZwKIIAlXeIkUpVRSkNYpUrY+trZDYJhgjVGqqiUatthRaq020CTatKWoTSyxEWlutqU2xxRo09VlbBQUREUGlFbsRTRF5RJ6//jFndod779yZWebOg/19khtmD+eend89u3vmnDlzjr8r3+iQzXt8WcyEOVRxr/T9odXgPii1JPlj+9ycc2JPYYfo0YCeHT2Hg2PylyTZs/9A7BmAaep+jUPMzX7S0tqrJy/MPTf1coss7jIbcT/v4NV93Cu8uM+DpHlD+qmbp9bdZhWSTZeefeYoPvk03tLwcXzvwrFccspxnDCo/gy2ek9FF4X/+xtVy0e2NLPq1vMjb+of3buZ5iZhblv9izhIftUe97mSyp6D36js3nuA9m31Vy9ohO7XONRbeK/4vxOl0Nk41M/XuQR4/XwH9RxSmMoalOZy+XFuZAYvSqLW6r/1a+kOMbY0N3FKyM5lQVkMWRwqf6mNOM8PRT0dDd7sw7fuaGvIkFrU7LAt23cDnVuKdrzPbxz27WfxC5tSP68oxf8pSFm95TN8fSOu/g43t1x8Ekui7g8k4K8z3yfic6xcJiBM8Oo+7nr6Ub/ka358Qaxy0rb8jQ86jrftSq9XkJfThveL9SRw2j5290/S/N6NutcS9xmnzVsP7h109Bz2Hcj8/hR0w55D2/jP8NnBrTWXGFZ3gVqC+3GpumrK8amWd8f08Zw9ZmCsq1SI/qUM9hzCVsv0bZp/cazv6Q8znBxj+YxGqffEdd7mtZ3IiAH15+cDPJrTzLnzThrE8P69uXrq6Fy+fxJxn56vfGK+l1tH6buPrKY9sJz56SPrb+eblm7XOIwY0Cd08Sx/uumMyeFr6JhofVuauXRi/QeEgob0qz8PPNg4DI75BG0cT940Nf4Ciw1Qb3/mvBX9j+6gvi08O+ecvE8jFU/dPJVX/r2VL487eDVav+fQXrHPxf2zTicLhWkcROQiYAHQBNyvqvOzPoeW5ibevP2i2N1Ac2hamptYMGNC5FPhcy46kRsfXgV0bSnyMCfEXNuoEaJmz5juY/TA1pr3qyqnwF5/3hhuSmn9szgK8VdQRJqAXwJtwDjgMhHJZcJ/S3NTqn+ATH3TJgyJHCqafmr0Es5lc1vKN5tN8Tw35xwevOoLXX5/5R7xcfa5SFNReg6TgY2q+g6AiDwETAPeyPWsTGE8fsOUXG58pu2Y1iP4aMeejnX5zeFrWP/eDOsffd8mTOUWsoNTWpcqrqI0DkOA9wJfbwa63uSaw06Rx+eT+MPVZ7BsTXu3mxFnkhMRNs2/mL37D/DuRzs5pjXb+2OFGFai9mMHVZPAReTbIrJCRFZ8+GE6O6sZk6XjB7Zy3bljSrFEhSmG5qYekbsvNkJRGofNQHCJw6HAfyszqepCVZ2kqpMGDoy3yY4xxpjkitI4vAyMEZFRInIEMANYmvM5GWNMt1WIgU9V3Sci1wFP4E1lXaSqa3M+LWOM6bYK0TgAqOoyYFne52GMMaY4w0rGGGMKxBoHY4wxVaxxMMYYU8UaB2OMMVVEozZXLSgR2Q6sr/Ffw4H/xCjiKGBbDvmS5LVYssuXJK/Fkl2+JHmLHkuSmBsZy1hVjX6qTlVL+QJWhKR/GPP9C/PIl7BMi8VisVgOk1gSxtywWML+dla+DsdhpY9j5nssp3xJ8los2eVLktdiyS5fkrxFjyVJzHnF0qHMw0orVHVS3PQysliKyWIpJosl3bLL3HNYmDC9jCyWYrJYisliSbHs0vYcjDHGNE6Zew7GGGMapPCNg4gsEpEtIvJ6IO0UEXlRRNaIyGMicqRLP0JEFrv01SLypcB7Jrr0jSJyt+SwoH6KsTwtIutFZJV7DcohlmEi8ncRWScia0XkBpfeX0SWi8gG9+/RgffMc5//ehG5MJCea92kHEuudZM0FhEZ4PLvEJF7K8oqVb1ExFK2ejlfRFa6z3+liJwbKCubeok7tSqvF3A2cBrweiDtZWCqO54N3O6OrwUWu+NBwEqgh/v6JeAMvI2FHgfaShzL08CknOvlWOA0d9wXeAtv/++fAnNd+lzgTnc8DlgN9AJGAW8DTUWom5RjybVuuhBLH+As4Brg3oqyylYv9WIpW72cChznjscD72ddL4XvOajqs8D/KpLHAs+64+XApe54HPCUe98WvOlgk0TkWOBIVX1RvU93CTC90edeKY1YMjjNWFS1XVVfccfbgXV4271OAx5w2R6g83OeBjykqrtV9V1gIzC5CHWTVixZnnOYpLGo6k5VfR74NFhOGeslLJYi6EIsr6qqv+HZWqBFRHplWS+FbxxCvA5c4o6/QecucquBaSLSU0RGARPd/w3B223Ot9mlFUHSWHyLXff4R1l39yuJyEi8K51/AYNVtR28Xwi8Xg/U3id8CAWrm0OMxVeIuokZS5gy1kuUstbLpcCrqrqbDOulrI3DbOBaEVmJ10Xb49IX4X1YK4BfAP8A9hFzj+qcJI0F4HJV/Twwxb2uyPSMA0SkFfgjcKOqflIva400rZOeuRRigYLUTYJYQouokVb0eqmnlPUiIp8D7gSu9pNqZGtIvZSycVDVN1X1AlWdCPweb8wXVd2nqt9R1QmqOg3oB2zA+yM7NFBEzT2q89CFWFDV992/24HfkdOQhog04/2gP6iqj7rkD1zX1x+a2OLSw/YJL0TdpBRLIeomYSxhylgvocpYLyIyFPgTMFNV33bJmdVLKRsHf6aBiPQAbgHuc1/3FpE+7vh8YJ+qvuG6a9tF5IuuOzkT+HM+Z3+wpLG4YaZjXHoz8FW8oamsz1uAXwPrVPWuwH8tBWa541l0fs5LgRlu3HQUMAZ4qQh1k1YsRaibLsRSU0nrJayc0tWLiPQD/grMU9UX/MyZ1ksj7nKn+cK7mm4H9uK1mt8CbsC72/8WMJ/Oh/lG4q3Uug54EhgRKGcS3g/E28C9/nvKFgvejIyVwGt4N6oW4GbKZBzLWXjd2deAVe71FWAA3o30De7f/oH3/NB9/usJzLDIu27SiqUIddPFWDbhTZTY4X4ux5W4XqpiKWO94F0o7gzkXQUMyrJe7AlpY4wxVUo5rGSMMaaxrHEwxhhTxRoHY4wxVaxxMMYYU8UaB2OMMVWscTCmAUTkGhGZmSD/SAms1mtM3nrmfQLGHG5EpKeq3pf3eRhzKKxxMKYGtzja3/AWRzsV7yHFmcBJwF1AK/AR8E1VbReRp/HWvzoTWCoifYEdqvozEZmA9+R7b7wHl2ar6lYRmYi3htYu4PnsojMmmg0rGRNuLLBQVU8GPsHbY+Me4OvqrYW1CPhJIH8/VZ2qqj+vKGcJ8H1XzhrgNpe+GLheVc9oZBDGdIX1HIwJ9552rmvzW+AHeBuvLHcrPjfhLYfie7iyABE5Cq/ReMYlPQA8UiP9N0Bb+iEY0zXWOBgTrnJtme3A2jpX+jsTlC01yjemMGxYyZhww0XEbwguA/4JDPTTRKTZrbcfSlW3AVtFZIpLugJ4RlU/BraJyFku/fL0T9+YrrOegzHh1gGzRORXeKtm3gM8AdzthoV64m3EtDainFnAfSLSG3gHuNKlXwksEpFdrlxjCsNWZTWmBjdb6S+qOj7nUzEmFzasZIwxpor1HIwxxlSxnoMxxpgq1jgYY4ypYo2DMcaYKtY4GGOMqWKNgzHGmCrWOBhjjKnyf3Koe5mlpfRuAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXmcHGd95/95qqqr+u45NTMazeiyZFu+ZFsY4wMDtjGEbGwIZE2ywUkcHI7NcmSzC/klr4RNeIVkE8iStWEhJEBMAoYQMGAcjLHBxgeWT0mWZd3SHJp7pu/qOp7fH1VPdfXd03dPP+/XS68Z1XTP1NRUPZ/nexNKKTgcDofDcSO0+wQ4HA6H03lwceBwOBxOAVwcOBwOh1MAFwcOh8PhFMDFgcPhcDgFcHHgcDgcTgFcHDgcDodTABcHDofD4RTAxYHD4XA4BUjtPoFaGRoaotu2bWv3aXA4HE5X8eyzzy5SSocrva5rxWHbtm3Yv39/u0+Dw+FwugpCyOlqXsfdShwOh8MpgIsDh8PhcArg4sDhcDicArg4cDgcDqcALg4cDofDKYCLA4fD4XAK4OLA4XA4nAK4ODSIBw7MYimutvs0OBwOpyFwcWgAsbSGD3ztOXzz2al2nwqHw+E0BC4ODSCW1gEAq0mtzWfC6XUSqg7NMNt9GpwNABeHBhBXLXGIpbk4cNrLrXf/HPc8crzdp8HZAHBxaABZcdDbfCacXmdmNYXp1WS7T4OzAeDi0ADi6fVbDk+fWMKhmbVmnRKnR1F1EymNu5U49cPFoQEkarAc/vg7B/E3/3GkWafE6UF0w4RhUqQyRrtPhbMB6NqW3Z1ELW6luWgaPlls1ilxehBVtyyGtMbFgVM/3HJoAEwcolW6ldKagWhax1I808zT4vQYXBw4jYRbDg2gWrfSE8cXsSnkhSJZmryc4OLAaRyqbolCiosDpwFwy6EBxGxxiKs6DJOWfN3H/u0APvPQq1iwK6lTmsH9w5yGodqBaC4OnEbAxaEBMMsByLqYipHMGDi7ksR8NNtmYznJrQdOY3DcSnzDwWkAXBwaQNzlTiqXzqoZJs4uJx3LAQCWedyB0yC4W4nTSCqKAyFkghDyCCHkMCHkECHkQ/bxAULIQ4SQo/bHftd7Pk4IOUYIOUIIucV1/EpCyAH7a58lhBD7uEII+YZ9/GlCyLbG/6rNI65mH8ZycYeMbmIlqeHUYsI5xi0HTqNglgMXB04jqMZy0AH8AaX0QgBXA/ggIWQPgI8BeJhSugvAw/b/YX/tdgAXAXgLgHsIISxn83MA7gKwy/73Fvv4nQBWKKXnAfgMgL9qwO/WMtxupWiqtOWQsXvePHdmxTm2nOCdXDmNgcUc0poJSkvHvjicaqgoDpTSWUrpc/bnMQCHAYwDuBXAV+yXfQXAbfbntwL4OqVUpZSeBHAMwFWEkDEAYUrpk9S6c7+a9x72vb4F4EZmVXQDcVXHYEAGUNpyMEzqBKsPTUcx3ucDAJ7OymkYzK1kfc6rpDn1sa6Yg+3uuRzA0wBGKKWzgCUgADbZLxsHcNb1tin72Lj9ef7xnPdQSnUAawAG13Nu7SSh6hjr8wIAYmpxyyHjelgzhokdwwGIAsEKdytxGoRbEHgWHKdeqhYHQkgQwL8B+DClNFrupUWO0TLHy70n/xzuIoTsJ4TsX1hYqHTKLSOm6hgNW5ZAKcshk7eTGwl70e+Xea0Dp2G4LQced+DUS1XiQAjxwBKGr1FKv20fnrNdRbA/ztvHpwBMuN6+BcCMfXxLkeM57yGESAAiAJbzz4NS+gVK6T5K6b7h4eFqTr0lJFQdm5nlUEoc8nrsD4cUDAZk7lbiNAzV1XCPiwOnXqrJViIAvgTgMKX0064v3Q/gDvvzOwB813X8djsDaTuswPMvbNdTjBBytf0935P3Hva93gngJ7RLImqGSZHMGOj3y5AloWQLjQJxCCroD3i4W4nTMLhbidNIqmmfcS2A3wRwgBDygn3sjwB8CsB9hJA7AZwB8C4AoJQeIoTcB+BlWJlOH6SUsjv1/QC+DMAH4If2P8ASn38mhByDZTHcXufv1TISGctSCHklhL1SRbdS2CshmtZty0HB4XPlPHQcTvXkBqS5OHDqo6I4UEofR/GYAADcWOI9nwTwySLH9wO4uMjxNGxx6TZYGmtAkRDyekqmsjJx2DEcxAtnV7EppGAgwGMOnMaR41bK8GwlTn3wCuk6YdXRQUVCqArL4bU7BhBUJOwYDmIgIGMtpUHnM385DSDHrcRjDpw64eJQJ6yXUlCREPZ6SrbPYDGHa3YO4eAnbsGwbTlQCqyWKZzjcKqllmylv/j+y3j/vc8265Q4XQwXhzqJ57iVKlsOHjHroRuwC+e4a4nTCNyWQ7XN9w7OrOHVuVizTonTxXBxqJOEWqVbybYc2CwHAOjzewAAq0luOXDqR9VMyKJ1f1VrOawmtYJMOg4H4OJQN7GcmEMZt5K9q5PF7GjQgGLlAyQz1Y8X5XBKoeoGIvaGo9ppcNGUVlCgeXopweNgHC4O9eJYDl4r5pDIGEUfLEccXJZDQGbi0JvBQ90wcdvdP8cjR+Yrv5hTEVU3EfFZ4lC15ZDSctxRZ5aSeNPf/hTfe2mmzLs4vQAXhzpJ2At7QBER8VmLfbSIayljWK9zi4NftqyIcgOCNjLRtI4Xzq7i+dMrlV/MqYiqm/B5RCiSUJU4ZHQTyYyRYzn86OVzMEyK6ZVUM0+V0wVwcaiTWFqHLApQJNEx6VeLVD1rulXwnWM5MLdSj4oDs7p4tlZjUHUDiiTAJ4tVBaTX7OvuFoeHXp4DAKzwOFjPU02FNKcEc9E0Hjgwi9GI1VeJmfRrRRY71SjMVgooluWQ6FG3ErOYil0vzvpRNROKR4DPI1ZlOaylrE2MblKYJsVaSsMzp6yWZrytC4dbDjWiGyZ+80tPYymu4u9u3wugvDiw3ZniCkjLogBJID0bkE6UEYfvPD+NI+d4iuV6UHUrW8nrEZHSKgeU3VlyGcPEI0fmYVLA6xF4Bh2Hi0OtzMVUvDoXxx+8+XxcMWlNSK1GHNxuJUII/LKIhNrblkP+QnRmKYmP3PcCvvzEyXacVtdiuZVEeD1iVdlK7vtU1Uw8fWIZ/X4Prpjs55YDh4tDrbCMJCYIABC2Py/WX6mYOABW3KFXLQeWpZV/ve59+jQo5cWB60XVmVtJqEoc3KKsGgbiqo6BgIz+gMwtBw4Xh1rRDCvA7HEt9mUtB8OAKBCIQm4PQ2455F6vVMbAN56xBgnyoOj6UDXTCUhX07LbnQiQ0U2kNQNej4h+P28lz+HiUDO6aQeYXYu9ZdILRcVBM6hTveomoEhO2+9ew52txMZ3PHhoFmspDSNhpWjWF6c0zK1UfUA6VxxU3RKXfr/VEJLNPOf0JlwcakS3LQcpb8GP+DwlYw7uTCWGXxaR7FHLgYmDYVInY+v0UhIAcMPuYSwnuOWwHtji7q1WHFzimzGylkOf32oIWaran9MbcHGoEdaPRspb8EuJg6qbkCWx4Hiwly0Hl+uDWQmxtA6/LGIwaFkOXTIQsCNgMQevp7o6h9W8gLSqm45bCeBuvV6Hi0ONMMvBI5S3HM6tpZHM6MjYu7p8/LLUs+0zEq7iP3bN4mkdIa+EAb8M3aSI9WiB4HrRDROGSR23Ulovncp6YGoNL02t5rqVbMuBuZUAXuvQ6/AiuBrRixS1AZY4TK+mnf//6ueewK17NyNjmAWZSoBVCJfo0QUwXkQcYqqGkNeT7Vib0BD2eoq+n5OF9UeqJiD9p/cfRFzV4ZclKJIAVTetgLTO3EqlK/05vQO3HGpEM4vHHMK+7KhQ06SYWUvhXDQNTTeLBqT9stSz4pBQdRBbW9dsF0YsrSOoSNlZF3yBqgq3OLCYQymX3NRKCq/OxXF2OYlNYQWAHZDWzFzLgcd8ehouDjWiFRneA+S6lWKqDkotV0lJy0EWkdQMmD2YGZLMGBgJWa1HnGtmu5X6uGujar70+Em8eHYVAKB4LLcSkDv8h6HqBuZjKgBgKZHBJvv6q7rhSmXl157DxaFmWCqrVCTmEFd16IbpWBBxVS+draRIoBRI670Xd4irOsb7fQCywdFYWrNiDgG2e+UL1D8/dRp/+M0Xi37NMCn+/Psv4x9/blWTK5JVBAegqGtp1uXyBIBNIctyUHUTaTugHfJKEAgfQtXrcHGoEVYEJ0uFlgNgtaOOpnPFoXjMwQr79GIhXELVMRSU4RFJNiCt6ggpHp4x4+KZk8t45MhC0a+xTLdX7D5UiiQ691S0SCrqzKrVipu5ON3iYCVNiBAEgj6/zC2HHoeLQ42UsxwAy03izsBRjeKprAF7pkMvttBIqAYCipTrirPdSmGvBwLhlgMAaIZZ8v6I27NDFmxXkSIJ2NLvB2DFFvKZssXhjRcMAwA2hb0538drWx19Pg+3HHocLg41wuYzFKtzAKxMj2jKeuBizHIoEZAGetRyyFjB57DPg7WkBt2whs8EvdKG3L0+f2YFt93983UXl2n2dSkWl8pPZlA8AiYGLFfdmeVkwetnVlMgBPiVy8YBAMO25cDE2WtvYPp4C42eh4tDjWisfUaRCmnAetiiLstBM4rXOWRnOvSi5aAjoEjosy0HJpAhO3XV3eNnOZHBpx96tatbOhycXsMLZ1fxxPGldb0vY7swi1U959eBKJKIsYgPkkCKisP0SgrDQQVvvmgEf/LLe3DThSMAsi4oxbYc+v0yd+n1OFwcasRpnyEUtxzWUprzwKU0A6mMUTTmkLUcekscVN2AZlAEZNFxK7HrFfJa16TfLzvplP/6izP47MNHS854ME3a8cLBsod+fmxxXe9jmXHF7pF4Ol8cBIgCwZZ+X444fPzbB/CXPzyM6dUUxvt98IgC7rxuO8L2tY4xt5JjOci8zqHH4eJQIxorgpOKWw5RV8wBsHa+xRvv2ZaDauBnry4433ejw6yEgGKlra6mMk5RXMgOqPYHsm6ln71qBWRLWVh/9r1DeO9X9zf7tOuCicPj6xUH+54oNjGwmFsJACYG/DjrEoenTyzhq0+cxomFBDb3+ZzjkihAINm26V47DTbskxzB4PQmXBxqRDeLt88IF3ErAZb14JEKU1kDtuXwyJF5vOcff4FPfO9Qs065o2CLmhOQTmrOYpTvVoqrOp49vQIgt6razZnlJE4vJVpw5rWj2m6hEwsJJ2uoGhxxKPK7M7cSM2AVe+c/mScO0bSOlGbgXDSNLS5xAKwZI45byd7seD1izmxpTu/BxaFGmKmfH5D2ekSEFAmL8QyieTsvWSySrWTvkp+0/dD3PnUG39x/thmn3FEwCyCoSOj3y4imdSwnrIyboNdlOSQ0PHl8yRHjUu43zTCLFn11Eu7zW49ricUciv3u7NjO4SCA7OI+MeDHSjLrqnMHwTfni4MoOMkTzHKQRQEZw+zJ4kyOBReHGnHaZwiF1sBYnxfTq6mC7qzFYw7Wwzi9msJo2ItLxiP48hOnGn/CHQZb1Pyy6BTCHTkXB5Abc8gYJr7xzJmC9+Wj6bQrxCHslTAcUnD/izNVd5zN2AWSxRo0spjDJeMRAFlxmByw0lnPLiedWQ3sXhvPEwfFIxYEpNnHTI+4OTmFcHGoEd2wKp4JKRSHzX0+zK6lEE1pTisDoLg4sAAiAFw8HsH2oUBPBKfjdswhqEjYYovD4dkogKw47Npk7YZ/fHger9nWn/O+fDKG6bhtOhVVN+CTRfze63fgsaOL+Mkr81W9jxVcFou3xFUdsiTgsok+eD2CY4m6xYFZDf/l6q24Yfcwrtjan/M9LMshN5WVuadUjYtDr8LFoUZ0kxYUwDE29/kws5rGWkpzdsUAiqayEkKcHd1Fm8NWu+UeeCCTrpgDE4dXztnioFgxhxsvHMGTH38T/uE9+/DpX9ub8z4AuG//WdxnjxTtCreSZlUg33HNNpy3KYj/9f2Xq/Lrs5hDsaFQVkW5hHdfNYn/+PDrHbfQhCMOKSeWc8FoCF/5nauc1iQMRRKcWA4rgmP3qtqDbV04FlwcaiSjmwXxBsbmiBfLiQzmoukc/26xbCUgG5TeszlstVvu8B1wI2CLUUCWMBr2QhQITi8nIQnEWaAAYCziw017RjAx4IcsCYi7ds/3PXMW9+3PFYdOHg7EJrV5RAHvv2EnTi8lcWIxXvF9TByKBePjdq2ILAnYOhhwjkd8HkR8HpxeThQE+vORJQEstKA4lgMTh84WXE7z4OJQI7ppFhTAMZggRNM6xvu8zvFijfcAwK9kLYdqRzx2K7ph4n9860XsP2VlHwUUEZIoYDTsBaWWS6mYqw6wp+a5Fsi0bjg+ceZ66WQfuaobji9/MGjt3qsZ9MSsi2ItNBKqVWVejNGwF/NR1XErMXddPm6L1rEcnM6uG/de5JSHD/upEd2gRYPRQG42yGjYZTkU6a0EWItexOfBeJ8PPjuF0DCpE4vYSEyvpnDf/inn/8xHvqXfh+nVlJOpVAxrMFJ2sVI107lGbAFV7eZxnYj73FjxY7mhPIxszKFIhXS6tDgMBGQsJ7JZc6XEwR0LY6LABKMXXJyc4nDLoUY0g5a2HCJZQegPeJzmesUC0gCwfSiA63YNgZCsS6Xbd2yUUjx3ZqUgFdKdwSUKxFmEWLM4Fm8oRkCWclwrqm66LAdbHDp4MWPDdAA4iQrVWA7l6hwSGb2koA4EZCwnM1nLocS1zREHKT/m0LnXk9NcuDjUiOVWKr6zH4kozoSziM/jPLylxOHv/vNe/P3tlwMAfLaQVLOj7GR+fHge77jnCTx1IrePEBOHkbCC4aDiuJBYULqc5VDgVtIMx2JwxKGDRVXVjaw4VNmN1zSpq8ajeCprJcshVslysDc5hLjFgbuVep2K4kAI+UdCyDwh5KDr2J8RQqYJIS/Y/37J9bWPE0KOEUKOEEJucR2/khBywP7aZ4m9KhBCFELIN+zjTxNCtjX2V2wOmmEWjAhlKJKI4aDV7TLs9TgPb6mANCEEgu0eYdkm+XGHR4/M49M/OtKQc28F9z51GgAws5Y7XIYVW33+v1yJb/ze1c5xJg7hsm6lXHFQddMRBeZ6UXUTjxyZxz88dqIBv0VjSWtut1J1mwDW4BEoLiQsIF2M/oCM1aTmDFIqJbyyy1pgYq14uOXQ61RjOXwZwFuKHP8MpXSv/e8BACCE7AFwO4CL7PfcQwhhDuDPAbgLwC77H/uedwJYoZSeB+AzAP6qxt+lpWhlYg4AMGbHHcI+D4J2lkixVNZ8mLshnScOdz9yDP/UJcVxZ5eT+NlRqxcSq3pmMMthLOLLya5hKb+lMmoAK+YQL2E5ZFxupW8/N43/8/DRBvwmjcUdkHbEoULyARM9oHjMIa7qJS2CQTtl9exyEj6PWNINquTVNlif20VwXBx6loqrFaX0ZwCWq/x+twL4OqVUpZSeBHAMwFWEkDEAYUrpk9TKNfwqgNtc7/mK/fm3ANxISqWrdBB6iZnQDJalFPFJTiO5Ug+nG8dyyGQfypVEBs+eXnE6wXY6X3/mDAismMJS3rAeJg6sQSFjwo45lHKRAFbMgfnoKaW25UBBKc1xK6UyOmJpvaBCvd2wVFbA7VaqIA6uxTk/5qAbJtKa6aRC59Nvi8PppURJAQGyloM7hTjrVuLi0KvUE3P4r4SQl2y3Eyu5HAfgbgw0ZR8btz/PP57zHkqpDmANwGAd59USrCK40hrGgtJhn8uttB7LweXrffTVeZg0O32uk6GU4jvPz+CG3cMYDipYjheKg0fMrWUAgNGIF16PgCHbHVeMgJINSDNLgWV2sfKGtGY6u/HpIpPQ2ok7W0m2u6FWdCsZpcWBxSBKuYsGHXFIVikOhZZDp1edc5pHreLwOQA7AewFMAvgb+3jxVZLWuZ4ufcUQAi5ixCynxCyf2Gh+EzdVmEVwZW+fPu2DWDHUAD9frliQNqNTy4cDv/wYavNAtsldzIHp6OYXk3hly4Zw0CgcJJbNK0h4vMU1DJ4RAH//oFr8dvXbSv5vVlAmlLqpFhmDDOntkHVDWc3PrVSOOymnahaNiBtVcZLFS0H9849/7UxlWUhlQ5IA8BSIlPWXcdiYUpOSiuPOfQ6NYkDpXSOUmpQSk0AXwRwlf2lKQATrpduATBjH99S5HjOewghEoAISrixKKVfoJTuo5TuGx4eruXUayaW1nLmGesmLZmtBABvuXgUP/nvb4BHFCoGpN3kB6Q1w8RPX80Kod7hXTJ/eHAWokBw04UjGAzKRd1KYV/xherCsTDCZWMOEkxqWQfuLBr3oqnqpiOs0+toi90KVN10Fl0AdjV8+WwlZjnkZ2oBuTMxiuFuk1HOcmDnlGs5cLdSr1OTONgxBMbbAbBMpvsB3G5nIG2HFXj+BaV0FkCMEHK1HU94D4Dvut5zh/35OwH8hHbg9vhPv3sIv3fvs87/dcMs2Vspn9B6LIe8gPS5tTRiaR07hgP2z+24S+NAKcWDB8/hdTsG0R+Q0e+3UindRFNaQbyhWoJ2JXlc1XPqGdzT0FQ961aa6iC3km6Y0E2aE/T1y2IVbiXr7x3xeQoa78XV8llI/f6sOJQTXaWY5cB7K/U8FSukCSH/CuANAIYIIVMA/hTAGwghe2G5f04B+D0AoJQeIoTcB+BlADqAD1JK2d31fliZTz4AP7T/AcCXAPwzIeQYLIvh9kb8Yo1mdi2d48MuVwSXD7McqslW8uaJAwuqDgUVnFhIIGOY8KEzK4CPzcdxYjGB375uOwA7z75IzMG9aK0HtkNOqHpO/MW9aKr2SFags9xKzPXlvgd8HrFyQNp+X3/Ag+nVVE7lPKtfYKKZjywJCCkSYmUymtjrgFIxB2459CoVxYFS+u4ih79U5vWfBPDJIsf3A7i4yPE0gHdVOo92k9SMnOyXckVw+WwKKxAFUjYTh+Hz5Oa/swWABRf1Du4ddNZejNlsgcGAjJiq28Vf1u8VTWnY5kphXQ9MHPIb0OW01OhQtxJbZHPEoYomi0xU+nzW3z+R0R0rwAlIl6kqHwjKVYuD26ohhEAWBe5W6mF4b6UqSWV0xFXdcieJglXnUKXl8MuXbsaFY2EntbAcToW0vZiwISzMf9zJMQeWE89iKwN2c7mVhIbRiPV7rdXhVmIpmwlVz+mIm2M5dKhbiS2yiifXrVRtKmuf37pmSdVwxKGSWwmw7hsrW6mMW4nVOeRlkCmSwN1KPQxvn1ElbJfGmphphglPlY3xPKKAC0bDVb2W7SzZAscsByYOnVyUxMZZsp3ooJMtYxXCUUoRTeu1i4OSrQ1wuzvcgdp4WoduUoQUCatJreTM6VbDFtlct1LlbCUWc2Di4BZCZ2BSiToHABiwXXhVuZXyGhYqHm459DJcHKqELdardmqmvo6Yw3pgzffSjjh0seUQsGoWWFA6ruowTFpHQDrrVspJ8XS5lVZT1s86b8SaItcptQ6O5eBagH2yWFAJn48Tc7AX+bWU5rg3WSA+UCLmAGTvm2pSWfNrTxRJ5DGHHoaLQ5WwvjbswdTN0sN+6sWaBmdbKqlcy0Hr4JgDOze2Ex0IWAsSEwd27cK+2ryZ7oC0e1F1WwerSetnsBGjnRKULhZz8HvEio33WMyBCer//NZLeNtnHwMALMZVRHyesu7NrDisL+bAzpW7lXoXHnOoAsPMFl2xBS6jlx72Uy8+j+gKSGvwy6Lz4HayODDLgQXq8y0HJnS1u5WyloM7s8btVmJFd+fZ4jDTIUFpx62UV+dQbbZSn205HJ23JsclVB2zaymMRbwl3wusTxzyLQdZ4m6lXoZbDlXgzijJWg7l22fUg9eVxRJNawh5JciS9bM6uc4h33Lo83kgkGKWQ60BaUsQEqqRYzm4G9Ixy2FLvx+EAAt5qbTtophbqbo6B5atlHvNZtdSmFnNHUNbDCYOZesciqSyAlbwnItD78LFoQrcQ+2jTBwMCk8VdQu14JVEV8zBSl1kBXcdbTkYzHKwzlUQCPr92SrpUk33qkUSBXg9AhKZ3JhDIsetZP2soCJhwC9jMa4WfJ92UDwgLUI3adkkA023NgNDIcsKu3SLlSY8tZLC7FoKoxUshzecvwl3vG4rdo+ESr7G3bLbjSIJyGwQt5KqG869wakOLg5V4Db911Ka1QXUrD5bab1YgUprwYilrRx1Ft/QOthyyA9IA7mFcNE6xQGwFv24XTvBcGfwsNkFPlnEUFDBUqeIA4s55LmVgPJtu1VbcDf3efFPv/UafPrX9gIAji8ksJLUsLmCOAyHFHzi1ovLVueXtBw2kFvpc48ex613/7zdp9FV8JhDFSTz3BasC2i1dQ7rxefJdSv1+2Vnwe3kzqyaYUISsoOLAEscTi4mMB9NOzUbtbqVgOzAn3SRVFZFEhzrxOcRMRiUsdjRbqXsHOl8wfyz+w9hIabiyq1Ww2NZFPDGCzY5FdLPnVkBYM3FqBdZzJ0bzVAkEUsdcv3qxbK00pVfyHHglkMVuDNK1lKak07arGwlb05AWkfYlZHS0W6lIkH6124fwJG5GF73qZ/gW89OQSDl8/IrEZAlxNN60cZ7Ia/ktO5mlkMnu5X8JUaFUkrxwIFZHJpZc/7e7LqKAsFo2ItnT9ni0FfecqgGFohmlgzDqnPYGG6lhKojo5sFM805peHiUAX5biXnga2y8d56cee/R1NWQNrTBW4lzaAF7ouPvvl8/OQPbsA1OwfxyrkYwj5PjmWxXsI+CdG0lmM5sFRWd3sSPxOHWKeIQ/H2GUBhK+6zyynMx1TEVaNAHABgvM+Hc1FrF9wIy2HncBB//LYL8aYLNuUc30huJXaPbJTfpxVwcagCtrMLKZJlOdgLdLW9ldaLVxJyKqQtcbDdSh0sDmqJ9N4dw0F88T37cO15g9g+VFtfJUbE58FaSsuNOTBxcKVr+jwihkIyEhmjYkZQK8jGHFxFcCVGwj5zyupYn8zoTtW5+17b7LIWKqWyVoMgEPzu9TsKCuUUaeNkKzFxqNTLipOFxxyqgO3sxvq8luVg+/2bFnOwLYe0ZiBjmHa2ErMcOvdh1QyzwG/N8HpE3Hvna+u2fBxx0EwEZBGJjOFqQOcSB9tyAKyqEL1kAAAgAElEQVRisYkBf10/t17Ku5VyF6z9p5ed46puQBaFnOFILH11ICAXBJEbiSIJG2YSHNtAVKpI52ThlkMVOOIQ8dlupeZaDiwg7QRwXZZDJ4uDFXMofU0IIVXNtCgHE4e0bjqWAstWYt1JRcHqKDpsi8NCB8QdVN2EQJBTG1PKrfSMHU8AgLWkVnBNx/stcWiE1VCOjdRbiW0guDhUDxeHKmBupbGI13Yr2ZZDk2IOXo+Vyhp1FY1lxaFz3UqaYda9+Fci4vMgrZmIpTWnsIv1VmJVwD6PCEIIBu2usJ0Qd2Dzo90WgJOt5JoGt5zI4Nh8HBMDlgCsJDMF9TTMcmhEvKEczK3UgbO31g3rUZbmvaKqhotDFbCd3WjEi2TGcPyWTSuCs10FCzErjdBd59DpqazNainCYCmf81EVfkWCQKziO1Egzk6cfcy6ldqfjqlqRkFLbBZzcFsOL02tAgCu32WNwV1JagXXdLyvRZaDfX9nOtharQZKqVNFn94g2VetgItDFSQz1mB41hmT5X43rQjOXkTmY1ZGSsjbHZaDqjffcmA1EvMxFYokONfFIxKn5TRbdB3LoY1upR+/PIcb//ZRxNJ6QTzGKYJzicOrczEAwD67vmE1mSmYPT7e54PXI2DncH3B/UpkR4V2tziougnDTmFNd0ByQrfAA9JVkMzoCCiS01OfLTbNDEgDwILtDgl7Pa5U1s59UFtpOSwlVOzZHHaaw3lEwdmZs0CvIokIe6W2isMr56I4vpCAIokFXU/9RcThlXMxjIQVx3W0ktScnlKMgCLhRx++oWLrjHrJGRXa3B/VVNxde7nlUD3ccqiCZMaAzyM6u1bmpmhmERwAzEWZ5SA58Y1OHhOa0UtnKzUKJg6UWosX21XLolC0DcRQqL2FcMwFeXQ+VnBtPKIAj0iQ1HIth/NHw07m1WoyU1RwJwf9TbfSmJh1eyEcm3sB8JjDeuDiUAVJ1UBAEbO7VnuxaVoRnL24zTPLwefpmiK4VlkOgCUCbIH0iIKzmPldO+2hgNLWmANbjDSDFsQcgNxqeMOkODoXx/kjQac9eSuuaSnY+Xa7WynHcuDZSlXDxaEKkpoBnyw5bZOZu6dpRXBMHKIqBGK1qiaEQBJIR7uVKqWyNgK3OOTEHCTi7Mx9OZZDezuzuhejfLcSwOZIW4vX6aUEVN3E+aPhHFdSsxIfKuEEpLtcHNxde3kRXPVwcaiCpKojIIsYtIfXzMVaE3OYj6URVCQn/VESSUePCbVSWZtXlAXkNu3zeoRcy6FIj6DBgNLW5nGpHHEovF/8soSUbV0cOWcFo88fCTmWAwDITRbcUmTdSl0uDhnuVqoFLg5VkMwY8Msiwj4JkkAwb8cCmlkEBwAnFxM5QUePKHS05aC2wHLwiIKzq1YkMS/mUOhWCnqliqM4m4l7BnMxcfB5RGdeyJG5GAgBdo0E4ZdFsJKItrmVnIB0d++2Y2nuVqoFLg5VkLLdSoQQ9AdkJ1DcrCK40YgXXo+A63cN47Pvvtw53uniUK59RiNhriVFEhyXi8cVkHa7lQKyCM0oP1CnmaQquJUCiujsbI+ci2HbYABeu4gvYBfJNTvwXIqNEnNg1dFA9wtdK+GprFWQsN1KADAYkPGKbf6z0Z2NZiio4PD/ektONS1gtV7o5MZ7rUhlBSzX0sxaGl6PCMVV5+CIg+zusWR9nszokCW56eeWT1ozQIidXVUkIB1QJGeM6lw07RS4WV8TEVf1NloOG8StlJPK2t2/SyvhlkMVpDKG48dmM3mB5lkOAAqEAWCWQ+eKQ0Y3Cwq2mkGu5WBdJ48oOIH8fMsBKOxf1CrSmoEddifaYlZVQJacbJqEauQ0D2Rxh1Zc02Jki+C6e7fNrm/E5+mIDr3dAheHClil97pj4ueIQ4sDhR6xs7KVTizEczKBtCbO1XbDxMHrccUcpKxbyR1z8CtZy6EdpDQT24cCCHklpzGgm4AiOr2h4qqeE4hmQtHsOE4pHMuhy4O4cVWHXxbhd81J4VSGu5UqoOomTJrNgBl0iUOrzX1JFDqmt1IsreHt9zyB12wbwD/csQ+UUmSMNlgOojtbyfobeeVCy8Htd24lqh2v+trvvrZoRTMbewpYi1hQcZ87Ewcec6iHhC26Xo/I3UrrgItDBZgZynaj/TlupVZbDgIyeme4lb729BmspTT89NV5rKU0x5XTiuBpjuUgFcYc/C63Eut8mmiT5ZDWDHglAZdu6Sv69YAsIZHRLQs1z3II2ELR7jqHjeBWCjJx4JZD1XC3UgXYosJ2cTmWQ4sfWo9I2mI53P6FJ/G9F2ec/6c1A196/CQmBnzQDIqHXp5zOne22nKQXZbDxIAfV0z24bKJiPNatsAm22Q5WJlupWs/AooEkwLRlA7dpHni0N6Yg08WIZD2Ni5sBAlHHAQuDuuAi0MFmOWQDUgrztea1T6jFB5RaHm2Uloz8NSJZTxyZN459qOX57AQU/GXb78U430+/OClGWg6m3XcfGsqYjdAVFxFcLIoIKhI+PYHrsV5m0LOa51pa21aFNKaWXZaGxOvuVi2j1b2a+2POVw20YfHjy215ec3CiuWI8IrccthPXBxqEAiz63UzoC0JJCW99ZfswcOnVxMOMdYEeClExH88qVjeOzoorO7bHaFNOByK0m5vZWKwdxKSbX1biVKKVK2W6kUzCJltTMBuVhAun2P6Q27h/HS1CpWEu2fiVErcTsLzLIceMyhWrg4VIAthMMhy2JgMwKA9sQcWt2VNVpEHFjFaUCWcMXWfugmxfGFuH2Ozb8mm0JWYDfi9+T0VipGwIk5tH7HyAK53rJuJdaBV7X/77Ic2hyQBixxoBR47Nhi286hXhI85lATXBwqMLWSApCdvsUsB49IitYiNBNPG3orMcthNak5u0eWGigKxHGDLCes17UiIH31jgH82/uvwUWbIxUtB2dOcxssB5YC6i1jTTExYIOdgkUC0u2qkAaAS7f0oc/vwU+PLLTtHOqFBfp9HpHPc1gHXBwqML2agtcjOKLQ75dBSHML4EohiULL20AwcQCAE7b1wHZiABCyc/eXE7ZbqQW7XEIIrrQnpTFRKPVzZTto3Y6YA2udUT7mYIuDYzm4+kK1OSANAKJAcN15Q/jZ0e4Vh5h9vyr2bHZOdXBxqMD0SgrjfT7HShAFgj6fp+XxBqC9lgOQdS3FVB1B22II5lkOrXaBKBUsB8CyHtphOTAXhk+uPubQSQFpxsXjESzE1Jw2FN2CZpjI6GY25sArpKuGi0MFpldT2NLvzzk2EJDb4gduR8whVxysuEI8rSNkL1xsd+tYDm1I77U+lluAxbbEHJgLo7xbKXfqX9EK6Ta6lQBgOGjF29gck26CCVq2CI6LQ7VUvOsIIf9ICJknhBx0HRsghDxECDlqf+x3fe3jhJBjhJAjhJBbXMevJIQcsL/2WWJvxQkhCiHkG/bxpwkh2xr7K9bH1EoS4/2+nGODAaXlwWjAcmW1urcSE4fJAb9jOcRdloMTc0i2x3KQKwSkAauFRjvaZ7A06LJuJZnFHAoD0ixDrp0BaSCbjLHQhfUOS3acrD/ggVeyOvQaHTwTpZOo5q77MoC35B37GICHKaW7ADxs/x+EkD0Abgdwkf2eewgh7Mn4HIC7AOyy/7HveSeAFUrpeQA+A+Cvav1lGk0yo2MlqeV0ygTaaTm0vrfSWkpDUJGwa1MQJxZscUhnYw5WCwvSPstBKh9zAGzLoQ1FcMy/va6Yg1zoVmpnzAFwiUMXWg4zq1ZCyVjEB6/dDuTlmSjufep0O0+rK6h411FKfwZgOe/wrQC+Yn/+FQC3uY5/nVKqUkpPAjgG4CpCyBiAMKX0SUopBfDVvPew7/UtADeSVqcBlWDazlTakmc5/MbVk3jfG3a2/Hw8otDymEM0pSPi82D7UACnlhIwTWq3I7AC0YQQBBUJy/a0tVYvZO4K6VL4Zakt3TjTTkC69LnJtrhmDBM+j5UBxtg5HMQtF404wfd20c3iMLtques2R3xO5toXHjuBP/7Owa5vC9Jsan2SRyilswBgf9xkHx8HcNb1uin72Lj9ef7xnPdQSnUAawAGazyvqqg213lqtbg4XL9rGL959daGn1clJJE4lcitYi2lIezzYEu/D2nNxEoyg1haywmcBr2SY7632nKolMoK5A7UaSXpKrKVgKyF4HYpAVYg/f/95j5MDPiLva1l9PtliALpSnGYWbOe4ZGI4sR+jpyLAgBWElrJ93EaH5AutuOnZY6Xe0/hNyfkLkLIfkLI/oWF2lLr7nn0GC74kwer2jVMOzUO7X04GR5RgNbi3krRlIaIT8KAHZRcTmSQyOTOHQgqHqfgq9WZNVnLofTP9clSW+Y5sOCnr5I4yLnxm05DFAgGA3JH9VhajKtYS1Ze3GdX0xgKKlAk0ekyy9yjy11c9d0KahWHOdtVBPsja7wzBWDC9botAGbs41uKHM95DyFEAhBBoRsLAEAp/QKldB+ldN/w8HBNJx72Wu6Q1SpurKmVFDwiwaaQUvG1rcCKObQ+IB32epyGg1OrKRh5DeJCrs/bZTmU+7lWzKEdAenKMQcgm7HkrnHoNIZDSkdZDh+49zn8yXcPVnzdzFoKm/usinr2d2CuWS4O5an1Sb4fwB3253cA+K7r+O12BtJ2WIHnX9iupxgh5Go7nvCevPew7/VOAD+x4xJNod9vLXIryco3xvRqCmMRH4Q2ZCYVQxIEGCZFEy9PAWspDRGfx7luZ5aSALL1DfmftzrmUKkIDujsmAPgcivJnWk5ALY4dJDlcHo5gXN2+m85ZtfSGIvkigNjuYo1oJepJpX1XwE8CeB8QsgUIeROAJ8CcDMh5CiAm+3/g1J6CMB9AF4G8CCAD1JK2VP5fgD/ACtIfRzAD+3jXwIwSAg5BuCjsDOfmkWfv3rL4eRiHJNt9ve6Ya6TVloPTBxYT6nTtji4rQW3O6TlqazriDm0UlSB6iqkgawoBJXOFYehYOdYDqZJsRTPVExPppRi1t7gAYXuveUOErtOpOLdSCl9d4kv3Vji9Z8E8Mkix/cDuLjI8TSAd1U6j0aRFYfyu4ZkRsfh2Rjed8OOVpxWVbAFUDPMlrhvMrqJlGbkWg7LtuVQpFgLaL1biWXSbAqXdv35ZWtmgqqXb5/daFTNACHFZ0e7ybqVOlcchkMKFuMqTJO23ZKOpjXoJq2YnhxN60hkDJdbKffvsFzFBrGX6bkK6axbqfyN8dLUGgyT4orJ9qYRupFscWjVTIdo2rpGEb8HsiQg5JVwZtkK5pVyK7Xactg5HMQv/ujGkpPWgOzi+9SJJfz9w0dbdWpI6ya8klixQaNjOXRoQBqwqqQ1g+ZUzLcLFhiPV4gjza5laxyArAWnSAIiPo9Tm8MpTufejU2iUszh6784gz6/x2kyd3kHiYPjVmpRxhJbCNj8hMGAXNRycLuY2tEHaFO4cDazG+ZO+PxPj+OpE8t43xt2tkTEUhmjYrwByFoMnexWcldJu0fltoNFu6amVJLBfCyNX/7s43jLxaMAkLUc7FTWrYN+6CblqawV6Ny7sUn4ZBGKJBSNOaylNPzp/Yfg9YjYMxbGjqFAznCfduN2K7UCJg4sw2sgIOMUizl4C91Ksii0vI15NbDF97kzqwCseRSt+LumNaMqN5afuZU6PCANAIsxFbtHQhVe3VyWbHFIZoyibq6nTyxjPqbiq09aVdBZy8F6frYOBrCSyGCJWw5l6dy7sYn0++WcyVZ/+M0XsW0ogIjPytdXdRNPnljCr16xpcx3aT2sn1Or3EqOOPiy4sDIiTnY4tHOuQPlYD2KWLvzWFprjTjoZsUaBwAIyqwIrrNTWYHO6K/krrdIZHSE7PuP8by9CQAAgcBJRWdDl7YO+EGQTa7gFKcnxaHP73FiDpRS3P/iDDKGidGwFxeMhjAcUvDY0UVcsbW0H7sdsMW3VZZDNM+tlCMO3sJspXa3li5FfqA3mmpNzUMqY0CpQhy6yq3UARlLS25xUI0CcXjuzApes60fAiGYi6adWF1QlvCOy8fx1kvG8M39Z/H82VVwStO5d2MT6ffLTrbSfEx1qntn19L43et3YN/WfpxcTOD1u2ortGsWbMBQK1JZF+Oq0wwuKw7WAiGLAhRXG2oWc+h0y4HBAu3NRtWrjTlY59fJAemQPUnt3Frl2oJmsxDPWv35bVFU3cDLM1H89rXb8IE3nJfztxYEgk//570AgIcPz2ElkQGltCNdoZ1A596NTaQ/4MGRczEA2dTMD9+0C784uYx3XD6O/oCMx//nm9p5ikWRnDqH5loOac3ADX/9iDMDwR2QBgrdH0HHcuhMccj35ccaJA7fenYKB6fX8Ge/clHO8bPLSTx3ZgWpjFGVW6lUb6VOghCCLf0+nF1pvysm13LIFYdDM1FkDBOXT/Yh4vcg4vfkvx2AZQXrJkU0rTv3NyeXznyam0yfX3YC0szveOvecfzLe69ueyZGOVgVcLM7s0ZTGhIZAxePh3H7ayYci4C5lfJ3uMEusRx2jwQBNM6t9MiRefzL02cKxPqfnzqND339BZxeTlYVkGaLU6cvUpMDfpxZTrX7NLAYV51nIT+dlcUbKmUZsnuZt9AoTWc+zU2m3+/BakoDpRRnlpMQCApmNnQirbIcYvYD997rd+BTv3qpc3zArpJm7boZTCzaPXegFEGvBFEguN52EzbKrZTOGMgYJk7Zac8M5npZiKlVuZWu2TmEu3/9Clw+0VkxrnwmBvw4u5xseaV5PkuJDLYMWM9rfiHcs6eXMRbxYqRCenM/F4eKdObT3GT6fDIM26Q8s5TAWMTXsbteN9mYQ3PFIZ62xCG/SyhzK4Xy3B8hWyw61a3klyX863uvxodu2gVCrMrZRsDaYxyZi+Ucn3P1/KnGchAFgrddOtbxvu+JAT/iql5V65lmshhTnbY2brdSWjPw0yMLuGF35VjhIBeHinTm09xk3C00Ti8nsXWwc/onlUOWWtNbiZnq+RYCKyDMdyt5PQIkgXS0wF61fQBhrwdBWWpYzIG1AWfxK8a8K6Onle06ms2EPdeExenaQSpjIJExsG0wACDXrfT40UUkMoZT/FYOpxiWi0NJOvdpbiLuFhpnl5Md1VyvHMxy0JvtVkozccizHBy3Uu5xQgiCXqljU1ndhH2ehsUcWNdVtzhQSnFuLe1sOLzSxhGHSft3amdQmtU4FLMcHjx0DiGvhGt2DlX8PuxeXuLiUJLeFIeAtSOeXklhMZ5xbvpOR2pRV1a2G8t3K/llK52xWMplUJEgd8FCGPJKDYs5OJaDy60UU3WkNAO37R1HQBYxFOrcBIf1MtFvPSfttBzYYp4vDpph4qGX53DzhSNVWbA+j4iwV+qI7KtOpXNz55pIn205vDRlZTZ0i+Ugt6h9RtxePIsVZf35bRfjwrHC9gmDQaVjJ5m5CXs9DXMrsZjDmeUkkhkdflnCvB1v2DEcwI8+eoPj294IBBQJgwEZZ9soDou2y24opFhDnGyBPjQTxVpKw40XjlT1fQghuGRLxFkDOIV0/tPcBJhb6dnTKwCArQOBdp5O1ThdWZvceI9ZDsXy7t95ZfGWIp/5tcuqqgZuN2GfhJnVxhRypTIGxvt8mF5N4ehcHJdN9GHOLhzcFPJ2RQbcerEyltqXzsoG9AwGZAQUybEcWOX2xED11/zSLX344s9OVN0Dq9foSbdSxOcBIcD+0ysYDinYualLxEFojVsppupQJGFdAeYdw8GuWAxDXg9iav2WA6UUKc3A3kkr/fTJE0sAsplKI2XmS3QzEwP+trqVWHbRQEBGUJGcjQxrv72enlmXbemDblIcno02/kQ3AD0pDqJAMGb3Ufr2+6+Bv4O7YbppVW+leFrvChdRLYS9UkMC0hnDhGFSXDgawnXnDeHunxzDfDSdtRwq5Nl3K5MDlqXUzKQIs8wo3JVEBrIkwC+LOZYDi0UMBqoX5csmIgCAF3mPpaL0pDgAwP2/fx2++1+vxUSXxBuA1nVljat6RzeBq4eQHXOot5ArnbEWR58s4c9vuxiqYeIvfnAYc9E0goq0Ya/fYECBYVIno60Z3PA3j+Dep04X/dpyIoMBvwxCiDX+1S6CW45n4POI8MnVu4dGw14MhxS8NLXWkPPeaPSsOAwFlZzmcd2A1LKAtN7RTeDqIeyzRoayQGatJDVrcfR5RGwfCuDO67bjey/N4Pmzq2VHlnY7rBVJUqvv+pXCMCnOLqdwOK92hLGSzDjVzblupcy627ATQnDZlj68yIPSRelZcehGArKIoaCCf3tuGqrenIcTsGIOG3XnywYXRescd5myxYUtlr9+1SQotVwUoxvUpQTA2Zmn6hTXUrD7eqnE3AhLBKy/oV+WnK6sS4mMU7uwHi4Zj+D4QqJpv083w8Whi5BEAZ96xyU4PBvFpx96tWk/J57WC6qjNwqs93+9bhGWxsqyXCYG/Lh6xwAAVOzr082wLrPNWkzTmmUVl2prsZLUnNbx7phDLZYDkC2Ga1R680aCi0OXcdOeEdy2dzP+6fFTTWuAFlc3cEDaZ/1e9RbC5VsOAPCuKycAYIO7lazrl8w0J+bAqs5LVS5bMQdL4IPumEON4sDaz9frZtyIcHHoQvZsDiNjmA25odnQIzcbPSAN1L9TZJaDOwD61ktGsWcsjNdsHajre3cyjlupSTEHNnhrKV54X+qGibWU5sQcAoqElGbAMCmWEmpNBYdM7PLnQnC4OHQlzOUTr9M1MruWwr6/+DF+9upCzvENHZC2f69601lZ6wz3MB+/LOGBD12Pm/ZUV6XbjTTfrWR937WUVpB4sWrHiQYCuT2+luIq0prpuJvWA/seXBwK4eLQhTCXT7273+mVFHST4uBMNpVP1a0ZBRvVcgj7GmM5pItYDr2Ak63UZHEArMwkN6yDKutwwCr4T9tFebVZDs39fboZLg5dCNvVx+rc7bC+/O5eOaVmOWwU2O9V70yHZJGYQy/gb5FbCbBcS5979DgeOTIPILc6GsiKwxl7mmNtMQfbcmhSDKWb2ZgrwAYn7FgO9d3Qa7aZ7m6HkJ3lsDFvDUUSIUtCwwLS1cyI3kg0O5XVbTnMRdP42x8dgSgQfP2uqx1LwpkrYgeTTy1Zk/gGakhlDXC3Ukm45dCFNCrmsFpEHErNcthIBGQRSbW+xa1YQLoXYGLYPLdS1nJ44ewqdJNCNyne+9VnMbViNfxjFsLuEas78EMvzwGoza0UsP9++eNGOVwcuhLmGonX2UBuzd6JzaymneCfYzlsULcSkJsfXyupjAGBdO7c7GYhiQJkUajLrXTkXAyP2q6ifNzFnc+cWgYAfOSmXViMq/jeS7MAsvNYtvT7sXskiFfsaupa3ErNTs3tZnrrzt4gBBvsVjJMilm7jbUTc9igRXCAZRXV62NOZgz4Zanj5z43A58sIlXH9bv7kWP47998sejXVJfl8PwZq63F7VdNwi+LePHsKoKKlNP25o0XbAJgiXQt1q4sWWIX55ZDAVwcupCg3BhxWHW1kDi9bPlte8Fy8Mti3W6EVA/PAPB5xLrcSsuJDBbjmZz4AiNtWw6iQJDMGBgKKhgKKnjdjkEAWauBceMFVtrwQECuWaj9isgthyJwcehCBIEgqEj1i0NSw1DQyg1ncYfYBg9IA7Zbqc7FIK0Z8Mm9+fj4ZbGuxnss6+jcWuHQJWY5sP5UO4atWSuv3z0MABjw57qOrpjsQ8TnqamvEiMgSzzmUITevLs3AFZHyvpiDqspDReMhiCLgiMOGz2VFWCLQb1uJR1+z8a9RuXwySLSdVgOK06sq3CiHLMmNvdZ4rBzOAggKw79eXEFSRTwwTfuxG17x2s+H6v1N7cc8unNu3sDEPLWbzlEUxomB/zY0u9zah3iqgZJIFDWMQWu27AC0vW6lUx4eyxTiVGvW8kRh7U0pldTOL2UwDU7hwBYbiVRINgUYuJgWQ7bBv24YDSE7UOFUxvvev3Oms8FyO3uysnCxaFLCXqzvexrZTWZQZ/Pg4kBP07bhURL8Yw9RnXjBloDilj3YpDK6PD3asxBFmsuIkxlDCdddWY1hb976FX84MAsDn3iFhBCoGomvJLguImY5UAIwbc/cA08TcgO45ZDcTbu9nCDE/J66qryNU2KtZSGiM+DXZuCODofR0Y3cWB6DReMhRp4pp1HQ1JZNaPnahwY/jrcSsuulhizaykcmokimTEwH7PmN6R1K9DP0lJZzMH6uVJzxEGWePuMInBx6FJCioR4HVW+8YwOkwJ9fg+u2NqPjG7i+TMrOHIuhku39DXwTDuPgCxCMygyeu0T9ZKZXhYHyZmEt15WXK24Ty8lcWw+DiCbEJHWTCiSgDfvGcUdr9uKif7mj/FtRILCRqQucSCEnCKEHCCEvEAI2W8fGyCEPEQIOWp/7He9/uOEkGOEkCOEkFtcx6+0v88xQshnyUb2aTSIemMOa3ZfpYjPgysmrT/R154+A92kuGxLpCHn2Kk0omVCOmP0XOsMhtcj1tw+g8UbBgIy9p9eQcYuvmRuTVU34fWI2LM5jE/cejEEoflLQSNSmzcijbAc3kgp3Usp3Wf//2MAHqaU7gLwsP1/EEL2ALgdwEUA3gLgHkIIe7o+B+AuALvsf29pwHltaNzzc2th1SUOoxEvNke8eOCAVYG68S2H+putJTWj55ruMfxy7eLA0lgv2hzOsdyyloMBpcWiG2yAm3Ej0gy30q0AvmJ//hUAt7mOf51SqlJKTwI4BuAqQsgYgDCl9ElqjTb7qus9nBKEvB4kMwZ0ozbXCKuO7rPzxq/Y2g/dpBgKKhiLbNwxl4Dbcqh9t5jqYcuB1TnUMomQuZUuHresU0USMBbxOtlyac2A19Nab7dflqDqZs3P0kal3r8CBfAjQsizhJC77GMjlNJZALA/buL+GmMAABjDSURBVLKPjwM463rvlH1s3P48/zinDKyCudYFbjVlPaR99shF5lq6bEtkQ2cqAe7RkLXtFg2TQtXNno05eD0iKM1tr10ty0kNhAAXjFpJDyw99bTdWVXVzZanUfNRocWpN5X1WkrpDCFkE4CHCCGvlHltsRWHljle+A0sAboLACYnJ9d7rhuK7FwCDRF/dX2QPv/T4wCA992wM8etBFiWAwBcssHjDUD9MQdn0E8PWw6AZT2tt4XIatJKlZ4YsALNF46FAQA/Pmw14lM1o6DQrdmw+yGZ0Z3ngVOn5UApnbE/zgP4dwBXAZizXUWwP7L2i1MAJlxv3wJgxj6+pcjxYj/vC5TSfZTSfcPDw/WcetcTUlhn1uoXuO+9OOPEFZhbiT0Ml4xH8KEbd+HX9k2UfP9GIVDn3GDWkbSXYw4AamqhsZzIYMAvY+uAH5JAcPlkHyYH/ViMq0ioOtKaCa/U2uvq5227i1KzOBBCAoSQEPscwJsBHARwP4A77JfdAeC79uf3A7idEKIQQrbDCjz/wnY9xQghV9tZSu9xvYdTgpCXjbusfoFbTmQci2EtpcHrEZydnygQfOTm3djc52v8yXYYjhuhxsWABWN7tfGe15kjvX5xXUlm0B+QMRhU8OOP3oB3XjmBSduKOLuSRFo3oLQ45lDvZmGjUo9baQTAv9v+aQnAv1BKHySEPAPgPkLInQDOAHgXAFBKDxFC7gPwMgAdwAcppezpfD+ALwPwAfih/Y9ThuA6ZzpQSrEUzzjBPmbe9yL1jobshYFI5WAzEFKZGmIOCQ3jdt+kbXYrDCYOZ5aSdoV0a0WXjwotTs13N6X0BIDLihxfAnBjifd8EsAnixzfD+DiWs+lFwmtc6ZDXNWRMUxkDBOGSbGa1HpXHOT6gvnzMaub6Kaw0rBz6iYct1INi+lqMoOLN4dzjjnisJy0K6TbFJDmbqUceIV0l8JiDtWKw7KrMjWa0rCcyGAw0JuLm9cjQCC1uxHmo1arB9Ycrtfw1RhzoJRaMYe8gHPE54EsCViIqW2pc+DT4IrTm3bxBmC9MYfFeFYc1mxxuDBvB9crEEKstt01LgbMchgO9aa4siyt9fRXOjwbxeHZKFTdLMhGIoRgOKhgIa5aFdItTmUNNqDuZSPCxaFL8XoEyJLgtCOohNtyWE1pWIyrGGpxymAnUU/zvfmYiojP07MB6axbqfrF9JM/OIzHjy0CAPqLpF4PBWXMrqZBKVpvOThuJW45uOHi0KUQQrB9MIATC/GqXr+cUJ3PF2IqomkdAz3qVgJYm+YaYw5RFZt61GoAanMrHZuPY+9EH7YO+nH9rsI09KGggiNzMQBoeREca73e6QHpZ04t4/h8HLdf1ZoaLx5z6GJ2bgrg+IJVWfrjl+dwcjFR8rVLLsvh5KIlKPWMVux2aunEyap452Ppng1GA+t3K8VVHeeiady8ZwT/5/bLi6ZLDwUVzNpjQ1ttkUmiAK9H6Pi23V/++Sn8yXcPtszC4eLQxewcDuLMchKxtIYPfO05fOJ7h0q+dskVczhhC8pgL7uV7FGh//s/XsGnfvgKzthdQUvx8kwUN/zvR/HEsUXMRVWM9GgwGnAHcKtbTE/a99vO4cIpboyhkAzDtBojtMNdF1Q8iKbqG7vbbBbjKjSD4ue2e67ZcHHoYnYOB2GYFA8cmEXGMPHY0UUsxtWir11OZJyh7Y44BHt39xtQRJxcTOLuR47j8z89jps+/VOnjXkxXrVdHs+eXsFCTMVwD1sOokAgS0LVMx1O2JbqDnuqWzGGXPdiO0bUDgbkHOu6E2Hn98iRhZb8PC4OXQwbofitZ62+hYZJ8f0Xi3YewVIig5GwgpAiOQ9rr7uVmJC+9/rtyBgmFuLpkq9nLaWfPrmMjGH2bBorI6RIOLdW+nq5OT4fh0CArYOlB/e4xaEdlsNgUM5J2uhEluz79dEj8zV1xF0vXBy6GDZC8ZlTK9jS78MFoyF854US4hBXMRCQEfZ5nLTWXnYrMdfIaNjrDLcvlxbMWkr/4uQyAPR0QBoA3nrJKB44MFuVQBxfSGBywA+lTOVzrji0wXIIKs7i24nohomVpIbxPh9m19JO8L6ZcHHoYgKK5MxeuHyyH7+ydzNeOLuKhVjhTW4VHylOi25JIAh7e7NCGgCCdvriGy8YrqranFkObHLZSLi3LYffe/1OmBT44mMnKr72+EK8rEsJAIZD2Y1KORFpFoMBOScu12mw2dtvv9yaZvCzV5vvWuLi0OUw19IVk304f8TqkT+9msp5DaUUS4kMhoKyIw79AbklIxg7FWY53LB7U1UFhVMrqZx2I71uOUwM+HHrZZvxL0+fKdsZ2DQpTi4mygajgQ6wHAIyYqoOVe/MjCUmXHs2h/GD/3Yd7rxuR9N/JheHLue8TZY4XD7Zj1Hbiji3lisOiYyBjG5iICCjz2ft0HrZpQRY120oqODa8wZdlkPxgHRGNzG7lsKNF2xyjvVyKivjlotHkdIMJxspH80w8f0Ds1B1s6LlEPF5INmblfbEHKy/Z6fGHZZcruCLNkcgtmBjx4vgupw3XrAJL89EsWcs7OzgZvP8wMyXymIOQG8HowHgP122GW+7ZCzHeiplOcyspmBS4Oqdg3jo5TlQZC2PXmaiP9swr9iQqI/e9yK+9+IMhoIKrjtvqOz3IoRgMChjLqq2J1vJfh6W4hmMRTqvbf2SXcTaygxDfod3OTfsHsYNu62K037RamCWHyRkKXBDwWzMoVeb7rlhwhCQJRBS2nI4u2LFGyYH/LhwLIzFROcGLlvJxIC1iLLrk8+x+Thet2MQX/mdqyBXseAPBRXMRdX2WA62Jd2p6awsiWSohZs6Lg4bCEIIxiLeAsthasVyM20KK+jjlkMBgkAQlCXESvjOWTB6csCPj//SBeuavreRCXk96PN7nEyufFYSGVwyHq5KGIBs3KHV8xyA7I68UzOWluJqy5NIuDhsMEbD3gLL4YUzq/B6BOweCeHg9BoAHnPIJ+SVSrqVzi6n4BEJRsLenpiUtx4m+v04u5IqOM7ac69nHjQTh1ZPggOym6VWxBwopVb32XVYSEtxq9V5K5NIeEB6gzEW8WI2mvuwPn92BZeMR+ARBSfjppero4sR8npKu5WWk9jS729JELDbmBzwF7UcEhkDGcPEgL96cRgJK5AEAlls/bIUUiR4RJLT2r5ZPHx4Hlf8+UNYWYcQLSXUlj+zXBw2GKMRH+bWVJh2nxpVN3BoJorLJ/sBAH32w5o/cKXXKWU5PHpkHv9x6Bwun+hrw1l1PlsGfJheSTn3G4MtfOuxHH7rmm344h372pJiTQjBYKA1hXAvTa8hmTFwokyjzHwW45mWxhsALg4bjrGIFxnDxLloGg8ePIeD01FkdNNZ3K6Y7Mcf3Lwbry/SNrmXKSYOpxYTeP+9z2H3SAifuPWiNp1ZZzPR70fGMDEXy3VlMvfMetyXm8JevPH8TZVf2CRa1UJjyg7g59cjlWMpobbcFczFYYPBah0++cBhvO/eZ/GH33oRALB30hIHWRLw+zfucnrycyyCRdxK335uCqpu4Eu/tc8plOPkMmHPfz67nLvQsYre9VgO7WYwqGCxBeIwbcdoplaSWE5k8NcPvgLNrrxnmCbFvz8/hYxuHV+KZ7hbiVMfrJ3GAwdmIRCrA+to2NuRududRDHL4YGD53DV9gF+7cow0W+ns+bFHZhbaT0xh3ZjtdBovluJZQ9Or6TwgwOzuOfR43hpai3nNT99dQEf+caL+Mkrc0hmdCQzRsszDLk4bDCY5UAp8PG3XogLx8J4/e7yBUgcWxxcKapH52I4Nh/HL10y1saz6nzG+30gJJvuy1iuIebQbgYDzXcr6bbLF7DcSkftBnqzeV0NnjqxBAA4uZjMqY5uJTyVdYMxFLAyPggB3nnlFvz2tdt4lk0VhL0eZHQTqm5AkUT88OA5AMAtF422+cw6G0USMRr2FlgOy4kMRIEg7O2eJWYwqCCZMZDM6E2rgD8XTTtDjaZXUkjbo1ZnVouLw5nlBE7bg6hYRXqr6J6/HKcqBIFgx3AAuzaFumrX1m7cnVmVoIifvDKPKyb7er77ajVsGwwUZN6sJDPo98sgpHs2JsN2M8X5qIptQ81ZGplLafdIEFMrKcdSmVnNBvRjaQ0H7Hqk00tJHJu3rAvWR61VcLfSBuRrv3s1/vqdl7b7NLqK/Lbdp5YS2LM53M5T6hqsWebxnAE0Vov47gris2FEp5aqTzFdLywY/drtg0hmDKddh9ty2H9qBSa14oenl5I4Oh9HyCs54tUquDhsQIZDCgIKNwrXQ1Bhbbs1xNIaVpMatrTYjO9Wdg4HEUvrOQVkKwmt62ppHHFYR/3BemGWw2u2DzjHvB4hp+XNUyeXIIsCfmXvZsyspXB4Nopdm4Itt8K4OHA4yLUcWP75OG+VURWsHffxhbhzbDmZ6TpxGA4qCMgiTi0V7xXVCKZXk9gUUrBjKDvf4pqdQzkB6adOLOOyiQjOHwmBUuCFs6stdykBXBw4HADImekwZefsb+nn4lANbJCPWxxWElbMoZsghGDbUKApbqWluIr7njmLY/NxjPf7nI1HyCth70QfFuMZpDUDsbSGg9NruHrHILYOWtfVpMCuTaGGn1MluO+BwwGcbpextI6EndLK3UrVsTnig9cj4MRCAk+dWILXI2KlCy0HwAquH5pZq/zCdXLPo8fxpcdPArBmifT5PQjIInaPhJxmjufW0ji5lIBhUlscsvdfOywHLg4cDnLdSrNrKSiS0PJeNt2KIBBsHwrixbOruG//WYgCgUnRdZYDAGwb8uPBQ+egGSY8DWoASCnFjw/P4bKJPuwcDuA/XbYZhBDccvEo9oyFsbnPyoibWUvhqRNL8IgEV0z2w+sREJBFJDIGFwcOp10Elaw4TK2ksKXf11VpmO1m53AA339pNudYt1oOhkkxvZLCtqHyc6+r5eh8HKeXkrjr9TvwG6/d6hz/9K/tBZANgM+spvH0iWXsnehz2ttMDgZwcjHelvgXjzlwOAAkUYDPI1oxh5UUdymtExaUvn7XEG660Gqe1411NkwQTtYRd0hrBj709efx+NFFAMBDL88BAG68YKTo61lXg6PzMRyYXsNrtw86X7t8sg+v2TbQlk613HLgcGxCXgmrKQ1TK0lcWmQmMqc0e8asmpD337AToxEvVN3ExV1YJ7LNDgKfWkwA51s9jgKyiH3bBuxZ4rTixuFLj5/Ed1+YwRPHl/CjD78e/3HoHC7dEnFEIB+vR8RQUMaXf34KhklxrWve9l/cejFo0Xc1H245cDg2l0/24cGD57DCaxzWzZv3jOBHH3k9rjlvCDuGg/jnO1/blQOlhoIygoqEU4sJmCbFR7/xAv7o3w8AAN5377O47e4nMJ/XntzNfCyNex45hsu2RLCcyOANf/MoXppaw9svHy/7cy/aHEHY58Gn3nEJrt6RrYEQBNK29jdcHDgcm4/cvBuJDMtU4mms60EQCHaPtD7dstEQQnDJeASPH1vES9NrWEpk8OpcHI8fXcRLU2tYjKv48NdfcPoj5fP5R09A1U383e2X4/ffdB4A4G/fdRl+65ptZX/uF9+zD0987E24/arJjol1cXHgcGwuGA3j7XutHd44F4ee5W2XjuH4QgKff/Q42Dr9/33Hsh5+/03n4YnjS7hv/1nMRdN49xeewo/tmIJumLj/xWncvGcE24cC+PBNu/H8n9yMX71yS8UFX5aEhmVHNQoec+BwXHzsly7AlgE/Lh3nMYde5a0Xj+JP7z+EBw+dw+WTfUiqBo7MxXDR5jA+evNuPHZ0EX//8FE8dnQBT55Ywv7Ty/i/v34FFEnAYjyD21wupHYEkhtFx0gVIeQthJAjhJBjhJCPtft8OL3JppAXH715N6QO28VxWsdgUME1O62MoTeevwk377GyjN568SgIIfjozbsxs5bGAwfO4bev3YaLxyP44Neew6d++ArCXglvOH9jjODtiCeAECICuBvAWwHsAfBuQsie9p4Vh8PpVVgA+cYLN+EdV4zjos1hxyK4ftcQXrdjEBMDPvyPWy7AV3/nKly6JYJXzsXwtkvHoEgbYwQvcbfZbdtJEPI6AH9GKb3F/v/HAYBS+pel3rNv3z66f//+Fp0hh8PpJSilOLGYwM7h4pXJac2AblJX8aSGux85jt947aQzV7tTIYT8/+3dW6hc1R3H8e8PT6qo8ZZL8RIbC2KMYk0NWC9tweJD+lJBoQmlifah3or2rbEU+tQHQyuieYihpmgrkopKk9YLKq2ora05GozHg5eI1KOhJhg1ibRU+PdhraHDmXOSM2f2PnvNzO8Dm5lZs8+ftf/Zmf+sPTNrjUbEyiPtV8pnDqcD77U9ngAubqgvZjbkJE1bGCD9NqHd/GPmsX7Vsrq7NaeKuKwETPWpTceQRtIPJe2QtGPv3r1z0C0zs+FUSnGYAJa0PT4D+GDyThGxOSJWRsTKRYsG40MfM7MSlVIcXgLOlnSWpC8Aq4FtDffJzGxoFfGZQ0R8LulHwJPAUcCWiBhruFtmZkOriOIAEBGPAY813Q8zMyvnspKZmRXExcHMzDq4OJiZWYcifiE9G5IOAG/0GGYhsK+C7rScCFS5OnmV8aruW0tVOSw5d3XFHabzr454zl93Wvn6UkQc+bcAEdGXG7CjhBiT4m0uNV7Vfas6hyXnrsZ/k6E5/2qK5/zVmC9fVqrW9oLjVd23qpWcu7mIW4XSc1hy7qD8453T/PXzZaUdMYPJo+qOMeycw9lz7nrj/HWn23z188hhcyExhp1zOHvOXW+cv+50la++HTmYmVl9+nnkYGZmNRmo4iBpiaQ/SxqXNCbp1tx+iqSnJL2Vb0/O7Qvy/gclbWyLM1/SzrZtn6Q7mzquuVRVDvNzayTtkvSqpCckLWzimOZKxbn7bs7bmKQNTRzPXJtF/q6UNJrPsVFJV7TFuii3vy3pLkn9u5hzU6r8qlXTG3Aq8NV8fz7wJmnZ0Q3A+ty+Hrg93z8OuBy4Adh4mLijwDeaPr5+yiFp3q4PgYX58QbSan+NH2Mf5G4B8E9gUX58H/Ctpo+vwPytAE7L988H3m+L9Q/gEtJaMY8Dq5o+vn7bBmrkEBF7IuLlfP8AME5aZe47pP9g5Nur8j6HIuJ54N/TxZR0NrAYeK7Grhejwhwqb8fld20nMMUaHYOkwtx9GXgzIlorWj0NXF1z9xs3i/y9EhGtc2oMOEbS0ZJOBU6IiL9FqhT3t/7GZm6gikM7SUtJ7yz+DnwxIvZAOgFJL/YztQbYmk+yodJLDiPiv8CNwC5SUVgO3Ftjd4vS4/n3NrBM0lJJI6QXtiVH+JuBMov8XQ28EhH/IRWUibbnJnKbdWEgi4Ok44GHgR9HxKc9hlsNPNh7r/pLrzmUNI9UHFYApwGvArdV2slC9Zq7iNhPyt1W0oj1XeDzKvtYsm7zJ+k84Hbg+lbTFLsN3Zu7Xg1cccgvSg8DD0TEI7n5X3moSb79cIaxvgKMRMRoLZ0tVEU5vBAgInbnUdfvgUtr6nIxqjr/ImJ7RFwcEZeQ5hB7q64+l6Tb/Ek6A3gUWBsRu3PzBGmp4ZYplx22wxuo4pCvbd8LjEfEHW1PbQPW5fvrgD/MMOQahmzUUGEO3weWS2pN8HUl6RrywKry/JO0ON+eDNwE/Lra3pan2/xJOgn4E3BbRLzQ2jlfejog6Ws55lpm/n/eWpr+RLzKjfTNjyBdwtiZt2+Tvv3xDOnd1zPAKW1/8y7wEXCQ9I5jedtz7wDLmj6ufs0h6Vs44znWdmBB08fXR7l7EHg9b6ubPrYS8wf8DDjUtu9OYHF+biXwGrAb2Ej+wa+3mW/+hbSZmXUYqMtKZmZWDRcHMzPr4OJgZmYdXBzMzKyDi4OZmXVwcTCrgaQbJK3tYv+lkl6rs09m3RhpugNmg0bSSERsarofZr1wcTCbQp747QnSxG8rSNNHrwXOBe4Ajgf2AddGxB5JfwH+ClwGbJM0HzgYEb+UdCGwCTiW9KOsH0TEfkkXAVuAz4Dn5+7ozI7Ml5XMpncOsDkiLgA+BW4G7gauiYjWC/sv2vY/KSK+GRG/mhTnfuAnOc4u4Oe5/TfALZHmTzIrikcOZtN7L/4/Z8/vgJ+SFpV5Ki8sdhSwp23/rZMDSDqRVDSezU33AQ9N0f5bYFX1h2A2Oy4OZtObPLfMAWDsMO/0D3URW1PENyuGLyuZTe9MSa1CsAZ4EVjUapM0L68lMK2I+ATYL+nruen7wLMR8THwiaTLc/v3qu++2ex55GA2vXFgnaR7SDOC3g08CdyVLwuNAHeSlqg8nHXAJknHkmb6vS63XwdskfRZjmtWDM/KajaF/G2lP0bE+Q13xawRvqxkZmYdPHIwM7MOHjmYmVkHFwczM+vg4mBmZh1cHMzMrIOLg5mZdXBxMDOzDv8Dw8BLFS0f/4wAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'][-200:].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "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": 14, + "metadata": {}, + "outputs": [], + "source": [ + "year = []\n", + "yearly_incidence = []\n", + "for week1, week2 in zip(first_sept_week[:-1],\n", + " first_sept_week[1:]):\n", + " one_year = sorted_data['inc'][week1:week2-1]\n", + "# assert abs(len(one_year)-52) < 2\n", + " yearly_incidence.append(one_year.sum())\n", + " year.append(week2.year)\n", + "yearly_incidence = pd.Series(data=yearly_incidence, index=year)" + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 15, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAD8CAYAAAC2PJlnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAF6FJREFUeJzt3X+QndV93/H3Vz8sYRvZEkgEJGyRsYZakDE2O4LUaRpMg+Q4Y9EJZFQIqA0dWpe0JJMpFoUOHUNnTCcTJww1NVNsC2OMKTEDLVZsWY4nNKUSK2MHC1mREoiQkdF6REF4BoHQt3/cs+hqWe3ee3R1f+y+XzM7e++5z3n23CPtfp7nnPM8NzITSZJqzOh1AyRJg8sQkSRVM0QkSdUMEUlSNUNEklTNEJEkVTNEJEnVDBFJUjVDRJJUbVavG9Bpp556ai5durTXzZCkgbJ169afZebCdutNuRBZunQpw8PDvW6GJA2UiPj7mnoOZ0mSqhkikqRqhogkqZohIkmqZohIkqoZIpJOiH2vvMZvf+EJ9h14rddN0QlkiEg6Ie7YtJMnn9vPHd/Z2eum6ASacteJaHz7XnmN3/vaU9x5xYdZdPLcXjdHU9jZN2/g4KHDbz2/b/Nu7tu8mzmzZrDjto/3sGU6ETwTmSY8KlS3PH7DRXzyvDOYO7vx52Xu7BmsPu8MHv/0RW/b1iGvweeZyBTnUaG6bdG8uZw8ZxYHDx1mzqwZHDx0mJPnzBr3DLj54Oa2f/pLPWitjpchMsU9fsNF3PbN7Xx720957Y3DzJ09g5Xn/AI3feKDvW6aprCfvXqQKy94P1eseB/3b9nNyJgzDQ9upg5DZIpr56hQ6pQvXDX01uPbLj33ba97cDN1GCLTwGRHhVK3eXAzdRgi08BkR4VSL3hwMzVEZva6DR01NDSU3gpektoTEVszc2jyLY/W0hLfiPiDiNgWET+KiK9FxNyIWBARGyNiZ/k+v2n7GyNiV0TsiIiVTeXnR8TT5bU7IiJK+ZyI+Hop3xwRS5vqrC0/Y2dErG33DUqSTpxJQyQiFgP/DhjKzHOBmcAaYB2wKTOXAZvKcyJieXn9HGAV8PmImFl2dxdwLbCsfK0q5dcAL2XmB4DPAbeXfS0AbgEuAFYAtzSHlSSpt1q92HAWcFJEzALeCbwArAbWl9fXA5eWx6uBBzLzYGY+C+wCVkTE6cC8zHwiG2No946pM7qvh4CLy1nKSmBjZu7PzJeAjRwJHklSj00aIpn5E+CPgN3AXuDlzPw2cFpm7i3b7AUWlSqLgeebdrGnlC0uj8eWH1UnMw8BLwOnTLCvo0TEtRExHBHDIyMjk70lSVKHtDKcNZ/GmcJZwBnAuyLidyaqMk5ZTlBeW+dIQebdmTmUmUMLF7b9OfOS1HHT5ZYurQxn/RPg2cwcycw3gG8A/xB4sQxRUb7vK9vvAc5sqr+ExvDXnvJ4bPlRdcqQ2XuA/RPsS5o2pssfo6lmutyvrpUQ2Q1cGBHvLPMUFwPbgUeB0dVSa4FHyuNHgTVlxdVZNCbQt5QhrwMRcWHZz9Vj6ozu6zLgu2Xe5FvAJRExv5wRXVLKpGljuvwxmirOvnkDS9c9xn2bd5PZuKXL0nWPcfbNG3rdtBNi0osNM3NzRDwEfB84BDwF3A28G3gwIq6hETSXl+23RcSDwDNl++sy882yu08BXwZOAjaUL4B7gK9ExC4aZyBryr72R8StwJNlu89k5v7jesfSgPD+UoNput3SxYsNpT6175XXjvnHyNuD9LebHn6a+7fs5h0zZ/D6m4e5csX7+v4uxbUXG3rbE6lPeX+pwTWdbuliiGja6LdPd2ylPdPpj9FUMp3uV2eIaNrotw9AaqU90+mPkQaTcyKa8sZOUI/q1QR1v7WnRr+d1en4ndAbMEqDrJ3P/J6O7anhsmONcjhLU16/TVD3W3va4bJjjeWZiKaF0Qnqh//NR7nygvcz8upB21NhKpxFqbM8E9G00G8T1P3WnlYN8lmUTgxDRFJbXHasZq7OkiS5OkuS1H2GiKS3eNt5tcsQkfQWr/9Qu5xYl+T1H6rmmYgkr/9QNUNE0rS5/sM5n84zRCQBg3sVfTuc8+k8rxORNOVNhTsnn2heJyJJx9DOnI9DXu0xRNTX/IVWJ7Qz5+OQV3tc4qu+1sqn//kBSWrFZPf8cplzHedE1JfaGcO++eGn+eqW3Vy54n198bG3Gkz7XnmN2765nW9v+ymvvXGYubNnsPKcX+CmT3xwWhyc1M6JeCaivvT4DRcd8xd6lEeO6qTpssy505wTUV9q5RfaC+TUadNhmXOneSaitnRz/mGyMWyPHNVpg/phYb1kiKgtrUx0d0orv9B+QJLUW06sqyXtXqzliilpsHixoU6oducfXGsvTQ+GiN4y0YV9rc4/nH3zBpaue4z7Nu8ms7Fiaum6xzj75g3dehuaYrzgtL8ZInrLZGcPraxcccWUOs2z2v7mnIg6fnO6mx5+mvu37OYdM2fw+puHvQiwTwzaPJU3Tewu50RUrdNnD661775WhnwG7Yjes9rB4BJfdfx6C9fad99ES68H9cr+dv5fDtpZ1lRiiAjweotB1UpAtHILmX7V6v/Lbl6/pKM5JyI1aeWItp+Oelu9aeBUnady3qRznBOROqCVeYN+mltodchnqs5TOW/Sew5nSbQ2LNSvcwutDPlM1Xkq75/Wey2diUTEeyPioYj4cURsj4hfjogFEbExInaW7/Obtr8xInZFxI6IWNlUfn5EPF1euyMiopTPiYivl/LNEbG0qc7a8jN2RsTazr116YhWjmj79aj3C1cNcdul57L8jHncdum5RwXGdNCLsywvgDyi1eGsPwX+PDP/AfAhYDuwDtiUmcuATeU5EbEcWAOcA6wCPh8RM8t+7gKuBZaVr1Wl/Brgpcz8APA54PayrwXALcAFwArgluawkjqllSNaj3r7Uy9CtJ+GNHtt0uGsiJgH/CrwzwEy83Xg9YhYDfxa2Ww98D3g08Bq4IHMPAg8GxG7gBUR8RwwLzOfKPu9F7gU2FDq/Keyr4eAO8tZykpgY2buL3U20gierx3He5bG1cqwkKvYprd+HdLspVbmRH4RGAG+FBEfArYC1wOnZeZegMzcGxGLyvaLgf/bVH9PKXujPB5bPlrn+bKvQxHxMnBKc/k4daSOamXeYKrOLag1g7xc+kRpZThrFvAR4K7M/DDwc8rQ1THEOGU5QXltnSM/MOLaiBiOiOGRkZEJmiZJ9RzSfLtWQmQPsCczN5fnD9EIlRcj4nSA8n1f0/ZnNtVfArxQypeMU35UnYiYBbwH2D/Bvo6SmXdn5lBmDi1cuLCFtyRJdfp1uXSvJvsnDZHM/CnwfEScXYouBp4BHgVGV0utBR4pjx8F1pQVV2fRmEDfUoa+DkTEhWW+4+oxdUb3dRnw3WxcBfkt4JKImF8m1C8pZZLUE/26Gq5Xk/2tXifyb4GvRsQ7gL8D/gWNAHowIq4BdgOXA2Tmtoh4kEbQHAKuy8w3y34+BXwZOInGhProh0zcA3ylTMLvp7G6i8zcHxG3Ak+W7T4zOskuSer9ZL+3PZGkAdbqrW8m421PNHC8YEs6fr2e7DdE1DNesCV1Ri8n+x3OUtd551Wp/zicNY0N2rBQv96DSlL7DJEpYNCGhXo9hiv1k0E7CBzLW8EPsF4v7Tse3oNKahj0T2V0TmSAdWppn6Tua2dusBufpumcyDTksJA0uNqZG+znIWuHswacw0LSYGrlIHAQhqwNkQHnrcmlwTXZQeAg3HreEJGkHpnsIHAQhqwNEUnqY/0+ZO3qLOkE6MZqGqmTXJ0l9ZF+Xk0jdZLDWVIHDcJqGqmTPBOROsj7gmm6MUSkDhqE1TRSJzmcJXVYv6+mkTrJ1VmSJFdnSZK6zxCRJFUzRCRJ1QwRSVI1Q0SSVM0QkSRVM0QkSdUMEUlSNUNEklTNEJEkVTNEJEnVDBFJUjVDRJJUzRCRJFUzRCRJ1QwRSVI1Q0SSVM0QkSRVM0QkSdVaDpGImBkRT0XE/yrPF0TExojYWb7Pb9r2xojYFRE7ImJlU/n5EfF0ee2OiIhSPicivl7KN0fE0qY6a8vP2BkRazvxpiVJndHOmcj1wPam5+uATZm5DNhUnhMRy4E1wDnAKuDzETGz1LkLuBZYVr5WlfJrgJcy8wPA54Dby74WALcAFwArgFuaw0qS1FsthUhELAE+Afz3puLVwPryeD1waVP5A5l5MDOfBXYBKyLidGBeZj6RmQncO6bO6L4eAi4uZykrgY2ZuT8zXwI2ciR4JEk91uqZyJ8ANwCHm8pOy8y9AOX7olK+GHi+abs9pWxxeTy2/Kg6mXkIeBk4ZYJ9SZL6wKQhEhG/CezLzK0t7jPGKcsJymvrNLfx2ogYjojhkZGRFpspSTperZyJfBT4ZEQ8BzwAfCwi7gNeLENUlO/7yvZ7gDOb6i8BXijlS8YpP6pORMwC3gPsn2BfR8nMuzNzKDOHFi5c2MJbkiR1wqQhkpk3ZuaSzFxKY8L8u5n5O8CjwOhqqbXAI+Xxo8CasuLqLBoT6FvKkNeBiLiwzHdcPabO6L4uKz8jgW8Bl0TE/DKhfkkpkyT1gVnHUfezwIMRcQ2wG7gcIDO3RcSDwDPAIeC6zHyz1PkU8GXgJGBD+QK4B/hKROyicQaypuxrf0TcCjxZtvtMZu4/jjZLkjooGgf8U8fQ0FAODw/3uhmSNFAiYmtmDrVbzyvWJUnVDBFJUjVDRJJUzRCRJFUzRCRJ1QwRSVI1Q0SSVM0QkSRVM0QkSdUMEUlSNUNEklTNEJEkVTNEJEnVDBFJUjVDRJJUzRCRJFUzRCRJ1QwRSVI1Q0SSVM0QkSRVM0QkSdUMEUlSNUNEklTNEJEkVTNEJEnVDBFJUjVDRJJUzRCRJFUzRCRJ1QwRSVI1Q0SSVM0QkSRVM0QkSdUMEUlSNUNEklTNEJEkVTNEJEnVDBFJUrVJQyQizoyIv4iI7RGxLSKuL+ULImJjROws3+c31bkxInZFxI6IWNlUfn5EPF1euyMiopTPiYivl/LNEbG0qc7a8jN2RsTaTr55SdLxaeVM5BDwh5n5QeBC4LqIWA6sAzZl5jJgU3lOeW0NcA6wCvh8RMws+7oLuBZYVr5WlfJrgJcy8wPA54Dby74WALcAFwArgFuaw0qS1FuThkhm7s3M75fHB4DtwGJgNbC+bLYeuLQ8Xg08kJkHM/NZYBewIiJOB+Zl5hOZmcC9Y+qM7ush4OJylrIS2JiZ+zPzJWAjR4JHktRjbc2JlGGmDwObgdMycy80ggZYVDZbDDzfVG1PKVtcHo8tP6pOZh4CXgZOmWBfkqQ+0HKIRMS7gT8Dfj8zX5lo03HKcoLy2jrNbbs2IoYjYnhkZGSCpkmSOqmlEImI2TQC5KuZ+Y1S/GIZoqJ831fK9wBnNlVfArxQypeMU35UnYiYBbwH2D/Bvo6SmXdn5lBmDi1cuLCVtyRJ6oBWVmcFcA+wPTP/uOmlR4HR1VJrgUeayteUFVdn0ZhA31KGvA5ExIVln1ePqTO6r8uA75Z5k28Bl0TE/DKhfkkpkyT1gVktbPNR4Crg6Yj4QSn7D8BngQcj4hpgN3A5QGZui4gHgWdorOy6LjPfLPU+BXwZOAnYUL6gEVJfiYhdNM5A1pR97Y+IW4Eny3afycz9le9VktRh0TjgnzqGhoZyeHi4182QpIESEVszc6jdel6xLkmqZohIkqoZIpKkaoaIJKmaISJJqmaISJKqGSKSpGqGiCSpmiEiSapmiEiSqhkikqRqhogkqZohIkmqZohIkqoZIpKkaoaIJKmaISJJqmaISJKqGSKSpGqGiCSpmiEiSapmiEiSqhkikqRqhogkqZohIkmqZohIkqoZIpKkaoaIJKmaISJJqmaISJKqGSKSpGqGiCSpmiEiSapmiEiSqhkikqRqhogkqZohIkmqZohIkqoNRIhExKqI2BERuyJiXa/bI0lq6PsQiYiZwH8FPg4sB/5ZRCw/ET9r3yuv8dtfeIJ9B17ri23a2U6SeqHvQwRYAezKzL/LzNeBB4DVJ+IH3bFpJ08+t587vrOzL7ZpZztJ6oXIzF63YUIRcRmwKjP/ZXl+FXBBZv7eeNsPDQ3l8PBwWz/j7Js3cPDQ4beVz5k1gx23fbzr27SznSR1QkRszcyhdusNwplIjFN2VPJFxLURMRwRwyMjI23/gMdvuIhPnncGc2c3umPu7BmsPu8MHv/0RT3Zpp3tJKmXBiFE9gBnNj1fArzQvEFm3p2ZQ5k5tHDhwrZ/wKJ5czl5ziwOHjrMnFkzOHjoMCfPmcWik+f2ZJt2tpOkXprV6wa04ElgWUScBfwEWANc0ekf8rNXD3LlBe/nihXv4/4tuxkZZyK7m9u0s50k9Urfz4kARMRvAH8CzAS+mJn/+Vjb1syJSNJ0VzsnMghnImTmN4Fv9rodkqSjDcKciCSpTxkikqRqhogkqZohIkmqZohIkqoNxBLfdkTECPD347x0KvCzLjenUwa17ba7u2x3d021dr8/M9u+WnvKhcixRMRwzRrofjCobbfd3WW7u8t2NzicJUmqZohIkqpNpxC5u9cNOA6D2nbb3V22u7tsN9NoTkSS1HnT6UxEktRhAx0iEfHFiNgXET9qKvtQRDwREU9HxP+MiHmlfHZErC/l2yPixqY634uIHRHxg/K1qI/a/Y6I+FIp/2FE/FpTnfNL+a6IuCMixvsAr35sd7f7+8yI+Ivy774tIq4v5QsiYmNE7Czf5zfVubH0646IWNlU3rU+73C7u9bn7bY7Ik4p278aEXeO2Vff9vck7e7n/v71iNha+nVrRHysaV/t93dmDuwX8KvAR4AfNZU9Cfzj8vh3gVvL4yuAB8rjdwLPAUvL8+8BQ33a7uuAL5XHi4CtwIzyfAvwyzQ+/XED8PEBaXe3+/t04CPl8cnA3wDLgf8CrCvl64Dby+PlwA+BOcBZwN8CM7vd5x1ud9f6vKLd7wJ+BfjXwJ1j9tXP/T1Ru/u5vz8MnFEenwv85Hj6e6DPRDLzL4H9Y4rPBv6yPN4I/Nbo5sC7ImIWcBLwOvBKN9o5VpvtXg5sKvX2Af8PGIqI04F5mflENv717wUu7fd2n8j2HUtm7s3M75fHB4DtwGJgNbC+bLaeI/23msYBx8HMfBbYBazodp93qt0nqn3H0m67M/Pnmfm/gaM+da3f+/tY7e62inY/lZmjnw67DZgbEXNq+3ugQ+QYfgR8sjy+nCMfrfsQ8HNgL7Ab+KPMbP6D+KVy2vkfT/Sw0DEcq90/BFZHxKxofLrj+eW1xTQ+OnjUnlLWbe22e1RP+jsiltI4EtsMnJaZe6Hxi0jjjAka/fh8U7XRvu1Znx9nu0d1vc9bbPex9Ht/T2YQ+vu3gKcy8yCV/T0VQ+R3gesiYiuNU7vXS/kK4E3gDBqn+n8YEb9YXrsyM38J+Efl66ruNhk4dru/SOMfc5jGpzv+H+AQjdPNsXqx1K7ddkOP+jsi3g38GfD7mTnRWeix+rYnfd6BdkMP+ryNdh9zF+OU9VN/T6Tv+zsizgFuB/7VaNE4m03a31MuRDLzx5l5SWaeD3yNxrgwNOZE/jwz3yjDK39FGV7JzJ+U7weA++nNEMC47c7MQ5n5B5l5XmauBt4L7KTxB3pJ0y6WAC+M3W8ftrsn/R0Rs2n8gn01M79Ril8sp/CjQyf7Svkejj5rGu3brvd5h9rd9T5vs93H0u/9fUz93t8RsQR4GLg6M0f/Rlb195QLkdFVEBExA7gZ+G/lpd3Ax6LhXcCFwI/LcMuppc5s4DdpDNH0Rbsj4p2lvUTErwOHMvOZcnp6ICIuLKfKVwOP9Hu7e9HfpX/uAbZn5h83vfQosLY8XsuR/nsUWFPGic8ClgFbut3nnWp3t/u8ot3jGoD+PtZ++rq/I+K9wGPAjZn5V6MbV/f3ZDPv/fxF48h3L/AGjRS9BriexuqEvwE+y5ELKt8N/A8aE0nPAP8+j6yw2Ar8dXntTykrWvqk3UuBHTQmy75D406bo/sZovGf82+BO0fr9HO7e9Tfv0LjtPyvgR+Ur98ATqEx+b+zfF/QVOem0q87aFqh0s0+71S7u93nle1+jsaijVfL/63lA9Lfb2t3v/c3jYO9nzdt+wNgUW1/e8W6JKnalBvOkiR1jyEiSapmiEiSqhkikqRqhogkqZohIkmqZohIkqoZIpKkav8fPVSNLoPTwy4AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.plot(style='*')" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "1986 0\n", + "1987 0\n", + "1988 0\n", + "1989 0\n", + "1990 0\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": 16, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "yearly_incidence.sort_values()" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 17, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEICAYAAABPgw/pAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAEbZJREFUeJzt3X+QXXV5x/H3Q6JtYDGiwcUGdSljqY5RkbX4Y4q7QVsRf7SOrVBUqD/ijFUzls40dXSoVcdMlY4Og+2kFaGFcavBFgzWwqiLv5EErAtGhULEBA1SJLoUCZGnf5yTugmwN3vOzZ67332/ZnZyz91z7nnuk3M/+73nnnNuZCaSpIXvkK4LkCT1h4EuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKsTS+VzZihUrcmRkpNGy99xzD4cddlh/CyqI/enNHs3O/syuy/5s2bLlzsw8std88xroIyMjbN68udGyk5OTjI2N9beggtif3uzR7OzP7LrsT0T84EDmc5eLJBXCQJekQhjoklQIA12SCmGgS1IhegZ6RFwQEXdExA0z7ntMRFwVETfV/x5xcMuUJPVyICP0C4EX73ffOuDzmflk4PP1tCSpQz0DPTO/BNy1392vAC6qb18E/EGf65IkzVEcyHeKRsQIsCkzn1ZP352Zj57x+59m5kPudomINcAagOHh4RMmJiYaFTo9Pc3Q0FCjZRcD+9ObPdrX1I5d+0wPL4Od987PuletXD4/K+qjLref8fHxLZk52mu+g36maGZuADYAjI6OZtMzrTyLbXb2pzd7tK+z1l2xz/TZq/Zw7tT8nDy+7YyxeVlPPy2E7afpUS47I+LxAPW/d/SvJElSE00D/XLgzPr2mcBl/SlHktTUgRy2+Ang68BxEbE9It4ArAdeFBE3AS+qpyVJHeq5wywzT3+YX53c51okSS14pqgkFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEK0CPSLeERE3RsQNEfGJiPj1fhUmSZqbxoEeESuBtwOjmfk0YAlwWr8KkyTNTdtdLkuBZRGxFDgUuL19SZKkJhoHembuAD4E3Ab8CNiVmVf2qzBJ0txEZjZbMOII4FLg1cDdwKeAjZl58X7zrQHWAAwPD58wMTHRaH3T09MMDQ01WnYxsD+92aN9Te3Ytc/08DLYee/8rHvVyuXzs6I+6nL7GR8f35KZo73mW9piHS8Ebs3MnwBExKeB5wH7BHpmbgA2AIyOjubY2FijlU1OTtJ02cXA/vRmj/Z11ror9pk+e9Uezp1qEwkHbtsZY/Oynn5aCNtPm33otwHPiYhDIyKAk4Gt/SlLkjRXbfahXwNsBK4DpurH2tCnuiRJc9Tq/VVmngOc06daJEkteKaoJBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhZif75uSpBlG9vv6u/mybf2pnax3vjhCl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEK0CvSIeHREbIyI70bE1oh4br8KkyTNTduvoPsI8LnMfFVEPBI4tA81SZIaaBzoEfEo4CTgLIDM3A3s7k9ZkqS5isxstmDEM4ENwHeAZwBbgLWZec9+860B1gAMDw+fMDEx0Wh909PTDA0NNVp2MbA/vQ1qj6Z27Oq6BACGl8HOe7uuYnC17c+qlcsbLzs+Pr4lM0d7zdcm0EeBbwDPz8xrIuIjwM8y890Pt8zo6Ghu3ry50fomJycZGxtrtOxiYH96G9Qejay7ousSADh71R7OnWq7F7Zcbfuzbf2pjZeNiAMK9DYfim4HtmfmNfX0RuBZLR5PktRC40DPzB8DP4yI4+q7Tqba/SJJ6kDb91dvAy6pj3C5BfjT9iVJkppoFeiZ+S2g534dSdLB55miklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFWLBfIHg1I5dnNXRdy+2+S5ASZovjtAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYVoHegRsSQiro+ITf0oSJLUTD9G6GuBrX14HElSC60CPSKOBk4F/qk/5UiSmorMbL5wxEbgA8DhwF9k5ksfYp41wBqA4eHhEyYmJhqt6467drHz3saltrJq5fJuVjwH09PTDA0NdV3GQBvUHk3t2NV1CQAML6Oz19hC0LY/bXJkfHx8S2aO9pqv8ZdER8RLgTsyc0tEjD3cfJm5AdgAMDo6mmNjDzvrrM675DLOnermO623nTHWyXrnYnJykqa9XSwGtUddffn5/s5etaez19hC0LY/85EjbXa5PB94eURsAyaA1RFxcV+qkiTNWeNAz8y/ysyjM3MEOA34Qma+pm+VSZLmxOPQJakQfdlhlpmTwGQ/HkuS1IwjdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKkTjQI+IJ0TEFyNia0TcGBFr+1mYJGlulrZYdg9wdmZeFxGHA1si4qrM/E6fapMkzUHjEXpm/igzr6tv/xzYCqzsV2GSpLmJzGz/IBEjwJeAp2Xmz/b73RpgDcDw8PAJExMTjdZxx1272HlvuzqbWrVyeTcrnoPp6WmGhoa6LmOgDWqPpnbs6roEAIaX0dlrbCFo2582OTI+Pr4lM0d7zdc60CNiCLgaeH9mfnq2eUdHR3Pz5s2N1nPeJZdx7lSbPUTNbVt/aifrnYvJyUnGxsa6LmOgDWqPRtZd0XUJAJy9ak9nr7GFoG1/2uRIRBxQoLc6yiUiHgFcClzSK8wlSQdXm6NcAvgYsDUz/65/JUmSmmgzQn8+8FpgdUR8q/55SZ/qkiTNUeMdQpn5FSD6WIskqQXPFJWkQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVolWgR8SLI+J7EXFzRKzrV1GSpLlrHOgRsQQ4HzgFeCpwekQ8tV+FSZLmps0I/XeAmzPzlszcDUwAr+hPWZKkuYrMbLZgxKuAF2fmG+vp1wInZuZb95tvDbCmnjwO+F7DWlcAdzZcdjGwP73Zo9nZn9l12Z8nZeaRvWZa2mIF8RD3PeivQ2ZuADa0WE+1sojNmTna9nFKZX96s0ezsz+zWwj9abPLZTvwhBnTRwO3tytHktRUm0C/FnhyRBwTEY8ETgMu709ZkqS5arzLJTP3RMRbgf8ElgAXZOaNfavswVrvtimc/enNHs3O/sxu4PvT+ENRSdJg8UxRSSqEgS5JhVgwgR4RMfNfSdK+BjrQI2JJRKyNiI3AnwGkO/37JiJWR8QxXdcxqOxPb/ZodvPdnzYnFs2HFwK/D3wYWFdfP+aSzPRsthbqa+5MAHcDD0TEOZl5dcdlDQz705s9ml1X/RnoETrwOuDKzLwSeBfweODl3Za08ETE0RHxqBl3vRq4NDNPotroTo+I53ZTXffm0p/FusvPbWh2g9KfQQ/0rwMj9e3rgZuA4+uRunqIiKdExGeBrwB/ExF7L572C+DQ+vYnqa5P8ZzFFlZN+rPYdvm5Dc1u0Poz6IF+CzAUEY/NzHvraagu16uHEBGHzZh8JrA9M0eALwAfqu+/C7gvIg7PzLuA7wPD/OqPZ7Fm6c8XeXB/hur+3MQi6Q+4DfUyyP0Z9ECfAu4DTq6nf0rVlJ90VtEAiogjIuLCiLgWWB8RR9YjgacDX61HlpcDd0fEqVQb1+HAqvohbqK6ktzuLuo/2Gbpzyp+1Z/LeHB/nl4/xPcpuD/gNtTLQunPoAf6DuAbwNvr6Z9QBfrPO6toMJ0E7AFeQnUVzHcCj6L6/z1qxm6Ci4A/Ab5J1cNTADLz68Bq4GfzW/a8ebj+LMH+7OU2NLuF0Z/MHPgf4BLgCuAO4G3A0q5r6qgPS4A3A1dTXWN+RX3/J4G317ePAdbXv3821b69JfXvhqj+KC4BVgJbgbcCHwc+Chza9XO0P/bI/jT/GfQR+l5nUf1FfFJmnpeZezqupysvpTrK5z3Ac4G/re+/CnheffuHwJeBUzLzWqrRxDhAZk4D1wDPzswdwGup3vH8GHhXZv7vPD2Pg8X+9GaPZreg+zPox6EDkJn3A//VdR3zYe+RFBHxbKq3bl8GrsjM+4DfAm7JzC9ExK3AByPi94AtwB9GxIrMvDMibgLuiYgnAucBr4mIx1Fds/5/qN4Okpmbgc3z/iRbsD8Hzh49tJK3oYUyQl8UZmxoJwEXUB369ELgA/UsDwDfj4hlmXkr1du5p1Ptq7ud6thXgF9SveU7BLiU6rKfZwAnABsy84F5ekp9VT/vjIgx7M9DiojH1v8+D7gQe7SPiHhkvQ29gGo3SFH9WRAj9JJFxKFUJ1CdBHwxIi4EngX8Q2aeHxFHAFdHxPFU3xJ1PNVbuG1Ub+1eTnX0z78DH4iIi4EnAcuB2+t3N1+LiJct4Bfh4VQnZ+wEXk/VA/vDPoOAPwLWAnsi4s1UX+L+UXsEEfEI4E1UH1Buj4hzqQ43/PvS+uMIvUMRcRSwCRgD/oXqQ5ZXUu272wOQmT8FLqM60ucq4HHAU+qH+BLwAmB3Zm4CPgZsBM6n+rT9/r3r6npDa2kZ8GvAsRGxAjiWanS06PtTh/ly4I+BD2fmWGZ+j2qkuHeeRd0jqutAraa6hEhSvcZWU42yy+rPwfzE1Z+en6gvA06cMX0W1YcwZwLfnHH/b1CNBKDaOCeAI+rlPwM8cca8K7p+XgehT2cCHwTeDbwBeBlwrf35/+f0FuC99e29X1rzSuAae5TUz+919e03UB0p9+oSX2OO0Lv1C+CbM04H/jZVwF9ENRo9CiAzbwdujIgTM/N84GbgX6neHn4lM2/b+4BZ0IXLZvTlEOC/gRuB1Zn5GeCYiBiGxdufGe4EfjcizgC2RMRFVKPP367f0Sz2Hv0HcGZEfAr4a+AZwA1U34l8JJTTH/ehdyj3Dql+FVxrgU/Xtz8LvAP4y4h4DPAD4Ef1794DHAfcnJm/mL+K59fe/lCdzHEa1TG+ayPiXVSHga0B3rtY+zPD9VS7pE6kOqv6FOrdBMDrI+JDwKNZpD3KzI9GxG6qww6Ppzou/FVU57W8KSLWU0h/HKEPgMzMiDgaOIpqNAHVSIKI2AR8nurEhdvq+e/PzBsW0obWVEQMUY1A/xH4HNURBycCpwNHRMRnWMT9qd1GFd5Ls9ofvKm+bxPVBaIuxx49FZisR9cfp/oM5hNUu1SK6Y8j9MFxPNW1a26NiDdSvdV7J9W+vu9m5nVdFtehPVShdD/VJZR/CbwzM6eAP4+I1wA3Zub1HdbYqcy8rx5lvq++6xCqD9bXZ+b1i71HUV2ddTvVkWT/DNxDdYbnWzLz1pL6s/cDFHUsIr4K/CbVoVK3A+/JzG93WtQAqk/keCUwkZk/7rqeQRIR76faho4HrgTOqUfsi15EHEt1rPhuqh79G/C+rM7sLIaBPgDq42TPAW4FLs7qjDXNUI+yHkg32FlFxHHADxbaroL5UA8Gngx8LavLcRfHQJekQvihqCQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5Jhfg/wpkCdy2JX6AAAAAASUVORK5CYII=\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 +2305,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