{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ " %matplotlib inline\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import isoweek" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "data_url = \"http://www.sentiweb.fr/datasets/incidence-PAY-7.csv\"" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
02022227212861628426288322440FRFrance
12022217198891589623882302436FRFrance
22022207235851900428166362943FRFrance
32022197185931418123005282135FRFrance
42022187178511396321739272133FRFrance
52022177203141600124627312438FRFrance
62022167196601486024460302337FRFrance
72022157177991371521883272133FRFrance
82022147170051316220848262032FRFrance
92022137154481165919237231729FRFrance
102022127147021079418610221628FRFrance
11202211711729834715111181323FRFrance
122022107133141003616592201525FRFrance
13202209710485760013370161220FRFrance
14202208712088874115435181323FRFrance
152022077140031078917217211626FRFrance
1620220679798704812548151119FRFrance
17202205710851779713905161121FRFrance
1820220479547672112373141018FRFrance
192022037139721068017264211626FRFrance
202022027849560261096413917FRFrance
212022017137931059716989211626FRFrance
22202152713239961116867201525FRFrance
23202151713326962917023201426FRFrance
242021507141281031217944211527FRFrance
252021497136741036916979211626FRFrance
26202148711549850314595171222FRFrance
27202147711419837614462171222FRFrance
282021467821657241070812816FRFrance
2920214578965646811462141018FRFrance
.................................
16141991267176081130423912312042FRFrance
16151991257161691070021638281838FRFrance
16161991247161711007122271281739FRFrance
1617199123711947767116223211329FRFrance
1618199122715452995320951271737FRFrance
1619199121714903897520831261636FRFrance
16201991207190531274225364342345FRFrance
16211991197167391124622232291939FRFrance
16221991187213851388228888382551FRFrance
1623199117713462887718047241632FRFrance
16241991167148571006819646261834FRFrance
1625199115713975978118169251832FRFrance
1626199114712265768416846221430FRFrance
162719911379567604113093171123FRFrance
1628199112710864733114397191325FRFrance
16291991117155741118419964271935FRFrance
16301991107166431137221914292038FRFrance
1631199109713741878018702241533FRFrance
1632199108713289881317765231531FRFrance
1633199107712337807716597221529FRFrance
1634199106710877701314741191226FRFrance
1635199105710442654414340181125FRFrance
16361991047791345631126314820FRFrance
16371991037153871048420290271836FRFrance
16381991027162771104621508292038FRFrance
16391991017155651027120859271836FRFrance
16401990527193751329525455342345FRFrance
16411990517190801380724353342543FRFrance
1642199050711079666015498201228FRFrance
16431990497114302610205FRFrance
\n", "

1644 rows × 10 columns

\n", "
" ], "text/plain": [ " week indicator inc inc_low inc_up inc100 inc100_low \\\n", "0 202222 7 21286 16284 26288 32 24 \n", "1 202221 7 19889 15896 23882 30 24 \n", "2 202220 7 23585 19004 28166 36 29 \n", "3 202219 7 18593 14181 23005 28 21 \n", "4 202218 7 17851 13963 21739 27 21 \n", "5 202217 7 20314 16001 24627 31 24 \n", "6 202216 7 19660 14860 24460 30 23 \n", "7 202215 7 17799 13715 21883 27 21 \n", "8 202214 7 17005 13162 20848 26 20 \n", "9 202213 7 15448 11659 19237 23 17 \n", "10 202212 7 14702 10794 18610 22 16 \n", "11 202211 7 11729 8347 15111 18 13 \n", "12 202210 7 13314 10036 16592 20 15 \n", "13 202209 7 10485 7600 13370 16 12 \n", "14 202208 7 12088 8741 15435 18 13 \n", "15 202207 7 14003 10789 17217 21 16 \n", "16 202206 7 9798 7048 12548 15 11 \n", "17 202205 7 10851 7797 13905 16 11 \n", "18 202204 7 9547 6721 12373 14 10 \n", "19 202203 7 13972 10680 17264 21 16 \n", "20 202202 7 8495 6026 10964 13 9 \n", "21 202201 7 13793 10597 16989 21 16 \n", "22 202152 7 13239 9611 16867 20 15 \n", "23 202151 7 13326 9629 17023 20 14 \n", "24 202150 7 14128 10312 17944 21 15 \n", "25 202149 7 13674 10369 16979 21 16 \n", "26 202148 7 11549 8503 14595 17 12 \n", "27 202147 7 11419 8376 14462 17 12 \n", "28 202146 7 8216 5724 10708 12 8 \n", "29 202145 7 8965 6468 11462 14 10 \n", "... ... ... ... ... ... ... ... \n", "1614 199126 7 17608 11304 23912 31 20 \n", "1615 199125 7 16169 10700 21638 28 18 \n", "1616 199124 7 16171 10071 22271 28 17 \n", "1617 199123 7 11947 7671 16223 21 13 \n", "1618 199122 7 15452 9953 20951 27 17 \n", "1619 199121 7 14903 8975 20831 26 16 \n", "1620 199120 7 19053 12742 25364 34 23 \n", "1621 199119 7 16739 11246 22232 29 19 \n", "1622 199118 7 21385 13882 28888 38 25 \n", "1623 199117 7 13462 8877 18047 24 16 \n", "1624 199116 7 14857 10068 19646 26 18 \n", "1625 199115 7 13975 9781 18169 25 18 \n", "1626 199114 7 12265 7684 16846 22 14 \n", "1627 199113 7 9567 6041 13093 17 11 \n", "1628 199112 7 10864 7331 14397 19 13 \n", "1629 199111 7 15574 11184 19964 27 19 \n", "1630 199110 7 16643 11372 21914 29 20 \n", "1631 199109 7 13741 8780 18702 24 15 \n", "1632 199108 7 13289 8813 17765 23 15 \n", "1633 199107 7 12337 8077 16597 22 15 \n", "1634 199106 7 10877 7013 14741 19 12 \n", "1635 199105 7 10442 6544 14340 18 11 \n", "1636 199104 7 7913 4563 11263 14 8 \n", "1637 199103 7 15387 10484 20290 27 18 \n", "1638 199102 7 16277 11046 21508 29 20 \n", "1639 199101 7 15565 10271 20859 27 18 \n", "1640 199052 7 19375 13295 25455 34 23 \n", "1641 199051 7 19080 13807 24353 34 25 \n", "1642 199050 7 11079 6660 15498 20 12 \n", "1643 199049 7 1143 0 2610 2 0 \n", "\n", " inc100_up geo_insee geo_name \n", "0 40 FR France \n", "1 36 FR France \n", "2 43 FR France \n", "3 35 FR France \n", "4 33 FR France \n", "5 38 FR France \n", "6 37 FR France \n", "7 33 FR France \n", "8 32 FR France \n", "9 29 FR France \n", "10 28 FR France \n", "11 23 FR France \n", "12 25 FR France \n", "13 20 FR France \n", "14 23 FR France \n", "15 26 FR France \n", "16 19 FR France \n", "17 21 FR France \n", "18 18 FR France \n", "19 26 FR France \n", "20 17 FR France \n", "21 26 FR France \n", "22 25 FR France \n", "23 26 FR France \n", "24 27 FR France \n", "25 26 FR France \n", "26 22 FR France \n", "27 22 FR France \n", "28 16 FR France \n", "29 18 FR France \n", "... ... ... ... \n", "1614 42 FR France \n", "1615 38 FR France \n", "1616 39 FR France \n", "1617 29 FR France \n", "1618 37 FR France \n", "1619 36 FR France \n", "1620 45 FR France \n", "1621 39 FR France \n", "1622 51 FR France \n", "1623 32 FR France \n", "1624 34 FR France \n", "1625 32 FR France \n", "1626 30 FR France \n", "1627 23 FR France \n", "1628 25 FR France \n", "1629 35 FR France \n", "1630 38 FR France \n", "1631 33 FR France \n", "1632 31 FR France \n", "1633 29 FR France \n", "1634 26 FR France \n", "1635 25 FR France \n", "1636 20 FR France \n", "1637 36 FR France \n", "1638 38 FR France \n", "1639 36 FR France \n", "1640 45 FR France \n", "1641 43 FR France \n", "1642 28 FR France \n", "1643 5 FR France \n", "\n", "[1644 rows x 10 columns]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "raw_data = pd.read_csv(data_url, encoding = 'iso-8859-1', skiprows=1)\n", "raw_data" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
\n", "
" ], "text/plain": [ "Empty DataFrame\n", "Columns: [week, indicator, inc, inc_low, inc_up, inc100, inc100_low, inc100_up, geo_insee, geo_name]\n", "Index: []" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "raw_data[raw_data.isnull().any(axis=1)]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**This csv file of chickenpox_data does not have any missing data, so No year entry has been deleted\n", "The code line that was used to remove the year entry has been commented out**" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# data = raw_data.dropna().copy()\n", "data = raw_data" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "def convert_week(year_and_week_int):\n", " year_and_week_str = str(year_and_week_int)\n", " year = int(year_and_week_str[:4])\n", " week = int(year_and_week_str[4:])\n", " w = isoweek.Week(year, week)\n", " return pd.Period(w.day(0), 'W')\n", "\n", "data['period'] = [convert_week(yw) for yw in data['week']]" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ " sorted_data = data.set_index('period').sort_index()" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ " periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", " delta = p2.to_timestamp() - p1.end_time\n", " if delta > pd.Timedelta('1s'):\n", " print(p1, p2)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXm8HUWZN/59zrk3N/u+koQkkLCELUAICKKySYBxcJ2JG4j4iyI6OvqOL6jj6CgO+qoojqgoAgIKjMoAssm+CATCmgWyEUhCQnITsid3O6d+f3RXd3X1U0uf7d7c29/PJznnVldX1emuep561iIhBHLkyJEjRw4Vhe4eQI4cOXLk6HnImUOOHDly5EghZw45cuTIkSOFnDnkyJEjR44UcuaQI0eOHDlSyJlDjhw5cuRIIWcOOXLkyJEjhZw55MiRI0eOFHLmkCNHjhw5Umjq7gFUitGjR4upU6d29zBy5MiRY5/Cc889t1kIMcZVb59lDlOnTsXChQu7exg5cuTIsU+BiN7wqZerlXLkyJEjRwo5c8iRI0eOHCnkzCFHjhw5cqSQM4ccOXLkyJFCzhxy5MiRI0cKOXPIkSNHjhwp5MwhR44cOXKkkDOHHL0CW3a1455FG7p7GDly9BrkzCFHr8CF1y/ERTc9j627O7p7KDly9ArkzCFHr8C6rXsAAF1l0c0jyZGjdyBnDjly5MiRI4WcOeTIkSNHjhSczIGI+hPRM0T0EhEtIaLvhOXfJqI3iejF8N/Zyj2XEtFKIlpGRGcq5ccS0aLw2pVERGF5CxHdEpYvIKKptf+pOfoCBHK1Uo4ctYCP5NAO4FQhxFEAZgGYS0QnhNeuEELMCv/dDQBENBPAPACHAZgL4CoiKob1fwlgPoAZ4b+5YfmFALYKIaYDuALAD6r/aTn6Fqi7B5AjR6+CkzmIALvCP5vDf7bt2bkAbhZCtAshVgNYCWAOEU0AMFQI8ZQQQgD4PYD3K/dcH37/E4DTpFSRI0eOHDkaDy+bAxEViehFAJsA3C+EWBBe+gIRvUxEvyOiEWHZRABrldvXhWUTw+96eeIeIUQXgO0ARjHjmE9EC4loYWtrq9cPzJEjR44c2eHFHIQQJSHELACTEEgBhyNQER2IQNW0AcCPw+rcjl9Yym336OO4WggxWwgxe8wY50FGOXI4seC1LdiwfW93DyNHjh6HTN5KQohtAB4BMFcIsTFkGmUAvwEwJ6y2DsBk5bZJANaH5ZOY8sQ9RNQEYBiAtzP9khw5KsA/X/00Tvvxo909jBw5ehx8vJXGENHw8PsAAKcDeDW0IUh8AMDi8PsdAOaFHkjTEBienxFCbACwk4hOCO0J5wG4Xbnn/PD7hwE8FNolcuTIhgpmzZ6OUu3HkSPHPg6fM6QnALg+9DgqALhVCPFXIrqBiGYhWI6vA/gsAAghlhDRrQCWAugCcLEQQq6+iwBcB2AAgHvCfwBwDYAbiGglAolhXg1+W44+hNx9IUeO2sLJHIQQLwM4min/pOWeywBcxpQvBHA4U94G4COuseTIkSNHjsYgj5DO0asgANy/dCP+6ddPIddM5shROXzUSjly7DMQArjoxufQVRboKgs0F3N9U44clSCXHHL0KuTpM3LkqA1y5pCjVyHXJOXIURvkzCFHr0JP5A17O0r4wh+ex1vb27p7KDlyeCNnDjl6FYToeYqlexZvwF9f3oAf3Ptqdw8lRw5v5MwhR69CrlbKkaM2yJlDjhw5cuRIIWcOOXoFpMNqORcdcuSoCXLmkKNXQeUNOZ/IkaNy5MwhR69Czg9y5KgNcuawD+DvKzejVM7Jng96YsqMHjikHDmcyJlDD8fjK1rx8d8uwFUPr+zuoewTEInvOVXOkaNS5Myhh2PjjnYAwOotu7t5JD0bMmV3T9yl5+nEc+yLyJlDjl4G0eNUSz1sODlyeCFnDjl6FXJvpb6HtW/vwW0vrOtxm4J9HXnK7hy9Cj2RPORqpfriO3cuwQOvbMKcaaMwcfiA7h5Or0EuOeToVeiJm8eeOKbehGdf3woA6CqVu3kkvQtO5kBE/YnoGSJ6iYiWENF3wvKRRHQ/Ea0IP0co91xKRCuJaBkRnamUH0tEi8JrVxIFeyoiaiGiW8LyBUQ0tfY/NUdfQE/2UGqkAHHBtc/gsG/d28Aeux85E64tfCSHdgCnCiGOAjALwFwiOgHAJQAeFELMAPBg+DeIaCaAeQAOAzAXwFVEVAzb+iWA+QBmhP/mhuUXAtgqhJgO4AoAP6jBb6sp7lvyFlZvzj2Gejp6MoFo5NAeXtaK3R2lBvbYfZBquzx1Sm3hZA4iwK7wz+bwnwBwLoDrw/LrAbw//H4ugJuFEO1CiNUAVgKYQ0QTAAwVQjwlAsvR77V7ZFt/AnCalCp6Cj57w3M45UePdPcwcjggRM+0O+SoH+K8Wt06jF4HL5sDERWJ6EUAmwDcL4RYAGCcEGIDAISfY8PqEwGsVW5fF5ZNDL/r5Yl7hBBdALYDGFXJD8rRt6GqlXraRrJH7XZ6EeJ9ZA974fs4vJiDEKIkhJgFYBICKeBwS3VuDQhLue2eZMNE84loIREtbG1tdQ07Rx+EED2XCOekqz7IJYf6IJO3khBiG4BHENgKNoaqIoSfm8Jq6wBMVm6bBGB9WD6JKU/cQ0RNAIYBeJvp/2ohxGwhxOwxY8ZkGXqOHoaNO9rQ0VUf75K+QCO27+nE35a81d3D6FHoaZLivg4fb6UxRDQ8/D4AwOkAXgVwB4Dzw2rnA7g9/H4HgHmhB9I0BIbnZ0LV004iOiG0J5yn3SPb+jCAh0QPjWj55v8uwp6Oru4exj6NzlIZx3//QXz1f17yqr+rvQsfuOrvWLlpp7NuIgiuh7GJWko0X/jj85h/w3P5udTIDdL1go/kMAHAw0T0MoBnEdgc/grgcgBnENEKAGeEf0MIsQTArQCWArgXwMVCCOk2cRGA3yIwUq8CcE9Yfg2AUUS0EsBXEHo+9UTc+PQa3PDUG1W301kq4+V122owon0PMsPsfZ4738eXt+KFNdvwo/uWG+tQSHp7O4FYuWknNmzfi3Vb9wJAvlEBgD7y7hsNZ4S0EOJlAEcz5VsAnGa45zIAlzHlCwGk7BVCiDYAH/EYb49AVw2Umz+891X85vHVuO/L78LB44fUYFT7IGq4lkshYVixaZej5r6N03/yGABg+tjBAGozF/d19OSki/sy8gjpbsKS9TsAAJt3tXfzSBqPaDHXkDu07gye4/9RVFUuYtFDNZdeaCoED7GrtO/+hloh8lXKH0VNkTOHClALotKX9aR98CfXHE3FYALlh0D17bVUT+TMoQLUYg4Wwhndl+dz1t/uI2m0NBW8292Xn32xECzdrrKfx9dP/rYMT7+2pZ5D6nbsw6+zRyJnDhWglpPQe7fTC2d+PX5S/+Zi9N3V/r78SCO1EiM5vNaatrtc+dBKzLv66bqPqzsgN1q55FBb5MyhAtRiDsqozr44neXzq4fO/9AJQ/pEimyTzWHlpp049cePdseQug2xzaEvrqb6IWcOFaBQA+JT8Iz4f3R5GAneCwme71L2IfbvnTkOAHD6oeMyqJUaQ0zq0YvJ5rB+W9+Ne8h5Q22RM4cK0FSs/rHFIf/2GX3nS+ut1/dFSNtBLRcz587Ym3eS0ubQ6Wlz6M2gSK3UzQPpZciZQwVoLla/je/LBul6/GYZBJfFPbZRj74eQl8xYobJX8FJWb2ZSaroK7+zUciZQwUo1ECp3Zfd7+r5i3vi46zHkKLdsiY4cL+/t7u7xmupe8fR25Azh25CXzZI1xOq947TW8nz4X/xjy/gtB8/kih7eNkmbNvTkW1wNUSWJNWlnsgxawgySFE5qkPOHCpALaZgX/awqPQ3+9z2/+5b5t+e55u886X1WNUanwK4s60TF1z7LD593bNe92eRM9dt3YNfPbqKvfbeK2IvJMrgvtlXzBJ9byXVF87cSjnSqAVB79M2hzq02Uj1nHQffc3z2NgsI7vg2mexYtMu/ONR+2G/4QMS15ZvjOMXnnsjyGjvMxd7veSQJ96rC3LJoZuQ60mT6CyV8aP7lmFXe2VZRjnCwNGK9dv2YmmY12rFxuqS9GWmRR4ihO/v37qn03sMuc0hRyXImUMFqE0QXNiW776yF0187vn97wtv4r8fXokfZVALqfAlgCde/hDOvvJxAMA//PyJivqq2B9BALc+uxa3v/imsYpkcr59uH71no4ulHs51ezLKtp6ImcOFaAW2USz+mbvqyLzvYs3pE8sY35KZ6iqae8qpS96oJGpq6t5FV/788v40s0vOtsmT0uFa17M/NZ9WL3FT/21r4L6sIq2nsiZgwf0HUltE+/5Nbavbv4+d2NwYpkKG3PlH4ebULJEsqc9Mw96L4fsKzn4zIsPXvWkX2P7KOSj2ra3Az97YEWvV6M1CrlB2gM63anF1CtE7nd+9fdVycEX1YaONPJcg2rUSs4qkeTg2WSGeVGLtC89Gd+5cym27enEoROG4L2Hje/u4ezzyCUHD+jLryY2h/DTl+j3JuZg+ymV/sysz6ez1Hj/zrsXb/ColY07ZPnZtQje7JEIf9aejkAl2d7VR3x36wwncyCiyUT0MBG9QkRLiOhLYfm3iehNInox/He2cs+lRLSSiJYR0ZlK+bFEtCi8diWFykIiaiGiW8LyBUQ0tfY/tXaohc0hqytrb/JV536yD9myPSrO5mB7T22dldk2Eu1nZEhtne6XuLs9GJevzSHLXKwVc3hp7bYeZfwl7bM3baS6Ez6SQxeArwohDgVwAoCLiWhmeO0KIcSs8N/dABBemwfgMABzAVxFRDLJ/i8BzAcwI/w3Nyy/EMBWIcR0AFcA+EH1P612qIfNAZH7nV9jtTxSsyej0t+ZVc/sQ6i7A3tDpuX7HLJsGmrBG+5etAHn/uLvuO0Fs8dVo0HaD8uZQ23gZA5CiA1CiOfD7zsBvAJgouWWcwHcLIRoF0KsBrASwBwimgBgqBDiKRFQ298DeL9yz/Xh9z8BOI30N96NqMdUK2RMn9GbbGzcrrPaQ+Kz2hxqITnUdYrWwRZVC8lhdRj4t2ITHyPS1lnC3J8+1tBT5yLJQW64eibf3+eQyeYQqnuOBrAgLPoCEb1MRL8johFh2UQAa5Xb1oVlE8PvenniHiFEF4DtAEZlGVs9kTJI1+IM6Yxt9SZfdV6tVB3h4iQH/dHe/Mya6HstmENPQDabQ/37XdW6C6++tRPfuXNp7TrzhJxDvT0ivFHwZg5ENBjAnwF8WQixA4GK6EAAswBsAPBjWZW5XVjKbffoY5hPRAuJaGFra6vv0KuGLuLX0pV1X4lzWLd1D/6wYI27YpWo9Fe6ns/L67bhkr8siv6uhVqpnnp335a/c+cSbNzhd8BPLbMJm9ReMl6lFmntfSHHFEsOOXOoBbyYAxE1I2AMNwkh/gIAQoiNQoiSEKIM4DcA5oTV1wGYrNw+CcD6sHwSU564h4iaAAwD8LY+DiHE1UKI2UKI2WPGjPH7hXVALaZec1Mwk329Zrp7vs+7+ml8/bZF2NNRWXoLFSxNtdASH5rmejw6M3jglY3uRk19NeBd+Paxu6OES/78cnyfpW4ttGCRhGfoqCucz6tbd6eDH+sEOSbJ/HLJoTbw8VYiANcAeEUI8ROlfIJS7QMAFoff7wAwL/RAmobA8PyMEGIDgJ1EdELY5nkAblfuOT/8/mEAD4ke5A6RVitV32ZzeJpch6fbXSMlh/9ZuBYX3/R8Ygcmd6e12H1mD4LzaJPLrWSp/7MHV1TWkaPdWiGLYd43OrxQA70S2XkDOkLmsLO9KxX8WGt0lco4/SePYtnGnYny7t5I9Rb4BMGdBOCTABYRkYz7/zqAjxLRLATz5HUAnwUAIcQSIroVwFIEnk4XCyGkgvciANcBGADgnvAfEDCfG4hoJQKJYV51Pys7Vm7aiY4ugZn7DXXWrYXnUL+QOXR6GlIbySr/7U/BTvTH/3QU+hcCRzM5zpqMg2kjPp+gsg56Gz2o9DnbyH9N1Erhp2nv1sjo5B1tXVjJGMZ70L5yn4aTOQghngA/5+623HMZgMuY8oUADmfK2wB8xDWWeuL0nzwGAHj98nNS1+ohOfRrksyh50kOtj7r5VLr4/nTU9Z8I4hPpT3Y7quFQdrlVZbFa6xUFlj21k6vDRk7FsPY8vQZtUEeIV0BamJzKPZ85uBaZM+9sZXdublQl1/CpVaq0zNriFqpDmOvheTgcqToyBB5ftXDK3H2lY9j0brtVY8LiJlFzhxqg5w5eKAeu+XI5tCDDdJcnyrN+tAvn8TpP3k0XckBju5FZKtSm4Pjej1CEur5SnqqWknCtCayMLXF6wOm8Oa2vRWOgUd3e/b1FuTMwYDtezsx9ZK7cO3fV6cXag0mn3T18zVIHzRucNV9ZkVWI281cBk6XWiknrkn0576q5Vqlx47lkIaExXfG3DrwrWY/vW7G5IbLGcOBkjvnD8sWJNOvFfDflwv+fhpIwEAIwe11LBXP/CSg/+vv2nBG2x5PSQxrsV6kY5GpDKpBwOqRUS3uwX/PqplDvpczHpGyr6I7/11KbrKAnva6x/EmTMHD9Q6t9LDyzZhW3jMo8uAF/XVoO3q5fe8Gn1no44929m0ow3fuG0xe41VK0WGzvRFr6R8joH1mFwsnsjCgHyJfsFjtT/86iasfXuPpS/fUfmMpzpibopn6M1qpUhya8AGJT/PwYFSWWD9tmQEajUvZvveTlxw7bPR3y61kpzojdoN/erRVdF3k5SweVc7BjQX2WsStkAkPgYuW66pdJsNJAj7KO3xsTlccN2z6N9cwKvfPYu97sp8moV5SDVXpRHN+m1yvvZi3lB1DrIsyCUHA+Qcf23zbpz508cS16p5MboayWWQlouwO7KymgzSs7/3QOqZcPUaCfYguDqNIWuzldhD6jF2X4O0LbWIy+bgyxs6uspYHxqiTYxGCIGP/OpJ3Gs4B0NnKrs7eke+LB80YnnlzKECVPpidrR14uyfPZ4oc9kc5Pzvjt0Qu/sPi9ZtzeZhoia6q9x4zN+3ZP32zOOpCer4TurRdE3SZxhyK23a0RalzvDBt25fjGdf3wrALBWXygLPvr4Vn7vxeeN17u/enFspa8LOapAzhwpQ6Xt5fPlmbNrZnihzRUhHonJlXVYFbpFVKsFc8cDyuA2rzSFbu6WywDlXPpHpnvFD+2frREPWMfraBFYoaSDqsfiLNUm8l5YcdrZ1Ys73H8R/3LHEeu+NT7+BqZfchR1tnYmU3ibJwZUjyXRfL+YNis2h/siZgwafRVkpgdzLpInuyZJDNX3qdGhnm1/CvqxdPr9mq6WtuDV1PMUqfTqjdr2P8vT7VWqa6yzPwbf9piK/3He3d2FnW6dXG7HNQb0/mNf3L91oHfe1f18NIJAyBvePzZ2mnb56LsPnb0rnaTKqo/ZVo1AGNEI6yg3SGoTwEL8rfC97mYymLs+K7rU5MJJDD1t3ldB5H68dL9T4WajvOMtz9vV5bzI8rCO/8zeUyoJNHaNj7VbpycQzXt9xD+ynMAfDPV0Kd7h7UTrDq+m+njZHawmOOdcLueSgQT5zG4Oo9L2w5xw7Gotu6SE2h0qHoT5O22/Oqk6pJOq32kjheqmVkm7N/p34BlI2Gc5YyBJM9utHXzNecylA1Suqisvokur4WaZx94XEe41w182Zgwafh371Y+YFYgOnznB11502BzZC2nNSZj3ZzaZLtRHXJpsYYBhqtbr3enkrqcRu+95Ob6Lvm4LFJDlUAtNP8vuplNh8mUZVqc2BKxVC4Cd/W+Z9MFJPRSOTC+bMQUM9GTK3Y3WpiyK1UrdkZU2X1WIU3G+ulGxVYj+o9lyDer0L1ePqQ798Cp+70e88hHbPU+2sjBTAk6s2e7UDJAlzQir0vD/BHAyvw0UATZIF93peWrcdVz60El+6+QXPEdrR3lWyBgvWD7VLX+JCn2cOVz+2CnMVn/166vYrkRy61ZW1TrsT62/J2GUlzKEWXju1xlvb2/CWtqt96NVNXveqkoONcZnUShIf+82CVNnTr23Ba63cmQl8G+n09nxFH8nSNf+yREjLtmpxPCwAfOXWl3DyDx9u+Fnk7WF/+gFH9UCfZw7fv/tVvPqW6kIov9WegHBEyUULy92oVqqXQZpXHXncx9xoYw7J6nG96iWHqm5n8fbujorv9VU/NRu8lWyYd/XTOPXH6cy7pkfge956Uq1ksIXUUK0k+3tx7Tbc+DSf8ysLHl0WnGHfiAR4Kna2B04tqitwvdDnmUMjwRIlp80h+VlP6BOOj4FLF/59ZVodUfmJblkN0tn7qIBGVoVrnljtrPOzB5c76/jAZp+pacpubm6IdHk1hlOXu6bRBdaRn+ub/8vn/NqXUEv7kQk5c9AQz6vG7NW9bQ7Wc5cFrn5sFaZeclcicZ6Ots4SVm/ebbw+7+qnE39v4XazzDCWvZUWcVPqBW28OqiButRaqZV8h/oq83xUbNrRhvuWbKx+QGicbSph2FVdWbV6PprJim0OZvGlYeguv6hanAfu7KPuPexjqKfNgff+sd8TG6TNdZZu2IHv3x0wBTVxno4v3fwCTvnRI956Us5A6asSympWqJReW/sxXGyUWml3exeeZKQqHf9+e6N2srWb20+s3BzPZ2VDlcpgbOhTlXBMb0N3/W7vSs5bE/NgJYc62Zm6y2u2R0gORDSZiB4moleIaAkRfSksH0lE9xPRivBzhHLPpUS0koiWEdGZSvmxRLQovHYlhW+MiFqI6JawfAERTa39T/VDPV+2JVWREdIjw7YjdPmDSzyxIiBUXLwFh18/+hre2t7mZCbcNE3l2vfqMfvzr+R9ZVlY2/ekI4d9NxC/efw1fOy3aSNvejzV7dF8n0Gt57YQgb1j215zdHWiT+W7zxvQifwvHk5ufExroqEEu9uYQ/339T49dAH4qhDiUAAnALiYiGYCuATAg0KIGQAeDP9GeG0egMMAzAVwFRHJ/M6/BDAfwIzw39yw/EIAW4UQ0wFcAeAHNfhtFaGe75ojKi41gE+cQz3nyStv7UgY3SolRMJyDVASitXQzTXRhtJIvyb/B/bNKnb1a9/2SwbY35H+3AVfZlXruS0QpLV47xWKt1/KWyl9HxHw6PJWZ/u6ZPDHZ9Ykr2fIrVSvfXYjz45Q12Ej7GbOLoQQG4QQz4ffdwJ4BcBEAOcCuD6sdj2A94ffzwVwsxCiXQixGsBKAHOIaAKAoUKIp0RA8X6v3SPb+hOA06hecqAD9cwJX5Hk4GGQrjZXkA0EjbBzxDvzq+LayNiEpa34Cn+tpcmfGO9uT6c8ke9iF3NNhe9veniZn8tqT4MQAg+8skn5O13Hh3iaauieQPrjNKbPaMB2Pt7MNA67lPxkxR4iOUQI1T1HA1gAYJwQYgMQMBAAY8NqEwGsVW5bF5ZNDL/r5Yl7hBBdALYDGMX0P5+IFhLRwtZW986jEkTEuA5tc2162xwsI6qlFwoHdYyVSg6+6TOyopK2svimD25Jpx/z7dL3rVTjxgpkeSe1ndXsfNZKOebw8wdXOOsASDlP6JKEyVupUiZVCRopOexR5m2PkBwkiGgwgD8D+LIQYoetKlMmLOW2e5IFQlwthJgthJg9ZswY15Arg8dOvVKwcQPOe9zjaeQE5Y196Xopf3fP9jPbHLJVBwAsWP22d91q+G69mXZW1HuWCNjViRL/++L6ZB3DwPTIbz1NiMkgzZXXK9tEI9eeKsV+/Pgpde/PizkQUTMCxnCTEOIvYfHGUFWE8FPKl+sATFZunwRgfVg+iSlP3ENETQCGAfBfwTWEJGrViKbHfPd+fOw3T6fK2XlUA5uDuhhqTY+IKNG5zV6gwppcj+/Jcq0ymMaQZUH7GNuN9zaIN/iM5qhJw+pikHaVCQ9nCdOw9PekB/uZz3PgmEOyrOro/wYe1ymxRznpbhAj0dYaPt5KBOAaAK8IIX6iXLoDwPnh9/MB3K6Uzws9kKYhMDw/E6qedhLRCWGb52n3yLY+DOAhUUeH7XP/+wlc/Rjv8lmLoLO3d3fgyVXpCEZeDLfDx5VV9VbyoUcpd0NL44TkwmIXFUMFbQTYZKR0oRK1XNR+BfdkGQfbZ08SHIjqYJBOt6indfDZZBm9jrS/dclBnWNnHT6eLZfQ522Xr4ufA8s37mTtUvWAbyR8reAjOZwE4JMATiWiF8N/ZwO4HMAZRLQCwBnh3xBCLAFwK4ClAO4FcLEQQrK8iwD8FoGRehWAe8LyawCMIqKVAL6C0POpHiiVBV5atz2KC9Ahp1B9DNKMWsnRTzypzRVVrw2TcfjhZZuiM3bVljbvase0S+/GTQvMKQXUxfbYCj9bTz04O5drqBIJb95xk92VQriM7c+9sRVTL7kLL67dxt2dcWT1A6EONgemOT1jsUqTTb37Snjp6Ov4u3qQEa9WqrHkEOKT1zyDz1y/sCZtudDoVB1O2UQI8QTMs/w0wz2XAbiMKV8I4HCmvA3AR1xjqQW27bEb/2I1TnryXPjOac5UCEvXm80xvukoVMhJbFvXCbUSc/35NVtxwbXPsuNYuSlIqna7pgeO2qPkInx53XbreLk+9HHZo729mmfrn37o2KT3jOGe0w4dh5ufXWu4moRLZfa78HSzh1/dhFmThyfv9coXVT2REh4EuLukGB8VnjFHkvfGCWhWPPY4uq8LCq7jeV1QH+dTDchzBPinZq8V+lyEtCuZl01ymDlhKHuP5OirWnfh7CsfN7ddgeTQ5cEcXAtwyy6NISrVpUg8qJ/ZvVMd932L0ydy8XTHbJC266rTF1XCNvWSu3DnSzwjmzJqEP7l1OnsNRWu7KTJzrnCeIx3vbwBAK+m8PEwbq+BqmDzrnZnHqJAcqi6qwTSMQ3Z57etjotxqvNedefmnoW+7k1tX3H/clz8h+et/XLYw5zyWGt09kC1Uu+CY7La5mOhAHz8+P0xenC/qOyJFZsx4xv34IU1W7F5Z3vmrn13R7bdtktE7tJ2HGpb0lffZuBSm9/J6FebGWKbnRD53/DFP8Y5+dV+wnB75/3NDh/xTTvbcNi37sXiN3kpifttXcxwraskAAAgAElEQVRO1Cct9d6O2qR8lhKMqcd6hA1VY09ItpOtXEJlDgm1EmeQ1t1gDY3/7MEVEcO3QX+e8hzteqJaaScr+hxzcD3eyFuJqVggQoEoMbEeeCVImPb8mm3ebqm+4xFC+EkOSsNcNT1dhtqWNHKZooYJ5JRMfIL7TAbhx1e0JnI4ZXdljW/Y3VHS+uFtMS7J4ZFlrdjdUcJ1T77ufaIdt3B96PGeGp0HIA+esc6nGluCfNxWfVT7JgbikobUPY+6QeE2S3rAomtD1SgjcxZ0lBp7dkTfYw5u7hB+pCsGzCG5Y5HJwFqaCs62sx676WPMA4Brn3xduYcRqS0LQV4xEUEit9qKPWtaK1K9t9Rn+8lrnsHHfrMg9hKz9pSG2s+ONnOOHxWcpGOCSuDf2h5kIuUZcBmX/uVl/E6xSfn00igPlLqolbzq+EgXlbWfkBwUaZBrb412aptLovHJiaWiEcd2dnblkkNd4U6RHdZjqhULBCJKXJOBOv2bixXtzPR+2jpLOPmHD+HhZZsSemwTge7oKuP+pXG6Z66a7uXAtWTb5bqICi8RJQvVaFcfSaMSdHSVvXbretIy2z3qJempZXrGf3xmLf7zr0uVdi0HEUUuyrVd8Ga1UnXM4aTLH0owPiA9dq59L8nB8Pb1ew+fmLT5qZKFat/hCLVe5rI9vsR6n8XQX62rvVogN0jXGc7dfRQEl0aBEKqVVMkheGEtTYWK7Bn6wnhz216sfXsvvnvn0uSENrT9n39dYu8U6YWRlSA5PU6Y6395/s1MfcRNVb7IukrlhARkakpf2Oqfy97aiSu19A4StgR5rM3BwnTkK9Hp2FGThplv8oDp6RHISIQHeCT+e3Pb3gTj4/razmRn9XmfxhxJ2r16xLl6n/rdR5Kt9Ua/1AB7QKNdWfsec3BdjyQHs1opobMPX1hzsZBJDJZYvnFXKk+9HGeXw5YAAM+/Yd/hAECnzhycd8QIguDsdbjrur+7C7VQK+m2lbsWuQ2LOs786WNYtzXIpkpIEnhJRDkiy6VBt6XPiOdC8r6zj5jgHOPn33Og8Vq7yYZhkRyGDqhftG3CzTajy6pebtv9q++EVd/qOZ9qzB1qFVRng2QOn3nntLr3BfRF5uDryspcKxChUEhKDrI9nWnY2tZxq+JzH2V7FCKxG6lmR71di+3I1JSHzSHr2CqJ9zC2pdynq5WuV2wxKnT7iqQTXAyMaiuJmAMzVE7kt2m44k1Istz3nOfvvj8VLgQA+NyNvBumnl1XRaV00ue1+2VlNTANrVxnDnsV91G1G06tlJYcqmMO+rtthM1B2qf+71mH1L0voE8yB9d1YawX2BwqP2nK1LeaQlptx0dy4NDRVcbl97waeVzox31mJcRuf/NMzbGQ51dzXfl6YXaVhVnn7tEeZxyWUgQA9G82LxfdXdjWDwC0daWj1QG/syYEgE+eMMVZT+Ivnz8xGAvzbG97YR1aHS7Y1oG4qlRRR59X+rrbqaSwTqR44TYf2r3Vztmt2iFQvgdoVYOO8Ic14hQ4wCNCuq/BFYxVIEpEW8r6wdpzEVH++tABzWy5uhvx1Z8DwC0L1+JXj66CEAIThvXHn59bl6yQYR4Hrqz2Oj7NjRvaotRP33Hj02tSZT5IqJVKfgZpE/QFrrcVSZXMD+ZtDubBnPXTx/H3S05NteWj/8+CY6eMwDH7jwilpWRnbZ0l/OstLxnvdUvZ7jf/6ls7MHX0IHs7nuomfXeuMgeX+kovqWan/+nrnk2VNcRbqVRGv2Khbkee6sglhwz1igVCc4HQVS4raTYC/PC+V/FshlTQKob0T/Po3R0l7GqPdydZJQcgMJZ/+86l2NGW9NnOOo1roVYaOUhhDpbq3CXf+h1VGAVLZZFa4Ju0HXVsF0n3w6qVLGv4zW172bbqmW3TpabR4dwUeDxuk5rLp5/02RDJ66qtTrU/+KiVKlHTlssCO9o62RxfjZAcOrvKmdywq0XfYw4O0vjKhh1hvTQKRGhpLqIs1LQWwefyjbtw5UMr7X0bJiQx31t3tuP0n6jHLxrutblhGq4lmvJQs72wxm70zh64VjuozyWQHLJ7KwGBQVFf4I8sSyYZtDFJLtJZtW1MHTWQvU+3Yw6wpDGxwbSblM/H5cr6wWMmpsp87XOV4othqpNKDdIJO0PJwRx0xlgBc/jpA8tx5Lf/xl4rNcgg3ZzhiNtq0feYgzInpl5yV+r6T+5fnqonQRS6rCI+TSzLFKvGBuZ7a2AYt9fmvHhMPObqx1/DpX9ZZG3Px7gnma6rfqX5eQD/3ZtKtCVNLZWF01XQZEQGgMGM9KfS6xMOGIVzjkx7IulEq6VOi187liOF/YYNwL+deXCizPQ4ZbxBtTEa7zgwOOzRaJDW2rcxh5FKShvW4UG3X3jQ8itCWiDxNyWeSEcjUlt0lMreDgu1QN9jDo7rksBwBKxIFC1eGd+QZX2YxWcPePajJiAzRT1/969LvY2QKlGvcmhRIFwt44V+8XAsrdmijVVCrX6Xxr3OknBGK9uIIXf0qGo3FALsg9KbtMVSmO6x1g0/CZQOWlO+NxUJIwb2067zHcnzi6t9jcXwRZjWhcsgLRDY1G64cA7+5dQZUTkb56D9fW2Yi8qGnz24Apt2tkV/m2yD3NjqgY4ugX45c6gfXLsdW4rsQoEiz6JKsmmad0hqHR6+Pvu+R1MuWe+XetsHvutiw/ZQz57V9dVyTU3R3VmBQVoyUy/JwXKNmw+V7PJsHlHVgJMcVILWXCykssiaXpNMj12tOrEgO6zAIL1tTwduXbgOG7a34eQZYzCgXxH3ffldGDGw2cvm8D+6k4YBqkQw3MIcapy2yjCWspc3W63Q57yVXO9QEghTbqWW5tqrldS+suxA3t7dgcVvJnf2RU83N72XauZ21l2TTftTzQYscGV1/36VgRw4ZjCWrN+BrlLZw0Br3ji4GYtgdXeyrZ9/9GgcNG4Ihg+0ECC1rQqQMsoqQ5Zu2rb6at1qxiERBZUaruvzSn09//zr9DG8B48fgqMmD8fbu9PxKrVIOjjMKjlU3bwTnaXcIF1XuIhPJFpzkgPFO/OyZpD26ttjTFnam3f1U6myo/cfHrvX1mAe+WYl9UKks7fYHJinxNXfviedqsE33736m2S8QFfZTT5s3kouPbdRrRQWDmgu4uDxQzBuaH+272cq9IST4AzWycR19usqoqy2VRLEAgXjMh72o/+t1NOPI43bNLRX4VjVPm3ModYZbzl05jaHesP+Ek+cPspYq1igiDlUMhV8zsrNknZi+cZdqbLBiiukF1kPK/3JU8zm4JuK4Bv/uziob6nu6/SxYlOaOHSWq1MrOQMkLdeypmOP6oSV1FyAB4xJxwX806/TGwEdfG4j83jUv5sZ/3nT+H1Vly4UC2RNz6Kul0PGD/GSUAtE4IQ4n3fhkv6Kll17IySH9q4exhyI6HdEtImIFitl3yaiN7UzpeW1S4loJREtI6IzlfJjiWhReO1KCmciEbUQ0S1h+QIimlrbn5iEbX41KTIuVy9YQMF3m4qhkr4lbl1YOZEGkpP0t7YjTWs4mbMbpM138G6Ifnj3QWOMDFHtkgh45P+8B3/4zPExsxduqc30zo+aPBwrGWalVhQAy61lmwkJrcJ3c9qPH00PIfwkbTxq30A4tw1jM6HaKVQgQrFAqc3Fpp1teHt3B4QIGMjL334vTjhglFegWbHAb1b0dztmSEuqTipYNAMaYZCWQXCNgk9P1wGYy5RfIYSYFf67GwCIaCaAeQAOC++5ioik+8UvAcwHMCP8J9u8EMBWIcR0AFcA+EGFv6VqFAoEgWByfeKadD73GWMHK7mPgk9bal/d+2WxwQhcy7TNmZsy1D8yzA5azRnAGbsE4O9/rm9e/+N9M/GzeUf73Qtg6uhBOHH66GjHXhYeigGR+IgwsLnIujKqJULE3jlsnep5AzbvMnuguQzSA5n4CtOrUBlqNSAKmJL+7OZc9iCO+e79KIsgHcrQ/s0oaKnyTSgWeLWSXtS6sx3rt+1NlF3CuGyvao1TzVtVrA2SHBppkHb2JIR4DICvwvNcADcLIdqFEKsBrAQwh4gmABgqhHhKBJTw9wDer9xzffj9TwBOI1NETw1ge4dFCtz9TASqSRG9hQiCrrhjMyVuWvBG4u/lb/F6Ute8OnziUC9DZTAu4aX/dNWRv9Nnt5aVSNiYic/CDpCcIvsNH4D+zUVz4J96p1InsiEJt1pJjm2Hpr4xOQGo7XWVy/jW+2amFndkH0qUVaHfMoCQfo7q3wP6FdNqJcM4ZDVuDnHR/uoYVBQLhKYiGTOaCsTvRz9ky9gHkZcrKxC4dLtwoZIqw56Cvf7cYW9HiWXi9UI1bOgLRPRyqHYaEZZNBLBWqbMuLJsYftfLE/cIIboAbAcwqopxWWF7hzKKVK3z5CWnpuoAwWTY6zjiUXdv7CwLnHPEBJwxc5w2KPuY+xUL3nr9shB4+NVWd0UHpHrVJ7CMIyK2YC7bO6g0R01MZLPtKyhiDm6GKfs473fPJMoLJuagtNfeWcbowS347rmHpWqp4/BBJYSIKH2egy45pNVKfFs2yWH2lBE4dsqI9AVuTCA0FdKSQ2J84aCCbMjuNouUVlNFbWkY2M/trOkbWNkA3oA9+whz+CWAAwHMArABwI/Dcm6GC0u57Z4UiGg+ES0kooWtrZURQBsBkKnJ1Imke4+oxjgXc9B/mBACwwc2Y/KIZCoFF1FqLrqPII3bAp56bYuzngvyd3LZRpP1+Jd1juVcAhtx45mD/8rzo7GU+iaEcHZj9N7xkBxkFladebGSg30YuPbvrztqpDtwSg6M1OVK92Iap4/jBRA8z+YimeeYiAMJifxUjkUDE+HyIalSjhrEePC4IWzbtqnVCMlhT0cXBtYx95aOipiDEGKjEKIkhCgD+A2AOeGldQAmK1UnAVgflk9iyhP3EFETgGEwqLGEEFcLIWYLIWaPGTOmkqE7JAdKSQ66ykD+VRaCzaejgju5qsDs4Fzzql9TwXvy+W68Xc3JsesHBekwuSLaFnItJIc0E+BUCXGZSrA4tZJgW3D1kGzDVl8eJ5uO8k2PyfVu5G72+Gkj7RUVcLmVEq6sirNFfN3UlpQc/JiACaWyCNVK/B0PvLIxYqZS5esCEZe9tRNvbNmTqtuiBByqnl6HT+RP4/M52a+eaOsso39TD5ccQhuCxAcASE+mOwDMCz2QpiEwPD8jhNgAYCcRnRDaE84DcLtyz/nh9w8DeEjU+mBdBVbmgICgvM0c+hLVUYI6nZJDarEJr0OBdDQXC96GWt9HZ/PZB+KxuyQHbkcK2Im8TVLytzlkryOhvhbJ+31sDqbrPpLDlDDxnv5YIldWNWFgpGqyj+d3nzrOXgEqoU42tqu9K3EcapEoLdU45obpefi+i1JZYO3be3HbC29iZ1vaDXdV6+7o/QTxC+42i8xmxTQXVQeB3YrtsF+T+30yV92DqxJd5XIcY9IAOGUUIvojgPcAGE1E6wD8B4D3ENEsBE/kdQCfBQAhxBIiuhXAUgBdAC4WQkgKehECz6cBAO4J/wHANQBuIKKVCCSGebX4YSZYVTgh4T7p8ocSxQ985d3YFU4edafZ2eXeVasolwWrW3ZKDsVCzSSCqJ6jvpSYXP0GklAaNubwzOqtxmuVpj6Od+DuXbyKyOZQ9rE5GIiMxebQr6mAaz91HI7ZP9DDp6N+Y7WPhMwwynk3qRjU0oSWpoJ3Khe158vueiXhNp0lQtoVhO/7BtV3/frmPTiCOT9bvh9fg3SxQKm5pzK9x792Ck7+4cNRXQn1FlMsgTU+p068oa2zhJ1tXRgzpAXlcu1iTHzgZA5CiI8yxddY6l8G4DKmfCGA1NmGQog2AB9xjaNWcEkOHKaPHZyqJN3sbBikGY+EkEE6ycXsmlfNTQVviaBWbqVyEu43rD/Wb28z1ILxGFEbc/jVo6vM42Fu41qqZomoDCST5GBsz3wDAThp+mhjG1GfShttIbE3JeE75eBYpeoiFpFNg5LMbfvepHTMMTjT85DE1nzdD4Na4t+3p4P3+pNtxSpffoMlYVNTAcDkkbG9L9lOfI8plsC2eaiXrmP+Dc/hseWt+PrZh6CjVEYDwxz6YoS0GeSh11Q9NVyEWF/ckVpJq+fqs7noJ1LLcWWBcYGH62a/4QOs9xekFV+DLb2xDb7MLe12GZYb6qvNJtVK8ft09Wwa28YdZubpMvJGKiRlVNI4akrCd96JU1PtP7bc7qChP5cOTeotFhhJ1/B7ZWyIiVj6pPe69oLjcMj4odHfz77Oe8tL19/Y5dje7pD+zdjZ1pl4ziaVbCprbgjTmQm2qVkvg7R8r9+/+1UAbmmyluhzzMHpyuq4X/VucU1UfQdTFoFL3i3Prk2UxyoevsF+xYK3odZ3ksq+jARAeiu5DNJI63h3W2I/XNANh52l9CE8HHz19Hod1TXZebiN4fKzr/NqMq56ys1S2dlLSDVRi8H4qBII+W0Ho7NXoRuk9VQRBUo7ARslJYdsYJJm1Dl8ysFjE9c2GKRTucEqKsGKNgzt34zOkkBbZ/z7XHM8qBOXmyQHmzs518dDr27Euq1pQ7jepnqinQsmt+l6oO8xB5crq4MOxQFA8f8m6OojqYrSdcSyTxMRlDpQH9VS1v2LWa0UfKqTnvOO4QzsH/9tOrq8Usz4xj34lz++4KxX6cYti7dSZpWdSGeJ1Rc3o1XCUaHu3RQNW0ioxSLFi3UshKSXnB541lRgvJUM8/HAMPeTSQVoImC28zJGDerHlst4mTgexf4OBoeqql3KBsXMHOLv6rMxPXe9lU+fNM06lk9ftxBzf/q4tc6lf1mEg795r7WOilxyqCPceuUsaiV7WzqxF7AHO5mkAzlZfYQHf9tEdAN7fXD/5tSYPvfuA1P1iPEiedGSUqRWMB1cY366cX3uJLiyj14pI4RISzL/NDvw9J45YWhUJxhHXPHaC+Zg+tjBfgZwB604ecbosH1NctDUSoWCXwDhhe+chmmSOTDX12/ba1QrcQciSZgefUsoOXin7IiYfVzRlMzRdKQsl8AQyB7ZDySZFIdbFq4Nx+g3+So5R6ZS9D3mYLmmE+4TDkjvlGWNshBGUViiS4v8FKHNQSLWKQf1TFkhpeTgl+fIWSUYW1meW8FjYmhrUBc0Z7QMVHGVUdU5U/399HWkXUKlWslN4Ng4B2FOO8JJUT4QSNPu/s1FHD9tZHSsqGxTfbTDBjTj8P2GGt+l+h5iySFd+fGvnYKvvjc4+lNXmXZoc43zVuJ+7+CWJsUgnb6+fOMuo/eWqurRIZt6YU1SRSfXiGzSpV6Nelaq+aiV1CpjBqeT8ul1gi4EbvrM8ey1rDBpDXQ36b8tfau6jjKg7zEHy1vU1Urnzkofuq7mVnKpO3ZrHhgyCE6qZ2J1UXDdNPGbIrdSD7WS5yyVQVmm6nJRqonHOB9rH1WcCdXkwM98mpxSXdXxqoZOl1tvJcZ+jlmRYsQXapkCLlgyHo/WlgEThw+Ixq5LBfpGpMjYHLjpSKSrVpk6BglEZ0gq5G99S9twSbVS0XMNcGMz3aM+O1nnm+ccivceNo6tn5JWBbD/SD5+JStMa18PyKvp+SoO9D3mYLmm7664DVAcAOSeDXe9nDzaU3oryeyh7ztqv8SYTDlmIrWSh0TpO0njM7D9diwAb2gsFPwiV4EgpXatYCRMHmtnaP84iWEkFVgeXMXMAby7s0r442eXrElExvfN2Ry4saXsG0olXarl4xzSjZp22yoqUYuXNWYpIQ3Sag4sGzg3W/NzTJftP3Kg0aBuj3NIawmywORRVWmusVqg7x0Tan3WybTA7K5PtuPR15L18RGeIvSjJyKMH9Yfr353Lt7Ysgd/WLCGlRwe+Mq7MG5of7yxZQ+eXLUZgKfk4DEuIN49m+oXC377BoIfQ5o5Yaj12EeJcUN5kV6HTsxl0+aAtBhjlXxZ5CE5NBUKAMqZ1UorNu5is/YSpYlhWnIwE5hRg1oS9dR2jNA2PvpJasFz011Z2WaUdeFhE/FEfLJislymi4iPm/aVHBSbQwa1kjxjgsMNT7+R+Pvg8UOUc7CTdbPS9JJhY6irmxpoj+57koNtGQUPPr5u2ikDlakY1DbV9NJyIqui/vSxQzCkfzMOnzgskVba3U82w5apOndWLbs7tag/VAzoV4yCmCQ4G4vvokp7hAYFRY+drQpVEjRVl5KbEMBVj6z0GyCAJ1Zu5vtELG09EMaD6E9bZSA69h+lJm4M7uTOTYZeK2yvdWf63IcCIzlwa6Wtq5Q60yTVltbQno4uLFrHn2UiITdG+hxvaebjHIYaUoNzY/NzZQ3qEPkzt3nHTU7YIFVk3fGb0pbrHo+NRJ9jDtY4B+06q1YKP9XJ8J6D0+qSjx+/f+JvzvCoN28ySvmI1MdNHYFTDxlrXAhfPn1G4m9paL5fCVYbr+yoOcbIkU+diJnUM4GEkYwNaWcMlP65oXjJIasfeMKV1aFiW7d1D3547zJrexMdQYNAUn1504I1YVly3Came4Smg/b9uaSkOeG8htJygy3xXvDpS/7+5Y8v4n3//YRVKpRqlZRaSZMcJNE9ZMJQcKjU5qBKcL4pKojMxwZnlTJNzKTSdDK1QN9jDpZrursfN0mImQxfPv0gXHxK0s1z5KB+msErbNOSpsC0S/DRi3eW7En9PnTMpMTfUnJQ03urE5FN8MW2nVTFtRkCeuTB7+oi4IJ/KpYcpJrG8ADMu1zZXkyK9VckmYOPqm2AZ759fTx6n5yLMJB+L75qhmDjY364LU3cYT+mtszS83nvmJJgMs1FwotrAw+kXW1mt85SpFYySA6RxB5KGMZNSLIekM0FnCxqJQ4mTzYf5vCzB+LEhyUhsGbLHlz8h+cT6yKVJyqPc6gf7JJDcrfGvQcyTAbdwCdzwUi8tnmXsS1ZzWSQ5jw19EUUnE5mTrOh/xbXrr2JIYTTxw1OlRU0VdweQxpzKWGov4Hz2XYtqqP3H26tl3WjxdkcdOIgD5b3oRk+RIFLVpgKljPYHHRHAf9drnlj1FQg9GvyP0Oa0+tLTB87OCodMbA54RbaUSrjqMnDccv8E1L3xcwhWd5iSJ9hfM6RmjDdtgTnYCCrEDIGmkX0IFnsUisJIXDFA8ujv7tKAv9++2Lc9fIGPLlyS6K8u9AHmUP6Yf/5ohOx+Dtn+kkOUUPJNlOGI62/88PTw5L63uQuxzShOLWSnmZi9OAWlqCoJ0ed944p0Xdu1z5vTnwUB+etNKQlfVQpUdIbhJvM137quOBZiuTz5Y6UVHeE3PncUs2g9yL/NqrVbubdjpM7v+Defz3joESdIvP8jfCoQ8x74lxZufmgM20Tc9BVnTZyd2Poq5/2VuLr29pKJjVMblY6SwJHThyG4w9IH/Qo35vOcOLcSsl6JprJjU2fE3d84aSwLxWx5ODpixGOi+FGcM8VjpkUGKarz4FDxvMHEdUDfY85MGWDWophcI/uysoYpJkITAFO/AuvhcV72kuJcv07YA6C8/HUuOoTx7CHsKu7oMP2i/W0XEDSV884GCPCs6o59Re3aFTd+OZd7Vj8ZtrweMohY1EoBItUddlTvW4k1PGf+4u/G8dgch00LUpTZtlEEFx474kHjk7UkTtNzt3wgDBaWMLrnGO4jeXFAp9d1Def/3UXzEmVxaq3dF9Aej66fgt3WY2X4AIkTakp5KZC16zKWCD9THOjWolx7dXHOSk8iVFdT/ctiR0DsqiVYhtkstwV8azTi/f86JFE+ng5Pj1b7aVnH+I9tmrR95gDp8eNFkcyiZwtzkGdxEKkvQ0i3Wf4t2xXZQC6Z4WcMPPfdUCiLUnASuFkWfbWTty3JBkpOXpwS6i6Me9I/2n2ZNz75ZNxwOhBrORQKBBmhEckct5KnLgdGJqDiTz7ew/gM79fmKojf0NZCHz3zvhQ94070wTbGeQUPrX/uvuVRLm8K6t/ufo+5Z1pm4M5t9Xcw8bjq4qk4dM7d56z/mibmHMJgrFVaHNQ+rzpmTcS10yBchyEiPtkvZ4Sm5+0mtPEHEwGaVm/qBF9o5Qtx2nZfadUrF0lXP3Ya+H4ySiNDWGO6DR5EroO5+LGL3NPySubd3Vgh2anMaVxrwf6XpwDs3xVQ6P6TjnjD++dEBtaPzpnf3xszv54ZFlwZm1ZCBQVW4CafExvX9ocdJWAqtb44h9ewIPMebhybPqvU3dBRIRDxg9FP8sBMXLHw+2eTAZ6IYCnX+NTLqv1yiLOJQOko2EBjyCncAjLN6ZtOEB2F8Kkt1KyTCKSHJi2L3rPgXhqVawjVufPyEH9MHwAo4pDepOiE+amYoGVHPSSbAbp4PuvH30tcc2kY2dP5UM8bzlJd+KIAdGYKGxDbcaV8VTvU9bXJUanPUS5nNowIVnnm7ctjq9Rcu6rDGH4oOZU3IqU5PRn4drkcMxDuj7LTQjn3pon3qsnrJIDnJKDRFqtEXzO3G8ojpg0LDVJZX0uM6VkWHIy6CdRxbtbgQWrLUSYkRw4gt5UJKPYKyctGyFtyq0khDU1QjCO9M6bOyvCFTNhIi6InrP19hRUBwNT2m/bqXgFogSTV5//iIHNOJRxudSdFbg+mwoUHdH6WmvMCCeNSD4z73QKoT2Nk37MaiW+qeh8cUbxf/KMeGOj2/AAj4ynWv20QdqTOShlaWk6qRpWPfZ0g/SowXG2WM6ZcFC/JhABu9qTkrgrPMEU9KaOnduMVBJgWCn6HHPgXkkxwRzicn6nHLaj6TTnhPmSZHCOPgGnjgp002fMHB+3pdwPxK6k+gRQUyTYJkeB0kSfJejMGQwSseRgnhpqGgzTMaHs2LSKN33meNx44fEYpuyuXcRdlapU1doxUwIvpqxqJc7mYJIcuGdGlDSCJuYFzB5vNvUfEH4Om58AACAASURBVDJwEbyPU3/8KADgoHGD8a1/mKmN3/DD9D5hZnBFRa2qwpwVNvg02cjUPlNBbRpzuPHC4xPj0uvrNoeyEOjoKhvTf3OurDoTK2hruEmTrtU1o+7w2XQiBcLgfk2pM7CdaiXL9dhdN32tkceEOpkDEf2OiDYR0WKlbCQR3U9EK8LPEcq1S4loJREtI6IzlfJjiWhReO1KCt82EbUQ0S1h+QIimlrbn5iE1eagTWaTATZoR5k0AD5w9ET8z+fegX8M8yXpTOSAMYNwwJhBOOfICdF9eh1plGvWPVLCP0tCsMzhD6G3SYHSC0HX1cp6OpG47fMnAgD+30eOwjlHTEgFW0k8felpuPq8Y/HROfvjt+fNjgid+wS99MIfN7Q/3jljdMJrydWOylw/e8Nz0ffpYwNbSaVqpbJQbQ7JZ3xwaIfhpC19556YF8Lh8aZgysikYVvOSVW1NHnEwJTO2dfvXUp4WXajZrdos1oJQCQtjRvWHwLAFiV6W5cc3jljNA4YM8joUJD2VgIO+uY9eH0Lf4gOJzno49QZZZPl7E2VQJuex+D+TakYDpdB2hQRDaS1DSp6muRwHYC5WtklAB4UQswA8GD4N4hoJoB5AA4L77mKiORs/iWA+QBmhP9kmxcC2CqEmA7gCgA/qPTH+IC3OcSSg7q4jRk1kQz2Cox0hOOmjozu0QOF2jpLGNQvaeLRjdYyCE73SFFFam5ynBieUUxITzreqB7s9h9fER8tefT+AX8/aNwQ/OLjx7AGaQAYP6w/WpqK+K8PHoHTZ46LPKRcJNkU1KX+PsAtOZiOzpRQzwj2QfIMaV6t9LW5B1vHptbX1RkmpwY5L/o3F/DZdx2QCp6TBEsl5lv3pFNkmNJLp/oMx8YRHPn802NN11UZnskH/0unzcCfLzoRR08eniKSnFpQtYeYvJuKyhrwgVpNlzJiBpJWoer9q/2Z+u7XVMhsc7CpnSJ6wL4ra7M1hZM5CCEeA6Arus8FcH34/XoA71fKbxZCtAshVgNYCWAOEU0AMFQI8ZQIVuDvtXtkW38CcBr5bocqAC85hKIrkovR5sq6fU8sRnK7XXUCCiHQuqs9RdhiySG4X+76TYFOQohEvzqWb9yFjTuSHiQRs1KalPp/k1FXvc8FQrgQHGvWlkhOus8C7kXl8tY46/DxmD42CNbzISOxzSEuKxYoOo0NUA9bSrfYXCSccMAoHDV5OGZNHp4iJrzkEHsOlYVhExJ+/vGZNVEZF5h4xbxZ5h+nthc2yEkOUV4hj8R7QEyoVTvTA195N26/OIgfaCoWcOyUEd42B+nJxvWZUis5dPlcUkCThCPHpm649PEm3yffZ1OB0KldNAW0StjVSmF/nKTak9RKBowTQmwAgPBTHgg7EYB6QPK6sGxi+F0vT9wjhOgCsB1AOkqmRuBeiYyA1Xe3ttxKe5T8NNyB5Kpe83d/fx2L39yBLbvsydHkwtXFXFX1YTP8Lt2ww3hNhXTZtYm+/jl7Qk8fBylWCYCOyz5wRPRdGk25Bf3kJaemjPUA8A4lqEpKcL4ghfHK4RGA27/wzlQdffzLv3cWmooFDGppwu0Xn4RDJwxJesmU+cXcVRbRORn6AVAS67ftBQD8519j11/u+amM1QUh+DakBKBLLx/51VPB57GTovQwRPHc6EokihyMoyYPT9zP2aO495eIz0nZHJJSjdvVOdnM3o4StoUbqi+dNiP6DSqarMwByrX4D6mGDcZYQKcmnbjsMXaDdHDNZbeoN2ptkOZIirCU2+5JN040n4gWEtHC1tZWrooTtpQEpF3nPEEifatyzOLR2qJQ7y0Lgb+FhtNNjF84EP9YOaHSkkPw6dKnX3qWX4BMgQLCZVtovl4wgeeNYCUyAPjUiVPDPpOMd95xcTS27oEjBPDzB1dAx37DB7CE9MCxg7QS/0WlquzkLj2dBC/41J9/ehdMKWLCjfeBVzaiVBboKpVRNtgl1m4NmMMARVLiiIX3ewqlFW7nLVWRx4SpSaaNTj7Pr809JBoHIXZycO2OI6lSgek0QZl3q0NrUzIuUzyBbhvTbQ7vuPxBXPKXRQDi/GK60drkbQYkd+/qNamGBdJqpRUbd+Lrty1K/U4VNsK/Y29X2Le1ibqjUuawMVQVIfyUjvfrAExW6k0CsD4sn8SUJ+4hoiYAw5BWYwEAhBBXCyFmCyFmjxlT2cExrOQgJ6xmqOUWtqzbUYolB5tt4hcPr4rcT1W3OLWOHFQsOWjEyeIto4JL+sae3kYhobA15ys5gFfVAcAHj56Ib//jYWGfyfF/8bQ4S6zuUVUWAm9u4yOauWddzYEosuunVm2J3pPegyQmusqObyuphuAIv3Ra6CoLo11CMvoTD4ylIu5nesc5hCoePsp7cFiHcPqhYxMMCQikE9WTK3JldVCv4CAorYxdK4GE8c+/fgrfVSQlANgbSuicO/FH50zGnV98Z6K+buvbpqhhm5uSKtZIUlSGZEugZ5plTVo0+6eufRYvrLGfo25LxS0Ziz6Wfz39IK563VApc7gDwPnh9/MB3K6Uzws9kKYhMDw/E6qedhLRCaE94TztHtnWhwE8JLL6I2YB07LUoeo7HY4QyV29yZVOx68eXRV9HzskaTzU3V2l3tLmyqrjs0o0NTfe6y6Yg8+/58BEKmnprWRjNr5qpcggzTSVjBlJEopmpQM9sKcsgD8/vw6+6OjSO/fXy3I++/pjlH8/8MpG2KDr2MtCsB5vM8M0JjJAjHtvM8YNwSHjhySYu+8h9MaxIc1IT5o+CoOVQC9O/ddULEREOalWckgOlJ5jJlWtEAIL39iautYvsjkEf+tziusT4NWcLU2x9BPUSUMvU5+X6fk3FwsJeuBDvhxaJwDADU+9kfj7Y9oxAPWGjyvrHwE8BeBgIlpHRBcCuBzAGUS0AsAZ4d8QQiwBcCuApQDuBXCxEEJusS8C8FsERupVAO4Jy68BMIqIVgL4CkLPp3pBnTSzp4zA8/9+RrRz9bE5FD2ZAzdxdTWEqh99dHkrbg7VGilXVmVhfPKEKYlrQ5UYAS56ctroQfja3EMSBEgeNmMjNllcJNV01yr0Z2livDozzJoLf+3buluj//1yGKr/va6q8d2d64TVZGwuagzJ5rte1piNDn/7JLE2h0vPOjQ5NkPajihAELGk5z5gyE9yKBTM0ueZh41P3KcSXlZFFX52lQR+qmQ9BeIklHr+JbWV6BCmr7wLALCjrQtvhvYf0xh1tZIPD7e5sgKBd6OaTQAAxgzx80yrFZzpM4QQHzVcOs1Q/zIAlzHlCwEczpS3AfiIaxz1QL+mAkYOilU9us2BCyCTE9KUfiJqi1m0prw4AnHWViA2kOv3lcoipTq64KSpSj3rkOJ6BaBcsus9fWmOeoiMDn2XVzYsbP25cMP6v3PT9pTmIqGzJLBtL0+k2gzpw1XIvjlDqYS/Xj/Jlkw2BzmvJBE2vbcCJc/nZgPYMjBxQOAPC9YkyvUD7AsFYudFVKSola578nVrnwVyOyoAsAZlymcVrwG1fbPkcOdL6/Hrx5JpQuQGIJYc0n1Kmi3jZgBg0brtmDh8gHGMulqJs1ukVKfdbE/wQd/LrZTYzaaJtSt9hq/kwC1Z3+jGdBBcrG9VicWxU0ZgYL+kSsAHkujYdjje5wRAevqkG7tEMZATUWJBqETNR3Lg3oWspleXf3PnN+vg1BWy7Lfnzcbujq4MnluUMmDy7tDBp/T2MZ1eF2Syjf/mJD3foChpG9ITNqbaIz61iiwpUDaG5LOLLpBb1uOy8fJzNCjjzhXR3bqj+aPU4cbRErky82NLqZW0611lgX7ae3JJDt3sqASgrzMH7YURkgTMZnNwSQ5symtdl60Zz6I+UpIDwnpJgq4TBn/jJKUYTbqOX1tSJaCO69oLjsPsKSMwpH+s8tLjHEjhf/pz4UbFEYL4DADt/gwLS5XKdJw+c5zxGgfSiJzJIC3fm9xtmp51gShBdFhvpRr7vRcVyWFI/ya8K8yVFHn2gDtrmocu+RhB5gBJtS0AmvTJNgXAfrymfpqjajzmxtvSVMA9izZExnEdzZpaSW/DFlsCALfMPwGfuvbZRPs+Ele90adzK6WINSUXIK8flczBrrLg1o+er8hkPDMFwZU1w69eL9suUmTW7fNtpVUCg1uaEowBSLuyqrtPncCxkgPz22R71fgvyK5tBuks+Yt0gzRHSOXv1d00uXoqc6jmfUnG5ZIIC4qE11yM1a6cZ4+7zzTR5xm/+x1Ga0BVK1lsDro30OUfPCJVl9tFcDylX1MBF930vHFszSm1UvI6JyWoxvzmpkLquVbhe1Az9D3m4CD+yevp+30lB24V6TudSPepSQM6sZR/lsoCv/v76kRdFf5qpWDyyc0Ol4HVf3coF7xCXLl6BX3XZ+5AsAnHzGNISQ4VxDmoCzj9/P2fq09upcggHfZpap2QDHqsxlupSIGh2fVTioV4p9tVKsdHaioGaV8e5ctHfNqU7//Xj8Xef+z5ItH7TDYoD/hR2yszDE99fzLBpN7N/f/6rsTfehCcj+SgbkKbCukzJPQ2TLnO6om+xxyU72nJQfdWMqsEnAZppizVXiQ5xGATo4X36bl1dKKeSdxHLDn4nt3AIlJR2cdBRIkgQFvzmROOVbHLKjDExOeX3/Olk1Nluo7dFMMgNwkLwjMwTM86SKRoVyv5oqW5iI6usgdziNVKpbKI5pj8XZzBetG338u2lckGpr3EYQOaI48hICb6j6/YbG0/lhyS7aU8BZU+1Z+j3vWZk6cBSO/i5YFYEs1NhUTwXlpyYJhDWWUO3Pnd8fcffujIVDxHI9D3mEOCiOk2B94wqUISaqfNwaKSivvzXDyFNAED0moqb7VSqDqIgu4s4rlzbJS2hXB3613YjJq25HAc9NqcZ5MJXCoIF0276uPHsOc06EQuUCuZd7cy2MnX5lBNEFRLUwHtXSUcq0T23vmFNMFR0753lUXkORflX6I04R3MnJAm6wLAuw6yB6wSpb13Fn7z9ITHELuBsUi8+lrh3MjlNFPdmBPSLWOP4k6Day6Q5sqa7JtLpZFgDkVKLRlVcvjAMRPRHehzzEElJcdPS+bg0Xd+NsnB6a3ELHgjgXPqW4NPXY9qsk24IL2yIldKdpFlsV8E50er7evQx+br2y9hVUNpz2+UZ6ZSIGbQWaKsjSOh5KsslfksujpjNEsOFKmVfvChI/DBYyax9XzQ0hQEsg0fGLtuDx2QJnQmyUEuG87GZJor8taBjmSJ+vG8QT98W4mxsow3+NRzF+nZYFXngbFD4/ly/LQ4Ij3yElTmxi8/cWyqz+ZiIbG50Jcz5zmlzrcio1aSl7/9vplWN+t6os8xB/nifv3JY3HhO6clrpGmMzaF+hcLVJlBOqXGCseklOlHhMo+AeasBkMOJhcKIRHrsui8s7hvCogEcxg3tD/TnjZWrYN7vnQy/u3MIDU2Z5y0MeNqbHfSa0rdabojf/mHoxqkg+SBIhEJLqFHTXMpToJ+4t/NZWTNAhkdrB9YnxpbaJsQQgSSgzw/O7quZy42t2XS/6fqIb279tlMsN5K4PvUA8jUd1UuB3mXXr/8nES96HhYZT4OakkzusBbKbbT7NJcqPc6mMMk5WhV/Xojs7Dq6HvMIfycMmogo1bSDvsxvJciEdo7K1ArMWosQPdASr+SOKI22Wf67IZsOl5JBLmlq0cw//yjRxvaChaXSrwnMsd/qkM77ZCxqeuHThgaecZktbu6VPE2TxjOldUlRZjVQLHh1pRhV+1TwiQVkWJzMDEQX7SE6eL14yx1FAuBWkk+giZt90wEHKAk5rPNudiRwp5aokCUIuYp8xzTDUs4mT6P2X94OrpYeVfBOSnppri5we3ig5TdQX+fZ7yabJLDPV86GS1NxRQ9kWvbd5NWD/Q95qCIxzoCfapS19BGsUDOM5PZ7XhqwlM4JnXymZmKvnvWJ526UL9x9qH43adms0OTbqVy935G6M9vwvtnTcT7wmRxOrrKAk+9tiWTDcZGXAHe5mBz47R5Jz34ykYr85B9qoZOn50uW66oJU1HvnJlnM0nGBtFz7Va1YK8f68iOUhpQh9bSYiIOMmxzpszGf2aCjjr8AkYO7Q/TjmY9+TRxw8ADkEsISHFZW4GajNIq+/w/zv5gFS95kJszymV+WBF7nhYY+JAEazjvy3dqNQNPjlpTUojsg89hkJu3HLJoYEwHSIPyLwy8SQ1qReaClRRhHSqP2bx3LM4HcEqx6pLDvpuQxU6jt5/OE49xED0Q5uD7P/yDx5pH6dllsgAoqXr7WdJJNe2aaecXowStt28jfhfeP3CxKl96XGlxzKQyW6rwrRgX1izDaWywOrNuyPiZGP2EiaVkWqQNjEQX3AH9Iwfxqv/OGeF6WOHYPn3zopO2pPqJhvxkldsGUhlny6GzL0nLjmjHI86X7gDoiaOGIA3w7ToQYJEs21IHf7ezjSh547iBWKGzLqyWjYPAHDP4g3h72EvNwR9jzlEkkMazYWCV9BRwcPm4KPikZ4geppiU1s6Q9JdCm3J7PT2pM1h7JAW9nQuIPbz9vGCcklS6vMwNSeLn1q1JXXtwDCtNACjFGPCzG/dZ7zGvab9GLVY4h5D+ZPhuJ9Y0RoZKDnCn2IONptDKbvkcM35aYkxCt50qEOLBYSSg514SaZhmxq26HMVqvrM3Fa6bON2Pq07kJQcWpj5PaR/M/Z0lLD4ze3YvKvDbtxW1tn+I/WzQ8wSrzSC21xZTV578oxsb5fyOqDvMYfwk3vmTWEiNwnTbkaVHOa/Ky2ymtrX4ZujRtIXnTmkXFu9mUMgAr+9u8NKdOQ1nwnqWtxD+seeMTbXTQC446X1ifJffeIYnKLYKfTDaKoJ9PZ1J07c47hFIHYe4CWH5N9mm0NsNM1ic+DVRcGnS/0ng+XuCyVYk8QSnZ5oeX6xzcFt4K/ECaCdmXOxtBK3d4CysZCQv/Mffv5E8LdF/SfH9q+nH8RmRlVzn6mQ74wLXnxx7bZEHR3ybPDutDn0wdxKEXtIXWsqFhI7YNOkLhbiYLnPZmAO+kLyjUuIVQLJ8ejienJ3bpccXt+yJ9qdmNAvPBzFZ5xthrwzEmOHxCoMo+tmSMB2tiVF95kTtNO+tPuqSZ+h/7SbPnO88x4fni519pxB2t/mYL7HBk5LZZI+0/cG9f7tTy8H/Ro2Dz6SA+et9I4D0ycAE4BNO81SQNBPuqOvMycfcgxpOHOUaqGQlAhMambAbRyW96ZPvQvVSsz8/H14VoNJpTg6PBisks1LrdDnJAcJVnIoENoVIjfG4C+vLlTTy/V5qVmC1gBGctCYhdqcbafpO92ySA5yR3rL/BPY69JbBnBLDjqx1+unUg04R2eG3pZ+IBMH17sVIn43PsFbejAjNzabhKenVrBG9jskPF2aNUoOEXOwS6iAFkTWP02oH13e6jxylJOy5x4+IVUWu7LGv9MkFSQy6DIbQf0gKGP23LCe7rIq57FNcjKt087IIG28te7oc8zBZnNoKhB2hLvW//rgEdh/1ECmVnKy6WcvSHAvVSWQAL8T+d77U0deGL0e9EmnTl6bysrXA0LqTG1MTJ4R3dZZxqzJw3H8Aemdod6nMU4gMtBruttU6gP7uG34xtnJw20KBUoYoE0EIDkAvniUcjZIFoO0qccEs7eM69gpIxJ/2zx7fCUHCZtqFbC/C3nJpTLygd7Pv515MMYNTTNyWU/tk1sLBUqmAdGlckBVK5WjezjI9v/l5hcS5ao3lAmm9yrVtLm3UgMReyulH7q6O7MlumpKSA52QqfiEk0M9k/oFtTTT97atreTrQfYiZyvhqK5yc0c5BnH7Z0lo2EbSBJAFzHUN7f6AUf6cLJolfRzvAFg9tQ4Ut5EANTYDNP4r//0HAABM5OEgdX/kx8BJk/JwSdgTJa5HCn0d/365t2Gem5vJd+zz32gzudRg/rh4lOm86lJws9EoJ6H5MDZzOIEiaHk4FArqe7QQCyl2ZiDaW11OfpsBPoec7BJDsouz6aWSezQPfTFAHDc1BEYyojUPpB9qD7UAHCSpr/13Wn6ekD081AryWttXSXWK4Tr05xuIblTk9DTL/ik+DaBW6dNHhLXN86JJQ7T+KWxsixERIS5Z6K/GjW6XIVazToftUs24+prrTyxl9Cf/XZtAxK3x/etwjdC2gdqP7Z5FgfBudVUquTAMYdYCrFLDqby7/zjYQDs89PE9DsdfTYCfdAgHXzyNof4RdnUMj7ERCfOnK81B9aQzZQ9felpGDGoWaunSA5WtZLXUCKViM2LUj6yzpKwH7Wp9GkiKJyOGkgbdatZL5zxWiWmLnsIYA+CC/qIbTCc5KDvZN9rCEJU+7Slz9DHzM3JSmPoPnA0n/RNjsc6z8LPLHmrTFB/U4tlLUmbg/sYX0pIqJzqS/dWMs+NdNknTtgfpxwcSJs2M49pcxmrlcz31htVSQ5E9DoRLSKiF4loYVg2kojuJ6IV4ecIpf6lRLSSiJYR0ZlK+bFhOyuJ6Eqqo6It9lVKd9HkIREE1+LHZlLfqIv5sP2G4scfOSrjSGNwC3D8sP4pwpMwlNsM0p6PVxJ7X/uFntxMhfqYzMTVb6epv7uDxw8x1EyDa9nnvdvOvNbHJYSI4gl0O5Pe1nnvmIKxTC4qHZztwjQefoPh987VeqcfOg4nTR/N1pO/wdasb5yDKauraVw2yUFWc3nPFQuaEZr5HfI3dpYdkgMzZz514rRo42RLtW6S8HtLhPQpQohZQggZeXMJgAeFEDMAPBj+DSKaCWAegMMAzAVwFRFJ6vZLAPMBzAj/za3BuFhExx0yz9xHXRRcc/ejEudLzzrUiwCYkOUQHwkbQffVY8ZGRz/CZLM5qKvPFDDHqZVW/9fZqXqqC+/PP3o0fvHxYyz9amDWadHjvauPwLWDFAC++j8vAnCrlWzv9l7lvGfOJTYeT7INn0ywPrDNc5+5IS+tedvuMv1bJmhPh69aSVYzHekpoWafDdpP/454Ptr1/9y908cOjjaIJYvoYHp+kbeS8c76ox42h3MBXB9+vx7A+5Xym4UQ7UKI1QBWAphDRBMADBVCPCUCyv175Z6aw7aH8fUrN7keqlBVLFUm1PS+398g7clsCvFO2NxW/H1wf/MOUK139yL+kHtZRxoA9xvWn108RytnEpw0fXQmWw6Xh2mDEmlrejQJ1ZOh7Sj9R1lg447AjsCqlRL2F/NYzzkidtXksrua2vDJ5+QDmyordmU13+/bpc98LCQkB4taSbqVOpiDem5F8He6TspbyeHKmro/8r6zDoVFT7A5VMscBIC/EdFzRDQ/LBsnhNgAAOGndPOYCGCtcu+6sGxi+F0vT4GI5hPRQiJa2NraWvmIYZAclELbYvLJc+Njl+Awe8rIVJk+QY6aPJy9N6mjNvfp8liRiCe3jTnE/Uw1uP4CfuKxbEsOncv9E7SVHqMvXKmvTe0l3oGH5CDBSw5+c+PE6bHDgU1y0F+1T0ZgH9g2GHEQnL/6Unpz6fCRxNWxcFHKcZ/Bp8tHIUi7H+/ordmQHemzjRJF2GQlx7v2hqysJwkhjgFwFoCLiehdlrrczxSW8nShEFcLIWYLIWaPGWM/XcoEmytrgujYJAePhdZUdNsldLzyn3NZ/blqLJszdST+ctGJ7P3q/Lb1uWVX7BJ79P48o1HbsBnU1Edh29F5PYKwjpQEfnOeOaushCnOZMbYdMqEk2eMxrmz0nmZEnEOJntCgjcY6oTlKi1w2Rxs72mAYnh15coytR/XMd6egNqUbYPhIz3rd7/bcCKc38Yh/j6acUc29WlCkQhvbtsb/X3p2Uy0tTzro2Qn1PqmJupDrh+GU7U0FfDZd/PZFYDYoL7PGqSFEOvDz00AbgMwB8DGUFWE8HNTWH0dgMnK7ZMArA/LJzHldYHNldV3R+ej5lEXlm139bN5s6Lvuj+/hBq1O6Bf0eI+6yc5qMP5xcfM+vqmiDmYuUOiT6sR3Hgp1VZHqYzJIwcYT3TzeU9fOHV6quwHHzqS3YGrnmQmYq32Y/otkpioarhq1Eoqc7AZpPUrXJS3OmeOmDgML/z7GYa2/OZtk3Z8KAdfacVns5WU6i1SlCcX1OuNZuZafI6KtDnYNw5ztJMlbQb5zlKZdeB4/GungAjYGsY07ZMGaSIaRERD5HcA7wWwGMAdAM4Pq50P4Pbw+x0A5hFRCxFNQ2B4fiZUPe0kohNCL6XzlHtqDlviveSuyfxofNzqVUJp18t6LIwC4YAxQbK5kYPMuybf3EoqORnUz2YnSO+EbX3aGVJ87cAx6cyWgJJDqqtsff4+tiHu95tsIgnm4NgdAhabQ/j5QphUDTColVQJz/Ke1HHZXVmDNmZNHo7XLz8HgxjvH5UYjh/WHyMM8yj5bI1dKmcdmOv42sp87CHqerKNy+Yxl+jTY91FKUccZ2pEXk2acSHeXCXLr35sVXCeCjOGof2bMWxAMx5eFqjN123dm6rTKFQT5zAOwG3hxGwC8AchxL1E9CyAW4noQgBrAHwEAIQQS4joVgBLAXQBuFgIIZXfFwG4DsAAAPeE/+oC22E/SYOuuQ3uZCcd6kSq1n4BxJPZzhyU+p6SA1l+p9wVujK8SlgJmPL9qo+nz+FV2+roKlvdG31Sk3PlJi+X5qK7PZW4mvT/cv7cvyQOVuT6LCaYONtU2E980So5eEyh/k0qo/Fry2fe2vTpvknjfMbf7OE+Drg85mL4SBh6VLnpHch6z72xNVFeNDCH79/9athuWiIvFIDhA5qxbU8QfGiKUG8EKmYOQojXAKSc94UQWwCcZrjnMgCXMeULAaSTCtUBtsN+fInroje3O/vxVSv5isHSe8FXcrDtjFznZEuULcZ7CZVQ26N442smnXGUYLBU9jbAZnEvNO0qbtpobgAAEv9JREFUE4THw+hoTpmSLuN+h/rMXtmwk20L0IPz/J6tCWOVPET2Zxu3ZY9bid+VCeqwuKNhJXwkB9+8Yb6SgzXKWhtXu+Mcb6OakQhEZtWbGpUtU64UC4RhA/sBYcZkX2ZXD+TpMxT4eiv5IOHKapnMvt42Mo3BKF+1kmX8vgfEy0ntk3kT8IviHdLSZLElBJ+BWsmPGJqIJvf+THWbm9yM3CcWIsoq68gTq96/RcuXpcLlWRX3G3zaCOPIgfG88VX/2SKRZQyGnl5dhZrk75pPHWesl9WTaq1F1dLsSUxfWrfNWUc+prtetp/K5lrf6npTg/PU53PeCVMABO982IDYNftz7z7QOc56oe8xB/mFe5+q5FClIajJQ1WRBXIicfpkCV8dr6oatUoOHsm/fBmqrDd9XNqLKKpTUCWH6nbKWR65jwrQ1/AOeLhRetppsqZ0t+0yCwXCp0+a5uzTN9jsTQ9d+KadfM6odJ/Z1sfzmvpGhc4gf/1JXoXpyjEFpDcTpgOtrMxBC7ZTo/8PUtbC188+FK9+dy76NRUwQPFwG9Til3anHuhzzEGuXJfNoVqCnky3UFVTAGJ3uCxpFEzIqlbySbznGpvPyKI4AeGvFze2leGhJ6U809jcBD2WHOzwdReV11zpJWS/LqYlu/JlvLZ8YEdOMmctlnClsNDHJTFphP2YVhvBVJnD+47aD2ceNp6tN8QSsGmCKQ2ILfizWCCU1NMlQwbzj0fth0+E0gIQzFf5vFUeVIuNZaXoc8zB5q2kvgcfFzKbyOdrQPPdNMl5mSV1swm+aqWPHj8Zx04ZgU8qk1iHrwEzSxAcYFep+PxOXfL7oCGBHAB84vj49/kEOplcKX3fZYI5WAh1lL/I0Z5kyi7nBjkPBzTbjP3xdxtx/9rcdFyADnn/cVNHWOvp88b1fgdYmJbqDm57HrY2TDCdS2E7ICqd/TX4ftzUEca59uzrb0fffVWL9UDfYw4Wm4PrNKosSLqymiepb7Zpudu3GRN9U1d3lf0kh7FD+uPPF51ozQtVaeZQZ1vWlOnutnSC8/lTzIzcdKiTqT2X5OB6D+OU52l7ZtEzcDw7OTbXLlOeVDaWOSRHYpYSfW9LmDfQEJOT6C9kDu87Kh14qCJ9HoW9XVM8EBCo1oaGO3lbM4dbzmsxodMQ72OT7AqFZJoO+Ux9ggiB2mgdKkUfZA7mCOk/PrMmW1sWBUJCj12DQJZYcrCpgfyYQ0KtVOXsU+c4Fw0c1Qufga+KqtrzKPQa08faM7de9oHDvdN/GJPzhZ/y8ZpdZwvKd7dayfVr5YbBtcvc3R4Yj20qlSMnxczBZsOwSbAS8oQ6FyHW7Wiu9+v6nfI0xztfNsfSfv8DR1jbkJivnBE/3rBJss173eYg7Ra+AaP7ZBBcb4Rv3nm5a3rfkeYdkUrcavF+JeG3eaRYApkTqEV+fYmEjtqSPkNEnk/mtpLqlurUSlkPmPn48VPwyL+d4lXX5co6JCR2N1x4vEdbNsJSCNt1EcuQ8ToYvZxDzZ67VhsDsDE1iffPmohnvnEajtnfrlYaMVA/l4Svd8OFQW4m302QTRNgkz5UqHPttEP5czdaiua2ioXYW+m11l14eV3gBl+tfa4R6HuH/YSf3Avw3eE//+9nYG9HyRhlCvi5PqpwGR3lerARTekCp58prEOumds+z+doygL1t6kuoTrimInqjds+wo7rsJdqYE7rHZRL9cPM/YY62/o/Zx5kvCa7GeQgZHI8LptDOZpDfvPc1p6P5EBE/3975x5jR1UG8N+3j7b0YZe2W6jdtiy0NKyllG4Ba1tEsEB5FVKJKKGFGpGg4jOK4DNIxAfEIibYaBvwhaI0FkVJNQJBUUotr1preRgtQWqjQlsCdeHzjzmzO733zsy527l3ZrrfL7nZ2XPPPfPNnJn5zvfNd77DxDHpqeorr4l3njClZr3+F/7ZjW1SCQ8zblY/VFtYHRFlFw1lPfXG+yPt+qdgz4uhpxwSJnb59smI9tbUld183BBRTuquzsYaZczwNna/2pf40Bw7sp1HP7s4NYX16/1+zwO/CKOj0KSHiY/lEL1fEvPneHTUvgYqh7TUHuGI1Wf2e/eE+NDe0GefFL4M9I940vYXujd8+31Zb1fsd8PbWjhn9qTEOoPhL9edGeuO63fbpcaDBcxNSCoZ5SvLZsd+Fw4Yp433Uw53XjGf3oilFFgO1b9JTMGeJGwTGXJupaSV4BoVNpbFQOCLF8ziyAmjmNyRHObXMXJYqnvhtf75CwcumO+Ld5+w2KjllnTzeCmH1/zCKAdDWihoeH59RtdJXTXl0JGMHt5WM2NolD5PZa91KodayehCRIRb3j23fynMrBjW2hI/cnbFvpbDqouO96p31uxJsd+F91JnwrnYz3pubdnv/mtpqe0GS/IAJKXIbyZD0HKID1dqlDmX5DIKHzRplsjSOZNZOic+HLMefGY++xJ9UCbl2ffZp28qDh+x3zg2WYkeCEmROqFsIn4P4aRrbtTwNp78whmx34eEWXPT9he+k8pzAZk0fBap8n10+iw/Csluu3DAkjSXIUql9VY5QzquXpS4kNlmM+SUQ0jaPIcsGZmQ+XTRjE4+8LbprFzY3Zid1+C1DN1K0RfkYxLcWf3KwXN2brJPNl2uk44cn15pkCQvjSmAer/0zYIwOiftAfa6h2uvyAxEg/k9PH1SaRw3pSOxP8OFftLe+4RU3lPtrS01Z1YnRTiF9Q9k3fksGHLKIWmeQ9Z8/tweNmx9IbFOa4vw8TNmNkGaAQZ8zwfeVpJ5HEX73UrxdXxfSPtGZU0YPZxde/xSOPhwSHtr6vKTodS+L32zILTY5nQl+9h9ggKKTNe4INQ4LmqokrSIqu3XL0l9DoRWme91XmkRHDKsteZkwqTUJKGb8MxZtWd3N4uhpxwSVoLL+p3DpQu6uXRB8ywCX17P8p2D5znzS+IXcSsljLz7fLWDtwPCj/s/cQovuQSIcYTHkJbW/foLZvGTTTsS6/hyYW8XR3WOorfGErNRwodlUdwW9TK54xA2f2bxftFAtXjvom7WbX4ucWVC8Hsn5Ps+J6Sy3oj21prXQqL16f76ZI5tJEPvhXSC5XCoy1y5NiGD5MFEJtFKniOq1z0sh+j9khSR4jtPY8ms+BeNg2HimBGpk+mSUlhHufikaay7ckEWYiEiqYoBBlJGvLwvPpNq0Tl01LBUy+fas3vYeO3bM9lfOJDyHQRVDmpGOsvh3wnZdytZd+UCPnjqdG9rpVEMPeXg/ta6vkLFMakjPTa7zHx7xQmsXNDN1HHpaSPSCBXMYQkpGWDAckic5xC5AR/cviu23uFj/frnc+f2eNXLkiwnGGZNOPErLSHeyUd3ej8Mi0pWrrN6LYfK3R7iLIdn/rVnv/IkC/rYrrF87PTmupprMeTcSktmHc7Rh42uaXJmGcVTZKZPHM1nM3xwfuuSXo5NSZHgk/67LzLqTooIGjOinb/dcHaqXHmPvIrGRxcfze5X+rhgbvLchNtXntgkiQa4/OQjGzo3ZbDU64KtfNEsAtt37uEdtz60X3nPpPQJknkz5JTDtPGjYie0nH/8ZL5677bEmGajmri0yFFGuId93EI/sP/Etes9c9+k8fA1p/VHnDQT39w9zWT86OHc/C6/2P9mc81Zx+QtQk3CBY/SZtzfdeVb2LX71aoZ4fc88c+quqsumpPrCm++FEY5iMiZwCqgFfi2qt7QbBmuPOUo3rOwO3XOgVE/pxzdyXXnz2LZ3Pi5GpMjefxnHp7s2/clKaNsIzkxZca7UQ7CtBlpL4fT8keFPPuls0oTLVYI5SAircA3gcXADmCjiKxX1T83WQ5TDA1CRBLXhYBgPsj0iaN5auce7xfdRWX6xPi0GEZ5uLB3Cm0tLZw3JzntuA+Js78LSFHuwBOBp1T1GVXdB9wBLM1ZJiMHfvy++ZkkBDSMLGhpEZb1dg16sPLRxQOJFd/ek22qkUZTCMsBmAz8I/L/DiA937Fx0DFu1DDGJWS7LTqrL+n1TiltHPxcddoMrjptBs/u2psa0Vc0iqIcatlaVXeYiFwOXA4wderURstkGHVzusfLeWPo0T0hPqtrUSmKW2kHEE3i3gVULeOkqqtVdZ6qzuvs7GyacIZhGEONoiiHjcAMEekWkWHARcD6nGUyDMMYshTCraSqfSLyAeBeglDWNaq6JWexDMMwhiyFUA4AqnoPcE/echiGYRjFcSsZhmEYBcKUg2EYhlGFKQfDMAyjClMOhmEYRhXiux5r0RCR3cC2Gl9NBf7u0cRY4MWM6mXZFvgdQ9b7zPI4s+yDPI7TrqHG1bNrqHH79G1rpqqmZ7ZU1VJ+gEdiyv/l+fvVWdXLsi3fY2jAPrM8zsz6IKfjtGso/+O0a6hxx1nz2Vn5ORjdSv/1rHd3hvWybAv8jiHrfWZ5nFn2QR7HaddQ4+rZNdS4ffq25UWZ3UqPqOo83/IyUfZjMPnzp+zHYPI3Dl/Zymw5rK6zvEyU/RhM/vwp+zGY/I3DS7bSWg6GYRhG4yiz5WAYhmE0iFIoBxFZIyI7ReTJSNlxIvKQiDwhIneLyBtc+TARWevKHxORUyK/6XXlT4nIzdKkNfsylP8+EdkmIo+6T1OWlhKRKSLyWxHZKiJbRORDrnyciGwQke3u76GR33zKnedtInJGpLzpfZCx/KXoAxEZ7+rvEZFbKtoqfB+kyN/0PhiE/ItFZJM7z5tE5NRIW7k8h+rGJ6Qp7w9wMjAXeDJSthF4q9teCVzntt8PrHXbE4FNQIv7/2FgPsHiQr8ElpRM/vuAeTmc/0nAXLc9Bvgr0AN8BbjalV8NfNlt9wCPAcOBbuBpoDWvPshY/rL0wShgIXAFcEtFW2XogyT5m94Hg5D/eOCNbnsW8Fye538wn1JYDqr6APDviuKZwANuewOwzG33AL9xv9tJEFI2T0QmAW9Q1Yc06KHbgfMbLbuT44Dlb4KYsajq86r6J7e9G9hKsLTrUuA2V+02Bs7nUuAOVX1VVZ8FngJOzKsPspK/0XImUe8xqOpeVX0QeCXaTln6IE7+vBiE/JtVNVywbAswQkSG5/kcqpdSKIcYngTOc9sXMrCS3GPAUhFpE5FuoNd9N5lgxbmQHa4sL+qVP2StM6U/k4c5KiJHEIyK/ggcpqrPQ3DzEFg6UHtN8MkUoA8OUP6QMvRBHGXpgzRy64NByL8M2Kyqr1KA8+9LmZXDSuD9IrKJwMzb58rXEJzwR4CvA78H+vBcp7qJ1Cs/wMWqeiywyH0uaabAIjIa+CnwYVV9KalqjTJNKG8KGcgP5emD2CZqlBWxD5LIrQ/qlV9E3gR8GXhfWFSjWiFDRkurHFT1L6p6uqr2Aj8k8Aujqn2q+hFVnaOqS4EOYDvBA7cr0kTNdaqbxSDkR1Wfc393Az+gia4OEWknuCm+r6p3ueIXnJkcuit2uvK4NcFz64OM5C9TH8RRlj6IJa8+qFd+EekC1gHLVfVpV1yo51ASpVUOYYSCiLQAnwZudf+PFJFRbnsx0Keqf3Ym324RebMzQ5cDP8tH+vrld26mCa68HTiHwDXVDFkF+A6wVVVviny1HljhtlcwcD7XAxc5H2s3MAN4OK8+yEr+kvVBTUrUB3Ht5NIH9covIh3AL4BPqervwspFew4lkteb8Ho+BCPr54H/EWje9wAfIogY+CtwAwMT+o4gyNa6Ffg1MC3SzjyCC+lp4JbwN2WQnyB6YxPwOMELrlW4CJomyL+QwPR9HHjUfc4CxhO8PN/u/o6L/OZad563EYnGyKMPspK/hH3wN4JAiD3uuuspWR9UyZ9XH9QrP8GAb2+k7qPAxLzO/2A+NkPaMAzDqKK0biXDMAyjcZhyMAzDMKow5WAYhmFUYcrBMAzDqMKUg2EYhlGFKQfDaAAicoWILK+j/hESydprGHnTlrcAhnGwISJtqnpr3nIYxoFgysEwauCSq/2KILna8QSTFZcDxwA3AaOBXcClqvq8iNxHkAdrAbBeRMYAe1T1ayIyh2AG/EiCiU8rVfU/ItJLkEvrZeDB5h2dYaRjbiXDiGcmsFpVZwMvEay18Q3gHRrkxFoDXB+p36Gqb1XVGyvauR34pGvnCeBzrnwtcJWqzm/kQRjGYDDLwTDi+YcO5MX5HnANwcItG1yW6FaCtCghP6psQETGEiiN+13RbcCdNcq/CyzJ/hAMY3CYcjCMeCpzy+wGtiSM9PfW0bbUaN8wCoO5lQwjnqkiEiqCdwF/ADrDMhFpd/n6Y1HVF4H/iMgiV3QJcL+q/hd4UUQWuvKLsxffMAaPWQ6GEc9WYIWIfIsg6+Y3gHuBm51bqI1gQaYtKe2sAG4VkZHAM8BlrvwyYI2IvOzaNYzCYFlZDaMGLlrp56o6K2dRDCMXzK1kGIZhVGGWg2EYhlGFWQ6GYRhGFaYcDMMwjCpMORiGYRhVmHIwDMMwqjDlYBiGYVRhysEwDMOo4v8qoarANEH8rQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ " sorted_data['inc'].plot()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXd8JHd99z+/7V1lV+0k3UnXfMX2ne/OxrgQcME2CTHNxCQQJw/ghIeQSgI8T3hSgCdACAQIONQHAwnNxBhiXLAxxvXsO5fz9aYraiettNL2mZ2Z3/PHzG92tmp3tX1/79dLL0mzM6vfjnbnM99OKKXgcDgcDseIqdEL4HA4HE7zwcWBw+FwODlwceBwOBxODlwcOBwOh5MDFwcOh8Ph5MDFgcPhcDg5cHHgcDgcTg5cHDgcDoeTAxcHDofD4eRgafQCKiUQCNCxsbFGL4PD4XBaiv379wcppX0r7dey4jA2NoZ9+/Y1ehkcDofTUhBCzpayH3crcTgcDicHLg4cDofDyYGLA4fD4XBy4OLA4XA4nBy4OHA4HA4nBy4OHA6Hw8mBiwOHw+FwcuDiwOFwOE0MpRT37J9EQpTr+ne5OHA4HE4TczoYwwd/9DIePjxb17/LxYHD4XCamJggAQDCSamuf5eLA4fD4TQxyZQCAIgLXBw4HA6Ho5FMqbGGGBcHDofD4TASTBx4QJrD4XA4DG45cDgcDicHQYs5cMuBw+FwODpJiVsOHA6Hw8mCFb9xceBwOByOTlJ3K3Fx4HA4HI5G2q3EYw4cDofD0eDZShwOh8PJgYsDh8PhcHLQ22ekZCgKrdvf5eLA4XA4TQyzHChNV0vXAy4OHA6H08QYBaGeGUtcHDgcDqeJSRrFoY4ZS1wc2oQXzoVwIZxs9DI4HE6VYTEHoL5BaS4ObcKd396Hrzx+utHL4HA4VSaZkuG2mQGo4rAQFUBp7QPTXBzahHBSwnIi1ehlcDicKpNMyfB77ACAyVACuz/+CL719Jma/10uDm2AolCIkoJEqr550BwOp/YkUwp63TYAwLELEQDAcLez5n+Xi0MbIEiqTzJR55a+HA6n9iQlGQGPKg5HZsIAgNFeV83/LheHNkDQeq/EuThwOG1HMiXD71bdSsdmVcuBiwOnJFg2Qz0LZDgcTu2hlKpuJc1ymIsI6HFZ4bFbav63uTi0ASwPmruVOJz2grmMvQ4LbBb1cl0PqwHg4tAWJLlbicNpS9iNn8Ni1tNZR3u4OHBKROBuJQ6nLWGfaYfVDLfmSmoay4EQMkoIeYwQcoQQcogQ8mfa9l5CyC8IISe07z2GYz5CCDlJCDlGCLnJsH03IeQV7bEvEEKItt1OCPmBtn0vIWSs+i+1feFuJQ6nffjar0/j60+oBa0snui0meC2MXGofRorUJrlIAH4K0rpVgBXAng/IWQbgA8DeJRSugnAo9rv0B67HcB2ADcD+DIhxKw9110A7gSwSfu6Wdv+bgAhSulGAJ8D8KkqvLaOISmlLYd6tvTlcDjVRVYovvSrk7j3xSkAWW4le5O5lSilM5TSF7SfIwCOABgGcCuAu7Xd7gbwJu3nWwF8n1IqUEonAJwEcAUhZAiAj1L6DFVrv7+ddQx7rnsAXM+sCs7KCAZ3Eos/cDic1uPg1DKW4iksxdVuB8lmdisZ0dw9lwHYC2CAUjoDqAICoF/bbRjAecNhk9q2Ye3n7O0Zx1BKJQDLAPzlrK2TYZYDwIPSHE4r8+vj8wCApbgIIB1zsFtVtxIhwJpuR13WUrI4EEI8AH4M4M8ppeFiu+bZRotsL3ZM9hruJITsI4Tsm5+fX2nJHYOxpS+PO3A4rcsTJ4IAgJgoQ5QUPdnEaTVjsMuB8YAbdou52FNUjZLEgRBihSoM/0Ep/S9t8wXNVQTt+5y2fRLAqOHwEQDT2vaRPNszjiGEWAB0AVjMXgel9KuU0j2U0j19fX2lLL0jEAyWA89Y4nBak0gyhRfOheDX+igtxcUMt9JfvX4zvvfeK+u2nlKylQiAbwA4Qin9rOGhnwK4Q/v5DgD3GbbfrmUgjUMNPD+nuZ4ihJArtef8/axj2HO9DcAvaT160rYJxpgDdytxOK3JS+eXICkUr98+CAAIxVN6DNFhNcPrsGLAVx+XElCa5XA1gHcBuI4Q8pL29QYAnwRwIyHkBIAbtd9BKT0E4IcADgN4EMD7KaXsivU+AF+HGqQ+BeABbfs3APgJIScB/CW0zCdOaXC3EofT+swsq8O6Lh3pAqBaDgkx7VaqNys26KCUPon8MQEAuL7AMZ8A8Ik82/cBuDjP9iSA21ZaCyc/xklRvG03h9OazEcEAMDmAQ8AzXLQ3Ur1r1fmFdJtgCCV71b68q9O4u9/eqhWS+JwOGUyHxHgdVgw2KUWuS3FxQy3Ur2pfWs/Ts3JsBxKFIdfH5/HYkys1ZI4HE6ZzEWS6Pfa0eOyAtAsB+3zbLfU/z6ei0MbkEzJsJlNEGWl5GylpXgqI8uJw+E0lrmwgH6vA06rGTazCUsJEZSqLqVG1ARzt1IbIEgKetzq3UapbqXFmKjnUDMePXIBUYHHLDicRjAXEdDvs4MQgm6XFUsxNebQCJcSwMWhLUimZHQ71dzoUsSBUqpZDul9Xz6/hHffvQ8/fWm6yJEcDqcWUEp1txIA9LhsCGl1Do46Fb1lw91KbUBSUuCwmeGwmjLSWgsRF2WIsgKLlDZV739lBoDq9+RwOPUlIkhIphT0e9U6hm6XFUvxFOxWM5w2bjlwKiSZkmG3mOCyWRAXV3YLhbS+LSzmQCnF/QdUcQjxIDWHU3fmwmoaa59mOXS7rFhKiFiKi3o31nrDxaENECQFDqsZTqu5JLdSKKZ2fJQVCklW8PLkMqaWEgCABS4OHE7dYRa70a20EBXx0rklXDrS3ZA1cbdSGyCkZDi8djht5pLcSsxyAABRVvDgwVlYzQTr/O6MxzgcTn1gBXD9PmY52PQbtVeN9zZkTdxyaAOSKRl2qxkuW4mWg0EAhJSCyVAcoz0ujPndWIhyceBw6k3araTGHFitAwBcwcWBUymCpMBhMZXsVmKDRNixyZTqlvK7bdxy4HAawHxUgN1igs+hOnN6XGr24dpeF4a66jMWNBsuDm0Ay4V2VeBWEiRZO96EXo8NizERvCEuh1Nf5sJJvcYBALo0y6FRVgPAxaEtUO/8TXCW6laKGcVBQTIlw2kzo9dlQ0qmiPBCOA6nrsxFBD2NFYA+04GLA6diKKUQJBl2ixlOq6Wk3koho1sppSApqYU2vdobkqezcjj15fR8DMPdaffRZWt78PE3XYzf3rGmYWvi2UotTkqmULT+K2pAuvBd/2NH56BQmuNWSoiqW4qJw0JMxDq/u+Zr53A4wNRSArPhJHatTaesmk0E77xyXQNXxcWh5TG29HXazEUb733m4WMIJ1PodtrgsVsQFSSIhoA0txw4nPqz/2wIALB7XeNcSPng4tDisAC0XSuCS6YUKAqFyZTbxfHcYhyRpISgVcRQtwPReUmPOTispgzLgcPh1IcXzobgtJqxdcjb6KVkwGMOTUwwKuD0fDTvY5RSvPGLT+I/954DAK19hlpmn896WI6nEElK+uNDXWrwi2UrOQ2WA5/zwOHUj/1nQ9g52g2Lubkux821Gk4G//Lwcbz32/vyPiZICl6ZWsYvDl8AAD2VFcjfmfV8KJ7xOxtUnkwpauM+7Xi7xcTdShxOnYiLEg7PhLF7XU+jl5IDF4cmJpxIFXTxhJNqxtHxCxEAgMNi0pt2zS7ndlY9v6iKw3hADTQPauIQFSTICtUHivS6bdytxOHUiZfPL0NWKHaPcXHoSOYjAr7x5AQUpbziMkFSEBOkvEVpUc1FlJLVx+xWM8a0C/+ZhVjO/uc0cbhtzwgAYFBzKy0nVJFhA0V63TZuOXA4deJ0UHUbbxlsrngDwMWhLjx4aBYf++/DeOFcqKzjUrKClEzzjvNk8QOGw2LCul5NHIK54nA+FEeX04q3XDaCnaPduHxMzYzIJw7ccuBw6sP0UgJmE8kogGsWeLZSHWBDwh87Noc9Y6Wnq4maKMQEKWdUYI44aKmsgz4HJgyWw9d+fRpehwXnFxMY7XVisMuBn7z/an0K3HI8VxyYlcHhcGrL9FISgz4HzHmyCxsNtxzqALsQP3Z0vqzjRFkVh3xznSPJVMbv7OI+FnBlWA7fe/4cPvngUZyaj2K0x6Vvt2mZEcxycGrH97i4W4nDqRfTS4mMyuhmgotDHUim1Iv84Zlw3mBxIVLFxEHbpvXpgt2i/ivHA26cXUjf+YcTEpbiKUyGEljbmxYHQgjsFpPBraQeb7ea8rqxOBxO9ZleTmCou/lcSgAXh7pg7JT6+PG5ko9Lu5VyU1OZW2lzvxrIYpbDOr8bCzER4WQKlFKEE2kLY8QgDgAyxIFZDnaLGYKk8M6sHE6NURSK2eUk1nDLoXNJSjJ6XFYMdzvx3WfPQS4xa4mJQ1RI5TzG3Eo7RrsApO/8x/zpoLQgKRBlBd1a+9/Rnsw3od1q1sXBrokDex5uPXA4tSUYFZCSKReHTkbQehd96JYteGVqGd96+kxJx6VjDvktB5fNjKs3BjDgs8NtV3MLWB3DRDCmWw3vuWYc77lmHK8a92c8Rz63ksNi1tfM4XBqB5vbvqaLu5U6FlaB/MZLh/C6i/rwmYeOYamEiWvGbKVsokkJHrsFv71jDZ79yPWwagFmFlc4uxDXL/xr/W787W9tg9OWmfFkt5j0eIbuVtJEgjX043A4tWF6SY0/csuhg0mmZNgtagXy71y+FomUjMlQYsXj9IB0Ml9AOgWvwwJCiD49CoCeznp2Ia5XUXc5rTnHA4DNkhYLFrPglgOHUx9mljXLgYtD5yJIiu7T92jun2KttRnpmEO+VFYJXkf+i/6Az475qIBwQj2OzaXNhmU4AWnLgYkEtxw4nNoytZSA22Yu+PlsNFwc6kAyJcOhXYhddvXim89VlE2xOodwUoK3wJvK77FjISrobiVfAcvBKA4OPVvJpK+Zw+HUjumlBNZ0OzMs/2aCi0MdEFKyfvEt1jnVCKVU75uUT0giyRR8BSwHv9uGhai4olvJbqi6ZqLA1smzlTic2jLTxGmsQAniQAj5JiFkjhBy0LDt7wkhU4SQl7SvNxge+wgh5CQh5Bgh5CbD9t2EkFe0x75ANLkkhNgJIT/Qtu8lhIxV9yU2HkFS9Iuv26be7a9kOTCrAchvObCAdD4CXjsWYoLeGqOQhcHWZLeY9OFALGuJWw4cTu0QJQXHL0T07MJmpBTL4VsAbs6z/XOU0p3a188BgBCyDcDtALZrx3yZEMJuT+8CcCeATdoXe853AwhRSjcC+ByAT1X4WpqWZB7LYaWYgygVF4dIMbeS24aUTHE+FIfDaoLdYs67X7a1oG7TYg48IM3h1IyD08tIphS8ary5RoMaWVEcKKW/BrBY4vPdCuD7lFKBUjoB4CSAKwghQwB8lNJnqFp6+20AbzIcc7f28z0ArifN6oSrEHVGs2Y52JnlUFwcmEtJ3VfKekxBIiUXDEizuQ4TwVhB1xOQFgK2NuPPAg9Iczg147kJ9ZJ6eSuLQxH+hBByQHM7sUkVwwDOG/aZ1LYNaz9nb884hlIqAVgGkFmt1eIkJVm/EKspreoEqGIYLYfsDqwstbWw5aCKw+n5WMF4AwDYNMvByS0HDqeuPDexiA19bgQ89kYvpSCVisNdADYA2AlgBsC/aNvz3fHTItuLHZMDIeROQsg+Qsi++fnyOpw2EsFgORBC4LZZSrAc0hfnWJaQMDeTp2C2kjoLeiEmFsxUAvK7lXjMgcOpLbJC8fyZRVwx3tz3wBWJA6X0AqVUppQqAL4G4ArtoUkAo4ZdRwBMa9tH8mzPOIYQYgHQhQJuLErpVymleyile/r6+ipZet2hlCIpyRkXYJfNjESquOXAsoW6nNYcIWFZSIXyo413I8VyqFk1dEbMgWcrcTg15ehsGJGk1NTxBqBCcdBiCIw3A2CZTD8FcLuWgTQONfD8HKV0BkCEEHKlFk/4fQD3GY65Q/v5bQB+SduoJagoK6AUOeKwkuXA3Eo9LmtOhXREdyvltwp6XFa9lXcxt1K+mAOvc+BwasvBqWUAwK61zTc32siKpXmEkO8BeC2AACFkEsDfAXgtIWQnVPfPGQB/BACU0kOEkB8COAxAAvB+Sim7yrwPauaTE8AD2hcAfAPAdwghJ6FaDLdX44U1C+wO3Fhw5rJZVow5MLdSj9uGMwtxiJKixwgiK8QcLGYTel22itxKLCYicHHgcGrCkpZizty/zcqK4kApfUeezd8osv8nAHwiz/Z9AC7Osz0J4LaV1tGqsDtwY8GZ216C5aCJQ69LfQPFBAk2i/oza+FdyHIA1DfeQkxcIVspNyDNhgBxtxKHUxvCyRTMJqKntTcrvEK6xrAGdo5sy2GFO/OUlLYcgMxaB2Y5FCqCA9IZS0XdSln9lPTtFjN3K3E4NYLVKDV7xj4XhxqTz3Jw2cyIr1AhLTDLoYg4FHIrAWqVNAD4nEUC0mbmVsp8GzisJp7KyuEA2H82hCMz4ao+ZzhRuPVNM8HFocYw90yO5bBCb6V0QFoVh39//BT++Dv7AQALURFOqznnjt+IXxOVom6lPNlK7HdeBMfpdARJxnvufh4f/cnBgvsUih3uPb2gZxVmU6y7QTPBxaHGMMvBkRVzKDUg3etWL+73vTSNhw7PQpBkTC3FMdxTvGEXq5Iunq2UXxzsFm45cDgPHbqAUDyFA1PLeW+W5iMCdv7DL/DMqYWM7QlRxu9+fS/+c++5vM8bLtI0s5ng4lBj2EXWeAF22syIlWk5AAClwFQogamlBIZX6OaoWw4lpLI681gOfJ4Dp9P53t5zIET9LB6aVl1LCVHGPz90FJFkCnORJERZwfnFeMZxUUGCrFDMLOUf6BVOSEXdvc0CF4caw+44jKmsbpsFoqRkVEFnw8SBWQBbBr0AgMlQAlOhxIqWw7Wb+3DrzjXY2O8puE/acjDlbG+XSXBxUcLscrLRy+C0GBPBGJ45vYA7Xj0GAHjhbAgA8LMD0/jSY6ew9/Si7jLO9gIktBu/uYiQ97kjyVTRTMNmgYtDjclnORSb6ZAQZcQESReOkR4Xvvqu3fjiOy4DABy/EEEonsLICuIw3O3E52+/rGhcoljMoV0shy89dhJvvevpRi+j5fjF4Qu461enGr2MhsHE4F2vXoeRHideOKf+/sArMwDUfmnsBiqRdSMV17ofzBcQh3BSagm3UvPbNi1OOuZgsBzYqFBRzokJfPCelxGKibhuSz8AtTne67cPQlYorGaCZ0+r/s2V3EqlkK6Qzk1lXYiKq37+ZmBmKYn5aP4PKSc/33xyAh+7/zCsZhP++DfWN33KZS1YiKnvmQGfA7vX9eBZLcD85MkgAPWmj91AJbIsB3bTl+99JysUUYEHpDlIz2I2zlRglkN2Qz1JVvD4sXlMLSX0lt02Ld3UbCIY7nZir9bqdyXLoRSYMBnjGoBqUbSL5RARJIiSAllpm44sNWV2OYl//O/D8NhV1+dKWXXtykJUhM1igttmxq61PbgQFvBPPz+qfy6TKaPlkHmOmFvJaDnMLCfw3WfP6q1wisUCmwUuDjVGL4KzZqayAkA8q0r6lallRAUJ0aSkxxxshljFaK9Lr3EY6XGtem2jvS785P1X61YKw2Ext03MIaKlE/KivtI4fiECALhx2wAAYDHWHhZkuQSjIgJuGwghuGn7IIa7nfjec+f0RpaCpOjxxGwBZb/HNRcxANz74hT+9icHcSoYBVC8RqlZ4OJQY9gdeEYqawHL4WktJS4iSBBlGWYTgdmUNumZtWAzm9BXpT7wO0e7M/6GulZT29Q5sOLBlSbvcVROz6sXrz3r1I6hoXhnisNiTIBf+4wNdjnwyw/+Bj5z2w58/nY19lfMcjAGqFlQejmh3qQcm1XFl8ccOHpAmrmHADWVFUibnwwWTxAlBTFBhtWcedFm1sKaboc+87kWqO0z2sVy0MShQ90j5TIRjMFjt+AiLTtuoUMth4WYmNEYz24x4227R8AaRgsp2RBzyO9WAlTX0njAjXCOODT/pZdbDjVGkGTYLKaMi7k+KtRwhyFIMp4/s6jXHCzGxAxBAVQ3EIAV01hXi9o+oz0upszH2y6vp9acDsawvs+t18mEOlUcoqLeusaIsTFlYcshUxwAtbYBMIgDjzlwhJSS0ToDMKSyGmIOB6fUgePXbgoA0MTBkplFNKqJQjUylYrhsJohKRRSkTqMViHC3UplcXo+hvUBt97wsRNjDpRSBKNCwRGeDqvamJLdcGTHHIzvtfmIWmPDWmkc02I6PObQwbx0fgl/+YOXEBOknFRRty3Xcjg4pVZgXr3RIA5ZbqW1muUwWoVgdDFYcVyrt+0WJFkP7HO30sokRBlTSwmMBzzwOSywmEhHikNclCFIim49ZcMaU7LPR7ZVGhclmAhgMRE9nZXFHNj55DGHDuaJ4/P4rxensGO0Wy82YzjzFMEdng6j123TK5oXY2JO5bLfY8ddv7cLV9R4vCATs2RK1l1grYhxgh63HFZmIhgDAKzvc4MQgh63rSMD0qzGx1/AcrBb1CJRJgrZNx5xUYbLZoHHbsFcmLmVMpvwFZr/3kxwy6FGsIvRkZkwHJbcxnZmE8nIajg8E8a2IZ9ubi7GRFjNuf+eWy4ZKvimrRbtYjlEDOLAYw4rYxQHQB001YmWQ1ArgCtmOQgGyyHHrSTKcNrM6PPadcshbHgvumzmvJ/tZqP5V9iiMHEQJSXHrUQIyZgjnZIVHJuNYNsanz7AR5SVjBqHemK0HFoZ4wwMbjmsDEtjHQ+o4tDjtiIUy992up1JWw6FxEG1HFi6d65bSYaLiUNEAKUUy4mUPte9FVxKABeHmmF8w9jzXORdtnTb7lPzUYiygu1rfBnmZuPEQf27rZ7OauynnxBb+7XUmmBUwL0vTWGd36UXafa6bXobiU5ikVkOhQLSFhaQzm85xEUZTqsZfR5VHOKiDFmhGPerotsKwWiAi0PNMPoh8zW/cxsG/hzSgtHbhnzw2tN3FY0yPVmrj1YvhOMxhzSPHZ3Dlx47iadPBvVcfYYgyXjn1/dieimBT731Un17r9uGULy2lsPTJ4P410eO1/RvlEuQWQ4F3Ep2qymjQjqRkjPOaSIl6ZZDMCrorjlWO9IKaawAF4eaYbwYZQeWAcBlN+ul9YdnwrBbTBgPuOGwmvSK5XwWRz2wt4nlYHQrtbqLbLV87L8P458fOobf/fpevDy5nPHYybkojs5G8Hdv3I4r1/v17b0uG5biYk37Uv3swDS+/KtTOYLVSBaiIjx2S8GOxqxI1Pj5MP7MAtKDXQ4oFDipueu2DPoAcMuh4zGamvY8bzKP3aLHHI7NRnDRoBcWswmEED3ukF0EVy/0mEOLWw7GgHQnp7IqCsXkUgKvvagPAHBGCzwz2HliqdKMHrcNCk2nYdaCcELtI9ZMNyILMSFvARxDDUjLGZa18WYwocUc1nQ7AADHtcI33XLgMYfOJpmS4dUu8vksAI/dqhdoheJiRq8kJg6Ncytp2UpN9IGtBGY52MymjnYrBaMCREnB1RvUGprJUObkMiYOnqy05d46FMKxuFAtBahcFqJiwWA0YCyCM1oLkuFnVRyGutRiVVYVvabbgT6vHf3e2mYbVgsuDjUikZKxfVg1I/OZpx6DWykqSBmBaGZ2NjpbqdVjDuFkCjaLCT6npaPF4XxIHVe5sd+DgMeOyVDm+ErWuTbb3cHEIV+twyuTywWH2ZQDy/9vFnEQJQVTSwn43YUv4Hr7DMnoVkq/v+KiDKfNgjWaOBwxNNv7wZ1X4gPXbarR6qsLF4cakRBl+D12/I+rx3F9VktsQC2CYXe20aSUcdfWaMtBF4dWtxySErya7zjZwW4lZimM9Dgx0uPMEQf2PsweXcnmfOQb/PQH/+85fPGXJ1a9Npb/3wzisBAV8Ja7nsJEMIbf2BwouB+zHISUrPdCM7qRE6IakPY5LXDZzDg1p8YcupxWrO/zoMvVGm6l1oiMtCDJlAKn1Yz/88ZteR/32K16Nk0kKWV8MD0NthyYW6kdYg5ehwXWDncrMTEY1sTh4FRmQJq5lUq1HARJxkJMxFSWyFQCs1qWmqAS+8FDszg4FcYX33EZ3rhjTcH9HFYTkprl0OOyIrEs6zEtSiniKdWtRAjBYJcDp+fVGE+rBKIZ3HKoEXFR0u8q8uGxmyHKCiLJFERZyXjjeIrEKuoBawy4VOM0xlrD3HUum7njxcHvtsFls2Ckx4WppQQUQwZSOJmCzWzKcX/6PTY4rCY8fmw+YzuLQVzQmsp966kJvPPre/F7X39Wd5WWAqVU71baDJYDW8sNWweK7me3mPVxn12adcXeX4KkgNJ0ixzmWnLbzLC0QFW0kdZabQuRSMn6GyQfTAAuhJMZvwPpO4zseQ71wmWzYMugV5+X26owd53Dau7obKXJUFwfFDXS40RKpvoQGiBtYWVjt5jxges24cFDs3js6Jy+nbmZWN+gzz96As+dWcRTJxdwdiGe8zyFECQFotb5txnEISqkYDGRvKnnRtjj4WQKPZqLiL2/mHvJpQntUJeasZQ9K74V4OJQAxSFIpnKbZthxKO5kWaWVXHIZzk0yq0EAK/fNoB9ZxZbup9/OJmC12GF02bu6DqHqVBCHxTFZoEYM5YKiQMAvPfa9djY78Hf/+yQXovABgAFowIiyRRC8RR2jHQBgH6xLwVjM7rsxnSNIJJULU1Cit+Usc81pem4DLMcWNYSqzIf0trrt0rhmxEuDjWAZTGs5FYC0uKQGZBW30g2c+Hja80N2wagUOCxY3Mr79ykRAU1IO20dq5bidU4MMuBzQSZWkrHCyKaiObDZjHhtt0jOLsQ1wPXC1ozOYUCB7SCOtaPSSjjPBvbmzSF5ZCVGFIIYyPNbs1yiCQl/O97X8GhabXbQdqtpFoOrVLbYISLQw1gFyJXUbeS+maZZeJgtByYW8nSGLcSAFy8pgsDPjseOXKhYWtYLexOsJPFgdU4jOiDolQLwpixFFnhotjvU9PSdXeWAAAgAElEQVQ6Weqqse7hhbMhAMD6PrXVfDmWw3IiHZ9YagJxCGclhhTC2IKfWQ4Hp5bxH3vP4fvPnQOQ/uwPMnHglgMHSItDUctBEwBmORjvLLwNrpAGAJOJ4LotA3j82HxTtTYoFUrVgKHXYYHDZu7YxnusxoG5lZw2MwIeW4ZbKVrErQQAfR71AsfEIWhIbX3x/BKAtOUgltHmnWUqEdIkloOQ0j97xbDnsRyOarUMrDWJbjnobqXWylQCuDjUhITmd3QUtRzUx2aXE9rvuZZDI2MOALB5wIOYKNe8+VotSKYUyAqFx26F09q5MYcZ7f01pLVyANQxs5mWQ2G3EgD0aRW9TBQWogIsWv+vF8+FYDUT3TIpRxxYjcOA19EU4lAs9mLEGLD22C2wmgmOa+M/mVWlxxx4QJpjhN2lFo85ZAakPfkC0g1OfRvwqW9sllHVSoQNVb/MrdSKFtBqYTUMRsvU77FnFLatdFFk4sDmIS/GRGzQ3EiheApDXU49SFtJQHq019kU4pDdqaAQRsvBbjXBaTXnDMZibiWvw4o3XzaM12zuq+5i68CKVx9CyDcJIXOEkIOGbb2EkF8QQk5o33sMj32EEHKSEHKMEHKTYftuQsgr2mNfIFpKACHETgj5gbZ9LyFkrLovsf6U41aaDefJVmoSy2FA8zW3ojhMawHXoS4HnDY1Lz2RkvVgaqfA6g6M4167nFZdPBWFIipK8BW5KHY7rRnzkIMxEQNdDgS0/kPD3c6KpgeyNYz2urDcBNbpSu41htFycFjMeVPWjZ/9z/3OTrzuotwuCc1OKVefbwG4OWvbhwE8SindBOBR7XcQQrYBuB3Adu2YLxNC2Fm6C8CdADZpX+w53w0gRCndCOBzAD5V6YtpFnRxsBU+vSwPeimuFiAZ70YGfQ6YTQT9Xkehw+sC+/tzVeihU2+MVcHsrvaLvzyJGz77eEYBWLvDMow8WeLA7tSjogRKc1tnGDGZCALa4BpAdSv53Tb0ae+P4R6nfiNTljgkJNjMJgz4VLdSoy07NTC/svvHmKJut5p0F9IWresqUDwZpVVYURwopb8GsJi1+VYAd2s/3w3gTYbt36eUCpTSCQAnAVxBCBkC4KOUPkPVd8C3s45hz3UPgOvJSonGTQ4riHFaC9+FmEzp1tzZpuyabiee+tB1uHqjP9+hdYNlqcy1oOWgi0O3U7+Le25iEaF4qilcGPUiJqiV+mxGCKBmzkSSEmSFFmydkQ0beQmobiW/26ZblsPdTti1tOtyA9I+pwVdTiskheZMVKsngiTndCooRIY4WMz676+9qF8vXGWC0cpU6rcYoJTOAID2ndlMwwDOG/ab1LYNaz9nb884hlIqAVgGkPeqSAi5kxCyjxCyb35+Pt8uTUEipX7gilVIA+m7uXxvyMEux4rFOLXGbjGjx2XFhXDzWA6ipEAqwa89tRRHl9OqFcGpb/MjM2oOeieNvowKUoZLCUgHRyPJlKEja/E75oDHhvmogLgoIa41lRzQLIcRg+VQbkDa57Dq62mkaJcqkkBmWxuH1aRbCRv63Fgf8ICQ/AO+Wo1qv4J8VzNaZHuxY3I3UvpVSukeSumevr7mDfCUEpAGALeWsVRK4U2j6Pc6mirm8N5v78NH7zu04n6ToXThV3bnzGCeLqPtSlSQ9cw4hvFiXK7lsGAYocksy+FKxSGRgtdpRXcecUjJCk5pE9TqQbTATIt8ZFsO7P21zu/G5kEvnFZzw2/sqkGl4nBBcxVB+87KaCcBjBr2GwEwrW0fybM94xhCiAVAF3LdWC1FKQFpIN1Co6nFwWfHhSaKORydDedMMsvHVCiBYS3HPLuNSbCDgtKxIpbDciKVviiWIA7BqKifO7/HpovvOr8bZhOBxUQgyuVVSPscFn09xkaP33hyArd8/omMITq1JC2SpcQcMi0H5iEY87vwh1eP4a9vuqg2i6wzlYrDTwHcof18B4D7DNtv1zKQxqEGnp/TXE8RQsiVWjzh97OOYc/1NgC/pI2OTK0SllPvKBKQBtLFbqW8IRvFgM/RNDEHWaEIRsWM2dD5oJRqloNW+JUlDvnmE7QrUSG3+tkoDixjqFi2EgD0eeyQFYqT2mwCv8eOW3cO4wd3XqmLsM1iKmsGSDiRgs9h1auHjZbDr4/PQ5SUutXYRAT175Ryo5aRyqpZDk6rGX1eO3at7cEfXj1es3XWkxXPBCHkewBeCyBACJkE8HcAPgngh4SQdwM4B+A2AKCUHiKE/BDAYQASgPdTStmtxPugZj45ATygfQHANwB8hxByEqrFcHtVXlkDSYgyzCayYp0Ccys1c5/3AZ8dcxEBikJhMjXWVA5pw+5XEofFmIhESk67lWzZ4tBZlsOgLzPrLb9bqfgNCstMYpXAfrcNDqsZr1qfDg/aLKby6hySkh6QBtJ1D8mUjH1aWw41JuIs+TmLoSgUf/Wjl/HOK9di97rejMeiZcQczCYCq5kgJVPYrSbccvEg1vldbeFKMrLimaCUvqPAQ9cX2P8TAD6RZ/s+ABfn2Z6EJi7tQkKbELXSm4WlzTW3ODggKxQLMVEvhmoUrEU0u6Dl45XJZUiKeoEazoo5eB0W2MwmzHeQ5bCSW6mcmAMAPHrkAqxmkve9YDObys9Wclj1FhRLCfX/8uK5Jf15wonquZWiooR7X5yC323LEYdyAtKAWt+QkiXYLWbccskQbrlkqGrrbBaa96rUwsRFuWi7bgZ7IzZ1zMGbrpJutDiwIqyokN/VcGQmjDf+25O4WJvdzSwH9r9Y3+dBUuysQrhi2Uos5mA2kRXjY+x/f2Yhjv9x9Xje97fNUro4CJKMZEqBz2mFx66KNmsF/szpBX2/arbyZi6vs4u5MycKjUothN1qRkSQ2iIrqRDt+8oaSDIlFy2AY+jZSk1tOWi1DpHGxx1Ynn0ypSCVx33xzCn1onJwSk1ZHelON5sDgA0BN/wem34R6gTUmEPmhdxhNcFqJrrl4C1hhgETB6/Dgg9ctzHvPnaLCUKJbiXjnTohRP2/aBbds6cW9BGlxrbeq4XFAs8u5CY0sJRet7204jWWztroFje1pH1fWQNJiPKKd2KA0a3U3AFpAE1R62AUqHzjKJ8/s4iAxw6n1Qyv3aJ3wmR3pluGvAh47B2TrSTJCpIpJafqlxCittBIpLSmeyvfnLhtZlw+1oOP3LIVPdqFOxubxVyy5RAXWFt79W8HPHYsRNXY1kvnl/R2E1W1HLS1nV2I61XyT50M4i9+8JLaqcCS2amgGA6rCXaLqe3iDEaa95a1hWExh5VgFkMpbYIbRZ/XDkLScycayXzWaMtuV/oiRSnF82cW8ZpNfbhmUwBngjH9g+uwmvGT91+N9X1ufOrBox2TrRTT6jry3Q37tBYaUUFGtzP/xd4IIQQ/+uOriu5js5hKbp8RT7GJaera/B6bmiobEyDKCrat8eHHLxSPL5VL0jDneS4iIJmS8b7v7kc4KeHSka4VM7aMOKzmklzHrUzzXpWamINTy5haSuCm7YMZ2xOijGBUWHF+NMPTAkVwVrMJQz4Hzufx09YbY4+n7IyliWAMwaiIy8d78ZZdI9mHYtsaNQ4R8NgRFSQkU6XFhVqZWJ6+SgzVcpAwEYxh97qenMcrwW42QZRKq3NgBYnsc+J323F8NqInHYz0qG1PqulWEgxrmwjG8LH/PqxX2x6YXMaY31Xyczms5oxK6XakvV9djfjaE6fxwR++nNMo7GtPnMbrP/drzEeEstxKzRxzANQip4k8ftp6Mx8RwKz4bHF4/oxaN3n5WPELnV9ziXSCaymapyMro8tpxcxyAlNLCWzs91Tl75UTkE5q4sAaUAa8NgRjom6hDvgc8DktVc1WMtZgPHBwBodnwvhfb9iKtb2qKJTj3rVbTG1/c8HFoQJigoSIIOH8YiJj+0QwhkRKxkQwVpLlcOX6Xrz7mnHsHO2u1VKrwljAjbMLjbcc5iOCnoEUzXI37D29iF63TZ8zUIiARw2sdoJrKV9HVkaX04pT86rgb6qSONjLqHPIthwCbjtEScFJrWXGoM8Bn8Na3YC0wXL40b5JEAK8ftsArhhX01rLseC55cDJS0wLph2eWc7YzqZuAbktG/LhdVjx0d/a1vR3IGN+FxZjYsO7mc5HBIwH1AtZxGA5CJKMR45cwG9s7lsxQOjXZhB0QvM93a2UxzI1TiZrhOUQz5qzHvCq/5dD02GYiNroz+uwVDXmwCwHQtS44O61PfB77Lo4lFNv9HuvWos7X7O+amtrRrg4VAB7Yx+eDmdsnzEEbUtxK7UKY9p84FJ6GtWKuCghKkhYr63FaDk8fmwe4aSEW3euWfF5mOUQjLS/5aAP+snTPpqJg9lEsM7vrsrfKycgzUbpOrW1+d3q/+XQ9DICHjssZhN8ztpYDus0N9KN2wYAAFeMaZZDGeJw/dYB3LZndOUdWxguDhXA3tiHDOJAKcXMclI30dtKHLSLx5kGxh1YptL6Pk0cDIVw9708Db/bhqs3BlZ8HmY5zHdEzEG9GBZyKwGqVVitiYPlVEgnsmIO7P8yEYxhUJu77HNYa1IEt3lAHcpzgyYO6/wubBvy4aIBb8FjO5HmjoQ2KWm3UlocQvEUREnBW3eP4FtPncG6QHXuxpqBdVoWx5lg4+IOTBzW+d0gJG05RAUJjxy+gN+5fBTWEgqSXDYLAh4bTs83PsBea6JFCrtYs7tquZSAytxKLObQp1l0lKar8n1OC8I1SGW946ox7Bjt1uNThBDc/6fXtHXNQiVwcaiAREqGiahupMWYiF63TZ9ZPOZ34YkPva6kC1Wr4LCaMdTlyFtZWi9YGmu/1w6P3aLHHA5OLUOQFFy/daDk57psbQ9eOBeqyTqbiXSdQ2HLoVHikBBlEJKuNDYW1g12sWps1XKglFblws1cXpet7c6xMrkw5NI+V7A6EhMkbB1S8+ZZ3IGl4A11OdtKGBhjDU5nZf2QAh47vHaLbjkwi2Koq/R527vX9WAiGGv7HktRQYLVTPJm1XTXQBzsFnPJ7TPiogyXoTml1WxCj9aAj02Y8znU8aHJMtqAF4M9T6lV0J1O+13FaoysUAiSgleN++GwmvCNJ0+r8YYwE4fSL1KtRKPTWYNREYQAPS4rPA6LnqYZNIhGqbCirxfOLVV/oU0E68ia765459pu/Ol1G3HjtsE8R1YGsxxKGceiFopmWjR+7X840JV2KwHV668kSDKsZpIxT5tTGC4OZcKmvA122fE3N23BY8fmcc/+ScwuJ2AxEf0N3m6M9jqxGBPz9jSqFv/++Ckcml7O+9hiTES30wqL2QSPPVMczCai3wmXwiXDXbCaCfafbW/XUlSQ8mYqAerd81++/qKqVuczCyVfrUM4mcqw1BJibnNKVqDI+nn5HJlzHlZLMqXAwa2GkuHiUCZxgfWEseAPrhrDFWO9+MTPj2AiGMOAz9G2dyWsejRWo7GNikLxyQeO4p79k3kfX4gJuvB6HFY9/30hqsZ8yhlE5LCasX1NF15oc3GI5ZkCV0tYh9J8cYcP3XMA7/n2Pv33uCjBZc1cW0Dr/MqGE7G6g2paDvY2brFdbfiZKhNW2emymWEyEfzFjZuxFE/hoUMX2talBKRTDlk3zWrDgoXBApXLwaio31l6syyHclxKjN3revDy5FLe1t/tgjrLoX53yiwlNlscKKV49vQCjs9GdJdTXMztPxbQ/r9MHFhGVbUylpIphccbyoCLQ5mkxUG9q7lyfS829nsgKxRD3dUZZ9iMsItMrSwHlmYYjOQPEi8YRMBjt+j99+ejIgKelbuKZrO21wVBUqqaR99sRAUZnjq2gy/kVjq7EEconkJMlLGozdJIiLJeHc24cr0fV23w67EG5lbad2YRf/XDl1ct5NxyKA9+psokLma2GiaE4J2vWgugfYPRQFoMWfFStWGxnEIN8RZiol4o5XGks5WCEUHPkS8Hlt650jzqViaWZ9BPLSlkObx0Ph34Px9SU77ztbW/5ZIh/Od7r9QD6EwkvvL4afz4hUmcW2VnYB5zKA8uDmVidCsx3rJ7BKO9TlzW5A30VgN7vbEaiUOyiDikZAVL8ZQ+HcxjtyAmypAVqrqVKhhfyi6a7SwOS/GUfvddD5g4ZLfQeNFQU8JavyfyuJWyYWuXtME8S/HVtTzhlkN58CK4Msl2KwHqm/iJv7muUUuqC2nLoTYXU2Y5hOIppGQlo1YkpLkiWECaBSovhJMQJEWPRZQDa5ceq1EMpdHICsVirLJ4TKUUCki/eH4Jl4504cDkMs6HVHGI53ErZWO3mNSWHJo7KRRbnQtQ4JZDWXAZLZNst1KnoMccanQxNRY6ZbfTZjOfAwbLAUg3AqzkAujWLYf2ijk8e3oBU0sJLMVFKBQVxWMqJZ/lkEzJODwdxtUbA/C7bbrlEBeljBusfBBCsG2ND2/SGiourtJySHLLoSz4mSqTfG6lToC5AOI1shyEVFp0glEBd3zzOXz6waOglOpikU5lVS8qrGK7ErcSsz6ibWQ5UErx3rv34UuPndSzvio5N5XCMoGMlsOJC1FICsWOkS6M9Lr0GSjJlFJSq/p7/+dV+Mc3XQygCm4lbjmUBXcrlYluOTTxaM9awIqp4jUOSANqZ87Hj8/j8ePzcNst+oAfFpBmfYFY65JK7o71gHQVG7s1mog2hGoqlKiocny12PJkK7G5GX1eB0Z7nHhlahmSrECUlZJusAgh8NotsJgIQvHVWXnccigPfqbKRJ9g1UYtuUuBvd7aBaTTF5RnTi8AUNtz//NDx/RK5oDW83/X2h7YLCb87OVpAKgoW4m5pmpZ8V1vWH+v2eWkQRzq51ay58lWYgOiupxWjPa6MBVK6EkApVrfhBB0u2zccqgzXBzKJCHKcFhNbVsJXQiTicBpNa8qIP35R07gXd/Ym/cxo+XwzClVHD52q+pOuPfFKVhMRE9tdNsteM2mAMJJCYRAz2IqB2YJtVO2EhOHmeVE2q3UAMtBMIzjZHUkPqcFoz0uSArFaS1WVMooXUaPy6rXSFQKtxzKg5+pMomVEEhrV9x286osh5cnl7B3YhGKktuYLZnlVnJYTXj1ej8uGvAikpTQ67ZlNJC7abvaMK7HZYOlgi64JhOBy2ZuS3EIJyWcX4zDYiIZ40BrTb5sJaPlsFabwHZ8NgKgPOu7x21btVtJKDHOwVHh4lAmcTG3eKdTcNrMem+pSliIChAlJW8tAxMHZgVs7PfAZCK4fms/AOQ0NLxh6wDMJlJRGivDY7e0l1spnB5Te3BqGX6Pra5zCthdebY4OKwm2C1mXRyOauJQTlJHj8u6KrcSpVS1HKo09a4T4GeqTBKiXNd+Nc2E22ZZVUCapaSyKlkjTBxY8HlzvzqykQ3xyfad97htuHHrgD5XoxKMQ4PaAeMM88Mz4bq6lACD5SBnigOzXtZ0O2AxERzRJihmt+wuRo/LhsVV1DmkZApKwS2HMuhM/8gqiIm5feg7BafNvDpx0Pzgk6E4xvwuzCwncfFwFwA15mA2EQx1OXBgchmbtHm+O0e70e+1Y01Xbt+qL/3eLqwm9ONxtJflcCGchFcTvLgo118c8gSkwwlJFweL2YTRXpcuDmVZDm41IF3pVLikFgfhlkPp8DNVJglR0juUdhqq5VDZxTQuSnrQeTKUwGcePoZ3fO1ZvUtnMqXAaTXrF7TNA+qEMrOJ4MfvuwofvmVLznOaTWRVbhO3rb3EYWY5iUtGuvTfGyUOgpTfcgDUeeSsy2pZMQeXOhWu0hiRwKbAdehntxK4OJRJTOhct5JrFZaDsep5MpTAi+eWEElK+mzoRErNAkuLg1fff7TXlTFjuFp4HBZ9LkQ7cCGcxFjArY/brGcaK5B2KxUThzG/W/+5nGylbpf6WpYqDEoztyW3HEqnM/0jqyDfeMNOwWUzV9yy25iGeHIughNzUQDA6Xl1SFIyJcNuMeONO4YgK1SPPdQStYFfe4hDMqW2wx7yOTDY5UQonqq75UAI0UeFMpYTKWwZSgv9mN+l/1xeQFoVh8WYiNFe1wp758IEi8ccSmdVMkoIOUMIeYUQ8hIhZJ+2rZcQ8gtCyAnte49h/48QQk4SQo4RQm4ybN+tPc9JQsgXSD1TLMok3sFuJZfdUnHLbiYOw91O7D8bgqyls57RWmAkU2qXzo39XnzwpovqkmXjtpvbpkJ6LqxaYANdDr11fMBbX8sBAOxmU1bMIcutFEhbDtmT4IrR61afI1RhxhKzHBzcciiZapyp11FKd1JK92i/fxjAo5TSTQAe1X4HIWQbgNsBbAdwM4AvE0LYVfYuAHcC2KR93VyFddWEuCDD1aluJau54sZ7LH310pEusDIHQtLN89ReO/X94Hrs1rbpysrSWIe6HBhk4tCAeeY2iwmirJ5TWaGICFJG2/Bxg1vJYSv9/71atxIrzOMxh9KpxafxVgB3az/fDeBNhu3fp5QKlNIJACcBXEEIGQLgo5Q+Q9Xo5LcNxzQVlFLEUyu3Gm5XXHYLEik5bxHbSjDLYYc286LPa8eGPo9eLZtoQP2Ix26GKCsZFb2tysyymh486HNgSBuz6Xc3SBw0yyFsKIBjDPc4YTYRmE1Ej1GUgtGtVAksIM0th9JZ7ZmiAB4mhOwnhNypbRuglM4AgPa9X9s+DOC84dhJbduw9nP29qZDlBXICu3cCmlNFI2tLkplISbCbjHhIi3QfMlwF8YD7rTlIMl19wen+yu1gziolsNglwM713bD77ZhtLf+Y2ttFpPu31/OIw5WswkjPU64rOayXIddTisIqbwza5JbDmWz2qvc1ZTSaUJIP4BfEEKOFtk33zuBFtme+wSqAN0JAGvXri13rasmLnRmu25GehqcpHc1LZWFqIiAx65fsC4Z7kIyJePxY/OQFYpEA/Ly3Ybme5X0Z2omzi7E0Ou2weuw4tpNfdj/0Rsbsg67xaTfpecTB0DNWCo3dmXWWoFU2kJDtxx4b6WSWdWZopROa9/nANwL4AoAFzRXEbTvc9rukwBGDYePAJjWto/k2Z7v732VUrqHUrqnr69vNUuviHiq08Wh8jnSCzEBvW4bxgMe/NFr1uNtu0cwFnBDlBVMLyUgSPXve8NmOrRDOutEMJaRCdQoRntceOjwLP7iBy/pcaYuV6Y4vGXXMN66eyTf4UUZ9Dn0BIZySRfBdeZntxIqFgdCiJsQ4mU/A3g9gIMAfgrgDm23OwDcp/38UwC3E0LshJBxqIHn5zTXU4QQcqWWpfT7hmOaCpbZ0qluJd1yqMANsxgT4ffYYDYRfOQNWzHa68K4lrlyZiGmdrutsz9YtxzaIJ31TDCOMUMmUKP47Nt34g+uGsO9L07h3henAORaDrfuHMaHbs4talyJK9f78fyZxYpiRNxyKJ/VnKkBAE8SQl4G8ByA+ymlDwL4JIAbCSEnANyo/Q5K6SEAPwRwGMCDAN5PKWX/5fcB+DrUIPUpAA+sYl01Yy6i+nX76zhdq5lgA44SqfIvpgtRMcd1w8RhIhhDUlp54Hy1adWBP5TSjC62CVHGbDiJ9U0gDl0uKz58yxbYLSY8fnxe3ValzrBXbwwgmVLwwtmlso6776UpXNBSfbnlUDoV3wJTSk8D2JFn+wKA6wsc8wkAn8izfR+AiytdS72YNQT9OhF3hZYDpRQLMSGng2qfxw6LiWBmOYlkqv4Baa+99WY6/Mfes/jHnx2GICn49Nsuxdv3jOqulmawHAD1ArxjpBvPnVkEUD1xeNX6XphNBE+fCuLVG/wF96OU4p79k/itS9cgFBfxZ99/Sc+E45ZD6fAzVQZMHAZ8nSkO6TnSpYvDS+eX8N1nzyKZUnLabptMBH1eOy6EkyXPFK4m7habBhdJpvCpB45iy5APAz47Hjo4CyBdK2JsTdFo9oypta82i6lq/1efw4odI1148mSw6H4n5qL463sO4IGDM5g3tGcBuOVQDlwcymA2nESPy9qxJfjpOdKlX0z/6edH8NH7DgHIP7Gt32vHpDZ0vu5FcI7Wshy+++w5hJMSPnbrdly/dQB7JxYhyQommsxyAIDLx3oBIKMArhpcvTGAA5PLWC6StcQC4XMRIaMuwmomHTfBcTVwcSiD2eUkBvO0ju4UWGV4qdPgKKU4fiGCazYG8N5rx3GDNpvBSJ/XgXOLcQD1n8vdKqNCzwRj+NSDR/GVX5/CtZsCuHSkG1dt8CMqSDgwtYwzwRj6vHa9bqMZ2LW2B4QAXc7qrumWi4egUIrPPHys4D4hbe5DMCLoQrGmy8HnR5dJx4qDJCslVfqmZAVff+I0kik16Dfo68xgNGBMZS3tYhqMigjFU3jdln7879/clt9y8Nn11g/1tsjMbFRokwek//WR4/jK46ewzu/G3/7mNgDAq9erPvdnTi3gTDCe0ZaiGehyWXHRgLfq9SPb1vjwh1eN4zvPnsVzE4t591nUCuUWYqJuOXzxd3fh429u+rBmU9E8txp1RFYobvjs4xAlBbdfsRYfuG5jwWrNZ04t4OP3H0HAY8fschKXjnTXebXNA7uzLzUgfeKCOg6SzWbIhzHzqxHBwi6nFUuJ1c0mrjUzy0nsXteDH/3xVfo2v8eOrUM+/OzlacwsJ3HT9lyrrNH889ty8lWqwgdv2oz7X5nG1544jSvGe3MeX9IEIRgVsBCzw2YxYdfabuxe15OzL6cwHWk5vHQ+hDMLcRBC8NlfHMe0YbxiNiwT5IVzISzERAx2aDAaUO+0HVZTyTGH45o4XGSYzZBNvzd9Phsxm7vPa9eDls3KfFTIOE+M12wK4OhsBIpCceO2wQasrDiXjHRlDB+qFi6bBZsHvHlnkQNpy2E+IqiV+e76ztJuFzpSHB4+dAFWM8Gf3bAJAIoGt84EVX/4o0fUQu+hDk1jZawPePDLo3N6y+1iHJ+LostpRV+RukvDFlAAABKESURBVBCj5dCIvjd9nhYQh7CQ9xx+4PpN+PH7Xo39H70RN25rPsuhlnS7bAgVaMLHOrcGoyIWY0JOlhynNDpOHCileOjQLF69IYCRbjW4HE4WFoezmuUwtaRm1Ax0uDh84LqNODUf06tfi3HiQgSbBzxF79r6DTGcRlgOAY+94B1oM5AQZUQEKa84eOwW7F7Xq4/n7CR6XJl9ls4vxrHn44/g2GxEjzMsxgTMR4WW75vVKDruXXVyLoozC3G8ftsAfFpxTriIz3liIQaLIf2t0y2Hmy8exMXDPvzrI8chyUrB/dRMpSg2FXEpAZlupUakCPd57ViIiRW1Ia8HzKopZn11Ij0uG8LJlP4efOpkEMGogJfOh/SBQApVJw366zwutV3oOHF4+PAFAMCN2wb0HOxwgWwVWaE4vxjHtZsC+rZOLYBjEELwB1eNYzKU0Gcx5OPA5DKWEyls7i8cjAbUOcfMsGhUzEFWaMUTxmrNfLSzW7YUosdlBaXpzq8vTy4DAKZCCYTiom5NxUU5pzKfUxodJw637RnBl39vFwZ8Dvi0HOxClsP0UgIpmeL6rQNwWE1w2czwOToywSuDTdoF/0wBcfjoTw7i1i89BbOJYNcKGSIWs0n/8DYiW4m1CZ9vUtcSG//JLYdMerT3DHMtHZhU+y1NhhIIxVLY0Je+KeExh8roOHHo9zrwhkuGAKSHvRSKOZxdUIPRG/o82Dbkw6DPwbMekG7TwM5PNvvOhrBztBtPfei6klJ/+zTXUqPcSgCaNig9p60rX7ZSJ9Ojjw0VkUzJODarZsZNLMQQFaSM9Gkec6iMjr4NtphN8NgtummaTbotgQsfecPWpi+WqhddLiu6XdaCvfUXogKu29JfcoPCfq8dR2YaIw4BzR/drEHp+YgAs4nwC1wWxrGhR2bCkBQKr8OCIzNhAMBmQ6wrwGMOFdFxlkM2XU4rwon8F/2zwRgcVhMGvA5cPtaL123pz7tfJ7LO785rOSgK1Wc3lArzpzfCrVRPyyElKzg4tVzWMXORJPxuG+8JlEW3NkBoKZ7CK9o5vXHbAJLa3IZ1fpeeSNKIWdrtQMeLg9dhyetWopTilallrOt1w8Q/mDmM+V15LYdwMgVJoWV9IMcCbvgclrIGzlcLj90Ch9WEYLT2Aen/emESb/y3JzGznCj5mPmIkJHuy1FJxxxEvHx+GQGPTW/2B6iuJHaDwq2uyuh4cfA5rTkBaUop/uFnh7F3YhG/vXNNg1bW3Iz53dp4z8xWGuwiW47l8O5rxvHAn7+mIfEcQggCdSqEOzkXBaXAxHzpoy7nIgL6eEA1B7fNDJvZhFA8hSMzYWxb04Xh7nRTzB6XTU824KmslcHFwWHNSWX9xeEL+NbTZ/Cea8bxP1+7oUEra27GAi4oVM0OMbKg+e7LsRwcVnPGB7ve1KuFBus+ez6kfi9lFvd8JH/rjE6HEIJulxULUQGng1Fs6vdgpCf9Hup1q+LgtJo7dqzvauHi4LTkWA7/9cIUAh47PnzLFp6dVIB1esZS5l3wQqx8y6HR9NWpSprFaM4txvHUySB2/MPDOS6myVAc1//Lr3BqPgpZoQhGuVupED0uGw5Nh5FMKdjY78Eaww1Gt8uK9X1ufRQtp3w6XlJVyyEtDsuJFH55bA6/e8VaWBrgA28VWDor6z3F0C2HFhKHgNeO/WdDNf0blKoFlQBwbjEBM1mAKCs4NhvBkGFGyMOHLuDUfAx7Ty/C67BAobzGoRDdLiue10aRbuz3wGE1I+CxIyFKsFvM+NDNWyBIhav4OcXpeHHocloRFSQoCoXJRPDQwVmIkoI3XTbc6KU1NT0uK7wOS05QmsUcel2tIw59HjsW4yIkWanZDcFiTNSHJJ1bjOsupWy33BMn5gEAp+ejmAqp6Zid3Am4GL1uG1jXk41a0dtIj1O3Ah1Wc8dObawGHS8OPqdahh9JSuhyWfGzA9MY87uwowathtsJQgjW93lwci6asX0hJqDHZW0pq6vPawelapX0UI0m/Z3VrIYBnx3nF+MIajEOozgIkoxnT6t3whPBmF7YddFg8f5UnUq3dgPid9v07KXXXtSnz3rnrA4uDo50lXSXy4ojMxHcsLWfxxpKYHO/B48dm8/YthAVW65dAXORTQRjNRMH5lK6ZmMffvzCpL59MpR2y71wdgmJlAyfw4KJYAxHZyNw2cwY7XHVZE2tTo9W62BslfHnN2xu1HLajta5vasRrDPrciKFmCAhGBWw1s8/jKXABq4Y++ovRMWWa3S2oV8Vh9NlpJiWyzktGH3VBr++zWYxZVgOT56ch9lE8JZdIzi3GMeh6WVsGvDyOpsCsPqFDSs0d+RUBhcHvTNrSk81XNvLxaEUNmn9a9jENwAIxgQ9v7xVGPQ54LKZcWo+uvLOFXJ2MY4Bnz2jrcPVG/wZ4vDEiSB2re3GJcNdkBSK/WdD2LJCy/NOhrmVNnJxqAlcHPTOrJIuDut6efpbKbAL3XFD3EF1K7WW5UAIwYY+D07VwHKYCMbwqQeP4qBWbc9uPHrdNuxa24NgVEAyJSMUE/HK1DKu2diH8T71/adQHm8oBnufbeLiUBM6PubQ5UxbDmxcKHcrlcZQlwNeuwUnLkTw3WfP6k0MW7GXzfo+N/adqX4665ceO4l79qsxhrfuGkGXywqfw4Itg16MakIxGUrg6GwYlALXbg5gvSE3fwsXh4JcszGAT7/1Uly9MbDyzpyy6XhxME6DO7sYQ7fLqgsGpziEEGwc8OBXx+bx3WfP6kH8VrMcADWoed9L00iIMpy26qQ/ipKChw/N4qbtA7hi3K8PjfrT6zdhzO/Wm8dNhuJ48kQQXocFlw53wWI2oddtw2JM5JZDEaxmE95++Wijl9G2dLw4eGwWEKJOgzu7EOfxhjLZ3O/FD86d16e4JRS5JVsks4yX08Eotq+pPI15ZjmBfq8DZhPBU6eCCCclvH3PKK7fOqDv855r1wOAnnJ5PpTAEyeCuHpDQE8BHg+4YSKk5TK/OO1Dx8ccTCYCr11toXF+kYtDubCg9B1XjeH9r1P7ULXiKFWWscTiDrJCy54rfWh6Gdd+6jF86McHAAA/PzADr92Cazbld3v0e+2wmU342cvTmFpKZOz3nmvG8Rc3bqrkpXA4VaHjLQdALYJ68fwSJkMJ/OalQ41eTktxw9YBPH9mEX/0mvXwOiy4eLgLO0dXnv7WbIz53SAEOKUF19/1jb0Y9Dnw2d/Zic8/cgJLCRH/57e2Fax/kRWK/3XvQSiU4p79k/DYLbj/lRncvH0Qdkt+N5XJRDDc48RzE4tY2+vCTdsH9cduuYS/DzmNhYsDgPdeux4f/q9XAPBMpXIZC7jxlXft0X9/7UWtORDJYTVj3O/GC+dCmF5K4OlTC7BZTPjQLVvw5V+dhCAp8Ltt+JPr8t/N/3Dfebx8fgmfuW0HvvPMGXzr6TPYtbYbf/n64kVZH/2trQhGRbxp5zBslo435DlNBBcHALftGcV3nj2LQ9NhPYOE03nccskg7vrVKfzH3rMA1IDy39xzAIKk4LK13fjMw8exZ6wXO0a6cdfjp/Dmy4YxHnCDUopvPjmBS0e68NZdw3jN5gBePLeEG7cOrFjAdt2WgaKPcziNgt+qADCbCP7vmy/BFeO92D7sa/RyOA3iLbtGoFDg3x8/jfUBNwZ9Djx+fB4jPU7853uuxHC3Ex+//zD+5eFj+MKjJ/DWu57GgcklvHh+CSfmonjHFWtBCEG/14Gbtg/yymZOS9M04kAIuZkQcowQcpIQ8uF6//0do9344R+9Wq+Y5nQeG/o8uGxtN2SF4qaLB3HLJWoM4Ld3rIHTZsZf33QRDk6F8fUnJ3DT9gG4bGa846vP4v/efwQumxlv3MGnBnLah6YQB0KIGcCXANwCYBuAdxBCtjV2VZxO5O171Lz5m7cP4u17RjHa68Rt2rbf3rEGO0a7Mehz4NNv24Efv+8qjPa6sO9sCL95yRA8du6l5bQPhNLy0vVqsghCXg3g7ymlN2m/fwQAKKX/VOiYPXv20H379tVphZxOQVEojs5GsG1NfvdiXJQgSore12c5kcKXHzuJd165jserOC0BIWQ/pXTPSvs1y63OMIDzht8nAbyqQWvhdDAmEykoDADgsllgnGPU5bTiI2/YWoeVcTj1pSncSgDyRe5yTBpCyJ2EkH2EkH3z8/N5DuFwOBxONWgWcZgEYGySMgJgOnsnSulXKaV7KKV7+vr66rY4DofD6TSaRRyeB7CJEDJOCLEBuB3ATxu8Jg6Hw+lYmiLmQCmVCCF/AuAhAGYA36SUHmrwsjgcDqdjaQpxAABK6c8B/LzR6+BwOBxO87iVOBwOh9NEcHHgcDgcTg5cHDgcDoeTQ1NUSFcCISQC4Fieh7oALNfgT1bzeQMAglV6LqD6r7nZn4+fv8rh5251tPL5Y2tfRylduRaAUtqSXwD2Fdj+1Rr9vao9b6G1N8PaWuT5+Pnj565Rz9ey56/ctbejW+lnLfa81aDaa2v256s2zf56m/n8NftrbeZzBzTx621lt9I+WkLzqGakldfeDPDzVzn83K2OVj5/5a69lS2HrzZ6AaugldfeDPDzVzn83K2OVj5/Za29ZS0HDofD4dSOVrYcOBwOh1MjuDhUAULIKCHkMULIEULIIULIn2nbewkhvyCEnNC+92jb/dr+UULIv2U91+8QQg5oz/PpRryeelPB+buRELKfEPKK9v06w3Pt1rafJIR8gRDS1oOcq3zuPkEIOU8IiTbq9dSbap0/QoiLEHI/IeSo9jyfbOTrqgrVTKPq1C8AQwB2aT97ARyHOu700wA+rG3/MIBPaT+7AVwD4I8B/JvhefwAzgHo036/G8D1jX59TXj+LgOwRvv5YgBThud6DsCroc4IeQDALY1+fS107q7Uni/a6NfVaucPgAvA67SfbQCeaPX3XsMX0I5fAO4DcCPUIr0hbdsQgGNZ+/1BljhcDuARw+/vAvDlRr+eZj1/2nYCYAGAXdvnqOGxdwD4SqNfTyucu6ztHSMOtTh/2mOfB/DeRr+e1Xxxt1KVIYSMQb272AtggFI6AwDa9/4VDj8JYAshZIwQYgHwJmQOQWp7Kjh/bwXwIqVUgDpudtLw2KS2rSNY5bnreKp1/ggh3QDeCODRWq631jRNy+52gBDiAfBjAH9OKQ2X6+6mlIYIIe8D8APg/7d3PyFWVXEAx78/0o0UlaAgRESbqCQMhYoKd0GtWxTRTEaLSIh2YQi1qIWhIupikjL6h0RQMBUUISRYtBTTZiEjgoEoUZmNS38tznn0mPte09P75j1nvh+4vOG8M4dzfzzmd8+Z+36XK8CPwJ2tT3RMDRq/iLgX2AE81mnq0W1Z3I7XQuyWtbbiVy/qDgF7M/P0kKa7KFw5tCQiVlI+XJ9k5ue1+XxErKvvrwMuLDROZn6ZmQ9k5kOUpe2pYc15nAwav4i4DfgCmMjM2dr8K+URsx09Hze71LQUu2Wr5fgdAE5l5p7hz3y4TA4tqHfEvAfMZOburremgcn68yRlP3OhsdbW11uBl4B3253t+Bk0fnXZ/jWwLTN/6HSuy/9LEfFgHXOC/xHz61lbsVuu2oxfRLxJKXz3yrDnvShG/U+PpXBQ7jxK4DhwrB5PUO4+Oky5+j8MrO76nTPA78DflCvee2r7IeCXejw16nMbx/gB24G5rr7HgLX1vU3ACWAW2E/9oudSPVqO3dv1s3ilvr4x6vO7XuJHWaUmMNPV/sKoz+9aDr8hLUlqcFtJktRgcpAkNZgcJEkNJgdJUoPJQZLUYHKQhiAiXoyIiQH63xERJ4Y5J2kQls+QWhYRKzJzatTzkK6FyUHqoRZh+4ZShO1+SinnCeBuYDdwI/Ab8FxmnouI7ym1sB4GpiPiJkp1050RsQGYopR1ngWez1JHayNwELgMHF28s5MW5raS1N9dwIHMvA/4C9gK7AOezMzOH/a3uvrfkpmbM3PXvHE+BF6t4/wMvF7b3wdezlJHSxorrhyk/s7mv/VzPgZeozzg5btatfMG4FxX/0/nDxARN1OSxpHa9AHwWY/2j4DH2z8F6eqYHKT+5teWuQSc/I8r/bkBxo4e40tjw20lqb/bI6KTCJ4GfgLWdNoiYmWt699XZl4E/oiIR2vTs8CRzPwTuBgRj9T2Z9qfvnT1XDlI/c0AkxHxDqU65z7gW2Bv3RZaAewBTi4wziQwFRGrgNPAltq+BTgYEZfruNLYsCqr1EO9W+mrzFw/4qlII+G2kiSpwZWDJKnBlYMkqcHkIElqMDlIkhpMDpKkBpODJKnB5CBJavgH8yH1B6I8C30AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sorted_data['inc'][-200:].plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**we need to start from Ist of september so below 8 has been replaced with 9\n", "Also the current data set is from 1990 to 2022, so we give 1991 the starting year to consider the year with full peak of the epidemic as described by the instructiors of the code\n", "--- Replaced 1985 by 1991**" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "first_august_week = [pd.Period(pd.Timestamp(y, 9, 1), 'W')\n", " for y in range(1991,\n", " sorted_data.index[-1].year)]" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "year = []\n", "yearly_incidence = []\n", "for week1, week2 in zip(first_august_week[:-1],\n", " first_august_week[1:]):\n", " one_year = sorted_data['inc'][week1:week2-1]\n", " assert abs(len(one_year)-52) < 2\n", " yearly_incidence.append(one_year.sum())\n", " year.append(week2.year)\n", "yearly_incidence = pd.Series(data=yearly_incidence, index=year)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAG61JREFUeJzt3X+MXeV95/H3x4wzA9QmYxgT/wDsKg6KISnEV7a7yXZL3NpmU2GjQjIbGkaLJTfApqRaCezACi24ElSrTWtFYbFCwUD44bpFeLt4yWA2qrd1bI8DKRjiziQQ4+Aww44DJpKnjP3dP+5z6zuT8cy58+vec+fzkq7uud85z3Of4+OZ731+nHMVEZiZmWUxrdoNMDOz/HDSMDOzzJw0zMwsMycNMzPLzEnDzMwyc9IwM7PMnDTMzCwzJw0zM8vMScPMzDJrqHYDxtsFF1wQCxYsqHYzzMxy5cCBA+9GRMtI+9Vd0liwYAEdHR3VboaZWa5I+lmW/Tw8ZWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmY1Znu90/wxQf30H38RLWbYnXIScOszmze1cn+N3vZ/EJntZtidajurtMwm6ouvWsnff2n/vX143sP8/jewzQ2TOPQpqur2DKrJ+5plHG33vJs9+1Xcc0Vc2maXvy1bpo+jTVXzGX3HVdVuWVWT5w0yrhbb3k2e2YTMxob6Os/RWPDNPr6TzGjsYHZM5rOWMYflKxSHp7C3XqrH+9+0McNyy7hy0sv5ol9h+kZIRmUf1DadO2nJqmVlmeKiGq3YVwVCoWo9N5T3e+fYNNzr/O9g7/gxIenaJo+jVWXfYw7v/DJYT+lmeXV4A9KJf6gNHVJOhARhZH28/AUo+vWm+WZ5z9stDw8lVTarTfLM39QstFy0kge/MrpXtmmtZdXsSVmk8MflGw0PKdhZmae0zAzs/HnpGFmZpllShqS/lTSQUmvSnpSUpOkWZLaJXWm5+ay/TdK6pJ0SNKqsvgSSa+kn22WpBRvlPR0iu+VtKCsTFt6j05JbeN36GZmVqkRk4akecCfAIWIuBw4C2gFNgC7ImIRsCu9RtLi9PPLgNXAtyWdlap7AFgPLEqP1Sm+DjgWER8Hvgncn+qaBdwNLAOWAneXJyczM5tcWYenGoCzJTUA5wBvA2uArennW4G1aXsN8FRE9EXEG0AXsFTSHGBmROyJ4uz7o4PKlOraDqxIvZBVQHtE9EbEMaCd04nGzMwm2YhJIyJ+Dvw34DBwFHgvIr4HXBgRR9M+R4HZqcg84K2yKo6k2Ly0PTg+oExE9APvAecPU5eZmVVBluGpZoo9gYXAXOBcSX80XJEhYjFMfLRlytu4XlKHpI6enp5hmmZmZmORZXjq94A3IqInIj4E/hb4N8A7aciJ9Nyd9j8CXFRWfj7F4awjaXtwfECZNAR2HtA7TF0DRMSWiChERKGlpSXDIZmZ2WhkSRqHgeWSzknzDCuA14EdQGk1UxvwbNreAbSmFVELKU5470tDWMclLU/13DioTKmu64AX07zH88BKSc2px7MyxczMxsS3hR+dEW8jEhF7JW0Hfgj0Ay8BW4DfALZJWkcxsVyf9j8oaRvwWtr/1og4maq7GXgEOBvYmR4ADwGPSeqi2MNoTXX1SroX2J/2uyciesd0xGZm+Lbwo+XbiNiU1/3+Cf7Tky/xrS9f6Rv2TQG+LfzQfBsRs4z8jY1Ti28LPza+y61NWf7GxqnJt4UfG/c0bMryJ86pq3Rb+Gdu+Sw3LLuEng/6qt2k3HBPw6Ysf+Kcuvz9OaPnpGFTmr+IyKwyXj1ldckroswq49VTNiHyckGUV0SZTQwPT1lFav2CKK+IMptYHp6yTPJyQVT3+yfY9NzrfO/gLzjx4Smapk9j1WUf484vfDL3w1QecrOJ5OEpG1d5WZ5azyuiPORmtcDDU5ZJnv4Y19uKKA+5WS1x0rDM8vLHuN7W4O++/aozDrmZTTYnDcus3v4Y50WeenlW/5w0zHIgL708q39ePWVmZl49ZWZm489Jw8zMMnPSMDOzzJw0zKooL/fyMitx0jCrIl/lbXnjJbdmVeCrvC2v3NMwq4K83MvLbDAnDbMq8FXellcjJg1Jl0p6uezxvqSvS5olqV1SZ3puLiuzUVKXpEOSVpXFl0h6Jf1ssySleKOkp1N8r6QFZWXa0nt0Smob38M3q57SVd7P3PJZblh2CT0f9FW7SWYjquiKcElnAT8HlgG3Ar0RcZ+kDUBzRNwhaTHwJLAUmAu8AHwiIk5K2gfcBvwAeA7YHBE7Jd0CfDoiviqpFbg2Ir4kaRbQARSAAA4ASyLi2Jna6CvCzfLB3w9SWybqivAVwE8i4mfAGmBrim8F1qbtNcBTEdEXEW8AXcBSSXOAmRGxJ4qZ6tFBZUp1bQdWpF7IKqA9InpTomgHVlfYZjOrQV45lk+Vrp5qpdiLALgwIo4CRMRRSbNTfB7FnkTJkRT7MG0PjpfKvJXq6pf0HnB+eXyIMmaWQ145lm+ZexqSPgJcA/z1SLsOEYth4qMtU9629ZI6JHX09PSM0DwzqyavHMu3SoanrgZ+GBHvpNfvpCEn0nN3ih8BLiorNx94O8XnDxEfUEZSA3Ae0DtMXQNExJaIKEREoaWlpYJDMrPJNpqVY75yvnZUkjT+A6eHpgB2AKXVTG3As2Xx1rQiaiGwCNiXhrKOS1qe5ituHFSmVNd1wItp3uN5YKWk5rQ6a2WKWQb+RbNaVenKMc9/1I5Mq6cknUNxbuE3I+K9FDsf2AZcDBwGro+I3vSzO4GbgH7g6xGxM8ULwCPA2cBO4GsREZKagMeAKyn2MFoj4qepzE3AN1JT/iwiHh6urV49ddpdz7zCd/cd5oalF7Pp2k9VuzlmFRs8/1Hi+Y/xl3X1lL+EqQ6N5hfNyx+tFnW/f+KM34/u/6fjy1/CNIWNZqLR3X+rRb5yvvb4hoV1qJJfNC9/tFrn70evLU4adSrrL9ru2686Y/ffrBY8+JXTIyab1l5exZYYOGnUray/aO7+m1klPKcxwfKw7NU3zjOzrLx6aoJ52auZ5UHW1VMenpognmA2s3rk4akJUumy1zwMY5mZOWlMkEonmH2dhJnlgYenJlCWZa8exjKzPPFE+CiM5y03fJsEywPfZqb++TYiE2g8h5J8nYTlgYdPrcTDUxWYqKEk3ybBapWHT20wD09VwENJllW9DOf4//zU4eGpCeChpKmtkmXR9TKc4//zNpiHpyrkoaSpqzwRnOnq/noczvH/eSvn4SmzEVTypVYezqmfobmpxsNTlhu1fjV8JVf3ezinfobmbGgenrIJk/UTZ5Zhn2qqNBFM1eGcehyas1/n4SmbMCPd4Xc032VeLX/8WActM5oGJILy7ywxD83lne9ya1WT9RNnnr410N8eNzIPzU0NntOwcZd1DsB/ZOqPv9Cr/rmnYeOukmQwVcf/61W1e2ReuTXxMvU0JH1U0nZJP5b0uqTfljRLUrukzvTcXLb/Rkldkg5JWlUWXyLplfSzzZKU4o2Snk7xvZIWlJVpS+/RKalt/A7dJlLWT5wPfqXAprWXs3juTDatvdzzBDYmXrk18TJNhEvaCuyOiO9I+ghwDvANoDci7pO0AWiOiDskLQaeBJYCc4EXgE9ExElJ+4DbgB8AzwGbI2KnpFuAT0fEVyW1AtdGxJckzQI6gAIQwAFgSUQcO1NbPRFuNvXkaVFFrRq36zQkzQR+B3gIICL+JSJ+CawBtqbdtgJr0/Ya4KmI6IuIN4AuYKmkOcDMiNgTxUz16KAypbq2AytSL2QV0B4RvSlRtAOrR2qzmU0tlX5Tpo1eluGp3wR6gIclvSTpO5LOBS6MiKMA6Xl22n8e8FZZ+SMpNi9tD44PKBMR/cB7wPnD1DWApPWSOiR19PT0ZDgkM6snXlQxebIkjQbgM8ADEXEl8CtgwzD7a4hYDBMfbZnTgYgtEVGIiEJLS8swTTOzeuWVW5Mjy+qpI8CRiNibXm+nmDTekTQnIo6moafusv0vKis/H3g7xecPES8vc0RSA3Ae0JvivzuozPczHZmZTSnVXrlViTyv8hqxpxERvwDeknRpCq0AXgN2AKXVTG3As2l7B9CaVkQtBBYB+9IQ1nFJy9N8xY2DypTqug54Mc17PA+slNScVmetTDEzs9zK8yqvrNdpfA34blo59VPgP1JMONskrQMOA9cDRMRBSdsoJpZ+4NaIOJnquRl4BDgb2JkeUJxkf0xSF8UeRmuqq1fSvcD+tN89EdE7ymM1M6uqerg/l+89ZWY2SWr5/ly+NbqZWY2ph1Vevo2Imdkkyvutczw8ZWZmHp6qZ7X+TXdmVr+cNHIoz8v1zCzfPKeRI/WwXM/M8s09jRzxTdnMrNqcNHKkHpbrmeWN5xAHctLIGd+UzWxyeQ5xIC+5NTMbwlT7YicvuTUzG4PRzCFOhaEsJw0zsyGMZg5xKgxlecmtmdkZZL3lx1RaDu85DTOzMarlu9dm5TkNM7NJMpWWw3t4ysxsHOT97rVZeXjKzMw8PGVmZuPPScNyZSqsgzerZU4alitTYR28WS3zRLjlwlRaB29Wy9zTsFzwbeHNaoOThuXCVFoHb1bLMiUNSW9KekXSy5I6UmyWpHZJnem5uWz/jZK6JB2StKosviTV0yVpsySleKOkp1N8r6QFZWXa0nt0SmobrwO3/PFt4c2qL9N1GpLeBAoR8W5Z7M+B3oi4T9IGoDki7pC0GHgSWArMBV4APhERJyXtA24DfgA8B2yOiJ2SbgE+HRFfldQKXBsRX5I0C+gACkAAB4AlEXHsTG31dRpmZpWbjOs01gBb0/ZWYG1Z/KmI6IuIN4AuYKmkOcDMiNgTxUz16KAypbq2AytSL2QV0B4RvSlRtAOrx9BmMzMbg6xJI4DvSTogaX2KXRgRRwHS8+wUnwe8VVb2SIrNS9uD4wPKREQ/8B5w/jB1DSBpvaQOSR09PT0ZD8nMzCqVdcntZyPibUmzgXZJPx5mXw0Ri2Hioy1zOhCxBdgCxeGpYdpmZmZjkKmnERFvp+du4BmK8xXvpCEn0nN32v0IcFFZ8fnA2yk+f4j4gDKSGoDzgN5h6jIzsyoYMWlIOlfSjNI2sBJ4FdgBlFYztQHPpu0dQGtaEbUQWATsS0NYxyUtT/MVNw4qU6rrOuDFNO/xPLBSUnNanbUyxczMrAqyDE9dCDyTVsc2AE9ExP+WtB/YJmkdcBi4HiAiDkraBrwG9AO3RsTJVNfNwCPA2cDO9AB4CHhMUhfFHkZrqqtX0r3A/rTfPRHRO4bjNTOzMfCt0c3MzLdGNzOz8eekYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZWWZOGmZmlpmThpmZZeakYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpllThqSzpL0kqS/S69nSWqX1Jmem8v23SipS9IhSavK4kskvZJ+tlmSUrxR0tMpvlfSgrIybek9OiW1jcdBm5nZ6FTS07gNeL3s9QZgV0QsAnal10haDLQClwGrgW9LOiuVeQBYDyxKj9Upvg44FhEfB74J3J/qmgXcDSwDlgJ3lycnMzObXJmShqT5wBeA75SF1wBb0/ZWYG1Z/KmI6IuIN4AuYKmkOcDMiNgTEQE8OqhMqa7twIrUC1kFtEdEb0QcA9o5nWjMzGySZe1p/AVwO3CqLHZhRBwFSM+zU3we8FbZfkdSbF7aHhwfUCYi+oH3gPOHqcvMzKpgxKQh6Q+A7og4kLFODRGLYeKjLVPexvWSOiR19PT0ZGymmZlVKktP47PANZLeBJ4CPi/pceCdNOREeu5O+x8BLiorPx94O8XnDxEfUEZSA3Ae0DtMXQNExJaIKEREoaWlJcMhmZnZaIyYNCJiY0TMj4gFFCe4X4yIPwJ2AKXVTG3As2l7B9CaVkQtpDjhvS8NYR2XtDzNV9w4qEypruvSewTwPLBSUnOaAF+ZYmZmVgUNYyh7H7BN0jrgMHA9QEQclLQNeA3oB26NiJOpzM3AI8DZwM70AHgIeExSF8UeRmuqq1fSvcD+tN89EdE7hjabmdkYqPiBvn4UCoXo6OiodjPMzHJF0oGIKIy0n68INzOrA93vn+CLD+6h+/iJCX0fJw0zszqweVcn+9/sZfMLnRP6PmOZ0zAzsyq79K6d9PWfvoTu8b2HeXzvYRobpnFo09Xj/n7uaZiZ5dju26/imivm0jS9+Oe8afo01lwxl913XDUh7+ekYWaWY7NnNjGjsYG+/lM0Nkyjr/8UMxobmD2jaULez8NTZmY59+4Hfdyw7BK+vPRinth3mJ4JnAz3klszM/OSWzMzG39OGmZmlpmThpmZZeakYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZWWZOGmZmlpmThpmZZeakYWZmmY2YNCQ1Sdon6UeSDkr6ryk+S1K7pM703FxWZqOkLkmHJK0qiy+R9Er62WZJSvFGSU+n+F5JC8rKtKX36JTUNp4Hb2ZmlcnS0+gDPh8RvwVcAayWtBzYAOyKiEXArvQaSYuBVuAyYDXwbUlnpboeANYDi9JjdYqvA45FxMeBbwL3p7pmAXcDy4ClwN3lycnMzCbXiEkjij5IL6enRwBrgK0pvhVYm7bXAE9FRF9EvAF0AUslzQFmRsSeKH7H7KODypTq2g6sSL2QVUB7RPRGxDGgndOJxszMJlmmOQ1JZ0l6Geim+Ed8L3BhRBwFSM+z0+7zgLfKih9JsXlpe3B8QJmI6AfeA84fpi4zM6uCTEkjIk5GxBXAfIq9hsuH2V1DVTFMfLRlTr+htF5Sh6SOnp6eYZpmZmZjUdHqqYj4JfB9ikNE76QhJ9Jzd9rtCHBRWbH5wNspPn+I+IAykhqA84DeYeoa3K4tEVGIiEJLS0slh2RmZhXIsnqqRdJH0/bZwO8BPwZ2AKXVTG3As2l7B9CaVkQtpDjhvS8NYR2XtDzNV9w4qEypruuAF9O8x/PASknNaQJ8ZYqZmVkVNGTYZw6wNa2AmgZsi4i/k7QH2CZpHXAYuB4gIg5K2ga8BvQDt0bEyVTXzcAjwNnAzvQAeAh4TFIXxR5Ga6qrV9K9wP603z0R0TuWAzYzs9FT8QN9/SgUCtHR0VHtZpiZ5YqkAxFRGGk/XxFuZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZWWZOGmZmlpmThpmZZeakYWZmmTlpmJlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpk5aZiZ1bDu90/wxQf30H38RLWbAjhpmJnVtM27Otn/Zi+bX+isdlMAaKh2A8zM7NddetdO+vpP/evrx/ce5vG9h2lsmMahTVdXrV3uaZiZ1aDdt1/FNVfMpWl68c900/RprLliLrvvuKqq7XLSMDOrQbNnNjGjsYG+/lM0Nkyjr/8UMxobmD2jqart8vCUmVmNeveDPm5YdglfXnoxT+w7TE8NTIYrIobfQboIeBT4GHAK2BIRfylpFvA0sAB4E/hiRBxLZTYC64CTwJ9ExPMpvgR4BDgbeA64LSJCUmN6jyXA/wO+FBFvpjJtwF2pOZsiYutw7S0UCtHR0ZH9X8DMzJB0ICIKI+2XZXiqH/jPEfFJYDlwq6TFwAZgV0QsAnal16SftQKXAauBb0s6K9X1ALAeWJQeq1N8HXAsIj4OfBO4P9U1C7gbWAYsBe6W1JyhzWZmNgFGTBoRcTQifpi2jwOvA/OANUDpU/9WYG3aXgM8FRF9EfEG0AUslTQHmBkRe6LYvXl0UJlSXduBFZIErALaI6I39WLaOZ1ozMxsklU0ES5pAXAlsBe4MCKOQjGxALPTbvOAt8qKHUmxeWl7cHxAmYjoB94Dzh+mLjMzq4LMSUPSbwB/A3w9It4fbtchYjFMfLRlytu2XlKHpI6enp5hmmZmZmORKWlImk4xYXw3Iv42hd9JQ06k5+4UPwJcVFZ8PvB2is8fIj6gjKQG4Dygd5i6BoiILRFRiIhCS0tLlkMyM7NRGDFppLmFh4DXI+K/l/1oB9CWttuAZ8virZIaJS2kOOG9Lw1hHZe0PNV546AypbquA15M8x7PAyslNacJ8JUpZmZmVZBlye3ngN3AKxSX3AJ8g+K8xjbgYuAwcH1E9KYydwI3UVx59fWI2JniBU4vud0JfC0tuW0CHqM4X9ILtEbET1OZm9L7AfxZRDw8Qnt7gJ9lPP5acwHwbrUbMc7q7Zjq7Xig/o6p3o4HJueYLomIEYdqRkwaNnkkdWRZJ50n9XZM9XY8UH/HVG/HA7V1TL6NiJmZZeakYWZmmTlp1JYt1W7ABKi3Y6q344H6O6Z6Ox6ooWPynIaZmWXmnoaZmWXmpDHBJP2VpG5Jr5bFfkvSHkmvSPqfkmam+EckPZziP5L0u2Vlvi/pkKSX02P2EG834SRdJOn/SHpd0kFJt6X4LEntkjrTc3NZmY2SulL7V5XFl6Rj7ZK0OV2/k+fjyeU5knR+2v8DSd8aVFfuztEIx5PXc/T7kg6kc3FA0ufL6prccxQRfkzgA/gd4DPAq2Wx/cC/S9s3Afem7VuBh9P2bOAAMC29/j5QqIHjmQN8Jm3PAP4ZWAz8ObAhxTcA96ftxcCPgEZgIfAT4Kz0s33Ab1O8XcxO4OqcH09ez9G5wOeArwLfGlRXHs/RcMeT13N0JTA3bV8O/Lxa58g9jQkWEX9P8YLFcpcCf5+224E/TNuLKd5mnojoBn4J1MTa7JKYnLseT5rxOp7JbfXwKj2miPhVRPxfYMA3/OT1HJ3peGrJKI7ppYgo3ULpINCk4l03Jv0cOWlUx6vANWn7ek7fX+tHwBpJDSregmUJA++99XDqUv+XagwTDKaJu+txVYzxeEryeI7OJK/naCR5P0d/CLwUEX1U4Rw5aVTHTRS/zOoAxa7pv6T4X1E86R3AXwD/SPFWLAA3RMSngH+bHl+Z1BYPoom96/GkG4fjgfyeozNWMUQsD+doOLk+R5Iuo/gldX9cCg2x24SeIyeNKoiIH0fEyohYAjxJcVyciOiPiD+NiCsiYg3wUaAz/ezn6fk48ARVHBLRxN/1eFKN0/Hk+RydSV7P0Rnl+RxJmg88A9wYET9J4Uk/R04aVVBasSFpGsXvP/8f6fU5ks5N278P9EfEa2m46oIUnw78AcUhrmq0fTLuejxpxut4cn6OhpTjc3SmenJ7jiR9FPhfwMaI+IfSzlU5RxM5y+5HQLEncRT4kOKngnXAbRRXS/wzcB+nL7JcAByiOCn2AsW7TkJxNcgB4J8oToL9JWnFThWO53MUu7//BLycHv+e4jct7qLYM9oFzCorcyfF3tQhylZ2UJzkfzX97Fulf4c8Hk8dnKM3KS7Y+CD9P12c83P0a8eT53NE8cPlr8r2fRmYXY1z5CvCzcwsMw9PmZlZZk4aZmaWmZOGmZll5qRhZmaZOWmYmVlmThpmZpaZk4aZmWXmpGFmZpn9f/2Yi6a8X2tzAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ " yearly_incidence.plot(style='*')" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "2020 221186\n", "2021 376290\n", "2002 516689\n", "2018 542312\n", "2017 551041\n", "1996 564901\n", "2019 584066\n", "2015 604382\n", "2000 617597\n", "2001 619041\n", "2012 624573\n", "2005 628464\n", "2006 632833\n", "2011 642368\n", "1993 643387\n", "1995 652478\n", "1994 661409\n", "1998 677775\n", "1997 683434\n", "2014 685769\n", "2013 698332\n", "2007 717352\n", "2008 749478\n", "1999 756456\n", "2003 758363\n", "2004 777388\n", "2016 782114\n", "2010 829911\n", "1992 832939\n", "2009 842373\n", "dtype: int64" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "yearly_incidence.sort_values()" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW8AAAEICAYAAACQzXX2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAEaZJREFUeJzt3XuQZGV5x/Hvwy7owgCiC0MFMEPQGImrImPwUsFZsCwV1BJNvCCBlGZNeSNmUxZJqWgSlcRgaSxIakWUEuOoYEoB46Wio+IF2VXLFQlqBBGIRGJAl6Cw8uSP9wyZXXd2zlx6pp/N91M1Nae7T/d5nj7dv3P67dPdkZlIkmrZa6ULkCTNn+EtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJU0OpB3fDatWtzbGxsh/PuvPNO9ttvv0EtclnYw3Cwh+FgD0tvy5Ytt2XmwXPNN7DwHhsbY/PmzTucNzU1xcTExKAWuSzsYTjYw3Cwh6UXET/oM5/DJpJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUN7EM6kn7V2FlX3De9cd12zphxepBuOOekZVmOlo973pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUUO/wjohXR8Q1EfGtiPhARNx/kIVJkmbXK7wj4jDgVcB4Zj4CWAU8f5CFSZJmN59hk9XAmohYDewL3DKYkiRJc4nM7DdjxJnAm4C7gE9l5qm7mGcDsAFgdHT02MnJyR0u37ZtGyMjI4uteUXZw3Co2sPWm++4b3p0Ddx61/Isd91hBw7kdquuh5mGrYf169dvyczxuebrFd4RcRBwKfA84Hbgw8AlmXnxbNcZHx/PzZs373De1NQUExMTcy5vmNnDcKjaw9hZV9w3vXHdds7dunpZlnvDOScN5HarroeZhq2HiOgV3n2HTZ4MXJ+ZP87Me4CPAE9YTIGSpIXrG943Ao+LiH0jIoATgWsHV5YkaXd6hXdmXgVcAnwN2Npdb9MA65Ik7UbvAbfMPBs4e4C1SJJ68hOWklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklRQ7/COiAdExCUR8W8RcW1EPH6QhUmSZrd6HvO+A/hEZj43IvYB9h1QTZKkOfQK74g4ADgeOAMgM+8G7h5cWZKk3YnMnHumiEcDm4BvA48CtgBnZuadO823AdgAMDo6euzk5OQOt7Nt2zZGRkaWpvIVYg/DoWoPW2++477p0TVw610rWMwS6NPDusMOXJ5iFmjYHkvr16/fkpnjc83XN7zHga8AT8zMqyLiHcBPM/N1s11nfHw8N2/evMN5U1NTTExMzLm8YWYPw6FqD2NnXXHf9MZ12zl363xGLodPnx5uOOekZapmYYbtsRQRvcK77xuWNwE3ZeZV3elLgMcstDhJ0uL0Cu/M/BHww4h4WHfWibQhFEnSCpjPa7ZXAu/vjjT5PvCHgylJkjSX3uGdmd8A5hyHkSQNnp+wlKSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSCDG9JKsjwlqSC5hXeEbEqIr4eEZcPqiBJ0tzmu+d9JnDtIAqRJPXXO7wj4nDgJOCCwZUjSeojMrPfjBGXAG8B9gf+LDNP3sU8G4ANAKOjo8dOTk7ucPm2bdsYGRlZbM0ryh6Gw2J62HrzHUtczcKMroFb71rpKhbHHnZt3WEHLvi669ev35KZ43PNt7rPjUXEycB/ZuaWiJiYbb7M3ARsAhgfH8+JiR1nnZqaYufzqrGH4bCYHs4464qlLWaBNq7bzrlbez0Fh5Y97NoNp04s6e3tSt9hkycCz4yIG4BJ4ISIuHhgVUmSdqtXeGfmn2fm4Zk5Bjwf+ExmvmiglUmSZuVx3pJU0LwHejJzCpha8kokSb255y1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBRneklSQ4S1JBfUK74g4IiI+GxHXRsQ1EXHmoAuTJM1udc/5tgMbM/NrEbE/sCUiPp2Z3x5gbZKkWfTa887M/8jMr3XTPwOuBQ4bZGGSpNlFZs7vChFjwOeBR2TmT3e6bAOwAWB0dPTYycnJHa67bds2RkZGFlHuyrOH4bCYHrbefMcSV7Mwo2vg1rtWuorFsYddW3fYgQu+7vr167dk5vhc880rvCNiBPgc8KbM/Mju5h0fH8/NmzfvcN7U1BQTExO9lzeM7GE4LKaHsbOuWNpiFmjjuu2cu7XvyOVwsoddu+GckxZ83YjoFd69jzaJiL2BS4H3zxXckqTB6nu0SQDvBq7NzLcNtiRJ0lz67nk/ETgNOCEivtH9PX2AdUmSdqPXQE9mXgnEgGuRJPXkJywlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqSDDW5IKMrwlqaCh/Nnnlfx178X86vNiLGfPG9dt54wZy6vY8849SP/fuOctSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQUZ3pJUkOEtSQX1Du+IeGpEXBcR34uIswZZlCRp93qFd0SsAs4DngYcDbwgIo4eZGGSpNn13fP+HeB7mfn9zLwbmASeNbiyJEm7E5k590wRzwWempkv6U6fBhyXma/Yab4NwIbu5MOA63a6qbXAbYsteoXZw3Cwh+FgD0vv1zPz4LlmWt3zxmIX5/1K6mfmJmDTrDcSsTkzx3sucyjZw3Cwh+FgDyun77DJTcARM04fDtyy9OVIkvroG95XAw+NiCMjYh/g+cDHBleWJGl3eg2bZOb2iHgF8ElgFXBhZl6zgOXNOqRSiD0MB3sYDvawQnq9YSlJGi5+wlKSCjK8Jakgw1uSChr68I6IEyLiyJWuYzGq91C9/ml7Qh/2MByGoYehfcOy++6USeB24F7g7Mz83MpWNT/Ve6he/7Q9oQ97GA7D1MPQ7HlHxOERccCMs54HXJqZx9PurBdExONXprp+qvcwn/ojYlefuh0K1dcD2MOwGOYeVjy8I+LhEfFx4ErgLyNi+guvfg7s201/iPbdA48bxtCo3sNC6s8hfMlWfT2APQyLCj2sSHhHxH4zTj4auCkzx4DPAH/Xnf8T4BcRsX9m/gT4DjAKjC1jqbOq3sNu6v8sv1r/SFf/dxmS+qdVXw9gD9jDgixbeEfEQRHx3oi4GjgnIg7utlaPBL7Y7c19DLg9Ik6i3Sn7A+u6m/gu7du/7l6umndWvYfd1L9uRv0f3UX9j+xu4jsrWf+06usB7KG7CXtYhOXc8z4e2A48nfYthX8BHNDVcOiMl+EXAS8Evgr8jPYDEGTml4ETgJ8uY807q97DbPWvokb906qvB7AHe1iszFzSP1oQvBT4HO27vdd2538IeFU3fSRwTnf5Y2njSqu6y0aAH3e3cxhwLfAK4D3A+cC+S13zntZD9fr3pD7swR4G9TeIPe+TgWcCbwQeD/xtd/6ngSd00z8EvgA8LTOvpm3x1gNk5jbgKuCxmXkzcBptTOlHwGsz838GUPOe1kP1+vekPuzBHgai748x7GD6aIOIeCztpcQXgCsy8xfAbwLfz8zPRMT1wFsj4inAFuDZEbE2M2+LiO8Cd0bEg4F3Ai+KiENo3xX+X7SXJ2TmZmDzIvvc43qoXv+e1kfXiz34nF42897znnEHHQ9cSDt05snAW7pZ7gW+ExFrMvN62suLR9LGiW6hHScJ8EvaS5C9gEtpX8t4KnAssCkz711wV3P3sKrr4Um0lz2leujqyoiYoOg66Po4oHofEfGg7v8TgPcW7eGQiHhQRIzTxnYr9rBP5ef0gvQZW6Ed1/jHwD8BfwTsDfwJ8PLu8oOAbwLH0O6Ec4Cx7rKTaXfA2m56K3AgbZD/48A+M5az16DGh4D9gJfQVshG2psS1XrYH7iC9n3qAK+uVP+Mx9LpwL/SPuxQqg/+71PJv0cbE52i/V5rtcfSfsAZtGGDO4CTCvawN/Ay4DLgH4CHAGdW6mExf3PueUfEocDlwATwPtpg/im0caPtAJn538BHgVd1D4ZDgId3N/F54EnA3Zl5OfBu4BLgPNpW/p7pZeXg9lT3o4XFCcC7gKcAz6G9KXFvhR46a4D7AUdFxFrgKNpeQon6I2Jv4BrgucBbM/M53UXHzFj+UPeRmRkRBwK/D7w9Mycy8zranlmJHiLiIbQhhROB1wI3AzfS3rAr83gCXk57Tr+d9pu6p3Snf1moh4XrsXVbQ/ul+OnTZ9AG+08Hvjrj/F8DbummX0776OhB3fUvAx48Y961y72VAh4wY/o1tBV6arEeTgfeCrwOeDHwDODqKvV3y/0IcOpO5z0PuKpKH7S9vb/qpqf3xE+p0gMtoO834/SFtA3qs6r00C3zMuAPuukXA6/sHktlntOL+esz5v1z4KvdgevQXoYcl5kX0fYADwXIzFuAayLiuMw8D/ge8EHajxdfmZk3Tt9gZt7WY7lLKjNvj4gDIuK9tGGTtbSVd1REjA5zDzPu+72Af6ftvZ6QmZcBRw57/Tu5EDg7Is6NiKmIeD3wFdpvpB7S1TbsfdwG/G5EnApsiYiLaHt7v9W9IhrqHjLzl9nexJsesw/aFy1dRlsPJZ7TwL8Ap0fEh4E3AI8CvkXr4eCurmHvYeHmsZWb3sO4CDizm34f8Dfd9AOBC+i2ZLTxqEcA91/pLdROfbyM9vJwE20M/EvA62kP4KHuAfgwba/pQNrxqq+lPVhfV2wdfJL2YYgjaE+kM4EvF1oPD+3q/XvaXtwLgbfRjgN+DW0jO9Q97NTP14FTuumLKz2nu+fwhbSdsTcAZwPXdY+vUuthvn+9jzbJzIyIw4FDaVs8ujuLiLicNqa8KrstWWbek5nfysyf913GcsjM87O923w+bZzyH2kvoT7GEPcQESO0Pb53AZ+gvVN+HPAC4KCIuIwhrn8nz87MN2fmD4E30w7jeicF1kPnRtrHoVdnG1e9vDvvctobshV6ICKmn/9X0t7sA3hTu6jM4+loYCrbXvN7aDs3H6DOY2nh5rmVewbtC1r2pm3xnkq7s14IPGalt0Tz7OUI4FPAg7rTLwKOWem6dlPv/Wmvei6gvdE0AXxqxuVDXf9u+vp12ptJD6zUB+3j0Vu66QfQXkEcU6mHrtZ9aRvOU3Y6/9Rh76HLnj8FLuhOr6XtWB5ZbT0s5G9eP8YQEV8EfgO4gXZs5Bsz85u9b2CFdUcJnEjb2BxNGzo5LzPv2e0Vh1D3IYJTgMnM/NFK1zMfEXE/2ob/NOC3aYd5nZ+Z21e0sHmKiDfRng/H0HYEzs62J15KRFwHvD4zPzj9OY6VrqmviDiK9jy+m7Yu/hn462yfiNyj9Q7v7jCvs4HrgYuze8OjkohYTft+g1/Qeij30ikiVgH3VnqC7UpEvJR2mOb7Kq6HaRHxMOAHFXuY8YG7R9PeBN9e8XHV7cg8FPhSZt610vUsl6H9GTRJ0uxW/Jd0JEnzZ3hLUkGGtyQVZHhLUkGGtyQVZHhLUkGGtyQV9L8Giii4LM1vYwAAAABJRU5ErkJggg==\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": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "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 }