diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb index 0bbbe371b01e359e381e43239412d77bf53fb1fb..46fa8fb25d7f7b814d1692ad11c584a63497c67d 100644 --- a/module3/exo2/exercice.ipynb +++ b/module3/exo2/exercice.ipynb @@ -1,5 +1,2270 @@ { - "cells": [], + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chicken pox analysis" + ] + }, + { + "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": 3, + "metadata": {}, + "outputs": [], + "source": [ + "data_url = \"https://www.sentiweb.fr/datasets/incidence-PAY-7.csv\"" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "The first line of the CSV file is a comment, which we ignore with skip=1." + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020201971590490001FRFrance
12020187824851563102FRFrance
220201772720658001FRFrance
32020167758781438102FRFrance
4202015719186753161315FRFrance
52020147387922275531639FRFrance
6202013773265236941611814FRFrance
72020127812357901045612816FRFrance
8202011710198756812828151119FRFrance
920201079011669111331141018FRFrance
102020097136311054416718211626FRFrance
11202008710424770813140161220FRFrance
1220200778959657411344141018FRFrance
1320200679264692511603141018FRFrance
1420200578505631410696131016FRFrance
152020047799158311015112915FRFrance
1620200375968410078369612FRFrance
17202002765344530853810713FRFrance
1820200179835701912651151119FRFrance
192019527794152461063612816FRFrance
2020195175823367579719612FRFrance
21201950764244276857210713FRFrance
22201949766214540870210713FRFrance
2320194875542338377018511FRFrance
242019477753650581001411715FRFrance
252019467263813163960426FRFrance
2620194574492261563697410FRFrance
2720194475728362778299612FRFrance
2820194374834275169177410FRFrance
29201942762793989856910713FRFrance
.................................
15061991267176081130423912312042FRFrance
15071991257161691070021638281838FRFrance
15081991247161711007122271281739FRFrance
1509199123711947767116223211329FRFrance
1510199122715452995320951271737FRFrance
1511199121714903897520831261636FRFrance
15121991207190531274225364342345FRFrance
15131991197167391124622232291939FRFrance
15141991187213851388228888382551FRFrance
1515199117713462887718047241632FRFrance
15161991167148571006819646261834FRFrance
1517199115713975978118169251832FRFrance
1518199114712265768416846221430FRFrance
151919911379567604113093171123FRFrance
1520199112710864733114397191325FRFrance
15211991117155741118419964271935FRFrance
15221991107166431137221914292038FRFrance
1523199109713741878018702241533FRFrance
1524199108713289881317765231531FRFrance
1525199107712337807716597221529FRFrance
1526199106710877701314741191226FRFrance
1527199105710442654414340181125FRFrance
15281991047791345631126314820FRFrance
15291991037153871048420290271836FRFrance
15301991027162771104621508292038FRFrance
15311991017155651027120859271836FRFrance
15321990527193751329525455342345FRFrance
15331990517190801380724353342543FRFrance
1534199050711079666015498201228FRFrance
15351990497114302610205FRFrance
\n", + "

1536 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202019 7 159 0 490 0 0 \n", + "1 202018 7 824 85 1563 1 0 \n", + "2 202017 7 272 0 658 0 0 \n", + "3 202016 7 758 78 1438 1 0 \n", + "4 202015 7 1918 675 3161 3 1 \n", + "5 202014 7 3879 2227 5531 6 3 \n", + "6 202013 7 7326 5236 9416 11 8 \n", + "7 202012 7 8123 5790 10456 12 8 \n", + "8 202011 7 10198 7568 12828 15 11 \n", + "9 202010 7 9011 6691 11331 14 10 \n", + "10 202009 7 13631 10544 16718 21 16 \n", + "11 202008 7 10424 7708 13140 16 12 \n", + "12 202007 7 8959 6574 11344 14 10 \n", + "13 202006 7 9264 6925 11603 14 10 \n", + "14 202005 7 8505 6314 10696 13 10 \n", + "15 202004 7 7991 5831 10151 12 9 \n", + "16 202003 7 5968 4100 7836 9 6 \n", + "17 202002 7 6534 4530 8538 10 7 \n", + "18 202001 7 9835 7019 12651 15 11 \n", + "19 201952 7 7941 5246 10636 12 8 \n", + "20 201951 7 5823 3675 7971 9 6 \n", + "21 201950 7 6424 4276 8572 10 7 \n", + "22 201949 7 6621 4540 8702 10 7 \n", + "23 201948 7 5542 3383 7701 8 5 \n", + "24 201947 7 7536 5058 10014 11 7 \n", + "25 201946 7 2638 1316 3960 4 2 \n", + "26 201945 7 4492 2615 6369 7 4 \n", + "27 201944 7 5728 3627 7829 9 6 \n", + "28 201943 7 4834 2751 6917 7 4 \n", + "29 201942 7 6279 3989 8569 10 7 \n", + "... ... ... ... ... ... ... ... \n", + "1506 199126 7 17608 11304 23912 31 20 \n", + "1507 199125 7 16169 10700 21638 28 18 \n", + "1508 199124 7 16171 10071 22271 28 17 \n", + "1509 199123 7 11947 7671 16223 21 13 \n", + "1510 199122 7 15452 9953 20951 27 17 \n", + "1511 199121 7 14903 8975 20831 26 16 \n", + "1512 199120 7 19053 12742 25364 34 23 \n", + "1513 199119 7 16739 11246 22232 29 19 \n", + "1514 199118 7 21385 13882 28888 38 25 \n", + "1515 199117 7 13462 8877 18047 24 16 \n", + "1516 199116 7 14857 10068 19646 26 18 \n", + "1517 199115 7 13975 9781 18169 25 18 \n", + "1518 199114 7 12265 7684 16846 22 14 \n", + "1519 199113 7 9567 6041 13093 17 11 \n", + "1520 199112 7 10864 7331 14397 19 13 \n", + "1521 199111 7 15574 11184 19964 27 19 \n", + "1522 199110 7 16643 11372 21914 29 20 \n", + "1523 199109 7 13741 8780 18702 24 15 \n", + "1524 199108 7 13289 8813 17765 23 15 \n", + "1525 199107 7 12337 8077 16597 22 15 \n", + "1526 199106 7 10877 7013 14741 19 12 \n", + "1527 199105 7 10442 6544 14340 18 11 \n", + "1528 199104 7 7913 4563 11263 14 8 \n", + "1529 199103 7 15387 10484 20290 27 18 \n", + "1530 199102 7 16277 11046 21508 29 20 \n", + "1531 199101 7 15565 10271 20859 27 18 \n", + "1532 199052 7 19375 13295 25455 34 23 \n", + "1533 199051 7 19080 13807 24353 34 25 \n", + "1534 199050 7 11079 6660 15498 20 12 \n", + "1535 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 1 FR France \n", + "1 2 FR France \n", + "2 1 FR France \n", + "3 2 FR France \n", + "4 5 FR France \n", + "5 9 FR France \n", + "6 14 FR France \n", + "7 16 FR France \n", + "8 19 FR France \n", + "9 18 FR France \n", + "10 26 FR France \n", + "11 20 FR France \n", + "12 18 FR France \n", + "13 18 FR France \n", + "14 16 FR France \n", + "15 15 FR France \n", + "16 12 FR France \n", + "17 13 FR France \n", + "18 19 FR France \n", + "19 16 FR France \n", + "20 12 FR France \n", + "21 13 FR France \n", + "22 13 FR France \n", + "23 11 FR France \n", + "24 15 FR France \n", + "25 6 FR France \n", + "26 10 FR France \n", + "27 12 FR France \n", + "28 10 FR France \n", + "29 13 FR France \n", + "... ... ... ... \n", + "1506 42 FR France \n", + "1507 38 FR France \n", + "1508 39 FR France \n", + "1509 29 FR France \n", + "1510 37 FR France \n", + "1511 36 FR France \n", + "1512 45 FR France \n", + "1513 39 FR France \n", + "1514 51 FR France \n", + "1515 32 FR France \n", + "1516 34 FR France \n", + "1517 32 FR France \n", + "1518 30 FR France \n", + "1519 23 FR France \n", + "1520 25 FR France \n", + "1521 35 FR France \n", + "1522 38 FR France \n", + "1523 33 FR France \n", + "1524 31 FR France \n", + "1525 29 FR France \n", + "1526 26 FR France \n", + "1527 25 FR France \n", + "1528 20 FR France \n", + "1529 36 FR France \n", + "1530 38 FR France \n", + "1531 36 FR France \n", + "1532 45 FR France \n", + "1533 43 FR France \n", + "1534 28 FR France \n", + "1535 5 FR France \n", + "\n", + "[1536 rows x 10 columns]" + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "raw_data = pd.read_csv(data_url, skiprows=1)\n", + "raw_data" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
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": 5, + "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
020201971590490001FRFrance
12020187824851563102FRFrance
220201772720658001FRFrance
32020167758781438102FRFrance
4202015719186753161315FRFrance
52020147387922275531639FRFrance
6202013773265236941611814FRFrance
72020127812357901045612816FRFrance
8202011710198756812828151119FRFrance
920201079011669111331141018FRFrance
102020097136311054416718211626FRFrance
11202008710424770813140161220FRFrance
1220200778959657411344141018FRFrance
1320200679264692511603141018FRFrance
1420200578505631410696131016FRFrance
152020047799158311015112915FRFrance
1620200375968410078369612FRFrance
17202002765344530853810713FRFrance
1820200179835701912651151119FRFrance
192019527794152461063612816FRFrance
2020195175823367579719612FRFrance
21201950764244276857210713FRFrance
22201949766214540870210713FRFrance
2320194875542338377018511FRFrance
242019477753650581001411715FRFrance
252019467263813163960426FRFrance
2620194574492261563697410FRFrance
2720194475728362778299612FRFrance
2820194374834275169177410FRFrance
29201942762793989856910713FRFrance
.................................
15061991267176081130423912312042FRFrance
15071991257161691070021638281838FRFrance
15081991247161711007122271281739FRFrance
1509199123711947767116223211329FRFrance
1510199122715452995320951271737FRFrance
1511199121714903897520831261636FRFrance
15121991207190531274225364342345FRFrance
15131991197167391124622232291939FRFrance
15141991187213851388228888382551FRFrance
1515199117713462887718047241632FRFrance
15161991167148571006819646261834FRFrance
1517199115713975978118169251832FRFrance
1518199114712265768416846221430FRFrance
151919911379567604113093171123FRFrance
1520199112710864733114397191325FRFrance
15211991117155741118419964271935FRFrance
15221991107166431137221914292038FRFrance
1523199109713741878018702241533FRFrance
1524199108713289881317765231531FRFrance
1525199107712337807716597221529FRFrance
1526199106710877701314741191226FRFrance
1527199105710442654414340181125FRFrance
15281991047791345631126314820FRFrance
15291991037153871048420290271836FRFrance
15301991027162771104621508292038FRFrance
15311991017155651027120859271836FRFrance
15321990527193751329525455342345FRFrance
15331990517190801380724353342543FRFrance
1534199050711079666015498201228FRFrance
15351990497114302610205FRFrance
\n", + "

1536 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202019 7 159 0 490 0 0 \n", + "1 202018 7 824 85 1563 1 0 \n", + "2 202017 7 272 0 658 0 0 \n", + "3 202016 7 758 78 1438 1 0 \n", + "4 202015 7 1918 675 3161 3 1 \n", + "5 202014 7 3879 2227 5531 6 3 \n", + "6 202013 7 7326 5236 9416 11 8 \n", + "7 202012 7 8123 5790 10456 12 8 \n", + "8 202011 7 10198 7568 12828 15 11 \n", + "9 202010 7 9011 6691 11331 14 10 \n", + "10 202009 7 13631 10544 16718 21 16 \n", + "11 202008 7 10424 7708 13140 16 12 \n", + "12 202007 7 8959 6574 11344 14 10 \n", + "13 202006 7 9264 6925 11603 14 10 \n", + "14 202005 7 8505 6314 10696 13 10 \n", + "15 202004 7 7991 5831 10151 12 9 \n", + "16 202003 7 5968 4100 7836 9 6 \n", + "17 202002 7 6534 4530 8538 10 7 \n", + "18 202001 7 9835 7019 12651 15 11 \n", + "19 201952 7 7941 5246 10636 12 8 \n", + "20 201951 7 5823 3675 7971 9 6 \n", + "21 201950 7 6424 4276 8572 10 7 \n", + "22 201949 7 6621 4540 8702 10 7 \n", + "23 201948 7 5542 3383 7701 8 5 \n", + "24 201947 7 7536 5058 10014 11 7 \n", + "25 201946 7 2638 1316 3960 4 2 \n", + "26 201945 7 4492 2615 6369 7 4 \n", + "27 201944 7 5728 3627 7829 9 6 \n", + "28 201943 7 4834 2751 6917 7 4 \n", + "29 201942 7 6279 3989 8569 10 7 \n", + "... ... ... ... ... ... ... ... \n", + "1506 199126 7 17608 11304 23912 31 20 \n", + "1507 199125 7 16169 10700 21638 28 18 \n", + "1508 199124 7 16171 10071 22271 28 17 \n", + "1509 199123 7 11947 7671 16223 21 13 \n", + "1510 199122 7 15452 9953 20951 27 17 \n", + "1511 199121 7 14903 8975 20831 26 16 \n", + "1512 199120 7 19053 12742 25364 34 23 \n", + "1513 199119 7 16739 11246 22232 29 19 \n", + "1514 199118 7 21385 13882 28888 38 25 \n", + "1515 199117 7 13462 8877 18047 24 16 \n", + "1516 199116 7 14857 10068 19646 26 18 \n", + "1517 199115 7 13975 9781 18169 25 18 \n", + "1518 199114 7 12265 7684 16846 22 14 \n", + "1519 199113 7 9567 6041 13093 17 11 \n", + "1520 199112 7 10864 7331 14397 19 13 \n", + "1521 199111 7 15574 11184 19964 27 19 \n", + "1522 199110 7 16643 11372 21914 29 20 \n", + "1523 199109 7 13741 8780 18702 24 15 \n", + "1524 199108 7 13289 8813 17765 23 15 \n", + "1525 199107 7 12337 8077 16597 22 15 \n", + "1526 199106 7 10877 7013 14741 19 12 \n", + "1527 199105 7 10442 6544 14340 18 11 \n", + "1528 199104 7 7913 4563 11263 14 8 \n", + "1529 199103 7 15387 10484 20290 27 18 \n", + "1530 199102 7 16277 11046 21508 29 20 \n", + "1531 199101 7 15565 10271 20859 27 18 \n", + "1532 199052 7 19375 13295 25455 34 23 \n", + "1533 199051 7 19080 13807 24353 34 25 \n", + "1534 199050 7 11079 6660 15498 20 12 \n", + "1535 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 1 FR France \n", + "1 2 FR France \n", + "2 1 FR France \n", + "3 2 FR France \n", + "4 5 FR France \n", + "5 9 FR France \n", + "6 14 FR France \n", + "7 16 FR France \n", + "8 19 FR France \n", + "9 18 FR France \n", + "10 26 FR France \n", + "11 20 FR France \n", + "12 18 FR France \n", + "13 18 FR France \n", + "14 16 FR France \n", + "15 15 FR France \n", + "16 12 FR France \n", + "17 13 FR France \n", + "18 19 FR France \n", + "19 16 FR France \n", + "20 12 FR France \n", + "21 13 FR France \n", + "22 13 FR France \n", + "23 11 FR France \n", + "24 15 FR France \n", + "25 6 FR France \n", + "26 10 FR France \n", + "27 12 FR France \n", + "28 10 FR France \n", + "29 13 FR France \n", + "... ... ... ... \n", + "1506 42 FR France \n", + "1507 38 FR France \n", + "1508 39 FR France \n", + "1509 29 FR France \n", + "1510 37 FR France \n", + "1511 36 FR France \n", + "1512 45 FR France \n", + "1513 39 FR France \n", + "1514 51 FR France \n", + "1515 32 FR France \n", + "1516 34 FR France \n", + "1517 32 FR France \n", + "1518 30 FR France \n", + "1519 23 FR France \n", + "1520 25 FR France \n", + "1521 35 FR France \n", + "1522 38 FR France \n", + "1523 33 FR France \n", + "1524 31 FR France \n", + "1525 29 FR France \n", + "1526 26 FR France \n", + "1527 25 FR France \n", + "1528 20 FR France \n", + "1529 36 FR France \n", + "1530 38 FR France \n", + "1531 36 FR France \n", + "1532 45 FR France \n", + "1533 43 FR France \n", + "1534 28 FR France \n", + "1535 5 FR France \n", + "\n", + "[1536 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+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm4HUWZ/vudc7ck5GbfE0ggYUkCBIjsKrJGVEABJ+oAjvjDYcBlGHFAR2VUFEaRGVxg2GQZlU0UEAHZ15BwCdlDyCULuUnIvueu59Tvj+46XV1dW5/lLrn1Ps99bp/q6qrq6qr66luLGGPw8PDw8PAQkenqBnh4eHh4dD944uDh4eHhkYAnDh4eHh4eCXji4OHh4eGRgCcOHh4eHh4JeOLg4eHh4ZGAJw4eHh4eHgl44uDh4eHhkYAnDh4eHh4eCVR1dQOKxdChQ9n48eO7uhkeHh4ePQpvv/32ZsbYMFu+Hkscxo8fj4aGhq5uhoeHh0ePAhGtdsnnxUoeHh4eHgl44uDh4eHhkYAnDh4eHh4eCXji4OHh4eGRgCcOHh4eHh4JeOLg4eHh4ZGAJw4eHh4eHgl44uDRK8AYw6Nzm9Dcluvqpnh49Ah44uDRKzBrxRZc9dB8/PjJJV3dFA+PHgFPHDx6BXa3dAAANu5s6eKWeHj0DHji4NErQERd3QQPjx4FTxw8PDw8PBKwEgciqiOiOUQ0n4gWE9F/hunXEdFaIpoX/p0tPHMtETUS0TIiOktIP4aIFob3bqFwO0dEtUT0YJg+m4jGl/9VPTw8PDxc4cI5tAI4lTF2JIBpAGYQ0fHhvZsZY9PCv78BABFNBjATwBQAMwD8loiyYf5bAVwGYFL4NyNMvxTANsbYRAA3A7ix9Ffz8PDw8CgWVuLAAuwOf1aHf8zwyLkAHmCMtTLGVgJoBHAsEY0CUM8Ym8UYYwDuA3Ce8My94fUjAE4jLyT28PDw6DI46RyIKEtE8wBsBPAsY2x2eOtKIlpARHcT0aAwbQyANcLjTWHamPBaTo89wxjrALADwBBFOy4jogYiati0aZPTC3p4eHh4pIcTcWCM5Rhj0wCMRcAFTEUgIjoIgahpPYCbwuyqHT8zpJuekdtxO2NsOmNs+rBh1oOMPDwSYCaeV4MVm3ZjQdP28jfGw6MbI5W1EmNsO4CXAMxgjG0IiUYewB0Ajg2zNQEYJzw2FsC6MH2sIj32DBFVARgAYGuqN/HwMKAUGeWpN72Mc379etna4uHRE+BirTSMiAaG130AnA7g3VCHwPFZAIvC68cBzAwtkCYgUDzPYYytB7CLiI4P9QkXA3hMeOaS8PoCAC+EegkPDw8Pjy6AyxnSowDcG1ocZQA8xBj7KxHdT0TTEIh/VgH4GgAwxhYT0UMAlgDoAHAFY4wHtLkcwD0A+gB4KvwDgLsA3E9EjQg4hplleDcPDw8PjyJhJQ6MsQUAjlKkX2R45noA1yvSGwBMVaS3ALjQ1hYPDw8Pj86B95D26JW46K7ZeHLB+q5uhodHt4UnDh69Eq8u34wr/jC3q5vh4dFt4YmDR6+Ct3Lw8HCDJw4evQLe397DIx08cfDw6EZoz+Vx1YPzsGLTbntmD48KwhMHD49uhPlrtuPRd9bi6kcWdHVTPHo5PHHw8PDw8EjAEwcPDw8PjwQ8cfDoVejuUVm6d+s8ehM8cfDoFfDWSh4e6eCJg4dHN4KnYR7dBZ449DA0rNqKvW0dXd0MDw+PfRyeOPQgbNzVggtum4WrHpzf1U3xqBC8zsGju8AThx6E5rYg8vmS9Tu7uCU9G91dKe3h0R3giYOHRzeC1zl4dBd44uDRq8BQ3DnSnYVu3LR9Gs1tOdz/5mrsbvX6PA5PHDx6BaiH7cl7Vmt7Ph6btxbf/8si3PvGqq5uSreBJw4evQ49YXfeE9q4L2FHc3vsv4cDcSCiOiKaQ0TziWgxEf1nmD6YiJ4louXh/0HCM9cSUSMRLSOis4T0Y4hoYXjvFqLANYmIaonowTB9NhGNL/+renh0f3Qlx3D1w/Mx/ponu7AFXQ9vrBDBhXNoBXAqY+xIANMAzCCi4wFcA+B5xtgkAM+Hv0FEkwHMBDAFwAwAvyWibFjWrQAuAzAp/JsRpl8KYBtjbCKAmwHcWIZ3KzueXLAea7bu7epmeJSI7rwAdGXLHn67qQtr71pwD/puPDQ6HVbiwALw4PLV4R8DcC6Ae8P0ewGcF16fC+ABxlgrY2wlgEYAxxLRKAD1jLFZLJid90nP8LIeAXAa5yq6E674w1yc8+vXuqx+P3B7D7rd4N/HwXVSfopFcNI5EFGWiOYB2AjgWcbYbAAjGGPrASD8PzzMPgbAGuHxpjBtTHgtp8eeYYx1ANgBYEgxL1RpbNvrZZI9GYz1jAWgJ7RxX4LnHJJwIg6MsRxjbBqAsQi4gKmG7KpNDzOkm56JF0x0GRE1EFHDpk2bbM328CiA+eXWwwF+nERIZa3EGNsO4CUEuoINoagI4f+NYbYmAOOEx8YCWBemj1Wkx54hoioAAwBsVdR/O2NsOmNs+rBhw9I0fZ9A9xO0VQZrtu4tu15ALK4n7A4r8alnvb/F68w04FLsnjA2Ogsu1krDiGhgeN0HwOkA3gXwOIBLwmyXAHgsvH4cwMzQAmkCAsXznFD0tIuIjg/1CRdLz/CyLgDwAuvGWsPv/2URtu9t6+pm7JNYvG4HPvpfL+Lu11c5P/Pa8s24+O45yOf1Q6b7jqbOwxfueBMf/a8Xu7oZ3RK9ZN+VCi6cwygALxLRAgBvIdA5/BXADQDOIKLlAM4If4MxthjAQwCWAHgawBWMsVxY1uUA7kSgpH4fwFNh+l0AhhBRI4CrEFo+dVfc/+Zq/PyZZWUr770Nu7CzxesyAOCDLcHOds7KLc7PXP5/b+OV9zZhtyFabU+jDeVq786WdrzzwbYylbbvItI59LSRUjlU2TIwxhYAOEqRvgXAaZpnrgdwvSK9AUBCX8EYawFwoUN7uwzyoMkZdqlpcebNr+DwMQPwxNdPtrShbFV2W1SKveff7+X3NvUqufJX72nAnFVb8d5PPtnVTenW8JxDEt5D2hHyYlXuxWvh2h3Oefdl3UNhB1fEs6Zv0tPIQbk+8bym7QCAfG/YWZQBvpcieOLgCHnQdMVk600Dt/ycQ+XK7s7IeBNNJ3iFdBKeODhCFit1xRjibdiXB3DlmKKe0Wnl/rbcuauYzcwL727Ab19qLG+DuikKRLSHjJPOgFXn4BFAHjJdsUD3hmEbOcaneFsHirIvE1QTMo5iusaNuxNpX7mnAQDwL6dMLHOruiE855CA5xyKRFdYNfSGgcvX+TLq+wH0HMJabn0SJ7YmzmFB03ac/suXy1txD0MRW5J9Hp44OEKeW5lM52uF54YmifuycrFSynbeZRnqHUSWo6Dgz+vzfOAd4wroTWPDBk8cHCHLIquznd9133lkAQBgV8u+e1pVSfbmRmul4Ga/mu4tSS334pRx4Bw8xE2J7ycOTxwcIc+t6mzX2ZPu06asFYqOyb8fQ+9SOvKxYiIOPe2UvErA90ESnjgUiaqM77qKoBTTS8P87inkoOw6h/C/SYfTm4ilDZ7BiuBXOEd0J86hNwzg8nMOQYntuXzF+u+CW9/AF+94M5a2sGkHGjfuci6jUmIlHxbCDE4gfTdF6N4C2G6MbBcopDn2ZflxwWqkQu/Y2pHHY/PW2TMWgYbVyRhGnwkPh1p1w6dSleXKQby9ehtWbt6DC44Zm7j3mxcbsWVPECByb1sucb9QlxepFOC5qAieODhCHjTlGkLFLIL7MG1ASQcAGvpFjIX13T8vLK54xkprX6q63PKdf+sbAKAlDhzn/uZ1fV1+QYx0Ur4rCvBiJUdUKrZSMeX0Ds4hnn73ayvx/qako5YrOnLufbZ2ezOWrNsJANiyu7WQ3hndXk7akxEK29Hso/6awKT/Hp44OCPhId2Fw2gfpg3KxTGfZ/jRX5fgPMPuFzB/k/a8wdBfwkk3vICzb3kVQHAGQlR+5SF+24ZVW/HLZ98ruqxi6czu1n3XVFqLXhCaJi08cXBEQvxTLs6hiGf2Bc5hx952/OSvS9DWEV+0I1PW6B35lW7R0nEbIooNsf7ehohb6cx+JwIuuG0Wbnl+eUllFIOpP3wGTy5YX3S9PRER59Dz51a54ImDI5KcQ9dhH6ANuOHpd3HnayvxxPy4clh10Hs5lNPtKcRKOnRmv5ejrlL0I1f8YW7pDehBEPv7gTkfYOn6nV3XmG4Cr5AuEuWypimmnNw+QB3acwHHIL+LiQuwvbbpdkfOXaykL79n9fu+7CxZbhTmIQOueTQwWEhrYbavwXMOjqiYQrqIZ8p5Cl1XQdt/itDJrm9rIrQdRfSZfCZ1Z9JklVlsWnjakB49f2aVD1biQETjiOhFIlpKRIuJ6Jth+nVEtJaI5oV/ZwvPXEtEjUS0jIjOEtKPIaKF4b1bKOR7iaiWiB4M02cT0fjyv2qJqIzKoddDt4CV+3Ce9iI4h9aO0rmNrkSmm7EOO/a2Y+XmPV3dDCX4ENsX9Hnlggvn0AHg3xhjhwE4HsAVRDQ5vHczY2xa+Pc3AAjvzQQwBcAMAL8lomyY/1YAlwGYFP7NCNMvBbCNMTYRwM0Abiz91cqLhJ9DF5qy7ssoePQKaa7iHFOuYritlva441h3+1a2d+osnwxXnH3Lq/jEL17q6mYo4f0ckrASB8bYesbY3PB6F4ClAMYYHjkXwAOMsVbG2EoAjQCOJaJRAOoZY7NYwP/fB+A84Zl7w+tHAJxG3WxkJ42V/CgqBbr+KyU4pmliF6OQbumQiEMnfHOVaEwnLuuwmOd2rxkU+I/Y8PU/voM7X13RCa2Jw8/mJFLpHEJxz1EAZodJVxLRAiK6m4gGhWljAKwRHmsK08aE13J67BnGWAeAHQCGpGlbpZGwViqbzqF3D0vdHiCmcyhDF+Uc/BxyeYZ731hV+N3SHn+mM3aVqip09dpeqQsjvBSNJ+avw0+eXNrp9RaO4O30mrsvnIkDEe0H4E8AvsUY24lARHQQgGkA1gO4iWdVPM4M6aZn5DZcRkQNRNSwadMm16aXBZWK9dPT2diHGtZgRQmeyzJKYRhNhNZFqvTgW2vww8cXF34nxEpFt8wdqvGgk4Pb5OO6mEliWA2POHyAwghOxIGIqhEQht8zxh4FAMbYBsZYjjGWB3AHgGPD7E0AxgmPjwWwLkwfq0iPPUNEVQAGANgqt4MxdjtjbDpjbPqwYcPc3rBC8IMo6IPvPLLA6rlcXNnueQsExfCMS3m7WuIhJp5aGHcE66pvriNsNpNmHefw82eWxX53x6F8xysrymJ+7ArxvA+PAC7WSgTgLgBLGWO/FNJHCdk+C2BReP04gJmhBdIEBIrnOYyx9QB2EdHxYZkXA3hMeOaS8PoCAC+wbrb6dicnuK7Ai8s24kt3vok9gpcyl+PvLOPJdCr2vhwjoRjx3S0vxHfYncI5KGrRtV02tZVh4sK62fRK4Pq/LcVfKhQ9V8YX73gTv3t9ZafU1ZPg4gR3EoCLACwkonlh2ncBfIGIpiGYM6sAfA0AGGOLieghAEsQWDpdwRjj/PnlAO4B0AfAU+EfEBCf+4moEQHHMLO010qP9lwezy3ZgBlTRyonVaX8HHoKrnt8MVZv2Yumbc04ZGR/AEBrqLCtKeORqYUwBkIHl8NaqSwEpjM2simc/+zWSvp7uTxDVXgmSXdRXMsEq7PMSt94f4u2Db0ZVuLAGHsNap3A3wzPXA/gekV6A4CpivQWABfa2lJJ3P7KCvz8mWW47R+PxoypoxL3kyG7y+UhXZZiKg7eTnFB4n4ANVVFEAfNe6vYe9c+KrUvbYtkp1grKdJ0i6RNrGR6n448Q1VoYN5dxqBsUda/1rw85fIMyz7chcmj68vWhhTxGfd5eA/pENvCQ1FWb9mrzlAxD+luMjMtUJ1FrCMOzy3ZgM1CqGtjuZr0ci9Y5dgRdtUiqmMQ7NZKeurQHZ29WiXT4b4W4nDzs+/h7Ftexbsfli8OUnfsl66CJw4hBvSpBqCPe9/bdQ4c4txpDa15agXi0NyWw1fva8DFd80prnx+XGMsLd2z6nulo6uslXSEzco5GO4VE06k0pDFZDZT3PlN2wEAH+5oKVsbumG3dBk8cQjBd7/ipHls3lqMv+ZJNG1LchNd4SE9bnAfAMDhYwaUp/IU4PM0p+AcqgWdA7+/eos5TIL2tSOlQ5TkKD4xZSuLzqETdpUqAqfnHGzWSvrVNVeGKLXFQtfutHo9flRvOT+L1zlE8MQhBB+v4nz6yztrAQDvbdilGICdP4imjg6IwqTh+3V63VxJL7Ld/Dqts1XDqq14+b1NYbnxe0onMEt5LvO5HOK7bsc5lKKQ7sJFUCe6SdsiTvzKGYjSi5Ui+JDdIaKFTm2pVInYSkvW7URzu7sZKJ8EnT2Af/NiYyFgmrjrM+7UDeVdcNssa53FvGHFrZU6oduVxLFIPweTKWu5I/u+8t4mjBvcFxOG9rPmzTGmXHhkImhrId+UlJPQdSFD1e3giUMIPjDF6cQnV3N7Dht3tkr5S6+TH0XJkc8zZAzbcD6fO3v8ik5TpjVl9ZY9GNyvpqS6VAHQbH0diZVK0znoPIqjMrpm5dDutLuRzuHiuwMdk8sZCLpmJ/R6jsSvFFFQgiB5zqEAL1YKwceEarN15R/ewbmSF3AlFoo2i0dowUGsC8dvTsM5vLliCz7+85fwUEMQPqtY0/lIIS1ShyILE8vtIRpp1eKk9ZAuxVqpCzWvemIn/TaUsbetAxt3toTl6fO9unwTPnXLq9qQ7XKdvfL8bA08cQjBx4hrDPxSF5tL73krkWY7P4BPqq7c2+ic05ZvDOIrvfNBukNqXmvcLJUf/5+uber09Tua8ery0mNxdYrOQZlWfp1DV1orye3euLMF7bm8QnSrb+OXf/cW5jftAGAWs37vz4uweN1OrNNEhJVFUu98sN3Y9t4ETxxC8AHmGvitlKnFGMPz725MpLdZiUP0fFdBnNcRtxUJZNLqQx6du1aZHhMrlbgsn/CzF9C0zR4uWvz0R+0/MHG/U3Q9SoW0Ouujc5vUNxzgEqW2UhDHUHNbDsf+9Hl899GFVl+iz/zqNZzy8xcBAHNWRqHXTESSm1nrNl77wqmKlYInDiH4GHG1vCllndjTllOm28RK3YFzUCn/GGORk1z4CsW2kUn/g/KLLKwEZB1CqHQWdETpztdWFl1mOTiH7XvbEpFrXSBubvjzf1+yITFm5CYuXLsDq0InVfHzmL5LbXWwxJ158yt4csH6xH1vnaSHJw4hIoW0K+dQ/KDaqXG0s8mBCwO5C8dzPiZWilAwKyxxskV6FXU95mdLqjoGlWFA54iVVMS3/PWIXGqxxU/70bP49K9eS/2cuFs3LfKmhXs/wXvalK+OxwgBcN0TixP3PeOghycOISKxkuMDJQwq3W7LtgjwXXlX7nbipqzCJA//84lfvEI6PXhd5TQSUHMOle93VRWV+N7NGu41LRo3pj/LQ7cgpzEXzwrE2yQa4sEFdfm8WEkPTxxC8IEobhhNC9yj76hl5S5wnRzJ57reWimmcxDSVbGXygXbolyJ7siqOIdOsVZyS5OxbU9bqvMPWgTOoTOCsorf0NlaybHDzaHJo2tVvWmttlo7cvjp35b2CqsmTxxC8CHSOUdXqwekbZxGUqXuYWkizjUujit5JxY+rhNfGR+tsFipq6Ba1OR+PurHz+L7jyXFJjqI3Gup3TZ/jd3CJ29ZpBlLjmrTUCqGwVdyDikHzUMNTbj9lRX4n+feS/WcDY0bd3U7H4teTxyWrNuJCdc+WbBmcTVlLQXaWDmWwdEdOIf4AI5RBwClK82LkblHYiUz0jjoZRXDoDt5SP+3YnH664L44TimjU4ximQdZB8gAFi7vRkvLYss8sSFOa/YVABJTsE0H8R3M85YkSgpJl5aTpdzZ3J48VIwe8UWnP7LV/D72R+UrcxyoNcTh9/PXg3GgjDTgCRWqhCd0HqI2nQO3cBayWrKGkor9hYp0y7Fz8GGNJ9TKVbqjPMclE5wybTXJf+QtGitsFhpxn+/gi//LvLlicXkErlPIQpvGie4eCQDfT7xm6m4hGItel9r3Iz731xd3MMSeGiahaHfRndBrycOsqy8M6RKukXGxlbmK7hw6vCGtAiJEyyuc0haK/EzMopBnD9xe+FysuW6GFuVxt2K4yrlDe/OlnbMVThrpRm6xViDpcEu6ejYOOegFk3KMH1P13ka0zkoCEFasRKvtnHjbnz/L4uMeXs6PHFAPOxvmpg+xUK3W7GJ66NJZc64fW8brn10AcZf8ySeXpS07Rbx/qbdRhHDF++cHfu9anMUijuucwggOifJC4QLou+g0Xybni3xvogqhVypM2jymyu2KlLjNV/14DxFniRM62dniyZFgiCemaKK8qv7HQcJ+fS5YjqHMiikexM8cZAiO3bGWNFyDo7WSrY2/vDxxfjjnDUAgL+8oz+kfU9rB0676WX828PzzQUK+PuSDwvX4mlv5eK4VOSvK6avmnOoXEu27WmLEVYR8vd+bmnSux5I6hhMre3sNVHcEJ3z60hHwYSLhJ+DQeQT94/Qv4x4L5dnqfQa6norJ1robqdCWokDEY0joheJaCkRLSaib4bpg4noWSJaHv4fJDxzLRE1EtEyIjpLSD+GiBaG926hsKeJqJaIHgzTZxPR+PK/qub9wv+RPD/6QJ2tc7DJPwvex5YBvac14gRM78A5hlnCAes2LFq7E39bGHAj//L7uYV01WJaTP+VFmHTfD/NQlCl0DmYFtQ1WzXHyzriJ08uxef/Vx3K3LXdafq7s31lVLv2lvZcjIuQc7jqHEyvovKyjrWrG3AOnWIgWQRcOIcOAP/GGDsMwPEAriCiyQCuAfA8Y2wSgOfD3wjvzQQwBcAMAL8lIu6meCuAywBMCv9mhOmXAtjGGJsI4GYAN5bh3VKhkopQXV0ynK2VLOULB7NVZOA93LAmkVb2elKI96LdnKX/rIHqopeoFTxrlY2ScP6tbxjLtuG9Dbu091zHpPwJzGKlTiYOir4/59ev4bSbXg7ao2iTaT5s3BVxrWl23Jt2xUPvp+ccUmXv0bASB8bYesbY3PB6F4ClAMYAOBfAvWG2ewGcF16fC+ABxlgrY2wlgEYAxxLRKAD1jLFZLBgF90nP8LIeAXAaUed8hk6qJoZi2UdXAlaViT6raziQNJD7TFdDMbtT/kTcz6E8C1maUnhMntjzhgK2WJTvtsXY5FYh9oUu9DSQbiyn6YtyEBLVWHhvQ9y7OsE5ONZrytYhseMykeoGjEO3RSqdQyjuOQrAbAAjGGPrgYCAABgeZhsDQNxaNoVpY8JrOT32DGOsA8AOAEMU9V9GRA1E1LBpU+khmFWID8jKEI6S/Rws5Wdc3bwrjGJY9lLWIduzacpW6hwM+fvXmc/NWm4LM+Ho5Sv7MhSLNIrYcjAZLmMhoXPQPHLT35c55QOAlZviZ5nL7Ug7RisxnWSjmO4CZ+JARPsB+BOAbzHGdpqyKtKYId30TDyBsdsZY9MZY9OHDRtma7IT5DnZKQppzSiw1c3ltrYdVbbCtMG1zOLk2kkC6FqM1VrJUpD4Xpt2tybumx639Ukp98V6d7fqLcsSYiUjN2JpUJnhNhbieXSP/OqFRukpfdmys5oc+bhUnUNZxHP8JMPSSyornIgDEVUjIAy/Z4w9GiZvCEVFCP9zE4omAOOEx8cCWBemj1Wkx54hoioAAwCoTTfKDFnsUi7qvWNvO8Zf8yT+T+Eoo6vCGkPIsW1ZUazUhULSUsJCF2HJakWa5qi+hWkRsvWz7TM4i5UsZ364Ik581e/1r6cfnMibuh5uYWdpNmMsVVTWeD739sjnOpjqUHJX0ofsykOTKg0XayUCcBeApYyxXwq3HgdwSXh9CYDHhPSZoQXSBASK5zmh6GkXER0flnmx9Awv6wIAL7AKasw+f9usQvgBedKWy1ppbXjylJI4FMk5RKasFs7BkR80lWISO6j6RRXlsxSxkvgdyiF7lsu0QdX2kjgHK/HQ3xe/t+nMjzTj1aVPXc7mttcT/HdxNtu4Sz6n3ZVldP8w8oFa4mcePaAOMz8S7Wtd2txRxjAa2/e2xfyIuhouy8hJAC4CcCoRzQv/zgZwA4AziGg5gDPC32CMLQbwEIAlAJ4GcAVjjK8clwO4E4GS+n0AT4XpdwEYQkSNAK5CaPlUCTDGMGfVVvz3c8s199XXurK09wwLke4x22TgAzHNXNWtF4vX7cBDodWRvBC++O5GHPjdv2HpepP0MI5rHl2YSCsmNAFvSZp33B062/1O4V0ca0+KMs+aMjKRZlyDLAszATjs+09rzVVNj4vtNp8W6E4dRIKjI0zlsHbiuWwbBQbgorviDpeu3ytNgD65/8R2VWUzse+oarNcnqzwLgXPLd2IU37xUtnKKxVmLRoAxthr0I+T0zTPXA/gekV6A4CpivQWABfa2lIOyKF25RezjcerzzoEP38mUIjlmTpAW3NbDss+DEwTVYpNfchuM/hgtc1LcXOpW7Q+dUt0SMsO6fChJ0Kl56K1pcV6KWbiqN7N9r58N/3AW2tww/lHGAoHzj58JP628EN9nhBTxwxQPF78LnHWii1obs9pHd3MAR8FziGFWMkYcyi2CbIv3DJcF+6gbHISEcllunOMxXOEYruqs5Q6QkI5OIfuah3b6zyk5fUqIVayjIiDR/SPypJMDPnAu+IPc3HVQ/OV5ZvqsFmQcPmmbZFSHcKTBnwnrrfAiZe6QsMKl+JoVeyj4695Em+vVi/AOcbws88aiEdJbTL39Pf+bInD46g8NpuymqsQMc8hzHYkVlK1KZ0+oJixUBbOQeoTWVQkEovqbCbWTpVYKWnAoq58b1sHTrrhBbzxfvoAibc8r5ZqdDZ6HXGQF9ZEyAEm3ks+n80E3AMQHxiTvvcULrgtcIR6c0XkcZzGJNJqrRRSNrtsPUIxCuk9bQFx6FdrZSyNSHH2TAGFKJ2in0PKdeXC29Sim1yeOVNLpfmc1I7/e3M1TrvppVRtS1Mfh7ga6jplAAAgAElEQVRpkBWqrmXIeGpRxD1pxUqOehAT+Pfctkd9NK4J7nXoIRucmDiHKolzUImV5CSdXqJx426s3d6Mn/3tXUPr1FAFX+wKlDb7eyDkb5kUK5kHJBEVFny5rHfCSJliusoKRXsSlqXugs7BmEs6o9eSV4WW9mABqqtWeQm7oxSx0rodLXh79VbsbcvhgMH9UpURDyteHAui5Piknv8PISpnqUZhJrGSWGu5FNKx8q1ipXTKeVU+bqBhy1cM0kRvlblzcYhWZTLxQ4kUxGG3FExSVzX/nssMnu869KvpHstyL+QcLPctGbJEhQVft8jnbeyHViFtrrtgNmfJ9/h8wVGqhEXLbbnSw5U2iKISseTzb52Fi+6aU5KsP7EZMPRH/J6C43NUfMo6HBnrtjfjtJtewvod0YJp9kmIKu6KKKKliJU4ijlcyJlzSKGQToiVEjoHs1hp9RazU12h3rDito485n6wTd9AJDm0frWlbcrKhd5HHCysg3hXNeiymYhzcJmnant5NeTJ0LRtL4647hms2BR41/KBaJo0G3a22BvliDTiL9Upa66T+++LN5jbkXI9FLk1+VFXWpmCpidgipMEAH+c8wHe37QHDzdEAQNclcem40vLHSqlHE503//LIlz3+GLc+FR68Yq7zkGfUebITGKlddtbYt9YRYhlvwYdcRAPi9qy2xxaRe7mzjiN0gW9jzhIv+UJJQ4WlViEKJo0PK98sHtsgCkGrusB60/MX4+dLR14MDQ55e0xzZnjfvp8vL3KHXBpOzJV+1WKUtXju1qSu2pxIpXDvUXcicnlleIUKJZ156srpHKj6z5FiOOMYiWmvgaA86aNVrYBqExcLQ5Xwv/w2024541V2NVqPtujWcFZlEPnIMfIShAH4ffa7c1WzkFOsomVAGBAn2pDC5PoDpFigd5IHKxypehS9ZGyos4hXA9lJWHcWshdPPHye5p4USwYxNFJcO6DR7XmmOTWsWpTOOupDvZRTe4PFKGt6xRB7mLtMN5NwuRt7ApVEfy9m7btxU+eXCrlj57oU2MmDmnpX7wfGfoJ5Y8f6qaPefRfTtTe0x3paor5w8pn3q+Faz+ZiIh8y2StRBQf20pHSGk0ujjKpd3wdHY4dR16H3GQPq78IcRfKhvmTCapc5AHiO3b6j7+Xa+tjP0WF3aRnU0zdFSL3FbH4zu1/hhltj8Xd9pqP4d0kyV2aH3inuk54VpBVflicd5vXk/c+7BEcZ6rZRBjcU4r72h8MHlUPYbulxT9AcC1CgfGoE368jpjASuHZ7x8SxYViRxLv5oqSb9jr0u3y4+ZxKbkBLoJbeh9xEEeLTIBEAekSqyUISrIfQvhLKSPL/4qRXYt5hcHWJrBQxSIfG546t2ColSOaa+v1z7wjc8rsqk4qZoq8zBctSVdSIH4KWH6e8YyFGlcfLjZIkO2LWoqEZyJ29kpcGUyceiI7XzjhYg/66qzeOLrJxvbJeJXXzgqqlMaB20dedz52gr5kbIjnaOdrgyZU4jfl40HbCJhuS5d3eJ8tcVf2rKntDMmKoVeRxzkbs/lZZFQdK36qNkMFSahy9Giyjnv+O3FZ9vzaoseexmEvy5Yh9tefh83Pv0unl+6Af/19DL7g4aKXK2QVP2iNIu03L/i9++4VRhClPeW8+jFds2HTi4Y5nJU8nUTzfrGH6P3Z2CxwIpplJfZhE+PvqGfOXJ0oU1ytttfeR+/efF953qLxX2zVjvlM3IO0j154d0pnWdt1TlIv3ViJXGOmDiHXzyzDD+VfCFcRFWdge5hUNuJSOwkDGIl1UfNEFAd7tw41yHm+8Fji6ymrGl3Brta2rFHVOileL66itDeEeRv68jj0nsbnJ/V1eLOOaTPV455YeQcHJW0qjVX3khwrJQ8xG07XpVZp6uiPM/igRVF6yVr5NeMTBzSmPZG2KnQL1UCmxVh01Uw9XeSc4j/5j49PK+YXWWt5CpWEtcVpbEGY9jZ3IFfv9iYuFfGcE0lofdxDomPG/8tyv3lWPBAsFPjzmF8kosD8L5Zqx0C9rm1lU/OP85Zg9PD4xSBdJxDdTZT2D277jH5Ym2yVlru4NxTTPgDpSgqpYLZtJt2FyslM6rGAwCcKnwbwE48xQWJQ1y3p40biCPGJmM7AaFYSXgJmRsQId+Rz8V2HUdyPpGYf/vMgx1LiVAOgwEAuOPi6UF7UgS5TG4GBfFPjlnDZyT7Ql2vSDRUps33v7kaR/7o75o2dw/OofcRB+H65BtfwB/nfKDNq9opZohQG8rIuZWS2c46mabLb7La2SNYlKThPKod43e7WGYU0hlwxs2vOJRq33kBgVMYRznkrcWuPeLuXbXmugZZK0asJLZ68uh6XP7xg9Rlg8Xamea88CTnYG5oZK2kz3fAkH74/qcnmysW8PGDh+FwRVDDYnBkSEBNm5DEudSyflD4OWlEf7u1koUTUeX7xd/fwzuSI9zrjfqYS91FrNT7iIPQ8U3bzC79J09MnjaXzSQ5B5NMURmjR5P3qHGDpGfVsz3N2KlWhY1V4Et3vulcT7kPYfnmA/OsdaaBuEi+8O5GfcaUcA0HYuccksQh5rjHDGImSRRkcoqTkdA5WPK7cFnVWcIQhQOkqcxyLX1EFLTRqJCOrkcPqEvoERmAmmwGj11xEu79p49gzMA+0bMqayXp969eUAfJk9eEz/72jdjvQX31fdZNaEPv0zmk6XjVIfPZTJTOxQNpZYSlfvzF63aiI5dHlQNXkKUomJhpsr+5IhnJVK9zcGgk1IukzaxP9YzODl8HccF8uGFN7F4pzqc6sVJaqIhDreQ4p1vzGeJiM3HBt+lTsgqdgwtM2aoyGRClCCOeol4bspngjY2cg9D6TIZinMOOve249aVAsX7kuIEAgGs+eSh2trTj0blr1Tt4Kem5perNh233P6Cv3jGumzhI9z7OIQ3U8m+Bc+hI6hxU+ZPlaljRRMRYfdtUC5UqhouYK63HrLadJZiy2sz6yrFuiG8pl1eMQnr8kL4A9AppGTbOQfXt+sQ2IkyrN8kzBiLgjWtOxVvfOx2jhZ2uDQni4Njb8uuIv6stZsgyiKhsFmQZCgilqTxxuGUofq7ERXfPTuSvq87inCMDr3N1VFZHrtkyVEycQ3dBryMOaXYtqkGXJSrs1vguxLRLSCNWKrVtn5NYVyC0wHAvNiw7gE7GXgrnYJtc5dA5bNvbrlWqOyukiTAo3N3d/A/TALhzDnaDhKSBgGieCgAZzcxkLHhu9MA+GNa/Fp8+YlT0jFBg48ZdmN8UP6wpacpqbqeLBVW1YNrtgnJyDpkMJbyaE4gRh/jcW9CkPswqip3mpjNTwSaCLDUcfmeg9xGHNEslQ0H5zJFRhOxOe96yboDJyXKIBhHuDkLmtpjw1fvUZq+lHKquEyutCcNqlCusTKkLEAHYf0gQmoLv8lw9Xe0EUJUaJQYEQKNvgqw4V1+f/sukwYBKP7Fxl4Nnt+F1qrKZVPwoURnFSkQgyylz4r2Ac3AoNxPf/IlwXT+s5zh0F8WCAVbiQER3E9FGIlokpF1HRGulM6X5vWuJqJGIlhHRWUL6MUS0MLx3C4UjmYhqiejBMH02EY0v7yvGYfomYyQWPc+S7H1NVXTObMFDOuV31g7mFOW4L1TR9QNvrdFnTIFJw/dzrNtd5/DOGn4WRvpJoyJ6vO6kWMm9nLsvmY7bLzqmEDgt7bnJMsYNDsbX4nXJHatYtMn/gDFWtnB6jAHHXv+89n7BCc4wMMX54AYqm0I6Q4SspEcAAr+Cxo2B+WieMXzt4wfi3R/PCLkMe+18zrsE3pPNg4FAl2E7x8Hsm2FtYqfAhXO4B8AMRfrNjLFp4d/fAICIJgOYCWBK+MxviYhr2m4FcBmASeEfL/NSANsYYxMB3AzgxiLfxQmmfhdlss8v3YA7Xl2ZMDscUV8bHZ8Ypq3bobd6emtVUg+wSLE4BOXZRwWPQ5TGwaxcm5QTDxoCAM5yblW9OuLAu76Ytp43bUwiTQ5SOH5IX8z9/hnOIhAiYMh+tThzyshUIdrFOmVwIqP0HBfS8oyhSpIrMYHY6V4hLdGwH2xlL6NfTTaVLivgHMqkc8gEp7fJ4r7rHl+M03/5CjbsDEJw12YzqKvOBvoJh7ojziF5T366I8+wZN3OWNoNT6u5htaOaC0xtaPH+Dkwxl4BoD6UN4lzATzAGGtljK0E0AjgWCIaBaCeMTaLBW9+H4DzhGfuDa8fAXAapRFipoTLqVGMMfxxjnqXTSErK5b1T797y1inHNJ78dqdynwuY+JjBw8FkE7EUc4QEhnH3VdQdzJNp5/hfVquuDJyOYP61WBwvxr38xyEnJRRl6kDzyZ/d/Npb1HZHXmGEw4agjMnj0iUyRhLJeN3aWcp+frWVqXiHExZbTG2ZGSIUJPNJOT7/Jjenc3tQdvDBhK5WRZyI0AXzgEALrs/Ln7dukft2f20cDTrvsI56HAlES0IxU7cQH8MAHFVbQrTxoTXcnrsGcZYB4AdAIaU0C4jXPpdZuuX/ijOOPFNnbtyKsn2fmT8IPzo3Cmp21ZTFXAOaXQOL5bJ1j9DgfLR1RlMKVbSPBuJ6tK3S3mgkqyITllmzJdA0jHZwN9BPEZULEcFseyW9hyyGcKPz5sa3RfyaUtJ+ZK214nESnoEnIM7TDqH4f1r8a3TJzmXlSUKOIeOeIFy8Zwrla2VdChwisrBmEzbT1Iuq8LXA/H5YA4z3j2oQ7HE4VYABwGYBmA9gJvCdJ1xjslox9mgh4guI6IGImrYtElz9oEFLv0uN1iOzx/tct3qlHf5jAX24ZNH1Uvp9gK5U5trdMg809tipwVRMNFsJp0F8VwqziF8pMycQxprJdJc246F1dX9hHhcK6D1XQDiXcX9IER5tmh9VS6+2uoh7VBRXXU2JeegNz3NCFy5U1kURABol8cji/0rlBmY0drBx29y3jLlXKqvi/ssiMShf11EOJw97J1yVR5FEQfG2AbGWI4xlgdwB4Bjw1tNAMYJWccCWBemj1Wkx54hoioAA6ARYzHGbmeMTWeMTR82LOm97Nh67R1xgaLCbiOZL1r73D6jPCjyjCGTcQ+2JqIm5HlNh5bIdZUTBLJaKxVktipTVu2Zu+l25yJUj+jeu5g+l3UOx00YbG5PmK+lI4VYSWhuc4E4RNMzOj6caRfQ9DoHx3zyRkS4DsZjWp2D/n6auEtEhOpsJhl2P/z/w8cWx8rMOOo7dArpvy5Yr8wvO8vuFE47HFlfV7gW52iaA4q6CkURh1CHwPFZAJx/fhzAzNACaQICxfMcxth6ALuI6PhQn3AxgMeEZy4Jry8A8AKrIF9lKrmwQAHYGJ55oDryMa0IRJaJcisoeSK4iZUyyro7LR6Lg86B73hVuXSERbYAKxXiYhqUn3LpFLLLbbO1UOXHAJhDXYgbjfGhCW2VEPqE3zdxDi7v+MrVnxDaac4rG14AQMOqrTHZOfc1cIUpfAYDSxUOBAgi4j4+f13sjBLe/7NC3QMvs1RT1u3NySNuxfwcrUJgRXFDIM7RNPGgugpWTwwi+iOAUwAMJaImAD8EcAoRTUPwnVcB+BoAMMYWE9FDAJYA6ABwBWOMq+gvR2D51AfAU+EfANwF4H4iakTAMcwsx4vp4NLt981ajXc+CEwruUjpr18/uRDELtrlun1EeVefYzwuTHxQuRTHOQd54OokPcUstrpHOEGzKcMLToIqzkFTOJ9Ef5m3TnnfBFWRshNcap2D8ERG/t6WLi2IM6RKTRFUwYCxg/rgv84/AkcfEKjwYsQhLDTP9ERAZVYpY//Q2zvWUA1U3NwFt81K5rPWKuY1WwwVKzJbuHY7Tj00UODrhmdaU1a5HLF7X7/mVJx0wwtBudLzIrEQiV1P4xysxIEx9gVF8l2G/NcDuF6R3gBgqiK9BcCFtnaUC0bOIfz/hhAxsTZUAE8VIknKCsoJQ/slYvqLkA+JYYwhS0kWOh3nIBOcyouViBzFStnkosKhIyy8K5auV1typYVcTWrGQaGQlrkRfd3BfflddV7PQZlBnSdOHFpIq1Y+oPdzqFNwuSaoo8NGcNUDpeLKDJwDgcwE1IA9rYKZqCYMDTlyDlyv15aL94+4YRB9omRxofjZxPBnovjL1KXdhDZ4D2kVxG/dV3FYPL8dOcGZy5QtdLhYKSE7dlJIa8RKWp2DtchUSMM5qBaVxzScgWnhtMGkc0hFG2NB7CLICmlbmeIuX4RZ55DUJYi7TrFuvVhJnc79K2T8+K9LtO0BhE2QMVdazkFfIANLBJPU6aj+dPmJsd8vLosUxfL34Rs8V51DfdhfO5vjVkc6xkwWhenO24ifBW7iHLoHeeh9xMHU75TMI1sqAUnOwRpoTdY55LlYSWqbcL1ordpRjnMO8gKtm0Smpp03bbT+pgKEoN06zoHb5Wc1OgfGGF5+T21lRiA0O0RflU/Vam7LKSeT/E3SBh2MPZvgHMywic5U4JyD9j6z59PRbJ24ab3BeRNCPbbxnU7nYLYYkkPMP7FAvZk45oB4ePu4ziGel+sNAyc4exs5MZXPl9a9p9y9bmIlff26e4/NW4tte8znl5cTnjgIUH37uiqTQjrkHCyONYmFnDHlDlxs2yvL1Ysoj/WUOHRE82KmXcjAlJEhAwfAJCHiLDafFNUavchDDfrwHUTAq5p3FrF6y97C9RuNm3HYD57G7JUK4zZHEZC+PdIuXth12nZ2urtGU1Zm3oGzQr5kSBduZ69rV1olb+E5xWZJBdWiKcckA4J2Esz9J3uG63wGZAzvX6e916cmKNPVgbM6m0FddQa7W+N168Rn8vcQiYPIOXRYdA48LI3q3pqte/HNB+bhyj/Otba/XOh9xCHlYlGnEitJVhy2xULe7TIWDCjxnIJTDx0e95LV2ETzXaBMDHScg2kyDNYc0qJ7ol/oDSsTNa5A5ZMkMmWNP7/sw93atrhC5OS4NYq4a+SIwmeECcUbKwGIO1Dp+ufz08eGdRbLOSTvjxpQFyszryAiL377FIwf0lf7rdPK8c8IOUBXr3UVV6Y6gXDWtacarZXWbG2OKeEBwfLNSpCj+2u3xzkivsEjRyc4QK041/Wi/F1FMV7MWkmYEKpmHL1/MI9E3wgOHnpj/Q6HQIllQu8jDg5jQ8zys88dnrgvWisxxrDO8sF0fg5726KdCUltk0MvcGQLu/J4ul4hbWxaKowd1EcpVuLDn+9Qdaasrk5gAPDZo5LxkuQyTN9SXshTWytJD4hmkLp6o7Aq6vtGU1amVjRfevIEAJJIS2rcsP61mDi8v5aDlU0tC3Uq0l79zifw6y8eBbEa65xRFK+qsqYqYw3ZXSMRFa6DaNPMBxfwDR7BfT6o/DF0xF1OrhWkDSIjFCcOyYb86LwpGD2gDiPqk1xQsVZ3paDXEQcTVDs3OVIrEGe3H5271lruHok9DQ5sIRwysr9Qd/wZ2cKJQ2cmqj/L1tq8BHQDsLYqgwwlj97kbeebvoxmt2daHGXOR5ezlLMkSoFoBqkr2eb/YtrBB4u+qsw4m6ojIiYzzTTK/hH1dYXFzdUxUdUe1SmFBLIe9qPjHNo6LMTB0Ma6qkjn4Co4IEWRun6UiUaeMRw6sj/m/eCMGGGOh89IllNblcWUMQOMY7yCYecS6HXEwUXn4Hrwep4Bm3erg2yJeFmSpedDsdLYQX3Rv7YKowuigyiPbrHn1lPyQqDbNRptym0Nl1BlOdglwTlIVW81KNMSr6tTujpSB7lulzklZpGDwMWUmZo+jXbabos0k36omqjyxle9i4krSyNWEpXXrpEAVGNCxa1kKMkhJ+uPdxIvp9VCHEwt5HMmk3HfNJBCea3jHBIm6SwQDQ3sWxN7RuS4dX2q08mUd6vjhu5/HFGZYT7S060McdK4PPO/L6/AtZ88LNYGvkF65wdBGOmv3f92bACIxOG5qz6O0QPrsHLzHmzc2RqWEa+jnH4OWoVqRm3zUxArJXQO8ZIeebsJOqjOzlC2LSXnEHkrpyOFfWviUyNDEWGykVsd/ZLfT1ws5q3ZnrCOiUqM61BU/WQKKnfMAeZwH7FyhJVO5oS271UTdyXnoCAOFDjKGL+hbK1U8DmwEQdDoQXikErn4K6flDniYF2IzwUgPqd/8+L7sWdOmjik0EZVE71YqRNg+tzOC4gwaUwLmu5QHHEhrMpmgoPSST/AJw7fD31rqjBl9IDChE14XUu/jxgbOO2lIQ02jolHZZUhp5nCZ8jQKRz1R5Tad19BPofKU0DWORwqiAQ5xHe55/WVyjJEiDv6tdubE9YxQFy/tW1PG15r3Kw0JTaFhvjp5xK+pwBQiAKgg+wZftn9byvzqaaAnnMwzzFZHJUhwnsbdlnNbk2fmzsHyk5wxhDhKXQOSbFScsMEmAPv3XHx9CC/hruJwsDom1xu9D7i4LBzsMn1Ch9cKkteMMYP7Ycjxg7A/oP7xtK5n0OsTmurAmQ1i6k8oG48/wgcOLRfqkWytSMPxhjmr4kWjc8JiuFs6Jshz/todxtUltFwDiroLJtkCy+OtJFRuSWTk1jJkCeuc2DKA49E/47rnkg6mMn91q/W7tEsqhz+NDfgvMTvo2qfjFqFOXYa8G+zXHO6mWrBV/Wl65nUIqqyhDNvfgXn35oM2yGi4HOkGPB9BIW0OG8OHNpPW57rwUBBXrkx0RohEgfT2OWcqs6iKuIcvM6hYjByDukYhwTn8PS3PhaviwWmiLKXNQv9HOS6VWNHZs91cV9kzoHIvGB8+cTxifdt7cgbw4BkMhTaiks3wnL4zoibMbrMrSoNIdGJEdKexrYq9IvQfVvnBSATLRY6nwQuAtErhiXxg0PVolJYZbJbKNvRwSsNZM6hpV39TVR9q7XoMnDIx44fnOAcsgotsIorMUUriMRK8flv1D9Scq1wdW4MDE6Ca3H+ygE4lfVq2lUgDp5zqBzU9sUDASBxvoIO4qQxfay2XD60l47SNu9uxbodLQkdgS7OvTwRdGKlBU3x3eTetlyoVIvnO+fIyCv6+APjZyq1tueMu5usQixwxuQRhTS+29c56qkwZlCwA5ft93Wmi7Yy//OcKYXyXOC6oMpiJbUoxUwU5TAWTqFcCu00RyxNc0KfK+TwGfKY432tapXe3Fe/QTt87ICEtZL8vf/fRyfgpW+fknjWdJ67aK0k9pEpkrFqkZb7l68XchQFBoFzSLnCZkjtQV7O0xyd29LpNXY5kp382aPH4p3vn4EjxwVEwh5oLPhv0zm0tOfCHUhU3o9CccNLy+IWTDLnwK91cmq5jf/+p4Wx35w7kV9FdLCZPj6uqGxpzydY/8PHRgEHuW5ExD9MH1dI48ShRhP/ScZxEwbjFxceqXwfG+ewq6Udb65IekZzb2HX8Bk8Vz7P8L0/L1LmAeKLbyj9TeTRcUEcl3/8oNhvFwImipVMUVczRM5Hx+pw6qHDlXUXAglK73XJiePDjMmydHNIxyEDwfiS/RzkvGdMHolxkpgWiLhW1SLKiSoRxaz6zMYpyc2avPG/558+AiD5XUTOQRx3LrRbR+QjzsGLlSoGnehmUL8ad7GSsKMyPdPSnkuw+9zTMeF9KZXDB0g/6QhCPtBt5zccvf8gpXWGkZh15GLzfPyQvvjyieMLnFUgVtI/364RK+1qacfbq7cl8p931JjCblpuZ4smYijPd+Uf3lGWqTsc3vZtbU5WFOMcgsn/D9PHxfKYDjkCkl7DTmIlwbFO58zG6y6VONz95Y/EfsvhM/QRdVViHnUdnENWEX+ipJ+D3Jc6JbJLUER54TWKlRT35Tk3vD4QGYvJ2/e24Z0PtkfjV3gdF85Op3Pgfd+JUqVeSBwUaXzSuXa8aM9uouTD+9cC0oDkE0m1yIpt46dLPfS142N5ZJ3DorU7NBNNrZuQF5jZ3z2tsHtvbc/H2pUJ/RomhlZXKnv50DoRgMA5SGHFD7/u7zj/1jeS7ytwXvk88NX7ooPaN+9Wm03yMt/9UB3am7/fL599T3k/2QZHnYMkKyckJ7uOMHEMksKVuJ0tELYTzMg5VGUzBeJcPkTErmnbXn0uRbP0J/EF333uBwrCTpTwc5CLUcVsAiKzYJupunjbRExVOgfVWJHFZN94YB4A4K1VwfvFFNLCuOBhUWxtLLQ1TPQ6hwrCqPR17HkxKqvuiVu/dDRuunCaVsaatFaK6wdyOYZ+NVkcOCxuDiuGj16zdS8+/avX8J9PLNa2U7bekInSiPo6jKivBRByDsJtTgz4pjqTIaUMVRYr8d0fg3vkT3lX9uFOdUgS23rK3++5pRvMGXl5KcqVw2brCK+Kq7vwmLGJNHmxOf2w4Yk8MRGmgThUZymm8DxCEAcWC5Fz+Mo9bxnymTc6IkxnKhw4bL+EoYYr58AXejPnEBcVmYlzUsGvYi5FZ7lNu1rxihR1WHwfcVzsV6t2MavOZJSbPT6PPXGoIFTUPzXnEP43OW598vBRGNC3OiFW4tlrJPYZ0k6lI880tuJ8d8oKTlNzNfbqmUxSXKKIaqA9UD3iUoL0LCXPoCACtu8N2sHtyaszmcJO23aoeuKUtRAjFfFlVG2UITtRcejWAZ5u28XHFNIIzl6QnylYKynaKJ8zLLeprjqT2AgAoljJwjlkorOU1+9oxvINpQc5FC2lVA56Ub5kmj74YHBPdfv8o5PxtOSu1HEOLt9R9gWxWSvJJE5F9EWLQKWfinAdD5+hrry+T5UyEi1vtzdlrSCMewXHfhetOH7wmF6JycsUJ0p9XSBjLyjzeD6pcbl88uATIO4XwHfosk8A35VkiBKLs2r3KRIBccz2De3w+YKsUkgDKCjyDxkR+HlUV1FBduq605cXgYf/+QQ8eNnx+MKx+8fSbYv4foLC/XXhRL8Lpyd37mnKjfk5sPhvDm6tpFxELHJ5rXGDsHs36bhjhbEAACAASURBVHuqBM7hhJ+9gObQGOLV73zC+F4miCKttOoMbXypDCHHmLK/VSJamcjInMODlx1faKOp3qD86Ju1tOeMeiaVzqFdJb4VrlXE+8pTJxauxU2Dbrj1r6tGc3suMadL1ScVAytxIKK7iWgjES0S0gYT0bNEtDz8P0i4dy0RNRLRMiI6S0g/hogWhvduoXAkEFEtET0Yps8movHlfcU4VB+FDxhXqixaccjeqnO+dxpevvqUKC+SuoT6uip8TbJckQ9B0XEOoikrH4yNG6Nd4sUnHIBn/vVjhbrlQaYLvQAEfSOGEf/tl44GAFzzyUMxY8pInDllROJ5AuGfP34Q3rjmVHzj9En4+MHD8PVTJxUspWwmeKKYTMTogX1w3IFDUN8nzn7b5kj/2shc9Et3zi5cnztNHeWVt8829UQOkCH4DvKEHd6/1qmNct2A2veF18vB+4hHahVRlUluBBiD0rLHFaJIy6a8lWEKW57PJ+cNByfuh4VmonIu2aHvuAOH4CPjB0XyfJtYKbx/6PefNvqNqGT/KsdMlbm4iInDI8fY+GZA/Qy3JtzdkgzWydvVWXDhHO4BMENKuwbA84yxSQCeD3+DiCYDmAlgSvjMb4mIf81bAVwGYFL4x8u8FMA2xthEADcDuLHYl3GBarHi39zdWiksS/F9h/evwwFDIs9L2WKopT2XsEACkt6buXxeuRMRPaRVC/0JBw4pRJIlokR0V5VSWVygxcNERg0Iyhk9sA9uu+gY9K2pUi4E2Qxh9MA+qK+rxr1fORYj6utCsYsiswAxBo02r5RuUyDzg11cUVjwLf5JMVNWxhL+K0B0OplKrKRW2savVd9TFGHy/N8+85BEvqpsBh155qxgd4Ecml6fL5lmii+VZ+o+AoL588SVJ+O/zj8iLCeeTyVWEs1ObQrpVOc5SINPZ2HFc1mPC46JldR5OGckczXd0lqJMfYKANmg/FwA94bX9wI4T0h/gDHWyhhbCaARwLFENApAPWNsFgtG2X3SM7ysRwCcxrmKikDFOaTseFEO7JA5toh8uLNVeRC8bB1h1TkwnadodD1vzfZEqAWVWEncIa7YpPeQFvPaELyPi1gp+K/ry3rJccxEcEbU1yr7VsZr/54UtaTTOQAg4OunTcSEofGNAKAWAZxyyDAAwPWfnYqLTzgg6EdJBq3iHMSNiGn3yJ/9w5wPjO8BBN7xLogIk42zUo1D9RPZMHaQ2EcvX30K/vwv0ZnQh48dUHAsk/tS9oMAOOE21xvkc/ciV3EOKjGU+Bl13BCHSBB17eTvJxMiLjKc36Q+PrgSKFbnMIIxth4Awv/czGIMAPEsyKYwbUx4LafHnmGMdQDYASDuultGqD5JWjMx2f7bBHGH+VDDGrzy3ias3pJcgGUZpyg2UtWd08jzbWIctTmqO7FL0G1Nn/HJZSsxMv+Mcp59+MjC9cUnHBDLzyfVhp1JkcCzV31cuXjI/Th2UFLUYmtnbNfJgtc+dGQ9XhS8dUVxo4j3fvJJnHpocLral447AD86d2poChsWF35L1Z6IE5yOfF7rGAkAH2wNTE1NjnwcR4V+KzZE9ah1BBxyeJhfPLMMu1o6cNCwfnjyGycnysxJYqUDhvTDUfvHz4TmVbvozERLP9N3TMc5xMvavretIHL9w1ePi9ojcC02vYDOx+I3Xzy6cM05B1mEZT3PogIot0JaKX40pJueSRZOdBkRNRBRw6ZN9vOGVVCNjYhlc6MOhR2iUNi9XzlWkzdadF8OzdxUY0heGGycA9Mo9cSyv3DsuMR9k1zbZd6Yzg2QyzSJI6aOqcfZU0cpxUr/evrBhev+dTLnoG9kfV21cvE4dFQygipHGmslxhgWNG3His17lAu5qA8SoTK/FBcV06K/cVdg0vvSsk0FAqr6BmtC4jAs1HuY4MqYixyl6RwNrh8YEvpx/PrFRgDA544eiymj4ya1QYwqt/4G7M6eQZnhRoQxrfMkL1NVnOpAL9FElTGGaT96Fne9thL9a6tw4sShQr5o7NqIg7jei+8vboa4o6Tos/L0ovX431dWGMuuBIolDhtCURHC//xosCYA4oo0FsC6MH2sIj32DBFVARiApBgLAMAYu50xNp0xNn3YsGFFNVy1sy4sYOFksJ0rkAmtdsQJM2agwakFwI1Pv4vnLbb3YttyOZZwCBLblsureQQxAqzKlpovVCLh0SmFVUiYsmryFWTLmiL/+vWPYlC/GmXdk0bYF/NkfcF/FWdkMqd1kVUH5QehF8759esA1O/N+2ajQdFZyKsQhagW/RMPChai+j7Vhfyq8XldGOfoSAf/BlfxqbhpkHtHrue8aaMTurSDFKa5GQoWfJsIJqMhtLp25hnDXa+txMk3vmgsU/7OXzlpAl6/5lRlfj42xIOGkoQ+IiK2toobpZhjrPA9I+IQ1fnP/zfXGl69EiiWODwO4JLw+hIAjwnpM0MLpAkIFM9zQtHTLiI6PtQnXCw9w8u6AMALrJxaNQmqkg8Ndz78E5mcjTiqMvJZyupnuD38rS+9X4hqOahvtSJfvG06zkE0ZVV107RxkchA9R5fOG5/fPnE8fjm6ZOifIrduy6csbPOAVyJ6rhDdDTxyTOmDB0dBTqLGsjb6nT+sKV6WSSh5sCC/3e9ttJanbgzLSz6ikIHCmPFpHOYMnoADh3Z37qx0T2vzCfUK3/GRKgNx/Ad2XAht53oxwmtjYhw5Bnw53fMR/aqfFNUfj9A2Edh1lYhGq2sEBf9Iaycg4NCunC4UQlnZpcLLqasfwQwC8AhRNRERJcCuAHAGUS0HMAZ4W8wxhYDeAjAEgBPA7iCMcb5vMsB3IlASf0+gKfC9LsADCGiRgBXIbR8qhTEb1JfV4W3vnc6PhIGoCuc3uQwe+RYNjp6wlleEXKMnaDyKN//PLcci9ftSMSZCdoY/E/jQyCib3UW150zpeBvIZfJcdnHDlSWmfDs1vQV1zm4Eoc0i8CXf5f01r3p80EIEJGg9gtj5K8yhCGXF2hTO8U8LocemcuLiLtp0RcdFHlMJ109pgN/Yu10NtmOOIdmSVwzsG88FEhVRhUTKNkYThR5jDF93cH/nMMiqYtkCgBPffOjUT7FXNQRU87xL1q7A/e/uaqQ3iZxoS4KaS5yFoNt6va/OoV0V8B6TChj7AuaW6dp8l8P4HpFegOAxJFUjLEWABfa2lEuiB8lk6GYjJYPE1OAM46qTCY2GPQLQ3LgKv0X+E4pl8fNzwVxgUYq4q+IHtLyWJTt39Vmq2Y9Bodusrkuf9z+27ab0snpddARmxPC8OPiO48f2heL1u40Lpj8lpWIZeL9o+oHV31M8Hw8HEfwvJ5TzIXf28QZ8DZOGV2PxevUsaeCMh3bGFYlf8dvnjYpMYazCU5ajTdXbAGQjCIsg3NRLmWavK4PE8Lw89AdeYd5y5Xcn/7Va7F0+cx40apJN4Y/fnBSBK7lHDQK6a5Ar/aQlicaF/G7iJUyBCfOgSi5S1Das4f5RNbT5OcQmBfGy/2cFH5AKfqwmMdyaGX7jiOGH6yiCkInWmekFStpnasyXKzk1j6OhxsC4zo3U1Z1nh98ejIeu+KkVKENxHDqJp1DlkTioDZ3ldto60vVATq68oL2ueWVRUUThiZ1Dk3bzLG2ovKC/y7jQhfJVFUmk/QdJrGSy4hMY60kwha51xOHroDwTVTevkBwqI0NVdlMbDDoFgbVmQpqziFoGjO0DxB22owlFl5ZgZ3eKkXkHNSD11UhzSesSOy+cdokzP3+GfjUEaMKaWkWAUBNbIBoERU5B5f5xVl929oS6FCi348K8u2vnDwBR44bCEoxm8QAdBFxSPYm/6acczB9U16mbbftyjmI4TNs4GExAODTR4xC/9oqHKI6Z9tlciGdziHyxrflSxJPrVgJrqbq0bh4e7XSjiYGxhgem7e2EI9MRiRWqpja1Rm9jjiIA12eJHycuCyqMhute4RcxUrhABcHrmoiiYf9yJNWfh8X8RggWitFabqJkcYXRH6f6gxhsBS22hTJVAXbUY3iO+ccjmXksnQrcXDYnaaQKhWcBIGI4Cn9HMJvynUOZs4hDOlioYouSmveRsCdc+DfOs+YUiQKmA8sipcX/HfnHOxl8oVcNHfV68zUeowF152ZzBdmXL7RHuwwz4BvhmG9VaipCtqj4xx0UWkrAavOYV+DOL91imcXY6mqDMUWHy1xoOSCpsrL3fXFvCr2X9Q5/HnuWuW96LfxFYT2pNA5JBTS2lJDU1YzAeXluesc1Om8q8T26c48FtEcOjbZxUr2trkuukB8ZxpFvU3mqxKC+ZmiAPP6GQN2KqJ6inDfNPBvE+9HneKci5Xac2pLOyDJ3erPNXDnHIK+dCMijDFc/ciCqN1GziFZZl+FBz4n8hM0Fn4i5HH+xJVxJ0GbWMlG+MuJXsc5mCxOrHF+BGQzFHNUMVmQJMRKBnGROA/l2PC8PADIMeBOyWRSnpAuuhOxzFg7NZPNdfnj8t2Y6M2wsN3u6ORjCugmw2YRAwicgyWfSedgasP/XXqcIid3Bot22jxNBv+mO5rb8fDbTcb+JwSLyra96oOS5DJt4EHuZMuZqaOTvhRVglgpiCisrkOse3j/2piHuYjIUMHRWkn6NH1rsrjiEwcp84lnfWi7QqNzkCMliyavJvPc78wI4mHJY+hwyV/ERhxsDonlRC/kHAR5oyxWCv+7KLeqMuRsbqY7MSzZNjermSBvMl+Sc0grVrJzDnLT9boWSoiVyoFA7p5MV/WpC+cgikJM4E5wJsjt+sl5U3HypKHqvBCd4KI0Gfy9bn3pfXPlYRt3t3YUFkpdDCXXcVEXnkEh9+Ppk0ckyxRMuzvyagdOIH7exn51VdpYWGlMnDNhvCZx03XD+UfgnCNHS2Umv7NuA0WAk0ZaVFybRKOiYYHYHhmcOMgmswDQv7YKu1o70JbLoy5jjyFWKnod5yB2+amHqA9UdxF/ZzPktDNV7Wq01kqID7DPHZUMM22y7kmYFzoThyTHpNc5uMur84zhwx0tsbRSoVvEVdE6TaEUPjl1ZKw8m1hCJR5U5TH9FqGyclGLEY1VJurjO87rPjO54DUtw1G9g9pw4XYZ51lB59CRU0cUDtoYpZtFZMH/nMPRp9ws+N0PI+dInaWePG10bXDhFHk+VuCYgrRshvAtwcmUpwHx+S07EgJmP4fB+wX6OtO4Lid6IecQ/L/z4umFSJkcfBfsstutymRiOyrd4kKULE/NOSRjzsgsp/isqolyua6LsdJaSfM+b6+On/2rq4M7JomhJA4fYw/tcP7RyUN5Xr76FDRta8aX7pytnbAqotVq4Oxkeb9tHRAVrjokLbn0H0BUovIFXXWKHRElHC5NbeShQlxFRybUVak5B2Xd4cLLTUV1YiWblzlHKp2DQgSkC3+e4BwMusK01kp5xlBXncG7P/5kIl9hUycQOzluGBAclAVEY2JvW6Q/GhBGKN7TlsPA4o/pcEav4xz4MBo9sI9afgg3sVJG4hy0O20ojt/UWCsBcfNUlSe1qY3yhEireHRRSIv4p5PG4/gD9QF084zFdkAfnWSPhzV1TH0i7YAh/QrB0fKMlZUDiYiEOT8XXRjzpOAciKL+joiDejo6f0fhWFiTL4OLaSoQcQ6yd7QKnEvNMx5RWF3/zI9EJ/uZOIdI5+CqaJbSdPmkNL0pa/I8hzsunm4sM5dnWm6df0NRl6COIBykccXz5B88U7jHoxrsVRxHWgn0Ws5BbTEUwFXnYNqZcmQy8YPfAbVVCrdeEetW7SRFayUZCYW08JI3XXik8pChWJlCkSobdRH9a6vww8+oxRZAIIqYvWKr9ixoHUxsPqD3c0gL/q5N25qxbU+b9ZvPX7NDeUawiISBgyGv6P/Cd8daJW4K8WDEORgyuoqVQs7BZTHi9eXygSltVhhrnzt6TEG38MXj9sf1Ty7BnjYzwUmlc1CI/NThz5McscnKUB4Spx06PJGvKkNoDYlnLs+0OgyeLh6+paKfKrNyDn4qoq3vyoVexznwPlcNiipFuFwdspJC2sQ5yB9atRPk1h7iTknFzmcVC7l8j0Os5rBR9ZgxdSRUELmRMQP74MSDhhSigepg66HNu9uwdnsztu4xW87o2mJqY3tZHISiMn785BKrcYGNMADpOIcNO1vxwFuBdzbfTep2266cA1G0ETFzDklMHpXk2Ph4EjdBZyiU0WJ9/OhccWPzy89Pw08/e7iQN7jnpHMo0lpJHRgxqXN4cuF6bbli1myGlAv/uMF9C2dp5JnehJf3pWiKulexyEdi4+RX4gS2s8xZex9xKFiGKBbobJL100HmHLSsumaQytjdmsP2ve1YIJz0tHZ7MtRAgd1WiZUy8u+oHpNnqqiIb8/lY2G/ZXz37EOD+h2tkFQTwAStSXCGi76AD3e2KPNwfEwRy0aG2H2Pzl2Lc3/zunsjAfz7jEMTaa5BCWXw3b6KUwzKcWtThiJFpqsXNMevvnhUsrxMkjgcoBkbvL6OfLDBMRE0vgkzRfFIwzmodAk6zkHOJ8dKEp8Xs6pEQEAQvHNvWw572zpw36zV2vDw/HFx4zl+SNIvwmROzzm5clsA6tD7iEO4iKsmXHXGPSJiNkMxqwHV6WKAOQSGiPc3Bd6V1/wpctDRHWiuYo95m3R125yngMAKYkdzu1b2DURycdcBmtayQtdK/mryYlFTlcHs78ZjQE4anozpI8P1RDAdxNDoOriqRsqmc6DIa9/IOSheXbf4yVZ5JgMEAFjQtB3vfrhLywXxMsVnTChW56C2+Ujm0zEm4Vak8HvsoOSBQEBkwsv9dHQcJn9X/q2/d/Zh2oOZMqQWG4sOkZ2B3qtzUNxLxTlkI87hf2ZO005g1wWiVnGwuG5i6MzsEgpp4bcpbAF/7idPLg3yGrgM7r7vOkBdlJmqtujSResNABjYpxojJL2Gy1pa7OZrQJ9qjKivxQkH2U+ydeUc+G5S1+9p/FX42DXpKfarS057k75DPM9A9058/H/xjtlhPn07q8O8pv5Jq3PYtCsZLVWVT8a1Zyc5QP68OMTlMcbhej4F7x/OWZj6J5tRz++0G7NS0Qs5hwBKziHLd9AunEOmoIgy7ZLUiu9kIm+XOCZOP0wt3w3MKlVtksUa+nvx8uK/dbtIIP0A5QvL3V9OWnoU6hN8FGw6B9kIQNm/DotpsdOrI5dPHH2pg8uSHph+mjkHd6vUyGvf9L1VXI8ufyZjNgnWPd9sECdms5xz0Jen8nM4eaJaD7Zw7Y7E4Tim0DMidHNMNo/VDSkeY80W6bZgrRR+a7OlVqAbWSeJlTnB8MShQojEMSqlcNAdq7fuwdD9ajD/h2cm8kR5I87BvPAm78kHsgcNC/6JOwbdZJDPFuCQd4tiu0xtlBdTE+egcjZTgVt2tHTkcNioepx6qHoSAvFFQtdM3o+yyE9FmF022ianN/GAGI4rPzExqD+XT+0/YkJrR76wm3RxHDNBfNwmipooid60ynCK69b0Yr/4HdEhTQavy9RClZ+DThSoqkt9lonYBsJ/fOowrYc2P88helbPMeXzrOCfoEOBEwq/tUsAxRNveCGWzr+DJw4VhtpaKeIcJgztV3A6UUGMymoS2aju3HjBEYk0rgvpiJm66cUMKosqOb84oM1ipfhvF52DDdyqpaU9ZyUoYjtNMaqA5PGJtdUK4uCwZzfpHFTfnSvp23Pm4HeiuaMp31lTRhTawWX6PJaRjLSe7oA9+ql8V885SJEANMXK9amMKeS6XKyV4o6Z2uzJ51WHLQr1TR5dj69+VH3aYZBXPvtFnY+HKq925Rxy9g1llpJnYwCImcx2BnodcTDpHMSFz7bzEids1mgJFL93woFDMHQ/tSLKFVkivLo8CsrXXyFDBuLEwhSET56kJgJgEjmp6m5uz1nDDMeIgzZP8F/mHFQLhosYxrTQqG6J3JRpsf7BZyYXrk1r+jEHDALAQ0gH71SnIHSAu1hJXJ9sQRflb64jJrtaOvDq8s2F39pYWik8sqsciAOfN6L+z9V5D7CfglinIcSF+mE/W4WnmwINys/z97HpW1Tr/2fCWFGlGlO4ovcppAvWSiqxkrjzsuwEsmJe/YeWzRP7qERKSLcrIkLMf+DOi6crwwXLbLSpPBE6k0ogOsbQBj4ZmtvyVs5BrE1/pq9arNSqUHi7iZX092wOhqby3ZXHFLZD4Bx0Ig6hzH89/WBtmWI+K+cg3S413Ia8GH/91In6vAWFtLnMDMV3yWk2zKrFWhzXKo4zBsnow6SIz+eZtb9lhbT5uNekQvqlb59S4Jpdz1svFSVxDkS0iogWEtE8ImoI0wYT0bNEtDz8P0jIfy0RNRLRMiI6S0g/JiynkYhuoTQntaeEK+dg2wm5EhIxRMfHDh6GG89PipSAdArSTCZ+0FBtdRbDFdYU4gBMoxcxvY8z5xAW2dqesz8jVK83lQz+y8ThohPGK/LGC5EPGALS777EPjH2ZYyI2HfGMc5BQ0TFz3GlYdGNPWOZQq6cQ/I5dbrYJ2MG9sG/nXmItoyCn4OljUTSudSaT6ZSKqvKFjlYnQivUDeAPYJZqk2sZFdIB/87CgppfV6VP8b4of1i54l3BsohVvoEY2waY4ybo1wD4HnG2CQAz4e/QUSTAcwEMAXADAC/JSL+hW4FcBmASeHfjDK0SwlT+IyY05jjTgAwK3DFcq6ZcajWtjkNskQxE1G9KMaNOMh9sZ8mzAYQnVRlQ8F3osMuVjK1RS5P1Dm8+p1P4PJTDkrkFfPc9o/H4OlvJRXMRrGS4l5MrGSRF3OYeorfY4zhh48vBqDnHOL+Kvoyn1n0YeHaxP0BxXMO2u/jaPwARHPChXOYszI6elMnVrrunMmJNFUbxM2fjXMgijtwanUyFIRyt5l2R+M35BwsmzWROHw0DPuuCvtdSVRC53AugHvD63sBnCekP8AYa2WMrQTQCOBYIhoFoJ4xNosF5gH3Cc+UHZGtkkomGV27DnBbXlc9hstJVhxb9rTFdtB6eWh0nYZzUNnBc9Rk3eLI8yLbc6wQMMwFuh2dylppnMZb99jxgwvXR+0/EMP7J7kqk/xadU/83qZPJW4gzeKnZFk6ziFGcAyFflYI8S4HlUzWHy+nVGY967gREe+7cA4idAYiqvqUnIPQJy46hz1tIuegEysFhCFn8Y3inGd7h92UVQ4t3q+mKqyrZxEHBuDvRPQ2EV0Wpo1gjK0HgPA/N98YA2CN8GxTmDYmvJbTEyCiy4iogYgaNm1KnpLm1GCmd0KJ7bQtA1dkI806B5E4uLYyHSYMUx9PGDNltSjAROgC9AGwmuypyhxvOz5RGOt6c87gv1OwQ4eFyijyUWwcRJPHx+evc6rbReEq7hB1Czovxra5FyPk2jiHYlUMOoW0q9k04M45iLfPmjICP7/gSHXdKv2hSucgEN8R9WYOnohivhpaP4dQIf0/zy83lpdJK1YSCEDk9d6ziMNJjLGjAXwSwBVE9DFDXlV3MEN6MpGx2xlj0xlj04cNs8fPUZZhuCcOAJMFEuDOOYiD1LRYXPaxyKxuSL8arLrhU8b6ORb/51laMZC7ziH+2zRx3HUOUaE6K5w0kMVK3O9ABbGbdcTml59XLzRXn3UI9h+S5EhEQwITl+cqVopMNQ2ZCnmDzLZFV/SfsZkcF8spHKYI0AfIOrjycA7i/S8edwAGKXRHgPpdVASjRpiLNvEuIR4HSfeZ+CZx2952AMBt/3i0Ol/BWslBIS2JlXL5uPlrZ4XPKGnWMsbWhf83AvgzgGMBbAhFRQj/bwyzNwEYJzw+FsC6MH2sIr0yMOkcUnEO0X2zX4DbAj1j6qjCxLtgevLAGxUyBPTRyKmB9MpRjkNHqhcAwN3PIWYplYJl0g17UcENQHv8JhDvZ51sVxcL6woN0XFdeON9rs1W6HMX80w+Lm0Lep3QxrTWSi74n5nT8KkjRinvxfrcUnhkymquT7yvskpT1V14VpEm6r6sYreMHHXZHFCP4xOKsN5iG7kpq03nIEqpZM6hs86QLpo4EFE/IurPrwGcCWARgMcBXBJmuwTAY+H14wBmElEtEU1AoHieE4qedhHR8aGV0sXCM2WHyZRVHNSL1u5I3BdR5TgZRCsX26RZun5n0EbLt+fHhw7uV2MZZOZyVLCd1ubqIS22qzpFQ4Zodofy5DKJTdI4g7lCJMLuXJiJIAf/f/1Co7XujONiKrbRRsR5H335xPHOXKpp0xDj1ixcN99tp+EcTOJE1UZOlRYT8Vrqrq3KSEYf6vzy/NNZ+skRbk1zIpOJ+3cUjpEN29xZpqyl+DmMAPDncJGtAvAHxtjTRPQWgIeI6FIAHwC4EAAYY4uJ6CEASwB0ALiCMcZ7/3IA9wDoA+Cp8K8iMJmyit+raZvewxOQd0r6fOIi5uoo1Lhxt1PdpTrTqWAz8XTnHIQFOgXnMF1QJoso+DlYzj4I8kbX5TguE4iLlUx9EOco7QTsvlmrrXXzIm0xv8T67MQh+J+GgzAt+q66lqBtnBMy1yeORFOkZFK8qs1ayTYkZfGpSedgq1fMx6MUG7lPIqzcvKfwmxMDvn50ls6haOLAGFsBICG4ZYxtAXBa8gmAMXY9gOsV6Q0AphbbljTg3aoOpU1CPvMHeHbJhkSZKqTZrRTyOcpsh+yn3mVzFDOGbAPP1SxVfAWbctS1mRmKFgnXhcq1z23oWxNNFee6y8TVuTrWxYwkrKasbnJ/Ea6OlLYxwhdIm5hMDH/NPcpVUHIOiramEX3J76DL7zofshLnYPo+WSK8tyGKF9VR2BB1LnHodeEz8kZrJfdy5BDBOoi7ZoufTPSMI3EY3M/MORTjZm97pBjOwXX3/rt/+ojxfnC2gDmCabF12+AuVhJFWgbuRuJdP62R5QdlurQwnQiNZ7UtblefFTmzXmf6EwAAEXZJREFUmTjAuAGC2Uz0xWWBpeFGxzn0jdMmGS3eVN9DdY5ImrEgi0914811oeb9c+tL7wMwG3YQxfszJ1sr9QSFdE+EWazkPnhs1kwcomzRdRdrE8PwRcDkrAYUp7iyi5Xc3iFuMWR+H67sG26zIKFISegaSLBczvYxcZHDQTaA2+l7HLrzAgB3caRY30BD0EgRNuuzc8J4PoBFTi4SB8fdtKsuao/17O5kmsoHRhwzaTjkyaPqCycgynA5UhhIEiYboRX1HT1OId1TUehWJefgLrMVB9qYgepTogCJc3AVD1iy8cWixpKxmDFkIw7u4aPdZO+xsi3RVEWxkqvFUBq4PuYqLjKLYdzb+M4H253yxcxJLeOXL442ziGuaHbT8+g8vWW4zgddYEkOm/K5kE90ZLTUKfbLTz93OIZo9Huu3LnL6ZIc8vt2lUK61xEHzjqoFiJxjFkHpDDQTGx0VREKade4TjYCVoxY6f1Ne6x5LjhmLG6/6BhjnjQK6YL82zIaMyLn4Cj3T4NnvmVy04ngejZGGp2DqxWYuV3uZfD1xVav+A7Kc0hCiH3+RuNmbT4AuOFzhzu0MMIlivhZIlzFRSIHa5sa/QQdk4nj150ZLUM+PrTD4FEtc5GiQprIcw4VQyF8hoVzOEBx+LcIm6iEoyZmIeE6iG06h6BMW4TUNCE50uAXFx6JM6eMNOZx3UEDgte6hXMQD54x9X+xaoZJI/ob7z9+5UkA3ENHm9sYb6QpENyRDudVB/W5vzjfOKQ5a8O0CdpfEOPssoiB0lrZ6ZzfOFy5MHGPYts4jRakAaZ5O3qgXhwoYuqYuBmwafcvh5sRRWBZIs85VAqupqy/0Ljqczgv9I7nAMTbYRevAHbOweEo7IpB5H5sCxAf63YLm8iU1dXPwRX8cCITJg7fD0P3q8F3zlLLn2Wk0TmYvMi/fOIBJdcnI1+EWMmEkQOiRbLWMi4LmxqHsqeM1vtWpIUYudi2vA7sGy3Qpn694Bg3h9X+ddWx6LGqSMEcclBAHnIDCOaVV0hXCFFspeQHFxeVAX3NCj3uKfoNSwjlmLOcY2/bRFp8aNh0DsUMois+kYxyWgzSWK/wXZyNeGYcZerFSJVcuqpvTRUa/uMMfOxgt9AtLjoHHqrkS8frCYC7Kav7i/OxkeYgJlfYIp666qAWXHcm/nT5ianr16G+rhpfOm7/4Iflg9c46gqJKMY1mSAO2akGZ1O+meIWcgJtQFVGfUpcJdALD/sJoPrcaZynvn3mIfjKSROsMVpix346TrSrztDHwgeixdSmmzh4RHBO8CUnuO08AeBqx12xDeJuy0Yc+Dy19fn2MH4NYF54izFfjR2FWSYYdRPh/741VThsVK3V8swFrqJOAGjvYIX6TSiKOFginu5sDr6j7Z3TRPMVcehIvXiwcMiSpYw0p0K2O7LovJzjJqgdPTl4/w2vr8XqLXsx8yNR1KHOFCv1PuJgiK0UUyY6xFZyOZuBT0KXMjl0p8Vx8NDFu1vMst1DR9Zj3g/OwMC+ZpltJSDuvOxipaCP0ohFXP0cXOES7TUtXNrY1pG3cop8XOriGkVlBv9dCA0Xz5mUzEBxToQ/tSicTzhoKD46aSiu+8yU1GXb8PA/n4DpBoc5/jq23beoz7P15/odLQCAI8eaQ8/wb36gJooyh8g5rPjp2YmAoF4hXSGYznMQUS7nqcmCzLRYE0sZZ08dhXGD+ziJN7qCMADprLS4iCNNn5fL+5gjjamhK8wWVcH/1o68c6gLm4FBv9oqjKyvw88crIH4+9qIgyo0hQ1HWOJzDehTjfsvPU4Z/bZU1GQzlnDsAeycQ5CztirjfECXeJ6GCrzOkfV603fg/7d35rFWVGcA/31vYXks78mmIFDQWsIiiyCVCNqauEDU10RrtFZQ2lgTrdbGtoptTdqYqGlNRZtSYqHazabpIlat1aZqXShK3UCKKDUVQa0BWbRC0a9/zBkY311m7ntn7sy8+/2SmzvvvHPPPd98M/eb7zvnfOdgqO9DVTdD6aMPreY5pIRWG5GO4Ms4+GonyrhhA/jb10/y3q5Poj94h8SM34QqqSUsUp1k5/yIYQPY7HLYpOE5VAuLhPf7vv0fxGdQdfIkWb2+eknZzDUlhJ5D/1b/YaWkC0TTIGnakLhzGXq+I2L2fYhyRmTBYDnCxXwdMfdD6DmUW6hXbn/ptGg4zyEk7ppP4Tc9lrsvPZ6bzi6/x3TRCFcRt/Vpjo1rh/gypEmn8P7BTU2F4EfaN9WeyqOJBOMHhYN3n78JO934zbCY/FzdUYmvTLjdIW7Fd3jfJx1ziIswwMGdBystlAsJjcPg/tXvh3ARYTnj0NokiVdl95QG9ByC9ziV+0q7AMHGNG/uej+23rQxHYnntOed8AnOV2K3Wkjqdkef7H16DqMP6c+WHf+teg2F52Xv/g8T773g84lx+pgO1ry6PXbNQbc8hwyNQ1yIbs4RQ1n5+KtMjRkfqEWGuy4+LtHMwPBHPW7A/oDnUKbNvq3NqYRAy9F4xqHKfg5pcdWp1Wcf9UYOPHnVcJp9/ah0J2ulzxtu1WVzeXtP9aRyB2LfGp+fZ84RwxjV3o8vn3SUpx7C7RfO4s2d78eOByXd+S9KtdxTaRMXVjpl8mGJJmmM6uiPCFx5cvw5b2oSmhJ4GAdSlsSc09B4fFDmuujb0pTKzLpyNJ5xSOg5GD0j6VaQUertOUTZ53HF4JABfaoucoKPrnl5LCbdRHtbK09ck2wsISmD+7UmmiranUkUviZedIdBCWRKMkljyIA+bqaQP1mS5rMK752tO0ujDX1bmmL39fBFw405VEuf0ch8p3MyP41JmV0L4QU+roYZKb48h/Eu9UkSY7Pm2uBHd2+dbriQJLHsohHdBz0rBsdMC64F39GFpMYhXAdSjr4tzeY5pMWCKSOZcNig2Lhfo7EwJrlZrbT3b+VH5x/DsTELfqL4uhnb21oTb3051O2JMbuGfvogKmpcqvKs+fbpkxjZHp9DaMmCiSxZMLEOPSrlprOm8ucX36hruLhWwjQYccah2v+bmmD15u2oauqyNpxxGDu0LXZ+dZI8O0Y884+uvmgrZHC/FnbFLOiL4jNq0dwkPPTVEz6SaK0eRMNtj3zt03X97lpZPHd81l2I5Zxjx3BOZCVxHgkzBcRdvmdOG8Xbe/aWlWf15u0A3L/uDRYkvL+6S26Mg4icBtwCNAO3q+oNWfRj0/XzvW0taSTj3svnsX7rrsT1Pxfmx/HEx0dUz8aaBtFLLG5FvNE7mDq6nfVbd8VOo21qEr44r3qIbsd7+/x1rFI/Uv+GBIhIM/BDYD4wCThPRCZl0ZfW5qZMB9QakTFD2jhtSvUU4ADLPh/sIXHk8IFpdyl1wtlcE0f6yzpq5JvrzpjM0vNmMMPDdPV6PMDmwjgAs4GXVXWzqu4D7gI6M+6TkTNOnXwody6eHbv5SxEIF5/Fpa8weg/9Wps5c9qoHo0V+Paaq5GXsNLhwGuRv7cAn8yoL0ZOEZHE6bLzzpHDB3LVKZ+gc3r1fDyGEWXJgols37OvLmNkeTEO5UxpSWhORC4GLgYYO7Z+FtQwfCMiXOZxUZvRGAzs28KymC16fZGXsNIWIDo0PxrY2rWSqi5X1VmqOmv48N7xBGkYhpFH8mIcngKOEpHxItIHOBdYlXGfDMMwGpZchJVUdb+IXAY8QDCVdYWqrs+4W4ZhGA1LLowDgKreB9yXdT8MwzCM/ISVDMMwjBxhxsEwDMMowYyDYRiGUYIZB8MwDKMESbrfbt4Qkd3AxjL/Ggv8O0ET7cBOj/XSaDOpLGl8dxpyN6o8WV5DtXy33Tv1rVdLXZ+6maCq8dkmVbWQL+DpCuX/Sfj55T7rpdFmUllS+u405G5IeTK+hmr5brt38iu3N91U+u3s+uqNYaV3Eta7x3O9NNpMKksa352G3I0qT5bXUC3fbfdOfevVUjeNa60qRQ4rPa2qs5KWF5HeJAuYPHmnN8nTm2QBv/IkbavInsPyGsuLSG+SBUyevNOb5OlNsoBfeRK1VVjPwTAMw0iPInsOhmEYRkrk3jiIyAoReUtE1kXKponIkyLygojcIyKDXXkfEVnpyp8TkU9FPjPTlb8sIkulJ9sx9QCP8jwsIhtF5Fn3GpGBLGNE5K8iskFE1ovIFa58iIg8KCKb3Pshkc9c43SwUUROjZRnrh/P8hROPyIy1NXfIyK3dWkrU/14lqWIujlZRNY6HawVkZMibaWjm6RTrrJ6AScAxwDrImVPASe648XAd93xpcBKdzwCWAs0ub/XAHMINha6H5hfcHkeBmZlrJuRwDHueBDwEsEe4DcBV7vyq4Eb3fEk4DmgLzAeeAVozot+PMtTRP0MAOYClwC3dWkrU/14lqWIupkBjHLHU4DX09ZN7j0HVX0U2N6leALwqDt+EDjLHU8C/uI+9xbB9K9ZIjISGKyqT2pwNu8EPpN238vhQ546dDMRqrpNVf/hjncDGwi2fO0E7nDV7uDgue4E7lLVvar6L+BlYHZe9ONLnvr2ujK1yqOq76rqY8D70XbyoB9fsuSFbsjzjKqGG6CtB/qJSN80dZN741CBdcCZ7vizHNxF7jmgU0RaRGQ8MNP973CC3eZCtriyvFCrPCErnVv8rSzCMFFEZBzB083fgUNVdRsENwGB1wPl9wo/nBzqp4fyhBRNP5XIlX56KEtIkXVzFvCMqu4lRd0U1TgsBi4VkbUELtk+V76C4OQ8DfwAeALYT8I9qjOkVnkAzlfVo4F57nVBXXscQUQGAr8FvqKqu6pVLVOmVcozwYM8UEz9VGyiTFkm+vEgCxRYNyIyGbgR+FJYVKaaF90U0jio6j9V9RRVnQn8iiDWi6ruV9UrVXW6qnYCHcAmgh/Y0ZEmyu5RnRXdkAdVfd297wZ+SUbhDBFpJbi4f6Gqv3PFbzp3NwxJvOXKK+0Vnhv9eJKnqPqpRC7040mWwupGREYDvwcWquorrjg13RTSOISzC0SkCfgmsMz93SYiA9zxycB+VX3RuWe7ReQ450IuBO7Opvel1CqPCzMNc+WtwOkEoal691uAnwAbVPXmyL9WAYvc8SIOnutVwLkuVjoeOApYkxf9+JKnwPopSx7040uWoupGRDqAe4FrVPXxsHKquvExqp3mi+BJehvwPwIr+QXgCoLR/ZeAGzi4mG8cQabWDcBDwMci7cwiuAheAW4LP1NEeQhmYqwFnicYnLoFN0umzrLMJXBhnweeda8FwFCCgfRN7n1I5DPXOh1sJDKrIg/68SVPwfXzKsGEiT3u+pyUB/34kqWouiF4aHw3UvdZYESaurEV0oZhGEYJhQwrGYZhGOlixsEwDMMowYyDYRiGUYIZB8MwDKMEMw6GYRhGCWYcDCMFROQSEVlYQ/1xEsnUaxhZ05J1BwyjtyEiLaq6LOt+GEZPMONgGGVwydD+RJAMbQbBAsWFwETgZmAg8DZwoapuE5GHCXJfHQ+sEpFBwB5V/Z6ITCdY9d5GsFBpsaruEJGZBPmz3gMeq590hhGPhZUMozITgOWqOhXYRbC/xq3A2RrkwVoBXB+p36GqJ6rq97u0cyfwDdfOC8B1rnwlcLmqzklTCMPoDuY5GEZlXtODeWx+Diwh2GjlQZfluZkgFUrIr7s2ICLtBEbjEVd0B/CbMuU/A+b7F8EwuocZB8OoTNfcMruB9VWe9N+toW0p075h5AYLKxlGZcaKSGgIzgNWA8PDMhFpdfn1K6KqO4EdIjLPFV0APKKq7wA7RWSuKz/ff/cNo/uY52AYldkALBKRHxNkybwVeABY6sJCLQSbMK2PaWcRsExE2oDNwEWu/CJghYi859o1jNxgWVkNowxuttIfVXVKxl0xjEywsJJhGIZRgnkOhmEYRgnmORiGYRglmHEwDMMwSjDjYBiGYZRgxsEwDMMowYyDYRiGUYIZB8MwDKOE/wOypFk5Txw8PQAAAABJRU5ErkJggg==\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": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXe4JGd95/t9K3U+fXKYHCWhUdZYyCAULBsJvAvYjzFiuUa2wVpjcX291+GC2We917vsru991mDsxdeywQRjgglG2CZKGEkgJEZxRmGCZkYzZ04OfTpXV3jvH1VvdXV3dTgdq7vfz/PMM2fqdJ+pU131ft9fJpRScDgcDofjRuj1CXA4HA7Hf3Bx4HA4HE4FXBw4HA6HUwEXBw6Hw+FUwMWBw+FwOBVwceBwOBxOBVwcOBwOh1MBFwcOh8PhVMDFgcPhcDgVSL0+gWaZnJyk+/bt6/VpcDgcTl/x1FNPrVFKp+q9rm/FYd++fTh27FivT4PD4XD6CkLIq428jruVOBwOh1MBFwcOh8PhVMDFgcPhcDgVcHHgcDgcTgVcHDgcDodTARcHDofD4VTAxYHD4XA4FXBxGEB+9Moazqyke30aHA6nj+HiMID8/j88j794+HSvT4PD4fQxXBwGkGReQyKn9fo0OENOXjOQKxi9Pg1Ok3BxGDAopcioOlJ5vdenwhlyfvcfnsN/+OKzvT4NTpP0bW8ljjc5zYBJgVSeWw6c3rKYyEE3aa9Pg9Mk3HIYMNK2xbAdy+HkUgo/OLXaqVPiDCmqbnK3Uh/DxWHASKvbF4ePPXwa/9eXn+/UKXGGFFU3kdO4OPQrXBwGjIxqPYxpVYfRoEm/ksxjPaOCUu4C4LQPVTeQ5+LQt3BxGDBSajHWkG7QelhNqdAM6lgdHE47UDXuVupnuDgMGMxyAKyU1mqcuLSFx19ZB2CJAwBsZAqdPTnOUKHqJvK6yS3SPoWLw4CRdlkOteIOH/3eaXzoa8eRUXVk7N3dOhcHThtRdQOGSaEZXBz6ES4OA0baZTnUSmfNawbmN3NYTuadYxtpLg6c9kAphaqbAMCD0n0KF4cBwx1nqGU5FHQTBcPECwtJ59hGlosDpz1oBgXzJvGgdH9SVxwIIbsJId8nhLxECHmBEPJ/2MfHCSHfJYSctv8ec73ng4SQM4SQk4SQu1zHbySEHLe/9zFCCLGPBwghX7SPP0EI2df+X3U4yLiCyu7gdDmqYe3qnr6w6RzjMQdOu1D1oiDwoHR/0ojloAP4XUrpawDcDOB+QsiVAD4A4CFK6WEAD9n/hv29ewAcAXA3gI8TQkT7Z/0lgPsAHLb/3G0ffw+ATUrpIQAfAfAnbfjdhhJ3xlEyV91y0HQmDgnnGBcHTrtgLiWAu5X6lbriQCldpJQ+bX+dAvASgJ0A3grg0/bLPg3gbfbXbwXwBUqpSik9B+AMgJsIIXMARiilj1MrfeEzZe9hP+vLAO5kVgVne6RVHZPRAIDaMYeCbTm8uLAFSSCYHQlincccOG2Ci0P/s62Yg+3uuR7AEwBmKKWLgCUgAKbtl+0EcNH1tnn72E776/LjJe+hlOoAtgBMbOfcGuUHp1bxga88D90w67+4D8moOsbCMgKSUDfmAFi+4cloABNRBZs85sBpE+44A4859CcNiwMhJArgKwB+h1KarPVSj2O0xvFa7yk/h/sIIccIIcdWV5vrBXR6OYUv/OTiwO5m0qqOaFBCLCgj2YA4AMBULIDxiMJTWTltQ9WK9xcXh/6kIXEghMiwhOFzlNKv2oeXbVcR7L9X7OPzAHa73r4LwIJ9fJfH8ZL3EEIkAHEAG+XnQSl9gFJ6lFJ6dGpqqpFTryAoW+GPgRaHgISRoNSQWwkApmMBTEQUbGTUbpwiZwgoDUgPppU+6DSSrUQAfALAS5TSP3V960EA99pf3wvg667j99gZSPthBZ6ftF1PKULIzfbPfHfZe9jP+iUAD9MOlVWGFUsc8gN6w6bzljjEglJDbiXAshzGIgqvc+C0DR5z6H8amefwegC/AuA4IYRN7vhDAP8DwJcIIe8BcAHA2wGAUvoCIeRLAF6Elel0P6WU3R3vA/ApACEA37T/AJb4fJYQcgaWxXBPi79XVUK25ZDVBrOPUEbVEQkwt1Jty2E0LCOR1TAdC0CRBGQKVqM0Zl1xOM3CxaH/qSsOlNLH4B0TAIA7q7znwwA+7HH8GICrPI7nYYtLpwnalsOg5l6nmFspJGHJVf3shlKKgm7iqh0jePpCAlOxAETBMiI3swXMxUPdPGXOAKK6A9ID+qwNOkNXIR0a4JgDGxEaDUiIBeSqMQfW6+ba3aPYORrCtbtHMR6RAYCns3LaQjOWg2nSklgFp7cMrTgMYgYFGxFqZStVjzmwYPRcPIgffuBncM2uUYxHrNoIXgjHaQfNiMPnnryA2/6ff+3QGXG2y/CJg+NWGryANKuOZjGHbMHwrOdgwWhFLH784xEFAHitA6ctuC2ARjdir6yksZTMD2wNUr8xfOIwwG4lNsshZmcrAfAc4OOIg1QMPMdDllspka0exOZwGoXVOSiS0LA4JOyNidvq4PSO4RMHZXDFgXVkjbjEwau/EhMHWSzmGUQD1uuzPHjIaQNsgR8NyQ0nfyRy1sbEnWZ9ejmFTe7q7AnDJw7McigMXiorsxKsbCXLEvBKZ2UxB0UqfvxBWQAhQHYAr0uj3PvJJ/HR753q9WkMBMytFA/JDW/EmNXK7s+8ZuAXPv4j/PnDZzpzkpyaNFLnMFA4FdIDHHOIBiSwtoVbOQ9xsHdmAZc4EEIQUaShniP9wsIWZHHo9ksdQdVNyCJBOCAhpzX2rCXte5W5pH70yhrSql4ykIrTPYZOHESBQJGEgXQrsVkO0aAEWbLUwVMcPCwHAIgERGTVwbsujZJRDWeB4rSGqpkISCJCstBwnYPjVjKs13/3xWUAPEmiVwzlNikkiwOXyrqV1fDAI2cRkkVMxQJOgLmW5aCIpZXQEUVCZkjdSoZJkdMMJHJ8IWoHqm4gIAkIyWJDGzHTpCUBadOk+N5LVru2TZ4k0ROGznIALHEYtArp93/+aZxZSeOv7z2KaECC0IBbqdxyCAfEoQ1IM1H0ul4PPrcAWSB409Vz3T6tvkXVTUsclMbEIV3QYdrd1Aq6iefmE1hNqYgFJEc0ON1lKC2HcIM3bD/x7IUE3vFTu3HbZVa32pAsQhZJFbeS9btXiMMQxxyYS648lTet6vjQV4/j//vBK704rb5F1U0EZRHBBq30Ldd1V3UTT71qja/9uSMz3K3UI4ZSHILy4O2QNdN0Os4CVoA5HpI96xYKurVFc6eyAlYge1izlViNiKqbJYvZ1565hJSqc9fGNlE1A4rtVmpIHFybmIJuOp/H/okI8po5cG7gfmAoxSGkDF7MQTMopLLFfiQkewZYWUA6UGE5DG9AOuOymNhCRSnFZ350HgAPim4XVTcRsC2HRly47k1MQTeh6gYkgWDCHnnLr3/3GU5xaDBI1i9QSmGYFJJQ+nHGQzIPSDeIlzg8N7+F0ytp7JsII5XXofG2Dg1THpCuN57FnQhQMEzkNcstNRa2Eis2M9xy6zZDKQ6N7mb6BdZltTyGMFpPHLwC0kNqObhjLWwXu5jIAQBef2iy5DinPu6AtElLJw96kSiJORjI6waCsoDRsGJ/n1sO3WYoxWHQ3Eq6aT14klDqVqpuOXgHpJnl0KEhfL7GHYNi14x1td0zHgbAXRvbgdU5sKLTepMXy2MO7P1jdit5HvPpPkOayjpYRXCaHWCWxNpupfNrGQC1iuAkmBTIa6bTg2pYKLUcLBFIqWXiwHv8NIyqGwjIQkmjyzjkite9sLCFFy4lK8Qhb79/LMy7BfeKoRSHsCINVLaSZlY20gMscUjmNZgmhSAQfOgfj4OA4OYD456vjwSsBzlT0IdOHLxiDmxY0m5uOWybolvJ2oBUs9T/9ofn8dWn5/Gmq+YQUURkCgZU23IISiJGw6xbML/23WYo3UrBAQtI6wZLTS39OEdCMigt7oDXUgWsZwoosBiFWFnnAGAo4w4ZVQchACHFHj+pvI6IImIiynav3LXRKJY4iAhKtbsgLyRyMCnwyOlVTI8EnfcyyyMgiQgrIr/2PWAoLYeQLKKgmzBMClGoNh67f2BZNF4xB8Ba7JiLSZYICroJRRRASJnloBQth2EjUzAQUSSIAnF6/KTyGqJByXFt8Cl59fmbR88iHpKhala2UrBOi/xLdtA/lddxxWwM59YylltJMxxhGQsr3GrrAUNpOdQzdfsNJg7llkP5AJ9kXkM6r1viIFV+9BF7pkNmCKukM6qOSEAsidOkVR2xoIygLCIki9y1AeBbJxbx0//9oZKZC26+8JOL+LsnLth1DsWYg1fzPdOkWEwUO66OhRUoooCCYdoV1tY9Wq2Yk9NZhlMcBmwanG6ygLS35bCV06AZJrIFA2lVR8EwqogDsxwG47psh7SqI6JIGA0XF6JUXneGJo1HFGzwXHu8sprB4lbeMwsOsET23GracSuNBKtnG62lVRQM07FYR8MyFEmAaldEB5jlEJG55dADhlIcijMdBmMRrGY5sBzxrZzmPMyaQZHK6xXxBsAdcxhWy0EqsRySed2ZkGeJxmAtUM2kLLN7rZp1mc7rSNopwAFJwO7xEADg4ma24rXztkvprqtmAVj3qyIJKBiGXQQnOMe55dB9hlIcBm1UaDEgXd1ycLfR2MgUnHkPbiK2OAyj5ZBRjUq3Ul5zdr7jEQUbAyQO85tZXPVH38ZzFxPbep8jDh5xKUop0q7jAUlALChjLCzjwkalOCzY4vC263aCEGAiYruV7IA028SNhbnl0AuGUhxYg7pBsxy82mcApZYDYImDp+Vgu5WGsflepmBZCW5xcLuVRsOKU+eQUXX89uefwUWPBa9fuLSZQ6Zg4FsvLG3rfawa3ysVPFsw4DZGAvbivmc8XHKtKKUo6CYubVricP2eUXz+N27GPTftQUAW7IC06fT+Ggsr2MppMMzhK87sJUMpDsEBizmwB7Y85hCUBSiiYFkO+eKCv5EpQJEq6xiYC2UY23ZnVB1hO+awldNAqeV+Y9fE2r1aovH9kyt48LkFPHluw/Nn5TXDycLxK6odUH7s9Nq23scC0V73SPkxtrjvHg+XWA5/8q2TeMtfPIZLiRxiQQmxoIybD0wgHpKhiILTGZc9p6NhBZSCT+nrMkMpDoMXkPaOORBCMGLvhN2Ww3qm4BmQDkgCBGLVOXzpJxcx7+EnHlTSquHEHAyTYiunIacZiNluJbZ71Q0Tj5xaBVA95fczj5/Hm//sUV+3IWHicGJha1uV38xK9aqFqSYOe8bDuLSZc3b+Z1fTeHkphW+dWMLO0VDJexRJsOscipYDs96GcdPSS4ZTHJTq6XX9SLUiOACIhyRs5Qolu66CbiLg8VpCCCKKhFdW0/iDrzyPX/3bnzhVwoNORtURDYgYDVlB/Hnb5cEWJtYdNJHT8ANbHKotVuvpgiUkPnaDqHZ/LUqBx8+uN/y+WjGHtG2dstgXyzbaMx6GblIsbhVrGgBgJaVi11ilOLDrytxSTCTqNe/jtJfhFIcBsxyqFcEBwFQsgJWkWpF66GU5AFbcgblLzqyk8YGvHG/z2foPNj86EpCcauizdh+qKBOHiHX8yXMbWE6qAKpn7LBFTK1SC+AHVK14bo+dady1xFyYXr87O3b5bAwAEJCLlgMAx7WUUov34o4yyyEgCc5GJuiIg1hxzpzOM9TiMCj9lbQalsOO0RAWt/JIllkA5ZlNjEhAwnqmAEKAXz66C/98fBGmj3fA7YAF4COKhJ32TvblxSQAYMSxHCxxeOCRswAsIc5UaTPCxFr18eaDCdf1e0bxzeOLJWM6a8FiDl7PDmvTcvXOUQClMQcATlA6ldcxYYttpVtJdCwL9n4mMsza4XSH4RQH5lby8cO7HZyW3R4L/s7REJaSeWxmCo6LBKhuObB01gOTEeyfjALw9w64HbBFPhKQnMXqJVscogHLnXT5bAw74kE8ezGBIztGMBULVHUrsS65fr5ubKH9/bsux1ZOw59+92RD7yvUqHNgbqVbD09CFolzLefiQUgCKVoOeR0/d+UM7r/jIN589VzJz1BEL8uBiYN/r+cgMpTiMHhFcN6N9ABgLh6CYVKcXkljOhZwXE9e2UpAMc33yI44QvaObVDcb9Vgi3wkICIWlDEalvHSYgpAMeYwMxLEDz/wM/jRB34Gn3vvaxEJSCUL5D8+cwnv//unAbgsBx8vZnnbRXP97jG867V78dkfv4ozK6m676tVBMfiEDftH8eLf3w39k5EAFit5HeOhXBhI2dngWkYDSv4/buucKwKRkASHAuEFcE5biUfX89BZCjFQRYFyCIZmEVPN6pbDjtGrU6XJ5dSiIdkx4fuJSRAsb/SkR0jA1csWA3mVmJpq7vGLGsLQIm1RQjBjtEQRsMKIgGpxHJ44twG/vWkFaguxhz8e91U18CnX339PpgUeO7iVt33FQPSHm6lPBNZqcLFuXvMSmdVdROaQUuuqxv3XHMmCo7lMOD3od8Yyq6sgGU9DE7MwbsIDij6dLMFAyMhGdGAhERWq+5WshfIq3bGsZa2Aq+DYmF58ZuffQoUluXF2ofsGg3jxCXbrVRlEYsGxJLds6objj+efR55HwdQVd2ELBKIAkHM/swb2QQwl5lXoWRa1SGLpGSBZ0zHAji3lnFiXyNVrqv7vmSWQ1DmbqVeMLTiMBkNYNVe/PodrUr7DACYcwX84rY4APB8gIFi2+4jO0bwhJ21NCixmXIopXj45RVnp++2HBisfUY5EUXCWqpYH6DqJgqGCUqp83n4eafLxnACrnYyDWwC2LVKewTjWX+q8lbwAGtcWHCsi1iV61oqDmXZSlwcuspQupUAYGYkgJVkvv4L+4BqRXCAteCxXdpIUHbM+WqWw2sPjOPNV89iNKwU2y37eJFrBEop/vGZSxWfd14zS3LnWVdaJg7VdsGAdV3dbiUmBJpB+yLmoOqG87ttJ3uvWATnHZBmAlvOeFRBTjOwYqcBV3Mrud2dTraSxLOVesEQi0PQyVfvd6q1z2CwXHK35VAtlfUXrt+Fj7/rRgCD06Dw+fkt/M4Xn8U/PDVfcpzVfly7K46IImIyFgAA7BqzgqTRKrtgwHK/uQvBmBAUDNNxL/lbHIoVyJJotVnJavUrkGvFHNJqDXGwU4EvbFj1I9UsB5a2CvA6h15TVxwIIZ8khKwQQk64jv1nQsglQsiz9p83u773QULIGULISULIXa7jNxJCjtvf+xixnzpCSIAQ8kX7+BOEkH3t/RW9mRkJYimZ93WLg0ZxWnZ7xByAYtxhJCQhaj+Uilh/RnSoSlbXj8+u4+c/9mjfWBSfe+JVAMBqqnQzwPzfv3HrATzzn97ouJB22W2mqy1gACqyldi1KOimy3Iw8KMza3jnAz92jvkFa5hO8R4IKWJDbqVaRXA1xcGua3h13UpnrW45FM+JTYIL8JhDT2jEcvgUgLs9jn+EUnqd/edfAIAQciWAewAcsd/zcUII+7T/EsB9AA7bf9jPfA+ATUrpIQAfAfAnTf4u22JmJIiCblYdWtJP6AaFQAChysjTOTtjyW05VHMruQlWSWX9wpMX8MJC0glY+5mtnIZvPLcIoHLMJ/vs4yG55HowMa22gAFWQFozqOPqUF3B6GLMwcTTFzbx+Nl1pz21X1C10oFP4QbFoVgEV0UcqlwzVnleVxzc2Ur2/cdcTdyt1F3qrhCU0kcAeLefrOStAL5AKVUppecAnAFwEyFkDsAIpfRxam3VPwPgba73fNr++ssA7iTVbPk2MjNiuRAGwbWkmaZnvIHB3EqNxBzcsJ2l25zXDRPft1M2WU8nP/P1Zy8hpxmIh+RKccgWxcENq3WotgsGilldrAGdt+VgOn581qvJL1hjPEsth2wDlmCxCM7brRSpcs1Yhfn59dpupZKAtG05CAJxurVyukcrMYf3E0Ket91OY/axnQAuul4zbx/baX9dfrzkPZRSHcAWgIkWzqshZkas3fTSAASldYPWFIedHjGHRsTBqwfVU69uOjtuFgj3Mw8+u4ArZmM4uncM61UsB6+MpKt2xLHPLuLygi2CLChdEnNwuZXYtbvkO3EwSoLtjVoOTPhymlExXyGd15202HImItZm7IJtOVQT3oCH5cCO85hDd2lWHP4SwEEA1wFYBPA/7eNeO35a43it91RACLmPEHKMEHJsdXV1e2dcxqwtDssDIA6aYVYNRgPAzQcmcPvlU7hyx0gxlbWGmDC8AtIPv7zifF3Q/W05rCTzeOrCJt501ZydSllqJbrdSuX8zb1H8V/edlXVn82uIwtKs4Wr3HLIOZaDv9qfuwPSgLURaGTIk+bavZe/vlbMYSQkQRQIUvZrxCou0BK3UplQcLdSd2lKHCily5RSg1JqAvhrADfZ35oHsNv10l0AFuzjuzyOl7yHECIBiKOKG4tS+gCl9Cil9OjU1FQzp+4wZWem9GM661paxY9eKXbS1AzqWQDHmBkJ4lO/dhNGw4rjE/YaE1oOM+vdO8qHX15xfMB+txy+/eIyKAXedPUsxqNWnr07AcGxHDzEISiLNa0rZjmwwGxed7mV9GLMwbduJVedAwCEFKnhgDQTU3fqq2FSZAtGVbcSIcRxLdWK5bjTV93e5YAkcrdSl2lKHOwYAuMXALBMpgcB3GNnIO2HFXh+klK6CCBFCLnZjie8G8DXXe+51/76lwA8TLuQQhSURYyG5b50K33m8Vdx7yefdLql6oYJpYbl4IaZ/Y1kKwkCgSIJJVlJFzay2D9puVs0n8ccvnViEQcmIzg8HcVERIFmUKdvD2BlK8Vq7GJrEbVrIlgxGLMcrIB0pVvJd+KgGyVum7As1k1ZppSiYJjObAt3nQezoGot/KwTa63XsI1HeX1JQOIxh25Tt0KaEPJ5ALcDmCSEzAP4IwC3E0Kug+X+OQ/g3wMApfQFQsiXALwIQAdwP6WU3XHvg5X5FALwTfsPAHwCwGcJIWdgWQz3tOMXa4TZPq11SOY0a6HL64iHZegmhdSAmwgotoNoJOYAWO4GtmiougFVN520RL+lZ7pJqzp+fHYD9916AIQQjNs+7410wYkxbOU0T6uhEdyWA6W0xHIo9KFbKazUbyfDhhfFwwqwni2ZBpdRi32VqjHuiEP1ax5wWmaUbl4USfB1xfkgUlccKKXv9Dj8iRqv/zCAD3scPwagwolLKc0DeHu98+gE0yPBvnQrMV9vIldAPCyjUCfm4GbveASKKGD3eKj+i2GJA7McWOuDcTst0c/ZSolsAYZJsd8OKrNd60a2gH2wjiVbEQelGJDWDApm66ouyyGvGc5ntZTMQzNqZ5V1E0sctlfnwH6vUfuauYsAWbvuWhle4w1ZDtY5lYtDQOZupW7jjzu1R8zEAn3pVmLVqU7WkGFWLYArZ89EGC//l7txxexIQ68PKSJytsuEicNEH1gOThtze3fMFqaNdDFjaSunIR5qrr1Y1GU5uAOlmu6qc9BN59qZFFja8s+9pmqV2Ur1LAcWS2FuJXchHHPXVatzABqzHBSplluJWw7dZKjFYTYexGpKrUjJ8zu5CnGgDQWYGdWK5bwIysUdJZsn3Q9uJVasxXbqjjhkysWhdbeSu/tqXjed+8lyK+mYjFourYs+ci1ZdQ6l2Uo5zajZMYC5y0btwLK7hQYTipYtB8nbrcRjDt1nqMVhLh6CSYvjC/sF9iAycdDM2tlKrRCSiwHpZK7UctB9LKpMuNhiwyp019skDopk9SNKq0bJjjZb1owvpxk4PG1N1PNLUJpS6uFWshbsWm3GmTiwa7aeVvGDU1ZKeapNbqWAIw7lloPI6xy6zFCLw+sOWrV2D7ly9/sBFiBO2BW+mt2bvxMEXQHpouVg7YT9bDmojuVgXZewIiEoCyW1Dsmc3rQ4AFYX13LLoaRTq10hfWAqAkL8Iw5skS93KwHebTEYrMaBuZU++r3TuPeTT2J+M+vUC03bKeJeMHGo1gYdcLuVymMO3K3UbYZaHPZNRnD5TAzfeWGp16eyLcotB900O2g5FAPSyQq3Uv9YDoBVpcssh4JuIqcZNReqerDme+5FK102AChXsFp3TEUDWNryhzgw4SwpglPqt+3WytxK7P67uJHD4lYeiiQ494YX20llrbQcuFup2wy1OADAXUdm8JPzG1jvgyZyDBYDYIPYNYNCbjA1dbsEFbflYLuVnGwl/z6sLObgng/ABs4ArurocPPiwGY6uC2HbIkf3oBuUoRkEZPRANbSBa8f03WYe8bdWyncQHt2ZnGElNICwUuJHBYSOczFg1VbnAPFwlMmLl44RXAVMQeerdRthl4c3nhkFibtL9dSRbaSaUJuopCrEUKyiLxLjAgp+pz97Fbyshw8xaElt5I106Ga5bCZtf6vkD0rwi9dbNn5ermVaqWzOhlgooCdoyG85dodAICFhGU5zMWDNf/fQ9NR/MW/ux5vvHKm6mvK50YXjw9OnUO2oOOhl5Z7fRp1GXpxOLJjBJNRBU+ea7TxbO9hD3Ax5kAbrnPYLu4iuKQ96Ys9uH52K5VnKwGWW2M9XSoOzdY5AJY4WAHpoki60zvZ5xNWJExGFayl/CIOlW6lYI1pcHnNwKvrmeLcEFHAP/7W6/E/f/laTEYDljgkctgRr107QwjBv7lmR0UmkpthyFb6p+cW8Z5PH/N9X7ehFwdCCGJBuW9uPM3V9bOYrWQ2XCG9XYKyUOJWGgnKzv/l595KhSqWw2pKxeefvIDVlPVgtmI5RO2AtHtHy1pZRxTRsSJCioAp263kh+FSjltJcruVrDhAzmMa3P/6/hn8/Mcec94niwTxsAxZFLBzNIiLm1ksp1RnbkgrOOJQHpC2xcEP169VEjlrg8LctH6luQqgAUMSiK/9527cO7uSOodOupU066FM5jXEgpKTAdQPloM75vDma+bwr6dW8cGvHi+Zq90s0YCEdF73tByiQclx/4VkCZPRAAqGiWS+tQypdsDcSu6gb7hGQPqxM2tIq7qTkOCOb+0YDeGx02swTIq6QbgeAAAgAElEQVS5OpZDI4gCQSwoORlRDBaDKBhmRSZTv8Gqyf0+SXHoLQfAuiH9vNC5yXmKQ+faMgTtRUPVTaTyVrsJVo3t75iD9Xm6r8sNe8bw3f9wKz76juucXX0rC/VIUMZWTit5yFkqqDvf34o5WEFYP8Qdim4lV51DFbdSrmDgxKUtAMUYiltwd46GnOroHW2wHADga7/1OvzaLftLjjEXWL9Y+LVIlw2I8ivccoC1gBg+dpG4Yf1sJqOKIw4Fo/HGe9vFPUc6mdOxYzQIQSAQiL/FoWDvjssbDBJC8Lbrd4IQ4OvPLtRMvaxHPCQjpxmOe0AgxYB01GWRhBURArEyddbTBRxsrdt8yzji4GE5lC9Yz15MOELLYihuwWVTBgG0xXIAgEPTsYpjjjhoJtAeDeoZTpt3nxf1cXEAIInE19W+bpjlMBsP4sSlJHTDtLKVOhiQBqwUx5SqIRa0HlxZFHzdeK9oOXhfl7detxNvvW6n5/cahaXBrtqB5mhAcmIO7oloIVl0LAlfWA5aZbZStTqHY+eLiRqbdqZXybztsaIg1AtItwKzcgahEC6tcrdS3yALgq93wW7YroPt0pJ5HXqdYT+tEHSJQzKnO756WRR87YrzCki3G+aSWrHFIRaUHcuuwq1k91fyhTh4uJVYALhcHJ50icOG7VZyCy4bQRtWRIw02cSwEZiVMxhupdIBUX6FiwOsmIOfd8Fu2BD4HXZO+VZOs1pBb6Px3nYIutxKqbzmdNSUROJrQXVSWTskmkAxDXY5mYcoEIQVsSQgzQjJIsYjCgQCX6SzeqWyCgKx0pZd7TMMk+KZCwlcsysOANiy3UqKh1upXgFcq5S4lfqcfnErcXFAf7mV2ICVOfuh3Mpp0E3asUWQuRvWMwWYFM7uUBYF36eyyiLZVgfa7eK2HAKSAFkUwG4jt+UQVkSIArFSaXtYJU0pxenlVLEIrqxFRXnb7vPrGaRVHa87OAmgGJB2xxzGwjJCslgSe+gEg+hWqjd5r9dwcYD/Fzo3LBuGVaNuZqyhNp0sggOKs7aZ5SD7PMPLakbY2ds77rIcApJQ4sJy9w9iAmu10Oid5fD0hU383EcewU/sgs/ylNCQUjoq9NRSCgBw494xAK6AtOv3JITgDYcncfOBiY6e+2BlK1nPsN8rvnlAGqzOwb8LnRu2s2M7NdZIrlMLoSMOtjtkxHEr+TtOUzDMjsYbgKI4pPI65uLBkv+PWQ4CKbphei0OqynrXnn6QgJAZYuKcNk0uJeXUiAEuHa35Vba9Ig5AMAD7z7asXNmDFLMIcMD0v2D3/3nbpg4zI5YlgNrGNi5lt3WLVK0HJhbyd+C2o2RnO4aiYA934HBYg5hRXJ88ZNRpafiwFwyF+z5JeXiEJJL3UqnllPYNxHBhN2inaVOdzKOUw3HreTzBbURMk6dg7/XHC4OACRB6J+YQ0EHIcD0iPXAsjTKTmcrOZaDvSDKPrccVN0sWaw7gSwKTn1AUBbL3Eqyc5wxGQ1gLdW7mIPbKlAkoSKAXD5H+uRSCpfPxCDawWqTWlZ2J+M41RgUt5KqG04mHY859AGSz3fBbrIFA2FZRECycudXO2w5MH95MV3Tncrq3wdVM2jH3UpA0XqwAtLFz4DVOTDxAIDJWAA5zShpztdN3G6McqsBsKycrFZ0eZxfz+CyWauuhY1F7cY19aIYkPbvPdcIzGoAuFupL5CF/gpIh+0HdTyiOJ0dO10hvZCwBtWMuFJZ/WxtaV2wHAC3OIhQXAFe5lYKuSyHcXuOAfPdd5u8a2H16k/kthzOrKRhUuAKWxyiAev1nXbVVaMYc/D3gloP98aAu5X6ALHfLAd7NzoWUbCSZJZDZ91Ki1t57BwNYTTMspX8bTkUOlj74Ya52QJyMeZASNFiCLksh3Cg/syETuL+fz0tB1fM4aSdqXTZTKnl0DNxsM+30OeWg3veh9+L4Hi2ElhaZn/cdBnVcNorT0QUnF62HuJOuZVEgeB/u3kP4iEZv/GGA87iIEvE1zsfzeiF5WB9BrIoOKLqthwi9ueW6ZE4uBej8hoHoDj2FABOraSgiAL2TYSd7wGA0qH7rB6D4lYqGSPrc7cSFwdYLhnDxy4SNzlNd3al4xHF2el1KiANAP/1bVdXHJMEAbrh3370ahfqHICiOARdloMiCs5O1x1zcNpi9yrmUDCcuQhebqVIQESmYIBSis1MAWOR4uwOlprbqXG09VAGpEKaiYMoEN8HpLk4wE5l7RNxyBYM50GdcHUU7VQRXDVkkaDgY1ecZpglVcqdotRysK0qkTiLb9AlDmz33TPLQTMxGpYhEOLpVooEJBgmhaqbyKhGyfXrtVtJFAhkkQxMzGE8ovja8ga4OACwA9J94lbKqgam7UHt7nbT3XChuLG6svrnmn3/5RVs5TS87Xqr02pBN6GEu2s5OC43UXDcNmGXW6nY+bRHloNuICiLuHbXqOf3mRhkVB1pVS8Rh14HpAFLgPverZRnLfcDvs9W4uIAa1diUsA0aU9yuLdDVtOdmMN4Dy0HSfRPbUiuYOB3/+E5qJqBu6+aRVAWrZhDV1NZ3ZaDt1vJiTmovQtIh2QRH3nHdfC6zcOu88uoumMtAMVz71XMAWCjQv29oNaDuZUmo4qTAehXeLYSisFcvyx2tciqxWyliahLHLpctSoLpCeZI2tpFYmyVNAv/uQCNjIFZAoG/vXkKgDLcuhmzCHgshwUu1paIEBIcTXgC/TacjARkK0mgF4dVJl1kLYth4iHW6lXdQ6ALQ4+d8XUg20MJvrArcTFAcUagX6odXCnso7bbQ2AzmUrVaNXzQrf93dP4Q+/dtz5d0E38cAjZ3HDnlGMRxT88/FFAN0vggtKomMtyKK1+P73X7wabz+6y3ltuMoozm6RLxgIeWQpMZjlkC1Y4hArcSv1NuYAANMjQbyymu7Z/98O0qqGoCwgEpB871bi4gCrJQAAX3cZBSy3V04rTWVldPuh7VVV+UIijxcXks6/TyxsYWErj1+/ZT/uvmoWD720jFzB6Fq2Ukmdg8utBADv+Kk9ODgVdV4r2e6mTI9jDtVg1kFa1SvdSj4Qh1sPT+LZiwlnrkQ/klYNRAMygrLIxaEfYOLgpwCrF2l7UYkEiqmsjO5nKwlOj5huspXTcHEz57i0knYzuB2jIfz81XPIFgw8fnYNmmF6ZuS0m6LlUBqQrkYkIDkzObpNrmA4E9+8iASKlk1GNcrEofcB6dsun4JJgcfOrPXsHFolo+qIBkQEZaGkYt2PcHFA0a3k91qHxYTVKoONCA0rbldGt7OVum856IaJtKrDMCkublqdRVmALxaQsG8yAsBqRmjFHDovmHPxIHaPh3DZTKykzqEaIVnsqeXgrtguhwWdN7MFFAzTiUEARbeS0oWq82pcu2sUI0EJj5xa7dk5tAqzyIKSCMOkvi6+5eKAor/e77UO8/aCyIa6E0Ic15LU5SwrqQcxh2S+uKieW80AKKYGRoOSs4Cl8npXWnYDliXw6B/8DF53aNIpEKvVtiMSEHtmOeQ102nB7gW7fst2Sxa/uZUkUcAbDk/hB6dWQam/n9VqpGxxYCLtZ9cSFwcUM3387la6ZKe+7XKNZBy3M5a6bzkI0Aza1YeUuZAA4NyaLQ5sZnOgKA7JvA7d7E5A2o3SgFsprEi9sxwKtWMOLJuKze6I+iwgDQBvODyJpWQe59ezPT2PZsnYgf4Am83OxcHfMH+93wPSlzZzUEQBk9FilhLLWOq6OAjdT//dconDWVscUrblEFEkiAJBRBGxkelsM8JqNOLiiwTE3mUr1QlIByQRskic9uxeFdK9TGUFgL0TlutwccvfNQLVKLqV/N8OhIsDXJaDz1NZ5xM57BwLlRTqOW6lHhTBAeiqz9SZRCYSnFuzUhrTqo6IIjrXJBqUsGGPTu1GQNqN3EDMIaxIPREH3TChGbSkEaAXYUXC0pZlOXgFpLtdiV/OVKx0yFW/kcrb4iAPgFuJEPJJQsgKIeSE69g4IeS7hJDT9t9jru99kBByhhBykhByl+v4jYSQ4/b3PkbsKhxCSIAQ8kX7+BOEkH3t/RXrwxZWv7ftnt/MYafLpQQUM5a6PbpR7oG1xcThyrmRolsprzuzEwBrt7vR4bna1XD3VqpGRBF7UgTHMmNqxRwA6/qtpCrFoehW6m0HgX4Wh4JuYiNbwHQs4BIH/25IG3l6PgXg7rJjHwDwEKX0MICH7H+DEHIlgHsAHLHf83FCCNuq/CWA+wActv+wn/keAJuU0kMAPgLgT5r9ZZqlXyqkL3mIw1w8aDV6q/PQtxu28HYzTsPE4brdo1hOqt49gIKyIw7ddoGw+6iWKIUUqSftM9gOtZZbCbAy4Nbt6xdziW5IFqGIQl3Lo9OMBCUokuBMQOwnlpN5UArsGA061/HZi5u46yOPVFT9+4G6Tw+l9BEAG2WH3wrg0/bXnwbwNtfxL1BKVUrpOQBnANxECJkDMEIpfZxaEczPlL2H/awvA7iTeNX2d5B+CEjnNQNradXJVGL8u9fuwZd/83V1H/p204s4jSMOe6zGcefWMkipOqL2dDrAWjw2Msz91CPLoYYo9cpyYIN+6t0nkYAElmPgthwIIfibe4/iXTfv7dg5NgIhBFPRQF9aDqyX0lw85FhwD728gpPLKV9Wfjf79MxQShcBwP572j6+E8BF1+vm7WM77a/Lj5e8h1KqA9gCMNHkeTVFryqkKaV44ux6Qxk/7MbaVSYOYUXCtbu9u2x2ErkHMYdkToMiCdgzbtczpFWk81pFmwc2hrPblgOLcdSMOQSsmIPZZSuVNayrJw4lVpirLxQA3HrZFGZGgu0/uW0yFetPcVi0Yzk7RoPO58Am7q2n+9By2CZeO35a43it91T+cELuI4QcI4QcW11tXyFMr4rg/vXkKt7xwI9x/NJW3deyNNZyt1Kv6IUrbiunYSQoO0H4jXTBo7W05HyO3e4g6m68V42I0psUxlzBEvH6AWn3/IneupCq4TdxeObCJr7+7KW6r1vYqrQcmGAwV6ifaFYclm1XEey/V+zj8wB2u163C8CCfXyXx/GS9xBCJABxVLqxAACU0gcopUcppUenpqaaPPVKHBdJl7OVXly0egSxoqNazG/a4jDmF3HogeWQ1xAPSU5tx0amUBmQdn3d9TqHBgLSYWfgT3ddS3nHcqgfkGavk3qcmVSNqVgAaz6KOfzNY+fw3/7lpbqvW0zkEQ/JJdlKjI1+jDlU4UEA99pf3wvg667j99gZSPthBZ6ftF1PKULIzXY84d1l72E/65cAPEy7XP4oOzGH7loOp+z5z5sN3BiXNnMQBYJZH5j1QDFO0+1U1nhIRiwgQRYJ1jMFK+bgshzcLqZeFAbW+38dy6HL6awNB6Rta6EbU/SaZSoawHqm4JsY4VpKbSjJYHErh7m49fxWiIMP3Up17wBCyOcB3A5gkhAyD+CPAPwPAF8ihLwHwAUAbwcASukLhJAvAXgRgA7gfkopu2rvg5X5FALwTfsPAHwCwGcJIWdgWQz3tOU32wZijxrvnV62glCNZCq8vJTE3omwb3Zzcg/Sf7dyGqaiARBCMB5RsJ5WrdbSLmsh5gpOdzsnv7wrqxfMbdPtjCUmRvXcSiwIHfGxOEzGAqDUshynfbBZWkuryBR0UEo952QwLiXy2GG7hfvBcqh7B1BK31nlW3dWef2HAXzY4/gxAFd5HM/DFpde0Qv/uWFSJ0Nhs04LYkopnnp1Ez/7mplunFpD9MKttJXTcMhugT0eCWB+MwdKy4KoLqGolTXUCaKKhJ+/eg6v3T9e9TXumQndpNE6B9Z8L6L4Vxym7A4BKynVF+KwnimAUqubbS1RXdzK4QY70y7oujcFMlgxh4GiXcN+tnIajs/XDy4DwMWNrDMPt57lcG4tg82shhv3jtV8XTfpSSprVnNaZI9HZFzYsPrrlBfBMbptOQgCwf961w04uq+6OLAg7w9OreKa//ztrvnO89tIZQV87lZihXA+iDtohomEvbnLqN6C//SFTSxu5ZDIao7lIImCkyV5+ewIFwe/0q5U1k88dg6//FePN5SayuINAgE2M96Ww4f/+UW85S8ew7HzmwCAG3wkDt22HEyTIqXqLnEIONkf1SyHXvcB8oJZDv9yfBHJvI6LG91pIJdvOJXVjjkE/SsO0z6qknYv6mkPccgWdNzzVz/GvZ98EoCVxsoIySLCiogrZmO+FAf/3gFdpF3tMxYTOeQ0AwXDRKDGUBUAOL1iuZSumB0pCUgnsgUUdBPRoIS/f+ICMgUDf5Y+jVhQclwqfkBuk7XVKKm8DkqLk9cmIopTrFUtIN3rPkBeMHfNK3bLcXcb8k7SaBEcEy9fxxyi/hEH9zl4xZGOz2+hYJg4ZccX2SwWAAjIIqZiAYxHFF+Kg/+enh7AMm+MFhc69gHnC/V/zunlFHbEg9g1FnLMUgD4w68dx5s/9hg+92NLGOIhGZcSOVy/Z6yk4V6v6XbhYDJvXaOi5VCcglciDq6AdLdjDo0QLqsdSOW7M/KS9fAJ1rkmUcet5M8aBwAIKSJiAckX4rBex3J45mICQFHQdrjEYSQo4cBkBOMRBdmC4bsmfP57enpAu5rIrTFx0Ot/yKdX0jg8E8NYWCmxHE4upbCWVvHhf3kJ+ybC+KN/eyUA4MY9/nEpAUWXTbfcSqx1hqc4VKtz8KHlEC6bxJbMdcdyyOsGZJHUzXZj5+fngDTgn0K4tRLLofKzfPrVTeybCOO/vu0q3Lh3DHMut9JH3nEdPvCmK5x72W/Wg7/vgC7RroD0uh0gq5fDbpgUZ1bSeN3BCQgCQSKrgVIKSoGLmzlMRq0in7cf3Y23XLsDS8k8fvH6XTV/Zrcpzt3uvOXwrROLePS0NTfY7VZixAJFa6GXAelGCEoiCIHjEku2yXLQDBOrKdUJeJaTqzPoh+EEpH0ccwCA2XjQiTn1EndCQXlhI6UUz1xM4A2HJnH3VbO4+6rZku+ztjducaj2+fUCf98BXaJdLhKm/PVaI8xvWplKh6dj2LDn9WYLBpJ5DQXdxP/+M4cwEpJw15FZSKKA37r9UEvn1Qm6GZD+42+8iAW7zcCkXR1d1XII+DsgLQgEYVmEqpugaJ9b6ctPzeP//sYLeOIPf9axrgCrPuZ3v/QcdoyGticOPo45AMCe8TC+99Jyr0+jplvpUiKH1ZSK6/fU7n3mV8vBf09PD2jHLjhb0J0hLvXEgQWnDs9EMRa2HuTNbAGv2qMP909G8AvX73KCg36kKA6dtxxSeR1vumoWf/UrN+LQdAwAMBEtioO7B5AoEMc10uvZA9UIByRcPhvDSFBqm1tpOZlHXjOdLDjGsfObeGEhiR+cXK1b4wBYmUCT0QAOz8Tacl6dYvd4GGvpQtX00W6xllIdMS4/l6cvWPGG6+u4hLk4+BhWId1KQNrdVbFeYOn0ivUAH5qOYjRs3RiJrIYLtjjsnQg3fR7dwsnw6nC2kmlSpAs6Dk9HcdeRolnOxqMqklCRGRYLShAIfFNNXs71u0fxxitnEQvKbXMrsQ3Jy0ul4sDmQRcMs6FZDJGAhGP/8Wdx22Xt613WCfaMW8/Ixc3ezpJeyxScc0mXZSs9/soaIoqIy2drC+142J/i4N+taRchhEAWCbQWKqTd5mVdcVhOY0c8iFhQxphbHDayEAXiK79jNWShO5ZDVjOsKugyH/hoSIZASlNXGdGAVJIB5jceePdRAMB3X1pyZmC3CotznSoThyVbHID6aaz9hCMOGzlcMTvSs/NYS6mYiwcRVsQSy8EwKb7zwjLuuGK6bo+vuH0v+00c/Lm16gGSILTUW4kNtQeK7ZGrcWo5hUO22V7iVtrIYsdosOsN45pBllicprOWQ9pePKOuoDNg+e7Hwopn4DQalH0ZjC5nJCgjmWuT5WCLw8kycVhOqmAZ0IMkDrttcbjQpSLCaqxnVExEFXs2eFEcnjy3gfVMAW+6aq7uz2D3st/6K/n/CeoSkkBa2gWvNehWYplKl01bBW1Ft1IBF9Yz2GsPsvE73Zqel1atxdNLBMYjimfKZSwg+TIYXU4sKLXNrZR13ErJkgr95WQeNx+YgCIKAyUOY2EZ0YDUtQpzL0yTYj1dwGQ0gGhALHErfevEIgKSgNsvb8w9xxpJ+gnuVrKRRNLSsB93zKFWQNrJVJph4sAsB8ut9Kar6+80/AAL9hY67FZibhcv99G+yYjnpKhoQOoL62skKLfNrcR6JyXzOpaTKmbt1tArKRU37h3DNbtGMTsSaMv/5QcIIdg9Hu6p5ZDMa9BNioloAJGA5LiVKKX41gtLuO2yqYYrzfdPRiqsvl7DxcFGEoWWgqsbGdXJX69lOZx2MpUst5IsCogFJFzYyGIzqzm+VL9DCIEkkC5YDrZbycNy+Mg7rvN8z9F9Y31iObTRraRZRW6aQfHyUhKz8SBU3cBGpoCZkSB++87Dbfl//MSe8RDO2m1IegGrcZiMKogEJOde3cxqWE6qeO2BxqcdX7t7FN95cRmJbMHxJvQa/z9BXUJu0a20ni5gJmbt1moVwZ1yZSoxRiMyHj1tjT3d2yfiAFjWVqfbnBdjDt6BZ6/j733DAXzsndd39LzawUhIQqZgtEVgswUDV85ZgVm2A2UVxH4ZENVudo9ZlkOXZ4M5bNgNM8cjCqIuy4HFHyejjS/y19kFcc832NW5G3BxsBHF1nbBa5kCZuJBSAKp2T7jzHIac/EgRlw9gMbCCpaTKm7aP47bL59u+hy6jSwKKOidtRxSanVx6HfYPeDVk2e75DUDs/EgDkxF8OWn5lHQTSzbmUrTA+ROcrNnIgxVN3vWRoNlF42FlRK3EnMxT0Qav+5X7YwDAJ6fT7T5LJuHi4ONLAgt7YLX0yomIgpCslgzW+nUSqrEagCAt163E++9ZT/+7j2vRUjpn6Ch3KIrrhGYT94tpoMCm2DXjkK4bMFAWJHwoTe/BqdX0vjEY+ec2eQzA2o57LLnqc8nOtdG48lzGzhxyXs3z3qiWZZDMSDNRMNdxV+PeEjGgakInr3ILQffIYmkpQrpjUwBExEFAVmsGpA2WaZSWfXpe27Zj//4b67sCz+5Gyvm0B23UsTHXUKbhfWJakfGUk6zeifd+ZoZ3HVkBn/20Cm8uJAEMLjiwDrwZjs4cvWP/+kF/L/fPun5PbcIRBSX5WAfn9iGWwkArt01iufmEz1zk5XTX6tRB5GE5nfBlFopbRPRAEKKUDUgPb+ZQ14zcXjaP3MZWkGRBMd10SnSqoaQLPq22rkVmDXUFnEoGE4F9O+98XLkNROf/fGrUETBqaUZNNjv28mRq1nVqDqtbzNTQFgREZRFRAIScpoBw6Ql7qbtcO2uOFZTaknhYi8ZvCeuSVimRzMkczoKhonJqOVWqiYOrG2G3/vWNMrbrtuJ759cxbdfWOrY/5FWdd93CG2WdrmVKKXIaYbTU+rwTAzX7h7FVk7D9Eig5tD7foa5YOv1MmuFvGaUpKm72cgWHAFgMbFMQcdGpoBYcPu1NmxdOL/W28I+BhcHG1Fovs6BtQ6ei4esmEOVm5U13CuPOfQrv33nYVy1cwQf/Orxjg2tSeV1zxqHQSDeJrdSwTBhmLQkXvX2G60W74PqUgKKsyeydVrkt0JeN7GRKXi6ejYyBSeuwOoZMqqOddvFvF2YwHTSEtoOXBxsJFFouhXEIhOH0aAVc6hys55eSWF2JFjSUrmfUSQB999+CBuZYkfZdjPIlgNzK7VaCMcmD7ob6/3ba3cgIAkDm8YKAGHZui/qzU9pBdUe+5vyyCjbzBQw5oiDde0zqo6NjLqtYDTD+Rkd/H22AxcHG7mFnP1LCctHuHPUshzyZemdiWwBf/qdk3j09JpTGT0oxMPt85t/7olXccaerc1I5/WBTGMFioV9rRbCZTVr4XJbDvGQjI+/6wbcf4f/ZoG0i6BiLV8ddSvZz/KGh2tpI1vAuH3/s3s0rVpuqPFtpLEy3NaHH+DiYNNK473FRA6SQDAZDVjiUKb833h+ER97+AxGghJ++ejudpyub3By9Vvc/eqGiQ997QQ++/j5kuNpdXDFQRQIogGpZcuB7ZzLR5De+ZoZXLmjdx1LO40iChAF0jE3jGa764DSrsuMzYzmshxadyux+S1+EYfBfOqaoJXGewuJHGZGghAFgqAsVOxkWEOtb/3OrX3R82c7sIW71QWOzYgu75WTyg+uWwmwhsy3anUxn/sgNdZrBEKsqXr1uiA3i+ryAKynVWiGCZNSBCQRqm4grerOLAb3c7CZKWw7jRUAIl2IoWyHwVqpWqCVxnsLW3nstGcwhJTKgHQiqyEW7I9mcNuFZdy0WuWbqCIOaXVwA9JAe/orsey4csthGAgqInJaZ3ba7qzDjUwB7//7p/GWP/8hknnNmRdSbjksbeWgm7SpmIMkCghIgm8sh8FbrZpEEgVoTdY5LG7lMDdqBf6CHqms7qyGQYPt6lvNVmKWw8XNHExbpCmlAx2QBqwgZPlg+u3CdpqNTHobNMKK2LGdtvs5Xs8U8NzFLZxcTuH+zz3ttOxg7iNWS3LCLjxsxnIALJFp9X5oF4P71G0TuclqX9OkWNrKO9PbvOocNn3UabHdBCQRiiS07layd2IF3cRKymo5ndcsn2/5oJ9BIhKQkGw15mDfb/3UeqVdhOROikNxs3hxI4ulZB6vmRvBo6fXsGf8AoCi5TAaVvCauRF864RV89NMQBqwxa6DFd/bgVsONmKTAem1tArNoNgRL1oOmkFL0mI3XVkNg0gsIHmm+m2HRK4Y8GOupVSNQT+DQjQgIdvitcsNseUQUqoXnTbC919ewd0ffcSzgaTqaqD5k/MbAID77ziIWFDCPz2/CKC0f9KdV0w77tVmAtKAdT+0oxFjO+DiYNNsKuulRLEADig+oF1ySI4AABxQSURBVO4bdjOjbbuUvp+IBVvPuHHPfH513erRX2y6N7jiEHb15GmWnBNzGNzrVI1W3UrPXNjEy0spLHg073NbDq/YcyMum4nhlkOTjhvU/VzfcUWxo3KzbuROusm2CxcHm2ZnEyxuWTUOzK0UVJg4lFoOYwMacwCsnX26xZgDEweBwBn9WGuWw6BgdfPkMYdmCclSS4spm9vsJQ6qLbrs/hMIsHcijFsvK47+HHV5BK7bPeqIQrPi4KeYAxcHG0lorkKa3VQ77IB0ueWg6gayBWNgm58BQCzQ+rjLrZyV0TUXDzlupfQAz3JgRALW4tZKJ878EMccwi26lTbtgT0LW5XN7lgqK8tE3D0eRkASHXEYKctAFAWCOy6fxlhYbjqtONIGS7JdcHGwabb99PxmDmFFdFpiBOXSqs3ylLdBJBZs3U+6ldMwGpaxd6I4F5gJziDHHCIBCbpJS3Lqt0u2oEMUiDPXe5iwAtLN33usg6q3W8l6htnG7+CU1d1g52gIB6cintbBf/z51+Bz77256fMJB0RkeEDaX0ii0FSdw4uLSVw2E3M6X5ZbDs227+0nom2JORQwGlKwZzyMCxvWg8oEJzbI2UptKHzKFUyEZXFgu6/WItSij54N7FncyuGpVzfxlr94zNm5s4mOzGV8cCrivO8P7r4Cv3V7ZWuSsYjSUlV6lLuV/Icskm3XOZgmxQuXtnDNrrhzjIkDyyBhN98gi8NIUG65ziGR0xAPydgzEcZaWkUyr+HcWhqiQAZ2zCUAhNvQTyen6U6sa9gIK9UbXTYC27xdSuTx3ReX8fz8Fk7b/b1Y3JCJw4GpYl+0u47M4pd/qv2tcMKKxFNZ/YYkCKAU27Iezq5lkCkYuHpnURyCZT3mmU9zLDK4u1+WfteK33wrpyEelp1r+dzFBJ6f38JlM7GBbgtRbNjWgjgUjKGsjgaszZhu0qbihZRSZ/O2kMjhpUWrgI25NVlAmg3nOtKFPlURRUTBMDs+m70RuDjYSLa/djs32fFL1jDwa3aNOseCUqlbyZkzO8CWQywowaStuUa2shpGQzKu2z0KQoBj5zdx/NIWrnVZZYNIpA09/LOuKXDDRqgFt1ymYEAzKESBYCGRw4tMHOxUataR9ZbDk/je/3lbyXPeKdpxP7SLlsSBEHKeEHKcEPIsIeSYfWycEPJdQshp++8x1+s/SAg5Qwg5SQi5y3X8RvvnnCGEfIz0wHkqCdZ/uR3L4fh8EkFZKPFFhspSWTdts3VQK6QBdwuN5m5oSikSdkA6FpRx+UwMDz63gERW68oD2UtYzCHdgishpxlDmakEFGs7mnEtsWfz8HQU2YLhtMRglgPb4AUksWsDuvw006EdlsMdlNLrKKVH7X9/AMBDlNLDAB6y/w1CyJUA7gFwBMDdAD5OCGF39F8CuA/AYfvP3W04r23BZhRvJ2Pp+KUEjuyIl8w3dmIOjuWgIRrY/sjAfoINek+rzcUd0qoOw6QYDVkCesPeMZxbs3Zv1wyL5dCiW2lYLYdwC6NCWbzBHUBWJMEZXKXqJmSRQBS6t1dtx/3QLjqxYr0VwKftrz8N4G2u41+glKqU0nMAzgC4iRAyB2CEUvo4tZzWn3G9p2uwNMBGg9KGSXHiUrIk3gB4B6RHB7jGAYDTNbXZHkEs3ZelA9+wxzI2A5KAy2cHY952NdoSc9CGN+bA4lHNuGFYAdyRHcVn+A2HJp0izLxmOG7ibhFR2tPluB20Kg4UwHcIIU8RQu6zj81QShcBwP6b1ZTvBHDR9d55+9hO++vy4xUQQu4jhBwjhBxbXV1t8dRLkYTtWQ5LyTxymlGxeAXK6hw2s4PbkZXhtO3ehjgksgWcWk4BKHZkZVPlbthjuZKu3DEykG3O3bBFvaVspYIx0EH7WjiWQwtupatsy2HnaAhX74pjMZmHqhvIayYCXb6u3ZiL3SitPnmvp5TeAOBNAO4nhNxa47VethmtcbzyIKUPUEqPUkqPTk1Neb2kaVjModGAdKJKimpAsqZTsYd9MzO4HVkZsSZmIf/5w2fwrr95AkBRHEZty2H/ZAR7xsN4/cHJNp+p/3AmiLVS5zDElkMriylzK102E4MkELxmbgR7xsOg1CpuVXUDgS67g/00KrSl0lNK6YL99woh5GsAbgKwTAiZo5Qu2i6jFfvl8wDcicG7ACzYx3d5HO8qLFup0YB0Mmc3hQuVXkJCCHaNhfCqbZpuZjXsn4xUvH+QaGamw0Iih9WUCsOkjluJiSghBN/+nVuHouK3fDPRDMOcrRSUm485bGYLEAWCeEjGb952EDfuHXOe5wvrWaia6XQ86BbFzULvxaHp35wQEiGExNjXAN4I4ASABwHca7/sXgBft79+EMA9hJAAIWQ/rMDzk7brKUUIudnOUnq36z1dI8BSUPXGbjLHFRKqjCccnIrilZU0DJNiJZXHRHRwi7iA5qbBsZm8yZzmtOt2x2ZCilgS6B9UCCGIKGLT4kAptdxKQ245NONW2shoGAvLEASC37vrctxxxTR2j4cBWBlLea377jonW8kHhXCtWA4zAL5mZ51KAP6eUvotQshPAHyJEPIeABcAvB0AKKUvEEK+BOBFADqA+yml7Aq8D8CnAIQAfNP+01Vm7XkMC4kcrpitX+zC5v6OBL3EIYIfnlnDubU08pqJKwY8qMqCaNsJSDOTPpHTKgLSw4bVMqG5xSCt6igYZtPzA/odlsrajFspkS1UuIWnogGEZBGvrmeR13sgDsoAuJUopWcBXOtxfB3AnVXe82EAH/Y4fgzAVc2eSzvYPW6VyF/cqGzA5UWyLIjq5uBUFKpu4jsvLgNAS71W+gFRIFaV9DbEYT1t5ZQnsgUksgUEZWF4g6qB5jtxrti5+dOxYDtPqW8ItZjKWt4QkxCCnWMhLCRyUDWz6zEH5h4clDqHgWAqGkBQFpw0tnokcxoIAaIeA1YO2gUzDz67AFkkODw92JYDYO1+G405GKZV9AZYlsN6uoCJJscqDgKRFiyH5aTVano6NpzXr5g6vn1xtSY0Vlpck1EFa2m1J5aDIFhuxkGtc+hLrEByGPObjVkOWzkNsYAEwaNAhrX2fXkphctmYgNdAMfYzjS4zWwBrA3TVlbDeqbQ9ED2QSAaaD7mwKp6p0eG03JQJAGSQLblVvrKU/N47X/7Hl5dz3q20p+MBixx6EFAGrAtyX4OSA8iu8ZCuLjZoOWQ1z1dSoA1BYoN97lybrBdSozZeNAZmVqP9XRxXvRWTsN6Rh1anznQ2qjQlSQTh+G0HIDtt+1+7Mwa1tIFqLqJPXYA2o0lDgWoeveL4ADYCQrcreQrdo+FG3YrbeU0z2A0g7X37UYnRz9wcCqKV1bTDXVmXc+ozteJrIaNdAHjQ+xWaqWH/3Iyj6AsOFXqw0hI3t40uLOrabzu4AR+/ME78eu37Kv4/lQsgLSqI5HVnKLWbmJNB+SWg6/YPR5CMq87aaq1SNrzB6rBmvEd2TnYvYEYB+3mZUvJPH7vH57Db3/+maqvZZlKgOViWssUMDnEbqWwvVNcT6uYb8ByNU2KD371eTx3MYGVlIqZkeBQDvphhLdhOVBK8cpqBgenopiNB50UdjfsXkzldc/vd5qI0vpkxXbAxcHF7jHLxGzEeqhnOVy7exTRgITXDIlb6ZBtKZ1aTuPbLyzhn55fcPzh5TBxCMkiLiVyKOjmkMccLLfSb/7dU7jlT76PX/nEEyUCWs5iMo/PP3kRX3vmEpaT+aENRjNCitSwOKykVKRVvaSTcjmTrrqkXmTQjYZlJ727l3BxcMEKYBrZvSXztS2He35qDx79gzucxmqDzsFp62H7zgtLSOV1mBT4xnPehe5rdsxh70QYZ1etqVvD7FaKBCSouoljr27iitkYHj29hmPnN6q+nm1eXlxMYjWlDm0aKyMkC8hpje20X7HvN/dUt3JKxaH7S+SEHfPoNVwcXBQth/qB1a2cVtE6w40oEM9MiEFlKhpALCjhQVsQJqMBfP3ZS56v3cioGAvLmIwGnPbIw2w5sCpfSoH777DmEtfK/GLzBl5aTFqWwxAHowHg8HQMT7+aqGltMV5ZtVrBH6wlDi5LrBdupYmIgo2MCrOJmfbthIuDi3hYRiwo1c1YYh0bh7Wi1wtCCA5ORZHK64iHZNx36348N7/lWAZuNjJWp9p4WIZuPwDDnK3ErMvJqIKfPjgBoHafqnlbHFJ5HZmCMfSWw2/cuh953cDf/vBc3de+spJGRBExU0NQ3fdibywHBSaFUwvUK7g4lHFwKorHzqzV7M5abLrHxcEN241dv2cUt11mdWp/YSFZ8bo1u+ht1HX9Br3/VC1Ys7XbLpt24li1ApIXymJiwx5zODQdw11XzuJTPzrvtLWpxiuraRycjtYM4Adl0ekX1ouYA3sWNjLeMbtuwcWhjPffcQhnVzP4zOOvVn0NuwG55VAKizvcsGcMc6PWbnZpK1/xOmY5uBvtDbXlYC9Et18+BUUSEJCEmm6li5s5XGvP2gaAmSEtgHPz67fsRyqv40dn1mq+7uxqBgca6JI8ZS/Q3W6fAQCT9rPQ67gDF4cy7nzNNG67bAof/e4pZxgIYE2FYrs5lupaK1tpGGENBo/uG0MsICGiiFisIg4TUcUZCxpRxKHtqwQArzs4gT9+6xHcfdUsAGs+Rq0mhhc2srh8Jor9E9YiN+wxBwC4fMa696rFC5+5sIlf+PgPcSmRayiDkAWle3Ffjtvxt3UuDv6CEIJ/f9sBpFQdxy9tOcf/8KvH8c4Hfgyg2HSPu5VKuf2yaXzuva/FTx+YACEEs/EglpKlD6thUmxmC5iwYw7AcLuUACvo+e6f3udMvRsJVu9TlSsYWE2p2DMedha5YXcrAVa8cCQoVbjcGF86No+XFpP4wJuuwK++fl/dnzcZsxbonsQc7My99R67lYYjz3KbsKwlt0vkpaUUXlpM4uWlpGuWA798bgSB4PWHitPb5uIhLCRKLYf1jApKrRYjLOYw6GNUt0utPlUszXr3eBg7RkO4sJHl7k2b3ePhqskkGxlLUH/ztoMN/SzHcuhBttJYWAYhvbcc+OrmATPTl5LFhW3B7hv04LMLmBu12ntzy6E2s/EgHjtd6gN+0Q5QXzYbc+Z2D3N1tBexoFzVcrjoEocb9ozhF2/Y5fm6YWT3WBinV1Ke39vMaBWzG2rBxKEX7TMkUcBYWOm55cDdSh4EJBGTUcXxl6fVYkuNbzy/gC17fjSPOdRmLh7ESioP3ZX59cyFBAQCXLNr1AlIc8uhlFqWwwW7LoRZt5wieyasrspe9QEb2e11/nXEoQeWA2A9E722HLg4VGE2HsTSlmUtLNpWw62XTeHiRg7/cnwJAWl4h9M0ylw8BJMCq+niDuiZiwlcNhNDNCA57pBhjzmUU00cDJPi4ZOriCgit7Y82D0WgqqbJfcbYzNTOfWtFjftH8O1u+LYM9EbEZ7g4uBfZkdCjuUwb4vDe2/Zj6t3xvHiYpK7lBpgzh69uriVB6UUpknx7IVNXL9nDIC1O9o/GcE1Q9KcsFGiAW+30n/6+gk8cmoVv3fX5UPdaK8au8a9e6OZdhLEdizUQ9MxfP39t/TMOzAZDWCNB6T9yWw8gGOvWv1tWLzh8EwUX3nf6/DXj56F6DHkh1MKm8v91PlN/Pbnn8Ev3rALybyO6/eMAgBkUcD3f+/2Hp6hP4kFrclwhkmd++zEpS187okL+I037MevvX5/j8/QnzjtbzazOLpv3Dm+ldNgUmzLcug1E1GloXYgnYSLQxXm4iEkshrymoGFRA6SQDAdC0IUiNP/hlMbZjk88OhZrKZUfOyh0wCAG2xx4HjDqnPTroFS3zyxCFEgeN/t/N6rxq4xK1Hkwnpp+vSGHSPsp/5d4xEFiawGzTCdFOduw91KVZgdKVb4LiTymI0HubWwTeIhGUFZwGpKxdU745iIKIgFJRyYrN70jFNMdGCV+JRSfPP4El67f5wH72sQlK2eSeXprKyYtb8sBysOt9lD64FbDlVw+8svbeaww05f5TQOIQRz8RDOrWXw7p/eiyt3jGAlpXrO3eYUYZYDC0qfWk7j7FoGv3YLdyfVw2ua47q9wPaTsLpbaPRqPjgXhyowf/lSModLiRxu2j9e5x0cL2ZHgljcyuFNV88hGpBwpNcn1AfEbMuBBaW/eWIRhAB3HZnp5Wn1BXsnInjszGrJMcdy6CNxYJbDmkfmVbfg4lAFJg6XNnNYSuaxY5Q3N2uG+249gLW0OjRDj9pBueVw4tIWLpuODX1r7kY4OB3BV56eR1rVnXuOxRzG+8ittNOOn9QbH9BJ+BNbhbAiYSQo4bn5LRgm5W6lJrnjiulen0Lf4QSk7UaP85s5Z0ohpzYsnnV2NY1rdlmJDxvpAkKyiJDSP3VJcyNBKJLgDMPqBTwgXYO5eAjffXEZQLHrI4fTadxuJUop5jdzTiYOpzaH7Lbxr7iGTG1ss8bBDwgCwd7xMM6tZXp2DtxyqMFbrtuBY+c38M6b9pTkTXM4nYRZDsm81bYlrepcHBpkz3gEokBwdrW4qG5mChiL9F/R6r7JCM5zcfAnvJ6B0wuCsghFtAb+zG9aOfu7eC+lhlAkAXvGw3hlNQ3TpDAoxUZWw3ik/1q07J+M4AenVmGatCcZflwcOBwfErNnOrAW3dxyaJyDUxG8spLBh/7xOJ44t4FcwcD+HvVIaoW9E2EUdBOLyTx29iDmycWBw/EhrPkesxx4F9bGOTj1/7d3/zFS1Gccx9/PccfBcccdHnf0OH63BLgiiNKCLbYRpS2NjTU2KdQU1DbWtqbtHyZiY9L+YxpJa4xiQolisG2UNrYpVNSoaYlYLEoh/PACApIeFJFf8rM2ok//mO+2W/YObu9md2buPq9kcruzs999vk/m9tnv7Ox3alm/+wh7jpzhwzBDa5ZOY83JXelv/9GziRQHfSEtkkK1/x05/Iu66kqG6sJS3fbxplo++NAxYPqoaFLHLJ3GmjMuXOs6qS+lVRxEUqiuuiqMHM7ROmywZmEtwoSm6E31phmt3H/T5QyosEyeCvyxoYOorqxI7EtpfRwRSaGW+kG81H6Y+poqJo0YmnQ4mTJtVAPfnjOeb10znpb6wWy459pM/oCwosIY21jD/mNRcXjslX1s6XiPR79xZXlevyyvIiJFuX3OeE69f56O4/qNQ7EGVlZw3w1ttNRHeWupH5zZSTMnNtex4+ApPvrIWbVxP89uO1S2qbxVHERSaGprPTdMawHI5CERice8thG8c+p9nn69g47j0ckJm94+VpbXVnEQSam7vzCJ1obBuv5FP3Z92wiqKyv42bp2IBoVvbbveFleW985iKTUuOFDeHXJ3KTDkATVVlcyd3Izz+14h8tb62moqWLj3n42cjCzL5nZLjPbY2ZLko5HRCQNvjJ9JADXTWlm9oRGdh0+zbEyTOWdiuJgZgOAR4H5QBuw0Mzako1KRCR5108ZwV3XfoJbZo1l9oRGADa9XfpDS2k5rPRpYI+77wMws6eBG4E3E41KRCRhAysruPuLkwBoqKli7uRmaspwfZS0FIdWoCPv/gFg1oUbmdkdwB0AY8aMKU9kIiIpUTWggpW3fqosr5WKw0pAZyche8EK9xXuPtPdZzY1NZUhLBGR/iktxeEAMDrv/ijgnwnFIiLS76WlOLwOTDSz8WY2EFgArEk4JhGRfisV3zm4+3kzuwt4ARgArHT3nQmHJSLSb6WiOAC4+zpgXdJxiIhIeg4riYhIiqg4iIhIARUHEREpYO4FPyfIBDM7Dewq4in1wMkYQ4i7veHA0RjbS3t/424P4sthFvpainbj3AfTnsM05w5K098qohjHuvulfyjm7plcgDeK3H5FzK8fd3tF9acP9DfW9uLMYRb6WqI4Y9sH057DNOeuVP0tNsb+dFhpbcrbi1va+5vm/GWlr8phetqLW+L9zfJhpTfcfWbSccSlr/UnCcph7yh/PZeF3BUbY5ZHDiuSDiBmfa0/SVAOe0f567ks5K6oGDM7chARkdLJ8shBRERKRMWhRMxstJn92czazWynmf0wrL/MzF40s7fC32FhfWPY/oyZLctrp87MtuYtR83soaT6VU5x5TA8ttDMtpvZNjN73syGJ9Gncoo5f18PudtpZkuT6E859SB388xsc9jHNpvZ3Ly2rgrr95jZw2bW2SUK0ifO06W0/N+pYy3AleF2HbCb6BKoS4ElYf0S4IFwewgwB7gTWHaRdjcDn0u6f1nKIdEcYu8Cw8P9pcBPk+5fhvLXCPwDaAr3VwHXJd2/lOVuBjAy3J4KHMxraxNwNdF1a54D5ifdv+4sGjmUiLsfcve/h9ungXaiK97dSPTPRfj71bDNWXffALzfVZtmNhFoBl4pYeipEWMOLSxDwqe2ofSD64XEmL8JwG53PxLuvwTcXOLwE9WD3G1x99w+tRMYZGbVZtYCDHX3jR5Viidzz0k7FYcyMLNxRJ8s/gaMcPdDEO2ARG/23bUQWB12sn6lNzl09w+A7wLbiYpCG/B4CcNNnV7ug3uAyWY2zswqid7cRl/iOX1GD3J3M7DF3f9NVFAO5D12IKxLPRWHEjOzWuAZ4EfufqqXzS0Anup9VNnS2xyaWRVRcZgBjAS2AffGGmSK9TZ/7n6CKH+riUat+4HzccaYVsXmzsw+CTwAfCe3qpPNMvHhTsWhhMKb0jPAb9z992H14TDUJPx9t5ttTQcq3X1zSYJNqZhyeAWAu+8No67fAp8pUcipEtc+6O5r3X2Wu19NNKfZW6WKOS2KzZ2ZjQL+ACxy971h9QGiyx7nZOYSyCoOJRKObT8OtLv7g3kPrQEWh9uLgT92s8mF9LNRQ4w5PAi0mVlusrF5RMeQ+7Q490Ezaw5/hwHfAx6LN9p0KTZ3ZtYAPAvc6+6v5jYOh55Om9ns0OYiuv8/n6ykvxHvqwvRWR9OdAhja1i+THTmx8tEn7xeBi7Le85+4DhwhugTR1veY/uAyUn3K6s5JDoDpz20tRZoTLp/GcvfU8CbYVmQdN/SljvgPuBs3rZbgebw2ExgB7AXWEb48XHaF/1CWkRECuiwkoiIFFBxEBGRAioOIiJSQMVBREQKqDiIiEgBFQeREjCzO81sURHbjzOzHaWMSaQYlUkHINLXmFmluy9POg6R3lBxEOlEmGzteaLJ1mYQTdm8CJgCPAjUAkeBW939kJn9Bfgr8FlgjZnVAWfc/edmdgWwHKgh+iHU7e5+wsyuAlYC54AN5eudyKXpsJJI1yYBK9x9GnAK+D7wCPA1d8+9sd+ft32Du3/e3X9xQTtPAveEdrYDPwnrnwB+4NF8RSKpopGDSNc6/H/z5Pwa+DHRhVxeDBfzGgAcytt+9YUNmFk9UdFYH1atAn7XyfpfAfPj74JIz6g4iHTtwrllTgM7L/JJ/2wRbVsn7Yukhg4riXRtjJnlCsFC4DWgKbfOzKrC/P1dcveTwAkzuyas+iaw3t3fA06a2Zyw/pb4wxfpOY0cRLrWDiw2s18SzcL5CPAC8HA4LFQJPER0WciLWQwsN7Maotl1bwvrbwNWmtm50K5IamhWVpFOhLOV/uTuUxMORSQROqwkIiIFNHIQEZECGjmIiEgBFQcRESmg4iAiIgVUHEREpICKg4iIFFBxEBGRAv8Bg3Vr5PFWp60AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'][-200:].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": {}, + "outputs": [], + "source": [ + "first_september_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": 15, + "metadata": {}, + "outputs": [], + "source": [ + "year = []\n", + "yearly_incidence = []\n", + "for week1, week2 in zip(first_september_week[:-1],\n", + " first_september_week[1:]):\n", + " one_year = sorted_data['inc'][week1:week2-1]\n", + " assert abs(len(one_year)-52) < 2\n", + " yearly_incidence.append(one_year.sum())\n", + " year.append(week2.year)\n", + "yearly_incidence = pd.Series(data=yearly_incidence, index=year)" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 16, + "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": 17, + "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": 17, + "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 +2281,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.3" + "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 } -