{ "cells": [ { "cell_type": "markdown", "id": "7b37cd74-dc44-405f-aaf0-52fe31a38c13", "metadata": {}, "source": [ "# Plot GRN" ] }, { "cell_type": "code", "execution_count": 1, "id": "4a4f718d-d97e-4227-92c0-c9d5fa62d4b4", "metadata": {}, "outputs": [], "source": [ "from grnet.clusters import CellClasses\n", "from grnet.models import PC\n", "from grnet.plot import grnplot\n", "from grnet.toydata import load_dataset, load_metadata\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "id": "6d8e71df-63a7-4d4f-a104-de44a3a72f15", "metadata": {}, "source": [ "## data loading" ] }, { "cell_type": "code", "execution_count": 2, "id": "3087432e-7a3f-44a5-8de4-b7d0aa5296e0", "metadata": {}, "outputs": [], "source": [ "data = load_dataset(\"prototype1\")\n", "meta = load_metadata(\"prototype1\")\n", "\n", "lst_data = [\n", " data.loc[meta.cluster == v, :].iloc[:, :5] for v in meta.cluster.unique()\n", "]" ] }, { "cell_type": "markdown", "id": "8d497e05-4ff7-4146-854f-f2090bf28182", "metadata": {}, "source": [ "## GRN formation with PC algorithm" ] }, { "cell_type": "code", "execution_count": 3, "id": "2d864ac2-b46f-4a6c-afc2-23fa8712cf2a", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "INFO:pgmpy: Datatype (N=numerical, C=Categorical Unordered, O=Categorical Ordered) inferred from data: \n", " {'gene_1': 'N', 'gene_2': 'N', 'gene_3': 'N', 'gene_4': 'N', 'gene_5': 'N'}\n", "INFO:pgmpy: Datatype (N=numerical, C=Categorical Unordered, O=Categorical Ordered) inferred from data: \n", " {'gene_1': 'N', 'gene_2': 'N', 'gene_3': 'N', 'gene_4': 'N', 'gene_5': 'N'}\n", "INFO:pgmpy: Datatype (N=numerical, C=Categorical Unordered, O=Categorical Ordered) inferred from data: \n", " {'gene_1': 'N', 'gene_2': 'N', 'gene_3': 'N', 'gene_4': 'N', 'gene_5': 'N'}\n", "INFO:pgmpy: Datatype (N=numerical, C=Categorical Unordered, O=Categorical Ordered) inferred from data: \n", " {'gene_1': 'N', 'gene_2': 'N', 'gene_3': 'N', 'gene_4': 'N', 'gene_5': 'N'}\n", "INFO:pgmpy: Datatype (N=numerical, C=Categorical Unordered, O=Categorical Ordered) inferred from data: \n", " {'gene_1': 'N', 'gene_2': 'N', 'gene_3': 'N', 'gene_4': 'N', 'gene_5': 'N'}\n" ] } ], "source": [ "models = [PC(v) for v in lst_data]\n", "\n", "for v in models:\n", " v.estimate()" ] }, { "cell_type": "markdown", "id": "df112136-4e11-4090-adca-92f06ac99a79", "metadata": {}, "source": [ "## Plotting GRN" ] }, { "cell_type": "code", "execution_count": 4, "id": "90fd76a4-9e6c-430a-a712-736caf1d3e94", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVMAAAFeCAYAAAA41fyqAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAMKhJREFUeJzt3XlYVPXiBvB3ZhiGRUDZERNRUJAdScENyyitDDU1t7xl2dVbarstppjZrcwbZd4W7/2VZeGeZZa3UMkQE1MQAUFUBkVRWQTZl5nv7w+vc0UGRBw4w/B+nsfncc585/CemeHlnDnLyIQQAkREdFvkUgcgIjIFLFMiIgNgmRIRGQDLlIjIAFimREQGwDIlIjIAlikRkQGwTImIDIBlSkRkACxTA0lISIBMJkNCQoLUUYza3/72N0RFRUkdo8uaOnUqpkyZInUMk8Qy7SSSkpIQExOD0tJSSX5+RUUFli5dijFjxsDe3h4ymQxffvnlLc0jNzcX//rXv/Daa6+1T8ibqKqqQkxMzG39wZP6dWjOihUr8NBDD8HFxQUymQwxMTF6xy1atAhbt27F0aNHOzZgF8Ay7SSSkpKwbNkyyX6Ji4qK8Oabb+L48eMICgpq0zw+/PBDeHp64q677jJwutapqqrCsmXLbrtMpXwdmrN48WIcOnQIISEhLY4LCQlBWFgYVq1a1UHJug6WaRdXVVXVqnFubm4oKChAXl4eVq5cecs/p76+Ht98843BNzFbm9/U5ebmoqCgAOvXr7/p2ClTpmDbtm2oqKjogGRdB8u0lc6dO4cnnngCPXv2hEqlgqenJ+bNm4e6urpmH9OnTx889thjTaaPGjUKo0aNajRt9erV8PPzg5WVFXr06IGwsDB8++23AICYmBi89NJLAABPT0/IZDLIZDKo1Wrd49evX49BgwbB0tIS9vb2mDp1Ks6ePdvk5/r7++Pw4cMYOXIkrKysWr3JrVKp4Orq2qqx+iQmJqKoqAj33HNPo+nXPmveuHEjXnvtNbi6usLa2hoPPfTQLeW/dOkSnnjiCbi4uMDCwgJBQUFYt26d7rFqtRpOTk4AgGXLlumew+s3h/fs2YMRI0bA2toa3bt3R3R0NI4fP667v6XXITIystk19gEDBuC+++7T5ZDJZHj//ffxwQcfwMPDA5aWloiMjER6enqTx2ZlZWHSpEmwt7eHhYUFwsLC8MMPPzQZ16dPn+ae+iaioqJQWVmJX3/9tdWPoZszkzpAZ3D+/HkMHjwYpaWleOqpp+Dj44Nz585hy5YtqKqqgrm5+W3Nf+3atViwYAEmTZqEhQsXoqamBmlpaTh48CCmT5+OiRMn4sSJE4iLi8MHH3wAR0dHANCVw4oVK/DGG29gypQpePLJJ1FYWIjVq1dj5MiRSElJQffu3XU/q7i4GGPHjsXUqVMxc+ZMuLi43Fb21kpKSoJMJmt2M3TFihWQyWRYtGgRLl26hNjYWNxzzz1ITU2FpaVli/mrq6sxatQonDx5Es888ww8PT2xefNmPPbYYygtLcXChQvh5OSETz75BPPmzcOECRMwceJEAEBgYCAAID4+HmPHjkXfvn0RExOD6upqrF69GsOGDcORI0fQp0+fFl+HRx99FHPmzEF6ejr8/f11eQ8dOoQTJ05g8eLFjZb3q6++Qnl5OZ5++mnU1NTgww8/xN13341jx47pXpOMjAwMGzYM7u7ueOWVV2BtbY1NmzZh/Pjx2Lp1KyZMmNCm12LgwIGwtLTE/v372zwP0kPQTc2aNUvI5XJx6NChJvdptVohhBB79+4VAMTevXt193l4eIi//OUvTR4TGRkpIiMjdbejo6OFn59fixlWrlwpAIjc3NxG09VqtVAoFGLFihWNph87dkyYmZk1mh4ZGSkAiE8//bTFn3Uzhw4dEgDEF1980erHzJw5Uzg4ODSZfu15c3d3F1euXNFN37RpkwAgPvzww5vmj42NFQDE+vXrddPq6upERESE6Natm26+hYWFAoBYunRpkxzBwcHC2dlZFBcX66YdPXpUyOVyMWvWLN205l6H0tJSYWFhIRYtWtRo+oIFC4S1tbWoqKgQQgiRm5srAAhLS0uRn5+vG3fw4EEBQDz33HO6aaNHjxYBAQGipqZGN02r1YqhQ4cKb2/vJstws2W8Xv/+/cXYsWNbHEO3hpv5N6HVarF9+3aMGzcOYWFhTe6XyWS3/TO6d++O/Px8HDp06JYfu23bNmi1WkyZMgVFRUW6f66urvD29sbevXsbjVepVHj88cdvO/OtKi4uRo8ePZq9f9asWbCxsdHdnjRpEtzc3PDTTz81Gqcv/08//QRXV1dMmzZNN02pVGLBggWoqKjAb7/91mK2goICpKam4rHHHoO9vb1uemBgIKKioppk0MfOzg7R0dGIi4uD+O/11jUaDTZu3Ijx48fD2tq60fjx48fD3d1dd3vw4MEYMmSI7meVlJRgz549mDJlCsrLy3Wva3FxMe677z7k5OTg3LlzN83VnB49eqCoqKjNj6emWKY3UVhYiCtXrjTadDO0RYsWoVu3bhg8eDC8vb3x9NNPY//+/a16bE5ODoQQ8Pb2hpOTU6N/x48fx6VLlxqNd3d3v+2PJdpKtPClDt7e3o1uy2QyeHl5NfpcGNCfPy8vD97e3pDLG7+dfX19dfe35Nr9AwYMaHKfr68vioqKUFlZ2eI8gKt/EM6cOYPff/8dwNWPDi5evIhHH320ydgblxcA+vfvr1vekydPQgiBN954o8nrunTpUgBo8treCiGEQVYE6H/4mWk7au7NqtFooFAodLd9fX2RnZ2NH3/8Ebt27cLWrVvxz3/+E0uWLMGyZcta/BlarRYymQw///xzo3le061bt0a3r//8sSM5ODjg8uXLtz0fqfK3xn333QcXFxesX78eI0eOxPr16+Hq6tpkp1traLVaAMCLL76o23l1Iy8vrzZnvXz5st5Cp7Zjmd6Ek5MTbG1t9e5pvZkePXroPR4xLy8Pffv2bTTN2toajzzyCB555BHU1dVh4sSJWLFiBV599VVYWFg0W8z9+vWDEAKenp7o37//LWfsKD4+Pvjmm29QVlYGOzu7Jvfn5OQ0ui2EwMmTJ3U7iFri4eGBtLQ0aLXaRmunWVlZuvuB5v+4Xbs/Ozu7yX1ZWVlwdHTUbaa3tDanUCgwffp0fPnll3j33Xexfft2zJkzR+8fuRuXFwBOnDih2yt/7f2hVCrbVMYtaWhowNmzZ/HQQw8ZdL5dHTfzb0Iul2P8+PHYsWMH/vzzzyb3t7Tp2q9fP/zxxx+NDp/68ccfmxzyU1xc3Oi2ubk5Bg4cCCEE6uvrAUD3y3xjOU+cOBEKhQLLli1rkkUI0WTeUomIiIAQAocPH9Z7/7W929ds2bIFBQUFGDt27E3nff/99+PChQvYuHGjblpDQwNWr16Nbt26ITIyEgBgZWUFoOlz6ObmhuDgYKxbt67Rfenp6fjll19w//3366Y19zpc8+ijj+Ly5cv461//ioqKCsycOVPvuO3btzf6zDM5ORkHDx7ULa+zszNGjRqFzz77DAUFBU0eX1hY2MyzcXOZmZmoqanB0KFD2zwPaoprpq3w9ttv45dffkFkZCSeeuop+Pr6oqCgAJs3b0ZiYmKjQ4+u9+STT2LLli0YM2YMpkyZglOnTmH9+vXo169fo3H33nsvXF1dMWzYMLi4uOD48eP4+OOP8cADD+h2ygwaNAgA8Prrr2Pq1KlQKpUYN24c+vXrh7feeguvvvoq1Go1xo8fDxsbG+Tm5uK7777DU089hRdffNEgz8PHH3+M0tJSnD9/HgCwY8cO5OfnAwDmz5+vd43zmuHDh8PBwQHx8fG4++67m9xvb2+P4cOH4/HHH8fFixcRGxsLLy8vzJkz56a5nnrqKXz22Wd47LHHcPjwYfTp0wdbtmzB/v37ERsbq3sOLS0tMXDgQGzcuBH9+/eHvb09/P394e/vj5UrV2Ls2LGIiIjAE088oTs0ys7OrtGxqM29DtdKNiQkBP7+/ti8eTN8fX0RGhqqN7OXlxeGDx+OefPmoba2FrGxsXBwcMDLL7+sG7NmzRoMHz4cAQEBmDNnDvr27YuLFy/iwIEDyM/Pb3RK6Ndff428vDzdSQz79u3DW2+9BeBqwV9b+waAX3/9FVZWVrxGgqFJcxBB55OXlydmzZolnJychEqlEn379hVPP/20qK2tFULoPzRKCCFWrVol3N3dhUqlEsOGDRN//vlnk0OjPvvsMzFy5Ejh4OAgVCqV6Nevn3jppZdEWVlZo3ktX75cuLu7C7lc3uTwnK1bt4rhw4cLa2trYW1tLXx8fMTTTz8tsrOzdWMiIyNveghWSzw8PAQAvf9uPFRInwULFggvL69G0649b3FxceLVV18Vzs7OwtLSUjzwwAMiLy+v0diW8l+8eFE8/vjjwtHRUZibm4uAgAC9h24lJSWJQYMGCXNz8yaHEMXHx4thw4YJS0tLYWtrK8aNGycyMzObzKOl10EIId577z0BQLz99ttNHnvt0KiVK1eKVatWiTvuuEOoVCoxYsQIcfTo0SbjT506JWbNmiVcXV2FUqkU7u7u4sEHHxRbtmxp8tw099rc+J4cMmSImDlzpt7nkdqOZUod5tSpU0KpVIr4+HjdtGtlunnzZgmTGVZsbKyQyWRN/hgI0bhMpZCSkiJkMplISUmR5OebMn5mSh2mb9++eOKJJ/DOO+9IHaXdCCHw73//G5GRkejdu7fUcZp45513MGnSJAQHB0sdxeTwM9MurqKi4qYXvHByctK7R7otPvnkE4PMx9hUVlbihx9+wN69e3Hs2DF8//33UkfSa8OGDVJHMFks0y7u/fffv+mxrLm5ubd0IY2uqLCwENOnT0f37t3x2muv8bCjLkgmRAvH9pDJO336NE6fPt3imOHDh8PCwqKDEhF1TixTIiID4A4oIiIDYJkSERkAy5SIyABYpkREBsAyJSIyAJYpEZEBsEyJiAyAZUpEZAAsUyIiA2CZEhEZAMuUiMgAWKZERAbAMiUiMgCWKRGRAbBMiYgMgGVKRGQALFMiIgNgmRIRGQDLlIjIAFimREQGwDIlktD+/fsRFhaGgIAARERE4NixY1JHojbit5MSSSgjIwN2dnbo1asX4uPjsWLFCuzdu1fqWNQGXDMlusHixYsxYMAAREVFISoqCgkJCdi5cyfCw8MRHByMuXPnQqvVQq1WIzQ0FDNmzICPjw/mzZunm4e+8fr4+fmhV69eAIBBgwbh3LlzHbKMZHgsU6LrJCcnY8+ePUhPT8e6detw8OBBFBUVITY2FgkJCUhNTYVCocC2bdsAAJmZmYiJiUFGRgYOHDiAnJycFse35KuvvsLo0aPbexGpnZhJHYDImCQlJWHChAlQKpXo2bMnRowYAZVKhbS0NISHhwMAqqur4eHhgbCwMPj6+sLb2xsAEBgYCLVajaysLL3jW3L48GF8+umnSExMbN8FpHbDMiW6jr5dCEIIREdH4/PPP280Xa1WQ6VS6W7L5XJoNJpmxzcnPz8f06dPx6ZNm+Dg4HB7C0CS4WY+0XWGDh2K7du3o6GhAQUFBUhMTERNTQ12796N/Px8AEBxcbHu//qEh4e3enxlZSWio6OxatUqBAUFGX6BqMOwTImuM2TIEERGRsLPzw9/+ctfEBQUBC8vL6xZswbR0dEIDAzEvffei0uXLjU7D2dn51aP/+KLL5CVlYXFixcjODgYQ4cOba9Fo3bGQ6OIblBZWQlra2sUFRUhIiICR44cgY2NjdSxyMjxM1OiG8yePRvZ2dmor6/H8uXLdUWq0Qok55bgUnkNnG0sMNjTHgq5TOK0ZCy4ZkrUCrvSC7BsRyYKymp009zsLLB03ECM8Xe76eOLi4ubHPbk6OiI+Ph4g2clabBMiW5iV3oB5q0/ght/Ua6tk34yM7RVhUqmjTugiFqg0Qos25HZpEgB6KYt25EJjZbrJF0dy5SoBcm5JY027W8kABSU1SA5t6TjQpFRYpkSteBSefNF2pZxZLpYpkQtcLaxMOg4Ml0sU6IWDPa0h721stn7Zbi6V3+wp33HhSKjxDIlakFpVR00+q+ep9ubv3TcQB5vSixTouZotQIvbj6Ksup6uNqq4GKranS/q50FD4siHZ4BRdSM/9ufi73ZhTA3k+OLxwejv4sNz4CiZrFMifQ4ll+Gd3dlAQDeeMAXvm62AICIfrxEHunHzXyiG1TUNmB+3BHUawTu83PBzPCWL+xMBLBMiZpYsj0d6uIq9LSzwLsPB0Im46Y83RzLlOg6Ww/nY1vKOchlQOzUEHS3Mpc6EnUSLFOi/zpdWIE3vk8HADx7T38eO0q3hGVKBKC2QYP5cSmoqtMgvK89nr7LS+pI1MmwTIkAvPtzNjLOX0EPKyViHwnhIU90y1im1OXtPn4R/7c/FwCwclIQXO14nj3dOpYpdWkXymrw4uajAIDHh/XBPQNdJE5EnRXLlLosjVbg2Y0puFxVD7+etnhlrI/UkagTY5lSl/XPvSfxx+kSWJkrsHpaCFRmCqkjUSfGMqUu6U91CWJ35wAA3oz2R1+nbhInos6OZUpdTmlVHRZuSIVGKzA+uCceDnWXOhKZAJYpdSlCCCzamoZzpdXo42CFtyYE8HRRMgiWKXUp6w+ewX8yLkKpkGH1tFB0U/HCaWQYLFPqMrIuXMHyHzMBAIvG+CCgl53EiciUsEypS6iu0+CZb1NQ16DFqAFOmD3MU+pIZGJYptQlvPljBk5eqoCTjQrvTw6CnKeLkoGxTMnk7Th6HnHJZyGTAbGPBMOxm+rmDyK6RSxTMmlnS6rw2rZjAIC/jeqHYV6OEiciU8UyJZNVr9FiflwKymsbENq7O569p7/UkciEsUzJZP3j1xNIPVsKGwszfDg1BEoF3+7UfvjuIpP0e04hPv3tFADg3YcDcYe9lcSJyNSxTMnkFJbX4rmNRyEEMH1Ib9wf4CZ1JOoCWKZkUrRagRc2H0VRRS36u3TDkgcHSh2JugiWKZmUfyWexr4ThVCZyfHx9FBYKHlZPeoYLFMyGUfPluK9XdkAgCXjBqK/i43EiagrYZmSSSivqcf8uBQ0aAXG+rti+uDeUkeiLoZlSp2eEAKvf5eOMyVVcO9uiXcmBvKyetThWKbU6W0+nI8fjp6HQi7DR9NCYGellDoSdUEsU+rUTl6qwNLvMwAAz0f1xyCPHhInoq6KZfpfeXl5CA0NRXBwMIKCgvDDDz9IHYluoqZeg/lxKaiu12BoPwfMjewndSTqwmRCCCF1CGNQV1cHIQRUKhUuXryIO++8E2fOnJE6FrUg5ocMfJmkhoO1OX5aOAIuthZSR6IuzKjXTBcvXowBAwYgKioKUVFRSEhIwM6dOxEeHo7g4GDMnTsXWq0WarUaoaGhmDFjBnx8fDBv3jzdPPSN18fc3Bwq1dVLs9XU1DQ7jozDLxkX8GWSGgDw/pQgFilJzmjLNDk5GXv27EF6ejrWrVuHgwcPoqioCLGxsUhISEBqaioUCgW2bdsGAMjMzERMTAwyMjJw4MAB5OTktDhen1OnTiEgIAB+fn74+OOPO2pR6RYVlFXj5a1pAIAnh3virgHOEiciAoz228SSkpIwYcIEKJVK9OzZEyNGjIBKpUJaWhrCw8MBANXV1fDw8EBYWBh8fX3h7e0NAAgMDIRarUZWVpbe8c3p168fjh07htOnT2PWrFkYM2YMLCy4xmNMNFqBhRtSUVpVjwB3O7w8xkfqSEQAjLhM9X2UK4RAdHQ0Pv/880bT1Wq1bhMdAORyOTQaTZPxGq1Acm4Jvk89B2cbCwz2tIdCz9dX9O3bF/b29jh27BjuvPNOAy8Z3Y7Ve3KQnFsCa3MFVk8LgbmZ0W5cURdjtO/EoUOHYvv27WhoaEBBQQESExNRU1OD3bt3Iz8/HwBQXFys+78+4eHhuvG70gswZOl3mLzqRyzckIppa//A8Hf3YFd6AQDg3LlzqK6uBgBcuHABaWlp8PTkl64Zk4Oni/HR7hwAwIoJAejjaC1xIqL/Mdo10yFDhiAyMhJ+fn7w8PBAUFAQvLy8sGbNGkRHR6O+vh5KpRJr166Fvb293nk4OztjzZo1uOve+5FXVA6ZXAH7MfNhZnv1qysulNVg3voj+GRmKBQFmXj++echl8shl8uxatUqODryKy6MxeXKOizckAqtAB4O7YXxIe5SRyJqxKgPjaqsrIS1tTWKiooQERGBI0eOwMbm1i5eodEKDH93DwrKavTeLwPgameBxEV3693kJ+kJITDnq8OIP34RfR2tsWP+cFirjHY9gLooo35Hzp49G9nZ2aivr8fy5ctvuUgBIDm3pNkiBQABoKCsBsm5JYjo53Abaam9fHUgD/HHL8JcIcdH00JYpGSUjPpduXHjxtuex6XyxkWqqb6CixtebzRNYWmLS1N33vbPIsPLPH8FK346DgB4ZawP/N3tJE5EpJ9Rl6khONs0PrRJYWmLno+vvuk4kl5VXQOeiTuCugYtRvs44/FhfaSORNQso92bbyiDPe3hZmeB5j4NlQFws7t6mBQZl6XfZ+B0YSVcbFVYOTmIl9Ujo2byZaqQy7B03NXvAbrxV/Ha7aXjBnLnk5H5PvUcNh/Oh0wGxD4SAntrc6kjEbXI5MsUAMb4u+GTmaFwtWu8Ke9qZ4FPZoZijD+/vdKY5BVX4vXv0gEA8+/y4o5B6hSM+tAoQ7t2BtSl8poWz4Ai6dQ1aDH50yQczS/DnX16IG5OOMwUXeJvPnVyJr8D6noKuYxrOUZu1S/ZOJpfBjtLJWKnhrBIqdPgO5WMxm8nCvHZvtMAgPcmBcK9u6XEiYhaj2VKRuFSeQ1e2JQKAHg03AP3+blKG4joFrFMSXJarcDzG4+iqKIOPq42eP0BX6kjEd0ylilJ7rN9p5F4sggWSjk+nh4CC6VC6khEt4xlSpJKOXMZq37JBgAse8gPXs63fv0FImPAMiXJlFXXY35cChq0Ag8GumFK2B1SRyJqM5YpSUIIgde+O4b8y9W4w94Sb08M4Omi1KmxTEkSGw+dxc60ApjJZfhoaghsLZRSRyK6LSxT6nA5F8sRsyMDAPDCvQMQ0ruHxImIbh/LlDpUTb0G8+NSUFOvxQhvR/x1ZF+pIxEZBMuUOtSKnceRdaEcjt3MsWpKEOS8NgKZCJYpdZhd6QX4+o88AMA/pgTzgtxkUlim1CHOlVbj5S1pAIC/juyLkf2dJE5EZFgsU2p3DRotFsal4EpNA4Lu6I4X7h0gdSQig2OZUrv7aHcO/sy7DBuVGVZPDYG5Gd92ZHr4rqZ2lXSqCKv3ngQArJgYgN4OVhInImofLFNqNyWVdXhuYyqEAKaE9cJDQT2ljkTUblim1C6EEHhp81FcvFKLfk7WiHnIT+pIRO2KZUrt4ov9auzOugRzMzlWTwuFlXmX+oYc6oJYpmRw6efK8M7PWQCAxQ/4YmBPW4kTEbU/likZVEVtA+bHpaBOo8W9A13waLiH1JGIOgTLlAxqyffpyC2qhJudBd6bFMjL6lGXwTIlg/kuJR/bjpyDXAZ8ODUE3a3MpY5E1GFYpmQQuUWVWPxdOgBg4ej+GOxpL3Eioo7FMqXbVtegxYK4FFTWaTDE0x7P3O0ldSSiDscypdv23q4sHDtXhu5WSsRODYaCl9WjLohlSrdlb9Yl/CsxFwCwclIQ3OwsJU5EJA2WKbXZxSs1eGHzUQDAY0P7IGqgi8SJiKTDMqU20WgFntuYipLKOgx0s8UrY32kjkQkKZYptcmnv51C0qliWJkrsHp6CCyUCqkjEUmKZUq37HBeCf7x6wkAwLKH/NDPqZvEiYikxzKlW1JWVY8FcanQaAWig3ti0qBeUkciMgosU2o1IQRe2ZaGc6XV8HCwwlvj/Xm6KNF/sUyp1b5NPoOf0y9AqZBh9bQQ2FgopY5EZDRYptQq2RfK8eaOTADAy/f5ILBXd2kDERkZlindVHWdBvPjjqC2QYvI/k54Yrin1JGIjA7LlG7qzR8zceJiBZxsVFg1JQhyni5K1ATLlFq0M60AcclnIJMBH0wJhmM3ldSRiIwSy5SadbakCq9sSwMAzIvsh+HejhInIjJeLFPSq16jxcINKSivaUBI7+54Lqq/1JGIjBrLlPSKjT+BI2dKYWNhho+mhkCp4FuFqCX8DaEm9p8swj8TTgEA3pkYiDvsrSRORGT8WKbUSFFFLZ7dmAohgGmDe+OBQDepIxF1CixT0tFqBV7cfBSF5bXo79INSx4cKHUkok6DZUo6/7c/FwnZhVCZybF6WigszXlZPaLWYpkSACAtvxTv7soCALzx4EAMcLWROBFR58IyJZTX1GN+XArqNQJj/V0xY0hvqSMRdTos0y5OCIE3tqcjr7gK7t0t8c7EQF5Wj6gNWKZd3NYj57A99TwUchk+mhYMOyteVo+oLVimXdjpwgos+T4dAPDcPd4Y5GEvcSKizotl2kXVNmjwzLcpqKrTIKKvA+aN8pI6ElGnxjLtov7+UxYyC67A3tocsVODoeBl9YhuC8u0C4rPvIgvk9QAgPcnB8LF1kLaQEQmgGXaxVwoq8FLW44CAJ4Y7om7fVwkTkRkGlimXYhGK7BwQwouV9XD390WL48ZIHUkIpPBMu1C1uw9iYO5JbA2V2D1tFCozHi6KJGhsEy7iOTcEsTGnwAALB/vD09Ha4kTEZkWlmkXUFpVh2c3pEArgImh7pgY2kvqSEQmh2Vq4oQQeHlLGs6X1cDT0RrLo/2ljkRkklimJm79H3n4JfMilAoZVk8LgbXKTOpIRCaJZWrCjhdcwfKdxwEAr4z1hb+7ncSJiEwXy9REVdU14Jlvj6CuQYu7fZwxe1gfqSMRmTSWqYla9kMmThVWwsVWhZWTeFk9ovbGD9BM0I6j57Hxz7OQyYAPHgmGQzeV1JGIWqWuQYuvD6iRV1IFD3srPBrRB+ZmnWOdTyaEEFKHIMM5W1KF+z/8HeW1DZh/txdeuJdnOVHn8PefMrH291xor2skuQyYM8ITr95v/F/uyDVTE1Kv0eKZuBSU1zYgzKMHFo72ljoSUav8/adMfLYvt8l0rYBuurEXaudYf6ZWef+XbBw9WwpbCzN8OC0EZgq+vGT86hq0WPt70yK93trfc1HXoG3XHNXV1Rg9ejS6deuGV1555ZYfz982E7HvRCE+++00AOC9SUFw724pcSKi1vn6gLrRpr0+WnF1XHsyMzPDkiVLsHLlyjY9nmVqAgrLa/H8pquX1ZsZ3htj/F0lTkTUenklVc3ed3nf1zi39q+4uGExVsyfgYSEBOzcuRPh4eEIDg7G3LlzodVqoVarERoaihkzZsDHxwfz5s3TzUPfeH2USiUiIyNhadm2FRGWaSen1Qo8vykVRRW18HG1weIHjPtzJaIbedhb6Z1eez4btXlH0XP2Gjg88BwKctJRVFSE2NhYJCQkIDU1FQqFAtu2bQMAZGZmIiYmBhkZGThw4ABycnJaHG9o3AHVya39/TR+zymChVKO1dNCYKHkZfWoc3k0og9W/HS8yaZ+7bksWPaPgExhBnNbB4yKHAGVSoW0tDSEh4cDuPo5p4eHB8LCwuDr6wtv76s7XQMDA6FWq5GVlaV3fHtgmXZiKWcuY+V/sgEAMeP84O1iI3EioltnbibHnBGeevbm/69d54zwRGKaDEIIREdH4/PPP280Uq1WQ6X63/HUcrkcGo2m2fHtgZv5ndSVmnos2JCCBq3AA4FueOTOO6SORNRmr94/EH8d6Ynrv9dR5e6L6pw/8OSwO/BYSA8kJiaipqYGu3fvRn5+PgCguLhY9399wsPDb2n87eCaaSckhMDr36XjbEk1evWwxN8nBvB0Uer0Xr1/IF641+e6M6B8caJnIb5+cTL2eXggKCgIXl5eWLNmDaKjo1FfXw+lUom1a9fC3t5e7zydnZ31ju/VS/81fQMDA5Gfn4+GhgasX78eKSkpcHJyalV+ngHVCW06dBYvb02DQi7D5rkRCO3dQ+pIRO2isrIS1tbWKCoqQkREBI4cOQIbG+P8OItrpp3MyUvlWPpDBgDghXv7s0jJpM2ePRvZ2dmor6/H8uXLjbZIAa6Zdio19RqMX7MfWRfKMdzLEV/NHgy5nJv3RLeiuLgYo0ePbjTN0dER8fHxtzVflmknsuT7dHx1IA+O3czx08IRcLaxkDoSEf0X9+Z3Ev/JuICvDuQBAFZNCWaREhkZlmkncL60Gi9vSQMAPDWyLyL7t27vIhF1HJapkWvQaPHshlSUVdcjqJcdXuT1SYmMEsvUyH205ySS1SXopjLDR9NCOs1Vx4m6Gv5mGrE/Thfj4z05AIAVE/zh4WAtcSIiag7L1EhdrqzDsxtSoRXA5EG9EB3sLnUkImoBy9QICSHw0pajuHClBn2drLEs2k/qSER0EyxTI7QuSY3445dgrrh6WT0rc56oRmTsWKZGJv1cGd7+KQsA8PoDvvDraSdxIiJqDZapEamsbcCCuBTUabSIGuiCWRHtcxFbIjI8lqkRWfpDBk4XVcLNzgLvPRzIy+oRdSIsUyPxfeo5bDmcD7kMiH0kGD2szaWORES3gGVqBNRFlXj9u3QAwILR3hjS10HiRER0q1imEqtr0GLBhhRU1DZgsKc95t/tLXUkImoDlqnEVv4nC2n5ZehupcSHU4Oh4PVJiTollqmE9mZfwtrfr34j43sPB8LNzlLiRETUVixTiVy6UoMXNx0FAPwlwgP3+rlKnIiIbgfLVAJarcBzm1JRXFkHXzdbvHq/r9SRiOg2sUwl8Mlvp7D/ZDEslQp8PD0EFkqF1JGI6DaxTDvY4bzL+MevJwAAy6L90M+pm8SJiMgQWKYdqKy6HgviUqDRCjwU1BOTB/WSOhIRGQjLtIMIIfDatmM4V1qN3vZWWDHBn6eLEpkQlmkH2XDoLHYeK4CZXIbV00JgY6GUOhIRGRDLtAOcuFiOmB8yAAAvjxmAoDu6SxuIiAyOZdrOauo1eObbI6ht0GJkfyc8Obyv1JGIqB2wTNvZ8h8zceJiBRy7qbBqchDkPF2UyCSxTNvRz8cK8M3BMwCADx4JgpONSuJERNReWKbtJP9yFRZtTQMAzBvVDyO8nSRORETtiWXaDho0WizckIorNQ0I6d0dz0f1lzoSEbUzlmk7iI3PweG8y7BRmeGjqSFQKvg0E5k6/pYbWNLJIqxJOAkA+PvDAbjD3kriRETUEVimBlRcUYtnN6ZCCGDqnXfgwcCeUkciog7CMjUQIQRe3HwUl8pr4eXcDUvH+UkdiYg6EMvUQP6dmIu92YVQmcnx8fQQWJrzsnpEXQnL1ACO5Zfh3V1ZAIDFDw6Ej6utxImIqKOxTG9TRW0D5scdQb1G4D4/F8wc0lvqSEQkAZbpbVqyPR3q4ir0tLPAuw8H8rJ6RF0Uy/Q2bD2cj20p56CQy/DRtBB0tzKXOhIRSYRl2kanCyvwxvfpAIBnR3sjrI+9xImISEos0zaobdBgflwKquo0CO9rj7/d5SV1JCKSGMu0Dd79ORsZ56+gh5USsY+EQMHL6hF1eSzTW7Qn6yL+b38uAGDVlCC42llInIiIjAHL9BZcKKvBi5uvXlZv9jBP3O3jInEiIjIWLNNW0mgFnt2YgpLKOvj1tMWisQOkjkRERoRl2kr/3HsSf5wugZW5AqunhUBlxtNFieh/WKat8Ke6BLG7cwAAy6P90depm8SJiMjYsExvoqyqHgs3pEKjFZgY4o6HB/WSOhIRGSGWaQuEEFi0NQ3nSqvRx8EKb473lzoSERkplmkL1h88g10ZF6BUyLB6Wii6qcykjkRERopl2oysC1ew/MdMAMCiMT4I6GUncSIiMmYsUz2q6zR45tsU1DVocdcAJzwx3FPqSERk5Fimerz5YwZOXqqAs40K708O4mX1iOimWKY3+DHtPOKSz0ImA2KnBsOhm0rqSETUCbBMr3O2pAqvbj0GAHh6lBeG9nOUOBERdRYs0/+q12gxPy4F5bUNGOTRA8/e4y11JCLqRFim//WPX08g9WwpbC3M8OHUYJgp+NQQUeuxMQAk5hTh099OAQDefTgQvXpYSZyIiDqbLl+mheW1eG5TKoQAZgzpjbEBblJHIqJOqEuXqVYr8MLmoygsr8UAFxu88eBAqSMRUSfVpcv0X4mnse9EISyUcqyeHgILJS+rR0Rt06VONtdoBZJzS3CpvAblNQ14b1c2AGDJg37o72IjcToi6sy6TJnuSi/Ash2ZKCiraTQ9tHd3TBt8h0SpiMhUdInN/F3pBZi3/kiTIgWAI2dK8Z+MCxKkIiJTYvJlqtEKLNuRCdHM/TIAy3ZkQqNtbgQR0c2ZfJkm55boXSO9RgAoKKtBcm5Jx4UiIpNj8mV6qbz5Im3LOCIifUy+TJ1tLAw6johIH5Mv08Ge9nCzs0BzVySVAXCzs8BgT/uOjEVEJsbky1Qhl2HpuKtnNt1YqNduLx03EAo5LwBNRG1n8mUKAGP83fDJzFC42jXelHe1s8AnM0Mxxp/n4xPR7ZEJIbrMMUHXnwHlbHN1055rpERkCF1izfQahVyGiH4OiA52R0Q/B71FWlBQAFtbW3z66acSJCSizqpLlWlrvPLKK4iKipI6BhF1MkZdposXL8aAAQMQFRWFqKgoJCQkYOfOnQgPD0dwcDDmzp0LrVYLtVqN0NBQzJgxAz4+Ppg3b55uHvrGN2ffvn2wtbVFQEBARyweEZkQoy3T5ORk7NmzB+np6Vi3bh0OHjyIoqIixMbGIiEhAampqVAoFNi2bRsAIDMzEzExMcjIyMCBAweQk5PT4vgbNTQ0YMmSJVi2bFlHLiYRmQijvWpUUlISJkyYAKVSiZ49e2LEiBFQqVRIS0tDeHg4AKC6uhoeHh4ICwuDr68vvL2vfgleYGAg1Go1srKy9I7XZ82aNZg8eTLs7Xm8KRHdOqMtU30HGQghEB0djc8//7zRdLVaDZXqf99vL5fLodFomh2vT3JyMhITE7Fy5UqUlpZCoVDAysoKs2bNuv2FISKTZ7Sb+UOHDsX27dvR0NCAgoICJCYmoqamBrt370Z+fj4AoLi4WPd/fcLDw1s9/ptvvkFeXh7UajWeffZZrFixgkVKRK1mtGumQ4YMQWRkJPz8/ODh4YGgoCB4eXlhzZo1iI6ORn19PZRKJdauXdvsprmzs7Pe8b169ergpSEiU2fUB+1XVlbC2toaRUVFiIiIwJEjR2Bjw68XISLjY7RrpgAwe/ZsZGdno76+HsuXL2eREpHRMuo10/ZQXFyM0aNHN5rm6OiI+Ph4iRIRkSnocmVKRNQejHZvPhFRZ8IyJSIyAJYpEZEBsEyJiAyAZUpEZAAsUyIiA2CZEhEZAMuUiMgAWKZERAbAMiUiMgCWKRGRAbBMiYgMgGVKRGQALFMiIgNgmRIRGQDLlIjIAFimREQGwDIlIjIAlikRkQGwTImIDIBlSkRkACxTIiIDYJkSERkAy5SIyADMpA5A1FlotALJuSW4VF4DZxsLDPa0h0IukzoWGQmWKVEr7EovwLIdmSgoq9FNc7OzwNJxAzHG303CZGQsZEIIIXUIImO2K70A89YfwY2/KNfWST+ZGcpCJX5mStQSjVZg2Y7MJkUKQDdt2Y5MaLRcJ+nqWKZELUjOLWm0aX8jAaCgrAbJuSUdF4qMEsuUqAWXypsv0raMI9PFMiVqgbONhUHHkelimRK1YLCnPdzsLNDcAVAyXN2rP9jTviNjkRFimRK1QCGXYem4gQDQpFCv3V46biCPNyWWKdHNjPF3wyczQ+Fq13hT3tXOgodFkQ6PMyVqJZ4BRS1hmRIRGQA384mIDIBlSkRkACxTIiIDYJkSERkAy5RIQmq1GtbW1ggODkZwcDCWLFkidSRqI17PlEhiAQEB+OOPP6SOQbeJa6ZEN1i8eDEGDBiAqKgoREVFISEhATt37kR4eDiCg4Mxd+5caLVaqNVqhIaGYsaMGfDx8cG8efN089A3nkwby5ToOsnJydizZw/S09Oxbt06HDx4EEVFRYiNjUVCQgJSU1OhUCiwbds2AEBmZiZiYmKQkZGBAwcOICcnp8Xx+mRmZiIkJARRUVFIS0vrqEUlA+NmPtF1kpKSMGHCBCiVSvTs2RMjRoyASqVCWloawsPDAQDV1dXw8PBAWFgYfH194e3tDQAIDAyEWq1GVlaW3vH6uLm5ITc3Fw4ODti3bx8mT56M7OzsjllYMiiWKdF19J0QKIRAdHQ0Pv/880bT1Wo1VCqV7rZcLodGo2l2vD4qlUo3j5EjR8Lc3BxFRUVwdHS8zSWhjsbNfKLrDB06FNu3b0dDQwMKCgqQmJiImpoa7N69G/n5+QCA4uJi3f/1CQ8Pb/X4wsJCaDQaAEBGRgYqKirg4OBg4KWijsA1U6LrDBkyBJGRkfDz84OHhweCgoLg5eWFNWvWIDo6GvX19VAqlVi7di3s7fVfw9TZ2Vnv+F69ejUZu2/fPixZsgRKpRIqlQpff/01ZDJePKUz4oVOiG5QWVkJa2trFBUVISIiAkeOHIGNjY3UscjIcc2U6AazZ89GdnY26uvrsXz5chYptQrXTIk6QHFxMUaPHt1omqOjI+Lj4yVKRIbGMiUiMgDuzSciMgCWKRGRAbBMiYgMgGVKRGQALFMiIgNgmRIRGQDLlIjIAFimREQGwDIlIjIAlikRkQGwTImIDIBlSkRkACxTIiIDYJkSERkAy5SIyABYpkREBsAyJSIyAJYpEZEBsEyJiAzg/wFMPMcR5qe+QwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(4, 4))\n", "\n", "grnplot(models[0], ax=ax)\n", "\n", "ax.set(title=\"cluster_1 (prototype1)\");" ] }, { "cell_type": "markdown", "id": "a55517a8-24a0-44e9-a011-8b89bfb89ab4", "metadata": {}, "source": [ "## Customize plot aesthetics" ] }, { "cell_type": "markdown", "id": "2bd2c75c-bd37-4ed2-b327-c69b9bb4ef19", "metadata": {}, "source": [ "### color" ] }, { "cell_type": "code", "execution_count": 5, "id": "a75d351f-af88-4dd9-a580-c54bc21eea07", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVMAAAFeCAYAAAA41fyqAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAALTlJREFUeJzt3XtcVHXixvHPgDgiioqI18ILpHiDlBLxgmWkdhE1c9ssNzVdzbba3VqrNcN1bWvNjdbsorWtZb+2vORaVlteyBBTy8wLSmYOhZIKeQkEBeb8/jhJIgMiDpxheN6vly/PnPnO4TkMPJwzc84Zm2EYBiIickl8rA4gIuINVKYiIm6gMhURcQOVqYiIG6hMRUTcQGUqIuIGKlMRETdQmYqIuIHKVETEDVSmbpKcnIzNZiM5OdnqKB7tnnvuIT4+3uoYddZtt93GmDFjrI7hlVSmtURqaiqJiYkcP37ckq+fm5vL448/ztChQwkKCsJms/Hvf//7opZx4MABXn75ZR599NHqCXkBp06dIjEx8ZL+4Fn9PJRnzpw5DB8+nJYtW2Kz2UhMTHQ5bvr06SxfvpyvvvqqZgPWASrTWiI1NZVZs2ZZ9kucnZ3NX/7yF/bs2UNkZGSVlvHss8/SoUMHrrnmGjenq5xTp04xa9asSy5TK5+H8syYMYOtW7dy5ZVXVjjuyiuvJDo6mnnz5tVQsrpDZVrHnTp1qlLjWrduTVZWFhkZGcydO/eiv05hYSFvvPGG23cxK5vf2x04cICsrCyWLFlywbFjxoxhxYoV5Obm1kCyukNlWkkHDx5k4sSJtGnTBrvdTocOHZg6dSpnzpwp9zHt27fnrrvuKjN/0KBBDBo0qNS8+fPn061bNxo2bEizZs2Ijo7m//7v/wBITEzkoYceAqBDhw7YbDZsNhsOh6Pk8UuWLKF37974+/sTFBTEbbfdxvfff1/m63bv3p0vvviCgQMH0rBhw0rvctvtdlq1alWpsa6kpKSQnZ3NddddV2r+2dea33rrLR599FFatWpFQEAAw4cPv6j8R44cYeLEibRs2ZIGDRoQGRnJ4sWLSx7rcDho0aIFALNmzSr5Hp67O7xu3ToGDBhAQEAATZs2JSEhgT179pTcX9HzEBcXV+4We+fOnRkyZEhJDpvNxtNPP80zzzxDaGgo/v7+xMXFsWvXrjKP3bt3L6NHjyYoKIgGDRoQHR3NqlWryoxr3759ed/6MuLj48nLy+Pjjz+u9GPkwupZHaA2OHToEFdffTXHjx9n8uTJdOnShYMHD7Js2TJOnTpF/fr1L2n5ixYt4r777mP06NHcf//9FBQUsGPHDjZv3sztt9/OqFGj+Prrr3nzzTd55plnCA4OBigphzlz5vDYY48xZswY7r77bo4ePcr8+fMZOHAgX375JU2bNi35Wjk5OQwbNozbbruNO+64g5YtW15S9spKTU3FZrOVuxs6Z84cbDYb06dP58iRIyQlJXHdddexfft2/P39K8yfn5/PoEGD+Oabb7j33nvp0KEDS5cu5a677uL48ePcf//9tGjRghdeeIGpU6cycuRIRo0aBUDPnj0BWLNmDcOGDaNjx44kJiaSn5/P/Pnz6devH9u2baN9+/YVPg933nknkyZNYteuXXTv3r0k79atW/n666+ZMWNGqfV97bXX+Omnn5g2bRoFBQU8++yzXHvttezcubPkOdm9ezf9+vWjbdu2PPzwwwQEBPD2228zYsQIli9fzsiRI6v0XHTt2hV/f382btxY5WWIC4Zc0Lhx4wwfHx9j69atZe5zOp2GYRjG+vXrDcBYv359yX2hoaHGb37zmzKPiYuLM+Li4kpuJyQkGN26dasww9y5cw3AOHDgQKn5DofD8PX1NebMmVNq/s6dO4169eqVmh8XF2cAxosvvljh17qQrVu3GoDx6quvVvoxd9xxh9G8efMy889+39q2bWucPHmyZP7bb79tAMazzz57wfxJSUkGYCxZsqRk3pkzZ4y+ffsajRo1Klnu0aNHDcB4/PHHy+SIiooyQkJCjJycnJJ5X331leHj42OMGzeuZF55z8Px48eNBg0aGNOnTy81/7777jMCAgKM3NxcwzAM48CBAwZg+Pv7G5mZmSXjNm/ebADG73//+5J5gwcPNnr06GEUFBSUzHM6nUZsbKwRHh5eZh0utI7nuuKKK4xhw4ZVOEYujnbzL8DpdLJy5UpuvvlmoqOjy9xvs9ku+Ws0bdqUzMxMtm7detGPXbFiBU6nkzFjxpCdnV3yr1WrVoSHh7N+/fpS4+12O+PHj7/kzBcrJyeHZs2alXv/uHHjaNy4ccnt0aNH07p1a95///1S41zlf//992nVqhW//vWvS+b5+flx3333kZubyyeffFJhtqysLLZv385dd91FUFBQyfyePXsSHx9fJoMrTZo0ISEhgTfffBPj5+utFxcX89ZbbzFixAgCAgJKjR8xYgRt27YtuX311VfTp0+fkq/1448/sm7dOsaMGcNPP/1U8rzm5OQwZMgQ9u3bx8GDBy+YqzzNmjUjOzu7yo+XslSmF3D06FFOnjxZatfN3aZPn06jRo24+uqrCQ8PZ9q0aWzcuLFSj923bx+GYRAeHk6LFi1K/duzZw9HjhwpNb5t27aX/LJEVRkVfKhDeHh4qds2m42wsLBSrwuD6/wZGRmEh4fj41P6xzkiIqLk/oqcvb9z585l7ouIiCA7O5u8vLwKlwHmH4TvvvuOTz/9FDBfOjh8+DB33nlnmbHnry/AFVdcUbK+33zzDYZh8Nhjj5V5Xh9//HGAMs/txTAMwy0bAvILvWZajcr7YS0uLsbX17fkdkREBOnp6bz33nt8+OGHLF++nOeff56ZM2cya9asCr+G0+nEZrPxwQcflFrmWY0aNSp1+9zXH2tS8+bNOXbs2CUvx6r8lTFkyBBatmzJkiVLGDhwIEuWLKFVq1Zl3nSrDKfTCcCDDz5Y8ubV+cLCwqqc9dixYy4LXapOZXoBLVq0IDAw0OU7rRfSrFkzl8cjZmRk0LFjx1LzAgIC+NWvfsWvfvUrzpw5w6hRo5gzZw6PPPIIDRo0KLeYO3XqhGEYdOjQgSuuuOKiM9aULl268MYbb3DixAmaNGlS5v59+/aVum0YBt98803JG0QVCQ0NZceOHTidzlJbp3v37i25H8r/43b2/vT09DL37d27l+Dg4JLd9Iq25nx9fbn99tv597//zVNPPcXKlSuZNGmSyz9y568vwNdff13yrvzZnw8/P78qlXFFioqK+P777xk+fLhbl1vXaTf/Anx8fBgxYgTvvvsun3/+eZn7K9p17dSpE5999lmpw6fee++9Mof85OTklLpdv359unbtimEYFBYWApT8Mp9fzqNGjcLX15dZs2aVyWIYRpllW6Vv374YhsEXX3zh8v6z726ftWzZMrKyshg2bNgFl33DDTfwww8/8NZbb5XMKyoqYv78+TRq1Ii4uDgAGjZsCJT9HrZu3ZqoqCgWL15c6r5du3bx0UcfccMNN5TMK+95OOvOO+/k2LFj/Pa3vyU3N5c77rjD5biVK1eWes1zy5YtbN68uWR9Q0JCGDRoEC+99BJZWVllHn/06NFyvhsXlpaWRkFBAbGxsVVehpSlLdNKeOKJJ/joo4+Ii4tj8uTJREREkJWVxdKlS0lJSSl16NG57r77bpYtW8bQoUMZM2YM+/fvZ8mSJXTq1KnUuOuvv55WrVrRr18/WrZsyZ49e3juuee48cYbS96U6d27NwB//vOfue222/Dz8+Pmm2+mU6dO/PWvf+WRRx7B4XAwYsQIGjduzIEDB3jnnXeYPHkyDz74oFu+D8899xzHjx/n0KFDALz77rtkZmYC8Lvf/c7lFudZ/fv3p3nz5qxZs4Zrr722zP1BQUH079+f8ePHc/jwYZKSkggLC2PSpEkXzDV58mReeukl7rrrLr744gvat2/PsmXL2LhxI0lJSSXfQ39/f7p27cpbb73FFVdcQVBQEN27d6d79+7MnTuXYcOG0bdvXyZOnFhyaFSTJk1KHYta3vNwtmSvvPJKunfvztKlS4mIiKBXr14uM4eFhdG/f3+mTp3K6dOnSUpKonnz5vzpT38qGbNgwQL69+9Pjx49mDRpEh07duTw4cNs2rSJzMzMUqeEvv7662RkZJScxLBhwwb++te/AmbBn936Bvj4449p2LChrpHgbtYcRFD7ZGRkGOPGjTNatGhh2O12o2PHjsa0adOM06dPG4bh+tAowzCMefPmGW3btjXsdrvRr18/4/PPPy9zaNRLL71kDBw40GjevLlht9uNTp06GQ899JBx4sSJUsuaPXu20bZtW8PHx6fM4TnLly83+vfvbwQEBBgBAQFGly5djGnTphnp6eklY+Li4i54CFZFQkNDDcDlv/MPFXLlvvvuM8LCwkrNO/t9e/PNN41HHnnECAkJMfz9/Y0bb7zRyMjIKDW2ovyHDx82xo8fbwQHBxv169c3evTo4fLQrdTUVKN3795G/fr1yxxCtGbNGqNfv36Gv7+/ERgYaNx8881GWlpamWVU9DwYhmH8/e9/NwDjiSeeKPPYs4dGzZ0715g3b55x2WWXGXa73RgwYIDx1VdflRm/f/9+Y9y4cUarVq0MPz8/o23btsZNN91kLFu2rMz3przn5vyfyT59+hh33HGHy++jVJ3KVGrM/v37DT8/P2PNmjUl886W6dKlSy1M5l5JSUmGzWYr88fAMEqXqRW+/PJLw2azGV9++aUlX9+b6TVTqTEdO3Zk4sSJPPnkk1ZHqTaGYfDKK68QFxfH5ZdfbnWcMp588klGjx5NVFSU1VG8jl4zreNyc3MveMGLFi1auHxHuipeeOEFtyzH0+Tl5bFq1SrWr1/Pzp07+e9//2t1JJf+85//WB3Ba6lM67inn376gseyHjhw4KIupFEXHT16lNtvv52mTZvy6KOP6rCjOshmGBUc2yNe79tvv+Xbb7+tcEz//v1p0KBBDSUSqZ1UpiIibqA3oERE3EBlKiLiBipTERE3UJmKiLiBylRExA1UpiIibqAyFRFxA5WpiIgbqExFRNxAZSoi4gYqUxERN1CZioi4gcpURMQNVKYiIm6gMhURcQOVqYiIG6hMRUTcQGUqIuIGKlMRETdQmYqIuIHKVMRCGzduJDo6mh49etC3b1927txpdSSpIn06qYiFdu/eTZMmTWjXrh1r1qxhzpw5rF+/3upYUgXaMhU5z4wZM+jcuTPx8fHEx8eTnJzM6tWriYmJISoqiilTpuB0OnE4HPTq1YuxY8fSpUsXpk6dWrIMV+Nd6datG+3atQOgd+/eHDx4sEbWUdxPZSpyji1btrBu3Tp27drF4sWL2bx5M9nZ2SQlJZGcnMz27dvx9fVlxYoVAKSlpZGYmMju3bvZtGkT+/btq3B8RV577TUGDx5c3aso1aSe1QFEPElqaiojR47Ez8+PNm3aMGDAAOx2Ozt27CAmJgaA/Px8QkNDiY6OJiIigvDwcAB69uyJw+Fg7969LsdX5IsvvuDFF18kJSWleldQqo3KVOQcrt5CMAyDhIQEFi5cWGq+w+HAbreX3Pbx8aG4uLjc8eXJzMzk9ttv5+2336Z58+aXtgJiGe3mi5wjNjaWlStXUlRURFZWFikpKRQUFLB27VoyMzMByMnJKZl2JSYmptLj8/LySEhIYN68eURGRrp/haTGqExFztGnTx/i4uLo1q0bv/nNb4iMjCQsLIwFCxaQkJBAz549uf766zly5Ei5ywgJCan0+FdffZW9e/cyY8YMoqKiiI2Nra5Vk2qmQ6NEzpOXl0dAQADZ2dn07duXbdu20bhxY6tjiYfTa6Yi55kwYQLp6ekUFhYye/bsX4rUWQwZqZB7GBq1hNBY8PG1Nqx4DG2ZilRG2ir4cDqcPPTLvMA2MPQp6Dr8gg/Pyckpc9hTcHAwa9ascXdSsYjKVORC0lbB2+OA839VbOZ/Y16rVKGKd9MbUCIVcRabW6RlipRf5n34sDlO6jSVqUhFMlJL79qXYcDJg+Y4qdNUpiIVyT3s3nHitVSmIhVp1NK948RrqUxFKhIaCw0rOsXTBoFtzXFSp6lMRSqSf6yCN5d+fjd/6JM63lRUpiLlcjph5VQoOA6N20DjVqXvD2yjw6KkhM6AEinP5hdg30fga4exSyEkQmdASbl00L6IK4e+hJfjwVkINzwNV0+yOpF4OO3mi5zv9E+wbIJZpF1ugqvutjqR1AIqU5HzrX4QfvwWAtvB8Plgs1mdSGoBlanIuba/CTv+AzYfuGURNAyyOpHUEipTkbOyv4HVfzSnBz2iY0floqhMRQCKTsOy8VCYB+0HwIA/Wp1IahmVqQjAmkT4YQf4B8GohTrkSS6aylQk/UP47HlzesTz5sH4IhdJZSp128lD5llOAH2mQudh1uaRWktlKnWXsxhWTIb8H6FVT4ifZXUiqcVUplJ3ffoPcHwKfgEw+lWoZ7c6kdRiKlOpm777DJL/Zk7f+DQEh1mbR2o9lanUPad+hOV3g1EMPcZA5K+tTiReQGUqdYthwKrfwYnvIagj3PQPnS4qbqEylbrl81dg73vg4wej/wX2xlYnEi+hMpW64/Bu+PBRc/q6RGhzpaVxxLuoTKVuOHMKlo6H4tMQFg8x91idSLyMylTqhg8fhux08wr5I14AH/3oi3vpJ0q8367lsG0xYDPPu2/UwupE4oVUpuLdjjng3QfM6QF/gI6DLAwj3kxlKt6ruBCWTYTTJ6Hd1eY1SkWqicpUvNf6OXDwc7A3gVteBl8/qxOJF1OZinfavw5Skszp4f+EZqGWxhHvpzIV75N7BFb8FjCg93joNsLqRFIHqEzFuzid8M4UyDsCLSJg6N+sTiR1hMpUvMum52D/WqjXAG59Ffz8rU4kdYTKVLzHwS9g7c8XeB76NwiJsDaP1CkqU/EOBSdh2QRwFkHEcPO1UpEapDKV2s8w4L3fmwfoN7ncfPdel9WTGqYyldpv+xuwaxnYfGH0K+DfzOpEUgepTKV2O/o1vP+QOX3No3DZ1dbmkTpLZfqzjIwMevXqRVRUFJGRkaxatcrqSHIhhQXm66SFp6DDQOj/e6sTSR1mMwzDsDqEJzhz5gyGYWC32zl8+DBXXXUV3333ndWxpCLv/wm2vAQNg2FKCgS2tjqR1GEevWU6Y8YMOnfuTHx8PPHx8SQnJ7N69WpiYmKIiopiypQpOJ1OHA4HvXr1YuzYsXTp0oWpU6eWLMPVeFfq16+P3W5+1G9BQUG548RD7F1tFinAyBdVpGI5jy3TLVu2sG7dOnbt2sXixYvZvHkz2dnZJCUlkZyczPbt2/H19WXFihUApKWlkZiYyO7du9m0aRP79u2rcLwr+/fvp0ePHnTr1o3nnnuuplZVLtaJg/DfaeZ033shPN7aPCJAPasDlCc1NZWRI0fi5+dHmzZtGDBgAHa7nR07dhATEwNAfn4+oaGhREdHExERQXh4OAA9e/bE4XCwd+9el+PL06lTJ3bu3Mm3337LuHHjGDp0KA0aNKj+lZXKcxbDikmQfwxaR8Hgx61OJAJ4cJm6einXMAwSEhJYuHBhqfkOh6NkFx3Ax8eH4uLisuOdxZCRCjuXmR9fERoLPr5lvk7Hjh0JCgpi586dXHXVVe5dMbk0G+ZCxkao38j8dNF69a1OJAJ48G5+bGwsK1eupKioiKysLFJSUigoKGDt2rVkZmYCkJOTUzLtSkxMzC/j01aRMyeCzPk3wPKJsPgmSOoOaea79gcPHiQ/Px+AH374gR07dtChQ4fqX1GpPMdG+OQpc/qmZ6B5J2vziJzDY7dM+/TpQ1xcHN26dSM0NJTIyEjCwsJYsGABCQkJFBYW4ufnx6JFiwgKCnK5jJCQEHP8kEEUZh/AzxcW3exPu8CfB5zMgrfHwZjXSDvozx/+8Ad8fHzw8fFh3rx5BAcH19wKS8VO/QjL7wbDCZG3Q88xVicSKcWjD43Ky8sjICCA7Oxs+vbty7Zt22jcuPHFLcRZbG6BnjxUzgAbBLaBB3a63OUXD2AY8J/bIf19aB4Gkz8BeyOrU4mU4rFbpgATJkwgPT2dwsJCZs+effFFCuZrpOUWKYABJw+a4zoMqHJWqUZbFplF6lvffJ1URSoeyKPL9K233rr0heQeLnUz55STwa+dKjUvuKGNNbeUHice4oed8NEMczr+L9A60to8IuXw6DJ1i0YtS91s3tCH7VNcbNmcN048wJk8WDoeik/DFUOhzxSrE4mUy2PfzXeb0FjzNVHKuySbDQLbmuPEs7z/J8jZB41bQ8LzuqyeeDTvL1MfXxj68+E0ZQr159tDn9SbT55m5zLYvgSwwahFENDc6kQiFfL+MgXoOhzGvFb2/O3ANub8rsOtySWu/fgtvPuAOT3wIb0xKLWCRx8a5XZnz4DKPVzhGVBioaIz8K8hcGgbXN4XfvMe+Hr/S/tS+9Wtn1IfX23leLp1s80ibdDU3L1XkUotUTd286V2+GYNpP7TnE5YAE0vszaPyEVQmYpn+OkwvPPzoU9X3Q0RN1mbR+QiqUzFek4nvDMZ8o5CSDe4/q9WJxK5aCpTsV7qs/BtMtTzh1tfBT9/qxOJXDSVqVgr83NY9/OW6A1/hxadrc0jUkUqU7FO/nFYNh6cRdBtFFx5p9WJRKpMZSrWMAx47wE4/h00DYWbk3S6qNRqKlOxxrbXYPc74FPPvKxegyZWJxK5JCpTqXlH9sIH083pa2dAu2hr84i4gcpUalZhPiybAEX50PEaiL3f6kQibqEylZr10Qw4shsCWsDIl8BHP4LiHfSTLDUnbRVsfdmcHvkSNNYFucV7qEylZhz/Hlbda07H3gdhg63NI+JmKlOpfsVF5sc0F5yAtr3h2sesTiTidipTqX6fPAXffwb2QLjlFahX3+pEIm6nMpXqdWADbJhrTt/0DAR1sDaPSDVRmUr1ycuBFZMBA668A3qMtjqRSLVRmUr1MAz47z3wUxYEXwHD/m51IpFqpTKV6rH5Rfj6Q/C1m6eL1g+wOpFItVKZivtlfQUfzzSnh8yBVj2szSNSA1Sm4l6nf4Kl46H4DHS5yfwIEpE6QGUq7vX+Q/DjfghsC8Pn67J6UmeoTMV9vnoLvnoTbD5wy8vQMMjqRCI1RmUq7pGzH1b/wZyOexhCY63NI1LDVKZy6YrOmJfVO5MLof1h4INWJxKpcSpTuXRrZ0HWdvBvBqMWgo+v1YlEapzKVC7N1x/BpufM6YTnoUlba/OIWERlKlV3MgtWTjGnr/4tdLnB2jwiFlKZStU4i+GdyXAqxzwoP/4vVicSsZTKVKom5RnzilB+ATD6VfBrYHUiEUupTOXifbcZ1j9hTt8wF4LDrc0j4gFUpnJx8o/B8olgFEOPWyHqdqsTiXgElalUnmHAqvvgxPfQrAPc+A+dLiryM5WpVN4Xr8KeVeDjZ15Wr0Gg1YlEPIbKVCrncBp8+Ig5fd3j0LaXtXlEPIzKVC7szCnzdNGiAgi7DmKmWZ1IxOOoTOXC/vcIHN0DjVrCiBfBRz82IufTb4VUbPc78MW/ARuMfAkatbA6kYhHUplK+Y5lwKr7zen+v4dO11ibR8SDqUzFteJCWH43nD4B7a6Cax61OpGIR1OZimvJf4PMLWBvAre8Ar5+VicS8WgqUynr22T49B/m9PBnoVmopXFEagOVqZSWexRWTAYM6H0XdBtpdSKRWkFlKr9wOmHlVMg9DC0iYMjfrE4kUmuoTOUXnz0P33wM9RqYp4vWb2h1IpFaQ2UqpoPbYE2iOT3kCWjZ1dI4IrWNylSg4KR5uqizECKGQ/QEqxOJ1Doq07rOMGD1H+HYAWhyGQz/py6rJ1IFKtO67qs3YefbYPM1jyf1b2Z1IpFaSWVal2V/A6sfNKeveQQu72NtHpFaTGVaVxWdhmV3QWEetB8A/f9gdSKRWk1lWld9PBN+2AkNm8OoReDja3UikVpNZVoXpX8Am180p0e8AIGtrc0j4gVUpnXNyUOw8h5zOmYaXDHE2jwiXkJlWpc4i2H5JMj/EVpHmp/lJCJuoTKtSz6dBxkpUL8RjH4V6tmtTiTiNVSmdUVGqnmNUoAb50HzTtbmEfEyKtO64NSP5u694YTIX0PkbVYnEvE6KlNvZxiw6ndwMhOCOsENT1udSMQrqUy93daXYe974ONnXlbP3sjqRCJeSWXqzX7YBf/7szkd/xdoE2VpHBFvpjL1VmfyYNl4KD4N4UMgZqrViUS8msrUW30wHbK/hsatYcTzuqyeSDWrZ3UAqQa7lsOXrwM2GLUQAoKtTiRSOUVnYOsiOOaAZu3hqklQr77VqSrFZhiGYXUIcaNjDnhxAJw+CQMfgmtnWJ1IpHI+egw2PWcewneWzQf63gvXz7YuVyVpy9SbFBeaHz9y+iRcFgNxD1udSKRyPnoMUv9Zdr7h/GW+hxeqXjP1Jutmw8EvoEETuOVl8NXfSqkFis6YW6QV2bTAHFeN8vPzGTx4MI0aNeLhhy9+Q0Rl6i2+WQsbnzWnExZA08uszSNSWVsXld61d8UoNsdVo3r16jFz5kzmzp1bpcerTL1B7hF4Z4o5HT0RIm62No/IxTjmKPeuGesK6PxcLvGv5xF/z1ySk5NZvXo1MTExREVFMWXKFJxOJw6Hg169ejF27Fi6dOnC1Km/HAroarwrfn5+xMXF4e/vX6XVUJnWdk4nvPNbyDsCId1gyByrE4lcnGbtXc7ecrCYdQeK2TU1gMUj/Nn89WGys7NJSkoiOTmZ7du34+vry4oVKwBIS0sjMTGR3bt3s2nTJvbt21fheHfTi2q13ab5sH8d1PM3Txf1q9pfVRHLXDUJPppRZlc/9fsiRnaph5+vjTaB9RgQ1x+73c6OHTuIiYkBzNc5Q0NDiY6OJiIigvDwcAB69uyJw+Fg7969LsdXB5VpbZb5Oaz9izk97CkI6WJtHpGqqFffPPzpvHfzSx202XcabNmGYRgkJCSwcOHCUmMdDgd2+y/X5/Xx8aG4uLjc8dVBu/m1VcEJ8zAoZxF0Gwm9xlmdSKTqrp8NsfeZx5X+LPYyX1amF1HU516yetxDSkoKBQUFrF27lszMTABycnJKpl2JiYm5qPGXQlumtZFhwHu/h+MZ0PRyuPlZnS4qtd/1s+Hax0rOgOozpD1xwd/R7YFlhIZuJTIykrCwMBYsWEBCQgKFhYX4+fmxaNEigoKCXC4yJCTE5fh27dq5HN+zZ08yMzMpKipiyZIlfPnll7Ro0aJS8XUGVG207XVYdS/YfGHC/+Cyq6xOJFIt8vLyCAgIIDs7m759+7Jt2zYaN25sdSyXtGVa2xxNhw/+ZE5fO0NFKl5twoQJpKenU1hYyOzZsz22SEFbprVLYQG8PBgO74KOg+COd8BHL3uLXIycnBwGDx5cal5wcDBr1qy5pOWqTGuT1Q+arycFtIApG6FxS6sTicjPtFlTW+x575fT6Ua+qCIV8TAq09rgRCb8d5o5Hfs7CLvO2jwiUobK1NMVF5kf01xwHNr0gmtnWp1IRFxQmXq6DX+H71KhfmMY/Uqtueq4SF2jMvVkjhTY8PPlwG5OgqCOlsYRkfKpTD3VqR/N3XvDCVF3QI/RVicSkQqoTD2RYcDKe+CnQ9A8HG74u9WJROQCVKaeaMtC+PoD8K1vXlavfoDViUTkAlSmnibrK/PajgDXz4HWPa3NIyKVojL1JKdzzcvqFZ+BzjfC1ZOsTiQilaQy9SQf/AlyvoHAtpDwnC6rJ1KLqEw9xY6lsP0N8+K4oxZBQ9fXZxQRz6Qy9QQ5+82LPQPETYf2/azNIyIXTWVqtaIzsHwinPkJQvvBwIesTiQiVaAytdraWXDoS/BvZu7e+/hanUhEqkBlaqV9H8Om58zphAXQpK21eUSkylSmVvnpB3hnijl99WTocqO1eUTkkqhMreB0worJcCobWvaA+NlWJxKRS6QytcLGZ+DAJ+DXEG59FfwaWJ1IRC6RyrSmfb8F1s0xp2+YC8Hh1uYREbdQmdak/OOwbCIYxdB9NESNtTqRiLiJyrSmGAa8ez+c+A6atYebntHpoiJeRGVaU7YthrSV4FPPvKxeg0CrE4mIG6lMa8KRPfDBdHN68OPQtre1eUTE7VSm1a0wH5aOh6IC6DQY+t5rdSIRqQYq0+r2v0fh6B4ICIGRL4KPvuUi3ki/2dUp7b/w+b/M6VEvQaMQa/OISLVRmVaX49/Bqt+Z0/1/D52utTaPiFQrlWl1KC6C5XdDwQlodxVc82erE4lINVOZVofkv8H3m8EeCLe8DL5+VicSkWqmMnW3bz+BT+eZ0zc/ax6gLyJeT2XqTnnZ5tWgMKDXOOg+yupEIlJDVKbuYhiwcirk/gDBnWHoU1YnEpEapDJ1l8+eh30fQb0G5mX16je0OpGI1CCVqTsc+hI+ftycHjIHWnazNo+I1DiV6aU6/RMsmwDOQuhyE0RPtDqRiFhAZXqpVj8IP34Lge1g+HxdVk+kjlKZXortb8KO/4DNF0a/Ag2DrE4kIhZRmVZV9jew+o/m9KBH4PIYa/OIiKVUplVRdBqWjYfCPGg/AAb8wepEImIxlWlVrEmEH3aAfxCMWgg+vlYnEhGLqUwv1tf/M48pBfP6pIFtrM0jIh5BZXoxTh4yz3ICiLkHrhhibR4R8Rgq08pyFpvn3Z/KgVY94bpEqxOJiAdRmVbWp/8Ax6fgFwCjX4V6dqsTiYgHUZlWxnefmdcoBbhxHgSHWZtHRDyOyvRC8o+ZV803iqHnbRD1a6sTiYgHUplWxDDMz3E68T0EdYQbn7Y6kYh4KJVpRT5/Bfa8Cz5+MPpfYG9sdSIR8VAq0/Ic3g0fPmpOx8+CNldam0dEPJrK1JUzp2DpeCg+DeHXm8eUiohUQGXqyocPQ3Y6NGoFI17QZfVE5IJUpufbtQK2LQZscMsiCAi2OpGI1AIq03Mdc8C795vTA/4IHQZaGkdEag+V6VnFhbBsIpw+CZf1Ma9RKiJSSSrTs9bPgYOfQ4MmcMvL4FvP6kQiUouoTAH2r4eUJHN6+HxoermlcUSk9lGZ5h6Bd34LGBA9AbomWJ1IRGqhul2mTie8MwVyD0NIVxjyhNWJRKSWqttluuk52L8W6vmbp4v6+VudSERqqbr1LouzGDJSzS3RghOwdpY5f+jfICTC2mwiUqvVnTJNWwUfTjc/euRc7a6G3ndZEklEvEfd2M1PWwVvjytbpACZW8wrQ4mIXALvL1NnsblFilHOAJt5Lr6zuCZTiYiX8f4yzUh1vUVawoCTB81xIiJV5P1lmnvYveNERFzw/jJt1NK940REXPD+Mg2NhcA2QHnXJLVBYFtznIhIFXl/mfr4wtCnfr5xfqH+fHvok+Y4EZEq8v4yBeg6HMa8BoGtS88PbGPO7zrcmlwi4jVshmGUd8yQ9zn3DKhGLc1de22Riogb1I0t07N8fKHDAOgx2vzfRZFmZWURGBjIiy++aEFAEamt6laZVsLDDz9MfHy81TFEpJbx6DKdMWMGnTt3Jj4+nvj4eJKTk1m9ejUxMTFERUUxZcoUnE4nDoeDXr16MXbsWLp06cLUqVNLluFqfHk2bNhAYGAgPXr0qInVExEv4rFlumXLFtatW8euXbtYvHgxmzdvJjs7m6SkJJKTk9m+fTu+vr6sWLECgLS0NBITE9m9ezebNm1i3759FY4/X1FRETNnzmTWrFk1uZoi4iU89qpRqampjBw5Ej8/P9q0acOAAQOw2+3s2LGDmJgYAPLz8wkNDSU6OpqIiAjCw8MB6NmzJw6Hg71797oc78qCBQu49dZbCQoKqpkVFBGv4rFl6uogA8MwSEhIYOHChaXmOxwO7HZ7yW0fHx+Ki4vLHe/Kli1bSElJYe7cuRw/fhxfX18aNmzIuHHjLn1lRMTreexufmxsLCtXrqSoqIisrCxSUlIoKChg7dq1ZGZmApCTk1My7UpMTEylx7/xxhtkZGTgcDh44IEHmDNnjopURCrNY7dM+/TpQ1xcHN26dSM0NJTIyEjCwsJYsGABCQkJFBYW4ufnx6JFi8rdNQ8JCXE5vl27djW8NiLi7Tz6oP28vDwCAgLIzs6mb9++bNu2jcaNG1sdS0SkDI/dMgWYMGEC6enpFBYWMnv2bBWpiHgsj94yrQ45OTkMHjy41Lzg4GDWrFljUSIR8QZ1rkxFRKqDx76bLyJSm6hMRUTcQGUqIuIGKlMRETdQmYqIuIHKVETEDVSmIiJuoDIVEXEDlamIiBuoTEVE3EBlKiLiBipTERE3UJmKiLiBylRExA1UpiIibqAyFRFxA5WpiIgbqExFRNxAZSoi4gYqUxERN1CZioi4gcpURMQNVKYiIm6gMhURcYN6VgcQqTWcxZCRCrmHoVFLCI0FH1+rU4mHUJmKVEbaKvhwOpw89Mu8wDYw9CnoOty6XOIxbIZhGFaHEPFoaavg7XHA+b8qNvO/Ma+pUEWvmYpUyFlsbpGWKVJ+mffhw+Y4qdNUpiIVyUgtvWtfhgEnD5rjpE5TmYpUJPewe8eJ11KZilSkUUv3jhOvpTIVqUhorPmu/dk3m8qwQWBbc5zUaSpTkYr4+JqHPwFlC/Xn20Of1PGmojIVuaCuw83DnwJbl54f2EaHRUkJHWcqUlk6A0oqoDIVEXED7eaLiLiBylRExA1UpiIibqAyFRFxA5WpiIUcDgcBAQFERUURFRXFzJkzrY4kVaTrmYpYrEePHnz22WdWx5BLpC1TkfPMmDGDzp07Ex8fT3x8PMnJyaxevZqYmBiioqKYMmUKTqcTh8NBr169GDt2LF26dGHq1Kkly3A1XrybylTkHFu2bGHdunXs2rWLxYsXs3nzZrKzs0lKSiI5OZnt27fj6+vLihUrAEhLSyMxMZHdu3ezadMm9u3bV+F4V9LS0rjyyiuJj49nx44dNbWq4mbazRc5R2pqKiNHjsTPz482bdowYMAA7HY7O3bsICYmBoD8/HxCQ0OJjo4mIiKC8PBwAHr27InD4WDv3r0ux7vSunVrDhw4QPPmzdmwYQO33nor6enpNbOy4lYqU5FzuDoh0DAMEhISWLhwYan5DocDu91ectvHx4fi4uJyx7tit9tLljFw4EDq169PdnY2wcHBl7gmUtO0my9yjtjYWFauXElRURFZWVmkpKRQUFDA2rVryczMBCAnJ6dk2pWYmJhKjz969CjFxeZHnuzevZvc3FyaN2/u5rWSmqAtU5Fz9OnTh7i4OLp160ZoaCiRkZGEhYWxYMECEhISKCwsxM/Pj0WLFhEUFORyGSEhIS7Ht2vXrszYDRs2MHPmTPz8/LDb7bz++uvYbOVdO1U8mS50InKevLw8AgICyM7Opm/fvmzbto3GjRtbHUs8nLZMRc4zYcIE0tPTKSwsZPbs2SpSqRRtmYrUgJycHAYPHlxqXnBwMGvWrLEokbibylRExA30br6IiBuoTEVE3EBlKiLiBipTERE3UJmKiLiBylRExA1UpiIibqAyFRFxA5WpiIgbqExFRNxAZSoi4gYqUxERN1CZioi4gcpURMQNVKYiIm6gMhURcQOVqYiIG6hMRUTcQGUqIuIG/w8Jh6k3twNdEgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(4, 4))\n", "\n", "grnplot(models[0], ax=ax, color=\"C1\")\n", "\n", "ax.set(title=\"cluster_1 (prototype1)\");" ] }, { "cell_type": "markdown", "id": "fb15ddc0-fe77-47d8-b653-32d7cc01ca6b", "metadata": {}, "source": [ "### edgecolor" ] }, { "cell_type": "code", "execution_count": 6, "id": "a2f87adc-50e5-4310-af8e-9d91a9f7685f", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVMAAAFeCAYAAAA41fyqAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAALqZJREFUeJzt3XtcVHXi//HXDMKAiCgiXgtRUPGKRgresIzUSlEzty1zU9O067Zba7Wul3Wtbc2NvmaW1rZu9uviJdMst7yQkaaWmSJqZoyFkQreERSY8/tjclZkQNSBM8D7+Xj4eHzmzGcO78Pgm3NmzhwshmEYiIjIVbGaHUBEpDpQmYqIeIDKVETEA1SmIiIeoDIVEfEAlamIiAeoTEVEPEBlKiLiASpTEREPUJl6SEpKChaLhZSUFLOjeLUHHniAxMREs2PUWHfeeScjRowwO0a1pDKtIjZu3Mi0adM4fvy4KV//9OnTTJ06lQEDBhASEoLFYuHf//73Za0jIyOD1157jaeffrpiQl7CmTNnmDZt2lX9wjP7eSjNzJkzGTx4MI0aNcJisTBt2jS38yZNmsTSpUv59ttvKzdgDaAyrSI2btzI9OnTTftPnJ2dzV//+ld2795N586dr2gdL774IhEREdxwww0eTlc+Z86cYfr06VddpmY+D6WZPHkyW7dupUuXLmXO69KlC7GxscyePbuSktUcKtMa7syZM+Wa16RJE7Kysjhw4ACzZs267K9TUFDAW2+95fFDzPLmr+4yMjLIyspi0aJFl5w7YsQIli1bxunTpyshWc2hMi2ngwcPMnbsWJo2bYrNZiMiIoKJEydy7ty5Uh/TokUL7r333hLL+/btS9++fYstmzNnDu3bt6d27drUr1+f2NhY/t//+38ATJs2jSeeeAKAiIgILBYLFosFu93uevyiRYu47rrrCAgIICQkhDvvvJOffvqpxNft0KEDX3/9NX369KF27drlPuS22Ww0bty4XHPdSU1NJTs7m5tuuqnY8vOvNb/77rs8/fTTNG7cmMDAQAYPHnxZ+Q8fPszYsWNp1KgR/v7+dO7cmYULF7oea7fbadiwIQDTp093fQ8vPBxet24dvXv3JjAwkHr16pGUlMTu3btd95f1PCQkJJS6x96mTRv69+/vymGxWHj++ed54YUXCA8PJyAggISEBNLS0ko8ds+ePQwfPpyQkBD8/f2JjY1lxYoVJea1aNGitG99CYmJieTm5vLpp5+W+zFyabXMDlAV/Pzzz3Tr1o3jx48zfvx42rZty8GDB1myZAlnzpzBz8/vqta/YMECHnnkEYYPH86jjz5Kfn4+O3bsYPPmzdx1110MGzaM7777jrfffpsXXniB0NBQAFc5zJw5k7/85S+MGDGC++67jyNHjjBnzhz69OnDN998Q7169VxfKycnh4EDB3LnnXcycuRIGjVqdFXZy2vjxo1YLJZSD0NnzpyJxWJh0qRJHD58mOTkZG666Sa2b99OQEBAmfnz8vLo27cv33//PQ899BAREREsXryYe++9l+PHj/Poo4/SsGFD5s2bx8SJExk6dCjDhg0DoFOnTgCsWbOGgQMH0rJlS6ZNm0ZeXh5z5syhZ8+ebNu2jRYtWpT5PNxzzz2MGzeOtLQ0OnTo4Mq7detWvvvuOyZPnlxse//zn/9w6tQpHnzwQfLz83nxxRe58cYb2blzp+s52bVrFz179qRZs2Y8+eSTBAYG8t577zFkyBCWLl3K0KFDr+i5aNeuHQEBAXzxxRdXvA5xw5BLGjVqlGG1Wo2tW7eWuM/hcBiGYRjr1683AGP9+vWu+8LDw43f/e53JR6TkJBgJCQkuG4nJSUZ7du3LzPDrFmzDMDIyMgottxutxs+Pj7GzJkziy3fuXOnUatWrWLLExISDMB45ZVXyvxal7J161YDMN54441yP2bkyJFGgwYNSiw//31r1qyZcfLkSdfy9957zwCMF1988ZL5k5OTDcBYtGiRa9m5c+eM+Ph4o06dOq71HjlyxACMqVOnlsgRExNjhIWFGTk5Oa5l3377rWG1Wo1Ro0a5lpX2PBw/ftzw9/c3Jk2aVGz5I488YgQGBhqnT582DMMwMjIyDMAICAgwMjMzXfM2b95sAMZjjz3mWtavXz+jY8eORn5+vmuZw+EwevToYURFRZXYhktt44Vat25tDBw4sMw5cnl0mH8JDoeD5cuXM2jQIGJjY0vcb7FYrvpr1KtXj8zMTLZu3XrZj122bBkOh4MRI0aQnZ3t+te4cWOioqJYv359sfk2m43Ro0dfdebLlZOTQ/369Uu9f9SoUQQFBbluDx8+nCZNmvDRRx8Vm+cu/0cffUTjxo357W9/61rm6+vLI488wunTp/nss8/KzJaVlcX27du59957CQkJcS3v1KkTiYmJJTK4ExwcTFJSEm+//TbGr9dbLyoq4t1332XIkCEEBgYWmz9kyBCaNWvmut2tWze6d+/u+lpHjx5l3bp1jBgxglOnTrme15ycHPr378++ffs4ePDgJXOVpn79+mRnZ1/x46UkleklHDlyhJMnTxY7dPO0SZMmUadOHbp160ZUVBQPPvggX3zxRbkeu2/fPgzDICoqioYNGxb7t3v3bg4fPlxsfrNmza76ZYkrZZTxRx2ioqKK3bZYLERGRhZ7XRjc5z9w4ABRUVFYrcV/nKOjo133l+X8/W3atClxX3R0NNnZ2eTm5pa5DnD+Qvjxxx/5/PPPAedLB4cOHeKee+4pMffi7QVo3bq1a3u///57DMPgL3/5S4nnderUqQAlntvLYRiGR3YE5H/0mmkFKu2HtaioCB8fH9ft6Oho9u7dy4cffsjq1atZunQpL7/8MlOmTGH69Ollfg2Hw4HFYuHjjz8uts7z6tSpU+z2ha8/VqYGDRpw7Nixq16PWfnLo3///jRq1IhFixbRp08fFi1aROPGjUu86VYeDocDgMcff9z15tXFIiMjrzjrsWPH3Ba6XDmV6SU0bNiQunXrun2n9VLq16/v9nzEAwcO0LJly2LLAgMD+c1vfsNvfvMbzp07x7Bhw5g5cyZPPfUU/v7+pRZzq1atMAyDiIgIWrdufdkZK0vbtm156623OHHiBMHBwSXu37dvX7HbhmHw/fffu94gKkt4eDg7duzA4XAU2zvds2eP634o/Zfb+fv37t1b4r49e/YQGhrqOkwva2/Ox8eHu+66i3//+98899xzLF++nHHjxrn9JXfx9gJ89913rnflz/98+Pr6XlEZl6WwsJCffvqJwYMHe3S9NZ0O8y/BarUyZMgQVq5cyVdffVXi/rIOXVu1asWXX35Z7PSpDz/8sMQpPzk5OcVu+/n50a5dOwzDoKCgAMD1n/nich42bBg+Pj5Mnz69RBbDMEqs2yzx8fEYhsHXX3/t9v7z726ft2TJErKyshg4cOAl133LLbfwyy+/8O6777qWFRYWMmfOHOrUqUNCQgIAtWvXBkp+D5s0aUJMTAwLFy4sdl9aWhqffPIJt9xyi2tZac/Deffccw/Hjh3j/vvv5/Tp04wcOdLtvOXLlxd7zXPLli1s3rzZtb1hYWH07duXV199laysrBKPP3LkSCnfjUtLT08nPz+fHj16XPE6pCTtmZbDM888wyeffEJCQgLjx48nOjqarKwsFi9eTGpqarFTjy503333sWTJEgYMGMCIESPYv38/ixYtolWrVsXm3XzzzTRu3JiePXvSqFEjdu/ezUsvvcStt97qelPmuuuuA+DPf/4zd955J76+vgwaNIhWrVrxt7/9jaeeegq73c6QIUMICgoiIyOD999/n/Hjx/P444975Pvw0ksvcfz4cX7++WcAVq5cSWZmJgAPP/yw2z3O83r16kWDBg1Ys2YNN954Y4n7Q0JC6NWrF6NHj+bQoUMkJycTGRnJuHHjLplr/PjxvPrqq9x77718/fXXtGjRgiVLlvDFF1+QnJzs+h4GBATQrl073n33XVq3bk1ISAgdOnSgQ4cOzJo1i4EDBxIfH8/YsWNdp0YFBwcXOxe1tOfhfMl26dKFDh06sHjxYqKjo+natavbzJGRkfTq1YuJEydy9uxZkpOTadCgAX/6059cc+bOnUuvXr3o2LEj48aNo2XLlhw6dIhNmzaRmZlZ7COhb775JgcOHHB9iGHDhg387W9/A5wFf37vG+DTTz+ldu3aukaCp5lzEkHVc+DAAWPUqFFGw4YNDZvNZrRs2dJ48MEHjbNnzxqG4f7UKMMwjNmzZxvNmjUzbDab0bNnT+Orr74qcWrUq6++avTp08do0KCBYbPZjFatWhlPPPGEceLEiWLrmjFjhtGsWTPDarWWOD1n6dKlRq9evYzAwEAjMDDQaNu2rfHggw8ae/fudc1JSEi45ClYZQkPDzcAt/8uPlXInUceecSIjIwstuz89+3tt982nnrqKSMsLMwICAgwbr31VuPAgQPF5paV/9ChQ8bo0aON0NBQw8/Pz+jYsaPbU7c2btxoXHfddYafn1+JU4jWrFlj9OzZ0wgICDDq1q1rDBo0yEhPTy+xjrKeB8MwjH/84x8GYDzzzDMlHnv+1KhZs2YZs2fPNq655hrDZrMZvXv3Nr799tsS8/fv32+MGjXKaNy4seHr62s0a9bMuO2224wlS5aU+N6U9txc/DPZvXt3Y+TIkW6/j3LlVKZSafbv32/4+voaa9ascS07X6aLFy82MZlnJScnGxaLpcQvA8MoXqZm+OabbwyLxWJ88803pnz96kyvmUqladmyJWPHjuXvf/+72VEqjGEYvP766yQkJHDttdeaHaeEv//97wwfPpyYmBizo1Q7es20hjt9+vQlL3jRsGFDt+9IX4l58+Z5ZD3eJjc3lxUrVrB+/Xp27tzJBx98YHYkt9555x2zI1RbKtMa7vnnn7/kuawZGRmXdSGNmujIkSPcdddd1KtXj6efflqnHdVAFsMo49weqfZ++OEHfvjhhzLn9OrVC39//0pKJFI1qUxFRDxAb0CJiHiAylRExANUpiIiHqAyFRHxAJWpiIgHqExFRDxAZSoi4gEqUxERD1CZioh4gMpURMQDVKYiIh6gMhUR8QCVqYiIB6hMRUQ8QGUqIuIBKlMREQ9QmYqIeIDKVETEA1SmIiIeoDIVEfEAlamIib744gtiY2Pp2LEj8fHx7Ny50+xIcoX010lFTLRr1y6Cg4Np3rw5a9asYebMmaxfv97sWHIFtGcqcpHJkyfTpk0bEhMTSUxMJCUlhVWrVhEXF0dMTAwTJkzA4XBgt9vp2rUrd999N23btmXixImudbib70779u1p3rw5ANdddx0HDx6slG0Uz1OZilxgy5YtrFu3jrS0NBYuXMjmzZvJzs4mOTmZlJQUtm/fjo+PD8uWLQMgPT2dadOmsWvXLjZt2sS+ffvKnF+W//znP/Tr16+iN1EqSC2zA4h4k40bNzJ06FB8fX1p2rQpvXv3xmazsWPHDuLi4gDIy8sjPDyc2NhYoqOjiYqKAqBTp07Y7Xb27Nnjdn5Zvv76a1555RVSU1MrdgOlwqhMRS7g7i0EwzBISkpi/vz5xZbb7XZsNpvrttVqpaioqNT5pcnMzOSuu+7ivffeo0GDBle3AWIaHeaLXKBHjx4sX76cwsJCsrKySE1NJT8/n7Vr15KZmQlATk6Oa+xOXFxcuefn5uaSlJTE7Nmz6dy5s+c3SCqNylTkAt27dychIYH27dvzu9/9js6dOxMZGcncuXNJSkqiU6dO3HzzzRw+fLjUdYSFhZV7/htvvMGePXuYPHkyMTEx9OjRo6I2TSqYTo0SuUhubi6BgYFkZ2cTHx/Ptm3bCAoKMjuWeDm9ZipykTFjxrB3714KCgqYMWOGq0iLHAZbMo5y+FQ+YUH+dIsIwcdqMTmteAvtmYqUw+q0LKavTCfrRL5rWZNgf6YOaseADk0u+ficnJwSpz2FhoayZs0aj2cVc6hMRS5hdVoWExdt4+L/KOf3SeeN7FquQpXqTW9AiZShyGEwfWW6q0i7WPZRhzMArmXTV6ZT5NA+SU2nMhUpw5aMo65D+2Yc4d9+z7HK72maW5zvzhtA1ol8tmQcNTGleAOVqUgZDp9yFqkPRbzoN5dgyxmOUpdfjBC386TmUpmKlCEsyB+A39daSqz1O04aATxc8BCFF50Ic36e1FwqU5EydIsI4bagfTzo8wEATxfcR6YR5rrfgvNd/W4RIaWsQWoKlalIGXzycnje5yWsFoN3C/vyoSPedd/5d/OnDmqn801FZSpSKocDlk/EP/8Ip4Na8krt8cXubhzsr9OixEXnmYqUZtNc+O/T4GODcesoCmuvT0BJqVSmIu78/A28lgiOArjleeg2zuxE4uV0mC9ysbOnYMkYZ5G2vQ2uv8/sRFIFqExFLrbqcTj6A9RtDoPngEWH8nJpKlORC21/G3a8AxYr3L4AauuUJykflanIednfw6o/Osd9n4JwXahZyk9lKgJQeBaWjIaCXGjRG3r/0exEUsWoTEUA1kyDX3ZAQAgMmw9WH7MTSRWjMhXZuxq+fNk5HvIy1G1qbh6pklSmUrOd/BmWT3SOu0+ENgPNzSNVlspUai5HESwbD3lHoXEnSJxudiKpwlSmUnN9/k+wfw6+gTD8DahlMzuRVGEqU6mZfvwSUp51jm99HkIjzc0jVZ7KVGqeM0dh6X1gFEHHEdD5t2YnkmpAZSo1i2HAiofhxE8Q0hJu+6c+LioeoTKVmuWr12HPh2D1heH/AluQ2YmkmlCZSs1xaBesfto5vmkaNO1iahypXlSmUjOcOwOLR0PRWYhMhLgHzE4k1YzKVGqG1U9C9l6o0wiGzAOrfvTFs/QTJdVf2lLYthCwOD93X6eh2YmkGlKZSvV2zA4rf+8c9/4DtOxrYhipzlSmUn0VFcCSsXD2JDTv5rxGqUgFUZlK9bV+Jhz8CmzBcPtr4ONrdiKpxlSmUj3tXwepyc7x4P+D+uGmxpHqT2Uq1c/pw7DsfsCA60ZD+yFmJ5IaQGUq1YvDAe9PgNzD0DAaBjxrdiKpIVSmUr1segn2r4Va/nDHG+AbYHYiqSFUplJ9HPwa1v56gecBz0JYtLl5pEZRmUr1kH8SlowBRyFED3a+VipSiVSmUvUZBnz4mPME/eBrne/e67J6UslUplL1bX8L0paAxQeGvw4B9c1OJDWQylSqtiPfwUdPOMc3PA3XdDM3j9RYKtNfHThwgK5duxITE0Pnzp1ZsWKF2ZHkUgryna+TFpyBiD7Q6zGzE0kNZjEMwzA7hDc4d+4chmFgs9k4dOgQ119/PT/++KPZsaQsH/0JtrwKtUNhQirUbWJ2IqnBvHrPdPLkybRp04bExEQSExNJSUlh1apVxMXFERMTw4QJE3A4HNjtdrp27crdd99N27ZtmThxomsd7ua74+fnh83m/FO/+fn5pc4TL7FnlbNIAYa+oiIV03ltmW7ZsoV169aRlpbGwoUL2bx5M9nZ2SQnJ5OSksL27dvx8fFh2bJlAKSnpzNt2jR27drFpk2b2LdvX5nz3dm/fz8dO3akffv2vPTSS5W1qXK5ThyEDx50juMfgqhEc/OIALXMDlCajRs3MnToUHx9fWnatCm9e/fGZrOxY8cO4uLiAMjLyyM8PJzY2Fiio6OJiooCoFOnTtjtdvbs2eN2fmlatWrFzp07+eGHHxg1ahQDBgzA39+/4jdWys9RBMvGQd4xaBID/aaanUgE8OIydfdSrmEYJCUlMX/+/GLL7Xa76xAdwGq1UlRUVGJ+kcNgS8ZRPth+kLAgf7pFhOBjLXk+YsuWLQkJCWHnzp1cf/31Ht4yuSobZsGBL8CvjvOvi9byMzuRCODFh/k9evRg+fLlFBYWkpWVRWpqKvn5+axdu5bMzEwAcnJyXGN34uLiXPNXp2XRfer73DH7Qx59Zzu/XfAlvZ5bx+q0LAAOHjxIXl4eAL/88gs7duwgIiKi4jdUys/+BXz2nHN82wvQoJW5eUQu4LV7pt27dychIYH27dsTHh5O586diYyMZO7cuSQlJVFQUICvry8LFiwgJCTE7TrCwsKYO3cuN9x8CweyT2Gx+hAy4GFq1Q0F4JcT+UxctI15I7vik5XOH/7wB6xWK1arldmzZxMaGlqZmyxlOXMUlt4HhgM63wWdRpidSKQYrz41Kjc3l8DAQLKzs4mPj2fbtm0EBQVd1jqKHAa9nltH1ol8t/dbgMbB/qROutHtIb94AcOAd+6CvR9Bg0gY/xnY6pidSqQYr90zBRgzZgx79+6loKCAGTNmXHaRAmzJOFpqkQIYQNaJfLZkHCW+VYOrSCsVZssCZ5H6+DlfJ1WRihfy6jJ99913r3odh08VL9KivJMceufPxZb5BNTl8J2rrvprSQX4ZSd8Mtk5TvwrNOlsbh6RUnh1mXpCWFDxU5t8AurSdPScS84TL3AuFxaPhqKz0HoAdJ9gdiKRUnntu/me0i0ihCbB/pT2aqgFaBLsPE1KvMxHf4KcfRDUBJJe1mX1xKtV+zL1sVqYOqgdQIlCPX976qB2evPJ2+xcAtsXARYYtgAC9Xq2eLdqX6YAAzo0Yd7IrjQOLn4o3zjYn3kjuzKggz7X7VWO/gArf+8c93kCInqbGkekPLz61ChPO/8JqMOn8sv8BJSYqPAc/Ks//LwNro2H330IPtX+pX2pBmrUT6mP1aLTn7zduhnOIvWv5zy8V5FKFVEjDvOlivh+DWz8P+c4aS7Uu8bcPCKXQWUq3uHUIXj/11Ofrr8Pom8zN4/IZVKZivkcDnh/POQegbD2cPPfzE4kctlUpmK+jS/CDylQKwDueAN8A8xOJHLZVKZirsyvYN2ve6K3/AMatjE3j8gVUpmKefKOw5LR4CiE9sOgyz1mJxK5YipTMYdhwIe/h+M/Qr1wGJSsj4tKlaYyFXNs+w/seh+stZyX1fMPNjuRyFVRmUrlO7wHPp7kHN84GZrHmptHxANUplK5CvJgyRgozIOWN0CPR81OJOIRKlOpXJ9MhsO7ILAhDH0VrPoRlOpBP8lSedJXwNbXnOOhr0JQI3PziHiQylQqx/GfYMVDznGPRyCyn7l5RDxMZSoVr6jQ+Wea809As+vgxr+YnUjE41SmUvE+ew5++hJsdeH216GWn9mJRDxOZSoVK2MDbJjlHN/2AoREmJtHpIKoTKXi5ObAsvGAAV1GQsfhZicSqTAqU6kYhgEfPACnsiC0NQz8h9mJRCqUylQqxuZX4LvV4GNzflzUL9DsRCIVSmUqnpf1LXw6xTnuPxMadzQ3j0glUJmKZ509BYtHQ9E5aHub80+QiNQAKlPxrI+egKP7oW4zGDxHl9WTGkNlKp7z7bvw7dtgscLtr0HtELMTiVQalal4Rs5+WPUH5zjhSQjvYW4ekUqmMpWrV3jOeVm9c6chvBf0edzsRCKVTmUqV2/tdMjaDgH1Ydh8sPqYnUik0qlM5ep89wlsesk5TnoZgpuZm0fEJCpTuXIns2D5BOe42/3Q9hZz84iYSGUqV8ZRBO+PhzM5zpPyE/9qdiIRU6lM5cqkvuC8IpRvIAx/A3z9zU4kYiqVqVy+HzfD+mec41tmQWiUuXlEvIDKVC5P3jFYOhaMIuh4B8TcZXYiEa+gMpXyMwxY8Qic+AnqR8Ct/9THRUV+pTKV8vv6Ddi9Aqy+zsvq+dc1O5GI11CZSvkcSofVTznHN02FZl3NzSPiZVSmcmnnzjg/LlqYD5E3QdyDZicS8ToqU7m0/z4FR3ZDnUYw5BWw6sdG5GL6XyFl2/U+fP1vwAJDX4U6Dc1OJOKVVKZSumMHYMWjznGvx6DVDebmEfFiKlNxr6gAlt4HZ09A8+vhhqfNTiTi1VSm4l7Ks5C5BWzBcPvr4ONrdiIRr6YylZJ+SIHP/+kcD34R6oebGkekKlCZSnGnj8Cy8YAB190L7YeanUikSlCZyv84HLB8Ipw+BA2jof+zZicSqTJUpvI/X74M338KtfydHxf1q212IpEqQ2UqTge3wZppznH/Z6BRO1PjiFQ1KlOB/JPOj4s6CiB6MMSOMTuRSJWjMq3pDANW/RGOZUDwNTD4/3RZPZEroDKt6b59G3a+BxYf5/mkAfXNTiRSJalMa7Ls72HV487xDU/Btd3NzSNShalMa6rCs7DkXijIhRa9odcfzE4kUqWpTGuqT6fALzuhdgMYtgCsPmYnEqnSVKY10d6PYfMrzvGQeVC3ibl5RKoBlWlNc/JnWP6Acxz3ILTub24ekWpCZVqTOIpg6TjIOwpNOjv/lpOIeITKtCb5fDYcSAW/OjD8DahlMzuRSLWhMq0pDmx0XqMU4NbZ0KCVuXlEqhmVaU1w5qjz8N5wQOffQuc7zU4kUu2oTKs7w4AVD8PJTAhpBbc8b3YikWpJZVrdbX0N9nwIVl/nZfVsdcxOJFItqUyrs1/S4L9/do4T/wpNY0yNI1KdqUyrq3O5sGQ0FJ2FqP4QN9HsRCLVmsq0uvp4EmR/B0FNYMjLuqyeSAWrZXYAqQBpS+GbNwELDJsPgaFmJxIpl3OFDt7cZOfA0TOEh9TmnvgW+NWqGvt8FsMwDLNDiAcds8MrveHsSejzBNw42exEIuXy7EfpLPg8A8cFjWS1wLjeETx1i/f/GR3tmVYnRQXOPz9y9iRcEwcJT5qdSKRcnv0onVc3ZJRY7jBwLff2Qq0a+89SPutmwMGvwT8Ybn8NfPS7UrzfuUIHCz4vWaQXWvB5BucKHRWaIy8vj379+lGnTh2efPLyd0RUptXF92vhixed46S5UO8ac/OIlNObm+zFDu3dcRjOeRWpVq1aTJkyhVmzZl3R41Wm1cHpw/D+BOc4dixEDzI3j8hlOHD0TKn3HdvwJgcX3M+hdyYz8+G7SUlJYdWqVcTFxRETE8OECRNwOBzY7Xa6du3K3XffTdu2bZk48X+nArqb746vry8JCQkEBARc0XaoTKs6hwPevx9yD0NYe+g/0+xEIpclPKS22+Vnf97L2QPf0nTMXBrc+hhZ+9LIzs4mOTmZlJQUtm/fjo+PD8uWLQMgPT2dadOmsWvXLjZt2sS+ffvKnO9pelGtqts0B/avg1oBzo+L+l7Zb1URs9wT34KZH+0ucah/9uAeAlrHY/GphV/dBvRN6I3NZmPHjh3ExcUBztc5w8PDiY2NJTo6mqioKAA6deqE3W5nz549budXBJVpVZb5Faz9q3M88DkIa2tuHpEr4FfLyrjeEW7ezf9fu47rHUHqDguGYZCUlMT8+fOLzbTb7dhs/7s+r9VqpaioqNT5FUGH+VVV/gnnaVCOQmg/FLqOMjuRyBV76pZ23N8nAusFH9SzNYsmb9+X3NfzGu7tUp/U1FTy8/NZu3YtmZmZAOTk5LjG7sTFxV3W/KuhPdOqyDDgw8fg+AGody0MelEfF5Uq76lb2vHHm9te8AmoaL5reoQ3H7+DDeHhdO7cmcjISObOnUtSUhIFBQX4+vqyYMECQkJC3K4zLCzM7fzmzZu7nd+pUycyMzMpLCxk0aJFfPPNNzRs2LBc+fUJqKpo25uw4iGw+MCY/8I115udSKRC5ObmEhgYSHZ2NvHx8Wzbto2goCCzY7mlPdOq5she+PhPzvGNk1WkUq2NGTOGvXv3UlBQwIwZM7y2SEF7plVLQT681g8OpUHLvjDyfbDqZW+Ry5GTk0O/fv2KLQsNDWXNmjVXtV6VaVWy6nHYugACG8KELyCokdmJRORX2q2pKnZ/6CxSgKGvqEhFvIzKtCo4kQkfPOgc93gYIm8yN4+IlKAy9XZFhc4/05x/HJp2hRunmJ1IRNxQmXq7Df+AHzeCXxAMfx1q+ZmdSETcUJl6M3sqbPj1cmCDkiGkpalxRKR0KlNvdeao8/DecEDMSOg43OxEIlIGlak3MgxY/gCc+hkaRMEt/zA7kYhcgsrUG22ZD999DD5+zsvq+QWanUhELkFl6m2yvoVPfv2LojfPhCadzM0jIuWiMvUmZ087L6tXdA7a3ArdxpmdSETKSWXqTT7+E+R8D3WbQdJLuqyeSBWiMvUWOxbD9rfAYoVhC6C2++szioh3Upl6g5z9zos9AyRMghY9zc0jIpdNZWq2wnOwdCycOwXhPaHPE2YnEpEroDI129rp8PM3EFDfeXhv9TE7kYhcAZWpmfZ9Cpteco6T5kJwM3PziMgVU5ma5dQv8P4E57jbeGh7q7l5ROSqqEzN4HDAsvFwJhsadYTEGWYnEpGrpDI1wxcvQMZn4Fsb7ngDfP3NTiQiV0llWtl+2gLrZjrHt8yC0Chz84iIR6hMK1PecVgyFowi6DAcYu42O5GIeIjKtLIYBqx8FE78CPVbwG0v6OOiItWIyrSybFsI6cvBWst5WT3/umYnEhEPUplWhsO74eNJznG/qdDsOnPziIjHqUwrWkEeLB4NhfnQqh/EP2R2IhGpACrTivbfp+HIbggMg6GvgFXfcpHqSP+zK1L6B/DVv5zjYa9CnTBz84hIhVGZVpTjP8KKh53jXo9BqxvNzSMiFUplWhGKCmHpfZB/AppfDzf82exEIlLBVKYVIeVZ+Gkz2OrC7a+Bj6/ZiUSkgqlMPe2Hz+Dz2c7xoBedJ+iLSLWnMvWk3Gzn1aAwoOso6DDM7EQiUklUpp5iGLB8Ipz+BULbwIDnzE4kIpVIZeopX74M+z6BWv7Oy+r51TY7kYhUIpWpJ/z8DXw61TnuPxMatTc3j4hUOpXp1Tp7CpaMAUcBtL0NYseanUhETKAyvVqrHoejP0Dd5jB4ji6rJ1JDqUyvxva3Ycc7YPGB4a9D7RCzE4mISVSmVyr7e1j1R+e471NwbZy5eUTEVCrTK1F4FpaMhoJcaNEbev/B7EQiYjKV6ZVYMw1+2QEBITBsPlh9zE4kIiZTmV6u7/7rPKcUnNcnrdvU3Dwi4hVUppfj5M/OTzkBxD0Arfubm0dEvIbKtLwcRc7P3Z/Jgcad4KZpZicSES+iMi2vz/8J9s/BNxCGvwG1bGYnEhEvojItjx+/dF6jFODW2RAaaW4eEfE6KtNLyTvmvGq+UQSd7oSY35qdSES8kMq0LIbh/DtOJ36CkJZw6/NmJxIRL6UyLctXr8PulWD1heH/AluQ2YlExEupTEtzaBesfto5TpwOTbuYm0dEvJrK1J1zZ2DxaCg6C1E3O88pFREpg8rUndVPQvZeqNMYhszTZfVE5JJUphdLWwbbFgIWuH0BBIaanUhEqgCV6YWO2WHlo85x7z9CRB9T44hI1aEyPa+oAJaMhbMn4ZruzmuUioiUk8r0vPUz4eBX4B8Mt78GPrXMTiQiVYjKFGD/ekhNdo4Hz4F615oaR0SqHpXp6cPw/v2AAbFjoF2S2YlEpAqq2WXqcMD7E+D0IQhrB/2fMTuRiFRRNbtMN70E+9dCrQDnx0V9A8xOJCJVVI16l6XIYbAl4yiHT+UTcXYvHddOxwIw4FkIizY7nohUYTWmTFenZTF9ZTpZJ/KpwxlW+T2NxVpIVvMBNLnuXrPjiUgVVyMO81enZTFx0TayTuQDBjN9/0W49TCZRij9v7+d1bt+MTuiiFRx1b5MixwG01emY/x6+w6fz0jy2UihYeWRcw9xikCmr0ynyGGUuR4RkbJU+zLdknH01z1SaGU5yPRaCwH4Z+EdbDNaYwBZJ/LZknHUxJQiUtVV+zI9fCrfNZ5S601qW86SWtSeeUWDSp0nInK5qn2ZhgX5u8Z/LJjI8qIePFbwAMZFm37hPBGRy2UxDKNav1hY5DDo9dw6fjmRj7sNtQCNg/1JnXQjPlZdt1RErky13zP1sVqYOqgd4CzOC52/PXVQOxWpiFyVal+mAAM6NGHeyK40Di5+KN842J95I7syoEMTk5KJSHVR7Q/zL3ThJ6DCgvzpFhGiPVIR8YgasWd6no/VQnyrBiTFNCO+VQO3RZqVlUXdunV55ZVXTEgoIlVVjSrT8njyySdJTEw0O4aIVDFeXaaTJ0+mTZs2JCYmkpiYSEpKCqtWrSIuLo6YmBgmTJiAw+HAbrfTtWtX7r77btq2bcvEiRNd63A3vzQbNmygbt26dOzYsTI2T0SqEa8t0y1btrBu3TrS0tJYuHAhmzdvJjs7m+TkZFJSUti+fTs+Pj4sW7YMgPT0dKZNm8auXbvYtGkT+/btK3P+xQoLC5kyZQrTp0+vzM0UkWrCa68atXHjRoYOHYqvry9Nmzald+/e2Gw2duzYQVxcHAB5eXmEh4cTGxtLdHQ0UVFRAHTq1Am73c6ePXvczndn7ty53HHHHYSEhFTOBopIteK1ZeruJAPDMEhKSmL+/PnFltvtdmw2m+u21WqlqKio1PnubNmyhdTUVGbNmsXx48fx8fGhdu3ajBo16uo3RkSqPa89zO/RowfLly+nsLCQrKwsUlNTyc/PZ+3atWRmZgKQk5PjGrsTFxdX7vlvvfUWBw4cwG638/vf/56ZM2eqSEWk3Lx2z7R79+4kJCTQvn17wsPD6dy5M5GRkcydO5ekpCQKCgrw9fVlwYIFpR6ah4WFuZ3fvHnzSt4aEanuvPqk/dzcXAIDA8nOziY+Pp5t27YRFBRkdiwRkRK8ds8UYMyYMezdu5eCggJmzJihIhURr+XVe6YVIScnh379+hVbFhoaypo1a0xKJCLVQY0rUxGRiuC17+aLiFQlKlMREQ9QmYqIeIDKVETEA1SmIiIeoDIVEfEAlamIiAeoTEVEPEBlKiLiASpTEREPUJmKiHiAylRExANUpiIiHqAyFRHxAJWpiIgHqExFRDxAZSoi4gEqUxERD1CZioh4gMpURMQDVKYiIh6gMhUR8QCVqYiIB6hMRUQ8oJbZAUSqiiKHwZaMoxw+lU9YkD/dIkLwsVrMjiVeQmUqUg6r07KYvjKdrBP5rmVNgv2ZOqgdAzo0MTGZeAuLYRiG2SFEvNnqtCwmLtrGxf9Rzu+TzhvZVYUqes1UpCxFDoPpK9NLFCngWjZ9ZTpFDu2T1HQqU5EybMk4WuzQ/mIGkHUiny0ZRysvlHgllalIGQ6fKr1Ir2SeVF8qU5EyhAX5e3SeVF8qU5EydIsIoUmwP6WdAGXB+a5+t4iQyowlXkhlKlIGH6uFqYPaAZQo1PO3pw5qp/NNRWUqcikDOjRh3siuNA4ufijfONhfp0WJi84zFSknfQJKyqIyFRHxAB3mi4h4gMpURMQDVKYiIh6gMhUR8QCVqYiJ7HY7gYGBxMTEEBMTw5QpU8yOJFdI1zMVMVnHjh358ssvzY4hV0l7piIXmTx5Mm3atCExMZHExERSUlJYtWoVcXFxxMTEMGHCBBwOB3a7na5du3L33XfTtm1bJk6c6FqHu/lSvalMRS6wZcsW1q1bR1paGgsXLmTz5s1kZ2eTnJxMSkoK27dvx8fHh2XLlgGQnp7OtGnT2LVrF5s2bWLfvn1lzncnPT2dLl26kJiYyI4dOyprU8XDdJgvcoGNGzcydOhQfH19adq0Kb1798Zms7Fjxw7i4uIAyMvLIzw8nNjYWKKjo4mKigKgU6dO2O129uzZ43a+O02aNCEjI4MGDRqwYcMG7rjjDvbu3Vs5GysepTIVuYC7DwQahkFSUhLz588vttxut2Oz2Vy3rVYrRUVFpc53x2azudbRp08f/Pz8yM7OJjQ09Cq3RCqbDvNFLtCjRw+WL19OYWEhWVlZpKamkp+fz9q1a8nMzAQgJyfHNXYnLi6u3POPHDlCUVERALt27eL06dM0aNDAw1sllUF7piIX6N69OwkJCbRv357w8HA6d+5MZGQkc+fOJSkpiYKCAnx9fVmwYAEhIe6vYRoWFuZ2fvPmzUvM3bBhA1OmTMHX1xebzcabb76JxaKLp1RFutCJyEVyc3MJDAwkOzub+Ph4tm3bRlBQkNmxxMtpz1TkImPGjGHv3r0UFBQwY8YMFamUi/ZMRSpBTk4O/fr1K7YsNDSUNWvWmJRIPE1lKiLiAXo3X0TEA1SmIiIeoDIVEfEAlamIiAeoTEVEPEBlKiLiASpTEREPUJmKiHiAylRExANUpiIiHqAyFRHxAJWpiIgHqExFRDxAZSoi4gEqUxERD1CZioh4gMpURMQDVKYiIh6gMhUR8YD/Dzx5Y35SXAg1AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(4, 4))\n", "\n", "grnplot(models[0], ax=ax, edgecolor=\"C1\")\n", "\n", "ax.set(title=\"cluster_1 (prototype1)\");" ] }, { "cell_type": "markdown", "id": "8a9d79cf-a198-4c17-a751-0aa0479ec830", "metadata": {}, "source": [ "### kwargs regarding annotations" ] }, { "cell_type": "code", "execution_count": 7, "id": "f53cf3ec-3831-4f18-9006-aa1a5fc055ff", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABPgAAAFeCAYAAAAYOqHOAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAkEJJREFUeJzs3XlcVXX+x/EX+y6IiBuKLC7IoqkpaqbmklaKUjnNL0tbxtGaaZvG3PetmtKZSkunvWbGFrTFJUMlc0krU0BAccEVZXFn597z+4ORiUBTUw5w38/Hw8dDzv3ec9+H7MPlc8/5HDvDMAxERERERERERESkVrI3O4CIiIiIiIiIiIhcOzX4REREREREREREajE1+ERERERERERERGoxNfhERERERERERERqMTX4REREREREREREajE1+ERERERERERERGoxNfhERERERERERERqMTX4REREREREREREajE1+ERERERERERERGoxNfjqiO+//57u3bvj4eGBnZ0dO3fuNDtSBaNGjaJly5ZmxxARqRPeeecd7OzsyMjIKN/Wu3dvevfubVomEZEbraraJyJia1QL5VLU4KsDSkpKuPfeezl16hQLFizg/fffJzAw8Lq/zty5c1mxYsV136+ISE11/Phxpk+fXuM+NBERqWlUL0VEVAvFXI5mB5Dfbv/+/Rw6dIilS5fy6KOP3rDXmTt3Lvfccw9Dhw696ucuXboUq9V6/UOJiNxAx48fZ8aMGbRs2ZIOHTqYHeey1q5da3YEEbFhtaleiojcKKqFYiY1+OqArKwsAHx8fMwNchlOTk5mRxARMV1eXh4eHh43ZN/Ozs43ZL8iIiIiIlLz6RLdWm7UqFH06tULgHvvvRc7O7vyGUzr16+nZ8+eeHh44OPjQ0xMDKmpqRWeP336dOzs7Ni3bx+jRo3Cx8cHb29vHnroIfLz88vX2dnZkZeXx7vvvoudnR12dnaMGjUKgPPnz/PUU0/RsmVLXFxc8Pf3p3///uzYsaNCzp/P4Ovdu3f5fn7555133ilfd+bMGZ566imaN2+Oi4sLoaGhPP/88zobUKQOOHToEI899hht2rTBzc2NBg0acO+991aaJ3JxzsjmzZt55plnaNiwIR4eHgwbNozs7OwKa1u2bMldd93Fpk2b6NKlC66urgQHB/Pee+9Vev0DBw5w77334uvri7u7O9HR0axcubL88YSEBG6++WYAHnrooSpr1OVczP3NN9/w2GOP4e/vT0BAwFUdO8Du3bu57bbbcHNzIyAggNmzZ1dZA385g+9S81kSEhKws7MjISGhfFt6ejp33303jRs3xtXVlYCAAO677z7Onj17RccqIjdWXa6XhmHQp08fGjZsWP6hNUBxcTGRkZGEhISQl5d32X0sWrSI8PBwXFxcaNq0KY8//jhnzpyptG7btm3ccccd1K9fHw8PD6Kiovj73//+qxlFpGaoy7XwwIED2NnZsWDBgkqPbdmyBTs7O/79739fdh+qhaIz+Gq5P/7xjzRr1oy5c+fyxBNPcPPNN9OoUSPi4+MZNGgQwcHBTJ8+nYKCAl555RV69OjBjh07Kt3wYvjw4QQFBTFv3jx27NjBP//5T/z9/Xn++ecBeP/993n00Ufp0qULo0ePBiAkJASAMWPG8Mknn/CnP/2Jdu3akZuby6ZNm0hNTaVjx45V5p40aVKly4k/+OADvvrqK/z9/QHIz8+nV69eHDt2jD/+8Y+0aNGCLVu2MGHCBDIzM1m4cOF1/E6KSHX7/vvv2bJlC/fddx8BAQFkZGSwePFievfuTUpKCu7u7hXW//nPf6Z+/fpMmzaNjIwMFi5cyJ/+9CeWLVtWYd2+ffu45557eOSRRxg5ciRvvfUWo0aNolOnToSHhwNw8uRJunfvTn5+Pk888QQNGjTg3XffZciQIXzyyScMGzaMsLAwZs6cydSpUxk9ejQ9e/YEoHv37ld1nI899hgNGzZk6tSp5b+kXumxnzhxgj59+lBaWsr48ePx8PBgyZIluLm5XdP3vCrFxcXcfvvtFBUV8ec//5nGjRtz7NgxvvzyS86cOYO3t/d1ey0RuTZ1uV7a2dnx1ltvERUVxZgxY4iLiwNg2rRp7N69m4SEhMue+Tx9+nRmzJhBv379GDt2LHv27GHx4sV8//33bN68ufwqkq+//pq77rqLJk2a8OSTT9K4cWNSU1P58ssvefLJJ6/8P4aImKYu18Lg4GB69OjBhx9+yNNPP13hsQ8//BAvLy9iYmIu+XzVQgHAkFpvw4YNBmB8/PHH5ds6dOhg+Pv7G7m5ueXbdu3aZdjb2xsPPvhg+bZp06YZgPHwww9X2OewYcOMBg0aVNjm4eFhjBw5stLre3t7G48//vhlM44cOdIIDAy85OObN282nJycKuSYNWuW4eHhYezdu7fC2vHjxxsODg7G4cOHL/uaIlKz5efnV9q2detWAzDee++98m1vv/22ARj9+vUzrFZr+fann37acHBwMM6cOVO+LTAw0ACMjRs3lm/LysoyXFxcjL/85S/l25566ikDML799tvybefPnzeCgoKMli1bGhaLxTAMw/j+++8NwHj77bev+vgu5r7llluM0tLSazr2izm3bdtW4Xi8vb0NwDh48GD59l69ehm9evWq9Po/X2MY//uZsWHDBsMwDOOnn36q9DNERGqWul4vDcMw3njjDQMwPvjgA+O7774zHBwcjKeeeqrCml/WtaysLMPZ2dkYMGBAeQ7DMIxXX33VAIy33nrLMAzDKC0tNYKCgozAwEDj9OnTFfb58++TiNRsdb0WXqyDqamp5duKi4sNPz+/Cr+HqxbKpegS3TooMzOTnTt3MmrUKHx9fcu3R0VF0b9/f1atWlXpOWPGjKnwdc+ePcnNzeXcuXO/+no+Pj5s27aN48ePX1PeEydOcM8999ChQwcWLVpUvv3jjz+mZ8+e1K9fn5ycnPI//fr1w2KxsHHjxmt6PRGpGX5+FlpJSQm5ubmEhobi4+NT4RL/i0aPHo2dnV351z179sRisXDo0KEK69q1a1f+iSlAw4YNadOmDQcOHCjftmrVKrp06cItt9xSvs3T05PRo0eTkZFBSkrKdTlGgD/84Q84ODhU2Halx75q1Sqio6Pp0qVLheO5//77r1u+i2foffXVVxVGM4hIzWEL9XL06NHcfvvt/PnPf+aBBx4gJCSEuXPnXvY58fHxFBcX89RTT2Fv/79fa/7whz9Qr1698kvnfvrpJw4ePMhTTz1VaWb1z79PIlKz1fVaOHz4cFxdXfnwww/Lt3311Vfk5OQwYsSISz5PtVAuUoOvDrpYsNq0aVPpsbCwMHJycirNMmnRokWFr+vXrw/A6dOnf/X1XnjhBZKTk2nevDldunRh+vTpFYrh5ZSWljJ8+HAsFgtxcXG4uLiUP5aens6aNWto2LBhhT/9+vUDqDCnRURqn4KCAqZOnVo+Y9PPz4+GDRty5syZKme/XWmd+uW6i2t/vu7QoUOXrJEXH79egoKCKm270mM/dOgQrVq1qvT8qrL/lnzPPPMM//znP/Hz8+P222/ntdde0/w9kRrEVurlm2++SX5+Punp6bzzzju/Oo7gUu95nZ2dCQ4OLn98//79AERERFy3rCJS/ep6LfTx8WHw4MH861//Kt/24Ycf0qxZM2677bZLPk+1UC7SDD4BqHR2yUWGYfzqc4cPH07Pnj1Zvnw5a9eu5cUXX+T5558nLi6OQYMGXfa5f/3rX9m6dSvx8fHlw+cvslqt9O/fn3HjxlX53NatW/9qNhGpuf785z/z9ttv89RTT9GtWze8vb2xs7Pjvvvuq/ImEldap35LPbsRqvoF9WqP/Vpc6pNYi8VSadtLL73EqFGj+Oyzz1i7di1PPPEE8+bN47vvvqtUm0Wk+tlKvUxISKCoqAiApKQkunXrZkoOEamZbKEWPvjgg3z88cds2bKFyMhIPv/8cx577LEKZ+aJXIoafHVQYGAgAHv27Kn0WFpaGn5+fpcdVnwplzttt0mTJjz22GM89thjZGVl0bFjR+bMmXPZBt9//vMfFi5cyMKFC8vvBPxzISEhXLhwofyMPRGpWz755BNGjhzJSy+9VL6tsLCwyrt9XW+BgYGXrJEXH4cbd7nClR57YGAg6enplZ5fVfZfuvgp9S/3ealPmCMjI4mMjGTy5Mls2bKFHj168PrrrzN79uxffS0RubFsoV5mZmby5z//mQEDBuDs7Myzzz7L7bffXr7/S2WDspoYHBxcvr24uJiDBw+Wv4e8eGO45ORkva8UqcVsoRYOHDiQhg0b8uGHH9K1a1fy8/N54IEHfjUbqBaKLtGtk5o0aUKHDh149913KxS75ORk1q5dyx133HFN+/Xw8KhUPC0WS6XTof39/WnatGn5J7BVSU5O5tFHH2XEiBGXvFvP8OHD2bp1K1999VWlx86cOUNpaenVH4SI1BgODg6VPhl95ZVXqjzD7Hq744472L59O1u3bi3flpeXx5IlS2jZsiXt2rUDKP8w5Hq/cbzSY7/jjjv47rvv2L59e/m27OzsCrNZLuXim7ifzyu1WCwsWbKkwrpz585VqqeRkZHY29tfto6LSPWxhXr5hz/8AavVyptvvsmSJUtwdHTkkUceuewZNP369cPZ2Zl//OMfFda9+eabnD17ljvvvBOAjh07EhQUxMKFCyvlM+tsRRG5erZQCx0dHfn973/PRx99xDvvvENkZCRRUVGXfY5qoVykM/jqqBdffJFBgwbRrVs3HnnkEQoKCnjllVfw9vZm+vTp17TPTp06ER8fz8svv0zTpk0JCgqiTZs2BAQEcM8999C+fXs8PT2Jj4/n+++/r/DJyi899NBDANx666188MEHFR7r3r07wcHB/PWvf+Xzzz/nrrvuKr9NeV5eHklJSXzyySdkZGTg5+d3TcciIua76667eP/99/H29qZdu3bll+s3aNDghr/2+PHj+fe//82gQYN44okn8PX15d133+XgwYN8+umn5ZdBhISE4OPjw+uvv46XlxceHh507dq1yrl6V+NKj33cuHG8//77DBw4kCeffBIPDw+WLFlCYGAgiYmJl32N8PBwoqOjmTBhAqdOncLX15f//Oc/lZp569ev509/+hP33nsvrVu3prS0lPfffx8HBwfuvvvu33ScInJ91PV6+fbbb7Ny5Ureeeed8rEAr7zyCiNGjGDx4sU89thjVT6vYcOGTJgwgRkzZjBw4ECGDBnCnj17WLRoETfffHP5UHp7e3sWL17M4MGD6dChAw899BBNmjQhLS2N3bt3V/lhsojUPHW9Fl704IMP8o9//IMNGzbw/PPP/+p61UIpV/037pXrbcOGDQZgfPzxxxW2x8fHGz169DDc3NyMevXqGYMHDzZSUlIqrJk2bZoBGNnZ2RW2//LW24ZhGGlpacatt95quLm5GYAxcuRIo6ioyPjrX/9qtG/f3vDy8jI8PDyM9u3bG4sWLaqwv5EjRxqBgYHlX1+8HXlVf35+S/Hz588bEyZMMEJDQw1nZ2fDz8/P6N69u/G3v/3NKC4u/m3fOBEx1enTp42HHnrI8PPzMzw9PY3bb7/dSEtLMwIDA42RI0eWr7tYj77//vsKz79Y+zZs2FC+LTAw0LjzzjsrvVavXr2MXr16Vdi2f/9+45577jF8fHwMV1dXo0uXLsaXX35Z6bmfffaZ0a5dO8PR0bFSjbqcS+W+mmM3DMNITEw0evXqZbi6uhrNmjUzZs2aZbz55puVavSljrFfv36Gi4uL0ahRI2PixInG119/XeH7duDAAePhhx82QkJCDFdXV8PX19fo06ePER8ff0XHKSI3Xl2ul0eOHDG8vb2NwYMHV3ps2LBhhoeHh3HgwIEKx/fz2mcYhvHqq68abdu2NZycnIxGjRoZY8eONU6fPl1pf5s2bTL69+9f/p41KirKeOWVV341o4jUDHW5Fv5SeHi4YW9vbxw9erTSY6qFcil2hqFzMUVEREREREREaoKbbroJX19f1q1bZ3YUqUU0g09EREREREREpAb44Ycf2LlzJw8++KDZUaSW0Rl8IiIiV6GgoKDSzYV+ydfXF2dn52pKJCJSM6leiohceS3cu3cvP/74Iy+99BI5OTkcOHAAV1fXakopdYFusiEiInIVli1bVn6joEvZsGEDvXv3rp5AIiI1lOqliMiV18KEhARmzpxJmzZt+Pe//63mnlw1ncEnIiJyFTIzM9m9e/dl13Tq1In69etXUyIRkZpJ9VJERLVQqo8afCIiIiIiIiIiIrWYbrIhIiIiIiIiIiJSi6nBJyIiIiIiIiIiUoupwSciIiIiIiIiIlKLqcEnIiIiIiIiIiJSi6nBJyIiIiIiIiIiUoupwSciIiIiIiIiIlKLqcEnIiIiIiIiIiJSi6nBJyIiIiIiIiIiUoupwSciIiIiIiIiIlKLqcEnIiIiIiIiIiJSi6nBJyIiIiIiIiIiUoupwSciIiIiIiIiIlKLqcEnIiIiIiIiIiJSi6nBJyIiIiIiIiIiUoupwSdiknXr1vHwww/TunVr3N3dCQ4O5tFHHyUzM9PsaCIiIiIiIiJSi9gZhmGYHULEFnXu3JlTp05x77330qpVKw4cOMCrr76Ku7s7O3fupHHjxmZHFBEREREREZFaQA0+EZNs3LiRW265BXt7+wrbevXqxaRJk5g9e7aJ6URERERERESkttAlulJrJSQk0LlzZ1xdXQkJCeGNN95g+vTp2NnZVVj3wQcf0KlTJ9zc3PD19eW+++7jyJEjFdb07t2biIgIUlJS6NOnD+7u7jRr1owXXnih0usWFRUxbdo0QkNDcXFxoXnz5owbN46ioqKryn/rrbdWaO5d3Obr60tqaupV7UtEREREREREbJej2QFErsVPP/3EwIEDadKkCTNmzMBisTBz5kwaNmxYYd2cOXOYMmUKw4cP59FHHyU7O5tXXnmFW2+9lZ9++gkfH5/ytadPn2bgwIHExsYyfPhwPvnkE5577jkiIyMZNGgQAFarlSFDhrBp0yZGjx5NWFgYSUlJLFiwgL1797JixYrfdFwXLlzgwoUL+Pn5/ab9iIiIiIiIiIjtUINPaqVp06bh4ODA5s2badq0KQDDhw8nLCysfM2hQ4eYNm0as2fPZuLEieXbY2Njuemmm1i0aFGF7cePH+e9997jgQceAOCRRx4hMDCQN998s7zB969//Yv4+Hi++eYbbrnllvLnRkREMGbMGLZs2UL37t2v+bgWLlxIcXExv/vd7655HyIiIiIiIiJiW3SJrtQ6FouF+Ph4hg4dWt7cAwgNDS1vxAHExcVhtVoZPnw4OTk55X8aN25Mq1at2LBhQ4X9enp6MmLEiPKvnZ2d6dKlCwcOHCjf9vHHHxMWFkbbtm0r7PO2224DqLTPq7Fx40ZmzJjB8OHDy/cnIiIiIvJLmzdvpnPnzkRGRtKtWzeSkpLMjiQiIibTGXxS62RlZVFQUEBoaGilx36+LT09HcMwaNWqVZX7cXJyqvB1QEBApfl99evXJzExscI+U1NTK10K/PNs1yItLY1hw4YRERHBP//5z2vah4iIiIjYBh8fH1asWEFAQADx8fE88cQTv+mDZhERqf3U4JM6y2q1Ymdnx+rVq3FwcKj0uKenZ4Wvq1oD8PMbTVutViIjI3n55ZerXNu8efOrznnkyBEGDBiAt7c3q1atwsvL66r3ISIiIiLmmjx5Mh9//DEtWrQAYNKkSeTl5TFr1iwKCwuJjo5m0aJFHD58mNjYWMLCwvjxxx/p06cPixcvBmDlypWV1v/ypmwA4eHh5X/v1KkTx44dq56DFBGRGksNPql1/P39cXV1Zd++fZUe+/m2kJAQDMMgKCiI1q1bX5fXDgkJYdeuXfTt27fS2X7XIjc3lwEDBlBUVMS6deto0qTJdUgpIiIiItVp+/btrF+/nuTkZLKzs8vHubzxxhskJCTg6urK448/TlxcHJ07dyYlJYVly5YRHBxMp06dSE9Pp379+ixcuLDS+nvuueeyr/3ee+/Rt2/fajpSERGpqTSDT2odBwcH+vXrx4oVKzh+/Hj59n379rF69eryr2NjY3FwcGDGjBkVzsKDsrPycnNzr/q1hw8fzrFjx1i6dGmlxwoKCsjLy7vifeXl5XHHHXdw7NgxVq1adclLiUVERESkZtuyZQvDhg3DycmJpk2b0rNnT1xcXEhMTCQ6OpoOHToQHx9fPts5LCyMVq1a4eDgQFRUFBkZGWzduvWS6y/lxx9/5PXXX2f27NnVcZgiIlKD6Qw+qZWmT5/O2rVr6dGjB2PHjsVisfDqq68SERHBzp07gbKz7WbPns2ECRPIyMhg6NCheHl5cfDgQZYvX87o0aN59tlnr+p1H3jgAT766CPGjBnDhg0b6NGjBxaLhbS0ND766CO++uorOnfufEX7uv/++9m+fTsPP/wwqamppKamlj/m6enJ0KFDryqbiIiIiJjjlx8mX9wWExPDkiVLKmzPyMjAxcWl/Gt7e3ssFssl11/K0aNH+b//+z8++ugjGjRo8NsOQEREaj01+KRW6tSpE6tXr+bZZ59lypQpNG/enJkzZ5KamkpaWlr5uvHjx9O6dWsWLFjAjBkzgLI5eQMGDGDIkCFX/br29vasWLGCBQsW8N5777F8+XLc3d0JDg7mySefvKpLgS82It966y3eeuutCo8FBgaqwSciIiJSS3Tv3p1nnnmGp59+muzsbDZt2sTIkSNZt24dR48eJSAggNzcXAoKCi65j+joaJ5++ulK6wMCAiqtzcvLIyYmhpdeeon27dvfyEMTEZFaws6o6uMmkVpq6NCh7N69m/T0dLOjiIiIiIgNmThxIp9++imBgYEUFhaycOFCsrKymDRpEiUlJTg5ObF06VJ8fX257777+O677wAYNWoU9913HwMHDmTNmjWV1nfs2LHSa7366qs899xz5SNe3N3d2bJlS7Uer4iI1Cxq8EmtVVBQgJubW/nX6enphIeHM3LkyCpn5ImIiIiI3Ch5eXl4eHiQk5NDt27d2LFjB15eXmbHEhERG6FLdKXWCg4OZtSoUbRsGcTWxD3EffgWjk7O/OXZv5qa6+zZs5e9/AKgcePG1ZRGRGyJxWqw/eApss4X4u/lSpcgXxzsf/sdv0VEahOzauHDDz/Mnj17KCkpYdasWWruiYjp9N7QtugMPqm1HnroIVatjSc76yTYO+HSrC0+tz5IyzYRTBvcjoERTUzJNWrUKN59993LrtH/diJyva1JzmTGFylkni0s39bE29XUeigiUt3qWi3Mzc2lb9++Fbb5+fkRHx9vUiIRqS3qWj2UX6cGn9Raa5IzGfvBDn75D/ji5xGLR3Q0pXClpKRw/Pjxy67p169fNaUREVtQU+uhiEh1Ui0UESmjemib1OCTWsliNbjl+fUVPo34OTugsbcrm567Tacgi0idpnooIqJaKCJykeqh7bI3O4DItdh+8NQlCxaAAWSeLWT7wVPVF0pExASqhyIiqoUiIhepHtouNfikVso6f+mCdS3rRERqK9VDERHVQhGRi1QPbZcafFIr+Xu5Xtd1IiK1leqhiIhqoYjIRaqHtksNPqmVugT50sT70gXJjrI7BHUJ8q2+UCIiJugS5Evjei6XfFz1UERsgd4bioiU6RLkSyO9N7RJavBJreRgb8cz/VtV+djFMaHTBrfT0FARqfMc7O1o39ynysdUD0XEVjjY2/FE39AqH1MtFBFb4mBvR7sm9ap8TPWwblODT2qtzftyAXD8RWFq7O2q236LiM1IOX6O+NQsAOq7O1V4TPVQRGzJhrRsAJwc9N5QRGzXjsOnSdhbVg999d7QpjiaHUDkWnyzN5sVO49jbwcfj+lGYYmVrPOF+HuVnWqsTyNExBZYrAYT4hKxWA0GRTTm1f/ryPaDp1QPRcTmrEk+wdqUkzja2/HZ4z04W1CqWigiNqfEYmXCp0kYBsR2bMaL97TXe0Mbogaf1Dr5xaVMWp4EwKjuQdzUor7JiUREzPHulgx2HT2Ll6sjM4aE42BvR7eQBmbHEhGpVucKS5j6WTIAY3qF0K6pt8mJRETMsWTjAfacPI+vhzOT72yn94Y2RpfoSq2z4Ou9HD1dQDMfN/4yoLXZcURETHH0dD5/W7sHgPGD2uJfT3dCExHb9PzqNLLOFxHk58Gfbqt6Dp+IXLlDhw7RsWNHOnToQPv27fn888/NjiRX4ED2Bf6+Lh2AKXeF4evhbHIiqW46g09qlaSjZ3lz00EAZg+NwMNF/4RFxPYYhsGUFcnkF1vo0tKX39/cwuxIIiKm+D7jFB9uOwzA3GGRuDo5mJxIpPZr0qQJW7duxcXFhZMnT3LzzTczZMgQs2PJZRiGwcTlSRSXWrm1dUOGdmhmdiQxgc7gk1qj1GJlfFwiVgMGt29Kn7b+ZkcSETHFF4mZbNiTjbODPXNjI7HXLBURsUFFpRbGf5oIwO86N9dlaFKnTZ48mTZt2tC/f3/69+9PQkICK1euJDo6mg4dOjBmzBisVisZGRl07NiR+++/n7Zt2zJ27NjyfVS1virOzs64uLgAUFhYeMl1UnN89MMRvjtwCjcnB+YMjcDOTu8NbZEafFJrvLX5ILuPn8PbzYmpd7UzO46IiCnO5Bcz84vdADzeJ5RQf0+TE4mImGNxwn72Z+fh5+nCxDvCzI4jcsNs376d9evXk5yczLvvvsu2bdvIyclh4cKFJCQksHPnThwcHIiLiwMgJSWF6dOns3v3brZu3Up6evpl11dl//79REZGEh4ezquvvlpdhyrXIPt8EXNWpgLwTP/WNPd1NzmRmEXXN0qtcDg3n5e/3gvApDvDaOjlYnIiERFzzF2VSs6FYlr5ezK2d4jZcURETLEv6zyLNuwHYPqQdni7O5mcSOTG2bJlC8OGDcPJyYmmTZvSs2dPXFxcSExMJDo6GoCCggICAwPp3LkzYWFhtGrVCoCoqCgyMjJIS0urcv2lhISEkJSUxIEDB3jwwQcZOHAgrq6a91sTzfhiN+cKS4loVo+HerQ0O46YSA0+qfEMw2DSiiQKS6x0C27AvZ0CzI4kImKKLftz+OiHowDMi43E2VEn4ouI7bFaDcZ/mkSxxcptbf25M7KJ2ZFEbijDMKrcFhMTw5IlSypsz8jIKL+8FsDe3h6LxXLJ9b8mODgYX19fkpKSuPnmm6/tAOSGWZ92ki8TM3Gwt2N+bBSODnpvaMv0X19qvOU/HePb9BycHctmTWmegIjYosISCxPjkgAYEd2Czi19TU4kImKOf20/zA+HTuPh7MAszZoSG9C9e3dWrFhBaWkpmZmZbNq0icLCQtatW8fRo2Uf/OXm5pb/vSrR0dFXvP7YsWMUFBQAcOLECRITEwkKCrrORyW/1YWiUiYvTwbgkVuCiGjmbXIiMZvO4JMaLfdCEbO+TAHgyb6tCPLzMDmRiIg5/rEunYzcfBrVc2HcwLZmxxERMcWJs4U8vzoNgGdvb0MzHzeTE4nceF27dqVXr16Eh4cTGBhI+/btCQ0N5bXXXiMmJoaSkhKcnJxYunQpvr5VfwDo7+9f5fqAgMpXR6WkpPDMM89gb2+Pvb09L730En5+fjf6MOUq/e2rPRw/W0hzXzee7tfa7DhSA9gZVZ3vK1JDPLNsJ3E/HaNtYy+++PMtOOmUYxGxQamZ5xj8yiZKrQZvPNCJ28Mbmx1JRMQUY97/kTW7T9C+uQ9xY7vjoLuIi43Iy8vDw8ODnJwcunXrxo4dO/Dy8jI7lphk55EzDFu0GcOA9x7uwq2tG5odSWoAncEnNdbGvdnE/XQMO7uyWVNq7omILbJYDcbHJVFqNbg9vJGaeyJis77afYI1u0/gaG/H/NhINffEpjz88MPs2bOHkpISZs2apeaeDSuxWBn/aSKGAcNuaqbmnpRTg09qpPziUiatKJs1NbJbS25qUd/kRCIi5nhvawa7jpzBy8WRmTERZscRETHFucISpn5WNmtq9K3BhDWpZ3Iikeq1bNmy677P3Nxc+vbtW2Gbn58f8fHx1/215PpZsvEAaSfOU9/dicl3hpkdR2oQNfikRloYn86RUwU09Xbl2dvbmB1HRMQUx84U8OJXewB4blBbGtVzNTmRiIg5XliTxslzRbRs4M4TfVuZHUekTmjQoAE7d+40O4ZchYM5efx9XToAU+5qRwNPl195htgSXfMoNU7ysbP889sDAMweFoGni/rQImJ7DMNgyopk8ost3NyyPv/XpYXZkURETPFDxik++O4wAHNjI3F1cjA5kYhI9TMMg4lxSRSXWunZyo9hNzUzO5LUMGrwSY1SarEyPi4RqwF3RTXhtraNzI4kImKKlUmZrE/LwsnBjnmxkdhr1pSI2KCiUgvj48rGtgzvHED3EN3JU0Rs08c/HmXrgVxcneyZMzQSOzu9N5SK1OCTGuXtzRkkHztHPVdHpg5uZ3YcERFTnM0vYfrnKQA81juUUH8N0hYR2/R6wgH2ZV3Az9OZiXdo1pSI2Kbs80XMWZkKwNP9WtOigbvJiaQmUoNPaowjp/J5+eu9AEy6Mwx/L82aEhHbNHdVKjkXigj19+SxPiFmxxERMcW+rPO8tmEfAFMHh+Pj7mxyIhERc8z8MoWzBSWEN63HI7cEmR1Haig1+KRGMAyDicuTKCixEB3sy/DOzc2OJCJiiq37c1n2wxEA5sVG4uKoWVMiYnusVoMJcUkUW6z0adOQwVFNzI4kImKKDWlZfLHrOPZ2MD82CkcHtXGkavqXITXCZzuP8216Ds6O9syLjdI8ARGxSYUlFiYuL5s1dX/XFtzc0tfkRCIi5vjP90f4PuM07s4OzB6mWVMiYpvyikqZvCIZgEduCSIywNvkRFKTqcEnpjuVV8zML8tmTT3ZtxVBfh4mJxIRMcer6/dxMCcPfy8XnhvU1uw4IiKmyDpXyLzVZbOmnh3QhmY+biYnEhExx0tr93LsTAEB9d14un9rs+NIDacGn5hu9soUTuUV06aRF3/oGWx2HBERU6SdOMfr3+wHYGZMOPVcnUxOJCJijmmf7+Z8YSntA7wZ2b2l2XFEREyx68gZ3tlyEIA5wyJxd3Y0OZHUdGrwiam+Tc8mbscx7Oxg/t2RODvqn6SI2B6L1WD8p0mUWg0GtGvEwAjNmhIR27R29wlWJ5/Awd6OebFRONjr0lwRsT0lFivPfZqI1YChHZrSq3VDsyNJLaBuipimoNjCpOVl8wRGdmvJTS3qm5xIRMQc72/NYOeRM3i5ODIzJsLsOCIipjhfWMLUz3YDMPrWYNo1rWdyIhERcyz99gBpJ85T392JKXe1MzuO1BJq8IlpFq7by+FT+TT1duXZ29uYHUdExBTHzxTw4ld7ABg3qC2NvV1NTiQiYo4Xv9rDiXOFtGzgzpN9W5kdR0TEFBk5efw9Ph2AyXe2o4Gni8mJpLZQg09MkXzsLP/8tmyewMyYCDxdNE9ARGyPYRhM/SyZvGILnQLrc3+XFmZHEhExxY+HTvP+d4cAmDssElcnB5MTiYhUP8MwmLQiiaJSK7eE+hHbsZnZkaQWUYNPql2pxcqEuCQsVoM7I5vQr10jsyOJiJhiVdIJ4lOzcHKwY35sJPaaNSUiNqi41MqEuEQMA+7pFED3UD+zI4mImOKTH4+yeV8urk72zBkWgZ2d3hvKlVOD7zrbuHEjQ4YMoXnz5ri6utK4cWMGDhzI5s2bzY5WY7yzJYOkY2ep5+rItCGaJyAitulsfgnTPi+bNTW2dyitGnmZnEhExByvf7OfvScv0MDDmUl3hJkdR0TEFDkXipizKhWAp/q1JrCBh8mJpLbRdZHX2d69e7G3t2fMmDE0btyY06dP88EHH3DrrbeycuVKBg4caHZEUx05lc9La/cCMPGOMPy9NGtKRGzTvNWp5FwoIrihB4/3CTE7joiIKfZlXeDV9fsAmDq4HfU9nE1OJCJijplfpHAmv4R2Terx6C1BZseRWsjOMAzD7BB1XX5+PsHBwXTo0IE1a9aYHcc0hmEw6u3v+WZvNl2DfPnP6GidciwiNum7A7nct+Q7AD76Yze6BPmanEhEpPpZrQb3LfmO7Rmn6N2mIW+PulnvDUXEJm3Yk8VDb3+PvR2seLwHUQE+ZkeSWqjGXaKbkJBA586dcXV1JSQkhDfeeIPp06dX+mH/wQcf0KlTJ9zc3PD19eW+++7jyJEjFdb07t2biIgIUlJS6NOnD+7u7jRr1owXXnih0usWFRUxbdo0QkNDcXFxoXnz5owbN46ioqLffEzu7u40bNiQM2fO/OZ91Waf7zrON3uzcXa0Z25spN7AiYhNKiyxMDEuCYDfd2mh5p6I2KxlPxxhe8Yp3JwcmBWjWVMiYpvyikqZvDwZgId6BKm5J9esRl2i+9NPPzFw4ECaNGnCjBkzsFgszJw5k4YNG1ZYN2fOHKZMmcLw4cN59NFHyc7O5pVXXuHWW2/lp59+wsfHp3zt6dOnGThwILGxsQwfPpxPPvmE5557jsjISAYNGgSA1WplyJAhbNq0idGjRxMWFkZSUhILFixg7969rFix4qqP5dy5cxQXF5OTk8N7771HcnIyEydO/C3fnlrtdF4xM79IAeDPfUIJaehpciIREXO8tmEfB3Ly8PdyYfygtmbHERExRda5Qub+d9bUXwa0prmvu8mJRETM8fLXezl2poBmPm4807+12XGkFqtRDb5p06bh4ODA5s2badq0KQDDhw8nLOx/w3YPHTrEtGnTmD17doWGWWxsLDfddBOLFi2qsP348eO89957PPDAAwA88sgjBAYG8uabb5Y3+P71r38RHx/PN998wy233FL+3IiICMaMGcOWLVvo3r37VR3L8OHD+eqrrwBwdnbmj3/8I1OmTLnK70jdMXtlKrl5xbRp5MUfe2nWlIjYpj0nzrM4YT8AM4aE4+3mZHIiEalpDh06xLBhw7BarRiGwaxZsxgyZIjZsa676V/s5nxhKVEB3jzUQ7OmRMQ27Tpyhrc3HwRgzrAIPFxqVItGapkac4muxWIhPj6eoUOHljf3AEJDQ8sbcQBxcXFYrVaGDx9OTk5O+Z/GjRvTqlUrNmzYUGG/np6ejBgxovxrZ2dnunTpwoEDB8q3ffzxx4SFhdG2bdsK+7ztttsAKu3zSsyfP5+1a9fy5ptvEh0dTXFxMaWlpVe9n7pg874cPt1xFDs7mHd3JM6ONeafnYhItbFaDcbHJVJqNejfrhEDIxqbHUlEaqAmTZqwdetWdu7cydq1a/nTn/5kdqTr7uuUk6xKOoGDvR3zY6NwsNeluSJie0osVsbHJWE1IKZDU3q38Tc7ktRyNaY9nJWVRUFBAaGhoZUe+/m29PR0DMOgVatWVe7Hyani2RABAQGV5nnUr1+fxMTECvtMTU2tdCnwz7NdrQ4dOpT/fcSIEXTs2JFRo0bxySefXPW+arPCEgsTl5fNmnowOpCOLeqbnEhExBwfbDvET4fP4OniyMyYcM2aEqlFJk+ezMcff0yLFi0AmDRpEnl5ecyaNYvCwkKio6NZtGgRhw8fJjY2lrCwMH788Uf69OnD4sWLAVi5cmWl9fb2lT/0dHb+311kCwsLsVqt1XOQ1eR8YQlTPyubNfWHnsG0a1rP5EQiIuZ4c9NBUjPP4ePuxJS72pkdR+qAGtPgu1JWqxU7OztWr16Ng4NDpcc9PSvOdqtqDZTd0fXn+4yMjOTll1+ucm3z5s1/Q+KyN2pDhgxh/vz5FBQU4Obm9pv2V5ssjE/nUG4+jeu58uztbcyOIyJiisyzBbywZg8A4wa2oYm37fwcEKnttm/fzvr160lOTiY7O7v8io833niDhIQEXF1defzxx4mLi6Nz586kpKSwbNkygoOD6dSpE+np6dSvX5+FCxdWWn/PPfdU+Zr79+9n6NChHDx4kA8++KCaj/jG+ttXe8g8W0gLX3ee7Fv1B/YiInXdodw8Fny9F4BJd4Th5+liciKpC2pMg8/f3x9XV1f27dtX6bGfbwsJCcEwDIKCgmjd+voMoAwJCWHXrl307dv3hp1RUVBQgGEYnD9/3mYafLuPn2Xpt2WXQs8aGoGXq2ZNiYjtMQyDKSt2c6GolI4tfBjRNdDsSCJyFbZs2cKwYcNwcnKiadOm9OzZExcXFxITE4mOjgbK3ucFBgbSuXNnwsLCyq80iYqKIiMjg7S0tCrXX0pISAhJSUkcOHCABx98kIEDB+Lq6nrjD/YG+/HQad777hAAc4dF4uZc9QfxIiJ1mWEYTFyeRFGplR6hDbinU4DZkaSOqDHD0BwcHOjXrx8rVqzg+PHj5dv37dvH6tWry7+OjY3FwcGBGTNmVDgLD8r+R8nKzmHr/lw+23mMcwUlV/Taw4cP59ixYyxdurTSYwUFBeTl5V3xcVR1Oe+ZM2f49NNPad68Of7+tnFdvcVqMCEuCYvV4I7IxvRv18jsSCI2xWI1ymvh1v25WKzGrz9JbojVySeITz2Jk4Md8++Owl6zpkSq1W+th798v3lxW0xMDDt37mTnzp3s2bOHcePGAeDi8r+zMOzt7bFYLJddfznBwcH4+vqSlJR0VZlrouJSKxPiEjEMuLtjALe08jM7kohN0XvDmuPTHcfYvC8XF0d75gyN1NgWuW5qzBl8ANOnT2ft2rX06NGDsWPHYrFYePXVV4mIiGDnzp1A2Seas2fPZsKECWRkZDB06FC8vLw4ePAgHy77BId2/bCLKrvT2InMc9gX5bEmOZOBEU0u+boPPPAAH330EWPGjGHDhg306NEDi8VCWloaH330EV999RWdO3e+omMYNGgQAQEBdO3aFX9/fw4fPszbb7/N8ePHWbZs2W/+HtUW72zJIPHoWbxcHZk+ONzsOCI2ZU1yJjO+SCHzbGH5tiberkwb3O6ytVCuv7MFJUz7fDcAY3uF0LqRl8mJRGzL9aiH3bt355lnnuHpp58mOzubTZs2MXLkSNatW8fRo0cJCAggNzeXgoKCS+4jOjqap59+utL6gIDKZ20cO3YMX19f3NzcOHHiBImJiQQF1f67zC7ZuJ+9Jy/QwMOZyXeGmR1HxKbovWHNkXuhiNkrUwB4ql9rWvp5mJxI6pIa1eDr1KkTq1ev5tlnn2XKlCk0b96cmTNnkpqaSlpaWvm68ePH07p1axYsWMCMGTMA8PVvwrkGYXgFdOLnF4KWWq2M/WAHi0d0vGTxsre3Z8WKFSxYsID33nuP5cuX4+7uTnBwME8++eRVXQr88MMP85///IcFCxZw5swZ6tevT3R0NP/617/o2bPnNX1fapujp/N5aW3ZrKmJd4ThX6/2X1IiUlusSc5k7Ac7+OVnsifOFv5qLZTrb/7qNLLPFxHc0IPH+lS+iZSI3DjXqx527dqVXr16ER4eTmBgIO3btyc0NJTXXnuNmJgYSkpKcHJyYunSpfj6+la5D39//yrXV9XgS0lJ4ZlnnsHe3h57e3teeukl/Pxq99luB7Iv8I/1ZSN3pg5uR30P5195hohcL3pvWLPM+jKFM/klhDWpx6M9a/+HN1Kz2BlVXXdQwwwdOpTdu3eTnp5e5eMWq8Etz6+v8InEz9kBjb1d2fTcbTjo0qgbyjAMHnrnexL2ZNOlpS//GR2ty9FEqolqYc2y7UAuv1vyHQDLRkfTNbiByYlEbMf1rod5eXl4eHiQk5NDt27d2LFjB15eOiP3SlitBr9f+h3bDp7i1tYNefehm3U5mkg10XvDmuWbvdmMfGs79naw/LEetG/uY3YkqWNqzAy+i355eUN6ejqrVq2id+/el3zO9oOnLlm0AAwg82wh2w+euk4p5VI+33WchD3ZODvYMzc2Us09kWqkWlhzFJZYmLC8bGbW77s0V3NPpJpd73r48MMP06FDB3r16sWsWbPU3LsKH/1whG0HT+Hm5MCcoRFq7olUI703rDnyi0uZ9N/3hqO6B6m5JzdEjbpEF8qGCY8aNYrg4GAOHTrE4sWLcXZ2vuwg4qzzly5a17KuKmfPnr3sbBWAxo0bX/P+64LTecXM/KJsnsCfbgsl1N/T5EQitqU6aqFcmUUb9nEgO4+GXi6MH6RZUyLV7XrXwxsxRzk3N5e+fftW2Obn50d8fPx1fy2zZJ0vZO6qVAD+MqA1zX3dTU4kYlv03rDmeHntXo6eLqCZjxt/GXDlI8BErkaNa/ANHDiQf//735w4cQIXFxe6devG3LlzadWq1SWf4+91ZTPernRdVZ588knefffdy66pBVc731BzV6WSm1dMK39PxvQKMTuOiM2pjloov27vyfMs/mY/ANMHh+Pt5vQrzxCR66021MMGDRqU30SurprxRQrnCkuJbObNqO4tzY4jYnNqQy20BUlHz/LW5oMAzB4agYdLjWvDSB1R4/5lvf3221f9nC5BvjTxduXE2cJKw0Phf7MFugRVPfj4SowbN44RI0Zc8/Prui37cvj4x6PY2cH8uyNxdqxxV3+L1HnVUQvl8qxWgwlxSZRYDPqF+XNHpG2f2S1iFtVD861LPcnKxEwc7O2YFxuJo4PeG4pUN9VC85VarIyPS8RqwOD2TenT1t/sSFKH1bgG37VwsLdj2uB2jP1gB3ZQoXhdnPIxbXC73zQ4tF27drRr1+63xKyzfj5rakTXQDoF6geEiBmqoxbK5X247RA/HjqNh7MDM2M0a0rELKqH5rpQVMrkFckAPHpLEBHNvE1OJGKbVAvN9+amg+w+fg5vNyem3qV+gtxYdeajtIERTVg8oiONvSueXtzY21W3/r7B/r4unUO5+TSu58q4gW3MjiNi01QLzZN5toDn1+wBYNzAtjT1cTM5kYhtUz00z9++2kPm2UJa+LrzVD/NmhIxk2qheQ7n5rMgfi8Ak+4Mo6GXi8mJpK6zM+rY4DiL1WD7wVNknS/E36vsdGN9InHjpGae465XNmGxGix5oBMDwnU5mkhNoFpY/Ua/9wNrU05yUwsfPhnTXd9vkRpC9bB6/XT4NLGLt2AY8P4jXejZqqHZkUQE1cLqZhgGD761nW/Tc+ge0oAPH+2qKzvkhqsTl+j+nIO9Hd1CGpgdwyZYrAbjP03EYjUYFNFYzT2RGkS1sHqtSc5kbcpJHP87a0pvmEVqDtXD6lNisTIhLgnDgNiOzdTcE6lBVAur1/KfjvFteg4ujvbMHRap5p5Uizpzia5Uv3e3ZLDr6Fm8XB2ZPiTc7DgiIqY4W1DC1M92AzCmVwhtG9czOZGIiDmWbDxA2onz+Ho4M/lOzZoSEduUe6GIWV+mAPBE31a09PMwOZHYCjX45JocPZ3P39aWzZoaP6gtjerp1uoiYpueX5NG1vkigv08+NNtoWbHERExxYHsC/x9XToAU+4Kw9fD2eREIiLmmL0yldP5JbRt7MXoW4PNjiM2RA0+uWqGYTBlRTL5xRa6tPTl9ze3MDuSiIgpth88xb+2HQZgbmwkrk4OJicSEal+hmEwcXkSxaVWbm3dkKEdmpkdSUTEFN/szWb5T8ews4P5d0fh5KCWi1Qf/WuTq/ZlYiYb9mTj7GDP3NhI7DVrSkRsUFGphQlxiQD8rnNzooM110ZEbNPHPxzluwOncHNyYM7QCM2aEhGblF9cyqTlSQCM6t6SDs19zA0kNkcNPrkqZ/KLmfFF2aypx/uEEurvaXIiERFzLNqwn/3Zefh5ujDxjjCz44iImCL7fBFzVqUC8Ez/1jT3dTc5kYiIORbGp3P0dAHNfNx4dkAbs+OIDVKDT67K3FWp5FwoJtTfkzG9NU9ARGxT+snzLErYB8D0Ie3wdncyOVFFmzdvpnPnzkRGRtKtWzeSkpLMjiQiddSML3ZztqCEiGb1eKhHS7PjiIiYIvnYWf757QEAZg0Nx8PF0eREYov0r06u2Jb9OXz0w1EA5sdG4uKoWVMiYnusVoPxcUmUWAz6tvXnzsgmZkeqxMfHhxUrVhAQEEB8fDxPPPEEGzZsMDuWiNQx69NO8mViJg72dsyPjcJRs6ZExAaVWqw892kiVgPuimrCbW0bmR1JbJR+CssVKSyxMDGu7AyQEdEt6NzS1+REIiLm+HD7YX48dBoPZwdmXcWsqcmTJ9OmTRv69+9P//79SUhIYOXKlURHR9OhQwfGjBmD1WolIyODjh07cv/999O2bVvGjh1bvo+q1lclPDycgIAAADp16sSxY8d++4GLiPzMhaJSJi9PBuCRW4KIaOZtciIREXO8tfkgu4+fw9vNiWmDw82OIzZMDT65Iq+sTycjN59G9VwYN7Ct2XFERExx4mwhL6xOA+DZ29vQ1Mftip63fft21q9fT3JyMu+++y7btm0jJyeHhQsXkpCQwM6dO3FwcCAuLg6AlJQUpk+fzu7du9m6dSvp6emXXX857733Hn379r32gxYRqcJLa/dw/GwhzX3deKpfK7PjiIiY4sipfF7+ei8Ak+4Io6GXi8mJxJbpEl35VamZ53jjm7J5AjOGRFDPtWbNmhIRqS7TPk/mfFEpHZr78GC3llf8vC1btjBs2DCcnJxo2rQpPXv2xMXFhcTERKKjowEoKCggMDCQzp07ExYWRqtWZb8wR0VFkZGRQVpaWpXrL+fHH3/k9ddfZ9OmTdd2wCIiVdh55AzvbMkAYM7QSNyd9SuFiNgewzCYuDyJwhIr0cG+3Ns5wOxIYuP001guy/LfWVOlVoPbwxsxMKKx2ZFEREyxJvkEX+0+iaO9HfPvjsTB/souzYWyN4BVbYuJiWHJkiUVtmdkZODi8r9Pf+3t7bFYLJdcfylHjx7l//7v//joo49o0KDBFWcVEbmcEouV8Z8mYhgw7KZm3Nq6odmRRERMsWLnMb5Nz8HZ0Z55sVFXPLZF5EbRJbpyWe9tzWDXkTN4uTgyMybC7DgiIqY4V1jC1M/KZk39sVcwbRvXu6rnd+/enRUrVlBaWkpmZiabNm2isLCQdevWcfRo2c2LcnNzy/9elejo6Cten5eXR0xMDC+99BLt27e/qqwiIpezZOMB0k6cp767E5PvDDM7joiIKU7lFTPry1QAnuzbiiA/D5MTiajBJ5dx7EwBL361B4DnBrWlUT1XkxOJiJjjhTVpZJ0vIsjPgz/fdvWzprp27UqvXr0IDw9n5MiRtG/fntDQUF577TViYmKIiopiwIABZGVlXXIf/v7+V7z+7bffJi0tjcmTJ9OhQwe6d+9+1ZlFRH7pYE4ef1+XDsCUu9rRwFOzpkTENs1emcKpvGLaNvZi9K3BZscRAcDOqOq6IbF5hmHw6Ls/sC4ti86B9fnoj92wv4rL0URE6oofMk5xz+tbAfjXH7rSPcTvmvaTl5eHh4cHOTk5dOvWjR07duDl5XU9o4qI3DCGYXD/P7exZX8uPVv58d7DXXQ5mojYpG/Ts3ngze3Y2UHc2O7c1KK+2ZFEAM3gk0tYmZTJurQsnBzsmBcbqeaeiNikolIL4+OSABjeOeCam3sADz/8MHv27KGkpIRZs2apuScitcrHPx5ly/5cXJ3smTM0Us09EbFJBcUWJi0vG9sysltLNfekRlGDTyo5m1/C9M9TAHisdyitGumXUBGxTYsT9rMv6wJ+ns5MvOO3zZpatmzZdUr1P7m5ufTt27fCNj8/P+Lj46/7a4mI7co+X8SclWWzpp7u15oWDdxNTiQiYo6F8Xs5fCqfpt6uPHt7G7PjiFSgBp9UMndVKjkXigj19+SxPiFmxxERMcW+rPMs2rAfgGmDw/FxdzY5UWUNGjRg586dZscQkTpu5pcpnC0oIbxpPR65JcjsOCIipkg+dpZ/bjoIwKyhEXi6qJ0iNYtusiEVbN2fy7IfjgAwLzYSF0cHkxOJiFQ/q9VgQlwSxRYrfdo05K6oJmZHEhExxYa0LL7YdRx7O5gfG4Wjg359EBHbU2qxMiEuCYvV4M6oJvQNa2R2JJFK9BNayhWWWJi0vGzW1P1dW3BzS1+TE4mImOPf3x/m+4zTuDs7MHuYZk2JiG3KKypl8oqyWVOP3BJEZIC3yYlERMzxzpYMko6dpZ6rI9MGtzM7jkiV1OCTcq+u38eBnDz8vVx4blBbs+OIiJji5LlC5q9KA+DZAW1o5uNmciIREXO8tHYvx84UEFDfjaf7tzY7joiIKY6cyueltXsBmHhHGP5eriYnEqmaGnwCQNqJc7z+TdmsqZkx4dRzdTI5kYiIOaZ9tpvzRaW0b+7DyO4tzY4jImKKXUfO8M6WsllTc4ZF4u6sWVMiYnsMw2DSimQKSix0DfLldzc3NzuSyCWpwSdYrAbjP02i1GowoF0jBkZo1pSI2Kavdp9gze4TONrbMT82Egd7XZorIranxGLluU8TsRowtENTerVuaHYkERFTfLbzOBv3ZuPsaM+8WI1tkZpNDT7hg+8OsfPIGTxdHJkZE2F2HBERU5wvLGHqZ2Wzpv5wazBhTeqZnEhExBz//PYgaSfO4+PuxJS7NGtKRGzT6bxiZn6ZAsATt4US3NDT5EQil6cGn407fqaAF9aUzZp6bmAbGntrnoCI2KYX1uzh5LkiWjZw58m+rcyOIyJiioycPBbGl82amnxnOxp4upicSETEHLNXpnIqr5g2jbwYfWuI2XFEfpUafDbMMAymfpZMXrGFToH1ub9roNmRRERM8eOhU3yw7RAAc4dF4urkYHIiEZHqVzZrKomiUiu3hPpxd8dmZkcSETHFpvQcPt1xFDs7mHd3JM6Oap1Izad/pTZsVdIJ4lOzcHIomzVlr1lTImKDikotjP80CcOAezsF0D3Uz+xIIiKm+OTHo2zel4urkz1zhkVo1pSI2KSCYgsTlycB8GB0IB1b1Dc5kciVUYPPRp3NL2Ha57sBGNs7lFaNvExOJCJijtcTDpCedYEGHs5MvCPM7DgiIqbIuVDEnFWpADzVrzWBDTxMTiQiYo6F6/Zy+FQ+Tbxd+evAtmbHEbliavDZqPlrUsm5UERwQw8e76N5AiJim/ZlXeC1DfsAmDq4HfU9nE1OJCJijllfpnAmv4R2Terx6C1BZscRETHF7uNn+ee3BwGYFROBp4ujyYlErpwafDbouwO5/Hv7EQDmx0bh4qhZUyJie6xWg4lxSRRbrPRu05Ah7ZuaHUlExBQb9mTx2c7j2NvB/LsjcXTQrwgiYnssVoMJcUlYrAZ3RDamX7tGZkcSuSr66W1jCkssTIwrmyfw+y4t6BLka3IiEZEb49ChQ3Ts2JEOHTrQvn17Pv/88wqP/+f7I2zPOIW7swOzh2rWlIjYpryiUiYvTwbgoR5BRAX4mBtIRMQkb28+SOLRs3i5OjJ9cLjZcUSums43tTGvbdjHgZw8/L1cGD9I8wREpO5q0qQJW7duxcXFhZMnT3LzzTczZMgQALLOFTJvddmsqb8MaENAfXczo4qImOblr/dy7EwBzXzceKZ/a7PjiIiY4sipfF5auxeAiXeE4V/P1eREIldPZ/DZkD0nzrM4YT8AM4aE4+3mZHIiEbE1kydPpk2bNvTv35/+/fuTkJDAypUriY6OpkOHDowZMwar1UpGRgYdO3bk/vvvp23btowdO7Z8H1Wtr4qzszMuLi4AFBYWVlg3/YvdnC8sJSrAm1HdW97QYxYRqakSj57h7c1ls6bmDIvAQ7OmRMQGGYbB5BXJFJRY6BLky+86Nzc7ksg1UYPPRlitBhPiEim1GvRv14iBEY3NjiQiNmb79u2sX7+e5ORk3n33XbZt20ZOTg4LFy4kISGBnTt34uDgQFxcHAApKSlMnz6d3bt3s3XrVtLT0y+7vir79+8nMjKS8PBwXn31VQC+TjnJqqQTONjbMT82Cgd7XZorIranxGJl/KdJWA2I6dCU3m38zY4kImKKz3cd55u92Tg72DMvNhJ7vTeUWkof09mID7YdYsfhM3i6ODIzJlyzpkSk2m3ZsoVhw4bh5ORE06ZN6dmzJy4uLiQmJhIdHQ1AQUEBgYGBdO7cmbCwMFq1agVAVFQUGRkZpKWlVbn+UkJCQkhKSuLAgQM8+OCD9OjdlykrymZN/aFnMO2a1rvBRy0iUjO9uekgKZnn8HF3Yspd7cyOIyJiitN5xcz8IgWAP90WSkhDT5MTiVw7NfhsQObZAl5YsweAcQPb0MTbzeREImKLDMOocltMTAxLliypsD0jI6P88loAe3t7LBbLJdf/muDgYHx9fZnwz5WcOOdFYAN3nurX6toORESkljuUm8eCr8tmTU26Iww/T5dfeYaISN00Z1UquXnFtG7kyZheIWbHEflNdIluHWcYBlNW7OZCUSkdW/gwouulz3QREbmRunfvzooVKygtLSUzM5NNmzZRWFjIunXrOHr0KAC5ubnlf69KdHT0Fa8/duwYBQUFAJw4cYIfduxk7ZGyOXxzhkbi6uRwPQ9PRKRWMAyDicuTKCq10iO0Afd0CjA7koiIKTbvy+GTH49iZwfzYqNwdlR7RGo3/Quu49YknyA+9SRODnbMvztK8wRExDRdu3alV69ehIeHM3LkSNq3b09oaCivvfYaMTExREVFMWDAALKysi65D39//yten5KSQpcuXWjfvj0DBw6iYb9HsXfz5u6OAdzSyu9GHaaISI0Wt+MYm/fl4uJoz5yhkRrbIiI2qbDEwsTlSQA8EB1Ip8D6JicSKRs/1LdvXzw9PRk/fvxVP9/OqOqaKakTzhaU0O/lb8g+X8QTt4XyzIA2ZkcSERuXl5eHh4cHOTk5dOvWjR07duDl5XXDX/eVdem89PVeGng4E/9ML+p7ON/w1xQRqWlyLxTR7+VvOJ1fwnMD2zK2ty5HExHb9PyaNBYn7KdxPVe+fuZWvFydzI4kQklJCVu2bCElJYVDhw4xf/78q3q+zuCrw+avTiP7fBHBDT14rE+o2XFERHj44Yfp0KEDvXr1YtasWdXS3NuffYFX1u8DYOrgdmruiYjNmvVlCqfzSwhrUo9HewaZHUdExBQpx8+xZOMBAGbGhKu5J5c1efJk2rRpQ//+/enfvz8JCQmsXLmS6OhoOnTowJgxY7BarWRkZNCxY0fuv/9+2rZty9ixY8v3UdX6qjg5OdGrVy/c3K7tvgm6yUYdte1ALv/efhiAecM0a0pEaoZly5Zd933m5ubSt2/fCtv8/PyIj4/HajWYEJdEscVKr9YNGdK+6XV/fRGR2uCbvdms2HkcezuYHxuJk4M+5xcR22OxGkyIS8RiNRgU0ZgB4Y3NjiQ12Pbt21m/fj3JyclkZ2fTtm1bcnJyeOONN0hISMDV1ZXHH3+cuLg4OnfuTEpKCsuWLSM4OJhOnTqRnp5O/fr1WbhwYaX199xzz3XPqwZfHVRYYmHCf+cJ/L5Lc7oGNzA5kYjIjdOgQQN27txZ5WPLfjjC9oOncHNyYPbQCM2aEhGblF9cyqT/vjcc1T2I9s19zA0kImKSd7ZksOvoWbxcHZkxJNzsOFLDbdmyhWHDhuHk5ETTpk3p2bMnLi4uJCYmEh0dDZTNzQsMDKRz586EhYXRqlUrAKKiosjIyCAtLa3K9TeCGnx10KKE/RzIzqOhlwvjB4WZHUdExBRZ5wqZuyoVgL8MaE1zX3eTE4mImGPB13s5erqAZj5u/GVAa7PjiIiY4ujpfF5auweACYPC8K/nanIiqemqumWFYRjExMSwZMmSCtszMjJwcXEp/9re3h6LxXLJ9TeCzs2vY/aePM/ihLJZU9MHh+PtpnkCImKbZnyRwvnCUiKbeTOqe0uz44iImCLp6Fne3HQQgNlDI/Bw0ef7ImJ7DMNgyopk8ostdGnpy303Nzc7ktQC3bt3Z8WKFZSWlpKZmcmmTZsoLCxk3bp1HD16FCgbF3Tx71WJjo6+qvW/hX7C1yEXZ02VWAz6hflzR6TmCYiIbYpPOcnKpEwc7O2Yf3ckjpo1JSI2qNRiZXxcIlYDBrdvSp+2/mZHEhExxReJmWzYk42zgz1zYyOxt9fYFvl1Xbt2pVevXoSHhxMYGEj79u0JDQ3ltddeIyYmhpKSEpycnFi6dCm+vr5V7sPf37/K9QEBAVWuj4qK4ujRo5SWlvLBBx/w008/0bBhwyvKa2dUdc6h1Ervb81gyme78XB24OtnetHU59ruvCIiUpudLyxhwIKNZJ4t5I+9gpmgUQUiYqPe+GY/81an4e3mRPwzvWjo5fLrTxIRqWPO5BfT7+VvyLlQzNP9WvNkv1ZmR5JaJC8vDw8PD3JycujWrRs7duzAy8vL7FhV0hl8dcSJs4U8v6ZsnsC4gW3V3BMRm/XS2r1kni2kha87T/XVrCkRsU2Hc/NZEL8XgEl3hqm5JyI2a+6qVHIuFNPK35OxvUPMjiO1zMMPP8yePXsoKSlh1qxZNba5B2rw1RlTP0vmQlEpN7XwYUT0jbkji4hITbfj8Gne3ZoBwJxhEbg5O5gbSETEBIZhMGlFEoUlVroFN+DeTlVfBiQiUtdt2Z/DRz+UzTubf3ckzo4a2yJXZ9myZdd9n7m5ufTt27fCNj8/P+Lj43/TftXgqwPWJGeyNuUkjvZ2zIuNxEHzBETEBhWXWpnwaRKGAbEdm9Gz1ZXNqhARqWuW/3SMb9NzcHYsmzVlZ6f3hiJiewpLLEyMSwJgRHQLOgVWPSNNpLo1aNCAnTt3Xvf9qn1dy50tKGHqZ7sBGNMrhLaN65mcSETEHEs27mfPyfP4ejgz+c52ZscRETFF7oUiZn2ZAsCTfVsR5OdhciIREXP8Y106Gbn5NKrnwriBbc2OI3LDqcFXyz2/Jo2s80UE+3nwp9tCzY4jImKKA9kX+Mf6fQBMuSsMXw9nkxOJiJhj9spUTueX0LaxF6NvDTY7joiIKVIzz7Fk4wEAZsZEUM/VyeREIjeeGny12PcZp/jXtsMAzI2NxNVJs6ZExPYYhsGEuCSKS63c2rohQzs0MzuSiIgpNu7NZvlPx7Czg/l3R+HkoLf6ImJ7LFaD8XFJlFoNBoY35vbwxmZHEqkW+qlfSxWVWhj/aSIA993cnOjgBiYnEhExx0c/HGHbwVO4OTkwZ2iEZk2JiE3KLy5l0oqyWVOjurekQ3MfcwOJiJjkva0Z7DpyBi8XR2bEhJsdR6TaqMFXSy3asJ/92Xn4ebowYVCY2XFEREyRdb6QOStTAXimf2ua+7qbnEhExBwL49M5cqqApt6u/GVAG7PjiIiY4tiZAl78ag8Azw1qS6N6riYnEqk+avDVQuknz7MooWzW1PQh7fB21zwBEbFNM75I4VxhKRHN6vFQj5ZmxxERMUXysbP889uyWVOzh0Xg6eJociIRkepnGAZTViSTX2zh5pb1+b8uLcyOJFKt1OCrZaz/nSdQYjHo29afOyObmB1JRMQU61JPsjIxEwd7O+bHRuGoWVMiYoNKLVae+zQRqwF3RTXhtraNzI4kImKKLxMzWZ+WhbODPfNiI7G319gWsS36baiW+df2w/x46DQezg7M0qwpEbFRF4pKmbIiGYBHbgkiopm3yYlERMzx9uYMdh8/Rz1XR6YN1qwpEbFNZ/NLmPHFbgAe6xNCqL+XyYlEqp8afLXIibOFPL86DYBnb29DUx83kxOJiJjjb1/t4fjZQpr7uvF0v9ZmxxERMcWRU/m8/PVeACbdGUZDLxeTE4mImGPuqlRyLhQT6u/J2N4hZscRMYUafLXItM+TOV9USofmPjzYraXZcURETPHT4dO8uzUDgDlDI3FzdjA3kIiICQzDYOLyJApKLEQH+zK8c3OzI4mImGLr/lyW/XAEgHmxkbg46r2h2CY1+GqJNckn+Gr3SRzt7Zh/dyQOmicgIjaoxGJlQlwShgHDbmrGra0bmh1JRMQUK3Ye49v0HJwd7ZkXG6WxLSJikwpLLExcngTA/V1bcHNLX5MTiZhHDb5a4FxhCdM+L5s19cdewbRtXM/kRCIi5liy8QBpJ85T392JyXeGmR1HRMQUp/KKmfVlKgBP9m1FkJ+HyYlERMzx6vp9HMzJo1E9F54b1NbsOCKmUoOvFnhhTRonzxUR5OfBn29rZXYcERFTHMzJ4+/r0gGYclc7Gnhq1pSI2KbZK1M4lVdMm0ZejL412Ow4IiKmSDtxjte/2Q/AjCER1HN1MjmRiLnU4Kvhfsg4xQffHQZgzrAIXJ00T0DkUjIzMxk/fjx9+vTBy8sLOzs7EhISzI4l14FhGEyMS6K41ErPVn4Mu6mZ2ZFEREzxbXo2cTuOYWcH8++OxMlBb+dFxPZYrAbjP02i1GowoF0jBkY0NjuSiOn0jqAGKyq1MD6ubJ7A8M4BdA/xMzmRSM22Z88enn/+eY4dO0ZkZKTZceQ6+viHo2w9kIurkz1zhkZq1pSI2KSCYguTlpeNbRnZrSU3tahvciIREXO8vzWDnUfO4OXiyMyYCLPjiNQIjmYHkEtbnLCffVkX8PN0ZuIdmjUl11dxqZX3t2Zw6FQ+gb7uPNCtJc6Otbvn36lTJ3Jzc/H19eWTTz7h3nvvNTuSXAfZ54uYs6ps1tTT/VrTooG7yYmkrqmL9VDqpoXxezl8Kp+m3q48e3sbs+NIHaNaKLXF8TMFvPjVHgDGDWpLY29XkxNJXVNb66EafDXUvqzzLNpQNk9g2uBwfNydTU4kdcm8VSks/fYgVuN/2+asSqWvdzY/ffIqycnJNGvWjHHjxpGZmcmMGTMwjP8t/uCDD1iwYAEpKSm4ubkxYMAAXnzxRZo3b16+pnfv3uTk5PDRRx/x+OOPs23bNurXr8+TTz7JuHHjKuQpKipi7ty5fPjhhxw5cgR/f39+//vfM2vWLFxcrnzOmpeX17V/U6TGmvllCmcLSghvWo9HbgkyO47UMZeqh3/oGcSEO9qZF0zkF5KPneWfmw4CMGtoBJ4uehsv149qodQWhmEwZUUyecUWOgfW5/4uLcyOJHVMba6HemdQA1mtBhPikii2WLmtrT93RTUxO5LUIfNWpfDGxoOVthee2M8/X3yW+n7+zJgxA4vFwsyZM2nYsGGFdXPmzGHKlCkMHz6cRx99lOzsbF555RVuvfVWfvrpJ3x8fMrXnj59moEDBxIbG8vw4cP55JNPeO6554iMjGTQoEEAWK1WhgwZwqZNmxg9ejRhYWEkJSWxYMEC9u7dy4oVK27kt0NquA1pWXyx6zj2djA/NgpHzZqS6+hS9dBqUL69pr+RE9tQarEyIS4Ji9Xgzqgm9A1rZHYkqUNUC6U2WZV0gnVpWTg52DEvNhJ7e41tkeunttdDNfhqoH9/f5jvM07j7uzArKERmjUl101xqZWl31YuWABnNn2Inb09nvfM5em/jMDZ0Z7hw4cTFva/y8MPHTrEtGnTmD17NhMnTizfHhsby0033cSiRYsqbD9+/DjvvfceDzzwAACPPPIIgYGBvPnmm+UNvn/961/Ex8fzzTffcMstt5Q/NyIigjFjxrBlyxa6d+9+Xb8PUjvkFZUyeUXZrKlHbgkiMsDb5ERyPRQUFHDXXXexbds2/vSnPzF//nxTclyuHl609NuD/GVA21pxSYbUbe9sySDp2FnquToybXDN/cVCah/VQqlNzuaXMO3z3QCM7R1Kq0a6ekeun7pQD2tmKht28lwh81elAfDsgDY083EzOZHUJe9vzahwqvFFhtVCYcYu3FpFY+/ZgPe3ZgAQGhpa3ogDiIuLw2q1Mnz4cHJycsr/NG7cmFatWrFhw4YK+/X09GTEiBHlXzs7O9OlSxcOHDhQvu3jjz8mLCyMtm3bVtjnbbfdBlBpn2I7/rZ2D8fOFBBQ342n+7c2O45cJ46OjkydOpUXX3zR1ByXqoc/ZzUor4ciZjlyKp+X1u4FYOIdYfh7adaUXD+qhVKbzFudSs6FIkIaevB4nxCz40gdUxfqoc7gq2Gmfbab80WltG/uw8juLc2OI3XMoVP5VW635J/FKC3CyadppXWhoaHlf09PT8cwDFq1alXlfpycnCp8HRAQUOkM1Pr165OYmFhhn6mpqZUuBb4oKyvrMkckddXOI2d4Z0sGAHOGReLurB9XN9LkyZP5+OOPadGibI7NpEmTyMvLY9asWRQWFhIdHc2iRYs4fPgwsbGxhIWF8eOPP9KnTx8WL14MwMqVKyutt7ev/Dmik5MTvXr14uDBy39CeqNdqh5e6zqRG8EwDCatSKagxELXIF9+d3PzX3+SyFVQLZTa4rsDufzn+yMAzL87ChdHB5MTSV1TF+qhfmOqQb7afYI1u0/gaG/H/NhIHDRPQK6zQN8ru/vopdZZrVbs7OxYvXo1Dg6Vf6h6enpW+LqqNUCFG3ZYrVYiIyN5+eWXq1z78xt3iG0osVgZ/2kihgFDOzSlV+uqm79yfWzfvp3169eTnJxMdnZ2+dm0b7zxBgkJCbi6uvL4448TFxdH586dSUlJYdmyZQQHB9OpUyfS09OpX78+CxcurLT+nnvuMfvwLum31kOR6vD5ruNs3JuNs6M982IjNbZFrjvVQqkNCkssTIxLAuD/urbg5pa+JieSuqgu1EM1+GqI84UlTPusbJ7AH24NJqxJPZMTSV30QLeWzFmVWunUYwd3b+wcnSk5U3Yzgwe6tSx/bN++feV/DwkJwTAMgoKCaN36+lwyGRISwq5du+jbt69+cREAln57gLQT56nv7sSUuzRr6kbbsmULw4YNw8nJiaZNm9KzZ09cXFxITEwkOjoaKJubFxgYSOfOnQkLCys/izcqKoqMjAzS0tKqXF+TXaoe/twv66FIdTqdV8zML1IA+HOfUIIbev7KM0Sunmqh1AavbdjHgZw8/L1ceG5gW7PjSB1VF+qhZvDVEC+s2cOJc4W0bODOk32rvvxR5LdydrTnDz2DKm23s3fANbA9Benf8bt2HuVDQ/ft28fq1avL18XGxuLg4MCMGTMqnIUHZWfl5ebmXnWm4cOHc+zYMZYuXVrpsYKCAvLy8q56n1J7ZeTk8ff4dAAm39mOBp4uJieq+375//LFbTExMezcuZOdO3eyZ88exo0bB4CLy//+m9jb22OxWC67vqa6VD38uT/0DKqxQ5Sl7pu9MpXcvGLaNPLij700a0puDNVCqen2nDjP4oT9AMwYEo63m9OvPEPk2tSFeqgz+GqAHw+d4oNthwCYOywSVyfNE5Ab5+JtvZd+e7DCpxP1e97PyQ928p+pD1E/cywWi4VXX32ViIgIdu7cCZSdbTd79mwmTJhARkYGQ4cOxcvLi4MHD7J8+XJGjx7Ns88+e1V5HnjgAT766CPGjBnDhg0b6NGjBxaLhbS0ND766CO++uorOnfufMX7mz17NgC7d5edEfv++++zadMmoGzOmNRchmEwcXkSRaVWbgn1I7ZjM7Mj2YTu3bvzzDPP8PTTT5Odnc2mTZsYOXIk69at4+jRowQEBJCbm0tBQcEl9xEdHc3TTz9daX1AQEA1HsnVu1Q9tLcrewN38XGR6rYpPYdPdxzFzg7m3R1Zo3+ZkNpPtVBqKovVYHxcIqVWg/7tGjEworHZkaSOq+310M6o6qN7qTbFpVbu/Me3pGdd4N5OAbx4b3uzI4mNKC618v7WDA6dyifQ150HurVk08YEnn32WXbv3k3z5s2ZMGECqampvPbaaxV+uY+Li2PBggX89NNPQNmcvL59+/LEE0+UX7rbu3dvcnJySE5OrvC6o0aNIiEhgYyMjPJtJSUlLFiwgPfee499+/bh7u5OcHAwQ4YM4amnnqJevSu/ZP1yl/mq3NVsH/9whL9+koirkz1fPXUrgQ08zI5kMyZOnMinn35KYGAghYWFLFy4kKysLCZNmkRJSQlOTk4sXboUX19f7rvvPr777jug7P/n++67j4EDB7JmzZpK6zt27Fjl60VFRXH06FFKS0upV68eP/300yVvtFMdqqqHaqiIWQqKLdy+cCOHT+UzslsgM2IizI4kNkK1UGqa97ZmMPWz3Xi6OBL/TC8ae+su4lI9ams9VIPPZP9Yl87LX++lgYcz6/7SCx93Z7MjiVQwdOhQdu/eTXp6utlRpA7LuVBEv5e/4Ux+CeMHtWWMLkerVnl5eXh4eJCTk0O3bt3YsWMHXl5eZscSsUnzV6fx+jf7aeLtytfP9MLTRRfciIjtyTxbQP+XN3KhqJRZMeE1eu6ZSE2hdwwm2pd1gVfXl93AYOrgdmruiekKCgpwc3Mr/zo9PZ1Vq1YxcuRIE1OJLZj5RQpn8kto16Qej95y+dkXcv09/PDD7Nmzh5KSEmbNmqXmnohJdh8/y9JvDwAwMyZCzT0RsUmGYTBlxW4uFJXSsYUP93et2TfuEqkp9K7BJFarwcS4JIotVnq3aciQ9k3NjiRCcHAwo0aNIjg4mEOHDrF48WKcnZ1NH5Z/9uzZy87/AmjcWDM5aqsNe7L4fFfZHZzn3x2Jo0PNP/29rlm2bNl132dubi59+/atsM3Pz4/4+Pjr/loidYHFajAhLgmL1eCOyMb0b9fI7EgiIqZYnXyC+NSTODnYMf/uKOztLz2CR0T+Rw0+k/zn+yNszziFu7MDs4dGXHZumEh1GThwIP/+9785ceIELi4udOvWjblz59Kqlbl3dn7yySd59913L7tG0wZqp7yiUiYvL5vT+FCPIKICfMwNJNdNgwYNym/QIyK/7u3NB0k8ehYvV0emDw43O46IiCnO5pcw7fOym+WN7RVC60a6qkDkSqnBZ4Ksc4XMW50KwF8GtCGgvrvJiUTKvP3222ZHqNK4ceMYMWKE2THkBnj5670cO1NAMx83nunf2uw4IiKmOHIqn5fW7gVg4h1h+NfTIHkRsU3z16SRfb6I4IYePNYn1Ow4IrWKGnwmmP7Fbs4XlhIV4M2o7i3NjiNS47Vr14527Wr2Lcnl6u06coa3Nx8EYM6wCDw0a0pEbJBhGEz5LJmCEgtdgnz5XefmZkcSETHFtgO5/Hv7YQDmDYvE1cnB5EQitYsGHVWzr1NOsirpBA72dsyPjcJB8wRExAaVWKyMj0vCakBMh6b0buNvdiQREVN8vus4CXuycXawZ+6wSM2aEhGbVFhiYcLyJAB+36U5XYMbmJxIpPZRg68anS8sYcqKsllTf+gZTLum9UxOJCJijn9+e5DUzHP4uDsx5S6dnSkitul0XjEzv0gB4E+3hRLq72lyIhERcyzasI8D2Xk09HJh/KAws+OI1Epq8FWjF7/aw4lzhQQ2cOepfubetEBExCwZOXksjC+bNTXpjjD8PF1MTiQiYo45q1LJzSumdSNPxvQKMTuOiIgp9p48z+Jv9gMwY0g43m5OJicSqZ3U4KsmPx46zfvfHQJgruYJiIiNMgyDSSuSKCq10iO0Afd0CjA7koiIKbbsy+GTH49iZwfzYqNwdtTbchGxPVarwfhPEymxGPQLa8SgiMZmRxKptfROohoUl1qZEJeIYcDdHQPoEepndiQREVN8uuMYm/fl4uJoz5yhkdjZadaUiNien8+aGtE1kE6B9U1OJCJijg+3HWLH4TN4ODswMyZc7w1FfgM1+KrBG9/sZ+/JCzTwcGbynZonICK2KedCEbNXls2aeqpfa1r6eZicSETEHH9fl86h3Hwa13Nl3MA2ZscRETFF5tkCnl+zB4BxA9vS1MfN5EQitZsafDfY/uwLvLJ+HwBTB7ejvoezyYlERMwx68sUzuSXENakHo/2DDI7joiIKVKOn2PJxgMAzIwJx8tVs6ZExPYYhsHUz3ZzoaiUm1r4MCI60OxIIrWeGnw3kNVqMCEuiWKLlV6tGzKkfVOzI4mImCJhTxaf7TyOvR3Mj43EyUE/fkTE9lisBhPiErFYDQZFNGZAuGZNiYht+mr3Cb5OOYmTgx3zY6NwsNeluSK/lX7DuoE++uEI2w+ews3JgdlDIzRPQERsUn5xKZNXJAMwqnsQ7Zv7mBtIRMQk727JYNfRs3i5OjJjSLjZcURETHG2oISpn+0GYEyvENo09jI5kUjdoAbfDZJ1vpC5q1IB+MuA1jT3dTc5kYiIOV5eu5ejpwto5uPGXwa0NjuOiIgpjp7O529ry2ZNjR/UFv96riYnEhExx/Nr0sg6X0SwnweP9wk1O45InaEG3w0y4/MUzhWWEtnMm1HdW5odR0TEFIlHz/DW5oMAzB4agYeLo8mJRESqn2EYTFmRTH6xhS4tffn9zS3MjiQiYortB0/xr22HAZgbG4mrk4PJiUTqDjX4boD4lJOsTMrEwd6O+XdH4qhZUyJig0osVsZ/moTVgMHtm9Knrb/ZkURETPFFYiYb9mTj7GDP3NhI7DVrSkRsUFGphQlxiQDcd3NzooMbmJxIpG5R5+k6u1BUypTPymZNPdoziPCm3iYnEhExx1ubDpKSeQ5vNyem3tXO7DgiIqY4k1/MzC/KZk093ieUUH9PkxOJiJhj0Yb97M/Ow8/ThQmDwsyOI1LnqMF3nf3tqz1kni2kha87T/XVrCkRsU2Hc/NZEL8XgEl3htHQy8XkRCIi5pi7KpWcC8W08vdkbO8Qs+OIiJgi/eR5FiXsA2DGkHC83Z1MTiRS96jBdx3tOHyad7dmADBnWARuzponICK2xzAMJi5PorDESrfgBtzbKcDsSCIiptiyP4ePfjgKwLzYSJwd9dZbRGyP1WowPi6JEotB37b+3BHZ2OxIInWS3mVcJ8WlViZ8moRhQGzHZvRs1dDsSCIipojbcYxN+3JwdiybNWVnp1lTImJ7CkssTIxLAmBEdAs6t/Q1OZGIiDk+3H6YHw+dxsPZgVlDI/TeUOQGUYPvOlmycT97Tp7H18OZyXdq1pSI2KbcC0XMXpkCwJN9WxHk52FyIhERc/xjXToZufk0qufCuIFtzY4jImKKE2cLeX51GgB/vb0NTX3cTE4kUnepwXcdHMi+wD/Wl80TmHpXO3w9nE1OJCJijtkrUzmdX0Lbxl6MvjXY7DgiIqZIzTzHko0HAJgZE0E9V82aEhHbNO3zZC4UldKhuQ8PdGtpdhyROk0Nvt/o4qyp4lIrt7ZuSEyHpmZHEhExxTd7s1n+0zHs7GD+3VE4OehHjIjYHst/Z02VWg0Ghjfm9nDNmhIR27Qm+QRf7T6Jo70d8++OxMFel+aK3Ej67es3+uiHI3x34BRuTg7M0TwBEbFR+cWlTFpeNmtqVPeWdGjuY24gERGTvLc1g11HzuDl4siMmHCz44iImOJcYQlTP0sG4I+9gmnbuJ7JiUTqPjX4foOs84XMWZkKwDP9W9Pc193kRCIi5ljw9V6Oni6gqbcrfxnQxuw4IiKmOHamgBe/2gPAc4Pa0qieq8mJRETM8fzqNLLOFxHk58Gfb2tldhwRm6AG328w44sUzhWWEtnMm4d6tDQ7joiIKZKOnuXNTQcBmD0sAk8XR5MTiYhUP8MwmLIimfxiCze3rM//dWlhdiQREVN8n3GKD7cdBmDusEhcnRxMTiRiG9Tgu0brUk+yMjETB3s75sVG4qhZUyJig0otVsbHJWI14K6oJtzWtpHZkURETLEyKZP1aVk4OZS9N7TXrCkRsUFFpRYmxJWNbfld5+Z0C2lgciIR26Gu1DW4UFTKlBVl8wQeuSWIiGbeJicSETHHW5sPsvv4Oeq5OjJtsGZNiYhtOptfwvTPUwB4rHcoof5eJicSETHH4oT97Mu6gJ+nMxPuaGt2HBGbogbfNfjbV3s4fraQ5r5uPN2vtdlxRERMcTg3n5e/3gvApDvDaOjlYnKiqmVmZlKvXj1ef/11s6OISB01d1UqOReKCPX35LE+IWbHERExxb6s8yzasB+AaYPD8XF3NjmRiG1Rg+8q/XT4NO9uzQDK5gm4OWuegIjYHsMwmLQiicISK9HBvgzv3NzsSJc0fvx4+vfvb3YMEamjtu7PZdkPRwCYFxuJi6PeG4qI7bFaDcZ/mkSxxcptbf25K6qJ2ZFEbI4afFehxGJlQlwShgGxNzWjZ6uGZkcSETHFip3H+DY9B2dHe+bFRmFnd2WzpiZPnkybNm3o378//fv3JyEhgZUrVxIdHU2HDh0YM2YMVquVjIwMOnbsyP3330/btm0ZO3Zs+T6qWn8pGzdupF69ekRGRv7mYxYR+aXCEgsTl5fNmrq/awtubulrciIREXP8+/vD/HDoNB7ODswaGnHF7w1F5PpRg+8qLNl4gLQT56nv7sTku9qZHUdExBSn8oqZ9WUqAE/2bUWQn8cVPW/79u2sX7+e5ORk3n33XbZt20ZOTg4LFy4kISGBnTt34uDgQFxcHAApKSlMnz6d3bt3s3XrVtLT0y+7/pdKS0uZOnUqM2bMuD4HLiLyC6+u38fBnDz8vVx4bpBmTYmIbTp5rpD5q9IAePb2NjTzcTM5kYhtcjQ7QG1xMCePv69LB2DKXe3w9dA8ARGxTbO/TOFUXjFtGnkx+tbgK37eli1bGDZsGE5OTjRt2pSePXvi4uJCYmIi0dHRABQUFBAYGEjnzp0JCwujVatWAERFRZGRkUFaWlqV66vy2muvce+99+LrqzNqROT6Sztxjte/KZs1NTMmnHquTiYnEhExx7TPdnO+qJT2zX14sFtLs+OI2Cw1+K6AYRhMjEuiuNRKz1Z+DLupmdmRRERMsXFvNnE/HcPODubfHYmTw5WfCG4YRpXbYmJiWLJkSYXtGRkZuLj876Yd9vb2WCyWS66vyvbt29m0aRMvvvgiZ86cwcHBAXd3dx588MErziwiUhXLf2dNlVoNBrRrxMAIzZoSEdv01e4TrNl9Akd7O+bHRuJgr0tzRcyiS3SvwMc/HGXrgVxcneyZMzRS8wRExCblF5cyaUXZrKmR3VpyU4v6V/X87t27s2LFCkpLS8nMzGTTpk0UFhaybt06jh49CkBubm7536sSHR19xes//PBDDh06REZGBk899RRz5sxRc09Erov3t2aw88gZvFwcmRkTYXYcERFTnCssYepnyQCMvjWYsCb1TE4kYtt0Bt+vyD5fxJxVZbOmnunfmhYN3E1OJCJijr/Hp3PkVAFNvV159vY2V/38rl270qtXL8LDwwkMDKR9+/aEhoby2muvERMTQ0lJCU5OTixduvSSl9X6+/tXuT4gIOC3Hp6IyBU5fqaAF7/aA8C4QW1p7O1qciIREXO8uGYPJ88V0bKBO0/0bWV2HBGbZ2dUdc2UlPvzv3/ii13HCW9aj88e74HjVVyOJiJSVyQfO0vMa5uxWA3++WBn+rVrdE37ycvLw8PDg5ycHLp168aOHTvw8vK6zmlFRG4MwzD4w3s/EJ+aRafA+nz8x27Y63I0EbFBPx46xT2vb8Uw4F9/6Er3ED+zI4nYPJ3Bdxkb0rL4Ytdx7O1gfmyUmnsiYpNKLVbGxyVisRrcGdnkmpt7AA8//DB79uyhpKSEWbNmqbknIrXKqqQTxKdm4eRQNmtKzT0RsUVFpRbGf5qEYcC9nQLU3BOpIdTgu4S8olImryibJ/DILUFEBnibnEhExBxvb84g+dg56rk6Mm1Iu9+0r2XLll2nVP+Tm5tL3759K2zz8/MjPj7+ur+WiNius/klTPt8NwBje4fSqpE+oBAR2/R6wgHSsy7g5+nMpDvDzI4jIv+lBt8l/G3tHo6dKSCgvhtP929tdhwREVMcOZXPy1/vBWDiHWH4e9W8WVMNGjRg586dZscQkTpu3upUci4UEdLQg8f7hJgdR0TEFPuyzvPahn0ATB0cjo+7s8mJROQiXXNahZ1HzvDOlgwA5gyLxN1ZfVARsT2GYTBpRTIFJRa6Bvnyu5ubmx1JRMQU3x3I5T/fHwFgXmwULo4OJicSEal+VqvBxLhkii1W+rRpyOCoJmZHEpGfUYPvF0osVsZ/mohhwNAOTenVuqHZkURETPHZzuNs3JuNs6M982IjsbPTrCkRsT2FJRYmxiUB8PsuLegSVPVdvkVE6rr/fH+E7RmncHd2YNbQCL03FKlh1OD7haXfHiDtxHnquzsx5a7fNmtKRKS2OpVXzMwvUwD4c59Qght6mpxIRMQcr23Yx4GcPPy9XBg/qK3ZcURETJF1rpB5q1MB+MuANgTUdzc5kYj8kq49/ZmMnDz+Hp8OwOQ729HA08XkRCIi5pi9MoVTecW0aeTFH3tp1pSI2KY9J86zOGE/ADOGhOPt5mRyIhGR6pOZmUmbNm144YUX2OV5M+cLS2kf4M2o7i3NjiYiVVCD778Mw2Di8iSKSq30bOVHbMdmZkcSETHFpvQc4nYcw84O5t0dibOjTvYWEdtjtRqMj0uk1GrQv10jBkY0NjuSiEi1Gj9+PP379yfl+FlWF5/Awd6OebFRONjr0lyRmki/tf3XJz8eZcv+XFyd7JkzVLOmRMQ2FRRbmLi8bNbUg9GBdGxR3+REIiLm+GDbIX46fAZPF0dmxoTrvaGImGry5Mm0adOG/v37079/fxISEli5ciXR0dF06NCBMWPGYLVaycjIoGPHjtx///20bduWsWPHlu+jqvWXsnHjRurVq0frsHZ8sSsTgNG3BtOuab0bfqwicm3U4ANyLhQxZ1XZPIGn+rWmRQPNExAR27Rw3V4On8qnibcrfx2oWVMiYpsyzxbwwpo9AIwb2IYm3m4mJxIRW7Z9+3bWr19PcnIy7777Ltu2bSMnJ4eFCxeSkJDAzp07cXBwIC4uDoCUlBSmT5/O7t272bp1K+np6Zdd/0ulpaVMnTqVGTNmsGVfLucKSwhs4M6TfVtV52GLyFXSJbrAzC9SOJNfQrsm9Xj0liCz44iImCL52Fn++e1BAGbGRODpoh8RImJ7DMNgyordXCgqpWMLH0Z0DTQ7kojYuC1btjBs2DCcnJxo2rQpPXv2xMXFhcTERKKjowEoKCggMDCQzp07ExYWRqtWZc24qKgoMjIySEtLq3J9VV577TXuvfdeDp63Y9fRMzh41GfusEhcnRyq54BF5JrY/G9vG/Zk8fmu49jbwfN3R+HooJMaRcT2lFqsTIhLwmI1uCOyMf3bNTI7koiIKVYnnyA+9SRODnbMvzsKe82aEhGTGYZR5baYmBiWLFlSYXtGRgYuLv+7WaS9vT0Wi+WS66uyfft2Nm3aRNaUWRTlX8DFyZH9W7rSI/TB334wInLD2HQ3K6+olMnLkwF4uEcQkQHeJicSETHHO1sySDp2Fi9XR6YPDjc7joiIKc4WlDDt890AjO0VQutGXiYnEhGB7t27s2LFCkpLS8nMzGTTpk0UFhaybt06jh49CkBubm7536sSHR19xes//PBDnn0rnkaj36RRt2HMmT2bBx9Uc0+kprPpM/he/novx84U0MzHjWcGtDY7joiIKY6cyueltXsBmHhHGP71XE1OJCJijvmr08g+X0RwQw8e6xNqdhwREQC6du1Kr169CA8PJzAwkPbt2xMaGsprr71GTEwMJSUlODk5sXTpUnx9favch7+/f5XrAwICKq3dn32BV9bvA6BXm4a4a2yLSK1gZ1R1vq8N2HXkDMMWbcZqwDsP3UzvNv5mRxIRqXaGYTDq7e/5Zm82XVr68p/R0bocrQ7IyMggPDy8fP7OkCFDmDlzpsmpRGq2bQdy+d2S7wBYNjqarsENTE4kIvI/eXl5eHh4kJOTQ7du3dixYwdeXtf/LGOr1eC+pd+x/eAperVuyDsP3ay7iIvUEjbZii+xWBkfl4TVgJgOTdXcExGb9fmu43yzNxtnB3vmxkaquVeHREZG8t1335kdQ6RWKCyxMGF5EgC/79JczT0RqXEefvhh9uzZQ0lJCbNmzbohzT2AZT8cYfvBU7g5OTB7aISaeyK1iE02+P757UFSM8/h4+7ElLvamR1HRMQUp/OKmflFCgB/ui2UUH9PkxPVbZMnT+bjjz+mRYsWAEyaNIm8vDxmzZpFYWEh0dHRLFq0iMOHDxMbG0tYWBg//vgjffr0YfHixQCsXLmy0np7e5sepytyXSzasI8D2Xk09HJh/KAws+OIiFSybNmy677P3Nxc+vbtW/51qdUg44I9fsNn85cBrWnu637dX1NEbhyb+63gUG4eC+PLZk1NvrMdfp4uv/IMEZG6ac6qVHLzimndyJMxvULMjlOnbd++nfXr15OcnMy7777Ltm3byMnJYeHChSQkJLBz504cHByIi4sDICUlhenTp7N79262bt1Kenr6ZddXJSUlhZtuuon+/fuTmJhYXYcqUuvsPXmexd/sB2DGkHC83ZxMTiQiUj0aNGjAzp07y//0eu5t/IbPJirAm4d6BJkdT0Sukk2dwWcYBhOXJ1FUaqVHaAPu7tjM7EgiIqbYvC+HT348ip0dzIuNwtnR5j7vqVZbtmxh2LBhODk50bRpU3r27ImLiwuJiYlER0cDUFBQQGBgIJ07dyYsLKx8fl5UVBQZGRmkpaVVub4qTZo04eDBgzRo0ICNGzdy7733smfPnuo5WJFaxGo1mBCXRInFoF+YP4MiGpsdSUTEFPEpJ1mZlImDvR3zYiNx0NgWkVrHphp8n+44xuZ9ubg42jNnaKTmCYiITSossTDxv7OmRnQNpFNgfZMT1X1V3c/KMAxiYmJYsmRJhe0ZGRm4uPzv7HJ7e3ssFssl11fFxcWlfB+33norzs7O5OTk4Ofn9xuPRKRu+XDbIX48dBoPZwdmxmjWlIjYpvOFJUz5LBmAR3sGEd7U2+REInItbOaUjZwLRcxeWTZr6ql+rWnp52FyIhERcyyMT+dQbj6N67kybmAbs+PYhO7du7NixQpKS0vJzMxk06ZNFBYWsm7dOo4ePQqUzcG5+PeqREdHX/H67OxsLBYLALt37+bChQs0aKCbBoj8XObZAp5fU3Zm67iBbWnq42ZyIhERc/ztqz1kni2kha87T/VtbXYcEblGNnMG3+wvUziTX0JYk3o82lPzBETENqUcP8fSbw8AMDMmHC9XzZqqDl27dqVXr16Eh4cTGBhI+/btCQ0N5bXXXiMmJoaSkhKcnJxYunQpvr6+Ve7D39+/yvUBAQGV1m7cuJGpU6fi5OSEi4sL77//vs5MEvmFaZ/t5kJRKTe18GFEdNWXu4uI1HU7Dp/mve8OATB3WCRuzg4mJxKRa2VnVHXdUB3zzd5sRr61HXs7WP5YD9o39zE7kohItbNYDWIXbWbX0bMMimjM4hGdzI5kU/Ly8vDw8CAnJ4du3bqxY8cOvLy8zI4lYpPWJGcy5oMdONrbsfKJnrRprP8XRcT2FJdaGfzKJvacPM/dHQN4aXh7syOJyG9Q58/gyy8uZdJ/Z02N6h6k5p6I2Kx3tmSw6+hZvFwdmT4k3Ow4Nufhhx9mz549lJSUMGvWLDX3RExytqCEqZ/tBmBMrxA190TEZi3ZuJ89J8/j6+HMpDvDzI4jIr9RnW/wvbx2L0dPF9DMx42/DNA8ARGxTUdP5/PS2rJZU+MHtaVRPVeTE9meZcuWXfd95ubm0rdv3wrb/Pz8iI+Pv+6vJVJXPL8mjazzRQT7efCn20LNjiMiYooD2Rf4x/p9AEy9qx2+Hs4mJxKR36rONfgsVoPtB0+Rdb6Q84WlvLX5IACzh0Xg4VLnDldEpEo/r4UNPV14Y+MB8ostdGnpy+9vbmF2PLlOGjRowM6dO82OIVKj/bwenrpQzL+2HQZgbmwkrk6aNSUituHntdDP04V/rEunuNTKra0bEtOhqdnxROQ6qFMdrzXJmcz4IoXMs4UVtt/csj592viblEpEpHpdqhY62tsxNzYSe3vdbEFEbMOl6mGP0AZEB+vO0iJiGy5VC50d7JkzNEI34hKpI+zNDnC9rEnOZOwHOyoVLYDvM06zJjnThFQiItXrcrWw1GqwL+u8CalERKrf5erh5n25em8oIjbhcrWw2GJl9/GzJqQSkRuhTjT4LFaDGV+kcKnbAdsBM75IwWKt8zcMFhEb9mu1EFQLRcQ26L2hiIjeG4rYmjrR4Nt+8FSVn0hcZACZZwvZfvBU9YUSEalmv1YLQbVQRGyD3huKiOi9oYitqRMNvqzzly9aV7tORKQ2Ui0UESmjeigiolooYmvqRIPP38v1uq67nv7whz9gZ2fHXXfdVe2vLSK2pSbXQhGR6qR6KCKiWihia+pEg69LkC9NvF251L1/7IAm3q50CfKtzlj88MMPvPPOO7i6qmCKyI1XU2uhiEh1Uz0UEVEtFLE1daLB52Bvx7TB7QAqFa+LX08b3A4H++q7/bdhGDzxxBM8+OCDNGrUqNpeV0RsV02shSIiZlA9FBFRLRSxNTWuwZeQkEDnzp1xdXUlJCSEN954g+nTp2NnV7HofPDBB3Tq1Ak3Nzd8fX15Z/bTzOzXiMbe/ztb7sS/xnPy7T8xvpsnz//5/3B3d6dZs2a88MILlV63qKiIadOmERoaiouLC82bN2fcuHEUFRVd03G8//77JCcnM2fOnGt6vojItRgY0YTFIzpWqIUAjb1dWTyiIwMjmpiUTESkeqkeioioForYEjvDMGrMPbF/+uknunXrRpMmTRgzZgwWi4XXXnuNhg0bsmvXLi5GnTNnDlOmTGH48OH06tWL7OxsXnnlFTw9Pfnhxx3sPW0l63whM/44nJNHM3BwcCA2NpY2bdrwySefsH79elatWsWgQYMAsFqtDBo0iE2bNjF69GjCwsJISkri9ddf584772TFihVXdRznz5+ndevWPPnkk4wfP56WLVsSERHBl19+eb2/ZSIiVbJYDbYfPEXW+UL8vcouvdCnsyJii1QPRURUC0VsQY1q8A0ZMoR169aRnp5O06ZNAdi3bx9hYWGUlpZiGAaHDh0iJCSEmTNnMnHixPLnJicnc9NNNzFjxozy7b179+abb77hvffe44EHHgCguLiYwMBAevTowSeffAKUnQ04cuRIvvnmG2655Zbyfb7xxhuMGTOGzZs307179ys+jr/+9a98+umnpKam4uLiogafiIiIiFyVzMxM2rRpwwsvvMCYMWPMjiMiIiI1XI25RNdisRAfH8/QoUPLm3sAoaGh5WfaAcTFxWG1Whk+fDg5OTnlfxo3bkyrVq3YsGFDhf16enoyYsSI8q+dnZ3p0qULBw4cKN/28ccfExYWRtu2bSvs87bbbgOotM/L2bt3L3//+9958cUXcXFxuervg4iIiIjI+PHj6d+/v9kxREREpJZwNDvARVlZWRQUFBAaGlrpsZ9vS09PxzAMWrVqVeV+nJycKnwdEBBQaX5f/fr1SUxMrLDP1NRUGjZseMlsV+rJJ5+ke/fu3H333Vf8HBERERGpuSZPnszHH39MixYtAJg0aRJ5eXnMmjWLwsJCoqOjWbRoEYcPHyY2NpawsDB+/PFH+vTpw+LFiwFYuXJlpfX29lV/1r5x40bq1atHUFBQtR2jiIiI1G41psF3paxWK3Z2dqxevRoHB4dKj3t6elb4uqo1AD+/MtlqtRIZGcnLL79c5drmzZtfUbb169ezZs0a4uLiyMjIKN9eWlpKQUEBGRkZ+Pr6Uq9evSvan4iIiIiYa/v27axfv57k5GSys7PLr/h44403SEhIwNXVlccff5y4uDg6d+5MSkoKy5YtIzg4mE6dOpGenk79+vVZuHBhpfX33HNPpdcrLS1l6tSpxMXF8Y9//MOEIxYREZHaqMY0+Pz9/XF1dWXfvn2VHvv5tpCQEAzDICgoiNatW1+X1w4JCWHXrl307du30tl+V+Pw4cMAxMbGVnrs2LFjBAUFsWDBAp566v/bu/uYKuv/j+MvzvkeD0J8TSADE5nclOQdKlNuvA1t9IeALm+2UtGW07WWNWaZN3GTrlXebGlOWTeW/RGmslqiDkSdqeFEI0HYcXIskyYc0ykDw8P5/eGPs4iDhfrlcMnzsblxfc7bz/X+MH3/8T6f6/osved7AAAAoOscO3ZM06dPl8ViUf/+/TV+/HhZrVaVl5crPj5ektTY2Kjw8HDFxcUpJibG/aTJ8OHDZbfbVVVV5THek82bN2vmzJkKDAzsmgUCAICHQrdp8JnNZk2ZMkUFBQW6fPlym0M2CgsL3XEzZszQ8uXLlZ2drR07drRpyLlcLl29elVBQUGduvesWbO0d+9e5eXladGiRW0+a2xsVEtLi/z9/f9xnmeeeUZ79uxpN75o0SKFh4drxYoVGjZsWKdyAwAAgPd4Oo/O5XIpLS1N27ZtazNut9vbvIPZZDLJ6XR2GO9JaWmpjh49qg8++EDXrl2T2WyWn5+f5s2bd/+LAQAAD61u0+CTpKysLB04cEBJSUlasmSJnE6nNm3apKFDh+rMmTOS7uy2e/fdd7V8+XLZ7Xalp6crICBANTU12rNnjxYtWqTMzMxO3Xfu3LnKz8/X4sWLVVJSoqSkJDmdTlVVVSk/P1/79+9XXFzcP84zcOBA97tZ/mrp0qV6/PHHlZ6e3qm8AAAA4F2JiYl644039Prrr6uurk5Hjx7V/PnzVVxcrEuXLmnAgAFyOBxqbGzscI74+Hi9/vrr7eIHDBjQLvarr75y/5yVlaWQkBCaewAA4B91qwbf6NGjVVhYqMzMTK1atUphYWHKycnRuXPnVFVV5Y5766239OSTT2rDhg3Kzs6WdOc9ec8++6xSU1M7fV+TyaSCggJt2LBBX3zxhfbs2SM/Pz9FRETotddee2CPAgMAAMBYxo4dq4kTJ2rIkCEKDw/XiBEjFBUVpc2bNystLU3Nzc2yWCzKy8vr8LHafv36eYz31OADAAC4Fz4uT88ddDPp6emqqKiQzWbzdioAAADoYRoaGuTv76/6+nolJCSorKxMAQEB3k4LAADArVvt4JPuvPOud+/e7mubzaa9e/dq/vz5XswKAAAAPdXChQtVXV2t5uZm5ebm0twDAADdTrfbwRcaGqqMjAxFRETo4sWL2rJli27duqXTp0+7TyTzhuvXr9/13SqSFBIS0kXZAAAAwMgcDoeSk5PbjAUHB6uoqMhLGQEAACPrdg2+BQsWqKSkRL///rusVqsSEhK0du1ajRo1yqt5ZWRkaPv27XeN6Wa/SgAAAAAAAPQA3a7B111VVlbq8uXLd42ZMmVKF2UDAAAAAAAA3EGDDwAAAAAAADAwk7cTAAAAAAAAAHDvaPABAAAAAAAABkaDDwAAAAAAADAwGnwAAAAAAACAgdHgAwAAAAAAAAyMBh8AAAAAAABgYDT4AAAAAAAAAAOjwQcAAAAAAIB/zW63y9/fX7GxsYqNjdXq1au9nVKP9x9vJwAAAAAAAABjGTZsmE6cOOHtNPD/2MEHAAAAAABgcCtXrtRTTz2lqVOnaurUqTp06JC+//57xcfHKzY2VosXL1ZLS4vsdrtGjRqlF154QYMHD9aSJUvcc3iKhzHQ4AMAAAAAADCw0tJSHTx4UGfPntX27dv1448/qr6+Xhs3btShQ4d05swZmc1m7d69W5JUWVmprKwsVVRU6Pjx47LZbHeN96SyslIjR47U1KlTVV5e3lVLRQd4RBcAAAAAAMDAjh07punTp8tisah///4aP368rFarysvLFR8fL0lqbGxUeHi44uLiFBMTo+joaEnS8OHDZbfbVVVV5THek9DQUNXU1CgoKEhHjhzRzJkzVV1d3TWLhUc0+AAAAAAAAAzM5XJ5HEtLS9O2bdvajNvtdlmtVve1yWSS0+nsMN4Tq9XqnmPChAnq1auX6uvrFRwcfJ8rwb3iEV0AAAAAAAADS0xMVEFBgW7fvq3a2lodPXpUTU1NKi4u1qVLlyRJDofD/bMn8fHx/zq+rq5OTqdTklRRUaGbN28qKCjoAa8KncEOPgAAAAAAAAMbO3asJk6cqCFDhig8PFwjRoxQVFSUNm/erLS0NDU3N8tisSgvL0+BgYEe5+jXr5/H+AEDBrSLPXLkiFavXi2LxSKr1aovv/xSPj4+/+tl4i58XJ72cQIAAAAAAMAwGhoa5O/vr/r6eiUkJKisrEwBAQHeTgtdhB18AAAAAAAABrdw4UJVV1erublZubm5NPd6GHbwAQAAAAAAoB2Hw6Hk5OQ2Y8HBwSoqKvJSRugIDT4AAAAAAADAwDhFFwAAAAAAADAwGnwAAAAAAACAgdHgAwAAAAAAAAyMU3RheM4Wl0prrurKjSb1C/DVmEGBMpt8vJ0WAHQ56iEAUAsBoBX1sGehwQdD23e2VtnfVar2epN7LLSPr96Z9rRShoZ6MTMA6FrUQwCgFgJAK+phz8MpujCsfWdrtWRHmf7+D7j1+4gtL46icAHoEaiHAEAtBIBW1MOeiXfwwZCcLS5lf1fZrmBJco9lf1cpZwv9awAPN+ohAFALAaAV9bDnosEHQyqtudpmq/HfuSTVXm9Sac3VrksKALyAeggA1EIAaEU97Llo8MGQrtzouGDdSxwAGBX1EACohQDQinrYc9HggyH1C/B9oHEAYFTUQwCgFgJAK+phz0WDD4Y0ZlCgQvv4qqMDvn1054SgMYMCuzItAOhy1EMAoBYCQCvqYc9Fgw+GZDb56J1pT0tSu8LVev3OtKdlNnVU1gDg4UA9BABqIQC0oh72XDT4YFgpQ0O15cVRCunTdmtxSB9fjv0G0KNQDwGAWggAraiHPZOPy+XibGQYmrPFpdKaq7pyo0n9Au5sNTbCtxGff/65FixY4PGz2tpahYSEdHFGAIzOqPUQAB4kaiEA3EE97Fn+4+0EgPtlNvkoITLI22ncs5ycHA0aNKjN2KOPPuqdZAAYmtHrIQA8CNRCALiDetiz0OADvOy5555TXFyct9MAAAAAAAAGxTv4YFiHDh1SXFycfH19FRkZqa1btyorK0s+Pm23HO/YsUOjR49W7969FRgYqDlz5ujXX39tEzNp0iQNHTpUlZWVmjx5svz8/PTEE0/o/fffb3ffW7du6Z133lFUVJSsVqvCwsK0bNky3bp1657XcuPGDTmdznv++wAAAAAAoOeiwQdDOn36tFJSUuRwOJSdna2XXnpJOTk5KigoaBO3Zs0azZs3T9HR0Vq/fr2WLl2q4uJiTZgwQdeuXWsT+8cffyglJUUjRozQunXrNHjwYL355psqLCx0x7S0tCg1NVUffvihpk2bpo8++kjp6enasGGDZs+efU9rmTx5sv773//Kz89Pqampstls9zQPAAAAAADomThkA4aUmpqq4uJi2Ww29e/fX5J0/vx5xcTE6Pbt23K5XLp48aIiIyOVk5Ojt99+2/13z549q5EjRyo7O9s9PmnSJB0+fFhffPGF5s6dK0n6888/FR4erqSkJH3zzTeS7uwGnD9/vg4fPqxx48a559y6dasWL16sH374QYmJif9qDfn5+SosLHQ3+E6dOqX169fLz89PZWVlCgsLeyC/KwAAAAAA8HBjBx8Mx+l0qqioSOnp6e7mniRFRUXpueeec1/v3r1bLS0tmjVrlurr691/QkJCFB0drZKSkjbzPvLII3rxxRfd17169dKYMWN04cIF99jOnTsVExOjwYMHt5nzmWeekaR2c97NrFmz9Nlnn2nevHlKT09Xbm6u9u/fL4fDoTVr1nT69wIAAICewW63y9/fX7GxsYqNjdXq1au9nRIAwMs4ZAOGc+XKFTU2NioqKqrdZ38ds9lscrlcio6O9jiPxWJpcz1gwIB27+/r27evysvL28x57tw5PfbYYx3mdj/GjRunsWPHqqio6L7mAQAAwMNt2LBhOnHihLfTAAB0EzT48NBqaWmRj4+PCgsLZTab233+yCOPtLn2FCNJf32KvaWlRcOGDdP69es9xj6Ix2rDwsJUXV193/MAAACg66xcuVI7d+7UwIEDJUkrVqxQQ0ODcnNz1dTUpPj4eH388cf65ZdfNGPGDMXExOjUqVOaPHmytmzZIkn6/vvv28WbTDx0BQD4ZzT4YDj9+vWTr6+vzp8/3+6zv45FRkbK5XJp0KBBevLJJx/IvSMjI/XTTz8pOTm53W6/B+XChQsd7hAEAABA91NaWqqDBw/q7Nmzqqurc7/OZevWrTp06JB8fX31yiuvaPfu3YqLi1NlZaW+/vprRUREaPTo0bLZbOrbt682btzYLv7555/3eM/KykqNHDlSwcHBWrdunYYPH97FqwYAdCd8HQTDMZvNmjJligoKCnT58mX3+Pnz59uceDtjxgyZzWZlZ2fr72fJuFwuORyOTt971qxZ+u2335SXl9fus8bGRjU0NPzruerq6tqN7d27V6dOnVJKSkqncwMAAIB3HDt2TNOnT5fFYlH//v01fvx4Wa1WlZeXKz4+XrGxsSoqKnK/2zkmJkbR0dEym80aPny47Ha7jh8/3mH834WGhqqmpkanT5/WqlWrNHPmzK5cLgCgG2IHHwwpKytLBw4cUFJSkpYsWSKn06lNmzZp6NChOnPmjKQ7u+3effddLV++XHa7Xenp6QoICFBNTY327NmjRYsWKTMzs1P3nTt3rvLz87V48WKVlJQoKSlJTqdTVVVVys/P1/79+xUXF/ev5kpMTNTIkSMVFxenPn36qKysTJ9++qnCwsLanPoLAACA7u3vXya3jqWlpWnbtm1txu12u6xWq/vaZDLJ6XR2GO+J1Wp1zzFhwgT16tVL9fX1Cg4Ovs+VAACMih18MKTRo0ersLBQffv21apVq/TJJ58oJydHycnJ8vX1dce99dZb2rVrl0wmk7Kzs5WZmalvv/1Wzz77rFJTUzt9X5PJpIKCAr333nv6+eeflZmZqezsbJ08eVKvvfZapx4Fnj17tmw2m9auXatXX31V+/bt08svv6yTJ0/q8ccf73RuAAAA8I7ExEQVFBTo9u3bqq2t1dGjR9XU1KTi4mJdunRJkuRwONw/exIfH/+v4+vq6uR0OiVJFRUVunnzpoKCgh7wqgAARuLj8vR1E2BQ6enpqqiokM1m83YqAAAA6EHefvtt7dq1S+Hh4WpqatLGjRt15coVrVixQs3NzbJYLMrLy1NgYKDmzJnjPgE3IyNDc+bMUUpKivbt29cuftSoUe3utWvXLq1evVoWi0VWq1Xr1q3TuHHjunrJAIBuhAYfDKuxsVG9e/d2X9tsNg0ZMkTz58/3+I48AAAA4H+loaFB/v7+qq+vV0JCgsrKyhQQEODttAAAPQTv4INhRUREKCMjQxEREbp48aK2bNmiXr16admyZV7N6/r162psbLxrTEhISBdlAwAAgK6wcOFCVVdXq7m5Wbm5uTT3AABdih18MKwFCxaopKREv//+u6xWqxISErR27VqPjzF0pYyMDG3fvv2uMfy3AwAAwD9xOBxKTk5uMxYcHKyioiIvZQQA6K5o8AEPWGVlpS5fvnzXmClTpnRRNgAAAAAA4GFHgw8AAAAAAAAwMJO3EwAAAAAAAABw72jwAQAAAAAAAAZGgw8AAAAAAAAwMBp8AAAAAAAAgIHR4AMAAAAAAAAMjAYfAAAAAAAAYGA0+AAAAAAAAAADo8EHAAAAAAAAGBgNPgAAAAAAAMDAaPABAAAAAAAABkaDDwAAAAAAADAwGnwAAAAAAACAgdHgAwAAAAAAAAyMBh8AAAAAAABgYDT4AAAAAAAAAAOjwQcAAAAAAAAYGA0+AAAAAAAAwMBo8AEAAAAAAAAG9n9z5sjlyAj4bAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1, 4, figsize=(16, 4))\n", "\n", "arguments = [\n", " {\"fontsize\": \"large\"}, # default :\"x-small\"\n", " {\"annot_radius\": 0.5}, # default: 1.1\n", " {\"annot_xloc\": 2}, # default: -1\n", " {\"annot_yloc\": 2} # default: -0.4\n", "]\n", "\n", "for i, v in enumerate(arguments):\n", " grnplot(models[0], ax=ax[i], **v)\n", " ax[i].set(title=list(v.keys())[0])" ] } ], "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.12.0" } }, "nbformat": 4, "nbformat_minor": 5 }