{ "cells": [ { "cell_type": "markdown", "id": "560a0699-b258-4d85-a984-785ebfacdbd7", "metadata": {}, "source": [ "# Cox (semiparametric proportional hazard)" ] }, { "cell_type": "markdown", "id": "bf8ce873-7445-4f10-a1f9-aa9a44593851", "metadata": {}, "source": [ "This demonstration requires 2 optional dependencies : [pandas](https://pypi.org/project/pandas/) and [lifelines](https://pypi.org/project/lifelines/). Install them if you wish to reproduce examples outputs. Otherwise, below examples may help you to understand how you can use Cox model with ReLife." ] }, { "cell_type": "markdown", "id": "dc178d4a7653235a", "metadata": {}, "source": [ "### Loading data" ] }, { "cell_type": "markdown", "id": "49c5c2fdc80b3b2a", "metadata": {}, "source": [ "Let's use a [lifelines](https://pypi.org/project/lifelines/) dataset." ] }, { "cell_type": "code", "execution_count": 1, "id": "9345ce34-5550-4eb0-9362-9f08c607ca64", "metadata": {}, "outputs": [], "source": [ "import lifelines\n", "import pandas as pd" ] }, { "cell_type": "code", "execution_count": 2, "id": "c43a4b58a2972191", "metadata": {}, "outputs": [], "source": [ "df = lifelines.datasets.load_canadian_senators()" ] }, { "cell_type": "markdown", "id": "04142090-c0c2-44af-b47b-452a2e3e01a9", "metadata": {}, "source": [ "Clean up the dataset" ] }, { "cell_type": "code", "execution_count": 3, "id": "25d7e066-fec3-4d78-bd6c-d5432558ffcb", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
NameProvince / Territoryreasondiff_daysobserved
0Abbott, John Joseph CaldwellquebecDeath2363True
1Adams, MichaelnewbrunswickDeath1090True
2Adams, WillienorthwestterritoriesRetirement11766True
3Aikins, James CoxontarioResignation5333True
4Aikins, James CoxontarioDeath3133True
\n", "
" ], "text/plain": [ " Name Province / Territory reason diff_days \\\n", "0 Abbott, John Joseph Caldwell quebec Death 2363 \n", "1 Adams, Michael newbrunswick Death 1090 \n", "2 Adams, Willie northwestterritories Retirement 11766 \n", "3 Aikins, James Cox ontario Resignation 5333 \n", "4 Aikins, James Cox ontario Death 3133 \n", "\n", " observed \n", "0 True \n", "1 True \n", "2 True \n", "3 True \n", "4 True " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "df = df[[\"Name\", \"Province / Territory\", \"reason\", \"diff_days\", \"observed\"]]\n", "df = df[df[\"reason\"] != \"Appointment declined\"]\n", "# standardize text (removing white spaces, non alphanumeric characters and put in lower case)\n", "df[\"Province / Territory\"] = (\n", " df[\"Province / Territory\"].str.replace(\"\\W+\", \"\", regex=True).str.lower()\n", ")\n", "display(df.head())" ] }, { "cell_type": "code", "execution_count": 4, "id": "f3292125-7173-4977-acd8-a1f4923c8fec", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Province / Territory\n", "quebec 246\n", "ontario 242\n", "novascotia 98\n", "newbrunswick 93\n", "britishcolumbia 44\n", "manitoba 44\n", "alberta 43\n", "princeedwardisland 39\n", "saskatchewan 32\n", "newfoundlandandlabrador 30\n", "northwestterritories 7\n", "yukon 3\n", "westernprovincesdivision 2\n", "maritimesdivision 2\n", "quebecdivision 2\n", "ontariodivision 2\n", "nunavut 1\n", "Name: count, dtype: int64" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df[\"Province / Territory\"].value_counts()" ] }, { "cell_type": "markdown", "id": "141d1d5b-423a-4262-8817-aade4e74f319", "metadata": {}, "source": [ "Dummy encoding of \"Province / Territory\" covariates" ] }, { "cell_type": "code", "execution_count": 5, "id": "c003398a-387c-4e7b-9cce-9b0836f56045", "metadata": {}, "outputs": [], "source": [ "# dummy encoding of \"Province / Territory\"\n", "df = pd.get_dummies(df, columns=[\"Province / Territory\"], prefix=\"covar\", dtype=int)\n", "# remove covariate columns corresponding to low count province\n", "covar_cols = [col for col in df.columns if \"covar\" in col]\n", "province_count = df[covar_cols].sum()\n", "province_with_low_count = province_count[province_count < 10]\n", "df = df.drop(columns=province_with_low_count.index)" ] }, { "cell_type": "code", "execution_count": 6, "id": "c936f4a2-33df-4c94-9546-8291c46eee47", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Index(['Name', 'reason', 'diff_days', 'observed', 'covar_alberta',\n", " 'covar_britishcolumbia', 'covar_manitoba', 'covar_newbrunswick',\n", " 'covar_newfoundlandandlabrador', 'covar_novascotia', 'covar_ontario',\n", " 'covar_princeedwardisland', 'covar_quebec', 'covar_saskatchewan'],\n", " dtype='object')\n" ] } ], "source": [ "print(df.columns)" ] }, { "cell_type": "markdown", "id": "edd09506cf6ec942", "metadata": {}, "source": [ "### Fitting" ] }, { "cell_type": "markdown", "id": "834746ac-407f-4593-9173-e061757fdd14", "metadata": {}, "source": [ "Lifelines version" ] }, { "cell_type": "code", "execution_count": 7, "id": "a94bda3d2ac9a025", "metadata": { "ExecuteTime": { "end_time": "2026-02-27T13:50:16.592449Z", "start_time": "2026-02-27T13:50:16.417480Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "covariate\n", "covar_alberta 0.345928\n", "covar_britishcolumbia 0.173319\n", "covar_manitoba 0.292570\n", "covar_newbrunswick 0.193955\n", "covar_newfoundlandandlabrador 0.444993\n", "covar_novascotia 0.266590\n", "covar_ontario 0.319706\n", "covar_princeedwardisland 0.397880\n", "covar_quebec 0.344936\n", "covar_saskatchewan -0.049557\n", "Name: coef, dtype: float64\n" ] } ], "source": [ "from lifelines.fitters.coxph_fitter import CoxPHFitter\n", "# Lifelines model fit\n", "\n", "covar_cols = [col for col in df.columns if \"covar\" in col]\n", "ll_model = CoxPHFitter()\n", "ll_model.fit(\n", " df=df,\n", " duration_col=\"diff_days\",\n", " event_col=\"observed\",\n", " formula=\"~ \" + \" + \".join(covar_cols),\n", ")\n", "print(ll_model.params_)" ] }, { "cell_type": "code", "execution_count": 8, "id": "7ad0d134-2103-49ea-bac1-5fbb31710c08", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAYbJJREFUeJzt3Xd4FOXexvHv7iabRgoQkhBICL13ISBNJVJERNQjKoqiYsNyxIoFjuUVu6hwRJFqAxvIEUQpIh2k954QBBIIkEJC2u68fwwEIyEkIcmm3J/r2mtnZ5+Z/U0Wk9uZZ57HYhiGgYiIiEgFYXV1ASIiIiLFSeFGREREKhSFGxEREalQFG5ERESkQlG4ERERkQpF4UZEREQqFIUbERERqVDcXF1AaXM6nRw5cgRfX18sFouryxEREZECMAyDlJQUQkNDsVrzPzdT6cLNkSNHCAsLc3UZIiIiUgSHDh2idu3a+bapdOHG19cXMH84fn5+Lq5GRERECiI5OZmwsLCcv+P5qXTh5tylKD8/P4UbERGRcqYgXUrUoVhEREQqFIUbERERqVAUbkRERKRCqXR9bkRERMorh8NBVlaWq8soMXa7/ZK3eReEwo2IiEgZZxgGcXFxJCYmurqUEmW1Wqlbty52u/2y9qNwIyIiUsadCzZBQUF4e3tXyEFozw2ye/ToUcLDwy/rGBVuREREyjCHw5ETbKpXr+7qckpUjRo1OHLkCNnZ2bi7uxd5P+pQLCIiUoad62Pj7e3t4kpK3rnLUQ6H47L2o3AjIiJSDlTES1H/VFzHqHAjIiIiFYpLw83SpUvp378/oaGhWCwWZs+efcltlixZQrt27fDw8KBBgwZMnTq1xOsUERGR8sOl4SY1NZXWrVszfvz4ArWPjo6mX79+XH311WzatIl///vf3H///fz6668lXKmIiIiUFy4NN3379uX1119n4MCBBWo/YcIE6taty3vvvUfTpk159NFHueWWW/jggw9KuNKCOXMqnqN7N7q6DBERkTJl/PjxRERE4OnpSWRkJGvXri3RzytXfW5WrVpFVFRUrnW9e/dm1apVF90mIyOD5OTkXI+SsGXxN3h92IjMGXeXyP5FRETKo5kzZzJixAhGjx7Nhg0baN26Nb179+bYsWMl9pnlKtzExcURHByca11wcDDJycmcOXMmz23GjBmDv79/ziMsLKxEaqtatz0AdRwHOTiuP+kxJZtKRUSk8jIMg7TMbJc8DMMoVK3vv/8+w4YNY+jQoTRr1owJEybg7e3N5MmTS+inUwkG8Rs5ciQjRozIeZ2cnFwiAad2REPmVL+P6xKmUCdhKVlT+5J281S8W/Yv9s8SEZHK7UyWg2ajXNPfdMervfG2Fyw+ZGZmsn79ekaOHJmzzmq1EhUVle9Vl8tVrs7chISEEB8fn2tdfHw8fn5+eHl55bmNh4cHfn5+uR4lwWKx0P/R9/il+4+sM5rgTjaeP9zFmT/GlsjniYiIlHUJCQk4HI48r7rExcWV2OeWqzM3nTt3Zt68ebnWLViwgM6dO7uootwsFgv9e17NhojZ/DD9YW62/I7H7//BWa8b1rD2ri5PREQqCC93Gzte7e2yzy7rXBpuTp8+zb59+3JeR0dHs2nTJqpVq0Z4eDgjR47k8OHDTJ8+HYCHHnqIcePG8eyzz3LvvfeyePFivv32W+bOneuqQ8hTu/o18XhwGnM/vZ1+lhWc/PFJqj22BIphGncRERGLxVLgS0OuFBgYiM1my/OqS0hISIl9rkv/2q5bt462bdvStm1bAEaMGEHbtm0ZNWoUAEePHiU2Njanfd26dZk7dy4LFiygdevWvPfee3z++ef07u2a9Jqf5qH+HO34EqmGB9VObcaxR2PxiIhI5WK322nfvj2LFi3KWed0Olm0aFGJXnWxGIXt9lzOJScn4+/vT1JSUon1v8n5rPQsFr55KzexmGz3KrjdNx9CWpboZ4qISMWSnp5OdHQ0devWxdPT09XlFNrMmTO5++67+fTTT+nYsSNjx47l22+/ZdeuXRf0xcnvWAvz97vsn9Mqx/w83VnX6Cla7NpNo6zDGJ9dhaX/h9D2TleXJiIiUioGDRrE8ePHGTVqFHFxcbRp04b58+dfEGyKkzqBlLB/97+C5+wjWetsjMWZjXPes5B4yNVliYiIlJpHH32UgwcPkpGRwZo1a4iMjCzRz1O4KWFBvp48clMvhvIf1jsbYs1KxZg74tIbioiISJEo3JSCa5sFM+ex7rxsPEymYcOy9zfY9qOryxIREamQFG5KSf0aVejdozsTHf0AMOY9DU6Hi6sSERGpeBRuStGw7nX50uN2Ug0PLGkn4JvbYd4zsHoCVK6b1kREREqM7pYqRd52Nx6JasZPc6/kDrffYe/fxr4JCIcm17muOBERkQpCZ25K2e0dw3nLci8vZt3LstChnKja2nxjxu0Qt821xYmIiFQACjelzM1mpVOjUL5yRHHXgWu5K24QJ40q5pufdoc1n7m2QBERkXJO4cYFXh3Qgutb1aRHoxqc9G1Cz4x3+csaCoYDVn6k/jciIiKXQeHGBYL9PBl3Rzum3duRx3s25BR+XJv2OpkWD0g6BP97HDJOu7pMERGRcknhxsVu7xjG070acQZP3sm82Vy5YTp80hmObHRtcSIiIpdh6dKl9O/fn9DQUCwWC7Nnzy6Vz1W4cTGLxcLwqxvwZFQjJjqu5wFexPCrDYmx5m3iIiIi5VRqaiqtW7dm/Pjxpfq5CjdlgMVi4dFrGhDk68Fv6c15vdr/mW/89SeknXRtcSIiIkXUt29fXn/9dQYOHFiqn6twU0bYrBae7dMEqwUm7XLjlMXffOPdRrDlO9cWJyIiZYthQGaqax7l4KYXDeJXhtzSvjZhVb149JuNPJ86lLc9p+DvTIKfhkNISwhq4uoSRUSkLMhKgzdCXfPZLxwBu49rPruAdOamjImsV51P72rPr86OdHV8ihHeGRwZMPFq2L/Y1eWJiIiUeTpzUwa1rOWPu81CSoaTj3yf5DG/WKzJh+H7e+GBJVA1wtUlioiIK7l7m2dQXPXZZZzCTRnkbrMS1TSYX7bF8cH6bLbVH89n9uewJOyGidfA0PlQo5GryxQREVexWMr8pSFX0mWpMuqDQW145Ybm2G1WFuxP49d2n0K1+pB2AmbcAelJri5RREQkX6dPn2bTpk1s2rQJgOjoaDZt2kRsbGyJfq7CTRnl6W7j7isjuCMyHICH5hzhp3aTwDsQTuyF9VNdW6CIiMglrFu3jrZt29K2bVsARowYQdu2bRk1alSJfq7CTRn3fN8mdKpXDYCXFhwjseO/zTfWT4XsDJfVJSIicilXXXUVhmFc8Jg6dWqJfq7CTRnn6W7j6/s70TzUj5SMbG5dUw+HTxCcPAA75ri6PBERkTJH4aYcsFotfDbkCoJ8PdhzCn63djbf2PZ9uRhMSUREpDQp3JQTtQK8+OK+SNysFj5OaI/TYoM982HLt64uTUREpExRuClHGof4cmenOmw2GvBL9bvNlb88Y16iEhEREUDhptwZ0MYcbvupI1eTGdjCvCX8p8fAkeXiykREpCQZlaAbQnEdo8JNOdM2vCqRdauR7rQxv8FLgAUOLoeF/3F1aSIiUgLc3d0BSEtLc3ElJS8zMxMAm812WfvRCMXlUPs6VVkTfZLp0f5c3+M5rH+8CavGQfUG0KgP+NV0dYkiIlJMbDYbAQEBHDt2DABvb28sFouLqyp+TqeT48eP4+3tjZvb5cUThZtyaFCHMKaujGHdwVPM7ngnNzXfA9t/hJ//bTa4eRK0vMWlNYqISPEJCQkByAk4FZXVaiU8PPyyw5vCTTlUp7oPj17TgLfn7+aNX3bT45F3qV4lCNZMMBtsmK5wIyJSgVgsFmrWrElQUBBZWRW3j6XdbsdqvfweMwo35dS9Xery9ZpY/jp1ho+WHeWVAW9BuyHwyZVwaC2kJoBPoKvLFBGRYmSz2S67P0ploA7F5ZSnu42Xr28GwLRVB1m65zgENYOabSD7DCx+3bUFioiIuIjCTTl2bdNgWtf2B+CFWVtJy3JAj+fMN9dPgc0zXVidiIiIayjclGNWq4UpQzsS4O3OX6fOMGz6OibGN8JR/1qzwU/DYeU4cDpcW6iIiEgpUrgp56r52Hm5XzMsFlix7wT/98tuxtUYBc0GgDMLfnsRPo8y++GIiIhUAgo3FcDN7Wvz48NX0josAIAlB1LgX9Og/4dgrwJHNsD0AXB4g2sLFRERKQUKNxVE2/CqjB3UBoDtR5JJyciG9vfAYxugbg/ISoNJveDnEZCR4tJaRURESpLCTQUSUd2bhkFVyMx2MmPtIXOlbzDcMgVqdzAvU62bBJ/2gFMxLq1VRESkpCjcVCAWi4V7u9YF4M35u9h2OMl8w6c63LcAhvwEPkFwcj+M76QzOCIiUiEp3FQwA9vWIqK6Nw6nwePfbORYcrr5hsUC9a6CO38wX2efgV3zXFaniIhISVG4qWA83W18cV8kNXw9OJCQynUfLWP53oTzDWq2gpCW5vKu/4HT6ZpCRURESojCTQUUVs2bbx/sTJMQXxJOZ3LnpDWM+mkbhmGYDToMM593/s8cC0dERKQCUbipoOoG+jB7eBfu7BQOwPRVB/l8WbQZcNreCd2eMhtu/hq2z3JhpSIiIsVL4aYC83S38fqNLXm2T2MA/m/eTm7+ZCXbjp6GnqOg86Nmw5+fhMw0F1YqIiJSfBRuKoFh3epxX9e6eLhZ2RCbyL8mrOKvU2kQ9QoE1IEzp2D5+64uU0REpFgo3FQC7jYrL1/fjGXPXk29Gj6cyXJw/cfLSXMAvV4zG62eYIYcERGRck7hphIJ8vPko9vaApCYlsXSPcehSX+o0QQyU+DHB1xcoYiIyOVTuKlkWtTyp26gDwD7j6eC1Qq9Xjff3PsbfHET/LXehRWKiIhcHoWbSqh/q5oAfLx4L6sPnICG10KP58Fig/2LYEofOLbTxVWKiIgUjcJNJXR7ZDih/p6kZzkZ/tUG0jKz4eqR8Ng6qNEUHJnmBJsa4E9ERMohhZtKqKa/F7880R1fTzdOpGbyf3PPnqWpVg/umAnuPhC7Erb94NpCRUREikDhppLy93bn0zvbAzDjz0PEn5uDqmod6HCfubxhmouqExERKTqFm0rsygaBXFGnKg6nwbjF+86/0W4IWKwQswxiVriuQBERkSJQuKnkRlzbCIAfNvxFSnqWuTKwoRlwAFaNd1FlIiIiRaNwU8l1rl+d+jV8SMt08NrPO86/0ekR83n3XFg32TXFiYiIFIHLw8348eOJiIjA09OTyMhI1q5dm2/7sWPH0rhxY7y8vAgLC+PJJ58kPT29lKqteCwWC4/3bAjAt+v+YtLyaDKznVCjMdTvaTZa85kLKxQRESkcl4abmTNnMmLECEaPHs2GDRto3bo1vXv35tixY3m2//rrr3n++ecZPXo0O3fuZNKkScycOZMXXnihlCuvWG5oHcpVjWsA8NrPO+j/8XLzEtV175gNju+ELd+6sEIREZGCc2m4ef/99xk2bBhDhw6lWbNmTJgwAW9vbyZPzvsyyMqVK+nSpQt33HEHERER9OrVi9tvv/2SZ3skfxaLhc+HXMEbA1tis1rYHZ/Cl6tjoXr9831v/vdvzRwuIiLlgsvCTWZmJuvXrycqKup8MVYrUVFRrFq1Ks9trrzyStavX58TZg4cOMC8efO47rrrLvo5GRkZJCcn53rIhdxsVu6IDOfNm1oCMG7xXo6lpEO/D8CrGmSlwjeDNLCfiIiUeS4LNwkJCTgcDoKDg3OtDw4OJi4uLs9t7rjjDl599VW6du2Ku7s79evX56qrrsr3stSYMWPw9/fPeYSFhRXrcVQ0N7erTbOafqRmOvhy1UGwuZ2feyp6Kfz1p2sLFBERuQSXdygujCVLlvDGG2/w3//+lw0bNvDjjz8yd+5cXnvttYtuM3LkSJKSknIehw4dKsWKyx+r1cKdneoA8NuOeAzDgLaDz3cuXvkROLJdWKGIiEj+XBZuAgMDsdlsxMfH51ofHx9PSEhIntu8/PLL3HXXXdx///20bNmSgQMH8sYbbzBmzBicF7lc4uHhgZ+fX66H5K9X82C87TZ2xaWweNfZzt1tB5vPu37WreEiIlKmuSzc2O122rdvz6JFi3LWOZ1OFi1aROfOnfPcJi0tDas1d8k2mw3APMMgxSKwigd3dAwHzKkZAGhx8/mxb355BuY9C06HiyoUERG5OJdelhoxYgQTJ05k2rRp7Ny5k4cffpjU1FSGDh0KwJAhQxg5cmRO+/79+/PJJ58wY8YMoqOjWbBgAS+//DL9+/fPCTlSPAZ1MPsmLdgRz6KdZ8+uXTUS6l1tLq/9FGY/AhmnXVShiIhI3txc+eGDBg3i+PHjjBo1iri4ONq0acP8+fNzOhnHxsbmOlPz0ksvYbFYeOmllzh8+DA1atSgf//+/N///Z+rDqHCahjsS9vwADbGJjLu9330bBoMnn4wZDZs+Q5+HAZbZsCB3+H+hRAQ7uqSRUREALAYlex6TnJyMv7+/iQlJan/zSVsiD3FTf9dCcCyZ68mrJr3+Tc3z4S5T0FmCnQYBv3edVGVIiJSGRTm73e5ultKSleLUH/qBvoA8NR3m3P3a2o9CG48O6nmnxNh4SuQrjGERETE9RRu5KLsblYm3X0FAGujT/JnzKncDRr1NR8Ay9+HD1tB/PZSrlJERCQ3hRvJV70aVehYtxoAD325nqQzWeffdLPDHTNg4KcQUAfOnIL5z4Mj6yJ7ExERKXkKN3JJ797SGoCTqZlsPpR4YYPWt8GdP4DVzRzFeMHo0i1QRETkbxRu5JLCq3vTr1VNAD5dup/M7DwGTAxsCNeeHSl6zy+lWJ2IiEhuCjdSILd1CMNqgRX7TvDN2ti8G7W+DazucPIAbPuhdAsUERE5S+FGCqRbwxo81KM+AP83byeTlkdf2Mi7GnR62Fz+/l7Yt7AUKxQRETEp3EiB3d+tHs1D/cjMdvLazzuYvirmwkY9R0H9a8zlPzUHlYiIlD6FGymwaj52/vdoVx67pgEAo37azpMzN5GR/bc5pmzu5/ve7FsAiZqFXURESpfCjRSK1WphxLWNuKV9bQBmbTzM9+v/yt0ouDnU6QqOTFg4WreGi4hIqVK4kUKzWCy8+6/WPN6zIQAT/thPcnrW3xtAl8fN5W0/wNeDICvdBZWKiEhlpHAjRTaoQxh2m5VDJ89ww8fL+W173Pk3G/WGAePNsW/2L4JZD0LlmsZMRERcROFGiqxWgBdThnagmo+dmBNpPPDFeuZv+1vAaXsn3D4TsMCO2bDyI1eVKiIilYjCjVyWLg0CWTiiBze0DgXgo0V7Sc3IPt+gYRT0eM5cXjAKlo+F7MzSL1RERCoNhRu5bNV87Lx8fTOqeLix42gyr8/dkbtB92cgopu5vHA0/Dqy9IsUEZFKQ+FGikUNXw9eHdAcgFX7T+R+0+YGd/4IXZ4wX//5Oaz8WJ2MRUSkRCjcSLG5qnEQVgvEnEhj/cFTud90s0PUKxAWab7+7SWYeh04si/ckYiIyGVQuJFiU83HzsC25vg39037k33HTuduYLGYZ3D6vAnu3nB4PWz73gWViohIRaZwI8VqVP9m1A30ITEti2krYy5s4FHFnH/qysfM11tm6hZxEREpVgo3Uqz8vdx5qV9TABbujMfhvEhwaXqD+bx/MSx7r5SqExGRykDhRopdlwaB+Hm6cTQpPffAfn8X0gL6vGUuL34NYleXXoEiIlKhKdxIsfN0tzG4Ux0A3py/68K+N+d0egjq9jCXf32xlKoTEZGKTuFGSsQdHcOp5mPn4Ik0hkxaw+HEM3k3HDAOsMDhdXAyulRrFBGRiknhRkpEWDVvfn6sK7UCvDiSlM5jX2/Iu/9NQDjU6WIur5+qzsUiInLZFG6kxIQGePHV/ZHY3axsiE3kw0V7827Y+jbzecVYmP1IqdUnIiIVk8KNlKiIQB9G9m0CmPNOfbM2Fuc/z+C0vROuuNdc3va9Ri4WEZHLonAjJW5I5wgi61YDYOSPW/l48b7cDSwW6Ps2eAeCIxM+j4K4bS6oVEREKgKFGylxNquFqUM7cluHMADGLtrDnviUfzRyh1smg7sPxG+F6QNg3WSdxRERkUJTuJFS4WW38eqAFgT7eWAY0HvsUj7+Zx+cej3g0bVQtS6kJcDPT8J7jWD+C5B6Iu8di4iI/IPCjZQau5uVL+6LpFvDQAwD3luwh22Hk3I38q8ND68w55/yD4P0JFg9HiZ0hdQE1xQuIiLlisKNlKpGwb58cV8kvZoFA/C/LUcubGT3MeefemIz3PEdBNSBlCOwYHQpVysiIuWRwo24RN+WIQD8tPEIJ1Mz825ktUGjXnDzJPP1pi8hMbaUKhQRkfJK4UZcIqppMLUCvIhLTueJGRsx8hu8L6wD1OlqLm/9vnQKFBGRckvhRlzC19Odyfd0wGqBZXsT+GHD4fw3aD3IfF72PiTsy7+tiIhUago34jKNQ3y5IzIcgAl/7M+/cZvBULMNZKbAmgklX5yIiJRbCjfiUs/2aYKb1cK+Y6fZf/wis4eD2f+m86Pm8s7/6dZwERG5KIUbcSk/T3c6RJijF/+6PS7/xg2vBd9QOB0H066H5DzutBIRkUpP4UZc7tydU2/P382HCy8yuSaAVwAMmQ2e/nBsByx5s1TqExGR8kXhRlzuzsg6XN24BgCTV0Tn37hG4/O3hm+fBZlpJVydiIiUNwo34nJWq4Xxg9thtUDSmSzGLtxDZrbz4hvU7wkB4ZCRDDvnlF6hIiJSLijcSJngbXfj0WsaAjB24V6enLmJ9CxH3o2tVmh9h7n8y7Nw4I9SqlJERMoDhRspM56Masjo/s0AmLv1KFe/u4QpK6JJTs+6sHHHB6BKiDn31PQbYM5jpVytiIiUVQo3UmZYLBaGdqnLh7e1oaa/J0eT0nnlfzvo/vbv/LTpH4P8+VSHB5ZARDfz9YbpsOQtcOZzOUtERCoFi5HvuPcVT3JyMv7+/iQlJeHn5+fqcuQi0rMcfLvuENNWxrD/eCoAPz5yJe3Cq+Zu6HTCp90hfqv5OrQdtPwXNLvBnGFcREQqhML8/daZGymTPN1tDOkcwa//7k6/VjUBeGPuzgvnoLJa4aFl0PdtsFjhyAb4dSSM66CB/kREKimFGynT3GxWXu7XDLublXUHT7E7PuXCRhYLRD4Ij22A3m+ATw3ISoPFr5V+wSIi4nIKN1Lmhfh70q1BIADztuYzinG1utB5ODS9wXy9fgrMexayzpRClSIiUlYo3Ei5cEObUAAmLTvAwROp+TeOGg2tbzeX134Kk3vrEpWISCWicCPlQv9WobSu7U9qpoOpK2Pyb+zpDwMnwJ0/gFdVOLoZlr9fKnWKiIjrKdxIuWC1WrizUx0Alu9NINtRgFu+G0TB1S+ay3FbSrA6EREpSxRupNzo2jAQL3cbe4+d5vkft5JVkIATFmk+Ry+FRepgLCJSGSjcSLlR09+Lsbe1AeD79X8x5VKTbALUbAWRD5nLy96FtJMlV6CIiJQJCjdSrvRuHsKTUY0A+GH94Uu0PrfRGHDzNJffawI/j4Cs9BKqUEREXE3hRsqde7pE4G6zsDs+hd1xeYx7809WKwz+Dmo0AUcGrJsEPz0CzotMzCkiIuWawo2UO/5e7vRoFATAN2tjC7ZR3e7wyGro/xFggW0/wLdDYN9CSE8uuWJFRKTUKdxIuXRL+1oATF0Zw7PfbyY9qwBnYSwWaH833PgJWGyw62f48mZ4ux6s+ayEKxYRkdKicCPlUq9mIQw4O7Dft+v+4rqPlhV84za3w9BfoNVt4B8Gziz47UU4vqeEqhURkdLk8nAzfvx4IiIi8PT0JDIykrVr1+bbPjExkeHDh1OzZk08PDxo1KgR8+bNK6VqpaywWi18eFtbXriuCQAHjqcy/KsNZGQXsB9NeCTc9Ck8uQ3qXQWOTPh+qKZqEBGpAFwabmbOnMmIESMYPXo0GzZsoHXr1vTu3Ztjx47l2T4zM5Nrr72WmJgYvv/+e3bv3s3EiROpVatWKVcuZcUD3evz76iGWC0wd+tRrnn3D95fsIcDx08XfCcDPwXvQIjfBt8NhX/OPC4iIuWKxTBc95s8MjKSDh06MG7cOACcTidhYWE89thjPP/88xe0nzBhAu+88w67du3C3d29QJ+RkZFBRkZGzuvk5GTCwsJISkrCz8+veA5EXG7xrniGf7WRM2f73thtVlY8fw01fD0KtoOY5TB9ADiz4drXoMvjJVitiIgUVnJyMv7+/gX6++2yMzeZmZmsX7+eqKio88VYrURFRbFq1ao8t5kzZw6dO3dm+PDhBAcH06JFC9544w0cjotfihgzZgz+/v45j7CwsGI/FnG9a5oEs/qFnnwwqDVBvh5kOpxMWl6AQf7OiegK3Z81l5e+CxmFOPMjIiJlisvCTUJCAg6Hg+Dg4Fzrg4ODiYuLy3ObAwcO8P333+NwOJg3bx4vv/wy7733Hq+//vpFP2fkyJEkJSXlPA4dOlSsxyFlh7+XOwPb1ubqxuZt4hP+2M/0VTEU+ORk92cgIBwykmDOo+DIKsFqRUSkpLi8Q3FhOJ1OgoKC+Oyzz2jfvj2DBg3ixRdfZMKECRfdxsPDAz8/v1wPqdie7t2YLg2qAzDqp+08/d2Wgk20abXCjRPAYoXts+DnJ0u4UhERKQkuCzeBgYHYbDbi4+NzrY+PjyckJCTPbWrWrEmjRo2w2Ww565o2bUpcXByZmZklWq+UHzV8Pfji3khe6tcUm9XCDxv+YnJB5qECiOgC/T80lzd+Adn6dyUiUt64LNzY7Xbat2/PokWLctY5nU4WLVpE586d89ymS5cu7Nu3D6fz/P+F79mzh5o1a2K320u8Zik/rFYL93erx/Cr6gOw/uCpgm/c9i6w+5rLO+eUQHUiIlKSXHpZasSIEUycOJFp06axc+dOHn74YVJTUxk6dCgAQ4YMYeTIkTntH374YU6ePMkTTzzBnj17mDt3Lm+88QbDhw931SFIGde2TlUA1kafLPjt4RaLOZs4wLynYf20EqpORERKgpsrP3zQoEEcP36cUaNGERcXR5s2bZg/f35OJ+PY2Fis1vP5KywsjF9//ZUnn3ySVq1aUatWLZ544gmee+45Vx2ClHGd61WndlUv/jp1hvunrWPRUz2wWCyX3vDmz2FafzixD/73OPjXhgY9S75gERG5bC4d58YVCnOfvFQMMQmpXPXuEgBG9m3CHZHh+HoWYJwkRxZ8d485B1X7e873xRERkVJXLsa5ESktEYE+dKpXDYAxv+ziqneWFOwSlc0d2gw2lzfPgLitJViliIgUF4UbqRS+vC+S//RvRoifJydSM7lv2jpOphbgTqjGfaHBtZCdDr+9XPKFiojIZVO4kUrBzWblni51mfFAJ4L9PIhOSOVfE1Yy/vd9bDqUePGB/iwWuO4dc/nAEji+u9RqFhGRolG4kUolItCHKfd0pIavB/uPp/LOr7u5cfwKXv15x8U3qlYXGvUFDPjlOXAWcOZxERFxCYUbqXSahfrxyxPd+E//ZkQ1Ne/Mm7oyhm2Hky6+0VVnJ3I98Lt595SIiJRZCjdSKQVW8eCeLnX5/O4r6NeqJoYBT87cRHrWRc7KhLaBfu+Zyxu/gv2LS61WEREpHIUbqfReG9CCwCoe7D12mqe+23zxhh3uh3ZDAAPmPA6ZaaVWo4iIFJzCjVR61XzsfHRbG6wWmLvlKD9u+OvijXu/YU7NkHQItn5XekWKiEiBKdyIAFc2CGRI5wgA3p6/G4fzIndPefhCh/vM5Z1zoHKNgSkiUi4o3Iic9XzfJvh6uBGXnM6k5Qcu3rDFTebzvoWw5tPSKU5ERApM4UbkLE93G8/1bQKYIxl/u+5Q3g1rtoYrHzOXf30Bts/SGRwRkTJE4UbkbwZHhnND61AMA975dffFB/fr/gwENQPDYc4/NetBBRwRkTJC4UbkbywWC+/8qxVWCxxPyWD7keS8G3r6wwNLoMfZGem3zIRjO0utThERuTiFG5F/8HCz4edlzho+dWXMxce+cfOAq1+AOl3M11u/BaezlKoUEZGLUbgRycOo65sB8P36v7hx/ApOZ2RfvHHLW8zn5R/AR23g0J8lX6CIiFyUwo1IHga2rcXbN7fCz9ONXXEpzFgbe/HG7e6Gbk+Dhz8kHoRJUXBsV+kVKyIiuSjciOTBYrFwa4cw7u1aF4BlexMu3thqg54vw5NboarZnrkj1MFYRMRFFG5E8nF9q5pYLfDHnuPM2Xwk/8ae/nD7N+bywRUw/3kFHBERF1C4EclHgyBfHrmqAQD//X3fpTcIagrXvWsur5kAE6+BpHymcxARkWKncCNyCfd3q4vdZmVXXApro09eeoMO98NVL5jLRzbAzDvhTGKJ1igiIucp3IhcQoC3nZva1QLgzV92kpF9kVvDz7FY4Krn4N5fwc0LjmyEX18shUpFRAQUbkQK5P5udbFYYENsIs98t6VgG4V3glsmm8tbZkDy0ZIrUEREchQ43MyZM4esrKySrEWkzGoQ5MvYQW0AmLP5CCv35XP31N817gtBzcGZDX+8qQ7GIiKloMDhZuDAgSQmJgJgs9k4duxYSdUkUiYNaFOLW6+oDcALs7aS5SjAaMQWC1zzkrm8fipsmFZyBYqICFCIcFOjRg1Wr14NgGEYWCyWEitKpKx6sV8zvNxtxJxI43+XujX8nCbXnZ9FfMt3JVeciIgAhQg3Dz30EAMGDMBms2GxWAgJCcFms+X5EKmo/L3c+dfZszfL8xvY75863G8+x66EM6dKoDIRETnHraAN//Of/3Dbbbexb98+brjhBqZMmUJAQEAJliZSNvVqFsL0VQeZtekw17euyTVNgi+9UdUICAiHxFhY8SFE/aekyxQRqbQKdbdUkyZNuP7663n66af517/+xYABAxgwYABt2rQhOjoaT09PBgwYUFK1ipQJXRpUp0/zEAwDXp69nWMp6QXbMKSV+bz8Azi8vuQKFBGp5Ip0K/jGjRuZPn06AImJiURGRvLee+9x44038sknnxRrgSJljcVi4dk+jfFws3I48Qy9PljK/G0FuM27/4fgE2QuH1xVskWKiFRiRQ433bp1A+D7778nODiYgwcPMn36dD766KNiLVCkLKpXowpzHu1K81A/EtOyePybTZw4nZH/Rj6Bf+tYPAOcBbjbSkRECq1I4SYtLQ1fX18AfvvtN2666SasViudOnXi4MGDxVqgSFnVOMSX2cO7EOLnSabDyYbYxEtv1HwgWN0hbits+77EaxQRqYyKFG4aNGjA7NmzOXToEL/++iu9evUC4NixY/j5+RVrgSJlmbvNSlQz81LTK//bzrd/HsLhzGegvoAw6PK4ubxgNGScLoUqRUQqlyKFm1GjRvH0008TERFBZGQknTt3BsyzOG3bti3WAkXKumHd6lHNx85fp87w7A9bGP7VBoz8RiLu/gwE1IGUI/DHW6VXqIhIJWEx8v0tfHFxcXEcPXqU1q1bY7WaGWnt2rX4+fnRpEmTYi2yOCUnJ+Pv709SUpLOMkmxOZWayddrY3nvt904DXhtQHPu6hxx8Q02z4BZD5rL3Z6CnqNKpU4RkfKqMH+/ixxuyiuFGylJz32/hZnrDhHg7c5vT3YnyNcz74aGAQtehpUfm68b9YUB48xOxyIicoHC/P3WrOAixeiVAc0Jr+ZNYloWt322mkMn0/JuaLHAta9B50fBYoM9v8DsR0q3WBGRCkrhRqQYebrbmDK0A17uNg4cT+Xqd5fwx57jeTe2WKD3/8F9CwAL7P0VTuwv1XpFRCoihRuRYla/RhVmPtiJiOreZDsNnpy5iT3xKRffoHZ7iOhqLv/ynHnJSkREikzhRqQEtKodwKxHulAv0IeTqZm8+r8d+W/Q732w2WHfAnPuKRERKTKFG5ESUtXHzoS72gOwfF8C+4/nM6ZNjUZw7avm8u9vQHIBpnMQEZE8KdyIlKBGwb5E1q0GwJh5u8hy5DPlQuRDENYJHBnwy7OlVKGISMWjcCNSwp7r2wS7zcrCnfFMWxlz8YYWiznmDcCe+ZASXyr1iYhUNAo3IiWsXXhVnohqCMAv2+Lyb9zwWqgSAo5MWD+lFKoTEal4FG5ESsHAtrUA2Bh7ilOpmRdvaLGc73uz6SvdOSUiUgQKNyKlIDTAi7BqXjgNmLYqJv/GTfqBuzckxsLKj0qlPhGRikThRqSUNA42hwsfu3Av2w4nXbyhRxWof425vGAU7F9cCtWJiFQcCjcipeSNm1pQzccOwJ8xJ/NvPPBTCGxsLi9+vYQrExGpWBRuREpJkK8n93aJAGD+pToWe1SBIbPN5cPr4fSxEq1NRKQiUbgRKUV9WtTEYoE10SdZerE5p87xC4Warc3l7bNKvjgRkQpC4UakFDUIqsK/2tcGYMwvu3A4L3E3VKvbzOdfX4S5T8Ff60u4QhGR8k/hRqSUPd+3KX6ebuw8msz36w/l37jjMGg+EJxZ8Ofn8HlP2DyjdAoVESmnFG5ESlk1HzsP9qgPwKiftvP7rnz609jc4ZYpMPgHqN8TMGDWg7DsPXBklU7BIiLljMKNiAsM7RJBZN1qZGQ7GTZ9HV+uPojzYpeoLBZoGAWDv4dOj5jrFr0KE7rBwZWlV7SISDmhcCPiAt52Nybf04HmoX5kOw1emr2NDxftzX8jqxV6vwH93gfv6nB8J0zpq87GIiL/oHAj4iI+Hm7MHt6FJ6MaAfDJkv2czsjOfyOLBTrcB4+ug6b9zXXLPyjhSkVEyheFGxEXcrdZeSKqIbWrepHpcDJ1RXTBNvSuBv0/ApsHHN1sjoUjIiJAGQk348ePJyIiAk9PTyIjI1m7dm2BtpsxYwYWi4Ubb7yxZAsUKWHnzt58vHgfO44kF2wj72rmnVQA80dCdkYJVSciUr64PNzMnDmTESNGMHr0aDZs2EDr1q3p3bs3x47lPyJrTEwMTz/9NN26dSulSkVKzk3tatGtYSAZ2U4e+2bDpce/OefqF8DuC4fWwISusG9RyRYqIlIOuDzcvP/++wwbNoyhQ4fSrFkzJkyYgLe3N5MnT77oNg6Hg8GDB/PKK69Qr169fPefkZFBcnJyrodIWWOxWPj49rb4ebqx/3gq36yNLdiGVevAv6aaHYwT9sCXN+kOKhGp9FwabjIzM1m/fj1RUVE566xWK1FRUaxateqi27366qsEBQVx3333XfIzxowZg7+/f84jLCysWGoXKW4B3nbu62qG9f/+vg/DKODZm4ZR8NgGqHe1+XrVeHA6SqhKEZGyz6XhJiEhAYfDQXBwcK71wcHBxMXlPbHg8uXLmTRpEhMnTizQZ4wcOZKkpKScx6FDlxgRVsSFBncKB+BIUjoLdsQXfEOvAIh80Fze9TP8NLz4ixMRKSdcflmqMFJSUrjrrruYOHEigYGBBdrGw8MDPz+/XA+RsiqwigcPdDfP3vyw4a/Cbdy4L1x/9rbwzd/ArnnFXJ2ISPng0nATGBiIzWYjPj73/6HGx8cTEhJyQfv9+/cTExND//79cXNzw83NjenTpzNnzhzc3NzYv39/aZUuUmIGtq0FwK/b49kbn1K4ja+4F9oNMZe/vxc2zwSns5grFBEp21wabux2O+3bt2fRovN3eDidThYtWkTnzp0vaN+kSRO2bt3Kpk2bch433HADV199NZs2bVJ/GqkQmtb0o34NHwBG/riVjOxC9p/p9wE0uBayz8CsB+C/neDwhhKoVESkbHL5ZakRI0YwceJEpk2bxs6dO3n44YdJTU1l6NChAAwZMoSRI0cC4OnpSYsWLXI9AgIC8PX1pUWLFtjtdlceikixeapXYwDWHTzFm7/sKtzGNje4dTp0fRI8/SFhN0zuDUvfgeSjJVCtiEjZ4vJwM2jQIN59911GjRpFmzZt2LRpE/Pnz8/pZBwbG8vRo/qFLJXLdS1r8mwfM+D8b/ORi0+qeTF2b4j6DzyxGep0AUcmLH4dxraEv9YVf8EiImWIxSjw/aYVQ3JyMv7+/iQlJalzsZRpmdlO2r22gNMZ2cx5tAutagcUbUeObNj+I/w4zHzd7m64fqw5EaeISDlRmL/f+u0mUkbZ3ax0a2jeFfjRor2kZxVx7BqbG7S6FW6eZL7eMA0+6wEn1AFfRComhRuRMmxY93rYrBYW7jzGFa8v5IVZW4secprfBFe/CPYqELcFPm4Hn18L8duLt2gRERdTuBEpw9qFV2XKPR2oFeDF6Yxsvl4Ty0NfrictM7vwO7NaocezcP9CCL/SXPfXWvjsKjhUsMlqRUTKA/W5ESkHnE6DeduO8vR3m0nPctI81I9XB7SgfZ2qRd9pwl6YfiMk/wU+QTB0HgQ2LLaaRUSKk/rciFQwVquF61uF8tX9nfD1dGP7kWRumbCSF2ZtLdpZHDCDzAO/Q/WGkHrMDDpHtxRr3SIirqBwI1KOtK9TlUUjenBzu9oYBny9JpaPF+8r+g6rBJlj4rh5mmdwJl4Dh/4svoJFRFxA4UaknAny8+S9W1vz9i2tAJixNrbwoxj/XXAzuH8RhHUCZxb88WYxVSoi4hoKNyLl1E1taxHi58mptCz+t/kyB7oMaQG93zCX9y2Eg6suv0ARERdRuBEpp9xsVu7qXAeA/5u7g82HEi9vh7XbQ/2e5vIP90Piocvbn4iIiyjciJRj93erS4tafpxKy+LG/65g5I9byHZcxizgfd8C/3Cz/82cx4qvUBGRUqRwI1KOebjZmHJPR25sE4phwDdrD/HjhsNF32FgQ7j7J7DY4MDvELOi+IoVESklCjci5VwNXw/G3taWYd3qAvDfJfuIPZFW9B1WqwfNbjCXV40vhgpFREqXwo1IBXFbx3D8vdyJOZHGfdP+xFHYmcT/rttT5vOe+ZAYWzwFioiUEoUbkQqifo0qzB7eBU93K3uPnWb0nG0UeQDykJYQ0Q0MB/w8ApyX0Y9HRKSUKdyIVCB1A31455bWAHy5OpZh09ez7XBS0XbW4X7zed8CGN/BvEVcRKQcULgRqWD6tw7llRuaA7BwZzz9xy3n/d92F/4sTvMboe/bYHWHE/vg60FwYn/xFywiUsw0caZIBbX1ryQmLN3P3C3mAH9NQny5vWM4gyPDcbMV4v9rTh+DGXfAX39CvavhrllgsZRQ1SIiedPEmSJCy9r+jL+jHaP7N8PHbmNXXAqj52xn7tZCjmZcJQgGjDfnnzrwOyx9FyrX/xOJSDmjcCNSwQ3tUpeVz/ekSYgvAJ8tPVD4fjg1GkOv183l31+Hb4dAdmYxVyoiUjwUbkQqAX9vd96/tQ12m5XtR5K5cfwK3vttN/uOnS74TjrcbwYcqzvsnAMLRpVcwSIil0HhRqSSaBbqx7wnutK1QSDZToOPF+8j6v0/eGv+LpwFGRPHYoErH4PrPzBfr/kE5jyuS1QiUuYo3IhUIg2CfJkytANv3tSSduEBAHyyZD8fLNxT8J20vg0a9TWXN0yDL26EhH3FXquISFEp3IhUMu42K7d1DOeHh6/kzk7hAHy8eB9frTlYsB3Y3OGOGXDDOLDZ4cASmBQFZ06VXNEiIoWgcCNSSVksFl4b0IIhnesA8OKsbbw8e1vBZxVvdxcMXwu+oWawWfx/kJVeghWLiBSMwo1IJWaxWHi+bxO6NggE4IvVB5myIqbgO6hWF66411z+cyJ8ciWknSz+QkVECkHhRqSS87a78eX9kTzdqxEAny07QGpGdsF30PXf5kjG3tXh5H749UVwOkqmWBGRAlC4EREA7uxUB18PN46nZPDRor0F39DmDpEPwsDPAAts/hrebwZL3iyxWkVE8qNwIyIABHjbeeTqBgB8uvQAX64uYAfjcxpGwS2TwO4Lp+NgyRhIPlIClYqI5E/hRkRyPNi9HsOvrg/AO7/uJjk9q3A7aHEzPLMX3LzM11P6wtqJkFGIwQJFRC6Two2I5LBaLYy4tjFh1bxIOpPFpGXRhd+Juxfc/rV5m/ipGJj3NPy3EyQdLvZ6RUTyonAjIrnYrBYe6mGevfl23SEyswt4a/jf1b/GvE086j+ABZIOwfhISChEXx4RkSJSuBGRC9zcrjZBvh4cTUpn5p+xRdtJtbrQ9Um4dz5UCYHMFJh6PaSeKN5iRUT+QeFGRC7g6W5jYNtaALz803b+u2Rf4W4P/7vwTnDXLHP5dBxMvwEObyimSkVELqRwIyJ5eqxnQ1rV9gfg7fm76T9uedEDTnAzuOM7czl+m2YUF5ESpXAjInmq4uHGtw925qV+TfHzdOPA8VQ+/WN/0XfYqBfcMtlcjt+mqRpEpMQo3IjIRXm627i/Wz1eu7EFAON+38fa6MuYXqHxdWb/mzOnYO5TYBjFVKmIyHkKNyJySQPa1GJg21o4DRgyeQ1fFHaAv3PcvaD/WLBYYdOXcHBlsdYpIgIKNyJSQP83sAXdGgaSnuXk5dnb+HlLEUcfbtwXWtxiLu/6ufgKFBE5S+FGRArE2+7GtKEdGdolAoBHv97I9FUxRdtZi5vM59WfwM8jIG4bOIswno6ISB4shlG5LnonJyfj7+9PUlISfn5+ri5HpNxJSc+iz9hlHE48A8DgyHCe7tWYqj72gu/E6YAf7oftP55f51UNIrpCRDeo2x1qNAaLpZirF5HyqjB/vxVuRKTQsh1O/m/eTqasiAGgVoAXC0Z0x9vuVrgdxayAFR9CzHLISs39XvOb4KbPzFnHRaTSU7jJh8KNSPFZtf8ED3+1nsS0LB7qUZ+nezXCzVaEq92OLHNgv5ilEL0MopcCBjS9wbx9XAFHpNIrzN9v9bkRkSLrXL86d0bWAWDCH/vp++EyEtMyC78jmzuER0L3Z+DuOXDHt+bEmzvnwNJ3irlqEanoFG5E5LI8ek0Dnu7VCD9PN/YeO81/5mwnI9txeTtt1Au6PW0u//EWHNl4+YWKSKWhcCMil8XT3caj1zRk3B3tsFpg9qYjXPHaQv4zZ/vl7bj1oPPLqz+5vH2JSKWicCMixaJ7oxpMuLM91X3spGRkM3VlDNEJqZfe8GKqRsCgr8zlLTNhwWjILsIlLxGpdBRuRKTY9GoewtoXowj28wCgz9ilvDR7KzFFDTlNr4cez5vLK8bCd/co4IjIJSnciEixslktTL6nA9V87GRkO/lydSw93/+DRTvji7bDq0fCLVPA6ga758Lk3pB0uHiLFpEKReFGRIpd81B/Vj5/DZPvuYIavh44nAb3TVvH7I1FDCUtbjIvUXkGwJENMOtBjWgsIhelcCMiJcLT3cY1TYKZ+1hXIutWA+Dp7zazav+Jou2wcR+4Zy5YbBCzDGYOhoyUYqxYRCoKhRsRKVFBfp5MuqcDdap7k+00uH3iaib8sb9oOwtpATdPBJsH7J4Hv71cvMWKSIWgcCMiJa6KhxszH+hMrQAvgKJfngJocTPcMslcXj8Fpg+AjV9C5mXcmSUiFYrCjYiUihB/T74eFgnArrgU1kafLPrOmlwPrc6Og3NgCfw0HL64CSrXbDIichEKNyJSaupU96F/61AA7pv6J8dS0ou2I4vFnFTz8U1w1QvmukOrYcZgOBVTLLWKSPmlcCMiperNm1rSJMSXlIxs/rf56OXtrFpduOo5uO7d87eKf3wFTLsBlrwFcduKp2gRKVcUbkSkVPl4uHFL+9oAfLJkP3O3HMXhvMzLSR2HwYNLod5V4MyC6D9gyRswoSv87wndVSVSyZSJcDN+/HgiIiLw9PQkMjKStWvXXrTtxIkT6datG1WrVqVq1apERUXl215Eyp5bO4RR09+ThNMZDP96A61f+Y0hk9cybvFe1hw4QXpWESbeDG4Od82GR1ZDv/egcT/AgPVT4ZMrYff8Yj4KESmrLIbh2h54M2fOZMiQIUyYMIHIyEjGjh3Ld999x+7duwkKCrqg/eDBg+nSpQtXXnklnp6evPXWW8yaNYvt27dTq1atS35ecnIy/v7+JCUl4efnVxKHJCIFkJSWxaQV0UxfFUNiWlau9+w2Kw9dVZ8noxpisViK/iExK8wB/5IOma8fWAKhbYu+PxFxmcL8/XZ5uImMjKRDhw6MGzcOAKfTSVhYGI899hjPP//8Jbd3OBxUrVqVcePGMWTIkEu2V7gRKVscToNdccn8GX2SP2NOsTbmJMdTMgB4vGdDRlzb6PI+IOM0TL8BDq+HxteZIx1by8RJaxEphML8/Xbpf+GZmZmsX7+eqKionHVWq5WoqChWrVpVoH2kpaWRlZVFtWrV8nw/IyOD5OTkXA8RKTtsVgvNQ/25p0tdxg9ux9oXejLq+mYAfLRoL58tLeKAf+d4VIGoV852OJ4Hk66FmOXFULmIlFUuDTcJCQk4HA6Cg4NzrQ8ODiYuLq5A+3juuecIDQ3NFZD+bsyYMfj7++c8wsLCLrtuESk5FouFe7vW5d4udQF4Y94u/oy5jDFxAOp2g5s/B3dvOLwOpvaDL2+GfYs0R5VIBVSuz82++eabzJgxg1mzZuHp6Zlnm5EjR5KUlJTzOHToUClXKSJFMax73ZzlDxfuvfw7qpoPNMfFueI+c36qfQvhy5vgv5Gw6RtwFqETs4iUSS4NN4GBgdhsNuLj43Otj4+PJyQkJN9t3333Xd58801+++03WrVqddF2Hh4e+Pn55XqISNlX09+LWY9cicUCy/clMOjTVcSeSLu8nfoGw/Xvw6N/QuRD4OEPCXtg9kPmNA4KOCIVgkvDjd1up3379ixatChnndPpZNGiRXTu3Pmi27399tu89tprzJ8/nyuuuKI0ShURF2gbXpVPBrfDx25j3cFT3D5xNZsPJV7+jqvXh75vwZPboOdoc13MMji84fL3LSIu5/LLUiNGjGDixIlMmzaNnTt38vDDD5OamsrQoUMBGDJkCCNHjsxp/9Zbb/Hyyy8zefJkIiIiiIuLIy4ujtOnT7vqEESkBPVpUZP5/+5O7apeHE48w02frLy8iTf/ztMPuo2AOl3N17MfgiMbi2ffIuIyLg83gwYN4t1332XUqFG0adOGTZs2MX/+/JxOxrGxsRw9en6I9k8++YTMzExuueUWatasmfN49913XXUIIlLCwqp58/NjXYlqGozDafDS7G2s2JdQfB/Q5w3wCYIT+2BiT1j2XvHtW0RKncvHuSltGudGpPzKyHZw22er2RibiMUCYwe1YUCbSw/eWSCpCTDvadg+y3x9x7fQqHfx7FtELlu5GedGRKQwPNxsfH1/J65vVRPDgJdmb+PTP/Zz6ORldjQG8AmEf02FNnear2cMhhUf6lZxkXJIZ25EpNxJzcim53t/EJecnrPuzk7hvH5jy8vfeWYq/DT8/Bmchr3MWcer1rn8fYtIkenMjYhUaD4ebsz/dzdev7EFneqZo5N/uTqW577fwtGkM5e3c7sP3DIFrv/AHA9n72/weU9I2FsMlYtIadCZGxEp996ev4v/LjGnabDbrLQJD6BdeFWublyDyHrVi77jYzvhh2EQvxWqBEO/96FBT3D3KqbKRaSgytXEmaVN4UakYloXc5K35+9m7T+maujRqAav3NCciECfou34ZDR8HgVpZ+/OslihekNocwd0/fflFS0iBaZwkw+FG5GKyzAM9h9PZcPBU6yOPsGcTUfIdhrYbVaGdo3gud5NsFothd/xmVOw4iPYMP18yAFoeyd0+TcENiy2YxCRvCnc5EPhRqTyiE5IZdRP21i21wwkH9/elv6tQ4u+Q8OA0/Ew+xHYf3ZkdXdveDYa3POe305Eioc6FIuIAHUDfZh+b0dublcbgNFztrNgRzxF/n86iwV8Q+CuH2HofHNdVhqM7whbvy+mqkXkcinciEiFZrFYeKlfU1rU8uNkaibDpq+j74fL+GnTYTKyL2OizDqdzbuqPPwg8SD8cB/89hJkpV96WxEpUbosJSKVQkp6FuN+38eXqw6SmmmGGn8vdwa2rcWdnerQIKhK0XZ8JhF+uB/2LTBfV6sPwxaDV0Cx1C0iJvW5yYfCjUjllpSWxbRVMXy9JjbXIIB3darDqwOaY7EUocOx0wGbvoI5j5mv+38E7e8upopFBBRu8qVwIyIADqfB8n0JfLEqhoU7jwHw0/AutA4LKPpOF74Cy983lzsMM+emqn8NWG2XX7BIJacOxSIil2CzWujRqAaf392BDhFVAXjwi/Us2hlf9J02vd68ewrgz4nw1S0wuQ/ErDDP7ohIqdCZGxGp9HYeTeb+aes4nGhO3dC/dShjB7XBVpQxcTJSIHqZ2Qdn80zISjXX+9SAxtdBuyFQ+4pirF6kctBlqXwo3IhIXtIys/lw0V4mLj2A04Bbr6jNqP7NqeLhVvSdxm2DVeNg9zxITzq/vmYbuOp5aNz3susWqSwUbvKhcCMi+ZmxNpbnf9wKQKi/J2/e3IrujWpc3k4dWRCzHLZ8C1u/A2eWOSnnzROhxc3FULVIxadwkw+FGxG5lGV7j/PCrK0cOmlepurTPIT3bm2Nz+WcxTknNQHmPgU7ZgMWePAPqNn68vcrUsGpQ7GIyGXo1rAG85/ozpDOdXC3WZi/PY47Jq5m86HEoo9ufI5PINwyGRr2Agz4eQRkZxZL3SJi0pkbEZF8LN1znHun/km20/xV6W23Ub9GFfq3rskD3esXfccn9sPEq82+ODeMg3Z3FVPFIhWTztyIiBST7o1q8L/HujKwbS3sNitpmQ62Hk7ijXm7GPz5av7Yc7xoO65eH66411ye+xT8PsacmFNELpvO3IiIFFCWw8nBE2m88+suftsRn5NFvhnWic71qxd+h2cS4dshEP2H+frGT6DNHcVWr0hFog7F+VC4EZHicOhkGn0/XMbpjGwA2oUH8NHtbald1btwOzIMmPMobPwSbB7QdjCEdYJa7SGwQQlULlI+KdzkQ+FGRIpL7Ik03l+wm9mbjgBgt1l5sEc9RlzbqHBzVGWkwA/DYM8v59dZrOYcVW3vhKLMdyVSwSjc5EPhRkSK26ZDiTz3/RZ2x6cAMOnuK+jZNLhwO3E6Yc98iFkGO/8HSYfM9VUj4MYJUKdz8RYtUs4o3ORD4UZESoJhGAybvi5nEs4hnesQWbc6PZsG4eleyIkzHdmw+FVY8aH52tMfBvwXmvTTWRyptBRu8qFwIyIlZW30SQZ/vposx/lfq+HVvLmrUx26NQqkcbBv4S5XnYyGKX0h5aj5uuW/YMB4cPMo5spFyj6Fm3wo3IhISTpxOoO10SdZE32Sn7ccIeH0+QH6gnw96NU8mAFtatE+vCrWgkzMmZECS9+FFWPN17Xaw9D54GYvmQMQKaMUbvKhcCMipeXE6QxmbzrCsr3HWX3gBOlZzpz32oYHMOnuDlTzKWBI2fk/+PFBc5bx6g3MOanCO5mTcHpXK5kDEClDFG7yoXAjIq6Qke1g9YGT/LTpMPO2HiU9y0ntql7MeqQLNXwLeJlp9y/w/X1mwDnH6ga3z4CG15ZM4SJlhEYoFhEpYzzcbPRoVIP3b23Dz491pVaAF3+dOsOTMzeR5XBeegcAjfvCiB3m3VPNbzLXObNh01clV7hIOaRwIyJSyhoE+TLujrZ4udtYvi+Bjxftxeks4El0rwBoczv8awrct8Bct30WfHcPJB3WFA4iKNyIiLhE2/CqPHqNOQLxR4v30X/cchLTCjk7eFhHiHzIXN4+Cz5oBm/Vgcl9YPMMBR2ptBRuRERc5IHu9Xi8Z0PsNivbjyRz56Q1LNoZz1+n0sgu6KWqPm9C50ehRhOw2MxZxmNXwawH4Yf7IWEvZBcyNImUc+pQLCLiYkv3HOeBL9blupsqsIoH793amu4NAws+Nk52BpzYZ95Z9cfbYDjOvmEBv1rmaMcd7jXvtBIpZ3S3VD4UbkSkLIpJSGXc7/vY+ldSzjQOADX9PXn9xhaFn87h0Fr49QWI2wbZZ86v9/SHp/aAu2cxVS5SOhRu8qFwIyJl3emMbF7/eQcz/jyUs+6mtrV4Iqohdar7FG5nhgGpx83LU1OvM9c1uR5u0x1WUr4o3ORD4UZEyosTpzN47JuNrNx/AgBfDzeev64JgyPrFG2HO36Cb4eYy60GQadHILRN8RQrUsIUbvKhcCMi5c3G2FM8/8PWnMtVvZoF061RDZrV9KVBkC/+Xu4F39miV2HZe+dfN+oDt36h6RykzFO4yYfCjYiUR9kOJ/83bydTVsTkWu9mtfDZkPZc06QQfXL+WgdrJsDW78zX3oFQ50rwDwP/2uYjLBJ8C9nPR6QEKdzkQ+FGRMqzLX8l8sfu46yNOcmyvQkANAnxZerQjoT4F7KT8PbZ8POTcObkhe95+MFj66FK0OUXLVIMFG7yoXAjIhXFoZNp9PpgKWeyzFu+O9erTkSgNxHVfejRuAZNQgrwOy47A6KXwsloSDoEibGwY/b5972qQkAdqFrHfA7vbE4DUdDb00WKicJNPhRuRKQi2RB7iiGT1nI6I/uC93o0qsE9V0ZQr4YPQb6eeNltBdvpgT9gzqNm0MlL7Q7QtD9UCQGf6uBTw7y05RMIbgWcBFSkkBRu8qFwIyIVTWpGNrvjU4hJSCXmRBo7jiSxeNcx/jldlZ+nG3dfGcFTvRoXbMfpyWbASTwIpw7Cib3mtA5ZaRffxsMPvKtDlWCzH0+jPmYYsmpAfLk8Cjf5ULgRkcrg4IlUPl8WzdK9x4lPTs81+vGEO9vTp0VI0XacfBQ2fwPx2yA1AdJOnH1OMGcoz0u1+tDxAWh3F9gLOU6PyFkKN/lQuBGRysYwDJLTsxn06Sp2xZm3kz98VX3u7FSHWgFexfUhkJ5oBp3UBDgVA/sWwt7fICPZbONXyww54Z2gZhuNkiyFonCTD4UbEamsMrIdDJm0ljXR5t1RVgtENQ3m9shwGgX7Ut3Hjqd7AfvlFPhDT8OWGbD0PUg5cn691d3st9P7DfCrWbyfKRWSwk0+FG5EpDIzDIMFO+KZujImZ+Tjv/P1cKN6FTtdGwby2oAWBZ+081JST8Cmr+DQGnPeq9Rj5np3H2g3xHwENyuez5IKSeEmHwo3IiKmbYeT+HL1QVbuP8HRpDNkOXL/Oaji4UatAC+C/DwI8vUk2M+DRsG+9G4eUvA7r/JiGHB4Pcy8K/fZnO7PQOvboXr9ou9bKiyFm3wo3IiIXOhcv5yE0xl8viyaHzb8RWa2M8+2PnYbjUN8qeptx8/LHT9Pt7PP7vh7uePn5YafpzsNgqoQ5JdPv5rMVNi/GH57yeyjc06T6+HqFyGoqcbTkRwKN/lQuBERubTTGdkcPnWG+OR0jqVkEJ+cTnxyOr/vPsahk2cKvJ/aVb3oXK86VzcJomfTIDzc8jjjk3UGts+CbT+YYcc4G6p8QyE80hw4sPVt4OlfTEcn5ZHCTT4UbkREis7pNNgVl0LMiVRS0rNIPpNN0pksktOzSD6TRXJ6NslnsjiZlklMQmqusXZq+HrQpX51Ggb70jCoCm3CAwjy/ceZneO7zck9d/8ChuP8ersvhHWE0LYQ0tIcMdk/HLyr6exOJaFwkw+FGxGR0pGSnsWG2ESW7jnOz1uOEJ+cket9m9XCjW1q0bdFCJH1quHr+bfZzTNT4fAGOLQaNs80BxDMi7sPBIRDRBdoe5c5XYRXVfDwVeipYBRu8qFwIyJS+tKzHCzbm8Ce+BT2HTvNzqPJOWPugNmP54Y2tehUrxrdGtagmo/9/MaObIjbDEc2wuGNkLDbHDn5dPzFP9BiA68A8Aw4H3i8AszRk8MioW53c1kBqNxQuMmHwo2ISNmweFc8v26LZ+HOeE6kZuast1ktRNatRvNQPxoEVaFBkC+tavvjbvvHFA5Z6ZD0Fyx9G45uMQcRPHMKstMLVoCbJ/iGmH17/GqC79mHfy0IbWeeEVL4KTMUbvKhcCMiUrZkZDtYuieBNQdOsGL/CXYeTb6gTRUPN+pU96ZWgBe1qnoR5OtJYBU7NXw9CKziQZMQX9zOhZ+sM3Am0Qw65wLPmVPmuqS/zE7LCbsvXZhXNfNsj90H7FXOPp9b/udrH/D4x3p3H7B7g7u3+drmfsmPlItTuMmHwo2ISNkWk5DK8n0J7Dt2mv3HT7P9SDIn/3ZmJy9uVgt1A33w9nDDx27D2+6Gj8fZZ7vt/HoPN0L8PKkXYKWa8yS+WQm4pcaZc2alnH2cPABxWy8+V1ZRWd3Phh0fcPc6v/z3AHTu2acGhLSAkNZQpUbx1lFOlbtwM378eN555x3i4uJo3bo1H3/8MR07drxo+++++46XX36ZmJgYGjZsyFtvvcV1111XoM9SuBERKV8cToN9x07z16k0Diee4fCpMxxPyeD46QwSTmfmeaanMKp4uOHv5Z7zCPB2J9THoKHbMXyt6XgZ6Xga6Xg60/BwnsHDOIPdcQZ3Zxru2Wdwd6Riy07LeVizUrFmnsaSnQaZaVj+ftdXUfjWNO8QC2llPoe2Ne8Wq2TKVbiZOXMmQ4YMYcKECURGRjJ27Fi+++47du/eTVBQ0AXtV65cSffu3RkzZgzXX389X3/9NW+99RYbNmygRYsWl/w8hRsRkYolPcvBrrgU0jKzSctwkJqZTVqmg9SMs89/W386PZvDiWeIPZlGSnoxn5nJk4GnJRt/tyz8bJn42bLPPmdSxXr2YcnAx5qBD5l4WzLwIoPqjmOEZ+6jWvohLOTxZzqwMTTpB/V6QI0mZy972c2zQ1brhe0rgHIVbiIjI+nQoQPjxo0DwOl0EhYWxmOPPcbzzz9/QftBgwaRmprKzz//nLOuU6dOtGnThgkTJlzy8xRuREQEINvhJCU9m8QzWSSmZZJ0JoukM1kkpmVxLCWdY8kZnMlykJntJMvhJNPhJDPbSabDMJ+zHWSdW3Y4ycp2knG2TXHxJp0mlliaW2NoZjlIC2sMTa2xuHHxs0FOixtOq3vOw/jbc86yzZ7zGpv97Gs3sLqZd5rZzi1bweqGxeoGVpu5zmozH5bz6/75vtUvmIA2A4rt5wCF+/vtVqyfXEiZmZmsX7+ekSNH5qyzWq1ERUWxatWqPLdZtWoVI0aMyLWud+/ezJ49O8/2GRkZZGScH1shOfnyTl+KiEjF4GazUtXHTlUfO+BTbPs1DINsp3E+FGU7yTgbgDKy/haSsp1kZDtywlHGuXbZTtKzHBxPySAuKZ2jySEsSGrFVynpOLPBj1Susm4myrae1pb91LEey/X5ViMbqyMbHAUfSbq47XJrWuzhpjBcGm4SEhJwOBwEBwfnWh8cHMyuXbvy3CYuLi7P9nFxcXm2HzNmDK+88krxFCwiInIJFosFd5vlwlvXL1O2w8nx0xkcTUrnSGJXDp86w5+nznDydBo4srEaWVgdWVidmeayMwubkYXFkYXVyMLt7DqrMxubYb5n+/uykY3FcGI1srHgxGo4sBoOLIYDG+ayFSc2w4GVs8s4sBnOnOVzz2medWlSrEdfOC4NN6Vh5MiRuc70JCcnExYW5sKKRERECs/NZqWmvxc1/b1oF17V1eWUaS4NN4GBgdhsNuLjc48yGR8fT0hISJ7bhISEFKq9h4cHHh4exVOwiIiIlHku7VJtt9tp3749ixYtylnndDpZtGgRnTt3znObzp0752oPsGDBgou2FxERkcrF5ZelRowYwd13380VV1xBx44dGTt2LKmpqQwdOhSAIUOGUKtWLcaMGQPAE088QY8ePXjvvffo168fM2bMYN26dXz22WeuPAwREREpI1webgYNGsTx48cZNWoUcXFxtGnThvnz5+d0Go6NjcX6t3v2r7zySr7++mteeuklXnjhBRo2bMjs2bMLNMaNiIiIVHwuH+emtGmcGxERkfKnMH+/K+YwhiIiIlJpKdyIiIhIhaJwIyIiIhWKwo2IiIhUKAo3IiIiUqEo3IiIiEiFonAjIiIiFYrCjYiIiFQoCjciIiJSobh8+oXSdm5A5uTkZBdXIiIiIgV17u92QSZWqHThJiUlBYCwsDAXVyIiIiKFlZKSgr+/f75tKt3cUk6nkyNHjuDr64vFYimWfSYnJxMWFsahQ4cq5XxVOn4df2U9/sp87KDj1/GX7vEbhkFKSgqhoaG5JtTOS6U7c2O1Wqldu3aJ7NvPz69S/gM/R8ev46+sx1+Zjx10/Dr+0jv+S52xOUcdikVERKRCUbgRERGRCkXhphh4eHgwevRoPDw8XF2KS+j4dfyV9fgr87GDjl/HX3aPv9J1KBYREZGKTWduREREpEJRuBEREZEKReFGREREKhSFGxEREalQFG6Kwfjx44mIiMDT05PIyEjWrl3r6pIKZcyYMXTo0AFfX1+CgoK48cYb2b17d642V111FRaLJdfjoYceytUmNjaWfv364e3tTVBQEM888wzZ2dm52ixZsoR27drh4eFBgwYNmDp1akkf3iX95z//ueDYmjRpkvN+eno6w4cPp3r16lSpUoWbb76Z+Pj4XPsor8cOEBERccHxWywWhg8fDlS8737p0qX079+f0NBQLBYLs2fPzvW+YRiMGjWKmjVr4uXlRVRUFHv37s3V5uTJkwwePBg/Pz8CAgK47777OH36dK42W7ZsoVu3bnh6ehIWFsbbb799QS3fffcdTZo0wdPTk5YtWzJv3rxiP95/yu/4s7KyeO6552jZsiU+Pj6EhoYyZMgQjhw5kmsfef2befPNN3O1KY/HD3DPPfdccGx9+vTJ1aa8fv+XOva8fg9YLBbeeeednDbl5rs35LLMmDHDsNvtxuTJk43t27cbw4YNMwICAoz4+HhXl1ZgvXv3NqZMmWJs27bN2LRpk3HdddcZ4eHhxunTp3Pa9OjRwxg2bJhx9OjRnEdSUlLO+9nZ2UaLFi2MqKgoY+PGjca8efOMwMBAY+TIkTltDhw4YHh7exsjRowwduzYYXz88ceGzWYz5s+fX6rH+0+jR482mjdvnuvYjh8/nvP+Qw89ZISFhRmLFi0y1q1bZ3Tq1Mm48sorc94vz8duGIZx7NixXMe+YMECAzB+//13wzAq3nc/b94848UXXzR+/PFHAzBmzZqV6/0333zT8Pf3N2bPnm1s3rzZuOGGG4y6desaZ86cyWnTp08fo3Xr1sbq1auNZcuWGQ0aNDBuv/32nPeTkpKM4OBgY/Dgwca2bduMb775xvDy8jI+/fTTnDYrVqwwbDab8fbbbxs7duwwXnrpJcPd3d3YunWry44/MTHRiIqKMmbOnGns2rXLWLVqldGxY0ejffv2ufZRp04d49VXX831b+Lvvy/K6/EbhmHcfffdRp8+fXId28mTJ3O1Ka/f/6WO/e/HfPToUWPy5MmGxWIx9u/fn9OmvHz3CjeXqWPHjsbw4cNzXjscDiM0NNQYM2aMC6u6PMeOHTMA448//shZ16NHD+OJJ5646Dbz5s0zrFarERcXl7Puk08+Mfz8/IyMjAzDMAzj2WefNZo3b55ru0GDBhm9e/cu3gMopNGjRxutW7fO873ExETD3d3d+O6773LW7dy50wCMVatWGYZRvo89L0888YRRv359w+l0GoZRsb/7f/6CdzqdRkhIiPHOO+/krEtMTDQ8PDyMb775xjAMw9ixY4cBGH/++WdOm19++cWwWCzG4cOHDcMwjP/+979G1apVc47fMAzjueeeMxo3bpzz+tZbbzX69euXq57IyEjjwQcfLNZjzE9ef+D+ae3atQZgHDx4MGddnTp1jA8++OCi25Tn47/77ruNAQMGXHSbivL9F+S7HzBggHHNNdfkWldevntdlroMmZmZrF+/nqioqJx1VquVqKgoVq1a5cLKLk9SUhIA1apVy7X+q6++IjAwkBYtWjBy5EjS0tJy3lu1ahUtW7YkODg4Z13v3r1JTk5m+/btOW3+/rM616Ys/Kz27t1LaGgo9erVY/DgwcTGxgKwfv16srKyctXdpEkTwsPDc+ou78f+d5mZmXz55Zfce++9uSaWrcjf/d9FR0cTFxeXq1Z/f38iIyNzfd8BAQFcccUVOW2ioqKwWq2sWbMmp0337t2x2+05bXr37s3u3bs5depUTpvy8DNJSkrCYrEQEBCQa/2bb75J9erVadu2Le+8806uy5Dl/fiXLFlCUFAQjRs35uGHH+bEiRM571WW7z8+Pp65c+dy3333XfBeefjuK93EmcUpISEBh8OR65c6QHBwMLt27XJRVZfH6XTy73//my5dutCiRYuc9XfccQd16tQhNDSULVu28Nxzz7F7925+/PFHAOLi4vL8OZx7L782ycnJnDlzBi8vr5I8tIuKjIxk6tSpNG7cmKNHj/LKK6/QrVs3tm3bRlxcHHa7/YJf7MHBwZc8rnPv5dfG1cf+T7NnzyYxMZF77rknZ11F/u7/6Vy9edX692MJCgrK9b6bmxvVqlXL1aZu3boX7OPce1WrVr3oz+TcPsqC9PR0nnvuOW6//fZcEyM+/vjjtGvXjmrVqrFy5UpGjhzJ0aNHef/994Hyffx9+vThpptuom7duuzfv58XXniBvn37smrVKmw2W6X5/qdNm4avry833XRTrvXl5btXuJFchg8fzrZt21i+fHmu9Q888EDOcsuWLalZsyY9e/Zk//791K9fv7TLLFZ9+/bNWW7VqhWRkZHUqVOHb7/9tsz80S0tkyZNom/fvoSGhuasq8jfvVxcVlYWt956K4Zh8Mknn+R6b8SIETnLrVq1wm638+CDDzJmzJgyORR/Ydx22205yy1btqRVq1bUr1+fJUuW0LNnTxdWVromT57M4MGD8fT0zLW+vHz3uix1GQIDA7HZbBfcORMfH09ISIiLqiq6Rx99lJ9//pnff/+d2rVr59s2MjISgH379gEQEhKS58/h3Hv5tfHz8ytTISIgIIBGjRqxb98+QkJCyMzMJDExMVebv3/HFeXYDx48yMKFC7n//vvzbVeRv/tz9eb333RISAjHjh3L9X52djYnT54sln8TZeF3x7lgc/DgQRYsWJDrrE1eIiMjyc7OJiYmBij/x/939erVIzAwMNe/94r+/S9btozdu3df8ncBlN3vXuHmMtjtdtq3b8+iRYty1jmdThYtWkTnzp1dWFnhGIbBo48+yqxZs1i8ePEFpxTzsmnTJgBq1qwJQOfOndm6dWuu/+jP/VJs1qxZTpu//6zOtSlrP6vTp0+zf/9+atasSfv27XF3d89V9+7du4mNjc2pu6Ic+5QpUwgKCqJfv375tqvI333dunUJCQnJVWtycjJr1qzJ9X0nJiayfv36nDaLFy/G6XTmBL/OnTuzdOlSsrKyctosWLCAxo0bU7Vq1Zw2ZfFnci7Y7N27l4ULF1K9evVLbrNp0yasVmvO5ZryfPz/9Ndff3HixIlc/94r8vcP5hnc9u3b07p160u2LbPffbF1Ta6kZsyYYXh4eBhTp041duzYYTzwwANGQEBArjtHyrqHH37Y8Pf3N5YsWZLr9r60tDTDMAxj3759xquvvmqsW7fOiI6ONn766SejXr16Rvfu3XP2ce524F69ehmbNm0y5s+fb9SoUSPP24GfeeYZY+fOncb48ePLxO3QTz31lLFkyRIjOjraWLFihREVFWUEBgYax44dMwzDvBU8PDzcWLx4sbFu3Tqjc+fORufOnXO2L8/Hfo7D4TDCw8ON5557Ltf6ivjdp6SkGBs3bjQ2btxoAMb7779vbNy4MeduoDfffNMICAgwfvrpJ2PLli3GgAED8rwVvG3btsaaNWuM5cuXGw0bNsx1K3BiYqIRHBxs3HXXXca2bduMGTNmGN7e3hfcDuvm5ma8++67xs6dO43Ro0eXyq3Q+R1/ZmamccMNNxi1a9c2Nm3alOv3wbm7X1auXGl88MEHxqZNm4z9+/cbX375pVGjRg1jyJAh5f74U1JSjKefftpYtWqVER0dbSxcuNBo166d0bBhQyM9PT1nH+X1+7/Uv33DMG/l9vb2Nj755JMLti9P373CTTH4+OOPjfDwcMNutxsdO3Y0Vq9e7eqSCgXI8zFlyhTDMAwjNjbW6N69u1GtWjXDw8PDaNCggfHMM8/kGuvEMAwjJibG6Nu3r+Hl5WUEBgYaTz31lJGVlZWrze+//260adPGsNvtRr169XI+w5UGDRpk1KxZ07Db7UatWrWMQYMGGfv27ct5/8yZM8YjjzxiVK1a1fD29jYGDhxoHD16NNc+yuuxn/Prr78agLF79+5c6yvid//777/n+e/97rvvNgzDvB385ZdfNoKDgw0PDw+jZ8+eF/xcTpw4Ydx+++1GlSpVDD8/P2Po0KFGSkpKrjabN282unbtanh4eBi1atUy3nzzzQtq+fbbb41GjRoZdrvdaN68uTF37twSO+5z8jv+6Ojoi/4+ODfu0fr1643IyEjD39/f8PT0NJo2bWq88cYbuf74l9fjT0tLM3r16mXUqFHDcHd3N+rUqWMMGzbsgv9ZLa/f/6X+7RuGYXz66aeGl5eXkZiYeMH25em7txiGYRTfeSARERER11KfGxEREalQFG5ERESkQlG4ERERkQpF4UZEREQqFIUbERERqVAUbkRERKRCUbgRERGRCkXhRkRERCoUhRsRKReWLFmCxWK5YBJTEZF/0gjFIlImXXXVVbRp04axY8cCkJmZycmTJwkODsZisbi2OBEp09xcXYCISEHY7XZCQkJcXYaIlAO6LCUiZc4999zDH3/8wYcffojFYsFisTB16tRcl6WmTp1KQEAAP//8M40bN8bb25tbbrmFtLQ0pk2bRkREBFWrVuXxxx/H4XDk7DsjI4Onn36aWrVq4ePjQ2RkJEuWLHHNgYpIidCZGxEpcz788EP27NlDixYtePXVVwHYvn37Be3S0tL46KOPmDFjBikpKdx0000MHDiQgIAA5s2bx4EDB7j55pvp0qULgwYNAuDRRx9lx44dzJgxg9DQUGbNmkWfPn3YunUrDRs2LNXjFJGSoXAjImWOv78/drsdb2/vnEtRu3btuqBdVlYWn3zyCfXr1wfglltu4YsvviA+Pp4qVarQrFkzrr76an7//XcGDRpEbGwsU6ZMITY2ltDQUACefvpp5s+fz5QpU3jjjTdK7yBFpMQo3IhIueXt7Z0TbACCg4OJiIigSpUqudYdO3YMgK1bt+JwOGjUqFGu/WRkZFC9evXSKVpESpzCjYiUW+7u7rleWyyWPNc5nU4ATp8+jc1mY/369dhstlzt/h6IRKR8U7gRkTLJbrfn6ghcHNq2bYvD4eDYsWN069atWPctImWH7pYSkTIpIiKCNWvWEBMTQ0JCQs7Zl8vRqFEjBg8ezJAhQ/jxxx+Jjo5m7dq1jBkzhrlz5xZD1SJSFijciEiZ9PTTT2Oz2WjWrBk1atQgNja2WPY7ZcoUhgwZwlNPPUXjxo258cYb+fPPPwkPDy+W/YuI62mEYhEREalQdOZGREREKhSFGxEREalQFG5ERESkQlG4ERERkQpF4UZEREQqFIUbERERqVAUbkRERKRCUbgRERGRCkXhRkRERCoUhRsRERGpUBRuREREpEL5f2pRJsUw41bQAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Lifelines sf\n", "max_offset = 0\n", "\n", "X = df.filter(regex=\"covar\").iloc[:2]\n", "\n", "sf_lifelines = ll_model.predict_survival_function(\n", " X=X,\n", " times=range(\n", " df[\"diff_days\"].min().astype(int),\n", " df[\"diff_days\"].max().astype(int) + max_offset,\n", " ),\n", ")\n", "\n", "sf_lifelines.columns = X.index\n", "sf_lifelines.plot(xlabel=\"time\", ylabel=\"sf\")" ] }, { "cell_type": "markdown", "id": "913c06f0-2d2b-49aa-bcc9-de420d3e01a6", "metadata": {}, "source": [ "ReLife version " ] }, { "cell_type": "code", "execution_count": 9, "id": "8f0742e4-ebf1-4d90-9f88-6d37f4700e4e", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from relife.lifetime_models import SemiParametricProportionalHazard\n", "\n", "cox_model = SemiParametricProportionalHazard()\n", "cox_model.fit(\n", " df[\"diff_days\"].to_numpy(),\n", " covar=df.filter(regex=\"covar\").to_numpy(),\n", " event=df[\"observed\"].to_numpy(dtype=bool),\n", ")" ] }, { "cell_type": "code", "execution_count": 10, "id": "2f78e90c-b1d4-4f2e-9ee3-dc74aa0740ae", "metadata": {}, "outputs": [], "source": [ "timeline, sf_values = cox_model.sf(covar=X, se=False)" ] }, { "cell_type": "code", "execution_count": 11, "id": "e98ff78d-a015-4c28-b5d3-72a8dfd1cb27", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAWZ9JREFUeJzt3Xd4FOXexvHv7iabRgoQUoDQewsQIUSkaRQRATsHOKKoqIgFsSA2XisqdkVRLOCxgCjYKIoBVCDSo9TQCSUJPQ3Sduf9YyAYSSAJSTbl/lzXXvvMzDOzv8licjvlGYthGAYiIiIiLmJ1dQEiIiJSvSmMiIiIiEspjIiIiIhLKYyIiIiISymMiIiIiEspjIiIiIhLKYyIiIiISymMiIiIiEu5ubqAonA6nRw4cABfX18sFouryxEREZEiMAyDtLQ06tati9Va+PGPShFGDhw4QFhYmKvLEBERkRLYu3cv9evXL3R5pQgjvr6+gLkzfn5+Lq5GREREiiI1NZWwsLC8v+OFqRRh5PSpGT8/P4URERGRSuZ8l1joAlYRERFxKYURERERcSmFEREREXEphRERERFxKYURERERcSmFEREREXEphRERERFxKYURERERcSmFEREREXGpYoeR33//nQEDBlC3bl0sFgvffffdeddZsmQJnTt3xsPDg2bNmjFt2rQSlCoiIiJVUbHDSEZGBuHh4UyePLlI/Xft2kX//v3p06cPcXFxjBkzhjvuuIOff/652MWKiIhI1VPsZ9P069ePfv36Fbn/lClTaNy4Ma+99hoArVu3ZunSpbzxxhv07du3uB8vIiIiVUyZXzMSGxtLdHR0vnl9+/YlNja20HWysrJITU3N9yoLa7+eyLrJw9m37e8y2b6IiIicX5mHkaSkJIKDg/PNCw4OJjU1lZMnTxa4zsSJE/H39897hYWFlUltPlvn0OnQ95xc8DQ4csrkM0REROTcKuTdNOPHjyclJSXvtXfv3jL5nGS/9gA0P7KY9BebkrF7VZl8joiIiBSuzMNISEgIycnJ+eYlJyfj5+eHl5dXget4eHjg5+eX71UWwm9/l1mBowGo4Ugh+/OhGNsWlslniYiISMHKPIxERUURExOTb97ChQuJiooq648+L38fL24Y/QLfXrqEJKMmNXMPYvniBjI3/OTq0kRERKqNYoeR9PR04uLiiIuLA8xbd+Pi4khISADMUyzDhw/P63/33Xezc+dOHn30UbZs2cJ7773H119/zYMPPlg6e3CBLBYL1/fsxMJes/neeQkAWT8+AjmZLq5MRESkeih2GFm9ejWdOnWiU6dOAIwdO5ZOnTrx9NNPA5CYmJgXTAAaN27M3LlzWbhwIeHh4bz22mt89NFHFe623psv7Yz9mjdJNGrhn3WAk0uLNo6KiIiIXBiLYRiGq4s4n9TUVPz9/UlJSSmz60cAnE6D9yY9xr0np5Bj9cD9zkUQ0q7MPk9ERKQqK+rf7wp5N42rWK0Wgi65le3Ourg7s3BMHwgJf7q6LBERkSpNYeRf+nZqxuOeT5JqeGM7eQRj+kA4FO/qskRERKoshZF/8fd255Gh/bgi+xXWOxthcWTh/PZODYomIiJSRhRGCtClUS2eu/kKHrA8xjGjBtakOFj7mavLEhERqZIURgpxeZtgRg24hMm5gwAw5o+D5I0urkpERKTqURg5h+s712dt6GDWOZthcebAlrmuLklERKTKURg5B6vVwsheLfjecbE5Y/EL8E4EfH497PrdtcWJiIhUEQoj53FluxD2N76BNc7m5owj22H7rzB9AGQccW1xIiIiVYDCyHlYLBa6tKjPDdkT6Jn9NjPrjT+z8KPLYPdS1xUnIiJSBSiMFMENEWHU8fUiwRnIuB3tGZH9CKlWfzi2C6b1hyM7XF2iiIhIpaUwUgS1fOz8PKYnbw/pxNUdQlns7MRVJ58502H5264rTkREpJJTGCmimj52BobX5Y3BHenZog77jCAey7nDXLhmGnw1BFL2ubRGERGRykhhpJjcbVY+u60r/+3WgK8dvZmcOxCnxQ3i58GHfeDEUVeXKCIiUqkojJTQ89e0J6JRbSbl/ocXwj7A8K4NGQdh12+uLk1ERKRSURi5ALdf0hiAj7d68Zvt1FgkKz7Ug/VERESKQWHkAvRtG8KT/Vtjt1n58kgzc2bCcngvCvaudG1xIiIilYTCyAWwWCzc0aMJc0ZfTKy9G/dm38eJGo3AcMDHl0PyJleXKCIiUuEpjJSCtnX96dMymJ+cUcxu/fqZBe9HwcHNritMRESkElAYKSVt6voB8OQfmbzf9H2M2qeGj//hft1hIyIicg4KI6Xk+s71aRnsC8DLG/2Z2exlsPvCvpXwfndITXRxhSIiIhWTwkgpqePrwYIxPRjRvREAj/2WSezFH4B3IKQdgB8fcG2BIiIiFZTCSCmyWCw8fEVLfD3dABiyAHb0eM1cuHOJTteIiIgUQGGklPl4uBHzUC8ualgTgMf/CsIZ1BYcWfDbyy6uTkREpOJRGCkDQb6eTLoxHC93Gyt2H+PlzGvNBRtmg2G4tjgREZEKRmGkjDQO9OHdoZ2wu1mZdbA+DoubOVz8n++5ujQREZEKRWGkDF3WOphxV7biKH68arnVnLnoeY09IiIi8g8KI2VsaNcGeLhZef9EH1IDO0POCZgxFJwOV5cmIiJSISiMlDEvu42hkQ0AC//ndj9Y3eHoTlj5oatLExERqRAURsrBNR3rAbAw0ZvcHo+YM39+HLb/6sKqREREKgaFkXLQLKgG/l7upGXl8k7ONdD2WjCcMHM4LH4Rkta7ukQRERGXURgpBz4ebjx3TTsAPvxjF/u7vwC1mkJOhjn2yEfRcPK4a4sUERFxEYWRcjKgQyjh9f05mePgwR8SYNQyGPgueNeG3Ew4sNbVJYqIiLiEwkg5sVgsPDvIPDqycvdR1hzIhM43Q4t+Zoc/XteAaCIiUi0pjJSj8LAAolsHA/DIrL9IzcyB3uPMO2x2/wGrPnJxhSIiIuVPYaScvXR9e0L9Pdl5OINnf9wEAQ2g+RXmwnkPw4E4l9YnIiJS3hRGyllgDQ+evroNAN+s2ceVb/7O9LpPQq0mZoepl0LMs7qgVUREqg2FERe4sl0I91/WHDerhS1JaUyYv4ukwfOgzSAwHPDHa/BOBCT+7epSRUREypzCiAtYLBbGXt6ClU9E07C2NwDrj1jhps/guqlQszGcOAyfXgXx83Vhq4iIVGkKIy5Uy8fOxU1rAzB/faI5s8NNcEcMeNWC7DT46j/waT9IP+TCSkVERMqOwoiL3RARBsDsdftZm3DMnOlTG0bGQNS95nRCLMx/xEUVioiIlC2FEReLaFiTjmEBANw2bRVZuaee5lurCfR9AW6dBxYbbJwDG79zWZ0iIiJlRWGkAnjtpnAAjp/IYfjHK0lMOXlmYaPu0P1+sx072QXViYiIlC2FkQqgaZ0afHFHJN52Gyt2HeXqt5fyxYo9Zzp0HGa+71sJS98AR65rChURESkDCiMVRPdmgcy9vwdN6/hwJCObJ+ZsYPry3ebC2s2gWbTZ/vX/YPUnripTRESk1CmMVCCNA32Ye38P7u7VFIAJP2zkrV+3keVwwpCZ0HqA2XHhU3B4uwsrFRERKT0KIxWMp7uNh69oQXh9fwDe+HUrV7zxO/tSs+HG6dC4p/mU3yUvurhSERGR0qEwUgG52azMvCuK128Kx9tuY8+RE7w0fwtYbXDpU2anDd/CkR2uLVRERKQUKIxUUJ7uNq7rXJ8Zd3bDYoGf/k4kdscRCOsKTXqbnX55SqOziohIpacwUsF1qB/A0K4NAHj/t1NHQjrfYr7Hz4WYZ8DpcFF1IiIiF05hpBIY3MUcpTV2x2GOZmRDu+sgfKi5cOkb5kP19sS6sEIREZGSUxipBJoH+RLg7U6Ow6DXpMXsPJQOA96E3o+bHY7tgjl3afwRERGplBRGKgEvu41pI7oCkJaZy0vzt2DY7NB7HIxZD/YacHwPxL7j4kpFRESKT2GkkugYFsDsey7GzWrhl03JLNl66im+AQ3gkjFm+7dJkJ3hshpFRERKQmGkEuncoCbDoxoB8P7iHRin76Tp/iB414acDFj9qesKFBERKQGFkUpmeFRD3KwWVu4+yraD6eZMmxv0Gme2YyeD0+m6AkVERIpJYaSSaRTowyXNAwH48PedZxZE3GpeO5J2AP6e4ZriRERESqBEYWTy5Mk0atQIT09PIiMjWbly5Tn7v/nmm7Rs2RIvLy/CwsJ48MEHyczMLFHBAld3qAvAN2v2MX99ojnTzQOaX262v7sHju5yUXUiIiLFU+wwMnPmTMaOHcuECRNYu3Yt4eHh9O3bl4MHDxbY/8svv+Sxxx5jwoQJbN68mY8//piZM2fy+OOPX3Dx1dX1nevRvVltAO77ah3P/7SJ1MwcuOq1Uz0M+Osr1xUoIiJSDMUOI6+//jojR45kxIgRtGnThilTpuDt7c0nnxT8WPvly5fTvXt3hg4dSqNGjbjiiisYMmTIeY+mSOEsFgtvDO5IZONa5DoNPlq6i/Gz14NPbej+gNnpj9fgULxrCxURESmCYoWR7Oxs1qxZQ3R09JkNWK1ER0cTG1vwCKAXX3wxa9asyQsfO3fuZN68eVx11VWFfk5WVhapqan5XpJfkK8nM+7sxuNXtQJg3vpEdhxKhz5PQFAbcObCnLshR6fDRESkYitWGDl8+DAOh4Pg4OB884ODg0lKSipwnaFDh/Lss89yySWX4O7uTtOmTendu/c5T9NMnDgRf3//vFdYWFhxyqw2LBYLd/ZsSrcmtTAMeGHuZvPakd7jzQ4H1sJvL7u2SBERkfMo87tplixZwosvvsh7773H2rVrmT17NnPnzuW5554rdJ3x48eTkpKS99q7d29Zl1mpPdLXPDqyaMtBklMzodXV0O4Gc+GORS6sTERE5PyKFUYCAwOx2WwkJyfnm5+cnExISEiB6zz11FPcfPPN3HHHHbRv355rr72WF198kYkTJ+IsZDwMDw8P/Pz88r2kcBENa9Iy2BeAT5ftBqsVLr7XXJgYB3+87rLaREREzqdYYcRutxMREUFMTEzePKfTSUxMDFFRUQWuc+LECazW/B9js9kAzowgKhfswcubAzBr9V6cTgPqdjLHHgGIeQY2zgH9vEVEpAIq9mmasWPHMnXqVKZPn87mzZsZNWoUGRkZjBgxAoDhw4czfvz4vP4DBgzg/fffZ8aMGezatYuFCxfy1FNPMWDAgLxQIhfustbB+NhtHMnIZun2w+bMfq+A/6nrbWbdCgseUyAREZEKx624KwwePJhDhw7x9NNPk5SURMeOHVmwYEHeRa0JCQn5joQ8+eSTWCwWnnzySfbv30+dOnUYMGAAL7zwQuntheBus9KtSW1ithxk5Ger+WNcH4J8PeGOXyH2XVj+DqyYAntXwA2fQK0mri5ZREQEAItRCc6VpKam4u/vT0pKiq4fOYe9R0/Q45XFALw3rDNXtQ89s/DnJ8xQAhDSAe76HSwWF1QpIiLVRVH/fuvZNFVIWC1v+p8KIC/N30KO4x8XCF/xPNw612wn/a0RWkVEpMJQGKliHry8BTarhYSjJ5ixMuHMAosFGl0CIe3N6e9Gwf+ug8PbXVOoiIjIKQojVUyzoBqM6tUUgOd+2syWpH+NXvvf2dDmGrO9Iwa+vBEKucVaRESkPCiMVEGDu4RhtUC2w8knS//19N4aQXDTdLjnT7DXgKM74ZcnXVOoiIgICiNVUlgtbyZeZ56OWbX7WMHjuQS1hsufNdt/TtbpGhERcRmFkSoqomFNAHYdzuCVnwt5em+X2yGwpdne/Xs5VSYiIpKfwkgV1SzIl/subQbAB7/tYG3CsYI7dvqv+T7/MVj/TTlVJyIicobCSBU2uk8z6vh64DTguveW8+WKhLM7dR0JAQ3BkQXzHwVHbvkXKiIi1ZrCSBXm6W7j67ui8sYeeXzOet6O2UbuP8cfcfc6NQCaFU4cge0LXVStiIhUVwojVVzjQB/eGdKJS5oFAvD6wq38348b83fyCoCLbjfby98p3wJFRKTaUxipBqxWC+8M6cQtUQ0B+PzPBE5mO/J3irzLfN+zDOaMAkdOOVcpIiLVlcJINVHTx84zg9rh62k+G/HGD5ZzMC3zTIfA5tBlpNn+60uYc7ee8CsiIuVCYaSaGdnDfFrvhv2pXPLSYr6P239mYf9X4cqXzPaGb2DZWy6oUEREqhuFkWrm/sua8+XISFqF+JLtcPLknA1k5vzjlE23UdB7vNn+dQLEz3dNoSIiUm0ojFRDFzcNZO79PQj19yQtK5eHZ/2Vf5TWXuMgtKPZ/uo/sHCCS+oUEZHqQWGkmrJZLTx/TTvcrBZ++juRbQfTzyy0WGDw/6B5X3N62Zuwf61L6hQRkapPYaQau6x1MM2DfQFYvy8l/8KABjDsa2jcy5ye1h/mPgwnj5dvkSIiUuUpjFRzfdsGA/DKz1uIT0o7u8O1UyCoDeScgFVTYfEL5VyhiIhUdQoj1dytFzfCz9ON5NQs7v2ygFMxfnXh7qVwyVhz+u+vwek8u5+IiEgJKYxUcwHedqbf1hWAbQfT2ZpcwNERqw36PA5uXpB5HP58r3yLFBGRKk1hROjUoCa9WtQB4Me/DhTcyeYOPR4y2wufhl1/lFN1IiJS1SmMCAADwusCMH35btKzCnlyb8+HofkVYDjgh/sg+0Q5VigiIlWVwogAcHWHUHw93UjNzKX3pCUkppw8u5PFAle/aZ6uObYL1n9d7nWKiEjVozAiAHi623hvWGcADqdnMfKz1WTlOs7u6F8Put9vtjf/VI4ViohIVaUwInl6NK/DT/ddgqe7lQ37U3ntl60Fd2w9wHzf/issfwecBYQWERGRIlIYkXza1fPnuUHtAPjw9538vvXQ2Z2C20GbQYABvzwJOxaXb5EiIlKlKIzIWW68KIyIhjUBGP7JSj76Y2f+DhYLXPcR1GpqTq+aChmHy7lKERGpKhRGpEBv3NSRro1rAfDCvM1nX9DqZod+L5vtrQvgnc5wYF05VykiIlWBwogUqEFtb76+K4pmQTUwDBj52WpidxzJ36n55XDjdPAJgswU+LA3fDcadv3ukppFRKRyUhiRc3r5+g7YbeYFrUOm/sk3a/bl79D2Ghg2Cxr3NKfjPofpA+CLGyE1sdzrFRGRykdhRM4pomFNfn6wJy1PPd33rOtHAOp2hFt+hOHfQ8f/gsUK236B2SPBMMq3YBERqXQURuS8Ggf68Nnt5vNr4pPTOJaRXXDHJr3hmslwR4w5MNruP2DjnPIrVEREKiWFESmSYD/PvOtH7p+xDqfzHEc86nWGrneY7W0Ly6dAERGptBRGpMjeHNwRDzcrf2w7zOcr9py7c+Pe5vvG2XBwS1mXJiIilZjCiBRZu3r+3HZJYwBe+2UrR9KzCu/c9FLztE1uJnxzm64dERGRQimMSLHcd2kzQvw8STmZw49/HSi8o9UKA9422wc3wuFChpYXEZFqT2FEisXb7sbNUQ0BWBxfwFDx/1SzIdS7yGz/OAYcOWVbnIiIVEoKI1JsvVrUAWD5jsP8sjHp3J0Hvg32GpCwHL4bVQ7ViYhIZaMwIsXWtq4f3ZvVJsdhcOf/1rBq99HCOwe3hes/MtvrZ+liVhEROYvCiBSbxWLh/f9GEODtDsCf/x4m/t9a9oOml5ntTd+VbXEiIlLpKIxIifh5unPfpc0BeG3hVh779u9zjz3S7jrzfdlbcOIcR1JERKTaURiREru5W0P6tQsBYMaqvQz96M+zn+57WrvrwTsQck7AWx1hw7flV6iIiFRoCiNSYnY3K+//N4LRfZoC8OfOo1zy8mLumL6KvUdP5O/s7gWDPzeHic9KMcce2TLXBVWLiEhFozAiF+yRvq2Yc8/FdAwLwOE0+HXzQS59bQkxm5Pzd2wYBY/ugMCW5vSMoTD/MchKL/+iRUSkwrAYRsUfGjM1NRV/f39SUlLw8/NzdTlyDluSUhk/ez3rEo4T4udJzEO98PFwy9/pxFGY/6h5dw2AzQ6Ne0LfiVCnRfkXLSIiZaKof791ZERKVasQP74a2Y2wWl4kpWYya/Xeszt51zJv973yZajZGBzZsP1XWPZmudcrIiKupzAipc7T3cbgi8IAWLXnWOEdu90N96+Da943p+Pnw6bv9RwbEZFqRmFEykREw1oA/LIxifiktMI7WizQeiAEtoCTR+Hr4TB7JORml1OlIiLiagojUia6NalF75Z1yHEYPDQr7tydPWrAnb9Br3FgdTevJZl1q46QiIhUEwojUiYsFgvPDmwHwIb9qRxMzTz3CnZv6PM4DJkBWCB+LmQcLvtCRUTE5RRGpMw0qO1N86AaALz2y9airdQ8Gvzqme0tP5ZRZSIiUpEojEiZeqxfKwBmrt7LzFUJRVup1VXm+08PQvyCMqpMREQqCoURKVOXtQ7mrl5NAJjww0bSMnPOv9IVz5u3/AL8/DhsmQdOZxlWKSIirqQwImVuXN9WhPh5kpnjZNGWg+dfwc0DrnzJbB/dATOGwLe3K5CIiFRRCiNS5qxWC4M61gXgh7gDRVup5ZVwxyLzAXsAG2fDwqfKqEIREXElhREpFwPCzTASs+UgP/5VxEBSPwKu/xi63WNOx74Lb4XDzJth269lVKmIiJQ3hREpF+3q+dO7ZR0A7vtqHY99+3fRVrRYoO+LEP0MWN3g2G7Y/AN8cQPsiS27gkVEpNyUKIxMnjyZRo0a4enpSWRkJCtXrjxn/+PHjzN69GhCQ0Px8PCgRYsWzJs3r0QFS+X15uCOXNEmGIAZq/byxsKt5DqKcB2IxQKXjIGHt8HwH6BuJ8CAeQ9DblaZ1iwiImWv2GFk5syZjB07lgkTJrB27VrCw8Pp27cvBw8WfGFidnY2l19+Obt37+abb74hPj6eqVOnUq9evQsuXiqXAG87Hw6/iP90MZ9b81bMNkZ9sZbk8w2Idpp3LWjSC4bOAu/akLwB5twFjtwyrFpERMqaxTCKN+Z2ZGQkXbp04d133wXA6XQSFhbGfffdx2OPPXZW/ylTpjBp0iS2bNmCu7t7iYos6iOIpXLIdTj58I+dvLIgHjAPfEQ0qMlz17SjdWgRv9/NP8HMYWa7x0Nw2dNlVK2IiJRUUf9+F+vISHZ2NmvWrCE6OvrMBqxWoqOjiY0t+Pz9Dz/8QFRUFKNHjyY4OJh27drx4osv4nA4Cv2crKwsUlNT872k6nCzWbmndzMmD+1M5wYBGAas3nOMqX/sLPpGWl8Nvceb7T9eg8QiXoMiIiIVTrHCyOHDh3E4HAQHB+ebHxwcTFJSUoHr7Ny5k2+++QaHw8G8efN46qmneO2113j++ecL/ZyJEyfi7++f9woLCytOmVJJ9O8Qyux7uvP2kE4AzFm3n40HUoq+gZ6PQMPuZjvmWT1YT0Skkirzu2mcTidBQUF8+OGHREREMHjwYJ544gmmTJlS6Drjx48nJSUl77V3796yLlNcqFeLOoT4eWIYcM3kZXzw246irWi1weXPmU/63b4QfrxfA6OJiFRCxQojgYGB2Gw2kpOT881PTk4mJCSkwHVCQ0Np0aIFNpstb17r1q1JSkoiOzu7wHU8PDzw8/PL95Kqy9/LnXkP9OCyVkHkOAwmzt/Cu4u2FW3l+hHQ59TpmrWfwc7FZVeoiIiUiWKFEbvdTkREBDExMXnznE4nMTExREVFFbhO9+7d2b59O85//B/r1q1bCQ0NxW63l7BsqWpq+dj5+NYuPNK3JQBv/LqN9Kwi3iXT4yFoe63Z3nvu28xFRKTiKfZpmrFjxzJ16lSmT5/O5s2bGTVqFBkZGYwYMQKA4cOHM378+Lz+o0aN4ujRozzwwANs3bqVuXPn8uKLLzJ69OjS2wupMkb3aUaAtzsOp0HCkRNFXzGsm/m+5lNdzCoiUskUO4wMHjyYV199laeffpqOHTsSFxfHggUL8i5qTUhIIDExMa9/WFgYP//8M6tWraJDhw7cf//9PPDAAwXeBiwC0DLYF4AxM9eRlFLEMUg6DQOfIEhPhg96wqbvy7BCEREpTcUeZ8QVNM5I9fL3vuMMfHcZAI/0bcnoPs2KtuLBzTDvEdj9BwQ0hDE6QiIi4kplMs6ISHnoUD+AEd0bATDp53ju/Gw1e48W4ZRNUGvzwXoAx/fA/rVlV6SIiJQahRGpkO7u1ZQO9f0B+GVTMoMmL+NYRsF3X+XjGwytB5jtGUMhtwjriIiISymMSIUU7OfJ96O787/bu+JmtXA0I5tRX6wh5WTO+VceNBk8/CEtEbYuKPtiRUTkgiiMSIVlsVjo0bwOU2+5CLvNyp87jxL+zC/0eXUJHy/dVfiKnv4QMdxs//p/epCeiEgFpzAiFV6flkF8OqILTQJ9ANh1OIPnftrEr5uSC1+p+xhw94ajO2D+oxoqXkSkAlMYkUqhe7NAFj3cm3VPXc4VbczbyF+ct5nMnEIeuOgTCJF3me3VH8PfM8upUhERKS6FEalUavrYmXRjOEG+Huw8nMFj3/6Nw1nIUY/LJkCnm832/HGQc7L8ChURkSJTGJFKx9/LneeuaQfAd3EH+KSw60csFuj/GvjVh8zjEPNc+RUpIiJFpjAilVLftiHc2bMJAK/8vIXdhzMK7ujmAb0eNdtxn+upviIiFZDCiFRaYy9vQasQX3IcBp8uO8fdNR2HgpsXZKbAn++VX4EiIlIkCiNSaXm623jgsuYATI/dw/z1iQV3tLlDj7Fm+5cnIGlDOVUoIiJFoTAilVq/9qH0bWveXXPvV+vYeSi94I69HoVaTc321Eth6Rs6ZSMiUkEojEil9/pNHWkV4ovDabBoy8HCOw6bBY17giPLHAxtzSflVqOIiBROYUQqPR8PN/q0CgLg730phXes3RSG/wBd7jCnt/1aDtWJiMj5KIxIldC3bQgAP/x1gAdnxpGRVcgQ8BaLeUErmM+tWTAeju0unyJFRKRACiNSJXQMC2BYZAMA5qzbz4QfNhbeuW5naNIbMMy7a6ZeBofiy6VOERE5m8KIVBkvXNueNwaHAzB/fSJGYc+jsVhg2Ddw0/8gqC2cOAzfj4bc7HKsVkRETlMYkSrlyrah1PBwIyPbwRcrEgrvaHOHNgNh8P/AzRP2rYJ5D5dfoSIikkdhRKoUL7uNUb3NW3if/WkTxzLOc7SjdlOI/j+zvXa6eQ1JViG3B4uISJlQGJEq557eTWlax4fsXCffxe0//wrdRkG7G8z2n+/BBz01BomISDlSGJEqx2KxcH1EfQAmL95BVq7j/Ctd8775UD2Aozvgy5sUSEREyonCiFRJt3VvjK+nG4fTs3js2/U4nYVczHqam90cf+S6j8zp7Qth129lX6iIiCiMSNXk6W7jnt7NAPNW3yKdrgHocCO0u95sr/2sjKoTEZF/UhiRKmtU76bc3K0hAM/9tKlop2sAuow03zfOhpVTy6g6ERE5TWFEqrQn+rcmyNeDYydy+PC3nUVbqWHUmaMjvzwFjpyyK1BERBRGpGrzdLdxZ88mAExbvrvwgdD+7bqPwN0Hck/CgXVlWKGIiCiMSJX3324NsVktHMnI5reth4q2ktUKTXqZ7XWfl11xIiKiMCJVn6e7jf7tQwG4839ryMwp4rUjLa8y39dOhy1zy6g6ERFRGJFqYdKNHbC7WcnOdXLde8vZlpx2/pWaXgru3mZ72dtlW6CISDWmMCLVgoebjTcHd6SGhxubElMZ8O5S1iYcO/dK/vXg9oVm+8A6KOr1JiIiUiwKI1JtXNU+lIVje9Kunh+ZOU6+OteD9E6r09J8kJ4jC5a9VfZFiohUQwojUq2E+ntxb5/mACzfcYQj6VnnXsHmDl1PjTuy9HXITC3jCkVEqh+FEal2Lm5WG2+7jf3HT9L1xRheXrDl3Lf8Xvo0+NWDzBSYO7b8ChURqSYURqTa8fN0Z8p/I2gd6ofDafD+kh1MW7678BXc7HDNe2Z7/SxY8WG51CkiUl0ojEi11LNFHeY/0IOHr2gBwKSf49l1OKPwFRr3gnY3mO35j0DsZMjJLIdKRUSqPoURqdbu7tWUViG+nMh2cMsnKwt/fo3FAtd/BM37mtM/Pw6fDYSck+VXrIhIFaUwItWam83KqzeGU8PDjYSjJ7jyzT9Ys6eQW34tFrjpM4h+xhwqfu8KiPuifAsWEamCFEak2mtXz5+Xrm+P1QK7Dmcw5MM/SUop5BSMuydcMgaiRpvTa6aBI7e8ShURqZIURkSAqzvU5Y9xl2KxQLbDyeAPY8nOdRa+QqdhYK8BSeth9cflV6iISBWkMCJySr0ALxY80BOAPUdOMHf9gcI712wE3ceY7fnjYO+qMq9PRKSqUhgR+YeWIb7c1r0xAN/HnSOMAPR4CNpeCxiw7M0yr01EpKpSGBH5lxsvqg/AkvhDvPpzfOEdrVboPd5sb5kLyRvLoToRkapHYUTkX1qH+nFnzyYAfPD7jsJv9wXz2TWtrgYM+PmJ8ilQRKSKURgRKcD4fq0I8HYnx2Gwfl/KuTtHP2O+7/odMg6XfXEiIlWMwohIASwWC10b1QLg27X7z905sBmEdADDAZt/KIfqRESqFoURkUJEtw4G4KuVCfy29dC5O3e4yXxf9DwcWFfGlYmIVC0KIyKFuLhZbew28z+RyYu3n7vzRbeBb104cQSmDYCMI+VQoYhI1aAwIlKI+jW9mXV3FACbE1MxDKPwznYfuG2BOUx8dhrsXFxOVYqIVH4KIyLn0CrUF5vVQlpmLt/FnefakZoNoesdZnvbwrIvTkSkilAYETkHDzcbA8PrAjBpQTwO5zmOjsCZp/qunwWbf4RzHU0RERFAYUTkvJ4d1BZfDzcOpGTy1q9bz925QRSEdTPvrJn5X3i9DexcUi51iohUVgojIufh6+nOk1e3BuDtRdtZm3Cs8M5WK9zyA3T6rzmddgC+/I85QquIiBRIYUSkCK7vXJ8WwTUAGDr1T7YfTCu8s5sHDJoMjyVAk96QexK+Hg4Ln4aTx8ulXhGRykRhRKQI3GxWZt19MWG1vMjMcXL5G7/z5vlO2Xj6w39nQ7vrwZkLy96C96Jg68/lU7SISCWhMCJSRP5e7nx6a1ciG9fCMODNX7exJP7guVey2uD6j2HIDAhoaJ62+eY2yM4on6JFRCoBhRGRYmgWVIOZd0VxXad6APwQd+D8K1ks0LIfjF4BPnUgOx3Wf1PGlYqIVB4lCiOTJ0+mUaNGeHp6EhkZycqVK4u03owZM7BYLFxzzTUl+ViRCmPQqTAye91+/t53vGgruXtB9wfM9vK3ISezbIoTEalkih1GZs6cydixY5kwYQJr164lPDycvn37cvDguQ9X7969m4cffpgePXqUuFiRiqJXizpc1ioIgLv+t4bjJ7KLtmLnW8xrSY5sh1m3wMlz3JkjIlJNFDuMvP7664wcOZIRI0bQpk0bpkyZgre3N5988kmh6zgcDoYNG8YzzzxDkyZNLqhgkYri+Wvb0ai2N4kpmXwWu6doK3n6wcB3zfbWBfBOBGyPKbsiRUQqgWKFkezsbNasWUN0dPSZDVitREdHExsbW+h6zz77LEFBQdx+++1F+pysrCxSU1PzvUQqmlB/Lx68vAUAn8XuJjPHUbQV2wyEm7+DOq3MB+vNH1d2RYqIVALFCiOHDx/G4XAQHBycb35wcDBJSUkFrrN06VI+/vhjpk6dWuTPmThxIv7+/nmvsLCw4pQpUm6uah9KqL8nh9Oz+fD3nUVfsWkfGP6D2T6yDX59RkPHi0i1VaZ306SlpXHzzTczdepUAgMDi7ze+PHjSUlJyXvt3bu3DKsUKTl3m5Ux0c0B+N+fe3Ce79k1/+QbDBfdZraXvg7rPi+DCkVEKj634nQODAzEZrORnJycb35ycjIhISFn9d+xYwe7d+9mwIABefOcTqf5wW5uxMfH07Rp07PW8/DwwMPDozilibjMNZ3q8fxPmzmUlsW6vceIaFir6Cv3fx1yTsJfX8HchyCwBTSILLtiRUQqoGIdGbHb7URERBATc+aCO6fTSUxMDFFRUWf1b9WqFevXrycuLi7vNXDgQPr06UNcXJxOv0iV4OFm49LW5p01X60s5lE8i8UMJPUiwJEFn1wBCx7XsPEiUq0U+zTN2LFjmTp1KtOnT2fz5s2MGjWKjIwMRowYAcDw4cMZP348AJ6enrRr1y7fKyAgAF9fX9q1a4fdbi/dvRFxkeFRDQH4Zs0+npizHqM413/YvWHoLAjtaE7/ORkmR8KBuFKvU0SkIip2GBk8eDCvvvoqTz/9NB07diQuLo4FCxbkXdSakJBAYmJiqRcqUpFFNKzFbd0bA/DFigQ2HijmHWA+tWHkIjOU+IZCehJ82BumD4Bdf5R+wSIiFYjFKNb/wrlGamoq/v7+pKSk4Ofn5+pyRAp188cr+GPbYcZd2YpRvc++HqpIUvbDt3dAwnJzOrg9jFpaekWKiJSTov791rNpREpRdGvzCOH8DYnFO1XzT/714Lb5Z279TV4P8x6BrLRSqlJEpGJRGBEpRVe0DcZigb/3pXDbtFUsiT9Y8lDSpBe0utpsr/wQpl0Nqz6GtORzryciUskojIiUolB/L54b1A6rBRbHH+LWT1fx+JwNxRt/5J9u+gwGvQdunpAYB3PHwpvtdHGriFQpCiMipey/3Roy484ormpvjr3z1coExsyMIzvXWfyNWW3QaRjcuQQuedCc58iGqZfCsd2lVrOIiCvpAlaRMvR93H4e+vovcp0G9Wt6MWFAWy5vE3z+FQtzeBu8390ck8TTHx7eBm4aIFBEKiZdwCpSAQzqWI/3hnUmsIYH+46dZORnqxk/e33JNxjYHO5bbbYzU2BSM9izvHSKFRFxEYURkTJ2RdsQ/ni0DyN7mOOQfLUygb1HT5R8gwENoN8k8K4NWakwawQ4ckupWhGR8qcwIlIOvOw2nujfhu7NagPw+Z97LmyDkXfCmA3g5mUOkPb3zFKoUkTENRRGRMrR6VFaP/9zD7sOZ1zYxuze0O46sz1/HKQfvMDqRERcQ2FEpBz1aRlE+3r+ZGQ7uOH95SzclFzycUgA+jwO7j6QnQZf3AC5WaVXrIhIOVEYESlHVquFyUM7E+LnyZGMbEZ+tppXf4kv+Qb968Mdv4KHPyT+BX99VXrFioiUE4URkXLWoLY3MQ/1yjtl81nsnpKNQXJacBvzGhKAJS/pYlYRqXQURkRcwMfDjYf7tsDbbiMtM5cB7ywlNTOn5Bvsdg9YbJCWCPHzSq9QEZFyoDAi4iLedjfeGNwRgPjkNF74afMFbKwWdBtltld8ABV/LEMRkTwKIyIu1LdtCK/fFA7AzNV7+b8fNpLrKOEpm+ZXmO97lsJXQyA1sZSqFBEpWwojIi52Xef63N2rKQDTlu+m96tLWLylBLfpNu4JbQaZ7a3z4Z0Ijc4qIpWCwohIBTDuypY80rclNb3d2XfsJCOmrWLSz1s4kl6MW3UtFrjhU+j/OgQ0hJwM+OUpcDrKrnARkVKgB+WJVCAZWbk8+u3fzP3bPMViscADlzVnTHSL4m3o+F6Y3BVyTsBlT0OPh8qgWhGRc9OD8kQqIR8PN97+TyeeGdiWmt7uGAZMX76bnOJeRxIQBle9arZjnoNFz0NudukXLCJSChRGRCoYm9XCLRc34vdH+2CzWjh2IoeoiYtYl3CseBsKHwJtrgEM+H0SzH+kLMoVEblgCiMiFZSvpzsvX98BH7uNw+lZXPf+cm54fzmz1+4r2gasVrhxGlw31ZxeMw3Wf1NW5YqIlJjCiEgFdkNEfZaOu5SIhjUxDFi95xhjv/6LNxZuLdoGLBbocBO06GdOf3s7LH0THBcwwJqISClTGBGp4Gr62Pnm7igWP9yblsG+ALwVs40vVuwp+kaueQ/a3WC2f50A73WDY7tLv1gRkRJQGBGpBCwWC40DfVgwpgcDw+sC8MScDXyzZh9OZxFuiPOuBdd/ZF7U6ukPR7bDb6+UcdUiIkWjMCJSiVgsFt76T0dC/T0BeHjWX7wwr4jDyFss0HUkDHzXnI77Ar4dCSeOllG1IiJFozAiUslYLBa+viuKOy4xn/r78dJdJKacLPoGWlwJLa8y2+u/hp/GlH6RIiLFoDAiUgmF1fLmyavb0CTQB4BJC+KLvrKbHYZ8BYO/MKc3fQ8LJ0DyxjKoVETk/BRGRCqxO3o0AWD2uv18vXpv8VZufTX0fNRsL3sT3r8YNs4p3QJFRIpAYUSkEhvcJYz+HUIBGD97PdsPphdvA30eh0HvnZle8DhkHCnFCkVEzk9hRKQSs1ktvDukE10b1cLhNHjtl2KcrgHzotZOw+C2X8zptAMwqQnMGqHh40Wk3CiMiFRyFouF+y5rBsDPG5PYd+xE8TfSIBL6ToSgNub0xtnwciPYv7b0ChURKYTCiEgV0KN5HTo3CMBpwLRlu0u2kah74J5YuOJ5czonA2beDBmHS61OEZGCKIyIVBHDIhsC8NHSXcTtPV7yDV18H4xYYLZT98GXN0FOMW4dFhEpJoURkSqiT6sgfD3cALhm8jLu+WINxzJKeN1Hwyi4+k2zvX8NvN0ZMlNLp1ARkX9RGBGpImr52Pn+3u60CjGfXzNvfRKv/FzMC1r/6aIRZwJJ2gFI/OvCixQRKYDCiEgV0qRODRaM6ckbg8MBmLkqgWXbL+Caj4tGQOOeZnv522AU4Tk4IiLFpDAiUgVd26k+N11UH6cBD86MY+WuC3j+TI+HzPdtv8C6/5VOgSIi/6AwIlJFPXV1G+oFeHEwLYubPojlq5UJJdtQk97Q8xGz/eeUUqtPROQ0hRGRKsrX050f77uE6zrXA+CJOevZmpxWso1ddBtY3eDgRvhtkk7XiEipUhgRqcJq+dh57cZw2tb1w2nA1W8vZdGW5OJvyK8utLrabC9+Hj7sDTHPwtZf4OSxUq1ZRKofi2FU/P/FSU1Nxd/fn5SUFPz8/FxdjkilE5+UxpiZcWxONG/P7dE8kLf/04maPvaib+TkMVj+Lix7C5w5+ZfVi4DwIdDlDnOIeRERiv73W2FEpJrIzHHw4rzNfP7nHpwG3BLVkGcGtSv+ho4nwI7FsHcFJPwJR3ecWdb1Tuj3igKJiAAKIyJSiK9X7+XRb/4GYFhkA8b1a4Wfp3vJN5iaCHFfwKLnAcO8++ayp0unWBGp1Ir691vXjIhUMwPD69KlUU0AvliRwBNzNlzYBv1CoefD5hERgGVvgyP3AqsUkepEYUSkmvF0t/H1XVF5A6P9+NcB/u+HjazefZQLOlB60W3g4WdeTzL7jlKqVkSqA4URkWrIYrFwbaf6jIluDsC05bu5YUoss9fuL/lGbW5njo5snAMrPiyFSkWkOlAYEanGHrisOZOHds6bfmjWX6zefZRch7NkG+w4BLo/YLbnPwJzH9KYJCJyXgojItWYxWKhf4dQFozpkTfvhimxXP7G7xxMzSzZRi99GiJHme1VH8GcuyGrhIOtiUi1oDAiIrQK8eODmyO4qn0IALsOZ9D1xRi+jyvBaRubG/R7CfpOBCzw9wyYdrUuahWRQimMiAgAfduG8N6wCD67rSt2N/NXw5iZcSzfUcKn/kbdA7f+ZA4jnxgHC8ZBzsnSK1hEqgyFERHJp2eLOqx8/DLsNiuGAaM+X0vM5hIMIQ/Q6BK4/Fmzveoj+Px6XUMiImdRGBGRswR421k4tiee7lZSTuZw+/TV7D16omQbixoNN/3PbO9ZBr88CZkppVesiFR6CiMiUqCGtX2Ieah33vQf20p4ugagzUDoMNhsx74LrzSFWbfqwlYRARRGROQc6gV4cVv3xgA899Mmdh/OKPnG+k6Eno9AYEtzYLSNc+DX/wNnCW8jFpEqQ2FERM5p7BUtaBFcg5M5Dl75eUvJN+RTGy59Eu5dCZc/Z85b9RG8FwkbvoWUfbqeRKSaUhgRkXOq4eHGC9e2B2De+iTik0rh1ErUvXDZBPD0h8Nb4Zvb4I22MLkrxC9QKBGpZkoURiZPnkyjRo3w9PQkMjKSlStXFtp36tSp9OjRg5o1a1KzZk2io6PP2V9EKp6LGtakbV3ziZt93/ydifM2s2z7YU5kl3DsEKsVeoyF+9ZClzsguB1YbGYw+WowfHs7pB4oxT0QkYrMYhTzyVgzZ85k+PDhTJkyhcjISN58801mzZpFfHw8QUFBZ/UfNmwY3bt35+KLL8bT05OXX36ZOXPmsHHjRurVq1ekzyzqI4hFpOys2XOUCT9sZMP+1Lx5blYLbev506VhTaLbBNOtSe2Sf0BmCvzxOix/BwwHuHnBkC+h6aWlUL2IuEJR/34XO4xERkbSpUsX3n33XQCcTidhYWHcd999PPbYY+dd3+FwULNmTd59912GDx9epM9UGBGpGAzD4Ie/DhCz+SCrdh8lMSX/kPGPXtmSe3o3u7AP2bMc5j0CyRsgNBzu+v3CticiLlPUv99uxdlodnY2a9asYfz48XnzrFYr0dHRxMbGFmkbJ06cICcnh1q1ahXaJysri6ysrLzp1NTUQvuKSPmxWCwM6liPQR3No5r7jp1g1e6jLNpyiB//OsArC+JxOg3uvbR5yT+k4cXwny/grY6Q+Bf8NgkuGQM291LZBxGpeIp1zcjhw4dxOBwEBwfnmx8cHExSUlKRtjFu3Djq1q1LdHR0oX0mTpyIv79/3issLKw4ZYpIOalf05trO9XnnSGdGBNtBpC3YraRmeO4sA3XbAS9xpntxc/DlB6wbaFuAxaposr1bpqXXnqJGTNmMGfOHDw9PQvtN378eFJSUvJee/fuLccqRaQk7u3TDD9PN3IcBrd+uhKn8wLviOn9GAx6D7xrw6HN8MUN8FYHWPIyZOpoqUhVUqwwEhgYiM1mIzk5/3MqkpOTCQkJOee6r776Ki+99BK//PILHTp0OGdfDw8P/Pz88r1EpGJzs1mZdGM4AH/uPMrjc9aTlXsBR0gsFug0DO5dDd3uMW8DTtkLS16EN9vBz09ATub5tyMiFV6xwojdbiciIoKYmJi8eU6nk5iYGKKiogpd75VXXuG5555jwYIFXHTRRSWvVkQqtL5tQ3jh2nYAzFi1l0HvLuNIetZ51joP71pw5UR4KB6umwq1m5t33sS+CyumlELVIuJqxT5NM3bsWKZOncr06dPZvHkzo0aNIiMjgxEjRgAwfPjwfBe4vvzyyzz11FN88sknNGrUiKSkJJKSkkhPTy+9vRCRCmNYZEPeGdKJGh5ubElKo99bf/Db1kMXvmF3L+hwE4xeCR3/a85b/w2cOHrh2xYRlyp2GBk8eDCvvvoqTz/9NB07diQuLo4FCxbkXdSakJBAYmJiXv/333+f7OxsbrjhBkJDQ/Ner776aunthYhUKAPC6/LVyG40DvThYFoWIz9bzbbkUnoontUKXUeC1Q2S18M7nWHVxxq1VaQSK/Y4I66gcUZEKqfMHAcD313K1uR0fOw2/m9gW268qJTujtuzHOY+DAc3mtPXvA8dh5bOtkWkVBT177eeTSMiZcbT3cbHt3QhvL4/GdkOHvnmbzbsTymdjTe82BwQrfOpwRMXToCtv5TOtkWkXOnIiIiUOafTIPr139h5OAOA8LAABobXZcTFjbBaLRe28cwU+OhyOBxvTne6Gfq/Dm72C6xaRC6UjoyISIVhtVqYdGM47eqZv4z+2nuc537axOTF27ng/x/y9IeRMebtvwDr/gcf9IT0gxdYtYiUFx0ZEZFylZyaycR5m/kuznwqb4f6/tzZswnt6/kTVtP7wo6UrP4UYp6Bk8egRrB5541XQOkULiLFVmYPynMFhRGRqiU718mL8zbz5coEsnPPDPFeL8CL+y5txg0R9XGzlfDA7eHt8F43cOaAVy0IHwIh7aB+Fwi8gGfmiEixKYyISIV3MDWTt2K2Ebf3ONuS08l2mMEkPCyAsZe3oFeLOiXb8K4/4McH4OiO/PP7vQKRd11g1SJSVAojIlKpZOY4+GJFAm8u3EpaVi4AT/ZvzR09mpRsg44c2DgH9q/JP1JraEe4abr5MD4RKVO6gFVEKhVPdxu3X9KYBQ/2pH+HUACm/LaDlJM5Jdugzd0csbXfyzB+H0Tcas5PjIO5D0FWKQ3CJiIXTGFERCqUegFevHJ9B4L9PDicns3lr//GjJUJ5Dqc51+5MB6+MOAtGDTZnN7+K7zaAtZML52iReSC6DSNiFRIGw+kcPfna9h79CQAof6eXN0hlEEd69Gunn/JNmoY8NcMmDsWck6Y87rcAb0egxolvD5FRAqla0ZEpNLLynXw+Z8JvLNoG8dPnDldE9WkNv/pGsbA8LpYLCW4FTg3G768CXYuNqfDh8C1egKwSGlTGBGRKiMzx8GS+EP8+NcBft6YRK7T/LX18vXtGdylQck2ahjw2yuw5EVzeuRiCOkANrdSqlpEFEZEpErad+wEj8z6m9idR2hax4fPbo+kXoBXyTaWmwUT64Mj25z2rm0+66b34xpOXqQU6G4aEamS6tf05on+rXGzWthxKIM+k5bwfz9sJCvXUfyNuXnA9R9BiyvBwx9OHIGlb8DrreHnJ2DX7+DILf2dEJF8dGRERCqluL3HeXn+FmJ3HgFgdJ+mPNK3Vck36MiF1R/DwqchN/PM/FpNIfJuaNUf/OtdYNUi1YtO04hItTD19528MG8zAL1b1mHCgLY0DvQp+QYzU2DHItj6C2ydbz7nBsDdG/77LTS8uBSqFqkeFEZEpFowDIOXFmzhoz924XAaeLhZeaRvS26/pHHJ7rT5p8xUWDsdVnwIKQlmIBk2CxpdUjrFi1RxCiMiUq3sOJTOyOmr2Xk4A4AB4XV5+z8dLzyQgDla65eDYc8yCG4Pd/wK7p4Xvl2RKk5hRESqnaxcB58t35N32qZlsC+DOtWlaZ0atK3rR/2a3iXfeFoSvBMB2enQ7R64cmIpVS1SdSmMiEi19fmfe3hx3mZOZOe/w2ZMdHPu6d0Mu1sJbySM+xK+GwVuXjB0BjTpfeHFilRhCiMiUq0dy8hm1pq9bDyQyvdxB/Lm1/R253+3R5ZsSPnsDHi3K6TuM6cHvmOOSyIiBVIYERE5JSMrl/eWbGfy4h0A1PKx82B0c67rXB8fj2KOuJqaCD89aN5p4+4Nd/0Ogc3LoGqRyk9hRETkXzbsT+GmD2LzTt80DvTh3aGdaFu3mEdJstJhyiVwbJc5XbcTtB4I4f8Bv7qlXLVI5aUwIiJSgMPpWby7aDvTlu/Om/f4Va24s2fT4m0o9QB8Pxp2LgHDac7z8INLxkD4UPALLa2SRSothRERkXP4e99x7v9qHbuPnADg5zE9aRZUA5u1mLcCpx+ELT/Byo/g4EZznr0GPLgRvAJKt2iRSkZhRETkPNKzcgl/5hccp54C7O/lzlXtQ3noihYE1vAo3sZyMuGvr+CnMea03Rcuvg96PQqlMdaJSCWkMCIiUgQzVyUwe+1+NuxPIeMftwJ3ahBAZOPaXN0htHh33uxdCXPugqM7zelWV8N1U8F+AWOciFRSCiMiIsXgcBrMW5/Iq7/Es+fUqZvTWgTXILp1MA9EN8fDzXb+jTmdsPApiH3XnA4fAle+pNM2Uu0ojIiIlNDeoydYsesoS+IPMn9DUt5pnKgmtXlnaKein8JZ+xn8cJ/ZtrpBUGuo2xk6DoUG3cqoepGKQ2FERKQUHMvIZva6/Uz6eQuZOU56tajDJ7d2KfqFrn/NhGVvnbm4FcAnCB6KB2sJR4IVqSSK+vdb/yWIiJxDTR87t1/SmM9ui8RmtfDb1kO8smBL0TcQPhhGLYMx6+GGT8x5GQdhzSdlU7BIJaQwIiJSBF0b1+Kp/q0B+OD3nUycv5mT/3r2TaEsFghoAO2uh863mPPmPgQf9IKYZ81n3hxPKKPKRSo+naYRESmiHIeTmz6IZV3CcQAubRXEx7dchKU4t+7mZsOHvfOftgGweZi3AXd/AGzupVaziCvpmhERkTKQ63DyWewenv1pEwBtQv3o2aIOlzQLpEvjmkW72ybnJOz8DY7ugCPb4UAcHFhrLgtuDy2ugJqNzFdQG/AJLKvdESlTCiMiImXof7G7efanTeQ4zvwKrRfgxc1RDbm6Qyj1axZjXBHDgL9nmg/gy8l/WzE2O9y9DOq0KKXKRcqPwoiISBk7mJbJ0m2HWbb9CEviD3IkIxsAqwWm39aVqCa1cbMV49K8w9sgfj4c3wPHdsP2X835ba6BG6dpJFepdBRGRETKUVpmDt+s2cczP27Km9c40Ic3B3ckPCygZBvdMhdmDDXbfV+EbvcokEilolt7RUTKka+nOyO6NyZ2/KUM6RoGwK7DGVz73jK+XbMPp7ME/9/Xqj90HGa2f34cpvaBNdPN0zoiVYiOjIiIlIHdhzMYMzOOuL3HAajtY6dliC/NgmowqGNdIhrWKtqGcrPM239XfQS5mea81gOh93jwrw+e+p0oFZdO04iIuFhWroPnf9rMnHX7Sc/KzZvv7+XOmieji3c9Sfoh+OVJ+HtG/vke/mYoqdMC+k4Ev9BSql7kwimMiIhUEFm5DjYdSGVrchrjvl0PgI/dxqs3htMsqAahAV7U8HAr2sY2zoE/p8DhrXDyaP5ldVpDxyHgVw98Q81g4lcP3Ir4LB2RUqYwIiJSAX23bj9jZsadNb9dPT8eu7I1lzQvxpgiWemQuh9WTIHVhQwv7+EHfZ6Ai0YolEi5UxgREamgFm1J5n+xezhwPJPElJOkZp45hdO5QQBN6tQgyNeDYD9PWob40q1J7fNvdMci2LEY0hIhNRHSDsDxveDMMZfba0CNIPCqab48A860vWqCbzDU7woBYWWz01ItKYyIiFQSh9OzeHfRdj7/cw+5Bdx18+mILvRpGVT8DTudsPpjWPqGeQSlKPwbQOfh0OV28C7iRbYihVAYERGpZBKOnGDl7qMcTMvkYGoW05bvBsDNauHD4RFc2iq4ZBt25MDRnXDy2L9ex8+0j+6ExL/AOPXwP6u7eWtxv1fMoyYiJaAwIiJSyR1Oz6LHy4s5mWMGhBbBNejfvi6jejfF7lYGw0RlpZsjwC57C5LXn5lfpxWEdIDWV0PLq/QgPykyhRERkSpg37ETPPPjJhZvOZh3CqeWj52+bUMIrGGnto+d2jU8aFjbm/b1/Iv3BOFz2bEIfhxjDk3/T/5hMPh/EBIOVo2bKeemMCIiUoWknMjh8xV7eGfRNjJznAX2uap9CFe1D6Wmt918+bhT28ej5EdRDAPSkiBpPexZBsvePLPMJwgGvg0t+5Vs21ItKIyIiFRBx09kE7P5IPuOneRwehZHMrKYtz6p0P7edhs3XRTG5W2Cqeltx8fDho+HGzU83PBwsxbvSMqORbD8XdgRc2ZeUBuIuhfCh+hIiZxFYUREpJrIcTiZtmw3a/Yc49iJbI6dyOZoRg5HM7I41yNxbFYL3nYbNTzcCPC20zyoBi1DfOnTMojWob6FB5WDW+DHB2DfSjBOHaUJ7QjtroOQ9ub1JT7FGC9FqiyFERGRas7pNPht2yFmr93P+n3Hych2kJGVy4lsx3nXrevvSa+WdfhPlwY0D66Bt72AEWJPHIV1/4PfX4Ws1PzL2l4LXe+Cep012Fo1pjAiIiIFcjgNTmSboSQ9K5eMrFwOpWURn5zGuoTj/BZ/iGxH/utSAmt4ENEwgMFdws6+xTj9EMR9DgfizOtLju44s8zmAbWbmUPT+54ant4vFILbQ/2Ist9ZcSmFERERKZET2bn8se0wHy/dRXxSGiknc/ItD6xhp01df0Zc3IgezQPPfuDfvtWw/G3YsxwyDhX+QUFtzNM7nn7g4WsOXe/he2r61MsrwAwzVlup76eUPYUREREpFSkncth+KJ2P/tjJ/A35L5b1dLdyWatgboioT4f6/tTysZ+51sQwzMHUju06NUR9IqQegDWfFq8A37rmXTv+9cA7EHzqnHqdatt9oLRuaZZSpTAiIiKlLi0zh0VbDvLFigQ2HUglPSs333JPdythNb2JbhNM18a1aBpYg2B/Dzzc/nFkIysNtv9qXnOSlWZeb5KVBpmp/2inmO/pByEn49xFedWCJr2hcQ/wrm0+h8fDDzxqnGrXALsv2Ir4ZGQpNQojIiJSppxOg02JqXyzZh+/bEziQEpmoX293G34e7mbL293fD3cGNSpHgFe7ni4WfF0t516nWq72fBwt+JBDpbtC2HfKsg4fOp16Mx77smiF+zmZQYTD18zpNhrgLsnuJ16uXuZF9u6eRUw/3Tb01zu5nFmvruXORicm70UfqpVS5mGkcmTJzNp0iSSkpIIDw/nnXfeoWvXroX2nzVrFk899RS7d++mefPmvPzyy1x11VVF/jyFERGRii8r10Hi8Uz+2necnzcmseNgBrsOZ5x1MWxx2KwWWgT70jrUFz9P93zjpPjY3fCzZRGWtYNGx5fjeWiDeTQlO/3UEZdTbUd2Ke5lIaxu5rUtQa2hTmsIamVeE1OzcbU+IlNmYWTmzJkMHz6cKVOmEBkZyZtvvsmsWbOIj48nKOjsp0ouX76cnj17MnHiRK6++mq+/PJLXn75ZdauXUu7du1KdWdERKRicToN0jJzSTmZQ8rJHI6fzObPnUeI23ucE9kOMnOcZOU4yMxxkJnrNN9zHOccH6UwAd7ueLvbsLtZ8Th1ZMVus+Ljlou/NRtfSyZ+1kx8LSfxsZzEmyw8ycFOFnYjB7uRhd3Ixk427s4s3I1s3I1s3JxZp17Z2JyZuDmzsDmysDqysDkyseWkY3UUclTIZod6EdD2Omg9wLyTqBopszASGRlJly5dePfddwFwOp2EhYVx33338dhjj53Vf/DgwWRkZPDTTz/lzevWrRsdO3ZkypQppbozIiJS+RmGQY7DIDPXQcqJHDbsT2HXkQwysnJJz8wlPcscLyUjO5f0rFwOpmax/3gxTteUfsXU5QgtrPto736Atu77ac4+whwJ2I2s/D3rtMbSuId5TYvN3TyiYnM3n5Kcb9rtPMvczSMuRV3mogt8i/r3u1jHjrKzs1mzZg3jx4/Pm2e1WomOjiY2NrbAdWJjYxk7dmy+eX379uW7774rzkeLiEg1YbFYsLtZsLtZ8fN0J6yW93nXSc/K5cDxk2TlOMnKdZCd6yQr12xn5bXNozDZDidZOU6yHU5yHU5yHAa5TicOpxmCch1Ocpzme67DyNfOdTrJ/Ue/XKdBVo6DIxlBLMkNZElWRziVPyw4aWA5yGXWdQywxRJu2YH10GY4tLlsf4AFcGDDabHhsLjhsLjhPPU607bhvOYDQlpFlnttUMwwcvjwYRwOB8HB+Qe8CQ4OZsuWLQWuk5SUVGD/pKTCn6WQlZVFVtaZNJmamlpoXxERkRoebrQI9nXZ5xuGQXpWLofTszmUlsXh9CwOpZmvXUc68UjSjaQeSaaLsZ721p14kIM7ubjhwN3iwA0HbuTihtOcd2qZm8WBe94yx6nluf9Y5x99Ty2zWc4+4WHDgc1w4G4Ufv3MlvQMQsryh3QOFfKqmokTJ/LMM8+4ugwREZEisVgs+Hq64+vpTuNAnwL7OJ0GB9MGkphyEofTwGmA0zBwGgZGXtt8zzUMsp35551uG6fbp5b/e13D6QBnLjhywJmLxZmDkdc2p0/Pt/zjvWdY0a7jLAvFCiOBgYHYbDaSk5PzzU9OTiYkpOA8FRISUqz+AOPHj893aic1NZWwsLDilCoiIlKhWK0WQvw9CfH3dHUpFU6xnvdst9uJiIggJubM46OdTicxMTFERUUVuE5UVFS+/gALFy4stD+Ah4cHfn5++V4iIiJSNRX7NM3YsWO55ZZbuOiii+jatStvvvkmGRkZjBgxAoDhw4dTr149Jk6cCMADDzxAr169eO211+jfvz8zZsxg9erVfPjhh6W7JyIiIlIpFTuMDB48mEOHDvH000+TlJREx44dWbBgQd5FqgkJCVitZw64XHzxxXz55Zc8+eSTPP744zRv3pzvvvuuyGOMiIiISNWm4eBFRESkTBT173exrhkRERERKW0KIyIiIuJSCiMiIiLiUgojIiIi4lIKIyIiIuJSCiMiIiLiUgojIiIi4lIKIyIiIuJSCiMiIiLiUsUeDt4VTg8Sm5qa6uJKREREpKhO/90+32DvlSKMpKWlARAWFubiSkRERKS40tLS8Pf3L3R5pXg2jdPp5MCBA/j6+mKxWEplm6mpqYSFhbF3795q97yb6rzvoP3X/mv/q+v+V+d9B9fsv2EYpKWlUbdu3XwP0f23SnFkxGq1Ur9+/TLZtp+fX7X8RwnVe99B+6/91/5X1/2vzvsO5b//5zoicpouYBURERGXUhgRERERl6q2YcTDw4MJEybg4eHh6lLKXXXed9D+a/+1/9V1/6vzvkPF3v9KcQGriIiIVF3V9siIiIiIVAwKIyIiIuJSCiMiIiLiUgojIiIi4lLVMoxMnjyZRo0a4enpSWRkJCtXrnR1ScU2ceJEunTpgq+vL0FBQVxzzTXEx8fn69O7d28sFku+1913352vT0JCAv3798fb25ugoCAeeeQRcnNz8/VZsmQJnTt3xsPDg2bNmjFt2rSy3r3z+r//+7+z9q1Vq1Z5yzMzMxk9ejS1a9emRo0aXH/99SQnJ+fbRmXdd4BGjRqdtf8Wi4XRo0cDVeu7//333xkwYAB169bFYrHw3Xff5VtuGAZPP/00oaGheHl5ER0dzbZt2/L1OXr0KMOGDcPPz4+AgABuv/120tPT8/X5+++/6dGjB56enoSFhfHKK6+cVcusWbNo1aoVnp6etG/fnnnz5pX6/v7bufY/JyeHcePG0b59e3x8fKhbty7Dhw/nwIED+bZR0L+Xl156KV+fyrj/ALfeeutZ+3bllVfm61NVv3+gwN8DFouFSZMm5fWpFN+/Uc3MmDHDsNvtxieffGJs3LjRGDlypBEQEGAkJye7urRi6du3r/Hpp58aGzZsMOLi4oyrrrrKaNCggZGenp7Xp1evXsbIkSONxMTEvFdKSkre8tzcXKNdu3ZGdHS0sW7dOmPevHlGYGCgMX78+Lw+O3fuNLy9vY2xY8camzZtMt555x3DZrMZCxYsKNf9/bcJEyYYbdu2zbdvhw4dylt+9913G2FhYUZMTIyxevVqo1u3bsbFF1+ct7wy77thGMbBgwfz7fvChQsNwFi8eLFhGFXru583b57xxBNPGLNnzzYAY86cOfmWv/TSS4a/v7/x3XffGX/99ZcxcOBAo3HjxsbJkyfz+lx55ZVGeHi48eeffxp//PGH0axZM2PIkCF5y1NSUozg4GBj2LBhxoYNG4yvvvrK8PLyMj744IO8PsuWLTNsNpvxyiuvGJs2bTKefPJJw93d3Vi/fr3L9v/48eNGdHS0MXPmTGPLli1GbGys0bVrVyMiIiLfNho2bGg8++yz+f49/PN3RWXdf8MwjFtuucW48sor8+3b0aNH8/Wpqt+/YRj59jsxMdH45JNPDIvFYuzYsSOvT2X4/qtdGOnatasxevTovGmHw2HUrVvXmDhxogurunAHDx40AOO3337Lm9erVy/jgQceKHSdefPmGVar1UhKSsqb9/777xt+fn5GVlaWYRiG8eijjxpt27bNt97gwYONvn37lu4OFNOECROM8PDwApcdP37ccHd3N2bNmpU3b/PmzQZgxMbGGoZRufe9IA888IDRtGlTw+l0GoZRdb/7f/8ydjqdRkhIiDFp0qS8ecePHzc8PDyMr776yjAMw9i0aZMBGKtWrcrrM3/+fMNisRj79+83DMMw3nvvPaNmzZp5+24YhjFu3DijZcuWedM33XST0b9//3z1REZGGnfddVep7uO5FPTH6N9WrlxpAMaePXvy5jVs2NB44403Cl2nMu//LbfcYgwaNKjQdarb9z9o0CDj0ksvzTevMnz/1eo0TXZ2NmvWrCE6OjpvntVqJTo6mtjYWBdWduFSUlIAqFWrVr75X3zxBYGBgbRr147x48dz4sSJvGWxsbG0b9+e4ODgvHl9+/YlNTWVjRs35vX558/rdJ+K8PPatm0bdevWpUmTJgwbNoyEhAQA1qxZQ05OTr66W7VqRYMGDfLqruz7/k/Z2dl8/vnn3HbbbfkeJFmVv/vTdu3aRVJSUr46/f39iYyMzPddBwQEcNFFF+X1iY6Oxmq1smLFirw+PXv2xG635/Xp27cv8fHxHDt2LK9PRf95gPm7wGKxEBAQkG/+Sy+9RO3atenUqROTJk3Kd0qusu//kiVLCAoKomXLlowaNYojR47kLatO339ycjJz587l9ttvP2tZRf/+K8WD8krL4cOHcTgc+X4BAwQHB7NlyxYXVXXhnE4nY8aMoXv37rRr1y5v/tChQ2nYsCF169bl77//Zty4ccTHxzN79mwAkpKSCvxZnF52rj6pqamcPHkSLy+vsty1QkVGRjJt2jRatmxJYmIizzzzDD169GDDhg0kJSVht9vP+mUcHBx83v06vexcfVy97//23Xffcfz4cW699da8eVX5u/+n07UWVOc/9yMoKCjfcjc3N2rVqpWvT+PGjc/axullNWvWLPTncXobFUFmZibjxo1jyJAh+R6Edv/999O5c2dq1arF8uXLGT9+PImJibz++utA5d7/K6+8kuuuu47GjRuzY8cOHn/8cfr160dsbCw2m61aff/Tp0/H19eX6667Lt/8yvD9V6swUlWNHj2aDRs2sHTp0nzz77zzzrx2+/btCQ0N5bLLLmPHjh00bdq0vMssVf369ctrd+jQgcjISBo2bMjXX39dIf5IlqePP/6Yfv36Ubdu3bx5Vfm7l4Ll5ORw0003YRgG77//fr5lY8eOzWt36NABu93OXXfdxcSJEyvk0ODF8Z///Cev3b59ezp06EDTpk1ZsmQJl112mQsrK3+ffPIJw4YNw9PTM9/8yvD9V6vTNIGBgdhstrPuqkhOTiYkJMRFVV2Ye++9l59++onFixdTv379c/aNjIwEYPv27QCEhIQU+LM4vexcffz8/CrUH/2AgABatGjB9u3bCQkJITs7m+PHj+fr88/vuars+549e/j111+54447ztmvqn73p2s913/TISEhHDx4MN/y3Nxcjh49Wir/HirC747TQWTPnj0sXLjwvI+Hj4yMJDc3l927dwOVf///qUmTJgQGBub7t17Vv3+AP/74g/j4+PP+LoCK+f1XqzBit9uJiIggJiYmb57T6SQmJoaoqCgXVlZ8hmFw7733MmfOHBYtWnTWIbaCxMXFARAaGgpAVFQU69evz/cf6ulfZG3atMnr88+f1+k+Fe3nlZ6ezo4dOwgNDSUiIgJ3d/d8dcfHx5OQkJBXd1XZ908//ZSgoCD69+9/zn5V9btv3LgxISEh+epMTU1lxYoV+b7r48ePs2bNmrw+ixYtwul05oW0qKgofv/9d3JycvL6LFy4kJYtW1KzZs28PhXx53E6iGzbto1ff/2V2rVrn3eduLg4rFZr3umLyrz//7Zv3z6OHDmS7996Vf7+T/v444+JiIggPDz8vH0r5PdfKpfBViIzZswwPDw8jGnTphmbNm0y7rzzTiMgICDfXQWVwahRowx/f39jyZIl+W7XOnHihGEYhrF9+3bj2WefNVavXm3s2rXL+P77740mTZoYPXv2zNvG6ds7r7jiCiMuLs5YsGCBUadOnQJv73zkkUeMzZs3G5MnT64Qt7c+9NBDxpIlS4xdu3YZy5YtM6Kjo43AwEDj4MGDhmGYt/Y2aNDAWLRokbF69WojKirKiIqKylu/Mu/7aQ6Hw2jQoIExbty4fPOr2neflpZmrFu3zli3bp0BGK+//rqxbt26vLtFXnrpJSMgIMD4/vvvjb///tsYNGhQgbf2durUyVixYoWxdOlSo3nz5vlu7Tx+/LgRHBxs3HzzzcaGDRuMGTNmGN7e3mfd2ujm5ma8+uqrxubNm40JEyaUy62d59r/7OxsY+DAgUb9+vWNuLi4fL8LTt8ZsXz5cuONN94w4uLijB07dhiff/65UadOHWP48OGVfv/T0tKMhx9+2IiNjTV27dpl/Prrr0bnzp2N5s2bG5mZmXnbqKrf/2kpKSmGt7e38f7775+1fmX5/qtdGDEMw3jnnXeMBg0aGHa73ejatavx559/urqkYgMKfH366aeGYRhGQkKC0bNnT6NWrVqGh4eH0axZM+ORRx7JN9aEYRjG7t27jX79+hleXl5GYGCg8dBDDxk5OTn5+ixevNjo2LGjYbfbjSZNmuR9hisNHjzYCA0NNex2u1GvXj1j8ODBxvbt2/OWnzx50rjnnnuMmjVrGt7e3sa1115rJCYm5ttGZd33037++WcDMOLj4/PNr2rf/eLFiwv8t37LLbcYhmHe3vvUU08ZwcHBhoeHh3HZZZed9TM5cuSIMWTIEKNGjRqGn5+fMWLECCMtLS1fn7/++su45JJLDA8PD6NevXrGSy+9dFYtX3/9tdGiRQvDbrcbbdu2NebOnVtm+33aufZ/165dhf4uOD3mzJo1a4zIyEjD39/f8PT0NFq3bm28+OKL+f5YV9b9P3HihHHFFVcYderUMdzd3Y2GDRsaI0eOPOt/Lqvq93/aBx98YHh5eRnHjx8/a/3K8v1bDMMwSucYi4iIiEjxVatrRkRERKTiURgRERERl1IYEREREZdSGBERERGXUhgRERERl1IYEREREZdSGBERERGXUhgRERERl1IYEREREZdSGBERERGXUhgRERERl1IYEREREZf6f68/NyBxR76UAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "\n", "plt.plot(timeline, sf_values[0])\n", "plt.plot(timeline, sf_values[1])" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.11.0rc1" } }, "nbformat": 4, "nbformat_minor": 5 }