deeso/survivor-analysis.ipynb

## survivor-analysis.ipynb
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "597b7325-20cb-4bff-9637-49376ac0cb2c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "count    21354.000000\n",
       "mean         0.108364\n",
       "std          0.310847\n",
       "min          0.000000\n",
       "25%          0.000000\n",
       "50%          0.000000\n",
       "75%          0.000000\n",
       "max          1.000000\n",
       "Name: status, dtype: float64"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEPCAYAAACjjWTcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAhNUlEQVR4nO3de5RdZZnn8e+PJIRLECJUMCTRIAYVUAKdhTrMOC1IE7Ht4Gpt4+pBdKRjK462Y48NLtcanWVap0fUHm2dhhaNSkvHC0280WIapFEkFhAgFzAhAVIkJAUkJJh71TN/vO/ps6lU1bnUOalTtX+ftc566uyz93veffbez373uy+liMDMzMa3I0a7AmZm1n5O9mZmJeBkb2ZWAk72ZmYl4GRvZlYCE0e7AgAnnXRSzJ49e7SrYWY2ptxzzz1PRURXPeN2RLKfPXs23d3do10NM7MxRdJj9Y7rbhwzsxJwsjczKwEnezOzEnCyNzMrASd7M7MScLI3MysBJ3szsxJwsjczKwEnezOzEuiYZD/7qh+PdhXMzMatjkn2ZmbWPk72ZmYl4GRvZlYCTvZmZiVQM9lLOkrSCkn3S1ot6VN5+CclPSFpZX5dUpjmaknrJT0s6eJ2zoCZmdVWz/Ps9wEXRMRzkiYBd0r6af7sCxHxueLIks4AFgJnAqcAP5d0ekT0tbLiZmZWv5ot+0iey28n5VcMM8kC4MaI2BcRG4H1wHkjrqmZmTWtrj57SRMkrQS2AbdGxN35ow9KekDS9ZKm5mEzgE2FyXvysIFlLpLULam7t7e3+TkwM7Oa6kr2EdEXEXOBmcB5ks4CvgqcBswFtgDX5NE1WBGDlHltRMyLiHldXXX9C0UzM2tSQ1fjRMQO4HZgfkRszTuBfuA6ql01PcCswmQzgc0jr6qZmTWrnqtxuiSdkP8+Gngj8JCk6YXR3gqsyn8vAxZKmizpVGAOsKKltTYzs4bUczXOdGCJpAmkncPSiPiRpG9JmkvqonkUeB9ARKyWtBRYAxwErvSVOGZmo6tmso+IB4BzBhl+2TDTLAYWj6xqZmbWKr6D1sysBJzszcxKwMnezKwEnOzNzErAyd7MrASc7Mehev/Fo/8VpFl5ONmbmZWAk72ZWQk42ZuZlYCTvZlZCTjZm5mVgJO9mVkJONmXwMBLLH3JpVn5ONmbmZWAk72ZWQk42ZuZlYCTvdWt0ccw+NyAWedwsrdDOEmbjT9O9mZmJVAz2Us6StIKSfdLWi3pU3n4CyXdKmldjlML01wtab2khyVd3M4ZMDOz2upp2e8DLoiIs4G5wHxJrwWuApZHxBxgeX6PpDOAhcCZwHzgK5ImtKHuZmZWp5rJPpLn8ttJ+RXAAmBJHr4EuDT/vQC4MSL2RcRGYD1wXisrbZ3Jff1mnauuPntJEyStBLYBt0bE3cDJEbEFIMdpefQZwKbC5D152MAyF0nqltTd29s7glkYH4ZKlIczgTpZm41fdSX7iOiLiLnATOA8SWcNM7oGK2KQMq+NiHkRMa+rq6uuypoTspk1p6GrcSJiB3A7qS9+q6TpADluy6P1ALMKk80ENo+0ota4du0YGi3XOyiz0VfP1Thdkk7Ifx8NvBF4CFgGXJ5Huxy4Of+9DFgoabKkU4E5wIoW17t0OjFhdmKdzGxwE+sYZzqwJF9RcwSwNCJ+JOkuYKmk9wKPA28HiIjVkpYCa4CDwJUR0dee6puZWT1qJvuIeAA4Z5DhTwMXDjHNYmDxiGtnZmYt4TtoO0wndY10Sp+/mY2ck705+ZqVgJP9GOdEbWb1cLIfR5z4zWwoTvZmZiXgZG81+YjBbOxzsjczKwEnezOzEnCyP8wa/T+uY0Gr6jqW5tlsrHGyH6PGYmIci3U2Gy+c7EeZn1dvZoeDk/1h4kRrZqPJyd7MrASc7M3MSsDJ3sysBJzszcxKwMnezKwEnOzNzErAyd7MrASc7M3MSqBmspc0S9JtktZKWi3pw3n4JyU9IWllfl1SmOZqSeslPSzp4nbOQKfyTVSt49/SbOQm1jHOQeCjEXGvpOOAeyTdmj/7QkR8rjiypDOAhcCZwCnAzyWdHhF9ray4mZnVr2bLPiK2RMS9+e9dwFpgxjCTLABujIh9EbERWA+c14rKmplZcxrqs5c0GzgHuDsP+qCkByRdL2lqHjYD2FSYrIdBdg6SFknqltTd29vbeM3HCHdBmFknqDvZS5oCfB/4i4jYCXwVOA2YC2wBrqmMOsjkcciAiGsjYl5EzOvq6mq03h3HSb15Q/12/k3NWqeuZC9pEinR3xARPwCIiK0R0RcR/cB1VLtqeoBZhclnAptbV2UbL5zMzQ6feq7GEfA1YG1EfL4wfHphtLcCq/Lfy4CFkiZLOhWYA6xoXZXNzKxR9VyNcz5wGfCgpJV52MeBd0qaS+qieRR4H0BErJa0FFhDupLnSl+JY2Y2umom+4i4k8H74X8yzDSLgcUjqJeZmbWQ76BtkPuZ28+/sVnrOdmbmZWAk72ZWQk42bdYvV0Q7qows8PJyd7MrASc7G3M8VGRWeOc7EfIicfMxgInezOzEnCyNzMrASd7GzPcZWbWPCd7M7MScLI3MysBJ3szsxJwsjczKwEnexs3fALXbGhO9mZmJeBkb2ZWAk72ZmYl4GTfJPcPm9lYUjPZS5ol6TZJayWtlvThPPyFkm6VtC7HqYVprpa0XtLDki5u5wyYmVlt9bTsDwIfjYhXAq8FrpR0BnAVsDwi5gDL83vyZwuBM4H5wFckTWhH5c0G46Mus0PVTPYRsSUi7s1/7wLWAjOABcCSPNoS4NL89wLgxojYFxEbgfXAeS2ut5mZNaChPntJs4FzgLuBkyNiC6QdAjAtjzYD2FSYrCcPG1jWIkndkrp7e3ubqLqVnVvwZvWrO9lLmgJ8H/iLiNg53KiDDItDBkRcGxHzImJeV1dXvdUwM7Mm1JXsJU0iJfobIuIHefBWSdPz59OBbXl4DzCrMPlMYHNrqmtmZs2o52ocAV8D1kbE5wsfLQMuz39fDtxcGL5Q0mRJpwJzgBWtq7KZmTWqnpb9+cBlwAWSVubXJcBngYskrQMuyu+JiNXAUmANcAtwZUT0taX2h5H7h8ceLzOzqom1RoiIOxm8Hx7gwiGmWQwsHkG9zMyshXwHrZlZCTjZm5mVgJO9mVkJONkP4JN6ZjYeOdmbmZWAk72ZWQk42VtpuIvOyszJ3sysBEqf7N3aM7MyKH2yt/HPO3QzJ3szs1JwsjfLfARg45mT/RAqG74TQOfzMjKrzcneSss7CSsTJ3szsxJwsjczKwEneys9d+dYGTjZZ97gzWw8c7I3MysBJ3szsxKomewlXS9pm6RVhWGflPSEpJX5dUnhs6slrZf0sKSL21Vxs2a5y87KqJ6W/TeA+YMM/0JEzM2vnwBIOgNYCJyZp/mKpAmtqqyZmTWnZrKPiDuAZ+osbwFwY0Tsi4iNwHrgvBHUz8zMWmAkffYflPRA7uaZmofNADYVxunJww4haZGkbkndvb29I6hGc3wob2Zl0myy/ypwGjAX2AJck4drkHFjsAIi4tqImBcR87q6upqshpmZ1aOpZB8RWyOiLyL6geuodtX0ALMKo84ENo+simZmNlJNJXtJ0wtv3wpUrtRZBiyUNFnSqcAcYMXIqmjWWdwFaGNRPZdefge4C3i5pB5J7wX+RtKDkh4A3gB8BCAiVgNLgTXALcCVEdHXttqbtZGTuo0nE2uNEBHvHGTw14YZfzGweCSVMjOz1vIdtGZ1ckvfxjInezOzEihdsnfrzMzKqHTJ3sysjMZ9sndL3sysBMnezMyc7M3MSqE0yd7dOVYvrys2Ho27ZO8N1czsUOMu2ZuZ2aGc7M3MSsDJ3sysBJzszcxKwMnerEV8cYB1Mid7sxqcxG08cLI3MyuBcZvs3RozM6sa88m+ktSd3G20eN2zsWDMJ3szM6utnn84fr2kbZJWFYa9UNKtktblOLXw2dWS1kt6WNLF7aq42Vjhlr91gnpa9t8A5g8YdhWwPCLmAMvzeySdASwEzszTfEXShJbV1szMmlIz2UfEHcAzAwYvAJbkv5cAlxaG3xgR+yJiI7AeOK81VTUzs2Y122d/ckRsAchxWh4+A9hUGK8nDzuEpEWSuiV19/b2NlkNs7HD3Tk2mlp9glaDDIvBRoyIayNiXkTM6+rqanE1zNrPydvGkmaT/VZJ0wFy3JaH9wCzCuPNBDY3Xz0zM2uFZpP9MuDy/PflwM2F4QslTZZ0KjAHWDGyKpqNLbVa/D4isNEwsdYIkr4D/D5wkqQe4H8CnwWWSnov8DjwdoCIWC1pKbAGOAhcGRF9baq7mZnVqWayj4h3DvHRhUOMvxhYPJJKmZlZa3X8HbQ+5DUzG7mOT/ZmZeGGjbXTmEv23iDMzBo3ZpK9k7yZWfPGTLIfyMnfxjqvw3Y4jdlkb2Zm9XOyNzMrASd7M7MScLI3a5N6++Tdd2+Hg5O9mVkJONmbmZVAxyV7H9KambVexyV7M3u+2Vf92I0gGzEnezOzEnCyN+tQbs1bK3VssveKbmbWOh2b7M3MrHWc7M3MSsDJ3qzDuAvT2mFEyV7So5IelLRSUnce9kJJt0pal+PU1lTVzCoG7hC8g7BaWtGyf0NEzI2Iefn9VcDyiJgDLM/vzawFnNStWe3oxlkALMl/LwEubcN3mJlZA0aa7AP4maR7JC3Kw06OiC0AOU4bbEJJiyR1S+ru7e0dYTXMzGw4I03250fEucCbgCslvb7eCSPi2oiYFxHzurq6RlgNMwN389jQRpTsI2JzjtuAm4DzgK2SpgPkuG2klTQzs5FpOtlLOlbScZW/gT8AVgHLgMvzaJcDN4+0kmY2OLfkrV4jadmfDNwp6X5gBfDjiLgF+CxwkaR1wEX5vZl1AO8cymtisxNGxAbg7EGGPw1cOJJKmZlZa/kOWrNxqNkWvFv+45eTvZlZCTjZm5WAH69gTvZm45iTulU42ZuZdwol4GRvVmJO8uXhZG9mVgJO9mbWNB8ZjB1O9mZ2CF+9M/442ZuZlYCTvZnVrdLCd0t/7HGyN7OW8U6gcznZm9mQRpq8nfw7h5O9mVkJONmb2YjVasH76p7R52RvZi03VDJvNul75zByTvZm1rH8XP7WcbI3s45T75FBu75nPHKyN7Mxa+B1/42eOygTJ3szG3XtvlmrzEm+om3JXtJ8SQ9LWi/pqnZ9j5mVT6ee2O3knUpbkr2kCcDfAW8CzgDeKemMdnyXmdlINdodNHC6saBdLfvzgPURsSEi9gM3Agva9F1mZg1pNJnX2hkM9b6e2Mi4w9WhFkVEQxPUVaj0NmB+RFyR318GvCYiPlgYZxGwKL99OfA08BRwUosiLSzLsTOjl/H4jV629cVjI6KLekREy1/A24F/KLy/DPhSjWm6WxnbUaZjZ8VOqIOjl+1ox3pf7erG6QFmFd7PBDa36bvMzKyGdiX73wBzJJ0q6UhgIbCsTd9lZmY1TGxHoRFxUNIHgX8BJgDXR8TqGpNd2+LYjjIdOyt2Qh0c2xM7oQ5jJdalLSdozcyss/gOWjOzEnCyNzMrASd7M7MSaMsJ2npIegUwAzgI7I6I3+RHKnwYeIh0F+5W4HfAbNJjFx4GPgN8DDiLdGPB94EbIuLZw1z/aRGxrfD+xIh4ulPKHFhWK+rY6nkuW3ntKtOsHqNyglbSh4ArgX7gpaSEv5mU1PtIO6EoxCDdYStgCrARmJTH3Z/H+0BE3F7jewfd0CQdD1wNXAp05e/pJ+1MngTWAy8DtgFvzMOnAMuBP8nFbQdOAH4JPAhcARybP+sD9uZXrTIn52mOKUx7K/DmPM+VYXuAJ4Yo6w7gAuAW4C3Aply3KXm+fkp6blGr6ufyDi1vNXAmcDfwAdK60gvML5Q1VJkHgbXA5yPiG5i1wGgl+weB1wF3kZ6Zs4q0UT1KSkaTSTdibSRtCJOA44EXkLqe/ivwb8BtwDRS4jsO2EFzG+9zwFTgSNKGthdYA5xK2viK4/bnOkwAns11uhc4B3gcOAXYTdp5VXZKL8n1f4a0MxmuzH7Szk25nOeAk4GbSMliOimR9+bhA+elH9gFvDDXb0p+fzywLte3i3SU1Kr6ubxDy9sCnEa1EbCdtI5BtQEzsMwP5OX1BtKO4rH8+rNcL4DTgQ3AfyTtUP4YuAc4QFr3XgPcHxHfk/QW4D/l73mcdLPjZtL6ejjKW5nn7YpOqGNEHIB/71XYQOpFuGNgOaTeggXA+Xn6Xw9Sn015/P3Aq0k75wuGqceWPM0K0g2nfaTHxEwH5ubvuRu4L5c5CXg3ads9hbReHDvIfK2IepN4I7fbtuoFrMnxvhyPIyXY/aSEdDNpI7sX+Dopca/MsT8vjP2klu1jwLmklepGUnL/HdUN6mAe92B+35f/7ss/5JY8fFEue3P+bEuebjfpqCJIG/i9pG6m/lz33+W4Oy/cStm7c5mRv2dXrteu4crMf/86L/R1pJ1UkHZs/cDOXP5the95Xlm5PpUjo0p5lfpVymt1/Vze88vbkD9bDuzJy2QvsG+YMleS1t8dediPcx36SAlpF2k935XL6qO6rvcV3u/JZeylug3syvO4r43lbcvlPZbH2ZHHebYD6thPSrbdpB3vZp6fEyrl7M3xOdK6sCvXY2B99he+q5+0zPcOUo9K46Gf1BjYnut0sDC8J5dZqc9nSAm9UvZzuZzifG0BbiflxD8Yzccl1PKkpLnAfknHRMQuUsI/QGqNv4G0BzuT1Fd/KnA2qQX1FPBaUqsrSHvJvwWOJrV0j87TPpO/6xHSkcP6PP79pAX+AOkIoNI18jHSOYKncj3OzeX8KiJOJC28Obk+6wBJugDYKemLpJbgItKK0ktaGP8jT/coaYFuIS2oIcsk7cFvAV5EtZugn7SQ95Bak/uBd+b5OaQsSR/P0/2kUN4k0nkQcjmtrN9Iyjuuw8trdn5fmsebCUyW9AtS0pk0TJkzgN9STWrzcnk78rSV1t0eUtLYTlrfKw8RfJa0XUA6gttC2sE/Tlpnjsjjtqu8o3J5q0nr5vV5umM6oI59+fc+N5c3ofDZwHKOBn5Iuil0D2mbHlifnbkuvaT1JPJ4A+uxk5SDHsl1qhxFPpGnWUpaJ/bkcXcCV+Uy+knr4dGFMn5C2pm+K8/fRaT8V9sotexn5h9+8iCfnU9KZNeRukbOzsPeRDqcOTPP6EPAl0iHT3NyfBFpj/3zXNae/IPtzQsvSBvTJ/PC2JwX9G7gzjxeZY/7WB4+K5d1KfCtvDC25PH+OS/wtaQdxWpSn/0DVFsTB0h75r2kHdAva5TZQ2rxbcjTrM/lFlsh+/J4G4Yp60Ae9jjpHMJO0oq5Pr9vZf1aWd66Di+v0fmdT0ooy0jr0+5c5r2DlFnZqfcD/zuP923Surg1DwvSOncf8CuqRx+3kNaT+0hJoVK3X5HWmSeoHv63q7yP5/J+lcc9hmoreLTr+Jf5t96Vf+OHhihnE9UW9yMMXZ+nSdv9vXnc1XmZD6zHblJj8978Hb+l2oNQOVJ4PJffQ1oXDuZ67M/THMjDK+Wsyt+1mtT1vL6evDtm76CVNJW0B/wvpL3lZNJhZA/wjojokXQpqX9vAWmv20VqzfeS+lhPIf2oe0n9ZxPy9J8j9dm/lbRAHwJekeMrSQv0fFJXynF5+KWkE6kX5TiVtLO6nXR+YivpGUFrSXvqPQPivlynSstlL2lj+V2Oe0g7pufy/L45Ii6U9M2IeFcxAu8DvhkRb68Mz7/ZwHG/FRGX5WmuI7X8+6m2VIqRQn0q71UY/8X5dz0pz8Pr8u95ch7eVfj86Tx8R/79Kp8/DZyY5/PUvGxmkjaGPuBDpPMvN+VlU4x/QjpqO2OIz7eSGhHb8nedCfxeLnfCgDgxz9eMXIfi8IM5Hsj1PT3Pw4Rc72I5+/P0x+blO4vqeaGB3zmB1H97Yn7/b8CrSDufjwKfjoi/lvQk6TyVSOvyi/Lv+EPSOnhUXmY7SInmRaQW7Lpcj+PycmxneU+SGmm78/xO7IA6biK1oj9J9eKOwepV6dIJ0nmBc/O0pw6oz0HS+vEc1e2il7S+FusxgbQjP4vqkYNI68/kPH3lHNGG/Ht15TJeksc/gmpjbztpx/D/gL8h7dCWRsRnqGHMJvvhSHpPRHy9GEn9+R+NiE8P8vkNpEPmY6iuSCfkeCIp2byY9CO/pPB+BSlhbCMl90lUF9yTpJVhO2ll66faZXQg/11vhGpiPYK00CuxOLwS95FWpGKcRDV5V6aFauIJ0so6jdRX+fJBYnGn94phPt9CSubrSFdYFX/PemNl4z6CtPFVokgbyJQG40FSd8sppCR8MtUT85UE3kis/IZ9pA10LylZ9BeG9+VhW0mJvDfPX7H7tLhT7aO689+Up/kZKUkti4i9AJJeRbpIYRrpJOFNpEbKVNJ6fDNpuf8+Kfk8R+r+O5m0Q7plQHlnAe9tYXnHA58gXWH0CVLCvWOEZY50nj8XEWskvSQvk78iNRz+aUA5PyBt1x8hNX4E/CK/ivU5IY+zIC+nz+XvrVxsUqnHEaT1bh+pFX8TaRvZS9oBnEQ68byD1FPwT6TtaT6pAfOCXLdn87xsz+WRh30zItZQj9HoxjkM3USPNxj3k5LCflLLsB/46xyfJO2lB4vdpIR8b47FaXeTkl6lK2ktaWN+qIkYeSUJ0pHC/rzyHMh/b8uxN8cv5Wl/U4i9eZrewrjPklb0g6SjlYdIK+vuJuK6QtyTy9qT53t3E/G+HNdTPaG1n+rheaOxv/B+HymZbiO1ug40EVcV3j9ISi4Pklp4K0lXmu0pxFfk755SGGdgvKsQ7yM1KO4b7e1prL+AaSOJLazHicPFdr/G7B20kh7Irz351Z9fAcxqME4i9dNOIh26CfhvOR5DtUtlYLyLtOf+ZY5XkBLSuaQN+FhSEgjSEcI+0obfaNxDahn1kxJ3ZfheUkK8P8eVOVZa1/+Yx709f7anMM5RuW47qLYo9zSzLLI1pCOElbm8nhy3kxJ1o/F40m9YObz9LenIqnI5aaNxJ9VW/pG5/JNIrbAjmohH5jgxD5swYJwjSevP5By353jsMOUeRWoNVs5VbQTmFtbtZl59eft4WNJnc4sUST9tNEp6gaRHJN0oaeuAuFHStyT9q6TPSFon6eYcl+VxNhTGrcTvStok6euS1kq6T9JeSdsl7WgyrpP02xx/BqyU9J+B+wrx9UO8vxK4V9LlwD2Sfp7rtVnSQw3ESj22S3pU0n7gQUl9wNqcdx7LcYOk/ZL+QdIqSZ8YJC6QdJukb0uaJenWPL+/kXTOwI1xUKO91x3BXnIr6frUXtLJ22eAy0nJ65kG4wHgv5MS9eU5/iHVQ/EtQ8RvkhL5DwfE38v16hskdpNabY3G40ktx5uonlSuxK8PEb9Mas1+d5DPdlE9ERWkBHuAat9go7GS3B/P5UTh1ddk3FcoexcpAd5H6uZoNH6ElDyfIR0pVC6dK54UbTRWjjT6c10PknZM/aQGQOXy4Gfz55XL8yrjDIzbqV5+vIF00cBjpKOuDU3E60iNl++QdsY/Je00v0za+TUaV1E96X2QtO0czL/nXqonZh8m7VQr460vjDMwbic1Mnbx/KPTyknvRmNQXacrcSSvg7nsnaT1sN64qxCfzvPVm5fnTlL+Kl5iWjmp25c/Gxgrv+M3SevU23IevBC4q56cOWb77CV9jZS43jNIfDHVfvV64hxS/9hXgfcPiNeTEsY5A2NE/FDpJPCvSZeDFuN2Uh9g94C4kdTKf7rBuIXUN7gl13UOaWOaSUos5w8VI+Ljkt5ca5z8ux5D6hvc2mT8HWlnd5B0omkTjfeHTySdO3mS1Mf9K+D4iPitpNObiXneTqHqBaRL43bQ+PmESpyR/z6SajI8imr//SpSi/+VhWHDxb8nbbw3R8TLJT1ciTTnZVSPQivnc1qhcu6ksqx6SUdNXyTdif6ZHKeRluEK0g1mtc6rTCP9Zq+gegdyo/HjpBOXnyddCn0tqSH3edJNa1+pER8lNQ62F8qtnA97JambsZ5YOb/3GOnE/KZCfDHPP/91PNW++Uoj4sgh4rOkPPAL4DsRca2k+yKiZut+zCZ7s/Emdzn8HLiEdD31JaTuuz9l8CukasUppORVbNw8QkouzTQ4djB4gp5CasFu5dAT+6eT7gD++zzOcQNisbzKDmkvKUlWGjd1x4iYJukp0pVyF+Tf749Il75eAPxrjXhkXhyV83i7qF4t1chFAT8knbz9MunquDVUn+d1fF4+T+S630m6C3c5qbF4J+mkbTG+hrSj+A3pzt9r8vL8BnBNRMyjhom1RjCzw+YdpMuJTwE+RUo8rx9hmSeSHzZIOuJcTGpxTiOdoG4kzmbw5PYq0qW2HyWdJ/oW6SbFSuwhJdvXkZJwMf4LKZFV4t053kE6au1pJEp6Gekc1dXAEqrP25pO2pnUik+RzuWQf7dVpMR/LI31GDxB6kqbQ+quOS2XeTxpZ3wk1SvO+nP93w18OyLeLenGAfFs0hHLuaQegveTLis+j7QzrW20+9798suv2i/gPWM5jtZ3k7pkPtFs7ITfrpHfeNh1aLRXYr/88qv2i8YvJ+6o2Al1GO+x1st99mYdQtID+c85OU6ufDQK1WlGMPw5BBu5ym+5l7R+rAZOj4jJw06F++zNOsnJwMWkx228C7iBdMno33LoXdKdGK8AvsfgV8iNlXno9HgF6WT3paS7fd9CulqtJid7s87xI9JJz2VUb61/hPRI70YvJz7sMSJ+IGljrvPAOCbmodNj/o3fFhG/lLQ8Ih6VdDt1cDeOmVkJHDHaFTAzs/ZzsjczKwEnextXJJ0g6QP571Mkfe8wfOftkublv3+i6sPGPpQfonWDpMn5oVorJb2j3XUyG8gnaG28OYH8rJOI2Ay87XB+eURcUnj7AeBNEbFR0muBSREx93DWx6zCyd7Gm88Cp0laSXpQ3Csj4ixJ7yZdrjaB9IySa0i3rF9Gup39koh4RtJpwN+RHuK2G/iziHio+AWSjiZdUngG1f88VvnsUdLt7J8GXgosk/Rt0i3tXblefxwRj7Rh3s2G5GRv481VwFkRMVfSbNLljBVnkZ5aehTpcbF/FRHnSPoC6br2L5KekvjnEbFO0mtIT0O8YMB3vB/YHRGvlvRq0j+veZ6I+HNJ84E3RMRTku4G/jIi/rCVM2tWLyd7K5PbImIXsEvSs6QnE0L6PwGvljQF+A/Ad6V/v+FzsDsTXw/8X4CIeKBw56tZx3KytzLZV/i7v/C+n7QtHAHsqLNf3Teo2Jjiq3FsvKk8K71hEbET2Cjp7QBKzh5k1DtIz5iv/LPuVzdZV7PDxsnexpWIeBr4paRVwP9poog/Bd4r6X7SQ6YWAEj6I0n/K4/zVWBK7r75GOk/MZl1ND8uwcysBNyyNzMrASd7M7MScLI3MysBJ3szsxJwsjczKwEnezOzEnCyNzMrgf8PjdzYCAZcHvgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "from dateutil.parser import parse\n",
    "from datetime import timedelta\n",
    "from lifelines import KaplanMeierFitter\n",
    "from lifelines.plotting import add_at_risk_counts\n",
    "\n",
    "tmm = pd.read_csv(\"../mdata/tmm_start_small.csv\")\n",
    "# tmm['first.seen'] = [parse(i) for i in tmm['first.seen']]\n",
    "# tmm['last.seen'] = [parse(i) for i in tmm['last.seen']]\n",
    "tmm['time.diff'] = [(parse(x) - parse(y) + timedelta(days=1)).days for x, y in zip(tmm['last.seen'], tmm['first.seen'])]\n",
    "tmm['status'] = [ (1 if parse(x) > parse(\"2020-12-31\") else 0)  for x in tmm['last.seen']]\n",
    "\n",
    "tmm[\"time.diff\"].groupby(tmm[\"time.diff\"]).count().plot(kind=\"bar\")\n",
    "tmm['status'].describe()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9a5e8801-250b-4f51-810f-e0e7b0cd33f9",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA70AAAIPCAYAAAC/qLlmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABJMElEQVR4nO3de5xddXno/8/DZGASiIFJAsUEkglFKRcNV61UDYo0KidojQLaFvuyai2g4o/TQ2sPVaivptWjrcdbqVLwFMELVKlFKCLBihcCCIRwSWBCYEAJIZgEyJALz++PvQI7wyQzzOyZtWbN5/16zWv2Wuu7v+uZPNkz+9nf7/quyEwkSZIkSaqjXcoOQJIkSZKkkWLRK0mSJEmqLYteSZIkSVJtWfRKkiRJkmrLoleSJEmSVFsWvZIkSZKk2ppQdgCDMW3atJw9e3bZYUiSJEmSKuqWW25Zk5nT++4fE0Xv7Nmzufnmm8sOQ5IkSZJUURGxqr/9Tm+WJEmSJNWWRa8kSZIkqbYseiVJkiRJtTUmruntz+bNm+np6aG3t7fsUCqjo6ODmTNn0t7eXnYokiRJklQJY7bo7enpYfLkycyePZuIKDuc0mUmjz/+OD09PXR1dZUdjiRJkiRVwpid3tzb28vUqVMteAsRwdSpUx35liRJkqQmY7boBSx4+/DfQ5IkSZK2N6aL3rHgyiuvZNGiRQB897vf5a677nru2LnnnssPf/jDF93nAw88wKGHHtqyGCVJkiSprsbsNb1jwZYtW1iwYAELFiwAGkXviSeeyMEHHwzAeeedV2Z4kiRJklR7Fr3D9PWvf53PfOYzRASveMUraGtro7Ozk1/+8pccccQRHHbYYdx88828+93v5sorr+SGG27gb//2b7n88ss5//zzOfHEE1m4cCFLlizhIx/5CE899RS77bYb1113HY8//jh/9Ed/xFNPPQXAF77wBV7zmteU/BNLkiRJ0tjR0qI3Ii4ETgRWZ+YL5t9G46LTfwLeAjwNvDczb21lDKNp2bJlfOpTn+LGG29k2rRprF27lo997GMsX76cH/7wh7S1tXHRRRcB8JrXvIYFCxY8V+Q227RpEyeffDLf/OY3Ofroo1m/fj0TJ05k77335tprr6Wjo4MVK1Zw6qmncvPNN5fwk0qSJEnS2NTqkd6LgC8AX9/B8TcDBxZfrwK+XHwflk/+xzLuemT9cLvZzsEvfQl/8z8O2WmbH/3oRyxcuJBp06YB0NnZCcA73/lO2traBn2ue++9l3333Zejjz4agJe85CUAPPXUU5xxxhncdttttLW1sXz58qH8KJIkSZI0brW06M3MH0fE7J00OQn4emYm8POI2DMi9s3MX7UyjtGSmf2umLz77ru3pJ/Pfe5z7LPPPtx+++08++yzdHR0DDlWSZIkSRqPRvua3hnAQ03bPcW+YRW9A43IjpQ3vvGNvP3tb+ess85i6tSprF27dqftJ0+ezIYNG16w/6CDDuKRRx5hyZIlHH300WzYsIGJEyeybt06Zs6cyS677MLFF1/M1q1bR+pHkSSpUj75H8uA8v7GS5IKPzin8f3Ni8qNYxhGu+jt70ay2W/DiA8AHwDYf//9RzKmITvkkEP4+Mc/zutf/3ra2to4/PDDd9r+lFNO4f3vfz+f//zn+c53vvPc/l133ZVvfvObnHnmmWzcuJGJEyfywx/+kD//8z/nHe94B9/+9rc57rjjXvQIsiRJY9Vdj6xnfe/mssOQJP16KfT+puwohiUaM41b2GFjevP3d7CQ1T8DizPz0mL7XmDeQNObjzrqqOy7gNPdd9/N7/zO77Qs7rrw30XSSLn+3tVs2GgRotHx+etWsMsuwX+d9fqyQ5Gk8e1f39ooej90Y9mRDCgibsnMo/ruH+2R3iuBMyLiMhoLWK0bq9fzStJ4s2HjZqZPdm0BjY72CbuwecuzZYchSaqBVt+y6FJgHjAtInqAvwHaATLzK8BVNG5XdB+NWxb9SSvPL0mSJElSs1av3nzqAMcTOL2V55QkSZIkaUdGe3qzJEmSVE0rroXedWVHIVXLU4/BLruUHcWwWPRKkiRJ0Ch499in7Cikamlrh61jeyHLsV2yS5IkSZK0Exa9w/DAAw9w6KEvuDMT5557Lj/84Q93+txPfOITfOYznxmp0CRJkiRJOL15RJx33nllhyBJkiRJwqJ32LZu3cr73/9+fvrTnzJjxgy+973v8aEPfYgTTzyRhQsXctVVV/Gxj32MadOmccQRR9Dd3c33v/99AO666y7mzZvHgw8+yEc/+lE+/OEPl/zTSJIk1dBgF6ja3DvysUgadU5vHqYVK1Zw+umns2zZMvbcc08uv/zy54719vbywQ9+kB/84Af85Cc/4bHHHtvuuffccw/XXHMNN910E5/85CfZvHlsXyAuSZJUSdsWqBroa69ZZUcqaQTUY6T3B+fAr5e2ts/fOgzevGjAZl1dXcydOxeAI488kgceeOC5Y/fccw9z5syhq6sLgFNPPZULLrjgueNvfetb2W233dhtt93Ye++9efTRR5k5c2ZLfwxJkiRJGs8c6R2m3Xbb7bnHbW1tbNmy5bntzBzycyVJkiRJw1ePkd5BjMiW4aCDDqK7u5sHHniA2bNn881vfrPskCSNE9ffu5oNG1t7yUTvlq0t7U+SJGk01KPoraiJEyfypS99ifnz5zNt2jSOOeaYskOSNE5s2LiZ6ZM7yg5DkiSpdBa9wzB79mzuvPPO57bPPvvsF7Q57rjjuOeee8hMTj/9dI466iigcZ/eZs39SJIk1cpgV08eKa7KLI1rFr0j7F/+5V+4+OKL2bRpE4cffjgf/OAHyw5JkiRpdG1bPVmSSmDRO8LOOusszjrrrLLDkKRhue7uR7nx/jVlh6FxZNXjT/PSKU7RlyQNn6s3S5IGdOP9a1j1+NNlh6FxZNbUSRw5a6+yw5Ak1cCYHunNTCKi7DAqY6BbJEnScMyaOolzTzyk7DA0jjy2weswJUnDN2ZHejs6Onj88cct9AqZyeOPP05Hh1PBJEmSJGmbMTvSO3PmTHp6enjsscfKDqUyOjo6mDlzZtlhSJIkSVJljNmit729na6urrLDkCRJGr8GeysibxkkqURjtuiVJElSyep0K6LlV0P34rKjkKpn7UqYMrZnk47Za3olSZKklule3HhzL2l7nV2w36vKjmJYHOmVJEmSoPHmfv6isqOQqufJR8uOYFgc6ZUkSZIk1ZZFryRJkiSptix6JUmSJEm1ZdErSZIkSaotF7KSJEkaLYO9r+1Y4f13JY0BFr2SJEmjpU73tZWkMcLpzZIkSZKk2rLolSRJkiTVltObJakCrr93NRs2bm5Zf71btrasL0mSpLHMoleSKmDDxs1Mn9xRdhiSJEm1Y9ErSZKk+lp+NXQvHrjd2pXQ2TXi4UgafV7TK0mSpPrqXtwoaAfS2QVz5o10NJJK4EivJEnScLyYe+96X9tydHbB/EVlRyGpJBa9kiRJw+G9dyWp0pzeLEmSJEmqLYteSZIkSVJtWfRKkiRJkmrLoleSJEmSVFsWvZIkSZKk2rLolSRJkiTVlrcskqQRcv29q9mwcfOg2vZu2TrC0UiSJI1PFr2SNEI2bNzM9MkdZYchjQ8rrm3cL7cMm3vLOa8kaVAseiVJ0tjXuw722KfsKNSf5VdD9+Lyzr92JXR2lXd+SaXzml5JkiSNnO7FjcKzLJ1dMGdeeeeXVDpHeiVJkjSyOrtg/qKyo5A0TjnSK0mSJEmqLYteSZIkSVJtWfRKkiRJkmrLoleSJEmSVFsuZCVJNXPd3Y9y4/1rWtrnqsefZtbUSS3tU5IkaTRY9EpSzdx4/5qWF6mzpk7i2AOmtaw/jW9THr6BCZvWD9iufeNmaNtrcJ1u7h1mVJKkurLolaQamjV1EueeeEjZYUj9mrBpPZsnTh+w3TO5CfbwwxZJ0vB4Ta8kSZIkqbYseiVJkiRJtdXyojci5kfEvRFxX0Sc08/xKRHxHxFxe0Qsi4g/aXUMkiRJkiRBi6/pjYg24IvAm4AeYElEXJmZdzU1Ox24KzP/R0RMB+6NiEsyc1MrY5EkSdIIWn41dC8euN3aldDZNeLhSNKOtHqk9xjgvszsLorYy4CT+rRJYHJEBLAHsBbY0uI4JEmSNJK6FzcK2oF0dsGceSMdjSTtUKtXb54BPNS03QO8qk+bLwBXAo8Ak4GTM/PZFschSZKkkdbZBfMXlR2FJO1Uq0d6o5992Wf794HbgJcCc4EvRMRLXtBRxAci4uaIuPmxxx5rcZiSJEmSpPGg1UVvD7Bf0/ZMGiO6zf4EuCIb7gNWAgf17SgzL8jMozLzqOnTB76XnyRJkiRJfbV6evMS4MCI6AIeBk4B3t2nzYPAG4H/joh9gJcD3S2OQ5IkjaIpD9/AhE3rB9V2l629IxyNJEnPa2nRm5lbIuIM4BqgDbgwM5dFxJ8Vx78CnA9cFBFLaUyH/l+ZuaaVcUiSpNE1YdN6Nk90ZpYkqXpaPdJLZl4FXNVn31eaHj8CnNDq80qSJEmS1Ferr+mVJEmSJKkyWj7SK0mSJI1JT6yC9o6yo5Cqp2NK2REMi0WvJEmSBI2C97CFZUchqcWc3ixJkiRJqi1HeiVJkkbD8quhe3HZUbTO2pXQ2VV2FJI0IEd6JUmSRkP34kahWBedXTBnXtlRSNKAHOmVJA3KQ088RceEtrLDUEW1b9zMM7mppX1O2q2Gb1M6u2D+orKjkKRxpYZ/TSRJI6FjQhsL5s4oOwyNthXXQu+6gdvNmAh7TRv5eCRJepEseiVJ0o71roM99ik7CkmShsxreiVJkiRJtWXRK0mSJEmqLYteSZIkSVJtWfRKkiRJkmrLoleSJEmSVFsWvZIkSZKk2rLolSRJkiTVlkWvJEmSJKm2JpQdgCRJ0pi1/GroXjy4tmtXQmfXiIajHXhiFbR3DNyuY8rIxyJp1Fn0SpIkDVX34sEXs51dMGfeSEek/rR3wGELy45CUkkseiXV2vX3rmbDxs2lnLt3y9ZSzitplHV2wfxFZUchSdoBi15JtbZh42amTx7ElDZJkiTVkgtZSZIkSZJqy6JXkiRJklRbFr2SJEmSpNryml5JKtl1dz/KjfevaVl/qx5/mllTJ7WsP0mSpLHMkV5JKtmN969h1eNPt6y/WVMncewB01rWnyRJ0ljmSK8kVcCsqZM498RDyg5DkiSpdhzplSRJkiTVliO9kiTVxYproXdda/vc3Nva/iRJGmUWvZIk1UXvOthjn7KjGH3Lr4buxeWce+1K6Owq59ySpEFxerMkSRrbuhc3is8ydHbBnHnlnFuSNCiO9EqSpLGvswvmLyo7CklSBVn0SpIkaWQ9sQraO8o7f8eU8s4tqXQWvZIkSRpZ7R1w2MKyo5A0TnlNryRJkiSpthzplTQk19+7mg0bN5cdxoB6t2wtOwRJkiSVyKJX0pBs2LiZ6ZNLvD5LkiRJGgSnN0uSJEmSasuiV5IkSZJUW05vliSp6lZcC73rBm63uXfkY5EkaYyx6JUkqep618Ee+5QdhSRJY5LTmyVJkiRJteVIryRJGj3Lr4buxa3tc+1K6OxqbZ+SpNqw6JWkEXDd3Y9y4/1rBtV21eNPM2vqpBGOSKqI7sWtL1I7u2DOvNb1J0mqFYteSRoBN96/ZtDF7Kypkzj2gGmjEJVUEZ1dMH9R2VFIksYJi15JGiGzpk7i3BMPKTsMSZKkcc2FrCRJkiRJteVIryRJkobmiVXQ3jFwu44pIx+LJO2ARa8kSZKGpr0DDltYdhSStFNOb5YkSZIk1ZYjvZIklWHFtdC7bnBtN/eObCySJNWYRa80Rl1/72o2bNxc2vl7t2wt7dxSLfSugz32KTsKSZJqz6JXGqM2bNzM9MmDWDxEkiRJGscseiWphh564ik6JrS1tM/JE9tb2p8kSdJosOiVpBrqmNDGgrkzyg5D48Xyq6F78eDarl0JnV0jGo4kSc1avnpzRMyPiHsj4r6IOGcHbeZFxG0RsSwibmh1DJIkaRR1L24Us4PR2QVz5o1kNJIkbaelI70R0QZ8EXgT0AMsiYgrM/OupjZ7Al8C5mfmgxGxdytjkCRJJejsgvmLyo5CkqQXaPVI7zHAfZnZnZmbgMuAk/q0eTdwRWY+CJCZq1scgyRJkiRJQOuL3hnAQ03bPcW+Zi8D9oqIxRFxS0T8cYtjkCRJkiQJaP1CVtHPvuznnEcCbwQmAj+LiJ9n5vLtOor4APABgP3337/FYUqSJEmSxoNWj/T2APs1bc8EHumnzdWZ+VRmrgF+DLyyb0eZeUFmHpWZR02fPr3FYUqSJEmSxoNWj/QuAQ6MiC7gYeAUGtfwNvse8IWImADsCrwK+FyL45AkSaqeJ1ZBe0fZUbROx5SyI5CkAbW06M3MLRFxBnAN0AZcmJnLIuLPiuNfycy7I+Jq4A7gWeCrmXlnK+OQJKmlVlwLveta2+fm3tb2p7GhvQMOW1h2FJI0rrR6pJfMvAq4qs++r/TZ/jTw6VafW5KkEdG7DvbYp+woJEnSELT6ml5JkiRJkirDoleSJEmSVFstn94sSZJqYvnV0L144HZrV0Jn14iHI0nSUDjSK0mS+te9uFHQDqSzC+bMG+loJEkaEkd6JUnSjnV2wfxFZUchSdKQOdIrSZIkSaoti15JkiRJUm1Z9EqSJEmSasuiV5IkSZJUWxa9kiRJkqTacvVmSdL4teJa6F03cLvNvSMfi8auJ1ZBe8fg2nZMGdlYJEkvYNErSRq/etfBHvuUHYXGuvYOOGxh2VFIknbA6c2SJEmSpNqy6JUkSZIk1ZZFryRJkiSptix6JUmSJEm1ZdErSZIkSaoti15JkiRJUm1Z9EqSJEmSasuiV5IkSZJUWxa9kiRJkqTasuiVJEmSJNWWRa8kSZIkqbYmlB2ANB5cf+9qNmzc3NI+e7dsbWl/Um2suBZ61w2u7ebekY1Fo+eJVdDeUc65O6aUc15J0qBY9EqjYMPGzUyfXNKbMWm86V0He+xTdhQabe0dcNjCsqOQJFWQ05slSZIkSbVl0StJkiRJqi2nN0uSVAfLr4buxa3tc+1K6OxqbZ+SJI0yR3olSaqD7sWNIrWVOrtgzrzW9ilJ0ihzpFeSpLro7IL5i8qOQpKkSnGkV5IkSZJUW470SlIFPPTEU3RMaGtZf5MntresL0mSpLHMoleSKqBjQhsL5s4oOwxJkqTacXqzJEmSJKm2LHolSZIkSbXl9GZJtXXd3Y9y4/1rSjn3qsefZtbUSaWcW5IkSc+z6JVUWzfev6a04nPW1Ekce8C0UT+vNCY8sQraO1rbZ8eU1vYnSaoNi15JtTZr6iTOPfGQssOQ1Ky9Aw5bWHYUkqRxwmt6JUmSJEm1ZdErSZIkSaotpzdLklRly6+G7sUDt1u7Ejq7RjwcSZLGGkd6JUmqsu7FjYJ2IJ1dMGfeSEcjSdKY40ivRtX1965mw8bNZYcx6nq3bC07BEljWWcXzF9UdhSSJI1JFr0aVRs2bmb65BbfpkKSJEmSdsDpzZIkSZKk2rLolSRJkiTVlkWvJEmSJKm2LHolSZIkSbXlQlaSJGn4nlgF7YNcqLBjysjGIklSE4teSZI0fO0dcNjCsqOQJOkFLHolaYQ89MRTdExoG1TbyRPbRzgaSZKk8cmiV5JGSMeENhbMnVF2GJIkSeOaRa8kSaNt+dXQvXhwbdeuhM6uEQ1HkqQ6c/VmSZJGW/fiRjE7GJ1dMGfeSEYjSVKtOdIr6UW77u5HufH+NWWHMaBVjz/NrKmTyg5D6l9nF8xfVHYUkiTVXstHeiNifkTcGxH3RcQ5O2l3dERsjQiXepTGmBvvX8Oqx58uO4wBzZo6iWMPmFZ2GJIkSSpRS0d6I6IN+CLwJqAHWBIRV2bmXf20+3vgmlaeX9LomTV1EueeeEjZYUiSJEk71eqR3mOA+zKzOzM3AZcBJ/XT7kzgcmB1i88vSZIkSdJzWl30zgAeatruKfY9JyJmAG8HvtLic0uSJEmStJ1WF73Rz77ss/2PwP/KzK077SjiAxFxc0Tc/Nhjj7UqPkmSJEnSONLq1Zt7gP2atmcCj/RpcxRwWUQATAPeEhFbMvO7zY0y8wLgAoCjjjqqb+EsSZJGwxOroL1j4HYdU0Y+FkmShqDVRe8S4MCI6AIeBk4B3t3cIDO7tj2OiIuA7/cteCVJUkW0d8Bh3mhBkjR2tbTozcwtEXEGjVWZ24ALM3NZRPxZcdzreCVJkiRJo6bVI71k5lXAVX329VvsZuZ7W31+SVIFrLgWeteVc+7NveWcV5IkVVLLi15JkuhdB3vsU3YUkiRJFr2SJO3U8quhe3Fr+1y7Ejq7Bm4nSZKGrdW3LJIkqV66FzeK1Fbq7II581rbpyRJ6pcjvZIkDaSzC+YvKjsKSZI0BI70SpIkSZJqy6JXkiRJklRbFr2SJEmSpNqy6JUkSZIk1ZZFryRJkiSptly9WRqDrrv7UW68f01p51/1+NPMmjqptPNLkiRJg+VIrzQG3Xj/GlY9/nRp5581dRLHHjCttPNLkiRJg+VIr4bt+ntXs2Hj5kG17d2ydYSjGT9mTZ3EuSceUnYYkiRJUqVZ9GrYNmzczPTJHWWHIQ3LlIdvYMKm9S3ts33jZmjbq6V9jhmbe8uOQJIkCbDolSQAJmxaz+aJ01va5zO5CfZwGnhlLb8auhcP3G7tSujsGvFwJEnSyPCaXknS+NS9uFHQDqSzC+bMG+loJEnSCHGkV5I0fnV2wfxFZUchSZJGkCO9kiRJkqTasuiVJEmSJNWWRa8kSZIkqba8pleSpLp4YhW0t/gWch1TWtufJEmjzKJXkqS6aO+AwxaWHYUkSZXi9GZJkiRJUm1Z9EqSJEmSasuiV5IkSZJUWxa9kiRJkqTaciErSSNqysM3MGHT+rLDGNAuW3vLDkGSJEkjwKJX0oiasGk9mydOLzsMjRfLr4buxYNru3YldHaNaDiSJKl8Tm+WJNVH9+JGMTsYnV0wZ95IRiNJkirAkV5JUr10dsH8RWVHIUmSKsKRXkmSJElSbTnSK0lS1T2xCto7Bm7XMWXkY5EkaYyx6JUkqeraO+CwhWVHIUnSmOT0ZkmSJElSbVn0SpIkSZJqy6JXkiRJklRbFr2SJEmSpNqy6JUkSZIk1ZarN0uqtdUbetm1rZzP9ybt5q9YSZKksvmOTFKt7dq2C8f+9rSyw5AkSVJJLHolDcmUh29gwqb1A7bbZWvvKEQjSZIk9c+iV9KQTNi0ns0Tp5cdhiRJkrRTFr2SJJXhiVXQ3jG4th1TRjYWSZJqzKJXkqQytHfAYQvLjkKSpNrzlkWSJEmSpNqy6JUkSZIk1ZZFryRJkiSptix6JUmSJEm15UJW0gi77u5HufH+NS3tc9XjTzNr6qSW9ilJkiTVkUWvNMJuvH9Ny4vUWVMncewB01rW3zZTHr6BCZvWD6rtLlt7W35+SZIkqdUseqVRMGvqJM498ZCywxjQhE3r2TxxetlhSJIkSS3jNb2SJEmSpNqy6JUkSZIk1ZbTmyUNyeoNvezaVv3PzSbt5q85tcATq6C9o7V9dkxpbX+SJKlfvhuUNCS7tu3Csb/d+sW0pEpq74DDFpYdhSRJGoLqD9NIkiRJkjREFr2SJEmSpNpqedEbEfMj4t6IuC8izunn+Hsi4o7i66cR8cpWxyBJkiRJErS46I2INuCLwJuBg4FTI+LgPs1WAq/PzFcA5wMXtDIGSZIkSZK2afVI7zHAfZnZnZmbgMuAk5obZOZPM/OJYvPnwMwWxyBJkiRJEtD6oncG8FDTdk+xb0feB/ygxTFIkiRJkgS0/pZF0c++7LdhxHE0it7f28HxDwAfANh///1bFZ8kSZIkaRxp9UhvD7Bf0/ZM4JG+jSLiFcBXgZMy8/H+OsrMCzLzqMw8avr06S0OU5IkSZI0HrR6pHcJcGBEdAEPA6cA725uEBH7A1cAf5SZy1t8fkmSBu+JVdDeMXC7jikjH4skSRoRLS16M3NLRJwBXAO0ARdm5rKI+LPi+FeAc4GpwJciAmBLZh7VyjgkSSVafjV0Ly7n3GtXQmfX4Nu3d8BhC0cuHkmSVLpWj/SSmVcBV/XZ95Wmx38K/GmrzytJqojuxS+++GyVzi6YM2/0zytJkiqr5UWvJEl0dsH8RWVHIUmSZNGr0XPd3Y9y4/1ryg5j1K16/GlmTZ1UdhiSJEnSuGTRq1Fz4/1rxmUBOGvqJI49YFrL+33oiafomNDW0j7bN27mmdw0qLaTdvPXhyRJkqrPd60aVbOmTuLcEw8pO4xa6JjQxoK5M1rbadtesEfrC3RJkiSpLK2+T68kSZIkSZVh0StJkiRJqi2nN0uS6uWJVY377w5Gx5SRjUWSJJXOoleSVC/tHXDYwrKjkCRJFeH0ZkmSJElSbVn0SpIkSZJqy+nNkrT8auheXHYU9bF2JXR2lR2FJEkS4EivJDUK3rUry46iPjq7YM68sqOQJEkCHOmVpIbOLpi/qOwoJEmS1GKO9EqSJEmSasuiV5IkSZJUWxa9kiRJkqTasuiVJEmSJNWWRa8kSZIkqbYseiVJkiRJtWXRK0mSJEmqLYteSZIkSVJtWfRKkiRJkmprQtkBSKqx5VdD9+KyoxjY2pXQ2VV2FJIkSRoBjvRKGjndixsFZdV1dsGceWVHIUmSpBHgSK+kkdXZBfMXlR2FJEmSxilHeiVJkiRJtWXRK0mSJEmqLYteSZIkSVJteU2vNAoeeuIpOia0tbTPyRPbW9qfJEmSVEcWvdIo6JjQxoK5M8oOQ5IkSRp3nN4sSZIkSaoti15JkiRJUm1Z9EqSJEmSastreiVJrffEKmjvKOfcHVPKOa8kSaoki15JL97yq6F78cDt1q6Ezq4RD0cV1N4Bhy0sOwpJkiSnN0sagu7FjYJ2IJ1dMGfeSEcjSZIk7ZAjvZKGprML5i8qOwpJkiRppxzplSRJkiTVlkWvJEmSJKm2nN6sYfnGLx7kwp900z5h4M9PVj3+NL/1kg4e29A7CpFVy+SJ7WWHIEmSJI1LFr0alu/d9jAP/2Yjs6ftPmDbWVMn8YoZU1gwd8YoRCa9SGXeYqeOvG2QJEmqCIteDduMPSdy7omHDKrteBzl1RjhLXYkSZJqyWt6JUmSJEm1ZdErSZIkSaoti15JkiRJUm1Z9EqSJEmSasuFrCQ13Pyv8PMvQ9sgbq+0diV0do18TJIkSdIwWfRKalj6HVj3IHQeMHDbzi6YM29w/Y6VWwF5ix1JkqRasuiV9Lwp+8P8Ra3t01sBSZIkqURe0ytJkiRJqi2LXkmSJElSbVn0SpIkSZJqy6JXkiRJklRbLmQlaWgGuyqzqyJLkiSpRBa9kobGVZklSZI0BrR8enNEzI+IeyPivog4p5/jERGfL47fERFHtDoGSZIkSZKgxSO9EdEGfBF4E9ADLImIKzPzrqZmbwYOLL5eBXy5+K4xapddgsc29A6q7eSJ7SMcjSRJkiQ9r9XTm48B7svMboCIuAw4CWguek8Cvp6ZCfw8IvaMiH0z81ctjkWjpG2XYMHcGWWHoVbYZRd48tHBtfVaXUmSJI0BrS56ZwAPNW338MJR3P7azADGbNH78y+9n8m/ubvsMEpx9qat7DZhF/jXPcsORcP166Ww535epytJkqRaaXXRG/3syyG0ISI+AHwAYP/99x9+ZCPs2S1byg6hFB27wJ7tW6H3N2WHouHacz+YM6/sKCRJkqSWanXR2wPs17Q9E3hkCG3IzAuACwCOOuqoFxTFVfLqP/+XskOQJEmSJPWj1as3LwEOjIiuiNgVOAW4sk+bK4E/LlZxfjWwzut5JUmSJEkjoaUjvZm5JSLOAK4B2oALM3NZRPxZcfwrwFXAW4D7gKeBP2llDJIkSZIkbdPq6c1k5lU0CtvmfV9pepzA6a0+ryRJkiRJfbV6erMkSZIkSZVh0StJkiRJqi2LXkmSJElSbVn0SpIkSZJqy6JXkiRJklRbFr2SJEmSpNqy6JUkSZIk1ZZFryRJkiSptix6JUmSJEm1ZdErSZIkSaoti15JkiRJUm1FZpYdw4Ai4jFgVdlxDGAasKbsILRT5qj6zFG1mZ/qM0fVZ46qzfxUnzmqvjJzNCszp/fdOSaK3rEgIm7OzKPKjkM7Zo6qzxxVm/mpPnNUfeao2sxP9Zmj6qtijpzeLEmSJEmqLYteSZIkSVJtWfS2zgVlB6ABmaPqM0fVZn6qzxxVnzmqNvNTfeao+iqXI6/plSRJkiTVliO9kiRJkqTasuhtgYiYHxH3RsR9EXFO2fEIIuLCiFgdEXc27euMiGsjYkXxfa8yYxzPImK/iLg+Iu6OiGUR8ZFivzmqiIjoiIibIuL2IkefLPabowqJiLaI+GVEfL/YNj8VEhEPRMTSiLgtIm4u9pmjComIPSPiOxFxT/E36XfNUTVExMuL1862r/UR8VHzUy0RcVbxPuHOiLi0eP9QuRxZ9A5TRLQBXwTeDBwMnBoRB5cblYCLgPl99p0DXJeZBwLXFdsqxxbg/8vM3wFeDZxevG7MUXU8A7whM18JzAXmR8SrMUdV8xHg7qZt81M9x2Xm3Kbbd5ijavkn4OrMPAh4JY3XkzmqgMy8t3jtzAWOBJ4G/h3zUxkRMQP4MHBUZh4KtAGnUMEcWfQO3zHAfZnZnZmbgMuAk0qOadzLzB8Da/vsPgm4uHh8MfC20YxJz8vMX2XmrcXjDTTeZMzAHFVGNjxZbLYXX4k5qoyImAm8Ffhq027zU33mqCIi4iXA64CvAWTmpsz8Deaoit4I3J+ZqzA/VTMBmBgRE4BJwCNUMEcWvcM3A3ioabun2Kfq2SczfwWNogvYu+R4BETEbOBw4BeYo0opps7eBqwGrs1Mc1Qt/wj8BfBs0z7zUy0J/FdE3BIRHyj2maPqmAM8BvxrcZnAVyNid8xRFZ0CXFo8Nj8VkZkPA58BHgR+BazLzP+igjmy6B2+6GefS2JLgxARewCXAx/NzPVlx6PtZebWYlrZTOCYiDi05JBUiIgTgdWZeUvZsWinjs3MI2hcAnV6RLyu7IC0nQnAEcCXM/Nw4CkqMA1T24uIXYEFwLfLjkXbK67VPQnoAl4K7B4Rf1huVP2z6B2+HmC/pu2ZNIb1VT2PRsS+AMX31SXHM65FRDuNgveSzLyi2G2OKqiY7reYxnXy5qgajgUWRMQDNC6reUNE/Bvmp1Iy85Hi+2oa1yIegzmqkh6gp5jFAvAdGkWwOaqWNwO3Zuajxbb5qY7jgZWZ+VhmbgauAF5DBXNk0Tt8S4ADI6Kr+CTqFODKkmNS/64ETisenwZ8r8RYxrWICBrXUN2dmZ9tOmSOKiIipkfEnsXjiTT+sN2DOaqEzPzLzJyZmbNp/N35UWb+IeanMiJi94iYvO0xcAJwJ+aoMjLz18BDEfHyYtcbgbswR1VzKs9PbQbzUyUPAq+OiEnFe7s30linpXI5ikxn4g5XRLyFxrVVbcCFmfmpciNSRFwKzAOmAY8CfwN8F/gWsD+NF+k7M7PvYlcaBRHxe8B/A0t5/nrEv6JxXa85qoCIeAWNxSfaaHxA+q3MPC8ipmKOKiUi5gFnZ+aJ5qc6ImIOjdFdaEyj/UZmfsocVUtEzKWxGNyuQDfwJxS/8zBHpYuISTTWzpmTmeuKfb6GKqS4peHJNO7M8UvgT4E9qFiOLHolSZIkSbXl9GZJkiRJUm1Z9EqSJEmSasuiV5IkSZJUWxa9kiRJkqTasuiVJEmSJNWWRa8kSaMgIvaMiD8vHr80Ir7Ton4/ERFnF4/Pi4jjW9GvJEl14S2LJEkaBRExG/h+Zh7a4n4/ATyZmZ9pZb+SJNWFI72SJI2ORcABEXFbRHw7Iu4EiIj3RsR3I+I/ImJlRJwRER+LiF9GxM8jorNod0BEXB0Rt0TEf0fEQX1PEBEXRcTC4vEDEfHJiLg1IpZuax8Ru0fEhRGxpDjHSaP4byBJ0qiz6JUkaXScA9yfmXOB/9nn2KHAu4FjgE8BT2fm4cDPgD8u2lwAnJmZRwJnA18axDnXZOYRwJeL5wB8HPhRZh4NHAd8OiJ2H/JPJUlSxU0oOwBJksT1mbkB2BAR64D/KPYvBV4REXsArwG+HRHbnrPbIPq9ovh+C/AHxeMTgAXbrgMGOoD9gbuH9yNIklRNFr2SJJXvmabHzzZtP0vjb/UuwG+KUeKh9LuV5//mB/COzLx3aKFKkjS2OL1ZkqTRsQGYPJQnZuZ6YGVEvBMgGl45xDiuAc6MYsg4Ig4fYj+SJI0JFr2SJI2CzHwcuLFYwOrTQ+jiPcD7IuJ2YBkw1AWozgfagTuKWM4fYj+SJI0J3rJIkiRJklRbjvRKkiRJkmrLoleSJEmSVFsWvZIkSZKk2rLolSRJkiTVlkWvJEmSJKm2LHolSZIkSbVl0StJkiRJqi2LXkmSJElSbVn0SpIkSZJqy6JXkiRJklRbFr2SJEmSpNqy6JUkSZIk1ZZFryRJkiSptix6JUmSJEm1NW6K3ohYEBHnFI/fFhEHNx07LyKOH0KfsyPizlbGKUmSJElqncjMsmMYcRExITO3NG1fBHw/M78zzH5nF/0cOrwIJUmSJEkjoTZFb0T8MXA2kMAdwFZgLXA4cCuwFDgK+AbwfWBd8fUO4H9TFMERcTTwT8DuwDPAG4GpwP8r9gGckZk/teiVJEmSpGqbUHYArRARhwAfB47NzDUR0Ql8FngZcHxmbo2I9wIUxeqVNI30RsS2fnYFvgmcnJlLIuIlwEZgNfCmzOyNiAOBS2kU0JIkSZKkCqtF0Qu8AfhOZq4ByMy1RSH77czc+iL6eTnwq8xcUvSzHiAidge+EBFzaYwgv6yFsUuSJEmSRkhdit6gMa25r6da1M9ZwKPAK2ks/tX7IvuVJEmSJJWgLqs3Xwe8KyKmAhTTm3dmAzC5n/33AC8truslIiZHxARgCo0R4GeBPwLaWha5JEmSJGnE1GKkNzOXRcSngBsiYivwywGechnwLxHxYWBhUz+bIuJk4P9GxEQa1/MeD3wJuDwi3glcz4sfQZYkSZIklaA2qzdLkiRJktRXXaY3S5IkSZL0Aha9kiRJkqTasuiVJEmSJNVWaUVvREyIiDUR8Xd99v/Vi+jjqxFx8E6OL46Io4YY3+si4taI2BIRC5v2HxcRtzV99UbE24ZyDg1PRFwYEasj4s5+jp0dERkR08qITQ395SgiOiPi2ohYUXzfq8wY1b+dvb5UDeao+iJifkTcGxH3RcQ5ZcejFzJH1WZ+qm+s5KjMkd4TgHtp3GoomvYPquiNiLbM/NPMvGtEooMHgfcC32jemZnXZ+bczJwLvAF4GvivEYpBO3cRML/vzojYD3gTjRyqXBfxwhydA1yXmQfSuN1YZX9BjnMX0c/rS5VyEeaosiKiDfgi8GbgYODUnX1Qr9FnjqrN/FTfWMpRmUXvqcA/0ShMXg0QEYuAicUI6iV9nxART0bEeRHxC+B3t43kRkRbRFwUEXdGxNKIOKvP83aJiIsj4m8HG1xmPpCZdwDP7qTZQuAHmfn0YPtV62Tmj4G1/Rz6HPAXgEuTl2wHOToJuLh4fDHwttGMSYOzk9eXKsIcVd4xwH2Z2Z2Zm2jcLvGkkmPS9sxRtZmf6hszOSql6C3ugftG4PvApTQKYDLzHGBjMZL6nn6eujtwZ2a+KjN/0rR/LjAjMw/NzMOAf206NgG4BFiemX/d4h/llCJ+VURELAAezszby45FO7RPZv4KoPi+d8nxSNJImAE81LTdU+xTdZijajM/1TdmclTWSO+JwPXFCOnlwNuL4fGBbC3a99UNzImI/xsR84H1Tcf+mUah/KnhBt0sIvYFDgOuaWW/GrqImAR8HDi37FgkSeNe9LPPGUjVYo6qzfxU35jJUVlF76nA8RHxAHALMBU4bhDP683MrX13ZuYTwCuBxcDpwFebDv8UOC4iOoYZc1/vAv49Mze3uF8N3QFAF3B78X9rJnBrRPxWqVGpr0eLD422fXi0uuR4JGkk9AD7NW3PBB4pKRb1zxxVm/mpvjGTo1EveiPiJcDvAftn5uzMnE2jUD21aLI5ItpfZJ/TgF0y83LgfwNHNB3+GnAV8O2ImDDc+JucilObKyUzl2bm3k3/r3qAIzLz1yWHpu1dCZxWPD4N+F6JsUjSSFkCHBgRXRGxK41Loq4sOSZtzxxVm/mpvjGTozJGev8A+FFmPtO073vAgojYDbgAuKO/hax2YgawOCJuo7Ga5V82H8zMzwK3Av8vIgb1M0fE0RHRA7wT+OeIWNZ0bDaNTzVueBExqsUi4lLgZ8DLI6InIt5Xdkza3g5ytAh4U0SsoLHK9qIyY1T/fH1VnzmqtszcApxB4zKou4FvZeaynT9Lo8kcVZv5qb6xlKPIrOS0a0mSJEmShq3MWxZJkiRJkjSiLHolSZIkSbVl0StJkiRJqq0hF70R0R4RiyJiRUTcGRE3RcSbWxlcK0XERRGxsEV9nVb83Csi4rSBn6HRFhHzI+LeiLgvIs4pOx69kDmqPnNUbean+sxR9ZmjajM/1TdWcjTkhawiYhGwL/CBzHwmIvYBXp+Z32plgAPE0NbffXt30PYi4PuZ+Z1hnrMTuBk4isbNl28BjizuFawKiIg2YDmNlYF7aCynfmpm3lVqYHqOOao+c1Rt5qf6zFH1maNqMz/VN5ZyNKSR3oiYBLwfOHPbrYcy89FtBW9EnBARP4uIWyPi2xGxR7H/gYj4ZLF/aUQcVOx/fUTcVnz9MiImR8Oni1HkpRFxctF2XkRcHxHfAJZGRFvRbklE3BERHyzaRUR8ISLuioj/BPYe7j9W4feBazNzbVHoXgvMb1Hfao1jgPsyszszNwGXASeVHJO2Z46qzxxVm/mpPnNUfeao2sxP9Y2ZHA11evNvAw9m5vq+ByJiGvDXwPGZeQSNUdGPNTVZU+z/MnB2se9s4PTMnAu8FthI436+c4FXAscDn46IfYv2xwAfz8yDgfcB6zLzaOBo4P0R0QW8HXg5cBiNAv01Q/xZ+5oBPNS03VPsU3WYo+ozR9VnjqrN/FSfOao+c1Rt5qf6xkyOJoxAn68GDgZujAiAXYGfNR2/ovh+C43CFuBG4LMRcQlwRWb2RMTvAZcW05cfjYgbaBS164GbMnNl8dwTgFc0Xa87BTgQeF3T8x+JiB+16OeLfvZ5s+NqMUfVZ46qzxxVm/mpPnNUfeao2sxP9Y2ZHA216L0P2D8iJmfmhj7Hgsb031N38Nxniu9bt50/MxcVU5DfAvw8Io6n/3/EbZ7qc74zM/Oa7YKIeAsj84/eA8xr2p4JLB6B82joeoD9mrZnAo+UFIv6Z46qzxxVm/mpPnNUfeao2sxP9Y2ZHA1penNmPg18Dfh8ROwKEBH7RsQfAj8Hjo2I3y72T4qIl+2sv4g4IDOXZubf05gOfRDwY+Dk4prd6TRGbm/q5+nXAB+KiPair5dFxO7F808pnr8vcNxQftYdnO+EiNgrIvaiMdJ8zQDP0ehaAhwYEV3F/89TgCtLjknbM0fVZ46qzfxUnzmqPnNUbean+sZMjoYzvfmvgb8F7oqIXhqjr+dm5mMR8V7g0ojYrant8p309dGIOI7G6O9dwA+ATcDvArfTGLH9i8z89bbFr5p8FZgN3BqN+dSPAW8D/h14A7C0OPcNw/hZn5OZayPifBpJBjgvM9e2om+1RmZuiYgzaHwY0QZcmJnLSg5LTcxR9ZmjajM/1WeOqs8cVZv5qb6xlKMh37JIkiRJkqSqG+rqzZIkSZIkVZ5FryRJkiSptka16I2IT0TEwxFxW9PXni0+x1+1sr8dnOO0iFhRfJ020ufTixMRHRFxU0TcHhHLIuKTZcek7UXEy/v8HlgfER8tO67xKCIujIjVEXFnn/1nRsS9xWvoH8qKb7wzP9XXX44i4ptNv98eiIjbSgxx3DNH1WZ+qq8OORrVa3oj4hPAk5n5mRE8x5OZuccI9t9JY4Xpo2gssHULcGRmPjFS59SLUyxotntmPlms6v0T4COZ+fOSQ1M/IqINeBh4VWauKjue8SYiXgc8CXw9Mw8t9h0HfBx4a2Y+ExF7Z+bqMuMcr8xP9fWXoz7H/w+wLjPPG/XgBJijqjM/1VeHHFVienNE/CIiDmnaXhwRR0bE7sUnC0si4pcRcVJx/L0RcUVEXF2Mtv5DsX8RMLH4xOGS4vn/WYz43RkRJ7cg3N+ncR/itUWhey0wvwX9qkWy4clis734csW26nojcL8Fbzky88dA3xXoPwQsysxnijYWVCUxP9W3gxwBz30I+y7g0lENStsxR9VmfqqvDjkqo+g9q2ko/Ppi32U0/rGIxj11X5qZt9D4JPtHmXk0jfvsfjoa9+AFmAucDBxG436++2XmOcDGzJybme+hUYw+kpmvLD6VuLoF8c8AHmra7in2qUKicX/m24DVND6k+EXJIWnHTqHivyjHoZcBry0+kLwhIo4uOyBtx/yMHa8FHs3MFWUHoh0yR9VmfqpvTOSojKL3c0VROjczjyv2fQt4Z/H4XcC3i8cnAOcUxctioAPYvzh2XWauy8xeGvf2ndXPuZYCx0fE30fEazNzXQvij372OYpYMZm5NTPnAjOBYyLiBVMxVL5o3Mh8Ac+/5lUNE4C9gFcD/xP4VvFJrqrB/Iwdp+KHelVnjqrN/FTfmMhRJaY3Z+bDwOMR8Qoao7eXFYcCeEdTkbx/Zt5dHHumqYutNN4E9O13OXAkjeL37yLi3BaE2wPs17Q9E3ikBf1qBGTmb2h8YOIU9Gp6M3BrZj5adiDaTg9wRXGpwE3As8C0kmPS88zPGBARE4A/AL5ZdizqnzmqNvNTfWMpR5UoeguXAX8BTMnMpcW+a4Azt32CHRGHD6KfzcXiRUTES4GnM/PfgM8AR7QgzmuAEyJir4jYi8Zo9DUt6FctEhHTo1gVPCImAscD95QalHZkTHw6OA59F3gDQES8DNgVWFNmQNrOdzE/Y8HxwD2Z2VN2INohc1Rt5qf6xkyOyr6m97aImF3s/w6Na/u+1dT2fBqLEN1RLJF9/iD6v6BofwmN631vKqZHfxz42+EGn5lriziWFF/nFftUHfsC10fEHTRydG1mfr/kmNRHREwC3gRcUXYs41lEXAr8DHh5RPRExPuAC4E5xe/dy4DTcjSX+tdzzE/17SBH4HoFlWGOqs38VF8dcjSqtyySJEmSJGk0VWl6syRJkiRJLWXRK0mSJEmqLYteSZIkSVJtWfRKkiRJkmrLoleSJEmSVFsWvZIkSZKk2rLolSRJkiTVlkWvJEmSJKm2LHolSZIkSbVl0StJkiRJqi2LXkmSJElSbVn0SpIkSZJqy6JXkiRJklRbFr2SJEmSpNoal0VvRMyOiDv72X9eRBw/wHM/ERFnj1x0kiRJkqRWmVB2AFWSmeeWHYMkSZIkqXXG5UhvoS0i/iUilkXEf0XExIi4KCIWAkTEWyLinoj4SUR8PiK+3/TcgyNicUR0R8SHS4pfkiRJkjSA8Vz0Hgh8MTMPAX4DvGPbgYjoAP4ZeHNm/h4wvc9zDwJ+HzgG+JuIaB+ViCVJkiRJL8p4LnpXZuZtxeNbgNlNxw4CujNzZbF9aZ/n/mdmPpOZa4DVwD4jGagkSZIkaWjGc9H7TNPjrWx/fXMM47mSJEmSpIoYz0XvztwDzImI2cX2ySXGIkmSJEkaIkco+5GZGyPiz4GrI2INcFPZMUmSJEmSXrzIzLJjqKSI2CMzn4yIAL4IrMjMz5UdlyRJkiRp8JzevGPvj4jbgGXAFBqrOUuSJEmSxhBHeiVJkiRJteVIryRJkiSptkoreiNiQkSsiYi/67P/r15EH1+NiIN3cnxxRBw1xPg+FhF3RcQdEXFdRMwq9h8XEbc1ffVGxNuGcg4NT0TsFxHXR8TdEbEsIj5S7P9mU34eKKapqwQ7ylHT8bMjIiNiWlkxjmc7eQ2dX/zuuy0i/isiXlp2rONVRFwYEasj4s6mfeZnjIiIPSPiOxFxT/E6+92yY9LzBvobpXKZn7EhIuZHxL0RcV9EnFN2PDtS2vTmiHgL8HHgt4DfziKQiHgyM/cYxPPbMnPrAG0WA2dn5s1DiO844BeZ+XREfAiYl5kn92nTCdwHzMzMp1/sOTQ8EbEvsG9m3hoRk4FbgLdl5l1Nbf4PsC4zzysrzvFsZzmKiP2ArwIHAUdm5poyYx2PdpQfoCcz1xdtPgwcnJl/Vl6k41dEvA54Evh6Zh5a7HuJ+RkbIuJi4L8z86sRsSswKTN/U3JYKgzmfYTKY36qLyLagOXAm4AeYAlwahVzVOb05lOBfwIeBF4NEBGLgInFp9eX9H1CRDwZEedFxC+A3902khsRbRFxUUTcGRFLI+KsPs/bJSIujoi/HWxwmXl9UyH7c2BmP80WAj+w4C1HZv4qM28tHm8A7gZmbDterLz9LuDSciLUADn6HPAXgAsLlGRH+dlWUBV2xxyVJjN/DKzts8/8jAER8RLgdcDXADJzkwVvtQz0PkLlMj9jwjHAfZnZnZmbgMuAk0qOqV+l3Kc3IiYCbwQ+COxJowD+WWaeExFnZObcHTx1d+DOzDy36Gfb/rk03qht+xR8z6bnTAAuKZ73qSGG/D7gB/3sPwX47BD7VAtFxGzgcOAXTbtfCzyamStKCUrbac5RRCwAHs7M25texypR39dQRHwK+GNgHXBceZGpP+ZnTJgDPAb8a0S8ksYo1Ucy86lyw1J/dvA+QhVhfiprBvBQ03YP8KqSYtmpskZ6TwS2jaReDry9GB4fyNaifV/dwJyI+L8RMR9o/hT8nxlGwRsRfwgcBXy6z/59gcOAa4bSr1onIvag8f/io31GQE7FUd5KaM4RsIXGpQ3nlhmTntffaygzP56Z+9H40PCMMuPTC5mfMWECcATw5cw8HHgKqOz1buPZTt5HqALMT6X1N3JRydlHZRW9pwLHR8QDND75nMrgPqnu7e863sx8AnglsBg4ncZ1gtv8FDguIjpebJARcTyNN+cLMvOZPoffBfx7Zm5+sf2qdSKincYvwksy84qm/ROAPwC+WVZsaugnRwcAXcDtxe+AmcCtEfFb5UU5fu3oNdTkG8A7RjcqvQjmp7p6aFwfv21k6js0imBVyCB+B6pE5qfyeoD9mrZnAo+UFMtOjXrRW1zj8nvA/pk5OzNn0yhUTy2abC7+g7+YPqcBu2Tm5cD/Zvs/Kl8DrgK+XRRCg+3zcBqjxAsyc3U/TRxFLFlxze7XgLszs+808+OBezKzZ/Qj0zb95Sgzl2bm3k2v/x7giMz8dYmhjks7eg1FxIFNzRYA94x2bNox8zM2FL/THoqIlxe73ghUbnGX8WyA9xEqmfkZE5YAB0ZEV7FY3ynAlSXH1K8yrun9A+BHfUZOvwf8Q0TsBlwA3BERt2bmewbZ5wwa18xsK+L/svlgZn42IqYA/y8i3pOZzw6iz08De9AolgEezMwF8Nx1BfsBNwwyPo2MY4E/ApbG87cl+qvMvIrGi84PJcq3sxypfP3mB3hf8Ub9WWAV4MrAJYmIS4F5wLSI6AH+BniL+RkzzgQuKd4MdgN/UnI82p5/o6rN/FRcZm6JiDNoXO7ZBlyYmctKDqtfpd2ySJIkSZKkkVbmLYskSZIkSRpRFr2SJEmSpNqy6JUkSZIk1daQi96IaI+IRRGxIiLujIibIuLNrQyulSLioohY2KK+Tit+7hURcVor+lRrRcT8iLg3Iu6LCO+LWEHmqPrMUbVFxIURsToi7iw7Fr2Q+ak+c1R95qjaxlJ+hjPSez6wL3BoZh4K/A9gckuiGqSIaBvN8xXn7KSxeuargGOAv4mIvUY7Du1Y8f/ii8CbgYOBUyPi4HKjUjNzVH3maEy4CJhfdhDaoYswP1V3Eeao6i7CHFXZRYyR/Ayp6I2IScD7gTO33XooMx/NzG8Vx0+IiJ9FxK0R8e2I2KPY/0BEfLLYvzQiDir2vz4ibiu+fhkRk6Ph08Uo8tKIOLloOy8iro+Ib9BYwrytaLckIu6IiA8W7SIivhARd0XEfwJ7D/cfq/D7wLWZuTYznwCuZYwkexw5BrgvM7szcxNwGXBSyTFpe+ao+sxRxWXmj4G1Zceh/pmf6jNH1WeOqm0s5WeoI72/TeO+tev7HoiIacBfA8dn5hHAzcDHmpqsKfZ/GTi72Hc2cHpmzgVeC2ykcT/fucArgeOBT0fEvkX7Y4CPZ+bBwPuAdZl5NHA08P6I6ALeDrwcOIxGgf6aIf6sfc0AHmra7in2qTrMUfWZo+ozR5IkqRYmjECfr6YxFe7GiADYFfhZ0/Eriu+30ChsAW4EPhsRlwBXZGZPRPwecGlmbgUejYgbaBS164GbMnNl8dwTgFc0Xa87BTgQeF3T8x+JiB+16OeLfvZ5s+NqMUfVZ46qzxxJkqRaGGrRex+wf0RMzswNfY4Fjem/p+7guc8U37duO39mLiqmIL8F+HlEHE//b7i2earP+c7MzGu2CyLiLYzMG7QeYF7T9kxg8QicR0PXA+zXtD0TeKSkWNQ/c1R95kiSJNXCkKY3Z+bTwNeAz0fErgARsW9E/CHwc+DYiPjtYv+kiHjZzvqLiAMyc2lm/j2N6dAHAT8GTi6u2Z1OY+T2pn6efg3woYhoL/p6WUTsXjz/lOL5+wLHDeVn3cH5ToiIvYoFrE4o9qk6lgAHRkRX8f/zFODKkmPS9sxR9ZkjSZJUC8NZvfmvgceAu4plqr8LPJaZjwHvBS6NiDtoFMEHDdDXR4sFq26ncT3vD4B/B+4Abgd+BPxFZv66n+d+FbgLuLWI459pjCD/O7ACWErj+uEbhv6jPi8z19JYuXpJ8XVesU8VkZlbgDNofBhxN/CtzFxWblRqZo6qzxxVX0RcSuPyoZdHRE9EvK/smPQ881N95qj6zFG1jaX8RKaXaEmSJEmS6mk4I72SJEmSJFWaRa8kSZIkqbZGteiNiE9ExMMRcVvT154tPsdftbK/HZzjtIhYUXydNtLn04sXEfMj4t6IuC8izik7Hm0vIjoi4qaIuD0ilkXEJ8uOSS9ULAT4y4j4ftmxjFcRcWFErC7WrNi27/yIuKP4G/pfEfHSMmMc7yJiv4i4PiLuLn6ffaTY3xkR1xbvFa4tFr/UKNtRfpqOnx0RGRHTyopxvDNH1bezHEXEmcV77mUR8Q9lxrkzo3pNb0R8AngyMz8zgud4MjP3GMH+O2msMH0UjVsi3QIcmZlPjNQ59eJERBuwHHgTjduuLAFOzcy7Sg1Mz4nGTbx3z8wni5XXfwJ8JDN/XnJoahIRH6Pxu+4lmXli2fGMRxHxOuBJ4OuZeWix7yWZub54/GHg4Mz8sxLDHNeKO0Tsm5m3RsRkGu8L3kZjUc+1xW0ZzwH2ysz/VV6k49OO8pOZd0XEfjQWRD2Ixnu5NWXGOl6Zo+rbye+5fYCPA2/NzGciYu/MXF1iqDtUienNEfGLiDikaXtxRBwZEbsXn3IvKUYbTiqOvzciroiIq4tPUP+h2L8ImFh8+n1J8fz/LEaT7oyIk1sQ7u/TuA/x2qLQvRaY34J+1TrHAPdlZndmbgIuA04qOSY1yYYni8324stV9SokImYCb6XxZkMlycwfA2v77FvftLk7vnZKlZm/ysxbi8cbaKx2PoPG352Li2YX03iDqFG2k/wAfA74C3wNlcocVd9OcvQhYFFmPlMcq2TBC+UUvWc1TW2+vth3GfAueO6ThJdm5i00Pjn4UWYeTeM+u5+Oxj14AeYCJwOH0bif736ZeQ6wMTPnZuZ7aBSjj2TmK4tPyK9uQfwzgIeatnt4/oWpajBHY0AxdfY2YDWND5J+UXJI2t4/0nij8WzJcagfEfGpiHgIeA9wbtnxqCEiZgOHA78A9snMX0HjDSOwd4mhie3zExELgIcz8/Zyo1Izc1R9fX7PvQx4bTGAeUNEHF1qcDtRRtH7uaIonZuZxxX7vgW8s3j8LuDbxeMTgHOKN8aLgQ5g/+LYdZm5LjN7adynd1Y/51oKHB8Rfx8Rr83MdS2IP/rZ56dP1WKOxoDM3JqZc4GZwDERcWjJIakQEScCq4sPH1VBmfnxzNwPuITG/ZRVsojYA7gc+Gif0XhVQHN+gC00Blb8wKhCzFH19fN7bgKwF/Bq4H8C3youYaucSkxvzsyHgccj4hU0Rm8vKw4F8I6mInn/zLy7OPZMUxdbafyj9+13OXAkjeL37yKiFS+cHmC/pu2ZwCMt6FetY47GkMz8DY0PtbxMoDqOBRZExAM0fh+/ISL+rdyQtAPfAN5RdhDjXbE2weXAJZl5RbH70WL22rZZbJWd9ld3/eTnAKALuL34PTcTuDUifqu8KMc3c1R9O/g91wNcUVy2dhON2WGVXHCsEkVv4TIaU+mmZObSYt81wJnbPjGIiMMH0c/mIilEY0XLpzPz34DPAEe0IM5rgBMiYq9iJcYTin2qjiXAgRHRFRG7AqcAV5Yck5pExPQoVm6PiInA8cA9pQal52TmX2bmzMycTeP186PM/MOSw1IhIg5s2lyAr51SFe9RvgbcnZmfbTp0JbDtDg+nAd8b7djUf34yc2lm7p2Zs4vfcz3AEZn56xJDHbfMUfXt5Pfcd4E3FG1eBuwKVHKxsReMjo6CsyKi+c3T2zLzAeA7wD8B5zcdO5/GdWV3FP/YDwADrSB6QdH+VuDrNK4DfhbYTONi62HJzLURcT6NwgrgvMxcu7PnaHRl5paIOIPGhxFtwIWZuazksLS9fYGLi5W2dwG+lZneFkfqIyIuBeYB0yKiB/gb4C0R8XIan6ivAly5uVzHAn8ELC0uxwL4K2ARjal+7wMe5PnLuDS6+s1PZl5VXkjqwxxV345+z10IXBiN2+ptAk7L0bw10IswqrcskiRJkiRpNFVperMkSZIkSS1l0StJkiRJqi2LXkmSJElSbVn0SpIkSZJqy6JXkiRJklRbFr2SJEmSpNqy6JUkSZIk1ZZFryRJkiSptix6JUmSJEm1ZdErSZIkSaoti15JkiRJUm1Z9EqSJEmSasuiV5IkSZJUWxa9kiRJkqTasuiVJEmSJNWWRa8kSZIkqbYseiVJkiRJtWXRK0mSJEmqLYteSZIkSVJtWfRKkiRJkmrr/wfIEJfzA8LMrwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 972x540 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "tmm_filter = lambda x:  (x['first.seen'] <= \"2020-10-01\") & (x['group'] == 1) & (x['idTeam'] == 1)\n",
    "tmm2_0 = pd.read_csv(\"../mdata/tmm_start.csv\")\n",
    "# tmm['first.seen'] = [parse(i) for i in tmm['first.seen']]\n",
    "# tmm['last.seen'] = [parse(i) for i in tmm['last.seen']]\n",
    "tmm2_0['time.diff'] = [(parse(x) - parse(y) + timedelta(days=1)).days for x, y in zip(tmm2_0['last.seen'], tmm2_0['first.seen'])]\n",
    "tmm2_0['status'] = [ (0 if parse(x) > parse(\"2020-12-31\") else 1)  for x in tmm2_0['last.seen']]\n",
    "tmm2 = tmm2_0[(tmm2_0.idTeam == 1) & (tmm2_0.group == 1)]\n",
    "tmm2 =tmm2 [['time.diff','status','severity']]\n",
    "\n",
    "\n",
    "fig, axes = plt.subplots(1,1, figsize=(13.5, 7.5))\n",
    "\n",
    "t_survs = []\n",
    "ax = axes\n",
    "for severity in tmm2['severity'].unique():\n",
    "    mask = tmm2['severity'] == severity\n",
    "    ds = tmm2[mask]\n",
    "    t_surv = KaplanMeierFitter()\n",
    "    t_surv.fit(ds['time.diff'], ds['status'], label=severity)\n",
    "    t_surv.plot_cumulative_density(ax=ax)\n",
    "    t_survs.append(t_surv)\n",
    "\n",
    "add_at_risk_counts(*t_survs, ax=ax)\n",
    "plt.tight_layout()\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "732b5598-0689-409e-9c77-055f2a7f2912",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bac74753-0632-4ab0-a55c-70884826a5ed",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}