{ "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 relife.data import load_circuit_breaker\n", "\n", "data = load_circuit_breaker()\n", "print(data[\"time\"])\n", "print(data[\"event\"])\n", "print(data[\"entry\"])" ] }, { "cell_type": "code", "execution_count": 2, "id": "5cd6ca17-3980-4e90-bf2c-1b8fa35843a1", "metadata": {}, "outputs": [], "source": [ "from relife.lifetime_model import KaplanMeier, NelsonAalen" ] }, { "cell_type": "code", "execution_count": 3, "id": "713b3e2f-7265-43a0-b3f9-a8e28dc26f35", "metadata": {}, "outputs": [], "source": [ "kaplan_meier = KaplanMeier()\n", "kaplan_meier.fit(\n", " data[\"time\"],\n", " event = data[\"event\"],\n", " entry = data[\"entry\"],\n", ")" ] }, { "cell_type": "code", "execution_count": 4, "id": "98e52f10-2452-4c4f-8bb4-1c0968a89422", "metadata": {}, "outputs": [], "source": [ "nelson_aalen = NelsonAalen()\n", "nelson_aalen.fit(\n", " data[\"time\"],\n", " event = data[\"event\"],\n", " entry = data[\"entry\"],\n", ")" ] }, { "cell_type": "code", "execution_count": 5, "id": "87dfc944-1c3e-4dea-b489-0c3672c519bf", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA04AAAFzCAYAAAAJ21nbAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAAcpdJREFUeJzt3XlcVGX7P/DPsIvsKouGLIL7hhuC5RaJS5Raaq6A6fMlUUTUlCcVcQHN3TQpF9TSXNI0tTBD0Ucjd0xDMRXDx1hcAgRiG87vD3/Mw8g2BwYODJ/36zWvZu4558w1J5x7rrnvc90yQRAEEBERERERUbm0pA6AiIiIiIiormPiREREREREVAkmTkRERERERJVg4kRERERERFQJJk5ERERERESVYOJERERERERUCSZORERERERElWDiREREREREVAkdqQOobUVFRfjrr79gbGwMmUwmdThERA2KIAh48eIFmjdvDi0t/nZXjH0TEZE0xPRLDS5x+uuvv2Brayt1GEREDdqjR4/w2muvSR1GncG+iYhIWqr0Sw0ucTI2Ngbw8uSYmJhIHA0RUcOSmZkJW1tbxWcxvcS+iYhIGmL6pQaXOBVPgTAxMWHnREQkEU5HU8a+iYhIWqr0S5xgTkREREREVAkmTkRERERERJVg4kRERERERFSJBneNExERUX0kCAIKCwshl8ulDoUaGG1tbejo6PDaRGrwmDgRERHVcfn5+UhOTkZOTo7UoVADZWhoCBsbG+jp6UkdCpFkmDgRERHVYUVFRUhMTIS2tjaaN28OPT09/vJPtUYQBOTn5+PJkydITEyEs7MzF6+mBouJExERUR2Wn5+PoqIi2NrawtDQUOpwqAFq1KgRdHV18eeffyI/Px8GBgZSh0QkCUl/Mjh37hy8vLzQvHlzyGQyHDlypNJ9YmJi0K1bN+jr68PJyQk7d+6s8TiJiIikxl/5SUr8+yOSOHHKzs5Gly5dsHnzZpW2T0xMxLBhwzBgwADExcUhMDAQU6ZMwcmTJ2s4UiIiagjCw8PRs2dPGBsbw9LSEsOHD0dCQkKF++zcuRMymUzpxl/kiYg0j6RT9YYMGYIhQ4aovH1ERAQcHBywZs0aAEC7du1w/vx5rFu3Dp6enqJeOye/EDr5haL2qYsa6WpzrjsRkZqcPXsW/v7+6NmzJwoLC/Hvf/8bgwYNQnx8PBo3blzufiYmJkoJFj+XiYg0T726xik2NhYeHh5KbZ6enggMDCx3n7y8POTl5SkeZ2ZmAgB6LY+Gln79nyvew84cB/3c2EkTEalBVFSU0uOdO3fC0tISV69eRd++fcvdTyaTwdrauqbDa7Ds7e0RGBhYYX9fFy1evBhHjhxBXFycZDHs3LkTgYGBSE9Pr3C7xYsXY8uWLUhLS8N3332H4cOH10p8RPVJvUqcUlJSYGVlpdRmZWWFzMxM/PPPP2jUqFGpfcLDwxEaGlpbIda6K3/+jZO/p0Jfp2qzLvV1tCpNutxaNanSsYmI6ruMjAwAgIWFRYXbZWVlwc7ODkVFRejWrRvCwsLQoUOHcrcv70c9TePj44Ndu3YhPDwc8+fPV7QfOXIEI0aMgCAIEkZXuf/+979wdHRE69atcevWLanDqTG3b99GaGgovvvuO/Tu3Rvm5uZSh0QaLPb+M7Ufs7a+q9arxKkqgoODERQUpHicmZkJW1tbXPrkTZiYmEgYWfXk5MvRY9nPAAC/r69W+Th2TQyx2Eu5c381mRL7B85Ei4g0QVFREQIDA9GnTx907Nix3O3atGmDHTt2oHPnzsjIyMDq1avh7u6O33//Ha+99lqZ+2j6j3olGRgYYOXKlfi///u/eveFfOfOnRg9ejTOnTuHixcvwtXVVeqQasT9+/cBAO+++y5nsBBVoF6VSLG2tkZqaqpSW2pqKkxMTMocbQIAfX19mJiYKN0AwFBPp17fmjTWQw+76ndAfz7Lge/Oy0q34O9uIrdArriJ/UUw9v4z0TciorrG398ft27dwr59+yrczs3NDZMmTULXrl3Rr18/HD58GM2aNcMXX3xR7j7BwcHIyMhQ3B49eiQqNkEQkJNfWOu3qowQeXh4wNraGuHh4eVuc/78ebzxxhto1KgRbG1tERAQgOzs7HLf++LFi9GyZUvo6+ujefPmCAgIUDz/999/Y9KkSTA3N4ehoSGGDBmCP/74Q/H8zp07YWZmhpMnT6Jdu3YwMjLC4MGDkZycXOp1IiMjMXHiRIwbNw7bt28vFcu8efPQunVrGBoawtHREQsXLkRBQUGF52Pbtm1o164dDAwM0LZtW3z++eeK5x4+fAiZTIbDhw9jwIABMDQ0RJcuXRAbG1vhMdPT0/F///d/sLKygoGBATp27Ijjx48rbVPe+128eDG8vLwAvKycx8SJqHz1asTJzc0NP/zwg1LbqVOn4ObmJlFE0pHJZDjo54aYhCdVPsbiY7/jz2elV6EvTqaKlTUqVRFVpv+9qrLkiaNYRFSbpk+fjuPHj+PcuXPljhqVR1dXFy4uLrh371652+jr60NfX7/K8f1TIEf7RbVfUTZ+iScM9cR9ddDW1kZYWBjGjRuHgICAUufz/v37GDx4MJYtW4YdO3bgyZMnmD59OqZPn47IyMhSxzt06BDWrVuHffv2oUOHDkhJScGNGzcUz/v4+OCPP/7A999/DxMTE8ybNw9Dhw5FfHw8dHV1AQA5OTlYvXo1vvrqK2hpaWHChAmYM2cO9uzZozjOmTNnkJOTAw8PD7Ro0QLu7u5Yt26dUpEQY2Nj7Ny5E82bN8fNmzcxdepUGBsb4+OPPy7zXOzZsweLFi3Cpk2b4OLiguvXr2Pq1Klo3LgxvL29Fdt98sknWL16NZydnfHJJ59g7NixuHfvHnR0Sp/7oqIiDBkyBC9evMDXX3+NVq1aIT4+Htra2optKnq/c+bMgb29PXx9fUslj0SkTNLEKSsrS6ljSUxMRFxcHCwsLNCyZUsEBwfj8ePH2L17NwDAz88PmzZtwscff4zJkyfj9OnTOHDgAE6cOCHVW5CUTCaDga525RuWI3xEJ+QVFim1lZVMvZpIVUZsogVUnmypMirF5IqIqksQBMyYMQPfffcdYmJi4ODgIPoYcrkcN2/exNChQ2sgwvppxIgR6Nq1K0JCQkqN3ISHh2P8+PGKwg/Ozs7YuHEj+vXrhy1btpQq7Z6UlARra2t4eHhAV1cXLVu2RK9evQBAkTBduHAB7u7uAF4mK7a2tjhy5AhGjRoFACgoKEBERARatWoF4GWivGTJEqXX2b59Oz744ANoa2ujY8eOcHR0xMGDB+Hj46PYZsGCBYr79vb2mDNnDvbt21du4hQSEoI1a9Zg5MiRAAAHBwfEx8fjiy++UEqc5syZg2HDhgEAQkND0aFDB9y7dw9t27Ytdcyff/4Zly5dwu3bt9G6dWsAgKOjo9I2Fb1fIyMjmJmZAQALnBBVQtLE6cqVKxgwYIDicfG1SN7e3ti5cyeSk5ORlJSkeN7BwQEnTpzArFmzsGHDBrz22mvYtm2b6FLkmkRsslAyASkr8Xo1mSpvVKoiYhMtoPJkS5VRLI5aEVF1+fv7Y+/evTh69CiMjY2RkpICADA1NVVMCZ80aRJatGihmHq2ZMkS9O7dG05OTkhPT8eqVavw559/YsqUKTUWZyNdbcQvqf2+r1E1fqxbuXIlBg4ciDlz5ii137hxA7/99pvSaI8gCCgqKkJiYiLatWuntP2oUaOwfv16ODo6YvDgwRg6dCi8vLygo6OD27dvQ0dHR+lapCZNmqBNmza4ffu2os3Q0FCRRACAjY0N0tLSFI/T09Nx+PBhnD9/XtE2YcIEbN++XSlx2r9/PzZu3Ij79+8jKysLhYWF5V4/nZ2djfv37+PDDz/E1KlTFe2FhYUwNTVV2rZz585KsQFAWlpamYlTXFwcXnvtNUXSVJbK3i8RqUbSxKl///4VzpfeuXNnmftcv369BqPSbKokDyUTkLJGpSpSlUQLqDzZUmUUq7qjVkysiGjLli0AXvY1JUVGRiq+MCclJUFL63+XCP/999+YOnUqUlJSYG5uju7du+OXX35B+/btayxOmUwmesqc1Pr27QtPT08EBwcrJR9ZWVn4v//7P6XrlIq1bNmyVJutrS0SEhLw888/49SpU5g2bRpWrVqFs2fPqhxL8ZS9YjKZTOn7yN69e5Gbm6uUgBUnc3fv3kXr1q0RGxuL8ePHIzQ0FJ6enjA1NcW+ffsUa02+KisrCwCwdevWUkUmSk6rezW+4n6tqKjsvri8a7zLO17xMet6RUOiuqh+fepSrXg1gRBTvEFsogWolmypMopV3VGrst4nkymihkWVL5MxMTFKj9etW4d169bVUESaZcWKFejatSvatGmjaOvWrRvi4+Ph5OSk8nEaNWoELy8veHl5wd/fH23btsXNmzfRrl07FBYW4uLFi4qpes+ePUNCQoKoRHb79u2YPXu2UoIHANOmTcOOHTuwYsUK/PLLL7Czs8Mnn3yieP7PP/8s95hWVlZo3rw5Hjx4gPHjx6scS2U6d+6M//73v4qEjohqDhMnqlRlyUNl0/8qU1mypeooVnVHrcpKrF5NpphIERFVXadOnTB+/Hhs3LhR0TZv3jz07t0b06dPx5QpU9C4cWPEx8fj1KlT2LRpU6lj7Ny5E3K5HK6urjA0NMTXX3+NRo0awc7ODk2aNMG7776LqVOn4osvvoCxsTHmz5+PFi1a4N1331Upxri4OFy7dg179uwpNTVu7NixWLJkCZYtWwZnZ2ckJSVh37596NmzJ06cOIHvvvuuwmOHhoYiICAApqamGDx4MPLy8nDlyhX8/fffSkunVOTx48d48803sXv3bvTq1Qv9+vVD37598d5772Ht2rVwcnLCnTt3IJPJMHjwYJWOSUSqYeJE1SZ2+t+rKku2VBnFUseoVRsrY4R4tWcpViKiGrRkyRLs379f8bhz5844e/YsPvnkE7zxxhsQBAGtWrXCmDFjytzfzMwMK1asQFBQEORyOTp16oRjx46hSZOXfVFkZCRmzpyJt99+G/n5+ejbty9++OGHUtPVyrN9+3a0b9++zOuJRowYgenTp+OHH37AO++8g1mzZmH69OnIy8vDsGHDsHDhQixevLjcY0+ZMgWGhoZYtWoV5s6di8aNG6NTp06KwhiqKCgoQEJCAnJy/tfnHTp0CHPmzMHYsWORnZ0NJycnrFixQuVjEpFqZEIDm+SamZkJU1NTZGRk1OsFcDVNdddyEgRBLaNWERO6Q1/nf9cuiC2tzhEpoorV1Gdwamoq5syZg+joaKSlpZWacieXy9X2WjWhovOSm5uLxMREODg4lKowR1Rb+HdI6lIT63dW5/uXmH6JI05UJ4iZDliW6oxa5RUWwe/rqwCg+G+xV6f3ib1OiokUUe3w8fFBUlISFi5cCBsbG44cExGR2jFxonqhJqcD6utooY2VMRJSX5R6rrLFgKuy2C8Rqd/58+fxn//8B127dpU6FCIi0lBMnEhjVLUaoEwmQ4hX+yotBlxZIsWFe4lqh62tLcsrExFRjWLiRBqrOosDA6otBswRKaK6Yf369Zg/fz6++OIL2NvbSx0OERFpICZOROV4dXqfOhIpoPJRKY5AEYk3ZswY5OTkoFWrVjA0NCxVQe358+cSRUZERJqCiRPR/1dWwlLRGlVVSaQA8dP7mEgRVW79+vVSh0BERBqOiRNRBSq6bqoqiRRQOpni+lFE1eft7S11CEREpOGYOBGpSWWJFFB2MpWQ+gKZuYWK9aM4AkVUNXK5HEeOHMHt27cBAB06dMA777wDbe3ylyogIiJSFRMnIhHErDdVVgn0kslUeetHVaVSH5Mpauju3buHoUOH4vHjx2jTpg0AIDw8HLa2tjhx4gRatWolcYRERFTfMXEiUqPKSqKXTKbKWz+KU/mIxAsICECrVq3w66+/wsLCAgDw7NkzTJgwAQEBAThx4oTEEdYMVZddUIeq/EDj4+OD9PR0HDlyBACQkpKC5cuX48SJE3j8+DEsLS3RtWtXBAYG4s033yz3ONeuXcO8efNw+fJlaGtr47333sPatWthZGSk2Kasz8hvvvkGH3zwAQDg+vXrmDx5Mv744w8MGDAAu3btUvytFBYWwtXVFVu2bEGvXr1Ev08iahi0pA6AqKEqXj8q0qen4mbXxLDUdsVT+XIL5MgtkJe5Vk3s/WdKN6KG5uzZs/j0008VX4QBoEmTJlixYgXOnj0rYWRU7OHDh+jevTtOnz6NVatW4ebNm4iKisKAAQPg7+9f7n5//fUXPDw84OTkhIsXLyIqKgq///47fHx8Sm0bGRmJ5ORkxW348OGK56ZMmYKBAwfi2rVryMjIQFhYmOK5NWvWoE+fPkyaiKhCHHEiqkHVqdSn6lQ+gOtFEenr6+PFixel2rOysqCnpydBRPSqadOmQSaT4dKlS2jcuLGivUOHDpg8eXK5+x0/fhy6urrYvHkztLRe/t4bERGBzp074969e3ByclJsa2ZmBmtr6zKPc/v2bezZswetW7fG2LFjcfz4cQDAgwcPsH37dly9erXM/YiIinHEiaiWubVqonQrqTiRMtDVhomBDtpYGZfav3gqX8lb6LH4MkeiiBqKt99+G//6179w8eJFCIIAQRDw66+/ws/PD++8847U4TV4z58/R1RUFPz9/ZWSpmJmZmbl7puXlwc9PT1F0gQAjRo1AgCcP39eaVt/f380bdoUvXr1wo4dO5Q+F7t06YJTp06hsLAQ0dHR6Ny5MwDAz88Pn376KYyNS3/eEhGVxMSJSGLlJVKqTuUDXk7nK1nBj1P3qKHZuHEjWrVqBTc3NxgYGMDAwAB9+vSBk5MTNmzYIHV4Dd69e/cgCALatm0ret+BAwciJSUFq1atQn5+Pv7++2/Mnz8fAJCcnKzYbsmSJThw4ABOnTqF9957D9OmTcNnn32meH7btm349ttv0apVK+jp6SE4OBhfffUVDA0N0bNnT3h6esLJyQkLFiyo/hsmIo3EqXpEdUxFBSZeLXFecjpfyXZO3aOGxszMDEePHsUff/yBO3fuAADatWunNI2LpKPqiLifnx++/vprxeOsrCx06NABu3btQlBQEIKDg6GtrY2AgABYWVkpjUItXLhQcd/FxQXZ2dlYtWoVAgICALycEljyerdnz54hJCQE586dw4wZM+Du7o7Dhw+jZ8+ecHV1hZeXV3XfNhFpGCZORPVIWSXOi5W8DoqV+KihcnZ2hrOzs9Rh0CucnZ0hk8kUSW15lixZgjlz5pRqHzduHMaNG4fU1FQ0btwYMpkMa9euhaOjY7nHcnV1xdKlS5GXlwd9ff1SzwcFBSEwMBCvvfYaYmJisGzZMjRu3BjDhg1DTEwMEyciKoWJE1EdV9EIVHklzYun7hUnWVxElzRRUFAQli5disaNGyMoKKjCbdeuXVtLUVFZLCws4Onpic2bNyMgIKDUdU7p6ekwMzODpaUlLC0tyz2OlZUVAGDHjh0wMDDAW2+9Ve62cXFxMDc3LzNpio6Oxu3btxEZGQng5eLJBQUFAKD4LxHRq5g4EdVjxddBlVWJr6Kpe0ykSBNcv35d8SX3+vXrEkdDldm8ebOi5PeSJUvQuXNnFBYW4tSpU9iyZQtu375d7r6bNm2Cu7s7jIyMcOrUKcydOxcrVqxQFJU4duwYUlNT0bt3bxgYGODUqVMICwsrc/QqNzcX06dPxzfffKOY6tenTx9s3rwZ/v7+OHToEBNtIioTEyeieqayEufFOHWPNN2ZM2fKvE91k6OjI65du4bly5dj9uzZSE5ORrNmzdC9e3ds2bKlwn0vXbqEkJAQZGVloW3btvjiiy8wceJExfPF5cpnzZoFQRDg5OSEtWvXYurUqaWOFRoaimHDhqFr166Kto0bN2LcuHHo27cvxo8fj/fee09t75uINIdMaGA1jDMzM2FqaoqMjAyYmJhIHQ6RWhQnToIgIPRYfKmpewAQMaE79HVe/rr66ggUR5yottTUZ/DkyZOxYcOGUiWls7OzMWPGDOzYsUNtr1UTKjovubm5SExMhIODAwwMDCSKkBo6/h2SutREtd/qfI8R0y8xcSLSACU/hARBKHPqXkmqjEAxmaKaUFOfwdra2khOTi51fczTp09hbW2NwsJCtb1WTWDiRHUd/w5JXepz4sSpekQaoLwCEhUVj8jMLSx3BIqovsjMzFQsePvixQulL3RyuRw//PBDhcUGiIiIVMXEiUiDVVQ8gtdAkSYwMzODTCaDTCZD69atSz0vk8kQGhoqQWRERKRpmDgRabiSaz+pWr4cYOU9qh/OnDkDQRAwcOBAHDp0CBYWForn9PT0YGdnh+bNm0sYIRERaQomTkQaqLype6qWLwc4fY/qh379+gEAEhMT0bJlS/7NEhFRjWHiRNTAlByBKunVIhKvTt/jCBTVZadPn4aRkRFGjRql1H7w4EHk5OTA29tbosjUp4HVcqI6hn9/RICW1AEQkXSKp+6VpXj6HlF9EB4ejqZNm5Zqt7S0RFhYmAQRqY+uri4AICcnR+JIqCEr/vsr/nskaog44kTUAKg6dQ8ov4Q5UV2WlJQEBweHUu12dnZISkqSICL10dbWhpmZGdLS0gAAhoaGnJJItUYQBOTk5CAtLQ1mZmbQ1i49Y4GooWDiRNTAlTd1D1C+7unVa544dY/qEktLS/z222+wt7dXar9x4waaNKn/f5vW1tYAoEieiGqbmZmZ4u+QqKFi4kTUAJU3AvUqMSXLmUiRlMaOHYuAgAAYGxujb9++AICzZ89i5syZ+OCDDySOrvpkMhlsbGxgaWmJgoICqcOhBkZXV5cjTURg4kREr+CiuVQfLV26FA8fPsSbb74JHZ2XXVtRUREmTZpU769xKklbW5tfYImIJCITGliZlMzMTJiamiIjIwMmJiZSh0NUJ7w6WiQIQpkly0uqyqK5HIWimv4Mvnv3Lm7cuIFGjRqhU6dOsLOzU/tr1AT2TUTUUJQ3y6U6qvP9QsznL0eciKjMD5ziDzYxi+YSSa1169Zo3bq11GEQEZEGYuJERBVSddFcTt0jKcnlcuzcuRPR0dFIS0tDUZFyKf3Tp09LFBkREWkKJk5EVClVFs1VZepeZcPznMpHVTVz5kzs3LkTw4YNQ8eOHZnEExGR2lU5ccrPzy/zV72WLVtWOygiqrtYPILqon379uHAgQMYOnSo1KEQEZGGEp04/fHHH5g8eTJ++eUXpXZBECCTySCXy9UWHBFJR9VFc0tO3RM7AkWkLnp6enBycpI6DCIi0mCiEycfHx/o6Ojg+PHjsLGx4Zciogao5NQ9VUegirflZwbVhNmzZ2PDhg3YtGkT/8aIiKhGiC5H3rhxY1y9ehVt27atqZhqFEu+EqlHyeuVVClfDogfheI1T5qnpj6DR4wYgTNnzsDCwgIdOnSArq6u0vOHDx9W22vVBPZNRNRQNKhy5O3bt8fTp0+rHBwRaR5VRqAAljCnmmNmZoYRI0ZIHQYREWkw0YnTypUr8fHHHyMsLAydOnUq9asefykjathevQYKqHgUikgdIiMjpQ6BiIg0nOjEycPDAwDw5ptvKrWzOAQRFSuvfDkRERFRfSU6cTpz5kxNxEFE9Ux5VfcqU3Ik6lWvFo949Zi85onK4+DgUOG1cw8ePKjFaIiISBOJTpz69etXE3EQUQNR0ZQ9ljCnqgoMDFR6XFBQgOvXryMqKgpz586VJigiItIoVVoANz09Hdu3b8ft27cBAB06dMDkyZNhamoq+libN2/GqlWrkJKSgi5duuCzzz5Dr169yt1+/fr12LJlC5KSktC0aVO8//77CA8Ph4GBQVXeChHVgooKRpTE4hFUVTNnziyzffPmzbhy5UotR0NERJpIq/JNlF25cgWtWrXCunXr8Pz5czx//hxr165Fq1atcO3aNVHH2r9/P4KCghASEoJr166hS5cu8PT0RFpaWpnb7927F/Pnz0dISAhu376N7du3Y//+/fj3v/8t9m0QkZq5tWqidCupuGBEpE/PMm8RE7pLFDVpuiFDhuDQoUMqbx8eHo6ePXvC2NgYlpaWGD58OBISEird7+DBg2jbti0MDAzQqVMn/PDDD9UJm4iI6iDRidOsWbPwzjvv4OHDhzh8+DAOHz6MxMREvP3226WmSlRm7dq1mDp1Knx9fdG+fXtERETA0NAQO3bsKHP7X375BX369MG4ceNgb2+PQYMGYezYsbh06ZLYt0FEtay4YERZt5KL5FYk9v6zUjeiinz77bewsLBQefuzZ8/C398fv/76K06dOoWCggIMGjQI2dnZ5e7zyy+/YOzYsfjwww9x/fp1DB8+HMOHD8etW7fU8RaIiKiOED1V78qVK9i6dSt0dP63q46ODj7++GP06NFD5ePk5+fj6tWrCA4OVrRpaWnBw8MDsbGxZe7j7u6Or7/+GpcuXUKvXr3w4MED/PDDD5g4cWK5r5OXl4e8vDzF48zMTJVjJKKqK6uQgyqJTsniEa8WiyAqj4uLi9LfiiAISElJwZMnT/D555+rfJyoqCilxzt37oSlpSWuXr2Kvn37lrnPhg0bMHjwYMW1VEuXLsWpU6ewadMmREREVOHdEBFRXSQ6cTIxMUFSUhLatm2r1P7o0SMYGxurfJynT59CLpfDyspKqd3Kygp37twpc59x48bh6dOneP311yEIAgoLC+Hn51fhVL3w8HCEhoaqHBcRSatk8QgWiyBVDR8+XOmxlpYWmjVrhv79+5fqr8TIyMgAgApHrWJjYxEUFKTU5unpiSNHjpS7D3/UIyKqf0QnTmPGjMGHH36I1atXw93dHQBw4cIFzJ07F2PHjlV7gCXFxMQgLCwMn3/+OVxdXXHv3j3MnDkTS5cuxcKFC8vcJzg4WKlDy8zMhK2tbY3GSUTilFc8or4Wi6jKFEKWWhcvKCgIS5cuRePGjTFgwAC4ubmVWpS9OoqKihAYGIg+ffqgY8eO5W6XkpJS5o+AKSkp5e7DH/WIiOof0YnT6tWrIZPJMGnSJBQWFgIAdHV18dFHH2HFihUqH6dp06bQ1tZGamqqUntqaiqsra3L3GfhwoWYOHEipkyZAgDo1KkTsrOz8a9//QuffPIJtLRKXyehr68PfX19leMiotpXXDyieJpeXmFRhWXLS+JaTw3XZ599hnnz5ikSp+TkZFhaWqrt+P7+/rh16xbOnz+vtmMW4496RET1j+jESU9PDxs2bEB4eDju378PAGjVqhUMDQ1FH6d79+6Ijo5WTLEoKipCdHQ0pk+fXuY+OTk5pZIjbe2Xv0QLgiDynRBRXVJcPIJIVfb29ti4cSMGDRoEQRAQGxsLc3PzMrct7/qk8kyfPh3Hjx/HuXPn8Nprr1W4rbW1tagfAQH+qEdEVB9VaR0nADA0NESnTp2q9eJBQUHw9vZGjx490KtXL6xfvx7Z2dnw9fUFAEyaNAktWrRAeHg4AMDLywtr166Fi4uLYqrewoUL4eXlpUigiKhhq26lPY5Y1R+rVq2Cn58fwsPDIZPJMGLEiDK3k8lkkMvlKh1TEATMmDED3333HWJiYuDg4FDpPm5uboiOjlaqLHvq1Cm4ubmp9JpERFQ/qJQ4jRw5Ejt37oSJiQlGjhxZ4baHDx9W+cXHjBmDJ0+eYNGiRUhJSUHXrl0RFRWlmCuelJSkNMK0YMECyGQyLFiwAI8fP0azZs3g5eWF5cuXq/yaRCSdypKSspKe2q6yxxLn9Udx2e+srCyYmJggISGh2lP1/P39sXfvXhw9ehTGxsaK65RMTU3RqFEjAKV/1Js5cyb69euHNWvWYNiwYdi3bx+uXLmCL7/8snpvkIiI6hSVEidTU1PFlxUTExO1fnGZPn16uVPzYmJilB7r6OggJCQEISEhant9IqrbWGWPKmNkZIQzZ87AwcFBaamMqtiyZQsAoH///krtkZGR8PHxAVD6Rz13d3fs3bsXCxYswL///W84OzvjyJEjFRaUICKi+kcmNLCLgzIzM2FqaoqMjAyYmJhIHQ4RlVA82iMIAkKPxZeqsgcAkT49NfJaqIYyRZCfwWXjeSGihqImZnZUpw8V8/kr+qe5gQMH4vDhwzAzMyv1osOHD8fp06fFHpKISEl1quwRERFR7Who09tFJ04xMTHIz88v1Z6bm4v//Oc/agmKiIhV9oiIiKguUTlx+u233xT34+PjlRb2k8vliIqKQosWLdQbHRHRK0oWiwBqp2AEERERkcqJU9euXSGTySCTyTBw4MBSzzdq1AifffaZWoMjInrVq1P2WDCCiIiIaoPKiVNiYiIEQYCjoyMuXbqEZs2aKZ7T09ODpaUl11Iiomp59eLO4rnT+jpaaGNlXGaxiITUF8grLOK0vgaosuUxShKzVAYREVFZVE6c7OzsAABFRUWVbElEpF6vFosAWDCCXi6VUUwQBHz33XcwNTVFjx49AABXr15Fenq6qASLiIioPKKLQ4SHh8PKygqTJ09Wat+xYweePHmCefPmqS04IqJiLBZBr4qMjFTcnzdvHkaPHo2IiAjF7Ae5XI5p06axvDcREamFVuWbKPviiy/Qtm3bUu0dOnRARESEWoIiIhIjr7AIuQVy5BbI0cCWpqP/b8eOHZgzZ47SlHFtbW0EBQVhx44dEkZGRESaQvSIU0pKCmxsbEq1N2vWDMnJyWoJiohIjJJT9lgsomEqLCzEnTt30KZNG6X2O3fucIo5ERGphejEydbWFhcuXICDg4NS+4ULF9C8eXO1BUZEVJHyCkawWETD5Ovriw8//BD3799Hr169AAAXL17EihUr4OvrK3F0RESkCUQnTlOnTkVgYCAKCgoUZcmjo6Px8ccfY/bs2WoPkIioLK8WjChZLOLVtZ7E4LpQ9dPq1athbW2NNWvWKGY/2NjYYO7cueybiIhILUQnTnPnzsWzZ88wbdo05OfnAwAMDAwwb948BAcHqz1AImq4yitPXqy8ghHVqbbHqX71T2FhIfbu3Qtvb298/PHHyMzMBAAWhSAiIrUSnTjJZDKsXLkSCxcuxO3bt9GoUSM4OztDX1+/JuIjIlJJRWs9iSHVVL9Xk8JXk0Yqn46ODvz8/HD79m0ATJiIiKhmiE6cihkZGaFnz57qjIWIqMrKWutJDK4LVb/16tUL169fV6w5SEREpG6iE6fs7GysWLEC0dHRSEtLK1Wt6MGDB2oLjohIDK711HBNmzYNs2fPxn//+190794djRs3Vnq+c+fOEkVGRESaQnTiNGXKFJw9exYTJ06EjY0NrwMgIiLJffDBBwCAgIAARZtMJoMgCJDJZJDL5VKFRkREGkJ04vTjjz/ixIkT6NOnT03EQ0REJFpiYqLUIRARkYYTnTiZm5vDwsKiJmIhIqpQWQUTXi2qQA0Tr20iIqKaJjpxWrp0KRYtWoRdu3bB0NCwJmIiIiKqkvj4eCQlJSmWyyj2zjvvSBQRERFpCtGJ05o1a3D//n1YWVnB3t4eurq6Ss9fu3ZNbcERERGp4sGDBxgxYgRu3rypuLYJgOI6XF7jRERE1SU6cRo+fHgNhEFERFR1M2fOhIODA6Kjo+Hg4IBLly7h2bNnmD17NlavXi11eEREpAFEJ04hISE1EQcREVGVxcbG4vTp02jatCm0tLSgpaWF119/HeHh4QgICMD169elDpGIiOq5Ki+AS0Skqaq6iG4xfR0tLtVQy+RyOYyNjQEATZs2xV9//YU2bdrAzs4OCQkJEkdHRESaQHTipKVV8RcCziMnovrO7+ur1dq/jZUxQrzaM3mqRR07dsSNGzfg4OAAV1dXfPrpp9DT08OXX34JR0dHqcMjIiINIDpx+u6775QeFxQU4Pr169i1axdCQ0PVFhgRkSpeLVFe1fLk+jpaaGNljITUF9WOKSH1BfIKi2Cgq13tY5FqFixYgOzsbADAkiVL8Pbbb+ONN95AkyZNsH//fomjIyIiTSA6cXr33XdLtb3//vvo0KED9u/fjw8//FAtgRER1SaZTIYQr/bVmqaXV1hU7dEqqhpPT0/FfScnJ9y5cwfPnz+Hubk5R/6IiEgttNR1oN69eyM6OlpdhyMiqnUymQwGutpVvunrqO0jlUTavXs34uPjldosLCyQl5eH3bt3SxQVERFpErX08v/88w82btyIFi1aqONwREREovj4+MDV1RWHDh1Sas/IyICvr69EURERkSYRPVXv1WkPgiDgxYsXMDQ0xNdff63W4IiIxHr1miexqnqNFEkvNDQUEydOxM2bN7F48WKpwyEiIg0jOnFav3690mMtLS00a9YMrq6uMDc3V1dcREREokyYMAHu7u4YMWIEbt26ha+++krqkIiISIOolDiNHDkSO3fuhImJCWQyGcaMGQN9ff2ajo2IiEglxTMhevfujYsXL+Kdd96Bu7s7IiIiJI6MiIg0hUrXOB0/flxR5tXX1xcZGRk1GhQREZEYgiAo7rds2RK//PIL7O3t8dZbb0kYFRERaRKVRpzatm2L4OBgDBgwAIIg4MCBAzAxMSlz20mTJqk1QCIiosqEhITAyMhI8djQ0BDfffcdQkJCcO7cOQkjIyIiTaFS4hQREYGgoCCcOHECMpkMCxYsKHNdDJlMxsSJiIhqXUhISJntXJidiIjURaXEyd3dHb/++iuAl8Ug7t69C0tLyxoNjIiISKz4+HgkJSUhPz9f0SaTyeDl5SVhVEREpAlEV9VLTExEs2bNaiIWIiKiKnnw4AFGjBiBmzdvQiaTKa55Kp4dIZfLpQyPiIg0gOgFcO3s7MqcpkdERCSVmTNnwsHBAWlpaTA0NMTvv/+Oc+fOoUePHoiJiZE6PCIi0gCiR5yIiKjmvboQb3UX9tV0sbGxOH36NJo2bQotLS1oaWnh9ddfR3h4OAICAnD9+nWpQyQionpO9IgTERFRXSOXy2FsbAwAaNq0Kf766y8AL2dJJCQkSBkaERFpCI44ERFRvdexY0fcuHEDDg4OcHV1xaeffgo9PT18+eWXcHR0lDo8IiLSAEyciIio3luwYIFiofYlS5bg7bffxhtvvIEmTZpg//79EkdHRESaQKXEycXFReWCENeuXatWQERERGJ5enoq7js5OeHOnTt4/vw5zM3NWdCIiIjUQqXEafjw4TUcBhERkXpZWFhIHQIREWkQlRKn8lZkJyIiqguys7OxYsUKREdHIy0tDUVFRUrPP3jwQKLIiIhIU/AaJyIiqvemTJmCs2fPYuLEibCxseH0PCIiUjvRiZNcLse6detw4MABJCUlIT8/X+n558+fqy04IiIiVfz44484ceIE+vTpI3UoRESkoUSv4xQaGoq1a9dizJgxyMjIQFBQEEaOHAktLS0sXry4BkIkIiKqmLm5Oa9pIiKiGiU6cdqzZw+2bt2K2bNnQ0dHB2PHjsW2bduwaNEi/Prrr6ID2Lx5M+zt7WFgYABXV1dcunSpwu3T09Ph7+8PGxsb6Ovro3Xr1vjhhx9Evy4REWmOpUuXYtGiRcjJyan2sc6dOwcvLy80b94cMpkMR44cqXD7mJgYyGSyUreUlJRqx0JERHWH6Kl6KSkp6NSpEwDAyMgIGRkZAIC3334bCxcuFHWs/fv3IygoCBEREXB1dcX69evh6emJhIQEWFpalto+Pz8fb731FiwtLfHtt9+iRYsW+PPPP2FmZib2bRARUT336lIZ9+7dg5WVFezt7aGrq6u0rZilMrKzs9GlSxdMnjwZI0eOVHm/hIQEmJiYKB6X1Y8REVH9JTpxeu2115CcnIyWLVuiVatW+Omnn9CtWzdcvnwZ+vr6oo61du1aTJ06Fb6+vgCAiIgInDhxAjt27MD8+fNLbb9jxw48f/4cv/zyi6JTtLe3F/sWiIhIA9TUUhlDhgzBkCFDRO9naWnJH/KISKPF3n8mdQiSEp04jRgxAtHR0XB1dcWMGTMwYcIEbN++HUlJSZg1a5bKx8nPz8fVq1cRHBysaNPS0oKHhwdiY2PL3Of777+Hm5sb/P39cfToUTRr1gzjxo3DvHnzoK2tXeY+eXl5yMvLUzzOzMxUOUYiIqq76tpSGV27dkVeXh46duyIxYsXV1iogn0TEVH9IzpxWrFiheL+mDFjYGdnh19++QXOzs7w8vJS+ThPnz6FXC6HlZWVUruVlRXu3LlT5j4PHjzA6dOnMX78ePzwww+4d+8epk2bhoKCgnI70PDwcISGhqocFxE1bG6tmlR535z8QjVGQvWFjY0NIiIi0KNHD+Tl5WHbtm3o378/Ll68iG7dupW5D/smImqIBEFAXmFR5RuKlJNfiEa62jW+FIXoxCk3NxcGBgaKx71790bv3r3VGlR5ioqKYGlpiS+//BLa2tro3r07Hj9+jFWrVpWbOAUHByMoKEjxODMzE7a2trUSLxERab42bdqgTZs2isfu7u64f/8+1q1bh6+++qrMfdg3EVFDIwgCFh/7HXdTs2rk+PFLPGGoV7NL1Io+uqWlJUaMGIEJEybgzTffhJaW6MJ8AICmTZtCW1sbqampSu2pqamwtrYucx8bGxvo6uoqTctr164dUlJSkJ+fDz09vVL76Ovri772ioiIqDp69eqF8+fPl/s8+yYiamjyCotqLGmqLaITp127dmHv3r149913YWpqijFjxmDChAno0aOHqOPo6emhe/fuiI6OVlzgW1RUhOjoaEyfPr3Mffr06YO9e/eiqKhIkbDdvXsXNjY2ZSZNREREUoiLi4ONjY3UYRAR1UkRE7pDX6dqgy9lcXW0QCPdsusdqJPoiEeMGIGDBw8iNTUVYWFhiI+PR+/evdG6dWssWbJE1LGCgoKwdetW7Nq1C7dv38ZHH32E7OxsRZW9SZMmKRWP+Oijj/D8+XPMnDkTd+/exYkTJxAWFgZ/f3+xb4OIiDRQfn4+EhISUFhY9evNsrKyEBcXh7i4OABAYmIi4uLikJSUBODlNLtJkyYptl+/fj2OHj2Ke/fu4datWwgMDMTp06fZNxERlUNfRwsGutpquxnq6dT49U1AFRKnYsbGxvD19cVPP/2E3377DY0bNxZ9oeuYMWOwevVqLFq0CF27dkVcXByioqIUBSOSkpKQnJys2N7W1hYnT57E5cuX0blzZwQEBGDmzJllli4nIpJSXmERcgvkKt8EQZA65HotJycHH374IQwNDdGhQwdFkjNjxgylokaquHLlClxcXODi4gLg5Y98Li4uWLRoEQAgOTlZcXzgZbI2e/ZsdOrUCf369cONGzfw888/480331TTuyMiorpAJlSxt87NzcX333+PvXv3KpKdsWPHiu6galtmZiZMTU2RkZGhtFAhEVF15eQXov2ik1Xat42VMUK82pf7i1l1qv3VJTX1GTxz5kxcuHAB69evx+DBg/Hbb7/B0dERR48exeLFi3H9+nW1vVZNYN9ERPVBddZxyi2Qw3fnZQBApE9PGKhxal11+kgxn7+ir3E6efIk9u7diyNHjkBHRwfvv/8+fvrpJ/Tt27fKARMRaYJGutroYWeOK3/+LXrfhNQXyCssUmtH0pAcOXIE+/fvR+/evZWSzw4dOuD+/fsSRkZERJqiSgvgvv3229i9ezeGDh0KXV3dmoiLiKjekclkOOjnhpiEJyrvk1dYBL+vr9ZgVA3DkydPYGlpWao9Ozu7Vua9ExGR5hOdOKWmpsLY2LgmYiEiqvdkMhlHjSTQo0cPnDhxAjNmzAAARbK0bds2uLm5SRkaERFpCJUSp8zMTMWcP0EQkJmZWe62nJtNRES1LSwsDEOGDEF8fDwKCwuxYcMGxMfH45dffsHZs2elDo+IiDSASlX1zM3NkZaWBgAwMzODubl5qVtxOxERUW17/fXXERcXh8LCQnTq1Ak//fQTLC0tERsbi+7du0sdHhERaQCVRpxOnz4NCwsLxX3OFyciorqmVatW2Lp1q9RhEBGRhlIpcerXr5/ifv/+/WsqFiIioirx8PDAhAkTMHLkSE4ZJyKiGiF6AVxnZ2csXrwYf/zxR03EQ0REJFqHDh0QHBwMa2trjBo1CkePHkVBQYHUYRERkQYRnThNmzYNJ06cQNu2bdGzZ09s2LABKSkpNREbERGRSjZs2IDHjx/jyJEjaNy4MSZNmgQrKyv861//YnEIIiJSC9HlyGfNmoVZs2bh7t272LNnDzZv3ow5c+ZgwIABmDBhAiZNmlQTcRIRNWhVWa29Oiup10daWloYNGgQBg0ahIiICBw7dgzLly/H9u3bIZfLpQ6PiIjqOdEjTsVat26N0NBQ3L17F//5z3/w5MkT+Pr6qjM2IiIi0VJSUhAREYGVK1fit99+Q8+ePaUOiYiINIDoEaeSLl26hL1792L//v3IzMzEqFGj1BUXERGRyjIzM3Ho0CHs3bsXMTExcHR0xPjx47F//360atVK6vCIiEgDiE6ciqfoffPNN0hMTMTAgQOxcuVKjBw5EkZGRjURIxERUYWsrKxgbm6OMWPGIDw8HD169JA6JCIi0jCiE6fiohD+/v744IMPYGVlVRNxERERqez777/Hm2++CS2tKs9AJyIiqpCoxEkul+OLL77A+++/D3Nz85qKiYiISJS33npL6hCIiEjDiUqctLW1MWPGDHh4eDBxIiIiSXXr1g3R0dEwNzeHi4sLZDJZudteu3atFiMjIiJNJHqqXseOHfHgwQM4ODjURDxEREQqeffdd6Gvr6+4X1HiREREVF2iE6dly5Zhzpw5WLp0Kbp3747GjRsrPW9iYqK24IiIiMoTEhKiuL948WLpAiEiogZBdOI0dOhQAMA777yj9OueIAiQyWRcZJCIiGqdo6MjLl++jCZNlBf9TU9PR7du3fDgwQOJIiMiIk0hOnE6c+ZMTcRBRERUZQ8fPizzh7u8vDz897//lSAiIiLSNKITp379+tVEHERERKJ9//33ivsnT56Eqamp4rFcLkd0dDSvySUiIrUQnTidO3euwuf79u1b5WCIiIjEGD58OABAJpPB29tb6TldXV3Y29tjzZo1EkRGRESaRnTi1L9//1JtJa914jVORERUW4qKigAADg4OuHz5Mpo2bSpxREREpKlEL7H+999/K93S0tIQFRWFnj174qeffqqJGImIiCqUmJjIpImIiGqU6BGnkvPHi7311lvQ09NDUFAQrl69qpbAiIiIxMjOzsbZs2eRlJSE/Px8pecCAgIkioqIiDSF6MSpPFZWVkhISFDX4YiIiFR2/fp1DB06FDk5OcjOzoaFhQWePn0KQ0NDWFpaMnEiIqJqE504/fbbb0qPBUFAcnIyVqxYga5du6orLiKiBievsEhxX19HS+n6UarYrFmz4OXlhYiICJiamuLXX3+Frq4uJkyYgJkzZ0odHhERaQDRiVPXrl0hk8kgCIJSe+/evbFjxw61BUZE1ND4ff2/qc5trIwR4tWeyZOK4uLi8MUXX0BLSwva2trIy8uDo6MjPv30U3h7e2PkyJFSh0hERPWc6MQpMTFR6bGWlhaaNWsGAwMDtQVFRNRQ6OtooY2VMRJSXyi1J6S+QF5hEQx0tSWKrH7R1dWFltbLekeWlpZISkpCu3btYGpqikePHkkcHRERaQLRiZOdnV2ptvT0dCZORERVIJPJEOLVXjFNL6+wSGnkiVTj4uKCy5cvw9nZGf369cOiRYvw9OlTfPXVV+jYsaPU4RERkQYQXY585cqV2L9/v+Lx6NGjYWFhgRYtWuDGjRtqDY6IqCGQyWQw0NWGga429HVEfywTgLCwMNjY2AAAli9fDnNzc3z00Ud48uQJvvzyS4mjIyIiTSB6xCkiIgJ79uwBAJw6dQqnTp1CVFQUDhw4gLlz53ItJyIiqnU9evRQ3Le0tERUVJSE0RARaYbY+8+kDqFOEZ04paSkwNbWFgBw/PhxjB49GoMGDYK9vT1cXV3VHiAREREREdUtgiAoVYOtjJht6yrRiZO5uTkePXoEW1tbREVFYdmyZQBenjy5XK72AImIiMri4uKictXBa9eu1XA0REQNhyAIWHzsd9xNzZI6lFolOnEaOXIkxo0bB2dnZzx79gxDhgwB8HLxQScnJ7UHSEREVJbhw4dLHQIRUYOUV1hU5aSpjZVxvb2eV3TitG7dOtjb2+PRo0f49NNPYWRkBABITk7GtGnT1B4gERFRWUJCQqQOgYiowYuY0F1UIlSfF3gXnTjp6upizpw5pdpnzZqlloCIiIiIiKh+0NfRajBrDopOnIiIiOoaLa2Kf8HkNbhERFRdTJyIiKje++6775QeFxQU4Pr169i1axdCQ0MlioqIiDQJEyciIqr33n333VJt77//Pjp06ID9+/fjww8/lCAqIiLSJPWzpAUREZEKevfujejoaKnDICIiDcDEiYiINNI///yDjRs3okWLFlKHQkREGkClqXrm5uYqlw18/vx5tQIiIiIS69V+ShAEvHjxAoaGhvj6668ljIyIiDSFSonT+vXrFfefPXuGZcuWwdPTE25ubgCA2NhYnDx5EgsXLqyRIImIiCpSsp8CXlbZa9asGVxdXWFubi5NUEREpFFUSpy8vb0V99977z0sWbIE06dPV7QFBARg06ZN+Pnnn7meExER1bqS/RQREVFNEF1V7+TJk1i5cmWp9sGDB2P+/PlqCYqIiKov9v4ztR/TrVUTtR9TXXJzc/Hbb78hLS0NRUVFSs+98847EkVFRESaQnTi1KRJExw9ehSzZ89Waj969CiaNKm7HSoREWmuqKgoTJw4Ec+elU4WZTIZF8AlIqJqE11VLzQ0FPPmzYOXlxeWLVuGZcuWwcvLC/Pnz6/yIoObN2+Gvb09DAwM4OrqikuXLqm03759+yCTyTB8+PAqvS4REWmGGTNmYPTo0UhOTkZRUZHSjUkTERGpg+jEycfHBxcuXICJiQkOHz6Mw4cPw8TEBOfPn4ePj4/oAPbv34+goCCEhITg2rVr6NKlCzw9PZGWllbhfg8fPsScOXPwxhtviH5NIiLSLKmpqQgKCoKVlVW1j3Xu3Dl4eXmhefPmkMlkOHLkSKX7xMTEoFu3btDX14eTkxN27txZ7TiIiKhuqdI6Tq6urtizZw+uXbuGa9euYc+ePXB1da1SAGvXrsXUqVPh6+uL9u3bIyIiAoaGhtixY0e5+8jlcowfPx6hoaFwdHSs0usSEZHmeP/99xETE6OWY2VnZ6NLly7YvHmzStsnJiZi2LBhGDBgAOLi4hAYGIgpU6bg5MmTaomHiIjqBtHXOAHA/fv3ERkZiQcPHmD9+vWwtLTEjz/+iJYtW6JDhw4qHyc/Px9Xr15FcHCwok1LSwseHh6IjY0td78lS5bA0tISH374If7zn/9U+Bp5eXnIy8tTPM7MzFQ5PiIiqh82bdqEUaNG4T//+Q86deoEXV1dpecDAgJUPtaQIUMwZMgQlbePiIiAg4MD1qxZAwBo164dzp8/j3Xr1sHT01Pl4xARUd0mOnE6e/YshgwZgj59+uDcuXNYtmwZLC0tcePGDWzfvh3ffvutysd6+vQp5HJ5qakVVlZWuHPnTpn7nD9/Htu3b0dcXJxKrxEeHl7la6+IiKh++Oabb/DTTz/BwMAAMTExSovhymQyUYmTWLGxsfDw8FBq8/T0RGBgYLn78Ec9IqL6R/RUvfnz52PZsmU4deoU9PT0FO0DBw7Er7/+qtbgXvXixQtMnDgRW7duRdOmTVXaJzg4GBkZGYrbo0ePajRGIiKqfZ988glCQ0ORkZGBhw8fIjExUXF78OBBjb52SkpKmT8AZmZm4p9//ilzn/DwcJiamiputra2NRojERFVn+gRp5s3b2Lv3r2l2i0tLfH06VNRx2ratCm0tbWRmpqq1J6amgpra+tS29+/fx8PHz6El5eXoq14rQ4dHR0kJCSgVatWSvvo6+tDX19fVFxERFS/5OfnY8yYMdDSqtKlu7UuODgYQUFBiseZmZlMnoiI6jjRPYyZmRmSk5NLtV+/fh0tWrQQdSw9PT10794d0dHRiraioiJER0fDzc2t1PZt27bFzZs3ERcXp7i98847igty2ekQkSbJKyxCboEcuQVyCIIgdTh1mre3N/bv3y/Ja1tbW5f5A6CJiQkaNWpU5j76+vowMTFRuhERUd0mesTpgw8+wLx583Dw4EHIZDIUFRXhwoULmDNnDiZNmiQ6gKCgIHh7e6NHjx7o1asX1q9fj+zsbPj6+gIAJk2ahBYtWiA8PBwGBgbo2LGj0v5mZmYAUKqdiKi+8/v6quJ+GytjhHi1V7p2h/5HLpfj008/xcmTJ9G5c+dSxSHWrl1bY6/t5uaGH374Qant1KlTZf4ASERE9ZfoxCksLAz+/v6wtbWFXC5H+/btIZfLMW7cOCxYsEB0AGPGjMGTJ0+waNEipKSkoGvXroiKilLMF09KSqo3Uy+IiKpLX0cLbayMkZD6Qqk9IfUF8gqLYKCrLVFkddvNmzfh4uICALh165bSc2KTzaysLNy7d0/xODExEXFxcbCwsEDLli0RHByMx48fY/fu3QAAPz8/bNq0CR9//DEmT56M06dP48CBAzhx4kQ13xUREdUlMqGK8z8ePXqEmzdvIisrCy4uLnB2dlZ3bDUiMzMTpqamyMjI4NQIIqoRsfefVWt/QRCQV/jy+s28wiLFyFOkT0/JEye3Vk2qtX99+AyOiYnBgAEDSrV7e3tj586d8PHxwcOHD5XWjYqJicGsWbMQHx+P1157DQsXLhS1KHx9OC9E1PCU15/lFsjhu/MygPrfN4n5/BU94rRkyRLMmTMHtra2StcU/fPPP1i1ahUWLVokPmIiIlKQyWSSd0INWf/+/Su8pmznzp1l7nP9+vUajIqIiKQmOnEKDQ2Fn58fDA0NldpzcnIQGhrKxImIiGrdgAEDKpySd/r06VqMhoiINJHoxEkQhDI7pxs3bsDCwkItQREREYnRtWtXpccFBQWIi4vDrVu34O3tLU1QRESkUVROnMzNzSGTySCTydC6dWul5EkulyMrKwt+fn41EiQREVFF1q1bV2b74sWLkZWVVcvREBGRJlI5cVq/fj0EQcDkyZMRGhoKU1NTxXN6enqwt7dn6VUiIqpTJkyYgF69emH16tVSh0JEVK+9WrioIVI5cSqe6uDg4AB3d/dSa2QQERHVNbGxsTAwMJA6DCKiek0QBCw+9jvupjbsEXzR1zj169dPcT83Nxf5+flKz7OMKhE1dFUpi1rdEuYN3ciRI5UeC4KA5ORkXLlyBQsXLpQoKiIizZBXWFRm0tTGyhj6Og1nvVXRiVNOTg4+/vhjHDhwAM+ele7o5XK5WgIjIiJSVcnp4wCgpaWFNm3aYMmSJRg0aJBEURER1R3q+oEuYkJ3RbKkr6MlepHx+kx04jR37lycOXMGW7ZswcSJE7F582Y8fvwYX3zxBVasWFETMRIREVUoMjJS6hCIiBoEfR2tBrvWoOjE6dixY9i9ezf69+8PX19fvPHGG3BycoKdnR327NmD8ePH10ScRERE5bp8+TKKiorg6uqq1H7x4kVoa2ujR48eEkVGRESaQvSkxOfPn8PR0RHAy+uZnj9/DgB4/fXXce7cOfVGR0REpAJ/f388evSoVPvjx4/h7+8vQURERKRpRCdOjo6OSExMBAC0bdsWBw4cAPByJMrMzEytwREREakiPj4e3bp1K9Xu4uKC+Ph4CSIiIiJNIzpx8vX1xY0bNwAA8+fPx+bNm2FgYIBZs2Zh7ty5ag+QiIioMvr6+khNTS3VnpycDB0d0bPSiYiIShHdm8yaNUtx38PDA3fu3MHVq1fh5OSEzp07qzU4IiIiVQwaNAjBwcE4evSoosJeeno6/v3vf+Ott96SODoiItIE1f4Zzs7ODnZ2duqIhYiIqEpWr16Nvn37ws7ODi4uLgCAuLg4WFlZ4auvvpI4OiIi0gRVSpwuX76MM2fOIC0tDUVFRUrPrV27Vi2BERERqapFixb47bffsGfPHty4cQONGjWCr68vxo4dC11dXanDIyIiDSA6cQoLC8OCBQvQpk0bWFlZKS161ZAWwCIiorqlcePG+Ne//iV1GEREpKFEJ04bNmzAjh074OPjUwPhEBERERER1T2iq+ppaWmhT58+NRELERERERFRnSQ6cZo1axY2b95cE7EQERERERHVSaKn6s2ZMwfDhg1Dq1at0L59+1IX3R4+fFhtwREREREREdUFohOngIAAnDlzBgMGDECTJk1YEIKIiIiIiDSe6MRp165dOHToEIYNG1YT8RAREanE3Nxc5R/vnj9/XsPREBGRphOdOFlYWKBVq1Y1EQsREZHK1q9fr7j/7NkzLFu2DJ6ennBzcwMAxMbG4uTJk1i4cKFEERIRkSYRnTgtXrwYISEhiIyMhKGhYU3EREREZcgrLKp8oxL0dbQ0ejq1t7e34v57772HJUuWYPr06Yq2gIAAbNq0CT///DNmzZolRYhERKRBRCdOGzduxP3792FlZQV7e/tSxSGuXbumtuCIiOh//L6+Kmr7NlbGCPFqr9HJU7GTJ09i5cqVpdoHDx6M+fPnSxARERFpGtGJ0/Dhw2sgDCIiKou+jhbaWBkjIfWF6H0TUl8gr7AIBrraNRBZ3dKkSRMcPXoUs2fPVmo/evQomjRpIlFURESkSUQnTiEhITURBxERlUEmkyHEq72oaXp5hUWiR6fqu9DQUEyZMgUxMTFwdXUFAFy8eBFRUVHYunWrxNEREdU/giAo+h6xU8U1lejEiYiIapdMJmsQo0bV4ePjg3bt2mHjxo2K9QTbtWuH8+fPKxIpIiJSjSAIWHzsd9xNzZI6lDpFpcTJwsICd+/eRdOmTSst/8qSr0REJAVXV1fs2bNH6jCIiOq9vMKiMpOmNlbG0NfRkiCiukGlxGndunUwNjZW3G8IFxoTEVH9cv/+fURGRuLBgwdYv349LC0t8eOPP6Jly5bo0KGD1OEREdVLERO6K5IlTa/WWhmVEqeSJV99fHxqKhYiIqIqOXv2LIYMGYI+ffrg3LlzWLZsGSwtLXHjxg1s374d3377rdQhEhHVS/o6Wpwu/v+JHmvT1tZGWlpaqfZnz55BW5snlYiIat/8+fOxbNkynDp1Cnp6eor2gQMH4tdff5UwMiIi0hSiEydBEMpsz8vLU+qsiIiIasvNmzcxYsSIUu2WlpZ4+vSpBBEREZGmUbmq3saNGwG8rO60bds2GBkZKZ6Ty+U4d+4c2rZtq/4IiYiIKmFmZobk5GQ4ODgotV+/fh0tWrSQKCoiItIkKidO69atA/ByxCkiIkJpWp6enh7s7e0RERGh/giJiIgq8cEHH2DevHk4ePAgZDIZioqKcOHCBcyZMweTJk2SOjwiItIAKidOiYmJAIABAwbg8OHDMDc3r7GgiIiIxAgLC4O/vz9sbW0hl8vRvn17yOVyjBs3DgsWLJA6PCIi0gCiF8A9c+aM0mO5XI6bN2/Czs6OyRQREUlCT08PW7duxaJFi3Dz5k1kZWXBxcUFzs7OUodGREQaQnRxiMDAQGzfvh3Ay6Spb9++6NatG2xtbRETE6Pu+IiIiCq1ZMkS5OTkwNbWFkOHDsXo0aPh7OyMf/75B0uWLJE6PCIi0gCiE6eDBw+iS5cuAIBjx47h4cOHuHPnDmbNmoVPPvlE7QESERFVJjQ0FFlZpVe5z8nJQWhoqAQRERGRphGdOD179gzW1tYAgB9++AGjRo1C69atMXnyZNy8eVPtARIREVVGEIQyV7O/ceMGLCwsJIiIiIg0jehrnKysrBAfHw8bGxtERUVhy5YtAF7+qscFcImIqDaZm5tDJpNBJpOhdevWSsmTXC5HVlYW/Pz8JIyQiIg0hejEydfXF6NHj4aNjQ1kMhk8PDwAABcvXuQ6TkREVKvWr18PQRAwefJkhIaGwtTUVPFc8VIZbm5uEkZIRESaQnTitHjxYnTs2BGPHj3CqFGjoK+vDwDQ1tbG/Pnz1R4gERFReby9vQEADg4OcHd3h66ursQRERGRphKdOAHA+++/X6qtuPMiIiLx3Fo1UXoce/+ZRJHUT/369VPcz83NRX5+vtLzJiYmtR0SERFpGJWLQwwdOhQZGRmKxytWrEB6erri8bNnz9C+fXu1BkdERKSKnJwcTJ8+HZaWlmjcuDHMzc2VbkRERNWlcuJ08uRJ5OXlKR6HhYXh+fPniseFhYVISEhQb3REREQqmDt3Lk6fPo0tW7ZAX18f27ZtQ2hoKJo3b47du3dLHR4RkaQEQUBugVzlW15hkdQh10kqT9UTBKHCx0RERFI5duwYdu/ejf79+8PX1xdvvPEGnJycYGdnhz179mD8+PFSh0hEJAlBEPB+RCyu/vm31KHUe6LXcaoJmzdvhr29PQwMDODq6opLly6Vu+3WrVvxxhtvKKZfeHh4VLg9ERFpvufPn8PR0RHAy+uZimdEvP766zh37pzo44npl3bu3KkoiV58MzAwqNobISJSs38K5FVOmtpYGUNfp06kC3WCyiNOxZ3Bq23VtX//fgQFBSEiIgKurq5Yv349PD09kZCQAEtLy1Lbx8TEYOzYsXB3d4eBgQFWrlyJQYMG4ffff0eLFi2qHQ8REdU/jo6OSExMRMuWLdG2bVscOHAAvXr1wrFjx2BmZibqWGL7JeBlslZyuro6+kciInWLmNBdVCKkr6PFz7MSRE3V8/HxUZQfz83NhZ+fHxo3bgwAStc/ibF27VpMnToVvr6+AICIiAicOHECO3bsKLO8+Z49e5Qeb9u2DYcOHUJ0dDQmTZpUpRiIiKh+8/X1xY0bN9CvXz/Mnz8fXl5e2LRpEwoKCrB27VpRxxLbLwEvEyVra+tqvw8iopqkr6MFA11tqcOot1ROnF4tNz5hwoRS24hNXPLz83H16lUEBwcr2rS0tODh4YHY2FiVjpGTk4OCggJYWFiU+XxeXp5SUpeZmSkqRiIiqvtmzZqluO/h4YE7d+7g6tWrcHJyQufOnVU+TlX7paysLNjZ2aGoqAjdunVDWFgYOnToUO727JuIiOoflROnyMhItb/406dPIZfLYWVlpdRuZWWFO3fuqHSMefPmoXnz5vDw8Cjz+fDwcISGhlY7ViKi+kjdlZFy8gvRSFe7zk/dsLOzg52dnej9qtIvtWnTBjt27EDnzp2RkZGB1atXw93dHb///jtee+21Mvdh30REYlRnbb/cArkaI2nYqrQAbl2xYsUK7Nu3DzExMeVeiBscHIygoCDF48zMTNja2tZWiEREkvL7+qrajxm/xBOGenWv+7h8+TLOnDmDtLQ0FBUpJ4xip+uJ4ebmBjc3N8Vjd3d3tGvXDl988QWWLl1a5j7sm4iI6h9Je76mTZtCW1sbqampSu2pqamVzhVfvXo1VqxYgZ9//rnCaRj6+vqK67KIiBoCfR0ttLEyRkLqC6lDqTVhYWFYsGAB2rRpAysrK6URMTGjY9Xpl4rp6urCxcUF9+7dK3cb9k1ERPWPpImTnp4eunfvjujoaAwfPhwAUFRUhOjoaEyfPr3c/T799FMsX74cJ0+eRI8ePWopWiKi+kEmkyHEq32NLGDo6miBRnXwwuINGzZgx44d8PHxqdZxqtovlSSXy3Hz5k0MHTq0WrEQEVHdIvlci6CgIHh7e6NHjx7o1asX1q9fj+zsbEU1o0mTJqFFixYIDw8HAKxcuRKLFi3C3r17YW9vj5SUFACAkZERjIyMJHsfRER1iUwmq5HKSXVxih7wsoBDnz591HIssf3SkiVL0Lt3bzg5OSE9PR2rVq3Cn3/+iSlTpqglHiIiqhsk7wHHjBmDJ0+eYNGiRUhJSUHXrl0RFRWluDA3KSkJWlr/qze/ZcsW5Ofn4/3331c6TkhICBYvXqy2uORyOQoKCtR2PGrYdHV1oa1d936lJ9IUs2bNwubNm7F+/fpqH0tsv/T3339j6tSpSElJgbm5Obp3745ffvkF7du3r3YsRESaxq1VE6lDqDKZIAiC1EHUpszMTJiamiIjIwMmJialnhcEASkpKUhPT6/94EijmZmZwdraus5XI6O6oToVlGpSdTu8yj6Dq6qoqAjDhg3D3bt30b59e+jq6io9f/jwYbW9Vk2oqfNCRJqhulX1fHdeBgBE+vSUfB2nupY4ifn8lXzEqa4pTposLS1haGjIL7lUbYIgICcnB2lpaQAAGxsbiSMi0jwBAQE4c+YMBgwYgCZNmvCzm4iI1I6JUwlyuVyRNDVpUreyYarfGjVqBABIS0uDpaUlp+0RqdmuXbtw6NAhDBs2TOpQiIhIQ2lVvknDUXxNk6GhocSRkCYq/rvitXNE6mdhYYFWrVpJHQYREWkwJk5l4BQPqgn8uyKqOYsXL0ZISAhycnKkDoWIiDQUp+oREVG9t3HjRty/fx9WVlawt7cvVRzi2rVrEkVGRESagokTqWzx4sU4cuQI4uLipA6lQv3790fXrl3VUpZYVUeOHMGcOXOQmJiIGTNm1OprExEUi9USERHVFCZOGsLHxwfp6ek4cuSIou3bb7/FhAkTsHz5csyePVu64Crw8OFDODg4QEtLC0lJSWjRooXiueTkZNja2kIulyMxMRH29vYqHfPw4cOlfm2uaf/3f/8HX19fBAQEwNjYuFZfm4heruVHRERUk5g4aaht27bB398fERERitXu67IWLVpg9+7dCA4OVrTt2rULLVq0QFJSkqhjWVhYVCsWuVwOmUymtMBlRbKyspCWlgZPT080b968Wq9NREREmqW66/IJgoC8wqIq71+dfUkZi0NUQhAE5OQXSnKr6trEn376KWbMmIF9+/Ypkqa1a9eiU6dOaNy4MWxtbTFt2jRkZWUp9tm5cyfMzMxw5MgRODs7w8DAAJ6ennj06FG5r3P58mW89dZbaNq0KUxNTdGvX79S1xHIZDJs27YNI0aMgKGhIZydnfH999+XOpa3tzciIyOV2iIjI+Ht7V1q21u3bmHIkCEwMjKClZUVJk6ciKdPnyqe79+/PwIDAxWP8/LyMGfOHLRo0QKNGzeGq6srYmJiSr3377//Hu3bt4e+vn6pZO3GjRsYMGAAjI2NYWJigu7du+PKlSuIiYlRjDANHDgQMplM6dhEVHMsLCwU//bNzc1hYWFR7o2IqD4SBAGLj/0O352Xq3zz+/qq1G9DY3DEqRL/FMjRftFJSV47foknDPXE/S+aN28ePv/8cxw/fhxvvvmmol1LSwsbN26Eg4MDHjx4gGnTpuHjjz/G559/rtgmJycHy5cvx+7du6Gnp4dp06bhgw8+wIULF8p8rRcvXsDb2xufffYZBEHAmjVrMHToUPzxxx9K09VCQ0Px6aefYtWqVfjss88wfvx4/Pnnn0pfZt555x1ERETg/PnzeP3113H+/Hn8/fff8PLywtKlSxXbpaenY+DAgZgyZQrWrVuHf/75B/PmzcPo0aNx+vTpMuOcPn064uPjsW/fPjRv3hzfffcdBg8ejJs3b8LZ2Vnx3leuXIlt27ahSZMmsLS0VDrG+PHj4eLigi1btkBbWxtxcXHQ1dWFu7s7EhIS0KZNGxw6dAju7u78kkZUS9atW6f4rFm3bh0rVxJRnVDdEaaS8gqLcDc1q/INVdDGyhj6OhwzqQ4mThrkxx9/xNGjRxEdHY2BAwcqPVdyBMbe3h7Lli2Dn5+fUuJUUFCATZs2wdXVFcDLqXLt2rXDpUuX0KtXr1Kv9+prfPnllzAzM8PZs2fx9ttvK9p9fHwwduxYAEBYWBg2btyIS5cuYfDgwYptdHV1MWHCBOzYsQOvv/46duzYgQkTJpS6VmnTpk1wcXFBWFiYom3Hjh2wtbXF3bt30bp1a6Xtk5KSEBkZiaSkJMU0ujlz5iAqKgqRkZGK4xQUFODzzz9Hly5dyjy3SUlJmDt3Ltq2bQsAioQLgCLJsrCwgLW1dZn7E5H6lRyR9vHxkS4QIqJaEDGhe7USH30dLbX/wOTWqolaj1fXMXGqRCNdbcQv8ZTstcXo3Lkznj59ipCQEPTq1QtGRkaK537++WeEh4fjzp07yMzMRGFhIXJzc5GTk6NYmFVHRwc9e/ZU7NO2bVuYmZnh9u3bZSZOqampWLBgAWJiYpCWlga5XI6cnJxS09w6d+6suN+4cWOYmJggLS2t1PEmT54Md3d3hIWF4eDBg4iNjUVhYaHSNjdu3MCZM2eU3lux+/fvl0qcbt68CblcXqo9Ly8PTZr87x+7np6eUpyvCgoKwpQpU/DVV1/Bw8MDo0aN4mKbRHWItrY2kpOTS40WP3v2DJaWlpDL5RJFRkSkHvo6WjAQ+d2Q1IuJUyVkMpno6XJSadGiBb799lsMGDAAgwcPxo8//ghjY2M8fPgQb7/9Nj766CMsX74cFhYWOH/+PD788EPk5+crEiexvL298ezZM2zYsAF2dnbQ19eHm5sb8vPzlbZ7ddRIJpOhqKj0hYqdOnVC27ZtMXbsWLRr1w4dO3YsVfo8KysLXl5eWLlyZan9bWxsSrVlZWVBW1sbV69ehba28odNyeSrUaNGFf4Ks3jxYowbNw4nTpzAjz/+iJCQEOzbtw8jRowodx8iqj3lXROal5cHPT29Wo6GiIg0Uf3ICEhldnZ2OHv2rCJ5ioqKwtWrV1FUVIQ1a9YoKsUdOHCg1L6FhYW4cuWKYnQpISEB6enpaNeuXZmvdeHCBXz++ecYOnQoAODRo0dKRRqqYvLkyZg2bRq2bNlS5vPdunXDoUOHYG9vDx2dyv98XVxcIJfLkZaWhjfeeKNasbVu3RqtW7fGrFmzMHbsWERGRjJxIpLYxo0bAfyvEE3JH0TkcjnOnTunmGJLRFQTKrqmiRXxNAsTJw1ka2uLmJgYDBgwAJ6entiyZQsKCgrw2WefwcvLCxcuXEBERESp/XR1dTFjxgxs3LgROjo6mD59Onr37l3mND3g5XU+X331FXr06IHMzEzMnTsXjRo1qlbsU6dOxahRo2BmZlbm8/7+/ti6dSvGjh2Ljz/+GBYWFrh37x727duHbdu2lRpVat26NcaPH49JkyZhzZo1cHFxwZMnTxAdHY3OnTtj2LBhZb5OcHAwHj9+jN27d+Off/7B3Llz8f7778PBwQH//e9/cfnyZbz33nvVeq9EVH3r1q0D8PLLSUREhNJngJ6eHuzt7cv8vCMiqmnFFfHUVdyBpMfSGhrqtddeQ0xMDJ4+fQo/Pz8sXrwYK1euRMeOHbFnzx6Eh4eX2sfQ0BDz5s3DuHHj0KdPHxgZGWH//v3lvsb27dvx999/o1u3bpg4cSICAgJKXV8glo6ODpo2bVruaFLz5s1x4cIFyOVyDBo0CJ06dUJgYCDMzMzKXXcpMjISkyZNwuzZs9GmTRsMHz4cly9fRsuWLcuNIzk5WXGtlra2Np49e4ZJkyahdevWGD16NIYMGYLQ0NBqvVciqr7ExEQkJiaiX79+uHHjhuJxYmIiEhIScPLkSUXBGyKimiYIAnIL5MgtkCMzt5AV8TSMTKjqYkH1VGZmJkxNTZGRkQETExOl53Jzc5GYmAgHBwcYGBhIFKE0du7cicDAQKSnp0sdisZqyH9fVH3qLG9bHdWtoFTRZ7A6yeVy3Lx5E3Z2djA3N6+x11GX2jovRFR95X0eVzTCVBcr4qmDJlTVE/P5y9SViIjqvcDAQGzfvh3Ay6Spb9++6Natm2LqMhFRTStvzaU2VsYwMdCBga52lW91MWlqiHiNExER1XsHDx7EhAkTAADHjh3Dw4cPcefOHXz11Vf45JNPyl3Im4ioOkoWfyhZyKHkCFNdHS0i8Zg4EYCXi0dyAUkiqq+ePXumWID6hx9+wKhRo9C6dWtMnjwZGzZskDg6ItJEFU3N09Q1lzRhal51MHEiIqJ6z8rKCvHx8bCxsUFUVJRiSYOcnJxS1TaJiMQoeU3TqyNM5U3NYyEHzcTEiYiI6j1fX1+MHj0aNjY2kMlk8PDwAABcvHiR6zgRkVqoWvyBU/M0FxMnIiKq9xYvXoyOHTvi0aNHGDVqFPT19QG8XE5g/vz5EkdHRJqgsuIPTJY0HxMnIiLSCO+//36pNm9vbwkiIaL6orKlHlj8gUpi4kRERPXW0KFD8c0338DU1BQAsGLFCvj5+cHMzAzAy6IRb7zxBuLj4yWMkojqI3UVf6grBRXqynqA9RkTJyIiqrdOnjyJvLw8xeOwsDCMHj1akTgVFhYiISFBouiIqD6ravGHupIokfoxcSIionpLEIQKHxMRlVRehbyycGoek8BXMXFSUW0Pb4r9Q/Xx8UF6ejqOHDkCAEhJScHy5ctx4sQJPH78GJaWlujatSsCAwPx5ptvlnuca9euYd68ebh8+TK0tbXx3nvvYe3atTAyMlJsU9aHxTfffIMPPvgAAHD9+nVMnjwZf/zxBwYMGIBdu3bBwsICwMtff11dXbFlyxb06tVL1HskIiIiUoeKpuGVpaKpeUwuGg4WmddADx8+RPfu3XH69GmsWrUKN2/eRFRUFAYMGAB/f/9y9/vrr7/g4eEBJycnXLx4EVFRUfj999/LXBg3MjISycnJitvw4cMVz02ZMgUDBw7EtWvXkJGRgbCwMMVza9asQZ8+fZg0EZFayGSyUj/mNJRfgomocrH3nyluv9x7itwCOXIL5MjMLVQ5aeK6TFSMI04aaNq0aZDJZLh06RIaN26saO/QoQMmT55c7n7Hjx+Hrq4uNm/eDC2tlx8QERER6Ny5M+7duwcnJyfFtmZmZrC2ti7zOLdv38aePXvQunVrjB07FsePHwcAPHjwANu3b8fVq1fV8TaJiCAIAnx8fBTlx3Nzc+Hn56f47Ct5/RMRab6KpuItPvY7/nyWU2qfktPwyvLq1DyOMDVcTJw0zPPnzxEVFYXly5crJU3Fii+YLkteXh709PQUSRMANGrUCABw/vx5pcTJ398fU6ZMgaOjI/z8/ODr66v4UOnSpQtOnToFJycnREdHo3PnzgAAPz8/fPrppzA2NlbHWyUiKlVufMKECaW2mTRpUm2FQ0S16NXLKFRNlEriGkwkBhMnDXPv3j0IgoC2bduK3nfgwIEICgrCqlWrMHPmTGRnZysWjkxOTlZst2TJEgwcOBCGhob46aefMG3aNGRlZSEgIAAAsG3bNkybNg2rV69Gnz59EBwcjK+++gqGhobo2bMnPD09cf/+fXzwwQdYtmyZet44ETVIkZGRUodARHWAqtcs2TUxxGKvDorHqhR60JQRJk15H1Ji4qRhVK0o5efnh6+//lrxOCsrCx06dMCuXbsQFBSE4OBgaGtrIyAgAFZWVkqjUAsXLlTcd3FxQXZ2NlatWqVInDp06ICzZ88qtnn27BlCQkJw7tw5zJgxA+7u7jh8+DB69uwJV1dXeHl5VfdtExERkYaraBpeeaXDK0uUmEyQGEycNIyzszNkMhnu3LlT4XZLlizBnDlzSrWPGzcO48aNQ2pqKho3bgyZTIa1a9fC0dGx3GO5urpi6dKlyMvLU1xnUFJQUBACAwPx2muvISYmBsuWLUPjxo0xbNgwxMTEMHEiIiKiUqpyvRLQcEuHU81j4qRhLCws4Onpic2bNyMgIKDUdU7p6ekwMzODpaUlLC0tyz2OlZUVAGDHjh0wMDDAW2+9Ve62cXFxMDc3LzNpio6Oxu3btxXTaeRyOQoKCgBA8V8iIiKikgRBQG6BXPFYleuVgMqvWeIIE1UHEycNtHnzZkXJ7yVLlqBz584oLCzEqVOnsGXLFty+fbvcfTdt2gR3d3cYGRnh1KlTmDt3LlasWKEoKnHs2DGkpqaid+/eMDAwwKlTpxAWFlbm6FVubi6mT5+Ob775RjHVr0+fPti8eTP8/f1x6NAhrF27tkbOAREREdUvxSNMVb1eCeBUPKpZTJxUVJ/+4Tk6OuLatWtYvnw5Zs+ejeTkZDRr1gzdu3fHli1bKtz30qVLCAkJQVZWFtq2bYsvvvgCEydOVDxfXK581qxZEAQBTk5OWLt2LaZOnVrqWKGhoRg2bBi6du2qaNu4cSPGjRuHvn37Yvz48XjvvffU9r6JiIiofhAEATEJT8p8rqrXKwH16/sa1T8yQdVqAhoiMzMTpqamyMjIgImJidJzubm5SExMhIODAwwMDCSKkDQV/76oOl4tuyuV6n4pqegzuCHjeSFN92qipOrUu4quV2KSROog5vOXI05EREREVGMEQcD7EbG4+uffovZ79XolJkokNSZORERERKQ2r44u5RUWlZk0lXWNUkn92zRjRTyqU5g4EREREZEoNVEqnCNKVNcxcSIiIiKiclVUyAGoeqlwJkpU3zBxKkMDq5dBtYR/V1QdVfmCUVcKShBR3SUIAv4psV5SsYsPnivuq5oYlVTWNDxOvaP6jolTCbq6ugCAnJwcNGrUSOJoSNPk5LzsdIr/zoiIiGpaWYlRdZOisryaKDFJIk3ExKkEbW1tmJmZIS0tDQBgaGjIf/RUbYIgICcnB2lpaTAzM4O2trbUIRERkYZQ1zS6yrCQAxETp1Ksra0BQJE8EamLmZmZ4u+LiIioMpUlRYB6EqPKkiKAiRERwMSpFJlMBhsbG1haWqKgoEDqcEhD6OrqcqSJiKgBUeXaocpwtIiobmHiVA5tbW1+0SUiImoAyktyXiVF0qMKJkZEtaNOJE6bN2/GqlWrkJKSgi5duuCzzz5Dr169yt3+4MGDWLhwIR4+fAhnZ2esXLkSQ4cOrcWIiYhIk7FfqjtUTWpeVVeTHLE4jY6o7pA8cdq/fz+CgoIQEREBV1dXrF+/Hp6enkhISIClpWWp7X/55ReMHTsW4eHhePvtt7F3714MHz4c165dQ8eOHSV4B0REpEnYL1VdVZOcV9VE1TcpqJL0vMrV0ULpcSNdbSZFRHWETJB4cRlXV1f07NkTmzZtAgAUFRXB1tYWM2bMwPz580ttP2bMGGRnZ+P48eOKtt69e6Nr166IiIio9PUyMzNhamqKjIwMmJiYqO+NEBHVMTWxjlN1F6ysD5/Btd0vAf87L8lPntXZ81IZQQBGRcQiPjlT6lCqRB1JzquY9BDVfWL6JUlHnPLz83H16lUEBwcr2rS0tODh4YHY2Ngy94mNjUVQUJBSm6enJ44cOVLm9nl5ecjLy1M8zsjIAPDyJBERabLsFy/UfszMzOqtQ1b82VtXF4SujX4JKL9v6rHoGLT0DavxDjSTrUUjBA9pJ2qfHvbmoravSpJTmFvxSNiLXFGHIyIJiOmXJE2cnj59CrlcDisrK6V2Kysr3Llzp8x9UlJSytw+JSWlzO3Dw8MRGhpaqt3W1raKURMRUXW9ePECpqamUodRSm30S0D5fdPjLT7ig24AHgHwWiR1FESkyVTplyS/xqmmBQcHK/0SmJ6eDjs7OyQlJdXJTru+yczMhK2tLR49elRvp5fUJTyf6sXzqV7qOJ+CIODFixdo3ry5mqOrX9g31Sz+21cvnk/14vlUr+qeTzH9kqSJU9OmTaGtrY3U1FSl9tTU1HIXCrW2tha1vb6+PvT19Uu1m5qa8o9VjUxMTHg+1YjnU714PtWruuezLicGtdEvAeybagv/7asXz6d68XyqV3XOp6r9klaVjq4menp66N69O6KjoxVtRUVFiI6OhpubW5n7uLm5KW0PAKdOnSp3eyIiIlWxXyIiovJIPlUvKCgI3t7e6NGjB3r16oX169cjOzsbvr6+AIBJkyahRYsWCA8PBwDMnDkT/fr1w5o1azBs2DDs27cPV65cwZdffinl2yAiIg3BfomIiMoieeI0ZswYPHnyBIsWLUJKSgq6du2KqKgoxYW2SUlJ0NL638CYu7s79u7diwULFuDf//43nJ2dceTIEZXXytDX10dISEiZUyRIPJ5P9eL5VC+eT/VqKOeztvsloOGc29rC86lePJ/qxfOpXrV5PiVfx4mIiIiIiKiuk/QaJyIiIiIiovqAiRMREREREVElmDgRERERERFVgokTERERERFRJRpc4rR582bY29vDwMAArq6uuHTpktQh1Xnh4eHo2bMnjI2NYWlpieHDhyMhIUFpm9zcXPj7+6NJkyYwMjLCe++9V2pBSCrbihUrIJPJEBgYqGjj+RTn8ePHmDBhApo0aYJGjRqhU6dOuHLliuJ5QRCwaNEi2NjYoFGjRvDw8MAff/whYcR1l1wux8KFC+Hg4IBGjRqhVatWWLp0KUrWEeL5VD/2TeKxb6pZ7Juqj32T+tSZvkloQPbt2yfo6ekJO3bsEH7//Xdh6tSpgpmZmZCamip1aHWap6enEBkZKdy6dUuIi4sThg4dKrRs2VLIyspSbOPn5yfY2toK0dHRwpUrV4TevXsL7u7uEkZdP1y6dEmwt7cXOnfuLMycOVPRzvOpuufPnwt2dnaCj4+PcPHiReHBgwfCyZMnhXv37im2WbFihWBqaiocOXJEuHHjhvDOO+8IDg4Owj///CNh5HXT8uXLhSZNmgjHjx8XEhMThYMHDwpGRkbChg0bFNvwfKoX+6aqYd9Uc9g3VR/7JvWqK31Tg0qcevXqJfj7+ysey+VyoXnz5kJ4eLiEUdU/aWlpAgDh7NmzgiAIQnp6uqCrqyscPHhQsc3t27cFAEJsbKxUYdZ5L168EJydnYVTp04J/fr1U3ROPJ/izJs3T3j99dfLfb6oqEiwtrYWVq1apWhLT08X9PX1hW+++aY2QqxXhg0bJkyePFmpbeTIkcL48eMFQeD5rAnsm9SDfZN6sG9SD/ZN6lVX+qYGM1UvPz8fV69ehYeHh6JNS0sLHh4eiI2NlTCy+icjIwMAYGFhAQC4evUqCgoKlM5t27Zt0bJlS57bCvj7+2PYsGFK5w3g+RTr+++/R48ePTBq1ChYWlrCxcUFW7duVTyfmJiIlJQUpfNpamoKV1dXns8yuLu7Izo6Gnfv3gUA3LhxA+fPn8eQIUMA8HyqG/sm9WHfpB7sm9SDfZN61ZW+SUdtR6rjnj59Crlcrlj5vZiVlRXu3LkjUVT1T1FREQIDA9GnTx907NgRAJCSkgI9PT2YmZkpbWtlZYWUlBQJoqz79u3bh2vXruHy5culnuP5FOfBgwfYsmULgoKC8O9//xuXL19GQEAA9PT04O3trThnZf3b5/ksbf78+cjMzETbtm2hra0NuVyO5cuXY/z48QDA86lm7JvUg32TerBvUh/2TepVV/qmBpM4kXr4+/vj1q1bOH/+vNSh1FuPHj3CzJkzcerUKRgYGEgdTr1XVFSEHj16ICwsDADg4uKCW7duISIiAt7e3hJHV/8cOHAAe/bswd69e9GhQwfExcUhMDAQzZs35/mkOot9U/Wxb1Iv9k3qVVf6pgYzVa9p06bQ1tYuVf0lNTUV1tbWEkVVv0yfPh3Hjx/HmTNn8Nprrynara2tkZ+fj/T0dKXteW7LdvXqVaSlpaFbt27Q0dGBjo4Ozp49i40bN0JHRwdWVlY8nyLY2Nigffv2Sm3t2rVDUlISACjOGf/tq2bu3LmYP38+PvjgA3Tq1AkTJ07ErFmzEB4eDoDnU93YN1Uf+yb1YN+kXuyb1Kuu9E0NJnHS09ND9+7dER0drWgrKipCdHQ03NzcJIys7hMEAdOnT8d3332H06dPw8HBQen57t27Q1dXV+ncJiQkICkpiee2DG+++SZu3ryJuLg4xa1Hjx4YP3684j7Pp+r69OlTqgTx3bt3YWdnBwBwcHCAtbW10vnMzMzExYsXeT7LkJOTAy0t5a5BW1sbRUVFAHg+1Y19U9Wxb1Iv9k3qxb5JvepM36S2MhP1wL59+wR9fX1h586dQnx8vPCvf/1LMDMzE1JSUqQOrU776KOPBFNTUyEmJkZITk5W3HJychTb+Pn5CS1bthROnz4tXLlyRXBzcxPc3NwkjLp+KVm5SBB4PsW4dOmSoKOjIyxfvlz4448/hD179giGhobC119/rdhmxYoVgpmZmXD06FHht99+E959912WfC2Ht7e30KJFC0XJ18OHDwtNmzYVPv74Y8U2PJ/qxb6patg31Tz2TVXHvkm96krf1KASJ0EQhM8++0xo2bKloKenJ/Tq1Uv49ddfpQ6pzgNQ5i0yMlKxzT///CNMmzZNMDc3FwwNDYURI0YIycnJ0gVdz7zaOfF8inPs2DGhY8eOgr6+vtC2bVvhyy+/VHq+qKhIWLhwoWBlZSXo6+sLb775ppCQkCBRtHVbZmamMHPmTKFly5aCgYGB4OjoKHzyySdCXl6eYhueT/Vj3yQe+6aax76petg3qU9d6ZtkglBiyV0iIiIiIiIqpcFc40RERERERFRVTJyIiIiIiIgqwcSJiIiIiIioEkyciIiIiIiIKsHEiYiIiIiIqBJMnIiIiIiIiCrBxImIiIiIiKgSTJyI6hgfHx8MHz5c6jCIiIgU2DcRATpSB0DUkMhksgqfDwkJwYYNG8B1qYmIqLawbyJSjUzgvwKiWpOSkqK4v3//fixatAgJCQmKNiMjIxgZGUkRGhERNVDsm4hUw6l6RLXI2tpacTM1NYVMJlNqMzIyKjUdon///pgxYwYCAwNhbm4OKysrbN26FdnZ2fD19YWxsTGcnJzw448/Kr3WrVu3MGTIEBgZGcHKygoTJ07E06dPa/kdExFRXce+iUg1TJyI6oFdu3ahadOmuHTpEmbMmIGPPvoIo0aNgru7O65du4ZBgwZh4sSJyMnJAQCkp6dj4MCBcHFxwZUrVxAVFYXU1FSMHj1a4ndCRESagn0TNTRMnIjqgS5dumDBggVwdnZGcHAwDAwM0LRpU0ydOhXOzs5YtGgRnj17ht9++w0AsGnTJri4uCAsLAxt27aFi4sLduzYgTNnzuDu3bsSvxsiItIE7JuooWFxCKJ6oHPnzor72traaNKkCTp16qRos7KyAgCkpaUBAG7cuIEzZ86UOSf9/v37aN26dQ1HTEREmo59EzU0TJyI6gFdXV2lxzKZTKmtuCJSUVERACArKwteXl5YuXJlqWPZ2NjUYKRERNRQsG+ihoaJE5EG6tatGw4dOgR7e3vo6PCfORERSY99E9V3vMaJSAP5+/vj+fPnGDt2LC5fvoz79+/j5MmT8PX1hVwulzo8IiJqgNg3UX3HxIlIAzVv3hwXLlyAXC7HoEGD0KlTJwQGBsLMzAxaWvxnT0REtY99E9V3XACXiIiIiIioEkzviYiIiIiIKsHEiYiIiIiIqBJMnIiIiIiIiCrBxImIiIiIiKgSTJyIiIiIiIgqwcSJiIiIiIioEkyciIiIiIiIKsHEiYiIiIiIqBJMnIiIiIiIiCrBxImIiIiIiKgSTJyIiIiIiIgqwcSJiIiIiIioEv8PK4xGcDGt0jkAAAAASUVORK5CYII=", "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", "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 }