diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb index 0bbbe371b01e359e381e43239412d77bf53fb1fb..d400aba7754470dfacd1878ece4d4fe267541657 100644 --- a/module3/exo2/exercice.ipynb +++ b/module3/exo2/exercice.ipynb @@ -1,5 +1,2264 @@ { - "cells": [], + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Incidence de la varicelle" + ] + }, + { + "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_local = \"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
020204774983294370238511FRFrance
12020467375219635541639FRFrance
22020457369620165376639FRFrance
320204474391237564077410FRFrance
420204374376250562477410FRFrance
52020427400019796021639FRFrance
62020417396120995823639FRFrance
7202040720786753481315FRFrance
8202039710492371861213FRFrance
9202038722537823724315FRFrance
10202037715844052763204FRFrance
1120203679191001738102FRFrance
12202035782801694102FRFrance
13202034722723714173306FRFrance
14202033712841772391204FRFrance
15202032726506894611417FRFrance
16202031713031002506204FRFrance
1720203071385752695204FRFrance
182020297841101672102FRFrance
19202028772801515102FRFrance
2020202779861491823102FRFrance
21202026769401454102FRFrance
2220202572280597001FRFrance
2320202473880959102FRFrance
24202023755811115102FRFrance
2520202272770633001FRFrance
262020217602361168102FRFrance
272020207824201628102FRFrance
2820201973100753001FRFrance
292020187849981600102FRFrance
.................................
15341991267176081130423912312042FRFrance
15351991257161691070021638281838FRFrance
15361991247161711007122271281739FRFrance
1537199123711947767116223211329FRFrance
1538199122715452995320951271737FRFrance
1539199121714903897520831261636FRFrance
15401991207190531274225364342345FRFrance
15411991197167391124622232291939FRFrance
15421991187213851388228888382551FRFrance
1543199117713462887718047241632FRFrance
15441991167148571006819646261834FRFrance
1545199115713975978118169251832FRFrance
1546199114712265768416846221430FRFrance
154719911379567604113093171123FRFrance
1548199112710864733114397191325FRFrance
15491991117155741118419964271935FRFrance
15501991107166431137221914292038FRFrance
1551199109713741878018702241533FRFrance
1552199108713289881317765231531FRFrance
1553199107712337807716597221529FRFrance
1554199106710877701314741191226FRFrance
1555199105710442654414340181125FRFrance
15561991047791345631126314820FRFrance
15571991037153871048420290271836FRFrance
15581991027162771104621508292038FRFrance
15591991017155651027120859271836FRFrance
15601990527193751329525455342345FRFrance
15611990517190801380724353342543FRFrance
1562199050711079666015498201228FRFrance
15631990497114302610205FRFrance
\n", + "

1564 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202047 7 4983 2943 7023 8 5 \n", + "1 202046 7 3752 1963 5541 6 3 \n", + "2 202045 7 3696 2016 5376 6 3 \n", + "3 202044 7 4391 2375 6407 7 4 \n", + "4 202043 7 4376 2505 6247 7 4 \n", + "5 202042 7 4000 1979 6021 6 3 \n", + "6 202041 7 3961 2099 5823 6 3 \n", + "7 202040 7 2078 675 3481 3 1 \n", + "8 202039 7 1049 237 1861 2 1 \n", + "9 202038 7 2253 782 3724 3 1 \n", + "10 202037 7 1584 405 2763 2 0 \n", + "11 202036 7 919 100 1738 1 0 \n", + "12 202035 7 828 0 1694 1 0 \n", + "13 202034 7 2272 371 4173 3 0 \n", + "14 202033 7 1284 177 2391 2 0 \n", + "15 202032 7 2650 689 4611 4 1 \n", + "16 202031 7 1303 100 2506 2 0 \n", + "17 202030 7 1385 75 2695 2 0 \n", + "18 202029 7 841 10 1672 1 0 \n", + "19 202028 7 728 0 1515 1 0 \n", + "20 202027 7 986 149 1823 1 0 \n", + "21 202026 7 694 0 1454 1 0 \n", + "22 202025 7 228 0 597 0 0 \n", + "23 202024 7 388 0 959 1 0 \n", + "24 202023 7 558 1 1115 1 0 \n", + "25 202022 7 277 0 633 0 0 \n", + "26 202021 7 602 36 1168 1 0 \n", + "27 202020 7 824 20 1628 1 0 \n", + "28 202019 7 310 0 753 0 0 \n", + "29 202018 7 849 98 1600 1 0 \n", + "... ... ... ... ... ... ... ... \n", + "1534 199126 7 17608 11304 23912 31 20 \n", + "1535 199125 7 16169 10700 21638 28 18 \n", + "1536 199124 7 16171 10071 22271 28 17 \n", + "1537 199123 7 11947 7671 16223 21 13 \n", + "1538 199122 7 15452 9953 20951 27 17 \n", + "1539 199121 7 14903 8975 20831 26 16 \n", + "1540 199120 7 19053 12742 25364 34 23 \n", + "1541 199119 7 16739 11246 22232 29 19 \n", + "1542 199118 7 21385 13882 28888 38 25 \n", + "1543 199117 7 13462 8877 18047 24 16 \n", + "1544 199116 7 14857 10068 19646 26 18 \n", + "1545 199115 7 13975 9781 18169 25 18 \n", + "1546 199114 7 12265 7684 16846 22 14 \n", + "1547 199113 7 9567 6041 13093 17 11 \n", + "1548 199112 7 10864 7331 14397 19 13 \n", + "1549 199111 7 15574 11184 19964 27 19 \n", + "1550 199110 7 16643 11372 21914 29 20 \n", + "1551 199109 7 13741 8780 18702 24 15 \n", + "1552 199108 7 13289 8813 17765 23 15 \n", + "1553 199107 7 12337 8077 16597 22 15 \n", + "1554 199106 7 10877 7013 14741 19 12 \n", + "1555 199105 7 10442 6544 14340 18 11 \n", + "1556 199104 7 7913 4563 11263 14 8 \n", + "1557 199103 7 15387 10484 20290 27 18 \n", + "1558 199102 7 16277 11046 21508 29 20 \n", + "1559 199101 7 15565 10271 20859 27 18 \n", + "1560 199052 7 19375 13295 25455 34 23 \n", + "1561 199051 7 19080 13807 24353 34 25 \n", + "1562 199050 7 11079 6660 15498 20 12 \n", + "1563 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 11 FR France \n", + "1 9 FR France \n", + "2 9 FR France \n", + "3 10 FR France \n", + "4 10 FR France \n", + "5 9 FR France \n", + "6 9 FR France \n", + "7 5 FR France \n", + "8 3 FR France \n", + "9 5 FR France \n", + "10 4 FR France \n", + "11 2 FR France \n", + "12 2 FR France \n", + "13 6 FR France \n", + "14 4 FR France \n", + "15 7 FR France \n", + "16 4 FR France \n", + "17 4 FR France \n", + "18 2 FR France \n", + "19 2 FR France \n", + "20 2 FR France \n", + "21 2 FR France \n", + "22 1 FR France \n", + "23 2 FR France \n", + "24 2 FR France \n", + "25 1 FR France \n", + "26 2 FR France \n", + "27 2 FR France \n", + "28 1 FR France \n", + "29 2 FR France \n", + "... ... ... ... \n", + "1534 42 FR France \n", + "1535 38 FR France \n", + "1536 39 FR France \n", + "1537 29 FR France \n", + "1538 37 FR France \n", + "1539 36 FR France \n", + "1540 45 FR France \n", + "1541 39 FR France \n", + "1542 51 FR France \n", + "1543 32 FR France \n", + "1544 34 FR France \n", + "1545 32 FR France \n", + "1546 30 FR France \n", + "1547 23 FR France \n", + "1548 25 FR France \n", + "1549 35 FR France \n", + "1550 38 FR France \n", + "1551 33 FR France \n", + "1552 31 FR France \n", + "1553 29 FR France \n", + "1554 26 FR France \n", + "1555 25 FR France \n", + "1556 20 FR France \n", + "1557 36 FR France \n", + "1558 38 FR France \n", + "1559 36 FR France \n", + "1560 45 FR France \n", + "1561 43 FR France \n", + "1562 28 FR France \n", + "1563 5 FR France \n", + "\n", + "[1564 rows x 10 columns]" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "raw_data = pd.read_csv(data_local, 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": 6, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020204774983294370238511FRFrance
12020467375219635541639FRFrance
22020457369620165376639FRFrance
320204474391237564077410FRFrance
420204374376250562477410FRFrance
52020427400019796021639FRFrance
62020417396120995823639FRFrance
7202040720786753481315FRFrance
8202039710492371861213FRFrance
9202038722537823724315FRFrance
10202037715844052763204FRFrance
1120203679191001738102FRFrance
12202035782801694102FRFrance
13202034722723714173306FRFrance
14202033712841772391204FRFrance
15202032726506894611417FRFrance
16202031713031002506204FRFrance
1720203071385752695204FRFrance
182020297841101672102FRFrance
19202028772801515102FRFrance
2020202779861491823102FRFrance
21202026769401454102FRFrance
2220202572280597001FRFrance
2320202473880959102FRFrance
24202023755811115102FRFrance
2520202272770633001FRFrance
262020217602361168102FRFrance
272020207824201628102FRFrance
2820201973100753001FRFrance
292020187849981600102FRFrance
.................................
15341991267176081130423912312042FRFrance
15351991257161691070021638281838FRFrance
15361991247161711007122271281739FRFrance
1537199123711947767116223211329FRFrance
1538199122715452995320951271737FRFrance
1539199121714903897520831261636FRFrance
15401991207190531274225364342345FRFrance
15411991197167391124622232291939FRFrance
15421991187213851388228888382551FRFrance
1543199117713462887718047241632FRFrance
15441991167148571006819646261834FRFrance
1545199115713975978118169251832FRFrance
1546199114712265768416846221430FRFrance
154719911379567604113093171123FRFrance
1548199112710864733114397191325FRFrance
15491991117155741118419964271935FRFrance
15501991107166431137221914292038FRFrance
1551199109713741878018702241533FRFrance
1552199108713289881317765231531FRFrance
1553199107712337807716597221529FRFrance
1554199106710877701314741191226FRFrance
1555199105710442654414340181125FRFrance
15561991047791345631126314820FRFrance
15571991037153871048420290271836FRFrance
15581991027162771104621508292038FRFrance
15591991017155651027120859271836FRFrance
15601990527193751329525455342345FRFrance
15611990517190801380724353342543FRFrance
1562199050711079666015498201228FRFrance
15631990497114302610205FRFrance
\n", + "

1564 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202047 7 4983 2943 7023 8 5 \n", + "1 202046 7 3752 1963 5541 6 3 \n", + "2 202045 7 3696 2016 5376 6 3 \n", + "3 202044 7 4391 2375 6407 7 4 \n", + "4 202043 7 4376 2505 6247 7 4 \n", + "5 202042 7 4000 1979 6021 6 3 \n", + "6 202041 7 3961 2099 5823 6 3 \n", + "7 202040 7 2078 675 3481 3 1 \n", + "8 202039 7 1049 237 1861 2 1 \n", + "9 202038 7 2253 782 3724 3 1 \n", + "10 202037 7 1584 405 2763 2 0 \n", + "11 202036 7 919 100 1738 1 0 \n", + "12 202035 7 828 0 1694 1 0 \n", + "13 202034 7 2272 371 4173 3 0 \n", + "14 202033 7 1284 177 2391 2 0 \n", + "15 202032 7 2650 689 4611 4 1 \n", + "16 202031 7 1303 100 2506 2 0 \n", + "17 202030 7 1385 75 2695 2 0 \n", + "18 202029 7 841 10 1672 1 0 \n", + "19 202028 7 728 0 1515 1 0 \n", + "20 202027 7 986 149 1823 1 0 \n", + "21 202026 7 694 0 1454 1 0 \n", + "22 202025 7 228 0 597 0 0 \n", + "23 202024 7 388 0 959 1 0 \n", + "24 202023 7 558 1 1115 1 0 \n", + "25 202022 7 277 0 633 0 0 \n", + "26 202021 7 602 36 1168 1 0 \n", + "27 202020 7 824 20 1628 1 0 \n", + "28 202019 7 310 0 753 0 0 \n", + "29 202018 7 849 98 1600 1 0 \n", + "... ... ... ... ... ... ... ... \n", + "1534 199126 7 17608 11304 23912 31 20 \n", + "1535 199125 7 16169 10700 21638 28 18 \n", + "1536 199124 7 16171 10071 22271 28 17 \n", + "1537 199123 7 11947 7671 16223 21 13 \n", + "1538 199122 7 15452 9953 20951 27 17 \n", + "1539 199121 7 14903 8975 20831 26 16 \n", + "1540 199120 7 19053 12742 25364 34 23 \n", + "1541 199119 7 16739 11246 22232 29 19 \n", + "1542 199118 7 21385 13882 28888 38 25 \n", + "1543 199117 7 13462 8877 18047 24 16 \n", + "1544 199116 7 14857 10068 19646 26 18 \n", + "1545 199115 7 13975 9781 18169 25 18 \n", + "1546 199114 7 12265 7684 16846 22 14 \n", + "1547 199113 7 9567 6041 13093 17 11 \n", + "1548 199112 7 10864 7331 14397 19 13 \n", + "1549 199111 7 15574 11184 19964 27 19 \n", + "1550 199110 7 16643 11372 21914 29 20 \n", + "1551 199109 7 13741 8780 18702 24 15 \n", + "1552 199108 7 13289 8813 17765 23 15 \n", + "1553 199107 7 12337 8077 16597 22 15 \n", + "1554 199106 7 10877 7013 14741 19 12 \n", + "1555 199105 7 10442 6544 14340 18 11 \n", + "1556 199104 7 7913 4563 11263 14 8 \n", + "1557 199103 7 15387 10484 20290 27 18 \n", + "1558 199102 7 16277 11046 21508 29 20 \n", + "1559 199101 7 15565 10271 20859 27 18 \n", + "1560 199052 7 19375 13295 25455 34 23 \n", + "1561 199051 7 19080 13807 24353 34 25 \n", + "1562 199050 7 11079 6660 15498 20 12 \n", + "1563 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 11 FR France \n", + "1 9 FR France \n", + "2 9 FR France \n", + "3 10 FR France \n", + "4 10 FR France \n", + "5 9 FR France \n", + "6 9 FR France \n", + "7 5 FR France \n", + "8 3 FR France \n", + "9 5 FR France \n", + "10 4 FR France \n", + "11 2 FR France \n", + "12 2 FR France \n", + "13 6 FR France \n", + "14 4 FR France \n", + "15 7 FR France \n", + "16 4 FR France \n", + "17 4 FR France \n", + "18 2 FR France \n", + "19 2 FR France \n", + "20 2 FR France \n", + "21 2 FR France \n", + "22 1 FR France \n", + "23 2 FR France \n", + "24 2 FR France \n", + "25 1 FR France \n", + "26 2 FR France \n", + "27 2 FR France \n", + "28 1 FR France \n", + "29 2 FR France \n", + "... ... ... ... \n", + "1534 42 FR France \n", + "1535 38 FR France \n", + "1536 39 FR France \n", + "1537 29 FR France \n", + "1538 37 FR France \n", + "1539 36 FR France \n", + "1540 45 FR France \n", + "1541 39 FR France \n", + "1542 51 FR France \n", + "1543 32 FR France \n", + "1544 34 FR France \n", + "1545 32 FR France \n", + "1546 30 FR France \n", + "1547 23 FR France \n", + "1548 25 FR France \n", + "1549 35 FR France \n", + "1550 38 FR France \n", + "1551 33 FR France \n", + "1552 31 FR France \n", + "1553 29 FR France \n", + "1554 26 FR France \n", + "1555 25 FR France \n", + "1556 20 FR France \n", + "1557 36 FR France \n", + "1558 38 FR France \n", + "1559 36 FR France \n", + "1560 45 FR France \n", + "1561 43 FR France \n", + "1562 28 FR France \n", + "1563 5 FR France \n", + "\n", + "[1564 rows x 10 columns]" + ] + }, + "execution_count": 6, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "data = raw_data.dropna().copy()\n", + "data" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [], + "source": [ + "def convert_week(year_and_week_int):\n", + " year_and_week_str = str(year_and_week_int)\n", + " year = int(year_and_week_str[:4])\n", + " week = int(year_and_week_str[4:])\n", + " w = isoweek.Week(year, week)\n", + " return pd.Period(w.day(0), 'W')\n", + "\n", + "data['period'] = [convert_week(yw) for yw in data['week']]" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [], + "source": [ + "sorted_data = data.set_index('period').sort_index()" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [], + "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": 10, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXmcHkWZ//d533cyMzkmd0JIAgESjnAkQIgRFkQCEq/FA9agq+iqCMKqP10V1luXFdcDZRWUBQVEBQQVVA4h3BASAoQcQEgggdwJua853/r90V3d1dVPHf3OO+/MZPqbTz7Tb3V1VXV3dT313CSEQI4cOXLkyKGi0N0DyJEjR44cPQ85cciRI0eOHCnkxCFHjhw5cqSQE4ccOXLkyJFCThxy5MiRI0cKOXHIkSNHjhwp5MQhR44cOXKkkBOHHDly5MiRQk4ccuTIkSNHCqXuHkClGDFihJgwYUJ3DyNHjhw5ehWeffbZN4UQI131ei1xmDBhAhYsWNDdw8iRI0eOXgUiet2nXi5WypEjR44cKeTEIUeOHDlypJAThxw5cuTIkUJOHHLkyJEjRwo5cciRI0eOHCnkxCFHjhw5cqSQE4ccOXLkyJFCThxy7NfYsbcNf31hXXcPI0eOXode6wSXI4cPPnfr83j0lc04btxgHDx8QHcPJ0eOXoOcc8ixX2Pd9n0AgJb2cjePJEeO3oWcOOToExCiu0eQI0fvQk4ccuzXIOruEeTI0TvhJA5E1EBE84noBSJaSkTfCcu/TURriWhh+P9dyjWXE9EKIlpGRGcr5ScS0eLw3NVEwadLRPVEdFtYPo+IJlT/VnPkyJEjhy98OIcWAGcIIaYAmApgFhHNCM9dJYSYGv6/BwCIaDKA2QCOBjALwDVEVAzrXwvgQgCTwv+zwvJPAtgmhJgI4CoAP+j8reXIkSNHjkrhJA4iwO7wZ1343ybBPQfArUKIFiHESgArAEwnojEAmoQQc4UQAsDNAN6nXHNTeHwHgJmSq8iRI0eOHLWHl86BiIpEtBDAJgAPCCHmhacuJaJFRPRrIhoalo0FsFq5fE1YNjY81ssT1wgh2gHsADCcGceFRLSAiBZs3rzZ6wZz5MiRI0d2eBEHIUSHEGIqgHEIuIBjEIiIDkMgaloP4MdhdW7HLyzltmv0cVwnhJgmhJg2cqQzkVGOHBVDCIG/LVqHto7cBDZH30QmayUhxHYAjwCYJYTYGBKNMoD/AzA9rLYGwHjlsnEA1oXl45jyxDVEVAIwGMDWTHeSI4cFwioJTeP+pRtw6e+fxzUPv9pFI8qRo2fDx1ppJBENCY8bAZwJ4OVQhyDxfgBLwuO7AcwOLZAOQaB4ni+EWA9gFxHNCPUJHwNwl3LNBeHxuQAeCvUSOXJ0CsQypW5s2dMKANiws7maw8mRo9fAJ3zGGAA3hRZHBQC3CyH+RkS/JaKpCMQ/qwB8BgCEEEuJ6HYALwJoB3CJEKIjbOtiADcCaARwb/gfAG4A8FsiWoGAY5hdhXvLkaMKyPcoOfomnMRBCLEIwPFM+Uct11wB4AqmfAGAY5jyZgDnucaSI0etIDmOnH/N0VeRe0jn6DNY+eYezPzxI9iyu8VZVxpS58QhR19FThxy9Blc99hreHXzHty3dIOzbu5kk6OvIycOOXLkyJEjhZw45MhhQVYT2Bw59hfkxCFHn0BW3UEtg7e0tpfxhVufx+tb9tSu0xw5HMiJQ479GslFPjsXUAuF9DOrtuIvC9fh8j8t7vrOcuTwRE4ccuRgUKnzXI4c+wty4pAjR44cOVLIiUOOPocsXEGujs7RV5EThxx9Apl1B7lUKUcfR04ccvQJVGqSmntI5+iryIlDL8LSdTvwpkfohxxpCJFtoZeMQy38HHIClKMnIicOvQjvvvoJvOOqx7p7GH0CeZbaHH0dOXHoZdga5hnIkQ3lSrfnNdjV53QoR09EThxy9Alk9pDummHkyNFrkBOHHH0ClTIAtVAH5DqH2uOexevxysZd3T2MHo2cOOToE1CzzvqIcXJRz/6L1vYyPvu75/DBa5/q7qH0aOTEIUefQE/enOeEqLaQ+qddze3dPJKeDSdxIKIGIppPRC8Q0VIi+k5YPoyIHiCi5eHfoco1lxPRCiJaRkRnK+UnEtHi8NzVFJqEEFE9Ed0Wls8jognVv9UcfRkV66NrIPPpLrHShh3NmHDZ3/G3Reu6ZwA5ejR8OIcWAGcIIaYAmApgFhHNAHAZgDlCiEkA5oS/QUSTAcwGcDSAWQCuIaJi2Na1AC4EMCn8Pyss/ySAbUKIiQCuAvCDKtxbVdHeUcaNT65EW0e5u4eSoyJkW4G7Yzdf6z5fWr8TAPDHBWtq23E3I9fx+MFJHESA3eHPuvC/AHAOgJvC8psAvC88PgfArUKIFiHESgArAEwnojEAmoQQc0WwHbtZu0a2dQeAmdTDDM1vefp1fPuvL+I3T67s7qHkqAA92JI17qubFq2e9aV1PfIETn7w0jkQUZGIFgLYBOABIcQ8AKOFEOsBIPw7Kqw+FsBq5fI1YdnY8FgvT1wjhGgHsAPA8EpuqKuwM5RP5nLK3gW5xyhnNmXd/1fMfJHMYYMXcRBCdAghpgIYh4ALOMZSnfuqhKXcdk2yYaILiWgBES3YvHmza9g5ckQQQvR4cUKtd/Dyeez/ZDCJnj4PegoyWSsJIbYDeASBrmBjKCpC+HdTWG0NgPHKZeMArAvLxzHliWuIqARgMICtTP/XCSGmCSGmjRw5MsvQc/QCbNjRjNb2rtHpZA7KGq6YvXkhaWnvwF0L1zqV6j1Mgtvl6MWvtKbwsVYaSURDwuNGAGcCeBnA3QAuCKtdAOCu8PhuALNDC6RDECie54eip11ENCPUJ3xMu0a2dS6Ah0QtzEQqwEvrd+GqB17p7mHsd2hu68CM78/BV+9c5H3N3S+swxdufd6rbl/UOfz0weX4/K0L8dDLm9jzPfML63r00KWlx8GHcxgD4GEiWgTgGQQ6h78BuBLAWUS0HMBZ4W8IIZYCuB3AiwDuA3CJEKIjbOtiANcjUFK/CuDesPwGAMOJaAWALyK0fOqJePCljfjZnOVobutwV/bAlt0tWPVmnli+NbQCe+DFjd7XfO4Pz+MvC/3MMFX5epZ9cm9cSJZv3IUNO5qxaWcQwdcUj0veWaFvMQ455+CJkquCEGIRgOOZ8i0AZhquuQLAFUz5AgApfYUQohnAeR7j7TFo6yijoa7orujAqf/zMPa2dmDVle+uwqh6L6IQ2VVejGV7S9buyDaebhC1VKvLs8LIvR+aFkh32w3a+PhZ9y3q0Avpfbcg95CuENWSje9trQ4H0tshF+Nqf7c797UBAP77nper3HLPR6kYPNN2g29On10j++yNZ0NOHCpEa+4MV1XEnEN1293XSfFfb15H6orB523mHIK/fUwfncMTOXGoEC1tOXHoDWiu8D3tD+tlMVQmdDicPOS97mxuw/+7bSF2Nrd18ci6F7l/hx9y4lAhcs6ha1DtD1flHCpquxevI6WQOHCcw/a9ramUs9c/vhJ/fn4tbnh8/44CkOsc/OBUSOfgUWvOoTdazWSBvLusnswujB/WiNVb9+HAwQ2ZrqulqKWrdrKSc+B0DlO/+0B0HPt0iMR1+yv27y+pesg5hwrRUePFet2O5pr2V2t0FfE7dVLgLHnutPGOmvsfSg6dgw4pftrPacN+v9GqFnLiUCFqPcF+/I9lNe2v21Dlxxq9pozJfuLhdP177qo4TpFYqcOlc0jGnyrs59QhJw1+yIlDhaj5BNvPZ3S0hlf9RkWifcBP5iwXzN68yYwU0s7wGcFfueEp5OZLOZATh4pR60WjF69RXpDPs9rPtdJ2axlbqau5E9976Dtipe4eQe9AThwqRm1n2H4vJxWJP9VrNmxQTdLklUPas/3fPv06Jlz294RT5Jpte/HMqlTcyJpDcgDuwHvB30istJ9zDrkpqx9y4lAhcs6ha1D18Bnhk/vVY69hTxd4o0td0O6WOM/HP/3gYZz3y7nebWTVOfz1hXV4+rUt7LmTrngwOi5Ei76vzqGPiJX6ysfUSeSmrBWi1vNr/2cc0rqBakC14lxQwW7etcusxjKadSf7738IItFy8bg274p9F+aGBCSrWCk3Zc0B5JxDxcg5h+rC9Dxvefp1LNuwq+J228sxddi4s8VSE3hqxZvYF3IXj75S+2RS1bZaWrouyBHttGTVOIz93lppf/+YqoScOFSIWusA9nedA2NxCgD4+l+W4OyfPlZxu742/qu37sWHr5+Hr965CB1lgVufWc2Ox4RqvB8BgceXb8YvH321020B/mKl9o4yWtvLis6hKt33WOQ6Bz/kxKFCVF9xam9xf5nOQghc9cAreGPL3pr0Z4pIqqOlPeAYlqzbgR37ahtbSOUYPnrDfFx5b3UiyMbmuPbZc//SjTj5yjkol/uIziGHF3LiUCG6yuTSfH7/IA9rtu3Dz+Ysx7/d9EyivKvuz+UAJtGvGOTmaG0vY9tePjmODZ3J/9BVO1nJAfi0/ubuVty2IOCWivs5cVCn2txXt+C+JRu6bzA9GLlCukJU+4MuC4GCRea8n9CGSMSh58PoqtvzFSvJ9bClvZyIm1XLx15tnYMkWC6xUvq6qg6jx0F9Guf/39MAeOV+X0fOOVSKKq8arjVsfyEOpvvoqvtTFdIS3CIsF9CWtg5niGsO1eB8nljxZqfbUFGpI99+b620v3xMXQwncSCi8UT0MBG9RERLiejzYfm3iWgtES0M/79LueZyIlpBRMuI6Gyl/EQiWhyeu5rCrQ0R1RPRbWH5PCKaUP1brS6qPb1cu7v9TYlWq91pm6dYSdKD1o5ygqC41pHuSCfqC925zRfV0jns2NvWI/Oj57TBDz6cQzuALwkhjgIwA8AlRDQ5PHeVEGJq+P8eAAjPzQZwNIBZAK4hIpls+VoAFwKYFP6fFZZ/EsA2IcREAFcB+EHnb61rUe0JlpX139/QVcTPVyEtuYWW9rL2LnrWezE5v3HwVUjrqJYp63t//gRO/9EjVWkrR+3hJA5CiPVCiOfC410AXgIw1nLJOQBuFUK0CCFWAlgBYDoRjQHQJISYK4LZejOA9ynX3BQe3wFgJvXkLRm6Qufg6K9nrVEVw3gbXXR/viIiuYAK4a/ETlyf+YrK8MLq7d51KxYrVenTe2Or3SJNCIH3X/Mk/rG0tgrh/eVb6mpk0jmE4p7jAcwLiy4lokVE9GsiGhqWjQWwWrlsTVg2NjzWyxPXCCHaAewAMDzL2GqNWnMO1U6C092oFeX3UUi/uG5nwulNJShOsVLFI6sMWUQ+cuxZudJaqRz2tXXg+Te243O3Pl+bDnNkgjdxIKKBAO4E8AUhxE4EIqLDAEwFsB7Aj2VV5nJhKbddo4/hQiJaQEQLNm+uvQeriqqHeXDuVnsWdVi9dS/+MP+Niq/X76ar7s6HqL7r6sfxX39/Kfpd60ROWVBJLoo/PrsmI8dRG+rQ1h6Mr65QW7uY/U1/11XweitEVIeAMPxOCPEnABBCbBRCdAghygD+D8D0sPoaAGrarXEA1oXl45jyxDVEVAIwGEAqEI4Q4johxDQhxLSRI0f63WEXodoWD20O2XhP4xz+5VdzcfmfFqO5rTrB7LpqPa7kPW3ZHfs5+F7dg+kJAOCcXzzpXbdWnIPMw76rpR23P7PaUbt66OnvqqfAx1qJANwA4CUhxE+U8jFKtfcDWBIe3w1gdmiBdAgCxfN8IcR6ALuIaEbY5scA3KVcc0F4fC6Ah0QfszdrabcTh+58HF//y2Jc+0gypMP6MG1ppZtM/bKetJv7yp2LvOtGcv0eNH4Jxoq3R6FV2RB95c5FNfFMv+mpVfjI9fPcFXN4OcGdAuCjABYT0cKw7D8BnE9EUxFsrlYB+AwACCGWEtHtAF5EYOl0iRBCbi8vBnAjgEYA94b/gYD4/JaIViDgGGZ37rayY1dzGxas2oa3HznKq361l4JWB+fQnUvPLU8H4qOLTz8sdS4rzcpC5HzrCiGMopBKaKruoOc3iOyXVIIebqeRCfpzrgXH8q27l3Z9J/sJnMRBCPEEeJ3APZZrrgBwBVO+AMAxTHkzgPNcY+lKXPanxfj7ovV49Mun4+DhA9wXVHkxqGhB6sXQF7nOMEZCmDkYdkefRW7vS6D8m2T66MTF1nYra7hWGxF9zrtEp5t3taCjLHDA4IYuHFUOidxDOsSGUEyyaZc9rLNEtcQI0hvVRRx6g5Dt9S17Mtnhq+Burxr33PXPrfO5prNcmmV3vS6c0z0VqTnveBAnXfEgZnx/TtcNKEcCOXEIMaghYKJ2eso9q7Xo9CsGr8AlVuqpTnLqsN72w0cw+7qn7fWN7aTPeCuDKzzX2baT9WrzfmohVKp+UEm+QT20SXfO8XJPs/joAciJQwi5SKvhFi67cxEmXPZ3tn615nFd0Y9zOHSEh6irG1DpouizyGXROXS2jc6ic2IxhjB262agNg6eenG3EoceuvnqTuTEIYScwKrs+laLeV21plK/Uhwq2oamxjoAwNTxQ6rUc3WQ5ZvatLMZtzz9eqfb6YnolM6BKzM02Bs3uCYvdf0eu/PWerJvS3chD9kdIZgcvh6ond3Z7W1tx+PL34w4B5cpq9zZ1HoKz/hvu4w3y3gu/O2zWJjBGatHiJVq8cCZPozity4dSNiHRyeL1+xAe7mM4w8a6qxr3pULz3pdj55u9tsdyIlDCFuKxG17WrG7pT1R1tlp/PW/LMGfnlsb/XbrHGTHtf2ANuy0KzWFENix109Ps3WPOYkOd1v+KTptJ/3a6Ozl1RYDlYVAkRG+1ULc5NPDe3/+BAC/PAi+Ydq7c/Oecw5p5MQhhNy1cIzD8d97IFXW2bmkp8l0iZW6i3PwwZTv/qO7h2BErZ5X56yV0hebdtHScKI3wbTwZtE5bHRsUiQu+PV8nDP1QHzghHHuygoqyeGxvyPXOYSQ89L/I698Mv3l+bVY8Pq2RJmvKWtP2+BUOpzXtDj/3ALpq+w21Vu8ZgdWMvkEslj8uHbq1fBJy8I1yTVswvD+ne84w3g6A33R37CjGeWyyMQ5vMUh3pR49JXN+OLtL2QdYm6txCAnDiHk1PCdI535gL5w28JUWWu7PUaR3Nn0NKuKLMPRF1I1LlO1xUqt7eVI9OGLo8Y0ZarvGkNnrjW198P7lwGoXkIedjxV5reEsu9Z9eYezPj+HFz76KspwqvObSEEJlz2d/zkgVeqOhYTetp31ROQE4cQcTx/391qdSB1HL5+Dj1uDldpUaz2bT308sbM1xQ78TV0RU5xDlJv0x1RNPa1dlQU/0gVK63dvg8A8MTyN9OReUX6+Oo5yzP3BwDnZNwY5DqHNHLiECKrTL9ac6kUhiv2TfbT06ZwTww4FyD76llJkht5RbU9pF3NdSnnYOj8zJ88iinfya5fUgmdOmqbWKmSx6lu7F5YsyPTtbm1Uho5cQghJ4e3KKNKi2Kc59feXjkjZ1MrZBIr6dcqz5B3BKt8DKUKorh1Jqhdp/wcmBtwzYdiF0apM/Usd/1Zwd2LCP+Z6lUyzzujVM45hzRy4hAinqh+k+TS31cne5XcAbrm5v4oE3XtFDtDgCtZPDna4E+gOOJWuYhSOHaypSJhy26/OGDdDeOuXLvxBHEI/2ah1z5Z/0zgFNJCCFz1wCtYVyFR7O3IiUMIOTdqvQbLRcy1kHTX+FzIMhx9Z14tgscRkYqIQ6fGwJRV2bxVVeAvWbsTJ/7Xg5hXYaBDa9+Ogd+3JFvOZ9N7TukcEmPI1IW1Hx9wXMfLG3bhZ3OW49LfP1dxuzpa28t4fUvagq4nos8Th7mvbsFR37gvCrhXa4s2uSC5dQ5SJ9KzqENnxFzqPVfbWqkSsRLbtuN5R8l+uPH79sFU5ObDh341N1W2ZN1Oz16qh4tueTZVtnVPa4JoqDvxxKKtvJa0zkHlHLLPq05xDhbR3r626ikkLv/TYrzth49gZ3PXJzbqLPo8cfjfh5ZjX1sHXt6wC0DtF1/b4qJC7mx6A+fgnZtCGH8Y2/ZFocY6B3b83i/LT+eQVclaS3zqpmdw0S3PRtZUqgyfEysJkf7WXJsFF3TRUBbfBRvX8dL6nfjlo68az2fB48s3Awgsv3o6+jxx0NeDWi++chFzK6ThVa+auOGJlc463HC+fIefE1JSAcm1XbnMviKFNNd2J7gX3zf16ydXOdtbsWkXe21XqKUrmWJvbA3k8u2hSXaHgXMgZcR6P75z2zQvdNFQWwYTJM6SXB3rlfe+7N3W/oKcOGifV6035rHns6e1kqO9rXta8cXbFmLCZX/HkyvetNZ9eYNZJLF66158728vamPlrU50zHlpk2OUAaqmc2DaqZVCWs6fzugc5q/cylybvPjMnzzm11gPgTp8ToQikH7/nJ+Dq20VOnHIYr2Uh89IIycOKc6htpPEd9GPhuWoeNmdi/Cn54OAfjbF4d8WrcOsnz6Ovy9az57n5LdcKAoOvstyQozAnO/Mm7ARB5P4SN8oZEEl02bHvjaWMAD+994V7g6ViVbDa8LxqGKlf/75k9GxOl69l6S1knkMpk2Fbo663TMgpKlN/dlWM8RGTxMPc3ASByIaT0QPE9FLRLSUiD4flg8jogeIaHn4d6hyzeVEtIKIlhHR2Ur5iUS0ODx3NYVfKRHVE9FtYfk8IppQ/VvtmZCTxNvPwdGeGj3WtnC8snE3AGC5QVzB4YwfPwoASRPKDOIUu58D046vSMevWkWwLVIdZRFFra0kNtSP7l+Gf2GUzIA/V9VTxEoSksByO/F9rR3Y3axEN04ppP3GYDql98mFqTGBG6/+/XRG4W1qsyfDh3NoB/AlIcRRAGYAuISIJgO4DMAcIcQkAHPC3wjPzQZwNIBZAK4homLY1rUALgQwKfw/Kyz/JIBtQoiJAK4C8IMq3JsXusq80he+iubYlNVeUd0xd9U83Kso0zrztJJ+DpW31BnCkgDzwGzihrsWxiHXKxmDjTD7rkOdU6J3Hbhd9j/94CF86uYF0W/9nZv8XvZo4fKNnIPWp4krY6/1cMLU05ru73ASByHEeiHEc+HxLgAvARgL4BwAN4XVbgLwvvD4HAC3CiFahBArAawAMJ2IxgBoEkLMFcEKd7N2jWzrDgAzqZtmfa3ZPTnRXYuBL+eghlWoxSPMpEjW2XTXw+7Eu6gSbUikjU3VdzzeXc3t1vM2MZb6DLfvNefB6ApU4xvgFtstaj4PYVdIy/snAL94eIXX+Noc8cls4OasTmyqwTlE/fUwk3QOmXQOobjneADzAIwWQqwHAgICYFRYbSwANb/mmrBsbHislyeuEUK0A9gBYDjT/4VEtICIFmzevDnL0M33pP2uNXHwVUjHgQHt7VXLvt+3FV6c4genn4NvS1V6Z9xib9stDqyvi4fAjOEnDyyz9lewfH1qe1/78xJjva7ROVRwjb7QeyykNuIQW+cBty9Yk6hn+gZWbNrt7NMEjq7om5cOy0bBF5EBQ8+nDf7EgYgGArgTwBeEEDbPG9Yi0FJuuyZZIMR1QohpQohpI0eOdA3ZCymFdFVa9UeHJ0fgG7K7Evt+DrZFx7nh913Tnaasnu10ysfAjrZ2C+fgGION6wiuNz9k9T2/2QvCZKSVy6766ScmDD/0+zdtGlypdm3gxIedMY01wTeWWk+AF3EgojoEhOF3Qog/hcUbQ1ERwr/SfnENgPHK5eMArAvLxzHliWuIqARgMAB/gWEnkOYcqvPSvv6XxZhw2d+d9SKxkuNr8g2fUUyIlZzdG2FXCNoXdRNsXFrVcwhkGEdcnj5jWxASljdMh67H70uAmy2LXjUZh6/MOiLsuzP6n+Ban0B2ej++HtKmT6UzCy7vIZ38XU1z115AG7yslQjADQBeEkL8RDl1N4ALwuMLANyllM8OLZAOQaB4nh+KnnYR0YywzY9p18i2zgXwkOhCm9LzfvkUfv5QECdel8tXq9Nbnn7Dq15sreSq5zeypEI6+9IhhIAQwvqhmRSHcRl/bXp36RClWc/y47GVucAt1gcPM2dccxFf93k/zqHZ5k3bwxTSctSuzY4QwKubk6bRvh7Spm8h6zv/n3OPi465hV+fn+1VECtJrHxzT811SVnhwzmcAuCjAM4gooXh/3cBuBLAWUS0HMBZ4W8IIZYCuB3AiwDuA3CJEELO7osBXI9ASf0qgHvD8hsADCeiFQC+iNDyqSsghMAzq7bhR//gM0x1V0hs1845tmryH59p3Zj32paIOOr49M0LcMjl91iJlXrqdcb3wXeI7thKfg09+FL2xD4cuOc1oN6cs9nm7SvPz75uLiZ/8z7D9Waoz2Zfm5k4VDPlaVQvQ5vpPoK/PrvsH9yX9DpWCYrtatO5rBv7+lK8/HGcTlcopOX7+tiv5+MdV/Vsx0ZntnIhxBMwz8GZhmuuAHAFU74AwDFMeTOA81xjqQZ2OixIXN/PL//1BFx0iz1K4yplwRRCeFkNufqVE9M1PU2JVaLzZYEPXfe08foHQ+9mO+cQn1u3I5343XRlV4nw/rpoPc6bNl4rjdt+75QD8dcX1sEFAuHYsYOxeG0cw8g6woRDV7rmQ8s2YfMus77APi3i9mzEIQuEsPcZEbtOvBb5HCrJj6BeYZsbpnDmWedTnZL6j82poROHTlhDcdhkmRs9AX3OQ1p/4WnHLDumjB/ClqtmdKf/6BGl3G/CukQskqV1+0MoxIFZCFzpSH3HI3HvYsbDugLOgYPvp/7YK5sx4bK/440te+NrlYuHD+iHH583xautlGe15xaWe1w2wgDYs7mpz8YmVsoiVdrahWIMubj6+uNw4PI5sH2ZxJYZu+ynEAfus9AJnGm+bt/bihn/PQeLKwiM+Nrmyi2suhp9jjjo7zf12zHBhvbvh+MPGoJTJ42Iyp59fSsmfe1ePMXEMvK1vXb1G4mVHEum2g7HsfhadNgMM9QRzHk5HUepVDSEp9D1O1WyVpL49l+XsmMM+nZfT5Q2BdaJ5NY9rTj2W/dj4ertScV8tqEG/VnOqf3axUr+1OF8C8cIAHVFGSeq8wppV+gKVlflqXMwLdJZx91PFSuLvbDDAAAgAElEQVR5WCuZRGVPv7YFG3Y24+cPu/Nd699Ad4Rc90XfIw667bL2wn3i9xeIEpP3sVcCovA0k3jFV4nl5BzC1ToT58Cc9w2n7auQ5qB+dPY+lDaV5/70a1vw4IsbM3/sqkxYHWNLu79YRucc9Hud99oW7Gppxy8feVVbzLIvqL7WStVyvlru8AMoRYmnsrctL5HXVpJS1Nd6zfSssz4mdZ5y831PS3LemL4JueBX4mcxsL7ortRN6HvEQfttiwzJoUiEAiWvk7vx+rr0i/YV47jmdeznYK/3oCMiqi8n46tz4FDvSRzU9Itqk7OvezoIs5DxYzettTv3tXtyDpTienxNKl1K8d0t7Xj31Y9j6TpV9OBnrWRDNY2VSsXqLQeuTYgzb3YFnENWU1YX57Bqi25Nxbcvzcdf3bwHa7btZeuY0FjnVPt2G/oecdDeb4pzcDmZEYGIkqaGIdvfyBAH33gsnM33O656NIrfE+8ezePbtFNTDnM6B2/OwTxOF4Ea1FDHluvD+dyt9jzcWXeC6qZfaMTbR/xCSOsB9OlgUpr+4mF7Mpgnlm/G0nU78T/3xV7TNn9F33WumoasUgbfGT7lR/9Yhq/c8QKueoC3BrShO3UO3MKvrw2m+ahOmd6QxMcXfY84aBMrLVayo1CQnENcJolDA0McbB62KnQa0lEWeGXjbvy/MLKkj0J6+n/PSfzmFkRfTsYkXxXCvUM7bOQAtlwXa6giN67FrDtBVZ6rXpklYJquc1DbWbJ2R5QRjCjbItrMpJr0FSsdOoJ/nq42skJyTabHPnZIo/Faec1dC9fh9gVrWCs2Fc+9sT3dBtOera90efKEi4Otd3AOro2khLqhcOUR0d9XT84j0eeIg/5FpzkH86XyvQc6h7iiXHDrGEWsr0nfbQtWJ37rfEKkc/BqLQC3cLR45sM1OhrBP2udCypb7xQzeCDJOcTHbR3lxLP4xl18rCIi4F3HjkmUqeN6z/8+gUWKRYp6zqVn4ZTKvuEzRg6qj+ZWXZG8xHafnznJWUeHS6zU2K9r5eO+HtJG4qD9ds2fpM6Ba0/nHEw6h/g462LfkyO99jnioL86ffG2iZXkrqBAlJxM4TFnmphlgVNzMaSCflXgBMdhyx4/22qbXNc1BN8hqoscd0n2W+UXW51722tg/QnAuSeOS5RJcd6StWkzxSVrY0uTkQPrrSNrZoiDNfCediznVoEocc5EYIiAaz9yAnvOZD5ZJxXShoW5q+MBqeuk3VrJMD6HdZHup6Du8jmPbr3I5PWttpPVeCDnHHoQ9Hmlv3Dbq5JiC9IU0vKY26lnyR6VCOSmXebrBKe3t7ulHVfe+3Kka3hzt5+te2fESr4EbOQg+4Ka1eIlwTkoT8o3YBqF+iQV8jnoOQWIgJ8roaQ//JaDrG2zYiUL56AnxZELUKmQlGfZdtjvPHYMjjxgUKp8MUPozjhyFKYeFPjwcK/vL8+vxdpt2S2QssBf5+BXXhbJuaiHUFfvk+Xw9Q2aYV6rm0KXdeIa7RlWMwx4tdH3iIOuc8hgrSQX2AIRu+gTpYlBlnevVo3yN4SFvkmB9PFcPWc5fvnoq7h9wWrcOv8N3Dx3ld9YjGIlt0Lad4f5vqljo2Puksv/tNirHYkE56a019perji3hbTuct2zS9bMcQ42bfK/3jAvOhYQUfvFAiWuMypJITcy6U64Z33lB4413kNLewe+cNtCq49MNbzdH14Wh+G3ekhnMGVVy7hc1hKsziHVL3+t+tRsYqIPMVn/cs6hByGtZEr+9lnY6oqU8HyWR7969DX848Vk3uYsrLhIcCNx+a7mtphzyNBeXbEQLUrtHWVc9qfFCZm5DaZdkhfn4D1C+1XZFdJ8a20dZS+rHr0OUbwT5JzhVLjG2sz4WviSq7LKORQLici7xgUr0o+lz3Fh3U3KfFsfiTruKk78YX4crLIaCmkgufjq3BvH/dv6MS3kajFXp62jjD0t7ZjHZKbLOYcehBTrqb2c6x57zdlGfamYcKySk+jlDbtScZey7Aw4zgEAjv32PyJ5aZb1sq5YiOqbds4m0YZp0ra0lVn5uwrfW05m/vK7xgZTOIr2DuHtIa3i1Ekj0RbejP48nn4t+aG7xs8ZAqjjPfrAJqM1kBAiIgjFAiVl5R6m1zqKTFnAkPAPSe3igyeMY+vYUK0EVN98z+TUeFRw5UlRlVlKwHMOfmIltXwlE4jykt89h6O/dT97bbXjNVUTfY84KC/yxO89gGUbk3l8fYJh1dcVErsQq3w0i1hJ8MdAZToHznpKx1UP8vbopqxXl/95Eb5511L2HBCaeHretCr75kUCGSmGwVrJ13xX38vXFQgdZSlW6pyehRMrqWv0ceOG4MLTDuXbRhwkTgizP0ei7fAvty5zRklJkZy2iCqzblBDCX/+7MlsnxxOnTQCx4wd7F1fgnveQ/rXGc+ZytVFP6VvNHDqpvpGcZZy8ZfvWIRlG5Jryj9eNDtIVjMMeLXRB4lDfLxlj105O34Yv5NrqCsm5K+2hSFTdMoEcdB2LZWE7PYUXNz57JpUmYlzeO71tH26irpCwZsg/m1RHLTPh613QV3gvve3F6PjQKzk4QSnVSkWKBYreeQnsIG1VtI6NHE3QgDDBvQDEIgYk5yDvV+XCEkdiyxOcddKQalAGMiFMTeKt6gikRP3POXzMrXHlXMhVb7zz0fj7ktPSfglce9XL/m6IV2rPnfP/ql/KO6enBGuzxGHLJgwnHc+qi8V0OIZRrnSl6/PVTnJ97R2JMJOuPqWuz6bWOXaR9PevZXaXxcLVFHgNu4x6ZYdLqjroMoR+oqVdNQVC5FC2iUbVk8f0NSQOs/5OST9Fcyh3QWA4QMD4jBtwlAvsVKsc2CIA1e/YNaBqH2UioVMyn0CKpIZcle4Umyy2dyUF/Orx4J5fsDgBhw3bggOHNKIn35oKgB+E6cXmRz7XKLjwY18xACgOrqarkKfIw7VINQNdcWEgtEuVsqiczAryNSwF5zM/4XV6R19YMoXHGddGyu1oigVyRrRVUe5Ao7IBNM9+oqV9OtLRfJ+Di6CKA0Y1HU14echeBFQcC6wVnrwi6fhVx+dhskHxmIac0iH0FqJHWsaBYuSW82fUFcka9iP9DgqWwC5+SDvqRKdw6I123HXwiCnhzr8WcccAMDkIe03ctcGcGh/M3HIOYcehCy7WtN7KxJpskyLWCnDQqk2Y5s0nKfqOb94Mt2e571y4zfloXAF7itl5Bzkjq0aRht/Wcgn9PG2VmLESv75OOznbYtdcN6sUBciWNAmjhqEgfUl/OLDx7PtPrwsHXSRa5MX2aiLr1nnUCpk5xwqWf9YzkGeM1rSmUVDqghZHb/kwnzESia4uMomC+fQg42V+iBxcLwMVZ5qWqALGXZDXBsm6w215q3zV7N1AD7AH9ue5yC5anoKR4ldLfZMeqViIdOEl7u5rtxBeVsrKSRkYH0JpQL5+5c4KvAKT2WDAWEkYEITOamBDdVuP/GbZ1LXcvf98oZ0DoECmbUyCZ2DgXMw3X2gc6iOmDEWK7mvmXbw0LCuSA1QHb+03OLFSn7jvvxOuz+O1QemN3MORPRrItpEREuUsm8T0Votp7Q8dzkRrSCiZUR0tlJ+IhEtDs9dTeFsJ6J6IrotLJ9HRBOqe4tJ2F7FiIH1kUUEYHlvWj4H2/vldiR1hhg26mT87dOvG9s0Xc/1LVv8hsXCiHsopixmLtO7uiyUE8BNT60yDaFqOP2IkcgiWLvxEyfh/v93WhiqQo7Mf/HnFuSFjNhPb9HFOfD92sVy+sJ018K1+OmD9qQ0eksq4a4vFQzciEH3gUrXP/OmShcBrdu+Dzv2tUXv4IVvvgMfCE1ubUQGiBX2LOfgMe4tu1ucGyYbejvncCOAWUz5VUKIqeH/ewCAiCYDmA3g6PCaa4hIbnOvBXAhgEnhf9nmJwFsE0JMBHAVgB9UeC9esO0GSoV40f/xP5ZhLpO8B4h3HrIt7sOX0COR7mxuM2b2UkfGTVYZKMx3PlVz3s06OpDNuiZzwDlk17NUwjlMHDXQWefnHz4eP/8wH2NIh1w0Tj9iFMYOaQxDs4fjdImNHE97x760d66+wTDtMIUw+6nI8Znei76Qv6yZWar1ImslrS313fTvly3/QGDanOkSdgxAINIC0sYSJ1/5EM78yaPRO2hqLKWU1+r70XkkPcpyNAZmXM++nvRvMZl1q9+vbVvSq3UOQojHAKRd+3icA+BWIUSLEGIlgBUAphPRGABNQoi5IlhRbwbwPuWam8LjOwDMJNOXUAXYFjf14/zfh4K4OYMa0h+DnFxlEciz11vCE3/r7uTk4VKJSiSccphJc/SBTWG/1VGUZYEp9WeqXoE3XTSZgsrvvBKFNOfMpWPc0P5oqCtW5ARXUHw2XKNTb882H1Qk4j91lPGeKWNw/EHpHOVlYRaLuTgH/TrTe/B1guvfr8iax5phNmW1EXfuGjkHOT3Q5l0t0TsgImUDxw4pgWKBvKyVAOCD1yZDYLy5m+ewVytJf2zLWW/nHEy4lIgWhWKnoWHZWACqsHxNWDY2PNbLE9cIIdoB7AAwvBPjcsD8NgqF9Ef2nuPG4MnLzkjWUzgHU4RPE+TEvu6jJzIji/vmPuKShQXmUBbA3QYlbbJfN1yxgyRKRWKJkklpJ0srsZzNsoXwUkindpQZOIfw/JMW4m+6BghCO9SXirjyA8cZxmZqgydestwUSFBHwlpJO5fkHIrZrZUMD++YA5vw3ikHsue4S2RyHpNRhFCIqLxvPUYZkOamTLHSOG5Q1/fpwfwkOjw5h2pY6XUVKiUO1wI4DMBUAOsB/DgsN1nO2SzqfK3tQEQXEtECIlqwefNmrooTtndRZBx2mhrrUnbKqmKMDahmgZysh3G7JoceoxBNeM/OhEiEATdXczfos0sHgGKBV0ibuBjfnXmlyMSCMpwDJ5bgIO9jPhM/x3iNcixNozkiHIiVTP0m/0rc+kywR9Obs2YzM3In8XFjvyLLYZiejnVhZMYXn2M2RyFxMHkVqxZfkcMcp3PQfhcLvMnyA4xnsy5J2NXCB/NLEIfeqY+ujDgIITYKITqEEGUA/wdgenhqDYDxStVxANaF5eOY8sQ1RFQCMBgGMZYQ4johxDQhxLSRI0dWMnTrJ14oUOpl9VMCnckPNzL5Q3bOQbZvUj5KcGxuZHbnLVbyG5OPuaavKKGuyAuZXf4CXS17rURSScqO0mmqGv61RS5NXZPgHIJ5xIU8CSyZ7DoHfTGVDoRTxg3R6pvET6pZbbKOuqsuFQoVcA6W84ZyVucgxUoGNvPx5Zuj9uRf7n5TJsuUFitt3NnMOmHqZuRGzqEbRL/VRkXEIdQhSLwfgLRkuhvA7NAC6RAEiuf5Qoj1AHYR0YxQn/AxAHcp11wQHp8L4CHRhbyWm3MQCba1rliIkrIcPjqIja8q7va2ZrNUkJPBZQ7ITRqTTbY5hLHfY/TxhvZdEEoFYhdS08fikpl3Fjb2VIdep6BYpbnGJ99JNhoUt3nwsMAbn8vGJoQ5MRAnNlHxuZmTcP70eL9muw9u7Ft2tyQ86EvFdM4LGwh2U1abhZYOKVbiTHYB4IU1O2LOoWBuJyVWKqTFSqaMiToHbaqXFCuZn1fPJQ2A0/SAiP4A4HQAI4hoDYBvATidiKYiuLdVAD4DAEKIpUR0O4AXAbQDuEQIIbfWFyOwfGoEcG/4HwBuAPBbIlqBgGOYXY0bM8E2UYsh5/Cl21+IyuqKBdSXivjNx0/CceMCz1SVZfV1kpKQc8blnMQtsJHOQTvnE0rYBp/gX0Vb2jIFpWIBrR1pgmkK5CfH+OQK3jLMBp9FKlMsqpQs2p/AylpZdtVCAKMG1eOqD03FCQcFajvOB6YszJyDcBCHYoEweUxT9Nu2o+VEMV/64wt4RMmzUCCDn4NFXGXlHIwhQzixko9VQfAnFsGmn096E5D+VtRhPXnZGTjv2qewbkcz6yjJIfFNWobdkzkHJ3EQQpzPFN9gqX8FgCuY8gUAjmHKmwGc5xpHtWDlHEJLm7tfiJW4ks1/+5GjojKVZc36cm1Z41wKablA6x+4MfeC577E5fUc9B38HTWo3hq5ts4QPsPFOdj8OjqDiHPwX1ciFAr+CumYI0y2ouf+SIwt1CWcMnFEVMYRB8ENLupX1kn28S/TFCmuMqZnVm4z3AEvitHzVpQKGTkHsu+OTS47NlNWa39Rv2b9nF6kxtCSUMWo0qwZ4JTZ/Dj8FdKWk92M3ENaQZHTOTDJ3NXokK6dqW6NIQyLiD42XqzEnzNJhXwnnk/CEclOyx2usV6hwC4Gr27i8xZ3hq/2WaLkM6jEOJoyKaT5MdkWNE6XwImVIHyc4JLl50+P05aq1+oh6lVwUU9bNR1KscBzDg11/H0SyO5b5OHQecrE4Xj8K29P6GNamORJQByvStb0MXhoaqxLZYnTb9FE/E2cg28Sn/3RWqnXwvaRB7uDpE05542sOti45kCHtnJbxUpMPRVyIqbCeZt25Z4T1Gciy52Uq25dgV8MPnTd02x9HwILpBeD5raOjKas2alDUudgr2syHdXXDj2/gH4PJs7B5QSnD89nl62Di3qa2lETzzkMH2jIB+7gHOq0+5XvWZ1nx4wdjPHD+ie+GVPEXhmG2xXeW8XgxrqUk6K+ATNtMkzEQX1utnnKfU4bdzbjviUb0idqjL5HHGzyT+Y8Txxi2axrAdbl+VaFdHjOpEOIxEqeqU199yQ+2aiKEZset/q+qWkbdY77skEI4RX59KX18Y73ry+sw5HfuA/LNxq4kWQP3mPhnOBcCl8g2K2axFf6QtqupZflzCp1CCGcJp86gVXb8SWinIe0LhILFNLpa5sYZ1EgvD/LszNxDjrHAiTvyfTdSWuigkbobBuQgfUl7G5Jbj706vI5cz4SHNT3rG9MLnrbYdEx9+1+6FdzcdEtz3Z7fuk+RxxskPJR9VVypoUqy+p6f/pOW/XilPjUPx0Sthf8NvlOxAppXaxkl+dzUOe0zxyUwd7U+/mMMskligXeCc4EAb/lu79iQvjgS4H9uU8o7mi4lYiVQIpYyYwB9SWj2EHvVrUMCziHZI0GJqhimRErPfjF06I2uPEliIPnzct66kKqL9KBQjrdHsepzPvPmdZkPxt2NKeUzLLrBKcYlqn31GFY9BtKwfOLdA7h8NU5rk9PYnaFqW9McvzabcokTDo6wzms2rI3HGdOHGoKH/GA+sGeHcYUUqG65rsS7+hseaxziMv6a5m1TLbysTWJJlYyeR87uKQsGDs0yIq3TzHd5UQgBWYxsE3yssY5vOvY9PMGYjPGSuFzv3dp3uSqFYvtHgqKRY6rn9YE58C3+dnTD0s8Wz0qKxCE7x7Qrxjnw9CmjO/jmnX0AVjw9TMBJOe1hD5/SwadA8fx1JcKITfO3+e8lVtTc0hW5TgHlSjJnbl09pNo6CeJQ/Dbh3MIhMlJGLl37T1Ipzjdc9pm5CEgMP8/Z4Y/bN+G8VRN0PeIg8PmWt2lHTJiALuTU13zv/THF1LnVehOcvJjVie6/iG4OAd93rksgTi4nPB0SA7qmVWxtQsrEiDzrovDyYeOSJoZGk0bK0NnNl9qfgN7M6TsLO3PNSHCY3QOQJrACoZzkPVM1krq+7W96lKRMCLUF3AWPjp3FiikmU0BNxWI3NZKWlty7iTS8Cp96/Ve0RTsDaFCWh+jzYSXKC0K1edsrHNIE7Njxjbh6ctnJspVcRz3/Ec1NYRctnFY3W7m2veIg1PnoMRnMdSLdlge/T37+rYEi6wqpC9/55EAgBFhCkiXWOngEf3DNnSxEt+33b7cNXK9PidK4DkH/cHYdlGjm+oTH64p14V6z1l1GkBlHtIJu39LnwVFsJ7qRvutPgvB1UfSSirq31FPfyZJsZIZ6g45ntdxma4zMxE/bg4S8Xo8W3vyflidg8o5hOPWd/j9DToHm/w+4BzsYiWO45f1hjT2w2At21tCrGR4A8T0k+zTeKom6HvEweO8U0bLKGdt+O3c16PjyM+hAFx42qFY9l+zMLR/SBzC0am7tUe/fDpe/t4s/PGit0YiLr3fSjiHrNY7XG1WWc/0+/17XjK2WxbJ+qZRVSp/lVeZiI6Kr73rqMRvdYGxW7nFi6P+XFM6B2WxXbJ2B9ZvT0dw5RZUNhc0xZZh+ugG1vuF1lZ1SBznoIeXLxk4B26uScumLB7Ssqa6oYoWZmW6SQ5c1+nFCunkvSSIrTYezlEvRRws45VFKkFWicMTWjDGSaMGRW3ZZnXOOdQYzkXGsEtTEc0BR1NTQo9qdRFVFZdEhPpSHE5aDk1dQA4eHoi2TpowjLUYAtIK6cNHDwzrWQaXcSPNKiEZZT2X68eUF6MujOAqkpIWFpXKX+Wj8okqq9eJEsEIx+43XAAfWbYJD7+cTtWpQlVgLt+0mzcNVp716q17sXb7PoNYKX5e6ryePmGY2bRUg26ppipnn3o1HWHW6CHNtO0zxfS2RDnImcDFLFPfj3xu+txvrCslOo85B9so0mIlndMw+SeVRUxUVc7GFjnhgycEAal1DlFHThxqDOt6qS1sprpqPgeJ9009MKWU+s45gUP46KaGqExeo06kiDiEv11hlXUWXq9/9fnHY9SgeqMV0ykThyc+XBkWZJvmDavCJBvnyvRJbQpOKJPpqJyP6XtI7vzMMOU/9uEc0uOL+5atHjJiQLoegvf68d88g/mr7FFZp45P52vQoYaEv+rBVwAAK9/cw9Tjrank+9Tvg4NOnOS9AMD2vemIo+3lMiui4zZdMvVoljVu0drt+OC1c/H9e+I0tfJ6db6ZOYekzkFeazM55/QiaVPWuK5ekZhzJlHq4Ma66PlxHIuKXCFdY9h1DgGb7lpGONnsp049FD86b0qiXmwGG9fjwmfoJoQmRzN5jSt8BiFg/U3RKy9622GJj0HqOL58xyK2ftBmGiwRozTx2mcgDjKZjs/C75vvwfR+faLKpv0c4gUmXqC468y+HfpC6vPBq5uPF9cF+Z650Otqpjr1Gfrm3gDS71BNjcrd0+DGOvYZGBkgywL4hTMnpco27gxCs6ii1aEht8VxDk0NSVl/Y11S5xD5DrmslRwWgJK46M9WKH2ZxEoq1MtNeSSitnPOodZIP3C5Gxw9uCEhQ3Rl11LfK8duFxWxhMSNT66K6uvtuTgHk4f0S+uTCeP3hd7DOmt72Mh41/vBE+LYO9IyRA8hoILbLQ7RlHAAz02YOIdiuPNNeOQazHhdLLb0FUnVC3/6cA56DU7nwNnzE5nnyshB9cl6HmYM6vyT75ybE2qmOrVZvaZNv6RvRFQLKG5c/fuVvDmHQCFt7vvQkXFOE2kSmtJXfeBYXHjaoQCS3LZc7IcPTPoZNIZpTHWdg1UhzRK7ZP1pE4YBQKQfVOtxcZdM3elrhm025JxDjcF9w5+bORHPf+MsHNBUH5oN2hcS1cRRolBIL6AxRxCXyaB1tsiWphDapmQ/l/7++cTvshAoEKXkyWrSos/NnBgdS85BHdOE4f2T3A0zXl+FNHc/D37xtGghUk+bYvTIJnfsbcOzjOhG7i5TJojhX9Nu2pa2VY3sKfv/PLPbtX3kv/3k9Oi4aOEwVKj8ZslKHOIkNepZ37ApA/oV8dVZR6Q6NxHi2z/zVmNbPBNJVsKZEK2Gf9W6/YoFnD/9oGieqdyfjPKrty05B1nTz8+BM2VNFvzw3CBDn67oD9aLAOoUM/sexeW5zqGHgXvcdcUChg7oF1lW6Dt5HfFEVsvSYQU48VN8TuUckheadQ728xLHjx/Ccg5qn/3r4kneHMakV4nij86bgme/flY8RmuPyT7U5zLvtS2pcYwY2A8TRw2KPg71I2g0JLCXdS74zXysY3I0q/oBFfKnKdbQf/55sdKGRtwVYixbnTYhHXjQ9pGPGdwYHRcMcad0FCIOMT7mRITFAkU7f2uzhpf30H+cjhMPHpbsm2Cc+DbfGO6SAsUOZpxosUDxM1bvWYILeikhOQf9U2gMNxdqiBvAzTmkTVmTdYYPrMeYwQ2J97x1TyuWrtsZicDUDYgwECX1l0sfkxOHGoN73nJ3Jl8Wt/ir4OLeC4hUfV38pO7o+O8sOC8/+M/NTO5SC5pY6ZWNu9hkQ0GCdbKGIR7cvw4Pfelt+PjJEyKzQXX9LBQoYVnj6zSnLpQ3PLGSDbgXye/DOEz//oeY89mwg/c4l20uWbuDPS/Hd+dzaxLl8qM3cQ7qMzESd0XpwHGVOkE0oWgR2XAoCxFzDoz1S12RIss2dXFLi5V4cOUyZIgQAr94eEXinE00x4uVJOcAPMRYcanPnuOK623EQVor6ZxDP5POQR2rPk7GlJV5UTrdvPiWZwEAjy8PuM9EeA8mbIfet2uz0M20oS8Sh/QTj3aVJO2Ww0XYsIXidqlCxBN26vgh+M3HTwKQXMxVr091l6oTI7kQnHXU6ES/8gPqKAu0tHfgHVc9hkt+95xxjLo4R/dGPXTkQDQ1lNDcVg4V8WSs60kbErF0vve3F611pbXNs6/HXtfrGa4AcMtf5Xf5tT8vSZTHnAN/A7bc2AWGc+Cqc6I0DqVC2ub/rYcOT7en9FGMOAdm3hYLsVhJOd2ZHacMBfL0a1vxoqbLsim6eT+H+B1zClrVik8XAwEOzsHAMclrVDNkwGGthLRYkFNgq4YHm3Y1Y56WLzypc+CJl7r+lAoFa3ywnHOoMbjHXSxKziGgDt6cg9quiInJ6KZ6vP3IUSm9gsnzWRdjSc4hZXevcCJybE++yvsQBJyDJlZi3raM9CpEclHyjT6pwyZjlpBnC5Rm91UFbuIaR5uc/iPZl4E4WBa8hEI6bIjdbZPBzl/nRLRwCUP712HS6IHQoeqq5Pi4yLmlQswdLloTc1TqohuMw3CPHKELuRtu0UhuJjgAACAASURBVLKKlZgHEHAOZm5JNesl5puycQ6yXjoabTKfQ+TnYJs/jFjQGHUg7JnLHc2F97CFkmlqKBlzUOt1uwN9jzjYxEpS9uhYB1nOAQJHHhB4Po4f2j+sR4l6zYYEJbriWi6YupNZvBuKF6u0E1AceCzl4MSKRNL3AjCLpuWZyPuW7flueLiF42ezj8eNnzgJnz09GfHV9aGoPia/n/dGdCwdAk2Lg7rgpcUN6vMWiTKJujCnMitW0X6XtHDmZWHyFQnHAxGd58ybS8VY53BRKOJ42+EjI8stF7j5IL8BbtdqS9Npej/Fgp9JJmchyHEO0srO5BleF33LMbEpl4XRYg4I35PWEKfjIUWuVMfstL57Tpzo0hTWROUcB7mIQzdTB58c0r8G8B4Am4QQx4RlwwDcBmACghzS/yKE2BaeuxzAJwF0APicEOL+sPxExDmk7wHweSGEIKJ6ADcDOBHAFgAfEkKsqtodauBERVK04KtzkNA5gYmjBuGuS07BEeFiqbcjlXI/+ZekPwQpiwEQL2RpzkFZrDT9BBA44v3H2UdEffvkt1bZ7407Y5HOUWOCe7jyA8dix742I22Y/7WZGFRfhxufWoX6UgEr39zjZIdlWwlTzBDDBvTD6UeMwuqtexPlrjbVAImqklkqhE3X2zmHmGhHYiWtzvHjh2L7vlYvgqjLmAMzyHQ9VVclg+J98z2TU/VKhXR6y3ccPdoru5raT6IM4Xxl7sfOOZiJb4cyX3UcPDzYSB0yYgA272pJLIgccbj0jIm487k1xp25vHd10/ONu5bgd8qGQR8JF1aci+2kcojco3jLobFyP9Y56GKl+HhQQx3rv8LV7Q74zKIbAczSyi4DMEcIMQnAnPA3iGgygNkAjg6vuYaI5Fd7LYALAUwK/8s2PwlgmxBiIoCrAPyg0pvxAvPA1RcuYGHDQ8iPZIfiQSpf5JTxQ9LZqMJz0ipIzU0AMDqHcNeiy8nj8Bn8Tu3kw0ZgXMi1cAppG8pC4JUwec7Jhw2PnsHs6QfhM287zLgwjBrUgMZ+RVx8+mH4t386hA2foUMV9fjKVbMQBw6mWEOZxUpa9es/Pi1S4urQ55FuyioMnIPqFNnYr4gRA+tx3rTxqXqlQmDKqi6ob+5Ke7kPrOefDXfnkoCxnINVIc2Xuxy93n3sGPzxorfiQ+H9JcVK6XHrYcX1piV3o4a3VwkDB0KauHHfTkHhEE2mxRIx8UrWUbvpVyqwREhvo7vgJA5CiMcA6Ibl5wC4KTy+CcD7lPJbhRAtQoiVAFYAmE5EYwA0CSHmiuDp3qxdI9u6A8BMqiSEpie4xy1fdOQh7SlWsjmNqfXkS96+N/hw67WFjLQJLy1Q9IWLEotV+k7U3dnyTbvx8gZzvmAJbnGypUZ1gRyLgd63qWpTY9LBzvadnHzYcKN/hMTopgb84IPHpsqt9u+qjkeWgTB9QrxDbGqoM4ZB+OiMgwEA139sGt5+xMiUyWvgj2LuV4TjM63J67bvw1OvbsHSdbHimBOfvWMynyOD+8wI5s2HjXMwLWTcO77vC6fins+dGo3hpAnDWIdRLoeHbtWkdyvFPdG789gfcTojlnNQ+uXEfOrTkd+A/lzUR1hXtG/gVm1Jh0ypJSrVOYwWQqwHgPDvqLB8LAA1+8aasGxseKyXJ64RQrQD2AEgbcJRJXBzWH5Q+iSxscpA0vPXJo+VrXz4+nkA0olBIqsmKVYyKqTDD6PMZ6CrRETJ6Ry47HeZFNKefdv8A/QkS7Ke/lEe0NSAGy44iRVBfOD4sYnfx45NxzXaoIjS9PedyOegaKRvv+itqXr6s3/lv94ZiYLOnDwav/nE9JSznHRWNEGUg8XNVEf6e/zjxTjfMBvjqECYPKYpVc5yDkShWXa6HU4OXy4L/H3ReizftBtEwJ0Xn5w4XyyEvizKpUce0ITJBybHw+kc7DoOKdPXxUoa52BsQekbyXVhT0t7pAtQI/Wq4ieXr1GUqU57ZH/+7CnRcV3Rbq30hdsWugffhfCL6+sP7m0KS7ntmnTjRBciEE3hoIMOqmR8rOyzHHEOms7B0IZcUNXdxRGjB6XrGRI/6CIQ/ds3WSsldleceEwpe9exB+Cexe4k5Zx9ucnz2Qcum//phwzDpW+fyNb9zj8fHR3rz8j0LR40vD8a+xVZk9QJWpA87t3v3NeunE9CPpuW9nIkmuDWaTWfgwRHrNKcA797lwlsnl+9zchdqBik5G82EVubbiM9Rv55q9GCp4wfghdWb0eHEPjOX5cCAI4+sAknHpx0EpRe3K7FVDfeAGL/gUR7yjfV3lFOp+GV3zKz6ZHQv1U9rPjR37o/On7/CfEGQxU/cfejPs9Y7JWspxLFfsWkzkgIge8q5t/drZCulHPYGIqKEP6VHi5rAKjC0XEA1oXl45jyxDVEVAIwGGkxFgBACHGdEGKaEGLayJEjKxo49+1E7yB8u9IczqzEDM6rVN/EogPA0nU78JU74oxxIwbyeWdT1kqaRYS6y+fGpkYMHdyYjnskiYtqpcJ9RKYMbz4gmJ8bEIRgOO3w4N0VtMVy9vS0XF3CLLYI/nK6A58hzzqGF7mobV//+GtRVFSuzZb2Ml7b7BYB6OInk8hopuLfYiIgAPDVWUGyKNXwwLQR5ZXPzDOLCHb6eY8fFnt7v2Py6LC/eOGfMi7NmUlHR6spKdK6BBPkiMtCYOLX7sUNT6xMnJfcvCml7qor340DBmumvpZ+E0Re4Yo54qC2YXLSU1FXLKCtPT6/YWczfhPGXgPcus+uRqXE4W4AF4THFwC4SymfTUT1RHQIAsXz/FD0tIuIZoT6hI9p18i2zgXwkOjCcIRcw9JiQr4KuTaadjtSMdfSZhdoynd709zXcfuCWKqWskPXrjNxDnHSdN72462HxdI4ThTxP+ceh3875RC2nsr+XnjqoalrfcVKcjHwras+YptDmhACL2/YmSqXz4iLuqqXceMaoywU+nl5z6q5IffBrti0G68xIbVT4yHdWol/rqNCX49I52D4Ss+bFuy31N2nkXPwpO5yd6w384ETxqK/Etok5mJFNF/5zIABYXBxDqa4YaZ6pmrRBkjRFzlh0BkBSV8LlUHkiB3v52Dutq5E2rvThtW9tMFNHIjoDwDmAjiCiNYQ0ScBXAngLCJaDuCs8DeEEEsB3A7gRQD3AbhECCEF8xcDuB6BkvpVAPeG5TcAGE5EKwB8EaHlU1dB/TjHDG7A/K/NxOEhmylfhtyxm0Jny0WnxeC3ENUzvF1dbKPGgblvyQbctXBtOI709TLvrMuSget7dFMDvvneyckUkgzncMxYP/k0BxfnoI/RN9R0uQycz4Ti+MrZwe6ZIyxbdpvzU3DQRy2bVK1mOvO9qgpNINh8cLec1C2Z9RKyXF1gTIsw24ShbyHSHtlSFCghn7fKORQZKlYMxUo2qxx1KO55ba93+hGSK42/KaeBSdj7lt0tqZAhiXeviJ86GP1LQ10R137kBACIxJC2fa6uc9DFSN1MG9w6ByHE+YZTM7lCIcQVAK5gyhcAOIYpbwZwnmsc1YL6+AtEGDVIdeGXu9DgNxfPBogXbd8J70Ks4xCRMxPAL5ZyJ8bt7PR6PuB2YqyIzFshbQ9DnOw7+fHY+lAT7qg4OpThcs9q295sxCE9vqDNesUSqjO7OVUhbcttrYbpNnEXQXvB36RYySR+495pup7Ui6hiUCAZXhtQ/GPKMXHipJGPhXqD7zpCqXCRjvmKwR/uNieOGhi14yIiyb6Dfv/jjy/g4WWbjfVUB089t7bEO48dk/ht4xx0nUO6v+4lD9VWSPd8KC9L3+jId2HzSgXij7fZQRy85fQGeSvPplPCCU7i/OlJBb1PcpugveCv747NBSLzB64G8gvaJO8kPqaPrBAtBukBpiJiMm2IxPlkBdlmncpphS/19s+8Fcs2uk2FE6C0uIHfAMQiFpOjHBAvqOomxfTsTaaZXN8CsU+OCVFAQCEw+cAmPP/G9lSEVwB4Q3NmNEE+BtM3p44PAPsy1ftRFdwEu9WSPL9jn8M0XfFncelQJJw6B4XI6HUnjkqHVqkl+hxxUBfVdPyg4K/KMnOIdQ6ViZV0qHbtKkwLhxDpb8Nk9uruO/6IfMbogil8xtOXz8TopmTcJCJ/JzijWXF439yz0q9Q331rexn9SgVr/zYfhOmHDMP0Q9KLoQ0Fxc5X9suKlQpxHWHhHGSxTW4toQfRC67nuQmfdyKfe3u5jOmHDMPzb2xnlfu+6VnJ8c1F9cK/XDX1diJu3GN6SUMBt0VV3N5zSrBIE4QQuP7xlcbzdWHgxI6ySIRfl9Atv2qNPh1bKR15VIqV4l0RB7kQtVRNrJRmqYsFYj9eGawuRRy0ur6pIrnw4+wYfYmNwVv4gMENqTYCayWvZgOrHct5TqRhu6cov4ClDhcMrjNQ9TFxulhGrKQslDZTVlbnkMGWg1VDGIi7aYzlciBiGdCP98L2n4fBX3/FdbqeSkRjcSnvE6QiiMqarnfzv01P1ZVVloXRBGx4Yc0O/PpJC3EoJd+ffu+uZ9HV6NPEwSR6cXMOwWNzEgff7TbDOdjyD5SFwCPLkvHxTd7ULsRy6+oQOp/wGXFdv+Q3gA9n48E5KAV729pT7aatldLlndY5iGSbbOA9xRLIrpAO/m5XwrhksY3nmi1SegfLIbLoE/HOl6+XLD+KccYD4ufgSxy46aD2JettDjMv2mDiHEYxnK7sd+JIt8hHnVtjhzTi0S+fnjgvPcDlt6d/gz6x0boSfY84KMepwKNSrOTY7Ui236WQ9pbTM2WmtuUCc9mfFqfKVdjMQhN9k72/uJ6/0sHblJWAOUwSGA6+wfxU2AjPnpaOsF1zm7a4RzZ895yj2XJVZGMTK8l3t2FHM+5futH47OX43twdL4BZdpvcvdSXimjVrPCmjBts7LtcFmgvl40h09Xyg4b1x58/ezJbTw7FSZgMOrJRg+rxy389Ma4W1vMmDkg/Oz22k8pH2ji0D7/lIIwY2C8R/mPc0EYcPDzplFkXEYegLV3JzVlE1RJ9T+egvFT945cfi4s4xJyDXefga1sej81dR4qVdHRW52Bz45f9+iDebftYifg/n0rEXvolx4yNFzlJDBNOadoVchOglruG/O7jxuBjb51gHKNsSb5Crj357r5/78sAzFydvHan4ofxoZPMjoSm61U01BXQ0l7G4Ma6SEF7EyNeUS2qfDmHwY11xgCJ/pxDumx0Uz3m/eeZiTKL3poB76ine7mrIjcbhybNd1UfiakHpR0E6zTOQSWM/UoFtOacQ22hPu5TJ41InIuslRwrobfOwfPpxgubDzvPy/T1jVtWU1a3Wa6/zgFwP5ugb68mAUirHfMF7BmL0p4L+aw/VtmfuqNzDdl2nqCEX4iCKzLB5bQiE3HQdQ5PfPXtmDYhm5JcR32piJa2Mt6iKNsHMBFtI+IgBNo6hFemPdv7jq2VXJuUUFykcEt8Tgz/yRVUFakFXw/8p/rlqITky2GYfIliGC1XxcVvS+YnAeIYZvLbU/OvDB/QD/uYFMC1RB/kHIK/1330RJxx5KjEOTmdLPG+AGTwkPYcU0QalPl01uTRfF2DhU9K2ZtREegWK3k1F7WnsvNq7J9k3bhR6aWu4vGvvB2rt+7Fh6+fV1FQQVMOAUBNM2mzVgqJQ1nlHOwPwna+UIjfsQxixwU51HfhLt8FUxTfSlBfV8Dulnb0L8c7fJNJNRDsoDvKIsqmqCPxOKzE3Y9zkE0sWBVbC9mIg8+0CYh2WqTV1Jict4qxGcplgYa6Al7+3jvZvssaJ8LF2tI5B9UvZ1BDCXssCYpqgT7HOcjXO25o/1QMIV+dQ8w5OMRKnTBl5RYNILTwYdZxXcfgHe7Ck3NQExuNGlSPX398Gt+e9AFR6v/h0zOsfQPAKRNHpM6PH9YfB4VEwxRPSoK7XZtIIYqaqdZPjS/4q+7oXE/VNnVUSy65oOvxs4I2ko2YRBiyVuSE1hlteQjJOagLGzePJcFoLwfhM7jMaEDS/8aLc3CIUuSz4WKHce356GA4ncNHZxzM6hzkY+koC+PzLhZCSzPlk+J0MjKCbHtZYMGqrbjoljgf/MD6UmxR103oc8RBvlzuvcrdi2tRly/VJTrxV0hTYmwAv2gAwQLAi5XI+Pvq849nzfKA+Dm47mWbYhHzkbccjDOO5DkbiS174l2Q6UPm7NJ1qAHUXEniddiIwzVhmASbWEniBSU/s2v9tZ1WLbnkgs6FpdaJgx7KQh+LKRZXJagvFdDS3uGW/Sd0DuVE3+dPPwjvOS7wFL7g5AmReMa2YYn8HDytGVRRG/upSI7YI+GVzOPiCgJZKsbJeTqEMHLnhUKQ/S4RzNLq7CgSOcABoH+/EvbmYqXaQr4u1ozP8+OSOwa3n4Pn7p0J9PfEinS4YiD4iLgPSJ+o6v1NnzAsFYkyuk67l2+9N52OEkharNji7G8K8yN8WQm/UG9IxMPZpZvqlIX9Qx/AZDuziZUefGkT9ra2W8VVUiG7Ugmq1xmxkkpk5OLGJbTxFSsREYjU8BXZiAPrZBnG7nIRh6KysLV1JBXS3/9AMqlSsUBAR3U5hwRxsIiV9DzqHCTnoIqVOMX5uKGNkTirbFPCU5DwKmkCbRPNpd9FkAioexXSfY84SM6BWbilKMf5YXh6SPty+ZJLuO6xV6OyrXv4uEDFAm8qahMr2RaNyJQ1/IjGGIjIqKYGDOlfh+1726w7wFfD0NWrt+6Lyrh0j0Byx2dqUg3vYQvpMHxgfcLCBnBbqkz+5v0JfYhOTCpxQvJdnqNIpiznoNV1xN+RC6WvnkmCE3VIazjfb0DWtW0Y5K7ZalCQ0VrJZSQQExEPsRKCuVJ2EIemhjrsbW1Hc1sHbpr7utXxr0MI7zDlZYYL4ZTatUafEyvZ7MulDsIlf9ed4Ljk74A/cZAfmi3ol0TgVcxxDlqbnsRBfkQyhalJnKW2actfvytMnaoq4EwpPAfVx+Imc4gIuQNM3vMBTUFEXRXHafb4Pp+WGo5bh84ZmhaDBDzeeXtHORYrsdZKfpwDEMyHWH/RebGSFFv6ipXWbt+Hh17eZNV3SGW1S+QG+CikfTmH4K+fn0MgVlI5h9GD6lP15IJ941OrAMCoMJa+SPJebvzEScZ6AJ8u1mSVWEv0OeJgEyvJj8vlLVzUdA4fP3kCW89XrJRFHCAza6Xa0H02lJ924hD8/eZdQTavOsaqQm/HxjnID0YVl3CiEyDtgWob3+6W5CI+YlC/RERdblyd/bZ0znDOl053XuPzzpvby9GCzlorafdhcwwjij2aqxHFs1Cg0LHNT6z0H38MxIeL1+4w1i15zJtoE6D0++lTD0mPL2xiqxpxl2lWvgeTeFaHQJJz4BJeRVZInlyV5PhM961mduTC3/h4qncl+h5xiFYMTqwUPA5paXPs2LRnKKCYsrZ3oEBmdt5frJSuOJLZuQAwJtOxKaTtYiVN1mmpK9uxtdc/3F1LzuFth480ihNUouHSOejcHJ/FTC/J9nEdpoVEGDEw+Q6GDrBbyPBjSGNva3tkz88tQr46ByDJAdvey9fffVTit6lqMdRpOf0NZHDAsJptHYsTMpnryOemegXrecSDekFFVWzJcS2cj9GMQ3kfECIAjBOcjmIh8MNwcWgRcXAQbVlcFgLLNyVjNRVCvUV3os8RBwmWcwh3cW9s3Ysp44fg9s+8NV0J8csPnH/MjzCrOamK/3pfKvVFWNckVjLrHGyTWe/bxjnIurZF6Jow2cnIcGH9mrYoqVDvwqxzCImDxs3ZQl1H7Wf8tvTFSI8y6vM+fZjA5tZyJAvniLHejZVzUIikTbSjh26wxe7qKLsVw3Lem0SGXF3b8+N8Snxfn89cOHXSCPxs9vHs9UHgPR+OoMAqj019R7ogx9ye++qWVLrTUoE3PKkl+hxxiBXSachJXBbAUQcMQqNJ4eSr7GXKfv7h9AS1KfN0+IqVCp6cg37KRkjkOG0f+dghQZ7h5va0eMkGU7dyB6hzDtyipDfR2U+LiDCkf8wtmBbfH583RRmD+dl86p8CMUlZiEgkyVlyZfFXkM+NLBxsMC6tD5OlTSGZ+tME+Vp9FL5y3lgdBBmdg+/O2ZQUS8VXZx2ZSs8rQYSUKSvfT8BdmOJIqfWA+NkYTV7D57GKSTFbKJCTQHc1+h5xCJcMbqKq8l/rglqgaHdnW0y5b+E9xx2YKsuqc1i9LWapZR/6hE2IG6yy3uRv28QvenAO8rlK0RznGWq7TodJrMTH89c5h85/XOr7NX3kHzxxHA4aFuYht7xKGdupLESkz+AsubLOB8BNUGwGC3r51j2tWLHJHpJa9rvPYbEHxPdjG6EkqipR8pWq+HAONg4nlCp5me+6LLPUvl06B5PQ4a2HDjf6M9USfY842DgHZWF0yRXl+V0tZmsXIvKyIGFlpoa6hQKwfkdMHH507hQ8ddkZqUU44UOQRedgIQ62xDrJvuNMYv7EwdwWkCYOnI+JPqxqfFqxEt5eT93BmyDPCcRZBFkOSGvkh+ce52zT9U7SuTTMYiUf6P2Zwr0AqkLaNr7gbzkhVvJ7gzYPbgmTObXsWwg3MZLPxlfn0BbpHAztGd7B7z/9FpSKvVwhTUSriGgxES0kogVh2TAieoCIlod/hyr1LyeiFUS0jIjOVspPDNtZQURXU5ZwnRlh85Cu89glSvju7uQuo7GuiKvP52WebIpLU79EidAUxQLhwFCUo9fzgd63bVdU9N2lKmO0EYeEk5CBHJp0DuedOC5VVx/WBE3OXgmkqNEdxl2edxNi4eAc1P5GN9XjvGnmSKuSCDuJg3G8Wr8VzptrQ10T22aFOgdv6s7ssHVFv8kREwhNWX2CXkamp456mrWSiytWmzt8dJAHe39RSL9dCDFVCCGD7VwGYI4QYhKAOeFvENFkALMBHA1gFoBriEh+GdcCuBDApPD/rCqMi0Vsq8SwosoH49wdeH5EMubMGUeNwj9PSYuUgjH5TwIiSjiDma71Ja/6bQ7oZ/aL9LFWAoLnKHf2vjoHE+R9qJzDk5edgY+85aBUXRniY0C/Im78xEm4/F1HdqpvFS4OX37otuceOz0BX75jEQDzoiXnl0sJ3hr5S/iJOqL2jTqHyjgH23V+TnDB34Wrt0dlpkd+6dsnJn77fD1WzgGe4fILvIgzVU9TSNv0Ozr+59xAf7W/KqTPAXBTeHwTgPcp5bcKIVqEECsBrAAwnYjGAGgSQswVgZD4ZuWaqkPKobl56iuKAfzloXInbiMm3Bww1V66LmlPbrI4zJomVGKgIYKq2qbr2cjT/UqFREx7Gxau5nPy6uE9gEDpzS00bzt8JIBgETr9iFHWBcEXa7cHIjxfu3/bo+GymLm8x30t3lxK0nS+dL5dX6ZdN5W2XRfpHDyejQoTwdPnn88aap2H5Edgip7EQbVmBCzGFsw9R3GoChSFdVfR1lHGc29sSyR46ip0ljgIAP8gomeJ6MKwbLQQYj0AhH9lXOyxAFYr164Jy8aGx3p5CkR0IREtIKIFmze7vYlNAzZBfYmunZgvVZfsbda4N6Z0irp1yEmG+P3+yX6Sv/sbkrEAis7BQ6wEABOG9/debJx+Dh4xcnwUn10FOXybtZI8oyoaTYpSed823wAAOGdqwI26lKT64zXREl8lfoI4ON6xj85BP3f5O4/E9EP4ua3358N5m5IMAeE787ht2e81j6zwqudUSHO6Rorb4NaYbXtb8YFrnsK9Sza4B9xJdJY4nCKEOAHAOwFcQkSnWepyT0hYytOFQlwnhJgmhJg2cuTI7KNVWubel7qQuSa8rzywLvowzO3pnsIrrngnxg9L5zfQ8eAX3xaFtNbhL1bSdpQe4gGXpEi2afsgdZh2drFCOpDRf37mJEu/4YHl3nVnMNn3Hy/ifVp84bM7jnUOcZlJ7OYrVmoMn7HLMztFHAzt6puPow/kNykljXOwwcezXt9EXHjaoRZZvbU7Y//mvv0IjGxHii9/NnsqPz7NCc5025EOQ1lLIuuzIm+yLnPI+HLknUGnehBCrAv/bgLwZwDTAWwMRUUI/8okwWsAqJq1cQDWheXjmPIugc2UNWH+6diy+ZqZlaJwxeY6/fuV8Llw0Zs0aiDrNSsxVLG7NyXRyTI+dVjnMkpeFb4Lljzt4r7UD9IkFpHvSYqV/mlSOu+DRMFjfO9l9D6fPvVQIwfmi0jnYKkTWeQ4ciUA8QLjetaSAG8Io+G6xifxGSYzGZAMHfP+48caHUEzWfZ5+Mfop+xJk7JRB1d1aaLqgk5kzjyKt9CSj6bVoXOQt9iWIA72Mck2ezRxIKIBRDRIHgN4B4AlAO4GcEFY7QIAd4XHdwOYTUT1RHQIAsXz/FD0tIuIZoRWSh9Trqk6bKas6uR1hQ/IrHNwLZThwHQ3eh0fnXFwdDy0fz9jvT0tfolC1A/tiNGDrHWjRd/TzttG5HQcbum7WKBIzuvj7Z1Vtm0jsr4oRJyDe3y3PvOGuz0ppnIsbL7cmXrfz379TFxgiAem+i0cNKw/myIUSPoEuRZrafGV9b2YoH9LndXb9isVvL5nLnIqW09uaNrMwRXVeq1K0jA5f2RkVx2X/C5ICFQL4tCZr2I0gD+HN1MC8HshxH1E9AyA24nokwDeAHAeAAghlhLR7QBeBNAO4BIhhHwqFwO4EUAjgHvD/12CyFqJFSvFx9WKpS6tlVwf0Gub016SHGQ7gxpKVjNR/9Ad6bZNkHO1waHolc2Ystlx+PyZdnGRJA5WPwySf839ch90NT40eat24hX8veVpN3GIw8LbNym+z1h9JDYCpoYHt7WtZn5zbXx8rJWyMAN63aXrdlrru0Ruvr44uijONBf1HPPGrI7h5fcv3RiXSc6hQFEYcfW7fHnDrkxj7gwqJg5CiNcATGHK0ARLMAAAEyRJREFUtwCYmb4CEEJcAeAKpnwBAD6YUJVhy+egLipZzEtt8LFWUus52wsnihrWoTNQP1jXByrFIbb4S0E7Iefg0qYqsC36RBSbbFqeUyyKMffDvYd+VbBqUmXFJmTYHEf5PKS1lKtfFxKbAMslKqdg4/zU9+AyV94bRuq1OsFlMCOoRvRZFb7m1mlHU76enGPS18dIRFjRdnK96BAChfDZbFNyvFTDEs+FvuchDbOSSC2rlgNKnYfOQa3nglwAXTl0fXUOPtnY9DZd90IRcfDjREwhz+NxQREr2YlIeGRui7m8k64YYbvue9Z3zpMNFmmAv9jSF2rftl3nv58Rc3C2e1Hnqyv43tzXtgAA9lijCcTHLuOArJZ/Lks33124zgGYOCFZLPM+mDY08vok4U5uMlS9wxLFjL0WnEPfIw6eOgcPy0kvyAnlEtn4Egf5wTY6ZM3+xEE59vTtcBERuQiYch/rOMhhmVUkX51D8Nc2vGrkWbb1bTNk0Hs+2GBpBvhzGb40RG3OtlM+YHADxg8LPO6tQRiVc767WF8lc1aHPp/IsDaoC+2ghhLuvvQUtp6vbuPN3cksjqbnrTpFSqimrECSOKjGAqMH8UEEq4m+RxzkATP/EmGuPcU8prwLUTsy/IKnLbgLPqEIAH/OR2XnXUMQEedgryiVmo+9Upkvio6CShw8wnvYRseNvb/FK9wXO/YFBNGWD0Pv2/YYvZWsYcUvWHQ2QJIougwFpC2GrZ7aXjWkPGoTLmKjb2I+PP1gQ00/qDqnL511OI4bN4St5xvrSOeQTM9xUEOa+4+IQyEWK0lIPegPzz3OaMJeTfQ54iA/Jl7nEB+7duYSP/jgsdbzMoNZtXascm10i4D82ssiSrvi/cfilInDceQYu1VTtUEEtHT4KKQp8ZcD9x70vA0qrv/YNOM5FS+tD5SiWXQO1ZAb28LBqMii/9nbGszZAfW2YHVpJakJMo2ubYTqzDvyAPv8SudLt1Z3Qt0cFC3zy3fDpRMHk0J6+IB+qfmoWivpfUou4thxfBKyaqPPEQe7tVJc6OIIovYc80WmT1yzza5Y9EXR0+O6Ep3DI44c1seMHYzffWpG1ZVhrpH6mrJSRDgtbWkv/sgDBlkJzpmWaKMcfILLSVTDSsoWSFJFFssx6eRVLdGFTzpYGYhwwvD+HkEvk7+z+j3oOHBIfJ82Dt93Tfig5i9UZzJlLRAGaqbCkhjo2eQAWPOOdwX6HnGw6hzi49knmSNhcu3VCpHoxCkC8muvRvOMhRRTudatpLjPzTlkcaD64bkpg7sUPn3qIXj7EXaP/JlHBlFibM89zTmY78WlpJfwtaqrhHNtchg9HDrCL+qtj8hU9vWuY8c46+piQNPrnu7p2DhE8ReySdxOd8wBiYOHD8Dho+OUszbipc+BDo04qJyDLe94V6APEgebh7S/XFbiOE8Wr1piJZ9k7YB/7Ce1nU+fdmjlA2PgImCRgtsZ3iA+b/swJKHbtrfVWEeHz3v52rsn4zefmG6tMyl04rNxbFL0I3egnz/zcGNdX3NNqaN03YevwUOWa1zEIwtGDKzHM187E//xjiOcdccOTYaoNz2r2z4zw6tvdU5xeUIkiAhTPL93OSaXY6ludSQ5hSg+U/j7yRVvYnMYbC+Lc2ln0HlNXC9DLKNNI4v99JLvnA0hBKtU4uC7WLs4Ft/kM1NCpdpXZtk/NtnO2CGNmHHocK8x+sKlt+mwEGoVagRKG0u96s29AGK7eh/UknOSz7pfsYCxQxoxbIDZw10+Epf8Xcq3dfGEjiypaCVcBMfXSmh7KKayRfwF/MU2elj50ybxO3rfoI+qNdG2PW3Wuq2ezrHy2Z1w8FBrPbnZG9RQwq7mdowYGMwJuWHqKAu0tHfgI9fPi66xGT1UE32POFgD7/m34/oYdXDhdzmYorFKDA8nj8ty4q2HDcdz3zjLugAFSCrAqglXaIdYrOTft41zcIU84fD/2zv3GCvKK4D/zu7ykPcb5LEsEkVXauWloDzUFBRpi621wVoWxQZobdTatJVqbVrbBJvWVLGJEgvVPnylNuKjGjX1WVuUAgoigo/4IlJTUYQILnv6x8zsDvcxM/cyd+7M3fNLbu7st3PnfmfO3Dlzvu9855Sy+CoI7/QFZTX154kKmux1+uXw1YkFExS34xmHYmkuPMoZpw4bvvCumZZpwdFCs5uHct/G9/h+BK8gCn5Dt+OXcw/7SdrvITUNCo4C8gIPRvbPL7DlZ8+njl7GhBzPMwLfPm0sS2eObT+nntFoU83L1pCU59D5hpXc90I3hZgXXgJw/kmOJxDVcwjrw4TG/owe2CM0SR4QwTDkh87FSajnEFJGsRCBdQPKUGBcYncEOgRNSDvvB1oPho7DS7uxCf7ellOa6NpQx8yAhIQQPVTaT9RrYk5z8WgvgIG9unHHkqmMKFCxsBz8T/px3Cj9xqFYQa5cwvbzVjMP6xsss6eXLnV1eTUywHkIzK0fUam1Onl9S+RbUoQGzEjHvSwfYOpRA7lj3duhkSne/8Oe8Poe0YUnf3B6bP1rK+MGHZXF08cEf3fEOYfIlKG/uFT+aXvZz/AV3AcOtoWO57en9w753hNH9ePVX8wN7V85w0rlzFMkQTmyBOH3kKIORbVMawr8v1dbvl/IvEz7iuic34DX3tamhyx+y+1vJUmn9hOg0DVQCePwsetehs1NXD77GJbNGhvJI4gTz6OpxNPIxaHGIdqiushE9M78NR3iKlfuGYcjugYXsgcnFUi4cXDe44qGK+dGn9QTaqnEHcpZyjXgVRsc1jdamG/YpL13jnNv+N7fnx3M9xzCEl/GRSf0HJz3sFDWuJjTPJSVj2/n4ulNgfv16tbAlXPjq3kclY6hneRvBJ4u4jrvUVew+gspxSV14wAnrHP0gOLhne1PgxqhchsdY85x4BmHqOGdUDw+P5ekL51Sn5zDIoZKYfWFU0qa2wobMegooHXoft5aogMH2/IiqGLztEPofMYhoNhPJepLDu3TnXVXfSH+A8dEW8RQyIp8d8yeQ5SCLXDo03hcCe6WzDyKiY39ODkg4ssv5b6QehsLpoziruffai8DerjU1wkPXjqdpoHR1iZAuAFLeo2PRynzDJt/dmYkYzJ2cE8mNgZHFoFzHuvroj+5hyXIK5aw0TMq+z87yL6AhIWVpPMZhwDPIaYs3ZnCu0EnFQHhx7uZx/XkGdVzOGpwxw0yrFh8VOrrJNAwwKFGcNv7wSknmgb1ZMM1c2Lpm8fxw0tLuxD1gSHpx4pSJtejRhU+dsWscrsTSFg68A1v7QY60uy0f84zDq1tkYpDVYJON+cQlD7Du7cMjbDcv1YYP6IvLdNGs3LBhMS/2zPUcXkt805wVteGPW0fM7Q3w90x4/2t0ddEHC5VGLk7LMJuwl+f7ETijR3SK3C/uKnERLmIxDb/5CdqipRN7+zO+ZzjnexvbaOMCO1Y6HSew9zxwzhmaK+C+YFinyDNAPV1ws/nx1tn6eHLZ4QuJoL4z/fExv68uWJepH3PGn8kq599I9ZVvmFk5bK69zun8MCmnaFG+5wJIzhnQvA6jEpQXydMaerPN30lc9NK1LoLuQtQu7kLDK+4a2N75BMc6vVWmk5nHEYP7MnoIuOu3o3ivIQjhmqNY4cFL+Tz8HLa9AiI8KkUy88+lq9MGMHYwck99fqN4HPLz0jse0tlYmP/SOPv1eSeZadUuwuRCPNynv7h6Ty9/YO8zAiex7EnZ7jp7qXBhZDiJDXGQUTOAm4A6oFbVXVF0n3o1a2BV649K5Hi3YYTUto8vE97eGAY3oLCOOhSX5dY6mMPv+cwrE/li7UY1WPMoJ688cHeUM9h1IAefOPkxrz2Pjmh74tPHcM1X2qOtY9hpOIuKCL1wO+AuUAzcL6IJHsmXLp3qa/I2KORT89uDSycOjr0fN+ycBJAok/5lcAfn27XWG1z99Jp3Noyuez5kdzhznknBK9CrwRp8RxOAnao6usAInInMB94uaq9MlLBnOah/OGiKUUTrGUFb+FU7xLzchnZY3DvbiXXAvGTO98zJIGyoLmk5SodAbzt+/sd4OQq9cVIGSLCaeOGVLsbh033LvVcPe84po2NN/utUZu8uWIeB9uU7bv2HLJwMynSYhwK+dh5QesisgRYAtDYmD9OZxhp51sz4q2ZYdQ29XUSOcAjblIx54DjKfhnG0cC7+XupKqrVHWyqk4ePDjbQwyGYRhpJi3G4XngaBEZIyJdgQXA2ir3yTAMo9OSimElVW0Vke8Cj+CEsq5W1S1V7pZhGEanJRXGAUBVHwIeqnY/DMMwjPQMKxmGYRgpwoyDYRiGkYcZB8MwDCMPMw6GYRhGHqLVKud0mIjIHmBbgX81AlGqY/QFPkrxflA7skSVoxLfXS1ZStFzrchSTZnT/lspZd9KyzJOVcNrp6pqJl/AC0Xa/xvx86vSvF8tyRJVjlqSpUQ914QsVZY51b+VNMlS7N6Z+6rFYaXd4bsAcH/K94PakSWqHJX47mrJUoqea0WWasqc9t9KKftWUy/tZHlY6QVVnRy1PYvUiiy1IgeYLGnFZIn/+Fn2HFaV2J5FakWWWpEDTJa0YrLEfPzMeg6GYRhG5ciy52AYhmFUiNQbBxFZLSK7RGSzr+3zIvKciLwkIveLSB+3vauIrHHbN4nIab7PTHLbd4jIjVKFOo0xyvKEiGwTkY3uK9FKOCIySkT+ISJbRWSLiFzmtg8QkUdFZLv73t/3meXuud8mImf62quql5hlyZReRGSgu/8nInJTzrEypZcQWbKml9kist49/+tF5AzfsZLTS9QwrGq9gJnARGCzr+15YJa7vRi41t2+BFjjbg8B1gN17t/rgGk4hYX+DszNsCxPAJOrqJMjgYnudm/gVZza378CrnTbrwSuc7ebgU1AN2AM8BpQnwa9xCxL1vTSE5gOLANuyjlW1vQSJEvW9DIBGO5ujwferYZeUu85qOpTwP9ymscBT7nbjwLnutvNwOPu53bhhIRNFpEjgT6q+pw6Z/h24JxK9z2XOGRJoJuhqOpOVf2Pu70H2IpT6nU+cJu72210nOP5wJ2qul9V3wB2ACelQS9xyZJkn4tRqiyquldVnwE+9R8ni3opJksaKEOWDarqFTvbAnQXkW5J6yX1xqEIm4Evu9vn0VFFbhMwX0QaRGQMMMn93wicanMe77htaaBUWTzWuC7yT5J2+f2ISBPOk86/gaGquhOcHwSOxwOFa4SPIGV6OUxZPLKkl2JkUS9hZFUv5wIbVHU/Ceslq8ZhMXCJiKzHcdMOuO2rcU7YC8BvgX8CrUSsUV0lSpUF4AJV/Rwww30tTLTHLiLSC/grcLmqfhy0a4E2DWhPnBhkgezppeghCrSlXS9BZFIvInI8cB2w1GsqsFvF9JJJ46Cqr6jqHFWdBNyBM+6Lqraq6vdU9URVnQ/0A7bj3GRH+g5RsEZ1NShDFlT1Xfd9D/AXqjCsISJdcC70P6vqvW7z+67r6w1N7HLbi9UIT4VeYpIli3opRhb1UpQs6kVERgJ/A1pU9TW3OVG9ZNI4eNEGIlIHXA3c7P7dQ0R6utuzgVZVfdl12faIyFTXpWwB7qtO7w+lVFncYaZBbnsX4Is4Q1NJ9lmA3wNbVfV637/WAovc7UV0nOO1wAJ33HQMcDSwLg16iUuWjOqlIBnVS7HjZE4vItIPeBBYrqrPejsnrpdKzXTH9cJ5mt4JfIZjOS8GLsOZ8X8VWEHHYr4mnEytW4HHgNG+40zGuSheA27yPpM1WXCiMtYDL+JMVt2AGy2ToBzTcdzZF4GN7utsYCDOJPp2932A7zNXued+G74Ii2rrJS5ZMqyXN3GCJD5xr8nmDOslT5Ys6gXnIXGvb9+NwJCk9WIrpA3DMIw8MjmsZBiGYVQWMw6GYRhGHmYcDMMwjDzMOBiGYRh5mHEwDMMw8jDjYBgVQESWiUhLCfs3iS9br2FUm4Zqd8Awag0RaVDVm6vdD8M4HMw4GEYB3ARpD+MkSJuAs0ixBTgOuB7oBXwAXKiqO0XkCZz8V6cCa0WkN/CJqv5aRE7EWfneA2fx0mJV/VBEJuHk0NoHPJOcdIYRjg0rGUZxxgGrVPUE4GOcGhsrga+pkwtrNfBL3/79VHWWqv4m5zi3Az9yj/MS8FO3fQ1wqapOq6QQhlEO5jkYRnHe1o7cNn8CfoxTfOVRN+tzPU46FI+7cg8gIn1xjMaTbtNtwD0F2v8IzI1fBMMoDzMOhlGc3Nwye4AtAU/6e0s4thQ4vmGkBhtWMoziNIqIZwjOB/4FDPbaRKSLm3O/KKr6EfChiMxwmxYCT6rqbuAjEZnutl8Qf/cNo3zMczCM4mwFFonILTiZM1cCjwA3usNCDTiFmLaEHGcRcLOI9ABeBy5y2y8CVovIPve4hpEaLCurYRTAjVZ6QFXHV7krhlEVbFjJMAzDyMM8B8MwDCMP8xwMwzCMPMw4GIZhGHmYcTAMwzDyMONgGIZh5GHGwTAMw8jDjINhGIaRx/8BWGCaLqGUPJMAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEKCAYAAAAMzhLIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXmUHOV97/19upbeZ3o2jUYzo10CJAQCZMArJpgAtm+wE+Pg3GuThIRcx/Fxltf39ZLFTkIc3yx+45sbX9uxY/CNjZeYYMcbGPBGBEKAQEILCC0zI83e09N7LV3P+0fVU1299/T0UjP9fM7R0aime1Rqdde3vr+VUErB4XA4HE4lPJ0+AQ6Hw+G4Gy4UHA6Hw6kKFwoOh8PhVIULBYfD4XCqwoWCw+FwOFXhQsHhcDicqnCh4HA4HE5VuFBwOBwOpypcKDgcDodTFbHTJ9Aog4ODdOvWrZ0+DQ6Hw1lTPPPMMwuU0qGVPGfNCsXWrVtx+PDhTp8Gh8PhrCkIIedX+hweeuJwOBxOVbhQcDgcDqcqXCg4HA6HUxUuFBwOh8OpChcKDofD4VSFCwWHw+FwqsKFgsPhcDhV4UJRJ8cvxnH4XLTTp8HhcDhthwtFnfzdw6fwx/9+rNOnweFwOG2HC0WdxLMaljNap0+D0+VoOQMpRe/0aXC6DC4UdZLI6khk+QeU01k++5NX8F/+8eedPg1Ol8GFok5Sqo6koiNn0E6fCqeLuRDL4sJSptOnwekyuFDUSdJyE8k6bf98QsF3nr/YylPidCGKnoOiGzD4DQunjXChqJOUkgMAJLL15Sm+9vQE3v/V57Cc5nkNTvNQdAMAkNVzHT4TTjfBhaIOFD0HNWd+QOvNU8zGFQDAYkpp2Xlxug9FM9+HGZULBad9cKGoA+YmgPqFYj5hCkQ0pbbknDjdiWI5iazlLDicdsCFog6c5YjVQk/Tyxn86PgsAGAukQXAhYLTXFjoiTsKTjvhQlEHThcRryIUX3lqAvd8+TCyWg7zSe4oOM3HzlFoXCg47YMLRR2kVKejqBx6Sqs5GBS4EMvYoadFLhScJqJYApHhQsFpI1wo6iCZrU8oNCvhfXI6gayVdFziQsFpIioPPXE6QE2hIISME0IeJ4ScIIS8SAj5gHW8nxDyCCHkZev3PsdzPkwIOU0IOUUIucVx/BpCyFHre58mhBDruJcQ8jXr+FOEkK3N/6c2jrN3olroiX2In51Yso/x0BOnmdg5Cu4oOG2kHkehA/gjSullAK4H8D5CyB4AHwLwKKV0F4BHrT/D+t6dAPYCuBXAPxFCBOtnfQbAPQB2Wb9utY7fDWCJUroTwKcAfLIJ/7amUZjMruwoygkFDz1xmold9cSFgtNGagoFpXSaUvqs9XUCwAkAowBuB3Cf9bD7ALzN+vp2AA9QShVK6VkApwFcSwgZAdBDKT1IKaUA7i96DvtZ3wRwE3MbboA5ioGgXFUoFCv09OKFOABgNOLHUpoLBad5sD4KLhScdrKiHIUVEroKwFMAhiml04ApJgA2WA8bBTDpeNqUdWzU+rr4eMFzKKU6gGUAAys5t1bChGJjr69qeSxzFKw575KNYSwmuVBwmgcvj+V0grqFghASAvBvAH6fUhqv9tAyx2iV49WeU3wO9xBCDhNCDs/Pz9c65aaRzOoIyAJ6/VJdoScAkASCbYNB7ig4TcMwqH0TktF4wx2nfdQlFIQQCaZI/Cul9FvW4VkrnATr9znr+BSAccfTxwBctI6PlTle8BxCiAigF0DJOjlK6ecopQcopQeGhobqOfWmkFJ1BL0ienxSVUfBqp4AYCjkRX9QRlrN8TABpymojvcXT2Zz2kk9VU8EwBcAnKCU/r3jW98GcJf19V0AHnIcv9OqZNoGM2l9yApPJQgh11s/8z1Fz2E/6x0AHrPyGK4gkdUR8ooI+0TEM/U5iqEeHwaCMgCe0OY0B8XhIvjNB6ediHU85rUA3g3gKCHkiHXsIwD+GsDXCSF3A5gAcAcAUEpfJIR8HcBxmBVT76OUsnf1ewF8CYAfwPetX4ApRF8mhJyG6STuXOW/q6mkFCYU1R2FmjPQ4xMRz+q2owDMXorRiL9dp8tZpyiOibE8R8FpJzWFglL6c5TPIQDATRWecy+Ae8scPwzg8jLHs7CExo2klByCXgFhn4iUmkPOoBA8pS+JqhvYPhTCkckYNvTkhYI7Ck4zUHQeeuJ0Bt6ZXQcJRUfIKyHsM3U1WSGhreoGRiN+bO4P4PJNvbZQRPmocU4TKHAUXCg4bYQLRR2YoScBPX4JQOXubDVnwCt68NP/cSN+7brNGAh6AQDRFF9exFk9WUeOQqlTKJ6bWMIVH/shFpP8ZoXTOFwo6iCpsKon01FUKpFVdQOymH9Jwz4RgodwR8FpCo2Enk7OJBDP6phezrbqtDhdABeKOkgqOkI+M5kNVN5JoeYKhcLjIej1S4jxdaicJsBCT7LgqTuZzd57ztJaDmelcKGogaobUHUDIVm0cxRVHYVQ+JIGvQKvUOE0BeYoegNS3Q13yxlLKBxuZCml8mGVnBXBhaIGbCBgyGc23AFVchRFoScACMpiwfTZbuPTj76Mjzx4tNOnsS5gfRQRv1R3H8VyxhQEZ9jqni8fxh9+/Uilp3A4JdTTR9HVsIt80Cui10pms7s0J4ZBoRu0RCgCsoB0FzuKp89FMRFNd/o01gUs9BQJSJiMZup6jh16soRiLp7F0+eWsGtDqDUnyVmXcEdRAyYUYa9oVz2VyzmwGLBUEnoSCzbkdRtJRec5miZhh578ct3J7OLQ049OmJN2lvj/CWcFcKGoAqUU9/3nOQDm5FjBQxD2imUdBRMKb5nQU1rpXkeRUnTEsxoMwzUTWdYsTCgigfpDT/lktvn4R47PWMdVuGhKDsflcKGowqcfPY0Hnp7E7924E1dtNhf49fglxMsJhfUhLgk9eYWudhQpJQdKzaZFJydn4njoyIUOndXahPVO9PolKLpRl/g6HUVK0fHEK4vwSR7oBu3q3BlnZXChqMLPXp7H/vEI/ugXd9vHev1SeUfBhKI49CSLBRvyug0mkstFoY6PPngMH/4WT3KvhHzoyQyBZvXaroK9VxXdwPOTMai6gTddNgygfAiVwykHF4oqaAZFj1+Cc9leTaEo6yi6O/QEFBYAHLuwjGfOL/ER7CuECQVr/KxVdq3lDNs1qLphu7qdViKb70rh1AsXiiroOQOyUDj8r6JQVEpmyyJU3SjYVdEtKHoOWs4Mjzhfs/sPnrO/5ne19aPoOXhFDwKyJRQ1RNb5miu6YYvyxh4fAJ7Q5tQPF4oq6DkK0VP4Eq3UUQS95oe6G0tknUn8mFXPn1J0PHTkIoZ7zDlY/K62fhTNnCXmkwUAtXdSON+nqm7YjmRjrykUMf7ac+qEC0UVNMOAWOQoIoHqjqK04c78UKe7MKHtTJay12w+oUDRDdx4iblifYl3CNeNohvwSgL8kvmeyqjVXarTrak5w06Gj/T6S77P4VSDC0UVtJxREkrqsSpOiu/mmKPwCsU5CtNRdGNC21ntxS5KbPzJeH8AAA9/rAQWevJJ5nusVjKbdWUDphthjoK7Oc5K4Z3ZVTBDT6U5CgCIZzT4JAFpVUc8o1cOPVmOItWFvRROcWQlxWyg4mZLKKL8YlU3im6GnvKOovx7Kp7V8PjJOeQc5bNqLl84ELDmlnFHwakX7iiqoOUoJLE0RwHkQyn/+Nhp/Mpn/rNy1ZOVeOzGXgqnOLLXi1XeMEcR46GnujFzFAJ8TCgq5CgefPYCPvDAETx5ZhGA2QSq6gaymgEPASSBoC8gc0fBqRsuFFXQDQNSBUfBLnwXYxnMxrMVq55CLJndxY6CkNLQU39ARtgrckexAhQ9B6/kgb9GMntqyZyt9dhJc1zHYMhrJbNz8IoCCCHoC0g87MepGx56qoKmGxCF6o5iOaNBN6gdUinXRwF0p6NgyezhsC/vKKzXKewTEQnyXR31cOhsFOcWUnWHni7EzIGBC0kVYZ8IvyxYeTXDzm9EuKPgrADuKKqgGbSk6qlYKOLWHfKiFUIp15kNdHeOYrTPj5j1erF94yGfiL6AzPciwHQAt/3DzzCXKL+F7v8+eR5/8d3jUDTTEfhrhJ4uLOUny0YCkiP0lLPDVqaj4K89pz64UFTBbLir7SgAIJo0P3TFQwGZo+jG8ljWkb4p4s8nsxUdPskDSfCgLyDzWn4AL80mcGI6jtOzybLfTyo6ElkdF5ezZsOd9Z6KZ8q/py7Esvb7NuKXIYseszzWciSA6ShifJc7p064UFQgZ1AYFCUNdz3FjoIJBXMUxUIhdXfVk+ghGAp5HaEn3V4p2xeQeI4CgKqb1UmVRr0wFzafUOCVBHhFAcM9XjsX4SSr5bCQVHDjpUMAzBsbWfDYJd3MUUQCEhKK3pUTAzgrhwtFBdgHqDj0xEaNs9g6uwCy0FNxMlsUzLr3rnQUim4vfEpaF6VEVkPYSvD3BdfXXW3OoHjnZw/i8VNzK3oee69Veo84J+8yRzDeFyi7EGp62Qxf3XTZMEJeEZGAZDoK3UDWatgDgL6ADIA33XHqgwtFBXSrBl0qEgogP2o8q+XsJiYW7y12FIA1QbYbhULN2RcrwHRfpqOwhCIgI6Hke1D+5YmzeGk20bHzXS1JRcehs1F88/DUip7HhKKS60yVEYrN/QFMlhEKlp8Y7wvg//vV/fjvN+yAVzST2WaOg4We2BIu7ug4teFCUQGdOQpP6UvE5j05d2cvJlUQgpIGPcCaINuloaeALBTkdZKKI/QUtO5qMyqmltL4+HeO49+erXyRdXuYhK0qfeKVhRUtasoLRfmbiWSBUJiOYLw/gOl41v47Hz81h7d8+md4ec4U2rE+P960ZxiXj/ZayewcsrrhSGabrz0vkeXUAxeKCrCpp8UNd4BDKBzJxGhKhSx4CkaSM9hOiuMX45herm/X8XogyUJP7O41oyGR1ezekj7r+FJKw09fWgBQ+WJ5+FwUl//ZDzGfUNpw5o2haOYFP5bWcHw6Xvfz1BzLUdQhFFLeUVCadxDPT8bw4sU4vvzkeRACDFsTYgHkk9laDj7r/cxybYksFwpObbhQVEA3rAa6Mg6BCYVzOGBGy5UNOwFAQDa33P3Glw7hnZ892DV2P6XoCFk5CsB0FMWhJ8AM2/30pXnrOeWd1+RSGopuYDZevoTUDbAwJAD8/PRC3c9jobdyE4YVPQdVNzuqAUfoacDsbJ+0hII1Mp6ZT2E47Ct4L8qCx54ey3IUrJ9C1d3t0jjugAtFBTSrEqW44Q4ABsMyFpJKQegJKC2NZQS9Is4tpDEbVzAZzeD9X32uK3ZIp9Ucgl4B/ZYgLCZVJLM6QkVCMZ9Q8MQr5oW10npO9v+huPjCpjiG9D2xAqFgoady/3YmnDuGzGVDduipzxQKltB2OoNNEZ/zR+ST2Q5HwX6Om19PjnvgQlEBjTmKMsnskV4/ltIaZpcL726LK54YQVm0u2Vv3jOMn728gDMLqSafsftIKjqCsogR68I1GU0jqTpzFObvDx25YN8RVwo9sREpSh3rPzsFu+ju3BDCU2ejdYcZNeYoygqFeWzfaC+A/M3IhrAXsuixE9qJrJkPAoBRS0QYXjFfHstCV+znuPn15LgHLhQV0FmOoszFfzRizvM/OWMmDlmnbMXQk9UgBQBvvWIEQHeMHWflsazu/9RMApTmV3kOBL3oD8r40Yk5iB6C/eORiq+LZguFe++AWY7id96wHQTAJ753sq7n2cnsMqEnJqBXbY5A9BAMhEwX5vEQjPf5MbGYF4pLN4bxa9dtxlv2bSz4GcxRKLoBn1j4XnXz68lxD1woKmD3UZTJUYxYG8JOWAlLZvWLu7gZbIzH1oEAhkLmLoBu2BWdUnP2hr+xvgBOzpivF0tmy6IHP/sfN+Jr91yPB+65HpsivoLwy9mFFD70by9Az+VXybKLsRthd+c7NoTwOzfswLefv4inz0VrPo8ls8v1UbAE97bBEH70hzfgLftG7O9t7g9gcikfegr7JPzV2/fh1stHCn6GLHqgGxQZR8Od7Shc/Hpy3AMXigpoFabBAuZICsB0FD7Jg36rzLOWo9i7qddeY1lr3/FaR8sZUHXD3scx1ufHeStMwkJPgJm/uW77AA5s7beqw/Kvy89fnscDT09iNqHYVWhuDpWwu3Ov6MF7b9gBSSB49ETt5rt8jqL038a6soNeAVsHgwU5s839gQJHwYoEimHvS0rzApHPUbj39eS4Bz49tgKs4a64Mxswdw4TYlbxDPd4C+6Qy8EcxZ5NPXaYaj07iq89PWF3CDNHMRrxg1r5+0oXtKBXLAg9Za27XVU37OocN4dK8kIhwC8LCMgiMnU0Wtqd2WXCbqwru9xrNtzrQ0LRkVFziDtGoxTDRAGA7SgkgYAQd7+eHPfAhaICWpWGO0nwYEPYi9m4gh6fhJD1Aa2YzLYulnsdQrGeHcVXDk3i+ckYgHyYacyRYA1VEIqQ1+xgp5SCEGLf7WrO0JOLL2xsJzW7aw/IQtmS12LyIzxKH8uEk72HnAxYTjaaVpHIanbupxjnDQwriyWE2EluDqcWPPRUAZbMlsVSRwHkw089fsm+GFYqj90+GETYJ+LKsYh9R5ddB7HhswspLCZLG+Dijv4SFnYb6/Pbxypd0IJeEQbNvzbsIqbqzhyFewXWdhTWxdgvC0jXcb5sKGC58lh7LHsZoWDlxTPLWSi6UfYxQOEed6e78IqCq19PjnvgQlEB1nBXzlEAwKZe88LX65fssEClZPaNl27AkT/9RfQF5ZpLZ9YKWS2HX/6nJ/DJH5RW9ixnNGy0OoMHreS9UyhC3vIhkpAlKuyCycJzim44chTuFVhn6Akwq+Hq+X+uNhSQhZ5Y+NIJq4CaiJql1rVyFEBexMzz5I6CUx81hYIQ8kVCyBwh5Jjj2McIIRcIIUesX292fO/DhJDThJBThJBbHMevIYQctb73aWLNuiCEeAkhX7OOP0UI2drcf2JjqHrlHAWQr3Tq8Yk1cxSAOXUWAHyy+Zji0FMsreL3H3iu7B26G/ne0WkspTU7F8GglCKe0fDLV4/iB7//ely3rR9A3oEB1XMUQD7cwi5iWs5w9FEYiKVV/L/ffKFic16nyJYJPdUjFKr976QlndIpRUdQFuApU33HHMW5hdIiASeFoSeHo5C4UHDqox5H8SUAt5Y5/ilK6X7r1/cAgBCyB8CdAPZaz/knQgh7Z34GwD0Adlm/2M+8G8ASpXQngE8B+GSD/5amYo/wqOAS2IWv1xF6qiYUDFnwwENKk9n/8cI0/v3IRTw/FVvNabeNrzw1AcDstnaSVnPQDYpev4RLN/bYs698koChsBceArsxrJiAdddc7ChU3bCb0hQ9h8PnlvC1w5N4+mzt0tN24qx6AgC/LNYVenIOOyx2Fc5O9mIGgqZbm7CrySqEnkRPha8FXvXEqYuaVzZK6U8B1PuJvB3AA5RShVJ6FsBpANcSQkYA9FBKD1JKKYD7AbzN8Zz7rK+/CeAmUm6yXpup1nAHmN3ZgJWjqBF6ckIIga9MSOLRE7MA8sMI3cxLswkcPr8EWfCUrNNk86/Y0DknY31+hLxi2cGJQD4OX85ROPso2MW33OKeTqLo5rwv9u8LSEJdVU+qQyiKXVJS1csmsgFTGAQPwblFFnpaoaMQPbyPglMXq8lR/B4h5AUrNNVnHRsFMOl4zJR1bNT6uvh4wXMopTqAZQADqzivplCt4Q7Id2f3+iV7EU+5SbPl8EsCso47ubSq44lXFgHkBcrNPPjcBQgegtv3b8JiSgWl+XNmQtFbRii2D4YwGPZW/LlBK0fBmswUzZnMzvdRZC2RnYq5axKvohkFd+z+FVY9AaWVT8msbr+/ivF4CPoCst1LUTFHIVQRCh564tRBo0LxGQA7AOwHMA3g76zj5a6qtMrxas8pgRByDyHkMCHk8Pz8/MrOeIVoNRzFtqEgRiN+7BnpWZGjAGA5ivwH9InTi444tbs/uJRS/ODYDF69fQA7N4Sg6kbB6Il4FaH40G2X4vPvOVDxZzNHwRrPmJiqzhyFZtjhmakllwmFbhRUFfnrzFE4XWTxCJOkUjn0BJglsmy7Yk8djoKHnjiN0JBQUEpnKaU5SqkB4PMArrW+NQVg3PHQMQAXreNjZY4XPIcQIgLoRYVQF6X0c5TSA5TSA0NDQ42cekUopQUbw+yqpwrJ7JBXxBMf+gW8ZudgzfLYYvyyUJCjeOzkrP2124Xi5bkkzi6kcMvlG+2O9KgjT1HNUQyFvfYU1HKUJLMrNNxlrOPuE4pcwXsgIAl19ctoOcMRdit8fMoarFgJNlgRqJajKG24A8xkNh8zzqmHhoTCyjkw3g6AVUR9G8CdViXTNphJ60OU0mkACULI9Vb+4T0AHnI85y7r63cAeIw6Yxlt4icvzeOGv3ncnvJqO4oK5bFO7PLYFYSenBeQkzMJbB8KAsh3hLuV7x+dASHALXuG7fLMxVS+UsvOUVS4u61GadVT3lFojumxLO5/wXVCYRSUnwZk8/+51ttZ1Q1bWIuXFyWqJLOBfEIbqNzIWNlR8NATpz5qdmYTQr4K4I0ABgkhUwD+DMAbCSH7YYaIzgH4HQCglL5ICPk6gOMAdADvo5SyK+J7YVZQ+QF83/oFAF8A8GVCyGmYTuLOZvzDVsrMchYGBaZjGYxG/PYqVKlCw50T1hdQf+jJU+Ao4hkNgyEvzsynXO8ofvjiDK7Z3IcNPT70W6Wx0VR9jqIWbC4Uu6tmjXeaXtiZzUR2IamYOxak8lVU7cbMUThDTyKo1UDor1DpBZiOoi8o4UIsU1L1lFL1io10QN5R+CWhYpjUWzGZLXCh4NRFTaGglL6rzOEvVHn8vQDuLXP8MIDLyxzPArij1nm0GpZEZBc6e9ZTHY6i1y9BFj2IWKGYWvgkoaC6JZHVsXs4DMD9VU+T0TR+5RozisgWEjmFIp7VQUjlMEg1RMEDn+TJJ7MLchRWMlszChK+U0sZ7NxQOZzVTopDT37LXaRVvYZQUIz4zdfSORiQUmqWx1YRin7LUVR7vcuN8ABY1RPPUXBqw2c9WbC7VCYULHZbbnFRMX5ZwHff/zqM9wdqPhYw7/6cu5/jWc2O9+sudxRKLl/Z0x8qIxQZDWGvWLZBrB6Csujoo8g3ojn7KJxhuwsxNwlFYdUT6wupladQdQMRa3+4czCgohvQDVozmQ1UF4riBLb9NW+449QJH+FhweLieUdhQPCQijX/xewaDtcdAnEms7Wcgaxm2ELh5tATpRRazrDvUIOyAFn0lISeyvVQ1ItzgqziSGAXhJ7UnF0y6qZeCudOagC2i6hV+aTlDPs1c1aQMcGsHnpiQlH5NWf/X5JA7AkBAA89ceqHC4UFC2fE0pZQ5GhdbqIRfGI+mc02mOWFwr2hJ92goDRfMkwIKSjPBEyhaCQ/wSgUivLTYzNaDlsGAxA9xFWVT4pWVPVkCUWtXgrVcmlBWUBa0W3htXdRVKl6qsdRsNyZTyy8kTGT2Tz0xKkNFwqLTFGOQsvRuiqeGsF0FOaFL+6oEhI8xC7LdSPsYu2MefcH5RJHsRqhCHnz+ZuyDXdaDmk1h6AsYlPE7yqhUPXShjugtlBoOQOy4EHAK+LxU3O45i8fwTPnl7Bgzf1i1WXlYPOeqlWZidbYGGdFFmA6Ci1HkXN5pR2n83ChsGAJ1LgtFEbFHorV4pNKHUWPX4LoIa52FPm8TaFQLBblKFbvKHIwDGo32al64VDAjJqDXxawKeLDRRd1Z5c03NW5pErLUUiC6ShemU+BUuDUTAIXraoy50DFYpiI1CoekEVPwbkBeeHgvRScWnChsCh2FLphFKydbCZ+SYCqG8gZFPGs+feFfSJkwePqHIVaxlEMBGUsFecoGuihYLDQkzN2Xi70FJAFDIV99l23G1D0XFEfhXnxruYocoZ5Ry8JHvvxAHAxlsG0JYJsR3s5mKOoJRReUSioeDKPeezz5nCqwaueLJijcIae6u2LWCnsA5vVckhk86EnUSCunvXE3I7scFr9QW9p6CmwitCTVfXkvHgVT4/NqDn4JRE9frFkem0nKZ71lM9RVB4MqDn6dSIBCcM9XlBqCgVbilUrUf1HN+/GDZdUn1RQ1lHYe7Pde3OyEuYSWWwIVxZVTuNwR2FR4ihaGHpiseusZu46Bsw7QtHtjkIv4yhCsn1hz2o5KI4u40Yo5yjMzmxTpLQcRVLR4Zc9GAx5kVR01+wfLzfrCaheHmvnfQQPPv5Le3H/b16HLQMBTMUymF7OVHUTjPfftAtXjEWqPka2elSc2I5iHUyQPTOfxHV/9SiOTK6NMf1rDS4UFnbVk8NRVJocu1p8jr3ZzmS2GXpys6Mon6MAgKfORPP/llUms1NqrqCklOUo2IUtntUQkEUMWdvznD0pnYLlVHxSqaMoVx57+FwUV//FI5iNm3kISfBg13AYl2wMYzTiN0NPy1mMVMlPrARvlRzFegg9XYxlQSkws+yenNV6goeeLIo7s7WcUXEkwmpxJjlZMjvkE83Qk4urnsolsy/f1Au/JOA9XzyEfaO9ACrvxK4HNu/JmSBnd+RhnwglqYJSU2wHw6ZILSSVupsdWwXL3xQM4BMrVz09fHwW0ZSKM/PmLgnna7op4sfMC9PIqDnsGelpyvl5JaGkO3w9hZ6KmzQ5zYULhQWLI6u6gaxmbmlrWejJ3pttIJ7VEPKaC2gkt4eeyiSz94314tBHb8KXnzyP//mDUwAam/PEYG7E6RKcjWcLVk4iIAv2QLwFF+QpWPjGmaPweAh8kqds6Onpc+aAZLb4ydmzsynih25QLKZUbKwj9FQPH7rt0pLGPXkdJbOLtyJymgsXCouU1e2bUHQsZ7SWOgpn6CmR1e2KlbVSHluc5A/7JPzuG3ciIAm493snsGUg2PDf0WsLRX4XN2vAc46y8EuCvQTJDZVP7GJb3KsQkMWSZHZGzeHo1DKAfIOnU3xH+/Lhpk29zQk93bC7NNm9nnIU7D1Sz1h3zsrhQgGzRFHVDWwdCCAxm8RyRjM7s1vWcFdY9cTKSd1YpTlDAAAgAElEQVTuKMo13Dn59dduw7uu21wSC18JTCjmLEcRlAV7mqzzjtgvC3ZX8oILchT5fdmF/3a/VLrl7shkzB46ucSEwiG+o468xEikdVU8+fJY977n6oWHnloLT2YjH3Zie7Bjaa3lDXcAS2bnHYXk8vLYcjmKYlYjEoBDKOLmxT/sk+zSZTbOHTAvwD5JQNgnFuQzOoXtKMRiRyGUhEMOn8vv5YrZoaf885yVTiNNchTlyOco1v5deIqHnloKFwrkk43sA7qc0aAZtKUNd4DlKJT8ED23l8c6SzlbBROKWSv0FPaJSFuOwpkkZxVFQyEv5l0QesqWyVEA5fdmHzoXxXi/KQB2jsLxvLBPsv+t9ZTHNkq+6sm977l6sR3FOhA9N8KFAk6hMD+8ZujJKGgsaya2o1BLHYWbhYLthJDrWObUKD0ljkK0k+jFOQoAGAx53RV6kqqHniileGFqGa/dMQhC8qGn4gGUo30B9PhEuwqsFaynHEXxfDBOc+FCgbxtdToKPUfrWlrUCIXlsZpDKDyuXoWaT2a3bqNc2CuCkHyOwtmVXJyjAIDBsNzRZHZGzeH1//Mxe+95udCTs49iLqFgOaPhspEeBGXRDj0Vu7SdG0LY0eI9G+sx9FRrpDunMXgyG/lKiQ09XhACu+qp1Z3ZaavqiSWzRc8aabhroaPweAh6fJK9h9s5w6jAUch5R/FEcrFl51OLhaSCyWgGj52cB1BOKERktHwT2MmZBADgko1hBL2Cw1EUPu8vb7/cdlKtYj2Fnuz1uetA9NwIdxTIh57CPhFhr4h4RoNmtK48ll1MllIqdIPad82uDz1VKI9tNr1+CdTSS6ejCDscRUAyvx4IerGc0To2AZXdZJyeMwWgpOqpyFG8xIRiOIygVyybzAaA3oCEIav8t1Wsp6qnBE9mtxQuFMivnwzIInoDEmJp1Qo9tebOmRACvyTY4ZUevyP05GKhyDuK1gsFo5Kj8Fklxqw7mzmQdpPfVGgqW3EfhZmjyPdRnJxJYEPYi76gjJBXzA9abKFLqwQT/PUgFCleHttSuFAg7ygCsoCIX0Yso5k7Alp4QfRJnoISUAAQBXc33CltdBTs7/E5/g+c5bFsJPegNe+pU1Nkiy9M5XIUzmT2S7MJXLIxbH+P0Sr3Wg1CyLrZcscb7loLFwrk+ygCsmhvbNMNA1KLHAUAy1HkS0ABQPKsjfLYVl/UmFB4JU/B31WQzLarnvLznjpB8YWpXOhJ0Q0Y1t6Jl2YT2D1sCoXz39MJoQCsdajr4C48X/XEhaIVcKFAoaMYCMpYTKrQ9NYtLgIAnyzg7II5EG5jj1ltJYnE9VVPgodAaKGAAvkSWa8oFHSBM0GVRY99Dsxl1Fo32iqKY+KlIzzyhQsT0TQU3bAdRdANQiEJaz70RCnloacWw6uekL/I+CXBdhQ52rqhgACwfTAIRTPwhzfvxmXWhFBxDTiKVoedAIejEAsdhV8WIHiI7SaA/IWYXSjaTYlQlKl6Asw82KmZOAAzkQ0UCkWlsSitZj2EnjJaDuz+ilc9tQYuFDBDT35JgMdD0B+S7XBCKy+Kn3/PAQBmnJjh9qonc7dz65OuTCh8kqfgAioLHnhFT0FsP79FrrOOghCA0tL3DAsvpdQcZq2cFBv65ww9tUOAy2EKhXvfc/WQdNwk8Kqn1sBDTzAvMuyCwwbNAWhZwx1gCoRTJABW9eTe0JOiG5BXOcupHnorhJ4kSyicjiJoX4g75SjMi+z2wSC8oqfk/9TpeJzj0gEgKDtDT+2vegLM13it5yhYD0UkIPGGuxbBhQKWUHjND3R/MF+73srQUzlEqzObUneIRVLRcWI6bv9Za+FYEycFjkJwCgWBVyxcwOMVPfCQznXkMvd55VikoJSXYTsKRUdK0SF6iB2eClrvOULQ8rxPJbzS2g89Ja3lX4MhL7Jr3B25FS4UMENPrIGr3+Eo2n2Xx6qs3FIi+6cPHcPt//gEotZ0VlU32hJLr+goRA+8UqGjIIQgKIv2XWW7YaGOD912Kb5w16tKvh9wOJ6UoiPoFW3XwUREEkqdSLtYT6GngaAM1aow4zQXLhQwHQW7S+1vU+ipHKxvww3rUCejaTx05CLUnIGHX5wB0Nr1sE6cjsL597EcRfFKT3NCa+dCT7LgwYYeH64cj5R8P2S5hqSSQ0LRC/ISTEQ6lZ8ArNDTGhcKVsjAemrW+r/HjXChgCkUQW+pULS6A7kYsUOO4gs/P4v3f/W5gmOf/9kZeAgw3OPFd49OA3CBoxA8uPGSDXjDrsJtbUGviFQHk9k+qfJr4qx6Mh1FXuSYiHQqPwGY3e5umL67Glh+ivXU8IR28+FVTzCFoi9gvsl6fKJVfURb2nBXDnZRbHfl01NnFvHjU/PIGRSCh0DRc/ja05N42/5RDIW9+OxPzyCaUqG22VF4HTkKjxXH//CbLyt5fEAW7DEs7cYUisoJfpZsTyo6UkquwFGwZHaneigA4Nqt/fjuC9M4t5DC1sHGV9h2kkS20FFktBz6OnlC6xDuKGDlKGSWWCS2aLSy4a4cLNTV7sqn5YwGNWfgYixj/1nRDVwxHsFbrhhBzqB4+MWZtjmKsM8cNe4TBXsGUrWLaVAWO1Yem6klFI7y3aSVo7C/5+28ULBd2j99eb5j57BaWOhpwBIK7iiaT9cLBaUUc3GlYFInCz+1OyTAqqza7Sji1h3ZGatTnFWRhL0i9oz0wCd58Mp8sm0Ndx4PwVDIi0hAsndfVPt7A95O5ihyBcn1YkQrr8KqnpyOgn3dqWY7ANg6GMSWgQB+cmptCwUhQF/AdKK8O7v5dL1QRFMqMlquYKH9QIgJRXtfHnYxbLtQZMydCGfnkwDyVSSsQifsk5BU9LY13AHAV377erz3jTvs3RfV8kUBWehgjsKomqMATEFIWn0U5R1F53IUgOkqDp5ZXLNlskklh5As2kUOvDu7+XS9UFywwi2sWxbI91K0asx4JZijaPe8p2UmFMxRFDWGhb0iElm9baEnwNzwFgnItnhWu5gGZLFjOYpaoSeAOR4z9OQ2RwEAb9g1hLSawzPnlzp6Ho2SVDQEvaL9/8BDT82HC8WSJRRORxHsjKNgOYp2Ogo9Z9jCUBJ6shrIQj7RchTtSWY7kWyhqJaj6JyjUOoQiqAsWsnsQqHwSWazYCdzFABw+WgvAODcQrqj59EoKSWHkI8LRSvhQmE5irECR8GS2e2uemp/eSzLTwCVHUXIKyKZ1a0RHu19y7Au5uo5CrFjndn1hJ6CXhGLSQUGLRwESAhB0Ct2XChYqHV+jZbJspAe+3/gOYrmU/MdSgj5IiFkjhByzHGsnxDyCCHkZev3Psf3PkwIOU0IOUUIucVx/BpCyFHre58mVisqIcRLCPmadfwpQsjW5v4TqzO1lEFQFgq2qtlC0eaGu3zVU/ve6Cw/sWUggAuxDLJaLi8UlqMIOxxFu5vDmDDVchRqzujIOtRMjWQ2YIoD22YY8hY+NiiLHW24A8zXtj8oYz6Z7eh5NEo0paIvIMEnckfRKup5h34JwK1Fxz4E4FFK6S4Aj1p/BiFkD4A7Aey1nvNPhBD2yfgMgHsA7LJ+sZ95N4AlSulOAJ8C8MlG/zGNcCGWwVhfoGCEwkDHq57a5yhYfmL/eASUAhPRtF2XnncUkpmjyLXfUdihpyqrQv1WP0InXEWtPgrAFDK2zdDpKMw/Cx1PZgPAUMi7Zh3FxVgGI71++/8hpebwvx8/jeW01uEzWz/U/NRTSn8KIFp0+HYA91lf3wfgbY7jD1BKFUrpWQCnAVxLCBkB0EMpPUjNiXf3Fz2H/axvAriJtHHwzYWlTEEiGwBes3MQ775+ix27bRedqHpyCgUAnJlPIano1gA+83zCPhGJrAZNb3+Ool5HAXRmgmw9yeygV4Rq/Z8WC8U7rhnHbZePtOz86mUovDaFIqvlsJhSsanXZzu7Q2ej+JsfnsLDx2c6fHbrh0Y7s4cppdMAQCmdJoRssI6PAnjS8bgp65hmfV18nD1n0vpZOiFkGcAAgIUGz21FTC2lcc2Wwj7OXr+Ev3jb5e346wtgDX7tnPXEhGLvJlMUL8YySGbNpKtzeF1S0SF4SNsdBas8qyYUbGZSJ3opFM2oy1EwQkVC8d437mjJea2UobAXh8+nOn0aK2Zm2QyXjUT89nZBtiCKDbPkrJ5mf+rLOQFa5Xi155T+cELuIYQcJoQcnp9ffYNQIqshntVLHEWnYBdFVW9/6Gm83w/BQxBNqWZ1jmNkdsgnwqBscVF7hYIQU5yqxfFtR9HmCbI5g0LN1ZfMZhQLhVsYDMmYTyiuGXGv6gYmo7WrsC4um8Uom3p91j4Q0xUDQDTNhaJZNPqpn7XCSbB+n7OOTwEYdzxuDMBF6/hYmeMFzyGEiAB6URrqAgBQSj9HKT1AKT0wNDRU7iErwu6hiLhDKOQOTI+NZ02h6AvI6AtIWEypSCh6wVId58WteNVnO/AKnpp9FED7t9yxpGk9yexyX7uJobAXWc0o2BbXSb5+eBI3f+onNfNO07G8oyDEDJeyPqRokgtFs2j0U/9tAHdZX98F4CHH8TutSqZtMJPWh6wwVYIQcr2Vf3hP0XPYz3oHgMdom25r7B4KlzmKds56Ws5okEUPfJKAvoCMaEpBMqsXLOFxft2JxKskeqqGvPLrUNt7kWNCsZrQk1tgI2zckqeYiKaR1Qzb8VZi2nIUI70+AIX/F0vcUTSNespjvwrgIIBLCCFThJC7Afw1gJsJIS8DuNn6MyilLwL4OoDjAH4A4H2UUnZL8F4A/wwzwf0KgO9bx78AYIAQchrAH8KqoGoHrG9gc3+gXX9lVVhYR21zeSwrDe4Pyoim1JIO4kKhaL+jkAVP9WS2N1/t0k4ytlBUf00CBY6i9atkG2EoZF5o3SIUbPR5rQKFi8tZ9AdlWyCc7m6R5yiaRs3bG0rpuyp866YKj78XwL1ljh8GUJIhppRmAdxR6zxawXMTMYxG/PZ44k7DLobtdhQ9lhAMhGScmknAoCgYOR3y5ntMOjFuYrjHiw1hX8XvO3c+tBPW2FXLUZQbLe42mKNYcEm4ZsG6yKdq/J9OxzK2mwCKHAUXiqbhzndtm3jm/BKu3dbf6dOwyc96am/VU7GjEDyesjOJgM44ivvvvq5qboRdfJfSGj7+nRfx3jfuqCoszaLu0JP1+gVlAZ4O7cauRT705I6mO+YoauVMppezGOvLRwTY+2Rjj487iibSdSM8Hj0xi9+67zCmltKYiWdx9ebS9ZWdQrI6s9vZYVwoFF7EMhriWa1ijqITyexev1T1Ysymhv7kpTn8yxPn8NOX2lJZXX8y2zo/tyayASDilyB4COaT7gg9Laas0FOFSrb3feVZfPAbz+NiLINNkVJHcdXmiD3IkrN63PvObRFpNYcfnZi15ypdvcU9u7BY93E7p8fGMzp2DoUAmB3plJpC5SZHUQtZNKuinp2IATDLnttBvaGnoN3h7t6Pm8dD7BLZTkMpxWKycuhJ1Q08cnzWFoGR3nwxCssX7R+P4PvHZhBLq9jQ03p3ud5x36e+xbx53wi2DwXxvaMz8EkeXDbS0+lTsunErCeno+hz7AsvEAqHo+j0XKJKBGTRvnDEM+3JVWTrTGaz0JibHQXgnu7s5Yxm3yyVCz2dmI4XOAWno/BLAobCXoxbBSo8/NQc3PmpbyGCh+B9b9wJALhiLOKqO2RWeqq2KZltGBTxrIYeSygGnEJRVOnELobVFgh1EmcJarscRabuPgqh4He3MhTyuiL05Eyol3MUz02YezPe/wvm53jrQL7w4rZ9I7jr1VvsdcY8od0c3H2L0yJu378JXzk0gbfs6/yMHSeEEIge0hZHkVFzmIlnQSkKktmMcNHdb8grIasp7nUUjvONty30tLJktrN6zI1sCPtw9EK806eBBYdYlSt5fm4yho09Pvzhzbtxy96N2LspHxV45wGz3/el2QQA7iiahTs/9S1GFDz4t/e+Bne9ZmunT6UEUSBtyVH88b8fw41/+2MAKOsoisMkLKEtV5ni2kkKHUVzQk8zy1nc+bmDJeEYw6C47z/P2Xff3hqhJ6/ogeAhJSPG3cZ4vx8LSaVjuz0YizUdRQxXbY6AEILLR3tRboao7Sh4011T6EqhcDOS4GlLpcbFWAZDYS9ev2sQ11klwn0VQk+AY22n4M6LHat82rUh1DRHcfTCMp48E8XBM4sFx0/MxPFn334RX3rinPl313AUhBD0+iVEAnLVx3UaFtefXOrspjvmKGTBUyIUC0kFE9E0rqpRrdgXMG9+Fl3SF7LW6crQk5uRBE9b+iiSio59o7344q+/quDv7vGJiGf1MqEn88/V9kJ0kr6AjO2DQYz2+Zt2cWA5iJdmEsCV+eOzcbPXgC0jqhV6AoDPvvuagi2KboRNKJhYTGP3cLhj57GYVOAhwEjEV5LMZnu9949Xr1YUBQ96/RJ3FE2COwqXYeYoWh96Kh7TwRiwutRLHIX1Zzcl/5185M2X4bPvvgY9PqlpyeyMNT7i5Eyi4PhsPB+KEj2krtfkVVv7C8o43YhbHMV8UkV/UEbYJ5Y4ikeOzyLsE+39KdUYCMo8R9EkuKNwGZLgacusp0RWLxEDwLTsZ1Fa82/nKFwqFOwiF7YcUTNgsXqWGGUwRwHU5ybWCgNBGQFZwEQd471byWJSwUDQi6AsFjTcaTmzf+Lmy4brGiXTF5R51VOTcOenvouRhPY4ikRWKwkvAWZ3NlA6k4g9thOznlZCj990FM0YQJyxGuomoumCO9vZeBZ9AQlhn1izh2ItQQjB5v4AJqOZjp7HQlLBYFi2F2YxnjyziOWMhlsv31jXz2EjaTirZ/28y9cJYhtyFKpuQCnqvmYMhswPaPFMopDLHQWjxydBy1G7a3o1ZByTS1+eS9pfz8YVbIr4cdOlG1wzULJZjPUF6loY1EoWU6rpKLxiwfTY7x+bQUAW8Ibd9e2i2djjw9RSBkYbJx2sV9z9qe9CzKqn1r6x2d1xudDTu1+9BR/7pb0lxweCXoge4vpQCwuRNaPyiSWzgfx6TcB0FMM9PvzF2y7Hfb957ar/HjexuT+AiWi6o5vuFhIKBkKyKRSO0NMTpxfwup2Ddb8H9432IqnoOLOw9la8ug0uFC5DEkjLHQWz8+Ucxd5NvXjHNWMlx3/1VeP4xn9/tV2G6lZYT0gzEtppNYe+gASf5MGpmUJHMdzjRdgnYXidzRHa3O9HRst1LAmc1XJIqTkMhrwIeYWCkN98QrFzUfVwpZXwfmEq1vTz7Da4ULiMdlQ9sYa0cBlHUYmgV8RVm90zQLES7N+03IR5Txkth6BXxCXDYTx1dhGGQaHlDCymlLaMMe8E7ELcqYQ2K2ftC5iOIqPlkDMosloOaTVXMD2gFjs3hBCQBbwwtdyq0+0auFC4jHZUPeUdhbtHSjRCj695jiKj5uCXBPy367fgxYtxfOOZSSwkFVCKdeckGKyXopV5CkppxbwBSz73ByXb8aZU3XY4AysQCsFjdm4fmeSOYrVwoXAZkuBp+aynpGJeRMvlKNY6PXaOojmOIiALeMc1Y7h2az8+8f2Tdk/FcM/6SmIzWB9NK8tK/+7hl/DOzx4s+72llPneZI4CMHNq0SQTkJV1t1851ovjRdNmOSuHC4XLEAUCrU2hJzfvR2gUlqOIZ5qTo/BJAggh+PjtexFLa/hfj74MYP06ioCVg0prrZv3dG4xVdLEyIim84LgFAq2yGggtEKhGI9A1Y2SXhjOyuBC4TK2DgRxfDqO51tol1noaSU5irVCPvS0ekeRtRwFAFw20oNrtvTZy5HWq1B4RQ8IAbItHAyY1QwkFd2evuuEOZm+oGwPekwqOUdIamVObt9oLwDg+MXOT8Vdy3ChcBl/cPNuDIe9+IOvHUFabc0CnuQ6dhQ+yQPRQ5pSHptWcwVVXqwaTPCQFcXK1xKEEPglAekWCoWimz+7XDNcNKWCEHM1a0HoKdVY6ImF0pab4DC7GS4ULqPXL+Het+/DmYUUfnxqviV/R1LRQUg+zLCeIITY3dmrJWOFnhhvvWIEPsmDDWFvSUPieiIgCy0NPTEnUU4oltIqev0SRMGTT2YrZjJbEoidg6oXNtk31aKbrm5h/d1SrgP2WItYmjH5MqPmQEjhTKJE1hwIWG6O/3og7BObsg414wg9mT9Xwl2v3rruJ5L6ZaGloSfFSiwvlNmmF02p6LfGsQcdVU/RpIq+gLzi96zgab1D6ga4ULgQdieVbEKc/bfvP4zhHh/+7p35OdlJpXSM+Hqixyc1pzPbKo918uE3X7bqn+t2Wn1hreUo2F4Utjo2qZgNgCsNOzGCRTOjOCtn/V4t1jABWYCHNCche3YhVbKhLVlhcux6occvrvq1MwyKjJaDX16/r1Ml/LLY4tCT6SjK5yg0jEbMceyhghyFsuKKJ0bQKyDNhWJV8ByFCyGElEzObJTljFYyu6fSLor1Qtgrrbo8loVHam2vW48EpFaHnsyfvZBUcehsFH//8Cm7AW8ppaI/aFau+SXzhonlKFZa8cQIyGLZ3duc+uFC4VLCTQifaDmzDDGj5bDg2PqWUHSEfOuvK5sR9q3eUbCKs/WY8K+FXxaQ1lp3B553FAr+9anz+PRjp/E3D58CpRRRR+iJEIKgbP5fRpNqw5VmxTOjOCtn/d5WrnHCPnHVOQrnXfVENI2hsHlHlsxqGIu4e9vaaigeT90IbHJsNzoKv9y+HMVMPAvBQ/CZH7+CSzeGoeqGncwGgN0bw3jyzCISit5wjiIgi4jx8thVwR2FSwn7Vh96cn44nLN71nvoKWjdQa5mVDbbbuf2abmtYLWhJ1U38OWD58qOoqGU2mG9+aSKM/Mp/Nq1mzHc48UXnzgHALajAIBfuHSD3cXdeDKb5yhWCxcKlxLyrj58EksXOgrGek9mB70iDJrPMzRC1zuKVSSzf356Hn/y0It46my05HvO/5OXZxNIqzns3hjG63cN2dMInCGmGy/ZYH/daOjJXKnKhWI1cKFwKWGftGpHURx6AoCcQZFSc+vbUViVSqt5/dJd7ChWG3patPJhF2KlK1UVKz8heIj9d+wYDBZsrXM6istGwhjpNcelrKY8liezVwcXCpcS8omr7i6OZcwP7GDIawsFi92vxzlPDNaolVYavzjYjqILhSIgiVB1A7kGV4iyhsTpWLbke6ziaaNjVtb2oRBet3MQrJfOmaMghODGS01X0Wh5bEAWWjYOp1vgQuFSmlG5w0JPV4712jmK9TzniZEfJtf462fnKLoy9GReFjINhp+i1qjwi2UcBat42hQxhSIoCxju8aI/KNsD/PqKnMN/vW4z3nTZ8Iq22zkJekVoOWqLFGflcKFwKWGvCEU3VjVHnwnF3tFezMSzyGq5/NKibnAUq7iLZELRneWxq3v92ATYi8sZGAbFQ0cuQLMS21nrYr3JqrrbPhSyx3LcevlGDIbkknlOezf14p/vOgCv2Nj/BbtxWI3D7Ha4ULiUsNXnsJq74uWMhrBPxLbBACgFppYy9nydXv/67aNgQrGqHEUXJ7MD1r8502Bcn4WeLsYyOHhmER944AgeOzkHIJ+jYEKxYyhoP+933rADP/7gjU2fQRZowvuh21m/t5VrHBYaSmS1hpN4yxkNkYCEXRvCAIAT03HMLJtx48tGeppzoi6EzQhKreIOMtvlyWyg8dCTnaNYzuLYBXNf9bmFFIDyjoIheEhLQqIh22FyR9EoXChcCks2ryZPEbNGNl+yMQyf5MGzE0uYTygYjfgxGFqfqzyBfNXTapru0l2do7BCNQ1eWNkMp7Saw5NnFgHkq+5Ys93OoRDede043nLFyGpPtyYsfMhHjTfOqkJPhJBzhJCjhJAjhJDD1rF+QsgjhJCXrd/7HI//MCHkNCHkFCHkFsfxa6yfc5oQ8mmyXudfr4BQM4QioyHilyEJHlw5FsGz55fwwtQyrhjrbdZpuhLnwptGyWg5yIIHotB90dnVh57yLviJVwqFgoWegl4Bn/jlK7DD4ShaRTPeD91OMz4FN1JK91NKD1h//hCARymluwA8av0ZhJA9AO4EsBfArQD+iRDCbtc+A+AeALusX7c24bzWND1NylH0Bsyfc/WWPhy7GMdENI0rxiJNOUe3wkJPqwk1ZLUcfFL3iQTgCD018PoZBkUsrWKvtVOFFWPYjsIKPfna6NRsh8mT2Q3Tik/C7QDus76+D8DbHMcfoJQqlNKzAE4DuJYQMgKgh1J6kJozF+53PKdrceYoGmU5rSFiJa2v3txn18Vfuc4dhSyY61BX13CnI9CFI8aBfKimke7seFaDQc1KJca2wSAuLGWg5wy7PNbXYAVTI+RvHLijaJTVCgUF8DAh5BlCyD3WsWFK6TQAWL+zHvxRAJOO505Zx0atr4uPdzUsR9HoxY5SilhGs6ubrtqcdxGXr3OhIIQg6BVXNd8noxldmcgG8uWxmQYurCw/sXs4BEnIl73qBsX0ctbuZfC20a0F5PUTelL0HP7ga0dw0ArptYvV/m+9llJ6NYDbALyPEPKGKo8tl3egVY6X/gBC7iGEHCaEHJ6fb80+abfQSI7ixYvL+I1/OQRFzyGl5pAzKCJW6Gkw5MWWgQC2DwbtsNZ6xtznsYrObFVva3jETfhXkaNgFU8DIS829voQlAW8bucgADP81AlHYS9AWgdVTxeWMnjwuQuYXi5tZmwlq/LWlNKL1u9zhJAHAVwLYJYQMkIpnbbCSnPWw6cAjDuePgbgonV8rMzxcn/f5wB8DgAOHDjQ+GjQNYBXFCCLnhUJxROnF/D4qXlcjGXtu7mIP19a+5EuWOPJCMir20FQvC+7m1hN6Il1ZfcFJOwYCmG8z8DWQbNX4vxi2q56aqej8EkeEIJ1MUH2vJXr2dxgly+eLIoAABhHSURBVHqjNCwUhJAgAA+lNGF9/YsA/hzAtwHcBeCvrd8fsp7ybQBfIYT8PYBNMJPWhyilOUJIghByPYCnALwHwP9q9LzWE2HvyuY9LVq2P5ZWIVnVOiyZDQC37N3Y3BN0MavdSZFWc3YStNvwiuaFtSFHYb0H+wIyPvXO/aAwmzslgWAimoZs3cB4xfYJBVuAtBqH6RYmFi2hGFgjQgFgGMCDViWrCOArlNIfEEKeBvB1QsjdACYA3AEAlNIXCSFfB3AcgA7gfZRS9j/3XgBfAuAH8H3rV9ez0p0UbGpnLKNBZkKxjjuwqxFc5VazZFbHhvD67TWpBiEEfklYVeipPyjbZakAMNYXwGQ0jbE+vyVE7a2AD3rXx2DAiWgafknAUJv7oBoWCkrpGQBXljm+COCmCs+5F8C9ZY4fBnB5o+eyXgmtcDAgSyQup7W8o+hWoZBFWzgbYS6h4PrtA008o7VFoMGdFNG0Cln0lITtxvsDmIimMRiSO5L7MR3F2hIKVTewmFLQ65fshPz5xTQ29wfaLrTd6a3XCGGvtKJ1qM7Qk+Ax30iNjmZe65jJ7MYuDFkth+WM1rWOAjB7KRoNPfUH5JIL2XDYi1MzcWQ1oyP9KQFva9e7NpuXZhO44/8cxHJGw1WbI3jwd18LwNxU2egU3dXQnR1Fa4SQT0R8BTmKaMoc+Lec0W3R6At0p1Cs5sIwnzBfx2HHzoRuY6Whp4yaw78+dR7nFtIlY8IBYDDsxWJSRUbLNTwFdjWstS13n/3JGai6gVv2DuPIZAzLGQ2UUkxE09jS5vwEwIXC1UT8kh3zrYd8jkLFYlJFJCDZIahuI7gKRzGXMAcnDvV0s6MQVxR6+s4LF/HRB4/h0Lko+gKl4c7BkBe6QTGXyHbEUay2uKGdzCWy+M7zF3HHgTHc9ZqtoBR45nwU8wkFGS3X9oongIeeXM22oSC+8cwUkopec6pmVsvZd9DLaQ2KbjQ8dXY9EJTNLW1azlixWM7FTUfRzaGngCSsqOHulbkkZMGDO68dx4Gt/SXfH7RCoFNLmY68LwOysCb2UTw/GcMXfn4WmmHgN167DRt7fJAEgkNnl+z+p3ZXPAFcKFwNG5h2Zj6JsE/CK3NJvGnPcNnHslATYFY9pRQdg8HuvdA516H2BuoTivmEAr8sYI6HnqzXof6w5yvzKWwdDODPby9fk8KqdKaXs/YO7HYSWgOO4pnzUfzKZw4CAH7tus3YZvWf7BvtxaGzi9i1wbwecEfBKYAJxSvzSTx+ch7fPTqNQx+5CQNlSuOijgqfWFpFIqu3ZTKnW2Fbzb5/bBqHzkXx9qtGrb3MlatF/us/P4mrxvswEJIhekjB7uZuwy+vLMdzZj6JSzaGK35/0HJnOYN2pOop5BURz7hbKNjgxG//3msLBne+als/vvjzszg5EwchwFifv+3n1p0B7DXCloEARA/BK3MpPDuxhJxB8d2j02Ufu2glskd6fYhlNCym1K6teALyjuLvH3kJ33r2At79hUP4+uHJio/XcgZOzyXx7MQS5hIKBkNeeDzdO+0+sIJktpYzMBFNY7tjW10xzv0nnUhm94dkZLRcw6PT2wHLMW4ZKHwdr93aDy1H8fmfncWW/kBHXj8uFC5GEjzYPBDAk2cWMbVkznZ58LkLZR/L3mQ7hkJYSqlYSqtlnUe3wHI6cwkFv/HarSAE9mtYjouxDAwKnFlIYTKaxnAXJ7IBs+JuKa3WdWE9v5iGbtCqDjbil+yS7U4ks1kYlt1QuZH5pAJZ8JTsDH/NjkG8/apRfPCWS/Cvv319R86NC4XL2TEUwuHzSwCAm/cM47mJGM4vpkoex5rttg8FsZTWQCkw0MXJbGfD1617NyLkrd68OBk1RSRnUDw3EcNQuHvzEwDw5n0jyGoGHnh6ouZjX5lPAkBVofB4iP1+7IijsP7u1TRhtprFpIrBUGkPil8W8Klf3Y/33bgTo5H2h50ALhSuh334JIHgw7ddCgB45PhsyeMWUyokgRTEL3noyZyXdfWWPvT4pOpCsZS2v1ZzBjZ0uaN41dZ+XLu1H5/76Rl7+VAlzsybNy7VQk9APvzUCUfBPgtudhQLScW1UQAuFC5nh/Xh27OpF9uHQgh5xbIhlGhKwUDQWzAttqvLYy2heP3uQUiCB2Ff9QGLE9E0RA+xnchwlzsKAPjdG3dgejmL7x4tO8zZ5pX5JDaEvQjXGF/PEtqdSGYzkXKzo1hIKnYZsdvgQuFydlglcVdbi4c29vows5wtedxiUkV/UC6YFjvo0ruTdjAU9qLXL+GXrjR3YIVrzM2ajKYx2ufHpVblTrc7CgC4YfcQ/JKAYxfiFR/z8mwCz04s1XQTQL6Xop2TYxl26CnlXqFYTLo3r8jLY13OpRvDuHI8grfsGwFgVjVNx8sIhVXlFHEMAexmRxHyinjuT262K5fCPsnuuC7H5FIG430BbBkI4NmJWFc32zEIMUOZk9F02e8/emIWd993GIKH4Ldet73mz8uHntrvKAKyAJ/kwWLSnaEnSqmVo3Dn+44LhcsJyCIeet9r7T9v7PHhpdnS7X7RlIqtAwFErNp/Qrp3zhPDWd4a8op4Zb66o7hl70bs2dQDoLub7Zxstqa+luP0nJnE/skH34ixvtpNYMxRdCJHQQjBQNDr2tBTPKtDzRmuDT1xoVhjjPT6MJdQCkZTqLqBmXgWwz0+e/Vpf0C2yxE51m6PCqGnlKIjmlIx3u/H7ftHoeoG9oz0tPkM3cl4fwBPnY2CUlpSjRNNmSPF663E6aSjAMyEtltDTwuW03Gro+A5ijXGxl4/KM1POAWAkzNxqLqBK8Yi9v6Jbg47lSNcpeqJVTyN9wUQ8or4jddu6+pmOydjfX4kFR2xdGkhQLTCSPFKsItgJ3IUgFku3q6qp6yWA6X1b2teSHCh4DQRNidn2pHQfm4iBgC4anMEPkmAV/R0dWlsOcI+EWrOsHc2Ozm/2Jk9xGsBtvvAWT7MWEqrK7ohYeG8YI0Bl61iINSe0FNWy+HVn3gU33hmqu7nMKfj1s8tF4o1xkZLKGYKhGIJwz1eW0T6AjIGunggYDlYt2uxq4ilVfzND08h7BPtCjNOnnEr98AaEp1EUysTit3DIfzDnfvxpsvKD7ZsNaajUFd0p98IF2IZLKU1HLuwXPdz3B564jmKNcamXjMePL2cwX+eXsDODSE8NxnD/vGIHQL4k7fuwWgHBoe5GVbjn8hqGHJUNL33/z6LiWga9//mtTVHuXcj4/3m+6i8o9AwWkcSm0EIwe37R5t2bitlICRD1Q0kFb1mz8dqYH1OF6qMjClmIamCEPeGjPknY43R4xfhlwQ8O7GEe793AtsGgji/mMa7rt1sP+YtV4x08AzdCRMBp6M4PZfEwTOL+MibL+3q/djVCPsk9AWkspVPi0kF/WWWFLkV5rIXk2pLhYIJxIWY+fsXfn4Wr9kxgMuKCiQOnY1iICRjx1AIC0nF1QUoPPS0xiCEYKTXh+8fmwG1htgBwFXjkRrP7G7CZUJPPzhmTuJlTXmc8oz3B0p6KbScgXhWL7v21K30h9rTdDdlua8LSxksZzT8xX8cx1cPFc7MMgyK377/MP72h6cAmMlst4adAC4Ua5KNvT5QCrxqax8+cNMuRAIS9o31dvq0XA27g0wq+eqd7x2dwdWbI3beh1Oe8b5AydgYVgXl1lBJOewJsi1uumNOIqHoeHbCHOh5MVbY7HliJo7ljIaz1o2eOefJva8lF4o1CLuwvf2qMfzBzbvx9EffhIDMo4jVYI4ibjmKicU0jk/H8eZ9PExXi/H+AC4sZZAz8klgtst9LQnFQNscRV5UHz85B8DMKTp58kwUgDljjFKKiWjaLhxwI1wo1iA7hkIIyII91mOlO6G7kR47mW0KxWMnzQm8t+zd2LFzWitsGwxAzRkFyVk21n4tbQEcCMkgBJgtMwKnmUwtpbF72Kyge/wUE4rCv/PJM4sAgLSaw+m5JBaSKna6uOqOX2HWIHe/bhse/aMbCgYAcqoTsnMUZsjk3GIaYZ9o9wlwKuNcyctgQrGWchReUcBIj8/um2kFip7DXELBddvM4ghWVhxN5ZdAGQbFobNRu5z9RydMMeFCwWkqPknASC8vf10JgjVCnDmKqaV0XfOJONWFYi2FngBzzei5Mou/msV0LAtKgX1jvfZMK/Y7Cz+dnElgOaPhjmvGAOTdLRcKDscFOHdSTC1lOrYtbK3RF5TRF5Dwynz+ArtkCUVkjbnarYNBnFtonVCwRPZ4XwCbrPfXa3YMAsgntFnY6VeuGYOHAM+cX4JfElz9fuRCweka2LwnSimmljIF2wA51dkxFCp0FGkVYa/YkbWmq2HbYABLaQ3LZWZXrQZKKR46cgE/fdmc7DzW57cd6w27hwAAFy1H8eSZRWzuD2DLQBCbIn4YFNixIejq+WK8VIbTNYR9IpKKjuWMhqSic6FYAduHgnjs5DzmElk8cXoB0ZS6pvITjC0D5oKlc4spXBloXu/R0QvL+MADRwAAHmJWJjKH8AYmFLEMDIPiqbNR3LJ32Dofs/R4Z5V9426ACwWnawj7JCxnNLt8keco6mfHUAhfPzyFD/3bUTx2cg6RgGRfdNcS2wYdQtHEJtUfHZ+FhwC//frt8HgIJMGDN+waxNmFJLYOBDAU9mI6lrXzE2wSwJaBIJ44vejq/ATAhYLTRYR9IqaW0nbnLHcU9cMS2o9ZfQGxtIarxtdWfgIwJwQTArvRrVk8fHwWB7b048Nvvsw+dtu+EdxmlbBv6vXh4nIGT5018xPXMaGwqu52bgg39XyaDc9RcLqGsNfcm80chZsbnNwG24ktCQQftS6GazH05JOaUyL75YPn8MkfnARgbkc8OZPAm/ZsqPj4TRE/LsYyePLMIsb7/XZY6qrNffCKHlzh8skK3FFwuoawT0Q8o+H8Yhphr4geP3/718vmfnOp01uvGMHdr9uGU7MJ/MKllS+MbmbrYNB2FNPLGTx+ch7vunYcWo7i+HQc+2uEpM4upPDn/3EcWo7iyrFeu1fi5j2VmzdHev14+Pgszi2m8c4DY/bxa7f149jHb3F90yz/pHC6hqs39+HzPzuLB5+7gLE+f92b2TiAKHjwH+9/HTb2+uDxEPztHVd2+pQaZutgEN87Og1KKT77kzP40n+ew1WbI3ji9AL+8rsn8A937q86Dv2vvncCsuDB9sEAPviNF5BSdewfj9j5j3JcPmpOjr3zVeP44C2XFHzP7SIB8NATp4u4Ze9G7N3UY1U88bDTStk6GOzYvutmcsVoL2JpDS9ejOPH1oiNR47P4rtHzWnCH/nW0YJSYCfPT8bwyPFZ/O6NO/F37zTF8o5rxvHlu6+t+ne+/apRHPvYLbj37fsQWUNjTxhcKDhdg8dD7Ls5nsjuXm7ZuxGih+AfHzuNc4tpEAJ845lJPDcRw3+7fjNk0YOPPngUAPCnDx3Dn3/nuD0Q8VvPTkEWPXj3q7fg8tFePP9nv4hPvuOKmvstCCHwy2tXZHnoidNV3LB7CH/8lsvw+l1DnT4VTofoC8p43a5B/ODFGQBmOOirhyYBAL/+mm3YMRTCx79zHPd+9zjuP3geABDLqPjEL+/Df7wwjZsvG7aHTLq5Sa6ZuMZREEJuJYScIoScJoR8qNPnw1mfEELwW6/fjks2ursckdNa/ssVmwCY1Vx3v24bAGDXhhB2bgjhXdduxkivD5//2VlsHwriAzftwreevYA7/s9BLKZU3L5/UydPvSO4QigIIQKA/w3gNgB7ALyLELKns2fF4XDWKzfvHUZAFnDznmHsGArhlr3D+E1LMHySgN9/0y54CPDxX9qLP7h5N/6fX9yNF6aW0euX8MZL1ma112oglNLaj2r1SRDyagAfo5TeYv35wwBAKf1EpeccOHCAHj58uE1nyOFw1htTS2kMhrwVE/SLSQUDjvWk3zg8Cb8s4K1XrG1HQQh5hlJ6YCXPcUuOYhTApOPPUwCu69C5cDicLqBW5dtA0Q7rOw6Mt/J0XI0rQk8AymWESqwOIeQeQshhQsjh+f+/vfuNkasq4zj+/YUFDVCg9I+BgFYTApZGqDQRBCGR8ALfaIKJbYit+Ab8E/GdYEz0TV/QICHQF7WBGkBDqkFjqwgBoo2AgK00lHVjSw2RYiM01tKWSCA+vDhnwrDZvezMnJl7787vk9zMzJm7557z7N155tydOef110fQLDMza0qiOAB0p+tzgH9N3ykiNkfEqohYtWSJP7ViZjYKTUkUfwHOk/RxSScBq4FtNbfJzMxoyP8oIuIdSd8GHgVOALZExGTNzTIzMxqSKAAi4mHg4brbYWZm79eUS09mZtZQThRmZlbJicLMzCo14pvZ/ZB0FPj7B+y2GDhU8LCnA0fGqD4oG8OS7RtGX4dR7zidg45ds+qbHr/O449FRG/fL4iIVm7AzhL79HjMzeNUX+kYlmzfMPo6pN/J2JyDjl3j6ttZ9biXzZeeerN9zOorrWT7htXXcYph6focu2bVV0ybLz3tjA+Y2Gou+1g1x3Awjl//HLvBTI/fIPFs84hic6F9rJpjOBjHr3+O3WCmx6/veLZ2RGFmZqPR5hGFmZmNQKsShaRzJf1B0pSkSUk35/IzJT0maV++XZjLF+X9j0na2FXPAkm7u7ZDku6sq1+jVCqG+bk1kvZIekHSI5IW19GnUSocv6/k2E1K2lBHf0apj9hdI2lXPsd2Sfp8V12X5PKXJN0lad4vXl04fuslvSLp2JwOXvLjWMPegLOAT+f7C4C9pKVTNwC35PJbgNvy/VOAK4CbgI0V9e4Crqy7f22KIWmesNeAxfnxBtIqhbX3sSXxWwT8E1iSH98HXF13/xoWu5XA2fn+CuDVrrqeAy4jrWXze+DauvvXsvhdmus7Npdjt2pEEREHI+Kv+f5RYIq0Ot4XSX9o5Nsv5X2OR8STwP9mq1PSecBS4E9DbHpjFIyh8nZKfjd3GjOsITLfFIzfJ4C9EdFZgetx4LohN79WfcTu+YjonFOTwIclfUjSWcBpEfHnSK9693d+Zj4rFb/83DMRcXCux25VougmaRkpYz4LfKTT6Xzby+rna4Ct+YQbK4PEMCLeBr4B7CEliOXAvUNsbuMMeA6+BFwgaZmkCdIf99istdlH7K4Dno+It0gvjge6njuQy8bGgPHrWSsThaRTgYeA70bEGwNWtxp4cPBWtcugMZR0IilRrATOBl4Abi3ayAYbNH4RcZgUv62k0ezLwDsl29hUvcZO0oXAbcCNnaIZdhubN3oF4tez1iWK/AL1EPDziPhVLv53Ho6Sb1+bY10XARMRsWsojW2oQjG8GCAi9ufR2C+Azw6pyY1S6hyMiO0R8ZmIuIw0b9m+YbW5KXqNnaRzgF8DayNify4+QFouuWPGpZPno0Lx61mrEkW+Fn4vMBURd3Q9tQ1Yl++vA34zxyrXMGajiYIxfBVYLqkzudg1pGum81rJc1DS0ny7EPgmcE/Z1jZLr7GTdAbwO+DWiHiqs3O+vHJU0qW5zrXM/W++tUrFry91/ye/l4306ZEgXebYnbcvkD5B8gTpHdkTwJldP/My8B/gGOmdyPKu5/4BXFB3v9oaQ9IneaZyXduBRXX3r2XxexD4W95W1923psUO+AFwvGvf3cDS/Nwq4EVgP7CR/OXh+bwVjt+GfC7+P9/+qOrY/ma2mZlVatWlJzMzGz0nCjMzq+REYWZmlZwozMyskhOFmZlVcqIwGwJJN0la28P+yyS9OMw2mfVrou4GmM03kiYiYlPd7TArxYnCbAZ50rVHSJOurSRN6bwW+CRwB3AqcAj4WkQclPRH4GngcmCbpAWkKZxvl3QxsAk4mfQFsa9HxGFJlwBbgDeBJ0fXO7Pe+NKT2ezOBzZHxKeAN4BvAXcDX46Izov8+q79z4iIqyLix9PquR/4Xq5nD/DDXP5T4DuR5noyayyPKMxm90q8N0fOz4DvkxaAeSwvqHYC0D2n/9bpFUg6nZRAduSi+4BfzlD+AHBt+S6YDc6Jwmx20+e3OQpMVowAjvdQt2ao36yRfOnJbHYfldRJCmuAZ4AlnTJJJ+a5/mcVEUeAw5I+l4u+CuyIiP8CRyRdkcuvL998szI8ojCb3RSwTtJPSDNz3g08CtyVLx1NAHeSlpmssg7YJOlk0ozFN+TyG4Atkt7M9Zo1kmePNZtB/tTTbyNiRc1NMaudLz2ZmVkljyjMzKySRxRmZlbJicLMzCo5UZiZWSUnCjMzq+REYWZmlZwozMys0rvYsIa5HbYGkgAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'][-200:].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": {}, + "outputs": [], + "source": [ + "first_sept_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n", + " for y in range(1991,\n", + " sorted_data.index[-1].year)]" + ] + }, + { + "cell_type": "code", + "execution_count": 20, + "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", + " #print(abs(len(one_year)-52))\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": 21, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 21, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAD8CAYAAACyyUlaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAHgNJREFUeJzt3X+QVeWd5/H3B5uAOmLAgBEQYSLjBswOBqox4242LhPAZEuwRjM9OkplqMIxmElSW5Vo6RaW8sc4lawbytKVxETUiLJsKNmNRFvc1Li7TANGE0HC0BkJIkh3qolgtujY8N0/znPldNt03+57u++P/ryqTt3T33uew3k8Vn/7+XGeo4jAzMysFKMqfQFmZlb7nEzMzKxkTiZmZlYyJxMzMyuZk4mZmZXMycTMzErmZGJmZiVzMjEzs5I5mZiZWckaKn0B5faxj30spk+fXunLMDOrKa+88spvI2LiYMvXXTKZPn06O3furPRlmJnVFEm/KaV8Ud1ckr4habekXZLWSxor6R5Jb0t6LW1fyB1/p6RWSXslLcrF50p6PX23RpJSfIykZ1K8RdL0XJllkvalbVkplTUzs6HRbzKRNAX4O2BeRFwOnAU0pa8fiIg5aXsuHT8rfT8bWAw8JOmsdPzDwApgZtoWp/hy4GhEXAo8ANyfzjUBWAXMBxqBVZLGl1ZlMzMrt2IH4BuAsyU1AOcAh/o4dgnwdER0RsSbQCvQKOkiYFxEbItsqeLHgaW5MuvS/kZgQWq1LAKaI6IjIo4CzZxOQGZmViX6TSYR8TbwbeAAcBh4NyJeSF/fLumXkn6QazFMAd7KneJgik1J+z3j3cpERBfwLnBBH+cyM7MqUkw313iylsMMYDJwrqS/Juuy+gQwhyzJfKdQpJfTRB/xwZbJX+MKSTsl7Wxvb++jNmZmNhSK6eb6c+DNiGiPiPeBHwN/FhFHIuJkRJwCvkc2pgFZ6+HiXPmpZN1iB9N+z3i3Mqkr7Xygo49zdRMRayNiXkTMmzhx0DPbzEastmMn+NIj22g7fqLSl2I1qphkcgC4UtI5aRxjAbAnjYEUXAfsSvubgaY0Q2sG2UD79og4DByXdGU6zy3As7kyhZla1wMvpXGV54GFksanFtLCFDOzMlqzdR879new5sV9lb4Uq1H9PmcSES2SNgI/B7qAV4G1wPclzSHrdtoP3JqO3y1pA/BGOn5lRJxMp7sNeAw4G9iSNoBHgScktZK1SJrSuTok3QfsSMfdGxEdpVTYzE677O4tdHad+uDnJ1sO8GTLAcY0jGLv6msqeGVWa1Rv74CfN29e+KFFs+K0HTvB6uf28MLudzjx/inGjh7Fotkf564vfpJJ542t9OXZMJL0SkTMG2x5r81lNoJNGjeW88Y00Nl1ijENo+jsOsV5YxqcSGzA6m45FTMbmN++18lN8y/hxsZpPLX9AO0ehLdBcDdXTtuxE9y+/lUevPEK/2VmZiOKu7nKyDNazMwGx91ceEaL2UC5FW89uWUCvPzNq7l2zmTGjs7+c4wdPYolcybz8reurvCVmVUnt+KtJ7dM8IwWs2K5FW9n4pZJUpjRsukrV3HT/Etof6+z0pdkVnXcirczccskeeTm05MYVi+9vIJXYla93Iq3M3EyMbMB8XMp1hs/Z2JmZn7OxMzMKs/JxMzMSuZkYmZmJXMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMicTMzMrmZOJmZmVzMnEzMxKVlQykfQNSbsl7ZK0XtJYSRMkNUvalz7H546/U1KrpL2SFuXicyW9nr5bI0kpPkbSMyneIml6rsyy9G/sk7SsfFU3M7Ny6TeZSJoC/B0wLyIuB84CmoA7gK0RMRPYmn5G0qz0/WxgMfCQpLPS6R4GVgAz07Y4xZcDRyPiUuAB4P50rgnAKmA+0AisyictMzOrDsV2czUAZ0tqAM4BDgFLgHXp+3XA0rS/BHg6Ijoj4k2gFWiUdBEwLiK2Rba65OM9yhTOtRFYkFoti4DmiOiIiKNAM6cTkJmZVYl+k0lEvA18GzgAHAbejYgXgAsj4nA65jAwKRWZAryVO8XBFJuS9nvGu5WJiC7gXeCCPs5lZmZVpJhurvFkLYcZwGTgXEl/3VeRXmLRR3ywZfLXuELSTkk729vb+7g0MzMbCsV0c/058GZEtEfE+8CPgT8DjqSuK9JnWzr+IHBxrvxUsm6xg2m/Z7xbmdSVdj7Q0ce5uomItRExLyLmTZw4sYgqmZlZORWTTA4AV0o6J41jLAD2AJuBwuyqZcCzaX8z0JRmaM0gG2jfnrrCjku6Mp3nlh5lCue6Hngpjas8DyyUND61kBammJmZVZF+X9sbES2SNgI/B7qAV4G1wB8BGyQtJ0s4N6Tjd0vaALyRjl8ZESfT6W4DHgPOBrakDeBR4AlJrWQtkqZ0rg5J9wE70nH3RkRHSTU2M7Oy82t7zczMr+01M7PKczIxM7OSOZmYmVnJnEzMzErQduwEX3pkG23HT1T6UirKycTMrARrtu5jx/4O1ry4r9KXUlH9Tg02M7MPu+zuLXR2nfrg5ydbDvBkywHGNIxi7+prKnhlleGWiVkR3JVhPb38zau5ds5kxo7Ofo2OHT2KJXMm8/K3rq7wlVWGk4lZEdyVYT1NGjeW88Y00Nl1ijENo+jsOsV5YxqYdN7YSl9aRbiby6wP7sqwvvz2vU5umn8JNzZO46ntB2gfwS1XPwFv1oe2YydY/dweXtj9DifeP8XY0aNYNPvj3PXFT47Yv0CtPvkJeLMh5K4Ms+K4m8usH+7KMOufu7nMzMzdXGZmVnlOJmZmVjInEzMzK5mTiZmZlczJxEYsL5FiVj5OJlYWtfiL2UukmJWPnzOxssj/Yl593acqfTl98hIpZuXn50ysJD1/MRdU8y9mL5Fi9mF+zsQqqhaX4fYSKWbl128ykXSZpNdy2zFJX5d0j6S3c/Ev5MrcKalV0l5Ji3LxuZJeT9+tkaQUHyPpmRRvkTQ9V2aZpH1pW1be6lupavUXc2GJlE1fuYqb5l9C+3udlb4ks5rW75hJROwF5gBIOgt4G9gEfBl4ICK+nT9e0iygCZgNTAZelPQnEXESeBhYAfwT8BywGNgCLAeORsSlkpqA+4G/lDQBWAXMAwJ4RdLmiDhacs2tbGpx7apHbj7dml+99PIKXsnQaTt2gtvXv8qDN15R9cndat9AB+AXAL+OiN+kRkVvlgBPR0Qn8KakVqBR0n5gXERsA5D0OLCULJksAe5J5TcCD6ZWyyKgOSI6UplmsgS0foDXbUNoJPxirkW1NCnCat9Ak0kT3X+R3y7pFmAn8B9Ti2EKWcuj4GCKvZ/2e8ZJn28BRESXpHeBC/LxXsp8QNIKshYP06ZNG2CVzOqLZ6tZJRQ9AC/pI8C1wH9LoYeBT5B1gR0GvlM4tJfi0Ud8sGVOByLWRsS8iJg3ceLEM9bBbCSoxUkRVvsGMpvrGuDnEXEEICKORMTJiDgFfA9oTMcdBC7OlZsKHErxqb3Eu5WR1ACcD3T0cS4zO4NanRRhtW0gyeSvyHVxSboo9911wK60vxloSjO0ZgAzge0RcRg4LunKNB5yC/Bsrkxhptb1wEuRPQDzPLBQ0nhJ44GFKWZmffBsNRtuRY2ZSDoH+Dxway78D5LmkHU77S98FxG7JW0A3gC6gJVpJhfAbcBjwNlkA+9bUvxR4Ik0WN9BNjZDRHRIug/YkY67tzAYb2Zn5kkRNtz8BLyZmfkJeDMzqzwnEzMzK5mTiZmZlczJxMzMSuZkYmZmJXMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMicTMzMrmZOJmZmVzMnErIa0HTvBlx7ZRlsNvBrZRhYnE7Makn8Vr1k1Gehre82sAvwqXqt2bpmY1QC/iteqnZOJWQ3wq3it2rmby6xGFF7Fe2PjNJ7afoB2D8JbFfGbFs3MzG9aNDOzynMyMTOzkjmZmJlZyfpNJpIuk/Rabjsm6euSJkhqlrQvfY7PlblTUqukvZIW5eJzJb2evlsjSSk+RtIzKd4iaXquzLL0b+yTtKy81Tczs3LoN5lExN6ImBMRc4C5wP8DNgF3AFsjYiawNf2MpFlAEzAbWAw8JOmsdLqHgRXAzLQtTvHlwNGIuBR4ALg/nWsCsAqYDzQCq/JJy8zMqsNAu7kWAL+OiN8AS4B1Kb4OWJr2lwBPR0RnRLwJtAKNki4CxkXEtsimkD3eo0zhXBuBBanVsghojoiOiDgKNHM6AZnZCOZ1yqrLQJNJE7A+7V8YEYcB0uekFJ8CvJUrczDFpqT9nvFuZSKiC3gXuKCPc5nZCOd1yqpL0Q8tSvoIcC1wZ3+H9hKLPuKDLZO/thVk3WdMmzatn8szs1rmdcqq00BaJtcAP4+II+nnI6nrivTZluIHgYtz5aYCh1J8ai/xbmUkNQDnAx19nKubiFgbEfMiYt7EiRMHUCUzqzVep6w6DSSZ/BWnu7gANgOF2VXLgGdz8aY0Q2sG2UD79tQVdlzSlWk85JYeZQrnuh54KY2rPA8slDQ+DbwvTDEzG6G8Tll1KqqbS9I5wOeBW3Phvwc2SFoOHABuAIiI3ZI2AG8AXcDKiDiZytwGPAacDWxJG8CjwBOSWslaJE3pXB2S7gN2pOPujYiOQdTTzOqI1ymrPl6by8zMvDaXmVkxPJV4aDmZmNmI4KnEQ8vvMxlB2o6d4Pb1r/LgjVd4sNJGDE8lHh5umYwg/svMRiJPJR4ebpmMAP7LzEYyTyUeHm6ZjACD+cvMg5VWTwpTiTd95Spumn8J7e91VvqS6o5bJiPAYP4yy3eJrb7uU8N4tWbl98jNp2e8rl56eQWvpH45mYwQxT7k5S4xMxsMP7Ro3bQdO8Hq5/bwwu53OPH+KcaOHsWi2R/nri9+0n3MZnXMDy1aWXmw0swGw91c9iFe98jMBsrdXGZm5m4uMzOrPCcTMzMrmZOJmZmVzMlkGPhpcjOrd04mw8ALLJpZvfPU4CHkp8nNbKRwy2QIeelrMxspnEyGkJ8mN7ORwslkiA106WsP1ptZLSoqmUj6qKSNkn4laY+kz0i6R9Lbkl5L2xdyx98pqVXSXkmLcvG5kl5P362RpBQfI+mZFG+RND1XZpmkfWlbVr6qD49Hbp7H6qWXM2vyOFYvvbzbUti98WC9mdWiopZTkbQOeDkivi/pI8A5wNeB9yLi2z2OnQWsBxqBycCLwJ9ExElJ24GvAf8EPAesiYgtkr4C/OuI+FtJTcB1EfGXkiYAO4F5QACvAHMj4uiZrrVWl1PpOVhf4MF6MxsOQ76ciqRxwGeBRwEi4g8R8bs+iiwBno6Izoh4E2gFGiVdBIyLiG2RZbDHgaW5MuvS/kZgQWq1LAKaI6IjJZBmYPGAa1kDPFhvZrWsmG6uPwbagR9KelXS9yWdm767XdIvJf1A0vgUmwK8lSt/MMWmpP2e8W5lIqILeBe4oI9z1R0P1ptZLSsmmTQAnwYejogrgN8DdwAPA58A5gCHge+k49XLOaKP+GDLfEDSCkk7Je1sb2/voyrVze+pNrNaVcxDiweBgxHRkn7eCNwREUcKB0j6HvA/c8dfnCs/FTiU4lN7iefLHJTUAJwPdKT453qU+VnPC4yItcBayMZMiqhTydqOneD29a/y4I1XlK314PdUm1mt6rdlEhHvAG9JuiyFFgBvpDGQguuAXWl/M9CUZmjNAGYC2yPiMHBc0pVpPOQW4NlcmcJMreuBl9K4yvPAQknjUzfawhSrOM+6MiuOp7uPDMUup/JV4EdpJte/AF8G1kiaQ9bttB+4FSAidkvaALwBdAErI+JkOs9twGPA2cCWtEE2uP+EpFayFklTOleHpPuAHem4eyOiY3BVLQ8vkWI2MPk/vFZf96lKX44NEb9pcYDajp1g9XN7eGH3O5x4/xRjR49i0eyPc9cXP+nBcrMcT3evLX7T4jDzrCuz4ni6+8jiVYMHoTDr6sbGaTy1/QDt7gs2+xD/4TWyOJkMgmddWbkMxazAauI/vEYOj5mYVdDdm17nR9sPcFPjNA9OW0WVOmbilolZBXhWoNUbD8CblVkxz1V4cNrqjZOJWZkV80CrB6et3riby6xMBtp15cHpD6v3CQn1zAPwZmXiB1pL5wkJleMBeLMq4a6rwfOEhNrnMROzMvJrBAbHExJqn1smVtVqrQ/dD7QOjlt1tc8tExt2A1mS3Ev9jxxu1dU2D8DbsCtmkNUrzpoNr1IH4J1MbNgMJEF4ZpTZ8PIS9FYzBjLI6j50s9riAXgbNgNNEH6oz6x2OJnYsBpIgvDMKLPa4TETMzPzmImZmVWek4mZWQ8DeRbKMk4mZmY9+GHZgSsqmUj6qKSNkn4laY+kz0iaIKlZ0r70OT53/J2SWiXtlbQoF58r6fX03RpJSvExkp5J8RZJ03NllqV/Y5+kZeWruplZd5fdvYXpd/yEJ1sOEJEtODn9jp9w2d1bKn1pVa/Ylsl3gZ9GxL8C/hTYA9wBbI2ImcDW9DOSZgFNwGxgMfCQpLPSeR4GVgAz07Y4xZcDRyPiUuAB4P50rgnAKmA+0AisyictM7Ny8oKTg9dvMpE0Dvgs8ChARPwhIn4HLAHWpcPWAUvT/hLg6YjojIg3gVagUdJFwLiI2BbZFLLHe5QpnGsjsCC1WhYBzRHRERFHgWZOJyAzs7Lyw7KDV0zL5I+BduCHkl6V9H1J5wIXRsRhgPQ5KR0/BXgrV/5gik1J+z3j3cpERBfwLnBBH+fqRtIKSTsl7Wxvby+iSmZmvfOCk4NTzEOLDcCnga9GRIuk75K6tM5AvcSij/hgy5wORKwF1kL2nEkf12Zm1ic/LDs4xbRMDgIHI6Il/byRLLkcSV1XpM+23PEX58pPBQ6l+NRe4t3KSGoAzgc6+jiXmZlVkX6TSUS8A7wl6bIUWgC8AWwGCrOrlgHPpv3NQFOaoTWDbKB9e+oKOy7pyjQeckuPMoVzXQ+8lMZVngcWShqfBt4XppiZmVWRYtfm+irwI0kfAf4F+DJZItogaTlwALgBICJ2S9pAlnC6gJURcTKd5zbgMeBsYEvaIBvcf0JSK1mLpCmdq0PSfcCOdNy9EdExyLqamdkQ8dpcZmbmtbnMzKzynEzMzOpApdcTczIxM6sDlV5PzC/HMjOrYZfdvYXOrlMf/PxkywGebDnAmIZR7F19zbBdh1smZmY1rFrWE3MyMTOrYdWynpi7uczMalxhPbEbG6fx1PYDtFdgEN7PmZiZmZ8zMTOzynMyMTOzkjmZmJlZyZxMzMysZE4mZmZWMieTGlfp9XjMzMDJpOZVej0eMzPwQ4s1q1rW4zEzA7dMala1rMdjZgZOJjWrWtbjMTMDd3PVtGpYj8fMDLw2l5mZ4bW5zMysChSVTCTtl/S6pNck7UyxeyS9nWKvSfpC7vg7JbVK2itpUS4+N52nVdIaSUrxMZKeSfEWSdNzZZZJ2pe2ZeWquJnZcKvn58IG0jK5OiLm9GgGPZBicyLiOQBJs4AmYDawGHhI0lnp+IeBFcDMtC1O8eXA0Yi4FHgAuD+dawKwCpgPNAKrJI0fRD3NzCqunp8LG4oB+CXA0xHRCbwpqRVolLQfGBcR2wAkPQ4sBbakMvek8huBB1OrZRHQHBEdqUwzWQJaPwTXbWY2JEbCc2HFtkwCeEHSK5JW5OK3S/qlpB/kWgxTgLdyxxxMsSlpv2e8W5mI6ALeBS7o41xmZjVjJDwXVmwyuSoiPg1cA6yU9FmyLqtPAHOAw8B30rHqpXz0ER9smQ9IWiFpp6Sd7e3tfVbEzGy4jYTnwopKJhFxKH22AZuAxog4EhEnI+IU8D2yMQ3IWg8X54pPBQ6l+NRe4t3KSGoAzgc6+jhXz+tbGxHzImLexIkTi6mSmdmwKjwXtukrV3HT/Etof6+z3zK1NGDfbzKRdK6k8wr7wEJgl6SLcoddB+xK+5uBpjRDawbZQPv2iDgMHJd0ZRoPuQV4NlemMFPreuClyB6AeR5YKGl86kZbmGJmZjXlkZvnsXrp5cyaPI7VSy/nkZv7f6SjlgbsixmAvxDYlGbxNgBPRcRPJT0haQ5Zt9N+4FaAiNgtaQPwBtAFrIyIk+lctwGPAWeTDbxvSfFHgSfSYH0H2WwwIqJD0n3AjnTcvYXBeDOzelWLA/Z+At7MrMq0HTvB6uf28MLudzjx/inGjh7Fotkf564vfnLIxln8BLyZWZ2pxQF7L/RoZlaFam0hV3dzmZmZu7nMzKzynEzMzKxkTiZmZlYyJxMzMyuZk4mZmZXMycTMzErmZGJmZiVzMjEzs5I5mZiZWcmcTKxu1NK7H8zqjZOJ1Y1aeveDWb3xQo9W82rx3Q9m9cYtE6t5L3/zaq6dM5mxo7P/nceOHsWSOZN5+VtXV/jKzEYOJxOrebX47gezeuNuLqsLtfbuB7N64/eZmJmZ32diZmaV52RiZmYlKyqZSNov6XVJr0namWITJDVL2pc+x+eOv1NSq6S9khbl4nPTeVolrZGkFB8j6ZkUb5E0PVdmWfo39klaVq6Km5lZ+QykZXJ1RMzJ9andAWyNiJnA1vQzkmYBTcBsYDHwkKSzUpmHgRXAzLQtTvHlwNGIuBR4ALg/nWsCsAqYDzQCq/JJy8zMqkMp3VxLgHVpfx2wNBd/OiI6I+JNoBVolHQRMC4itkU26v94jzKFc20EFqRWyyKgOSI6IuIo0MzpBGRmZlWi2GQSwAuSXpG0IsUujIjDAOlzUopPAd7KlT2YYlPSfs94tzIR0QW8C1zQx7nMzKyKFPucyVURcUjSJKBZ0q/6OFa9xKKP+GDLnP4HswRXSHLvSdrbx/XVgo8Bv630RQyxeq9jvdcP6r+OI61+l5RysqKSSUQcSp9tkjaRjV8ckXRRRBxOXVht6fCDwMW54lOBQyk+tZd4vsxBSQ3A+UBHin+uR5mf9XJ9a4G1xdSlFkjaWcp871pQ73Ws9/pB/dfR9RuYfru5JJ0r6bzCPrAQ2AVsBgqzq5YBz6b9zUBTmqE1g2ygfXvqCjsu6co0HnJLjzKFc10PvJTGVZ4HFkoanwbeF6aYmZlVkWJaJhcCm9Is3gbgqYj4qaQdwAZJy4EDwA0AEbFb0gbgDaALWBkRJ9O5bgMeA84GtqQN4FHgCUmtZC2SpnSuDkn3ATvScfdGREcJ9TUzsyFQd8up1ANJK1LXXd2q9zrWe/2g/uvo+g3wfE4mZmZWKi+nYmZmJXMyGSaSfiCpTdKuXOxPJW1LS8z8D0njUvwjkn6Y4r+Q9LlcmZ+lZWpeS9ukXv65YSfpYkn/S9IeSbslfS3Fy7bsTiWVuX51cQ8lXZCOf0/Sgz3OVfP3sJ/6Vd09HET9Pq/s2cHX0+e/z51r4PcvIrwNwwZ8Fvg0sCsX2wH8u7T/N8B9aX8l8MO0Pwl4BRiVfv4ZMK/S9emlfhcBn0775wH/DMwC/gG4I8XvAO5P+7OAXwBjgBnAr4Gz0nfbgc+QPWe0BbimzupXL/fwXODfAH8LPNjjXPVwD/uqX9Xdw0HU7wpgctq/HHi7lPvnlskwiYh/JJuplncZ8I9pvxn4i7Q/i2y9MyKiDfgdUNXz3SPicET8PO0fB/aQrVZQzmV3KqZc9Rveqx6YgdYxIn4fEf8b6PYmsnq5h2eqX7UaRP1ejfQMIbAbGKvskY5B3T8nk8raBVyb9m/g9MOevwCWSGpQ9qzOXLo/CPrD1LT+T9XQfdCTslWfrwBaKO+yO1WhxPoV1MM9PJN6uYf9qdp7OIj6/QXwakR0Msj752RSWX8DrJT0Clmz9A8p/gOyG7gT+C/A/yV7Zgfgpoj4FPBv03bzsF5xPyT9EfDfga9HxLG+Du0lVvQSOpVShvpB/dzDM56il1gt3sO+VO09HGj9JM0mW6n91kKol8P6vX9OJhUUEb+KiIURMRdYT9avTkR0RcQ3IlvyfwnwUWBf+u7t9HkceIoq6jqRNJrsf+IfRcSPU/hIajYXuj9KWXanospUv3q6h2dSL/fwjKr1Hg60fpKmApuAWyLi1yk8qPvnZFJBhRkgkkYBdwP/Nf18jrKla5D0eaArIt5I3V4fS/HRwH8g6yqruNTMfxTYExH/OfdVOZfdqZhy1a/O7mGv6ugenuk8VXkPB1o/SR8FfgLcGRH/p3DwoO9fpWYejLSNrOVxGHifLPMvB75GNuPin4G/5/RDpNOBvWQDaC8Cl6T4uWQzu35JNmD2XdIMoUpvZLNeIl3ba2n7AtmrBLaStay2AhNyZe4ia43tJTdbhGyywa703YOF/y71UL86vIf7ySaWvJf+v55VZ/fwQ/Wr1ns40PqR/QH7+9yxrwGTBnv//AS8mZmVzN1cZmZWMicTMzMrmZOJmZmVzMnEzMxK5mRiZmYlczIxM7OSOZmYmVnJnEzMzKxk/x+qftJYVEr3+gAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.plot(style='*')" + ] + }, + { + "cell_type": "code", + "execution_count": 22, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "2002 516689\n", + "2018 542312\n", + "2017 551041\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": 22, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "yearly_incidence.sort_values()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], "metadata": { "kernelspec": { "display_name": "Python 3", @@ -16,10 +2275,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.3" + "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 } -