{ "cells": [ { "cell_type": "markdown", "id": "9dafcc21-b668-4831-a68a-c3b7096b81ce", "metadata": {}, "source": [ "# Analysis of the Error per Age Bin" ] }, { "cell_type": "code", "execution_count": 40, "id": "952a9552-f4a6-41af-a63c-1f59ff1b16ae", "metadata": { "tags": [ "hide-input" ] }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
subjectIDgenderagey_hatmodality
0sub-474134.0140.339008DF+GM
1sub-292123.7331.652921DF+GM
2sub-382165.9146.452941DF+GM
3sub-434167.2456.620143DF+GM
\n", "
" ], "text/plain": [ " subjectID gender age y_hat modality\n", "0 sub-474 1 34.01 40.339008 DF+GM\n", "1 sub-292 1 23.73 31.652921 DF+GM\n", "2 sub-382 1 65.91 46.452941 DF+GM\n", "3 sub-434 1 67.24 56.620143 DF+GM" ] }, "execution_count": 40, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import os\n", "import math\n", "import pandas as pd\n", "import seaborn as sns\n", "import matplotlib.pyplot as plt\n", "from pathlib import Path\n", "\n", "path_files_test = os.path.join('..', '..', '..', 'resources')\n", "\n", "all_data = []\n", "for path in Path(path_files_test).glob(\"*test.csv\"):\n", " if \"gm_wm\" in path.name:\n", " continue\n", " df_mod = pd.read_csv(path)\n", " df_mod[\"modality\"] = \"+\".join(path.name.split(\"_\")[:-1]).upper()\n", " all_data.append(df_mod)\n", "df = pd.concat(all_data)\n", "df.head(4)" ] }, { "cell_type": "code", "execution_count": 25, "id": "dd77cad1-e113-40e7-9d92-850547d29889", "metadata": { "tags": [ "hide-input" ] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEiCAYAAAD+lrlwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAACQdElEQVR4nOydd3hUxfeH39matum99wRC6L1JEZAu2BXs8hMs2L4WREVF7F0siB0VlSKIgqAogvReQiiBQDrpPVvn98eGQEghQEKR+z5PHrKzd+aeu2Q/d+6ZM+cIKSUKCgoKCpcPqgttgIKCgoLC+UURfgUFBYXLDEX4FRQUFC4zFOFXUFBQuMxQhF9BQUHhMkMRfgUFBYXLDM2FNqA+tmzZ4qvRaGYDbVBuTgoKCgqNYQN2WyyWuzt16nSsKR0uSuHXaDSz/f39W/n4+BSqVCplo4GCgoJCA9hsNpGbm9s6Ozt7NjCqKX0u1tl0Gx8fnxJF9BUUFBQaR6VSSR8fn2LsHpKm9WlBe84FlSL6CgoKCk2jWi+brOcXq/BfcNRqdaf4+PjW0dHRCXFxca2fe+45P6vVCsCSJUsMBoOhfXx8fOv4+PjWPXv2jG3quLt27dL3798/OiQkpE1CQkKrbt26xS5dutQF4L333vMSQnT6+eefDceP/+abb9yFEJ2++OILj2a/yAvE0aNHNSNGjIg8/hlcccUV0Tt37tTffvvtITExMQmxsbGt27Rp0yo5OVkHEBQUlBgbG9v6+Oe9YsUK5wt9Dc3BXXfdFfLCCy/4Hn/du3fvmBtuuCHs+Ot77rkneNq0aX5CiE4PPvhg4PH2rKwsjUaj6XjrrbeGnm+bW4qW+r5d7KSlpWlGjhwZERwcnJiQkNCqffv28V9//bX7kiVLDEKITm+99Zb38WPXrl3rKITo9Oyzz/qd63kvSh//xYBer7clJycnAWRkZGiuu+66yJKSEvXbb7+dCdC5c+eyv/7662BD/bt27Rr3zTffHI6LizMdb6uoqBAjR46Meemll9JuueWWYoBNmzY5rFu3znno0KFlADExMZXff/+959VXX10KMHfuXM+4uLjKlrzW84nNZmPUqFHRN998c/6SJUsOAaxbt87xm2++8czOztYmJyfvUavVpKSkaF1dXW3H+61atWp/QECA5cJZ3vz07t277KeffvIAjlmtVgoLCzVlZWXq4+9v2rTJ5cYbb0wLCgoyrVixwh3IBPj66689oqOjqy6Q2S3CuX7fLkVsNhsjR46Mvvnmm/N/+eWXwwD79+/X/fTTT+6enp6VMTExlfPnz/d45JFH8gC++eabZtOC/8SMf862Ms8uH2Qnhr+a2anLB9mJc7aVeTbn+EFBQZbZs2enfvHFF742m+30HRrgk08+8erYsWPZcdEH6NKlS9WDDz6Yf/x1t27dyrZt2+ZsNBpFcXGxKjU1VZ+QkFBxjpdwdmTO8mR9eCL/OHZifXgimbPO+XNdsmSJQaPRyMcffzz3eFuPHj0qnZ2dbX5+fma12q57UVFRZh8fH+u5nq+5+HjpYs/A265JVI3q3ynwtmsSP166+Jw/i/79+5dt3brVBWDLli2OcXFxlc7Oztbc3Fx1ZWWlSElJcfDx8bE4OjraoqOjK//55x8ngPnz53teffXVBed6/rNl555fPGd/c13iu58M7DT7m+sSd+755aL8vjUnKZ9/7rmkVavEeV5enZa0apWY8vnn53zNv/zyi0Gr1db6LsTGxpqefvrpYwBBQUEmo9GoSktL09hsNlauXOk2cODA4oZHbDqX/Ix/zrYyzxf+LAkzWu03sWPlNt0Lf5aEAYzr4NJsX47WrVubrFYrGRkZGoDNmze7xMfHtwYYPXp0wauvvpp9ujH27Nnj0KFDh0ZFXAhB3759SxYsWOBaVFSkvuqqq4pSU1P1zXMVZ0DmLE8OPR6Grco+OTBl6zj0uN0NETjhrD/XnTt3OrZr167OZzB+/PiCvn37xsfHxxv69OlTcvvtt+f36tWrZnZzxRVXxKpUKnQ6nW3nzp3JZ3v+s+HjpYs9H549M6zKbFIBZBUW6B6ePTMM4N6ho876swgPDzer1Wp54MAB3apVq5y7d+9enpGRoV25cqWLh4eHJTY2tlKv10uAG2+8sWDOnDmegYGBZrVaLQMDA82ZmZm65rnCprNzzy+e/6z7MMxqtX8W5RUFun/WfRgG0DZh5EX1fWsuUj7/3HPH00+H2YxGFUBVTo5ux9NPhwFE3XnnWV/zrl27HNu2bduoHlx99dWF33zzjUfnzp0rEhMTK47/PZwrl/yM/91/y4KOi/5xjFZU7/5bFtSS5+3cuXNZcnJyUnJyctLxP8J3333X67gfcvfu3U7Dhg2LiY+Pbz1o0KCo+sYYNGhQVExMTMLgwYNrvX/LLbcUzJ071/Onn37yvO222y7MzO7ojKAa0T+OrUrF0Rkt8rlGRUWZDx48uPuFF15IV6lUDBs2LG7RokU1ax2rVq3an5ycnHS+RR/ghblfBR0X/eNUmU2qF+Z+dc6fRadOncr++usv53Xr1rn06dOnrGfPnuX//vuv8+rVq126detWdvy4a665pmTVqlWuX331lec111xzwWb7G7d+HXRc9I9jtZpUG7d+fd6/b+eLva+/HnRc9I9jMxpVe19/vVmvefz48aFxcXGt27Rp0+p426233lqwcOFCzzlz5njdfPPNzfb/fskLf265rd5ZT0PtZ0tSUpJOrVYTFBTUoJ958uTJ+cf/ONu0aVPx22+/HUhOTk5asWJFCkBCQkLVtm3bnI4fv2LFipTPPvvscFFRUa0nr/79+1ckJyc7FhQUaNq2bWtszutoMqac+j+/htqbSGJiYuWOHTuc6nvP0dFRXn/99SWffPJJ+uTJk7MWLFjgfi7nai6yCwvqveaG2s+Enj17lq1du9YlOTnZsUuXLpX9+vUr27Rpk8v69etdevXqVSP8Dg4Osm3bthUfffSR//jx4wvP9bxnS3lF/dfcUPvZ0pTv2/mi6tixeq+tofamkpiYWLlz586a78I333xz9O+//95fWFhYowehoaEWrVYr//nnH9dRo0aVnMv5TuaSF34fZ5XpTNrPhszMTM0999wTdscddxxTqc7+I7vnnnvyN2/e7PLtt9+6HW8rLy+vd8AXX3wx/cUXX8w465OdKzq/+j+/htqbyMiRI0tNJpN44403aqIVNmzY4Pjrr7+6pKamagGsViu7du1yDAsLa7b/w3PB38OzXjsaaj8T+vbtW/bHH3+4u7u7WzUaDX5+ftaSkhL1tm3bXAYMGFB+8rFPPPFE9nPPPZfu5+d3wdY+nJ3qv+aG2s+G5vq+NRcOvr71XltD7U1l5MiRpUajUbz66qs+x9vKysrqXPDzzz+f8eKLL6ZrNM3nmb/wn+o5MrmXS4ZeTa0VIL0a2+ReLuckmkajUXU8vKx///6xAwcOLHnjjTcyz2VMFxcXuWjRooOzZs3yCQ4OTmzfvn389OnTA6ZMmVJn3Ouvv75k5MiRpedyvnMidEoGKofaK2sqBxuhU87pc1WpVCxevDhl5cqVriEhIW2io6MTnnjiiaDt27c7DR8+PDomJiYhPj4+QaPR8OSTTzZp+3lL8+yNt2U4aHW1PgsHrc727I23nfONuWvXrpVFRUWazp0718zu4+PjK11cXKynRjF17ty56oEHHsivO8r5o2vHWzPU6tqfhVqts3XteOtF931rLlr9738ZKr2+1jWr9Hpbq//975y/C7/88kvK6tWrDUFBQYmJiYmtxo0bFz5t2rT0k48bNGhQ+fjx44vO5VynIi7G0os7duxIbdeuXV5Tj5+zrczz3X/LgnLLbTofZ5Vpci+XjOZc2L1syZzlydEZQZhydOj8TIROyTiXhd1LmY+XLvZ8Ye5XQdmFBTp/D0/TszfelnEuC7uXMjv3/OK5cevXQeUVBTpnJ09T1463ZjTnwu7FSMrnn3vuff31oKpjx3QOvr6mVv/7X8a5LOy2BDt27PBu165deFOO/U8Iv4KCgsLlzpkI/yXv6lFQUFBQODMU4VdQUFC4zLhYhd9ms9nEhTZCQUFB4VKgWi+bvM35YhX+3bm5uW6K+CsoKCg0TnU+fjdgd1P7XJQpGywWy93Z2dmzs7OzlQpcCgoKCo1TU4GrqR0uyqgeBQUFBYWWQ5lNKygoKFxmKMKvoKCgcJmhCL+CgoLCZYYi/AoKCgqXGYrwKygoKFxmKMKvoKCgcJmhCL+CgoLCZcZFuYHrVLy9vWV4ePiFNkNBQUHhkmLLli15UkqfU9svCeEPDw9n8+bNF9oMBQUFhUsKIcSR+toVV4+CgoLCZYYi/AoKCgqXGYrwKygoKFxmKMKvoKCgcJmhCL+CgoLCZUaLCr8QYrIQYrcQYo8Q4qHqNk8hxAohxIHqfz1a0gYFBQWFRpESynZAzreQOx8qUi60RS1Oiwm/EKINcA/QFWgHjBBCRANPAn9KKWOAP6tfKygoKFwYStbDtr6w7y7Yewvsugoq9l9oq1qUlpzxtwI2SCkrpJQWYBUwFhgNfFV9zFfA1S1og4KCgkLDWI1w9HWQxhNtxjQo+vuCmVSD6Zj9BmQubvahW1L4dwN9hBBeQggnYBgQAvhJKbOqj8kG/FrQBgUFBYWGsVVBVT2unaqj59+W40gJhX/B9itgc1vYPcruimpGWkz4pZR7gVeB5cAyYDtgPeUYCdRb+1EIMUEIsVkIsTk3N7elzFRQULic0bpBQD2laj0GnH9bjlO5H/aMharD9telG2DvbWDOa7ZTtOjirpTyMyllJyllX6AQ2A/kCCECAKr/PdZA31lSys5Sys4+PnVSTSgoKCg0Dz7XQMiToHIGnT/EzQZD1wtnT+UhsFWe0pbcrE8hLZqrRwjhK6U8JoQIxe7f7w5EALcBr1T/u6glbVBQUFBoFH0ghD9rn/mrtKC7wN5nrWfdNpUzqN2a7RQtHcc/XwiRBPwC3CelLMIu+IOEEAeAK6tfKygoKFw4hAocgi+86AM4tYLASbXbot4Ex8hmO0WLzvillH3qacsHBrbkeRUUFBQuWTSuEPYMeF8NpmxwiADnNiBE852i2UZSUFBQUGgetB7g3rfFhldSNigoKChcZijCr6CgoHCZobh6FBQUFC4ySqqs7D1mIbfcRoi7mnhvLXqt4uNXUFBQ+E9SZrLx3toyPt1UXtP25jB3rmnjiGimBV7F1aOgoKBwEbE/11JL9AGe/aOYo0XWBnqcOYrwKygoKFwgLBYLX3zxBenp6TVthZW2OseVmyQlxrrtZ4si/AoKCgpAuclGlbn5xLUxpJQsWLCAtm3bcuedd/LFF1/UvBfqrkZ/ihM+ylNNoKu62c6vCL+CgsJlTWGllfm7K7hmTh63zytgTWoVJmu9uSObjbKyMiZMmADA/PnzmTp1as170V4aPhvrSaDBLvSJfhreHemBl1PzCb+yuKugoHBZszLFyCO/FtW8Xn+0gHm3eNE5WN+s59m8eTOffvopH374IQaDgVWrVhEXF4dGU1uGhRB0D4EfxxRSWlGAh6s/fl5ezWqLMuNXUFC4bCk32pi1sQwPRxU3JlgZGSvRqWF1qvH0nZtIcnIy1113HV26dGHBggXs32+v7pWQkFBH9AHMliq271rIgsX3seKPZ5j3870cSl3bbPaAMuNXUFC4jFEJGBZlJFa3iezDC9BonRg9YDzF2sRzHrukpIRHHnmEL774AicnJ5577jkeeeQRXF1dG+1XUHCEtRtn17y22az8seoNfLxjcHP1P2e7QJnxKygoXMY46lT08tzB/p0zKSnNoqAghaRN02jnlnrWY1qt9rBLJycnNm/ezIMPPsihQ4eYNm3aaUUfoLyyoE6b0VRGlbH5SjAqM34FBYXLFqOxnOS9C+u05+ZsIya0/RmNVVZWxttvv83XX3/N1q1bMRgMbN68uV53TmO4GvxRqdQ4O3nhZggkr/AQWo0DLk7eZzROYygzfgUFhcsWtVqDo6N7nXZHh6YXPTEajbz33ntERkby7LPPkpiYSGlpKcAZiz6Ap3soYzs9R+LeIDw/SqJrXkeGd3ocZ+fmW+BVZvwKCgqXLRqNnq4dx5GRtRMp7TH8DnpXQgI7NKl/bm4uXbp04ciRI/Tr149XXnmFbt26nZNNVdk57Pq/hyhLsReBz9+wAVPyUbzefReNk9M5jX2cli69+DBwN/aC6ruAO4AAYC7gBWwBxkspTS1ph4KCgkJDBAW05brR75GVvRut1pFA/zZ4eYY3eLyUkqSkJBISEvDx8eHqq69m2LBhDBo0qFly6ZTs21cj+sdJmz+fuIcewj0h4ZzHhxZ09QghgoAHgc5SyjaAGrgReBV4W0oZjb0A+10tZYOCgoLC6VCp1AT4taJju+tIbD2iUdFftWoVPXv2pGPHjhw9ai9+/s477zB48OBmS6AmVPXIshDNNj60vI9fAzgKITSAE5AFDADmVb//FXB1C9ugoKCgcE5s27aNoUOH0q9fP9LS0pg5cyaBgYEtci7XuDjcEmuHk4aPH49z5CVQc1dKmSGEeAM4ClQCy7G7doqklJbqw9KBoPr6CyEmABMAQkNDW8pMBQUFhUbJzs6ma9euuLq68sYbbzBp0iQcHR1b7HyO/v50//xzspYvJ3/DBgIGD8a3Xz80Dg7Ndg4hZcvkpBBCeADzgRuAIuAn7DP9adVuHoQQIcDSaldQg3Tu3Flu3ry5RexUUFBQOJWMjAx++eUX7r33XgAWLlzIgAEDcHNrerTPxYAQYouUsvOp7S25uHslcFhKmVttwAKgF+AuhNBUz/qDgYwWtEFBQUHhtJgKCyk9dIiSykpmzp/PzI8/xmq1Mnz4cEJCQhgzZsyFNrFZaUnhPwp0F0I4YXf1DAQ2A38B12KP7LkNWNSCNigoKCg0SumBA6y+917mrF3LIqORSim56dprmf7aa4SEhFxo81qElvTxbxBCzAO2AhZgGzAL+BWYK4SYXt32WUvZoKCgoNAY0mbj0FdfkbNlCz9XVZGg1XKjoyM3jBtHcETEhTavxWjROH4p5XPAc6c0HwK6tuR5FRQUFBrDZrPx/fffs+Cnn7g7PR2DSsW7bm54VIdSFmzZQvDo0RfYypZDSdmgoKBw2SClZMmSJXTo0IFx48ZxKDUVdadOADWiD+Detu2FMvG8oAi/gsJxqtKhbBeY8i60JQotQHp6On369GHkyJFUVFTw/fffs2XrVro99FCtGPnAYcPw7tHjAlra8ii5ehQUpBXyl8KBiWDOBcd4iP8cDB0vtGUKzUBpaSkGgwFfX18APv74Y+688060Wi0Abq1a0W/JEkoPHkSt02GIjUXn7n4BLW55WiyOvzlR4vgVWpTy3bC1O9TsKwQcYqD9StD5XDi7FM6JlJQUnn32WVavXs2+fftadNPVxUpDcfyKq0dBofJwbdEHqDoAJmWLyaVIVlYWkyZNIj4+noULFzJu3Lia4igKdhRXj4KCzrdum8YdNB7n3RSFc2Pfvn106NABs9nMhAkTmDp1KgEBARfarIsOZcavoODUGoIfPalBBTEzwSHsgpmkcBpM+VCZApYSKioqWL16NQCxsbE8/vjjJCcnM3PmTEX0G0CZ8SsoaAwQ+gR4jQDTMXCMAKdWF9qqy5b8CivpRVacdIIITw0a1SnpiIvWwIH7MJcm8/nKcJ7/rJCS0grS9izAQ3OAafdGgavi2mkMRfgVFAA0ruD23w7huxRIPmbmvsWFHMy3oFXBY30NjGvvjIu+2jlRcRDbrtH8+EcJz3xm5WDGfnq1c+Llz1/H48jVgL2KFrpAaLsMnGIv1KVc1CiuHgUFoMpYRkbWTvanrCIndz8W6yVcFM5SCtaKC23FGVNhsvHyqhIO5tsX2s02ePnvUnbnmE8cVJnC3kOl3PSCBUc9LHlFw+p3LfRpY6JG9AFMmVCy7vxewCWEMuNXuOwxmsrYsPkrtu9eUN0iuGrAFOJiBlxQu84YcxEU/A7pb4LKEUKfBPd+oL40whgLKm38c9hIuIeKweFVZJVrWHZQxdEiC7b0LaxevZonJvYnIULFP+9p6ZUoUKkEqBzAVlV3QEvh+b+ISwRF+BUue/ILjpwk+gCSlWvewd+vFW6ul9DiYNGfsO+2E6/3jLG7O9z7XTCTzgQ3B8ELV1ThVvEH2Ud/J9TJl7ZhfXj/0W9Z88dvBAQEMPGecbgGTqIPH57oGPkKdaVMgOHcip7/l1GEX+Gyp7KqqE6byVSO0VR2/o05W2wmyPigbnvuz5eM8DtrIdC0lN37F1BcbOSXBWvZsvFjXFxceOmll5g8eTLOzs5gmAreo8CYBQ7h4NwWrCX2n4yZoPWGiBfA0OlCX9JFiyL8Cpc97q6BqFVarLYTvmRP9zAMzvXE91+0qEDjXbdZ63n+TTlLysvzSEpeXPN67+58+g8KZfoLb9Gj2zUnDtR61r2ZaZwh5BHwGwdCD9pLq1LW+UZZ3FW47PH0CGPEkBdwdrILp7dXNEMGPIWj4yUkHioNBD8IqE9qc7KHqF4CFBcX89JLrzPrg11IKXFz0/Pcyz0ZNTYaH58zSJuh81VEvwm02IxfCBEH/HBSUyTwLPB1dXs4kApcL6VUVmEULhgSQZq1LeUhr2HQVpBldiPX6s2lNN8HwK2nPb9Q4Ur7gqd7fzC0v9BWNUplZSUzZ87k5ZdfpqCggGFDe+BjaEWMZxuqZAXZVSn4uhoutJn/ORpM0iaEaEpqQrOUctdpTyKEGntt3W7AfUCBlPIVIcSTgIeU8onG+itJ2hRakn25ZoZ/mYv5pGjAaC8NP93shaeTuuGOlwjbM00sP1BFqdFKrzAHeobpcHW48Ne1Y8cOhg8fTkZGBkOGDGHGcw8QrTnCrvfWkPXb7+i9vWn74uME9/ZAHTj2Qpt7SXI2xdZXAZsA0cgxEdhn7qdjIJAipTwihBgN9Ktu/wr4G2hU+BUUWpKjRZZaog9wMN9CVqn1khf+5GMmqoq2MMF/EzpbJemiDbsyutIr6sJkHbXZbGRmZhIQ4I+PrwsdOiTyzTff0L9/f6zl+Wx5eC5Zvy4DwJiby6Z7/4fLku/xCrwg5v5naUz4N0kpGw1kFkKsbOJ5bgS+r/7dT0qZVf17NuDXxDEUFFoET6e6S12ueoGr/tJfAvOQh/DK2MOuT/6gKqeQ6FtNRPUSFFVehbvj+bs+KSUrVqxgypQp5Bfk8clnD7Pv4G8MHaPFxSOPqqpSrAUVpP38e52+panH8Op53ky9LGjwf/50ot/UY4QQOmAU8FM9/SVQr69JCDFBCLFZCLE5Nzf3dKdRUDhrYr213NnZuea1AKYPdiPE/RIOerNWQdEqNAe28s8tU8hZuYbiPXvY8sTblK/aik4WnfGQpUYrB/LMZJY0kAdHSqhKA2PtdNYbNmxg4MCBDBkyhLzcLO64/Ur2JC/CZrNgNleydsPHHM7YicbZuVYlrOPoPJWaCM3NaW/5QoheQgjn6t/HCSHeEkKcSdrCocBWKWVO9escIURA9XgBwLH6OkkpZ0kpO0spO5/Rqr6Cwhli0Kt4uJcLP93sxczRHiy5zZuhsZfGbtcGKf4Hdo+hYNtepKV2rYF9H/+Euqz8jIbbn2vmjp8KuPKzXIZ/mcuvyZWYrCfN2Uw5cHQGbG4Hmzva4+nNhaxevZru3buze/du3nv9KXb/PICoVtU7bk/iSNpm9J6edHj5ZYTmxA3X94orzqj+bVVVKWZz5Rld2+VIU6Y0HwHthBDtgEeB2dgjc65o4jlu4oSbB2AxcBvwSvW/i5psrYJCC+HqoKZryKXtz6/BZoK0t0Ba0Oi1dd7WGFwQetcmD1dusjHtz2I2Zdj3ORRU2rhvUSGLb/Omrb/OflD+b3DkRQCOZEt2rX6IEXeH0KvXCGbOnMn48eMxlM3DlrkSb/eB5BUcqnUOB+cQwC70A/78k9L9+9G5u+PWpg2Ofqf3BpdXFHLg0N9s37UAJ0cPune+jaCAdqjVl/BTWwvSFCefpdolMxr4QEo5E2hSfFX1k8Ig4OT98K8Ag4QQB4Arq18rKCg0F9IGQgvShGdrLTqP2gVlEp56Ep1b02Pdc0qt/HukdtI6CRwuqH6SsJkh+0tyiyQPf2AhdpyJu1+3YMpcgEqlYtKkSRgMBrAZUZVvp314BDrdCdeawRCIj589iFCoVHgkJhJ6zTX4DxzYJNEH2H9wJav+/YDikkyycvbw829Pcixvf5Ov8XKjKbfDUiHEU8A4oK8QQgXUnUbUg5SyHPA6pS0fe5SPgoJCM2KzWcnJ3UdaxlY04lqC/a/Co/IPOs35gJx/NmIuKMS7Xx8MXc+siLxBr8LfoMLHSXJVpJHDRRoW71fhVb0oXlJawVufV/DmZyYqjHDnUBXP3q5B59669kCuXUBo8Ts2hRt6Pka+yQWrNpgiTQcCvc8+bKeysphtu+bVapPSRlbOHgL8WjfQ6/KmKcJ/A3AzcJeUMlsIEQq83rJmKSgonClZ2XuYv+RRpLTHpmo0Doy88kkWrp+G3s+ALtiR0qMf0t23nG5etzZ5XB8XNd8Pq8C0axMV6/Zyha8PEwd2xcvXvsVtT1ISz3+wmWv76XjxTkl8mAq0fnV3Dbt0gLa/Q/o7eJQvwuB/L6VO3Ylz80GIxqLGG0el0uCgd6W0rPZyoU7r3EAPhUaFv3rj1fdSyv7H26SUR7H7+BUUFC4SbDYrW3fOqxF9AIuliiOZe3B28qK8Ih+jsRSAXUm/kthqJE5OTaspbLEYqVr+K0lPPgeAVUo2+vmgGjSQ19//iB49erB//35iAq1QvguEGlzagmN07YGEyr672NAFIS1o1Y40RyYhvd6ZHl3vZPHSpzkeJOjo4E6Af0IzjP7fpFHhl1JahRA2IYSblLL4fBmloKBwZkhpo9JY9ytqNJWj1TjUanM1+KHR6ps8dunRwxx4+S2klGw0m/m+spKMoiI6CRUWiwWNRkNMTIz9YOf40w+o0tJEb3GTCQnqyLWj3iYjaycOeleCAtvi5aHUTG6Iprh6yoBdQogVQE0MmJTywRazSkFB4YxQq7W0bzOGrOzdtdqjI/qQkbmj5rVKpaZHlzvQaZ2aPrjZxsHCQj4uK+Og1UqQSsVjzs7cO+0ZNJqLI2pGo9bi5JaAo4hDpxK4GC4Ouy5WmvLpLKB2VI6CgsJFSGhwZ4YMmMKW7XNRa3R06zieoIB2XD38VY7lHcBsrsLLMxxf7+jTD1aNyWTCNTyKyGFDKfvhRyY5OdFXp0Pr6Ih764vHlXIwz8zDvxayM9uCSsDdXZyZ2M3lkk+50VI0mKTtYkJJ0qZwuWKxmsjK3sPOpMVIm422CaMI9E9Eo9E12MdoKkclVGi1Z78JLTk5mWeeeYaysjKWLl1KyaEUDsyeTeaChbjERNP6qSfx69n7rMa2mUzkb9lC9h9/oHF2xn/gQDzatTtrW81WydTlRczdWXvj1ufXeDIw2qGBXpcHZ5Ok7XjHGOBloDVQ8ylKKevurVZQUGhWsnKSWLDkfxxftExJXcOY4a8RGtxwdSm97uyjWdLS0nj++ef54osvcHJy4tFHH8VqteIaGUXH6S+R8NDDaFxc0DidgavoFHLXr2f1mDH2FA9A8ltv0e/XX89I/I8WWSiosOFvUKNRwcoUY51jknPNl73wN0RTXD1fAM8BbwP9gTtQCrgonAs2M5gy7QXBdZdc1vvzyt59v3NqOqudexY3Kvxny++//87o0aORUvLggw8yZcqUWkVQhEqFg++5/X9ZTSb2vfdejegDWCsqyF6+vEnCb7VJlh+o4n+/FVFqkvg6q5g91oPOwTp+21e74Hqkp+Lnb4imCLijlPJP7G6hI1LKacDwljVL4T9L5SE48CBsSoBtPSFvsT3FwAWmPC2NIz/+yO7p08lctgxjQcGFNgkAlagrXkLVfH7rsrIy9u7dC0CPHj2466672L9/P2+//faZVb5qKjYb5uK60Uem0pImdU/Jt/DA4kJKTfYbx7FyG3cvKGRCVxf8XE7I2fB4BzoENm/k0H+JptwSjdW7dQ8IIe7HXlDFpWXNUvhPYjND+tuQ84X9tTEdkm6A9n+Da7cLZlZVbi6bJk4kb+3amrZWjz9Oq8ceQ3WBo1ZiIvuQtH/ZSfH5goTYwec8rtFo5JNPPmH69On4+vqyc+dOXF1dmTlz5jmP3RiZlRqiJ4xj0/9tOdEoBIFX9mpS/4wSK7HuNm73ysKpLI8itwA+zPBDBfw83ptDBRYcNSqivdW4XQTFZs6WwtJSth8+SFZhAZH+AbQLj8RR33xuq6b8VU8GnIAHgRexu3tuazYLFC4fTNmQ880pjRIq9l1Q4S/Zu7eW6AMkv/02oddcg+F4fPoFws28nWv63cTe9CNIKWkVHIGLZTPQ/azGs1qtfPvttzz77LMcOXKEfv368corr6BStbz3NrXQwrbUdLomuNPhwzdJnT0HjbMD8fcMxTNoJ3DVaccIcjDzQuUvpE54hiqbDSe9ng/f+Rgv5+EEumoIdL303Tsl5WU8893nzPz155q2T+9/jLsGDTunHc4nc9pPSUq5CUAIYZNS3tEsZ1W4PFE5gT4EKg/Ubtdc2OLY1qqqOm3SbMZmuvAuKOfyZbgX/UaQIc6+8zVnL7j1BnkfnIUIzJ8/n9tuu42OHTsya9YsBg0a1HQxsRnt6ZfVLqCtu+e2pDSHgsJUhFDj5RGOi4t3rfdzi3KoPPwO83J3odE4EP5wL7q27oZP2g2gfbZJJngeO8DO56fWrBHYjEaOPPkQUV3bgltE067jImdP2pFaog/w0OwP6JfYjuiA4GY5R1OienoAn2F374RWp2f+PynlpGaxQOHyQecFUW/C7jFAdTEPQxdwObOkYc2NIS4OnZcXpvz8mjb/wYNxCg+/cEZVY/UYjrboN6jcV9NW5T4GhzMQ/VWrVpGbm8u1117L2LFjWbRoESNGjDizWX7FATgyA/J+AocoiH7LXsxd2MfIyz/MoqVPUVZuL5rk5RHB8MHT8HA/IVSyfB+5ufYS3RZLFQfT/8VoLWJExHB0nkOaZEZ55tFaC8MA5uJiynPSMUT8N4S/oJ71jvKqKorLz6yGQmM05X/+HWAIkA8gpdwB9G02CxQuLzwGQIfVEPcZtJ4Hrb4Hh5ALapJLWBh95s8neOxYnEJDiX3wQdrNmIHW+cIn+dJ4DcbsM47jpa8tXteh9mlabMXWrVu56qqr6NevHzNmzEBKiUajYdSoUWcm+tZKODwVcr8HabHfhHaPtuflqSZp37Ia0QfILzxMatqG2sMY69Zcyi/KxBwyDVyaFsqp8/epVagFQO/ri8bbvenXc5ETFRCEo652So3WoeGE+TZfldomOcSklGmnPA42UHtNQeE0CA0YOtp/LiI82raly4cfYikrQ+fh0Wy+1HNF4xQCse8hQyYjkWgco0Dd+A0pJSWFp59+mh9++AFPT0/eeOMNJk2adPbXZMyA/MW126QFKvaDSzusVjOZObvrdMs5llzrdaBvDAYXX1p59EQjNaQZD+DhG4mja2yTTXGKiiD+jec5MPVlLGVl6L29iXtjGi7B/528PHFBIfz67MvcM/NNUrIy6BmfwMx7H8Lb1b3ZztEU4U8TQvQEpBBCi32xd29TBhdCuGOv2NUGezDyncA+4AcgHEgFrpdSFp6p4QqXJmVGG1szTfxz2EiQq5re4XpivC+OsDu1TofasznyRTYzaieESyJNle0jR47wyy+/MGXKU9wz4Wbc3T3Q6xve6duU86P1BvMpta+r12bUai1xUf3rCH14aO0FaB/nUHqU9GDfU69jKSvDf/hVtJp6J6ozCE91cw/EOGwwMtQDU2EhDv5++LTqjIND0yuKXewIIejftgPrXvuAwrIy/Nw9cGvmp8/TpmwQQngD72KvliWA5cDk6oIqp+v7FbBaSjm7uui6EzAFKJBSviKEeBLwkFI+0dg4SsqG/w4/J1VwOCebIIdjWIUja7J8eOwKHyI8Lw7xvxQpKCjg1VdfRa1WM2PGDABSjyRz4PAv7N2/HJVKTZcON5PYehROju5nd5JjP0Hy+BOv3QdA3Oeg9wfsC7vrN3/J3v0rEELQvs1YOrW/AWenEzfSnL//ZvXYsbWGjf6//6Pd9OkI9ZmFXprNVRhNZTg4uKFRK387DXHWKRsAKaW85SxO6IZ9LeD26kFMgEkIMRroV33YV8DfQKPCr/DfILfcipPlEC5p08iqtD/kXRk1mpzim4nw9D5Nb4VTKS8v59133+W1116jpKSEu+66CyklQgjyCreTtG8ZAFarjfWbv8LLM5LoiLPLr4P3KGi/yh56q/WyF1WpFn0AtN7oQiYSZLgalRBovAJRnZIBtGh3XXdQ2oIFxD/0EA5NLLFYczqtA1qtko7hbGmK8K8XQmwHPgeWyaZndYsAcoEvqiOBtmB3E/lJKbOqj8kGmm/FQuGixmapID/lM6qMJfh4x1BVVczRlEVEhHYDFOE/E5YvX86tt95KTk4OI0eO5KWXXiIxMREAs7mSfQf+rNPn8NFNZy/8Kr19r0UD+y3WpJqYuKiME+W4S5hzg4Y+4SfE2SkoqE4/t9at0bgo+0HPN00R/ljsbp47gfeFED8CX0opT1fJWAN0BB6QUm4QQrwLPHnyAVJKKYSo90YihJgATAAIDQ1tgpkKFzuOqnIMOJCgH0nxz+vQhwSi630lwqIs8TRE6rEKNAd2Yko9RGlqKk4dOhB55UDCwsJo3bo1CxYsoGfPnieOL7RQWCExuMeSV3Co1lgOLuEtYqPFKvlqa91Qw1/3VtUSfq8uXfDq0YP8desAUDs7kzBlCpqLIHrqcqMpG7gksAJYIYToD8wBJgkhdgBPSinXNdA1HUiXUh6P6ZqHXfhzhBABUsosIUQAUDfGy37eWcAssPv4z+SiFC5OHPRuhGX6sO+Z52ratHPd6bbguwto1cXL4SN5lH7/OQffe4dNZWXM02gI9fHh+9mzievTh5UrV9Y6PinHxM0/FNAhQMtNccNxzFxHZZU9L46raxhmx7NPfdwYKhV4O9UND/U8pc0pOJgen39O0Z49WCsqMMTG4hrb9IiexrBUVJC/cSMZv/2Go68vAVddhXubNs0y9n+Rpmzg8gLGAeOBHOABYDHQHvgJu0unDtWF2dOEEHFSyn3AQCCp+uc24JXqfxed+2UoXApYC4pJfffjWm3moiKMB9Og/dmlIDhrzMVQeRBUGnCMsu9Gvcio3LqRRS+9yLeVleyzWvFTqRjUuzdHFyzAt0+fWsfapGTO9goKK22sTzPRPjAYfeRrhGszkELNvtIgYp0CW8ROlRDc3tmZZQeqsFSnFHLUCq6Kq+uDd/Dzw/8M/flNIfvPP1l/24lMMvvef5/+S5fi1rp1s5/rv0BTXD3rgG+Aq6WU6Se1bxZCfNxAn+M8AHxbHdFziBMpnX8UQtwFHAGuP3OzFS5JpERa69kCcr6f5ypS4OD9UPSX/bXvLRDxAujr+qAvJHM++YCXy8rwEIIJTk701+nwyM3FqV+/OseaLJKd2WYAKsz231v5+PLWTgNOWsEjvQ10Cmq56JdOQTrm3+LNuqNGNGpBz1AdCX6NFIsxm0k9lo1GpSbCz/+ccgWZiovZ8/LLtdospaXkb9yoCH8DNEX44xpa0JVSvtpYRynldqBOKBH22b/CZYZjQADxjz3GzqefrmnTGAzn/5E894cTog9w7FvwGAh+N59fO+ohJSWFkpISOnTowMihV5H2z99cpdejr9585RIZiW/vugu0DloVYxMc2VUt/n8cNLIp3cRHoz1oG6DFoG/ZTJUqIWgfqKN94On3CxzNzWH6D9/w2R9L0Wk0PH39OCZeNQov17PL2SStVmz15FuyGusWZ1Gw0+BtVggxDWp8/I0eo6DQVMKuu44un3yCd+/eRN5xB1csWoRbq1bnzwBrFeT/Urf95BvBBSArK4tJkyYRHx/P5MmTAQi5ohO3tGtbI/paNzci774D98T6b5RXxTpwW0cnNCrwcRZ83rmAoM2/UrJiKWVHjjSrvZVmG6tTq3h8aSGvriphV3bTE9r9sPovPl3+KzabjSqTiWfmfM7qpF2n79gAek9P4h58sFab0Grx6tLlrMf8r9PYjP9uIURj1REEcCMwrVktUvhPo/f2Juy66wi95hpEc6QCNueDpcReyes0qQwAe1ii+wAo21a73fU8rzFUU1hYyGuvvca7776L2WxmwoQJTJ06FYCU8h04Tb0W/yI9GrTIQAN/FSzgGlWPer+4ga4anhngxh2dnHFI2cmGa8eSWloKgGNwMH3mzWu2xdTVh42sP3yMK4PKMdnUvLfGhcl9fGjTiHsHoKyykm/+XlGnffm2TVzd/SxDTYGgUaNQOzmRMns2Dv7+xE6ahEf79mc93n+dxoT/U04E5TZ2jILCGdMsol+0Gg5MhsokcL8Sol4B59O4jYQA//FQ8BtUVGcecbvC7uq5AHz//fe88sor3HzzzbzwwgtERUXVvOfrE8Py3a+cOHg/tE0Y3ejGJa1aEOYq2Dx7FpZq0QeoTE8nZ+XKZhH+MqMNl6ojXFv0L3lfrcDg5cb/Rg3lWEk38AtotK+DTkfn6Dh2pdYONU0IPbfMmnpPT8Kuv57gUaMQGs0FL6BzsdPgpyOlfP58GqKgcEZUVGeItFXYXxf9AfvuhsTf6s0VXwuneGi71D6G0IBjHOjOzwYys9nM559/jru7OzfccAN33303vXr1ol099WaDA9oSHtKtJsulu1swbVuPOm1uG5vJRPHeuum0SlNSmuUapLTiuXs9u+//X01b1rzFdPzpW4hpXPg1ajUPDB/Dog3/1qQfbhUSyuAO9S0FnjlqB2U3b1NQbosKlyaVB0+I/nHKtoMx7fTCD6Dzt/+cJ2w2Gz/++CPPPPMMBw8e5Prrr+eGG25Ap9PVK/oABj0MaRNBQYQ/VinxcFDj4mA57bk0Tk6EjxvHjidr7ZfEf2DzPNWoKgo5+tHsWm02o5HSjZug94DT9u8QFcOGNz5k95HDaDUa2oZFEuJzbkXcFc4MRfgVLk00HnXbVE6galo8fklJNvmFR1Cp1Hh5huPi3HIz/tWrV/Pggw+yfft2EhMTWbJkCcOGDTt9x+J/cUh9gFrR9+IYhD9/2upbQSNGUJmVxcGPP0at15Pw9NN4d2+edQyNRl9vplC1TU3SW28hbDb8rrwSj3btGkwFHR0QRHTAxRU+ezmhCL/CpYlTK8zBz5FvdqHCZMNVZ8PT3R+VY+Rpu+blH+Ln356ivCIPAG/PSIYPnoa7W/MKkc1mQ6VSUVRURHFxMXPmzOGmm25qesx68b912/IWQfCjoG089NEpMJDEqVOJuuMOhFpdb56cMyG7sIBlWzfyzV8r6BHfmjsnP8DWCRNr3lfp9ai0WnY/b/cQ733rLfotWYJnx4ur7oKCnabs3I0FPsKeXK2NEKItMEpKOb3FrVP472PKgYpkkDZwbnVa94vFYqKwOI2qqhIqbd35Y80bmC1VqFRqhl45lejTzISllOxO/rVG9AHyCg5xJG1Tswn/rl27ePrpp2nbti3Tp09nxIgRDBkyBJ3uDHPiu3So2+beFzRNe6oRajXOzZDnymq18sGvC3npxzkArNy5lZQuPXj1q09Im7sAvY8X3l17sPvFE5Jgq6oiY/FiRfgvUpoy9fgUeAowA0gpd2IP47z0sBRD8VrInQ+l2+zFoxUuHBUHYdcY2DkEdg2FncPsC64NYDSWsXn793w//14WLHmMVes+pEfXu1Cp1NhsVv5Y9RYlpdmNntJqNZGVvadOe07u6XIOnp7Dhw9z66230q5dO1atWoWXlxdgL6xxxqIP4N7HXtf2OPowCPw/EC27GetUjuYe442FP9Rqm7tpHVud8uj1WBodJ1Ry5IfvqMrJqXWMqbi4zlgZ+Xn8tnk9P635mz1HD7eo3QoN0xRXj5OUcuMpvrrTrzBdbFhK4ejrkP5GdYOAuNnge/Np/aUKLUTBb1C+9cTriiTI/QnCptZ7+LG8g2zY8vWJwysKSN6/nMiwnhw8vBqjsZSqqhJcDQ0/NWg0ekLD+3Es70Ctdr+Ac9vs8+mnn3LfffehVqv53//+xxNPPIHnuVbzcgiD+G/sYafSDE5xFyathKBeX72UNijdjKpsJ9G3PEvu6tr5GkPGjKn1OjUnmxtef56N++2Vuhx1ela88Aa9WivJ1M43TZnx5wkhoqjOqCKEuBbIarzLRUjF3pNEH0DCwclQdajBLgotTPEazA5DKTA/S4HxGUyOo6FwJTSwWby0rG4i12N5B/Bwt7szXA0BODv7NHpKk1Vy1NaN4NCBgEClUhMRcy2HjHFnbn5xMTnVs9yuXbtyxx13cPDgQV599dVzF/3j6LztM3+PAWcs+maLhQOZ6RzOyaLpZTTqEubjxxPX3FSrLcDTi7bhsaAPBad4fHt3pvsXn+PRqRPePXvS+8cf8Txl5+z6fUk1og9QaTIy9dvPKK+qPGvbFM6Opsz478OeHjleCJEBHMaerfPSwpxXt81aBpai827Khab08GHydu5GWix4tmmNe9yZi15zUK66nR1vzSJzuf2G7Nu3Ox2mP4ShgSew+kTd1zuGwqKjGFz8uWrgFJyd6on2OQkBLD1kwGS5nasSx2C1qfj4gCt9hFOj/U6msrKSmTNn8vLLLzN06FDmzJlDu3bt+OSTT5o8RkuTlnuMV+Z/x6zfl6DTaHnmhvHcM3gEXq5nXptWpVIxaehoYgKDmfvPSjpERnN97/5EhUVA2FpAi1brRvBo8B80GKFS1RtPn56fW6ctKe0IpZUVODs41rSVV1WSnpeLk4MDId5KmGdL0JR8/IeAK4UQzoBKSll6uj4XJQ4RIPQgT/LrO0SCPuTC2XQByEtKZv0N11GVkQHYk6R1n/8z/p3rWUhsYTL+TSdz+T81r4/9s570PwbRqoEnf4Ne0DVxMBt3rQAkTo4e9Ol6M+XSi6g2PhjcTx+SqVUL7u7qwvgfC1iXbo+M0agkLw1xPE1PsFgsfPnll0ybNo2MjAyGDBnCI4880qRrPd/8+O/ffPibPeO5xWrlqa8/JcDTi87RsWe+S7byIL62FG5p784tvZ4CzUkb+rW1P3ONU8M30I6RMXXaxvW7El+3EzfrfelpPPblRyzZuA4vgyvvTniAa3r0xeFs1kgUGuS0rh4hxCNCiEeA/wPuqX59lxCifYtb15w4xUPCPNBVPy47tYFW39hzvFxGZC5fXiP6YE9fe/iLz5E223m3JXvl33Xbfl/eoC3ayj100f/Czf2GcU2fEdzQJYbgI2MxWrUM+cbCT7sqsDXBpdEtRM+fI838GLWLeXF7WXE1tAs4fcriF198kXvuuYeQkBD++usvli1bRseLMGqlwljF1yuX12lfumUDN772Imv31q192yAlG2BrL/su6R1XwMGHwFRv7aTT0jU2nln3PYqbszMqlYpx/a7kvmFX14S3Gk0mXvjha5ZstK8V5JeWMO7Nl9h26EBjwyqcBU1x9XSu/jme0nAEsBO4Vwjxk5TytZYyrlkRAjwHQYc1YCkEnZ+9aPRlRuXButErFXv3YKwy4uDUhFlveTmFu3ZRmnoEpwB/PNq2Re/RuHulITx7dufYsmW12jz69W1wsV1oPdGUb8anfHNNm00XRkqxHoBXVpXSL1JPhGfjIl51IJl9N1xPZfUNsLhtW3y/+AKXiNozYSklK1aswMvLi06dOjFx4kQ6duzIqFGjGtyYdDGg12jpFBXDztTaKRqCvXxYumUjz3z7BYunvoTz6dIbWEog5XGwnhSdc+xb8L0ePIecsV0ujk7cM2QEQzp2wWg2E+rti/6kmXx2USE/rqmbJXVfRho94hPO+HwKDdOUxd1goKOU8lEp5aNAJ8AX6Avc3oK2tQz6AHBufVmKPoBnv7q7N33HDkOjP/2MV1qtHPxmDquGDWPrpImsGTOGXS+/irms7KxsEW0Dce9+YgHQtX0imu72DVjG/HzyNmwgf+NGTAUFAKhd22F0HQKOsdgMV4DajeKQd5j2t108qiySCnPjM34pJYe//rpG9AGKdu4k+5Qyhhs2bGDgwIEMGTKE119/HQB/f39Gjx59UYs+gFqt5v4RY/BwOeGSifQPRK/VUlJRTlJaKmVVFY2MUI2lqG4WUwBj5jnZF+rjR0xgcC3RB3BxdCQmMLjO8d5nmadfoWGaMuP3BU4OeDdj38xVKYRoNBBeCJEKlAJWwCKl7CyE8AR+AMKBVOB6KaVSbfs8UREMkVMe4uj7n2EzmQi+82aqWrsjrFVQlQ7SAo6RoK7rqy04kELS89NqtaXOnkXA2GsI6n7m4ZB5mlwsE3sQeddokDYKHIrItB0h/PBhNt1/P/nr1wPg3bs3nd97D5fwcAp4hj1z51C8dz9BN8xgaUEMRUa72LcP0BLs1viftM1oJK963JMp3GYXuKSkJKZOncrChQvx8fHhvffeY8KECWd8bReajlGxrH/9Q9bvS+JQTiZF5WW8Ov97AMb1G4SPq/vpB9H6gMcQKDilfkETdkefDV4GV96f8ABDn38Ss8UeMT6kY1c6RTVPKmmFEzRF+L8FNgghjtfGHQl8V73Ym9SE/v2llCeH1DwJ/CmlfEUI8WT16yfOxGiFs8fLz53f/ZOI//gRVKhJLtpEzwBH1DlfwOEnASv4XAsR08EhvFbfisJCbPVUNTIW1I3WaArxMQOYt/cR9lns1ZPUai3XjX6XjO+X1Ig+QN6aNWQtW4bfwIGsHnMd5uqNQfnr19PloceJDridGG8dj/Yx4ObQ+EOs2sGBkDFjKNqxo1a73wB7crElS5bwxx9/8MILL/DQQw9hMJwuM/nFS2xQMAEennz/z5+8s3geEhjffzATh45qWtoItSNEvEC5sSvlGVVoXHS4RgahqW9HcTPRP7EDm9/6mOT0NNydnWkXEYWfezOFxirUIJoS3yuE6Az0qn75r5Ryc2PHn9QvFeh8svALIfYB/aSUWUKIAOBvKWWj8YSdO3eWmzc36ZQKp8FSvJmMrD1s3r8Fs8VMp9hOhPr4od9/k322f5zw6RD6WK2+mSlJ7LjxVspTTux90BgMdFr0IyHtu52VPcfyDpKVvRubtBHo3wZf7xhWjRpF3r+189T4XnklETfdxIa77qrVrnZyoteqtbiFBaPXNC0HTvmRI+x6/nnSf/6ZEiFYGRnJkLvvZvz//R+VlZWUl5fj7X1+0jSfL47m5tTrVz8dhdu3s+bGGzEesy/oRk+YQKvHH0ffXPsUFFoUIcQWKWWdnNdNStJWLfSbq2f5Y4UQz0sphzelK7BcCCGBT6SUs7C7iY5vAMsG/BoweAIwASC0GfKNnC+qcnMx5uWh9/HB4SIUD41re8IwEqTPRtpMaB0sUP5vbdEHyJsHQffZZ33VaN0dCZ4+mWPvzqFw/SZc4uMIfWoi+DUtd0x9+HpH4+sdXavNoVdPOEX49d271Vu8RajVaEVlk0UfwDksjLiXX2aRoyMffPklFVu3EpVp91s7Ojri6Hj6Re5LjVCfer9mjWIuK2PHs8/WiD7AwVmz8B88GP8Bp0+/rHDx0pRwTp0QYowQ4ifsO3YHAB83cfzeUsqOwFDgPiFE35PfrK7n21Ah91lSys5Sys4+Po3vxrxYyFu3jr+uuooVvXrx11VX1etLvuAIDbj1QhPxBFqPbpA6rf7jXHuBqnbUh4dbMIWulVRO6EjodzMQU4ZyUJuCt0fdBblz4UhUBIaTUggb+vTmaHgIzq1i0fvWDr8Nf2ACVbYKpNXa5PHnzJlDbJs2vPLhh1w1bBh79uzh+eeVukOnYiosJH/DhjrtJy+MH0dKydaU/cz89WdmLfulToUthYuLBmf8QojBwE3AYOAv4Gugi5TyjqYOLqXMqP73mBBiIdAVyBFCBJzk6jm7oOCLjPIjR1g7bhymQvs6dfmhQ6wdN46BK1c2S4bEplJZVYLJVI6jozs67Wlmrlqv6lA9FTglQEV18jJ9CATcWSesUqPR0rn99WTnJJGTd4AI9z4E+Cfg6NC8URf6kGC+69mRK0cPRwjBz9np3BgcjDbQi5B3nsC8YR9Vh9IIGDyEyqNpbB1zK0EjRhB1110YoqPrHdNisWC1WtHr9Wi1Wjp06MCMGTPo3Ll5Kj/9F9F7euLTty/H/vyzVrtTPX/P6/cl0W/Kw5gsZgAMjk6smvEOHaLqbtpSuPA05upZBqzGPms/DCCEeLepA5+807f698HAC8Bi4Dbglep/FzU8SvNSVp5PzrFkikuz8PIIx883Dgd98yzelael1Yj+cUwFBZSnpZ298JfvgcI/7amLPa60FwRXNyzm6Zk7+GvNexQUphIW0pU+3f8PL8/whsd3ToTWP8DBh+1x2QF3gy4Io60N5fuL0TjvxyUyslb9UidHdyLDexIZ3vPsrukkLFYr6/cl8dWfv2OxWbltwBB6xremf2J7KoxVvPjDN6hUKp65YTxXRDninPUqVe4h7IpMoeuAa9g74XmM1blyDn7yCSUHDtDjyy/RupxwPUkpWbhwIVOnTmXcuHFMmTKlpvqVQuNonJ1p+9xzrD14kIojR0ClIv6RR/A4pWKYzWbj/SULa0QfoLSyggXrVyvCf5HSmPB3xJ5++Q8hxCFgLnAm+WD9gIXVMc8a4Dsp5TIhxCbgRyHEXcAR4PqzsvwMqawq4e8175GSuqamrUfnO+jc4abT1jBtCjpPT4RaXcvlINTqs97cREWyPV3x8RxD6W/aRdp7dL2HFxQeZdHSp7BY7FE3R9I2UllVzJjhr+Ggb8AHr9KB9yhw7QbWKtD5U7z/EBvvvYPiXbtQ6XQkPPUUEXfcge4scrycjnXJe+j39MPYqnfqfrXyd/548Q3aRrWjRN+NQYPjAUGPiFwMe/uBrYo2XvegjeuJSCmsEf3jHFu5kvLUVNzb2HM+rFy5kqeeeoqNGzcSHx9PYmIiWCsQ1nLQXRruwwuNe5s29F+2jPLUVDQGA4aoKNR6fa1jrDZbvXl4MvLOLtpLoeVp0McvpdwupXxSShkFPAe0B7RCiKXVC6+NIqU8JKVsV/2TIKV8qbo9X0o5UEoZI6W8UkpZ0FwX0xgFham1RB9gw9ZvKCqu6688GwzhviQ+/UCttsSpD2IIP/NFNcC+Vf7kxHIqJ4r2ZbLn5Rmsu/120hctwlhw4qMrKthXI/rHOZa7j9Ky2uJYL5ZSyJuHJfkpdk9/nuJduwB70e5dzz9fJ/SxufjmrxU1og/22fmHvy1iV7aZN9eUUmZWUWZS2TcR2ewhny75n9LJ9i6+XnWffIRWi6palJ588kkGDhxIVlYWn332Gbt27WJkXy/Ycy1s6wWp06HqaItc138NrY8v5ladkBHxdUQfQKvRMHHoqDrt1/Xudx6saz5yi4soKrs0U5GdKU2N6lkLrBVCTAauxP4kMKslDWtuzOaqOm02mwWLtXmKsagth4js/gdePzxFxbFKnPwccXP5BZXlKuAsdglba6eqLRWT+ef2tzDl5wOQsXgxbV96idiJE6FiPzpz3XwmWo0D2tP5+U25sO8uKN2AyeERclauqnNI2eHD+Pbpc+bXcBqstroLslarjUN5Jfw2ciMR5e8jUSFd74SKjlBWnbtfH4Z7hCO+fftw7J/VNX0dbr2VUgcHXIGxY8fi5+fHxIkTcXBwgLLdsGtYzQ2Eo9PBkg+Rr4Hq9LuWLxWslZWUpaYi1GpcIiJQac/t2o4WWfhySxmb081oVXB/TwO9w/Vo1bXXf4Z06MKn9z/GjJ++Ra/V8vzNt9OndeI5nft8cayoiB/WrOSNhT/i5uzM9FvuZHCHLv/pxHBnVHNXSmkDllf/XFK4GIJw0LtSZSypafPzbY2rS+Ol/pqMtKExbcRLvxGv4wk/Tfb2s8K1iz0CpzrMsuiIrkb0j5P0yisEjx6Nk34/3qZVxEV0Zt/hE/sd+nS/GzdDQOPnqdwPpfbIDa3Yh0fbVuRv2l7rEMeAk8YwF0DpJvv6g2MkGLrZ02CcBbcOGMwXfy6rlSt+0vDRtPPcjG/qSUXeDq62C3TZNgj5H1TsxSHrDjo9NZn8m8dwcPt+Ptu+nR/efpv/q6hg5syZdO3ala5du54Yo2LPCdE/TtZsCJ5cZ6PapUr50aPsnj6dtPnzESoV0f/3f8Q98AAOfmf31GmxSZbtzcehMhld9r/4ewWw/XA3vJyiaRtQWxQ9Da7cPXg4Y7r3RiUEHobmdw22FD+vX8ODs963v8iF0S9N5e8Z73BFm3aNd7yEuWyKradX+hDc7gVci9ejN0sqHfVk6XqQXelI+GlyVTUJxxhw7QElJ1Uhcu0FTvVHmZwWlw6QuBTSXgdjBmjr3qCkzWYvWqJ2waFkBX3bfEqr8PaUV5Xh7uKFj1dME/LKnHhfW7mMtk++yZo7D9fsjvW+5nosMdUzN5sJ0t+DtFdOdPcaA7EzQXvmG3p6xCWwcvpbfPjbIixWKxOHjaJP67bo9t+CdIii1O16QGIo+hFKNiLiv0UW/IrIt6cQqMx5i/cW2PjgZxVSwgMPPMCUKVPqP5mqnicftSuI/86sLn3RItLmzQPseZUOfPghnh07EjJ27FmNl1duJSNtHS/MORHT4WX4mYSn3qZtQP1pG7wusbw6JRXlvPvL/Drtf2zfogj/f4ESow2HoyVUfvATubt24t5/IPpx7SgNaaZ0xDovcoM+BoeF+FT9Tq7DEPAei8/ZJoMTKnvlJUNXkGbcnJLQurnVCDJAqwfvxtHfB6w6iJiB05GHCDNXL6gVARUjweVLUDs3fB7HWHu0UMl6wIYMcyb406mInDKEo56qAF+O2jSEAFQetN+ITiZ/IVQ8AG5nHuWj02rpl9iefonta7VXOPdnV2k3Nq9ZCQi6JNxFG2cDKuceOOy/p+a4p2dbmP2rjZvH9mT6G98RFhbW8Mlc2ttTc1ecqABF5MugDzxjuy8E1uqgAbW6/kAEa1UVaQsW1GnPXrnyrIW/orKQj5d8Xastv7SEo9mHoE3L5Os532g1GgI8PUlKS63V7uN2ad3AzpTG4vgbncKdr0XZc8FS/WXRqNVEV6az8cGbsZTaF28Kfl+KITOT4B/ncVY++FNIL7Yw4ns3tKrbifa6i4P5Viw2WHKbhaDTJA5rFLUe0OPqfYS+Xz7IkcV7KTqQQfiYHvi3OoAwZ9lrs+oDwHxKFEXBL1B5GFwaqWmq84bYzyD/F0zGAv7Zs520zO3298qAXLhiQAjgDbZK7Pn2TsHahEyPp1BYVoqDToejru5i4RFjMOt3nphlrtvxG26ej+JZJfnoZ2f6tVHToU0gz9y6j8nXWgka8BhuwY2IPoBDKCTMh+J1YEy3RzIZujbe5yLAaDazes9O3v1lPlablQdHXMMVbdriqK/9mKrS6/Hu1q3OQvypoZdngl4t0KrVPHvdNUR6O2GRauas3oxGnP/aDS2Fo07P09eP4+9d27FWBxp4GVwZ0Pbiq7PQnDSmSFuw76oVQChQWP27O3AUOMMyPucPk9nMP3t28s7ieVSZTTw06lraF5TUiP5xSnftgOx08GsO4bdSWGn/w8mrsNEzUOCklmSWWM9a+I35+RQnJ2OtqMAQ7IOHZhIeN/ojtf6IsrcAtxO7a9UnV0XyoUzeRkW+Hn1qJYZ4S61Y/Do4RYHTQ1SVHSPzr1vrvC2N1XvsHCLBpTOUbQa1C+j87aLv1PRY7Yz8XH5Y/RcfL/uFSL8Apt4wnl6t2tRySe09uLpWH5tN8tGsj/l53k6OHEmj38j+hOf3YHT7e+kRexizUxO/pI5R9p9LiC0p+1mzdxcdo2JIy8tl+ItP8duzrzCkY+1sqEIIIm69lfQlS6iqTj/hlpiI38CBZ33uIC9vFj3+ENs2v0vRUfv60oP9BtEm7tL6DE9H39Zt+ffV91m3LwlnvSO9WiXQOjT8QpvVojSoBlLKCAAhxKfAQinlb9WvhwJXnxfrzpJ1+/Yw+Ln/1Swa/rljK3sfeKrOcSq9Ho1zI26QM8DNQYVKQKSr5BnvvVjnfIitpBhP9URMhivRubuf0XgVmZlsfeQRspfb19G17u70+eo1PCsnISqrI3iiZtkLygA4J9gXWiv3k53/BOvvex1LaSlC+wEd33iD0BtuQN1IlEJpWS6lJdkE+CWQnrW91nuebtVpErQeEDfbvo5hzoPKQ+BeN7VDQ9hsNj5Z9gsv/vANAAcy0/lr13bWvz6z1kYfH69o0jLsETxJu/JYvDCF7Mxy2rRNxHNwH/6WNli1li9XwZt33ssjbU8z279EycjP5dV537F441oAWoeE89yNt/HR0kV1hB/ArXVr+i9dSklyMiqNBtdWrXD0P/vgBZO5gsMHF1BecSKo4GjqCjq2vvKsx7wYUavVdItrTbe41hfalPNGU6ai3aWUNY5VKeVSIcRFXXXr+39W1ooUAZifdohB11/P0R9/rGlrM3VqnapLZ0ukp4YXrnShU/4Wkm68wb7oCmy7dwLaTz4h9LrranewWUDV8MdfsGVLjegDmIuK2PPer/R4bzFq236Ec3swnDTT1flhjv2SyrQDbLztvpqnG2k2s+XhhzHExOB9Uv6bkyksSmfJ8mcxmSro3vl2CovTKa+w7yFoH98fX7cTERpm4UZ2mZ60HAvODvEEixy8yt+F8Gn2DWGNkFmQxzuLay+kmSxmdh05VEv442OvZM++pRiNpRxJLcFmFXz68TsMaNWKyuJiStzduOXrLzmck8Xzc7/h+t79CW5iUW6r1YLVakKna3px9QvFP3t21og+QFJaKgezMgjzbVjMnUNCcA5pnjrSVVWlZGXXLdNYVn7xZVkxmc1sOpDMP3t24mlwpW9CW1qF/DcnBM1BU4Q/UwgxFZhT/foW4NxK8LQw9fmNN2dl8PC0afiMGklVVjYukRH4d+mKaGCx7EzRqSU9PHZRuPC3GtE/zv6ZMwkcPtxeiLryCBVFOxDlu3HUacH76npdJRVH624uKty+gxWHY/gurQNjEpzoG6nDu1q/Kkw2Pt3uTWBmFuKUsE9sNvI2bsStdWu0p+zAlVYrZYcP0ddtDPn/rqNs+2/0H3QDllAn3DXpeJbORVsVA672zNmpGTv4beVnNf2dHD24tltXPCpTwLlVo5+RVqPF1cmJ0sraawIOJ/1/bd26lSlTpjBu/LUM7BPDgD4mvHV+HP3wU7Y+ZU+k5hgUxPKP3iXh7VcwWyxYbadPLQ6QnbOXrTt/oqDoKAlxVxEddQUG54t3B+/65LrlLjYdSOarh548L+d30LsSFNiOtPQttdoNLme5KbEZkVJitJhx0NonG3/s2MKIF6fUTPh83T34+6W3FfFvgKbksr0J8AEWAguqf7+pJY06V67v3Q/tST5tIQT3jxjDgr27SPj8fVovmkOH2e+y5nDzFXEuLE7n7zXvIhzq3nQ0BgNCpaKyLJsdycv5/s9v+WFTEvuK3TGlvlFv8Wq31nUfOw2DhvF+ko5/Uk08/GsRP+yoqPlD35dn4a01pWRqPNGfmg5apcJmNFKWmlqr2VJezsHPPqPy351suX0Chz74mLTP57Dt5nvgQD6elUvRlq2p2UtQWVXC1t1L6BN8I1eoh3OFYQwxnl3INrrXTetcD37uHrx++7212kJ8fOkYGc2BAwe48cYb6dSpE5s2bULa9MREXUFi60FUJaeQ9uW3NZ9lyNixWHfs5d/bJjLrzomEeJ9evPMLUlmw5DEOHFpFfsFh/ln3ETt2LcRWzyayi4Xu8XX/Bq5s14n24efHx67TOdK72wQMLsefpgSd29+Er/eJiliFZaWUVVbWP0ALkXQ0lf998TG9Hn+Al3/6luT0I0yd81mtp/xjRYX8eyZF5S8zTjvjr47emSyEcJZSlp8Hm86ZbrGtWDXjXeavXYXRbOK63v3xdDEw5Ln/1azc55eWcPMb09n81scEeZ37rK+iooDyinw07SNQOztjLa/+qIQg/uGHUTs4kLp3I39vnFvTZ9nauYzucyPhlftBV9tV4dGhA22eeYak117DZjTi1bMXWUPvYs9Jk8AP1pUxJsGRQFcNuWV2AfvsqBuL33iT7Q/cj6W0FJVOR9zkyaT9/DPBI2uXUCjatYv977+P34ABtStrScmxOQuJeVqCyhmc7ZEhNquFdpoeJN31KDaTCQDXDm1xeuPpJpfjG9WtF39Of5O/dm0nxMuHK9q0Y86s2UyfPh29Xs/UqVN57LHHcDspnK4iLd3+i0pFwpQp7H3ttZqEeDE33oipS08cTpO6O6/gMGZL7Q1c23cvoG3CKFwNdtdJUVkZxRVl+Lp54FhPaoLzTd+EttzYpz9zV9sLkLePjGbSsNHozuOOUl/vaK6/+gOKSzLR6ZzwcAtGo9GTW1zI/LWreXvxPLwMBp678Xb6J7ZHd447hU9HRl4uo1+aysEse6qVrSn7WZu8hzZhkWw7dLDWsSUVl4RcXRBOK/xCiJ7AbMAFCBVCtAP+T0o5qaWNO1tUKhU94lvT46QZ05JN62pE/zjZhQVkFRQ0i/A7O3mh0ejZkLuE7h9Pw7w1BWtZBWEjxuDTrQdWq5kdyXU3PKdkpxMeUjc1sM7dnYj7HsB50HCslZWUewVz97za6SVO3jUf4q5BJaC4SmKJ9qHN1CepyslHqFSkzZ9P8IheuATUfsArP3IElUZTbzlFa0Ul0jAQIp+uCQdVV1WS/s5nNaIPULJtJ/JwHnRo2iK5s4MDA9p2pH1wOM7Ozuj1eqKiopg4cSJTp07Fr55dpg6hdp+1/8CBpC1YUCsLatrcuYSOGUPAoEGNnlet0iCEilCf9rho3TlavBuTqQIh7K6+1Xt2MvnT99mReogRnbsz49Z7SLjAkR1BXj58MulRHrn6ekwWCzGBQfi6nWXSv3PAxdkLF+fakW8L1q1h4kdv17we+vwT/PPyu/Ru4TQNe9OP1Ij+cZZsWsdPT0zjm79OfL9UKhXdL6PF2jOlKa6et4EhQD6AlHIH0LfRHhchgZ5edXaxehlc8XFzb5bx3d2CGTJgClarib/Tv2NXZCp+992CX+++qHU6hFDhaqi7AGlw9gT9CT+k2VxFWsY2Vq39mLXbfuLfcgvD1/jy0F9WXhjkxqPdKpnWu5jhMZKHerkQ6Gq/d0d7aXhvpDsGnaDAWkxgwm7cIlQIWx7tptxAzFWZqNS13TEO/v6UHzmCW0JCndz7fnfciyr8ITB0qmmrKquk/HBqnWsw5+XVaWuI8vJyZsyYQWRkJB999BEA48eP5/33369X9AGMEXGETbgLQ2wsxbvrPr7nHTrEpgPJGM3menrbMbhFMzj4Tlw/34Pxqbm0SQ5hUNuHMLj4sC/9KEOff4Jthw5is9lYvHEtd7z3KoUXQcIuV2dnusTE06tVmwsi+vVRVFbG24t+qtUmpeTv3dtb/Nwadd25qhCCMB8/3pvwAFEBgfSIT2DZtFfpEhPf4vZcqjQ1SVvaKaJ58TpGG6B1SDjvTXiAhz79AKvNhqNOz1cPP0WYb/MsVAkhiArvxU3XfEJFZQFOTp64uwbV3GxUKjXt21zDodS1WK12gdLrDUQEt4P8n8FjADiEczR9C0uWP1szrl7nwvcjp7Kv2JNAsYWSwx9RYiqnn387+kY+CNjj97VqwchWTrQL0KG3JeJYsZ+QwCLoEwjSBKaAOjHs7omJhN9xB6nffUfbF14ge8UKrGYzTrf8H849+tQpdaj1DSPkmtEc+uyrWu3OcY0v6gKYTCZmz57NCy+8QE5ODiNHjmRgE2PMXQJjOHrd3cRaCvA5eJDs33+v9f4hYePqRyfy1UNPcuuAIfWOUZaSR9Jdj9Q83RydORtLmQ2P2J7sz0ynvKq2G2jT/mSOHsvBw+XSLbbeUmg1antqhoy0Wu3uTmdfgrOptA4No0dcAuv27alpu3vwcNqERdAlNp6br7gSvUaLy3+wfGZz0hThT6t290ghhBaYDOxtWbOaHwedjgmDR9A3oS05RYWE+vgRG9i8JQOFEHi4B+PhXs+40kKAhy/Xj36XnNz9qGUVfrpsvFKvBWsZuF2BMfpr1m/+slY3o6mM0vwdRDq4sWb9iYqXWdk72LTlc64a+DRajd6+uCoh1F0DpTlw5AWwVs9Y1QZos7hOEZcci5H8gT749rkds1lLxPOPIT3ccfaKx8u57p+Gs6OWgDvuoaqolMyFP6Nzd6fVtGfRt6k/TPRkxo0bx08//USfPn1YsGABPXs2PcVDqLsGW2QU+3NDiXjkaSrS0ihJSkJoNHjeczfvpu4H4P5P3qN360Qi/eumYag6kFzHpZX13TdETbofN6e6bionvYMiHg3g7ODIczfeylXTnqhZUPU0uNZJvdES+Lp5MOeRKfy+bRPr9yUxqENnBiR2qFmT8bqEksNdSJoi/PcC7wJBQAb2zJxN9u8LuxN1M5AhpRwhhIjAXtTFC/vu4PFSSlNjY5wNVceOUZKcjLTZMMTG4hQYiE6rpe15ioiwlJdTsH075YcP4+DthkdIJg7WX/D1vhpf+ScgwbEHaH3twl+8CpspC7OlEg/XYKLcO2G0VbDv2FqsNgtmc93IiUOpa6koz8PNtg8y3rdnnwx6AFmRjLCe5KawlmI5tgCNW49a/fMLUjiYtRrfrneSdGg1xzbPI8AvgW6d7wLnuo/JUkqWZR7hSK+O3Hj3OCyYKXBxwL8efZRS8ttvv9G1a1d8fHx45JFHuOOOO7jqqquakDiuLr5VufgZwDG2DaGLFnFox3b+ObiPF3dsYkd1npXSygqKy+tf0NO71DVS6+qKk5OONv4R3NR3IN//c6LE4Ot33FvvDUTBTr827Vnzynv8vXsH7s7O9E1oR5uw87OZPzIgkIkBo5k4rP6iRAqnpynCHyelvOXkBiFEL+DfJp7j+BPC8Vvxq8DbUsq5QoiPgbuAj5o4VpMoO3SI9XffTdH27QC4REXRc84cXOPizngsm9mM1WRCW98OX3M+lCeBrcKenbM6skVKSer337P98cdrDg0a3p+Oz41Hf/C2E/3zf4GI6XDYvqvYUZRzRfxd5P6wmIxvX0fn7k6f/90LTkGUW4rqnN7DPQSd7RjsGg5I0HiSfvQfChyuJCjkc7xECmCFY3OxVaRQVGnF3VHNjsMprE7aSYcASGw1nM3b51JRYU+9lJaxlZLSHK4b/Q7OTrXTNR3MyiA7P4duoVX8vv3FmnY9j9A+4UTE0Jo1a3jqqadYs2YN06ZN47nnnqN7A5vHTkdVXh5H5s4l+c03EWo1CU8+Scg110BMNJNnvlrLr986NJxQn/o3cvl2aIdLXDxl+04kaYuZOpWksiKiDc68ffckxvcfRHZhAdEBQXSMakpm08sXnVZLz1Zt6NmqkTxQChctTVncfb+JbXUQQgQDw7FHBSHs36QBwLzqQ76iBdI/ZK1YUSP6AGUpKRz54YczHid/82Y23HMPf191FQdnz6YyO/vEm8ZMSH0BCv+wV8s69DSU2je6lB0+zK5p02qNlfHrX5QcOjWvndWeJ94hHBxiwTGG4mUbOTrrK6zl5VRmZJD80DOI1HwKi44SHtqtpqdaraV/7wdxzPsakFh1YSw2TeeaOcVM+W4Z76wuYq+xC8nFfmT5vkmO+wOkFVnZmnKAPk8+yAOfvMe+7FK0Wqca0T9OcUkGJSVZdT6P0ooKesUEsP/Ar7Xa1234mKLiTHbs2MGIESPo06cPKSkpfPTRRw2nSW4iOX/+ya5nn8VcXIypoIBtjz9O7po1xAeHsnjqDEKqd+x2jW3FnEemNJgW2Obpjua5p3F75mk875uE85uvsSvIjz5PTea6V6dRYTQytFM37rhyKH0S2uLsoLh5FP67NJadswfQE/ARQjxy0luuNL327jvA4xxfgbS7d4qkrNntk47dhVTf+ScAEwBCz7BYef7GjXXacteswWY2N7kiUdGePfwzejTW6s0p2x9/HGN+Pq0ff9w+EyyvXuZIf9u+eOo+AEo2g1MrrBUVWCvqZqw0l5ng1LQ2unAIuxKs5VRkHSDjuzl1+lUlHaKqnZne3SfQse11GI1leLgF4anJgzL7efbq7qN0n2SGWQ1rt+Hl5kXegT1szfgNgD49JuEfJvj077U1O2ef+fEXlj11f53zCaGqt3KXwcmJw0frlnI0mSswmcqYMWMG//77L6+88goPPPAATk7nlhZB2mwcnlP380hbuJCgkSMZ3KEzG9/8kKLycvzdPXBvZCF226GDDPrwDbQaDc56B4p2/Evb8Ciu7tabn/79mz+2b+GeISPOyd6LASklmfl5qNUq/D3OPfmgwn+Txmb8Ouyx+xrswn38pwS49nQDCyFGAMeklFtOd2x9SClnSSk7Syk7+5xmc86pBAyqm0QqeNTIJol+RWURR9I2k7t9U43oH2f/Bx9QWZ35EFMG5P1MhW4CZdonsRor7KUBzcU4eRbh1bVTrb5qZ2cMUadsctKHg84d9t0OBx/AWvYHjsF1/cruIREMHTgVb88IQoI6EB3ZBy9NNmLHFfY0Ce4DkaZEXF9/m8I531G4bRsHZ7wOC7YT7G1/FF+36TMcySGn6MTsPrMgjye+W0yr+FF0ChlJL98xdAgZRteO43F3q5vvxcPZQJlFi1pt/xyLi43M+34f5SUuGFx8efvttzl06BBPPPHEOYs+gFCpcI2NrdNuiDmR4sLfw4v44NBGRR8gLc++O9pssVBUXgbAztQUYoPsC/H/hV2eOYUFvDr/exIeuIN2k+/h0+VLKK6+VgWFk2ms2PoqKeXz2JO0PX/Sz1tSyqbkOugFjBJCpGJfzB2AfZHYXQhx/EkjGPuCcbPi0zWM8JtH18SmBw0fSNCgtqftV2Us5d8Nn/Lzb09QYSqu877G0RFRnQrCUl5C6t7/Y8VNC1g28gO2zAqisKwr/x44xPzt24l47lGCRg5BpdPh3r49feZ+jsGvCFr/AF4jwfsaiHoNDj8HKmc2O89k7NxcbLfcgjjpBuUSEYpPogd6/SlrDIUrAQk534LnVTjn5mHMrJ1CKWfRb4Tq7Yu0FosRk6mCMT1q187ddvgwkaURFEyeSepdz1Ly2Gf4FIUiRd2HQW83N1ycA/EJvpGVywt56ZkNrFuThZDtcXR0JzAwEA+P5o01jxg/Ho3hhKjrPDwIGjkSaTWSl72FA8kLSU1ZSllR43+S4fUkNusUHcueo6kADGh36edf/2XTOp76+lOKy8s5VlTIhA/eZHXSrgttlsJFSFMWd2cLIa6TUhYBCCE8gLlSyvoDpquRUj4FPFXdpx/wmJTyFiHET9ifGOYCtwGLztr6BnCyfEf7W1OJuf4hpA1cnNajMb4Gspe9jm0D5BccIWnfMgAq/dToAwMwZp7wdbd59lkcqzcZFRwNZvNjd9a8d3T+UmzuAUwoW0VSdXzzN/feyagHu+CgyUCX0AM0LnD0NSzCgFEXg96Yi8ZWzhHD4wx7azG5xUXsy8jh1VdeIrisnBBfB3xCU3EOrSdETVPd5nEVlcZyqkx1Z3ZCpUIKe7idh3sori6+9IoPZeGUF3l+7ldUGI18d9Nd7Lr5lhrXVFV2NnvuvRv9b38QHBteZ8yk1Rt4/aXplJaVcmVwMI/cey/+g2+pc1xz4dG+PQN+/52i3bsRKhVuiYm4xsSQnraGn3+fXrMnIsA3mqH9HsXgUfcJAezpDl4afxfPfvsFVpuNYG8fxvUbxGOff8T4/oMZ2LZDi13D+aDKZGLW70vqtC9a/y8juvSop4fC5UxThN/7uOgDSCkLhRBNy4FbP08Ac4UQ04FtwGenOf7MUTuhqVqHG+vszzRVgPNoTreWbTpJPDdnLqH76w+iSsqBY2UEDh6C10nFuwuT0ur0P/bzIkY/MKFG+G+f9SVbpw6jrUgGYZ/F5xPLpvQCjmRtJNA3lu6t55ByNIvc4iIA0vLzuHmhPSHZuqeuJjw0uNbu2Ro8rqTA9gYZX6RSWZiH3zVtcWkVT9neE1ErYXfeBn4e9A7+PwL8EtDrDahUaq7u3psBbTtgtdnI+3dznfUIU0GB/emhWvgtFgtqtRohBEcPHCJOreJ6g4GI8nLK33yToiOZZMx4jSDv5qltcCqucbG4hjsCKtAHY6w8xppN39SIPkDWsYPk5B1oUPhdnZx57OobGNm1JyXlFQR4elJUVsrmtz4hJjDokl/M1arVxAeHsulAcq32qAAlJFWhLk0RfpsQIlRKeRRACBGGvTJXk5FS/g38Xf37IaBFa96ZPK/BbDRTQCQ2KfFUZaD3HY5ONC787m5BaLWOmM2VWK1m/j36I4HRbRh1/3T0+to+ZI133cqUTpHhuHmecHVYbTaKqlQQ/Tio9VRUFvLb+sUUFB4h0rszocVhpM/fR0R0PBO79+Wj9f/U9FWrVNjcrsDsm4BWY79pZObnse3QQYrKy+itd2LrLS/XrEMc/f5Heiz4gfytmynZnYyhZzu8e/Xlt40vU1FZiFqlpW/PSbSKHYJWq8e1etNSlZ83Qq1GWk9sxlY5OGDw9cZms/Hjjz/yzDPPMHPmTAYPHszDV48laeGPcNKu3vz5c7HeOxG8WyC0z5QDmbMg/U1ADWFPYzKMpLC4btRRRWVRo0PptFoSw/4btWJPRa1Wc//wq1mw7p+aXci+7h7KbF+hXpoi/E8Da4QQq7CXXuxDdbTNxUpyiT+7kirJy/8ZADfXQDp0D6TdaYIcPNxDGD30ZVb9+z65+YcIC+5E707XoS/+HZxiwbmNvQg6IKL9cW3XhpId9kVBlV5P0AO382fuibq3Pm7uhIcPALcEAIpLsigoPIKPexT+22wc+HBqzbE3jB6F+9DRlJWVseboYbq2SmDy7I8Z0aUHD4wYQ5XJxLi3X+KvndvxcDHwU0z7GtFXOzrS+oknSPv6eyoyMvAdfRWO7RJYse0dKirtCc2sNjN/rXkXP584/HxP7GfwbR1HqxdnkPT0k/Y6AioVCa++zpbUg0y55Qa2b99OYmIiDg72cCQnff1/Mo7aFop5L1gGR1868frwFJwSEoiL7MnufX/i4RZClamUiooCvDzCW8aGS4Susa1Y99pMdqSmoFGr6RgZQ2xQ8xRlUfhv0ZS0zMuEEB2B4ztwHpJSNj0r1wUgL2cTefknFvuKSzJJP7qS1qFRaNWNC1RQQCJjR7yJsSoXx9zP0O3vA0j72kCbReBhzy9T7u6B7tERRJdcD1VmrH7OZLuXs3e7XfjbhEYw676HCPXQYSnKpDApjYpDB+jnOhangBCOfPQWGoPBnjpZryegU2dGbt5KwcZN3HDFFRQHhvHp8l/ZfHAfXWPjsVpt/LVzO2B/GsBij4h1DAwk8fnn2TFlCsbqm07Bxo0Ejbka1+FeFJ9SM6ekLKeW8Kt1OuJuH4931y6UZ2bhEhLMfdNf4Kd584iIiGDOnDncdNNNqKpn+L6JrXCJiaHswInPN+Smm/GIaYEd0dIG2d/UaVbnfEXH2OcIygsib/lf6AOj8R9yJb5+7ZvfhkuMxPBIEsP/m081Cs1HY3H88VLK5GrRhxNVt0KrXT9bW968s6OsqG4qobKCnaiw0pSHHAcHAw7l/0DWmycapQUOPgLt/wKtJ4crw9H7xLMudzYmcwVB4grKVTfxzr3eTLmpGB99KaHHJiF353F44w3seO7dmqGCrr4a/yFDsFVVIa1WpMXCgU8+oTLDHuBU+d13uOzcyYThA/n437/Zduggkf4BNf3zSoopiY3GJSaGsOuvp/TgwRrRP07Gz4uIvv4N0nJrR3W4OJ9SpAVQ61Tk6rXEDrkSrVbLiJEj6de/P3fffXed3O9OAQH0nDOHtEWLyV+3lsCRowgcdCUah6bV3T0jhArpkogoWXOiTeVAptN4ildsYM9DD9U0Z3z1Az0W/4pvm9MnjFNQuNxpzOn9aPW/b9bz80YL23VOhId0q9MWEd4XdT0pXRvEfOouW6AqxZ5XB/B0cWJXdgLdPe+jn9f96LVjSS33xklVQYAqiTDbMtSunSmzXs+ulz+pNUzGzz+j0mjY//775G/ahEtUVI3oH6ds926u8LaHIAZ5ehPi5VMrhcCj/ywn7qXpJL3ySr2pBYRGg7NnK3Raezy9ECp6dr0LL8+T8qlU7OPwqvu5dVQwiYmJfPWp/b/11ltvZdKkSQ0W/HCNiSHhsUfpO38+0bffhlNQvXvwmoUS91tBc3zdRLDN+X0+XrqPw2/U/hM0FxWRt217i9mhoPBfokElPF5gXUrZ//yZ0zxEhnQkJ24YSfuWApKoiL7ERfU5bb9aGDqC2g2sJ8Xze18LOns4Z4I8RtVHD5CycT0AWnd37vvxOzwOPYS6Yrv9eH0YFtVL2E5J+Qv2IugABZs2EX3v/9U9vxBEBAXSJ6Et/ybvZu7qlcz93zM8+vlHZBbk061VArK8BGm1UpmVhWurVpTsPfGkEz5xEhuK8mjT4X48HUz4qNPxdCxGI+wuopyju3jp8YF8PD8XtQr+d6OasYl/gXkSaOtPe3AheGljGD19f6d9wDFUGjemfDaPNm6e2OrJvW81n778o4KCQuOunrGNdZRSLmh+c5oHFxdv+ve+nw6JY5HShptrIDpdE8P1yvdB3nwoWgWhT9ijSjLeA69REPY0qOzpX0s2b6CoWvTBPuM88P77dL2v+jy6IBAanN1ycG+bQNHOE/nDNS4uyJOqgam93fEdOohjS1fUtAXfdiPaIHcCPb34tDo+O9DTmw2vf0iV2UyQpxfl21cBcPjrr4m9/34CBg2iIiMD70GDeO3Abma9b3dV6bVa/n7yenyz7wfX9uDWgzHX3sTGrbncNUzFs7dpCPJ3J88pkpyCdPz8Lg7ht0mJq86Eoyhk/fbFqFRqHhjUg/eXb2TY+FuoeuOtmmPVjo64JJ5+k56CgkLjDu+R1f/6Ys/Zs7L6dX9gLfbC6xctGo0eb68zTBNrzISk66Fyn/118SrwvRk67wJ9UK189mUHD9bpXrhtFxZ1P7Y53cjCnQUUlRu5rosXnd6YTNKbP5D9x9+4tm5N+M03s/e11wAwJLSiwsWE5cb2xAzqhSk1E110MMfcCtmalMQP1fVWARZv/Jfp4+6s8fdXhEUR/eTDHP7gU6TVisbVleBrx7JJq2bWnNUnLsts5rN/DrE+w5PbAw/i7taDd1+djFvmfcSGqCh16MO3Bdcx9cu/MJqf48lrbuamvgOarTrZ2aISglGhKfz11wsnGrN3cu/ACXz59w7GP/M0Ytly9CGheI27G9+2LVv2T0Hhv0Jjrp47AIQQy4HWUsqs6tcBwJfnxbrzTfneE6J/nGPfQ/AjdYqYGNq3r9M9YNQI9qtb03f6HCpN9qIfH65YxS9PT2Popx9SkV9FbnEl5St+wSUyEpceHRA9Ilm373viYgayJmcumnA9VtNuukQ8xNuvfVhr/H5t2uPhYq9ylFmQT9KRXeSFpdPzi0/YMemhmgXe4PHjeOXq63D1dsRqtTJz9lzmvPA+VSVluEYf5M5o6NJrDCT9CCVr+LfqesZ9+F3NeSZ/+j5uzs7c1kA1q/OFlDYOHai7G9VWnsyWY1nM372dNx+cxJAuVxDm46akUVZQaCJNWe0MOS761eQAZ5Yu81LhZOHQelOpuRabzYCj1NZZBU/39sBr4r0Ufv4FNqMRtwH9ye3ZnUPpZTWif5wZ834k0GsyHaNjsVVZKfK8D8/hfdixfSbOpZW00/ZCu7WCq/wnkuvrzy6zG/MPOzCsUw+++NOeXTPM15+pN4zHQWd3Nf29exdeWguOejcOvP1eraie9G/mcO2Vn/BD0io+/nAeaUdzSGzXiveeHUO/q6vTJOu8IW4WlGxk4Y9189x8tHQRN/Tuj0MDC7znB1E3RxHgbvDki8k3oFVriAsKxlHfAhFFCgr/YZoi/H8KIX4Hvq9+fQPwR8uZdAFxagVOrbBY3MjYdxU7Z3yBubiYqLssxEy8D6fgEyUVf03ezc8F6dz34rM4qdR8c2gf2+d/x/Rxd9YZtlVwEJiz2LhjLwYXP9y0DqQXp5DYeiyO/6ax5+kna451HTycvaOeJ1M48e7tE5k0bCTlVVVEBwYT5HUiFFMtK/BwlKhc40jb9m2dc5ozjrH06z/QmSSzXn6OPmN7Ex8zsPbNzTECHCOIDPwe+LNW/3BffzTqpmbfbhmEELRNGMWBQ6uQ0r4molZpiY8ZSICfUkhbQeFsacoGrvuFEGOAvtVNs6SUC1vWrAuEPgBrqx/IX7eVTQ+fqC554KNP0Lm70up/T9W0RQUEcf+wAQS5mJE2EzeHdyb0QAhxgSGoVCps1Yu3kX4B3NkrgtX/2P3Uvbrdw99bv8VkrqBbyBiOTH+5lgkly39l3K1349mjF35uGvzc6he4MHcb6zd8wlV9nsK7d29yVqzggMXCD5WV3O3kRJsqM3dmZOMkBLr3P8W5Vy+Ire0KKa8oJOfYXvpHqln46D18t3YXP61bj16rZfLIay648AME+CVw3ah3OHR0A2qVmvDQbvj5nHklNQUFhRM0NbB9K1AqpfxDCOEkhDBIKUtP2+sSw2iWLE4JIn7Lj3XeOzxnLpG334nexx7O2SnMhw1rP+Vwpj3cUwgVd/Z5igAfX5Y88wrv/vITZZWVvD7+JrZumA6Ai7MPRcUZmMz2pGg6m65Ozn8AP8oJ8mj8v6aiNA0pbZSa8sgb2J63ly9nbWkprkKgGjqUvLVrMVTvtrVWVpL9828E976yxg9uNJazZsMnJO8/EUl0b59ruanfEMJ9g2gfGX2mH1+LoFKpCfBPIMA/4UKboqDwn+G0pReFEPdgL5V4fBdSEPBzC9p0wTiQb+aN1aU4+rnXec8lIohK1YkQzKqyVCqrTsT4S2lj376FeBmcSIhoy1Njb+Cegf1QCYnFYvf563XOVBlLavpo/L1w71K3YIs2xO+0tgb4tcLNNYipU17nmnunsN1UxX03j2H5z7PoHhpC9h+1vXHGrJyaYiQABUVHaok+wO6kBXQI9aV9ZLSyUKqg8B+mKTV378NeVKUEoLoIy7mkZb5oyS61klNmQ5eYgEtUeE27Sq8n+qFJfJXkiNlqT0xaZaq7KctYVUxZZRmlx/5i//aXKUv7HGNJEp3a3QRAQdFR/H1OuG7KbCV4PXITfiOvQuXggFuH9sR9NINyt8Y3IhUXFxMSmEhUeE9cDW6MHzGYuUOHcr1Z4Fwlcb+ibmFz9ZBBNamfAczmuk8aUtr4e+cmFm34t9HzKygoXNo0RfiNUkrT8RfV1bPOKC3zpUKAqwYh4L30MDp8+T7tPnmf1m+/Rux3XzBl207UlmyOldnTFxs822BPVnqCwIirKS5J549Vr1NRWYjFUsW2nXPR6V0ICemFlJLc/MMM6PMI7q5BWG1mtuWtIGukF36fPYnpwR78mTUHR8e6KZ8BSkpKmDZtGsHBwWzbmkS3Trfx1sQnGfnvZowr/+bYH3+RfO8TCKkh7p03cYmOxiU6mqCXXyI90I9QnxNPEu5uwTg51q6W5e4WxsZDWdz69sukZDV7YTQFBYWLhKb4+FcJIaYAjkKIQcAk4JfTdRJCOAD/APrq88yTUj4nhIjAXn3LC9gCjD/5xnIhifHS8OpVbpSWlPLr9pmUltkLi6sOari2ywQ+Wvkt13WaAmhwco2lddfnOZYyB5OpFL/wq/EO6E3usT/p3P5GhFBx+OgG8vJT2HdgOf36vUhM/PU46F0JVOuIGNAek6MzNpcu7Nsynaz8JNRqHbHtHqREBnFyocCqqio+/vhjXnrpJfLy8rj22msJCAhAp3Mkc9Fv9nTKJ3H0y695s10cvmOHAfDjyiW8H/torQ1ZrgY/Rg99mX83fkbOsb14eidQqUvk0x8+w2K1klVYQFRAy+XgUVBQuHA0RfifAO4GdgH/B/wGzG5CPyMwQEpZJoTQYs/pvxR4BHhbSjlXCPExcBfw0VlZ38zoNIJr2jixc+9B/tmTU9Nus1koOfY3WpUzVVXF4OpInK8jRmsXcm1xCKsZFx93wj2r2JNdxvZdC5DSSmzUAPx84qioLMTLxQlftQNH583jr3feQaXTEfPEFL619SA+5HX8HAqptLnw4lZXnnUXxFY706SU9OjRg+3bt3PllVcyY8YMunTpAtgLiJdb6uasEWoVt/Tqyo3vnMgI6uVed93A1yeGdh3v4/PlP/PtmmTW7J1FpH8gE4YMZ1/GUfJLS+gYGUOIz3/Ss6egcNnSqPALIdTAHillPPDpmQwspZTA8VqG2uofib3o+s3V7V8B07hIhB9ArRJUVhcvORljVQGdIhNxc7LvnBVC0C5QR7vAE7H1+1M2smX73JrXyQdW0LHdDfh6x7B20+dEZQSw46kTIaE77p/ILTO/5bZN7ZAcT7ss8XFRsXTpUoYMGYJKpeLxxx/H19eXgQMH1rJp95HD+AwaSP7cH5CWE+sCnjcMpUh1YhH5/pE3Uo4/u7NNtPGvvSErxN2bka278OeOPTjq9DwwYgyPff4R1upw1M7Rccx/6vlabqKLhZyiQlKyMnBxcCQ2KLhmc5uCgkLjNCr8UkqrEGLfyaUXz4TqG8cWIBqYCaQARVLK4yqVjj1K6KLC0a11nTZv/57EB/bGvTplQn0cOrK2TltG1g5KSrIwVZVz+Mvf657r39+4bVQXvtxmQgBjXLYy4eoX2bxpI/PmzeOaa67hpptuqt9OnZ6/qwoY+vnrlCxfj81kwm1IT7ZVrSbecyyfPPgMjg4erMn24fHlRtwdTCy61Yfw6lDR7O27OPDuOxRvXM/bAwfi9+R0RnzxQY3oA2w+uI/NB/dfdMK/K/UQ1746jf0ZaQgh+N+YG3h87E14udZTmF5BQaEWTXH1eAB7hBAbgfLjjVLKUafrKKW0Au2FEO7AQqDJ2y2FEBOoLvEYGnp+M0QcrIigVacp5B35keCADnh7x1GhCiE2qPHKRt6eEZyS6Qd31yCO5e1Hq3bAPyQYNm+u9b7W0YFR1g0ERhj5+sOPeGfVH4SEhPDZZ58xevToRs+XEBrOJ78vIbaDG4U9K1Cp1BQc+w61WouzWxwzlrrUWoUvqpIcyDMT7qGhIDWdzeNvpqq6DkD6nDmU7N5N/2EDSUo7Uus8ReUX15aNSqORaXO/ZH91UXspJa8tmMvAdh0Z3KHLBbZOQeHipylRPc8AI4AXqF2MpclIKYuAv4AegHt1ZBBAMFBv+IiUcpaUsrOUsrOPj8+ZnO6cifJ2Yv6RBCJjrmPfwT9Y8ddLZBycg9WY2Wi/yLCeuBpOLMs6Orjh7RVJYVEaBcVHibznLlQnVarSurvjGBBA4ZzveHTS3STv3sZbb73F/v37ufPOO9FoGr8v+7i58/QNd+Ps1ZXQ8MEIoSI4uDtd+ryKziUcXT0bb4+XnizYf6BG9I9Tsn0710bW3hUrhKB1SHidcVKyM9i0P5nMgvNfhTO/tIQV27bUad+fmX7ebVFQuBRpLB+/A3AvdjfNLuCzk1w0p0UI4QOYpZRFQghHYBDwKvYbwLXYI3tuAxadvfnNT3laGp7r1nHv7l2oK/wIMsRxoHI9qUdW46DTM7DvY2g02nr7enqEcs3It8jLP4TZUklVVSlrNszC0z2M/n0m4xfYjg5vvkl5Sgo5ZWV8uWYNw956C98O7Wk9ehhT77yHUb2uOCN724T4E+TlS2ZJewwBw3HU6QjzcsBRI7i/p4E3V5+Yrbf21RDvY7ddONRTn0ClIiw4nHsGD+eLP5cR4OHFexMeoGNkTM0hZouFhetXc/f7b1BaWUGQlzc/Pv4cPVu1OSO7zwVPFwN9Etry2+b1tdqj/APPmw0KCpcyjU0pvwLMwGpgKNAamHwGYwcAX1X7+VXAj1LKJUKIJGCuEGI6sA347KwsbwGMeXlsuu8+8tacqPEaeu/tFLaLIK/oMPtT/qJHlztxNTTs73Y1+NW8bzJXEh7aFZ3OCUcHN5KOprK/MJ/PZ8xgWWUlEgh3cSFw5Ah2LFuAt6fXWdnt4aTCw0mPPXL2BOM7ONHKR8O6oyZivTX0CNXjb7A/BrjFxeIxbCSFv52IzPW/4//wbtWamR3a8fT143DU6fF1rx3rvzftCDe/Mb1mHSAjP4+b3pjOhjc+xN+j/v0H54rJYuZwdhZCCCL8AnBycGD6uDvZdugAWQX5AEwYMoJOUbEtcn4Fhf8ajQl/ayllIoAQ4jNg45kMLKXcCXSop/0Q0PVMxjpfFO/dW0v0AdJmzyHmi2nkFR3GzTUQrbbxSl7msjLKDh5E2my4REbi5m6P1jGbzdz94P1sXLocabYw0M+PcUFBtLn/PpL9fBjRpTutQsKa9Xo8HNUMinFkUExdm739PIl9/iXyho+h/MB+DAkJeHfpjMFgr9Eb5utfpw/AkdycWou/AEdzc8gsyGsR4U/Py+Xled/yybJfUKlUPDTyWh65+lo6RMaw4Y0POZiZgYujI/HBoRgcnZr9/AoK/0UaE/6aAHEppeVyyN1iMxrrtEmLBZVVoFKp6d19Eo4ODUeNlKens+OZZ8lc9DMAXn360PmttzBERWGRNnZt2YrV0w1dYixuvfryu0rFsvxMbm7TmvcmPICHi6GlLq1eQqKCCY4MwmiROGibstxDveLuaXDFy9Ay0TS/bFzLh7/ZvYFWm43XF86lbUQk4/oNIsTblxBvZY+BgsKZ0ti3vZ0QoqT6pxRoe/x3IURJI/0uWQyxseh9awuJd79+GH1D6NT9acJDOjXQ0076H3/ViL5VShb88Qfte/YkIyMDR52eydOfhy5tMDnomLdlPd9tWkuUfyARfv4EeZ3fBezjCCGaLPoACSHhvHzrPTWvtRoNnz/weINPCOeC1Wplzqq6pR9+Xr+mnqMVFBSaSmOlFy98MvbzjHNoKH3mzWPfBx9QsHEjDv37sTchjvQD2Uwc2hWVqvGPJPefv5FSstls5vvKStJsNmLUao5lZxMUFMQdg4ex+fABft+6CSEEY3v0oWtsPLGBwY2OezHh5ODA5JFjubJ9J3KKConw9ScuKKRFzqVWq+kZn8DavbtrtXeJVoqwKCicC03Nx3/Z4N6mDV3efx9LWRklKkGM2USgpzcq1elnxfqOnXn2m69JtlgIUKl4xNmZ6x9+lDbtOwL24i0/PPYs21NTMFssBHl5ER0QjPY0YZsXG456BzpHn59iKLcNGML3/6wkI99eWjLKP5DR3Xqel3MrKPxXubQU5zyh0mrReXjgffpDAUhNTSU8PBxL2wQSQ0K4IieH/jodrq1bY+nTsyZ2HsDNxYUr2rRrGcP/g7QJi2DNK++x++hhVELQJiySUCV3kILCOaEI/zlw4MABnnnmGX788Uc2bdrE+xv/wfXu8Vzh5QsSlpcVsW7ZQv7p1w9dA7H/F5ryqkr2ph+lqLyser0h4PSdzjPhfv6E+zX/GoKCwuWKIvxnQWZmJi+88AKzZ89Gr9fz9NNPEx0dTa/sBO7/5L1aGxOevPbmFhf91JwsdqYewmqzkRgWSXSgPf2RqaSEwi1bKNyxA+fQULy6dq1VML6wtJQZ877ljYU/APbonCXPzKBHvFLmUEHhv4wi/GeI0WikQ4cOFBYWMnHiRKZOnYqfn33D1lWdutoXI5P3ANAqJJRb+w9uUXv2ph1h6PNPcuRYNgBeBlf+nP4mbcOjSP32W3Y+/XTNsV49etD9889xrLZ32+GDNaIPUFBawn0fv8Of098676GlCgoK5w9F+OvBarVSYazC4OQMQHl5Od9//z133XUXer2ejz/+mHbt2hEZWTtpW5R/EAunvMje9CNYrTZahYQRcJa7cZvKog3/1og+2PPYzF7+Ky8PGc2el16qdWz+unWUJCXVCP/xBdOT2XboIIWlpYrwKyj8h1GE/xR2Hk5h5m8/syZpN2O790Lz/+3deZAUZZrH8e+PbhSbGw/uBlrw4FYRcREVBEeRFW8kXAePWZwNB3VWXVHDGETxjpnVnRk3PMEdUMeDGXQ3VARR8QAPblDxQESBdhQQARHoZ/9434Lqthtauuvqej4RHZ35VmXmU9XZT2W9mfm8q0q5/977WLt2LYcffjj9+/fnzDPPrHL5g5o1/0mZg5TGu+LTn7S98/GHbD/hB3Zs3vyTx7YntXWopNRy/8O7c0DTprUbpHMuq1T/zp08sLJ0LUPHX88DLzzH0rfmcOuvxzDuhhs5uHNn3njjDfr375/pEH/ijH4/jWnUwF/QqLiYNqedVq69sHFjGh+66zLMI0q6cMcvR1NYEO5P6HBQS+4bPYYm8ZuOc65u8iP+JMtWfR66P8oMPvgU6teHvj3402MT6V3SOdPhVWpgj96Mv+ASbn9qMtvLdjBm2FkMP6Y/hUVF9Lz5Zorat2fV1Kk07d6drtddR5POu15H46Ii/v2McxnW91g2bNpEp5atUt415ZzLPFmFgbqzUZ8+fezdCgOY1LbZs2dz0y3jmVX4AxQUwOYtsF8DCgoKWHjfw3Qt7pjS7ddEWVkZK0rXUGZGx4Na7TyCB7AdO9j6zTcUNmpEYZEXMXMun0h6z8z6VGzP+66ehQsXMmzYMAYMGMCyhYvodWCs6V60H0hcefrZdG6ddaNDllOvXj1KWrWhc+u25ZI+gAoKaHDQQZ70nXM75W1Xz6ZNm7jsssuYMmUKTZs25Y477mDMmDGs3biBVxbNZ96nyzmhey8GdO3BPvWz8+Yr55zbG3mX+Ldu3cq+++5LUVERq1evZuzYsVx77bU0bx6uxOlUVBTvXj01s4E651yKpKyrR1J7Sa9IWippiaQrY3sLSdMlLY+/03Lt4/r167nhhhvo0KEDpaWlSOLll1/mtttu25n0nXMuH6Syj387cLWZdQX6AZdL6gqMBWaYWRdgRpxPmc2bN3PXXXdRUlLC7bffzqBBg9i+PQwdnA+DyzjnXEUp6+oxs9XA6ji9UdIyoC0wHDgxPm0SMAu4LhUxbNiwgW7duvHll18ydOhQJkyYQO/evVOxKeecyxlp6eOX1JEw/u4coGX8UABYA1Q6crmk0cBogOLi4r3abtOmTRk9ejQDBw5kwIABe7UO55yra1J+Hb+kRsCrwAQze1bSejNrlvT4OjPbbSd7Oq7jd865uiYj1/FLqg88A0w2s2dj81pJrePjrYHSVMaQC7Zs3crnpWtY9/3GTIfinMsDqbyqR8DDwDIz+33SQ9OAUXF6FPD3VMWQC5auXMGIu8dTMvoCBt74W15bsiDTITnn6rhUHvH3By4EBkmaH3+GAncAQyQtBwbH+by0/vuNXHLf3Tw3903KyspY8NknnDpuLB988Xmtbuf7LVvYvmNHra7TOZe7UnlVz2ygquslT0rVdnPJyq9LmfPR0nJtm7f+wEdfreKw9h1qvP4Va9fwxOsz+Musl+nb5TDG/PNZHFHSpcbrdc7ltry7czebNGzQgIYNGrDphx/KtTdr2KjG697644+Me3wik2a+CMCSlSuYNvdN5tzzZw7O8tpDzrnUyvsibZlU0qoN91z8b+Xazh8wiO7FnWq87s9K1/DYKy+Va/tm43csWbmixut2zuU2P+LPIElcOHAI3Tt0YvlXq2jVvAVHHtyFFk2a1HjdhfUK2KewkK3btpVrr1/of3Ln8p1ngQxr2GA/juvag+O69qjV9XZq2Yrrz7mAcY9P3NnWo0MJPTuUVL2Qcy4veOKvowoKCrj8tDPo1elgXl7wHt2LOzK411G0PeDATIfmnMswH4HLOefqKB+ByznnHOCJ3znn8o4nfuecyzOe+J1zLs944nfOuTzjid855/KMJ37nnMsznvidcy7PeOJ3zrk8k8oRuB6RVCppcVJbC0nTJS2Pv3c71q5zzrnal8oj/onAKRXaxgIzzKwLMCPOO+ecS6OUJX4zew34tkLzcGBSnJ4EnJGq7bvqKSsrY9OWLeRCzSbnXO1Id3XOlma2Ok6vAVpW9URJo4HRAMXFxWkILf8sWbmCB154jpmL5jH8mOMYNehkurRpl+mwnHMplrGyzGZmkqo8zDSzB4AHIFTnTFtgeeKrb//BWbfdxEdfrQJg8eefMXvpIv52w3iaNWqc4eicc6mU7qt61kpqDRB/l6Z5+y76cNUXO5N+wquL5/Px6i8zFJFzLl3SnfinAaPi9Cjg72nevouqGoKxfmH9NEfinEu3VF7O+TjwFnCopFWSLgXuAIZIWg4MjvMuAw5v14EhvY8q13bpkKHex+9cHkhZH7+ZjazioZNStU1Xffs3acKDv7mGmQvnMXf5BxzfrScndOtF0b77Zjo051yK+dCLzjlXR/nQi8455wBP/M45l3c88TvnXJ7xxO+cc3nGE79zzuUZT/zOOZdncuJyTklfA5/v5eIHAP+oxXDSxeNOL487vXIx7lyMuYOZHVixMScSf01Ierey61izncedXh53euVi3LkYc1W8q8c55/KMJ37nnMsz+ZD4H8h0AHvJ404vjzu9cjHuXIy5UnW+j98551x5+XDE75xzLkmdSvyS2kt6RdJSSUskXRnbW0iaLml5/N0807EmSGogaa6kBTHmm2N7J0lzJH0s6UlJ+2Q61spIKpA0T9LzcT7r45a0QtIiSfMlvRvbsnYfSZDUTNLTkj6QtEzSsdket6RD4/uc+PlO0lXZHjeApN/G/8nFkh6P/6tZv39XR51K/MB24Goz6wr0Ay6X1BUYC8wwsy7AjDifLbYCg8ysF9AbOEVSP+BO4A9m1hlYB1yauRB360pgWdJ8rsQ90Mx6J12el837SMK9wAtmdhjQi/C+Z3XcZvZhfJ97A0cBm4GpZHncktoCVwB9zKw7UACcT+7s37tnZnX2hzC04xDgQ6B1bGsNfJjp2KqItwh4HziGcKNIYWw/Fngx0/FVEm87wj/tIOB5QDkS9wrggAptWb2PAE2Bz4jn5XIl7gqxngy8kQtxA22BL4AWhAGrngd+kQv7d3V+6toR/06SOgJHAHOAlma2Oj60BmiZqbgqE7tL5hMGn58OfAKsN7Pt8SmrCDtitvlP4D+Asji/P7kRtwEvSXpP0ujYltX7CNAJ+Bp4NHatPSSpIdkfd7LzgcfjdFbHbWZfAvcAK4HVwAbgPXJj/96jOpn4JTUCngGuMrPvkh+z8FGdVZcymdkOC1+F2wF9gcMyG9GeSRoGlJrZe5mOZS8cZ2ZHAqcSugOPT34wG/cRwlHnkcD9ZnYEsIkK3SNZGjcAsS/8dOCpio9lY9zxnMNwwgduG6AhcEpGg6pFdS7xS6pPSPqTzezZ2LxWUuv4eGvCkXXWMbP1wCuEr5DNJCXGRG4HfJmpuKrQHzhd0grgCUJ3z71kf9yJoznMrJTQ39yX7N9HVgGrzGxOnH+a8EGQ7XEnnAq8b2Zr43y2xz0Y+MzMvjazbcCzhH0+6/fv6qhTiV+SgIeBZWb2+6SHpgGj4vQoQt9/VpB0oKRmcXo/wjmJZYQPgHPi07IqZgAzu97M2plZR8JX+JlmdgFZHrekhpIaJ6YJ/c6LyeJ9BMDM1gBfSDo0Np0ELCXL404ykl3dPJD9ca8E+kkqinkl8X5n9f5dXXXqBi5JxwGvA4vY1e98A6Gf/69AMaHK53lm9m1GgqxAUk9gEuGqgXrAX81svKQSwpF0C2Ae8C9mtjVzkVZN0onANWY2LNvjjvFNjbOFwBQzmyBpf7J0H0mQ1Bt4CNgH+BS4mLjPkN1xNyQk0hIz2xDbcuH9vhkYQbhacB7wK0Kfftbu39VVpxK/c865PatTXT3OOef2zBO/c87lGU/8zjmXZzzxO+dcnvHE75xzecYTv6sxSa0kPSHpk1gG4f8kHSLpxETVzjTGMk7SNXu5bB9J9+3lsrMkVToeq6QDJG2T9Ou9WXcV67xI0teSHqqtdf7M7U+W9K2kc/b8bJdtPPG7Gok3t0wFZpnZwWZ2FHA9P6P2iqSCVMX3c5jZu2Z2RQpWfS7wNuEmptr0pJn9qpbXuVPSHao/EW/Wm5aqbbvU8sTvamogsM3M/jvRYGYLzOz1ONsoqYb85PhBkaiJf6ek94FzJY1UqJG/WNKdiXVJ+l7SBIXxCt6W1DK2d5Q0U9JCSTMkFVcMTFLvuMxCSVMTNd8lHR3b5ku6W9Li2L7zG4qkRpIejTEtlHR2bL9f0rtKGjuhGkYCVwNtJbVLiu9SSR8pjMfwoKQ/xvYDJT0j6Z34039PG5D0WrzBKzE/W1KveKfyI3Eb8yQNT3r/Xpf0fvz5p6T34HVJ04Clcfn/je//YkkjqvmaXRbzxO9qqjuhamFVjgCuAroCJYR6JwnfxGJprxHqnA8ijElwtKQz4nMaAm9bGK/gNeBfY/t/AZPMrCcwGaisi+Yx4Lr4nEXA72L7o8BlsTDejirivgnYYGY94vIzY/uNFmr49wROiHdeV0lSe0L54bmEO1VHxPY2cRv94nuSXJjvXkLN96OBswl36+7Jw8BFcd2HAA3MbAFwI6GcRl/Ch/Td8U7aUmBIfP9HUP79OxK40swOIRQm+8rMelmoS/9CNWJxWc4Tv0u1uWa2yszKgPlAx6THnoy/jyZ0FX0dS95OBhIVM38k1EKH8AGTWP5YYEqc/h/guOSNSmoKNDOzV2PTJOB4hbpIjc3srdg+hcoNBv6UmDGzdXHyvPgtZR7QjfCBtjsjCAkfwq3+ie6evsCrZvZtLAKWXLVyMPBHhVLd04AmChVnd+cpYJhCkcJLgImx/WRgbFzXLKABoUxCfeBBSYvissmvY66ZfRanFwFD4rezAYmSCy63VdmH51w1LWFX0arKJNcx2UH5fW5TNda/zXbVFam4fFpJ6gRcAxxtZuskTSQk0t0ZCbSSdEGcbyOpyx6WqQf0M7MfqhubmW2WNJ1QSvg8wmhXEAbHOdvMPqzwWsYBawkjedUDkre18+9iZh9JOhIYCtwqaYaZja9uXC47+RG/q6mZwL7aNaAJknpKGvAz1jGX0G1yQDzROxJ4dQ/LvEmoCgpwAaE4307xyHRdUhwXEo6w1wMbJR0T28+nctOByxMz8fxAE0JS3BDPNZy6uwBjl0sjM2trZh1jJdPb4+t7J77m5vEk6tlJi74EjElaT+/dbSfJQ4Qum3eSvqG8CIxJOrdyRGxvCqyO38QuJBQJrOw1tAE2m9lfgLsJ3UAux3nidzUSj8bPBAYrXM65hJDc1vyMdawmDCryCrAAeM/M9lTudgxwsaSFhMR1ZSXPGUXo015IOHeQOFK9lNDNMZ9wDqGy7otbgebxhOYCwhi9CwhdPB8Quoje2EOMI9lVCTThGWBkHBPgNsKH3huE4SATcVwB9IknlZcC1boMNA6K8x3hHEbCLYRunYXxb3NLbP8zMCq+tsOo+ttXD2BufK9+R3hfXI7z6pwu70hqZGbfx+mxhJOvlX1wpCWOeMQ/FXjEzCp+UFS17EWEgcB/k9TWhtCPf1g8kk+p2NX1vJk9neptudrlR/wuH50WL+VcDAwgc0ex4+KR9GLCQOp/+xnLbgFOVbyBS9IvCeNO3JimpD8ZOIHy5wZcjvAjfuecyzN+xO+cc3nGE79zzuUZT/zOOZdnPPE751ye8cTvnHN5xhO/c87lmf8HMJXFxMI0oPAAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "colors = [\"#1E88E5\", \"#FFC107\", \"#004D40\",\"#909842\", \"#AB1111\"]\n", "sns.scatterplot(x=\"age\", y=\"y_hat\", hue=\"modality\", data=df, \n", " palette=colors)\n", "plt.plot([20, 85], [20, 85], 'k--')\n", "plt.xlabel(\"Chronological Age [years]\")\n", "plt.ylabel(\"Predicted Age [years]\")\n", "plt.legend(ncol=5, bbox_to_anchor=(1, 1.15))\n", "plt.savefig(f\"predicted_age_vs_age.png\")" ] }, { "cell_type": "code", "execution_count": 26, "id": "346ce6b5-9cbf-4079-9f73-48ea24d153e4", "metadata": { "tags": [ "hide-input" ] }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 26, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEiCAYAAAD+lrlwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABmBklEQVR4nO3dd3hUxdfA8e/sbnohCamkU0MIPfQiRVRQsTdUwAJi98WGCtJUQP2JYFcUFRWpKiBFOqH33iEJoaSTnmy2zPvHbiIgJJtGQjKf58mT7M0t527g3Nm5c88IKSWKoihK3aGp7gAURVGU60slfkVRlDpGJX5FUZQ6RiV+RVGUOkYlfkVRlDpGJX5FUZQ6RlfdAVzNrl27fHU63QwgCnVxUhRFKYkZOGg0Gp9u3759si0b1MjEr9PpZvj7+zf38fG5qNFo1IMGiqIo12A2m0VKSkpkYmLiDGCgLdvU1NZ0lI+PT5ZK+oqiKCXTaDTSx8cnE0sPiW3bVGE8FaEpa9If+FNKs4E/pTSrqoBqrD3dm7Gne907b6uOrz7brOOrz9bZ87/U7wufa/b7wufq/Huxum/fZqv79q1T74M1X9qcz2tq4q92Wq22fURERGTjxo1bNGvWLHLs2LF+JpMJgCVLlri5ubm1iYiIiIyIiIjs2rVr02oOt9IkJCTo7rzzzvCgoKCWLVq0aN6mTZuIn3/+2WPJkiVuQoj2n3zyiXfRups3b3YSQrR/9913/aoz5sry1FNPBU+YMMG36HX37t2bPPTQQ6FFr4cNGxY0btw4PyFE+5deeqlB0fILFy7odDpdu8GDB4dc75ir0pkzZ3R33HFHw+Dg4KgWLVo0v+mmmxrv37/fYejQocFNmjRp0bRp08ioqKjmR48etQcIDAxs2bRp08ii/xcrV650qe5zqAxVlQsOHDjg0Lt378ZF72+nTp2aLlu2zBVg+vTp9YUQ7f/880+3ovVnzZrlIYRoP3PmTM+KnlON7OOvCRwcHMxHjx49DHDu3DndAw880DArK0s7derU8wDR0dE5a9euPVm9UVYus9nMnXfe2XjQoEFpixcvjgU4fvy4/bx58zy8vLzymzRpkr9gwQLPkSNHpgLMmjXLq1mzZvnVG3Xl6d69e868efM8gWSTycTFixd1OTk52qLf79ixw/Xhhx9OCAwMLFy5cqUHcB7g559/9mzcuHFBNYVdJcxmMwMHDmw8aNCgtCVLlpwG2LJli9OsWbO8EhMT7Y4ePXpIq9Vy6tQpO3d3d3PRduvXrz8eEBBgrL7IK19Fc0HHjh2bzZo1K7ZZs2aFRcvy8vLEnXfe2eT9999PePTRRzMBduzY4bhlyxaX/v375wA0adIkf/bs2V533313NsDvv/9eaf/fVIvfBoGBgcYZM2bEzZw509dsNpe+wQ1q8eLFbnZ2dvKNN95IKVrWtGnTwnfeeScZIDAwsFCv12sSEhJ0ZrOZNWvW1Ovbt29m9UVcuXr37p2ze/duV4Bdu3Y5NWvWLN/FxcWUkpKizc/PF6dOnXL08fExOjk5mRs3bpy/YcMGZ4AFCxZ43X333enVG33lWrJkiZtOp7vs30KXLl3yXVxczH5+fgat1nI9bNSokcHHx8dUbYFeZ5WVC7755pv67dq1yylK+gAdOnQoeOmll9KKXnfq1Clnz549Lnq9XmRmZmri4uIcWrRokVfBUwBugBb/60svBh9LMTqXtt7JNKMTWPr6S1u3mY8u76MBnglliSMyMrLQZDJx7tw5HcDOnTtdIyIiIgHuuuuu9ClTpiSWZX+lOjY8mLzDpZ43eUedAGzq53eOzKPZt9c87wMHDji1atWqxH9Yd99998VZs2Z5RkdH57Vs2TLPwcGhSm7APzn9w+CD8bGlnv+RhHgnsPT1l7ZuVGh43g8vvXHN8w8LCzNotVp54sQJ+/Xr17t07tw599y5c3Zr1qxx9fT0NDZt2jS/6Hwffvjh9F9++cWrQYMGBq1WKxs0aGA4f/68fVnO0VYr130UnJZe+nuRnnHGCSx9/aWtW98rPK9fr9ev+V7s37/fqXXr1v/5t/D444+n9+zZMyIiIsKtR48eWUOHDk3r1q1bcSv0pptuaqrRaLC3tzfv37//aGlxlMXOF18Mzjxc+v+J7OPHncDS11/auvUiI/OiP/vsuueCQ4cOObZt27bE/2tCCHr27Jm1cOFC94yMDO1tt92WERcX51CWWK+lxif+mqo2dvVc6fHHHw/Zvn27q52dnZw8efJZgMGDB6ffd999jY4ePeo0aNCg9I0bN7pWd5yVqX379jlr16512bJli+vrr7+edObMGftNmza51KtXz9SpU6ecovXuu+++rAkTJgT6+fkZ7rvvvlrV2i9Jo0aNDCdPnjy4ePFit9WrV7sPGDCg2c8//3zqrrvuyoba2dVTmqvlgmnTptX/6quv/ADOnDnjMGDAgCZ2dnYyODhYv3LlylNX7qNfv36N4uLiHMPDwwv++eef4t8/+uij6Z9++qlfdna29tNPP00YP358QGXEXOMTv60t86KW/qIhPseqIo7Dhw/ba7VaAgMDjfv27auKQ1yuhJb5ZYpa+m03Vvi8W7Zsmf/XX38V3ziaNWvWmQsXLuiio6ObFy0LCQkx2tnZyQ0bNrj/8MMPZ6oq8ZfUMr9UUUt/+/++qpS/e9euXXM2b97sevToUacOHTrkN2zYsPDTTz/1c3V1NQ0dOjS1aD1HR0fZqlWrvK+++sr/0KFDB+fMmeNRGce/mpJa5pcqauk/fO+XlfJv4c8//7zqTUQnJyf54IMPZj344INZfn5+hoULF3oUJf6qZGvLvKil33f16mrLBS+//HLayy+/nAZX7+Nv0aJFQUxMTPH/nZUrV57asGGD82uvvRZ86X569+6dN2LECCcnJydzq1at9JV1DqqP3wbnz5/XDRs2LPSJJ55I1mhq71t25513Zuv1ejFlyhSfomU5OTn/OeHx48efmzhx4lmdrsa3G8qsZ8+eOatWrfLw8PAw6XQ6/Pz8TFlZWdo9e/a49unTJ/fSdd98883EsWPHnvXz86t1fdx33nlndmFhofj444+LR3Ft27bN6e+//3aNi4uzAzCZTBw4cMApNDS08Np7ql0qKxcMGzYsbefOna6//vprvaJlubm5V93hxIkTz06cOPFcuQ92FbXvf24l0ev1moiIiEij0Si0Wq186KGH0saOHZtU3XFVJY1Gw+LFi089//zzwdOnT/f38vIyOjs7m8aNG3f20vX69euXe6193Og6duyYn5GRobv33nuLb7JFRETk5+bmagMCAoxZWVnF/fjR0dEF0dHRtWo0TxGNRsOiRYtOPffcc8HTpk3zd3BwkEFBQfp+/fplvfHGG8GFhYUagDZt2uSOGjXKpjIBN6qqyAWurq7yr7/+OvnKK68EvfnmmyHe3t4GFxcX09tvv33+ynUffPDBrIoc62pETZx6cd++fXGtW7dOLX3Nf1V1V0+NVYldPTeiyu7quZFVZlfPjayqu3pqqn379nm3bt06zJZ1a03iVxRFqcvKkvhrb4e1oiiKclUq8SuKotQxNTXxm81ms6juIBRFUW4E1nxp86PENTXxH0xJSamnkr+iKErJrPX46wEHbd2mRg7nNBqNTycmJs5ITExUM3ApiqKUrHgGLls3qJGjehRFUZSqo1rTiqIodYxK/IqiKHWMSvyKoih1jEr8iqIodYxK/IqiKHWMSvyKoih1jEr8iqIodUyNfIDrSt7e3jIsLKy6w1AURbmh7Nq1K1VK6XPl8hsi8YeFhbFz587qDkNRFOWGIoSIv9py1dWjKIpSx6jEryiKUseoxK8oilLHqMSvKIpSx6jEryiKUseoxK8oilLHqMSvKIpSx6jEryiKUseoxK8oilIDzV80kvmLRlbJvlXiVxRFqWNU4lcURaljVOJXFEWpY1TiVxRFqWNU4lcURaljVOJXFEWpY1TiVxRFqWNU4lcURaljVOJXFEWpY1TiVxRFqWNU4lcURaljVOJXFEWpY1TiVxRFqWNU4lcURaljVOJXFEWpY6o08Qsh/k8IcUgIcVAIMVsI4SiECBdCbBNCnBRCzBFC2FdlDIqiKDekScssX1WgyhK/ECIQeAmIllJGAVrgYWAKMFVK2Ri4CDxVVTEoiqIo/1XVXT06wEkIoQOcgQtAH2C+9fc/AXdXcQyKoijKJaos8UspzwEfA2ewJPxMYBeQIaU0Wlc7CwRebXshxHAhxE4hxM6UlJSqClNRFKXOqcquHk/gLiAcaAC4ALfZur2U8lspZbSUMtrHx6eKolQURamZcg2SXIOskn1XZVfPzUCslDJFSmkAFgLdAA9r1w9AEHCuCmNQFEVRrlCVif8M0FkI4SyEEEBf4DCwFrjfus4Q4K8qjEFRFEW5QlX28W/DchN3N3DAeqxvgTeBkUKIk0B94PuqikFRFEX5L13pq5SflHIsMPaKxaeBjlV5XEVRFOXa1JO7iqIodYxK/IqiKHWMSvyKoih1jEr8iqIodYxK/IqiKHVMlY7qURRFUcrHUeRW2b5Vi19RFKWOUS1+RVGUGsjgZU9ulCdGUyE6beVOW6Ja/IqiKDVIQUEWazdOJ+3uUPQhLqRfjK/0Y6gWv6IoSg1gNps4dHQpm3f8gF6fQ9beVGLXnuOll5tU+rFU4leUG9z8RSMBuH/gJ9UciVJe5xMPsnbjZ6SmnUSDL4/c/zGrvr0Jra5qUrRK/IqiKNUkJzeFmC3fcvzUGrKzjCyccxQ/n7a8+EwjNELQyKVqpiRXiV9RFOU6M5oK2XtgAVu2/4TBWMjqFfEknnVl3NjvuPPOO6v8+CrxK4qiXEex8VtZv/kLMrPOU8+tGV9OW89rr07i/vvvR6O5PuNtVOJXFEW5DjIyz7Fw0QSy805SWGDHA/dMIiykI4MfNl+3hF9EJX5FUZQqVGjI54/FUziXGIPBYGRrTDq33/YiYSGWaUmud9IHlfgVRVGqhJSSYyfXsHzV/xAaPYf2X6RV84eY/cvLODo6VmtsKvEriqJUsq3bVnA8dgEXM0/hUS+U5HO+TP/fGFxcXKo7NEAlfkW54R1ONlR3CIrVocN7+G3Ou3j752Iyarjt5pFENrsNjUZb3aFdRiV+RVGUCoqPj+OLr9/AwzuR+n468rP9GPr4ZBoEhJZ/p2aJyDYipUQIUXnBohK/oig11I3yRPK5C/uZv+gdAsPyyM/1YMDNo4lo2rHc+5NScnbhQjRn8hEmSfquXdSPjq7EiFXiVxRFKbPU1FQ++ngiTSMN5OQfo56HNy2aDaVzh3sr1DrPPHyYPW++SeqmTRidHNnY6RHuatu+EiO3uGbiF0K0s2F7g5TyQCXGoyiKUmNlZGTwv/99xO7987ipTwOy83R0av8Y0W0ewc6u/CN1CjMzOTxpEqe+/x6tmztxT77Hl5790UojR1MMtPSv3NINJbX41wM7gJIuX+FA2NV+IYRoBsy5ZFFD4F3gZ+vyMCAOeFBKedHWgBVFUarDN998ww8/fkC/AQ24ZUAwvvVbMeCWN6jnHlDufUqzmfjZszkwfjz6tDQybnmUzyJHkOvoQZdjv9PnyExajq78tnVJiX+HlLJPSRsLIdZc63dSymNAG+t6WuAc8AcwClgtpZwshBhlff1mGeNWFEWpcvn5+Wi1WnLzkigkhkFDG+Ps6MstfV4lNLhi/e7pe/aw9403SN+1C0Pzdvxwz1ec8GrOPS2cGNnDjc03T62ks/ivayb+0pK+retY9QVOSSnjhRB3Ab2sy38C1qESv6IoNUhhYSHff/89H374Af/3+u0Iu1i0dvZ0ix5B6xb3oNWW//aoPi2NgxMnEjtrFtLDm6X3TmJt49vp3diJ6T3dae5rV4lncnWlRi+E6AbslVLmCiEeA9oB06SUZZkW5mFgtvVnPynlBevPiYDfNY47HBgOEBISUoZDKUrdEvr1SssPT1dvHLWB0Whk1qxZjB8/Hm+/Ap55KRKpPUFE41vp1ulpXJy9yr1vaTJx+scfOfj++xiys9nbfQjz2j9DZLgXc3u50zHYoRLPpGS2XLa+AloLIVoDrwIzsPTT32TLAYQQ9sBA4K0rfyellEIIebXtpJTfAt8CREdHX3UdRVGUyvTYY4+xcfPfDHmqHb7+Ovx8mtGr24v4+zWv0H5Tt25l75tvknHgABeadGbWvaNwj4jgs5vc6NfYsdLH6ZfGlsRvtCbou4DPpZTfCyGeKsMx+gO7pZRJ1tdJQogAKeUFIUQAkFzWoBVFUSqDlJK//vqL7t274+Jix533htKld0ecHD3o1ulpIpvdihDlL6KWn5jIgfHjOTNnDnmeAcy76xNSo2/jtR7u3BflhFZzfRN+EVsSf7YQ4i3gMaCnsLwLZemEeoR/u3kAFgFDgMnW73+VYV+KoigVJqVkxYoVjB49mj17djPpw2dwr3+BwsJc2rS8j87tB+Pg4Fru/ZsNBk5+8w0HP/wQY0Ehq7sMZ1evYQy/yY/B7Vxw1FVPwi9iS+J/CBgEPCWlTBRChAAf2bJzIYQL0A945pLFk4G51k8N8cCDZQtZURSl/NatW8fo0aPZtGkT3XpE8OlX92KWx/Gp35Ze3V6gvldYhfaftG4du94YRd7J4xxp1JNlt47inlsi+aiTK+4O178E89WUmPitwzBnSyl7Fy2TUp7B0sdfKillLlD/imVpWEb5KIqiXFNVFZ/7/PPPSU1LYPpXgzGLs7i42NGjy1gah/eoUF97bkICu98eTdLfi0nzDGbR/Z/T5v4B/NHVDT/XG6hIm5TSJIQwCyHqSSkzr1dQiqIolWXPnj2MHTuWKVOm0KRJI557qR9HTuQhZTKd2gymfeuHKvTUramggMPTP+Po1KkYTbCyx4u4Pj6Cz/v4EO5V/mGf9QKMNO2RD+ZC0Fy/J3eL5AAHhBArgdyihVLKlyo1EkVRlEp0+PBhxo4dy/z58/Hw8GDPvmVs27ubzKwLNArvQc8uI3B386/QMc4uW87W19+C8/Hsb3YLiYPf4f/uaVaxEgsF8RA3gVtHXqSwQEDeYXBtU6E4r2RL4l9o/VIURbluyvt8gpSSZ555hhkzZuDi4sK48a/RorXkXOLfeHmEcM/tHxISVLHCZ1knT7Lm/97CuGk1ifUbsu/ZHxky7Fa6h1VgLL4hHRI+hHNfAXB0nRNH1jhz74Y2FYr1akpN/FLKnyr9qIqiKJUsKSkJPz8/hBD4+Pjw+uv/R99bgzl2aikp6Y707PIcrVrcVaGnbo25uawZ/xEXZ36FQWvPztvf5OY3n+W5Fm5oynt/wJQP57+EhI/AmAl+j0PoGPaP7FruOEtT6i1mIUQTIcR8IcRhIcTpoq8qi0hRFKUMEhMTeemllwgJCWHNmjVIaeaxIb1oHHWGIycWE9GkH0Me/om2re4rd9KXUrL95wXMbtWBrBnTORLVH/Hrej768Q3uiHIvX9KXJkj8GXZEQew74N4F2u+EZt+CY3C54rSVLe/CTGAsMBXoDTyBDRcMRVGUqpSWlsaHH37IZ599RmFhIU8++STevnbM++tlLiQdxs83gjtvnYi/b0SFjnN8x0FiXnmDeke2kuEfQf5Hn/PG471wti9nGpQS0pdD7GjIOwRuHSBiJnj0rFCcZWFL4neSUq4WQghrfZ5xQohdWEosK4qiXHdms5nOnTtz6tQpHn30UUa9NZLElLWs3TQeJ6d69Ov1Os2b3lKhp26TEy/yx+sf4LH0R+wcXEgaNpFHxwzHy7UCRdSydkDs25AZA46NoPmv4H0v1MCSDXrr07onhBAvYCmvXP5H2hRFUcohNzeXmTNnMmLECHQ6HZ9++imhYSEY5Ulito3DYMynbav76NTu8Qo9dZtbYOS3j2Zh9+0HeOamk9L7Ye7431hCw3zLH3z+KYh9F1IXgJ0PNJ4G/k+CpuorcV6NLYn/ZcAZeAmYiKW7Z0hVBqUoilKkoKCAr7/+mkmTJpGcnEzjxo257bbbaNWmAes3TyMtPZbgwHb06vYCXp7ln9zcYJL8Pn8baR+8TVDCXjIatqbpR7/xYO8O5Q++MBnOfAAXZoBwgJB3IOgV0LmVf5+VwJZRPTsAhBBmKeUTVR+SoigKmKRkQ2om7zRuzLlz5+jbty8TJ06kRVRDlq6cwInT63F38+f2W8bTKKxbuZ+6NUvJ4q3n2TXhPaK2zcXbzQvviZ9y37OPITTl7Coy5cDZaXB2qmXUTsBTEPo22Nv+3IBbFVZptqUefxfgeyzdOyHW8szPSCmfq7qwFEWxhcks2drpHszBOhLnzyTc240mvm64O7tib+eMvb2z9bsL9nZO6HTXvwRweTmJXFYnpxHariOzZs2iR49u7No3h1lzJyKlpHP0ENq3fgidrnwZUkrJhlN5/Pnh97RZMpWowhycHn6Kuz54GwcPj/IFbTZA4kyIfx8MSeB9D4RNAOcm5dtfFbGlq+dT4FYsVTWRUu4TQly/28+KolxVjt7I+3N/omnLHWiEhLTNXEiDC8euvY0QGuztnLGzd8bBemGwu/QCUfSz9UJhb++CnZ0TDnYul69r74xOW7llBMxmM/PmzWP69OksXboUjRC828yHp2JiiI3fzKy5T5KVfYHGDXvSo/MI3N2uOoeTTfZdKOTbH9bR5JfxdEs6Am060/fzj/GMjCzfDqWE1D8h7l3IPwHu3aDFXHDvVO4Yq5JNg1qllAlXtBJMVROOoii2OJ2Sxg9/vI+/3IfcaMB30ynajn2f82Et2W10Yf+5TOLSctDKfJy1+TTyNNConoFANz2e9nqMxjwKDfkUFuaiL8whJyeZQkPRsjyg9LmPNBqd9YJRdJG45IJxlQvFvxeWfy8q9vbO2Omc+PvvpYwZM4b9+/fTokULzp49C4CztyN/LXuLM2d3Ut8zjHvv+JjgwLblft9OpRmZvuQUjt9PpvfBvzB7+9Pu2xmE33dP+T8JZW6E0+9A9jZwjoQWC8Gr/3UfqVMWtiT+BCFEV0AKIeyw3Ow9UrVhKYpyLesO7Gbzpg/wNGcTuCcAli5FCjg4ciQAkSEh3NStG+6du5PQsCPbCr3ZHF/IomMGJOBsJ+gQZE/XUAe6hdoT6Wt32YQgUpoxGAosF4LCvCu+52Iw5KO3fi8szLVeMCzr5OVdJMNwjkJDHobCPAzGApvOqbDQxENDPBnmejdeXr4cPf096bcGUhjojH3yEW7q+jytWtyFRlO+KpeJ2SamrU/n/I/f02/TlziY9IS/8DKt33gVnWs5RwDlHrGMxU//G+wDoenXlqduRc2qxHk1tiT+EcA0IBDLUM5/gOerMihFUf5LSjOz/vmVtNifweBJ83XO5K1disHNkUJvN+6e9QcpGzeSsnEjF5YvJ362Zf6jduHh3NqtGy6duhEb2oGt+fXZHK9n0rosANwdBF1CHOgSak+3EAeaeOuKW+m4VCxms9lkuUAUXTwKcyk05JNwNhZPTzcKDbnMmzeboGB/mrYKx2AswGC9yJidtDgdy+SxSQtwdvIs1/EzC8x8tTWH9QtWc8eKD2ifegqvXn3o8OFk3Bo3Lt9J6c9C3ERImgVaNwibCIHPg9a5fPurBkLKkj/SCSHqW2voV5vo6Gi5c+fO6gxBUapVXn4G3y78AJGzi4zcKDr+uYfcQ4doNXEiW76fBkIwZOfR4vWl2Uzm4cPFF4KUTZswZFoqq7s2aoRP9+44RHflZEgHtuR4sPmMnjMZlh5cb2cNXUMd6BpiT5dQB0I9tJV2Q3jbtm2MGTOGdevWceLECUJDrz38ck4HS9mCh3YklPk4BQbJj7tz+GXFCXov+4jWx1bgEBxC+8mTCLjttvKdjzEDEj6Gc5+DNEODERDyJtjVL3XT8lja3XL+AzaW/fyLCCF2SSmjr1xuS4t/qxBiL/ADsFyWdqVQFKVSxZ09yPxlExCmTAxZ/Wj348/oc3LoNns2AbfcwpYfpv9nG6HR4BEVhUdUFE1GjECaTGQcOkRKTAwpGzeSsHAhxp8s9Rf7NGnCQ927o2nXleNBHdic7cbmeD2LjuQD0MBNS9dQe+vFwIEA97J3Zezbt48xY8awePFivL29mTRpEj4+PhV7Y67CaJbMO5DHZ+vSiFj7Ay9s/Q47IWk+ahTNXnwRrZNT2Xdq1sP5b+DMZDCmg+/DEDYOHMMqO/zrxpbE3xS4GXgS+EwIMRf4UUp5vEojU5Q6Tkoz63fMYc+eH8g1eROUMQCvLz5E6+1Nz+XLqWcdgeIockvZEwitFs9WrfBs1Yqmzz+PNJm4uH+/5UKwaRNn5s3DOHMmGmBAs2YM7t4d2aYrhwOj2ZzhyOpTeuYftFwIwj21xReBziH2eLuUfCFITk6mQ4cOODs789577/HSSy/h5la5DzBJKVl+vICPNmRjt2MNT6+bgntqPIF33EGr997DJSSkHDs1Q/IciBsH+njwvBnC36v02vjX4mxXdTeHS+3quWxlIXoDv2Dp+dsHjJJSbqmi2Iqprh6lrikoyGLB8smkJm3jlCGavgmuGL6bileHDnSdNQtH33/LB1SkS6SI2Wjk4t69pGzaREpMDKnbtmHKtVxQ6kVG4t29O/qWXTjk355NF53YllBITqEld0T46Oga4kDXUAc6BttTz1HD6dOnWbRoEa+88goAf/31Fz179sTT0/a+elvPa8sZPVPWZXHm8CkeifmI8CNrcW3cmDaTJ+Pfp0853g3g4irLjducvZZEH/4+eF7nGWP39bN8b72y3Lu4VlePTX38wGPA40ASloe5FgFtgHlSyvByR2UjlfiVuuRC0mEWLpuAviCdfYaHeHTnNrIWzyf4/vuJnj4drePl0wRWRuK/ktlg4OKePSRv3EhKTAxp27djys8HIfCIiqJ+127kRnVhv29bNqU5suNcIXojaITEreA8Z7f9hSl+GzuXzKRJWFC5YijtvA4lGZiyPovNxzMYuOt7umyeiZ29juavv06TESPQ2JfjOYOcvZahmRmrwSEUwseDz4NQgWJv5VbNif84MAuYKaU8e8Xv3pRSTil3VDZSiV+pC6SU7N4/n41bvyPT5EWcaSgPL5pK1q4dtHj7bSJeffWqNyWrIvFfyVxYSPru3SRbu4bStm/HXFBguRC0aoVL+2jmn83m9zRnzI174hDaHim02GmgbQP74nsEbQLscdDZ1oVxrfM6k2Hk45hs/jqUR+fYldy7/n9ok88R8sADtBw3DqeAgLKfYH6spUsnZQ7o6kPIKGgwHDRVWDehNFWY+G3p4292rRu61yPpK0pdUKDPZsWaKcSd2cKJwmhcuJX7ZowkJyWZzj/8QNDdd1drfBp7e7w7d8a7c2d4/XVMBQWk79pl+USwcSPnf5lFl8JCOgOuBTvw8+5DatNO7PZuw6ZkyfTNOXy6KQdHXdEzBPZ0DXEgyt8Onca2C0FKronPNufw295c/NNPM3HLFJz3b6Jeixa0+eEbfLqWY8YqQyqcmWK5eSu0EPwGBL8Kunpl39cN5JqJXwgxTko5rqRRPEXrVElkilJHJCYfZck/E8jOTWV93mPciicBU4ciXV3ptWQJXu3aVXeI/5FbWMhXq1ezatUq1q1bh9TrSdiwgVzrxSD226+Rhs9ootXSsW1b3Dt3I7lJJ3Z4tmZTkokp67OBbNzsBZ1C7IvvETTz0f1nNqtsvZnvtufw3Y5cyM3m/47PwH/Fj+hcXGgxZQoNn3gCja6MM2uZ8izDMhM+thRU8x8CoaPBIbDy3qQarKR362khRFYJvxfAw8C4a64ghAcwA4jC8gz4k8AxYA4QBsQBD0opL5YhZkWpFaSU7D2wkJit35IjPViaM5pXs/dh+vx1XKOi6PrrrzgH1qxElJubyxdffMGUKVNIT0/nnnvuISsrCw8PD8JuvRVuvZUWgDEvj7Tt24ufI4j/+gukcRpROh0927XDpVM3LjTuxDb3lmxOMrLqpB4ALycNnUMs3UIX3cI54deRyd8kk55n4smLK2j954cYU1MIffRRosaMwbGsQ0Kl0TLdYfx7UHge6t9heQDLpXnlv1k1WEmJ/zugtDFX35Xy+2lYxv7fL4Swx1LX/21gtZRyshBiFDAKeNPWgBWlNijQZ7Nq3UecittEnLEdO4xP8e7xr8ma9zMNbr+djl9/jc6lgo/NVrJDhw7Rt29fkpKS6N+/PxMnTqR9+/ZXXVfn7Ixfr1749eoFgDEnh9Rt24pHDSV8OR1pMhFtb88t7dvj2KEbZxt2YKtbSzaeL2TpsQK49XcABuiPc/vKDyjYsx23du1o+9uveF3juNckpaW0QuwYyDsCbp2g+c9Qr3tF3pIbVpmGc5Zpx0LUA/YCDS/tLhJCHAN6SSkvCCECgHVSymYl7Uvd3FVqk8TkoyxbNZGsnBQ25D2EQduDF/95ncxNMTR75RWiRo+2vQ58YSIJ3zTFO9yAU0gvcG0Lrq0t350aV3g0isFg4MSJE0RGRmIwGHjqqad45pln6NatW8X2m5VluRBYPxFc3LcPzGY0Dg7U79ABu/ZdWLlhJw3OH8UzJQl7T09avvsuYY8+WvYa+Vnb4PTbkLUJnJpA+ESof1eNLqIGVO+ongocsA3wLXAYaA3swlLg7ZyU0sO6jgAuFr2+YvvhwHCAkJCQ9vHx8VUSp6JcL1JK9h36k5gtX2MQHszNeJ4Ozi70nzmC/DNnaP/pp4Q98oitO7OMQDk5ElN+OmcPOhDaowXkHgRZaFlH6wouLS+/GDg3t2m6P5PJxG+//ca4cePIy8vj9OnTOJXnqVcbGbKySNmypfjJ4owDB0BKJNBk+HAiR43Cvqw18vOOW8okp/4Jdn6WPnz/odU23WGZVfOonvLSAe2AF6WU24QQ07B06xSTUkohxLVGDH2L5cJBdHS0KhOh3ND0+hxWrf+Yk7ExZOja8VvaMEY4xdLw46EYNRp6/vknPl262LazwkQ48SKkLQa3Tqx4/zjZyTpCX95smQgk74hlPHrOXsjZA4k/gdn6dK+wB5eofy8Erm0sr60FxsxmMwsXLuTdd9/lyJEjtGnThunTp+N4xbMDlc3O3Z0Gt95Kg1tvtZxiRgZ/dG0OdhraTJ5ctp0VJlomQrnwA2idIPRdCHrJciFUgKpN/GeBs1LKbdbX87Ek/iQhRMAlXT3JVRiDolS7pJTjLFs1gazsJA6KR1mTfguTjMvRTnwbx4YN6fb777iGhZW+Iykh+Xc4NdIyKiV8EgS9RHbyJdtq7MC1leWLwdbtzJB/8pKLwV5I/csyU5RlI3COANfWnE5y5/P3P8PduSnz58/nnnvuQVPe6QcrwN7DA1zKmJ6M2ZapDs9OA6m3jMMPeQvsKzBJei1ly9SLTYGvAD8pZZQQohUwUEr5XknbSSkThRAJQohmUspjQF8s3T6HsUzWPtn6/a+KnoSi1ERSSvYf+ouYLV+jta/HEv0Yko0N+eLcl+T8/DW+ffrQ+YcfsHN3L31nl7by3TtD02/AucRbY/8SGnBuavnyfbAoONAnILP3EHvwD7RJBwg1rqexOM+6afZAHDiOhqMLLvl00Brsyz/rVZUxF8KF7y2TmhtSwOd+CBsPTo2qO7KKqUAXT2lsuaR+B7wOfAMgpdwvhPgNKDHxW70I/God0XMaeALQAHOFEE8B8cCD5QlcUWoyvT6H1Rs+4cTp9Th4RPN5wpMEOtrx0abXyFzzD42HD6fVe++VPv780la+OR8aTobAFys+2YcQbNx5htGjp7J+/XqaNWvGwYPH0ZnTIGffFZ8OFv67nX0D64Wgzb8XA4eQ6rlRKiWkLoDYd6HgNNS7CRq+D27/6dJWrmBL4neWUm6/4lFxoy07l1LuBa72V7jO1Y4U5fpJTjnO0lUTycpORO/9OJ+cvJlb3NK579enyDp+jLYffUSjp54qfUf6C3DyRUhbYm3lf2tptVfQkSNHGDlyJMuXL8fPz4/p06czfPhwdDod4Adet1i+ihgz/3sxSF8BmC2/13ld/qmgkkYUlShjPcS+A9k7Lfcoov4Cz1tq/kidMuj19isArPvg00rfty2JP1UI0QjrJJxCiPuBC5UeiaLc4KSUHDi8mA2bv8TBsR4n3cfx18lwnnE9RstpI9Dr9XSfOxe/3r1L21GZWvkFLzSxKT6TyYRWq8VsNrNr1y4+/PBDnn/+eZydS5k5SlcPPHpavop3lmcZQXTpxeDc5/+OKNK4WO81lH1EUYlyD1qnO1wODkHQ9DvwG3RDTHdYk9iS+J/HMromQghxDojFUq1TURQrfWEuazZ8wvFT6/ALiOa31GEcuODM+5o1OL/3OroGDei2aBHuzUrpl6+CVv6xY8cYN24cALNnz6ZFixYkJCTg4FCBAmRaZ3DvaPkqUuYRRa0tw01tmbKwIAHiJ0DSL5YLUfgH0OBZy6gdpcxKTfxSytPAzUIIF0Ajpcyu+rAUxXbzF1kmGb9/4CfVcvyU1JMsXTWBzKwLhEUMZeKBvugLJV+kziBvxqd4de1Kl59/xsHL69o7kRKSZ8OpV8vclx9vaHHV5XFxcUyYMIGffvoJJycnXnnlFaSUCCEqlvSvpdwjippZ7xkUfbUGnQcAdk5my8NX576wrB70CgS/DnYlvJdKqWwZ1TPyitcAmcAuax++otRJUkoOHlnC+s1f4Ojgjl/LD3h9UyCBDoVM2jWOzGWLCBs0iHaffFJybfgKtvIjff/bfTJ//nwGDRqERqPh5ZdfZtSoUfj6VsOwxhJGFF12MchYb7nwFXEMp9vQTHwaGSxDNP0etYzHdyzHTFrKf9jS1RNt/VpsfX0HsB8YIYSYJ6X8sKqCU5SaqrAwj9UxUzl+cg0hQe2Jc3mBiRt09HTL4LF5L5K5fy8tx42j6YsvXnti7/+08qdA4Avl7q9OTk4mPT2diIgIevTowTPPPMOoUaMIrGGF3hDCksAdQ8B74L/LC5Muu4lcz/8MqbF2BD6+CVxbVlu4tZEtiT8IaCelzAEQQowF/gZ6YinDoBK/UqekpJ1i6coJZGadp0O7J/jt3G0s3lbIYNfTdPpyBHmZmXT95Rca9O9/7Z1UYl9+Tk4Bi/7YxdBBDWnbti0xMTH4+fnx2WeflfMMq4n95SOKlg6xTsTyrEr6lc2WxO8L6C95bcDyMFe+EEJ/jW0UpdaRUnLo6FLWbfocBwdXevWdwphNwey7UMhY+814TnoF4elJ76VL8Wh5jWQlJST/Zm3lF1SolZ+dnc2nn37K5Ck/kJer5+GHHy6+iasoJbEl8f8KbBNCFD1heyfwm/Vm7+Eqi0xRapBCQz5rYz7l6IlVBAe2I6zl6wz/20xGvoHpObMp/PID3Nu2pduvv+Lod42nW/UX4MQLlvLA7p0tQxGdbRuKeTU//vgj7777Lh06NeKhQV149aVZ5d6XUrfYMqpnohBiGVBUh3WElLKoRvKjVRaZotjocLKhSvefmh7L0pUTuJiRQOfoIWS63cejC7Lw1Bn47PBkMv/4naB77qHD55+jvVoFy/+08j+EwOfL3MrX6/V89913+Pv7c//99zNs2DC6dOlC3PnfKulMaxa3apzutraz6dE6KeVOKeU0LLNpNRdC/F21YSlK9bN07Sxjzh/Po9dnc8/tH7LHeA/D/8gkyjGLcUtHkPnH7zR/4w06zZhx9aSvPw+H7oNjT1kKobXbYakUWYakbzAYmDFjBk2aNOHFF19k0aJFADg6OhIdrcoTKGVny3BOe+B2YBBwK7AA+LqK41KUamUw5LN24zSOHF9JUIM29O31Nu/FaJh3IJsH3c/S+9sR5CReoNN33xF8333/3UEltfKXLl3Kyy+/zMmTJ+nYsSM//PADffuqiidKxZQ02fotwCPALcBa4Gegg5TyiesUm6LYJPRraxXDpytnf2npcSxdOZ70jAQ6tR9M44hHGL4oi20J+bzptJvAKS9gcnTkpkWLqN+hw393oD9v7ctfCu5doNl3lto1NpJSYjAYsLe3Jz8/HxcXFxYtWsQdd9xx7aGhilIGJbX4lwMxQHcpZSyAdTIVRam1Dh9bwdqN07Gzc+Se26egd2zFPb+mk5RtYmrhn8iPx+IcEUG32bNxDgq6fGMpIflXOPVauVr5UkqWLVvGmDFjGDhwIGPHjuXee+8ttSZ+dT2xrNy4Skr87YCHgVVCiNPA74CqhKTUSgZDAWs3TufI8RUENWjNbX3eYVeyC8/PT8VRmPj89FSyfp9JwG230fGbb7Bzc7t8BxVs5a9Zs4bRo0ezZcsWwsPDiYiIACxPyqtWvlLZrpn4reUY9gKjhBBdsXT72FlH+PxhnRpRUW546RfjWbpyAmkX4+nY7jE6tR/Mz3sKmLA6nRYu+bz0z+tkxqyj6Qsv0HLsWIT2kvbPZa18PTT8CAKfK1Nf/ttvv82kSZMICgrim2++4YknnsDO7gaZF1a5Idk0t5mUcjOwWQjxMnAzlk8CKvErNYKjyC33tkeOr2RNzKfY6Ry5e8BkAgPbM3Z1Jj/vzmOgeyIDfnyOrNjTtJ82jfDHH7984wq08nft2oWvry/BwcHcd999+Pn58cwzz1T53LaKAmWcc1dKaQb+sX4pyg3LYChg3abPOHxsOYEBrbit7zuYtF48MS+dDXF6XnY5SOP/PYdBSnosXIhv9+7/biylpTzw6dfL3Mo/ePAgY8eOZeHChTz//PN8/vnntG/fnvbt21fh2d6YjqZYvg+o3jBqpaqcbF1RaqT0i2dYumoCaemxdGg7iM7RQ0nIlDwxP5X4DCMfsgLdxFE4hIXRbfZsXBs2/Hdj/Xk48TykLwP3rtDsW5ta+SdOnGDcuHHMnj0bV1dXxo0bxyuvvFJ1J6koJVCJX6lTjp5YzZoNn6DTOXBX/0mEhXRkW4KeZ/64iDSb+OLC1+T8+AU+N91E55kzsffwsGx4ZSu/0cfQ4Dmbpxf8+OOP+fPPP3nzzTd57bXXqF+/ftWdpKKUoqRx/CXOdCClTK/8cBSlahiNetZv+oKDR/+mgX8Ut/UdjZurD3P35/H2igwaOul5Y+M7ZKxaTsMnn6TNpEloim6w6s9Z+/Jtb+WfP3+eDz74gMcee4zOnTszYcIEJkyYgN+16vgo/6FVg5mqTEkt/l1Y5tkVQAhw0fqzB3AGCK/q4BSlMlzMSGDpygmkpp8mus0jdOnwBBINH6zN5JvtufRzS+OB2S+QceQQbSZPptGwYZYhlEWt/FOvWeaStaGVn5KSwpQpU/jiiy8wGo1ERETQuXNnlfDL4Y9bLBfXl6s5jtqopOGc4QBCiO+wDN9can3dH7j7ukSnKBV07OQaVm/4BK1Gx8D+HxAe0oncQjMvL77IypMFDHc/QatpI8jPy6Xb7NkE9Otn2fCyVn43ayu/UYnH+vDDD5kwYQL5+fk89thjjB07loaX3h9QlDKQUqI3GDCbzSU+wFcetvTxd5ZSDrskmGVCCDX5ilKjGY2FbNjyJQcOLybArwX9bx6Nm6sv57NMPLUgjaMpRt6324DzxJFofX3p+cdC6jVvbm3lz4JTr9vUys/JycHZ2RmNRoPRaOT2229n3LhxNG/e/DqfsXIjS8nM4GB8LAfPxFq+x8ey7fgRTGYzsUkXaBRQubOo2ZL4zwshRgO/WF8/Cpy3ZedCiDggGzABRilltPXewRwgDIgDHpRSXixb2IpybRczz1q6dtJO0b71g3Tp8BRarY695wt5emE6BYVmvsj4ibyvP8KzUye6zpqFg7e3tZX/PKQvL7WVn5+fz1dffcWkSZP48ssveeCBB3jrrbfUU7ZKibLycjl0Ju6KJB9Hcua/KdDLzZ2WoeH4eXji4uBIPRfXSo/DlsT/CDAW+ANLn/8G6zJb9ZZSpl7yehSwWko5WQgxyvr6zTLsT1Gu6fipdaxe/z+ERsudt71Hw9AuACw6nM9ryy4SYG9k8r4JZCxeSOjDD9Nu6lS09vaQ+PMlrfz/QYNnr9rKLywsZMaMGbz33ntcuHCBfv360aiR5eKgkr5SJF+v5+jZM/8m9zOWBH8mJal4HRdHR6JCwrmzYxeiQsKJCrV8+Xl4IoSg19uvAODtXq/S47NlIpZ04GUhhIuUsvyPSP7rLqCX9eefgHWoxK9UkNQK1sZMY//hRfj7Nqf/zWNwd/NDSsmnm7L5dFMO3d2zGLLgJTJ27yRqzBiavfIKovA8HLKtlQ/Qv39/1qxZQ/fu3Zk9ezY33XTTdTxLpaYxmkycOH/28hb8mThOXjiH2WwGwF5nR/PgEHpEtiQqdKAlwYeEE+LjW+l997aypR5/VywTsLgCIUKI1sAzUsrnbNi/BP4RQkjgG2t9Hz8p5QXr7xOBqw53EEIMB4YDhISE2HAopa4yutuR0SeApMOLaNfqAbp2fBqtVkeBQfL6sgwWHclnsHs8nb96hty0NLr89BOBd9wBST/DqTdKbOWbzWbmz5/PHXfcgbOzM6+99hpvvvkm/fr1Uy38OsRsNhOfnHRFCz6Wo2cTKDRaZoDTaDQ0CQikZWg4j/ToU9yCbxwQiE5bs+pb2tLVMxXLBCyLAKSU+4QQPW3cf3cp5TkhhC+wUghx9NJfSiml9aLwH9aLxLcA0dHRV11HqZuklFzMOMPp+M2cjt9K6gNhiEIzd946kYZhXQFIzjExfGE6ey8YeNdpO/U/eBHp7k6vpUvxjPCGQ/dYWvn1ukPTb/7TypdS8ueff/Luu+9y8OBBvvvuO55++mn69+9fHaesXCdSShIvpl/Rgo/l0Jk4cgsKitcL9fUjKiSc/u07FXfTRASF4GhvX43R287WIm0JV7RuTDZud876PVkI8QfQEUgSQgRIKS8IIQKA5DLGrNRBZrOJ84kHOB1nSfaZWecA8PFujMuedJyPZdLwJUvSP5xs4Kn56VzMN/FZ/lz0H07ErXVrus6ahZN2Fex8A6QBGn0CDUZc1sqXUrJixQpGjx7Nrl27aNq0Kb///jsPPPBAtZy3UnUu5mQXj6C5tJsmPTureB0/D0+iQsN5ut/txV00kSGhuDu7VGPkFWdL4k+wdvdIIYQdlucpjpS2kRDCBdBIKbOtP98CTMDyyWEIMNn6/a/yBq/Ubnp9DvEJOzgdv4W4hO3o9dloNXYEBbahXav7CQ/tjJurL3MmBxdvs/JEAS8tvmiZCP3YR2TO/5XAO++kw6ej0Z19Di6uuGYrv8j7779PWloaM2fO5LHHHkOnU5VNbmS5BfkcToi/IsnHcT793zEn9VxciAoJ54FuNxW34FuEhOFTz6P6Aq9CtvyLHgFMAwKBc1gqc9rSv+8H/GH9pKADfpNSLhdC7ADmCiGeAuKBB8sTuFI7ZWZdIDZ+C6fjt3Duwj7MZhNOjvVoGNqFhqFdCQlqj72983+2k8DX23KYvC6LaPdcnlkykswtm4gYOZIWTwYiDne/Zit/69atvPfee8yYMQN/f39mz56Nr68v9jfIx3bFotBg4Ni5hP+04GOTLiClpbfY0d6eFiFh9GvTvrgFHxUaTmB97zp1z8aWxN9MSvnopQuEEN2ATSVtJKU8DbS+yvI0QM0WrQAgpZnE5KPFyT4tPRYAL48Q2rZ6gIahXfD3bY5Gc+2bY0ahY2G7Uexcl8UD9c7T57tnyD53jg6ffUBoq2Vw8r2rtvL37NnDmDFj+Pvvv/Hx8eHIkSP4+/sTdOWUikr1cP1P+gDAZDJxOunCf7ppjp8/i9Fk6YXWabU0CwymQ5NmPNH3tuIkH+7nj7aG3WitDrYk/s+wTMNY2jJFsYnBkM+Zc7uJjd9CbPxW8vIvIoSGwICW9OjyLA1Du+BRr/QnFfMNZpYdK+CL3jM459WcN1z2ETTlOUx2dvSc+Qze7u9CpvE/rXyTycSgQYOYO3cuHh4efPDBB7z44ou4ulb+gzJK+RWVLFi2a9tlSf5wQjwFhYWA5dmJhn4BRIWGc0+XHsUt+KYNgrBXs5hdU0nVObsAXQEfIcTIS37ljpp7t9rMX2T5U9xoE2zn5KYSG7+V2PgtnDm3G5OpEHt7F8KCOxAe2pWw4A44OrqXuh8pJfsTDczZn8eiw/lkF0rq27vy0sJH8T99EKfGYXQb74mL4xRw7QFNvy5u5ScnJ+Pr64tWq8XT05MxY8YwcuRIPIpKLyvVymgysef0CWIO7WfDof1sPnoIo8nEgPGjAAis701USDjPD2hb3IJvHhyCi6NTNUd+4ympxW+PZey+Drh0Zuks4P6qDEq58UkpSU0/zem4zcTGbyEp5RgA7m7+tGx+O+GhXQkMaIlWa1urLD3PxB+H8pl7II+jKUYcdYIBzRy5t14iFx67E022Eb+ezen0/HHsnE5B+FRo8AwIDWfOnOG9997jxx9/ZOfOnbRq1Yqvv/66Kk9fsUG+Xs+240fYcGgfMYcOsOXYoeIhk438G1DfzR13J2d+ePlNWoSE4enqVsoeFVuVVJ1zPbBeCPGjlDL+Osak3KCMpkLOnd9vHV+/hZycZEDg7xtBlw5P0jCsK/U9w2y+iWYyS2Li9Mzdn8c/JwowmKF1gB0f9HYg+vQaEr+dRdLGjQggoLuJri/vQnj928q/cOECkyZN4ptvvgHg2WefVeWRq9HFnGw2HTlIzKH9xBw+wM6TxzAYjQghaBXWkCf69qdHi5b0iGxFgFf94pIF3SNbVm/gtZAtffwzhBAPSCkzAIQQnsDvUspbqzQy5YaQX5BJ3JntnI7fTHzCDgyGfHQ6B0KC2tOp/eOEh3TGxbnEOX3+40yGkbn785h/MI8L2WY8nTQMbufCQOcEWPwbZ8bP4cDFi7gEeRE1LITgTkexdwPR9N9WfkFBAa1ateLixYs8+eSTjBkzhuDg4NIPrlSa82mpxBzeT8yhA2w4tJ+DZ2KRUmKn09GhcTNG3vUAPSJb0q15SzzU/ZXrypbE712U9AGklBetT+Iq1SFnX3VHwMWMBOuDVFu4kHQIKc24ONenWeM+hId2ISSwHTqdQ5n2WWCQLDtu6crZHF+IAHqGOzC6qyDi4DISPprFye3bETotgd3dCe9diG/LswjnRuxf7MSxTU7c8scgZn/zLc888wyOjo58+eWXtG3blsaNS58TV6kYKSUnL5wr7p+POXyA04mWIr4ujo50jYjigW430TOqFR2bNMfJoWz/PpTKZUviNwshQqSUZwCEEKFYhk0rdYTZbOJC4qHiLpyMzLMAeNdvRIe2g2gY1hVf7yYIG+efLSKl5GCS5UbtX4fzydJLgutpebWHG/3FKXIW/kr863PYm52DW7ADrYYWEtrbiENAEHgPAe97wCWKTU8FsTo1h2fCw7l48SIdOnSgffv26mnbKmQymdgfd7q4RR9zeD9JGZbSwt7u9ege2ZLnB9xFzxatadOwcY2rVVPX2ZL43wE2CiHWY5l6sQfW4mlK7aUvzOVMwk5Ox28m7sw2CvTZaDQ6ghq0oU3UPYSHdsHdrXz95Rn5Zv44nMfc/XkcTjbioIP+TZ24v5GZBtsXE/fOd+zbdwSNvSCoq4Hwfka8OzZB+NxjSfbOzSwx6vV8NW0aY48mkmU0c8cddzBhwgTatm1bmW+FAugNhew4cay4f37TkYNk5VmK9Yb4+HFz6/b0bNGKHi1aEREUUqcehroR2VKWebkQoh3Q2brolSvq6yu1RFZ2kvVBqs2cPb8Ps9mIo4MbYSGdaRjWhZCgaBzsy1ejxCwlm+L0zNmfx4oTBRSaoKW/HRNvdqdXwWGSf/2MhL/+ISnfQL1QM22eNhJyR2vsw+8D77vB6b9TGJrNZj766COCHO24x9+d9xYvruA7oBTJzstj89GDxBw+wIaD+9l+4gh6g6UKZfPgUB7p2YcekS3p0aIVIT7qhvmNpqRx/BFSyqPWpA//zroVYu362V314SlXyjVUXi+blGaSUo5bkn3cZlLTTwPg6RFMm6h7aRjWhQC/FiU+NVuas5lG5h3IY96BfM5lmajnKBjUxoX7Q/Q4rfic2Gd/Y/vJNLSOkuDuJsLvicSrx6MI77vB8fKbsUajkV9//ZWff/6ZZcuW4eTkxO7du1l6a3RF3gYFSM64eFm3zd7YU5jNZrQaDe0aNeX5AXfTo0Uruke2rJKJQZTrq6QW/6vAMOB/V/mdBPpUSURKiRw/P2H5YUj5tjcYCkg4t5vY+K2cPrOFvLx0hNDQwD+K7p2foWFoFzw9Kjb6pcAo+ed4AXP257Ip3vKEZfcwB0b1dCH63GzO/jiDE6tOYC4Ez0Zm2o1sSPCDQ7ELexDs/f+zP7PZzLx58xg7dizHjh2jffv2XLhwgdDQUDU8E4qHPa774FOb1pdSEpeUaEn0hw8Qc2g/x84lAJZaNp2bRfLOA4/SI7IVXSJa4OqkHpCqbUoaxz/M+r339QtHqQq5eenF5RHOnNuF0ajH3s6Z0OAONAztQmhIR5wcK96KO2S9Ufvn4TwyCyRB7lr+r6sz93uuJ3fBl8R+sIvNCWZ0zpKwAQ0If+wRPHuMADvva+4zMTGRW265hQMHDtCiRQsWLlzI3XffrfqQy8BsNnM4Id7aP28ZdXMuzdJb6+HiSrfmUTxxc396RLYkunEzVeqgDiipq+fekjaUUi6s/HCUyiClJC09ltPxW4iN30JisqWKtpurLy2a9Sc8rAtBAa1tfmq2JJkFZv48nM+c/XkcSjLgoIX+TXU83WgHvntnEvu/DWzbbMBsFNSPqkf0s3cS9MgodPUalBh/fHw8YWFh+Pr6EhERwVtvvcWDDz6oCmzZwGA0svvU8eJhlZuOHCyuMR/gVZ8ekS0tN2IjWxEVGl5t0/8p1aekrp47rd99sdTsWWN93RvYDNT4xH+j1rUpD5PJwLkL+4uTfVZ2IgB+Ps3o0uEJwkO74O3VsFJaymYp2RxfyNwDeSw/lo/eBK19TXzfaycdChdwYcE/xI4q5FiiBjs3OxoOupnwYaOo16J9qfvesGEDo0ePZt++fcTGxuLl5cXcuXMrHHNtZjSZyM7PY9xvPxJzeD9bjx0hT28pfdA4IJC7O3WjR4tW9IhsSUP/BurTklJiV88TAEKIf4DIonlyrbNm/XhdolNKVFCQRVzCdk7HbSH+7A4KC3PRau0JCWxHh7aDLE/NutSvtOOdzzJZb9TmkZBpwtexgPFtt3BbvWUUbFpO7PhCVu3QIU3g0ymSFuNeIHDgvWgdHUvd9/bt2xk9ejQrV64kICCASZMm4eJyY89yVNkKCgs5khB/2ZSAB+JiSUi1TGJ3ID6W1mGNeKpff3pEWm7EBnhV3t9fub5svWdTHraM4w++ZHJ0gCRAzX5eDaSUGN3t0Ie4MH/RSM4nHkBKM85OnjRp2JOGoV0JDmyHnV3pidZWeqNk1ckC5uzPY0OsHmdNDs82ieHeNivwuLCK+AWFbFntQF6KxMHLiybPPU744KG4Nbr67FZXc/z4cTp16oS3tzf/+9//ePbZZ3GqwzcUjSYTpy6cK07sRYn+xIVzmM1mAOx0OpoHhdCjRUu2HD2Ei6MTMZOmq9IHik1sSfyrhRArgNnW1w8Bq6ouJKWIXp9DYvJREpOPkJh8hKTkY+Q/GA6AvT6L6DYP0zC0K36+zcr81GxpjiQbmLs/jz8O52MqvMgDDdbybpd/CDetJXGHkdgvnUncLUA64NfrJloNHkKD/v3R2Dhr1ZEjR4iJiWH48OHF89oOGDAAN7e6U4FRSklCavJlteYPxMVy5Gx88Zh5IQSN/BvQMrQhD3bvVVyOuEmDIOysU0IWjepRSV+xlS0PcL0ghLgH6Gld9K2U8o+qDavuMZmMpKaftib4IyQmH+ViRkLx7708QggL6cS5H37G/mwuj65dXekxZOnNLLLeqD2bmkh/71X81mIFEXYbyUs0ETfbk2WrXSlIzcfR34eI/3uU8McewyU01OZjnD59mvHjx/PLL7/g7u7Oww8/jLu7Ow899FCln09NkpqVycH4WA7En74k0ccVP/0K/9ab79u6HS1DLROKNA8Oxdmh8j7BKQrY1uIH2A1kSylXCSGchRBuUsrsqgysNpNSkp2TRGLSEWtr/ijJqScwmSxj3p2cPPD3bU5Ek374+0bg59MMBwdLa27OG59VeixbEwqZsz+Pnafj6O2+gon+/9AqfCsYzJzbE8DGVeEkbz8HmkIC+vUjfPBg/Pv1Q1OGSciTkpJ49913+eGHH9DpdIwcOZI33ngDd/fSJ1+5keTk53PokvleixJ9UR0bAE9XN1qGhvNYr5uLZ4yKCg1X9eaV66bU/7lCiGFYavN4AY2wTLr+NTfCvLk1oJIlWOreJCUfs7bmLV03efmWRKDV2uPr3YRWkXfi79scf7/muLn6VfnIi8Rsy43a9UdO0MrubwZ7r+CT1rvQCElWakMOLutC/NJY9GkZOAe7E/nWW4QNGoRzYOlTIl5KSokQAr1ez+zZsxkxYgRvvfUWDRpcezjnjaBoYu/iFvyZOA7GxxKb9O/tMCd7B1qEhDGgfafi5B4VEk6AV/1K/ftW5U1ApXaypcn2PNAR2AYgpTyhyjJfm9lsIi09trgln5h8hPSLZygqaOrpEUxocDR+vs3x922Ot1dDtFrbW84VUWiy3Khde+gQ9fMWcZvXcl5ssh8Ag7YFCfvuI3bxeVK37UXoUmlw222EDxmCX69eiDKOn09LS+Ojjz7i6NGj/Pnnn4SEhHD27NkbroVvNps5nXjhipE0p686sXfHphE8eXP/4m6acL8ANUa+AtQFrerYknH0UsrCohaKEEKHKssMWFqzObkpliSfZEnyyaknMBotY6idHOvh5xtB00a98fdtjp9vMxwdrv/H+eMpBlbt3wOpf3KT+zI+8rY80KV3akdm9kvELs0mfuEKDBlLcW3YkKixYwl75BEcfct+fc/MzGTq1Kl88skn5OTkMGjQIPR6PQ4ODlWW9C/mV3wfUkoupKf9ZyTNoTNx5Bfqi9cL9wugZWg4d3fuXtxN0ywwWD3tqtxQbEn864UQbwNOQoh+wHNAnSyDWFiYR1LKMRKTjxbfgM3NSwNAq7HDx7sxURED8Pdrjr9vBO5uAdX2sEx2gYmYQzvIObeAdvZLec75FARAhn1nCj0ncm6zPbGz/yZ95ww09vYE3nkn4YMH49OtG6KcrdTNmzdz5513kp6ezn333cf48eNp0aJFJZ9ZxV3Myf7PSJqDZ2K5mPPvbSt/Ty+iQsIZ0f/O4gQfGRym6tYotYItif9N4GngAPAMsBSYYesBhBBaYCdwTkp5hxAiHPgdqA/sAh6XUhaWNfCqZjabSLsYV9wnn5h8lLT0OIo+7HjUCyQ4sK21JR+BT/1GlVICoSKk2cyhU1tIiV9AI/MSBjiewVRPQ5K2KzlBz2NIbkLs78s4M/8zjDk5uDVtSqv33iP0oYdwqF++B30KCgpISEigSZMmtGzZkltuuYXXXnuN9u1Lf0q3quXpCziSEM+BoiRvTfRFdWoA6rm4EBUSbhkqGRJOy7BwWoSEqwqUSq1WYuK3Ju1DUsoI4LtyHuNl4AhQ9Dl/CjBVSvm7EOJr4Cngq3Luu9Lk5KYUd9ckJh8lOeUYBmuXjaODG36+zWkc3hN/3wj8fSNwdKwhfdXSRHriJhJOzcdfv4Qou/MY7HTEmntwxuc1/H1uJn/xeg6P/omM/fvROjkRdPfdhD/+OPU7dSr3JxKDwcDMmTOZOHEirq6uHDx4EDc3N2bPnl36xpXMBKTY2zF341rrkElLkj+VeB4pLRdqBzs7IoPD6Nuq3WU3WoO8fVQJA6XOKTHxSylNQohjl069WBZCiCDgduB9YKSw/A/rAwyyrvITMI4qSvzXql1faMgnOeV4cZJPSj5CTq6lFajR6PDxbkxkRH9rkm9OPfeaVd9EaCTGtDWcPz0f95wleGmTcTbbc9B4E2c8RhMZMZD6h04T+8XP7PpzAqa8POpFRdHmww8JeeAB7OuVvzVrMpn49ddfGT9+PKdPn6ZLly6899571714mpSSnSePMXPVMn5sHE6+VgsfTkCr0dCkQRBtwhvzeO9+xd00jfwbqAJvimJlS1ePJ3BICLEdKH7aREo50IZtPwXeAIruaNYHMqSURuvrs1iGh/6HEGI41ikeQ0LKXyFCSki1jrJJsrbo0y7GIaXl0fd67g0IDGhdnOS9vRuh09r29Ol1Yy6EvMPIrF1EPSYJaZ6H7tAAvE1ObM7pTZ77XbRqMZBWOoifO5dNLw4k6+hRdK6uhNx/P+GDB+PZtm2lXLwWLlzIkCFDaNu2LX///Tf9+/e/rhfFpIvp/LJuFTNXL+PQmTgc7e1pkVdI8zw9r/0+j2aBITja+PSwotRVtiT+MeXZsRDiDiBZSrlLCNGrrNtLKb8FvgWIjo4u1yiiPAMUmuDXeU8D4ODghr9PBA3DulluwPpE4ORUw/pyrUme7F2QswdD5m40eQfRUogAglq6s/piXy4Y7iKiyQBuauRBxtYtxI4exe5FizDr9Xi2a0f7Tz8l6J57sKtgCQQpJUuWLCE7O5tBgwZx7733snjxYm6//fbrlvALDQb+3rmVmauXs3TnVkxmM10iWvDNcyN5qEdvfuzYAYDW4Y2vSzyKcqMrqR6/IzACaIzlxu73l7TUbdENGCiEGAA4YunjnwZ4CCF01n0FAefKG3xp7LekonOxo/sLk/H3jcCjXlCN6rK5MsmTvQeZewBhvdedY3JnX04LDuQM5VRhS+w92pM+6yciz65j2OJbiZ/9C6t/+YWckyexc3cn/PHHCR88GI+oqAqHJqVk9erVjB49mm3bttGlSxceeeQRtFotd9xxR4X3b4v9saeYuXo5v6xbSWpWJgFe9XntnocY2vc2IoL+/RT4x123ApabSYqilK6kFv9PgAGIAfoDkZTh/5aU8i3gLQBri/81KeWjQoh5wP1YRvYMAf4qT+C2cN1veTq2+fR+VXUI2xUn+d2Qsxuy90DuAbAm+QLpzrGClmy9OJT9OVEc17ciwKcxXUMd6dbBnmG+dmAo5Nc35mKXa+TvqCikwUD9zp1pPnIkgQMHonN2rpRQd+/ezciRI1m/fj0hISHMmDGDwYMHX5eLZlpWJr9tWM2Pq5ez+9QJ7HV2DOzYlSduvo1b2nZAp/rpFaXCSkr8kVLKlgBCiO+B7ZV0zDeB34UQ7wF7gO8rab81RylJ3qipxwXZij1ZT7ImKZLd2VGcKwylTYAD3ULtGRzmQJsAe+wwkbFvH8m/x7BpwwbStm3DOV+PWQONnxlO+OOP4x4RUXlhm81oNBqysrI4duwYn332GcOGDcPBwaHSjnE1RpOJf/bsYOaq5SzavplCo4G2DZswffiLDOrZl/pqaKWiVKqSEr+h6AcppbEirT0p5TpgnfXn01hKQFS5PEPp61TYZUl+j+X7JUleauuR59CGU9rhbE6P5K8zzTmcHQIIIn11dAtzYFyoAx2D7HGxg8wjR0hZsoGdMTGkbNqEMdvyUFG9yEjCBw9m/4IfMTlpaf3++5V2CgcOHODdd98lJCSEadOm0atXL2JjY3G0YQKVijh29gwzVy/n57X/cCE9DW/3ejw34C6G9r1V9dcrShUqKfG3FkJkWX8WWJ7czbL+LKWUNWQg+3VUSpJHWw9c25Dl/RwHcqJYlRTJ4lh/UvIsvw7z1NKtkQMvhDrQJcQeTycNOadOkbwhhoMbNpC8cSOFaZYngV0bNiT43nvx7dEDnx49cPTxAWDNwl+goHJO59ixY4wbN445c+bg5ubG22+/Xfy7qkr6mbk5zN24jpmrl7Pl6CG0Gg0DojvzRN/buD26syp9oCjXQUlTL9btztQrk3zOHsjZ/58kT+DzZNm1YUdGFCvPBrDxgIGETEvxLl9XDd3DHega4kC3MHsC3XXknT1L8oalnJwWQ8qGDeRfsFRzdGrQgIB+/fDp0QPfHj1wDgqq0tP79ttvi2e6euutt3j11Vfx8vKqkmOZzWbWHtjDzFXLWbglhvxCPZHBoXz0xAge69UPf8+qOa6iKFd3fcpC1nRmA+QdsinJ49qWHIe2bE0OZFO8gc279BxNsQx2cnfQ0yXEgac7uNIt1J7G9XXoU1JIiVlD4g8x7I+JITc2FgAHb298unfHt2dPfHr0wLVh5UyEXpKzZ89iMBgIDw+nV69evPTSS7z11lv4lqMYmy1iEy/w45rl/LRmBfHJSXi4uDK0760M7XsbHZpE1KwRVopSh9S9xH/VJH8ApLUC4xVJHrf2FOjC2H3eyKZ4PZvi9ey/YMAkM3DQQccgB+6KdKJbqANRfnaYsjJJ2bSOlHkxrNywgayjRwHQubnh0707jYcNw7dnT9wjIspdDK2skpOTmTx5Ml9++SW33347CxYsoGnTpkydOrXSj5VbkM/8Tev5cc0K1h3YixCCfm3aM3nwcO7u3F09XKUoNUCtTvxCI/FsYIQLP5SS5J+zJvl24NgQkxQcSDSw6ZQl0e88l4TeCFoBrQPseK6zK93CHGjXwB5tQS6pW7eS/HUM62JiuLhvH0iJ1skJ786dCXnoIXx79MCjVasyzVhVGdLT0/n444+ZNm0aBQUFDBkyhHfffbfSjyOlZNORg8xctYy5m9aRk59P44BA3nvsKQb3voVgHzV9g6LUJLU68d/6XCYBTQ1w4jnQuluS+xVJHqFBSsmJNCObDuvZHJ/B1jN6svSWh4Wb++h4rI0L3UId6Bhsj7MsJG3nTlJ+2cDmmBjSd+1CGo0IOzvqd+hA5Btv4NOzJ17t2qGt4mGQpZk6dSqTJ0/m4YcfZty4cTRt2rRS9382NYWf16zgxzUrOHH+LC6OjjzYrRdP3Nyf7pEtVVeOotRQtTrxH1rnxLHNjvT6Yktxki9yNtPIphMFbIrXszleT0qupXZPiIeW2yMsXTddQuzxcpBc3LOH5EUx7ImJIXXbNswFBaDR4NmmDU1feAHfHj2o36lTpT1AVV55eXl88cUXtG3blptvvpn/+7//48EHH6Rly5aVdoyCwkL+2raRmauWs3LfLsxmMzdFteadBx7lvq43qXr1inIDqNWJP+GgpcXdy6kxaXkmtsTnF/fTx2dYRt74uGjoEmJ5cKpbmANBbhoyDx8mec16jsbEkLp5M8acHADqtWhBw6FDLTdku3bFroZMI6jX6/n22295//33SUpK4v/+7/+4+eab8fLyqpSROpdWwpy9YQ0ZuTmE+PjxzgOPMqTPrTQKKNs8vIqiVK9anfhP+XXgRIOuzPkhmSPWkTdu9oLOIQ480d7SfdO4vpbcU6dIjonh7NQN7Nm4kcL0dABcGzcm5IEHLGPpu3fHwdu7Ok/nqubOnctrr71GQkICN910E/Pnz6d79+6Vsu+rVcK8r0tPhva9jT6t2qr5ZBXlBlWrE/+2pg9yKqATnZw1vN7TjW6hDrT0t0N/7iwpG5aRPCOGpTExFBSNpQ8MJODWWy0t+u7dcQ6smS1Zs5QYjUZ0Oh1paWk0aNCAmTNn0qdPnwr3qxcaDCzdtY2Zq5bxt7USZudmkcWVMOu5uFbSWSiKUl1qdeK/Y+eHOBVmMXjpJpI3/kPK7xtYGRNDblwcAA4+PsVPxvr27IlLWFiNviFpNpvZn53H8rQsXGfMYMSIEQwfPpwRI0ZUOO6iSpi/rl9FSmYG/p5evHr3gwztexvNg0Mr6QwURakJanXid0k4h7PeyJLmzQGwc3fHp3t3mowYgU+PHpax9DU40ReRUrJ06VLGjBnDngvp+NrrCA4OBqjQrFLp2Vn8tn41M1cvY/epE9jpdJZKmH1v49Z2HVUlTEWppWp14teaJSaNhpZjxuDToweerVohbsBkNnz4cGbMmEHDhg15xN+Tdm7O3H777eXal8lk4p+9O5m5ahl/bVOVMBWlLqrVif+im2VUT7OXXqrmSMpu8+bNRERE4OXlxcMPP0yHDh144okn+CyqYbn2V1QJc9balZxPT8XbvR7P9h/IEzffpiphKkodU6sT/41o165djBkzhmXLljFu3DjGjh1L37596du3b5n3lZWXy5yYtZdVwuzfvhOfDX+ROzp0UZUwFaWOUom/hjh48CBjx45l4cKFeHl5MWXKFJ5//vn/rKct5ZaE2Wxm3YG9zFy9nAWbN6hKmIqi/IdK/DXE6NGjWbt2LePGjeP//u//cC/jw2FXVsKs5+LCkD638sTNtb8S5roPPq3uEBTlhqISfzWJi4tj4sSJjBo1iiZNmjB9+nRcXFyoX79+yRtq/x1Hn1uQz4LNG5i5enlxJcybW7dn0uPDuLtzd5yquVaQoig1k0r819n58+d5//33+e6779BoNPTp04cmTZoQEhJi0/YSiHW04+nPPmLuxnVk5+fRyL+BqoSpKIrNhJSyumMoVXR0tNy5c2eZt3uhTRMydToG/O9zzNKMyWzGLCVmc9HP1u9meY2fL9nmkmVmKa/42YTZLIv3aZYSk8l0xXaSvfv3cfz4caSUhIaF0bhxY+wdHTCZLtnPZfH9N9ajZ8+gNxhUJUxFUUolhNglpYy+cnmtbvGv9vLkqKsLv/zvvUrZn0ajQavRoBECrUaLRiPQiBKWaQRCgk6nQ6PRkJqXg3sDf/x8/XBydCQ9PxetvgCNEJftW6PRYKfTXbafop/Pp6Xi5evOzk++UZUwFUUpl1qd+O9LSqEwJZWnlq1DU1LC1motCfeSRKvVaC9LyGVtUWdnZzN9+nQ+/vhj/vzzT2666SaklBVumfd6+xUAlfQVRSm3Wp34vYyWipzNgmzrP68M+fn5fPnll0yePJnU1FQGDhxYPKet6o5RFKUmqLLEL4RwBDYADtbjzJdSjhVChAO/A/WBXcDjUhbNan5jk1LSqVMnDhw4wC233MLEiRPp2LFjdYelKIpymaosqK4H+kgpWwNtgNuEEJ2BKcBUKWVj4CLwVBXGUOWMRiO///47JpMJIQRjxoxh/fr1rFixQiV9RVFqpCpL/NIix/rSzvolgT7AfOvyn4C7qyqGAulCgXSpkn2bTCZ+++03IiMjeeSRR1iyZAkADzzwAD179qySYyqKolSGKp1CSQihFULsBZKBlcApIENKabSucha46mwnQojhQoidQoidKSkp5Tq+i53Axa5y+9WllPzxxx+0bt2aRx99FCcnJ/766y8GDhxYqcdRFEWpKlV6c1dKaQLaCCE8gD+AiDJs+y3wLVjG8VdJgOVgNpt5++23kVIyZ84c7r//fjUFoaIoN5TrMqpHSpkhhFgLdAE8hBA6a6s/CDhXVcf9465bAXi5gvtZt24dn3zyCb/++itubm4sW7aMoKAgdLpaPShKUZRaqsqaqkIIH2tLHyGEE9APOAKsBe63rjYE+KuqYqiorVu3cvPNN9O7d292797NiRMnAAgLC1NJX1GUG1ZV9lEEAGuFEPuBHcBKKeUS4E1gpBDiJJYhnd9XYQzlkpeXxx133EGXLl04cOAAU6dO5eTJk7Rr1666Q1MURamwKmu2Sin3A22vsvw0UCPHOaampuLt7Y2zszPOzs5MmjSJF154AVdX19I3VhRFuUGo/grg5MmTjB8/ngULFnDs2DGCg4OZO3dudYelKIpSJer0cJQzZ84wbNgwIiIiWLBgAS+++CIuLlUz7l9RFKWmqLMt/vT0dCIiIjCZTDz//PO89dZb+Pv7V3dYiqIoVa5OJf7U1FSWLFnC0KFD8fLy4quvvqJ37942T4KiKIpSG9SJrp6MjAzGjBlDeHg4Tz31FLGxsQAMGTJEJX1FUeqcWt3iNxUaOLtzD+Hh4WRkZHD//fczfvx4wsPDqzs0RVGUalOrW/xmo5GErbvo0aMHe/bsYd68eURGRlZ3WIqiKNWqVrf47Zyd6Dh8MIumf1fdoSiKotQYtbrFD2DvqoZnKoqiXKrWJ35FURTlcrW6q2fdB59WdwiKoig1jmrxK4qi1DEq8SuKotQxtbqrpzZS3VeKolSUavEriqLUMSrxK4qi1DEq8SuKotQxKvEriqLUMSrxK4qi1DEq8SuKotQxKvEriqLUMSrxK4qi1DEq8SuKotQxQkpZ3TGUSgiRAsRXdxzX4A2kVncQ1UCdd92izvvGFCql9Lly4Q2R+GsyIcROKWV0dcdxvanzrlvUedcuqqtHURSljlGJX1EUpY5Rib/ivq3uAKqJOu+6RZ13LaL6+BVFUeoY1eJXFEWpY1TiLwMhhKMQYrsQYp8Q4pAQYrx1ebgQYpsQ4qQQYo4Qwr66Y61sQgitEGKPEGKJ9XWtP2cAIUScEOKAEGKvEGKndZmXEGKlEOKE9btndcdZ2YQQHkKI+UKIo0KII0KILrX9vIUQzax/56KvLCHEK7XxvFXiLxs90EdK2RpoA9wmhOgMTAGmSikbAxeBp6ovxCrzMnDkktd14ZyL9JZStrlkWN8oYLWUsgmw2vq6tpkGLJdSRgCtsfzta/V5SymPWf/ObYD2QB7wB7XwvFXiLwNpkWN9aWf9kkAfYL51+U/A3dc/uqojhAgCbgdmWF8Lavk5l+IuLOcMtfDchRD1gJ7A9wBSykIpZQa1/Lyv0Bc4JaWMpxaet0r8ZWTt8tgLJAMrgVNAhpTSaF3lLBBYTeFVlU+BNwCz9XV9av85F5HAP0KIXUKI4dZlflLKC9afEwG/6gmtyoQDKcBMa/feDCGEC7X/vC/1MDDb+nOtO2+V+MtISmmyfhQMAjoCEdUbUdUSQtwBJEspd1V3LNWku5SyHdAfeF4I0fPSX0rLsLjaNjROB7QDvpJStgVyuaJ7o5aeNwDW+1UDgXlX/q62nLdK/OVk/ei7FugCeAghdNZfBQHnqiuuKtANGCiEiAN+x9LFM43afc7FpJTnrN+TsfT3dgSShBABANbvydUXYZU4C5yVUm6zvp6P5UJQ28+7SH9gt5Qyyfq61p23SvxlIITwEUJ4WH92Avphuem1FrjfutoQ4K9qCbAKSCnfklIGSSnDsHz8XSOlfJRafM5FhBAuQgi3op+BW4CDwCIs5wy18NyllIlAghCimXVRX+Awtfy8L/EI/3bzQC08b/UAVxkIIVphubmjxXLRnCulnCCEaIilNewF7AEek1Lqqy/SqiGE6AW8JqW8oy6cs/Uc/7C+1AG/SSnfF0LUB+YCIViqxj4opUyvpjCrhBCiDZab+fbAaeAJrP/mqd3n7QKcARpKKTOty2rd31slfkVRlDpGdfUoiqLUMSrxK4qi1DEq8SuKotQxKvEriqLUMSrxK4qi1DEq8SsVJoTwF0L8LoQ4ZS1tsFQI0VQI0auomud1jGWcEOK1cm4bLYSYXs5t1wkhrjo3qxDCWwhhEEKMKM++r7HPoUKIFCHEjMraZxmP/6sQIl0IcX/pays1jUr8SoVYC7b9AayTUjaSUrYH3qIM9UyEENqqiq8spJQ7pZQvVcGuHwC2YnkwqDLNkVI+Xcn7LHbJk9n/YX2Ib1FVHVupWirxKxXVGzBIKb8uWiCl3CeljLG+dL2krvuv1gtFUZ37KUKI3cADQohHrHXvDwohphTtSwiRI4R4X1jmQNgqhPCzLg8TQqwRQuwXQqwWQoRcGZgQoo11m/1CiD+K6qgLITpYl+0VQnwkhDhoXV78CUUI4SqEmGmNab8Q4j7r8q+EEDvFJfMx2OAR4FUg0FrptCi+p4QQx4VljofvhBCfW5f7CCEWCCF2WL+6lXYAIcQG60NXRa83CiFaW58+/sF6jD1CiLsuef9ihBC7rV9dL3kPYoQQi4DD1u3/tr7/B4UQD9l4zkoNphK/UlFRQEkF3NoCrwCRQEMstX+KpFkLoG3AUt+/D5Z5DjoIIe62ruMCbLXOgbABGGZd/hnwk5SyFfArcLUump+BN63rHADGWpfPBJ6xFtszXSPuMUCmlLKldfs11uXvWOvytwJusj7NfU1CiGAgQEq5HcvTnw9ZlzewHqOz9T25tNjfNCxzHXQA7sNaDrsU3wNDrftuCjhKKfcB72Aps9ERy0X6I+vTqclAP+v7/xCXv3/tgJellE2B24DzUsrWUsooYLkNsSg1nEr8SlXbLqU8K6U0A3uBsEt+N8f6vQOWrqIUa6nnX7HUgwcoBIruE+y6ZPsuwG/Wn2cB3S89qLDUlPeQUq63LvoJ6GmtteQmpdxiXf4bV3cz8EXRCynlReuPD1o/pewBWmC5oJXkISwJHywlLoq6ezoC66WU6VJKA5dXgrwZ+FxYyn8vAtyFEK6lHGcecIcQwg54EvjRuvwWYJR1X+sARyylB+yA74QQB6zbXnoe26WUsdafDwD9rJ/OehSVMVBubNfsw1MUGx3i32JtV3Np/R4Tl/+by7Vh/wb5b12RK7e/roQQ4cBrQAcp5UUhxI9YEmlJHgH8hRCPWl83EEI0KWUbDdBZSllga2xSyjwhxEosk4Y8iGUGKQAB3CelPHbFuYwDkrDMrqUBLj1W8d9FSnlcCNEOGAC8J4RYLaWcYGtcSs2kWvxKRa0BHMS/k5QghGglhOhRhn1sx9Jt4m290fsIsL6UbTZjqRYK8CgQc+kvrS3Ti5fE8TiWFnYGkC2E6GRd/jBXtxJ4vuiF9f6AO5akmGm919C/pACtXS6uUspAKWWYtcLpJOv57bCes6f1Jup9l2z6D/DiJftpU9JxLjEDS5fNjks+oawAXrzk3kpb6/J6wAXrJ7HHsRQevNo5NADypJS/AB9h6QZSbnAq8SsVYm2N3wPcLCzDOQ9hSW6JZdjHBSwTfawF9gG7pJSllb59EXhCCLEfS+J6+SrrDMHSp70fy72DopbqU1i6OfZiuYdwte6L9wBP6w3NfVjm3d2HpYvnKJYuok2lxPgI/1b3LLIAeMRa5/8DLBe9TUDcJXG8BERbbyofBmwaBmqdLCcLyz2MIhOxdOvst/5tJlqXfwkMsZ5bBNf+9NUS2G59r8ZieV+UG5yqzqnUOUII16K5k4UQo7DcfL3aheO6xGFt8f8B/CClvPJCca1thwLRUsoXLlnWAEs/foS1JV+lrF1dS6SU80tbV6lZVItfqYtutw7lPAj0oPpaseOsLemDQCzwZxm2zQf6C+sDXEKIwcA2LKOOrkfS/xW4icvvDSg3CNXiVxRFqWNUi19RFKWOUYlfURSljlGJX1EUpY5RiV9RFKWOUYlfURSljlGJX1EUpY75f+o+axuKH7CrAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "df[\"age_bin\"] = df[\"age\"].round(1).apply(lambda x: math.floor(x/10)*10)#.groupby(\"age\")\n", "df[\"age_bin\"] = df[\"age_bin\"] - (df[\"age_bin\"]==80).astype(int)*10\n", "\n", "df_plot_metrics = df.groupby([\"age_bin\", \"modality\"]).apply(lambda x: pd.Series([x[\"age\"].mean(), x[\"y_hat\"].mean(), x[\"y_hat\"].std()],\n", " index=[\"age_mean\", \"y_hat_mean\", \"y_hat_std\"])).reset_index()\n", "\n", "\n", "for i, mod in enumerate([\"DF\", \"GM\", \"WM\", \"CSF\", \"DF+GM\"]):\n", " df_mod = df_plot_metrics.loc[df_plot_metrics.modality==mod]\n", " plt.errorbar(df_mod[\"age_mean\"], df_mod[\"y_hat_mean\"], yerr=df_mod[\"y_hat_std\"], \n", " label=mod, color=colors[i])\n", "plt.plot([df_plot_metrics[\"age_mean\"].min(), df_plot_metrics[\"age_mean\"].max()], [df_plot_metrics[\"age_mean\"].min(), \n", " df_plot_metrics[\"age_mean\"].max()], 'k--')\n", "plt.xlabel(\"Chronological Age [years]\")\n", "plt.ylabel(\"Predicted Age [years]\")\n", "plt.legend(ncol=5, bbox_to_anchor=(1, 1.15))" ] }, { "cell_type": "code", "execution_count": 39, "id": "c0962057-7ada-43db-9558-9ed693b975ab", "metadata": { "tags": [ "hide-input" ] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "########################## CSF ##########################\n", "\n", " Age bin Mean age [years] Mean predictions [years] Difference [years]\n", "0 20 26.09 39.43 -13.34\n", "1 30 34.56 42.96 -8.40\n", "2 40 45.74 49.52 -3.78\n", "3 50 56.12 64.64 -8.52\n", "4 60 64.79 63.75 1.04\n", "5 70 75.19 75.59 -0.40\n", "\n", "\n", "\n", "########################## DF ##########################\n", "\n", " Age bin Mean age [years] Mean predictions [years] Difference [years]\n", "6 20 26.09 35.23 -9.14\n", "7 30 34.56 39.52 -4.96\n", "8 40 45.74 47.86 -2.12\n", "9 50 56.12 64.74 -8.62\n", "10 60 64.79 60.42 4.37\n", "11 70 75.19 71.20 3.99\n", "\n", "\n", "\n", "########################## DF+GM ##########################\n", "\n", " Age bin Mean age [years] Mean predictions [years] Difference [years]\n", "12 20 26.09 34.39 -8.30\n", "13 30 34.56 38.62 -4.06\n", "14 40 45.74 47.17 -1.43\n", "15 50 56.12 64.06 -7.94\n", "16 60 64.79 59.82 4.97\n", "17 70 75.19 71.77 3.42\n", "\n", "\n", "\n", "########################## GM ##########################\n", "\n", " Age bin Mean age [years] Mean predictions [years] Difference [years]\n", "18 20 26.09 36.38 -10.29\n", "19 30 34.56 39.66 -5.10\n", "20 40 45.74 46.70 -0.96\n", "21 50 56.12 60.09 -3.97\n", "22 60 64.79 56.84 7.95\n", "23 70 75.19 67.77 7.42\n", "\n", "\n", "\n", "########################## WM ##########################\n", "\n", " Age bin Mean age [years] Mean predictions [years] Difference [years]\n", "24 20 26.09 31.78 -5.69\n", "25 30 34.56 31.93 2.63\n", "26 40 45.74 40.12 5.62\n", "27 50 56.12 43.59 12.53\n", "28 60 64.79 47.60 17.19\n", "29 70 75.19 51.07 24.12\n", "\n", "\n" ] } ], "source": [ "df_groups = df.groupby([\"modality\", \"age_bin\"]).mean()[[\"age\", \"y_hat\"]].round(2).reset_index()\n", "\n", "df_groups[\"Difference [years]\"] = (df_groups[\"age\"] - df_groups[\"y_hat\"]).round(2)\n", "df_groups.rename(columns={\"y_hat\": \"Mean predictions [years]\", \"age\": \"Mean age [years]\",\n", " \"age_bin\": \"Age bin\"}, inplace=True)\n", "\n", "for el in df_groups.modality.unique():\n", " print(f\"\\n########################## {el} ##########################\\n\")\n", " print(df_groups.loc[df_groups.modality==el].loc[:, df_groups.columns!=\"modality\"])\n", " print(\"\\n\")" ] } ], "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.9.12" } }, "nbformat": 4, "nbformat_minor": 5 }