{ "cells": [ { "cell_type": "markdown", "id": "42e833db-bf17-48c6-a4ea-7c5f1d310360", "metadata": {}, "source": [ "# Non-parametric lifetime models" ] }, { "cell_type": "markdown", "id": "e828b3dd-0bec-4b88-8731-0b857ba98176", "metadata": {}, "source": [ "A built-in dataset of circuit breaker lifetimes." ] }, { "cell_type": "code", "execution_count": 1, "id": "264693a3-5de5-4457-9632-48ee20a0e452", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[34. 28. 12. ... 42. 42. 37.]\n", "[ True True True ... False False False]\n", "[33. 27. 11. ... 31. 31. 26.]\n" ] } ], "source": [ "from assetlife.datasets import load_circuit_breaker\n", "\n", "dataset = load_circuit_breaker()\n", "print(dataset[\"time\"])\n", "print(dataset[\"event\"])\n", "print(dataset[\"entry\"])" ] }, { "cell_type": "code", "execution_count": 2, "id": "5cd6ca17-3980-4e90-bf2c-1b8fa35843a1", "metadata": {}, "outputs": [], "source": [ "from assetlife.lifetime_models import KaplanMeier, NelsonAalen" ] }, { "cell_type": "code", "execution_count": 3, "id": "713b3e2f-7265-43a0-b3f9-a8e28dc26f35", "metadata": {}, "outputs": [], "source": [ "kaplan_meier = KaplanMeier(\n", " dataset[\"time\"],\n", " event=dataset[\"event\"],\n", " entry=dataset[\"entry\"],\n", ")" ] }, { "cell_type": "code", "execution_count": 4, "id": "98e52f10-2452-4c4f-8bb4-1c0968a89422", "metadata": {}, "outputs": [], "source": [ "nelson_aalen = NelsonAalen(\n", " dataset[\"time\"],\n", " event=dataset[\"event\"],\n", " entry=dataset[\"entry\"],\n", ")" ] }, { "cell_type": "code", "execution_count": 5, "id": "87dfc944-1c3e-4dea-b489-0c3672c519bf", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0MAAAF0CAYAAADhIOLTAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvctoD+AAAAAlwSFlzAAAPYQAAD2EBqD+naQAAY6pJREFUeJzt3Xd8k1X7P/BP2nRC96SllNlC2VBAloxHRpFalgyliiIIPKIyBAEZglrgQUVBRGUPQYbsimxlFcoQSkF2yy6l0N2mI+f3R3+9v02btAmkSdp83q9XX5qTc5Irg9y5cs65bpkQQoCIiIiIiMjMWBg7ACIiIiIiImNgMkRERERERGaJyRAREREREZklJkNERERERGSWmAwREREREZFZYjJERERERERmickQERERERGZJSZDRERERERklpgMEf1/MTExcHZ2xs6dO40dSrkxh8dIRFRRjBkzBg0bNjR2GEbh7u6Ozz//3NhhEDEZItPSrl07tG7dukR7TEwM6tati6CgIFy/fr1c7js/Px8pKSnIyckpl9svy507d+Ds7AxnZ2dMnz5dbZ/vv/9e6vPPP//ofB/GfoxERBXNlClT4OzsjMDAQGRlZZW4vmPHjujYseNz3XZmZiZSUlJeNESDyM/PR2BgIJydnbF58+YXvr3k5GS1zyeRoTEZIpOSmpqK1NRUlbZjx47h5ZdfhoWFBfbs2YN69eoZKbrypVQqkZKSAoVCgeXLlyM/P79En19++QUKhQIpKSnIy8vT+T4aN26MZ8+eISwsTB8hExFVellZWUhJScG1a9ewcOHCEtenpaUhLS3N8IEZWGRkJK5duwZra2v89NNPxg6HSG+YDJFJ2717N7p37446derg2LFjqFWrlrFDKndhYWF4+PAh9u3bp9J+5swZXLp0CX379n3u27a0tISzszOsrKxeNEwiIrPSpUsXzJ07F0+ePDF2KEbxyy+/oHHjxvjqq69w6NAh3Lp1y9ghEekFkyEyWatWrULfvn3Rrl07HD58GJ6entJ1jx8/lpaLOTs7w8vLCy+99BIWLVoEpVIp9Tt27BicnZ1x8OBBrFy5Ek2aNEG1atXw6quv4sKFC2XG8Dz389tvv6FFixbw9vZGjx49cOnSJZ0ed/369dGmTRusWrVKpX3lypUICAhAu3bt1I7LycnBvHnzEBwcDE9PTwQGBmLChAkqM22a9gxpM7boY1y3bh1atmwJNzc3xMTE6PT4iIgqorlz5yIzMxOzZ8/Wqn9sbCyGDh2K2rVrw9vbGy+//DJ27NhR5rgbN25g2LBhqF+/Pnx9fdG5c2esXr1a5ZgDAIcOHUJISAh8fHxQvXp1hIWF4dSpUyp93nzzTbRu3RqpqakYMWIE/Pz8ULduXUydOlWn1QUPHjxAZGQkRo8ejTfeeANOTk5Yvny52r4NGzaUjpkeHh5o0qQJpkyZgvT0dK3ua/369ejcuTOqVasGf39/vPXWW7h9+7Z0fUZGBpydnbFgwQJERUWhS5cu8Pb2Rps2bRAZGan1YyIqxGSITNL//vc/vPPOO+jbty8iIyPh4OCgcr2Hhwfi4uKkv6ioKLzzzjuYPHky5s6dK/XLy8tDSkoKVqxYgQsXLmDLli3Ys2cPFAoFOnbsiMuXL5cah6738+uvv+L06dPYsGED/vjjDyQmJqJXr14679F55513sGPHDjx79gwAoFAosGHDBgwbNkxt/9zcXPTo0QPffvstJkyYgLNnz2LZsmXYv38/unbtCoVCAUD9niFtxxY+xlWrVuH06dNYv349li9fDgsLfowQUeUXEBCAESNGYOnSpbhx40apff/++2+0atUKGRkZ2Lx5M06dOoU+ffrg9ddfx7JlyzSOy87ORteuXfHw4UNs2LABZ8+exRdffIG//vpLZbXAb7/9hm7dusHHxwcHDhxAZGQkbGxs0LFjR5V+GRkZSE5Oxvvvv4/evXvjxIkTmDx5MubPn48FCxZo/dhXrlwJe3t7hIeHw97eHm+//TZWrlypNqE6efKkdMz8559/MHv2bKxfvx5vv/12mffz0Ucf4b333kNoaCiOHTuGXbt2ITk5GW3atMHdu3cBAEIIpKSk4Ny5c/jmm2/wzTff4NixY/D390dYWJhK4kSkFUFkQho2bCgsLCwEANGrVy+Rn5+v0/iJEycKLy8v6fLhw4cFAPHKK6+o9EtNTRUuLi6iT58+Utv58+cFALF58+bnvp9XX31VpV9h+9atW8u8zdu3bwsAYubMmSI5OVnY2tqKJUuWCCGE2Lhxo7CwsBD37t0TixYtEgBEdHS0NHbx4sUCgDh06JDKbV69elVYWlqKn3/+WeNj1HZs4WPp2rVrmY+FiKiy+OijjwQA8ezZM5GQkCAcHBzEgAEDpOubNm0qmjZtKl1WKpWiXr16onnz5iIvL0/ltsaOHStcXV1FVlaWEEKIt99+W/j6+krXnz17VgAQkZGRGuPJzc0V3t7eonXr1irt+fn5IiAgQAQEBEhtYWFhAoDYt2+fSt+wsDDh5+en1eNXKpWiVq1aYsyYMVLb1atXhUwmE9u3b9fqNjZu3CgAiLt370ptlpaWYvLkydLlqKgoAUDMnz9fZWxmZqbw8fERI0eOFEIIkZaWJgAIPz8/oVAopH5JSUnC2tpaTJkyRauYiArxJ10yOV5eXggKCsKBAwewa9cujf3+/PNPvPbaa6hZsyZcXV3h7OyMxYsXIyEhoUQRhv79+6tcdnBwQEhICPbv3w8hRKnx6HI/ffr0UbncvHlzAFD5FdHb21tl6V3RGaZCTk5O6NOnj7RUbuXKlejevTt8fX3Vxrht2zZ4e3ujS5cuKu0BAQGoXbs2Dhw4oPHx6Tq2+HNJRGQuPD098cknn2DLli0llqQVio2NxfXr1zFo0CBYWlqqXNezZ088ffoU58+fVzu2Ro0asLOzw+eff479+/dLM/NFXbx4EY8ePcLgwYNV2i0sLDBo0CBcu3YNcXFxUnuVKlXQrVs3lb7NmzfHvXv3pNtfuXKlynHJ2dlZqnJ38OBB3L59G2PGjJHGBwQEoGvXrmpnuS5duoTw8HAEBATAzc0Nzs7OGD58OACUWg1227ZtAIA33nhDpd3Ozg6dOnUqcSzq1asXrK2tpcuurq6oUaNGmbN2RMUxGSKT4+joiCNHjqB+/foYMGCA2hKeO3fuREhICPz8/LB9+3ZcvXoVcXFxmDhxIgCUWJZWrVq1ErdRrVo1ZGRkIDMzU2Msut6Pj4+PymUnJycAkJa7AQXlRFNSUqS/7Oxstfc9bNgwnD59Gvv378f+/fs1LpEDCtZzP3nyBO7u7nB3d4ebmxtcXV3h4uKCW7duITExUW9jq1evrvG2iIgquwkTJqBatWrScaC4Bw8eAAC++OIL6TO18HO1MIHR9Jns7u6OHTt2IC8vD927d4eTkxO6du2KlStXSnuGCscWP94UbXv8+HGJtqKcnJwghEBycjIASFVKi/4V/lD4yy+/AADat2+vkiwdO3YMf/zxB+7fvy/d7vXr1/HSSy8hMTERq1atQmxsLOLi4rBmzRrpfjQpfN6aNm2qcixydXXF77//XuI50/S4ih5vibQhN3YAROp4eHjg8OHD6NatG4YMGYLc3FyVX4t++eUX1KxZEz/88IPKuKIHgKLUHXgSExNhY2MDe3t7jXHoej+a9s8UnX1KSEhQuWxra6t2TLdu3eDr64uhQ4dKM0WaODs7o06dOoiKilJ7vVyu+Z+6rmOL/hJHRGRu7O3tMXv2bIwYMQLbt28vcb2zszOAgmRI0z6ZKlWqaLz9bt26oVu3bnj8+DGOHz+ODRs24N1338X9+/fx2WefSbev7rhWeGxycXGR2krb11l4LHr33XdLzDQ5OTnhyZMn2L59O3799VeEhISUGB8cHIwVK1ZI58b79ddfkZGRgfXr18PNzU3ql5SUpDGGQs7OzrCwsMCFCxfUPj8ymUzlsjbHWyJtcGaITJarqysOHTqEVq1aITw8XKW6Wm5uLhwdHVX6Z2VlaazU8+eff6pczsnJwcGDB9GxY8cSH7BF6Xo/2nByclL5dU1TMmRhYYHw8HA8fvwYgwcPho2NjcbbDAkJwfXr15GQkFBiqYOzszOqVq1aLmOJiMzRO++8g4YNG+LTTz8tUUSgefPm8PLywr59+9R+pmp7egNPT0/07dsXmzZtQr169bB//34ABTMnTk5O2LNnT4kxe/bsgY+PD+rWravT47G2ti4Ro0wmw5o1ayCEQO/evdU+jpCQEKxYsUJKQHJzc2FhYVHiuLFhw4YyYwgJCYFSqcTff/+t9r4KV1oQ6RuTITJpTk5O2LdvH9q1a4d3330XP//8M4CCddcXLlzA+vXroVQqcf/+fQwZMkTao1NcWloaFi1ahKysLDx9+hQjR47Eo0ePMGvWrFLvX9f70bfZs2fj2bNnak/0V9THH3+MoKAg9OnTB4cOHUJOTg7y8vLw77//Ytq0afjtt9/KZSwRkTmytLTE3LlzcfXqVcTGxqpcZ2VlhR9++AF79+7F2LFjpWVkycnJ2Lt3b6knvd6/fz8++ugjXLx4ETk5ORBC4MCBA7h79y7atGkDoGA1wcyZM7F371589dVXSE9PR2pqKiZNmoTTp08jIiKi1B/5dLFs2TK0a9euREXXQiEhIYiLi5MStR49ekCpVGL69OnIyclBcnIyPv30U40/+hW/rf79++ODDz7A+vXrkZGRASEE7ty5g0WLFuGLL77Qy2MiKo7JEJk8BwcH7N27F126dMGoUaOwePFijB07Fp988gk++OAD2NraomPHjhg0aBD+85//qL2Njz76CAkJCfD19YWbmxtOnTqFXbt2oX379qXet673o29WVlZwdnYuc2mag4MDjh8/jl69euGtt96CnZ0dnJyc0L9/fzg4OKBHjx7lMpaIyFz17t0bnTt3Vntd//79ceTIEVy9ehV169aFjY0N6tati8WLF+ODDz7QeJvt27dHvXr18Pbbb8PZ2Rn29vYYOXIkpk6dqpIMjBs3DsuWLcOaNWvg7OwMFxcX7N27F5s2bcJbb72ll8d37NgxXLlyBT179tTYp0uXLrC1tZX2FXXs2BE///wzfv31V1SpUgX16tWDvb09PvnkE63u87fffsO0adPwxRdfwMHBAVWqVEGnTp1w69YthIeH6+VxERUnE1xcSSYkLS0NANT+CpWfn4+0tDTIZDKV6XKFQiEtIVMoFMjKyoKTkxNkMhmOHDmCLl264I8//pA+0LOzs9X+SlV4+1WqVFG7hKG0+8nLy0N6ejqqVq1aYo9NcnIybG1ty/xlTKlUIjU1tcy+OTk5yMzMhIODQ4lKRUVjlcvlJa4v6zGWNra0x0hEVFllZWVBoVBIn/dFFR4LLC0tNc6eCCGQnZ0NOzu7EtdlZmYiLy+vxHJsoODzWghR5udtYSEfdT+aZWRkQKlUloit+DFMncI+pR0vgILjthCixGMoesxUd/xISUmBjY2NxuNdXl4ehBBq7zs5ORl2dnYllo+np6dDJpOVuieLqDgmQ1SpqUuGiIiIiIgALpMjIiIiIiIzxWSIiIiIiIjMEpfJUaXGfS5EREREpAmTISIiIiIiMktcJkdERERERGapQqwbUiqVePDgARwcHPR2IjEiIiqbEAJpaWnw8fGBhQV/PyvE4xIRkfHo89hUIZKhBw8ewM/Pz9hhEBGZrbt376J69erGDsNk8LhERGR8+jg2VYhkqPBkYXfv3lV7YjIiIiofqamp8PPz03hCSXPF4xIRkfHo89hUIZKhwiUIjo6OPOgQERkBl4Kp4nGJiMj49HFs4gJwIiIiIiIyS0yGiIiIiIjILDEZIiIiIiIis1Qh9gwRUeWVn5+P3NxcY4dhtqysrGBpaWnsMCotvr+Nz9rammXhiUgjJkNEZBRCCDx69AjJycnGDsXsOTs7w9vbm0US9Ijvb9NhYWGBWrVqwdra2tihEJEJYjJEREZR+EXR09MT9vb2/CJuBEIIZGZm4vHjxwCAatWqGTmiyoPvb9NQeHLchw8fokaNGnwdiKgEJkNEZHD5+fnSF0U3Nzdjh2PW7OzsAACPHz+Gp6cnl8zpAd/fpsXDwwMPHjxAXl4erKysjB0OEZkYnRfRJiQk4KuvvsKAAQNw9OhRrcbcvXsXU6dOxZtvvolZs2bhyZMnOgdKRJVH4R4Ke3t7I0dCwP+9Dtzboh98f5uWwuVx+fn5Ro6EiEyRTsnQunXrEBwcjNTUVGzduhXx8fFljomLi0OLFi1w5coVdO7cGUePHkVwcDCSkpKeO2giqhy4ZMU0VPTXITExEXPnzsXAgQPx3nvvYevWrWWOmTRpEnr27Kny9+WXX+o1ror+vFYWfB2IqDQ6LZPr3Lkzbt68CWtra8ybN0+rMXPmzEH16tWxdetWWFhYIDw8HHXr1sW3336LL7744rmCJiIiAoDbt2+ja9euGDJkCAYMGID79+9j5MiROHz4MBYvXqxx3OnTp+Hj44O33npLavP29jZEyEREZEJ0SoaqV6+u8x1ERkZizJgxUllLW1tbhIaGIjIyUudkKDMnD/KcPJ1jMFV2Vpb8xYqoAsnPz8fDhw/h4eEBGxubEtdnZGQgNzcXzs7OOt1uSkoKLC0tUbVqVZ3GWFtbS3t+irp//z6EECptzs7OOt1+ReHl5YXY2FiVJWnu7u54++238cUXX5T6WtSuXRs9e/Y0QJQVR2JiIqysrNQ+b7m5uUhOToabm5tOpaoVCgXS0tLg5uam9TFPoVAgMzMTLi4uJa579uwZMjIyVNpsbGzg4eGhdUxERIXKtYBCVlYWHj16hBo1aqi016hRAxs2bNA4TqFQQKFQSJdTU1MBAK2/PAgLm8qzBjuomiM2j2oLXfIhJlBU2Z28abgltG3r6La5/e7du6hVqxYOHz6Mzp07AyioyLZ06VJ89913iI+Pl76U9+vXD9OnTy/x+VfUjRs3EB4ejvPnz0OpVKJTp05Ys2ZNqVXd1q5di6lTpyI1NRUKhQKvvfYali1bBkdHR6mPv78/XFxcVBK2zz//HMOHD9fp8VYE6vbleHp6QgiBjIyMUpOhyMhIxMTEoFq1aggJCUFYWFg5RlrAkO9vQPf3eGhoKJo1a4alS5dKbadPn8Znn32Go0ePomrVqsjKykLLli0xbdo0dO/eXeNt5eXl4cMPP8SKFStgaWkJd3d3LF26FCEhIRrHXL58GaNHj8bp06chl8tRo0YNrFixAm3atJH6jB07Ftu2bVNJlFq0aIGdO3fq9FiJiIByToYKE5riB6uqVasiOztb47iIiAh8/vnn5RmaSbj8MBUNZ/6p0xh/N3vMCm0oXbaRW5SZHOl6MCQi7Y0ZMwabNm3CL7/8gtdeew1yuRwpKSn4/fffsW/fPrz33ntqx+Xn5yMsLAz16tXD06dPkZubi1dffRVDhgzBkSNH1I7566+/MGzYMKxevRpDhw5FcnIyBgwYgPfeew+bNm1S6bt27VqznPUQQuCbb75B06ZN4evrq7Gfk5MT2rVrhxYtWiA2Nhbh4eEYNmwYvv/+e7X9Nf1IV9nt378fvXv3xkcffYSNGzfC1dUVOTk5OHPmDH799ddSk6E5c+Zg69atOH/+PAICAjB//nz069cPV65cQc2aNUv0VygU6N27N4KDg/H06VPY2trim2++QUhICK5duwZ3d3ep7+uvv45Vq1aVwyMmIn0rjx+BMtLS9HZb5ZoMVa1aFZaWlnj27JlKe1JSUqm/1k2ZMgXjx4+XLqempsLPzw+np/1H5dfPikoI4PWlJ3H5oe4H0/ikTLyzKlq6XDw5UufEjSd6mU1iUkWk6vjx41i6dCnWr1+Pfv36Se1OTk545513Sh17+PBhXL58Gdu2bZN+MJozZw66du2K2NhYNGxY8t/1zp07Ub9+fQwdOhRAwdK3yZMno0ePHnjw4AF8fHykvnl5eXjy5InKF0hzMHXqVJw4cQLHjh0rtd/69eulZYMDBgxAkyZN0K9fP4wcORKNGjUq0d9cfqQrKj8/HyNHjkSPHj0wf/58qd3a2hrt2rVDu3btNI5VKpX48ccf8cEHH6BBgwYAgMmTJ2PRokVYvnw55syZU2LMxYsXcfv2bezYsUNa/jlhwgR89dVXWLduHT7++GOV/klJSXB0dGS5bCJ6IeWaDMnlcjRo0AAXL15Uab948SKaNGmicZyNjY3a9fj21nLYW1eOUyPt+bADjlxN1GnMrF2xiE/KVGkrnhypo03CpM0MU1mZPZMlMjdbtmyBs7MzhgwZovPY06dPw93dHQEBAVJbhw4dpOvUJUN2dnbIyMiAEEL695qWlgYhBE6fPo0+ffpIffv27Qt7e3sIITBs2DB89dVXlXLPUFFz5szBokWLEBkZiaZNm5bat/hz0atXL8hkMpw/f15tMqTpR7rK7MyZM4iLi8OSJUt0HhsXF4fExES0b99earOwsEC7du1w+vRptWMKE6D09HSpLTc3FwqFAlFRUSp9165di23btiErKwvt2rXDokWL0LhxY53jJCLSe2axadMm7NmzB6tXrwYADB06FF9//TU+/fRTVK9eHbGxsfjjjz/w008/6fuuKxSZTAZbK91ObhjRtzEUeUrpsrrkSB19JUzFFU+gnncalEkUVVS3b99G7dq1n2vmNTExscSsjZWVFZycnJCYqP6HksGDB2PBggUYN24c3n//fTx8+BAzZsyATCbD48ePpX6ffvopJkyYABcXF5w5cwZ9+/ZFSkqK9LlcGX355ZeYO3cuIiMj8fLLL+s8/tmzZxBCSOekKU7Tj3SV2e3btwEAdevW1Xls4Xu4+Hvc3d0dt27dUjumQYMGaN26NcaOHYvvv/8ejo6OmD9/PnJzc1Xe3+3bt8enn36KRo0a4enTpxgxYgS6d++O2NhYuLq66hwrEZk3nZKhmJgYlWUC33//PbZv345XXnkFo0aNAlCw+XHHjh1Sn3HjxuHEiRNo2rQpmjVrhujoaAwZMgTh4eF6eggV1/MkAUUTjuLJkTr6TJiKK68ZJyZHVFHI5XKVfSSaPH78GDk5OdIYb29vWFhYIC+vZHXM3NxcWFqq/6GkUaNGOHz4MCIiIhAWFgYPDw8sWLAAISEhKkuFilbqDA4OxuzZs/Hee+/hxx9/rJQnAp07dy4iIiIQGRmJTp06qe0zceJEuLi4YNq0aYiPj8c///wjFUzIycnBpEmT4OjoiFdeecWQoZs0ubzgK0JZ7/G0tDSkpKRIlz09PaVqc8Xf46W9vy0tLbF7927Mnj0bY8aMQU5ODt58802EhIQgKytL6jd69Gjp/11dXbFy5Uq4ublh586dGDZsmE6PkYhIp2TIy8sLgwcPBgDpvwBQq1Yt6f8HDhyIFi1aSJetra2xY8cOXLhwAXfu3EHdunWl9cOku+KJQlmzMfpMmIp7nhknfSzHU4cJFBlDUFAQIiMjkZWVpbbEdaHw8HDExsYCAPz8/HDy5ElUr14djx8/Vlnylp6ejszMzFI3/rdt21alalZ0dMG/wcDAQI1jateuDaVSiTt37qB+/fo6PUZTd+7cOUyZMgX+/v6IiIhARESEdN3XX38tLTeMioqSziPk5OSEDRs24IMPPkCtWrVw7do1uLi4YPfu3XBz42dJoaCgIABAbGys2qWDhdauXYuvvvpKurxr1y7pVByPHj1S6ZuQkFDq+9vDwwOLFi1SaatTpw5effVVjWMcHR3h7u6OuLg4jX2IiDTRKRny9PTEgAEDSu0TFBQkfYAW1bRp0zLXcJPu1CUBRZMJbZbjaZMwFfe8M076mk0qTl0CxQSJylt4eDgiIiLw7bffYurUqSWuz8/Ph6WlJf78s2TVyJdffhmpqak4deoUXnrpJQDA3r17IZPJpL1DAHDv3j2VcwQplUqVc7ysWrUK/v7+Uunh4tcDBYmAtbX1c50rztTVrl0bf/zxh9rrin7p/vrrr6Vlbs7Ozti4cSMeP36Mq1evwtvbG7Vr19Y4Y2GugoKC0KZNG8ydOxcDBgwo8fwUvr/HjBmDMWPGlBhft25d7N+/X0pksrKy8Ndff6msMHn69Cny8/OlcwQVf/8ePXoUt27dwqBBgwBA5ceDQnFxcUhISEDt2rX188CJyKxUjmoEpELX2SN97F9Sx9QKPjA5In2rV68eFi1ahA8//BApKSkYOHAgnJyccPPmTWzbtg1+fn6YNm2a2rEtW7ZEr169MGLECCxduhTZ2dkYP3483n33XZWkxc/PD99++61USatXr16YNGkSfH19sXHjRixfvhy7du2SvqiuWLECly5dQv/+/eHh4YEDBw5gzpw5+PjjjytlAQVnZ2etSogXPU9NIU9PT3h6epZHWJXG6tWr0aVLF3Tr1g2TJ09GvXr18OjRI0RFRWHPnj04ePCgxrEzZszAe++9h2bNmqF58+b48ssv4ejoqHK+qzFjxiAuLk4qkDBlyhTUr18f7du3R2xsLMaOHYtRo0ZJhRiePXuGPn36YMKECahfvz5u376NSZMmITAwEK+//nr5PhlEVCkxGTID2iQBui5Ne54ZJ0MXfCirwAOTI9KVXC6Hr6+vykb6UaNGoUWLFvjhhx/w9ttvIzc3F3Xq1EG/fv3w9ttvl3p7GzduxPTp0zF8+HDI5XK89dZbmD59ukofX19fODg4SJdnzZqFmTNn4tatW2jcuDH+/vtvtG7dWrr+nXfewcqVKzFjxgw8fPgQtWrVwvLly6Vf1olK4+npqXIy08DAQFy4cAHfffcdZs6cicTERPj4+KBdu3ZYt25dqbcVHh4OpVKJJUuWICkpCc2bN8dff/0FJycnqY+bmxsyM//vuDB58mRMnToVc+fOhbu7O6ZPn46RI0dK17u6umLBggVYsGABLl68CDc3N4SFheGTTz4pdakqEZEmMiGEMHYQZUlNTYWTkxNSUlIqxXmGKgtdEyghRLntX1KHJ6g1XdnZ2bh9+zZq1aoFW1tbY4dj9kp7Pfj5q15pzwvf36aFrweRcZXXSVe7Na+ll2MTZ4bouek646Sv/Uv62q+kLjni7BERERGR+WAyROWqPPYv6Wu/UqCXA2aGBpU6W8TCDERERESVF5MhMih97F/S136lqwlpSM3Og43covhwCWePiIiIiCovJkNkcsojYSqaHCnylBi17iwASP/VRF2hBhZmICIiIqocmAxRhfQi+5Vs5BYI9HLA1YS0Mm9DXWU7bZbXEREREZHpYzJEZkcmk2FmaNBzF2oovryurJkigLNFmlSAYpZmga9D+eDzahr4OhBRaZgMUaVVWvGG59l3pGl5HavU6c7KygoAkJmZyXODmIDC87wUvi70Yvj+Ni05OTkAIJ2YmIioKCZDZDZetLKdpuV1ZZXwLhzL2aP/Y2lpCWdnZzx+/BgAYG9vz2WHRiCEQGZmJh4/fgxnZ2d+WdQTvr9Nh1KpRGJiIuzt7SGX8ysPEZXETwYyW+qSj9ISJHXL67Qp4Q083wlgKztvb28AkL4wkvE4OztLrwfpB9/fpsPCwgI1atQw+89cIlJPJirAYlqeAZ2MpazZIyFEmcmROtoUYTCXmaL8/Hzk5uYaOwyzZWVlVeqMED9/1dP2eeH72/isra1hYaH5FApEVL7K+i71PDLS0tCteS29HJs4M0RUirJmj0or4V1I23McmetskaWlJZdnUaXF9zcRkWljMkSkI10LM2h7jqPiS+lO3HiikhyZy0wRERERkaEwGSJ6QWUVZtD2HEfF9xrxfEZERERE5YvJEJEBaTrHkaaldIo8ZZklwImIiIjo+TAZItIzXfcZAZqX0hVNmoovm9N0X0RERESkHSZDRAagy1K6ooruK+KyOSIiIiL9YjJEZGI07StSt2yueFLFmSIiIiIi7TEZIjKC0pbSFd9XpGnZHGC+5biJiIiI9IHJEJGJKGspHVCyHHfxpXPqxnC2iIiIiEg9npKZyMQVLptTp3DpHBERERHpjjNDRCZOXTluTUvnuGyOiIiISHtMhohMlDbL5gDVpXP+bvaYFdpQ5friJbm5bI6IiIioAJMhogpIU8W5+KRMvLMqWqWNJbmJiIiI1GMyRFQBqVs6N2tXLOKTMkv0VVeSm4iIiIiYDBFVWMVP1BrRt7HGfUVF8dxERERERAWYDBFVEKWdmwgomRwVxSILRERERCUxGSIyA0VniHhuIiIiIqICTIaIKrDSKs5pKrLAPUREREREBZgMEVVSxYss6HJuIu4rIiIiInPAZIioEilrX1Gh0pbNEREREZkLJkNEZqK0ZXOp2XmwkVuo9C2aHHFfEREREVVGTIaIzERpy+aKl+DWZraIS+mIiIioorMouwsRVRaF5bdtrSzhaCtHoJeD2n6FRRaIiIiIKjPODBFVcpoqzhWfKQI0n6iViIiIqDJiMkRkxniiViIiIjJnTIaISK2iM0T+bvaYFdpQ5XptiiwUx31FREREZEqYDBGRRFPFufikTLyzKlqljSW5iYiIqKJjMkRkZjTtIQLU7yOatSsW8UmZJW6nsMiCpmV2RERERKaOyRARqSi+jyiib2ONRRa4r4iIiIgqMiZDRGZO3T6e4rNFmmZ/iu4r4rmJiIiIqKLheYaISCeF+4qK47mJiIiIqKLhzBAR6aT4viKem4iIiIgqKiZDRKQzTUvnypoZ4r4iIiIiMiVMhoiohNIqzpWmrBmi4vuK1N0u9xERERGRoei8Z2jNmjUICAiAjY0NGjVqhF27dpXaPysrC+PGjUONGjVga2uL2rVrY/r06cjPz3/uoInIdGjaQ6QO9xVRecnKysK5c+dw48YNKJXav8euXLmC06dPIzOzZPl4IiKq/HSaGdq7dy+GDx+O5cuXIzQ0FCtXrkS/fv0QFRWFli1bqh3z6aefYuvWrdixYwcaNmyIqKgohIWFoUqVKvj000/18iCIqHyVVnFO3bmJiuO+Iiov6enpmDJlCn799Vf4+/sjISEBDg4OWLNmDVq3bq1xXEJCAnr37o3bt2/Dw8MDDx8+xMqVK9G3b18DRk9ERMam08zQggULEBoairfeegsuLi4YP348mjVrhoULF2occ+7cOfTs2RMtW7aEra0tOnfujA4dOuDsWX4xIqosCvcQafqzkWv/UXPyZpLKH1FpkpKSEBgYiPv37+PcuXOIj49HmzZt0Ldv31JniEaOHAkLCwvcvXsXV65cwbRp0/Dmm2/iwYMHBoyeiIiMTetvKEIIREVFoXPnzirt//nPf3DixAmN4wYPHow///wT58+fh0KhwN9//42oqCgMGTLkuYMmIiICAH9/f3zwwQewtbUFAMjlcrz77rt48OAB4uLi1I558uQJdu/ejQkTJsDOzg4A8OGHH8LKygq//faboUInIiIToPUyubS0NGRkZMDDw0Ol3dPTE48ePdI47r///S9u3ryJFi1aAAAsLCzw1VdfoV+/fhrHKBQKKBQK6XJqaqq2YRKRgTxvkYWiy+lYXY7Kw+nTp2Fra4vq1aurvf7ChQtQKpUqy7ttbGzQsGFDnD9/Xu0YHpeIiCqnF64mJ4Qo9cvM5MmTsXXrVkRFRaFx48aIiorC66+/Djs7O3z44Ydqx0RERODzzz9/0dCIyAQV3TtUvLoc0Yu6dOkSZs+ejSlTpsDa2lptn2fPngEAXF1dVdrd3d2l64rjcYmIqHLSOhlycHBAlSpVkJiYqNKemJgILy8vtWOUSiUWLVqEiIgItGnTBgDQtWtXjBw5Et9++63GZGjKlCkYP368dDk1NRV+fn7ahkpEJqaw4tzVhDSV9sLqcurOWVTI1MtvP8++JlOKvzK5desWevbsidDQUHz22Wca+xUmSVlZWXBxcZHaMzMz4eTkpHYMj0tERJWT1smQTCbDSy+9hMOHD6skMYcOHUK7du00jrG0tCzxq68QApaWmr/82NjYwMbGRtvQiMjEFa84x+pypG+3bt1C586d0b59e6xduxYWFpq3xPr7+wMA7t+/Dx8fH6n9/v37aNKkidoxPC4REVVOOi2TmzhxIkJDQ7FmzRqEhoZi1apVOH/+PH788Uepz6xZs7Bw4UIkJydDJpOhT58++Prrr9GmTRs0adIEUVFR+OWXXzB8+HC9PxgiMl2FFef0QR9V5jg7U3nExcWhS5cuaNu2LdavX6/2x7bz58/D2toaDRs2ROPGjeHt7Y2dO3eiVatWAICrV6/i33//xTfffGPo8ImIyIh0SoZ69uyJ5cuXY86cORgxYgTq1auHrVu3ajzHEAAsWbIEM2bMwKBBg5CQkABfX1+MGTOm1CUMRETljWW7K4fHjx+jS5cuqFq1KkaOHIljx45J1zVv3lxa9jZ27Fh4e3tjy5YtsLCwQEREBEaOHAlHR0fUqFEDc+bMQZcuXdCzZ09jPRQiIjICmRBCGDuIsqSmpsLJyQkpKSlwdHQ0djhE9JwKE5Ds3Hy8syoaALB0aEuV8xCZQ4W5ijQrZeqfvzExMRg7dqza6xYtWoTGjRsDKEiGXF1dVYog7N69G2vWrEFaWhrat2+P8ePHw97eXqv7NfXnhYjIVJTHj48ZaWno1ryWXj6DX7iaHBHRiyi+d4gV5kgXjRs3xpEjR8rst2jRohJtvXv3Ru/evcshKiIiqii0Py08EZGeFFaXU6ewwhwRERFReePMEBEZXPHqcgArzBEREZmyyrrXlskQERlFadXliiZJ5rCHiIiIiIyDyRARmZyiM0TcQ0RERETlhXuGiMgkaNpHxD1EREREVF44M0REBlO8pHTR9cfF9xFxDxERERGVNyZDRGQySttHRERERKRvTIaIyOQVXybHogpERESkD0yGiMjk8cSsREREVB5YQIGITBJPzEpERETljTNDRGSSeGJWIiIiKm9MhojIZGl7YtbnwX1HRERExGSIiIymeKltQLXcdmledIaI+46IiIiIyRARVRiF+4iuJqS98G0V7jsyRinv4gmfuqSQiIiIyh+TISKqMNTtI9IV9x0RERFRISZDRFSh8MSsREREpC8srU1ERERERGaJyRAREREREZklLpMjIpNSvJiAttXliIiIiHTFmSEiIiIiIjJLTIaIiIiIiMgsMRkiIiIiIiKzxGSIiIiIiIjMEgsoEJHZepGTtxaykVtAJpPpIRoiIiIyNCZDRGS2Rq07+8K3EejlgJmhQUyIiIiIKiAmQ0Rk0oqX2gZerNy2jdwCgV4OuJqQ9iJhSa4mpEGRp4StlaVebo+IiIgMh8kQEZkVmUyGmaFBL7xETpGn1MvMEhERERkPkyEiMjsymYwzOURERMRqckREREREZJ44M0REFY66fUS6epF9R0RERFQ5cGaIiIiIiIjMEpMhIiIiIiIyS0yGiIiIiIjILDEZIiIiIiIis8RkiIiIiIiIzBKTISIiIiIiMktMhoiIiIiIyCwxGSIiIiIiIrPEk64SERmZuhPA6uPEskRERFQ6zgwREREREZFZYjJERERERERmickQERERERGZJSZDRERERERklpgMERERERGRWWIyREREREREZonJEBERERERmaXnOs/QxYsXER8fj3r16qF+/fpajcnKysKpU6cgk8nQpk0b2NraPs9dExERERER6YVOyVBOTg4GDhyIo0ePolmzZoiOjsbrr7+OZcuWQSaTaRy3detWjBw5ErVq1YKHhwfu37+P3377DQ0aNHjhB0BEROZNCIEDBw5g2bJluH//Pnbu3AlXV9dSx4wdOxbnz59XaevSpQvmzJlTnqESEZGJ0SkZWrhwIY4fP44LFy6gevXqiI2NRXBwMDp37ozw8HC1Y86ePYvBgwfjhx9+wMiRIwEAcXFxePr06YtHT0REZm/w4MFISkpCs2bNsGnTJuTk5JQ5JiYmBnXq1MGIESOkNnd39/IMk4iITJBOydDatWsxaNAgVK9eHQDQsGFDhISEYO3atRqToXnz5qFZs2ZSIgQANWvWRM2aNZ8/aiIiov9v6dKlcHFxwYEDB/D1119rPc7Pzw8dOnQox8iIiMjUaZ0M5eXl4cqVKxg7dqxKe5MmTbB06VKN4/766y+MHDkSjx8/RnR0NDw9PdG0aVNYW1trHKNQKKBQKKTLqamp2oZJRERmxsXF5bnG7dixAydPnkS1atUQEhKCN998U8+RERFVXCdvJhk7BIPQOhlKT09Hfn5+iYOOm5sbkpOT1Y4RQiAxMRExMTFo1aoVgoKCcO3aNVhYWGD79u1o2LCh2nERERH4/PPPtX8UREREOnB3d0f37t3RokULxMbG4sMPP8ShQ4ewfPlytf35Ix0RUeWkdTJkY2MDAMjMzFRpT09P11gZTiaTwdraGlFRUYiJiYGHhwfy8vLw6quvYuTIkTh+/LjacVOmTMH48eOly6mpqfDz89M2VCKiMrWt4/ZC4zNz8vQUCRnDunXrpGNXz549Ub9+ffTu3RsffvghmjZtWqI/f6QjItKdEAKKPKXebzc7L19vt6V1MmRnZwdvb2/cuXNHpf3OnTuoXbu2xnF16tRBkyZN4OHhUXCHcjn69u2Ljz76CEIItVXobGxspOSLiIhI34r/iNetWzfIZDL8888/apMh/khHRKQbIQRm7YrFtYR0vd+2UpFZdict6XTS1ZCQEPz+++9QKgsyvOzsbOzatQshISFSn8uXL2Pnzp3S5d69e+PatWsqt3Pt2jX4+vqWWo6biIjIUB4/fgwhBOzt7dVeb2NjA0dHR5U/IiLSTJGnLJdESN90qiY3Y8YMtGrVCv3790evXr2wceNGyOVylV/LNm3ahIULF0r7iD755BNs2rQJgwYNwquvvorY2FgsWbIEK1as0OsDISIi0mTMmDFwc3PDnDlzcOvWLZw8eVIqmJCZmYnx48fD1dUV3bp1M3KkRESVz9KhLWEj12kOplQZ6WkIXaif29Ipqpo1a+LcuXOoX78+jhw5go4dOyI6Ohpubv+39j4oKAhhYWHSZXd3d5w5cwaNGjXCvn37oFQqcfz4cbzxxhv6eQREREamyFMiOzdf6z8hhLFDrlSWLl2KDh06YNy4cQCAsLAwdOjQAQcPHpT6XLx4EVeuXAEAeHh44PDhw3B3d0eLFi3g4+ODGzduYO/evXB2djbGQyAiqtRs5BawtbLU35/cUm+xyUQFOCqnpqbCyckJKSkpXJpARCYhMycPQTP+fK6xgV4OmBkaVOpS4Rct8KAvFeHzNz4+Hnfv3i3RXr9+felEqjExMbCyskL9+vWl69PS0nDjxg14e3ujWrVqOt1nRXheiIhexIuW1s7Ozcc7q6IBACuHtYKtlf4SmIy0NHRrXksvn8E6LZMjIqICdlaWCPZ3wZn4ZzqPvZqQBkWeUq8HBnPm7+8Pf3//Uvs0bty4RJuDgwOaN29eXmEREVEFwGSIiOg5yGQybB7VFkeuJmo9RpGnxKh1Z8sxKiIiItIFkyEiouckk8k4u0NERFSB6a+sAxERERERUQXCZIiIiIiIiMwSkyEiIiIiIjJLTIaIiIiIiMgsMRkiIiIiIiKzxGSIiIiIiIjMEpMhIiIiIiIySzzPEBGRCTp5M0nnMW3ruJVDJERERJUXZ4aIiIiIiMgsMRkiIiIiIiKzxGSIiIiIiIjMEpMhIiIiIiIyS0yGiIiIiIjILDEZIiIiIiIis8RkiIiIiIiIzBKTISIiIiIiMktMhoiIiIiIyCwxGSIiIiIiIrPEZIiIiIiIiMwSkyEiIiIiIjJLTIaIiIiIiMgsMRkiIiIiIiKzxGSIiIiIiIjMEpMhIiIiIiIyS3JjB0BEZI4UeUrp/23kFpDJZEaMhoiIyDwxGSIiMoJR685K/x/o5YCZoUFMiIiIiAyMy+SIiAzERm6BQC+HEu1XE9JUZoqIiIjIMDgzRERkIDKZDDNDg6TER5GnVJkhIiIiIsNiMkREZEAymQy2VpbGDoOIiEjFyZtJxg7BKLhMjoiIiIiIzBJnhoiIiIiISCMhhM57WyvKXlgmQ0REREREpJYQArN2xeJaQrqxQykXXCZHRERERERqKfKUL5QIBXo5wEZuuikHZ4aIiIiIiKhMS4e21DmxMfUTizMZIiIiIiKiMtnILSpdRVTTnbMiIiIiIiIqR0yGiIiIiIjILDEZIiIiIiIis8RkiIiIiIiIzBKTISIiIiIiMktMhoiIiIiIyCwxGSIiIiIiIrPEZIiIiIiIiMwST7pKRFRJnLyZpPfbzEhL0/ttEhERmQqdZ4bu3r2LqVOn4s0338SsWbPw5MkTrcdGR0djwIAB+Prrr3W9WyIiolKdP38eu3fvhkKh0Kq/UqlEdHQ0Dh06hGfPnpVzdEREZIp0Sobi4uLQokULXLlyBZ07d8bRo0cRHByMpKSyf41MTU3Fm2++idOnT+PkyZPPHTAREVFRW7ZsQcuWLREaGorQ0FCtEpt79+6hcePG6NevHyZPngw/Pz+sX7/eANESEZEp0SkZmjNnDqpXr46tW7dixIgR2LNnD/Ly8vDtt9+WOfb999/H66+/jiZNmjx3sERERMU9efIEP/30E1atWqX1mBEjRsDV1RU3b95EdHQ05s6di+HDh+Pu3bvlFygREZkcnZKhyMhI9OvXDxYWBcNsbW0RGhqKyMjIUsctW7YMN27cwKxZs547UCIiInVGjRqF4OBgrfsnJCTgzz//xLhx42BtbQ2g4Ac7Ozs7/Pbbb+UVJhERmSCtCyhkZWXh0aNHqFGjhkp7jRo1sGHDBo3jrly5gqlTp+LYsWOwsrLS6r4UCoXKmu/U1FRtwyQiIipVTEwMhBBo2rSp1GZlZYWgoCBcuHBB7Rgel4iIKietZ4YKDwL29vYq7VWrVkV2drbaMdnZ2Rg0aBC+/PJLBAQEaB1UREQEnJycpD8/Pz+txxIREZUmOTkZAODq6qrS7ubmJl1XHI9LRESVk9bJUNWqVWFpaVliY2pSUhKcnZ3Vjtm3bx+uXr2KP//8EwMGDMCAAQNw9uxZREVFYcCAAUhMTFQ7bsqUKUhJSZH+uIabiIj0pXBpXGZmpkp7ZmamdF1xPC4REVVOWi+Tk8vlaNCgAS5evKjSfvHiRY1FEYKDg0tU54mLi4OVlRUGDx6MKlWqqB1nY2MDGxsbbUMjIiLSWq1atQAUVJTz9fWV2u/evYtmzZqpHcPjEhFR5aRTAYWhQ4di06ZNuHfvHgAgNjYWf/zxB4YOHSr12bRpE95++20AgI+PjzQjVPjn7e0NX19fDBgwoMSSOyIic6XIUyI7N1/6E0IYO6RK5cSJEzh79iwAoFGjRqhevTp+//136fqYmBhcu3YNPXv2NFaIRERkBFrPDAHAuHHjcOLECTRt2hTNmjVDdHQ0hgwZgvDwcKnP5cuXsWPHDr0HSkRUmY1ad1blcqCXA2aGBkEmkxkpoorj33//xY0bN6TiBwcOHICzszNatGgBHx8fAMCkSZPg7e2NLVu2QCaTYcGCBRg6dChsbW1Ro0YNzJs3D7169cIrr7xizIdCREQGplMyZG1tjR07duDChQu4c+cO6tatiwYNGqj0GThwIFq0aKHxNqZMmQJLS8vni5aIqBKxkVsg0MsBVxPSSlx3NSENijwlbK34eVmWU6dOYfPmzQCAV199FRs3bgQATJ48WUqG2rdvDxcXF2nMoEGD4OnpibVr1+Ly5csYM2YMxowZY/jgiYjIqGSiAqzFSE1NhZOTE1JSUuDo6GjscIiIJCdvJr3QeCEEFHlK6bIiTynNEq0c1sroyVBGWhq6Na/Fz99ieFwiospG0/EsOzcf76yKBmAaxyVAv8cmnWaGiIhIv2QymUkcWIiIiMyRTgUUiIiIiIiIKgsmQ0REREREZJa4TI6IiIiIiCRF97MW3ddaGTEZIiIiIiIiAAWJ0KxdsbiWkG7sUAyCy+SIiIiIiAhAwUyQukQo0MsBNvLKlzpwZoiI6AW0reOm85gXLcdNRERkCEuHtpQSIBu5RaU8ETiTISIiIiKiSkKfP7jZyC0q/ekfKt9cFxERERERkRaYDBERERERkVliMkRERERERGaJyRAREREREZklJkNERERERGSWmAwREREREZFZYjJERERERERmickQERERERGZJSZDRERERERklpgMERERERGRWWIyREREREREZonJEBERERERmSUmQ0REREREZJaYDBERERERkVliMkRERERERGaJyRAREREREZklubEDICIi9RR5Sp3628gtIJPJyikaIiKiyofJEBGRiRq17qxO/QO9HDAzNIgJERERkZa4TI6IyITYyC0Q6OXwXGOvJqTpPJtERERkzjgzRERkQmQyGWaGBumU1CjylDrPIhERERUSQkjHHXP7UY3JEBGRiZHJZLC1sjR2GEREZAaEEJi1KxbXEtKNHYpRcJkcEREREZGZUuQp1SZCgV4OsJFX/lSBM0NERERERISlQ1tKCZC5VChlMkRERERERLCRW5jdMu3KP/dFRERERESkBpMhIiIiIiIyS0yGiIiIiIjILDEZIiIiIiIis8RkiIiIiIiIzBKTISIiIiIiMktMhoiIiIiIyCwxGSIiIiIiIrPEZIiIiIiIiMwSkyEiIiIiIjJLTIaIiIiIiMgsyY0dABGRuWlbx03l8smbSUaKhIiIyLxxZoiIiIiIiMwSkyEiIiIiIjJLXCZHRERERFQJCCGQnZuv0xhFnrKcoqkYdE6G1qxZgy+++ALx8fGoV68eIiIiEBoaqrH/uXPnMG/ePPz111/Izc1Fq1atEBERgebNm79Q4ERERIXWr1+PNWvWIC0tDe3bt8e0adPg7Oyssf97772HM2fOqLR169YN//vf/8o5UiKi8iGEwIClJ3E2/pmxQ6lQdFomt3fvXgwfPhyfffYZHj16hHfffRf9+vXD2bNnNY6ZMWMGBgwYgAsXLuDff/+Fn58f/vOf/+D+/fsvHDwREdGSJUswYsQI9O/fHzNnzsRff/2FkJAQKJWaf+28ceMGWrVqhVWrVkl///3vfw0YNRGRfmXl5r9QIhTo5QAbufntoNFpZmjBggUIDQ3FW2+9BQAYP348NmzYgIULF2Lt2rVqx+zevVvl8uLFi7Fq1Srs378fw4YNe76oiYiIAOTl5WHmzJmYOnUqRo4cCQAICAhA7dq1sXv3brz22msax3p5eaFZs2YGipSIyHCWDm2pc2JjI7eATCYrp4hMl9bPkhACUVFR6Ny5s0r7f/7zH5w4cULrO3z27Bny8/Ph5OSk9RgiIiJ1Ll68iCdPnqBXr15SW61atdCwYUMcPHiw1LGbN29GmzZt0KdPH/z000+lziQREVUkNnIL2FpZ6vRnjokQoMPMUFpaGjIyMuDh4aHS7unpiUePHml9hx999BH8/PzQo0cPjX0UCgUUCoV0OTU1VevbJyIi83H37l0AQLVq1VTafXx8pOvU8fHxQb9+/dCiRQvExsZi5syZOHLkCDZs2KC2P49LRESV0wtXkxNCaJ1JTps2DZGRkTh06BDs7e019ouIiMDnn3/+oqEREZkdfVcFys7TrSqRoeXm5gIAbGxsVNptbGyk69RZvXo1rKysAAAdOnRA7dq10b17d0ycOBEtW7Ys0Z/HJSIyhBc5CbeuVeSogNbJkIODA6pUqYLExESV9sTERHh5eZU5ftasWfj+++/xxx9/oFWrVqX2nTJlCsaPHy9dTk1NhZ+fn7ahEhGZrVHrNBe0eR5KRaZeb0/f3NzcAABJSUlwdXWV2pOSkhAQEKBxXGEiVKhz586QyWSIiYlRmwzxuEREVDlpvWdIJpPhpZdewuHDh1XaDx06hHbt2pU69vPPP8fXX3+NyMhIdOjQocz7srGxgaOjo8ofERGpZyO3QKCXg7HDMIpmzZrB0tISp0+fltqysrJw8eJFtUmNJg8ePIAQAg4O6p9HHpeIiConnZbJTZw4EaGhoVizZg1CQ0OxatUqnD9/Hj/++KPUZ9asWVi4cCGSk5MBAF9++SUWLFiAyMhIdOzYUa/BExFRwY9VM0ODyuXEeRnpaQhdqPeb1RsXFxcMGDAA//vf/9C7d284OTkhIiICFhYWGDx4sNTvnXfegbu7O/73v//hxo0bOHjwIIYPHw65XI6UlBR88MEH8PT0RLdu3Yz4aIiIyNB0SoZ69uyJ5cuXY86cORgxYgTq1auHrVu3lvrr25w5c6BQKPDyyy+rtM+cOROzZs16rqCJiEiVTCaDrZWl3m83X67/29S3JUuWYODAgahWrZpUqXTLli0qBX+uX7+OtLQ0AAXFFq5cuQJ3d3d4enri7t27CA4OxoEDBzjjQ0RkZmRCCGHsIMqSmpoKJycnpKSk8EBFRJXOi2yYLW8ZaWno1rxWhfj8ffDgAdLS0lCnTh3I5aq/9d24cQNyuRw1a9aU2nJzcxEXFwdvb2+Ny+M04XGJiMrDixZQeGdVNABg5bBW5fID2fNoW8dN77epz8/gF64mR0REZAp8fHw0Xle3bt0SbVZWVqhXr155hkRERCZOt1PTEhERERERVRJMhoiIiIiIyCwxGSIiIiIiIrPEZIiIiIiIiMwSkyEiIiIiIjJLTIaIiIiIiMgssbQ2EREREdELMOXzxVHpmAwRERERERmZEAKKPOVzj3+RseaMyRARERERkREJITBrVyyuJaQbOxSzw2SIiIiIiEgH+l4Wp8hT6i0RCvRygI1c/2UB2tZx0/ttmgImQ0REREREJmLp0JYvlMzYyC0gk8n0GFHlxmSIiIiIiMhE2MgtYGtlaewwzAaTISIiIiIiDVgprnJjMkRERERE9JxetAocwEpwxsRkiIiIiIjoObAKXMXHZIiIiIiISEtFZ4L0WQUOKL9KcKQZkyEiIiNTV66Ua9SJiIyjtM/f0maCXrQKHMBKcMbAZIiIiIiISAuaZoICvRzgaCtnIlMBMRkiIiIiIlKjeHGEov9fdCaIMzoVF5MhIiIiIjJbmpbFlVUcobKfD0jdEu7KiMkQERERERG0L47AQgeVB5MhIiIiIjJ7uhRH4LK4yoPJEBERERGZPRZHME9MhoiIiIioUirrNAXFl8UVYnEE88FkiIiIiIjMTmnL4nQtjmBKxQZ4njrdMBkiIiIiIrNT2rK4soojmFLyQy+GyRARERERVXjFZ0SKnyOoOC6LU2WuCR6TISIiIiKqcEpbDlbWOYKKK21ZnLkmCeaCyRARERERmbyykp+iMz2lnSOouOLL4pj8mBcmQ0RERERUoRRPfmbtikV8UqbavsXPEVScuS6LowJMhoiIiIjIpJS1/6e05Kcobc4RVNlmgirb4ylvTIaIiIiIyKiKJj/qCh9ok/z4u9ljVmhDlTbO+lBZmAwRERERkdEIIZCdmy9d1nbWp3jyoy7x4SwJlYXJEBEREREZTPFZIG2rvmmT/BDpiskQERERERlE8VkgTVXfnmfJG2eB6HkwGSIiIiKiciGEwJGridJlbau+cckbGQqTISIiIiLSC12Sn6KKV31j4kOGwmSIiIiIiHRWPPEBnr/qW+dAD+7/IaNgMkREREREZXreWZ/iyQ8THzIlTIaIiIiIqFRCCAxYehJn45+V2besqm9MhMiUMBkiIiIiMnPqlrwVpchTqk2EtKn6xv0/ZMqYDBERERFVYmUlOoD2S94A1apvXPJGFR2TISIiE/Q8v6QWPZEhEZmv593bo43iVd+YCFFFx2SIiIiIqILQ9yxPUeqWvBXHmSCqbJgMERERERmBEAJZufkqbaduPS11DBMdIv1iMkRERESkZ9okOvpcvlYcy1kTaee5kqGLFy8iPj4e9erVQ/369cttDBERkTZyc3MRFRWFtLQ0BAcHw9PTs1zGUOWjLmlRp6wZm+IMmei0qe1aoo+dlSWTHyIt6JQM5eTkYODAgTh69CiaNWuG6OhovP7661i2bJnGf3DPM4aIiEhbcXFx6N69O/Lz8+Hj44Nz585h0aJFePfdd/U6hgxL2ySlLOW17ExftFm+VjzZYaJDpD86JUMLFy7E8ePHceHCBVSvXh2xsbEIDg5G586dER4errcxRERE2ho5ciR8fX2xf/9+yOVy/Pzzzxg9ejS6du2KmjVr6m0MaUdfMy3GTlL0gft0iEyfTAghtO3cuHFjdOrUCYsXL5ba+vXrh/T0dOzbt09vY4pLTU2Fk5MTUlJS4OjoqG24RERmpTxKa2ekpaFb81om+/n76NEj+Pj44Pfff0efPn0AAHl5efDy8sLkyZMxadIkvYwprvC49DAxySSfF0MrmthUxCRGm6SlOHVL04rjDA5R+dBnbqD1zFBeXh6uXLmCsWPHqrQ3adIES5cu1dsYAFAoFFAoFNLllJQUAAUPnIiI1MtIS9P/baYX3KYOv5sZVExMDIQQaNKkidQml8vRoEEDXLx4UW9jNB2XgmfsgoWNvT4eCqnh52qHKSENXug2gmu6lNnneZKWvOyyE760bJ1ukoi0VJgT6OPYpHUylJ6ejvz8fLi4qH6ouLm5ITk5WW9jACAiIgKff/55iXY/Pz9twyUiIj1KSkqCk5OTscMooTAp0eU48zxjNB2X7v84TLeASSd3AYTOMHYURGSq9HFs0joZsrGxAQBkZqr+EpKeng5bW1u9jQGAKVOmYPz48dLl5ORk+Pv7486dOyZ5MC4uNTUVfn5+uHv3boVYPsF4yxfjLV+Mt3ylpKSgRo0acHUte0mQMRQ9zhRNbtLT00skOy8yhsclw6toMTPe8sV4y1dFi1efxyatkyE7Ozt4e3vjzp07Ku137txB7dq19TYGKDhQFR6sinJycqoQL1AhR0dHxluOGG/5Yrzlq6LFa2FhYewQ1Co8lty5cwe+vr5S+507d9CyZUu9jeFxyXgqWsyMt3wx3vJV0eLVx7FJp1sICQnB77//DqVSCQDIzs7Grl27EBISIvW5fPkydu7cqdMYIiKi5xEUFAR/f39s2bJFajt//jxu3Lihcpw5cuQIoqKidBpDRESVn06ltWfMmIFWrVqhf//+6NWrFzZu3Ai5XK6ydGDTpk1YuHChtO5amzFERETPQyaT4ZtvvsGgQYMgl8tRo0YNfP311+jTpw+6dOki9fvss8/g7e2NLVu2aD2GiIgqP51mhmrWrIlz586hfv36OHLkCDp27Ijo6Gi4ublJfYKCghAWFqbTmLLY2Nhg5syZapcomCLGW74Yb/livOWL8epfv3798NdffyEtLQ1RUVGYNGkSNm/erNKnS5cuaNu2rU5jSlMRnpeiKlq8QMWLmfGWL8Zbvsw5Xp3OM0RERERERFRZmOaOWCIiIiIionLGZIiIiIiIiMwSkyEiIiIiIjJLOlWTM4abN28iJiYGnp6eeOmll0zuXBfHjh3DvXv3EBYWBjs7uxLX5+Tk4MSJE0hJSUGrVq3g4+NjhCgLJCQk4Pz587C3t0fz5s3h4OBQok9eXh5OnDiBZ8+eoUWLFvDz8zNCpAXS09Nx5swZZGdno1GjRqhevXqJPvn5+Th58iSePHmCZs2aoWbNmoYPtJi7d+/i+PHjCAoKQpMmTVSuy8/Px6lTp/D48WM0adKk1PNtlaczZ87gxo0bKm3Ozs7o2bOnSptSqcSpU6eQkJCARo0aoW7duoYMs4SUlBRERUXBxsYG7dq1g7W1tcr1QgicPn0aDx8+RFBQEAICAowS586dO0ucbBoA/Pz80L59e+myEAJnzpzBvXv30KBBA9SvX9+QYarIzc1FdHQ0Hj9+DB8fHwQHB6v9vD137hzi4+NRr149NGrUyAiRGp8pvW7qKBQK7Nu3D5aWlujVq5faPs+ePcOJEycgl8vRoUMHVKlSxcBR/p+rV6/i2rVrqFatGlq0aKH2fZecnIzjx4/D0tISHTp0QNWqVY0QaYEHDx7gn3/+QZUqVdCiRQu1x9LU1FQcP34cANC+fXuTOG9LdHQ0bt68ia5du8LT01PluvT0dBw7dgz5+flo3749nJ2djRLjnj17kJaWptIWGBiI5s2bq7RlZmbi2LFjyMnJQbt27Yx+Uug7d+7gn3/+QY0aNdCsWbMS12dlZeHYsWPIzs5G27Zt4e7ubvAY09LSsGfPHrXXtWrVCnXq1JEuZ2dn4/jx48jIyMBLL71U4v1iSE+fPsX58+eRmZmJwMBAtcd1hUKB48ePIz09Ha1bt4a3t7dudyJM2LRp04S9vb145ZVXhK+vr2jVqpV4+vSpscMSQgixfv160aBBA1G3bl0BQNy9e7dEn9u3b4u6deuKunXris6dOws7Ozvx448/GjzWjIwMMXToUOHr6yt69OghgoODhaurq9i2bZtKv/v374ugoCBRq1Yt0bVrV2FnZycWLFhg8HiFEGLhwoWiRo0a4pVXXhH/+c9/hJ2dnfjwww9V+iQkJIgmTZoIf39/qc9XX31llHgLKRQKERwcLKytrcXkyZNVrnvy5Ilo2bKlqF69unjllVeEvb29mDFjhlHifP/990X16tXFoEGDpL+JEyeq9ElOThYvvfSS8PHxkeIt/pgM6eeffxYODg6iQ4cOIiQkRDRv3lzExcVJ16empoqOHTuKatWqiW7duokqVaqIjz/+2CixjhkzRuW57devnwAgxo0bJ/XJyMgQXbt2FV5eXqJ79+7CwcFBjB492ijxxsTECD8/PxEQECDCwsJEjRo1RGBgoIiPj5f6ZGdni1dffVW4u7uL7t27CycnJ/HWW2+J/Px8o8RsLKb0uqkzadIk4ePjI72G6uzatUs4ODiItm3biqZNmwpPT09x+vRpA0cqxKVLl0Tbtm1FYGCg6N27t/D39xeNGjUSN2/eVOn3559/CicnJ9GmTRvRvHlz4ebmJo4dO2bweLOzs8XQoUNFzZo1Re/evUXLli2Fk5OT+O2331T6HT58WLi4uIjg4GDRqlUr4eLiIg4fPmzweIu6evWqcHV1FQBKxHLs2DHh5uYmWrRoIdq0aSOcnJzEn3/+aZQ469SpI4KDg1U+P1evXq3S58yZM8LLy0s0adJEtGvXTjg4OIjt27cbJd7c3FwxYsQIUaVKFdGrVy/Rrl07ERYWJpRKpdTnwoULwsfHRzRs2FC0b99eVKlSRWzatMngsd6/f1/leR00aJBo166dACD27Nkj9bt8+bLw8/MTDRo0EB07dhT29vZi3bp1Bo9XCCFWrFghqlSpItq3by969+4tqlatKl5//XWRm5sr9bl+/bqoWbOmCAwMFC+//LKwt7cXy5cv1+l+TDYZOnz4sAAg/v77byGEECkpKSIgIECMGjXKyJEVWL16tYiNjZXiVJcM9ejRQ3Tu3Fl60VauXCnkcrm4fv26QWN9+vSpWLt2rcjLy5PaZsyYIezt7UVqaqrU1r9/f9G6dWuRnZ0thBBi8+bNQiaTiYsXLxo0XiGE2LRpk8jIyJAu//333yrvByGEePPNN0Xz5s1FZmamEKLgAA9AREdHGzzeQuPGjRPDhw8XgYGBJRKH4cOHi4YNG4q0tDQhhBD79+8XAMTRo0cNHuf7778v+vfvX2qf//73vyIgIEAkJycLIYQ4evSoACAOHDhgiBBV/Pnnn0Imk6kc8K5evSquXLkiXZ4wYYKoVauWSEpKEkIIcerUKWFhYSF2795t8HiL27x5swAgYmJipLZp06aJ6tWri8ePHwshhDh37pyQy+Viy5YtBo/vtddeEx06dJASm+zsbBEYGChGjhwp9YmIiBCenp7i/v37QoiCA6adnZ1YtWqVweM1JlN63dRZsGCBSExMFNOmTVObDKWkpAgXFxcxc+ZMqW3o0KEiICBA5QucIZw8eVKcPHlSuqxQKETHjh1Ft27dpLaMjAzh4eGh8nk6fPhwUbNmTZUvRIaQlpYmtm7dqvI8TZ06VTg4OEjH1+zsbOHj4yM++ugjqc9///tf4evrKxQKhUHjLZSdnS2aN28u5s+fXyIZys3NFf7+/uL999+X2iZMmCA8PT1VjsGGUqdOnVJ/NFYqlaJ+/frizTfflNpmzpwpnJ2dpWOVIX3yySfCy8tL3Lp1S2rbtWuXyo9ETZs2FQMGDJDeNxEREaJq1aoiMTHR4PEW9/777wtfX1+V74dt2rQRoaGh0mP49ttvhZ2dnXjw4IFBY8vNzRX29vZi9uzZUtuVK1cEAPH7779LbZ06dRLdu3eXHsOPP/4orK2tVX7MK4vJJkPDhw8XrVq1UmmbP3++cHR0NKlfIjUlQwkJCUImk4mtW7dKbXl5ecLDw0N88cUXhg6zhAsXLggA4uzZs0KIgl/V5XJ5iV9g/P39xaeffmqMEFUkJCSo/HqRlZUlbGxsxE8//aTSLyAgwGizAXv27BH16tUTaWlpJZKh3NxcUaVKFfHdd9+pjGnSpIlREvz3339fdO3aVezYsUMcOXJEPHv2TOV6pVIpnJ2dxbx581TaW7duLYYNG2bASAt06tRJ9OrVq9Q+Xl5e4vPPP1dpe/nll8XgwYPLMzSt9OzZU7Rt21alzd/fX0yZMkWlrXv37qJPnz6GDE0IUfDDzRtvvKHS1qVLF5XXumHDhmLs2LEqffr16ydeeeUVg8RoKkzpdSuNpmRow4YNwtLSUvrRQIiChA6AiIqKMmSIas2bN094eHhIl7dt2yZkMpl4+PCh1Hb58mUBQBw5csQYIapYvXq1sLKyEllZWUIIIf744w8BQGXW+tatWwKA2Lt3r1FiHDt2rBg2bJi4fft2iWSo8DvM1atXpbZ79+6V+PHJUOrUqSMmTpwofv/9dxEdHV0igYyOjlb57iKEEElJSUIulxt89iI1NVXY2dmJb775RmOfmJgYAUAcP35cZZyNjY1YtmyZIcLUKCMjQzg6Oorp06dLbdevXxcAxMGDB6W2zMxMUaVKFbF48WKDxpednS2sra3FmjVrVGKWy+Vi48aNQggh7t69W2JmKycnRzg7O+u0sslk9wzFxMSUWI/euHFjpKam4s6dOyaxN6Q0sbGxEEKoPAZLS0sEBQUhJibGiJEVOHDgAGxsbFCvXj0AwL///ou8vLwSz3mjRo2MFu+dO3ek/UurVq3CkCFDpD0t169fh0KhUPseMUa8Dx48wHvvvYdt27apXct++/ZtZGRkmEy8AHDhwgUsXboU9+/fR1xcHL777jsMGzYMAHDv3j0kJyerjfeff/4xaJx5eXk4efIk5s+fjzt37uD8+fPS3gK5vOAjLDExUdrXVDzeI0eOGDTe4u7du4d9+/Zh2bJlUltaWhri4+PVxrt9+3YDRwh89dVXGDhwICZPnoygoCBpn9gvv/wCoOA1uHLlCj788MMS8S5dutTg8RqLqb1uzyMmJgY+Pj4qeywaN24MmUyGmJgYtGnTxojRAQcPHlR5fmNiYuDu7q6yB6BBgwawsrJCTEwMOnXqZPAYT506hVu3biEuLg4///wzvvvuO9ja2krxOjg4wN/fX+pfq1YtVKlSBTExMejRo4dBY925cyciIyNx/vx5JCUllbg+JiYG1tbWKvswfH194erqipiYGISFhRkyXADAtm3bcOXKFVy4cAF2dnbYsGEDWrZsKcULQOU94urqCl9fX4MfS6Ojo5GVlYXu3bvj3LlzuHfvHurVq4cGDRpIfdTF6+DggJo1axr9u+CWLVuQnp6Od999V2pTF6+dnR3q1q1r8HhtbGywaNEizJ49G4mJiXBycsKvv/6KsLAw9OvXT2O8VlZWCAwM1Clek02GUlJSSmyIc3NzA1CwkdLUpaSkAIDax2Ds+C9cuIAZM2Zg2rRp0sbP0uK9efOmwWMEgPv372P79u1ISEjAnTt3MGjQIGljbWnxxsXFGTROpVKJoUOHYvTo0Rq/SJQWr6GTCwAYPHgwvv32W6nox8KFCzFixAi0bNkSjRs3Nqn3b3JyMnJycvDXX3/h22+/RZMmTRATEwM7Ozvs3r0btWvXNql4i1u5ciWqVq2KgQMHSm2mFq+Pjw+aNm2K33//HVeuXMHFixfRpUsXaZNveno6lEqlycRrLKb2uj0PdcdWuVwOBwcHoz+GH374AYcOHcJff/0ltamLFwBcXFyMFu+5c+dw6NAhXL9+HVWrVlVJJDTFa4z3yL179zBixAhs374dDg4OapOhlJQUuLi4lGg31nt60aJFCAkJAVBQgGrw4MEYOHAgrly5Amtra6SkpKBKlSoliucYI97Hjx8DAGbPno3Lly/D398fx44dQ5cuXbBp0yZYWVkhJSUFlpaWJQpomMJnxvLly9GtWzeVyQVT+4xr0KABqlatis2bN8PFxQU3btzA2LFjYWVlBUB/8ZpWabYibGxskJ6ertJWeLnwFxhTZmNjAwBqH4Mx47969Sp69uyJ119/HZ999pnUborxtm3bFhs3bsThw4exa9cufPrpp9i8eTMA04p3zZo1uHjxImrXro2NGzdi48aNSEtLw5UrV7Bp0yaTixcAOnfurFL98OOPP4aTkxP++OMPAKYVb+H9Xbp0CZcuXcLOnTtx9epVODo64qOPPjK5eIsSQmDlypV44403VKp1mVq8AwcORHZ2Nv7991/p+b18+TJGjRplkvEaS2V4HtQdW4UQyMzMNOpj2LBhA8aNG4eVK1eiXbt2Uru6eAHjPuejR4/G5s2b8c8//yA8PByvvfYaHj16BMC04p04cSIaNGiA+Ph4bNy4Ebt27QIAHDp0CEePHjW5eAFIiRAAWFtbY/r06bh16xYuXboEoCDerKwsKJVKlXHGPDbZ2tri4sWL2L17Ny5cuICDBw9iyZIlUrz5+fnIzs42erxFXb9+HX///TdGjBih0m5Kn3GJiYkICQnBW2+9hZMnTyIyMhIHDx7EzJkzsXr1ar3Ga7LJUJ06dXDnzh2Vtvj4eFhaWqpMP5uqwhKF6h6DscopX7t2DV27dkW3bt2wfPlyyGQy6TpTjLeo4OBgNGjQQPoAL4zJFOJ1d3fHK6+8gl27dmH79u3Yvn070tLScP36dezYsQNAwTIJmUxmEvFq4uDggMTERABAjRo1IJfLTSLeqlWrwsvLCz169JCWIFpbWyM0NBRnz54FUDCzYWtraxLxFnXo0CHcvn27xAHH3d0djo6OJhFvfn4+jh8/jn79+sHS0hJAwQEmNDQUhw8fBlCwTKJatWomEa8xmdLr9rzq1KmDhw8fIjc3V2q7f/8+8vLyjPYYfvvtNwwbNgy//PILhg4dqnJdnTp1kJiYqPJl8smTJ8jMzDSJ5zw8PByZmZk4d+4cgIJ4nz17plIeOjU1FcnJyQaPt3HjxvD29paOS/v27QNQcEqQkydPSvFmZGTg6dOn0risrCwkJiaaxPNbOKNSeGyqU6cOlEol7t27J/XJy8vDgwcPDB5v4femfv36Sd+n/P390bJlS+nYpO67VWH8xnx+V6xYAU9PT7z22msq7Zq+C965c8fg8Z45cwYZGRkqqyrq1q2L5s2bS8cmvcWrp31OevfLL78IW1tblWobPXv2NLnNuqVVk6tbt6744IMPpMv//POP0apxXbt2Tfj4+Ijw8HCNBSiaNm2qsmH6+vXrRtlEmZWVVWJDf1JSknB0dBTz58+X2tq0aaOyOT4+Pl5YWlpKG+uMSV01uU6dOqlssn7w4IGwsrIyeDWu/Px8kZCQoNJ25swZIZPJVErE9ujRQ/Ts2VO6nJiYKOzs7MTSpUsNFmuh4cOHl/i3/+abb4rWrVtLl8PCwkTnzp2ly8+ePRMODg5i4cKFBouzuCFDhogWLVqovW7QoEGibdu2UoWh1NRU4erqKubOnWvIEIUQQvj6+ooJEyaotL3xxhuiZcuW0uXhw4eLpk2bShV7MjMzhY+Pj/jss88MGquxmdLrVhpNBRRu3rwpLCwsVD7XFyxYIKpWrSpVujSkTZs2CWtra7Fy5Uq119+9e1dlw7QQQixevFjY2dmVOE6Ut4SEhBLHz8IqpoVVVxMSEoS1tbXK5/qyZcuEtbV1ic9dQ1NXQOHp06fC1tZW5XN93bp1Qi6Xi3v37hk0vidPnoicnByVtvnz5wsrKyvx6NEjIUTBBnoHBweVzfHbt28XMplMXLt2zaDxKpVKUatWLZV/+zk5OcLPz09MnTpVCFFQIdHV1VWlcNbevXtLVBc1pLy8PFGtWjUxadKkEtfl5uYKLy8vlc/1I0eOCADizJkzhgxTKj6xf/9+qa2wWuO0adOEEAWvgZ+fn8qpQaKiogQAncrvm2wypFAoRKtWrUSzZs3E4sWLRXh4uLC3tzf4i6HJxYsXxYYNG8T06dMFAPHDDz+IDRs2qFSQ2bVrl5DL5eLjjz8W3333nahVq5YICwszeKxJSUmievXqolatWmL9+vViw4YN0l9hmVwhhDh48KCwsrISo0ePFt9//70ICAgQr7zyisHLrT558kQ0aNBATJ48WaxYsULMmzdPBAYGikaNGqmUzjx69KiwtrYWI0aMEIsWLRJBQUHi5ZdfVikRaSzqkqFTp04JW1tb8c4774jFixeLxo0bizZt2hi8PGxOTo6oX7++GDdunFi2bJmYMWOGcHFxEb1791Z57s6fPy+qVKkihg4dKn744QfRvHlz0aJFC6n0uiHdvXtXeHt7i/DwcLF8+XIxZswYYW1trXIujEuXLgkHBwcxePBg8cMPP4hWrVqJRo0aSaXXDa3wS4amMrFXr14Vzs7Oon///mLJkiWibdu2on79+irl7g1l6dKlwsrKSnz88cdixYoVYuTIkcLS0lLlXGTx8fHCw8ND9O7dWyxZskR06tRJpZS5uTCl102d/fv3iw0bNoi+ffsKHx8f6bO+sNqZEEKMHz9euLq6innz5omZM2cKGxsbsWjRIoPHevDgQSGXy8Vrr72mclzasGFDifLVTk5OIiIiQsyePVvY2dmJ//3vfwaPd9euXaJ169biiy++ECtXrhSTJk0STk5OYsSIESr9Zs+eLapWrSq+/PJL8eWXX4qqVauKOXPmGDze4tQlQ0IUJByFJYwjIiKEo6Oj9GXTkE6ePCmaNGkiZs2aJZYvXy5GjhwprK2tVX4EFUKIJUuWCBsbGzFjxgwxf/584ebmplLK3JB2794tHBwcxLRp08SyZctEt27dhI+Pj5S8CVFwrhwrKysxbdo0sWDBAuHp6alSytzQdu7cKQBoTB7Xr18v5HK5mDx5svjmm2+Ej4+PePvttw0b5P/Xt29f4e3tLebNmyd+/vln0alTJ+Hh4aGSqG/ZskXI5XIxceJE8e233wo/Pz8xaNAgne5HJoQQzzmDVe4yMzOxdOlSXLhwAZ6ennjvvfcQGBho7LAAAJs2bcLvv/9eov2DDz5Ahw4dpMvR0dFYt24dUlNT0bZtW7z77rtSBSxDuXfvHiZOnKj2usmTJ6uc2fmff/7B6tWr8ezZM7Rq1QojRowosVHREFJSUrBq1SrExMTA0dERrVq1wuuvv17iubt06RJWrFiBpKQktGzZEu+//760htSYJkyYgODgYAwZMkSl/cqVK1i+fDkSExPRrFkzjBo1SmXvjqFkZGRg9erVOH/+PFxcXNCxY0eEhoaW6Hft2jUsW7YMCQkJaNy4MUaPHm20M9UnJCRg6dKluH37NqpXr47w8PASnwc3b97EL7/8gocPH6Jhw4YYPXq02rPDG8KxY8fwww8/4KefftJ49vn4+Hj89NNPuH//PurXr48xY8bAycnJwJEWOHXqFLZt24aEhAT4+vpi8ODBJaqm3b9/H0uXLsWdO3dQr149jB49WipsY05M6XUrbvr06bh+/XqJ9p9++kklxs2bN+PPP/+EXC5H//790a1bN0OGCQDYvXs31q1bp/a6X3/9VSqYAxRUGIuMjISFhQX69OmjsrfEkK5evYoNGzYgPj4ePj4+6NWrF9q3b1+i386dO7F7924AQO/evUssRzKGxMREjB07FjNmzEBQUJDKdZGRkdixYweUSiV69eqFvn37GiXG+Ph4rF27Vvqc79+/P5o0aVKi34EDB7B161bk5OSge/fuGDRokBGiLXD+/HmsWbMGqampqF+/PkaMGAFnZ2eVPkeOHMGmTZugUCjQtWtXvPHGGypbFQyp8Dg6b948jX2OHTuGjRs3IjMzE506dUJ4eLjKv0dDUSqV2LRpE06cOIHMzEwEBARg+PDhJY47UVFRWL9+PdLT09GxY0e8/fbb0rJvbZh0MkRERERERFReTLaAAhERERERUXliMkRERERERGaJyRAREREREZklJkNERERERGSWmAwREREREZFZYjJERERERERmickQERERERGZJSZDRERERERklpgMERERERGRWWIyREREREREZonJEBERERERmSUmQ0REREREZJb+H2Yq2mTJkTBzAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "\n", "fig, axs = plt.subplots(ncols=2, nrows=1, figsize=(10, 4))\n", "kaplan_meier.plot(\"sf\", ax=axs[0])\n", "nelson_aalen.plot(\"chf\", ax=axs[1])\n", "axs[0].set_title(\"Kaplan-Meier\")\n", "axs[1].set_title(\"Nelson-Aalen\")\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.0rc1" } }, "nbformat": 4, "nbformat_minor": 5 }