{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Incidence de la varicelle " ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import isoweek" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "data_url = \"https://www.sentiweb.fr/datasets/all/inc-7-PAY.csv\"\n", "\n", "data_file = \"inc-7-PAY.csv\"\n", "\n", "import os\n", "import urllib.request\n", "if not os.path.exists(data_file):\n", " urllib.request.urlretrieve(data_url, data_file)" ] }, { "cell_type": "code", "execution_count": 10, "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
0202443727017254677417FRFrance
12024427265912604058426FRFrance
2202441720353813689315FRFrance
3202440721257253525315FRFrance
42024397289813334463426FRFrance
5202438775101513102FRFrance
62024377916281804102FRFrance
7202436722358703600315FRFrance
8202435716202852955204FRFrance
9202434725606224498417FRFrance
10202433719715363406315FRFrance
1120243274399194468547311FRFrance
1220243174500221367877410FRFrance
13202430770044278973011715FRFrance
1420242979270630312237141018FRFrance
1520242879364649812230141018FRFrance
16202427710247709013404151020FRFrance
172024267143681039918337221628FRFrance
18202425711174803914309171222FRFrance
19202424712621935715885191424FRFrance
202024237146571133917975221727FRFrance
21202422711628836114895171222FRFrance
2220242179701685112551151119FRFrance
232024207136611020917113201525FRFrance
2420241971008364131375315921FRFrance
25202418713438951417362201426FRFrance
262024177153031121919387231729FRFrance
272024167181381354022736272034FRFrance
282024157249291731532543372648FRFrance
292024147161811254419818241929FRFrance
.................................
17391991267176081130423912312042FRFrance
17401991257161691070021638281838FRFrance
17411991247161711007122271281739FRFrance
1742199123711947767116223211329FRFrance
1743199122715452995320951271737FRFrance
1744199121714903897520831261636FRFrance
17451991207190531274225364342345FRFrance
17461991197167391124622232291939FRFrance
17471991187213851388228888382551FRFrance
1748199117713462887718047241632FRFrance
17491991167148571006819646261834FRFrance
1750199115713975978118169251832FRFrance
1751199114712265768416846221430FRFrance
175219911379567604113093171123FRFrance
1753199112710864733114397191325FRFrance
17541991117155741118419964271935FRFrance
17551991107166431137221914292038FRFrance
1756199109713741878018702241533FRFrance
1757199108713289881317765231531FRFrance
1758199107712337807716597221529FRFrance
1759199106710877701314741191226FRFrance
1760199105710442654414340181125FRFrance
17611991047791345631126314820FRFrance
17621991037153871048420290271836FRFrance
17631991027162771104621508292038FRFrance
17641991017155651027120859271836FRFrance
17651990527193751329525455342345FRFrance
17661990517190801380724353342543FRFrance
1767199050711079666015498201228FRFrance
17681990497114302610205FRFrance
\n", "

1769 rows × 10 columns

\n", "
" ], "text/plain": [ " week indicator inc inc_low inc_up inc100 inc100_low \\\n", "0 202443 7 2701 725 4677 4 1 \n", "1 202442 7 2659 1260 4058 4 2 \n", "2 202441 7 2035 381 3689 3 1 \n", "3 202440 7 2125 725 3525 3 1 \n", "4 202439 7 2898 1333 4463 4 2 \n", "5 202438 7 751 0 1513 1 0 \n", "6 202437 7 916 28 1804 1 0 \n", "7 202436 7 2235 870 3600 3 1 \n", "8 202435 7 1620 285 2955 2 0 \n", "9 202434 7 2560 622 4498 4 1 \n", "10 202433 7 1971 536 3406 3 1 \n", "11 202432 7 4399 1944 6854 7 3 \n", "12 202431 7 4500 2213 6787 7 4 \n", "13 202430 7 7004 4278 9730 11 7 \n", "14 202429 7 9270 6303 12237 14 10 \n", "15 202428 7 9364 6498 12230 14 10 \n", "16 202427 7 10247 7090 13404 15 10 \n", "17 202426 7 14368 10399 18337 22 16 \n", "18 202425 7 11174 8039 14309 17 12 \n", "19 202424 7 12621 9357 15885 19 14 \n", "20 202423 7 14657 11339 17975 22 17 \n", "21 202422 7 11628 8361 14895 17 12 \n", "22 202421 7 9701 6851 12551 15 11 \n", "23 202420 7 13661 10209 17113 20 15 \n", "24 202419 7 10083 6413 13753 15 9 \n", "25 202418 7 13438 9514 17362 20 14 \n", "26 202417 7 15303 11219 19387 23 17 \n", "27 202416 7 18138 13540 22736 27 20 \n", "28 202415 7 24929 17315 32543 37 26 \n", "29 202414 7 16181 12544 19818 24 19 \n", "... ... ... ... ... ... ... ... \n", "1739 199126 7 17608 11304 23912 31 20 \n", "1740 199125 7 16169 10700 21638 28 18 \n", "1741 199124 7 16171 10071 22271 28 17 \n", "1742 199123 7 11947 7671 16223 21 13 \n", "1743 199122 7 15452 9953 20951 27 17 \n", "1744 199121 7 14903 8975 20831 26 16 \n", "1745 199120 7 19053 12742 25364 34 23 \n", "1746 199119 7 16739 11246 22232 29 19 \n", "1747 199118 7 21385 13882 28888 38 25 \n", "1748 199117 7 13462 8877 18047 24 16 \n", "1749 199116 7 14857 10068 19646 26 18 \n", "1750 199115 7 13975 9781 18169 25 18 \n", "1751 199114 7 12265 7684 16846 22 14 \n", "1752 199113 7 9567 6041 13093 17 11 \n", "1753 199112 7 10864 7331 14397 19 13 \n", "1754 199111 7 15574 11184 19964 27 19 \n", "1755 199110 7 16643 11372 21914 29 20 \n", "1756 199109 7 13741 8780 18702 24 15 \n", "1757 199108 7 13289 8813 17765 23 15 \n", "1758 199107 7 12337 8077 16597 22 15 \n", "1759 199106 7 10877 7013 14741 19 12 \n", "1760 199105 7 10442 6544 14340 18 11 \n", "1761 199104 7 7913 4563 11263 14 8 \n", "1762 199103 7 15387 10484 20290 27 18 \n", "1763 199102 7 16277 11046 21508 29 20 \n", "1764 199101 7 15565 10271 20859 27 18 \n", "1765 199052 7 19375 13295 25455 34 23 \n", "1766 199051 7 19080 13807 24353 34 25 \n", "1767 199050 7 11079 6660 15498 20 12 \n", "1768 199049 7 1143 0 2610 2 0 \n", "\n", " inc100_up geo_insee geo_name \n", "0 7 FR France \n", "1 6 FR France \n", "2 5 FR France \n", "3 5 FR France \n", "4 6 FR France \n", "5 2 FR France \n", "6 2 FR France \n", "7 5 FR France \n", "8 4 FR France \n", "9 7 FR France \n", "10 5 FR France \n", "11 11 FR France \n", "12 10 FR France \n", "13 15 FR France \n", "14 18 FR France \n", "15 18 FR France \n", "16 20 FR France \n", "17 28 FR France \n", "18 22 FR France \n", "19 24 FR France \n", "20 27 FR France \n", "21 22 FR France \n", "22 19 FR France \n", "23 25 FR France \n", "24 21 FR France \n", "25 26 FR France \n", "26 29 FR France \n", "27 34 FR France \n", "28 48 FR France \n", "29 29 FR France \n", "... ... ... ... \n", "1739 42 FR France \n", "1740 38 FR France \n", "1741 39 FR France \n", "1742 29 FR France \n", "1743 37 FR France \n", "1744 36 FR France \n", "1745 45 FR France \n", "1746 39 FR France \n", "1747 51 FR France \n", "1748 32 FR France \n", "1749 34 FR France \n", "1750 32 FR France \n", "1751 30 FR France \n", "1752 23 FR France \n", "1753 25 FR France \n", "1754 35 FR France \n", "1755 38 FR France \n", "1756 33 FR France \n", "1757 31 FR France \n", "1758 29 FR France \n", "1759 26 FR France \n", "1760 25 FR France \n", "1761 20 FR France \n", "1762 36 FR France \n", "1763 38 FR France \n", "1764 36 FR France \n", "1765 45 FR France \n", "1766 43 FR France \n", "1767 28 FR France \n", "1768 5 FR France \n", "\n", "[1769 rows x 10 columns]" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data_csv_path = \"inc-7-PAY.csv\"\n", "raw_data = pd.read_csv(data_csv_path, skiprows=1)\n", "raw_data" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "raw_data[raw_data.isnull().any(axis=1)]\n", "data = raw_data.dropna().copy()" ] }, { "cell_type": "code", "execution_count": 14, "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": 16, "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" ] }, { "cell_type": "code", "execution_count": 17, "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": 18, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXmYXUWZP/55773dnX1PICSBRBIghCVAjCCLQATC4oAKGlxgFIRhQJ2vywjqjKhEQUZQ/AkOi8MyyCLIgGyCIRCBQAiQhYQkNCSQfV86S3e6+9bvj3PqnDq1nTp36ZtO1+d5+ul769apqlOnTr317sQYg4eHh4eHh4hcrQfg4eHh4bHnwRMHDw8PDw8Fnjh4eHh4eCjwxMHDw8PDQ4EnDh4eHh4eCjxx8PDw8PBQ4ImDh4eHh4cCTxw8PDw8PBR44uDh4eHhoaBQ6wGUikGDBrGRI0fWehgeHh4enQpvvvnmBsbY4LR6nZY4jBw5ErNnz671MDw8PDw6FYjoQ5d6Xqzk4eHh4aHAEwcPDw8PDwWeOHh4eHh4KPDEwcPDw8NDgScOHh4eHh4KPHHw8PDw8FDgiYOHh4eHhwJPHDw6Bd78cDMWrtpW62F4eHQZdFonOI+uhc/f9ioAYNn1Z9d4JB4eXQOec/Dw8PDwUOCJg4eHh4eHAk8cPDw8PDwUpBIHIupGRLOIaC4RLSCin4bl1xLRSiKaE/6dJVxzDRE1EtFiIjpDKD+GiOaHv91CRBSWNxDRQ2H560Q0svK36uHh4eHhChfOoQXAqYyxIwGMBzCZiI4Nf7uZMTY+/HsaAIjoUABTAIwDMBnArUSUD+vfBuAyAGPCv8lh+SUANjPGRgO4GcAN5d+ah4eHh0epSCUOLMD28Gtd+Mcsl5wL4EHGWAtjbCmARgATiWgogD6MsZmMMQbgXgDnCdfcE35+BMAkzlV4eHh4eHQ8nHQORJQnojkA1gF4njH2evjTVUQ0j4j+SET9w7JhAJYLl68Iy4aFn+XyxDWMsTYAWwEM1IzjMiKaTUSz169f73SDHh4eHh7Z4UQcGGPtjLHxAIYj4AIOQyAiOhCBqGk1gF+H1XUnfmYpt10jj+N2xtgExtiEwYNTExl5dGEwxvDXuavQXrQxuR4eHiZkslZijG0B8CKAyYyxtSHRKAK4A8DEsNoKACOEy4YDWBWWD9eUJ64hogKAvgA2ZboTDw8Bj7y5At984G3c/eqyWg/Fw6NTwsVaaTAR9Qs/dwfwaQCLQh0Cx2cBvBN+fgLAlNACaRQCxfMsxthqAE1EdGyoT7gIwOPCNReHn88H8EKol/DwKAkbtu8GAKzb1lzjkXh4dE64hM8YCuCe0OIoB+BhxtiTRHQfEY1HIP5ZBuByAGCMLSCihwEsBNAG4ErGWHvY1hUA7gbQHcAz4R8A3AXgPiJqRMAxTKnAvXl0YXhzBg+P8pBKHBhj8wAcpSn/quWaqQCmaspnAzhMU94M4IK0sXh4ZIVnPz08SoP3kPboVGgvMnzu1lfw0hK7tZpnHDw8yoMnDh6dCpt27MZbH23Bdx+eU+uheHjs1fDEwWOvhrdr8PAoDZ44eHQqMEctgldIe3iUB08cPPZqVItxWNfUjCvvfws7Wtqq04GHR43hiYPHXgmqskr6pueW4Kn5q/HE3FXplT08OiE8cfDopHDb/L3GwcOjNHji4NFJYd/2vc7Bw6M8eOLgsVfDGyt5eJQGTxw8Oik8a+DhUU144uCxV8PV9NXDwyMJTxz2QLQXGV5p3FDrYXh4eHRheOKwB+K2Fxvx5Ttfx4yU+EEeHh4e1YInDnsgPtiwAwCwrqmlxiPpvPApyD08yoMnDh57Nby1kodHafDEwaNzwXGzrzbf4IlOZfBq4wa8+eHmWg/DQwOXTHAeHh4eVcGX7nwdALDs+rNrPBIPGZ5z8NgrUW2Vg1dpeOztSCUORNSNiGYR0VwiWkBEPw3LBxDR80T0Xvi/v3DNNUTUSESLiegMofwYIpof/nYLhVpDImogoofC8teJaGTlb9Vjb0KtN+eOFCsdf/0L+Oytr3Rchx4ecOMcWgCcyhg7EsB4AJOJ6FgAVwOYxhgbA2Ba+B1EdCiAKQDGAZgM4FYiyodt3QbgMgBjwr/JYfklADYzxkYDuBnADRW4t4rhhUVr8d7aploPw0OA6+a8NyT7WbllF97+aEuth+HRxZBKHFiA7eHXuvCPATgXwD1h+T0Azgs/nwvgQcZYC2NsKYBGABOJaCiAPoyxmSx4Y++VruFtPQJgEu1Btohfv3s2Trt5Rq2H4ZEB1V48e87q9PCoDpx0DkSUJ6I5ANYBeJ4x9jqAfRhjqwEg/D8krD4MwHLh8hVh2bDws1yeuIYx1gZgK4CBpdyQR9eA6+ZcLb5hL2BIPDyscCIOjLF2xth4AMMRcAGHWarrXltmKbddk2yY6DIimk1Es9ev997DXRGue/IexHh6eHRKZLJWYoxtAfAiAl3B2lBUhPD/urDaCgAjhMuGA1gVlg/XlCeuIaICgL4ANmn6v50xNoExNmHw4MFZhu5RA2zY3oJdu9trOobOcsKfvmgdNm73HvEeew5crJUGE1G/8HN3AJ8GsAjAEwAuDqtdDODx8PMTAKaEFkijECieZ4WipyYiOjbUJ1wkXcPbOh/AC2wP1CRe+8QCtLYXaz2MToMJ1/3d2crm6kfn4eE3lqfWc+UHOopxqEQ3LW3t+Nrdb+Ard82qQGseHpWBC+cwFMB0IpoH4A0EOocnAVwP4DQieg/AaeF3MMYWAHgYwEIAzwK4kjHGj49XALgTgZL6fQDPhOV3ARhIRI0AvoPQ8mlPw92vLsMLi9alV7RgxeadWLutuUIj2vOxaI2bldeDbyzHvz86L7Ve1hNDtUN2l9M6Y0H0XX4MalznLeI89hykekgzxuYBOEpTvhHAJMM1UwFM1ZTPBqDoKxhjzQAucBhvzVEslrfZnHDDdADeI7RUuPKTnUHj8PDs5fjBo/Nxw+cPBwC0lbm2PDwqCe8hnRH+9a0ttuzanal+tYWT5RChFZt3AQA+2rQTQOfRj3h0DXjikBH+Ba4tCrlgyaZuyh2kdChnOfAReobBY0+EJw4Z4dNO7hlwfQp79NMKCViWA8eqLbvwvT/Pxe42bxjhUV144pARnnOoNdweQEfpHMrph19rM8x7f/32xPcf/987eOTNFT5LoEfV4YlDRnjaUFtwEYzrplxtYl6WWCm8iaJhkPNWbMGkX7+UKON187nOoHL36MzwxCEjOuKV/L+3V3ZAL50TztZKnWDvzIWDNOkcPty4UymLiGMnuD+Pzg1PHDKiI15Kr6A0w1Xn01Hiv0qIlUycgw5cBOXDg3hUG544ZITXOdQWe9P88/09yz1xQuKlSh7VhicOGbEX7U2dEnuTWIkisZL7qiqGRkq5znCDHp0anjhkxB4Y8qlLIctGGsBc/5q/zMcnfvH36PvutiKemrc60zOupkJah5hz8MTBo7rwxMGjU+FbD7ztVI8ctAEPzPoIa7fFkVB/8/cluPJPb2H64vLiZ4nY1tyKXz+3GG1SwMaN21vwq2cXAwCamtu01+pIBvMKaY8OgicOGeEZh9rigw07AGRI9pPhea3eGgRE3LKz1fmatGH88ul38bsXGvH0O2sS5fw+AODxOavky4zwnINHR8ETh4zwHtKlY+WWXbjtxfczi+Zc6y/ftBOLwyiwr76/IfP44v4y1E35fWeYz6K9mOQcSt3aeX9eIe1RbXji4NFhuPSe2bjh2UVYvmlXpuvaHW17T/zVdJzxmyDX95PzVgPoACc4FuT5eG+tPtx2JAaSyEEppqjNre3RXHhTVo9qwxOHjPBiJRW7drfj508uxM7detk5B/89q1I5ayhrV2JSCazYvBN3v7oMX7/nDe3vUR5caS932dvlKof8x7MRF+U5B49qwxOHjPDEQcX/vLoUd728FLfP+MCpftYpzLrZb9sV6wyyiAFL2W/TWjeJxFz60l05d8VWAHufzuFvC9Zg5vsbaz0MDwGpyX48kvC0QUVbezAraZt4qdtZVs5hdw1SuZqso2IdQfliJRF7G3G4/L43AfgkWHsSPOeQEd7PQUW1pyQr5yDWr/bYXm20K77jcBfJ8r1ra/fYG5FKHIhoBBFNJ6J3iWgBEX07LL+WiFYS0Zzw7yzhmmuIqJGIFhPRGUL5MUQ0P/ztFgqPT0TUQEQPheWvE9HIyt9qZeBJgxnukVKz6hyycQKck+kIcDGP6SDf0hqMXeYsyj35l2s19/767WhqdjfZ9eh6cOEc2gB8lzE2FsCxAK4kokPD325mjI0P/54GgPC3KQDGAZgM4FYiyof1bwNwGYAx4d/ksPwSAJsZY6MB3AzghvJvrUrw1KFkcFFKtXUOO1tjxXitH9e0RYFDXSUU0pXEpF+/hC/+92tV7MGjsyOVODDGVjPG3go/NwF4F8AwyyXnAniQMdbCGFsKoBHARCIaCqAPY2wmC46O9wI4T7jmnvDzIwAmUblC2SrB+zmUjji5TbbrsnICO1ras3UgoZQnnLZYK63sroS4bOHqbdrytz7ajDNunoFdu8ubR4/OjUw6h1DccxSA18Oiq4hoHhH9kYj6h2XDACwXLlsRlg0LP8vliWsYY20AtgIYmGVsHrWDM8Eskdy7KKQfezteWqJJbYeF7k45y8g/ly9Wqh5+9teFWLy2yUg8PDoOH23cia27aiP+cyYORNQLwKMA/o0xtg2BiOhAAOMBrAbwa15VczmzlNuukcdwGRHNJqLZ69fXJk1iV9BH/3n2ciwVwjs4o0rMXpqOYsP2Fvy/h+ZG38vlHEpBKuegWCvp6017d61Tf9U0jOA6nrr8Hsm8dymcdON0fOZ3L9ekbyfiQER1CAjD/YyxvwAAY2wtY6ydMVYEcAeAiWH1FQBGCJcPB7AqLB+uKU9cQ0QFAH0BbJLHwRi7nTE2gTE2YfDgwW53WGHs7bSBMYbvPzIP5/5/1VyQGcNnaMpEBa8sdio5dEYV90K5aRPncMk9s53aq+Y65PP55LzVWLFZzUbn0bH4aFNtnoGLtRIBuAvAu4yxm4TyoUK1zwJ4J/z8BIApoQXSKASK51mMsdUAmojo2LDNiwA8Llxzcfj5fAAvsD3UZnTPHFXpuGPGB/jF0+9G37kIZ5shUmg5KFXnkBX3zvww+txhOqIUwuLKOewJ4Gvg9hkf4Ct3vp5SuzQs37QTE677e3pFj5rBxQnueABfBTCfiOaEZT8EcCERjUdwiFkG4HIAYIwtIKKHASxEYOl0JWOM8/lXALgbQHcAz4R/QEB87iOiRgQcw5Tybis7/r5wLY4fPQjd6/PWenubQnpqSBh+eNZYAEFOAwAo7EHxGfYGglzp2azmnIjhxbdUSd79yJsrsGF7S3pFj5ohlTgwxl6Gfm0/bblmKoCpmvLZAA7TlDcDuCBtLNXCkrVNuPTe2fjsUcNw8xfHW+vuDRuVDZw41OXdbRXcs7OVZsqadoX1FF6l5yUfEtIUzKpCuvwRVAutgpiuR535sNTaXsQH63fg4H17Z+5jT+acPAJ4D2nEdvQLV6VbZ1TqlZSTv2j7qgEl4qEnRGVk47omvPXR5tRrq2HOCXQOgpyukK5sf1XlHASnQxsnPfWpd3HGb2ZgeQkycZdkTB61hScOABoKwTQ4eYxW6K1sdbDdr8WeyDmH+kK8ND590wx87tZXazCaAOXMQynXuhDlrJtbpd12qrk2RMthG0fEDwwbd+zO3IfnHPZ8eOKA+GVoFd6Ksf/xLM77/StV69MlOFy+Bm9QazvXOVRvaWSlrx3mq5Bhw5fFSmmPSt5ka8kNpRE/11ng95Q9r3fXiC014bq/47YX36/1MEqGJw6IXxZRDryrtR1zlm9R61aoz1YLceDK4FqcrviLnncUiq/Z2ow/zfrIqS6/n6xKfV39asxNlnFl3Q/l4Za7jsohLmnXinNrq8qXSCniz67AOWzY3oIbnl1U62GUDE8cEHMOupOjvPBLfSlb2trx9PzV0XcbcSi3L1d8789zjb+5vvBfu/sNrG8KrE7SXng+v5W4L7GNWuwz8i2kcR2V1znYJ/Gpeaux0WANVKllFXMO2a/dQ6PjeAjwxAGxQlp3WJZTWpaqJL7x2cX41/vfir63tpnbYdL/auGRN1coZeLt7drdjnVNzdY2Ngvy5iwn0iwo75ScJdmP+wCL0o6Ydm/yMMolkLbLN25vwZV/eguX3qt3qEsXK8U3Y6vL71meC4+9Az7ZD2JRiu40c9KN0xPfS30NVm5JEhmbzoG/kLWwVhJ7nHL7zCgkdS2hm4ZqHjxdZj2rnL3SfiO27rmxw6ot+lzdlVhVre1FrNwctG+iDbfPeB+LVjfhJo15+N6WrGhvhOccIBKH9Lql7NfvrNyKZ95ZkyjjVkHaPqT/HQl+f0SUmTDM/CBO82gjbFnn8KUlpcfRqtYclntYLteZspzrTfO/eusuMMacdA4/f3IhVm1tDtvT1/rF04vwl7dXan/ztGHPhycOEHQOLsShhPZ/9teFSplN58BqSB3K2XRmLVXCYZWNtduaO1yp57JvyZyDjusUTaP3JMGL7hm/v347jvvlC/hvSx7wpRt2YOTVT+GlJesTz7qUe+tKtMGk+9nT4YkDRJ1D+pItRdTT3KZGCXVSSNdgS6mUJMvWTpb7+tFj89MrWR5bKffjcolLu9PeXWesX/Y8a67ftbsdW3cKBMnQh658RSgieqVxQ3I6hbpvLAsIwl/nrkLPhlginSZim3Dd8wqnvLeLlcR94tZOas7qiQPiB1mt5apLmuIilugUnsGGSdMNPQqfkeG+XJwFRZQj2s+yXyl+DhnrV0MhffYt/8CRP3su9VqXg0lapz0Ez+m0TH0btu/GphIc5TozxOe7ZWfnTMfqiQNEsVJ1yIPu5XHhQGqpc6gGOuKsKPtnZLmdLPcup7XOaq1ULnTtfSDl4DCN6ZZp7yllYlXxPdASeSTn2UU5X4qjXGeGeL+dlUnyxAHx5u3yEK976t30ShJ07e6pnEOlRFm1irjeUeKKrB7S25pbE6KVjlBImxXPdvPk9L6TcDEBlg9IrvfPGMPNzy8xWl7tqRDvrpPSBk8cgOqLlXTexk4vdwfyDrH5bIXa05SVsm9nHU45xCHLpbJkRrdBbt4Zi1Ku+tPbuOSeN5zaLpewvtxoT3Zk1QcpuhG9t2Hibh3mTeYcXG9x0Zom/Hbae7jqT2+lV07BuqbmDku5uTdwDl3ez+Erd74exZWv1qlT167Ly9GRh2/GOm4RZ7mvrBulIlaqkkbapd2fSlZq/3gv3rRtl7cXGQopKTpt13PPd/4829qLeOztlfj80cORy1FFwpGQi72rgFJNf/kmu1Ojt8uKiVOnoXdDAfN/ekbZbaUhSVM7J3Xo8pzDy40bsGhNE4DqbY66F92JOFRhLCYUHTgH3YZomjKb41ol5M8JGbnwrcMInPS9kv22OeykWSyq7n51Gb7/yDw8+MZyAKq+xLkfPRPhxOGqYiVDf7KJcNjTojVNeHyO3mciC5paKp/hUAfxNh6avbxD+qw0ujxxEFEtzkEXXdVpg6wi6/Ca4LAGuJ3s7np5qXP72tMpj63k3IqtfT0qEuDOYRm8sGid9ffnF64tpWcArlyl+51t2B6It7bsCv7bNnMGlkroCFJwvhLG62pmK/bz7QfnoLNgb1DAe+LQAdBZQTmd/Cy/PTlvFUZe/RS+epc5x+/yTTuNOSqm3P5a4jtP8GLbOGYvS0/444JqKqsrYnFWwvASJ2nG8A1DXKO4jq17zsXZNvHs4M3pDgIm7qtSZ5h2x2cu1+us8vrOTxo8cUigWqasOoW0y8nCVuWqP70NICnHlnHir6bjs45Jet5dHWTBsynsMvkBWDe/8mEaijLGDJ253F5zazte1s250PH6Mj1ineTzZUyiSHRmaEKTuMnIs5mybt6RXFemQ4jied5J5fVdgnMgohFENJ2I3iWiBUT07bB8ABE9T0Tvhf/7C9dcQ0SNRLSYiM4Qyo8hovnhb7dQuBsTUQMRPRSWv05EIyt/q+mocGy0CNpmKyw6MKFx3Xanep+/bSYA4Kt3zSq7TxOifA5VfG+UpDqGiW7ReK27DOu+mR/iKxpuzaQDMcEq2qnQBMmEMtb5xGWL1mxDa3sxEjkpYzGMM6M+GpdKllqmW8zqP7KngpXoZ5hoo8YExoVzaAPwXcbYWADHAriSiA4FcDWAaYyxMQCmhd8R/jYFwDgAkwHcSkTcnfI2AJcBGBP+TQ7LLwGwmTE2GsDNAG6owL1lRrUWoq5dN1PW2kG3MHU6mSzclt21So9PHTQ4tY44l3WS8r/N4GF9+0vmGEI2yNF1dSiXK3RjHNz7UL2zBTNLEL71wNsRJwrYvN7j60SdissetsPR2qhUk9c9Db+ZtqTsNmp976nEgTG2mjH2Vvi5CcC7AIYBOBfAPWG1ewCcF34+F8CDjLEWxthSAI0AJhLRUAB9GGMzWbA675Wu4W09AmASVUvGY0G1FNK6k6SLxUgtF4e263LFSiUkh+ndLZu1tZze1BTaQQyZLvu5LAxFbDqYlohYzuMUlQp+6rS9AqWsDb4O5SmRIwan9SMPq5InXFnn0FnFM//zyrKy26j1nWfSOYTinqMAvA5gH8bYaiAgIACGhNWGARBtt1aEZcPCz3J54hrGWBuArQAGavq/jIhmE9Hs9etLD+NsQkdSoyqLlctGNd/JircttCef7FsNxEEMHCeP5+5Xlxm7chEZXX6fXRkN2J+tk0I6yxwy81d1o09/Dx6enUwSVZpyXH+VLI5Ji9u0N6MziJUAAETUC8CjAP6NMWY+WunXFrOU265JFjB2O2NsAmNswuDB6eKGrKgas6INn+EiFqjd4tCNL8vsVEqmXu4UtBmCzCXt9MMyhxs0cg7CZ246WiqcQquU0G6s86ns/OvqpPVhiv0ncw6dlHGoCGp9607EgYjqEBCG+xljfwmL14aiIoT/ueH3CgAjhMuHA1gVlg/XlCeuIaICgL4AKp8cIAXl0oYfPTYfI69+Sm1XU3dvX/RacQT/LUs7ZY7DpHNImmvyk3p6eyajhawHCztXYP7toH16ZepH3378WTfuzPeieUppc2k6HMnlriaweyNqfesu1koE4C4A7zLGbhJ+egLAxeHniwE8LpRPCS2QRiFQPM8KRU9NRHRs2OZF0jW8rfMBvMCqeGw+7aaXcO/MZUp5uXzD/a9/5FzXyFYL5dWYAddp1asMKsNZVeK+0iKHcgwf0F1/veZpuwzLNAdZZ8bWl41ziBwJs5z+le+iQlr6jQVJfbTtmCyMNOVpnLE4/s8dPSz6LOejrobOYcnapoq3WQ3UIp+LCBfO4XgAXwVwKhHNCf/OAnA9gNOI6D0Ap4XfwRhbAOBhAAsBPAvgSsYYN1W4AsCdCJTU7wN4Jiy/C8BAImoE8B2Elk/VQEtbO95btx3/+fgC5beqKaS11kp6iO9GqYtj685W/OWtFdrfRl3zNL778NzUNsqlzRZ9dKb7MnvS6n+QfUp61qcrtHlL4kZ0+s0v4fjrX1DqOvtXlAHb/MRzmB18iKIxREV8BjXPIm18IkdQn4+3IZnQyMSiEjj95hlobi0/VlO1UWvOIfXNYYy9DPM7MclwzVQAUzXlswEcpilvBnBB2lgqAZ50pL6g0kXbi3L+McPxyJv6DRcA1mQMg2w6EbUVRSuaTE1G+Nrds/DWR1uUcv4SP2ogHMm6alkmnYNFZ5Ep8J6wzZxy8GBMX6waIojtXXrCqESqy0wmn0LVJWsN/iEVIgL2GFblXZ96rYVzUPspraO0E7+459cJxEEWI1VLIe0Sv6qro8t5SMcpQdXfbGKTw4f11bbF2zv2l9OUPqJ2dSIMw9oUrzUt37QX9p2VenuB3RkygG2sYuYu2/BHXv0Upj6l5twGgMs/dWD02fSsckpUVn0/SSeuoJLLRmSyVqqkJ+9GJ4V26ZubeJvyPLpwLWp7JegchEGIgSlVsZK9HSBw5Js49e9RdGUX/FyT190jiS5HHHSnRA7b650j4OSDB+PI4TGROOM3M3Dwj59R6sppGLOIlRLEIYOMN9G/wYliR0t5rLR4wkuDbYhpp/k7/rE0rituZA79ykEO5Z6+8IeZ+P30xmSdsJKLfNt4fshMG8x9nfv7l42/6UKxpPYk3df25jgyqXw/8tpyIUGlMBfiXIvPTH4GLs/kjhlLsa6pBdNTgiGK6AyRUmstVupyxIFDN+82nQNRcDYUX57Gddu17KlcpiUOhief5BzS6+hgWlQ8ZITsReyKbnXuojhxDH9bsAazl20qKYe0WNXELYjzlMY5zFq2CTf+bbE+GKLDuCpGGyyw5c0uhATaSfQkfee3vK7JLAIVs9Up7Rn1P2pZFrGSSPDktb3DIcQ2vzyraHdPR2dQSO9ViNasjnOwvOH5HCFH+kQpMmTb+ixipTYHzsE1XpKpz1JFIC56Gp0C8fL73sT5f5hZkilrVqjh0dN74x7Nug3t+YVr8eU7X4u9qDtA52BDXbgTup3o9bVEAiDfTotEHEohQi7XiUQgl8hHnazn4m3OD3W/fr78kBV7El7S6Nc6El0uExzf3LVOXpYXP0fBydUl7IV88tPnkE7W+cIfZuLsI4bizMP2TW3/vtc+TB+EBuVuyt3q8ql15q/ciiNH9LMSXzV+jpt21oVLkSVfRp2D8Hnau2sxekgvLUfGQ28zFvRv4i6zEo1SnwWXz2dRWvP/N/5tMeau2IptCbFScuC6gIRp0OocUscm6BwSxMEuVmKMKWPO7aVH3FnLOtzVK4G9dFrToVvQaWKlHLnJQNscKIhO3PGTJxZInIPa15qtzXhglrs/hYjoVF/i6Vcn7l6+KXmyizYv2/ZQqljJob5JrFQsMvzymXe113CiZ5PW8Z9cxtDdgYiWCq730c3v6CFJBzl5rba2Mzw1b3WiTF7yH6xP+jg4iTYcxEoNEtcpznXOonOQXwH9M+rIwDcdh2qZ1jv3X9Pea4DoNJXxuhwRiNxObC4WFi8s1ivP0hTS2wzJe7Kg1CXncl8NhfSNUY0SaqmbsKyJP5uio35sUM/Ed77ZPD4PPtcyAAAgAElEQVR3Jf5biMQqthXrUswDiQi10QkuLh87tHf0efI4PSco3te3Jo3R1mlpVQ8Z/JStm7Me9cm5t0hQI6SJGN3EShrOQVZsS99Fk1WbzkFuW3+oSx9jZ0Stb6vrEQf+Xyv2MD+OfA7OOgcXxx35BMfRlmLKWs6CKdf6weX6OovYI/buLXUEJoV0jH37Jj2i+Sb0/x5KOv7NX7k1+twjdJSzKfrTmC5x6SSUrQ7KfxMB+fdH5ylleYsMRRHJVMCWnxk+J+polQ7Jr/LYmoUQ3iJxSCUqmnuq9QlbhyG9G3DhxBHpFS2o9W11PeJg2ZkKliNIjgKFtMv7pmSzyvCU0zgHXfe7drfj+mcWpXp9lmv9kMWhTAcuQ1U2AGt76ToHEXIV04b/l7fiZPWDezcASBMr2RXSYgY9sRnTmhLva0DPevzkM4eaOxdwxckHKn0AQUY32b/FyfqqAhuQS/gM2bmtSbBCEjd3lXNIQndPeyLn0NpeVMLHZ0WtiV4XVEibYVtkFIqVXE5jcpUsj9jFlFXGH156H3946X0M6FmHVVvM5nwdwTnwTcG64Wfp06WOhYC4eMK6+DmkWXotWLUNT81bjbOPGJqYKFnWHrUnfXfZCH527jj07V4XjifZwkV/VDP4lc83AOubWrRK4GQ/6QppxpLKZNHXQjQiSNc5qH2lHb5a2tpx8I+ftdapNFrbWSa/IC0859CxKDVsQY6CU6Bus/mZ5G1ZTrCwJkGn4Oqoxz2ftze3WXMRlLtZuIRSiMV2NhGNvAFUYhvTo93iMyD377I2bPvQax9sDOoKZS7K6cAKKrVaav8yXOa1xSHGUCl+NWle06JVlFUh7aBzsM1Ja3sRK8tMvlQKWtuLqCuUt7vXOn92lyMOti1ymsXDMk+EhkJea+r3x1eWJntQxEruozv/DzPjdhyvcd1b0zaLddvsTkRu5pMOnEMGWiDWzRreon+POjfOIfxv5Ry4WEkoGzmwR6IO3+DFZkzmv0pXDouEsfT4VMePjnNk8Vu3zfegUKRWCs4Yt0/Yvp0QxOPR69PyFj8HhXPQGAKKvi2yaPVr//MGTv31S+pFVUZrezERULAU1Fpc1uWIg7jYdLkXTCAi1BdyVg9SjnLESgk4npKizGGGZng4j7RtcuIvpll/d9G3cNm7nXuyy5VNv8lhSaI6QiVxfvK5XCazYhvx4XVEU9nPHLlfok4uclCL29E5DsogZOEc7BX37dMdZx8+FEDlwl2biMwJowcF5ZprdKKmRMRh4XPBYq2ktmF/J4649rlEGy83brC2Vw20FxmKLEhZe9T+/RQrMlfUWufQ9YhDidflKJAfyx6kOrgopI85oH9qO1nHatoLvnTHa9bfndt3GNHnb5vJK5vbyTCOnz8Zi+xMiXtMyOfcuA3+vI47UMlMq9RJJglK1uEvs1huOkzIc5lVhGB6FnV5ihTsLKUukE2PJCNnMat14Rw+MWoA7vn6RIzZp7e2TtBO8vuLS1TuXtxEd7cXtblaOhL8EFNXIIwZ0gt9utWV1I63VupgpL0MprSSuSycg1RFdwJwkQf/7oXG1DpBY/yfvs03lm1OViwRWTZ1F4eyrO26cAEiCrmcE0HhNWw+GrxOncUCJa/ZLEdJfhdKgwDgqHNgjKWKlRKB+VzEgOlVIsjLOPa5cBMr8bI1W5sxd/kWbN3Vik8dNBifGDUAFx13AAD13ZHb+Z9XlmnGlRzY5ipGFHYBJw71+Zyz+bsOtTbC6nLEIQ07DQq6XC7YPNqEMN0mOOWHFj9n2HVnLDGzyWnNlG+t5N5AuadVHWwB6TjEfaJXQ8GJc4j1JOljtoa05nUBfHrsPnjiquPxxY+72bq7ihDiPND63+sSiXOYtW7wmztnJYP7XOimWC9WCsp4ePtFa4KMbESEc8fvp+1LbkW3BmyEtZre6ibwMRZyoYVjieu9V7faGpN2OeKQRsVN1i05oig0RHbrjTRiYv05wpK1TfjZk+Y49GntlCuBznK9bSwuxHPF5p1KmUnnIEIUzwzp01BSKBN9HXWjVURDFJ+kiYAjhvcDESnmrIwx/HXeqvg6kLMIIXIkNPyez8VtOR1SnMRKpr5gHEvW1KF87mR/CPmSrBGFe9dgg43FSjkAVNJhqKm5Fb94elFlB5YRXY84OG6gsvMSD58BpL90aQtcLns9NIFMg+hoxXGAYDGTdgrkP7voTWzXu9W1nMKV72rdVzSKRBNRtolYnDgHqBu/qQ+mKYv7iz+Lq+fo/ZP6pVlLN+GBWXE+AV1AP12oagaRc9APtpCnVAKSbDMDZ2XUsejESjrOwdxH3tCWPD6dU+GtL74vXWOvX21ExCGfC7ma7NSBc1W1RCpxIKI/EtE6InpHKLuWiFZKOaX5b9cQUSMRLSaiM4TyY4hofvjbLRQeFYiogYgeCstfJ6KRlb3FJNJFL1zEkERdKD90aSONNZbLnl2wxt5gNDZDefhfTI+pr1ce75DF+qVcUYYOw/v3SK1DBPzflcfjqW+dEMp70+HEOWjqypeJ60OXaY5j525VdCmrMrgRQdZx1uVyqaKnrG2m0Vd5XSxe0xS1+7Nzx+HHZ48N+zI3xOdOZg7lS0YNSgYXXN9kz/6WJTpBpcDFSvX5XMlipVKSOlUaLpzD3QAma8pvZoyND/+eBgAiOhTAFADjwmtuJSIu9LsNwGUAxoR/vM1LAGxmjI0GcDOAG0q8Fye4bJDfuHe2cuI8cEjPOEE7Y9bQxss3JUUii1Zr0nYKq9606Q7oWS9d4qb4M8FUd2jfbm7Xu3dV9rhk652/XnWC0cIrkROZgPEj+mHcfn3D5EwuMvX0Mcae33Fl40kaLDF+lcNI3htBvd+5K7bCBtOQu9fnM+b6dqmjr8Q3XlEP8NyCNTjjNzPw+JwgPEmP+kKkB7HNMyeOaQerR99akXj3rvmLGn+qFuarIrjRSl0+h0Iuh1YDp/7Oyq3YuVufzKgWHI+MVOLAGJsBwDWw+LkAHmSMtTDGlgJoBDCRiIYC6MMYm8mClXYvgPOEa+4JPz8CYBJVkdy7iJWeX7hWKe9RXxBefuCXFnngtx+cE33e0dKGVZoMVeIwdC9N/x51GCOFYDYN3fUk7mThUsL1uvG4hKLIggOHGKx+LHDNv8HnzxTpNaijlrUXi3j48uOi7ybOQYbuxXdZ8aJYybQYutflM5lAlqOX4LchWoS9FyaiejcUi+QESyxbX3zulPA0mmtEkaMu1/nbH20x9tMR4Bt+j/o8ejUUsGN3m/KObm9pwzm/exnfeuBtbRu19nEAytM5XEVE80KxEz/SDQMgJmddEZYNCz/L5YlrGGNtALYCMBucVxniM7z+c4dHnxsKuYTOYclaN5kg99D9+vGjLH2qL0BdPuco7mCY9q5b7lwT1+S6ELMQIVtNJ/FUie+GeALPkavcPcCfXjfnydDpJVraipg4akD0XVTQJsVKSSicQxjU0QUkcCc6dK/PZxKl6Fr5xonJtWp6XpwItgsUOFojjI83HrMtLhmfE5cc1j3rYyWzSxrRjsauUGzYvT6Png0FFBmwS7KC5NxEbGaeRGcRK+lwG4ADAYwHsBrAr8Ny3R0xS7ntGgVEdBkRzSai2evXl5ZCL51ziCsM7BWHFijk4heYFeEcVIu/LMP7d5fK48+6021dPpcahIzjgw079D8oY9GXV5xzgFsQOxvkEdmcxMT2xH0xR+RE0Bhjzsp8EXK+hYS1EszUQTffzsTBMh5Ab7ppz7SnFskJkyIPaalye5GhkCO0Cjt6UXqrCWTc+BN9hvVlYw6diKhnQ0wcxAB+OtTiAM51St3r8ujZEDyPHS16MbSJYMprxMVSr9IoiTgwxtYyxtoZY0UAdwCYGP60AoBo2D0cwKqwfLimPHENERUA9IVBjMUYu50xNoExNmHw4MGlDD1d5yD8LL5oPBMcEGx87sQh+C/vB+I4dBtpfUElDrb2y4ELbSByV2gXGbMrpOXvLsSihJc8UAYGjR+yb29jPcbczZPFjVbWO0W5LAAkaUOybZ3OYZh0eDAhbR7yuWy+1rpnKufhNj33gDjkEo6jfHp4vhIxqGApYiWdiKggmLNuDzkH0xquBXFoDtdFt7q84BiZvK/W8ERomhN52GN+9ExlB+mAkohDqEPg+CwAbsn0BIApoQXSKASK51mMsdUAmojo2FCfcBGAx4VrLg4/nw/gBVbFMJ1pLYtmnt3qcgkb9Vx0AmLO9tb84ZvSVwJ6NqmQU3NHlJ2PoQzOoS6XU8ZjOvUUU07ilX68YmvinQScg9v1adV0YiV5Q4/i9zO7VEyO888QKNG/c9pB9jEInYvjvefVZdHnfI4yieR0j1BeD6Y5bGcMu1rbccc/lgrjSlYOQt3H740JOYc6uvHsdowb1pHgB428IG2Q55nXMZ1JdOVrNLrLaiLVQ4SIHgBwMoBBRLQCwE8AnExE4xE8k2UALgcAxtgCInoYwEIAbQCuZIzx49UVCCyfugN4JvwDgLsA3EdEjQg4himVuLFScfl9b0afu9XlMe27n4ryJJPwoPv1cIuXEmUQk09jCbGSSefgJlZyhTFGjsPxKp9TN1pZBMDBGLDUIupSOAfNq12qTYJ4mcg52OYujdMJ6gT/ebU//vMEnDgmyb3yE22RscScyl698q1xDmTcfn3sg4CaTW99Uwt+8sSC6PcckSKCs92a7r7l9WBaNyK3VSwy5QAUjFcfc0oGJ0gmDu4f/34KTvzVdKWdtLY7Muz1tuZWMBbPV0Acgt/kOeRKfNM7pJvzC+94DdO/d3LlBpyCVOLAGLtQU3yXpf5UAFM15bMBHKYpbwZwQdo4KoW0/XWhYHbarS6H4f17RPb1cUhmZo3DM0QIg8w3eNvhXLcQ6gq5iqR5FGFqzYU4FHKq/N70IhcZw2UCkXUeiAWliZVigpYWGiPtxCo7gn3qoCHKCZtzBEwab5rTYb/u9eF4rdUSdfj9yKIt0UPaBbp5kbli0zLs3yM2td7Z2o5eDQVlk3YVK/G5NIVYHzFAcPYUxiyLwGqJI659DgAwrF/AUeYFrkl+VyLOwfIOyejovBRd0EPafWeSCQA/hRSZvZ2zDo+lbvzZyxuweLVufdTnyUkhnc2fwMA5OKyCfF51KDM5IKWfwivLESXvK2mt5BRfyKkPXpdJvcSIA9Elf+cbuC7fAxBYtbhCVkjLgQXzOXVspmx0urEAqrHFLo0t/q8+f0SUzwEAPli/PWwv2WC94DxqO+vwOEJpCmZ5zDpuRUQtaAc3iRY93+V55kTQNCVa0+nqSdu16HrEIUNd+aXNCac219hBUahneRxCnSfmroKMgkbGr0O2YHh6uJy+dJzDtx/U22inncKffif2CGeMaSPdZrFWStSTrZV4P5ZrmINYSXciVjvnfSXTanKrJp7bwcTFON2jVEWOHSWGeeEQrXtk6EYic0Q/DTMdiqfcM8btm7jHZRt3attrEPwubOuiV30BRHImxPS1naYvy0IbXly8Dh9tVGN6iWhubXd+56xipfC5EYL7VHJn64hDhSUJaeh6xCHD/MpZvMQTkG2hi6wx/5TVqaWuoNE5aF7lLMvFNOS00xegJ1ZvGZyN0tbwjCWxGfJv/v4ejvr586n926bP1F1S52AXK6Up+13iL/EmZM7hsGGBLqFXQ130e7ngTcj6jMBaKV2h/OmxQ8Lf1B/ltbpuW8Ah7hb7kp5HnSG3Q0NBDDtjvvFcjtCjLo8dQmgRU3WdzsEVZx2+r/G3f/6fN3DSjdONv7+3tgmH/MezeGr+aqe+kgppvc6BCLjp+SU48IdPJ7LYVSpRUznocsQhy3baTWLHoxNQ0c45iJFd+WlLXMMTDuifukHoxEpZIl3qne70dXtZTpYc+ZxbnCIgGzfz59nLteXyO++6BVDic+whbd/T7c8TEMVKfHxxT7d++eionaiuMJCbvjAeBwzsgT6h6MTYVcpNHrJvn3jjDwck27/ruEDd8/jl544Qm0lAdtjiTYomq/J5gos8Nu1MeiyLxCFtjkUdUVBff4FIyAf1qtfWEdsUMW6/vvZBWPBSeKh5Y6lbwAgi0cIx+RvnAgiE+0PnS+7Qt76pBVf/ZX7J46wUuhxxyEKQC5LsVZQfvv2R3rMRiG2Yxf7EE46LjLmQyznmItCX60xtTXVdYivlNWIlEzqY+01A3AxyJGyMNjFgMZ2gRRuVpt74Ef2Un8TTe8+GAo4c3s+JizHh5R+cghPGDBIU0gFkziGnUUgrop5Czmh/DwAPvpH0FOdrtzVBHJKd8LUqe5k3FPJOCmkgFLEk/H/09cRmBoWOqsakShIuOWEUhvbtpoSmcQE3m+3mqCPKC75RijGHZS5+9uRCvKuLx9bB6HrEwbHex0eqQd7E4GC2kLq7dqvsoaj0dckOVcirpqNZwyDLMFWVxWeu4zGhEiyxwjlYxAdid2I4BeeorEhfF6LuQh2bVIcxpY4YysMoBrP0zy3mZIW0zDkEpqzS2JV1JOrPVMg5TfjaFcVKruauDXU5IXyGtkoMgiPnIH5m+NignvjrN08wNZlAXT6Ho/fvX5JfBJcCuAbFs/k58HdZd4+mbJQdja5HHBxXhc3+O23ze0ZQuEbEgQiXnfQxAG5xk/I5nVhJd5G+Id1maupTdsrSIUu6w7R7O+mgdO92RQHs1DPQX4hkm9A5WPtyV0jL+oRgbEl5uyRVAhDMn4vlVBrEEB2AmqO6rb2oUC810mkc3kMvqkx+13EO8vIy7WeBWEk/DhnynJljOjHhMzCwV71ZNCo1yp0ESznAcF2iq/lszqKQ5l91DqMmc96ORhckDo4bnKaMDKcAHRas2qrUvXryIVh83WTUF/Qn2unfOxlv/vjTeObbJ4abSfJ33cnLNBbd4caWICYNOelU9+aHZrlr2hyLfiCVQdwftzEHwqis0aZuHhNj6WPmvy/duEOZc9n3gDGNAp3E5+dO0JU60vd3ViVDe+9uL6qcg/S9tT2mcLr7lkUeUWjuNrNYyWSvH4iV3A5VJMXCMoqVxH6ZFMdKwgfrVWfMHJH9tGAAv8e85jBlUuybvMP5vRUZsHlnYKFl8omoFboccSgHMUuf/vBi2/KYc8jlCA2FPAh6+f2oQT0xsFcDxg7tkzj1cujklJUwZXVhkwlJTkb30nFkTQ6jg5PpaApEgmbrcXd7Ef/13GL7eML/PGZQYmy8DgPWbWvGmm1qmANxMxXv7YFvHGvtN208v3o2Oe7dbUVlrnTzbZtPeXPKETBn+RbsEPwdFM7BIlZyCbzH25Q3fh3EYln5b8NnjwoCQZOlbRv4PerCqmnDkBDB5Odg61+cfxcuu1rocsTBdUnYTPzEXz4mKcK+PWlM+Cn5QiROWg6LWRcXSB7TxccdYHzhxNSU/OV8Z6U+gYwTcZA4B5vXb7q3cWp3ShuuOofkNfGcMQYcOVxvqfLY2yvxv6+Zw3Xb+hDHxgDcHmbje/PDpMFCwiFPKBej9brscRGXYhhPwDmkaKSRJGgydKFhzvv9K4nUpiaFtAw51D2gvjPymKJhG5cYS3xyPTeMHBj0qwvl7nLI4qIznem3Lld5Lmf2DnclDhlTZlcUXY84OFIHXTXdg37u/52UqMND9PIWYp1DevtyXwrnIK0/mQ0XsU+f2AKJE4frnnpXW1fHJssgSblbDnEQfzfVrMvnEor9UpCQ80Mf+wdw2xhs9xRv2CxKdiOnAhU5L7EpFx8TuR3APG9HDOvnxDmIWetk/NcFRya+Hz5MJaqu+oH6vOrnYJttJ4W0WKwT4RkQOaSSm5m4iHdWbsUfXgpyVRMITc2tuPy+2diwPfAB0RHHHNkU0ua+xLZ0eeM7Cl2POJShVNVZXcgnKJmN5CdI8eRLwUCsyGkWsBzbXj7NixjaLyYOu9uK1g3QRedASG6itjSpaXuti0i1kKfIdLBUiHqbwEJHf58u1lpWzkGoE4tQJNFMTlRqiyfDjEdDi64AAPYf2CPVWgkQfHY0vw1U0tPGn+vzOfz3V49RiBrf0ORrSbNBmrgMkoweXK2VnJNVRf1kE/MAwI1/i8V3Rcbw0BvL8bcFa/H76Y0A9ErkvOCtnolzEH4r9x0oB12OOOg25bk/OR03feHItGpazkFel5H9ePj9R4+9k7g2uCbdxFJUpnI8MEu1P2fQhxAY1Cup9JUX2TlHxPGfZJ+Izxy5n2Y8GlmvhDiRvPqbCL6xbWtuxWpDGOIic8szAQDXP6NP2SqLlUz+HPWa2EPcsS0as+WJiRZEcXRRtY5OrCQ+O9Me90/C87Dtg1/+xP7aOrqRyxZWItQ4YHGl4QO644xxqpcx3/B1B43ovYnCVBuIgzQeWY/xn+ccqoy5mIFz4Gsh0Pnpx2/C+2HsKF6XpEOgbP4LcGulJNckjts0PnEscuysjkSXIw66qe7bvU5xeNNBp1ySZeH8waoKVYlzSO3LLu44b/x+kWLNxbpBFgP9dspR0WcXsZLoM7B8004lJPeRw/tiYOitKr78nxDSaHLwny+7d7axvyyK9mmL9GlS5U3AFEm3XvPsDxws5e924RwQcwLy+As50q4Nl1wat1x4VGodwMyR2RTSOqJnygQH6OfK1EfcV5JzMFk2yQppucmDw4RNfG6bW9vx5oebsW2XW6pQ3p7u3UpbbiLBvOn5JVi+KRmDSfYql69TLA81Heo4q1qatXY94mCYa+UV1clpBSc4E0w5ftXFGJ6iwoe/n3SqFUUiH23ciU1SIvWxQ/sgp8mxIOKNH306UpDLKS3zOYpsw10V0kUWBMk78VfT8cibK5Q6/EUQY8+s365GbuXzt2Cl2Qv0508uLMlRKTkeO5dnK99fCBENpCmk4zomsVI+R0IkTr1YySXwnijCUg0MwpOxfPK3iJWaNXodRT8mXG/KgGiTfsgct3m/s4fPEIkwEMTlAoD5BkMLGVFEXdKEs0ihDvKcPPrWisQYlxnyl8hcUzQWTX86Ql2L9KAc6UF19jIYI2K6sOIp9tq//9LRkZmbbLwghjkQT0h8UX4pFAlwiBvbSTdOR0/JZb/I0k3yBvduwH6h7kEnuzSJAnRbFI97c99rHxr74/MjymfXbzMTBxs2bN9dlrMYEJyAXfwcdJBDnPB1M2ZIL6yTQpWLSmJTLB2RcxAXV9bsfuLh45zfvZwco6Ep3b3zMbdqdmo1D7VAzAwHCZuDnyKCMYmVJGWcMjTJUosriE0wcQeFfE6xLkpbkzI31RSGFueX6eYRMPtGab3FNU2cO34Ybpn2nnVs1YLnHEK4RLLkNUwnn7OPGBrVkl/6hHep0D6XKcqiHSJKsJc7pBNekTEn3QVvV8fKRxmrHJ2viowlIkfq6sgQFeMczvt0mcSBED8HBrM4z2U8fPr6dK9TrXcEJTGnsyrnkIs4hzkr4mi2OQedg6Yr7Zh5mdtBJ/ivE0n262FWSJu4TGvEVUmEZaor6xzkNUuGd8sEuRv+tVdDIco9bepLhknpzcdizrLIx5L8XVefSf8BoIdwSPlc6KfRUeh6xMG5nkasFC6QJWvNcZVMdugH7dNbqBNbZfATjPzS6dJyiujbvS6hcDWOJ/yvl3FywqQu/CMknwAyLHIOBv0L9L+XfAL3fn2i1K95vIMF7+n7Z5m5FBckPaRLc6TjEGPhyO2I33ORQjp5k4Ucoa29iLb2YuS49s1TR6NPN7d0s3JfuimMExHJOgMd5xDAReEpXi+foH87ZXxYJzk+3ZiLLLCcM5lBy0YPsrw9atvhJWaMGS2EetYX0NxaTMQwShPtmySvvAsTcTFxknqFtH0MpabOLRWpxIGI/khE64joHaFsABE9T0Tvhf/7C79dQ0SNRLSYiM4Qyo8hovnhb7dQeKdE1EBED4XlrxPRyMreohtc5p1vfo+9tTJRPuP7p+C1ayYF7UjXHLV/P5wwehAOFfIDi3VMoh2dnwPHD886BBdO3F8JaaEdM895L9T7bpjI/v5Lj8UXJgxXrHWIgAcvOxZv/OjT+Nm54/Dd0w6KnPJs/eleoCF9uilenjaCdt15cSZZ2fs3K0TFoxhPiKN/6OzllFQp/B9YUenFcIzFqTOv/adxiTqFfECoRPHip6R5KffVN3IOFlFPe2o0PEmBLjX+mSP2M/Yh9BbWYTjox89EIhm1VtKUVZa3p9GGjw2OnesY0zzX8Dv3RRK58TSjDjPnECAtjI36Lus4BzMRB7KLRsuFC+dwN4DJUtnVAKYxxsYAmBZ+BxEdCmAKgHHhNbcSEeeLbgNwGYAx4R9v8xIAmxljowHcDOCGUm/GBTaWNllPrcMfdJPEku4/sAf2DRXKskVTS2tRa0cfiZUMkR51fg4c5x8zIkrqkh7MTNWTnBkmPJk4agB+df6R2oXfo76Awb0bcNFxI/HNSWMisZKtN9eTjW3Mqry7dIhe5jrO4f5Ljw1/S3/pYpNYHecQizsaCjkUcoQvTBiRqMOfryiWk3ttcLh3mwkqi+roy0Xw5WaSlb92zSQ8+28nAkhufIp+yrj5qX1l3drkoIKy7sKGdgvnwA9DIvFJWwNG4hBeZvPdEPuOx2Juy3SDHZ0AKJU4MMZmAJCjrJ0L4J7w8z0AzhPKH2SMtTDGlgJoBDCRiIYC6MMYm8mCp3CvdA1v6xEAk8h1lykBpul1kveGlba3mL0WxZelWGTYuKMFDXXSNAsnfpvOIe00snprs1O8GiB532pf9jZ4nVI4Bx1sY95PCJynwytXn+rWCWILKyDUORBw+ac+Fv0+oKdqegsAIwcGlko/mHwIxg4NOD6RyKiOj/FvJvEav+e/vB1znbIo4uj9+zndU3A/mpOngXPQ2eBHnIPw21+vOgHPhx7/+/btFoWbEDdOhWuKCGNyDNo6aQcZiRNWOAfLvQNJoqg33w1qcIsrsX0TkeQw/xofGkT8+V+OA6A3fzeNL23r72iXh1J1DvswxlYDQPh/SFg+DICY2igeR8UAACAASURBVGtFWDYs/CyXJ65hjLUB2ApgYInjSodxgtMV0nwTkM1CE60Im/FNzy/B2m0tSn1RvGHSOdj8HPj1plhJyXZUGbgsGpDFLbo9nutJbLyDzVP1h2cdEn3mLcgc2G+njMdoTRKWK085MPo8LIV4iJB1DgDhmjPHRr9HDmvCLU35+AhM++7JAIArTj4QPzprbGLMRcYUIsjnb8uuVq1OAkDkF/LzJxdGZWp0V9KGqtBBr5AOCkU9xtadrXgozLZ3/6WfUK4R5fqHD++LMYJuTFdHp58iWZMsQSSeNhDieW5rL2LD9sB8WwyY59IOryNvwFedOhpA/K6J+pZScyhwqZxsgcXv2RxbSTfm+CDDIR6WdluiElQDlVZI63YHZim3XaM2TnQZEc0motnr16/XVUlFOaassXLJvDpjtp/h8bnBKVG28BEVyWadg5lz4PfwzUmjjeMQ+wKSLLrs8+ZqJWMaz5ghvXDtP42ztnPk8PhUbCJ6suMZx9C+7gRBBH8x3129DRu2tyjj457hiZNxjrRey7F3ryo+4/N524vvazkLAFi5eRcAJEySS5Eh2+aYt3ZBKNL69Nh9sFgwnjh+9CClLbs4iBPPuM4LGodDvjZa24tWb2x5/Vw4MWm+LaYJ/d6f5+IboZMkr2dTxssoMjX1KyeaMuewcXsLvv3gHIdWVfDYR3JfXLeUluwnWaa2L0YyaLYcSquBUonD2lBUhPA/XzErAIjC1uEAVoXlwzXliWuIqACgL1QxFgCAMXY7Y2wCY2zC4MGlhbI1m7IqfSl1+J7BT1Jc/KBr6M6Xl2L5pmBD6C/HmxE+RwlENJyDMZhZKDM1easm2wnaFf0c5L5cpEEUHut0Q3r8quNx9P79lU3xe6cfpG0rq+y0VCEj35TO/O0/wu9JcGK/XvBbkIcmK0GZhnPgntcTRw5Asaj+DgA/PTdQUB8mcAY6e/+0e7U5yv3zJ0cCCJ7vIfv2Ro704aU5ckRWJ6vI3DVFnsETQY350TOJuZTbEd+pK085EL/83OFKXX7w+b85q6KyOLxLTMzTCGt7UdU5cPCDGH/3rn9mEeYs36Ktm4ZnFwSJvWQR4dh9A3GkTiezvaUNLy52O9yKBxGbGXk1UCpxeALAxeHniwE8LpRPCS2QRiFQPM8KRU9NRHRsqE+4SLqGt3U+gBdYFdXyRuLgtAMlRTS6E2LkJSzE/dclQY/Z55Bz0OgcAhm/OuAe9QVj/3/85wn41flHCOMJ/ouJWuS4S/K96+aCbwBajiqsL18mKuLFNovFIO+BK1w8h3VI04Fw8YKY8lWeb1kJakouc+Dgnhjcu0FrzQQEHu1D+3ZLzJ+OE3O9U926OHJEzJ1xztMWGoVgt9LRcQ7adsiuRxJNWeW2lXqadvhhSOQcdP2JIpii4SADxJwD56adrNUsc8AYUxzy+oaWcLrYSt//81w89nbS4jGok/zP8fDlx2HkwB5otkRCrgZcTFkfADATwMFEtIKILgFwPYDTiOg9AKeF38EYWwDgYQALATwL4ErGGCd3VwC4E4GS+n0Az4TldwEYSESNAL6D0PKpWhDnfcrHR2D+tacDcHspZVd4XbhlXTvyCZ8r3jZsb8G1TywAoOMc9IqsRB3N0zv1kH0kS5mgHdGEUg6B4KqQLjL7gEzKSiD5cjAw/PSvCyFDp28A7AYANsjPR75Pnfxc4RykUy9j+nnnaV2LjBkXU44oId7T2sanPAyTH83PBRNgIBhjkTFraJS2IjMGPgTi20gz8xTFQSL+9m+BcjteF4L4zkAcGIDH5yQ3Tn4Yiq5g+jH97sKjopAws5dtMm7odRLnUO5ZdM7yLcac8jqxUuO67dq6JpH3xFEDcOh+fbBRE4qmmkgNn8EYu9Dw0yRD/akApmrKZwM4TFPeDOCCtHFUCuJC6NVQQO9QDuliGx7l0o04B7WO3glIFuMEp/Brn1iAmR9sBKBGRjUpsr52/Ehjuzrwdna3Z/NsVuuYc0fEfbmde4tMfREKOTKGzjbZxWeFqnjXEAdpXPK2xj3TZYimx7Z5aE04XqnzmRY9Xd1mA3z12AOkdsgpIKNNnBKJldKIA/T3wgPlOXMO4RqT5f/cS1g0Gdb1169HPfp2r8P2ljY8PX8NjhISXongXDpXQsstFYvm3B86yNMjRv91zZ8NCJyDhkgcOLgXnp6/JjSlLlHOmhFdMLZSDFvoApuHNF9Uus1FVybLfTnnkLAgkjkHg2flaYfuo4zHBr6QdreZF6dKvHR1grnTtcJTorq+T9pYP5ZrTe9V47rtRm4jGE+yUb4xP/WtEzBjyQbteOW++HOwmbICCMOd6K2Z4rZk4qDWkZME6frRjVPtKyBW5UT1dM1pTGQPEKfzzdDpQvgak8HjXEWE0cA5APE4CjmzD1CskA5+l+u1FRnqMxAHW1VdbCXTWrfNMh9ze5E55V+pBLpc+AwRlPicnHC9vXbwv2ih8LoHr5iOhl/FTUbVOfC+5BN2XM9l/cacg/nldWmHm4WKw3ntmklYct2ZMXFxXLNZ9ysTu924zhzGBFDv68E3ApPOcfv1xRUnH6g9gSnmpVF5vJHo5iufQxSywUS0AwVw3IFugzOJJ+TxpCEfcjLcA/rOiyY4XplEjty8h21hOPR5UPQiWd1715OLlQSRmhylmCNyKs2TEsmAg3PpnJDIfepSfop1Tj44aQxjO6SZYitx/EkwL+Z13tFEK47NrqumjlXQ9TiHBAXPRoGdTuq6Ms11TKoscw756MSRXAxiPZfxRNZKFmWWqykrWFIlPbh3g/N4EgppHedg2fZM70OaaV/a/Oh+VRXSYXn4nUfD1fXVHppPmtYVIfkcypF1pwWfy4ViJb5py1FmXSEHgNTWQQrnoCEOLnksAODUQ4ZEdcXot43r9XJ7Po66fA5Tn47T4oomxDx3CycCOs7BBlP2R1tdExcjcr62bqNQ8B2ok+6CnIOoFBOKFbGSClkRqecu9Io2qUQRT6g6h+C7nFQnK3Hg9+XC9svXJMcTipVYsixZx7bBCwpp5tZnVN9QbktVCjiYhWp+l73ZZe9eU1pKbh2kM3UV66SJldLgqgfI5YKNhJ80XTdjpR1yyXWQ5Ihk8Dm86+WlUZkuErAuyvAYYfMUjQMaDDa6V5wcOEweJDnzTf/eydFnrqTnY5b7TDPdTZvJi0OTYkAgDobXL7lfJPvt2z12ZuTzpeNqqoUuRxySm5ugc1DqmXUOUR1N+66WPwBLirXkvTL8/q0H3k6UFxLEIb2vSIlu8wB1Fislz6uKM5jjfq/PL2BGqZxDGmeo+/1qwYNaHFdkylo0mfoGdugPvrHcuAaIkuK9UkQERIT6Qi5VN5HPBZxM5GSZsljEzTPRHyhdb6HROYjhpXnP81bEHv1aww2oa0NrUQbz3J0e6uR6dYuFIkP7dsOQPrGSWI6tJPeZdr/y2MUN+7VrJuFfPnWgUpdzDrvbiliyNuZ6xNtT9F3CbznPOVQfZoW0g87BoX3XOgHnEJftaEm+7LqwF4DEOThQh0jnYBMrpbbCw3nYFaG2zTgpVtK1b+E6DLxDuizc+rOCJ795QuK0BiQtZABoA+8FfRFmLd0ExoC1mgRHvC3R3+TjI/XWNGloyKcTh0isZDG7FtG7m17CTOSW60DcIM8/Zjhu+uL4xO8ytJwMqQcunbc6Y7YoqPy9icczYWQyVa0cPsNJ55A8FiV+E9fhvnJGR8mg4fYZSX8IcW7elizHkjrJ8L46UOfQ9YiDMLfjDGG0AXvgrrit9DpB2yTVUeueIIU2kL2xObJyDrxvu0JaGp9ec4Iig9XWWh5Pvx7xRttf+KybW1MOXgBGuVLaa1KaICUJOS5QW1HvO+BKqLko44mrji85LEhDXQ47d9vNe4OItCwSkaRxDqZfidxO0u8Jp2G5LRMx1Y1BjnGU5O4jPi5BsLgoCYiJiagg/85pSU99OXyGonPQiJV4FSJ1nfP5ufXLRyvX8bp845fNssX7+/dH5iX2FHGK+PryYqUqgp8A7rhoAs4Yt29ULq9VHYWWF0WpYiV+re7hy9/lU5uN29HBhXOQNxqT/wZjDJt36q1EdOM5b3wsWhizT2888i/H4VMHDQ5etAw7d6lnpSy26oBJ1BE+h3AQ7UWmPfW6Eurdkall6a9efT6ncJoyIrGSo87BZmGVxjkQgPfWiaISWdzoRkxJo7sQkz8lOYeYyP5gchzUMTI3F8bcIOUrqYsU0kEYDh4fiUPHkfK+f/PF8epeEdbXzSAfzwOzPgq+y46Z0jIQOULScA5erFRF8D1fjoukZM/SPAQ36yC1TNEnhM4+LtY9ds7BaUcCYNc5LN2wM/o8rF93XHPWWKVOFIbDwWSRj1N+ESaMHICGQq5icenTQ0BnJA46v5VoQwr6am0vKh7mgJvCVxyObICQBfWFHHa1pnMO7UUhsGMKMTLqSZB+Ws0RYaAQP8xl2g1SJWWdfvHjsbe/aDlmCmETnbCFdmS/gLxQ58+zV+C1D5Kh3HScUpEBZx2+L84dP0xZJ3zMLsYoaUYcWwRCJSrtvUK6A8Afu0kBzKE7PdQVZHMlXQ9um0Ta4ZkvIqvOQWjg0hNG4YFvHKtpJ6hk29R3CZzD1M8eFuU5EMFYYIPPx6ML+icudFNvYgIeV5iIwG/+bk+83pIxUJmWcwjL0jkHl4OD8LKXaD0EBCffXcIJc0jvBqVOPhdwnXFgx6D8xDGDIjHM5SfFuS1MhNQWHTi+VsqRbtiw5XbT2vn68aMShJgTuN1tRaM3eryJxoOukwhjJCoE8OISNcqsbgMOnnsyxlNUP3y3XMRnplwgHNsE4qCTLHjOoYqINxr9aYJDd3qQ2dNSxUpcWehiHy2Lt0ymrMePHoTjDlTTYMh+Dg9ephKQs4SwwKZN67mFawEALzduAKA/+Yq3Y8xFQYFIIM1cUISJmGxvaUNTsznu0v/N0TtBmWATK/1+eiOAIHSKzkPVhUtJcg6lv3r5HCV0SL+ZMl5bR3SC4xvbfZd8IhLD6JS96qDTT6tESYsmlVNWMV0b+jspVpLNinlYilVbd0UESxfNGEjqDXTh8IGAeOqem07n0F5kUWgT+b29//UPteWALlaXnViYIq96hXQHQn6O8manE33IuZZdTTLlsq27WtHU3IYPNugdeQCBc2h3Iw4m+TqvwjcTnuNYxMGCTbjO/lwHhYuCO+dQZCwKdawbw+ePHp74zfY6HH7tc8bf5D3tmjMP0VcMYYu1ND9MrNRuUEi7SIlK4RyuO08JRwY5gJ+Y3IeDiJuyBt+1Y3YQUQY6B/sYCckNVW5J17TOw5koKQ6SD2M9GwrI5wi7drfjf1/7MLwHabyR4jYejxyzSwyGp1vvOqlBuxBvSd4bFode7boZlA1BFHNv6XuL8FxFyznet0vO70qhyxIHGXJoY90CkUUpt1x4lFLH5QQ5a2kg43xj2WZjnSiujS18Rk4sN73cwf/N4cuYdup1VeLqZNg8IiZgPu1zkZqIYf2645Erjou+HyDpg0o9LMkv8eHD7VnWdI9O9mgOiIOdMJqQ0Mk46hzkDRIInrvNwAAIw2cUGW57MeB4dMTIJKIUQQ6cQ44IrUIdF4W06QQscg66lLH5XMBdvPr+RmtfnMj8+OyxahTi8KspFIosNSgWGVZu2RUtXPl0P2JAsF51ap19JJGfTIzk8fPnWp/P4Q4h5Am/rsSEdSWhyxGHyCRNKpc3V52FhviiHjq0D44Yrub8dXnleziEMogXucQ5iBu500k0KH/mneCknmbS6HqirddsboM1sm8ZOp3DiAHdo+i4gDqHaaEiTBCJw7B+3fHJAwdZauufnXiSi+L2lGjKKlIfV2slrY+AxDno6BK3Vnp//Y7ouwzxPkzhS7him+PHZ6vGCu2MSWFpkr/rWv7WqWO0/Yk6h8G91PVUyFHi9GzSObQaFNZimSmIomxOyzM68vzfMjdwaGgSr5vDQj6HSYcMiczmf/38ksTvhGTucL7eHr3ikxHRAYTYSmUEUsyKrkccwJVHyQcpn+QmjR0CGUQUcQ+mk5/LCdJlA+aHnVbp1JbPi8QhvU25uJAi63bVk9799YlKme6Uq2tfPtGb/EA4SuUcxgyJRVWHDetjqanvF0ie0KPETFplfPp4XDiHqyXRl1b5nSPsdsjOJookTe1w2KyVxOelI7By9jdbXg8OnX4sMGU1ZyzkZbac1pzmcgKiey68rMgYFq5Wg9zJG/CWnUm9lsxI2RTSQbnZCCNHhH89OU73G3EO0rvkiUMHwIVzOHxYX/zq/CO11/PFZ9qMdQtEls+7+SfEpxvTOMmBc5Bf1LoS7d3FTeukgwYrsWvk8ZjAw3CIkDfKUuPVn39MUlfx43PiU65brgm1jhiwjotXtJxDRmslEwc3cmDP5DWGvsRk87oTa47i07Opv4ILcaBk4D0XLk5JQyu1fc2Zh2hzdxCSYiUTtyOOR/EPkjhum6XUq+9v1EbBlcVK8jsozwFfF2aLr4AQuYSq5/HC5IOWV0h3AERPRxGizmHs0N4K5Y7qhRe6Km6BOLdv3Ff6NaaFZlIimjYbuZk0zsFEZI4SUlDqREquIFJfNvlFUDkHtxdCNunsKehATGIfMYSFbsonHBD8fuKYQcZTna19EWINk1hJbka3zvKUtFYyiZXEE32a+a2JuBElT6suj8KU1RAAutflcbkQe0juK51zyCUIiDJfYUFru15CwPsB1KCWHGk6FvnwHnEOhvp8zT88e7nmt+RVXKwkW2p5hXQHIDZklcRKwiqziX34Qnc9qX98ZH+jtYQNpiEkdQ7pbcrlaYpQUzvi/ZoIpwt0Ogc1h0JyDK5nJbleXSL3hf6+fnT2oUK/KogIhw3rg/p8Ln5xUyy1THAh5ibTTBHbmlsTQQd1XU9fnDQVTdc56JGj7KIMmy2/vOmJkImDidtJEBCDgtcuVqKoPxH8ECfr+eS7NwXqs1p8MYYfPDpf+7sIUSEtwiukOwD8waonajfiwDeItz/Sp1eUN1+e/zY5hvRxygvtuI8NxMxrTk2cUBObjWHTV0x2UxShxgUu9FuOjb5O5yBb3pTqHybPa51BP2Pqy5yHIXi5bZyDOP0ypyg0FPVp4jRc5PVpCYEAKCEhdE6L+Xw68SRQZs7BpufSjUPsS4yxpZujfC6pjFfFSsF/V4W0iCkTA29smRjKxECmlbGHtNJVVO4qDorEStKBstPpHIhoGRHNJ6I5RDQ7LBtARM8T0Xvh//5C/WuIqJGIFhPRGUL5MWE7jUR0C1UxSappakU23yYy4uy8KZCdeNoZ1q87brzgCM0Y0h+w+oKREqhNfDFMY1Y4pBTOwYUjKoc48I1WxHlCeGdAfcmGaUwadbCFezbNjy1se1wnWDf8mWvFSkI73zx1tPI7EFvA2SyVVHGni7hKL3ri+I9zDtWKE8U6rhtbKToH8Z1I4xxEjkjLOeQpkceju2ET5RZHNoW0iG+eOjqO1pqyASucg4NC2popT5gSI+fQ2YhDiFMYY+MZY9wo92oA0xhjYwBMC7+DiA4FMAXAOACTAdxKRPzJ3gbgMgBjwr/JFRiXHgadQ0JEU+rRFckX/+ozD8GQ3t2UOi6HCJcQ4k7WStITTjdlNZRXiDjkcsl7eeKq4xVFsihTvu3LR+PSEz8GF8hTRInNz4E4WITGRQb8NgzX0VDQKFMd+pr9YeDXYiPQimmmk35KLRPXg+mZm4wb5PGIG5Jub/rhWUkLK8UIQrgJ29qRR2CyVhLNi2WRrRyTTK9zUMsuOm5k9O7KOgd5Q1Z0DlwhbTheENyTbbW0FUGkOuXyaevsCulzAdwTfr4HwHlC+YOMsRbG2FIAjQAmEtFQAH0YYzNZQJLvFa6pGmxu7FmUzTJcxFOlPF/diS0pw9Y/StVM1H5vxk1CaN4U+98FsrXS4cNUx7QjBf+R0w7dx9n3wqa4NoXZEO/L9nIzxvDU/NUADGIlocg0XH6vNgKtWvq4cA72dowirBIOQS6WWkrKW1EkaRNrprTD+xfFSkZzV26tZLEgkq/RhfsGdNZLye+tUV/ariBnAFQgXLdlZyt6NxSU584NZtIi5FYS5RIHBuA5InqTiC4Ly/ZhjK0GgPA/dxgYBkBU168Iy4aFn+VyBUR0GRHNJqLZ69evL3HA+skVn0U5QdGy5lsAgMtOUk/Gb36Y9J7+0sQDlDrimE3vXNZbcRG/jOhfWh4CPp5kHB7dZiOMJ8MN6E70HDw2lNpXOucgDyFNIW3a0E89JHgVbJuyIlYq8aDiIlJL4yKB5Lr6xomjEjlQojoOmzpH/55qqA8T9CE/cpi7XK/vi+oI3I5twxYPKYT4YCdzCvx58/fUyDlYCNGG7XG4EPlAJM7f6q27MLx/MkIAv6egr85DHI5njB0N4EwAVxLRSZa6upljlnK1kLHbGWMTGGMTBg8enH20MPs5lBJCQoeCg5Lv8k/FxODi4w7ADzUhsjcLsWd+MPkQnC0Ex9O1b+QcMm4uphdbLLdtwmnIUZwN7cgRqoc5kC6ieexfP6mUDevXHVeeopf128eTXkfmdtJMWU1tchHILksWN5fQEzLk7HWApEMziQoz6jMuPfFj2ueR5tcgQhfbK+5LblcdeCFH2BHOn0m5LVo92SyIEocrEjgHOYdKWPFfw2i2cpOckzFNpzxnt33laFxywqjITFq8rqWtmPCticYQOfd1EuLAGFsV/l8H4DEAEwGsDUVFCP9zm7oVAEYIlw8HsCosH64prwoiU1bpQYoPyOVEZYJLMLNzhSQ4E0epnqJAkpvo1aDfjF3CZ4hDEImSCWYrmvhzOZwVIVbsnnXYvto6ac0ftb+aXvP6zx+ufanS4Mo5mHKPi3VsvwNAt1AZ22KJiyRfe9C+vYx1AeDOiyYk8iNzFPLp6zBrDgpXzsrG7diInU2RHZUJ92U6DIhe1DZFe0I6RLHISw6fwQkN15eMGhQ4Kv7nOYeiW10OG8ODnK0veXz/cc6heOSKT4ZdxxV2txW1EY85se8UxIGIehJRb/4ZwOkA3gHwBICLw2oXA3g8/PwEgClE1EBEoxAonmeFoqcmIjo2tFK6SLim4og5B7Nsd83W5pLbT5hPOsyuSTkpxinqZzhtOSmkhfs6YEBPbR0RLmIl16BxOojZvswELXv7OpNhF7iIgwiUMA3V1XJxKNN5BcsQD8PLrj9ba9AgYsw+euKRyINgeF4uxMHFjFWVj1vatfwkb4q6+RIJhvG+KNZLmJ0NkxZEOYpD04icw+62In75zKJEf3yt1Rdy6Nu9Dh+E8at27dYTfXmvkd8x8etuQzIpnmNlxeadym/VQumaRWAfAI+FC6MA4E+MsWeJ6A0ADxPRJQA+AnABADDGFhDRwwAWAmgDcCVjjPPXVwC4G0B3AM+Ef1VBHFspWS6uZ104YVe4cA4iTBnBxBO8iRV3C5+h/2yCKT5SQobt0NChQ/WxjBLKSYNYoBTGxCWuk74vUTRnOvohiMrJv2qquXAgbif1yogBXTIGunDIRWmjdBnD++vNoehtkNdDN43Zq8v7RRSP1RwDDVizrVn4HouZRIIoHhQ5Z8HXWktbEX261WHttsAT3ZRCV55m+RmLX1vbi1px2eDeDSjkCNss+UsqjZKJA2PsAwBKACLG2EYAkwzXTAUwVVM+G4AauL4KMOkcXHIRuMDFK1eE6WQjvrj9euiVeC4btnhqcRmPyQ7daRMVcNbh6SKjSnIOJRMH4TIbgRVPzTpi7cI5OCmAM967cf04OJ656NZWCZujLkMgoG5+z76zBr/8nOrfkwZ5nGmclo0Ic18I4+FLs0HHpqzxw16xJT6p8/nia213WzGZc8GB6+bjM/3+zspt6DlKvy3nctSpFNKdDtHUli4ZsSKX8YTtIg7q7/BSGp3gHGTGIlw2Epdw01ecrJcHuxCZUjiHUkN6OMWnkhaLGEqZw8WU1WU9uNz7f3/1mNQ2xU3RNDfrtrmLT7956uhErCoRMjHfYVG4v/2hOYeJPE4dZynmezZN1eadrZi7IkjOZOMuROSIYic4Qdy0ROONfuKYwBjm6P37JYjDmQYdWpoFmjzC15dugg48R0dHocsRB34ElF/4rCc2F7g0aaoivvQ9DYpWccymDUC2yEiDKTBfItm5g87BhSsw9VXKsyjVgipx4rc4EqYF/8s53JcLUXW5908K4a5dHNxMc7N5p7uIwmZhJW92UzXZ6zhWWfR5IjH49iR9voesMIbK1lhY5XIUBhpUc3iIOGHMILz7s8n4xMcGJjhtczZGqa+c/XcTgoi0TlUrgi5HHCplrfTVY1W/AxnlEJy8w2YjNu/CgZRjZVQpay6XdkqZt1IslVzHI8cX0mGdlNNABz7/3N9BP55soieTTD0hVjIcHA7QcEAm9NGYy0bjkZrXWZO5QCQO//bpyhAHkximp2TAwA+LBUl0wz///TtJK32+3mxhMaK2UzgHedk98A01zzsQEB8flbUDYNM5mNhnEVecrA87LKKszTihuC19ExXr2Da433/paHzuKK3voTKe8oiMg56khOZ1tv4AcPLBdn8Yl/sK7Obtm8CT89Ktr3n7tttzsXAT59DEFYjiQZM+ZrJBDKLDRceZD0NK5N8S10f3+lzYXml6Jx1Mm+kgKbw7H3Ihl0u8J/zz/gZLv5cbN6SOQTH1lQpka8SRg/REm2f36yh0OeJgy23MMbSv3XwQcN2YXUdlhynkgMsYxCqvWBby2UcMxU1fHG/8vVLhRcQ5KYfoibhc42HO8ceLP269VuzJrLchY6BFjp/+U7o9BW/ddntZn6mJKzhTMAgwixzd59lkTg2U5rinw76hv4arWN2lGxNRl32H+JiDkOAC58Cz/5Xl26OKsESMHtILvYUDaY86/eE078VK1UUcslt+YPH3qwxRNbPC6SVxqOLimGZsvkIEStw45XDCWZA0vzUpBtwuegAAD7JJREFUv7O1eckJo4y/5XKEYf2644bPH679faCQp9hs8ZWOE8fY81MDblZwlXqmX5q4f/TZJWe5Db1TOGm5fdv4bAcvMY+4CTeen24F9d3TDoo+D9LkoQZUwws+5nyeJJ1DMdJH6OCyNsRL+3QraA97xwp6pG71Zr3fA7M+wj2vLnPotXyU4+fQKREn+zGjm4Ny02ReKsKFONji26chS7pRALjIlGfAAeJ6LtVsVB5PpXQOafPwytWnltWOOMx/NYgTXRwD43g/5rpZRYUmiPfSy0FMasL8a09PFSPuI3lo64Y34/un4KQbp1uDTrrM4QED0x05xY1cl6saULkpijiHXKRneO2DjVi9tdkaLNA1NS4AHDCwB176/in6MVP837Qn8OfwkycWYNSgnjjpoNJCCLmi6xEHQ8huETb773nXno48kZO3q8tLfMwBpSnvRAw0mLoCSSI43hDLyAUullEuEGMU2RyUsqAKhmZS+3EHJmWyix5mV2sbAPtJvlLEQYSLDs0El9O8fD86PxD+rOVcHiJcQsGLa8YURZfPjy0PSL3hACjmqJ5y+2sA1JwRWRFxJZbnxtdPt7q8keCIjphNzW1ljckFXY84hP9NCysNfRxeFg4X8Ui5ire5Pzm97M3GBS75AVywfrs9r3GAjJxDyaPJ3r7JcoyfLm1Ts6MlMAftYdmsq6HLsnF6Xzt+JMbuq/dmd4Voots49Ux9YqFw0DbRmkkHJUI8VX/h4yO0dXiVYZbowbyvrx0/Ej8WUsXmQ52DaLps0ze5iZWCWvZovMFvroRo9BB7zK1KoOsRB5OLdBVQ7sb8z58ciVXCaUEHk5UOBx9CuTTCxQPYBaITjznuTbY2q+GjIiLhm2EYHPf9sG1+pxwyBPkc4SufMFv+OKmpKih2+8lnxmVqSwfxNJ9mdm0TK/HnaNPfiLdiWvu8HdtBjhPMunwuefAJdQ5tGosl/YDMP3Hw9l0CErpIJACgpyEYZyXR5YgDRzX3k0G96rFh+27Fljorrv2n8l9cvjGUY2EEJDfsgb3MYqw0iC+aOe911s2v5OFkbt/oV+BA0Yb16473f3GWtQ7fSMrhzjoa2QIxmjdaTjhsllE7Lc54HHwObZwI50B0ZrhtRWZPziPCwcqA7wO2dcqH6uqv4yLuKxdd0Fop+F/NV297SyAP7GdJbHL31z6OP3zl6CqOIkCs6CrvjsXTVakRUIGkmaJpA3QJS3yd4IVbqojQFQmxkoHbicQmFTJDtzmd7WnI4vld7vxwPcIvPqu3PgNE4mAeF9ebyWuwkMuhrZ2htc1toJzDuOXCo4x1+CnftvHnMoqV0iQGlUCX4xziqKzV21D+64Ij8btpjVYTwJMPNnvJVhIUyTsr0065EE9kJp2DC3EQTSLlcASVhkhYTadRa/rLDOBzckoHrY9KwIVz4KfnczRJqziiHAyWdvbr1x2Lfj7ZKn7ha9VmOMF/06U0bSuyVL+WeMxBvUkWr3dOpPazKMj5mF3Mjss1TXZFlyMOHNU8a55zxH4454j9qthDdpQrVqoUxLDPppOdC3EQT6Dl3tnLPzgFtqgETmFKcoSD9+nt5Dlvwz59uuHxK4/HWEPI8z0RLoSxe30ec//zdPSy5B8/ev/Amu6LBkUzR2q0Vj4uC+fA3wf5teA6BzF/h70vAsCsSn/exwCLuIwPdXBvvV8GAPz/7Z17jFT1Fcc/h31QebgqwrouykJVgmIB2ajoKrVKK7QRTWvQWFmFtNhipaZ/VPpIY1pTbOpbEiUVYmsfamtTLBaiTTWx1cqSijyWh1hSMasilYeaqiSnf9zfuHd35t65M3PvnTvM+SSTuXvmN7/7nbN35tzf6/zuuHIK33l8Y0Wzz0qh7oJDiqvPM0Fu6mAlKS/iJK6Wgz/FQKWtmkJ79voZ2HII/hFYd3PYLrnRCdo+NatEvbZaiqwNmjB6BP/+6ZyK/5+56ydszOGN/V4CwJ1vDdx7IjfmECUdCsAT3zyPtZvfDByIh/5gFTaNN3eNhQWH3A1DsSSQcVF/wcE9Z+RGOnFGDm1k0YUTuDwkb1IpVLKgCgZOCwxa7HOySwgX1i1wypiRFekohTWb+j45TnpmVJw0NwzhnAnHJX6eKFNQoxJH9+XgbT0LsXaz9z9ds6mP5T57bswhaGOjwUxub2Fye0tomZyOMD25wBEWHJobc2tFIkmrmLoLDrMnn8BprSPKTvFca4gIS+dMiqWuh7o7mXhCZT/K/i9d0BehZVgTu5d9sWhdl0xq5Znet2LLYRWFtPp7i3HrZWcUHZTccdvsVLSICJ3jjuXakOR8afJhkV3gwMsc+/yr7+RNmz34v4/Z9uYhXnht3ye2sPGEKFw+rZ3evoPcfMlpgWVyXWVhGROaG7wykWdSVUjdBYdxo4aHLsG/qEgWz3rm4kmtFdfhbzlEndMdxH1XT2P3vvcrmj1VCiKVrTaOk+4KUqEkwe+/cV61JXxCriUTlk57+TVn8dj617nu/I4B9m0FNvcJS0gZhU81NXDr3PDEjG0t3mB1WGu5tcW7mRo3Knqq9UrIxpUOiMilwD1AA/ALVV2WtobtP7k00rQ8I5hHFp4TmnfqxotOYeHDPbGc66jmhlQHbrtndKR2LqN8TnJjSGFdoC1HNfG1kGy+4O2D/tSSC2LVFsSCrg6GSPhg/NDGBh5bNCMwpXfcZCI4iEgDsByYBewB1ovIalXdmqaOeulqSpKuItlJ42h9VItvxZSt10iWSyefwF3zpjDnzOBps1E4d0LhpH1JMLSxgUUzi890O3t88mNIObJym3w28KqqvqaqHwG/A+ZWWZOREGtu6uLFpRdXW0ZkjncrwitJOGikh4hwxbSxZd3sXeGbuDF+dPEMsEcymWg5AO3A676/9wDnVEmLkTBnnBg+uyNrPLpoBms3v5lKygKjutw1byp3zZtKb99BJramNyMui2TlVqjQtIK80SQR+bqI9IhIz969e1OQZRjw6dEjWHyRdSnVE5Pajg7NoloPZCU47AH8IzFjgbxVKKq6QlU7VbVz9GibVWQYhpEUWQkO64FTRWS8iDQDVwGrq6zJMAyjbsnEmIOqHhaRG4F1eFNZV6rqlirLMgzDqFsyERwAVPUp4Klq6zAMwzCy061kGIZhZAgLDoZhGEYeFhwMwzCMPCw4GIZhGHlIWhtHxI2IHAK2B7x8MvCfIlW0AAcyVMY0p1OmXjVHOU/UcnFojus8UcrZdTGQiapafPm3qtbkA+gJeW1vhPevyFgZ02yaE9Mc5Txpao7rPHF99nq6LsJ+O/2PI7VbaX+EMk9mrIxpTqdMvWqOcp6o5eLQHNd5opSz66IMarlbqUdVO0t9LauY5nQwzelQa5prTS+Urznq+2q55bCizNeyimlOB9OcDrWmudb0QvmaI72vZlsOhmEYRnLUcsvBMAzDSIiaCA4islJE3haRzT7bFBF5QUQ2iciTInK0szeLyCpn3ygin/W9Z7qzvyoi94pIYgnbY9T8rIhsF5GX3WNMgppPEpG/iUiviGwRkSXOfpyIPC0iO93zsb73LHX+3C4iX/DZU/F1zJpT8XWpmkVklCv/nojcP6iuTPq5iObE/VyG3lkissH5coOIfM5XV1Z9HKa5ch9HmdJU7QdwIXAWsNlnWw/MdMcLgB+748XAKnc8BtgADHF/vwTMwNtc6C/A7BrQ/CzQmZKf24Cz3PFIYAdwOvAz4BZnvwW43R2fDmwEhgLjgV1AQ5q+jllzKr4uQ/NwoAu4Abh/UF1Z9XOY5sT9XIbeacCJ7ngy8EYN+DhMc8U+TvRLELPjOhj4Q3uQ/jGTk4Ct7ng58FVfub/i7VHdBmzz2a8GHsyy5rj+yRXo/xMwC2+xYZuztQHb3fFSYKmv/Dr3JUrd15Vqrqavi2n2lbsO3w9tlv0cpLlafo6q19kF2Id3A5F5Hw/WHJePa6JbKYDNwGXu+Er6d5LbCMwVkUYRGQ9Md6+14+04l2OPs6VJqZpzrHJNwx8m1aQdjIh04N2Z/BNoVdU+APeca6IW2vu7nSr5ukLNOVL1dUTNQWTZz8VIzc9l6P0y8C9V/ZDa8bFfc46KfFzLwWEBsFhENuA1wT5y9pV4/8Ae4G7gH8BhIu5TnTClaga4RlXPBC5wj2uTFikiI4A/AN9W1YNhRQvYNMSeGDFohpR9XYLmwCoK2LLi5zBS83OpekXkDOB2YFHOVKBYpnxcQDPE4OOaDQ6quk1VP6+q04Hf4vUdo6qHVfVmVZ2qqnOBY4CdeD++Y31VFNynOmOaUdU33PMh4Dd4XWSJISJNeBfmr1X1CWd+S0Ta3OttwNvOHrT3d6q+jklzqr4uUXMQWfZzIGn5uVS9IjIW+CMwX1V3OXOmfRygORYf12xwyI2+i8gQ4AfAA+7vYSIy3B3PAg6r6lbXHDskIue6JtZ8vD69zGp23UzHO3sT8CW8rqmk9AnwENCrqnf6XloNdLvjbvr9thq4SkSGuu6wU4GX0vR1XJrT9HUZmguScT8H1ZOKn0vVKyLHAGvwxqP+niucZR8HaY7Nx2kMrFT6wLvL7gM+xovkC4EleKP5O4Bl9A/0duAN4PQCzwDjfPV0OiftAu7PvSermvFmfGwAXgG2APfgZtYkpLkLr8n8CvCye8wBRuENku90z8f53vN958/t+GZxpOXruDSn6esyNe8G/gu8566n02vAz3ma0/JzqXrxbtbe95V9GRiTZR8HaY7Lx7ZC2jAMw8ijZruVDMMwjOSw4GAYhmHkYcHBMAzDyMOCg2EYhpGHBQfDMAwjDwsOhpEAInKDiMwvoXyH+DL4Gka1aay2AMM40hCRRlV9oNo6DKMSLDgYRgFc4rO1eInPpuEtXJwPTALuBEYA7wDXqWqfiDyLlxPrfGC1iIwE3lPVn4vIVLzV8MPwFlItUNV3RWQ6Xl6tD4Dn0/t0hlEc61YyjGAmAitU9TN46dYXA/cBX1EvP9ZK4DZf+WNUdaaq3jGonl8C33X1bAJ+5OyrgJtUdUaSH8IwysFaDoYRzOvan7PmEeB7eJuqPO0yIDfgpUjJ8ejgCkSkBS9oPOdMDwOPF7D/Cpgd/0cwjPKw4GAYwQzOLXMI2BJyp/9+CXVLgfoNIzNYt5JhBHOyiOQCwdXAi8DonE1Emlwu/UBU9QDwrohc4EzXAs+p6n7ggIh0Ofs18cs3jPKxloNhBNMLdIvIg3gZMe/D21b0Xtct1Ii3OdOWIvV0Aw+IyDDgNeB6Z78eWCkiH7h6DSMzWFZWwyiAm630Z1WdXGUphlEVrFvJMAzDyMNaDoZhGEYe1nIwDMMw8rDgYBiGYeRhwcEwDMPIw4KDYRiGkYcFB8MwDCMPCw6GYRhGHv8HIxW7SzTeDPYAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sorted_data['inc'].plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Incidence annuelle" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [], "source": [ "first_august_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": 26, "metadata": {}, "outputs": [], "source": [ "year = []\n", "yearly_incidence = []\n", "for week1, week2 in zip(first_august_week[:-1],\n", " first_august_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": 27, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 27, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAHARJREFUeJzt3X9wndV95/H3x8iRDLWJDDLxD8DuxGFiSGriO7a7yXZL3Npm08FmCok2NGgWz7gBNiWdnQE7sMMsuDPQ2dm0nkxYPKFgIPxw3TJ4u3hBmM3U2zq25UAKhrhSAjEODhIrB5vMWEX2d/+4R/WVkKXnSlf3lz6vmTv3uV895+g8fnT9vc855zlXEYGZmVkWUyrdADMzqx1OGmZmlpmThpmZZeakYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmDZVuQKldeOGFMX/+/Eo3w8ysphw4cOC9iGgZbb+6Sxrz58+no6Oj0s0wM6spkn6eZT93T5mZWWZOGmZmlpmThpmZZeakYWZmmTlpmJlZZk4aZjWg+/hJvvzgHrpPnKx0U2ySc9IwqwGbd3Wy/61eNr/YWemm2CRXd/dpmNWTy+7aSV//6X99/fjewzy+9zCNDVM4tOnqCrbMJitfaRRwF4BVm923X8U1i+fQNDX/Vm2aOoU1i+ew+46rKtwym6ycNAq4C8CqzawZTUxvbKCv/zSNDVPo6z/N9MYGZk1v+si+/tBj5eDuKdwFYNXtvQ/6uGHZpXx16SU8se8wPWdJCoUfejZd+5kyt9ImC0VEpdtQUrlcLopde6r7+Ek2PfcGLxz8JSc/PE3T1CmsuvwT3PmlTw/7ic6smgz90DPAH3qsGJIORERutP3cPUVxXQBm1cbjHlZO7p5KsnYBmFUbf+ixcnLSSB782pmrsk1rr6hgS8yK5w89Vi4e0zAzM49pmJlZ6TlpmJlZZpmShqQ/lXRQ0muSnpTUJGmmpHZJnem5uWD/jZK6JB2StKogvkTSq+lnmyUpxRslPZ3ieyXNLyjTln5Hp6S20h26mZkVa9SkIWku8CdALiKuAM4BWoENwK6IWAjsSq+RtCj9/HJgNfBdSeek6h4A1gML02N1iq8DjkXEJ4FvA/enumYCdwPLgKXA3YXJyczMyitr91QDME1SA3Au8A6wBtiafr4VWJu21wBPRURfRLwJdAFLJc0GZkTEnsiPvj86pMxAXduBFekqZBXQHhG9EXEMaOdMojEzszIbNWlExC+A/wYcBo4C70fEC8BFEXE07XMUmJWKzAXeLqjiSIrNTdtD44PKREQ/8D5wwQh1mZlZBWTpnmomfyWwAJgDnCfpj0YqMkwsRoiPtUxhG9dL6pDU0dPTM0LTzMxsPLJ0T/0e8GZE9ETEh8DfAv8GeDd1OZGeu9P+R4CLC8rPI9+ddSRtD40PKpO6wM4Hekeoa5CI2BIRuYjItbS0ZDgkMzMbiyxJ4zCwXNK5aZxhBfAGsAMYmM3UBjybtncArWlG1ALyA977UhfWCUnLUz03DikzUNd1wEtp3ON5YKWk5nTFszLFzGwS8zLwlTPqMiIRsVfSduBHQD/wMrAF+A1gm6R15BPL9Wn/g5K2Aa+n/W+NiFOpupuBR4BpwM70AHgIeExSF/krjNZUV6+ke4H9ab97IqJ3XEdsZjXPy8BXjpcRsbrWffwk/+nJl/nOV6/0An51wMvATxwvI2KGv42x3ngZ+MrzKrdWl/xtjPXJy8BXnq80rC75E2n9GlgG/plbPs8Nyy6l54O+SjdpUvGVhtUlfyKtX/7um8py0rC65S8mMis9z56ymuWZUWal49lTNm7VfgOVZ0aZlZ+7p+ysqvUGKs+MMqscd0/ZR1T7DVTdx0+y6bk3eOHgLzn54Wmapk5h1eWf4M4vfbrmuqncxWbVwt1TNmbVPl21nmZGuYvNao27p+wjauE/5VqfGeUuNqtVTho2rGr/T7nW5+rvvv2qs3axmVUzJw0bVq3/p1ztauFqzmw4ThpmFVLtV3Nmw/HsKTMz8+wpMzMrPScNMzPLzEnDzMwyc9IwmwDVvm6X2Vg5aZhNAN/pbfXKU27NSsh3elu985WGWQlV+7pdZuPlpGFWQr7T2+rdqElD0mWSXil4HJf0TUkzJbVL6kzPzQVlNkrqknRI0qqC+BJJr6afbZakFG+U9HSK75U0v6BMW/odnZLaSnv4ZqU3cKf3M7d8nhuWXUrPB32VbpJZyRR1R7ikc4BfAMuAW4HeiLhP0gagOSLukLQIeBJYCswBXgQ+FRGnJO0DbgN+CDwHbI6InZJuAT4bEV+X1ApcGxFfkTQT6AByQAAHgCURcexsbfQd4Wbl5e8EqQ8TdUf4CuCnEfFzYA2wNcW3AmvT9hrgqYjoi4g3gS5gqaTZwIyI2BP5TPXokDIDdW0HVqSrkFVAe0T0pkTRDqwuss1mNoE8U2xyKXb2VCv5qwiAiyLiKEBEHJU0K8Xnkr+SGHAkxT5M20PjA2XeTnX1S3ofuKAwPkwZM6sgzxSbnDJfaUj6GHAN8Nej7TpMLEaIj7VMYdvWS+qQ1NHT0zNK88ysFDxTbHIqpnvqauBHEfFuev1u6nIiPXen+BHg4oJy84B3UnzeMPFBZSQ1AOcDvSPUNUhEbImIXETkWlpaijgkMxurYmeK+S75+lBM0vgPnOmaAtgBDMxmagOeLYi3phlRC4CFwL7UlXVC0vI0XnHjkDIDdV0HvJTGPZ4HVkpqTrOzVqaYnYXfmFZOxcwU89hHfcg0e0rSueTHFn4zIt5PsQuAbcAlwGHg+ojoTT+7E7gJ6Ae+GRE7UzwHPAJMA3YC34iIkNQEPAZcSf4KozUifpbK3AR8KzXlzyLi4ZHaOtlnT931zKt8f99hblh6CZuu/Uylm2P2kbGPAR77qC5ZZ0/5S5jqRLFvTE+TtHLpPn7yrN+H7r+96uEvYZpkih2UdFeBlYvvkq8vXrCwTmR9Y3qapFWCvw+9fjhp1JEsb8zdt1911q4Cs4ny4NfO9HpsWntFBVti4+WkUUeyvDHdVWBm4+ExjQlSzVNfvaCemY2VZ09NEE99NbNaknX2lLunSswDzWZWz9w9VWLFTn2t5m4sM7OhnDRKrNiBZt8vYWa1xN1TEyDL1Fd3Y5lZLfKVRpGydCc9+LUcm9ZewaI5M9i09opBU2EHeFlpqxR3idaWajtfThpFKlV3ku+XsEpxl2htqbbz5Sm3GU3ESp1//FgHLdObBnVjDXdVYlYKXm22tpT7fHmV2xLzSp0Gtb06sP+Ga0u5z5dXuS0xdyfVt6z9xtXWVVAM/w3Xlmo9X549VQSv1Fm/CpPBcHfw18tsN/8N15ZqPF/unrJJLWu/8WTs2qnlrjgrnrunrOpU29RByD71uVq7CiZSLXfF2cRx95SNW9ZPpKN1AVVCMcmgGrsKJkK9dMXZxHD3lI3baCv6VvtUT099HmwydsWZV7m1Msj6ibTavy3Q3yo32GTsirPsPKZhY+bxgPrlL+qys/GVho2ZxwPqV6Wuvjxjq/plutKQ9HFJ2yX9RNIbkn5b0kxJ7ZI603Nzwf4bJXVJOiRpVUF8iaRX0882S1KKN0p6OsX3SppfUKYt/Y5OSW2lO3QrhayfSLMs4mjmGVvVL9NAuKStwO6I+J6kjwHnAt8CeiPiPkkbgOaIuEPSIuBJYCkwB3gR+FREnJK0D7gN+CHwHLA5InZKugX4bER8XVIrcG1EfEXSTKADyAEBHACWRMSxs7XVA+FmtafaJ0tMBiW7T0PSDOB3gIcAIuJfIuJXwBpga9ptK7A2ba8BnoqIvoh4E+gClkqaDcyIiD2Rz1SPDikzUNd2YEW6ClkFtEdEb0oU7cDq0dpsZrXFXxVQO7J0T/0m0AM8LOllSd+TdB5wUUQcBUjPs9L+c4G3C8ofSbG5aXtofFCZiOgH3gcuGKGuQSStl9QhqaOnpyfDIZlZNfFkidqRJWk0AJ8DHoiIK4FfAxtG2F/DxGKE+FjLnAlEbImIXETkWlpaRmiamVUrz9iqDVlmTx0BjkTE3vR6O/mk8a6k2RFxNHU9dRfsf3FB+XnAOyk+b5h4YZkjkhqA84HeFP/dIWV+kOnIzKym1ML9Mp7dleFKIyJ+Cbwt6bIUWgG8DuwABmYztQHPpu0dQGuaEbUAWAjsS11YJyQtT+MVNw4pM1DXdcBLadzjeWClpOY0O2tlipmZlZ1nd2W/T+MbwPfTzKmfAf+RfMLZJmkdcBi4HiAiDkraRj6x9AO3RsSpVM/NwCPANGBnekB+kP0xSV3krzBaU129ku4F9qf97omI3jEeq5nZmHg9rjO89pSZ2Sgmw3pcXhrdzKxEPLvrDC8jYmaWgZfCyXP3lJmZuXuq3lTjt96Z2eTjpFEjPNXPzKqBxzSqnKf6mVk18ZVGlfNCbmZWTZw0qpyn+pmNnccCS89JowZ4ITezsfFYYOl5yq2Z1R1/qVPxPOXWzCatYscC3Y2VnZOGmdWdYscC3Y2VnafcmlldyrLsh6e0F89jGmY2aU2G1Wuz8piGmdkoPKW9eO6eMrNJzavXFsfdU2Zm5u4pMzMrPScNqzqeM29WvZw0rOp4zrxZ9fJAuFUNz5k3q36+0rCq4WXgzaqfk4ZVDc+ZN6t+mZKGpLckvSrpFUkdKTZTUrukzvTcXLD/Rkldkg5JWlUQX5Lq6ZK0WZJSvFHS0ym+V9L8gjJt6Xd0Smor1YFbdfIy8GbVLdN9GpLeAnIR8V5B7M+B3oi4T9IGoDki7pC0CHgSWArMAV4EPhURpyTtA24Dfgg8B2yOiJ2SbgE+GxFfl9QKXBsRX5E0E+gAckAAB4AlEXHsbG31fRpmZsUrx30aa4CtaXsrsLYg/lRE9EXEm0AXsFTSbGBGROyJfKZ6dEiZgbq2AyvSVcgqoD0ielOiaAdWj6PNZmY2DlmTRgAvSDogaX2KXRQRRwHS86wUnwu8XVD2SIrNTdtD44PKREQ/8D5wwQh1DSJpvaQOSR09PT0ZD8nMzIqVdcrt5yPiHUmzgHZJPxlhXw0TixHiYy1zJhCxBdgC+e6pEdpmZmbjkOlKIyLeSc/dwDPkxyveTV1OpOfutPsR4OKC4vOAd1J83jDxQWUkNQDnA70j1GVmZhUwatKQdJ6k6QPbwErgNWAHMDCbqQ14Nm3vAFrTjKgFwEJgX+rCOiFpeRqvuHFImYG6rgNeSuMezwMrJTWn2VkrU8zMzCogS/fURcAzaXZsA/BERPxvSfuBbZLWAYeB6wEi4qCkbcDrQD9wa0ScSnXdDDwCTAN2pgfAQ8BjkrrIX2G0prp6Jd0L7E/73RMRveM4XjMzGwcvjW5mZl4a3czMSs9Jw8zMMnPSMDOzzJw0zMwsMycNMzPLzEnDzMwyc9IwM7PMnDTMzCwzJw0zM8vMScPMzDJz0jAzs8ycNMzMLDMnDTMzy8xJw8zMMnPSMDOzzJw0zMwsMycNMzPLzEnDzMwyc9IwM7PMnDTMzCwzJw0zM8vMScPMzDLLnDQknSPpZUl/l17PlNQuqTM9Nxfsu1FSl6RDklYVxJdIejX9bLMkpXijpKdTfK+k+QVl2tLv6JTUVoqDNjOzsSnmSuM24I2C1xuAXRGxENiVXiNpEdAKXA6sBr4r6ZxU5gFgPbAwPVan+DrgWER8Evg2cH+qayZwN7AMWArcXZiczMysvDIlDUnzgC8B3ysIrwG2pu2twNqC+FMR0RcRbwJdwFJJs4EZEbEnIgJ4dEiZgbq2AyvSVcgqoD0ieiPiGNDOmURjZmZllvVK4y+A24HTBbGLIuIoQHqeleJzgbcL9juSYnPT9tD4oDIR0Q+8D1wwQl1mZlYBoyYNSX8AdEfEgYx1aphYjBAfa5nCNq6X1CGpo6enJ2MzzcysWFmuND4PXCPpLeAp4IuSHgfeTV1OpOfutP8R4OKC8vOAd1J83jDxQWUkNQDnA70j1DVIRGyJiFxE5FpaWjIckpmZjcWoSSMiNkbEvIiYT36A+6WI+CNgBzAwm6kNeDZt7wBa04yoBeQHvPelLqwTkpan8Yobh5QZqOu69DsCeB5YKak5DYCvTDEzM6uAhnGUvQ/YJmkdcBi4HiAiDkraBrwO9AO3RsSpVOZm4BFgGrAzPQAeAh6T1EX+CqM11dUr6V5gf9rvnojoHUebzcxsHJT/QF8/crlcdHR0VLoZZmY1RdKBiMiNtp/vCDczq4Du4yf58oN76D5xstJNKYqThplZBWze1cn+t3rZ/GJnpZtSlPGMaZiZWZEuu2snff1nbnl7fO9hHt97mMaGKRzadHUFW5aNrzTMzMpo9+1Xcc3iOTRNzf/32zR1CmsWz2H3HVeNq95ydXc5aZiZldGsGU1Mb2ygr/80jQ1T6Os/zfTGBmZNbxpXveXq7nL3lJlZmb33QR83LLuUry69hCf2HaZnHFcH5e7u8pRbM7Ma1n38JJuee4MXDv6Skx+epmnqFFZd/gnu/NKni7p68ZRbM7NJYKK6u87G3VNmZjWulN1do3H3lJmZuXvKzMxKz0nDzMwyc9IwM7PMnDTMzCwzJw0zM8vMScPMzDJz0jAzs8ycNMzMLDMnDTMzy8xJw8zMMnPSMDOzzJw0zMwsMycNMzPLbNSkIalJ0j5JP5Z0UNJ/TfGZktoldabn5oIyGyV1STokaVVBfImkV9PPNktSijdKejrF90qaX1CmLf2OTkltpTx4MzMrTpYrjT7gixHxW8BiYLWk5cAGYFdELAR2pddIWgS0ApcDq4HvSjon1fUAsB5YmB6rU3wdcCwiPgl8G7g/1TUTuBtYBiwF7i5MTmZmVl6jJo3I+yC9nJoeAawBtqb4VmBt2l4DPBURfRHxJtAFLJU0G5gREXsi/yUejw4pM1DXdmBFugpZBbRHRG9EHAPaOZNozMyszDKNaUg6R9IrQDf5/8T3AhdFxFGA9Dwr7T4XeLug+JEUm5u2h8YHlYmIfuB94IIR6jIzswrIlDQi4lRELAbmkb9quGKE3TVcFSPEx1rmzC+U1kvqkNTR09MzQtPMzGw8ipo9FRG/An5Avovo3dTlRHruTrsdAS4uKDYPeCfF5w0TH1RGUgNwPtA7Ql1D27UlInIRkWtpaSnmkMzMrAhZZk+1SPp42p4G/B7wE2AHMDCbqQ14Nm3vAFrTjKgF5Ae896UurBOSlqfxihuHlBmo6zrgpTTu8TywUlJzGgBfmWJmZlYBDRn2mQ1sTTOgpgDbIuLvJO0BtklaBxwGrgeIiIOStgGvA/3ArRFxKtV1M/AIMA3YmR4ADwGPSeoif4XRmurqlXQvsD/td09E9I7ngM3MbOyU/0BfP3K5XHR0dFS6GWZmNUXSgYjIjbaf7wg3M7PMnDTMzCwzJw0zM8vMScPMzDJz0jAzs8ycNMzMLDMnDTMzy8xJw8zMMnPSMDOzzJw0zMwsMycNMzPLzEnDzMwyc9IwM7PMnDTMzCwzJw0zM8vMScPMzDJz0jAzs8ycNMzMSqz7+Em+/OAeuk+crHRTSs5Jw8ysxDbv6mT/W71sfrGz0k0puYZKN8DMrF5cdtdO+vpP/+vrx/ce5vG9h2lsmMKhTVdXsGWl4ysNM7MS2X37VVyzeA5NU/P/tTZNncKaxXPYfcdVFW5Z6ThpmJmVyKwZTUxvbKCv/zSNDVPo6z/N9MYGZk1vqnTTSsbdU2ZmJfTeB33csOxSvrr0Ep7Yd5ieOhsMV0SMvIN0MfAo8AngNLAlIv5S0kzgaWA+8Bbw5Yg4lspsBNYBp4A/iYjnU3wJ8AgwDXgOuC0iQlJj+h1LgP8HfCUi3kpl2oC7UnM2RcTWkdqby+Wio6Mj+7+AmZkh6UBE5EbbL0v3VD/wnyPi08By4FZJi4ANwK6IWAjsSq9JP2sFLgdWA9+VdE6q6wFgPbAwPVan+DrgWER8Evg2cH+qayZwN7AMWArcLak5Q5vNzGwCjJo0IuJoRPwobZ8A3gDmAmuAgU/9W4G1aXsN8FRE9EXEm0AXsFTSbGBGROyJ/OXNo0PKDNS1HVghScAqoD0ietNVTDtnEo2ZmZVZUQPhkuYDVwJ7gYsi4ijkEwswK+02F3i7oNiRFJubtofGB5WJiH7gfeCCEeoyM7MKyJw0JP0G8DfANyPi+Ei7DhOLEeJjLVPYtvWSOiR19PT0jNA0MzMbj0xJQ9JU8gnj+xHxtyn8bupyIj13p/gR4OKC4vOAd1J83jDxQWUkNQDnA70j1DVIRGyJiFxE5FpaWrIckpmZjcGoSSONLTwEvBER/73gRzuAtrTdBjxbEG+V1ChpAfkB732pC+uEpOWpzhuHlBmo6zrgpTTu8TywUlJzGgBfmWJmZlYBWabcfgHYDbxKfsotwLfIj2tsAy4BDgPXR0RvKnMncBP5mVffjIidKZ7jzJTbncA30pTbJuAx8uMlvUBrRPwslbkp/T6AP4uIh0dpbw/w84zHXy0uBN6rdCPGycdQHXwM1aPWjuPSiBi1q2bUpGETT1JHlvnR1czHUB18DNWjXo5jKC8jYmZmmTlpmJlZZk4a1WFLpRtQAj6G6uBjqB71chyDeEzDzMwy85WGmZll5qQxAST9laRuSa8VxH5L0h5Jr0r6n5JmpPjHJD2c4j+W9LsFZX4g6ZCkV9Jj1jC/bqKO4WJJ/0fSG5IOSrotxWdKapfUmZ6bC8pslNSV2ryqIL4kHV+XpM3pPp1aO4aKnItij0HSBWn/DyR9Z0hdNXEeRjmGmnlPSPp9SQfSv/kBSV8sqKsi56IkIsKPEj+A3wE+B7xWENsP/Lu0fRNwb9q+FXg4bc8CDgBT0usfALkKHcNs4HNpezrwz8Ai4M+BDSm+Abg/bS8Cfgw0AguAnwLnpJ/tA36b/LIwO4Gra/AYKnIuxnAM5wFfAL4OfGdIXbVyHkY6hlp6T1wJzEnbVwC/qPS5KMXDVxoTICL+nvxNioUuA/4+bbcDf5i2F5FfWp6I6AZ+BVR8bneUZ3XjmjiGcrT1bIo9hoj4dUT8X2DQN//U0nk42zFU2hiO4+WIGFj26CDQpPxKGRU7F6XgpFE+rwHXpO3rObOm1o+BNZIalF92ZQmD19t6OF2G/5dKXcJq4lY3LptxHsOAip6LjMdwNrV0HkZTK++JQn8IvBwRfVTJuRgrJ43yuYn8F1gdIH9p+y8p/lfk/2g6gL8A/pH88isAN0TEZ4B/mx5fK2uLmfDVjcuiBMcAFT4XRRzDWasYJlat52EktfSeGNj/cvJfLPfHA6FhdquZaaxOGmUSET+JiJURsQR4knx/ORHRHxF/GhGLI2IN8HGgM/3sF+n5BPAEZe4q0cSvbjzhSnQMFT0XRR7D2dTSeTirGntPIGke8AxwY0T8NIUrei7Gy0mjTAZmeUiaQv47z/9Hen2upPPS9u8D/RHxeuquujDFpwJ/QL6Lq1ztLcfqxjVxDJU8F2M4hmHV2Hk4Wz019Z6Q9HHgfwEbI+IfBnau5LkoiUqPxNfjg/yVxFHgQ/KfKtYBt5GfbfHPwH2cubFyPnCI/KDai+RXmoT8DJIDwD+RH0T7S9JMnjIdwxfIXzL/E/BKevx78t+ouIv81dAuYGZBmTvJX0EdomA2CPmB/dfSz74zcOy1cgyVPBdjPIa3yE/E+CD9/S2qwfPwkWOotfcE+Q+Hvy7Y9xVgViXPRSkeviPczMwyc/eUmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZWWZOGmZmltn/B3v28ZbrUIyxAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "yearly_incidence.plot(style='*')" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "2020 221186\n", "2023 366227\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 641397\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": 28, "metadata": {}, "output_type": "execute_result" } ], "source": [ "yearly_incidence.sort_values()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }