{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Lifelines examples\n", "##### Gehan Survival Dataset\n", "\n", "The data show the length of remission in weeks for two groups of leukemia patients, treated and control, and were analyzed by Cox in his original proportional hazards paper. The data are available in a file containing three columns:\n", "\n", "- Treatment: coded Treated (drug) or Control (placebo),\n", "- Time: weeks of remission,\n", "- Failure: coded 1 if a failure (relapse), 0 if censored\n", "\n", "Thus, the third and fourth observations, 6 and 6+, corresponding to a death and a censored observation at six weeks, are coded 6, 1 and 6, 0, respectively." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "WARNING: pylab import has clobbered these variables: ['datetime']\n", "`%matplotlib` prevents importing * from pylab and numpy\n" ] } ], "source": [ "from lifelines.estimation import KaplanMeierFitter, NelsonAalenFitter\n", "import pandas as pd\n", "import numpy as np\n", "%pylab inline\n", "figsize(12.5,6)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/camerondavidson-pilon/.virtualenvs/data/lib/python2.7/site-packages/IPython/kernel/__main__.py:1: ParserWarning: Falling back to the 'python' engine because the 'c' engine does not support regex separators; you can avoid this warning by specifying engine='python'.\n", " if __name__ == '__main__':\n" ] } ], "source": [ "data = pd.read_csv('../lifelines/datasets/gehan.dat', sep=\"\\s{1,3}\", header=None)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "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", "
012
371121
381151
391171
401221
411231
\n", "
" ], "text/latex": [ "\\begin{tabular}{lrrr}\n", "\\toprule\n", "{} & 0 & 1 & 2 \\\\\n", "\\midrule\n", "37 & 1 & 12 & 1 \\\\\n", "38 & 1 & 15 & 1 \\\\\n", "39 & 1 & 17 & 1 \\\\\n", "40 & 1 & 22 & 1 \\\\\n", "41 & 1 & 23 & 1 \\\\\n", "\\bottomrule\n", "\\end{tabular}\n" ], "text/plain": [ " 0 1 2\n", "37 1 12 1\n", "38 1 15 1\n", "39 1 17 1\n", "40 1 22 1\n", "41 1 23 1" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data.tail()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " 0 1 2\n", "0 2 6 1\n", "1 2 6 1\n", "2 2 6 1\n", "3 2 6 0\n", "4 2 7 1\n" ] } ], "source": [ "data = data.dropna()\n", "print data.head()" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false }, "outputs": [], "source": [ "treatment = data[0] == 2" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": false }, "outputs": [], "source": [ "T = data[1]\n", "E = data[2]\n", "t = np.linspace(0,40,150)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 20, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(0, 1.05)" ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuAAAAF9CAYAAABFzt0JAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X24XWV97vv75i1VMCLiMYJA3KDY7K0sbMHUooZMqohu\n6e6pbegRxLqV3XMo51RAG7WAuOtLwcrmpFpq0Qj2ktaqDW7BKIsEwRoMLRPqBpQASYAYxBJISgSl\n+e0/5lgwWaz3NX5jjDnW93Nd60rGnGM+41k3I4NnPvM3n+GIEAAAAIBq7FZ3BwAAAIC5hAE4AAAA\nUCEG4AAAAECFGIADAAAAFWIADgAAAFSIATgAAABQoT2qOtAnP/nJGBoaqupwc0q32xXZ5iDbXOSb\nh2zzkG0ess1Dtnm63a7OPPNMT+c1rmod8FNOOSUuu+yyaXUOU2P7vIg4r+5+tBHZ5iLfPGSbh2zz\nkG0ess0zkzFuZSUoW7durepQc9HCujvQYgvr7kDLLay7Ay22sO4OtNjCujvQYgvr7kCLLay7A3gK\nNeAAAABAhSobgL/xjW+s6lBz0cq6O9BiK+vuQMutrLsDLbay7g602Mq6O9BiK+vuQIutrLsDbXXE\nEUdM+zWVDcAp/M8TEWvr7kNbkW0u8s1DtnnINg/Z5iHbPDMZ41a2Ckq321Wn06nqcHOK7SX8w8pB\ntrnINw/Z5iHbPDPNdnh4eDdJfyRpvqRqVpcYMD/60Y8WvuxlL9tYdz8GlCVtl/SpTqezq4wGKxuA\nAwAAJPkjSdd0Op1b6u5IUx133HG8cZyF4eHhI9Q7zz5ZRnuUoLQA/6DykG0u8s1DtnnINs8ssp3P\n4HtinLezU5xf88tqj1VQAADAoKPsBFUo7TyrbADe7XarOtScY3tJ3X1oK7LNRb55yDYP2eYh2zxk\n2yyTDsBtf872A7b/ZYJ9LrZ9p+1bbB9ZbhcBAADaz/Zrbd8xwfMLbe+yTQXDgJvKf8DPSzp+vCdt\nnyDpsIh4qaT3SPrMWPtRA56Huq48ZJuLfPOQbR6yzdO2bG0vt33VqMfuHOex34mI6yPi5X2Pb7S9\ndBbHX2n7I1JOtraX2L637Hancfy1tt9V1/FnY9IBeERcL2nbBLu8VdIXin1vlLSv7ReW0z0AAICB\ndZ2k19i2JNl+kXor0A2NzGIXjx0q6TtjvD7UWwIvhe1BXw1vYGv/y/gI40BJ/e9+7pP04tE7UQOe\nh7quPGSbi3zzkG0ess3TwmxvkrSnpJEygNdKWiPpR6Me2xARW/tnlG1fLulgSV+3vcP2WX3tvt32\nJtsP2v7AWAe2/R5JvyfpfcXrbyge32j7fbZvlbTD9m62F9v+R9vbbHdtv76vnXfavs32dtt3Fe3K\n9t6SrpZ0QNH+dtsvsn2e7S/bvrx47FbbLy0+DXig6Pdv9LX/XNuX2t5i+z7bH+l7c3Kq7RtsX2D7\nIdt32z6+eO5Pi+xWFMe/eKb/kepQ1juf0e/Oxn5HYpfyTuWGl79af/T7Hy2jqVZ42Wmf1FEXXFt2\ns1etP3vpm8tuFACAuSIifm77Rkmvl3SzpNdJul7SluLv/1z8+YzZ74g42fYxkt4VEddKvRrw4ulf\nl/QySYdL+r7tr0bEHaNe/1e2f03SvRFxzqg3N8skvUnSTyW9SNL/lPT2iPim7eMkfcX24RHxr5Ie\nkPTmiLjH9uskXW17fUTcXAyGvxgRB400XEz2v0W9ColTJX1O0rclXSLpAEnvLP7+H4qXrJS0Vb1P\nAfYp+nKvpL8qnj9avXLo50s6TdKlkg6MiA/afo2kyyPic2P/F2iuMgbg90s6qG/7xcVjT7Nhwwad\nKmlhsb2vem/9lhTba4s/p7J9zB03avtdvRn1+Yf23kDO5e35hw5ltH+C3Xly0f6Rf7hss8324GyP\naEp/2rI98lhT+tOm7YhYO5PXf/rTn144crft0c8fdcG1pZUprD97qWdwvmyS9FuSLpJ0jKR/lPSI\negPvi9T7nt3f9h1mXv/5JekI27v6zz9J34qIxyXdavtuSSdJOnec4y8c1d684vX3F/ufI+mfIuKb\nxfNPSLpb0pslXSZpp6RDJN0TEd+x/c/qDaJvVm8CdnR/F0r6l4j4dtH+HZLeJunjERG27y/6NF/S\ns4rjvCUifibpZ7ZXS/pvemoA/hNJdxWvvUzSp23/ZkT8Q/H8y2v49zik3lBWkhZeeOGF077buyMm\nPy/de8f19Yh4xRjPnSDp9Ig4wfZiSRdFxOLR+w0PD0en05l9HdPILHpEWk3UXDdysVp/9lIyBgA0\n3vDw8HmdTue8sZ4rewA+3dfYPla9Afbhkn4QEQcWg88fSfplSQ9KOjQiNhWDu8ujmFG2fY+eOQN+\nt6Q9ImJX8dgajTMLbPvzku6LiD/pe+weSf81IoaL7U+rN6B+rO+le0j6SET8me03qTe4f6l6pcvP\nVm8wfe7o/hbtnVf8PicX28dJ+mxEvKTY3kPSz9WbsH2xpO+pd5v3EbtJ2hwRr7B9avH7v7av/V3q\nLf5x90S/e4bxzrOZjHEnnQG3/SX1PjrZ3726pHPVq2dSRFwSEVfZPsH2BkmPqvcf8Rm63e603x1g\naka980SJyDYX+eYh2zxkmycj2wZMJq2T9FxJ75b0XUmKiO22t6i3etyWiNg0zmtn++bhydePyra/\n3c3qDWLfM/rFtudJ+oqkt0taFRH/bvtr0pOlx2P1bzp9vlfS45KeP/KGYpoG9kuYkw7AI+KkKexz\nejndAQAAaI+I+JntmyS9V9J/73vqhuKxb03w8gfUq42e7Ite473JeEBP1VqP54uS1tt+g6Rh9SZZ\nF0u6U72Z6b3UqxXfVcyGv0HSyL1hHpD0fNvzI2JkFnvKb3gi4se2vyXpz23/iXoTuS9Rr8Z7rFVh\nxvr9Dp3q8ZqksoXcWQc8DzMxecg2F/nmIds8ZJunxdleJ+kF6g26R1wvaX898wuY/bO6H5P0IfdW\nJ3nvGM+P9Zp+l0paZHubpDPG2iEi7pN0oqQPqFdvvVnSmeqVKe8oXvd3kh5Sr9Z8Vd9r75D0JUl3\nu7dKyYuKvozuz0Tbp6g3yL+tOMaXJS3o22+i1/4PSb9dHPuisX6/phr09R8BAAAaLSI+oN4At/+x\nL6s32Ox/bK16Sw+ObF8p6cpRze0+6jXHTnDcDZKOHPXYS8bY7/t6at2L0c99WtKnJzjGuyS9q++h\nD496/hr1zcJHxBP9v0Mxc/5/Fz+j2/6CinvN9D3W/9p16tXWD5zKZsBZBzzP6BUPUB6yzUW+ecg2\nD9nmIds8ZNsslQ3AAQAAAFRYglJ6Dbi9rIRWtilidQnt1KrFNXO1I9tc5JuHbPOQbR6yzUO2zTLI\nNeBbS2hjweS7AAAAAOWhBrwFqOvKQ7a5yDcP2eYh2zxkm4dsm4UacAAAAKBCrAPeAtR15SHbXOSb\nh2zzkG0ess1Dts3CDDgAAABQIWrAW4C6rjxkm4t885BtHrLNM9eztf1a23dM8PxC27tsT3v8Ntez\nbZpBXgUFyY664NpSlnpcf/bSgV/qEQCA6bK9XNJrI+KEvsfulHTnGI99MCL+TtLL+x7fKOn3I+La\nCvq6UtK9EfEnE+yzS9JhEXF3XX3IYnuhpLsl7RERu7KPN7jrgONJiXVdc36pR2rmcpFvHrLNQ7Z5\nWpjtdZLeb9sREbZfpN7Ya8j2bhGxq3jsUEnfGeP1IclldKTEbMftj+09ilvND7JS8p4MNeAAAAA5\nbpK0p6SRWcjXSloj6UejHtsQEVttL7F9ryTZvlzSwZK+bnuH7bP62n277U22H7T9gZEHbc+zfZHt\n+4ufT9neq3juVNvX93euKGc51PZ7JP2epPcVx1o1+hexPfIG4ZZin7cV/b3P9vts/1jSpe75Y9sb\nbP/U9t/afl5fO1+2/WPbD9u+zvai4vEx+2B7o+2zbN9aPH6p7Rfavtr2I7a/bXvfvvYX2/5H29ts\nd22/vu+5tbbPt32D7e22V9t+fvH0yO/3cHGcV0/w33XWqAFvAeq68pBtLvLNQ7Z5yDZP27KNiJ9L\nulHSyCDwdZKul3RD8feRx54x+x0RJ0vaLOktEfGciLiw7+lfl/QySR1J59g+vHj8g5KOlnRE8XO0\npA8Vz71cY4uI+CtJfyPpE8WxThxjp5H+vrLY58vF9gslPU+9NwunSTpD0luL3+tFkrZJ+ou+pr4h\n6TBJL5D0z8VxNUEfQtJvFb/r4ZLeIulqSX8s6f9Qbyx7hiTZPlDS/5R0fkQ8T9JZkr7SN8iWpJMk\nnVq8dq9iH6n3RkiSnlsc/8Zx8ioFM+AAAKC97CjtZ2au01OD7WPUG2xf3/fYa4t9puPDEfF4RNwq\n6Rb1BttSbwb5/Ij4aUT8VNKHJZ08jXZnUn6xS9K5EfGLiHhMvUH4hyJiS0T8oujDb498cTQiVkbE\no33PHWH7OZP04f+PiAcjYot62X0vIm6JiMclfU3SkcV+b5d0VUR8szjWNep9CvHm4vmQ9PmI2FD0\n9e/01CcRlZSejGAd8BZoYc1cY5BtLvLNQ7Z5yDZPS7P9jqRjijKMF0TEXZK+J+k1xWP/UWPXf0+k\n/ztaOyXtU/z9AEmb+p7bXDwmSeOurjJLDxYz/SMWSvpaUQKyTdJtkp6Q9ELbu9v+eFGe8oike4rX\n7D/JMR7o+/vPRm0/pqd+/0MkvW3k2MXxf11P/z5af3Y/63ttpVgFBQAAtFdEpTObY1gn6bmS3i3p\nu5IUEdttb5H0HklbImLTOK+d7qz7FvUGwLcX2wcXj0nSo5KePbKj7dGLJMx0hn/06zZLemdEfG/0\njrZPVq88pRMRm4ra7Yf01OzzVPsw3n/TzZIuj4j3TLGdfjP9/WeEGvAWaFvNXJOQbS7yzUO2ecg2\nTxuzjYifqVcG8V49fab7huKxicpPHlBvhZSp+pKkD9ne3/b+ks6RdHnx3DxJ/9H2EbZ/SdJ5Yxzr\nP0zS/lT685eSPmr7YEmy/QLbby2e20fS45Iesr23pI/OoA8T+aKk/2z7DcVs+y8VXxQ9sG+f8Qbv\nD6pXTjOdvGeMGnAAAIBc16n3pcMb+h67Xr3Si9HlJ/0zsR9Tb0C9zfZ7x3h+tP+u3mD/1uLnpuIx\nSbpP0vmSrpH0w+L4/W1dKmlRcayvjtP+eZK+UOzz28XrR/fnf0i6UtK3bG9Xr9zm6OK5y9Qrkblf\n0g+K56bbB416zZN9iIj7JJ0o6QOSfqLejPiZevqge7zX7pT0p5K+Wxz/aCVyRDUz7sPDw9HpdGb/\nMdBTX4I4dtZt9WqFHi+hHUnapohW3HDmqAuuLTPjBevPXnpFCe0AADCm4eHh8zqdznl19wPtNt55\nNpMx7lyvAR+v5momBvqGM4nmlXRHzabgzp4AAGBWKhuAd7tddTqdqg43p9he0uBvjpf5Jqdym1et\nGDr4xNP7v8DAG60SNfzcHWhkm4ds85BtHrJtFmrAAQAAgAqxDngL8I42z6jZb5SMczcP2eYh2zxk\nm4dsm2WQa8DXlNDGOknLS2gHAADUp+61vjE3lHaeDeI64FeV1ZCkxSW2VZs2rpvaFJtXreCjm0Sc\nu3nINg/Z5plFtnsMDw8zCJ8A5+3sFOdXaRPXgzcDHvFm2cv09FuJzkQZM+gAAKB+V0l6h6SVNfcD\n7fUOSd8oq7HKBuDUgOehrisPNeC5OHfzkG0ess0z02w7nc53h4eHDx8eHj5fvbsZYpRrrrlGw8PD\nS+rux4DaTdLGTqfzj2U1OHgz4AAAAKN0Op3P1d0HYKpYB7w884rSmOma9R00WdszzxjrgKNEnLt5\nyDYP2eYh2zxk2yzMgJdnpjec4cYug6VNd/bkrp4AANSAGvAW4B1tnjFmvwf6zp6j1P7mj3M3D9nm\nIds8ZJuHbJuFO2ECAAAAFRrEdcAxCmt75mEd8Fycu3nINg/Z5iHbPGTbLMyAAwAAABWqbABODXge\n6rrysAJKLs7dPGSbh2zzkG0esm0WZsABAACAClED3gLUdeWhBjwX524ess1DtnnINg/ZNgsz4AAA\nAECFWAe8BRLrutaU0MY6SctLaKcW1IDnoiYxD9nmIds8ZJuHbJuFGXCM5aoS21pcYlsAAAADjxrw\nFii7rmv92UvfLOkkScfO8mfgUQOei5rEPGSbh2zzkG0esm0WZsABAACAClED3gLUdeVpeQ34vKMu\nuHZZnR341T8bVol92Lb+7KWrS2pr4HFdyEO2ecg2D9k2S2UDcACNs6nuDpRsQd0dAABgKqgBbwHq\nuvJQA56LfPNwXchDtnnINg/ZNgs14AAAAECFKhuAUwOeh7quPC2vAa8d+ebhupCHbPOQbR6ybRZm\nwAEAAIAKUQPeAtR15aFGORf55uG6kIds85BtHrJtFmbAAQAAgAqxDnj95sme1TrIIaloY5siWAe5\nRNQo5yLfPNR75iHbPGSbh2ybhXXA61fmWsysgwwAANBw1IC3wBkSHy8koUY5F/nmod4zD9nmIds8\nZNsskw7AbR9v+w7bd9p+/xjP72/7m7a7tn9g+9SUngIAAAAtMGEJiu3dJa2QdJyk+yWtt31lRNze\nt9vpkm6OiOW295f0Q9tfjIgn+ttqcA34mpLaWSdpeUltTcvFEh8vJKFGORf55qHeMw/Z5iHbPGTb\nLJPVgB8taUNEbJQk21dIOlFS/wD8x5JeWfx9vqR/HT34bqh1khaX2F6ZbbVNWW9ymqC2N1oAAKAd\nJhuAHyjp3r7t+yS9etQ+n5V0re0tkp4j6XfGaqjb7arT6cy0nxnKHETVOsA8Qxpq6Cx42W9yKrf9\nrq7mH/q0T28G+vdpms2rVgyVOAs+76gLrp3VikJtsuWayxYdcNwpt9Xdj8K29Wcvbc0KTbaXMJuY\ng2zzkG2zTDYAjym08QFJ3YhYYvtQSd+2fURE7Jh99zDgBn6m+OEf3DA0/9ChkQFim2by26jMFYUG\n3hOPbl8gaWvd/SiwQhMA9JlsAH6/pIP6tg9Sbxa832sk/akkRcRdtu+RdLikm/p32rBhg2yvlLSx\neOhh9Qbua6Wnvp071e2RlT9GZn7r3l4r6at9M9FVHv9iqXuGNLS/tN85RbjTzXP09pZrLlv0xKPb\nF4zMTo6sVjGXt/f9T8c8ORvehP6wzfZk2yPq7s+Way5b5Pd1lsz0etS07ZHHmtKfNm1HxNom9Ydt\ntsfZHpK0r3oWXnjhhdOu8nDE+JPctveQ9ENJHUlbJH1f0kn9X8K0/eeSHomID9t+oaR/kvTKiHio\nv63h4eHodDqeVu/G79gyNWdmR3pqZvTYWnshLVDEFWU0VHyU36SMm6Ap/52BQbNg/dlLS7k2AUDT\nzGSMO+EyhMWXKU+XtFrSbZL+NiJut32a7dOK3T4q6Vdt3yLpGknvGz34llgHPBPrgOdhnepc5JuH\nbPOwnnIess1Dts0y6Z0wI+JqSVePeuySvr//VNJ/Lr9rAAAAQPtUdifMBq8DPvAaugJKK7BOdS7y\nzUO2eVhJIg/Z5iHbZqlsAA4AAACgwgE4NeB5qAHPQx1tLvLNQ7Z5qKXNQ7Z5yLZZmAEHAAAAKkQN\neAtQA56HOtpc5JuHbPNQS5uHbPOQbbMwAw4AAABUiBrwFqAGPA91tLnINw/Z5qGWNg/Z5iHbZpl0\nHXBMy5rJd3mGdZKWl90RAAAANBM14OVYN4vXLp7twakBz0MdbS7yzUO2eailzUO2eci2WZgBL8dM\nZ7BnMmMOAACAAVbZALzb7arT6VR1uDnlDGmIWfAcm1etGBpjJrEtb5xqL38aJ1+UgGzz2F7CbGIO\nss1Dts3CKijA1M2m1KiJZl3+BAAAps8RUcmBhoeHo9PpuJTG7GWStpbSVr1GZlKPLam9QyQ9XkZD\nXzv6hF977s7tj5TR1kw88qz5Oz76tjPX13X8OaDscw+YyIL1Zy+9ou5OAECGmYxxqQFvl01lNfTA\nvi945Od7znuorPam6wWPPLhfXccGULp5R11w7bKS2tq2/uylq0tqCwBqwTrgLcA64HlYSzkX+eZp\nWLab1PvUsoyf51Xc92dgPeU8ZJuHbJuFGnAAAACgQqwD3gKsgJKHVSRykW8ess3DShJ5yDYP2TYL\nM+AAAABAhagBbwFqwPM0rI62dcg3D9nmoZY2D9nmIdtmYQYcAAAAqBA14C1ADXge6mhzkW8ess1D\nLW0ess1Dts3CDDgAAABQIWrAW4Aa8DzU0eYi3zxkm4da2jxkm4dsm4UZcAAAAKBC1IC3ADXgeaij\nzUW+ecg2D7W0ecg2D9k2CzPgAAAAQIWoAW8BasDzUEebi3zzkG0eamnzkG0esm0WZsABAACAClED\n3gLUgOehjjYX+eYh2zzU0uYh2zxk2yzMgAMAAAAV2qOqA3W7XXU6naoON6ecIQ0xC55j86oVQy2f\nSVxT58G339XV/ENL+3RsnaTlZTU26ObAuVsb20uYTcxBtnnItlmYAQfmpnV1dyDB4ro7AADAVFQ2\nA04NeB5mv/O0eAaxETPFJc5+1zqT30QtPndrxyxiHrLNQ7bNwgw4AAAAUCFqwFuAGvA81NHmIt88\nLc523lEXXLuszg5sueayRQccd8ptdfahrcg2D9nm+firpv+aygbgAACUYFPdHXji0e0LJG2tux9t\nRLZ5yLZZWAe8BZj9ztPSGcTGIN88ZJuHbPOQbR6ybZZBnQHfJmlBCe3MUwNmU1TeF8hYhg0AAKDh\nKpsB73ZLfOMVsVoRV8z6R3q8vE7NSClLwa196q8sw1ayzatW8NFNIvLNQ7Z5yDYP2eYh22YZ1Bnw\ntihltvqr0tAS6VNltAUAAIBc1IC3ADXgeaiZy0W+ecg2D9nmIds8ZNssrAMOAAAAVGgwa8DxNGdI\nfLyQhJq5XOSbh2zzkG0ess1Dts3CDDgAAABQIWrAW4Aa8DzUzOUi3zxkm4ds85BtHrJtFmbAAQAA\ngApRA94C1IDnoWYuF/nmIds8ZJuHbPOQbbMwAw4AAABUiBrwFqAGPA81c7nINw/Z5iHbPGSbh2yb\nhRlwAAAAoEKV3Yq+2+2q0+lUdbg55Qxp6OK6O9FSm1etGGLWIE9CvmtKbGugbb+rq/mHNuaTx3WS\nltfdibJwXchDtnnItlmYAQfQBuvq7gAmtLjuDgBAk1Q2A04NeB5qwPMwW5CrxHxbM7talgbNfrfu\nUwmuC3nINg/ZNktlA3AMlr0f27lj297z95ttO/++2x577r/joQfK6BMAAEAbUAPeAhk14Cd/58vr\ny2jnL9/wjqVltFMXauZykW8ess1DtnnINg/ZNsukNeC2j7d9h+07bb9/nH2W2L7Z9g9sry29lwAA\nAEBLTDgDbnt3SSskHSfpfknrbV8ZEbf37bOvpL+Q9MaIuM/2/mO1RQ14HmrA8zBbkIt885BtHrLN\nQ7Z5yLZZJpsBP1rShojYGBG/kHSFpBNH7fN7kr4SEfdJUkT8tPxuAgAAAO0w2QD8QEn39m3fVzzW\n76WS9rO9xvZNtk8eq6FulzdeWc6Q+HghyeZVK8g2EfnmIds8ZJuHbPOQbbNM9iXMmEIbe0p6laSO\npGdL+p7tdRFxZ/9O1113nc4666yVkjYWDz0sqRsRa6VeHbkkVbn9YWnROdJW6alB7Eg5x6Btry1C\nXVL8WXd/RrYXFf25esuGwyTpTQcctmEq2994aOsh/V8YGblwVL09oq7jt32bfPO2d/b+TTWiP9vv\n6urhH9xQ+7/nsrZ3btlw2OZVKxrTH7bZnsr2iKb0Z5C3d27ZcNiunz+2jyQ9sXPHgu4Hz5r2QiOO\nGH+MbXuxpPMi4vhie7mkXRHxib593i/pWRFxXrH915K+GRF/39/W8PBwdDodT6t32exlKgbgLTCy\n1u6xtfZilL98wzuWPu/R7Q9N93UveOTB/d5/yoevzegTgEo18toEAGX5+KtizXTHuJOVoNwk6aW2\nF9reS9LvSrpy1D6rJB1je3fbz5b0akm3TacTAAAAwFwx4QA8Ip6QdLqk1eoNqv82Im63fZrt04p9\n7pD0TUm3SrpR0mcj4hkDcGrA81ADnmf0R3coF/nmIds8ZJuHbPOQbbNMeiOeiLha0tWjHrtk1PaF\nki4st2sAAABA+0x6I56ysA54HtYBzzPypQvkIN88ZJuHbPOQbR6ybZbKBuAAAAAAKhyAUwOehxrw\nPNTM5SLfPGSbh2zzkG0esm2WSWvAgTrs9cQv9vzEZecurbsf33ho6yFvvuzc/WbbziPPmr/jo287\nc30ZfQIAAINtwnXAy8Q64OnWTL7L4Nj+rH1uu/pVv/HZuvtRFtY1xxzGOuAAWm0m64DP9RnwbZIW\nlNTWPEmbSmprJtZJWlzj8Us1/2f/tmjyvQAAAAZPZQPwbrc77dt0potYXVpbvdn0WpwhDV0sLa/r\n+BOZyZ0wf/e7X/tUVn+m6+otGw570wGHbai7H221edWKIb6Zn4Ns85BtHrLNQ7bNMtdnwJFs78d2\n7ti29/wZ1VDP9HVleuJZ+zynvx//vtsee+6/46EH6uwTAAAYbHO7BrxM7aonr1tja0ZnMqMvUQOO\nOa2x/54BoAwzqQFnHXAAAACgQqwD3gKsA56HbHOxLm0ess1DtnnINg/ZNgsz4AAAAECFKhuADw3x\nxivLxRIfLyQh21x8Iz8P2eYh2zxkm4dsm4UZcAAAAKBC1IC3AHXKecg2FzWJeRqY7ZqSfj5WdcdH\na2C2rUG2eci2WZgBBwBkWldye6254y+AuauyG/FQA56HOuU8ZJuLmsQ8Dcq2zLv0rpl8l3wNyrZ1\nyDYP2TYLM+AAAABAhagBbwHqlPOQbS5qEvOQbR6yzUO2eci2WZgBBwAAACrEOuAtQJ1yHrLNRU1i\nHrLNQ7Z5yDYP2TZLZV/CnAO2SVowg9fNk7Sp5L4AAACgoagBL0vEakVcMe0f6fHZHpo65Txkm4ua\nxDxkm4ds85BtHrJtFmrAAQAAgApRA94C1CnnIdtc1CTmIds8ZJuHbPOQbbMwAw4AAABUiBrwFqBO\nOc/obPd+bOeObXvP32+6P488e/5z1fuSbpN+DkmIbFqoScxDtnnINg/Z5iHbZmEVFGAaTv7Ol9fP\n8KUL/susQAU+AAAYkUlEQVSN37ii1M7M0lEXXLus7j4AADAXOSIqOdDw8HB0Oh1XcrBBYi+TtLXu\nbjTMmuLPY2vtRbkWFKveNEYxAOfcwyBp47UBwID7+KtizXTHuNSAAwAAABWiBrwFqAHPQ7a5qEnM\nQ7Z5yDYP2eYh22ZhBhwAAACoEDXgdaMGfCxtrPOkBhyYvTZeGwAMuJnUgLMKClCNecWbrTJsU8Tq\nktoCAAAVowa8BahTzlNitpvUm20u4+d5JfWpdtQk5iHbPGSbh2zzkG2zUAMOAAAAVIga8LpRAz4W\n6jwnVko9OTXgGEBcGwA0DuuAAwAAAA1HDXgLUAOeh2xzUZOYh2zzkG0ess1Dts3CDDgAAABQIWrA\n60YN+Fio85wYNeCYq7g2AGgcasABAACAhqMGvAWoU85DtrmoScxDtnnINg/Z5iHbZuFOmACAQbNm\n8l3y7Pufjqnz8K1GtnnItlkqmwEfGuKNV5aLJT5eSEK2uQ4+8XTyTdLSbNfV3QFJmn8o/z/LQrZ5\nyLZZmAEHAAyK5XV3AADGMO1P5agBbwHqlPOQbS5qEvOQbR6yzUO2eci2WVgFBQAAAKgQ64DXjXXA\nx8JavxNjHXAAABpiJuuAUwMOzF3bJC2ouxMlmidpU92dAABgMpUNwLvdrjqdTlWHm1POkIZYrSNH\nm7Ndf/bS1XX3wfaSiFhbRlvFjD4Km1etGGrpSii1I9s8ZJuHbJuFGfD6zXoWcn9pv6INZgABAAAa\njhrwNmlPPTk14BMrpQa8bahpBwDUYSY14KyCAgAAAFRo0gG47eNt32H7Ttvvn2C/o2w/Yfu3xnqe\ndcDz2F5Sdx/ainXAc3Hu5mHN3zxkm4ds85Bts0w4ALe9u6QVko6XtEjSSbZ/eZz9PiHpm5IoMwEA\nAADGMdkM+NGSNkTExoj4haQrJJ04xn5/KOnvJT04XkNDQ7zxylLWKhJ4praugNIUnLt5WO0gD9nm\nIds8ZNsskw3AD5R0b9/2fcVjT7J9oHqD8s8UD1XzrU4AAABgAE02AJ/KYPoiSX8cveVUrHFKUKgB\nz0MdbR5qwHNx7uah3jMP2eYh2zxk2yyTrQN+v6SD+rYPUm8WvN+vSLrCtiTtL+lNtn8REVf273Td\nddfprLPOWilpY/HQw5K6Ix9Bj/yPmO2Zb39YWnROsQzbyMBxpIRi0LbXSvpq301w6urPiLrzGL3d\nhPOtjO0RZbR3wBveseiA407ZKj31P5qRj1zn4vbOLRsOU3G+NKE/bdreuWXDYZtXrWhMf9hmeyrb\nI5rSn0He3rllw2G7fv7YPpL0xM4dC7ofPGvaN5uccB1w23tI+qGkjqQtkr4v6aSIuH2c/T8v6esR\n8dXRz7EOeAVYB3yuYB3wMbAOOACgDjNZB3zCGfCIeML26ZJWS9pd0qURcbvt04rnL5lxbwEAAIA5\naNJ1wCPi6og4PCIOi4iPFY9dMtbgOyLeOdbst0QNeCbqaPNQA56LczcP9Z55yDYP2eYh22bhTpgA\nAABAhSasAS8TNeAVaF8NeJusk7S8pLaoAR8DNeAAgDrMpAacGXA00bq6O5Bgcd0dAAAAzTDZMoSl\n6Xa7016iBVNje0nL7ihY1kzxrJ3RtxTiLJQ9oz+v+LRj4J0lHXRhxAV196ONNq9aMcSd73KQbR6y\nzUO2zVLZABxAaTbV3YGyzJcW1d0HAACqRg14m7SnBrxtWNd8fKXVs1MDDgCoAzXgAAAAQMNRA94C\nLawBb4ySasAxjvOlReeU19w2SQvKa26wbbnmskUHHHfKbXX3ozBPLSqdopY2D9nmIdtmoQYcQCus\nP3vp6rr70CR+X2fJ/atXrq27H9KT5UEAgEJlJShDQ9yAKQuz33mY/c51jtSUGdrW4bqQh1nEPGSb\nh2ybhRpwAAAAoEKVDcC7Xd54ZbG9pO4+tNUZEh/dJDqfZQjTcF3Is3nVCq4LScg2D9k2CzPgAAAA\nQIWoAW8Baj3zUAOeixrwPFwX8lBLm4ds85BtszADDgAAAFSIGvAWoNYzDzXguagBz8N1IQ+1tHnI\nNg/ZNgvrgLdLm25E0qobdwAAAIxwRFRyoOHh4eh0Oq7kYBh89jJJW+vuRknWFH8eW2svmmmBIq6o\nuxPIVdyIpy3/ngHgaT7+qlgz3TEuNeAAAABAhagBbwFqPfNQA56LGvA8XBfyUEubh2zzkG2zMAMO\nAAAAVIgacDQTNeBzxSGSHq+7E2PYpojVdXeiLagBB9BmM6kBZxUUAHVq6ko3bVlNCADQQNSAtwC1\nnnmoAc9Fvnm4LuShljYP2eYh22ahBhwAAACoUGUD8KEh3nhliYi1dfehrS6W+OgmEfnm4bqQ5+AT\nT+e8TUK2eci2WZgBBwAAACpEDXgLUOuZhxrlXOSbh+tCHmpp85BtHrJtFmbAAQAAgApRA94C1Hrm\noUY5F/nm4bqQh1raPGSbh2ybhRlwAAAAoELUgLcAtZ55qFHORb55uC7koZY2D9nmIdtm4U6YQHXW\nTL7LlKyTtLyktgAAQMWoAW8Baj3zlFSjvK6ENvotLrm92lADnofrQh5qafOQbR6ybRZmwIF8Zc5W\nlzWLDgAAakINeAtQ65mHGuVc5JuH60IeamnzkG0esm0WVkEBAAAAKkQNeAtQ65mHGuVc5JuH60Ie\namnzkG0esm0WasABANm2SVpQUlvzJG0qqS0AqEVlA/But6tOp1PV4eYU20uY7cpxhjTELG0e8s3T\npOvC+rOXri6rraMuuHZZWW3N1OZVK4aYTcxBtnnItlmYAQeAZ5onu/aB3mx8WFoku6xZ59napojS\nBuEAMOgqG4BTA56nKbNcbcTsbK4G5zvwJQ7nSFvr7kOfprwRKAWziHnINg/ZNguroAAAAAAVYh3w\nFmC93zysU52LfPOQbR7WU85DtnnItlmYAQcAAAAqxDrgLUANeJ4G1yi3AvnmIds81NLmIds8ZNss\nzIADAAAAFaIGvAWoAc9DHW0u8s1Dtnmopc1DtnnItlmYAQcAAAAqRA14C1ADnoc62lzkm4ds81BL\nm4ds85BtszADDgAAAFSIGvAWoAY8D3W0ucg3D9nmoZY2D9nmIdtmYQYcAAAAqBA14C1ADXge6mhz\nkW8ess1DLW0ess1Dts3CDDgAAABQoSkNwG0fb/sO23fafv8Yz/9ftm+xfavt79p+5eh9qAHPQw14\nHupoc5FvHrLNQy1tHrLNQ7bNMukA3PbuklZIOl7SIkkn2f7lUbvdLel1EfFKSR+R9FdldxQAAABo\ng6nMgB8taUNEbIyIX0i6QtKJ/TtExPci4pFi80ZJLx7dCDXgeagBz0MdbS7yzUO2eailzUO2eci2\nWaYyAD9Q0r192/cVj43nXZKumk2nAAAAgLaaygA8ptqY7WMl/b6kZ9SJUwOehxrwPNTR5iLfPGSb\nh1raPGSbh2ybZY8p7HO/pIP6tg9Sbxb8aYovXn5W0vERsW3089ddd53OOuuslZI2Fg89LKk7Uj4x\nMohkm21JOl9a9FNpwcjH6CODiaq3R9R1/DG2m9aftuXbmu2udJga0p/zpUXn2kvKuj6MDCRGPlKv\nenvnlg2HbV61orbjs832TLZHNKU/g7y9c8uGw3b9/LF9JOmJnTsWdD94ljqdjqbDERNPcNveQ9IP\nJXUkbZH0fUknRcTtffscLOlaSW+PiHVjtTM8PBydTsfT6h3mLnuZpK11d6OB1hR/HltrL4DpOUTS\n42U0tHLJslf+xQnv/lYZbQFAGT7+qlgz3THupDPgEfGE7dMlrZa0u6RLI+J226cVz18i6RxJz5P0\nGduS9IuIOHq6vwCAKVsz+S4DYZ2k5XV3Auk2ldXQc3du/7Wy2gKAukxpHfCIuDoiDo+IwyLiY8Vj\nlxSDb0XEf42I50fEkcXPMwbf1IDnoQY8TwPraMf8hGlQrZUW192Htmrgudsa1NLmIds8ZNssU6kB\nB9AcbZotbsssPir07Md3PippQZ192GPv+fvV3Ye2Its8ZNssk9aAl4UacEyL/Ub1ypraYp5K/Bi+\nJahnx0wsUMQVdXcCAEbMZIzLDDiaKWJ13V0oVe9LpQAAAFOrAS8DNeB5qAHPQ7a51tbdgRajBjwP\n14U8ZJuHbJulsgE4AAAAgAoH4ENDTMZkGbk5BcpHtrmW1N2BFhu5CQ7Kx3UhD9nmIdtmYQYcAAAA\nqBA14C1AXVcess21tu4OtBg14Hm4LuQh2zxk2yzMgAMAAAAVoga8BajrykO2uZbU3YEWowY8D9eF\nPGSbh2ybhRlwAAAAoELUgLcAdV15yDbX2ro70GLUgOfhupCHbPOQbbMwAw4AAABUqLJb0VMDnoe6\nrjxkm2tJ7481tXZibOskLa+7E7NBDXgergt5yDYP2TYLM+AA6rKu7g5MYHHdHQAAtFdlM+Ddbled\nTqeqw80ptpfwzjYH2aZafoY01MCZ2ibOyE9bQ7NtBa4Lecg2D9k2CzPgAAAAQIVYB7wFeEebh2xz\nMUObh2zzcF3IQ7Z5yLZZmAEHAAAAKsQ64C3A2p55yDYXa1XnIds8XBfykG0esm0WZsABAACAClED\n3gLUdeUh21zUKech2zxcF/KQbR6ybRZmwAEAAIAKUQPeAtR15Skx222SFpT0c0hJfaoddcp5yDYP\n19w8ZJuHbJulshvxAHNaxOrS2rKXldYWAACoHDXgLUBdVx6yzUWdch6yzcN1IQ/Z5iHbZqEGHAAA\nAKgQNeAtQF1XHrLNRZ1yHrLNw3UhD9nmIdtmoQYcAMa2pu4OzMZv1d2Bp1snaXndnQCApqAGvAWo\n68pDtrkaWqe8ru4OlGFJ3R14usV1d6BMXBfykG0esm0WZsAB4OmYqS3XQH+SAAAZqAFvAeq68pBt\nLuqU85BtHq4Lecg2D9k2C6ugAAAAABWiBrwFqOvKQ7a5GloD3gpkm4frQh6yzUO2zcIMOAAAAFAh\nasBbgLquPGSbizrlPGSbh+tCHrLNQ7bNwgw4AAAAUCFqwFuAuq48ZJuLOuU8ZJuH60Iess1Dts3C\nDDgAAABQIWrAW4C6rjxkm4s65Txkm4frQh6yzUO2zcIMOAAAAFChym5FTw14Huq68jQ0222SFtTd\niTJc3PtdkIAa8DwNvS60AtnmIdtmqWwADqAkEavr7kJp7GV1dwEAgKpRA94C1HXlIdtc50uL6u5D\nW1EDnofrQh6yzUO2zcIMOACgCmtKaudmSVeU1BYA1IJ1wFuAuq48ZJvrHOm2uvvQVg2qAV9XcntH\nltzetHFdyEO2eci2WZgBBwBkWl5iW2XNogNAragBbwHquvKQbS5qwPNQA56H60Iess1Dts3CDDiA\n2myXdqiZSyrOk7Sp7k4AANqJdcBbgLquPGSb68KIC+ruw5hasDxig2rAW4frQh6yzUO2zcKdMAEA\nAIAKUQPeAtR15SHbXOSbhxrwPJy3ecg2D9k2CzPgAAAAQIVYB7wFqOvKQ7a5yDcPNeB5OG/zkG0e\nsm0WZsABAACACk26Cort4yVdJGl3SX8dEZ8YY5+LJb1J0k5Jp0bEzaP36Xa76nQ6s+8xnsH2Et7Z\n5iDbXA3Od5uauTzilJ0vLWrQnUZbtaxjg8/bgUe2eci2WSYcgNveXdIKScdJul/SettXRsTtffuc\nIOmwiHip7VdL+oykxaPb2rBhQ6kdx9MMSVpbdydaimxzNTPfiNV1d2G2zrX/v3Mirqi7H5Jasazj\nKM08b9uBbPOQbZKZTDJPVoJytKQNEbExIn4h6QpJJ47a562SviBJEXGjpH1tv3B0Q48++ui0OoZp\n2bfuDrQY2eYi3zxkm4ds85BtHrJNcsstt0z7NZMNwA+UdG/f9n3FY5Pt8+Jp9wQAAACYAyarAY8p\ntuPJXrd169YpNoUZWFh3B1psYd0daLmFdXegxRbW3YE+5dbU21P9f1OKd/T6cG6dfWgrss1DtolO\nPnnaL3HE+Ncx24slnRcRxxfbyyXt6v8ipu2/lLQ2ilpD23dIen1EPNDf1h/8wR9EfxnKEUccwdKE\nJel2u2SZhGxzkW8ess1DtnnINg/Zlqfb7T6t7GTvvffWZz7zmdGT0ROabAC+h6QfSupI2iLp+5JO\nGuNLmKdHxAnFgP2iiHjGlzABAAAATFKCEhFP2D5d0mr1liG8NCJut31a8fwlEXGV7RNsb5D0qKR3\npvcaAAAAGFATzoADAAAAKFf6nTBtH2/7Dtt32n5/9vHmGtsbbd9q+2bb36+7P4PM9udsP2D7X/oe\n28/2t23/yPa3bLOM0wyMk+15tu8rzt2bi5t+YZpsH2R7je3/ZfsHts8oHufcnaUJsuXcnSXbv2T7\nRttd27fZ/ljxOOdtCSbIl3O3JLZ3LzL8erE9rXM3dQa8uJHPD9V3Ix+NqiHH7Ni+R9KvRMRDdfdl\n0Nl+raR/k3RZRLyieOzPJP00Iv6seAP5vIj44zr7OYjGyfZcSTsi4s9r7dyAs71A0oKI6NreR9I/\nSfpN9coBOXdnYYJsf0ecu7Nm+9kRsbP4vtkNks5S794inLclGCffjjh3S2H7vZJ+RdJzIuKt0x0v\nZM+AT+VGPpi9aX3zFmOLiOvVWy6t35M3mir+/M1KO9US42Qrce7OWkRsjYhu8fd/k3S7evdn4Nyd\npQmylTh3Zy0idhZ/3Uu975ltE+dtacbJV+LcnTXbL5Z0gqS/1lN5TuvczR6AT+VGPpidkHSN7Zts\nv7vuzrTQC/uW1HxA0jPu8opZ+UPbt9i+lI+aZ8/2QklHSrpRnLul6st2XfEQ5+4s2d7Ndle983NN\nRPwvcd6WZpx8Jc7dMnxK0tmSdvU9Nq1zN3sAzjc88/16RBwp6U2S/p/io34kiF69Fud0eT4j6SWS\nhiT9WNIn6+3OYCtKJL4i6f+NiB39z3Huzk6R7d+rl+2/iXO3FBGxKyKG1Lt79utsHzvqec7bWRgj\n3yXi3J0122+R9JOIuFnjfJowlXM3ewB+v6SD+rYPUm8WHCWJiB8Xfz4o6Wvqlf2gPA8UdaCy/SJJ\nP6m5P60RET+Jgnof43HuzpDtPdUbfF8eEf9QPMy5W4K+bL84ki3nbrki4hFJ31CvnpbztmR9+f4q\n524pXiPprcV38L4kaantyzXNczd7AH6TpJfaXmh7L0m/K+nK5GPOGbafbfs5xd/3lvQGSf8y8asw\nTVequINv8ec/TLAvpqG4QI34L+LcnRHblnSppNsi4qK+pzh3Z2m8bDl3Z8/2/iPlD7afJek3JN0s\nzttSjJfvyACxwLk7AxHxgYg4KCJeImmZpGsj4mRN89xNXwfc9pskXaSnbuTzsdQDziG2X6LerLfU\nu6nS35DvzNn+kqTXS9pfvfqtcyStkvR3kg6WtFHS70TEw3X1cVCNke25kpao9zFoSLpH0ml99XOY\nItvHSPqOpFv11Eeey9W7czHn7iyMk+0HJJ0kzt1Zsf0K9b6otlvxc3lEXGB7P3HeztoE+V4mzt3S\n2H69pDOLVVCmde5yIx4AAACgQuk34gEAAADwFAbgAAAAQIUYgAMAAAAVYgAOAAAAVIgBOAAAAFAh\nBuAAAABAhRiAA0AD2D7Y9o7i5i9lt73S9keKv7/W9h1lHwMAMHUMwAGgJrY32l4qSRGxOSKeEzk3\nZ4jiRxFxfUS8POEYAIApYgAOAPUJSaXPeI+jquMAACbBABwAamD7cvVuWfz1ovTkbNu7bO9WPL/W\n9kdsf7d4/krb+9v+G9uP2P6+7UP62nu57W/b/lfbd9h+2zjHXWL73r7tjbbPtH2L7YdtX2F7Xt/z\nb7Hdtb2t6Msr8lIBgLmBATgA1CAiTpa0WdJbIuI5kr48xm6/K+ntkg6UdKik70m6VNJ+km6XdK4k\n2d5b0rclfVHSCyQtk/Rp2788la5IepukN0p6iaRXSjq1aPfI4njvLo55iaQrbe817V8YAPAkBuAA\n0Ewh6fMRcU9EbJd0taQfRcS1EfHv6g3Yjyz2fYukeyLiCxGxKyK6kr6q3sB6Ki6OiK0RsU3S1yUN\nFY+/R9IlEbE+ei6T9LikxeX8igAwN+1RdwcAAON6oO/vj0n6yajtfYq/HyLp1ba39T2/h6TLpnic\nrX1//5mkA/raPcX2H/Y9v6ekF02xXQDAGBiAA0B9prPiyUT7bpZ0XUS8oeRjbZb0pxHx0Wm8FgAw\nCUpQAKA+D6hX2z0ej/P30b4h6WW23257z+LnKNsv73vtdFZBGdn3s5L+m+2j3bO37Tfb3meiFwMA\nJsYAHADq8zFJH7L9kKT/U8+cpY5Rfx/z+YjYIekN6n358n5JPy7a3muc1040G96/Zvg/qfcFzBWS\nHpJ0p6RTpvB7AQAm4Jx7PgAAAAAYCzPgAAAAQIUYgAMAAAAVYgAOAAAAVIgBOAAAAFAhBuAAAABA\nhRiAAwAAABViAA4AAABUiAE4AAAAUCEG4AAAAECF/jd15AZ+z4u6xAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "kmf = KaplanMeierFitter()\n", "kmf.fit(T[treatment], timeline=t, event_observed=E[treatment], label='With treatment')\n", "ax = kmf.plot()\n", "\n", "kmf.fit(T[~treatment], timeline=t, event_observed=E[~treatment], label=\"Without treatment\")\n", "kmf.plot(ax=ax, c=\"r\")\n", "ylim(0,1.05)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Hazards, too" ] }, { "cell_type": "code", "execution_count": 21, "metadata": { "collapsed": false }, "outputs": [], "source": [ "naf = NelsonAalenFitter()" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuAAAAGCCAYAAAChCnbFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcXHWZ7/Hvk4QUEEBQViEYZZHCUWA0gCsxxVWMgDNz\n0bG8iMtc4ToybhBHnbkal1HHgAIuDF7RAGorwnVB4ap0IIALRqGDmkaICGFfZCfQGPLcP+p0KDq9\n13nq/PrXn/fr1a+qU3VO1dNfj+RXv37qd8zdBQAAAKA7ZlRdAAAAADCdMAAHAAAAuogBOAAAANBF\nDMABAACALmIADgAAAHQRA3AAAACgi8Y1ADezmWZ2tZldMMLzp5nZ9Wa2yswOKLdEAAAAIB/jnQF/\nj6TVkjZZNNzMFkna0933knSspNPLKw8AAADIy5gDcDPbTdIiSV+VZMPscqSksyTJ3a+UtK2Z7VRm\nkQAAAEAuxjMD/nlJiyVtGOH5XSXd3LZ9i6TdOqwLAAAAyNKoA3AzO1zSXe5+tYaf/d6465Btrm8P\nAAAADGPWGM+/RNKRRZ/35pK2MbOz3f2Ytn1ulTS3bXu34rGnOPLII/2xxx7TzjvvLEmaM2eO9txz\nT+2///6SpL6+PkliexLbg/dTqSen7cHHUqknt+3Bx1KpJ6ftNWvW6Kijjkqmnpy2zzvvPP79Ctrm\n3zP+eztVtiVp1apVuuOOOyRJe+yxh04//fTRJqufwtzHN1ltZodIOtHdjxjy+CJJx7v7IjM7WNIp\n7n7w0OOPOeYYP/vss8ddGMbPzJa4+5Kq68gR2cYi3zhkG4ds45BtHLKNNdFx7lgz4EO5JJnZcZLk\n7me4+4VmtsjM1kh6RNLbhjtw8BMCQsyruoCMzau6gMzNq7qAjM2ruoCMzau6gIzNq7qAjM2rugA8\nadwDcHdfIWlFcf+MIc8dX3JdAAAAQJa6diXMV7/61d16q+loWdUFZGxZ1QVkblnVBWRsWdUFZGxZ\n1QVkbFnVBWRsWdUF5Gy//fab0P5dG4APNq+jfO5+adU15IpsY5FvHLKNQ7ZxyDYO2caa6Dh3oj3g\nk9bX16dGo7HJ4729vTMkvU/SNmL5wkm57rrr5u299943Vl1HokzSg5I+32g0RlrLfuSDzRbwH604\n5BuHbOOQbRyyjUO2aenaAHwU75N0caPRWFV1IVPVoYceyv+pRtHb27ufWufZyVXXAgAAkEILyjYM\nvjvD4Ht0xfm1zWSOJdtY5BuHbOOQbRyyjUO2aenaAHwUtJ2gGzjPAABAEro2AG+/chDKZWYLqq4h\nV2Qbi3zjkG0cso1DtnHINi0pzIBnwcxebmbXjvL8PDPbYGZkDgAAMI2l0AOeJDP7kJldOOSx60d4\n7A3ufrm779P2+I1mtrCD919mZp8Yz76T6esyswVmdvOECyuJmV1qZv9U1fuPFz1zscg3DtnGIds4\nZBuHbNPCbOzIVkh6iZmZJJnZLmqtGrP/4Cx28dgeki4b5nhXawm8EGaWwgo2naAnGwAATEv0gI/s\nN5I2kzQ4df9ySZdIum7IY2vc/Y72GWUzO0fS7pIuMLOHzOzEttc92sxuMrO7zezDw72xmR0r6U2S\nPlAc/4Pi8RvN7ANmdo2kh8xshpkdbGa/N7P7zKzPzA5pe523mdlqM3vQzP5UvK7MbI6kiyQ9s3j9\nB81sFzNbYmbfNbNziseuMbO9ir8G3FnU/d/aXv9pZnammd1mZreY2SfaPpy81cyuMLOlZnavmd1g\nZocVz/1Hkd0Xi/c/bbL/I0WjZy4W+cYh2zhkG4ds45BtWpgBH4G7Py7pSkmDA9pXSLpc0hXF/cHH\nNpn9dvc3S1or6XB339rdT2p7+qWS9pbUkPQRM9tnmOO/Iumbkv6zOP51bU+/UdJrJG0raRdJP5J0\ntrtvJ+lESeeb2TOKfe+U9Fp330bS2yR93swOcPdHJB0m6bbi9bdx99uLYw6XdLak7SRdLelnxePP\nlPQJSWe01bJM0uNq/RXgAEmvkvQ/254/UNK1kp4h6bOSzix+v38rsnxX8f7vHpoBAABArrrWxjCZ\nHvD5S5eX1qawcvHCybSDrFBrkH2KpJdJ+ryk2yQdVzz2ckknjXj08D7m7gOSrjGzVZL2U2uQOpyh\nNbuk09z9Vkkys6MlXejun5Ukd7/YzH4j6bVqDco39qu7+2Vm9tOi5quHee1Bl7n7z4rXP0/SP0j6\njLu7mX1H0lfMbBtJW6j4IODuj0l61MxOkfQOSV8pXusmdz+zeK2zJX3ZzHZ097tG+P2SQ89cLPKN\nQ7ZxyDYO2cYh27RM9T7iaJdJepeZbSdpB3f/k5ndLems4rHnafj+79Hc0XZ/naQ5Ezy+/YuTz5L0\nejM7ou2xWZKWS5KZvUbSRyXtpdZfO7aUdM0Yr39X2/1HJd3j7t62LUlbSdpNrRad24s2eRXvsbbt\n+I2/q7uvK/bbqu096AMHAADTTtcG4H19fWo0GhM6ZpKz1mX6laSnqTWr+3NJcvcHzew2Sceq1cJx\n0wjHdjq4HOn49sfXSjpH0reGfrI1s5qk8yUdLekH7v6EmX1PT846D/f6E6n5ZkkDkp7h7hsmcNxk\n3qsyZraAWYM45BuHbOOQbRyyjUO2aaEHfBTu/qhaX8Z8v546031F8diKUQ6/U63e6LGM9CHjTknP\nGePYb0g6QtKLzGymmW1efBl0V0mzi597JG0oZsNfNeT1n1G0k4xVyyaKnvGfSvqcmW1dfCF0DzN7\nxVjHtr3/ePIBAADICuuAj22FpB3UGnQPulzS9tq0/aR9VvfTkv69WJ3k/cM8P9wx7c6UtG9x/P8d\nbgd3v0XS69T64uRdas2InyDJ3P0hSe+WdK6keyU1Jf2g7dhrJfVIuqFYpWSXopah9Yy2fYxag/zV\nxXt8V9LObfuNduypko4q3vuUYRNIALMFscg3DtnGIds4ZBuHbNNiT7b3xurt7fVGo7HJDGtvb++S\nRqOxpCtFYNriPAMAAFFGGueOhHXAM8DannHINhb5xiHbOGQbh2zjkG1a6AEHAAAAuoge8AzQ1xWH\nbGORbxyyjUO2ccg2DtmmhRlwAAAAoIvoAc8AfV1xyDYW+cYh2zhkG4ds45BtWpgBBwAAALqoa1fC\npAc8Dn1dccg2FvnGyS3bnlp9hqRXSNqu6lq+NXsf9dTqf191HTki2zhkG2vHC784of27NgAHAKAD\nr5M07EXJAGCq6doAvK+vT41Go1tv13Vm9nJJ/8fd9xnh+XmSbpA0y903lPzeC3Kb7UoF2cYi3zgZ\nZrtHcbtG0u+rLOTKDQ9uf9CMbe6psoZckW0csg33dxPZmRnwEZjZhyS93N0XtT12vaTrh3ns39z9\nXEn7tD1+o6S3u/vyLpT7QTNruPv/HmkHM9sgaU93vyGiADNbJunm0WqIEvnhBkAytiluv9Ec6P9Y\nlYW8yWyB+62XVllDrsg2DtnG6u3tndCl5VkHfGQrJL3EzEySzGwXtT6w7G9mM9oe20PSZcMc75LG\nfUnSDt0xzv1GrMfMcvgwVnremc0gJod842SY7dOK2wcrrUJZZpsMso1DtmlhFZSR/UbSZpIGPzm8\nXNIlkq4b8tgad7/DzBaY2c2SZGbnSNpd0gVm9pCZndj2ukeb2U1mdreZfXjwQTOrmdkpZnZr8fN5\nM5tdPPdWM7u8vTgz22Bme5jZsZLeJOkDxXv9YOgvYmaDHxBWFfu8vqj3FjP7gJndLulMa/mgma0x\ns3vM7Dtmtl3b63zXzG43s/vNbIWZ7Vs8PmwNZnajmZ1oZtcUj59pZjuZ2UVm9oCZ/czMtm17/YPN\n7Bdmdp+Z9ZnZIW3PXWpmHzezK8zsQTP7iZk9o3h68Pe7v3ifg0b53xXA1DQ4A/5ApVUAQAlYB3wE\n7v64pCslDQ4CXyHpcklXFPcHH9tk9tvd3yxpraTD3X1rdz+p7emXStpbUkPSR8zsucXj/ybpQEn7\nFT8HSvr3scv0r0jqlfSfxXu9bpidBut9QbHPd4vtndRaUWB3ScdJerekI4vfaxdJ90n6UttL/VjS\nnpJ2kHSVpG8Wr/+V4v7QGlzSPxS/63MlHS7pIkkflLSjWuffuyXJzHaV9CNJH3f37SSdKOn8tkG2\nJDUlvbU4dnaxj9T6ICRJTyve/8oxchs31k2NRb5xMsx2cABe+Qx4htkmg2zjkG1akm476KnVJ9RP\nM5rmQP9k2hNWqDUYPUXSyyR9XtJtag1WT1Fr4HfSiEcP72PuPiDpGjNbpdZg+49qzSAf7+73SJKZ\nfUzSGZI+Ms7Xnczvt0HSR939r5L+ambHFTXc1lbDTWZ2tLtvcPdlG9+s9dx7zGxrd39olBq+4O53\nF8dcLulOd19VbH9PrcG5JB0t6UJ3/3+S5O4Xm9lvJL1W0tlqDea/7u5rimPPVevDwmR/dwBTSzIt\nKADQKXrAR3eZpJcVbRg7uPufJP1Srd7w7SQ9T8P3f4+mvV97naStivvPlHRT23Nri8cm+poTcXcx\n0z9onqTvFS0g90laLWm9pJ3MbKaZfaZoT3lA0p+LY7Yf4z3ubLv/6JDtx/Tk7/8sSa8ffO/i/V8q\naee2/dt/z0fbjg1Dz1ws8o2TYbbJtKBkmG0yyDYO2aZlzBlwM9tcrZngmlp/9v+Bu39oyD4LJP1A\nrZUoJOl8d/9kp8VNcta6TL9Sa9blHZJ+Lknu/qCZ3SbpWEm3uftNIxw70dn729QaAPcX27sXj0nS\nI5K2HNzRzHbWU032LwVDj1sr6W3u/suhO5rZm9WacW64+01F7/a9enL2ebw1jPS/6VpJ57j7seN8\nnXal/aUEQLKSaUEBgE6NOQPu7o9JeqW77y/pBZJeaWYvG2bXFe5+QPGzyeB7qvWAS5K7P6rWlzHf\nr6fOdF9RPLZilMPv1JPr1o5Hj6R/N7PtzWx7tVpPzimeWyXpeWa2X/GBaMmQY2uSnjPG64+nnv+S\n9Ckz212SzGwHMxts89hK0oCke81sjqRPDfP6Y9Uwmm9IOsLMXlXMtm9efFF017Z9Rhq8361WO81E\n8h4XeuZikW+cDLNNpgUlw2yTQbZxyDYt42pBcfd1xd3ZkmaqNfM5VNWz1VFWqPWlwyvaHrtcrdaL\noe0n7TOxn1ZrQH2fmb1/mOeH+qRag/1rip/fFI/J3a+T9HFJF6vVL375kNe6UNK+xXuNdKW4JZLO\nKvY5qjh+aD2nSvqhpJ+a2YNqtdscWDx3tlotMreqdRGMXw45/sxx1KAhx2yswd1vUetKdx+WdJda\nM+In6Knn1UjHrpP0H5J+Xrz/gQKQm2RaUACgU+Y+9l/vrbXu9VVqzTCe7u4fGPL8IWpdIvgWtQZo\nJ7r76vZ9ent7vdFobDJI7+3tXdJoNJZM+jcAxoHzDJi6emr1WZL+qtaH7pnNgX7azgAkZaRx7kjG\nOwO+oWhB2U3SK4b5M8ZVkua6+36SviDp++MtAACAMWzs/2bwDSAHE1oFxd0fUGst6BcNefyhwTYV\nd79I0mZm9vT2fU499VSZ2TIzW1L8vLd9IF/0+7I9ie3B+6nUk+L2ddddN28yxw/NOJXfJ5dt8g3d\nfm9i9Ux6+9N/Xfuq1RsekYr2k6rrMf79CtsevJ9KPTltD8246nqm+nbxs8RaY9tlE/2u45gtKNb6\nQuB6d7/fzLaQ9BO11rLubdtnJ0l3ubtbq//2XHef1/46J598sp9wwgm0oAQwswUsLzS6yZ5nZBuL\nfOPklG1Prf4Ctb6M/vvmQP/zq64np2xTQ7ZxyDbWRFtQxnMhnl3U+vLeDLVmzM9x915rXbRF7n6G\npKMkvdPM1qu1tvUbh77IFF0HfErg/1BxyDYW+cbJLNtkVkCRsss2KWQbh2zTMuYA3N1/J+lvh3n8\njLb7X9JTL1k+EbmunoK0cJ4BUxcroADISteuhDlKb8ys3t5eBkcdaO9NwqaK82s8f+3ZBNnGIt84\nmWWb1EV4Mss2KWQbh2zT0rUB+CgulPSWqotA1t6i1peHAUxNSbWgAECnJjUrOBkj9YA3Go2f9/b2\nPre3t/fjal3NEBN08cUXq7e3d0HVdSRqhqQbG43GLyZzMD1zscg3TmbZJjUDnlm2SSHbOGSblq4N\nwEfTaDS+VnUNAIBk0QMOICsp9ICjQ/R1xSHbWOQbJ7Nsk2pBySzbpJBtHLJNSwo94AAAjCapFhQA\n6FTXBuCsAx6Hvq44ZBuLfONklm1SLSiZZZsUso1DtmlhBhwAkDpmwAFkhR7wDNDXFYdsY5FvnMyy\npQd8miDbOGSbFmbAAQCpS6oFBQA6RQ94BujrikO2scg3TmbZJtWCklm2SSHbOGSbFmbAAQCpS6oF\nBQA6RQ94BujrikO2scg3Ti7Z9tTqsyVtLukJSesqLkdSPtmmiGzjkG1amAEHAKRs6+L2weZAv1da\nCQCUhB7wDNDXFYdsY5FvnIyyTa79JKNsk0O2ccg2LcyAAwBSxgooALJDD3gG6OuKQ7axyDdORtkm\ntQKKlFW2ySHbOGSbFmbAAQApS64FBQA6RQ94BujrikO2scg3TkbZJteCklG2ySHbOGSbFmbAAQAp\nS64FBQA6RQ94BujrikO2scg3TkbZJteCklG2ySHbOGSbFmbAAQApS64FBQA6RQ94BujrikO2scg3\nTkbZJteCklG2ySHbOGSbFmbAAQApS64FBQA6RQ94BujrikO2scg3TkbZJjcDnlG2ySHbOGSbFmbA\nAQApowccQHboAc8AfV1xyDYW+cbJKNvkWlAyyjY5ZBuHbNPCDDgAIGXJtaAAQKfoAc8AfV1xyDYW\n+cbJKNvkWlAyyjY5ZBuHbNPCDDgAIGXJtaAAQKfoAc8AfV1xyDYW+cbJIdueWr0mabak9ZIeq7ic\njXLINlVkG4ds08IMOAAgVRvbT5oD/V5pJQBQInrAM0BfVxyyjUW+cTLJNskvYGaSbZLINg7ZpmXU\nAbiZbW5mV5pZn5mtNrNPj7DfaWZ2vZmtMrMDYkoFAEwz9H8DyNKs0Z5098fM7JXuvs7MZkm6wsxe\n5u5XDO5jZosk7enue5nZQZJOl3Tw0NeiBzwOfV1xyDYW+cbJJNvkVkCRssk2SWQbh2zTMmYLiruv\nK+7OljRT0r1DdjlS0lnFvldK2tbMdiqzSADAtJRkCwoAdGrMAbiZzTCzPkl3SrrE3VcP2WVXSTe3\nbd8iabehr0MPeBz6uuKQbSzyjZNJtkm2oGSSbZLINg7ZpmU8M+Ab3H1/tQbVrxjhf0AbelgJtQEA\nprckW1AAoFOj9oC3c/cHzOzHkl4k6dK2p26VNLdte7fisadYs2aNzGyZpBuLh+6X1DfYkzQ4sGd7\n4tvufmlK9bDNNttpbA9KpZ6Jbn9r9j7bSNKlT9y/7ZvMFlRdT3sPrSVUT07b/HvG9lTZLiyQNE+S\nTjrpJDUaDY2XuY88WW1m20ta7+73m9kWkn4i6WPu3tu2zyJJx7v7IjM7WNIp7r7JlzB7e3u90WgM\nnSkHAGBYPbX6ZyT9q6R/aw70f6rqegBgJBMd547VgrKLpOXW6gG/UtIF7t5rZseZ2XGS5O4XSrrB\nzNZIOkPSPw/3QvSAxxk624XykG0s8o2TSbZJtqBkkm2SyDYO2aZlrGUIfyfpb4d5/Iwh28eXXBcA\nAKyCAiBLXbsSJuuAx2nvS0S5yDYW+cbJJNskV0HJJNskkW0csk1L1wbgAABMUJItKADQqa4NwOkB\nj0NfVxyyjUW+cTLJNskWlEyyTRLZxiHbtDADDgBIVZItKADQKXrAM0BfVxyyjUW+cTLJNskWlEyy\nTRLZxiHbtDADDgBITk+tbkq0BQUAOkUPeAbo64pDtrHIN04G2W4uaTNJjzcH+geqLqZdBtkmi2zj\nkG1amAEHAKSI2W8A2aIHPAP0dcUh21jkGyeDbJPs/5ayyDZZZBuHbNPCDDgAIEWsgAIgW/SAZ4C+\nrjhkG4t842SQbbItKBlkmyyyjUO2aWEGHACQomRbUACgU/SAZ4C+rjhkG4t842SQbbItKBlkmyyy\njUO2aWEGHACQomRbUACgU/SAZ4C+rjhkG4t842SQbbItKBlkmyyyjUO2aWEGHACQImbAAWSLHvAM\n0NcVh2xjkW+cDLKlB3waIts4ZJsWZsABAClKtgUFADpFD3gG6OuKQ7axyDdOBtkm24KSQbbJIts4\nZJsWZsABAClKtgUFADpFD3gG6OuKQ7axyDdOBtkm24KSQbbJIts4ZJsWZsABAClKtgUFADpFD3gG\n6OuKQ7axyDdOBtkm24KSQbbJIts4ZJsWZsABAEnpqdVNzIADyNisbr0RPeBx6OuKQ7axyDdOKtn2\n1Orvk/T8CR42s/h5rDnQ/3j5VXUmlWxzRLZxyDYtXRuAAwCml55a/TmSPtfBS6wtqxYASEnXBuB9\nfX1qNBrdertpxcwW8Mk2BtnGIt84iWT77OK2X9JJkzj+shJrKU0i2WaJbOOQbVqYAQcARJlb3F7V\nHOj/WqWVAEBCWAc8A3yijUO2scg3TiLZ7lbc3lxpFSVLJNsskW0csk0Lq6AAAKIMzoDfUmkVAJAY\n1gHPAGt7xiHbWOQbJ5Fss5wBTyTbLJFtHLJNCzPgAIAogzPgWQ3AAaBT9IBngL6uOGQbi3zjJJJt\nli0oiWSbJbKNQ7ZpGXMAbmZzzewSM/uDmf3ezN49zD4LzOwBM7u6+Pn3mHIBAFNBT62+laRtJT0m\n6Z6KywGApIxnBvyvkt7n7s+TdLCkd5lZfZj9Vrj7AcXPJ4c+SQ94HPq64pBtLPKNk0C2g/3ftzQH\n+r3SSkqWQLbZIts4ZJuWMQfg7n6Hu/cV9x9W64IKzxxmVyu5NgDA1JVl+wkAlGFCPeBmNk/SAZKu\nHPKUS3qJma0yswvNbN+hx9IDHoe+rjhkG4t84ySQbZYroEhJZJstso1DtmkZ95UwzWwrSedJek8x\nE97uKklz3X2dmb1G0vcl7V1emQCAKYYVUABgBOMagJvZZpLOl/QNd//+0Ofd/aG2+xeZ2ZfN7Onu\nfu/g46eeeqoOPfTQZZJuLB66X1Lf4Ceywd4ktie+3d7XlUI9OW0PPpZKPbltDz6WSj2Zbe/v7qdU\n9f7vmrnLgS+d+TRJuiWRPMrcfq/494t/z6bY9uBjqdQz1bcLCyTNk6STTjpJjUZD42Xuo383xsxM\n0lmS/uLu7xthn50k3eXubmYHSjrX3ee173PyySf7CSecQJ94ADNbMHhioFxkG4t841SdbU+tfpGk\nwyQd0Rzo/1FVdUSoOtuckW0cso3V29vrjUZj3OPc8cyAv1TS0ZKuMbOri8c+LGl3SXL3MyQdJemd\nZrZe0jpJbxz6IvSAx+H/UHHINhb5xkkg22y/hJlAttki2zhkm5YxB+DufoXG+LKmu39J0pfKKgoA\nMOVl+yVMAOhU166EyTrgcYb0I6FEZBuLfONUmW1Prb61pKdJelTSvWPsPuVw3sYh2zhkm5auDcAB\nANPGxvaT3C7CAwBl6NoAnB7wOPR1xSHbWOQbp+Jss24/4byNQ7ZxyDYtzIADAMrGGuAAMAp6wDNA\nX1ccso1FvnEqzjbbFVCkyrPNGtnGIdu0MAMOAChb1i0oANApesAzQF9XHLKNRb5xKs426xYUzts4\nZBuHbNPCDDgAoGxZt6AAQKfoAc8AfV1xyDYW+capONusW1A4b+OQbRyyTQsz4ACA0vTU6k+TtLWk\ndZLur7gcAEgSPeAZoK8rDtnGIt84FWa7sf8714vwcN7GIds4ZJsWZsABAGXKuv0EAMpAD3gG6OuK\nQ7axyDdOhdlm/wVMzts4ZBuHbNPCDDgAoEzMgAPAGOgBzwB9XXHINhb5xkmhB7yi9w/HeRuHbOOQ\nbVqYAQcAlCn7FhQA6BQ94BmgrysO2cYi3zgVZpt9CwrnbRyyjUO2aWEGHABQip5a3TQNWlAAoFP0\ngGeAvq44ZBuLfONUlO22kuZIeljSgxW8f1dw3sYh2zhkmxZmwAEAZdnYfpLrRXgAoAz0gGeAvq44\nZBuLfONUlO20+AIm520cso1DtmlhBhwAUBb6vwFgHOgBzwB9XXHINhb5xqko2+xXQJE4byORbRyy\nTQsz4ACAskyLFhQA6NSsbr1RX1+fGo1Gt95uWjGzBXyyjUG2sVLLt6dWf6akMyRtV3Utnfrthoee\n9sIZWz/Q5bfdt7jNegY8tfM2J2Qbh2zT0rUBOABMAU1Jh1ddRBm2qO4PnH+V9Puq3hwApoKuDcDp\nAY/DJ9o4ZBsrwXz/prj9rKQLqiykU/vOmFPVW69tDvTfWtWbd0OC5202yDYO2aaFGXAAeNLzitsL\nmwP9V1RaCQAgW6wDngHW9oxDtrFSyrenVp+hJ3uY/1BlLWVIKdvckG0cso1DtmlhFRQAaJmr1mXU\n72oO9N9TdTEAgHyxDngG6OuKQ7axEst3sP1kys9+S8llmxWyjUO2ccg2LcyAA0BLVgNwAEC66AHP\nAH1dccg2VmL5ZjUATyzbrJBtHLKNQ7ZpGXMAbmZzzewSM/uDmf3ezN49wn6nmdn1ZrbKzA4ov1QA\nCJXVABwAkK7xLEP4V0nvc/c+M9tK0m/N7Gfu3j+4g5ktkrSnu+9lZgdJOl3Swe0vQg94HPq64pBt\nrFTyzW0FFCmdbHNEtnHINg7ZpmXMGXB3v8Pd+4r7D0vql/TMIbsdKemsYp8rJW1rZjuVXCsARHmW\npC0l3dEc6L+36mIAAHmbUA+4mc2TdICkK4c8taukm9u2b5G0W/sO9IDHoa8rDtnGSijf7NpPEso2\nO2Qbh2zjkG1axj0AL9pPzpP0nmImfJNdhmx7+8aKFStkZsvMbEnx8972k8HMFrDNNttsV7H90yfu\nO2L1hkekYgBedT1lbEvav5Pj2R55W9L+KdXDNttsd3+7+FlirbHtsolONJu7j72T2WaSfiTpInc/\nZZjn/0vSpe7+7WL7WkmHuPudg/v09vZ6o9EYOkgHgMr11OpnS3qzpOOaA/1fqboeAMDUMtFx7nhW\nQTFJZ0paPdzgu/BDSccU+x8s6f72wTcAJC67FhQAQLrG04LyUklHS3qlmV1d/LzGzI4zs+Mkyd0v\nlHSDma3V+dIAAAAgAElEQVSRdIakfx76IvSAx2n/0wjKRbaxUsi3WAGlXmxmMwBPIdtckW0cso1D\ntmkZcxlCd79C41st5fhSKgKA7nq2pC0k3dYc6L+/6mIAAPnr2pUwWQc8Dmt7xiHbWInkm2X7SSLZ\nZols45BtHLJNS9cG4ACQqCwH4ACAdHVtAE4PeBz6uuKQbaxE8s1yAJ5Itlki2zhkG4ds08IMOIDp\nLssBOAAgXfSAZ4C+rjhkG6vqfHtq9ZmS9ik2V1dZS9mqzjZnZBuHbOOQbVqYAQcwnT1H0uaSbmkO\n9D9QdTEAgOmBHvAM0NcVh2xjJZBvtu0nCWSbLbKNQ7ZxyDYtzIADmM6yHYADANJFD3gG6OuKQ7ax\nEsg32wF4Atlmi2zjkG0csk0LM+AAprNsB+AAgHSNeSn6svT19anRaHTr7aYVM1vAJ9sYZBunp1b/\n0KoNDx+/34ytBios41nFbVYroEicu5HINg7ZxiHbtHRtAA4Ag4rl/z64mWybqmuR9MvmQP9DVRcB\nAJg+ujYApwc8Dp9o45BtmH0lbbPvjDk3S1pQcS1rK37/EJy7ccg2DtnGIdu0MAMOoAovLm5/0Rzo\nv6HSSgAA6DLWAc8Aa3vGIdswL5aki564966qC8kV524cso1DtnHINi2sggKgCi+WpLU+wOojAIBp\nh3XAM0BfVxyyLV9Prf4MSc+V9Nhxs3b5etX15IpzNw7ZxiHbOGSbFmbAAXTbwcXtyuZA/+OVVgIA\nQAXoAc8AfV1xyDbE4Bcwf0m+ccg2DtnGIds4ZJsWZsABdNvGAXilVQAAUBHWAc8AfV1xyLZcPbX6\nLEkHFZu/dPc7q6wnZ5y7ccg2DtnGIdu0MAMOoJv+RtIcSTc0B/oZfAMApiV6wDNAX1ccsi3dU9pP\nyDcO2cYh2zhkG4ds08IMOIBuov8bADDtsQ54BujrikO2pdt4CXqJfCORbRyyjUO2ccg2LcyAA+iK\nnlp9R0l7SnpE0u8qLgcAgMrQA54B+rrikG2p2i/As14i30hkG4ds45BtHLJNCzPgALqF/m8AAEQP\neBbo64pDtqV6Sv+3RL6RyDYO2cYh2zhkmxZmwAGE66nVN5N0YLH5qyprAQCgavSAZ4C+rjhkW5oX\nSNpC0vXNgf57Bh8k3zhkG4ds45BtHLJNS9cuRQ+gMz21+jxJF0vaoeJSJmOz4vYXo+4FAMA0MOYA\n3My+Jum1ku5y9+cP8/wCST+QdEPx0Pnu/smh+9EDHoe+rjiJZbtI0h5VF9GB9ZK+2/5AYvlmhWzj\nkG0cso1DtmkZzwz41yV9QdLZo+yzwt2PLKckACN4bnH7EUmnVVnIJD3eHOh/tOoiAACo2pgDcHe/\n3MzmjbGbjfU6fX19ajQa4ywLE2FmC/hkGyOxbAcH4Nc0B/ofqLSSkiSWb1bINg7ZxiHbOGSbljK+\nhOmSXmJmq8zsQjPbt4TXBLCpwQH4HyutAgAAdKSML2FeJWmuu68zs9dI+r6kvYfuRA94HD7Rxkkl\n255afQtJz1Krj/pPFZdTmlTyzRHZxiHbOGQbh2zT0vEMuLs/5O7rivsXSdrMzJ4+dL/zzjtPZrbM\nzJYUP+9tXxLHzBawzTbbw29/cf1tzdUbHjFJNzQH+v9adT1ss80222yzPZ23i58l1hrbLpvoctvm\n7mPv1OoBv2CEVVB2UmuFFDezAyWd6+7zhu538skn+wknnDBmrzgmzoy+riipZNtTq79e0rmSLmgO\n9GfzhedU8s0R2cYh2zhkG4dsY/X29nqj0Rj3OHc8yxD2SDpE0vZmdrOkj6pY09fdz5B0lKR3mtl6\nSeskvXEyhQMYFf3fAABkYjyroDTHeP5Lkr401uvQAx6HT7RxEso2ywF4Qvlmh2zjkG0cso1Dtmnp\n2qXoAXQkywE4AADTUdcG4BNtTsf4tX85AOVKIdueWt2U6QA8hXxzRbZxyDYO2cYh27QwAw6kbydJ\n20i6X9LdFdcCAAA61LUBOD3gcejripNIthtnv5sD/WMvWzSFJJJvlsg2DtnGIds4ZJsWZsCB9GXZ\nfgIAwHRFD3gG6OuKk0i22Q7AE8k3S2Qbh2zjkG0csk0LM+BA+rIdgAMAMB3RA54B+rriJJJttgPw\nRPLNEtnGIds4ZBuHbNPCDDiQsJ5afbakZ0tySWsqLgcAAJSAHvAM0NcVJ4Fs95A0U9KfmwP9j1Vc\nS+kSyDdbZBuHbOOQbRyyTQsz4EDasm0/AQBguqIHPAP0dcVJINt9itssB+AJ5Jstso1DtnHINg7Z\npoUZcCBtzIADAJAZesAzQF9XnASyzXoAnkC+2SLbOGQbh2zjkG1amAEH0pb1ABwAgOmIHvAM0NcV\np8pse2r17SU9XdLDkm6vqo5InLtxyDYO2cYh2zhkmxZmwIF0bZz9bg70e6WVAACA0tADngH6uuJU\nnG327Secu3HINg7ZxiHbOGSbFmbAgXRlPwAHAGA6ogc8A/R1xak42+wH4Jy7ccg2DtnGIds4ZJuW\nWVUXAOSup1bfVdJekzj0BcVttgNwAACmo64NwPv6+tRoNLr1dtOKmS3gk22MTrPtqdXnSPqdpO0m\n+RIu6frJvn/qOHfjkG0cso1DtnHINi3MgAOx/latwff9klZN4vj/1xzof6TckgAAQJW6NgCnBzwO\nn2jjlJDt/OL23OZA/3EdvlZ2OHfjkG0cso1DtnHINi2sggLEOrC4/XWlVQAAgGSwDngGWNszTgnZ\nDs6Ar+zwdbLEuRuHbOOQbRyyjUO2aWEGHAhSXEr+OZLWSVpdcTkAACARrAOeAfq64nSY7YuK26ua\nA/3rSygnO5y7ccg2DtnGIds4ZJsWZsCBOPR/AwCATdADngH6uuJ0mC3932Pg3I1DtnHINg7ZxiHb\ntDADDgToqdVNzIADAIBh0AOeAfq64nSQ7VxJO0r6i6Q/l1ZQZjh345BtHLKNQ7ZxyDYtzIADMTa2\nnzQH+r3SSgAAQFLGHICb2dfM7E4z+90o+5xmZteb2SozO2C4fegBj0NfV5wOsh1sP6H/exScu3HI\nNg7ZxiHbOGSblvHMgH9d0mEjPWlmiyTt6e57STpW0ukl1QZMZXwBEwAADGvMAbi7Xy7pvlF2OVLS\nWcW+V0ra1sx2GroTPeBx6OuKM5lse2r1GXpyDXAG4KPg3I1DtnHINg7ZxiHbtJTRA76rpJvbtm+R\ntFsJrwtMVc+VtLWkm5sD/XdUXQwAAEhLWV/CtCHbm3zpjB7wOPR1xZlktiw/OE6cu3HINg7ZxiHb\nOGSbljIG4LeqteTaoN2Kx55ixYoVMrNlZrak+Hlv+8lgZgvYZjuT7fmrNzyinz5x318SqYftabgt\naf9Ojmd75G1J+6dUD9tss9397eJnibXGtssmOtFs7mOvkGZm8yRd4O7PH+a5RZKOd/dFZnawpFPc\n/eCh+/X29nqj0Rg6Uw5kp6dW/7VaX8Jc2Bzov6TqegAAQKyJjnNnjbWDmfVIOkTS9mZ2s6SPStpM\nktz9DHe/0MwWmdkaSY9IetvkSgemvp5avSZpP7XasH5bcTkAACBB41kFpenuz3T32e4+192/Vgy8\nz2jb53h339Pd93P3q4Z7HXrA47T/aQTlmkS2L5A0W9K1zYH+B8uvKC+cu3HINg7ZxiHbOGSbFq6E\nCZSL9b8BAMCoxmxBKQvrgMfJbW3Pnlr9JZJOk7RF1bV8a/Y+6qnVJ3LIzsUtK6CMQ27nbkrINg7Z\nxiHbOGSblq4NwIEJ+CdJL6y6iA6sl3Rx1UUAAIA0dW0A3tfXp0aj0a23m1bMbEFmn2z3Km7fropn\nkk9df+v898zadaLtJHc3B/rvCikoMxmeu8kg2zhkG4ds45BtWpgBR4oGB+C9zYH+tVUW8iazHZoD\n/X+osgYAAJCXrn0Jkx7wODl9ou2p1bdWq4/6MUm3VFxOVtmmiHzjkG0cso1DtnHINi2sgoLU7Fnc\n/qk50L+h0koAAAACdG0AzjrgcTJb23Ow/eT6SqsoZJZtcsg3DtnGIds4ZBuHbNPCDDhSs3dxm8QA\nHAAAoGz0gGcgs76upGbAM8s2OeQbh2zjkG0cso1DtmlhBhypSWoADgAAUDZ6wDOQWV9XUgPwzLJN\nDvnGIds4ZBuHbOOQbVqYAUcyemr1bSVtL2mdpNsqLgcAACAEPeAZyKiva+Psd3Og3yutpJBRtkki\n3zhkG4ds45BtHLJNCzPgSAkroAAAgOx17VL0fX19ajQa3Xq7acXMFmTyyTap/m8pq2yTRL5xcst2\n/tLlu0s6Q9I2Vddy/+pfbLPtvi95sOo6ckS2ccg21mf+dmL7d20ADoxDcgNwAMl4j6TDqi5CkmbU\ntqy6hGyRbRyyTUvXBuD0gMfJaJYruQF4RtkmiXzj5JTt/KXLZ0h6fbH5NlX834ht9uDfsyhkG4ds\nw10xkZ2ZAUcSemp1U4IDcABJOEjSXEk3Szp75eKFGyquBwCeore3d0L7sw54BjJZ2/MZkraV9KCk\nuyquZaNMsk0W+cbJLNs3FLffTWHwnVm2SSHbOGSbFlZBQSo2roCSyhKEAKo3pP3ku1XWAgBlYR3w\nDGTS65lk+0km2SaLfONklO2LJe0qaa2kKyuuRVJW2SaHbOOQbVqYAUcqkhyAA6hce/sJfx0DkAV6\nwDOQSV9XkgPwTLJNFvnGySHbIe0n51ZZS7scsk0V2cYh27QwA45UJDkAB1Cpl0raRdJNklZWXAsA\nlIYe8AxM9b6ulJcgnOrZpo5842SS7cbZ75TaTzLJNklkG4ds08IMOFKws6StJN3bHOj/S9XFAKje\n/KXLZ0o6qthMpv0EAMpAD3gGMujrSnL2W8oi26SRb5wMsh1sP/mzpN9WXMtTZJBtssg2DtmmhRlw\npCDZATiAygyufpJU+wkAlKFrl6KnBzxOBn1dyQ7AM8g2aeQbJ4Vs5y9dPlvSOZLqkzh8z+I2ufaT\nFLLNFdnGIdu0dG0ADowi2QE4gI4cqSdnsidjlaSrS6oFAJJBD3gGMujrSnYAnkG2SSPfOIlk+z+L\n249L2m8SPy9Osf0kkWyzRLZxyDYtzICjUj21+gwlPAAHMDnzly6fJ+lVkgYknbpy8cJ7q60IANIx\nrgG4mR0m6RRJMyV91d3/c8jzCyT9QNINxUPnu/sn2/ehBzyOu1/aU6sfI+nAqmuZhM2Ln7ubA/0P\nVF3MUPTMxSLfOAlk+zZJJun83AbfCWSbLbKNQ7ZpGXMAbmYzJX1R0qGSbpW00sx+6O79Q3Zd4e5H\nBtSIMfTU6ntLOqvqOjq0quoCAJSjWMP77cXmV6usBQBSNJ4Z8AMlrXH3GyXJzL4t6XWShg7AbbQX\n6evrU6PRmEyNGMMlT9z/kVfO3FaSeiV9v+JyJmODpB9VXcRwzGwBswZxyDdOxdm+StJuktZIqqqG\nMJy3ccg2DtmmZTwD8F0l3dy2fYukg4bs45JeYmar1JolP9HdV5dTIkbTU6vPfIZt9upi8+PNgf7L\nKi0IAKR3FLdnpvglSgCo2ngG4OP5j+dVkua6+zoze41as7B7t+9AD3iYQ18wY872avXfX151Mblh\ntiAW+capKtv5S5fvLOkISU9o6rfGDYvzNg7ZxiHbtIxnGcJbJc1t256r1iz4Ru7+kLuvK+5fJGkz\nM3t6+z7nnXeezGyZmS0pft7bviSOmS1ge1Lbb5Wki5+479I3PX7tIQnUwzbbbE/v7WMe/FPfrPuu\nuexXKxcvvD2Bethmm222S98ufpZYa2y7bKLLbZv76BPcZjZL0h8lNSTdJunXkprtX8I0s50k3eXu\nbmYHSjrX3ee1v87JJ5/sJ5xwwqh94piYnlp9O0m3r97wyOx9Z8yZ1xzoX1t1Tbkxo2cuEvnGqSLb\n+UuXm1r/Xuwl6fCVixf+uJvv3y2ct3HINg7Zxurt7fVGozHuce6YLSjuvt7Mjpf0E7WWITzT3fvN\n7Lji+TMkHSXpnWa2XtI6SW+cVPWYqH+UVFunDVcx+AaQgFeoNfi+Va1/MwAAwxhzBrwsE/1kgLH1\n1OpXqrVKzdHNgf5vVl0PUIX5S5dvIek8SftUXQu0XfHzyZWLF/7vqosBgG4pfQYcaeqp1fdVa/D9\noKTvVVwOUKUPS1pUdRHY6BGx9jcAjKprA3DWAS/dW4vb77zp8WsPbGa41m4K6JmL1Wm+85cu30fS\nvxabR2jT6xNMW9d/7UMH7fX2T19ZwVvfs3LxwuSualsm/rsQh2zjkG1amAGfgnpq9VmS3lxsLpM0\nu7pqgGoUX/g7XdJmkr66cvHCJC/mVBX7wK/nrly88E9V1wEA2FTXBuCsA76pYiC9pzT6VUSHcbCk\nndVabeCX3q1G/mmI2YJYHeb7PyQtkHSPpA+WUU9OOHfjkG0cso1DtmlhBrxa35H0Dx0cv6w50M/g\nG9PO/KXLt5P0uWJz8crFC/9SZT0AAEzEeC7EU4qJLlCeu55afR+1Bt9/lXTtJH4uVfFFp/YF4lEu\nso3VQb6fkrSDpMuU6dUWO8W5G4ds45BtHLJNCzPg1fmX4nZZc6D/2EorwbQzf+ny2ZIOUsX/Ddjx\npX+/X9HLPRG7SjpO0npJ/7xy8UL+CgQAmFLoAa9AT62+raS3FJundfp69HXFyTjbr6nVQ12p3V93\nfCeHn7Ry8cI/lFVLbjI+dytHtnHINg7ZpoUZ8Gq8XdIcSb3Ngf7fV10Mppf5S5f/d7UG349KqmKZ\nujLcLOkTVRcBAMBksA54l/XU6jP1ZPvJqWW8Jmt7xskt2/lLl++g1tJ9UuvLi1+qsp7c8k0J2cYh\n2zhkG4ds09K1L2FioyMkzZN0g6QLqy0F09CX1Pry4iV6ciAOAAC6qGsDcHrAN3pPcfuF5kD/E2W8\nIJ9o4+SU7fyly98g6fWSHpb09pWLF26ouKSs8k0N2cYh2zhkG4ds08IMeBf11OovUOvCIQ9L+nq1\n1WA6mb90+U6Svlxsnrhy8cIbKywHAIBpjR7w7np3cfv15kD/A2W9KH1dcdqznb90+VxJz6u2okk7\nXtIzJF0s6SsV17IR524cso1DtnHINg7ZpoVVUErQU6vvKGmhRr+k/Cw9uezbF8KLQqnmL12+QNJP\nJW1WcSmdeEjSP7FuNgAA1WId8A711Op7S7pc0o7jPOTC5kD/9WXWwCfaOO5+6fyly/eUdL5ag++V\nku6ttqpJ2SDpyysXL1xbdSHtOHfjkG0cso1DtnHINi3MgHegp1Z/llp/0t9R0ipJ/WMcMiDp09F1\noTzzly7fVtIFkp4u6UeS/m7l4oWlfHkWAABMT/SAT1JPrb6LpF5JcyX9XNKrmwP9j1RRC31dMeYv\nXT7rvt9f8bPt/uZl+0j6naQ3MfguF+duHLKNQ7ZxyDYO2aaFVVAmoadW316tme89JF0l6bVVDb4R\n6tSZW2z1Ikl3STpi5eKFD1VdEAAAmPqmbQ94T62+jaQ3SqpN4vC3StpX0mq1Zr5LW9FkMlL9RDt/\n6fKDJZ0gafOqa5mELSUt3GaP/R9Xq+3kpqoLylGq524OyDYO2cYh2zhkm5bp3AN+mqS3dHD8nyQd\n2hzov6ekerJSDL4vljSn6lo69E8rFy/8ZdVFAACAfEzLHvCeWr0u6c2S1qu1JvJErwj4kKQvNQf6\nby+7tslIra9r/tLl+0u6SK3B97clfavaiibtT7/5QGNHLWbVviipnbs5Ids4ZBuHbOOQbVqm6wz4\nx9Tqf/9Kc6D/XVUXk5P5S5fX1Vove1tJ35P05pWLF66vtqrJsw+Me3lJAACAcenalzBT6QHvqdX3\nl/R6tZYE/I+KyylFKp9o5y9d/hy12k52kPQTSc2pPPiW0sk2V+Qbh2zjkG0cso1DtmmZjjPgHy9u\nv9wc6L+l0koSNH/p8u0lfV3S9pM4fA+1Bt+XSfqHlYsXDpRZGwAAQA6mVQ94T61+sKQjJK2T9JlK\niylRyX1d/yzp8A6O/7VaS/atK6meStEzF4t845BtHLKNQ7ZxyDYt020G/BPF7anNgf67Kq0kQfOX\nLjc9uTLM/5J0zQRfYr2kq6d62wkAAECkabMOeE+tvkDSoZIekLS00mJKVuIn2pdJeo6kWyV9las+\n0jMXjXzjkG0cso1DtnHINi1Tbga8p1Z/taSdJ3Ho8cXtyc2B/vtKLCkng7Pf5zD4BgAAiDGlesB7\navV/VGtd6cn6i6RTOyoiQWX0dc1funxLSW8oNs/quKhM0DMXi3zjkG0cso1DtnHINi1TZga8p1bf\nQk+2jlwoaaJXoHRJ32gO9D9YamH5+HtJW0u6cuXihddWXQwAAECuplIP+ImS5krqk3Rkc6CfFolC\nSZ9oB9tPmP1uw2xBLPKNQ7ZxyDYO2cYh27R07UI8neip1XeV9MFi870Mvss1f+ny3dT6gurjkr5T\ncTkAAABZG3MAbmaHmdm1Zna9mf3rCPucVjy/yswOGG6fvr6+Tur8jKQtJZ3fHOhf0ckL5cjMFnT4\nEkdLMkk/XLl44b2dV5SPErLFKMg3DtnGIds4ZBuHbNMy6gDczGZK+qKkwyTtK6lpZvUh+yyStKe7\n7yXpWEmnD/daa9asmVSBPbX6QWoNEAckLZ7Ui+Rv0v09xdrfby02l5VRTGaqXT8zf+Qbh2zjkG0c\nso1DtoEmOtE81gz4gZLWuPuN7v5XtVYged2QfY5U0Tfs7ldK2tbMdhr6Qo888siECpOknlrdJJ1S\nbH6uOdD/5wm/yPSwbQfHHijpuZLulPSTcsrJSifZYmzkG4ds45BtHLKNQ7aBVq1aNaH9x/oS5q6S\nbm7bvkXSQePYZze1BnRPceIr3/U/JlLcrs/c/YCdb1t78BMzZt7znWPff9Hnli5/4USOny423/FZ\nu8yffDb/Utx+kytYAgAAxBtrAO7jfB0b67g77rhDr/nFym+M8/We4uLXNbe/Y7d5l03m2Olgy7nP\nlVrtP51g9ZPhzau6gMzNq7qAjM2ruoCMzau6gIzNq7qAjM2rugA8ydxHHmOb2cGSlrj7YcX2hyRt\ncPf/bNvnvyRd6u7fLravlXSIuz9lBvyd73ynt7eh7LfffpVfnj4XfX19ZBmEbGORbxyyjUO2ccg2\nDtmWq6+v7yltJ3PmzNHpp58+dEJ6RGMNwGdJ+qOkhqTbJP1aUtPd+9v2WSTpeHdfVAzYT3H3gyf8\nmwAAAADTwKgtKO6+3syOV+vLeTMlnenu/WZ2XPH8Ge5+oZktMrM1kh6R9LbwqgEAAIApatQZcAAA\nAADl6sqVMMdzMR9MjpndaGbXmNnVZvbrquuZyszsa2Z2p5n9ru2xp5vZz8zsOjP7qZmxjNMkjJDt\nEjO7pTh3rzazw6qscaoys7lmdomZ/cHMfm9m7y4e59zt0CjZcu6WwMw2N7MrzazPzFab2aeLxzl3\nOzRKtpy7JTGzmUWGFxTbEzpvw2fAi4v5/FGtS53fKmmlhvSRY/LM7M+SXujuXMGyQ2b2ckkPSzrb\n3Z9fPPZZSfe4+2eLD4/bufsHq6xzKhoh249KesjdP1dpcVOcme0saWd37zOzrST9VtLfqdUOyLnb\ngVGyfYM4d0thZlu6+7riO2dXSDpRreuLcO52aIRsG+LcLYWZvV/SCyVt7e5HTnS80I0Z8PFczAed\nGfe3bjEyd79c0n1DHt54oani9u+6WlQmRshW4tztmLvf4e59xf2HJfWrdX0Gzt0OjZKtxLlbCndf\nV9ydrdZ3ze4T524pRshW4tztmJntJmmRpK/qyTwndN52YwA+3IV6dh1hX0ycS7rYzH5jZu+oupgM\n7dS2pOadkja5yis68i9mtsrMzuTPzJ0zs3mSDpB0pTh3S9WW7a+Khzh3S2BmM8ysT61z9BJ3/4M4\nd0sxQrYS524ZPi9psaQNbY9N6LztxgCcb3nGeqm7HyDpNZLeVfypHwG81a/F+Vye0yU9W9L+km6X\ndHK15UxtRYvE+ZLe4+4PtT/HuduZItvz1Mr2YXHulsbdN7j7/mpdQfsVZvbKIc9z7k7SMNkuEOdu\nx8zscEl3ufvVGuGvCeM5b7sxAL9V0ty27blqzYKjBO5+e3F7t6TvqdXyg/LcWfSBysx2kXRXxfVk\nw93v8oJaf8bj3J0kM9tMrcH3Oe7+/eJhzt0StGX7jcFsOXfL5+4PSPqxWj21nLslasv2RZy7pXiJ\npCOL7+D1SFpoZudogudtNwbgv5G0l5nNM7PZkv5R0g+78L7ZM7MtzWzr4v4cSa+S9LvRj8IE/VDS\nW4r7b5H0/VH2xQQU/4Ea9Pfi3J0UMzNJZ0pa7e6ntD3FuduhkbLl3C2HmW0/2AJhZltI+m+Srhbn\nbsdGynZwgFjg3J0Ed/+wu89192dLeqOk5e7+Zk3wvO3KOuBm9hpJp+jJi/l8OvxNpwEze7Zas95S\n66JK3yTbyTOzHkmHSNperf6tj0j6gaRzJe0u6UZJb3D3+6uqcaoaJtuPSlqg1p9BXdKfJR3X1j+H\ncTKzl0m6TNI1evJPnh9S68rFnLsdGCHbD0tqinO3Y2b2fLW+rDaj+DnH3Zea2dPFuduRUbI9W5y7\npTGzQySdUKyCMqHzlgvxAAAAAF3UlQvxAAAAAGhhAA4AAAB0EQNwAAAAoIsYgAMAAABdxAAcAAAA\n6CIG4AAAAEAXMQAHgASY2e5m9lBx8ZeyX3uZmX2iuP9yM7u27PcAAIwfA3AAqIiZ3WhmCyXJ3de6\n+9Yec3EGL37k7pe7+z4B7wEAGCcG4ABQHZdU+oz3CLr1PgCAMTAAB4AKmNk5al2y+IKi9WSxmW0w\nsxnF85ea2SfM7OfF8z80s+3N7Jtm9oCZ/drMntX2evuY2c/M7C9mdq2ZvX6E911gZje3bd9oZieY\n2Sozu9/Mvm1mtbbnDzezPjO7r6jl+XGpAMD0wAAcACrg7m+WtFbS4e6+taTvDrPbP0o6WtKukvaQ\n9EtJZ0p6uqR+SR+VJDObI+lnkr4haQdJb5T0ZTOrj6cUSa+X9GpJz5b0AklvLV73gOL93lG85xmS\nfo/q+OQAAAGKSURBVGhmsyf8CwMANmIADgBpcklfd/c/u/uDki6SdJ27L3f3J9QasB9Q7Hu4pD+7\n+1nuvsHd+yT9X7UG1uNxmrvf4e73SbpA0v7F48dKOsPdV3rL2ZIGJB1czq8IANPTrKoLAACM6M62\n+49JumvI9lbF/WdJOsjM7mt7fpaks8f5Pne03X9U0jPbXvcYM/uXtuc3k7TLOF8XADAMBuAAUJ2J\nrHgy2r5rJa1w91eV/F5rJf2Hu39qAscCAMZACwoAVOdOtXq7R2Ij3B/qx5L2NrOjzWyz4me+me3T\nduxEVkEZ3Pf/SPpfZnagtcwxs9ea2VajHQwAGB0DcACozqcl/buZ3Svpv2vTWWofcn/Y5939IUmv\nUuvLl7dKur147dkjHDvabHj7muG/VesLmF+UdK+k6yUdM47fCwAwCou55gMAAACA4TADDgAAAHQR\nA3AAAACgixiAAwAAAF3EABwAAADoIgbgAAAAQBcxAAcAAAC6iAE4AAAA0EUMwAEAAIAuYgAOAAAA\ndNH/B3q66ha7Vyi1AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "naf.fit(T[treatment],timeline=t, event_observed=E[treatment], label=\"With treatment\")\n", "ax = naf.cumulative_hazard_.plot()\n", "naf.fit(T[~treatment], timeline=t, event_observed=E[~treatment], label=\"Without treatment\")\n", "ax = naf.cumulative_hazard_.plot(ax=ax)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Try a more complicated model?" ] }, { "cell_type": "code", "execution_count": 26, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from lifelines import AalenAdditiveFitter" ] }, { "cell_type": "code", "execution_count": 34, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " [-----------------100%-----------------] 30 of 30 complete in 0.0 sec\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 34, "metadata": {}, "output_type": "execute_result" } ], "source": [ "aaf = AalenAdditiveFitter()\n", "aaf.fit(data, duration_col=1, event_col=2, timeline=t)" ] }, { "cell_type": "code", "execution_count": 35, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 35, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAt4AAAFwCAYAAACPanxcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYbFV97//3Ogco5CigMYoB9BhwKBLCES3iGA7dBtGf\n4uXGyG0TNGpM1MRExf45JdHrEEwKEpM4xCQaxKFNHBPjBJ4SnBItlSYojXJUjKBAFJmlGc66f1Q1\ndOBMXbW/vVftfr+ep5/qXV177y+fZ3N61epvrZ1yzkiSJEmKta7uAiRJkqS1wIG3JEmStAoceEuS\nJEmrwIG3JEmStAoceEuSJEmrwIG3JEmStAoqGXinlNanlM5NKX20iuNJkiRJTVPVjPcfAhcALgou\nSZIkbcfYA++U0kHAE4B/ANLYFUmSJEkNVMWM918Cs8C2Co4lSZIkNdJYA++U0hOBK3LO5+JstyRJ\nkrRDKefR27JTSn8KnATcAuwN7At8MOf89KXXHH/88fnGG2/kgAMOAGDDhg0ceuihbNq0CYD5+XkA\nt0fYXvq+lHqatL30XCn1NGl769atPOUpTymmnqZtm2/c9gc+8AF/fwVt+/vM32eTsL1161auv/56\nAC677DIOOeQQ3vrWt65o4nmsgff/OFBKRwMvyTk/afnzT3/60/MZZ5zhbHiAlNKrc86vrruOJjLb\nOGYby3zjmG0cs41jtnFGGeNWvY73nUbxl112WcWn0DIb6y6gwTbWXUCDbay7gIbbWHcBDbax7gIa\nbGPdBTTYxroL0O32qOpAOedzgHOqOp4kSZLUJOF3rnzc4x4XfYq17PS6C2iw0+suoMFOr7uAhju9\n7gIa7PS6C2iw0+suoMFOr7uApjriiCNWvE/4wHupKV3VyzmfXXcNTWW2ccw2lvnGMds4ZhvHbOOM\nMsatrNVkR+bn55menr7T81u2bFkHvIjBSije8XIE3/rWtzY+8IEPvHgXL0vANcBfTk9Pu9b6bkop\nbfYfqxhmG8t845htHLONY7ZlCR9478SLgE9PT0+fV2MNE+2xj33sbv3PtGXLliMY5H1aeFGSJEna\nrjpbTfZ10D2e3X0HO8x539hqmsXZgThmG8t845htHLONY7ZlCR9474TtJavLvCVJkmoUPvBefuck\nVSultLnuGprKbOOYbSzzjWO2ccw2jtmWpc4Zb0mSJGnNcDnB7Ugp3SOl9OGU0nUppYtTSjN117Q9\n9m3FMds4ZhvLfOOYbRyzjWO2ZalzVZOSvRm4EbgX8BDgYyml83LOF9RbliRJkiaVPd53kFLaAPxv\n4I9zzjfknL8A/AtwUr2V3Zl9W3HMNo7ZxjLfOGYbx2zjmG1Z7PG+swcCt+Scty577jzgF2qqR5Ik\nSQ0Q3moySo93p9urbOm7/uxUWuEud2Vwp8flrgXuVk1F1bFvK47ZxjHbWOYbx2zjmG0csy2LM953\ndh13vtnMfgwG35IkSdJIwme85+fnmZ6eXtE+I8xSV+lbwB4ppUOXtZscAXy9xpq2K6W0W7eM18qZ\nbRyzjWW+ccw2jtnGMduyOON9Bznn64EPAa9JKe2TUno08CTgXfVWJkmSpEnmOt7b93zgLsAVwLuB\n5+acF+ot6c58BxvHbOOYbSzzjWO2ccw2jtmWxXW8tyPn/BPghLrrkCRJUnO4jvcEc23OOGYbx2xj\nmW8cs41jtnHMtiz2eEuSJEmroMh1vLV77NuKY7ZxzDZWE/Oda7XvAuxZdx3v3evBX5trte+43Kwq\nYLZxzDbOvT7+phXvY4+3JKlYc632CcA/4+8rSQ1Q5Dre2j2uzRnHbOOYbawG5vtEBr+rbgRurrOQ\nr2+7fv0vrttwa501NJXZxjHbUCu+q7kzCJKkkj14+PiEmcWFz9RZSAPf1BTDbOOYbZwtW7bkle7j\nOt4TzP+R4phtHLON1cB8lwbe36y1ChqZbTHMNo7ZlsVVTSRJRZprte8J3AO4FvhhzeVI0thcx/sO\nUkoXp5RWrSk9pbQtpfTzw+/fmlL6oxXsuzmssDXObOOYbayG5bs0233hzOLCiv+kW7WGZVsUs41j\ntmWxx/vO8vBr9U+c8/PqOK8kFeq2gXetVUhSRezxnmD2bcUx2zhmG6th+RY18G5YtkUx2zhmWxZ7\nvLfvqJTSN1JKV6aU3pFSaqWU7p5S+reU0hXD5z+aUjpwaYeU0m+llL6dUrompfSdlNLTlv3sWSml\nC4b7fTKldN/tnTSldHpK6bXD7zenlC5JKb04pXR5SukHKaXfWvbaVkrp1JTS91JKlw3bVPYOzESS\nVltRA29JGpc93neWgKcBxwKHAA8E/mj4/NuB+w6/fgq8CSCltAH4K+C4nPO+wCOA+eHPngy8HDgB\nuCfwOWBuB+e+Y5vLvYF9gZ8Dng28OaW03/BnbwAeDhwBHAocCPzJWP/luo09cXHMNlbD8i1mRRNo\nXLZFMds4ZluWInu851rtynqsZxYX0gp3ycCbcs6XAqSUXg/8Tc75j4EPL70opfSnQG/ZftuAw1NK\nl+ScLwcuHz7/XOCUnPM3h/udArwipXRwzvn72zn/8npvBl6Tc94GfCKldB3woJRSH3gO8Jyc81XL\njvse4BUr/O+VpOLMtdot4P4M/m3dWnM5klQJe7y3b/mA+L+An0sp3SWl9LbhqidXA+cA+6WUUs75\neuBEBoPsHwxbUh403P9+wF+llH6SUvoJ8OPh8weyaz8eDrqX3ADcFfhZYB/gLcuO+wkGM+qqgD1x\nccw2VoPyPZTB76jvzCwuLNZdDDQq2+KYbRyzLUuRM94jzFJX7b53+P4HwMkM2k6OyjlfkVLaBHyN\nwQx1zjmfCZyZUmoBrwf+HvgVBgP31+acd9Recke7M9v/IwatLoflnF3bVlIT2d8tqXHs8b6zBPxe\nSunAlNI9gFcC7wPuxmCwe/Xw+VfdtkNK90opPXnY630zcD1w6/DHf8ugteSw4Wv3Syn9+k7Ovcs3\nHcNZ8L8H3pNS+tnhcQ9MKR278v9cbY89cXHMNlaD8i1u4N2gbItjtnHMtiyuanJnmUGv9JnAt4GL\ngNcBbwTuwmC2+YsMWjuWZqfXAS8CLmXQSvIY4HkAOeePAH8GvG/YonI+8Lg7nG/593fc3pGXDs/3\nH8PjnsVgRl6SmqC4gbckjSu81WTSerxzzvcffvtnd/jRD4Fj7vDc3w0fLwM27+SY7wbevYOfrV/2\n/TOXfX82/7PlZXlt5JwXgZN2dE6Nx564OGYbq0H5FrWiCTQq2+KYbRyzLYsz3pKkosy12glY+oC6\nM96SGsMe7wlm31Ycs41jtrEaku99GHyu5scziws/qruYJQ3JtkhmG8dsy+KMtySpNPZ3S2ok1/Ge\nYPZtxTHbOGYbqyH5Fjnwbki2RTLbOGZbljpnvOteq3utMW9Jk6LIgbckjavOHu89tmzZ4mBwDLvb\ntzXMucibJZXKnrg4ZhurIfkufbCymBVNoDHZFsls45htWeqc8f448Iwaz7+WPAP4WN1FSNJucsZb\nUiOlnHfnDuWj27JlS56ent7uzPaWLVueBWwEtoUWsbatAy6enp5+R92FSNKuzLXaG4DrGNwFeJ+Z\nxYVbai5JkrZrZ2PcHam1/cDBoCTpDpbuwHuRg25JTeM63hPMvq04ZhvHbGM1IN9i20wakG2xzDaO\n2ZZlrIF3SmnvlNKXUkrzKaULUkqnVFWYJGlNKnbgLUnjGqvVJOd8Y0rpmJzzDSmlPYDPp5QenXP+\n/NJrXMc7jmtzxjHbOGYbqwH5FrmiCTQi22KZbRyzLcvYrSY55xuG3+4FrAeuHPeYkqQ1yxlvSY01\n9sA7pbQupTQPXA58Jud8wfKf2+Mdx76tOGYbx2xjTXK+c632Ogqe8Z7kbEtntnHMtixjr2qSc94G\nbEop7Qd8KqW02T9rSNLaNddq3xP4AnDwCndNwN7AD2cWF66uvDBJqlllywnmnK9OKX0MeBhw9tLz\nW7duJaV0OnDx8KmrgPmlwfnSOzG3V76dcz67pHrcdnt3t5eUUk/TtpfUdf737vXgg4EHXrDtegAO\nW7cBgBVsf7CkPJdPJi2fXKq7niZt+/vM7QnZ3gTsz8DGU089lenpaVZirBvopJTuCdySc74qpXQX\n4FPA/805b1l6zSiLi0uSJtdcq/1m4PnAHwN/scLd88ziwk+rr0qSqlXHDXTuA7wzpbSOQb/4u5YP\numHQ473SdwPaPctnXlQts41jtrEKyffhw8fPzSwu3LDTV06QQrJtJLONY7ZlGXc5wfOBIyuqRZI0\n4eZa7X2AI4BtwFdqLkeSihJ+50rX8Y7jO9g4ZhvHbGMVkO+RDJaWPX9mceH6mmupVAHZNpbZxjHb\nsoQPvCVJa8pSm8l/1FqFJBUofODtOt5x7riCgapjtnHMNlYB+TZ24F1Ato1ltnHMtizOeEuSqrQ0\n8P5SrVVIUoHs8Z5g9m3FMds4ZhurznznWu2DgAOBqynwzpPj8tqNY7ZxzLYsznhLkqryy8PHL80s\nLmyrtRJJKpA93hPMvq04ZhvHbGPVnG9j+7uh9mwbzWzjmG1ZnPGWJFXF/m5J2gl7vCeYfVtxzDaO\n2caqK9+5VntP4GHDzUYOvL1245htHLMtizPekqQqHA7sDVw0s7jw47qLkaQS2eM9wezbimO2ccw2\nVo35Nrq/G7x2I5ltHLMtizPekqQq2N8tSbtgj/cEs28rjtnGMdtYNebb+Blvr904ZhvHbMvijLck\naSxzrfbPAA8AbgT+s+ZyJKlY9nhPMPu24phtHLONVVO+Rw0fvzqzuHBzDedfFV67ccw2jtmWxRlv\nSdK4Gt9mIklV2CP6BPZ4x7FvK47Zxikx27lWe1/gGcCGumsZ13v3ejBzrfbDd/3KSv3a8LHRA+8S\nr92mMNs4ZluW8IG3JE2Ak4E/qbuICZdp+MBbksYVPvCen59neno6+jRrUkpps+9kY5htnEKznRo+\nvge4pM5CxnX2rVfdd/P6/f+rhlOfO7O4MNHZ7Uqh124jmG0csy2LM96S1rS5VntvBh8OzMALZhYX\nflJzSWN5Wkqb33bLD8+uuw5J0p25jvcE8x1sHLONU2C2HWAv4OuTPuiGIvNtDLONY7ZxzLYsrmoi\naa17zPDxc7VWIUlqPNfxnmCuzRnHbOMUmG2jBt4F5tsYZhvHbOOYbVmc8Za0Zs212uuBRw43GzHw\nliSVyx7vCWbfVhyzjVNYtr8E7At8d2Zx4dK6i6lCYfk2itnGMds4ZlsWZ7wlrWWNajORJJXNHu8J\nZt9WHLONU1i2jRt4F5Zvo5htHLONY7ZlccZb0po012onGjjwliSVyx7vCWbfVhyzjVNQtocC9wau\nAL5Vcy2VKSjfxjHbOGYbx2zL4oy3pLVqabb78zOLC7nWSiRJa4I93hPMvq04ZhunoGwb2WZSUL6N\nY7ZxzDaO2ZbFGW9Ja1UjB96SpHLZ4z3B7NuKY7ZxSsh2rtW+D3AIcB1wXs3lVKqEfJvKbOOYbRyz\nLcsedRcgaW0ZriYy6r89t1TUj7002/3FmcWFWyo4niRJu2SP9wSzbyuO2cbpb7v2P4CbRvy6dK7V\nftlcq73fmGU0ts3EazeO2cYx2zhmWxZ7vCWtmrlW+5ANrDtquHnLCr9uBe4DnAL811yr/Ya5VvuA\nEUtp7MBbklSu8FYTe7zj2LcVx2zDHHfYug0A/zyzuHDiSnYctqj8KvAy4BjgpcAL51rtrSPUcRhw\nM/DlEfYtmtduHLONY7ZxzLYs9nhLWk3HDR8/udIdh73dZwJnzrXaRzEYeJ8A/MKItXxyZnHhpyPu\nK0nSitnjPcHs24pjttWba7X3BqYu2HY9jDDwXm5mceHLM4sLvwYcCBw+wtcvAr82Tg2l8tqNY7Zx\nzDaO2ZbFGW9Jq+XRwD6L5G/PLC78sIoDDo9TybEkSYrmOt4TzL6tOGYb4jiAh6y76wfqLqTJvHbj\nmG0cs41jtmVxVRNJq+Xxw8ex2kwkSZpU9nhPMPu24phtteZa7fsyWEnk2t++6Vu2uAXy2o1jtnHM\nNo7ZlsUZb0mr4XHDxy03sM07RUqS1iR7vCeYfVtxzLZyt7WZmG0s841jtnHMNo7ZlsUZb0mh5lrt\nPYHHDjft75YkrVn2eE8w+7bimG2lHgHcDViYWVz4ntnGMt84ZhvHbOOYbVnGGninlA5OKX0mpfSN\nlNLXU0p/UFVhkhrD1UwkSWL8G+jcDLwo5zyfUror8NWU0lk554WlF9jjHce+rTglZTvXaj8IeDuw\n7wi73wJ8CPjrmcWFayotbPct3Sb+E1BWtk1kvnHMNo7ZxjHbsow18M45XwZcNvz+upTSAvBzwMJO\nd5S0Ek8HHjXG/g8BXjTXap8K/M3M4sJ11ZS1a3Ot9n2ATcBPgc+t1nklSSpRZevpppQ2MvgF/6Xl\nz8/PzzM9PV3VabRMSmmz72RjFJbt0p+NXgxsWeG+BwIvBx4D/Cnw4rlW+23AldWVt1O/OHz8zMzi\nwo1QXLaNY75xzDaO2cYx27JUMvAetpl8APjDnPOqzaZJa8QRw8d/m1lcuGiF+/7nXKv9SWAaeC3w\ncOCVVRa3mz5ewzklSSrK2APvlNKewAeBd+ecP3LHn2/dupWU0unAxcOnrgLml959LX3a1u2Vb+ec\nzy6pHrer3964bu8nn7T+Xgcetm7D9cC3xzjep+da7S3vueWK2YPSXg89ev3+lwKcc+tVBwEcvX7/\nS6K2byZf+9j1d3/H8nqW1J1vU7eXlFJPU7aXniulniZt+/vM7QnZ3gTsz8DGU089dcVdHSnnvKId\n/sfOKSXgncCPc84v2t5rtmzZkqenp9PIJ5HWsLlWexr4NPDvM4sLj6y7HkmSNDDKGHfcdbwfBfwm\ncExK6dzh13HLX+A63nHuOLul6hSU7VJ/93m1VlGhgrJtJPONY7ZxzDaO2ZZl3FVNPo93v5QiLfV3\n+w5WkqQJFz5odh3vOMv7DlWtgrJt3Ix3Qdk2kvnGMds4ZhvHbMvibLVUqLlWuwW0gQycX3M5kiRp\nTOEDb3u849i3FaeQbA9j0A520cziwvV1F1OVQrJtLPONY7ZxzDaO2ZbFGW+pXPZ3S5LUIPZ4TzD7\ntuIUkm3j+ruhmGwby3zjmG0cs41jtmVxxlsq19KMd6MG3pIkrVX2eE8w+7bi1J3tXKudaGirSd3Z\nNp35xjHbOGYbx2zL4oy3VKaDgbsDPwZ+UHMtkiSpAvZ4TzD7tuIUkO1ts90ziwu51koqVkC2jWa+\nccw2jtnGMduyOOMtlamRH6yUJGkts8d7gtm3FaeAbBvZ3w1FZNto5hvHbOOYbRyzLYsz3lKZnPGW\nJKlh7PGeYPZtxakz27lW+27AIcBNwIV11RHF6zaW+cYx2zhmG8dsy+KMt1SeXxo+XjCzuHBTrZVI\nkqTK2OM9wezbilNzto3t74bas208841jtnHMNo7ZlsUZb6k89ndLktRA9nhPMPu24tScbaNnvL1u\nY5lvHLONY7ZxzLYse9RdgNREw1u+Pw94wAi7L/V4O+MtSVKDhA+85+fnmZ6ejj7NmpRS2uw72RgV\nZHsk8OYx9t86s7jwkzH2L5bXbSzzjWO2ccw2jtmWxRlvKcZjho+fAz48wv5nVliLJEkqQPjA2x7v\nOL6DjVNBto8aPr5zZnHh7WMeq1G8bmOZbxyzjWO2ccy2LK5qIlVs2N+9NPD+Qp21SJKkcriO9wRz\nbc44Y2Z7P+A+wJXANyspqEG8bmOZbxyzjWO2ccy2LM54S9Vbmu3+4sziQq61EkmSVAzX8Z5g9m3F\nGTNb20x2wus2lvnGMds4ZhvHbMvijLdUvUcOHx14S5Kk29jjPcHs24ozarZzrfa+wOHAzcBXqqyp\nKbxuY5lvHLONY7ZxzLYsznhL1fplBv9ffW1mceGndRcjSZLKYY/3BLNvK84Y2drfvQtet7HMN47Z\nxjHbOGZbFme8pWrdtqJJrVVIkqTi2OM9wezbijNKtnOt9nrg4cNNZ7x3wOs2lvnGMds4ZhvHbMvi\njLdUncOBuwLfmVlcuKzuYiRJUlns8Z5g9m3FGTFb20x2g9dtLPONY7ZxzDaO2ZbFGW+pOq7fLUmS\ndsge7wlm31acEbN1RZPd4HUby3zjmG0cs41jtmXZo+4CpCVzrfaRwAuBPWsu5XuHpL2/tJId5lrt\nA4H7AdcAF4RUJUmSJlr4wNse7zgN7Nt6FXB83UUAvHbPjdfPtdqnAqfNLC5cuxu7LLWZ/PvM4sKt\ngaVNvAZet0Ux3zhmG8ds45htWZzxVkkOGz6+GKhrVZB1wInAkxi8EXj+XKv9GuDjQN7JfscOH20z\nkSRJ2xU+8J6fn2d6ejr6NGtSSmlzU97JzrXaewM/D2wD3jKzuLBYYznvOXb93V/wzD0OmAEeAfzN\n8Gt3uKLJLjTpui2R+cYx2zhmG8dsy+KMt0pxKIPZ5q01D7oBOGvbVec/kwPeBDwZeBlw793Y7ULg\n86GFSZKkiWWP9wRr2DvYBw8fF2qtYmhZth8ZfqkiDbtui2O+ccw2jtnGMduyuI63StEePl5YaxWS\nJElBXMd7gjVsbc6lGe8iBt4Ny7YoZhvLfOOYbRyzjWO2ZXHGW6UoqtVEkiSpauEDb3u84zSlb2uu\n1V7H7QPvb9ZZy5KmZFsis41lvnHMNo7ZxjHbsjjjrRIcBOwDXDGzuHBl3cVIkiRFsMd7gjWob2vp\ng5XFtJk0KNvimG0s841jtnHMNo7ZlmXsgXdK6R0ppctTSudXUZDWpKI+WClJkhShihnvfwSO29EP\n7fGO06C+reI+WNmgbItjtrHMN47ZxjHbOGZblrEH3jnnzwE/qaAWrV2u4S1JkhrPHu8J1qC+reJa\nTRqUbXHMNpb5xjHbOGYbx2zLEj7wPuecc0gpnZ5SevXw64XLL4KU0ma31+72Qan1xAu2XX9v4Abg\n+3XX43b8NrBpnP3dNt+6toFNJdXjtttur/r2C9Pt49nTR5lcTjnnFe90p4OktBH4aM758Dv+bMuW\nLXl6ejqNfRI10lyr/Qjgi8C5M4sLR9ZdjyRJ0u4YZYzrOt6qW3EfrJQkSYpQxXKCcwxmLB+YUvp+\nSumZy39uj3ec5X8CmWDF9XdDY7ItktnGMt84ZhvHbOOYbVn2GPcAOeeZKgrRmuWKJpIkaU0Ye+C9\nK67jHSfnfPZcq30S8MARD/GvM4sL/SprGkGRrSauexrHbGOZbxyzjWO2ccy2LOEDb8WZa7UfDpwx\nxiFm51rtR84sLnytqppWYq7VbgE/D2wDttZRgyRJ0moJH3jPz88zPT0dfZo16XO3Xv3Sx6zfD+As\n4HMr3P0o4InAB+Za7YfOLC7UcROkQ4H1wLdnFhdurOH8O5RS2uwsQQyzjWW+ccw2jtnGMduyOOM9\noeZa7T33T3scM9x86cziwrkr3H9v4AvAkcAZc632k2cWF7ZVXecuFNlmIkmSFCF8OUF7vMMce/i6\nDfsBFwArXjpmOMP8FOAnDGa+X1ptebulyBVNwJ64SGYby3zjmG0cs41jtmVxHe/J9RvDx/fMLC6M\ndBekmcWF7wInDTdfN9dqH7Oz1wdYWtHEGW9JktR49nhPoLlW+67Aky/Ydj2Hrdvw3nGONbO48LG5\nVvv1wCuB98212qcAKx3I3wJ8bGZx4eIV7lfsjLc9cXHMNpb5xjHbOGYbx2zLYo/3ZPpfwD4/ZdvX\nRxjsbs+rgIcD08BfjniMv5xrtf8WeP3M4sLlu3rxXKudKHjgLUmSVDXX8Z5MvwHw0HV3e0sVB5tZ\nXLh1rtX+deD3gXuOcIiDgBOAFwDPmmu1/wI4bWZx4epd7LMBuGJmceHKEc4ZytmBOGYby3zjmG0c\ns41jtmVxxnvCzLXa9wJ+lUF7x/urOu5wOcHXjlHX4cDrgScBfwy8bK7Vvmknu6wfPjrbLUmS1oTw\nD1fOz694wQ3t3IkMBq2feNpNF/5i3cUsmVlcOH9mceF44NHAZ4E9Gcxo7+hrbwa95B+ppeBdSClt\nrruGpjLbWOYbx2zjmG0csy2LM96T57bVTGqtYgdmFhe+ABw912rvw67f2N06s7jw01UoS5IkqXb2\neE+QuVb7UOCXgeuAj+acb6i5pB2aWVwotrbdYU9cHLONZb5xzDaO2cYx27I4412DuVb754FHjbDr\nscPHD036wFaSJGmtcR3vVTbXau8BfAa47xiHeQ+4Nmcks41jtrHMN47ZxjHbOGZbFme8V98TGAy6\nLwM+PcL+F424nyRJkmpkj/fq+53h46kziwunjXMg38HGMds4ZhvLfOOYbRyzjWO2ZQlfTlC3m2u1\n7ws8HrgJeGfN5UiSJGkVuY736no2g8w/NLO48KNxD+banHHMNo7ZxjLfOGYbx2zjmG1ZnPFeJcMP\nVT57uPm2OmuRJEnS6gsfeNvjfZvHAwcy+HDkOVUc0L6tOGYbx2xjmW8cs41jtnHMtizOeK+e3x0+\n/t3M4kKutRJJkiStOnu8V0HUhyrt24pjtnHMNpb5xjHbOGYbx2zL4oz36ngWt3+o8r/rLkaSJEmr\nzx7vYMMPVf72cPPvqjy2fVtxzDaO2cYy3zhmG8ds45htWbxz5RjmWu0Z4Em7eNn+3P6hyrOja5Ik\nSVKZwgfe8/PzTE9PR59m1c212r8JvGsFu7y16g9VppQ2+042htnGMdtY5hvHbOOYbRyzLYsz3iOY\na7UfBbx9uPkG4Pxd7HId8PHQoiRJklS08IF303q851rt+wMfBvYC3jSzuPDyumrxHWwcs41jtrHM\nN47ZxjHbOGZbFlc1WYG5Vns/4N+AnwU+Bbyo3ookSZI0KdZkj/dcq/0zrPy/PQH/CBwGXACcOLO4\ncEvVta2oIPu2wphtHLONZb5xzDaO2cYx27KsuR7vuVb7NcAfj3GIHwFPmllcuLqikiRJkrQGrKke\n7+FM98nDzStGOMR/A78zs7jwneqqGp3vYOOYbRyzjWW+ccw2jtnGMduyrLUZ798D9gE+NbO4cFzd\nxUiSJGntCP9w5fz8fPQpdstcq70P8ILh5p/VWUtVUkqb666hqcw2jtnGMt84ZhvHbOOYbVnW0qom\nzwLuCfQlUI6EAAAX80lEQVTxDpKSJElaZeED7xJ6vOda7T24vbf7z6q+g2Rd7NuKY7ZxzDaW+cYx\n2zhmG8dsy7JWZryfCmwELgI+Um8pkiRJWosa3+M912on4P8fbnZnFhdurbOeKtm3Fcds45htLPON\nY7ZxzDaO2ZZlLcx4HwscAVwGvKvmWiRJkrRGTcQ63nOt9jpGf5Pw0uHjG2cWF24cu5iC2LcVx2zj\nmG0s841jtnHMNo7ZlmVS1vE+CTh9jP2vAf62mlIkSZKklZuUHu8M3Dri1yLwmibe4t2+rThmG8ds\nY5lvHLONY7ZxzLYsEzHjPbO4cAZwRt11SJIkSaNaE+t4N5V9W3HMNo7ZxjLfOGYbx2zjmG1Z1sKq\nJpIkSVLtxh54p5SOSyldmFK6KKX00jv+vO51vJvMvq04ZhvHbGOZbxyzjWO2ccy2LGMNvFNK64E3\nAccBhwEzKaV2FYVJkiRJTTLujPdRwNac88U555uB9wFPXv4Ce7zj2LcVx2zjmG0s841jtnHMNo7Z\nlmXcgfeBwPeXbV8yfE6SJEnSMuMOvPOuXmCPdxz7tuKYbRyzjWW+ccw2jtnGMduyjDvwvhQ4eNn2\nwQxmvW9zzjnnkFI6PaX06uHXC5dfBCmlzW677fba2QY2jbO/2+Zb1zawqaR63Hbb7VXffmG6fTx7\n+iiTyynnXU5a73jnlPYAvglMAz8AvgzM5JwXll6zZcuWPD09nUY+iSRJklSYUca4Y925Mud8S0rp\n94FPAeuBty8fdEuSJEkaGPuW8TnnTwCf2NHP5+fnmZ6eHuscnW7v54FHjLj7zcAn+7NT14xVRIFS\nSpv9tHIMs41jtrHMN47ZxjHbOGZblrEH3qvkMcDpY+x/YafbO7o/O3VFRfVIkiRJKxI+8K5oHe/v\nAu8dcd+jgAcDZ3W6vWP6s1NXVlFQCXwHG8ds45htLPONY7ZxzDaO2ZZlIma8+7NTnwU+O8q+nW7v\n3sA5wC8Bn+p0e4/tz05dXWV9kiRJ0q6Mu5zgLtW9jnd/dupyBquufAd4GPDxTrd311qLqsjyZW5U\nLbONY7axzDeO2cYx2zhmW5aJmPEeV3926tJOtzfNYNb8kcAnO93eWSMc6jrgjP7s1H9XWqAkSZIa\nb6x1vHdHSet4d7q9BzAYfB8wxmGuBl4DvKk/O3VTJYVJkiRpoqz6Ot6Tpj87dVGn23s48JvAXiMc\n4uHAscBpwHM73d7JwL/1Z6di371IkiRp4oUPvKtYx7tK/dmp7wGvH3X/Trf3BOAvgAcB/wpc1On2\nfrqL3a4EXtyfnTp31PNuj2tzxjHbOGYby3zjmG0cs41jtmUJ/3Bl0/Rnpz4OHA68ELgKeACDFVN2\n9rUZ+Fin2zuwhpIlSZJUgDXV41214eooh+ziZQl4I3A08BXgV/qzU7uaIZ9onW5vE/BEYP0uXnoL\n8O7hXyEkSZImhj3eq6w/O3UdcN6uXtfp9p4CfJnBcobv6HR7T2tiX3in27sHgzae32XwhmN3PILB\nIF2SJKnRGr+Odwn6s1M/Ap7EYDnC/wO8oorjlrI2Z6fbW9fp9n4b+BbwXGAb8A8MVn/Z0dfrhq87\nttPt3b2OunemlGybyGxjmW8cs41jtnHMtizOeK+S/uzUNzrd3gyDD2S+rtPtXcpgFnylLunPTl1T\nbXXQ6fbuC7wB2DDC7hsZ9LIDnA38fn926hu7cc5HAccAxwPvHOG8kiRJEyN84L1p06boU0yM/uzU\nv3W6vZcBfwb844iHubrT7R3Tn506t6pPKXe6vfXAe4BHj3GYHwAnA/+0gjaa9zMYeP86hQ28/QR4\nHLONZb5xzDaO2cYx27I44736usA+wFNH2PeuwMHAP3W6vYf2Z6euraimFzEYdF8GPJ9BC8hK3Ax8\nboR6PgS8mUG7yX792amrV7i/JEnSxFhz63jXbTgb/Orh14p0ur27AF9isJzhW9K69W/P2249e5x6\nOt3eL3D7uubPHi6XuCr6s1OXd7q9zzJY8eV44F2rde5dcd3TOGYby3zjmG0cs41jtmVxHe8JMlyG\n8ETgBuA3DzzuWceNc7xOt7cncAaDu3j+w2oOupd5//DxKTWcW5IkadWED7zt8a5Wf3ZqAfh9gPsc\nM/OCTrfXHuNwrwSOBL7HoD+7Dh8CMvC4Tre3b0013ImzA3HMNpb5xjHbOGYbx2zL4oz3ZDqdwYch\n92HQ732XlR6g0+09FPij4eZvRayUsjv6s1M/BD4PtBgsuShJktRI9nhPoP7sVO50e8+7+ptfPnq/\nBx11OLDQ6fZWOnA+kMGdJf+qPzt1duVFrsz7gccwaDd5T821APbERTLbWOYbx2zjmG0csy2Lq5pM\nqP7s1LU/89BT/u9+DzrqNOB+Ix7mAuDlFZY1qg8Cfw08vtPt3a3C1VrG1un2EoMlD39uN17+nf7s\n1BeDS5IkSRMq5Rx75/JR7mOv3dfp9vYH7jvi7t/uz05dX2U9o+p0e58HHgXM9Gen3ld3PQCdbu9+\nwNuAx+3mLhno9GenvhpXlSRJKsEoY1xnvCdcf3bqKuCquuuowPsZDLx/Hah14D28odDzgVMY3Mnz\nJ8AnGAysd+RQ4JcZfGD1f0fXKEmSJo893hOsYX1bHwTeCDyh0+2dyMpv4lOV9cAfXPPt+Ufse8gm\nGLwheEF/duryne3U6fYOAC4GTuh0e7/Yn536enilE6ph121xzDeO2cYx2zhmWxZnvFWE/uzUJZ1u\n79+BR1DzjDdAvvWWHwO/3Z+d+sjuvL4/O3VZp9v7ewZLPb4CeFpkfZIkafKED7xdxztOA9/Bvgh4\nIfW/Ifz2fg982BuGbTwr8efA7wIndrq9V/dnp74VUNvEa+B1WxTzjWO2ccw2jtmWpe4BjnSb/uzU\nl4CZuusYVX926vudbu904DkMVot5Zr0VSZKkkoTfQGd+fj76FGtWSmlz3TU01RjZvgG4FTip0+1t\nrKygBvG6jWW+ccw2jtnGMduyeOdKqUL92anvMLgJ0HrgZTWXI0mSChI+8LbHO459W3HGzPYUBksP\nPrPT7R1UTUXN4XUby3zjmG0cs41jtmWxx1uqWH926sJOt/d+4KnApzrd3iUjHOb9/dmpf6i4NEmS\nVCN7vCeYfVtxKsj2dcAtwGHAsSN8vanT7d1jzBqK5HUby3zjmG0cs41jtmVxxlsK0J+dOr/T7R0B\nHDzC7q8EHgOcBPxVpYVJkqTapJx3dhfs8Y1yH3tpLet0e78GfAD4BnB4f3Yq9n9SSZK0YqOMcV3V\nRCrPR4ErgF8AHl5zLZIkqSL2eE8w+7bi1Jltf3bqJuD04eZz6qojitdtLPONY7ZxzDaO2ZbFGW+p\nTEsrmpzY6fb2rbUSSZJUCdfxnmCuzRmn7mz7s1MXAWcD+wBPq7OWqtWdbdOZbxyzjWO2ccy2LK5q\nIpXr74HNwO8Af1vFATvd3gbgYQzurLlaru7PTn11Fc8nSVKRwgfe8/PzTE9PR59mTUopbfadbIxC\nsv0QcCXwkE6399BxB6+dbu+RwHuB+1VR3ArP/dv92am3QzHZNpb5xjHbOGYbx2zL4oy3VKj+7NSN\nnW7vDOCFDD5kOdLAu9PtrQdeDryawUz3VuD7FZW5K3sDjwD+vNPt/Ut/dupHq3ReSZKKEz7wtsc7\nju9g4xSU7d8zGHg/rdPtvaQ/O3XdSnbudHsHAe8Gjh4+1QX+aLhySrhOt5eATwNTwJ8Cv1NQto1k\nvnHMNo7ZxjHbsngDHalwnW7vC8AjgR8CN65w93sBG4DLgaf3Z6fOrLi8Xep0ew8G/pPBG/1H9Gen\nvlThsdcBe46yb392arGqOiRJa0+RN9BxHe84rs0Zp7Bs/2L4eB/g/iv82gB8AvilOgbdAP3ZqQuB\n04AEvGXdnq2pKo7b6faewuANxY2jfHW6vXdWUUdpCrt2G8Vs45htHLMtiz3eUuH6s1Mf7HR7BzAY\nRK/UTf3ZqUuqrmkEr2OwLOKRBx737OOB3qgH6nR7ewOnAr83fOpmYKV/utsLeHqn23tLlTPwkiTt\njK0mklZFp9s7gcFKLVcDD+rPTl0+wjEOBf4ZeAiDAffJwJv6s1Mr+oes0+29HngFcFZ/durYldYh\nSdIoY1xnvCWtlo8waHt5PPBPnW7viyvcfy8Ga5rfDfgu8NT+7NRXRqzlNOD3gV/tdHu/0p+d+uyI\nx5EkabeN3OOdUvr1lNI3Ukq3ppSO3NHr7PGOY99WHLOt3nBW+gXXbD33ZgarrLx8hV8nMxh0fxA4\ncoxBN/3ZqSsZDL4BXjdcfaURvHbjmG0cs41jtmUZZ8b7fOAE4G0V1SKp4fqzU9++16P++iX7HvqQ\nu454iAuBD6+0tWQH3gj8IfAY4LHAWRUcU5KkHRq7xzul9Bng5Jzz17b3c3u8JZWq0+29FHgD8GXg\n4RUN6CVJa0CRywlKUsHeBFwBHAX8fzXXIklquJ0OvFNKZ6WUzt/O15N29wT2eMexbyuO2cYpKdv+\n7NT1wCnDzdcOb8gz0UrKt2nMNo7ZxjHbsuy0xzvn/KvjnuCcc87hJS95yenAxcOnrgLml25hunRB\nuO12SdtLSqmnYdubgGLq2WPD/hdsetUHLwU2XXPR165+0O98ddu+D3joLQDXXPTVPQAmafugJzxn\nfafbW6zh/P3z//wZf7P4o0uuL+x6q2wb2JRSKqYet93ene0lpdQz4dubgP2HkW489dRTmZ6eZiWq\n6vF+Sc75q9v7uT3ekkrX6fZOAs6ou44G+Bjw5P7s1K11FyJJ0VZ1He+U0gnAXwP3BD6WUjo35/z4\nUY8nSXXpz069q9PtfYTBWuFaufsA5zDok38N8Mp6y5GkMoXfufK0007LJ598sjPeAVJKm5f9CVQV\nMts4Zhurrnw73d408ClgPXBif3bqn1e7hmheu3HMNo7ZxnFVE0lSLfqzU1uAlww3/7HT7R1RZz2S\nVKLwW8Zv2rQp+hRrlu9g45htHLONVXO+f8Xgw0fPAP6l0+0dD9ywwmNk4Hv92albqi5uXF67ccw2\njtmWxRlvSVIlhjcgei7QB+4HnAdctMKvrcBHO92eLYqSGid84O063nHuuFSQqmO2ccw2Vt359men\nbgROAM4Evj3C1yJwHHDsate+K3Vn22RmG8dsyxLeaiJJWlv6s1OXAo8bZd9Ot/cSoAu8odPtndWf\nndpWaXGSVKPwGW97vOPYtxXHbOOYbawG5Ptm4BIGveJPrbmW/6EB2RbLbOOYbVns8ZYkFaM/O/VT\n4NXDzdd1uj3XVpfUGPZ4TzD7tuKYbRyzjdWQfN8JXAgcAjy75lpu05Bsi2S2ccy2LM54S5KKMlxK\ncOnul3/S6fY21FmPJFXFHu8JZt9WHLONY7axGpTvh4EvAwcAf1hzLUCjsi2O2cYx27K4qokkqTj9\n2anc6fZeBvSAl3a6vUsAVziRVIw3HLnyfcIH3vPz80xPT0efZk1KKW32nWwMs41jtrGalG9/duoz\nnW7vTAZrer+z7nqu+fY8+x7iX3EjmG0csy2LM96SpJI9l0G/9951F3LLtT+5N3B53XU0kdnGMdtQ\nv7HSHVLOOaKQ22zZsiVPT097619JkiQ1xihjXFc1kSRJklaB63hPMNfmjGO2ccw2lvnGMds4ZhvH\nbMvijLckSZK0ClzHe4I1ZeWCEpltHLONZb5xzDaO2cYx27I44y1JkiStAnu8J5h9W3HMNo7ZxjLf\nOGYbx2zjmG1ZnPGWJEmSVoE93hPMvq04ZhvHbGOZbxyzjWO2ccy2LM54S5IkSavAHu8JZt9WHLON\nY7axzDeO2cYx2zhmWxZnvCVJkqRVYI/3BLNvK47ZxjHbWOYbx2zjmG0csy2LM96SJEnSKrDHe4LZ\ntxXHbOOYbSzzjWO2ccw2jtmWxRlvSZIkaRXY4z3B7NuKY7ZxzDaW+cYx2zhmG8dsy+KMtyRJkrQK\n7PGeYPZtxTHbOGYby3zjmG0cs41jtmVxxluSJElaBfZ4TzD7tuKYbRyzjWW+ccw2jtnGMduyOOMt\nSZIkrQJ7vCeYfVtxzDaO2cYy3zhmG8ds45htWZzxliRJklaBPd4TzL6tOGYbx2xjmW8cs41jtnHM\ntizOeEuSJEmrwB7vCWbfVhyzjWO2scw3jtnGMds4ZlsWZ7wlSZKkVWCP9wSzbyuO2cYx21jmG8ds\n45htHLMtizPekiRJ0iqwx3uC2bcVx2zjmG0s841jtnHMNo7ZlsUZb0mSJGkV2OM9wezbimO2ccw2\nlvnGMds4ZhvHbMvijLckSZK0CkYeeKeUuimlhZTSeSmlD6WU9tve6+zxjmPfVhyzjWO2scw3jtnG\nMds4ZluWcWa8zwR+Ied8BPAt4OXbe9HWrVvHOIV2wT6eOGYbx2xjmW8cs41jtnHMNsgok8sjD7xz\nzmflnLcNN78EHLS9111//fWjnkK7tn/dBTSY2cYx21jmG8ds45htHLMNct555614n6p6vJ8FfLyi\nY0mSJEmNs8fOfphSOgs4YDs/ekXO+aPD17wSuCnn/N7tHeOyyy4bu0jt0Ma6C2iwjXUX0GAb6y6g\n4TbWXUCDbay7gAbbWHcBDbax7gJ0u5RzHn3nlH4LeA4wnXO+cXuved7znpeXt5scccQRLjFYkfn5\nebMMYrZxzDaW+cYx2zhmG8dsqzM/P/8/2ks2bNjAW9/61rSSY4w88E4pHQecBhydc/7RSAeRJEmS\n1ohxBt4XAXsBVw6f+vec8/OrKkySJElqkrFaTSRJkiTtnrA7V6aUjkspXZhSuiil9NKo86xVKaWL\nU0r/mVI6N6X05brrmWQppXeklC5PKZ2/7Ll7pJTOSil9K6V0ZkrJ5ZhGsINsX51SumR47Z47bFvT\nCqWUDk4pfSal9I2U0tdTSn8wfN5rd0w7ydZrd0wppb1TSl9KKc2nlC5IKZ0yfN7rtgI7yddrtyIp\npfXDDJcWGVnRtRsy451SWg98E3gscCnQB2ZyzguVn2yNSil9F3hozvnKXb5YO5VSegxwHXBGzvnw\n4XN/Dvwo5/znwzeOd885v6zOOifRDrJ9FXBtzvkvai1uwqWUDgAOyDnPp5TuCnwV+F/AM/HaHctO\nsn0qXrtjSyntk3O+IaW0B/B54CXA8XjdVmIH+U7jtVuJlNKLgYcCd8s5H7/S8ULUjPdRwNac88U5\n55uB9wFPDjrXWraiT9Jq+3LOnwN+coenjwfeOfz+nQx+6WqFdpAteO2OLed8Wc55fvj9dcACcCBe\nu2PbSbbgtTu2nPMNw2/3AtYz+DfC67YiO8gXvHbHllI6CHgC8A/cnueKrt2ogfeBwPeXbV/C7f9o\nqRoZ+HRK6SsppefUXUwD3TvnfPnw+8uBe9dZTAO9IKV0Xkrp7f5JeXwppY3AQxjcRdhrt0LLsv2P\n4VNeu2NKKa1LKc0zuD4/k3P+Bl63ldlBvuC1W4W/BGaBbcueW9G1GzXw9hOb8R6Vc34I8Hjg94Z/\n0leAPOjH8pquzluB+wObgB8yWJZUIxq2QnwQ+MOc87XLf+a1O55hth9gkO11eO1WIue8Lee8CTgI\n+JWU0jF3+LnX7Ri2k+9mvHbHllJ6InBFzvlcdvDXg925dqMG3pcCBy/bPpjBrLcqknP+4fDxv4EP\nM2jvUXUuH/Z5klK6D3BFzfU0Rs75ijzE4M91XrsjSintyWDQ/a6c80eGT3vtVmBZtu9eytZrt1o5\n56uBjzHol/W6rdiyfB/mtVuJRwLHDz9jNwdMpZTexQqv3aiB91eAB6SUNqaU9gJOBP416FxrTkpp\nn5TS3YbfbwCOBc7f+V5aoX8FnjH8/hnAR3byWq3A8B+mJSfgtTuSlFIC3g5ckHN+47Ifee2OaUfZ\neu2OL6V0z6U2h5TSXYBfBc7F67YSO8p3aWA45LU7gpzzK3LOB+ec7w/8H6CXcz6JFV67Yet4p5Qe\nD7yRQWP/23POp4ScaA1KKd2fwSw3wB7Ae8x3dCmlOeBo4J4M+rP+BPgX4J+B+wIXA0/NOV9VV42T\najvZvgrYzODPnRn4LvC7y/rjtJtSSo8GPgv8J7f/afPlwJfx2h3LDrJ9BTCD1+5YUkqHM/gA2rrh\n17tyzt2U0j3wuh3bTvI9A6/dyqSUjgZOHq5qsqJr1xvoSJIkSasg7AY6kiRJkm7nwFuSJElaBQ68\nJUmSpFXgwFuSJElaBQ68JUmSpFXgwFuSJElaBQ68JUmSpFXgwFuSJElaBf8PDk7BLbQy3rYAAAAA\nSUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "aaf.cumulative_hazards_.plot()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.9" } }, "nbformat": 4, "nbformat_minor": 0 }