{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Analyse Varicelle" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import isoweek" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "data_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": 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
02024427310413324876528FRFrance
1202441720103613659315FRFrance
2202440721257253525315FRFrance
32024397289813334463426FRFrance
4202438775101513102FRFrance
52024377916281804102FRFrance
6202436722358703600315FRFrance
7202435716202852955204FRFrance
8202434725606224498417FRFrance
9202433719715363406315FRFrance
1020243274399194468547311FRFrance
1120243174500221367877410FRFrance
12202430770044278973011715FRFrance
1320242979270630312237141018FRFrance
1420242879364649812230141018FRFrance
15202427710247709013404151020FRFrance
162024267143681039918337221628FRFrance
17202425711174803914309171222FRFrance
18202424712621935715885191424FRFrance
192024237146571133917975221727FRFrance
20202422711628836114895171222FRFrance
2120242179701685112551151119FRFrance
222024207136611020917113201525FRFrance
2320241971008364131375315921FRFrance
24202418713438951417362201426FRFrance
252024177153031121919387231729FRFrance
262024167181381354022736272034FRFrance
272024157249291731532543372648FRFrance
282024147161811254419818241929FRFrance
292024137183221420622438272133FRFrance
.................................
17381991267176081130423912312042FRFrance
17391991257161691070021638281838FRFrance
17401991247161711007122271281739FRFrance
1741199123711947767116223211329FRFrance
1742199122715452995320951271737FRFrance
1743199121714903897520831261636FRFrance
17441991207190531274225364342345FRFrance
17451991197167391124622232291939FRFrance
17461991187213851388228888382551FRFrance
1747199117713462887718047241632FRFrance
17481991167148571006819646261834FRFrance
1749199115713975978118169251832FRFrance
1750199114712265768416846221430FRFrance
175119911379567604113093171123FRFrance
1752199112710864733114397191325FRFrance
17531991117155741118419964271935FRFrance
17541991107166431137221914292038FRFrance
1755199109713741878018702241533FRFrance
1756199108713289881317765231531FRFrance
1757199107712337807716597221529FRFrance
1758199106710877701314741191226FRFrance
1759199105710442654414340181125FRFrance
17601991047791345631126314820FRFrance
17611991037153871048420290271836FRFrance
17621991027162771104621508292038FRFrance
17631991017155651027120859271836FRFrance
17641990527193751329525455342345FRFrance
17651990517190801380724353342543FRFrance
1766199050711079666015498201228FRFrance
17671990497114302610205FRFrance
\n", "

1768 rows × 10 columns

\n", "
" ], "text/plain": [ " week indicator inc inc_low inc_up inc100 inc100_low \\\n", "0 202442 7 3104 1332 4876 5 2 \n", "1 202441 7 2010 361 3659 3 1 \n", "2 202440 7 2125 725 3525 3 1 \n", "3 202439 7 2898 1333 4463 4 2 \n", "4 202438 7 751 0 1513 1 0 \n", "5 202437 7 916 28 1804 1 0 \n", "6 202436 7 2235 870 3600 3 1 \n", "7 202435 7 1620 285 2955 2 0 \n", "8 202434 7 2560 622 4498 4 1 \n", "9 202433 7 1971 536 3406 3 1 \n", "10 202432 7 4399 1944 6854 7 3 \n", "11 202431 7 4500 2213 6787 7 4 \n", "12 202430 7 7004 4278 9730 11 7 \n", "13 202429 7 9270 6303 12237 14 10 \n", "14 202428 7 9364 6498 12230 14 10 \n", "15 202427 7 10247 7090 13404 15 10 \n", "16 202426 7 14368 10399 18337 22 16 \n", "17 202425 7 11174 8039 14309 17 12 \n", "18 202424 7 12621 9357 15885 19 14 \n", "19 202423 7 14657 11339 17975 22 17 \n", "20 202422 7 11628 8361 14895 17 12 \n", "21 202421 7 9701 6851 12551 15 11 \n", "22 202420 7 13661 10209 17113 20 15 \n", "23 202419 7 10083 6413 13753 15 9 \n", "24 202418 7 13438 9514 17362 20 14 \n", "25 202417 7 15303 11219 19387 23 17 \n", "26 202416 7 18138 13540 22736 27 20 \n", "27 202415 7 24929 17315 32543 37 26 \n", "28 202414 7 16181 12544 19818 24 19 \n", "29 202413 7 18322 14206 22438 27 21 \n", "... ... ... ... ... ... ... ... \n", "1738 199126 7 17608 11304 23912 31 20 \n", "1739 199125 7 16169 10700 21638 28 18 \n", "1740 199124 7 16171 10071 22271 28 17 \n", "1741 199123 7 11947 7671 16223 21 13 \n", "1742 199122 7 15452 9953 20951 27 17 \n", "1743 199121 7 14903 8975 20831 26 16 \n", "1744 199120 7 19053 12742 25364 34 23 \n", "1745 199119 7 16739 11246 22232 29 19 \n", "1746 199118 7 21385 13882 28888 38 25 \n", "1747 199117 7 13462 8877 18047 24 16 \n", "1748 199116 7 14857 10068 19646 26 18 \n", "1749 199115 7 13975 9781 18169 25 18 \n", "1750 199114 7 12265 7684 16846 22 14 \n", "1751 199113 7 9567 6041 13093 17 11 \n", "1752 199112 7 10864 7331 14397 19 13 \n", "1753 199111 7 15574 11184 19964 27 19 \n", "1754 199110 7 16643 11372 21914 29 20 \n", "1755 199109 7 13741 8780 18702 24 15 \n", "1756 199108 7 13289 8813 17765 23 15 \n", "1757 199107 7 12337 8077 16597 22 15 \n", "1758 199106 7 10877 7013 14741 19 12 \n", "1759 199105 7 10442 6544 14340 18 11 \n", "1760 199104 7 7913 4563 11263 14 8 \n", "1761 199103 7 15387 10484 20290 27 18 \n", "1762 199102 7 16277 11046 21508 29 20 \n", "1763 199101 7 15565 10271 20859 27 18 \n", "1764 199052 7 19375 13295 25455 34 23 \n", "1765 199051 7 19080 13807 24353 34 25 \n", "1766 199050 7 11079 6660 15498 20 12 \n", "1767 199049 7 1143 0 2610 2 0 \n", "\n", " inc100_up geo_insee geo_name \n", "0 8 FR France \n", "1 5 FR France \n", "2 5 FR France \n", "3 6 FR France \n", "4 2 FR France \n", "5 2 FR France \n", "6 5 FR France \n", "7 4 FR France \n", "8 7 FR France \n", "9 5 FR France \n", "10 11 FR France \n", "11 10 FR France \n", "12 15 FR France \n", "13 18 FR France \n", "14 18 FR France \n", "15 20 FR France \n", "16 28 FR France \n", "17 22 FR France \n", "18 24 FR France \n", "19 27 FR France \n", "20 22 FR France \n", "21 19 FR France \n", "22 25 FR France \n", "23 21 FR France \n", "24 26 FR France \n", "25 29 FR France \n", "26 34 FR France \n", "27 48 FR France \n", "28 29 FR France \n", "29 33 FR France \n", "... ... ... ... \n", "1738 42 FR France \n", "1739 38 FR France \n", "1740 39 FR France \n", "1741 29 FR France \n", "1742 37 FR France \n", "1743 36 FR France \n", "1744 45 FR France \n", "1745 39 FR France \n", "1746 51 FR France \n", "1747 32 FR France \n", "1748 34 FR France \n", "1749 32 FR France \n", "1750 30 FR France \n", "1751 23 FR France \n", "1752 25 FR France \n", "1753 35 FR France \n", "1754 38 FR France \n", "1755 33 FR France \n", "1756 31 FR France \n", "1757 29 FR France \n", "1758 26 FR France \n", "1759 25 FR France \n", "1760 20 FR France \n", "1761 36 FR France \n", "1762 38 FR France \n", "1763 36 FR France \n", "1764 45 FR France \n", "1765 43 FR France \n", "1766 28 FR France \n", "1767 5 FR France \n", "\n", "[1768 rows x 10 columns]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "raw_data = pd.read_csv(data_url, skiprows=1)\n", "raw_data" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On a enlevé la première ligne du fichier csv car c'est un commentaire" ] }, { "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": "markdown", "metadata": {}, "source": [ "Signifie qu'il n'y a pas de ligne manquante." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "data = raw_data\n", "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": "markdown", "metadata": {}, "source": [ "On met la date au bon format pour que pandas puisse la lire." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" ] }, { "cell_type": "markdown", "metadata": { "hideOutput": true }, "source": [ "Tri par ordre chornologique" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", " delta = p2.to_timestamp() - p1.end_time\n", " if delta > pd.Timedelta('1s'):\n", " print(p1, p2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Juste pour tester qu'il ne nous manque pas une semaine, on regarde si les dates sont bien espacées de 7 jours" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm4HUWZ/vudc+692feFLEACiUDCEiCEQARZggRkZFGcuIEOGkUc3B1wfiozmhFkFEWFUQQBRRYBBWWRLShLSAgQspGQC9kTsu/JXU/9/uiu7qrqqurqs9wlt97nuc/tU11dVd1dXV99OzHG4OHh4eHhISLX3gPw8PDw8Oh48MTBw8PDwyMBTxw8PDw8PBLwxMHDw8PDIwFPHDw8PDw8EvDEwcPDw8MjAU8cPDw8PDwS8MTBw8PDwyMBTxw8PDw8PBIotPcASsWgQYPYqFGj2nsYHh4eHp0Kr7322hbG2OC0ep2WOIwaNQrz5s1r72F4eHh4dCoQ0SqXel6s5OHh4eGRgCcOHh4eHh4JeOLg4eHh4ZGAJw4eHh4eHgl44uDh4eHhkYAnDh4eHh4eCXji4OHh4eGRgCcOHp0Ci9btxJtrdrT3MDw8ugw6rROcR9fCBb94EQCw8voPtfNIPDy6Bjzn4OHh4eGRgCcOHh4eHh4JeOLg4eHh4ZFAKnEgom5ENJeI3iSixUT0X2H5dUS0jojmh3/nC9dcS0T1RLSMiM4Vyk8kooXhuZuJiMLyOiK6PyyfQ0SjKn+rHh4eHh6ucOEcGgGcxRg7DsAEANOIaHJ47ibG2ITw73EAIKJxAKYDGA9gGoBbiCgf1r8VwAwAY8O/aWH5FQC2M8bGALgJwA3l35qHh4eHR6lIJQ4swJ7wZ034xyyXXAjgPsZYI2NsBYB6AJOIaBiAPoyx2YwxBuBuABcJ19wVHj8I4GzOVXh4eHh4tD2cdA5ElCei+QA2AXiaMTYnPPVlIlpARHcQUf+wbASANcLla8OyEeGxWi5dwxhrAbATwEDNOGYQ0Twimrd582anG/Tw8PDwyA4n4sAYa2WMTQAwEgEXcDQCEdHhCERNGwD8JKyu2/EzS7ntGnUcv2GMTWSMTRw8ODWRkUcXxxMLN6CxpbW9h+Hh0SmRyVqJMbYDwPMApjHGNoZEowjgNgCTwmprARwsXDYSwPqwfKSmXLqGiAoA+gLYlulOPDwEvPzOFlx5z+u44Yll7T0UD49OCRdrpcFE1C887g5gKoCloQ6B42IAi8LjRwFMDy2QRiNQPM9ljG0AsJuIJof6hMsAPCJcc3l4/FEAz4V6CQ+PkrBzXzMAYN2Ofe08Eg+PzgmX8BnDANwVWhzlADzAGPsbEf2eiCYgEP+sBPAFAGCMLSaiBwAsAdAC4CrGGOftrwRwJ4DuAJ4I/wDgdgC/J6J6BBzD9Arcm0cXhjdn8PAoD6nEgTG2AMDxmvJPW66ZCWCmpnwegKM15Q0ALk0bi4dHVnj+08OjNHgPaY9Oh0/fPgePLdiQUsuzDh4e5cATB49OhxeWb8FVf3y9vYfh4XFAwxMHjwMaXqrk4VEaPHHwOCDhFdIeHuXBEwePAxrVUkjvbmjGVfe8jq17GqvTgYdHO8MTB48DEtVmHO5/dQ0eW7gBv5r1TpV78vBoH3ji4HGAw2sdPDxKgScOHgckfFBfD4/y4ImDxwEN7wTn4VEaPHHwOCDh+QYPj/LgiYPHAQ3POHh4lAZPHDogikWGF5dvgQ9M6+Hh0V7wxKED4o6XVuBTt8/B00s2tvdQPDw8uig8ceiAWLU1yEHw3q6Gdh5J54U3VvLwKA+eOHgc0PCiOQ+P0uCJg0enguti7zmHzoEFa3fgH29vbu9heGjgiYPHAQ3PN3RsfOZ3r+LyO+a29zA8NPDEwaNTwVVKRN7ToVNg296m9h6ChwGpxIGIuhHRXCJ6k4gWE9F/heUDiOhpIloe/u8vXHMtEdUT0TIiOlcoP5GIFobnbqYwxgER1RHR/WH5HCIaVflb9fDonPjY/83GSTOfae9heHQxuHAOjQDOYowdB2ACgGlENBnANQCeZYyNBfBs+BtENA7AdADjAUwDcAsR5cO2bgUwA8DY8G9aWH4FgO2MsTEAbgJwQwXurWKYtXQT3t64u72H4YHsYqIDQR89d+U2bN7tQ4N7tC1SiQMLsCf8WRP+MQAXArgrLL8LwEXh8YUA7mOMNTLGVgCoBzCJiIYB6MMYm80CreLdyjW8rQcBnE0dKHLaZ+98FR+86Z/tPQyPLOgws8fDo3PCSedARHkimg9gE4CnGWNzAAxljG0AgPD/kLD6CABrhMvXhmUjwmO1XLqGMdYCYCeAgaXckIeHiAOAcfDwaBc4EQfGWCtjbAKAkQi4gKMt1XV7NmYpt10jN0w0g4jmEdG8zZu9+VtXhLMpa5XH4eFxoCOTtRJjbAeA5xHoCjaGoiKE/zeF1dYCOFi4bCSA9WH5SE25dA0RFQD0BbBN0/9vGGMTGWMTBw8enGXoHu2AbXubsLexpV3H0Fmc4F5+ZwvW7djf3sPw8IjgYq00mIj6hcfdAUwFsBTAowAuD6tdDuCR8PhRANNDC6TRCBTPc0PR024imhzqEy5TruFtfRTAc6wDftXXPboYza3F9h5Gp8EJP3gaZ/3keae6P/jbEtz50orUeq6Toq1UVpXq5hO3zcEZN86qTGMeHhVAwaHOMAB3hRZHOQAPMMb+RkSzATxARFcAWA3gUgBgjC0mogcALAHQAuAqxlhr2NaVAO4E0B3AE+EfANwO4PdEVI+AY5heiZurNO58eSUmjR6A848ZVnIb63bsR46AYX27V3BkHRcbd7lZ2dz+YkAYPjNldDWHU3GUu4V5+Z0tmDRqAACgubXD7Yc8ujBSiQNjbAGA4zXlWwGcbbhmJoCZmvJ5ABL6CsZYA0Li0tHRWizvA55y/XMAgJXXf6gSw+lycH3+nUHn8OrKbfjEbXPw5TPHtPdQPDwS8B7SGeH3du2Lxet3ZqpfbeFkOWKlLaHvQv2mPSk1PTzaHp44ZEQHVIV0KdQV8umV0HaB98qZDnyMzG85/HfVAeGJg0enQreabFO2Yy+8AXXIsi7uamjG1+6fj537m6s0pvaBpw0dD544ZISfxO0LVyuktgq8Vw6Hwq+1qVFWbd0r/b79hRX48xvrIgX+gQL/WXU8eOKQER17J3rgIytxrjYxL0usFLeiPb91TyM+cOPzSn9B3XzHiS5TERT9rqvDwROHjMi1wUf51wXr0yt1URxIyX74XDLd0q6GpAMh5zJyneD+ssATh44HTxwyoi2cq3bsO7DkyZWE6xLSVmtNZcRK7oPldTsD8csCTxs6HjxxyAhvVdG+OJAef2yt5A7OOXSgoMUVwYH0Xg8UeOLg0angusvuDGsnlWCtxDcnbSHebEt4XV7HgycOGeF3OO2LrB7qtvf1k6eWYex/Ph79LhYZHn1zPYoZ+qiMn4M7OHHMH2BfbpmBBzyqgANsilUffofTvvjEba841XPZV//iuXopntH989bg6nvfwB/mrCpxdEk0trTixr8vxf6mVqm8ubWIz/zuVQAw+izoRJiRWKlTBAhxhxfXdjx44pARfg63L3QWPDZkIeY8FWeWlJxp0p17XlmNX816B7c8Xy+VNzTHxOLNNTuc+ztQFdKec+h48MQhIzxxKB1b9jTi588sr9ou8b2dDVi0Loi9NG/V9pLbyaYDsJ9vbAlCvDe1yKHeS1Uo8/4ONJ2DZ8g7Hjxx8GgzfOtPb+KmZ94ua+G24Yz/nYULfvEiAOCnT78NoG2I+f/+fRleX62/p4hzUdbyUvwU9je1RpyD93PwqDY8ccgIP4WTaGopYuZjS1Lj/ewL5e4tVcpb0NDc9omYiIBfzqrHJbe8rD3PDDoCF52Byl0c9b0no2RTuQOMOvjvquPBE4eM8IqzJP4yfx1ue2EFbvz7Uqf61VbqixZNWXqq5nKrSoFcnoFurt07d03Y3oFFHF55dyueXLShvYfhIcAlE5yHAE8akuCLcRpHUOp6lpUgq/L9aiJtaKxKYqADjHHAl+55HYBPgtWR4DmHrPDUIYFqM1NZfRv2NQkWTVUe23NLN1rPm0xPy31mB5xC2qPDIZU4ENHBRDSLiN4iosVE9JWw/DoiWkdE88O/84VrriWieiJaRkTnCuUnEtHC8NzNFPLGRFRHRPeH5XOIaFTlb7Uy8H4OZjivVxkfYUtm4tCaXqlCWLl1n/U81xFUei0vl3NYs20ftu5xN9n16Hpw4RxaAHyDMXYUgMkAriKiceG5mxhjE8K/xwEgPDcdwHgA0wDcQkQ8fdetAGYAGBv+TQvLrwCwnTE2BsBNAG4o/9aqA69yKB1RuIiM12XnHGLi0N7E/BfPBf4N6lpe7qjK1Tmc9uNZOPl/ni1zFB4HMlKJA2NsA2Ps9fB4N4C3AIywXHIhgPsYY42MsRUA6gFMIqJhAPowxmazQBB7N4CLhGvuCo8fBHA2dVCNm6cNpSMKF5GVc8ho3bSnMZujnIqqEBRSxUrl9VEJsZKJI1u9dR/OuHEWNu5qKLsPj86LTDqHUNxzPIA5YdGXiWgBEd1BRP3DshEA1giXrQ3LRoTHarl0DWOsBcBOAAOzjM2j/eC6mJa6nrUU0xXMjy2ILV1EnUNH4fQqvdOppkL6rtkrsXLrPjw63+cVaW9s2t2ATe1EpJ2JAxH1AvAQgK8yxnYhEBEdDmACgA0AfsKrai5nlnLbNeoYZhDRPCKat3nzZtehVxQdZbGpJv7yxjos37i7hCurs2K5PPKr/vh6dNyWOgdXqDt90z09Mn+dU3vV5KtbQj1JId8hmfcuhUkzn8WkdhL/OREHIqpBQBjuYYw9DACMsY2MsVbGWBHAbQAmhdXXAjhYuHwkgPVh+UhNuXQNERUA9AWwTR0HY+w3jLGJjLGJgwcPdrvDCqO9Zdhtga/ePx/n3PTPEq50ezZZn2FWgvzcW5uyXRCimguua9tfuW++U71qWis1h+KmWcs2Y+l7u6rWj0fHhou1EgG4HcBbjLGfCuXDhGoXA1gUHj8KYHpogTQageJ5LmNsA4DdRDQ5bPMyAI8I11weHn8UwHOsg3qbteWo2qKvu2evxHWPLq5+Rygtf0EpuH9eLNXsKJMooZDuwKasraGO559vb8a0n71QlT72NLZgyvXPVaVtj8rAxQluCoBPA1hIRHxb8x0AHyeiCQi+v5UAvgAAjLHFRPQAgCUILJ2uYoxxPv9KAHcC6A7gifAPCIjP74moHgHHML2828qOZ5ZsxJQxg9C9Nm+t15aLTVvQx+89EhCG6z48HgAy5TJIolpipeo+B8ZY23scd2Di0Oyg4ykXi9ftxLod+6vej0fpSCUOjLEXof/qH9eU8WtmApipKZ8H4GhNeQOAS9PGUi0s37gbn7t7Hi6aMBw/m368vXKFFuzWIkM+RavYHmGMm1qzLwyuj6SU5DalXSBc6jA4xrKLlLJOg0q/y+rqHBzFg4xh8fpdOHpE38x9dFBjRA8B3kMa8Ye7eH26fLVS33izwyLcHpEqGzWhJ1Zs2YtXVyZUQAlU63uv9lNoi6esvstyuaFqTg0X6zAAeGDeGlzwixdTvcR18LSh48MTBwC1heAx7HZIJFOpj9Jlh94eWpco6qfw8Z75v8/j0v+b3faDCVHOcyjlUpf+snMaHUX7kQ5XqdKy9/YAAN7dvDdzH542dHx44oB4VyfumMZ/70lc+KuXqtZns0NwuB51dv1HNdAcmTFWb2pkXSjbQucAZBN1ZF3r1erl0opqPhPXx8CnSCkcrqmPzkRE0/DJ376Cr93vZn3WEeGJA/SLw96mVm36xkpN3maLXLdPt0AVVGiH0Js8VEXecYXYuqcRd7280qkuf75Zn2C11wvefDUXpqRYqTxU85m4EgeuFC9Nn6Lv5EBKF/pS/Vb8+Q03v5WOCE8cIEbOTEepc7eppSh58broHKq9KH73L4uM51x3pl+9fz6Wb9rjVDd6vh1sASjlOasLaNqCWsqiZ7skrblZSzdhzTZ7UMByQRFxqBznkDWOlkf14PM5IJ6QOvNA9QMrdcG++dnl+OWsOMm8TefAlP/Vwu9fWZXsW+i0obk1Nbvbjn328yJKzueQdt7yUlzeFyeE1RQrqSbC5XIpadd/9s5X0a0mh6U/OC9z2y5Z6oDSY2UFfejh04V2HHjOAfGE1K0Np/14lvS71Km7XrHptnIOYSft/Z18/u55HSJyZ7Xl0GrzpfRmWuz4nDpqWJ8SWi0P1UybyhjDqq2BItrkG/PXN9fj4795RXvORIg9ceg48JwDYusMF8eiUhaqJet34WFF9mjLVhZzDm3/oYg5j19YviXTtf982y3eVdb7erl+q709S3PtvdSM6Ncda7fv16QJ7XjYtLsB/XvUOslX7527Bo8vfA+AWWT27/e+AUDvZGj6jrxUqePAcw6o/m7lR0+8lSizcQ78w2mPTVQ5BGntdrvHaylSpZ37mvHthxaUNqASUco4dTvhppZi6jPRITKQsNbJ3KwVDc2tmDTzWVzz0EJjnT2NLRh1zWO4e/ZKKeZS2pzRcTAmInAg6hw6a3wqTxwAtPI8vw5Po5SPsqE5GSW0qcUiJ1f+twXakiBl6ePnzy6vemfVEivVC4r6RB+WTnieBbtCOnm2qaWIbXubLFeZwXVgTy1+z0iUNu8OMsfd8eIK9KyLhQ6m9ZxHADjqe09G13KYOIcDyZSV44o757X3EEqCJw4Qd2rVMR3dryEOTkrUNvxQSt2xZVEyR6asFb4tsbm0kCT668sfUNpzUPuw9ekSvkL3DL90z2s44QdPO1ybvjFJu66nEIPMpHOoFXxlXqqXRZSm+XYgcg5pRh0dFZ44QDBldVhXSllIdB+7y66wLT8T/gyq2WcppDerL57qn+FyPx2NW+LOmFnH9YxjqPK/L34vUSY+tTSrLSJCTiDCrYaBin46qujWRAMOQNrQab3BPXGA3ZRVxf88vjRz+7p2XfQcbclhFyOxUvU7zdJDViMBF9FgOePhSPg5pCwBexpbZPGiTaxUIufgim17y9vJqnPEdOdiLZUjMM1/XfmdL63A4vU7swyxY6GTUgdPHCCYslapfZ2ow2phE5myth11YG3BObTBR1LNUNYiEq9G0+2uhngR/t4jizHtZ24JlOKQ2aWJf55YuMFy1g7GmMUsNz4jEkOXUBhJzsGdOFz31yX40M0vGkbljh37mrBlT2N6xQqjk9IGb8r6qd/OwdZQiVetdSWnIQ42zqE9FNIx5xD8ruYam4XoZX0GCbGSixOcGtrC5Rqk756nKzb+K7fGDpW2Lvgu276BMJ+88p7Xpd+MMfxl/jqcd/QwdKvJt5mJtNiLKi5qD2ulCf8d6GNWXv+hqvWhQ2cNT97lOYcX67fgrQ2BqVm1XqIuRpL1E4g4h6oMR4uYWJk71X242Z5YabGV0iC2pxJil4WwlPGo76aS+SCcxEoZ+nr5na342v1vYuZjSZPqLO1J4jvhfk0iNfEeE2IlAxGwPZc7XlyRPsgOCK+QPgBQrTh3uiB2douR6iukX1sl52dQv1Xd8Mo1K20bsZL8uxT5vcs4b09ZqOauSM9/kTYelw2EC3aH4q1Nuxuk9k19pt0/EUl1jGIlYZDqfDdxzjbO4b//tsQ+MI+KwhMHAdUyZdV9POXqHJ5ctAGjrnkMH7PkWVi3Y79x1/KRW+XrGltaU8c1+51sHtOVgGsmNw5Vv9PiIqZwIIwqNil2++rc+div7fkvbBxNtDmotMmvligkKVGlvgIr52C0VnLjMDyqD08cBFRrZ6tTSNvmvMvn8MU/BHLluZYMbVOuf85ZCbpoXWANsrsxSHikexaVIp7VFJepCukWl+i3UeC99PZbWouYtczNZLRUuNGz0h+ieOXvXlqZaM9FvCoppB36VImpiUNwVVx7VB+pxIGIDiaiWUT0FhEtJqKvhOUDiOhpIloe/u8vXHMtEdUT0TIiOlcoP5GIFobnbqZwFhJRHRHdH5bPIaJRlb/VdFRL56DnHCw7xwp6K2/Y2eBU799CL84vhURnX1PSca9ctIVaTtXvmPJm6BYnl+f998Ub8dnfvZooF99xVm4neY6LFW1zJLWLBHTz8NWV21AsMmzZ01RyuyaITd3y/DvyOcfYSiYfCo/qw4VzaAHwDcbYUQAmA7iKiMYBuAbAs4yxsQCeDX8jPDcdwHgA0wDcQkTcnfJWADMAjA3/poXlVwDYzhgbA+AmADdU4N4yo1qLl27H7aQAbIfwbO/tMhMTLe3MQFDjqu73NX54evJ68TnxlK8cpnzIj8yPAyFmWX927NeHpxCfgovFja2G29xwR+L+mMglAD975m1M/ek/4jKHPuessAdDTINp0VefXWelDaXk1e5oSCUOjLENjLHXw+PdAN4CMALAhQDuCqvdBeCi8PhCAPcxxhoZYysA1AOYRETDAPRhjM1mwbbhbuUa3taDAM6mam3jLSjFgcoFmTmHqE51xlMqyn0jnEhmua+6QraXohIH00ItjoEf8vt7fdV2Y/sm0Zo4XbdnyHGhH5uLKWv2dk1jd/WsjttRrjFMDBuRdNU5JIlFB/soDHjo9c6bAY4j05cXinuOBzAHwFDG2AYgICAAhoTVRgBYI1y2NiwbER6r5dI1jLEWADsBDNT0P4OI5hHRvM2b3cJDZ0G1FNI6OCmk22YonR7is1QDz5nESmLgOHXBselxXAjkdY8uTq3jElurUutggnEQjtU5z8BSWeh3t+yVfpuq24iDUaykMHppxKKjolETT62zwZk4EFEvAA8B+CpjzBaDVjdXmKXcdo1cwNhvGGMTGWMTBw8enDbkzGjLlM1OSscO9h1UTCFdpboAItk5h0khLekIeJnD/Rm9h4XjUiOjclRLIa0XC2ZupiJwVkgrr8/J+qwDoJqJltoKTsSBiGoQEIZ7GGMPh8UbQ1ERwv+cz1wL4GDh8pEA1oflIzXl0jVEVADQF0DphuKloky5yfcfWYRR1zymaVanc3CRS3eOD8EVpTzecglks2ExEYcSc2rpnbmE5+BmwTZYOcfIzyVZqX+PmtTrs0BPL7K9KJ2TZxqMHtIp1kqdhXOwJfPqLHCxViIAtwN4izH2U+HUowAuD48vB/CIUD49tEAajUDxPDcUPe0moslhm5cp1/C2PgrgOVZF4eI5P/0H7p69MlFe7ibqrtmrnOu2F+dQzmPVm7d2bBx/cD9teckqLcNlYnONZS4MfLese1VR2PMM7dlySSQztAHzVun3Zaapk8+7PUvRZ0Fc9D9yQrxnVOenSiwqwTlws+1q4kCwsnLhHKYA+DSAs4hofvh3PoDrAZxDRMsBnBP+BmNsMYAHACwB8CSAqxhjfCt1JYDfIlBSvwPgibD8dgADiagewNcRWj5VA40trVi+aQ++90hSLlwtsZJWZlbFybO7oRl/eUOvEBt97eO48g+vVa1vG1wS0icUkBm5pzFDekm/h/TplnpN5OcgvKlP3z4HY//z8URdl6B0uuROWWC759jgqzJiJdKUrxJiQDm167hFaBJEfCKhqC3E16tSwGo4xV3wixcxz6JXqgQOBP+M1MB7jLEXYf4mzjZcMxPATE35PABHa8obAFyaNpZKgMuDVasWwL6bvOSEEXjYYoGwUTABLRaZNtieCNPcESd/qQTkqj++Yc3n/MSiZDx/F2QZjm7skbWSLdFNsYh8Lk4kIzZz2thB2rzWYp1LThiBHz+5zDqOZAPJIlP+bBeOw0XeXKolUiWWHPGZVMIm0HUhFOuJa3yNkLQjTefQbDBNzgpX/59ScSB4dnc5D+k4d0PynG09P2ZE0t6+tcii9k7+n2ejcnUCa01ZDZ+5yI6WuvlYsl5vL1Aut9IWu6Ej/t+TuDpMTK/i05MPTb3eNdkPOdRJu04ql8RK5XEOLgtXpd6Eej+2dsslJOJ6Kc5ziTik5X3QDHDb3iZMmvlMppwP/26YY5XCAUAbuh5xsJkJ2tjjHBHOOnKIRCTO//kLGKMRPajB3vRiJX0/olil1PllcvzSpSvNgpoMadnKsdF/9M310bFkdumwOiVClSh9zbh7HmY+tkTxaE5tVhhDep1GF87B8nY/f3fgra4bF1f+ZuLilL72hCFSAM0zLWHSuW46itLGx0AclKZUMaNOlv/C8s3YtLsRv/7Hu07jaAt0FsW5DV02n4Pu1dk+/BwFi7w4wZdt3K2tqyrNdIuaSzz7UjfqrQbbfpdFywa9KE5fVxzBEws3oG+PmpK015IIxFRH6E21JlIXxqeWBJ6rpxw+0FjHBtP9isUNZXIOHLpx8YW0FC6Qb3427jInvLHmGTGJQh2HwoTpJ+kcBIW2uvirYVx0ffHva9W2bLqSauJA0Dl0Wc5Btx5YiUOOQEROC7ZqW59FId0icQ76OmtSPgKTRUe507VWwzmo96bzK7jyntfxidvmRHWr+dkkGAeHzpZuCIi87v3PXbENF/3qpcg00cVD2uR4l3VcOhTy2a2VVDQL70i9m1IWNddLxLbFR2TTOahzXSfL5+/8zTU7OsyOvVwuHUC7ZK0T0fWIQ/hZ6T6CNLFSjtw+niTnoBmH0swnbnsFv33hXUUhrW//T6+t1Z8IYQ4ZYR972vnutXnreQB4ftlmY1uRGWaJK6M5HWV8rIqVTF2J5fe9utpY9xt/mo/5a3Zg9ba91jFkZYpcnoBNrJSpr7CdxxZuwOV3zMV9r8YBDNTcDPa11bTpKE+sJHKk6uKf5vcAyNyiSaTa1tgriO6A0ua8Jw7tBO0kszyNHAWT0OUdN2cIE83x8jtb8cPH3lI4hyR27GvCzSlJd0w21mmbqrR76xc6YIl4ffUO6TdfnCu1f5Nt8tPrmzLBMcZw09NvR+Xis+hWk5fqilBFcS56j2NHpgcLLBWFHBcrJc+ddeSQZKGCfyhWbFmUzOWKlcR64vdn0zmofer6El95R+Ec9iSIQ/Y2RKLXHjGluhxxsMUtsnEOfIflwjkkqyTbNQU7kya3pi+XUNpGziFlyU67M5f5mTVQXhqkPAMO+/P+PWql3/xRLF6/S8pkJ35s3QohcdDcn+rQ5mKtNFTwrThpVH9Nbbn/H1zfFwp3AAAgAElEQVSUsO4GoBcP2sRKurwhprqlwNiO48Il3rO4f6oR5kxy7qZzDiRxDh2DOIwb1kf6Xcqocs4cXXXQ9YgD/59R55AnCjgHhz5cCMjTS/QhfUW2OKvSPBWpnEN5YicAqKsx72zLNql3UH6LAfWAWExxwS9elMrfFowJ+nQvSHVFqGapbuKguJZpwRZx5EG9teXfeXhhoqxgUUiXsrt0SYYUt28oN9TvXiOLIcXHK8rkRYW0eg8JziElj7nJGKOt0a0mL20MSjIgED729lBwdz3iYHnItg85l4Mz56DufrIs6KLIVNeVrqyhuRXXP7EU+5pakifFa9P6TiMeKdcDQG3erJf464LARDWbM1187OJYlFCQG67536diEdOJhwYfsa6q6tBmmj879zdrQ20XDLJKsZVuhTwe+MIpiTo6a7gbPnKMtr3XVm3DK+9m9/rNEiDOxHma3qc670Vx556GeK6KYiVVJKo2rRcrma2d2gvNrUWruMwFeUmsVIlRZUOXM2W1PWPVgUpELuQcXBYotUqWHbPMObjNiD+8sgr/9493UJMn7LeIndImWFp/WXYvurZikV6GdoSqJl2OzePXKfFOWKWc3VlzK8NtL7yLGacfLt2dzvxX7JPDRc/8sYkjMaJfd+05NSe4K1ysalpaiyjkcxadgxvREL+dPY1xzgtxU5auc7DrCk3ve+IPn9aWVwtNrQw9auOBlRJEM+c5h7ZFqaEJiAiFHGl3ov/1VzlOUznKo13CjspV9MVNJ5tai/jtiyuMbafqHFLFTvbzYh+2uqUalLiYiKpwsV7hr7Tc9J5cjyS2o4pWdCDKHgiwUmuFSxyoeWHyo6xiJXVBE3+aghOqmy91zmbVORSLDLsbmhOh3KuNltYiavKV2/m3B+fQ5YiDjQQ8t9ScEStPhLqanHZS8yTtHAnOIcN3f9GvXoqOXedD9AGVubinmc656VvS62SZ52Jdl4VeVFrnyDVlZ1Anax4FlSvgG2CxmRpjxNLsokfGxKis+sHqLKVsRO+MI4akKvpj/xS5nUmjBqC2YOYokiKhuEQSF0qK6jSdg3l8ALBVmcM3/H0pjrnuKf0AqwhVrFTK4i4+7/YI39/liIP4knS5F0zIUeAE5hKnXd3dlJokR6/UTbZls8AC4t1h2vR6/w2zUsaTPkH5x2n3tM0iVorrmjgHsVRcZAu5XCaHtKzy6rOOkE1H85rQFq4cgUueCEBYqA1D7du9Bt+7YJxTW4CNeAl9UvK+AODso4aAYHmfFpGQeErUyyS4DaVJ3bwSiz78y5ewuyEWWT1ZYpDJctHcylCTz+HyUw4FUNriLpv+Vmpk7uh6xKHE64gIdTV5p6BqCeKg+f6Oc7CFzzqhTB8pTx5frq20y+VX3GWOCxS1k6HPbz24IDpWY1alIZ8jJ86Bv6+R/fXyfBHifanvJxft6mPsbtAbCZSicwCEsOeG8/kcoX/PpD+KCS7v1ES3WhkDkbtYSbbSZjhsUE/cfvlEnDNuqPEadc7erhGbqtf8x0PxnOnTzf1ZVBJNLQHnMCKcU6Us7uK9e51DGyDtGZtM+0TOId3k02Ec6VUyBxIz9bt2+37r+XLb18E6mUsch0msJO3UhfK8QUdkgs0ggbdsu62IODCGHqE3+REH9dLWVbkdF86BIZ0LNVlHmdssfVIUiyx0DE3n6IB4sdu5vxl/W7AB727Zi7OPGopuNXn8z8WBFVba6xKDMnKoG4D1O+KottxMua3RUiyitkDSnMiKFJenqqPLEYc07DMo6HJEqCvkUGTpjjY2RVy5eKlen2cAKJ0rckU2ayXbudJG6qSQFtbOwb3r0Oqgp0gTy4lnbXXE9f2Ig3rjb//+fnx16vtS+1evdYHpVRRylCAgVuOADI9Ubae1yINRmsao5xwuu2Nuou7Uo4aEbWb/dmz30F6cQ3MrU8Rl2dsQjQW8h3QbIG1hMjnR5HMUpUNME1UkzPEsZp1ZsHrrPnzjT29a+i2fo7Fen6WuZU12+VDe1tj4u4QlERfG8cP7OHEOmYieQakKiJxDsGgePaIvavI5rZf0U4tjWTghuaCbEIuV9GMu5KkiITFc+mwtFq3BKNVHz5/zm2t2JOqa4m65bCTU9yeGruhV1z6cQ3MoVqI0OaAFouOm1zm0AdJt/QOoijqR9U9rw2X3I076N1ZvtzcYYrdgG84xsGccLqJcP4Y0uK6hLa1Fa18u7fxtwYZEmZEoG4pddQ4R52CtmtQnqNVzgiWRqIgeP1zWL23Z0yg54REl43rpdFuSvsMw1iw5N+KxptXSz/uWIkuEsY/HlyyzEWGuzC+Fc/j+o7IpubjjLjjmt640motF1BRIsGAr99vrgJwDEd1BRJuIaJFQdh0RrVNySvNz1xJRPREtI6JzhfITiWhheO5mCr8eIqojovvD8jlENKqytygj3ZafaevV5nPRi07foau7H/s4ZllMaE3X6Pq78+WVJV3vCtcJ2tBStO50Sv1Q1PzQOhABj199Gv78pVOdAyWy6L9dGCZVhlmpzDmHuJ664072o3IO1zyUDJ3hspAXwtDyrsgyJ9SqrSxIh6veX/2m3RHH9sUPHI5bPnlCal/xt2Uf34CecuyshuZWbN5tM8FuJ+LQylCbz0W9l7vz76icw50ApmnKb2KMTQj/HgcAIhoHYDqA8eE1txAR9wK6FcAMAGPDP97mFQC2M8bGALgJwA0l3osTXBamGXfPS4gjDh/SK/qAi4xZTVpXKzHol76XTNsp23nr2+nhECI7k6jHUD68bzfDmdL6spo3wlGUofy+53Mn49zxBxnGFTdIAMYN74PjD+lv3NWqyBJMkSm8g8hhiqasUqY5pS01TEvAlcp1XjToltLETz1q8xVfDiPJiErkWoP7bxYm8MotezH1p//Ej59cCiBwAuTBGG3PmUfTTTNl3ba3SfJluHfu6kRb3ACjvdDSWkRrMTBl5bGwdCLRdzfvcQ7L3SE5B8bYPwG4Bm25EMB9jLFGxtgKAPUAJhHRMAB9GGOzWXCXdwO4SLjmrvD4QQBnU5atT0a4iJWe0gTF616bjz6SIkM0+XX4+gOxXqChuRVrttknq+6jGdCzNhGQzcw5WJsX6hn0KY6sd5YJauccsuPwwelcgwrX5Ex8QLYMabpmWooMf7lqSvRbEitZlmidVZQ65bVhWiSuRX9j3YR5Gl9mfgi6fhLEK9k9gIBzKORykoXf5nCxe3VlICrloe4B+5yI66hipeRFv5xVHx3v2JcUtUpjbwfGgRu19KjNRzoPNYQ3AJz1k3/gtBTfIo72iBlVjs7hy0S0IBQ7cY3bCABrhDprw7IR4bFaLl3DGGsBsBPAQLQTxHfAzeuAIBR1TlAu1W/e49QeFyF85tRRcj/Cse6jKeTIzWOXsUSMfmNdQ7ndhDOGe9x+lqJzqN5EFxfZHLn1xcdqy5Oh00s0NhcxfnjfaAHIiU5wlLw27k8ZMyjBOZh22fEuXj/O7jX5TE6XumYunDBc6VOvcygWGQp52VxYtfwKQoMExzb9D5+D6gZbd0XP2ljJrCbV6Qjg8c261+ajKMGmcbpmjMvq41MJlEocbgVwOIAJADYA+ElYrpuVzFJuuyYBIppBRPOIaN7mzW4LYqLhVM4hrjCoVyzfVHUOupSZOvCPXA2YJlu9JAdVk885W27w2DdpMN27miDHeL1TraCe1c0hgx6Aw0a/JD8HoV6O3AisWx2WGBdXGvMuxfAZks5BuRvde1X9HExjMu3iOXSxnMp9F/HI5MotRRbEGxMWLv6c+DVEJIjbzJ2R8G2JmK+xbOrVTSAOKZGI20PjwHOudK/Jo2coGt7bWF7aUJcYWJVGScSBMbaRMdbKGCsCuA3ApPDUWgAHC1VHAlgflo/UlEvXEFEBQF8YxFiMsd8wxiYyxiYOHjy4lKGn6xyE02JaTCKS5KI1jkltxB2UCbpdYm3oU+Havhv0tV0csMhxFw4ArJi2IDmxRHL/Tj3LEEOsTz5sgGtX1joS5xDqnXhRZNeu6hzUHbfyu7G5iN7dZJNL0zMy7eI58rmMpqwOM4jXSN4HQyGfkyMJh3X4oi5a+dnmMycgqpjrj3OSOgVxY7YnZdFtD7ESX8i71eSdCKMLzrnpn1i3o211KSURh1CHwHExAG7J9CiA6aEF0mgEiue5jLENAHYT0eRQn3AZgEeEay4Pjz8K4DlWRblDWstiYL1uNXlJKUzCJK9x3W2HzdkWYN2YArGSuuN06jJTP4CbWKkmZw6wpqLImD22klszzhDbE++EhORMVmLltEAy6T8AjB7UM2w7KKuJsrTZdQ5qf/uaWjGwVx1+9YkTojKTyiHmHOIKDwo5xfO5ZM+2ey+ydMKr45qAQExUv2kPHl/4XhTPSL03Ahm5AhEuBIRD1Ik0pYSzKTWuWTng4rN8TtxQZmtDtwQ+rjHvriZSPUSI6F4AZwAYRERrAXwfwBlENAHBfFkJ4AsAwBhbTEQPAFgCoAXAVYwx/vauRGD51B3AE+EfANwO4PdEVI+AY5heiRsrFV/8w2vRcbdCHk9//QNYvTWwPopNFRl6O3pe8g9CpSXiu9cpm2rySc6h3EXVKKpw+H7yOXI2QS0yhg07zbucUsRKrt+4LFaKPzJbl266Hfn/z6dPwNSjhkp18sJCIPotJOTFys/9zYFo5AjBACGLY943BcdIHedgbSmLpZZSVdQ1rNq6D0ePSMYLkzkHG3EI/psUrw9deSo+cuvLUt3guH1MVXVoaG7F3saW6D6D7JHBuayxkXRzcubjb+Hzpx9W7jCdkUocGGMf1xTfbqk/E8BMTfk8AIlkuYyxBgCXpo2jUkh7RYvXx2an3WpyGNGve6QvEHc3piQuANCvR0w4IuJgSHwP6L/PmkIOjYqcsezAeYa7d0llWciRcx6GIgO+9acFxvMud6Heaik7wECs5NBpJlPWAKccPjCRkrRGSOEpjjct1Sh3khPXOVOu8DSFdDBHVf2FSUTl+i6Y9J9D1KNFitXEe3NzHuW6CVNsM7EvyeggZe62Je3411/Pxptrd+KsI4NQIHnB5ySrtVF7BNpT0QU9pN0fejdFuSfuAmwv78PHxdYefHFSTRVlPweNziGvEStp+soyh+wLih35fJJz2Llfb0bIGEOzhZJksSBygdyeaK1EAudQnpgrzpnBwl6SatpIvgx5UVJzgKi3P7RPt2i81jGw2POaN6FaAGl1DiaOUTMWHXS2+FeecTiuPnts9Hv+6kDHoO54awvuzqO96gpak08gaWjAkSYSbUu+4s21OwHEeWHkqArZFvsOQBu6IHHIUFclDiSwx7aX1yqZ9vHFxDyOu2evSrRRyLkppLPATBzSry3kkj4DP/jbEm3dIoPVt+Oh19cJY2JaS4wE51CSWCm2VrLK3UtJ/aoZT6QPUOomiINhFmZayMJOVOeqPCV5LNOinCODqFAp+uIfXgcgz+upRw2VQnWsD8WIant1hbykq7OhT/eCMcS5eE/ifE3jerO4TL2+ejsWrdtprdPQ3Oq80OdzglhJ2Supc07lmDzn0A7I8sy71ciPR/z4bS9PIg7h/6yy0ZpCzskhKFPiHJO1kpNYKUmsREWoSz8cb22IRXd3z16FI7/7ZGr/thGaehOtlaxSpdTe7boLLmIUY6yJi5Ia9qMcDo73w5tQCU9OEz5D191xI/vKYjcHNEkLmHwhj/Gl3ltdBs6hZ23B7Ldg4Byyio2uPONw47lLbnlZCnanorm1iCO/+yRueHKZU1+BzkGvbxE561lLN2HMfz4hEaYOQBu6HnHIwjskxUrxi7ZNdFFRp1NIn3ho/9RFvVaT0D3Lh/zZKaMSZaYuXSJXBjs0twFkGedfNfH5oenJdQcoWSsh5nbKDecRqS64mFA494crTgYA9OtRG1USz3/zg0fg2JF9I3NVGzGzYeKoAYm+1TAuOvcb3Vy9/wunxM9H6XizIaSD2JfaJP+9T/E7qKvJZRKtGJ8N9AShrpAeYkaEa6gYHXhYHNOcVUEScZDPiRtILoZ6PQzA2dRSxLUPm3V2bYUuRxyyUGQ1wmVkxs70zjkcIovI54C42+lRm09dZmvyyYiiOg9TUzuuTnoAcFCf9A8mrxErmeAipklDqTsnkTC7Wys5LFqKUlYkVmOHKpwB5AWstpDDxEMHpBIqG3F45uun45MnHyKMJ/jfpIgjchqxkslUGqTn8l5Yro/rpPYl9RH+v/re+VJ5XSEfL5ApBg22cCcmnUMa06u+26NH9MWUMQOdM++J4J7PqkTBhMCUNTh2Cb7I8dzSjfjLfDcCVE10PeLgWO/EQ5Mx+EXOYdG6ZDA9DtHShC+UrkpCjsCUVa5UrhzSdHmdw2QvaBTkWfvJgqS9vFt/6wVHoZxjCBKn8fCFXTMevhyLEX3V8UqEyrgAmu9yZP8e0XkSdAUq57C3sdXJlDXHzSyZZUAKbPk0TMSqrpBz8nMAdEHR1XPhsfCDhX288d1zrOMSMeHgfiWZwHIRnknPcejAHtJvmwOgbV5myV5YTXQ94uC6+9VUdFWsPbVkY2IhyBHhS6G8UxVX6VDIJ3dR2p1XhoXGmCDGIbVkIUfOhDVtETh6RJ/0RpQmXL/lQwf2jK8RxmIbUkliJWE8qnmpms8BcCNUtt2suCCJVkYqcVixZU+SOGjnckDUsmw4bJGITXNLjEmWdv+BF77pnGyFFvXLGPp0r0F/JZQ3v2ftd5zxvjn4/ZsIi47TN+pbmHjozlW0JbogcSj9wYtOcGlYEJq1RbFmCPjWuUdg2Q+noVuNXqw065tn4PXvnoPHrz5N6yGts5U2jUS30JiG7eLnkFc8pOs3JTO1caR9eCP79bCezwrx4xo3LCY8koe05XqnsN7hB8s5E0kGrvSh4xxcCJVrJFdRIV2/SQ4AGYiVzGbTcRsUL8aOlNdGHEzrWV1NPhKtpHIOZH5PphHqnjUQp5S9Rwm/MaBnbcDFWUeiB+ecTN9LUkfIIFo4quc4/vBKMEbequg0KfpMtTW6HHFwhW4eZ3Hx56IaXjVHgRVJXSFvzHcwelBPDOhZi3HD+4RmmOlipWzWSnoUXJ3ghL5sMfPTRuQk41d+uzrByd6zwvOx9FlkwM+eedt4XhzPb19cwQcUj01w8tq5rxmL1+9KiAZkQhWf++nHjtOOXUVO4hzi46v++LpckZBYLY1OcHBfJN/asAub9zRFv9XrTHqmukIuImx8HCYdV6AgT9fHqA6kOtrWomGziULO0qLbsIETR9P3ogt3Y3IAtFo7CufOO1qfw6Qt0OWIQzkMm45FHKuYKX516tiwrvxBqGKItHEQJcUQ6gf4+dNGGwnVlDGDEmVLN+j1JLqdkPrBqeauqgmliDSC5RQFNeFYYKksVCVlhx35OQA4bWzymQDA3BVb8bNnzOG6geTHrYbpCPpgeOj1wLxXDaNu0jlIcuoMYnDjI9bspG2iGpdFsq6Qw3k/f0EKaa5eZ/IADnQO8qbqCCVPSTwe8xgkp0NFJKMT8zRrQlwfO7IfgGwSABFcn2Iy/U5w+haxku0bEMVKVUxtk4quRxwc54Oumo5FfOprp0t1eKA+1TJFss12GIcuF4E6oYxOTABG9k+GCL/m4WTqScBAHJTfqhOcjTikLf5SuHJDHVWG7Apx3CL3Je7iVLgoJ+3hpuPF771dDYY6MqGK2xXH7rgQGKyMAGDquKHJxEGW3TgDSxDirwiezwDw2SmjE9eqocFNnIPoIW3TA3AYz4icg1DJFDhQG4aDb9QgEysXrN+xH1ff+wYAoH+PWhSLDN/805tYvD72TVC7LDJmlDbY5pNIHFzNZquBrkccHHkH3cvTsYgJxWMU3iCo9EYYVkB23ElfBHS5CF55d6tcEC44Oj8FdXG1LeYuOgdRZg4gEfdJRJrYyEWspOYKdmX5ZA5N8HMAM96ni4FAUsylq6TfxQJKKA8DcXT2AtcNKMRJowY4MyB8k6LGftJZ3Yj46ceOwzEj5SB7fEG74NhhUrlkyppCHGybJtVCKTpm+u/JlkivFM7hodfWRqLCfj1qsG7Hfjz42lopUGciT4cg8srCOYh1TR7jbYEuRxx0H9Wb3/+gJPs1gU9Bl9DDvMq3HlwQlqvDsE/MHCX7uZ3Lu8W+WGDZpKKPEjVWNTG85IQR0bHqz3Hq4QMTH5xqSaJ7BD+4cDyAdHt2fum+pha8ZkhUlIwrZX5evxLSRqqZ4IKxMqPiEoCUB5rjexeMk/tPiJU0zwdmeTQXEzLGLJyDHlPGyIkRTWLJY0YkA/jpxi6OibFkPnSRiOaU996tJodLThgJFVyspD6XIElWSBzCbozzg8wWcQTgj58/Obwfkbjqr9At/KqINwvnsHVvrG+Rw+Mk2+fo37PW6ABodaJth6xvOnQ54qB77H2710SJwKN6moqiE5wJXBFmk1E7iZVyZHUmO/OIwdFu3sXprLFZ/iJ/cqmoCE0XK4m29et37Mc7SprU4X27YUioaBQn/oeUnWRwPvj/7QfNXqBJNtxYFXdpYlMBsvjApLg0Yfxw2dw2zTiAv1OTPJoTDT4WDpfd6z2fm6z0ZbCWMSyUaVFZV4Yh6ePyuH01I6HJuZLPQfV+agrJfA4m/YTJUIOP6eD+AUfDaxSLDA+/vk6bBEfLObC4raAd90VY/EYemb8eL9ZvkdosFlmCOxf9KZJ+Dpq+DZZN7YWuRxxMuyiHa01maSKMsVTE3QAlxzFQEaPkBJHImm37sFUJaXD0iL6hzsG8A3rt/03Ft849AkBSdEBEURhk3W5XXUgp7KeltYhTr38Ov/7nu4lr+L2Lylh13EC8ACx7z2wO+92/LDKec0VSfGBSJCbLRqg6m0TbKucQEE9TlFC+Gw82Dyaxkjv1Yoxh5Za9SlnYjmrKamhDnGNyuVxH3CmrXCYHX/DV5mrzOSmDIh+7DnaFtHA+vNyWHpePp2/3mIOOiYP82wXqJ/Jff10cthE0smVPozbMukkh3UHWfyu6HnEweWA6iH3STFl/+YnjccgAvQ2/yLrLlhdBY5+afKjSVzyhTvvxLJw08xnpfGBDHfw3Oc0M7FUXEQCdjXqUlCSnLnTJBYaLsB5baM5GxZu58e9xYLJNu5PEoZSdUSnfkpiFy3a9bic7sr/8Hvl74hZPqp6Hcw6mqCWcALcWmcI5WAZmABfxXXzLS8oY4/MibDJ+3TnR1JRIfj4mvU0ktVSaEzPT8WZM8zVN56Du+D/269n6yog5mf6CnwBvuiYUAdg8vlWo993QLKeIbTbdUxbOQYMfXpRIgdNm6HrEwcg5KBYemnmTpsi64Njhks27CNNE5EoudfeeI5I+ymTgrljhaptoeWXXJoKX6XQWqo6CEyJVPKWOWcWQ3nWJMtcFUZItl7HVCkKsM+PO1EUsx7vvVpPHUcOSHt5cRGMSK8WcA8MSwaRYJDKZFMkAtu+T82mY+CPbgqs7xQP8AaFhhPB8zD4xQZ2EuI3i4HOROMg4HrP1neja5zIVInGPRoTHkzTtTck/rY7N2o/hpkxrhm3KiVVHCR7/owZW1nk0DV2POBjKXTh6Psnf3rjHUof3I/fEs33xvvhkifLNKgu0zs9BxICeNdHHYnO3tynf+Pqv2w2+X/GT4AuJUQEI/TO85ZMn4nefPUkqc901/fmNOO+DmnXNBQnRj6Gem1NeLBLRrY88wqlJrMQX1ZZWhq/cFwSnu2jC8ITVjwtM/gnRApRQSBtnfepCS1CCGSo3/8fPBUpivpmyiNJRZAzNrUWj5ZzKOcj2/oI4yD7kqC/xvzi2nnWBddqeRn2yKh1MNJEZiGJ8HVfGq8TB7RsQI9u4piauFFKJAxHdQUSbiGiRUDaAiJ4mouXh//7CuWuJqJ6IlhHRuUL5iUS0MDx3M4WkmIjqiOj+sHwOEY2q7C26IbHb0tUJKz365jqp/IVvn4nZ154l1eFz4aRR/TH5sAHSIsB3fkD84dUo8Y10fg4c3/zg+/BvU0ZHE88WqEtnfnvF+0cDAH57+URcfPwIDOol7+4JhNsum4i53zkbP/7IsfjymWMiU0zbnNZxDgN61uLMI4ZIZTaid8NHjomOv/5AnBvZxdw0OR7eHzMqpHOUtE/XQdyJas1VKVgouGHDF5Rcv7xc9Nw980j5uWQyZbWeV3QOFs7BtNRy0RLnGDnUjcSpYwZhWN9uMfHUtCfu+M+48Xkpn4c6HnGsKrcdBzjU389nTh0VHetClfCx9SqBczCJ04pM/q/C7Ofg2K8wKdpaUe3COdwJYJpSdg2AZxljYwE8G/4GEY0DMB3A+PCaW4iIf9W3ApgBYGz4x9u8AsB2xtgYADcBuKHUm3FBKcowDv6i1VSGBw/ogWF9A9m+GqGzsaWYcBgSP4LWVr3cX+fnwHHJCSNRyOccxxz8FyfW+ccEFkQTDu6Hm/51gna32702jyF9uuFjJx2Mb4ZK7TTZvWukS5uIqJIJ48UPk4ElFs17Pz85MudMQ7wI6MVT6jr7pTPGSOdFnQPHDkUs5OJvwqFbhE06B9PdqWaqIp74yml4+mung4gk00qTw6RtkRRTm+osi0zQiTaDdlzeV3JHz8fGlepZdA5pxgJGPYohrpQ1v4hwfyKn1tbx+FKJA2PsnwC2KcUXArgrPL4LwEVC+X2MsUbG2AoA9QAmEdEwAH0YY7NZ8FTuVq7hbT0I4GxKexNlwCKAketpXl5EHGyOKQLrWywybN3TlEhIIi5SPCOUKvfX+TlEYwv/b9Yoe00Dsu3+HIKyRrGBbOuo61uzTfJRg3qaTwKY+52z3ToBJHEG5xy+Pe2I6Pzg3nVapSyP1z/z4qOjCLJSOG7Njeai5xPUI+WZ8o9cTAmr9ussNjAs6pEdv1JuVgDL8bKe+foH8MhVUwAENvpjh/YOOCtx7ujMnklMqqTpJ0VXJ45HXBibleRCWXQOuvSwfOycOJjyU2jHafoWo4s9X4UAACAASURBVH7k85eeGPiCmGMrGYcuQdwsVSJPShaUqnMYyhjbAADhf84fjwCwRqi3NiwbER6r5dI1jLEWADsByF4/lYSVxbaDr6k2b2NxAt/83HKs27Efe5uSxIR/BJHOIWExFHxwttSgJvZcN2axGXXzp+6oTeIXMGbdtdmeoehkyMe/XIko+h/TjsRJgjKU4+OT4iQ3QxwSE8XjiT9MFo7vS2eMweiQAOWIK/7ja04a1R+LrgukoZ88+VDc8JFjgzbC80WDzmF/cyu27G6MnrPKAW0Lnahuff6dqEz3rZ+liJq092Uo5811r403I8Uiw8/DmEi3XTZRbkcRqY0Z0gvHHdxPqUPatLdqO6JORoVpgUw2JC+2G3cFm5/DBvcM7qlMncOvP30igNjp0eRspgsqaeoz0h0qN8d1ZObYShbOQTjF9SOAPdlSNVBphbRu3jJLue2aZONEM4hoHhHN27x5s65KKlyTu9t2QDYFcLwoMPxtQWD2qYqhRLGSWedg/qB42Rc+YM6HG/eVbCchunEgjKLoQMWhA3vgJ5ceZxUJHXlQbOFjakcNYshx8IDu2vI08A/z3c17UL9pj+AUx8LzFC5IAhtPJDlE6hL5mO7z4TfWCW3L59bpFpwSZcgmURgv4wHmxgzpJS0o54wbKreDdDl2oJOJ66ghwoPxBM+lpbWo1X/xR6H2NUnZCIi6uN++sALn3/wCgFhEF21iXMSA3BtbqMq5UlX/09RSxDf/FOu3dBZnpufEPadNXt+ZnOA0GNgz1gc2WELWVAOlEoeNoagI4f9NYflaAAcL9UYCWB+Wj9SUS9cQUQFAXyTFWAAAxthvGGMTGWMTBw8eXNLAzco5RaykU6wpDm7DNPloeTO3v7gi+pAG9KhN1OGtm3QOeYOsEogVarqwDyp4s6IiNMGlpLYSmxnqnt/9M07BqWMGJRbNiyYM1/aZdVF0DkiXuC7Axbe8LBWIYRSaWorYsDMOlqeOTHWYMnEOADCiX/doEVCfxRc+ECioJwg7c90mw+1d6MuvOe+o6PgD7xuMnrV56yJERA5y92S6WhXcWGHMfz6RiEbLzwPyt/fh44bjgS+eoowH0Qt4asl7UTmf56VYK+nmGtf/cMfU55ZuxIOvxYIN3fu1PaY12/YlnvOlE0dKY5Z0H0WGJxe9BxeIY+ksxOFRAJeHx5cDeEQonx5aII1GoHieG4qedhPR5FCfcJlyDW/rowCeY+UYtafASBwcrlWtg/RhJ4KyxxfGL1+1BhJ7M+kcbFnneFA9Xf93/9sk/M/FscUPryN6aI9W5PpugQBj8YwK0265V7fY/FQ8V2RB3gNXlKqBMvkcqKIfcTesTr2kjX5SsQ0Eu+CDB3RPxO/hGNm/B44Z0Re1hfiT0ytvjbcj34OmTOQM8mHmOZslG5GdCwaSnIO2HaTI0DULpNEcOLwzcb7Whc8sTecgRhmIiIOmHn8HnDCqY9eJm2xE9r1dDXhg3hqpjJuu62Ir3T9vDX7xXD1UxDKHGPkc4cmvnobJhw2IHO/aCi6mrPcCmA3gCCJaS0RXALgewDlEtBzAOeFvMMYWA3gAwBIATwK4ijHGyd2VAH6LQEn9DoAnwvLbAQwkonoAX0do+VQtiA9++kkHY+F1HwSgsfDQzAU+ofnHolPk6j7umkKykDFg+94mXPdo4Iavs1YC3EJ1iDj9fYPxCSERPZ9x4g5RNQt1M5IJrad04gxeg/T3AMgLNQPDL2cl8ydMPlyvatpoCIOdBnU8nEuzbT3Uc+rOz2gSm+Nxk0LioCEgOZI91fXvNv1lkGacn5p8iPSbL+qtliBuO/c3Y/W2fcbzgCsB0QfM43om3fzSEW4unnp+2SbMX7MjKu9eWwjPJxdaEV884/AoSOG9c4PFWveMRZ8TXR3d/dqC4eVzJBkaiNCJlURONQ1EhCMP6oPJhw3E/uZW7G1suyitqZ5FjLGPG05pzUYYYzMBzNSUzwOQ8AVnjDUAuDRtHJWCOLF61RUiCxEX8z/VocUlYJ2uHrcv/9ETb+GF5UEALzWPs6hILuQo2gFOP+ngRB0beN+2FI8uYhsi80fJn4dNlSHZaxf1C7Qu9DiQNPl0hXpXDWF8Kb5z1C1QCbESL2fxf505ZxCDqChwJboByWIcnfWJKfyG3EzSk/irU9+XGE+RMWOQOyAIB714vd2oQfXU1w9Ivwgff0j/aLyAyjkYiAOAz/zuVak8EiuFv02jqcnncPrYwXipfivunbsaP7rkGO0zVk1Z1aHruC3bJs2WSdGYQ9oCnfHI4YMDfdya7fsk/V010aU9pG0svCkpOxBPHpNZnwqdpy5TWH5d+AyAy7jjc9OEtIGu4iAgxcJKXdS19xD8103x4WH8Jl0wOl0fWZ15TETp1ZVa1VQEdTz7w8Bod3zmJFx91hgM1+iM1L5EG30gfB+ar4YHqDPpHIIymYOzea3boOMcVORzwXjKTVZPcBMr2RQBOnGQ9tuBXtHO35uq/9FBXdh1dWsihbQ752B7Bjb/FJ14OIuUlM8jPuZy32cWdDniIELOzia/MhvnwCeVViigefO6wHZM6T8ZPiP438pkpyuRw3DhHPh92RSPLpOV28SL39Fz3/gAlv1wWnR/CRNZkq/nyKpRMlWfv3qH4Qy045mzIiAmowf1xNc/eISWCKrfnmqjrxLrqK9Qxm/SOQTjIUmWrtuRP/PWRuP9qGOyIRgPi+bqdf8yLuUKU1+OCmmtU57MYYvNmESyup64Ej+yHIMh0xvkeT5r2Sbs1ohhuH7P9E3oysV3xSMMcKQ5btqiHYgRAaKyJ5dGx/zbimKktaHaocsRB/Edkbx6pYJPAmssI+11ah1K1DVxDqwoLwYioXHxqOVVrGIlh9WGy4PFST6od53k4Kd+JHL2u7g8O+egL1fDkKsoxds6wTkoYzCHSZAD/OmeKUF+D+XYXaRdyT3suc6hV4lxeYjsSm1eR7doMeF88NsuVjJhIDfoEAi1LtovEOsGiIDPKuIpjjgqawbOoVUcu3wu3RyYjHXU6L8quH4w8rCvnq1OAl2OOAD6l5yYqhaFdFRFK3rS7yoTzVs4ArGvNdv3SVYKIkFwW9QdOAeLriAeT+gBLJSpogFbO+KjYpq6Npg+B5uoLDEARySV9Xy3GnOLuoUtHy4AxthL4JxDmkI6HWqkVB3yFCqkwz5scvFy+7LlMgfid/3Ku7EY0LSxUR+JmHTJZYPBxa7/cuxwqXzWN8+IjguRE5y7zkFclBM5upVp+KFj5ARXtlA40v7U8or4GtLahqxD9lCXnRwmZ7Ckn0MSLkpr3QvWeSCri0zChDKcDGq2NDWFYxp4F+WKlbjSUZIbG7id6LcxWFlldA5pdt+lcA43f/x46XdkrRRFHTVEZaUgBtF9r642h6sg2cu1VIfX2kIuwQmq8btykc4hqJfGZf7hipO15S6cA5DkqE4a1R+HhSbTfP4vWBuLAfUK6aTVk0jURC7OtEYecVBvAEmHStF8W1VIJ01ZNblPxAixyjnRh+hHlxwjefQDIWcl+F3c9+rq6JzrHOUbsbZ0ku5ynIM4D3KaicehW8CyKIDl65Tf4CaR8Ykd+5uVa/QirIJEHFzGE9SxK6Qd22HK87OIkQD5mYqnKqVzSPtQsm6Wf3/FJAxVwnOIcm6AB97TcIcELH1vN7bsaUqci9oiOV7QvxyXTKHqgtpCDo3KzashzTknExlPpDyMfj30YiebSESsI1Y55bCB+NMXT42ek8s3AXDFtn6TFFwTc3GmMUWLqGXMNZHOIYNYSZJHK+eE+h89MZlfW3w+i9fvisKC8HMcTy0265vyOf16UE10PeIgPFsTywroiYO6+Onnn17WLP0O2xGbm3qUHNqAT/IWZYuUz0gceBVbXBabIjkqQ/BMxLSfOkW7iIG9YqckKd5PmdSBm7ymRefMyjjoTHpVhXRLK9MutG76n1ghfcdnJko5PrKgrpC36pB4X0XGIhl8loivIgjpCe+JgJ37m6Tf8nm3vnUKaVF0KXIO4mJ9nmDBJ2b/4+AxlcTx5HMUf1sOYiUpp4VyP/zcN855nzaN6v7m1sjQQP0Oxdei8y7n8MShDcAXlNsum4hzxwtmocrCYMsEp7YlwnVBYoxJPaqhMGL7aPk6WeeQ3g+v3txinlTJCatfJBmA7fuahHr6vjj+bcro6HhQrzr8+UunYtr4g4yRTU1Qn7PrlVnFSlqiGBGH4H9rkWnl966Emj9rVceUBbWFXKq+JRArxYQ4TedgGr8oEjGBCFi4bqexLV3wR12TpCnnZtK8HwBS9NtffPx43PopefHPkSwG6q3xn6nJx6HIdypcu24B5u19depYY9Rb2xR4d3OQ7zupp3MUK3mFdPXBn62ack99R7oJ4rQA6Mo0harOweRdrHIOWcVKvN2mVrN8XhSFDO5dhx9dkjSv46yxaIqpjln9XVB2Uccf0h896uzxflygLtjmehmJg6UNTqBaiixxX4BZvyLVEcajS83qirp8Dk0pllr5HDKJlUyPiohSdQ45IvSsFdKdGjiHtLSvOue+H14c+82KyX74fkZrHJCTHfd076uQy6G5lWHRup2Y+fhb0jmtKWuR4ciDeuOrU9+XuD9e31k8K/1OvQSAyDm0ndKh6xGH8H+ahY6OQqsso25xchb1MPvE4M2o4Q9MpqyfOXWUVqnIazRbRAP7BFvwH1w4PiF351i9bZ/Rvhxw14Fkliop9flHeOfLK63moJnj32uGr3JwLcWigXNwaF4kDmVwDjUFShUr5Ul2guNz5dITR+Li44No+d8XfB/MFlbpz5EgL6haM16lqEkzH1XOYepRQ9FHMMHl99DUUhTuSzdmWU+i1x8GBH/pe7sT50xOcLx/VcqgOunZoNZJiKgM35dXSLcB4sVEfinqzko3QeoK6cTBZYLkiULnNnPlyKeCmYmDOOTJhw3A+8fKeZ/FdvhiouZzBoApQr5o0yLB8zk/8uZ67Xl1PCZw3UUW2WmSOMTHuo+bY26KB3VybJpFLSzjMbBaWssQKwnH5XAOqjPd15TQGYDgBKfoHG689Djc9K8TgjE4iCgJ6ZyDyl1oOTDI39S9c1cn65D8rutq5O+ttpDDoF612LBzvxx2XYFqfquKjXgdZtigmUxZTc6e335ogXEsKtIs/EziQq9zaEOo71FlPXUvobaQ/rhsCwzH7oYW7GtqxfJN5oWNbyzVcZj8HPKGnWhEHMIthxo+HAjSnOra18G263flHIqM4fYXV0jlYkTNzykeqKqoQeznvJ+/YBmrfN1HTkhakojQDV8NPxKIlUojDuKjdfU7+NeJByfK8jmZcxjcW436Gy+Qsc4hOT/yDt72RA7hM0hWWut1N+lERiUr6mYMCLLl7WtqxeMLN4Tj1ouVxK7GDUvGIuJOi7r5ropygeAZRH1R8ntOjl4P20YH8MShQ0L9WHWstEocfvXJExJ1XDiHN0N775fqtxrrmLyxTeEzTAHb+Hh4iOy0HauL7NwE1ZxS337yA6nN5/DcN86Ifo/oLyf3SXxQjuNRv/FTDVFfre0qhS2tRe1CmyUIIuAuVhKtvMR2miQxTvI6Lne/Z86q6HeyTnxs4mJdw2c0Cw9bt2AzxtKdFiEbDg3vm0zyxGNG3RJm1NPmtA4J2oCetfj05EOlzY84RpNZss4664XlW7C7IfiG1Ct4Ey4bBFXEpV4jEv1HvzwlOubvyiukqwj+bNXXmBAraV5CIUfRInDkQb2lxC1ZoPvgVfBJq05UcU2RYjMZFhte5bFwp5W2KOkCoumgSzTU32ArL48n6S06vF839BWuVVtW34Srnll8hz1q87jkhBGW2m4KxZZi6aasYvOuYiU9NyNzDrqWAs4hziuiH7OY8U6PHMnP8YRD9HPe4gYAIGl1d8GxSR+PgAOxc0QF0QQV+s1MPse91c2Jmfg81J1XieGmMGT8O6HFkfpOLj9llLac4wunHxZxQd99ZJEyDuA75x8Z/eYhYX5y6XFRRr/gnoLr2zKPdNcjDuBmZ/KbVJXNakpFfg3nHkwft5PcMYMIQmXFZc5BIA5GhaJcnso5OK68D185JVHmnjRIMU1NsXpK0mm3MR4q7BiPGdE3dXy606oiv6WVaQljllAmgFms9KtPyNyo7r3mc3DgHORdqq4/F8s3InlB4spsEa+t2i79duE+P6YRlxFkxbaJCIuLt66rKJRJ0azX4/Nw1dZkPgv1m1M5HvVZxaasZu6Lt/iGEiwyRySZ1HOir+pbYr8nTxyqBhfO4egRfXDjR4/TXs9fUtpOXYRKeFyV1kCcg0A3TnEIpp2r2peaq1qF6fT1gnnrSaP645iRpTlw8eiuIpJKOvUqlZjo2554aH/p9xfPiHNsl2qGrC7OgdVKqWKl+NjEwR2khBE35Y5wcYITd/zadhwU0mo7LktT2ubnIyeMlIwgxDGIxFg3ZjG3iamvQBQGa77vHBF27GvGjX9fljhns8gDkvOEczI27stkVZcjWSfJCVGtsmbwqM2ec6giIuKgvElxF3XUQX2Mymc+2UwbcN1cnHH6YZnHGduGm8fpYjeflXMwfdhHj4iJgc4L1BU5Si4w6oegPkRXnYO6sNbkc9GCbCKelwg7Yd2t8zbHDOkFxhiaWovaueEkVkL6+0oQSsOiLouV9AQkTTztxDlANqd2EXmnPYvvnH+kXlcAcuIcRFGrXqyESBlvGkqOgL1Nel+RtJwQ6qMSI8HqQBSI1eo1BihEJF3HxUp1Sqwszzm0AWJDVvOiaZvcfDKadn7qR3bsyL5OOoZkO/pycWEW+7LZqYtIJQ6GjsX2Xay2TCBKRvlUOQl1BK6fg64et0IzfbhyBjV9pdPGDkLvboVoV6ezonGyVhIuM4mV0rzOgWAXul8MOqips2Hnfum3Lr2km55E5hxcHBjTxEqm+RNwDjFx0IvCcrLOQSd2IyFCrpGjTpZ/dsooAEmdQ7Ni2aBea8spz8sZY5j6039qzkEhDvo5FlkveoV09RDl+E1wDvGjsH00fMc2b5Xehl69toeGMLi8Xx2Refmas6QPy800Ui5PEyuZ5KbifZXHOSR3tKpHatpCa9w8aco5V2KTqeuO5TqB8pLL+dOIw2kafxNA5RzcNhe6cat2+7phq3LyUUJUUg6RSzQvovJu1W3u2s+LOUDUvkzh6cUycfHWPUZO0AJrJMMYc0kRzTXnBYphdXeuzk+1zdhDWt8XwZYHhKR3HImVlDlW6GwKaSJaSUQLiWg+Ec0LywYQ0dNEtDz831+ofy0R1RPRMiI6Vyg/MWynnohupqxxDzLA9GhdOQe+YzO9I3HxHda3G34+/fhEHZfXq67h+RxJcWYA1c+h2pxDfJwQA2VAwGLLT0AX4ljEoQNVU0T9E9Qmkw/v13hfotxdWyO8d8ZiZaGWOMTH370gPeuaC4em+w24ZQMTF/4Zpx+GQb2Slj9i9kGjVQ/kBakSOgedQj/oiySOSKtzyMs6BzX/Br+uWGSpOgdRqT+iX/doAVYtBNXfqtQhFiu5cyniOMT7NM2xzipWOpMxNoExNjH8fQ2AZxljYwE8G/4GEY0DMB3AeADTANxCRPzN3gpgBoCx4d+0CoxLD4POwUVE4wLxg/uPaUdqQ1G4ZABLJBTRXCL7ObhNTN3HJMJo9SRxDqU/H1XncOsnT0ikXRRb/99Lj4t2dBymx6cr3xXusk2LnxT10yJ3LzLgD6+sAqDnnFyIDDcnDtowPWdlfJqBJ8RwWoum9I2DHPHUvIjKnEPyIasWVmliJeNimeDm0zkHNY8FwB3ckKJzkPUbN156LPI5Cn0kZOqrBqZUnTLTFNK25YQIWs5B5a74vOg0nIMBFwK4Kzy+C8BFQvl9jLFGxtgKAPUAJhHRMAB9GGOzWTDz7hauqRqSge7i41LDGwNAjWRNpG+npNebEkLcVYadJhIyTWRxIXFxdjNBjXszcdSAxLsYJ4RSP+vIIYkPxST31kXJ5d/S4vW7DOOJj80fdxAQ7mfPLI/uIdlOOpH5kGDb76qz0r3WNB0NoBAHw3hkhbS2SmI+6DYX6rWuvjKJvhLtajgHRSGtKm6B4LvjYiWbOFFqh5uo55Ke3CrnkBSLpuscTMjlZM5hZxj1uHc3+Rvj86Uz6RwYgKeI6DUimhGWDWWMbQCA8P+QsHwEgDXCtWvDshHhsVqeABHNIKJ5RDRv82Zz7HP7gPUPV9w5lUMcsmZqA4KgeSre2iAvZpdr6uQcCFFWLshF3KHzOHWF6gRnMtV0PS9Ct4vk2LCzwTie+Fh/bWCKKF6kr6M7FnHKYbGHtrNCWlPPJYSCy9zIO1BG8fmcNnYQ/vWkpH9COV71NuhEoPkc4b1d8bvU+oGECuAiMxPqpDlwLH406RyGh5Zr6u494hws80cFF80S5PewbkcDavKEIYoDoCmcTjVRLnGYwhg7AcB5AK4iotMtdXWPjlnKk4WM/YYxNpExNnHw4MHZRwuznwMZPI+zQlQ0mnZQn3t/bNo6/aSDcd2Hxyfq8HAXAPD1c96HSzRxgdwU0tlgls3HxzqZuytEm+9DBvTAgJ7JWE/iY9ON54+fn5wo69u9RvscXcYT9WvmHSSClrZTN80fcdftSsy1Ogfl60jqZGQO1jQPXcYsll7x/tF6kZpqemzhjU0Z54DkQq7zJynkcti2Nw4xP7SPPq5US7Sb1/eVI1nRzOsVcrkEp8Dr/TIMlyMOs1ddISIyZu5LPvF/nzoR1314PAb1qsOQ3nXSe2hsaUW3mnzyWYS/Ow1xYIytD/9vAvBnAJMAbAxFRQj/bwqrrwUgbjtGAlgflo/UlFcFkSmrukNzENG4QI50mS5eOMUQ7+fTpxwaHZvEOGk7bEBehD41+RBtHVObpvJyng+3+SYCLpwwXFsnzfOb5wkWMfPio9FPE1QwDS7Jk0SCFtRLVhTLzMQh/XNTr/3A+5KbINGU878vHI+JowZYx2OSJEr37jAek+Ndlulg23ipjEKNpmGRmxjcu077LnLC7t+mS2nShBkv5Cmhc+BiI77bHz0oyE990qj+OHRgD+zazwPvmUVYInrU5vGJkw/BvP83FYV8Tprjza1FrcEHf1edQiFNRD2JqDc/BvBBAIsAPArg8rDa5QAeCY8fBTCdiOqIaDQCxfPcUPS0m4gmh1ZKlwnXVBwx5yC/MfGXSQThgoKDBYhU3yB7FnfUpphFTsRBKD5sUC9tHRFG5aWk4CzPlLU1tCRx8ep27apUPYhb2AtZzq+7QuJADE12M5hwyu3EFy/74TQcpYkoKpp7HnlQ8rzajul9mUKxiBCLTTlBsoiVbDVVrkSnT5BMqo1zNSZkJus8IsWZLuIckmKlz989T2qrR10wru61BfTtXoNlGwPnNjWMSNRXSnoA8fU0tzAtd0ZEGNCzFuu270+cqxZK1ywCQwH8Ofy4CgD+yBh7koheBfAAEV0BYDWASwGAMbaYiB4AsARAC4CrGGPcbu1KAHcC6A7gifCvKohjK8nl4sexdW8jSoWrv0RU38Gzub9hR5wmfgHkienyDZtERlmzmJnEB+JibFKOu8SMUlFXonmtJJqzmFhKHsmaai4+A3mH5ya2bdqJijJvFwJrejSmII5SHZFzMGQTVK+dtaw0fWCNMvd0Fl0i12p6nqI+wWgVRpB0F1HkA0XhLXKMfL7yb6SxuVX6NhcIqVLVvkSoj1pcM5pbi6gp6Md8UJ9u2twU1ULJxIEx9i6ARAAixthWAGcbrpkJYKamfB6Ao5NXVB5GnQMl65QCFzmuVN9BHmxaaN0C7wnHDtTBTBziYxexkil3Ql7a0aaP2dU4QA1U5gqxfaOSOCfvmrUhoB0MEVyem4tBg+hx76LYNhFhiWs19LV2exyY7viD+2vrqHNv8+7SNleqOMWUq4HDZvHFw1C4WoWJbYqcg3gv3IeJfyNNrUX07R5/mz0NkRDUvtTnJb6rh8OkWjoEZrw+TWjVEL32FGpeCbgsxi6moybOIaufg4sIxRTawMVuXoTqm8DhQmSyjhkAavPZQ5QE40kXv6gxf07SyPhd/GRcnpvMoenH85OPxXuycsKdiM2bHBtXhlFLjzyot9bLOuhLW6zFVkGZrEIdp+4x3v9qbPBomj/z1+zAm2uDXfzmPXpCZdrNqzoHMdMg39HzzcH5Rw+TQs1/71/0zo8JKYXSuescD0x0napWBF2OOHC2QGXZy7FQMqEcAyJxAnEZpwpxYbDFq8kyHlNog5zDjk2EUWTkslMv4V2UyjlI4heL3F11hFIh3oppoXV5bi7vSNQzGEV8QrHpnYq6izRRoU7+z6EudmcdOcRQ0w5RBHTmEYNxlEafkmYGreJtQxpZU1j4fI7QLHQiKv95Pushvbvhrf+ehs+dNlritPt1N4l/5b7UcByuKOSSccmqiXJ0Dp0SLtZKLh7AH5+UtPdWUSmCY4qHJE4604KUVX5vltFm0zmYIC3GxjFnb1cXw8oF0vOxKC+bU0JkixE+Vdl51H5YPKxv0mte7CsLXMQmpo1DDwfxFMehFt8WdZ5PPSqZC8UFIvd33YfHp3LeLvPQ1EYvxYCBi5IKOZKi0HL9w8+nT5CMHrhoT+KojTorGWnh1v9jmp7rzhNpU5hWC12PcwihvjBxDvWoTaeZXzpjTGqdihEHB78CkwhCHEGzZdfx28sm4kPHDDMH3quQKaubV3f29kf21y9eOgdDES73RTBb6nC8JgRiNBFqFyuvrA6Ypoi/4n2Z9EiHDY6t19KIkklMCGgc90qcHt1rsxlzuDxPUzNqljkuSlJ1DtyvQGc+Dcj6SVfHxrR0qWKEAKmdnFtcrUqhyxEHk7JZ/DiG9TPv7Dhc9AmVchyt1GL890XvGetNHTdUmxObQzYvLWc88bFZ55CtzbMtYow0xzhZ7Gb+uNPESledGW8WTNwXv1/b7WV9tAM1ToSA7D9TToh1DjXoJHGu+AAAD2dJREFUowh1PpS6KRLjkJmIg6j0NZmyis/ENJaeiqiWv15V5yByFDqkhQ8HkoRXRxzeLyQ/MnHBhVwOm/c04pjv/x0PvLpGW6eS6ILEwR5BEQCuPmusczs2ZM0rbEI5IbJFz+9yvCsl/UaZIbs5ykm1KuI/LLtaABg/vA+uNdRxSd3pIuoRLXlsIRvS2ssqVjIFUhQzrZXj0a6L5qoisZhpbuGMI9IjGnCZPmAWgd7/hVOiY9Mm5fvChuB9Q/U7ftWAgX8bavgMvvibuBSxrivncPLopEHDsUJmRVMYmFyOsGLLXuxubLF6oVcKXVfnYKmTFrkUcPtonKyDHKxsyon1JF75lanpRM84BgdrHKfxOFgHZW0/rf5jV5/m1I5Z5xAfXyB4t0vXZvBpsQ23nHdt7LcMp8UXvn1mqpz7ICXysO593HH5STjsO49b25EcSA3PQRQHGcV3Qv//fpZe/KtyU5w41AjhMxat24n6TXuCsZk4BzGtqUVnxfH2D8/TcnLiMzOtP2LzleAG09D1iIMhZLcI2/e54LoPIkfkREBcvvMpY/ThM7LANlHESTf5sNL7cnG4c2pHODaJBbI2XylDM6O1kjDqfzlOH/LDRfTHF5KeFp1WFWhDWQYEgU7DPtdVvcdxmvziLqLIGoe4ZHJ4Gn07Is0w6eLUb2b04J7htXFI8At+8aLQjok42DPXqeN0CXNjCrMiOheaLNAqia5HHML/5iBr9h2/yPqmwTUNYzlYcN0HrdFIK6UUd4nV4wIxDIXx+WQlDqUPR4KLt7E5D0P6KPY3BzF4TKbJQHVMqsdpQnBwfOvcI8omriJnUj/zPOOCnAbx2Zo2X6Jxxg8u1PvNujzD2rCv948ZhDs/e1I05kKe0NCS9AQ3zQ1RrOTiZW6aJmK57XvmKEfU7IquRxxMLtJVQLkf+udPG413Nu+11kkjVlVx7iujUVGxWymdQ6UWVBeZcTkimrFDe6N7TR7fOOcIY51qvC8blysq0kuF5LVcxqIlPluTFZZo1m12ygvGYzNJ55xDa5HJkZRzpNXNmd67KFYyLdguCaXEOeySc95F51kuuhxx4KjGR8gxol93rNuxv6ykOADwnx9KTzeZhmrc5yFl5HNobhEVeJXROVTqHs3mwOlKaxf06VaDt35gT3JYaqKcAwF8s2CywBLr2MCJlW13zfUVavIcXchuW7/Ngj7GaGHlsA6I17oEaGwLX7guaK0U/K/mJ7i7IQiO1a+7eVf/+ysmJdIrVgPVEFMMdFDGmyAGb3O17tDhns+dHB3bRISVgPgIy9kZu6Aa76uzwMUKjs8ZW/j5nAtxCBdg1eO4YOAcTBZfvO6phtD7QNJsVgcXU3FRT3L2UaV5oWdBl+Mc4qis1fsIb/jIsfjp029LQblUnDa2tGRFWdHRlhrRmcw177UOYqasaq+nkrlrGcpdt76C/8cf0q+q/bQXJmniUnG46uiW/XCaMWoAEHNfduKg5xzyeZK4gai+oS3OZXx8kplY8ZSftvvjmwJbzo8PHzccD74WJM30Oocqopqf+HnHDMN5x+hNHtsaHW0nKoYOKGcXXkmues53zsaexhbjeUkhbVmUTjlsIKaOKy10RNwX4fGrT8MhmuxunR0Lr/ug1cpmUEjwv3muWScDpFvq8FdU66BzMHEOqkw/bcNiW6x71QWbRF3WQw5O0Ab2NHPl3z73iIg4tAW6HHFow/zcHQIdjDbICmnDTsoluJj4oaWFI0jD0D7dYFvSxWHanue9M5LpS0uBKXxCZ0fvFOOJPt1qsOJH55fN1Uecw/9v7+5j5KrKOI5/n3a7xXahQEtf2HZbIi/SFi20FKoFCklLQbQaNdkGaaEmiAEBNQYQlUTRgFEiFBJotA2CKDFiLCISMEKCitDKaynlLQ0WCG8CbalSio9/3HPpZWfu3dmdmTtn2N8nmczdM2dunn12ds495957TsEl3ulnLJ29NZWec9jwwtZqb6twyaemM2bUiMLJBtNhpaKTyOmvPLYrvwGpZzh3MIZe4xCeY/vSbBYz45zjD2TRjPqOaBsl23PIO9qq5QTeuK6RjOvq5NXtOwc9y2Wt1m3evcLXjp3VF7yJ0YHju/qd5C02jRjureWEdN5qj8lKcP9j82vFVwmmxu+1Bz/87GGFddKrxYp6POk5jb4TAmY14wbJIkOucThp5kQOntBVyk0kseivm16r61fMLZydsxbZL/K+k5+lOjuGsfmyT/a7r8UzJ3Ljfc8V/kM1QnqXLED3PvlzDJXpyt5ZvLGjeFWwu75+XEnRJDcHzoikx5Oe0C26siw9D9CX47y09W3OuenBhsUzdnQnXz3hwNw10wEmh/+r/7wTz8HHkGscpo4dzdSx1a+PhuoLukuiEbnJHsmOK+hC1+K7p8yg98ieqiuzNdLIjmHsCj2G7oIJ6Mq0ZFZ3q0N4n5VLD291CO9J53oq6lH2HjmFbf/dxbJ5U99Xfss/K1diu+3c+XXFY2Z8Y1HxAdr+Y5LPVdEVjmWLpnEws8XAlST36v/M3S8rO4ZNly6u6yYngZvPPLpwOo/TPz6N28PssPUOIXR2DGNmd+VUDY32oc4O3tr5LnOrTJgm8UnnPZuwV/7syh3Dh/GVBR+uKO/7kaylB9sIB0/o4ntLZrB45sTCerefd0wp8ypBJI2DmQ0HrgEWAluAB8xsrbs/XmYcQ2moqVmO6mf+pv5ej1F6eeGlnyllmXOp07iukaw6bTZzCi6bzbPHiOG8825y5VrRzXiNZmYsmzet33qHFkyF0mixHCbPBZ5292fdfSfwa2BJi2OSJrnj/GO594LjWx1GzdKV26rdOStxWjRjYuGlo3kuPvnQ97YXHNL8G81iFkXPAegGsqtXbAGOyqkrbS5vVa1YrVx6BDfct5mPtFncMnC9c3vondvDs69sL1zgaCiIpXGoNvhccZhmZmcCZwL09OTfkSjSSBPH7ME3TyxeUEg+WLJLqA5VsQwrbQGmZH6eDLzQt5K7r3L3Oe4+Z7/9dFWRiEizxNI4PAAcZGYHmFkn0AusbXFMIiJDVhTDSu6+y8zOAe4guZR1tbtvaHFYIiJDVhSNA4C7/xEoXmRWRERKEcuwkoiIRESNg4iIVFDjICIiFdQ4iIhIBStagCJmZrYN2JTzcg/wXD+7GAO8GVEdxVxOHcVcTh3FXN8+aq03mJgPcff+b/d397Z8AOsKXnulhveviqyOYlbMinkIxFzLPpoZc9F3Z/bxQR1WeqOGOrdGVkcxl1NHMZdTRzHXt49a6zUq5grtPKy0zt3nDPS1WCnmcijmcijmcgwm5lrf0849h1WDfC1Wirkcirkcirkcg4m5pve0bc9BRESap517DiIi0iRt0TiY2Woze9nMHsuUfczM/m5mj5rZrWa2VyjvNLM1ofxhM1uQec/sUP60mV1l9S5iXE7Md5vZJjN7KDyatjyVmU0xs7+Y2UYz22Bm54Xyfc3sTjN7Kjzvk3nPRSGfm8zsxEx5KblucMyl5HqgMZvZ2FB/u5ld3WdfUea5n5hjzfNCM1sf8rnezE7I7CvWPBfFXF+ea7mkqdUP4FjgCOCxTNkDwHFhewXw/bB9NrAmbI8H1gPDws/3A/NIFhe6HTipDWK+G5hTUp4nAUeE7T2BJ4HpwI+AC0P5hcDlYXs68DAwEjgAeAYYXmauGxxzKbkeRMyjgfnAWcDVffYVa56LYo41z4cD+4ftmcDzbZDnopjrynNT/zgNTto03v9Fu5Xd50ymAI+H7WuAL2bq/ZlkjepJwBOZ8qXAdTHH3Ig/cJ3x/x5YSHKz4aRQNgnYFLYvAi7K1L8j/AOVnut6Y25lrvuLOVPvdDJftDHnOS/mdshzKDfgNZKDiOjz3DfmRuS5LYaVcjwGfDpsf4HdK8k9DCwxsw4zOwCYHV7rJllxLrUllJVpoDGn1oRu4Xea1Z3ty8ymkRyV/AOY4O4vAoTntHtabe3vblqU6zpjTpWa6xpjzhNznvsTe54/Bzzo7m/TPnnOxpwadJ7buXFYAZxtZutJul87Q/lqkj/eOuCnwN+AXdS4TnWTDTRmgFPd/TDgmPA4rdlBmlkX8FvgfHffWlS1SpkXlDdNA2KGknM9gJhzd1GlLJY8F4k6z2Y2A7gc+HJaVKVaVHmuEjPUmee2bRzc/Ql3X+Tus4FfkYwd4+673P1r7j7L3ZcAewNPkXz5Ts7souo61ZHFjLs/H563ATeRDJE1jZmNIPlQ/tLdbwnFL5nZpPD6JODlUJ639nepuW5QzKXmeoAx54k5z7lizrOZTQZ+Byxz92dCcdR5zom57jy3beOQnnk3s2HAt4Frw8+jzGx02F4I7HL3x0NXbJuZHR26V8tIxvOijTkMM40L5SOAU0iGppoVnwE/Bza6+xWZl9YCy8P2cnbnbS3Qa2Yjw3DYQcD9Zea6UTGXmetBxFxV5HnO20+0eTazvYHbSM5J/TWtHHOe82JuSJ7LOKlS74PkKPtF4B2SVvxLwHkkZ/KfBC5j94neaSQnbzYCdwFTM/uZExL0DHB1+p5YYya54mM98AiwAbiScGVNk2KeT9JdfgR4KDxOBsaSnCR/Kjzvm3nPxSGfm8hcwVFWrhsVc5m5HmTMm4F/A9vD52l6G+S5IuaY80xywPZWpu5DwPiY85wXcyPyrDukRUSkQtsOK4mISPOocRARkQpqHEREpIIaBxERqaDGQUREKqhxEGkCMzvLzJYNoP40y8zgK9JqHa0OQOSDxsw63P3aVschUg81DiJVhEnP/kQy6dnhJDcuLgMOBa4AuoBXgdPd/UUzu5tkTqxPAGvNbE9gu7v/2MxmkdwNP4rkJqoV7v66mc0mmVdrB3Bveb+dSP80rCSS7xBglbt/lGS69bOBlcDnPZkfazXwg0z9vd39OHf/SZ/9/AK4IOznUeCSUL4GONfd5zXzlxAZDPUcRPL9y3fPV3Mj8C2SBVXuDLMfDyeZIiV1c98dmNkYkkbjnlB0PfCbKuU3ACc1/lcQGRw1DiL5+s4tsw3YUHCk/9YA9m1V9i8SDQ0rieTrMbO0IVgK3Afsl5aZ2Ygwj34ud38TeN3MjglFpwH3uPsbwJtmNj+Un9r48EUGTz0HkXwbgeVmdh3JbJgrSZYVvSoMC3WQLM60oZ/9LAeuNbNRwLPAGaH8DGC1me0I+xWJhmZlFakiXK30B3ef2eJQRFpCw0oiIlJBPQcREamgnoOIiFRQ4yAiIhXUOIiISAU1DiIiUkGNg4iIVFDjICIiFf4PiCFO+dEhlugAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sorted_data['inc'].plot()" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXt8HGd97/955rp3SaubbVm2fE0cm8SJHcdJSAINadJACddfA20SOLShHNrSA6enTXvacjn8Cj2nDU0LobRQQsohBGibEAg0NyCE4MTOzfElvsSWZUvWXdpd7c7uXJ7zx8wzO7vam6Td1e7qeb9eekmenVmPVjPzfT7fK6GUgsPhcDgcL8JynwCHw+FwGg9uHDgcDoczD24cOBwOhzMPbhw4HA6HMw9uHDgcDoczD24cOBwOhzMPbhw4HA6HMw9uHDgcDoczD24cOBwOhzMPablPYLF0dXXRgYGB5T4NDofDaSoOHDgwQSntLrdf0xqHgYEB7N+/f7lPg8PhcJoKQshgJftxtxKHw+Fw5sGNA4fD4XDmwY0Dh8PhcObBjQOHw+Fw5lHWOBBC+gkhTxFCjhBCDhFCPuZs/yQh5Bwh5CXn62bPMXcRQk4QQl4jhNzo2b6LEHLQee0eQghxtquEkG872/cRQgaq/6tyOBwOp1IqUQ4GgE9QSrcB2Avgo4SQi5zX7qaU7nS+fggAzmu3AtgO4CYAXyKEiM7+9wK4E8AW5+smZ/uHAExTSjcDuBvA55f+q3E4HA5nsZQ1DpTSEUrpC87PcQBHAPSVOOQWAA9QStOU0lMATgDYQwhZDSBCKX2W2uPnvgHgHZ5j7nN+/i6A65mq4HA4HE79WVDMwXH3XApgn7Pp9wghrxBCvkYI6XC29QEY8hx21tnW5/ycvz3nGEqpAWAWQOdCzo3D4XCWQlzT8dBL55b7NBqGio0DISQE4HsA/pBSGoPtItoEYCeAEQB/w3YtcDgtsb3UMfnncCchZD8hZP/4+Hilp87hcDhlefTgeXzsgZdwflZb7lNpCCoyDoQQGbZh+Cal9N8AgFI6Sik1KaUWgH8CsMfZ/SyAfs/hawEMO9vXFtiecwwhRALQBmAq/zwopV+hlO6mlO7u7i5b/c3hcDgVoxkmACClm8t8Jo1BJdlKBMBXARyhlP6tZ/tqz27vBPCq8/PDAG51MpA2wA48P0cpHQEQJ4Tsdd7zdgAPeY65w/n5PQCedOISHA6HUxd0037kZAxrmc+kMaikt9LVAG4DcJAQ8pKz7U8BvI8QshO2++c0gA8DAKX0ECHkQQCHYWc6fZRSykzxRwB8HYAfwKPOF2Abn/sJISdgK4Zbl/ZrcTgczsLQTdsocONgU9Y4UEp/jsIxgR+WOOazAD5bYPt+ADsKbNcAvLfcuXA4HE6tMJhxMLlbCeAV0hwOhwMAyDhupTRXDgC4ceBwOBwAHuXAjQMAbhw4HA4HAI855MONA4fD4cCTrWRy4wBw48DhcDgAssohrXPjAHDjwGky9r0+iVfPzS73aXBaEIMrhxy4ceA0FZ/6/mH87WPHlvs0OC0Ijznkwo0Dp6lIZgzEUvpynwanBdEtXiHthRsHTlOR0k3ENWO5T4PTgugGK4LjxgGorH0Gh9MwaLqFhMCNA6f6GJYTkObKAQA3DpwmI6WbsHhPRk4N4I33cuHGgdM0mBZFxrCgmxYsi0IQ+LBATvXgAelceMyB0zSknX77lAJJ3nOfU2Wyqaz82gK4ceA0EalM9qaNazxjiVNdMlw55MCNA6dp8E7oys9YGo1pOD0xV+9T4rQQLCDNjYMNNw6cpkHztDXINw6f/v5h3HzP09h/et50WQ6nInSDV0h74caB0zRoenG30uRcGsmMiQ/8y/MYnOQKgrNwdK4ccuDGgdM0aCXcSnNpE33tfiTSBl44M13vU+O0AG7jPW4cAHDjwGkivDGHRDrfOBhY3eYDkOt+4nAqxeCT4HLgxoHTNJTKVkqkDXSGFABAmqe5chYBr3PIhRsHTtOgeW7axDy3koHOkDpvPw6nUniFdC7cOHCaBs2jHGIe42BZFHMZE51BWzloXDlwFoGrHHi2EgBuHDhNhOZUSPtkIScgzaqlwz4JiijwmANnURhcOeTAjQOnaWAxh+6wikQ6G3OYc4LTQVWCKgtcOXAWDKWUV0jnwY0Dp2lg2UrdITVHObDMpZAqwSeLbg8mDqdSTCvb6Ze7lWy4ceA0DZpuQREFtPnlHOPgKgdFgk+23UqabuKmL/wMz/OKaU4FsGA0wJUDgxsHTtOg6SZ8soCwT86pc0h43UqSrRzGYmkcPR/H4eHYcp0up4lg1dGKKHDj4MCNA6dpsI2DiJBPyqlzmEvbbiTbrWQrh6RuGwzuYuJUAhsRGlRFZEwLlA+U4saB0zykdBN+RUTYJ+WksmYD0iJ8kghNN12DwTOXOJVgODGHgGLPP+NxB24cOE1EKmPCL4uI+GRkDMtVBfkBaU033cwmrhw4lcBcSSFVyvn3SoYbB07ToBkWVFl0b2BWJe1NZXXdShl7G1cOnEpgyiGoigC4cQC4ceA0EVrGhF8WEPY5xiGdNQ6EAAFFhCqL0AzTTXvlyoFTCaw6OqhytxKjrHEghPQTQp4ihBwhhBwihHzM2R4lhDxGCDnufO/wHHMXIeQEIeQ1QsiNnu27CCEHndfuIYQQZ7tKCPm2s30fIWSg+r8qp9k4PBzD795/wF3FaYbtVmLKgaWzJtImgooEQgh8koi0biHJ3EpcOXAqwDUOCncrMSpRDgaAT1BKtwHYC+CjhJCLAPwJgCcopVsAPOH8G85rtwLYDuAmAF8ihIjOe90L4E4AW5yvm5ztHwIwTSndDOBuAJ+vwu/GaXK+uW8QPzp0HkPTSQB2zMEni+iJ2K25z89qAGzlwNwBPqdCmhkH3oSPUwmszoEpB962uwLjQCkdoZS+4PwcB3AEQB+AWwDc5+x2H4B3OD/fAuABSmmaUnoKwAkAewghqwFEKKXPUjtP7Bt5x7D3+i6A65mq4KxMKKX4yWvjAICZZAaAk60ki1gfDQAATjsT3xIZw72p7ToHCykn5sDbd3MqwTCzqawAVw7AAmMOjrvnUgD7APRSSkcA24AA6HF26wMw5DnsrLOtz/k5f3vOMZRSA8AsgM6FnBuntTgxlsC5mRQAYGrOrmnQdDsg3R6QEfFJGJy0FcVc2nBdTfnKga8AOZXAlcN8KjYOhJAQgO8B+ENKaamy00Irflpie6lj8s/hTkLIfkLI/vHx8XKnzGlimGoAgGlHOWiOciCEYKAriMGprHFgvmKfLMKwqBuP4E34OJWQjTlw5cCoyDgQQmTYhuGblNJ/czaPOq4iON/HnO1nAfR7Dl8LYNjZvrbA9pxjCCESgDYA85riUEq/QindTSnd3d3dXcmpc5qUnxwbQ1+7HwAwPedxKyn2JbsuGsAgcyulTXfF55Pt15lB4StATiUYFs9WyqeSbCUC4KsAjlBK/9bz0sMA7nB+vgPAQ57ttzoZSBtgB56fc1xPcULIXuc9b887hr3XewA8SXn9+oqFUooXz8zg+m09kEWC6aQO3bRgWhQ+yV7ZDXQGcXY6Bd20HLcSC0jb371qg8MpR8Zw3Eo8W8lFqmCfqwHcBuAgIeQlZ9ufAvgcgAcJIR8CcAbAewGAUnqIEPIggMOwM50+Silld+hHAHwdgB/Ao84XYBuf+wkhJ2ArhluX+HtxmphE2kAyY6Kv3Y+OgILpuYxbt+B3ZP/6zgBMi2J4JuVkKznKwTEekwnbOPCbnFMJ85QDv27KGwdK6c9ROCYAANcXOeazAD5bYPt+ADsKbNfgGBcOZzSWBgD0Rny2cUhmXAXAlMH6ziAA4PRkEglPQFrlbiXOImAxhwDLVjK54uQV0pyGYyxm1y/0RFR0BGXMJHVoGfvmZcZhoNNOZz05lkDasHJSWQFg2s1w4jc5pzxuthJ3K7lw48BpOMbiucphKulxKznGoTuswi+LODxiJ87lB6RZQJErB04l6LzOYR7cOHAajlGmHMIq2gMKZnLcSvYlSwjB+s4ADgxOA8C8gDSDKwdOJRh5yoEvKrhx4DQgo7E0AordQykalDGd1N2iNr/n4f/OS/twasJOZ80qh1zjYFjUrX7lcIrBG+/NhxsHTsMxFtfQG/GBEIKOgALTohhP2K4mn5J9+N957Ubcee1GAEBXSLVfl+df0vxG55QjWyHN3UqMSlJZOZy6MhZLoydsP+w7AgoAYNhppcFSVQHbtXTXr12Id17ahwt6w/Neb/PLmE3p0HQLzttwOAVhykERBUgC4cYBXDlwGpDRuOZ2Xu0IygCA505NgRC4VdMMQgi2rY5AEOxsa69bKRq0LQKf6cAph2FaIAQQBQJFEnjMAVw5cBoMSilGYxp685TD08fHcdHqCNoCcsnjVSm73ml39uXT4DjFGItrePrYBDImhSwIIMQ2Dlw5cOXAaTDiaQOabqGXKQfHOOgmxd6N5Rv15iiHAFcOnNI89OIwPvGdlzEW1yCLtvpURG4cAG4cOA2GtwAOyBoHALiyAuOQqxwc48CVA6cILAtuZEaDJNrXjioLfEEBbhw4DYa3dQYAhH0SRIFAIMDlG6JljxccnzEARIPMrcRvdE5hNMcInI9pkB3jEFJldz75SoYbB05D4S2AA+yHfbtfxvY1bWjzl443MHyOcXCVA3cRcIrAFg4jsynXrdTmlzCb0pfztBoCHpDmNBRDU3bKKlMOAPCbe9djU3ew4vfwySJimuG6pLhx4BSDJStouuUqhza/7BZXrmS4ceA0FAfOTOOC3rBbqQoAH79h64LegwWluVuJUw5vbEFylYPMlQO4W4nTQJgWxQuD09g90LGk91G5W4lTId5kBcVRDhEfNw4ANw6cBuLo+RgSaQOXD5QPPJeCKQfmVuLKgVMM77XhVQ6abq34jCVuHDgNw/On7LHhlWQllYL1V2LV1Vw5cIrhvTbcmINTPBlLreyMJW4cOA3D84PTWNPmm9ciY6H4ZBECsd0DAC+C4xTHqxxkIRuQBrDiXUvcOHAahv2np7B7iS4lwJ4GF1AkN/bA22dwiqF5Fg6yZLuVItw4AODGgdMgxDUdo7E0LloTWfJ7+WQBfkUEIQSqxKtdOcXxLhykPOUQW+HGgaeychqCc05L7qW6lABgz4aoG5RWJYG3z+AUJW2YIASgFDl1DgBXDtw4cBqCc9OOcehYunG4/coB92dVFrly4BRF0y30hn1O+4xsthLAjQN3K3EaAqYc1lZBOXjxyVw5cIqj6SbWOgsSrhxy4caB0xCcm0lBEQV33Ge1UCUxJ+jI4XhJG5arVlmdgywKCCgiNw7LfQIcDmC7lVa3+9yJbtWCKwdOMSyLImNYWOOoVVYhDfAWGgA3DpwG4dxMqirB6HxUSeRFcJyCsOsi4pPREZDhV7KDongLDR6Q5jQI56ZTuG5rd9XfV5UE3j6Dk0Nc03FmKok1bfZixCcL+Oc7Ls9ZnLT55RWfysqVA2fZSRsmxuLpqmQq5eOTuXJYyZgWxV3/9gpOjifcbV/9+Sm8+95fIOUsGlRJxK71HVjVlm0TH+FuJW4cOMvP+Vl7wE9t3EpcOaxkzsc0fOu5Ifz40Hl32+BkEppuYTKRAZDtxeWFKwduHDgNgFvjUCPjwJXDyiWVsZvnsWsMAIadtOnxhL0oYQWTXnhAmhsHTgNwdqZ6BXD5+HgR3IomlbEXBqyOBgBGHKU6HrfnlRdTDnMZE7q5chcW3Dhwlp2hqSQIAVa3Vd84hFR7HnCGq4cVSdJRDkwtUEpdNyYzDqpUSDnYuTor2bVU1jgQQr5GCBkjhLzq2fZJQsg5QshLztfNntfuIoScIIS8Rgi50bN9FyHkoPPaPYQQ4mxXCSHfdrbvI4QMVPdX5DQ6+09PY/uaCBSp+muVvRs7oekWnnNmRbQqsykdPz8+sdyn0XCwoPO56RQopZicyyDjqIGSyiHAq6QruRu/DuCmAtvvppTudL5+CACEkIsA3Apgu3PMlwghzCzfC+BOAFucL/aeHwIwTSndDOBuAJ9f5O/CaUI03cSBM9PYu6GzJu9/9eYuqJKAx4+M1uT9G4UHnjuD2762D5OJ9HKfSkPBkhHmMiZmUzpGZjT3tfFEceXQGbQr9SecoPVKpKxxoJT+DECly65bADxAKU1TSk8BOAFgDyFkNYAIpfRZSikF8A0A7/Acc5/z83cBXM9UBaf1eWloBhnDwt6NtTEOfkXENVu68PiRUdiXXmsyMquBUuDk+Nxyn0pDkcxk403nZlIYmc3GHkoph96IndY6FtfmvbZSWIqO/z1CyCuO24lNhO8DMOTZ56yzrc/5OX97zjGUUgPALICCTwpCyJ2EkP2EkP3j4+NLOHVOo/DL1ydByNJHg5bi+m29ODudwrHRRPmdmxT2EPPm83OybiXAdi2xYDSQVQWFlENvxFYOo7GVq8QWaxzuBbAJwE4AIwD+xtleaMVPS2wvdcz8jZR+hVK6m1K6u7u7+tW0nPrz7MlJbF8TcTth1oI3X9ADAPjFydb1yY85D7GTY9w4eEnlKYfh2RRkkUARBY9yKJzKqkgCxmJcOSwISukopdSklFoA/gnAHuelswD6PbuuBTDsbF9bYHvOMYQQCUAbKndjcZoY3bTw4tBMzeINjK6QAgCIa607MH7MedC9PsHdSl6YcVAkAcMzKYzMaFjV5kPELyGRtq+HQm4lQgh6wipGSxgHy6L4hyePYzbZmkHrRRkHJ4bAeCcAlsn0MIBbnQykDbADz89RSkcAxAkhe514wu0AHvIcc4fz83sAPElb2TnMcZlJ2imm6zsDNf1/JFGAJJCWrXeglHK3UhGSuglZJFjb4ce5mRTOz2pY3eZH2JdVqoXcSoAdd2BGtxCvjcbxf/7zWE71dStRtvEeIeRbAN4EoIsQchbAXwJ4EyFkJ2z3z2kAHwYASukhQsiDAA4DMAB8lFLK7siPwM588gN41PkCgK8CuJ8QcgK2Yri1Gr8Yp/GJafaKK1JDlxLDJ4st27o7njag6RbCqoShqSQ03SzoKlmJpDL2Z9HX7se56RQm5zLYtb7DzWISCNwJcPn0RlS8dj5e9L1ZDcS4J0MslTER13T0RHzFDmsayhoHSun7Cmz+aon9PwvgswW27wewo8B2DcB7y50Hp/VgN1c9jEMrt9Fg8YY9G6J44ugYBieTuGBVeJnPqjHQdBMBRcTajgAePDkEi1K87eI1nr5KIoolR/aEfXj6WPE4VcxxU4571MW9PzmB771wDs/8ya9U8bdYHniFNKciMoaF46PFV1GLgd1cEV99jEOrNuBjLiWWDvw6dy25JDMm/LKI37lmA37rinW4elMXbrioB2GfvS5WSxRe9kRUxNOGW2WdD1vcTHiUw9mZVMk4RTPBjQOnIu5+/BhuvufpBRdZzSb1nJsn5zXn5mKtCmqJ2sKtu9nK9YqNdjowjztkSekm/IqEjd0hfOqWHfjX374Cu9ZHXeNQyv3WG3ZqHYqks8Ydt6hXOcQ1A4YzYa7Z4caBUxZNN/HAc2egmxQHBqcL7vOJB1/GX//o6Lztn/z+IXz4/gMFj3HdSnVSDq0akGYPp/WdQXSFVJz1dCBd6aQyJvwFspFYQLqkcXDiBsWUgOtWSniNg31NF1MbzQQ3DpyyPPrqCKaddL39HuPw40PncWh4FgDw8tkZHDw3O+/Y0ZiG00XSK+sZkG5l5TAWT0OVBER8EjoCvNW0F1s5zDcAlbqVAGC0SMaS61bKUw6A3a6j2eHGgVOW+58dxMauIC5b147nT9slKIm0gd//1ov456dPAbDVRarADZHMmJicyxRsfRxLGVBEoeQNWi1aOuYQ09ATUUEIQXtAxkyL5t0vBjvmMN9tyZSDWpFbyVYOh4djeNWzAGKLm5hmuNcWMw4prhxan3MzKXzy4UOYSa7MBlwZw8ILZ2bw1otX4/INUbx6bhaabuKJI6PIGJYrn9OGldPHhsFeLxR3iGk6In65aLZINWnpbKV4Gj3Og6zNr2CmysrhxFgC33j2dFXfs15oZZSDr8TCJOKXoEqCW+vwmUcO49OPHHZfj6WyBmByzn4+MLfSXLr5FyLcOJTg4NlZ3PIPz+DrvziNX76+Mou2p5yLflWbD5evj0I3KV4emsEjr4wAADSndkDTzZw+NgxmMAoF9WIpHZE6BKMBu9CpVescxuJpdIdsF0h7QMZslRcy//7iWfzFQ4easnFhKmMiUEAdRCqIORBC0BvxuTGHmZSeM98hns7+PB5Pg1LqcStx5dDSfPqRQ647hD0kVxpsxd8VUrFrvd1f8b5nT+Onr9mND5mcTutWwSCcaxyc1RelFH/3+HGcn9UQ04y6BKMBu0WC1qIBaeZWAoB2v1x15cCMqm42jnH42bHxoskRXpIZo6ByiFQQcwCAaFBx7/25tJGjjmMpw23NMhFPQ9MtGJb9GRVysTYb3DgUgVKKo+fjuHF7LwBgam5ldmccd42Dgo6ggtuvXI8fHjyPjGkhGlSgGRZMiyJjlnYrsVz8s9Mp3P34MTz88jlHOdTHOLSqcjBMCzHNQDRoP6TaAzKSGbOqmVlsOE4jjcz81PcP4YtPnSi7n6ZbBdVBJdlK9n7ZHkyJvJqHmKZjY1cIgH2fMJcSwAPSLY39xzawbXUEIVVyfYorDVZJ2uW4LT59yw48/vHr8PfvuxSXD3QgrZtuTnf+asm0qOt2Ym4llkkzPKPZMQdfndxKcmvGHJgbg3W1bQvYRqKaGUvs79soufuUUgzPaGUTDAzTQsa0ECgVcyiQ5pq/H/uME2kjJ5YQS+nY0BUEYCuHmKexIw9ItzAnx+z0y809oRxpudJgbqVOxzgA9mfy65esgU8Woemme5MaFs1ZXXpjEMytxB5a52ZSdVYOrVnnMJtXK9LufJ7V7BTKjEKjKIfZlI6UbpY19uz68xdUDsytVFo5hFQJCc1AxrCQMSykdBOmRUEpRUwz0BVW0OaX5ysHHpBuXU44Vaabule2cZhMpOGTBQQLrL58kghNt3J8+V7XkleCjztuJdc4TKcQS9Uz5tCabiWWTsmUQ7sz+3i6isYh7RiFTIMYh2Fn1Gc5Y8+UbOFsJeZWKv0IDKmyoxg8qsBJvjAtiohPRndYxYTjaWDwIrgW5uRYAgFFxOo2HzqDiuteWWlMJDLoCqkF0019soCUbrquIyDXtZRMF1cOpyfnkDGtmg758aJKAjKmBctqnKBqNWDplEyBtfttt1I1U6+zyqExPrvzMbsCvJybq5RyUCQBN21fhcsHSk8gZDGHnAd/2nA/97BPRldIwXg83zg0v3Koj8O3CTk5nsCm7hAIIYgGFRwaji33KS0LE4l0jkvJC3MrpY3CaoHdIGGfNC/mwF6rZyorYK9+fULrtLPOVpnbnyNTDpVmLL1ydgaJtIGrNnUV3afR3EpZ5VChcSigHADgy7ftKvt/MffTqGeWdCJtuFlJEb+E7rAPB8/O5LiVWsE4cOVQhJNjCWzqtoNN0ZDtVmrGPO+lMpHIoNtJ18uHtaTIUQuen1O6bSgGOoOYSKRhWXReoLRebiWWsthqrqVs88Jct1IlMYe0YeJ37z+AP/v3V0vu12gB6ZFZWzmU+1smS7iVKiWkSs7/mTUOyYyZ0xdsTbsPwzOa+7cIq1KOG6pZ4cahAHNpA8OzGjZ122lqnUEFGdNyU9pWEhOJNDqDxZSDffnkZGl4gtAsKDfQFYRhUUwlM/ONQ53cSixlsdVqHfKbF4ZUCaJAMJMq71b63oFzGJ7VcHY6CbOEuy3TYDEH9qAuF3PQMsXdSpXCYhOjHuMwlzZy+oL1dwSQMS28Pj4HQoDusIpkC7Rq4cahAEed6U+be2zjEHUejistKG1ZFFNzGXSFCysHn+Oq8fq3C6mIAWcM6GjMXl15UwvrlsraosohpukQBeJ+poQQuxCujHLQTQtffOoEBGLHEkrNIHDdSo2iHBy3Uikl89LQjLuYK5TKWikhXzHlwGaRSOiP2tf3oZFZhBQJQVVCsgUWkivOOGi6WdJ3OjSVxMceeBFhVcJlTkVwp1NgVOtah0YzPjMpHaZFSygH+6bzqoFC2UoDnbZ7biyeRiylY3NPCIIT365bKqujclotnTWWMtCW15+qLVC+SvqXr0/i3EwK779iHQD7ui9GowWkXbdSEeMwNJXEO774DL79/BCApSkH5lbyGs+5jOHGF8I+Gf0dfgDAsfMJhH0SAorIi+CakW89dwbb/+LHeOs9T+Mnr43Ne/2Pv/cKYikd//rbV7j93Fn16VQNM5ZOjCWw6389hpeHZmr2fywUNtinK1zareT1b7M4A5A1FBuc2M3orK0cokEFq5zPtn4xB/sB0WqFcLOp+YWE7X65bMzhwOA0BALcerljHErMgGikCmlKqbuKNywKo8A5sar+Z1+fBFC+CroUYVc5ZD+fZNp0Xalhn4S+Dj8IsT+nsE9GQBF5+4xm5JL+dnzw6gGcmUziP148N+/1oekkrt/Wi0v6291trnGo4cp+aDoJSoHBEiu4euNtnVEIdtPNlFEOG7uCIAQYdoxDm19Gn7Paql+2kn2pt1rbbtbZ1kt7QCkbczgwOI2tvWFs7Q2DEOBMBcqhEQzrdFJH2rDcxUWhOEjCeXCza3EpbiVmHM7P5iqHWEqHKgnwySJUSXTbe4d9EgKqxBvvNSOXrevAXTdvw/a+SMHVUlwz3AuC0RmqvVuJBRar3VFzKUzktc7IhykHr3+7UMwh7JPRG/ZheCblGoc17X74ZKFshWq1YIasER5w1SSW0uepr3IxB8uieOnMDC5b3wFFErA64sPZEsYh3UCprMMz9j273oljFYo75CeOVCNbaSyehuT4QpMZc55R7o/ai52wT0JQEXNqfJqVFWccGP0dgXl+VtZyN984BBQJPlmoafM9ZhwaaVDLpKcjayFYQHrWs0rNSWXNmFAlAaJAsKbdh7PTScQc43DLzjW4/cqB2p18Hm5AusViDszYemkLlHYrHR9LIJ42cNk6O6a2NhrA0HQp5WB/Zo1gHNgKnsWxChl7b70BkL1OF0NQkUCI7cJqD8iQBGJnK6VynxP9Hbaxst32mvZEAAAgAElEQVRKEq+Qbmb6owGMxdM5boZkJlsSn09nUC2oHP71l4O47av7lnw+zIdZ7XbLi+XloRl849lBhFTJ7deTj5oXkBYFkpvKmjEQdFZea9r9OD6agEXtOMOvXNiLP715W41/C++5tmq2kjHPNdfuVxBPG0Uf5i+csVtdsxbs9kIphe/sH8LnHp0/B7yRYg5shkKv06K80N/TW6nskwUIwuKHSQkCQUixP9+QamcizaUNJ8U7625dG2XGwQ5I8yK4JobJQO8w9pgnAyGfYv2Vnj89hZ+fmCgYGFsIsw2kHDTdxG/98z5ouol/un130ZvL61YSBYKQKs1zK7FMkb52v2tc69UywwtzX7VincM8txIrhCuy0HhhcBrRoOKmGPdH/RiNa/jMI4fx7efPzNvfLYJrgGylVMY+F+bSKaQEmVtpVcS3pEwlBktnDaq2y2guY2LcM30PgJuxFPbJCKoSDIs2TNHgYlm5xsGRgV45zVYchYKkXSEFpybm5hmByUQGlNrB29mkjgefH8KXf3oS52aKZ38Uwo05VFC8VGtmUzriaQO/9yubceWmzqL7eQPSPkmAXxZz22ekTQRVe5817X53e73SV734WlA5aE5n0vzPs8NZ0U4XiZGdGE/gwlVhN/21vyMASm0VwtKXGYZpgf2zER52TJm2O63JC7mVEpqBoCLi8g1RdAQKJ1MsBBZ3CKl2sDmZMezpe54svn6PcmAGqdldSyvXODh/TG8gLl5CObx3dz8GJ5P4l2dO52xnLa3Pz2r4ytMn8T++9wo+9+hRPPDc/BVYKRpJOZT6HLy4dQ5JHT5ZnCenk7oJv5J1KzGWUzm0UkDaW6XrpcsxDhNFUq+n5zI5/bLYvSAJBJTmFjV6s4Eawa3E5iS0u8qhsFsp5JPwl79+Eb5y++4l/59MOYQc5TARzyCRNtzpe4CdkScQoCesuguiZq91WLHGoTukQpGEnIylbKfF+crh13aswlu29eJvHnsNZz1qg92Ao7E0Tk8ksb4zsKg8Z3ajN0LMwZvDXQo2nD1jWlAlAf683zuZNtz5vWvasxJ8eYxD6wWkvVW6XtiDf7JIAsXkXCbHX759TQTXX9iD371uE4DclG2vWmiECumUbkIUiBvLKuZWCqkSukKq2+VgKbBFUsgnIaBIOD1pz3rxupV6Ij48/HtvxNt3rnEXRM0+8GfFGgdBIFjb4c/JWHJXYgVWzIQQ/PFNF0DTLfzihF1cY1rUzWAajdk9atZFA85gmYXdSLOuW2n5jYPrXitnHDz+XK9yODA4hf2np5DMZN1KfV7lEFg+46C1kFupmHJgqdeFYmS6aSGuGTnulqAq4asfuNx1IU4WMw4NoRwsBGTRU/FeQDmkDYSqWFwZVj0xB1V028935xWH7uhrgyqJ7uyTZh/4s6Jbdvd35Kbwxco8FJn8ZsVhM8mM64+1jUMKv7qmDcdG4wv2z7JV4GxSB6W04PyEesHcSiG1MrcSYGcu+RUJsZSOT33/MHSTQvO4ldr8sms8lkM5SKIASSAtphxyO7IyOgIKCCnsVpp2XEbRAoWNhYo9vQ/fhghI6yZ8iuga+0L3WbzK42e9MYeEkn3fniKdAwLOPs1eCLdilQNgZ2kMTWXdSuV87T5ZdGYT2LnW3pvv1MQcJucy6I/67WH2C3wIsVVgxrRy0kGXg3iFbiVRIJBF24j5ZAF+WUAqY+LU+BxOjMURS+muW4kQgjXtftslsISipKWgSkJLBaTzR4QyRIGgI6BgMpHGgcFpXPKp/3SvWfbgjxYI1HYWMA6NFnPQdDsDrlQMKaEZ7gO9GuTEHNRKjIN9bs3eQmNlG4eOAGZTuvtgjmsGZJGUHB3YE1ZdWcmKxAQC7B+0c8fXdgSgOFPHKsWyKGIp3S02W+6gdNZIlr/BWIGRKgkIKBLOzaScHHuKybkMAmrWEKxp989rEldPVFlsqVTWWInsOja98MUz05hN6Tg0Yg+rco1DcL5xYBlAxWIODZGtlGHGgWWfFY85VItwXiorYAfvi2VC8YB0C8CyDdgI0FhKR9hX+uHVE/a5xoG5lzZ2hzDubFvb4bdHUi7gRprLGLBotiXA8hsHA4TY1aHlYIVwPlmEXxHntS7w9rW5dksXrtlSfOJYrWk15ZA/y8FLpzOgitXxnJ6wg6iljIMiCQj7pBzjoDeYcsh3KxVVDjVwK4WdVFbAjjcUq/9hbqVmb9td1jgQQr5GCBkjhLzq2RYlhDxGCDnufO/wvHYXIeQEIeQ1QsiNnu27CCEHndfuIc4TmBCiEkK+7WzfRwgZqO6vWBzmq2XyPK4ZZX2VPREVY87IQGZUdqyJuK/3O8phIQFp9v+vc2IalQxqqSVxR5ZXUlnqV+xLyCeJrgvJS8BjYH77mo34u1svrd6JLpDFJAo0Mok0U7rzP/fOkIqJubSbWTc4aX+fLmEcAEdxFFMOjWAcMib8np5c+Yswy6JIZIyyadgLoZByKOZSArILomavkq5EOXwdwE152/4EwBOU0i0AnnD+DULIRQBuBbDdOeZLhBB25d4L4E4AW5wv9p4fAjBNKd0M4G4An1/sL7NQ2vKGscc0vexF1RNWMR5Pg1KKiYTdjGtLbxiA/fDpCilQxIU9hFgwmgW8KxnxWEti2vyq22Iwt5JPFtybQhQINjptupfSEbPa+OSFx4IambRuFe0bxNxKrnKYZMrBvrbai2SMRYNKTvFcbrZSYwSk/SWyleYyBijNZhhVA/ZMCKqiRzn4iu4fUiWokrDgQthGo6xxoJT+DMBU3uZbANzn/HwfgHd4tj9AKU1TSk8BOAFgDyFkNYAIpfRZag9i/kbeMey9vgvgelInp3Qh5VDOz94dVqHpFuJpA5OJDKJBBavb7AtlbYcfhBB3tnKlsP9/vasclt+tVEm8AchmLKmSCJ9jCNZ2+PGGvjYAjWUcVEloqVRWzTBdt14+nUEVsyndVQzs+9RcGhGfBFksfOtH83qIpc3Gq3MIKBIUsXDdCnNrVtOtxAo4V7f5s8ohUlw5SKKAHX1tDTWbZTEsNubQSykdAQDne4+zvQ/AkGe/s862Pufn/O05x1BKDQCzAAr2bCCE3EkI2U8I2T8+Pr7IU8+S34MmXsGKmRW+jMXSmEik0RVS3d7ya52WHIq4sJgDC4iva5iYg74A4yC435lbaaAziAtW2WoqUEHcol4sJouskdF00/W958NqHVK6iZAqYWgqCcO0MJXUc6qj84kG5Zzuw+w6FkjjuJV8sgjByZTLX4SxTLtqBqR39rfj6f/xZlywKuxez6XcSgBwydp2vDo82xBxmsVS7YB0oRU/LbG91DHzN1L6FUrpbkrp7u7u7kWeYhamHNjDOL8NbyHYRTEW1zAxl0FnSEGPYxxYMz9VFhb0EGLGqTfsgyIK82IOvzgxUde0OFs5VOhW8gSk2Y2zoSuIbavsOExQbSDlILdWzCGtW0Uz67wDmvZujMKwKIZnNEzNpdFRoggxGlQxNZeBLfCzxiGoSg3xoLNrZ+zfWZXEeYuwStOwFwpz+bLrOb8ALp+d69qh6RaOjcareh71ZLHGYdRxFcH5zuZtngXQ79lvLYBhZ/vaAttzjiGESADaMN+NVRNkUUBQEXOVQ5kCLSYnx+NpTDrKYU27Dz5ZwFYWe1iocvAUM7UFZPffAHBsNI73//M+PPzy/Kl1tWIhbiUWGFSdbCUAGOgM4I1buvCnN1+IqzYtX3ZSPqokNk220mxSd/t2FSNtmEVHYHrVwdWb7b/Bqck5TM3piBaZCQ7YsQrdpIg77hl2HYdVqTFSWfVsp187waCwW6naxoHBAvnMS1CMnWvtSZIvNbFrabHG4WEAdzg/3wHgIc/2W50MpA2wA8/POa6nOCFkrxNPuD3vGPZe7wHwJGXLljrQHlAwk9RhmBbmMmYFMYd8t5KCgCLh8Y9fh/ftsefxLnSFGnNSR8M+ad4Ur6ePTwCo7RS6fBbjVlIlwV1Vre8KQhYF3HntpiXN7602qiw0TZ3DJ79/CB/51wMl99F0q+jn681GYunDg5NzmJ7LIBospRxyO7oyV1IjKAdK6Xzj4Bh706K4/5eDmHBSystV9y+W7Wva8O079+LaMinZ/VE/okGlqeMOZZ8AhJBvAXgTgC5CyFkAfwngcwAeJIR8CMAZAO8FAErpIULIgwAOAzAAfJRSyu7Gj8DOfPIDeNT5AoCvArifEHICtmK4tSq/WYVE/DJmUxnPiqP0RRXxSVAkAacm56DplrtC864kFhxzSOlu6mh7INc4/OKEbRyq3XPp8z86iqGpJP7h/ZflbM9Ow1u4W2nvxk78wfVbcFWJNt/LSTPVOYzMpjA8o5XcR9PN4m4lRx2EfRI2dYfgl0WcmpjD1FympHJgxmFyLoP1ncE8t9LyZiulDQuUwk188CZ+vHhmGn/+H69ipzP7vZoB6Xyu2Fj++iaE4JK1bXjhzAxMi0JcwsCh5aLsJ0gpfV+Rl64vsv9nAXy2wPb9AHYU2K7BMS7LQbtftucXVNhsjhCCnrCKR162vWIXOK4kL+oCUya9A1va/Iqbm66bFn75+qSzT3ULag4NxwquajTdgmHRxQWkFQkfv2FrVc+zmtgB6eYwDom0keNeLIRmmEXdoBG/BEkgWNsRACF2avHTxyeQMa2KlMOUU8PDjENIldxCz+WCxd2YcvAuwpiyZm6cWrmVFsI1W7rx6UcO46Yv/Az3/tZl2Nwz/1nRyKzoCmnA9vPPJHV3ZV7JirknrCKmGbh8oANvumB+YHyhysE7B7g7rLg9m14emnFL8GNadZWDppuYTenzBsJUOsuB4dY5LGFOb70o5KNuVBKagUTGgGUVX61rJQLShBB0hhSsdSaU3XHlAE6MJQCg5ACc/OZ7zK0UagC3Eus5FnCVQ/bv6Z1BAVRW3V9rPnDVAP7h/ZdicDKJbz03VP6ABmPFG4f2QJ5yKNCnJh+Wzvrnb7uoYKsNRRJgUVQ8OtQbAO4KqZiaS8O0KJ45MQlC7EEi5VaRC4X1pDnlFEcB9uSqcp1p83HrHEr0o2oUfHLzBKTjml3MFS/RgiFtmCWN8qfevgO/9+bNAIB371qLrb32bIPOAh1ZGSwLh7WGYUoroIjLnsrKjIMvJyBtn5PXFRtUxIZw4wgCwdsuXoPeNrVg+/RGp/Hv6BrT5pcx42m+V0ll8H954wb81bvegIudjIR8SvV9KcRcxnC7PXaFVFjUbq18eGQWG7uC6I8Gqm4cWDEY67mz7/VJ7Pz0Y/i/++wJdgt2KzWBcgj7JGRMC9oyd72tBGYUSv3dNd0qWgQHADftWIVLHB+8KBD82VsvgioJ2NhVfAAO6zzMXEgZw4IiCVAkYfmVQ55byesmnE7qYPaglvGGxZBfWNgscOMQkJExLLw+zqY7lc5fBoA9G6JuZlIhlBK95guRyphuGijrzDoeT+P8rIY17X5E/LK7oq8WTI6fnpjD0FQSH/nmC8gYFr7/ih1LWUxAutHpdMdnLq/vvBxpw3SvnVLuxFJFcIW4bms3Xv3UjRjoCpbcz+48bAfDM4YFVRQgi8KyBaRPjCXwiQdfdtU9u1cUT4PLmaQ9+vSi1ZGq9lWqBp1BJaewsFloLBO7DLQ7/ZX2n55Cm18uW9xSCQudV5zMmG51MStemkikcT6mYWtvGKosVD1biSmHU5NJfO7Ro9ANC7vWd+CA03q84joH1zg0/jqDGd7JRKZsnvpykvAsBEolIqRLpLIWo1jbDC89YR/GYo5yME0okmMclimY//DLw/jeC2ex3WlwWajOYTqZQUdAxv982zbXiDQK0aCCI07L9GZixRsHFgjePziNrb2hqswaWKhySGYMN8jGjNP5WQ3j8TRWtflgOvMeqjkhjuX7Hx2J4ex0Cu/e1YcLV0U8xqHSgHS2WrXRYb72YrOVGwVv2/NiysG0KDJm8YD0UugOq27WD3MryRLJ6bNUTw4PzwKw6zSA4jGHdr/SUEWXDNbpdrknPC6UFW8cvP2VthRIS10MCx1mn9JNt9tjl2McjozEYVFgVZsPcc2AYVG36Vg1YH73404Gy43bV2F1W3bOc6XKYXWbPd2tK1w8yNkoMOUwEW9s/288RzkUNg5s4VELd5638zAzDqpoxxyW4wF3aNhedZ92mgcyt5K34n0mqbvzUBqNaFBBxrCLbKvZ86nWNL4voMZ45+9u7SkeqFsIygIC0rppQTep61YKq3aR3avOamlVxDeve+xSoZRC0y3XAER8EvZu7MSm7iC6QioIAUIVGqGrN3fi2T/5lRzD0qgw5TDRVMqhsIuEGXffAmIOldITUZHSTSTSBjKmBcWJOVBqK5Z6MjWXwcisHf9gyiHgjTmYLCCdKZmiu5zk1440C9w4eI1DlZTDQowDGwjCVkOEEHSHVBx2Vkur2nxuBlW1CuHYebHmeNdv64UsCiCE4IqNUUR8ckWDftj5ssaDjU5AkRBQxIZQDppu4t9eOItCnWLylcMvTk64WWTu8Y4qLZWttFjczsPxtMetZF/T9Q5Ks/sAgDubIrd9hglKqe1WKlHct5y4sa4GX5Tkw42Dp0Nltd1KlcQckhn7QeB1F3WFFHf1WAvlwIzDpevbEVREvPuybE/EP77xQtzzvuWb1lZrukJqQ9ykjx0exccffNl1mXhJpLN/55im4+vPnMb//8MjOYaEJRTUKuYA2P3D0sw4OIFstlKfTep1KSg85CjozqACw1Et3tqatGEhpZvImFbDK4dJrhyai7AqQRQIOgJyTpvjpeAahwoCeEw5eIfisJtTEQVEg4pbmFetWgdWALcuGsCrn7oRb/Q0EVvXGcB1W5feDr1R6QwpDXGTsopethr2wrKVVElALGXg3EwKibSB87FsryX2YK5FfYm3LX3GsN1KimgrSbbgeee9z+CLT52s+v+dz6HhGPra/djsuHwJyd5fqiTCsKj792wv01F5ucivOm8WVrxxIIQg4pOwpSdctUCbm8paQbFVqoBxYDK0t011zs9xK1WphQZbdaqS2FTZE9WgM6hiIpHG6Yk53PuTkwXdOvWAxRIKjZJkr/W1+xHTdAw7+7D2F4BXOdTOrTQeT9sxB6cIDoBbCDc8k8JIHcZgvjo8i4vWRFzXpV/OXrPsnEYdo9neoMohmyXHjUPTcf22Xtz8hlVVez9lUcrB61ayjcPqiB3krbZbifmrm6E2odqw3lXf3DeIz//oKJ51GhvWG6YChws8YBNpA7JI0BVSMRrTMO20hsg1Dk7MoQYB6YjfTooYd2IOqsetpJsWLMtOaEjVuNL8yEgMr4/PYe/GTvQ6asbvMYbsd2eKqtQQo+UkoEjwyULTFcKtvKdDAf7Pey/BB67eULX3c1NZK+jjw2IO/hzlYK80VjmzqVlWUbUC0tlMl8avTag2nUG7dxXL4/9mXqC3XjAVWNA4aAZCqoSIX8Zr57OTxAoahxooB5YUkROQ9hgHFrOqdRuS/7vvDBRJwLsv63OHbPlyjIP986hTsNcRbEzlANjXHVcOnJLK4fnTU9j1mcfc3jUF3UrOKokZB8mZWFdtt1IztLyoNp0hBRYFXjwzA1kk+PGr55elFTVTgYXcSom0PU8j4pfcB7EiCUXcSrW5hXsidgsNbyorAGQM6i5oaqkckhkD//HiObz1DavRHlDQ67iVvPeJOs+t1JjKAbDjDjzmwCkZc3j88Cgm5zJuOX3BgLTjVlrlSRFtc+ZOVANXOaxAtxJz2RkWxZ3XboRhUXzvhbN1Pw+mAgsph7hmD3/yNoHcMxDFyfGscXAD0jUy8KwQzi2C8yx4mFFI1nCu+Q8Pnkc8beD9V9g9zFiShldhz4s5+BtXOUSDjZEIsRBW3tOhDpRSDi+csdtTDE7Z1Z5JPbfOAQA2dAfhkwXs6Gtzt0WcuRNffOqE20l1saRrWF3b6HjbVb/z0j50hVS3uKqeMBU4kcjMc8/ENQMhn+QO8hEFgqs2d2IikXGznNI1Vn+rIj6cm04hmTHnuZXY+aZqaBx+cXICnUEFu9d3AICrHHyFYg6zGoKK6N53jUhniCsHDuwUVGB+zCFjWHj5rJ23PcSMg1PP4B1O0hP24dVP3og9G6LutohfxjMnJvC/f/waHnam0C2WWgYzGx2mykKqhI1dIfgVwXXR1JNYSndnDuSrh0TaQFiV3JkaqyI+XLjKrsFhriW3CK5Gf8OrN3dhLmMPhFJEEbKTyqobFlKZ2sccDgxOY9f6DjczqadAQJq1ojh6Pt6wmUoMu78SD0iveGSRgJD5yuHwSMzNEz/j9IlJ5vWoZ0h53TMjPtmV89PJpa1ANL22LolGhs383tEXgSAQ+CSxpivgYsQ0A5u67dbZ+bOi2fAnphzWtPuwqdvO82et5Wv9N7x2azfCzsPXWyGdMa2axxzG42kMTiaxe6DD3RZS7ep2731y+YYo3rKtF7MpHR0NWh3NiAZVaHr2s2sGuHGoAYSQgqNCWcfT7Wsirlsp5QyJL9euwjuhLn+050LRnPNqhult1abdLyPsk7BnwFZlfkV0V+H1glK7y+6FTvuSczPJnNcTacet5GPGwe+6VdicBTcgXSPl4JNF3HBRLwDbOCiuW4nWPOZwYHAKALBrfVY5E0Jw4aow+qPZHl6yKODe37oMt+1djxsvql4qei1gabbN5FpqnhaBTYa3nTDjhcFp9LX7cflAFN89YPfVsdt1l/8zbFsVwabuGSiSiKnk0gLT6RWsHASB4Id/cI0b4PRJYt0nw6V0E4ZFsbU3BIEA5/KUg53KKrsLgr52vz2hTZXc+eJpw4QkkHkKs5q89eLV+LcXz0H1FMFlDAtMENfqc9t/ehqKJGBHXyRn+7fu3Asxr2hTFgV85h07anIe1YSl2c4kdaztKLNzg7Dylo51QvGMMGS8OjyLnf3tWBcNIJE2MJ3UkcyY81xKhfidazfisf92HVZF1KUrhxVc5wAA/dFAdoKdIiJVh5jDFx4/hqdeGwOQzVSKBlX0OoFfhub0CQrnKQfAaf3h/O013ap5zOiaLd148wXduGxdR05AmikH3aQ1GR26f3Aal6xtmzcjRJXEmhrDWsL6Pi3VJVxPmvOTbgK8U6oYU4kMeiIq1kXtvvODk3NIZcycNNZSCAJBRxXypdOGBYHADTKuZHxOZ89aYpgW/uHJE7jreweh6WZ2XrlfwoWrwth3ahKW01SONVwM+yRs6Q3hfXv6cf22HgB2vGTCqcnQdLPmyk+RBPzLB/fgyk2d7rVixxyyn1e11YNhWjg0PItL1zXJ8rpCmFtpeomqv55w41AjVCk35mCYFuJpA21+2R1KcmYqaY8IXcAAkGhAqUpAeiX2VSqETxZr3gZiNJ6GYVGcj2n41nNn3NYZEZ+Mt+9cg7PTKewfnMapiTn84qTdziOkSlAlEX/1rovdWRldoWzGi7aIEaFLQfEqB49xqHYw/9xMCrpJsbm7OrNVGgWWTTXTRMqBxxxqhJIXc2DN1Nr9sju/+Mxk0o45LOAm7wgqSGbMJa0c7QcLXxcAdpZYrWMOZ53kgza/jC8+dRKfuWU7ADs9efdABwLKq/jyT0/ipaEZVxUWmhjWGVKx/7Sd1JA2zLomFLiN9wwr5/OqtmE95dTwDHQFq/q+yw2r3p6e48phxZOvHFh1c1tAhl8R0RNWMciUQ4VuJSDb/ncp6qEeLolmwScLVV/9nhhL4P5nT7v/Zm25P3zdRkwk0th3ys7GifgkBBQJN+1YhSePjkE3LHzmHTvwgasGsHdT57z37QoqmEpmYDqN7+oZM5ILZCsB1TcOg06K90BXY478XCyyKCCsSjzmwGHKIXvjMDnJSvwHOoM4M5lEKmPmVEeXgwW2lhJ30Iz6uiQaGZ8iuqm91eI7B4bw5w8dwqzjX2bG4S3b7NTQ5087xsGpY/jNK9ajPSDj7t/Yidv2rscn3749p3UGoyusglL7b19v5eAd9pOsoVvp1MQcgoroFiu2Eh3BpbuE6wk3DjVClcSCyoE9EAa6Anh9Ym7xymEJ8tSOOfA/PWBnbGUMyw0IL4RHD47gC48fm7edZSMdG7M7qp6dTqI3omJzdwh+WXT7arFuu7vWd+DA/7wBb3HqCorRGcyOm9R0s87KITvsJ8etVGXjcHpyDus7gy0ZD+sIyDwgzckdfg5kjQPzPW7oCmEikcbUXKaiOgdG1KkEneJuparAVNtiCuH+/cVz+OrPT83bzrKRjo0y45DC2o4ABIFg66owLGq7s7ypmmIFM7vdoTGJTN3jRoQQqJIATTdzA9I1cCttaLF4A6M9oDRVQJobhxphDz8vEHPwM+Ng3wAZ01qQcnDzpZfgVkobPCDNYBXGi+mvNJFII64Z81qpx53kg2POLIazM0ms7bAzjrY5PZLaFjHSks35mEikkTbqb+DbAzKmkxmkdLsAD6iucTBMC0NTyZaLNzA6nM+vWeBPiBqRrxxmkoWNA4AFGYc2vwxClhZzSDuprJyscljMQ45VK+c3zmOpqsdGEzBMCyMzmmscLnCMQ6GYQjlYu/EJVznU92/YEVDcwk1W8VtNt9LZ6RQMi2J9Zwsrh5WSrUQIOU0IOUgIeYkQst/ZFiWEPEYIOe587/Dsfxch5AQh5DVCyI2e7buc9zlBCLmHtIDDsZByCCqiG9hb3xkA+y39C3ArSaKANv/SViA8lTULe8AuJp11MmHXHMwzDh63EqtxYOnLrnFYhHKI+GRIAsFkIr0scaP2gIyZpN1iPOoo2Goqh9NO6/RWdSt1BBTE00ZNqsprQTWurjdTSndSSnc7//4TAE9QSrcAeML5NwghFwG4FcB2ADcB+BIhhC197gVwJ4AtztdNVTivZaWQcvC6EnyyiDVOcdNClANgF8ItLVuJxxwYTEEtdAWcypiYc47J740U1wwIxB4o/7IzjpQpB9Zsj7XjXgiCQNyhMcsRN4oGbeWQ0k23C2o1laQ6GyMAABYeSURBVMNJp+PsQIsqBxYvbBbXUi2WHrcAuM/5+T4A7/Bsf4BSmqaUngJwAsAeQshqABFK6bOUUgrgG55jmhZFFHPaMsymdLTl9ZxnK6SFGoelpsTVO9OlkWFupfxWJ+WYSGR78xdyKzEj8MODIwDgKodoUEFfu991ES2UrpBqZysZVt276rKAaipjurGvaimHM5NJfPGpE9jSE3JjK61Gtkq6OVxLS726KID/JIQcIITc6WzrpZSOAIDzvcfZ3gdgyHPsWWdbn/Nz/vamRpXzs5UyaPPnrhaZcaik8Z6XjoCCqSWlsnK3EmOxAenxIsYhbZhIG5Y7i+CRV0awtsOPvvZsq+n7/sse/NFNFyzqfDtDiju+s94GnqViJjMmgqoERRKqYhwsi+LO+/fDtCj+8bZdLZnGClQnmaSeLPUJcTWl9DIAvwbgo4SQa0vsW+gvTktsn/8GhNxJCNlPCNk/Pj6+8LOtI4ooQDepmz8/m9LnzbhlLQKCC+itBNjydGnZStytxHAD0mXcIzFNx3/95gGcn7VdSKwBXliVcowDy1Ta3BPCzW9YhQ9cNYAf/P41OSMsN/eE0BP2YTF0hVQcGrbrJIJq/QPSpkUxkUgjoIgIKNUZlHRuJoWj5+P4xK9uxcYW66nkpb3Jmu8tyThQSoed72MA/h3AHgCjjqsIzvcxZ/ezAPo9h68FMOxsX1tge6H/7yuU0t2U0t3d3d1LOfWawx46X3vmFNKGPW4xP31x+xrb9cBGIFZKNKhiai6zqMAWpbQu7Z6bBTcgXcat9OrZWfzw4Hl83xnRyjKVdvS15Uxyi3lSlr/0m7vwybdvR1ugelPKrtzYiXWdAXz42o14z67+8gdUEeYW0U0Kv2xPZauGcWCB6C094SW/VyOTnenQ4sqBEBIkhITZzwB+FcCrAB4GcIez2x0AHnJ+fhjArYQQlRCyAXbg+TnH9RQnhOx1spRu9xzTtLzrsj5csSGK//WDI7j7seOYSeruyoGxd2MnfvpHb8KW3oXdFJesbUPGtNxg50JIu1PguHIAsi69cg85Vqfy8xMTALIxh4vXtuF8TIPhGGrWYDG8iIBzJfx/l/fjyU+8CXfdvM2tlq8XHYHchAp/lTranm7Rfkr5NFvb7qUsH3sB/JwQ8jKA5wD8gFL6IwCfA3ADIeQ4gBucf4NSegjAgwAOA/gRgI9SStmV9REA/ww7SH0SwKNLOK+GoCfsw7c/fCWu3tyJHx86j7RhFUxfXExO91WbuiAQ4GfHJxZ0nKabbjCMu5VsWFC3XH8lZhyeOzWFjGFhIpFGm1/GQFcQpkUx5riZvO24W412T0KFXxHhq1JH28GJOfhkAb2LdLU1C35ZhCoJODIScxcTjcyilzeU0tcBXFJg+ySA64sc81kAny2wfT+Axp/1twj2bujE3zxm99/JVw6LpS0g45L+djx9fBwfv2FryX1nkzpEkSCkSviDb72IE2MJAOABaQemHLQKYg6AnZ3z4plpTCTS6Awp7pS24ZkU1rT73ZjDYuoYGh2vcmAxh2rMkT49mcT6aLDsHPVmhxCCm3aswkMvDePEWAL//tGrGroYlT8hasyeDdkh6YtpmVCMa7Z04+WhGbfzZzHuvH8//uI/XgUAvD4xh9edfvk8ldWm0iK42ZQOgQACAZ45MYGJRAZdIRV97fZqd2jado0wI1Irt9Jy0uFRDj5ZhF+plltpzh2A1ep84Td24o9uvACHR2IYmkqVP2AZ4cahxlzS3+5O0crPVloK127pgkWBZ06Wdi0NTiZdgzDpSb/kbiUbWRQgCqRsQHo2paM9oOCS/nY8cXQME4k0ukMq1kWD6Awq+OHB8wBa260U8ctgi3u/bLuVlhqQNi2KMy3cbC8fQgje0NcGoPGL4bhxqDE+WcTO/nYA1VUOl/S3oyuk4Ms/PVnUf0kpxVQyg/F4GoZpYTqpI+ykzfJspSx21k25mIOBiE/Cuy5bi0PDMZyamENXSIEiCXjfnnV44sgohqaSiGsGRIEsuLCxGRAF4l7DLFtpMTGH6bkMDg/HQKk9OjVjWi3bT6kQLJFgqbPgaw1/QtQB5lqqpnGQRQGfevsOvHJ2Fv/4s9cL7pPMmMgYFsad1uAA8NvXbMQ7dq7BznXtVTuXZscnCxUphza/jPfuWovOoAJKs43w3n/FOhBC8M19ZxDTdER8UssXcvkV2zgsJubwZ/9xEDff8zTe9vc/xyNOanCrZyp56Qg2RzEcNw514Df3rsPHrt+C/qi//M4L4K0Xr8bNb1iFv3v8eMEVHDMIGcNyZ/Nu7Q3hC7deuuj2Da2ITxbLBqRnUzoifhk+WcTtVw4AsCezAcCadj9u2NaL7+wfwkxSR7gFXUoMllThVxYXc9B0E08dHcfu9R2Ymsvgrx49CqB1+ykVgjUtbPSUVm4c6sDqNj/+2w1ba7KavG5rNzKmhckCqxCvbD3qzBboWmDB3UrAJ4tllUPMU8R4x1XrceP2Xly5MTvn+Ve29WByLoMXh6YR8bdeMJrhKgcnIL1Qt9IvTk4gpZv4/eu34Fu/sxddIQU+WcCqSGunsXqx04CFho85tO5VvELw9mvx9u8BcqfFHT1vt1zorHPhVDNQSaWv1zi0BxT84227c17fM2C7DoemUujf1LouknavcZBF6CaFblpuK/pyPHZ4FCFVwt6NUaiSiO/87lU4PTHX8mms+XQssbNyPeDGockpFNz64L88h13rO9wcfAA4MmIrh07uTpqHTxZKNt6jlBZsf+JlfWcAXSEVE4l0S6axMjq8biVPGnAlxsGyKB4/Mobrtna7+f0buoIrJlPJS0dA4TEHTm1xg1uOSohpOn5ybBzPvj6ZYzCOjcahiMKi5gi0Or4ybSCSGROGRUsWthFCcLnTibUV01gZbHERcGIOQOUzHV44M43xeBo3XNRbs/NrFqJLbLtfD7hxaHJYcIsZgoNnZ0Gp7d6YmstAEghEgSCZMdEZUlo2i2YplGsDkT//uxiXO66lVqyOZvzG5f348m9dhrBPdj+PSgOr33vhHPyyyI0D2EyWhQekKaV4yqmzqTXcODQ5rDCJSdSXnGZ8wzMpTCTSiAYVd3hKZ4sOUVkqPll0GxIWYsHGoYWVQzSo4KYdqwFkU3knK3hQabqJH7wyjJt2rFpwi/pWpCMgLyrmcD6m4YNff94dIlVLuHFockSBoD2guMFnZhwMi+LISNwxDvZN3Bnk8YZC+GWhpGukUuOwbXUYb79kDd64pbPkfq0CW3SMV2Acnjw6hphm4F2XNf0cr6rQEVAwm9LLNuA7MDiNJ4+Ouv9+zck63LrATs6LgRuHFqAjIGN6TgelFC8NzbjzIY6ej6EjoKDb+TevbShMuVTWSo2DJAq4532XYtf6aMn9WoWscii/Av7BwRH0hFVctamr1qfVFLBEEnZtFePvnzyO//6dV2BPUAaOj9qNM7lx4FRENGinxY3MahiPp3HzG2zZr5sU0aCC7hAzDtytVIhyqayVGoeVRptfhigQTM6VVw5HhmO4bF0HxBWWslqM/ESSYozH7e4GJ8dto/DaaBzdYbUuszy4cWgBOgJ25gMb/vPWi1eDxZ2jQcUtfOMxh8KoTsyBjXTNp5Wb6S0FQSCIBhVMxEs/4DTdxOnJOWxd1dqT3hYCSyQ5MDiN377vecS1wgqCqbJ9p6YA2FmHF9RBNQDcOLQETDkcOR+HQIA39LW5g1M6PMqBxxwKw/L1iwWlYykdhLRmG+6l0hVSyyqHk+MJWNRu3cKxYW1IvvSTk3j8yBiecx7+Xiil7mf73KkpWBbFsdF4XVxKADcOLUFHUMFMUsfx0TjWdwbhk0Ws7bAL4KIBORtz4K0zCsIGHxVLZ51N2d1sV1oVbyV0hRR3nnYx6uknbxaYW2jQGZF68NwsgOw8EPtnA7ppq9l9r09haDoJTbdwwar6GFluHFqAaEBBxrTw4pkZbOmxLxxmHDqCCi5b34Hd6zuwY01kOU+zYWEtSEZmtYKvz6Z0tFVpil+rwarCS3FsNA5JICuquV45vIOTAODVc7PY9/okdn7qP92MJPa5XrauHedjGp44MgYAC545v1i4cWgBWHDrfEzDFke690ft/j6dQRV97X589yNX8dYZRbhio51d9MyJwoOTyrXOWMl0BpWy2UrHRuPY2B2EwmeIuLD2IwIBrt3ajYPnZvHIKyOwqG0ogGy84e2XrAEhwF//2O5gyxaAtYb/tVqAaDD74NrSY68q+jsCzms8CF2O1W1+bOkJ4WfHx+e9dmw0jmdOTuKCXq66CtEZUpHSTcyljaL7HBtN1G2120x0hhTsXh/FdVu7MRpL42FntsXgZO7kxis2duIrt+1GNKDgwlXhurWE5xG2FsArUZlyuPni1UhmDGxbzW/KSrhmSze+uW8Qmm7CJ4uglOL1iTn84QMvIaxKuOvmC5f7FBsSlh49mcgUrHxOZgycmUri3ZetrfepNTx//Z6L0R1S3TYaLGX6tBOHYG6lzpCCbasjuHZrlxuDqAdcObQATB0QAmzqto1DSJXwgas38F5KFXLN1i6kDQvPn7azRj72wEu4/m9+iqPnY/jcuy/mBYRFYJ/LRIGMpQOD0/jA154HAL5IKcBVm7qwpTeM7WsiIMS+fy9aHcFpRzlMJDIgJJv2qkoiQnVsPcKVQwvAYg7rogH45NabXVwPrtgQhSIK+Olr47i4rx2PvDKMd+xcg4/fcAHWdbbufIal4hqHeK5xoJTiD7/9ItK6hT+9+UK8ZRtvtleMoCpha08YEb+EC1aF8dBLw24aa0dAgVThrIxqw41DCxBWJUgCceMNnIUTUCRcvbkTPzg4gov722FR4LYr13PDUAZWWJk/ifDFoRkMTaXwv99zMd67u385Tq2p+PJtu6BIAh49OIK4ZmAmqWMinlnW4VzcrdQCEEJw445VuGnHquU+labmXZetxcishi88dgxtfhk7+zuW+5QaHubSzO/M+vBLw1AkATfya7IiNnQF0dfud9N9T03OYXIuvaxdDbhyaBG++P7LlvsUmp4bLupFWJXw+sQc3nbxat4HqAJ8soiwT8ophDMtih8cHMGbL+jmLUcWyECXrVQHJ+cwmcjgomWsTeLKgcNx8Mmi27TwTRf0LPPZNA9bekJ46KVzODFmV0L//MQExuNp/Pola5b5zJqP/mgAhACnJpKYSKSXNRGCGwcOx8MH3ziAPRuieMs2bhwq5e7f2AlRILjja89hLK7hyz85id6Iyie+LQJVErGmzY/jo3HENIPHHDicRuHCVRE8+OEr0R7gxYOVsr4ziK9/cA8m59L/r717j5GqPOM4/v2FXbAUu1xEslwsNVUjRSKFQElpeklMi4khkSaCRpAGL4nG8l+pMWkTJS2mGqsx0k1K09ZeTGuNtCSgMWq8YcSIXIqIEBNoiUVjEbCtUR//OO+2k51d6e6cM+fMzO+TTGb2ncOb533Y7DPnnZnnsKJvO88feps1i89lTJc/OTcSF03r4dG/Zhf4KbMfmouDmTVs9rQeNiybw8Fjp+j5VDdXLjyn7JBa1o8uv4i5M8YD/Lejchn8hrSZ5WLpxdP41/sfMmncGF8nugETPj2aB9YsZMuuoyw+r7wr5/l/0Mxys3yBzxjycEb3KJbNK7flSGW2lSR9S9J+Sa9LWld2PGZmnawSxUHSKOA+YAkwC1ghaVa5UZmZda5KFAdgAfB6RByKiPeB3wNLS47JzKxjVaU4TAMO1/x8JI2ZmVkJqlIcButTUNe4XNJ1knZI2nHsWP2FWczMLB9VKQ5HgNrWjdOBvw88KCL6ImJ+RMyfPHly04IzM+s0VSkOLwLnSfqcpNHAcmBzyTGZmXWsSnzPISI+kHQTsA0YBWyKiL0lh2Vm1rEU0bxrkuZJ0glgf4PTnAW8lUM4/XqA4xWdL+/YOil3Rczn/DXOORye/nx9NiJOvy8fES15A3ZUYY4B8/VVdb4CYuuY3Dl/1ZvPOSw+X1V5z6Fd/LnC8+UdW96qnLsi5stb1ddb9fxB9dfc1By28rbSjoiYX/Ycncq5a4zz1zjncHiGm69WPnPoq8gcncq5a4zz1zjncHiGla+WPXMwM7PitPKZg5mZFaStioOkGZKekLRP0l5J303jEyU9JulAup+Qxi+R9JKk3en+GzVzrZd0WNLJstbTTHnlTtJYSVskvZrm+XGZ62qWnH/3tkp6Jc2zMXUtbnt55rBmzs2S9jR7LW0h74+XlXkDeoEvpsdnAq+RtQC/A1iXxtcBG9LjucDU9Hg28Leaub6U5jtZ9rpaKXfAWODr6fFo4GlgSdnra5X8pZ8/k+4FPAQsL3t9rZbDNHY58FtgT9lra8Vb6QEUujh4BLiE7MtyvWmsF9g/yLEC3gbGDBjviOJQRO7Scz8Fri17Pa2YP6Cb7OOLV5S9nlbLITAOeCYVFxeHEdzaaluplqSZZK8sXgCmRMRRgHR/9iD/ZBnwckT8p1kxVlVeuZM0HrgMeLzIeKsmj/xJ2gb8AzgB/LHgkCsnhxzeBtwJvFd4sG2qLYuDpHFkp+NrI+Ld/+P4LwAbgOuLjq3q8sqdpC7gd8A9EXGoiFirKK/8RcQ3yV4ljwHq9tLbWaM5lHQx8PmIeLjQQNtc2xUHSd1kv1i/iYg/peE3JfWm53vJXpH1Hz8deBhYGREHmx1vleScuz7gQETcXXzk1ZD3715E/JusO3HHXBUxpxwuAuZJeoNsa+l8SU82ZwXto62KgyQBPwf2RcRdNU9tBlalx6vI9jL7tz22AN+PiGebGWvV5Jk7SbeTNQlbW3TcVZFX/iSNq/lD2AVcCrxa/ArKl1cOI+L+iJgaETOBxcBrEfG14lfQZsp+0yPPG9kvQgC7gJ3pdikwiWzf+0C6n5iOvxU4VXPsTuDs9NwdZBch+ijd/7Ds9bVC7sgu1BTAvprxNWWvr4XyN4Xs+ia7gL3AvUBX2etrpRwOmHMmfkN6RDd/Q9rMzOq01baSmZnlw8XBzMzquDiYmVkdFwczM6vj4mBmZnVcHMwKIOkGSSuHcfxMdw+1KukqOwCzdiOpKyI2lh2HWSNcHMwGkRq/bSVr/DaXrH30SuBC4C6yrp9vAddExNHUnuE54MvAZklnknX0/Unq9bORrJ35QeA7EfGOpHnAJrLmcM80b3Vmp+dtJbOhXQD0RcQc4F3gRrJvLH87Ivr/sK+vOX58RHw1Iu4cMM+vgO+leXYDP0jjvwBujohFRS7CbCR85mA2tMPxv549DwC3kF1U5rGsDRCjgKM1xz84cAJJPWRF46k09EvgD4OM/xpYkv8SzEbGxcFsaAN7y5wA9n7CK/1Tw5hbg8xvVhneVjIb2jmS+gvBCmA7MLl/TFJ3upbAkCLiOPCOpK+koauBpyLin8BxSYvT+FX5h282cj5zMBvaPmCVpJ+RdQS9F9gG3JO2hbqAu8m6p36SVcBGSWOBQ8DqNL4a2CTpvTSvWWW4K6vZINKnlf4SEbNLDsWsFN5WMjOzOj5zMDOzOj5zMDOzOi4OZmZWx8XBzMzquDiYmVkdFwczM6vj4mBmZnU+BjTRa2B41TZsAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sorted_data['inc'][-200:].plot()" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "first_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": 13, "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": 14, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 14, "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": "markdown", "metadata": {}, "source": [ "Ce sont les incidences annuelles" ] }, { "cell_type": "code", "execution_count": 15, "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": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "yearly_incidence.sort_values()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Les années qui ont eu les plus grandes incidences (plus grosses épidémies) sont à la fin et celles qui ont les plus faibles incidences (plus petites épidémies) sont au début." ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEICAYAAABPgw/pAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAEuVJREFUeJzt3X2wXHV9x/H3l0Q0cCGigUsbtKFIqZSolGvxYUrvBXVQqFa0VUQLHW3s+ERtOg51QLT1Ia3iaB1oJyqCol4V7MiD9WHUq+IDkoBjwIhQQSQoikr0UgQC3/7xO5feXBOSu2dzd89v3q+ZTM6e3T37/d6z+9mzv3PObmQmkqTu223QBUiS+sNAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFVi8UI+2LJly3LFihVbzbvzzjvZc889F7KMvrOH4WAPw8Ee+m/9+vW3Z+a+O7rdggb6ihUrWLdu3VbzpqamGB8fX8gy+s4ehoM9DAd76L+I+OHO3M4hF0mqhIEuSZUw0CWpEga6JFXCQJekSuww0CPi3Ij4aURcM2veIyLi8xFxffP/Pru2TEnSjuzMFvp5wLFz5p0GfCEzDwa+0FyWJA3QDgM9M78C/GLO7OcA5zfT5wN/0ee6JEnzFDvzm6IRsQK4NDMPay7fkZkPn3X9LzNzm8MuEbEKWAUwOjp6xOTk5FbXT09PMzIy0mv9Q8EehkNXe9iwafMD06NL4La7FuZxVy5fukuW29X1MNuw9TAxMbE+M8d2dLtdfqZoZq4F1gKMjY3l3LOvhu2MrF7Yw3Doag+nnHbZA9OrV27hrA0LcwL3TSeN75LldnU9zNbVHno9yuW2iPgdgOb/n/avJElSL3oN9IuBk5vpk4FP9accSVKvduawxY8C3wAOiYhbIuKlwBrg6RFxPfD05rIkaYB2OFiXmSdu56pj+lyLJKkFzxSVpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVaJVoEfEayPi2oi4JiI+GhEP61dhkqT56TnQI2I58BpgLDMPAxYBL+xXYZKk+Wk75LIYWBIRi4E9gFvblyRJ6kXPgZ6Zm4B3ADcDPwY2Z+bn+lWYJGl+IjN7u2PEPsBFwAuAO4BPABdm5gVzbrcKWAUwOjp6xOTk5FbLmZ6eZmRkpKcahoU9DIeu9rBh0+YHpkeXwG13Lczjrly+dJcst6vrYbZh62FiYmJ9Zo7t6HaLWzzG04AbM/NnABHxSeApwFaBnplrgbUAY2NjOT4+vtVCpqammDuva+xhOHS1h1NOu+yB6dUrt3DWhjYvy51300nju2S5XV0Ps3W1hzZj6DcDT4qIPSIigGOAjf0pS5I0X23G0K8ALgSuAjY0y1rbp7okSfPU6rNdZp4JnNmnWiRJLXimqCRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVaJVoEfEwyPiwoj4XkRsjIgn96swSdL8LG55/3cDn8nM50fE7sAefahJktSDngM9IvYGjgJOAcjMe4B7+lOWJGm+IjN7u2PEE4C1wHeBxwPrgVMz8845t1sFrAIYHR09YnJycqvlTE9PMzIy0lMNw8IehkNXe9iwafMD06NL4La7BlhMH+xMDyuXL12YYno0bM+liYmJ9Zk5tqPbtQn0MeCbwFMz84qIeDfwq8w8Y3v3GRsby3Xr1m01b2pqivHx8Z5qGBb2MBy62sOK0y57YHr1yi2ctaHtSOhg7UwPN605boGq6c2wPZciYqcCvc1O0VuAWzLziubyhcAft1ieJKmFngM9M38C/CgiDmlmHUMZfpEkDUDbz3avBj7cHOHyA+Bv2pckSepFq0DPzG8DOxzXkSTtep4pKkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SapE60CPiEURcXVEXNqPgiRJvenHFvqpwMY+LEeS1EKrQI+IA4DjgPf1pxxJUq8iM3u/c8SFwNuAvYB/zMzjt3GbVcAqgNHR0SMmJye3un56epqRkZGeaxgG9jAc2vSwYdPmPlfTm9ElcNtdg66iHXvYtpXLl/Z834mJifWZObaj2y3u9QEi4njgp5m5PiLGt3e7zFwLrAUYGxvL8fGtbzo1NcXceV1jD8OhTQ+nnHZZf4vp0eqVWzhrQ88vy6FgD9t200njfV3etrQZcnkq8OyIuAmYBI6OiAv6UpUkad56DvTM/KfMPCAzVwAvBL6YmS/uW2WSpHnxOHRJqkRfBokycwqY6seyJEm9cQtdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEga6JFXCQJekShjoklQJA12SKmGgS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIq0XOgR8SjIuJLEbExIq6NiFP7WZgkaX4Wt7jvFmB1Zl4VEXsB6yPi85n53T7VJkmah5630DPzx5l5VTP9a2AjsLxfhUmS5icys/1CIlYAXwEOy8xfzbluFbAKYHR09IjJycmt7js9Pc3IyEjrGgZpV/WwYdPmvi9ze0aXwG13lemVy5cu2OP2U5v1sJB/6wczez10lT1sW5vX1cTExPrMHNvR7VoHekSMAF8G3pKZn3yw246NjeW6deu2mjc1NcX4+HirGgZtV/Ww4rTL+r7M7Vm9cgtnbSgjcDetOW7BHref2qyHhfxbP5jZ66Gr7GHb2ryuImKnAr3VUS4R8RDgIuDDOwpzSdKu1eYolwDeD2zMzHf2ryRJUi/abKE/FXgJcHREfLv596w+1SVJmqeeB4ky83Ig+liLJKkFzxSVpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIqYaBLUiUMdEmqhIEuSZUw0CWpEp354b9B/uZjV39js4varOfVK7dwypD8Nqg0CG6hS1IlDHRJqoSBLkmVMNAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJQx0SaqEgS5JlTDQJakSBrokVcJAl6RKGOiSVAkDXZIq0SrQI+LYiLguIm6IiNP6VZQkaf56DvSIWAScDTwTOBQ4MSIO7VdhkqT5abOF/ifADZn5g8y8B5gEntOfsiRJ8xWZ2dsdI54PHJuZL2suvwQ4MjNfNed2q4BVzcVDgOvmLGoZcHtPRQwPexgO9jAc7KH/fi8z993Rjdr8SHRsY95vvTtk5lpg7XYXErEuM8da1DFw9jAc7GE42MPgtBlyuQV41KzLBwC3titHktSrNoF+JXBwRBwYEbsDLwQu7k9ZkqT56nnIJTO3RMSrgM8Ci4BzM/PaHha13eGYDrGH4WAPw8EeBqTnnaKSpOHimaKSVAkDXZIqYaBLUiU6GegRcXREHDjoOtroeg9dr39GDX3Yw+ANS/2d2inafFfMJHAHcD9wZmZ+ebBVzU/Xe+h6/TNq6MMeBm/Y6h/qLfSIOCAi9p416wXARZl5FOWPeGJEPHkw1e2crvcwn/ojYltnDw+Frq8HsIdhMOz1D2WgR8RjI+LTwOXAP0fEzJd+/QbYo5n+OOW7Fp40jEHS9R56qT+H8ONe19cD2MMw6Er9QxPoEbHnrItPAG7JzBXAF4F3NPN/AdwdEXtl5i+A7wOjwIoFLHW7ut7Dg9T/JX67/pGm/usZkvpndH09gD0wBD10sf6BBnpE7BMR50XElcCaiNi3eWd7HPC1ZqvvYuCOiDiO8sfaC1jZLOJ6yrei3TOI+qH7PTxI/Stn1f+pbdT/uGYR3x9k/TO6vh7AHppFDOvroRP1D3oL/ShgC/Asyrc3vh7Ym1LX/rM+wp8PvAj4FvBryo9qkJnfAI4GfrWwZW+l6z1sr/5FdKP+GV1fD2APw9BDt+vPzF3+jxIOLwe+TPlu9GXN/I8Dr2mmDwTWNNc/kTJWtai5bgT4WbOc5cBG4FXAB4BzgD3soe76a+rDHgbfQ9fr396/hdpCPx54NvAm4MnAvzXzPw88pZn+EfBV4JmZeSXl3XECIDOngSuAJ2bmJuAllHGqnwCnZ+b/2kP19dfUhz0Mvoeu179NbX7gYiszRzlExBMpH0W+ClyWmXcDfwD8IDO/GBE3Am+PiGcA64HnRsSyzLw9Iq4H7oyIRwPvAV4cEftRvmv955SPN2TmOmBdv2qvpYeu119bH00v9uDrYcH0ZQt91h/uKOBcyqE8TwPe1tzkfuD7EbEkM2+kfDx5HGXs6VbKsZwA91E+wuwGXET5CsuTgCOAtZl5fz/q3U4Pi5oe/ozysalTPTR1ZUSM09F10PSxd9f7iIhHNv8/BTivoz3sFxGPjIgxynhxp3qIiN27/HruWa9jNZRjL/8O+Ajwt8BDgL8HXtlcvw/wHeBwyh9nDbCiue54yh9mWTO9AVhK2ZnwaWD3WY+zW7/HmWYte0/gZZQVtZqy86NrPewFXEb5PnqA13ap/lnPpZOBL1BO0uhUH/z/Gdd/SRlnnaL8fm7Xnkt7AqdQhh02A8d1qQdKBr0CuAT4D+AxwKldqb8f/3raQo+I/YFLgXHgQ5SdBidQxqK2AGTmL4FPAa9pniD7AY9tFvEV4M+AezLzUuD9wIXA2ZStgXtnHit33Tv4npQAORp4L/AM4HmUnR/3d6GHxhLgocBBEbEMOIiyRdGJ+iPiIcC1wPOBt2fm85qrDp/1+EPdR2ZmRCwF/gp4V2aOZ+Z1lK24TvQQEY+hDEkcA5wObAJupuwY7Mrz6ZWU1/O7KL9vfEJz+b6O1N9ej++ES4AjZ10+hbJT4WTgW7Pm/y5wazP9Ssqpsfs0978EePSs2y5b6Hcz4OGzpl9HWdEndayHk4G3A2cALwX+HLiyK/U3j/tJ4KQ5814AXNGVPihbhv/STM9ssZ/QlR4oof3QWZfPpbzJPqdDPVwC/HUz/VLg1c3zqDOv57b/eh1D/w3wreaAeygfY47MzPMpW4r7A2TmrcC1EXFkZp4N3AB8jPID05dn5s0zC8zM23uspWeZeUdE7B0R51GGXJZRVupBETE6zD3M+tvvBvwPZSv36My8BDhw2Ouf41zgzIg4KyKmIuINwDcpv1m7X1PbsPdxO/CnEXESsD4izqdsGf5h88lpqHvIzPuy7Cyc2QcQlC+cuoSyHrrwmv5v4OSI+ATwRuDxwDVN/ft2oP72Wr4jzmyJnA+c2kx/CPjXZvoRwPto3vUoY1yHAQ8b9DvZnD5eQflouZYypv514A2UJ/VQ9wB8grJ1tZRyTO3plCfxGR1bB5+lnMTxKMoL7FTgGx1aDwc39f47ZYvvRcA7Kccqv47yxjvUPczp52rghGb6gq68ppvX77mUjbM3AmcC1zXPrU6tg17+tTrKJTMzIg4A9qe8O9L8EYmISylj1IuyedfLzHsz85rM/E2bx+23zDwny57ucyjjnv9J+Qh2MUPcQ0SMULYM3wt8hrKX/kjgRGCfiLiEIa5/judm5lsz80fAWymHlb2HDqyHxs2U070XZxmrvbSZdyllp28XeiAiZjLhcspORYC3lKs68Xw6FJjKsnX9AcrGzkfpzvOolX4ch344ZY/wjRHxMspHl9dTxq6+l5lX9eExFsrPKTsVT8/MD0bEi4FrM/PqAde1PVsoYXEvZcv8PuD1mbkB+IcO1P+A3PpEjDsogX5GZn6kC31k5t0RsQZ4czNrN8pBAmsy8+ou9ABlh19EzHx74A3NvI3A65rhpO8Oaw8RsYiSP0cBHwTupBzk8IrMvLEr66CN1j9wERFfA34fuIly/OabMvM77UtbGM3RCcdQPiIfShl2OTsz733QOw6h5uSHE4DJzPzJoOuZj4h4KHAs5Yy7P6IcdnZOZm4ZaGHzFBFvobweDgc+R/nBg18Otqr5i4jrgDdk5sdmzjMZdE07IyIOoryG76Gsh/8C3pzlzM7qtQr05pCzM4EbgQuy2anSJRGxmPKdDndTeujcR69my+T+rrzoticiXk45ZPRDXVwPMyLiEOCHXexh1kmCT6DsaN/StedVs2FzMPD1zLxr0PUspE79BJ0kafsG/fW5kqQ+MdAlqRIGuiRVwkCXpEoY6JJUCQNdkiphoEtSJf4Pm1ZPCOJC1cIAAAAASUVORK5CYII=\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", "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 }