{ "cells": [ { "cell_type": "markdown", "id": "cb5649fb-fbec-4281-aa8b-5f5c802d8e13", "metadata": {}, "source": [ "# Lifetime regressions" ] }, { "cell_type": "markdown", "id": "d27b45aa-cfc2-462f-ba23-5b6d2e0e7a4e", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:51.641448Z", "start_time": "2025-09-02T12:31:51.638613Z" } }, "source": [ "A built-in dataset of insulator string lifetimes and covariates." ] }, { "cell_type": "code", "execution_count": 1, "id": "613e4f41-4c7b-4998-92ba-c32b1266356c", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:52.179542Z", "start_time": "2025-09-02T12:31:52.173357Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "('time', 'event', 'entry', 'pHCl', 'pH2SO4', 'HNO3')\n" ] } ], "source": [ "from relife.datasets import load_insulator_string\n", "\n", "dataset = load_insulator_string()\n", "print(dataset.dtype.names)" ] }, { "cell_type": "code", "execution_count": 2, "id": "0e83fb9f-d750-4109-a1c6-d7a2639d47d3", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:52.567088Z", "start_time": "2025-09-02T12:31:52.563513Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[0.49 0.76 0.43 ... 1.12 1.19 0.35]\n" ] } ], "source": [ "print(dataset[\"pHCl\"])" ] }, { "cell_type": "markdown", "id": "211ec5be-cfc0-4897-9987-44b6b8fbf50b", "metadata": {}, "source": [ "Covariates must be grouped in one ``np.ndarray``." ] }, { "cell_type": "code", "execution_count": 3, "id": "65a1c5b3-dedc-4799-bc9d-301987024177", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:53.015236Z", "start_time": "2025-09-02T12:31:53.011333Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "(12000, 3)\n" ] } ], "source": [ "import numpy as np\n", "\n", "covar = np.column_stack((dataset[\"pHCl\"], dataset[\"pH2SO4\"], dataset[\"HNO3\"]))\n", "print(covar.shape)" ] }, { "cell_type": "code", "execution_count": 4, "id": "7eae8b87-3cf6-4850-90ab-f445bfbfc7a1", "metadata": {}, "outputs": [], "source": [ "from relife.lifetime_models import ParametricProportionalHazard, Gompertz" ] }, { "cell_type": "code", "execution_count": 5, "id": "2826f516-aa7a-4bad-9eb7-86e296de1f21", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:54.407376Z", "start_time": "2025-09-02T12:31:53.549002Z" } }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "regression = ParametricProportionalHazard(Gompertz())\n", "regression.fit(\n", " dataset[\"time\"],\n", " covar,\n", " event=dataset[\"event\"],\n", " entry=dataset[\"entry\"]\n", ")" ] }, { "cell_type": "code", "execution_count": 6, "id": "704e8ca0-4165-40a2-a473-560233bb1c57", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:55.538299Z", "start_time": "2025-09-02T12:31:55.535144Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 4.10887984 -2.66050611 3.23411769 0.21810962 0.02945802] ('coef_1', 'coef_2', 'coef_3', 'shape', 'rate')\n" ] } ], "source": [ "print(regression.get_params(), regression.get_params_names())" ] }, { "cell_type": "markdown", "id": "da36f246-9f28-4252-88a9-0a641ebd9f9b", "metadata": {}, "source": [ "Notice that the regression parameters include parameters from the baseline (Gompertz) and the covariate effect." ] }, { "cell_type": "code", "execution_count": 7, "id": "e4be562b-f824-478e-9ec9-3ed867801e3c", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:56.225005Z", "start_time": "2025-09-02T12:31:56.065229Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGdCAYAAAAMm0nCAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAkiBJREFUeJzs/Xl8VOX5+P+/zuyTTHaysoQtsi8B2bWisrmgVO1H+7Vvpa325/sNWqW0ldqi0FpsFXervuunLq1WP7YulVoqgoACgoJYQQgSCIGQkH0mM5PZz++PkwwJJJDlTBa4no/HPGbmzDn33CdGc3nf133diqqqKkIIIYQQvYihuzsghBBCCNFeEsAIIYQQoteRAEYIIYQQvY4EMEIIIYTodSSAEUIIIUSvIwGMEEIIIXodCWCEEEII0etIACOEEEKIXsfU3R1oi0gkwvHjx0lISEBRlO7ujhBCCCHaQFVV6urqyMnJwWDQd8ykVwQwx48fp3///t3dDSGEEEJ0wNGjR+nXr5+ubfaKACYhIQHQfgCJiYnd3BshhBBCtIXL5aJ///7Rv+N66hUBTOO0UWJiogQwQgghRC8Ti/QPSeIVQgghRK8jAYwQQggheh0JYIQQQgjR6/SKHJi2CIfDBIPB7u6GEF3CaDRiMpmkrIAQ4rx1TgQwbrebY8eOoapqd3dFiC4TFxdHdnY2Foulu7sihBBdrtcHMOFwmGPHjhEXF0d6err8H6k456mqSiAQoKKigsOHD5OXl6d7gSghhOjpen0AEwwGUVWV9PR07HZ7d3dHiC5ht9sxm80cOXKEQCCAzWbr7i4JIUSXOmf+t01GXsT5RkZdhBDnM/kvoBBCCCF6HQlghBBCCNHrSADTTWbOnMndd999xnPKysqYPXs28fHxJCcnd/i79u/fz9SpU7HZbIwfP77N123cuBFFUaitre3wd/cECxcuZMGCBd3dDSGEEDqSAKYHe+yxxygtLWX37t0cOHCgw+3cf//9xMfHU1BQwPr161s8py0BlV727t3L9ddfz8CBA1EUhccff/y0c1atWsWkSZNISEggIyODBQsWUFBQcMZ2i4qKUBSF3bt3x6bjXeRcuQ8hhIglCWB6sMLCQiZOnEheXh4ZGRmdaueiiy4iNzeXtLQ0HXvYMV6vl8GDB/PQQw+RlZXV4jmbNm1i0aJFfPrpp6xbt45gMMicOXPweDxd3FshhBA9UbsCmGeffZaxY8dGd4WeNm0a//rXv854zZtvvsnw4cOx2WyMGTOG999/v1MdPhtVVfEGQt3yaG8hvUgkws9+9jNSU1PJysrigQceiH42cOBA/v73v/PKK6+gKAoLFy5stY2VK1fSr18/rFYr48ePZ+3atdHPFUVh586drFy5EkVRmn1Ho4ULF7Jp0yaeeOIJFEVBURSKioqin+/cuZMLL7yQuLg4pk+fftpIyLvvvsuECROw2WwMHjyYFStWEAqFWr3vSZMm8fDDD3PTTTdhtVpbPGft2rUsXLiQUaNGMW7cOF566SWKi4vZuXNnq+0OGjQIgPz8fBRFYebMmc0+f+SRR8jOziYtLY1FixadsXJzYWEh1157LZmZmTgcDiZNmsSHH37Y7Jw//OEP5OXlYbPZyMzM5IYbboh+9re//Y0xY8Zgt9tJS0tj1qxZzYKvF154gREjRmCz2Rg+fDh/+MMf2nwfQggh2lkHpl+/fjz00EPk5eWhqiovv/wy1157LV988QWjRo067fytW7fy3e9+l1WrVnH11Vfz2muvsWDBAnbt2sXo0aN1u4mm6oNhRi7/d0zaPpuvV84lztL2H+nLL7/MkiVL2L59O9u2bWPhwoXMmDGD2bNn89lnn3HLLbeQmJjIE0880WqNmyeeeILVq1fz/PPPk5+fz5/+9CeuueYa9u7dS15eHqWlpcyaNYt58+axdOlSHA5Hi20cOHCA0aNHs3LlSgDS09OjQcx9993H6tWrSU9P54477uAHP/gBW7ZsAeDjjz/mlltu4cknn+Tiiy+msLCQH/3oR4A2daUXp9MJQGpqaqvn7Nixg8mTJ/Phhx8yatSoZhVqP/roI7Kzs/noo484ePAgN954I+PHj+f2229vsS23282VV17Jgw8+iNVq5ZVXXmH+/PkUFBQwYMAAPv/8c+666y7+/Oc/M336dKqrq/n4448BKC0t5bvf/S6///3v+fa3v01dXR0ff/xxNMB99dVXWb58OU8//TT5+fl88cUX3H777cTHx3Prrbee8T6EEKI3CYYjMWtbUTtZfz81NZWHH36YH/7wh6d9duONN+LxeFizZk302NSpUxk/fjzPPfdcm7/D5XKRlJSE0+kkMTGx2Wc+n4/Dhw8zaNAgbDYb3kCoVwQwM2fOJBwOR//oAUyePJnLLruMhx56CIAFCxaQnJzMSy+91Go7ffv2ZdGiRfziF79o1s6kSZN45plnABg/fjwLFixocfSlaX/Gjx/fLB9l48aNXHrppXz44YdcfvnlALz//vtcddVV1NfXY7PZmDVrFpdffjnLli2LXveXv/yFn/3sZxw/fvysP4eBAwdy9913nzH/JhKJcM0111BbW8snn3zS6nlFRUUMGjSIL774olmy8sKFC9m4cSOFhYUYjUYA/s//+T8YDAZef/31s/ax0ejRo7njjjtYvHgxb731Ft///vc5duwYCQkJzc7btWsXEydOpKioiNzc3NPaGTp0KL/+9a/57ne/Gz32m9/8hvfff5+tW7e2eh+nOvV3Xwghepqi0koG5aS3+Pe7szpciTccDvPmm2/i8XiYNm1ai+ds27aNJUuWNDs2d+5c3nnnnTO27ff78fv90fcul6vN/bKbjXy9cm6bz9eT3Wxs1/ljx45t9j47O5vy8vI2X+9yuTh+/DgzZsxodnzGjBl8+eWX7erLmTTtZ3Z2NgDl5eUMGDCAL7/8ki1btvDggw9GzwmHw/h8PrxeL3FxcZ3+/kWLFrFnz55mwcsdd9zBX/7yl+h7t9t9xjZGjRoVDV4a7+Orr75q9Xy3280DDzzAP//5T0pLSwmFQtTX11NcXAzA7Nmzyc3NZfDgwcybN4958+bx7W9/m7i4OMaNG8fll1/OmDFjmDt3LnPmzOGGG24gJSUFj8dDYWEhP/zhD5uN/oRCIZKSktr9sxFCiJ7M7Y/dJsvtDmC++uorpk2bhs/nw+Fw8PbbbzNy5MgWzy0rKyMzM7PZsczMTMrKys74HatWrWLFihXt7Rqg5Xy0ZxqnO5nN5mbvFUUhEondcFtHNe1nY8Xjxn663W5WrFjBddddd9p1eowKLF68mDVr1rB582b69esXPb5y5UqWLl3a5nba+7NeunQp69at45FHHmHo0KHY7XZuuOEGAoEAAAkJCezatYuNGzfywQcfsHz5ch544AE+++wzkpOTWbduHVu3buWDDz7gqaee4r777mP79u3RgO6Pf/wjU6ZMafadTQMsIYTo7YLhCPWBcMzab/cqpGHDhrF79262b9/Of//3f3Prrbfy9ddf69qpZcuW4XQ6o4+jR4/q2v65IjExkZycnGg+SqMtW7a0GlS2xmKxEA63/xdtwoQJFBQUMHTo0NMenSl1r6oqixcv5u2332bDhg3RxNZGGRkZzb6r8R6ADt3HqbZs2cLChQv59re/zZgxY8jKymqW2AxgMpmYNWsWv//97/nPf/5DUVERGzZsALQAacaMGaxYsYIvvvgCi8XC22+/TWZmJjk5ORw6dOi0n1fjPep5H0II0V3qfK0v5tBDu4cqLBZL9A/GxIkT+eyzz3jiiSd4/vnnTzs3KyuLEydONDt24sSJVpfONrJara2uThHN/fSnP+X+++9nyJAhjB8/nhdffJHdu3fz6quvtqudgQMHsn37doqKinA4HGdMlm1q+fLlXH311QwYMIAbbrgBg8HAl19+yZ49e/jNb37T4jWBQCAa9AYCAUpKSti9ezcOhyP6u7Vo0SJee+013n33XRISEqKjdklJSa0mNGdkZGC321m7di39+vXDZrN1eFomLy+Pt956i/nz56MoCr/61a+ajdisWbOGQ4cO8a1vfYuUlBTef/99IpEIw4YNY/v27axfv545c+aQkZHB9u3bqaioYMSIEQCsWLGCu+66i6SkJObNm4ff7+fzzz+npqaGJUuW6HofQgjRXVz1sZs+Ah3qwEQikWb5Kk1NmzbttMJp69atazVnRrTfXXfdxZIlS/jJT37CmDFjWLt2Lf/4xz/Iy8trVztLly7FaDQycuRI0tPTo7keZzN37lzWrFnDBx98wKRJk5g6dSqPPfZYi8mrjY4fP05+fj75+fmUlpbyyCOPkJ+fz2233RY959lnn8XpdDJz5kyys7OjjzfeeKPVdk0mE08++STPP/88OTk5XHvttW3/AZzi0UcfJSUlhenTpzN//nzmzp3LhAkTop8nJyfz1ltvcdlllzFixAiee+45/vrXvzJq1CgSExPZvHkzV155JRdccAG//OUvWb16NVdccQUAt912Gy+88AIvvvgiY8aM4ZJLLuGll16KjsDoeR9CCNFdYj0C065VSMuWLeOKK65gwIAB1NXV8dprr/G73/2Of//738yePZtbbrmFvn37smrVKkBbRn3JJZfw0EMPcdVVV/H666/z29/+tt3LqNuzCkmI84X87gsheqpQOMLnR2owhuqZPKx/969CKi8v55ZbbqG0tJSkpCTGjh0bDV4AiouLm+U9TJ8+nddee41f/vKX/OIXvyAvL4933nknZjVghBBCCNH96nwhOlek5ew6XQemK8gIjBCnk999IURPVVzlpaS2HlM4diMwsheSEEIIIXTl8sU2gRckgBFCCCGEjsIRFbc/tgm8IAGMEEIIIXTk7oL8F5AARgghhBA66orpI5AARgghhBA6kgBGCCGEEL1KJKLijnEBu0YSwAghhBBCF+5AiEgXFWfpHds2d8C2wqou/b5pQ9Ladf7ChQupra3lnXfeAbSdux988EH++c9/UlJSQkZGBuPHj+fuu+/m8ssvb7Wd6upq7rzzTt577z0MBgPXX389TzzxBA6Ho9VrCgsLWbp0KZ988gl+v5958+bx1FNPNds5fODAgRw5cqTZdatWreLee+9t130KIYQ4f8R6/6OmZASmBygqKmLixIls2LCBhx9+mK+++oq1a9dy6aWXsmjRojNee/PNN7N3717WrVvHmjVr2Lx5Mz/60Y9aPd/j8TBnzhwURWHDhg1s2bKFQCDA/Pnzm21WCLBy5UpKS0ujjzvvvFOX+xVCCHFuivX+R02dsyMwvcn//M//oCgKO3bsID4+Pnp81KhR/OAHP2j1un379rF27Vo+++wzLrzwQgCeeuoprrzySh555BFycnJOu2bLli0UFRXxxRdfRKsivvzyy6SkpLBhwwZmzZoVPTchIeGsO4cLIYQQAKqqdmkAIyMw3ay6upq1a9eyaNGiZsFLo+Tk5Fav3bZtG8nJydHgBWDWrFkYDAa2b9/e4jV+vx9FUbBardFjNpsNg8HAJ5980uzchx56iLS0NPLz83n44YcJhbruF1MIIUTv4gmECXdVAgwyAtPtDh48iKqqDB8+vN3XlpWVkZGR0eyYyWQiNTWVsrKyFq+ZOnUq8fHx/PznP+e3v/0tqqpy7733Eg6HKS0tjZ531113MWHCBFJTU9m6dSvLli2jtLSURx99tN39FEIIce7ryvwXkBGYbtfWvTTvuOMOHA5H9NFR6enpvPnmm7z33ns4HA6SkpKora1lwoQJzXYSX7JkCTNnzmTs2LHccccdrF69mqeeegq/39/h7xZCCHHu6qr6L41kBKab5eXloSgK+/fvP+N5K1euZOnSpc2OZWVlUV5e3uxYKBSiurr6jLkrc+bMobCwkMrKSkwmE8nJyWRlZTF48OBWr5kyZQqhUIiioiKGDRvWhjsTQghxPunK/BeQEZhul5qayty5c3nmmWfweDynfV5bWwtARkYGQ4cOjT4Apk2bRm1tLTt37oyev2HDBiKRCFOmTDnrd/fp04fk5GQ2bNhAeXk511xzTavn7t69G4PBcNqUlRBCCOHxhwiFuy7/BSSA6RGeeeYZwuEwkydP5u9//zvffPMN+/bt48knn2TatGmtXjdixAjmzZvH7bffzo4dO9iyZQuLFy/mpptuiq5AKikpYfjw4ezYsSN63Ysvvsinn35KYWEhf/nLX/jOd77DPffcEx1Z2bZtG48//jhffvklhw4d4tVXX+Wee+7he9/7HikpKbH9YQghhOh1unr0BWQKqUcYPHgwu3bt4sEHH+QnP/kJpaWlpKenM3HiRJ599tkzXvvqq6+yePFiLr/88mghuyeffDL6eTAYpKCgAK/XGz1WUFDAsmXLqK6uZuDAgdx3333cc8890c+tViuvv/46DzzwAH6/n0GDBnHPPfewZMkS/W9eCCFEr9fV+S8AitrWLNJu5HK5SEpKwul0RmuXNPL5fBw+fJhBgwZhs9m6qYdCdD353RdC9BQ7j1QTCJ0eTpjC9Uwe1r/Fv9+dJVNIQgghhOiw+kC4xeAl1iSAEUIIIUSHdcf0EUgAI4QQQohOqJMARgghhBC9jbO+e7aZkQBGCCGEEB3iC4YJhCLd8t0SwAghhBCiQ7or/wUkgBFCCCFEB3VHAbtGEsAIIYQQokO6egfqpiSAEUIIIUS7+UNhfMHuyX8BCWC6zcyZM7n77rvPeE5ZWRmzZ88mPj6e5OTkDn/X/v37mTp1KjabjfHjx7f5uo0bN6IoSnRDyd5q4cKFLFiwoLu7IYQQ55TunD4CCWB6tMcee4zS0lJ2797NgQMHOtzO/fffT3x8PAUFBaxfv77Fc9oSUOll7969XH/99QwcOBBFUXj88cdPO2fVqlVMmjSJhIQEMjIyWLBgAQUFBWdst6ioCEVR2L17d2w63kXOlfsQQpzbunP6CCSA6dEKCwuZOHEieXl5ZGRkdKqdiy66iNzcXNLS0nTsYcd4vV4GDx7MQw89RFZWVovnbNq0iUWLFvHpp5+ybt06gsEgc+bMwePxdHFvhRBCtERGYM5jkUiEn/3sZ6SmppKVlcUDDzwQ/WzgwIH8/e9/55VXXkFRFBYuXNhqGytXrqRfv35YrVbGjx/P2rVro58risLOnTtZuXIliqI0+45GCxcuZNOmTTzxxBMoioKiKBQVFUU/37lzJxdeeCFxcXFMnz79tJGQd999lwkTJmCz2Rg8eDArVqwgFGr9F3vSpEk8/PDD3HTTTVit1hbPWbt2LQsXLmTUqFGMGzeOl156ieLiYnbu3Nlqu4MGDQIgPz8fRVGYOXNms88feeQRsrOzSUtLY9GiRQSDrf/fQ2FhIddeey2ZmZk4HA4mTZrEhx9+2OycP/zhD+Tl5WGz2cjMzOSGG26Ifva3v/2NMWPGYLfbSUtLY9asWc2CrxdeeIERI0Zgs9kYPnw4f/jDH9p8H0II0d2C4QjeQLhb+2Dq1m+PBVWFoLd7vtscB4rS5tNffvlllixZwvbt29m2bRsLFy5kxowZzJ49m88++4xbbrmFxMREnnjiCex2e4ttPPHEE6xevZrnn3+e/Px8/vSnP3HNNdewd+9e8vLyKC0tZdasWcybN4+lS5ficDhabOPAgQOMHj2alStXApCenh4NYu677z5Wr15Neno6d9xxBz/4wQ/YsmULAB9//DG33HILTz75JBdffDGFhYX86Ec/ArSpK704nU4AUlNTWz1nx44dTJ48mQ8//JBRo0ZhsViin3300UdkZ2fz0UcfcfDgQW688UbGjx/P7bff3mJbbrebK6+8kgcffBCr1corr7zC/PnzKSgoYMCAAXz++efcdddd/PnPf2b69OlUV1fz8ccfA1BaWsp3v/tdfv/73/Ptb3+buro6Pv74Yxo3fn/11VdZvnw5Tz/9NPn5+XzxxRfcfvvtxMfHc+utt57xPoQQoifo7tEXOBcDmKAXfpvTPd/9i+NgiW/z6WPHjo3+kc/Ly+Ppp59m/fr1zJ49m/T0dKxWK3a7vdVpFtBGFX7+859z0003AfC73/2Ojz76iMcff5xnnnmGrKwsTCYTDoej1XaSkpKwWCzExcW1eM6DDz7IJZdcAsC9997LVVddhc/nw2azsWLFCu69915uvfVWAAYPHsyvf/1rfvazn+kWwEQiEe6++25mzJjB6NGjWz0vPT0dgLS0tNPuIyUlhaeffhqj0cjw4cO56qqrWL9+fasBzLhx4xg3blz0/a9//Wvefvtt/vGPf7B48WKKi4uJj4/n6quvJiEhgdzcXPLz8wEtgAmFQlx33XXk5uYCMGbMmGhb999/P6tXr+a6664DtBGXr7/+mueff55bb731jPchhBA9gbOb819AppC61dixY5u9z87Opry8vM3Xu1wujh8/zowZM5odnzFjBvv27dOlj9C8n9nZ2QDRfn755ZesXLkSh8MRfdx+++2Ulpbi9eozErZo0SL27NnD66+/Hj12xx13NPvOsxk1ahRGo7HZfZzpZ+12u1m6dCkjRowgOTkZh8PBvn37KC4uBmD27Nnk5uYyePBg/uu//otXX301er/jxo3j8ssvZ8yYMXznO9/hj3/8IzU1NQB4PB4KCwv54Q9/2Kz/v/nNbygsLOzQz0cIIbparTfQ3V04B0dgzHHaSEh3fXd7Tjebm71XFIVIpPvW1LemaT+Vhimyxn663W5WrFgRHU1oymazdfq7Fy9ezJo1a9i8eTP9+vWLHl+5ciVLly5tczvt/VkvXbqUdevW8cgjjzB06FDsdjs33HADgYD2L21CQgK7du1i48aNfPDBByxfvpwHHniAzz77jOTkZNatW8fWrVv54IMPeOqpp7jvvvvYvn07cXHa78gf//hHpkyZ0uw7mwZYQgjRU/mC3Vv/pdG5F8AoSrumcXqzxMREcnJy2LJlS3SKB2DLli1Mnjy5XW1ZLBbC4fYnZE2YMIGCggKGDh3a7mvPRFVV7rzzTt5++202btwYTWxtlJGRcdrKrMZckY7cx6m2bNnCwoUL+fa3vw1ogVrTxGYAk8nErFmzmDVrFvfffz/Jycls2LCB6667DkVRmDFjBjNmzGD58uXk5uby9ttvs2TJEnJycjh06BA333xzi9+t530IIYTeesL0EZyLAcx55qc//Sn3338/Q4YMYfz48bz44ovs3r2bV199tV3tDBw4kO3bt1NUVITD4ThjsmxTy5cv5+qrr2bAgAHccMMNGAwGvvzyS/bs2cNvfvObFq8JBAJ8/fXX0dclJSXs3r0bh8MRDYQWLVrEa6+9xrvvvktCQgJlZWWAlq/TWkJzRkYGdrudtWvX0q9fP2w2G0lJSe36OTTKy8vjrbfeYv78+SiKwq9+9atmIzZr1qzh0KFDfOtb3yIlJYX333+fSCTCsGHD2L59O+vXr2fOnDlkZGSwfft2KioqGDFiBAArVqzgrrvuIikpiXnz5uH3+/n888+pqalhyZIlut6HEELordbbMwIYyYHp5e666y6WLFnCT37yE8aMGcPatWv5xz/+QV5eXrvaWbp0KUajkZEjR5Kenh7N9TibuXPnsmbNGj744AMmTZrE1KlTeeyxx6LJqy05fvw4+fn55OfnU1payiOPPEJ+fj633XZb9Jxnn30Wp9PJzJkzyc7Ojj7eeOONVts1mUw8+eSTPP/88+Tk5HDttde2/QdwikcffZSUlBSmT5/O/PnzmTt3LhMmTIh+npyczFtvvcVll13GiBEjeO655/jrX//KqFGjSExMZPPmzVx55ZVccMEF/PKXv2T16tVcccUVANx222288MILvPjii4wZM4ZLLrmEl156KTrKpOd9CCGEnlRV7dYdqJtS1Ma1nT2Yy+UiKSkJp9NJYmJis898Ph+HDx9m0KBBuuRcCNFbyO++EKKruXxB9pa42ny+KVzP5GH9W/z73VkyAiOEEEKINnH2kOkjkABGCCGEEG3UU/JfQAIYIYQQQrRBMBzBE+j+CryNJIARQgghxFk564P0pKxZCWCEEEIIcVY9afoIJIARQgghRBv0lAJ2jSSAEUIIIcQZeQMhAqHu3z6gKQlghBBCCHFGPW36CNoZwKxatYpJkyaRkJBARkYGCxYsoKCg4IzXvPTSSyiK0uwhRbeEEEKI3qOnTR9BOwOYTZs2sWjRIj799FPWrVtHMBhkzpw5eDyeM16XmJhIaWlp9HHkyJFOdVoIIYQQXSMSUXH1wACmXZs5rl27ttn7l156iYyMDHbu3Mm3vvWtVq9TFIWsrKyO9bCjDn/ctd836OJ2nb5w4UJqa2t55513ACgrK+PBBx/kn//8JyUlJWRkZDB+/HjuvvtuLr/88lbbqa6u5s477+S9997DYDBw/fXX88QTT+BwOFq9prCwkKVLl/LJJ5/g9/uZN28eTz31FJmZmdFzBg4ceFqguWrVKu6999523acQQojezeULEulBy6cbdSoHxul0Apx152K3201ubi79+/fn2muvZe/evWc83+/343K5mj3OZUVFRUycOJENGzbw8MMP89VXX7F27VouvfRSFi1adMZrb775Zvbu3cu6detYs2YNmzdv5kc/+lGr53s8HubMmYOiKGzYsIEtW7YQCASYP39+s92WAVauXNls5OzOO+/U5X6FEEL0Hj1x+gjaOQLTVCQS4e6772bGjBmMHj261fOGDRvGn/70J8aOHYvT6eSRRx5h+vTp7N27l379+rV4zapVq1ixYkVHu9br/M///A+KorBjxw7i4+Ojx0eNGsUPfvCDVq/bt28fa9eu5bPPPuPCCy8E4KmnnuLKK6/kkUceIScn57RrtmzZQlFREV988UV0Y62XX36ZlJQUNmzYwKxZs6LnJiQkdP3ImRBCiB6lJybwQidGYBYtWsSePXt4/fXXz3jetGnTuOWWWxg/fjyXXHIJb731Funp6Tz//POtXrNs2TKcTmf0cfTo0Y52s8errq5m7dq1LFq0qFnw0ig5ObnVa7dt20ZycnI0eAGYNWsWBoOB7du3t3iN3+9HURSsVmv0mM1mw2Aw8MknnzQ796GHHiItLY38/HwefvhhQqGeU0JaCCFE7PlDYbyBcHd3o0UdGoFZvHhxdLqitVGU1pjNZvLz8zl48GCr51it1mZ/YM9lBw8eRFVVhg8f3u5ry8rKyMjIaHbMZDKRmppKWVlZi9dMnTqV+Ph4fv7zn/Pb3/4WVVW59957CYfDlJaWRs+76667mDBhAqmpqWzdupVly5ZRWlrKo48+2u5+CiGE6J166vQRtHMERlVVFi9ezNtvv82GDRsYNGhQu78wHA7z1VdfkZ2d3e5rz0VqGzeWuOOOO3A4HNFHR6Wnp/Pmm2/y3nvv4XA4SEpKora2lgkTJmAwnPx1WLJkCTNnzmTs2LHccccdrF69mqeeegq/39/h7xZCCNG7OHvo9BG0cwRm0aJFvPbaa7z77rskJCRE/y8/KSkJu90OwC233ELfvn1ZtWoVoCWCTp06laFDh1JbW8vDDz/MkSNHuO2223S+ld4pLy8PRVHYv3//Gc9buXIlS5cubXYsKyuL8vLyZsdCoRDV1dVnzF2ZM2cOhYWFVFZWYjKZSE5OJisri8GDB7d6zZQpUwiFQhQVFTFs2LA23JkQQojeTFXVc2cE5tlnn8XpdDJz5kyys7OjjzfeeCN6TnFxcbOpiJqaGm6//XZGjBjBlVdeicvlYuvWrYwcOVK/u+jFUlNTmTt3Ls8880yL9XRqa2sByMjIYOjQodEHaPlFtbW17Ny5M3r+hg0biEQiTJky5azf3adPH5KTk9mwYQPl5eVcc801rZ67e/duDAbDaVNWQgghzk2eQJhguAeun27QrhGYtkx3bNy4sdn7xx57jMcee6xdnTrfPPPMM8yYMYPJkyezcuVKxo4dSygUYt26dTz77LPs27evxetGjBjBvHnzuP3223nuuecIBoMsXryYm266KboCqaSkhMsvv5xXXnmFyZMnA/Diiy8yYsQI0tPT2bZtGz/+8Y+55557oiMr27ZtY/v27Vx66aUkJCSwbds27rnnHr73ve+RkpLSNT8UIYQQ3arWG+juLpxRh5dRC/0MHjyYXbt28eCDD/KTn/yE0tJS0tPTmThxIs8+++wZr3311VdZvHgxl19+ebSQ3ZNPPhn9PBgMUlBQgNfrjR4rKChg2bJlVFdXM3DgQO677z7uueee6OdWq5XXX3+dBx54AL/fz6BBg7jnnntYsmSJ/jcvhBCiR+qpy6cbKWpbs0i7kcvlIikpCafTGa1d0sjn83H48GEGDRokeyyJ84r87gshYiUcUfmsqJrORgimcD2Th/Vv8e93Z8lu1EIIIYRoxlkf7HTwEmsSwAghhBCimZ6e/wISwAghhBDiFD15+XQjCWCEEEIIEeULhvEFI2c/sZtJACOEEEKIqJ6++qiRBDBCCCGEiOoN00cgAYwQQgghGvT07QOakgBGCCGEEAC4fCHCkR6+frqBBDBCCCGEAMDVS0ZfQAKYbjNz5kzuvvvuM55TVlbG7NmziY+PJzk5ucPftX//fqZOnYrNZmP8+PFtvm7jxo0oihLdULK3WrhwIQsWLOjubgghRI9X7en59V8aSQDTgz322GOUlpaye/duDhw40OF27r//fuLj4ykoKGD9+vUtntOWgEove/fu5frrr2fgwIEoisLjjz9+2jmrVq1i0qRJJCQkkJGRwYIFCygoKDhju0VFRSiKwu7du2PT8S5yrtyHEKJ38QXDeAPh7u5Gm0kA04MVFhYyceJE8vLyyMjI6FQ7F110Ebm5uaSlpenYw47xer0MHjyYhx56iKysrBbP2bRpE4sWLeLTTz9l3bp1BINB5syZg8fj6eLeCiHE+UHv0ZeIqnKoInb/zZYAphtFIhF+9rOfkZqaSlZWFg888ED0s4EDB/L3v/+dV155BUVRWLhwYattrFy5kn79+mG1Whk/fjxr166Nfq4oCjt37mTlypUoitLsOxotXLiQTZs28cQTT6AoCoqiUFRUFP18586dXHjhhcTFxTF9+vTTRkLeffddJkyYgM1mY/DgwaxYsYJQKNTqfU+aNImHH36Ym266CavV2uI5a9euZeHChYwaNYpx48bx0ksvUVxczM6dO1ttd9CgQQDk5+ejKAozZ85s9vkjjzxCdnY2aWlpLFq0iGCw9bnewsJCrr32WjIzM3E4HEyaNIkPP/yw2Tl/+MMfyMvLw2azkZmZyQ033BD97G9/+xtjxozBbreTlpbGrFmzmgVfL7zwAiNGjMBmszF8+HD+8Ic/tPk+hBAiFvQOYA6Wu/np3/6ja5tNmWLWcjdRVZX6UH23fLfdZEdRlDaf//LLL7NkyRK2b9/Otm3bWLhwITNmzGD27Nl89tln3HLLLSQmJvLEE09gt9tbbOOJJ55g9erVPP/88+Tn5/OnP/2Ja665hr1795KXl0dpaSmzZs1i3rx5LF26FIfD0WIbBw4cYPTo0axcuRKA9PT0aBBz3333sXr1atLT07njjjv4wQ9+wJYtWwD4+OOPueWWW3jyySe5+OKLKSws5Ec/+hGgTV3pxel0ApCamtrqOTt27GDy5Ml8+OGHjBo1CovFEv3so48+Ijs7m48++oiDBw9y4403Mn78eG6//fYW23K73Vx55ZU8+OCDWK1WXnnlFebPn09BQQEDBgzg888/56677uLPf/4z06dPp7q6mo8//hiA0tJSvvvd7/L73/+eb3/729TV1fHxxx/TuPH7q6++yvLly3n66afJz8/niy++4Pbbbyc+Pp5bb731jPchhBCxEAhFcPtb/x/Pjth+qErX9k51zgUw9aF6prw2pVu+e/v/t504c1ybzx87dmz0j3xeXh5PP/0069evZ/bs2aSnp2O1WrHb7a1Os4A2qvDzn/+cm266CYDf/e53fPTRRzz++OM888wzZGVlYTKZcDgcrbaTlJSExWIhLi6uxXMefPBBLrnkEgDuvfderrrqKnw+HzabjRUrVnDvvfdy6623AjB48GB+/etf87Of/Uy3ACYSiXD33XczY8YMRo8e3ep56enpAKSlpZ12HykpKTz99NMYjUaGDx/OVVddxfr161sNYMaNG8e4ceOi73/961/z9ttv849//IPFixdTXFxMfHw8V199NQkJCeTm5pKfnw9oAUwoFOK6664jNzcXgDFjxkTbuv/++1m9ejXXXXcdoI24fP311zz//PPceuutZ7wPIYSIhVpvQNfdpyOqyvbD1fo12AKZQupGY8eObfY+Ozub8vLyNl/vcrk4fvw4M2bMaHZ8xowZ7Nu3T5c+QvN+ZmdnA0T7+eWXX7Jy5UocDkf0cfvtt1NaWorX69Xl+xctWsSePXt4/fXXo8fuuOOOZt95NqNGjcJoNDa7jzP9rN1uN0uXLmXEiBEkJyfjcDjYt28fxcXFAMyePZvc3FwGDx7Mf/3Xf/Hqq69G73fcuHFcfvnljBkzhu985zv88Y9/pKamBgCPx0NhYSE//OEPm/X/N7/5DYWFhR36+QghRGdV67z7dGG5mypPAJs5dmHGOTcCYzfZ2f7/be+2724Ps9nc7L2iKEQiPW8Drab9bJwia+yn2+1mxYoV0dGEpmw2W6e/e/HixaxZs4bNmzfTr1+/6PGVK1eydOnSNrfT3p/10qVLWbduHY888ghDhw7Fbrdzww03EAho/5InJCSwa9cuNm7cyAcffMDy5ct54IEH+Oyzz0hOTmbdunVs3bqVDz74gKeeeor77ruP7du3ExenjdD98Y9/ZMqU5iOFTQMsIYToKuGIilPn/Y8aR18uzE3hG11bPumcC2AURWnXNE5vlpiYSE5ODlu2bIlO8QBs2bKFyZMnt6sti8VCONz+5XMTJkygoKCAoUOHtvvaM1FVlTvvvJO3336bjRs3RhNbG2VkZJy2MqsxV6Qj93GqLVu2sHDhQr797W8DWqDWNLEZwGQyMWvWLGbNmsX9999PcnIyGzZs4LrrrkNRFGbMmMGMGTNYvnw5ubm5vP322yxZsoScnBwOHTrEzTff3OJ363kfQghxNrXeAHoW31VVle2HtfyXaUPS+Kt+TTdzzgUw55uf/vSn3H///QwZMoTx48fz4osvsnv3bl599dV2tTNw4EC2b99OUVERDofjjMmyTS1fvpyrr76aAQMGcMMNN2AwGPjyyy/Zs2cPv/nNb1q8JhAI8PXXX0dfl5SUsHv3bhwORzQQWrRoEa+99hrvvvsuCQkJlJWVAVq+TmsJzRkZGdjtdtauXUu/fv2w2WwkJSW16+fQKC8vj7feeov58+ejKAq/+tWvmo3YrFmzhkOHDvGtb32LlJQU3n//fSKRCMOGDWP79u2sX7+eOXPmkJGRwfbt26moqGDEiBEArFixgrvuuoukpCTmzZuH3+/n888/p6amhiVLluh6H0IIcTY1ek8fVXiodGvTR/kDknVtuynJgenl7rrrLpYsWcJPfvITxowZw9q1a/nHP/5BXl5eu9pZunQpRqORkSNHkp6eHs31OJu5c+eyZs0aPvjgAyZNmsTUqVN57LHHosmrLTl+/Dj5+fnk5+dTWlrKI488Qn5+Prfddlv0nGeffRan08nMmTPJzs6OPt54441W2zWZTDz55JM8//zz5OTkcO2117b9B3CKRx99lJSUFKZPn878+fOZO3cuEyZMiH6enJzMW2+9xWWXXcaIESN47rnn+Otf/8qoUaNITExk8+bNXHnllVxwwQX88pe/ZPXq1VxxxRUA3Hbbbbzwwgu8+OKLjBkzhksuuYSXXnopOsqk530IIcSZqKpKje7TR9roS/6AFKym2E2NK6qqZ95xbLhcLpKSknA6nSQmJjb7zOfzcfjwYQYNGqRLzoUQvYX87gshOsvpDfJ1qUu39lRV5a7Xv6DSHeCeWRcwfYCdycP6t/j3u7NkBEYIIYQ4T+m9+uhQpTZ9ZDUZGN8/Wde2TyUBjBBCCHGe0rv67qcNxesmDEjBYoptiCEBjBBCCHEecvtDBEL6le5QVZXth7Tl01MGtW0hSGdIACOEEEKch6rd+k8fVbj92vRRDFcfNZIARgghhDgP6Z3/0rj30fj+yTFdfdRIAhghhBDiPFMfCFMf0K9Yptpk76Opg9N0a/dMJIARQgghzjN6j74crvRQXufHYoz96qNGEsAIIYQQ55kanVcfNY6+5A9Ixmbumn3dJIARQgghziOBUAS3P6Rbe6qqRpdPTxnUNdNHIAGMEEIIcV6p8QbQswZ/UZU3On0Uy72PTnXObub4WdlnXfp9k7Imtev8hQsXUltbyzvvvANAWVkZDz74IP/85z8pKSkhIyOD8ePHc/fdd3P55Ze32k51dTV33nkn7733HgaDgeuvv54nnngCh8PR6jWFhYUsXbqUTz75BL/fz7x583jqqafIzMyMnjNw4ECOHDnS7LpVq1Zx7733tus+hRBC9CyxKl43vgunj0BGYHqEoqIiJk6cyIYNG3j44Yf56quvWLt2LZdeeimLFi0647U333wze/fuZd26daxZs4bNmzfzox/9qNXzPR4Pc+bMQVEUNmzYwJYtWwgEAsyfP7/ZbssAK1eupLS0NPq48847dblfIYQQ3SMcUXHV67d5o7b6SAtgpnZB8bqmztkRmN7kf/7nf1AUhR07dhAfHx89PmrUKH7wgx+0et2+fftYu3Ytn332GRdeeCEATz31FFdeeSWPPPIIOTk5p12zZcsWioqK+OKLL6Iba7388sukpKSwYcMGZs2aFT03ISGBrKwsvW5TCCFEN6vxBojoPH10wuXHbFTIH5CiX8NtICMw3ay6upq1a9eyaNGiZsFLo+Tk5Fav3bZtG8nJydHgBWDWrFkYDAa2b9/e4jV+vx9FUbBardFjNpsNg8HAJ5980uzchx56iLS0NPLz83n44YcJhfRL+hJCCNH19F99pI2+5PdP6dLpI5ARmG538OBBVFVl+PDh7b62rKyMjIyMZsdMJhOpqamUlZW1eM3UqVOJj4/n5z//Ob/97W9RVZV7772XcDhMaWlp9Ly77rqLCRMmkJqaytatW1m2bBmlpaU8+uij7e6nEEKI7qeqKrV6Tx817n00uGunj0BGYLqd2sZU8DvuuAOHwxF9dFR6ejpvvvkm7733Hg6Hg6SkJGpra5kwYQIGw8lfhyVLljBz5kzGjh3LHXfcwerVq3nqqafw+/0d/m4hhBDdx1kfJBTWb/7oSLWXMpcPs1FhQhdPH4GMwHS7vLw8FEVh//79Zzxv5cqVLF26tNmxrKwsysvLmx0LhUJUV1efMXdlzpw5FBYWUllZiclkIjk5maysLAYPHtzqNVOmTCEUClFUVMSwYcPacGdCCCF6Er1XHzWOvozv37WrjxrJCEw3S01NZe7cuTzzzDN4PJ7TPq+trQUgIyODoUOHRh8A06ZNo7a2lp07d0bP37BhA5FIhClTppz1u/v06UNycjIbNmygvLyca665ptVzd+/ejcFgOG3KSgghRO9Q443N6qOuLF7XlAQwPcAzzzxDOBxm8uTJ/P3vf+ebb75h3759PPnkk0ybNq3V60aMGMG8efO4/fbb2bFjB1u2bGHx4sXcdNNN0RVIJSUlDB8+nB07dkSve/HFF/n0008pLCzkL3/5C9/5zne45557oiMr27Zt4/HHH+fLL7/k0KFDvPrqq9xzzz1873vfIyWl64cJhRBCdE6dL0ggFDn7iW1UXO2l1Nl900cgU0g9wuDBg9m1axcPPvggP/nJTygtLSU9PZ2JEyfy7LPPnvHaV199lcWLF3P55ZdHC9k9+eST0c+DwSAFBQV4vd7osYKCApYtW0Z1dTUDBw7kvvvu45577ol+brVaef3113nggQfw+/0MGjSIe+65hyVLluh/80IIIWKuxqPf6Auc3PtoXL9k7Jaunz4CUNS2ZpF2I5fLRVJSEk6nM1q7pJHP5+Pw4cMMGjQIm83WTT0UouvJ774Qoq2+KK7BF9RnBEZVVZa++SXHnT4WXTqUi4b2afVcU7ieycP6t/j3u7NkCkkIIYQ4h9X5groFLwBHa+o5Hp0+Stat3faSAEYIIYQ4h1W69V19tLWwEoCx/ZKJs3RfJooEMEIIIcQ5SlVVqj361e+KRFQ+/kYLYC4+w9RRV2hXALNq1SomTZpEQkICGRkZLFiwgIKCgrNe9+abbzJ8+HBsNhtjxozh/fff73CHhRBCCNE2zvoggZB+qa57jjup9gSItxqZkNu9q1LbFcBs2rSJRYsW8emnn7Ju3TqCwSBz5sxpsX5Jo61bt/Ld736XH/7wh3zxxRcsWLCABQsWsGfPnk53XgghhBCtq3TrWz1984EKAKYP6YPZ2L2TOJ1ahVRRUUFGRgabNm3iW9/6Vovn3HjjjXg8HtasWRM9NnXqVMaPH89zzz3Xpu9pyyqkgQMHYrfbO3orQvQ6Xq+XI0eOyCokIUSLwhGVnUdqCOu0/bQ3EOKOv+wkGFb5zYLRDEk/+7Y2sVyF1KnsG6fTCWjVZFuzbdu20+qHzJ07l3feeafVa/x+f7M9d1wuV6vnms1mFEWhoqKC9PR0FEVpY++F6J1UVSUQCFBRUYHBYMBisXR3l4QQPVCNN6Bb8AKw7VAVwbBK32Q7g/vE69ZuR3U4gIlEItx9993MmDGD0aNHt3peWVkZmZmZzY5lZma2ulsyaLk2K1asaFM/jEYj/fr149ixYxQVFbXpGiHOBXFxcQwYMKDZJpxCCNEoVtNHl1zQMwYLOhzALFq0iD179vDJJ5/o2R8Ali1b1mzUxuVy0b9//1bPdzgc5OXlEQzqW2lQiJ7KaDRiMpl6xH9EhBA9TzAcoVbHvY9KnfUcOOFGUeCivO5dfdSoQwHM4sWLWbNmDZs3b6Zfv35nPDcrK4sTJ040O3bixIkz7pZstVqxWq3t6pPRaMRo7J5yxkIIIURPUu0JoGed/c0HtKXT4/olkxLXM6at2zX2rKoqixcv5u2332bDhg0MGjTorNdMmzaN9evXNzu2bt26M25SKIQQQoiOq6jTt/bL5m+06aNv5aXr1m5ntWsEZtGiRbz22mu8++67JCQkRPNYkpKSoiuAbrnlFvr27cuqVasA+PGPf8wll1zC6tWrueqqq3j99df5/PPP+d///V+db0UIIYQQvmCYOl9It/b2lrq02i8WIxO7ufZLU+0agXn22WdxOp3MnDmT7Ozs6OONN96InlNcXExpaWn0/fTp03nttdf43//9X8aNG8ff/vY33nnnnTMm/gohhBCiY/RO3t3UkLw7bUgfLKaes2igXSMwbSkZs3HjxtOOfec73+E73/lOe75KCCGEEB1QpePeR95AiM8OVwNwyQU9I3m3Uc8JpYQQQgjRKR5/CG8grFt72w9VEwhHyEm2talwXVeSAEYIIYQ4R8Rq+uiSvJ5R+6UpCWCEEEKIc4CqqlTqOH1U5vRRcKKuofZLz1l91EgCGCGEEOIc4PKFCIQiurXXuHR6bN8kUuN7Ru2XpiSAEUIIIc4Bek4fRVSVjxtrv1zQ80ZfQAIYIYQQoteLRFSqPfpNH3193EWlO0CcxciFua1v2NydJIARQggherna+iChsH57BzRu3DhtcFqPqv3SVM/slRBCCCHaTM/pI28gxPZo7ZeeOX0EEsAIIYQQvVooHKFGx+mj7Ycbar8k2Ria0bNqvzQlAYwQQgjRi1V7A0R03Xlamz66+IKeV/ulKQlghBBCiF6ssk6/0ZcTLh/7y+pQgIuH9qytA04lAYwQQgjRSwVCEVy+oG7tNdZ+GdM3iTSHVbd2Y0ECGCGEEKKXqvL4acM+y20SUdXo9FFPrf3SlAQwQgghRC+l5/TR/lKt9ovdbGTSwJ5Z+6UpCWCEEEKIXsjtD+H2h3Rrb8P+cgCmDem5tV+a6vk9FEIIIcRpTrh8urVV6w3waUPtl8uHZ+jWbixJACOEEEL0MuGISpWOO09v2F9OOKKSl+FgcHrPrf3SlAQwQgghRC9T6fYT1qn4SygS4cN9JwCYMypLlza7ggQwQgghRC+j5/TR50U11HiDJNrNTBnU85N3G0kAI4QQQvQidb4gHn9Yt/b+vbcM0HJfzMbeExb0np4KIYQQgvI6/TZuPFLlYX9ZHQal9yTvNpIARgghhOglQuGIrsm7677Wcl8mDUzt8ZV3TyUBjBBCCNFLVLoDuiXvuv0hPjlYCcDcXpS820gCGCGEEKKX0DN5d/OBCvyhCP1T4xielaBbu11FAhghhBCiF6jzBfEG9EnejagqH3ytJe/OHZmJoii6tNuVJIARQggheoETLv2Sd/9zrJYTLj9xFiMzhvbRrd2uJAGMEEII0cNpybv6BTD/3qsl7868IB2b2ahbu11JAhghhBCih6tw+9Epd5cyp48vj9YCMHtk70vebSQBjBBCCNHDles4fbRu3wlUYHz/ZLKSbLq129UkgBFCCCF6MJeOybu+YJhNBeUAzBmZqUub3UUCGCGEEKIHK9dx6fSWwko8gTAZCVbG9U/Wrd3uIAGMEEII0UMFday8q6oqHzQk784ZmYWhFy6dbkoCGCGEEKKHqtQxebegrI7iai8Wo4FLhqXr02g3kgBGCCGE6KH0rP3y74bCdRfl9cFhNenWbneRAEYIIYTogZz1Qep1St6t9gT47HAN0PuTdxtJACOEEEL0QHom767ff4KwqjI8K4HctHjd2u1OEsAIIYQQPUwwHKHao0/ybigcYf2+xqXTvbdw3akkgBFCCCF6mIo6/ZJ3tx+uxlkfJDnOzKRBKfo02gNIACOEEEL0MCd0mj5SVZV/7SkF4PLhmZgM586f/XPnToQQQohzQLUngC8Y0aWtPcddFFZ4sBgN3ZK8G8tSMxLACCGEED3I8dp63dp654sSAC4bnkGi3axbu22V7ojdXku9fyG4EEIIcY6o8wWp84V0aevAiTq+LnVhVBSuHputS5vtkWAzkRIfuyEYGYERQggheohSp35Lp9/dfRyAi/P6kOaw6tZuW/VPiYtp+xLACCGEED2ALxjWben0kSoPu4prUIBrxuXo0mZ7JNnNJMXFdspKAhghhBCiBzheW4+q09Lpf3ypjb5MHpRKdrJdn0bboX9q7L9TAhghhBCimwXDESrq9Nn36ITLx7ZDVQBcO76vLm22R0q8mQRb7BOGJYARQgghulmZ06db4br3vjyOqsL4/skM6tP12wbEOvelUbsDmM2bNzN//nxycnJQFIV33nnnjOdv3LgRRVFOe5SVlXW0z0IIIcQ5IxJRdStcV+0JsOlABQDXju/63Jc+DgvxXbTTdbsDGI/Hw7hx43jmmWfadV1BQQGlpaXRR0ZGRnu/WgghhDjnlNf5CYb1GX7551elhCLapo3DsxJ1abOtFAX6ddHoC3SgDswVV1zBFVdc0e4vysjIIDk5ud3XCSGEEOcqVVUpdepTuK7OF2T9vhNA94y+pCdYsVuMXfZ9XZYDM378eLKzs5k9ezZbtmzpqq8VQggheiw9tw1Yu7cMfyhCbloc4/ol69JmWxkU6NvFq51iPlGVnZ3Nc889x4UXXojf7+eFF15g5syZbN++nQkTJrR4jd/vx+8/mY3tcrli3U0hhBCiy+lVuK4+EObfe7Xc0gXj+6LEchOiFmQk2rCZu270BboggBk2bBjDhg2Lvp8+fTqFhYU89thj/PnPf27xmlWrVrFixYpYd00IIYToNi4dtw1Yv/8EHn+Y7CQbkwem6tJmWxkNSpePvkA3LaOePHkyBw8ebPXzZcuW4XQ6o4+jR492Ye+EEEKI2Cut1Wf0JRCK8M//lAIwf1wOBkPXjr5kJdqwmLo+nOiWzRx3795NdnbrG0tZrVas1q7ft0EIIYToCvUB/bYN2HSggtr6IKnxFi4e2keXNtvKZFTISY7djtNn/O72XuB2u5uNnhw+fJjdu3eTmprKgAEDWLZsGSUlJbzyyisAPP744wwaNIhRo0bh8/l44YUX2LBhAx988IF+dyGEEEL0Isd1WnkUjqi817BtwPyx2ZiMXTsSkp1k6/LvbNTuAObzzz/n0ksvjb5fsmQJALfeeisvvfQSpaWlFBcXRz8PBAL85Cc/oaSkhLi4OMaOHcuHH37YrA0hhBDifBEIRajUaduArYWVVLj9JNhMXDq8a+urmY0K2Uldn/vSSFFVvbaOih2Xy0VSUhJOp5PExK4tzCOEEELo6Wi1l2M1nR+BiagqP//7fzhWU8+NF/ZnQX7X7nuUmxZHzlmSd2P591v2QhJCCCG6SFjHbQN2HanhWE09drOR2SMzdWmzrSwmA1mJ3ZP70kgCGCGEEKKLlNf5dNk2QFVV3vqiBIDZIzO7bP+hRv1T7V2+2ulUEsAIIYQQXUDbNkCf0Zfth6s5XOnBajJw5ZjWV/XGgsNqIiOhe0dfQAIYIYQQoktU1Pnx67BtQCgS4Y3PtPpoV4/NJslu7nSb7TGwT9dt2HgmEsAIIYQQMRaJqByr1Wfp9Ef7Kyhz+Ui0mbhqTNdu2pieYCHB1rUBU2skgBFCCCFirFyn0RdfMMxbu44B8O38fl26+7PRoNA/tWeMvoAEMEIIIURMRSIqJbVeXdr6154yauuDZCRYmTWia+u+9E2xYzV17YaNZyIBjBBCCBFDZS4fgVDnVx65fMFo1d3vXNi/Syvg2swGsrt52fSpJIARQgghYiQcUTmuU+7Lu1+UUB8Mk5sax/Qhabq02Va5afHdvmz6VBLACCGEEDFS6qzXpe5LRZ2fD74+AcBNkwdgULoumEiym0mNt3Ts4nBI3840IQGMEEIIEQOhcES3ui9v7jxKKKIyMjuRcf2SdGmzLRSlk8umncf068wpJIARQgghYuB4rY+QDqMvxdVePvmmEoDvTh6A0oWjL1mJNuIsHazy666AoD7Jyy2RAEYIIYTQWSAUoUynPY/e+KwYFZgyKJWhGQ5d2mwLs1GhX0oHd5tWVag9om+HTiEBjBBCCKGz47X1hCOdH33ZX+piV3EtBgVuvLC/Dj1ru/6pcR1f6VRXCkF9kpdbIwGMEEIIoSN/KKzLjtOqqvLajmIALh2WQXZyB0dDOiDeaiQjwdqxiyNhqC3Wt0MtkABGCCGE0FFJTT06DL7w+ZEavil3YzEauG5Cv8432A65afEdz7VxHoNwUN8OtUACGCGEEEInvmCY8jp/p9sJR9Toho1XjMnq+DLmDujjsHR8g8hQAFwl+naoFRLACCGEEDo5VlOPqsPoy+ZvKiiprcdhNTF/bNdt2GhQ6Nx+R7XF2hRSF5AARgghhNBBfSBMpbvzoy+BUIS/7dTqp1w7Pod4aweXMXdATrIdm7mD+x0F68Fdpm+HzkACGCGEEEIHR2u8uoy+/HtvGdWeAGnxFuaMzOp8g20UZzHStzOJwjVF6PIDaCMJYIQQQohO8vhDVLkDnW7HWR/knd1aDsl3LuyHxdQ1f6YVBQand2K/I38deCr17dRZSAAjhBBCdNLRGn0qzr62/QjeQJiBaXFcPDRdlzbbIivRRoKtg4m7ANWH9etMG0kAI4QQQnSCyxekxtP5ZcP7y1xsbtgy4AczBnXZ7s9Ws6FzibveavA59etQG0kAI4QQQnSQqqocqez86Es4ovKnLUWAVrQuLzOh02221ZA+DoydCZZqinTrS3tIACOEEEJ0UHmdH7c/1Ol2/r23jKPVXhxWEzdN7rotAzISrSTFdWLqyF0OAY9+HWoHCWCEEEKIDgiFIxyt7vzoS7UnEF02fdPk/iR2JhelHSwmhdzOTB1FIlAT2w0bz0QCGCGEEKIDjtbUEwx3ftnwX7YfoT4YZmiGg0uHZejQs7YZ1MfR8c0aQduwMaTPjtsdIQGMEEII0U7eQEiXDRv3lDjZVliFojQk7nZ0/6F2SnNYOrc9QTgEzqP6dagDJIARQggh2ulwpafTNdtC4QgvbtWWH88ekcmgPvE69OzsTEaFgWmd/C7n0S7ZsPFMJIARQggh2qHK7cdV3/nE3X9+VcrxWh+JdjP/58KuS9zNTYvrXIG8gKfLNmw8EwlghBBCiDaKRFSO6JC4W+n28/YXWhBw8+QBXbbfUXKcmYwEW+caqSrs0i0DWiMBjBBCCNFGJbX1+IORTrfz521H8IciDM9K4OK8Pjr07OyMBqXz01R1J7qlaF1LJIARQggh2sAXDHO8tr7T7ew+WsuOomoMCnx/xiCULkrc7Z/aiZ2mQUvc7aaidS2RAEYIIYRog+JqL5FOzpwEQhFeakjcnTc6mwGdqcPSDgk2E1mJnZw6qimCcOc3rNSLBDBCCCHEWTi9QV12m37vP8c54fKTEmfmhgn9dOjZ2RkNCkPSHZ0b6fG5wF2mX6d0IAGMEEIIcQaqqlJU1fly+SdcPt7drSXu/tfUXOyWTkzntMPAtLjOfZeqQnXPSNxtSgIYIYQQ4gzKXD68gXCn2lBVlZe3FhEMq4zum8TUwWk69e7M0hwWMjo7dVRXCn63Ph3SkQQwQgghRCuC4QjHajqfuPvJwUq+OFqL0aCwcPrALknctZoNDO7sqqNQoFv3OzoTCWCEEEKIVhRXewl1cr+jKrefl7YWAXD9hH70Tbbr0LMzUxQYmtHJvY4Aqg9BpPNF+2JBAhghhBCiBW5/iIo6f6faUFWV/918CG8gzJD0eK4Zl6NT786sb7K987ta19eAp0KfDsWABDBCCCHEKSIRlcJyd6fzVj/cV85/SpyYjQr/PXMoRkPsp44SbCb6pXRylCcSgapD+nQoRiSAEUIIIU5RUlvf6cTdEy4fr27X8ke+O3lAl0wdmYwKQzM6uWQawHUMgp3fMiGWJIARQgghmnD7Q5R0suJuJKLy7MZC/KEII7MTmTsqS6fendngPvGdq7YLEKyH2qP6dCiGJIARQgghGug1dfT+nlIKTtRhNxu545LBGLpg1VFGopU0h7XzDVUVgtr5/Z5irWu2vxRCCCF6AT2mjo5We3njM20E47+m5pLe2d2f28BuMTIwrZNLpgE8lVryrh4iYfj0j/q01QIZgRFCCCHQZ+ooFInw7KZCQhGV/P7JzByWrlPvWmdQIC/D0fkE4XBQG33Ry96/w9fv6NfeKSSAEUIIcd5TVX2mjt75ooTDlR4cVhO3f2twlxSsG5AWR7xVhwmVym/026yxYj988Rd92mqFBDBCCCHOe8dqOj91VFjh5u0vtL2OfjBjIClxFj26dkbJcWayk3RY3VR3ArxVnW8HIOCGzQ9reTSDZ+rTZgvaHcBs3ryZ+fPnk5OTg6IovPPOO2e9ZuPGjUyYMAGr1crQoUN56aWXOtBVIYQQQn8eHaaOAqEIz24sJKLCtMFpTBvSR6fetc5qNjA0w9H5hoI+reKuHlQVtj0D7hPgyIRpP9an3Ra0O4DxeDyMGzeOZ555pk3nHz58mKuuuopLL72U3bt3c/fdd3Pbbbfx73//u92dFUIIIfSkqiqFFZ2fOvp/nx+lpLaeZLuZ788YqEvfzsSgwAWZCZg7u1WAqkLlAf22Czi4Doo+BsUI3/oZWHVILG5FuyfNrrjiCq644oo2n//cc88xaNAgVq9eDcCIESP45JNPeOyxx5g7d257v14IIYTQzbGaejz+zk0d7St18f5XpQDc/q3BJHS2hH8bDE534NAj78VVAj5n59sBrXbMjue11/n/BenDIIbbKMU8B2bbtm3MmjWr2bG5c+eybdu2Vq/x+/24XK5mDyGEEEJPekwdeQMhnttUiApcOiydCQNS9OncGeQk20hP0KHeS8Cj307T4QBs/h2E/JA9HkZfp0+7ZxDzAKasrIzMzMxmxzIzM3G5XNTXt/yLs2rVKpKSkqKP/v37x7qbQgghziN6TB2pqsrzmw5RXuenj8PC96bm6tfBViTHmRmQGtf5hiIRqCjQr2Dd5y9CTRHYkuCiJaDEfo1Qj1yFtGzZMpxOZ/Rx9GjPL2kshBCi99Bj6uj9r8rYUVSN0aDw48svIM4S29qwNrOBPD32OQKoPaKNwOiheDvsf097PeNuiEvVp92ziHkl3qysLE6cONHs2IkTJ0hMTMRub3npl9VqxWrVYXhMCCGEOIUeU0f7y1y8tkObfrllWq4+q4HOwGhQGJaVgKmzSbug5by4SjrfDmiVe7c+rr0eeS30m6RPu20Q8xGYadOmsX79+mbH1q1bx7Rp02L91UIIIUQz4YjKN50sWFfrDfDE+m+IqDBjSBqzR2Se/aJOGpIer88ITyQMFQfo9LKrxrY+fgT8dZA6BCYs7Hyb7dDuAMbtdrN79252794NaMukd+/eTXFxMaBN/9xyyy3R8++44w4OHTrEz372M/bv388f/vAH/t//+3/cc889+tyBEEII0UaHK93Ud6JgXTii8tSGg9R6g/RNtnPbxbGvttsvxa7PJo2gbRUQ8unT1ldvwok9YLJpS6aNsV991VS7A5jPP/+c/Px88vPzAViyZAn5+fksX74cgNLS0mgwAzBo0CD++c9/sm7dOsaNG8fq1at54YUXZAm1EEKILlXu8lFR17lS+W/uPMrXpS5sZgP3zL4Am9moU+9alhJvpr8eSbsAniqtwJweyr+GL1/TXk/9b0jqq0+77aCoqh7jSLHlcrlISkrC6XSSmJjY3d0RQgjRy3j8IfaUOIl04i/eziM1PPJBAQB3XZbHtCFpOvWuZXaLkdE5ifrkvYQCcHyXtmFjZ/nd8N6d4KnQtgq46CfQyiiUK2Qm6YJpMfn73SNXIQkhhBB6CUdUDpyo61TwcsLl49mNBwGYNyor5sGLyagwLFOnpF2Aqm/0CV4iYdj8ey14SciGKf/TavASaxLACCGEOKcdqnDjC3a83kkgFOGJ9d/gCYTJy3Bw85QBOvbudIoCQ9Md2C06TU/VFoO3Wp+2dr6kjeQYrXDJvWDRaXqrAySAEUIIcc4qc/qodHcu7+XlbUUcrvSQYDPx48vz9BsVacWgPvGkxOu0k7W3Wgtg9HBwPXz9tvb6onsgbYg+7XaQBDBCCCHOSW5/iCNVnSvWtulAORv2l6MAiy8dqt9qoFb0S7GTmWjTp7FgfUO1XR1SXcv3wbantNdjb4KBF3W+zU6SAEYIIcQ5JxSO8E0n816OVHn4v58cBuCGif0Y2y9Zn861IiPRqt+Ko0hECzr02GXaUwkfPai1NWAajP//Ot+mDiSAEUIIcc45VOnpVN6Lxx/i8Q+/IRhWGd8/mQX5sV0mnBJvZnCfeP0arPpGn60CQj7Y8Gvw1ULKwC7b56gtekYvhBBCCJ2UOuup6kTeSygc4dF1Byhz+ejjsPA/M4dgiOFKmwSbiQsyEvQriOc6Du7yzrejqrDlCaguBGsiXPYrMLe8BVB3kABGCCHEOaPOF6S4ytvh61VV5X8/PhQtVrd0zjASbLGrMGu3GBmWlYDBoFPw4nNC9SF92vrqTSj6GBQjzPwFOGK/ZUJ7SAAjhBDinBAKR/im3N2pvJe/7yrh428qMShw9+UXkJum47TOKSwmA8OzEjDrtaopFIDy/fok7RZ/Cl+8or2e+t+QNbrzbepMAhghhBC9nqqqHKxw4+9E3sumAxX8fdcxAH5w0SDG9U/WqXenMxkVRmQn6LcVgapCxT4Id27JOAA1RfDxau31sKvggnmdbzMGJIARQgjR6xVVeanxdLzS7J4SJ3/crE29XDs+h8uHx266xKDABZkJ+uwu3aj6EPhcnW/H59SSdkP1kDUWJt/e+TZjRAIYIYQQvVqps54yZ8d3WD5a7eWxDw8QVlWmDUnj/1zYX8feNacoMDTDQZJdx7wad7mWuNtZkRBsXKVt+JiQrVXaNegYZOlMAhghhBC9Vo0nwJFOJO3WegP8/t/78QbCDMtM4I5vxXbF0cC0eH2L4fndUHWw8+2oEdj6NJzYo600uvSXYOvZmydLACOEEKJX8vhDfFPu7nDOqi8Y5vf/LqDSHSA7ycZP5lyAxRS7P4sD0uLIStKpyi5oSbsV+7UNFjtDVeHzP0Hhh1qNl4t/Cim5+vQxhnru2JAQQgjRCn8ozP6yOsIdXHIUiag8teFgdI+jn88bHtPl0rlpceQk61hDJRKG8r3adgGd9dWb8PU72uvpP4b+kzvfZheQERghhBC9SjiicqDMTSDUsRVHqqry8rYidhXXYDYqLJ0zTL/9h1qge/Ciqto2AX5359sqeP/kcukLb4Ohl3e+zS4iAYwQQoheQ1VVDpa7cfs7vsfPv/aU8cHXJ1CARZcO5YLMBP06eArdgxeAym+gvqbz7RzeDJ8+q70ecyOMWtD5NruQBDBCCCF6jSNVXqo9Ha91suVgJX/59AgAN0/JZcqgNL26dpqYBC81Rdoqoc4q2QmfPAqoMOxKyP9e59vsYpIDI4QQolcoc/oo7cRy6U8PVfGHjQdRgTkjM7lyTJZ+nTvFwD5xZCfpHLy4SqH2aOfbKd8HH/1WWzY98Fsw5Q5tfXcvIyMwQggherxab4Ciqo7vrvxZUTVPbzhIRIVLLkjn1ukD9ds88RQxCV48Vdqmip1VfRjWPwBhP/SdCBfd02N2l26v3tlrIYQQ5w2PP8SBEx1fLr3rSA1PrP+GsKpy0dA+/OjiwTGr9RKT4MXn1JZLd3aPI1cpfLgcAh5IHwEzl4ExdiuvYk0CGCGEED2WL9i55dK7j9ZqVXYjWpXdOy4Zot/Oz6eISfAS8GpTPmrH93gCwFsN636lJf+mDITL7wdT7FZedQUJYIQQQvRIvmCYr0tdHV4u/VWJk0fXFRCKqEwelMqimUMx9qbgJRSAE3sh3PE9ngBtufW6X4G7TNsiYPavwerQp49nE8PpKUniFUII0eP4Q2H2lbo6vLv018edPPLvAoJhlQtzU7jzstgEL4oCg/vEk6F3HZlwSCvrH+p40jKgbfD44XKoPQL2VC14safo08ezsTrAFrt9pSSAEUII0aMEQhH2ldbh62Dwsr/Mxe//XUAgHCG/fzJ3XZ6HyaD/SIDRoJCX4SAl3qJvw5EIVOzTclU6o74GPvilFrxYE2H2SkiI3cqrZuwpkDEC3J28hzOQAEYIIUSPEQxH2Ffqoj7Qsf19Dpyo43dr9+MPRRjbN4m7Z12A2ah/8GI2KgzLStB/+4FIBMq/hvrazrXjqdCCF1eJNvIy50FIjt1oSDOOTOiTF/Ol2RLACCGE6BFCDcGLt4PBy8FyNw/9az++YITROYn8ZM6wmGzOaDUbGJGViN1i1LfhSETb36izwUtdGfz7F+Aph/gMLXhJzNali2eV3F9LEu4CEsAIIYTodlrwUofH37HgpbDCzUP/2kd9MMyI7ASWzo1N8OKwmhiWlaB/25GwPiMvzqPw7/ugvhoSc2D2b8CRoUsXz0hRIHVI1wVKSAAjhBCim4UjKvvL6jq8v9F/jtXy6LoD+EMRhmUm8LO5w7GadB4dAZLjzFyQmaB/MnAkrK028jk71071IW21kc8JyQO0kZeuSNhVDJA+HOJjty1DSySAEUII0W0iEZX9ZS7qfB0LXrYcrOTZTYWEIyqjcxJZMnsYNrP+wUt6goUh6Q79q/fqFbxUFMCH90PArY2EzF4JtiR9+ngmRrOWrNsV33UKCWCEEEJ0i0jDyIurvmPBy7/2lPLKNm1jxmmD0/jvmUNikrDbN9nOgLQ43dvVgpc92lLnzjixB9avgGC9NhIy6wGwdEGdF5MNMkeBJQY/m7Z8fbd8qxBCiPNaJKJyoLwOZ337i7Spqsrrnx3lH18eB2DuqCxumZar+/YAigID0+LJSopBxVq9gpfjX8CG32h7G2WNhct+BWadC+q1xBIPmaPBpPMS8naQAEYIIUSXCoUj7C+r69C0UTii8sePD7HpQAUAN07qz7XjcnSf2jEZFYamx6DGC5wsUuev61w7xZ/Cpoe0XaX7XqjtbWSy6tPHM4lP15ZJG/SfqmsPCWCEEEJ0GX8ozP7Sug4tlfaHwjy5/ht2FddiUOC2iwdz6TD9V9jYLUaGZSbov0wa9AleVBW+fhc+/7+ACrnT4eKfxn5jRkXRlkgn9Yvt97SRBDBCCCG6RH0gzL6yjm0P4PaF+P2/9/NNuRuzUeHHl1/AxFz9V9ikxlsYkh6PKQa5NISDWsJuZ4KXSAi2Pw8H/qW9z5sLU/8n9qMhRjOkD+u6bQjaQAIYIYQQMef2h9hf6iIYbv+u0lVuP6v+tZ+S2nrirUZ+Omc4w7ISdO9jvxQ7/VNjlJAa9GnBS9Db8TYCbtj4Oyj9AlDgwh/AyAUxr3iLJR4yRoK5Z+1eLQGMEEKImKr1Bjhwwk040v7g5Wi1l4fW7qfaEyA13sK984brHmTENN8FtETd8n0QDnS8jboybaWR86i2+ufin8KAKfr1sTU9JN+lJRLACCGEiJmKOj+FFW7U9scufHa4mj9sOogvGKFvsp17rxhOH4e+SaoxzXcB8FRBZYG26qijyr/WVhr5XRCXBpcth7Qh+vWxJT0s36UlEsAIIYSIiVJnPUWV7Z8yiagqf995jLe+KAFgZHYid8/K033jxJjmuwA4S6DmMB2K3hod2ghbHtdyX1KHaMuk4/vo1cOW9cB8l5ZIACOEEEJ3xVVeSmrr232dNxDimY8Osqu4FoArRmdx85Rc3cv3xzTfRVW1sv6u451r48vX4Mu/au/7T4WLl8Y+D6WH5ru0RAIYIYQQulFVlcIKDxV1/nZfW1Jbz6MfFHDc6cNsVLj94sFcnJeua/8sJoUh6Q6S42KU7xKJaFNGnsqOtxEOwJYn4PAm7f3o62HCrdqeQ7GUmAMpg8AQ4+/RiQQwQgghdBEIRThwomMF6nYeqeGZjw5SHwyTGm9hyewLGJKubzn85DgzQ9IdMdmlGtBnmbS7HDb9TguCFKO2RPqCufr1sSVGi5aoG5ca2+/RmQQwQgghOs3tD1FQVkcg1L4aLxFV5e0vSvjbzmMADM9K4MeX5+k6QmJQYEBaHNlJMSyxH6xvWCbd/mmzqKM74JNHteXSlniY+QvIHqdfH1tiT4E+F3TrlgAdJQGMEEKITimv83G4wkN7V0l7AyGe3VjI50dqAJgzMpP/mpaLSccpjDiLkaEZDuKtMfxzV18LFfu1EZiOiIRg1yuw9y3tfVoeXPJzSMjSrYunUQwNq4z6xu47YkwCGCGEEB2iqipFVV7KnL52X1taW8/qdQcoqa3HZFD4wUWDdN8WIDPRysC0eAw6JwA3U3sUao90fKWRpwI2/R4q9mnvR1wDE78f220BLHHartWW+Nh9RxeQAEYIIUS7BcNavourvn35LqqqsvFABS9vLcIfipASZ2bJ7AsYmqFfZV2zUWFwuoPUWBWmA21Po8oC8FZ3vI1jn2lTRv46MMfDjB9r+xrFUi9L1D0TCWCEEEK0S0fzXdy+EC98cojth7U/+iOyE7jrMn3zXZLsZoZmxDBRF8Dv1qaMOprvEgnBF3+BPX/T3qcNhUvuje2UkdGs5br0skTdM+nQP+FnnnmGgQMHYrPZmDJlCjt27Gj13JdeeglFUZo9bLaev75cCCHE6Srq/OwtcbY7eNl73MnP3/oP2w9XY1QUbprUn19eOVK34MWgQG5aHCNzEmMbvNSVQemXHQ9ePJXw71+cDF6GXw1XPBzb4CW+D+RMOKeCF+jACMwbb7zBkiVLeO6555gyZQqPP/44c+fOpaCggIyMlucvExMTKSgoiL5XYr3xlBBCCF2pqsqRKi+l7cx3CYUjvLnzGO99eRwVyEq0sfiyoboukU6ymxmcHo/NHMP9eiIRqDoI7hMdb6NkJ3y8WtsSwBwH0++CgRfp18dTmaxa9d74tNh9RzdqdwDz6KOPcvvtt/P9738fgOeee45//vOf/OlPf+Lee+9t8RpFUcjKimF0KYQQImbqA2EOlrtx+9uX73K8tp6nPzrI4UoPAJcNz+C/pubqFmiYjAq5qXFkJMZ4VD9Yr23GGPB07PqAF3b+CQ6s1d6nDtGmjBKz9evjqRKytVVGxnM3U6RddxYIBNi5cyfLli2LHjMYDMyaNYtt27a1ep3b7SY3N5dIJMKECRP47W9/y6hRo1o93+/34/efrOLocrna000hhBA6OeHycaTK266dpFVVZUNBOX/edgR/KILDauJHFw9m0iD9pjDSHBYGpsXHdroIGjZjPKDlrXTE8d2w9QlttRHA8Plw4fe14nGxYI7TNnq0J8em/R6kXQFMZWUl4XCYzMzMZsczMzPZv39/i9cMGzaMP/3pT4wdOxan08kjjzzC9OnT2bt3L/36tbzL5apVq1ixYkV7uiaEEEJHwXCEQxUeqj2Bdl3n8gX54+ZD0douo3MS+e+ZQ3VbEWQxGRjUJz62K4xAmzKqOdzx/YyCXvi8yaiLI1NbZZQ1Vr8+NqUokNgXknPPiRVGbRHzsaVp06Yxbdq06Pvp06czYsQInn/+eX7961+3eM2yZctYsmRJ9L3L5aJ///6x7qoQQgig1hugsMJNINS+UZdth6p4edsRXPVBjAYtUffKMdkYdMp7zEy0MiA1Lna7RzfyuaDyGy0I6YjTRl2u1vYyMseoErDVoRW/s+q79UJP164Apk+fPhiNRk6caJ7EdOLEiTbnuJjNZvLz8zl48GCr51itVqxWa3u6JoQQopMiEZUj1e0vTFdR5+P/fnKYL485AeibbGfRpUMZ1EefQml2i5FBfeJJssewuBtooy61R8BV0rHCdF096mIwQvIAbeTlPFwc064AxmKxMHHiRNavX8+CBQsAiEQirF+/nsWLF7epjXA4zFdffcWVV17Z7s4KIYSIDY8/xMFyN95AuM3XhCMq/9pTyt92HsMfimAyKHw7vy/zx+Vg1mGUxGhQ6JtiJzvRFttquqAVk6s8oCXcdkRXj7o4MrUk3V64h5Fe2j2FtGTJEm699VYuvPBCJk+ezOOPP47H44muSrrlllvo27cvq1atAmDlypVMnTqVoUOHUltby8MPP8yRI0e47bbb9L0TIYQQHXK8tp6j1d527WVUWOHmhY8PUVSl/cEfkZ3AbRcNJie583+wFQXSE6z0T4mLfZKuqkJtMTiPdmzUxV8Hu17uulEXWyKkDgarfpWLe6t2BzA33ngjFRUVLF++nLKyMsaPH8/atWujib3FxcUYmiQQ1dTUcPvtt1NWVkZKSgoTJ05k69atjBw5Ur+7EEII0W7eQIhDFR7qfG1fYeMLhvl/nx9l7d4yVBXirUZunpLLzAvSdanxlWAzMbBPPI5Ybr7YKODRRl387vZfGwnDwXXaJoz+hpWyw66CiQtjM+pismpbADjS9W+7l1JUtaM7UHUdl8tFUlISTqeTxMTE7u6OEEL0auGISklNPced9e0adNhVXMOLWw5T6dZWJs0YksZ/TRuoS26K1WxgQGocfRxdkP+oquA8po28qO2rKAxo2whsf04rbAdaHsqU/4asMfr2E7Q8l8S+kNRPe93LxPLv97lb4UYIIcRpar0BDlV68Afb/oe73OXj1R3F7GjYwyjdYeUHFw1ifP/kTvfHaFDISbaRk2SPfZ4LaFM+VQc7NupSXwM7X4bCD7X35jgYfzMMvwoMMfhzGp8OqYO00RdxGglghBDiPBAIRSiq8lDlbntdF48/xDu7S1i7p4xQRMWgwJVjsrl+Qj9dqummJ1jonxqH1dQFIwvhINQUaVsBtHfiIRKCfWvgy9dOLq0eOktL0rWn6N5VbIlagq4tSf+2u1g40vak8PaSAEYIIc5hqqpywuXnaI2XULhtf7hDkQjr95Xzt53HotsHjO6bxPemDCA3rfNLo1PjLfRPtRNn6YI/QaqqbcBYe0QLYtqr9EvY/jw4i7X3aXkw5f8H6cP17SdodVySc8+JTRfDkTDHPcepcFbE7DskgBFCiHOUxx/icGXbk3RVVWVXcS2vbT/C8YZaMH2T7dw8ZQDj+yd3Okk3Oc5M/9S4rknQBa0gXXVhx6aLXMe11UVHtmjvrYnaiEvebFB0XhlliYPkgefMpouV9ZUcqztGMBLEoMZuFZkEMEIIcY4JhiMcq6nnhMvX5tmSw5Ue/vLpEb4u1VbUJNpM3DCxP5cNz8DYydyURLuJ/qlxJNpiXIiuUeN0UV1Z+6/1VMJ/XodvPtASfBUDDLsSxn9P/0q35jgtAfgcWVnkDrg5UncEb0crGLeTBDBCCHGOiERUSl0+jtfWt3m6qNoT4P99fpTNBypQAbNR4YrR2Vw7PqfTUzwJNhP9U+JIiuuiwEVVoa4Uao60f/NFnxO++hvsXwORhqmmvhdqoy6pg/Ttp8nWELhknBMVdAPhAMfqjlHlq+rS75UARgghejlVValw+zlaXU8g1LbVRTXeAO99eZwP950g2BDsTB+Sxk2TBpCe0LlVL/FWI/1T4kiJ9YaLTXmrtVGXgKd91wW88PU78PXbEKzXjmWMggm3QOYofftosmnLoROyzonAJaJGKPOUUeopJdKR5eidJAGMEEL0YrXeAMXVXjz+tq32qPZogcv6/ScDl2GZCXxv6gCGZnSuumui3UROkr1rAxefUwtcfK72XRfyQ8H78NWbJwvRpQ7RApecCfoGGJY4SOqvLYs+BwIXgKr6Ko65jxEIt2+3cj1JACOEEL2Qxx+iuNpLrbdtK2uqPQHe3V3CRwXlzQKX6yb0ZUzfpA4n6CqKtqooO8lGQlfluICWmFtTpNVmaY9wEArXw5d/BW/DlEdiP8j/HuRO1zdB15akjbicA6uKGtX4ajjuPo431DV5LmciAYwQQvQi/lCYo9X1VLr9bUrQrXL7effL43y0v5xQw2ZHw7MSuH5CP0blJHY4cDE07FeUk2zXpSZMmwXrtRwXTzuX5wa9cODf2nRRY+ASnw7jvgtDLtevyq2iQFyaVj3Xdu5Ujq/11VLiKemyBN22kABGCCF6AV8wTEltPZV1/jZtuljp9vPu7hI2FlQ0C1xumNiPkdkdD1xMRoWsRBuZibbYb7TYVMivlf5vbyG6+hqtCF3BmpP5MfZUGH29trrIqNOokWLQNnJM6hu7Hai7gdPvpMRdgifYztyiBoFI7KaYJIARQogerD4QpqTWS6U70Ka/24UVbv61p4xPC6sIN1wwMjuR6yf0ZWROxyu7Ws0GspNsZCTYOr2sul1CfnCVgKu0ffsWuUph71vadFFjnkZiPxh9HQy+VL/AxWjRknITssHUhbk/MeYKuCipK8Ed7EANHbTA56OjH/FhwYc69+wkCWCEEKIH8vhDlNTWU+05e+ASiah8fqSG978qpeBEXfT4qJxErp/QjxHZHZvKUBSt+Fxmgo3kOLMuu023WcCrBS7u8vYFLlWFsOdvWgG6xuv6DNNGXAZM1S/HxZYEiTnadNE5kpgLUBeoo8RdQl2g7uwnt+CI6wjrjqxje+l2wmqYcFC2EhBCiPNCnS9ISW09NZ6zJ+d6AyE2FlSwdk8ZFW4/oG2OOG1wGleMzmJwescKr5mNChkJNjISrV2b3wLaaiLnMaivbvtUUSQMx3fB1/+A0i9OHu87EUbfAJmj9QkyDCatdktCtray6Byhqio1/hrKPGUdmiqKqBF2l+/mgyMfcKDmQPT4kOQhzBk2h0Us0rO7URLACCFED+CsD1JSU4+z/uyBywmXj7V7y9hUUEF9w//hOqwmZo3IZPbITFI7uIw5wWYiM9FGWryla3aGbspbrQUuPmfbr6mvhYPr4MBaLTcGtBGWQd+CUddB6mB9+maJ14IWR4Z+yb49QDgSpqK+ghPeEx1aDl0fqufjYx/zYfGHVNZXAmBUjFyYeSGzcmcxJHkIBp9sJSCEEOecSESl0u2nzOU7ax2XiKqyp8TJh/tO8PmRmujgRN9kO1eMzuKivD4d2tXZaFDo47CQmWgjvqv2KGqkqtpqIuexthegU1Wo2Af7/6lNEzVW3LU4tB2iR8zXkmk7y2DUpocSss6JXaGbCoQDnPCeoMJbQVht/xRPqbuUjcc28vGxj/GFtT2z4s3xzOw3k0sHXEqqrWuWjUsAI4QQXcwXDFPu8lNe54vWZGlNpdvPxoIKNh0op9J98v+Sx/ZL4srR2Yzt1/4aLooCiTYz6QlWUuMtXZuUC1otlroy7RHyte2aoBcObdSKz9UUnTze5wJtNdHAi8HUuQrCgBasODIhvs85NdoC4A16KfOUdajkvz/kZ8eJHXx87GMO1h6MHs+Oz2Z27mym5UzDatTh598OEsAIIUQXcdYHKXP6qPGeOTE3FI6ws7iGj/aX859jThpPjbcYuSgvncuHZ9A/tf05GHaLkT4OC+kJ1g6N1nRafa0WtHir2paYq6pQfUjbWPHQhpOl/o1WGHyJFrikDe18v8x2bXooPgPMts6314M05rdUeCtwBdpXrVhVVQ47D7O5ZDM7SndER1sMioGxfcZy2YDLGJk2EoPeu3O3kQQwQggRQ+HGaSKnD2/gzMP1JbX1bCwoZ/OBCly+k5sRjsxO5LLhGUwamNru2itmo0JqvBa0dGml3EbhEHjKtU0WA20sguaphMMboXCDVvulUWJfLWgZcnnnd4Y2mLRRFkfGOTdFBOAL+aior6CqvopgpG3Vmhu5A262lW5j87HNlLhLoscz7Blc3O9iZvSdQbI1Wecet58EMEIIEQN1viAVdX6qPIEz7gxd5wuy43A1H39T2WwJdLLdzCXD0pl5QQZZSe0bFTAokBxnoY/DQkpcNyTkglbqv65Uy3GJtCHPIlgPxdu0oKX0S2gcdzKYof8UGDYPssZ1bjWRYtDK+sf10fJbDN0zchArETVCja+GivqKdi+DDkVC7K3ay7bj29h1YhchVQugzQYzF2ZeyMX9LmZYyrCuXUp/FhLACCGETvyhMJXuABV1furPMNpSHwjz+ZFqthZW8dUxZ7TgnEGB8f1TuHR4Ovn9U9qVm2JQICXeQmq8FrR0eV4LNIy2VGgrgvxt+AMaCUPZV1qxueJtzfNhMkbBkEth4EVagm5HGYxgT2kIWlLPubwW0FYDVXgrqPJVEYqEzn5Bg4gaYX/1fnaU7mDniZ14QicTqQckDOBb/b7F1OypxJl75pJxCWCEEKITIhGVaq8WtDjrg63mtgRCEb44WsO2wip2Fdc0S94dmBbHtCF9uGhon3YtgTYaFJLjzN0btEQiWs0Wd7lWtv9suS1qBCr2ayuIDn+sXdsoIRuGXKZVyk3I6nifGoOW+HTt+RwMWoKRIDW+Gqrqq9pVLTeiRjhUe4jtZdv5rOyzZnkxSZYkLsy6kIv6XkRuYm4suq0rCWCEEKID3P4Q5S7fGaeIQuEIe4472VpYxedFNdGaLQA5STamDenD9CFp5CS3fe+cxqAlLd5CcncFLaDVa3FXNEwRneX/+iNhOLEHjmzVRlqaBi2WeBj4LS1wSR/e8Skig6khaOlzzgYt4UiYWn8tVb4qnP6218tRVZXiumJ2lO5gR9mOZquQ4s3xTMycyJSsKQxLHaZrQm6cOQ47sdsXSgIYIYRoI7c/RLU7QJXHjy/Y8kiDxx9i99FadhbX8OXR2maJu30cFqYNTmPakD4MTItrcz6B1Wwg2W4mJc5Ckt3cPTktoCXhesq1wOVsy5/DQSj7jzbSUvwp+JusgDHHQb9JkDtDe+7ovkRme8P0UCrYks+pkv6NVFXF6XdS5aui1l9LpI3bKoQiIb6p+YYvyr9gd8XuaKE5AKvRyoTMCUzJmsLItJGYDPqEAgbFQIIlgRRrCknWJCxGCy5X+1Y+tYcEMEIIcQaNQUulx4+/laClos7PziM17DxSzb7SumhOC0CS3cyUQalMH9KHvEwHhjb8kVUUrbJuSryFlDgzcZZu/E91wKOtCvJWnb3YXMALpbu1UZajO6BpWXprAvSfCrnTIXt8x4IWRQFrohaw2FPPqXL+TamqSl2wjqp6LWhpa16LN+hlT+Uedlfs5j8V/8EbOrnqy2wwMy59HJOzJjM2fSwWoz4bT5oNZpKtySRbk0m0JnbpkmoJYIQQ4hR1viDVngBVnkCLQYuqqhyu9DQELTUcqW6+PLhvsp2JuSlcmJvCkIy2BS0mo0Ky3UxynIXkODNmYzeukPG5wFuplfdvrL3SElXVljmXfK49TnwNTSu72lNgwDQtaMkc07FpHaMF7MlawGJPAeO5+WcrHAnjCrio8dfg9DvbHLRU1Vexu3w3uyt2s796f7PKugmWBMaljyM/I5+RqSOx6lDoT0HBYXGQaEkk2ZrcrQm+5+ZvghBCtEMkouLyBanxBqnxthy01HgDfHXMyVcl2qPpnkWKAsMyE5iYm8LE3BSyk84+729QwGEzkWQ3k2Q347Caum+JqqqCr1YLWLxVEPK3fm6wHkr/AyWfQclOLQemqcS+0PdCLWhJH97+oMVg1EZZ7MnatFBn6730YMFIkFpfLbX+Wpx+Jypn37zSH/ZzoPoAe6v2srdqb7M6LaBVxh2fPp78jHwGJw/WZUTEbrKTaEkkyZqEw+zA2EPyiySAEUKcl/yhMLUNAYurPkQ4op72+f7SOv7TELAcPWWUxWoyMK5fMhNzUxg/IJnENhSJi7caowFLgs3cfQm4oOWo1NecfIRbKXYWCUPNYS1oOb5LS8ZtOjpgtEDWGC1o6XshJGa3rx+Koi2TjgYsiedcfZam6kP1OP1Oanw1bVo9FFEjFLuKowHLwZqD0RotoI2I5KXkMT59POMzxpMV34nVWw3MBjOJ1kQtaLEkYe5ojlKMSQAjhDgvqKpKnT9ErSdIbX3gtM0TwxGVoioPe4+7+OpYLQUn6potdVaAgX3iGdsvibF9k8jLTDjrNI/dYiShyShLt04LqaqWSNsYsPhb+eOpRqDmiJaAW/YVnPjq9NwXRxb0awhYssa0bw8iRdFWHlkTtYDFlnTOTguBlkxbF6jD6XfiCrjwh88wuoX2e1rlq+Lrqq+jj1MDnTRbGqPSRjGyz0hGpo7E0Zk6OYDFaCHBnIDD4iDBkoDdFLuVQ3o6d39rhBDnvfpAGJcviLM+iKs+2CwgCYYjFFa42V9ax74yFwdO1J22sig13sLYvkmM7ZfEqL5JZxxlURSIt5hIsDU+zO0u+6+7kL/JKEtty8udVVXbDbrsPyeDFv8pK0fMcZA5Sku+7XshJOa0fcWPYtASeG2JDUFL0jm5xLmRqqq4g25cARcuv+usoyyqqnLcc5wDNQf4puYbCmoKqPHVNDvHZrQxPHU4o/qMYlTaKDLjMjs13Wgz2kiwnAxYunoTRr1IACOEOGf4glrA4qoP4qwPEQhFmn124EQd+8vq2FfqorDCfdpO0HEWI8OzEhjTN4kx/ZLJSbK1+oeiMYcl0WaOBizdOiUEEApo9VkaH8EW9h4KB6DqIJTv0wrKle/T8l+aMtkgYyRkjYXssZA6pO1Bh8HUJFhJBEvCOT0lBNq0kMvvwhVwUReoa5ZIe6pwJExxXXE0WPmm5pvTghyjYmRg4sBowDIoaVCHlzorKMSZ43CYHdGgxWzomVNC7SUBjBCi1wqEIs1GWBpHUCKqSpnTx8FyNwcr3Bwsd3OkysMpaS4k2s2MyEpgeFYiI7IT6J8S12qNFbvFiMNq0h42E/EWY/fvCxPynxKwtLBiyFsNFfsaApZ9UFV4+kiMwQwZI04GLGl5bVvmrBgapoMSTj7MvWP6oaNUVdUCloALd8CNO+g+42aJTr+TQ85DHHYe5pDzEIW1hadNI1kMFoYkDyEvJY8LUi5gSNKQDq8YshltxJvjiTfH47A4iDO1vd5QbyMBjBCiV1BVFW8gTJ0vRJ0vSJ0/FF0t5KwPUtgkWDlU4cbTwl5EfRwWRmQlMjw7kRFZCWS1MsJiMSk4rGbirUYSGp5N3Zm/0ijg1fYY8rtaDlgCXqgu1EZYqg5qIyzuE6e3Y0vWVghljNAeaUO1ZNyzMdu1IMXiOPl8jo+uRNQInqAHd0CbFvIEPa2OsPjDfo64jnDIeYhDtVrQ0rTqbaM4U1w0WLkg5QJyE3M7NMJiNpiJM8dpwYrZQbw5XreidL3B+XOnQoheJRSO4PaHGgKWEG6/tlKo1hugqMpLUZWHokoPhyo8VLhPT4w0GxUG93EwJMPB0HQHeZkO+jhO/79ai8lAvNVIvMVEvNVEnMWIzdwDcjTCIS1Q8ddpj4C7+UqhU4OVqkJwlcCpS3EVAyTnaoFKekPA4sg8cw6LomjBisXR8IjXns/hZNtG/rAfd8CNN+jFHXTjCXpaXN7sC/k45j5GsauY4rpiDjsPU+IuOa1SroJCjiOHQUmDGJw0mMFJg+mX0K/dy5stRgvxpnjizHHEmbSgpaeuDuoq5/5voxCixwtHVDyBEB6/9nD7w3j8IU64fBRVejlSrQUrR6q81Na3PFzfN9nO0AwHQ9IdDM1w0D/VjqnJ6ICigM1sJN5iJN5qIt5iIs5q7N6VQY1UVVvp0xis+OtO5q+oqlabpaao4XG49WAFtA0M04Zqj/Rh0OcCLQm3NQaj9nljoGJ1gDn+nB9ZAa0OS9NAxRv0tjgd5PQ7OVp3lCOuIxytO0pxXTEnPCdaDGySrEnRQGVw0mAGJg1s16oeBQW7yY7dZNeClYaA5XwaWWkr+YkIIbrUqcFKnU8LVI5W13OsxsuxmnqO1ngprva2uN+QAmQn28hNi2dgWjyD+sQzJD2+Wbl9i8lAnMVInMWI3WIkzmLCbjZ2f5ItaHVVAh5tRCX67NWWLwfrofZIk2Cl4RFoZSVLfLqWYNsn72TQYktq+VzFAGabFpxY4k4+m2zn5B5CpwqGg3hDXrxBL96QF0/Qc1ouii/k47j7OCWeEu3ZXcLRuqOtbpyYZE1iQMIABiQOIDchl8HJg0m1pba5T1ajNRqsNH2cqzkrepMARggRM75gGG8gjLchYCl1+iiscHOsWgtSjtXUc6ymHre/5bLpZqPCgNQ4BqbFk5umPfdPjYtO8VhMCnazqSFIaQhWzD0kXwW0VUFBT8PoSkPAEqrXSvU7j2rLl53HtNe1R8Fd1nI7igES+0FKLqQM1IKWtKFa8beWzjXbtFEVs71JwBJ3XgQqqqriC/uigYo36KU+VN9sZMUf8nPcowUox93Ho8FKS/kqoI2KZMZnasFKQ8DSP6E/SdZWgsVT2Iw2bCbt0TRQ6cp9g85FEsAIITotGI7gDYSpD4Sp9Pg5VO7mUKWHkpp6ypw+jju155YSa0H7u5qZYKNfip1+KXH0S7EzIDWOnGQ7ZqOCvSEvxW5ueLYYsZkMPSdQCQe14CTo1UZRAh4tf8VZAq7jJ4MU5zFwHdMScFtjT20IVAadfE7qd3qSrdHSEKA0BioNj/NkRAW0fJX6YL0WsIS0QKU+WI+KSigSoqq+ijJvGSc8J5o9n1pnpalESyJ9HX3JceTQ19GXfgn96O/of9ZVQSaDKRqo2E12rEarFrQYW1+KLzpHAhghRJuoqoo/FKE+EKaizs+hSjdHqrwcq/FyvNZHqdNHmctHtSfQahsKkJ5gjQYp/VLs9E+No1+ynUS7GavZgM2kBSk2s6HhuQck1DYK+RuCFJ/2XF8D1Ye1aR9XKdSVQt1xqCvTVv+oLe9eDWjTP0n9tEdif+05Jbf5FJDJqgUkJps2qmKyn3w+DxJqoeH3LuynPlSPL+TTgpSw9tof8lPpq6TCW0FFfQUV3gpOeE9Q5imjsr7yjPVYEi2J0SAlx5FDTrz2+kxVbc0GM1ajFavJis1oiwYpFqPlnKmt0pucH/8GCCHazB8K46wPUlzt5WiVl6IqL0drvByvrafc5ae8zt/qlE8jh9VETrKN7CQ7WUk2cpLsDEi1MyAtjkSbGZvZqAUrZiNWkwGrqQcFKeGQFpyEfNpISc0RLUBxHj0ZmLjLtWdPJS0m0jYyWiAh+2SgktQPkvprGx6a7VoCbdMgxWRrGEWxakHKeZBIC1qQEogE8IV8WmAS9kengU54T1BVXxUNUCrrK6Ova/ytj6SAVl8lIy6DrPgssuKzyIzLjD63FKgYFANWoxWL0aIFKqc8esomhkIjAYwQ55lQKEJZnY+iSg/F1fWU1Hgbpnj8nKjzUekO4GplpU9TiXYzmQlW0hOsZCZq0z8D0+IY2CeejARbNDCxmg1YjIZWC8R1OVXVgpNgvRaQ1B5pmNop0UZRGgMUT7m2+udMoyigBRwJOVqgkpitPUcfmU2CE+vpr8+TURTQKtD6w/7oIxAOUOOrocRdQqmnlCpfFdX11dqzr5qq+ipq/DWnLUs+ldVoJd2eTnpcOn3sfciKyyIzPpOsuCySbcnN8kzMBjNmoxmLwXJaoGI2mmUUpZc5f/7tEeI84PGHKKnxUlJbT0mNlntywuXTRk7cfqrcfmq8wdN2Xm6J1WSgj8NKZqI1OorSP1Wb8slNiyPZbsFiMmgPo6HnzPNHwtrUTmOCbF2pFqg0PrvLwVupjZ6coYJqlNEC8RngaPrI1EZRUnIhPlOb1jFatMAk+mw9b0ZQQNu0MBAOEAgH8IQ8lLpLKfOUccJ7gsr6Sqrqq6j111Ljr9GefTVn3dgQtFGRFGsK6XHppNu1IKXxdXpcOgnmBBRFaRacWIwWLAYLZqM5GqiYDWZJmj3HSAAjRA/nC4apdPupdAcoa0iGPeHyU94wWlLl9lPlCVDjDbS47LglCpASbyHdYSUj0UpWoo3sJBs5ySfzUtLiLVhMxu7fkLBRyN8kECmFuhOnTOdUaIGJt7rlPYBaY0uG+D4Q10d7dmRCUl9tqic5V3tvsjYEJWYtQDFazptEWdDqpdQH66msr+SE9wTl3nIq6iuoqq+KjpQ07rbs9DvPuoFhU/HmeNJsaaTaUkmzNzzb0ki1p9LH1oc+cX20ERKDNkLSGIw0fW0ymHpOAC26jAQwQnShcESlzhekxhuk1hug2qMFIOVuP5V12utqT4Bqb4Bab5Bab5D6YOuJiC2xm42kOSz0cWjTOxkNUzzZSTaykmz0TbaTnWTDbunGf/1VVVul463WHo0BiKdce++pbAhGqrRHfY22sqc9jFawp0BcKsSlaUmzjdM8iQ0BSmKOVrzNaGkITs6PKQRVVfEEPVoA4tOCkEqfNkpS46vRRkl8tTgDTlx+F3XBOjzB9v38TYqJZFsyyVbtkWJLIdmaTJotjT72PvSx94nmojRO35gMpmhwYjaYJedEnJEEMEK0k6qq1AfDDRsIhqIbCTrrg9rGgt4g1V5tRKTGE6S2XgtGXPVB6nyhM6V8tspkVEiJs5AW3/BoCE4aR1AyE61kJdrJTLLhsHbhv9bRzQRdWpDhq4X62pMjId4q7b2vpuHZqZ3nd2pTPe2lGLXaJ/YUbblxXNrJURNH5skclMS+2uemc3ukJBAORHdAdvqd1PprqfXX4vQ7o+8bX7sCrmYbEJ5phU5rFBQSLAkkWhJJsiZFA5MUawqp9lRSbana1I49nVRbKmajFpREH4qMlAj9SAAjzhuNy4C9gXBDufpQ9PnknjtB3L4QrobXLl+Iunpt48C6JsfakkNyJnEWI0l2M0l2sxaYOCykxlvoE2+lT4I2epLmsJIWbyHVYSHBqvN/+MMhCDSWrXc37LXTpIx94zGfUwtKfLVakOJzNuzP49LOCbe+ZLpNTDZtCseerD3HpZ4MThzpJwMTRxYkZmnnnAN/AEORULR0vSfowRPyaM9NNg1sDEzqAnXRoMMT9OAOnnwd6OTP32wwE2+OJ8GSQJIlKRqUJNuSSbYkk2JPaTa9k2JNwWayYVSMMjoiup0EMKLHCUdUfMEw9UGtMFqz54bXjdVdvYEwXr/27AmEqQ+E8EQrvzYPVLyBMKFOBh5NGQ0KCVYTiXYzCTYTiTYzSXFmkmwmUuItpMRpj+Q4M8lxFlIanpPs5rbllURXy9SBy6uVmw96Gp7rm7z2NBRP856s+tp0X52mZeuDXu116OzJk+1icWj1S2yJDc8pTUZKmoyWxDV9naqt4OmhGuuP+EK+aKG0xmW+jfVIosXTQvXRiq/uoPtkyfom1WAbz/UGvQQinQz8mlBQmu1I7DA7cFgcpwUlSdak6EhJsjWZNHuaFFkTvVqHAphnnnmGhx9+mLKyMsaNG8dTTz3F5MmTWz3/zTff5Fe/+hVFRUXk5eXxu9/9jiuvvLLDnRaxF4moBMIR/KEIgVAEfyhMIBQhENbeB6LHtc/8oQj+YJPXoQj+oPbaFwzjC0bwhcINgYl2zN8QkPga34ci1Ae174k1m9kQ3X043mrEYdUCkATbyYDEYTWTYDWQZFFJtKgkmCIkmSMkmMM4jBHsShAlHNCCjLAPQrVaYBDya8dCfvDVg9uvBRghv1ZGPtTkfWO9kYBX+yzY8Aj5tEesGS0NOw0naJv4WRxgTWh4JGrP9pTmwYk1sfl7S0LMVtuoqkowEowuuw1GgtGVLoFIIHqs6dLcps+nvm6sMRINTBqCk+jy3pBWf6TxeKyZDCZt0z5THHaTHYfFEQ1EEiwJJ4MRc0I0IGkMThKtidFzZHWNOB+1O4B54403WLJkCc899xxTpkzh8ccfZ+7cuRQUFJCRkXHa+Vu3buW73/0uq1at4uqrr+a1115jwYIF7Nq1i9GjR+tyEz1JJKISVlXCEZVQRCUcVglFIoQbjofCDccjEYLhk+eFwpGG4yrBcKThWbs2FNaONZ7XeDwYbmxPCywazws2nBNsOBaIHosQDJ18HwidPDfQ5H0gFNFhpELFSAQTYQxEMBPGSBgTEe1ZaXgmjJ0wCZx8b1bCGJUwZsLEGVXspkjDs4rdqGI3RqLHbYYIdkMYq0F7bTWEsSgRLIrWjlUJYyYUfRgJYoqEUCJBrfx7OKA9ggGo9zccawhCwsG2LbPtCgazNlrRuAmfOe7k/jaWeO3RuO+N1YFqjiNsiSNsjiNsthMy2wib7ITNNsImKyGTlbDZRshgIBQJEY6ECathQpGQ9r7hdVgNE4wEo8dDkQChQBkh3zFCqnas6eeNr8/0HAwHteeGYKTxdTAcJBDR3ofCoejrnsBsMEfLwltNVuxGOzaTTQs8zFoAEmeKi+4e3DQoiTPHEW+Kj34Wb44n3hRPvDke83mSNCxELCiqqrbrL9WUKVOYNGkSTz/9NACRSIT+/ftz5513cu+99552/o033ojH42HNmjXRY1OnTmX8+PE899xzbfpOl8tFUlISj7y2BpvNjqqqqKpKJBJBVSOoKlqxIzVCJKISUVXUhs8iEVU7J/oc0T5XT35Ow7HGNtRou2qT183baPycpueqERQ1ggIYiKCgYkBteA0Gpekx9bRzDKgN7yPRc5SGQEBRwi2cG2lybgRFiTS0qR3TrlOj7xvPMUav0x6NwYZBUaNtaI/wyT40XqdEMCpqw700thuO3otCGKWhD9qVWp1SVSH6WhtfOflZRDlZyzTScDyiNL5vvEYB5eT7k9cqRDh5bgRQFSX6PRHtqxp6dPL7ot9zSjsnP6fhJ3yyfxEUwgYDqsFExGAkYjARUYwNrxseiiH6HDYYURueI4qBsKJo7xUDEUUhrCjaM43P2neFG3564Yb+hdF+N8NqmHAkTESNEFJD2rGGwCN6LNJwXgcSNHsDk8F0ss5HwxLaU6ulNq390Xis6evGIKSxFPypZeEbdwhuuqeNySCz7UJ0ROPfb6fTSWJioq5tt+vfykAgwM6dO1m2bFn0mMFgYNasWWzbtq3Fa7Zt28aSJUuaHZs7dy7vvPNOq9/j9/vx+0/O0Tud2sZn/zh6F8Y4LXGsadSlnvJMC59pr5XTjynAKbloZ2r75HPzeWNVOf37W2sn2g/lTOdof4h7FkPDQ2hUUENaQmwvihdMigmjwaglYirG6OoQo8EYPd64aqTxddNjZsWMwWBotrKk6eumK0+aLodtfN14TnS5bNPXTZbSRut9GM1dX4hMBYIQDobx0o6aMkKIZlwuF6BNB+utXQFMZWUl4XCYzMzMZsczMzPZv39/i9eUlZW1eH5ZWSvbxgOrVq1ixYoVpx3f/PND7emuEEIIIXqAqqoqkpKSzn5iO/TIcdFly5Y1G7Wpra0lNzeX4uJi3X8APZnL5aJ///4cPXpU96G3nkzuW+77fCD3Lfd9PnA6nQwYMIDU1FTd225XANOnTx+MRiMnTpxodvzEiRNkZWW1eE1WVla7zgewWq1YrdbTjiclJZ1X/+AbJSYmyn2fR+S+zy9y3+eX8/W+DTFYqdiuFi0WCxMnTmT9+vXRY5FIhPXr1zNt2rQWr5k2bVqz8wHWrVvX6vlCCCGEEGfT7imkJUuWcOutt3LhhRcyefJkHn/8cTweD9///vcBuOWWW+jbty+rVq0C4Mc//jGXXHIJq1ev5qqrruL111/n888/53//93/1vRMhhBBCnDfaHcDceOONVFRUsHz5csrKyhg/fjxr166NJuoWFxc3GyqaPn06r732Gr/85S/5xS9+QV5eHu+88067asBYrVbuv//+FqeVzmVy33Lf5wO5b7nv84Hct/733e46MEIIIYQQ3U2KegghhBCi15EARgghhBC9jgQwQgghhOh1JIARQgghRK/T4wOYZ555hoEDB2Kz2ZgyZQo7duzo7i7patWqVUyaNImEhAQyMjJYsGABBQUFzc6ZOXMmiqI0e9xxxx3d1GN9PPDAA6fd0/Dhw6Of+3w+Fi1aRFpaGg6Hg+uvv/60goi90cCBA0+7b0VRWLRoEXDu/LPevHkz8+fPJycnB0VRTtv7TFVVli9fTnZ2Nna7nVmzZvHNN980O6e6upqbb76ZxMREkpOT+eEPf4jb7e7Cu2i/M913MBjk5z//OWPGjCE+Pp6cnBxuueUWjh8/3qyNln5HHnrooS6+k/Y52z/vhQsXnnZP8+bNa3bOufbPG2jx33VFUXj44Yej5/TGf95t+bvVlv+GFxcXc9VVVxEXF0dGRgY//elPCYVCbe5Hjw5g3njjDZYsWcL999/Prl27GDduHHPnzqW8vLy7u6abTZs2sWjRIj799FPWrVtHMBhkzpw5eDyeZufdfvvtlJaWRh+///3vu6nH+hk1alSze/rkk0+in91zzz289957vPnmm2zatInjx49z3XXXdWNv9fHZZ581u+d169YB8J3vfCd6zrnwz9rj8TBu3DieeeaZFj///e9/z5NPPslzzz3H9u3biY+PZ+7cufh8vug5N998M3v37mXdunWsWbOGzZs386Mf/airbqFDznTfXq+XXbt28atf/Ypdu3bx1ltvUVBQwDXXXHPauStXrmz2O3DnnXd2Rfc77Gz/vAHmzZvX7J7++te/Nvv8XPvnDTS739LSUv70pz+hKArXX399s/N62z/vtvzdOtt/w8PhMFdddRWBQICtW7fy8ssv89JLL7F8+fK2d0TtwSZPnqwuWrQo+j4cDqs5OTnqqlWrurFXsVVeXq4C6qZNm6LHLrnkEvXHP/5x93UqBv7/7d1tSFPvGwfwrzPnA6VLp04TRc2s8AE1GiPyRYopwU/qRWZC9lyWVGQhBhH1ogRBX/RCemEaGEVBJvRIPkHlshSHmTV0TKVwRsamoebDrt+L/3+H/5rp/qSbZ10fEOZ9n7Ndt9c59325nW2XLl2ixMTEOfuMRiN5eHjQ/fv3hbaPHz8SAFKr1Q6K0DFOnz5N0dHRZDabicg1cw2A6urqhN/NZjMpFAoqKysT2oxGI3l6etKdO3eIiKinp4cA0Lt374Rtnj59Sm5ubvTlyxeHxf4nfh33XN6+fUsAaGBgQGiLiIigioqKpQ1uCc017vz8fMrOzv7tPn9LvrOzs2nbtm1WbWLPN5HtumXPHP7kyROSSCRkMBiEbSorK8nX15d+/vxp1+Mu22dgpqam0NHRgfT0dKFNIpEgPT0darXaiZEtLZPJBAA2X3x1+/ZtyOVyxMXFoaSkBOPj484Ib1H19vYiNDQUUVFRyMvLw+DgIACgo6MD09PTVrlfv349wsPDXSr3U1NTqK2txcGDB+Hm5ia0u2Ku/5der4fBYLDKr5+fH5RKpZBftVoNmUyGTZs2Cdukp6dDIpGgra3N4TEvFZPJBDc3N8hkMqv20tJSBAQEICkpCWVlZf/X0+rLVUtLC4KCghAbG4uCggKMjIwIfX9DvoeHh/H48WMcOnTIpk/s+f513bJnDler1YiPjxc+BBcAtm/fjtHRUXz48MGux12W30YNAN++fcPs7KzV4AAgODgYnz59clJUS8tsNuPMmTPYsmWL1ScV7927FxEREQgNDUVXVxeKi4uh1Wrx4MEDJ0b7Z5RKJWpqahAbG4uhoSFcvnwZW7duRXd3NwwGA6RSqc2kHhwcDIPB4JyAl8DDhw9hNBqxf/9+oc0Vc/0rSw7nOrctfQaDAUFBQVb9K1asgL+/v8scA5OTkyguLkZubq7Vl/udOnUKycnJ8Pf3R2trK0pKSjA0NITy8nInRvtnMjMzsWvXLkRGRkKn0+HChQvIysqCWq2Gu7v7X5HvW7duYdWqVTYvhYs933OtW/bM4QaDYc45wNJnj2VbwPyNTp48ie7ubqtrQQBYvQ4cHx+PkJAQpKWlQafTITo62tFhLoqsrCzhdkJCApRKJSIiInDv3j14e3s7MTLHqaqqQlZWFkJDQ4U2V8w1szU9PY3du3eDiFBZWWnVd/bsWeF2QkICpFIpjh07hmvXron2Y+j37Nkj3I6Pj0dCQgKio6PR0tKCtLQ0J0bmODdv3kReXh68vLys2sWe79+tW46wbF9CksvlcHd3t7lqeXh4GAqFwklRLZ3CwkI8evQIzc3NCAsLm3dbpVIJAOjr63NEaA4hk8mwbt069PX1QaFQYGpqCkaj0WobV8r9wMAAGhoacPjw4Xm3c8VcW3I437mtUChsLtafmZnB9+/fRX8MWIqXgYEBvHjxwurZl7kolUrMzMygv7/fMQE6QFRUFORyuXBcu3K+AeDly5fQarULnu+AuPL9u3XLnjlcoVDMOQdY+uyxbAsYqVSKlJQUNDY2Cm1msxmNjY1QqVROjGxxEREKCwtRV1eHpqYmREZGLriPRqMBAISEhCxxdI7z48cP6HQ6hISEICUlBR4eHla512q1GBwcdJncV1dXIygoCDt27Jh3O1fMdWRkJBQKhVV+R0dH0dbWJuRXpVLBaDSio6ND2KapqQlms1ko6sTIUrz09vaioaEBAQEBC+6j0WggkUhsXmIRs8+fP2NkZEQ4rl013xZVVVVISUlBYmLigtuKId8LrVv2zOEqlQrv37+3KlwtBf3GjRvtDmTZunv3Lnl6elJNTQ319PTQ0aNHSSaTWV21LHYFBQXk5+dHLS0tNDQ0JPyMj48TEVFfXx9duXKF2tvbSa/XU319PUVFRVFqaqqTI/8zRUVF1NLSQnq9nl6/fk3p6ekkl8vp69evRER0/PhxCg8Pp6amJmpvbyeVSkUqlcrJUS+O2dlZCg8Pp+LiYqt2V8r12NgYdXZ2UmdnJwGg8vJy6uzsFN5tU1paSjKZjOrr66mrq4uys7MpMjKSJiYmhPvIzMykpKQkamtro1evXlFMTAzl5uY6a0h2mW/cU1NT9M8//1BYWBhpNBqr893yrovW1laqqKggjUZDOp2OamtrKTAwkPbt2+fkkc1vvnGPjY3RuXPnSK1Wk16vp4aGBkpOTqaYmBianJwU7sPV8m1hMpnIx8eHKisrbfYXa74XWreIFp7DZ2ZmKC4ujjIyMkij0dCzZ88oMDCQSkpK7I5jWRcwRETXr1+n8PBwkkqltHnzZnrz5o2zQ1pUAOb8qa6uJiKiwcFBSk1NJX9/f/L09KS1a9fS+fPnyWQyOTfwP5STk0MhISEklUppzZo1lJOTQ319fUL/xMQEnThxglavXk0+Pj60c+dOGhoacmLEi+f58+cEgLRarVW7K+W6ubl5zuM6Pz+fiP7zVuqLFy9ScHAweXp6Ulpams3fY2RkhHJzc2nlypXk6+tLBw4coLGxMSeMxn7zjVuv1//2fG9ubiYioo6ODlIqleTn50deXl60YcMGunr1qtVCvxzNN+7x8XHKyMigwMBA8vDwoIiICDpy5IjNP6Kulm+LGzdukLe3NxmNRpv9xZrvhdYtIvvm8P7+fsrKyiJvb2+Sy+VUVFRE09PTdsfh9t9gGGOMMcZEY9leA8MYY4wx9jtcwDDGGGNMdLiAYYwxxpjocAHDGGOMMdHhAoYxxhhjosMFDGOMMcZEhwsYxhhjjIkOFzCMMcYYEx0uYBhjjDEmOlzAMMYYY0x0uIBhjDHGmOhwAcMYY4wx0fkX2bd3E169n0YAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "\n", "i, j, k = 12, 1073, 7999\n", "\n", "timeline = np.linspace(0, 200)\n", "regression.plot(\"hf\", timeline, covar[i], label=f\"hf of the {i}-th asset\")\n", "regression.plot(\"hf\", timeline, covar[j], label=f\"hf of the {i}-th asset\")\n", "regression.plot(\"hf\", timeline, covar[k], label=f\"hf of the {i}-th asset\")\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 }