{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "from sklearn.linear_model import LogisticRegression"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Data are available in the MOOC repository. The CSV file contains data for the 1314 women that were polled in Whickham, England, in 1972-1974 and were categorized as \"currently smoking\" or \"never smoked\". Each line is related to a person and contains whether she smokes or not, whether alive or dead twenty year after the survey and her age at the time of the survey."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of rows: 1314\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Smoker</th>\n",
       "      <th>Status</th>\n",
       "      <th>Age</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>21.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>19.3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>No</td>\n",
       "      <td>Dead</td>\n",
       "      <td>57.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>No</td>\n",
       "      <td>Alive</td>\n",
       "      <td>47.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>81.4</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  Smoker Status   Age\n",
       "0    Yes  Alive  21.0\n",
       "1    Yes  Alive  19.3\n",
       "2     No   Dead  57.5\n",
       "3     No  Alive  47.1\n",
       "4    Yes  Alive  81.4"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data_url = \"https://gitlab.inria.fr/learninglab/mooc-rr/mooc-rr-ressources/-/raw/master/module3/Practical_session/Subject6_smoking.csv\"\n",
    "data = pd.read_csv(data_url)\n",
    "print(\"Number of rows:\", len(data))\n",
    "data.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The CSV file does not contain any missing data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of rows with missing values: 0\n"
     ]
    }
   ],
   "source": [
    "print(\"Number of rows with missing values:\", data.isnull().any(axis=1).sum())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "\n",
    "Let's visualize the number of women alive and dead after twenty years, according to their smoking habits. A heatmap is effective in this case."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAi8AAAHwCAYAAABuVI8jAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3XWYFdXjx/H32V526ZLu7hKQVMAAUSxsvyoGKip2d3ciovITCxFEDMRACaWREESQDsmltoPdnd8fZ3a5d3fvBizuDnxez8PD3Zk5M2fmTnzmTFzjOA4iIiIiXhFU0hUQERERKQqFFxEREfEUhRcRERHxFIUXERER8RSFFxEREfEUhRcRERHxlFIXXowx44wxT5d0PQrDGHO1MWZOSdfjeGWMuckYs9sYk2CMqVzEstnrkTGmlzHmn2NTS79p1jfGOMaYkGMwbs9sF8XNGHOeMWabux50KOn6yLFhjJlljLnO/Xy5Mebnkq5TSclvezfGvGuMeeQYTPOo9l9u2cYB+hX791lgeDHGbHYPIFE+3a4zxswqzopIySitB0VjTCjwKnC64zjRjuPsO9JxOY7zu+M4zYqxbu8ZY24orvFJgV4GRrjrwTJ3n9S/uEZujAk3xow1xmwxxsQbY5YZY87KMUw/Y8waY0ySMWamMaZecU1fcnMc5zPHcU4/FuMu7vXnv+Y4znDHcZ4q6XoURc7vM7+gU1iFbXkJAW4/mgmVBGNMcEnXoTgcizN5D6gORACrSroieTgTmFbSlTiB1KOY1gNj5dzvhQDbgD5AeeARYKIxpr5bpgrwldu9EvAH8EVx1OdIBZgPTzhB92dS3BzHyfcfsBm4H9gPVHC7XQfMcj/XBxwgxKfMLOA69/PVwFzgNeAgsBE4xe2+DdgD/M+n7DjgXWA6EA/MBur59G/u9tsP/AMMzVF2NPbAkgj0z2N+rgFWu+PeCNzo068v8C9wl1uvncA1Pv0rA98CccAi4ClgToDllrVcbgB2uOO6y6d/kLtcNwD7gIlApRxlhwFbgd/c7j2Bee5y3AZc7XYPx56dbgV2u8svsqB5cut2CEgDEoDv3O5Z9YoH/gbO86l3MPAKsBfYBIzw/f6xO/+x7nS2A08DwQGWUTjwurt8drifw4Gm7vfnuPWaEaD8JGAXEAv8BrTKsS487bsMfObtyxzjeQN4szD1B9oCK3yWxcvustgI3FLYZQE0Ama43/1e4DPc7cvt3wFY6n4HXwATsuYnj+WwBejkfr7CrUNLn2316/yWd4715F6f9WQIMBBYi93eHizi+vs/7Dq5F3gon33MIGAZdrvaBjzuU98Ed1yJ7rQ+ATKBZLffve6w3Ti8bfwJ9M2xP3oGux9KBhoXYr+3ArjAZzuZ59Mvyh1P8zzKXQQsydHtrhzfQaBttSIwFYgBDrifa+c3H9j96EZ3PdkEXF7QvBW0H3T7nwssd7+TDcCZbvdKwIfu+nMga77cftcD69115Vugpk8/B7t9rAM2ud0GAGuw2+/b2H2973FjTo7yw93yB4BRgCnMPinHfAVaf87BBuSD7nJuEWC5GeyxbI9b7xVAa599zjvAD+645wInYbezA+68dvAZVwt3WgfdaZ8TYP9VFpgJvOlO37dfXwo+Zn3nfo+Lsfuggo5ZeW63wMnAfLe+O93vLCzHd3Qbdn3aC7wEBOX8PrH76qxtOgG4GKiCXd8PYtef37PKBlyHC7GSbwb6Y888shZYUcNLOnZjCXYX3lbsyhcOnI7dgKJ9vrR4oLfb/w2fmY7C7tyuwZ4tdXQXUiufsrFAD+zONSLAjrKRuxL0AZKAjj4rQjrwJBCK3XEnARXd/hOwO+kooDX2gFTQivC5O3wb7E6pv9t/JLAAqO3O5xjg8xxlP3bLRgJ13eVyqVu3ykB7d/jXsTuLStgV/TvguULO0zhyHBSxO+Ca7jK8GLuS1XD7DccGmtrYne0v+B+wv3bnJQqohg15NwZYRk+6y6AaUBV78Hkq0HqVR/lr3fnNOigvD7Dx9+VweKnnzn85nx3fTqBbYeqPPWA/57Ms1gB13GU/s7DLAnvgGeDWvSp2g37d7ReGDSR3uN/ZhdiQGSi8fIwbjIH3sAebm3z63VGI5Z21njzqTvN67Po63l3GrYAUoGER1t/3setuOyCVwAeEvtjtIwgbDncDQ3LsFBv7/L0ZnxMToBY2QA10xzHA/buqz/5oqzsPIUBoAfu86u68Nnf/fgMYnWOYv3DDTY7u4didbwufbss4HITy21YrAxcAZdx+k/APBznnozz2oNTM7V8DnwBfwDzmtx88GbsfHeAuz1o+y+J7bJiu6K4nfdzup2H3xR3dZfAW7kmXz3c43Z3vSOzBKg67bodi1/V08g8vU4EK2H1hDIcDVb77pEDHNJ+/s06WBrh1uRcbwsLyKHsGsMSth8EGkKx94zh3GXTCthrPwIapqzh87JvpDhvqTuNB7PZ+Gnb/3sxnXE+768QifLZ9cu/bCjpmTcCuUy2xx8+Cjll5brfufHXDrnv1seF3ZI7vaKb7HdfFnvTk9336btPPYYN8qPuvF244DbgOF2Il34wNL62xK3RVih5e1vn0a+MOX92n2z4OH4jHARN8+kUDGdgDxMXA7znqNwZ4zKfsx4XZeH3Kfw3c7rMiJOeYlz3uFxaMPYA09+n3bCFWBN/hXwTGup9XA/18+tVwxx/iU7ahT/8HgCl5TMdgN7xGPt26c/jsJuA85dwQ8llGy4Fz3c8z8D+Y98/6/rE7/VTcM0m3/6W4G2we490ADMyxY9gcaL0qoI4V3OHLB9jA//UZdg5wlft5ALDB/Vxg/bFnBL18lsVwn36nH8WyGAIscz/3xp7ZGp/+8wJ9T9gWum991qvrcLchbAjKOijlt7yz1pOslqGy7rx09Rl+CW6ooHDrr2+rwSLgkkJ+l68Dr/n8XVB4uQ/4JMc4fsJt0cXuj54s5LRDsQe/MT7dxgLP5xhuLm7LZx7jGA08435uhT3rDqeAbTWP8bQHDvj87Tcf2FB8EBt4Igszf/nMt+9+cIzv8s/xHWfiHhhz9BsLvOjzd7S7PtT3+Q5P8+l/FbDA52+DbUHI72DX0+fvicD9PtthnvukAPOac/15BJjo83cQ9sS0bx5lT8MelLuRo2UAu8953+fvW4HVPn+3AQ66n3thW42DfPp/zuFWx3HA/2FD8j15TMd331bQMauZT7/CtLwUarvFnsBM8fnbwQ2U7t83A7/m8336btNPAt9QiFbRrH+FvmbqOM5f2OR7f2HL+Njt8znZHV/ObtE+f2/zmW4C9kymJvasuasx5mDWP+BybNNcrrJ5McacZYxZYIzZ75YfiD0LyLLPcZx0n7+T3LpV5fC18Sxb8ptWHvXZ4s4H7rxM8ZmP1diQVj1A2TrYg09OVbGpeonPuH50uxc0T3kyxlxljFnuM77WHF5GNXPUy/dzPezOf6dP2THYM/281MR/Gfoun3wZY4KNMc8bYzYYY+KwOyTw/y4DGY8NEgCXuX8XWH9jTAXsZct5PvUPtD4UNK5qxpgJxpjtbv0/xX8Zb3fcrTqPcec0G+hljDkJu8P6Aujh3rNRHhs+s8ab3/Le5zhOhvs52f0/0HZamPV3l8/ngOucMaarexNsjDEmFnsmXZjvMUs94KIc+4We2INtlnz3C249grCXFdKwlx6yJADlcgxeDnumnJePgMuMMQa4EntgTKWAbdUYU8YYM8a9cTgO2xpXIce9e777xkTsCd1w7Hr2vTGmeUHz6U4rv/1goH1NHWC/4zgH8ujnt265++192FabXHUnx7bjrusFfUeB1qf89kmFkbPume44auUc0HGcGdjLJaOA3e7N+77rRs7tJdD2UxPY5k4ry5Yc0xyEbQF5t4D6F+WYVZhlk+dyNsY0NcZMNcbsctfPZ8m9nQY63hXkJWxL1M/GmI3GmAJzRlFv+HoM25zsu4AT3f/L+HTzDRNHok7WB2NMNLYZagd2wcx2HKeCz79ox3Fu8inrEIAxJhyYjL3mXN1xnArY+2NMIeoUg22eq+PTrW5R5sUdfof7eRtwVo55iXAcZ7vP8L7zsg3bzJvTXuxG0cpnPOUdxwkYTnLwW17uUxTvY3feld1l9BeHl9FObPNsXvO3DdvaUMWnLuUcx2kVYNo7sAeeLL7LpyCXYa/L98ceoOtnzUIhyk4C+hpjagPncTi8FFT/M7BnElkH+J0EXh8KGtdz2GXf1nGccth7VXyXcS334JfXuP04jrMeu5O5DdtUH4/dAd2APdvJ2kEezfLOqTDrb2GNx15KqeM4Tnnszjq/7zHnNr4N2/LiW5cox3Gez6eMH3dZj8WGrwscxznk03sVtgk9a9go7LaY503EjuMswAagXtj19BO3V0Hb6l1AM2xrVzlsCxz4Lwu/+XAc5yfHcQZgg9oa7Labr0LsBwPta7YBldwQn5PfuuUuo8rYFoy86u637bjL33dbKor89kl5ybku5Kx7Vl3yXJcdx3nTcZxO2Fa1psA9Ra2wO806OW66rptjmu9jw+0036d9iyDrmFWUZZOf0dh1rIm7fj5I7u000PEuX47jxDuOc5fjOA2BwcCdxph++ZUpUnhxd5JfYHeSWd1isAv8Cvds+FryXvGLYqAxpqcxJgx7U+xCx3G2YVt+mhpjrjTGhLr/uhhjWhRyvGHY5tsYIN19HLJQj+O5B6yvgMfdM6SW2BubCvKIO3wr7L06WU8pvAs8k/XIpTGmqjHm3HzG8xnQ3xgz1BgTYoypbIxp7x6Y3gdeM8ZkndXXMsacUZj5wp4ZNPT5Owq7cce447oG2/KSZSJwuzuNCtgmewAcx9kJ/Ay8YowpZ4wJMsY0Msb0CTDtz4GH3Xmvgr3f4tNC1rssNhzswwbnZwtZLmudnYW98XCT4zirC1n/Qfg/ZTQRuM0YU9sYUxGfVslCjKss9oz+oDGmFv47wPnYnc5t7nd9PvY+hPzMxgbO2e7fs3L8DUe3vHMq6vqbn7LYM/oUY8zJ2AN+fnKus58Cg40xZ7j7oAhjTFY4LazR2PsXBjuOk5yj3xSgtTHmAmNMBHa5rXAcZ00+4/sYe4ae7jjOHMg+o89vWy2LDTcHjTGVsCeLARljqhtjznEPbKnY9SnD7Zf1zo76eRQtaD84FrjG2MfDg9w6NnfX6R+Ad4wxFd39b1bAGu+Wae+Go2ex++3NAar/PdDKGHO+sU8f3caRn/QG3CcFkHP9mQgMcuc3FBsiUzncwprNPd50dYdLxN4blZFzuEJY6Ja/112OfbEH7Qk5hhuBfTBlqjEmsigTyOOY1Rx7ue5IlcXep5TgjuumPIa5x1036mCfUA70VJ7fd2CMOdsY09gNjnHYZZrvcj2SR+2exB7gfF2P3fnuw6bRXF96EY3Hbrj7sTcJXQ42nWE3skuwiW4X8AJ2QyyQW/427Mp6ALuT/LYI9RqBbULbhb3u+GEhyszGNof9CrzsOE7Wi3recKf9szEmHnvzY9d86r4V27R7F3a5LOfw2eB97jQWuM15v2DP4ApjLNDS2Gbsrx3H+Rt75/587ArWBnt9P8v72IPyCuyNiNOwB9qsFe0q7M7xb+wy/hL/5ntfT2MfO10BrMQ+XVPYd858jG2W3O5Oa0Ehy2UZj221GZ+je571dzeqAdgzoSzvY++t+NOt+1eFGZfb7wnszY2x2B15dlnHcdKA87HXiQ9gLw3kHHdOs7E7l98C/A1Ht7xzKtL6W4CbgSfd8TyK3T7z8xw2hB00xtztnticiz0TjMG2ENxDIfdvbgC7EXuPyS5jX4aXYIzJ2u/EYO8reQb7fXTF7oPy8wk29H+So3t+2+rr2MsEe7HL80fyF4TdH+zA7hP6YJcl2DPgrO3DT0H7QcdxFmFPtF7Drp+zOdwycSX2Poo12HsrRrplfsXeOzIZ2xLSiHyWkeM4e7EPBjyPPW40wX8/UxQF7ZNyyrn+/INt+XwLu+wHY0NsWh5ly7nTO4BdvvuwLVhF4o77HOAsd5rvYO/DW5NjOAfbgroN+MYNz0UxAtsyvQu7Ln6ODWZH4m7suhKPXQZ5BZNvsPfGLcfu18YGGNfjwEfudzAU+/3/gg3g84F3HMeZlV9lsh41k2LmnvFswj7ZkJ7/0N7lnrW96zhOvQIH9jC3ReBtx3EKagERwT1L3oO9WXpdCUz/YSDGcZwx//W0S9qJsk86EsaYF4CTHMcpzFWDUk0vC5IicXfKp2LPdKpjW8imlGil/jv5NuOL+LgJWFwSwQXAcZxS99bsY+UE3yfly728E4Ztae2CfTLxuhKtVDFRy8sxcry2vBhjymCbkZtjr89/j33EMq5EKyZSShhjNmNvZBziOM6yEq7OcU/7pMCMMV2wl4pqYlsCx2Af+/f8gV/hRURERDzFk7+NISIiIicuhRcRERHxlBP2ht3IDiN0vUzkGJox6YS5Z1SkxHRvXKEwL+Y87qjlRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8JaSkKyDHpysGd+X9J68scLiBw99i5sJ/8uxXrVJZ7ry6P2f1ak2dkyqSnHqI1Rt28unUhYybMj/gOBvWqcLZfdrSu3MT2jStRfXKZUnPyGTHnljmLlvPexN/Z9nqbUc8byKlQUJcLMsW/sbfy/9gy4Y17N2zi8yMDMqWr0D9Ji3o2W8QnU7pm2fZNSuX8tfShWxev5qYXduJj40lNSWJMtHlqFW3AZ1O6UufM84lLDwi3zrs2fkv0778hL+WLeTg/n1ElomiXqOm9DlzCF16nHYM5lrEMo7jlHQdSkRkhxEn5oz/R7LCS0ZGJjEHEgIPd+9Y5i7bkKt7hxZ1+HbULVSpGA1AfGIKEWGhhIYGAzB93mouHDmGtEPpfuW6t2vIjHF3+nWLS0gmPCyE8LBQADIyMnlh7E88Nfr7o5pHyd+MSU+XdBWOa8POOYWMjIzsv0PDwgkKCiI1JTm7W9vO3bnlgecJj/APIa89fid/Lp6b/Xd4RCSAX9mq1Wty11NvcFKtunlO/8/Fcxn13IOkpaYAEFkmipSUZJzMTAB6DTiba29/GGPMUc6p5Kd74won5AJWy4scU//uPkDzQY8VqUy56AgmvzGcKhWjWbNxF8Me+Zilf28lNCSYa8/vwYt3n8+AU1rw4t3nM/K5iX5lQ0KCSU/PYNpvfzHhh8XMXryO/bGJBAUZOjSvwwt3nU+Pjo158Iaz2LpzPx99HbgFR6Q0y8jIoGHTVvTsP4jWHbtRrUYtAGJ27+C7CR/y28/fsuKP+Yx7+zluvPsJv7It23ehdcduNG3Vjmo1ahNZJgqwrTnzZ/3EpHFvE7N7B28+fS9PjxpPUJD/HQYxu3bwzvMPkZaaQpOWbRk28hFOqlWXlOQkfpj8Kd98Ppbfp0+lRu36DLyw4BZYkaLSPS9S6oy8qh81qpYnKTmNIbeOZunfWwE4lJ7BmIm/8dS70wAYdn4PGtet5ld247YY2l/wNBff9T5TflnO/thEADIzHZb8vZWzbnyLFWv/BeCea07/D+dKpHjd9+woHn3t/zht0AXZwQVsi8m1tz9E37POA2D+zB/ZF7Pbr+wZQy5lwDlDqdeoWXZwAYguV54B5wzlsuvvAGDH1k1sWLMy17S/+nQMqSnJlK9YmZGPvZrdOhMRWYbzrriBvmcOAeC7Lz4kMT6ueGdcBIUXKYUuP7srAJN+WsKWHfty9R/9+SziE1MICQnmkoGd/fpt33OQDVtjAo77UHoGE75fDECjulWpUDayGGsu8t9p0a5zvv17n35O9ufN61YXadyNmrfO/rx/7x6/fqkpySyZNxOA0waeT1R02VzlBw39HwDJSYksXTC7SNMWKQyFFylVmtSrRt0alQD4ee7feQ6TmJyWfZ9M/+4tijyNlLTD98kEB2sTkONTaFhY9ufMzIx8hsztn1XLsz9Xq1Hbr9/aVX+SlpoKQJvOp+RZvmr1mtSsUx+Av5YuLNK0RQpD97zIMVWlYjRzP7uXpvWrExxk2LU3jgV/buTDKfP5fcm6XMO3alwz+/OqDTsCjvfv9Ts4s2crmjc4qch16t25CQA7Y2LZdzCxyOVFvGDNiqXZn2vXb1zg8GmpKezfu4fFc2bwzedjAWjWugMNmvifIPy75fAN9rXqNgw4vlr1GrFj22a2b91U1KqLFEjhRY6pqMhwOrasy/7YRKIiw2hQuwoNalfh0kEn89HX87nl6c/JyMjMHr5G1fLZn3fsiQ043qx+5ctGEhUZRmJyWqHq07VtAwb3bQvAh1PmHcksiZR6iQnxfD/pIwCatmpPjdr18hzu4P59jLxyYJ792p/ci+vufCSPMnsBiIoul+spJl8VK1d1hw98GVfkSCm8yDGxMyaWp9+dxjczlrN28x7SDqUTFGQ4uU19Hh4+iH7dmvO/Id1JSknjzhcmZZeLLhOe/TkpJXAg8e1XNiqiUOGlSsVoPnruaoKDg1i3ZQ+vjpt+hHMnUnplZmby3iuPc3D/XkJCw7hi+F0Bhw0KDqJcBXuZNjkpkUNp9nJQl579OO+KG4guWz5XmZQk21oZFh6eq5+vrHfEpCQlHdF8iORH4UWOiV8XrOHXBWv8umVmOiz4cxODbx7FF69cx+BT23HDRb0Y9fmsfG+yLQ5RkWF8+fqN1KtZmbiEZC6/d2yhW2tEvOSzMa/y56I5AFx18z3Ubdg04LDlylfkzc9+AMBxHA7s28PMaVP4ccp4li6YzZXD785+akmkNNHdivKfcxyH+1+bAtgbZgf1bpPdLyEpNftzmYiwXGXz6hefmJLv9MpEhDHlrZvo2rYB8YkpnHfraFau3X6k1RcptSZ88Aa/TrUtmZdeP9LviaOCGGOoVKU6F1w1nBvvfoKM9HQ+eudFtm5c6zdchPtoddZNu4FkvbwuokyZosyCSKEovEiJ2LhtLzEH4gFoUKtydvedMYfvc6lZLXeTdc5+sfHJ+bagZAWXXp2akJCUynm3jWbe8o1HW32RUueL/3uLH6eMB+Dia2/ljCGXHvG4Ovc4lcrVauBkZvLbz9/59atQqQoAiQlxpKYEPnE4sC/GHb7qEddDJBCFFylVVq0//IRRq0Y1Aw7X0n0qac2mXQGHyQouvTs3ITHZBpe5S3P/FIGI100Y+yY/TP4UgKHXjuCsC6446nFWdEPKnp3+vwNWu16j7M/btwY+EdjuPpVUq26Do66LSE4KL1IiGtSuQtWK9uVWm31eRLduyx627twPwIAeeb/DpUxEGD062B3oL/PzfvlWmYgwvn7bBpeEpFSG3DqaOUvWF+csiJQKEz54gx+/+gywwWXgBUf/On7HcYjZbU8kIiKj/Po1bdUu+2bdlUvy/nmNvXt2smPbZgBad+x61PURyUnhRUrEcyPt68MzMjKZ9ttffv0+m2pfanXRGZ2yX1jna/jFvSkbFUF6egYTpv2Rq39WcMm6VDTk1ncUXOS4NOGDNw5fKhp2W6GCS0ZGeoHD/D79O2IP2JOK5m06+vULj4ik0ymnAjBz2lckJeb+4dVpkz4B7M8FdOzWp8DpiRSVwosUu7o1KvH7J3cz7IIe1Pe5n8UY+6j0N2/fzLn92gPwweQ5rNvi//rx1z/+lZ0xsURFhjPlrZvo0KIOAKEhwVx/UU8evXkQAGO/msv6rf5lIyNC+erN4fTq1IT4xBSGjHhHl4rkuDTx/97ODi6XXjeSs86/vFDl1q76k2fvvZG5M6axf6//bx7t2r6ViR+O4qO3nwfs23V79h+UaxznX3Ej4RGRHNy/l9efuItd2+3vj6WmJPPN+A+Y+cNXAJxzybVElS13xPMoEohxHKek61AiIjuMODFn/D9Qt0Yl/pn2ZPbfKamHiE9KoWyZCCLCQ7O75/WSuiwdWtTh21G3UKViNABxCclEhIcSFmqf7p8+bzUXjhxD2iH/s8jLzj6ZsU9dBUByShqxCfk/iXTp3e+z4E+9AfRYmDHp6ZKuwnFr355d3HXNuQCYoCDKlquQ7/BnnX959n0wq1cs4YUHbs7uFxoWTkREJKmpyX5PENVp0ITbHnmRqtXzvvfsz8VzGfXcg9lPFUVGRZOanJz9UwQ9+w9i2MhHMMYc+YxKgbo3rnBCLmC950WK3Z798dzx/ES6tm1A22a1qVIxmoply5CSdojN2/exYMVGPv56AfP/DHyz37LV2+h04TPcdc0AzurVmtrVK5CYnMbilZv5dOpCPvp6AXkF7yCfHWVkRBiR+TxuDRAaok1AvCfTORz4ncxM4g7uz3f4lJTk7M/1Gzfn+rseY82KpWxev4bYA/tIjI8lJDSMajVqU69RMzr3OJUuPU4jKDg44DjbdenB06M+4/tJH7Nq+SIO7ttLmeho6jVqRt+zzqMvksy/AAAgAElEQVRLj9OOfkZFAlDLi4gcE2p5ETn2TtSWF93zIiIiIp6i8CIiIiKeovAiIiIinqLwIiIiIp7i6fBijKltjJlijIkxxuw2xkw2xtQu6XqJiIjIsePp8AJ8CHwL1ABqAd+53UREROQ45fXwUtVxnA8dx0l3/40DAv6EqTHmBmPMH8aYP9L3rvrvaikiIiLFxuvhZa8x5gpjTLD77wpgX6CBHcd5z3Gczo7jdA6p0uo/rKaIiIgUF6+Hl2uBocAuYCdwodtNREREjlOefje64zhbgXNKuh5e1L55bQb2bkOHlnVpUrcqVSqWpVxUBHGJKazdvIsf5/zN+5N+50BcUsBxnNe/PZef3ZWOLetSuUIUh9Iz2L77IHOWrmfMF7+xYu32I6pbwzpVOLtPW3p3bkKbprWoXrks6RmZ7NgTy9xl63lv4u8sW72tUPN42xWn0btzE6pUjGZ/bBKLV27mnQmzmb14bcByHVvW5ZGbBtG9XUNCQ4L5e8MOXvvoV776ZVnAMr07N+Gn92/n+9kruXDkmCOabzn+pKak8M9f9jX8W9b/w+b1a9gXswuAcy+7jvMuvz7PcjG7d3DPtecVejo9+w/iujseLVLd3n/1Seb++n2Bw439di7BwYEPFZvXr+GnKZ+zZuUS4mMPElW2HI2at6b/4KG0bNc5YLlN61Yz5dMxrPt7BRkZ6dSq14izzr+Ck3v1C1hm9Z9/8MKDt9Du5J7c8dgrBdZdjl+eDC/GmPy2UsdxnKf+s8p41P/O7c7wSw7/VH1yShrJqWlUrhBF9/aN6N6+ESMu78tFI99j4Qr/Hy4MCw3hs5eGcXafNtnd4hNTCAsNpmn96jStX53/ndudB16bwlufzSxSvbq3a8iMcXf6dYtLSCY8LIQm9arRpF41rhzcjRfG/sRTowPveK8+rztvPnAJoaH2t1kOxidRvXJZzjmtHeec1o6n353GM2Om5SrXpXU9fnr/diIjwkhPz+BQegadW9fns5eGMfK5iYyZ+FuuMmGhIbz10CXEJ6Zwx/MTizS/cnzbuHYVrz52R5HLBQUFU65CpXyHOXQojeTEBAAaNGl5RPUD+8OMkWWiAvY3BH77/OyfvuHjUS+QkWF/jDEyKpq4g/tZOn82S+fPDhjQNqz5i+cfuJlDaakEBQUTHBLCprV/887zDxIXezf9z74oV5lDh9L4aNQLRESW4aqb7jmCOZXjiSfDC5CYR7coYBhQGVB4KcDiVVvY8uoU5i3fwD+bdhObYH+4LSoyjCH9O/DsyCFUq1SWL169nrZDniTO59eZ7x12enZwefeL33hp7E/siInFGEP75rV56e4L6NGxMc/feR5zl21g6d9bC12vkJBg0tMzmPbbX0z4YTGzF69jf2wiQUGGDs3r8MJd59OjY2MevOEstu7cz0dfz881jq5tG/DWg5cQEhLMtzP+5M4XJrF9z0EqlY/i8RGDuf7Cnjw8fCBrNu5k8nT/1pRn7ziPyIgwxk9dxK3PTCAl7RC3XNqXF+++gCdvHcxnUxeSkJTqV+a+686gaf3q3PPSl2zbdaDQ8yonhqjoctRr1Ix6jZtRr1EzPn//dWIPBLw1D4DKVavz5mc/5DvMJ6Nf5tepkwgLD6db3zOOuH4n9+rP9XcWrdUGYP3qlXz09gtkZmbQsXsfrhh+F5WqVCchLpYvPx7NrB+m8M34D6hVtwEn9+rvV/aL/3uLQ2mpdD/1TK4ecT+hoWFM/3Yin3/wOl+Oe4cepw3MFai+m/Ahu7Zv5dLrR1K52klHPL9yfPDkPS+O47yS9Q94D4gErgEmAA1LtHIeMX7qIl7/5FcWrdycHVwAEpPT+Oy7hVz70EcAVK9cjoG9WvuVvfzskwH47Y913PH8RHbExALgOA7LVm/j/NvfJT4xhaCgIIb0a1+kem3cFkP7C57m4rveZ8ovy9kfa3NqZqbDkr+3ctaNb7Fi7b8A3HPN6XmO45nbzyUkJJiVa7dz+X1j2b7nIAD7YxO57ZkJ/Dz3bzvcyCEEBfn+CnUo3ds1JD09g5HPTyQpJY3MTIe3PpvJstXbKBcdSde2Dfym1axBde66uj9L/t7KOxNmF2le5fjXrFV7Rn0xnXuffZuLr72Vbn1OJyQ0/186L4y0tFQWzPoJgE6nnEpUdNmjHmdRTfzwLTIzM6hdvzE33/8slapUByC6XHmuHnE/rTt2A+CL/3ubTLdlBuyltHWrVxAUFMxVN99DeEQkQcHBnHHepdRr1IyU5CQ2rFnpN60dWzcx7ctPqN+4OQMGD/3vZlJKLU+GFwBjTCVjzNPACmwLUkfHce5zHGdPCVftuLBo5ebsz7WqV/Drd1KV8gABW1TiElJYt8V+DdFlwos03e17DrJha0zA/ofSM5jw/WIAGtWtSoWykX7969eqTI+OjQF4/ZNfSU/PzDWOlz78GYB6NSvT0x0WoGLZMgQHB7H3YCLxiSl+ZTZstfNTtWK0X/e3H76U4KAgRjw1nsxM/VC5+AsKDj4m410ybxaJCXEA9Dnj3GMyjfzs2bmdtav+BOCs8y8nJCR3I/7ZQ/8HwL49O/nnr8MtnEkJcTiZmUSXK09kGf/tqXrNOgDExR7M7uY4DuPefp7MzEyuvvWBY7ZMxVs8GV6MMS8Bi4F4oI3jOI87jqP2+mLUo0Oj7M8b/93r12/Tdvt3hxZ18ixbLjqCJvWqAbCkCJeMCislLT37c3Cw/yrcr1vz7M/T3RaWnOYt20Cc29rUv3uL7O4H4pPIyMikSoUoykZF+JVpUMe+PijmQEJ2t2vOO4WeHRsz6vNZLF/z7xHOjUjR/fbzt4A92Ddv0/E/n/6q5QuzP7fp1C3PYZq2bEdEZBkA/lp2ePgy0eUwQUEkxMWSnJTgV2bPTrsdlSt/+IRp9k/fsHbVcgacezH1GzdHBDwaXoC7gJrAw8AOY0yc+y/eGBNXwnXzrLDQEOrWqMTwi3sz9umrAFi/dQ/fz/7Lb7j3J80BoE+Xprx2/1BqVi2f3a9989p89cZwykZFsHDFJiZMW1zs9ezduQkAO2Ni2XfQ//anVo1rArB7X5xf0PCVmemwdvNuAFo0qpHdPTnlEAv+3EhISDCv3X8RkRGhBAUZbr60D51a1iU+MSW7RapapbI8ffu5bN25nyffKfiJDZHismfndtasWAJA79OP/mHL1X8u5r7rL+S6Ib0YfuGpPHzzZXz23qvs2h74xOPfzRsBKFehYsAbi4OCg6lRpz4A27ccvuk/PCKCJi3akJmZwSejXyY1JYXMjAymf/sFm9evISKyDI2a20vVsQf2MenDUVSuehLnX3HDUc+rHD88ecOu4zheDV2l0oEFrxERHpqr+7xlG7j6wXGkHUr36/7uF79Rq1oFRl7Vj+EX92b4xb2znzYKDwtlZ0wsL/3fzzz73g9kZOS+bHM0urZtwOC+bQH4cMq8XP1ruEFqx57YfMeT1b+GT/ACeOiNb/jxvdu4/OyuXHxmZw6lZxAZYe9ReGLU1OzLSS/dcwGVykcx7OGPSUpJO7qZEimC36d/h+M4BAcH06PfoKMe3/69ewgKCiayTBTJSYn8u2UD/27ZwMzvv+KyG+7gtEEX5CpzcL+9tFuxcrV8x12xclU2+QyfZeg1t/L8Azczb8YPLJj1M8EhIRxKszfCn3/ljdmXk8a/9xqJCXFcf9djhEdE5hy9nMA8GV6keO3eF0d4WCjRZcKz71GZtegfHnrjmzyfnnEch0fe+pbVm3bx2n0XUTYqwu8yS0R4KOWiI4iKDCMl9VCx1bNKxWg+eu5qgoODWLdlD6+Om55rmKz6JxcQKLICR9kc9+QsXLGJAcNe59Gbz6Zr2/qEhgSz5O+tvPHxr0z6yZ7tDjilBUPP7Mzkn5fy4xz7MxM3XdKH6y/qScPaVdh3MJEpvyzj8VFTcz2ZJHI0MjMymPPLVADadulBhUqVj3hc9Rs3o0HTFrQ/uSeVKlcjKDiY1JQUVi6Zz8QP32bPzn/5+J0XKVu+Al16+r97JSXZvv8pLDz/e9rCwiP8hs/SuEUbHnzhXb76dAzrV68kIyOD+o2bc+b5l9Otj70Rf8Uf81n423S69OxH+5N7AjD924nMnDaZPTu3E12+Al16nMr5Vw7P91FvOT4pvAjNBz2W/blqxWguO/tk7h12Br9/cjfPf5D7fSqVK0Tx2YvD6NOlKb/MX82zY35g1YYdRIaH0rVtA56+fQg3Du3NgO4tGDDs9eynkY5GVGQYX75+I/VqViYuIZnL7x1LYvKxafH4Y9UWzrllVJ79IiNCeeOBizkYn8TdL30JwLMjh3DH//qzffcBJv64hLbNanPLZafSsWVdBlz3RrG3PsmJa8WS+RzYZ1sxjvZG3QHnXJyrW3hEBJ17nErzNh15fOT/2Lt7JxPGvkXnHqdhTOD3vRyJhs1acfdTb+bZLzUlhU/eeZHIqGguv9G+92nC2Df58avPqFi5Kl37nM7WjWuZ/u1ENq1bwwMvjM73RXpy/NHlF/ETcyCBNz6Zwbm3vIPjwIM3nMVZOR6Vfv/JK+nTpSm//bGOwTePYv6fG4lLSGH3vni+nbmCfte8SsyBeBrWqcpTtx39kxBlIsKY8tZNdG3bgPjEFM67dTQrA7y9N6ulI+tST37jBIgvYsvIQzcMpEHtKjz8xjfs2htHk3rVuP3K09i1N47ul77ADY99Ss8rXmT24rV0b9+Iq87J+2ZGkSORdaNuxcpVadup+zGbTnS58gy++BrAPi20ZcM/fv2zbsRNS81/+0lLTfEbvrC+Gf8BMbt3MPSaW6hQqQo7/93CT1PGU65CJZ5482Ouv/NRHnt9HM3bdGT96hX8Pn1qkcYv3qfwInn6Y9UW5i3fAMC1F/TI7t6sQfXsMPPGJ7/mWTbmQALjpy4C4Nx+7Y6qHlnBpVenJiQkpXLebaOZt3xjwOF3uq08NauVDziMb/+dRWgVat2kJrddcRrzl29g7OS5AAzu25agoCAmTFucfYNwenomb4+3bxYefGrbQo9fJD9xB/fz52K73vXsf/Yxf2S4cfPDb9CO2eV/slChkn367sC+/N9MkdVKlDV8YWzbtI6fvh5P4xZt6Xum/YmEZQt+w3EcTjn1zOwbhENCQjj93EsAWDpf71g60Si8SEBZN7U2qlMlu1uLhoefzsn5CLWv9e67WqIiw6lW6cheoJUVXHp3bkJisg0uc5duyLfMqvU7APtyvSo53smSJSjI0LS+faHW6g07C1UXYwyjHr4UB4dbnvo8u3v92nbZbPzX/4bErPlvUKsKIsVhzq/fk5GejjGGXgMGl2hdate37wKNO3iAuNi831KRmZHBzm2bAahVr0Gew+Qqk5nJh289B8Zwza0PZF+q2uOGp2o1avsNX72WfV1DzK4dRZ4H8TaFFwmoQS17M2BC4uGm4czMw/dv1K0R+LdXfAPLkdy0WiYijK/ftsElISmVIbeOZs6S9QWW+3XBmuzPA05pkecw3ds3pFy0fXLhl/mrC1WfG4f24uS2DXh13C+s3rgrV//I8LAcf9untxz04jopHr/9/B0Azdt2olqNWsd8euvXHH5FQtXqNf36tWrfNfvzyiW5f6IDYN3qFdk36rbu0DXPYXKa8f1kNv6zioEXXEmterlflp6W5r8vOeRetirm23HEAxReTkC+r8QPpO/JTencuh5gfwYgy/I1h3/N+fqLeuVZtkxEWPZPCKxY+2+RHyXOCi5Zl4qG3PpOoYILwObt+5i71A57+5X9CAnJvYrf7f6swJYd+5iztODx1qxansdvGcz6rXt4/oMf/fptcV/Yl7WsspzcpkF2fUSO1tpVy9n17xYA+hTDu10cJ/9QnRAfy9SJ4wCoVKUadRs18+tfrUYtmrayl4R//Go86enpOUfB95PsT4xUrlaDZq07FFinA3v3MPnj0VSvWZvBl1zj1y8rPG1cu8qv+/p/bMCqkiNcyfFP4eUEVLt6RRZMuJ9hF/Sgfq3KOfpV4O5rBjDptRsJCgpi38FE3vpsRnb/rTsPMHW2/d2Rs/u0YexTV9HAvXQSEhJEt3YN+PmD22novpH2jU9mkNN7T1xB8rK3SV72dq5+kRGhfPXmcHp1akJ8YgpDRrxT4KWinB564xvS0zNo16w2nzx/bfZL9CqWK8PrDwzlzJ6t7HCvf12oV/q/et9FlC8bya3PTCA1zX8nnbUszj2tHeeeZnfmLRqexD3XDrD9Z/n/RouceBLj44iPPZj9z3Fs62Vaaopf95yPE/ua/ZO9UTeqbDk69Ti1UNN9/9UnuXpQV64elLvVY97MH3jr6ftYPHcGcQf3Z3dPS01hyfzZPHXnsOz7XC4edhtBQbkPFUOvuZWgoGC2bVrH6Bce4sBee/9LQnwsH496kRV/2BaZi68dUaj7cz599xWSkxL53y33Exbm/wh2h672RGnJvFn8MXcmjuOwfctGvp9oA1LHbr0LsUTkeKJny05Q7ZrV5u2HLwUgNe0QcYkpRIaH+f0W0aZ/93Lp3R+we1+8X9nhj3/KN6NuoVPLulx29slcdvbJJCanEhYSQmjo4Z3Uax/9kn3jbmGd178Dfbo0BSAkOIhPXxyW7/CX3v0+C/7c5Ndt4YpN3PrsBN584BKG9GvPkH7tORCXRPnoiOyd8NPvTsv1i9J5ObtvG87t155Pv1vIrEVrc/X/Z9NuRo2fyS2XncqEV64nKTmNMpH2EtKiFZv4+Nu8m9TlxPHobVexb0/ue6t+mPwpP0z+NPvvHv0G5fnrzslJiSyeY2+OP+XUMwkthh92zMzIZMn8WSyZPwuA8IhIQsPCSEpIIDPT/ohiSGgYl153O117D8hzHI1btOF/I+7j41EvsGTeLJbMm0WZqLIkJyVkt+yce9l1uX5ROi9L589myfxZ9Og3kJbtu+TqX7NuAwacM5Tp307k7WfvJyw8PPtJp4bNWtGzhO8Bkv+ewssJaGdMLJffM5ZenZvQpU19alQpR+UKUWRkOGzduZ8Va7czddYKvvjhjzxfMrfvYCJ9rnqZKwZ35fz+HWjbrDaVypchPSODbbv2s+DPTYydPCffp4ICCTK+v/IcVuAjz6F5/CAcwLgp81m+ehu3X9mPXp0aU6ViNHv2J7BoxSbemTCb2YtzB5GcoiLDePXei9h7IIH7X/0q4HB3vzSZLTv2M+yCHjSoXYWdMbF8/etyHh/1XZ4/DClSFAtn/5z9yHHvYvoRxhbtOnHBVTexfs1Kdm7bbH9nKDGByDJRVKtZmxZtO3PqWedR9aT8L8f0OeNc6jVqxo9TxvPPyqXExx6kbPmKNG7Rhv6Dh9KyXecC65KSnMSnY14hulx5Lhl2e8DhLrvhTqpUq8nMH6cQs2s75StWpnOP07jgquF5/jCkHN9MQdc+j1eRHUacmDMu8h+ZMenpkq6CyHGve+MKJ+TtyrrnRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDxF4UVEREQ8ReFFREREPEXhRURERDylxMOLMSbYGPNLSddDREREvKHEw4vjOBlAkjGmfEnXRUREREq/kJKugCsFWGmMmQ4kZnV0HOe2kquSiIiIlEalJbx87/4TERERyVepCC+O43xkjIkE6jqO809J10dERERKrxK/5wXAGDMYWA786P7d3hjzbcnWSkREREqjUhFegMeBk4GDAI7jLAcalGSFREREpHQqLeEl3XGc2BzdnBKpiYiIiJRqpeKeF+AvY8xlQLAxpglwGzCvhOskIiIipVBpaXm5FWgFpALjgThgZInWSEREREql0tLyUs1xnIeAh7I6GGO6AItLrkoiIiJSGpWWlpevjDG1sv4wxvQG/q8E6yMiIiKlVGkJLzcCXxtjTjLGDATeBAaWcJ1ERESkFCoVl40cx1lsjLkN+Bn7UwEDHMeJKeFqiYiISClUouHFGPMd/o9ElwFigbHGGBzHOadkaiYiIiKlVUm3vLxcwtMXERERjynR8OI4zuysz8aY6kAX989FjuPsOZbTDm/Z7ViOXuSE16F+hZKugogcp0rFDbvGmKHAIuAiYCiw0BhzYcnWSkREREqjkr5slOUhoEtWa4sxpirwC/BlidZKRERESp1S0fICBOW4TLSP0lM3ERERKUVKS8vLj8aYn4DP3b8vBn4owfqIiIhIKVUqwovjOPcYYy4AegAGeM9xnCklXC0REREphUpFeAFwHGeyMWY6bp2MMZUcx9lfwtUSERGRUqZUhBdjzI3Ak0AykIltfXGAhiVZLxERESl9SkV4Ae4GWjmOs7ekKyIiIiKlW2l5omcDkFTSlRAREZHSr7S0vDwAzDPGLARSszo6jnNbyVVJRERESqPSEl7GADOAldh7XkRERETyVFrCS7rjOHeWdCVERESk9Cst97zMNMbcYIypYYyplPWvpCslIiIipU9paXm5zP3//hzd9ai0iIiI+CnRlhdjTBdjzEmO4zRwHKcB8ATwFzAV6FySdRMREZHSqaQvG40B0gCMMb2B54CPgFjgvRKsl4iIiJRSJX3ZKNjnJwAuxv6m0WRgsjFmeQnWS0REREqpkm55CTbGZAWoftjHpbOUdLASERGRUqikA8LnwGxjzF7s7xr9DmCMaYy9dCQiIiLip0TDi+M4zxhjfgVqAD87juO4vYKAW0uuZiIiIlJalXTLC47jLMij29qSqIuIiIiUfiV9z4uIiIhIkSi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKcovIiIiIinKLyIiIiIpyi8iIiIiKeElHQF5MQxcnArHr+kQ/bfFS7/tNBlX732ZK7t1xSArTEJtB35dZ7DndmhFj1aVKd9g0rUrhxF5XLhRIQGsy8+lb+2HuDrBVuYMGcTGZnO0c2MSAlLTk5myR+L+HvVKlav/pvVq1axc+cOAIbfPIKbbrk1YNk/Fi9i/ry5rPrrL7b/u40DBw+QnJRE2XLlaNSoMaf1H8AFFw4lIiIi4DgyMzP5YdpUvv3ma/5Z/TcJCQlUqFiRDh06csllV9Cpc5din2eRLAov8p9oXKMc953f5ojK9mxRnatPbVKoYR+9uAMt61TI/jsuOY2MTIcaFctQo2IZBrSrxfWnN+OiF2cSE5dyRPURKQ3+WrmCW4bfcERlx304lt9nz8r+OzKyDKGhYRzYv58/9i/ij8WL+OyTj3hnzAfUr98gV/mkpCTuvP1W5s+bA0BwcDBRUdHs27uXn3/6kek//8SNN92Sb4ASORoKL3LMGQNvXd+NyLAQFq6NoWvTqoUuGxkWzJvXdSM9M5O/Nh2kY6PK+Q7/7aKtvPvTGhaujWHLngRSDmUAcFKFSK46tTH3nd+G9g0qM3r4KVz44oyjmi+RklauXHlatGxJixYtad6yFS+/8Bx798YUWK5bt+6c0qMnHTp2om7dukRFRQNw8OABpn0/lTdefZnt//7LHbeNYPLX3xEU5H+HwVNPPMr8eXMICgpixG0jueSyy4mKiiYuLo6PPxzL+++9y7vvvE29evUZePbgYzLvcmJTeJFj7sbTm9G9WTW+mLOJTbvjixReHhnanoYnleWlr1dSq1KZAsPL81+tyLP7roPJvDhlJeGhwdx1bmv6t6tJzUpl2LE/qUjzIlJadOzUmd/nL/Lr9uZrrxSq7BVXXZ1n9woVKnLZ5VcSFhrGU088ysYN6/lz+TI6dOyUPcy6dWuZNvU7AC67/EqGXX9jdr9y5cox4vY72L17N99+M4VXX3mRAaefQWhYWBHnTiR/umFXjql6VaN4eGh79sWn8OCnfxSpbOfGVbjxjGas2xHHy1+vLJb6/LF+b/bnGhUji2WcIiUhODj4mI27bbt22Z93797t18/3ctP/rh2WZ/mrryEm8PUAACAASURBVL0OgJg9e5g3b27xV1BOeAovcky9cV03oiNCeejTJeyLTy10ubCQIN6+oRsGwx3/t5DUQ5nFUp/uzaplf968J6FYxilyvFm6ZEn25zp16vj127nD3hRctmxZqlWrnmf5uvXqERJiG/bnzvn9GNVSTmS6bCTHzFWnNqZv6xrMXLmTCXM2Fansvee1oXmtCnw0cx1zVu8uuEA+osJDqFs1ikt6NWTEwBYAfP7bhiKFKZHjXUpKCrt372L6Tz8yZvQoADp17kKr1nnfaJ+RkRFwXI7jkJlpTzjWr1tb/JWVE57CixwTNSpG8tSlHUlKTWfk2IVFKtu2XkVuP7sVuw8m8+j4ZUc0/c6Nq/DLE2fm6p6ekclnszZwz7jFRzRekePJ3pgY+vXtmWe/Pn1P5alnns/VvWatWoB94mjHju3UrFkr1zAbN2zIDi8xe/YUY41FLF02kmPitWFdKR8VxvNfrWBLTOEvzwQHGd6+oTuhIUHc9/FiYpPSjmj6h9Iz2X0wmd0Hk0lLP3yG+OGMdbzw1Yrsp5BETmRBwcFUrlyFypWrEB4ent399DPO5I677qF8hQq5yvTs3Sf78/vvjs5zvB+8d7h7QqIuz0rxU8uLFLuhPRpwZofarNi8n1HTVhep7B3ntKJt/Ur8uPRfvl649Yjr8Ofm/TS7ZTJgH9WuXy2am89qwbX9mnBJz4bc8M5cflj67xGPX+R4UKlSJWb8Zm+odRyHPbt3M+mLz/n4ow+Z8euvPPDQI1w49GK/Mk2aNOXMswbx4w/f89XkSZQpU4bLrryKatWqs3PHDsZ9OJbpP/9ESEgo6emHCDI6R5bip7VKilWVcuE8d2Un0jMyue2DBUV6k22zWuW5Z0gb4pMPcdeHiwouUEiOA5t2J3DPuMU8On4ZZSNDee/mHlSvoKeNRLIYY6h+0kmMuP0Onn3hZdLTD/HMU4/zz5o1uYZ97Imn6HZKDwA+/eQjBp7ej87tWzN44OlMnvQF7dp34KyBgwD7+LRIcVN4kWL1xCUdqVw2gnEz1rFuRxxR4SF+/0JDDq9y2d2CbbeXr+5CeGgwr3zzFwcT03KVDXZflGXM4bIhwaZI9fvgl39IScugbGQoF3avX2zzLXI86T/gdGrWrEVmZiZTvvoyV/8yUf/f3p3H6VQ9cBz/nNlX+8yUdeyM7GtS2qyFRIVUSimlffkpftosrRIivxKyRJIWqZAipcJYspTs+zIYZl/P749n5mHMDDMM0x3f9+vl9Xrm3nPuc66XZ3yfc88SyPgJH/LmyFFcf2MbKoWHU7ZsOZo0bcbzg4bw0ZRpREW5liWoFB5+kVsvlwI9NpICVSnEtVLn/W1qcn+bmmcsu/ejHgCM/3YTz09b5a77Uo+GWfZAOl2FMkHuus9PXcn477J/M8xNUko6x+KSuNwngMphQXmuJ3KpCQkNZd++vezetTPH8x4eHrRt14G27TpkO5eSnMyf69YA0KBhowvaTrk0qedFLilBfl6UCXZtNhebmFrIrRH5d7LWsnePa0xYQGBgvusvWPAdsbGxeHl50eEmbQ8gBU89L1Kgbh628IznB95aj4Hd6gHZd5XObafoTOMevJJe11TNdVdpTw9z1jE2j90U4X50db7rx4g4UWpqqnsBudx8MXeOe4+kpk2b5ev6hw8fYtTItwC45dZuhIXlvJCdyPlQeJEi4/arKtO5aUVmLN3K8s2HiDrhWoTOGIgoX4IH2takz/Wu3amX/32IRWv3FWZzRc7biePHSUs/Oe0/c22VxIQEjh076j7u6+Pr7kFZHbmKcWNHc2v322jWrAVhl13mLrdz5w7mzvmMqVMmAVChQkU633JrtvdduuQndu/ayTWtr6NsuXJ4enoSHx/Pkh8X8+47b3Po4EEqV6nCU888d0HuW0ThRYoMY6BD4/J0aFwegNjEFPfgXF/vk/vALFm/n3tGa8lycb47undl37692Y5PnjSRyZMmun/u3KUrrw4/ueBc5KqVRK5y7TXm6+tLQEAACQkJJCYmusvUrFmLd8a8h5+fX7br79q5gzdfH8Ebrw3Hy8uLgIBAYmJOYK2r57Nuvfq8O2ace7dqkYKm8CJFxver9/LYh7/RqnYYV1QsSWhxP0oE+pCQnMaOQ7Gs3naEOct3sFA9LnIJi6hTh6EjXmflij/YtGEDUVFRHD8ejbe3DxUqVKR2RAQ3tGlHm7btct38sUXLq+jV+y5WR0ZyYP8+YmJiKVOmDLUi6tCh40106HgzHh4aUikXjslMypeaEndOuzRvXOQiOTCld2E3QaTI8/Mif+tFFBGKxiIiIuIoCi8iIiLiKAovIiIi4igKLyIiIuIojg0vxpgBxphiGa8nGGP+MMbcUNjtEhERkQvLseEF6GetPWGMaQuUA/oDbxRym0REROQCc3J4yZzq3AGYZK1dxVnuxxjTzxiz0hizMnnL4gveQBERESl4Tg4va40x84FOwLfGmCBOBpocWWv/Z61tYq1t4lPt+ovSSBERESlYTl5h916gMbDFWhtvjCkD9C3kNomIiMgF5tjwYq1NM8ZUAdoAwwB/nN2TdFHVDy9F+0blaBBemqqXB1Mm2I9gf29iElLYvP84C9fsY+KizUTHJWere1WtUK6rezkNq5QmPDSI0sG+BPp6Ex2XzF97o5m3cjdTFm8hMSUth3c+s1a1w5g3uE2ey4+Ys5bXP/8zy7GKZQJpUTOUBpVLUT+8FHXDS1LM3weAeo/PZVdU3Bmv2aByKV7oXp/mNULw9vRg055oRs/byJd/7Mq1ztURYXw9qA3fRe6hx9s/5bn9UrQlJCSwauUfbNywgU2bNrJpwwb273dtT/HQwwPo/8ijZ73GkagoJn30IUuX/MiB/fvx9fWjarVqdO7Sla7dumPMuS2w+t8XBvLVl3PPWm7V2g257kKdmJjI7Fkz+WHRArZt3UpcXCx+fn5UqFiJq1pdTc9evSkTEpJj3Q3r/2Tc2NGsWR1JamoqVatVp899fWnbrkOubfnj99944L57uKb1dYwZ937eblSKJMeGF2PMWMAbuAZXeIkD3geaFma7nOLO1lXp17am++eE5FQSUlIpFexLi+BQWtQIpX/7WvR8+ydWbInKUvfRmyNo37C8++fYxBSSUtMIKe5HSPHLuDriMvq3r0W31xez9UBMvtqVnJrGweiEM5YJ8PUi2N8bgMitR7KdH9itHr2uqZqv983UuGpp5g1ug7+PF6lp6aSkpdO4ahmmPH4Nz0z+gw8Xbs5Wx8fLg5H3NicmIYVnJv9xTu8rRdP6P9fxyEP9zrn+xg3r6d+vL9HR0QAEBAQQHx/H6shVrI5cxYLvv2P0e+Px8fE55/fw9fUlKCg41/O5haN9+/by4P33sWvnDvex4OBg4uLi2LRxA5s2bmDWJ9N5d+x4GjfJ+mt53bq13N/nLpKSkvD09MTLy4sN6//k2aee4Oigo/TodWe290tOTmboKy8SEBDAC4OHnNvNSpHh2PACtLTWNjLGrAaw1h41xpz7J/gSE7k1isHTY/lt82H+2Xec4/EpAAT6etG5WUVe6dmIkOJ+TH+qNU2e/ooTCSnuukvWH2Dxuv389vchth2MITYxFYCSQT7c1rIyL/VoSHhoMNOebE3LgfPIz/ZZf/wTRc1H5pyxzMynr6V9o/LsPRrHD+v2Zzufnm7ZdiCGtTuOsmb7EYwxvNSjYZ7e/9VejfD38WLWsm08+dHvJCan81C7mgy/qwlD7mjAzJ+3ue830zNdrqB62WI8P3Ule47E5/1m5ZJQrFhxakdEULt2BLUi6vDW6yOIijp81noxMTE8+vBDREdHU7lKFYaNeIM6V9QlJTmZOZ/N5s3XR7D812W8+dpwBg156Zzb1659xyw7TufV4Of/w66dO/D29uaZ556nc5dbCAgMJCU5mV+W/cywoS9z6OBBnn3qCeYv+CHL7tTvvPUGSUlJ3NSpM4OHvIyvry8zpk3lrTdGMHrU23Tq0iXbjtQfTBjPzh07ePY/z3N52bLnfL9SNDg5vKQYYzzIGKRrjCkNpBduk5xj5rLtOR6PS0rlk5+3cSA6gbkDbyC0uD/tGpZj9q873GXGf/dXjnWPxSbzvwV/k5SSxrv3t6B2+RI0qx7C75vP/os6ry4r4c+N9V2/uGYs2UZ6DsnosQ9/z3K8Ve2wPF3b38eT5jVCSE1L55nJK4hPcj32GvfdX9zRqgr1K5eiafUQfvzzZGCqUbYYj3eqw+ptR5jw/d/nc2tSBDVq3ISfl2ftjRv9ztt5qjtl0kSiog7j5+fH2PH/o3z5CgB4+/jQo9edxMXFMnrUSOZ89il33n0P4eGVC7z9udm3by+rVq4AoO8DD2bpKfH28eHa62/APyCAfn37cORIFKtWruCqVlcDrkdpa1ZH4unpyQuDXT0pAHfd04d5X3/JX5s2snbNGlpe1cp9zW1btzJp4gdE1KlDzzvvumj3Kf9eTh4j8h4wBwgxxrwMLANeL9wmFR0rT3lUVLZUwEWreza9rqmCl6cH6emWaUu25Fgmp0CTFyUCffD08OBITBIxp/Q0AWw9eAKAMsG+WY6P6tscTw/DExN/P+f3laLL09PznOvO++pLANp16OgOLqfq2as3AQEBpKWlMX/e1+f8Puci6vDJLyQRda7IscwVdeu6XyfEn+yRPHHiBOnp6ZQoUZKgoKy9KxUrVQLg2LGj7mPWWl59eQjp6ekMeenV8/o7laLDceHFGDPfGBNurf0YGAy8BRwDbrPWzizc1hUdV9YMdb/efig2f3VrnVL3YP7GvJxN79bVAFiy4QA7D5954G1+Rcclk5aeTulgX/eYmkyVQ11jAqJiktzH7r6uGi1rhfH+93+xdsdRRArKju3b3AN7W119TY5lAgIDadS4CQDLf/3lorUNoNwpYWrjhvU5lln/p2sgvYeHB7VqR7iPFytWDA8PD6KjjxEbm/V3y57drkHxJUuWch/7/LPZRK5aSa/ed1M7ok6B3YM4m+PCCzAZWGCMGQRstta+a60dZa3N+RMkeebj5UHFMoE80KYGE/q3BGDrgRN8F7nnrHX9vD2pEhbMU53rMLRXYwB+2XSQNdsL7j/1qyPCqHKZK0R8/FPOvS7nIyE5jd83R+Hl6cGb9zTF38cTD2N4sF1NGlYpTUxCirtXKaSYHy/3aMjuqFiGf7a2wNsil7Yt//zjfl2tWo1cy1WrVh2AbVvP/fPw++/L6dSxHU0b1qVls0Z0u6UTb4wYxs5TBuKernTp0txwY1sAJn4wgZkzphMf5/oykZKSwk+Lf2DwC/8B4O577qV8hZNhx9/fn/oNGpKWlsaIYa+QkJBAWloa06d9zMYNGwgICKBe/QaAa6bVqJFvcfnlZXl4wGPnfI9S9DhuzIu19lNjzDfAEGClMWYqp4x1sdaOLLTGOdSBST3x88neFbv870M88N4yklNzHkoUWtyPzeO653ju21V76D/h1wJtZ2avy5GYRL5ZubtAr53pxU8i+XpQG3pcXYXuLcNJSUvH38f1MRk2e437cdKIu5pQMsiXB8f/6h4bI1JQDh0+5H4dGpb7mK3Mc7GxscTHxREQGJjv9zp44ACenp4EBgYRFxfLln82s+WfzXw66xOeG/gCt/folWO9l14ZSkJCPL/+sowRw15hxLBX3LON0tPTqV6jBv0ffpRbu9+Wre6TTz/L/ffezbyvvuTbb+bh5eVFUpKrV3PAY0+4Hye98dpwTpw4zrDX3nCPjREBB4aXDCm4pkb7AsFooO55OXQ8AV9vTwL9vAjycz0uWbrhAEM+iTzj7Jm0dOue1lwswNv9n/zc33Yy/LO1Oa4Rc66KB3jTuZnr29uny7bnGqjO14otUXR8dQGDbqtP0+pl8Pb0YPW2I4z5ZiOf/7YTgBvqXU73luHM/W0nC9bsBaBf25rcd2N1KocGcyQmia/+2MXQ2WuyzUwSyYvMXgwgyyyd0/n5+btfx8XnL7zUioigzhV1uebaawkLuwxPT08SEhL4ZdnPjHr7TXbv3sWwV1+mZMlStGnXPlv9YsWLM/LdsYwbO5qpUyZhrSUm5uRj4vj4eI4dO0ZaWlq2cSr1GzTkoynTeW/Mu6xbu5rU1FQi6tTh7j596dDxJgB++Xkp3337DW3bteea1tcCMGP6VGbP/ITdu3dRomRJ2rRtx4DHnsg2M0mKPseFF2NMe2Ak8BXQyFqruannqd4TX7hflynmS49WVXi6yxUsfqUDb33xJ8PnrMux3pGYpCzTmsuWCuDeG6ozoGNtbmpSnmcnr2DKjwXzeOe2qyq7w9GFeGR0qshtR+j2es57X/n7ePL2vc04HpfMwKkrAXilZyMeuzmCvUfjmLN8B3UrleSh9rVoUKUUN726kLR0DeSVf587e9+d7Zi/vz83tmlLk6ZN6Xl7N/bt3cvIt97gxrbtsq33smnTRp4Y8DCHDh3k9h69uO2OHpQvX4EjUVEs+Wkx498bw+hRbxO5agVjxk3AwyPrKIW69erx/gcTc2xbQkICw4a+THBwMM89PwiAkW+9zpRJHxEaFkaHjjfz91+bmDFtKhvXr2filGm5LqQnRZMTx7wMwjU4d6CCS8GLOpHE2Pmb6Pb6YiyW526tR7uG5fJUd9/ReIbNXku/937Bx8uTkfc144qKJQqkXXdlPDJaseUwm/YcL5Brnov/3FqP8NBgXpy5moPRCVS7vBgDOtbmYHQCrQfN5+EJy7nuv9/y88YDtKgRyp2tz22xPLm0ndqDkpiYmGu5xMSTCzoGBuT/kVFuSpQoyQP9+gOuadF/bdqY5XxcXCyPPPgABw7s54EH+/PC4CFUr14Df39/yleowJ133cPbo8ZgjGHZz0v5Yu6Z12463fvjxrJ3zx4ef+oZQkJC2bF9Gx9PnkTp0mWYOXsurw5/jemzPqNJ02asWbOaL+d+XmD3Ls7guPBirb3aWruhsNtR1EVuO8Jvf7umQ/a5rnq+6n69cje7Dsfi6eHBXddWO++21A8vRf3KrtkHH/+49byvd67qVCjBIx1q89vmQ0xe7BpQ2bFxeTw8DJ/+sp2oE65n9qlp1r0Wzk2Ns09xFTmb0JCTM/YOHTyYa7nMc0FBQec03uVM6jVo4H69Z3fWMWbffP0VR464Bq/f3ee+HOs3b3Gle5bRDwsX5Pl9N//9F9M+nkyDBg3pftsdAPy4+AestdzUqTOlS5cGwNvbm95335NxflGery9Fg+PCi1w8+465OrYqX5b/58n7j7m+EVYOy33Z8bzqfa2r9yI2MYXPl+847+udC2Ng1P3NsVie+PB39/HwUNffzelTwjO3Rcg8L5If1aqf/MKwZUv2LSlOnnOF6CpVz/9LQn5s3er6ElGyVKlsa7WcKnPdlr17zz5jESA9PZ1XXhoCGP778qvuR1V79rjCU4UKWb8MVKwYnuW8XDoUXiRXmf/xxibkf9BppZCMuokpZyl5Zn7entzWMhxwDQSOSyqcAbD331iDptVCeHfeRv7am/2xVeZ4nEx+3q4BihaNd5H8C69chcsvd60k/cuyn3MsEx8fT+Qq17irK1teVeBtWLd2jft1ufLls5zz8HCFiuhjx4iPz/3p/ZEoV+9MQB4fac36ZAZ/rlvLvff1dU8DP1ViYlKWn5OSXI/UDOe2OaU4l8LLJcgjD7vQXlPnMhpXKQPAsk0nu609Pc5et3frqlxW0jULYtnG3Lu886Jzs4qUCHStavtxAQ3+za/LS/oz+PYGbD1wgre+yLqD9c6MBfwaVS2d5XjT6mUyzhfsQnpy6bi5cxcAvvt2fo49F7M+mU58fDyenp50vLlTvq5tz7Ia9PHoaCZ+MAGAsMsuy7LIHEDt2nXc1/ns01k5XuOffzazZnUkAPVPeQSVm4MHDzJ29DtUrFiJBx56OMu5cuVc4Wn9+qyTB9atda2xdHq4kqJP4eUSVL50AD8P70if66u7e0gylSsVwBOd6jDjqdZ4eBiOxiQx7ttN7vNX1gxl/n/bcEerytmW/q8SFsyLdzTgnfuaA7DtQAwzlm7L9v7jHryS6Om9iZ7e+6xtvTtjzMymPdHZdrfOjZenoVSQr/tPsYCTq+WWCPTJcs7L8+xh7I17mlI8wIcnP/qDpJSsU7S/zVjAr1PTCnRq4urSrlWuOE92dv1yn79K3dmXuhPHj3Ps2FH3n/R017+hxISELMdPnR4NcM+9fSlTJoTEhAQG9H/QvZJtSnIyn86cwXtj3gWgW/fbc9zX6L8vDKR+nZrUr1Mz27l5X3/Jk48PYNGC7zly5OTO7ImJiSz+YRG9e93hHufy1DP/yTZT6Ma27QgJCQFgzLsj+WDCeKKjjwGuHqFv539D/wfuIzU1FW9vb3r0Ovtn/bXhrxIbG8vgF10bNZ7q2utuAOCHhQtZtHAB1lq2bPmHiR9OyHJeLh2aW3aJqlupFKP6ukJGUkoaMQkp+Pl4utd5AdhxKIa7Ry3l0PGssx1a1gqjZS3X4lgJyanEJaYS4OtFgO/Jf05/7jzKnSOXkJhy7gu4VQ4LomXGVgNT89Hr0qJGKPMGt8nx3NLhN2X5+eahC7P0LJ2uY+PydGpakU+WbmXphgPZzm/ed4L3v/uLh9rXYuqTrYlPSnX/PazYcpjpSwtvgLH8O9zRvSv79u3NdnzypIlMnnRyqnDnLl2z7O4cHBzMmHHv079fX7Zt3ULP27sRGBhIUlIyqamux7FXtmzFswNfyHeb0tPSWbxoIYsXLQTA3z8AX18fYmJiSEtzfWZ9fHx45rmBtO/QMVv9gIAA3hk9jkcfeZBjR48ydvQoxo4eRWBgIPHx8e6eHT9/f4YOe41KlcLP2J4fFy9i8aKFdOp8C81bXJntfJWqVenV+y5mTJvK0088ip+fn3sWVt169enS9dZ8/x2Isym8XIL2H0vgnneX0ioijCZVyxBWwp/Swb6kpVt2R8Wyfmc081ftZvavO7KFjzXbj/Dg+F9oVTuMBpVLEVrcn1JBviSlprHtQAxrdxzlqxW7+PL3Xee9UWHv1tXw8DAkpaTlugv2hRTo68XrdzflSEwig6ZH5lpu4NSV7DocS58bqhMeGsSBYwl89cdOhs5eS2qaxrzIuYuocwWff/kNH038gKVLfuLggf34B/hTrVo9Onfpyi23dsvWK5IXTZs359HHn2TtmjVs37aV6OPRxMbGEhgYRIWKFWnWvAXdb78jxw0hM9WtV48vvp7PpzM/4eelS9ixfTtxcbGu6dLlK9CsxZX07NU7y9YAOYmPi+O1YUMpUaIETz/3n1zLPTdwEGXLluez2TPZs3sPZcqEcGPGInXe3t651pOiyZzt2WdRVeLOaZfmjYtcJAemnP1RgYicHz+vS3O0ssa8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoyi8iIiIiKMovIiIiIijKLyIiIiIoxhrbWG3QSRPjDH9rLX/K+x2iBRV+oyJU6jnRZykX2E3QKSI02dMHEHhRURERBxF4UVEREQcReFFnETP4kUuLH3GxBE0YFdEREQcRT0vIiIi4igKL1LojDFdjTHWGFMr4+dwY8z6jNdNjDGjC7eFIs5kjEkzxqwxxmwwxqw1xjxljCmQ3/vGmJeMMc8UxLVE8kvhRf4NegLLgB6nn7DWrrTWPnbxmyRSJCRYaxtYa+sAbYCOwIuF3CaR86bwIoXKGBMEXAX0JYfwYoy51hgzzxjjYYzZYYwpccq5LcaYMGNMiDFmjjFmRcafqy7iLYg4grX2EK51XAYYF09jzJsZn5l1xpgHwfWZNMb8YIyJNMb8aYzpknkNY8wgY8zfxphFQM1CuhURvAq7AXLJuwX4zlq72Rhz1BjTCDh6eiFrbbox5kugKzDJGNMc2GGtPWiMmQG8Y61dZoypCHwP1L6YNyHiBNbabRmPjUKBLsBxa21TY4wv8IsxZgGwG+hqrT1hjCkD/GaM+QpohOsLRkNc/3dEAqsK5UbkkqfwIoWtJzAq4/XMjJ/fy6XsLGAIMAnXL9FZGcdvBCKMMZnlihljgq21MRekxSLOlvlBaQvUM8Z0z/i5OFAd2AMMN8ZcA6QD5YAw4GpgrrU2HiAj0IgUCoUXKTTGmNLA9cAVxhgLeAIWGJdLleVANWNMCK4em6EZxz2AK621CRe4ySKOZoypAqQBh3CFmEettd+fVqYPEAI0ttamGGN2AH4Zp7W2hvwr4du87gAAAvdJREFUaMyLFKbuwMfW2krW2nBrbQVgO1A+p8LWtSjRXGAksMlaeyTj1AJgQGY5Y0yDC9tsEefJCP3vA2MzPkvfA/2NMd4Z52sYYwJx9cAcyggu1wGVMi6xFOhqjPE3xgQDnS7+XYi4qOdFClNP4LXTjs0BXjhDnVnACqDPKcceA94zxqzD9W96KfBQwTVTxLH8jTFrAG8gFZiKK/wDfAiEA5HG9cz1MK4ezenA18aYlcAa4C8Aa22kMWZWxrGdwM8X8T5EstAKuyIiIuIoemwkIiIijqLwIiIiIo6i8CIiIiKOovAiIiIijqLwIiIiIo6i8CIiZ5Sxn82GjP1v1hhjmhtjnjDGBOShbp7KiYjkh6ZKi0iujDFX4loX5FprbVLGXjc+wK9AE2tt1Fnq78hLORGR/FDPi4icyeVAlLU2CSAjhHQHygI/GmN+BDDGjDfGrMzooXk549hjOZSLzbywMaa7MWZyxuvbjDHrjTFrjTFLL+L9iYgDqedFRHJljAkClgEBwCJglrV2yek9KsaYUtbao8YYT+AH4DFr7bocysVaa4MyXncHbrbW9jHG/Am0t9buNcaUsNZGX+x7FRHnUM+LiOTKWhsLNAb64Vo+flbGxn2nu90YEwmsBuoAEfl8q1+AycaYB3Bt0CkikivtbSQiZ2StTQN+An7K6CG559TzxpjKwDNAU2vtsYxHQX6nXyfzcqe8dpex1j5kjGkO3ASsMcY0OGXjTRGRLNTzIiK5MsbUNMZUP+VQA1yb8sUAwRnHigFxwHFjTBjQ4ZTyp5YDOGiMqW2M8QC6nvI+Va21v1trhwBRQIWCvxsRKSrU8yIiZxIEjDHGlMC1K/EWXI+QegLfGmP2W2uvM8asBjYA23A9Asr0v1PLAQOBecBuYH3G9QHezAhJBteYmbUX/tZExKk0YFdEREQcRY+NRERExFEUXkRERMRRFF5ERETEURReRERExFEUXkRERMRRFF5ERETEURReRERExFEUXkRERMRR/g8i0hcDYy/L1wAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 720x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "count = np.array(data.groupby(['Smoker', 'Status']).count())\n",
    "count = np.reshape(count, (2, 2))\n",
    "annots = np.array([f\"{v}\\n{v/len(data):.2%}\" for v in count.flatten()]).reshape(2,2)\n",
    "\n",
    "plt.figure(figsize=(10,8))\n",
    "sns.heatmap(count, annot=annots, fmt=\"\", cmap='Blues', cbar=False, square=True,\n",
    "           xticklabels=['Alive', 'Dead'], yticklabels=['No', 'Yes'], annot_kws={\"fontsize\": 25})\n",
    "plt.title(\"Number and percentage of alive/dead women after 20 years, according to smoking habits\")\n",
    "plt.xlabel(\"Status\")\n",
    "plt.ylabel(\"Smoker\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It is possible to see that the fraction of smokers and non smokers is quite balanced (in total, 582 smokers and 732 non smokers). As expected, there are less dead than alive people (369 versus 945).\n",
    "\n",
    "We can then compute the mortality rate for the two groups. For a population proportion $p$, confidence intervals are computed as $\\hat{p} \\pm z \\cdot \\sqrt{\\frac{\\hat{p}(1-\\hat{p})}{n}}$, where $\\hat{p}$ is the sample proportion, $n$ is the sample size and $z$ is the value derived from the standard normal distribution. For 95% confidence intervals, $z=1.96$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Mortality rate for smokers:\t23.88% ± 3.46%\n",
      "Mortality rate for non smokers:\t31.42% ± 3.36%\n"
     ]
    }
   ],
   "source": [
    "z = 1.96\n",
    "\n",
    "num_smokers = sum(data['Smoker'] == \"Yes\")\n",
    "num_dead_smokers = sum(np.logical_and(data['Smoker'] == \"Yes\", data['Status'] == \"Dead\"))\n",
    "rate_smokers = num_dead_smokers / num_smokers\n",
    "ci_smokers = z * (rate_smokers * (1 - rate_smokers) / num_smokers) ** 0.5\n",
    "print(f\"Mortality rate for smokers:\\t{rate_smokers:.2%} \" + u\"\\u00B1\" + f\" {ci_smokers:.2%}\")\n",
    "\n",
    "num_non_smokers = len(data) - num_smokers\n",
    "num_dead_non_smokers = sum(np.logical_and(data['Smoker'] == \"No\", data['Status'] == \"Dead\"))\n",
    "rate_non_smokers = num_dead_non_smokers / num_non_smokers\n",
    "ci_non_smokers = z * (rate_non_smokers * (1 - rate_non_smokers) / num_non_smokers) ** 0.5\n",
    "print(f\"Mortality rate for non smokers:\\t{rate_non_smokers:.2%} \" + u\"\\u00B1\" + f\" {ci_non_smokers:.2%}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Surprisingly, the mortality rate is sensibly higher for women categorized as non smokers. However, we are not taking into account an important information: the age of those people at the time of the poll. This result can be expected, for example, if the average age of polled non smokers was higher than the one of smokers.\n",
    "\n",
    "---\n",
    "\n",
    "Let's now include the age in the analysis. The following age classes are considered: 18-34 years, 35-54 years, 55-64 years, over 65 years."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Smoker</th>\n",
       "      <th>Status</th>\n",
       "      <th>Age</th>\n",
       "      <th>Age group</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>21.0</td>\n",
       "      <td>18-34 years</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>19.3</td>\n",
       "      <td>18-34 years</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>No</td>\n",
       "      <td>Dead</td>\n",
       "      <td>57.5</td>\n",
       "      <td>55-64 years</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>No</td>\n",
       "      <td>Alive</td>\n",
       "      <td>47.1</td>\n",
       "      <td>35-54 years</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>81.4</td>\n",
       "      <td>Over 65 years</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  Smoker Status   Age      Age group\n",
       "0    Yes  Alive  21.0    18-34 years\n",
       "1    Yes  Alive  19.3    18-34 years\n",
       "2     No   Dead  57.5    55-64 years\n",
       "3     No  Alive  47.1    35-54 years\n",
       "4    Yes  Alive  81.4  Over 65 years"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def bin_age(age):\n",
    "    if age < 18:\n",
    "        return None\n",
    "    if age < 35:\n",
    "        return \"18-34 years\"\n",
    "    elif age < 55:\n",
    "        return \"35-54 years\"\n",
    "    elif age < 65:\n",
    "        return \"55-64 years\"\n",
    "    else:\n",
    "        return \"Over 65 years\"\n",
    "\n",
    "data['Age group'] = data['Age'].apply(bin_age)\n",
    "data.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Again, let's check that no missing are present, to ensure that no women under 18 was polled."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of rows with missing values: 0\n"
     ]
    }
   ],
   "source": [
    "print(\"Number of rows with missing values:\", data.isnull().any(axis=1).sum())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's visualize the number on women alive and dead after twenty years, according to their smoking habits and age. Different colors correspond to different couples of smoking habits and status."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3EAAAKZCAYAAAAWFKIXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3XmYFdWZ+PHvC4gQQUFtFUUFd5RFtAUyLoiiomJQRwQmLkQTE7fEJKODSyIa/WmcxKjJZNQYhTiK4DaJiInJiEmIiqLgCnEBBZQgEFFQkYDn90cV7aXp5Ta0tpX+fp6nn7636tSpt07VrXvfe07VjZQSkiRJkqRiaNHUAUiSJEmSymcSJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEid9SiJiTERc0UTrjoi4LSLeiYgn12P51yNiYP74ooi4pfGjXGedIyNiyqdUd9X2NDcRcUVELI6IvzV1LPr0RESKiF3yxzdGxPeaOqZyRMToiPif9Vz2Mzk3NVRd55uIeDEiDv4U1rne58+IODgi5tcx/3N5PDXn87oEJnFqRvIT/sKI2KRk2lcj4tEmDOvTcgBwGNA5pdRnQypKKf2/lNJXGycsiIiXI2K3xqpPtYuI7YHvAnumlLaJiC75h/1WjbiO3SLi1xGxKCL+HhG/i4jdq5X5dkT8LSLejYhbI2Ljxlq/1pVS+kZK6QeNXe+GHj/1JQsN1djnps9CSmmvlNKjTR1HQ5QeT429DyWtP5M4NTetgG81dRANFREtG7jIjsDrKaX3P4141ldE7Ay0SCm93NSxNBM7AktSSm83RmW1fHjvAPwG2B3YGngS+HXJMkcAo4BDgS7ATsBljRHP+mrMJPazVuTYP0u2k6R/diZxam7+E/j3iOhQfUZN3zJHxKMR8dX88ciI+EtE/CQilkbE7Ij4l3z6vIh4OyJOrVbtlhHx+4hYFhF/jIgdS+reI5/394j4a0ScWDJvTET8d0RMioj3gQE1xLttRPwmX/7ViPhaPv104BbgixGxPCLW+cAcETtHxCMRsSQfandHTW2Sl60a7hQRv42Ic6rNfzYijq9vm3JHA5Pyslvk8b+XD/ncuVq9dbXP0RExPV92XkSMrrbsyRHxRr59F9e0XXm5rvm+bJE/vyUi3i6Z/z8RcV5d7V3SRnfn5ZdFxPN5D9WF+XExLyIOLym/WUT8MiIWRMSbkQ15bJnPGxkRUyLiR5ENh50TEUfWsQ2jIuK1fL0vRcRx+fSBwO+BbfPjYAzwp3yxpfm0L+ZlT4uImfn6flftOE0RcXZEvAK8Un39KaUnU0q/TCn9PaX0D+AnwO4RsUVe5FTglymlF1NK7wA/AEbWsi3/FRE/rjbtgWr74N7Iev3mRMQ3S8r1iYjH8/25ICJ+FhGta9uOyPwk3z/vRsRzEdG9tnauFtP1+T59LyKejogDS+a1jGyY35p98nRkPaJExF4lx/TCiLgon75xRFwXEW/lf9dF3lsZec9HRPxHZENib8unn59v51sRcVq1+KqGcpcs/918WxdExFdKym6Rt/F7EfFUfizWNixvneMnIlpExCWRvd7ejohfRcRmNbTZJsBDfHI8Lo+IbfPZrfPllkU23LCyZLm69nnpuWnN+fv0iJgLPFJDDFtGxMT8GPl7RPw5Pnntv5636XMR8X5kr8+tI+KhPK4/RETHkrq+lMe6NLL3iW41NVhk57E5ETG8ZD1rhqqPjogJdWz7PpGd55ZFdn4ZH/UM0Y9azhsR8ZXIXuPLInvv+noNy14U2fvB6xHx5ZLpY/LjosZ9GNlrb1p+DC2MiGtria1j3v6L8hgnRkTnkvmPRsQPInufXRYRD0fEliXzyzqv52Xre484paSu71XbLy3ik/PqknwfbV7X+qQmkVLyz79m8Qe8DgwE7gOuyKd9FXg0f9wFSECrkmUeBb6aPx4JrAK+ArQErgDmAv8FbAwcDiwD2uXlx+TPD8rnXw9MyedtAszL62oF7AMsBvYqWfZdYH+yL1va1LA9fwR+DrQB9gYWAYeWxDqljrbYhWy45cZABdmHs+uqt1X+eDTwP/njU4C/lJTbE1ia11PnNuXlfwsckT++C5iQL9cdeLMB7XMw0CNvm57AQuDYkpiWl7T7tfl+G1hLW8wF9s0f/xWYDXQrmde7jPYeDawAjsjj/RUwB7gY2Aj4GjCnZJ3/C9yUb+dWZL1XXy/Zd//Il2kJnAm8BUQt8Q8Fts3bYhjwPtCppJ3ml5TtwrrH+LHAq0C3PPZLgMdK5ieyZHBzoG0Zr7NjgQUlz58FhpU83zKvc4salu2Tb2uLkrIfkPXwtQCeBr4PtCbr0ZvNJ8fTvkC/fBu6ADOB82rbjnxfPU3Wkxj59ncq81xyErBFvq7vAn8jf40C5wPPk/VMBtArL9seWJCXb5M/75svcznwRH4sVACPAT8o2YergB+SHc9tgUFkx3z3/Bi6M9++XUrOH1dUW/5ysmPxqLxNO5a8Du8CvkD22plHLeeOWo6f0/LjZyegHdn59fZalj+YkuOx2mvnKLLj/SrgiXxefft8NJ+cm9bE9qu8TdY5VvO6b8zbYSPgQPLXFdk57wmyY2074G3gGaB33u6PAJfmZXcje50dltdzQd4Grau91+xDdg4ZXMe5tbZtbw28QTZyZCPgeGDlmv1aw7aNpI7zBtkXaDuTHZP982Ngn2rHyLX5tvbPt2/3Wo6n6vvwceDk/HE7oF8tMW4B/CvZsdYeuBv435L5jwKv5e3bNn9+9Xqe1w+m/veIA/J2/lHedmv2y3n5sdA5X9dNwLhyzg3++fdZ/jV5AP7591n9lbyxdidLkCpoeBL3Ssm8Hnn5rUumLQH2zh+PAe4qmdcOWA1sT/Zh+8/V4ruJTz4kjAF+Vce2bJ/X1b5k2lXAmJJYa03iaqjvWGB69bbKH4/mkw9K7cne3HfMn18J3Jo/rm+bvpC3TxuyDxn/APYoKfv/+CSJq7OuGuK/DvhJ/vj71dp9E7IPP7W92d8OfAfYhiyJuwb4BtCVLEFtUUZ7jwZ+XzLvGLIPCS1L2i2RJQxbAx9R8iETGAFMLtl3r5bM+0K+7DZl7ssZwJD88cHUn8Q9BJxe8rwF2Qe8Nfs4AYeUue7OZMn4iJJprwGDSp5vlNfZpZY6ZgKH5Y/PASblj/sCc6uVvRC4rZZ6zgPuL3m+1nYAhwAvkyV+Lcp9rdSyrneAXvnjv65p/2plRlDyGqs27zXgqJLnR5ANh16zD1dS8kUOcCv5h9v8+W7UncR9WG2fv51v95rX4e4l866gYUnc/wFnlTzfPa+zVQ3Lr3U8lrx2/lDyfE/gw3L2OTUncTvVsZ8uJxvqu0sN814Hvlzy/F7gv0uen0uecADfAyZUe828CRxcUtdlwHxgQA3rGVjGth+U1xkl86dQdxJX9nmD7Iukb5Xsl1XAJiXzJwDfq+V4qr4P/5Rv75YNfN3sDbxT8vxR4JKS52cBv80fN+i8XsO6qr9HjCuZ94XSusjOQYeWzO9ELce0f/415Z/DKdXspJReACaSXafTUAtLHn+Y11d9WruS5/NK1rsc+DtZr8mOQN98KM7SiFgKfJkskVhn2RpsC/w9pbSsZNobZN8g1ysitoqIuyIbyvce8D9kvR51ytf3IDA8nzQcuCN/XN82HUrWw7OCLIFuVW0b3yh5XGddEdE3Iibnw3LeJUu61sS/LWu3+/tkyWNt/kj2weQgsg8jj5J9E92fLJH8mPLau/pxsDiltLrkOWTHxo5kicyCkm27iawXZo2qO0mmlD4oWXYd+bCgGSV1daeMfVliR+D6kuX/TvZtfem21XUsromjAngY+HlKaVzJrOXApiXP1zwubctSY8l6usj/314S57bVjomLyJLiNTdYmRjZDVTeI/tSoHo7lB4XjwA/I+tJXxgRN0fEppQhsqGJMyMbhrkU2KxkXduTJWXV1TYdsuOr9Ph/I5+2xqL8dVNavrbXTk2WpJRWlTz/gOx4qul1WO++rqam2FuR75cyld459QOgTWTD2uvc57WoK/7/JOsxezgfUlj9PaD6a7i2c/ta25yfI+ax9mvmG2Tnu8l1xAO1b/u2wJsppVQyv759U+t5IyKOjIgnIhtGupSs96/09fFOWvsa6urHYF1OJ/siYVZkQ3IH11QoIr4QETflwxjfIzvfdoi1r/mu3h6lbV72eb2B7xEfVKtrR+D+kmNuJtmXeA05pqVPnUmcmqtLyYadlL7prnkD+0LJtNKkan1sv+ZBRLQjG8r1FtkbyB9TSh1K/tqllM4sWTZRu7eAzSOifcm0Hci+uS3HVXn9PVNKm5J9WI4ylx0HjIjseqq2wJoPKfVt01FkCSBkQxFXUdI+efyUWdedZDfT2D6ltBnZEKk18S9g7Xb/Atkwntr8kWxY1cH54ylkw1j7589hw9u71DyynrgtS7Zt05TSXg2tKLJr135B1mO1RUqpA/ACte/Lmo6peWRDOUvbum1K6bF6liuNoyNZAveblNKV1Wa/SDakcI1ewMKUUm0fwP4HGBIRvciGOP5vSZxzqsXZPqV0VD7/v4FZwK75MX0R67bDWtuRUrohpbQvsBfZh9Dz69rOfFsPBP4DOJFsSGIHsp79NeuaR7XrO+uZDtnxtWPJ8x3yaTXGTbVjnLVfOw2x5nXYuWTa9rWUrSkOqDn2VaydANW1fF3q2+flxpjNSGlZSum7KaWdyHrLvxMRhzYwJqi2zRERZO1Wej74BrBDRPxkPeqHbB9vl9e9Rl37plaRXV95L9mwwa3zY3YSa78+OkbJnZtZ9xhcY532TSm9klIaQfZF1A+Be6rVtcZ3yXpq++av0YPWhFjGZjT0vF7fe0TptXhtq9U1Dziy2nHXJqW0Pud76VNjEqdmKaX0KjAe+GbJtEVkb8InRXZzgtOo/UNXuY6KiAMiu8HCD4CpKaV5ZD2Bu+UXam+U/+0XtVwcX0P888ium7kqItpERE+yb0PvqHvJKu3JekiWRsR2lPHhtcQksg8wlwPj82+hKWObjsyXJe+hug8YnX87uyfZDTAos672ZD1jKyKiD/BvJcveAwwuaffLqeNcl1J6hexb9pOAP6WU3iP7APqv5ElcI7R36foWkCU8P46ITfOL6HeOiP4NrYtsSFEi+zBOZDesqOvmHIuAj8muLVrjRuDCiNgrr2OziBhabgB579XvyK6VrKl3+1fA6RGxZ57sXUI2PKtGKaX5wFNkPXD3ppTW9GI+CbwX2Q0+2uav0e4RsV8+vz3wHrA8IvYguyaorrj3y7+t34jsC5wVZN+2r7m5zOu1LNqeLElZBLSKiO+zdk/jLcAPImLXyPSM7CYvE4FtIuK8yG5k0j4i+ubLjAMuiYiKyG7k8H2yZLY2E4CReZt+gexLqQar4XW4B9l1r7Wp6fgZB3w7spsEtSPrAR1fredvjYXAFlHDjU9qUd8+b5CIGBwRu+SJ0Xtk+3t1PYvVZAJwdEQcmh8/3yX7Yqb0i49lZNcuHhQRV6/HOh7PYzsnIlpFxBCya0bXR2uya7sWAasiu+HJ4TWUuywiWudfVAwmu2atunX2YUScFBEV+XvB0nxyTe3anuxcuzSyG4U05Lht0Hmd+t8jjonsxmStyYaCliaSNwJX5l+Skb8uhzQgVukzYRKn5uxysg/Bpb5GltAsIft2/rHqCzXQnWRvVH8nu/HCl6FqWOLhZMMR3yIbQrLmxgXlGkF2HchbwP1k14v9vsxlLyO76P5dst6x+8pdaUrpo7z8QLLtWzO91m2K7K5/y1NKc0uqOodsqMzfyD7U31ZOXXmRs4DLI2IZ2QfeCSXLvgicnce2gOx6pfp+1+iPZEPO5pY8D2B6SZkNae/qTiH7YPVSHt89ZNddNEhK6SXgx2Qf+BaSXaf5lzrKf0B2HeNf8qFC/VJK95O17V2RDXF6gSzhLtdxwH7AV+KTO9Ytj4gd8nX+luw6w8lkQ7TeoP4Pb2PzbVkzlHJNwnEM2XU0c8hudHML2VBGgH8n+6C2jKx3cnw969g0L/dOHtMSsp4KyL7xr60df0d2HeHL+XIrWHuY27Vkx+PDZInCL8muf1xGdiOMY8iO51f45K6zVwDTgOfIboryTD6tRimlh8iu8XmEbHjgOndibIBzyNrwb2TtPY4sIalpvescP2TX591ONjRuDll7nFvL8rPy+mfny9c5XK+Mfd5QuwJ/IPsC63Gyob+PNrSSlNJfyb70+Wke0zHAMSmlldXKLSXb50dGRIN+ty+v63iyL4uW5uubSC37pp66lpF9YTmB7Hj/N7JeqlJ/y+e9Rfbl1Dfy/VW9rpr24SDgxYhYTnYDr+HVhv+ucR3Z6I3FZDcO+W0DtqGh5/X63iPOJbuhzwKyc8bbfNK215O1z8P58k+QXZ8pfa6suWuRJH1qIuICsuGDFzR1LPr8i4iDyHqiupT09H6W63+Y7KYPMz/rdTe1iPgh2c0wTq23sD5TETEVuDGldFu9hVW2vAd5KdlQ7DlNHY9ULnviJH0WXqekp02qTT487VvALU2RwAGklA5vLglcZL9j1jMf+tmHrOfn/qaOSxAR/SNim3w45alkt8ovu/dKtYuIY/IhxJuQ9cA/T/Y+JRVGq/qLSNKGSSlNqL+Umrv8msdpZL8t95V6iqtxtCcbHrct2ZCyH5Pdhl9Nb3eyYYDtyO5sekJ+Ta023BCyYcBBds4ZnhyapoJxOKUkSZIkFYjDKSVJkiSpQEziJEmSJKlAPhfXxG255ZapS5cuTR2GJEmSJDWJp59+enFKqaKcsp+LJK5Lly5MmzatqcOQJEmSpCYREW+UW9bhlJIkSZJUICZxkiRJklQgJnGSJEmSVCCfi2viJEmSJDWNf/zjH8yfP58VK1Y0dSjNQps2bejcuTMbbbTRetdhEidJkiQ1Y/Pnz6d9+/Z06dKFiGjqcP6ppZRYsmQJ8+fPp2vXrutdj8MpJUmSpGZsxYoVbLHFFiZwn4GIYIstttjgXk+TOEmSJKmZM4H77DRGW5vESZIkSSrblVdeyV577UXPnj3Ze++9mTp16gbV9+ijjzJ48OBGiq72+K677jo++OCDepcvt1xT8po4SZIkSWV5/PHHmThxIs888wwbb7wxixcvZuXKlU0Wz6pVq2jV6pOUpq74rrvuOk466SS+8IUv1FlnueWakj1xkiRJksqyYMECttxySzbeeGMAttxyS7bddlu6dOnCRRddxBe/+EUqKyt55plnOOKII9h555258cYbgeymHueffz7du3enR48ejB8/fp36n3rqKXr37s3s2bN5//33Oe2009hvv/3o3bs3v/71rwEYM2YMQ4cO5ZhjjuHwww8vK74bbriBt956iwEDBjBgwAAAzjzzTCorK9lrr7249NJLAWos165du6r677nnHkaOHAnA3XffTffu3enVqxcHHXRQYzVxeVJKTf637777JkmSJEmfvZdeeqnsssuWLUu9evVKu+66azrzzDPTo48+mlJKaccdd0w///nPU0opnXfeealHjx7pvffeS2+//XaqqKhIKaV0zz33pIEDB6ZVq1alv/3tb2n77bdPb731Vpo8eXI6+uij01/+8pe0zz77pDfeeCOllNKFF16Ybr/99pRSSu+8807adddd0/Lly9Ntt92Wtttuu7RkyZKy41sT46JFi6qer1l+1apVqX///unZZ5+tsdwmm2xS9fjuu+9Op556akoppe7du6f58+dXxdcQNbU5MC2VmT/ZEydJkiSpLO3atePpp5/m5ptvpqKigmHDhjFmzBgAvvSlLwHQo0cP+vbtS/v27amoqKBNmzYsXbqUKVOmMGLECFq2bMnWW29N//79eeqppwCYOXMmZ5xxBg888AA77LADAA8//DBXX301e++9NwcffDArVqxg7ty5ABx22GFsvvnmDYqvugkTJrDPPvvQu3dvXnzxRV566aUGtcX+++/PyJEj+cUvfsHq1asbtOyG8po4SZIkSWVr2bIlBx98MAcffDA9evRg7NixAFVDGFu0aFH1eM3zVatWkXU21axTp06sWLGC6dOns+222wLZiMF7772X3Xfffa2yU6dOZZNNNmlQfGuGQK4xZ84cfvSjH/HUU0/RsWNHRo4cWett/0vvJlla5sYbb2Tq1Kk8+OCD7L333syYMYMtttii1rgakz1xkiRJksry17/+lVdeeaXq+YwZM9hxxx3LWvaggw5i/PjxrF69mkWLFvGnP/2JPn36ANChQwcefPBBLrroIh599FEAjjjiCH76059WJX/Tp0+vte499tij3vjat2/PsmXLAHjvvffYZJNN2GyzzVi4cCEPPfRQ1TKl5QC23nprZs6cyccff8z9999fNf21116jb9++XH755Wy55ZbMmzevrHZoDPbESZIkSSrL8uXLOffcc1m6dCmtWrVil1124eabb2bixIn1Lnvcccfx+OOP06tXLyKCa665hm222YZZs2YBWbL0wAMPcOSRR3Lrrbfyve99j/POO4+ePXuSUqJLly41rmfx4sVViV5t8QGcccYZHHnkkXTq1InJkyfTu3dv9tprL3baaSf233//qvqql7v66qsZPHgw22+/Pd27d2f58uUAnH/++bzyyiuklDj00EPp1avXBrdvuaKubs3PSmVlZZo2bVpThyFJkiQ1OzNnzqRbt25NHcZ6mzhxIrNnz+ab3/xmU4dStpraPCKeTilVlrO8PXGSJEmSCqsxfyi8KLwmTpIkSZIKxCROkiRJkgrEJE6SJEmSCsQkTpIkSZIKxCROkiRJkgrEJE6SJEmSCsSfGJAkSZJUpcuoBxu1vtevPrreMh9++CGDBg3ikUceYd68eXTt2pUbbriBc889F4BzzjmHyspKRo4cud5xDBkyhLfffpvHH3+8atro0aNp164d//7v/873v/99DjroIAYOHNjguocPH84PfvADdt111/WOryFM4iR9Ls3co7g/OtpcdJs1s6lDkCT9k7j11ls5/vjjadmyJQBbbbUV119/PV//+tdp3br1Bte/dOlSnnnmGdq1a8ecOXPo2rXrOmUuv/zy9a7/zDPP5JprruEXv/jFhoRZNodTSpIkSWpSd9xxB0OGDKl6XlFRwaGHHsrYsWPXKTtjxgz69etHz549Oe6443jnnXfqrf/ee+/lmGOOYfjw4dx11101lhk5ciT33HMPDz30ECeeeGLV9EcffZRjjjkGgIcffpgvfvGL7LPPPgwdOpTly5cDcOCBB/KHP/yBVatWNWi715dJnCRJkqQms3LlSmbPnk2XLl3Wmj5q1Ch+/OMfs3r16rWmn3LKKfzwhz/kueeeo0ePHlx22WX1rmPcuHGMGDGCESNGMG7cuDrLHnbYYTzxxBO8//77AIwfP55hw4axePFirrjiCv7whz/wzDPPUFlZybXXXgtAixYt2GWXXXj22WcbsOXrzyROkiRJUpNZvHgxHTp0WGd6165d6dOnD3feeWfVtHfffZelS5fSv39/AE499VT+9Kc/1Vn/woULefXVVznggAPYbbfdaNWqFS+88EKt5Vu1asWgQYN44IEHWLVqFQ8++CBDhgzhiSee4KWXXmL//fdn7733ZuzYsbzxxhtVy2211Va89dZbDd389eI1cZIkSZKaTNu2bVmxYkWN8y666CJOOOEEDjrooPWuf/z48bzzzjtV18G999573HXXXVxxxRW1LjNs2DD+67/+i80335z99tuP9u3bk1LisMMOq7Unb8WKFbRt23a942wIk7hPiTdl+HzzhgySJEmfDx07dmT16tWsWLGCNm3arDVvjz32YM8992TixIn06dOHzTbbjI4dO/LnP/+ZAw88kNtvv72qV+5nP/sZkN3JstS4ceP47W9/yxe/+EUA5syZw2GHHVZnEnfwwQdz+umn84tf/IJhw4YB0K9fP84++2xeffVVdtllFz744APmz5/PbrvtBsDLL7/MXnvt1TiNUg+TOEmSJElVyvlJgMZ2+OGHM2XKlBpv73/xxRfTu3fvqudjx47lG9/4Bh988AE77bQTt912GwCzZs1i//33X2vZ119/nblz59KvX7+qaV27dmXTTTdl6tSptcbTsmVLBg8ezJgxY6purlJRUcGYMWMYMWIEH330EQBXXHEFu+22GwsXLqRt27Z06tRp/RuhASKl9JmsqC6VlZVp2rRpTR1Go7In7vPNnrjPP19Dn3++jiTpn8PMmTPp1q1p33enT5/Otddey+23377edQwePJj77ruvUX6SoKF+8pOfsOmmm3L66aeXVb6mNo+Ip1NKleUsb0+cJEmSpCbVu3dvBgwYwOrVq6t+K66hJk6c2MhRla9Dhw6cfPLJn9n6TOIkSZIkNbnTTjutqUNYb1/5ylc+0/X5EwOSJEmSVCAmcZIkSZJUICZxkiRJklQgXhP3KTnxQpv28+z5pg5AkiRJWk9mGpIkSZI+MXqzRq7v3XqLfPjhhwwaNIhHHnmEefPm0bVrV2644QbOPfdcIPsB78rKSkaOHNng1Y8ZM4bzzz+fzp07s3z5cnbaaScuvfRS/uVf/qXBddWkS5cuTJs2jU033ZSBAwfyyCOP0KrVp5tmOZxSkiRJUpO69dZbOf7446t+XmCrrbbi+uuvZ+XKlY1S/7Bhw5g+fTqvvPIKo0aN4vjjj2fmzMb9vdPWrVtz6KGHMn78+EattyYmcZIkSZKa1B133MGQIUOqnldUVHDooYcyduzYdcrOmDGDfv360bNnT4477jjeeeedBq1rwIABnHHGGdx8880AvPbaawwaNIh9992XAw88kFmzZgHwwAMP0LdvX3r37s3AgQNZuHAhAEuWLOHwww+nd+/efP3rXyelVFX3scceyx133NHg7W8okzhJkiRJTWblypXMnj2bLl26rDV91KhR/PjHP2b16tVrTT/llFP44Q9/yHPPPUePHj247LLLGrzOffbZpypZO+OMM/jpT3/K008/zY9+9CPOOussAA444ACeeOIJpk+fzvDhw7nmmmsAuOyyyzjggAOYPn06X/rSl5g7d25Vvd27d+epp55qcDwN5TVxkiRJkprM4sWL6dChwzrTu3btSp8+fbjzzjuVIQQOAAAgAElEQVSrpr377rssXbqU/v37A3DqqacydOjQBq9zTe/Z8uXLeeyxx9aq46OPPgJg/vz5DBs2jAULFrBy5Uq6du0KwJ/+9Cfuu+8+AI4++mg6duxYtWzLli1p3bo1y5Yto3379g2Oq1z2xEmSJElqMm3btmXFihU1zrvooov44Q9/yMcff9yo65w+fTrdunXj448/pkOHDsyYMaPqb821cueeey7nnHMOzz//PDfddNNaMUZErXV/9NFHtGnTplHjrc4kTpIkSVKT6dixI6tXr64xkdtjjz3Yc889mThxIgCbbbYZHTt25M9//jMAt99+e1Wv3M9+9jN+9rOf1bu+P/7xj9x888187WtfY9NNN6Vr167cfffdQNZD9+yzzwJZr992220HsNa1eQcddFDVdW8PPfTQWtfkLVmyhIqKCjbaaKMGt0NDlD2cMiJaAtOAN1NKgyNic2A80AV4HTgxpfROXvZC4HRgNfDNlNLvGjluSZIkSZ+GMn4SoLEdfvjhTJkyhYEDB64z7+KLL6Z3795Vz8eOHcs3vvENPvjgA3baaSduu+02AGbNmsX+++9fY/3jx49nypQpfPDBB3Tt2pV7772Xbt26AdlNVc4880yuuOIK/vGPfzB8+HB69erF6NGjGTp0KNtttx39+vVjzpw5AFx66aWMGDGCffbZh/79+7PDDjtUrWfy5MkcddRRjdYutYnSu6nUWTDiO0AlsGmexF0D/D2ldHVEjAI6ppT+IyL2BMYBfYBtgT8Au6WUVtdWd2VlZZo2bdqGbsvnSo+xPZo6BNXh+VP9ue/Pu5l7dGvqEFSPbrMa99bMkqSmMXPmzKqEpqlMnz6da6+9lttvv3296xg8eDD33XcfrVu3bsTIGub444/nqquuYvfdd6+zXE1tHhFPp5Qqy1lPWcMpI6IzcDRwS8nkIcCafsWxwLEl0+9KKX2UUpoDvEqW0EmSJEnSOnr37s2AAQPWuRNlQ0ycOLFJE7iVK1dy7LHH1pvANYZyr4m7DrgAKL2icOuU0gKA/P9W+fTtgHkl5ebn09YSEWdExLSImLZo0aIGBy5JkiTpn8dpp51W9WPfRdS6dWtOOeWUz2Rd9SZxETEYeDul9HSZddZ0q5Z1xmymlG5OKVWmlCorKirKrFqSJEmSmrdybmyyP/CliDgKaANsGhH/AyyMiE4ppQUR0Ql4Oy8/H9i+ZPnOwFuNGbQkSZIkNVf19sSllC5MKXVOKXUBhgOPpJROAn4DnJoXOxX4df74N8DwiNg4IroCuwJPNnrkkiRJktQMlf0TAzW4GpgQEacDc4GhACmlFyNiAvASsAo4u647U0qSJEmSytegJC6l9CjwaP54CXBoLeWuBK7cwNgkSZIkfcYa+6eyyvlppw8//JBBgwYxadIk9t13X+6++2569MjiuOaaa5g9ezY33nhjo8bVvXt3evfuvdbPGpx00kmccMIJHHvssXzlK19h1KhR63W3yUMPPZT77ruPzTbbrDFDrlLu3SklSZIk6VNx6623cvzxx7PJJptw3XXXcdZZZ5FS4s033+Smm27iqquuatT1Pffcc7Rq1YpHHnmEDz/8sMYyt91223r/XMC//du/NXrSWcokTpIkSVKTuuOOOxgyZAgAgwYNolOnTvzqV7/i29/+NqNHj6Zjx44AXH311fTp04eePXty+eWXA7Bs2TKOPPJIevXqRffu3bnnnnvqXd+4ceM45ZRTOOSQQ5g4cWKNZQ444ABmzJjBT3/6Uy666KKq6bfccgvf/va3ARg7dix9+vRh77335qyzzuLjj7NfZBsyZAh33nnn+jdIPUziJEmSJDWZlStXMnv2bLp06VI17brrruPiiy9m0aJFnHzyyQBMmjSJuXPnMnXqVGbMmMFjjz3GY489xqRJk+jSpQvPPvssL7zwAocddli965wwYQLDhg1jxIgRjBs3rs6yQ4cOXSsxHD9+PMOGDeOFF17g/vvv57HHHmPGjBmsWrWKu+66C4Att9ySZcuWsXTp0vVokfptyI1NJEmSJGmDLF68mA4dOqw1bdttt+WQQw5h8ODBVdMefvhhHnroIXr37g3A8uXLefnll+nbty+jRo1i1KhRHHPMMey///51ru/xxx+nc+fObLfddmy11VZ87Wtf49133631+rVtttmGzp07M23aNHbYYQfmzJlD3759uf7663nqqaeorKwEsuv6tt/+k19aq6ioYMGCBetsW2MwiZMkSZLUZNq2bcuKFSvWmd6iRQtatPhk4GBKiUsuuYTTTz99nbLTpk1j0qRJnH/++QwePHit4Y/VjRs3jhdeeKGq5++9997j/vvvZ+TIkbUuM2zYMCZMmECXLl3413/9VyKClBKnnXYaP/jBD2pcZsWKFbRt27bWOjeEwyklSZIkNZmOHTuyevXqGhO5UkcccQS//OUvef/99wGYP38+ixcv5s0336Rdu3acfPLJfOc73+GZZ54B4IILLuCBBx5Yq47Vq1dz77338tJLL/H666/z+uuvc99999U7pPKEE07gvvvu46677mLYsGEADBw4kAkTJrB48WIAlixZwty5cwH4+OOPWbx48Vo9c43JnjhJkiRJVcr5SYDGdvjhhzNlyhQGDhxYa5mjjjqKWbNm0a9fPwDat2/PnXfeyUsvvcSoUaNo0aIFrVu3rror5HPPPcfQoUPXqmPy5Ml07dqVrbfeumragAEDOOmkk1i4cGGt695iiy3YZZddeO2119hnn30A6NGjB5deeikDBw7k448/ZqONNuLGG29khx124Mknn+SAAw6gZcuW690mdYmU0qdScUNUVlamadOmNXUYjaqxf19DjaspTk5qmJl7dGvqEFSPbrNmNnUIkqRGMHPmTLp1a9r33enTp3Pttdeu9ZttGyKlxKBBg/jd737XKPU11Nlnn82JJ55I//79a5xfU5tHxNMppcpy6nc4pSRJkqQm1bt3bwYMGMDq1asbpb6IaLIEDrLtqS2BawwOp5QkSZLU5E477bSmDqHRfPWrX/1U67cnTpIkSZIKxCROkiRJkgrEJE6SJEmSCsQkTpIkSZIKxBubSJIkSarS2D/zU85P0nz44YcMGjSISZMmse+++3L33XfTo0f2k13XXHMNs2fPrvr9tw11ySWXcNttt1FRUcH7779Pz549ufLKK9ljjz02uO5Vq1ax5ZZbsnTpUv72t79x+umn8+CDDzZC1GuzJ06SJElSk7r11ls5/vjj2WSTTbjuuus466yzSCnx5ptvctNNN3HVVVc16vrOP/98ZsyYwSuvvMIJJ5zAgAEDWLJkSaOuY5tttmHzzTdn6tSpjVovmMRJkiRJamJ33HEHQ4YMAWDQoEF06tSJX/3qV3z7299m9OjRdOzYEYCrr76aPn360LNnTy6//HIAli1bxpFHHkmvXr3o3r0799xzT4PWPWLECAYMGMBdd90FwFNPPUX//v3Zd999OfLII1m4cCEAN954I/vttx+9evVi6NChfPjhhwC89tpr9O3bl/3224/Ro0evVfexxx7LHXfcsd7tUhuTOEmSJElNZuXKlcyePZsuXbpUTbvuuuu4+OKLWbRoESeffDIAkyZNYu7cuUydOpUZM2bw2GOP8dhjjzFp0iS6dOnCs88+ywsvvMBhhx3W4Bj22WcfZs2axUcffcS3vvUt7r33Xp5++mlOOukkvve97wEwdOhQnnrqKZ599ll23nlnxowZA8C5557Lt771LZ566ikqKirWqreyspI///nP69cwdfCaOEmfSyde6Onp8+75pg5AkvRPYfHixXTo0GGtadtuuy2HHHIIgwcPrpr28MMP89BDD9G7d28Ali9fzssvv0zfvn0ZNWoUo0aN4phjjmH//fdvcAwpJQBmzpzJiy++yMCBAwFYvXo1nTt3BuC5557j+9//PkuXLmXZsmVVsT3++OM88MADAJx88slceumlVfVutdVWvPXWWw2Opz5+SpIkSZLUZNq2bcuKFSvWmd6iRQtatPhk4GBKiUsuuYTTTz99nbLTpk1j0qRJnH/++QwePJiLLrqoQTFMnz6dAw44gJQSPXv2rLH37JRTTuGhhx6ie/fu3HLLLTzxxBMARAQRUWO9K1asoG3btg2KpRwOp5QkSZLUZDp27Mjq1atrTORKHXHEEfzyl7/k/fffB2D+/PksXryYN998k3bt2nHyySfzne98h2eeeQaACy64oKqHrC4TJkxg8uTJDBs2jD333JM333yTJ598EsiGer744osAvP/++2yzzTb84x//4M4776xavl+/fkyYMAFgnevfXn75Zbp3715mS5TPnjhJkiRJVcr5SYDGdvjhhzNlypSqYYw1Oeqoo5g1axb9+vUDoH379tx555289NJLjBo1ihYtWtC6deuqnyJ47rnnGDp0aI11/ed//idjxozh/fffp0ePHkyePJktttgCgHvuuYdvfvObLFu2jFWrVvHd736Xvfbai8svv5w+ffqwww470L1796qk84YbbuDLX/4y1157Lccdd9xa65k8eTJHH330BrdPdbFm/GdTqqysTNOmTWvqMBpVj7E9mjoE1eH5U72a5/PO19Dnn68jSfrnMHPmTLp1a9zfhmuo6dOnc+2113L77bc3Sn0pJQYNGsTvfve7RqlvfWM48MADefDBB9lss83WmldTm0fE0ymlynLqdjilJEmSpCbVu3dvBgwYwOrVqxulvoho0gQO4O233+aCCy5YJ4FrDA6nlCRJktTkTjvttKYOoVFtvfXWfOlLX/pU6rYnTpIkSWrmPg+XWDUXjdHWJnGSJElSM9amTRuWLFliIvcZSCmxZMkS2rRps0H1OJxSkiRJasY6d+7M/PnzWbRoUVOH0iy0adOm6gfE15dJnCRJktSMbbTRRnTt2rWpw1ADOJxSkiRJkgrEJE6SJEmSCsQkTpIkSZIKxCROkiRJkgrEJE6SJEmSCsQkTpIkSZIKxCROkiRJkgrEJE6SJEmSCsQkTpIkSZIKxCROkiRJkgrEJE6SJEmSCsQkTpIkSZIKxCROkiRJkgrEJE6SJEmSCsQkTpIkSZIKxCROkiRJkgrEJE6SJEmSCsQkTpIkSZIKxCROkiRJkgrEJE6SJEmSCsQkTpIkSZIKxCROkiRJkgrEJE6SJEmSCqTeJC4i2kTEkxHxbES8GBGX5dNHR8SbETEj/zuqZJkLI+LViPhrRBzxaW6AJEmSJDUnrcoo8xFwSEppeURsBEyJiIfyeT9JKf2otHBE7AkMB/YCtgX+EBG7pZRWN2bgkiRJktQc1dsTlzLL86cb5X+pjkWGAHellD5KKc0BXgX6bHCkkiRJkqTyromLiJYRMQN4G/h9SmlqPuuciHguIm6NiI75tO2AeSWLz8+nSZIkSZI2UFlJXEppdUppb6Az0CciugP/DewM7A0sAH6cF4+aqqg+ISLOiIhpETFt0aJF6xW8JEmSJDU3Dbo7ZUppKfAoMCiltDBP7j4GfsEnQybnA9uXLNYZeKuGum5OKVWmlCorKirWK3hJkiRJam7KuTtlRUR0yB+3BQYCsyKiU0mx44AX8se/AYZHxMYR0RXYFXiyccOWJEmSpOapnLtTdgLGRkRLsqRvQkppYkTcHhF7kw2VfB34OkBK6cWImAC8BKwCzvbOlJIkSZLUOOpN4lJKzwG9a5h+ch3LXAlcuWGhSZIkSZKqa9A1cZIkSZKkpmUSJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFUm8SFxFtIuLJiHg2Il6MiMvy6ZtHxO8j4pX8f8eSZS6MiFcj4q8RccSnuQGSJEmS1JyU0xP3EXBISqkXsDcwKCL6AaOA/0sp7Qr8X/6ciNgTGA7sBQwCfh4RLT+N4CVJkiSpuak3iUuZ5fnTjfK/BAwBxubTxwLH5o+HAHellD5KKc0BXgX6NGrUkiRJktRMlXVNXES0jIgZwNvA71NKU4GtU0oLAPL/W+XFtwPmlSw+P58mSZIkSdpAZSVxKaXVKaW9gc5An4joXkfxqKmKdQpFnBER0yJi2qJFi8qLVpIkSZKauQbdnTKltBR4lOxat4UR0Qkg//92Xmw+sH3JYp2Bt2qo6+aUUmVKqbKiomI9QpckSZKk5qecu1NWRESH/HFbYCAwC/gNcGpe7FTg1/nj3wDDI2LjiOgK7Ao82diBS5IkSVJz1KqMMp2AsfkdJlsAE1JKEyPicWBCRJwOzAWGAqSUXoyICcBLwCrg7JTS6k8nfEmSJElqXupN4lJKzwG9a5i+BDi0lmWuBK7c4OgkSZIkSWtp0DVxkiRJkqSmZRInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQVSbxIXEdtHxOSImBkRL0bEt/LpoyPizYiYkf8dVbLMhRHxakT8NSKO+DQ3QJIkSZKak1ZllFkFfDel9ExEtAeejojf5/N+klL6UWnhiNgTGA7sBWwL/CEidksprW7MwCVJkiSpOaq3Jy6ltCCl9Ez+eBkwE9iujkWGAHellD5KKc0BXgX6NEawkiRJktTcNeiauIjoAvQGpuaTzomI5yLi1ojomE/bDphXsth8akj6IuKMiJgWEdMWLVrU4MAlSZIkqTkqO4mLiHbAvcB5KaX3gP8Gdgb2BhYAP15TtIbF0zoTUro5pVSZUqqsqKhocOCSJEmS1ByVlcRFxEZkCdwdKaX7AFJKC1NKq1NKHwO/4JMhk/OB7UsW7wy81XghS5IkSVLzVc7dKQP4JTAzpXRtyfROJcWOA17IH/8GGB4RG0dEV2BX4MnGC1mSJEmSmq9y7k65P3Ay8HxEzMinXQSMiIi9yYZKvg58HSCl9GJETABeIruz5dnemVKSJEmSGke9SVxKaQo1X+c2qY5lrgSu3IC4JEmSJEk1aNDdKSVJkiRJTcskTpIkSZIKxCROkiRJkgrEJE6SJEmSCqScu1NKkiRJzcrMPbo1dQiqR7dZM5s6hCZjT5wkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYi/EydJ0j8hf+Pq8685/8aVpA1jT5wkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFUirpg7gn9Xzc+Y2dQiSJEmS/gnZEydJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQXiTwx8SrqsuLOpQ1AdXm/qACRJkqT1ZE+cJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQVSL1JXERsHxGTI2JmRLwYEd/Kp28eEb+PiFfy/x1LlrkwIl6NiL9GxBGf5gZIkiRJUnNSTk/cKuC7KaVuQD/g7IjYExgF/F9KaVfg//Ln5POGA3sBg4CfR0TLTyN4SZIkSWpu6k3iUkoLUkrP5I+XATOB7YAhwNi82Fjg2PzxEOCulNJHKaU5wKtAn8YOXJIkSZKaowZdExcRXYDewFRg65TSAsgSPWCrvNh2wLySxebn0yRJkiRJG6hVuQUjoh1wL3BeSum9iKi1aA3TUg31nQGcAbDDDjuUG4YkSZL0qTvxwrI/JquJPN/UATShsnriImIjsgTujpTSffnkhRHRKZ/fCXg7nz4f2L5k8c7AW9XrTCndnFKqTClVVlRUrG/8kiRJktSslHN3ygB+CcxMKV1bMus3wKn541OBX5dMHx4RG0dEV2BX4MnGC1mSJEmSmq9y+on3B04Gno+IGfm0i4CrgQkRcTowFxgKkFJ6MSImAC+R3dny7JTS6kaPXJIkSZKaoXqTuJTSFGq+zg3g0FqWuRK4cgPikiRJkiTVoEF3p5QkSZIkNS2TOEmSJEkqEJM4SZIkSSoQkzhJkiRJKhB/xVCSJEmq5vk5c5s6BKlW9sRJkiRJUoHYEyfpc8lvQCVJkmpmT5wkSZIkFYhJnCRJkiQViEmcJEmSJBWISZwkSZIkFYhJnCRJkiQViEmcJEmSJBWIPzEgSdI/oRMv9C3+8+75pg5AUmHZEydJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEidJkiRJBWISJ0mSJEkFYhInSZIkSQViEifp/7d359GWnWWZwJ+XFARCIAwZiBCSGIN2EDuQIoDQDRINiAMBGRIbRKZ0t9BLGl0ubOlmkLgio4CiMsfGEEFFkGYxSINIY0zCICTMEAgxQMIsMwlv/3FOya3LrVRR+9zs+qp+v7Vq3X32mZ5krTep5+69vw0AwECUOAAAgIEocQAAAAPZMncAgI0c9c2z547ATnxi7gAAsI9yJA4AAGAgShwAAMBAlDgAAICBKHEAAAADUeIAAAAGosQBAAAMRIkDAAAYiBIHAAAwECUOAABgIEocAADAQJQ4AACAgShxAAAAA1HiAAAABqLEAQAADESJAwAAGIgSBwAAMBAlDgAAYCBKHAAAwECUOAAAgIEocQAAAANR4gAAAAaixAEAAAxEiQMAABjITktcVb24qi6vqgvX7HtiVf1LVb1n+edea5777ar6aFV9qKrusVnBAQAA9kW7ciTupUnuucH+Z3X38cs/r0uSqjouyalJbr18z/Oqar9VhQUAANjX7bTEdffbknxhFz/v3knO6e5vdffFST6a5MQJ+QAAAFhjyjVxj66q9y5Pt7zxct/Nk3xqzWsuXe4DAABgBXa3xP1xkmOSHJ/k00mesdxfG7y2N/qAqjq9qi6oqguuuOKK3YwBAACwb9mtEtfdn+3uq7r7u0lekO+dMnlpkiPWvPQWSS7bwWc8v7u3dvfWQw45ZHdiAAAA7HN2q8RV1eFrHt4nybaVK1+T5NSq2r+qjk5ybJLzpkUEAABgmy07e0FVvTzJ3ZIcXFWXJnlCkrtV1fFZnCr5iST/OUm6+6KqekWS9ye5MsmjuvuqzYkOAACw79lpievu0zbY/aKref0ZSc6YEgoAAICNTVmdEgAAgGuYEgcAADAQJQ4AAGAgShwAAMBAlDgAAICBKHEAAAADUeIAAAAGosQBAAAMRIkDAAAYiBIHAAAwECUOAABgIEocAADAQJQ4AACAgShxAAAAA1HiAAAABqLEAQAADESJAwAAGIgSBwAAMBAlDgAAYCBKHAAAwECUOAAAgIEocQAAAANR4gAAAAaixAEAAAxEiQMAABiIEgcAADAQJQ4AAGAgShwAAMBAlDgAAICBKHEAAAAD2TJ3AABg9d538SVzRwBgkzgSBwAAMBAlDgAAYCBKHAAAwECUOAAAgIEocQAAAANR4gAAAAaixAEAAAxEiQMAABiIEgcAADAQJQ4AAGAgShwAAMBAlDgAAICBKHEAAAADUeIAAAAGosQBAAAMRIkDAAAYiBIHAAAwECUOAABgIEocAADAQJQ4AACAgShxAAAAA1HiAAAABrLTEldVL66qy6vqwjX7blJVb6qqjyx/3njNc79dVR+tqg9V1T02KzgAAMC+aFeOxL00yT3X7Xtckjd397FJ3rx8nKo6LsmpSW69fM/zqmq/laUFAADYx+20xHX325J8Yd3ueyc5a7l9VpJT1uw/p7u/1d0XJ/lokhNXlBUAAGCft7vXxB3W3Z9OkuXPQ5f7b57kU2ted+ly3/epqtOr6t4lCnwAABTcSURBVIKquuCKK67YzRgAAAD7llUvbFIb7OuNXtjdz+/urd299ZBDDllxDAAAgL3T7pa4z1bV4Umy/Hn5cv+lSY5Y87pbJLls9+MBAACw1u6WuNckechy+yFJXr1m/6lVtX9VHZ3k2CTnTYsIAADANlt29oKqenmSuyU5uKouTfKEJGcmeUVVPTzJJUnunyTdfVFVvSLJ+5NcmeRR3X3VJmUHAADY5+y0xHX3aTt46qQdvP6MJGdMCQUAAMDGVr2wCQAAAJtIiQMAABiIEgcAADAQJQ4AAGAgShwAAMBAlDgAAICBKHEAAAADUeIAAAAGosQBAAAMRIkDAAAYiBIHAAAwECUOAABgIEocAADAQJQ4AACAgShxAAAAA1HiAAAABqLEAQAADESJAwAAGMiWuQMAAKt31DfPnjsCO/GJuQMAw3IkDgAAYCBKHAAAwECUOAAAgIEocQAAAANR4gAAAAaixAEAAAxEiQMAABiIEgcAADAQJQ4AAGAgShwAAMBAlDgAAICBKHEAAAADUeIAAAAGosQBAAAMRIkDAAAYiBIHAAAwECUOAABgIEocAADAQJQ4AACAgShxAAAAA1HiAAAABqLEAQAADESJAwAAGIgSBwAAMBAlDgAAYCBKHAAAwECUOAAAgIEocQAAAANR4gAAAAaixAEAAAxEiQMAABiIEgcAADAQJQ4AAGAgW6a8uao+keRfk1yV5Mru3lpVN0nyF0mOSvKJJA/o7i9OiwkAAECymiNxP9Xdx3f31uXjxyV5c3cfm+TNy8cAAACswGacTnnvJGctt89KcsomfAcAAMA+aWqJ6yRvrKp3VtXpy32Hdfenk2T589CJ3wEAAMDSpGvikty5uy+rqkOTvKmqPrirb1yWvtOT5Ja3vOXEGAAAAPuGSUfiuvuy5c/Lk7wqyYlJPltVhyfJ8uflO3jv87t7a3dvPeSQQ6bEAAAA2GfsdomrqutX1Q22bSc5OcmFSV6T5CHLlz0kyaunhgQAAGBhyumUhyV5VVVt+5yzu/v1VXV+kldU1cOTXJLk/tNjAgAAkEwocd398ST/foP9n09y0pRQAAAAbGwzbjEAAADAJlHiAAAABqLEAQAADESJAwAAGIgSBwAAMBAlDgAAYCBKHAAAwECUOAAAgIEocQAAAANR4gAAAAaixAEAAAxEiQMAABiIEgcAADAQJQ4AAGAgShwAAMBAlDgAAICBKHEAAAADUeIAAAAGosQBAAAMRIkDAAAYiBIHAAAwECUOAABgIFvmDgAAAHuao7559twR2IlPzB1gRo7EAQAADESJAwAAGIgSBwAAMBAlDgAAYCBKHAAAwECUOAAAgIEocQAAAANR4gAAAAaixAEAAAxEiQMAABiIEgcAADAQJQ4AAGAgShwAAMBAlDgAAICBKHEAAAADUeIAAAAGosQBAAAMRIkDAAAYiBIHAAAwECUOAABgIEocAADAQJQ4AACAgShxAAAAA1HiAAAABqLEAQAADESJAwAAGIgSBwAAMBAlDgAAYCBKHAAAwECUOAAAgIFsWomrqntW1Yeq6qNV9bjN+h4AAIB9yaaUuKraL8kfJfnZJMclOa2qjtuM7wIAANiXbNaRuBOTfLS7P97d305yTpJ7b9J3AQAA7DOqu1f/oVX3S3LP7n7E8vGDk9yhux+95jWnJzl9+fBHk3xo5UFYpYOTfG7uEDAwMwTTmCGYzhzt2Y7s7kN25YVbNilAbbBvu7bY3c9P8vxN+n5WrKou6O6tc+eAUZkhmMYMwXTmaO+xWadTXprkiDWPb5Hksk36LgAAgH3GZpW485McW1VHV9V1kpya5DWb9F0AAAD7jE05nbK7r6yqRyd5Q5L9kry4uy/ajO/iGuPUV5jGDME0ZgimM0d7iU1Z2AQAAIDNsWk3+wYAAGD1lDgAAICBKHEAAAADUeLYqaq6cVX9xNw5YFRmCKYxQ7B7qurXq+qGtfCiqnpXVZ08dy6mU+LYUFW9dTn0N0nyz0leUlXPnDsXjMIMwTRmCFbiYd39lSQnJzkkyUOTnDlvJFZBiWNHDloO/X2TvKS7T0jy0zNngpGYIZjGDMF0tfx5ryzm6J/X7GNgShw7sqWqDk/ygCSvnTsMDMgMwTRmCKZ7Z1W9MYsS94aqukGS786ciRXYlJt9s1d4UhY3a397d59fVT+c5CMzZ4KRmCGYxgzBBFVVSf5XFqdRfry7v15VN83ilEoGp8TxfapqvyRHdPe/XUTe3R9P8kvzpYJxmCGYxgzBdN3dVfU3y1ORt+37fJLPzxiLFXE6Jd+nu69K8otz54BRmSGYxgzBypxbVbefOwSrV909dwb2QFV1RpKDkvxFkq9t29/d75otFAzEDME0Zgimq6r3J7lVkk9mMUeVxUE6t+wYnBLHhqrqLRvs7u6++zUeBgZkhmAaMwTTVdWRG+3v7k9e01lYLSUOAAD2YlV1aJLrbnvc3ZfMGIcVsLAJO1RVP5fk1tl+6J88XyIYixmCacwQTFNVv5jkGUl+KMnlSY5M8oEs5oqBWdiEDVXVnyR5YJL/lsX50/fPYvCBXWCGYBozBCvxu0numOTD3X10kpOS/L95I7EKShw78pPd/StJvtjdT0pypyRHzJwJRmKGYBozBNN9Z3lbgWtV1bW6+y1Jjp87FNM5nZId+cby59er6oeyuKfI0TPmgdGYIZjGDMF0X6qqA5P8Q5I/r6rLk1w5cyZWQIljR15bVTdK8rQk70rSSV44byQYihmCacwQTHfvLH4h8pgk/ymL23a4rnQvYHVKdqqq9k9y3e7+8txZYERmCKYxQ7D7lrcZOLa7/66qDkiyX3f/69y5mMY1cWyoqg6oqv9ZVS/o7m8lObSqfn7uXDAKMwTTmCGYrqoemeQvk/zpctfNk/zNfIlYFSWOHXlJkm9lcSF5klya5CnzxYHhmCGYxgzBdI9KcuckX0mS7v5IkkNnTcRKKHHsyDHd/dQk30mS7v5GFks8A7vGDME0Zgim+1Z3f3vbg6raksX1pQxOiWNHvl1V18ty0KvqmCx+IwrsGjME05ghmO7vq+p/JLleVf1Mklcm+duZM7ECFjZhQ1V1cpLfSXJckjdmcSj+V7v7rXPmglGYIZjGDMF0VXWtJA9PcnIWR7LfkOSFrQAMT4ljh6rqpknumMXQn9vdn5s5EgzFDME0ZgimWS4G9Lru/u7cWVgtJY4NVdVfJnlxktcbfPjBmSGYxgzBdFX1siwWB/qrJC/p7g/MHIkVcU0cO/InWdwU8iNVdWZV/djcgWAwZgimMUMwUXc/KMltk3wsyUuq6h+r6vSqusHM0ZjIkTiuVlUdlOS0LK5L+FSSFyR5WXd/Z9ZgMAgzBNOYIZiuqg5O8qAkj0nygSQ/kuQ53f3cWYOx25Q4dmh5LcKDkjw4yWVJ/jzJXZLcprvvNmM0GIIZgmnMEExTVb+Q5GFJjknyv5Oc1d2XV9UBST7Q3UfOGpDdpsSxoar66yQ/lsXAv7S7P73muQu6e+ts4WAAZgimMUMwXVX9WRarUb5tg+dO6u43zxCLFVDi2FBV3b27/+/cOWBUZgimMUMAO6bEAQAADMTqlAAAAANR4gAAAAaixLFTVfVrc2eAkZkhmMYMwQ+mqu65ZvugqnpRVb23qs6uqsPmzMZqbJk7AHuWqnrs+l1Jfruqrpsk3f3Maz4VjMMMwTRmCFbi95K8frn9jCSfTvILSe6b5E+TnDJTLlZEiWO9JyV5XZKLsvgfZ5Lsl+QGsyWCsZghmMYMwWpt7e7jl9vPqqqHzJqGlbA6JdupqlsmeWaSjyV5Und/vao+3t0/PHM0GIIZgmnMEExXVZdmMUeV5FFJjunlX/qr6r3d/RNz5mM618Sxne6+pLvvl+QdSd5UVfebOxOMxAzBNGYIVuIFWRy9PjDJWUkOTpKqulmS98yYixVxJI4dqqoDsjit5Q7d/R/nzgOjMUMwjRkC2JgSBwDs0arqpt39+blzwCiq6g5JPtDdX6mq6yV5XJLbJXl/kt/r7i/PGpDJnE7JdqrqZlX1x1X1R1V106p64nJJ2ldU1eFz54M93bplnW9UVS+0rDPsuqo6s6q2nfq1tao+nuTcqvpkVd115ngwihcn+fpy+9lJDkry+8t9L5krFKujxLHeS7P4Lc2nkrwlyTeS/HySf0jyJ/PFgmH83prtpyf5TBbLOp+fxbLOwNX7ue7+3HL7aUke2N3HJvmZLJZKB3buWt195XJ7a3c/prvf3t1PSmKRoL2AEsd6h3X3c7v7zCQ36u7fX15k/twkR84dDgaztbsf392f7O5nJTlq7kAwgGtX1bZbIF2vu89Pku7+cJL954sFQ7mwqh663P7nqtqaJFV1qyTfmS8Wq+I+cay3ttj/2dU8B2zs0OXNiivJDauq+nsXH5sh2Lk/SvK6qjozyeur6g+S/HWSk2JVPdhVj0jy7Kp6fJLPJfnHqvpUFmdaPWLWZKyEEsd6r66qA7v7q939+G07q+pHknx4xlwwim3LOiffW9b5Css6w67p7udW1fuS/Nckt8ri7yo/muRVSZ4yZzYYxXLhkl+tqhtkcfrkliSXdvdn503GqlidEgAAYCBO7WGnquq1c2eAkZkhmMYMAWxPiWNX3HzuADA4MwTTmCGANZQ4dsW75w4AgzNDMI0Zgh9QVe1XVX83dw42h2viAABgL1RVr0ny4OVCJ+xFHIljO1W1tareUlUvq6ojqupNVfWlqjq/qm47dz7Y01XVPddsH1RVL6qq91bV2VV12JzZYARmCFbqm0net5yj52z7M3copnMkju1U1XlJnpDkRkmemuS/d/dfVtVJSZ7S3XeaNSDs4arqXd19u+X2C5N8JovbDtw3yV27+5Q588GezgzB6lTVQzba391nXdNZWC0lju1U1bu7+7bL7Uu6+5YbPQdsbN1fQN/T3ceveW67x8D3M0OwWlV1vSS37O4PzZ2F1XGzb9b7ZlWdnOSgJF1Vp3T331TVXZNcNXM2GMGhVfXYJJXkhlVV/b3fljmFHXbODMGKVNUvJHl6kuskObqqjk/y5O7+xXmTMZX/GLLef0nyG0keluQeSX6qqr6U5HlJfn3OYDCIFyS5QZIDk5yV5OAkqaqbJXnPjLlgFGYIVueJSU5M8qUk6e73JDl6zkCshtMpAQBgL1RV/9Tdd1h3ucx7u/sn5s7GNI7Escuq6qFzZ4ARVNWJVXX75fZxVfXYqrrX3LlgRFV1l+UMnTx3FhjQhVX1y0n2q6pjq+q5Sd4xdyimcySOXbZ+oRPg+1XVE5L8bBbXHL8pyR2SvDXJTyd5Q3efMV862PNV1XndfeJy+5FJHpXkVUlOTvK33X3mnPlgJFV1QJLfyWJ+kuQNWaw2/s35UrEKShzbqar37uipJLfq7v2vyTwwmqp6X5Ljk+yfxdLot+juryxXB/snp7DA1Vt32tf5Se7V3VdU1fWTnNvdt5k3IYyjqm7b3e+eOwerZ3VK1jssiwVNvrhuf8Xhd9gVV3b3VUm+XlUf6+6vJEl3f6OqvjtzNhjBtarqxllc8lHdfUWSdPfXqurKeaPBcJ5ZVYcneWWSc7r7orkDsRpKHOu9NsmBy9WLtlNVb73m48Bwvl1VB3T315OcsG1nVR2URImDnTsoyTuz+OVhV9XNuvszVXXgch+wi7r7p5Yruz4gyfOr6oZJ/qK7nzJzNCZyOiXAClXV/t39rQ32H5zk8O5+3wyxYHjLa3sO6+6L584CI6qq2yT5rSQP7O7rzJ2HaZQ4AADYC1XVv0vywCT3T/K5JOck+avuvnzWYEymxAEAwF6oqs5N8vIkr+zuy+bOw+oocQAAsBdarox8TJJO8jG3Fth7uNk3AADsRapqS1U9NcklSc5K8rIkn6qqp1bVtedNxyoocQAAsHd5WpKbJPnh7j5hee/FY5LcKMnTZ03GSjidEgAA9iJV9ZEkt+p1f9Gvqv2SfLC7j50nGaviSBwAAOxden2BW+68Kovr4xicEgcAAHuX91fVr6zfWVUPSvLBGfKwYk6nBACAvUhV3TzJXyf5RpJ3ZnH07fZJrpfkPt39LzPGYwWUOAAA2AtV1d2T3DpJJbmou988cyRWRIkDAAAYiGviAAAABqLEAQAADESJA2AIVXWfquqq+rG5swDAnJQ4AEZxWpK3Jzn1mvzS5c1xAWCPocQBsMerqgOT3DnJw7OmxFXVtarqeVV1UVW9tqpeV1X3Wz53QlX9fVW9s6reUFWHb/C5x1TVuVV1flU9uaq+utx/t6p6S1WdneR9y32PraoLl38es9x3VFVduObzfrOqnrjcfmtV/UFVvWP5nhM37V8QAPsUJQ6AEZyS5PXd/eEkX6iq2y333zfJUUluk+QRSe6UJFV17STPTXK/7j4hyYuTnLHB5z47ybO7+/ZJLlv33IlJfqe7j6uqE5I8NMkdktwxySOr6ra7kPv63f2TSX5tmQEAJlPiABjBaUnOWW6fs3ycJHdJ8sru/m53fybJW5b7fzTJjyd5U1W9J8njk9xig8+9U5JXLrfPXvfced198ZrveVV3f627v5rFTXT/wy7kfnmSdPfbktywqm60C+8BgKu1Ze4AAHB1quqmSe6e5MerqpPsl6Sr6reyuIHthm/L4sa2d5rw1V9b93kbuTLb/0L0uuueX38zVjdnBWAyR+IA2NPdL8mfdfeR3X1Udx+R5OIsjo69PckvLa+NOyzJ3Zbv+VCSQ6rq306vrKpbb/DZ5yb5peX21S2Y8rYkp1TVAVV1/ST3SfIPST6b5NCqumlV7Z/k59e974HL779Lki9395d/oH9yANiAI3EA7OlOS3Lmun1/leSXkzwqyUlJLkzy4ST/lEVZ+vZygZPnVNVBWfz/7g+SXLTucx6T5GVV9RtJ/k+SDUtWd7+rql6a5Lzlrhd297uTpKqevPzei5N8cN1bv1hV70hywyQP+0H+oQFgR6rbmR0AjKuqDuzury5PuzwvyZ2X18ftynsPSPKN7u6qOjXJad197xXlemuS3+zuC1bxeQCwjSNxAIzutcsFQ66T5Hd3tcAtnZDkD6uqknwpjpYBMABH4gAAAAZiYRMAAICBKHEAAAADUeIAAAAGosQBAAAMRIkDAAAYiBIHAAAwkP8PK3JKkUPrKpcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1080x720 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "_, ax = plt.subplots(figsize=(15, 10))\n",
    "data.pivot_table(index=['Age group'], columns=['Smoker', 'Status'], aggfunc='size').plot(kind='bar', stacked=True, ax=ax)\n",
    "ax.set_title(\"Number of alive/dead women after 20 years, according to their smoking habits and age\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One can see that most women that were over 65 in 1972 are dead twenty years after and that, at that a great majority of polled older women were non-smokers. For the other age groups results look similar but are difficult to interpret.\n",
    "\n",
    "Let's therefore compute the mortality rates for smokers and non-smokers in different age groups."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Age group: 18-34 years\n",
      "\tMortality rate for smokers:\t3.70% ± 2.69%\n",
      "\tMortality rate for non smokers:\t2.64% ± 2.09%\n",
      "\n",
      "Age group: 35-54 years\n",
      "\tMortality rate for smokers:\t17.03% ± 4.87%\n",
      "\tMortality rate for non smokers:\t9.95% ± 4.24%\n",
      "\n",
      "Age group: 55-64 years\n",
      "\tMortality rate for smokers:\t44.35% ± 9.08%\n",
      "\tMortality rate for non smokers:\t33.06% ± 8.38%\n",
      "\n",
      "Age group: Over 65 years\n",
      "\tMortality rate for smokers:\t85.71% ± 9.80%\n",
      "\tMortality rate for non smokers:\t85.49% ± 4.97%\n",
      "\n"
     ]
    }
   ],
   "source": [
    "def ci_per_age(age_group):\n",
    "    pop = data[data['Age group'] == group]\n",
    "    print(f\"Age group: {group}\")\n",
    "    \n",
    "    num_smokers = sum(pop['Smoker'] == \"Yes\")\n",
    "    num_dead_smokers = sum(np.logical_and(pop['Smoker'] == \"Yes\", pop['Status'] == \"Dead\"))\n",
    "    rate_smokers = num_dead_smokers / num_smokers\n",
    "    ci_smokers = z * (rate_smokers * (1 - rate_smokers) / num_smokers) ** 0.5\n",
    "    print(f\"\\tMortality rate for smokers:\\t{rate_smokers:.2%} \" + u\"\\u00B1\" + f\" {ci_smokers:.2%}\")\n",
    "\n",
    "    num_non_smokers = len(pop) - num_smokers\n",
    "    num_dead_non_smokers = sum(np.logical_and(pop['Smoker'] == \"No\", pop['Status'] == \"Dead\"))\n",
    "    rate_non_smokers = num_dead_non_smokers / num_non_smokers\n",
    "    ci_non_smokers = z * (rate_non_smokers * (1 - rate_non_smokers) / num_non_smokers) ** 0.5\n",
    "    print(f\"\\tMortality rate for non smokers:\\t{rate_non_smokers:.2%} \" + u\"\\u00B1\" + f\" {ci_non_smokers:.2%}\")\n",
    "    \n",
    "for group in sorted(data['Age group'].unique()):\n",
    "    ci_per_age(group)\n",
    "    print()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, the mortality rate is considerably higher for smokers, especially for people between 35 and 65 years old. This might seem like a contradiction, as before the rate was higher for non-smokers. However, from the previous bar chart it is clear that the percentage of polled smokers/non-smokers is different in different age groups, in particular for older women as mentioned. In addition, the fact that most polled women over 65 are non-smokers can be an argument in favor of the hypothesis that smoking is dangerous for health, but this can't be proven through statistics.\n",
    "\n",
    "---\n",
    "\n",
    "The age groups are fixed a-priori. In order to have more flexible results and reduce the introduced bias it is possible to try to perform a logistic regression, studying the probability of death in the two groups according to the age."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Smoker</th>\n",
       "      <th>Status</th>\n",
       "      <th>Age</th>\n",
       "      <th>Age group</th>\n",
       "      <th>Death</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>21.0</td>\n",
       "      <td>18-34 years</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>19.3</td>\n",
       "      <td>18-34 years</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>No</td>\n",
       "      <td>Dead</td>\n",
       "      <td>57.5</td>\n",
       "      <td>55-64 years</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>No</td>\n",
       "      <td>Alive</td>\n",
       "      <td>47.1</td>\n",
       "      <td>35-54 years</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Yes</td>\n",
       "      <td>Alive</td>\n",
       "      <td>81.4</td>\n",
       "      <td>Over 65 years</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  Smoker Status   Age      Age group  Death\n",
       "0    Yes  Alive  21.0    18-34 years      0\n",
       "1    Yes  Alive  19.3    18-34 years      0\n",
       "2     No   Dead  57.5    55-64 years      1\n",
       "3     No  Alive  47.1    35-54 years      0\n",
       "4    Yes  Alive  81.4  Over 65 years      0"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "death = []\n",
    "for status in data['Status']:\n",
    "    if status == \"Alive\":\n",
    "        death.append(0)\n",
    "    elif status == \"Dead\":\n",
    "        death.append(1)\n",
    "assert len(death) == len(data['Age'])\n",
    "data['Death'] = death\n",
    "data.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It will be done using statsmodels: sklearn does not allow to derive prediction confidence intervals and uses regularization by default."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "smokers = data[data['Smoker'] == \"Yes\"]\n",
    "non_smokers = data[data['Smoker'] == \"No\"]\n",
    "log_smokers = LogisticRegression().fit(smokers['Age'].values.reshape([-1, 1]), smokers['Death'])\n",
    "log_non_smokers = LogisticRegression().fit(non_smokers['Age'].values.reshape([-1, 1]), non_smokers['Death'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "At this point, we can plot the predicted probabilities of death as a function of the age."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlYAAAF1CAYAAAAqdaQaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xd0VNXax/HvTiGBAKGDEHqT3kJTmkqzgKD0DipivdjvtRe8Yq+vF1Q6UgQRaQp2UJCO9F5DTYAEkhBI2e8fZ4KTmApJJiS/z1pZmZmzZ59nTptn9tlnH2OtRURERESunpenAxARERHJK5RYiYiIiGQRJVYiIiIiWUSJlYiIiEgWUWIlIiIikkWUWImIiIhkESVW+YQxpq0xZpen40iJMaaDMSYkh+ZVxRhjjTE+V/h+a4ypkcq0gcaYZSmVNcaMM8a8cGVRZzrGB4wxJ40xkcaYkhkof9AY0zGbYsm2uuXKGWOGGWN+d3seaYyp5smYPCWt4092HjeNMb8aY+69wvdONsaMSWN6vl2fuYESq1zO9cV0wbWjJP59koH3JUkArLUrrLW1synGNHfy/MJa+6W1tnMq00ZZa1+D7E0kjTG+wHtAZ2ttYWvt6eyYTyrz1nZwjXJtK/uzul5jzMvGmOlZXW9Oyc7jZnZyX5/aL3PeFf1qlxzXzVr7o6eDuBYYY3ystXGejsODygL+wDZPB+JJ2g7+ZozxttbGezoOkfxCLVbXMGNMDWPMb8aYCGNMmDFmtuv15a4if7lauPombyVxtYQ9ZYzZbIyJMsZMMMaUNcZ8Z4w5b4z50RhT3K38HGPMCde8lhtj6rleHwkMBJ52zWuh6/XyxpivjTGhxpgDxphH3eoq6PoVddYYsx1ons7ntMaYR40x+12f821jjJdr2jBjzB/GmPeNMWeAl40xXsaY540xh4wxp4wxU40xgcmqHWGMOWaMOW6MecJtXi2MMauMMeGuaZ8YYwoke+9tacTyOylI/NVojAkAvgPKu7VAljfGRLuftjPGNHMtO98U6vIzxnzgiv+Y67GfMaYWkHjaItwY83MqsQx2LZvTxpjnkk3zMsb82xizzzX9K2NMCbfpmdoOXBq7trMIY8xsY4x/KnElrsuPXWV3GmNucZse6NpOjxtjjrqWp3ey917eDlKoP811a4zpbIzZ5Zr3p8bZt+51mz7CGLPDtd0uNcZUTulzpDDf6saYn13LM8wY86Uxppjb9IrGmHmu9X3auLVIG2Puc83zvDFmuzGmqev1OsY5lRRujNlmjOnu9p7Jxpj/GWOWGGOigJuMMSWNMQuMMeeMMWuA6slidD9tPdkY83/GmMWu+a42xlR3K5vmcnIr1xV4Fujr2ib+cr1e3hXLGWPMXmPMfWksu9tcn/u8a50/6Xq9gzEmxBjztHH28ePGmB6u8rtddT/rVk+K+0wq83zUNc8gk/Jx88nUtmdXPMdd87jXpNF1wKWya7s9b4xZZowp5VZXivuam1LGmB9c7/3NfXtMnK9J/fj8jGt5nnety1uQrGOt1V8u/gMOAh1TmTYTeA4nQfYH2rhNs0ANt+cdgJBk9f6J08JRATgFbACaAH7Az8BLbuVHAEVc0z4ANrlNmwyMcXvuBawHXgQKANWA/UAX1/SxwAqgBFAR2OoeWwqf0wK/uMpXAnYD97qmDQPigEdwWmALumLd65pvYWAeMM1VvoqrvplAANAACE1cxkAzoJWrrirADmB0JmL5PaV14L6Mkq8L12tLgAfcnr8PfJzK8njVte7KAKWBlcBryT6fTyrvrQtEAu1c6/I91/JL/PyjXXUHuaaPB2ZeyXbgtp2tAcq7ltkOYFQqsSWuy8cAX6AvEAGUcE2f74onwPXZ1wD3p7YdpFB/qusWKAWcA+5yTf8XEOu2bnvgbFN1XNOfB1ZmcB+uAXRyLbPSwHLgA9c0b+Av1/oOwG0/BnoDR3F+eBhXPZVdy2YvTtJSALgZOA/UdlsPEcCN/H1smAV85ZpHfVe9aW2rZ4AWrs/6JTArI8sphc/+MjA92Wu/AZ+64mqMs//dksr7jwNtXY+LA03d9qE4nGOML3Cfq54ZONtnPSAGqJaBfaYDrv0ReAHnOFg6jeNmitsz0BU44Zp3IWAayY7DyT7br8A+oBbOcetXYGwm9rXz/L0ff5jO+nQ/PtcGjgDl3Y4Z1TOyLesvY38eD0B/6awgZ0eOBMLd/u5zTZsKfAYEpfC+jCRWA92efw38z+35I8D8VGIq5qo/0PU8+Y7bEjic7D3/ASa5Hu8HurpNG0n6iZV7+QeBn1yPh6Uwr5+AB92e18Y5+Cd+oVrgerfpbwETUpn3aOCbTMRypYlVX+AP12NvnAN0i1Ri2gfc5va8C3DQ9Tjx86WWWL2I60vS9TwAuMTfidUO3L7kgOsSl11mtwO37WxQsmU9LpXYhgHHAOP22hpgMM4PgIu4JUxAf+CX1LaDDOxbl9ctMARY5TbN4Hz5JCZW3wH3uE33AqKBypmZp+u9PYCNrsetcRKClJbvUuBfKbze1rV9eLm9NhN42W09THWb5u1ah+7b/H/T2Va/cJt2G7AzI8sphVhfxi2xwvkhFQ8UcXvtDWByKu8/DNwPFE32egfgAuDtel7E9RlaupVZD/TIwD7TASfRfA/4Hdf2nNK+ShrbMzAReMNtWg3ST6yed3v+IPB9KmVT2tfc9+PCruVaMZX1OSZZXKeAjoBvZrdf/aX/p1OB14Ye1tpibn+fu15/GufAtsZ1OmBEJus96fb4QgrPC4PTR8MYM9Y4p4fO4RxcwPn1mpLKOKe6whP/cH5dl3VNL49zME50KAOxJi9fPpVpifW713kIJ6kq6/ZaivUZY2oZYxa5muDP4XwBJf+cacVypb4F6hrnSp5OQIS1dk0qZVP6fBmNIcmyt9ZGAe4d3CsD37ittx04B+yyV7AdJDrh9jga13aViqPWdfR3SfxsiS01x91iG4/TApEo+XaQRDrrNvlysYD7BQaVgQ/d5n0GZ9+rkNY8XfMtY4yZ5Tr1cg6Y7jbfisAhm3J/sIo4CUFy5YEj1toEt9cOJYvFfVmUxtn+M7PPpbbO0ltO6SkPnLHWnk8WS2rL8W6cxO6Q63RXa7dpp+3ffccuuP6neAwj/X2mGM4PvDestRHpfIYMLRvS2R7TqiuD+5r7eojE2SbTPQ5Ya/fi/Kh4GTjl2jaz4hgmLkqsrmHW2hPW2vusteVxftV9ms75/Cs1ALgT5xdOIE6rCDhfLOD8OnJ3BDiQLBksYq29zTX9OM6XRqJKGYghefljbs+Tz/8Yzhehe/k4kh50U6vvf8BOoKa1tihOQmhIKq1YMiJ5vFhrY3BO1QzEaaGZlsb7U/p8GY0hybI3xhQC3IdkOALcmmzd+Vtrj5L57eBKVDDGuC/vxM92BKfFqpRbXEWtte79TtKbf1rr9jjO6U8AXDEEub33CM5pR/flUtBauzIDn+kNV2wNXfMd5DbfI0Alk/LwH0dI1hfK5RhQ0bj69rlUwml1SeS+LEJxtv/M7nMpSW85JZfSvlnCGFMkWSxHSYG1dq219k6cBHo+zj5yJdLbZ84CdwCTjDE3XuE8kiwbki7vzEpvX0tSvzGmMM6pyZSOAykdb2ZYa9vgLBMLvHkVsUoySqyuYcaY3saYxB35LM4OkvgL7iROH6OsUATnS+00Tt+B/yabnnxea4Bzrg6SBV2/vuobYxI7qX8F/McYU9wV/yMZiOEpV/mKOP06ZqdRdibwmDGmquuA819gdrJWgReMMYVcHUKHu9VXBKcPSaQx5nrggauMJSUngZLmnx3qp+Kc0uqO06qRmpnA88aY0q7Ori+mU97dXOAOY0wb43TcfpWkx4FxwOuJHWFd87jTNS2z28GVKAM8aozxNcb0xunTtMRaexxYBrxrjClqnE721Y0x7TNRd1rrdjHQwDgdoH2Ah4BybtPH4WyziZ31A13x4Xr+qzHm5TTmG4lzQUEF4Cm3aWtwvpDHGmMCjDH+bl/sXwBPGudCBuPqjFwZWA1E4XRI9jXGdAC64fSj+gdXq848nAs7Chlj6gJD01pQaUhvOSV3EqiSmARaa4/g9G96w/VZGwL34PTjSsIYU8A4Y8MFWmtjcdbdlV7dmO4+Y639FeeHzTfGmJZXMI+vgOHGubCgkGseVyq9fQ2ci2gS9+PXgNWu5Ztckv3SGFPbGHOzcTrvx+C07Omq0SykxOrasNAkHcfqG9frzYHVxphIYAFOf4wDrmkvA1Ncpy76XOX8p+I0nR8FtuN0AnU3Aec0VrgxZr7rQN4Np2PqASAM50siMZF4xVXfAZwvy7RaZxJ9i9NnYhPOwX1CGmUnuupc7ppHDP9M3n7D6QD8E/COtTZxYM8ncX4tngc+J+WkKTOx/IO1difOgX6/a5mVd73+B5AAbLDWHkyjijHAOmAzsAWns22Gxqmx1m7D+TKcgfOFfpakp3I+xNmWlhljzuOs68QvmUxtBxmJJwWrgZo428zrQC/791hcQ3A6a293xT0Xpw9YRqW6bq21YTidxd/C+TKri7OML7qmf4Pzq36W69TMVuBWt7orAn+kMt9XgKY4HcoX4yQ5ifNN3Fdq4PQnCsHpb4e1do5rGcxwxTwfpyP/JZzk+1ac5fQpMMS1XaXmYZzTTCdw+txMSqNsqtJbTimY4/p/2hizwfW4P04LzDHgG5yLZH5I5f2DgYOuZT4Kp7XvSmRon3HFMRxYYIxplpkZWGu/Az7CubhlL7DKNSm1ZZOW9PY1cLaLl3BOATbDSQpTkny/9MO5gCgMZ3sog9N6K1nEJO3OIJL7GGMszumbvZ6OJbsZZ4iEGdbaLzwdS04zxgzD6QTdJhfE4oWT5Ay01v6STtkgYI61tnVa5fKizCyn/MYYUwcnAfdLpQ+d5FFqsRLJJVynSpuS+VOLkgWMMV2MMcVcp0gS+1+l1FKQhLU2JD8lVVe6nPIDY0xP1ynM4jgtnAuVVOU/SqxEcgFjzBTgR5xxlc6nV16yRWucq/DCcE7P9bDWXkj7LfmSllPq7se5WGAfTr+llPpoSh6nU4EiIiIiWUQtViIiIiJZRImViIiISBZJaVC6JIwxE3EGTjtlra2fwnSDc4n2bTgjxw6z1m5IXi65UqVK2SpVqmQ6YBEREZGctn79+jBrben0yqWbWOGMefIJzrgaKbkVZ9yZmjjj3fyPv8e9SVWVKlVYt25dBmYvIiIi4lnGmIzcfi39U4HW2uU4A5Cl5k6cG35aa+2fQDFjTGYG7RMRERHJE7Kij1UFkt5sMoRUbqhpjBlpjFlnjFkXGhqaBbMWERERyT2yIrFKfoNaSOVmqNbaz6y1wdba4NKl0z1NKSIiInJNyUgfq/SEkPQu3kGkfIftdMXGxhISEkJMTEwWhCVZwd/fn6CgIHx9fT0dioiISK6XFYnVAuBhY8wsnE7rEa470WdaSEgIRYoUoUqVKjgXG4onWWs5ffo0ISEhVK1a1dPhiIiI5HoZGW5hJtABKGWMCcG5m7YvgLV2HLAEZ6iFvTjDLQy/0mBiYmKUVOUixhhKliyJ+sOJiIhkTLqJlbW2fzrTLfBQVgWkpCp30foQERHJOI28noLXX3+devXq0bBhQxo3bszq1auvqr5ff/2VO+64I4uiExERkdwqK/pY5SmrVq1i0aJFbNiwAT8/P8LCwrh06ZLH4omLi8PHR6tJRETkWqAWq2SOHz9OqVKl8PPzA6BUqVKUL1+eKlWq8Oyzz9K6dWuCg4PZsGEDXbp0oXr16owbNw5wOns/9dRT1K9fnwYNGjB79ux/1L927VqaNGnC/v37iYqKYsSIETRv3pwmTZrw7bffAjB58mR69+5Nt27d6Ny5M8ePH6ddu3Y0btyY+vXrs2LFipxbICIiIpJhubYp5JWF29h+7FyW1lm3fFFe6lYvzTKdO3fm1VdfpVatWnTs2JG+ffvSvn17ACpWrMiqVat47LHHGDZsGH/88QcxMTHUq1ePUaNGMW/ePDZt2sRff/1FWFgYzZs3p127dpfrXrlyJY888gjffvstlSpV4tlnn+Xmm29m4sSJhIeH06JFCzp27Ag4LWebN2+mRIkSvPvuu3Tp0oXnnnuO+Ph4oqOjs3S5iIiISNbItYmVpxQuXJj169ezYsUKfvnlF/r27cvYsWMB6N69OwANGjQgMjKSIkWKUKRIEfz9/QkPD+f333+nf//+eHt7U7ZsWdq3b8/atWspWrQoO3bsYOTIkSxbtozy5csDsGzZMhYsWMA777wDOFdFHj58GIBOnTpRokQJAJo3b86IESOIjY2lR48eNG7cOKcXi4iISO4TfQZO7QD/QChX39PRALk4sUqvZSk7eXt706FDBzp06ECDBg2YMmUKwOXTg15eXpcfJz6Pi4vDuUAyZddddx0xMTFs3LjxcmJlreXrr7+mdu3aScquXr2agICAy8/btWvH8uXLWbx4MYMHD+app55iyJAhWfZ5RUREcrXYGAjdASe3wYmtcHIrhO6EKNdwQE2HQvePPBujS65NrDxl165deHl5UbNmTQA2bdpE5cqV2bJlS7rvbdeuHePHj2fo0KGcOXOG5cuX8/bbb7Nz506KFSvGhAkT6Ny5MwEBAXTo0IEuXbrw8ccf8/HHH2OMYePGjTRp0uQf9R46dIgKFSpw3333ERUVxYYNG5RYiYhI3mMtnD/hJE4nt/6dRIXtARvvlPEtBGXqQK0uULoOlL4+17RWgRKrf4iMjOSRRx4hPDwcHx8fatSowWeffcaiRYvSfW/Pnj1ZtWoVjRo1whjDW2+9Rbly5di5cycAZcuWZeHChdx6661MnDiRF154gdGjR9OwYUOstVSpUiXF+fz666+8/fbb+Pr6UrhwYaZOnZrln1tERCTHRZ6Coxvg6Hrn7/hfEB329/TAilC2PtTp5vwv1wCKVwEvb4+FnB6T1umr7BQcHGzXrVuX5LUdO3ZQp04dj8QjqdN6ERGRq3YxEo5v+juJOroRIpx+xRgvKFMXyjeGsg2cFqiy9aBgcc/G7MYYs95aG5xeObVYiYiISNYLPwyHVsKhPyBkndMnyiY404pVhqBgaDkSKjSD6xpBgYC067tGKLESERGRq2MtnN7rJFGHVjp/EUecaX6BULE51OnuJFEVmkJAKc/Gm42UWImIiEjmWOtcoXfoj7+TqcQr9AJKQ+Ub4IZHnP9l6ubqPlFZTYmViIiIpO/ccdj/C+z7Bfb/ClGnnNcDK0H1W5wkqvKNULI6GOPRUD1JiZWIiIj8U2wMHPwd9v3kJFOhO5zXA0pDtQ5Q7Sao2haKVfJYiBHRsaw/fIYSAX40rljMY3G4U2IlIiIijnPHYc8y2L3UaZ2KjQYff6jUGhr3h+o3Q5l64OWZWw2fOh/D6v1nWH3gNGsPnGXXyfMA3NW0Ao0r5o67kiixSsYYw+OPP867774LwDvvvENkZCQvv/yyZwNLRYcOHXjnnXcIDk73ClAREZGkrHXGjtq1BHZ/7zwGZ/yoxgOdQTirtAHfgh4JLyzyIqv3n2HV/jD+3H+GvaciASjs50PTysW5o+F1BFcpkWtaq0CJ1T/4+fkxb948/vOf/1CqVN69agEgPj4eb+/806FQRESAhAQIWQs7Fjh/4YedcaSCWsAtL0Gtrs7I5h7oJ3Um6hKr959m1f7T/Ln/NLtPOolUQAFvmlctQa9mQbSuVpJ65Yvi4+2ZVrP0KLFKxsfHh5EjR/L+++/z+uuvJ5l26NAhRowYQWhoKKVLl2bSpElUqlSJYcOGUbRoUdatW8eJEyd466236NWr1z/qnjNnDq+88gre3t4EBgayfPlyJk+ezPz584mPj2fr1q088cQTXLp0iWnTpuHn58eSJUsoUaIEmzZtYtSoUURHR1O9enUmTpxI8eJ/D5yWkJDA8OHDqVixImPGjGHZsmW89NJLXLx4kerVqzNp0iQKFy5MlSpVGDFiBMuWLePhhx/m1KlTjBs3Dh8fH+rWrcusWbOyfRmLiEgOS4h3+kvtWAA7FkHkCfDyheo3QbunofZtEFAyx8OKjU9gw6Gz/LY7lN92h7Lt2DkAChXwJrhKCXo0qUDraiWpXyEQ31yaSCWXexOr7/4NJ9K/P1+mlGsAt45Nt9hDDz1Ew4YNefrpp5O8/vDDDzNkyBCGDh3KxIkTefTRR5k/fz4Ax48f5/fff2fnzp107949xcTq1VdfZenSpVSoUIHw8PDLr2/dupWNGzcSExNDjRo1ePPNN9m4cSOPPfYYU6dOZfTo0QwZMoSPP/6Y9u3b8+KLL/LKK6/wwQcfABAXF8fAgQOpX78+zz33HGFhYYwZM4Yff/yRgIAA3nzzTd577z1efPFFAPz9/fn9998BKF++PAcOHMDPzy9JTCIico2z1rldzJY5sPVr5yo+n4JQsyPUuRNqdQb/wBwP68iZaJbvCeW3XaGs3HeayItx+HgZmlYuzpOda9G6eikaBl07iVRyuTex8qCiRYsyZMgQPvroIwoW/Pu88qpVq5g3bx4AgwcPTpJ49ejRAy8vL+rWrcvJkydTrPfGG29k2LBh9OnTh7vuuuvy6zfddBNFihShSJEiBAYG0q1bNwAaNGjA5s2biYiIIDw8nPbt2wMwdOhQevfuffn9999/P3369OG5554D4M8//2T79u3ceOONAFy6dInWrVtfLt+3b9/Ljxs2bMjAgQPp0aMHPXr0uLIFJiIiuUfobtg610mozuwH7wJOX6n6vaBmpxwf4TwmNp7VB87w265Qftt9in2hUQBUKFaQ7o3L075WaW6oXpIi/r45Gld2yb2JVQZalrLT6NGjadq0KcOHD0+1jHE7/+zn53f5ceL9F5977jkWL14MwKZNmxg3bhyrV69m8eLFNG7cmE2bNv3jvV5eXpefe3l5ERcXl26sN9xwA7/88gtPPPEE/v7+WGvp1KkTM2fOTLF8QMDfO9XixYtZvnw5CxYs4LXXXmPbtm34+OTezUJERFIQfQa2zIVN010d0A1UbQdtHnduYFwwZzt3n426xM87T7Fs+wmW7w7jQmw8fj5etKxWkgEtK9O+Vmmqlw5I8j2aV+gbNBUlSpSgT58+TJgwgREjRgBOAjNr1iwGDx7Ml19+SZs2bdKs4/XXX0/ST2vfvn20bNmSli1bsnDhQo4cOZKhWAIDAylevDgrVqygbdu2TJs27XLrFcA999zD8uXL6d27N9988w2tWrXioYceYu/evdSoUYPo6GhCQkKoVatWknoTEhI4cuQIN910E23atGHGjBlERkZSrFjuubpCRERSkRDvDNS5cTrsXATxl6BcQ+jyX6h3FxS9LkfDOXImmh+2n2TZ9hOsPXiW+ARLuaL+9GoWxC11ytCqWkn8ffP+BVNKrNLwxBNP8Mknn1x+/tFHHzFixAjefvvty53XM+Opp55iz549WGu55ZZbaNSo0eVWq/RMmTLlcuf1atWq/WPejz/+OBEREZeTvsmTJ9O/f38uXrwIwJgxY/6RWMXHxzNo0CAiIiKw1vLYY48pqRIRye3OHoJNX8KmGc79+AoWh+ARzvAI1zXMsTCstWw7do5l20/yw/aT7DjudDyvXbYID7SvTud6ZWlQITBPtkqlxSSetsppwcHBdt26dUle27FjB3Xq1PFIPJI6rRcREQ9LSIC9P8Laz2HPD85r1W+GJoOcK/p8/XMkDGst24+fY+Ffx1m0+RghZy9gDDSvXIJOdcvSqW5ZqpTK2T5cOcUYs95am+6gkWqxEhERya2izzin+tZNgLMHIaAMtHsKmg6BYhVzLIy9pyJZ+NcxFm4+xv7QKHy8DG1qluLRm2tyc50ylCrsl34l+YQSKxERkdzm+F+w+jPn6r64GKh0A9zyIlzfDXwK5EgIR85Es3DzMRb+dZwdx89hDLSqWpJ721Sja/1ylAjImTiuNUqsREREcoPE030rP4KDK8A3ABr1h+b3Qrn6ORJCWORFFmxyWqY2HnbGNmxaqRgvdavL7Q2uo0zRnDnleC3LdYmVtTbfdXTLzTzVB09EJN+IjYEtX8HKTyBsFxStAJ1eg2ZDc2QAz0txCfy88yRz1x/l112niEuw1L2uKP++9Xpub3AdFUsUyvYY8pJclVj5+/tz+vRpSpYsqeQqF7DWcvr0afz99QtFRCTLRZ9x+k6t/swZFb1cA7jrc6jXE7yzd7DMxCv65q4P4dtNRzkbHUuZIn7c07YqvZoGUbNskWydf16WqxKroKAgQkJCCA0N9XQo4uLv709QUJCnwxARyTuiwmDVJ7Dmc7gUCTU6wg2PQNX22X7j44joWL7ZGMKstUfYeeI8BXy86Fy3LL2aBdGmRqlce2Pja0muSqx8fX2pWrWqp8MQERHJeudPOv2n1k2E2AtOy1S7J6FsvWydrbWW9YfOMmPNYRZvPs7FuAQaBQUypkd9ujUsT2ChvHErmdwiVyVWIiIiec65Y/DHh7B+sjM6eoPe0PZJKF0r3bdejfDoS8zbcJSZaw6z51Qkhf186B0cRL/mlahfIedvvpxfKLESERHJDpGhsOIdp4XKJkDDftD2cShZPVtnu+P4OaasPMj8TUeJiU2gUcVivHV3Q25veB0Bfvraz25awiIiIlkpJsK5wm/V/zljUDUe4JzyK14l22YZG5/Asm0nmbLyIGsOnsHf14ueTSowqFVl6pVX61ROUmIlIiKSFWIvOB3Sf38PLpx1+lDd9DyUqpFtszwdeZGZaw4z/c/DnDgXQ8USBXnutjr0Ca6ovlMeosRKRETkasTHOTdF/nUsnD8G1W+BW16A8k2ybZb7QiP5YsUBvt4QwqW4BNrWLMXrPevToXYZvL00XJEnKbESERG5Uvt+gaXPwqntENQc7voMqrbNlllZa1l94AxfrNjPjztO4efjxd1Ng7inTVVqlCmcLfOUzFNiJSIiklmn98HS52D3d1CsMvSZCnW6Z8s4VHHxCSzZeoIvVuxnc0gEJQIKMLpjTQa3qkxJ3fw411FiJSIiklEXwmH527B6PPj4Q8dXoOUo8M36O1RcjIvn6/VHGffbPg6fiaZaqQBe71mfu5sG4e/rneXzk6yhxEpERCQ9CfGwfhL8/LrTMb3pELjUqeZkAAAgAElEQVT5eShcJstnFX0pjhmrD/P5iv2cPHeRRkGBPHd7MzrVKYuX+k/lekqsRERE0nJ0Ayx6DI5vgiptoct/4bqGWT6biAuxTF15kIl/HOBsdCytqpXg3d6NubGG7p97LVFiJSIikpIL4fDzGFj7BRQuC70mQr27srwfVcSFWCb8foBJvx/g/MU4br6+DA/dVJ1mlUtk6XwkZyixEhERcWctbJnjdE6PDoOW98NNz4J/1g60eT4mlkl/HOTzFfs5HxNH13rleOSWGhrQ8xqnxEpERCRR2B5Y/DgcWA4VmsHAOVC+cZbOIupiHJNXOglVeHQsneqWZXTHmkqo8gglViIiIvFxsPJDZ5BP34Jw+3vQbBh4Zd3VdzGx8UxZeZDxy/dzJuoSt1xfhtEda9EgSAlVXqLESkRE8rcTW+Hbh5zO6XXvhNveydKr/eLiE5i7PoQPftzDiXMxtKtVmsc71aJxxWJZNg/JPZRYiYhI/hR3ybmv3/J3oGAx6D0F6vXIsuqttSzbfpK3vt/JvtAomlQqxof9GtOyWsksm4fkPkqsREQk/zm2Eb59GE5uhQZ9oOtYCMi6hGf1/tO8+f1ONhwOp3rpAMYNakaXemU1bEI+oMRKRETyj/hY+O1NWPGec7qv/yyofWuWVb/3VCRvLNnBTztPUbaoH2PvakCvZkH4eHtl2Twkd1NiJSIi+UPYHph3n9Na1XigM9Bnwazp53Q26hIf/rSH6X8ewt/Xm6e71mb4DVUpWEC3nslvlFiJiEjeZi2smwBLn3fu6ddnGtTtniVVx8YnMG3VIT78aQ/nY2Lp16ISj3eqRSndHDnfUmIlIiJ5V+Qppy/VnqVQ/Wa481Moet1VV2ut5acdp/jvkh3sD4uibc1SPHd7Ha4vVzQLgpZrmRIrERHJm3YugQWPwKVIuPUtaH4feF19X6e9pyJ5ZeE2VuwJo1rpACYOC+am2mXUMV2ADCZWxpiuwIeAN/CFtXZssumBwHSgkqvOd6y1k7I4VhERkfTFxsCy55x7/JVrAHd9DmXqXHW1URfj+OjnPUz8/QD+vt68eEddBreujK86poubdBMrY4w38H9AJyAEWGuMWWCt3e5W7CFgu7W2mzGmNLDLGPOltfZStkQtIiKSktP7YM4wOLEZWj8Mt7wIPlfX38lay6LNx3l98Q5OnIuhV7Mgnul6PaWLqB+V/FNGWqxaAHuttfsBjDGzgDsB98TKAkWM0w5aGDgDxGVxrCIiIqnb+jUs+JdzG5osGkZhz8nzvLRgGyv3naZe+aL838AmNKtcIguClbwqI4lVBeCI2/MQoGWyMp8AC4BjQBGgr7U2IUsiFBERSUtsDCx91rnyL6g59JoExSpeVZUXLsXz4U97+GLFfgoV8Oa1O+sxoGVlvL3Uj0rSlpHEKqWtyCZ73gXYBNwMVAd+MMassNaeS1KRMSOBkQCVKlXKfLQiIiLuTu+DOUPhxBa44VHn1J+371VV+euuUzw/fyshZy/Qu1kQ/771ekpq+ATJoIwkViGAe+ofhNMy5W44MNZaa4G9xpgDwPXAGvdC1trPgM8AgoODkydnIiIiGbd9Acx/ELx9YMBXUKvLVVV36nwMry3awcK/jlGtdACzRraile7rJ5mUkcRqLVDTGFMVOAr0AwYkK3MYuAVYYYwpC9QG9mdloCIiIgAkxMPPY5wbKFdoBn2mQmDQlVeXYJm19ghjv9tBTGwCozvW5IEO1fHz0ajpknnpJlbW2jhjzMPAUpzhFiZaa7cZY0a5po8DXgMmG2O24Jw6fMZaG5aNcYuISH4UfQa+vhf2/QRNh8Jtb1/VVX97T0Xyn3mbWXvwLK2qleD1ng2oXrpwFgYs+U2GxrGy1i4BliR7bZzb42NA56wNTURExM2JrTB7IEQchTs+gODhV1xVXHwCn684wPs/7qagrzdv9WpI72ZBGuRTrppGXhcRkdxvy1xnFHX/QBi+BCq2uOKqdp44x9NzN7M5JIKu9crxao96lCnin4XBSn6mxEpERHKvhHj44UVY9QlUag29p0CRsldUVWx8Ap/+so9PftlDUX9f/m9AU25rUE6tVJKllFiJiEjuFHMO5o6AvT9A83uhyxvgU+CKqtp6NIKn5m5mx/FzdG9Unpe61dUQCpItlFiJiEjuc/YgzOgHYbvhjvcheMQVVZPYSvXxz3soHlCAzwY3o3O9clkbq4gbJVYiIpK7HF4NswZAQiwM+hqq33RF1ewLjeTx2Zv4KySCOxuX55Xu9ShW6MpavEQySomViIjkHpu/gm8fcsalGvAVlKqZ6SoSEixTVx3kje92UrCAN58MaMIdDctnfawiKVBiJSIinpeQAL+8DivegSptnUE/C2X+ZsfHwi/w1Ny/+GPvaW6qXZo3725ImaK64k9yjhIrERHxrNgY+OZ+2D4fmgyG29+7ok7q8zce5YVvtxKfYHnjrgb0a15RV/xJjlNiJSIinhN9xulPdXgVdHoNbngEMpkMnYuJ5cX5W5m/6RjBlYvzbp9GVC4ZkE0Bi6RNiZWIiHhG+GGY3gvOHoBeE6H+3ZmuYv2hs4yevZFj4TE83qkWD91UA28vtVKJ5yixEhGRnHdii5NUxV6AQfOgattMvT0+wfLpL3v54Kc9XBfoz1f3t6ZZ5eLZFKxIximxEhGRnLXvF5g9GPyLwojvoWzdTL39aPgFHpu1iTUHz3Bn4/K81qM+Rf19sylYkcxRYiUiIjnnr9nw7YNQqjYMnAOBFTL19u+2HOeZrzeTYOH9vo3o2SQomwIVuTJKrEREJGf88aFz378qbaHfl84NlTMoJjae/y7ZwdRVh2hUsRgf9WusDuqSKymxEhGR7GUt/PiSk1jVuwt6jgOfjN+n72BYFA/N2MC2Y+e4r21VnupyPQV8vLIxYJErp8RKRESyT0I8LHoMNkyB4HvgtnfAK+NJ0aLNx/j311vw9jJ8MSSYjnXLZmOwIldPiZWIiGSPuEsw7z5n4M+2T8LNz2d4jKqY2HjGLN7O9D8P06RSMT4Z0JQKxQpmc8AiV0+JlYiIZL1LUc6Vf/t+gs5jnIE/M+hgWBQPfrmB7cfPcX+7ajzZpTa+3jr1J9cGJVYiIpK1LpyFGX0hZC10/wSaDs7wW3/cfpLHvtqEt5dh4rBgbr5ep/7k2qLESkREsk7kKZjWE8J2Q+/JUPfODL0tPsHywY+7+fjnvTSoEMj/BjUlqHih7I1VJBsosRIRkaxx7hhM6eb8HzAbqt+cobedjbrEo7M2smJPGH2Cg3j1zvr4+3pnc7Ai2UOJlYiIXL3wI05SFRUGg7+BSq0y9LYtIRGMmr6e0PMXeeOuBvRvUSmbAxXJXkqsRETk6pw5AFO6Q0wEDJkPQcEZettXa4/w/LdbKRVQgDmjWtOoYrFsDlQk+ymxEhGRK3d6n9NSFRsNQ7+F8k3SfcvFuHheXrCdmWsOc2ONknzUrwklC2d8wFCR3EyJlYiIXJnQXU5LVUIsDF0I5Rqk+5aj4Rd4cPp6/gqJ4IEO1XmiUy18NJSC5CFKrEREJPNOboep3QEDwxZDmTrpvmXl3jAenrmRS3EJjBvUjK71y2V/nCI5TImViIhkzoktTkuVj5/TUlWqZprFrbVMWXmQ1xbvoGqpAMYPbkb10oVzKFiRnKXESkREMu7EViep8i0IwxZBiWppFr8Ul8BLC7Yyc80ROtYpw/t9G1PE3zeHghXJeUqsREQkYxJP//n4ZyipOh15kQe+3MCaA2d4oEN1nuxcG2+vjN0rUORapcRKRETSd2qnk1R5+WYoqdpx/Bz3TV3HqfMX+aBvY3o0qZBDgYp4lhIrERFJW+huZ0gF4+UkVSWrp1l86bYTPDZ7E4X9fJhzv8ankvxFiZWIiKQubK+TVGFh6KI0O6pba/nk5728+8NuGgUF8tmQYMoW9c+5WEVyASVWIiKSstP7YModkBDntFSVrp1q0QuX4nlq7l8s2nycOxuX5827G+p+f5IvKbESEZF/OnvQaamKu+gkVWmMU3U84gIjp65n67EInul6PaPaV8MYdVKX/EmJlYiIJHXuGEy9Ey5FOeNUla2XatEtIRGMmLKW6ItxfD44mI51y+ZgoCK5jxIrERH5W1SYk1RFhcGQBXBdw1SLLtt2gn/N2kSJgAJMf/BGapcrkoOBiuROSqxERMRxIRym9YDwwzDoawhqlmIxay0T/zjImMXbaVghkM+HBlOmiDqpi4ASKxERAbgYCV/2csar6j8LqrRJsVhcfAKvLtrO1FWH6FKvLB/0bULBAuqkLpJIiZWISH4XewFm9oOjG6DPFKjZMcVikRfjeGTGBn7ZFcrIdtX4d9fr8dJI6iJJKLESEcnP4i7BV0Ph4O/QczzU6ZZiseMRFxgxeR27T55nTI/6DGpVOYcDFbk2KLESEcmvEuJh3n2wZync8T406ptisa1HI7hnylqiLsYzcVhz2tcqncOBilw7lFiJiORH1sLiJ2D7fOj0GgSPSLHYTztO8sjMjRQr6MvcB1pzfbmiORyoyLVFiZWISH70y+uwfhK0eQxufDTFIlNWHuSVhduoVz6QCUODKaPb04ikS4mViEh+8+c4WP42NBkMt7z0j8kJCZax3+/ks+X76VS3LB/2a0yhAvq6EMkI7SkiIvnJ5jnw/TNw/R1wxweQ7NYzF+PieWrOZhb8dYwhrSvzUrd6eOvKP5EMU2IlIpJf7PkB5o+CKm3h7gngnfQr4FxMLKOmrWflvtO655/IFVJiJSKSHxxZA7MHQ5m60G8G+CbtL3XyXAxDJ65h76lI3uvTiLuaBnkoUJFrmxIrEZG87tQO+LI3FL3OuVWNf9Ir+/aeOs/QiWsJj77EpOHNaVtTwymIXCklViIieVn4YZh2F/j4w+BvoHCZJJPXHjzDvVPWUcDHi9n3t6Z+hUAPBSqSNyixEhHJq6LCYFpPiI2C4d9B8SpJJn+/9TiPztpEULGCTBnRgoolCnkmTpE8RImViEhedPG8c1PliBAY8i2UrZdk8pSVB3l54TaaVCzGF0ObUyKggIcCFclblFiJiOQ1cRdh1kA4vhn6z4RKrS5Pstby1tJd/O/XfXSsU5aP+zehYAFvDwYrkrcosRIRyUsSEmD+A3DgN+gxDmp1uTwpLj6Bf8/bwtz1IQxoWYlXu9fDx9vLg8GK5D1KrERE8pIfX4KtX0PHl6Fx/8svx8TG89CXG/hp5ylGd6zJv26pqTGqRLJBhn6qGGO6GmN2GWP2GmP+nUqZDsaYTcaYbcaY37I2TBERSdfq8bDyI2h+L9w4+vLLERdiGTxhNT/vOsVrPeozumMtJVUi2STdFitjjDfwf0AnIARYa4xZYK3d7lamGPAp0NVae9gYUybl2kREJFvsWAjfPQO1b4db37p8q5pT52IYMnEN+0Ij+bh/E+5oWN7DgYrkbRk5FdgC2Gut3Q9gjJkF3AlsdyszAJhnrT0MYK09ldWBiohIKg6vhq/vhaBguPsL8HI6ox8Mi2LwxNWcjrzExGEa+FMkJ2TkVGAF4Ijb8xDXa+5qAcWNMb8aY9YbY4akVJExZqQxZp0xZl1oaOiVRSwiIn8L2wMz+0LRCtB/NhRwxqLadiyCXuNWERkTx8z7WimpEskhGUmsUjoRb5M99wGaAbcDXYAXjDG1/vEmaz+z1gZba4NLl9ZOLiJyVc6fhOl3gZePc6uagJIA/Ln/NP3G/0kBb8OcUTfQqGIxDwcqkn9k5FRgCFDR7XkQcCyFMmHW2iggyhizHGgE7M6SKEVEJKmLkTCjtzO6+rBFUKIqAEu3neCRmRupVKIQU0e0oHyxgh4OVCR/yUiL1VqgpjGmqjGmANAPWJCszLdAW2OMjzGmENAS2JG1oYqICADxsTBnKJzYCr2nQIVmAHy19ggPTF9P3euKMuf+1kqqRDwg3RYra22cMeZhYCngDUy01m4zxoxyTR9nrd1hjPke2AwkAF9Ya7dmZ+AiIvmStbBoNOz9Ebp9BLU6AzD+t3288d1O2tYsxbhBzQjw0zCFIp6QoT3PWrsEWJLstXHJnr8NvJ11oYmIyD/8OhY2Tof2z0CzoVhreXvpLj79dR93NLyO9/o0poCPRlMX8RT9pBERuVZsnA6/jYXGA6HDf0hIsLy8cBtTVx2if4uKjOnRAG8vDfwp4klKrERErgX7f4OF/4JqN0G3D4lLsDw9dzPzNh5lZLtq/OfW6zWaukguoMRKRCS3C90FswdDyZrQZwoxCV48OnMDy7af5MnOtXjophpKqkRyCSVWIiK5WWQofNkbfArAgNlEmQBGTlnLH3tP83K3ugy7saqnIxQRN0qsRERyq9gYmDUAIk/CsCVE+JVn+ITVbDoSzru9G3F3syBPRygiySixEhHJjRISYP4DELIG+kwlNLA+gz9bxf7QKD4d2Iyu9ct5OkIRSYESKxGR3OiXMbBtHnR8haPlOzNo/CpORMQwYViw7vsnkospsRIRyW02TocV70LToeyvdQ+D/reS8xfjmHZPC4KrlPB0dCKSBiVWIiK5iduwCtuavsjQz/7EWpg1shX1ygd6OjoRSYeG5xURyS1Cd8NXg6FkDTa2/oB+X6zH19uLr0a1VlIlco1QYiUikhtEhcGXvcC7AGtvGM+AqTspGVCAOaNaU710YU9HJyIZpFOBIiKeFhsDM/tD5ElWtZ3K0LnHqVY6gKn3tKBMEX9PRycimaDESkTEk9yGVfiz2XsMWhpPw6BAJg9rQWAhX09HJyKZpMRKRMSTfnkdts1jfa3R9PujHDdUL8HnQ4IJ8NPhWeRapD1XRMRTNn4JK95hW7ke3L25OZ3qluXj/k3w9/X2dGQicoWUWImIeMKB5diF/+JQYAvuPHgXPRpX4O3ejfD11jVFItcyJVYiIjktdDd29iBCCwTR/eR99G5Rjdd71MfLy3g6MhG5SkqsRERyUlQYdkZvImO9uCt6NH3a1Oe52+tgjJIqkbxAiZWISE6JjSFhZn/iwo8xOOZ57r75BkZ3rKmkSiQPUWIlIpITEhKI/+YBvEPW8K9L/+K2W+9gZLvqno5KRLKYEisRkRxw6ccxFNg+j7Gx/bih+z0MblXZ0yGJSDZQYiUiks2i10yj0Mp3mR3fgZo9n+fu4IqeDklEsokSKxGRbHRuxy8UXDKaPxLqE9jrY7o2UlIlkpdpwBQRkWwSdnArzB7EIVuWhN5T6NqokqdDEpFsphYrEZFscOzYERKm3I2xXpzrOYO2DWp4OiQRyQFqsRIRyWIHjocR+nkvStnThHWbTNPGjT0dkojkECVWIiJZaOfxCHZ9NpRGdiehHT+idvAtng5JRHKQEisRkSzy15Fwfhv/GF3t74S1+g8V2wzwdEgiksOUWImIZIHV+08z+4u3uJ+viaw7gFJdnvF0SCLiAeq8LiJylX7bHcqEaVOY4D2eixXbUvjuj0C3qRHJl5RYiYhchaXbTvDejMXMLfA+png1/AZMB29fT4clIh6ixEpE5ArN33iUMXNWsLDg2wT4+eM1eA4ULObpsETEg5RYiYhcgRmrD/PK/PV8W+RDyiWcxfRfBMWreDosEfEwJVYiIpn0xYr9vL54GzNLTOb66O3QezJUbO7psEQkF1BiJSKSQdZaPv55L+/9sJtPy39PqzO/QseXoV5PD0cmIrmFEisRkQyw1jL2u52MX76fsdU2c9ux6dB0CNw42tOhiUguosRKRCQdCQmWFxdsZfqfh3mxXhh9D7wD1W6C29/TsAoikoQSKxGRNMTFJ/D015uZt+Eoz7YwDN/1PKZkTegzRcMqiMg/KLESEUnFpbgE/jVrI99tPcHz7Utyz677MD7+MPAr8A/0dHgikgspsRIRSUFMbDyjpq/n112hvNS1KsP3PAKRoTB8MRSr5OnwRCSXUmIlIpLM+ZhY7pu6jtUHzvBGz3r0P/gCHF0PfadDhWaeDk9EcjElViIibs5EXWLYpDVsO3aOD/o25s7Q8bBjAXT5L9S5w9PhiUgup8RKRMTlREQMgyas5siZaD4b3IxbopbAHx9C83uh1YOeDk9ErgFKrEREgINhUQyasJrw6FimjGhBq4SNsPgJqNEJur6pYRVEJEOUWIlIvrfj+DkGT1hDfEICM+5rSUPfozBhGJSpC70ngbcOlSKSMV6eDkBExJPWHzpL3/Gr8PEyfHV/axoWvQBf9gG/wjBgNvgV8XSIInIN0c8wEcm3VuwJZeTU9ZQp6sf0e1pSsbCFSX3hwlkY8R0EVvB0iCJyjVFiJSL50vdbj/PozE1UKx3A1HtaUCbAF2YPghNboN9MuK6Rp0MUkWuQEisRyXfmrDvCM19vpnHFYkwa1oLAgj7w3dOwawnc+jbU7urpEEXkGqXESkTylQm/H+C1RdtpU6MU4wc3I8DPB/74CNZ8Bq0fhpYjPR2iiFzDlFiJSL5greX9H/fw0U976FqvHB/2b4yfjzdsmQs/vAD1ekKn1zwdpohc45RYiUiel5BgeXXRdiavPEivZkGMvasBPt5ecGAFzH8AKt8IPcaBly6UFpGro8RKRPK0uPgEnv56M/M2HGXEjVV5/vY6eHkZOLUDZg2E4lWh35fg6+/pUEUkD8jQzzNjTFdjzC5jzF5jzL/TKNfcGBNvjOmVdSGKiFyZmNh4Rk3fwLwNR3msYy1euMOVVJ07BtN7gW9BGDQXChb3dKgikkek22JljPEG/g/oBIQAa40xC6y121Mo9yawNDsCFRHJjIjoWO6dupZ1h87ySvd6DL2hijMh5hx82RtiwmH4d1CskkfjFJG8JSMtVi2Avdba/dbaS8As4M4Uyj0CfA2cysL4REQy7eS5GPqMX8WmI+F81K/J30lV3CX4ajCE7oQ+U+G6hh6NU0TynowkVhWAI27PQ1yvXWaMqQD0BMalVZExZqQxZp0xZl1oaGhmYxURSdf+0Eju+nQlIWejmTSsBd0alXcmWAsLH4X9v0L3j6HGLR6NU0TypowkVind0t0me/4B8Iy1Nj6tiqy1n1lrg621waVLl85ojCIiGbI5JJxe41YRExvPzJGtaFOz1N8Tf3kd/poJNz0PjQd4LkgRydMyclVgCFDR7XkQcCxZmWBgljEGoBRwmzEmzlo7P0uiFBFJx+97wrh/2jqKFSrAtHtaUK104b8nrpsEy9+GpkOh3ZOeC1JE8ryMJFZrgZrGmKrAUaAfkOTnnrW2auJjY8xkYJGSKhHJKYs2H+Ox2ZuoXrowU0a0oGxRt6ETdn0Pix+Hmp3h9vfApNQILyKSNdJNrKy1ccaYh3Gu9vMGJlprtxljRrmmp9mvSkQkO01ddZCXFmwjuHJxvhjanMCCvn9PPLQK5gx1bqjcaxJ4a+g+EcleGTrKWGuXAEuSvZZiQmWtHXb1YYmIpM39FjUd65TlkwFN8Pf1/rvAyW0wsy8EVoSBc8GvcOqViYhkEf18E5FrTlx8Ai98u42Zaw7TJziI//Z03aIm0dlDMO0u8A2AwfMgoFTqlYmIZCElViJyTYm+FMcjMzby085TPNihOk91qY1x7zcVGQrTekJcDIz4XgOAikiOUmIlIteMsMiL3DN5LVuORvBaj/oMblU5aYGL5+HLXs4ta4Z8C2XqeCZQEcm3lFiJyDXhQFgUQyeu4dT5GMYPDqZT3bJJC8RddG6qfGIL9J8FlVp6JlARydeUWIlIrrfh8FnumbwWYwwz72tFk0rJbpqcEA/z7oMDv0HP8VCrs2cCFZF8T4mViORqy7ad4JGZGykX6M+U4S2oUiogaQFrYcmTsP1b6PJfaNTPM4GKiKDESkRysWmrnDGqGgQVY8LQYEoV9vtnoV/HwrqJ0OYxaP1QjscoIuJOiZWI5DoJCZa3lu5i3G/76FinDB/1b0KhAikcrtZ8Dr+NhSaD4JaXcj5QEZFklFiJSK5yMS6ep+du5ttNxxjYshKvdK+XdIyqRFu/hiVPQe3b4I4PdasaEckVlFiJSK4RER3LqOnrWbX/NE93rc0D7asnHaMq0e6lMG8kVGoNvSbqVjUikmvoaCQiucLBsChGTFnLkTPRvN+3ET2bBKVc8MBymD0YyjWAAbPBt2DOBioikgYlViLicWsOnOH+aeuwwPR7WtKyWsmUCx5ZCzP6QYlqMGge+BfN0ThFRNKjxEpEPOqbjSE8M3cLQcULMnFY838Op5DoxBb48m4oUhaGzIdCJXI2UBGRDFBiJSIeYa3l/R9289HPe2ldrST/G9SUYoUKpFw4bI9z/78ChZ1b1RQpl7PBiohkkBIrEclxMbHOlX8L/jpGn+AgxvRoQAGfFK78Azh7CKbe6TweskA3VRaRXE2JlYjkqLDIi4ycuo4Nh8N5puv1jGpfLeUr/wDOn3CSqkuRMGwJlKqRs8GKiGSSEisRyTF7Tp5n+OS1hEVe5H8Dm3Jrg+tSLxx5CqZ0g6hQ5/Rfufo5F6iIyBVSYiUiOWLFnlAenL4B/wLezB7ZmkYVi6VeOCoMpnSHiBAY9DUEBedcoCIiV0GJlYhkK2stU1Ye5LXFO6hZpjAThjWnQrE0xp6KPuOc/jt7EAZ+BZVvyLFYRUSulhIrEck2F+PieXH+NmavO0LHOmX5oF9jCvulcdi5cNZJqk7vhf6zoGq7nAtWRCQLKLESkWxx6nwMo6atZ8PhcB69uQajO9bCyyuN+/ldCHeGVAjdyf+3d9/xWZX3/8dfVwIZJCGDDDIhYe89RIsCDhQV66ggrg791mqtrR1f7bDt79uHdtcObXEVK0oRaUErUkerCJUlIBtCgCRk7z3v6/fHSTBFkEDukdx5Px+PPO517nM+PC6Jb865zudi8UswZI73ihURcRMFKxFxu49yK7j7+e1U1jfzxJLJXPVpk9QBGqrghRugYA8sWg5DL/VOoSIibqZgJSJu9bcdufzvK7uJDQ/mlXtmMTrpLMvONFbD8hshfyd87nkYfoV3ChUR8QAFK9+bOAQAABvhSURBVBFxi1aX5advHGDpe1nMSI/hiSWTGRAe/OlfaqiEF26EE9vhpudg5ALvFCsi4iEKViLSZZV1zdy/YgfvHirm9gsG8f2rR9M38Ayd1NvVV8AL10P+LrjpzzD6Wq/UKiLiSQpWItIlmUXV3PX8dnLL63j0+nEsnt6JJWfqyuAv10HhPvjcX2DkVZ4vVETECxSsROS8vbmvkK//dSchfQN48a6ZTBscc/Yv1ZbA89dBySFY9CIMv9zzhYqIeImClYics1aX5VdvHuQP/zrCuORI/njblE9v+tmuphievxbKspyWCkPneb5YEREvUrASkXNSVtvE/S/t4P3MEhZPT+WRa8YQ0jfw7F+sLmhbpiYHblkJGRd7vlgRES9TsBKRTtuZU8FXXthOSW0TP71hHDdP68R8KnDW/Ht+IVTlw5JVMPhCzxYqIuIjClYiclbWWl7aksMP1+4lLiKYV748i3EpkZ37ckmmM1G9oRJu+xukzfBssSIiPqRgJSKfqr6plR+s2cPL23OZPTyOx2+eSHRYUOe+XLDbWabGWrjzNUic4NliRUR8TMFKRM4os6iGe5d/yMHCar7att5f4Ket99dR9mZ48SYIioDb/w6xwzxbrIhIN6BgJSKntWbnCR5avZuQvoEs+8J0Lh4e1/kvH3kHViyBiES4fQ1EpXquUBGRbkTBSkT+S0NzKz96dS8vbclh2uBofrt4EomRnWil0G7fWnjlixA7Am5bDeHxnitWRKSbUbASkZOyimv4yvIPOVBQzT2XDOHBy4bT52xL03S0YzmsvQ+Sp8KSlRAa7bliRUS6IQUrEQGcS38Pr95NUJ8Anvv8NOaMOIczTdbC+7+Gt38EGXNg0XIICvNcsSIi3ZSClUgvV9vYwiNr97Jqey5TBkXzu8WTSOpMF/V2rlZY9x3Y+hSMuwkWPgF9OnnXoIiIn1GwEunFdudWcv+KHRwrreWrc4dy/7xh9D2XS3/NDbD6S7D/VZh1P1z6Iwg4h++LiPgZBSuRXsjlsjz9fhY/X3+Q2PBgXrprJjMzBpzbTurL4aVbIPs/MP8xmHmPZ4oVEelBFKxEepmi6gYeXLmLDYdLuGJMAj+9YTxR/c7x0l1FDiy/0VlM+cZnYez1nilWRKSHUbAS6UXeOVDIt17+iNqmFn7y2bHcMj0NYzrZ8LNd4V544UZoqoVbV0P6ZzxTrIhID6RgJdIL1Da28JPX9/Pi5mxGDoxgxeKZDEuIOPcdHX4LXr4TgsPhC+sgYYzbaxUR6ckUrET83Pbj5Xxj5U6yy+r4n9kZfOPy4QT3CTz3HW1eCm98xwlTi/8KkcnuL1ZEpIdTsBLxU00tLn779mGe+HcmiZGhrLhrJjPOdYI6QGsLrH8ItiyFEVfB9U85Z6xEROQTFKxE/NDhwmq+vnIne05UceOUFB65ZjQRIX3PfUcNVbDqC5D5JlxwH1z2Ywg4j7NdIiK9hIKViB9pdVme23iUn60/SHhwH/546xTmjx14fjuryIYXb4aSQ3DN4zDlTrfWKiLijxSsRPzEkeIavvXyLj7MruDSUfE8ev144iKCz29nOVthxWJoaYJbX4GMS9xZqoiI31KwEunhWl2Wpzdk8cs3DxHaN5Df3DyRhROTzr2NQrvty+AfD0JkCty5EuKGu7dgERE/pmAl0oMdLqzmm6s+YldOBZePTuD/PjuW+IiQ89tZS5Nz19+2Z2HIXLjhGegX496CRUT8nIKVSA/U0upi6YYsfvPmYcKCA/nt4klcMz7x/M9SVRfAytshZzNc+ADM+4EmqYuInAcFK5EeZldOBQ+t3s2+/CquHDuQHy8ce/5zqcCZT7XyNmiohBuf0/I0IiJdoGAl0kPUNLbwy38eZNmmY8SGB/PkkslcOS6xazv98HlnPlVEInzxTRg41j3Fioj0Up0KVsaY+cDjQCDwtLX2sVM+XwJ8p+1lDXCPtXaXOwsV6c3e2lfID9bsIb+qgSUz0vj2/JH0P5++VO2a6+H1b8GOv0DGHGchZc2nEhHpsrMGK2NMIPAH4DIgF9hqjFlrrd3XYbOjwMXW2nJjzJXAUmCGJwoW6U2Kqhr44at7eX13AcMTwll1ywVMGdTFAFSSCS/fAYV74DMPwiUPQ6BOXouIuENnfptOBzKttVkAxpgVwELgZLCy1m7qsP0HQIo7ixTpbVpaXbzwwXF++eYhGltcfPPy4dw9ewhBfQK6tuPdq+DVr0FgECx5BYZd6p6CRUQE6FywSgZyOrzO5dPPRn0RWNeVokR6s+3Hy/je3/eyP7+Ki4bG8v+uG0t6bFjXdtrcAOsfhm3PQOoM59JfpP79IyLibp0JVqe7f9uedkNj5uAEq4vO8PndwN0AaWlpnSxRpHcoqWnksXUHWLU9l8TIEJ5YMpkrxw48/xYK7cqy4OU7IX8XzLrfaaUQ2IX5WSIickadCVa5QGqH1ylA3qkbGWPGA08DV1prS0+3I2vtUpz5V0ydOvW04Uykt2l1WZZvPs4v1h+kvrmVey4Zwn1zhhIW7IZ5Tx+9DP/4BpgAWLwCRlzZ9X2KiMgZdeY391ZgmDEmHTgBLAJu6biBMSYNWA3cZq095PYqRfzU5qxSfvzaPvbmOZf9fnjtGIbGh3d9xw2VThuF3S87l/6ufwqiB3V9vyIi8qnOGqystS3GmPuA9TjtFp611u41xny57fM/Aj8ABgBPtF22aLHWTvVc2SI9W3ZpHY+u28+6PQUkRYbw+1smsWBcFzqnd3R8E6z+H6g6AXO+Cxd9Q3f9iYh4ibHWN1fkpk6dardt2+aTY4v4SnVDM0/8+wjPbDhKYIDhK5cM4a7ZGYT0dcPyMa3N8O9H4f1fQ9QguOFpSNG/b0RE3MEYs70zJ430z1gRL2h1WVZtz+Hn6w9RUtPI9ZOT+fYVIxkYeZ4LJp+q9Ai88iXI+xAm3QrzH4PgCPfsW0REOk3BSsSDrLW8d7iEx9YdYH9+FVMGRfPMHVOZkBrlngO4XLBlKbz1Q+gTDJ97HkYvdM++RUTknClYiXjIR7kVPLbuAJuOlJIaE8pvF0/imvFumkcFzlmqNfdB9iYYdjlc8zj0T3LPvkVE5LwoWIm42fHSWn6+/iCvfZRPTFgQj1wzmiUzBnW9a3o7Vyts/hO8/WOng/p1T8KExeCuwCYiIudNwUrETUpqGvnd24dZvjmbvoEBfHXuUO6enUFEVxZL/sRBMmHNvZDzAQyfD1f/WmepRES6EQUrkS4qr23iT+9lsWzTMZpaXdw8LZUH5g0jvr+bJqaDc8ffB0/Cv34CfULgs0th/Od0lkpEpJtRsBI5T5X1zTyzIYtnNx6jtqmFayck8bV5w8iIc0ODz45ytsJrD0DhHhixAK7+FUQMdO8xRETELRSsRM5RTWMLz71/lKc2ZFHV0MJV4wbywKXDGZ7g5vYG9eXw1o9g+5+dy303L4eRC3SWSkSkG1OwEumkqoZm/vKf4zy9IYvyumYuHZXA1y8bxpikSPceyFpnKZr1D0NdGVxwL1zyEAS7+UyYiIi4nYKVyFmU1Tbx3Maj/HnTMaobWrhkRBxfv3S4+3pRdVRy2Fnj7+i7kDwFbl0NiePdfxwREfEIBSuRMyiqauCpDVks35xNXVMr88cM5L65Qxmb7OYzVAD1FfDuz2DLn6BvGCz4JUz5PAS4YakbERHxGgUrkVNkl9bx1IYs/roth5ZWF9dOSOIrc4a6fw4VOD2pPlwG7/yfc9lv8u0w9/sQHuf+Y4mIiMcpWIm02ZVTwdL3sli3J5/AAMMNk1P48sVDGBwb5pkDHn0P3njIudtv0IUw/1FInOCZY4mIiFcoWEmv5nJZ/nWwiD+9l8WWo2VEhPTh7tlDuHPWYPctkHyq0iPw1iOw/1WITIObljnr++luPxGRHk/BSnqlhuZW1uw8wVMbjpJZVENSZAjfWzCKRdPTCA/20F+L6kJ496fOpb/AIJjzPZh1H/QN9czxRETE6xSspFc5UVHPCx8cZ8WWbMrrmhmd2J/HF03kqnGJ9A1001p+p2qohI2/hQ+egNYmmHwHXPxtNfkUEfFDClbi96y1bD5axrJNx1i/twCAy0cP5I5Zg5mZEYPx1CW45gbY+jRs+IXT7HPsDTDnuzBgiGeOJyIiPqdgJX6rtrGFtbvyWLbpGAcKqonq15e7Zw/h1plppET389yBWxphxwuw4VdQlQtD5sK8RyBpoueOKSIi3YKClfidvXmVvLg5mzU786hpbGHkwAgeu34cCycmExrkwb5QzQ2w4y/w/q+h6gSkTIPrnoCMiz13TBER6VYUrMQv1DW18NqufJZvyWZXTgXBfQJYMD6RJTPSmJwW7bnLfeAEqg+Xwfu/geo8SJ0JC38PGXN0p5+ISC+jYCU9lrWWXbmVvLwth7U786hubGFofDg/uHo0109OJqpfkGcLaKyBD5+HjY9DTQGkzYLPPgnpFytQiYj0UgpW0uMUVTXwtx0nWLU9l8NFNQT3CeCqcYksnp7GtMEePjsFUFMEm//kTExvqIBBF8ENT8HgzyhQiYj0cgpW0iM0trTy9v4iVm3P5d1DxbS6LFMGRfPo9eNYMD6R/iF9PV9EyWHY9DvYtcJpmzByAVz4NUid7vlji4hIj6BgJd1Wq8uyOauUNTvzeH1PPtUNLST0D+bu2RncOCWFIXHhni/CWsj+D2z6PRx83WnsOfEWuOA+iB3q+eOLiEiPomAl3Yq1lr15Vfx9xwle/SiPwqpGwoICuWLMQBZOSuaiobEEBnjhcltTHex+GbY8BYW7ITQaZn8Lpt+tBZJFROSMFKzE56y17Muv4o09Bfxjdz5ZxbX0DTRcPDye7y1I4tJRCZ5tk9BRWRZsfcZpm9BQCfFj4JrHYdxNEOShxZhFRMRvKFiJT1hr2ZlTwRt7Cli3p4DssjoCDExPj+FLF2Vw1biBnr+rr11rCxx52wlUh/8JJgBGX+ucnUq7QBPSRUSk0xSsxGtaXZbtx8tZtyef9XsKyKtsoE+AYdbQWO65ZAiXjU4gNjzYewWVHXU6pO9cDtX5EBbvrOE35U7on+S9OkRExG8oWIlHNTS38kFWKW/tL+SNPYWU1DQS1CeA2cPiePDyEVw6KoHIfl64o69dcz3sf81p6Hlsg3N2auhlcNXPYfh8CPRiLSIi4ncUrMTt8ivreedAEf86UMTGzFLqm1sJ7RvInJFxzB+byNyR8YQHe/E/PZcLcj6Aj1bC3tXO3KmoQTD3ezDhFohM9l4tIiLi1xSspMtaXZYd2eW8c6CIdw4UcaCgGoCU6FBumprCnJHxXJAxgJC+XpqA3q5ovxOmdq+Cymzo28/pPTXpNqeZZ0CAd+sRERG/p2Al56WoqoGNR0r498Fi3j1UTEVdM4EBhqmDonnoypHMHRnP0Phwz3dBP1VFNuz9O+xeCQW7wQTCkDkw7/sw4ioI9kLvKxER6bUUrKRTqhqa2ZxVxsbMEjZmlnC4qAaAmLAg5o6IZ+6oeD4zLI7IUB/MUSo9AvvWwP61kLfDeS95Klz5MxjzWQiP935NIiLSKylYyWk1trSy/Xg5mzJLeT+zhI9yK3BZCOkbwPT0Adw4JYULh8YyOrE/Ad5o2NmRtVB80AlS+9Y6DTwBkqfAZT+GUddCTLp3axIREUHBStrUN7WyI6ecrUfL2XqsjG3Hy2hodhEYYBifEsm9c4Yya0gskwdFEdzHy3OlAFqaIHsTHHwDDq2D8mPO+6kz4YpHYdQ1EJXq/bpEREQ6ULDqpcprm9oCVDlbjpax50QlLS6LMTAiIYJF09K4cGgsMzJivLPA8enUlkLmm3BwHRx5BxqrIDAY0mfDrK86c6bUb0pERLoRBatewFrLsdI6duaUs/VYOVuPlp2cIxUUGMD4lEjump3BtMHRTEmL8W5fqY5amiB3ixOijrwDeTsBC+EJMOY6p89UxiVaWkZERLotBSs/VFLTyK6cCnblVLAzt5JdORVU1jcDEBHch8mDorluUjLTBscwPiXS+20Q2rXPlTr6rhOkjm6A5lrnTr7U6TDnYRg6DxInqTWCiIj0CApWPVxtYwv78qvYlVPBjrYwlVteD0CAgeEJEVw5diATU6OYkBrF8IQIAr092bydywVFe+HYRjj+PhzfBHWlzmfR6TBhkROkBl8EIZG+qVFERKQLFKx6kOLqRvblV7E3r5K9eVXsy6viWGkt1jqfJ0eFMjE1itsvGMSElCjGpUTSL8iHQ9zaAgUfwfGNTpjK3uR0PQeITINhl8OgC2HwhRCT4bs6RURE3ETBqhtqdVmOl9ZyoKCafXkfB6mi6saT26REhzImqT/XTUxmTFJ/xqdGEh8R4sOqgdoSyN0GJ7bDiW2QsxWanC7sxAyB0QudIDVoFkSl+bZWERERD1Cw8iGXy5JbXs+hwmoOFVVzqKCaQ4U1ZBbX0NTiAiAwwDA0LpyLhsYyOqk/o5P6MyYx0ncTzNs110P+LidEtYepiuPOZyYA4sfA+M85Z6PSZkH/RN/WKyIi4gUKVl7Q1OIiu6yOrOIaskpqySyq4VBhNYcLa6hvbj25XVJkCMMSIrhoWCzD4sMZMTCC4QkRvptc3q6xBor2OUvEFOx2QlTRPnC1OJ9HpjrNOad9CVKmQuIE3bknIiK9koKVm1hrKa5pJKu4tu3HCVFHS2rJLquj1WVPbhsXEez0ipqeyvAEJzwNSwj3Xb+odtZCVZ4TngrbQlTBHijLAtrqD4mEpElw4QNOmEqeAhEJPi1bRESku1CwOgeNLa3klteTXVZHblkd2Sd/6skpq6OmseXktsF9AkiPDWNUYgQLxiWSERdGRlw46bFhvllPryNroaYISg467Q6KD0LxASjcA/XlH28XnQ4Dxzp36yWMdZ5HpoK3F1YWERHpIRSsOmh1WYqrGzlRUUdOWX2H4FRHTlkdBVUNJ+/AAyc8pcX0Iy2mHzPSYxg8oB8ZceFkxIWRFBnq/TX0TuVqhYpsKDnkhKeSg1B8yHlsvzsPICgC4kY4a+wNHOf8xI+GkP6+q11ERKQH6jXBylpLaW0T+RUN5FXWk19RT35lA3mVDSefF1Y10NLhkh3AwP4hpMX0Y9aQWNJi+pEaE3oyTMVFBGN8ffampckJT2VZUH7UeSzLgrKjznp6ruaPtw2Lg7iRMPZGJ0jFDnceIxJ1FkpERMQN/DZYZZfW8Zu3D5Ff0UB+pROcGtvutGsXFBhAYlQIiZEhzEiPaXseSnJUKKkx/UiJDvX9xPHWFqjOh8ocqMiBymyozHVCU9lR533b4c8VFA4x6RA/CkYucPpDtYeofjE++2OIiIj0Bn4brFqt5YMjpSRGhTI2OZLLxwwkMTKEpKhQkiJDSYwKYUBYkG/POLlaobYYqgvafvKg8kSHEJXjTCa3rf/9vbA4pw9U6nRn/lN0uhOgYtKdz3T2SURExCf8Nlilx4ax6aF5vjm4ywX1Zc6ZpuqCMz/WFH0yNJlA6J8MUalOM82oVIhMcSaNR6U5z/uG+ubPJSIiIp/Kb4OVW7lczt1ydSVOd/Ha4rbnpR2et/3UlUBd2ScDE0C/Ac58poiBkDDm4+ftj+EDITwBAjUsIiIiPVHv+T+4tdBcB/UV0FDR4bH8NO+d8lhX+t/zmDoKiXQuv/WLhQFDnMtzYXHOT3to6p/oBKY+wd79M4uIiIhX+W+wKj4Ea77y3yGp4x1yn2CckBQaBSFRzmNkivO834C2sBT7yeeBPu5JJSIiIt1Gp4KVMWY+8DgQCDxtrX3slM9N2+dXAXXAndbaD91c67npE+wsq9I/+b/DUkgUhEZ/8r3g/hAQ4NOSRUREpGc7a7AyxgQCfwAuA3KBrcaYtdbafR02uxIY1vYzA3iy7dF3ogfB7Wt8WoKIiIj0Lp05RTMdyLTWZllrm4AVwMJTtlkIPG8dHwBRxphEN9cqIiIi0q11JlglAzkdXue2vXeu24iIiIj4tc4Eq9N1m7TnsQ3GmLuNMduMMduKi4s7U5+IiIhIj9GZYJULpHZ4nQLkncc2WGuXWmunWmunxsXFnWutIiIiIt1aZ4LVVmCYMSbdGBMELALWnrLNWuB245gJVFpr891cq4iIiEi3dta7Aq21LcaY+4D1OO0WnrXW7jXGfLnt8z8Cr+O0WsjEabfwec+VLCIiItI9daqPlbX2dZzw1PG9P3Z4boF73VuaiIiISM+ijpgiIiIibqJgJSIiIuImClYiIiIibqJgJSIiIuImClYiIiIibqJgJSIiIuImxumU4IMDG1MMHPfJwT0vFijxdRHiMRpf/6Wx9V8aW//mjfEdZK0967IxPgtW/swYs81aO9XXdYhnaHz9l8bWf2ls/Vt3Gl9dChQRERFxEwUrERERETdRsPKMpb4uQDxK4+u/NLb+S2Pr37rN+GqOlYiIiIib6IyViIiIiJsoWHWRMSbVGPMvY8x+Y8xeY8zX2t6PMca8aYw53PYY7eta5fwYYwKNMTuMMa+1vdbY+gFjTJQxZpUx5kDb398LNLb+wxjz9bbfyXuMMS8ZY0I0vj2TMeZZY0yRMWZPh/fOOJbGmIeMMZnGmIPGmCu8Xa+CVde1AA9aa0cBM4F7jTGjgf8F3rbWDgPebnstPdPXgP0dXmts/cPjwBvW2pHABJwx1tj6AWNMMnA/MNVaOxYIBBah8e2p/gzMP+W9045l2/9/FwFj2r7zhDEm0HulKlh1mbU231r7YdvzapxfzsnAQmBZ22bLgOt8U6F0hTEmBVgAPN3hbY1tD2eM6Q/MBp4BsNY2WWsr0Nj6kz5AqDGmD9APyEPj2yNZa98Dyk55+0xjuRBYYa1ttNYeBTKB6V4ptI2ClRsZYwYDk4DNQIK1Nh+c8AXE+64y6YLfAN8GXB3e09j2fBlAMfBc22Xep40xYWhs/YK19gTwCyAbyAcqrbX/ROPrT840lslAToftctve8xoFKzcxxoQDrwAPWGurfF2PdJ0x5mqgyFq73de1iNv1ASYDT1prJwG16LKQ32ibb7MQSAeSgDBjzK2+rUq8xJzmPa+2P1CwcgNjTF+cULXcWru67e1CY0xi2+eJQJGv6pPzdiFwrTHmGLACmGuMeQGNrT/IBXKttZvbXq/CCVoaW/9wKXDUWltsrW0GVgOz0Pj6kzONZS6Q2mG7FJzLwF6jYNVFxhiDM09jv7X2Vx0+Wgvc0fb8DmCNt2uTrrHWPmStTbHWDsaZDPmOtfZWNLY9nrW2AMgxxoxoe2sesA+Nrb/IBmYaY/q1/Y6ehzP/VePrP840lmuBRcaYYGNMOjAM2OLNwtQgtIuMMRcBG4DdfDwP52GceVYrgTScv+Q3WWtPnXwnPYQx5hLgm9baq40xA9DY9njGmIk4NyUEAVnA53H+samx9QPGmB8BN+Pcub0D+BIQjsa3xzHGvARcAsQChcAjwN85w1gaY74LfAFn7B+w1q7zar0KViIiIiLuoUuBIiIiIm6iYCUiIiLiJgpWIiIiIm6iYCUiIiLiJgpWIiIiIm6iYCUiIiLiJgpWIiIiIm6iYCUiIiLiJv8foblOLb3QKzYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 720x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(10,6))\n",
    "data_pred = np.linspace(start=15, stop=100, num=150)\n",
    "smoker_proba = log_smokers.predict_proba(data_pred.reshape([-1, 1]))\n",
    "plt.plot(data_pred, smoker_proba[:, 1], label=\"Smokers\")\n",
    "non_smoker_proba = log_non_smokers.predict_proba(data_pred.reshape([-1, 1]))\n",
    "plt.plot(data_pred, non_smoker_proba[:, 1], label=\"Non-smokers\")\n",
    "plt.title(\"Estimated probability of death per age, according to smoking habits\")\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}