{ "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" ] }, { "cell_type": "code", "execution_count": 1, "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": [ "/Users/camerondavidson-pilon/.virtualenvs/data/lib/python2.7/site-packages/IPython/core/magics/pylab.py:161: UserWarning: pylab import has clobbered these variables: ['datetime']\n", "`%matplotlib` prevents importing * from pylab and numpy\n", " \"\\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": 34, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/camerondavidson-pilon/.virtualenvs/data/lib/python2.7/site-packages/ipykernel/__main__.py:2: ParserWarning: Falling back to the 'python' engine because the 'c' engine does not support regex separators (separators > 1 char and different from '\\s+' are interpreted as regex); you can avoid this warning by specifying engine='python'.\n", " from ipykernel import kernelapp as app\n" ] } ], "source": [ "data = pd.read_csv('../lifelines/datasets/gehan.dat', sep=\"\\s{1,3}\", \n", " header=None)\n", "data.columns = ['treatment', 'time', 'failure']" ] }, { "cell_type": "code", "execution_count": 35, "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", " \n", "
treatmenttimefailure
371121
381151
391171
401221
411231
\n", "
" ], "text/plain": [ " treatment time failure\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": 35, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data.tail()" ] }, { "cell_type": "code", "execution_count": 36, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " treatment time failure\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": 37, "metadata": { "collapsed": false }, "outputs": [], "source": [ "treatment = data['treatment'] == 2" ] }, { "cell_type": "code", "execution_count": 38, "metadata": { "collapsed": false }, "outputs": [], "source": [ "T = data['time']\n", "E = data['failure']\n", "t = np.linspace(0,40,150)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 39, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(0, 1.05)" ] }, "execution_count": 39, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuUAAAF3CAYAAAACfDUrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3XuUV/V97//nmwEZAiOiIDlykYtg\n0MGMiMolOcESdFIjnno0xpi22Bijuaz+1LTBc34LKU26YoImzTKtscaQk/YYI7+kIR5iaBNtraOe\njJZKANOSMEZGg1zEMIHh+vn9MV/ocJ3bHj4zX56PtVjf7/5+P/u937PdC17u+ey9I6WEJEmSpHz6\n5G5AkiRJOtkZyiVJkqTMDOWSJElSZoZySZIkKTNDuSRJkpSZoVySJEnKzFAuSZIkZWYolyRJkjIz\nlEuSJEmZGcolSZKkzPrm2vDQoUPTmDFjcm1ekiRJ6hYvvPDC5pTSsI6sky2Ujxkzhvr6+lyblyRJ\nkrpFRLzS0XWcviJJkiRlZiiXJEmSMjOUS5IkSZllm1MuSZLU2+3Zs4cNGzbQ3NycuxVlUFlZyciR\nI+nXr1+XaxnKJUmSOmnDhg1UVVUxZswYIiJ3OzqBUkps2bKFDRs2MHbs2C7Xc/qKJElSJzU3N3PG\nGWcYyE9CEcEZZ5xR2G9JDOWSJEldYCA/eRX5395QLkmS1EvdfvvtfPnLXz64fMUVV3DzzTcfXL7z\nzju57777eO2117j22msBWLlyJcuXLz84ZuHChSxevLjNbf3FX/xFgZ3DkiVLeO211wqtmXM7XWUo\nlyRJ6qVmzpxJXV0dAPv372fz5s2sXr364Pd1dXXMmDGDs846i6VLlwJHhvL2OlYoTymxf//+Dtcz\nlB+qzVAeEQ9HxBsR8bNjfB8R8ZWIWBcRL0XElOLblCRJ0uFmzJjBs88+C8Dq1auprq6mqqqKN998\nk127drF27VqmTJlCQ0MD1dXV7N69mwULFvDoo49SU1PDo48+CsCaNWuYNWsW48aN4ytf+coR25k/\nfz47d+6kpqaGG2+8kYaGBs4991z+4A/+gOrqal599VVWrFjB9OnTmTJlCtdddx1NTU0ALFq0iIsv\nvpjq6mpuueUWUkosXbqU+vp6brzxRmpqati5cydjxozhrrvuoqamhqlTp/Liiy9yxRVXMH78eB54\n4IGDvXzxi1/k4osv5oILLuDuu+8GoKGhgUmTJvHRj36U888/n8svv5ydO3cedTs9VaSUjj8g4r8C\nTcD/SilVH+X73wU+BfwucCnwlymlS9va8NSpU1N9fX2nmpYkSeoJ1q5dy6RJkwD4sx+sZs1rvym0\n/nlnncrdV51/3DFjx47ln/7pn/jhD39ISonGxkamT5/O4MGDmT9/Pk8//TQNDQ28//3v52c/+xlL\nliyhvr6e+++/H2iZvrJixQqefPJJtm/fzrnnnsuvf/3rI27zN2jQoINBu6GhgXHjxlFXV8e0adPY\nvHkz11xzDT/84Q8ZOHAg99xzD7t27WLBggVs3bqV008/HYDf//3f5wMf+ABXXXUVs2bNYvHixUyd\nOhWAMWPG8JnPfIbbbruN22+/nR//+Mc888wzNDc3U11dzcaNG1mxYgVLly7la1/7Gikl5s6dy5/+\n6Z8yevRozjnnHOrr66mpqeEDH/gAc+fO5cMf/vAR2yla62PggIh4IaXUoQ22eUvElNI/R8SY4wy5\nmpbAnoDnIuK0iPgvKaXXO9KIJEmSOm7GjBnU1dVRV1fHHXfcQWNjI3V1dQwePJiZM2e2q8aVV15J\n//796d+/P2eeeSYbN25k5MiRx13n7LPPZtq0aQA899xzrFmz5uD2du/ezfTp0wF48skn+cIXvsCO\nHTvYunUr559/PlddddVRa86dOxeAyZMn09TURFVVFVVVVfTv359t27axYsUKVqxYwYUXXghAU1MT\n//Ef/8Ho0aMZO3YsNTU1AFx00UU0NDS062fvKYq4T/kI4NVWyxtKnx03lO98/WVW/8W7Ctg8PDPg\nMn78tt8tpJaO7uqaEXzo0tG525Akqcdq64x2dzkwr3zVqlVUV1czatQo7r33Xk499VRuuummdtXo\n37//wfcVFRXs3bu3zXUGDhx48H1KiTlz5vDII48cMqa5uZmPf/zj1NfXM2rUKBYuXHjcWwge6KNP\nnz6H9NSnTx/27t1LSom77rqLj33sY4es19DQcMTP0JOnqhzNCb3QMyJuiYj6iKhva9pMe43Z80tm\n7nyykFo6ujWv/4bvr2zM3YYkSTqKGTNm8Pjjj3P66adTUVHB6aefzrZt23j22WeZMWPGEeOrqqrY\nvn17h7fTr18/9uzZc9Tvpk2bxjPPPMO6desA+O1vf8u///u/HwzgQ4cOpamp6eDFpp3t44orruDh\nhx8+OI2msbGRN95447jrdPbnPdGKOFPeCIxqtTyy9NkRUkoPAg9Cy5zy8//Hv3R969+4kvOBR2+a\n3vVaOqrrv/Zs7hYkSdIxTJ48mc2bN/OhD33okM+ampoYOnToEeMvu+wyPv/5z1NTU8Ndd93V7u3c\ncsstXHDBBUyZMoXPfe5zh3w3bNgwlixZwg033MCuXbsA+OxnP8vEiRP56Ec/SnV1NW9/+9u5+OKL\nD64zb948br31VgYMGHDwYtW2XH755axdu/bg1JhBgwbxt3/7t1RUVBxzncO3M2DAgHb/zCdSmxd6\nApTmlD9+jAs9rwQ+yX9e6PmVlNIlbdUs7ELPb1zZ8nrT/+l6LR3VgVD+6Mf8Hx9Jklo72kV+Ormc\nsAs9I+IRYBYwNCI2AHcD/QBSSg8Ay2kJ5OuAHUD7Ji9JkiRJAtp395Ub2vg+AZ8orCNJkiTpJOMT\nPSVJkqTMDOWSJElSZoZySZIkKbMibomY386tsGpp2+Pao3IwTJhTTC1JkiSpHcrjTPn+fTBoeDF/\nmt/K/dNIkiS1y+23386Xv/zlg8tXXHEFN99888HlO++8k/vuu4/XXnuNa6+9FoCVK1eyfPnyg2MW\nLlzI4sWLC+lnyZIlvPbaax3+rjOeeuop6urqCquXezvlEcolSZJOQjNnzjwYGPfv38/mzZtZvXr1\nwe/r6uqYMWMGZ5111sGnaR4eyovU2VC+b9++Dm/LUC5JkqQeYcaMGQefhrl69Wqqq6upqqrizTff\nZNeuXaxdu5YpU6bQ0NBAdXU1u3fvZsGCBTz66KPU1NTw6KOPArBmzRpmzZrFuHHj+MpXvnKw/n33\n3Ud1dTXV1dUHz8gfqHXA4sWLWbhwIUuXLqW+vp4bb7yRmpoadu7ceXDM0b4bM2YMn/nMZ5gyZQqP\nPfYYv/jFL6itreWiiy7i3e9+Ny+//DIAP/jBD7j00ku58MILee9738vGjRtpaGjggQce4Etf+hI1\nNTU8/fTTzJs3j9tuu41p06Yxbtw4nnrqKf7oj/6ISZMmMW/evIO9rFixgunTpzNlyhSuu+46mpqa\nABgzZgx33303U6ZMYfLkybz88stH3U53KY855ZIkSbn9cD78elWxNd8+Gd73+WN+fdZZZ9G3b19+\n9atfUVdXx/Tp02lsbOTZZ59l8ODBTJ48mVNOOeXg+FNOOYVFixZRX1/P/fffD7RMX3n55Zd58skn\n2b59O+eeey633XYbL730Et/4xjd4/vnnSSlx6aWX8p73vIchQ4YctZdrr72W+++/n8WLFzN16tR2\nfXfGGWfw4osvAjB79mweeOABJkyYwPPPP8/HP/5xfvKTn/Cud72L5557jojgoYce4gtf+AL33nsv\nt956K4MGDeLTn/40AF//+td58803efbZZ1m2bBlz587lmWee4aGHHuLiiy9m5cqVjBw5ks9+9rP8\n4z/+IwMHDuSee+7hvvvuY8GCBQAMHTqUF198kb/6q79i8eLFPPTQQ0dsp7sYyiVJknqxGTNmUFdX\nR11dHXfccQeNjY3U1dUxePBgZs6c2a4aV155Jf3796d///6ceeaZbNy4kX/5l3/h937v9xg4cCAA\n11xzDU8//TRz584trPfrr78egKamJurq6rjuuusOfrdr1y4ANmzYwPXXX8/rr7/O7t27GTt27DHr\nXXXVVUQEkydPZvjw4UyePBmA888/n4aGBjZs2MCaNWsO7pfdu3czffr0g+tfc801AFx00UV897vf\nLeznbA9DuSRJUhGOc0a7Ox2YV75q1Sqqq6sZNWoU9957L6eeeio33XRTu2r079//4PuKigr27t17\nzLF9+/Zl//79B5ebm5s73fuBwL9//35OO+00Vq5cecSYT33qU9xxxx3MnTuXp556ioULFx6z3oGf\no0+fPof8TH369GHv3r1UVFQwZ84cHnnkkeOu39Y+6A6GcrXLtp17WLaysct1qgb047JzzyygI0mS\nBC1nyhcvXsy4ceOoqKjg9NNPZ9u2baxevZq/+Zu/OWJ8VVUV27dvb7Puu9/9bubNm8f8+fNJKfG9\n732Pb33rWwwfPpw33niDLVu2MGjQIB5//HFqa2vbrH2870499VTGjh3LY489xnXXXUdKiZdeeol3\nvvOdvPXWW4wYMQKAb37zm4fU+81vftPmz9HatGnT+MQnPsG6des455xz+O1vf0tjYyMTJ0485jqd\n2U5neKGn2mXfvv0Mq6rs8p/tO/fk/lEkSSorkydPZvPmzUybNu2QzwYPHszQoUOPGH/ZZZexZs2a\nQy70PJopU6Ywb948LrnkEi699FJuvvlmLrzwQvr168eCBQu45JJLmDNnDu94xzsOrjNv3jxuvfXW\nIy70bOs7gL/7u7/j61//Ou985zs5//zz+f73vw+0zHm/7rrruOiiiw75ea666iq+973vdegCzGHD\nhrFkyRJuuOEGLrjgAqZPn37wgtJj6cx2OiNSSt1W/HimTp2a6uvru17oG1fCbzfBlfd2vRZA00aY\nfG0xtcrE9V97li1Nu/jz/za5y7U2bW9mbs2IArqSJCm/tWvXMmnSpNxtKKOjHQMR8UJKaeoxVjkq\nz5RLkiRJmRnKJUmSpMwM5ZIkSVJmhnJJkqQuyHV9nvIr8r+9oVySJKmTKisr2bJli8H8JJRSYsuW\nLVRWVhZSz/uUS5IkddLIkSPZsGEDmzZtyt2KMqisrGTkyJGF1DKUH25PM6xa2vU6lYNhwpyu1ykz\nzXv3FfIQop7ChyFJ0smtX79+x33su9RehvLDDTm7mDpNG4upU2ZGDRmYu4VCbdre+UcLS5IkHeCc\nckmSJCkzQ7kkSZKUmaFckiRJyqw85pS/9St4Yn4xtcbNgom1xdSSJEmS2qH3h/LJ18JvC7oN0db1\nLa+GckmSJJ1AvT+UT70J+lfBoOFdr1XU2XZJkiSpA5xTLkmSJGXW+8+U91SdfQiRDx2SJEk66RjK\nu0tnH0LkQ4ckSZJOOoZyqQua9+5j2crG3G0UompAPy4798zcbUiSdFIylEtdMGrIwNwtFGbT9ubc\nLUiSdNLyQk9JkiQpM0O5JEmSlJmhXJIkScrMUC5JkiRlZiiXJEmSMjOUS5IkSZkZyiVJkqTMvE+5\n2qVx204WPb66y3Vmjh/K7EnDC+hIkiSpfHimXG26umYEI04b0OU6r2zZwTO/2FxAR5IkSeXFM+Vq\n04cuHc2g/hUMq6rsUp0izrRLkiSVI0O5JACa9+5j2crG3G0UpmpAPy4798zcbUiS1C6GckkAjBoy\nMHcLhdq0vTl3C5IktZtzyiVJkqTMDOWSJElSZoZySZIkKTNDuSRJkpSZoVySJEnKzFAuSZIkZeYt\nEXuaPc2wamkxtSoHw4Q5xdSSJElStzGU9zRDzi6uVtPG4mpJkiSp2zh9RZIkScqsXWfKI6IW+Eug\nAngopfT5w74fDXwTOK00Zn5KaXnBvZ4YW9fDE/O7XmfcLJhY2/U6kiRJKntthvKIqAC+CswBNgA/\njYhlKaU1rYb9v8B3Ukp/HRHnAcuBMd3Qb/caN6uYOlvXt7wayo/wypYdLHp8de42CjNz/FBmTxqe\nuw1JktTLtedM+SXAupTSLwEi4tvA1UDrUJ6AU0vvBwOvFdnkCTOxtpggXcSZ9jI0c/xQYHPuNgrz\nypYdwGZDuSRJ6rL2hPIRwKutljcAlx42ZiGwIiI+BQwE3ltIdyorsycNL6sAW05n/MtR8959LFvZ\nmLsNHUXVgH5cdu6ZuduQpB6lqLuv3AAsSSndGxHTgW9FRHVKaX/rQRFxC3ALwOjRowvatCQdadSQ\ngblb0DFs2t6cuwVJ6nHac/eVRmBUq+WRpc9a+wjwHYCU0rNAJTD08EIppQdTSlNTSlOHDRvWuY4l\nSZKkMtOeUP5TYEJEjI2IU4APAssOG/MrYDZAREyiJZRvKrJRSZIkqVy1GcpTSnuBTwI/AtbScpeV\n1RGxKCLmlobdCXw0Iv4NeASYl1JK3dW0JEmSVE7aNae8dM/x5Yd9tqDV+zXAzGJbkyRJkk4OPtFT\nkiRJysxQLkmSJGVmKJckSZIyM5RLkiRJmRnKJUmSpMwM5ZIkSVJm7bolojph63p4Yn7H1xs3CybW\nFt2NJEmSejBDeXcYN6tz621d3/JqKJckSTqpGMq7w8TazgXrzpxZlyRJUq9nKJe64JUtO1j0+Orc\nbRRi5vihzJ40PHcbkiSdlAzlUifNHD8U2Jy7jUK8smUHsNlQLklSJobycranGVYtLaTUWa++yakD\n+hVSq7P2nnIqb414T9YeWps9aXjZhNhyOdsvSVJvZSgvZ0POLqzUrsoK9gw4pbB6ndFv56as25dU\njOa9+1i2srGQWlUD+nHZuWcWUkuScjKUS5JOqFFDBhZWa9P25sJqSVJOPjxIkiRJysxQLkmSJGVm\nKJckSZIyM5RLkiRJmRnKJUmSpMwM5ZIkSVJmhnJJkiQpM0O5JEmSlJmhXJIkScrMUC5JkiRlZiiX\nJEmSMjOUS5IkSZkZyiVJkqTMDOWSJElSZn1zNyCpZ3hlyw4WPb46dxuFmTl+KLMnDc/dhiRJ7WIo\nl8TM8UOBzbnbKMwrW3YAmw3lkqRew1AuidmThpdVgC2nM/6SpJODc8olSZKkzAzlkiRJUmZOX5Ek\n9VrNe/exbGVj7jYk6RB9Bpw6pKPrGMolSb3WqCEDc7cgSUdI+/fu7eg6Tl+RJEmSMiuPM+WVg6Fp\nYzG19jTDkLOLqdUZW9fDE/OLqTVuFkysLaaWJEmSuk15hPIJc4qrtWppcbU6atys4mptXd/yaiiX\nJEnq8cojlJeLibXFheiizrZLkiSp2zmnXJIkScrMUC5JkiRlZiiXJEmSMjOUS5IkSZkZyiVJkqTM\nDOWSJElSZoZySZIkKTNDuSRJkpSZoVySJEnKzCd6SipLr2zZwaLHV+duQ0cxc/xQZk8anrsNSepR\nDOWSys7M8UOBzbnb0FG8smUHsNlQLkmHMZRLKjuzJw039PVQ/vZCko7OUK52eVv/vmzbsbvLdXbv\n28+ZVZUFdCRJklQ+DOVqlwtHnVZInWfWOaVAkiTpcO26+0pE1EbEzyNiXUTMP8aYD0TEmohYHRH/\nu9g2JUmSpPLV5pnyiKgAvgrMATYAP42IZSmlNa3GTADuAmamlN6MiDO7q2FJkiSp3LTnTPklwLqU\n0i9TSruBbwNXHzbmo8BXU0pvAqSU3ii2TUmSJKl8tSeUjwBebbW8ofRZaxOBiRHxTEQ8FxG1RTUo\nSZIklbuiLvTsC0wAZgEjgX+OiMkppW2tB0XELcAtAKNHjy5o05IkSVLv1p4z5Y3AqFbLI0uftbYB\nWJZS2pNSWg/8Oy0h/RAppQdTSlNTSlOHDRvW2Z4lSZKkstKeUP5TYEJEjI2IU4APAssOG/P3tJwl\nJyKG0jKd5ZcF9ilJkiSVrTZDeUppL/BJ4EfAWuA7KaXVEbEoIuaWhv0I2BIRa4AngT9JKW3prqYl\nSZKkctKuOeUppeXA8sM+W9DqfQLuKP2RJEmS1AHteniQJEmSpO5jKJckSZIyK+qWiFK367OvmTPW\n/yB3G4XYe8qpvDXiPbnbkCRJPYShvJxtXQ9PzM/dxSGqd+6hb0XnfkHz1ttn8ObI2QV3lEe/nZty\ntyBJknoQQ3m5GjcrdweFqtz+CkDZhHJJkqTWDOWHqxwMTRu7XmdPMww5u+t1OmtibcufHuZn6zZz\n2ttO6fB6Y+r/vBu6kSRJ6hkM5YebMKeYOquWFlNHkiRJZc9QrhPqbf37sm3H7g6vt3fffoBOrdud\ndu/bz5lVlbnbkCRJvZyhXCfUhaNO69yK6/oBMPOcoQV203XPrNucuwVJklQGvE+5JEmSlJmhXJIk\nScrMUC5JkiRlZiiXJEmSMjOUS5IkSZl59xVJ0gn1ypYdLHp8dSG1Zo4fyuxJwwupJUk5GcolSSfM\nzPFDgWJuJfrKlh3AZkO5pLJgKJcknTCzJw0vLEQXdbZdknoC55RLkiRJmRnKJUmSpMwM5ZIkSVJm\nzinvLpWDoWljx9fb0wxDzi6+H0mSJPVYhvLuMmFO59ZbtbTYPiRJktTjOX1FkiRJysxQLkmSJGVm\nKJckSZIyM5RLkiRJmXmhp9QFb+vfl207dnd4vf7Ne9i0vbkbOuq85r37GDVkYO42JEk6KRnKpS64\ncNRpnVuxaR9MHlFsM120bGVj7hYkSTppOX1FkiRJysxQLkmSJGVmKJckSZIyM5RLkiRJmXmhp5TD\nnmZYtbSYWpWDYcKcYmpJkqQsDOVSDkPOLq5W08biakmSpCycviJJkiRlZiiXJEmSMjOUS5IkSZkZ\nyiVJkqTMDOWSJElSZoZySZIkKTNDuSRJkpSZoVySJEnKzIcHSZJ6rVe27GDR46tztyFJXWYolyT1\nSjPHDwU2525DkgphKJck9UqzJw1n9qThuduQpCP86DMdX8c55ZIkSVJmhnJJkiQpM0O5JEmSlJlz\nyiUBUDWgH5u2N+duozDNe/cxasjA3G1IktQuhnJJAFx27pm5WyjUspWNuVuQJKndDOU9TeVgaNpY\nTK09zTDk7GJqSZIkqdsYynuaCXOKq7VqaXG1JEmS1G280FOSJEnKzFAuSZIkZdauUB4RtRHx84hY\nFxHzjzPuv0dEioipxbUoSZIklbc2Q3lEVABfBd4HnAfcEBHnHWVcFfDHwPNFNylJkiSVs/acKb8E\nWJdS+mVKaTfwbeDqo4z7c+AeoHxudCxJkiSdAO0J5SOAV1stbyh9dlBETAFGpZT+T4G9SZIkSSeF\nLl/oGRF9gPuAO9sx9paIqI+I+k2bNnV105IkSVJZaE8obwRGtVoeWfrsgCqgGngqIhqAacCyo13s\nmVJ6MKU0NaU0ddiwYZ3vWpIkSSoj7QnlPwUmRMTYiDgF+CCw7MCXKaW3UkpDU0pjUkpjgOeAuSml\n+m7pWJIkSSozbYbylNJe4JPAj4C1wHdSSqsjYlFEzO3uBiVJkqRy17c9g1JKy4Hlh3224BhjZ3W9\nLekotq6HJ455m/zeZdwsmFibuwtJktRDtCuUS9mNm5W7g+JsXd/yaiiXJEklhnL1DhNryyfElsvZ\nfkmSVBhDudTb7WmGVUtzd1GMysEwYU7uLiRJOuEM5VJvN+Ts3B0Up2lj7g4kScqiyw8PkiRJktQ1\nnimXVJaqBvRj0/bm3G3oKJr37mPUkIG525CkHsVQLqksXXbumblb0DEsW9nY9iBJOsk4fUWSJEnK\nzFAuSZIkZWYolyRJkjIzlEuSJEmZGcolSZKkzAzlkiRJUmaGckmSJCkz71NezioHl9djy/c0l9cj\n5SVJkkoM5eVswpzcHRRr1dLcHUiSJHULp69IkiRJmRnKJUmSpMwM5ZIkSVJmzimX1HPsae6Z1w5U\nDi6/azQkST2KoVxSz9FT765TTncxkiT1SE5fkSRJkjIzlEuSJEmZGcolSZKkzAzlkiRJUmaGckmS\nJCkzQ7kkSZKUmaFckiRJysxQLkmSJGXmw4OkHLauhyfmF1Nr3CyYWFtMLUmSlIWhXDrRxs0qrtbW\n9S2vhnJJkno1Q7l0ok2sLS5EF3W2XZIkZeWcckmSJCkzQ7kkSZKUmaFckiRJysxQLkmSJGXmhZ6S\npBOqakA/Nm1vLqRW8959jBoysJBakpSToVySdEJddu6ZhdVatrKxsFqSlJOhXJLasqcZVi3N3UX5\nqBwME+bk7kKSehRDuSS1ZcjZuTsoL00bc3cgST2OF3pKkiRJmRnKJUmSpMwM5ZIkSVJmhnJJkiQp\nM0O5JEmSlJmhXJIkScrMUC5JkiRlZiiXJEmSMjOUS5IkSZkZyiVJkqTMDOWSJElSZoZySZIkKTND\nuSRJkpRZu0J5RNRGxM8jYl1EzD/K93dExJqIeCkifhwRZxffqiRJklSe2gzlEVEBfBV4H3AecENE\nnHfYsH8FpqaULgCWAl8oulFJkiSpXLXnTPklwLqU0i9TSruBbwNXtx6QUnoypbSjtPgcMLLYNiVJ\nkqTy1Z5QPgJ4tdXyhtJnx/IR4IddaUqSJEk6mfQtslhEfBiYCrznGN/fAtwCMHr06CI3LUmSJPVa\n7QnljcCoVssjS58dIiLeC/xP4D0ppV1HK5RSehB4EGDq1Kmpw91Kknq/Pc2wamkhpYb9eg9UXV5I\nLUnKqT2h/KfAhIgYS0sY/yDwodYDIuJC4GtAbUrpjcK7lHRsW9fDE0fcFKl3GjcLJtbm7kLdbUhx\nN+jqt+dl9hRWTZLyaTOUp5T2RsQngR8BFcDDKaXVEbEIqE8pLQO+CAwCHosIgF+llOZ2Y9+SoCXE\nlout61teDeXqgAGnVPDa9ubcbUjSIaJP3w5PEW/XCiml5cDywz5b0Or9ezu6YUkFmFhbPiG2XM72\n64Q6/6xTOX/y8e49IEkn3v6dv3mzo+sUeqGn1K0qB0PTxtxdFGNPc6G/wpckSb2boVy9x4Q5uTso\nTkEXuUmSpPLQnvuUS5IkSepGhnJJkiQpM0O5JEmSlJmhXJIkScrMUC5JkiRlZiiXJEmSMjOUS5Ik\nSZkZyiVJkqTMDOWSJElSZj7RU1LPsXU9PDE/dxdHGjcLJtbm7kKSVMYM5ZJ6hnGzcndwdFvXt7wa\nyiVJ3chQLqlnmFjbM4NvTzxzL0kqO84plyRJkjIzlEuSJEmZGcolSZKkzAzlkiRJUmaGckmSJCkz\nQ7kkSZKUmaFckiRJysz7lEvydcBlAAAJgUlEQVQ5VA6Gpo3F1NrTDEPOLqaWJEnKwlAu5TBhTnG1\nVi0trpYkScrC6SuSJElSZoZySZIkKTNDuSRJkpSZc8olqS1b18MT83N3UT7GzYKJtbm7kKQexVAu\nScczblbuDsrL1vUtr4ZySTqEoVySjmdirQGySP7GQZKOyjnlkiRJUmaGckmSJCkzQ7kkSZKUmaFc\nkiRJysxQLkmSJGVmKJckSZIyM5RLkiRJmRnKJUmSpMwM5ZIkSVJmPtFT6u0qB0PTxtxdFGNPMww5\nO3cXkiSdcIZyqbebMCd3B8VZtTR3B5IkZeH0FUmSJCkzz5RLkk6srevhifnF1DprCky+tphakpSR\noVySdOKMm1Vcra3rYd+e4upJUkaGcknSiTOxtuVPEZ6YbyiXVDacUy5JkiRlZiiXJEmSMnP6iqSe\no6fec937p0uSupmhXFLP0VPvue790yVJ3czpK5IkSVJmhnJJkiQpM0O5JEmSlJmhXJIkScrMCz0l\nqS099a4wvZV3s5GkI7QrlEdELfCXQAXwUErp84d93x/4X8BFwBbg+pRSQ7GtSlImPfWuML2Vd7OR\npCO0OX0lIiqArwLvA84DboiI8w4b9hHgzZTSOcCXgHuKblSSJEkqV+2ZU34JsC6l9MuU0m7g28DV\nh425Gvhm6f1SYHZERHFtSpIkSeWrPdNXRgCvtlreAFx6rDEppb0R8RZwBrC5iCYlSWWkqDn6+/bA\nW6/CN67sei1JyuyEXugZEbcAt5QWd0XEz07k9k8iQ/F/iLqL+7Z7uX+7T1nu28H9qYLlWXvYtY/B\n/St4K2sTZcp9233ct92rIjh8qneb2hPKG4FRrZZHlj472pgNEdEXGEzLBZ+HSCk9CDwIEBH1KaWp\nHW1YbXPfdh/3bfdy/3Yf9233iYj6nXvct93Bfdt93LfdKyLqO7pOe+aU/xSYEBFjI+IU4IPAssPG\nLAP+sPT+WuAnKaXU0WYkSZKkk1GbZ8pLc8Q/CfyIllsiPpxSWh0Ri4D6lNIy4OvAtyJiHbCVluAu\nSZIkqR3aNac8pbScwybtpZQWtHrfDFzXwW0/2MHxaj/3bfdx33Yv92/3cd92H/dt93Hfdh/3bffq\n8P4NZ5lIkiRJebVnTrkkSZKkbpQllEdEbUT8PCLWRcT8HD2Uq4hoiIhVEbGyM1f+6j9FxMMR8Ubr\nW3dGxOkR8Q8R8R+l1yE5e+ytjrFvF0ZEY+nYXRkRv5uzx94qIkZFxJMRsSYiVkfEH5c+99jtouPs\nW4/dAkREZUT834j4t9L+/bPS52Mj4vlSZni0dNMJdcBx9u2SiFjf6tityd1rbxURFRHxrxHxeGm5\nw8ftCQ/lEVEBfBV4H3AecENEdPhejjquy1JKNd7+rMuWALWHfTYf+HFKaQLw49KyOm4JR+5bgC+V\njt2a0rUs6ri9wJ0ppfOAacAnSn/Heux23bH2LXjsFmEX8DsppXcCNUBtREwD7qFl/54DvAl8JGOP\nvdWx9i3An7Q6dlfma7HX+2NgbavlDh+3Oc6UXwKsSyn9MqW0G/g2cHWGPqTjSin9My13E2rtauCb\npfffBP7bCW2qTBxj36oAKaXXU0ovlt5vp+UfiRF47HbZcfatCpBaNJUW+5X+JOB3gKWlzz12O+E4\n+1YFiIiRwJXAQ6XloBPHbY5QPgJ4tdXyBvxLrUgJWBERL5SeoKpiDU8pvV56/2tgeM5mytAnI+Kl\n0vQWp1d0UUSMAS4Ensdjt1CH7Vvw2C1EaQrASuAN4B+AXwDbUkp7S0PMDJ10+L5NKR04dj9XOna/\nFBH9M7bYm30Z+FNgf2n5DDpx3HqhZ/l5V0ppCi3Tgz4REf81d0PlqvSALM80FOevgfG0/Gr1deDe\nvO30bhExCPj/gP8npfSb1t957HbNUfatx25BUkr7Uko1tDw9/BLgHZlbKhuH79uIqAbuomUfXwyc\nDnwmY4u9UkS8H3gjpfRCV2vlCOWNwKhWyyNLn6kAKaXG0usbwPdo+UtNxdkYEf8FoPT6RuZ+ykZK\naWPpH439wN/gsdtpEdGPltD4dyml75Y+9tgtwNH2rcdu8VJK24AngenAaRFx4LkqZoYuarVva0tT\nslJKaRfwDTx2O2MmMDciGmiZkv07wF/SieM2Ryj/KTChdFXqKbQ8/XNZhj7KTkQMjIiqA++By4Gf\nHX8tddAy4A9L7/8Q+H7GXsrKgcBY8nt47HZKaS7j14G1KaX7Wn3lsdtFx9q3HrvFiIhhEXFa6f0A\nYA4t8/afBK4tDfPY7YRj7NuXW/2PetAy59ljt4NSSnellEamlMbQkml/klK6kU4ct1keHlS6XdSX\ngQrg4ZTS5054E2UoIsbRcnYcWp7W+r/dt50XEY8As4ChwEbgbuDvge8Ao4FXgA+klLxgsYOOsW9n\n0fLr/wQ0AB9rNQda7RQR7wKeBlbxn/Mb/wctc589drvgOPv2Bjx2uywiLqDlgrgKWk4afieltKj0\nb9u3aZle8a/Ah0tndtVOx9m3PwGGAQGsBG5tdUGoOigiZgGfTim9vzPHrU/0lCRJkjLzQk9JkiQp\nM0O5JEmSlJmhXJIkScrMUC5JkiRlZiiXJEmSMjOUS1IPFBGnRcTHS+/PioilBdVdGBGfLr1fFBHv\nLaKuJKlrvCWiJPVAETEGeDylVF1w3YVAU0ppcZF1JUld45lySeqZPg+Mj4iVEfFYRPwMICLmRcTf\nR8Q/RERDRHwyIu6IiH+NiOci4vTSuPER8UREvBART0fEOw7fQEQsiYhrS+8bIuLPIuLFiFh1YHzp\nScEPR8T/LW3j6hO4DyTppGEol6SeaT7wi5RSDfAnh31XDVwDXAx8DtiRUroQeBb4g9KYB4FPpZQu\nAj4N/FU7trk5pTQF+OvSOgD/k5bHRl8CXAZ8MSIGdv7HkiQdTd/cDUiSOuzJlNJ2YHtEvAX8oPT5\nKuCCiBgEzAAei4gD6/RvR93vll5foCX0A1wOzD0wDx2oBEYDa7v2I0iSWjOUS1Lvs6vV+/2tlvfT\n8vd6H2Bb6Sx7Z+ru4z//fQjgv6eUft7JXiVJ7eD0FUnqmbYDVZ1ZMaX0G2B9RFwHEC3e2ck+fgR8\nKkqn3CPiwk7WkSQdh6FcknqglNIW4JnSBZ5f7ESJG4GPRMS/AauBzl6g+edAP+CliFhdWpYkFcxb\nIkqSJEmZeaZckiRJysxQLkmSJGVmKJckSZIyM5RLkiRJmRnKJUmSpMwM5ZIkSVJmhnJJkiQpM0O5\nJEmSlNn/D8juRkxcLw/dAAAAAElFTkSuQmCC\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)\n", "ylim(0,1.05)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Hazards, too" ] }, { "cell_type": "code", "execution_count": 40, "metadata": { "collapsed": false }, "outputs": [], "source": [ "naf = NelsonAalenFitter()" ] }, { "cell_type": "code", "execution_count": 41, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuUAAAF3CAYAAAACfDUrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3XuclXW99//Xh9OgnJSDR0DwkAcG\nHREV0O4s85jSzq2WuSvbmWmHvdPapf3ujeau7iyydj93mWlh7TKLbaXeutNKdxZCjYQa4AGLYACB\n4TSDMAPDfO8/1mJCTnNgrbnWWvN6Ph48rnW4Dp+5vB7w9juf63tFSglJkiRJ2emVdQGSJElST2co\nlyRJkjJmKJckSZIyZiiXJEmSMmYolyRJkjJmKJckSZIyZiiXJEmSMmYolyRJkjJmKJckSZIyZiiX\nJEmSMtYnqwMPHz48jRkzJqvDS5IkSUXxzDPP1KeURnRmm8xC+ZgxY6itrc3q8JIkSVJRRMRfO7uN\n7SuSJElSxgzlkiRJUsYM5ZIkSVLGMusp352tW7dSV1dHU1NT1qUoA/3792fkyJH07ds361IkSZK6\nVUmF8rq6OgYNGsSYMWOIiKzLUTdKKbFmzRrq6uoYO3Zs1uVIkiR1q5JqX2lqamLYsGEG8h4oIhg2\nbJi/JZEkST1SSYVywEDeg/nfXpIk9VQlF8qzdP311/O1r32t7f15553H1Vdf3fb+E5/4BLfffjvL\nly/n0ksvBWDevHk88sgjbevccsstTJ8+vd1jfeELXyhg5TBjxgyWL19e0H1meRxJkqSexFC+gzPO\nOINZs2YB0NraSn19PfPnz2/7ftasWUyZMoXDDjuMmTNnAruG8o7aUyhPKdHa2trp/RnKJUmSypeh\nfAdTpkzh6aefBmD+/PlUV1czaNAg1q1bR3NzMwsXLmTChAksXryY6upqtmzZwrRp07j//vupqanh\n/vvvB2DBggWcddZZHHnkkXz961/f5Tg33ngjmzdvpqamhiuvvJLFixdz7LHH8t73vpfq6mqWLl3K\nY489xuTJk5kwYQKXXXYZGzduBODWW2/l1FNPpbq6mmuuuYaUEjNnzqS2tpYrr7ySmpoaNm/ezJgx\nY7jpppuoqalh4sSJzJ07l/POO4+jjjqKO++8s62WL3/5y5x66qmceOKJ3HzzzQAsXryY448/ng9+\n8IOMGzeOc889l82bN+/2OJIkSdp3kVLK5MATJ05MtbW1r/ts4cKFHH/88QB89qH5LFjeUNBjnnDY\nYG6+eNxe1xk7diz/8z//w6OPPkpKiWXLljF58mSGDBnCjTfeyFNPPcXixYu56KKL+NOf/sSMGTOo\nra3ljjvuAHLtK4899hhPPPEEjY2NHHvssbz66qu7TPM3cODAtqC9ePFijjzySGbNmsWkSZOor6/n\nkksu4dFHH2XAgAHcdtttNDc3M23aNNauXcvQoUMBeM973sPll1/OxRdfzFlnncX06dOZOHEiAGPG\njOHTn/401113Hddffz2/+tWv+N3vfkdTUxPV1dWsXLmSxx57jJkzZ/Ktb32LlBJTp07lU5/6FKNH\nj+boo4+mtraWmpoaLr/8cqZOnco//MM/7HKcQtvxGpAkSSpHEfFMSqlTYamkpkQsBVOmTGHWrFnM\nmjWLG264gWXLljFr1iyGDBnCGWec0aF9vO1tb6OqqoqqqioOOuggVq5cyciRI/e6zRFHHMGkSZMA\nmD17NgsWLGg73pYtW5g8eTIATzzxBF/60pfYtGkTa9euZdy4cVx88cW73efUqVMBGD9+PBs3bmTQ\noEEMGjSIqqoq1q9fz2OPPcZjjz3GySefDMDGjRt5+eWXGT16NGPHjqWmpgaAU045hcWLF3foZ5ek\nbrP6RVi/JOsqJKkgSjaUtzeiXSzb+8qff/55qqurGTVqFF/5ylcYPHgw73//+zu0j6qqqrbXvXv3\npqWlpd1tBgwY0PY6pcQ555zDfffd97p1mpqa+PCHP0xtbS2jRo3illtu2esUgtvr6NWr1+tq6tWr\nFy0tLaSUuOmmm/jQhz70uu0WL168y89gq4qkknP3OdC8IesqJKkgSjaUZ2XKlClMnz6dI488kt69\nezN06FDWr1/P/Pnz+fa3v73L+oMGDaKxsbHTx+nbty9bt27d7dMrJ02axEc+8hEWLVrE0UcfzWuv\nvcayZcs46KCDABg+fDgbN25k5syZbbPAdKWO8847j3/913/lyiuvZODAgSxbtqzdp2l29eeVpILa\n2pQL5BM/ADXvzroaSXq9z57a6U0M5TsZP3489fX1vPvd737dZxs3bmT48OG7rP/mN7+ZL37xi9TU\n1HDTTTd1+DjXXHMNJ554IhMmTODzn//8674bMWIEM2bM4IorrqC5uRmAz33uc7zhDW/ggx/8INXV\n1RxyyCGceurf/oNfddVVXHvttey3335tN6u259xzz2XhwoVtrTEDBw7kP//zP+ndu/cet9n5OPvt\nt1+Hf2ZJKpjm/ODAQcfDyOLc4yJJ3alkb/RUz+Q1IKlD1rwC//8EeMddcNI7s65Gkl6nKzd6OiWi\nJKn8NOV7yfsPzrYOSSoQQ7kkqfw056fMrRqUbR2SVCCGcklS+dneU17lSLmkymAolySVn6b8SLnt\nK5IqRLuhPCL6R8TvI+LZiJgfEZ/dzTpXRcTqiJiX/3N1ccqVJIkd2lcM5ZIqQ0emRGwG3pJS2hgR\nfYHfRsSjKaXZO613f0rpo4UvUZKknTTZUy6psrQ7Up5yNubf9s3/yWYexSK7/vrr+drXvtb2/rzz\nzuPqq/826P+JT3yC22+/neXLl7c9tGfevHk88sgjbevccsstTJ8+vSD1zJgxg+XLl3f6u6548skn\nmTVrVsH2l/VxJFW45gbouz/03vsDzySpXHSopzwiekfEPGAV8HhKac5uVvv7iHguImZGxKiCVtlN\nzjjjjLbA2NraSn19PfPnz2/7ftasWUyZMoXDDjuMmTNnAruG8kLqaijftm1bp49lKJdUVpobbF2R\nVFE6FMpTSttSSjXASOC0iKjeaZWHgDEppROBx4F7d7efiLgmImojonb16tX7UndRTJkype1pmPPn\nz6e6uppBgwaxbt06mpubWbhwIRMmTGDx4sVUV1ezZcsWpk2bxv33309NTQ33338/AAsWLOCss87i\nyCOP5Otf/3rb/m+//Xaqq6uprq5uG5Hfvq/tpk+fzi233MLMmTOpra3lyiuvpKamhs2bN7ets7vv\nxowZw6c//WkmTJjAT37yE1555RXOP/98TjnlFN74xjfywgsvAPDQQw9x+umnc/LJJ/PWt76VlStX\nsnjxYu68806++tWvUlNTw1NPPcVVV13Fddddx6RJkzjyyCN58skn+cd//EeOP/54rrrqqrZaHnvs\nMSZPnsyECRO47LLL2Lgx90uVMWPGcPPNNzNhwgTGjx/PCy+8sNvjSFKXNDXYuiKponSkp7xNSml9\nRDwBnA/8aYfP1+yw2t3Al/aw/V3AXZB7oudeD/bojfDq850pr32HjIcLvrjHrw877DD69OnDkiVL\nmDVrFpMnT2bZsmU8/fTTDBkyhPHjx9OvX7+29fv168ett95KbW0td9xxB5BrX3nhhRd44oknaGxs\n5Nhjj+W6667jueee47vf/S5z5swhpcTpp5/Om970Jg488MDd1nLppZdyxx13MH36dCZOnNih74YN\nG8bcuXMBOPvss7nzzjs55phjmDNnDh/+8If59a9/zZlnnsns2bOJCO6++26+9KUv8ZWvfIVrr72W\ngQMH8slPfhKAe+65h3Xr1vH000/z4IMPMnXqVH73u99x9913c+qppzJv3jxGjhzJ5z73OX75y18y\nYMAAbrvtNm6//XamTZsGwPDhw5k7dy7f+MY3mD59Onffffcux5GkLmlucOYVSRWl3VAeESOArflA\nvh9wDnDbTuscmlJakX87FVhY8Eq7yZQpU5g1axazZs3ihhtuYNmyZcyaNYshQ4ZwxhlndGgfb3vb\n26iqqqKqqoqDDjqIlStX8tvf/pZ3vOMdDBgwAIBLLrmEp556iqlTpxas9ne+M/eo6Y0bNzJr1iwu\nu+yytu+am5sBqKur453vfCcrVqxgy5YtjB07do/7u/jii4kIxo8fz8EHH8z48eMBGDduHIsXL6au\nro4FCxa0nZctW7YwefLktu0vueQSAE455RQeeOCBgv2ckkRzo+0rkipKR0bKDwXujYje5NpdfpxS\nejgibgVqU0oPAv8UEVOBFmAtcNU+V7aXEe1i2t5X/vzzz1NdXc2oUaP4yle+wuDBg3n/+9/foX1U\nVVW1ve7duzctLS17XLdPnz60tra2vW9qaupy7dsDf2trKwcccADz5s3bZZ2Pfexj3HDDDUydOpUn\nn3ySW265ZY/72/5z9OrV63U/U69evWhpaaF3796cc8453HfffXvdvr1zIEmd1tQAgw/LugpJKpiO\nzL7yXErp5JTSiSml6pTSrfnPp+UDOSmlm1JK41JKJ6WU3pxSeqHYhRfLlClTePjhhxk6dCi9e/dm\n6NChrF+/nqeffpopU6bssv6gQYNobGxsd79vfOMb+dnPfsamTZt47bXX+OlPf8ob3/hGDj74YFat\nWsWaNWtobm7m4Ycf7tC+9/bd4MGDGTt2LD/5yU8ASCnx7LPPArBhwwYOP/xwAO6992+t/x39OXY0\nadIkfve737Fo0SIAXnvtNV566aW9btOV40jSLprtKZdUWXyi507Gjx9PfX09kyZNet1nQ4YMYfjw\n4bus/+Y3v5kFCxa87kbP3ZkwYQJXXXUVp512GqeffjpXX301J598Mn379mXatGmcdtppnHPOORx3\n3HFt21x11VVce+21u9zo2d53AD/4wQ+45557OOmkkxg3bhw///nPgVzP+2WXXcYpp5zyup/n4osv\n5qc//WmnbsAcMWIEM2bM4IorruDEE09k8uTJbTeU7klXjiNJu2huhKohWVchSQUTKWUz5fjEiRNT\nbW3t6z5buHAhxx9/fCb1qDR4DUhqV+s2uHUonHUTnHVj1tVI0i4i4pmU0sT21/wbR8olSeWlefvT\nPL3RU1LlMJRLkspL0/ZQbk+5pMphKJcklZfm/M3izlMuqYKUXCjPqsdd2fO/vaQOsX1FUgUqqVDe\nv39/1qxZYzjrgVJKrFmzhv79+2ddiqRS12Qol1R5OvLwoG4zcuRI6urqWL16ddalKAP9+/dn5MiR\nWZchqdRtHym3fUVSBSmpUN63b9+9PvZdkiTbVyRVopJqX5EkqV1NjpRLqjyGcklSeWlugF59oI/3\noEiqHIZySVJ5aWrIta5EZF2JJBWMoVySVF6aG21dkVRxDOWSpPLS3OBNnpIqjqFcklRemgzlkiqP\noVySVF5sX5FUgQzlkqTy0rzBkXJJFcdQLkkqL00NjpRLqjiGcklS+Ugp175SNSjrSiSpoAzlkqTy\nsXUTpG22r0iqOIZySVL5aGrILW1fkVRhDOWSpPLRnA/ljpRLqjCGcklS+WgylEuqTIZySVL5aLZ9\nRVJlMpRLksqH7SuSKpShXJJUPtraV5wSUVJlMZRLksqH7SuSKpShXJJUPpobgYB+jpRLqiyGcklS\n+WhqyLWu9PKfL0mVxb/VJEnlo7nBfnJJFclQLkkqH80NzrwiqSIZyiVJ5aOpwZs8JVWkdkN5RPSP\niN9HxLMRMT8iPrubdaoi4v6IWBQRcyJiTDGKlST1cI6US6pQHRkpbwbeklI6CagBzo+ISTut8wFg\nXUrpaOCrwG2FLVOSJP52o6ckVZh2Q3nK2Zh/2zf/J+202tuBe/OvZwJnR0QUrEpJkiA3JaLtK5Iq\nUId6yiOid0TMA1YBj6eU5uy0yuHAUoCUUguwARhWyEIlSbJ9RVKl6lAoTyltSynVACOB0yKiuisH\ni4hrIqI2ImpXr17dlV1Iknqqli3Q0mQol1SROjX7SkppPfAEcP5OXy0DRgFERB9gCLBmN9vflVKa\nmFKaOGLEiK5VLEnqmZobckvbVyRVoI7MvjIiIg7Iv94POAd4YafVHgTel399KfDrlNLOfeeSJHXd\n9lDuSLmkCtSnA+scCtwbEb3Jhfgfp5QejohbgdqU0oPAPcD3I2IRsBZ4V9EqliT1TE2OlEuqXO2G\n8pTSc8DJu/l82g6vm4DLCluaJEk7aBspd0pESZXHJ3pKkspDk+0rkiqXoVySVB6aG3NL21ckVSBD\nuSSpPLS1rwzJtg5JKgJDuSSpPDTZUy6pchnKJUnlobkB+vSHPv2yrkSSCs5QLkkqD80N3uQpqWIZ\nyiVJ5aGpwZs8JVUsQ7kkqTw0N9hPLqliGcolSeWhudH2FUkVy1AuSSoPtq9IqmCGcklSefBGT0kV\nzFAuSSoPTYZySZXLUC5JKn2trbCl0fYVSRXLUC5JKn1bGnNLR8olVShDuSSp9DU15JZOiSipQhnK\nJUmlrzkfym1fkVShDOWSpNLXbPuKpMpmKJcklb7t7Sv9h2RbhyQViaFcklT6mu0pl1TZDOWSpNLX\nFsptX5FUmQzlkqTS1+SNnpIqm6FcklT6mhsgekPf/bOuRJKKwlAuSSp9TQ25fvKIrCuRpKLok3UB\nkqQe5IkvwNPf6Px2LZth8OGFr0eSSoShXJLUfRb9CvYfCsdd1PltR08qfD2SVCIM5ZKk7tOwDI46\nG87/QtaVSFJJsadcktQ9WrZA46swZGTWlUhSyTGUS5K6R+MKIMEQe8MlaWeGcklS99hQl1s6Ui5J\nuzCUS5K6R8Oy3HKwoVySdmYolyR1j7aRcttXJGlnhnJJUvfYUAf7HQj9BmRdiSSVnHZDeUSMiogn\nImJBRMyPiH/ezTpnRcSGiJiX/zOtOOVKkspWwzJbVyRpDzoyT3kL8ImU0tyIGAQ8ExGPp5QW7LTe\nUymlLjwNQpLUI2yogyGjsq5CkkpSuyPlKaUVKaW5+deNwELAhkBJUudsqLOfXJL2oFM95RExBjgZ\nmLObrydHxLMR8WhEjCtAbZKkStG8EZrWOx2iJO1BR9pXAIiIgcB/AR9PKTXs9PVc4IiU0saIuBD4\nGXDMbvZxDXANwOjRo7tctCSpzDgdoiTtVYdGyiOiL7lA/oOU0gM7f59Sakgpbcy/fgToGxHDd7Pe\nXSmliSmliSNGjNjH0iVJZWPD0tzSkXJJ2q2OzL4SwD3AwpTS7XtY55D8ekTEafn9rilkoZKkMrYh\nP1JuT7kk7VZH2lfOAN4DPB8R8/KffQYYDZBSuhO4FLguIlqAzcC7UkqpCPVKkspRwzKIXjDo0Kwr\nkaSS1G4oTyn9Foh21rkDuKNQRUmSKsyGOhh4CPTum3UlklSSfKKnJKn4NtTZTy5Je2EolyQVn3OU\nS9JeGcolScWVUq6n3JFySdojQ7kkqbg2rYGWJucol6S9MJRLkoprQ11u6Ui5JO2RoVySVFxtodye\ncknaE0O5JKm4GvIPDrJ9RZL2yFAuSSquDUuhdxUMGJ51JZJUsgzlkqTi2rAs17oSe30OnST1aIZy\nSVJxOR2iJLXLUC5JKq4NdfaTS1I7DOWSpOLZ1gKNKxwpl6R2GMolScXTuAJSq9MhSlI7DOWSpOLZ\nPh2iI+WStFeGcklS8Wx/cJA95ZK0V4ZySVLx+DRPSeoQQ7kkqXg21EH/IVA1KOtKJKmkGcolScXT\nsAyGjMq6CkkqeYZySVLxbKiDwbauSFJ7DOWSpOLZUOfMK5LUAYZySVJxbNkEm9d6k6ckdYChXJJU\nHG1zlNtTLknt6ZN1AZJU0upfhvV/zbqK8rRyfm5pT7kktctQLkl70roNvn02NG/IupIyFjD0yKyL\nkKSSZyiXpD1ZtzgXyN/0aTj6rVlXU572OxAGH5p1FZJU8gzlkrQnqxbklm84Dw4/JdtaJEkVzRs9\nJWlPVi3MLUccl20dkqSKZyiXpD1ZtQAOHAP9BmRdiSSpwhnKJWlPVi2Eg07IugpJUg9gKJek3Wlp\nzk2HaCiXJHUDQ7kk7U79y5C2wUHHZ12JJKkHaDeUR8SoiHgiIhZExPyI+OfdrBMR8fWIWBQRz0XE\nhOKUK0ndZPtNno6US5K6QUemRGwBPpFSmhsRg4BnIuLxlNKCHda5ADgm/+d04Jv5pSSVp1ULoFcf\nGHZ01pVIknqAdkfKU0orUkpz868bgYXAzs9MfjvwvZQzGzggInxahKTytWohDDsG+vTLuhJJUg/Q\nqZ7yiBgDnAzM2emrw4GlO7yvY9fgLknlY9V8ONjWFUlS9+hwKI+IgcB/AR9PKTV05WARcU1E1EZE\n7erVq7uyC0kqvuZGWL/EmzwlSd2mQ6E8IvqSC+Q/SCk9sJtVlgGjdng/Mv/Z66SU7kopTUwpTRwx\nYkRX6pWk4lv9Ym7pTZ6SpG7SkdlXArgHWJhSun0Pqz0IvDc/C8skYENKaUUB65Sk7rMqfx+7I+WS\npG7SkdlXzgDeAzwfEfPyn30GGA2QUroTeAS4EFgEbALeX/hSJambrFoIffeHA8ZkXYkkqYdoN5Sn\nlH4LRDvrJOAjhSpKkjK1cj6MOA56+Xw1SVL38F8cSdrZqoX2k0uSupWhXJJ29Fo9vLbKfnJJUrcy\nlEvSjlYtzC0N5ZKkbmQol6Qdtc28YvuKJKn7GMolaUerFsB+B8KgQ7KuRJLUgxjKJWlH22/yjL1O\nOiVJUkEZyiVpu5Tyodx+cklS9zKUS9J2DcugucFQLknqdoZySdpupTd5SpKy0e4TPSWpW6QEa/8M\n27ZmV8Pip3JLR8olSd3MUC6pNDx7H/zsuqyrgCGjc7OvSJLUjQzlkkrDn5+E/YfDhV/Otg5HySVJ\nGTCUSyoNS2bDEZOh+pKsK5Ekqdt5o6ek7DW+Cuv/CqMmZV2JJEmZMJRLyt7SObnlqNOzrUOSpIwY\nyiVlb8kc6NMfDj0p60okScqEoVxS9pbOhsMmQJ9+WVciSVImDOWSsrVlE6x4FkbbuiJJ6rkM5ZKy\ntXwutLZ4k6ckqUczlEvK1pLZueWo07KtQ5KkDBnKJWVr6RwYfizsPzTrSiRJyoyhXFJ2Wltzodx+\ncklSD2col5Sd+hehaYP95JKkHs9QLik72/vJRxvKJUk9m6FcUnaWzoH9h8PQI7OuRJKkTBnKJWVn\nyezcKHlE1pVIkpQpQ7mkbGxcBev+AqO8yVOSJEO5pGzYTy5JUhtDuaRsLJ0Dvavg0JOyrkSSpMwZ\nyiVlY+kcOOxk6FOVdSWSJGWuT9YFSNoHD338b20g5ab+RZjysayrkCSpJBjKpXK1bSv88fsw7BgY\nfkzW1XTeQcfDye/JugpJkkpCu6E8Ir4DXASsSilV7+b7s4CfA3/Jf/RASunWQhYpaTfW/gVaW+DM\nj8NJ78q6GkmStA86MlI+A7gD+N5e1nkqpXRRQSqS1DH1L+WW5ThKLkmSXqfdGz1TSr8B1nZDLZI6\nY3soH2YolySp3BVq9pXJEfFsRDwaEeMKtE9Je1P/Mgw6FPoPzroSSZK0jwpxo+dc4IiU0saIuBD4\nGbDbobuIuAa4BmD06NEFOLTUg9W/ZOuKJEkVYp9HylNKDSmljfnXjwB9I2L4Hta9K6U0MaU0ccSI\nEft6aKnnSikfyo/NuhJJklQA+xzKI+KQiIj869Py+1yzr/uVtBcbV0JzAwx/Q9aVSJKkAujIlIj3\nAWcBwyOiDrgZ6AuQUroTuBS4LiJagM3Au1JKqWgVS3LmFUmSKky7oTyldEU7399BbspESd2lLZQ7\nUi5JUiUo1OwrkrpT/cvQbyAMPizrSiRJUgEYyqVytPrFXOtK7nYOSZJU5gzlUjmqf9nWFUmSKoih\nXCo3zRuhoc6bPCVJqiCGcqncrFmUWzpSLklSxTCUS+Wm/uXc0lAuSVLFMJRL5ab+JYheMPTIrCuR\nJEkFYiiXyk39i3DgWOhTlXUlkiSpQAzlUrlx5hVJkiqOoVwqJ63bcjd6OvOKJEkVxVAulZP1f4Vt\nWxwplySpwhjKpXLizCuSJFUkQ7lUTla/mFvaviJJUkUxlEvlpP4lGDAC9h+adSWSJKmADOVSOXHm\nFUmSKpKhXCon9S/ZuiJJUgUylEvl4rU1sHmtI+WSJFUgQ7lULupfyi0N5ZIkVZw+WRcg9UgpdX6b\n+u0zrxjKJUmqNIZyqbs9eRs8+YWubdtnPxgyqrD1SJKkzBnKpe72wsO50e7qv+/8tgdXQy+7ziRJ\nqjSGcqk7bdkEK+fDmR+Hs27MuhpJklQiHHKTutOKZyFtg8MnZl2JJEkqIYZyqTstq80tRxrKJUnS\n3xjKpe5UVwtDRsPAg7KuRJIklRBDudSdlj0DI0/JugpJklRiDOVSd2lcCRuW2k8uSZJ2YSiXuov9\n5JIkaQ8M5VJ3qauFXn3g0JOyrkSSJJUYQ7nUXZbVwsHjoO9+WVciSZJKjKFc6g6t22DZH+0nlyRJ\nu2Uol7pD/UuwpdF+ckmStFvthvKI+E5ErIqIP+3h+4iIr0fEooh4LiImFL5MqczV5W/ydKRckiTt\nRkdGymcA5+/l+wuAY/J/rgG+ue9lSRVmWS1UDYFhR2ddiSRJKkHthvKU0m+AtXtZ5e3A91LObOCA\niDi0UAVKFaHuGTh8AvSyY0ySJO2qEAnhcGDpDu/r8p9JAtjyGqyabz+5JEnao24dtouIayKiNiJq\nV69e3Z2HlrKzfB6kVvvJJUnSHhUilC8DRu3wfmT+s12klO5KKU1MKU0cMWJEAQ4tlQGf5ClJktpR\niFD+IPDe/Cwsk4ANKaUVBdivVBnqauGAI2DA8KwrkSRJJapPeytExH3AWcDwiKgDbgb6AqSU7gQe\nAS4EFgGbgPcXq1ipLC17BkZPyroKSZJUwtoN5SmlK9r5PgEfKVhF0u6seQU27W0SoBLVvAEaltlP\nLkmS9qrdUC5lrvFVuGNi7mbJcuVIuSRJ2gtDuUrfqgW5QH7u52HEcVlX03lVg3JzlEuSJO2BoVyl\nb80rueX4S2HQIdnWIkmSVAQ+XlClb80i6DcQBh6cdSWSJElFYShX6at/GYYdBRFZVyJJklQUhnKV\nvjWLYNgxWVchSZJUNIZylbaWZli/BIYdnXUlkiRJRWMoV2lb+xcgGcolSVJFM5SrtK15ObccdlS2\ndUiSJBWRoVylbc2i3NKRckl8nVYmAAATLklEQVSSVMEM5SptaxblpkLsPzjrSiRJkorGUK7StuYV\nR8klSVLFM5SrtK1ZZD+5JEmqeH2yLkDao83r4bXVzlEuabd+OGcJr6zemHUZklQQhnKVrjWv5Ja2\nr0jaydK1m/jMT5+nqk8v+vb2l76Syp+hXKXLmVck7cH/fX4FAI9f/yZGD9s/42ok6fXi1s5v4/CC\nSteaRRC94MAxWVciqcQ8OG85NaMOMJBLqhiGcpWuNS/DAUdAn35ZVyKphCxatZEFKxq4+KTDsi5F\nkgrGUK7StWYRDPcmT0mv99Czy4mAi048NOtSJKlgDOUqTSk5R7mkXaSUeOi55Zw+digHD+6fdTmS\nVDCGcpWmxhWwdZNzlEt6nfnLG/jz6tdsXZFUcQzlKk31L+eWjpRL2sFDzy2nT6/ggmpbVyRVFkO5\nSlPbdIj2lEvKSSnx8LMrOPOY4Qwd4A3gkiqLoVylac0r0Hd/GORomKScuUvWsWz9ZqbauiKpAhnK\nVZrWLIKhR0EvL1FJOQ89u4J+fXpxzgkHZ12KJBWciUelac3L3uQpqU3LtlYefm4Fbzn2IAb175t1\nOZJUcIZylZ6WLbDur85RLqnNnL+spX5jM1NrbF2RVJn6ZF2AtIv1f4W0zZlXpAr0rz/7E9+f/dcu\nbTugX2/efOxBBa5IkkqDoVylp23mFUO5VEnWbGzmR39YwpSjhjFxzNBOb18zagj79etdhMokKXuG\ncpWe7aF86JHZ1iGpoGY+U8fWbYnPTh3HMQcPyrocSSop9pSr9NS/DPsPg/07P5ImqTS1tibu+/0S\nTh1zoIFcknajQyPlEXE+8O9Ab+DulNIXd/r+KuDLwLL8R3eklO4uYJ3qrD/cDb/996yr6JrXVsOh\nJ2ZdhaQCevrPa1i8ZhP//FZv4Jak3Wk3lEdEb+A/gHOAOuAPEfFgSmnBTqven1L6aBFqVGdta4Hf\nTId+A2HkqVlX0zXj3pF1BZIK6IdzlnDA/n25oNoHgknS7nRkpPw0YFFK6c8AEfEj4O3AzqFcpeKV\nX0HjCnjnf8LxF2ddjaQebnVjM7+Y/yrvmzKG/n29UVOSdqcjPeWHA0t3eF+X/2xnfx8Rz0XEzIgY\nVZDq1DVzvwcDRsAbzs+6Eknix7VLaWlNvPv00VmXIkklq1A3ej4EjEkpnQg8Dty7u5Ui4pqIqI2I\n2tWrVxfo0Hqdjavgpf+Gk66A3j71TlK2WlsTP/rDEiYdOZSjRgzMuhxJKlkdCeXLgB1Hvkfytxs6\nAUgprUkpNeff3g2csrsdpZTuSilNTClNHDFiRFfqVXuevQ9aW+Dk92RdiSTx1KJ6lq7dzLtPPyLr\nUiSppHUklP8BOCYixkZEP+BdwIM7rhARO965MxVYWLgS1WEpwdzvw6hJMOINWVcjSfxwzl8ZOqAf\n5407OOtSJKmktXujZ0qpJSI+CvyC3JSI30kpzY+IW4HalNKDwD9FxFSgBVgLXFXEmrUnS2bDmpfh\nzOuzrkTK1KsbmnjqZVvksrZ1W+KXC1dx9ZljqerjDZ6StDcdmqc8pfQI8MhOn03b4fVNwE2FLU2d\n9sfvQ79BMO7vsq5EyszmLdu4/FtPs2TtpqxLEdC3d3DFad7gKUnt6VAoVxloaoD5P4Xxl0G/AVlX\nI2Xm33/1MkvWbuJb7zmFEw4dnHU5Pd7Aqj4cOKBf1mVIUskzlFeK+Q/A1k0w4X1ZVyJlZsHyBr79\n1J+5fOJIzht3SNblSJLUYYWaElFZm/s9OOgEOHxC1pVImdjWmrjpgec4cP++fObC47MuR5KkTnGk\nvJSsWgg//yhs29LJDRO8+jyc938goiilSaXue08v5tm6DXz9ipM5YH/bJSRJ5cVQXkp+8+VcMB/7\nvzq/7YjjoObdha9JKgPL1m/my794kbOOHcHFJx7a/gaSJJUYQ3mpWL8U5v8MJl0H530+62qkspFS\nYtrP/kRK8G9vryb8bZEkqQwZykvF77+VW55+bbZ1qEdauKKB933n96x9rbOtU9lL5PrJ//fbjmfU\n0P2zLkeSpC4xlJeC5kZ45t7c/OIHjMq6GvUwzS3buP7+ebSmxIfedGTW5XTJwYP7827nwpYklTFD\neSn4439CcwNM+kjWlagHuv3xl3jh1Ubued9Ezj7eR6FLkpQFp0TMWus2mP0NGD0ZRp6SdTXqYf6w\neC13/ebPvOvUUQZySZIyZCjP2sKHYP0SmOwoubrXxuYWbvjxPEYeuB//+6ITsi5HkqQezfaVrD39\nH3DgGDj2wqwrUQ/zuYcXULduMz/+0GQGVvlXgSRJWfJf4iwt/T3U/R4u+BL06p11Neqkpq3beGll\nY9ZldMnCFQ386A9LufZNR3HqmKFZlyNJUo9nKC+Gli2w+Kn2n8z5+29D1RCoubJ76lLBNDZt5fJv\nzWbhioasS+my4w4ZxPXnHJN1GZIkCUN54bU0w/3vgZd/0bH1z7wBqgYWtyYVVMu2Vj76wz/y0spG\nPvd31Rw6pH/WJXXJqWOHUtXH39BIklQKDOWFtG0r/OT9uUB+7udhzBl7Xz96wUHeYFdOUkpMe3A+\n//PSar54yXje5dzYkiSpAAzlhbKtBf7rA/Di/4ULp8NpH8y6IhXBXb/5Mz+cs4TrzjrKQC5JkgrG\nKRELoXUb/PQaWPBzOO8LBvIK9cjzK/g/j77ARSceyr+ce2zW5UiSpAriSPl2ra3wx+/n5gzvrFef\nz7WsvPWzzje+F0vXbmLaz//Euk1bsy6lSxasaOCUIw5k+mUn0atXZF2OJEmqIIby7eY/AA/9U67P\nm04Grt594eyb4cyPF6W0SrB8/Wbeffds1m/aysmjD8y6nC65sPoQpl08jv59vTlSkiQVlqEccjdo\n/vpzcHA1fOgp6GVXTyG9uqGJK749m/WvbeUHHzydE0cekHVJkiRJJcVQDjD3e7DuL/DuHxvIC2xV\nYxPvvns29Y3NfP9qA7kkSdLuGMq3bIL/uQ1GT4Zjzs26mopSv7GZK789h1c3NHHvP57GhDJtW5Ek\nSSo2Q/mcO2HjSrj8exDevLezlBKX3vk0z/x1XZe279+3F9+96jQf5S5JkrQXPTuUb14Hv/savOEC\nGD0p62pK0twl63nmr+uYetJhjBk+oNPbn33cQZw0ypYVSZKkvenZofy3X4WmBjj7X7OupGT9pHYp\n+/frzRcuGc/Aqp59uUiSJBVLz72rsWE5zPkWnPhOOHhc1tWUpE1bWnjo2eW8bfyhBnJJkqQiKv+k\nNfubsOiXnd9uQ13uSZxv/kzha6oQjzz/Kq9t2cblp47KuhRJkqSKVt6hfMls+O8bYdjR0H9I57bt\nNwAu/BIceERxaqsAP65dytjhA5h4hLOmSJIkFVP5hvLWbfDIJ2Hw4fCh3+RCtgrmL/Wv8fu/rOVT\n5x9LOCuNJElSUZVvKH9mBrz6PFz6HQN5Ecx8Zim9Av5+wsisS5EkSap45Xmj56a18Ot/gyPOhHGX\nZF1NxWnZ1srMZ+o469iDOHhw/6zLkSRJqngdCuURcX5EvBgRiyLixt18XxUR9+e/nxMRYwpd6Os8\n8Xlo2gAX3OYDf4rgqZfrWdnQzOUTvcFTkiSpO7QbyiOiN/AfwAXACcAVEXHCTqt9AFiXUjoa+Cpw\nW6ELbfPq81D7HTj1ajikumiH6cl+XLuUYQP68ZbjDsq6FEmSpB6hIyPlpwGLUkp/TiltAX4EvH2n\ndd4O3Jt/PRM4O4pxd2BK8MinYL8DncqwSNZsbOaXC1fyjpMPp1+f8uxukiRJKjcdudHzcGDpDu/r\ngNP3tE5KqSUiNgDDgPo97bRpxUJe/NzOu2mv2K0c1fIK3xj0Tzz+3QWd2lYds2HzVrZuS1xm64ok\nSVK36dbZVyLiGuAagOMOHUBzn87NmtIMPNr/EmYPuYCB0bsIFWpgVR/eevzBHHvIoKxLkSRJ6jE6\nEsqXATsOm47Mf7a7deoiog8wBFiz845SSncBdwFMnDgxnXjjrztd8InkmtslSZKkStGRpuE/AMdE\nxNiI6Ae8C3hwp3UeBN6Xf30p8OuUUipcmZIkSVLlanekPN8j/lHgF0Bv4DsppfkRcStQm1J6ELgH\n+H5ELALWkgvukiRJkjqgQz3lKaVHgEd2+mzaDq+bgMsKW5okSZLUMzjnnSRJkpQxQ7kkSZKUMUO5\nJEmSlDFDuSRJkpQxQ7kkSZKUMUO5JEmSlDFDuSRJkpQxQ7kkSZKUMUO5JEmSlDFDuSRJkpSxSCll\nc+CIRuDFTA5e+YYD9VkXUcE8v8XjuS0ez23xeG6Lx3NbXJ7f4jk2pTSoMxv0KVYlHfBiSmlihsev\nWBFR67ktHs9v8Xhui8dzWzye2+Lx3BaX57d4IqK2s9vYviJJkiRlzFAuSZIkZSzLUH5XhseudJ7b\n4vL8Fo/ntng8t8XjuS0ez21xeX6Lp9PnNrMbPSVJkiTl2L4iSZIkZSyTUB4R50fEixGxKCJuzKKG\nShURiyPi+YiY15U7f/U3EfGdiFgVEX/a4bOhEfF4RLycXx6YZY3lbA/n95aIWJa/fudFxIVZ1liO\nImJURDwREQsiYn5E/HP+c6/dAtjL+fXa3UcR0T8ifh8Rz+bP7Wfzn4+NiDn5zHB/RPTLutZys5dz\nOyMi/rLDdVuTda3lKiJ6R8QfI+Lh/PtOX7fdHsojojfwH8AFwAnAFRFxQnfXUeHenFKqcZqjfTYD\nOH+nz24EfpVSOgb4Vf69umYGu55fgK/mr9+alNIj3VxTJWgBPpFSOgGYBHwk/3es125h7On8gtfu\nvmoG3pJSOgmoAc6PiEnAbeTO7dHAOuADGdZYrvZ0bgH+ZYfrdl52JZa9fwYW7vC+09dtFiPlpwGL\nUkp/TiltAX4EvD2DOqS9Sin9Bli708dvB+7Nv74X+LtuLaqC7OH8ah+llFaklObmXzeS+0ficLx2\nC2Iv51f7KOVszL/tm/+TgLcAM/Ofe+12wV7OrQogIkYCbwPuzr8PunDdZhHKDweW7vC+Dv9CK6QE\nPBYRz0TENVkXU4EOTimtyL9+FTg4y2Iq1Ecj4rl8e4stFvsgIsYAJwNz8NotuJ3OL3jt7rN8C8A8\nYBXwOPAKsD6l1JJfxczQRTuf25TS9uv28/nr9qsRUZVhieXsa8CngNb8+2F04br1Rs/Kc2ZKaQK5\n9qCPRMT/yrqgSpVyUxc50lBY3wSOIvfr1RXAV7Itp3xFxEDgv4CPp5QadvzOa3ff7eb8eu0WQEpp\nW0qpBhhJ7jfrx2VcUsXY+dxGRDVwE7lzfCowFPh0hiWWpYi4CFiVUnpmX/eVRShfBoza4f3I/Gcq\ngJTSsvxyFfBTcn+pqXBWRsShAPnlqozrqSgppZX5fzhagW/j9dslEdGXXGD8QUrpgfzHXrsFsrvz\n67VbWCml9cATwGTggIjok//KzLCPdji35+fbsVJKqRn4Ll63XXEGMDUiFpNryX4L8O904brNIpT/\nATgmf1dqP+BdwIMZ1FFxImJARAza/ho4F/jT3rdSJz0IvC//+n3AzzOspeJsD41578Drt9PyvYz3\nAAtTSrfv8JXXbgHs6fx67e67iBgREQfkX+8HnEOuZ/8J4NL8al67XbCHc/vCDv+jHuR6nr1uOyml\ndFNKaWRKaQy5TPvrlNKVdOG6zeThQfmpor4G9Aa+k1L6fLcXUYEi4khyo+MAfYAfem67LiLuA84C\nhgMrgZuBnwE/BkYDfwUuTyl5s2IX7OH8nkXu1/8JWAx8aIc+aHVARJwJPAU8z9/6Gz9Dru/Za3cf\n7eX8XoHX7j6JiBPJ3RDXm9yg4Y9TSrfm/237Ebn2ij8C/5Af2VUH7eXc/hoYAQQwD7h2hxtC1UkR\ncRbwyZTSRV25bn2ipyRJkpQxb/SUJEmSMmYolyRJkjJmKJckSZIyZiiXJEmSMmYolyRJkjJmKJek\nEhQRB0TEh/OvD4uImQXa7y0R8cn861sj4q2F2K8kad84JaIklaCIGAM8nFKqLvB+bwE2ppSmF3K/\nkqR940i5JJWmLwJHRcS8iPhJRPwJICKuioifRcTjEbE4Ij4aETdExB8jYnZEDM2vd1RE/HdEPBMR\nT0XEcTsfICJmRMSl+deLI+KzETE3Ip7fvn7+ScHfiYjf54/x9m48B5LUYxjKJak03Qi8klKqAf5l\np++qgUuAU4HPA5tSSicDTwPvza9zF/CxlNIpwCeBb3TgmPUppQnAN/PbAPx/5B4bfRrwZuDLETGg\n6z+WJGl3+mRdgCSp055IKTUCjRGxAXgo//nzwIkRMRCYAvwkIrZvU9WB/T6QXz5DLvQDnAtM3d6H\nDvQHRgML9+1HkCTtyFAuSeWneYfXrTu8byX393ovYH1+lL0r+93G3/59CODvU0ovdrFWSVIH2L4i\nSaWpERjUlQ1TSg3AXyLiMoDIOamLdfwC+Fjkh9wj4uQu7keStBeGckkqQSmlNcDv8jd4frkLu7gS\n+EBEPAvMB7p6g+a/AX2B5yJifv69JKnAnBJRkiRJypgj5ZIkSVLGDOWSJElSxgzlkiRJUsYM5ZIk\nSVLGDOWSJElSxgzlkiRJUsYM5ZIkSVLGDOWSJElSxv4foUL6QJdf/0kAAAAASUVORK5CYII=\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 regression models" ] }, { "cell_type": "code", "execution_count": 42, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from lifelines import AalenAdditiveFitter" ] }, { "cell_type": "code", "execution_count": 43, "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": 43, "metadata": {}, "output_type": "execute_result" } ], "source": [ "aaf = AalenAdditiveFitter()\n", "aaf.fit(data, duration_col='time', event_col='failure', timeline=t)" ] }, { "cell_type": "code", "execution_count": 49, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 49, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuoAAAFpCAYAAADKqw1NAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xu0XHV99/H3N8lJgkm4JSECSUhQ\nVG4x0JOIBjGAQEAIuhZttVbBS4Nalu16Wi3UPg+V1qdYbasoiqnlsdZ6TUWjog1FEG1rJSjKJSgC\nSZNAE0hIQggHck5+zx8zCSeTmXPmsmdmz8z7tdZZZ/Zl9v6dffZKPvM73/37RUoJSZIkSfkypt0N\nkCRJknQgg7okSZKUQwZ1SZIkKYcM6pIkSVIOGdQlSZKkHDKoS5IkSTlkUJckSZJyyKAuSZIk5ZBB\nXZIkScohg7okSZKUQ+Pa3YCRTJs2Lc2ZM6fdzZAkSZIyc9dddz2RUpo+2n65Dupz5sxh9erV7W6G\nJEmSlJmIWFfNfpa+SJIkSTlkUJckSZJyyKAuSZIk5VCua9QlSZLUGrt372bDhg0MDAy0uyldY+LE\nicycOZO+vr663m9QlyRJEhs2bGDKlCnMmTOHiGh3czpeSoktW7awYcMG5s6dW9cxLH2RJEkSAwMD\nTJ061ZCekYhg6tSpDf2FwqAuSZIkAEN6xhq9ngZ1SZIktd22bdv41Kc+ldnxPvaxj7Fr167MjteO\n8xjUJUmS1HaVgvrg4GBdxzOoS5IkSRm48soreeihh5g/fz4LFizg1a9+NUuXLuWEE04A4Atf+AIL\nFy5k/vz5XH755QwNDQHw7ne/m/7+fk488USuvvpqAK677joeffRRzjzzTM4880wAJk+ezPve9z5O\nPPFEXvva1/KTn/yExYsXc+yxx7Jy5UoAhoaGeN/73seCBQuYN28en/nMZwC4/fbbWbx4MZdccgkv\ne9nLePOb30xKqex5shQppcwPmpX+/v60evXqdjdDkiSp661Zs4bjjz8egA9+6z7uf3RHpsc/4aiD\nufqiEytuX7t2LRdeeCH33nsvt99+O6973eu49957mTt3LmvWrOH9738/X//61+nr6+M973kPp512\nGm9961vZunUrhx9+OENDQ5x99tlcd911zJs3jzlz5rB69WqmTZsGFOrFb775Zs4//3ze8IY38PTT\nT/Od73yH+++/n0svvZS7776b5cuXs3nzZv7sz/6MZ599lkWLFvG1r32NdevWcfHFF3Pfffdx1FFH\nsWjRIj7ykY9w+umnH3Ceka7rXhFxV0qpf7Rr5vCMkiRJyp2FCxfuG9bw1ltv5a677mLBggUAPPPM\nMxxxxBEAfPWrX2X58uUMDg7y2GOPcf/99zNv3rwDjjd+/HiWLFkCwMknn8yECRPo6+vj5JNPZu3a\ntQCsWrWKX/ziF6xYsQKA7du38+CDDzJ+/HgWLlzIzJkzAZg/fz5r167l9NNPb+o1MKhLkiRpPyP1\nfLfKpEmT9r1OKXHppZfyV3/1V/vt88gjj/DRj36UO++8k8MOO4zLLrus/HCIAzvo6+sjnnkSgDF7\ndjNhwsGF12PG7KuDTynxiU98gvPOO2+/t99+++1MmDBh3/LYsWPrrp2vhTXqkiRJarspU6bw1FNP\nld129tlns2LFCjZv3gzA1q1bWbduHTt27GDSpEkccsghbNq0ie9+97uFNwzsYMqkF/DU5vWwaysM\n7S6sH9tX+Ep7yp7nvPPO49Of/jS7dxf2/9WvfsXTTz9dd7sbZY+6JEmS2m7q1KksWrSIk046iYMO\nOogZM2bs23bCCSfwl3/5l5x77rns2bOHvr4+rr/+ek477TROmXcSL3vpS5h19FEsOm0BPLsThnaz\n7J1vY8nrf4ujjjyS21Z9p6o2vPOd72Tt2rWceuqppJSYPn063/jGN0Z8z7Jly1iyZAlHHXUUt912\nW0PXoJQPk0qSJKnsQ4+5NLAD9gwrO0kJxo2v7RhDu+EFh2fbrgp8mFSSJEndKYtg3qEM6pIkScqH\n0lAOPRXMSxnUJUmS1B493FteDYO6JEmSWsNgXhODuiRJkrJnGUvDDOqSJElqnL3lmat6wqOIuDEi\nNkfEvcPWfSQiHoiIX0TETRFxaIX3ro2IeyLi7ohwvEVJkqRON7CjMJnQ3q+h3c9PKDS2r66Qvnbt\nOk469RVNaCzc/oMfcuEbfhOAld/5Ltdee21TzpOlWnrUPwd8Evj8sHW3AFellAYj4sPAVcCfVHj/\nmSmlJ+pqpSRJktqny8pYlr7ufJa2aBz1RlTdo55SugPYWrJuVUpp72/tx8DMDNsmSZKkdhitt7zO\nHvNqDA4O8uZL38HxL+/nkje9hV27dnHNh65lwaLXcNKpr2DZe97L3gk7r7v+05wwfwHz+l/JG99y\nGQBPP/00b1/2HhaevphTXnE63/zWgbOSfu6fvsgVV1wBwGWXXcZ73/teXvWqV3HssceyYsWKfft9\n5CMfYcGCBcybN4+rr766KT/vSLKsUX878JUK2xKwKiIS8JmU0vIMzytJkqRGPHgLPHd4IZQD3PpB\nePyBbM9xxIlwzgdH3e2Xv3qQf7jheha96jTevuw9fOozn+WKdy/j/3zgSgDe8rbf49s3f4+LXnc+\n137073jkgXuYMGEC27ZtA+BD136UsxafwY3LP8W2bdtYePqZvPasxSOe87HHHuNHP/oRDzzwAEuX\nLuWSSy5h1apVPPjgg/zkJz8hpcTSpUu54447OOOMMxq9ElXLJKhHxAeAQeCfK+xyekppY0QcAdwS\nEQ8Ue+jLHWsZsAxg9uzZWTRPkiRJez14Cwxs33/d7gHom1boKQeIsa1vV9GsmTNZ9KrTAPjdN/02\n133qBubOOYa//tuPsWvXM2x98klOPOF4Lnrd+cw76UTefNk7ef1Fr+P1Sy8EYNWt32fld27mox/7\nBAADA8/y3+s3jHjO17/+9YwZM4YTTjiBTZs2FY6zahWrVq3ilFNOAWDnzp08+OCDnRXUI+Iy4ELg\n7LT37xAlUkobi983R8RNwEKgbFAv9rYvB+jv7y97PEmSJFWpNJjvHoDDjjlwv53DXlfR890sEXHA\n8nv+4H+x+t9/wKxZM/nzv/i/DAwMAPCdb6zgjh/+O9+6+bt86MMf5Z67fkxKiX/58hd46UuO2+84\nmzZtrnjOCRMm7Hu9N86mlLjqqqu4/PLLs/rRalZ1jXo5EbEEeD+wNKW0q8I+kyJiyt7XwLnAveX2\nlSRJUoMevAXuWfH811ObYPKM57/KhfQc+e/16/nPH/8XAF/8ytc4/VWvBGDatKns3LmTFTd9E4A9\ne/awfv0Gzlx8Bh/+0DVs37GDnTt3ct5rz+YTn7phX+D+2d0/r6sd5513HjfeeCM7dxY+wWzcuJHN\nmyuH/Waoukc9Ir4ELAamRcQG4GoKo7xMoFDOAvDjlNK7IuIo4LMppQuAGcBNxe3jgC+mlL6X6U8h\nSZLUiyqVseQ8jI/kpS85jutv+Hvefvnvc8LxL+Pdy97Bk08+yUmnvoIXzpjBgt84FYChoSF+922/\nx/YdO0gp8d73vItDDz2U//2n7+cP//hK5vW/kj179jB3zjF8+6av1dyOc889lzVr1vDKVxY+KEye\nPJkvfOELHHHEEZn+vCOJCtUqudDf359Wr3bYdUmSJKD6MpY6rNk5meNf8uJMjpV7Q7uhRcMzrlmz\nhuOPP36/dRFxV0qpf7T3OjOpJElSXjUxmCv/DOqSJEl50IVlLGqMQV2SJKkd7C3XKAzqkiRJ7TCw\nvTAKS46klA4YHlH1a/RZUIO6JElSs1Uqa8mRiWP2sOXJ7Uw97BDDegZSSmzZsoWJEyfWfQyDuiRJ\nUtY6sKxl5sRdbNi6mcefeKLdTWm+PXtg/Kamn2bixInMnDmz7vcb1CVJkhrVgcG8VN8YmPuCsvNX\ndp+dm+D4S9rdilEZ1CVJkmrh6CxqEYO6JEnSSLqgt1ydyaAuSZI0nMFcOWFQlyRJvcsyFuWYQV2S\nJPUOe8vVQQzqkiSpd+RwkiGpEoO6JEnqXuV60KUOYVCXJEndwXpzdRmDuiRJ6kzWm6vLGdQlSVJn\nMJirxxjUJUlS/ljGIhnUJUlSDjk6i2RQlyRJOeDoLNIBDOqSJKn1rDeXRmVQlyRJzWW9uVQXg7ok\nScqWveVSJgzqkiQpWz4IKmXCoC5Jkhrjg6BSUxjUJUlS9aw3l1rGoC5Jkiqz3lxqG4O6JEmqzHpz\nqW1qCuoRcSNwIbA5pXRScd3hwFeAOcBa4LdSSk+Wee+lwJ8VF/8ypfSP9TdbkiQ1hfXmUm7U2qP+\nOeCTwOeHrbsSuDWldG1EXFlc/pPhbyqG+auBfiABd0XEynKBXpIktYj15lKu1RTUU0p3RMScktUX\nA4uLr/8RuJ2SoA6cB9ySUtoKEBG3AEuAL9XUWkmSVD/rzaWOkkWN+oyU0mPF1/8DlCtkOxpYP2x5\nQ3GdJElqFevNpY6S6cOkKaUUEamRY0TEMmAZwOzZszNplyRJPcl6c6mjZRHUN0XEkSmlxyLiSGBz\nmX028nx5DMBMCiUyB0gpLQeWA/T39zcU+iVJ6hnWm0tdZ0wGx1gJXFp8fSnwzTL7/CtwbkQcFhGH\nAecW10mSpCzsLWsZ/mVIlzparcMzfolCz/i0iNhAYSSXa4GvRsQ7gHXAbxX37QfelVJ6Z0ppa0T8\nBXBn8VDX7H2wVJIk1cGyFqnr1Trqy5sqbDq7zL6rgXcOW74RuLGm1kmSpAJHbJF6jjOTSpLUCRyx\nReo5BnVJkvKm0oOhknqKQV2SpHazrEVSGQZ1SZLazbIWSWUY1CVJajVHbJFUBYO6JEnN5EREkupk\nUJckqZksa5FUJ4O6JElZsqxFUkYM6pIkZckedEkZMahLktQIe9AlNYlBXZKkavlgqKQWMqhLklQt\ny1oktZBBXZKkSixrkdRGBnVJkiqxB11SGxnUJUmCyvXnktQmBnVJUm8qV9biQ6GScsSgLknqTZa1\nSMo5g7okqTf4YKikDmNQlyT1BnvQJXUYg7okqfv4YKikLmBQlyR1H3vPJXUBg7okqfNZfy6pCxnU\nJUmdzx50SV3IoC5J6jz2oEvqAQZ1SVLnsQddUg8wqEuS8s0RXCT1KIO6JCnf7D2X1KPGtLsBkiRJ\nkg5kUJckSZJyqOGgHhEvjYi7h33tiIg/LNlncURsH7bP/2n0vJIkSVI3a7hGPaX0S2A+QESMBTYC\nN5XZ9YcppQsbPZ8kSZLUC7IufTkbeCiltC7j40qSJEk9Jeug/kbgSxW2vTIifh4R342IEzM+ryRJ\nktRVMgvqETEeWAp8rczmnwLHpJReDnwC+MYIx1kWEasjYvXjjz+eVfMkSZKkjpJlj/r5wE9TSptK\nN6SUdqSUdhZf3wz0RcS0cgdJKS1PKfWnlPqnT5+eYfMkSZKkzpHlhEdvokLZS0S8ENiUUkoRsZDC\nB4QtGZ5bktQtSmcidRZSST0qk6AeEZOAc4DLh617F0BK6QbgEuDdETEIPAO8MaWUsji3JKnDlQvm\nhx3TvvZIUk5kEtRTSk8DU0vW3TDs9SeBT2ZxLklSlxnYDpNntLsVkpQ7WZa+SJI0stLec7C0RZIq\nMKhLkprHshZJqptBXZLUPJa1SFLdsp7wSJIkSVIGDOqSJElSDhnUJUmSpBwyqEuSJEk5ZFCXJEmS\ncshRXyRJ2Sk3HKMkqS4GdUlSfSpNXuQ46ZKUCYO6JKk6Tl4kSS1lUJcklWcwl6S2MqhLkixjkaQc\nMqhLUi+yt1yScs+gLkm9wGAuSR3HoC5J3cYyFknqCgZ1Sep09pZLUlcyqEtSpzGYS1JPMKhLUp5Z\nxiJJPcugLkl5Ym+5JKnIoC5J7WQwlyRVYFCXpFaxjEWSVAODuiQ1i73lkqQGGNQlKSsGc0lShgzq\nklQPy1gkSU1mUJekegxsh8kz2t0KSVIXM6hLUjXKlbVIktREBnVJqoY96JKkFhvT7gZIkiRJOlBm\nQT0i1kbEPRFxd0SsLrM9IuK6iPh1RPwiIk7N6tySJElSt8m69OXMlNITFbadDxxX/HoF8Onid0mS\nJEklWln6cjHw+VTwY+DQiDiyheeXJEmSOkaWPeoJWBURCfhMSml5yfajgfXDljcU1z2WYRsk9bpy\n45tnwVFeJEktlmVQPz2ltDEijgBuiYgHUkp31HqQiFgGLAOYPXt2hs2T1JWcDVSS1KUyC+oppY3F\n75sj4iZgITA8qG8EZg1bnllcV3qc5cBygP7+/pRV+yR1AWcDlST1kEyCekRMAsaklJ4qvj4XuKZk\nt5XAFRHxZQoPkW5PKVn2Iqkye8slST0sqx71GcBNEbH3mF9MKX0vIt4FkFK6AbgZuAD4NbALeFtG\n55bUrZxkSJLUwzIJ6imlh4GXl1l/w7DXCfj9LM4nSZIkdTtnJpUkSZJyyKAuSZIk5VDWM5NK6nTN\nGoe8Ho5dLknqYQZ1qdc5sookSblkUJd6jcFckqSOYFCXupkTBEmS1LEM6lI3cxxySZI6lqO+SJIk\nSTlkUJckSZJyyKAuSZIk5ZBBXZIkScohg7okSZKUQwZ1SZIkKYcM6pIkSVIOOY661MnKTWg03O6B\n1rVFkiRlyqAudQpnGZUkqacY1KW8Kg3mhnL1ol99Dx6+vd2tkNTpDj8WFi5rdytqZlCX8mpgO0ye\n0e5WSO318O2w9RE4fG67WyJJLWdQlyTl2+FzYcm17W6FJLWcQV2SJEld7Wfrt7Hr2cF9y4fu2cGJ\nJ7exQVUyqEuSJKljlYbwcp4b2sMRUybuW35m61Czm5UJg7okSZJya7QgXhrCu4lBXWqG0hFbJh4C\nx53TvvZIkpQD1fR+l+rmID4ag7rUqGrGN39yHdyzorbjOlmRJKmD1FOCopEZ1KVGVTOMouOfqxdl\nMQa6QzNKuWAIbw+DulSrchMRSTpQFmOgHz4Xjl2cUYMkVdLLdeB5ZlCXRuMMoVL9HANdajt7wzuX\nQV0arpp6c0mScsIQ3t0aDuoRMQv4PDADSMDylNLHS/ZZDHwTeKS46usppWsaPbfUMHvLJUk5ZklK\nb8uiR30Q+KOU0k8jYgpwV0TcklK6v2S/H6aULszgfFJ2qnkQVJKkJrA3XKNpOKinlB4DHiu+fioi\n1gBHA6VBXWo/HwRVPbIYvaQXOWKLepghXFnItEY9IuYApwD/VWbzKyPi58CjwB+nlO7L8tzSAaw3\nV1ayGL2kFzlii7qUIVytkllQj4jJwL8Af5hS2lGy+afAMSmlnRFxAfAN4LgKx1kGLAOYPXt2Vs1T\nL7KsRVly9BKpZ1gXrrzIJKhHRB+FkP7PKaWvl24fHtxTSjdHxKciYlpK6Yky+y4HlgP09/enLNon\nSZIE9oars2Qx6ksA/wCsSSn9bYV9XghsSimliFgIjAG2NHpuSV0sq7rwYxfDS5Y0fhxJuWcIV7fJ\nokd9EfAW4J6IuLu47k+B2QAppRuAS4B3R8Qg8AzwxpSSveWSKsuiLnzTvYUvp7GXuoIlKeo1WYz6\n8iMgRtnnk8AnGz2XpB7TaF14Vr3yPhQptVy5UG4QV69xZlJJ3eslSyx7kXLK3nFpdAZ1qVtlOfb3\nsYsNvJKqZq24lA2DurJVbuzyiYfAcee0pz29LKuxv7Oq866VdeFSxygN5oZwKRsGdTWm3EyfpRMK\n7dxU2zHKMezXJ4uxv9s1K6d14VIu2DsutY9BXZVVE6Crmelz9wDcs6KxYzy5buRjVDquGmedt9S1\nDOFSvhnUVVlWM3uOFsJbdQxJ6nGWqEidxaAuSVIXsHdc6j4GdUmSOpC941L3M6hLeZTFA5yOmiJ1\nLHvHJYFBvbeN9rCoD2O2TxZDKzpqitQx7B2XVI5BvVeUC+XVjLai9sliaEVJbWfvuKR6GdS7VTXj\nm0uSMmfvuKSsGNS7VVZDK0qSRmQwl9QsBnV1p3bNplnJsYudNEjqQJatSGong7q6UxYPY2Zl072F\nr1o+OOSl7VIPKRfKDeGS2smg3olGG60FHLEF8vMwZj29+47YIjWdJSuS8s6g3gl8MLSzvWSJZS9S\nDhjMJXUag3on8MFQSaqJZSySuoFBvd0sY5Gk/ZSG7BdMGMcpsw4dcZ9ShnJJ3cCg3mqdUsaS5agp\nxy629EPqQtWMiFIasusZRWXzUwP8+6+fGHEfSepGBvVm65RgXiqrUVPqGfGknGMXG/alFsui17o0\nZNcTsA3kknqVQb3ZOrm+PItRU7LomXd4Q6lh1fRkl8qi19qQLUn1M6h3oyzCcVZBN4sRTxzeUBqR\nk/JIUncyqHejLMpW8hR0Hd5QPcwQLkm9y6CepXIjuLRrxJa8TPYjaUSOXiJJqsSg3ohWPChaT9mH\n9dm5U099cLkh6VS7eq59OfX8PuwNlyQ1wqBei3aM4FJPGUueylZ6VBYzIG7b9VzD580q7GcVdmtV\nTftb1SNd+vswhEuSms2gXkmlMpZagnmWD3X2eBlLnnqkWxXQnhvac8DY0dW8Z/h5swj75Y7bKuXG\nzy7VqraV/j4M4ZKkZjOoV5LFsIrd9lBnkzQr+GYVUrNoSz2adY5OqonOSzsgX22RJPWGTIJ6RCwB\nPg6MBT6bUrq2ZPsE4PPAbwBbgN9OKa3N4txN8d0r4aHvw9i+xo5jbziQ72CYRYlKnpXrle+2n1GS\npG7VcFCPiLHA9cA5wAbgzohYmVK6f9hu7wCeTCm9OCLeCHwY+O1Gz517PdAbXo1dzw5y6AvGt/y8\n1ZSOdHto7eafTZKkbpdFj/pC4NcppYcBIuLLwMXA8KB+MfDnxdcrgE9GRKSUUgbnz97518I9K0Yt\nfammTOIFE8ZxSo3vyUJpfXazaryrLRVpB0OqJEnqZFkE9aOB9cOWNwCvqLRPSmkwIrYDU4HanpRr\ns3rKJEofhmtVD24W583Tg3ySJEm9JncPk0bEMmAZwOzZs9vShg9+6z5W3zeWQbbtt35PgvHjaivh\neM3REzm/ySM4lpNFeO70AD5h53r2jN3/ZxgzNMCzk2e1qUWSJEnVyyKobwSGJ5+ZxXXl9tkQEeOA\nQyg8VHqAlNJyYDlAf39/20pjBtMY+qKkpCP2kKg+vN6zpfD1g9KrUaPXHE1bwn6n2zN2IlvmXrTf\nuqmPfKtNrZEkSapNFkH9TuC4iJhLIZC/Efidkn1WApcC/wlcAnw/t/XpwNUXncjKWYcyvaRH+ZCN\nP2Dcc4/vWx6td/a76xoP6fWG/bOnbeWCWUP7luvpSe7GHunB8QfT90z1v0NJkqR2aTioF2vOrwD+\nlcLwjDemlO6LiGuA1SmllcA/AP8UEb8GtlII87k25aA+Hn9qYL9161/Qz6yjJ+1bPmTjD0YMfecf\nA6+femDYHa4ZYf/h7RBDk3nFGWdUbGs15y7XI13uOLUcsx71fuA5a/pETisZwn770a/Zb3m032Gn\n68YPW5Ik9YrIccc2/f39afXq1e1uxj63/XIzTz2zu+L2Qx+9nReO33/74PiDDwiHwxV66Xfst640\nSJULW6WGv+fK/4CHt+9h9vRDRnzPuGe3E3uGl/fsIQ07z1nTn+K0V58z4jFKlf4839kwkdv+p2+/\n49bqnmKR1MlTq3/Pw9vhRZOf46rfPGP0nYeZ+si32H3Q9Jrek2d9zzxe9sPW8N9RnoN76b3fzrbm\nqS2SpMY8s/VR+i/6vbadPyLuSin1j7Zf7h4mzbMzX3rEiNtvO+gcHigJ8gODQ4z0X3m5EF/ay/vc\nQUeMGPZL33PW9ImQJjHaYIyDE/YP8uOe3U4MFWbz/MW28dyzZQqrtt83ylFKTSt+Fax57CmgtpBd\n6uSptdfpX/kfwNCoux2gtDQGui+Qld5Lea7bL/2rTj1/GerGtkiSeoNBPUPlgvxtv9y8XwnNwOAQ\nsw6bdMB+w40Wykd7z2lz4bSaj7C/W9ds4t8fanz0zOOPnMJrp27lohduG33nEYwZGuDZET/yZKOa\nD061hrF6S3eyeM6gGp1Ut9+M3083tEWS1J0M6k1WGt5Lg3s1qgn3WTv7+BmcffzIEz7VouwQPzUY\nrTa+VAwdwkM7+7jm27X+RaCc/f9KUCgZem7Y9j0jlvbUW7qTxXMG1Ritbr9e7QrM7QzLeWqLJKnz\nGdRbbLTymXLq6ZXvNrX+lWHhwCZ2Z/AXgXJGKhl63vPhvd7SnYe2jyv5oLH/B4Z95x7hQ0O5B2pH\nU89fdMrJIvAPjj+45vc064NHs9pieJckVWJQ7wD19Mr3YpgfLuu/CNSq7FCeNZTuvOZoiKFBKj+6\nXDDSh4aHdo6DNKnhMqh6ZRX4G5WXdoDlMpKk2hjUO1A1vfIr725wAHc1pNFe3aVHwAUzD2b70SfW\n3YZrvn3fqA8Uq/066eFeSVJrGdS7VLlx4EfT673wzZSnXt1aZfVgcT0WvWhaW/8y0g6d9HCvJKm5\nDOpdylp4AazbsqvhB2r3DrF5/JFTsmhS1dZt2QU80VBQz/JDRqs+NFjXLknay6CufUrDfVblM+uf\nfJqJ48buW67nA0DpMcrxg8X+Fr1oGpDNEJvt6Nm+5tv3NfxBI6sPGVl8aKiXde2S1LsM6qqomvKZ\nasLxxHFjWTr/6H3L9QxROX3KxNEnnPIvAvtp9wO1jcrig0ZWHzKy+NCQrcaGDG2lWkc8kiQ9z6Cu\niqopn6kmdE85qK/m49ajWX8RUHvk6YNGVn+daJZahwxtlYe3F743GtQn7FzPnpK2+1cESb3AoK6G\nNCt0S3mSpw8N9So7ZGiTg+6V/1EI61f+R2PHOWv6oZz26nP2W+foOJJ6gUFdknpAO2ZNLcwHMABD\nY4atra1nvzBL74QD5gNwdBxJvcCgrq5Vrsa+1+vWpb1aMbrM+cfA0iOeYsvci/Y7z/Ce/dHOU36W\nXiit0+97Zgtp5OfNJfWwYw+GZSe1uxW1M6ira5Uryxmtpt4gr17VqtFlaj1PtbP0pjHjSury8/NA\nrSTVy6CunlLryDFgeFfvqmeG3XrC/EizsxZ65bezZe4ZNR2z1p57Scojg7o0TLN64R0HXt2gmhl2\nS8P84PiDaz5Paf15Pcco11YfQJXUaQzq0ijq6YUv5Tjw6hXVhPlWHKMcH0CV1GkM6lKDshqi0nHg\npeYaqcRGkvLIoC5J6kn2sEsIkO+iAAAMbUlEQVTKO4O6lFPlhpccjeUyUvXsYZeUdwZ1KafqKamx\nzl2qX2kPO9jLLqm9DOpSF7HOXaqfI8VIyhuDuiRJFVjHLqmdDOqSJFVQzaRPhndJzWJQl7pYNQ+k\n1jNhk7Xv6lWWx0hqJYO61MWqeSC1ngmbqnlPFvxAoE5geYykZjGoSz2untFlsprkaTQ+DKtO4DCP\nkprFoC4pt0pLd+xhVyewh11SVhoK6hHxEeAi4DngIeBtKaVtZfZbCzwFDAGDKaX+Rs4rqTeU9txX\nU3JTGuZL6+vr5YcEVcsHUCVlpdEe9VuAq1JKgxHxYeAq4E8q7HtmSumJBs8nqYfVU3NfWl9fL8tw\nVC8fQJVUr4aCekpp1bDFHwOXNNYcSWpMq+rnpUZYHiOpGlnWqL8d+EqFbQlYFREJ+ExKaXmG55Wk\nprNeXlkarTzG4C4JqgjqEfFvwAvLbPpASumbxX0+AAwC/1zhMKenlDZGxBHALRHxQErpjgrnWwYs\nA5g9e3YVP4IkNV8W9fJSJY4cI6mcUYN6Sum1I22PiMuAC4GzU0qpwjE2Fr9vjoibgIVA2aBe7G1f\nDtDf31/2eJLUblmNUW+YVzmlpTFgL7vUixod9WUJ8H7gNSmlXRX2mQSMSSk9VXx9LnBNI+eVpE6Q\nRZg3yPemcg+gWh4j9Z5Ga9Q/CUygUM4C8OOU0rsi4ijgsymlC4AZwE3F7eOAL6aUvtfgeSWpK4wW\n5uuZBdZw350sj5F6T6Ojvry4wvpHgQuKrx8GXt7IeSSpV9Uzio0lN73BkWOk7ufMpJLUZayf7w1O\nrCR1P4O6JPWgrMJ8KcN9+zixktR9DOqSpLKyKLsxuLeX5TFSZzOoS5IyUxruV969sU0tETixktTp\nDOqSpKZxRtd8sa5d6iwGdUlS0zija75Z1y7lm0FdktQy9TzEanBvLevapfwwqEuScsU69/ayrl3K\nD4O6JEmqyLp2qX0M6pKkXCt9IBUsh2kn69ql1jGoS5JyrVxdu3Xs+WJdu9QcBnVJUscZbTQZg3tr\nWdcuNYdBXZLU8aoZBtLw3jrWtUvZMKhLkrpONeUypQzyzWNdu1Qfg7okqSeMNoa7vfCtZV27NDqD\nuiRJ+NBqq1nXLo3OoC5JUgXWvreOde3SgQzqkiRVydr31ilX126vu3qNQV2SpAZY+946lsuo1xjU\nJUlqInvhm6c0uDuSjLqNQV2SpBazF745SkeSAXvZ1dkM6pIk5Ywj0NTHunZ1G4O6JEkdoJoRaEoZ\n5q1rV2czqEuS1IFGK58Ba+HLcRhIdRKDuiRJXcpa+NFZLqM8M6hLktSjHJGmPMtllBcGdUmStI+9\n8AdyGEi1i0FdkiRVzV54h4FU6zQU1CPiz4HfA/berX+aUrq5zH5LgI8DY4HPppSubeS8kiQpP+rp\nhS/VSWHeuna1ShY96n+XUvpopY0RMRa4HjgH2ADcGRErU0r3Z3BuSZKUc1mMUAP5DvPWtasZWlH6\nshD4dUrpYYCI+DJwMWBQlyRJQH1hvpODOxjeNbosgvoVEfFWYDXwRymlJ0u2Hw2sH7a8AXhFBueV\nJEk9pJMnfbJcRvUYNahHxL8BLyyz6QPAp4G/AFLx+98Ab2+kQRGxDFgGMHv27EYOJUmSulinT/pk\nr7tGM2pQTym9tpoDRcTfA98us2kjMPwOm1lcV+l8y4HlAP39/amac0uSJJXTSQ+62uuuUo2O+nJk\nSumx4uIbgHvL7HYncFxEzKUQ0N8I/E4j55UkScpC3h909SHV3tZojfpfR8R8CqUva4HLASLiKArD\nMF6QUhqMiCuAf6UwPOONKaX7GjyvJElSS+SpxMZymd7SUFBPKb2lwvpHgQuGLd8MHDC+uiRJUjdo\nV4mN5TLdzZlJJUmSmqyVJTaWy3QPg7okSVIONCvMWy7TuQzqkiRJHSKLML/+Bf3MOnr/Xnl73fPJ\noC5JktRF6qmXf/zg/eeiPPTR23mhwb3tDOqSJEk9pKpe+YPO4YFndu9bnv4/d9C39dH99tmzexeT\nZrw48/bpeQZ1SZIk7efAMP+mA/b5+W1f45lh4X3s0AC7Js3ct/zc0B6OmDKxWU3sCQZ1SZIk1ezl\nZ/7m/isevAUGtu9bXLNhM1t3HTniMQzzIzOoS5IkqXHHnbPf4vET9w/uAOwegMOO2bf4s/Xb2Lbr\nuZpO00vh3qAuSZKk7JUEd6DQ675z077FUyYPwKxjDtxvBNWE+24J8wZ1SZIktUZpeC8J7qU97uWc\nMuvQUU8zWpg/dPzYUY+RBwZ1SZIktcdowR2qCu+lRg3zO4dqOl67GNQlSZKUD1WUy9QT3DuVQV2S\nJEn5lUG5TKcyqEuSJKlzNKlcJo8M6pIkSepcXVwuY1CXJElSdxmt133iIa1tT50M6pIkSepu5Xrd\nO8CYdjdAkiRJ0oEM6pIkSVIOGdQlSZKkHDKoS5IkSTlkUJckSZJyyKAuSZIk5ZBBXZIkScohg7ok\nSZKUQwZ1SZIkKYcM6pIkSVIOGdQlSZKkHDKoS5IkSTlkUJckSZJyKFJK7W5DRRHxOLCujrdOA57I\nuDm9xmvYOK9h47yGjfH6Nc5r2DivYeO8ho3L2zU8JqU0fbSdch3U6xURq1NK/e1uRyfzGjbOa9g4\nr2FjvH6N8xo2zmvYOK9h4zr1Glr6IkmSJOWQQV2SJEnKoW4N6svb3YAu4DVsnNewcV7Dxnj9Guc1\nbJzXsHFew8Z15DXsyhp1SZIkqdN1a4+6JEmS1NG6KqhHxJKI+GVE/Doirmx3ezpRRKyNiHsi4u6I\nWN3u9nSCiLgxIjZHxL3D1h0eEbdExIPF74e1s415V+Ea/nlEbCzei3dHxAXtbGPeRcSsiLgtIu6P\niPsi4g+K670XqzTCNfRerFJETIyIn0TEz4vX8IPF9XMj4r+K/z9/JSLGt7uteTXCNfxcRDwy7D6c\n3+625llEjI2In0XEt4vLHXkPdk1Qj4ixwPXA+cAJwJsi4oT2tqpjnZlSmt+Jwxi1yeeAJSXrrgRu\nTSkdB9xaXFZln+PAawjwd8V7cX5K6eYWt6nTDAJ/lFI6ATgN+P3iv4Hei9WrdA3Be7FazwJnpZRe\nDswHlkTEacCHKVzDFwNPAu9oYxvzrtI1BHjfsPvw7vY1sSP8AbBm2HJH3oNdE9SBhcCvU0oPp5Se\nA74MXNzmNqkHpJTuALaWrL4Y+Mfi638EXt/SRnWYCtdQNUgpPZZS+mnx9VMU/oM6Gu/Fqo1wDVWl\nVLCzuNhX/ErAWcCK4nrvwxGMcA1VpYiYCbwO+GxxOejQe7CbgvrRwPphyxvwH9h6JGBVRNwVEcva\n3ZgONiOl9Fjx9f8AM9rZmA52RUT8olgaY8lGlSJiDnAK8F94L9al5BqC92LViiUHdwObgVuAh4Bt\nKaXB4i7+/zyK0muYUtp7H36oeB/+XURMaGMT8+5jwPuBPcXlqXToPdhNQV3ZOD2ldCqFEqLfj4gz\n2t2gTpcKQyvZG1K7TwMvovCn38eAv2lvczpDREwG/gX4w5TSjuHbvBerU+Yaei/WIKU0lFKaD8yk\n8Nful7W5SR2n9BpGxEnAVRSu5QLgcOBP2tjE3IqIC4HNKaW72t2WLHRTUN8IzBq2PLO4TjVIKW0s\nft8M3EThH1nVblNEHAlQ/L65ze3pOCmlTcX/rPYAf4/34qgioo9CwPznlNLXi6u9F2tQ7hp6L9Yn\npbQNuA14JXBoRIwrbvL/5yoNu4ZLiqVZKaX0LPD/8D6sZBGwNCLWUiiDPgv4OB16D3ZTUL8TOK74\nVO944I3Ayja3qaNExKSImLL3NXAucO/I71IFK4FLi68vBb7ZxrZ0pL3hsugNeC+OqFiD+Q/AmpTS\n3w7b5L1YpUrX0HuxehExPSIOLb4+CDiHQq3/bcAlxd28D0dQ4Ro+MOwDd1Cor/Y+LCOldFVKaWZK\naQ6FLPj9lNKb6dB7sKsmPCoOmfUxYCxwY0rpQ21uUkeJiGMp9KIDjAO+6DUcXUR8CVgMTAM2AVcD\n3wC+CswG1gG/lVLyYckKKlzDxRRKDRKwFrh8WK21SkTE6cAPgXt4vi7zTynUWHsvVmGEa/gmvBer\nEhHzKDyoN5ZCZ+BXU0rXFP9/+TKFko2fAb9b7BlWiRGu4feB6UAAdwPvGvbQqcqIiMXAH6eULuzU\ne7CrgrokSZLULbqp9EWSJEnqGgZ1SZIkKYcM6pIkSVIOGdQlSZKkHDKoS5IkSTlkUJckSZJyyKAu\nSZIk5ZBBXZIkScqh/w9rNnYot73EhAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "aaf.plot()" ] }, { "cell_type": "code", "execution_count": 45, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from lifelines import CoxPHFitter" ] }, { "cell_type": "code", "execution_count": 46, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 46, "metadata": {}, "output_type": "execute_result" } ], "source": [ "cp = CoxPHFitter()\n", "cp.fit(data, duration_col='time', event_col='failure')" ] }, { "cell_type": "code", "execution_count": 47, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "n=42, number of events=30\n", "\n", " coef exp(coef) se(coef) z p lower 0.95 upper 0.95 \n", "treatment -1.5092 0.2211 0.4096 -3.6849 0.0002 -2.3121 -0.7063 ***\n", "---\n", "Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 \n", "\n", "Concordance = 0.690\n" ] } ], "source": [ "cp.print_summary()" ] }, { "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.11" } }, "nbformat": 4, "nbformat_minor": 0 }