{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Dr. Semmelweis and the Discovery of Handwashing\n", "> A Summary of project in datacamp\n", "\n", "- toc: true \n", "- badges: true\n", "- comments: true\n", "- author: Chanseok Kang\n", "- categories: [Python, Datacamp]\n", "- image: images/Ignaz_Semmelweis_1860.jpg" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "13f090f9f0" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 1. Meet Dr. Ignaz Semmelweis\n", "

\n", "\n", "

This is Dr. Ignaz Semmelweis, a Hungarian physician born in 1818 and active at the Vienna General Hospital. If Dr. Semmelweis looks troubled it's probably because he's thinking about childbed fever: A deadly disease affecting women that just have given birth. He is thinking about it because in the early 1840s at the Vienna General Hospital as many as 10% of the women giving birth die from it. He is thinking about it because he knows the cause of childbed fever: It's the contaminated hands of the doctors delivering the babies. And they won't listen to him and wash their hands!

\n", "

In this notebook, we're going to reanalyze the data that made Semmelweis discover the importance of handwashing. Let's start by looking at the data that made Semmelweis realize that something was wrong with the procedures at Vienna General Hospital.

" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "dc": { "key": "13f090f9f0" }, "tags": [ "sample_code" ] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " year births deaths clinic\n", "0 1841 3036 237 clinic 1\n", "1 1842 3287 518 clinic 1\n", "2 1843 3060 274 clinic 1\n", "3 1844 3157 260 clinic 1\n", "4 1845 3492 241 clinic 1\n", "5 1846 4010 459 clinic 1\n", "6 1841 2442 86 clinic 2\n", "7 1842 2659 202 clinic 2\n", "8 1843 2739 164 clinic 2\n", "9 1844 2956 68 clinic 2\n", "10 1845 3241 66 clinic 2\n", "11 1846 3754 105 clinic 2\n" ] } ], "source": [ "# importing modules\n", "import pandas as pd\n", "\n", "# Read datasets/yearly_deaths_by_clinic.csv into yearly\n", "yearly = pd.read_csv('./dataset/yearly_deaths_by_clinic.csv')\n", "\n", "# Print out yearly\n", "print(yearly)" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "45ea098e15" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 2. The alarming number of deaths\n", "

The table above shows the number of women giving birth at the two clinics at the Vienna General Hospital for the years 1841 to 1846. You'll notice that giving birth was very dangerous; an alarming number of women died as the result of childbirth, most of them from childbed fever.

\n", "

We see this more clearly if we look at the proportion of deaths out of the number of women giving birth. Let's zoom in on the proportion of deaths at Clinic 1.

" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "dc": { "key": "45ea098e15" }, "tags": [ "sample_code" ] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " year births deaths clinic proportion_deaths\n", "0 1841 3036 237 clinic 1 0.078063\n", "1 1842 3287 518 clinic 1 0.157591\n", "2 1843 3060 274 clinic 1 0.089542\n", "3 1844 3157 260 clinic 1 0.082357\n", "4 1845 3492 241 clinic 1 0.069015\n", "5 1846 4010 459 clinic 1 0.114464\n" ] } ], "source": [ "# Calculate proportion of deaths per no. births\n", "yearly['proportion_deaths'] = yearly['deaths'] / yearly['births']\n", "\n", "# Extract clinic 1 data into yearly1 and clinic 2 data into yearly2\n", "yearly1 = yearly[yearly['clinic'] == 'clinic 1']\n", "yearly2 = yearly[yearly['clinic'] == 'clinic 2']\n", "\n", "# Print out yearly1\n", "print(yearly1)" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "2bc9206960" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 3. Death at the clinics\n", "

If we now plot the proportion of deaths at both clinic 1 and clinic 2 we'll see a curious pattern...

" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "dc": { "key": "2bc9206960" }, "tags": [ "sample_code" ] }, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Proportion deaths')" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEGCAYAAAB/+QKOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3dd3hUZfbA8e9JpwZIIJTQIXRpAUFXLEixgdjAgoKFVddV17a6utbVXbvu6voTBQW7oii6IgjYRSUUhQCB0AOSQIAQCCHt/P64g8QQyCSZyZ2ZnM/zzMPkzr1zzzVxzty3nFdUFWOMMaasMLcDMMYYE5gsQRhjjCmXJQhjjDHlsgRhjDGmXJYgjDHGlCvC7QB8JT4+Xtu1a+d2GMYYE1QWL168U1WblvdayCSIdu3akZKS4nYYxhgTVERk09FesyYmY4wx5bIEYYwxplx+TRAiMlJE0kQkXUTuLOf1ISKyRESKROSCMq+1EZG5IrJKRFaKSDt/xmqMMeb3/NYHISLhwPPAMCADWCQis1R1ZandNgMTgNvKeYvpwMOq+rmI1AdK/BWrMSZ0FBYWkpGRQX5+vtuhBJSYmBgSExOJjIz0+hh/dlIPBNJVdT2AiLwNjAZ+SxCqutHz2u8+/EWkOxChqp979tvnxziNMSEkIyODBg0a0K5dO0TE7XACgqqSnZ1NRkYG7du39/o4fzYxtQK2lPo5w7PNG0nAHhH5QESWisjjnjuS3xGRSSKSIiIpO3bs8EHIxphgl5+fT1xcnCWHUkSEuLi4St9V+TNBlPfb8bZ0bARwEk7T0wCgA05T1O/fTHWyqiaranLTpuUO4zXG1EKWHI5Ulf8m/kwQGUDrUj8nAtsqcexSVV2vqkXAh0A/H8dXaxUUlfDBkgz2HyxyOxRjTADzZ4JYBHQWkfYiEgWMA2ZV4tjGInLotuA0SvVdmOqZvnAjt7z7M/+cvcrtUIypNe6//36eeOIJAO69917mzZt3zP1nzZrFv/71L6/f/+6776Z169bUr1+/WnGW5rcE4fnmfwMwB1gFvKuqqSLyoIiMAhCRASKSAVwIvCgiqZ5ji3Gal+aLyHKc5qqX/BVrbbI3v5Dnv0gnOiKMN37czNLNu90OyZha58EHH+T0008/5j6jRo3izjuPmB1wVOeccw4//fRTdUP7Hb/Og1DVT1U1SVU7qurDnm33quosz/NFqpqoqvVUNU5Ve5Q69nNVPU5Ve6nqBFUt8GestcXLX69nd14h064cSLMG0fxt5gqKim0EsTG+NH36dI477jh69+7N+PHjj3h9woQJzJgxA3DKBN13333069ePXr16sXr1agBeffVVbrjhBgAyMzMZM2YMvXv3pnfv3nz//fdHvOegQYNo0aKFT68jZGoxmYrtyD3Iy99u4OzjWjCoQxz3n9OD695YwivfbeSaIR3cDs8Yn3vg41RWbtvr0/fs3rIh953T46ivp6am8vDDD/Pdd98RHx/Prl27KnzP+Ph4lixZwn//+1+eeOIJXn755d+9fuONN3LyySczc+ZMiouL2bevZkb+W6mNWuS5BWs5WFTCrcO7ADCyZ3NO69qMpz5fw9Y9B1yOzpjQsGDBAi644ALi4+MBaNKkSYXHnHfeeQD079+fjRs3lvue1113HQDh4eHExsb6LuBjsDuIWmJzdh5v/rSZsQNa0z6+HuAMe3tgVA+GPf0V932UystXJLscpTG+daxv+v6iqpUeUhodHQ04H/5FRYEzutDuIGqJp+etIUyEm4Z2/t321k3qcvPpScxblcmc1O0uRWdM6Bg6dCjvvvsu2dnZAF41MXnzni+88AIAxcXF7N3r22azo7EEUQus+nUvHy7bysQT25PQMOaI16/6Q3u6Nm/A/bNSbW6EMdXUo0cP7r77bk4++WR69+7NLbfcUu33fPbZZ/niiy/o1asX/fv3JzU19Yh97rjjDhITE8nLyyMxMZH777+/2ucVVW8nNwe25ORktQWDynfVq4tYtHEX39xxGrF1yy/UtXjTLs5/YSFX/6E995zdvYYjNMZ3Vq1aRbdu3dwOIyCV999GRBararnty3YHEeIWbdzF/NVZXHtKx6MmB4D+bZtw8cA2vPL9RlK35dRghMaYQGUJIoSpKo/OXk2zBtFMPKHiCo53juxK47qR/G3mCopLQuPO0hhTdZYgQtgXaVmkbNrNjUM7UyfqiGK4R4itG8k9Z3Xn5y17ePPHoy5Ta4ypJSxBhKiSEuWxz9JoG1eXsQNaV3yAx+g+LflDp3ge+yyNrL224IoxtZkliBA16+dtrN6ey63DuxAZ7v2vWUR46NyeHCwu4cFPrD6iMbWZJYgQVFBUwpOfp9G9RUPO7lX52izt4+vxp1M68ckvv/LVGluIyZjayhJECHp70Wa27DrAHSO7EBZWtYVTrj2lAx2a1uPvH64gv7DYxxEaU/v4s9x3Xl4eZ511Fl27dqVHjx6VqgJ7LJYgQsz+g0X8e346x7dvwslJVV9lLzoinH+c25PNu/L4z4K1PozQGOOPct+33XYbq1evZunSpXz33XfMnj27umFaggg1r3y3gZ37DnLHyK7VXnbxhI7xnNevFZO/Xs/azFwfRWhM6Kvpct9169bl1FNPBSAqKop+/fqRkZFR7euwYn0hZPf+Al78aj3DuifQv21jn7zn3Wd2Y8HqLO6euYK3Jw2qcpOVMa6YfSdsX+7b92zeC844etOP2+W+9+zZw8cff8xNN93k/TUdhd1BhJAXvlrHvoIibh/RxWfvGVc/mrvO6MpPG3cxY3H1v5EYE+rcLPddVFTExRdfzI033kiHDtVf48XuIELErzkHePX7jZzXN5GkhAY+fe8L+7dmxuIMHpm9iqHdmhFXP9qn72+M3xzjm76/uFnue9KkSXTu3Jmbb765yu9Rmt1BhIhn560FhZtP71zxzpUUFiY8PKYX+/KLeOTT1T5/f2NCiVvlvu+55x5ycnJ45plnqn2+Q/yaIERkpIikiUi6iBzRHS8iQ0RkiYgUicgF5bzeUES2ishz/owz2KVn7ePdlC1cOqgNrZvU9cs5khIaMGlIB95fksHCddl+OYcxocCNct8ZGRk8/PDDrFy5kn79+tGnT58j+jGqwm/lvkUkHFgDDAMygEXAxaq6stQ+7YCGwG3ALFWdUeY9ngWaArtU9YZjna82l/u+/o3FfJW2g6/uOJV4Pzb/HCgoZvgzXxEZHsbsm04iOqLi+k7G1DQr9310gVTueyCQrqrrVbUAeBsYXXoHVd2oqr8AJWUPFpH+QAIw148xBr2ft+zh0+XbuWZIB78mB4A6UeE8NLon63fs58Wv1vv1XMYY9/kzQbQCtpT6OcOzrUIiEgY8Cdzuh7hCyuNz0mhSL4qrT6r+iAVvnNKlGWcd14Lnvkhnw879NXJOY4w7/JkgyuvG97Y963rgU1XdcqydRGSSiKSISMqOHbWvZtC3a3fybfpObji1E/Wja25A2n1ndyc6PIy/f7iCUFmR0IQW+7s8UlX+m/gzQWQApetMJwLbvDx2MHCDiGwEngAuF5Ejxqup6mRVTVbV5KZNq15WIhipKo/NWU2rRnW4dFCbGj13s4Yx3D6yC9+m72TWz97+So2pGTExMWRnZ1uSKEVVyc7OJibmyDXpj8WfXzsXAZ1FpD2wFRgHXOLNgap66aHnIjIBSFZV31SfChGfrdjOLxk5PHFhb1c6iy89vi3vL87goU9WckpSs2MuZ2pMTUpMTCQjI4Pa2KpwLDExMSQmJlbqGL8lCFUtEpEbgDlAODBVVVNF5EEgRVVnicgAYCbQGDhHRB5Q1R7+iilUFBWX8PjcNDo3q8+Yvl516/hcuGduxKjnvuXROat5ZEwvV+IwpqzIyEjat694iV1TMb82XKvqp8CnZbbdW+r5Ipymp2O9x6vAq34IL2i9vySD9Tv2M3l8f8JdrI3Us1UsE05oz9TvNnB+v0Sf1X8yxgQGm0kdZPILi3lm3lr6tmnEsO4JbofDLcOTaBEbw90zl1NYfMRoZWNMELMEEWReW7iJX3Py+asPynn7Qv3oCO47pwert+cy9dsNbodjjPEhSxBBZG9+Ic9/mc7JSU0Z1CHO7XB+M6JHAqd3a8Yz89aSsTvP7XCMMT5iCSKITP5qPXvyCn1aztsXRIT7RzljC+77KNWGFxoTIixBBIms3HymfLuBc3q3pGer8mvBuymxcV3+Mqwz81dnMSd1u9vhGGN8wBJEkHhuQTqFxSXcOizJ7VCOauKJ7enavAH3z1rJvoNVr2lvjAkMliCCwObsPN78cTNjB7SmXXw9t8M5qsjwMB45rxeZufk8OTfN7XCMMdVkCSIIPPV5GhHhwo1Dfb8YkK/1a9OYSwa2Ydr3G1mxNcftcIwx1WAJIsCt3LaXj37exsQT25PQsHJ1VNxyx8iuNKkXzd9mLqe4xDqsjQlWliAC3BNz02gQHcG1Qzq6HYrXYutE8vezu/FLRg6vLdzodjjGmCqyBBHAftqwiwWrs7julE5BVwxvVO+WnNQ5nifmrmF7Tr7b4RhjqsASRIBSVR77bDXNGkQz4YR2bodTaSLCQ6N7UlBcwoOfpFZ8gDEm4FiCCFALVmeRsmk3N53emTpRwbn2c7v4evz51E58unw7X6zOcjscY0wlWYIIQMUlymOfpdEuri4XJbeu+IAANunkDnRsWo+/f7SCAwXFbodjjKkESxABaNbPW0nLzOXW4V2IDA/uX1F0RDgPj+lFxu4D/HvBWrfDMcZUQnB/+oSggqISnpy7hh4tG3JWrxZuh+MTgzrEcUH/RF76ej1p23PdDscY4yVLEAHmrZ82k7H7AHeM7EqYi4sB+drfzuxGg5gI7p65nBKbG2FMULAEEUD2HyziPwvWMqhDE4Z0jnc7HJ9qUi+Ku87sRsqm3bybssXtcIwxXrAEEUCmfruBnfsKuCNAFgPytQv7JzKwfRP+OXs1O/cddDscY0wFLEEEiF37C5j89XqGd0+gX5vQXNtZRHhkTE/yCop45H+r3A7HGFMBvyYIERkpImkiki4id5bz+hARWSIiRSJyQantfURkoYikisgvIjLWn3EGghe+TGd/QRG3BdhiQL7WqVkD/jikIx8s3cr36TvdDscYcwx+SxAiEg48D5wBdAcuFpHuZXbbDEwA3iyzPQ+4XFV7ACOBZ0Skkb9iddu2PQeYtnAT5/VLJCmhgdvh+N0Np3WibVxd7vlwBQeLbG6EMYHKn3cQA4F0VV2vqgXA28Do0juo6kZV/QUoKbN9jaqu9TzfBmQBTf0Yq6uenbcWFG4+PfDLeftCTGQ4D43uyfqd+3nhy3Vuh2OMOQp/JohWQOnhKhmebZUiIgOBKOCITxIRmSQiKSKSsmPHjioH6qb0rH28t3gLlw1qS2Ljum6HU2OGJDXlnN4t+e8X61i/Y5/b4RhjyuHPBFHeMJxKDYAXkRbAa8BEVS0p+7qqTlbVZFVNbto0OG8wnpybRp3IcP50avCU8/aVv5/djejIMO75cAWqNjfCmEBTYYIQkXoiEuZ5niQio0TEm9rTGUDpQkKJwDZvAxORhsD/gHtU9QdvjwsmP2/Zw+wV27lmSAfi6ke7HU6Na9YghjtGduX7ddl8uGyr2+EYY8rw5g7iayBGRFoB84GJwKteHLcI6Cwi7UUkChgHzPImKM/+M4HpqvqeN8cEo8fmrKZJvSiuPqmD26G45tKBbejTuhH/+GQVe/IK3A7HGFOKNwlCVDUPOA/4j6qOwRmVdEyqWgTcAMwBVgHvqmqqiDwoIqMARGSAiGQAFwIvisihhQMuAoYAE0RkmefRp9JXF8C+XbuT79KzueHUTtSPjnA7HNeEhQmPjOnFngOFPPrZarfDMcaU4s0nk4jIYOBS4KpKHIeqfgp8WmbbvaWeL8Jpeip73OvA696cIxipKo9+tppWjepw6aA2bofjuu4tG3Llie146ZsNnN8vkeR2TdwOyRiDd3cQNwF3ATM9dwAdgC/8G1Zom71iO8u35vCXYUlERwTnYkC+dvPpSbSMjeHumSsoLD5iPIIxxgUVJghV/VpVR6nqo56f16vqjf4PLTQVFZfwxJw0Ojerz5i+lR71G7LqRUfwwOiepGXm8vI3G9wOx5igkbJxF9+v809VAm9GMSWJyGQRmSsiCw49/BJNLTBjcQbrd+7n9hFdCA+hct6+MKx7AsO7J/Ds/DVs2ZXndjjGBLyV2/Yy8dVFPPjxSor9UEbfmyam94ClwD3A7aUeppLyC4t5Zt5a+rZpxLDuCW6HE5DuH9WDcBHu/cjmRhhzLBt27ufyqT9RPzqCKRMG+OULpzcJokhVX1DVn1R18aGHzyOpBaYv3Mj2vfn8NUTLeftCy0Z1+MuwJL5I28HsFdvdDseYgLQ9J5/LXv6RElVeu+p4WjWq45fzHDVBiEgTEWkCfCwi14tIi0PbPNtNJeQcKOT5L9ZxclJTBnWIczucgDbhhHZ0b9GQBz5OJTe/0O1wjAkou/YXcNmUH8k5UMi0iQPp1Ky+3851rDuIxUAKcAVOk9L3nm2HtptKmPz1OnIOFHJ7iJfz9oWI8DAeOa8XWbkHeXLuGrfDMSZg7DtYxMRXfmLLrjxeviKZXomxfj3fUeczqGp7ABGJUdX80q+JSIxfowoxWXvzmfrtRs7p3ZKerfz7Cw0VfVo3YvygtkxbuJHz+rXiuMSQrfZujFfyC4uZND2FFdv2Mnl8/xppifCmD+J7L7eZo/jPgnQKi0u4dViS26EEldtGdKFp/Wj+NnM5RTY3wtRiRcUl3PjWUr5fl82TF/ZmaLeaGeRyrD6I5iLSH6gjIn1FpJ/ncQpQe+pSV9Om7P289dNmxg5oTbv4em6HE1QaxkRy7zndWbF1L9MXbnI7HGNcUVKi3PnBcuauzOSBUT04twbnTx2rZMYInNXeEoGnSm3PBf7mx5hCylOfryEiXLhxaO1YDMjXzurVgveSMnhybhpn9GpOi1j/jNYwJhCpKv/43ypmLM7glmFJXHFCuxo9/1HvIFR1mqqeCkxQ1VNLPUap6gc1GGPQSt2Ww0fLtnHlie1JaGjdNlUhIjw0uidFJcoDs1a6HY4xNeq5BelM/W4DV57Ynj+f1qnGz19h0T1VfV9EzgJ6ADGltj/oz8BCwRNz0oitE8kfT659iwH5Upu4utw4tDOPz0lj/qrMGmt/NcZN0xdu5MnP13B+v0TuOaubK3OnvCm18X/AWODPOKvEXQi09XNcQe/H9dl8kbaD607pSGwdb9ZXMsdyzUkd6NysPvd+lEpeQZHb4RjjVx8u3cq9H6UyrHsCj57fizCXyvJ4M4rpBFW9HNitqg8Ag/n9SnGmDFXlsTlpJDSM5orB7dwOJyRERThzI7buOcCz89a6HY4xfjN/VSa3vvczgzvE8Z+L+xIR7s+VoY/NmzMf8PybJyItgUKgvf9CCn7zV2WxeNNubhqaRJ0oK+ftKwPaNWFscmte/nYDq37d63Y4xvjcD+uzuf6NJfRs2ZCXrkgmJtLdzw9vEsQnItIIeBxYAmwE3vZnUMGsuER5fE4a7ePrcWHyEWshmWq684yuxNaJ5G8zl1Pih+qVxrhlxdYcrp6WQusmdXll4sCAWGnSm/UgHlLVPar6Pk7fQ1dV/bv/QwtOHy3bSlpmLrcOTyLSxVvDUNW4XhR3n9mNpZv38NaizW6HY4xPpGft4/KpPxFbJ5LXrhpIk3pRbocEeNdJXVdE/i4iL6nqQaCZiJxdA7EFnYNFxTz1+Rp6tmrImT1buB1OyDqvXysGd4jj0dmr2ZF70O1wjKmWrXsOcPmUHwkTeP3q4wNqro83X3FfAQ7idE4DZAD/8FtEQeytHzeTsfsAd4zo6tqog9pARPjHmJ7kF5bw8P9sboQJXjv3HWT8yz+Se7CIaVcOpH2AVVvwJkF0VNXHcDqnUdUDOMNdKyQiI0UkTUTSReTOcl4fIiJLRKRIRC4o89oVIrLW87jCm/O5ad/BIv6zIJ3BHeI4qXO82+GEvI5N63PtKR35cNk2vl3rn+UWjfGnvfmFXDH1J7blHGDqhAH0aBl4hTy9SRAFIlIHUAAR6YhzR3FMIhIOPA+cAXQHLhaR7mV224xTzuPNMsc2Ae4DjgcGAveJSGMvYnXN1G83kL2/gDtGdrHFgGrI9ad0pH18Pf7+0QryC4vdDscYr+UXFnP1tBTStufywmX9GdAuMJfY8SZB3Ad8BrQWkTeA+cAdXhw3EEhX1fWqWoAz8ml06R1UdaOq/gKULdU5AvhcVXep6m7gc2CkF+d0xa79BUz+ej0jeiTQt01A57GQEhMZzkOje7Jh537+++U6t8MxxiuFxSX86Y0lLNq4i6fG9uHULs3cDumovBnF9DlwHs43/beAZFX90ov3bgVsKfVzhmebN7w6VkQmiUiKiKTs2LHDy7f2vf9+kU5eQRG3DbfFgGraHzrHc26flvzfl+tYt2Of2+EYc0wlJcrt7/3M/NVZPDS6J6N6t3Q7pGM6VrnvQ+W9++EMb/0V2Aa08WyrSHntLN4OXPfqWFWdrKrJqprctGlTL9/at7buOcD0HzZxfr9EOic0cCWG2u7us7oTExnG3TOXo2pzI0xgUlUe+DiVD5dt4/YRXbhsUOBXLDrWTIwnPf/GAMnAzzgf3McBPwJ/qOC9M/h9SY5EnATjjQzglDLHfunlsTXq2XlrQOFmWwzINU0bRHPnGd3428zlfLBkK+f3twmKJvA8PW8t0xZuYtKQDlx/SnAU8DxWue9TPeW+NwH9PN/U+wN9gXQv3nsR0FlE2otIFDAOmOVlXHOA4SLS2NM5PdyzLaCkZ+UyY3EG4we3pVWjwBm7XBuNG9Ca/m0b8/Cnq9i9v8DtcIz5nSnfbuDf89cyNrk1d53RNWgGsnjTSd1VVZcf+kFVVwB9KjpIVYuAG3A+2FcB76pqqog8KCKjAERkgIhk4FSIfVFEUj3H7gIewkkyi4AHPdsCyhNz1lA3KiJovg2EsrAw4eExPdl7oJB/zV7tdjjG/GbG4gwe+mQlZ/RsziPn9Qqa5ABerAcBrBKRl4HXcfoBLsP5wK+Qqn4KfFpm272lni/CaT4q79ipwFRvzuOGZVv28Fnqdv5yehJx9aPdDscAXZs35KqT2vPiV+s5v38iA9sH5tBBU3vMSd3OX9//hT90iueZcX0ID7IJtN7cQUwEUoGbgJuBlZ5ttZaq8ujs1cTVi+Kqk6ywbSC5aWhnWjWqw90zl1NQVHb0tDE15/v0nfz5zaX0ahXLi+P7Ex0RfJWdvRnmmq+qT6vqGM/jaVXNr4ngAtW36TtZuD6bG07rFBAVF81hdaMieOjcHqzN2sdL36x3OxxTS/28ZQ/XTE+hfXw9Xp04gHpB+jlh5UYrqaREeeyzNFo1qsMlx7dxOxxTjtO6JnBGz+b8e/5aNmfnuR2OqWXWZuYy4ZWfaFI/iulXDaRR3cCozFoVliAqafaK7SzfmsMtw5KC8paxtrjvnB5Ehofx949W2NwIU2O27Mpj/JSfiAgP4/WrjiehYYzbIVWLJYhKKCwu4Ym5aSQl1Ofcvt5OCjduaB4bw63Dk/hqzQ7+t/xXt8MxtcCO3IOMn/IjeQVFvHbVQNrGBVZl1qrwZj2IJBF5SUTmisiCQ4+aCC7QzFicwYad+7l9RNegG41QG10+uB29WsXywMcr2Ztf6HY4JoTlHCjk8qk/kbn3IK9MHEjX5g3dDsknvLmDeA9nqdF7gNtLPWqV/MJinpm3hn5tGnF6t8AtrmUOCw8THhnTi+x9B3liTprb4ZgQdaCgmKteXUR6Vi4vju9P/7ahU7DTm671IlV9we+RBLhXv99I5t6D/Htc36Ca6FLb9UqM5fLB7Zi2cCPn9UukT+tGbodkQkhBUQnXvbGYJZt389wl/RiS5E5NOH/x5g7iYxG5XkRaiEiTQw+/RxZAcvIK+e8X6ZzSpSnHd4hzOxxTSbcOT6JZg2jufP8X5qRuZ8PO/RSXWMe1qZ7iEuWWd5fxZdoOHhnTizN7hd4yw97cQRxaza10s5ICHXwfTmB68et17M0v4vYRVs47GDWIieThc3tx3RuL+eNriwGIjgijU7P6dEloQOeEBnRpXp/OzRrQqlEdWy7WVEhVufejFXzyy6/cdUZXxg0MzSHvFSYIVa3VU4Wz9uYz9bsNjOrdMiCXBDTeOb17AsvuHU561j7SMnNZsz2XNVn7+H5dNh8s3frbfvWiwumc0ICkhPokJTQgKaEBXZo3oFmDaGtaNL95Ym4ab/y4metO6cgfTw7dWmwVJggRiQSuA4Z4Nn0JvKiqtWJYyL8XrKWoWLnFynkHvXrREfRu3YjeZfohcg4UsjYzlzWZ+1iTmcuazFwWrM7i3ZSM3/ZpGBNBl+aeu42EBnROcO4+rA5X7TP563U8/8U6Lh7YhjtCvFXBmyamF4BI4L+en8d7tl3tr6ACxcad+3n7py2MG9iadvHBP6bZlC+2TiTJ7ZqQXGZd4Ox9B3+XNNZk5vK/X37lzQObf9snvn4UnZs18CSPw01WsXUia/oyTA14Z9FmHvl0NWcd14J/nNsz5O8qvUkQA1S1d6mfF4jIz/4KKJA89fkaIsKFG0/r7HYoxgVx9aMZXD+awR0PD0xQVbJyD7ImM5e07bmszXSarN5L2cL+guLf9mveMOa3hJGU0ICk5g3o3Kx+0NbkMTB7+a/c9cFyTk5qytMXBV9l1qrw5q+1WEQ6quo6ABHpABRXcEzQS92Ww6yft3H9KR1pFuTT5Y3viAgJDWNIaBjDSZ0PD2lUVbbuOfBbwjh0x/HaD5s4WKqqbGLjOkd0jHdqVp+YSCvbEsi+WbuDm95eRt82jXnhsn5ERdSOIhTeJIjbgS9EZD3OkqNtqQXlvh+fk0ZsnciQ7oAyviMiJDauS2Ljupza9fBEyuISZcuuPNIyc1mbmUta5j7WZuby9dodFBY7Q23DBNrG1TuiY7xdXL1a80EUyJZs3s0fX1tMh6b1mHrFAOpG1Z67QG9GMc0Xkc5AF5wEsVpVD/o9Mhf9sD6bL9N2cOcZXa0t2VRLeJjQLr4e7eLrMaJH89+2FxaXsCl7P2nbD/dxpGXmMm9V1tvZ0Z8AAB/fSURBVG9zNCLChA5N6/3WMX4ogbSNq1crmjcCQdr2XCa+soimDaKZftVAYuvWrs+DoyYIETlNVReIyHllXuooIqjqB36OzRWqymOfrSahYTRXDG7ndjgmREWGh9GpWQM6NWvAWRyeYJVfWMz6HftZm+X0cazJ3MfyjBw+Xf4rh4rSRkWE0alpfSdhNG9AkqeT3OZw+Nbm7DzGT/mRmEinMmuzBrWvqflYdxAnAwuAc8p5TYGQTBDzVmWxZPMeHhnTizpR1i5salZMZDjdWzake8vfF3vLKygiPWvf70ZV/bRhFx8u2/bbPnWjwunc7HAzVVLzBnRr0aBWfrBVV9befC6b8iMFxSW8+8fBtG5S1+2QXCEV1coXkfaquqGibW5LTk7WlJSUar1HcYlyxrNfU1iszP3LECLDrf3XBLa9+YWsLTMUd03mPnbkHm4F7p0Yy/AezRnRozmdmtV3MdrgsCevgLEv/sCW3Xm8ec2gkK/fJSKLVTW5vNe86W15H+hXZtsMoL8XJx4JPAuEAy+r6r/KvB4NTPe8VzYwVlU3eibnvew5bwQwXVX/6UWs1fLh0q2sydzHc5f0teRggkLDmEj6t218RAXR3fsLWJOZS8qm3cxN3c7jc9J4fE4aHZo6fSEjejTnuFax1iRVxv6DRUx8dREbdu7nlYkDQj45VORYfRBdgR5AbJl+iIZAhfesIhIOPA8MAzKARSIyS1VXltrtKmC3qnYSkXHAo8BY4EIgWlV7iUhdYKWIvKWqGyt3ed47WFTMU5+voWerhpzZM/SKbpnapXG9KI7vEMfxHeL406md+DXnAJ+vzGRO6nYmf72eF75cR/OGMQzvkcDw7s05vkOTWv+l6GBRMde+vpift+zhv5f258RO8W6H5Lpj3UF0Ac4GGvH7fohc4Bov3nsgkK6q6wFE5G1gNFA6QYwG7vc8nwE8J87URAXqiUgEUAcoAPZ6cc4qe/PHzWzdc4B/ntfLvlWZkNMitg6XD27H5YPbsSevgPmrspi7cjvvpmxh+sJNxNaJZGjXZgzv0ZyTk5rWuv634hLl5reX8c3anTx2wXGM7Nm84oNqgaMmCFX9SEQ+Af6qqo9U4b1bAVtK/ZwBHH+0fVS1SERygDicZDEa+BWoC/xFVXeVPYGITAImAbRpU/VqivsOFvHcgnQGd4jjpM72rcGEtkZ1ozi/fyLn90/kQEExX6/dwZzU7cxflcUHS7cSExnGkM5NGdGjOUO7NaNR3Si3Q/YrVeVvHyxn9ort3HNWNy5Kbu12SAHjmH0QqlosIsOAqiSI8r6Gl+0RP9o+A3Fma7cEGgPfiMi8Q3cjpeKbDEwGp5O6CjECMOWbDWTvL+COkV1CvraKMaXViQr/rU+isLiEnzbsYm7qduauzGTuykzCw4Tj2zdhRI/mDO+RQIvYOm6H7FOqyj9nr+adlC38+bROXH1SrVnFwCvedFJ/LyLPAe8A+w9tVNUlFRyXAZROxYnAtqPsk+FpTooFdgGXAJ95KsZmich3QDKwHh/L3neQl75Zz4geCfRtEzpLBRpTWZHhYZzYKZ4TO8Vz/6ge/JKRw5zU7cxJ3c59s1K5b1ZqqRFRCXRq1sDtkKvtha/WMfnr9Ywf1NYqNpfDm2GuX5SzWVX1tAqOiwDWAEOBrcAi4BJVTS21z5+AXqp6raeT+jxVvUhE/gp0Ba7EaWJaBIxT1V+Odr6qDnPdk1fAs/PXcsnANnROCP4/eGP8IT1rH3NXbmdOaiY/b9kDEPQjot74cRN3z1zB6D4tefqiPkEXv68ca5hrhQmimic+E3gGZ5jrVFV9WEQeBFJUdZaIxACvAX1x7hzGqep6EakPvAJ0x2mGekVVHz/WuXwxD8IYU7HSI6J+WL+L4hKlecMYhnVPYESP4BgR9fHP27jx7aWc2qUZL47vH/Dx+lO1EoSIxAL3cXjBoK+AB1U1x6dRVpMliEpY/yV88hdoMxj6XAptTwDrezFVsCevgAWrs5iTup2v1uwgv7DkdyOihiTFB1xxuy/Tsrh6Wgr92jZm+pUDa30l3eomiPeBFcA0z6bxQG9VLVujyVWWILy0Zi68cxnUbwYH9kBBLjTpAH0ugd6XQGwrtyM0QarsiKicA4XERIZxkmdE1OkBMCIqZeMuLpvyIx2b1uetSYNoGFO7iu+Vp7oJYpmq9qlom9ssQXhh1cfw3kRI6A6XzYTIGGfb0tdh4zcgYdDhVOh7GXQ9CyJsOU1TNYXFJSzasIs5nhFRv+bk/zYianj3BIb3aE7LRjU7Imrltr2MnbyQpvWjeffawcTbcrFA9RPEQuB2Vf3W8/OJwBOqOtjnkVaDJYgKLJ8BH0yCVv3g0hlQp0wJgV0bYNmbzmNvBsQ0guMucpJFi97lv6cxXlBVfsnI+a2TOz1rHwDHJcZ6Orn9PyJqw879XPh/C4kMF2ZcdwKtajg5BbLqJog+OM1LsTgdxruAK441osgNliCOYekb8NGfnL6GS96B6GP8z1hSDBu+cu4qVn0CxQchoRf0vRR6XQT14o5+rDFeWLdjn2f47O9HRA3v7iSL3omNfDqiaHtOPue/8D0HCot594+DrWBhGT4ZxSQiDQFU1a8lL6rKEsRRLHoZ/ner03Q07k2IqkTZ4gO7nTuPZW/AtqUQFgldzoC+46HjaRAeWJ2PJvhsz8nnc8+dxQ/rsykqURIaRnuSRfVHRO3aX8BFLy5ke04+b10ziF6JsT6MPjRU9w4iDmcU0x9wZjl/izOKKdvXgVaHJYhyLHwe5vwNkkbChdOcPoeqykx17kR+eQfydkL95tDnYuhzGcR38l3MptbKyStk/urM342IahgTwdBuCYzokcCQpKaVGhG172ARl770A6u25zL9yoEM6mB3v+WpboL4HPgaeN2z6VLgFFU93adRVpMliDK+fgIWPATdRsH5UyDCR6NHigpg7RwnWaydC1oMrQc5TVA9xhy7+coYLx0oKOabtTuYk5rJ/NWZ7MkrJDoijCFJTRnePYHTuyXQuN7R/6bzC4u58tVF/LhhFy9e1p/TuyfUYPTBpboJYrGq9i+zLeVob+gWSxAeqvDFw/D1406fwbkv+K8pKHe7c0ex9HXYuQYi60L3c52ObZtbYXykyFMjquyIqIHtmjCix5EjooqKS7j+jSXMXZnJ02N7M6ZvoovRB77qJogngBTgXc+mC4AeqnqfT6OsJksQOMlh7j2w8Dnodzmc/QyE1cAkIFXISIGlr8GKD5y5FY3bO3cVNrfC+JCqsnxrzm+d3IdGRPVqFftbspj89XpmLM7g/nO6M+HE9i5HHPiqmyBygXpAiWdTGIeL9qmqNiz3wBpW6xNESQnMvt3plB44CUY+CmEulA8oyINVsw7PrUCcDu2+l0GXM6vXD2JMGYdGRM1NzWSZZ0QUwF9OT+Km0zu7GFnwcK0WU02q1QmipBg+vtH5UD7hRhj2YGA07+zaAD+/5cytyNnizK3odeHhuRWBEKMJGYdGRAFcNqitle73UrUThIiM4nAtpi9V9RMfxucTtTZBFBfBh9fC8vfg5L/CKXcF3gdvSUmpuRUfe+ZW9HQShc2tMMZV1W1i+hcwAHjDs+liYLGq3unTKKupViaIogJ4/0rnQ3fovXDSrW5HVLEDu2HF+84oqG1LSs2tuAw6DrW5FcbUsOomiF+APqpa4vk5HFiqqsf5PNJqqHUJojAf3r3cGXI64p8w+Hq3I6q8zJXOJLyf3z48t6L3OCdZxFv7sTE14VgJwttezNKFe2wqotsK9sNbY53kcPbTwZkcwCkaOOJhuGUVjH0DWvaF7/8DzyXDlOGwZDoczHU7SmNqLW/u5/8JLPWsLCc4fRF3+TUqc3QHc+GNi2DLD84chz6XuB1R9UVEQbeznUduJvzyttMENevPMPuvNrfCGJccs4lJnGEAiUARTj+EAD+q6vaaCc97taKJ6cAeeP18py7S+S9Bz/Pdjsh/Ds2tWPY6LH//8NyKPpc6JT5ibfKTMb7g85nUgSjkE8T+bHjtXMhaBRe+6nzbri0K8jzrVrxWam6FZ92KLmfZ3ApjquFYCcKbJqYfRGSAqi7ycVzGW7mZMH007N4AF78NnQOqDJb/RdWF3mOdx+6Nh9etmHFlqbkVl0KLPtYEZYwPeXMHsRLoAmzEmUEtODOobRRTTcjZCtNHwd5tTnLocLLbEQWGQ3Mrlr3h3F0U5TtzK/pc6ix0VC/e7QiNCQrVbWJqW952Vd3kxYlHAs8C4cDLqvqvMq9HA9OB/kA2MFZVN3peOw54EWiIU+ZjgKrmH+1cIZkgdm+CaedA3i64bAa0GeR2RIHpwB7P3IrXS82tGOlZt8LmVhhzLFVKECISA1wLdAKWA1NUtagSJw0H1gDDgAxgEXCxqq4stc/1wHGqeq2IjAPGqOpYEYkAlgDjVfVnz5oUe1S1+GjnC7kEkb3OSQ4F+2H8B9Aq4LuBAoPNrTCmUqo6D2IakIyTHM4AnqzkeQcC6aq6XlULgLeB0WX2Ge05D8AMYKhn5NRw4BdV/RlAVbOPlRxCTtZqeOUMp9lkwieWHCrj0NyKW1c7K+i16nd4bsXUM2DXercjNCZoHCtBdFfVy1T1RZwS3ydV8r1bAVtK/Zzh2VbuPp67kxwgDkgCVETmiMgSEbmjvBOIyCQRSRGRlB07dlQyvAD16y/w6pmAwIRPoXkvtyMKTuGR0PUsuPgtZyLesIdgxyp4eZgzfNYYU6FjJYjCQ08q07RUSnnDScq2Zx1tnwicJU4v9fw7RkSGHrGj6mRVTVbV5KZNm1YhxACTsRimnQ0RdWDip9Csq9sRhYYGCXDijXDVPIiuD6+eDasCrt6kMQHnWAmit4js9TxygeMOPReRvV68dwbQutTPicC2o+3j6XeIBXZ5tn+lqjtVNQ/4FOjn3SUFqU0LnaGsMY2c5BDX0e2IQk98J7h6PiT0gHcugx9ecDsiYwLaUROEqoarakPPo4GqRpR67s0iQYuAziLSXkSigHHArDL7zAKu8Dy/AFigTq/5HJyEVNeTOE4GVhKq1n8Fr5/nfNO98jNoXO7AMeML9eLhio+d5qfP7oTP7nLW0zDGHMFvS455mqVuwPmwXwW8q6qpIvKgZ30JgClAnIikA7cAd3qO3Q08hZNklgFLVPV//orVVWs/hzcuhMbtYOJsaNjS7YhCX1RduGg6DLoefvivUxW3IM/tqIwJOLainJtWfQzvTYRm3WD8h7Zwjht+eMG5i2jVHy55xybYmVrHF+W+ja8tnwHvXgEt+zhNHpYc3DHoOhj7GmSugJdPh53pbkdkTMCwBOGGZW/CB9c4M6PHz4Q6jSo+xvhPt3Pgik/g4F6YMgw2/+B2RMYEBEsQNW3RFPjwOmh/Mlw6A6IbuB2RAWg9AK6eB3Uaw7RRkDrT7YiMcZ0liJq08L/wv1sgaaRTeC+qrtsRmdKadHCSRMu+8N4EZwZ2iPTRGVMVliBqyjdPwpy7oNsouOg1W8MgUNVtApd/5KxiN/ce+PR2GwZrai0rc+lvqvDFI/D1Y9DrImeZUKsuGtgiY+CCV2BeG/j+35CTARdMgah6bkdmTI2yOwh/UoXP/+4kh77jYcz/WXIIFmFhMPwhOPMJWDsHXj3LWbjJmFrEEoS/lJQ4zRPf/wcGXAPn/BvCwt2OylTWwGucqrA70mDK6c6/xtQSliD8oaQYPr4RFr0EJ/wZznzc+UZqglOXM2DC/6Aw3xkGu/E7tyMypkbYp5avFRfBzD/C0tdgyB1OmWlbJzn4ternjHCqnwCvnetMdDQmxFmC8KWiApgxEZa/B0PvhdPutuQQShq3havmQuJAeP8q+OYpGwZrQpolCF8pzHdKSK+aBSP+CSfd6nZExh/qNHaWgO11Icx/AD652blrNCYE2ZAaXyjIg7cvgfVfwFlPwYCr3I7I+FNENIyZDI3aOPNbcrbCha/YrHgTcuwOoroO5sIbF8CGr5w5DpYcaoewMKcZ8exnYN0CeOVM2Pur21EZ41OWIKrjwB54bYxT3O38l6HPJW5HZGpa8kSnTPiu9U412MzQXdfK1D6WIKpqfzZMOwe2LXMWn+l5vtsRGbd0HuYsE1tSBFNHOCsEGhMCLEFUxb4smHY27FwDF78F3c52OyLjtha9nWGwDVvB6+fDz2+7HZEx1WYJorJytsIrZ8DujXDJu863R2MAGrV21hRvO9iZC/PVYzYM1gQ1SxCVsXuTkxxyM52Ffjqc7HZEJtDUaQSXvg+9L4YvHoaPboDiQrejMqZK/JogRGSkiKSJSLqI3FnO69Ei8o7n9R9FpF2Z19uIyD4Ruc2fcXole50zUiV/j1MOus0gtyMygSoiyhnRdvKdsOx1eONCyN/rdlTGVJrfEoSIhAPPA2cA3YGLRaR7md2uAnaraifgaeDRMq8/Dcz2V4xey1rt3DkUHXCWpkzs73ZEJtCJwKl3wejnYeM3zt9Pzla3ozKmUvx5BzEQSFfV9apaALwNjC6zz2hgmuf5DGCoiFObQkTOBdYDqX6MsWK//gKvnuk8n/AptDjO1XBMkOl7GVz6ntM8+fLpsH252xEZ4zV/JohWwJZSP2d4tpW7j6oWATlAnIjUA/4KPODH+Cq2dbEzWimiDkycDc26uhqOCVIdT3M6rwGmngHp892Nxxgv+TNBlFelruyQjqPt8wDwtKruO+YJRCaJSIqIpOzYsaOKYR7F5h9g2miIaeSMcY/r6Nv3N7VL857OMNjGbZ0+iSWvuR2RMRXyZ4LIAFqX+jkR2Ha0fUQkAogFdgHHA4+JyEbgZuBvInJD2ROo6mRVTVbV5KZNm/ou8vVfOTOkGyQ4dw6N2/ruvU3tFdvK+XvqcDLMugEWPGzDYE1A82eCWAR0FpH2IhIFjANmldlnFnCF5/kFwAJ1nKSq7VS1HfAM8IiqPufHWA9bOw/evAgatXX6HGLLtooZUw0xDZ35M33HO0vRzrzWKRNvTADyWzVXVS3yfOufA4QDU1U1VUQeBFJUdRYwBXhNRNJx7hzG+Sser6z6BN6bAM26wfgPoV6cq+GYEBUeCaP+43wJ+eIfsHcrjH3dmUNhTAARDZFb3OTkZE1JSan6G6x4H96/Blr2hcvet/9ZTc34+R346E8Q18kZ7dSodcXHGONDIrJYVZPLe81mUgMsexPev9qZ/Hb5h5YcTM3pPdZZgGjvNnh5qFP80ZgAYQlixxr48HpoPwQunWGLvpia134IXDUHwqOc2fpr5rodkTGAJQhomgTj3oCL34Goum5HY2qrZt2cYbDxneCtcZAy1e2IjLEEAUDXsyAyxu0oTG3XoLkzcq7TUPjkLzDvfigpcTsqU4tZgjAmkETXh3FvQf+J8O3T8MHVUHTQ7ahMLeW3Ya7GmCoKj4Czn4bG7WDefZC73RkGW7eJ25GZWsbuIIwJRCLwh5vh/CmQsQimDHcWqTKmBlmCMCaQ9brAWX9k/w6nGuzWxW5HZGoRSxDGBLq2J8BVn0NkXXj1bFj9qdsRmVrCEoQxwaBpkjMMtmlXeOdS+OkltyMytYAlCGOCRf1mMOETSBoJn94Gc+62YbDGr2wUkzHBJKqeM6Lpszth4XOQswXGvAiRddyOzNSU4kLITofMVMhc4fxbrxmc+7zPT2UJwphgExYOZzzmVIOde48zDHbcW1Z9OBTtyzqcBA4lhB1pUOwpER8WCU27QFxnv5zeEoQxwUgETrgBYhPhg0kwZZhTDdZWPgxORQedD/7SdwWZqbA/6/A+9ZtDQg9nCduEns7zuM4QEeW3sCxBGBPMepwLDVo49ZumDIOL34bWA92OyhyNqlO5t2wi2LkGtNjZJyLGGYzQebiTBA496sXXeLiWIIwJdm2Od0Y4vXEBTDsHznsJuo9yOypTkAc7VpVqHvIkhQO7D+8T28b58O96licR9IQmHZzZ9AEgMKIwxlRPXEdnrsRb4+Ddy2HEIzD4erejqh1UYc+m3yeBzFTIXgd4FmSLrAcJ3aH76MPNQ826B/zaM5YgjAkV9eLhio+dPok5dzkfWiMecTq1jW/k74WsVZ4kcKiJaCUU5B7ep0kHJwH0uvBw81CjdhAWfLMKLEEYE0oi68CF0+DzvzvDYLcvh+bHQUzssR/RDYPyA8xvSoph14YjRxDt2XR4n+hY58O/97jDzUPNujkVeUOEJQhjQk1YGIx42KkG+/2/nSRxcG8FB4mTJCpKJKGYYPJ2QdZKJwlsX+78m7UKig44r0uYM1qoVX/od/nhJqLYRGc0WQjza4IQkZHAs0A48LKq/qvM69HAdKA/kA2MVdWNIjIM+BcQBRQAt6vqAn/GakzIGXiN8wDnG/HBvZCf4/1jz+bDzw/mVHCyIEgw5U0wy0yFvVsP71OnCTTvCclXHm4eatql1k5E9FuCEJFw4HlgGJABLBKRWaq6stRuVwG7VbWTiIwDHgXGAjuBc1R1m4j0BOYArfwVqzEhLywc6jR2HlVRUgwHc4MnwezbUaaf4CgTzNr9odRQ0p5QPyHk7woqw593EAOBdFVdDyAibwOjgdIJYjRwv+f5DOA5ERFVXVpqn1QgRkSiVdWW1jLGDWHhzoibqo66qckEU3TAKY9+iAsTzEKFPxNEK2BLqZ8zgOOPto+qFolIDhCHcwdxyPnA0vKSg4hMAiYBtGnTxneRG2N8qyYTTFgYNHN3glmo8GeCKO8+TSuzj4j0wGl2Gl7eCVR1MjAZIDk5uex7G2NCRXUTjKkSf/YKZQCtS/2cCGw72j4iEgHEArs8PycCM4HLVXWdH+M0xhhTDn8miEVAZxFpLyJRwDhgVpl9ZgFXeJ5fACxQVRWRRsD/gLtU9Ts/xmiMMeYo/JYgVLUIuAFnBNIq4F1VTRWRB0XkUKGYKUCciKQDtwB3erbfAHQC/i4iyzyPZv6K1RhjzJFENTSa7pOTkzUlJcXtMIwxJqiIyGJVTS7vtSCd+miMMcbfLEEYY4wplyUIY4wx5bIEYYwxplwh00ktIjuATRXueHTx/H4Gd21Q2665tl0v2DXXFtW55raq2rS8F0ImQVSXiKQcrSc/VNW2a65t1wt2zbWFv67ZmpiMMcaUyxKEMcaYclmCOGyy2wG4oLZdc227XrBrri38cs3WB2GMMaZcdgdhjDGmXJYgjDHGlCtkE4SITBWRLBFZUWpbHxH5wVMdNkVEBpY5ZoCIFIvIBaW2fSYie0Tkk5qMvyp8cc2e/ReKSKqI/CIiY2v6OirDR9fcVkQWe/ZPFZFra/o6KsNXf9ue7Q1FZKuIPFdT8VeFD/9/Li5VIbrs8gMBxYfX3EZE5orIKhFZKSLtvA5CVUPyAQwB+gErSm2bC5zheX4m8GWp18KBBcCnwAWltg8FzgE+cfuaauKagSSgs+d5S+BXoJHb1+bna44Coj3P6wMbgZZuX5s/r7nUa88CbwLPuX1dNXHNwD63r8WFa/4SGOZ5Xh+o620MIXsHoapf41mdrvRmoKHneSy/X+Huz8D7QFaZ95kP5PopTJ/yxTWr6hpVXet5vs3zWrmzLAOBj665QA+veR5NgN9Z++pvW0T6Awk4HzoBzVfXHEx8cc0i0h2IUNXPPe+5T1XzvI3Bn2tSB6KbgTki8gTOh8AJACLSChgDnAYMcC88v6jyNXtuX6OAYFvytdLXLCKtcVYx7ATc7kmOwaRS1ywiYcCTwHicu+RgVJW/7RgRSQGKgH+p6oc1GK8vVPaak4A9IvIB0B6YB9ypqsXenCygvyn5wXXAX1S1NfAXnBXtAJ4B/urtf7QgU6VrFpEWwGvARFUtqZFIfafS16yqW1T1OJwEcYWIJNRYtL5R2Wu+HvhUVbfUYIy+VpW/7TbqlKS4BHhGRDrWTKg+U9lrjgBOAm7DSRwdgAlen83tdjY/t+G14/ftdzkcnvshwF7P8w047c4bgX04t2jnljruFIKgD8JX14xzC7sEuNDt66nJ33Op41+hTFt9oD2qe83AG8Bmz/adwF6cb9SuX1sN/p5frQW/50H8vp9iPPC8t+evbU1M24CTcTptTgMOtbW3P7SDiLyKkwyC7dbzaCp1zSISBcwEpqvqezUerW9U9poTgWxVPSAijYETgadqOuhqquzf9oeltk8AklX1ToJLZX/PjYE8VT0oIvE4v+fHajroaqrsNYcDjUWkqaru8Bzj9drMIZsgROQtnG/+8SKSAdwHXAM8KyIRQD4wyYv3+QboCtT3vM9VqjrHb4FXg4+u+SKc0RNxng8OgAmquswvQVeTj665G/CkiCjOt7InVHW5/6KuHl/9bQcTH/6eXxSREpzm9X+p6kr/RV09vrhmVS0WkduA+SIiwGLgJa9j8Nx2GGOMMb9T2zqpjTHGeMkShDHGmHJZgjDGGFMuSxDGGGPKZQnCGGNMuSxBGGOMKZclCGMCiGdikzEBwRKEMVUkIg+JyE2lfn5YRG4UkdtFZJE462k8UOr1Dz3rTqSKyKRS2/eJyIMi8iMwuIYvw5ijsgRhTNVNAa6A36qjjgMygc7AQKAP0F9Ehnj2v1JV+wPJwI0iEufZXg+n3s7xqvptTV6AMccSsqU2jPE3Vd0oItki0hdnXYWlOBUzh3ueg7NAS2fga5ykMMazvbVnezZQjFPH35iAYgnCmOp5Gad8cnNgKs7aCv9U1RdL7yQipwCnA4NVNU9EvgRiPC/na2iWmjdBzpqYjKmemcBInDuHOZ7HlSJSH5yFXESkGc7qX7s9yaErThlmYwKa3UEYUw2qWiAiXwB7PHcBc0WkG7DQKZ7JPuAy4DPgWhH5BUgDfnArZmO8ZdVcjakGT+f0ocWV1rodjzG+ZE1MxlSRZ0H4dGC+JQcTiuwOwhhjTLnsDsIYY0y5LEEYY4wplyUIY4wx5bIEYYwxplyWIIwxxpTr/wE7d0yWiTm/WwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# This makes plots appear in the notebook\n", "%matplotlib inline\n", "\n", "# Plot yearly proportion of deaths at the two clinics\n", "ax = yearly1.plot(x='year', y='proportion_deaths', label='clinic 1')\n", "yearly2.plot(x='year', y='proportion_deaths', label='clinic 2', ax=ax)\n", "ax.set_ylabel('Proportion deaths')" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "0c9fdbf550" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 4. The handwashing begins\n", "

Why is the proportion of deaths constantly so much higher in Clinic 1? Semmelweis saw the same pattern and was puzzled and distressed. The only difference between the clinics was that many medical students served at Clinic 1, while mostly midwife students served at Clinic 2. While the midwives only tended to the women giving birth, the medical students also spent time in the autopsy rooms examining corpses.

\n", "

Semmelweis started to suspect that something on the corpses, spread from the hands of the medical students, caused childbed fever. So in a desperate attempt to stop the high mortality rates, he decreed: Wash your hands! This was an unorthodox and controversial request, nobody in Vienna knew about bacteria at this point in time.

\n", "

Let's load in monthly data from Clinic 1 to see if the handwashing had any effect.

" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "dc": { "key": "0c9fdbf550" }, "tags": [ "sample_code" ] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " date births deaths proportion_deaths\n", "0 1841-01-01 254 37 0.145669\n", "1 1841-02-01 239 18 0.075314\n", "2 1841-03-01 277 12 0.043321\n", "3 1841-04-01 255 4 0.015686\n", "4 1841-05-01 255 2 0.007843\n" ] } ], "source": [ "# Read datasets/monthly_deaths.csv into monthly\n", "monthly = pd.read_csv('./dataset/monthly_deaths.csv', parse_dates=['date'])\n", "\n", "# Calculate proportion of deaths per no. births\n", "monthly['proportion_deaths'] = monthly['deaths'] / monthly['births']\n", "\n", "# Print out the first rows in monthly\n", "print(monthly.head())" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "2da2a84119" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 5. The effect of handwashing\n", "

With the data loaded we can now look at the proportion of deaths over time. In the plot below we haven't marked where obligatory handwashing started, but it reduced the proportion of deaths to such a degree that you should be able to spot it!

" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "dc": { "key": "2da2a84119" }, "tags": [ "sample_code" ] }, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Proportion deaths')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAD8CAYAAABthzNFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOy9eZxjZZn3/b2yV1JJdW2973s3dLM1IKvA6MDMaKuMDjju+I6POqjz+IwOzszjgqPjMpu+oogPKPKqiIDaw4OjjoDSAjYNNL1BrzTd1Wt17VXZk/v945yTSqWS1MlyUt1V9/fzyaeSk3OSuyqpc53r+l2LKKXQaDQajaYQ12QvQKPRaDRnJtpAaDQajaYo2kBoNBqNpijaQGg0Go2mKNpAaDQajaYo2kBoNBqNpiieyV5Avejo6FCLFy+e7GVoNBrNWcVzzz13WinVWey5KWMgFi9ezNatWyd7GRqNRnNWISKvlnpOh5g0Go1GUxRtIDQajUZTFG0gNBqNRlOUKaNBaDQae6RSKbq6uojH45O9FE0DCQQCzJ8/H6/Xa/sYbSA0mmlGV1cX4XCYxYsXIyKTvRxNA1BK0dPTQ1dXF0uWLLF9nA4xacoST2WIpzKTvQxNHYnH47S3t2vjMI0QEdrb2yv2GrWB0JTl7x/ewYd/8PxkL0NTZ7RxmH5U85nrEJOmLK/0jDCSSE/2MjQazSSgPQhNWQZjKaJJHWLSTC22bdvGo48+mnu8adMmvvSlLzn2fu9973t58MEHqzr2Zz/7Gbt37849vuaaaxpWFKwNhKYsA7G0NhCaSSGTceZ7l06nxxmIjRs3cttttznyfrVSaCAaiQ4xacoyGE+ho9VTl8/95y52Hxus62uunRvhM288p+w+hw4d4oYbbuDSSy/lhRdeYOXKlXz/+99n7dq13HLLLfzqV7/i1ltvZfXq1Xzwgx8kGo2ybNky7rnnHlpbW7nmmms4//zz2bJlC4ODg9xzzz1ccskl9Pb2csstt3Dw4EGCwSB33XUX69ev57Of/SzHjh3j0KFDdHR0sHnzZmKxGJs3b+ZTn/oUsViMrVu38o1vfINXX32VW265he7ubjo7O/nud7/LwoULee9730skEmHr1q2cOHGCr3zlK7z1rW8t+vsppfjIRz7CY489xpIlS8gf7fzcc8/x8Y9/nOHhYTo6Ovje977HnDlz+M53vsNdd91FMplk+fLl3HfffWzbto1Nmzbx29/+ln/6p3/ioYceAuAnP/kJH/7wh+nv7+fuu+/mqquuYteuXbzvfe8jmUySzWZ56KGHWLFiRU2fpaMehIjcICJ7RGS/iIwzzyLyQRHZISLbRGSziKzNe+5T5nF7ROR6J9epKU48lSGZzpJIZ8lk9exyTX3Zs2cPH/jAB9i+fTuRSIRvfvObgJGvv3nzZm6++Wbe/e538+Uvf5nt27ezbt06Pve5z+WOHxkZ4amnnuKb3/wmt9xyCwCf+cxnuOCCC9i+fTtf/OIXefe7353b/7nnnuPnP/85P/zhD7n99tu56aab2LZtGzfddNOYdd166628+93vZvv27bzjHe/gox/9aO6548ePs3nzZh555JGyHsdPf/pT9uzZw44dO/jOd77DU089BRg1KB/5yEd48MEHee6557jlllv4h3/4BwBuvPFGnn32WV588UXWrFnD3XffzeWXX87GjRv56le/yrZt21i2bBlgeEFbtmzhP/7jP3J/kzvvvJOPfexjbNu2ja1btzJ//vyqPxsLxzwIEXEDdwCvB7qAZ0Vkk1Iq31f6oVLqTnP/jcC/ATeYhuJm4BxgLvDfIrJSKaVjHQ1kMJbK3Y8m04QD9gtsNGcHE13pO8mCBQu44oorAHjnO9/J17/+dYDcCXtgYID+/n5e+9rXAvCe97yHt73tbbnj3/72twNw9dVXMzg4SH9/P5s3b85dZV933XX09PQwMDAAGGGkpqamCdf19NNP8/DDDwPwrne9i09+8pO559785jfjcrlYu3YtJ0+eLPkav/vd73j729+O2+1m7ty5XHfddYBhFHfu3MnrX/96wAijzZkzB4CdO3fyj//4j/T39zM8PMz115e+Lr7xxhsBuOiiizh06BAAl112GV/4whfo6urixhtvrNl7AGc9iEuA/Uqpg0qpJHA/8Kb8HZRS+b5tCLAuU98E3K+USiilXgH2m6+naSCD8XwDoW2zpr4Upl1aj0OhUNXH54dyqn3dcu/j9/tz94u9V7n1Wcecc845bNu2jW3btrFjxw5+9atfAYaQ/Y1vfIMdO3bwmc98pmzNgrUOt9tNOm1kGf7lX/4lmzZtoqmpieuvv57HHnvM/i9ZAicNxDzgSN7jLnPbGETkr0XkAPAV4KMVHvsBEdkqIlu7u7vrtnCNwUBMGwiNcxw+fJinn34agB/96EdceeWVY55vaWmhtbWVJ598EoD77rsv500A/PjHPwZg8+bNtLS00NLSwtVXX80PfvADAJ544gk6OjqIRCLj3jscDjM0NFR0XZdffjn3338/AD/4wQ/GrcsOV199Nffffz+ZTIbjx4/z+OOPA7Bq1Sq6u7tzv3cqlWLXrl0ADA0NMWfOHFKpVO53mGit+Rw8eJClS5fy0Y9+lI0bN7J9+/aK112IkwaimLY5zuQqpe5QSi0D/g74xwqPvUsptUEptaGzs+i8C00NDMZG6x90LYSm3qxZs4Z7772X9evX09vby4c+9KFx+9x777184hOfYP369Wzbto1Pf/rTuedaW1u5/PLL+eAHP8jdd98NwGc/+1m2bt3K+vXrue2227j33nuLvve1117L7t27Of/883OGxuLrX/863/3ud1m/fj333XcfX/va1yr+3d7ylrewYsUK1q1bx4c+9KGcYfP5fDz44IP83d/9Heeddx7nn39+Tp/4/Oc/z6WXXsrrX/96Vq9enXutm2++ma9+9atccMEFHDhwoOR7/vjHP+bcc8/l/PPP5+WXXx6jv1SLTOQmVf3CIpcBn1VKXW8+/hSAUuqfS+zvAvqUUi2F+4rIL83XerrU+23YsEHpgUH15WcvHOVvfrwNgAf+x2VcsqRtklekqQcvvfQSa9asmdQ1HDp0iDe84Q3s3LmzquOvueYa/uVf/oUNGzbUeWVTm2KfvYg8p5Qq+od00oN4FlghIktExIchOm8qWFi+ivJnwD7z/ibgZhHxi8gSYAWwxcG1Tnt+9sJR9p0c68aO1SC0B6HRTDccy2JSSqVF5Fbgl4AbuEcptUtEbge2KqU2AbeKyOuAFNAHvMc8dpeIPADsBtLAX+sMJme57eHt/PmF8/nCW9bltg1qDULjEIsXL67aewBDXzgT2LFjB+9617vGbPP7/fzhD3+YpBXVF0cL5ZRSjwKPFmz7dN79j5U59gvAF5xbncYinckST2XpGU6O2Z4vUmsNQqMZz7p169i2bdtkL8MxdKsNDSOmd3B6ODFm+2Asjc9tfEViuuX3lMIp7VFz5lLNZ64NhCanL/SMjPcgZrUY+dYjCW0gpgqBQICenh5tJKYR1sCgQCBQ0XG6F5MmFz4a50HEU8wMB+jqi2mRegoxf/58urq60LVD0wtr5GglaAOhYdj0DobiaeKpDAGvGxg1ECGfR4vUUwiv11vR2EnN9EWHmDRE8wTo3rww00AsRUuTlyafW3sQGs00RBsIDcN5BiI/zDQYSxMJeAj53FqD0GimIdpAaMaEj6xU12xWMRg3PIigDjFpNNMSbSA0YzyIbtODGE6mUQoiTV6COsSk0UxLtIHQjCmCszyIgahRJBcJeAn6tQeh0UxHtIHQ5ArlAl5XToOw+jBFmryEtAeh0UxLdJqrhpFEmqDPTXuzjx7LQJitviNNHpq0SK3RTEu0B6EhmkwT8nvoaPbnqqmtPkwtTV5CPo9utaHRTEO0gdAwnMgQ8rlpD/npHioIMQW8BP1u3axPo5mGaAOhYSRheBCdYV/Og7BafbcEvQS9HhLpLJms7t2j0UwntIHQ5AxEe8hP70jSqIGIpRCBZp+HkN9ovaGFao1meqENhIaRZNoIMTX7yGQV/bEUA7EUYb8Hl0sI+oxcBp3qqtFML7SB0BBNZHIiNRjtNgbjaVqCXgCCPsOD0DqERjO90AZCw3AiTbPfQ3uzDzANRCxFJDDWQGgPQqOZXmgDoTHrIDx0mh5Ez3Ay18kVIOTXISaNZjriqIEQkRtEZI+I7BeR24o8/3ER2S0i20XkNyKyKO+5jIhsM2+bnFzndCGTVdzx+P5cCisYTfmiqQzNfjftY0JMox5EkxVi0iK1RjOtcMxAiIgbuAP4E2At8HYRWVuw2wvABqXUeuBB4Ct5z8WUUuebt41OrXM68WJXP1/95R5+89LJ3LZYKoNSEPR7mNHkxe2S8R6EKVLHtAeh0UwrnPQgLgH2K6UOKqWSwP3Am/J3UEo9rpSKmg+fASqbh6epiKN9MQB6R0Y9CEt4DpkZS20hn6lBpIk0GYZBi9QazfTESQMxDziS97jL3FaK9wO/yHscEJGtIvKMiLzZiQVON471GwaiL29qnNWor9msdWgP+Tg+ECeWymiRWqOZ5jjZrE+KbCtaiisi7wQ2AK/N27xQKXVMRJYCj4nIDqXUgYLjPgB8AGDhwoX1WfUU5qhpIHqjeQbC9AqsWofOsJ+Dp4cBcmmuWqTWaKYnTnoQXcCCvMfzgWOFO4nI64B/ADYqpXLzLpVSx8yfB4EngAsKj1VK3aWU2qCU2tDZ2Vnf1U9BLA+id3i8gWg2jUB7yEeXGYqyPAi/x4WIrqTWaKYbThqIZ4EVIrJERHzAzcCYbCQRuQD4NoZxOJW3vVVE/Ob9DuAKYLeDa50WWCf+MR5EclSDAGhv9qNMP88SqUWEkM+jW35rNNMMx0JMSqm0iNwK/BJwA/copXaJyO3AVqXUJuCrQDPwExEBOGxmLK0Bvi0iWQwj9iWllDYQNVJMgxg2T/ohU2ewqqmBnEgNhg4RS2kPQqOZTjg6MEgp9SjwaMG2T+fdf12J454C1jm5tunGUDzFYNw4wffleRDRRKEH4cs9Z3kQYBgI7UFoNNMLXUk9TTjWHwdgcXuQvmiKrNm6e9gyEJZIne9BBPINhEdrEBrNNEMbiGmCFV46d14LmazKVVNbmUlWS+98DyKS50GE/G6dxaTRTDO0gZgmdJkGYt28FgB6TR1iJJHG73HhcRtfBavdhs/jIuB1545v8nlyNRMajWZ6oA3ENOFYfwyvW1g5OwyM6hDD5rAgi/aQ4UHkh5fAELFjOsSk0UwrtIGYJhztizGnpYmOkOEhWO02oslMLrwEEPC6Cfs9tDSNzV8I6jRXjWbaMaGBEJGQiLjM+ytFZKOIeCc6TnNmcaw/xtwZAVpDxkdnpboOJ9I5gdqiI+wfoz+AkcWkRWqNZnphx4P4HUZfpHnAb4D3Ad9zclGa+nO0P8a8GUHazBBSj2kgosmxISaA9fNbWDsnMmZbUIvUGs20w04dhCiloiLyfuD/VUp9RURecHph04VUJovX7WykL5XJcnIwzrwZAZq8bvweV54GkRlT7wDwtZvHdTUh6PWQSGdJZ7I5QVuj0Uxt7Pyni4hcBrwD+L/mNkcL7KYLyXSWS7/4G775xH5H3+fEQJysgnmtTYgYLb3zs5isKupyWDpFNKW9CI1mumDHQHwM+BTwU7NVxlLgcWeXNT0YiKXoHUny77/ey76TQ469j1UDMXdGEwBtIV9Og4gmxoeYihHUQ4M0mmnHhAZCKfU7pdRGpdSXzccHlVIfdX5pUx9L9E1lFH/30HYy2aLd0GvGavM9L89A9OaluTbbMhB6aJBGM92wk8W0UkTuEpFfichj1q0Ri5vqWGmj158zi+cP93Pf04cceZ9CD6I1aHgQSimiyUzu5F8OPTRIo5l+2NESfgLcCfwfQJ8d6ojVavsdly4ilsrylV/u4XVrZzG/NVjX9znaH6Oj2ZerjG4L+egZSRqic1bZCjHpoUEazfTDjgaRVkp9Sym1RSn1nHVzfGXTgPx50F98y7kA/OPPdqJUfUNNR/vjOe8BDA9iKJ5mIGYUy9kRqZusEJOuhdBopg0lDYSItIlIG/CfIvJhEZljbTO3a2okv1He/NYgn7h+FU/s6ebn28YN3quJo33RnP4A0GYWy1kDhGx5EKZIHdXV1BrNtKHcmeE5jBnS1mzpT+Q9p4ClTi1qulDYavvdly1m04vH+Nx/7uKqFR25xnm1oJTiWH+ca1bNzG1rM9ttdPVFASoSqXU1tUYzfSjpQSilliillgJrzPu5G7C2cUucuhQO63G7hC//+XqGE2luf6Q+A/T6oiliqcwYD6K1wIMIVmQgtAeh0UwX7GgQT9ncpqkQq312fhbRyllhPnzNcn6+7RiPv3yq1KG2KcxgAnLtNo70Wh6EnUI5j7lm7UFoNNOFchrEbBG5CGgSkQtE5ELzdg1Q3zSbacpIIo3HJfg9Yz+GD1+7jLktAX645XDN71FYAwHQFjQMxGHTQNjRIPweFy7RhXIazXSi3JnheuC9wHzg3/K2DwF/7+Capg1WDYKIjNnu97hZ2B6kP292dLUMm3OoI3ntu2eYBiInUvsmNhAiolt+azTTjHIaxL1KqWuB9yqlrs27bVRKPWznxUXkBhHZIyL7ReS2Is9/XER2i8h2EfmNiCzKe+49IrLPvL2nqt/uDKdcFXNLk5fBWO3hnJjZO6kpbzqcz+Mi7Pfkwk92PAjQLb81munGhGcGpdRDIvJnwDlAIG/77eWOExE3cAfweqALeFZENiml8tXXF4ANZrfYDwFfAW4y02g/A2zAyJh6zjy2r7Jf78wmmkyXFIhbmry5OoVaiJsGIlBQ69Aa8jFkhpjsVFJb+2mRWqOZPthptXEncBPwEYyU17cBi8oeZHAJsN/s3ZQE7gfelL+DUupxpVTUfPgMRjgLjPDWr5VSvaZR+DVwg433PKsYSWRKFqlFAl4G4/UzEPkeBIwK1cU0kFIEfR7tQWg00wg7Z4bLlVLvBvqUUp8DLgMW2DhuHnAk73GXua0U7wd+UcmxIvIBEdkqIlu7u7ttLOnMYqRMJ9WWJi/RZIZUJlvTe8RSGdwuGTdzwjIQIb9nnAZSipAeGqTRTCvsGIiY+TMqInOBFLDExnHFzjpFe0iIyDsxwklfreRYpdRdSqkNSqkNnZ2dNpZ0ZjGSzOTaaBdijfwcrDHMFE9lx3kPYLTbAHttNiyafJ5caq5Go5n62DEQj4jIDIyT9/PAIYxw0UR0MdbTmA+M6yEhIq8D/gHYqJRKVHLs2c5IIl2yBsGa8larDhFLZXJN+vKx2m3YFajBMCZR3e5bo5k22BGpP2/efUhEHgECSqkBG6/9LLBCRJYAR4Gbgb/M30FELgC+DdyglMqvCvsl8EURaTUf/zHG0KIpRTmR2kpLrdVAxJMZmnzjrwNazRCTnSpqC0OD0B6ERjNdsCNSB0Xkf4vId8wr/Jki8oaJjlNKpYFbMU72LwEPmBPpbheRjeZuXwWagZ+IyDYR2WQe2wt8HsPIPAvcbm6bUpQTqS0PYjBe2xV7LJUh4CniQZghJjtV1BY6zVWjmV7YuXz8LkbjvsvMx10YMyIemehApdSjwKMF2z6dd/91ZY69B7jHxvrOSjJZRSyVKRniiQTqE2KKpzK5Vt355ERqG0VyFu3NPvpjKWLJ4q+pOfvZ0TXAD7e8yhfevA6Xy17ygmbqYkeDWKaU+gqGOI1SKkZxEVlTAdaVeKkTdEudROrSGsRoFpNdVs0KoxTsO+Xc/GzN5PLzbUf50ZYjnB5OTLyzZspjx0AkRaQJM4tIRJYB+ttTI1YsP1gixBOpm0hdIospZyDsewKr50QAePmENhBTlQPdwwCcGtL/4hp7IabPAP8FLBCRHwBXYPRo0tSANQuiVKuNgNeNz+OquVgunswQiIyfK2FpEJV4EAvbggS8Ll4+rg3EVOVA9wgAp4biQMvkLkYz6UzoQSilfg3ciGEUfoTRGuMJZ5c19bEms5WqgwCrH1ONBiKdKepBtDR5WdoRYtWssO3XcruElbPC7Dk5WNOaNGcm8VSGI+YQqVODZ44HEU2m+esfPs+JgfhkL2XaUa7dt9Xe+0KM1hrHMWoRFprbNDWQmyZXJsQTCXhqbthXSlB2uYTH/vYabrxwfpGjSrN6dpg9OsR01rPtSD+f/vnY+eevnB7BelhJiOmBrUc43BOdeMcq2XtymP+7/ThPHTjt2HtoilPOg/hX83YH8AfgLuA75v2vO7+0qc1EIjXUp2FfKZG6WlbNjnB6OEm3jlGf1Tzy4jG+//SrHOmN5bZZ+gNYIaaJGUmk+eSD2/neU4fqvcQc1gwSLZw3nnLtvq81232/ClxotrS4CLgA2N+oBU5VrJYVZT2IOhiIeJ0NxOrZRkhKexFnN8cGDMOw+/hozev+U8OIwKL2oO0Q03Ez7HPw9PAEe1aP1XDy9HDt81E0lWEni2m1UmqH9UAptRM437klTQ9GCuZRF6OlqbaOrulMllRGFdUgqsUyEC+fcF6H+Ppv9vG3P3nR8feZjhw1h0XtOjb6OR7oHmFBa5AFrUHbIabjpqE5aIrbTpAzENprbTh2DMRLIvJ/ROQaEXmtiHwHozJaUwOWgSgnUkcCtXkQ8bTRCbaeBqK92U9Hs78hqa4vHO5jyyu1F9Arpdje1Z/TfTRwtN+48t+dbyBODbOsM8TMsN92CPG4+TpdfVESaWfasFhDr06PaA+i0dgxEO8DdgEfA/4G2G1u09SAVQdRrpuqlcWULyRWghW7LRwWVCuNEqrjqWxNJ3WlFI+9fJI33fF7Nn7j99zrYJz8bCKeyuTi+ZYHkc0qDp4eZvnMZjojhoGw872zQlVZBa86JFTHtAcxadhp1hcH/t28aerESCKN3+PC4y5toyNNHrLKyHgKm603KiE3Tc7mQCC7rJ4d5r5nXiWTVbgdbMcQT2dyM7UrPjaV4d13b2HLoV7mtzbhc7vo0TFsYFQ3WD07zMsnhugZThBNZoinsizrbGYkmSGZydIfTeUKKkuRn3p6sHuYlRWkTdtFi9STR33PHBrbjCRLDwuyqLVhX26aXJ09iFWzwyTSWQ71OBd3BsODSGayVYUu9pwYYsuhXj5y3XIe/9traA15daNBE2sW+evXzgJg9/HBXAbTspnNzAwbhZV2dIhjA3GWdoaA0SK7epMwQ6U9I0my2eq8aU11aAMxSYwkMhO2ucjNhIhWp0PESowbrZXVs42WG06HmRLm+kcSlRsI62rzj9bMwut2EdLDjnJYArVlIHYdG2T/KdNAdOYbiIlTXY/3x1je2cysiN8xodryIDJZRX8d5rRr7KMNxCQxkkhP2EnV6uhabSaT9Y9VbwOxYlYzLoGXjzubyWR5QNWEmSyRtaPZmnuhhx1ZHO2PIWIY+nkzmth9bJAD3SO0hXy0hXzMjAQAe9XUxwfizJ3RxNKOZsdSXa0LHdBhpkYzoQYhIiuBT2BUU+f2V0pd5+C6pjzRZOlW3xa1Nuyz/rH8dTYQAa+bxR0hxzOZrCysaoRq60TS0WxcDQd9HkZ0iAkwDMTMsB+fx8WaORF2HRugPeRnmRkqshtiGoynGE6kmdMSIJXJ8sj24yilbM84t8sYAzGUcETn0BTHTqe2nwB3YlRRax+9ThjCs00NokoDEU/VP83VYvXsMDuPNsiDqMpAJAn7PbkiwZDPTY9OkwQMDWLujCYAzpkb4Tcvn+TUUII/WzcHMGpzQj73hCEmS6CeM6MJj9vFQCxF70iS9ubxzSFrIZ5nILq1B9FQ7ISY0kqpbymltiilnrNujq9sihNN2ggx1ehBOCVSgxGeONwbzdVz1BulVJ6BqPz37x5O0BkePVEF/R7H1nq2caw/xjzTQKydG0EpGIqnWT6zObfPzEhgQg/CErvntgRyQvXB0/XXIeKpTM4T1NXUjcWOgfhPEfmwiMwRkTbr5vjKpjiGSF3eQIT9HkSq9yAs1zzgrb/UtMI8mbziwAkBIJVRWAkrw1WI1N1DidxJBQwPQs/TNuodjvXHcwbinLmR3HPLOkcNRGfYT/cEGoSVLju7JcCyDuPYg9311yFiyQyzW/x4XKI1iAZjJ8T0HvPnJ/K2KWBp/ZczfTDSXMtf2btcQtjvqTrN1SmRGmBRu3HFeKhnhHPn1X9uQDwvtbUakfr0cCLXFgSMsIn2IOD0SIJkJsu8VsNAzJvRlGsKmW8gZob97Dw6UOplACODSQRmRQK4RPB5XI5kMsVSGYJeD+3NPl0s12DszINYUuRmyziIyA0iskdE9ovIbUWev1pEnheRtIi8teC5jIhsM2+b7P9KZwdRGx4E1NawzzrJ1rNZn8Wi9iDgXPVsfty5mhDT6XEehIdoMlN1VfpU4ZjZGmNui2EgRIS1cyL4Pa6c0QCYGZ44xHR8IM7MsB+v24XbJSxuDzpSCxFPZQn43HQ0+7UH0WDsZDF5gQ8BV5ubngC+rZQq+18rIm6MVuGvB7qAZ0Vkk1Jqd95uhzEGEf1tkZeIKaWmZFPAZNooACvXZsOilqFB8WQGEfDXuZIajCvymWE/hxwKMSVMgR0mDjENxFI5QR8gkc4wGE/T2ZyvQbhJZxWJdNYRg3m2YNVAWCI1wDtes5ALF80YUxU/M+InmswwnEiXnHp4fCDOnJbR11na0cxeB+aVx1MZZkX8poHQGkQjsXPm+BZwEfBN83aRuW0iLgH2K6UOKqWSwP3Am/J3UEodUkptB7LFXmCqYlX0lmvUZ1HLTIhYKkPA46572qHF4vZQYzyIMiGmF4/0c8Htv2LfydETk3US6QiP9SCAaa9DWMJyvrfwhvVz+cT1q8fsl0t1HSydyXRsIMbcGYHc46WdIQ73RElnjH/nenlr1kwT7UE0HjsG4mKl1HuUUo+Zt/cBF9s4bh5wJO9xl7nNLgER2Soiz4jImys47ozHqugtdWWWTyRQfcvvWKr4NLl6sag96Fi7jfgYD6L07//K6RGyCnYeG42Xnx4aWwMBEDT/DlNVh9h9bDCnOZXjaH+MZr+HyAQp1jPDZrFciTCTUorj/XFmR0YNzZKOEOms4khfjOFEmpu+/Qz/8NMdRY+vhFjSGJvbEcRVXy0AACAASURBVPbRM5yc9mHCRmLHQGREZJn1QESWYq8eothlayWf7EKl1AbgL4H/yF9D3lo+YBqRrd3d3RW89OSSa/U9gUgNtXkQ8VTWEYHaYnFHiFNDCUdOuvkidblWG71mbcOh06OezGiR3GijOUvvmYoexFA8xZvu2Mx9zxyacN+jZorrRF7lzEj5YrmBWIpYKlPgQRgi9+5jg7z/e8+y5VAvf6hDu3Zr6FVns59kJlt10oamcuwYiE8Aj4vIEyLyW+Ax4H/ZOK4LWJD3eD7GTGtbKKWOmT8PYugeFxTZ5y5z0t2Gzs5Ouy89htPDCW785u/5xY7jVR1fDblhQTZCTJGm6udSG665c91UnBSq80NMQ2UMUF/UNBB5noxlIMbUQVgexBSspj7SGyOVUew9OXGKqVEkF5hwv4lCTFaKa74GYVVi3/bQdrYc6mXNnAhHeqM1N9iLp7I0mSI16HYbjcROFtNvgBXAR83bKqXU4zZe+1lghYgsEREfcDNgKxtJRFpFxG/e7wCuwJhDUXeavG6eP9zPq73ODV0vJDcLwkaIqaXJSyyVIZmuXKaJJ+s7brSQxWaq66sOhJmsEFMk4GG4TIht1IMYXUN3kRBTzoOooqbiTKerz/ju2kkYOJpXRV2OliYvPo+r5OAga5LcnDxjMyNo9HIaSqT54lvW8ZeXLCCRztZU/ZzOGAkdAU+egdCprg2j5BlKRK5TSj0mIjcWPLVMRFBKPVzuhZVSaRG5Ffgl4AbuUUrtEpHbga1KqU0icjHwU6AVeKOIfE4pdQ6wBvi2iGQxjNiXCrKf6kbQ58bncdHXwDYMw7lpchOfvPOrqfOviO0QS2UcDTFZHsQhBz2IjrC/bIjJ8iBeOT2S6wNU2GYDprYH0WVmJk2kB40k0vRHU2ME6lKICJ3N/pIhpsJ0WYsPX7OMSMDLX1y8gCf2nALgcG+UWZGJvZZi5KYi+lx0hI2Qoc5kahzlLmFfixFOemOR5xRQ1kAAKKUeBR4t2PbpvPvPYoSeCo97Clg30evXAxGhPeTLXYk2AiuLyY5IPToTonIDEU/Zq7WolnDAS0ezzyEPwjQQzf5camYxrCFAg/E0fdEUbSEf3cOJMRlMkJ/FNHUNxOnhJEPxVMnhUrkMJhseBBg6RKl+TMcHYrhdMu47+f9cNVoitSjnYUa5eHF1zRfieS3rdYip8ZQ8eyilPmPevV0p9Ur+cyKyxNFVNZjWYGMNhJXXb0ektlp+VyNUx1JZ2pudzflf1B5ypN2GdeXY0ewr21a8L5qkyesmlspwqMdoWd09lBhTAwGjf+tqZkuc6VghJjBOxqUq249WaiDCpWc8HB+IMyvsLztR0BDDDQ+iWnJjc71uWoM+XKINRCOxo2A+VGTbg/VeyGTSFvLRG22gB1GRSF19R1cr+8NJnKqFSOR5ECNlKqB7R1Ksn2+cEK0Y/OnhRC4cYTGVPYgjfaPN98oZ61xYyLaBKF1Nfbw/zpwJXsfncTG3pYkjNRiI3Nhcrxu3S2gL6VqIRlLSQIjIahH5c6BFRG7Mu70XqC6geIbSFvI1VIMYMSuc7egDLU3Gia0qDyKZocnBLCaAxe1BTgzGbeXgV0J+iCmTVWPqIiyUUvRFk5y3YAYuyTMQBW02YPRvPVU9iMuXtQPlEwaO9kdxuySXoTQRM8N+BmKpMRllFscHYsxpmfg0sKCtqTYPomAqYkezj+4hrUE0inJnj1XAG4AZGDqEdbsQ+Cvnl9Y42kK+hs4KGEmkCXrduMq45xaRGuZSx9POitQAizqMOHMtJ4FixFNZRKA1aPz+Q0WK5QbjaTJZxcywn3mtTbzSE8212Sg0EC6XEPS5p5wHMRBLMRRPs3JWmFkRP6+cLv05dPXFmB0J4HHbu2iYbRoAS7uwUErlJslNxMK2YE3fjdxME59lILQH0UjKaRA/F5FHgL9TSn2xgWtqOK1BH0PxNKlMFq/Nf55aiCbTtsXj3NjRKj2IgIOV1GB4EGBk0KyaXb9JX3GzTYgluI4kMlDw8pbX1xbysbg9xKHTI7kMl2KCfnAKzqW29IcFbU3G36CMB/HS8cGKPqMV5uS2faeGc0VwYKQWJ9JZZtvITFrYFqR7KGF4s1V8F2OpsQ0nO5p9jlXva8ZT9myolMpgNNub0rSZFbd9DdIh7MyCsAh43fg9rooNRNZqTOdxXqQGezn4lRBPG0V+1t+pWD8mSzdqDflY0mGcHIu12bAITcG51FYG0/zWIEs6QiVDTLFkhv2nhsfMf5gIa+bH3oLRspZHMN9GuuyCNuMC4khfdV7EqEhtnKosD0K322gMdi6XnxKRb4jIVSJyoXVzfGUNpC1oGIhGZTKNJNK2aiAsqmm3kUiPdc2doqXJS1vIV/daiHjK6LpqpQIXCzHlPIig4UEMxdPsNZv25bfZsJiKHsSRvJP1ovZQLtW1kJdPDJJVcM5c+7M7Qn4P81ub2HtqbIX2S8eNv/Hq2RMbm4WmgThc5fcjXqhBhP3EU9kp9zmeqdi5jL3c/Hl73jYFXFf/5UwOrSEjjNEwA1FBiAkMHaLShn2F4p6TLGoP1r0WwsrAsuZ2FxOXe/NDTB3GiejZQ0bvn6IexBTUILr6jOZ7LU1elph/g0Ono6ybP9YQ7DpmpApX4kEArJwVHtMpF4xQVbNpPCYiVwtRpQ5RODY3v5raTh2RpjYm/Asrpa5txEImk/aQ8aXrG6muKV6ljCQyRa9wS1GNB9FIA7G4PcTTB3rIZpUt4d0O8VQWvycvxFTMg4iO1SAAth7qA0poEH5P1bM1zlS6+mLMbzWa7+VP+StmIFqavLZO6vmsmNXMk/u6x+hzLx0fZPXssK3PujXopdnvqTrVtVgWExipzIvNBAmNc0wYYhKRFhH5N6trqoj8q4jUf8bkJJLzIBqlQVTqQQQ8lRsI0wX3O5zmCnDpkjZODMa59UfPF02JrIZEOjMmxFRUgxhJ4fO4CPrcLGgL4nYJB0+PjGuzYTE1PYgo81sNz2FxGT1o97EB1s6JVDwbZNWsMKmMynmI2azi5RNDrJljzxMRERbUkMk0XqS2qql1qmsjsHP2uAcYAv7CvA0C33VyUY2m1dIgGvSliyYytorkLGZFAnT1xSrqilkYu3WSmy5ewD/86Rp+sfMEN931TMkGb5UQNzvR5gxE0RBTgragDxHB63blro4L22xYBH2eKVUHoZTKeRBghGFmRwK8UhDuS2eyvHxiqOLwEhghJiDXKbbLnPVg10AALKyhFiJuXeiYUxHbQo1NKJnu2DEQy5RSnzEnwx1USn0OsDWT+mzB63YRCXgamMWUttVmw2LD4jb6oyn2nLQ/zrEwduskIsJfXb2UO995EXtPDHHjt35fc+GcJVIHvMa842Ihpt6RFK2h0VCddQVdKnwX8runVLO+gViK4UR6TNjI0IPGnowPdI+QSGc5Z17lBmJZZzMisMfMZHrphKFlrJljP112YVuw6rbf8bQx08TyfCxNqtyUQU39sGMgYiJypfVARK4ASndPO0tpa1DDPqUUI8nSc36LcZlZJfvMwR7bxzRSg7C4/pzZfOEt53KkN8aB7olnE5TDqoMQEZr9nqInhL5okrbQaGM6qyajVFPDoM8zpdp956e4WizpCI0LMe0yp+2dW0EGk0WTz82itiD7zFnTLx0fRISK6ikWtgWrbvtdWD8R8nkQKT8jRFM/7BiIDwF3iMghEXkV+AbwP5xdVuNpbZCBiKeyZJW9edQW82Y0sbAtyNMHKjAQybGx20ZhCYe1hpmsOggwut4WCzH1jSRz4cH89y6WwQSGBpHMZKuarXEmYhXJ5XsQiztC9Iwkx2S97Tw6SMDrGlPsVgkrZoVzIaaXjg+ypD1U0ffXqoWoJsxkzFUfPU25XEKzr/gFg6b+2BkYtE0pdR6wHlinlLpAKbXd+aU1lka1/LZCHKEKQkwAr1naxh9e6bXtplvdUBttIKwuqrUMiYHREBNYBqJIiCmazMWkwYaBML22eveNmiwsD2JBngdheVGv5rXc2HVsgNWzI2U7r5Zj5axmXjk9QiKd4aXj9gVqi/y235USS43vBtAc8BSt9dDUHztZTO0i8nWMsZ+Pi8jXRKTd8ZU1mNagryEahBXiqESkBiPMNBBLsbtM6+t8LHGvERpEPtbJuWYPIq8TbXPAkxuyZJHOZBmIpcZ4ECtmGvHyUu2sLaM8VXSII71RwgEPLcG8MJtpJC2hWinF7uODVQnUFitnhclkFTuPDnC4N1qR/gC1tf1OFBl6ZVwwTI3P8EzHTojpfqAb+HPgreb9Hzu5qMnAatjndAm/9cWu1IO4bGkHYF+HmAwNAgyDFPZ7ajYQiVQ2l6IbKhJiGoilUIoxHsT81iA/+/AVvPG8uUVfM1jnlt/pTJZ/emQ3L9k02vXGyGAKjtm2qM28Wjd1iCO9MYbi6YoqqAuxMpk2bTNGytupoM6nlrbfxaYiFrtg0DiDHQPRppT6vFLqFfP2TxgdXqcUbSEfyXQ2Ny/aKUZDTJV5ELNbAizpCI0xEEopUpni8fTR/HHn6yAK6QzX1nEzk1W5OcQAYf/4udR9eX2Y8jlvwQx8nuK/c6jOQ4MOdI/wfza/wtvufJrf7e2uy2tWQn6Kq0WTz83CtiAPPHeEA93DOYG6Fg9iaWcIt0t4ZPtxANZU8VrVtv2OFZmrHg54GdIaREOwc/Z4XERuFhGXefsL4P86vbBGY51onNYhrEreUmMhy2HpEJmsYjCe4u3feYa/+PbTRffNDVpxuFlfMTqa/TV5EIn0WIG9WEih16x6bwvar0i3PIh6hZgsI+VxC7d871ke2HqkLq9rB6MGIlq0MvrfbzqfaCLDW+74PT969ghul9TUadfvcbO4PUjPSJJIwMNcG3MgClnQGqzSg8iONxB+rUE0CjsG4n8APwSS5u1+4OMiMiQiZX1rEblBRPaIyH4Rua3I81eLyPMikhaRtxY89x4R2Wfe3mP/V6oO60TjtA5hXflEApX3kXnN0naG4ml+t7ebm7/9DM8c7OXAqeLppLFUBr/HVbfWF5XQGfbXJFJbMwACeSGmwqt+y5C3huwb2txUuTp5EP1R4yT17XdexGXL2vnkg9v5zUsn6/Ladt57JJkpqrdctKiVn/31FcxuCfC7vd2smNlcc7KCFWZaU0U1Nhge8OnhBJkKayESqfFtwrUG0TjsZDGFlVIupZTHvLnMbWGlVElfU0TcwB3AnwBrgbeLyNqC3Q4D78UwQPnHtgGfAS4FLgE+IyKtlfxilWK1/HZ6cJB15VONB3HZUiM34K++v5WDp4e5fFk7g/E06SJhpniV/ffrQWe4Ng8iXtBewYo552dwWYbc6qNlh2CdReqBmLGG+W1B7n7PxYR8bp7cd7ourz0RVghvVomZDAvagjz0ocu58cJ53HzxgprfL99AVMPMSICsgp4KLxwMDWLsaao5oNNcG4Wty1gR2QhcbT58Qin1iI3DLgH2K6UOmq9xP/AmYLe1g1LqkPlc4RnueuDXSqle8/lfAzcAP7Kz3mrIeRBOh5jML3a4Cg9iZiTAqllhjg/E+O77LmZH1wBPHehhMJ4eI9aClT8+OQaio9kYwFTtTOx4gX4S9o+GhizDankQM4JVeBB10pksD6I16MXncbFmTiQX83eanrxOtqUIB7z821+cX5f3swzE2moNhFm8eHIwwUwbg4YsYkW+Q+GA0bY9k1VVp+5q7GEnzfVLwMcwTuy7gY+Z2yZiHpAflO0yt9mhlmOrolEaxFA8jc/tqtrlv+d9F/Prj7+Wixa1McM0av1FwmLxVHZSPQigaqE6F2IyDZwl6OeHmXpHkoR87or+jjkPok7hib5oCp/blcuyWTs3wkvHh6pqKVHxe9swEPXkiuXt/Nm6OVy7emZVx1uezsnBeEXHGXPVx4eYAB1magB2NIg/BV6vlLpHKXUPxpX8n9o4rphpt/ufY+tYEfmA1WW2u7u2LJJIwIPHJQ0wEKmqvAeLeTOacv9s1tVzX3S8YFfsyqtRWAai2jBT3BKpfaMhJhjb8rtvJDkug2kigubfo14exEAsSUvQm4vJr50TYTiRrnp6WiXY8SDqyYygjzvecWHJNiYTYXkQpyr4TuSmIhbxIEAbiEZgNwcyP63VbkJ1F5Af/JwPHKvnsUqpu5RSG5RSGzo7O22+dHFEhNaQ88VyQ/F0TQYiH8uDsGLh+cSLxG4bRa0tmQszsKwQU35qY2EVtR08bhd+j6tuGkR/NEVrXojLqjWwhvM4ieVBtFaQxTWZdOZCTPY9iESJbgBWmFHrEM5j5wzyz8ALIvI9EbkXeA74oo3jngVWiMgSEfEBNwObbK7rl8Afi0irKU7/sbnNUdqCzrfbMDyIygXqYsxoMl6nv5gHUSR/vFHU6kEkimQxwdgQU2EfJruE/PVr2NcXTTKjKa+Se1YzbpewuwEGomckSdjvKVnzcabhdbvoaPZV5EGMFnsWiNS5Cwad6uo0Zb9dYvjOm4HXAA+bt8uUUvdP9MJKqTRwK8aJ/SXgAaXULhG53RS9EZGLRaQLeBvwbRHZZR7bC3wew8g8C9xuCdZO0hryOl8HUUcPojWXmjv+HyWeHh+7bRRWZlHVIabCLKYiU+Wq8SAAgr76tfzuj6bGtLkIeN0s72y23Q6lFvqiyVzm3dlCZzjAqQo8iFiJlvVWyFF3dHWesmcqpZQSkZ8ppS7C/tV//vGPAo8WbPt03v1nMcJHxY69B2NYUcNoD/l5+YSz/9xD8RSdHdV11SwkHPDgEhgoEhaLJcc3OWsUPo+L1qC3epG6oFDOMqj5Iaa+kVR1HkQdW34PxFKsL8iiOmduhKcq6LpbLb1VelCTyayIvyIPovBCwSJcZsqgpr7Y8U+fEZGLHV/JGUBryFv0arye1FODcLmElqbia46nspPmQUBt1dTFCuVgNPsokc4wnEiPmQVhl2Adhwb1RZM5Hchi7dwIJwbjFef7V0rvSJL2BgnU9WJm2F+RBlGqZX1Og9AehOPYMRDXYhiJAyKyXUR2iMiUa/cNhgbRF01WXO1ZCYaBqI8GAYZQ3V9kXnUslZmUPkwWtVRTF4rUVg8l64SQqz+o4gQZ8nnqksUUT2WIp7K0NI39LK06AafDTL1VZHFNNrMilVVTlxqbmwsxaQ3Ccexcyv6J46s4Q2gL+VDKCB04kT6YySqGE/XzIABamrwl6iAmT4MAw0C8cLi/qmNHPQhj/X6PG5/HlYs5WzpRJX2YLII+d02NBC0GYlaR3HgPAmD3sUGuWlFbZl0plFJnrQdhVVPbKZYrpUEEvW5EdIipEZS8xBSRgIj8DfAJjNqHo0qpV61bw1bYQGotlrt/y2E+/sC2ks9bV8D1NBCtQe+4LCalVNE2yY2ko7n6jq7WlaM/L0On2e/JhZhyKZ7VeBB+T11CTFY6dGEl94ygj3kzmhxNdY0mMyTS2bPOg5iZK5az972wLhQKv8fWVDktUjtPuRjEvcAGYAeGF/GvDVnRJGJ5DdXWQjy68wS/3Hmi5POWSxype4hp7HoT6SxKMWkiNRgeRDSZqapqOZ7O4CtoNJg/l7o3Wn2RWNDnrotIbRnlGU3jP8s1cyKOhph6G1wkVy+sAs9TQ/Z0iHIt68O6H1NDKHcpu1YptQ5ARO4GtjRmSZOH9Q/XU2WB14FTw4wkMyV7EA3G6u9BGCGmsR7EZLb6tujMmyxX6eyLRCo7Zg4xWEODjL+fdQVabR1EPTyInIEosoa1cyM89vJJo02EA0a6lhDbZJLfj8kO8TJz1Y2xo9pAOE05DyJ31jFrGqY8tXgQ0WSao/2xssfX0sm1FK1BozFefkfXnGs+iR5ERw39mIoZ2LBpIOKpDPdsfoVVs8J0VFEHEPS5iaeyNSci9JcIMYGR6ppVOJYynTMQZ10dhNVuozIPolioVLf8bgzlDMR5IjJo3oaA9db9ieZAnK1YV6TVaBAHu0dy90sdn5sF0VQ/D8I6QQ3kZTJN1rjRfDprmE1dzEBYLb+/9cQBjvbH+OzGc6qaSxCq09hRK3OsmIFwOpPpbPUgvG4X7SFfBRpE6bnqzQGv1iAaQMkzlVJq8s4uk0TA6ybkc1cVYjrQPTq4p2+kePrdUKL+HkR+w75286Q8mj8+uWmuQFWprvFUdtzaQ34PR16N8a3fHuCN583lsmXtVa3L6ugaTWZq+hz6oyl8HldRIzy/tYlIwOOYUH22ehBgCNV2q6ljZUKl4YCHrgY0RZzunB2NXBrIvNamqrpx5k926y0ZYqq/BlGsYV+sRAVqI2kL+XBJlR5EuogH4fcwEEvhcQl//6erq16X5UHU2vK7P5pkRpO3qBcjIpwzt8U5AxFN4nVLrqL4bGJm2H41dSw1PlnBIuzXInUj0AaigGWdzWO8Absc6B7JFU2VGjrkiIEo0rAvcQaEmNwuoS1UXaprvMiwI+tv9pHrVjCnZfyYTbsEffVp+d0fTZUdVnTuvAgvHR8kVWTa30TsPjbIv/96L0oV10l6h402G9WE2CabWRH71dTxIrMgLLQG0Ri0gShgWWczr/ZESaYr+8fef2qYCxbOQKT02NLBuBGW8Ncxu6hYw75SBUaNptrRo/FUFn9BiOnK5R386brZvP/KJTWtqdlfJw8iNr7NRj7nzmshmc6yv8TM8HLc98yrfO03+/jlruIp09U2KjwTqKSauly7mHDAS9ScKqdxDm0gClg2M0QmqzjcOzLxziaZrOKV0yOsmhU2eiOV8SAidfQegFw30fxq6jMhxAS1GIjxIaarV3byzXdcVHN766C/PmNH+6OpojUQFufOM2ZD7Dxa+QjSHUeNCvR/+dXeovPGe0fOXgORX009EeXaxeSGSOkwk6NoA1HA8k5j9u7+U/YNxJHeKMlMlmUzm2kL+UpqEIOx+s2CsAj7jY6u+SEmS6SezBATGLOpqxkaVGyKWL0ImV5VrbUQE4WYlrSHCPncFesQiXSGPSeGWD07zP5Twzz8wtFx+/SdzQYiVyxn10CU8CCsmRAJ3Y/JSbSBKGBpZwigIh3C2ndZZ7PR8K+MB1FP/QGMtgOF1dRWiCtS5gq3EVgeRKlYeikMDcKZr2bOg6ixmro/Vr7dtsslrJ0bqdiD2HNiiFRG8ZHrVnDe/Ba+9t/7cumeFj1ns4GoYLJcPFW60LBZjx1tCNpAFBDye5jTEhiTlTQRloFY3tlMa6j0VLpa51GXYkZBNfW+k8PMDPvHdRptNJ3NfpKZbK6C3C6lKtHrgV0P4kdbDnPjN39f9LlcJ9cyHgQYI0h3Hx+sKE6+vcswKOvnt/DJG1ZztD/GD/5wOPd8OpN1rJlkI5hVQT+mcg0nwzrE1BC0gSjC8pmVZTLtPzVMR7OflqC37NhSQ4Oo/0m7paBh3/5TQ6yYVZ+hRLVQaeWsRbE6iHoR9NnTIB7YeoTnD/cXbSk92oep/En63HktRJMZXjltP1y5o2uA1qCX+a1NXLG8gyuXd3DH4/tJmEOUrGSEs9VAWPPK7XwnyoWYmovMKdfUH20gimCkuo7YDo0c6B5hmRmaams2ZkoUO9aJEBMYmUxWiCmbVew7NcyKmeG6v0+lLOs0jJR1VWwHpVTROoh64fO48LqlbBZT70iSbUcMobhYKMT6W7dO4EGcO8+oqK4kzLT96ADr5s/IpbC+8zWL6B1JsvOooWX01dCo8EzA57FfTR0rk+Ya1mNHG4I2EEVY1hliOJG29SVWSrH/1DDLZhonw7agj1RGFY2NGiGm+nsQ+SGmYwMxosnMGeFBrJ0ToT3k43f7um0fk8yYnWgdFNib/Z6iQ5YsntzXjWXfjw+MNxBWpfxEIablnc34PS7bBiKeyrD35BDrzQwogAsXzgDghcN9wGgjybOtzUY+MyMBum14EIYnWcqDMKfKaQ/CUbSBKIJ15WsnzNQzkmQglmK5eYzVo7+w3UY6k2UkmXFGgwj6cgZi30ljzStnTb4H4XIJV63oYPO+02RtTxEz0jr9DonUYHy++0+W/myf2NON121cwRczEFbV+kQhJo/bxZo5EXYes2cgLL1i3fxRAzEzEmDejCZeMD2anAdxFrbZsDBGj9rUIHzFvwdhPVWuIThqIETkBhHZIyL7ReS2Is/7ReTH5vN/EJHF5vbFIhITkW3m7U4n11nIctMbsFPkZInZOQ/CnJPcMzL2H2B0WJADHkTQy3AiTSqTZd+pIQBWzJx8DwLgqhWd9IwkbTeuSzSghmPV7DAvnxgsGgbMZhW/3dvNH6+dDcCJIgZidOTpxJ/lufMi7Do6aMtA7sgTqPO5cFErL7xqehBnaaO+fGZF/PY1iBJFpUGfOVVOh5gcxTEDISJu4A6MYUNrgbeLyNqC3d4P9CmllgP/Dnw577kDSqnzzdsHnVpnMTrDfsJ+jy0PYr+VwWSekEcrm8cK1U602bCYERxtt7H35DCdYX/ZKt9GctXKDgB+u9demKlw3KgTrJ4dZjCeLuodbD86QO9IktevnUVHs694iMmmSA1w7twWhhJpW/29tncN0NHsZ3bBOM4LFszg2ECcEwPxmqbpnSksag9xcjBRtlguNxWxRJqriNDs1zMhnMZJD+ISYL9S6qBSKgncD7ypYJ83YUyuA3gQ+CM5AxrMiAhLbWYyHTg1QpPXzRzzn7o9ZGRp9BaEmAZz0+ScCTGBEfowBOozw3sAmBkOsGZOhN/ZNRBp5zvRrpptiMd7TgyNe+6JPacQMSq3Z7cEODEQG7dPfyyJz+OytcbRiuqJPaidRwdYP79lXI+lC0wdYtuRPnpHkoQDHrzuszc6/JqlRifeZw72ltwnNxWxzIVCJODVHgSGl1tYK1MvnPyWzQOO5D3uMrcV3cccSjQAWH2cl4jICyLyWxG5qtgbiMgHRGSriGzt7rYvhNpheWczB2xUU79yepjFHaFcx0kr7FBYLDfqQTgjUoNxZbv/5NAZvs+KqgAAHehJREFUoT/kc/XKDp57tc/WP3MjpuGtmm38fV4uaiC6OW/+DNpCPmZHmoprENEUrcHinVwLWTGrGa9b2DGBUB1Nptl3aoh181rGPbd2bgSf28Xzh/vpHUnSfhZ7DwDnzW+h2e/hqQOnS+6TKDGPOh/Dg9AaxN//dAdvvqN4zU6tOGkgiv33FAZiS+1zHFiolLoA+DjwQxGJjNtRqbuUUhuUUhs6OztrXnA+y2aGODEYn/ALeLg3yqK2YO5xs9+D1y3j2m3khgU5pEEAvHR8kJFkJhfuOlN47YpO0lnF0wd6Jty3ESGmliYvc1sC4ya+9Y4kebGrn2tWGd+lOS0BThRLc42mbIWXAPweN6tnR/jpC138evfJkqnTu48NklXj9QfrNc6dF+GFw4YHcTaHl8AQ7y9Z0lb2+2Cnn5g1RGoislmVaz8z1UikMzx9oIdLlrQ58vpOGoguYEHe4/nAsVL7iIgHaAF6lVIJpVQPgFLqOeAAsNLBtY7DymTKnxRXSDarONIXY2H7qIEQEVqDPnqHCw2ENSzImToIgC2vGC77meZBXLS4lSav21aYKedBODzsaNXs8LgQ0+/2Gumt16yaCcDslgD90dS4k0tfNDlhims+t7/pHMIBL3/1/a28+54tbO/qH7ePVStSzIMAuGBhK9u7Bjg1FD/rPQiAy5e1c/D0CMeLhPAgvyNx6e9Bs82ZEP/x33u58suPjStYPD2c4KHnumxn2J2JPP9qP7FUhqtW1PcC2cLJ/8JngRUiskREfMDNwKaCfTYB7zHvvxV4TCmlRKTTFLkRkaXACuCgg2sdh2UgymUydQ8nSKazLMjzIICiDfucFKmtk9WzhwwDcSZpEGBcAb9maRtP2qiHiDeoE+2q2REOdA+PmdfwxJ5TtIV8uTqEOS2GrlToRQzEUhMWyeVzwcJWfvGxq/j0G9ay7Ug/G7/xe/7s609y3zOv8l87T/DJB1/k64/tY05LINfMbvxrzCCRzrL35HDZHlBnC9ZEwKf2F/ci4jZmmoQDHluFck8d6KFnJMn7vrsl1+XgSG+Ut37rKf7XT17kib2nKl3+GcOT+7rxuITXLD3LPAhTU7gV+CXwEvCAUmqXiNwuIhvN3e4G2kVkP0YoyUqFvRrYLiIvYojXH1RKlVa0HGBRexCfx8VLZdIzD/camSkLixiI8RpE/ceNWoT9Htwu4eRggo5m3xkZgrh6ZSeHeqIc7imfzRNPWyEmZz2INXPCpDIq5yEm01kee/kU166amdOTZpsGovAqt5IQk4XX7eKWK5fw+9uu4/NvOoesgv/9s5188P97jl/sPMFVKzq54x0Xljz+goWtuftncw2ExZrZEVqDXp4qEWayE2IKBybOYkpnsuw8NsBrlrZxbCDOB76/lR1dA7z1zqdygv/PtxUGNs4entx3mgsXtjpyXoEyM6nrgVLqUeDRgm2fzrsfB95W5LiHgIecXNtEeN0u1s6JlBUXrZNdoYFoDfnGGZbBeBq/x1XzPINiiAgzmrz0jCTPiBYbxbhqhZHu+swrPWNCcoVYV471HKpUjFGhepBVs8M8c7CHwXiaPzl3dm4fa3JdYS1EXzRZttV3OSIBL++6bDHvfM0idh4dJJpMc+Gi1gmzkua2BMxpbImzugbCwuUSLlvWztMHTqOUGif4x5M2NAgbIaa9J4eJp7K8/ZKFvOs1i/nrHz7PG7+xmZlhPw988DK+//Sr/OyFo0ST6VyfrrOFnuEEO48N8PHXORd9P3tz5RrA+vkt7Dw6UDJGebg3igjMmzF2BGaxlt9OtdmwsMJMZ0KLjWIs6TCyeSZKHW5EoRzA0o5mPC7JZTL9164TBH1urjQNGZCrR8jPZIqnMiTS2ZrrTESEdfNbuHRpu62UVRHhggWGF3G29mEq5LJlHRwbiPNqEa8yZiPE1Oz3Ektlig5VsrD0nvXzZ/Bn6+fw2Teu5cKFM3joQ5ezenaEjefNJZrM8N8vnX1hpt8f6EEpuGqlM/oDaANRlnXzWhhJZjhYohvnkd4oc1uaxnkFrSEf/bHUmDbPgw5MkxvznuYJa8UZJlBbuF3CovZQWdEf8rOYnP1q+jwulnU2s+fEEJms4le7TnLtqpljDFOTz82MoHeMB5Hr5FqlB1ELVj3EVDEQV5g6xO+LpLta34NyY3PDNmZCvNg1QCTgYbHptb73iiU8/OErcrrhJYvbmB0JsGnb2MFMp4cTZ7x4/eTebmYEvSUTG+qBNhBlsHrilGq2drg3yvzWpnHb20M+lBo7BtSpTq4WVi3EmSZQ57O0I8TBCTyIRonUAKvnhHn5+CAvHO7j9HCC6/PCSxazI4ExHoRVIV9u3KhTXLt6Jh3N/jMuS61alnSEmB0JFNUh9p0aQsSYKVKK5sDELb+3d/WzPq87biEul7Dx/Ln8dm937v/1ke3HuOQL/82mF89cbUIpxZP7TnPF8g7cLudqi7WBKMPyzmYCXlfJdtWHe6Pj9AfIa9g3xkA0JsR0Jp88lnY2c7g3WjYkEE9ncLukIZXCq2aHOTYQ54GtR/C5XVy7aryrbtRCjIrUox5E46/iV84Ks/UfXzcua+5sRUS4fFk7zxzoGXe1/tjLp7hgwYyyCRfW2NFSHkQ8ZYxvPW9B+SvsjefNJZVR/GLnCX67t5v/+eNtZBVsfbWheTEVsf/UMCcG41ydFxJ1Am0gyuBxuzhnbktuiHw+8VSGU0OJogbCEhHz220MxdNEmpzzIM6Z28KaOZEzOvywtDNEKqPo6iue+w5mi2cHO7nms9oUqh9+/ihXLG8vasBntzRxvH/Ugzjca4TIZkZKX9lq7HPVyg56RpJsOTR6Mj45GGd71wB/tGZW2WOtz6uUgdh1bJB0VrF+/oyyr3PO3AhLO0N858mDfPC+51g+M2zqj5XNE28kVm+zKx2qf7DQBmIC1s1rYdex8WMju8zma8Uycqx2G/mT5YbiKcJ+5zyI91+5hF98rGhHkjMGa6jSwdOlw0xOjhstxOrJlM4qbigSXgLDg+gZSeZCX0/uO83MsJ+lHaGGrHGqc8M5c5gR9PK93x/KbXv8ZUMw/qM1M8se2zxBy29LoD5vAgMhImw8by4Hu0eYFfHz/Vsu4ZLFbbx0fLCst+sEmazi7s2vsP/U+DYwFn0jSe787QHOXzBjXIJMvdEGYgLWmWMjC2PnVg1EMXffuorPNxCDMWc1iLOBpR0TV6eXGxJTb+a2BAgHPLgEXlfiatWqhTg1aIiWTx3o4coVHbb6MGkmpsnn5u2XLORXu0/kLrr++6VTzJvRxKoJwqUTjR3d3jXAzLA/9xmW4x2XLuLtlyzgvvdfSmfYz7nzWkiksxyYIKminmSyik/85EU+/8huvvSLl0vud/sju+mPpvjiW9Y5viZtICbA6o1TqEOUqoGA8S2/U5kssVTGUQ3ibKA15KM16C37TxdPZ/A7nMFkISJsWNTK1Ss7aS8hhs7JK5bbfXyQ3pFkrqZDUx/e9ZpFiAj3Pf0q8VSGzfu7ed2amRMa4YmymF7s6ue8BeW9B4vOsJ9/vnF97oLPGhc7UZPFepHJKj754HYefuEoSztC/HZvNwPR8Z7RYy+f5KcvHOXD1y5n7dxx7enqjjYQE7C0s5mgzz3ui3K4N0bQ5y7aFyfgdRPyuXMexLCDbTbONpZ2NpfNZEqUGRLjBN9650Xc+c6LSj6f327jyX1GOuYVy7WBqCdzZzRxwzmz+dGWw/zmpVPEU1mum0B/gDwDUcSDGIynONg9wnlFmh/aYUlHM01ed0XzxKtFKcWnHt7OQ8938TevW8HXbr6AVEbxX7uOj9lvMJ7i7x/eyapZYW69drnj6wJtICbE7RLOmTu+otrKYCp1ldOa127DyT5MZxtLOkIl60oyWcXR/jjBMrnv9SbgdZcNac02q6mPD8TZvL+b1bPDzAxPHLLQVMZ7r1jMYDzNZ/9zFyGf21ZvoSavG5cwru8ZwM7cdD57HkQhbpewdm6EXTbHxdbC0wd6eGBrFx++Zhl/87qVnDsvwuL2IP/54lgD8c+PvsypoThfeet6RzoyFEMbCBusmzeDXccGxghWR3qjzG8tnW7YFvLlxkMOOtiH6WxjaWeI7qFEUWHx33+9l5eOD/K2DfMnYWXFafZ7CPs9HDo9wrOH+rhSew+OsGFRK+fOi9A9lOCqFZ22Wq2ICJcsaePh54+O67i7LVdBXX0R2blzI+w6Zm9cbC3c/+wRWpq8fPSPVgDG7/XG8+by1IHTudGszx7q5UdbDvP+K5fYDpvVA20gbLBufoR4KpsbL6qUKlkDYdEW8uU0iCOmoO1kJfXZQimh+le7TvCNx/dz04YF/MWGBcUOnTRmtwT4r10nSKazY1pxaOqHiPDey5cAcN0E2Uv5/K8/XkX3UIJ7nz6U26aUYvO+0yxqD9ZUr3KOmaDySo9zQnXfSJL/2nmCt1wwb4wnu/G8uWQV/GKH8b37+4d3MG9GE//z9Q2deqANhB3WzTMstjVU/vRwklgqw8K20ilmbUEfvSNJnj7Qw9/+5EUWtDVxbg1XM1OFYqmuB7qH+fgDL7J+fgufe9M5Z1yGkDUXwud2cemS9okP0FTFWy6Yx9duPp83n184eLI0Fy9u47UrO7nztwdynvp3njzIUwd6eMelC2taz7p55Tsp1IOfbTtKMpPlpovHXhStmPX/t3fmUVJVdx7/fLvZoQHZtyAgi+yMEFABNRJWwxEjGdBgRoOC+xaSqMyZGI1jMs4cJRyi4jFiMkPUEJPgNjhqHHE0ILghiKgoi6waQWRv+jd/3FdQdqqhm66q9wp+n3PqSL37qvz0q/fe797fvffdEk5uVcL8tzbwwMLVvL/lS247t2feHyjoAaISdGpWnwa1a/Die2FyyrrDzIFIcUL9WmzavoeLH1pM68Z1+f3U03Oymlyh0b5pPYp0qAVReqCMq/7rdWrVKOLeSf3zNsS1KqQ6qvufeMJhnw3kVI/iInFuv7ZVzq9PG9GNbbv28+DCj3hp1VZ+/sxKzundmsuGdqqWT+cWDahVoyhnAcLMeGTxOvq2C5NcyzO2bxuWrvmcGc+/z5jerY44cTAXeICoBEVF4uLTO/DUso0sWL7pYMroSCmm0jKjS8sGPDb1tEqNxT4eqF2jmK81qXewo3re0vWs3LSDO8b1yvmkn6Ml9dhvTy8lk97tGjGqZysefPkjrvndG3RtWcJd3+lT7ZZozeIiurcqydmM6jfXbeO9zTuYODBzS2dsnzYA1Cou4idje+bE4Uh4UrySXDusCy+s3MItjy9jbN/wwx2uk3p4j5Zs/mIP00Z285ZDOcJD+3aye98B7n5uFae0b1zhTOYkkHog4xk5fqyBc/TcOKIrC1ZsomGdmsy+aEDWUjE92zbiibc2ZFyzIhObv9jDs8s3UaO4iAsquPGnePS1ddSrVXzwflKe9k3rMfXMTvRp25iWFaw0mGs8QFSSWjWKuHtCP8bOfJk5r3xMy4a1D5sO6dqyhNvO7ZVHw8KhU/MGvLr6M379fx+x+Yu9zLzglMT1O6Qztm8bmpfUPvh0Xyd5dG1ZwqwLT6F9k3qHTf1WlV5tGjF30VrW/W13hd9beqCMR5esY97S9byx9tBz25rUr8XInpkrPl/uLWX+Wxv4Vp/WB2eEZ+Lm0d2r9wdUE08xVYFurUqYNjKMIjhcesk5PB2b1WfP/jJmPP8+3+zegoEdc7OebraoU7OYs7pVfmSNEw9jeremV5bXRkjNqH4nw3wIM+OFlZsZNWMh0//4Dnv2lzFtRFeeuW4ovds24qY/vM2WcuuZp7jvxQ/Zte8Ak049Mau+2cYDRBWZPKQT5/RpXeGze5wj0ykayVR6oIwfjTo5ZhvHqZiuLUuoWSyee3fzV7aXHijj6rlv8P05SzhQZsy+qD9PXzuEq8/uQvfWDbl7Qj927z/AtHlv/908irWf7WL2wtWc9w9tj3oiX77IaYCQNErSe5I+kHRThvLakh6NyhdJ6pBWdnO0/T1JI3PpWRWKi8SsC09h6pknxa1SsHSOFjUa379dotevcJw6NYv5/uCOPP76J19ZQOj2J1fw1LKN/GB4VxZcfwYjerb6Spq0c4sGTD+nBy+t2sqcVz7+ynfe8fQKahSJHxdA5ShnfRCSioFZwHBgPfCapPlmtiJtt8nA52bWWdJE4BfABEk9gIlAT6AN8Jykrmb21emSTkHSoqQOcy8dlNcZoY5ztEwb2Y2laz7n5j+8Tc82DVm4aisPv7qGKWd04ppo9nMmJg1qz4srt3D7UyvYuH03PxgRvmfB8s38cGS3ghjZKLPcTCOXdBpwq5mNjN7fDGBmd6btsyDa51VJNYBNQHPgpvR90/er6P83YMAAW7JkSU7+Fsdxjm82bt/NOb98mbo1i9m4fTfDurfkvkn9j7jc5869pdzx9LvMXbT24Boi+8vK+J8bzkzMnB9JS81sQKayXKaY2gLr0t6vj7Zl3MfMSoHtQNNKftZxHCcvtG5Ul3sm9GPD9t10b92QGRP7VWot6Pq1a/Cv5/XmPycPYm9pGas/3cn0Md0TExyORC6HuWY6euWbKxXtU5nPImkKMAWgffvqTat3HMc5HGd0bc6frhxMh2b1qzzPYkiXZiy44QyWf7I98aP20sllC2I9kP6AkXbAhor2iVJMjYC/VfKzmNlsMxtgZgOaN/dJTI7j5Ja+X2tMo7pHN/G1Qe0aDOrUNNFzfsqTywDxGtBFUkdJtQidzvPL7TMf+Kfo3+OBFyx0iswHJkajnDoCXYDFOXR1HMdxypGzFJOZlUq6GlgAFAO/NrPlkm4DlpjZfOBB4LeSPiC0HCZGn10u6TFgBVAKXOUjmBzHcfJLzkYx5RsfxeQ4jlN14hrF5DiO4xQwHiAcx3GcjHiAcBzHcTLiAcJxHMfJyDHTSS1pK7CmGl/RDPg0SzrZwp0qRxKdIJle7lR5kuiVC6cTzSzjRLJjJkBUF0lLKurJjwt3qhxJdIJkerlT5UmiV76dPMXkOI7jZMQDhOM4jpMRDxCHmB23QAbcqXIk0QmS6eVOlSeJXnl18j4Ix3EcJyPegnAcx3Ey4gHCcRzHyYgHCKegUAIfpi/Jr6NKIOnoFlLIMUk8p5LCcXFiS2oc/Tcxf6+kDpJaxu1RHkmdJZ0Ut0c6kvpKGg5gCek0k9Rf0pUAZlYWtw+ApFMl3SWpQdwu6UgaIOkh4OS4XVJI6inpG5Coc6qjpFZxe6STmBtmLpDUWNJvgCcltTezsriDhKQTJN0P/BmYK2lynD4pJDWRNAt4HPhV6uYXN9Hv9QwwTVL3tG1x+TSVNBO4F9gZbYu9Bhqts3I/sAnYJyn2RY+jc2omMJOwMFgSfr96kh4AHgFukfQTSbXj8inn9GfgPkkXSKoblcV6bh3TAQK4HmgB/BW4CuKt7UU1u7uBfWbWF/g3YKqkRnE5RV61Iq9SM+sD3AZcl4ALpxhoCrwIvApcCOE3jPHCeRQYbmYDzezhyCfWGmi0XG9TYIyZ/YeZ7Yt7ga2oJvwgsN/MTgMuBoZC7C2u84CGZtYbuAYYBcQd4McBDaJr73aC04WSiuM+t465ACFpoKTUc0VmAJcBc4EekgZH++T1746cGpvZl8CvgGsBzGwB8CXRhZNvIq8mZrYPuN7MrouKuhJuyHlPCUROJwBEN7k9QF3gPaC5pNOisrxdOJFTi+jtLUD9aPtISVMk5f33S/120ds2wBAz+0TSaEkPSPpeHKnCyKshsBm4xMxujIqaAruifeK4/ppFb/dxaL37voQW19BUjT0mp1qE+wBmtpRwXz4fOD2fTpk4pgKEpLMJrYURkmqZ2edmto5wc/kLMAXyW4NJczonqhEvMTOTVCypPmHZ1zfz5ZPBa7Skmmb2uaQiSROBnwIrgdmSzs/XxZPmNCZq1QC0Bj4ys98BbwBXSPrnVBDJo9NwSXXNbDGwSFIZcANwAvCQpIkxHKfRkuqY2drI6TfAJcBLwDDgRknt8uFUzmssUNPMtqV1TL8NfBtiu/5GRtffy8Cbkp4F7iHcF6YD10hqGoOTgC2EtOCUqHJbBKwHusadEj+mAgRwIvA6cBJptV8z20nIY9eRNAEONsvz7dSr3MVhhBry9jy5VOTVAw5euE+YWXsz+zkhd/zdmJy6R9s+BWpGta1hwHigbSqgxeR0CTDdzEaZ2S+AWyOvfKV10p26RNteBAYBD5jZbwnpy+K08nx7nQxgZvujsr8Aq6ObYz5JOXUGepjZRkIL/hNgqJn9EvgZ0BtoUuG35MapS/R6Hvhv4JvAn6J//xU4K+Z0amEHCEmtyt0kNgDzgA7A6eUO7BrgYWC8pH8BJucix34Ep9NSTlH6ZDCwxsx2SLpK0rhcnQyV8CqKvHam7TcfaAjUicMp2taOcHG/CewAfggUS2qbi5roEZwGSqphZjvM7M603+pJoDFQL9s+lXBKpbeeA5YQAhVmtpyQSvkkF06V8Cp//bWIynMaRA/jdCIwJPr9SgnHZWC0z/OE82xXnp3aE4LCPjN7ApgMjIgC/CvAfklFcfZDFGSAkDRM0kJgFiGnn2IwsIrQ99CHkI44Q5LMbBfhJB1D6AR6ycz2xuUUlZ0KdJD0FKGjalm2T4YqeF0uaWh0AZVJGgk8ASwm3Jjz7dRX0qVAK0LteKSZTQaeBZYSWl75duoDTImOk6JUYeo4LSKe49Rb0uWE1s0thOP2U0mvAluBLdmudBzl9bcJaAmMzqZLFZ16A5dKGgZsAyZJ+jEhJbcK2J7NY1VJp17AlQr9WDujCto4YA6wKPYh1GZWUC9CB+oiQk2pBfA0YVQJwARCjaoT8D7wBSEiQ0ijvAaMT4DTmKjsHuBjYFRCjtVwQp/IVMJN+NyYnXYC3yj3+aIkHKeo7CJCv0jcx2kHMDoqa0+4CY9LyDk1Mu2zZ6aux5jPqSGEjuFvEW7e307AcUrdp04H3gXOy8XvV9VXvvLw1SIt/VEG9AMWm9k8hdESOwkHlKjsJsKol+cIo032Rc20FcDXE+K0MyqbaWbXZ8spC177zKxU0uNmdn9CnCxVW0/7jrid9kVlT1hIB2SFajrtjo7TWmBttpyy4LU39fuZ2f8mwKkBoW+m1MyeJKQH43ZKv0+9wqG+rthJfIpJ0iWEHv3bo01vA/0VJpYsIzRb75J0D2Es/2JC59MVwOqoPKsTTrLg1Dq6aD7MllMWvYrMbGuCnFqS5XHqWTxO2xLklPXjlFSvajp9SEhbJu2ekPX7VFaIuwlzhGZaA0Kv/nWEXv9u0fbmhM7KK6L3dYDPgK+X+3zt48EpqV7uVLhOSfVyp/y+Et2CsDCx7Fozm0HonLw1KvqMMIzunWi/PYTJcKkonGrqZa0TOslOSfVyp8J1SqqXO+WZuCNUFaJ0K0KzLNUZdzMhp9eNMHrjTaDj8e6UVC93KlynpHq5U+5fBbWinKSpwCQzGxq9/3fCTNsi4EcWZk0f905J9XKnwnVKqpc75ZaCCRBRp2CZpHmEqem7gMcIcwd2u1OyvdypcJ2S6uVOuSfRfRDpRAe9HmFM8T8Ca81scZwHPYlOSfVyp8J1SqqXO+WegpgHkcaVhFECwy05HTtJdIJkerlT5UiiEyTTy51ySMGkmOBQ8y1uj3SS6ATJ9HKnypFEJ0imlzvlloIKEI7jOE7+KJg+CMdxHCe/eIBwHMdxMuIBwnEcx8mIBwjHyRKSbpU07TDl4yT1yKeT41QHDxCOkz/GES3v6jiFgI9icpxqIGk68D1gHWEFt6WENcanEBal+YCwsFA/wtoD26PX+dFXzCI89XMXcJmZrcynv+McDg8QjnOUSOpPWBpyEGHS6evAfcBDZvZZtM/PgM1mNlPSHOBJM5sXlT0PXG5m70saBNxpZmfn/y9xnMwU2kxqx0kSQ4E/WljvHEnzo+29osDQmLBWwILyH5TUgLC85O/T1oipnXNjx6kCHiAcp3pkaoLPIawJ/Zaki4GzMuxTBGwzs365U3Oc6uGd1I5z9LwEnCeprqQSYGy0vQTYKKkm8N20/XdEZZjZF8BHkr4DYalJSX3zp+44R8b7IBynGqR1Uq8hrEm8grBA/Y+ibcuAEjO7WNJg4AFgLzAeKAPuJawVUBN4xMxuy/sf4TgV4AHCcRzHyYinmBzHcZyMeIBwHMdxMuIBwnEcx8mIBwjHcRwnIx4gHMdxnIx4gHAcx3Ey4gHCcRzHyYgHCMdxHCcj/w9shw0fRgBSdQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot monthly proportion of deaths\n", "ax = monthly.plot(x='date', y='proportion_deaths')\n", "ax.set_ylabel('Proportion deaths')" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "518e95acc5" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 6. The effect of handwashing highlighted\n", "

Starting from the summer of 1847 the proportion of deaths is drastically reduced and, yes, this was when Semmelweis made handwashing obligatory.

\n", "

The effect of handwashing is made even more clear if we highlight this in the graph.

" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "dc": { "key": "518e95acc5" }, "tags": [ "sample_code" ] }, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Proportion deaths')" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAD8CAYAAABthzNFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOy9eXhjZ3n3/7klWZIlyx6v49m3TCazJjOZ7CQkadKkFBLWElraBGgppaT0xwstS4E2bXkLpewJhEII7VsSQkJCCKEhK5BlSCbbZPY9M57V493WLj2/P845sixL8tFyZHv8fK5L10hH5xw9Y8vnPvf2vUUphUaj0Wg0ubgmewEajUajmZpoA6HRaDSavGgDodFoNJq8aAOh0Wg0mrxoA6HRaDSavGgDodFoNJq8eCZ7AdWira1NLV68eLKXodFoNNOKF1988ZRSqj3fe6eNgVi8eDGbN2+e7GVoNBrNtEJEXi/0ng4xaTQajSYv2kBoNBqNJi/aQGg0Go0mL6dNDkKj0WhySSQSdHV1EY1GJ3spk47f72f+/PnU1dXZPkYbCI1Gc9rS1dVFKBRi8eLFiMhkL2fSUErR09NDV1cXS5YssX2cDjFpihJNpIgmUpO9DI2mLKLRKK2trTPaOACICK2trSV7UtpAaIry6Z++xof/56XJXoZGUzYz3ThYlPNz0AZCU5QDPSN09YUnexkazbTl4MGDrFmzxvb+O3fu5JxzzmH9+vXs27fPwZVNjDYQmqIMRhKE4zrEpNHUigceeIDrr7+el19+mWXLlk24v1KKdDrtyFq0gdAUZSCS1AZCo6mQZDLJjTfeyLp163jnO99JOBzmxRdf5I1vfCPnnnsu11xzDceOHePhhx/ma1/7Gt/73ve44oorAPjKV77CmjVrWLNmDV/72tcAwytZuXIlH/7wh9mwYQOHDx/mV7/6FRdddBEbNmzgXe96F8PDwxWvW1cxaYoyGE2gI7ia04F/+vk2th8drOo5V81t5PNvWT3hfrt27eL73/8+l1xyCe9///u59dZbuf/++/nZz35Ge3s7P/7xj/nMZz7DHXfcwYc+9CEaGhr4+Mc/zosvvsgPfvADfve736GU4oILLuCNb3wjzc3N7Nq1ix/84AfcdtttnDp1in/5l3/hscceIxgM8sUvfpGvfOUrfO5zn6vo/+eogRCRa4GvA27ge0qpf8t5/0PAXwMpYBj4oFJqu/nep4APmO/9jVLqESfXqhlPNJEinjRc11Ra4XZpU6HRlMOCBQu45JJLAHjve9/LF77wBbZu3crVV18NQCqVYs6cOeOOe/rpp3nb295GMBgE4O1vfzu//e1vue6661i0aBEXXnghAJs2bWL79u2Zz4jH41x00UUVr9sxAyEibuBW4GqgC3hBRB60DIDJj5RS3zH3vw74CnCtiKwCbgBWA3OBx0TkTKWUjnXUkMFIIvM8HE8S8ttvsNFophp27vSdIreCKBQKsXr1ap577rmixymlCr5nGQ1rv6uvvpq77rqrsoXm4GQO4nxgr1Jqv1IqDtwNXJ+9g1Iq298LAtZP43rgbqVUTCl1ANhrnk9TQwaj2QZC22aNplwOHTqUMQZ33XUXF154Id3d3ZltiUSCbdu2jTvusssu44EHHiAcDjMyMsL999/PpZdeOm6/Cy+8kGeeeYa9e/cCEA6H2b17d8XrdtJAzAMOZ73uMreNQUT+WkT2AV8C/qbEYz8oIptFZHN3d3fVFq4xGIhoA6HRVIOVK1fywx/+kHXr1tHb28vNN9/Mvffey9///d9z9tlnc8455/Dss8+OO27Dhg3cdNNNnH/++VxwwQX8+Z//OevXrx+3X3t7O3feeSfvec97WLduHRdeeCE7d+6seN1SzIWp6MQi7wKuUUr9ufn6T4HzlVI3F9j/j839bxSRW4HnlFL/z3zv+8DDSqn7Cn3exo0blZ4HUV2e3HmS9935AgAP3fwG1sxrmuQVaTSlsWPHDlauXDnZy5gy5Pt5iMiLSqmN+fZ30oPoAhZkvZ4PHC2y/93AW8s8VuMA2oPQaGY2ThqIF4DlIrJERLwYSecHs3cQkeVZL/8Q2GM+fxC4QUR8IrIEWA487+BaZzwPvHyEPSeGxmwbm4NI1npJGo1mknHMQCilksBHgEeAHcA9SqltInKLWbEE8BER2SYirwAfA240j90G3ANsB/4X+GtdweQsn/zpFu589uCYbYPag9BoZjSO9kEopR4GHs7Z9rms5x8tcuy/Av/q3Oo0FslUmmgiTc9wfMz27BDTSEx7EBrNTENLbWgYMb2DU8OxMdsHI0m8buMrEtGS3xrNjEMbCE0mv9AzMt6DmN3kA2Akpg2ERjPT0AZCkwkfjfMgogk6Qn5EdJJao6k2P/nJT1i5ciVXXHEFr7zyCg8//PDEB9UYbSA0DJvewVA0OWZ63GA0QVN9HUGvRyepNZoq8/3vf5/bbruNJ598siwDkUw6f9Om1Vw1hLMS0L0jcebOqgeMENPyjhD1Xrf2IDSaCnjrW9/K4cOHiUajfPSjH+X48eM8/fTTHDhwgDe96U3cd999RCIRnn76aT71qU/x5je/mZtvvpnXXnuNZDLJP/7jP3L99ddz55138otf/IJoNMrIyAhPPPGEo+vWBkLDcJaBODUcyxiIwUiSRr+HoNetcxCa6c8vPwnHX6vuOTvXwh/824S73XHHHbS0tBCJRDjvvPP49a9/zRNPPMGXv/xlNm7cyNlnn83mzZv51re+BcCnP/1prrzySu644w76+/s5//zzueqqqwB47rnn2LJlCy0tLdX9v+RBGwjNmPCRVeqaTqtMiCmgQ0waTUV84xvf4P777wfg8OHD7Nmzp+j+v/rVr3jwwQf58pe/DEA0GuXQoUMAXH311TUxDqANhIaxHkS3magejidRChrr6wjoEJPmdMDGnb4TPPXUUzz22GM899xzBAIBLr/8cqLRaNFjlFLcd999rFixYsz23/3ud2Nkvp1GJ6k1Y5rgLA9iIGw0yTX66wj4tAeh0ZTLwMAAzc3NBAIBdu7cyaZNm8btEwqFGBoalbq55ppr+OY3v5mZB/Hyyy/XbL3ZaAOhyTTK+etcmVJXS4epsb6OoPYgNJqyufbaa0kmk6xbt47PfvazmSlw2VxxxRVs376dc845hx//+Md89rOfJZFIsG7dOtasWcNnP/vZSVi5DjFpMDyIgNdNa4OXHstARAyD0FjvoV4nqTWasvH5fPzyl78ct/2pp57KPG9paeGFF14Y8/7tt98+7pibbrqJm266qdpLLIj2IDSE40mCPg9tDb5MN7Wlw2T1QWipDY1m5qENhIbhWIqg101r0Ef3UE6IyV9HwOfWYn0azQxEGwgNIzHDg2gPeTMehCX13RSoI1DnIZZMk0o7M31Qo9FMTbSB0GQMRGvQR+9I3OiBiCQQgQavh6DPDWg9Js30xKmxytONcn4O2kBoGIknjRBTg5dUWtEfSTAQSRDyeXC5hIDXqGXQpa6a6Ybf76enp2fGGwmlFD09Pfj9/pKO01VMGsKxFMFWI0kNhtzGYDRJU6AOgIDX8CB0HkIz3Zg/fz5dXV10d3dP9lImHb/fz/z580s6RhsIDcOxJA0+D60NXsA0EJEEjf6xBkJ7EJrpRl1dHUuWLJnsZUxbdIhJY/ZBeGg3PYie4TgDEUOHCSDo0yEmjWYm4qiBEJFrRWSXiOwVkU/mef9jIrJdRLaIyOMisijrvZSIvGI+HnRynTOFVFpx65N7MyWsYIjyhRMpGnxuWseEmEY9iHorxKST1BrNjMIxAyEibuBW4A+AVcB7RGRVzm4vAxuVUuuAe4EvZb0XUUqdYz6uc2qdM4lXu/r590d28fiOE5ltkUQKpSDg8zCrvg63S8Z7EGaSOqI9CI1mRuGkB3E+sFcptV8pFQfuBq7P3kEp9aRSKmy+3ASUlkHRlMSRvggAvSOjHoSVeA6aFUstQa+Zg0jSWG8YBp2k1mhmJk4aiHnA4azXXea2QnwAyBYs8YvIZhHZJCJvdWKBM42j/YaB6DOb4WBUqK/B7HVoDXo5NhAlkkjpJLVGM8NxsopJ8mzLW4wsIu8FNgJvzNq8UCl1VESWAk+IyGtKqX05x30Q+CDAwoULq7Pq05gjpoHoDWcZCNMrsHod2kM+9p8aBsiUueoktUYzM3HSg+gCFmS9ng8czd1JRK4CPgNcp5SKWduVUkfNf/cDTwHrc49VSn1XKbVRKbWxvb29uqs/DbE8iN7h8QaiwTQCrUEvXWYoyvIgfB4XIrqTWqOZaThpIF4AlovIEhHxAjcAY6qRRGQ9cDuGcTiZtb1ZRHzm8zbgEmC7g2udEVgX/jEeRHw0BwHQ2uDDajq1ktQiQtDr0ZLfGs0Mw7EQk1IqKSIfAR4B3MAdSqltInILsFkp9SDw70AD8BMRAThkViytBG4XkTSGEfs3pZQ2EBWSLwcxbF70g2aeweqmBjJJajDyEJGE9iA0mpmEo53USqmHgYdztn0u6/lVBY57Fljr5NpmGkPRBINR4wLfl+VBhGO5HoQ3857lQYBhILQHodHMLHQn9QzhaL8xJH1xa4C+cIK0Kd09bBkIK0md7UH4sw2ER+cgNJoZhjYQMwQrvLRmXhOptMp0U1uVSZakd7YH0ZjlQQR9bl3FpNHMMLSBmCF0mQZi7bwmAHrNPMRILInP48LjNr4KltyG1+PCX+fOHF/v9WR6JjQazcxAG4gZwtH+CHVu4czOEDCahxg2hwVZtAYNDyI7vARGEjuiQ0wazYxCG4gZwpG+CHOa6mkLGh6CJbcRjqcy4SUAf52bkM9DU/3Y+oWALnPVaGYcExoIEQmKiMt8fqaIXCcidRMdp5laHO2PMHeWn+ag8auzSl2HY8lMgtqiLeQbk38Ao4pJJ6k1mpmFHQ/iNxi6SPOAx4H3AXc6uShN9TnSH2HerAAtZgipxzQQ4fjYEBPAuvlNrJrTOGZbQCepNZoZh50+CFFKhUXkA8A3lVJfEpGXnV7YTCGRSlPndjbSl0ilOTEYZd4sP/V1bnweV1YOIjWm3wHg6zeMUzUhUOchlkyTTKUzCW2NRnN6Y+cvXUTkIuBPgF+Y2/So0ioQT6a54AuPc9tTex39nOMDUdIK5jXXI2JIemdXMVld1MWw8hThhPYiNJqZgh0D8VHgU8D9plTGUuBJZ5c1MxiIJOgdifPVR3ez58SQY59j9UDMnVUPQEvQm8lBhGPjQ0z5COihQRrNjGNCA6GU+o1S6jql1BfN1/uVUn/j/NJOf6ykbyKl+Pv7tpBK51VDrxhL5nteloHozSpzbbBlIPTQII1mpmGniulMEfmuiPxKRJ6wHrVY3OmOVTZ6zerZvHSon/9+7qAjn5PrQTQHDA9CKUU4nspc/IuhhwZpNDMPO7mEnwDfAb4H6KtDFbGktv/kgkVEEmm+9Mgurlo1m/nNgap+zpH+CG0N3kxndEvQS89I3Eg6p5WtEJMeGqTRzDzs5CCSSqlvK6WeV0q9aD0cX9kMIHse9BfetgaAf3hgK0pVN9R0pD+a8R7A8CCGokkGIkaznJ0kdb0VYtK9EBrNjKGggRCRFhFpAX4uIh8WkTnWNnO7pkKyhfLmNwf4xDUreGpXNz97ZdzgvYo40hfO5B8AWsxmOWuAkC0PwkxSh3U3tUYzYyh2ZXgRY4a0NVv6E1nvKWCpU4uaKeRKbf/ZRYt58NWj/NPPt3Hp8raMcF4lKKU42h/l8hUdmW0tptxGV18YoKQkte6m1mhmDgU9CKXUEqXUUmCl+TzzAFbVbomnL7nDetwu4YvvWMdwLMktD1VngF5fOEEkkRrjQTTneBCBkgyE9iA0mpmCnRzEsza3aUrEks/OriI6c3aID19+Bj975ShP7jxZ6FDb5FYwARm5jcO9lgdhp1HOY65ZexAazUyhWA6iU0TOBepFZL2IbDAflwPVLbOZoYzEknhcgs8z9tfw4SuWMbfJz4+eP1TxZ+T2QAC0BAwDccg0EHZyED6PC5foRjmNZiZR7MpwDXATMB/4Stb2IeDTDq5pxmD1IIjImO0+j5uFrQH6s2ZHl8uwOYe6MUu+e5ZpIDJJau/EBkJEtOS3RjPDKJaD+KFS6grgJqXUFVmP65RSP7VzchG5VkR2icheEflknvc/JiLbRWSLiDwuIouy3rtRRPaYjxvL+t9NcYp1MTfV1zEYqTycEzG1k+qzpsN5PS5CPk8m/GTHgwAt+a3RzDQmvDIope4TkT8EVgP+rO23FDtORNzArcDVQBfwgog8qJTKzr6+DGw01WL/CvgS8G6zjPbzwEaMiqkXzWP7SvvvTW3C8WTBBHFTfV2mT6ESoqaB8Of0OjQHvQyZISY7ndTWfjpJrdHMHOxIbXwHeDdwM0bJ67uARUUPMjgf2GtqN8WBu4Hrs3dQSj2plAqbLzdhhLPACG89qpTqNY3Co8C1Nj5zWjESSxVsUmv01zEYrZ6ByPYgYDRRnS8HUoiA16M9CI1mBmHnynCxUurPgD6l1D8BFwELbBw3Dzic9brL3FaIDwC/LOVYEfmgiGwWkc3d3d02ljS1GCmipNpUX0c4niKRSlf0GZFECrdLxs2csAxE0OcZlwMpRFAPDdJoZhR2DETE/DcsInOBBLDExnH5rjp5NSRE5L0Y4aR/L+VYpdR3lVIblVIb29vbbSxpajEST2VktHOxRn4OVhhmiibS47wHMOQ2wJ7MhkW915MpzdVoNKc/dgzEQyIyC+Pi/RJwECNcNBFdjPU05gPjNCRE5CrgM8B1SqlYKcdOd0ZiyYI9CNaUt0rzEJFEKiPSl40lt2E3QQ2GMQlruW+NZsZgJ0n9z+bT+0TkIcCvlBqwce4XgOUisgQ4AtwA/HH2DiKyHrgduFYpld0V9gjwBRFpNl//PsbQotOKYklqqyy1UgMRjaeo946/D2g2Q0x2uqgtjByE9iA0mpmCnSR1QEQ+KyL/ad7hd4jImyc6TimVBD6CcbHfAdxjTqS7RUSuM3f7d6AB+ImIvCIiD5rH9gL/jGFkXgBuMbedVhRLUlsexGC0sjv2SCKF35PHgzBDTHa6qC10matGM7Owc/v4AwzhvovM110YMyIemuhApdTDwMM52z6X9fyqIsfeAdxhY33TklRaEUmkCoZ4Gv3VCTFFE6mMVHc2mSS1jSY5i9YGL/2RBJF4/nNqpj+vdQ3wo+df51/fuhaXy17xgub0xU4OYplS6ksYyWmUUhHyJ5E1JWDdiRe6QDdVKUldOAcxWsVklxWzQygFe046Nz9bM7n87JUj3PX8YU4NxybeWXPaY8dAxEWkHrOKSESWAfrbUyFWLD9QIMTTWLUkdYEqpoyBsO8JnDWnEYCdx7WBOF3Z1z0MwMkh/SeusRdi+jzwv8ACEfkf4BIMjSZNBVizIApJbfjr3Hg9roqb5aLxFP7G8XMlrBxEKR7EwpYA/joXO49pA3G6sq97BICTQ1GgaXIXo5l0JvQglFKPAm/HMAp3YUhjPOXssk5/rMlshfogwNJjqtBAJFN5PYim+jqWtgVZMTtk+1xul3Dm7BC7TgxWtCbN1CSaSHHYHCJ1cnDqeBDheJK//tFLHB+ITvZSZhzF5L4tee8NGNIaxzB6ERaa2zQVkJkmVyTE0+j3VCzYVyih7HIJT3z8ct6+YX6eowpzVmeIXTrENO155XA/n/vZ2PnnB06NYL0sJcR0z+bDHOoJT7xjmew+Mcwvthzj2X2nHPsMTX6KeRD/YT5uBX4HfBf4T/P5N5xf2unNRElqqI5gX6Ekdbms6Gzk1HCcbh2jntY89OpR/uu51zncG8lss/IPYIWYJmYkluTv7t3Cnc8erPYSM1gzSHTivPYUk/u+wpT7fh3YYEpanAusB/bWaoGnK5ZkRVEPogoGIlplA3FWpxGS0l7E9ObogGEYth8b7Xnde3IYEVjUGrAdYjpmhn32nxqeYM/ysQQnTw1XPh9FUxp2qpjOUkq9Zr1QSm0FznFuSTODkZx51Ploqq9M0TWZSpNIqbw5iHKxDMTO487nIb7x+B4+/pNXHf+cmcgRc1jUtqOjv8d93SMsaA6woDlgO8R0zDQ0+83kthNkDIT2WmuOHQOxQ0S+JyKXi8gbReQ/MTqjNRVgGYhiSepGf2UeRDRpKMFW00C0Nvhoa/DVpNT15UN9PH+g8gZ6pRRbuvozeR8NHOk37vy3ZxuIk8Msaw/SEfLZDiEeM8/T1RcmlnRGhsUaetWtQ0w1x46BeB+wDfgo8LfAdnObpgKsPohiaqpWFVN2IrEUrNht7rCgSqlVojqaSFd0UVdK8cTOE1x/6zNc961n+KGDcfLpRDSRysTzLQ8inVbsPzXMGR0NtDcaBsLO984KVaUVvO5QotoyED06xFRz7Ij1RYGvmg9NlRiJJfF5XHjchW10Y72HtDIqnkKm9EYpZKbJ2RwIZJezOkP896bXSaUVbgflGKLJVGamdsnHJlL82fef5/mDvcxvrsfrdukLjImVNzirM8TO40P0DMcIx1NEE2mWtTcwEk8RT6XpDycyDZWFyC493d89zJkllE3bRSepJ4/qXjk0thmJFx4WZFGpYF9mmlyVPYgVnSFiyTQHe5yLO4PhQcRT6bJCF7uOD/H8wV5uvvIMnvz45TQH67TQoIk1i/zqVbMB2H5sMFPBtKyjgY6Q0VhpJw9xdCDK0vYgMNpkV21iZqi0ZyROOl2eN60pD20gJomRWGpCmYvMTIhweXmISIFxo5VyVqchueF0mClmrn8kVrqBsO42f2/lbOrcLoJ62FEGK0FtGYhtRwfZe9I0EO3ZBmLiUtdj/RHOaG9gdqPPsUS15UGk0or+Ksxp19hHG4hJYiSWnFBJ1VJ0LbeSyfrDqraBWD67AZfAzmPOVjJZHlA5YSYrydrWYM290MOOLI70RxAxDP28WfVsPzrIvu4RWoJeWoJeOhr9gL1u6mMDUebOqmdpW4Njpa7WjQ7oMFOtmTAHISJnAp/A6KbO7K+UutLBdZ32hOOFpb4tKhXss/6wfFU2EP46N4vbgo5XMllVWOUkqq0LSVuDcTcc8HoY0SEmwDAQHSEfXo+LlXMa2XZ0gNagj2VmqMhuiGkwmmA4lmROk59EKs1DW46hlLI949wuYwzEUMyRPIcmP3aU2n4CfAeji1r76FXCSDzbzEGUaSCiieqXuVqc1Rli65EaeRBlGYg4IZ8n0yQY9LrpGdFJajByEHNn1QOwem4jj+88wcmhGH+4dg5g9OYEve4JQ0xWgnrOrHo8bhcDkQS9I3FaG8aLQ1ZCNMtA6FLX2mInxJRUSn1bKfW8UupF6+H4yk5zwnEbIaYKPQinktRghCcO9YYz/RzVRimVZSBK//93D8doD41eqAI+j2NrnW4c7Y8wzzQQq+Y2ohQMRZOc0dGQ2aej0T+hB2Elu+c2+TOJ6v2nqp+HiCZSGU9Qd1PXFjsG4uci8mERmSMiLdbD8ZWd5hhJ6uIGIuTzIFK+B2G55v666qealpsXkwMOXBAAEimFVbAyXEaSunsolrmogOFB6HnaRr/D0f5oxkCsntuYeW9Z+6iBaA/56J4gB2GVy3Y2+VnWZhy7v7v6eYhIPEVnkw+PS3QOosbYCTHdaP77iaxtClha/eXMHIwy1+J39i6XEPJ5yi5zdSpJDbCo1bhjPNgzwpp51Z8bEM0qbS0nSX1qOJaRBQEjbKI9CDg1EiOeSjOv2TAQ82bVZ0Qhsw1ER8jH1iMDhU4DGBVMIjC70Y9LBK/H5UglUySRIlDnobXBq+U2aoydeRBL8jxsGQcRuVZEdonIXhH5ZJ73LxORl0QkKSLvzHkvJSKvmI8H7f+XpgdhGx4EVCbYZ11kqynWZ7GoNQA41z2bHXcuJ8R0apwH4SEcT5XdlX66cNSUxpjbZBgIEWHVnEZ8HlfGaAB0hCYOMR0biNIR8lHnduF2CYtbA470QkQTafxeN20NPu1B1Bg7VUx1wF8Bl5mbngJuV0oV/asVETeGVPjVQBfwgog8qJTanrXbIYxBRB/Pc4qIUuq0FAWMJ40GsGIyGxaVDA2KxlOIgK/KndRg3JF3hHwcdCjEFDMT7DBxiGkgksgk9AFiyRSD0STtDdk5CDfJtCKWTDtiMKcLVg+ElaQG+JMLF7Jh0awxXfEdjT7C8RTDsWTBqYfHBqLMaRo9z9K2BnY7MK88mkgxu9FnGgidg6gldq4c3wbOBW4zH+ea2ybifGCvUmq/UioO3A1cn72DUuqgUmoLkM53gtMVq6O3mFCfRSUzISKJFH6Pu+plhxaLW4O18SCKhJhePdzP+lt+xZ4Toxcm6yLSFhrrQQAzPg9hJZazvYU3r5vLJ645a8x+mVLXwcKVTEcHIsyd5c+8Xtoe5FBPmGTK+HOulrdmzTTRHkTtsWMgzlNK3aiUesJ8vA84z8Zx84DDWa+7zG128YvIZhHZJCJvLeG4KY/V0VvoziybRn/5kt+RRP5pctViUWvAMbmN6BgPovD//8CpEdIKth4djZefGhrbAwEQMH8Op2seYvvRwUzOqRhH+iM0+Dw0TlBi3REym+UKhJmUUhzrj9LZOGpolrQFSaYVh/siDMeSvPv2TXzm/tfyHl8KkbgxNrct5KVnOD7jw4S1xI6BSInIMuuFiCzFXj9EvtvWUn6zC5VSG4E/Br6WvYastXzQNCKbu7u7Szj15JKR+p4gSQ2VeRDRRNqRBLXF4rYgJ4dijlx0s5PUxaQ2es3ehoOnRj2Z0Sa5UaE5K99zOnoQQ9EE19/6NP+96eCE+x4xS1wn8io7Gos3yw1EEkQSqRwPwkhybz86yAfufIHnD/byuyrItVtDr9obfMRT6YrH8GrsY8dAfAJ4UkSeEpFfA08A/8fGcV3AgqzX8zFmWttCKXXU/Hc/Rt5jfZ59vmtOutvY3t5u99RjODUc4+23PcMvXztW1vHlkBkWZCPE1Fhf/lxqwzV3Tk3FyUR1dohpqIgB6gubBiLLk7EMxJg+CMuDOA27qQ/3RkikFLtPTFxiajTJ+Sfcb6IQk1Ximp2DsDqxP3nfFp4/2MvKOY0c7g1XLLAXTaSpN5PUYFRiaWqDnSqmx4HlwN+YjxVKqSdtnPsFYLmILBERL8JKlsgAACAASURBVHADYKsaSUSaRcRnPm8DLsGYQ1F16uvcvHSon9d7nRu6nktmFoSNEFNTfR2RRIp4svQ0TTRe3XGjuSw2S11fdyDMZIWYGv0ehouE2EY9iNE1dOcJMWU8iDJ6KqY6XX3Gd9dOwcCRrC7qYjTV1+H1uAoODrImyc3JMjazAoaW01AsyRfetpY/Pn8BsWS6ou7nZMoo6PB7sgyELnWtGQWvUCJypVLqCRF5e85by0QEpdRPi51YKZUUkY8AjwBu4A6l1DYRuQXYrJR6UETOA+4HmoG3iMg/KaVWAyuB20UkjWHE/i2n+qlqBLxuvB4XfTWUYRjOTJOb+OKd3U2dfUdsh0gi5WiIyfIgDjroQbSFfEVDTJYHceDUSEYHKFdmA05vD6LLrEyaKB80EkvSH06MSVAXQkRob/AVDDHllstafPjyZTT66/ij8xbw1K6TABzqDTO7cWKvJR+ZqYheF20hI2SoK5lqR7Fb2DdihJPekuc9BRQ1EABKqYeBh3O2fS7r+QsYoafc454F1k50/mogIrQGvZk70VpgVTHZSVKPzoQo3UBEE/Z6Lcol5K+jrcHrkAdhGogGX6Y0Mx/WEKDBaJK+cIKWoJfu4diYCibIrmI6fQ3EqeE4Q9FEweFSmQomGx4EGHmIQnpMxwYiuF0y7jv555eOtkgtyniYYc5bXJ74QjRLsn5UbkN7ELWi4NVDKfV58+ktSqkD2e+JyBJHV1VjmgO1NRBWXb+dJLUl+V1OojqSSNPa4GzN/6LWoCNyG9adY1uDt6iseF84Tn2dm0gixcEeQ7K6eyg2pgcCRn/W5cyWmOpYISYwLsaFOtuPlGogQoVnPBwbiDI75Cs6UdBIhhseRLlkxubWuWkOeHGJNhC1xE4G87482+6t9kImk5agl95wDT2IkpLU5Su6WtUfTuJUL0Qsy4MYKdIB3TuSYN1844JoxeBPDccy4QiL09mDONw3Kr5XzFhnwkK2DUThbupj/VHmTHAer8fF3KZ6DldgIDJjc+vcuF1CS1D3QtSSggZCRM4SkXcATSLy9qzHTUB5AcUpSkvQW9McxIjZ4WwnP9BUb1zYyvIg4inqHaxiAljcGuD4YNRWDX4pZIeYUmk1pi/CQilFXzjO2Qtm4ZIsA5EjswGjP+vT1YO4eFkrULxg4Eh/GLdLMhVKE9ER8jEQSYypKLM4NhBhTtPEl4EFLfWVeRA5UxHbGrx0D+kcRK0odvVYAbwZmIWRh7AeG4C/cH5ptaMl6K3prICRWJJAnRtXEffcorGCudTRpLNJaoBFbUacuZKLQD6iiTQi0Bww/v9DeZrlBqNJUmlFR8jHvOZ6DvSEMzIbuQbC5RICXvdp50EMRBIMRZOcOTvE7EYfB04V/j109UXobPTjcdu7aeg0DYCVu7BQSmUmyU3EwpZARd+NzEwTs8igPaQ9iFpSLAfxMxF5CPh7pdQXarimmtMc8DIUTZJIpamz+cdTCeF40nbyODN2tEwPwu9gJzUYHgQYFTQrOqs36StqyoRYCdeRWApyTm95fS1BL4tbgxw8NZKpcMmX0A+chnOprfzDgpZ642dQxIPYcWywpN/RcnNy256Tw5kmODBKi2PJNJ02KpMWtgToHooZ3mwZ38VIYqzgZFuDzzGJec14il4NlVIpDLG905oWs+O2r0Z5CDuzICz8dW58HlfJBiJtCdN5nE9Sg70a/FKIJo0mP+vnlE+PycobNQe9LGkzLo75ZDYsgqfhXGqrgml+c4AlbcGCIaZIPMXek8Nj5j9MhDXzY3fOaFnLI5hvo1x2QYtxA3G4rzwvYjRJbVyq2hq8nBqOabmNGmHndvlZEfmWiFwqIhush+MrqyEtAcNA1KqSaSSWtNUDYVGO3EYsOdY1d4qm+jpagt6q90JEE4bqqlUKnC/ElPEgAoYHMRRNstsU7cuW2bA4HT2Iw1kX60WtwUypay47jw+SVrB6rv3ZHUGfh/nN9ew+ObZDe8cx42d8VufExmahaSAOlfn9iI7LQfiIJtKn3e9xqmLnNvZi899bsrYp4MrqL2dyaA4aYYyaGYgSQkxg5CFKFezLTe45yaLWQNV7IawKLGtud77kcm92iKnNuBC9cNDQ/snrQZyGOYiuPkN8r6m+jiXmz+DgqTBr5481BNuOGqXCpXgQAGfODo1RygUjVNVgGo+JyPRClJmHyB2bm91NbaePSFMZE/6ElVJX1GIhk0lr0PjS9Y2UJ4pXKiOxVN473EKU40HU0kAsbg3y3L4e0mllK/Fuh2gijc+TFWLK50GEx+YgADYf7AMK5CB8nrJna0xVuvoizG82xPeyp/zlMxBN9XW2LurZLJ/dwG/3dI/Jz+04NshZnSFbv+vmQB0NPk/Zpa7jqphCo81yi80CCY1zTBhiEpEmEfmKpZoqIv8hItWfMTmJZDyIWuUgSvUg/J7SDYTpgvscLnMFuGBJC8cHo3zkrpfylkSWQyyZGhNiypuDGEng9bgIeN0saAngdgn7T42Mk9mwOD09iDDzmw3PYXGRfND2owOsmtNY8myQFbNDJFIq4yGm04qdx4dYOceeJyIiLKigkml8ktqS29CVTLXAztXjDmAI+CPzMQj8wMlF1ZpmKwdRI42XcCxlq0nOYnajn66+SEmqmLmxWyd593kL+MybVvLLrcd593c3FRR4K4WoqUSbMRB5Q0wxWgJeRIQ6tytzd5wrs2ER8HpOqz4IpVTGgwAjDNPZ6OdATrgvmUqz8/hQyeElMEJMQEYptsuc9WDXQAAsrKAXImrd6JhTEdtDxmS5REonqWuBHQOxTCn1eXMy3H6l1D8BtmZSTxfq3C4a/Z4aVjElbclsWGxc3EJ/OMGuE/bHOebGbp1ERPiLy5bynfeey+7jQ7z9289U3DhnJan9dca843whpt6RBM3B0VCddQddKHwX9LlPK7G+gUiC4VhyTNjIyAeNvRjv6x4hlkyzel7pBmJZewMisMusZNpx3MhlrJxjv1x2YUugbNnvaNKYaWJ5Ph0hP5v/4Srecvbcks+lKR07BiIiIm+wXojIJUBh9bRpSkuNBPuUUozEC8/5zcdFZpfspv09to+pZQ7C4prVnfzr29ZwuDfCvu6JZxMUw+qDEBEafJ68Iaa+cJyW4KgwndWTUUjUMOD1nFZy39klrhZL2oLjQkzbzGl7a0qoYLKo97pZ1BJgjzlresexQUQoqZ9iYUugbNnvcvsnNNXBjoH4K+BWETkoIq8D3wL+0tll1Z7mGhmIaCJNWtmbR20xb1Y9C1sCPLevBAMRHxu7rRVW4rDSMJPVBwGG6m2+EFPfSDwTHsz+7HwVTGDkIOKpdFmzNaYiVpNctgexuC1Iz0h8TNXb1iOD+OtcY5rdSmH57FAmxLTj2CBLWoMlfX+tXohywkzGXHXn82ia/NgZGPSKUupsYB2wVim1Xim1xfml1ZZaSX5bIY5gCSEmgAuXtvC7A7223XRLDbXWBsJSUa1kSAyMhpjAMhB5QkzhOC3BEgyE6bVVWzdqsrA8iAVZHoTlRb2eJbmx7egAZ3U2FlVeLcaZsxs4cGqEWDLFjmP2E9QW2bLfpRJJOK8GoCmMnSqmVhH5BsbYzydF5Osi0ur4ympMc8BbkxyEFeIoJUkNRphpIJJgexHp62ys5F6t3XPr4lyxB5GlRNvg92SGLFkkU2kGIokxHsTyDiNeXkjO2jLKp0se4nBvmJDfQ1MgK8xmGkkrUa2UYvuxwbIS1BZnzg6RSiu2HhngUG+4pPwDVCb7HXN46JWmOHZ8t7uBbuAdwDvN5z92clGTgSXY53QLv3WhK9WDuGhpG2A/DzEZOQgwDFLI56nYQMQS6UyJbjBPiGkgkkApxngQ85sDPPDhSwomMANVlvxOptL8y0Pb2WHTaFcbo4IpMGbbohbzbt3MQxzujTAUTZbUQZ2LVcn04CvGSHk7HdTZVCL77fRURE1x7BiIFqXUPyulDpiPf8FQeD2taAl6iSfTmXnRTjEaYirNg+hs8rOkLTjGQCilSKTyx9NH68drH7+tVHEzlVaZOcQAId/4udR9WTpM2Zy9YBbeAjHrYJWHBu3rHuF7Tx/gXd95jt/s7q7KOUshu8TVot7rZmFLgHtePMy+7uFMgroSD2JpexC3S3hoyzEAVpZxrnJlvyMOz1XXFMfO1eNJEblBRFzm44+AXzi9sFpjXWiczkNYnbyFxkIWw8pDpNKKwWiC9/znJv7o9ufy7psZtOKwWF8+2hp8FXkQseTYBLuRgxh7199rdr23BOx3pFseRLVCTJaR8riF99/5AvdsPlyV89rB6IEI5+2M/uq7zyEcS/G2W5/hrhcO43ZJRUq7Po+bxa0BekbiNPo9zLUxByKXBc2BMj2ItDYQk4gdA/GXwI+AuPm4G/iYiAyJSFHfWkSuFZFdIrJXRD6Z5/3LROQlEUmKyDtz3rtRRPaYjxvt/5fKw7rQOJ2HGDLLNRv9pevIXLi0laFokt/s7uaG2zexaX8v+07mLyeNJFL4PK6qSV+UQnvIV1GS2poB4M8KMeXe9VuGvDlo39BmpspVyYPoDxtG6vb3nstFy1r5u3u38PiOE1U5t53PHomn8uZbzl3UzAN/fQmdTX5+s7ub5R0NFV9krTDTyjK6scHwgE8Nx0iV2AsRS+gy18nEThVTSCnlUkp5zIfL3BZSShX0NUXEDdwK/AGwCniPiKzK2e0QcBOGAco+tgX4PHABcD7weRFpLuU/ViqW5LfTg4Mspc1yPIiLlhq1AX/xX5vZf2qYi5e1MhhNkswTZopOYv14e6gyDyKaI69gJamzK7gsQ27paNkhUOUk9UDEWMP8lgDfv/E8gl43v91zqirnnggrhDe7wEyGBS0B7vuri3n7hnnccN6Cij8v20CUQ0ejn7SCnhJvHIwchC5znSxs3caKyHXAZebLp5RSD9k47Hxgr1Jqv3mOu4Hrge3WDkqpg+Z7uVe4a4BHlVK95vuPAtcCd9lZbzlkPAinQ0ymBxEqw4PoaPSzYnaIYwMRfvC+83ita4Bn9/UwGE2OSdaCVT8+OQaircEYwFTuTOxoTv4k5BsNDVmG1fIgZgXK8CCqlGeyPIjmQB1ej4uVcxozMX+n6clSsi1EyF/HV/7onKp8nmUgVpVrIMzmxRODMTpsDBqyiNRgrrqmMHbKXP8N+CjGhX078FFz20TMA7KDsl3mNjtUcmxZ1CoHMRRN4nW7yv7S3/G+83j0Y2/k3EUtzDKNWn+esFg0kZ5UDwLKF1TLhJhMA2cl9LPDTL0jcYJed0k/x4wHUaWhQX3hBF63K1Nls2puIzuODZUlKVHyZ9swENXkkjNa+cO1c7jirI6yjrc8nROD0ZKOM+aqawMxWdjx3d4EXK2UukMpdQfGnfybbByXL1Bp9y/H1rEi8kFLZba7u7Iqkka/B49LamAgEmV5DxbzZtVn/tisu+e+8Pgmssm887IMRLlhpqiVpPaOhphgrOR330h8XAXTRATMn0e1PIiBSJymQF0mJr9qTiPDsWTZ09NKwY4HUU1mBbzc+icbCsqYTITlQZws4TuRmYqoDcSkYTe4l13WaregugvIDn7OB45W81il1HeVUhuVUhvb29ttnjo/IkJz0PlmuaFosiIDkY3lQVix8Gyikxi7zQx1KVMdN7cCywoxDWXpMeV2UdvB43bh87iqloPoDydozgpxWb0G1nAeJ7E8iOYSqrgmk/ZMiMm+BxGbJDUAzSh2riD/F3hZRO4UkR8CLwJfsHHcC8ByEVkiIl7gBuBBm+t6BPh9EWk2k9O/b25zlJaA83IbhgdReoI6H7PqjfP05/MgJrF+vFIPIpanignGhphydZjsEvRVT7CvLxxnVn1WJ/fsBtwuYXsNDETPSJyQz1Ow52OqUed20dbgLcmDGG32nB7/x9ORoj95MXznp4ELgZ+aj4uUUndPdGKlVBL4CMaFfQdwj1Jqm4jcYia9EZHzRKQLeBdwu4hsM4/tBf4Zw8i8ANxiJaydpDlY53wfRBU9iOZMae54AxFNTl7s1qosKjvElFvFlGeqXDkeBEDAWz3J7/5wYozMhb/OzRntDbblUCqhLxzPVN5NF9pDfk6W4EFEaihZr8lP0SuVUkqJyANKqXOxf/efffzDwMM52z6X9fwFjPBRvmPvwBhWVDNagz52Hnf2j3somqC9rTxVzVxCfg8ugYE8YbFIfPJEzrweF82BuvKT1DmNcpZBzQ4x9Y0kyvMgqij5PRBJsC6nimr13EaeLUF1t1x6y/SgJpPZjb6SPIjcGwVN7bHju20SkfMcX8kUoDlYl/duvJpUMwfhcglN9fnXHE2kJ7X6o5Ju6nyNcjBafRRLphiOJcfMgrBLoIpDg/rC8UweyGLV3EaOD0ZLrvcvld6ROK01SlBXi46Qr6QcxGRJ1mtGsWMgrsAwEvtEZIuIvCYip53cNxg5iL5wvORuz1IwDER1chBgJKr788yrjiRSk6LDZFFJN3VuktrSULLkNjL9B2VcIINeT1WqmKKJFNFEmqb6sb9Lq0/A6TBTbxlVXJPN7MbSuqlrOTZXkx87t7J/4PgqpggtQS9KGaEDJ8oHU2nFcKx6HgRAU31dgT6Iya0fbw/5ePlQf1nHjnoQxvp9Hjdej4sh00BYeaJSdJgsAl53VQbeD0SsJrnxHgTA9qODXLq8ssq6Qiilpq0HYXVT22mW0zmIyafgLaaI+EXkb4FPYPQ+HFFKvW49arbCGlJps9zdzx/iY/e8UvB96w64mgaiOVA3ropJKTXpMsltDeUrulp3jr6sCp0GnycTYsqUeJbjQfg8VQkxWeXQuZ3cswJe5s2qd7TUNRxPEUump50H0ZFplrP3vbBuFLQHMXkUi0H8ENgIvIbhRfxHTVY0iVheQ7m9EA9vPc4jW48XfN/SYWqseohp7HpjyTRKMamTuNpDPsLxVFldy9FkCm+O0GD2XOrecPlNYgGvuypJassoz6of/7tcOafR0RBTb42b5KqF1eB5csheHmIyJes1BsVuZVcppdYCiMj3gedrs6TJw/qD6ymzwWvfyWFG4qmCGkSDkep7EEaIaawHMZlS3xbtWZPlSp19EUukx80hDmZJflt3oOX2QVTDg8gYiDxrWDW3kSd2njBkIhww0pWE2CaTbD0mO0R1knrSKWaaM1cds6fhtKcSDyIcT3KkP1L0+EqUXAvRHDCE8bIVXTOu+SR6EG0V6DHlM7Ah00BEEynuePoAK2aHaCujDyDgdRNNpCsuROgvEGICo9Q1rXCsZDpjIKZdH4Qlt1GaB6FDTJNHMQNxtogMmo8hYJ31fKI5ENMV6460nBzE/u6RzPNCx2dmQdRXz4OwLlADWZVMU+EPq72C2dT5DIQl+f3tp/ZxpD/CP163uqy5BMEqjR21KsfyGQinK5mmqwdR53bRGvSWkIPQSerJpuCVSik1434r/jo3Qa+7rBDTvu7RwT19I/l7KYZi1fcgsgX7Ws2L8mj9+OSWuQJllbpGE+lxaw/6PBx+PcK3f72Pt5w9l4uWtZa1LkvRNRxPVfR76A8n8HpceY3w/OZ6Gv0exxLV09WDACNRbbebOjIFQqUzHZ39yWFec31ZapzZk916C4aYqp+DyCfYF5kCHagtQS8uKdODSObxIHweBiIJPC7h0286q+x1WR5EpZLf/eE4s+rr8noxIsLquU3OGYhwnDq3ZEQMpxMdIfvd1JHE+GIFTW3RBiKHZe0NY7wBu+zrHsk0TRUaOuSIgcgj2BebAiEmt0toCZZX6hrNM+zI+pndfOVy5jSNH7Npl4C3OpLf/eFE0WFFa+Y1suPYIIk80/4mYvvRQb766G6Uyp8n6R02ZDbKCbFNNrMb7XdTR/UsiElHG4gclrU38HpPmHiytD/svSeHWb9wFiKFx5YORo2whK+KLnM+wb6p0mBU7ujRaCKNLyfE9IYz2njT2k4+8IYlFa2pwVclDyIyXmYjmzXzmogn0+wtMDO8GP+96XW+/vgeHtmWv2S6XKHCqUAp3dSTLRej0QZiHMs6gqTSikO9IxPvbJJKKw6cGmHF7JChjVTEg2isovcAZNREs7upp0KICSoxEONDTJed2c5tf3JuxfLWAV91xo72hxN5eyAs1swzZkNsPVL6CNLXjhgd6F/+1e6888Z7R6avgcjupp6IyZaL0WgDMY4z2o3Zu3tP2jcQh3vDxFNplnU00BL0FsxBDEaqNwvCIuQzFF2zQ0xWknqy777aGrxlDQ1ycopY0PSqKu2FmCjEtKQ1SNDrLjkPEUum2HV8iLM6Q+w9OcxPXz4ybp++6WwgMs1ydg2E9iAmE20gcljaHgQoKQ9h7busvcEQ/CviQVQz/wCGomtuN7UV4moscodbCywPolAsvRBGDsKZr2bGg6iwm7o/Ulxu2+USVs1tLNmD2HV8iERKcfOVyzl7fhNff2xPptzTomc6G4gSJstFE840Gmrsow1EDkGfhzlN/jFVSRNhGYgz2htoDhaeSlfpPOpCzMrppt5zYpiOkG+c0mitaW/wEU+lMx3kdinUiV4N7HoQdz1/iLff9kze9zJKrkU8CDBGkG4/NlhSU96WLsOgrJvfxN9dexZH+iP8z+8OZd5PptKOiUnWgtkl6DFNtuCkRhuIvJzRUVol096Tw7Q1+GgK1BUdW2rkIKp/0W7KEezbe3KI5bOrM5SoEkrtnLXI1wdRLQJeezmIezYf5qVD/Znu92xGdZiKX6TXzGsiHE9x4JT9cOVrXQM0B+qY31zPJWe08YYz2rj1yb3EzCFKVjHCdDUQ1rxyO98JHWKafLSByINR6jpiOzSyr3uEZWZoqqXBmCmR71gnQkxgVDJZIaZ0WrHn5DDLO0JV/5xSWdZuGCnrrtgOSqm8fRDVwutxUeeWolVMvSNxXjlsJIrzhUKsn3XzBB7EmnlGR3UpYaYtRwZYO39WpoT1vRcuonckztYjRi6jrwKhwqmA12O/mzqiy1wnHW0g8rCsPchwLGnrS6yUYu/JYZZ1GBfDloCXREplhOWyMUJM1fcgskNMRwcihOOpKeFBrJrTSGvQy2/2dNs+Jp4ylWgdvDA0+Dx5hyxZ/HZPN5Z9PzYw3kBYnfIThZjOaG/A53HZNhDRRIrdJ4ZYZ1ZAAWxYOAuAlw/1AaNCktNNZiObjkY/3TY8CMOT1AZiMtEGIg/Wna+dMFPPSJyBSIIzzGMsjf5cuY1kKs1IPOVMDiLgzRiIPSeMNZ85e/I9CJdLuHR5G0/vOUXa9hQxo6zT51CSGozf794ThX+3T+3qps5t3MHnMxBW1/pEISaP28XKOY1sPWrPQFj5irXzRw1ER6OfebPqedn0aDIexDSU2bAwRo/azEF49SVqMnH0py8i14rILhHZKyKfzPO+T0R+bL7/OxFZbG5fLCIREXnFfHzHyXXmcobpDdhpcrKS2RkPwpyT3DMy9g9gdFiQAx5EoI7hWJJEKs2ek0MALO+YfA8C4NLl7fSMxG0L18Vq0MOxojPEzuODecOA6bTi17u7+f1VnQAcz2MgRkeeTvy7XDOvkW1HBm0ZyNeyEtTZbFjUzMuvmx7ENBXqy2Z2o89+DkLrME0qjhkIEXEDt2IMG1oFvEdEVuXs9gGgTyl1BvBV4ItZ7+1TSp1jPj7k1Drz0R7yEfJ5bHkQe60KJvOCPNrZPDZR7YTMhsWswKjcxu4Tw7SHfEW7fGvJpWe2AfDr3fbCTLnjRp3grM4Qg9FkXu9gy5EBekfiXL1qNm0N3vwhJptJaoA1c5sYiiVt6Xtt6RqgrcFHZ844zvULZnF0IMrxgWhF0/SmCotag5wYjBVtlstMRdRlrpOKkx7E+cBepdR+pVQcuBu4Pmef6zEm1wHcC/yeTAGBGRFhqc1Kpn0nR6ivczPH/KNuDRpVGr05IabBzDQ5Z0JMYIQ+jAT11PAeADpCflbOaeQ3dg1E0nkl2hWdRvJ41/Ghce89teskIkbndmeTn+MDkXH79EfieD0uW2sc7aie2IPaemSAdfObxmksrTfzEK8c7qN3JE7I76HOPX1DLxcuNZR4N+3vLbhPZiqizkFMzOBRSJRWKWgXJ79l84DDWa+7zG159zGHEg0Alo7zEhF5WUR+LSKX5vsAEfmgiGwWkc3d3fYToXY4o72BfTa6qQ+cGmZxWzCjOGmFHXKb5UY9CGeS1GDc2e49MTQl8g/ZXHZmGy++3pc3cZ9LLabhreg0fj478xqIbs6eP4uWoJfOxvr8OYhwguZAfiXXXJbPbqDOLbw2QaI6HE+y5+QQa+c1jXtv1dxGvG4XLx3qp3ckTus09h4Azp7fRIPPw7P7ThXcJ6bnUdvn538L3/s9R07tpIHI99eTG4gttM8xYKFSaj3wMeBHItI4bkelvquU2qiU2tje3l7xgrNZ1hHk+GA0bx18Nod6wyxqCWReN/g81LllnNxGZliQQzkIgB3HBhmJpzLhrqnCG5e3k0wrntvXM+G+tQgxNdXXMbfJP27iW+9InFe7+rl8hfFdmtPk53i+MtdwwlZ4CcDncXNWZyP3v9zFo9tPFCyd3n50kLQan3+wzrFmXiMvHzI8iOkcXgIjeX/+kpai34eq6oml0xAvXcJ/WpCMwcHfwsKLHDm9kwaiC1iQ9Xo+cLTQPiLiAZqAXqVUTCnVA6CUehHYB5zp4FrHYVUyZU+KyyWdVhzui7CwddRAiAjNAS+9w7kGwhoW5EwfBMDzBwyXfap5EOcubqa+zm0rzJTxIBwWaVvRGRoXYvrNbqO89fIVHQB0NvnpDycy2lYWfeH4hCWu2dxy/WpC/jr+4r8282d3PM+Wrv5x+1i9Ivk8CID1C5vZ0jXAyaHotPcgAC5e1sr+UyMcyxPCg2xF4ip8D576v/C1tdCzb+z24W545S7DgExXDj8PiTAsu9KR0zv5V/gCsFxEloiIF7gBeDBnnweBG83n7wSeUEopEWk34sppiQAAGiRJREFUk9yIyFJgObDfwbWOwzIQxSqZuodjxJNpFmR5EEBewT4nk9TWxeqFg4aBmEo5CDDugC9c2sJvbfRDRGukRLuis5F93cNj5jU8teskLUFvpg9hTpORV8r1IgYiiQmb5LJZv7CZX370Uj735lW8crif6771DH/4jd/y35te53+3Hufv7n2VbzyxhzlN/oyY3fhzzCKWTLP7xHBRDajpgjUR8Nm9+b2IaDVnmhz4DYRPwf+8E0bMz+t7He74fXjgQ7D30co/Y7LY9wS4PLD4DY6c3jEDYeYUPgI8AuwA7lFKbRORW0TkOnO37wOtIrIXI5RklcJeBmwRkVcxktcfUkoVzmg5wKLWAF6Pix1FyjMP9Rpu68I8BmJ8DqL640YtQj4PbpdwYjBGW4N3SoYgLjuznYM9YQ71FHf1o0krxOSsB7FyTohESmU8xHgyzRM7T3LFio5MPqnTNBC5d7mlhJgs6twu3v+GJTzzySv55+tXk1bw2Qe28qH/9yK/3HqcS5e3c+ufbCh4/PqFzZnn07kHwmJlZyPNgTqeLRBmqlqIKZWEY6/C4kth4Ajc/cdw9GW44xoI94CvCV77SWWfMZnsewLmnw/+cRH4quDozEKl1MPAwznbPpf1PAq8K89x9wH3Obm2iahzu1g1p7FoctG62OUaiOagd5xhGYwm8XlcFc8zyIeIMKu+jp6R+JSQ2MjHpcuNctdNB3rGhORyse4cqzlUKR+jiepBVnSG2LS/h8Fokj9Y05nZx5pcl9sL0ReOF5X6Lkajv44/vWgx771wEVuPDBKOJ9mwqHnCqqS5TX5zGltsWvdAWLhcwkXLWnlu3ymUUuMS/tF4lQxE9w5IRmDDjXDeB+AnN8F3L4eGTnjfL+H5/4Qt90B8BLzByj6r1oycMozfFZ9x7COmb61cDVg3v4mtRwYKNjkd6g0jAvNmjR2BmU/y2ymZDQsrzDQVJDbysaTNqOaZqHS4Fo1yAEvbGvC4JFPJ9L/bjhPwunmDaciATD9CdiVTNJEilkxX3GciIqyd38QFS1ttlayKCOsXGF7EdNVhyuWiZW0cHYjyeh6vMlKtENORl4x/522A1W+DP/iSccf9gUdg9mpY+05IjMCuX1b2OZPB/qcA5Vj+AbSBKMraeU2MxFPsL6DGebg3zNym+nFeQXPQS38kMUbmedCBaXJjPtO8YC2fYglqC7dLWNQaLJr0h+wqJme/ml6Pi2XtDew6PkQqrfjVthNcsaJjjGGq97qZFagb40FklFzL9CAqweqHOF0MxCVmHuKZPOWu1veg4ka5oy+BvwlalhqvL/hL+PNHoXmx8XrhxRCaC6/dO/a44e6pn7ze9yTUN8Pccxz7CG0gimBp4hQSWzvUG2Z+c/247a1BL0qNHQPqlJKrhdULMdUS1NksbQuyfwIPolZJaoCz5oTYeWyQlw/1cWo4xjVZ4SWLzkb/GA/C6pAvNm7UKa44q4O2Bt+Uq1IrlyVtQTob/XnzEHtODiFizBSpiCMvwdz1UKhnxeWCte+AvY9B2Exzbv0p/MeZsPXe/MdMBZQy8g9LLweXc38r2kAU4Yz2Bvx1roJy1Yd6w+PyD5Al2DfGQNQmxDSVLx5L2xs41BvOO2fZIppM4XZJTTqFV3SGODoQ5Z7Nh/G6XVyxYnwvjdELMZqkHvUgan8Xf+bsEJv/4apxVXPTFRHh4mWtbNrXMy6M+8TOk6xfMKuygotEBE5uh3nnFt9vzTshnYAdDxqG4qcfBJWGQ5vK/2yn6d4FQ0cdDS+BNhBF8bhdrJ7blBkin000keLkUCyvgbCSiNlyG0PRJI31znkQq+c2sXJO45QOPyxtD5JIKbr68te+gynx7KCSazZnmYnqn750hEvOaM1rwDub6jnWP+pBHOo1QmQdjRXe2WoAQ6urZyTO8wdHixRPDEbZ0jXA762cXdnJj78G6STMLVwdBsCcs6F1OTz7Lfjxn0L7WYbXcezVyj7fSfY+Zvy79ApHP0YbiAlYO6+JbUfHj43sMsXX8lXkWHIb2ZPlhqIJQj7nPIgPvGEJv/xoXkWSKYM1VGn/qcJhJifHjeZiaTIl04pr84SXwPAgekbimdDXb/ecoiPkY2nbNKt4maJcu3oOswJ13PnMwcy2J3eeBOD3VnZUdvLsBHUxRIxkdc8eCHXCn/4UFl0CJ7YaZbK1JJ2C524zPIRChHvhma/BvI0wa0Hh/aqANhATsNYcG5kbO7d6IPK5+9ZdfLaBGIw4m4OYDixtm7g7vZZDYuY2+Qn5PbgEripwt2r1QpwcjJFOK57d18MblrfZ0mHSTEy91817zl/Ir7Yfz9x0PbbjJPNm1bOi0nDp0ZeMctbGuRPvu/H9Rinsn94PDR2GV5GMwqndla2hFNIpeODD8Min4NHPF97vfz8JkT54y9ccX5I2EBNgaePk5iEK9UDAeMnvRCpNJJFyNAcxHWgOemkO1LGvmIFIpvA5XMFkISJsXNTMZWe201ogGTonq1lu+7FBekfimZ4OTXX40wsXISL893OvE02keHpvN1et7KjcCB95aeL8g0VDB1z3jdHqpjlnG/8ee6WyNdglnYKffQS23A2tZxghpEjf+P12PwJbfgyX/h/oXOv4srSBmICl7Q0EvO5xDXOHeiMEvO68ujj+OjdBrzvjQQw7KLMx3Vja3lC0kilW4yEx337vuXznvYUvItlyG7/dY5RjXnKGNhDVZO6seq5d3cldzx/i8R0niSbSXGk3/6CUcXHNJTpghIzmrS9vUa1nQF2gNnkIpeDnfwOv/ggu/xS843tm0vznY/eLDhjKrR2r4NKPO78utIGYELdLWD13fEe1VcFU6C6nOUtuw0kdpunGkrZgwb6SVFpxpD9KoIZDYvx17qIhrU6zm/rYQJSn93ZzVmeIjlB+vSRN+dx0yWIGo0n+8efbCHoN7a4JifQbMtcvfH/8e0fNO/+JEtSFcLmNO/RaGIgDv4GX/x+84WNw+SdhzjlG38bWHDGJRz8Hw8fh+m+BpzbFKNpA2GDtvFlsOzowpjzzcG+Y+c2Fyw1bgt7MeMhBB3WYphtL24N0D8Xyyqh/9dHd7Dg2yLs2zp+EleWnwech5PNw8NQILxzs4w3ae3CEjYuaWTOvke6hGJcub7cnteJvAo8ffvvl8XLeR140/p1bpgcBRpjp2BbnG+Ze+i/wz4I3/r3xWgTWvMMwHEMnjG2vPwcv3gkXfth+2KwKaANhg7XzG4km0pnxokqpgj0QFi1BbyYHcdhMaDvZST1dKJSo/tW243zryb28e+MC/mijs5UZpdLZ5Od/tx0nnkyPkeLQVA8R4aaLlwBwpd3qJRG48rMwfAKe/+7odqVg/5PQvAQCNjyRQsw525Dh6N038b7lEu41+i/WvRvqsjzTNe80ejG2/wyScXjob6FpAVzxaefWkgdtIGywdp4hcWANlT81HCeSSLGwZXwXtUVLwEvvSJzn9vXw8Z+8yoKWetbkGQYz08hX6rqve5iP3fMq6+Y38U/Xr55yFULWXAiv28UFS1onPkBTFm9bP4+v33AObz0nd/BkERZdBGdcZZR9Rs0w8LPfNO6+N76/sgXNMSUsjjqYqN5yD6TisOHPxm7vOAs6Vhvd3M9+A7p3wpu+XHNBQW0gbLC0LUiDz8NTu4x5BoeL9EBYNAe9HB+IctMPnmfOrHp+8pcXOzJNbrqxsDWAS0Y9iGQqzV//z0t4PS6+/d5zp+QMYitRfe6i5sq1gTQFcbuE68+ZV7ri8ZX/YFT8PHcb7H0cHvs8rHorXHxzZQtqXwFun3OVTErBSz808iSda8a/v+btcPh38OsvwarrYcW1zqyjCNpA2MDlEm66eDG/eO0Yj2w7ngkZTRRiSqYVy2c3cM9fXpSpp5/p+DxuFrQEMonqe1/sYufxIf71rWvGqeJOFSzZbx1emqLMXQ8r3wLP3Qr3vh/aV8Jbbyusv2QXd52h+OpUovrIi4YUyLk35n9/zTvMdXjh2i86s4YJ0EFxm/zN7y3niZ0n+fRPX+MtZxuNN8WS1Fevms2JwSgfv2aF9hxyMET7RojEU3z1sd1sWDirYCfzVMASZLxseXXnnmuqyBWfgR0PGYnrG/6neqGYOWcb4n1K2TM4g8dg50OGcTn3puL7vvRDqAuOGoJcWpbAJR81DGDjnJKXXg20gbCJ1+Piq+8+h7d882nufPYgsxt9RcMhZ84Occv1edxGDUvbG3hufw93PHOAE4MxvvmeDVMu75DNW86eS3vIl1H31UxBOlbCu+40LqotS6p33jlnw4s/gL6Dhc+bSsLL/wWv/Ai6XhjdHmiDlW/Of0xsCF67D9a8DXxFOsavvqXspVcDHWIqgRWdIT5+zZlA8fCSpjhL2oJEE2m+/vgerlrZwflLKqg0qQH+OjeXr6hQF0jjPKvfOtoBXS0yHdV5wkxKGZ3N374YHvr/IBE18iEfesZIcD94Mwwdz3/ep79mVEht/EB111tl/v/2zjxIiuqO45/vIgIKiiCHYjjkUFHBCnhxeCFyGCsoRIhXaaGIGI9YakSNwaPKGJOKhkJFyuBRhUcog2ASMaJGUyorGNRIvI2C9wUiKMvKL3+8NzpsZt1Zdqa7B36fqimZfj3jZ3u6+/fe+73XzwNEI5k4ZHeO7rdLvc/ucRpm9ziSqfabjVw8cs+UbRzne+jYF6qaw6sPbbr9m9qwfOmc48MTYyfMgclPwiEXhYTzcbPC48bnTfn/eRSfvRVGWvUb3/CDBFOmrAFC0khJr0h6XdIlBcpbSLo3li+W1D2vbGrc/oqkEeX0bAzNqsSME37ImYf2TFulYukVFzUaN2C3TK9f4Tg0bwkHnQXP373pqnMLp8LyeXD45TDlGdjz6E1zFB36wIhr4I1FUD1z0+98+HKo2gaOnJbEX9AkypaDkNQMmAEMB1YCz0qab2bL83abCHxuZr0kTQCuA8ZL6gtMAPYGdgUekdTHzAo8dMWpNDq2acmc0w+k/w/apq3iOA0z7ApYUQ0LzgtdTm88GibmDToHDr2o/s8NnAiv/R0emgqrV4bupxWLQxL7iF8W95TZlJGZNbzX5nyxdDAwzcxGxPdTAczs2rx9FsZ9npa0DfAB0AG4JH/f/P3q+/8NHDjQlixZUpa/xXGcrZzV78LMoWHU0Rcroc8oGH9Xw8t9rv8ytBiWzg4PAAT4ZgOcXb3pzOkUkbTUzAYWKitnF1MXYEXe+5VxW8F9zKwWWA20L/KzjuM4ybBjl5BXWL0COu0DY2cVtxZ0i9Zh3YaT50Htevj0dTjqmswEh4Yo5zDXQuMW6zZX6tunmM8iaRIwCaBr166N9XMcxymeXsPgjEXQrmfj51n0PBymPB0e/tdtUHn8ykA5WxArgfynru0GvFffPrGLaUfgsyI/i5ndamYDzWxghw4+iclxnDLTZQC02szcWYs20H1w02d4J0g5A8SzQG9JPSRtS0g6z6+zz3wgN898HPCohaTIfGBCHOXUA+gNVJfR1XEcx6lD2bqYzKxW0s+AhUAz4I9m9pKkq4AlZjYfuA24S9LrhJbDhPjZlyTdBywHaoGzfQST4zhOspRtFFPS+Cgmx3GcxpPWKCbHcRyngvEA4TiO4xTEA4TjOI5TEA8QjuM4TkG2mCS1pI+Bt5vwFTsDn5RIp1S4U3Fk0Qmy6eVOxZNFr3I4dTOzghPJtpgA0VQkLakvk58W7lQcWXSCbHq5U/Fk0StpJ+9ichzHcQriAcJxHMcpiAeI77g1bYECuFNxZNEJsunlTsWTRa9EnTwH4TiO4xTEWxCO4zhOQTxAOI7jOAXxAOFUFFL2HqYvya+jIpDUPG2HQmTxnMoKW8WJLalt/G9m/l5J3SV1StujLpJ6SeqZtkc+kvpLGg5gGUmaSRogaQqAmW1M2wdA0kGSrpfUOm2XfCQNlDQb2DNtlxyS9pZ0OGTqnOohqXPaHvlk5oZZDiS1lXQn8KCkrma2Me0gIWknSTOBB4A5kiam6ZNDUjtJM4D7gZtyN7+0ib/X34ALJe2Vty0tn/aSpgM3A2vjttRroHGdlZnAB0CNpCIWTC4v8ZyaDkwnLAyWhd9vO0mzgHuASyX9SlKLtHzqOD0A3CLpp5JaxbJUz60tOkAA5wMdgWeAsyHd2l6s2f0eqDGz/sBvgDMl7ZiWU/TaNnrVmlk/4CrgvAxcOM2A9sDjwNPACRB+wxQvnHuB4WZ2gJndEX1SrYHG5XrbA6PN7HdmVpP2AluxJnwbsMHMDgZOBYZC6i2uY4EdzGxf4BxgJJB2gB8DtI7X3tUEpxMkNUv73NriAoSkAyTlnityI3AGMAfoK2lw3CfRvzs6tTWzL4GbgHMBzGwh8CXxwkma6NXOzGqA883svFjUh3BDTrxLIDrtBBBvcl8DrYBXgA6SDo5liV040aljfHspsH3cPkLSJEmJ/3653y6+3RUYYmbvSholaZakU9LoKoxeOwAfAqeZ2QWxqD2wLu6TxvW3c3xbw3fr3fcntLiG5mrsKTltS7gPYGZLCfflscCgJJ0KsUUFCElHEFoLR0na1sw+N7MVhJvLY8AkSLYGk+d0dKwRLzEzk9RM0vaEZV+XJeVTwGuUpOZm9rmkKkkTgCuBl4FbJY1N6uLJcxodWzUAuwBvmdndwL+AsyRdngsiCToNl9TKzKqBxZI2Aj8HdgJmS5qQwnEaJamlmb0Tne4ETgOeAIYBF0jaLQmnOl7HAM3NbFVeYvoF4DhI7fobEa+/fwLLJD0M3EC4L1wGnCOpfQpOAj4idAtOipXbKmAl0CftLvEtKkAA3YDngJ7k1X7NbC2hH7ulpPHwbbM8aad96lwcRqghr07IpT6vvvDthbvAzLqa2a8JfccnpuS0V9z2CdA81raGAeOALrmAlpLTacBlZjbSzK4DpkWvpLp18p16x22PAwcCs8zsLkL3ZbO88qS99gQwsw2x7DHgzXhzTJKcUy+gr5m9T2jBvwsMNbM/ANcA+wLt6v2W8jj1jq9FwEPAkcC8+O9ngMNS7k6t7AAhqXOdm8R7wFygOzCozoF9G7gDGCfpCmBiOfrYG3A6OOcUu08GA2+b2RpJZ0saU66ToQivqui1Nm+/+cAOQMs0nOK23QgX9zJgDXAR0ExSl3LURBtwOkDSNma2xsyuzfutHgTaAtuV2qcIp1z31iPAEkKgwsxeInSlvFsOpyK86l5/HWN5WYPo9zh1A4bE36+WcFwOiPssIpxn6xJ26koICjVmtgCYCBwVA/xTwAZJVWnmISoyQEgaJulJYAahTz/HYOBVQu6hH6E74hBJMrN1hJN0NCEJ9ISZrU/LKZYdBHSX9BdCourFUp8MjfCaLGlovIA2ShoBLACqCTfmpJ36Szod6EyoHY8ws4nAw8BSQssraad+wKR4nBS7CnPHaTHpHKd9JU0mtG4uJRy3KyU9DXwMfFTqSsdmXn8fAJ2AUaV0aaTTvsDpkoYBq4CTJP2C0CX3KrC6lMeqSKd9gCkKeay1sYI2BrgdWJz6EGozq6gXIYG6mFBT6gj8lTCqBGA8oUa1O/Aa8AUhIkPoRnkWGJcBp9Gx7Abgv8DIjByr4YScyJmEm/CPU3ZaCxxe5/NVWThOsexkQl4k7eO0BhgVy7oSbsJjMnJOjcj77KG56zHlc2oIITH8I8LN+7gMHKfcfWoQ8B/g2HL8fo19JdUP3yTyuj82AvsB1WY2V2G0xFrCASWWXUIY9fIIYbRJTWymLQf2z4jT2lg23czOL5VTCbxqzKxW0v1mNjMjTparred9R9pONbFsgYXugJLQRKev4nF6B3inVE4l8Fqf+/3M7B8ZcGpNyM3UmtmDhO7BtJ3y71NP8V2uK3Uy38Uk6TRCRv/quOkFYIDCxJIXCc3W6yXdQBjLX01IPp0FvBnLSzrhpAROu8SL5o1SOZXQq8rMPs6QUydKPE69hMdpVYacSn6csurVRKc3CN2WWbsnlPw+VRLSbsI00ExrTcjqn0fI+u8Rt3cgJCvPiu9bAp8C+9f5fIutwSmrXu5UuU5Z9XKnZF+ZbkFYmFh2rpndSEhOTotFnxKG0f077vc1YTJcLgrnmnolS0Jn2SmrXu5UuU5Z9XKnhEk7QjUiSncmNMtyybiphD69PQijN5YBPbZ2p6x6uVPlOmXVy53K/6qoFeUknQmcZGZD4/vfEmbaVgEXW5g1vdU7ZdXLnSrXKate7lReKiZAxKTgRklzCVPT1wH3EeYOfOVO2fZyp8p1yqqXO5WfTOcg8okHfTvCmOLjgXfMrDrNg55Fp6x6uVPlOmXVy53KT0XMg8hjCmGUwHDLTmIni06QTS93Ko4sOkE2vdypjFRMFxN813xL2yOfLDpBNr3cqTiy6ATZ9HKn8lJRAcJxHMdJjorJQTiO4zjJ4gHCcRzHKYgHCMdxHKcgHiAcp0RImibpwu8pHyOpb5JOjtMUPEA4TnKMIS7v6jiVgI9icpwmIOky4BRgBWEFt6WENcYnERaleZ2wsNB+hLUHVsfX2PgVMwhP/VwHnGFmLyfp7zjfhwcIx9lMJA0gLA15IGHS6XPALcBsM/s07nMN8KGZTZd0O/Cgmc2NZYuAyWb2mqQDgWvN7Ijk/xLHKUylzaR2nCwxFPizhfXOkTQ/bt8nBoa2hLUCFtb9oKTWhOUl/5S3RkyLshs7TiPwAOE4TaNQE/x2wprQz0s6FTiswD5VwCoz2698ao7TNDxJ7TibzxPAsZJaSWoDHBO3twHel9QcODFv/zWxDDP7AnhL0k8gLDUpqX9y6o7TMJ6DcJwmkJekfpuwJvFywgL1F8dtLwJtzOxUSYOBWcB6YBywEbiZsFZAc+AeM7sq8T/CcerBA4TjOI5TEO9ichzHcQriAcJxHMcpiAcIx3EcpyAeIBzHcZyCeIBwHMdxCuIBwnEcxymIBwjHcRynIB4gHMdxnIL8DzsH2Wx00smdAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Date when handwashing was made mandatory\n", "handwashing_start = pd.to_datetime('1847-06-01')\n", "\n", "# Split monthly into before and after handwashing_start\n", "before_washing = monthly[monthly['date'] < handwashing_start]\n", "after_washing = monthly[monthly['date'] >= handwashing_start]\n", "\n", "# Plot monthly proportion of deaths before and after handwashing\n", "ax = before_washing.plot(x='date', y='proportion_deaths', label='before')\n", "after_washing.plot(x='date', y='proportion_deaths', label='after', ax=ax)\n", "ax.set_ylabel('Proportion deaths')" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "586a9f9803" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 7. More handwashing, fewer deaths?\n", "

Again, the graph shows that handwashing had a huge effect. How much did it reduce the monthly proportion of deaths on average?

" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "dc": { "key": "586a9f9803" }, "tags": [ "sample_code" ] }, "outputs": [ { "data": { "text/plain": [ "-0.0839566075118334" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Difference in mean monthly proportion of deaths due to handwashing\n", "before_proportion = before_washing.proportion_deaths\n", "after_proportion = after_washing.proportion_deaths\n", "mean_diff = after_proportion.mean() - before_proportion.mean()\n", "mean_diff" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "d8ff65292a" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 8. A Bootstrap analysis of Semmelweis handwashing data\n", "

It reduced the proportion of deaths by around 8 percentage points! From 10% on average to just 2% (which is still a high number by modern standards).

\n", "

To get a feeling for the uncertainty around how much handwashing reduces mortalities we could look at a confidence interval (here calculated using the bootstrap method).

" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "dc": { "key": "d8ff65292a" }, "tags": [ "sample_code" ] }, "outputs": [ { "data": { "text/plain": [ "0.025 -0.101488\n", "0.975 -0.067532\n", "dtype: float64" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# A bootstrap analysis of the reduction of deaths due to handwashing\n", "boot_mean_diff = []\n", "for i in range(3000):\n", " boot_before = before_proportion.sample(frac=1, replace=True)\n", " boot_after = after_proportion.sample(frac=1, replace=True)\n", " boot_mean_diff.append(boot_after.mean() - boot_before.mean())\n", "\n", "# Calculating a 95% confidence interval from boot_mean_diff \n", "confidence_interval = pd.Series(boot_mean_diff).quantile([0.025, 0.975])\n", "confidence_interval\n" ] }, { "cell_type": "markdown", "metadata": { "dc": { "key": "0645423069" }, "deletable": false, "editable": false, "run_control": { "frozen": true }, "tags": [ "context" ] }, "source": [ "## 9. The fate of Dr. Semmelweis\n", "

So handwashing reduced the proportion of deaths by between 6.7 and 10 percentage points, according to a 95% confidence interval. All in all, it would seem that Semmelweis had solid evidence that handwashing was a simple but highly effective procedure that could save many lives.

\n", "

The tragedy is that, despite the evidence, Semmelweis' theory — that childbed fever was caused by some \"substance\" (what we today know as bacteria) from autopsy room corpses — was ridiculed by contemporary scientists. The medical community largely rejected his discovery and in 1849 he was forced to leave the Vienna General Hospital for good.

\n", "

One reason for this was that statistics and statistical arguments were uncommon in medical science in the 1800s. Semmelweis only published his data as long tables of raw data, but he didn't show any graphs nor confidence intervals. If he would have had access to the analysis we've just put together he might have been more successful in getting the Viennese doctors to wash their hands.

" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "dc": { "key": "0645423069" }, "tags": [ "sample_code" ] }, "outputs": [], "source": [ "# The data Semmelweis collected points to that:\n", "doctors_should_wash_their_hands = True" ] } ], "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.7.6" } }, "nbformat": 4, "nbformat_minor": 2 }