diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb index 0bbbe371b01e359e381e43239412d77bf53fb1fb..64d6880597cce0dce2c61e56cf287d79c0689f97 100644 --- a/module3/exo2/exercice.ipynb +++ b/module3/exo2/exercice.ipynb @@ -1,5 +1,2291 @@ { - "cells": [], + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + "%matplotlib inline\n", + "import matplotlib.pyplot as plt\n", + "import pandas as pd\n", + "import isoweek" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [], + "source": [ + "data_url = \"https://www.sentiweb.fr/datasets/incidence-PAY-7.csv\"" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020225275239272177578412FRFrance
120225176226382286309513FRFrance
22022507659031001008010515FRFrance
320224975095321269788511FRFrance
420224874985304369278511FRFrance
520224776087373384419513FRFrance
62022467303313924674537FRFrance
72022457382717205934639FRFrance
82022447427122316311639FRFrance
920224375863330284249513FRFrance
102022427377019505590639FRFrance
112022417417722196135639FRFrance
1220224074883147282947212FRFrance
13202239720413313751306FRFrance
14202238717714193123315FRFrance
15202237717254992951315FRFrance
16202236710691781960213FRFrance
17202235715814002762204FRFrance
18202234722667883744315FRFrance
192022337734001739911026FRFrance
202022327780140861151612618FRFrance
21202231768964170962210614FRFrance
222022307903957701230814919FRFrance
232022297148511006019642221529FRFrance
242022287154711102819914231630FRFrance
252022277211911619826184322440FRFrance
262022267168541280620902251931FRFrance
272022257222461801126481342840FRFrance
282022247224581810526811342741FRFrance
292022237187721487522669282234FRFrance
.................................
16441991267176081130423912312042FRFrance
16451991257161691070021638281838FRFrance
16461991247161711007122271281739FRFrance
1647199123711947767116223211329FRFrance
1648199122715452995320951271737FRFrance
1649199121714903897520831261636FRFrance
16501991207190531274225364342345FRFrance
16511991197167391124622232291939FRFrance
16521991187213851388228888382551FRFrance
1653199117713462887718047241632FRFrance
16541991167148571006819646261834FRFrance
1655199115713975978118169251832FRFrance
1656199114712265768416846221430FRFrance
165719911379567604113093171123FRFrance
1658199112710864733114397191325FRFrance
16591991117155741118419964271935FRFrance
16601991107166431137221914292038FRFrance
1661199109713741878018702241533FRFrance
1662199108713289881317765231531FRFrance
1663199107712337807716597221529FRFrance
1664199106710877701314741191226FRFrance
1665199105710442654414340181125FRFrance
16661991047791345631126314820FRFrance
16671991037153871048420290271836FRFrance
16681991027162771104621508292038FRFrance
16691991017155651027120859271836FRFrance
16701990527193751329525455342345FRFrance
16711990517190801380724353342543FRFrance
1672199050711079666015498201228FRFrance
16731990497114302610205FRFrance
\n", + "

1674 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202252 7 5239 2721 7757 8 4 \n", + "1 202251 7 6226 3822 8630 9 5 \n", + "2 202250 7 6590 3100 10080 10 5 \n", + "3 202249 7 5095 3212 6978 8 5 \n", + "4 202248 7 4985 3043 6927 8 5 \n", + "5 202247 7 6087 3733 8441 9 5 \n", + "6 202246 7 3033 1392 4674 5 3 \n", + "7 202245 7 3827 1720 5934 6 3 \n", + "8 202244 7 4271 2231 6311 6 3 \n", + "9 202243 7 5863 3302 8424 9 5 \n", + "10 202242 7 3770 1950 5590 6 3 \n", + "11 202241 7 4177 2219 6135 6 3 \n", + "12 202240 7 4883 1472 8294 7 2 \n", + "13 202239 7 2041 331 3751 3 0 \n", + "14 202238 7 1771 419 3123 3 1 \n", + "15 202237 7 1725 499 2951 3 1 \n", + "16 202236 7 1069 178 1960 2 1 \n", + "17 202235 7 1581 400 2762 2 0 \n", + "18 202234 7 2266 788 3744 3 1 \n", + "19 202233 7 7340 0 17399 11 0 \n", + "20 202232 7 7801 4086 11516 12 6 \n", + "21 202231 7 6896 4170 9622 10 6 \n", + "22 202230 7 9039 5770 12308 14 9 \n", + "23 202229 7 14851 10060 19642 22 15 \n", + "24 202228 7 15471 11028 19914 23 16 \n", + "25 202227 7 21191 16198 26184 32 24 \n", + "26 202226 7 16854 12806 20902 25 19 \n", + "27 202225 7 22246 18011 26481 34 28 \n", + "28 202224 7 22458 18105 26811 34 27 \n", + "29 202223 7 18772 14875 22669 28 22 \n", + "... ... ... ... ... ... ... ... \n", + "1644 199126 7 17608 11304 23912 31 20 \n", + "1645 199125 7 16169 10700 21638 28 18 \n", + "1646 199124 7 16171 10071 22271 28 17 \n", + "1647 199123 7 11947 7671 16223 21 13 \n", + "1648 199122 7 15452 9953 20951 27 17 \n", + "1649 199121 7 14903 8975 20831 26 16 \n", + "1650 199120 7 19053 12742 25364 34 23 \n", + "1651 199119 7 16739 11246 22232 29 19 \n", + "1652 199118 7 21385 13882 28888 38 25 \n", + "1653 199117 7 13462 8877 18047 24 16 \n", + "1654 199116 7 14857 10068 19646 26 18 \n", + "1655 199115 7 13975 9781 18169 25 18 \n", + "1656 199114 7 12265 7684 16846 22 14 \n", + "1657 199113 7 9567 6041 13093 17 11 \n", + "1658 199112 7 10864 7331 14397 19 13 \n", + "1659 199111 7 15574 11184 19964 27 19 \n", + "1660 199110 7 16643 11372 21914 29 20 \n", + "1661 199109 7 13741 8780 18702 24 15 \n", + "1662 199108 7 13289 8813 17765 23 15 \n", + "1663 199107 7 12337 8077 16597 22 15 \n", + "1664 199106 7 10877 7013 14741 19 12 \n", + "1665 199105 7 10442 6544 14340 18 11 \n", + "1666 199104 7 7913 4563 11263 14 8 \n", + "1667 199103 7 15387 10484 20290 27 18 \n", + "1668 199102 7 16277 11046 21508 29 20 \n", + "1669 199101 7 15565 10271 20859 27 18 \n", + "1670 199052 7 19375 13295 25455 34 23 \n", + "1671 199051 7 19080 13807 24353 34 25 \n", + "1672 199050 7 11079 6660 15498 20 12 \n", + "1673 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 12 FR France \n", + "1 13 FR France \n", + "2 15 FR France \n", + "3 11 FR France \n", + "4 11 FR France \n", + "5 13 FR France \n", + "6 7 FR France \n", + "7 9 FR France \n", + "8 9 FR France \n", + "9 13 FR France \n", + "10 9 FR France \n", + "11 9 FR France \n", + "12 12 FR France \n", + "13 6 FR France \n", + "14 5 FR France \n", + "15 5 FR France \n", + "16 3 FR France \n", + "17 4 FR France \n", + "18 5 FR France \n", + "19 26 FR France \n", + "20 18 FR France \n", + "21 14 FR France \n", + "22 19 FR France \n", + "23 29 FR France \n", + "24 30 FR France \n", + "25 40 FR France \n", + "26 31 FR France \n", + "27 40 FR France \n", + "28 41 FR France \n", + "29 34 FR France \n", + "... ... ... ... \n", + "1644 42 FR France \n", + "1645 38 FR France \n", + "1646 39 FR France \n", + "1647 29 FR France \n", + "1648 37 FR France \n", + "1649 36 FR France \n", + "1650 45 FR France \n", + "1651 39 FR France \n", + "1652 51 FR France \n", + "1653 32 FR France \n", + "1654 34 FR France \n", + "1655 32 FR France \n", + "1656 30 FR France \n", + "1657 23 FR France \n", + "1658 25 FR France \n", + "1659 35 FR France \n", + "1660 38 FR France \n", + "1661 33 FR France \n", + "1662 31 FR France \n", + "1663 29 FR France \n", + "1664 26 FR France \n", + "1665 25 FR France \n", + "1666 20 FR France \n", + "1667 36 FR France \n", + "1668 38 FR France \n", + "1669 36 FR France \n", + "1670 45 FR France \n", + "1671 43 FR France \n", + "1672 28 FR France \n", + "1673 5 FR France \n", + "\n", + "[1674 rows x 10 columns]" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + " raw_data = pd.read_csv(data_url, encoding = 'iso-8859-1', skiprows=1)\n", + "raw_data" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
\n", + "
" + ], + "text/plain": [ + "Empty DataFrame\n", + "Columns: [week, indicator, inc, inc_low, inc_up, inc100, inc100_low, inc100_up, geo_insee, geo_name]\n", + "Index: []" + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "raw_data[raw_data.isnull().any(axis=1)]" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
020225275239272177578412FRFrance
120225176226382286309513FRFrance
22022507659031001008010515FRFrance
320224975095321269788511FRFrance
420224874985304369278511FRFrance
520224776087373384419513FRFrance
62022467303313924674537FRFrance
72022457382717205934639FRFrance
82022447427122316311639FRFrance
920224375863330284249513FRFrance
102022427377019505590639FRFrance
112022417417722196135639FRFrance
1220224074883147282947212FRFrance
13202239720413313751306FRFrance
14202238717714193123315FRFrance
15202237717254992951315FRFrance
16202236710691781960213FRFrance
17202235715814002762204FRFrance
18202234722667883744315FRFrance
192022337734001739911026FRFrance
202022327780140861151612618FRFrance
21202231768964170962210614FRFrance
222022307903957701230814919FRFrance
232022297148511006019642221529FRFrance
242022287154711102819914231630FRFrance
252022277211911619826184322440FRFrance
262022267168541280620902251931FRFrance
272022257222461801126481342840FRFrance
282022247224581810526811342741FRFrance
292022237187721487522669282234FRFrance
.................................
16441991267176081130423912312042FRFrance
16451991257161691070021638281838FRFrance
16461991247161711007122271281739FRFrance
1647199123711947767116223211329FRFrance
1648199122715452995320951271737FRFrance
1649199121714903897520831261636FRFrance
16501991207190531274225364342345FRFrance
16511991197167391124622232291939FRFrance
16521991187213851388228888382551FRFrance
1653199117713462887718047241632FRFrance
16541991167148571006819646261834FRFrance
1655199115713975978118169251832FRFrance
1656199114712265768416846221430FRFrance
165719911379567604113093171123FRFrance
1658199112710864733114397191325FRFrance
16591991117155741118419964271935FRFrance
16601991107166431137221914292038FRFrance
1661199109713741878018702241533FRFrance
1662199108713289881317765231531FRFrance
1663199107712337807716597221529FRFrance
1664199106710877701314741191226FRFrance
1665199105710442654414340181125FRFrance
16661991047791345631126314820FRFrance
16671991037153871048420290271836FRFrance
16681991027162771104621508292038FRFrance
16691991017155651027120859271836FRFrance
16701990527193751329525455342345FRFrance
16711990517190801380724353342543FRFrance
1672199050711079666015498201228FRFrance
16731990497114302610205FRFrance
\n", + "

1674 rows × 10 columns

\n", + "
" + ], + "text/plain": [ + " week indicator inc inc_low inc_up inc100 inc100_low \\\n", + "0 202252 7 5239 2721 7757 8 4 \n", + "1 202251 7 6226 3822 8630 9 5 \n", + "2 202250 7 6590 3100 10080 10 5 \n", + "3 202249 7 5095 3212 6978 8 5 \n", + "4 202248 7 4985 3043 6927 8 5 \n", + "5 202247 7 6087 3733 8441 9 5 \n", + "6 202246 7 3033 1392 4674 5 3 \n", + "7 202245 7 3827 1720 5934 6 3 \n", + "8 202244 7 4271 2231 6311 6 3 \n", + "9 202243 7 5863 3302 8424 9 5 \n", + "10 202242 7 3770 1950 5590 6 3 \n", + "11 202241 7 4177 2219 6135 6 3 \n", + "12 202240 7 4883 1472 8294 7 2 \n", + "13 202239 7 2041 331 3751 3 0 \n", + "14 202238 7 1771 419 3123 3 1 \n", + "15 202237 7 1725 499 2951 3 1 \n", + "16 202236 7 1069 178 1960 2 1 \n", + "17 202235 7 1581 400 2762 2 0 \n", + "18 202234 7 2266 788 3744 3 1 \n", + "19 202233 7 7340 0 17399 11 0 \n", + "20 202232 7 7801 4086 11516 12 6 \n", + "21 202231 7 6896 4170 9622 10 6 \n", + "22 202230 7 9039 5770 12308 14 9 \n", + "23 202229 7 14851 10060 19642 22 15 \n", + "24 202228 7 15471 11028 19914 23 16 \n", + "25 202227 7 21191 16198 26184 32 24 \n", + "26 202226 7 16854 12806 20902 25 19 \n", + "27 202225 7 22246 18011 26481 34 28 \n", + "28 202224 7 22458 18105 26811 34 27 \n", + "29 202223 7 18772 14875 22669 28 22 \n", + "... ... ... ... ... ... ... ... \n", + "1644 199126 7 17608 11304 23912 31 20 \n", + "1645 199125 7 16169 10700 21638 28 18 \n", + "1646 199124 7 16171 10071 22271 28 17 \n", + "1647 199123 7 11947 7671 16223 21 13 \n", + "1648 199122 7 15452 9953 20951 27 17 \n", + "1649 199121 7 14903 8975 20831 26 16 \n", + "1650 199120 7 19053 12742 25364 34 23 \n", + "1651 199119 7 16739 11246 22232 29 19 \n", + "1652 199118 7 21385 13882 28888 38 25 \n", + "1653 199117 7 13462 8877 18047 24 16 \n", + "1654 199116 7 14857 10068 19646 26 18 \n", + "1655 199115 7 13975 9781 18169 25 18 \n", + "1656 199114 7 12265 7684 16846 22 14 \n", + "1657 199113 7 9567 6041 13093 17 11 \n", + "1658 199112 7 10864 7331 14397 19 13 \n", + "1659 199111 7 15574 11184 19964 27 19 \n", + "1660 199110 7 16643 11372 21914 29 20 \n", + "1661 199109 7 13741 8780 18702 24 15 \n", + "1662 199108 7 13289 8813 17765 23 15 \n", + "1663 199107 7 12337 8077 16597 22 15 \n", + "1664 199106 7 10877 7013 14741 19 12 \n", + "1665 199105 7 10442 6544 14340 18 11 \n", + "1666 199104 7 7913 4563 11263 14 8 \n", + "1667 199103 7 15387 10484 20290 27 18 \n", + "1668 199102 7 16277 11046 21508 29 20 \n", + "1669 199101 7 15565 10271 20859 27 18 \n", + "1670 199052 7 19375 13295 25455 34 23 \n", + "1671 199051 7 19080 13807 24353 34 25 \n", + "1672 199050 7 11079 6660 15498 20 12 \n", + "1673 199049 7 1143 0 2610 2 0 \n", + "\n", + " inc100_up geo_insee geo_name \n", + "0 12 FR France \n", + "1 13 FR France \n", + "2 15 FR France \n", + "3 11 FR France \n", + "4 11 FR France \n", + "5 13 FR France \n", + "6 7 FR France \n", + "7 9 FR France \n", + "8 9 FR France \n", + "9 13 FR France \n", + "10 9 FR France \n", + "11 9 FR France \n", + "12 12 FR France \n", + "13 6 FR France \n", + "14 5 FR France \n", + "15 5 FR France \n", + "16 3 FR France \n", + "17 4 FR France \n", + "18 5 FR France \n", + "19 26 FR France \n", + "20 18 FR France \n", + "21 14 FR France \n", + "22 19 FR France \n", + "23 29 FR France \n", + "24 30 FR France \n", + "25 40 FR France \n", + "26 31 FR France \n", + "27 40 FR France \n", + "28 41 FR France \n", + "29 34 FR France \n", + "... ... ... ... \n", + "1644 42 FR France \n", + "1645 38 FR France \n", + "1646 39 FR France \n", + "1647 29 FR France \n", + "1648 37 FR France \n", + "1649 36 FR France \n", + "1650 45 FR France \n", + "1651 39 FR France \n", + "1652 51 FR France \n", + "1653 32 FR France \n", + "1654 34 FR France \n", + "1655 32 FR France \n", + "1656 30 FR France \n", + "1657 23 FR France \n", + "1658 25 FR France \n", + "1659 35 FR France \n", + "1660 38 FR France \n", + "1661 33 FR France \n", + "1662 31 FR France \n", + "1663 29 FR France \n", + "1664 26 FR France \n", + "1665 25 FR France \n", + "1666 20 FR France \n", + "1667 36 FR France \n", + "1668 38 FR France \n", + "1669 36 FR France \n", + "1670 45 FR France \n", + "1671 43 FR France \n", + "1672 28 FR France \n", + "1673 5 FR France \n", + "\n", + "[1674 rows x 10 columns]" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "data = raw_data.dropna().copy()\n", + "data" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": {}, + "outputs": [], + "source": [ + "def convert_week(year_and_week_int):\n", + " year_and_week_str = str(year_and_week_int)\n", + " year = int(year_and_week_str[:4])\n", + " week = int(year_and_week_str[4:])\n", + " w = isoweek.Week(year, week)\n", + " return pd.Period(w.day(0), 'W')\n", + "\n", + "data['period'] = [convert_week(yw) for yw in data['week']]" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [], + "source": [ + "sorted_data = data.set_index('period').sort_index()" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [], + "source": [ + "periods = sorted_data.index\n", + "for p1, p2 in zip(periods[:-1], periods[1:]):\n", + " delta = p2.to_timestamp() - p1.end_time\n", + " if delta > pd.Timedelta('1s'):\n", + " print(p1, p2)" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm8HUWZ9vOec+692ROyQUgCCRCWhCVACFEEWSWICCJo1BkYZAZ1YIRxmQ8GUVBRcMNBBEVQFkcBkRFkE2SVLSRACFkIBLInZN+Te3PvPfX90V2nq6vfWvpsd+vn90vOuXWqq6q7q+qtdychBDJkyJAhQwYVuY4eQIYMGTJk6HzIiEOGDBkyZEggIw4ZMmTIkCGBjDhkyJAhQ4YEMuKQIUOGDBkSyIhDhgwZMmRIICMOGTJkyJAhgYw4ZMiQIUOGBDLikCFDhgwZEih09ADKxdChQ8WYMWM6ehgZMmTI0KXw2muvrRNCDHPV67LEYcyYMZg5c2ZHDyNDhgwZuhSIaIlPvUyslCFDhgwZEsiIQ4YMGTJkSCAjDhkyZMiQIYGMOGTIkCFDhgQy4pAhQ4YMGRLIiEOGDBkyZEggIw4ZMmTIkCGBjDhk6NJYtmEHnn9nbUcPI0OGbocu6wSXIQMAfPTHz6AogMXXnd7RQ8mQoVsh4xwydGkURUePIEOG7omMOGTIkCFDhgQy4pAhQ4YMGRJwEgci6kVErxLRm0Q0l4iuCcuvJqIVRDQr/Pdx5ZoriGghES0golOV8iOJ6K3wtxuJiMLyJiK6NyyfTkRjqn+rGTJkyJDBFz6cQwuAE4UQhwGYCGAqEU0Jf7tBCDEx/PcoABDReADTAEwAMBXAzUSUD+vfAuAiAOPCf1PD8gsBbBRC7AfgBgDXV35rGTJkyJChXDiJgwiwLfyzIfxnUwOeCeAeIUSLEGIRgIUAJhPRCAADhBAvCyEEgLsAnKVcc2f4/X4AJ0muIkOGDBky1B9eOgciyhPRLABrADwphJge/nQJEc0mot8S0W5h2UgAy5TLl4dlI8PvennsGiFEG4DNAIYw47iIiGYS0cy1azPb9gwZMmSoFbyIgxCiXQgxEcAoBFzAwQhERPsiEDWtAvDTsDp34heWcts1+jhuFUJMEkJMGjbMmcgoQwYWbe1FPDx7JQIGNkOGDBxSWSsJITYBeBbAVCHE6pBoFAH8BsDksNpyAKOVy0YBWBmWj2LKY9cQUQHAQAAbUt1Jhgye+PXz7+OSP7yBh2ev6uihZMjQaeFjrTSMiAaF33sDOBnA26EOQeJTAOaE3x8CMC20QBqLQPH8qhBiFYCtRDQl1CecB+BB5Zrzw+/nAHhaZMe6DDXCB5ubAQAbtu/q4JFkyNB54RM+YwSAO0OLoxyA+4QQDxPR3UQ0EYH4ZzGALwGAEGIuEd0HYB6ANgAXCyHaw7a+AuAOAL0BPBb+A4DbAdxNRAsRcAzTqnBvGTKwyEwdMmRww0kchBCzARzOlP+z5ZprAVzLlM8EcDBT3gzgXNdYMmTIkCFDfZB5SGfoNrj1+fdw+Z9ne9fPJJcZMpiREYcM3QY/ePRt3DNjmbNeJlXKkMGNjDhk6LHI+IYMGczIiEOGHofM+T5DBjcy4pChx6KWKgchBK5+aC5eW5K562TomsiIQ4Zugc6oXL7jpcX49C0vd/QwMmQoCxlxyNAtUA5tqCU56YS0KkOGVMiIQ4ZugTR7cT1UDhltyNDVkRGHDN0CxeyoniFDVZERhwzdAmWJlWpIUDqjDiRDhjTIiEOGbgGRQpBDdXCDy0hDhq6OjDh0csxevgmbdmTRQ13IDuoZMlQXGXHo5PjkTS/is79+paOH0a1QF4V0RqwydHFkxKELYMHqrR09hE6P8nQO1R9Hqe1MsJShiyMjDhm6BdLpHGoPSXiySB0Zuioy4pChW6A8J7jan+4z8VL18fTbqzFnxeaOHka3h08muAwZOj06mxNchtrhi3fMBAAsvu70Dh5J90bGOXRiZLby/uhszyoTK2Xo6nASByLqRUSvEtGbRDSXiK4JywcT0ZNE9G74uZtyzRVEtJCIFhDRqUr5kUT0VvjbjRTGTiaiJiK6NyyfTkRjqn+rXQ+dbL/r1CjnUdXj+dbzHY65/BFc99jb9eswQ7eGD+fQAuBEIcRhACYCmEpEUwBcDuApIcQ4AE+Ff4OIxgOYBmACgKkAbiaifNjWLQAuAjAu/Dc1LL8QwEYhxH4AbgBwfRXurWoQQuDulxdj5672+vZb1966NtJswvXI59BR1kq/eu69Duk3Q/eDkziIANvCPxvCfwLAmQDuDMvvBHBW+P1MAPcIIVqEEIsALAQwmYhGABgghHhZBDKAu7RrZFv3AziJOlFGlifnrcZVD87F9Y/X91TW2UQlnRqdNCpr55nFGTKkg5fOgYjyRDQLwBoATwohpgPYXQixCgDCz+Fh9ZEA1ES+y8OykeF3vTx2jRCiDcBmAEPKuaFaYEfIMWyss6dyRhr80dlMWSXqRd+zg0SGasOLOAgh2oUQEwGMQsAFHGypzq09YSm3XRNvmOgiIppJRDPXrl3rGnaXR7be/dH5nODqi2I2VzJUGamslYQQmwA8i0BXsDoUFSH8XBNWWw5gtHLZKAArw/JRTHnsGiIqABgIIJFfUQhxqxBikhBi0rBhw9IMvUuiO3rZbt7Ris07W6vebqonVZfwGcGIqilWWrxuO2Ys5tOOZpxDhmrDx1ppGBENCr/3BnAygLcBPATg/LDa+QAeDL8/BGBaaIE0FoHi+dVQ9LSViKaE+oTztGtkW+cAeFp0wtn+5rJNuPnZhXXrr/M9gcpx2HefwGHXPOFVt7m1Hef+6iUvh6dypkstiW8tWj7+J8/i3F/xaUczziFDteHDOYwA8AwRzQYwA4HO4WEA1wE4hYjeBXBK+DeEEHMB3AdgHoDHAVwshJBmPl8BcBsCJfV7AB4Ly28HMISIFgL4GkLLp86Gxet34EePL6hKW++s3oqtzdU/QXcnzFq2CTMWb8R3H57nrJtmc6xHyO5qYvbyTdi8wz5XuiOXmaFj4fSQFkLMBnA4U74ewEmGa64FcC1TPhNAQl8hhGgGcK7HeDsFhBAVm0N+7IbncdjoQXjw4mMs/VTURY/C+u0tHT2EGKr57j5504s4eOSAuvWXIQOQeUh7QT+V7WovVqXdN5dtStVvBjN+8VR6cV9NN9Qqtz1nxRZ7d9lUyVBlZMShDOxqqw5xcCFb8P5oKvhP5e7oe5Dl0M5QbWTEwQO6jLpuxKEuvXQBeDyIpoYUxEE2W8sc0nV+ezpxWLV5J77xpzfrNlczdD9kxKEMVEus5EInNNjqtGgq5N2VQnS3THC72opYsn5HrOyqv8zB/a8tx7ML1hiuypDBjixktwcSOoeMc+h02HNQr9TXdAXayx0Q2osC+VxE4b71l7dw38zlsTrSekutlyFDGmScQxlor5NR+VLtNNjTUBL/WMjkgF7B+WavwX1TtFuPwHu1Q6vGuf7j3XWJOlLMlOtmCpaMm64fMuLggY6yi/ex7+/pKOSDKVyeE1ztUK1NjGumzXI4kbRAVulmtKFLcHvdBRlxKAN1m5/ZQnAi4i5SXFMPnUMN27ZxrrqyvdtxDh09gB6EjDiUgUpPL76nyszPwR/d9UTJ3ZZt/kjnzO4qVspQP2TEoSzUZyfqrhteLdBWTG8kYHu+yzbswJjLH8Gri6JAd4vXbcfrSzdW3HYacITAR+UluYvupo/OdA71Q0YcykDlnINnvcq66fJIE6Lk0ntm+bfrUeel9wIl759mRqlJjv/Jszj75pe8+kjD9S3bsAO3Ps9ncNvvyscSZbpYSZ1P8t5klVw3ow49fU3UE5kpaxnIJmjnQSWewZ1FbPfPt0/H4vU78KnDR2FY/yZnfZ/Tc7fVOXSOV9YjkHEOZaBizsG7n56xEm77x/tYsn57Wdf6mhW/uWwT1mxtBgA8ncIxrOy9NcWr29bSlqpp/ZbVMerWSt2Mceg0BL0nICMOHtAnZKUT1F8h3f2xpbkV339kPqbd+kpZ1/u6nJz5yxfxsRuex7aWtlIQO5/XUC59Vi97+u3VuP2FRca6ac1O24VZrNRWFGhpay9xVJ0oFXuGLoaMOJSB+uUFrk8/tcCKTTvx0ycWOAmh/Dnt6VkijSJ6045WbNzulwe8mr4tX7xjJr5n8VkpZY3zbK9ooYhCAAd86/Huyzl04TXR1ZARBw/oG0XdxEqVddOhuOQPr+MXTy/EvFX2UNMlWG7W9rzTGimpKUpr6wSXom746XvK19vmLpPh4LubzkHFM2+vwfPvdP9c8h2FTCFdBuom9+zCx6SW1mDXdt1CpXuXLmJxIW3QxHLfQJo5Im/B91HYxEo6uhttUO/1gjtmAAAWX3d6B42meyPjHDyQ0DlkpqxOpB17ufeaNs6Vt5K/wk1VduOn10h3Dz05d0OmkK4fnMSBiEYT0TNENJ+I5hLRpWH51US0gohmhf8+rlxzBREtJKIFRHSqUn4kEb0V/nYjhXw0ETUR0b1h+XQiGlP9W+166A57gOvkasutUItTb4xx6CQPeMeuIMV6LazYqnGLC9dsK1snlKHrwodzaAPwdSHEQQCmALiYiMaHv90ghJgY/nsUAMLfpgGYAGAqgJuJSAbbvwXARQDGhf+mhuUXAtgohNgPwA0Arq/81qqHagfe8z39dOUTou8GVm9rmnpF1E3Tiy2QHoc01atxtyf/7Dl84bbpVWipcnThJdHl4CQOQohVQojXw+9bAcwHMNJyyZkA7hFCtAghFgFYCGAyEY0AMEAI8bIIdo67AJylXHNn+P1+ACdRvXeNFMislaqHevtyqMShs0Vl9b0mDYGr1vO15Tv//G9ewYOzVlSlHxd6wJLoNEilcwjFPYcDkMeIS4hoNhH9loh2C8tGAlimXLY8LBsZftfLY9cIIdoAbAYwJM3Yaonq+zn49tt1kZa21+peN+9sxT2vLi39nVaBXU/a5dtVGo6yHsN/6b31qcKXVIKe4hjaGeBNHIioH4A/A7hMCLEFgYhoXwATAawC8FNZlblcWMpt1+hjuIiIZhLRzLVrO86ErX6cQ+dYCI/PWYXZy80nRw5pHf1qdav/df+buPyBt0p/qz4CtXy85bTte801D81Dm6fVVa2nUL3EdBKdY0X0DHgRByJqQEAY/lcI8QAACCFWCyHahRBFAL8BMDmsvhzAaOXyUQBWhuWjmPLYNURUADAQwAZoEELcKoSYJISYNGzYML87rAF62gT98u9fxydverGsa136Gi9rnrJ6DrBuW9zpbekGv+x6HSHT9OVIX128gc3+1hFQs9LVS7SUoT7wsVYiALcDmC+E+JlSPkKp9ikAc8LvDwGYFlogjUWgeH5VCLEKwFYimhK2eR6AB5Vrzg+/nwPgadFZjs3gnODqJFbqgCfw/+6fjd+9aA71UHVIk886kdzvPDRX6drdZ11NJ1N05X9ir+34W5R86rUWLb20cB2O+9EzNe0jQwQfJ7hjAPwzgLeISL79/wbwOSKaiGD2LQbwJQAQQswlovsAzENg6XSxEKI9vO4rAO4A0BvAY+E/ICA+dxPRQgQcw7TKbis95qzYjL5NBYwdmsxFnNQ51AcdYdN9bxii+oJjxlalPV/VQ+c5CgSo1B6iLLFSirq+eodaJ6ba1ZY+j0a5+N4j87FpR6u7YoaqwEkchBAvgOeyH7Vccy2Aa5nymQAOZsqbAZzrGkst8YlfvADAz9uy8vAZ9VnYXQG2Z9GR5moVc4dlEPZyQm5Uq57xekcDrZruQwhhJKy72opYtG47Dtijf1lj6bTmi90UmYd0Weh8tvK1ArdJtrS148FZK6qiMO9IAtjpFNJlhNyoxThi1zt+1zkHW39X/3UuTv358/hgc3NZY+m8xu3dExlx8EDVA+95L+yOJw+cg9YNT76LS++Zhaff9s+L4EKt7rTc/aQTu9mE8OU+K9WP2a/X54dN3PXa4iDF6qadfpFxdXT6V9LNkBEHA155fz3GXP4IZi/f1GHxXA4bPQgAUOjAuMu62AAAVm8JTn5qhFMdqX05OuAR19QJrpxr0oiV6mVOnbKGrX4pEVGZaoruHGGWw0d//Ax++sSCDus/Iw4GyFPxy++tT/xWsRzXs97w/r0AAMeOG1phj+WDs4pxnSZfX7oRC1ZvBeA+7XUG7siKModXlod0irq+xkq11jnov9s4B7m5lxsWpmeRBmDJ+h34xdMLO6z/jDiUgUr2s6Xrd+C1JRs9+wk6qtf2+dEfJ80Ey+n77JtfSn2NjTurhIDYiJM11HXZPZaPVAH1qmTUsHLTTvx93uqy+9F/tfWXy/mNyYTOL+rrXsjyOTiwvaUN67dVz8/hOGYDNkGesOp1uF6yPukgpvYthMCyDTurukg7Od9QNjqLWMm1uX/yphewbtsuo5VeWs7BShzCeZM2hIlERhvqi4w4GCDn4Y0MW1evDU2KDjpLdNbfv7IEVz04F/swviDlokOtlVJL1FO0XeP78p4Tjmq6B3nichdx8Iw71tzaXrJSMo19zZZmfPbWV3DHBUdh7yHJOZbRhvoiEyuVgXIX/s///k6izObp2imIgjIEKQ57f91278vnrvBLE8rdquuk6Cue6ypI87p9PaQr14/xLWzcvgs7d7UzOge+na/+8Q2s2doStGm40YfeXIlF67bjjpcWs79nYqX6IiMOZaBc66Wf//3dRBlnDVTqR8Q/OwLqvZazOL/+pze92k97i82t7fj0Lel1G1rnTpS/Hfnd0buh4j64wv8ptLXXy0OaLz/8e0/i9F/8w6qQ/t/pSzDm8kewpbkVzym5nk0EpKkhSPvSYvC67qmkYZlnPLBqIyMO5aCKm7Utp3FJ59CBknl18ddkcZZ5a39+fbm7EipP1FRrsdJ3H56X+hoAaFXsQcuZH9Ic2QVby++v3W5NofvbFxYBCMRF/XtFEuyigTo0FYLt6A/Tl+IOJr5XTzNllfjWX+a4K9UAGXEoA9XcqoXF5rveCmkOsa5ruDbTKvmrnZ0v0X4NmnfdY5on4M05MK0+/fZqHP2Dp/CMhxNjS2u79fekQjoqUH/q06gQBxPnUIi2o6v/Oo+vxI6hE4hfa4gN28tzGqwUGXEwoUwTyLSw6RXkIurYEBO17Vxon77IV2Hm1tsJjtsUfTd5Hao40kYoudc3a9nm8NOdo+OHj72dalz8dCHkFUdOk74kX6azZ71zStQbHcUwZcShDHz/Ef9TjQu2aS06g1hJ+V6L03q5+RzyufpM3XKJI3cZ15a6sW3ascuqgzJdZx0HU5bmLa7ctNPevqcTnNpnufPZdB1nGvviwnXdJr9ERwnTMuJQBt7+YKu7kiesnEO4T3SW4HS1PMGkvcdqhBSxbfyV3iu3kXH7+Qpl8/3UzS/hGw4FvgQX84odh+Uef/F00kBCh8tiTr9P47CU52lqMq1Pha38C7dNr1p+ifaiwKIUFnpVRwexDj2eOPzo8bfx2V+/3GH92xZEpJDuHKhFiKdyT5G5GsebqgVB1u917daWGHEAgAdnrYQPVA7D6l1uaYPbyD/Y3BzTRbhoUELnYOhRfVsmguN65Kbfay1W+skTC3DCT57FUsZJtDujxxOHm599D9MXJTKS1lzhKWE72ZXmfGcxZdWeiS3wnrNdL2W7+R3kfU9TddId+bStl21pLv/5eesqfKuFgzvzly/ggjtmJMp9m68k7Wu5Ijybx/V3HpxT0XMGgOnvB/HV1m4rL9R4OVCfxZseuqFaoMcTh45GZ9E53DdjmWEQ0Vd9P74mhUVJotkKb6lc5aUvKuXkH5vzQaJMv+dfPlN+UDUb0XaBuzc5ttVbWmLlbs7BbMoa75OUOuXqcQwch0VNc+fLS3Aj419UHuon3ukMBlgZcehg2K2V/ExZf/vCIoy5/BFc89e5xjpbmluxfKOZLf6vP89my2MK6SqujUqZomoQh1quvxufSm5I6rvesH0XHni9fIVpLOaVVaxU2V2m5RyK8YGVEFNIG5p0iYeMYiXHGMuN5eQ/guqjM0RHyIhDB8OucwjrWK5ftXlnyZHqdy8uNtY77ef/wEeuf8Z7XNK2eucu1c69etQhInzlLYJqmLL6IO3ohBB4ZgHvP6C29aPH05mI6vAO2V2h9Ent529z3dyQUR/toZDWiYMvV9IdTVk7wy05lxgRjSaiZ4hoPhHNJaJLw/LBRPQkEb0bfu6mXHMFES0kogVEdKpSfiQRvRX+diOFvCYRNRHRvWH5dCIaU/1b7Zzw4xwsJ0PPSaQrPl044ntPAgD+n4GjqBSd4GBkxT/eXVfWdfNXbcUFv5vB/qa+60K+Us9tP7GSN3HwUBLPXbEZQgis3dpiqGH2flbHaJrzevmMxelC25tQLf1hPedsR5qvS/icv9oAfF0IcRCAKQAuJqLxAC4H8JQQYhyAp8K/Ef42DcAEAFMB3ExE+bCtWwBcBGBc+G9qWH4hgI1CiP0A3ADg+ircW0Wol/WYbcKVYitZrq+17H3eyihwXjWfSRrvb65O7wa/gMK2Idv6liKftBuCjdirP/VuyBvr+cBfrFRu+3LyxZVOf3h1KY669u/sOBJ/Kw9/QSyGFA/dxWPx+rj5aLlipUrnbUcE/LupA5P8SDiJgxBilRDi9fD7VgDzAYwEcCaAO8NqdwI4K/x+JoB7hBAtQohFABYCmExEIwAMEEK8LIKZd5d2jWzrfgAnUUe8kQ6AlymrBwGpbAx+m0sHZitNoCNnR0ubPaSECepzfmFhMsMg4H9fvuI473qJ64JPlRHIUeBcZm/H3Z9pSAkHQM+57RIrdaJp642OzAAnkUpyG4p7DgcwHcDuQohVQEBAAAwPq40EoJq+LA/LRobf9fLYNUKINgCbAQxh+r+IiGYS0cy1a9fqP3dJeImVLNdXQ9nmQ6CA6pr3VqpwqwpRLPNcbbTsgp8OCQDmr+JDmfs+4UrShHLv0aQ7UJ8RF/guqZBO9qeHqjcRLN3RLDG3TaIvx7PwdRh0oeMFPfWFN3Egon4A/gzgMiGELUg/N7+Fpdx2TbxAiFuFEJOEEJOGDRvmGnKXgG3CRX4OFgJShYlvIzAx4lBNayVGasHUsvxSW6Jog+2kWqt0pzqqbc2S9HQO/laCv4KQJCzJ/Ts5rodnr9L64qGH6tbv0SyOsj8L35AkJnRFzqMa8CIORNSAgDD8rxDigbB4dSgqQvgpTTSWAxitXD4KwMqwfBRTHruGiAoABgJIeqbVEeVOiNeWbMCYyx/BW8s3e9W3y6g9OIcqEAdfOXk1F0lnCEde7gZr88725RxM8PUHcA39c5NHG+v5EPlIrBQ/HCRDdOtERf3R3nay3K7cNl3neo+VEoeeCh9rJQJwO4D5QoifKT89BOD88Pv5AB5UyqeFFkhjESieXw1FT1uJaErY5nnaNbKtcwA8LWoYDvSMX7yA2/7xfk3afiJM1v6PhX5iLy9TVkudSsVKQgi7TkP53txavUXmM2rrrXnetm0jbCzTHtaq5Lb+5h60LwH2Xx6+ugn9qqS+i1MD6lZwPuMyhs/Qim3ioItP2NfZnoQaLrwSLFq7vaKoAF0NPqvjGAD/DOBEIpoV/vs4gOsAnEJE7wI4JfwbQoi5AO4DMA/A4wAuFkJIDd5XANyGQEn9HoDHwvLbAQwhooUAvobQ8qkWaC8KvLViM77/yPzadBDOU1/5vG0xSa7AHpzPvjCEEHgyJFg6Hn1rFcZe8SjeX2sOKqaO796ZZll7WogUQQW50NLVODn4xGdi+yHCmMsfwWX3vJGsXwWzY5/63qShzAdl4hz0ef21++KBAr2IvqFcn+d6iBCVuDYohN3FGIwc1NtjVGZImvhff56NT938YkVtdSX4WCu9IIQgIcShQoiJ4b9HhRDrhRAnCSHGhZ8blGuuFULsK4Q4QAjxmFI+UwhxcPjbJZI7EEI0CyHOFULsJ4SYLISozbEelSXO+JcPj7H+vrW5Fe+FG6231YnlN0kcKuEc/vTacvzbXTPZ3x4Is6m9s9ocZbZWzjg+6UGtvyk/nnzQ7l59nn7ICK2N8m7uxdAH4i9MkDzfMZsQcxaz1HOLxMjYRprbVqv6HHh8RHU+PhWAnXNoVBIDOcVKxepxvLaDVHdDj/OQ9pUzc5v71IP3AADsObBX1F5RoC08unz+N9Px9/nBKd3f6sQ8nrZwUleic1ix0ez8tr0lYOh6WWzua+XGn5bojLn8EWN8/qPHDsYnDh3B/qZC9wkp99YeZzyFfdr0eZa+HOdMbwexZBl3oNCLSgpp5QcfU+Zy/VZsY+B+b4xxDi7umS9/5f31+Mj1T2PHrjbr9fo76ajMbPVGjyMOlex1owf3wchBvfGhfYeWyj5/2yvY78qAOXprRaSE9s13a+UKSpyDTaxkb7/NUmF7uCh6NZinQa00P1FUVn8xjBqfP+YhTH6OSnoOCL3nn/xtAc7/7avOdsqF7G9Xm99J1vZs3v5gq3WTarB4YHOiSFMu6K3N0cbpM6W9OAfDcUe/32Q4jei7KlZy9WkSvf7w0flYvnEnFqTMz+JKgNRdUB1NTRdCJdYxeSLkcvFJ/Mr7vFGVr1jJzjm4x+oSK9nakJtUrb2sOTw25wMM7N2A8XsOKOt6X6Klnvp0HYPexk0VRElVWjX+Ijepnba8zCleRbOlHflOufnuM69kDTVUBoGc4/N5L6bzin6pzZQ1Llay92daI/JAoUeidaF3Y2Xe7V0FGeeQAjkKOIJqilr8OAdznT+/ttz8I+xx/2W7HRG47Ft/mYP/+OMbWLhmm7GOr+d2UQivPVXPAeFzUPjrm/bkOw/OWhHT6fhMDW/OwasWDyl24cbDPVcff4U0JrCAefy+CmmbsYUtJ7VONE3NyCa+/PvXjP0A6LGODj2QcygfRME5tJy91KiE04q/du8sDO3fhP/++EGljd20iRWLAne/ssTab5uHjXetwwO3F4WRO9Edn8pBa7vw2rhcnEM50FNR+iiRbUTJJ7S1D2RgP64JTsek1/vS3a9h/Ig4V0eBo4MVFSmktWL9xK9epxL6BFHxJDK+ol8d3TEKLIceyDl4KqSZ40I+R96cg01eGivXVtsDb6zArc9PB2TdAAAgAElEQVS/H2vDdO3Nz7rFIF6iqTL3Z85PgDsVl53cxfab0qavk5NNevbbFxZ5jsoOq/mp4OuMG96v9L1aXuiFnOQckgMa2q+RGVu83kvvrcdt2jPxUkh7jM3Xmc32WtWoti5nOdN69SUOei3vLHxdHD2QOLh+N1fIkfQSdUOfd6Zr3lhqTgHYXjpp8pizwhbFJIDPRC73JNSnKXkC3cqkZLQ1X666Q22ytb0YW8DvreVFVSccMDz2t9w01mxpLuXEqBT2tK88sf/EoXvybSl3+eNzDjW2u2ZrUmZuVUiXubdR6T8zKlJIO9pS/1Q39gSHkWiHH0e5hDjjHHooSqc7ZgITBZyDz0k4GYOGv+Y7D5mzt/lYK7mwaafb7K7cye47LLsYpcxcBMpvbe0iZq300nt8xNOTx8f9IeR9f+G26bFym6LXBbtYSdaJ12ooRGNXn4funXzeh/bm+2Q6LVh0DtwG7vMqvZwGPRoym7LqnEP8b1WRr4opdQV3MqRHdcVKNgvAamFov0Z84ei9at6PDT2OODg5h/CTY2lLYqUy5kY5+6+Pn4OO5tZ2XPfY26UMbhu3u939yw3B4aursFVbusGcutQXrRp3pC552/qXXNVyzRfkF09HKT5181cXbPcq34leRxXPmcZLAL575sHxMkPdA3bvjxMPHM7+tnT9DlbJ7uWg567iGT7DdG3874RoVlkJMeLg4hxMOgfP3W+ZNkfrkcJzV1sxZq7bEeiBCmn7iw0mN69XkGIln800YTdehiq8XW56KdjiO15ajF899x76NOaRzxFeXeyOX1h2ZFf2VMpUszR/vSVdpm9uZF2s5Hs3Uh+jb0Jbdka2/SMG9UIa2MZ8xk0vYPF1pydq+JhGpjnk/uicQ0sbiz6e03/xj5jvQhoUhZtAzFyyEZPGDLbWMRpYWBTQQNznIm8TK+lEpgLO4Vt/eQsrNzfHyuqhc2htFzFz3Y5Axjnov4ef3ATIEaEhn2PFMF/XYszoVdIeNjbvaC1tXmkubQmD47W2F/Hjvy3wuqZssRJTxhHVWpy01CbbisWyzA0lZ7ZLYxNVbib10L3EKvFKvQoRcaiGPjpHkXBKH7+RMHiM20fcdt1jEbE3cRHenINW0KIEfsxbFdJ2YiHhcpzcvLMVv39laaK8HjqH1vaiVW9UD/Q44uDC6i3BKaGdkR3liNBUyLGZwP78etzfoNIN8bDvPuHUOdjmtpcIINwRTDLUNVub2XIJ7h65hVP2k7BcqHZtEyuZsO+wvkZLrufeiSLqpn2NevXRg5NB3/Q2mxpUsZJJ5+A/BqKovm38JxwQ5UTx4WxH7lZZALsvfXQf66D00mRiOEWsFDNl1eo5xFMSqsRQJ3yPz/kAh13zBHtdtZIHmVAsCrQVRYeLlXoccVBf65jLH0n8/n9h7mBuAuQoWMgtHqGr9dOMiVgcNnqQsQ1poplmKroW+WtLkmImfXOVmHztU/a+mMu4hShl7bVCW3sxdZa6Qi4Xie2qCP2ZXHbS/s5rVM4h1pbjXZruOeAczH4OQEA8fnfBZOfYqonJobjJNKY0+RxsTnAJnYOHWOnI7z0ZM4l+fak5dlWtFdIyUGBDPodPHT7SUbt26HnEwTPcxNihfRO/BZxD3stxy1estO+wZD8SUtxRFhNiOGp++paXE2U+jnIcfEVIF//h9bLat922auVlc4IzledzVNVFHiVm8lHIxuuosmUTiTMRApMnc8Q58ONJxJmqvaSktKGbdFzqGAb2bmBNVD83eTTuvWgK+veK1KUuXcX/Tk+KhoA457B9Vztufua9WP8m1NpYSR7WGvM57D2kD4DKLBbLRc8jDo7f5SmE8+jN5cxiJR2+KQ5tA5LJdZZu2OGMHFkJymWTuas4zuG1JX4RRNNATTKj6wx8RDANefKLMZTSIkuvzo1Fr2OKiluuWah6IjY1oStjfe6y0v0pivfEQ10zTYVcjIjMWLwBu9qK6NtYwNH7DMGhowbhspPHJa4ztc8T0fgzWL898hexEYdab9PSkbQhH3GAHeFa0fOIg+Mhl5TAhnpNhZxXRrQk5+AnZzVh0464Seq2lja8zNjz+yxg/eRmEiu5wN1TNRXSvpe1adZKPsjnyMvqxHfkQvv0qXvd2Yfg8cuOxQglBLx6I2pbLJEx8Ck5chPItCa6+njKQYk4mBTSyveGfC42b879VcDxquHSpfOgTzQC7tBiewQDbJxDjU/xUrzVUMiVxtgRnEOPM2V1TXEZvtpUK5/jrZV0JHUOhtF4vnRdOfW1e2dho0Ywjtgr0l/Yln5RCOSUGuWKlfhF6FevXGxpbsWAXvGF21b0i62kopDPeYmVfMdeFAL5lNn/ejXkceAeWvwi5ftzC+ypZoXwOBEbxp+MM+XHRfmERjdBhvQwba7qGinkiV0zKsdTElMlOAfugOIeX0ynYbnPWu/TEeeQK83rjHOoA1wvdnDfxrCe2cKhrDy5Hqcl3dkmXi/ewAIme9ugPo1epztdlvvTJ9/xuIobE9N2FWcxt8jfY6K46rGVXMrpHAWLv6pjLYmV/LkRdf8Z2q8pUc9HV6PncJbtyk3cpANJJD5y9pQeepvyfOPjIZ03xDBTOR75NUHjmfa5tv42l0+fGzRhfiK1PsXL+dxUyDnfYy3hJA5E9FsiWkNEc5Syq4lohZZTWv52BREtJKIFRHSqUn4kEb0V/nYjhXdNRE1EdG9YPp2IxlT3FuNwPWLXew8C77n7SYiVPF6umizINS5unOqk/Z+n3k1WkGOrgkKNqH4+DfF+kxv/pL0HGwmCXn73hZPx3DdPCPNyuPvzXZQ+aU9LddlKQaHpZM6VCwAf/fGziXKbn4NEOQrpSl+t5H71w8myDTuwvaUNQgTWe3OvORW5HE8c8jHiQGx7vM4h/jeX4Ed917Y1XuttWop5Vc6hA6RKXpzDHQCmMuU3qDmlAYCIxgOYBmBCeM3NRCS1bbcAuAjAuPCfbPNCABuFEPsBuAHA9WXeixecTnBCYNmGHfj+I/PZ33O5aAO0RQN95u01sb9XbeZ9BtKILWLXlck6q22Z+j501ED0tqQOBYLNxVe2K5H2xMVV17fIW75wBP7zlP29xUrHjhuG0YP7eEfX9R1yqZ5Wn/fCTRICeb3Lyijep5m7LbVjuLbJYD5r7a/CbVFaZelizGN/9Az+6fbpEBBoyBH6NhWQI/4QU1DEqybrJ+6x3P/astjfn7zphUSdp+ZHa9YngGKtIA1eGvK5UmgV3eiiHnASByHE8wDcMRgCnAngHiFEixBiEYCFACYT0QgAA4QQL4vgqd8F4CzlmjvD7/cDOIlMx6cqwBk+A3brGlI4h59ZxDG62OeV9/lgcOpouDk3IDTZ8zGNFaYfNLjCf/hkhouUi/G2bG3XQm46ZmjfsjLZkScHaKqid1kSK2lXnHbIHthPCcmt1lXb0PvRx8aSGJM+QfVzMNTRQ3b4bPyV7olyo+MMIN5YugnFYkTUTMS7wHAOPlz6VQ/GA1xy5ujqAc52r7WW/+8I/YL6NubRrylY/9tbametaEIlOodLiGh2KHbaLSwbCUAl0cvDspHhd708do0Qog3AZgBDKhiXFW7OIS4LnvfdU2O/qzqH9w2hoTlIy5gffToeetl1QpQnJZ/TkRACf5+/JvmDBlcspTyRc7OQCjt9HLa2ayGrle8q7XGCPMcjOaFnF8Sfq84RmLixpkIeP/jUIVpdOYbkoCNxUPnPSvVzMCGRN7wCsdLZno5ahVzO6l8iECm8TeJbNY+DDJzniq1UDmzcQa11DtJptHdjHn1D4rCtzHhYlaBc4nALgH0BTASwCsBPw3ITD23mre2/xUBEFxHRTCKauXat3ZLDBCdxQHzh92mMG3SpJxqTfToHec0howYm+rMhOh25J2RRCMxb5c7x4Dr5BPJedx1uXDaxUtoTF1dd3/SiDdYkq+fbzlHUvjydcWgNT5gztACGCWsf+ckMOqn8lWIlpUzERU16M7OXJ/N+mKZETFxlmGEmr+xy4Mvo53LByd9kQqwG9tNztUs0KKFU5QFFP5DoUXZtGN4/aQgA2PeJWsv/t4c+TX2bCqW5ua2rcA5CiNVCiHYhRBHAbwBIP/zlAEYrVUcBWBmWj2LKY9cQUQHAQBjEWEKIW4UQk4QQk4YNG8ZVcY/dJVYSwhqtUT3RNKWImiivSdiXC/Urp4CLX2+D76R1WenkPTIamRya7GKlyleVftrmXpVXylCFyH94XzOjKmW9eoa7pFgp5ByYNvQAaiX9glqmXaM/Kk5nZZrLBDgV0ibiZoOpji/XJkQgR1fFSioBWPDBVqdYiVNI6/U+8+tkFAATuEgIXJsqam05JP2oehXyRnPdeqAs4hDqECQ+BUBaMj0EYFpogTQWgeL5VSHEKgBbiWhKqE84D8CDyjXnh9/PAfC0qCHf5tPyhh3mBDmqlU4apZ68xidhioq8YQH4OqDZxmLsM+cWKxVKOget7SrqzdJMg9RiJSUvh60XaXSgHxh0O/hSIh9lzNLLNsE5SOJATJn8WxuVlJHPuPLkxDU6SKEOapWn5kemmyadiQ1CCJY/M03pxNwQAtta2vDbFxeVnpN6Ttm8s7X0nIkIHIMRFyuF1koVKAFMV1p1DjXWDcugn/k8daifg9MJjoj+COB4AEOJaDmA7wA4nogmIni2iwF8CQCEEHOJ6D4A8wC0AbhYCCFjTXwFgeVTbwCPhf8A4HYAdxPRQgQcw7Rq3Fi5WLu1Bb8OczhL3HvRFOw5KIhIGWSCC8oTclsL5GJImBAq05ObkDmD4pebK9XiHHzESqYTjY1zWL/dnZUuLTihkg+dyJFZT6Ci9Bx0cZYhD6ysft3Zh5QS7ugOjNE7j9qQ7yRnILqSOAxTxCC2ja3EYSkNXXjnzGj8ZQQHN/Xnm1FNnXctbUX0asgn5o9sKk/84UB97q5T9VWfGI+Wtnb86HFz6Hpz5AIb51AbrNvWgn5NhZIjaT7MPAnUXpTFwUkchBCfY4pvt9S/FsC1TPlMAAcz5c0AznWNo1pwPeS125L5eI/eJxI7+G4qOuS6SMS0cbRRKJ2O4uW8Qtp3LA7iYFiY8XHxhNGmkFZj/fvA53bKtWsjUp+XuSe5AbnEWfqcOGTUQAwf0Cu8Ng6Oc5Aho00Hjs9OGp0o22jgcIXStq8oyEtU4imiMvWhBhjc1tKGXg35xEFFPmeTWEltUq4NU/iXgb0b0NJWnlrVdjiqlYhn0vf/jhwBu4fzJpeL7rfLiJW6MiqVF6qT1jaB1KiRQd2QczDIn01II3P0Fis5xCmBtZIduZIuJF5TzYWgo9wwHSoSmxqz0fogp1hk2R7bQSP6h/Xj5WYP4+BbLBe0oW21BRnTSyqK9TF9ZL+hieu/dPdrpmG7dQ4pDiml+D7gw5T4xGm664uTMWq3PqW/563cwvYrCYhPOt6mQuAHYEpgFEjX7GMz3bddIV39jVoeqooi0i/liSIDhQ7gHHoecXA8ZNdkUu3jbZuxLkowcg7K1vEff3wj0U6a+DG+88fl55AzOLgB0eKVnINe746XFhvbfWzOB8bfgraixtrai6kC46nvzV8hHW+DHxPfpilfB1c/qS+KxqBDPl/9/XKM2tqtSS5XwmVBZJCKsRu9ukFx88LERapcwXH7xw1IlqzfHtTRGpSck6rbU+eFWpuIMKB3AVua+TzpPjmiTe/ebsqaLHt8zipncqxdbUXjAYkjcPlcpHPIQnbXAa5H7NpYfGMr6ZMgCFpmPvma+zMofpnrqqaQ5jatUmAcOa5kW+8y8Z4kvJ6ZclOHf/dJ/NtdMy21Zbt+ZTpMm48OM1Hn+ywRK8s8kv2qdQ7coz/bnoRpA3bB7EXND7Bfr6SkWSacMT0lUzpLW/TiptAMXJ+L0sjDl0Pv21QwOogJ4ZYUmNaC/twm7b1b6bs+nubWdnz596/jC7+Zbu3rQz98Cif/7Dn2tx2tyXsgNQyKteXaoOcRB8fO4Tp05jw5By6MsL4g9Q3B1B/XF5uO0zCccZqHrssJjtsgjts/LtbgTFnXWE6yaU1xtxoWvOmZp03coxoWmIY2ZZ/BCf+DEoT+p8Y5KDNp/+Haxh9+qk3ec9EUjBzU2zimNF7g/RoLzk0leUgJxZ5aP/sM7Yuvf2z/sI6pLX5strwncg7q3KHkHPKKUYTNgIIs42ovCudB4d3VvCOrlZsEvxZd/hXrt+/C4vV8cE2OS87nyOhPVA/0OOLggjrPrzjtwMTvqkLalkqwlRE75ChpL+965VJHoS8QLq6TifD11Zy8XPkb9lTzC5QQ3wA4jsbntGyDXxY1/hr1nnzESnHOga+TV8RrepMmURE3hlyOcPYRIzEytHgrERyl1UF9GjFhzwHKaTneoG/+hVeuOAkD+zQ4A7aZlOQ6EerTlI8darhnaxqaLZ2uFCfd/MzCWLnkHNT3Y5s7ZNGPBZwDj90HBFZfU/bhfVy4A9RDlxwTjkfrpzQW4zCd4AhgntwBFGuJHkccnGIlZdlMCnPexn4PT5xCCMxZYfZG1uWLRRFcO6hPY6ItG0zxY1jiYGhDn7SuIF4+ic0l0VIJks2ksVpB7hImveGftiCIHHw4B1Vp7S9W4ltjFdQJ81ilHa0ZX/+YPQZKC6k4Z6dvPi8bYn3pIsV4Vjn+3kxcjW2eyc1X9+hvUhXSHhx6wDnwv7eblCQhTjxwOG76/OH8+BgCIM3Z9TYjYl8+2Jz1ObPHfD3Q84iD4ynHTnzsKcnfeuD9dduVfgPOoV9TASMH9UZjIQef175PmGNaXwAcG+qTSB0AWlrtaU5tXsdyg8gzCulmS7s+z8vPIotv10d5rUI9mZoQ4xwcuqKEQprpTzjqqMSo4t1AU2TqHKtJH5DX9AekjDOtWMkGSaz0TVHqIlTdntUvh8yPqr0ojOJMIYLQGcYUrXxXAJJzUEo0bc/B5Vtk5BwY3V69kGWC039VfuZOwpwi1oRpt76C1686pVRftvfcN48HAJz5yxfZBffoV4/FqMG9seCDrWhpLeKB11ckJiR30jCJ3XUiJx2qTPfATXK9JM88B5tpZTkbPwfjKVG5OG65FHzXQ5DHOAdDm2rCmYQYRptH8i9pQaQ/QoI680RsbOqYIguqyjYDvf/lG82JpABFrKRd2C4E622tooyguKX5q2+KKucgf1Pntf6uyDKwohAWJzf7Zv69h+fF/j589CCz/s+Dc9CJsw5OZ5bPRU5wHcE6ZJyDBvXFs8ShpCAyt7HX4MCeW5UTF0XUXiGfQyFM5ME1M2pwbwzo1YCjxgyOIk967JxGsZI2beUJ39SmlXMILymZsip1uDDIEn4JksrgHOCWSwNAk+ZcpuqOAKBPY/IE2ZDPQSAIuvfw7FX2cYRtfbe0qWgncEVk9Pw765gaZj3ItKOSDnAu6G2ffmMyfwEHXUQ0Z8WWmLc1twHq66S5tR2zliUDBaqQ96lb9cmTvOql79Y5mA8Mpiul9aAv/t9pB1osBw2soAIXceDWomqtlCmk6wDXI1bfkW2TVF/WYVqk1TMnBonPP31kFGuwyExGk0+FSlSi0NjmkQ/t14iPjd/dWOeCY8bE/pab+BKD5QQfSjp+ajI5wZngpXPwUBuYdA4u4pPYiEGxjZizHCvkA+7i3F+9nEzLahBvldrXOlRFRjJLn76pqjootblxu7ut2kyQ7bn0TJG4kHn3yuHV523/9wNv4axfvljKRcJBDkf3c4g4B0Ws5NQ58L8Vhfk3hSHSyvkLGvKRu7I+l13WfwDQ0m4X5XKSACCdGLva6HnEwfGQ1cnBW2YkX9aDl3wk0QcXiZPjRLjJGIs86cGptBdFTEyi41jNAamlrYg1W8wOO5yYICKKwWeD5vfgastn469E9BR7jh4nQjVNqJpHQEWgc+A7dFkrJcQflBx7koCoeglFCFXGzuDKPaxHI5VdcHM0rcXMm2F48S2WHATy+Y0Y2DtWHldIS7GSjXOwEAcb5wD+Xk2bdDAmvtwlVmpubcfF/xvPB/7w7JW46ekola+Li884hzrAtdCcYiWPl9XaXgQhvlnf8dJibN4Z9+Q0iZVUhyef/rY2t8XSl+rQ76Klrd1KbHzYbekU55qyg/sGJoOVhv/Y16CYjzgHe/v65q96uptOkQQzMTKFR5HY1aZfmCTenELax2vbB64N/a4vTmbLufqkcK/6mI9UnMOAZGyoCXsOwJ++/KFEm3Iz3H/3uA9OpJBWxUrsUIOxwSJWEgJD+zWyv3GcvCw39mXQObgYhwUfbMWMxXGz90v+8AZ+8kSUSdJkUNGBKoceSByYsr9ddhzmXhNkfFNfNH+K4ieIiubW9nDjt598TCyxekIxTUgVew7qHROTJNsj7KOcFFtai4lT0OjB0QlOFytNHjs4sZAk52Ab16S9d8OXP7pPop5pwcoa85mERTKtpVnnkGxv5652PB/GekpuxHHrIV6EaN54vvaxA9ixS+hiHOI0p7qYkaJNU32svofGq88Yr/Vn3lRGD+4T+1vWswW749oyhSOXOGz0IBzFmIS3MfcJRIcOomjdWMVKFs6hvShw+iEj0L9XgfVA58VKxq4iLlirI8dpUnC7wtXIsXKI0r1mnEPNwT3jfr0KJUexOOeQrKuLVzg0txZjk9aoqDU48KiTzCcT3L1fmmLkQuSYP6QktGlpK8ZWxtFjB+Ph/zg2Vj92PZIEo6EQN2XlcmR/bvJeJTGB+rwuPXl/dpzyHk/7n38k78FAlG2cw/wPzH4oqthCQPB6FsvG07/JzjnofhfcQUDvs5CLj0nCx3Lp+AOG4V+OGZto23dPsSUrUufDQ2+ujP1WyMXvwrTZ61ADzcXqF4LW8jk/sZJpzMFYBAr5HD6875DYuJZt2IHNO1vZ6+ze2LyIlwuHEvtdu+CKB2Yn6sjsbzqkECHTOdQB3EJTFcAx0TUrVnJT8kF9GkKxUlDHFEaAHO0AairEYOK+tXwzXlsST5Q3YmBvp0Pdd86YgL9/7bjSeNQlnSMqJaYBDCcgrahR4RzueHERpt36CtuvbOt5JVrrTsNCsI1fDunHf+Nj87NZwyyevarISAiwx0hdNKji9ENHxP7W6+liJ45462PK56IsaWk5h4TCXSlYucmdNtPGOajj0DdGnXNIxkoyBOVjiCAANOY5sVJ5O6Okz7ro6YybAsuttz9IxgKTfU0eyznA8uORBMUkjdX1GH98dVmiziLFJyrWp4Eg1QM9jjhwx4VSzH5ycw6RDoBv/mefOQz/ecr+odgiKDM5HOmbQ5/GPP71I2PZOkUhcMuzC3HGTS/gX+9MBqSzBQQkCqJ97je8PxryhJa2Ymyx6AvcttFIyBOeEMDVf52XrFAaV1Dvh4/NL5Wt2cLHYPKR9761YnOsvLSpMY/YFo8oIKbRaZmAhHWN6ZmOGdKH8SIP6n36iMBC7cA9BmhtJZXb+uga8lTKAhYjDsa7UNoyHFsFRCzJjxNMZzbR5vgR8fvUazQaiEPRoFyR9dW1aD3NWw5FJUWxshYXrduOTTta42NQxxUWcc5xUf6I+GSLOAf+HfhYM5mitarh0uuNHkccuEcsX3pwmlSJA8M5OAJhnX3EKPRqyMcUnjbPYRVt7QEbHOuPopPD/FXBSWfjjmSIYoJZtqlyCfkcBVYcMQ5Jq8/M8S2aMl1ujq5DnZzcqjKe8ykA7KejRkPkzyjdZHSxrBl/f/HrY9ZdIrhnfcPnLIwAfuOQ9Qo5wh4DkrGpOOW23nQ+R6VTpioSG9KX19Ho7XMQIvnuTPUAfl6bZO0A8I1T47oX/XqTWMnEqUgrP/X9uMJnmMhnpAuIaqxXknlxYbLlNdyhoJDPoXdDPhEFVr5X0zvw0znw5T5i7Fqh5xEHK+cQT4/JKqTDsrkr4/Js7vQtqb1J56DLoduKxUSAtUjmKKwn4RyR0+Ih6JNCz9GojPMejo+TcPjoQQCi03WJODhONKWcwMrS+dfj9sFNnz8cl508LlbXdsI66aDdS98fnLWi9F06HHIn7fjzircdy+gX6hwSKSsNljC8SXJ4D0IYdFVJzkE/ERdy0Tu84HczAADHjhta4kZssBF4P7GMUP6Pw5ZwqiF05iy1YuAEJH71T0cEbZkU0oqfg3w+21uiw1VCb2NTSCsnevns1TnRyrCbrmB//XoVEkTF5aBqM4/V+9WhWorVG07iQES/JaI1RDRHKRtMRE8S0bvh527Kb1cQ0UIiWkBEpyrlRxLRW+FvN1J410TURET3huXTiWhMdW8xDm6x5xXOwe3nEHzeO2NprHz6f5+Ep77+0aggJlYKJvdt502KXaNawxTDODB6pji5EbUbiMPvLjiq1JYp+Jy6mUlxl+qx+cOzDwEA/OicQ/HNUw9g7/urJ43Di5efiD/82xScfNDuOCC0/rDNWQFeTttUyOETh+6JEVr0V1tbn1Bk/JfeM6v0fUi/wFSWO53Z4geqHrjSWklfw6ovhAqbCaQMsMhBb+ogTSRTyOcSG81H9x/mHXQvNkbFysVnX7HJ903BHyX2GhJYwsmos6p/g04cph48Av2aCgrnEG9LchpSIb184w6c+csXjeO2GWKooU/kbcXyYoRlqhNr6TkYTvL9mwrYluAcIg6FHUclxEEOtZNyDncAmKqVXQ7gKSHEOABPhX+DiMYDmAZgQnjNzUQkj6W3ALgIwLjwn2zzQgAbhRD7AbgBwPXl3owP7JyDZsrKLEq5UPTTw9B+Tdh3WL9EPSAiDr01cYo6aeXpIsE5KAuTIw4nHDC8NHbTCUW9Sir6vvNQidaXNtjPTBqNi0/Yj9U5FPI5jBzUGwePHIjbzp8UU0jbIMevjs0UvdXWlilAmsSEPaONlphvSUshlXNI6puCHBhk3BB1yEuFEGwGMiKUNqPdBzRh2lGjMVgTFxVylIixs24bnyea6SHZXzguH3k1Z0Krt6U+H0kIAOCMQ0fgT1/+EM6dFO5VPvcAACAASURBVHA4W5XMbFyE3/iz18RQJZ1D8OzfXxtX1CYNC8yOirJYDeutvhtuvrk4h8ZCLqlzcDh4lsM5yINDyQCmM+ochBDPA9igFZ8J4M7w+50AzlLK7xFCtAghFgFYCGAyEY0AMEAI8bII3uRd2jWyrfsBnESmo1cVwD1ieZrQRQvcIOT+bMpbW7pW2WykjkB3EFJZYrkpJHUO4biFsMb0X7hmW0nRxo1FH9cLC9fZB++Ar3MOu1GGn/2aGmLltrZcxOHK0w9i2uG+hWNQNozAuYvws89MLP1+4IgB4bNPjuqsiUFmtH/58JiSVYsa48nkH6OelrkpLqOI/lUxF+U2Pi7PCGeNJUfvQ+BKUVK5pDMli7notwf+/cO4+8LJYd+Eo8YMjlnWSXAKaXXe28RK3L2zYiXDPRUYiYBL0e/K8VHIUyIfihRPmQi5D+egM/2/Lz1b2YaziaqjXJ3D7kKIVQAQfg4Py0cCUO20lodlI8PvennsGiFEG4DNAPgMHFUAN+FKey7F5Yc2nYMtyFzYVGmCffn3QbRSmcgkqhO17+Yc7OIBW6AzdSOS8tdywgRwbbpkodwzlBzQSQcNj5Wbomi+ePmJbCrKTx62Z+m7/myDsdnGFbftJwJOGb87fv7ZiaXfOVHTgu9Pxb8eOxYAcPUnJ5SsyyKdg5kzkuIIGb5dhzwlf/+RyPqLk2cfwMSBSnJ7Ktfkd3L9yxsr8MGWZuRzhHsumlL6jRMrDevXhGPHxcOycBwGp5DOKWFJ9LFJTkP1Q7FB19tt2L4LQ0NO+OIT9iuNKyLMyoGBaZ4Tr6ni4IZ8knNodUZcNd8HZ1ABxD3FAT+LtWqj2gppblUIS7ntmmTjRBcR0Uwimrl27VquihNcw3Ih6RPNls+h3UHKOe9anXMIxhM/sZmIQ3vRnnbxW8rJ2YZcuFBs6053zrLpXkztHLhHf5x80HD2hCzLdG5ACJGIfgoEIgwur/UYLT6QDtuBTYrXlqzfjjeWblKuEaXfOZFFUyHPOimqCmnueT33zhoAwJotzSHnkKyzLAyrLTc30z3wz5S5yXBctnf9tVMCh8T2osC3/hKIGo/ca7eYPkS27bK64UxeOdFsIFYKzDf1jImNWshuvUcpulIHJ+tMf389jvjek1i3rQWfP3qvkmMrJ76NjVN5eJFzXlRPPcQ05OLEYc6KzfjOQ3MT96jCRuRK4fO155AM0dEJxUoGrA5FRQg/14TlywGo8YVHAVgZlo9iymPXEFEBwEAkxVgAACHErUKISUKIScOGDeOquGHbFDW5M+8EF3zKSXbIyIGJOrLesg07cckfooBbMs5Q1EFSrJRnzCmBYFy2yJousUs0LveJzEeo53LOefyy4zCoT2OCwH78kD2MbRYFsG4b7wPBRgt1jDHmZazds7RWOvXnzwdtaadjCjkHzmw4NgbttCwEPy4pitq+q90oerru7EMBIKa7soWzsI4r/BSMtPo/FQ/1w0IrtKIQpTk4ca9BmiiSEmIebo5wZpemA5aAwCFXP4E/TI8bdjQqnIN+63dfOBlnThwZK1NFRm8o3LPKseiGHxKl9kWSYKhdq/tAQyEuVvrMr19mnelU2KyZPv+bwHlUJ7x9lNDliQHVCeUSh4cAnB9+Px/Ag0r5tNACaSwCxfOroehpKxFNCfUJ52nXyLbOAfC06Ai7LSR1DjbOQZ70zZYkhOfeWRs7Ce/Wp0GrEb1z2V6D1p7cFDmrE9WLk9tsHv6Pj+AbH4uHqtBFJZxMWL8lU2gJwN+UVaK/pmdQURQCf359OftbmtNyqT01SQxzsRBJB0WVc/jza/xYuDGoYgvuXew1JDC5lUnvuTqHjBqIkYN6xzYKljhwG7NFIa23caliQiw5svZi9LxiSWZC5Ihip25b9AAX5Bzcyfj/RGKl5Lg57lEtUjdh1TPbxDnIbWYPxWrOFeyvIZ/T0v+6tyqbWOn1pZuwcfsu/O7FxbFyqXuUt9cpOQci+iOAlwEcQETLiehCANcBOIWI3gVwSvg3hBBzAdwHYB6AxwFcLISQM+ArAG5DoKR+D8BjYfntAIYQ0UIAX0No+VQrqJvZseOGYta3Tyn9nbBWskzG0kQ0maCxJytmAYen5atD1lQ/Iavy3o+MGxr7bZgifuBo1MEjB+KSE+O+BLpH6SWhXDZWx+Ns6hIrmcbFKaglikVY83LrcIWFsC0ofVzSokyUiENwyndB172YNn6VyJtETxKqDJvnHNyEUs093KIQQD0mlJqXo8S9hpyCiraiwPKN9uftRxqSc1Bin6F9S89JmhrHuRWeexQiECmpivyYZaDCoRdjYqXg88fnHlYyZ31oVtCG6XxayOWwq10lMOb7LPXjUEhf8cBbsb9PVsRYvoYftYAzTagQ4nOGn04y1L8WwLVM+UwABzPlzQDOdY2jWlBfZlMhj0F9InNClUUFTCG7gzIpdzS9ND/WnyBQxLcfnIMn5q0GkDT9U0N265vCNWdOiOp52sLbwmwoA4v/yTadlC97NWZ5MlxbF5+wr7H+hu28yMfk76FCf7c7Q0IgNyDOBJNDdLKTn/zGr/qrmAgIEGzWqviQuxWjqSwHIWIn9B+fexg7LulnE5TxBOivWtC9xLh8OQfwc3DqwZHIURfxAqZERIHI6LNabK/p70eSafVeVK5MHhQH9GrA+D0H4M3lm/H76Utw6cnjLKasFOMcfIiDy5RVD7qnmjj7BN6sFXq0h3RCfEJxu3Zink6COBjeme9CAeLmgwlP69jijeqNHNQ7prj07U/XOXDj92nLxjl8RlEaJjgHS9Ppx8I//Cv/L1Cs2u5TH8eOcIGeNXEkLjhmDP7r1KS5KIdofFKRaT9UtIfv0fQcckQxB0U2GRR7gua5CdeWIudbuxAlXxFd5+AL32tMiamEVke30OPotSoyUrFJCRlCyoGojeEcAJT0CFKsa7I3KWgK6UrzlAD2gwinx6kXeh5xUL5zIS+csZXCIrmhm2TuXkpdSoohdJPNKHSBfey+TrSE+ETjxu/F9WhWOhKvXHESfnROdDpNpsI0t8mOxXKBa12qvyczs8XbfS80I20s5PCdMyZgYB+zbiTeTvAZSRl5J7hIrGTWSwDB82qNcQ4McWDjc3gVJZ6/SrSk2PLEA3dPdbiJ2k6jc7C/vHx4UHM5T8o1pEP1BQmltwCA15SkO+p1khso5HN45f31mMfkFAGkKatdJ6RDPfzd96Vk4iNbWBzpg+XjK1FtOMVK3Q1WsVFC55C8Xi6AViZ6ptaUE/Jkpw4jn7OIlUT82ng934UZP7XxppLudkzRIvfQQmLoG6VtnNxYuFOyhGu52KO8Oi72BGfKajtUSEc300ZKQEymzZmPpgmn4dq71LhJ7e0CfUNZfTnPx3dYgbVSEjpXXxQiFuTO7FyYbO04JTWuSkDeUXKBq3uBNKkt5MkYfh4IxEoq8VZ7NjmpqvOQ89exWTPJ9+PjZV1t9GjOgaENcQUYs8WrixywEAeP1SXt6NVJr8dWUk0sb3n2vdhYY+Py1Tnk3DoHX65HjsteT+dw0ukc7GIoe+fC8N01jjSIdA6RWIm15FE2YZMTHBA817hYKVmHFysx9XLk9E1QrZXailH8rrI4hxR1XVZY8hBz2b1RHC1fjkmHKko9UPHfUEdw7pGBKPSjWr71By8+Jva3LlaKcR9F3olT3fy5ubFTM3w47ZAojpgkOC6/qlqgB3IO0fekyEMLn8FMvEjnUJ0NVucc2rV21fAZqzY3x8bK1XP2qUcfZTdknkDp7chx2ZBmk2G91y035jpL2cZWRiw7XPlxxtGw9H6CT5M+QT6H2cs3o7XdbK0U6ByizcJXrMS9o8Z8LkZogCiKbam/XNRPezEKGW97Pr/7l6P4sac4oOgvb48BvfDvx0fGBza9jQ6nfYWie4n7OUTfjwhzYasxo4DID0QiMGU1d1gUgM4cxEVjyWtUg4HF150e+y3jHOoK84siiiui+HwOptaSbflCXdirtzbHflNlwrFr9HF5K6TtJ2pf+JrYcc/YBFaslGIX/59pE2N/2zaNtCfjyWMG49+O28fYjhRtCGE4xYf1pDetTeegWit9/ui9km0xz4QjhE0N+YRXvR4JVhUrqZyDjfM16WPkFR/axx79Rj+g7DW4D17575PQv1fULktgTQppa29BLdld3IckqpE3rDMdDXmyOqNylnIqQeKkEUeH/kpnHz4y8VuJc3AcRmuBHkcc1DV0pJb4nKArpJPX64vGmH3Nk8kOFNLR37oHqKqQtiGNvFdti2fvPURiihWVqz8VqoWVDq4tU4YsILn5q5FZAe2ZaXVTx3Y0iYG0sQQ5OTwcC006BwJa24LGbj9/EqYwG60vwWwqJDmHRFvKptjO5BPhYCJsUViU6P55CzRgm5KjgScEyUKjE5wH51CyJlMmxdWfHB+1rVhtSXCcIhdbSYUuIgKgOQ8mr5Gh66/6xPjEb3JcT85fbeyzVuh5YqXw89f/fCQ+Nn732G+6sjYNa5us564jI4OqVXtpHsuRbD++AvREJem8U+2rad3WeAgLXkwSfG5VFIZ6CGrZn4ovaafv5795AqYvWo9v3j+bJYArNzcnC0Mk74OMv9tr8hjevwlrwmdhND1VrJCAwDKF9TrXGjBtwmpeDj1Cb6kt5l1zb7SxkHMGiMwZOAfrNRZ9CeAWgRARFq3bFvs72YdZbxNrK/QVUtHXEhq/vSjQpzGPed+NZyHQDQsAYDdmPjfkcyiKyCRZx47Wduymlcl6Jx+0O0sc1JAtOuSzmb18c/LHGqPHcg5jhvRNeixDj62UvF6fn5UppIMG1IVg8pDWORT9hJLOWkmVuybrqDFq9hnatxTzJz72oD/V+/bOCyYn6unj0je8vYb0wehSNrfkYPSFboP+CGymrD7EO2a8YCAnslTOG9MGqz8H0yaco0ghrYdSUev4oKmQi70fDnHOwR4WXsL0LPTQMsbrKZ50x2edqe3rbenT5rFLj0vWCb8HOgHu/QSfqlhp565kWH5pMNJWLOLf7krm5uauaSsKNOQJt50/iX12sscaZiooCz2POISvgqfScbbSpOTzgc9rjhTSpJTxxEE/jO3QiYMyrCtOOxB3fjG5UQftxRfToaOSgQN3KNzANWdOSJinynYAYOmGKBnLgN5JRtTPoU7K7ZNIo4dT7+vtD7ZYFdLNjhO13p7pNnTdS1uxyDo16UWmTThHhBYH5+ALlwgkGFdEHHTCts/Qvvj+WUFQg/+aGuWKNivTg09n9FbEZfNcc2wIEhNxUP4+cI/+pThWEoVcJF4zebCrz0GCy49Set8CeHZBFBlamqjq6xIIRFmldLlM33KedjLa0AOJg2ThmN/yOXMeZokmPcy0SeDp8aLzFIQltk0K7kQDRLGASt0pjUwaMzhhkhe1FygDJ44ehH2G9cXUg0ck6qhKUKOfQVj8g0ffLhX5BPHjxxR8cmy6TUGYrB4VnHnTi6VNlsNf3lhh/I1rz/yOgh9+8/z7AIJTM7fx60TfxDkQQREr8XW4JyIzAup9qBu1mv9CQhUrtbeL2In+6W8cj3+asncwFmW8Rp0DkhsshxxRTCxqi2EWHyvfp3oI0NcnEFggbWtpw+YdrWg3cHZElBC56ulA1bHqc1UeGjmRWpvCkdnyj1fLvLpa6HnEIfzk3oNuw8yhSdsATYcknxe9aP12zF25BRu3m1NBykmrp4/U56C6ibtEA4HnaRFjh/D5EGLRXg1tcYpOLuFOGuX20g07Er+NV6xrTp0Q1xHphFl9Fy1tRVzwuxlK3Tj0RXz8AUliGuMcDNRePgeZWa+tKNhNXSeyNp2D7NdUh5tz5xw5KlGWo/hhZ9RuvZk6waf0RjYRrbxDDKSWu4gDUVz0xIuV/PQQOueg6+wAoH+vgKPdvLMVd7+yJBGJVyJPcSuqkzWdpKwDJNef5PK4e28vitI64nUOIefAjqrj0POIQ+nlJ18FlwJQh346vuGzE9l6Pi9aZv566u011npBSAV/qyCbUjHYMIrYtKPVeDJVxSKmtrg0qTrhBIABvdw2D7KLGYviaTx+8KlDYsld9lHyHADJTTKNCEpfxNzJO6ZzMDxS3VS0rd1kreSvc5Dg2kmOzIx8jrytldZtbcHfLRYx6lRxWSt5cQ7KfLZFsVXBhfgG4vNA91MAogPOvTOXWtsJEgwB+w7ri9MPHYGjNGtGwGwgIsVK3L0/s2BNab1wT05ekXEOnQQ850CpOIfJYwbjYEOyn2q+Z1VJaasTfbcQhxzwzIK1WL5xp1GmrRIHU1vbGZabEysN75/UVyTGFPaxRSM44/ccENfHaNcliYO/CEp9z70acjjvQ3sz17g3Yf29tLbznIOvf4wtzpYEx6FxyCv6C8AkqgkKf/rkOwBgjCmUj80Jvj9OBHqcQbxpSupk62PPgcmNX80FDgDDByTnmySALpFxLic92O1+KAAgtOVYMjlniMOS9RFHzHHSRYfO4SvH71sSW903YxnGXP6INStktdDjiINN59CQzzkfurowbSd0Xz8HHxDZHW9kHQkTRwBoHIZhNqr35WtTP2HPAawYhCMYpjHppy69OX3B2sRKOo7ZL54LY9zwiAs5ZORAdtHGOQf+OeiK7faiQANz4tefo2mjUmuZiPdozcvZBNIOFdyctMWuUqG+2zRipYG9kw5zeuY0X4/v3ozlGgGxF889slL7jluVOkB71Nzgc0dr/CAjb9mdSjVZVtqTLH1KAnLj0+8CAD6wmHhXCz2POJSslXhWVsojH/j3D7PXqxPZtglXk3PIE7En9VgdD4UhEF8fPmM0bR6fDmXccvP/9ifG82E2PPowseq6WMXVlo1zuOnzh8f+/vYZSi4MQ8Nxz1YeuqltW7GIPMc5aH2YRC+xOFsWwjxe83TmEIiVosOOjXPwaStqx36q1vVjadrW2wKA3g153H7+JPZaXeegB66UdfQ2OeRykjhYOIdwrJf+cVasXBJhp76FKStllbQo+uXclvsTZ0lVbfQ84mDhHNST2oF79GevVzdL6yas/XbLF47wH6SGHLlzGavry8rRpKRaJrF3Qz4Xe0YNBg4hDXHQ0Ud3ZtI5hxQ6Bj3Htrqpm97jXRce7RzjiQcOj322tBVZ3Yveh+7EyPVj5QA9Vm4+5+Y4fTlDL2slg7LWBVuudiCwsjrpoKRyGJBmsVGHHEGNggnax5HPkTNTnxzrq4vj+jEpcdAPKEn/mmTDkpiahidD3hSLoiSO40xmq42eSxwMOgcJo6mh8sRsJzv9JzXSYlrkiPD8O2tjZQePjJ8cyfPEqf7ksy24CI08MZn8P9L4OegLK0EctOt04pCGWPiIziaOHlTa9E13QUQ4cI/+pTYC4pAUf+hdrN/GW6jFOQdLEhiPt9faXsSyDVFqT+4KX7FSjHMwjilAWs6Bg8rR2ESTqkc5wL9LeY9OzoHImanP9Ly+GSaH0mmxbhXHXd7q4hxCC7bP/PrlUllzHXQOPS58hg0x4uBgnQFg4dptbB29rUqhD2XGlScnNs6YQrqKnINt81C74U7LgB8BMukcdBmza3GnSaXo8tDV+7Q9N5nZTAiBXQbOQd+0vsAE1NPHYlJIB32axyyxcpMml+ZO6NpQj9xbD/4QwEdsWSLyVYgurfZgmluynhotgFt3JTNSR5/SWslH56Di4hP2LQUb1ENr6wYuHFGXdVy6nJlLokRFLs/3aqAizoGIFhPRW0Q0i4hmhmWDiehJIno3/NxNqX8FES0kogVEdKpSfmTYzkIiupFq6EceGbImu1DZeKO9tzI01Qoh2Vb0aH+q5e1NC32zH9a/CX21RPF5D8IGaJyDz6neYRYrYUp1mCblqE4c+jTG71FvauzQuJ9GGmmGuinaTX/jnxyIAsIg4xg1NdjFSqcfOiJhliuhEjirh7THc/Uhlvpc+c15vGw/jUI6NefgGKedc4ibpto4h1aHvCufi9LxugigimlH7aV4WMd/k0EU1fHqaCsRB/8+u4q10glCiIlCCDmrLgfwlBBiHICnwr9BROMBTAMwAcBUADcTkTwa3gLgIgDjwn/xqFhVhM1VXb6EHJlflC/Zkovp8L0GlZS35cKH9ff1c6DYhu5u15fQmBawn84hqKRmQFv0w48n7kO1vPnVPx2Ja86cEPvdx/RUQuUcbAQsGoODcwDwp9eWG9vLeZy8AeCV9yNZdhqxJQfdrJIVK2kNVeIEJy9dvcVupnrDZ+OHJddbc5nuujgHeU8+Ph9F4cjUxxSPHtynZCWlWyvpOh9WrFQ0cyoA/671CAm1QC10DmcCuDP8fieAs5Tye4QQLUKIRQAWAphMRCMADBBCvCyC1X2Xck3VYZuIctG6nMh8IE/SvjJdG9ZbPKgl1G58TsJAUknLwZfQ9DM4u3lxJ2GVmNklc90Re0eJV44/YFhi/GkUoep9cZ7Z0dikWMncFlFwSr/qL3MA8GGb434o5raOGhOJdezEwYdzSI5Th29YDx+FtK+Pb1pnL44TU/uMcQ4Mt5XzJA4U6hyKls3aLVKLP/SEWIlTSLcXreuE+82mj6oWKu1BAHiCiF4joovCst2FEKsAIPyUrqcjASxTrl0elo0Mv+vlCRDRRUQ0k4hmrl27lqviN2KYOAf56XOStKOcdItfPGass45J/urjvwAAC9dEehKfIIJ2sVLwuVufBgzoxSeA8UGUXa8Ya1eHukjSJAHioHJN6jNJjM1DXq2Heue4KF8rtxMOjDy17X40brhs7jn4+L6YN8j438aMcSmtzoZZcoDoQ+HGL8vufmWJtZ98jtxiJcOS4QL3AUl/Fu6VtrXbOQduKJVKI3xQKXE4RghxBIDTAFxMRMdZ6nK3LyzlyUIhbhVCTBJCTBo2jPe8dMHm5yDLKmXngWjzsRH4r564X+z7t89IJvtQMaRvI+Zccyr7W2zxWvpcp1jJmOTesXYdiljA32PXBNmFPNk9843jrf2ZxjVRS+kocf2nD2H69D3lxsfIgRCX77ucumxt9VKepUsJ7oIuZjOZZ6swi5Xc49fHpBI6nz5MGNrfQhy0v7lAmHp3MtIsN652EYb1NozRxTmYxEo3fu7wUh9qfwDQ0l60Wp/pfX7tlP2NdauJioiDEGJl+LkGwP8BmAxgdSgqQvgpAwctBzBauXwUgJVh+SimvCaw+TnI92ZTXPluKpLtsy3i/1Re8hEGKxEVvRryFsVvsm8Xph012lnHLgt1i+F8oHIOOQL2NgQEdFlk5XOEqRP2SJR/9ijeMijN2GyLV3/FXE1fJ0XOC5jvNPqqK+YlVDPKqz4xno3Aq6MS4uDLJOtdcBu6WmI7rOnj5ZauPlc+NoH3mchRqJC2rH/57g7fK34QMYXPkNxwY3hYVLmgG6cFBGP5xp1WfxT97rkkRLVA2cSBiPoSUX/5HcDHAMwB8BCA88Nq5wN4MPz+EIBpRNRERGMRKJ5fDUVPW4loSmildJ5yTdVRslZi55ufbNIH0vLJXzlsfhWnHxosai7TGteWrzjSxzvWPv7g00ex7TOOXW180LpSPY/dp9p2bj46h8AE0h5lNC5WMrfVyypfV/sMPi8+YV8jp6XO44NGuLkG29h8dA7+WRLdYqUcufsDkuuGG39C6W4RnRWFQGuxaDSw0EWgEvL56H4N3/jTm7Fxqut0SL9gPb+pJNey9ekaf7VRiZ/D7gD+L7zZAoA/CCEeJ6IZAO4jogsBLAVwLgAIIeYS0X0A5gFoA3CxEEJqkr4C4A4AvQE8Fv6rCSLOobYPWE4WX07DRhzkZLARB1+dQ1pw8fH1Pn05B9MJt6SQbi9WHJKk+sQh+HTJ/+PJaxixUg7W3yV8uT4fjkYVddnmlwrTfPVxgvNdUvpmN3dlMtifj0MqkNTvcM8jkY3Q8IwlkTfl5AjqBJ9zVsTHrObFUCHvjXv+g/r46ekSepU6uS6XTRyEEO8DSBjwCyHWAzjJcM21AK5lymcC4AWBVYYtE1w6S3k7CiVrJd/6tg3DTRxUglCpmEcFFx9fQq4x343nDCbZDBAtXpMDmV6vnpDP0mU5pEY/3WdYkgiq7+fROatw/TnJ1Ktqfy7I5mzV+zUV0Ny6y7tdWxUf01+1XPfgV+E1lhinYq6nEweurj6nuNhXwbVUyqVt8jExEk+pc1A4B9XclCUOvf3EQ/rz3uWILlst9NzwGZbfqoEGD7FSrL5NnBK2sVsfm1hJ/V69TdTmiJWWc1AV8CpKCun2old/NlSbI6TSPdoc0oBWRYTDxQFSx87lwpDwfZbSCsbG2akObT4e+9ZwHcrlpjhaKu7/Mh+4EvATZ7r0SxK6xR1XVycOpmcxb9UW/H1+oCJdun47W8do4spYK72zOoo+q4peB/VpQP+mAhuxloM+7bfsrH3QPaAHhs8ovTrmJVeTHstF7i1WKpjryVNJ3yaLiKeK3IIvJAttihOkw7Txl2LkiypYilX5Mcghu8blCnDnqwfyDbsi+7NxWofvpfpMuAdg417VqMAmHdO2ZrWOjcjH/+b0LDFO2LKGVM7h8L0GsalQdWs6HwKshqpQYYy5xIiV1ICAQxVF9IwrT4YQ/ro6vc96OMABPZA4SPaAO2Gm8bAFgO9YTE/lYvSPl+8+taVZcPXAjEXBAvpgS2Wx5X3DVHulHE20Xe6oAsiToEkUEfTpThLlKxKTm8whhiRSEiXOweMUD9g3fglbpM8ByinXxOW6wliUftOexTTGmszXo1zdYH949iGsY6fuROezJk3jNzl7RmKlqEyGxfjW6QdhjKJv8xXDSuhDqRdx6HlipfCTmx9pOYcvHJ3MHiYh2+fM9DjYThHyNOKbxKca+J9pE52mrmkC3dkQy39dIQHUCcjfLjO73px9BOtrGcMfXw38NpfZvKhzfqEZfCAPCa7qkmj5+pj09vCGt+EgJX+EiUuVxMElLtHfkdO6y7JLqfG3TAcLH9GTDtN63GQltgAAEAVJREFUGmbwuZBjVNeEfEcuQu+ENpaWKlhT+qDnEQeLzsFXPCJh29DvfnkJAOBvc815eVXYZOXyNGLzaK6mvgQAzpw4Etd9mleaVhuxMOhVJoDjdjebcP7sM3z+bw7W8AYenIN6uW1v8tU5lIiDp+mraVOrJppDzuPTR9i9d/V7dOWQtj2TEQOjtKCm+eH7jPorXIHpdfdv4jkHSdRVj2jpL+XDtdmg3349IrICPZI4mD2k31qxGQCfpJyDbcPQ7Z3dbZl/k2O2saPVOsWnQbWYFV+xkp9COoLPe5yw5wBccdqBxt+v+sT4RLuJPsl9miNP7khuJK47lZ64vmIlWxyt/zAYCqTF+D0D7mLKPoOt9Vx5OgBgf8Wb2/be1bAtJiLSp7GAMUPcqVVVz+mD9+RP+2b/h+CzPcY5BHOi0jhI+v1/zhDuvdrocTqHKGS3GTYzvFrBxur6iJVc6Ql1mHwO0iDgdionSr7e3T5rTH0OPvUf+eqx1t9H7RYQGBvxVRPOSIdFG6x6FaVNG9KKlezjqc4Z8Zj9huLVK0/C8P69rPX21Ig2N/dHDuqNMUP6YPH6HXadg7JZ2ziMD+07FIvXL7WOS73+mk9OYOuYnjcRlTysAeCt5Zvx3prtznH5QL16xpUn14ULBHoicbAE3pNwmUO+cdUp/mEOPGE75cr9zsY5SIWh7tZvwmOX2jdFH+RyAKqgG4txDhWKldQTfFU83T3NLmVfZ0006zF6N+Sxs7UdP7Hk95BKYZtlGhDlTPDlHOoFF2EAkpzDR/fn46RFRMNMmBs8xU8++7Oq59DzpUjonINqhRTEZgrGesZNL0RjrNBrzTccf7XR84hD+JkmVo6OesU2kZCy0BaLlcLA3g2Y9e1T0N8zOmo1MtVVSwlOMc6hMuKgmpRWJwxKsLBtUjtCRMBtxK2xkMPO1nYrAd+xKzAH7dtoX5rtHn4OabHHAPvG/s1TD6jKO1ff8cJrT/MycTZB5Rz6NJif2aJ1vN9CrD9lXCbxkbo533PRFEweE4nQckRsXKaKOQfl8nqFzgB6sM7BxhxU4/n/6p+OAFB53CEA+NzkvbDnwF44cm+7LHdQn8bUIcUrQbXUHNW0VtqlZMiqBnGQcMVWkrA5M8pqtmd/5N6D0ZjP4cvH72sdT1sxnc7BB59xWKddfMJ++IpjXD5Q79/2vmUtm8RUJTQDLeEoXnpvPQDee700rpROnQ35XIyg5HPEind99gDboUg9cNQhjUPUb/266lywipWqQB0OGRWcDntVQSZ88MiBeOkKNiJJ2ajGPUo5/If3HVJRO74KaR9adNioQaVsaq1VCDMgbcptsv0Y52P1hwhgO/0N69+Ed649zTmuYkqFdGeC79zzMQeXIpshDm6+kCO0FQWuPSsZvl1CEoW+niJjfa7midgcGj7E5srTDzL+1jdmrlu/9931ZlYdUA3GTbKX/Q1OMxK//PwRuOuLk6vQY/0hl4FNhu6D+OZqnpJ64hQO3zj1gNJ3l9eyDyRxsEVL9U29Ws2wJlt2BuKnQZaQKr7w9cWpN7zESp7y/AuOGQMgnmkv0Z8kDgZ9gw590y/keZNmn4yLZ1vMf1XHu2rrOm3occTBQ6pUlUU8fEATxg7ti2vPNp9UgMC65TiDQq7TI3yWtgTwPiBPzsHHIqtS5Z8OqXCUZpoc1BFbPd3Dz2psxYeNDkwtXSfmroypBwe5OYZaMsH5+hBccdpBWPD9qdbDh+TofNvU6/VuyKOZ8UHw4e5sjoP9m8rPslgJepxYyZYJTmKQZ0AsG5oKeWOc/e6Coof/RVrYFJ7lpL2sFMfsNzSheNTha231g7MPwQ8fne8dcM2G284/Ch9sbq5qTK36qTr98NUTx+FfPjzGyh1Js06VY+SQyxGacvZTt6TrvqIb/SDTuzGPnbvaE4cYH87BBsk5VMOIJA16HnHw4Byu+LjZKSpDhFrIvbc0myNOSqciF/Yc2Ov/t3fvMVJVdwDHvz+WRWAXFoQF5A0VSJEWkI1AWYsxroBtxMRaJVao2FATTKk1abWPNH2l2IdRxITSFoK1FtO0jdhajSVFY0urSwVhpbxNAanQighaUdJf/7hnupedmTuPPXsf8Pskk7175s7Z3/wY5sw5c+45vHaic+s9hc0YG/2dylkrlkY0lHMuGcKcAjvVVaOhV62XRibNunWTksNmfXvWcuC713jp7ed6DuW+B3dcqbdXj2Cqcuurb5xVHvWm/u3rJjGoxHUL9W6YK+6PRuffsJL7GfVa6l1iGmGWPXLb9Mgrgisx8sLgqtOoZT0q9cKBN4reN3ZgsOd1x3nyHW286wpv8ZTjqbZ//v84rkXRfMpdZzB7QjaHN319l9O+d0d5r+eOF0b2qg16Dvs7TJuNiu9TM0ZxdYkPDLnvLaP2VukK5+67YBHzJg1h/OB6L1eWZlHzuIE0jxvopa71S2ay9eCbsS0X3r+uB68u/1jJ83Jf2kVNW/Qp/B5RzkVgcbrvk5PPWm67kKkj+5eVV5+unTw08nucJORex+W+njt+SHn3/f+y/fBxNu//t9e4etbWcOdV42mZWHjv665y3jUOowbUFd3A/trJQ9mw7bWYI8quIQ09mdvgZ5gk5wst473Us/Gu2Qysi2eZgbAhDelqHKJmwSRpxYKpSYeQJzd8c/pMdO/v8aWzOHLiP1zUcPaqBrm12cKeuKPZS2zLrhrnpZ5KpKZxEJG5wANADfATVV0edwz33ziF798Qz0qk55N1iy8rOa6aM7OT10zkfKCx3ks9lbhhWjrfiE15Gt2sqMaI2VEAk0f0Y/KI6GVq6nrU0PbNud5iS0IqGgcRqQEeAlqAQ8CLIrJBVV+JM45yZjSYyhVbOyfsxqYRPNZ6MHIr1LT74lybyJBl/et68OOFTWWvTxZl1sV+hm6TlJYvpC8D9qrqflV9D1gPzE84JhOjb8y/hEc/M52LB8X/ib+zcr2izl7vYZLXMnFw5HUVUcJblIY3R8qqVPQcgGHAwdDvh4DpCcViEtCztoaPZPTT1volM3hy+xH6lrga3pzbViyYyooFU9n9+kkvS+InLS0fdQpND8ib1isiS0SkVURajx07FkNYxpQ2trGeO64c53V5DJNd4wf38X6lfhLS8gwOAeElIYcDedOGVHW1qjapalNjYzbnZBtjTBakpXF4ERgnImNEpAdwE7Ah4ZiMMea8lYpBUlU9IyJ3AE8TTGVdo6ptCYdljDHnrVQ0DgCq+iTwZNJxGGOMSc+wkjHGmBSxxsEYY0weaxyMMcbkscbBGGNMHtEEdtfyQUROAruK3D0S+EeJKhqA/GUUKz/Hd13lxO7zb/p8jnHH7rsuX68b33FZ7pOL61zM/QRV7VPyL6hqJm9Aa8R9x8p4/Gof53RBXSVj9/k3PT/HWGNPIvcJxWW5t9x7qyvqvTN8O1eHld4s45wnPJ3ju65yYvf5N30+x7hj912Xr9dNuedZ7ttZ7ss/x3ddBWV5WKlVVZsqvS/tLPbkZDn+LMcO2Y4/a7GXG2+Wew6rq7wv7Sz25GQ5/izHDtmOP2uxlxVvZnsOxhhjuk6Wew7GGGO6SCYaBxFZIyJHRWRHqGyyiGwWke0i8oSI9HXlPURkrSvfJiJXhB4zzZXvFZEVEsMC/B5j3yQiu0Rkq7sNiiH2ESLyRxHZKSJtIrLMlV8oIs+IyB73s3/oMfe4/O4SkTmh8iRy7zP+WPNfaewiMsCdf0pEVnaoK/W5LxF/2nPfIiJbXI63iMiVobpiz7035UxpSvoGfBS4FNgRKnsRmO2OFwPfcsdLgbXueBCwBejmfn8BmEmwudDvgXkZin0T0BRz3i8CLnXHfYDdwETge8Ddrvxu4F53PBHYBlwAjAH2ATUJ5t5n/LHmv4rY64Bm4HZgZYe6spD7qPjTnvupwFB3PAk4nGTufd0y0XNQ1eeANzoUTwCec8fPANe744nARve4owTTzJpE5CKgr6pu1uBf7WHguizE3tUxFqOqR1T1b+74JLCTYEvX+cA6d9o62vM4H1ivqqdV9QCwF7gswdx7ib+r4yyk0thV9W1VfR54N1xPVnJfLP4kVBH7S6qa25ysDegpIhcklXtfMtE4FLEDuNYd30D7TnLbgPki0l1ExgDT3H3DCHacyznkypJQaew5a123+mtxd09FZDTBJ6S/AoNV9QgE/5EIejlQeC/wYaQg952MPyeR/JcZezFZyX0pWcn99cBLqnqaFOS+M7LcOCwGlorIFoKu33uufA3BP0IrcD/wZ+AMZe5THZNKYwe4WVU/BFzubrfEFayI1AO/Aj6vqm9FnVqgTCPKY+Ehfkgo/xXEXrSKAmVpzH2UTOReRC4B7gU+mysqcFpmpodmtnFQ1b+r6tWqOg34BcH4MKp6RlXvVNUpqjof6AfsIXjTHR6qouA+1XGoInZU9bD7eRJ4lJiGO0SkluA/yM9V9deu+HXXZc4NWxx15cX2Ak8s957iTyT/FcZeTFZyX1QWci8iw4HfAAtVdZ8rTs17TjUy2zjkZiyISDfgq8Aq93tvEalzxy3AGVV9xXUDT4rIDNctXQg8noXY3TDTQFdeC3ycYGiqq+MU4KfATlW9L3TXBmCRO15Eex43ADe58dYxwDjghaRy7yv+JPJfRewFZSj3xepJfe5FpB/wO+AeVf1T7uQ0vedUJalvwiu5EXy6PgK8T9Aa3wYsI5hFsBtYTvsFfaMJVmvdCfwBGBWqp4nghbUPWJl7TNpjJ5jJsQV4meALrwdws2i6OPZmgm7wy8BWd7sGGEDwxfke9/PC0GO+4vK7i9DMjIRy7yX+JPJfZeyvEkx+OOVeaxMzlvu8+LOQe4IPeG+Hzt0KDEoq975udoW0McaYPJkdVjLGGNN1rHEwxhiTxxoHY4wxeaxxMMYYk8caB2OMMXmscTCmC4jI7SKysILzR0to5V5jktY96QCMOdeISHdVXZV0HMZ0hjUOxhTgFlx7imDBtakEFywuBD4I3AfUA/8CPq2qR0RkE8FaWLOADSLSBzilqj8QkSkEV8H3JrgYarGqHheRaQTrab0DPB/fszOmNBtWMqa4CcBqVf0w8BbBfhsPAp/QYF2sNcB3Quf3U9XZqvrDDvU8DHzJ1bMd+LorXwt8TlVnduWTMKYa1nMwpriD2r5WziPAlwk2c3nGrRpdQ7A0Ss5jHSsQkQaCRuNZV7QO+GWB8p8B8/w/BWOqY42DMcV1XFvmJNAW8Un/7QrqlgL1G5MaNqxkTHEjRSTXECwA/gI05spEpNat4V+Uqp4AjovI5a7oFuBZVX0TOCEiza78Zv/hG1M96zkYU9xOYJGI/IhgJc4HgaeBFW5YqDvBpkxtJepZBKwSkd7AfuBWV34rsEZE3nH1GpMatiqrMQW42Uq/VdVJCYdiTCJsWMkYY0we6zkYY4zJYz0HY4wxeaxxMMYYk8caB2OMMXmscTDGGJPHGgdjjDF5rHEwxhiT53+7qZdYkh/tNQAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + " sorted_data['inc'].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXeYXGd99/25p7ftTdKupFW1JXdLbhgCxrgQQiAJJiZPAuEigRCSJ5UE3uR5Q96EBNJIeAgQJxBq6CGYYEPAFOMqJBfZsmRpVXellbbvTq/3+8c595kzu7Nldqfu3J/r2mtnz8w5e+bMmfM7318VUko0Go1Go7HjqPUOaDQajab+0MZBo9FoNAvQxkGj0Wg0C9DGQaPRaDQL0MZBo9FoNAvQxkGj0Wg0C9DGQaPRaDQL0MZBo9FoNAvQxkGj0Wg0C3DVegdWS3d3txwcHKz1bmg0Gk1DcejQoQkpZc9yr2tY4zA4OMjBgwdrvRsajUbTUAghzq7kddqtpNFoNJoFaOOg0Wg0mgVo46DRaDSaBWjjoNFoNJoFaOOg0Wg0mgVo46DRaDSaBWjjoNFoNJoFaOOg0Wiahm8/P8pEJFnr3WgItHHQaDRNQSyV4Tc+9xRfPjhc611pCLRx0Gg0TUEslQVgLp6p8Z40Bto4aDSapiCRNoxDJJmu8Z40Bto4aDSapiCRzgEQSWjlsBK0cdBoNE1BXjlo47AStHHQaDRNgTIOc1o5rAhtHDQaTVOg3UqloY2DRqNpCrRbqTS0cdBoNE1BXBuHktDGQaPRNAWWcijiVvqVTzzJXz94tNq7VNc07JhQjUajKYVExog5pLI5kpksXpfTeu7gmWl+fGKCDa0+3nrrtlrtYl2hlYNGo2kKkqZyAAjb1EM2J4mns3icDv7iv19gJpaqxe7VHdo4aDSapiCeyhsHu2spljIe7+wNkZMwE9MV1KCNg0ajaRISGZtxsAWlo0ljeVfIs+B1zYw2DhqNpilQdQ5Q6FaKmsqhM2gYh6Ttdc2MNg4ajaYpiKcXUw6FxiGR1soBtHHQaDRNQqLAOKRtjw3j0B3yApDMaOUA2jgs4OjoHGcmorXeDY1GU2aS6RztATcwLyBtxhy0cihEG4d5/PHXDvOBB4/Vejc0Gk2ZSaSzljqYWyLmkNDKAdDGYQHhRIbZuE5l02jWG/F0lja/G7dTFM9WsgLSWaSUDE/FarKf9YI2DvNIpLPEtKzUaNYdiXQWv9tJyOsqcCstCEhncjwyNMHL//YHnJ+J12Rf6wFtHOaRSGdJpLRx0GjWG4l0Dp/bQcjnKlQOC1JZs1ycTZCTMBVp3mppbRzmkczkiKV110aNZr2RyGTxup20eN2FdQ7JDH63E7/H6LWUzOSsoHSyiQvidOM9G1JKEulsQZm9RqNZHyRSplvJ5ypIZY2msgS9TjxOB0IY3gOXQwDNndaqlYONdFaSk2jjoNGsQxIZw63U4nUtUA5BrwshBD6X0/AepLRy0MbBhuqpEjOzFTQazfohkc7iczkXxhySWQIew4nidTuMuKNyKzVxKw3tVrKhTgQpDTnpczuXWUOj0TQCUhptuX2LZCuFvMZ33edykkznUPeG2q20BEKIzUKIHwghjgohjgghfsdc3imE+K4Q4oT5u8O2znuFEENCiBeFEHfZlu8TQjxnPvdhIYQwl3uFEF8ylz8phBgs/1tdHntlpHYtaTSNzdefHrGmu6WyxgXf7zGUQziRsbwD0ZThVgJTOWSy2q3EytxKGeAPpJR7gJuBdwkh9gLvAR6SUu4CHjL/xnzuXuAK4G7go0IIdQv+MeDtwC7z525z+duAaSnlTuBDwAfL8N5Kxn4i6FoHjaaxue/h03z5J8NAviOr1+WgJ+Qllc1Zxa7RZIag6VbyuZxGUoqZsaiVwxJIKUellE+Zj8PAUaAfeB3wafNlnwZebz5+HfBFKWVSSnkaGAJuFEJsBFqllI9Lw2R/Zt46altfBW5XqqKa2Fv6xlM6nVWjaVTGwgmOjs4xG0+Ty0lrCpzP7WSgIwDA8JRR4BZNGtlKYCiHZCZneQ6auc9SSQFp091zHfAk0CelHAXDgAC95sv6gWHbaiPmsn7z8fzlBetIKTPALNBVyr6VgwLloN1KGk3D8uPjEwDkpNESR7Xr9rudbO70AzA8bbTHiKYyVkBaKQfLrdTEAekVGwchRAj4GvC7Usq5pV5aZJlcYvlS68zfh7cLIQ4KIQ6Oj48vt8slU6gctHHQaBqVH5/IXx+mYynru+1zO9ncaSiHc1MxpJRmQDofc0hmcpYx0W6lZRBCuDEMw+ellP9pLr5kuoowf4+Zy0eAzbbVB4AL5vKBIssL1hFCuIA2YGr+fkgp75NS7pdS7u/p6VnJrpeEjjloNI1PLif58YkJq5GeYRyUW8lBq89Nm9/N8FSMZCZHTkJAuZVcThLpvFtJB6SXwPT9fwI4KqX8B9tT9wNvMR+/BfiGbfm9ZgbSNozA8wHT9RQWQtxsbvPN89ZR23oD8H1Zg0IDrRw0msbn1ESUyWiKO6/YAMBMLG0zDoYR2NzpZ3g6btU7KOXgcztI2t1KWjksya3ArwCvFEI8Y/78NPAB4A4hxAngDvNvpJRHgC8DLwDfBt4lpVRX2ncC/4YRpD4JPGgu/wTQJYQYAn4fM/Op2uhUVo2m8bk0lwDg6oE2AGbiKctNpIzDls4AI1Mxa9CPVQRnVkjHdRHc8kVwUspHKB4TALh9kXXeD7y/yPKDwJVFlieAe5bbl0pjVw7araTRNCZjYcM47O5rAWA6msbvzisDgM0dAb73whhzCSOd1SqCMyuktVtJt88owH4irDSV9WuHRnjX55+q1C5pNJoSGZtLArCzN4QQMBNLWd9tpRwGOgOksjnOTBojga0iOJeTeDqrA9Jo41BAYcxhZSfFoycn+NHx8mdOaTSa1TEWTuJ3O2n1uWjzu5mOpS0l4Fcxhw4jnfXFi2Eg71byuR0FaezaOGiAfMzB7RQrnukwE0s3tfTUaOqNsXCS3lYvQgja/W5m4sUC0kY669FRIys/H5Au7KfWzN9tbRxsJDM5PE4HQa9rxQHpqWiKdFaSzeWTqx4dmmA8nKzUbmo0miUYm0vQ2+IFoD3gYSaWIpFRdQ7GJa+/3Y9DwMGz0wAEPCqVtfCS2MwBaW0cbCTSWbxuB363c8UV0jMxY4xgyjz5Ls0l+OVPPMmnHztTqd3UaDRLMB5O0tviA6Aj4C6sc3CpwLOTd991OTMxFZAurhwSTawcdMtuG8mM0dLX73FaAanlmDZPrmQmi9/j5MHnRpEynzGh0Wiqy3g4yU/tNpRDR8DD8UsR4uksHpcDhyOfePnOV+xgd1+IH7w4RnvADRQqB6PmoXmVgzYONhLpHF6Xg4DHuSK3UsbW2VEFrh547iIAU9H0outpNJrKEE9lCScz9JhupbaAm9l4mplo2lIHdm7f08fte/qsv+3Kod3v0QFpjYGlHNxOYitIZVWGAQy30qW5BD85a3T9mIrqmINGU22UYlcxh46Ah0gywxOnJ62iuKWwK4f2gFsHpDUGibQxY9bvcRFfgZxULiUwDMt3jlxESti7sbXgOY1GUx3GzESQ3tZ8zAHg7GSMm7Yt3+jZrhw6Alo5aEwS6Sxel5OA27miIrhpMxhtrJtjZDqO1+XghsEOJiNaOWg01UYVwNmzlRQ3butcdn2ve55yaOKYgzYONoy50Q78npVlK01H88YhmcmRSBtB6c6gl7lEhnS2eU8sjaYWFHMrgRFcvqp/JW4lW8wh4CGZyVKDHqB1gTYONhLpLD6Xka20kglQM/PcSmr9zlC+VbBGo6ke4+EkLoewjILKQtq3tQOPa/nLnW+ecshJyOS0cWh6EmkjIB1YYZ2D/eJvdHLMGcrBPDGnoto4aDTVZCycpKfFa6WsdpozHW4cXNlgyQLl4DcMS7PGHXQqq41kxkhlVXUOUkqWGmU9ZTcO6ZwZs3BYJ6Q2DhpNdTk1HmHA7JsEsKndz9++4Wru3LthResr5aDcywDJdLZoGux6p/ne8RIk0jm8ZhGclMbf6gSxc+jsFOemYsxEi7iV3E5tHDSaGpBIZ3n+/BxvvXWwYPk9+zcXX6EISjn43U4rrVUrBw3JdBaf20HATGeLmwHm+Xz8R6f48Ylxbt7eRcjrIpLMkFIBaZtxmNbGQaOpGkcuzJLK5rh+a8eqt6GUQ8DjsgzFSuKP6xEdc7CRsLXPABYthBueipFI53jy1BS9rUZWhJGtZGQ7qSDYpDYOGk3VOGQ20bt+y1qMg6kcPE7LUGjl0MT8+mcOclV/G+msNGMOxmEp1kJDSsnwVMx4Pp1lY5uPU+NRa7Sg3+PE7XTQ5ndrt5JGU0UOnplma1fAap2xGlwOgUMot5IZc9DGoXk5dHaa0dk4gJWtBBTNWJqOpYnalveZlZj2VFaArqBHGweNpkpIKXnq3DQ/tbtnTdsRQljeAyvmoN1KzYmUknAibU2E8rkc1p3H6OzCzqpKNWzrDgI245DOWQFtgA5tHDSaqnFuKsZEJMW+NcQbFKr5prfJ3UpNbxySmRzprCSdNQpdvG4ng+aF/6w5X9bOOdM43LN/AIDOgAePy5GvkDaNQ6c2DhpN1Tg1YXxXL9/QsuZtqeabze5WanrjEE4UBp19biNe0Bn0WMPH7QxPG8bhDfsGuGlbJzds68TrcthSWY1D2hnQxkGjqRYXZgy38KZ2/zKvXJ6X7OjmhsFOWyprc7qVmj7mEE4Udk9VMYPBrgCnJ/LG4WuHRrg4l2BkOk5n0ENvi48vveMWwJCh0WSGTE5a2Q6dIQ/TsdSyhXQajWbtXJiJ43QIawLcWvj7N14DwLlJ40awWZvvaeOwQDko4xDksZOT1vKvP32eJ05NsrM3ZA0nV3hdTubixnaUW6kj4CadlURTzVldqdFUk9GZBBtafTgd5bsR0zGHJieSNC7q6uZeScnB7iAX5xJWOutcIk0mJzl2MczmjkLp6nU5rME/+fL7fOm9RqOpLOdn4mxqX7tqsNPsbqWmNw7KrbS71whkqWwjKyg9ZbiW7FPf5isHT4FxMNZXJ1aiSe86NJpqcmE2XpZ4g518hXRzfoeb3jjMmW6l67a0A/k7/21dhnE4Y8Yd5uJpq/J5c8c8t5LbucA4qN/NWnqv0VSLXE5ycTZRduPgaXLl0PTOcBVzeNWePv778CgbzLqFrd2GATg9EUNKyVwiw6++ZBCnQ3DH3r6CbXiLKgflVmrOuw6NplpMRJKks5JNbeV1KzkdArdTNG3MoemNQ8Q0Dq+4rIfn3nenlVnU6nPTFfRwdjJKLJUlm5P0tnh5x8t3LNiG1+WwYhcqIK2CWYkmvevQaKrF+TKmsc7H53I27Q1e0xuHcCJNwOPE5VzoYdvaFeDsZMxSBa3m8I/5eG0TpqyAtFYOGk1VUJ0MKmEcvG5H07qVmj7mEE5kaPEVt5F9rT7GI0nmzKB126LGId/WOx9z0MpBo6kGVgFcWwWMg8vZtG4lbRySaVp8xS/6XSEPk5GkVcPQusjrCpXD/JiDNg4aTSU5PxMn6HHS6i+/I8RrtsZpRrRxWEI5dAW9TMfSVhuMxU4+r7uIW6nJC2g0mmoxOmNkKlWiE4HH5Sjaur8ZWNY4CCE+KYQYE0I8b1v2PiHEeSHEM+bPT9uee68QYkgI8aIQ4i7b8n1CiOfM5z4szE9SCOEVQnzJXP6kEGKwvG9xacKJzKIVzN1md1bVY2klbiW/TmXVaKrK0YtzVpfkcrOjN8ThkRmklBXZfj2zEuXwKeDuIss/JKW81vx5AEAIsRe4F7jCXOejQgh15fwY8HZgl/mjtvk2YFpKuRP4EPDBVb6XVRFOpBd1F3Wb4z5PjUeAUt1KZsxBB6Q1mopxcTbB2ckYN27rrMj2b7usl7FwkiMX5iqy/XpmWeMgpXwYmFrh9l4HfFFKmZRSngaGgBuFEBuBVinl49IwwZ8BXm9b59Pm468Ct4sqdqpbyq2klINqwLfY64oZB6t9hg5IazQV48AZ49J007auimz/FZcZw4N+cGysItuvZ9YSc/gtIcRh0+2kJmz0A8O214yYy/rNx/OXF6wjpcwAs0BlPukiLB1zUMohSnCRdFfIV1J6nA6r8ZdWDhpN5TlwepKQ18WejWuf41CM7pCXawba+P6L2jislI8BO4BrgVHg783lxe745RLLl1pnAUKItwshDgohDo6Pj5e2x0XIZI25zyHvYtlKhnKYjKYWjTdAPuZgD0y7nA5cDqFjDhpNBTlweop9WzsWvXErB7dd3sszwzNMRpIV+x/1yKqOqJTykpQyK6XMAf8K3Gg+NQJstr10ALhgLh8osrxgHSGEC2hjETeWlPI+KeV+KeX+np61zYqFfEfWxZRDq8+FxzzpFiuAg7xRUMFohc/dvDnSGk2lmYqmOH4pUrF4g+L6LR1ICUNjkYr+n3pjVcbBjCEofg5QmUz3A/eaGUjbMALPB6SUo0BYCHGzGU94M/AN2zpvMR+/Afi+rFJqgOqrtJhxEELQFTJcS4sFoyHvQvLNMw5el0MrB42mQhy5MAsYF+9KEvAY3+tUtrlu9JatGhFCfAF4BdAthBgB/gx4hRDiWgz3zxngHQBSyiNCiC8DLwAZ4F1SSnV1fCdG5pMfeND8AfgE8FkhxBCGYri3HG9sJajK58WK4MDwOY7OJpYssFFuJZ+70Nb63E4dc9BoKsRMzPj+dps3cJWiWZtoLmscpJRvKrL4E0u8/v3A+4ssPwhcWWR5Arhnuf2oBKrpXusiygHIK4clYw7F3UrN3JdFo6k06uZuqe9mOWjWiXBNXSGt3EqhpYxD0AhKL+lWMk8e7wK3klYOGk2lWM4tXC6adSJccxuH5ArcSi3LKwePs7C2QeHTykGjAeDYxTkOnV1pudTKmIuncTnEAsVebiy3klYO65tPPXqaq9/3HdLZHJMRo2dSR2AJ42AphyViDla20ryYQxP3gtdoFFJKfvs/nua3/+PpRV8TT2XJ5RbmoTw7PMNYOFF0HVWjVOmaWUs5NFlySdMZh4DHxVwiw+hMwurmuFQNg1IOS9c5LJKt5Hbolt2apufQ2WlOjEW4MJtgdDa+4Pl0NsdLPvAQX3/6/ILnfu0zB/nI94eKbncuka54vAF0zKFpGOg0er4PT8c4Px2nv2Ppbo5WzGEFRXAL6hxcTp3Kqml6/uPAOdRX7KmzMwBkc5IPP3SCCzNxoskM07E0w9OxBevOxtJcmCmuHObii/dFKyeq1kkbh3XO5g5jNvTwVIzzM3H6l5kede2Wdn7m6o3s27p4LvVSymG9nFCpTI7hqYVfXo1mKWbjab51eJQ37tuMz+3g0NlpAJ44Nck/fPc4Dx29ZH1H5rfGzmRzpLI5xhepTF6q9U05Ud0Omi1+2HTGYWObD6dDGMphxlAOS9Hqc/ORX7qebrOVRjHy2UoLYw7rRTl85dAwd37o4XXzfqrFgdNTfODBY03Z8hnghQtzJDM5fuaajVw90M6hc4ZxeOC5UcDoPabicvF551bM/Ht8bhHlsERH5XLjdTmaLn7YdMbB5XSwqd3HsdEwM7E0A6aSWAuLupXWkXIYnUkQT2etliOa5fnvwxf45X97ko//6CST5sCoZmMyatz197X62Le1gxcuzBJLZfjOkYuAMe9ExeVi85SDUhLjkWRR41ot5QBGmvp6+S6vlKYzDmC4llSr3+XcSish4HHicTroDBZWanrd60c5KKPQrFOxSiWVyfEHX37WUpPTzWoczIzAzqCH67d0kM5K/vJbR5kwlycy2cWVg3mupbOS2biRdh5OpLnv4ZNkc9KIOVQhIA1qXGhznftNaRwGOvxWAc1ybqWV4HM7+a933co9+zYXLnc5SKRz68KloKpR14uxqzTnpmIkMznu2NsHYI2abTYmI0mEgI6Ah1t2dHH5hhb+48lzeF0OPKarRimH+TcesVRepY6HDQXy/WNj/NUDx3hmeIZoKls15eBpwlnSTWkcNttcSQNlMA4Aeze14vfMD0ivn4ZdqtXI/Ls7TXHU9MD9W42OodOxJjUO0RSdAQ9OhyDkdfGt//0yPvq/ruef7r2OFq+rUDks4laCvHFQCuL4pTCwdOeCcqJjDk3C5k7DOHhcDqvIrRKsp4E/Smlpt9LKUNMD9w8aWW5NG3OIpKz+ZABOh+Cnr9rI3VdusBpTKjUaW8StBFgZS7Nms70XLxrGoWoxB5dTu5Wagc1mrUN/ux+Ho3LVldao0HVwt23FHNbBe6kGp8ajdIc8bDFvRJo25hBNLojFKVRLe+WuSSxwKy1UDsq9qYxDdWMOjX+TVwrNaRxMt1I5gtFLkZ8j3fgnlTIOOuawMk5NRNjeHcLndhL0OJmKpmu9SzXBUA7F1bl3gXIozISL2/4eU8Yhbix78VKVlcM6yjxcKU1pHHpavPjdTktBVIq8W6nxL6hh845NK4eVcWo8yvaeIACdIQ9T0eYaMamYjKboXkQ5qMaU+SK4wouvUg4el2NBzEEF+KsXc2g+t1J1zG6dIYTg47+yj+3dwYr+H6Uc1lfMofHfS6WZjaWZjKbyxiHgYSrWfMohlckxG08vrhxUtpJ5wzH/JkrFt7Z0Bha4lRQ6IF05mtI4ALx899pnUC/HeukDn8rk8nd3Wjksy6kJI1Npe3cIgI6gx8r3byZUhtZiMQef28lUNGWdW7FUBiml1etMKYetnQFGpo2GfQuMwxITGsuJ1+VYF1mHpdCUbqVqsV6Ug70qej24yCrNqXEjU2mbUg5BT1PWOUyYGUaLjfFU7WXUzVNOFqZ9x1JZPE4HG9p8+WyleBp7n8yQt4rZSg3+PS4VbRwqiJop3egX1LDtbk2nsi5NNJnhPw6cw+92WplKnYHmNA7qPS/mVlLtZew3T/bzK57K4Pc46WnxMhVNkc7mmItn2NZlGN2gx4nLWZ1LWDOO/NXGoYKslwlSKt4A2q108MwUH37oBA8dvUR23nAaKSW/8blDPDM8w9/dcw1u88LVEfQQT2craliPjs7xvvuPFB2YUyuUK61r0VTWQuUAhedXNJUlYBoHMDKW5hJpLt/YAlQvjdXYV52tpCkj60U52N1KzW4c/vF7J/iH7x7nbZ8+aDWPU8zG0/z4xATvesUOXnP1Rmu5ujhWskr6ey9c4lOPnamrSmzlVlpKORiprIWuJEU8lcXvcbKxzQfAiUthpITL+lqB6qWxgspW0sZBUyasmEODy1G7cphfqNRsjEzHeMVlRjKDqoJWqOM00FnY6bfDNA6VdC2pQO1MvH6yoiajKdxOseiIXZ+7iHKwnV+xVIaAx8nGNiPlXBW+bWz30epzVS1TCQzlkM1JMk0UlNbGoYLkZ8829gkVSRoXHI/T0dTKIZeTnJ+Js2djK11Bj5VBo1AX6PkXw85qGAezOGy2noxDxKiOXmzSomqDbVfW9vMrlsoScLvYNM84tPrcbOsJ0dtaudY3C/e1+abBNW0qazVYb8qhp8Xb1MZhLJwknZUMdPgZ6PAzMm+spTpOLfPuaDuXcCsNjYXxupxWv6/VogxTvRiHXE5ydjJmjdkthrp5KohppQoNRUfAQ6vfRcDj5KjVMsPFP//SdXhc1bu3tccPK9iOra7QyqGCrBfloL683S3eps5WUsZgoCPAQEeA8/OUQ8QyDvOUQ8AwDsVqHX7vS8/y5998Yc37pozDXB0Yh0Q6y5s/eYAnT0/x8ssWrydSN0+z8bQ1KCuWmqccPE6EEGxs83FyzKgfafO7GegI0Nviq+C7KGS91CyVgjYOFUQIgcflWBfKweN00O53N3xwfS0oN9JAh5/+Dj/nZ+IFszrCpvttvnJo87txiOLKYXQ2wfmZ+ILlpaLcSjN1UIn9+MlJHhma4L2vvpw/uuuyRV+nEjZmYmk6AsYxK3AxmQFpgE3tfqsGopqxBoXlVmrwG71S0MahwgQ8TqvNcKMSSaYJ+Vz43c6mdisp5dDfbriVkpmcVZwFdrdSoXJwOAT9HX4ePj5ekGqay0mmYynGzBnJ33jmPG/+5AHe+PHHuVCiwQjXkVtJqZjb9/QtGm+AvKtmNp6mzVRXhTEHIyANWBlLUN0UVoXHuT7S0ktBG4cKs29LB48MTTT0NLhIIkPI68LvaXbjEKenxYvP7bSGRNmD0osZB4Dfv2M3z47M8oWfnLOWzcbTZHOSSbPA676HT/HkqUkOnJnihQtzJe3bXKJ+AtLqOCyWpaRQymEukVcOC91KxjZUxpIQ0FKlqmg72q2kKTuv2tvHyHTcajHciKhB7j63s6kb741Mxy2jMGC2fbcbh7lEGo/TYd0R23n9tf3csr2LDz54zBp/OWnr1DoeTnJ+Js41A+1AadMDpZRWrKEe3ErKOISWMw7mcZLSGCMKebdSNidJZnJWLEIphxavq6IzWBajGbOVtHGoMLdf3gsYRUqNSjhpKgczL71ZGZmOWTNA1G97xpIyosUQQvBLN21hLpFheMowKBO2APXpiSgzsTTbzE7BpdyhxtNZMqa7qh6UQySZxukQ1oV9MXy250NeF06HIJbK8MFvH+ORoQmAvFvJPN61cCmBLVtJxxw05aK31ce1m9v57tGxWu/KqlEXPb/HqHNoZBfZalE1DkoxBL0uOufVOixlHAB6zTYQqv20ve7hqbPTAFab71QJd6gqGG08rr1xCJtuyKXiDZC/GwfDxRRwOxmbS/KxH57kU4+eBvLGYZOpHNpqZhy0W0lTAW6/vJdnh2cadlRkJJmmxefG73aSzUnS2eYzDvYaB8VAh78gnTWSSC/IVLKjegSNR4wA9KQtmP308AyApRxKMg5mAFiIOlEOyxhJhc/mfvO6nfg8To6ZtQzPjswC4FcxB6UcapCpZOyfditpKsCuPqNRWDlSFmuBCkgrN0AzBqUvzBqf3ab2fNZMf7t/xW4lyBuHibBxk6DcSg4BT58rVA6lXISUWuhr8TETr/0NyJx5viyHz64cXA4CHifHzdicUlVKOYS8Llp8rhoqB5Wt1DznvjYOVaCvVXWVTNR4T0onl5M2t5KaT9E8XxBFPgMnf3HqChW24l7OOBgG1mGlv05FU3QE3PS0eJmOGcHs/nbDbVWRxm8FAAAgAElEQVRKQFoph82d/vpQDsn0ipSDd55y8LsXNrdT5xzAvTds5vY9veXb0RJYLwWtpbCscRBCfFIIMSaEeN62rFMI8V0hxAnzd4ftufcKIYaEEC8KIe6yLd8nhHjOfO7DwnRICiG8QogvmcufFEIMlvct1p6+VuNu89Jc480RnogkyeQkG9t8BVWsKje/WYianWmDtjviNr+buUTGisGEl3ErCSHoafFaMYfJaJKukNc6Pza1+6yLUCluJWW4NncGSKRzNb+7jSQzSx4HhV05eF2OAkOgCNiC1n/ymr3cs39zeXayRPIxB20c7HwKuHvesvcAD0kpdwEPmX8jhNgL3AtcYa7zUSGE+nQ/Brwd2GX+qG2+DZiWUu4EPgR8cLVvpl5R7oRLDXhBHTZ96v0dfss4fPngMC/9mx80bAxlNai25aF5xiGbk9Zz4RW4U3pCeeMwEUnRGfRYger+Dj8Oh8DtFKtyK202g+W1Vg8rOQ5gqAX7Y3V+Xb6hxVqu6hxqjdrXWhvearKscZBSPgxMzVv8OuDT5uNPA6+3Lf+ilDIppTwNDAE3CiE2Aq1SyselcZv1mXnrqG19FbhdLJfm0GC4nQ66Q56GVA72fkI+887uyVOTpDK5ho2hrIboIsYBjItxLieJpDLLFn7ZlcNUNEV3yEOvqRxUeqzH6SgxIJ1XDkDNK/JXHJCerxzMC/DeTa3WsSimJmqBdiutnD4p5SiA+Vs5AvuBYdvrRsxl/ebj+csL1pFSZoBZoKvYPxVCvF0IcVAIcXB8fHyVu14belt8deWKSWdzK7oAKQPQ355XDkdHjaDhZDMph0RxtxIYxiGSyiDlwr5K8+kOea2Yw2QkSVfQS5/ZQE6lyXpcJRqHeBqvy2Ep1LpQDiswDh6nw5oH7XM7LUOwtTPI7r4QkA9I1xqXQ+AQpcWCGp1yB6SL3fHLJZYvtc7ChVLeJ6XcL6Xc39OzeLfHeqS31culOgpI/7/feJ5f/8zBZV83Mh2nI+AmaBbBQT5byZ6Kud6JpIzmg/Y20a0247BU6ww7ah5yIp1lOpY23EpmwoKlHEo1Dok0rX437bb9UUgprQygapDMZEllcytKORVCWHfkduUw2B1gt+laqhfjYOxrc02DW61xuGS6ijB/qwqvEcAeMRoALpjLB4osL1hHCOEC2ljoxmp4+lp8deVWOjoa5sxkdNnXnZ/OF37Nl/gTTWQcoskMQW/h+1fKYS6etrXrXvqiqO7uh8z2090hj1U7MdhtUw6lZCvFMwVpnnbj8M3Do9z5oYc5NxlbbPWyYrXOWGH/I5Ue7XM7LUOwtSvIL+7fzO/cvqtmqavF8LodJJsoU2+1xuF+4C3m47cA37Atv9fMQNqGEXg+YLqewkKIm814wpvnraO29Qbg+3IdluD2tXqNzJ86kaXj4aR1QVuKkemYdfGa3w6h2HyC9Uo0mV3gKskbh4zVFXVZ5WDOUz46ajTW6wp5uXVHN597201cv8VI+jPuUFd+EZpLpGn1ua39sfdX+tGLhvu1Wqp1sZkWi2FXDj7LrRRge0+I37tj97JV1tXE63I0lXJY9hMUQnwBeAXQLYQYAf4M+ADwZSHE24BzwD0AUsojQogvAy8AGeBdUkp1lr8TI/PJDzxo/gB8AvisEGIIQzHcW5Z3Vmf0tvqQ0shQ2dBWvSElxZBSGkHRZb53UhotI267zAgp+dzzlUPzGIdwIkPQU9w42N1Ky/nalXJQlcCdQQ8Oh+Clu7qt16wmIN3mdxe4ucD4/J44NWnuf3njEH/x3y9w47ZO7rpiQ8HyYlldS2FXDq+8rJdIIkN7oH7Ugp1mcyst+wlKKd+0yFO3L/L69wPvL7L8IHBlkeUJTOOynsnXOiRqbhzm4hnLbZHMZIt2ET06Okeb300incsrB5tbaXt3sOncSvMveKpZ3Gw8vej86Pko4/DQUaMRo4oz2PGUeIcajqfZ3OHH6RC0+FyWcRieilsJBfb+S+Xgyz8ZZnQ2vsA4zFkKamUXeNVCw+tycO3mdm7aXjQXpS4wlEPzuJXqI4m4CVBV0vVQ62Cv1I4mFxqHyUiS1/7fR7h8oxEUVDEHn+kC6GnxsrUrUDDoZr0TTWWsttIKIQSt5sV4sfnR8+k23UpnJmPcfcWGorOjVxuQBugKeqwsssdOThS8ppwkMznOFoljlOpWUums81VpPeJxOXQqq6b8WMohXPsL6rhtH4rFHX5yZopMTvL8ecMv3m8qB5fTgcfpYHt3kK6Qt6liDpEiygGMjKVSspV8bietPhcuh+CP7i4+QtO7yoA0YH4uxuf7+KlJa4hOObu1ZnOSVNYwDvPDgys9DgqvTTnUO80Wc6j/T2Sd0BX04BDURa1DwWjL5MKLxoHT03hdjoLKXUVn0MOeja10m8ZhHeYOFKWYWwmMuINhHFY2wwDg5Zf18pu37WR7T6jo86XEHNLZHKlszoqHdIc8ltF+6tw0L9nRjdflsC7a5UDtWySZsXpLPX9+lnd+7hAXzfN7pTEH1e20EYxDwOMimiqve66e0W6lKuFyOugOebk4W3vjMDa3vHK4dnM7v3nbTv7nyMWCnPXP//pNdAe9fOXQMKlszgqGrnciiUxBAZxCGYfpWIp2v3tF2TX/903XLfm8173yO1Q1VlOlgXaFvBw8M42UkouzCX76Kj8tPndZ3Ur2xotnJmMIIXjHZw9xfiYf41hJERwUBqTrnY6gh+Hp6qQE1wPaOFSRzZ0Bzk3V/uSyKweVXWL/+8iFWd51205evruHl+8uLDbcYd7tdoUM//tkJLnujUMuJ4mmsoS8Cy9grX4356fjnByLWrMY1kopykFdqFWyQHfQw1QsxXjEmD+xodVHq99V1oB0whaUPTcV5WM/PMl4JInH6eDwyCweV/FRqcVQRqERlENX0MNUE7lS6/8TWUds7QqsqPCs0ozNJay2BfONw1Nnp8lJuGGwc8ltdAXN2QRN8GWJmRfgpZTD0HiEXX3F3USlUkpAer5y6G7xIiW8aKbK9rX6aC2zcrAHZQ+emeahY5d4x09t5/qtxvzrlhW6lMAwCi6HwOWs/0tRV9BDOJlpmoyl+v9E1hHbuoJcmktaA+ZrxXgkaaWnzvdFHzg9hUPA9Vs7iq1qobJumqGFRrF23Yo2v5vJaIqpaMpSVWullAppdS6pWIcy2kcuGMkEfa0+Ws3W4uXCrhz+86nzSAl3X7mBG7cZaagrDUaDka3UCKoBDJcdFI53Xc9ot1IVGTTdDmcnY+zZ2Fqz/RgPJ9nWHWJ4Kl6gHKSUPPD8KPsHO5cNKHabbqWJJviiLJWBY3ep7ewtk3FwOlesHOIp5VbKB6TBCBADbGjz0epzMVJGd6ZSDkIYfbb62/3s3dhqVWavNN4A8HPX9TPYVR53XKXJu1JTbGxbWJ+y3mgMk71OUF+CMxO1dS2NhZNs6TSKpuwB6SMX5jg1HuX11/YvsbZBR9A0DnWQmltpLOVQZLaA3TiocbBrxVNCsZVqgmgPSAO8YLbn6Al5KxaQ3mrWaLxqTy9CCK7b0o7LIWjxrjwGtW9rJ7/2su1l27dK0mWe883SjVgbhyqiGqudqVITtGIkM1lmYml6W3yEvK4C5XD/sxdwOwWvvnLDElswcDsddATcTVElvZxbCYyL86YyVb57XQ7SWUkut3yasIo5KLeSUg6nJ6J0hzx4XA4jIF1Gt5LKpNptGsNX7e0DjFTPV1zWw2UbymMk642uJnKlgnYrVZUWn5vukKemykHlwPe2eAl5XZbLJJeT3P/MBV6+u8dSBcuxuTPAqfHaB9grzVL9gpRx2NkbKluTONUWPJXN4XMsnfUTnxeQbvO7cTkEmZyk15wT0epzk8rkSKSzZUkZVcrhF/YNsKM3xM22lhf/+ub9ddUsr5x0BvNupWZAK4cqs7UrWNOMpTHTDdTT4qXF5yJiFsGdn4lzcS7B7Xv6Vryt67d08MzwTN10mq0UqvBpfstusBmHMgWjIZ/WuZKgtKUcTOMghLB846qHl2qtUS7XklIOO3qC/PHdl+O2ZRqtV8MARt8st1Not5KmMgzW2DgoSdwd8ha4lZTRKKUp4PVbO4ins1aH0fVKZImOq8o47ChTMBpsymEFQWkr5uDO75vKJFP9vFQzwHLVOijlsNJahvWCEIKuoNGe5J9/MMS7v/JsrXepomjjUGUGuwI1TWdVkrgr5CHodVkXPhU7UPMGVsI+M9310NnpMu9lfRFJGhfDYm6lgQ4/777rMt6wb2DBc6vF4yzBOKhUVlvH3C7LOOTdSlC+tt1KOXjdzXf56Ax6mIqm+PrT5/nKoRFO1zi5pJI036dbY7Z0GUHpC2abgUrwhQPneGxoouhzShJ3Bb2EfC7CyULjoFwSK2FTm48Nrb51bxyiyQwOsXDYERh3k++6bad1IS4HSjkUa6ERT2ULOvvGUllcDlEwvrTb9I1bxsFvKocyBaWVcmiElhflpitktNA4OW5M8vvST4ZrvEeVQxuHKqNy5aPJylVZfvihE3zm8bNFn5uMJAl4jGHuLTblYCmK4MqVgxCCfVs71r1xiCSNQT/V8qcrd00x5fB3//MiP/uRR6y/Y6nsAqPVbTZM3DBPOZSrM6ulHBqkeK2cdAU9HL8UQUroCLj56qER0us05tZ8n26N8Zu+4Up2d4yns4ummE5GU5Y6sMccJsweSZ4Sv/DXb+3g/Ey8LrrNVopIsnjTvUqxVMzhiVOTXJpLMmsWnCXS2QWzvbsWKAflVirPOZdMZxEi7/5qJrpsbtd333U5E5Ekjy6i0hud5vt0a4zKeIlVUDkkljAOE5GkpQ5CPhexVJZsTjIRSVo58qWwxSyEujS3fnO/o8lMSVW/ayWfylp4jsRSGSv4rxo4xlJZK41VsX+wkyv7W626GqUcjl8K8ztffHrNWUvJTA6vy7GuM5MWQ91Y9bZ4eZk52nV8nRaCauNQZQJmlW0sXRnjIKUkkc4t2hBvMpKyjIAKsEaSGSbCKSvLpRSUsSs2F2K9UHXl4Cwec3huZJasWRinWkfHUlmrdYZi39YO/vu3X2adaz630dzuc0+c5RvPXODw8Oya9q9c9RKNiFJlV/W3WQH5xDodAKSNQ5VRd3mxZGXcSknbIJZEEQM0GU1axTwq/hFJZgzl0FK6cVAGppIxlFozG08vOxu6nCwWkH5meMZ6PGwqh0R6oXKYjxCCVr+bjGlYpmNry9NPpHNNGW+AfEzuyv42y0AmK3SjV2ua8xOuIao/TzRVmRPKbhDmy10pJVPRlOU3DZk9cCKJDOORpJXlUgp547B+J2RNRlIlpfiuFe8iMYenz82wpTNAm99tUw6ZFU2fsxu3mTUah2SmeZXDVjPb8MZtnfjMxIF4hb7LtUYbhyqjgofxCgWk43bjMC/uMJfIkM5KSxorP/pkJEk4kVmVW0kZh/A6Mw6HR2YYGgsjpTQM5ypU1WopZhyklDx1bprrtrSzpTPA8JSRCm24lZa/UF+xqY3XXbsJgKno2lyAzawcdvW18Mgf38atO7txOwUOUdjCfD3RnJ9wDfG4HLidooLKIX9BmQgn+eOvHubdX3mWbE4WVEdD/sKuGgGuyq3kW5/K4d1fOcxffuso4WSGVCa3qmD9aimWrTQZTTEWTnL1QDubO/2WcoivwK0E8M//63r+8RevpcXnWrNbqZmVA8BAh6EehBD43M6C71wlSGVy/H/ffMEawVotdOO9GhDwuComRe3bnYikeOD5UcKJDB6Xg9eZrbhVxkV7wHArvTBqBChXoxz8bicOUXwWdSNzYSZOTkqr/mM1x2a1eIr0VlIDZnpbvGzuCPC9F8bI5STxItlKiyGEoCPgWbNbqZmVw3wM41BZ5XDkwiyffPQ0k9Ek/3Tv0vPHy4n+hGtAwOOs2J22XeIeuTBLOJFhe0+Qzz95jm8dvgDkg2rbuoJsaPVx/zPG8tXcHQshjDYc60g5RJIZwskMF2cTtsrxKhqHIu0zZs0Ctja/m4HOAKlsjrFwkniqtLv4joCbqdhaU1mbWznY8bkcFVcOSjHc/+wFhsYiFf1fdrRxqAEBj9PqprkavnxwmDv+4UdIubDff8K23cdOTgLwp6/Zg8fl4D+fPg/klYPDIbjzij6rrcJq745DXte6citdnDUK+sLJjNVevRZuJfvAH1X01up3s9kc8To8HSO2QreSol0rh7Li8zgrHnM4P20YB4/TwUe+f6Ki/8uO/oRrQNDrWlPjvcMjM5wYixStZbCfqKop2NUD7dy6o8uqkO0I5C90d1+RH+yzWuOw3pSDMg6Qn8Vc3Wylhe0z7MpBFR6eGo+QzUmrnmEldAY9ZYk5eLVyAMDncpY9lfXkeKRAIZyfidPidfHz1w/w3RcurWgIVDnQxqEG+N3ONQWklR98ZHrhRDklcVUtQ2fQQ3fIa81pmN8i48ZtnbQH3ATNfkurYf5EuUbnoq0ViJrFvNIBSOXA7TQqj+3GQVU1t/nd9Hf4EQKOjhrV0itJZVW0B9xM62ylsuFzl9+t9GffOMIf2NqBn5+O09/h55qBNqKprJWMUGn0J1wD1qocVGfVkemF2QsqID1guh52mXMGbt/TCyzsuupyOvj56wbYs7F11fuz7ozDbP64vjA6R0fAXTDQptIIYXRZTWYXKodWnwuvy8mmNj/HLhqqphS3UkfAQ8TMwFotyUxOxxxM1hKQfuTEBC8WmYUyGU1x/GLYUgjnZ+L0t/ut7+hRcz54pdHGoQb41xhzUCmpI9NxoskMh85OWc8pt9JmM91Ozfnd2GbcefS3+xds709fs4ev/MYtq96fdRdzmEtYab6xVLaqmUoKr9OxwK0U8rpwmUZqsDvAC6bLqxTFpxTQTHz1rqVkOquVg4nPvfqYwx9/7TD/+L3jC5aHE2ni6awViD4/YyiH3X0tOGyKsdLoT7gGBD3ONTXem4rm3UqfeOQ093z8cSvIOF857O7LTyi77837+bt7rlmwPYdDrKmJWtDrWlftMy7OJhjo8FuuuVoYB4/LUdA+YzaetqbOgTFuViUSlOJW6jDTl9fiWtLKIY/P7VhVWrqUkvFwsmjtgmqtPjQeYS6RJpzI0N/ux+9xMtgdrJpy0HUONSDgWb1bKZPNMW1mroxMxxmdTZCTcGoiyvVbPNYFZcAMWu4ylQNQ1oE0dlp8rrJNGasHLs4l2NjmwyGE2W6kevEGhcdVqBzm4mmr9TYYaciKUgLSKhlhtUHpbE6SyuqYg8LnWl0RXDiZIZXNcWGmsNW9lNJy0Q5dilgzOfrNm709G1s5PDJDNdCfcA1QqazFUlGXY9qWoz48HeM5M2B6etzITIqnsjgE3HVFH2+9dZDrtrSXZ6eXIOg1AuyreT/1yMXZJBvafGw052nXxK3kWuhWavPnjYDq8QOluZVU4eNq01nVPmnlYOB1OwtSjleKSiqZiCQL1o+msqhkpBNjYSuNVbmD925sZXgqXpWbsTUZByHEGSHEc0KIZ4QQB81lnUKI7wohTpi/O2yvf68QYkgI8aIQ4i7b8n3mdoaEEB8W67xRfNDrImPegZXKZNSIN/S3+zkzEbWa652ZNIyDaqfc2+Ljz157RVWGwAe9LrI5WfFioGqQyuSYiCTpa/WxwTQOPVXsq6RYqBwyBW6lwe68cijFraRcZavtr2SNCNXKAVh9ttKkre/Zpdn8Y/tF/8RYxHI75ZWD4Qk4ViSQXW7K8QnfJqW8Vkq53/z7PcBDUspdwEPm3wgh9gL3AlcAdwMfFUKos/pjwNuBXebP3WXYr7pFfZlX46ucMu84rh5os+4whMjXNMTTC8dGVpoW21yIRmcsbMj8jQXKoUZupXnZSmpoDxhDltQtVKnZSrB6t5IKvuo6BwP/KrOV7DVKF2zZcaoWqSvoYeiSYRw8LgfdZleDamYsVcL8vw74tPn408Drbcu/KKVMSilPA0PAjUKIjUCrlPJxafglPmNbZ12iBuSsptZhwgxGX7PZcBc5HYIbBjst45BIVz9YGFxHbbtVAVxfq4+NbcbdWilztcuFp0i2kl05+NxONpr+6FKMg8/txOd2rNqtlEwrt5JWDmAcz0xOkinRC6A8AACjBcbBUA7Xb+0gnMzwzWcvMNDhx+Ew7gQ2tPr491+9gddctbEMe780a/2EJfA/QohDQoi3m8v6pJSjAObvXnN5PzBsW3fEXNZvPp6/fN1iTYNbxcVUydFrBgzjsKs3xN6NrZyZiBpT4DJZa0JVtQitI+VwYVYpBz9XDbTR4nVx2YaWZdYqP0a2knHzkMrkiKezBcYBjIwlKC3mANAZ8KzeraSUQxXclY2Ab4XT4DLZHJdsxZUTYZtysAWl5+LGd+jGwU7ASAB432uvsJ4XQnDb5b1V6fW11mylW6WUF4QQvcB3hRDHlnhtsTiCXGL5wg0YBujtAFu2bCl1X+sGaxrcatxK0RQOAVf0G/Lyqv42BrsCRFNZxsNJEqnqu5XWk3E4ayqwLZ0B/B4nz/35XcusURm8Lqfl+rFaZwQKjcNgd5DHT02W/Hmvpb+SVg6FKJWeSGet70ExvnpohPd98wiPved2OoMeJqNJ2gNupCxUDqoS/rbLe9jZG2L/YActPvdim60oa/qEpZQXzN9jwNeBG4FLpqsI8/eY+fIRYLNt9QHggrl8oMjyYv/vPinlfinl/p6enrXsek0JWNPgSr+YTkRSdAY9tPrc/OGdu3nzLYNs6zFqGU5PREnUoGOmmumwHtp2n56MsqHVt+pWIuVic6efwyOzvPXfD1jtEuYrh5+5eiNv3D9gFcatlI1tPk6bCQylovzrWjkYqGlwy8Udzs/ESaRzPHHKaIY5GUnRFfSwsc3HqE05qJhDq8/NbZf31swwwBqMgxAiKIRoUY+BO4HngfuBt5gvewvwDfPx/cC9QgivEGIbRuD5gOl6CgshbjazlN5sW2ddEvCsISBtmwH9W6/cxVUDbVbO+5nJKPEaKAcr5lCh6XbV5MxElMHuwPIvrDB/8tN7efddl/GDF8f55COnAQrqHABu3dnN37xhYVHjcty8vYtT49GCBoMrJZnRysGOcuEul7Gk1N/jZqfkiUiSrpCXTe1+y5UJeeNQS6OgWMsn3Ac8IoR4FjgAfEtK+W3gA8AdQogTwB3m30gpjwBfBl4Avg28S0qpro7vBP4NI0h9EnhwDftV96wlIG3ccRT6G/s7/LidglMTUTMgXZuYQ7gBlcP8O77TE1G2dYcWeXX18HucvOu2nfS3+/nR8XFgoXJYLS/Z2QXAo0MTJa33P0cuWinTWjkY2N1KS6Gqnh87aRzzyWiK7pCHTe2+BQFpl0PUhfFddcxBSnkKWHDbIqWcBG5fZJ33A+8vsvwgcOVq96XRWEtAeiqaYs+mwiZ5Toegr9XHxdmEVedQTUINmK30P0cu8vtffpZIMsPvvmoXv/uq3czG0kzH0myrA+WguGGwg/8yhzGVyzjs2dBKZ9DDoycn+IV9A0u+9oHnRo3OvX437/z8U9Ygonq4eNUDKzYO5o3TyfEoY3MJJiNJurZ3saHNx0wsTSyVIeBxEU5kaPG51tTOplzoT7gGrCYgPTQW4bOPn2E8nKS7SPvo3hYvY3PJmhiHgMeJEI1jHDLZHO9/4Ci9rV529ob49vMXASw//KCtNUWt2W9mrUD5jIPDIbhlRxePDk0sWdU+l0jzm59/ii8eOMd0LE02J4nrmEMBqhhwObfSXDxt9bX68YkJpmNpukIeq/L5lNnhIJxI14VLCbRxqAmWcijBR//RHwzxf75xhHAyQ2eRvPveFh9j4URNiuCEEIQ8LsINYhy+9dwoZydj/NFdl/Fz1/Vz7GKYyUjSmvq2rbuejIPVYKCgCG6t3Lqjm0tzSU4sMXZy2qypGQsnrWaPimqnS9crKnFhOeUwG09zw6AxO+VLB42M/q6Ql1t3duNyCL5uTmmcM5VDPaA/4RrgcTlwO0VJMYcXL4W5qr+Nd7x8O6+7dtOC5/tavUYqaw1iDqA6s9a3cRgLJ/iH7x7ngw8eY1dviDv3buCWHYb//fFTk5yeiCIEbO6sH7fS7t4WWnwu/G5nwZCmtXL7nl78bicfePDYoupBGYSJSNKqr1F3urq3koHlVlqmv9JcIk1HwMMb92/mwGmjxX530ENPi5e7rtzAVw+NkEhnTeWgjUNT43c7V5ytlM1JhsYi3LStk/e+ek9BXx1Fb6uPuUSmJsoBjHTWeq9z+NSjZ/jwQydo9bv5i9dficMhuLq/jZDXxWMnJzkzGWVTm7+uLnwOh2D/1g4rQ61c9LX6+IM7d/P9Y2Pc/2zRzHGrzmIikrIGTP3Vz1/FX77+yrKqmEYmn8q6nFspQ6vfxVteMojTrHZWhWy/fNNWZuNpvnV41Iw51MexrQ8T1YSUcqc9PBUjmclZg3uKYW8OV4u+N21+NzOx+m7bfXE2QX+7n2//7k9Zy1xOBzdt6+SHx8ZwOkVduZQUf/ozexmbSy7/whJ5663b+M+nznPfw6d43bULmxKomQ925XDlplZevrtxa4zKjVUhvYRbyV7h3t/u59VXbuC/D49areBv3t7Jtu4g//XMeSsgXQ9o5VAjSpkGd/yS0YFxV9/iKZa9NuNQC+XQE/JaHWLrlfFIsmiH1Zft6ubCbILxcJJXX7WhBnu2NDt6Qpb7q5w4HYKr+tuYiBT/3CzlYMYcHMKortbk8a4gW0lVPas6lT+88zLeeuuglfgghOClO7t56uz0ggaLtaQ+TFQTsr07xOOnJokkM0uW3QNW0HDXEsqhtyU/yKcWbpGeFi9PnJ6s+v8thbG5JFu6FsYTfummrezd1MZV/W01r4yuNu1BN9PRNFLKBemTyjjMJTKMziboCHgsl4jGQCmH5BK9leas+d/GRX+wO8if2folAezb2sFnnzhrvq4+LstaOdSI33rlTqaiKav6dSmOXwrT3+5f0oj0tqkNIV0AABIcSURBVNqUg6f6H2t3yMtMLL2mwfWVZjySLFBYCo/LwY3bOpvOMIDRwjuVzVkqdi6R5tYPfJ+Hj48XNOc7PhapyUS8esfjdCDEcsrBbInhX/z7u29rPiutXmIO2jjUiGs3t3PH3j7+9eFTVmn9Yhy/FFnSpQRGp02XeVfnq0EOunLX2FsR1xOpTI6paKomg3vqmc558x2eOTfD+Zk4B89OW6msACcuhcseFF8PCCHwuZZOLrEaJy5RpzLQ4bduXHTMQcNbbx0knMzw1LnpRV8zMh3j5HhkyWA0GFktapxlLdxKaiCOvRVxPaGMlt39psmPDVXBZzWf+Px0nOlYykqfjaWyVWkT3Yj43I4lU1nnu5WKIYSw1INWDhp29RoXfNUmej4f/eEQL/3gD0hlctxgq5RdDOVaqlXMAWA8Unozt2qgsn20ciikI1ioHJ4dMWaSj0zHmI6l2NGTV6xdWjkUxZgGt0TMYV5AejHyxkErh6anO+Qh6HFyZjJW9Pmnzk6zpTPAj979Cu7Y27fs9pQsrUURnFIt9ZqxNBZWykEbBzvzx4Y+ZxqH8zNxpmNpdvfZjYM+dsXwLTMqVA3wWS4L6VV7+tizsZXLazBcqhj1YaKaFCEEW7uCnF2kt/5EJMXWroA18Ws5ekyXSS0Cq+qO3D4bt55QRsseuNdg9fuZjqYYm0twcS5Bi8/F6GwCAWxq9xMw0647dUC6KN5llMNsPI3H6Vj2pm2wO8iDv/Oycu/eqtHKocYMdgc4u4hymIwmrTvylWAphxoEpH1uJy1eV1WUg5SSp5eI0xRjLGy4u/TdbyEqSDodS3PYVA137O0jm5NkcpLOgMc6B4s1fNQYSj25VMwhkabVXx+dVktBG4cas7UryPB0rOiAcjUtaqVs6w7idgoryFhtelq8jC9SUFVOHhma4Oc++lhJBmI8bAxJKmd/ovWAy+mg1ediJpbi8MgMDgF37s0XAnYEPVYKq85WKo7PtZxbqX4K20pBu5VqzGBXgHRWMjqbKGj4FktlSs4Qee01m7h2c3vNqli7W6pTJT1kFgWeHI9y3ZaOZV5tMBZO0qOzbYrSGfQwHUtzfibB9p4QO3vzcYaOgNtSDjpbqTg+t4OJyOKtcOYSGVrK1G69mujbqBqz1Tbi086k6bsvRTk4HaJoU75q0RPyMlEF43BuKlbweyUNDMfDSR1vWIT2gIfpWIqT4xF29YYY6PBbz3UE824lna1UnOUC0rPxdNlmcVQTbRxqzKBlHArjDqoLZiNVpVbLrXTOPFbDUzFOXApz1fu+Y+XnK+YSaV79Tz+2BrqPh4v3VdIY6mBsLsnZySg7e0P43E6rbqUz4GFHT5DOoKchL3DVwO9xLjkiNxxP101LjFLQxqHG9LZ48bkdC2odVBfMRpLy3SEP4URm2cEna8WuHA6enSaTkzx/fq7gNY8NTXB0dM6adqaNw+J0BDycGAuTk1guJTW3oSPg4S0vGeR7v/9yHLqvUlGu3NTGxbkEI9PFE0uMgHTjGVZtHGqMwyHY2hksi1up1liFcBV0LUkpC4zD0VHDKMz/Yj58whjkfmo8yng4SSqbY2Orro4uRkfQQ86c96OK3gY6AjgdghafC7fToYPRS/DSXd0APDa0sPHk8FTMGAnagMdPG4c6YHtP0AqyKiaiSjk0zkmljMNYBY3DWDhJMpOzJt89fc5wJ41Mxwte94hpHE6ORzh20Wh5vrtOiovqjQ5bdpsyDi/d1c3LdnVrtbACdvWG6Gnx8sjQRMHyXE7y7q8+S8Dt5E03bqnR3q0ebRzqgF19LZydihW4YyYjKQIepzVvuhGwguuLtAMpB0o1vHSnMXDmufP5dg+Ks5NRzk3FaPW5ODMZ5dhFQ11cvqG1YvvVyKjstv52v1VA+aYbt/Cpt95Yy91qGIQQ3Lqji8dOTlgjV79ycJh7//UJnjg1xf/5mb1savcvs5X6QxuHOmB3XwgpKVAPk5FkQ6kGgC2dAVwOwcnxxYfWrxUVjH7prvzwG4/LUaAcfmyqhjfu30wineOHL47T0+LVrpFFUMfFnsKqKY2X7OxmIpLixUthzk5GefdXDzMRTvL7d+zmnv0Dtd69VaGNQx2gOq6eGAtbyyajqYar5nU7HWztClTUOJydiuEQcPP2vHG4dUcXY+GkpbweOTFBf7ufV+7pBeCJU5N106+mHlFFk9o4rJ6X7jTiDt8/NmbFuz7xqzfwv2/f1XCV0QptHOqAwa4gLofg+KX8RXUikrLSCRuJHT0hTo6X3600EUny1w8c5YlTk2xs87Oh1UfA40QIeOUeoynhhZk4mWyOR09O8NKd3Zb/PCfhsmVanjczqo5hlzYOq2ZTu5/rtrRz/zMXePj4OAMdfgaLTB1sJLRxqAM8LgfbuoOcuBTmgedG+ezjZwy3UoMpB4AdvSHOTkZJF2kHshY+8/hZ/uXhUxw4PcWWzgBCCLZ0BtjaGbAu/CPTcQ6fnyWcyPCy3d30tngJmj70y7RyWJRdvSH+9g1X87pr+2u9Kw3N667ZxLGLYX744hg/tbunYRWDonGineuc3X0tHDgzxeMnJ4mmsgjRWJlKih09IdJZyfBUjO095bsTfeC5UfZt7eDnrutnz0bjQv9rL9uOlNKq6B2ZjjMeTiIE3LqjGyEE23qCPH9+Tgejl0AIwT37N9d6Nxqe11y9ib/41lHSWclPmemtjYxWDnXCrr4Q4+EksXSWjoAbKRurAE6xo8fIWDq1RtfS2FzCmkd9/FKYobEIr792E79881b2bTUGH71h3wD37N9MX6sPl0MwMh3jkaFxrupvs4bYbOsO4RAsO2ZVo1krPS1ebt3ZjUPALTu0cdCUCRWU/tlrNvGnr9kLwIYGLNpSakEFpXM5SVZVWK2Q0dk4r/z7H/HWTx0gl5N86/AoQsBdV24o+nqnQ7Cp3c/jpyZ5+tyMFRwE+KUbt/AHd15Wk+l4mubj/7xmD/9473XrotWIdivVCbds7+KOvX384Z2XGcPGW73cuG350aD1RpvfTU+L1zIOf/S1wwxPxfjSO27hCwfO8cSpSf72Ddcs2Tr7z+9/gVgqw6NDk7z3P5/je0cvceNg55Lznwc6/Dx2cpLukIefvz6fOnjLji5u2dG16HoaTTnZ1dfCrnWS/KCNQ53QEfTwr2/eb/39sl09NdybtXFZXwtPnZshkc7yrcOjxNNZzkxE+afvneDiXIKAx8lf/dxVRQN2P3xxjG8fuci777qMp8/N8KWDw+zqDfG+n71iyf/5e3fs5u4r53jDvoGGKhzUaOoV/S3SlJ1XX7WBP/n683zshyeJm7UHf3b/ES7OJbhuSztfODDMTdu6eO01m/j4j07ysl3dXD3QDsC/P3qGDa0+fv1l20lksjx8fJw7925YdkjPDYOd3DDYeEpLo6lXdMxBU3Z+5upNeFwOPvKDIVp8Lq7sb+VHx8cJeJx87m03cc1AGx/89jHue/gUf/udF/nFf3mCHx0f58JMnIdPjPPG/QN4XA5afW5rWxqNprrUzbdOCHG3EOJFIcSQEOI9td4fzepp87utOcSvvLyX15v583ddsYGg18X/89N7GJ1N8MFvH+MlO7rY1h3kbZ/6Ce/+6rNIiU6r1GjqgLowDkIIJ/DPwKuBvcCbhBB7a7tXmrXwRvMC/+orN/LaazaxtSvAL99sdKa8aXsXd1+xgZDXxd+84Wq++I6b2T/YwaNDk7x0Z3fBuFSNRlMbhOoiWNOdEOIW4H1SyrvMv98LIKX868XW2b9/vzx48GCV9lCzGo5cmGXvxtaigedUJsdMPGVlICXSWT76gyHuvGIDV/a3VXtXNZqmQQhxSEq5f7nX1UtAuh8Ytv09AtxUo33RlIkrNi1+kfe4HAWpqT63k9+/87Jq7JZGo1kBdeFWAoo1IVkgaYQQbxdCHBRCHBwfH6/Cbmk0Gk1zUi/GYQSwRyEHgAvzXySlvE9KuV9Kub+np3HrADQajabeqRfj8BNglxBimxDCA9wL3F/jfdJoNJqmpS5iDlLKjBDit4DvAE7gk1LKIzXeLY1Go2la6sI4AEgpHwAeqPV+aDQajaZ+3EoajUajqSO0cdBoNBrNArRx0Gg0Gs0C6qJCejUIIcLAi0A3MFGmzbYBs2XaViNsr5zHDur//dbz8av391rPxw7q//3Ww/FT62yVUi5fCyClbMgf4KD9d5m2eV+Z97Het1e2Y9cg77duj18DvNe6PXYN8n5rfvxKXUe7lQr5ZpNtr9zU+/ut5+NX7++1no8d1P/7rffjt4BGdisdlFLuV79rvT+NiD52a0Mfv9Wjj93aWM3xK3WdRlYO9837rSkdfezWhj5+q0cfu7WxmuNX0joNqxw0Go1GUzkaWTloNBqNpkLUnXEQQmwWQvxACHFUCHFECPE75vJOIcR3hRAnzN8d5vIu8/URIcRH5m3rF4UQh83t/E0t3k81WcWxu0MIcUgI8Zz5+5W2be0zlw8JIT4sik3sWWeU+fi9XwgxLISI1Or9VJNyHTshREAI8S0hxDFzOx+o5fuqFmU+974thHjW3M7HzUmbpVPO9KoypWhtBK43H7cAxzFGh/4N8B5z+XuAD5qPg8BLgd8APmLbThdwDugx//40cHut31+dHbvrgE3m4yuB87ZtHQBuwZi18SDw6lq/vwY7fjeb24vU+n010rEDAsBt5mMP8GN97pV87rWavwXwNeDeVe1TrQ/KCg7aN4A7MAreNtoO5IvzXver84zDDcD3bH//CvDRWr+fejx2thNpEvCarzlme+5NwL/U+v00yvGbt7wpjEMljp353D8Bv17r99OIxw9wY6TQ/uJq9qHu3Ep2hBCDGBbySaBPSjkKYP7uXWb1IeByIcSgEMIFvJ7CgULrmlUcu18AnpZSJjHGto7YnhsxlzUNazx+TU25jp0Qoh14LfBQJfe33ijH8RNCfAcYA8LAV1ezH3VrHIQQIQxJ9LtSyrlS15dSTgPvBL6EIU3PAJly7mO9UuqxE0JcAXwQeIdaVORlTZPWVobj17SU69iZN3RfAD4spTxViX2tR8p1/KSUd2EoDS/wyiKrLktdGgchhBvjAH1eSvmf5uJLQoiN5vMbMazikkgpvymlvElKeQuGPDtRqX2uF0o9dkKIgf+/vfsHbSqK4jj+PRgXqfgHrJu4iSKi6CBYcbaLi4uIiRUHURA3UQQddBEtxTrUohX/gIiDUBQUKVhQEVyKtlaQiqDQRay1tjh5HO4Nhr6kNfHFJO3vA4+Um5fLeZfQ8+59L+cB94Gsu4/E5s+ER7XmFX1s61yU0vjNSymPXTfw3t07qh95fUj7u+fuPwlP1NxVSTx1lxziXTHXgGF3by94qxfIxb9zhDW52fpqjq/LgMPA1XSjrS/ljl2ctj8ETrj78/zOcfo6YWZbY59Z/mK8G11a4zcfpTl2ZnaWUKjuWLXjrhdpjZ+ZNRUkkwzQCryrKKhaX3gpcnGlhbCE8RoYiFsr4e6jPsLZfx+wvOAzH4GvwA/CWe+62H4HeBu3iq7YN9JW7tgBp4DJgn0HgOb43hZgEBgBLhN/MDmXt5TH73z8Lv6Kr2dqfXyNMHaEWaoDwwXtB2t9fA00fiuBV7GfIaATyFQSk34hLSIiCXW3rCQiIrWn5CAiIglKDiIikqDkICIiCUoOIiKSoOQgUgVmdsjMsmXsv9rMBqsZk0g5MrUOQGSuMbOMu3fVOg6Rf6HkIFJELH72iFD8bBOhhHIWWAu0A03AF2C/u4+a2VPgBbAN6DWzxYSKrBfMbCPQRShHPQIccPcxM9sM9ABTwLP/d3Qis9Oykkhpa4Bud98AfAeOEH5xutvd8//YzxXsv9Tdd7j7xWn93ASOx37eAKdj+3XgqIfaXyJ1RTMHkdI++Z+6NbeBk4QHqzwJpXBYAIwW7H93egdmtoSQNPpj0w3gXpH2W8DO9A9BpDJKDiKlTa8tMwEMzXCmP1lG31akf5G6oWUlkdJWmVk+EewBXgIr8m1mtjDW0y/J3ceBMTPbHpv2Af3u/g0YN7OW2L43/fBFKqeZg0hpw0DOzK4QqmJ2Ao+BS3FZKAN0EKpfziQHdJnZIuAD0Bbb24AeM5uK/YrUDVVlFSki3q30wN3X1zgUkZrQspKIiCRo5iAiIgmaOYiISIKSg4iIJCg5iIhIgpKDiIgkKDmIiEiCkoOIiCT8Bjq95xcmTML7AAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "sorted_data['inc'][-200:].plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "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": 12, + "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": 13, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 13, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAG+dJREFUeJzt3X+QldWd5/H3B5t0qwOm0cbwS2EqxAqaDIZbwGyys2OYAdxMCdZo0hsndq1UMVE3Y6a2SiG6Za0yVTq1tZmhUnGl4ihq/MGyY8nOypoWNzXsLgGaaEbRMN2JBAnEbreJYqroseG7f9zTw+22+/Zzm+77qz+vqlv3uec+5/Q5PN187/nxnKuIwMzMLIspla6AmZnVDgcNMzPLzEHDzMwyc9AwM7PMHDTMzCwzBw0zM8vMQcPMzDJz0DAzs8wcNMzMLLOGSldgvF1yySUxf/78SlfDzKymHDhw4N2IaBntvLoLGvPnz6ejo6PS1TAzqymSfpHlPA9PmZlZZg4aZmaWmYOGmZll5qBhZmaZOWiYmVlmDhpmNar7/VN8+eE9dJ88Vemq2CTioGFWozbv6mT/4V42v9RZ6arYJFJ392mY1bsr7tlJX/+Zf3795N4jPLn3CI0NUzi06doK1swmA/c0Cri7b7Vg953XcN3i2TRNzf/5Nk2dwprFs9l91zUVrplNBg4aBdzdt1owc3oT0xob6Os/Q2PDFPr6zzCtsYGZ05qGPd8fhmw8eXgKd/et9rz7QR83Lbucry69jKf2HaGnSEAo/DC06frPlLGWVo8UEZWuw7jK5XJR6t5T3e+fYtMLb/KDg7/i1IdnaJo6hVVXfoK7v/TpET+9mVW7oR+GBvjDkA1H0oGIyI12noenKL27b1YLPPdhE8HDU0kp3X2zWuAPQzYRHDSSh792tle2ae1VFayJ2fjxhyEbb57TMDMzz2mYmdn4c9AwM7PMMgUNSX8u6aCk1yU9LalJ0gxJ7ZI603NzwfkbJXVJOiRpVUH6Ekmvpfc2S1JKb5T0bErfK2l+QZ629DM6JbWNX9PNzKxUowYNSXOAPwNyEXEVcB7QCmwAdkXEQmBXeo2kRen9K4HVwHclnZeKewhYDyxMj9UpfR1wIiI+CXwbeDCVNQO4F1gGLAXuLQxOZmZWXlmHpxqA8yU1ABcAx4A1wNb0/lZgbTpeAzwTEX0R8RbQBSyVNAuYHhF7Ij/7/viQPANlbQdWpF7IKqA9Inoj4gTQztlAY2ZmZTZq0IiIXwL/CTgCHAfei4gfAJdGxPF0znFgZsoyB3i7oIijKW1OOh6aPihPRPQD7wEXFynLzMwqIMvwVDP5nsACYDZwoaQ/KZZlmLQokj7WPIV1XC+pQ1JHT09PkaqZmdm5yDI89QfAWxHRExEfAn8L/AvgnTTkRHruTucfBeYV5J9LfjjraDoemj4oTxoCuwjoLVLWIBGxJSJyEZFraWnJ0CQzMxuLLEHjCLBc0gVpnmEF8CawAxhYzdQGPJ+OdwCtaUXUAvIT3vvSENZJSctTOTcPyTNQ1g3Ay2ne40VgpaTm1ONZmdLMzEbk7eAnzqjbiETEXknbgR8D/cArwBbgt4BtktaRDyw3pvMPStoGvJHOvz0iTqfibgUeA84HdqYHwCPAE5K6yPcwWlNZvZLuB/an8+6LiN5zarGZ1T1vBz9xvI2ITSrd75/i3z39Ct/56tXeuK8OeTv4sfM2ImbD8Lcz1jdvBz/xvMutTQr+dsbJwdvBTzz3NGxS8CfQyWNgO/jnbvs8Ny27nJ4P+ipdpbrinoZNCv4EOnn4u3EmloOGTRr+QiKzc+fVU1YXvCrK7Nx49ZSNq2q/WcqroszKw8NTlkm13izlVVFm5eXhKSuq2m+W6n7/FJteeJMfHPwVpz48Q9PUKay68hPc/aVP1+wwlYfarBI8PGXjotqXqtbjqigPtVk18/CUFVUL/ynXy6ooD7VZLXDQsFFV+3/K9bIuf/ed14w41GZWLRw0bFT18p9ytauFXp2Zg4ZZFan2Xp2ZV0+ZmZlXT5mZ2fhz0DAzs8wcNMzMLDMHDbMyqfb9u8yycNAwKxPf6W31wEtuzSaY7/S2euKehtkEq/b9u8xK4aBhNsF8p7fVk1GDhqQrJL1a8Hhf0jclzZDULqkzPTcX5NkoqUvSIUmrCtKXSHotvbdZklJ6o6RnU/peSfML8rSln9EpqW18m29WHgN3ej932+e5adnl9HzQV+kqmY1JSXeESzoP+CWwDLgd6I2IByRtAJoj4i5Ji4CngaXAbOAl4FMRcVrSPuAO4EfAC8DmiNgp6TbgsxHxdUmtwPUR8RVJM4AOIAcEcABYEhEnRqqj7wg3qx7+bpDaMVF3hK8AfhYRvwDWAFtT+lZgbTpeAzwTEX0R8RbQBSyVNAuYHhF7Ih+pHh+SZ6Cs7cCK1AtZBbRHRG8KFO3A6hLrbGYV4hVj9afU1VOt5HsRAJdGxHGAiDguaWZKn0O+JzHgaEr7MB0PTR/I83Yqq1/Se8DFhenD5DGzKuUVY/Urc09D0seA64D/Otqpw6RFkfSx5ims23pJHZI6enp6RqmemU00rxirX6UMT10L/Dgi3kmv30lDTqTn7pR+FJhXkG8ucCylzx0mfVAeSQ3ARUBvkbIGiYgtEZGLiFxLS0sJTTKziVDqijHfLV87Sgka/4azQ1MAO4CB1UxtwPMF6a1pRdQCYCGwLw1lnZS0PM1X3Dwkz0BZNwAvp3mPF4GVkprT6qyVKc2K8B+gVYNSVox57qN2ZFo9JekC8nMLvx0R76W0i4FtwGXAEeDGiOhN790N3AL0A9+MiJ0pPQc8BpwP7AS+EREhqQl4AriafA+jNSJ+nvLcAnwrVeUvIuLRYnX16im457nX+P6+I9y09DI2Xf+ZSlfHbERD5z4GeO6j/LKunvKXMNWRsfwBekmkVVL3+6dG/F50/z6Wl7+EaRIay+SjhwWskny3fO3xhoV1pJQ/QC+JtGrh70WvLQ4adSbrH+DuO68ZcVjArJwe/trZEZFNa6+qYE0sCweNOpP1D9DDAmY2Fp7TmEDVvvTVm+iZWam8emoCeemrmdWKrKunPDw1ATzJbGb1ysNTE6DUpa/VPoxlZjbAQWMClDrJ7HslzKxWeHhqgmRZ+uphLDOrNe5pjEGW4aSHv5Zj09qrWDR7OpvWXjVoKewAbx9t1cTDpLWtXNfPQWMMxms4yfdKWDXxMGltK9f185LbEkzEjpx/+kQHLdOaBg1jDdcrMZso3mm2to3X9fMutxPAO3LacGp9p2D/Xte28bp+3uV2Ang4aXLJOkZc68M6/r2ubeW+fl49VSLvyDl5FAaD4e7or6fVb/69rm3lvH4enjIbIusY8WQe1qn1ITn7KA9PWVWr5uWdWZdCT+ZhnVofkrOx8/CUjausn0BHG/qppFKCwWQb1qmnITkbGw9P2bgabWffWlne6aXQw5vMQ3L1zrvcWlll/QRaK98Y6G+TG95kHpKzPM9p2LjwPMDk4S/vmtzc07Bx4XmAyaOSvTCv2qq8TD0NSR+XtF3STyW9Kel3Jc2Q1C6pMz03F5y/UVKXpEOSVhWkL5H0WnpvsySl9EZJz6b0vZLmF+RpSz+jU1Lb+DXdxlvWT6BZNnM0G45XbVVepolwSVuB3RHxPUkfAy4AvgX0RsQDkjYAzRFxl6RFwNPAUmA28BLwqYg4LWkfcAfwI+AFYHNE7JR0G/DZiPi6pFbg+oj4iqQZQAeQAwI4ACyJiBMj1dUT4Wb1p1YWUNSycbtPQ9J04PeARwAi4p8i4tfAGmBrOm0rsDYdrwGeiYi+iHgL6AKWSpoFTI+IPZGPVI8PyTNQ1nZgReqFrALaI6I3BYp2YPVodTaz+uKvEageWYanfhvoAR6V9Iqk70m6ELg0Io4DpOeZ6fw5wNsF+Y+mtDnpeGj6oDwR0Q+8B1xcpKxBJK2X1CGpo6enJ0OTzKyWeAFF9cgSNBqAzwEPRcTVwG+ADUXO1zBpUSR9rHnOJkRsiYhcRORaWlqKVM3MapVXbVWHLKunjgJHI2Jver2dfNB4R9KsiDiehp66C86fV5B/LnAspc8dJr0wz1FJDcBFQG9K//0heX6YqWVmVldq5d6Zel/hNWpPIyJ+Bbwt6YqUtAJ4A9gBDKxmagOeT8c7gNa0ImoBsBDYl4awTkpanuYrbh6SZ6CsG4CX07zHi8BKSc1pddbKlGZmVpXqfYVX1vs0vgF8P62c+jnwb8kHnG2S1gFHgBsBIuKgpG3kA0s/cHtEnE7l3Ao8BpwP7EwPyE+yPyGpi3wPozWV1SvpfmB/Ou++iOgdY1vNzCbMZNmXy3tPmZmNg1rfl8tbo5uZldFkWeHlbUTMzMbJZNgix8NTZmbm4al6VM3fdmdmk4ODRg2p96V8Zlb9PKdRAybLUj4zq37uadQAb9ZmZtXCQaMGTJalfGbl5nnC0jlo1Ahv1mY2/jxPWDovuTWzScdf6vRRXnJrZjaCUucJPYx1loOGmU06pc4TehjrLC+5NbNJKcuWH17u/lGe0zAzG0Gt71xbCs9pmJmdIy93/ygPT5mZFTEZdq4thYenzMzMw1NmZjb+HDSsqnl9vFl1cdCwqub18WbVxRPhVpW8Pt6sOrmnYVXJ28GbVScHDatKXh9vVp0yBQ1JhyW9JulVSR0pbYakdkmd6bm54PyNkrokHZK0qiB9SSqnS9JmSUrpjZKeTel7Jc0vyNOWfkanpLbxarhVP28Hb1Z9Mt2nIekwkIuIdwvS/hLojYgHJG0AmiPiLkmLgKeBpcBs4CXgUxFxWtI+4A7gR8ALwOaI2CnpNuCzEfF1Sa3A9RHxFUkzgA4gBwRwAFgSESdGqqvv0zAzK1057tNYA2xNx1uBtQXpz0REX0S8BXQBSyXNAqZHxJ7IR6rHh+QZKGs7sCL1QlYB7RHRmwJFO7D6HOpsZmbnIGvQCOAHkg5IWp/SLo2I4wDpeWZKnwO8XZD3aEqbk46Hpg/KExH9wHvAxUXKGkTSekkdkjp6enoyNsnMzEqVdcnt5yPimKSZQLuknxY5V8OkRZH0seY5mxCxBdgC+eGpInUzM7NzkKmnERHH0nM38Bz5+Yp30pAT6bk7nX4UmFeQfS5wLKXPHSZ9UB5JDcBFQG+RsszMrAJGDRqSLpQ0beAYWAm8DuwABlYztQHPp+MdQGtaEbUAWAjsS0NYJyUtT/MVNw/JM1DWDcDLad7jRWClpOa0OmtlSjMzswrIMjx1KfBcWh3bADwVEf9T0n5gm6R1wBHgRoCIOChpG/AG0A/cHhGnU1m3Ao8B5wM70wPgEeAJSV3kexitqaxeSfcD+9N590VE7zm018zMzoG3RjczM2+NbmZm489Bw8zMMnPQMDOzzBw0zMwsMwcNMzPLzEHDzMwyc9AwM7PMHDTMzCwzBw0zM8vMQcPMzDJz0DAzs8wcNMzMLDMHDTMzy8xBw8zMMnPQMDOzzBw0zMwsMwcNMzPLzEHDzMwyc9AwM7PMHDTMzCwzBw0zM8vMQcPMzDLLHDQknSfpFUl/l17PkNQuqTM9Nxecu1FSl6RDklYVpC+R9Fp6b7MkpfRGSc+m9L2S5hfkaUs/o1NS23g02szMxqaUnsYdwJsFrzcAuyJiIbArvUbSIqAVuBJYDXxX0nkpz0PAemBheqxO6euAExHxSeDbwIOprBnAvcAyYClwb2FwMjOz8soUNCTNBb4EfK8geQ2wNR1vBdYWpD8TEX0R8RbQBSyVNAuYHhF7IiKAx4fkGShrO7Ai9UJWAe0R0RsRJ4B2zgYaMzMrs6w9jb8C7gTOFKRdGhHHAdLzzJQ+B3i74LyjKW1OOh6aPihPRPQD7wEXFynLzMwqYNSgIemPgO6IOJCxTA2TFkXSx5qnsI7rJXVI6ujp6clYTTMzK1WWnsbngeskHQaeAb4o6UngnTTkRHruTucfBeYV5J8LHEvpc4dJH5RHUgNwEdBbpKxBImJLROQiItfS0pKhSWZmNhajBo2I2BgRcyNiPvkJ7pcj4k+AHcDAaqY24Pl0vANoTSuiFpCf8N6XhrBOSlqe5ituHpJnoKwb0s8I4EVgpaTmNAG+MqWZmVkFNJxD3geAbZLWAUeAGwEi4qCkbcAbQD9we0ScTnluBR4Dzgd2pgfAI8ATkrrI9zBaU1m9ku4H9qfz7ouI3nOos5mZnQPlP9DXj1wuFx0dHZWuhplZTZF0ICJyo53nO8LNzKpY9/un+PLDe+g+earSVQEcNMzMqtrmXZ3sP9zL5pc6K10V4NzmNMzMbIJccc9O+vrP3hr35N4jPLn3CI0NUzi06dqK1cs9DTOzKrT7zmu4bvFsmqbm/5tumjqFNYtns/uuaypaLwcNM7MqNHN6E9MaG+jrP0NjwxT6+s8wrbGBmdOaKlovD0+ZmVWpdz/o46Zll/PVpZfx1L4j9FTBZLiX3JqZmZfcmpnZ+HPQMDOzzBw0zMwsMwcNMzPLzEHDzMwyc9AwM7PMHDTMzCwzBw0zM8vMQcPMzDJz0DAzs8wcNMzMLDMHDTMzy8xBw8zMMnPQMDOzzBw0zMwsMwcNMzPLbNSgIalJ0j5JP5F0UNJ/TOkzJLVL6kzPzQV5NkrqknRI0qqC9CWSXkvvbZaklN4o6dmUvlfS/II8belndEpqG8/Gm5lZabL0NPqAL0bE7wCLgdWSlgMbgF0RsRDYlV4jaRHQClwJrAa+K+m8VNZDwHpgYXqsTunrgBMR8Ung28CDqawZwL3AMmApcG9hcDIzs/IaNWhE3gfp5dT0CGANsDWlbwXWpuM1wDMR0RcRbwFdwFJJs4DpEbEn8t8x+/iQPANlbQdWpF7IKqA9Inoj4gTQztlAY2ZmZZZpTkPSeZJeBbrJ/ye+F7g0Io4DpOeZ6fQ5wNsF2Y+mtDnpeGj6oDwR0Q+8B1xcpCwzM6uATEEjIk5HxGJgLvlew1VFTtdwRRRJH2uesz9QWi+pQ1JHT09PkaqZmdm5KGn1VET8Gvgh+SGid9KQE+m5O512FJhXkG0ucCylzx0mfVAeSQ3ARUBvkbKG1mtLROQiItfS0lJKk8zMrARZVk+1SPp4Oj4f+APgp8AOYGA1UxvwfDreAbSmFVELyE9470tDWCclLU/zFTcPyTNQ1g3Ay2ne40VgpaTmNAG+MqWZmVkFNGQ4ZxawNa2AmgJsi4i/k7QH2CZpHXAEuBEgIg5K2ga8AfQDt0fE6VTWrcBjwPnAzvQAeAR4QlIX+R5GayqrV9L9wP503n0R0XsuDTYzs7FT/gN9/cjlctHR0VHpapiZ1RRJByIiN9p5viPczMwyc9AwM7PMHDTMzCwzBw0zM8vMQcPMzDJz0DAzs8wcNMzMLDMHDTMzy8xBw8zMMnPQMDOzzBw0zMwsMwcNMzPLzEHDzMwyc9AwM7PMHDTMzCwzBw0zM8vMQcPMzDJz0DAzq4Du90/x5Yf30H3yVKWrUhIHDTOzCti8q5P9h3vZ/FJnpatSkoZKV8DMbDK54p6d9PWf+efXT+49wpN7j9DYMIVDm66tYM2ycU/DzKyMdt95Ddctnk3T1Px/v01Tp7Bm8Wx233VNhWuWjYOGmVkZzZzexLTGBvr6z9DYMIW+/jNMa2xg5rSmSlctEw9PmZmV2bsf9HHTssv56tLLeGrfEXpqaDJcEVH8BGke8DjwCeAMsCUi/lrSDOBZYD5wGPhyRJxIeTYC64DTwJ9FxIspfQnwGHA+8AJwR0SEpMb0M5YA/w/4SkQcTnnagHtSdTZFxNZi9c3lctHR0ZH9X8DMzJB0ICJyo52XZXiqH/j3EfFpYDlwu6RFwAZgV0QsBHal16T3WoErgdXAdyWdl8p6CFgPLEyP1Sl9HXAiIj4JfBt4MJU1A7gXWAYsBe6V1JyhzmZmNgFGDRoRcTwifpyOTwJvAnOANcDAp/6twNp0vAZ4JiL6IuItoAtYKmkWMD0i9kS+e/P4kDwDZW0HVkgSsApoj4je1Itp52ygMTOzMitpIlzSfOBqYC9waUQch3xgAWam0+YAbxdkO5rS5qTjoemD8kREP/AecHGRsszMrAIyBw1JvwX8N+CbEfF+sVOHSYsi6WPNU1i39ZI6JHX09PQUqZqZmZ2LTEFD0lTyAeP7EfG3KfmdNOREeu5O6UeBeQXZ5wLHUvrcYdIH5ZHUAFwE9BYpa5CI2BIRuYjItbS0ZGmSmZmNwahBI80tPAK8GRH/ueCtHUBbOm4Dni9Ib5XUKGkB+QnvfWkI66Sk5anMm4fkGSjrBuDlNO/xIrBSUnOaAF+Z0szMrAKyLLn9ArAbeI38kluAb5Gf19gGXAYcAW6MiN6U527gFvIrr74ZETtTeo6zS253At9IS26bgCfIz5f0Aq0R8fOU55b08wD+IiIeHaW+PcAvMra/mlwCvFvpSowTt6X61Es7wG2ZKJdHxKhDNaMGDSsPSR1Z1kjXArel+tRLO8BtqTRvI2JmZpk5aJiZWWYOGtVjS6UrMI7clupTL+0At6WiPKdhZmaZuadhZmaZOWhMEEl/I6lb0usFab8jaY+k1yT9d0nTU/rHJD2a0n8i6fcL8vxQ0iFJr6bHzGF+3ES3ZZ6k/yXpTUkHJd2R0mdIapfUmZ6bC/JslNSV6r6qIH1JameXpM3pnp1abUvFrk2p7ZB0cTr/A0nfGVJWTV2TUdpS0b+XMbTlDyUdSP/+ByR9saCsil6XEUWEHxPwAH4P+BzwekHafuBfpeNbgPvT8e3Ao+l4JnAAmJJe/xDIVbgts4DPpeNpwD8Ci4C/BDak9A3Ag+l4EfAToBFYAPwMOC+9tw/4XfJbxOwErq3htlTs2oyhHRcCXwC+DnxnSFm1dk2KtaWify9jaMvVwOx0fBXwy2q5LiM93NOYIBHx9+RvVCx0BfD36bgd+ON0vIj89vJERDfwa6Bq1m5HeXY6Lovxaks56zycUtsREb+JiP8NDPq2n1q8JiO1pRqMoS2vRMTA1kgHgSbld9Oo+HUZiYNGeb0OXJeOb+Tsvlo/AdZIalB+65UlDN5z69HU1f4Ple6iauJ2Oi67c2zLgIpfm4ztGEktXpPRVPyawJja8sfAKxHRR5Vdl0IOGuV1C/kvsTpAvuv6Tyn9b8j/UnQAfwX8X/JbsADcFBGfAf5lenytrDUuoInd6bisxqEtUAXXpoR2jFjEMGnVfk2Kqfg1gdLbIulK8l8+96cDScOcVhVLXR00yigifhoRKyNiCfA0+fFxIqI/Iv48IhZHxBrg40Bneu+X6fkk8BQVGhrRxO90XDbj1JaKX5sS2zGSWrwmI6r0NYHS2yJpLvAccHNE/CwlV8V1GY6DRhkNrOSQNIX8957/l/T6AkkXpuM/BPoj4o00XHVJSp8K/BH5Ia5y17scOx2XxXi1pdLXZgztGFaNXpORyqn430upbZH0ceB/ABsj4v8MnFwN12VElZ6Jr9cH+Z7EceBD8p8a1gF3kF9N8Y/AA5y9uXI+cIj8pNlL5HebhPwqkQPAP5CfJPtr0sqdMrflC+S7xv8AvJoe/5r8tyvuIt8r2gXMKMhzN/me1CEKVn2Qn+B/Pb33nYF/g1prS6WvzRjbcZj84owP0u/kohq+Jh9pS6WvyVjaQv7D428Kzn0VmFkN12Wkh+8INzOzzDw8ZWZmmTlomJlZZg4aZmaWmYOGmZll5qBhZmaZOWiYmVlmDhpmZpaZg4aZmWX2/wFIULHVnmtEQwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + " yearly_incidence.plot(style='*')" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "2020 221186\n", + "2021 376290\n", + "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", + "2022 641328\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": 14, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + " yearly_incidence.sort_values()" + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 15, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEICAYAAABPgw/pAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAEtNJREFUeJzt3X2wXHV9x/H3l0Q0cCGigUsbtKFIqZSolGvxYYr3gjooVCvaKqKFjjZ2fKI2HYc6INr6kBZxtA60ExVBUa8KduTB+jDqVfEBScAxYESoIBIURQW9FIHAt3/8zqU314Tk7tnc3fPr+zWTydmzu2e/33t2P3v2d87ZjcxEktR9uwy6AElSfxjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEosXsgHW7ZsWa5YsWKLeXfeeSe77777QpbRd/YwHOxhONhD/61fv/62zNx7e7db0EBfsWIF69at22Le1NQU4+PjC1lG39nDcLCH4WAP/RcRP9yR2znkIkmVMNAlqRIGuiRVwkCXpEoY6JJUie0GekScExE/jYirZ817RER8PiKua/7fa+eWKUnanh3ZQj8XOHrOvFOAL2TmgcAXmsuSpAHabqBn5leAX8yZ/VzgvGb6PODP+1yXJGmeYkd+UzQiVgCXZOYhzeXbM/Phs67/ZWZuddglIlYBqwBGR0cPm5yc3OL66elpRkZGeq1/KNjDcOhqDxs23fHA9OgSuPWuhXnclcuX7pTldnU9zDZsPUxMTKzPzLHt3W6nnymamWuBtQBjY2M59+yrYTsjqxf2MBy62sNJp1z6wPTqlZs5c8PCnMB94wnjO2W5XV0Ps3W1h16Pcrk1In4HoPn/p/0rSZLUi14D/SLgxGb6ROBT/SlHktSrHTls8aPAN4CDIuLmiHgZsAZ4RkRcBzyjuSxJGqDtDtZl5vHbuOqoPtciSWrBM0UlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlWgV6BHxuoi4JiKujoiPRsTD+lWYJGl+eg70iFgOvBYYy8xDgEXAi/pVmCRpftoOuSwGlkTEYmA34Jb2JUmSetFzoGfmJuAdwE3Aj4E7MvNz/SpMkjQ/kZm93TFiL+BC4IXA7cAngAsy8/w5t1sFrAIYHR09bHJycovlTE9PMzIy0lMNw8IehkNXe9iw6Y4HpkeXwK13Lczjrly+dKcst6vrYbZh62FiYmJ9Zo5t73aLWzzG04EbMvNnABHxSeApwBaBnplrgbUAY2NjOT4+vsVCpqammDuva+xhOHS1h5NOufSB6dUrN3PmhjYvyx134wnjO2W5XV0Ps3W1hzZj6DcBT4qI3SIigKOAjf0pS5I0X23G0C8HLgCuBDY0y1rbp7okSfPU6rNdZp4OnN6nWiRJLXimqCRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVaJVoEfEwyPigoj4XkRsjIgn96swSdL8LG55/3cDn8nMF0TErsBufahJktSDngM9IvYEjgBOAsjMe4B7+lOWJGm+IjN7u2PEE4C1wHeBxwPrgZMz8845t1sFrAIYHR09bHJycovlTE9PMzIy0lMNw8IehkNXe9iw6Y4HpkeXwK13DbCYPtiRHlYuX7owxfRo2J5LExMT6zNzbHu3axPoY8A3gadm5uUR8W7gV5l52rbuMzY2luvWrdti3tTUFOPj4z3VMCzsYTh0tYcVp1z6wPTqlZs5c0PbkdDB2pEeblxzzAJV05they5FxA4FepudojcDN2fm5c3lC4A/brE8SVILPQd6Zv4E+FFEHNTMOooy/CJJGoC2n+1eA3y4OcLlB8Bfty9JktSLVoGemd8GtjuuI0na+TxTVJIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklSJ1oEeEYsi4qqIuKQfBUmSetOPLfSTgY19WI4kqYVWgR4R+wHHAO/rTzmSpF5FZvZ+54gLgLcDewD/kJnHbuU2q4BVAKOjo4dNTk5ucf309DQjIyM91zAM7GE4tOlhw6Y7+lxNb0aXwK13DbqKduxh61YuX9rzfScmJtZn5tj2bre41weIiGOBn2bm+ogY39btMnMtsBZgbGwsx8e3vOnU1BRz53WNPQyHNj2cdMql/S2mR6tXbubMDT2/LIeCPWzdjSeM93V5W9NmyOWpwHMi4kZgEjgyIs7vS1WSpHnrOdAz8x8zc7/MXAG8CPhiZr6kb5VJkubF49AlqRJ9GSTKzClgqh/LkiT1xi10SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqRM+BHhGPiogvRcTGiLgmIk7uZ2GSpPlZ3OK+m4HVmXllROwBrI+Iz2fmd/tUmyRpHnreQs/MH2fmlc30r4GNwPJ+FSZJmp/IzPYLiVgBfAU4JDN/Nee6VcAqgNHR0cMmJye3uO/09DQjIyOtaxgkexgObXrYsOmOPlfTm9ElcOtdg66iHXvYupXLl/Z834mJifWZOba927UO9IgYAb4MvDUzP/lgtx0bG8t169ZtMW9qaorx8fFWNQyaPQyHNj2sOOXS/hbTo9UrN3PmhjYjoYNnD1t345pjer5vROxQoLc6yiUiHgJcCHx4e2EuSdq52hzlEsD7gY2Z+c7+lSRJ6kWbLfSnAi8FjoyIbzf/nt2nuiRJ89TzIFFmXgZEH2uRJLXgmaKSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVojM//DfI33xs81uAbSxkz6tXbuak5vEG1S+063l2D9L/R26hS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIq0SrQI+LoiLg2Iq6PiFP6VZQkaf56DvSIWAScBTwLOBg4PiIO7ldhkqT5abOF/ifA9Zn5g8y8B5gEntufsiRJ8xWZ2dsdI14AHJ2ZL28uvxQ4PDNfPed2q4BVzcWDgGvnLGoZcFtPRQwPexgO9jAc7KH/fi8z997ejdr8SHRsZd5vvTtk5lpg7TYXErEuM8da1DFw9jAc7GE42MPgtBlyuRl41KzL+wG3tCtHktSrNoF+BXBgROwfEbsCLwIu6k9ZkqT56nnIJTM3R8Srgc8Ci4BzMvOaHha1zeGYDrGH4WAPw8EeBqTnnaKSpOHimaKSVAkDXZIqYaBLUiU6GegRcWRE7D/oOtroeg9dr39GDX3Yw+ANS/2d2inafFfMJHA7cD9wemZ+ebBVzU/Xe+h6/TNq6MMeBm/Y6h/qLfSI2C8i9pw164XAhZl5BOWPeHxEPHkw1e2Yrvcwn/ojYmtnDw+Frq8HsIdhMOz1D2WgR8RjI+LTwGXAP0XEzJd+/QbYrZn+OOW7Fp40jEHS9R56qT+H8ONe19cD2MMw6Er9QxPoEbH7rItPAG7OzBXAF4F3NPN/AdwdEXtk5i+A7wOjwIoFLHWbut7Dg9T/JX67/pGm/usYkvpndH09gD0wBD10sf6BBnpE7BUR50bEFcCaiNi7eWd7HPC1ZqvvIuD2iDiG8sfaA1jZLOI6yrei3TOI+qH7PTxI/Stn1f+prdT/uGYR3x9k/TO6vh7AHppFDOvroRP1D3oL/QhgM/Bsyrc3vgHYk1LXvrM+wp8HvBj4FvBryo9qkJnfAI4EfrWwZW+h6z1sq/5FdKP+GV1fD2APw9BDt+vPzJ3+jxIOrwC+TPlu9GXN/I8Dr22m9wfWNNc/kTJWtai5bgT4WbOc5cBG4NXAB4Czgd3soe76a+rDHgbfQ9fr39a/hdpCPxZ4DvBm4MnAvzbzPw88pZn+EfBV4FmZeQXl3XECIDOngcuBJ2bmJuCllHGqnwCnZub/2EP19dfUhz0Mvoeu179VbX7gYgszRzlExBMpH0W+ClyamXcDfwD8IDO/GBE3AGdExDOB9cDzImJZZt4WEdcBd0bEo4H3AC+JiH0o37X+c8rHGzJzHbCuX7XX0kPX66+tj6YXe/D1sGD6soU+6w93BHAO5VCepwNvb25yP/D9iFiSmTdQPp48jjL2dAvlWE6A+ygfYXYBLqR8heUJwGHA2sy8vx/1bqOHRU0PT6N8bOpUD01dGRHjdHQdNH3s2fU+IuKRzf9PAc7taA/7RMQjI2KMMl7cqR4iYtcuv5571utYDeXYy78FPgL8DfAQ4O+AVzXX7wV8BziU8sdZA6xorjuW8odZ1kxvAJZSdiZ8Gth11uPs0u9xplnL3h14OWVFrabs/OhaD3sAl1K+jx7gdV2qf9Zz6UTgC5STNDrVB/93xvVfUMZZpyi/n9u159LuwEmUYYc7gGO61AMlg14JXAz8O/AY4OSu1N+Pfz1toUfEvsAlwDjwIcpOg+MoY1GbATLzl8CngNc2T5B9gMc2i/gK8DTgnsy8BHg/cAFwFmVr4N6Zx8qd9w6+OyVAjgTeCzwTeD5l58f9XeihsQR4KHBARCwDDqBsUXSi/oh4CHAN8ALgjMx8fnPVobMef6j7yMyMiKXAXwLvyszxzLyWshXXiR4i4jGUIYmjgFOBTcBNlB2DXXk+vYryen4X5feNj2su39eR+tvr8Z1wCXD4rMsnUXYqnAh8a9b83wVuaaZfRTk1dq/m/hcDj55122UL/W4GPHzW9OspK/qEjvVwInAGcBrwMuDPgCu6Un/zuJ8ETpgz74XA5V3pg7Jl+M/N9MwW+3Fd6YES2g+ddfkcypvsczvUw8XAXzXTLwNe0zyPOvN6bvuv1zH03wDfag64h/Ix5vDMPI+ypbgvQGbeAlwTEYdn5lnA9cDHKD8wfVlm3jSzwMy8rcdaepaZt0fEnhFxLmXIZRllpR4QEaPD3MOsv/0uwH9TtnKPzMyLgf2Hvf45zgFOj4gzI2IqIt4IfJPym7X7NLUNex+3AX8aEScA6yPiPMqW4R82n5yGuofMvC/LzsKZfQBB+cKpiynroQuv6f8CToyITwBvAh4PXN3Uv3cH6m+v5TvizJbIecDJzfSHgH9pph8BvI/mXY8yxnUI8LBBv5PN6eOVlI+Waylj6l8H3kh5Ug91D8AnKFtXSynH1J5KeRKf1rF18FnKSRyPorzATga+0aH1cGBT779RtvheDLyTcqzy6ylvvEPdw5x+rgKOa6bP78prunn9nkPZOHsTcDpwbfPc6tQ66OVfq6NcMjMjYj9gX8q7I80fkYi4hDJGvSibd73MvDczr87M37R53H7LzLOz7Ok+mzLu+R+Uj2AXMcQ9RMQIZcvwvcBnKHvpDweOB/aKiIsZ4vrneF5mvi0zfwS8jXJY2XvowHpo3EQ53XtxlrHaS5p5l1B2+nahByJiJhMuo+xUBHhruaoTz6eDgaksW9cfoGzsfJTuPI9a6cdx6IdS9gjfEBEvp3x0eQNl7Op7mXllHx5jofycslPx1Mz8YES8BLgmM68acF3bspkSFvdStszvA96QmRuAv+9A/Q/ILU/EuJ0S6Kdl5ke60Edm3h0Ra4C3NLN2oRwksCYzr+pCD1B2+EXEzLcHXt/M2wi8vhlO+u6w9hARiyj5cwTwQeBOykEOr8zMG7qyDtpo/QMXEfE14PeBGynHb745M7/TvrSF0RydcBTlI/LBlGGXszLz3ge94xBqTn44DpjMzJ8Mup75iIiHAkdTzrj7I8phZ2dn5uaBFjZPEfFWyuvhUOBzlB88+OVgq5q/iLgWeGNmfmzmPJNB17QjIuIAymv4Hsp6+E/gLVnO7Kxeq0BvDjk7HbgBOD+bnSpdEhGLKd/pcDelh8599Gq2TO7vyotuWyLiFZRDRj/UxfUwIyIOAn7YxR5mnST4BMqO9s1de141GzYHAl/PzLsGXc9C6tRP0EmStm3QX58rSeoTA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRV4n8Bc3dQMZDbe6kAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + " yearly_incidence.hist(xrot=20)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], "metadata": { "kernelspec": { "display_name": "Python 3", @@ -16,10 +2302,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.3" + "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 } -