{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import isoweek" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "data_url = \"https://www.sentiweb.fr/datasets/incidence-PAY-7.csv\"" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
weekindicatorincinc_lowinc_upinc100inc100_lowinc100_upgeo_inseegeo_name
02022467294111004782417FRFrance
12022457389517626028639FRFrance
22022447427122316311639FRFrance
320224375863330284249513FRFrance
42022427377019505590639FRFrance
52022417417722196135639FRFrance
620224074883147282947212FRFrance
7202239720413313751306FRFrance
8202238717714193123315FRFrance
9202237717254992951315FRFrance
10202236710691781960213FRFrance
11202235715814002762204FRFrance
12202234722667883744315FRFrance
132022337734001739911026FRFrance
142022327780140861151612618FRFrance
15202231768964170962210614FRFrance
162022307903957701230814919FRFrance
172022297148511006019642221529FRFrance
182022287154711102819914231630FRFrance
192022277211911619826184322440FRFrance
202022267168541280620902251931FRFrance
212022257222461801126481342840FRFrance
222022247224581810526811342741FRFrance
232022237187721487522669282234FRFrance
242022227189161494122891292335FRFrance
252022217203101630724313312537FRFrance
262022207235851900428166362943FRFrance
272022197185931418123005282135FRFrance
282022187178511396321739272133FRFrance
292022177203141600124627312438FRFrance
.................................
16381991267176081130423912312042FRFrance
16391991257161691070021638281838FRFrance
16401991247161711007122271281739FRFrance
1641199123711947767116223211329FRFrance
1642199122715452995320951271737FRFrance
1643199121714903897520831261636FRFrance
16441991207190531274225364342345FRFrance
16451991197167391124622232291939FRFrance
16461991187213851388228888382551FRFrance
1647199117713462887718047241632FRFrance
16481991167148571006819646261834FRFrance
1649199115713975978118169251832FRFrance
1650199114712265768416846221430FRFrance
165119911379567604113093171123FRFrance
1652199112710864733114397191325FRFrance
16531991117155741118419964271935FRFrance
16541991107166431137221914292038FRFrance
1655199109713741878018702241533FRFrance
1656199108713289881317765231531FRFrance
1657199107712337807716597221529FRFrance
1658199106710877701314741191226FRFrance
1659199105710442654414340181125FRFrance
16601991047791345631126314820FRFrance
16611991037153871048420290271836FRFrance
16621991027162771104621508292038FRFrance
16631991017155651027120859271836FRFrance
16641990527193751329525455342345FRFrance
16651990517190801380724353342543FRFrance
1666199050711079666015498201228FRFrance
16671990497114302610205FRFrance
\n", "

1668 rows × 10 columns

\n", "
" ], "text/plain": [ " week indicator inc inc_low inc_up inc100 inc100_low \\\n", "0 202246 7 2941 1100 4782 4 1 \n", "1 202245 7 3895 1762 6028 6 3 \n", "2 202244 7 4271 2231 6311 6 3 \n", "3 202243 7 5863 3302 8424 9 5 \n", "4 202242 7 3770 1950 5590 6 3 \n", "5 202241 7 4177 2219 6135 6 3 \n", "6 202240 7 4883 1472 8294 7 2 \n", "7 202239 7 2041 331 3751 3 0 \n", "8 202238 7 1771 419 3123 3 1 \n", "9 202237 7 1725 499 2951 3 1 \n", "10 202236 7 1069 178 1960 2 1 \n", "11 202235 7 1581 400 2762 2 0 \n", "12 202234 7 2266 788 3744 3 1 \n", "13 202233 7 7340 0 17399 11 0 \n", "14 202232 7 7801 4086 11516 12 6 \n", "15 202231 7 6896 4170 9622 10 6 \n", "16 202230 7 9039 5770 12308 14 9 \n", "17 202229 7 14851 10060 19642 22 15 \n", "18 202228 7 15471 11028 19914 23 16 \n", "19 202227 7 21191 16198 26184 32 24 \n", "20 202226 7 16854 12806 20902 25 19 \n", "21 202225 7 22246 18011 26481 34 28 \n", "22 202224 7 22458 18105 26811 34 27 \n", "23 202223 7 18772 14875 22669 28 22 \n", "24 202222 7 18916 14941 22891 29 23 \n", "25 202221 7 20310 16307 24313 31 25 \n", "26 202220 7 23585 19004 28166 36 29 \n", "27 202219 7 18593 14181 23005 28 21 \n", "28 202218 7 17851 13963 21739 27 21 \n", "29 202217 7 20314 16001 24627 31 24 \n", "... ... ... ... ... ... ... ... \n", "1638 199126 7 17608 11304 23912 31 20 \n", "1639 199125 7 16169 10700 21638 28 18 \n", "1640 199124 7 16171 10071 22271 28 17 \n", "1641 199123 7 11947 7671 16223 21 13 \n", "1642 199122 7 15452 9953 20951 27 17 \n", "1643 199121 7 14903 8975 20831 26 16 \n", "1644 199120 7 19053 12742 25364 34 23 \n", "1645 199119 7 16739 11246 22232 29 19 \n", "1646 199118 7 21385 13882 28888 38 25 \n", "1647 199117 7 13462 8877 18047 24 16 \n", "1648 199116 7 14857 10068 19646 26 18 \n", "1649 199115 7 13975 9781 18169 25 18 \n", "1650 199114 7 12265 7684 16846 22 14 \n", "1651 199113 7 9567 6041 13093 17 11 \n", "1652 199112 7 10864 7331 14397 19 13 \n", "1653 199111 7 15574 11184 19964 27 19 \n", "1654 199110 7 16643 11372 21914 29 20 \n", "1655 199109 7 13741 8780 18702 24 15 \n", "1656 199108 7 13289 8813 17765 23 15 \n", "1657 199107 7 12337 8077 16597 22 15 \n", "1658 199106 7 10877 7013 14741 19 12 \n", "1659 199105 7 10442 6544 14340 18 11 \n", "1660 199104 7 7913 4563 11263 14 8 \n", "1661 199103 7 15387 10484 20290 27 18 \n", "1662 199102 7 16277 11046 21508 29 20 \n", "1663 199101 7 15565 10271 20859 27 18 \n", "1664 199052 7 19375 13295 25455 34 23 \n", "1665 199051 7 19080 13807 24353 34 25 \n", "1666 199050 7 11079 6660 15498 20 12 \n", "1667 199049 7 1143 0 2610 2 0 \n", "\n", " inc100_up geo_insee geo_name \n", "0 7 FR France \n", "1 9 FR France \n", "2 9 FR France \n", "3 13 FR France \n", "4 9 FR France \n", "5 9 FR France \n", "6 12 FR France \n", "7 6 FR France \n", "8 5 FR France \n", "9 5 FR France \n", "10 3 FR France \n", "11 4 FR France \n", "12 5 FR France \n", "13 26 FR France \n", "14 18 FR France \n", "15 14 FR France \n", "16 19 FR France \n", "17 29 FR France \n", "18 30 FR France \n", "19 40 FR France \n", "20 31 FR France \n", "21 40 FR France \n", "22 41 FR France \n", "23 34 FR France \n", "24 35 FR France \n", "25 37 FR France \n", "26 43 FR France \n", "27 35 FR France \n", "28 33 FR France \n", "29 38 FR France \n", "... ... ... ... \n", "1638 42 FR France \n", "1639 38 FR France \n", "1640 39 FR France \n", "1641 29 FR France \n", "1642 37 FR France \n", "1643 36 FR France \n", "1644 45 FR France \n", "1645 39 FR France \n", "1646 51 FR France \n", "1647 32 FR France \n", "1648 34 FR France \n", "1649 32 FR France \n", "1650 30 FR France \n", "1651 23 FR France \n", "1652 25 FR France \n", "1653 35 FR France \n", "1654 38 FR France \n", "1655 33 FR France \n", "1656 31 FR France \n", "1657 29 FR France \n", "1658 26 FR France \n", "1659 25 FR France \n", "1660 20 FR France \n", "1661 36 FR France \n", "1662 38 FR France \n", "1663 36 FR France \n", "1664 45 FR France \n", "1665 43 FR France \n", "1666 28 FR France \n", "1667 5 FR France \n", "\n", "[1668 rows x 10 columns]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "raw_data = pd.read_csv(data_url, skiprows=1)\n", "raw_data" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "def convert_week(year_and_week_int):\n", " year_and_week_str = str(year_and_week_int)\n", " year = int(year_and_week_str[:4])\n", " week = int(year_and_week_str[4:])\n", " w = isoweek.Week(year, week)\n", " return pd.Period(w.day(0), 'W')\n", "\n", "data['period'] = [convert_week(yw) for yw in data['week']]" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [], "source": [ "sorted_data = data.set_index('period').sort_index()" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "periods = sorted_data.index\n", "for p1, p2 in zip(periods[:-1], periods[1:]):\n", " delta = p2.to_timestamp() - p1.end_time\n", " if delta > pd.Timedelta('1s'):\n", " print(p1, p2)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXm4HUWZ9+89527ZdyAkQICENUIgAQKIyCIEnRH0g2+CI+CIooi4fp+iMy6MMuqooOgniqLsAooKwyabCEIIXJKQECBkIWTf9+Su59T3R3d1V1fX1uf0Offm3vo9z31un+rqquru6nrr3YkxBg8PDw8PDxGFnh6Ah4eHh0fvgycOHh4eHh4peOLg4eHh4ZGCJw4eHh4eHil44uDh4eHhkYInDh4eHh4eKXji4OHh4eGRgicOHh4eHh4peOLg4eHh4ZFCQ08PoFKMHj2aTZgwoaeH4eHh4bFX4ZVXXtnEGBtjq7fXEocJEyagtbW1p4fh4eHhsVeBiN5xqefFSh4eHh4eKXji4OHh4eGRgicOHh4eHh4peOLg4eHh4ZGCJw4eHh4eHil44uDh4eHhkYInDh4eHh4eKXji4LFXgzGG+19ZhbbOUk8PxcOjT8ETB4+9GrOWbcaX//AqvvPw6z09FA+PPgVPHDz2auxq7wYAbNjR0cMj8fDoW/DEwcPDw8MjBU8cPDw8PDxSsBIHImohopeI6FUiWkhE14bl3yai1UQ0L/x7v3DN14hoCREtIqJzhfKpRLQgPHcjEVFY3kxE94bls4loQv636tG3wXp6AB4efQounEMHgDMZY8cCmAJgBhFND8/dwBibEv49AgBEdBSAmQCOBjADwC+IqBjWvwnAFQAmhX8zwvLLAWxljE0EcAOAH1R/ax79AeH+wsPDI2dYiQMLsCv82Rj+mbZp5wO4hzHWwRh7G8ASACcS0VgAQxljsxhjDMDtAC4QrrktPP4jgLPIf/UeHh4ePQYnnQMRFYloHoANAJ5gjM0OT32WiOYT0W+JaERYNg7ASuHyVWHZuPBYLk9cwxjrBrAdwCjFOK4golYiat24caPTDXp4eHh4ZIcTcWCMlRhjUwCMR8AFTEYgIjoUgahpLYAfh9VVO35mKDddI4/jZsbYNMbYtDFjrImMPDy0mL9qG5Zt3GWv6OHRT5HJWokxtg3AMwBmMMbWh0SjDODXAE4Mq60CcIBw2XgAa8Ly8YryxDVE1ABgGIAtme7EwyMDPvjz53Hmj//e08Pw8Oi1cLFWGkNEw8PjAQDOBvBmqEPg+BCA18LjBwHMDC2QDkageH6JMbYWwE4imh7qEy4F8IBwzWXh8YUAng71Eh4eTvCzxcMjX7jkkB4L4LbQ4qgA4D7G2ENEdAcRTUEg/lkO4FMAwBhbSET3AXgdQDeAqxhjPPDNlQBuBTAAwKPhHwDcAuAOIlqCgGOYmcO9efQDeKsFD4/awEocGGPzARynKL/EcM11AK5TlLcCmKwobwdwkW0sHh4eHh71gfeQ9ugz2LCjHWf++Bms3LKnp4fi4bHXwxMHjz6D++esxrKNu3Hn7Hd6eigeHns9PHHw8PDw8EjBEwcPDw8PjxQ8cfDwqBEWrtmOa+6fj3LZ29l67H3wxMGjz4D1ssisH7/1Zdzz8kps2OkTEXnsffDEwcPDw8MjBU8cPPoEehfP4OGx98MTB4+9GqrA7uT9pj08qoYnDh4eHh4eKXji4OHh4eGRgicOvRxLNuzC6m1tPT2MXg8fxNfDI1+4RGX16EGcfX2Qc2D59z/QwyPx8PDoT/Ccg8deDZ9p3MOjNvDEwaPPoLdJlnrbeDw8ssATBw8Pj70OD8xbjXc27+7pYfRpeOLg4VEjeJFXbbBtTyc+f888fPzWl3t6KH0anjh49Dn0lkXZi5Vqg65S8GCXbvScQy1hJQ5E1EJELxHRq0S0kIiuDctHEtETRLQ4/D9CuOZrRLSEiBYR0blC+VQiWhCeu5Eo+IyJqJmI7g3LZxPRhPxv1aO/oL8uyk+9sR4TrnkYKzb37Ux4vS3AYl+FC+fQAeBMxtixAKYAmEFE0wFcA+ApxtgkAE+Fv0FERwGYCeBoADMA/IKIimFbNwG4AsCk8G9GWH45gK2MsYkAbgDwgxzuLTcwxnD37BXY0d7V00Px2ItQbw7mT3NWAwBeXbWtvh3XG5421AVW4sAC7Ap/NoZ/DMD5AG4Ly28DcEF4fD6AexhjHYyxtwEsAXAiEY0FMJQxNosFHku3S9fwtv4I4CzOVfQGzFmxDV//8wJ8/U8LenooHg7oLTOnpziY3nL/Hns3nHQORFQkonkANgB4gjE2G8C+jLG1ABD+3yesPg7ASuHyVWHZuPBYLk9cwxjrBrAdwKhKbqgWaO8qAQC27O7s4ZF4eOjRX8Qt/eMuex5OxIExVmKMTQEwHgEXMNlQXbVvYYZy0zXJhomuIKJWImrduHGjbdgeHj2Keu/gOafS16PS9ledUr2RyVqJMbYNwDMIdAXrQ1ERwv8bwmqrABwgXDYewJqwfLyiPHENETUAGAZgi6L/mxlj0xhj08aMGZNl6B69BN2lMtb0k1hRtVjEnlm0ARs1meUi4tC3aUO/4ZB6Gi7WSmOIaHh4PADA2QDeBPAggMvCapcBeCA8fhDAzNAC6WAEiueXQtHTTiKaHuoTLpWu4W1dCOBp1gsjqb25bieuf+Ktnh7GXo3vPfomTvn+09oFTsbPn17s9Mx73WSpET72u5dx0S9fMNbp47TBcw51ggvnMBbA34hoPoCXEegcHgLwfQDvI6LFAN4X/gZjbCGA+wC8DuAxAFcxxkphW1cC+A0CJfVSAI+G5bcAGEVESwB8CaHlU2/Dlt2duPGpxZEOohos37Qbm3f1v9zCzywKGMxte9z0Nz96/C3c+NRi7fm9YaGodqfb2V3Gi8s2R5Fnl2tMVf2O2iNPWKOyMsbmAzhOUb4ZwFmaa64DcJ2ivBVASl/BGGsHcJHDeHsFukpltDQW7RUNeO+PnkFTQwFvffe8nEa1d4AboeW9jD2zaCNOmDAy6CPntnsa//XIG7j1heV45HOnOdXv+2Ilj3rAe0hXgM7ucq9qZ28CX7fy2vGr2ulti0e19/rW+p0AgM27zZzm3sBF5YFeKHHuk/DEoQJ0lvrfop4X8t7V7g3LRLVj5M+sbGkoPt23WQdPG+oDTxwqQH/c8ecFbmaZl3xctYvsa0sjf2Zlx1WRE5NHFqzFrc+/XathefRx+ExwFcATh8rBF67cxEr5NFNTVCsGiZ+Zvp1NuzpSSv7P3DUHAPCxUw+uqn+P/glPHCpAhycOVaOWOofehmrHyJX4ZcO0m/bdJ+P61XXX67E3vPO+AC9WqgCu7H212N3RXZd+egL5mV0G7bQ0FvqsorIQrvalPnp/WeFNdusDTxx6MR59bV1PDyF3RKasOXMOA5v6LhMcW3i56hz6Nu/gaWR94IlDBfCTs3LkvWztDa+i2vlS4GIlx3b6NmnYO955X4AnDhWgXpOzL4pJcjdlVTyivrZxjk1Zs1kreXhUA08cKkBfXLTrhfytlYKGugTfk7zavuL2Vrz/p88lyh57bR3aOrOFT6leRp6Nc+jr8N9ffeCJQwXwU7N65OfnEPzf2d6NF5elAvlWhcdfX4/X1+6Ifs9ftQ2fvvMVfOvB1zK147KWPTBvNV5enh7/S29vwZNvrAcA7HE0UOjrnIP//uqDvqvFqyH8xqVy5J1rQBS1/GPJpqCPGi2OO9uDxXnV1vxDjn/+nnkAgOXf/0Ci/NHX1kbH1/hMhAD891cveM6hIvjZWSlUYqVZSzfj4flr1RdY0F1yexfb27qinfnqKvNJZF2cqpktlRDTvp7sx39/9YEnDhXA71wqR2SWKZRd/OsXcdXdcypqr+QoiL/i9lZc9MtZaOss4YwfPVNRX5WizBi+/eBCLNmwy15ZQiVcUHtXCWWvoPCoEp44VAD/2dnBGMP1T7yFlVs0uQdyorBdJrdhAYvCyKZtXaWqw59kXbCXb9qNW19Yjk/e3pq5r0IFxOHKu+bgv/+6KPuFewn85qw+8MShAvjJaceyTbtx41OLccUdryRP5JzPwZVzaCoGUz2PuFiZxUpV3GylDm2//PvSyjvt5RAf55/nrsKrK7f12Fj6MrxCugJ4Uzo7+DPq7E6afeadz8FV59DUEBCHju7Ks/hVKsmvSufQ19UHFUCcO1+891UAaUW+R/XwnEMF8KTBDt3in/di58w5hMShLYcUr67gI5uzYiuAyohL31cuZ4ePrVQfWIkDER1ARH8jojeIaCERfT4s/zYRrSaieeHf+4VrvkZES4hoERGdK5RPJaIF4bkbKeSZiaiZiO4Ny2cT0YT8b9WjJyCLReJf9dU5cLFSe1f9I+re9Ewo4qlgne8tnMP8Vds8x9zP4MI5dAP4MmPsSADTAVxFREeF525gjE0J/x4BgPDcTABHA5gB4BdExBMu3wTgCgCTwr8ZYfnlALYyxiYCuAHAD6q/tdrBfyN26B5R3oH3So5ipeYGThzqxznI67punVc5v3FUopDOG48uWIsP/vx53D9ndU8PBYD//uoFK3FgjK1ljM0Jj3cCeAPAOMMl5wO4hzHWwRh7G8ASACcS0VgAQxljs1iwBbkdwAXCNbeFx38EcBbJW85eBM/WVo+8nmC3QqykEsX0pFiJQzelTQrV3iBWenvzbgDA4g07lec37uzAe3/4N7y9aXddxuOJQ32QSecQinuOAzA7LPosEc0not8S0YiwbByAlcJlq8KyceGxXJ64hjHWDWA7gFFZxlZX+MlphW5Jy1sh7aJzeHzhOqzZ1g4A6MiBOFS6OdA+E8M+qDdskSICpbntRxasxfLNe/C7OqUk9Zuz+sCZOBDRYAD3A/gCY2wHAhHRoQCmAFgL4Me8quJyZig3XSOP4QoiaiWi1o0bN7oOPXf0t6n5xtod+J9X12S6Ri9WCs/nRB1cIpVecccrkVd0VTqHjAu1a3VTPR3h+Mvc+ol4bASKBz1sKNTHvsVzDvWB09skokYEhOEuxtifAIAxtp4xVmKMlQH8GsCJYfVVAA4QLh8PYE1YPl5RnriGiBoADAOQEsQyxm5mjE1jjE0bM2aM2x3WAP1tcp730+dw9e/nVnRtWu6e71Y466u4f84qe6WcOkuLldT1TIuv7tQX7p2XbTBVQOXVLoKL9uau3IoXl22uy5g8ag8XayUCcAuANxhj1wvlY4VqHwLAQ1U+CGBmaIF0MALF80uMsbUAdhLR9LDNSwE8IFxzWXh8IYCnWS82jejLbO0ds5bjS/dVv/Bo355ZQpFfPxo8t3hTTj3bsXFnR+K3jjCaOYccB1QhonwSGhFed8g5zF2xDTNvfrGmY7n35RX48E0v1LQPjwAuTnCnArgEwAIi4qvG1wFcTERTEHznywF8CgAYYwuJ6D4AryOwdLqKMcYFvVcCuBXAAACPhn9AQHzuIKIlCDiGmdXdVnYsXLMdLY1FHDpmsLVu7yVb1eMbDywEAFz/v6fk0p68uOWtc6groa5yoa5koS/0AurAiZruSXc5Wozlga/e7yPT1gtW4sAY+wfUn8UjhmuuA3CdorwVwGRFeTuAi2xjqSU+cOM/ALh5WvZh2pA7dESgL3NfWWFUSNdxHDrYEjR1O/qaAMD2PV3Y0d6FA0YOzGFkHrWE95CuAL1Y4pUbbPf41BvrsXlXh7GOCmQTYGdFD7yKvF+/iTko9AZHhxA6gu4awgQAzr7h7zjtv/+W15A8aghPHCpA3XJI16kfFUwWom2dJVx+Wys+9ruXre2kxUr5Bt7bm8i0jkPo7XuNgsVx0TW3NZDWw3j0XnjioMH8Vdsw4ZqHa259YYq735N7RtMHXwrPLTc4Pel2mZFyM6cVUcXh1FpMn3f7Ji6tF6gccjc/9nDHBf/veXz9zz2jZ/HEQYPnlwRE4W+LNqRP5viNdJb08lreTU9IFkyLN18kdDW2t3XhlucChyidhU5uCum9SKyke42m5nqDh7RtBL04mMFej3krt+Hu2St6pG8fsrsCVKNM3bizA4vWxWEIukpltDQWlXX5Ilyvj+/bDy6MjjPoGFP45gOv4YF5aqe53DmHXFpxQ7ULte41mh5FPV79y8u3YHBzA44cO9RYz/MN/QueOFjQ1lnCJknxWs269i83z8KyjbE4xpR8hkuc6rUvu/WF5ULf6ptcs60NAzTEjGN7W5f2XO46h71oxdISB9M1NRlJEhf9chYAvaWeLVii5xv6JrxYyYLbZ72Dz9+TdAqrZkESCQNgthEvR5xD5f1VChVxaOss4ZTvP42v3j+/4nbzll/vTSaxOs5j+IBG7TW9ws8hclys7llv36PfNHC88s4WzPjJs2jrrF+ARA81PHHQwPRNVvqJqGSHbpxD/RcIlViJZ1F74o31zu3w3M0y8qAN29u68MTr7mPpaejmFFfwHzJ6kPM1QP0UxJxzkG0n1m1vDwwqHKfn8d99IjrWjf27D7+BN9ftxOtrt1c0Vo/84IlDBaj0o1RZHRgV0qzOciUBKs4hEglVsSbpFppKcOy1j2PV1jbn+qMGNVXfaQ1w41OLg4OM71kVrrwWUHm1r9i8B9O/9xT+39+WGDcvn7nrFRz270EgBDGCrm7sLQ2ByLKjB5Iy9Wa8uW5H3fv0xKEC5PlJGjmH8APqCcGCUueQ40B6wiwyL4eySkeu650TuKwipCzOZzI6ukvYurvTqW48rLi/NduDMdtiVT2yYJ1yA6QLtd7cGCxJH/nNbDy/pH5xsHo7PnFba9379MShAuS5rpmsdiKxUo9wDorCHO6b30qdNr0JFKt8kFW/B0sDWWmXa4pUFf7tdy/juO88Ya+ImHiJU1UcaiXPRcc58HSuAHD9E29lb7iPYosjIc8TnjhUhPxWNhOhiRTSPaFzUAwsD/NT1S40b+ieVrWMQ62ZnWLGfAiuKVJVeGGpu3NnlAM7R3RrxKkNxfgluSRy6i/oCcMETxw0ML2KL9/3am79mCxAWI9yDjUiDlFbVTeVGT0dp8jWe7GQ3iGafFyq4Rwqger1M7CKti7agIxCebXz7bHX1uLpN/cegwUTemLmeuJQAXbnaGZnWiTzchSrBKpxuS7o8kRW6Rd64taKVRKHaom06npR5/Ta6h04/jtPYN329qjMpJupRufA8dB89wx/4nys1jHTZeTVzv9P3zkHH781H1l9V6mMdzbXJ0d2b0G/Jw7XP74IF/Zg8hDTx19vJ7hE3wpKUKkSOSGrjqyV6k8dqtU5VAtV76rESht2titqppEHcfjs3ekMf3NXbFVax+h6q+Sx6t6/WGww5Ks7vvPQ6zj9h884v5vc0QNTt98ThxufXoLWd7b2WP8unENPxK5Rfbxyya6O7sxt5R2xOwtyEytVGltJ8R7/unBd5ms4aiVW+tAvXsCMnzyXKtd7SGd/ri55PkxBKbPg2w8uNFoFuoDraFwc+WqBne1u31qe6PfEoedhD3BXD9LwwtKk2aBarJQuVH108oLGEufCsj7OOSifi/R76+5OpYe862Jb70eYpze67v0nOIecbvDWF5ZnEp95BPDEQYN6rSNmziE8sIzloflrMOGah/H+n6Z3exzdpTIWa7yVAeAjv56d+L1DER9JNVaXLGBJolK9I12lyI1zcGjGxXP7k7er5eHuc8/+EBtyVMJv2JHOxcBYZd+KbuTiHMtT9JiXAUR/sp/yxKGH4WbKagaXG7++Vu9F+cPHF+F9NzyLZRt3OY1LJWpT6yHsbSV1DmFbPcE55DXbLUNfs60NixQye3kRdRFnmt69y4JXrRJexD8Ep7RqN0/61x+fWLllT3WdCOj5CFV7H6yfCxEdQER/I6I3iGghEX0+LB9JRE8Q0eLw/wjhmq8R0RIiWkRE5wrlU4loQXjuRgrlD0TUTET3huWziWhC/rfaO+HiBJcH5oQL0aZdbs4033nodezq6MaCVXGMG9VQXVh/pQlkj1gruVOHbXvSz8l1gbnkltm48ekliuvzXaKcCHOuPQbo6C5FnCVDZQuvi0K6q8ScvbhtyEsSsDdFAa4WLl9LN4AvM8aOBDAdwFVEdBSAawA8xRibBOCp8DfCczMBHA1gBoBfEBGP8XwTgCsATAr/ZoTllwPYyhibCOAGAD/I4d72Cpgmmy2pTq2xcWcHNu+ORQlK3wcFBZO/Q5VCuqacg2YlaM7AOlQTeXbpRo3JY87bV6dnmPNjZozhI7+ejcszhnOQdQyuw9rRno8CuGoz5FxGkQ33vbyyB3qNYf1aGGNrGWNzwuOdAN4AMA7A+QBuC6vdBuCC8Ph8APcwxjoYY28DWALgRCIaC2AoY2wWC2bK7dI1vK0/AjiLTGYafQhmU9bwXA9RB5K6Vi1GLl6sSfv44L8yOkeNt2U8bo8LTNYhLy3fUvVYjOlhHWe+G23I95kyBrwii8McBrx8c1JEpHvXcliNrpzsWfPi2uq5Kon5VXoCmaSwobjnOACzAezLGFsLBAQEwD5htXEARJK3KiwbFx7L5YlrGGPdALYDGKXo/woiaiWi1o0bN2YZembUK2SF6dMt9yxtSHWuGoeTWEk4jiO71t8xrjED5zC4Of88WOKMes4QVM517rlwDnk/U7k5xtw8pH/8+CLpOnU92dGsw9EE1bax2BvFSu1dPZvTwvlrIaLBAO4H8AXGmCl+rOo16ESTJnuc1GtgjN3MGJvGGJs2ZswY25D3Cph1DqFYqQcFneLOUzUOlbGS/CGycvqcOhRDbZHF1r0WO0SxzbZON7t1cz6HKgeUE1ye1UPz1yZ+68Yum/aakmGJqHUcpp6QY7TtDcSBiBoREIa7GGN/CovXh6IihP83hOWrABwgXD4ewJqwfLyiPHENETUAGAaget69B/DO5t2YcM3DeMoxIY5Z5xD+z2FclUD+IFTfnwvnoBIrqdqqNRFUmd3q+qw15+i6IzbBJjL61OmH5D53VM+rIic4x5G5EvS+GKMvjzlSDVyslQjALQDeYIxdL5x6EMBl4fFlAB4QymeGFkgHI1A8vxSKnnYS0fSwzUula3hbFwJ4mtVwpfjnn/0Dv3luWU3anrtiGwDggXluTjdGzqHMOYfqx6WD7TGLp1XhGly8WJX6BUVpT3AO9dx9i4uo647YBNujLxDlTnBVrc1dmT3CgOtibtI5fHT6gUJ79XmRyzbuws6clOQ29GRsNcCNczgVwCUAziSieeHf+wF8H8D7iGgxgPeFv8EYWwjgPgCvA3gMwFWMMc4fXQngNwiU1EsBPBqW3wJgFBEtAfAlhJZPtUCpzLBg9XZ89+E3atI+X/RczctNr5/vyqtVKs5auhkvL09/wFt3d+Lgrz2C2wyKL3F+/uzpxekxZlVIQ58JzvQxrN9RfUybEw4emSrT9agSI4h173t5JSZc8zA27Uo7hukgtmla9FwXBdvCX6D8Ca6qy2cWZdf/aTk26bmbOIcGwTTZNg+rXWf5vL3yrjm4+NcvVteYI3o6ZLlV68YY+wf0llxnaa65DsB1ivJWAJMV5e0ALrKNJQ9sVdivq6BaHK5876HG2PbtXSW8uW5neL1rCAT9BCjlwDm0dZa0k3l5qPz709zVyvMESiwuT7+5IVVHLVaSwmcwxamMCumlG3Zh36Et2vOjBzdbF+qDR6VzNAfPP/2ubK/vhieDRDRvb9qN0YObzZUVbZqIg+v7tq0dRaIaKKQrM0mVoRu7/NxNqVBF728bQc3L6gkIoufWAzKnzhira5y1fuchXU0wr/dMSivBy2UWJS750n3z8Ku/B+Iq11domtP8wzCN2LZ77OhOKrXEubUnDD0+qKkIFeR5qJKBOomVlH4O1ssS+MhvZmPCNQ9r7/emjx5vbcPGDSTqKt6g2PXaMKx2lkVHbNO0I3blFFdttXgQ12AhyY/YuOl6TLtn0fvbFsVFN+7HF67DWT9+Rpt8KBqX9Ci3K8LL5A1ZR1ZvRqLfEQfb8zUttqMHN+GI/YbgnKP2jco+f+88TAwTqL+4LNahu1J40wuPMn0Z6tgmjGnnxaOqDmzSM5A24uOmkHZr12XhWbNdLV4izXGyjnnB37yrA+/61l8xb+U25fWqRTtL2Owk52CyUhOuMbT3+XvS4b4T/TmOqxpUSizcuSM34mCbh7p2vnr/fCzduDvzYp+X57YJ8rdbb6vF/kccLM/XdL5QIBSIEh/v/7yqVjy7OzLpO+STw/SB2OSSpsWL716bG/TTwDYdXeSi4qLKieZts97BowvWJjKfucjaqw29bBrbi8u2YGdHN25+dqmjcbVb4EEVTPeRtyIyz0Ulr6ZcxUquxEGuJ3N0OuLB02+6hpXhaDR8M3lBHrLnHGoMG8tuOlsgApH+Y1OJT2wwcg7hwmMa01uGSKuARbbNDzSDdVkIXPwcNgsfHj/19qbduPKuOYnIpC5zX3zGQ1pijseJGKvESrpFSjjebNBldHYzXHzzi/ibQh9jgomoJOZRFaIhk09JpcjL49q1HdPmQ8yrLIs3125Lcph6YhS0ce5PnnUaj66/eiBvb3cb+h9xsDxf006lQKF5oK7tykaU+HXDE2/hU3cEC2akczCM6dcWk9xOB5m4bvlhYFU9L47/eiS2DJPXOlFu7rLDFb/JpoTHs30RzbLMiosyj56qGt3O9i7MWrYZn7lrTqY2XfxbguPKF4TIG73C60/776fx62eT80seTqVt626rIE2QSsVK8m/dYl5p0FqTuLZWqLdla/4xAno57DoH/bkCEQrktiC6JiqR59hPn4rNRUsWhfSCVdut/hTVWGkE3WfntORcBuLuTw4hLS4Gbk8sruUSDsO2C+en3960OxXigWNAYzFRV0SWRULsvR7ReJOJlZL3fvIhozBr2Wbj9Su3tOG6R5Im3/LQXtXoZ2zQ3X9jUVZI69tIipWS52SiqutPJkY6yHOnVKMsfCbUmzj0Q87BttgZOIcCgSSdQ7Lt+NjVa/TvBhtxmynr2u1t1varkdEzlg/n0CTIZxulsNkJ4uAwVLE7cXEQv12RwNrAx3/Gj57Bsk2BaS9Jb2+AxpoLyEZ8meZYNyYguSj96KJjnfsCBKthxbmGYoVb5pygmzYpObuBUorfslxPvkrXTG9WWi1tAAAgAElEQVTmHEYPbkr8rrdTXD8kDpWfD8RKegJTycu7t1UfljcPJ5gtVVhVuOkA7HVE8Y+cjU1c1F1kqlrZsaZ+kmCnUWJM+T7FtKktDUXt+Dq6MhAHoR/Tc9PNo7OP3EdZ3qqJEmvSOVS60NTaYkZu3cSBi3NB/lZciUylOp0sVmqVYtiARnzgmLHR73oLsvodcagGReLWSrptT3yYh4m5uDvJ8lHe9MzSKLplVisMEYzZl2uXcY0c1KQ9lyAODrdYjVJO9U66S0y5CxSfm2mRzRI5c9ueLiFHh0npoC4ePrAJt/7bCalylff7f55/tHbhe/atjZgvJHHKgqxGAzroviG5XFz0ZV8EsWq6PUexkuMKKEeLrccuvqvEEhsrzznUGNVwDhQSB5cdfR5MeylBHNTjkbF1dyd+8Nib+MivZ+ORBWvxi7+lM5K5gjGH5+XQzkmHiNHXk1eIeROycipZrXqUxKFcTouGMphTZomcuWD1dtzDE7gYOQdhKNJYVDJyFaG59OQJ2vOX/vYlY74KE1zWJ5cNibYdg+6A++XEVQWxktSg3L7NlNWEL9/3auQwylEPsVJXqZzQwXidQ41hN2UNzqsWhAIFsloVcfjCPXOxSwjDXC3nsHlXR8Lc0XVe8I+grauEz9w1J5KjVwIGh+flsts3iFMS5ohV7Dj1FldmlMoMXd3JWht3JE1XeZeqtrKGVX5u8UbruMwWc+k7LVomW66LikNb33/0TWudSjgH2UM/KVaS20n+drWOSl7DsL2tC/fPWZU6V4+4RwFxiJdo7wRXY9ie75rQPlq1MygWCM0NBWUYib/MW1ORQlqHqd99Eu9sNpt5mnpw6d022dq7Sli4xhxHxinhTKLP5LlCVrGSyDnYqyegeifdJZYy95UzvZnyatjESv/r+PGJ33wxMj171+fFIVuAcdg2KKdOTOXTssJFrGeq8x8fODKso7s2CVG0n5YcGTgH6ImMCPEZyWKr++esxrHXPq68rh46h85umTjUvMsE+h9xEI4nXPNw6vzvX1oBQAhdIYCI0NxQdFJC5qFzWCuEinCdF7YJ9PvZK1JlOpHMN/7ymjHQoEt/QBA0j0P+iJOmrNlmf4IYax35mLFOd5lZLY5Mo2rrNBOHy045KPE7Jg76MYvPSCYiqnel2/3GWffU/QxsasBtHz9RfVLC4fsOMbaVgKHOofsMNrZjUiSbgv5ZFdIOYqVzbng20d+spXpT30o947Ogq8TQ1FDAp04/BIDXOdQc1lhB4eQYpYi0WaAgD7EczE4FV+JwoiKMtAqu84J/QLr+/zR3NZZu3KU+KeFVB6WlC6v7K8GRSq6dIA4O98g/kNXb2rBOCONdKadWUukcJGQVK4nPXh4X3+SXGTBIE9PK9ExVTILOHDPO161ur7FIzulQOXdSJW2IRGCuFn+lBKEM/n/urEl46Op3Y5Rg6GDTOfzkSbt587JNu/GHV2LrweEDG7V16+EDx3UO44cPAOCtlXocfEegWlxNYqU0XIMruVVbvMEcJiMLnn0rv/zbWSes/BGLC3wWsdJ8R+crW5NdJQfOwUAd2hRcpGhholMoM8S5l1P27EKTKWKqoARasRJvQ/MQGgoFZzt/7hdRrdw7JlhqpBTJwir85zC0/LjhLZg8bhguOXkCzp+yPwCFjsHVpVLqcOueOADfsAF64lBr+X+5HFjRNRYL0UPznEON4aKgBNQfVIGLlXIM/uY6iZ+XEtIzxvBUxng+8bUVXaZE1glrqu4mz1bDhVMTRTJ8AS+V1aasLn0CQKeCi2xpjJ3mUsShIIh6CHj2/56Bx794Oo7Yb4iyPxedg04saHsmDUVytvOPiZodZs7HLOqSr+W/y2WGH/418GDfFi7gxQLhQ8eNA2AXK2nHajhnJg5u7VcKrgdrLAoE3OscagvbS+Xhq1ULFVEwId0sFVx3Lk7V0NyQ9NL9Q+uqSD+i7rY+HrBZPxKjItahrd89/3a2DgUQ4mB9370gyDnVXXbwAo/CmKQrqkTPLY0C56ARK/HxHDhqIEYOasLtl8eyf05wd3d045EFa5P3oLJWsmz/dbfXmIFz4PXcfFH04N3pNhVyKf/URPGSeL/8uNJdtYmoDDR4xtc6CB7nZpuKBWP2xFqi3xEH26I9dliQbUw116xOcAJc9VViS6aEI01SiODV29KhM8TQ26YNoTx6XdhxF7h+lK+t3h7Wr64tHkvKPSR68veBIwcCiIlEqVx29uVQ1VMpJsV6spNVJLtnyaxeIwbGoiVOQL/xl9eioH9Re4ob10XmjRXS6hssFqni2EImGJ9npHNx0znw3+IiLmaA4+NPhc9wXEhFi0AZpiZqrY/m+T4aixQT5t4WlZWIfktEG4joNaHs20S0Wsopzc99jYiWENEiIjpXKJ9KRAvCczdSONuIqJmI7g3LZxPRhHxvMYlq2E0eeM+lDdfAeyJM4ipxNwrYTQE37nTPbVwNXO+Sh7Q2EYBqWHU3sRJw28dPxK8umYrBIXFgzP7R6cZ16JhBypSR+w2L05nKnAMJ4hlxzEl/j+D/WkViI3mn/9D8Nfjd88uV47PJ95scAhdGbYG3Vd0CVYwW82T54vU7UQq5uA8fPw5vfmdGWE8v5gXi5yZ/by7jnLsi7VmeaMPR3LgW4JxDY0Mheo+9kXO4FcAMRfkNjLEp4d8jAEBERwGYCeDo8JpfEBHnzW4CcAWASeEfb/NyAFsZYxMB3ADgBxXeixNsz5exIJn9N/7yWupcoUAoFGLOwRQU7KH5yd24LtNUIniYYTKOGSzlT3bceVWCY8cPc6/s2J1toaoGrf9xttZaKZFoCITRg5tx7tH7xTtOZidKUcgLqd6abe3KkOiiFVJaIc3b1Gev4/2oguPJ97lond1QQXd/A5qKzpxDIcMLNFVpCAlSl0AdFq/fiffd8Cx++uRbYAhEqC2NRRQLpBQriYjEStJrUFV/fOG6xO8P/eKFVB0xpla1XG414ObyTcVC5OvQlXOiKxusxIEx9iwAdWSvNM4HcA9jrIMx9jaAJQBOJKKxAIYyxmax4Eu7HcAFwjW3hcd/BHAWZeFhM8IuQmB4fa3e8Ysonhh3zn5HW69dsmJ5QVIox/3FMLGqrs5mbhY/5ko2GbYI148k2jHXgHMYrTA7Vo8hfZzFiU8OYqjzjhYJUoGAuz5xUvQ7MuUESyzM4tj4mFTvIUumNLIs6AMbi+7+OIamxg5LblxM75iHgxCdyHjq17krt4XitqBcDI+fCCUjtFfQvEfVCK644xXtuDieWxx/p0bjiRoThz1dQbSFQc0NkbmxHD6k1qhG5/BZIpofip1GhGXjAIhhRleFZePCY7k8cQ1jrBvAdgDZXTcd4SJCED/ahdeemzgvpgldtdUeMpujK7zopzOnpPrjMDnWpNlmBZg9M5yubxGuu0lTG5XUr2Y3lnU7ESv53MOSf+X++U5tJ98pw6kTRwvjjDmWJLESxUpBAw0K4iC/G5MVru2RqDgH3aJnUkjPmLyfpacYfBcs6tfk8Cq8LzE8vjZZT0EjVsph8a6VCNQFPJbTgKbiXkccbgJwKIApANYC+HFYrpqPzFBuuiYFIrqCiFqJqHXjxsps9V2Uj+I3OUhyEhJDdrdkyCPLrzly7NBUfxwmK6i0iZ+6nktGsrhv8wfnAtcFPRbj1EaOq/eQFuokxhN3mkf8qER94VjmIGOFtDyiNBoUIUPlIpdFUHd/AXFIlulDosccT2pMMoExjIVzDl0KToBCYw/el/itJXJcCO3FOoxkr4s3uDl6AsDE0Gtbhuk+ai3/39MREIdBTQ3RGrR7byAOjLH1jLESY6wM4NcAuB3eKgAHCFXHA1gTlo9XlCeuIaIGAMOgEWMxxm5mjE1jjE0bM2ZMJUN3kC+bd84i59DcqDd1k6EVFQgDMtnbp2SqiqnLwJxCSPMu9ZyDtQmhTze4KNVqzaonFMCRCaTbhsHWnu6CDul9xHoEpr0+misKncPwAe4JYEyhxoGAM3EVU3GipDqdGqXheXKCl7DMC+s/+9ZGbN3TJYiVYrPxhKe0OC6BExPxlT+6cXkAcNi+GuJg3MjUdq7y77ilsRCtGfUI9ieiIuIQ6hA4PgSAa28fBDAztEA6GIHi+SXG2FoAO4loeqhPuBTAA8I1l4XHFwJ4mtVwlXCJyrp1jz7kMN/dABll8+G3YIqgaXr5tvAAcT37WPgzyEes5KhziOqbxlU53MJniCaQwf8ys3/mugVTfJfJxSs+ljlPnglPVkgn+gvnikqstN+wFvzm0mnC2PTjFpijsE+GX/09jpVVoLQTnO79mPJRZ+E0OXH80n2vxn1KrfL5V0yIldTtcaJVzcKpU3GaDE5qvU7zjWJRIOD1tlayBlYhot8DeC+A0US0CsC3ALyXiKYgmCvLAXwKABhjC4noPgCvA+gGcBVjjG+drkRg+TQAwKPhHwDcAuAOIlqCgGOYmceN6WBby3a0deOzd89NlN1/5ckYOShQenJTVsZYpgnJFxgTC27kHFx0Do6IOAfN+VroHHiTT76xXlunmj1BZp2DIOay9qsjooHGFEAyHStjwAkTRuDqMydh8rik5VcchsIeeE+3+ThibOxNbZqDshHAq6u243tCOO1iIU1StSHRIy4kfV5uw0RuVXm/tQRJUEgngxHGdWxOcD+dOQXPL9mE+1rTYbc5dPPdNCtqtX/dsrsTLY2FxByICWN9qYOVODDGLlYU32Kofx2A6xTlrQAmK8rbAVxkG0e9sLM9bXI69aA4OJ7o/m/aWcjgNdMy4/jYlLTcJT59duWweSFwgesjcNnZ1+J7S+gcRLGS8B7z4BwSfSJwSHzPYWnRJ18cg9hK5kVJxTkAycVsriHGlGystKczKbMuFtyd4GxhL0SY6qjuSa4fiZUKFM1RHRHkYipdfKyRg5owwCL+1aXTNotAjU1WjOO/8wSAWA9SJAIKnEDWpk8d+p2HdLUvVRRHZHlZfJKnOQezzkFvqlf5jTDpvwy+G5s8bqimhn4c2/do/Dkcxpvtjsjwy34FP3axVtKd15qyMv3C3yAopG1hxrWhuIXiV03EIWov+N8lhaFXEQd9aGttNynoLOYe/OypCSU7l6vLPfJoAKJ+TzeuoaEzoy6zHYGs80rLORgmRq39HJaECvVCgayh12uF/kccLFPFtpESTRFtE0TkLPhh2nQwPv6yIIflKAqKUxuyzh292SK3AHFpI/n7Ww+mnQeBYPwrt+hDFQDJRYUxN+U6R9bAe8mdsE0PlQ3MMJ6GiHPQL0r8Xetk4c5iP6me7ERVpLRCWi/icacOSzeqsw8eM354gnNeHuZllr+jljCOWIFiXY6OcxgaBsfboeD4eRvWRdXB0s3l3GOvrcWGnWmvdhFdpbI2TI6qvEi9OHxGX4NtotjEH6JJpu1ViUrKWOegH48q61pEjNIxiauGzVpJ/GinH6LOOyF/2H+Zp47TxBjwvUffMI7n2w++Hh3f+eI7OOIbjxnru0D3lkQnONucyDPyrI4TVPWnzdOQaTTxM5D9aAoZrJUOy5LsxwBRj6KLgNzcGHMOsSmrur3mhsCaR2vmSQ5GKJrT8rM4+8h9tOfau0r49J1z8NHfzDb2dcaPnsGpP3haea5d4QHdkwrp/kccLOdtGyTRGcg1cRCg3g2+93C7OW6DRuFWzTyxKaS5043Y52VCwvpq+9VhQFM8Hf/n1bWGmqr3ZF8yxRqJ8BmW6yrxc5B32ndeflKyLYVY6fYwKxuL5oq6fdsunp+OaoXtdcpiJZW1kqK9S6YfhEk8g5uihsm6T0YihlS4Fsp6vphziJWxekV5sJ0zWVnZ3p9ODGa6TD7Hv3WbY+yqrW1Yv0Md90yVfTIg4EnDgnqh/xEHywMWT1995sTUeZFzWLDanCmtS+EFKu4GGwoF68JUFBax5DgrnyhxTmT1+bFh5ilRpq6znHHdVTMwa92WBne/ERkuTnDiIhbb7bvoHLI/a3k47wrjVUXPHukxH7X/0ER/LjoHGdecdwTmf+ucRD0++pRYSYj4ycEUEo/BLQ1Gn4n7Wlc56yQSnEPY2P+VfBJizsEePgMI7lP3hkznOI49YLiyXGU+/sdPnxwe57dZ41BFSCiQmxl4LdD/iIP1fFzj5EPSUTxEccQzi8xe2nuE/MJcLCTH07EtPAUN56CC6yLGiZbWcza07li5Jd4FqUwQgz6dusSoQU1WtpiH0a4V5qyIlbfJwHv56hxUDyUt208rrWMleXKMMkw6h0HNDRjS0hi2l1RkphYflhajKvOYIE1osoxJV0+nR+CbBCKy6hyC8em5AxNX0dRQwAVT9se3/vko5XnVdeNHDFSe499nVpGfCNU98jQBgNc51Bx2+XJ8bErm7iL/W7gm5izED56Lk1wm0pRwV+MSPsN16nCbfN2zUHEJes7Brc/GYsG6CJ8/ZZzxvAm6Zyn2KEabFa2VdGgsJhdXV6i4Ann3p6pTkMQHWrGSoW9RqRkv6EF7slipravkFD6jQKLFjPphuBMHsS8NcQg3J4WC8LxM78CgVyBBb5ECA/YbNiCVSCs6bSDyKU6+HPdXKZTWiqLOob5BWfsfcbBapggTwpTM3WWX/vFbW6PjaGdRAG657AQsvu68kHNIX/fLj05F63+cjb9cdSq+c37gGiIbMqh6d13E2jQmhBwqW3RV+Oigz7iVWUs3a/sMrLvM48piLilDFYdIxplHxApFMa+CbliVZhmbv2p7KkQ7Sbs/xtKLvLzxEBdcMfKpaSEWd58yQWpdnoxK01UqKzyk1QuijXNQDUndVlyxWyFjB+KkVUUhGoEpPwoZBqb7xoJLmHHO/ejxtxK/3zV+mDaab8mBc7D5Rek4B5N3ei3R74iDlXMQFmFVWAAxJo8OKvEIEz74YoHQGKb/U+14Rg1uwujBzZhywHC0hEpa1/AZLog+Sk0jKi7BRay0dKM+2FngF2IetPhMs7LQDrQh8eWKRF4rdijGz/4nT76lriSBL4hzVyT9D8SFuqO7hAdfXYNdHUlTXQrvQWWt9OjnT0vVUyFBHKTX+IBkSdbRXVZYK6XbFEOV6y3ckg1t2NluDX5n5RwEPwezJ7iBaEE/l+TIuDZcfeYk7UIdjc/Qnir3hwgr5+AV0rWFkTul5AswJXM3mXmecmigq+DJz8X6YpMuE1MWNejwxbMP056T8x3wxf+pMDsbkMyXq+IcdGIl8cMzZbJjzO40mGXyt3VKAe001EF8bqJ8PWmSnF6Mgdgn4aW3t+AnTy42juf979ovbE99viBwKq8sDzKQbdqVtFqRd/rigjtcSCNqFCspBqB7qp3dZaNTZmJsUT23dzT9v57COTc8a6yjW/C5Qlr8HpNhvZPXBToHdVsMBs5BofcxIXAaVI+BKb5vGabvA1BHSBBNWevNOvQ/4mB4wAQXnUPwX1zI7rni5FQf+w9rSSyovLbJCU4FnY4jHaxM38aXz0kSDj5JRYW6eD9FBZegC+Ugjt/ktOai+M1ix/2Fe+clfrtwDqrwGeUyopcjExj+/rLEmrLFJiozhg7NDlJWPFbiBJcUK5k3FicdPDIdF0mrc9Cf/9FFxyZ+H3vAcKd3qRMVNYse0jwqq41z0JwulfU8KEN2UaYujIgtLXBndxmfvTsZTv+J19fjx48vin6rCLtXSNcRpgVKVl6pPkKXGDNclivWufPFd8I+xP7UmwGxVz555Y/j/leSgcS4uGvqQSOistGDk+GdOTq604u42HyjQ5KZBoV4zbQzcgtTYa4wsKmIDxwzVnlOyzkIx6rMa6Izo9wEv8cmh7wdfOg27ocxfUymtM5B3YaJViUWmGiXq647ad8hzpsVk85h/+GxPmTkoCYMbk4qeGccvR+e/vLpqev4wn/apNGJcq4gLhbcdQ66s+Uyw4RRgYWR/B4ZQza5EqDXOZTNBH3Zpl2JLHMA8MnbW/Gzp5dEv1U6mEJBrwSvNfodcVDhkc+dhoXXnhtyDhaxUvjftAi0dZXChT+uw51jEguUgR2OxqBRiu7QxJIRFwceWgAAvvfhd0XHHV3pRfyjJx0UHbtYK4kB5DhU93Lfp04Ox+8QbiQ8v72tCy8vTyd/HzusRbsK6MReNz0Th6hWifREscPVZ05KXMvflYo4HDJmkHLstsVVPlbV4QuN7p50nMOQlgZcNDVOm2Ja9rgYUtZfaDkfQ2vieBqLlLKqOfPIfXDImHTOBF0UYv68k5ngtN2nNmIiSozh8ncfAgA4cUIs/nURA+n6Cq5PlvPfunerU74nxqqMrUZOG9JaoN8RB9XzHdISZFsKZJxxueqDcDFlbe8qa1ndxPeu4xyEGeYyMa4574ionmjKKH60E0bFi5mKc/iPDxwZHasWJXlB4mae/Dks2bALq7elvUMPDRdRJiikD993SKpeUCf4/8nbW5XnC6RW4OvGDEjEQSFWYoLOYZrAdQHJ8NoyTpcirtpMLkURj64O51SiXai6mhbPX3MmDhg5MO7TYe6kxEqWPlTEQ5SLNxYLqXfUrOG8dNY7fG7pMsHJsCmdBzQVMfWgEYk6i0Kv6KyWp1bOQXOdXP8L98xN1dndmd7wBYYr6jZqjf5HHBTPly8sBEq8ANXEiW2O9S9qnyHNWq4gyTkg8TXuMyRQHB8neGzKOo75q7Zh486kIvOwfQdH4xJ3KOJ6efKho/DYF07DcQcOV4p/CgWKQiQ0KHQOclG0m2YMm3d14Ozr/44/zVmtvd/2rhJeXBaYUo4Y1JiqF9xj8H+hxvN88+5OPLJgXcJ6hkOnExGhIrpiJjjZOi1O6Zl+jycdPArjRwyIfvMa9nwIhhDgYX/dFhGF1nPa8Pvvb6kdNlNRWZUK7Thr3Tub08ETxcfWVCykNk5NGks3nahIFZXVKNM36BzKwqIt1pnxk+cABPlbdDhqbDoqsU3noHtnMpekikH29iZ1wELR7Lqe6H/EQfGIo4WFZLFSdp3Dzz9yHH540bHBrof3qSE48kQaNbgZZx+5b2KREhexOSu24oM/fx6X/valxHWlcvyBdpXVnAMAHLHfUDQ3FJRiJSCe4KqFVh5rbOZpTsnIx7B2e3uqTIZtZ8SJwr0vr9T2Y0KSc4j7ZFIZh7yTFzFj8n6RVRqQ3uF+/NSDleNjMPkKEBoKFFmt6D2k9derwMBwmTRn4ms0g9H0qRKPjBAsqRqLhdQN6nQ2WmulohCVNaxj2oyZbiHWBcTfrOiPs6sjHc2V9zWoOe0cx+eEbJZqE1O55H7RiZ7ijUV9yUNt4xX0Rhg5B9nPIV2Xl+kWsn86Zv+gLWHXI+7UTYHOukvliKWO68f9rQsX2DfWJqO3lsqxSZ74walELQ2FAtpKaquiskHWLX/IjQ083SXTxtIHYpm2+Li0ZrGOk1+1oDlxDok2RM6BKWrYRYii2DFWSAf/xwlchdiyLfNcsUDoLjGUywzPL92krOPqhStyK7pzrvkceJ+q+D+iPqGxgVJtuCj0VfULgkLatLaavKBj57RYJLliS7xDV3EOZcZQACn7bG4ooLGYjgJrcWEwZnmM21DX8TqHOkH1fLmlS9rPQc858Fj0HAeNGpj4HbCxQVu6nbpYBwgmh7xwRoH3FOc4Bjc3KMVKqvryPYowcQ483hKPN9UocA4mEzvd7vfuT5yEv1x1aqLM1RpDDKDH9RcueYzFRVXcjTGpjKMh0quYd3RiHVWAxWR/ZuVqY7GA7jLDrS8sx0tvB2K4P3/mlEQd11DeJgsjve2/uozPw06NRRrvu1gopImDJFb66PQDAeiV+LHOIVY0i1ns5Po6qz9AECsJnENR2PWpiJ0pwRARYXBzQ2pDJHIopnGYoDcGMJ+vFazEgYh+S0QbiOg1oWwkET1BRIvD/yOEc18joiVEtIiIzhXKpxLRgvDcjRR+qUTUTET3huWziWhCvreYhFJJHD4FknYLJie4vy5M5kJ+8LPvxmNfOE2oGE/Y9lAB/J0LkllS5UndVS6nPJHF3avKBPIn/zIFp04cFdXrLOm5FF4mPwMeafKmf52K86fsjwNHDkxdd8DIgXjuK2fgd/92As48Yp/I+onBvKNRhyAhnDJxNI7ePynTNU3+86fsHx3/7h/Lo+Ovvf8IfeeGsSR2Y2G38lAP22eIdVwcsc4h2T6HKDfm7R2pkGkHnEM5sfk47sCkolzHOcjFJj8HnehKmwkunJa6HfC00BKouZiONCxzDp887ZBwXOo+uc6LR2V9+s31+Pw9Sb8WEbI+QUQcKj+uI3LnvEw0+7ZxK4OaG1KcQ3wPbjoHETalu4kDrCVcOIdbAcyQyq4B8BRjbBKAp8LfIKKjAMwEcHR4zS+IiAvubgJwBYBJ4R9v83IAWxljEwHcAOAHld6MC9Q6h5hzSOoHVJxD8F+OQT9sQCOO2G+oUC9e+blz2EApl608qUulNHcgipWKivhGFxw3DkSxtYjImupiQzHG8Mo7cZwd/mFPHjcMP515nHYXfsDIgWhpLOK3HzsBh+/Hk7+YXXNUHt68eZkQmjZXoh+DKO4zOUfJSHpI8z5jayX5ff+fcw+P6ijbS3AOSNTVMjIs9hT/seQ8BgQLV3eZVRTdMxXl1cA58CB8qZDdurYV1nAifvGvx+Ohq9+NlqZi6j3K75mPk9fTm/8GYqVn39oklafr6a2V0mIl8RtTcS8qr2wRTQ2FlM7BtoEwKdQjvYqGg+q1CmnG2LMAtkjF5wO4LTy+DcAFQvk9jLEOxtjbAJYAOJGIxgIYyhibxYInfrt0DW/rjwDOIt3WKAeYrZXcxUq7dJmnQoht8SBsLTJxkCZ1V5kpdA4UiYJMcnXeh2jJo+I0uAXI25v0KTtdHj6vw5j+IwLUC5RWcWpo5xghoqourMTvPzldP2AkF5WkziEsQ9KJkC9qOnnyladPxORxQzFh1MDU7k/NtXFOK1yYFV/fpl2duGv2Cq0fC4eaq1XXVWs5ivYAACAASURBVD1Wvsi7BN5jiJ+5mIf6kc+dFsV8GtTcgMnjhgXjktqQTVllc1BT3mqVCC7LDpqPW+QcxO9I1VS8WKs7aiwUErlagPi5yCFROExiJV1Y8oeuPi0aO1B/hXSlOod9GWNrASD8z8NdjgMgmpKsCsvGhcdyeeIaxlg3gO0A0okUcoLq8fLJIiqRAf3OG9DrETjEyfjBnz8PQPGRSNeodA5AHJ1S51kLJHMvRGNV6RwQTHrTZHWx/HF16Y93iXE9VRIlsc5uKW7SxSceiEtPPij6LT4HcTHn3Ix2LAmxUvBf9DsgAu6/8pRUffGj3HdoMxZfdx4A4MBRA/HQ1adh5KCYQPGqqmfIWGCOqxM9iVhmCGIIAB88dn/j+WD8+va5mEOeIv/0s38AAM45al/8nzDsCiF+5uKieNT+Q1OisWB+JdtMiUqjh5/4lwL3kLZNR5EeMcYSG6Rzj943qCOIcEWdg/juOVSOd3d94qT4fhooQSQBpIiFDJNYSWeRxR9TX1FIq14jM5Sbrkk3TnQFEbUSUevGjeZEOzoofQ8ycA78g7OJM1QRV2XOIRhPfNxVKivDQPDdvk4hDQAXn3iA4jrFuEKdgzEcgQPrEO/+1C9r8rih+NUlU4XwH/E5UYYuLi66Rzppn8FJH4UKrZ1EsYsYQPFPc1elzge/+djjdpsaCmkxCVFqF6x7VXfPXiGIOvSwBWnTcSZqpJ8LV8QSEUYJxK093PR87qxJkYd94KXLrzM/Y5WjYkpUGv4vM4auUlm7USlQ2vIJAE6ZmNw7igv/vS+vxPHfeQIA8KX3HRbpL0S/I3E0qvZVOoBTJ8YhPhqLSc7htdXbce3/xPnPVTBtxnhSMLkKF12Jz6ueqJQ4rA9FRQj/8/CeqwCIq9R4AGvC8vGK8sQ1RNQAYBjSYiwAAGPsZsbYNMbYtDFj7PmXlW2YTlJysVZ9bJxgxE4vmqYI2LqnC1fe+UpUJu4wVf2VymrRERcrmYiDynGtiWfUSow/5BxMxMFBsCTuZlRN3Xn5STj36P2EAHf2ie06903PwQQV5/DO5j14ZMG61PmgTlCwURAVKA0ahPcoKkBVOGTMIKPoiXtp28I7ZxX9ybj4xAOj438+dn8MlcLM7zO0OWHpw8eqs1aK+qR0+AyZmIs+H5P+/VG0vpMOlRK1xZLz8R9fPQNH7z9Mrhnd498WxZGGRSIuEhBx7sccR9yabeMni5U+fNMLKfNyGSaien4oWZA3bAeFUQ3E51VPVEocHgRwWXh8GYAHhPKZoQXSwQgUzy+FoqedRDQ91CdcKl3D27oQwNOslsI1Q8tuOofgP3/ZOlEPEeGlt7fg0dfWRWVjhiRDZ8uLcHeJKRd5Hp3SJCZQrZdfnXE4Ljv5IHzo+Dh0OBd3me3G9edkiE5kprZsUSuBYMe2coteF8KRJUpqYiyKNlyegxi9VnUbomd9vKCmx3jcgcMxbnjs/6B6Z1efFcR3si3CKuqQRSEtLrCiJzJHc0MxwdpzgmwTn/CNjAidWa9tJ1xQtqXjmNJye1mMy5tK5qMOjkVvd5MpK5AWK7ksV6Z7Xb2tDbs6uvHzp5Nh4Qc3BwTb9XnlDRdT1t8DmAXgcCJaRUSXA/g+gPcR0WIA7wt/gzG2EMB9AF4H8BiAqxhjXIB8JYDfIFBSLwXwaFh+C4BRRLQEwJcQWj7VCuJSdurEUZj7jfdFv0XxAGBWSEemgJpdrKpUNunj1lHtXSV8+8GF6CyVlZxD8JEEo9dBNdaRg5pw7fmTE2kQOdvvokQ2IdGfoq14Bx0q1DWrsMziy6G4gWBnK6Kl0R6BVQVVVFaRFqusYFwg6pdUfXFw3VEselIsdOH/A0akzYmT9exiJTmHNMf1/ztpJVUspBee5oaCII+PxUpb96Q9ikWoPgd5AxXrctJ1xXzOXKwkXq703QnbeuL19QnLJtHDOTD+CCBuVDiXc8vHTsABIwMC8WwYakS3cWgsFhJWWy5rti3w3g8fezPRnxjY0fS8agmrhzRj7GLNqbM09a8DcJ2ivBXAZEV5O4CLbOPIC+IDbmkoYoQg6gkWsfi8ySKEzw0955AuS+2gECxof3hlFW59YTkAdTpO7ikqTw5RSWYiZHJZmdni49sXxWg3U1ZzDjKRdeEckt7KMWSOSxfIzYakWMmuO1I9BrVjVDxvTDqHAhFKZRbVNb0zMfmSCqpwRXJrvPk125PGCh8+fnziNx+XiOaGgiD+iufE719aYRyXHJ8suF7WOeif/fRDYn1CQSWi0nAOjKUDNm7aFSumA+YizV3w2Tt6cDP2GdKClVva8Oe5q/G/po7X7tQbCoVEPm6XHb1t/rdJuVD2Eea8LS9HrdDvwmckdQppRVniRZt0DuGszRJW2SVshJpzIKUo6FgxQJ8q1IfmQwp2r8qhhNfpz0XthP91zYwUTEwLpJe5iqWubLNu7OLl5TJLcXWqwHuieIB/hLdcNg27OrqV3J+SOIDAUE6MTbnwFwIiIi66qTphmU3nYDKWiMcVgIeL17ZVMC/oBHKaE0Fb6R2uPO+jOGCKexSv5WIl8Wo156D2c+gQFlxR5yASB3EucW4gDriYahIA0NRAibG7uNrY9BgqcTJH7JNj7ydP9D/iIBzL7yOLWCmOnKnux6TMFuswlvwQ1ToHdUwekZC4EiNurWQOZGZfCRJpNoWmWv/jbAxtaUwszESUq0Ja96HJKUsHGHbffHSqnMtnHRmYP8pesLoxFgoAF56aFn7u/RyZuxoCHNqsle5RBB/UcQ4mE2h+nj+HQ8cMigwnYrNcdz2P/A3x6+U6gBsBTIvrVH2q38v5QppeLnoCpAi1wnV8o9BYJDy3eKM2SqosVnKBOM/++OmTceEvZyXbNFBf/h1ncfjMA/0vtpJh8Q84B2jPB2XBf7spaxrp/oIdj1iq4xxK5fRuv2ghDiYPaaO1kgvnoJGDjh7cnNKtmDgHEa6cg00pCiRj8ajAx9+tIA6634CaUxJFKXFsJc17TOgc9ONyuUdXGDalwTgKFIn0CkSRGE8kdC5iQSC5CIvtJ37zZ6+Qw4sEvlAI+hXzHOh8d1SjO1QICCg6nC5atzMqF+ccN+9tKBRwyS3qKLb8vOznEJ9TfzzieqEyaS8a8twSEYqFtOiv1uh/xEE4Vu/kReKRvj7l56B5Xyq5vSo0BmPJcSiJQ/jx/mVuMl+CuCNUEgINoSkze9pFG0SXfpsjnJgL2NSXK3HQ5qoWLpdluKl+Q0W5mNRd5phUC/wAxYctiiyMYqXw2dvqAJURB7k5V+epKLgj446YPMsfbzf9/mRPfrFPuTvVpgiw3yNXSP/+pZhLUnFBqnhhIwYmc4aI3tYniBnhhDo8zLoorgWAh65+d+K3LFYS0V1WG3uIC7vsJxOcT7Z30dSk31KxQE4brDzR/8RKCZ2DfDZp0meKrWQTkzgppImHJojLNisS2fDd/r2tSVGCKu9D8jr1Lkul3E6Oy1WEEHIhlnVsT2cJT725wVwJ7jJV28IPAG2d9joFSpokuuzM7rj8xFRZIp2lUSEdzJs/vBK8RzV3Gfy3mrIqIL+35pCQ7bE8C75WlcoM3YKvTWytlH42r37rHM0Y0kQ+Za0U9ifvvqceNAJHSvHJbPoLDpvjHaC2VhIX8tPC7H4jpWRUk8cl/SpkJzgZQZDMZJlKfClCnNPLv/+B1Hkxz0e90O+Ig7hXUHEOZQvnwBdkTsV1y4lqCqcXXYr65RATkYjjtNuEKzgVwy7LqHNwow2CojyfHY1rM+2a0CWiVZNNng0E71f8aOXwzfJzuPPykyLHJLkdPnjehE7nsEDIcqeLmgvExEHeAWcBt+oSRWynTRqdqifmKRdDuIhB68QFtbmhgIFN6qVDtaCnxHXhf3GBPeeofXHzpdOkttJzS/VcSSFXUnPp/B2pFdL8e7FNHZNYCQjuq1hIcpji81N9XoeFoeePP3C44mzPcA79T6wkPN9pE6RQyMiuc9CZl7nuvmWF9FfPS4egDnQO5nZUIksd98Is1kqOtCHiQkSTwayYPC7eKTJmdvSzjke41mWTRZQUD8gfn2rzoGwH8bzh80IlQ5bvzaxzYBg9uBlzv6neobuA637EWFWjZC99iIsiQ3c57WtDBMnEW/+O+JwQPa7T1krBb9GrWC+GM3MhfHzydE5xK8KxuCH4wtmT4v44ByX0yUPTi2g0iJUANdea5BzS98ADSP505nHKNne2d+MBRWrRWqLfcQ78Ff3qkqk456h9E+fSOgf9zs7GOTiZg1LQgthPKsQGuHlgBZxDhToHV8sULl/evFsdiVKFr8w4PPH7P8+fjH854UB89u45KbPFauCiQCUkP1pZdJLi83TEQVB28jStKn8VecHSxdECAmslnVzfFRHnIFhdmXRhJYlzEGMRlSwctdg+Y2aRJ//5gsAlqzmCtJ+D1pRVet8HSDlJRIumEmMYN3wAnr/mzGTbXPcizIPhA9KcW5MgVlIRiT1dJYyQyvjcOnXiKGf/GRmq3Om1RL/lHCaMGqR0zhFfknrnHe+yxPZkuJiDcmmEON8btYH3shMHtWhLbRYr1+FoaSzg9o+n5ey8A1F/oUpeI+OQ0UmxTEtjEVMPGhERrbzgoj8oUJJVl61ndIrUdDvxPOBtKN+jVKTKzyH6AKgITBZwz3iRczDN6bJG51AgwqR9kpY/OnDOlKQyVX9HCFF01ZsblVjJjXO46aNTk3UEXwjG1ERGFqcFx6lqYbTY4Hl9/NaXU+fbFJZyfJ7deflJ6lnEn3WFccNqgf5HHKLELulzsuJNNWF5ykOtxUzUmH0ssUI6rqz0kHZYOMU59fFTD8YPLzxG6TPBwwjwyTrj6P1SdcRncNkpE/Cew9RBDju7y3hNkKHLnsxqqB+MSpFZDX7z3DJrnbauUpSKE1BwDrKsXPtO4/fTbXCOlOeTzmwZCJ6tirPIAj5XRZ2DGNuJQ7SjDxJOBdd9cMr+aGoo4MPHj8P+wwdE+goTYxm8x+QiriMmiV23hmi5zAm+yeI44/AxKQ68oUiRHkcXGl8VcXlHezpciBgI77nFcbiOWMeTXhvKZYYCIQxFku7bJVIvkK+Jsw39jzhwCwzFOVnpoyLizWFcn1ispNE5OIyFOx+Jc0VFHFwWTnHHcdCogbhoWjqENxDv7LgTz08vnpKqM1iUF1tETM8viUUDTQ47XV1zXJFpUrodYcnXIOLR19ZZfR1kDJeUvzqPYxlPvrEeb6zdge5SOeYcNERehCm6bJcmzlYWcNoiWj59+vRDFfVisVJ3mUVz8KBRg/DWd8/DIaG/AB+vSezIY3eZRi47kuralPOr6ECU9AEQY4lxjBs+AGu2taNcZigxpvy2VZzDLkXSJZ2vEyfGqjlcYiyRcVIGv8Im0q2nr0O/Iw4cqnfQWCygS/iQVC9KF+lRhovcfkd7N7a3dWHJhjixSx5iJRNrynd2XSWGIS0Nyg+puaGIY8PMa1nCY6vst03jTJYHH+W8ldukscRt/mRmmpCZcNQ3/5qpvmyyKMOmLN/dWYodqRTPQn6Wukx9AMJFurrPk/cnelur3hEfFmP6XTUgJsXS9xlzDuY6ABLfmqo638jYQEhap6kCMw5paUBnqYztbV14eP5atbhIYa0k544AktZdIjhRVS3gpTJL5KpXnQfsloL1jMzabxXSqulYLBC6LPbIsvfvbRp5vItO99VwIfzZ00uiMp0i02Z9I37Pph1nsENn2Lqn07iYx0ne3YmDy2Kma00lOttvaEuUhhIABmnMJ+sF26Po6C5FJo66vBwilE6KwiOsViHNF7uOxIZHX2/W0s3ayMDBtXbOQeWQpmuny8KlF4jcLOEoKW4ZNyItOuPP+oYn3wIAZWgM/uzLjGHssBacNmm0IneE3rkwTiubfgDPLNoQmWCrHl+ctjV9Lmg78MmppzVrv+McIgsMxUtoKCRN1FQ7RXGnfcR+QxJRJEW4rKmq2D+qBTurWMkkCioQYeueLjwwb43R+oG3kY1zsNfVidFV97jfsJZE1Nxa4qcOXInt7to7RbGSgnNwmBQ2b3mOYZIVjVLHUYj1F6r25Xo8XLqWcyhy4qAdlrDbNxCQ8P+2PULUVMW4Nuxs13cktSd+t/sObUnViQilIb2vaK1UKjMjlwsoOIdCfL2Mt9bH0gElcQiv0fX51RmBiXs9I7P2O+LAoXoFjcWCNTe0KOYwLZwu1kqq69Vsv4tC2l2s5AK+iLsQB/5Mvv7+I611dc+FKG2O6LLTzgv7KxS1MmzPrr27ZFRIuzxLsYZJIf2rS5LWOCrCw8tE4qC6h7SJrXqcse+Ged7bNjKRWKkkcg7pNtftcDORJiJ0dZvbKjqIxCJLxNAPSCdG5PVkxXPk62K9/3S7/FnohhentTU2nSv6HXGIQwKoJxB3Y7/7kyelzvM6HKadncsirKqi3gG6Zc1yGpejJ4GL8pFbvhSI8MnTDsbowWprJbEJvUI6zaannLEsYzZhwqh08pyJgnmmm/jMXKets2RUSLvlyRDml4ETGzEwyVGZQrZ3Wrhhec5t0yT0aYjmhHZYwXsEYEq9pB6Dui2OscNa8JerTlW3h+Q9qkV6lPivQkHY+TPGtAELeRtX3TUnUc7HUElQTi5WsnErnnOoISJTVsU58YM+SmOznwh2ZyQOyXMqT0vX11wku+u8K+fguvnm7ZnUCOdP2R8NBUJnqZzSxYgQn5lpNyZ/VHL0yiz6DxnjFZnVjhkvpspUXyeGM9B1/41/CrKXlRmLFgiVot9Fv5wg8g4x/qO2leKi4L8tTpM8X95Yp86H7LJh4EH6su5wbaa/HztlAqYcoA4tQZQUKym5NuJt2sfBI+faFuqXlidT3fPnXIlPUpdFIe0aRDFP9D/iEHEO6XPipNKJABJmpxkW4ZknqE1LXUBEaJUmYro/YewOuyMbGhwWAu5EViozNBX1uRPENkxss2zDbfM7yALVB9vg8L6v/WCcvFDX/aFhSscyixPMqLLVuRA3sY4pxn86FIeeI1hhycstz5cdbWoTYF7PZonEkN2qxhYbzJT9j5CcO0qjDof5zPsslYN3qaurmytfPucwAOnYTDZxKSAqpM3rzufumYsTr3tSN/xc0f+slQxzVpT364lDXC6bXYpIx6dJt+f6/RQoKd+8/8qTccjowak6HGZdiBvcbNrj48YGNxGbaTcmBzOTiUU1xEH1rMX4R9pxJXJM67mesJfIMkjFSYl97D8srTQN+oiPTWIlJ8snxwcmMygff/cEdb2iy5ywB3ZUj0HP+QBAk4ITEyEGY1TFteLPx/ZIuAhXzl0tQjUPvnD2JJx0cGCcIkdPlee1qtmuSKxk7lN0uqs1quIciGg5ES0gonlE1BqWjSSiJ4hocfh/hFD/a0S0hIgWEdG5QvnUsJ0lRHQjVRN9zQL+mlSyd1vyHBmmyIyiMvHaDx6tGYvbBySPZdiAxpQVT4JzyLDjtPXp2laTUQQiipX0dWTxh5zr2FVfooLqWYtD1hMtcU6o2+ZVyiw2G7VxDj/RBFhLesu7PVPAzSPbpa1jxw/Dh44br6zn5OcAe0h42xjithw5B2nzpOLoORdieyaFUCwm51lJ1kmXXTTtAMHTPHkuHe1XIVbiCukMfdYaeYiVzmCMTWGM8Xi71wB4ijE2CcBT4W8Q0VEAZgI4GsAMAL8gIr4duAnAFQAmhX8zchiXEiZTVpGNrdY7le/6xg0fgMtOmVBVW+lJqviQEpyDqS23PiPOweIzwWH6eMUmdK0VKKlUfPJL70mJ4oYIntv/etKBePyL79H2KUO1WIkEXGcYJC66+thKsTz4+icCO3q1E5x4jbo/8XlnEVu6Bl5UIbEpMlzjZPFTIOzpLGGnIsWqiE+edrB2DHFb8XGzwrFNhOgNbzLrtaFIos5BXUf1rMcNHxDnxZAmm2hJBehMWS2cQ24hKd1RC53D+QBuC49vA3CBUH4PY6yDMfY2gCUATiSisQCGMsZmsWDlvl24JneYNjQi1c5i368CF1GZ2nHdXYmhjQH15EqO3X3Hqa0XjtvsMxEfiyE3TH3qPnKSOIeJ+wxJ7aJaGouR9+shYwZHMfBdoHrULnqaBGHTfrgBsigidTtEUdxkMmWVFwubtY9xXI4EyUUP5frVZOV8VAr+qE8iZ87BGt24wBXqdoW0bryySE3OLaK6vDsyZe07nAMD8DgRvUJEV4Rl+zLG1gJA+H+fsHwcADGV2aqwbFx4LJenQERXEFErEbVu3LhRVcVtxFB/6PwFELmLX3Tgkz0LkTlV4aoPILULUwVPS/Sdg1QudoLT1xEXlQNHppPgRBCGo/vIZc7BhqycnSq9pyjTt+sT9IhSptp8UZwMHggHhuGmTU6F8rB0iZ1ckEw3a+Ic7F7zcp+ih7uIdF5pM3EYPVjvDClfafIxuW3WO9p2eL1SmN1QrxxWl4sBDEXIujOltVJJL9HQXVNrVEscTmWMHQ/gPABXEZGJz1fdnc6VUvmZMcZuZoxNY4xNGzNGHSnUhjgqq140U61ICYg/bNM7veI9h0THZxw+Bnd9Yrq13Ze+fpYyQbkIUzBPkeXdxxBF1c1sMT42W5PE0I1dtFb6wLvGatuKwkhr3tE3Q7NSEfsMacaPLjo2Ve5ireSiL+Hl9pwb6mMZnDsyKaRTi2sVc1Z8DqY++UbB1JP8jHRh3FOmuIrpI74TnQ+Nqk8VxyXX0YW9KUTWStnESmJ5SqwUzusfh3NQFaq90xI+owekStURB8bYmvD/BgB/BnAigPWhqAjhf548eBUAUYg8HsCasHy8orwmMEVljWK+GBTNruAT1LSLP3/KuGhRPfPIfbX19h0afxgm8Y3ctwqi6OOOy9WOfoAgRzf0Yws1ziG2oU9MH9uJn33UPso6QPz+dAT8vHelQ5B/54LJynDiol7AZSHQenfzsambiOC6Q+ehFh5fuF5bR75a3GjoICe34Ug+BzvnYFqonHVaslhJyTmIYzSNK3nOZMrKceqhai69WOCZEvViJf645JSeUVRXDefQGH7rQ1vi0Ce/+7cTAABrtrUBcDOMAOrjDFcxcSCiQUQ0hB8DOAfAawAeBHBZWO0yAA+Exw8CmElEzUR0MALF80uh6GknEU0PrZQuFa7JHZG1klKslB95bnAw+wNi6xaTTfvMEw4EEOzOVeIRGSZRkDhxTYSGt2EySUx8vI65B1w4B5POhHN+OqKrKh/SrL5PMcS4mymrekx81/7Lvy9VV5DqmdoSsWGnPnyEON4nvvgeXHryBGt7OnGki38P4MZRu35DKT8Ni87BtMmSQ86ovpEUMTJwijzAnY4gqUKOA/HzkcuvuX8BAHVI+zEhR/TKO1sTbaf7TP6uRz7pavwc9gXw5/AlNwC4mzH2GBG9DOA+IrocwAoAFwEAY2whEd0H4HUA3QCuYoxxLdKVAG4FMADAo+FfTRBzDnqxUh6IlHeOWynTq+YTefjARic5sukDFcWfpqGVLeIbuR9XUZwqKBoQfLwuhNI2LlW5Tgme2DE7iJV0XA8vtdmgJ81iq5ts4uXVztuEpZuRc1DvjBNtCcdGr3npeS9ck/bKFhdn0zyUzahVGxD5epOOiYt4dFGL+bXzV21PlOtCebeGC7+qPTmAolaqJJ1o7yo5hcivBhUTB8bYMgApQS5jbDOAszTXXAfgOkV5K4DJ6SvyhzkTXJ6cA1feudU3VeMfkjyRtH07ipVMCxSv56p8NKa0dNjkFAox+22y7+fstC2stAid93ZC52BRMuraBjIofjP60ZiQJA7VtSWOxWVBN4ldxR3t1WdMNPSZ/P3E62kRmivnII9ZxRWY7ku+lnu56+eY5lqucxCeQacll8awVIIpddvyfNmyuxNDWtzWg0rRf8NnKM7lyDhEu99qTWLFNoY6TgajQtqSr4KD7w5Nw88qVrrjck0uaiDhBGfiQvjodTtJ1SKi4xzEBcPm4BaMS7eTTP7+0HFKYzvJU1zdnytcfUxkubi6rfjY1NZf5q4GAKwO5eMqzFmxNTq++qxJ+j4ziqhMxFRedFVNm+5LxIote/B4SKg2asR6ehFkmjiIibzEeT0+zDkhizxdLaQ2G8Lt54X+Rxz4geId5CnFM8nNM7cVTgxV/gcVTB+CaElh2o1F4htDHVeFNIfZ4iQOS25UblsU0qrHrnseSSc4O7FRWZkA6WekswIrOu7QOT5vWFzFHk0+AKq0oDJc/Qne2WyO0QSo8yfb+gSAf1eEe0865+nbEt/vhVPHY8KotFm1K3EQ8ea6ncpyLeegECuJ3tGjhPn/t//zXrz13fOcuT65Vpvjc64G/Y448NVFpXPIGizs8ncfrD3HFzjXNk1zJPZMdZtIpg9clBe7iJVMXSY5h+p0IeLlLtxWlqBogzQK6YSfg0M/unuUS3VvXIz/YyIOvMtzjtZbsIn3afIedvPTiI9NbemiooqwRYDlkDcmJx0yMlXHNSSMqAv65j8fpSH02Vk13YZhcLOag1elGeUiuK/MOByHCznQG4sFZ1EXkH6P7V2eOOQOk7VSVtbhyvfqd2V8Idnd4fYSTTQk9jlwG5dp0pUcicOXzzkcJ0wYgTOP0JuVZvUoNxOa+KSLos3FL4FDx7GcLZgP6x6/ywIl96mLM/Xqqm3WOkA8F1zSuAJAi2Ez4MLAihulVyUlq4hvfTDtQyKDL1oqv5JEn9KjVIf/MJ/nGCCkj9UR8EpyL+sYWJVZNBA/a7Ev/r0dO95OWE2Qb6vdkpQsD/Q/4mDQOWzdk02OZ/p4/7EksFxR5arNikJEHFw5BzexEhne/sGjB+EPnz7FqPRK7qrtU8k0+qxciKtC2iSaGTmoKfK81S0eoihJyzlIxbrnL24mXHaNLqEsAJsndXBupCHd6n5CyI7dhphILjm8ebIsnbc/h0xoVVPbVYEvhhzR1XPV14nQtTWoWWfgELzTbkFh3x2Z7bMp1wAAD3VJREFUZ1enZFJZK9Ua/ZA46D2kX1i6GYC7ZYPp43WVvXKY1n2XDFwi8hArucBVkckXXtcIry6ERqsjkMr/7dQJxnau/eBkjB7cjOED1IvnYEEk5ZLjA9CLZkSxjMscM3IOCZ8Ju15It6ABScJh4mhcOLpzjgqcEG2WdTKhUYl5eVuAee4MFfrS1Tt8vyFO3891H4qNJo8Yq47dZQoBAyQ3YNx6yyW/ugnyO27vrj1x6H/5HML/pld1rJAhzATXkBEuMIXEcA03zGFaeEQT9ZYKlHQiEtYkDr4JpkU/ERHX4UNyNTO06Wk+cMxYfOAYfbiOZJ8aayXpbX/kpIM01wsKaYeF1rQYu+5EeTXXudNi0Dm4LHDf+uej8IWzJ2GghcsYOzzp76J6tGcIIk3T7SbysBjuc/oho6INoA7i9V88+zBlHd33RUQoULwBe2319shaqVoDFfGuLpiyv5P+p1r0P+JgCLzHYVtQXv3mOWhsIKM9ftZN+T8ds7/2nGsWK7m+CnxXc+PFxxnHX20/IiKfCUN329pikV41ym353eXp2KjVOYT31VAIguYN1ijAxbG5GBcMHaD/PLPm5XAXSeo3KS7zpaFYSFjl6DBQsrzbT+McyWG6X1dnuW6HsDiJCLWa+5U3heLYeeA+APinn/0jbqtKsZL4/q46YyImZYhIXCn6HXHgMMVHt71G2XEla/sqmBbakaFcXE4aIuO1a8+17u74riaP4IKumyHmwDm8uCxOg+qyCLnvnPOjDnprpaC8u8ycuB5XuIRKsSEmDm71vzLjcO25Ic0NmHH0frj0ZDVnlAXiXFhy3XlVbVRcuDAAmBcaBJiiGnPOwcR9i/398dMn4/gDo3xm4MmCZFQ7L8RPp5ogi5n6rEsvvQgmP4foVA7P/r8+nJ/D9/EHjMCBIwemMqPJGNzcYNz5AbH1RB6LpqsBSMmBcxCRB+HiyPMz0ofYiI9dY0yZ8O/vPxKT9hmci8d+VrGSKUdGoUD45SVTccrE0VWPS1wsq+VgeVvTDhphrMfNbE8/XB/RmW86xhi4H3EeNBYLqVDscshuoPp5IYql8gjJ74J+xzmYMsFx5LFwHqRwxKkUwwY24tmvnJFLW3FYjDzaytan6wfisstSfYAq1CMOvthFtYpHAPjkew7BJx2irLogq6VbHh79LsgzLhBfLE1OliI+YzBB58/LpMBP9C09L55JTobLhsfksDhYGE+93lG/4xxckMd64rr7vfsTJ+HGi9X5hGsBvqjmIf7ghPbDx6vDRcT1gv+uk9qFiJhCFl8yPRZ71GOTlYwx1bs+Kf7IXZ9DvRaePLlDVw8Gnidk/IiB2jqc0NgU6hzyd9TYUEgl9wHMBiccJqda0fFO52eRN/oh5xD8N03NPPK1EhGOHT8MH51uls/mwaJnAdfJ5SJWQra23ImDA+dgIA4HjdJ//LWAONo8Fz0X6MKRc1BGhXS9RBaunMN7DhuDZ98yZ31ULcYq3Hjxcbj+X8zOeXxYrhyg/L4HNBaVDmou4TtMi74YXt+F0OSB/kccDJngOIY7KJxd8MBn351LO3minKvOgYcicYPrwlmtWEmkG3nc56yvnWn0SBX7aLM4J51+2BhlqIhK8NSXT8eIgXrnNiD2nrYZM3DUS9npyrn+5tJpVpv+Y0Lv44+cdKCxXrFAKBYsWRSzElOJy21pLKCts5RSSle7oA9xSPKVN/ofcXDgHP7z/LpED+9R5CE+cAnOJ8K1nsuu0qQMP3pcnJoyj43w2GHmnN1iH3KMfxm69JSV4NAxg611uFltVqfMWsNV/9TUYI9BdPDoQXj7e+/PRYHPvwtX4pXiHJqKaOsqYe7KbYly0z388MJjrMTDxTs9b/QuAWkdYIytFGK4Y96EvRE/uPCY/9/e3cfYVZRxHP/+tlsKbSmlpSB9o4VgQykvpStQbMWQgBRftoiYEqQVNECAgCYmgmJMACMYJbzUCA22Ad8wRNAiKEEjIIoClfJSaimtRgoNL0phS6Ww5vGPMzd72ft27u3ZM+ewzye52bOz584+++zZnTtz58zwuWOnc8zMXX/1Wnk1k3YMNHXPocl5lfnxc5ssRX3cQQNDdXm8IR1j8/e0Kiv5tlqLJ4vpqe2opCztHiWt68vmd9Buz2Gw0SO72fFOP5te3f6e8mYvxk7vmcYnj2h8n1Or5w+VYddzWDTnA3xwv7FNp3zm1bWOYcr4Pbhq8WGZ1HXKnP25+tP9nNriDemKtHlt9qbus1ecnKqO+QdO5JHN/46xL3uh7DN2N5bOP4DP9kxret4VvXNy7TGP32Mkh04ex1dOanxfRQxplnqpNvgV/9v9/+OpLW+8576drFz+8UM4Ioc7oyuGXeNwwMQxDaeZnj5vKnes2ZJzROXV1SWWHN18nBeS93C27Xg3h4gG3LKshxe3/TeXhr76e/wxoynHWZFUyGHS7hFd3HPxwthh1Ki8aNzZoqd15wXH8Vrfzppec71hxbsuOC6T2L64MJvpzWkVpnGQdDJwPTACuMXMrs47hmtOO5wrFxfvD6ns7r5oQcux+POOP5CbH9zc1hr3zYwZ1d30hq4sVTc/lR2+XDlV1nzae0zz4a7qu6IbGdXdxYarFmUSVwyFeM9B0gjg+8AiYDZwhqTWi8dnrKtLuU0TG06mTRjdcnG7yiJneSwolrXqaybLfchd/g6aNJYfnHkU3z718F2ua0HO09SzVpSew9HA82a2GUDS7UAv8GzUqFxudh85gjvOn/+e3bLKYu+Mpj67Ylh0WLpVeutZfORkfrn2JQBmTx7X4uxiK0rjMAV4oerzLcAxkWJxkXxoRjbz//MmiSt7D2XOlHRLvbv3r+uWzOW6JXPZ+HIfM/bJbgmdGIrSONTri9fMZJd0LnAuwPTprd8IdS4vZ82fETsEVyB5LKk91ArxngNJT6F6rt1U4KXBJ5nZCjPrMbOeSZMar6zonHNu1xSlcXgMOFjSTEm7AUuA1ZFjcs65YasQw0pm1i/pIuA+kqmsK81sXeSwnHNu2CpE4wBgZvcC98aOwznnXHGGlZxzzhWINw7OOedqeOPgnHOuhjcOzjnnaqjZXrxFJqkP2NDgy9OBf7WoYi+g+Wpw6c/Lsq40sceIK21dnvvs4ypq7tOeV9T4h+u1M8vMWt+lZ2alfACPN/naqymevyLl92l5XsZ1tYw9Ulxp6/LcD5Pclz3+4XrtNPvfWf14vw4rbWt9CnenrCvNeVnWlSb2tHXF+Bk99+nPybquvHOf9ryixu/XThNlHlZ63Mx62v1a0ZU5dih3/GWOHTz+mMoUe9pYy9xzWNHh14quzLFDueMvc+zg8cdUpthTxVranoNzzrmhU+aeg3POuSFSisZB0kpJr0h6pqrsCEmPSHpa0t2SxoXy3SStCuVPSvpo1XPmhfLnJd2gHPZ0zDD2ByRtkLQ2PPYd6tjD950m6Q+S1ktaJ+mSUD5B0v2SNoaPe1c957KQ4w2SPlZVnmv+M4499/y3G7+kieH87ZKWD6orxrWfZfy55r+D2E+UtCbkeI2kE6rqyj33mUgzpSn2A/gIcBTwTFXZY8Dx4fgc4MpwfCGwKhzvC6wBusLnjwLzSTYX+g2wqESxPwD0RMj9/sBR4XhP4DmSfb6/A1wayi8FrgnHs4EngVHATGATMCJG/jOOPff8dxD/GGABcD6wfFBdMa79LOPPNf8dxD4XmByO5wAvxsx9Fo9S9BzM7CHgP4OKZwEPheP7gdPC8Wzg9+F5r5BMMeuRtD8wzsweseQ3dhuwuAyxD3WMzZjZVjP7WzjuA9aTbOvaC9waTruVgVz2Areb2U4z+wfwPHB0jPxnFftQxthMu/Gb2Vtm9jDwdnU9Ea/9TOKPoYPYnzCzygZl64DdJY2KlfsslKJxaOAZ4FPh+HQGdpJ7EuiV1C1pJjAvfG0KyY5zFVtCWQztxl6xKnSpvxGjayppBskrpL8C+5nZVkj+kEh6OlB/P/ApRM7/LsZeES3/KeNvJPq1v4vxV0TJfwexnwY8YWY7KUDuO1XmxuEc4EJJa0i6fe+E8pUkv4DHgeuAPwP9pNynOiftxg5wppkdBiwMj7PyDFjSWOAXwJfM7M1mp9YpsyblQy6D2CFi/tuIv2EVdcpyu/YziB8i5b/d2CUdClwDnFcpqnNaKaaIlrZxMLO/m9lJZjYP+BnJ+DBm1m9mXzazI82sFxgPbCT5pzu1qoq6+1TnoYPYMbMXw8c+4KfkONwhaSTJH8hPzOzOUPxy6DJXhi1eCeWN9gOPkv+MYo+W/zbjbyTatZ9R/FHy327skqYCdwFLzWxTKC7M/512lbZxqMxWkNQFXA7cFD4fLWlMOD4R6DezZ0MXsE/SsaFLuhT4VRliD8NM+4TykcAnSIam8ohVwA+B9WZ2bdWXVgPLwvEyBnK5GlgSxltnAgcDj8bIf1axx8p/B/HXFevazyr+GPlvN3ZJ44F7gMvM7E+Vk4v0f6dtsd4Jb+dB8up6K/AuSUv8BeASkhkEzwFXM3BD3wyS1VrXA78DDqiqp4fkotoELK88p+ixk8ziWAM8RfJm1/WEWTQ5xL+ApBv8FLA2PE4BJpK8eb4xfJxQ9ZyvhxxvoGpmRt75zyr2WPnvMP5/kkyA2B6ut9kRr/1M4o+R/3ZjJ3mR91bVuWuBfWPlPouH3yHtnHOuRmmHlZxzzg0dbxycc87V8MbBOedcDW8cnHPO1fDGwTnnXA1vHJwbApLOl7S0jfNnqGrlXudi644dgHPvN5K6zeym2HE4tyu8cXCujrDY2m9JFlubS3LD4lLgEOBaYCzwGvB5M9sq6QGStbA+DKyWtCew3cy+K+lIkrvgR5PcCHWOmb0uaR7Jelo7gIfz++mca82HlZxrbBawwswOB94k2W/jRuAzlqyLtRL4VtX5483seDP73qB6bgO+Gup5GvhmKF8FXGxm84fyh3CuE95zcK6xF2xgnZwfA18j2cjl/rBi9AiSpVEqfj64Akl7kTQaD4aiW4E76pT/CFiU/Y/gXGe8cXCuscFry/QB65q80n+rjbpVp37nCsOHlZxrbLqkSkNwBvAXYFKlTNLIsH5/Q2b2BvC6pIWh6CzgQTPbBrwhaUEoPzP78J3rnPccnGtsPbBM0s0kq3DeCNwH3BCGhbpJNmVa16KeZcBNkkYDm4GzQ/nZwEpJO0K9zhWGr8rqXB1httKvzWxO5FCci8KHlZxzztXwnoNzzrka3nNwzjlXwxsH55xzNbxxcM45V8MbB+ecczW8cXDOOVfDGwfnnHM1/g/MfPZoFzeh3AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sorted_data['inc'].plot()" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 19, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEKCAYAAAAMzhLIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXl4Y2d977/v0dE52m1JlpexPfbsk5msM5OdQMhCAqUJpdCGUghc2rTctOW2vb0PdLnltjdlaQstDwUalktSKJBCWmgTCMmEkJBMFs9MkpnJ7Is9Hnu8yLK1S0fSe/845z06kmUttna9n+fxY/lIOj6Sjs73/e2EUgoOh8PhcFZCaPQBcDgcDqe54ULB4XA4nKJwoeBwOBxOUbhQcDgcDqcoXCg4HA6HUxQuFBwOh8MpChcKDofD4RSFCwWHw+FwisKFgsPhcDhFERt9AKulp6eHjo6ONvowOBwOp6XYv3//PKXUV8lzWlYoRkdHMTY21ujD4HA4nJaCEDJe6XO464nD4XA4ReFCweFwOJyicKHgcDgcTlG4UHA4HA6nKFwoOBwOh1MULhQcDofDKQoXCg6Hw+EUhQsFh8PhGNh32o9Ts6FGH0ZTwYWCw+FwDHz80dfxhb2nGn0YTQUXCg6HwzEQSaQQjCuNPoymggsFh8PhGIgrGYTjqUYfRlPBhYLD4XAMxJU0wgkuFEa4UHA4HI6Gks4glaEIcYsiBy4UHA6HoxFX0gDALYo8uFBwOByORlzJAFCFglLa4KNpHrhQcDgcjgazKNIZiph2m8OFgsPhcHTiBnHgmU9ZuFBwOByOBnM9AUAoL07x1//1Bu57eAxKOpP/tLaHCwWHw+FoxFMrWxT7xwP46Rsz+Iv/ONxx8QsuFBwOh6OR43rKsygiiRQkk4DvvnIeY+OBeh9aQ+FCweFwOBqxZFYo8msposk0Nvc6AAD+cKKux9VouFBwOByORjyVjT/kWxThRApehwQASKQ6K07BhYLD4XA0crOechsDRpMpeOzSssd1AlwoOBwOR8MoAEbXUyKVhpKmBqHgFgWHw+F0JCsFsyMJdXuPQwagCkcnwYViFUwGojg0udTow+BwOFWGWQrdNnNOHUVEu80tCk7Z/ONTJ/EH3z3Y6MPgcDhVJq6kIQoE3VZzTh1FJKne7rKaYRIItyg4pQnFU1iK8QlYHE67EVPSsJhNcFjEgq4nm2SCRRR0i+L8QrQjiu+4UKyCeCqNaJL3geFw2o24klGFQhZzLQpNNByyCNlsQiKVxvRSDG/525/hmRNzjTrcusGFYhXElTTiSgaZTPuvJDicTiKhpGExC3DIuTEKtjC0SaJuUcyHkshQYC7U/sV3XChWASu2iXeYn5LDaXfiKdX15LKICCey7uWw5npiFkVcSettyDuh+I4LxSpg/slokgsFh9NOxJJpWLUYhbGOQrcoZBNkUUAildG3JTqg+I4LxSpgJ0aMCwWH01aoMQpBj1GwQHXYEKOwMIsiyS0KThFYUQ6fgMXhtBfM9eSwiEhlqC4C0UQaAgFkUdAtCu564hRFP3m4RcHhtBWxZBqyaIJTFgFk23iEEynYZRGEEFjMJiSUtP7974SaipJCQQgZJoT8jBBylBByhBDyMW27hxDyJCHkpPbbbXjOJwghpwghxwkhdxi27yaEHNLu+wIhhGjbZULI97TtLxFCRqv/UqtHnLueOJy2YP/4Au7/1wNIaxmMiVQGVkm1KAAgpDUGjCRScGjioVsUTCg6oEq7HIsiBeCPKaWXALgOwP2EkB0APg5gL6V0C4C92t/Q7rsHwE4AdwL4EiHEpO3rywDuA7BF+7lT2/4RAAFK6WYAnwfwmSq8tprBWhHHFF5LweG0Mt96cQKPvT6tz5eIK2lYREHv6TSrpb5Gk2nYJPUyxmIUUR6jyEIpnaaUHtBuhwAcBTAI4G4AD2kPewjAu7TbdwP4LqU0QSk9C+AUgGsIIQMAXJTSfVSNED2c9xy2r+8DuJVZG81GKp3RVx/c9cThtC6ZDMVzJ9ViuUWt00Jcq8wedtsAqJXXgOp6WmZRsBhFB8QqK4pRaC6hqwC8BKCPUjoNqGICoFd72CCA84anTWrbBrXb+dtznkMpTQFYAuCt5NjqhXGwCRcKDqd1eWM6iPlwEgAQiKi/Y0oaVsmEdd1WEJIVimgyBZukCkU260lLj+UWRRZCiAPADwD8D0ppsNhDC2yjRbYXe07+MdxHCBkjhIzNzTWmbN7YhrjThpdwOO3Ecyfn9duBqAJKqZoeKwqQRAEDLgvOB2IA1II7u8yEQq3M5sHsPAghZqgi8W1K6aPa5hnNnQTt96y2fRLAsOHpQwCmtO1DBbbnPIcQIgLoArCQfxyU0gcppXsopXt8Pl85h151Etyi4HDagmdPzMGrtQ1fjCb177ZsVmMRwx5bjkVhl9Xtsqj2eory9NgsWqzg6wCOUko/Z7jrRwDu1W7fC+CHhu33aJlMG6AGrV/W3FMhQsh12j4/mPcctq/3AHiaNmlLRqMVwYWCw2lNMhmK/eMBvG1nPwDVomDZSxajUARUoYho6bHq/QIyNJs62wlZT2IZj7kRwAcAHCKEvKpt+1MAnwbwCCHkIwAmALwXACilRwghjwB4A2rG1P2UUnZF/SiAbwKwAvix9gOoQvQvhJBTUC2Je9b4umoGdz1xOK3PQjSJZDqDbX0OSCYBi7GkHpy2MqFw2zATTCCupBFJpGGXshYFkI1rdILrqaRQUEp/gcIxBAC4dYXnPADggQLbxwBcWmB7HJrQNDvGyVa81TiH05qwjq+9Lgu6bWYsRhR94Wcxq46WYY8VgBrQjim5MQoACESZULS/RcErsyvEuHoo1/V0YCKAd/3T89wC4XCaBFYf0euU4bZJCESTejdo5npa71FTZI9dDAGAIT1WvX8xqqbUcqHgLMPojyz3wj92bgGvnl/EfLj9+9ZzOK3AbDAOAOh1ahZFVNErra2GGAUAHNeEgqXHyppFwRoF8joKzjKYOEgmoWyLItBBKw8OpxXQLQqXrApFLKm7lZkQ+BwyJFHA0Wm1GoBlPTGLg9EJ32suFBXCTgq33Vy+ULCgl8EaOTS5hNNz4eofIIfDKclcKAGnRW0ZrrqelGWuJ0EgWO+xYWw8AACwS9nKbCNcKDjLYBaF2yaV3RQwG/RSH5/OUHz4m6/gU48fq81BcjicosyG4uh1qv2cum0SFqNJxLXvs0XMWgz/821b9e98Npida1F0QuyRC0WFsJOi22Yuex5Fvuvp5bMLmA8nMBeK1+YgORxOUWaDCfQ6LQAAt80MJU3h1yx/q5QVgjsvHcAPPnoD3n3VIHascwHItSgsZgGpDEUq3d5WBReKCmG9njz2CiyKSG4a3eOHpgGoudwcDqf+zIYS8OkWhRkAcHY+AiAbi2BcOtiFz/36leiyqo8zWhTdVrWyO8mFgmOExRm6bVLZdRTMokim1M6zPz58EQCwEOZCweHUG0rpMtcTAOw9OoOBLgt8WovxlTBaFExk2r06mwtFhcRTaZhNBA5ZLMv1RCnFoiFGsX88gPlwAjsGXIgk0x3h3+RwmolQIoW4kkGvSxUEtyYU5/xRXLvBg1ITDowWBXtuuwe0uVBUSFxRRyVazSbElQwymeItqUKJFFJsepaSwaTWO+amrT0AsoFuDodTH2aDrNguG6NgXLOh9HSDHNcTsyjavI0HF4oKSaQysJgFfdpVKatiMaLkPJflag91q+0BFiJcKDicejJnqMoGsq4nALhmg6fk83NdT9yi4BSAWRTlCoXRYkik0vrjB7q4UHA4jWBWyzZkricWpO5xSNjks5d8fifGKMrpHssxkFBUi4KZn6Uyn3KFIjtGdR23KDichsAsCp/mepJEAU6LiGvKiE8AgGgSIAoEqQzV3Vbt7nriQlEhiRSzKNS3rlR1do5QKBmkMxkQAvR3qScpFwoOp76cnovAKYtwWbKXv3/49Suxoae0NcGwmE0IJ1Id43riQlEhcSU3RrFSiuyEP4onjlyESciuUBKpNJR0BhbRhG6rGQLhQsHh1JsD4wFcNeLOsR5uvaSvon3IooBwAui2doZFwWMUFRJX0rCYTVnX0woxin/bfx4PPH4UY+MLIETtSJnUgtlWyQRBIHDbJC4UHE4dWYopODEbwu717jXth33/Pdoo1TiPUXCMxFNpuKzmbDB7BdfThDZr92fH5tBlNUMgBIlUBjElDYsWDHPbuVBwOPXk1fOLoBTYM7o2oWAdZjslPZYLhUY6Q/Ha5CJ2rbDS+OxPjuHCYgwJJQNZNLqeCp8gbCh7TEljoMuCmJJGIqUW2Fmk7GrEz4WCw6kb+8cDEAhwxXD3mvYjiyZIJkGPVbZ71hN3PWk8fWwW7/7SCxj3Rwref3BiEU8fm0U8pbqerCUsivOBmH6722aGLAp6HQXrTum1S3ofKA6HU3sOjAewvd+lT6tbLRazAKtk0lNl2z2YzYVCYymmFsattMIPJRSE4ilcXIrDYhbQo/WDmV5a3gE2lkxjLpTQsyjcNgmyaEJCyWgxDu564nDqTTpDcXAisGa3EwDdqyBrsYp2dz1xodBIaiuCSKJwFlM4rm5X0hSyqAaz13VZClogrE3He/cMAVCrN2WzoLuemDXitauzeku1AeFwOGtnNhRHJJnGtn7nmvdlMattfHSLgrueOgNFaxPMBCGfkGE7C2SN9thxtoBQnNeE4toNXrzz8gG8eWuP7npSg9mqULhtEjI0a81wOJzaMbWouoNZsetauHrUg+s3eSEKBAJpf9cTD2ZrMIsivIJFYRQKdqEf8drxk8PT+vYXz/jxn69NYUuvAwAw7LHii7+xCwDwyNh5xHXXk2ZRONTUuoVoEm57tt8Mh8OpPhcWVTfxYBWE4v63btZvy6Kp7V1PXCg02OCRQkIRV9I5g0nYhX5Djw2BqIKlqIIumxlPvTGDb780gcuHumAxCzl97WXRhKWYohXsqc9nVZ2LvIMsh1NzpjWLYkDrilAtVLdye1sU3PWkUSxGwcSDFXIyv+SoVw1Wn9PcT8G46kJ6fXIJw25bTuWnLApawV02mG3pEP8mh9MMTC3G4LSIcFrMpR9cAbIotP13mAuFRtaiWG5CMrfT1l41CMYsgtGeXKEwxhqGPbacfWTTY9Owas9nGRPxNjdbOZxm4MJivCpup3xk0dT232EuFBrZGMXywHJIsxSu1Ip0mEWw3mMDIdlZu8FYSq/UHHbnnpCyaEJcUduMM6Fh+2n38n8OpxmYXopVJZCdTydYFDxGocGyniJFLIobt/Tg8UPTuiWhpshacW4+63ravd6NXSNu3LzNl7MP2SwgkkgjQ2FwPXVGDjaH0wxMLcb0xV41Yanv7QwXCo1iWU9MKDb57Hj9k2/LiT2MeG0451fTYZdiCrb1OXMyIhhqt0l1Pxbd9cQtCg6nHsSSaQSiSk0sCoto4sHsTiFZpI6CuZ5cFvOywSYjXrveADAYU+CyFg6USYapWLrriVkUJabkcTictTG1pGY81SRGwbOeOgc966nAfAlmUTgtyw2wPpeMhUgSyVQGoURqRaGQxexAdqseo2DB7PY+yTicRjNVo9RYoDPqKLhQaOiup4IWhbrNXqCRmFerlZhYiIBS5EzNMiIXsCjYtji3KDicmlLNqux8OiGYzYVCQylScBeKK7CaTTCblr9dPVpF9Zk5NaC9skVhFAr1tiAQSKb2N1s5nEYztRjPGUFcTWRRWHGAWbtQUigIId8ghMwSQg4btn2SEHKBEPKq9vMOw32fIIScIoQcJ4TcYdi+mxBySLvvC0Rz9hNCZELI97TtLxFCRqv7EssjmS5ecFfI7QQAPU7VojijZT51rSQU5uWuJ3W7wC0KDqfGHJ0OYr3HVnCxt1Y2+RyYDMQwF0pUfd/NQjnv2jcB3Flg++cppVdqP48DACFkB4B7AOzUnvMlQgi7Kn4ZwH0Atmg/bJ8fARCglG4G8HkAn1nla1kT2RhFelk311B8ZaHw6hZFGIAa8C6E0aIwioZaX8EtCg6nVmQyFK+cW8A1o56a7P+t23sBAD8/MVeT/TcDJYWCUvosgIUy93c3gO9SShOU0rMATgG4hhAyAMBFKd1HKaUAHgbwLsNzHtJufx/ArSQ/tagOJNNZccgPaAfjChwrCACzKFjRncu6UowiKw7M9cRut3sgjMNpJKfmwghEFVyzoTZCsXOdC71OGT87PluT/TcDa7HDfo8Q8rrmmmKTQAYBnDc8ZlLbNqjdzt+e8xxKaQrAEgDvGo5rVSQNcYL8OEUonloxSO2URUgmIRujWEFQjOmxOa6nDgiEcTiN5KWz6jr32g21uawQQnDzNh+ePTGnxzrbjdUKxZcBbAJwJYBpAH+vbS9kCdAi24s9ZxmEkPsIIWOEkLG5ueqaecYPOD9OUSxGQQiB15Gdfd1lKyeYbcq5zWMUHE7tePnsAvpdFgx7qp/xxLhley9C8RT2jwdq9j8ayaqEglI6QylNU0ozAL4K4BrtrkkAw4aHDgGY0rYPFdie8xxCiAigCyu4uiilD1JK91BK9/h8vkIPWTXJVEYPROc3BgzFFTjllTtOsrGohAAOqXR6rDVPKHjWE4dTGyilePmsH9ds8Cwrlq0mu9arTpUTM6Ga/Y9Gsiqh0GIOjF8BwDKifgTgHi2TaQPUoPXLlNJpACFCyHVa/OGDAH5oeM692u33AHhai2PUlWQqA7dmDeTXUoTiKThWsCiA7AAipyxCEAqfjHKeOOjbRZ71xOHUirlwAjPBBHatr36PJyNsvHGyTRd9JXs9EUK+A+BmAD2EkEkAfwngZkLIlVBdROcA/A4AUEqPEEIeAfAGgBSA+yml7Cr4UagZVFYAP9Z+AODrAP6FEHIKqiVxTzVeWKUo6QwGui0454/mxChS6QyiyfSKricga1GsVEMB5GU95bmh2BwLDodTXZai6nfLaxgiVgtkvcFnhwoFpfR9BTZ/vcjjHwDwQIHtYwAuLbA9DuC9pY6j1iRTGXi0iXPGGAXrJlts2AmzKFaqoQCy4iCLQo7VYTHzYDaHUyvYIqzYIq4amE0EhLRv3zZema2RTGf0udVGi4KdaEUtCrtmURQREzmvv5O+vQOGnnA4jSJYpE9bNSGE6MPJ2hEuFFADXsl0Bp4CQqE3BCzQ54nR41Sft1INBQBIWkWosYaC/c0L7jgcYCYYx96jM1XdZzCW7fxca+Q2bjfOhQJAKkNBKeCQRYgCwSvnFvBLX3gOs8E4FrS0127NLVUIb1kWhfpWWwtYFO1qrnI4lfB/HzuK33p4TI8r5JNIpZEqUKdwZi6M01pnhHzYQm+lOqhqoloU7fld5kKBbA2FJAqwyyKeOT6HI1NBHDy/iAuL6qyJIffKOdgsmF1OjGKZ68ks8DbjnI7HH07gicMXQSlw4HzhWoT3f/UlfPaJ48u2/+WPjuATjx4q+Jx6xSgAbS5Fm3oHuFAgm9ImmQQ4DC6m8wtRXAjEIJToOtnjYK6nlU/GrOspVygsognJVGZZfykOp5N49MAFJNMZEAIcMBStfevFcbx2fhEAcD4QxYQ2TdLIYlTBtDaYKJ9gLAXJJORkGtYK7npqc3ShEAV028wY9lhhl0yYDMQwuRhDn8tStOtkj0PG+65Zj1u05mCFYMGu/BgFc0kl26T0n/W84nDKhVKK77w8gT0jbuwYcOHAhCoUs8E4/uKHh/G9MbUrUCKVQbSAmzaaTGEulECh8qtQXIHTIta02I7BXU9tCqUU/nBCv0hLJgGffvfl+H8fuhrDHhsmA6pFUcztBKhzJT717stw6WBX0cepQrHcogDaY3jR/vEFvPXvnmnb6tRacXY+gv/9w8MFZ6F0AotRBWfmI7jz0n7sHnHj1YlFpNIZPHFEdUWx70ZCySCeXP49iSXTiCvqhMl8gvGVp05WG5711KY8euACrv/005gPqwFrSRRw2VAXNvc6MeS24fxCDJOBWNXm7Mpm07JgNhOOdjjBphbjANDWffmrzcGJAN79pefx8L6si6XT8EfU86XXZcHuETciyTSOz4Tw48MXAagCQSlFPJVGVFkuBszKKHTeMYuiHqiJKa3/PS5ERwvF3mMzSKYymNbGJBrdS8MeKyYWorgYjGOwhEVRLt1W87LsqXYah8pWxLECqz5OYf7s3w8jpcWnWIZdp8EWal67pPdMeuiFc3jxjB+A+t1Q0mpmYqFzK5rMFYpMhuKrz55BMK4gGFPqkhoLaMHsNnU91Udqm5BMhuLFM2rvQdb51dgKfNht08cbDrltVfmfX/3gnmWrG2ZRtEMtRUjLMGn3sZDVIp2hODUXxjsvH8CjBy4gEO1MofAzoXBIGHJbcdslfXhkTJ1K4LKISKQyelFqvlCkM1SPMTKhODEbwgOPH4XLKiIUT6HPVf3xp4VoZ9dTxwrF8ZmQvoJbKCAUxrhEtVxPoz32ZdtYcLsdViKsmSIXivK4EIghmcpg13o3Hj1woWMtigXN9eS1yyCE4Ksf3I0XzyzgzHwYj70+jbiS1l06+edW1DBkjAkFq8OYWowjGK+jRcGzntqPfaf9+m32BTWbspkRw56sFVEqmL0WZLF9LArWLqEd3Gj14My8WiS2tc+JLqu5Y4ViPpwEIdC7NxNCcP0mL95/7Yg6ryWV1s+paJ5FYbQw5sKaUMSYUMSKjjGuNuoQsvY89ztWKF447dcDy+wLasy1NgrFuipZFIVgFkU7XFx5jKIy2FTEjT47vHapY4XCH0mg22qGWCAFnU2AZCv1RF7NkVE4mEXBFiznA1FEk+n6ZT2Z29f11JFCwYaZ3LSlB4DB9WTKZiQ5ZBFumxk9DnlZSms1aaesJ+56qowz82G4LCK8dgluu9TRMYqV2oDnWxRA7vllFIpZJhSaRXFiRrXY6pr11Abf40J0pFBMBmIIxlO4cXOuUJjF3KKc9V57TccnAu2V9RRK8GB2JZyZi2CDzwFCCNw2SQ/qdhr+SBJee+FeaqxppvECHMsRDXVxIolCNkahCQX7XtcvRsGzntqKI1NBAMAVw92QRcFgUeTq5l/fvbPmx5LNemr9E4xZFIWKojjLOTMXwQ2bvADU1NBDFzq0jiKcwPZ+V8H7ZNGkBbMN4pBcblGs99gMrqfcpoL1tCiUNEUmQ1ecdNmqdKRQvDEdhECAbX1OOC2iIZidKxSXD9V2fCKQtSjawWQNcddT2UQSKVwMxrHRp2bCue0SAhEFlNK6tJtoJvyRpN7iPx/m91/JomBCMeKx4czcLNIZimAstyivXjEKljWZTGdgEWrnrm4EHel6emNqCRt9DlglExyyqLfwqEfjsHzkNrIoWAuFWBtkcNUa1hNro88BAPDYzUimM4h0mDWmpDNYjCr6lMh8WNPM2AoWBbu93mtDhqqB8aWYAqPW1tP1BKAtq7M7VCiC2LlONXUdBrNUaoBQZOsoWv/k0gvuOuxiVymZDMU3nj8LANjW7wQAuLWK/YUOi1OwAH6xYDaQDVADuQFso0UBqJlPwbiCDd5szVLdXE9tVBOVT8cJRSCSxNRSHDsGNKEwtBUv1iG2VkgmoS1m7SrpjF4L0g7W0Vo4Ox/B5588gf98bSqnIIzxwONH8eiBC/gft23BJs2iYCvqhRpmPgUiSfzRI6/qgt4MsAB+z0quJ23xtmQQiniO60l9f0e0YtbZYALBmIKNPrteF1W/poDtk8GYT8cJxdFpNZC9g1kUcvYkaoRFwdqPt/rwooihc2enxyi+9eI4/nHvSfz+dw7iKz8/s+z+/3p9Crfv6MPHbt2ib2MWRaCGtRSvnFvAowcu4PXJpZr9j0rJtu8oblEsrWBRMOt1U48quNNLcQRjCrqskt66w1FkjHE1ycYb2+/877hgNst4YhYFM0sJAcQGZSpYzKaWX4WzQDbAXU+TgSg2+uxQ0hmcKzCfIxRPYb3HlhO0ZsFcfw2FghWiLa4warQRsM6xKwWzmWt20SAUOcFsJQ2ziWBdtwUCAaaX1NT3LqsZ67qsWIoqMNXpe51NdW/tRV8hOk4o3rrdB7ss6isYttowm4SGZZuw6tNWhgmFZBJaXvTWymQghhGPDcl0BpOB3IlsqXQG0WR62SqXXShraVEwP79xdd5oWOfYnpWC2QUsilgyd1FiNZsgmgT0uSw4vxBFOJGCyypiQ48d85H6tbyX26h4Np+OE4rNvU5s7nXqf7NgttyA+ASDVZ+2Mszv7XPKHe96mgzEsGu9G4lUGs8cn8u5j7U5yQ+wOmQRZhMpGKO4uBTHQiSpu0tXC6svaCahGPdHIApkxcwktkoPrmRRJFOwSep7OdBlwXGtGttlMeODbx+tazymnV1PHRejyEe3KBoQn2BY2mDgCbsA9nS4UATjCpZiCobcVgy5bZgNJXIsLGZ55V8YWXV2oaynv/vpcfzOt8bWfmxafUEzCAWlFH/+H4fw8L5xvHmrb8UCNaNFwdxQ+VlPNkl9zEC3FadnVaHosprhsUsY8S7v2Fwr2qkmKp+OFwq2ssuvyq4nslloA4tCvQj1OuWOjlFcCKhDsIbcNr3r8JQ2GAvIvk+FUjY9dmlFi2JqMZ7TDG81hHSLovEpuBeDcXzrxQm8d/cQ/vkDu1d8nB6jiCqwSyIkUVhWU2FlQuGy6DVR9cp0MqJnPbX4oq8QHed6yodZFI3IeGJYzaamCjCuBlZs1+uU9Q6f7dbGoBwmdaGw6ivLyUBML6xjF2tnAVfLkNuGgxMBLMUUdBkudPPhBNIZCn8kiXP+CL7yzGkEokn86TsuwZ5RT9nH1kyuJ2bd3Lytt2haOrv4LsUU2CUTUhlTTouYfIuC4apT7YQRXkfRxmSD2Y27qO0ecePQhSV94Eorwvo8+ZxqkkCrW0irhQWvVdeTVdtWnkXxsVu3YCGSxN//9HjOdtZiZjYUx3densBzp+ZxYGIRL51dqOjYmsn1lBXM4hd0ZlEE4wpkswk2yZTrelLSsGoxinVd2Ul2XbZGWBTc9dS2sGC2JDauN8vtO/qQzlA8c2K2YcewVkJxBaJA9HqATnU/TQZisJpN8NjVPH5RIDmZT6zDbqEL5GVDXfjg9aP4lxfHcWo2BECt4taFIpjAhUAMlw12AainqdtFAAAgAElEQVT8gtRMFgWzQB0lhIJZFJSqF2Kr2ZTnekrBZi5kUTTQ9cSFov1wagV3jXQ9XTHUjR6HjCffmGnYMayVcCIFh0XU/cWdGtCeDEQx5LaCEAKTQDDQbcGFgjGKwhey91+7HpRm632CcQUpLTYxE4zjwmIMQ24rpFW0tG4qoWDvQ4liOOMsGNlsgsVsQiyZxr+8OI4fvTaV43oyWhQNiVEw11Mbnvs8RqEHsxvnehIEgtsu6cVjr08jmco0VLRWSyiegkMW9amBnVpLMRmI5YzOHeq2le16AoBep3qxYy2zjQV4U0txXFyKY7DbCtkkIFmpRRFrnoK7cAnBZLCLLwBYRAFmgSCmpPGVZ07D65BygtlehwxRIKAA7FL9PQTc9dTGNEMwGwBuvaQPoUQKByYCDT2O1aLOJjbrQhFLtt+XpRxUocidt57jeoqnYDaRFTsVu6wiJJOgz382DjM6fGEJqQzFoGZRVCIUlFKE4mpX1VA8hfQaM6jWSlhzwZVyPVnEXIvCKqmJH1NLMZyYCSGSTOkWhUkg6HNZ4LKIDSmeZZmTXCjaELaya0RDQCNb+9SsGOPqs5UIJxQ45c52PUWTKSzFFAx0Z10gg24rZoIJ3U0UiitwWswrXsgIIfA5ZcyHVIHwa4IhCgQHtUXEkNtWsVBEkmlkKNCnWSyNbgwYiqdASOmVv9lEwJLnLFqM4vRcGJSqrTLiSkYPZgPAum5LTsZYPWF923jWUxsiiwJEgTS0jgLIuhxmgvGGHsdqCcZScFpE3afciUJRqJiOtYph7h7V8iq+iu5xylmLQnM9be51IKDtY7BbsyjS5QsFq2xmo30bHadgrspSK3/14queU8yiyF+x2wxic9cV6/DLV6yr/gGXSTu04ylEyasjIeQbhJBZQshhwzYPIeRJQshJ7bfbcN8nCCGnCCHHCSF3GLbvJoQc0u77AtHOEEKITAj5nrb9JULIaHVfYsnXB4dFbLjrySqZ4LKImG1RoZheiqG/y2JwPaVbVvRWC6tON/ZxYqtbdqFWLYriQuFzyNkYheZ6umQg275jsNsKuUKLggWyh7W5DY2OU4TiqbIzk1iKrEUUckSBYdz2getH8cdv21adg1wFsnm5kLUD5VwdvwngzrxtHwewl1K6BcBe7W8QQnYAuAfATu05XyKEsE/xywDuA7BF+2H7/AiAAKV0M4DPA/jMal/Mahnstuor+kbS57JgNlS/JmbVIpxIIRBVMOi26q6ng+cDuPZv9rZszGU1RIoIxVLMYFHIxS+QPqdBKCIJtROq5s7y2iVYJZOW9VSJRaEe27AWP2m0RRFOKGW3/2ZWqmwW9NsOWcSglg5rNTcutT2fjnU9UUqfBZBf2XM3gIe02w8BeJdh+3cppQlK6VkApwBcQwgZAOCilO6jlFIAD+c9h+3r+wBuJaXs0Srz7d+6Fn9yR+NWIYw+l6UlV+HGthXsSzt2ThWIs3PL22y3K8yisJcSilIWhVPGQiRbje11ZGcrDGoZVVKFWU+hPIui0UIRiqdKBrIZulCIJv38GvHa9OmANql5kjflCgW8VVitv6WPUjoNANrvXm37IIDzhsdNatsGtdv523OeQylNAVgC4C30Twkh9xFCxgghY3Nzc4Uesiq6bZK+Em4kvS4ZM8HmsSgyGVpWmuuFxWw1MvsiswFR/jq2eW40LOWzuEWhlEwJ9Tkkff6zP5xAj13WLV6WeltpMFt3PbmbI0YRTpQWTAbLELOYs66nEa8NW/uYUDT+u8uQtRnf7Ua1HfOFLAFaZHux5yzfSOmDlNI9lNI9Pp9vlYfYvPQ6LZgNxaEaXY3nX14cx81/+0zJ49H7G3VbYZFyO3z6O2gGdCTJLIrshWu1FgWg1lL4w0l47BJ6Xeo25m6RRBMSFQWzNdfTChbFqdnQmpsOVgILZpeDbLAomHWx3mPHtn41U7AZFnkM2cwtCiMzmjsJ2m/We2ISwLDhcUMAprTtQwW25zyHECIC6MJyV1dH0OeSoaSpnt3SaI5dDOJiMF7yxL8QiEESBfQ4ZEgmAcZegCx7pxMIJ1RxNF4AWXO6pZiCTIYinKxMKBY01xMTCNY2u+Jgdiw7L0QWhRyhOH4xhNs+9yz2HqtfCxlWd1MOlhyLQn3vRr023HpJHz7ypg24cri7ZsdZKWrWUwfGKFbgRwDu1W7fC+CHhu33aJlMG6AGrV/W3FMhQsh1Wvzhg3nPYft6D4CnabMsqesM80M3S5xiVnODGcecFmIyEMNQtxWCQEAIyQkudpRFUaB/kWgS4JBFLMUURJIpUFq6EZ7PkT0PFqJJeB0y+lwWfO++6/Ce3ep6q9IWHsG4AqvZBLNJQJfVnNOA8rmTqhv3Yh3Pu3KyvxhGi8KqWazrvTa4LGb8xTt35LT5aDSy2KFZT4SQ7wDYB2AbIWSSEPIRAJ8GcDsh5CSA27W/QSk9AuARAG8A+AmA+yml7Gz+KICvQQ1wnwbwY2371wF4CSGnAPwRtAyqTqRPcy80i1Awa4AFaVdiMhDVg6xAriugk2IUkUQKAlmehdNlNWMpppTs88TocaqNFU/MqIVlXm1M6rUbvdnAboXB7GBMHQ9qPB7GC6f9AKpfhPfN58/ia8+dWbY9mcogkcqU7PPEMFoUV4968I7L+nHFUPNYEUbaNZhd8pOilL5vhbtuXeHxDwB4oMD2MQCXFtgeB/DeUsfRCbCA5WyTBLTZcYRXsCjG/RG4LGZcWIzljOlkF7ONPXa9wrgTCCdSsEvLi8hcVjOCOUJR/Gtnk0Q4ZBF7j6pNItcZuqIyKg1mhxKKXrfQbTNjURtelEpn8LLWrpzFMarFf70+jYVoEr9108ac7eEyO8cyjFlPQ24bvvT+lQcdNZrVNGtsBZonr4yjByybwaLIZCjmNYuCtcY2QinF+x58EW67hPlwMqe/EVtRXz3qwaMHJ0EpbUjvnXoTjqdyUmMZXVZRsyhWHlqUj88p4+x8BJt8drx12/LEjcors1N6R1WPXcLZeTVt+fBUUL9wB6tsUSRSGUwuxJDOUJgMgatyGwIy9II7c/M3kpDbYKxxIZr/ne8gZNEEt82MmVDjhWIxlm1vXciiGPdHMbUU19thD3bnup66bWZs6XNASVMES8Q42oVIMpWT8cRwWfJdT6XXZz6t9cfH334JxALtZVZTmc3+r9ch67GjfZrbyW0z6wHvahFX0kimM5heyu1fxgSp7Kwn0ZTzu5nhWU+cuqAW3TXe9TRrEKtCMYqXz6nuivVauqWxtbbHLmF7vxM92sVuvkMyn8KJNBwFVsksJsAukOWM6bxxcw9+6bIB3HZJb8H7K24KqLnFAKDHIWMhmkQ6Q3FgIoBNPjvWe+1VF3R2wRz3q3U2s6E4Pvqt/Th0YQlAeYIJZC2JlTruNhM2swnRZPstjLjrqcnobZLq7DlDK5FCQvHK2QW4bWZ87d49+Me9J3NiFJ9692UAgFOzYQBq5tOm9it7WUYkkYKjgEXBhCKgNfjrskol9/Wx27YUvV8ymZDK0GVunZUwzm3ocUigVB2xOr0Uw7DHhnSGVj2YzYo1x/1RXLMhg/u/fQCvnAvog5zKF4psU8Bmx22XEE2mEVfSTZWNtVaaX6I7jGG3FRML0dIPrDHGgHqh9NhXzi1gz6gHW/uc+Kff2JXTRmGgy4qBLiu8dtWi8HeIRWFctRvpspoRVzI4djEEpyyix1FaKErBmliWa1VElewkOP1ziSQwE0yg32WBy1Ib1xOgJj3841Mn8cq5AJyyiNcnVYui4l5PLWBRsM/WOHCqHWj+d77DGPXasRhVsBht7InGUmMJWW5RzIbiOOeP4upRd6Gn6rAvTae4nlaqNu6yqe6o/eMBbOp1VCWwX6lQ5FsUAHBxKY75cAK9LgtcVrFmrqez8xE8MnYet+/ow6/uztbdlhvMzrbwaP4VuqdNF0dcKJqM0R618vacv7FWxWwwAbtkgtsmLQtmv3RGjU9cPeopug+PnQlFe62uVkINZhe2KADg5GwYW3odVflf+tjNdOlUzHSGIpHK6NlobEbGsYshUIqaWBSUUl0onjs5j9lQAm+/tB/XbMieM5UX3DX/5crbphYFj1E0GaNeNTh8bj7S0NYEc+EEfE4ZGbrcovjJ4Yvw2iVcNthVdB+iSYDbZu6YortIonBHVJdh4trmKglFJRYFGyLFXE8so4plrPV3qd1qE6kMEql0VbKLmEgQov5/k0Dw1m29eiZdsXGw+dyyvRcTN23IyaxrVnp0i6K9hKL5JbrDGPbYQAj0PPdGMRuMo9dpgUMWc2IUobiCp47O4J2XDxRM28zH65A7ougukUpDSdPCrqcaCIVuUZQhFCwLh40MdVlFiALBkSk1VtDrtOhiVqpdS7mwWoIRLStuz4gbbrsEn1PGxh57WdPtGIPdVvzZL+2AUEbQvtF4NItioc0WR1womgyL2YR1XVaM+xsrFMyicFhEhA0Fdz89MoNEKoO7rhws8uwsPQ6pI2IUEa0hYKEZ0Eah2NLrrMr/Y6N7y7Eo4kn1MTbNhUMIgdeRLbrr77LoVdvVcj+x6mTWCvz2HX36fXde2o/t/a6Cz2t17JIJsii0nUXBXU9NyGiPDWcbHKOYCyXw5i0y4ko6pwDwh69NYchtxa715bnFht027D022/bV2SyOUyxGIYtCTk+stVCJ6ymqMIsiK2I9DnX2idlE4LFJeh+oagW045pF8eatPqz32PCru7JB7D+5Y1vbnguEEHi1bgXtBLcompARr72hFkVcSSMUT2UtCsPFY+zcAm67pK/sL/quETcWIsmGB+drTaF52QwmFBt9jrJqHspBF4oy2niw2SBGoWAB7V6nBYJA9AykalsU3TYz/vydO+C2Z1OC21UkGF6HzF1PnNqzocEpssxV5HPIcMiifhGMJFKIJtPo7yp/vvjuETWFdv94e8/OZkOLCgWzzSYBdslUtfgEkG1nUVYwWxMKm9loUagXbtZfjLmeqhWjYBZFK7TdqDZehwR/JInHD03jngf3QamgJ1ezwoWiCWl0iizzr3odEhyWbDCbCQhrzVEOm30OOC1i2wtFoXnZRv7yl3fivrwuqmuhoqwnJhSGYkD2GfZrM1CyrqfqWhSt0Miv2njsEvzhJB47NI0Xzyzg6ToOhKoVnfcptgDDHtWPPRmonVD87NgsHj0wWfC+BS0H3GOX4JRFJFIZJFMZXSi8FVQWCwLBrvVuHGhzoYgUcT0BwK9dPYzLhoqnE1cCC2YXammdyVBMGBYZUYW5nrJfdzbjgg3LqnYwm1kUrVAkV216HDL8kQQOaz2tvvvyRIOPaO1woWhCmL84UmJg0Fp4eN85fGHvyYL3GS0HduGLJFJ6gM5XgUUBqO6nE7OhZXOa24lICYui2khF0mO/f2ASt/z9M3pjx1heeiyQtSiYUNgkE0wCqbpF0QpFctXGY5cQVzIY90fhtpnx8xNzmFqMlX5iE9N5n2ILwHzJLOWyFsSU9IqZGayqVHU9qaIVTqRW5XoCVKGgFHjt/OIajri5Ye45R4FeT7VALuJ6evGMH6kMxdk5NSEiWiBGwazC/i71sySEwGkRqx6j6ESLwmsI3P/x27YhQ4H/fG2qgUe0drhQNCE2rQNpLdsVx5UMwomU3rjNiD+cgNVs0ietAarvmhXOVeJ6ArKtyJuhK26t0OsoCnSPrQVykaynVydUQWbNJQtlPe1Y58LWPgeuGs7263JZzJgNJvAn//YaTs2G1nR8nWxRGBdS77hsADbJlNONuRXhdRRNiGQSIApE/4LXAiYQc6EEhj22nPv84aQuBqwfTziuWhTdNjPMZVRkG2HumFKzt1uZSDIFi1koq1q9GqwUzF6MJnFGK6Q7H1DdHXElDUJyL9q9Tgt++odvyXmuyyriqaMzSGUotvU7sXkNxYGdbFGwHmeD3VZ47BIsZhPiLT4etfPkvgUghMAqmeoiFIWqpucjSd18dhgu8vPhRMVuJyC7yq5lzKXRLEaTekC4HqwUo3jV4N6bNFgUNrOpZP2CUzbrvZgCa0zNZudXJ1oUbJHFeqFZREEXzlal8z7FFsEuiTV1PcV0oVh+QfCHE3pBFqsLYEJh9L+WiyyaIJkEhGsYc2k08+HkqkR0tazUwuPgxCIIAS4ddOF8ICsU1jJiJyxFFgAC0bUFtZmAdaJF0eOQYZdMuFrrlGsxm/TvW6vCXU9Nik02IVJTi0L9IhfynS5EktgxoPbicWoWRSiuZj0ZJ9lVgl025fSMagfG/RHMhxPYPeKBP5xAj7N+QiGaBAhkuVC8en4R2/qc2N7vwnMn5wCoWU+2Aj2o8tkx0IWLS3GE4il9Gt9qYRaFVCdXXDNhMZuw949v1osaZbMJiRYXis77FFsEm2RCtIauGqPr6WvPncFvfPVFhOIKKKVajKKARRFKVJway3BYxJpmcTWCv//pCfzevx4EwCyKtU+uqwRZNC0LZh+ZCuLyoS4Mu22YCSYQV9KIKWl9FkUxPnbbFvzH/TfC65DW7HpKpDKQRKElOr7Wgv4uix6vspi564lTI2ySWLMYRSaTHSozH07giSMX8cJpP+57eD/mw0kk0xn9omc1q26jiYUoQonUqi+Gdql6qZfNwtRiDBeDcSjpjNptt46uJ0CNUxgtCkopAtEkep0WvWjzwmJMcz2V5wIihKDbJmFxja6nuJLuyPhEIaxmU8HswlaCf5JNiq2GwWxjAHQ+nMCJmTA2+uzYd8aPLz6tFuGxgBwhBNdt8uK/tDzw1frhnRax7YLZ00txUAqcngsjmcpUnDa8ViRRyPksw4kU0hmKLqtZz2Q7vxBFLJkuy/XEcNvMVbEoOjE+UQie9cSpGWsNZh+cCGD3Xz9ZMAZhDKwdnVYrpj9w3Qi29Drw6MELALKzfwHgjp19evvp1QqFXRb1xnntQCZD9crnwxfUSXH1DGYDqv/f2MKDVb67rCKG3VmhiFYsFBICEdUNuVoS3KLQ4a4nTs1Yq0VxeCoIfySJ4xeXF04ZzWA2vGZrnxO3XtKnu4eM2U237+gDy6xcbcDWLovLZm+3Mv5IEkpavZCynj71Fgo5z/XEhKLLakavU4YkCjgfiCGmpCta3bvtEpLpzJrOP25RZLGI3PXEqRE2ybQmV41fq48o1FiQWRQegxhs6XPgtkt69b+NF71epwV7tHbhq0mPBdTsqXYquDNWmTOhaITrqZBQuKxmCALBsFudlLga1xOwtloKHqPIIptN3KLg1AabLK4p95p1gJ0MLG9GxlY3Q9q0tW6bGT6HjKvWu3Xx8OQJwnt3D8Nrl/T5BZVibzOhmF7KCsUb06rrqd7BbFkUcrKeggaLAgBGvXaM+6OIJlM5LcZL0W1TP/tAZPUBbW5RZFFdT6v7Lh+ZWsK+0/4qH1HlcKFoUmxmE5Q0LWveQCHYTInJQBSZDMUvTs7rPmd20jI/9tZeJwghMAkEt13SC5/mtjDy3j1DePnPblv1IBqHrGZxZTKr93s3Exc1i8KpvS5Clotrrcm3KIIxVYiZUIxoQhFTys96ArKvg1sU1cGyhqynv//pCfzZvx+q8hFVDv8kmxSbVugWW6Wf2B9hrqcYnjo6g9/8+kv68CBmBg9pKZRb+rKT1/78nTvwyO9cv2x/TEhWi96uvE0C2heXYjAJRC9AdNukuvV5YuRnPS3lWxQ9NsSUNJQ0LauOglEN1xO3KLJYRBNSGYrUKibdzYUSuLAYW1NiQTXgQtGk2LUV4GovrFmLIqb3/zk5GwaQFR/doujLNn9zWczYoE3YqybGwr124OJSAr1OGYOa+261sZu1IJmWxyhMAtFFedSb/RwriVEw19Naaim4RZGFDYyKr8I74A8nkEhldFdyo+CfZJPCXAWrTZFlJ9ZMKI4xzZJgGU4sp/uq9d340A2jePtl/Ws93JLYDQOQ2oGZYBz9XRYMaPPD653xBGiV2XlC4bKIevM/o1BU4nrq1iyStVycuEWRhb0PlbqfKKWY1z4DY0ysEaxJKAgh5wghhwghrxJCxrRtHkLIk4SQk9pvt+HxnyCEnCKEHCeE3GHYvlvbzylCyBdIqTaXHYBdCz6uJkUxnVErdAe7raAUGDu3ACArFMyicFnM+ORdO9HrtFTpqFfGoXWQbZfq7OmlGPpdFvR3qRZFPfs8MaS8YPZSTNHdTgCwrtsCUXMXVmJRiCYBLouIxTXGKDpxXnYhLOLqhCKcSOkLgUZPyKvGJ/lWSumVlNI92t8fB7CXUroFwF7tbxBCdgC4B8BOAHcC+BIhhJ29XwZwH4At2s+dVTiuloZ9sVfTH2kxmkSGApdrM5ozFCAEOKdbFOrJJ9fxi+yQ2XjX1s4nZ8wEE6pF4WIWRQNcTwXSY10GoRBNgj40qpIYBaDWUqylg6zqeuIWBZD9nlWaIus3dHZuaYtiBe4G8JB2+yEA7zJs/y6lNEEpPQvgFIBrCCEDAFyU0n1Ujdg8bHhOx6IHs5XKV+DMZXDFcLe+7YZNXoz7o0hnKOJs4lkdXQNsJkU7xChCcQXhREqzKBrneioUzDZaFAAw4tWEosIRrd22tTUGTKQydV2INDOrdT2xhBQAmFpqbYuCAvgpIWQ/IeQ+bVsfpXQaALTfrIprEMB5w3MntW2D2u387R2NfQ0WBZsxsXOdCyaBQDIJuGNnP5LpDKYWY/oJW08fslPOzt5uddjqrr/LghGvDT0OGZdqQ2rqSX4Lj2CeRQGoKbJAZa4nAPCsod8TpWrTSW5RqFQiFBcMLqa5kMGiWGxti+JGSukuAG8HcD8h5M1FHlso7kCLbF++A0LuI4SMEULG5ubmKj/aFmItwWxmUfQ6LVjXbcH2Aaee2XR2PoKYkoYokIpHmq6Fdppyx1x4I147nBYzxv78Nrxlq6/uxyGbVYuCpU4WsihYBlvFriet39NqyA4t4hYFoE64A0q7ng5OBHDjp5/GgQk1+YRZFKNeG6Zb2aKglE5pv2cB/DuAawDMaO4kaL9ntYdPAhg2PH0IwJS2fajA9kL/70FK6R5K6R6fr/5fzHqylmA2O8G8Dgkfu3Urfv+WLfoF45w/grhS/4yUdkqPPedXhWKDt/ppxJUw7LYhmcrgri8+jzNz4YJCcfM2H27Z3ouNvsqOtb/LgplgfFULlYR2QeQWhUq5FgWzJp45ri6CWYxi52AXplrVoiCE2AkhTnYbwNsAHAbwIwD3ag+7F8APtds/AnAPIUQmhGyAGrR+WXNPhQgh12nZTh80PKdjyVoUqxCKcBKEqKvC9+wewu07+tDrlGGTTLpFUe/VniyaYDaRthCKs/NRuG1mdNnqNyO7EO+7Zj3+9j2X4/RcGJ978gRSWotxIyNeO77xoasrauEBANdt9CKVoXjlXKDi42LuMG5RqLDvcqlW46xg8kWtZYc/nECX1Yz1HhtmgnGkG9jVYC2jUPsA/LuWySoC+FdK6U8IIa8AeIQQ8hEAEwDeCwCU0iOEkEcAvAEgBeB+Sil75z4K4JsArAB+rP10NLIowCSQVa3o/JEE3DYpp5KaEIJRrx1n5yPw2KSG5Li3agfZRCoNs5Cd1nZ2PlyTosRKMQkE790zjCeOXMTPT6ir0HyhWC1Xj3ogmQQ8f2q+Irfa/vEFzAZVi5ZbFCrZ9NjirifWguXg+QBiyTTmI0l4HRLWdVuRylDMhxPoc9U+lb0QqxYKSukZAFcU2O4HcOsKz3kAwAMFto8BuHS1x9KOEEK0DrKVWxQLkWTBvkPruq2YDERhk0wNEQqH3FrDi96YCuLD33wZM8EE3n5pP778m7sBAOfmo7hhs7fBR5dl94gHTx1VPbzVEgqrZMKukW48f2q+5GNfOD2PXqeMzb1O/Nm/H8a4X+1YzC0KFYueHlv8uxyMqxaFkqYYG19Q57DbZazTMuumFmMNEwr+STYxNslUUa+nhUgSX3vuDM7ORwu2lOh1yZgNJRBLljdDudo4WqyD7OeePIG4ksG1Gzx4+tisOn86mcbFYLzh8QkjV4/qNa1VEwoAeNPmHhyZCpas0P6Tf3sd//CUOhlxLpTQux5zi0JFLjNGEYwpcMgiRIHghdN+bXa9alEA2RY8jYALRRNjlyqbCvdvY+fxfx87iqPTwYKzEXqdMhYiSYTiqYas9lpJKI5dDOKpozP48I2j+O2bNiKRyuDgxKIeyB5tAtcT47KhLr3br8tSPaG4YXMPAOAXJawKfySB2VBC7wjA4BaFCnsfEiV6PS3FFPicMq5a340n35jBfDgBr0PC1j4nNvrs+H/Pn2tYc0D+STYxNrmyKXfHZ0Lw2iV87NYt+NANG5bdz1p1TAZiDYtRNLvrKZZM48FnT+MPv/ca7JIJH7phFNds9EAgqouFpcY2Q4yCIYsmXKFV4VfTorh8sAvDHiv+4ckTK66GY8k04koG8+GE3hGAuUp4rycVySSAkNKdoIPxFFwWEfdcvR6nZsMIRBV47TJMAsF/v3kzjk4H8fSx2aL7qBVcKJoYm7myudknZ8K4ZMCFP7x9K67Z4Fl2f6/Wj+hiMN6YGIVFRKjJheKnb1zE3zx+DJFECp+8aye6bRJcFjMuH+rGC6f9ONuEFgUA7BlVP+9qZmKJJgF/8yuX4cx8BF/Ye7LgYxY0C8IfTuouqj+4dQv+9B3bcaWhM0AnQwgpaxwqK5j85SvWwad9V1lrmLuvXIchtxX//PMzNT/eQnChaGIqsSgyGYpTs+Gc2RL5GKfTNUIouqxmLK2hf1A9uKhVXT/2B2/Ce/dky35u3OzFq+cXsffoLHocst7Ku1m476aN+PL7d1XVogCAm7b4cPeV6/C1584WHDoV0MRhKaboFevDHhvue/OmuhZ0NjtWyVQyPTYYV4VCEgXce/0IAMCrtYYxmwTcdcU6HJgI5FTj1wv+STYxlczNngzEEFPS2GaYLZGPsUustQH+Y59DxkI0CWUVA1zqxVwoAavZtL9/HqAAABMNSURBVEwIbtriQzpD8er5Rdx95boGHd3KuO0S3n7ZQE32fcVQN5LpjJ6VY8QYk2DB1nrPDm8FLKJQRnqsoseYPnD9KD54/Qiu35jNrtu5rgupDMXJmfoHtZtrWcTJYWOPAz85fBGnZkPY3LuyAADAiZkQAGBLEaHocUggBKC0MRaFzymDUjU7q1FpfqWYDSXgc8rI73R/7QYPfvDR67HJ59AH+3QKbnt2PkX+azd2mD1+UZ0d7rXXv0Fis1NqHCqlFMFYSrcIu6xm/NXduRUDO7VpikemlureW4xbFE3Mh28chdVswuefLOwfNnJilgnFyq4n0SToabONSI9lHVbnQokSj2wcc6GEHssxQgjB7hFPx4kEkJ14ZxSFd3/peXz9F2d11xMAHNdWuu4GV6w3I7LZVNSiSKQySKYzcFlXXruv99jgkEW8MRWsxSEWhQtFE+N1yPhvb9qAxw5Nlzw5Ts6EMdBlKZke6dPcT3KDLAoAmAs3r1DMhuL6cXJUPPpoVG1qYjCOAxOLePmsP6fG4uRMCG6bue6zw1sBi1koGltg7TuKfX8FgeCSASeOcKHg5PPhG9U012dPrtwtdyGSxKELS0XdTgy2Wm6ERcH+93wLWhSdjFsTCiYKr08uAVCb2C1Gk3oNRzSZLtgRgIOSWU9BJhQlkhF2DLhwdDpYMLGglnChaHI8dgleu4RxLS0znyffmMHVDzyFU7NhXFsgJTYfdhFsRDGU7npqUosirqQRjKe4RZFHtxajWNRcT69PLgJQEygWogrWdVn0hYe3AQOcWgGLWdAr1gvBEgVKZa3tWOdCJJnG+EK0qsdXCh7MbgFGvDacmy98Yrx6Xu3u+V+//yY92FUMliLbiGC2VVKziZo1RsGOqx4zxFsJp9ZWgmU4MYtiMargQiAKt11CmlKcX4gVbB3DYcHslWMUrCGgy1L8krxznRrEPjK1VNeiT25RtACjXvuKFoU/nITXLuHSwa5lmTqFYBfBRrieADXzaj68+hGbtWRWEwqfi6+KjRBC9NGolFK8PrkIp3ZBe2M6CLdN0q1FnhpbGGuJrKelMl1PW/ocuH6jt+7fXy4ULcCI146ppXjBE20+nKzI3G+k6wlQA9pzofoMYTl8YQnJEv11jLDj8nH3yTLcNjMCEQWTgRgCUQW37+gDoLbOzhEKnhpbkFJZT+W6nmTRhO/cdx1uvaSvqsdXCi4ULcBojw0AcL6AX9IfSehl/uXtSzVXG+WH9znlulgUF5fiuOuLv8CjByZLP1hDdz1xi2IZbs2iYG6nd1yaLe7z2M36OcgtisJYzAISZQSznSVcT42CC0ULMOplY0wLCIXmeiqXSwZcePqP34Jd692lH1wDehxyXWIUZ+bDyFDgVAWtmWdDCQiEr4oL4babsRhVcHQ6CJNA8KYtPZC0NNhublGUxGIu3sIjqHV0btbW7FwoWgAmFIXiFP5wAp4Kv5wbfY6y4hm1wOeQsRRTat6vhllf5wPq77iSLtmieS6UgNch50wG5Ki4bRIWokmcmg1jxGODxWzCum413uWxZ4WCp8cWxiKaoKTpiuNMl6JKVVvEVxsuFC1Al82MbptZn4XAiCXTiCTTLWXuM5eXv8buJzZlbWIhhmgyhes+tRff35/rhqKU4t5vvIx/fWkCgGpR8BqKwnTbJCxGkzg1F8amXrX6f9CtDtRx2yRs9NlhEgjWe22NPMymxabNzQ4V6JcFqDGKajd0rCZcKFqEEa9dv/gx/BHVhVNJjKLR1KuNxwSzKBaiOHYxhMWogsMXlnIec84fxc9PzOk9/nlV9sq4bWYoaYozc2FsZkLRbdXve9PmHuz7xC36Nk4uOwfV1PX944GC97POsc0KF4oWYdRrw9n5XIuCrcpbyS+st/GosVAw11M4kcK+034AaoGYkee0avcz82FQSnEhEMNAF6+hKIRbcyllKLDZpwrFkFu1Hjx2CYQQXn9ShF3r3ZBFoeC0wGgyhVOz4aZ223GhaBE29jhwYTGWMyWLWRSt6HqaqXGK7PhCFH1a9tITRy4CKCQU6pd2wh/F9FIcgaiCrWW0QelE3IZmiMyiuH6TF1v7HLoLirMyFrMJ12zw4IVT/mX3febHxzAbSuC3b9rYgCMrDy4ULcLWPgdoXhYPSzPtaaG8/z6XBbIo6CNFa8FSTMFiVMGbNvsAZCuJJwNRPaCtpDPYd9oPl0VEKkOx9+gMAGBbPxeKQhg7wm70qckVV4968NM/fAtsUnOmdDYbN2zqwfGZEGa1RdKzJ+Zw38NjeGjfOD58w4aCUymbBS4ULQJr+MfmTgAG11MLWRQmgWBDjx1n5monFMztdOPm7NAXSRQQSab1fkWvnV9EOJHCr2lT7H58WLU6tveXboPSiTDXU7/LAmcTZ+c0M2/a3AMA2HfaDyWdwf3fPoCD5xfxwetH8Cd3bGvw0RWHC0WLMOq1QTIJ+twJQE2NtZpNLbei29TrwOm52k3pYoHsbf1O3dq6cZMqGsz99NzJeQgEeP916sjJF8/44XPKTe0nbiTM9cTcTpzK2bHOBbfNjKePzeLV84sIJVL4q7t24q/uvhRWqTnrJxhcKFoE0SRgo8+eMwbRH0m2lDXB2ORzYGIhWvVaikQqjb974jgeOzQNQB30st6j+s9v39EPQHU/AWog+7KhbmzosaPbZkaGAtu522lFuqxmmATChWINmASCOy8dwJNvzOCJwxchENUd1QpwoWghtvQ5cWImhIMTAfztE8f0ArFWY5PPjgzFsnTftbL36Cy++LNTeOz1aXjsEpwWM9Z7bCAEuGV7LwDVoliKKXhtcglv3qJ+SVkXzmLzxjsdk0Dwld/cjd95S/MGXFuBu69ch2gyjYf2ncOVw93oapFpgK3ls+hwtvY68J+vTeEPv/cqzvmjMAkEN2/1NfqwKmaTll55ejZc1Syjxw5No8ch4c9/aYduyv/a1cPY6HOgv8sCp0XEZCCKfaf9SGeo7jPe2OPAwYlFHsguAWsEyFk914x6MNBlwfRSHDdtaZ3vLrcoWggW0D7nj8Jrl5DO0JZ0PbEV/FrjFEsxRe+6GUum8fTRWdyxsx/vumoQd+xUXU03bOrBH9y6BYCa9z8ZiOEXp+Zgl0y4Sut3xbJ4eCCbU2sEgeCuK9YBAN7cQos8blG0EFv71JX49n4nPvuey3H3Pz2PflfrFTnZZRHruiw4rWU+Uar2wKlk1nJcSeNXvvQ8lHQGj//BTfjFyXnElDR+6bKBFZ8z7LbiyFQQR6eDuG6jVx/h+c7LBzAbjGP7ALcoOLXnvjdvxLpuK64a7m70oZQNtyhaiBGvHe+6ch3+6u5LcflQN77/uzfoM7VbDWPm09d/cRY3fuZpxJU0Xjg9jw98/SXMlxiX+uVnTuPMXASTgRj+6JHX8LknT8Brl4rmog+5bbiwqM5T+M3rR/TtI147/s/dl8JcgVBxOKvF65Bx7w2jEFqo+SS3KFoIk0DwD/dcpf+9e6QxrcKrwdY+J7714jjCiRR+cOACZoIJPHtiDt96aQLPnZzHf//WAXzrt67VV/1GJvxRfPmZ07jrinVY77Hhiz87BZ9Txt/92hVFrZIPXD+CXpeM9+weaqkiRQ6n0XCh4DSEd1zWj6//4iwefPYMjk4HAQAP7xvHvjN+XDncjZfPLeBzT57Ax9++Hd95eQJ9Lhm3bFeDqd9+aRxpSvGn77gEXoeE7QNO3LTFV7L75oYeO373LZtq/to4nHaDCwWnIexa78aGHjv+6WenAAA3bPLqDdM++57L8eCzZ/CNX5zFxh47PvHoIRAC/J+7duJ916zHDw5M4pbtvejXGvi98/J1DXsdHE4n0DROWULInYSQ44SQU4SQjzf6eDi1hRCCd181iHSGYnu/Ex+9WV3pb+93YmufE//zbdsgCMD/+sHrGPHacOv2PvzvHx7Bbz00hvlwEvdcPdzgV8DhdA5NIRSEEBOAfwLwdgA7ALyPELKjsUfFqTW/smtQq1btx/UbvbhiqAv/TQvO93dZ8Dtv3gRCgE+9+zJ85Td34df2DOHnJ+bQ65TxlhZKLeRwWh1SajxkXQ6CkOsBfJJSeof29ycAgFL6qZWes2fPHjo2NlanI+TUipMzIQxrozXzoZTiYjCOgS6r/vdDL5zDoNvGi784nFVCCNlPKd1TyXOaJUYxCOC84e9JANc26Fg4dWRLkcpsQoguEuzvD7VoOjCH08o0hesJQKGE4mWmDiHkPkLIGCFkbG5urg6HxeFwOJxmEYpJAMbo5BCAqfwHUUofpJTuoZTu8fm4j5rD4XDqQbMIxSsAthBCNhBCJAD3APhRg4+Jw+FwOGiSGAWlNEUI+T0ATwAwAfgGpfRIgw+Lw+FwOGgSoQAASunjAB5v9HFwOBwOJ5dmcT1xOBwOp0nhQsHhcDiconCh4HA4HE5RmqIyezUQQkIAjpf58B4A81X8910AlppwX7XYXzO/d62wv056//h711z7y3//2N8jlNLK6gsopS35A2CsFo8tc38PNuO+arS/pn3vWmR/HfP+8feu6fY3VuzvSn6462l1/GeT7qsW+6s2zf56+fvXHPuqBc383tVif1WjlV1PY7TMxlaVPJaTC3/v1gZ//1YPf+/WRv77t5b3s5Utigdr9FhOLvy9Wxv8/Vs9/L1bG/nv36rfz5a1KDgcDodTH1rZouBwOBxOHWhJoSCEDBNCfkYIOUoIOUII+Zi23UMIeZIQclL77da2e7XHhwkhX8zb168TQl7X9vPZRryeerKK9+52Qsh+Qsgh7fcthn3t1rafIoR8gRBSqF18W1Hl9+8BQsh5Qki4Ua+nnlTrvSOE2AghjxFCjmn7+XQjX1e9qPK59xNCyGvafr6iTRldmWqmY9XrB8AAgF3abSeAE1BHqH4WwMe17R8H8Bntth3AmwD8LoAvGvbjBTABwKf9/RCAWxv9+prsvbsKwDrt9qUALhj29TKA66HOE/kxgLc3+vW12Pt3nba/cKNfVyu9dwBsAN6q3ZYAPMfPvYrPPZf2mwD4AYB7iv7vRr/4Kr2BPwRwO9QCvAHDm3o873EfyhOKqwE8Zfj7AwC+1OjX04zvneGk8gOQtcccM9z3PgD/3OjX0yrvX972jhCKWrx32n3/COC3G/16WvH9A2CGmpb768X+V0u6nowQQkahKudLAPoopdMAoP3uLfH0UwC2E0JGCSEigHchd4BSW7OK9+5XARyklCagjq+dNNz3/9u7vxApqzCO498fjQSxogUpQkR30h+iKKjICLoJheimm4h2M4L+QXgXRVAXeVGUiHaxiRlpIVIUWEESCwkVQTdS2gayIiRIEW66KV1ETxfnTA3LzKnRM847O78PvMxy5p3znvdh2GfOmXee90RuGxsXGL+xVit2klYC9wEzgxxv09SIn6QDwC/AAvBB6XgjnSgkTZCmTZsi4ky/r4+IeeBJYB9p+noc+LPmGJuq39hJuh54BXi83dRlt7G5hK5C/MZWrdjlD3d7gW0RcWwQY22iWvGLiHtJM5BLgXu6vPQfI5soJC0jBeu9iPgwN/8saU1+fg0pWxZFxMcRcVtE3EGawh0d1Jibot/YSboK+AiYjIi53HyCdMvatq63r12KKsVvLFWO3Q7gaERsHfzIm6H2ey8i/iDdTfT+0nFHMlHkq2veAmYjYkvHU/uBqfz3FGkN77/6WpUfLweeAnbWHW2z9Bu7PLX/FHguIr5q75ynuAuSbs99TvI/4j3qasVvHNWMnaSXSUX0Ng163E1RK36SJjoSSwvYAPxYPPiwv5A5zy9x1pGWOb4DDuVtA+kqphnSrGAGuKLjNceBU8DvpE/D1+X2vcAPeSt+878Utn5jB7wAnO3Y9xCwKj93K3AYmAPeIP+AcylvleP3an4v/pUfXxr2+Y1C7Eiz1wBmO9ofG/b5jVD8VgPf5n6OANuBVunY/mW2mZkVjeTSk5mZXTxOFGZmVuREYWZmRU4UZmZW5ERhZmZFThRmAyDpCUmTfex/jaTDgxyT2flqDXsAZkuNpFZETA97HGa1OFGYdZGLrn1GKrp2M6mk8yRwLbAFmAB+BR6JiJOSvgC+Bu4E9ktaTqoK+5qkm4BpUnnsOeDRiJiXdAuwCzgHfHnxzs6sP156MuttLbAjIm4EzgBPk37F+kBEtP/Jb+7Yf2VE3B0Rry/qZzfwbO7ne+DF3P428EykOmNmjeUZhVlvP8W/NXLeBZ4n3QDm81R2h0uAkx3771vcgaQVpARyMDe9A7zfpX0PsL7+KZhdOCcKs94W17dZAI4UZgBn++hbXfo3ayQvPZn1drWkdlJ4EPgGuLLdJmlZrvXfU0ScBuYl3ZWbHgYORsRvwGlJ63L7Q/WHb1aHZxRmvc0CU5LeJFXm3A4cALblpaMWsJVUgbNkCpiWdBlwDNiY2zcCuySdy/2aNZKrx5p1ka96+iQibhjyUMyGzktPZmZW5BmFmZkVeUZhZmZFThRmZlbkRGFmZkVOFGZmVuREYWZmRU4UZmZW9DcLOMgjUY2PgwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sorted_data['inc'][-200:].plot()" ] }, { "cell_type": "code", "execution_count": 47, "metadata": {}, "outputs": [], "source": [ "first_sept_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": 48, "metadata": {}, "outputs": [], "source": [ "year = []\n", "yearly_incidence = []\n", "for week1, week2 in zip(first_sept_week[:-1],\n", " first_sept_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": 49, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 49, "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": 50, "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": 50, "metadata": {}, "output_type": "execute_result" } ], "source": [ "yearly_incidence.sort_values()" ] }, { "cell_type": "code", "execution_count": 51, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 51, "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)" ] } ], "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 }