{ "cells": [ { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "from matplotlib import pyplot as plt\n", "import seaborn as sns\n", "import glob" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "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", " \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", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
Meansampleconcentrationtimelabel
01187.830125uM-180min-11250001801
11141.036125uM-180min-11250001801
21144.027125uM-180min-11250001801
31116.973125uM-180min-11250001801
41080.741125uM-180min-11250001801
..................
69538.214DMSO-1001
69658.696DMSO-1001
69744.679DMSO-1001
69867.089DMSO-1001
69928.902DMSO-1001
\n", "

700 rows × 5 columns

\n", "
" ], "text/plain": [ " Mean sample concentration time label\n", "0 1187.830 125uM-180min-1 125000 180 1\n", "1 1141.036 125uM-180min-1 125000 180 1\n", "2 1144.027 125uM-180min-1 125000 180 1\n", "3 1116.973 125uM-180min-1 125000 180 1\n", "4 1080.741 125uM-180min-1 125000 180 1\n", ".. ... ... ... ... ...\n", "695 38.214 DMSO-1 0 0 1\n", "696 58.696 DMSO-1 0 0 1\n", "697 44.679 DMSO-1 0 0 1\n", "698 67.089 DMSO-1 0 0 1\n", "699 28.902 DMSO-1 0 0 1\n", "\n", "[700 rows x 5 columns]" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "['125uM-180min-1', '1uM-180min-1', '200nM-180min-1', '25uM-180min-1', '40nM-180min-1', '500uM-10min-9', '500uM-120min-1', '500uM-180min-1', '500uM-30min-1', '500uM-60min-2', '5uM-180min-1', '625uM-180min-1', '8nM-180min-1', 'DMSO-1']\n" ] } ], "source": [ "flist = glob.glob('quantification/*.csv')\n", "\n", "concs = [125000, 1000, 200, 25000, 40, 500000, 500000, 500000, 500000,500000, 5000, 625000, 8, 0]\n", "times = [180, 180, 180, 180, 180, 10, 120, 180, 30, 60, 180, 180, 180, 0]\n", "\n", "df=pd.DataFrame()\n", "sample_list =[]\n", "for i,file in enumerate(flist):\n", " df_tmp = pd.read_csv(file, index_col=0).drop('Slice',axis=1)\n", " df_tmp['sample'] = file[15:-4]\n", " df_tmp['concentration'] = concs[i]\n", " df_tmp['time'] = times[i]\n", " df = pd.concat([df,df_tmp])\n", " sample_list.append(file[15:-4])\n", "df = df.reset_index(drop=True)\n", "df['label'] = np.full(len(df),1)\n", "display(df)\n", "print(sample_list)" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARwAAAGTCAYAAAD6Eq3tAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABEJklEQVR4nO3de3xV1Z3//9c+l9zvkCsnkUC4JOESIHKxVlEEwSp0lCLUKiIOMzKtts63/dLx12lnOpbot9NWW9sZpgjRIpna6RCrBUSUilWMAS+FiARIyP1G7rdz2Wf9/kBOxSQQyD47OeHzfDx8PJKdfdZaJ+a8WXutvdbWlFIKIYQwgWW4GyCEuHpI4AghTCOBI4QwjQSOEMI0EjhCCNNI4AghTOO3wHnggQdISEhg2rRpfX724x//GE3TaGpq8h3bvHkzGRkZTJkyhb179/qOHz58mOnTp5ORkcHDDz+MzOILEbj8Fjj3338/e/bs6XO8srKSffv2kZaW5jtWUlJCQUEBx44dY8+ePWzcuBFd1wF46KGH2LJlC6WlpZSWlvZbphAiMPgtcG644Qbi4uL6HP/Wt77Fk08+iaZpvmOFhYWsXr2a4OBg0tPTycjIoKioiNraWtrb21mwYAGapnHfffexa9cufzVZCOFnNjMre+mllxg3bhwzZ8684Hh1dTXz58/3fe9wOKiursZut+NwOPocH8iWLVvYsmULAMePH2fq1KkGvwMhxKWUl5dfMFzyWaYFTnd3N48//jivvvpqn5/1Ny6jadqAxweyYcMGNmzYAEBubi7FxcVDaLEQ4krk5uYO+DPTAufUqVOUlZX5ejdVVVXMnj2boqIiHA4HlZWVvnOrqqpISUnB4XBQVVXV57gQIjCZNi0+ffp0GhoaKC8vp7y8HIfDwZEjR0hKSmL58uUUFBTgdDopKyujtLSUuXPnkpycTGRkJIcOHUIpxXPPPceKFSvMarIQwmB+C5w1a9awYMECPvnkExwOB1u3bh3w3OzsbFatWkVWVhZLly7lmWeewWq1AvCrX/2KBx98kIyMDCZOnMiyZcv81WQhhJ9po3V7ChnDEWJ4XOyzJ3caCyFMI4EjhDCNBI4QwjQSOEII00jgCCFMI4EjhDCNBI4QwjQSOEII05i6WlwIcXG7d++mrq6uz/HK1iC8BHFNTGefnyUlJQXMHfjSwxEiALSTQpuWdukTRzjp4QgxggzUU9m2bRsA69atM7M5hpMejhDCNBI4QgjTSOAIIUwjgSOEMI0EjhDCNBI4QgjTyLS4EAGg1jMejwoa7mYMmfRwhAgAFrxY8A53M4ZMejhCBIBEW8WnX31xWNsxVNLDEUKYRgJHCGEaCRwhhGkkcIQQppHAEUKYRgJHiBGus9dLmWsaHhX4k8oSOEKMcN1ORas3nl4VPtxNGTIJHCFGuIRoKxGWNs64s4e7KUMW+H00Ia4C19iP4VUWYMpwN2VIJHCEGMG8SlF0woUVN0EWWdoghPCjlk4v29/oIlabiTbcjTGAjOEIMYKNibTy8wdjibdWEWutH+7mDJn0cIQY4ew2jVhr43A3wxASOEKMYB09Xv77rW469Tismj7czRkyCRwhRjDdC+09Xuo9E/BiH+7mDJmM4QgxgsWEW3h0eRSTgw4zJejd4W7OkPktcB544AESEhKYNm2a79i3v/1tpk6dyowZM/ibv/kbWltbfT/bvHkzGRkZTJkyhb179/qOHz58mOnTp5ORkcHDDz+MUspfTRZixLJoCqsW+NPifguc+++/nz179lxwbPHixRw9epSPPvqIyZMns3nzZgBKSkooKCjg2LFj7Nmzh40bN6Lr565XH3roIbZs2UJpaSmlpaV9yhRCBA6/Bc4NN9xAXFzcBceWLFmCzXZu2Gj+/PlUVVUBUFhYyOrVqwkODiY9PZ2MjAyKioqora2lvb2dBQsWoGka9913H7t27fJXk4UQfjZsYzjPPvus78Ht1dXVpKam+n7mcDiorq6muroah8PR57gQV5sObyyt+tjhbsaQDcss1eOPP47NZuOee+4B6HdcRtO0AY8PZMuWLWzZsgWAxsbRcd+CEAANnlRcKmS4mzFkpvdw8vPzefnll9mxY4cvPBwOB5WVlb5zqqqqSElJweFw+C67Pnt8IBs2bKC4uJji4mLi4+P99yaEMInbozhY4iTZepp4SwXferaFhrbAvR/H1MDZs2cPTzzxBC+99BJhYWG+48uXL6egoACn00lZWRmlpaXMnTuX5ORkIiMjOXToEEopnnvuOVasWGFmk4UYVnWtOi8c7OKMJ5M6PZ2ls0KICQ/cu1n8dkm1Zs0aDhw4QFNTEw6Hg3/5l39h8+bNOJ1OFi9eDJwbOP6P//gPsrOzWbVqFVlZWdhsNp555hmsVisAv/rVr7j//vvp6elh2bJlvnEfIQLd7t27qauru+g5SsEYNR1b22E0pfNm0Xwq3z96ybuOk5KSRuRnRVOj9MaW3NxciouLh7sZQgxo27Zt1NTUEB0dPeA5Xiyc1m4CdzsKO5o9jPHqTWy4BnxNW1sbKSkprFu3zh/NvqSLffZkaYMQwyg6OpqFCxde9Jyb0aloicWlW8gY2wtcd9HzDxw4YFj7jCaBI0QASIsduEcTSAJ39EkIEXAkcIQQppHAEUKYRgJHCGEaCRwhAkBNm50zzcEAlDaG0NZjHeYWXRmZpRJiBFMKdh+PxesFl67xfnU4LT02Zo3rImdc13A377JJ4Agxwlk1BRaweDVC7V6+lNWAPTA7OBI4Qoxkmga3Tm2lujUIl66RPsY53E0aEgkcIQLAuBi58U8IIS6LBI4QI5jTo/GfbydSUhdKU1fgX5BI4AgxggVZFfOv6aC8OZg3Tw28qjxQBH5kCjGKaRrMHNfN2Ag3bk/g9w8kcIQIAMfrw+hyWRkvs1RCCH+7eVLbcDfBEBI4QgSAizysJKAE/kWhECJgSOAIEQA+qgnj3TMRw92MIZPAESIAKHXuv0AnYzhCBID6jiA6XVagc7ibMiQSOEIMk+bmZnp6egb1lAVFOGFYOXCg/ZLntra24nSOzOlzuaQSIgAE00Uv0XSQNNxNGRLp4QgxTOLi4ujq6rroc6ncusZfasMYH9fLn8uiiA31cP2EjouWe+DAAcLDww1urTEkcIQYwTqcVooqIqluDcLttVwybEY6CRwhRrC4MA9/f10dbl3DOwpmqWQMR4gAoHvB4w38240lcIQIAH86Fc2e47HD3Ywhk0sqIQLA9ent6Ep6OEIIExRXRvDW6ajhbsaQSQ9HiBHOqyAsSCc1dmTezHc5pIcjxAh3tsvG2+VRRIXow92UIZPAEWKEi4/wMCO5i+LKwF8tLpdUQgSAjLG9uPTA7x9I4AgRAFJjR8eD8CRwhLiK7d69m7q6uj7HK1uDsCgP42K9fX6WlJTEsmXLrqg+v/XRHnjgARISEpg2bZrvWHNzM4sXL2bSpEksXryYlpYW3882b95MRkYGU6ZMYe/evb7jhw8fZvr06WRkZPDwww+jRsMuREJcpuMNoXxYE2ZafWctmTRrEw0v1289nPvvv5+vf/3r3Hfffb5jeXl5LFq0iE2bNpGXl0deXh5PPPEEJSUlFBQUcOzYMWpqarjllls4ceIEVquVhx56iC1btjB//nxuu+029uzZc8XpKkSg6nJa6HFbDS93oM/Ss89uQ9Ng3bp1htbntx7ODTfcQFxc3AXHCgsLWbt2LQBr165l165dvuOrV68mODiY9PR0MjIyKCoqora2lvb2dhYsWICmadx3332+1whxNZmT2sX1Ey69+ZZR/PWUCFOHvevr60lOTgYgOTmZhoYGAKqrq0lNTfWd53A4qK6uprq6GofD0ef4QLZs2UJubi65ubk0Njb66V0IIa7UiJhn629cRtO0AY8PZMOGDRQXF1NcXEx8fLyhbRRCDJ2pgZOYmEhtbS0AtbW1JCQkAOd6LpWVlb7zqqqqSElJweFwUFVV1ee4EFeTbpeFXnfgL9wEkwNn+fLl5OfnA5Cfn8+KFSt8xwsKCnA6nZSVlVFaWsrcuXNJTk4mMjKSQ4cOoZTiueee871GiKvFq5/EsL80Brce+KHjt8BZs2YNCxYs4JNPPsHhcLB161Y2bdrEvn37mDRpEvv27WPTpk0AZGdns2rVKrKysli6dCnPPPMMVuu5Eflf/epXPPjgg2RkZDBx4kSZoRJXnfPPFX/5mOyHM6CdO3f2e3z//v39Hn/sscd47LHH+hzPzc3l6NGjhrZNiEASFaKTldiN3dr3JrxAMyIGjYUQF/enU9FUtIQMdzOGTAJHiAAwNaEbPfA7OLKWSohAEB2q0+sO/P6BBI4Qw6itrW1Qj/qt6x2LVwum7eTAN75+tkx5EJ4Q4gJJSYN7bG+bPhanKwivFkJ4eOslzw8PDx902WaTwBFimAzmFo/GNp1/2tFGQng5VjyGL6b8vM5eL0/+bwdxehwRljbDyw/8i0IhRrH4aCt590YTY2kk1NLp9/qC7Rozx9upcGdy2j3D8PKlhyPECDcm0kq7dywu5f9pcbtV464FYdSVfIQFHcgxtHwJHCECQLy1Eh0bkG1KfWGWDr+UK4EjRACItLZc+qQAIIEjhOijWU/EisfwciVwhAgQurLS61aE2P2/arxRdxBEr+HlyiyVECPc0y93UObK5qR7Fr/e5/+ZKoApQYdJDzpmeLkSOEKMcDdND8ZNMF5l4a755j25wR/kkkqIEW76NUEkWivwYCM5zpxZKn+RwBEiAERbm4a7CYaQSyohhGmkhyNEAOjyRuIdBR/XwH8HQlwFaj0TTVna4G8SOEIEgIn2Dzn3lLbMYW7J0EjgCBEA6vTxuFXwcDdjyCRwhAgAVjx4TZzjcalgNIzfRFkCR4gAEG+tRKEBN5hSX6lrDkGa8UsbJHCECACn3DmmDhpPtL+PRfMCMw0tVwJHiBGs163YebCLWEsdFk3HrEHjEEuPX8qVwBFiBNN1xdl2Lw16mqn34SgF5e5pVDd7GBdnXL0SOEKMILt376auru6CY2OAlo4gFMFs2/Z2n9ckJSUNakP2wWjv9vLOJ0669Eg6vLF09ihDyj1PAkeIAOAKTsejhQHv+7We+jad1z7qxeuZQailiynj7IaWL4EjxAhysZ6KUgpNy/Fr/ZOS7fy/tbH8v/96FytuZNBYiKuUpvl/p7/zmnSHTIsLIcwxNfi9T7/KNbRcCRwhhM/ZDp3/3NtJuDeGUK3b8PIlcIQQPiF2jbR4G0cas/1yL44EjhDCJzzEwtduDKfz1GEsGH+nsez4J4ToI1jrxaVCaGrX+bDcZVi50sMR4irR302Fn+dWQVS4MwlqbqMu4np+vPMYZ/VxzA7ed9FLrMHefDgsPZyf/vSnZGdnM23aNNasWUNvby/Nzc0sXryYSZMmsXjxYlpa/vpo082bN5ORkcGUKVPYu3fvcDRZiIBXV1dHTU0NXV1dA/7X0e2iwxNBY/gtBOvNpLjfItv7e/SepgFfU1NTc8kgO8/0Hk51dTVPP/00JSUlhIaGsmrVKgoKCigpKWHRokVs2rSJvLw88vLyeOKJJygpKaGgoIBjx45RU1PDLbfcwokTJ7BarWY3XYiAFx0dzcKFCy9xlov6Djt2azRxYZfeDuPAgQODrn9Yejgej4eenh48Hg/d3d2kpKRQWFjI2rVrAVi7di27du0CoLCwkNWrVxMcHEx6ejoZGRkUFRUNR7OFuGokRrqJCzP+2eKmB864ceP4P//n/5CWlkZycjLR0dEsWbKE+vp6kpOTAUhOTqahoQE41yNKTU31vd7hcFBdXd1v2Vu2bCE3N5fc3FwaGxv9/2aEEJfF9MBpaWmhsLCQsrIy3/Xkb37zmwHPV6rvatWBbvHesGEDxcXFFBcXEx8fb1ibhRDGMD1wXnvtNdLT04mPj8dut3PnnXfy9ttvk5iYSG1tLQC1tbUkJCQA53o0lZWVvtdXVVWRkpJidrOFEAYwPXDS0tI4dOgQ3d3dKKXYv38/mZmZLF++nPz8fADy8/NZsWIFAMuXL6egoACn00lZWRmlpaXMnTvX7GYLcVVx6xoe4/dQN3+Wat68eaxcuZLZs2djs9mYNWsWGzZsoLOzk1WrVrF161bS0tJ48cUXAcjOzmbVqlVkZWVhs9l45plnZIZKCD97+Vgs4cFelkxpNbRcTfU3SDIK5ObmUlxcPNzNEGLE2LZtG11dXYOYFofGThs2iyI2TL/kuQcOHCA8PJx169YBF//syZ3GQog+QmxeLH4YcJG1VEKIPl49Ectbp6MML1d6OEIIn+ZuG2+djmLq2C7Cgo0fbZEejhDCR0NhsypOtYRypsX4Z5lL4AghfGLDdG7LbCElys3Y8FGwtEEIMfLpXuh2WThwMoo3TkYbVq4EjhDiAi8fi8Vu9XLoTCQ2iyItxmlY2TJoLIS4QGqMk5r2IKYldTF/fAc2A7slEjhCiAvMHNcNGoTYlKFhAxI4Qlw1mpub6enpueiGWS7CadCyiFFHseGmtqTtkuW2trbidA7uskvGcIQQPjp2eoihXptOozbV8PIH3cN5++23KS8vx+P561TZfffdZ3iDhBD+ERcXN8i1VG5q220EWSIYE3Gpc/+6lmowBhU49957L6dOnSInJ8e3UlvTNAkcIUapQ2eiCA8yfrX4oAKnuLiYkpISUx+mLoQYPrdntaBpw7S0Ydq0aYN+DIQQIrC9Xx3OayeiDZ+hgkH2cJqamsjKymLu3LkEB/91fcVLL71kfIuEEMNqTJibhg47p88GM2GMcTf9wSAD5wc/+IGhlQohRq60WBcnm0Jp6AgansC58cYbDa1UCDGyxYZ5CLEZv6nxoK7SDh06xLXXXktERARBQUFYrVaioozfnEcIMTJ0uyz0uI0fxBlUiV//+tfZuXMnkyZNoqenh1//+td8/etfN7wxQoiRocdtwekxPnAGfeNfRkYGuq5jtVpZt24d1113neGNEUKMDI4YF8FW4y+pBhU4YWFhuFwucnJy+M53vkNycjJdXV2GN0YIMTJMTejxS7mD6jM9//zzeL1efvGLXxAeHk5lZSX/8z//45cGCSGG39tlkbxfNbjlCpdjUD2ca665hp6eHmpra/n+979veCOEECNLaJB3+Gap/vCHP5CTk8PSpUsB+OCDD1i+fLnhjRFCDC+PF6pagxgf14sj2th7cGCQgfODH/yAoqIiYmJiAMjJyaG8vNzwxgghhteZ5hB+/9EYXvskhrfLh+m5VDabjeho4zZSFkKMTBPH9rJuXj1KgW04F2++8MIL6LpOaWkp3/jGN2RaXIhRKjLYy/4TMbxVbnwnY1CB8/Of/5xjx44RHBzMmjVriIqK4mc/+5nhjRFCjAzXT2hnblqH4eUO+j6cxx9/nMcff9zwBgghRpYDJ6OYNLaX+AiX4WVfNHAuNRMl21MIMfr0uK2U1IfS2mslO8nYGwAvGjjvvPMOqamprFmzhnnz5qGU8YNIQoiRZVlmC4fORODxGr/D50UDp66ujn379rFz505eeOEFvvSlL7FmzRqys7MNb4gQYuSYf02nX8q9aOBYrVaWLl3K0qVLcTqd7Ny5k4ULF/LP//zPfOMb3/BLg4QQ/tPW1nbR51L5zut0oqGIiggZVJmGPbXB6XTyyiuvsHPnTsrLy3n44Ye58847B1W4EGLkSEpKGvS5J/QbsKouksOPX/Lc8PDwQZd90cBZu3YtR48eZdmyZXz/+99n2rRpg2utEGLEWbZs2aDP7d76OyzorFu3ztA2XPQ+nOeff54TJ07w1FNPcd111xEVFUVUVBSRkZFD2vGvtbWVlStXMnXqVDIzM3nnnXdobm5m8eLFTJo0icWLF9PS0uI7f/PmzWRkZDBlyhT27t17xfUKIQZHx4aO1fByLxo4Xq+Xjo4OOjo6aG9v9/13/vsr9cgjj7B06VKOHz/Ohx9+SGZmJnl5eSxatIjS0lIWLVpEXl4eACUlJRQUFHDs2DH27NnDxo0b0XX9iusWwggflLn4/aHu4W6G31S6p1DrmWh4uaY/W7y9vZ0333yT9evXAxAUFERMTAyFhYWsXbsWOHcpt2vXLgAKCwtZvXo1wcHBpKenk5GRQVFRkdnNFuICTo+i1zV6bxOZGvQuE+0fGl6u6YFz+vRp4uPjWbduHbNmzeLBBx+kq6uL+vp6kpOTAUhOTqahoQGA6upqUlNTfa93OBxUV1f3W/aWLVvIzc0lNzeXxsZG/78ZcdWaNymYr95g/AZVI4VFU8P35E0jeTwejhw5wkMPPcT7779PeHi47/KpP/3dbDjQI4c3bNhAcXExxcXFxMfHG9ZmIYQxTA8ch8OBw+Fg3rx5AKxcuZIjR46QmJhIbW0tALW1tSQkJPjOr6ys9L2+qqqKlJQUs5stxKh3tkOnsKgbr9d/l4qmB05SUhKpqal88sknAOzfv5+srCyWL19Ofn4+APn5+axYsQI4t56roKAAp9NJWVkZpaWlzJ071+xmCzHqne3wUnzSxRtHe3F57XTosfz0D+30uo0LoEE/JsZIP//5z7nnnntwuVxMmDCBbdu24fV6WbVqFVu3biUtLY0XX3wRgOzsbFatWkVWVhY2m41nnnkGq9X46TohrnaTU+x87cZwfv5KB0H6XGyam6lhFqwGLqnS1ChdkZmbm0txcfFwN0OIgKOU4lfPvogFL3+/fvVlv/5in71h6eEIIUYuTdOo9mQQpBm/iboEjhCij2hLEzbNbXi5EjhCiD7ibf3f6zZUps9SCSGuXhI4QlyBA0d7+Y+9xm8yPtrJJZUQVyA51opu/JNwRz0JHCGuwJRxdqaMsw93MwzX1u3lzx87afGMxWYxflcGCRwhhE99q86fjjnxeiYSbDH2iQ0gYzhCiM+YnGLniftiyAx+l4lBHxlevgSOEKKPATZkGDK5pBJC9PGJK5cgeg0vVwJHCNFHrKUOm2byo36FEFenRj2NIE16OEIIE0y0v49F8wIzDS1XAkcI0UeIH6bEQWaphBAmkh6OEJewe/du6urqLjhW1+LGrUWQGtN3z5ikpKTLesrl1UR6OEJcgWbSqbflDnczAo70cIS4hP56K88+uw1FteHP3h4plPLPzX8SOEJcAU0DjVG5HTgAx11z/TItLpdUQog+UmynSbBVGF6u9HCEEH1EW5v8Uq70cIQQppHAEUKYRgJHCGEaCRwhhGkkcIQQppHAEaOC26OobzV+029hLJkWFwGnv7VNJ9uSaLLnMD9sT5/z/bG2qU0fS6c32tAyrwYSOGJUCOs5SoqzHMJCTKmvV4XTpWJMqWs0kcARAae/3sq2bdsATFvblGg7QyJngBmm1DdaSOAIcRXr7/IU4MO2aVhVty/IP2sol6gSOEKIPmLVaSy4DS9XAkeIq5jZG4XJtLgQwjQSOEII0wxb4Oi6zqxZs7j99tsBaG5uZvHixUyaNInFixfT0tLiO3fz5s1kZGQwZcoU9u7dO1xNFkIM0bAFzlNPPUVmZqbv+7y8PBYtWkRpaSmLFi0iLy8PgJKSEgoKCjh27Bh79uxh48aN6LrcUSqGh0dX/O7tbrq8kfR6Q4e7OQFnWAKnqqqKV155hQcffNB3rLCwkLVr1wKwdu1adu3a5Tu+evVqgoODSU9PJyMjg6KiouFothB4dCipclPtmcTHrgXD3ZyAMyyB881vfpMnn3wSi+Wv1dfX15OcnAxAcnIyDQ0NAFRXV5Oamuo7z+FwUF1d3W+5W7ZsITc3l9zcXBobG/34DsTVSgFRoRaitEbircZvwTnamR44L7/8MgkJCcyZM2dQ5yvVd6NqbYDt5Dds2EBxcTHFxcXEx8cPqZ1C9MeiQUSIhpMw2r1jh7s5Acf0+3D+/Oc/89JLL/HHP/6R3t5e2tvb+drXvkZiYiK1tbUkJydTW1tLQkICcK5HU1lZ6Xt9VVUVKSkpZjdbCABcHsXJOg/h6IRZ2oe7OQHH9B7O5s2bqaqqory8nIKCAm6++WZ+85vfsHz5cvLz8wHIz89nxYoVACxfvpyCggKcTidlZWWUlpYyd+5cs5stBAChQRrXZwZjw4VSclfJ5Roxdxpv2rSJVatWsXXrVtLS0njxxRcByM7OZtWqVWRlZWGz2XjmmWewWq3D3FoxkpyocdPrDSPE0u33umxWjdtzQyl5X0cpPzwpbpQb1sBZuHAhCxcuBGDMmDHs37+/3/Mee+wxHnvsMRNbJgLJb//cTZd7Knat73O+xcgifUIR8L7zN1GMtVbRq8JNq7NLRdPmlYmJyzViLqmEuFJBNo04WwNxtgYg+4rLGWirhs/SlZWT7hwsbeWEWhvYtu31QZXtj10HA5EEjhg1Kt2TqW3RSY69sjG+uro6ampqiI4eeOtQHRsuTcNiT8BpTaCr68wly21ra7ui9oxGEjgi4Hm9iirXRM56U2hou/LAAYiOjvaNKw7k9NlgkqMS8CqN6taljItxER7kHfD8AwcOXHF7RhsJHBHwelyKWj2DBOsZUuLS/FpXp9PCy8fiyBjbQ1OXnW63lYUT25ia2OPXekcLGTQWAS88xEJuyF4a9VQ+qfb4ta6IYC/Lprbg9FiYf00Ha3PrJWwugwSOGBVqPBPRNMX1mcF+r8vt1YgN87C/NIZPGsP8Xt9oIoEjRoUwSxvRWpPf6+l1a+w7EYPuhaRIF9OTu/xe52giYzhiVNBVEN0qyu/1hNgV9+U20OWy0txtwyr/ZF8WCRwxKoy11TDWVgNM9HtdsWE60aE6CZHGP9VgtJN8FgHP5fnrFiYHS3r5uMpFYZF/11UdKo/ktx+c257CrWt4Bp4VF58hPRwR0NwexaPPthCrptOqx1PyXg9zJwVxosbDHdcqLAPsndSf5uZmenp6BrxvppN46rVphNKMTdUSy1n+cMBKg5aJHSdp6p1+X9fa2orTKeu8QAJHBDi7TWPdogj+/EYlYZYOvrv2Jr/VZcVFMB2AwkkUXdpY2kklRR0hkosviRDnSOCIEWswa5sAlIKTXXOIcX7Q76Np+9Pf2qa4uDi6uroueadxZWsQHl2jojWY6lYPd8wYR4h94E3hDhw4QHi4eQtLRzIJHDFiDWZtE5zbZ9jubaPTPp6yLhceLYQU9cGA5w9lbVOvW+N/PxpDZmI3LT02kqPc2Cx9t8EV/ZPAESPaYNY2VbcF4akLpaIlhInjOtE0mO2IGfD8oaxtCrErHphXT6fLSku3jUy5y/iySOCIgOf0aHS6rOSmdjBznP93/YsI9vJxfRhlzSESOJdJpsVFwIsN9XBNrJO/1Jo3ThId4iEpUmaeLpcEjgh4B09HU9MWRHKUi1dKYk2ps7XXRl2H/9dtjTYSOCLg3Tq1heZuG0pBVlIX1W1B6H6+EW9uWiercvy/dmu0kcARAU/3anxxQjs9bgtHqiL43YdjqGz1b++jrDmYd89E+LWO0UgGjUVAUwq2FyWQldSNW9f44oQOgm1tRIfqfqvzw+owetwWOpxWul0WiioiyEzsIVHWVl2SBI4IaJoGy6c1o3tBI8SUBZWlTaGMDXcTavfy3HsJJEW56PXIxcJgSOCIgOeIcVHZEoTTgA99W1vbJe/TGQscb76FEE8NqfYGws82UXYWyi5SptxpfI4EjhixLrWYEsBDEG7C0AmiRUvn1bqPsOHEwsCjxgMtpkxKShp028KaT9Njv4aE8DPAxcMkPDz8ssoezSRwREBrYjIt2gRS1BFchFOqLWWs+oRESi67rME8N6qr10t5g07Hq3+gV1Vz+533ceS0i8QYKznpQVfyFq4qEjhixBrMYkqPF373gZ246FmEeqzMGtdBVEgiwbaEAV8zlMWUR067+e8/d6G5ZxNpaeWxHW3MnxJEsF2eMz4YEjgioNkskJnUQ3uvldNnQ1gypdWv9V2bYedkbRBnTtagofi/d0YyMcnu1zpHEwkcEfBmpnTT1GUjxOb/Vdu6F+pbvaBpdHmjJWwuk8zliVGhqcvOicZQv9cTHmJh011RWNHxIGM2l0t6OGJUmJrQw9QEc1Zu7/hTFxY8JFtPA5NMqXO0kB6OCHgVLUF0OS00ddr4/UdjqGjxf8+jyxtNnT7e7/WMNhI4IqB5vPDS0TEcPB3Fbz8Yy5hP7wD2p3tuDCfZXkai9Yxf6xmN5JJKBDSbBa6f0EaPy8JX5zQS48c1VJ8VYWkjwnLlW5VeraSHIwKezQLtThuHK2X19kgnPRwR8KYldxMZolNS5/9ZqvPa9DGmPFp4tDG9h1NZWclNN91EZmYm2dnZPPXUU8C5dTOLFy9m0qRJLF68mJaWFt9rNm/eTEZGBlOmTGHv3r1mN1kEgDFhbibH95pWX7eKpE0fa1p9o4XpPRybzca///u/M3v2bDo6OpgzZw6LFy9m+/btLFq0iE2bNpGXl0deXh5PPPEEJSUlFBQUcOzYMWpqarjllls4ceIEVqvV7KaLYTCY1duntJvQ3C10264hjY+IpBaNgW8CNGL1drKtnGRbOTBtSOVcbUwPnOTkZJKTkwGIjIwkMzOT6upqCgsLfX9Ya9euZeHChTzxxBMUFhayevVqgoODSU9PJyMjg6KiIhYsWGB204XJBrvCOtbdSrM7iFD3acqCb2Ra0FtEWZsHPN+I1du93lBcyrxLuNFiWMdwysvLef/995k3bx719fW+IEpOTqahoQGA6upq5s+f73uNw+Ggurp6WNorzDWY1dvnPflf+2jzjuGacXYeuX2FH1t1Tr0+nhY90e/1jDbDFjidnZ3cdddd/OxnPyMqauDBN6X6do21AR5Qv2XLFrZs2QJAY2OjMQ0VI5ruVWz+n3Y0gtCxMnO8OcsN0mwf47CdANJNqW+0GJZpcbfbzV133cU999zDnXfeCUBiYiK1tbUA1NbWkpBwbnsBh8NBZWWl77VVVVWkpPT/HOcNGzZQXFxMcXEx8fHxfn4XYiSwaDD9GjsaOl5l56ZpIabUq2lg1cy552c0Mb2Ho5Ri/fr1ZGZm8uijj/qOL1++nPz8fDZt2kR+fj4rVqzwHf/qV7/Ko48+Sk1NDaWlpcydO9fsZosRStM0VswNo+ovZwH/TCTs3r2burq6C46dafTgssSybdu2PucnJSVd1uXg1cT0wPnzn//M888/z/Tp08nJyQHgRz/6EZs2bWLVqlVs3bqVtLQ0XnzxRQCys7NZtWoVWVlZ2Gw2nnnmGZmhEn10eaOp96SZVp8rOI0uyzjgPdPqHA001d8gySiQm5tLcXHxcDfjqtBfD6CyNZg6+wJmhx3sc+nhjx7A+Z7GunXrDC1XXL6LffZkaYPwC83ZRITrEyzIOIf4K1naIIasv97KuR5HrSk9jh1/6qJVH0uQ1vdJDGJkkcARhvuk2o1LhWDBY0p9Hb2Ks3oKTbrDlPrElZNLKmG47W90Ue7Oorh3Kb0u/w8R/v2tEaTbjzIz+IDf6xJDI4EjDPeN2yJIsx5lUlAxIUHmPD7FonkJsXSbUpe4chI4wnA//UMHlXoWZa4ZON3mTYLqykphUTc9JvSqxJWRwBGG6Ojx0tB2bkbqm3dEMs5WSor9JEGfjhJ6dMWbJb2GB5DuVfz0D+3Uea6hzJ3NO5+4aOv27xaj4srJoLEwxB+KezhZ4yE7zc7SWSE068m0e8dS06zT3uPlgzIXB0tcXBNv45p44/7sNA3io6ycqbXgVqE8eW+MYWUL40ngCEN8eW4oda062/Z30d6to5SGnV5++GI712cG096t+MXfxmKxGDumY9E0vnZjONtOl5FCGZBlaPnCWBI4YtD6u6P48yYAJ0/m4OrqoS14ItODD+A81U6opsjPv3j5g70Dub92VDT00mMbJ2ubRjgJHDFodXV11NTUEB0d3ednXixYODd2ksyfafbY8FjCQdXSg0IBGqAA9Zlzz2trG9oTEPTgJDpt2cCbQypH+JcEjrgs0dHRLFy48IJjHh22vJPEnNROgmwKt0cjKcTDx/VhjB83ljMtITR320mMcOPSNRo67aye1XRBGZfaRvSzLt5bmTj4NyNMJ7NUYshsVlg8pRW7VVHZEsSpsyGcbg6hx23hlY/HEGr3Mj62h/LmYBwxTq5Pbx/uJothIj0cMWjNzc309PT06Y20k0IoZwELvYynWZtKdPsBnFo2iaqe98sno7ASy3FOHKtDw8thQojmr1vFtra24nTKWqjRTgJHDIkCzmjXEaWqaNdSSVQfEqmqqNLmobDiooMIVU8U1XQxllOWxcSo0+iEEK1kb+qrjQSOGLS4uDi6urr6jOF8QT+LUqE0drbyx49nkZXYhULjL7XhWMMmkx7Xi1tPQe+yszajnpjQ89uA/rWcAwcODPnRLWLkk8AZRq1dXmLCA2sY7VLPiUojmsNVXyKh61ViLXa6PFM50dZFszaRbPV7Pjjb/2WTEc+KEiOfBM4waWzT+acdbfzfv4kkI9k+3M0ZlEs9y6lNH0uc5SyZ3rf5WN1CjPMIITY3SdZqUqkgxBKMVev/T86IZ0WJkU8CZ5jER1t55PYI0hMD53/BxaajG9t1/uk3bdw1P5TbJwbxdMH7eK3xTJk8iUUzpvH2cSdHKzz865q+9/CIq0fg/LWPQqFBFgZ4xFbAiY+y8sS90Tz9SifNnV50Zcei6ZysdXOixk12mp2v3Rg23M0UwyywBhD8pL3by4kat6l1NrTp5P2+nZO15uyKZ4a4SCv/31eiONPoQaHh8oawYEow3/mbKFYuCGNySmBcOgr/kcAB3v7EyX/t6+Q/93ZS2WROACREW/nhmmhTP4Sdvf7ftsFm1Xjk9khSbKewajpLckKIDrMQZBslXTkxJHJJBSyZGcKCKUG88GY3/ti6affu3VTUNnPMeT2pto8JtvQC4GotQ8NLZGwSZ/VkEq1n0LShLzbsb3FjVauNiqDbmB3yWp+d8fxR39lmL10hkJ//HBbtwqCTxZRXLwkcwGLRiA6z8tDSSL/VYcNDnLWWFj0Jlx5Km3csid4eXNYYWtwJNOnjiLdWYzV443FdWelREXhd9cR7Xic4tJsK91TirHVEWFoNreuzgm2KENrR/BLhIlBdlQ/CG2ibhZoWcGpRtIXMYYL9Q6KszYCx/yJ7vQqvgpeLe3j3/RNYNZ1F12Vz03T/PBP7jaO9/O7tbiyeDuJsdfztXXP4zZ+6uTUnhFkTgvxSp7i6yYPwPqfTHYxb/fXDpisbxT1LaLJk0WydQqy1jmCtG6VAV8b+iiwWjaMVbty6wq45qdfTSYr136OLv5gZzMNfimSc7QQtngReeLObxBgLsyYEcbDESdXZ0TNoLUa+UX1Jdb4nc8o1gwhLK1bNzRhLLR85byBEi8dFKKm2jwnROkmwldPgiSfKdZS0iDYqPVPxKgtNXgfWulfZtm0bPd5w3MpOs3cc4+3H+vR8LrZBVYMnFS9WgrQeOvRYnCqMJu8kxnQf4K29XbyqouhREWh4SbBV9VvG5fa0Pq5y80pxD1FhGvX6NYRaOpkyLoHIUAvP7u+k6qyO7g3GMWZU/xmIEWRU/6Wd3zCKmAl0KhsN2gyqOEsQrbixYKOFk3o24TTTQjoR3hO4LbEU9cwnmHbGqJOEYaOyK4F2bRxdJBClqnBqkVQ31vRbX1VNIyHRKYTSCoDi3OxMFxouLYSzTCdOncROCzGuRtqDp3G8t4t27RrGqFJs9BLu7LqgXDchdLad7VPfQAGnlIamKbq8UbR40inTxxLsrKfdlsqh908RqnXQ5o1nUtBhThf1cvTdKKx4DB9MFuLzRnXgNDc3o5QiSf0FgFhVRgcJ6Nhp19JwYyGUVlpJIUl9QLM9nShVS4yqwE4XNVoOQfQQqlqw4iJRfUQUVQSpXlqVorm5uU99HUFTqdLmMoZSxqhT1GizsaDjJoQ09Q5xnKKDZFxE0Bo8kRR1mBDVQrhqxIpOOI00k06zNpHx6k+cZQqtWir20ArCmo9cUF9/O/A1MJWzWgYaEKPKUVo3EdShqx5iXB9x1jKfLPU27doCKrodnNUmE0IrIbSSqt7zlTPUHfiE6M+oDpzPa9Km0Ek8TqKJU6dwa6HEqlOE0cxZLYNQmujUElHYiFIVjFWltGmp6NhwEUarNh6rchPEmQHriHQeJSG0mQrtC+iEYFfduLRzj73tJoYWbSI6drq0JNK8f6JC+yKxqpweLRbQUAp6icWuuqjVcugghTT1Ds6eSvjcuPL5QD2vgyS82AhXDYTShgU3Z5nEGHWcs8GZ9GAhQR3jmLaSeHWcYNqJVyW4tFAS1dELylb9BKoQQzWqAycuLu7cJdWnQtRZXFoIIbRh49zlQ7l2E4nqGHHqJN1uG3bNSZc9HbcWTifxxFJBCG1Ylc543vJNW2uaRlxcXJ/6emtqCKGTyWovDUwFNLoZQzBtNGrZeLESp04Rqpqp0L5ItPMvRAW1fLrPr5sWLR034cSpk+iEkKUKAdDpJS4u5aLvt5conFokXcRjVW5iOYPCQoc2jl4iGatO0KnFk6iOEkYTjVomGl7chKFG95+CGCFG9V/Z51cfu3XA20uvCkNXcVg0GKNqsRCCh2i6VTQRruPEhTRgBSK0KjpVHG3edMbZS4myBQPBQP+rmz//fTqVn35VQY8eRo83ggZvGkqLwe2NIsxZhscSQqc1A68WTJSlniDVSYs3nAnB5Z+uswofsL7z+9PceONCPmkIZUakC5fHQlJUO//z0QRCQlP5UnwPuoqgpq2Xk42TGRfhZnxsCG+cmsb8tA4yk7oJtfcAuReULfvTCH8Y1YFzsQHPv5xxcbbDy8ESJ/HRFrp6FW01HdiDkpg4cSK1LV6+kBlMt1PhGGMldWwycRELr7i+37/TzdFKNx0tOjlZwXgVvHesDqcKZdbUKN476cIVnsiUcTYcY2wsnLZuUO+xra2NfX8q4rh2B9FU0EsMcaqCsfyJE9oyKmq6aNPGE999AKfNQWtXD682ZONQRVSd6qL+VC+9xOJFQwPiKPOVK4EjjDaqA2cgn53dcQCdHdFYvHbszgSCnEewVH7AOKD80zFUS1ISM4d4639rbR0pWIm3BXH642TavXF0e8cS3fsXOk61M1ELIaK3BecpnZ7uRJh26fo+2+OZp/4IKNq9Y/jEdS1OSxJJlgo0peOwVFIS+gUiXCfosqaSaXuHFj2JZpVJmKWNs7qDRGsZNs1NuL0BGPr+NB5dUdOsk2bgUzZF4JO/BiDC0oauWem2zyVGnQG6Lvmay6VpYEXHqvWQYjlNkirjdHsi0Zwg3BoFdFx2mQP1qHbv3k1RZTQ6Nmr0TEL0DkLdZ+gImkyM1kyNJwM3IUyyFxNm6eQLjjpu/9LST1993ZW/yc94v8zFr/d18dT6WEKCZOGmOOeqDJyBPqj36wqbdZVp9X36U8PrA0izfwKAw1ZKm3cs3b1N9LgbGBfSgUWDIK3Xd67VYvzqljkTgkhdbZOwERe4KgNnIDbr6PhwjISb9SwWza9LNkRgCpi1VHv27GHKlClkZGSQl5c33M0RQlyBgAgcXdf5h3/4B3bv3k1JSQk7d+6kpKRkuJslhLhMARE4RUVFZGRkMGHCBIKCgli9ejWFhYXD3SwhxGUKiMCprq4mNTXV973D4aC6uu9TG7ds2UJubi65ubk0Njaa2UQhxCAEROD0t0eY1s/jDjZs2EBxcTHFxcXEx8eb0TQhxGUIiMBxOBxUVlb6vq+qqiIl5eLrioQQI09ABM61115LaWkpZWVluFwuCgoKWL58+XA3SwhxmQLiPhybzcYvfvELbr31VnRd54EHHiA7O3u4myWEuEwBETgAt912G7fddttwN0MIMQQBcUklhBgdJHCEEKaRwBFCmEYCRwhhmlH75M2xY8cyfvz4y35dY2OjqTcNmlnfaH5vUt/Iqa+8vJympqb+f6jEBebMmTNq6xvN703qC4z65JJKCGEaCRwhhGkkcD5nw4YNo7a+0fzepL7AqG/UDhoLIUYe6eEIIUwjgSOEMI0EjhDCNAGzWtxo5/fVSUlJ4ZZbbuGFF17g7bffJjMzkw0bNmC324e7iUOmlKKoqIjq6mo0TSMlJYW5c+f2u1uiGLzm5mY0TSM2Nna4mxJwrtpB43vuuQePx0N3dzcxMTF0dnZy5513sn//fpRS5OfnG1pfW1sbmzdvZteuXb79lhMSElixYgWbNm0iJibG0PpeffVVNm7cyKRJkxg3bhxwbqfEkydP8stf/pIlS5YYWt9w2Lt3L7t27bogUFesWMHSpUsv/eLLVFFRwXe+8x32799PTEwMSina29u5+eabycvLu6K72keatrY29uzZc8Hv89ZbbzX0b/OqDZwZM2bw0Ucf4fF4GDduHDU1NVitVpRSzJw5k48++sjQ+m699VZuvvlm1q5d63tmd11dHfn5+bz22mvs27fP0PoyMzPZvXt3nw9CWVkZt912Gx9//LGh9YG5AfDNb36TEydOcN999+FwOIBzgfrcc88xadIknnrqKUPrW7BgAd/85jdZuXIlVuu5B/zpus6LL77Iz372Mw4dOmRofWBOAJz33HPP8S//8i8sWbLkgn+g9u3bx/e//33uu+8+Yyoy/N7lAJGdna2cTqdqbm5WERER6uzZs0oppXp6etTUqVMNr2/y5MlX9LMrlZGRodxud5/jTqdTTZw40fD6HnnkEbVs2TK1c+dOdfDgQXXw4EG1c+dOtWzZMvXwww8bXt+kSZP6Pe71elVGRobh9V2sTH/Ul5+fryZMmKD+/u//Xv3whz9UP/zhD9Xf/d3fqQkTJqj8/HzD65s8ebJqaWnpc7y5uXnA3/WVuGrHcNavX8/UqVPRdZ3HH3+cr3zlK0yYMIFDhw6xevVqw+u75pprePLJJ1m7di2JiYkA1NfXs3379gsegWOUBx54gGuvvZbVq1f7yq+srKSgoID169cbXt8f//hHTpw40ef43XffzeTJkw3vcYSEhFBUVMTcuXMvOP7ee+8REhJiaF0Ac+bMYePGjaxdu/aC32d+fj6zZs0yvL7HH3+cw4cP9+nNtLS0MG/ePON6HJ9SSvU7tmexWPp9asqVumovqQBqamoASElJobW1lddee420tLQ+f8RGaGlpIS8vj8LCQurr6wFISkrijjvuYNOmTcTFxRleZ0lJCS+99BLV1dUopXA4HCxfvpysrCzD65oxYwa//vWv+/zuioqKWL9+PX/5y18Mre/IkSM89NBDdHR0+C6pKisriYqK4pe//CVz5swxtD6Xy8XWrVspLCzs8/tcv349wcHBhtY3efJk3nvvPaKjoy843tbWRm5uLqWlpYbWl5+fz7/+67+yZMkSX6BWVFSwb98+vve973H//fcbUs9VHTjD7d577+X5558f7mYYwuwAOK+uru6CADg/PhbozAqAz2ppaWHv3r0X/D5vvfVWQ2fjJHBM0t9jbV5//XVuvvlmAF566SVD6zN7Vuw8MwNAmTztb+agOJgTAP3x57T/VTuGY7aqqiqysrJ48MEH0TQNpRTvvfce//iP/+iX+latWsXNN9/MgQMHLpgV2759O1/5ylcMnxWDcwFw5swZ3wdS13USExP9EgBmT/sPNCv29NNPs3v3bsPHqABiY2NZvXq1Kff9nJ/2f/3114mOjvbftL9hw8/ionRdVz/5yU/ULbfcot5//32llFLp6el+q8/sWbG9e/eqiRMnqqVLl6r169er9evXq1tvvVVNnDhR7d271/D6pk6dqsrKyvocP336tF9mGc2eFTtz5oy6++67VXx8vMrIyFATJ05U8fHx6u677+73fQ/V/PnzVUFBgfJ4PL5jHo9H7dy5U82bN8+weiRwTFZZWalWrlyp/uEf/kGlpqb6rZ7FixerJ554QtXV1fmO1dXVqby8PLVo0SLD6zM7AMye9p8+fbp69913+xx/99131bRp0wyvz6wAOM+saX+5pDKZw+HgxRdf5JVXXiEqKspv9fz3f/83eXl53HjjjdTX16NpGomJiSxfvpzf/va3htfn8Xh8lxqfNW7cONxut+H1mT3tv3379gEHxbdv3254fU1NTdx9990XHLNaraxevZrvfe97htdn1rS/DBqPYsePH6eqqor58+cTERHhO75nzx7DBzo3b97Mb3/7234DYNWqVXz3u981tD4wd9r/PLMGxVevXk1cXFy/AdDU1GT4PxpmTftL4IxSTz/9NM888wyZmZl88MEHPPXUU6xYsQKA2bNnc+TIEcPrHI4AMJMycVbM7Pt+zCKBM0pNnz6dd955h4iICMrLy1m5ciX33nsvjzzyCLNmzeL9998f7iYOiSyGNZ4Z0/4SOKNUVlYWJSUlvu87OztZuXIlWVlZvP7663zwwQeG1md2AAy0GHb79u3s379fFsNeJtMWwxo2/CxGlJtuusk3/X6e2+1W9957r7JYLIbXt2TJEpWXl6dqa2t9x2pra9XmzZvVLbfcYnh9shjWWGZN+0vgjFKVlZUXfPg/66233jK8PrMDwOxp/x/96EcqJydH5eXlqR07dqgdO3aovLw8lZOTo370ox8ZXp/Z9/2YNe0v0+KjVH9T1Od94QtfMLw+s1fDmz3t/93vfpcvf/nLFBYW8s477/gGcXfs2OGXQXGzV8ObNe0vYzjCEJ9dDd/Q0ADgC4BNmzb5ff3PwYMHKSoqYvr06aNiAHfULoY1rK8kxACeffZZw8u89tprfV//13/9l8rJyVE/+MEP1HXXXac2b95seH27d+/2fd3a2qrWr1+vpk+frtasWXPBZZ3RamtrVXFxsXrvvfcGvET2h46ODnX48OF+N+UaCgkc4Xf+WMKRk5Pj+zo3N1c1NDQopZTq7Oz0y1KDWbNm+b5ev369euyxx1R5ebn6yU9+olasWGF4fZ/nrwA476GHHvJ9ffDgQZWamqoWLlyoHA6HeuWVVwyrR8ZwhCFmzJjR73GllG/DMSN5vV5aWlrwer0opYiPjwcgPDwcm82/f9bFxcW+2wq+9a1vGb7hPsDGjRv55S9/CcBbb73FV7/6VSZOnMjJkyf5z//8T2677TZD6/vsnszf+9732LVrF7Nnz+b06dOsWrXKsPokcIQh6uvr2bt3b5+xGqUU1113neH1tbW1MWfOHN/WmHV1dSQlJdHZ2WnolpjnNTQ08JOf/MS3bcP5euFc+BnNrADoT3t7O7NnzwZgwoQJ6LpuWNkSOMIQt99+O52dneTk5PT52cKFCw2vr7y8vN/jFouF//3f/zW8vr/927+lo6MDgLVr19LU1ER8fDx1dXX9vmcj+TMAzjt+/DgzZsxAKUV5eTktLS3Exsbi9XoNXXwrs1RCDNLx48eprq5m3rx5fl8MGxYWRkZGhi8AKioqfAEwY8YMjh49amh9Z86cueD7lJQU7HY7TU1NvPnmm9x5552G1COBI8Qg/PznP+cXv/iFaYthzQoAs0ngCDEIo30xbHt7O5s3b6aqqoply5bx1a9+1fezzw5gD5XFkFKEGOV0XfddRo0fP54DBw6we/duHn30Ub8MUre3t/Pd736Xe++9lxdeeOGCn23cuNHw+tatW4dSirvuuouCggLuuusunE4ngKFPFZXAEWIQkpKSLlhhHxERwcsvv0xTU5Phz9wC8wLgvFOnTpGXl8eXv/xlXnrpJWbPns3NN9/M2bNnja3IsDt6hBjFzF4MO3PmzAu+/7d/+zd13XXXqaampgtuQjTK1KlTla7rFxzbvn27ysrKUmlpaYbVI9PiQgyC2YthnU4nXq8Xi+XcRchjjz2Gw+HghhtuoLOz0/D67rjjDl5//XVuueUW37HzC3G/8Y1vGFaPXFIJMQKdD4DPWrt2Lf/+7/9OUFCQ4fXdddddXHvttQD09PTwz//8z9xxxx288cYbFBcXG1aPzFIJEWC2bdvGunXrDC0zOzubDz/8EJvNxoYNGwgLC2PlypXs37+fDz/8kN///veG1COBI0SASUtLo6KiwtAyMzMzfdukfv6+opycHMO2pJUxHCFGILMXw06bNs3Xc5o5cybFxcXk5uZy4sQJ7Ha7YfVID0eIESgxMfGii2FramoMra+trY1HHnmEgwcPMnbsWI4cOUJqaiqpqak8/fTTzJw505B6pIcjxAhk9mLY6Ohotm/fTkdHB6dPn/Y9SfX8drFGkR6OEMI0Mi0uhDCNBI4QwjQSOGJEKi8v77NocTBaW1svWNlcU1PDypUrjWyaGAIJHDEiXSxwPB7PgK/7fOCkpKTwu9/9zvD2iSsjgSMu8NxzzzFjxgxmzpzJvffey5kzZ1i0aBEzZsxg0aJFvhvO7r//fh5++GGuu+46JkyYcMGH+sknn2T69OnMnDmTTZs2AedWIy9dupQ5c+bwxS9+kePHj1+0nE2bNnHw4EFycnL46U9/yvbt2/nKV77CHXfcwZIlS+js7GTRokXMnj2b6dOnU1hY6HvdqVOnyMnJ4dvf/jbl5eVMmzYNgN7eXtatW8f06dOZNWsWb7zxBnDuIXB33nknS5cuZdKkSXznO98x55d9NTJsGagIeEePHlWTJ09WjY2NSimlzp49q26//Xa1fft2pZRSW7du9T0SZe3atWrlypVK13V17Ngx3/O1//jHP6oFCxaorq4uXxlKKXXzzTerEydOKKWUOnTokLrpppsuWs4bb7yhvvSlL/natm3bNjVu3DhfeW63W7W1tSmllGpsbFQTJ05UXq9XlZWVqezsbN/rPvv9j3/8Y3X//fcrpZT6+OOPVWpqqurp6VHbtm1T6enpqrW1VfX09Ki0tDRVUVFh5K9WfEruwxE+r7/+OitXrmTs2LEAxMXF8c477/jW0dx7770X/Ov/5S9/GYvFQlZWlu/u19dee41169YRFhbmK6Ozs5O3336br3zlK77Xnt/bZaBy+rN48WLi4uKAczfA/dM//RNvvvkmFouF6urqS96B+9Zbb/lWPk+dOpVrrrmGEydOALBo0SKio6MByMrK4syZM355RPHVTgJH+KjPPPpkIJ/9eXBw8AWvHagMr9dLTEzMgOtx+iunP+Hh4b6vd+zYQWNjI4cPH8ZutzN+/Hh6e3sv2vaLlf3ZNlit1ouOE4krJ2M4wmfRokX89re/9e3y1tzczHXXXUdBQQFw7kN+/fXXX7SMJUuW8Oyzz9Ld3e0rIyoqivT0dF588UXg3Af/ww8/vGg5kZGRvsey9KetrY2EhATsdjtvvPGGb9Pxi73uhhtuYMeOHQCcOHGCiooKpkyZctF2CGNJ4Aif7OxsHnvsMW688UZmzpzJo48+ytNPP822bduYMWMGzz//PE899dRFy1i6dCnLly8nNzeXnJwcfvzjHwPnwmrr1q3MnDmT7Oxs3yDvQGbMmIHNZmPmzJn89Kc/7fPze+65x7fAcMeOHUydOhWAMWPG8IUvfIFp06bx7W9/+4LXbNy4EV3XmT59OnfffTfbt2+/oGcj/E+WNgghTCM9HCGEaSRwhBCmkcARQphGAkcIYRoJHCGEaSRwhBCmkcARQpjm/wcXoYmgc5DMOgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(facecolor = 'w',figsize=(4,6))\n", "plt.rcParams['pdf.fonttype'] = 42\n", "sns.boxplot(data=df[(df['time']==180)&(df['concentration']!=500000)], \n", " x='concentration', y='Mean', color ='lightgrey',fliersize=0)\n", "sns.swarmplot(x=\"concentration\", y=\"Mean\", data=df[(df['time']==180)&(df['concentration']!=500000)]\n", " , color=\"cornflowerblue\", size=1.5)\n", "plt.xticks(rotation =90)\n", "\n", "plt.savefig('dose_titration.pdf')" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAOQAAAF3CAYAAAC10ufUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAypElEQVR4nO3daXRU1533+++pKlWV5gE0lCgBEsUoJoMswEkcDJGHDKLtYBq3OybYueplp1eSlQ5ur+u12p3n2i25b3c8JO50EztEdhLriTt5wAOGxOTiITYti8EYCYMMktE8z6q59n2Bo9hGAiGdU1WS/p9XcOrU3ruk+mmfYZ+9NaWUQggRFUyRboAQ4i8kkEJEEQmkEFFEAilEFJFAChFFJJBCRBHDAnn33XeTkZHB8uXLR7bt2rWLJUuWsHLlSm699VZ6e3tHXistLcXlcrF48WIOHjw4sv3o0aOsWLECl8vFd77zHeQujZjWlEFef/11dfToUZWfnz+y7eDBg8rv9yullLr//vvV/fffr5RSqrq6Wq1cuVJ5PB51/vx5lZeXpwKBgFJKqWuvvVa9/fbbKhQKqZtvvlnt37/fqCYLEXGG9ZDXX389aWlpn9p24403YrFYAFi/fj2NjY0A7Nu3j+3bt2Oz2cjNzcXlclFZWUlLSwv9/f1s2LABTdO466672Lt3r1FNFiLiInYO+fOf/5xbbrkFgKamJnJyckZeczqdNDU10dTUhNPpvGS7ENOVJRKVPvLII1gsFu68806AUc8LNU0bc/tYdu/eze7duwH44IMPWLJkiU4tFkI/9fX1dHZ2jvpa2ANZXl7Oyy+/zKFDh0bC5XQ6aWhoGNmnsbGR7OxsnE7nyGHtJ7ePpaSkhJKSEgAKCgqoqqoy6FMIMXEFBQVjvhbWQ9YDBw7w6KOP8uKLLxIXFzeyvbi4mIqKCrxeL3V1ddTW1lJYWIjD4SAxMZEjR46glOLZZ59ly5Yt4WyyEGFlWA95xx13cPjwYTo7O3E6nfzwhz+ktLQUr9dLUVERcPHCzn/+53+Sn5/Ptm3bWLZsGRaLhaeeegqz2QzAT3/6U775zW/idru55ZZbRs47hZiONDXaido0IIesIlpd7rspI3WEiCISSCGiiARSiCgigRQiikgghYgiEkghoogEUogoIoEUIopIIIWIIhJIIaKIBFKIKCKBFCKKSCCFiCISSCGiiARSiCgigRQiikgghYgiEkghoogEUogoIoEUIopIIIWIIhJIIaKIBFKIKCKBFCKKSCCFiCISSCGiiARSiCgigRQiikgghYgiEkghoogEUogoIoEUIopIIIWIIhJIIaKIBFKICXj1mJtHf9eve7kW3UsUYgbInxvD7CT9+zMJpBATMHe2hbmz9Y+PHLIKEUUkkEJEEQmkEFFEAilEFJFAChFFJJBCRBEJpBBRRAIpRBSRQAoRRQwL5N13301GRgbLly8f2dbd3U1RURELFy6kqKiInp6ekddKS0txuVwsXryYgwcPjmw/evQoK1aswOVy8Z3vfAellFFNFiLiDAvkN7/5TQ4cOPCpbWVlZWzevJna2lo2b95MWVkZADU1NVRUVFBdXc2BAwe47777CAaDANx7773s3r2b2tpaamtrLylTiOnEsEBef/31pKWlfWrbvn372LFjBwA7duxg7969I9u3b9+OzWYjNzcXl8tFZWUlLS0t9Pf3s2HDBjRN46677hp5jxDTUVjPIdva2nA4HAA4HA7a29sBaGpqIicnZ2Q/p9NJU1MTTU1NOJ3OS7aPZffu3RQUFFBQUEBHR4dBn0II40TFRZ3Rzgs1TRtz+1hKSkqoqqqiqqqK9PR0XdsoRDiENZCZmZm0tLQA0NLSQkZGBnCx52toaBjZr7GxkezsbJxOJ42NjZdsF2K6Cmsgi4uLKS8vB6C8vJwtW7aMbK+oqMDr9VJXV0dtbS2FhYU4HA4SExM5cuQISimeffbZkfcIMS0pg2zfvl1lZWUpi8Wi5syZo55++mnV2dmpNm3apFwul9q0aZPq6uoa2f/hhx9WeXl5atGiRWr//v0j2999912Vn5+v8vLy1Le//W0VCoXGVf/atWt1/0xC6OFy301Nqel5Y6+goICqqqpIN0OIS1zuuxkVF3WEmGqaugIcO+/TvVwJpBATcLzOz4Hjbt3LlUmuhJiArxbE8tWCWN3LlR5SiCgigRQiikgghYgiEkghoogEUogoIoEUIopIIIWYgL7hEBc6ArqXK4EUYgJee8/D068N6V6uDAwQYgKKC2O56Rq77uVKIIWYgBizRox57IflJ0oOWYWIIhJIIaKIHLIKcQWvvvoqra2tn9p2vm82QyYHKxLfv2T/rKwsbrnllgnVJYEUYgI0X7ch55ASSCGuYLTebs+ePUAHO3fu1LUuOYcUIopIIIWIIhJIIaKIBFKIKCKBFCKKSCCFiCISSCGiiARSiCgigRQiikgghYgiEkghoogEUogoIoEUIorI0x5CTEBrYD59wdm6lyuBFGICYrUBlFn/A0w5ZBViApLNXTgs53UvVwIpRBSRQAoRRSSQQkQRCaQQUUSusgoxAd3BLAZDKbqXKz2kEBMQVGb8yqp7uRJIISYg3dLEAutJ3cuVQAoRReQcUogJaPQvpCeUqXu5EkghJiDV3IbdNAgs1LVcCaQQExBv6ieeft3Ljcg55GOPPUZ+fj7Lly/njjvuwOPx0N3dTVFREQsXLqSoqIienp6R/UtLS3G5XCxevJiDBw9GoslChEXYA9nU1MSTTz5JVVUVp06dIhgMUlFRQVlZGZs3b6a2tpbNmzdTVlYGQE1NDRUVFVRXV3PgwAHuu+8+gsFguJstxCVCSqOmwa9rmRHpIQOBAG63m0AgwPDwMNnZ2ezbt48dO3YAsGPHDvbu3QvAvn372L59OzabjdzcXFwuF5WVlZFothAEQ4rX3vPQH0yjMbCIx18aoGcwpFv5YQ/knDlz+MEPfsDcuXNxOBwkJydz44030tbWhsPhAMDhcNDe3g5c7FFzcnJG3u90Omlqahq17N27d1NQUEBBQQEdHR3Gfxgx4wy4FS9VuekIOukLZfDoXSmkJugXo7AHsqenh3379lFXV0dzczNDQ0P88pe/HHN/pdQl2zRt9IUyS0pKqKqqoqqqivT0dN3aLMSfpcSbeOKeVBZYT7LC9qauYYQIBPK1114jNzeX9PR0YmJiuO2223j77bfJzMykpaUFgJaWFjIyMoCLPWJDQ8PI+xsbG8nOzg53s4UIi7AHcu7cuRw5coTh4WGUUhw6dIilS5dSXFxMeXk5AOXl5WzZsgWA4uJiKioq8Hq91NXVUVtbS2FhYbibLcSndAez+Mi/VPdyw34fct26dWzdupU1a9ZgsVi45pprKCkpYXBwkG3btvHMM88wd+5cXnjhBQDy8/PZtm0by5Ytw2Kx8NRTT2E2m8PdbBHlfAFF90CIrNTwfDdCykxQ6R8fTY12kjYNFBQUUFVVFelmiDB57T0PL77r5slvpYalvj179gCwc+fOq37v5b6bMlJHGKK5O8jvT7j52y/GYzGPfhFOT9fn21g+N8bweowmT3sIQ3j9it4hRbgOv6wWLWyHq0aSHlIYIjfTwve+lhjpZkw50kMKEUUkkEJEEQmkEFFEziGFmID+YBrDKkn3ciWQQkzAkEqmN6j/eGk5ZBViAhyWOpba9H8MUAIpRBSRQ1ZhmCFPiHj71Pmb/+qrr9La2nrF/TyhWHraWwlqcfzTf1UyL+Y0Ju3yDylnZWVxyy23XLFsCaSYtNG+yPW9ibRYv0hh7AHMWuBTr433yxlura2tNDc3k5ycPOY+QWI4pX2NJPspgqY4goFZpPm9WPCN+Z6+vr5xt0ECKQxh9dSRGfBijgtceecokpyczMaNGy+7T2zdIN3DC0mLC3Bd7gBw3WX3P3z48Ljrl0CKSRutt7v4NIRvQk9DRLuU2CC9HgtjTFwxKRJIIa7Ssiw3Jg2a+/VfbEcCKcRVGvCaWJLpZkmmW/eyJZDCMCFl/HOQ4eb2a+z5n0xWZg/hC5q4cXGvruVPnWvSYkrpDmZR6fkyHv/0mpAiNkZx68ou5qZ4SYnV/4KV9JDCEEmmLhbEnMAec2Okm6K7nJSLtzjy8OpetvSQwhAWzU+6ZfQJraeDdy8k8JsTs3UvVwIpxFUIhuCt80kk2/0sSh/WvXwJpBBXIRjS+KjHxvmuWN5vide9fDmHFOIqWC2KO9d24AtoeAL692fSQwoxAa0DMdR22HUvVwIpxDgNek0cb4pHKej3WGgf1H8eWDlkFeJj3d3duN3uMQeD9+OgQVvHqQ+HSFJnCBDLT9rXME+9RSJjP7bV29uL1zu+WyTSQwoxTkm0sFjtx0Mqg2TRreWSpU4QS49udUgPKcTH0tLSGBoauuLjV+Dj/RYHjb02blnqBJyX3fvw4cPEx4/viqwEUogJSLEH8CfqP1ZXAinEBOSk+shJHXuWgImSc0ghJuBYYzxtA3KVVYiocK7LjoYiM9Gva7nSQwoxAY5EH6da9R86J4EU4iqdaokD4CvL9Lvd8WcSSCEmYMBrpqYtVvdy5RxSTDmjzQNb15dCl2UlBfFvXLK/3vPALncMYzWHaDFgkivpIYUh+oKzOO65AW+YpvAwe1tJ9J4yvJ5hn4mf/08G7YMx9Hv1X0JdekhhCJs2zCxzMzGWebqXPfY8sN2GzwNrs4TIzxomPcFPkj2oe/njDuTbb79NfX09gcBfJva56667dG+QmB7sJjdzTWcwaZef1XuqMZtg3bxBw8ofVyC/8Y1vcO7cOVavXo3ZfLGb1jRNAimEzsYVyKqqKmpqatCMmDtdCB10BObQHHBFuhmTNq6LOsuXLx/XMl1CREq8qZ90S0OkmzFp4+ohOzs7WbZsGYWFhdhstpHtL774omENE+JqxJkGiDMNAF80tB5fQOP3Z1JwpbvxBUyszNZ35rlxBfKf//mfda1UiKlMAT1DFjqHYyITyC9+0di/OkJMFVaL4mv5PTT3xZBgv/yqyRMxrnPII0eOcO2115KQkIDVasVsNpOUlDThSnt7e9m6dStLlixh6dKlvPPOO3R3d1NUVMTChQspKiqip+cv4wRLS0txuVwsXryYgwcPTrheET5DoSTO+NYSCE6vtT3+rH3Qyvku/WedG1cP+fd///dUVFRw++23U1VVxbPPPkttbe2EK/3ud7/LzTffzH//93/j8/kYHh7mX/7lX9i8eTMPPPAAZWVllJWV8eijj1JTU0NFRQXV1dU0NzfzpS99ibNnz47cfhHRSaGhVPgGgoWUCZ+afED6+vrGteLx4OAgZuBwZ8K4ytR9Cg+Xy0UwGMRsNrNz506uu25iN3z7+/t54403+MUvfgGA1WrFarWyb9++kR/Ejh072LhxI48++ij79u1j+/bt2Gw2cnNzcblcVFZWsmHDhgnVL8IjwdTHQusx3L5lJMYaf7usNTifBv+SSZWRlZU17n27hiz4TclkxnuuuG98fPy4yx5XIOPi4vD5fKxevZr7778fh8PB0NDQuCr4rPPnz5Oens7OnTt57733WLt2LU888QRtbW04HA4AHA4H7e3tADQ1NbF+/fqR9zudTpqapu8iLtPBhY4APmXlgn8p/3FgkH+8deKnN+OVaf6IZFMHsHXCZVzNAPT/tftPDAWz2LlzwYTrG824jimee+45QqEQP/nJT4iPj6ehoYHf/va3E6owEAhw7Ngx7r33Xo4fP058fDxlZWVj7q/UpecgYw1Q2L17NwUFBRQUFNDR0TGh9onJ2/2HQep9y+kLprP983FhqbMvlM5H/vyw1AWQE3OWlfZLnyyZrHEFct68eSilaGlp4aGHHuJHP/oRLtfERkU4nU6cTifr1q0DYOvWrRw7dozMzExaWloAaGlpISMjY2T/hoa/3PBtbGwkOzt71LJLSkqoqqqiqqqK9PT0CbVPTN4/3ppEjuU06ZYGcmaH51w/RvMQb+oLS11GGlcgX3rpJVavXs3NN98MwIkTJyguLp5QhVlZWeTk5HDmzBkADh06xLJlyyguLqa8vByA8vJytmzZAkBxcTEVFRV4vV7q6uqora2lsLBwQnWL8EiMNTFMMs0BFz59p5wZu05TL/NiToenMqAnmMEF3yKefHmA0436fchxDwyorKwcmUB29erV1NfXT7jSH//4x9x55534fD7y8vLYs2cPoVCIbdu28cwzzzB37lxeeOEFAPLz89m2bRvLli3DYrHw1FNPyRXWKSDZ1InLegy79eaw1NcWmEtDYDHGPnz1F35lw62SSNUg1qrfRatxBdJisZCcnKxbpatXr6aqquqS7YcOHRp1/wcffJAHH3xQt/qF8SxagNnmlrDVl2zqRLMorjSL+GQppTh6zk+C1s2Alkp7X5D5Gfo9VjyukpYvX86vf/1rgsEgtbW1PPnkkxO+7SGEEeymYewm/Vc0/qyeIcXTrw2Spi1iKJTMj7br11HBOM8hf/zjH1NdXY3NZuOOO+4gKSmJxx9/XNeGCDEVpCWYePzuVFwxx1lhexOzSd97rOO+D/nII4/wyCOP6Fq5EFOR3aqhaWAmzFN4XOlKqjx+JWaq7mAmQ6EU3cu9bCDfeecdcnJyuOOOO1i3bt2oN+mFmIkCKgavDmNnP+uygWxtbeUPf/gDzz//PL/+9a/5yle+wh133EF+fvhGRIipZ9gb4n3v58mLeS/STTFMhqWRDBqBNbqWe9mLOmazmZtvvpny8nKOHDmCy+Vi48aN/PjHP9a1EWJ6sZg1bAzTFphHaJoeVSl18QkTvV3xoo7X6+WVV17h+eefp76+nu985zvcdtttujdETB9Wi0am5SPq/CvxB8Cm/6ptn3KizsdwKPHjKTzCoyGwhO7g+J8OGa/LBnLHjh2cOnWKW265hYceeojly5fr3gAxPSWbu1ht/v+wxRg/dmb/UTcDgYUElcHJ/4QM8wWSTZ1Arq7lXrbPfe655zh79ixPPPEE1113HUlJSSQlJZGYmDipGQOE0NP/vTUZh6UOuzbIGzUefAHjD5PtpmGSzZ26l3vZHjIU0n/OECH0tqu8h8RgOgoTz785zPwMC3NnT81VMmSxHWGI3mA6Rz1fCstiO7euiyNG8+BWiTz1f6WGJYxtgbmc9el7hRUkkMIgsaYBssz1xISho7puiY1Z5laclrOYdB7KNhab5iZO69e93KnZr4uIGG1dxrG0dIbotqeyZ8+zmLTxnfpMZh3HzuCcSc+pczUGQqnhv8oqxCe1trbS3Nw8rkfxTJixhAYYHh5E48qHrX19k3vaP8P0EYnWTiYzp87VmGVu+XiGAn3n1JFAiquSnJw88qD65dR22DnwQSqf25CE1XLlQI429eLV9Mgf9qbTY1/78TqRVzbZVZXjTAPEof99TwmkMERumoe/WtE1rjCO5Wp6ZEvIhjXQPq7ZECfbGxtJAikMYTFDTopv0uWMt0du6rNypj2WjQuvvO94JkKOFLnKKnTX3GdlyGui32Pm9XNJ1HfbrvymSbJZQiTbA1feMcpJDyl0t/90Ks4UL7UdsSxOdxMMGX8ronvIwtmOWNbmTGwC72ghgRS6W5oxhN2iWJo/zLy0yR+2jofFrLDHhGdk2dlmP55QnCFz+Mghq9BdrFXR57XwUvUsfIHw3KiPs4bISAjPJLD/+61hLviX0BqYr3vZ0kMK3a1xDuELaGQl+iZ1lfVqZCX6yUoMTyDvvzWJf9vTQ08wQ/eypYcUhviox8ah2pSw9ZAt/TG8cS48TyDZYjSyY+pYaqvUvWwJpDDE/DQvxcu7w9ZDevwm+jxTf0Z7OWQVhogxK+alesNWX2aiH5Nm/ETJRpMeUkwLtR12Dp5JiXQzJk16SDEtLM10k5UUpqW2DCQ9pDBE+0AML1WnEgjTpBPVrXH87uSs8FRmIOkhhSFMJkWMSRGea6ywwjHE/DRPmGozjvSQM0T3YIin/zDIsNf4Lqu1PwabOcTNS3sxh+kbptAIqXDF3zgSyBkiEFT0u0MEw3AIeag2hbfqkniuKh1/MDwheb8ljt+cmB2Wuowkh6wzREayme8Xh+fG+dZVnbh9Jk63x2ExGX8fsvJCAmfaY9l+TYfhdRlNekihO5tFkRIXZMP8AbQwdJAL093kpHj53UnpIYWIuNTYIGucQ2Eby2ok6SHFlFffbSPGFCJv9tS/yio9pNDdsM/EB+2xXDNnyPBD1pC6+EB03iwPF3pslGxo072O0Sbb+rArDo85a9RJtSYzgZYEUuiua9hCVUMCdkuIZVluQ+syaXDPujYCIY3OwfAtthNjNhMw4PhSAinGrbu7G7fbPa5JojJx8tqZ62j54A+YufJcN729vXi9ExuMbrMobChUvDHnkJfv7b6oa10SSGEIO73M4kNMBCdcxtX8ARgijQ9NN+MKHSSersvuO5nwG00CKcYtLS2NoaGhcU3L2Nxn5Y3zSXxhZRKWcTymePjwYeLj4yfUrmHS+FC7CWfoCHF0T6iMaCGBFIbITvax/ZrJrZ843j8AwRAk1g/S3LeGPstqFsz2sMIx9rORkwm/0SSQYsp663wSg14T6+cPsCTDjc2iMGmKBGuA022xLM009oKSESSQQldNfVaOfJRIRoKPWXEBQ6+yzkv1cKHXRvm7GdywoI+zHbEEghqZiT7quu0smO3Bag7PFCJ6kUAKXdksIeyWEMGgRlOflcoLiXyzsN2QunJSfWQn+3D7TZzrsrM8a5j0BD/J9gA3Lg7fkyZ6moJNFtFsdnyAAa+ZYb+JziEL187Vf4WoTzKbYGmGm+TYAO81xzMn2UeCLTQlwwgRDGQwGOSaa67hq1/9KnDxEndRURELFy6kqKiInp6ekX1LS0txuVwsXryYgwcPRqrJYpxuW9FF4dwBHEn+sJzHzUnxcY1ziLXOQcPrMlrEAvnEE0+wdOnSkf+XlZWxefNmamtr2bx5M2VlZQDU1NRQUVFBdXU1Bw4c4L777iMYnPi9LWE8q0UxOyHIRlc/YVphnNTYIEszZda5CWlsbOSVV17hW9/61si2ffv2sWPHDgB27NjB3r17R7Zv374dm81Gbm4uLpeLykr9J6gV+mvstTLgCc9X7P3mOH5emRmWuowUkUB+73vf41//9V8xmf5SfVtbGw6HAwCHw0F7+8ULAU1NTeTk5Izs53Q6aWpqGrXc3bt3U1BQQEFBAR0dU/9h1alq0Gvi92dSOFSbQnVbnOH1tQ/E4A1q3OCK3oVYxyvsV1lffvllMjIyWLt27biGRCl16WVrbYxHCEpKSigpKQGgoKBgUu0UExcIafQOm3HNHubanMmd1/X19V3xe9KJi87QfJL8p2m0jf7H+rNlysCAj/3pT3/ixRdfZP/+/Xg8Hvr7+/nbv/1bMjMzaWlpweFw0NLSQkbGxYVMnE4nDQ0NI+9vbGwkOzs73M0WVyElNsi6eYP84WwKhXOHMDOxe4FZWVnj2i+eFnxdQ/TY1rAwvvfK+8fHj7vscAt7IEtLSyktLQUuDmH6t3/7N375y1+ya9cuysvLeeCBBygvL2fLli0AFBcX8zd/8zd8//vfp7m5mdraWgoLC8PdbHGVggp8QW2CUbzoap4p/Mkz/01csJGdO3dOosbIi5qBAQ888ADbtm3jmWeeYe7cubzwwgsA5Ofns23bNpYtW4bFYuGpp57CbJ76i6pMd9lJPja5+sI2UsYdSqQr6AhLXUaKaCA3btw4MnB41qxZHDp0aNT9HnzwQR588MEwtkxMlj1GsSSMY0lTzO3YTMOAK2x1GmGKjmcQ0a7PbeZPdYmEwjSUtD04lxrvhvBUZiAJpDDEoNfMR912Ah+P4TB6jQ+HuY5r7H80tpIwiJpzSDE1jOc2RIO2jqCvg3j/W/zXn75MivqIfm0O+WrvZcud/K2IqfVkx2gkkNPQaLOkdXV349VmkZ166Zd2vLOkjfdWQXqwi1avheb4YvKsJ0k0deNV3cSbxw7cZG9FtARzafQvnvD7o4UEcoboD6bRklBMWugP2E0TG/M53tsQ/3VwEO9HVcSrD3ngWzeMOZBDTxnmBhK0HuCvDK/LSBLIaWi04OzZs4eM0B+5956/Nrz+7DQzHR+Z6Azm4AuAzeDZGbsGggRUDDaTm7dOe/n8UpuxFRpILurMIPEmY59N/LOi1XZSzO0UxP4eW4zxveMLb7s5519Fiz+P3x0ZJhSuS7sGkEAK3b1R7aHa+3kGQ8lhqe+ujXGkmNrpCObwo52pmML1zJcBJJBCdzcst5NjOc1J7/V4fMb3VnE2E3NizrHaNvrAkqlEziGF7mIsGlmWelLMHdit28JSZ2tgPg1ylVVMFQEVQ0sgF39AEWMx/pBO0yBWGzK8nj9LM7dg1dzA1H4SSA5ZZwiviqUtOI8hb3gueASVmd5g+BZQtWpe0sxt1LcH8AXkoo6YAsJ5qaMnlMFp3wbDzyGDIcWu8h4u+BfxoXcVj/6un3drfYbWaSQJ5Aww4A5h1dykmy8QbzM+lv/0fB+DwVTitR7Dp2M0mzS+sjaWWG2IkGbhf/1NMusXW42t1EASyGnu6Dkf95f3goK5MWcY9Cj+91tDnKgzrhf58ho7yaYOEkx9mMLwDdu43E66pYlF1qOkJ5kxy20PEa3y58bwlQI7R7030e6fwz8+14vbpwgZ+PTF+sU2Ui0d5FnfD1s4uoIOTnmvC0tdRpKrrNOcPUZj88pYjh07zrnAau75UjzrFk3doWVjsWnDJJkuvy7kVCCBnAG6BoLMMjdhIkiB6wuRbo4hEkx9JJj6gKndS0ogZ4Bf/HGIAf9SukPZdA+ESE82fk4ipcCjonOqxWgm55AzwPe+lkiu5SSZ5noudPg5es742wJdwWyOezeFZejcdCI95BQ22oPInxRUJloCeTgsdbR2e2mPd/J/DtUAipNvnLxi+eN9cHk0aeYWVprewG69bULvn6kkkFNYa2srzc3NJCeP/lSFlwSatTz6/AmYLO0k+07iIZ0Y3Az5Lz+sra9vctPymzRFgjb1p/YPNwnkFJecnDwyleZobqKHk83JhFQyJ5oSWDevnyR7HHOSx34PMK5lHoT+5Bxymuv3mFmZPcy8VC/2mBDzUn3MSZ66Q8umOwnkNNbvMbOnMpPjTfG4/SYSbEFqO22cbbdHumliDHLIOo0l2YN8fWUnb55PoqXfyrnOWOyWEIk24xe8DSoLfaHwPe0xXUgPOc05U3zcvqqTuJgga5wDrJ83QLw1RGu/sTNP9YbSOeMrkNseV0l6yCmsu7sbt9t9xQswg6TTzxxCmDne4CKeDpJVI7M5O+Z7ent78Xq9E27bLHML19oPYLd+Y8JlzEQSyBmgQ1sKKAbJZK56mxQarvgePVi0QFjqmU4kkFNYWloaQ0NDl73tAfCFEPiDGjVtbtLjl9IxtJI1zsvfhzx8+HDUrjI8nUkgp7l+j5lnqzJYM2eQkNL4yGeibcBKalyA3LSJH5J+0mgjhppbO/CaM9izZ88l+09mBNB0Jxd1prkEW5ClGcMopej3mDjamEiCLciR+kRD6/XZF9AS/1WCSv7mXw35aU1zJg0GvGaCSqO2M46vr+zAmeJH6Xjxc7TeTilFz5AiLUEu6lwNCeQUd7nl4RRwXruBdPUOfe4+YuzX8/6Jd3kfRT85ZHLqsuVO5hxS0zTSEqbuVBqRIoGcwq60fJtSkOB3068W0m+zYQ4Ncz7mBlJNHQypZHJt9Wja6F3lZJeHExMjgZzCxnth5EJngF/97nXqAiv5+68ksyDTQbzdBOQa20Bx1eSizgzwf4646Q45SNHaWZBp+TiMIhrJb2YGuGtjPDmWDwjDuqlikiSQM0BqgolEcx+Lre9yuNpL/7CBc0CKSZFzyBng8CkPXYEs+kOz+KDay5I5FpLi5G9xNJJAzgBnmwMMhFLpCWXx/3w9iZR4CWO0kt/MDFByYwLzradZYz8kYYxy8tsRIopIIIWIIhLIGSSoTLT3BXmvXia5ilZhv6jT0NDAXXfdRWtrKyaTiZKSEr773e/S3d3NX//1X1NfX8/8+fP5zW9+Q2pqKgClpaU888wzmM1mnnzySW666aZwN3tK8wcUfYFZnPavh0o3DV1BVs2fumsojva4V0PbEB6LY8o/7hX2HtJisfDv//7vnD59miNHjvDUU09RU1NDWVkZmzdvpra2ls2bN1NWVgZATU0NFRUVVFdXc+DAAe677z6CQeMnaZpOjtf5qPFfh8P8IXd8IY5/2pYU6Sbpzm+fS4+9INLNmLSw95AOhwOHwwFAYmIiS5cupampiX379o08tbBjxw42btzIo48+yr59+9i+fTs2m43c3FxcLheVlZVs2LAh3E2fsgpcVl4/XMVZfyFDXjXlh85dvrfbGbZ2GCGi9yHr6+s5fvw469ato62tbSSoDoeD9vZ2AJqamli/fv3Ie5xOJ01NTaOWt3v3bnbv3g1AR0eHwa2PXqMd0rm7hrDFOtn72/3EaP5PvTaVDummu4j9qRwcHOTrX/86jz/+OElJYx9CqVGepNXGGJRZUlJCVVUVVVVVpKen69bW6cBqtWDVPGjItIzRLCI9pN/v5+tf/zp33nknt912cXWkzMxMWlpacDgctLS0kJGRAVzsERsa/jJLWmNjI9nZ2ZFo9oR9tscaDiWgAQO9HVhwk5aW9qn9J9tjXf69U/88azoLew+plOKee+5h6dKlfP/73x/ZXlxcTHl5OQDl5eVs2bJlZHtFRQVer5e6ujpqa2spLCwMd7N1dcG/jDr/cj6M3U5/MDXSzRHRRIXZm2++qQC1YsUKtWrVKrVq1Sr1yiuvqM7OTrVp0yblcrnUpk2bVFdX18h7Hn74YZWXl6cWLVqk9u/fP6561q5da9RHmDSPL6Q8vpD695+9rJ55Zk+kmyPC7HLfzbAfsn7+858f9bwQ4NChQ6Nuf/DBB3nwwQeNbFZY2WIungOnmtsj3BIRbab29e8pKhBUPFTRR08wI9JNEVFGAhkBZhOsnh9Dsz8PTyg20s0RUUQCGQGapnHjajsmLYRbxdPSIyOPxEUSyAh44e1h9h9zE8RCc2Ahe/9nmD6ZVkMggQyr2hY/zd0BAkFFVoqZGHzE0cOsRBM/+EUvTd2yWtRMJ1N46Gi0IWufdNLzBeymITqDTlwxVXQFC0j01XDy/dNkmDUO7j152ZnhZIjb9CeB1FFrayvNzc0kJyeP+nouv0cLKmzk8aHvWma536TbXogl0ITSzAz7x14irq+vz6hmiygigdRZcnLyFddr7HWbONc5wFt117N2zgBLsxIwaZAaN/b7rrRKspge5BwyAg6dTaW5z8oqxwBHmxKpvJDI3vdnAeD2m3jzfBL+oMxqPBNJD6mj7u5u3G73qL1Zs3YNVjVIABuxnGaQdC5oLuar1wm0+YljFr87PMSA5sSNHc+FI8TgHnl/b28vXq8+C6yK6CWBDBMTAUJY6NbysDEL0Ahgp0fLxY8dLynMUmcw4WeJ2h/p5ooIkUDqKC0tjaGhocueQ37Q5uf1c1nMT/VgjxlmWZaNc512TreZuGPdrI/3uvT9hw8fntR6jWJqkHPIMBr0mjjbYWf9vAHibCG8QRN/OJNKr8fC5/P6I908EQUkkGEUUhr+oImu4RgWzh6mvttGvmMIb0Cjxy0HK0IOWXV3uSXGAeJJprb3c5xr1iDQwOEP55MZOkl7ZxuHz489D9BklxgXU4MEUkdXWgI8oCx0BVzkm4/gU3ZOh1aT6q5kQVrLx3uMHThZYnxmkEDq6ErD2lp7gzz58gDDSfmYTRoJjY3MjR3E49jKuoU21iyYupMXC33IOWQYZaWYuedLCQx7FSvmWvCoeD7yL8Ni1kiIlYEAQgIZdgl2jbwsCy9VuUk1taE0E5tW2FmUHRPppokoIIEMs8wUM/PTLWSlWAgqEwqNBVly5iAukkBGQGaKGVeWhfZQLnMtH0S6OSKKyJ/mCFiQZWFBloWPav5Ioqkn0s0RUUQCGQajPbislEZHbwpuf8clS6jJg8gzlxyyRohX2emKvQ5lk6kgxV9IDxkGY/V2IaUwaTvC3BoRzWZkIEc7hAwqM0fdN+DwvUVOiudTrxl1CGm63AQ6YkaSQ9aPmQiS5HkPzSvT+4vImZE95Fi93RNP7yPBlMLOnXIYKSJDesiPeXyKat/n6A7JAG4ROTOyh/wzf0Dx2kkPBQusaBost76FSfk4ds6H26foGQrx1mkvZd9IiXRTxQwxowPZNxzitfc81LUF6BsOUef7PLPMjeytHEYp2LDYxl+tk8VwRPjMuED6A4oYy8Wrm7OTzPy/O1Jw+xT+IJQ//xIJpl5uvmk5L73rZvNK+8hajkKEw7QO5Gdvb/iUnaOeL5FvfZskczcdASd1/hUkal0kmzvpHkomJqT44DdnMaF47rmjhDBT71/OvJhTWDXfp8qXETVCb9M6kJ+d2l8xhINj1HivJU3VMaylkaCaMONhMGjBa07DHPIQCnoZJo1zwy6GtFmYCTLk9+JneKRsmdpfGGFaBxIundo/pOB/PvKRkZDBsM/Euw2zyJ3t5qMeOwnArPhEOgZjWJbqJSclHV/QxJJMP1D4qXJlan9hhGkdyLFmEq/VbsRODznqXUza56hrVCjAghvPYDtBLZfq4TRamprJVO/TcNpGDJ+eNVxmEhdGmNaBHMtc9SdqtSJicGNWPhJppktbhJ84Alo8IWJIVvV0ai78xNKrzWe5+g1yeUcYbVoPDEhLS0P7zHjREGYGcJClTmLGi0dLplHbQKzqJhgM4QnaiFHDDGnZZKsTzOYDHOooAeyfKkfTNNLS0sL5ccQMMK17yNGmTRwKJdHiXUOi1oVZC2JFkag1YlYxmLw+rMEu4uNjsdEHWjL1wdWY8aPMacyPqRkpR6ZlFEaY1oEc65aEUooLHWl4/CFer/aRGm+iZyhE24dDJNsDBJNczE+3sNRpYdgL6xbFYLfOI8a8LsyfQMw00zqQY9E0jXkZFz/64jlWBtwhfvvOMKusr3E+sIZV2TGcuuDnni8lRLilYqaZkYH87IABbyiWs741BLyL0Riio/YQWcCePV2ADAAQ4TOtL+qMl83kJt/2DmazhYzQKVLMXSSbuyLdLDEDzcge8vK93epwNUOIS0gPKUQUmTKBPHDgAIsXL8blclFWVhbp5ghhiCkRyGAwyLe//W1effVVampqeP7556mpqbnyG4WYYqZEICsrK3G5XOTl5WG1Wtm+fTv79u2LdLOE0N2UCGRTUxM5OTkj/3c6nTQ1NUWwRUIYY0pcZVVKXbLts2NUAXbv3s3u3bsB+OCDDygoKLjqujo6OkhPT7/6Rk5QOOubzp9tKtVXX18/5mtTIpBOp5OGhoaR/zc2NpKdnX3JfiUlJZSUlEyqroKCAqqqqiZVRrTWN50/23Spb0ocsl577bXU1tZSV1eHz+ejoqKC4uLiSDdLCN1NiR7SYrHwk5/8hJtuuolgMMjdd99Nfn5+pJslhO6mRCABvvzlL/PlL3/Z8Home8gbzfVN5882XerT1GhXTIQQETElziGFmCkkkB8zemje3XffTUZGBsuXLx/Z1t3dTVFREQsXLqSoqIieHv2WN/d4PBQWFrJq1Sry8/N56KGHDK+zt7eXrVu3smTJEpYuXco777yja32j/Qx37drFkiVLWLlyJbfeeiu9vb0jr5WWluJyuVi8eDEHDx7Upb4TJ06wfv16Vq9eTUFBAZWVlbrVB4ASKhAIqLy8PHXu3Dnl9XrVypUrVXV1ta51vP766+ro0aMqPz9/ZNuuXbtUaWmpUkqp0tJSdf/99+tWXygUUgMDA0oppXw+nyosLFTvvPOOoXXedddd6mc/+5lSSimv16t6enp0rW+0n+HBgweV3+9XSil1//33j5RfXV2tVq5cqTwejzp//rzKy8tTgUBg0vUVFRWp/fv3K6WUeuWVV9QXv/hF3epTSinpIQnP0Lzrr7/+kkmx9u3bx44dF5e+27FjB3v37tWtPk3TSEi4OOOB3+/H7/ejaZphdfb39/PGG29wzz33AGC1WklJSdG1vtF+hjfeeCMWy8Vrk+vXr6exsRG4+LPdvn07NpuN3NxcXC7Xp3qzidanaRr9/f3Axcmy/3w/XI/6QA5ZgcgNzWtra8PhcADgcDhob9d3sdhgMMjq1avJyMigqKiIdevWGVbn+fPnSU9PZ+fOnVxzzTV861vfYmhoyPDP+Ek///nPR551Nep3+vjjj7Nr1y5ycnL4wQ9+QGlpqa71SSAZ/9C8qcZsNnPixAkaGxuprKzk1KlThtUVCAQ4duwY9957L8ePHyc+Pj6sj8k98sgjWCwW7rzzTsC43+lPf/pTHnvsMRoaGnjsscdGjgj0qk8CyfiH5uktMzOTlpYWAFpaWsjIyDCknpSUFDZu3MiBAwcMq9PpdOJ0Olm37uLMfFu3buXYsWNh+Yzl5eW8/PLL/OpXvxoJgVG/0/Lycm677TYAbr/99pHDUr3qk0ASuaF5xcXFlJeXAxd/0Vu2bNGt7I6OjpErjm63m9dee40lS5YYVmdWVhY5OTmcOXMGgEOHDrFs2TJDPyNcvDr+6KOP8uKLLxIXFzeyvbi4mIqKCrxeL3V1ddTW1lJYWHiZksYnOzub119/HYA//vGPLFy4UN/6rvoy0DT1yiuvqIULF6q8vDz18MMP617+9u3bVVZWlrJYLGrOnDnq6aefVp2dnWrTpk3K5XKpTZs2qa6uLt3qe++999Tq1avVihUrVH5+vvrhD3+olFKG1nn8+HG1du1atWLFCrVlyxbV3d2ta32j/QwXLFignE6nWrVqlVq1apX6u7/7u5H9H374YZWXl6cWLVo0cmV0svW9+eabas2aNWrlypWqsLBQVVVV6VafUkrJSB0hoogcsgoRRSSQQkQRCaQQUUQCKUQUkUAKEUUkkDNUb28v//Ef/wFAc3MzW7dujXCLBMgDyjNWfX09X/3qVw0dTieu3pSZwkPo64EHHuDcuXOsXr2ahQsXcvr0aU6dOsUvfvEL9u7dSzAY5NSpU/zDP/wDPp+P5557DpvNxv79+0lLS+PcuXN8+9vfpqOjg7i4OH72s5+xZMmSSH+sqW9CwwnElFdXVzfynN8n/71nzx61YMEC1d/fr9rb21VSUpL66U9/qpRS6nvf+5567LHHlFJKbdq0SZ09e1YppdSRI0fUDTfcEP4PMQ1JDykuccMNN5CYmEhiYiLJycl87WtfA2DFihWcPHmSwcFB3n77bW6//faR93i93kg1d1qRQIpL2Gy2kX+bTKaR/5tMJgKBAKFQiJSUFE6cOBGhFk5fcpV1hkpMTGRgYGBC701KSiI3N5cXXngBuPgs4Hvvvadn82YsCeQMNWvWLD73uc+xfPlydu3addXv/9WvfsUzzzwzMomWrEamD7ntIUQUkR5SiCgigRQiikgghYgiEkghoogEUogoIoEUIopIIIWIIhJIIaLI/w9TT/z8n319gwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(facecolor = 'w',figsize=(3,6))\n", "plt.rcParams['pdf.fonttype'] = 42\n", "sns.boxplot(data=df[(df['concentration']==500000)|(df['concentration']==0)], \n", " x='time', y='Mean', color ='lightgrey',fliersize=0)\n", "sns.swarmplot(x=\"time\", y=\"Mean\", data=df[(df['concentration']==500000)|(df['concentration']==0)], \n", " color=\"cornflowerblue\", size=1.5)\n", "plt.ylim(0,1200)\n", "plt.savefig('time.pdf')" ] }, { "cell_type": "code", "execution_count": 29, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAF3CAYAAABUsGfpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABBeklEQVR4nO3deVzU950/8NdcDPcpIDioKIhcSpRobLq5XJKYJpirRms2JDa1q2mz/TWbNK1tjm6MZttukyY2XRvbkMRIN2kj5lDbZHM1x1I8YgAPRFAGhnuAYe6Z7+f3x8BEFBjEOeH1fDzyCH6Z4+0X/L7m+zllQggBIiKiMcgDXQAREQU/hgUREXnEsCAiIo8YFkRE5BHDgoiIPGJYEBGRRz4Li3Xr1iElJQUFBQXuYw8++CDmz5+PBQsW4JZbbkFvb6/7e1u2bEFWVhZycnKwf/9+9/EDBw6gsLAQWVlZuP/++8GRvkRE/uezsLj77ruxb9++YcdKSkpQU1ODI0eOYN68ediyZQsAoK6uDhUVFaitrcW+ffuwceNGOJ1OAMCGDRuwfft21NfXo76+/rzXJCIi3/NZWFxxxRVITEwcduzaa6+FUqkEAFx22WXQarUAgMrKSqxevRpqtRqZmZnIyspCVVUVdDod+vv7sWzZMshkMtx1113YvXu3r0omIqJRBKzP4g9/+ANWrFgBAGhpaUFGRob7exqNBi0tLWhpaYFGoznvOBER+ZcyEG+6efNmKJVKrF27FgBG7IeQyWSjHh/N9u3bsX37dgDAsWPHMH/+fC9VTEQ0NTQ1NaGrq+u8434Pi/Lycrz11lt477333Bd+jUaD5uZm92O0Wi3S09Oh0WjcTVVnHx/N+vXrsX79egBAcXExqqurffS3ICKanIqLi0c87tdmqH379uGpp57Cnj17EBkZ6T5eWlqKiooKWK1WNDY2or6+HkuWLEFaWhpiYmLw+eefQwiBl156CStXrvRnyUREBB/eWaxZswYffPABurq6oNFo8Pjjj2PLli2wWq0oKSkB4Ork/t3vfof8/HysWrUKeXl5UCqV2LZtGxQKBQDg+eefx9133w2z2YwVK1a4+zmIiMh/ZJN1iXI2QxERXbjRrp0B6eAmIpqs7HY7tFotLBZLoEsZU3h4ODQaDVQq1bgez7AgIvIirVaLmJgYzJ49e8zRm4EkhEB3dze0Wi0yMzPH9RyuDUVE5EUWiwVJSUlBGxSAawpCUlLSBd39MCyIiLwsmINiyIXWyLAgIgqA6OjoMb/f1NQ0bCHW8bj77rvx+uuvX0xZo2JYEBGRRwwLIqIAGhgYwPLly7Fo0SIUFhaisrLS/T2Hw4GysjIsWLAAt99+O0wmEwDX1g1XXnklFi9ejOuuuw46nc7ndTIsiIgCKDw8HG+88QYOHjyI999/Hw888IB7Xbzjx49j/fr1OHLkCGJjY/Hb3/4Wdrsd3//+9/H666/jwIEDWLduHTZt2uTzOjl0logogIQQ+MlPfoKPPvoIcrkcLS0taG9vBwBkZGTg8ssvBwDceeed+M1vfoPrr78eNTU17pUwnE4n0tLSfF4nw4KIKIB27tyJzs5OHDhwACqVCrNnz3YPaT13xNLQatz5+fn47LPP/Fonm6GIiAKor68PKSkpUKlUeP/993H69Gn3986cOeMOhV27duHrX/86cnJy0NnZ6T5ut9tRW1vr8zoZFkREAbR27VpUV1ejuLgYO3fuHLYPT25uLsrLy7FgwQL09PRgw4YNCAsLw+uvv44f/ehHWLhwIYqKivDpp5/6vE42QxERBcDAwAAAYNq0aaM2KdXV1Y14vKioCB999NF5x1988UWv1Xcu3lkQEZFHDAsiIvKIYUFERB4xLIiIyCOGBRERecSwICIijxgWRESTzLp165CSknLBS5yPhWFBRBRAe/fuxY033ohLL70UN954I/bu3XvRr3n33Xdj3759XqjuK5yUR0QUIHv37sXmzZvda0G1tbVh8+bNAIAVK1ZM+HWvuOIKNDU1eaNEN4YFEZGP/OpXv8Lx48dH/f6XX34Ju90+7JjFYsHPf/5zvPHGGyM+JycnBw888IBX6xwPNkMREQXIuUHh6Xgg8c6CiMhHPN0B3HjjjWhrazvv+PTp07F9+3ZflTUhvLMgIgqQ++67D+Hh4cOOhYeH47777gtQRaNjWBARBciKFSuwadMmTJ8+HTKZDNOnT8emTZsuqnMbANasWYNly5bh+PHj0Gg02LFjx0XXymYoIqIAWrFixUWHw7l27drl1dcDeGdBRETjwLAgIiKPGBZEROQRw4KIyMuEEIEuwaMLrZFhQUTkReHh4eju7g7qwBBCoLu7+7xhu2PhaCgiIi/SaDTQarXo7OwMdCljCg8Ph0ajGffjGRZERF6kUqmQmZkZ6DK8js1QRETkEcOCiIg8YlgQEZFHDAsiIvKIYUFERB4xLIiIyCOGBRERecSwICIijxgWRETkEcOCiIg8YlgQEZFHPguLdevWISUlBQUFBe5jPT09KCkpQXZ2NkpKSqDX693f27JlC7KyspCTk4P9+/e7jx84cACFhYXIysrC/fffH9QrORIRTVY+C4u7774b+/btG3Zs69atWL58Oerr67F8+XJs3boVAFBXV4eKigrU1tZi37592LhxI5xOJwBgw4YN2L59O+rr61FfX3/eaxIRke/5LCyuuOIKJCYmDjtWWVmJsrIyAEBZWRl2797tPr569Wqo1WpkZmYiKysLVVVV0Ol06O/vx7JlyyCTyXDXXXe5n0NERP7j1z6L9vZ2pKWlAQDS0tLQ0dEBAGhpaUFGRob7cRqNBi0tLWhpaRm23vrQ8dFs374dxcXFKC4uDvq15ImIQklQdHCP1A8hk8lGPT6a9evXo7q6GtXV1UhOTvZqjUREU5lfwyI1NRU6nQ4AoNPpkJKSAsB1x9Dc3Ox+nFarRXp6unvHqXOPExGRf/k1LEpLS1FeXg4AKC8vx8qVK93HKyoqYLVa0djYiPr6eixZsgRpaWmIiYnB559/DiEEXnrpJfdziIjIf3y2reqaNWvwwQcfoKurCxqNBo8//jgefvhhrFq1Cjt27MDMmTPx2muvAQDy8/OxatUq5OXlQalUYtu2bVAoFACA559/HnfffTfMZjNWrFiBFStW+KpkIiIahUxM0okLxcXFqK6uDnQZREQhZbRrZ1B0cBMRUXBjWBARkUcMCyIi8ohhQUREHjEsiIjII4YFERF5xLAgIiKPGBZEROQRw4KIiDxiWBARkUcMCyIi8ohhQUREHjEsiIjII4YFERF5xLAgIiKPGBZEROQRw4KIiDxiWBARkUcMCyIi8ohhQUREHjEsiIjII4YFERF5xLAgIiKPGBZEROQRw4KIiDxiWBARkUcMCyIi8ohhQUREHjEsiIjII4YFERF5xLAgIiKPGBZEROQRw4KIiDxiWBARkUcMCyIi8ohhQUREHjEsiIjII4YFERF5xLAgIiKPGBZEROQRw4KIiDxiWBARkUcMCyIi8ohhQUREHjEsiIgmCYPFjsPNephsDq+/dkDC4te//jXy8/NRUFCANWvWwGKxoKenByUlJcjOzkZJSQn0er378Vu2bEFWVhZycnKwf//+QJRMRBS0JEmguceEA6f16DTY4JSE19/D72HR0tKC3/zmN6iurkZNTQ2cTicqKiqwdetWLF++HPX19Vi+fDm2bt0KAKirq0NFRQVqa2uxb98+bNy4EU6n099lExEFJbPNiS+0vTjZMYD4iDCoFDKfvE9A7iwcDgfMZjMcDgdMJhPS09NRWVmJsrIyAEBZWRl2794NAKisrMTq1auhVquRmZmJrKwsVFVVBaJsIqKgIYRAW68ZVY3dMNucmBathkLum6AAAhAWM2bMwL//+79j5syZSEtLQ1xcHK699lq0t7cjLS0NAJCWloaOjg4ArjuRjIwM9/M1Gg1aWlpGfO3t27ejuLgYxcXF6Ozs9P1fhogoAKwOJ+pa+3G0zYCYcBViwlU+f0+/h4Ver0dlZSUaGxvR2toKo9GIV155ZdTHC3F+25tMNnJ6rl+/HtXV1aiurkZycrLXaiYiChZdBgv+0dgDvcmGadFqqBT+uYz7PSzeffddZGZmIjk5GSqVCrfeeis+/fRTpKamQqfTAQB0Oh1SUlIAuO4kmpub3c/XarVIT0/3d9lERAFld0o40d6PL7R9iFApERcR5tf393tYzJw5E59//jlMJhOEEHjvvfeQm5uL0tJSlJeXAwDKy8uxcuVKAEBpaSkqKipgtVrR2NiI+vp6LFmyxN9lExEFTJ/ZjgNNPdD1WpAcrUaY0v/dzUp/v+HSpUtx++23Y9GiRVAqlbjkkkuwfv16DAwMYNWqVdixYwdmzpyJ1157DQCQn5+PVatWIS8vD0qlEtu2bYNCofB32UREfueUBM70GNHYaURMuAqJUb7vmxiNTIzUKTAJFBcXo7q6OtBlEBFNyIDVgWO6fgxYHUiIDIN8lL7ac3UbrVg8K2HCnd6jXTv9fmdBRESjE0KgpdeM+vYBRKgUSIpSB7okAAwLIqKgYbE7cbzNgG6jFYmRvp03caEYFkREASaEQKfBimNt/VDI5EiODg90SedhWBARBZDNIeFkhwFtfRbER4b5bd7EhWJYEBEFiN5oQ52uH5IQmBatHnXCcTBgWBAR+ZnDKaGpy4jTehPiwlVQK4N/OgDDgojIj/otdhxt7YfF7kRyVHDfTZyNYUFE5AdDe040dA0gKkyJxCAZEjteDAsiIh8z2Rw4quuHweIIuiGx48WwICLyESEE2vosON5ugFoZPBPsJoJhQUTkAxa7E/UdBnQabEiIUEEZpENix4thQUTkZZ0GC461GSCXyZAcHbp3E2djWBAReYnNIeFU1wBae82D+2GH9t3E2RgWRERe0GtyTbBzOAWmhdCQ2PFiWBARXQSnJHC624imbhNi1ErEqIN/gt1EMCyIiCbIYLHjmK4fRpsTSVHj33MiFDEsiIgukCS59pxo6AyuPSd8iWFBRHQBzDYnjrX1o9dsQ0JEaE6wmwiGBRHROAgh0D44wU6lkGNaVPDtOeFLDAsiIg+sDidOtg+g3WCZdENix4thQUQ0hu4BK47q+gEgKHew8xeGBRHRCOxOCY1dA2jucU2wC1MG991EuEqOMKUckWEKGK0ORIUpIfdifwrDgojoHH1mO47q+mG1O5Ec5DvYAa6gsDgk3PPiP6DVm6FJiMDv7ypGTmqM1wIjuKOSiMiPnJJAY9cADjT1QA4ZEkNkJnaYUo4NrxyAVm8GAGj1ZnznpWp0G21eew/eWRARATBav9pzIilaHVIT7Mw2pzsohmj1ZtgcTq+9B+8siGhKE0JAqzehqrHHta5TCAVF94AVT75zFKd7TNAkRAz7niYhAmFe3NubYUFEU5bF7sQRbR9OtBuQEBmGKHVoNLZIQuDtL3XYsPMgDpzRo7nbiOfXLnIHxlCfRVJUmNfeMzTODBGRFwkh0Gmw4rh7z4nQGRJ7utuI594/iWNtBhRlxGPjVXORFheBcJUcO+9dCodTIEqtQEpMOEdDERFNlM0hoaFzALpeM+IjQ2eCnc0h4U/VzfjzQS2iwhT4Yck8XDUv2d0Bb7FLsNgldButWDwrwatBATAsiGgK0Rtde05IkqtvIhRGOgHAEW0vtr1/Eq19FlwzPwXrLs9EXITKrzUwLIho0nM4JTR1GXFGb0JsuApqL3b8+lK/2Y4/ftqId492IC0uHE+sLMDCjPiA1MKwIKJJrd9ix9HWfljsTiRFhcZIJyEEPjzRiRf+3ogBqwPfXKzBHZdmBDTkGBZENClJkmtIbEOnEZFhCiSGyJ4Tbf0WPP/BSRw804uc1Bjcd3UWMqdFBboshgURTT4m21cT7BIiw0JizwmnJFB5uAU7q85AIZPhu1fMwYqCtKCpnWFBRJOGEAJtg3tOqBWhs4NdfbsBz31wEqc6jViamYjvXjEXyTHBVTvDgogmBYvdifoOAzoNNiREqKAMgSGxZpsTO//vNN480or4iDD8eMV8LJuTFJSjtBgWRBTyOg0WHHNPsAuuT+SjqW7qwW8/bECnwYoVBdNRtmx2UM8gD97KiIg8sDtdE+xae82ICw/+PScA11yP3//9FD6u70JGYiSeum0B8tJiA12WR+MOi08//RRNTU1wOBzuY3fddZdPiiIi8qTPZEetrs+1+F8ILCUuCYG/1bXjj582wmqXcOfSmbh1kSZkZpCPKyz+5V/+BQ0NDSgqKoJC4RrnK5PJGBZE5HdOSeB0txFN3SbEqJWIUQf/BLtmvQnb3j+J2tZ+FKTH4r6rs6BJiAx0WRdkXGFRXV2Nurq6oE9uIprcDBY7jusMMNqcSIoKC/oJdnanhNcPaPE/1c0IVylw/zVZ+Ofc1JC8lo4rLAoKCtDW1oa0tDRf10NEdB5JEmjpNeNkx8DgBDvvLb3tK7Wtfdj2/kk06824IjsZ9/5TJhIig7/u0YwrLLq6upCXl4clS5ZArf5qpMGePXt8VhgREeAaXnq83QC9yYqESHXQTFIbzYDVgfJPm7Cvtg0pMWo8elMeimclBrqsizausHjsscd8XAYR0XBCCLQPTrBTKeSYFhXce04IIfBpQzf++6MG9JntuLloBtYunYlwVfD3qYzHuMLiyiuv9HUdRERuVocTJ9sH0G6wID4i+Pec6DRY8bsPG1DV1IO5yVF45MZ8ZKVEB7osrxrXT+Dzzz/HpZdeiujoaISFhUGhUCA2duLjgnt7e3H77bdj/vz5yM3NxWeffYaenh6UlJQgOzsbJSUl0Ov17sdv2bIFWVlZyMnJwf79+yf8vkQU/LoHrPhHUw/0JhuSo8ODOiicksCeL1px36sH8YW2F9++PBO/+mbRpAsKYJxh8b3vfQ+7du1CdnY2zGYzXnjhBXzve9+b8Jv+27/9G66//nocO3YMX3zxBXJzc7F161YsX74c9fX1WL58ObZu3QoAqKurQ0VFBWpra7Fv3z5s3LgRTqdzwu9NRMHJ7pRwor0fX2h7Ea5UIC4iuDuDG7uMeOjPX+D3H59Cblostn1rEW6+ZEbQ96lM1Lgn5WVlZcHpdEKhUOCee+7B1772tQm9YX9/Pz766CO8+OKLAICwsDCEhYWhsrISH3zwAQCgrKwMV111FZ566ilUVlZi9erVUKvVyMzMRFZWFqqqqrBs2bIJvT8RBZ8+sx1Hda49J4J9gp3V4cSuqma8cUiL2HAVHrw2B/+UPS2oa/aGcYVFZGQkbDYbioqK8NBDDyEtLQ1Go3FCb3jq1CkkJyfjnnvuwRdffIHFixfjmWeeQXt7u3toblpaGjo6OgAALS0tuOyyy9zP12g0aGlpmdB7E1FwcUoCzT1GnOo0IlqtCvpVYg839+K3H5yErs+CkrxU3PO12YgJ9+/2poEyrmaol19+GZIk4bnnnkNUVBSam5vx5z//eUJv6HA4cPDgQWzYsAGHDh1CVFSUu8lpJEKI846NluDbt29HcXExiouL0dnZOaH6iMg/jFYHDjXr0dhlQmKUGhFhwTtqqM9sx3/97Th+VlkDuUyGJ28uwP3XZE+ZoADGeWcxa9YsmM1m6HQ6PProoxf1hhqNBhqNBkuXLgUA3H777di6dStSU1Oh0+mQlpYGnU6HlJQU9+Obm5vdz9dqtUhPTx/xtdevX4/169cDAIqLiy+qTiLyDSEEWnstONFuQLhKgWlBvEqsEALvH+/AC39vhNnmxB2XZmDV4oyQWLDQ28b1N37zzTdRVFSE66+/HgBw+PBhlJaWTugNp0+fjoyMDBw/fhwA8N577yEvLw+lpaUoLy8HAJSXl2PlypUAgNLSUlRUVMBqtaKxsRH19fVYsmTJhN6biALLYnfiiLYPJ9oNSIgMQ3QQL8mt6zPjkT21+PW79dDER+DpO4pw59JZUzIogAuYlFdVVYWrrroKAFBUVISmpqYJv+mzzz6LtWvXwmazYc6cOfjjH/8ISZKwatUq7NixAzNnzsRrr70GAMjPz8eqVauQl5cHpVKJbdu2uRczJKLQ0dFvwfE2A2QyWVDfTTicEnYfbsWuqjNQKmTYcOVcXF8wPejXofK1cYWFUqlEXFyc1960qKgI1dXV5x1/7733Rnz8pk2bsGnTJq+9PxH5j83h2nOirc+MuCCfYHei3YBn/7ceTd0mLJuThO9eMQdJQRxs/jTuhQRfffVVOJ1O1NfX4ze/+c2Eh84STYTBYofZ5oRcLoNCJnP93/01oJC5/jzZhy+GGr3RhjpdPyRJICmIh8SabA688vlpvHVEh8SoMGy6IReXzUkKdFlBZVxh8eyzz2Lz5s1Qq9VYs2YNrrvuOvzsZz/zdW1EAFzt3IfP9MIpCWDwWiMDIOD+I4bGzMllgEohh0ohh1Ihg0ouh0oph1Iug0ohR5hSDrnsnKCRy845FpwXtFDicEpo6jbhTI8RMWoVwsODt+m4qrEbz3/YgO4BG75RmIZ/WTYLkWHB25cSKDIx0tjUSaC4uHjEpi4KLZIkcKSlD0aLA7ERnocpSkJAkgQkMfi1EJAkwOn+moHja/0WO462uibYxUcG754TPUYbtn/UgE8aujErMRLfuyYL86cH//amnnQbrVg8K2HCw3pHu3aOGZ+eRjxxiXLyteYeE/RG27g7ROUyGeSKi7s4nR04doeAVTggmRk4nkiSgFZvQkOncXDPieBs65eEwP7aNpR/2gSbU8Jdl83CLZfMgDKI+1KCwZhh8dlnnyEjIwNr1qzB0qVLR5wgR+QrfSY7GjoH/H7RYeBcOJPNgWO6fvRbHEiIDAva9ZHO9Jjw3PsncVTXjwWaONx3VRbS4yMCXZZX2J0S+sx2xIQroVZ6v9lvzLBoa2vD3/72N+zatQuvvvoqvvGNb2DNmjXIz8/3eiFEZ7M6nKhp7UNMuCpoLzxj8UbgCPFVc5o7cCwYDCHX9wBX6MgEIAZTRyaTQRr8YOfrwBFCoG1wzwm1QhG0y3XYHBJeO9CM1w9oEaFS4AfLs3HN/JSg7XC/EE5JoM9ig0IuQ15aDJJjwn3yIWHMsFAoFLj++utx/fXXw2q1YteuXbjqqqvwyCOP4Pvf/77XiyECXBegk+0DEEJMmo1jJkImk0EhAxSQYaKn4WICx/0tmRg1cAasDnQP2BAfoQraZpyalj489/5JtPSacVVOMu79+hzEjaP/K9gJIdBvccAhSZiVFAlNQqRPhyV77PK3Wq14++23sWvXLjQ1NeH+++/Hrbfe6rOCiFp7LWg3WJHM8e0XzdeBI5cH7wS7AYsDf/y0EX+ta0dqrBqPl+Zj0cyEQJflFSabA0abA9Njw5E5Ldov62qNGRZlZWWoqanBihUr8Oijj6KgoMDnBdHUZrDYUd9hQGIIb2w/2XgjcPxJCIG/n+zC9o9Pod9sx22LZmD1pZNje1N3v4RaicUzExEX6b87pDGHzsrlckRFRbkeeFbbnhACMpkM/f39vq9wgjh0NvTYnRIOnnbtkMhx7jQRHf0WPP9hA6pP65GVEo3vX52FOcmhv2vd2f0S81JiMC1a7bPBCxMaOitJkk+KIRpJQ+cALHZn0A65pODllATePNKKVz4/DZkM+M4/ZeIbhekhOTjibP7ulxgLP75RUOjot6BVbw7a9m8KXg2dA3juf0/iZOcAimclYMNVc5ESEx7osi5aIPolxsKwoIAz2Ryo0/UjPjJsUgxlJP+w2J14teoMKg+3IDZChR9dPx+Xz00K+d+hQPZLjIVhQQHllATqdP1QKxVBvRopBZeDp/XY9sFJdBisuC5/Ou5eNhvR4aF9OTu7XyI/Pdan/RITEdpnl0JeY5cRAxZH0E7mouDSa7Lhhb834sMTndAkRGDrrYXIT/fe9gmBEEz9EmNhWFDAdA9YcbrHiGkMCvJACIH3jnbgD580wmx34ltLZuL2xZqgvKheCJPNgQGbA2lB0i8xFoYFBYTF7nT1U4QH76qkFBxae83Y9v5JHGnpQ15aLL53dRYyEiMDXdZFObtfojiI+iXGwrAgv5MkgeNtBsghm7L7GZNndqeEvxxqwZ/+cQZhCjm+d3UWSvJSQ/rDRbD3S4yFYUF+p9Wb0HMBy47T1HNM14/n3j+J0z0mfD1rGr7zT3OQGBW6s/pDpV9iLAwL8qs+kx0nO/y/7DiFBqPVgZc+P429X+qQFK3Gz76RhyWZiYEu66KEUr/EWBgW5Dc2h4RaXeguO06+9dmpbvzuwwb0mmy4aWE67lw6K2QvrEBo9kuMhWFBfiGEQH27AU5JIEYduhcA8r7uASv++6NT+OxUNzKnRWHTDbmYlxoT6LImLJT7JcbCsCC/0HHZcTqHJAT21ri2N3UKgXu+NhulC9ODdl8MT87ul5idFIUZCREh1y8xFoYF+ZzBYsfxdi47Tl853W3Ec++fxLE2A4oy4nHfVVmYHhe66zkF2zpOvsCwIJ+yOyXUtfYjMkzBfgqCzSHhT9XN+PNBLaLCFHigZB6unJccsus5Bes6Tr7AsCCfaugcgNnu5HIehCPaXmx7/yRa+yy4Zn4K1l2eGbLbmzolgT6zDQrF5OqXGAvDgnymo98CXa+ZQTHF9Zvt+OOnjXj3aAfS4sLxxMoCLMyID3RZEzLZ+yXGwrAgnzDZHDiq60dcBJcdn6qEEPjwRCde+HsjBqwOfHOxBndcmgG1MjTb86dCv8RYGBbkdU5J4KiuH2EKLjs+VbX1W/D8Bydx8EwvclJjcN/VWcicFhXosiZkKvVLjIVhQV53utsIg5eWHQ9XyRGmlEMSAnKZDDaHBIud2/0GK6ckUHm4BTurzkAhk+G7V8zBioK0kBzcMBX7JcbCsCCv6h6worHbO8uOh6vksDgk3PPiP6DVm6FJiMDzdy52HWdgBJ36dgOe++AkTnUasTQzEf965dyQXP9rKvdLjIVhQV5jsTtx1IvLjocp5e6gAACt3owNrxzA9n9ZjCPaXijkMihkMsjlMvfXisGv5ef82fU4DHsc+1Im7uw7PgjgjYNaPPdBA+IjwvDjFfOxbE5obm861fslxsKwIK8YWnZcJvPesuOSEO6gGKLVm2GwOLBpd81Fv75cBsjPCRSF4pzQOff77iDCyI8ZK8CGvRZGDrUx33P480cLSeVYzxt8fbkME76Yj3TH99RtC2CXBK6cl4IodehdVtz9EuFTu19iLKH3U6WgpNWb0GOyYlqUd2fhahIihgWGJiECsREqPHlLISRJwCkJOMXg/yUB6ayvh45Lwx4DOIX46rmjPu781z739W1OAckuTaAGAUl49TRNyPAwwjhDUoZN38jFT974ctgd34/+fAQ7712KfrMjwH+rC8N+ifFjWNBF6zPb0dA5gIRI77VPm2wO/OGTU3jqtgX40Z+PDOuzUMplKJwR+vsuSwKeA2vMYBt6/lBgffV6Y77OqMGGr547xuOiw5Uj3vEFQwCOF/slLhzDgi6KzSGhtrUP0WrvLTtusjnw2J5aHG834PK507Dz3qWQhKvZaLKMhpLJvmqKCjWxEcoR7/hC5a/CfomJYZTShAkhcLLDAKdTIFzlnX9wA1YHHqmsxYmOATx03XzkTI9Fv9mBAYsD/WbHpAiKUGdzSHj+zsXQJEQAgPuOz+YI7p+N3Smha8AKhUKGxTMTkZcex6C4ALyzoAlr67Ogvd/qteGRBosdj+ypRVOXEQ9fPx+XzUnyyuuSd1nsEsJV8pC542O/hHcwLGhCBqwOHGszIN5LC8H1m+34WWUNzvSY8OMVuSG/leZkZ7EHbzgMYb+EdzEs6II5nBLqWvsQGabwykY1fWY7frr7S7T0mvHTb+Rh8awEL1RJUxn7JbyPYUEX7FTXAEw27yw7rjfZ8NPdNWjrs+CRG/NRFKKrkVJwcM+XiOB8CW9jWNAF6ei3QKs3e2U5jx6jDZt2f4lOgxWP3pSHBZr4iy+QpiT2S/gew4LGzWxz4lhbP+K9sOx494AVm3bXoNtoxWM35aMgxOdNUGCwX8J/GBY0Lk5JoE7XB5X84pcd7zRYsWn3l+g12fF4aQHy0mK9VCVNJeyX8C+GBY3L6W4j+s2Oix4m295vwabdX8JgceDnK/MxfzqDgi4M+yUCg2FBHumNNjR1Gy+6Q7utz4Kf7P4SJpsDT6wsQHZqjJcqpKmA/RKBFbDGPafTiUsuuQQ33ngjAKCnpwclJSXIzs5GSUkJ9Hq9+7FbtmxBVlYWcnJysH///kCVPCVZ7E7UtvYhNlx1UcuOt/aa8eM3jsBic+KJlYUMCho3IQT6zHb0mm2YnRSFpZlJSIkNZ1D4WcDC4plnnkFubq77z1u3bsXy5ctRX1+P5cuXY+vWrQCAuro6VFRUoLa2Fvv27cPGjRvhdDoDVfaUcvay4xezb7JWb8KP3/gSNoeEzbcUICsl2otV0mRmsjnQZbQiMUqFpZlJmDUtih3YARKQs67VavH222/j3nvvdR+rrKxEWVkZAKCsrAy7d+92H1+9ejXUajUyMzORlZWFqqqqQJQ95bT0mtFjsiI2fOJtwmd6TPjJG19CkgSevKUQmdMYFOTZsHWcZnEdp2AQkLD4wQ9+gP/8z/+EXP7V27e3tyMtLQ0AkJaWho6ODgBAS0sLMjIy3I/TaDRoaWkZ8XW3b9+O4uJiFBcXo7Oz04d/g8mvz2xHfYcBCRET76c43W3ET974EgDw5C2FmJUU5a3yaJJySgI9RitMdgfy02OxKCMBcV5aUoYujt/D4q233kJKSgoWL148rscLcf4i+aON8V+/fj2qq6tRXV2N5OTki6pzKrM5XMt5xFzEsuONXQP4yRtfQiGX4clbCpGRGOnlKmkyYb9E8PP7aKhPPvkEe/bswTvvvAOLxYL+/n7ceeedSE1NhU6nQ1paGnQ6HVJSUgC47iSam5vdz9dqtUhPT/d32VPG0LLjDqdAtHpit/0nOwbwSGUN1Co5Nt9ciPT4CC9XSZPJ0HyJtLgIzE6KYnNTkPL7ncWWLVug1WrR1NSEiooKXHPNNXjllVdQWlqK8vJyAEB5eTlWrlwJACgtLUVFRQWsVisaGxtRX1+PJUuW+LvsKWNo2fH4yLAJPf9EuwE/rfwS4WEKbLllAYOCRnVuv0RuWiyDIogFzTyLhx9+GKtWrcKOHTswc+ZMvPbaawCA/Px8rFq1Cnl5eVAqldi2bRsUCv5C+cKA1YHj7RNfdvxYWz8e3VOLmHAlNt9ciNRY7+7HTZPD0HwJpULO+RIhRCZG6hSYBIqLi1FdXR3oMkKGwynh0Bk9JAFEhl34Z4g6XT8e21OL+EgVNt9ciOQY7+3HTZPD2es4ZSZFIZ3rOAWl0a6dQXNnQYE1tOx44gRmade09OHxt2qRFKXG5psLkOSlnfNo8mC/ROhjWBA6DRZo9RZMi7rwfooj2l78/K06JMeosfnmQiRO4DVo8hq2jtOsRA6DDWEMiynObHPiqK4f8RGqC152/HBzL/7j7TpMjw3HEzcXIGGCneI0+ZzbL5Eco77oZe0psBgWU5hTEjiq64NSLr/gtuMDp/XY/E4dZsRH4ImbC/mJkQCwX2IyY1hMYWd6jOibwLLj/2jqwZPvHMXMpEj8R2kBYhkUBPZLTHYMiylKb7ThVKfxgoPi81PdeGrfMcyeFoX/KC1AdDh/haY69ktMDfyXPgVZ7E7U6voQF3Fhy45/crILv/jrcWQlR+Ox0nxEq/nrM5kIISCJr74WAIQApMHR9UIAAsL1/8GvbU4JKvZLTAn81z7FCCFQ32GADBe27PjH9Z345V+PIyc1Bo+V5k9oLgZ95eyL8bCvz7kYu/4/+Jizvz7ruTjn+iwb/E8M/ic777uu55w7w0ouA+QyGeRymev/ABRyORQKGRRDxwEo5EOPAcJVCqTGhrNfYgrgv/gpRqs3o9NgRXL0+GdXf3C8A79+9wRy02LxyI15IRkU0uDFdujr0Y8BGLxQS2d/b+hCjeFXWJn7sgz3d2Sy4VdiMXjs7PmvMrguujKZDIqzL9KDF2G5zHWBlskBpVwOGQC5HFDI5JDLBi/YMplroUeZqw65zPU+MrieDxnOOyaTwfUfXF/Lh77HGdTkQej9q6cJ67fYcbJjAImR4++n+N9j7XjmvXrkp8fhkRvzEK7yf6ela4SNHQ5JDP+UPHhNlp194LzL+eAnbRkgl8mhkLsunkq5DDLIoFB8daFWyAE5ZIMXcddxpfuCLoNCIRt8reEXWpls5GNDTXzui/TZX7O5hkIMw2KKsDsl1LX0ISpMOe5lx/9W14Zn//ckFmbEY9MNuQEJCtcIGyc0CeGYFh3u+qQM1yfuoQu3/KxPykOfmoc+Uct5YSbyCobFFOBadnwANqdAQuT4Lvj7atqw7YOTWDQzHj+5IfeitlWdCLtTQp/FhphwFYpnJ1zUbn1EdPEYFlNAe58Fuj4zpo1z3ae3j7Tidx+dQvGsBPx4RS7ClP7rvBzaBAcyIDc1lhvgEAUJhsUkZ7Q6cKzdgISIsHE1x+z5ogW//7gRSzMT8aPr5/t1lIvJ5oDJ5sSMhHDMSory+90MEY2OYTGJOZwSjrb2I1ypgHIcF/03Dmnxh0+asGxOEh68LsdvQWF3Sug12xAbocJiNjkRBSWGxSR2qmsARptjXMuOv3agGS99dhpfz5qGB0rmjStcLtbZTU5509nkRBTMGBaTVFuvGdoe87iW86j4xxns/L8zuHJeMv7fP88b92ipi2GyOWC0OpCRGImZSZFsciIKcgyLSWjA6sDRNgMSIsfupxBC4NWqM6j4RzOuyUnB/cuzfR4UZ68jVJyZyCYnohDBsJhk7E4JNS19iAwbu59CCIGXPz+N1w5oUZKbivuuzvJpUEiDTU4yGZA7PYZNTkQhhmExiQghUN9ugN0hIX6MjYiEEHjx0yb85VALrsufjo1Xzb2gBQUvlNHqgMnmanKalRTl16G4ROQdDItJRKs3o63fiuQx+imEEHjh743Y80UrbihMw3evmOOzoBhqcoqNUKJ4BpuciEIZw2KS6DPZcbLDMOa6T0IIbP/oFN76UofShem49+uZPlkKY6jJSc4mJ6JJg2ExCVjsTnzZ2ototWrUfgdJCPzuwwbsrWnDzUUzsO7y2T4JCjY5EU1ODIsQJ0kCx9tc+1OMttCfJAS2vX8Sf61rx+2LNLhr2SyvB8XZTU6XzkhEDJuciCYVhkWIa+o2Qm+yIemciXfhKjnClHJIkkB7vxWdBivuuDQDa5fM9GpQnN3kxN3SiCYvhkUI6zJY0NhtPG+BwHCVHBaHhHte/Ae0ejM0CRF4+o4ixEWoYLFLXnt/o9UBs92BjIRIzGSTE9Gkxn/dIcpkc6BO14+EiLDzRjOFKeXY8MoBaPVmAK5RUj/402GvXcztTgmdA1aEqeQonp2IrNQYBgXRJMc7ixDkcEqoa+2HSqEYcbE/SQh3UAzR6s3ubUInShICvSYbFHIZCtjkRDSlMCxC0KmuARitoy8Q2D1ggyYhYlhgaBIicDGjV91NTolRmJkYyTsJoimG/+JDzNACgQmjzNB++0sdNr99FL9eVQRNQgQAV1A8f+di2BwX3l8x1OSkHmpySolmUBBNQbyzCCEGi33MBQI/OdmF//6wAZfOTkRMuBI7710KSQByGWBzSBfUuc0mJyI6G8MiRNgcEmpb+0ddIPCIthe//OtxzJ8egwevy4HdKWA3Oyb0XgODTU4z2eRERIMYFiHA0wKBpzoH8MTbR5EWH4Gf3Zg36uQ8T1w71tmREKlCASfWEdFZGBYhQKs3o90w8gKBbX0WPPZmLaLUSvy8NH9CF3h3k5OCTU5ENDKGRZDrM9lR324YceRTr8mGR/bUwOEU2Hxz/rh2xTsXm5yIaDwYFkFsaIHAmPDzFwg02Rx4/M06dBtt2LyyABmJkRf02nanhD6LHfERKhRqkhCt5q8CEY2OV4ggJUkCx0ZZINDulLBl7zGc6hrAT7+Rh/lpseN/3bOanPLT2OREROPDsAhSTd1G6I2285qWJCHw9LsncLi5Fz9Yno1LZyeO+zUHrA5Y7E5kJEayyYmILgjDIgiNtkCgEAIvfHwKH9V34e6vzcby3NRxvZ7NIaHfYkd8pAqFmjg2ORHRBeNVI8iMtUDg6we1ePOIDisXpuPWS2Z4fC1JCOhNNigVMi4fTkQXhWERRBxOCXUtIy8Q+G5dO1767DSunJeMdePYDtUpCfSYrJiZEIlZ06JGXHCQiGi8GBZBpKFzAEbb+QsEVjV249n363FJRjz+bXn2eXcc57I7JehNNuRMj4Em4cJGSRERjYRhESTaes1o0ZvP69A+quvHU/uOY25yNH68ItfjHcJQ/0ThjDikxIb7smQimkIYFkFgtAUCT3cb8fO36jAtOgyP3pSPiLCxl/Gw2J0w2Z0oyohHQtTIq9ISEU0EwyLARlsgsNNgxWNv1kKlkOHxlQWIixh7GQ+j1QGHJOGSmfGI5ZpORORlfu/1bG5uxtVXX43c3Fzk5+fjmWeeAQD09PSgpKQE2dnZKCkpgV6vdz9ny5YtyMrKQk5ODvbv3+/vkn3m7AUCI8O+yu1+sx2P7qmByebE46X5mO6hOanfbAcgsGhWAoOCiHzC72GhVCrxq1/9CkePHsXnn3+Obdu2oa6uDlu3bsXy5ctRX1+P5cuXY+vWrQCAuro6VFRUoLa2Fvv27cPGjRvhdDr9XbZPuBYItAxbSdZid+Lnb9Whrd+Cn34jD5nTosd8jV6zDWqlHEUzE4YFDhGRN/k9LNLS0rBo0SIAQExMDHJzc9HS0oLKykqUlZUBAMrKyrB7924AQGVlJVavXg21Wo3MzExkZWWhqqrK32V7Xa/J5logMPKrDm2HU8JT+47hRLsB/35tDgpnxI35Gj1GK2LCVViQET/hZcmJiMYjoIPvm5qacOjQISxduhTt7e1IS0sD4AqUjo4OAEBLSwsyMjLcz9FoNGhpaRnx9bZv347i4mIUFxejs7PT93+BCbLYnahp7Ru2QKAQAs+9fxLVp/XYcNVcfG3utFGfL4RAt9GKadFqFKTHctkOIvK5gF1lBgYGcNttt+Hpp59GbOzoC+EJIc47NtqEtPXr16O6uhrV1dVITk72Wq3eNNoCgS99dhrvHevAt5bMxIqCtNGfLwS6jDakx4cjNy12xF3ziIi8LSBXGrvdjttuuw1r167FrbfeCgBITU2FTqcDAOh0OqSkpABw3Uk0Nze7n6vVapGenu7/or1kaIHAszuiKw+34PWDWqwomI7Vl2aM+lyn5LqjmDMtEtkpMZDLuXQHEfmH38NCCIFvf/vbyM3NxQ9/+EP38dLSUpSXlwMAysvLsXLlSvfxiooKWK1WNDY2or6+HkuWLPF32V7RabCgqduIxLPmQHx4ohMv/L0Ry+Yk4btXzB31rsnulNBttCInNRazp0VzjSci8iu/D5/55JNP8PLLL6OwsBBFRUUAgCeffBIPP/wwVq1ahR07dmDmzJl47bXXAAD5+flYtWoV8vLyoFQqsW3bNigUodeZa7I5cFTXj/izFgg8eEaPp989gYL0WPz7tTnnbXA0hLOyiSjQZGKkToFJoLi4GNXV1YEuA4BrlNPhM72wS8K9PPiJdgM27f4SaXER2HJLIaJGWTbcYnfCaHOgKCN+2BBbIiJfGO3ayd5RHxNC4GTnAIx2hzsoWvRmPP5mLWLDVXjspvxRg8JodcDqcGLRrAQGBREFFMPCx9r6LND1mpEQ4brY9xhteGRPDQDg56UFw/ovzsZZ2UQUTDjl14cMFjuOtRkQH+FaINBodeDRPTXot9ix+eZCzEiIGPF5vWYbIlQKFMyI42Q7IgoKvLPwEZtDQm1LH6LClFAq5LA5JDzxdh2a9Wb8eEUu5qXGjPg896xsDWdlE1HwYFj4gHuBQEkgIkwBpyTwy78eR01rP36wPBuLZiaM+JyuASuSYzgrm4iCD5uhfKC9z4J2gxXJ0WoIIfC7Dxvw2alu3Pv1TFyVk3Le4yUh0G20QZMQjqxkTrYjouDDsPAyq8OJ+o4BxA/uP1Hxj2bsq23DbYs0WFk047zHD+2VPWdaFGYlRXGyHREFJYaFl53qNAIyQKWQY2+NDq9WncE181NQtmzWeY8d2it7/vTYUTu7iYiCAcPCi3pNNhgsdsxMiITZ7sTCjHisvlSDO4pnnnfHYHU4YbA4OCubiEICw8JLHE4JZ3pMiFApsHbH/0GrN0OTEIHn1y5CuEoBi11yP9Zsc8Jkd+CSmZyVTUShgUNuvKRFb0ZkmBIbdh6EVm8G4NoJb8POg8NGNhmtDticTizmrGwiCiEMCy8wWh041WWEQgZ3UAzR6s2QBlffMljsgAxYNCsBMZyVTUQhhGFxkYQQONFuQLhKgT6LA5pzOqo1CRGQy1z9GWFKOYoy4rlXNhGFHIbFRWrvs0BvsmPA4sCTb9fhV99c6A4MTUIEnr9zMboMVsRGqLCQe2UTUYjiR9yLMDSnIkatxGNv1qKxywiZTIad9y6FJFxJ3GGwIkqtQM702FH3qyAiCnYMi4swNKdib40OtYNLeUSoFOg3O1yzsgds0CRyVjYRhT6GxQT1mmzQ9Zlhsjjx8uensTQzEdfMdy3lMTQre25yNGYmRXJWNhGFPPZZTIDDKeFYmwERSgV+/d4JRKgU+N7VWe5Q0JutmDstGrOmcfkOIpocGBYT0KI3w2J3Ys8XrWjoNOK+q7Pccyb6zDZMi1YjIzEywFUSEXkPm6Eu0NCciu4BK/5U3YyrcpLxtbnTALj2y5bJZJiXyj4KIppceGdxAYbmVMhlMjzzXj0SIsPw3SvmAnD1UxisdhSkx0Gt5PBYIppcGBYXYGhOxRuHtGjWm3H/8mxEq103Z3qTDdkpMYiL5MxsIpp8GBbjNDSnokVvQuXhVqwomO7e8a7XbENKjPq82dtERJMFw2KcTnUaYbY78ez/nsT0uHCsuzwTgKufQiGXITs1hiOfiGjSYliMw9CcitcPaNFpsOIH/zwP4SrX3toDVgfy0+O4ZzYRTWq8wnkwNKeivn0Af61rx62LZiAvLRYA0GOyITs1GnER7KcgosmNYeFBi96MrgEr/vvDU5iVGIm1S13bo/aabEiNUWNGPPspiGjyY1iMYWhOxZ+qmtFnseP/lcyDSiGH2eaESsF+CiKaOhgWoxiaU3HoTC8+PtmFNZdmYG5yNJySgMnuQN4M9lMQ0dTBq90o2vssaOwy4o+fNGJeajRuX5wB4Kv5FLHc6Y6IphCGxQisDidOtBuw8//OwOqQ8IN/ngeFXAa9yYbUWDXS48MDXSIRkV8xLEZwqtOID+u7cOC0HmVfm4WMhEiYbA6EKeXspyCiKYkLCZ5FkgS6BqxQyGW4bE4i+kw23LggHQ6nBJPNieLZCVApmK9ENPUwLAZJksDxdgO+81I1tHozNAkR+O23FiFcJYe214zc6bGIYT8FEU1R/Jg8qNtocwcFAGj1Zmx89SAUcjnS4iIwPY79FEQ0dTEsBtkcTndQDNHqzZDJgKyUaPZTENGUxrAYFKZUnLdqrCYhApFhSvZTENGUx6vgoKSoMPz+rmJ3YGgSIvDfdy5GSow6wJUREQUeO7gHyeUy5KTG4C8bvgaDxQGVUoYZcRHcHpWICLyzGEYul2FatBoCAmlxEVCw+YmICADvLM4jl8uQlRIT6DKIiIIKPzoTEZFHDAsiIvKIYUFERB4xLIiIyKOQCYt9+/YhJycHWVlZ2Lp1a6DLISKaUkIiLJxOJ+677z7s3bsXdXV12LVrF+rq6gJdFhHRlBESYVFVVYWsrCzMmTMHYWFhWL16NSorKwNdFhHRlBESYdHS0oKMjAz3nzUaDVpaWgJYERHR1BISk/KEEOcdG2kV2O3bt2P79u0AgGPHjqG4uHhC79fZ2Ynk5OQJPdefWKf3hUqtrNP7QqVWX9fZ1NQ04vGQCAuNRoPm5mb3n7VaLdLT08973Pr167F+/fqLfr/i4mJUV1df9Ov4Guv0vlCplXV6X6jUGqg6Q6IZ6tJLL0V9fT0aGxths9lQUVGB0tLSQJdFRDRlhMSdhVKpxHPPPYfrrrsOTqcT69atQ35+fqDLIiKaMkIiLADghhtuwA033OCX9/JGU5Y/sE7vC5VaWaf3hUqtgapTJkbqPSYiIjpLSPRZEBFRYDEszhLMS4o0Nzfj6quvRm5uLvLz8/HMM88AAB577DHMmDEDRUVFKCoqwjvvvBPgSoHZs2ejsLAQRUVF7uHLPT09KCkpQXZ2NkpKSqDX6wNa4/Hjx93nrKioCLGxsXj66aeD5nyuW7cOKSkpKCgocB8b6xxu2bIFWVlZyMnJwf79+wNa54MPPoj58+djwYIFuOWWW9Db2wvANSQzIiLCfW7/9V//NaB1jvWzDqbzeccdd7hrnD17NoqKigAE4HwKEkII4XA4xJw5c0RDQ4OwWq1iwYIFora2NtBlubW2tooDBw4IIYTo7+8X2dnZora2Vjz66KPiF7/4RYCrG27WrFmis7Nz2LEHH3xQbNmyRQghxJYtW8RDDz0UiNJG5HA4RGpqqmhqagqa8/nhhx+KAwcOiPz8fPex0c5hbW2tWLBggbBYLOLUqVNizpw5wuFwBKzO/fv3C7vdLoQQ4qGHHnLX2djYOOxx/jRSnaP9rIPtfJ7thz/8oXj88ceFEP4/n7yzGBTsS4qkpaVh0aJFAICYmBjk5uaG1Cz2yspKlJWVAQDKysqwe/fuwBZ0lvfeew9z587FrFmzAl2K2xVXXIHExMRhx0Y7h5WVlVi9ejXUajUyMzORlZWFqqqqgNV57bXXQql0jZ257LLLoNVq/VLLWEaqczTBdj6HCCHwP//zP1izZo1fajkXw2JQKC0p0tTUhEOHDmHp0qUAgOeeew4LFizAunXrAt68A7hm11977bVYvHixe0Z9e3s70tLSALiCr6OjI5AlDlNRUTHsH2Cwnc8ho53DYP7d/cMf/oAVK1a4/9zY2IhLLrkEV155JT7++OMAVuYy0s86WM/nxx9/jNTUVGRnZ7uP+fN8MiwGiXEuKRJoAwMDuO222/D0008jNjYWGzZsQENDAw4fPoy0tDQ88MADgS4Rn3zyCQ4ePIi9e/di27Zt+OijjwJd0qhsNhv27NmDb37zmwAQlOfTk2D93d28eTOUSiXWrl0LwBVwZ86cwaFDh/Bf//Vf+Na3voX+/v6A1TfazzpYz+euXbuGfajx9/lkWAwa75IigWS323Hbbbdh7dq1uPXWWwEAqampUCgUkMvl+M53vuO32+WxDJ23lJQU3HLLLaiqqkJqaip0Oh0AQKfTISUlJZAluu3duxeLFi1CamoqgOA8n0NGO4fB+LtbXl6Ot956Czt37nRfaNVqNZKSkgAAixcvxty5c3HixImA1TjazzoYz6fD4cBf/vIX3HHHHe5j/j6fDItBwb6kiBAC3/72t5Gbm4sf/vCH7uNDFw8AeOONN4aNoggEo9EIg8Hg/vqvf/0rCgoKUFpaivLycgCuC8nKlSsDWabbuZ/Wgu18nm20c1haWoqKigpYrVY0Njaivr4eS5YsCVid+/btw1NPPYU9e/YgMjLSfbyzsxNOpxMAcOrUKdTX12POnDmBKnPUn3WwnU8AePfddzF//nxoNBr3Mb+fT791pYeAt99+W2RnZ4s5c+aIJ554ItDlDPPxxx8LAKKwsFAsXLhQLFy4ULz99tvizjvvFAUFBaKwsFDcdNNNorW1NaB1NjQ0iAULFogFCxaIvLw893ns6uoS11xzjcjKyhLXXHON6O7uDmidQghhNBpFYmKi6O3tdR8LlvO5evVqMX36dKFUKsWMGTPECy+8MOY5fOKJJ8ScOXPEvHnzxDvvvBPQOufOnSs0Go379/S73/2uEEKI119/XeTl5YkFCxaISy65ROzZsyegdY71sw6m8ymEEGVlZeL5558f9lh/n0/O4CYiIo/YDEVERB4xLIiIyCOGBRERecSwICIijxgWRETkEcOCyAt6e3vx29/+FgDQ2tqK22+/PcAVEXkXh84SeUFTUxNuvPFG1NTUBLoUIp8ImW1ViYLZww8/jIaGBhQVFSE7OxtHjx5FTU0NXnzxRezevRtOpxM1NTV44IEHYLPZ8PLLL0OtVuOdd95BYmIiGhoacN9996GzsxORkZH4/e9/j/nz5wf6r0XkxmYoIi/YunUr5s6di8OHD+MXv/jFsO/V1NTg1VdfRVVVFTZt2oTIyEgcOnQIy5Ytw0svvQTAta/ys88+iwMHDuCXv/wlNm7cGIi/BtGoeGdB5GNXX301YmJiEBMTg7i4ONx0000AgMLCQhw5cgQDAwP49NNP3SvfAoDVag1UuUQjYlgQ+ZharXZ/LZfL3X+Wy+VwOByQJAnx8fE4fPhwgCok8ozNUEReEBMT415t90LFxsYiMzMTr732GgDXCsNffPGFN8sjumgMCyIvSEpKwuWXX46CggI8+OCDF/z8nTt3YseOHVi4cCHy8/ODaktfIoBDZ4mIaBx4Z0FERB4xLIiIyCOGBRERecSwICIijxgWRETkEcOCiIg8YlgQEZFHDAsiIvLo/wNAD//Uw6+yCQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(facecolor = 'w',figsize=(6,6))\n", "sns.lineplot(data=df[(df['concentration']==500000)|(df['concentration']==0)], \n", " x='time', y='Mean', style=\"label\",ci='sd', markers=True, dashes=False)\n", "plt.ylim(0,1200)\n", "plt.savefig('time_scale.pdf')" ] }, { "cell_type": "code", "execution_count": 30, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgEAAAF7CAYAAABLt1hbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABbAUlEQVR4nO3de3xU1bk//s+e+0wyuZKEhAEJJHIJCEJE8VZbigdbC61FCrUVb+Uotra1l+P5+vr+er59HQvtaY+1HqoHSxWtmlPbHrGtxWqrrYAUg1dAIAQCud8z99vee/3+mGRIyMzktpOZST7vfyA7MzuLLbKeWc9azyMJIQSIiIhoytElewBERESUHAwCiIiIpigGAURERFMUgwAiIqIpikEAERHRFMUggIiIaIoatyDgjjvuQGFhIRYtWjToez/+8Y8hSRI6Ojqi17Zt24aysjLMmzcPr7zySvT64cOHsXjxYpSVleG+++4DTzQSERFpY9yCgNtuuw179+4ddL2+vh6vvvoqZs2aFb127NgxVFVV4ejRo9i7dy+2bt0KRVEAAPfccw927tyJmpoa1NTUxLwnERERjdy4BQHXXnst8vLyBl3/5je/iR/96EeQJCl6bc+ePdi4cSPMZjNKS0tRVlaGQ4cOobm5GS6XCytXroQkSbj11lvx4osvjteQiYiIppQJ3RPw0ksvYcaMGViyZMmA642NjZg5c2b0a4fDgcbGRjQ2NsLhcAy6TkRERGNnmKgf5PP58NBDD+HPf/7zoO/FyvNLkhT3ejw7d+7Ezp07AQDHjx/H/PnzxzBiIiKi9FFXVzdgr91wTFgQUFtbizNnzkRXARoaGrBs2TIcOnQIDocD9fX10dc2NDSgpKQEDocDDQ0Ng67Hs2XLFmzZsgUAUFlZierq6nH60xAREaWWysrKEb9nwtIBixcvRltbG+rq6lBXVweHw4F33nkH06dPx9q1a1FVVYVgMIgzZ86gpqYGK1asQHFxMex2Ow4ePAghBJ5++mmsW7duooZMREQ0qY1bELBp0yasXLkSJ06cgMPhwK5du+K+tqKiAhs2bMDChQuxZs0a7NixA3q9HgDw2GOP4a677kJZWRnmzp2LG264YbyGTERENKVIk7WVMNMBREQ0lYxm3puwPQFEREQTKRwOo6GhAYFAINlD0ZTFYoHD4YDRaBzzvRgEEBHRpNTQ0AC73Y7Zs2cnPFmWToQQ6OzsRENDA0pLS8d8P/YOICKiSSkQCCA/P3/SBABA5Jh8fn6+ZqsbDAKIiGjSmkwBQB8t/0wMAoiIiOLo6enBz3/+cwBAU1MT1q9fn+QRaYtBABERURz9g4CSkhL85je/SfKItMWNgURERHE88MADqK2txdKlS1FeXo6PPvoIR44cwVNPPYUXX3wRiqLgyJEj+Na3voVQKIRnnnkGZrMZL7/8MvLy8lBbW4t7770X7e3tsNlseOKJJ1KqpD1XAoiIiOLYvn075s6di/feew//8R//MeB7R44cwXPPPYdDhw7hwQcfhM1mw7vvvouVK1fi6aefBhApZ//oo4/i8OHD+PGPf4ytW7cm448RF1cCiIiIRuHjH/847HY77HY7srOz8ZnPfAZApEz+Bx98AI/HgwMHDuDmm2+OvicYDCZruDExCCAiIhoFs9kc/b1Op4t+rdPpIMsyVFVFTk4O3nvvvSSNcGhMBxAREcVht9vhdrtH9d6srCyUlpbihRdeABAp9PP+++9rObwxYxBAREQUR35+Pq666iosWrQI3/nOd0b8/meffRa7du3CkiVLUFFRgT179ozDKEePDYSIiGhS+uijj7BgwYJkD2NcxPqzjWbe40oAERHRFMUggIiIaIpiEEBERDRFMQggIiKaohgEEBERTVEMAoiIiKYoBgFERETj5I477kBhYSEWLVoUvfbee+/hiiuuwNKlS1FZWYlDhw5Fv7dt2zaUlZVh3rx5eOWVV8Z9fAwCiIiIAPzpT3/CjTfeiMsuuww33ngj/vSnP435nrfddhv27t074Np3v/tdfO9738N7772H73//+/jud78LADh27Biqqqpw9OhR7N27F1u3boWiKGMeQyIMAoiIaMr705/+hIceeggtLS0QQqClpQUPPfTQmAOBa6+9Fnl5eQOuSZIEl8sFAHA6nSgpKQEA7NmzBxs3boTZbEZpaSnKysoGrBKMBzYQIiKiSe8nP/kJTpw4Eff7H374IcLh8IBrgUAA3//+9/G///u/Md8zb948fOtb3xrxWH7605/in/7pn/Dtb38bqqriwIEDAIDGxkZcccUV0dc5HA40NjaO+P4jwZUAIiKa8i4MAIa6PhaPPfYYHn74YdTX1+Phhx/GnXfeCSDSYOhCkiRp/vP740oAERFNekN9Yr/xxhvR0tIy6Pr06dOxc+dOTceye/duPPLIIwCAm2++GXfddReAyCf/+vr66OsaGhqiqYLxwpUAIiKa8u69915YLJYB1ywWC+69917Nf1ZJSQn+9re/AQD++te/ory8HACwdu1aVFVVIRgM4syZM6ipqcGKFSs0//n9cSWAiIimvBtuuAEAsGPHDrS2tqKoqAj33ntv9Ppobdq0CW+88QY6OjrgcDjw//7f/8MTTzyBr3/965BlGRaLJbrSUFFRgQ0bNmDhwoUwGAzYsWMH9Hr9mP9sibCVMBERTUpsJTw0pgOIiIimKAYBREREUxSDACIioimKQQAREU1ak3Hbm5Z/JgYBREQ0KVksFnR2dk6qQEAIgc7OzkHHGUeLRwSJiGhScjgcaGhoQHt7e7KHoimLxQKHw6HJvRgEEBHRpGQ0GlFaWprsYaQ0pgOIiIimKAYBREREUxSDACIioilq3IKAO+64A4WFhVi0aFH02ne+8x3Mnz8fl1xyCT73uc+hp6cn+r1t27ahrKwM8+bNwyuvvBK9fvjwYSxevBhlZWW47777JtUuTyIiomQatyDgtttuw969ewdcW716NY4cOYIPPvgAF198MbZt2wYAOHbsGKqqqnD06FHs3bsXW7duhaIoAIB77rkHO3fuRE1NDWpqagbdk4iIiEZn3IKAa6+9Fnl5eQOuXX/99TAYIgcSrrjiCjQ0NAAA9uzZg40bN8JsNqO0tBRlZWU4dOgQmpub4XK5sHLlSkiShFtvvRUvvvjieA2ZiIhSQFBWkj2EKSNpewJ++ctfRls0NjY2YubMmdHvORwONDY2orGxccBZyL7r8ezcuROVlZWorKycdOdCiYimAlcgjHfP9SAQZiAwEZISBDz00EMwGAy45ZZbAMQugShJUtzr8WzZsgXV1dWorq5GQUGBdgMmIqIJ0e0JoccXwokWN1SVe8DG24QXC9q9ezf+8Ic/4C9/+Ut0Qnc4HKivr4++pqGhASUlJdFqTxdeJyKiyUcIgWZXAIV2C7p8QTQ5/XDk2pI9rEltQlcC9u7dix/+8Id46aWXYLOd/w+7du1aVFVVIRgM4syZM6ipqcGKFStQXFwMu92OgwcPQgiBp59+GuvWrZvIIRMR0QTxhRQEQgqMeh1yrWbUtHrgCcrJHtakNm5BwKZNm7By5UqcOHECDocDu3btwle/+lW43W6sXr0aS5cuxd133w0AqKiowIYNG7Bw4UKsWbMGO3bsgF6vBwA89thjuOuuu1BWVoa5c+dG9xEQEdHk0uMLoS/jq9dJsJn0ON7kgsK0wLiRxCQ9eF9ZWYnq6upkD4OIiIapuq4LQgAWoz56rdMbxKw8G+YUZCZxZOlhNPMeKwYSEVHSBcIKPEF5QAAAALk2E+o6fejxhZI0ssmNQQARESWdM84kr5Mk2M0GHGt2ISSrEzyqyY9BABERJV2LKwDrBasAfSxGPWRFoLbdM8GjmvwYBBARUVKFZBXdvnDcIAAAcqxGtDj9aHMFJnBkkx+DACIiSipXIAwgcTE4SZKQbTXhRIub1QQ1xCCAiIiSqs0VgNkQfxWgj1GvgyRJONHiZkdZjTAIICKipJEVFR2eEGymoYMAAMi2GtHlDaGph2kBLTAIICKipHEHZKhCQJcgFXChXJsJJ1vdrCaoAQYBRESUNO2eIIy6kU1Fep0Eq1GP482sJjhWDAKIiCgpVFWgzRVAhnnkvewyzAZ4gjLOdXnHYWRTB4MAIiJKCndQhqwK6HXDTwX0l2sz4Uy7F05fWOORTR0MAoiIKCm6vEHoR7AX4EI6SYLdYsSxZifCCqsJjgaDACIimnBCCLQ4R5cK6M9i1CPMaoKjxiCAiIgmnDekICirMOrHPg3lWI1o6gmg3c1jgyPFIICIiCZctzc0omOBiUiShByrEcdZTXDEGAQQEdGEa3H6h10gaDiMeh10koSaNlYTHAkGAURENKH8IQXeoDKsUsEjkWUxot0dQjOrCQ4bgwAiIppQTn8I0CYTMEiu1YiTbW54WU1wWBgEEBHRhGp2BmAzju1UQDwGvQ5mgx7HW1hNcDgYBBAR0YQJygp6/GFYjOM3/WSaDXD5ZTR0+cbtZ0wWDAKIiGjCuPwyICI7+sdTrs2E2nYPnH5WE0yEQQAREU2YVlcAVmPsDYEWow5ZVgMyLXpkWQ1jWi3Q6yLVBD9qdrGaYALjk5QhIiK6gKyo6PQGkWs1DfqexahDQFZx+1Nvo6HbD0euFY99aXnkenh0k7jFqEeXN4gzHR5cXJQ11uFPSlwJICKiCeEKyBBxUgEmgw73/OowGrr9AICGbj/u+dVhmAxjm6ZybSbUd/nRwWqCMTEIICKiCdHuDsAUp0ywKkQ0AOjT0O3HWDf4R6oJmlhNMA4GAURENO4UVaDdHYTNFDsLrZMkOHKtA645cq0YZZfhAUwGHSRJwilWExyEQQAREY07T0CGrAro48zqIVnFT25eEg0EHLlWPLrpUoRkbTb19VUTbHEyLdAfNwYSEdG46/AEYdDF/9zZ4Qnhx68cx39uWIqCTBPOdvnwH68cxz9fW4a8jMEbCUcjx2rEiVY3sm3GuCsSUw1XAoiIaFwJIdDiCiAjQcOgQ2c6cfhcD+q7fPAEFfiCCqrrevDIX2o0W8I36HUw6/X4qNkFldUEATAIICKiceYJygjLKgxxNgUCwJs1HSi0m1FemAkAmJlnwx1XzcY757rx8ofNmo0l02KAOyCjntUEATAIICKicdbtDcXdCwBEgoT36ntwVdm0AccHP7W4GMtm5eKX++s0nbRzrCbUdnjgCrCaIIMAIiIaV83OQMIc/D9Od0JWBa4umzbguiRJ+MaqcliMOvz41ROaVf7T6yRkmoz4qMkFeYpXE2QQQERE48YXkuEPKwmL/uw7NTAV0F9uhglf/UQ5Trd78dw/zmk2LqtJj4Cs4EyHV7N7piMGAURENG56vImX3D2BSCrg6gtSAf2tnJOP6xcW4bfvNOBIo1OzseVaTajv9qHTE9TsnumGQQAREY2bZlcAGQlSAQfPRFIBV12QCrjQXVfPwfRsC/7ztZPwBmVNxiZJErIsRnzU4kJQnprVBBkEEBHRuAiEFbj8YVjidA0EEqcC+rOa9Lh/9cXo9ATx+N9rNRuj2aAHBHCq1TMlqwkyCCAionHh8ocRZ4UfwPlUwDXl8VMB/c2fnoUvVM7EGyfa8WZNu2bjzLaa0OYOonUKVhMctyDgjjvuQGFhIRYtWhS91tXVhdWrV6O8vByrV69Gd3d39Hvbtm1DWVkZ5s2bh1deeSV6/fDhw1i8eDHKyspw3333TclIjYgoHbW4ArAmWAU4eLoTiipw1dzEqYD+NlTOxLwiO3a8cQodGubys3urCfpDUystMG5BwG233Ya9e/cOuLZ9+3asWrUKNTU1WLVqFbZv3w4AOHbsGKqqqnD06FHs3bsXW7duhaJE/kPcc8892LlzJ2pqalBTUzPonkRElHrCiooubyhhEPDmqQ4UZZlRNkQqoD+DXof7V18MRRV4+LWTUDX6YGjU62DS63G8ZWpVExy3IODaa69FXl7egGt79uzB5s2bAQCbN2/Giy++GL2+ceNGmM1mlJaWoqysDIcOHUJzczNcLhdWrlwJSZJw6623Rt9DRESpy+WPnAqIt8zvDoTxfkPiUwHxlORYcdfVc/BBgxMvvdc05rH2ybQY0OMLo6F76lQTnNA9Aa2trSguLgYAFBcXo62tDQDQ2NiImTNnRl/ncDjQ2NiIxsZGOByOQdeJiCi1tbmDMOuHTgVcXVYwqvtfv7AIl5fmYfdbdZqe9c+1mVDb7p0y1QRTYmNgrDy/JElxr8ezc+dOVFZWorKyEu3t2m0aISKi4VNUgQ5PEDZzolMBnZieZcHcgoxR/QxJkvC1T5Qj02LAT/58QrOWw3qdBJtJP2WqCU5oEFBUVITm5kgjiObmZhQWFgKIfMKvr6+Pvq6hoQElJSVwOBxoaGgYdD2eLVu2oLq6GtXV1SgoGF10SUREY+MOhKGoAro4H9pc/kgq4MJeASOVbTXi66vKcbbLh6ffqhv1fS5kMxkQCCuomwLVBCc0CFi7di12794NANi9ezfWrVsXvV5VVYVgMIgzZ86gpqYGK1asQHFxMex2Ow4ePAghBJ5++unoe4iIKDV1eIIw6OJPLwfP9KUChn8qIJ7Ki/Lw6cXF2PN+E96r7xnz/frk2kw42+1Dtzek2T1T0bgFAZs2bcLKlStx4sQJOBwO7Nq1Cw888ABeffVVlJeX49VXX8UDDzwAAKioqMCGDRuwcOFCrFmzBjt27IC+N5f02GOP4a677kJZWRnmzp2LG264YbyGTEREY6SqAi3OIDLN8asE7j/VgeLs0acCLnTblbPhyLXip6+dhFujXL4kSci2GHG02TmpqwlKYpIevK+srER1dXWyh0FENKW4AmG8c7Yb+Rnm2N/3h/HlX/4DN13qwOYrZ2v2c0+1efDt37yPK+bk41/+ad6Y0gz99fhDyM8wYUFxlmb3HC+jmfdSYmMgERFNDl2eUNy9AEAkFaAK4OrysacC+isrzMQtK2Zh/6kOvH5Cu43hOVYTWl0BtLsnZ5MhBgFERKQJIQRanIGEqYB9NZFUwJxp2qQC+rtpmQMLi7Pw+N9q0erSrgRwttWE4y2uSVlNkEEAERFpwhdSEAgrMOpjTy1O/+gLBA2HXifh/tUXAwD+89WTUDSq/GfU62DQ6XCi1T3pqgkyCCAiIk10+0IJGwYdPN2bCtDgVEA8RVkW3P2xuTjW7MLv3mkY+g3DZLcY0e0LorHHr9k9UwGDACIi0kSLM4CMRKmA3lMBpeOQCujv4/MKcHXZNDx76BxOtXk0u2+u1Yzado9mJxBSAYMAIiIas0BYgScow2yIXSXQ6Q/jg3FMBfQnSRK2XjcXOVYjfvLqCQTC2uTy9ToJVqMex5snTzVBBgFERDRmTl/iojoTkQroz24x4pufvBgN3X48daBOs/vaTAZ4QwrqOidHkyEGAURENGYtrgBsxsSpgJIJSAX0t2RmDtYtKcEfP2xGdV2XZvfNtZlwttOLniECn3TAIICIiMYkJKvo9oVhMcY/FfBBQw+uLi8YVipA1bCG3a0rZ+OiPBse+WsNnH5tcvk6SUKWxYijTS7NGhclC4MAIiIak762u/Em+Ldq+1IB+UPeyx9S0OwMIKxRzt1k0OFb18+DJyDj0b/WxOxOOxoWox6qKlDbrt3Gw2RgEEBERGPS5grAEmdDIADsO9WOGTlWzM4fOhXgDckozbfB6Q9pNmGXTsvA5pWz8Y8zXfjzsVZN7glEuhg29/jRpmFhoonGIICIiEZNVlR0eEKwmuKfCviw0TmsUwGqEJAkoLQgE0VZFs2W7wFg7dISXOLIxhNvnkaTRmf9JUlCji1STVCrEwgTjUEAERGNmjsgQxUibr+AA7UdUAVw1TBOBfiCCgrtZpgMOpQV2qHTSZpNrjpJwjc/eTEMeknzaoJ6SYcTLelZTZBBABERjVq7JwhTnDLBQORUQCQVYBvyXkFFQXG2FUAkl19RkgV3UNZso+C0TDPuva4MJ1rd+HV1vSb3BIAsqxFdviCanOlXTZBBABERjYqqCrS5ArCZYh8N7PGFcGSYqQBFFTDoJGRZjdFrOTYTZufb0K3hUbxrygvw8XkFqHr7HI63uDS7b67VjJpWDzxBWbN7TgQGAURENCruoAxZFdDr4pwKGEGBIE9QxvRsy6B7XZSfgUyTAV4NJ9d/vnYupmWa8Z+vntSsM6BeJ8Fm0uN4k0uzVMNEYBBARESj0ukJwqBLnApw5Fpx0TBSAWFFRYHdMui6XidhQUkW/GFFs8k1w2zANz95MVqcAfxi32lN7glEqgl6QjLOdno1u+d4YxBAREQjJoSIVAmMcyqguzcVcNUwUgFhRYXZqEOWJXZaIcNswLwiO7o0TAssmpGNzy9z4M/HWvHW6U7N7ptrM6Gu05c21QQZBBAR0Yh5QwpCYRXGOJsC+woEXTOMVIA3KMORY00YLBTnWFBgN0ULE2nhi5fPwpyCDDz61xp0ebWZtHWSBLvZgGPN6VFNkEEAERGNWLc3FHcvAADs700FzMobOhUgqwL5meaEr5EkCeWFdgghNJtcjXodvr16HoJhFY/8RdtqgrKSHtUEGQQQEdGINff4454K6PaGcKRpeKcCAmEFmWY9Mszxmw/1sRj1mD/djh4NqwnOzLPhjqtm451z3Xj5w2ZN7gkAOVYjmp1+tLtTu5oggwAiIhoRf0iBL6zAZIg9hRwYwakAX0jGjJyhVwv6TLNbMDPPim6/djn3Ty0uxrJZufjl/jrUd2nTIliSJORYTTje7E7paoIMAoiIaESG2vS2/1QHZuZacdEQvQKEEFAFkJdpGtHPnzMtExaDHr6QNscGJUnCN1aVw2LU4SevntCseZFRr4MkSTjR4tZs5UJrDAKIiGhEml0BZCRKBfQWCBqKP6wgL8MIizF+86FYDHodFpRkwRfS7thgboYJX/1EOWrbvXj+0DlN7glEmgx1eUNo6knNtACDACIiGragrMDpC8eduA+c7oTA8HoF+MPnywSPVJbFiLkFGZpWE1w5Jx+rFxbhN4cbcLTJqdl9c20mnGx1a1rwSCsMAoiIaNicvjAS7fXbV9OOmXm2IVMBff0AcjNGlgroz5FrQ26GCW4Njw1+5eo5mJ5twX++elKzSVuvk2A16vFRc+pVE2QQQEREw9bqDsIaZxWgyxvC0SbXsGoD+IIKirIscesMDIdOJ2FekR2yqmqWx7ea9Lh/9cXo8ATx33+v1eSeQKTgkSco41xXalUTZBBARETDElZUdHriBwFv1XYMOxUQVBRMzxpcJnikrCY95k/PQo9Pu2OD86dn4QuVM/H6iXa8WdOuyT2BSFrgTLsXTp92KxdjxSCAiIiGxR2ILI/HO/v/5qkOzMqzDVkgKFbHwLEozLKgOMcKp1+7yXVD5UzMK7Jjxxun0OEJanJPnSTBbjHiWLNTs5WLsWIQQEREw9LmCsCsj58KONbkGlPHwLGYW5AJvV7S7Ey+Qa/D/asvhqIK/PS1k9E9DGNlMeoRktWUqSbIIICIiIakqALtniCscRoGHehNBQwnCJBVFYUapAL6Mxl0qCjOhjsoazZhl+RYcdfVc/B+gxMvvd+kyT2BSFqgqSeQEtUEGQQQEdGQPAEZiirifnrfd6oDF+XZMHOIVEBYUWE26GAfRpngkcq2GVE6zaZZMyAAuH5hES4vzcPuA3Wo69BmU58kScg0GdDY49fkfmPBIICIiIbU7gnAoIs9ZXR6gjjW5BrWhkBPUMaMIToGjsWsvAzYrZGd+FqQJAlf+0Q5Mi0G/OTVE5o1L5IkIBWKCDIIICKihIQQaHEGkRE3FdA57FSAogpMsyfuGDgWep2EBdOzEJQVyBptvsu2GvH1T5SjrtOHZw7WaXLPVMEggIiIEvIEZciqCkOcM/37TnVgdv7QqYBAWIHdYojbfVArGWYD5hXZ0aPhaYHK2Xn41OJivPheE96v79HsvsnGIICIiBLq8oSgi7N83+kJ4qPm4aUCIh0DR1cmeKSmZ1tQaDfDqWG3wduvnI0ZOVY8/NpJTasUJhODACIiSqjZFUBmnI18+2uH1yugr2PgWMoEj4QkSSgrygSkSL8DLViMenz7+nno8Yfx8zdqU7Yz4EgwCCAiorh8IRmBkBK3vO/+vlRAbuJUgD+sIC/TNOKOgWNhNuixsDgbzkBYswm7rDATt6yYhX2nOvDGSe2qCSZLUoKAhx9+GBUVFVi0aBE2bdqEQCCArq4urF69GuXl5Vi9ejW6u7ujr9+2bRvKysowb948vPLKK8kYMhHRlNTtDcVtGNTpCeJY8/AKBPnDyoSlAvrLyzDholybpt0Gb1rmwMLiLDz+t1q0upJ/1n8sJjwIaGxsxM9+9jNUV1fjyJEjUBQFVVVV2L59O1atWoWamhqsWrUK27dvBwAcO3YMVVVVOHr0KPbu3YutW7dCUbRZ2iEiosSanYG4G/n213YAGDoVoAoBSYrssk+G2dMyYDXp4Qtp1xXw/tUXQwjg4ddOplxnwJFIykqALMvw+/2QZRk+nw8lJSXYs2cPNm/eDADYvHkzXnzxRQDAnj17sHHjRpjNZpSWlqKsrAyHDh1KxrCJiKaUQFiBJyjHXcLfd6oTs/NtcAyRCvAGZRTax9YxcCwMeh0WFGfBH1Y0m7CLsiy4+2NzcLTJhd+926DJPZNhwv+LzJgxA9/+9rcxa9YsFBcXIzs7G9dffz1aW1tRXFwMACguLkZbWxuAyMrBzJkzo+93OBxobGyc6GETEU05zgRL6B29pwKuLi8Y8j4hRdWkY+BY2C1GzC3IRLdfm2ZAAPDxeYW4qmwanv3HOZxqS41eACM14UFAd3c39uzZgzNnzqCpqQlerxe/+tWv4r4+1maOeJWmdu7cicrKSlRWVqK9Pf03bBARJVOLKxC3bfD+U5FUwNVzE6cC+joGJisV0N+MHCvybGa4NKofIEkS7r1uLrKtRvzk1ROaNS+aSBMeBLz22msoLS1FQUEBjEYjbrrpJhw4cABFRUVobm4GADQ3N6OwsBBA5JN/fX199P0NDQ0oKSmJee8tW7aguroa1dXVKCgYOjolIqLYQrKKbl84YRBQOi0DM3ITb/bzBGUUZ1uh07Bj4GjpdBLmTbdDFUKzVr52ixHf/OTFaOj246kDdZrccyJNeBAwa9YsHDx4ED6fD0II/OUvf8GCBQuwdu1a7N69GwCwe/durFu3DgCwdu1aVFVVIRgM4syZM6ipqcGKFSsmethERFOKq7cYTqyV13Z3EB+1uId1KiCsqCjIGr8ywSNlMeoxv9iOHn9Is2ODS2fmYN2SEvzxw2ZUn+3S5J4TZXxrN8Zw+eWXY/369Vi2bBkMBgMuvfRSbNmyBR6PBxs2bMCuXbswa9YsvPDCCwCAiooKbNiwAQsXLoTBYMCOHTugj9PPmoiItNHmCsBsiN82GBi6V0BYUWE16selY+BYFNgtKMkJoc0VRK5Nm+JFt66cjffqe/Czv9Tg0U3LUiL9MRySmAwlj2KorKxEdXV1sodBRJR2ZEXFgdpOZFuNMcsFf+c37yMkq3hk46UJ79PtC6F0WsaQPQWSIayoOFzXBb1Op1kBozMdXtz/6/dQOTsX/+eGBQk7JQbCCgx6CZfOytXkZwOjm/dYMZCIiAZwB2SoQsQMANrdQRwfZipAFQL5mRNTJnikjHodFpRkwx0Ia3ZssHRaBm5deREOnu7Cqx+1anLP8cYggIiIBmj3BGHUxSkTPMwCQYGwgkzz+HcMHItsa++xQQ2rCa5bOgOXzMjGE2+eRlOPX7P7jhcGAUREFKWqAm2uADLi5PH31XRgTkEGSoYoAewNyXAMcXIgFczMsyHLaoAnoE01QZ0k4RufvBh6nYT/fDX1qwkyCCAioih3UIasCuhjHOlrcwdwotU9ZG2Avq1mE9UxcCx0OgkLirMQVBTIGh0bLLCbce91ZTjR6savq+uHfkMSMQggIqKoLm8Q+jgb2g6c6gQwdCrAH1aQl2GKe7og1dhMBswvsqPbr11a4JryAlw3rwBVb5/D8RaXZvfVGoMAIiICEPkE3+JMkAo41YG5w0gF+MPKkK9JNUXZFkzPsqBHw0Dg7mvnIj/TjP989ST8odSsJsgggIiIAADekIKgrMZs9NPmiqQChtMxUCdJyEmTc/J9JEnC3MJMSBIQlLWZsDPMBtz/yYvR4gxg177TmtxTa8PetnngwAHU1dVBls9vnrj11lvHZVBERDTxur2hmMcCgfOnAoY6GugNyijKMsOQpI6BY2E26FFRnI1367uRn6GL+yxGYtGMbHx+mQO/eacBlbPzcMWcfA1Gqp1hBQFf/vKXUVtbi6VLl0ar9UmSxCCAiGgSaXH6YTPFaxvcgbKCTBRnJ17mDykqipLcMXAscjNMmJWXgYZuH/IztCl3/MXLZ+Gd+m48+tcazCuyp9SGyWEFAdXV1Th27FjC6kdERJS+/CEF3qCC/MzBE1+rK4CTrR5sXjk74T36OgZmWdIrFXCh0mkZ6PGG4AvJmtQ5MOp1+Nbqefjm/7yHR/5ag+/duFCDUWpjWOs1ixYtQktLy3iPhYiIksTpDyHe57xo2+AhUgHuQDhlOgaOhV4nYX5JFnwhRbNz/rPybLj9qtk4fLYbLx9Jnfl0WCFOR0cHFi5ciBUrVsBsPh8lvvTSS+M2MCIimjjNzgCsxthTwv7aSCpgenbiZX5FFShMoY6BY5FpNqC8KBMnW90oyNQmvfHpxcV4u64bh+u68LlLS2DS69DuDiI/w5S0wGlYQcC//du/jfMwiIgoWYKygh5/GPkxOur1pQJuu3J2wnuEFRUWkx6ZKdYxcCxm5FjR5Q3BFQhrkuKQJAkPfmo+goqKO3dXo6HbD0euFU/cWol5RfakBALD+q/1sY99bLzHQURESeLyy4BAzH1ffamAoY4GugNhzCnInFR7xyRJwsVFdrxd14WQrMJkGPuJh2l2M275xT/Q0B3pK9DQ7cdXnq7G/269CgX2iV9FGdaf6ODBg7jsssuQmZkJk8kEvV6PrKys8R4bERFNgFZXANY47XT3nepAWWEmpg+x418VAtNibCpMdxajHvOn2+EMhKPlkMdCFSIaAPRp6PYjpFFtgpEaVhDw1a9+Fc8//zzKy8vh9/vxi1/8Al/96lfHe2xERDTOZEVFpzcIa4yjgS2uAGraPLhmGB0Ds6ymmPeYDArsFjhyLZp0G9RJ0qDGSo5cK0xJKrE87LWNsrIyKIoCvV6P22+/HW+88cY4DouIiCaCKyBDCMQsjNOXCrhyqAJBIRmOnPStDTAcpdMyYTbqx1z+NySreOxLy6OBQN+egPwk1Q4Y1p4Am82GUCiEpUuX4rvf/S6Ki4vh9XrHe2xERDTO2t0BmOJU99t3qgPlQ6QChBAQAHJSqADOeDDqdVhYkoXqui6YDLqYXRaHIxBWYTHq8MydKwABZFqMST0dMKyVgGeeeQaqquK//uu/kJGRgfr6evz2t78d77EREdE4UlWBdncwZkGcFmcAp9o8Q9YG8IUUTEujjoFjkWUxoqwgc8xpgUBYRYc7hB5/GAV2c1LrKgxrJeCiiy6C3+9Hc3Mzvve97433mIiIaAK4AzJkVcT8VNvXK2CoUwEBWUF5Uea4jC8VOXJt6PSG4A6EYU/zyojAMFcCfv/732Pp0qVYs2YNAOC9997D2rVrx3VgREQ0vjo8QRh0cVIBNR24uCgzYR8ARY10DMxOs46BY6HTSZg/PQthRUVYUZM9nDEbVhDwb//2bzh06BBycnIAAEuXLkVdXd04DouIiMaTEAItrgAyYp0KcAZwqn3oVEA6dwwcC6tJj/nTs+D0hzQ5NphMw/ovZzAYkJ2dPd5jISKiCeIJygjLaswJfF9fgaC5iYOAdO8YOBaFWWYUZVng9IeTPZQxGXYDoeeeew6KoqCmpgZf+9rXcOWVV4732IiIaJx0e0Nxd7jvO9WOeUV2FCaY4GVFhVGvS/uOgaMlSRLKCiOlfgPh5BT60cKwgoBHH30UR48ehdlsxqZNm5CVlYWf/vSn4zw0IiIaL83OQMxTAc1OP2rbvbiqLD/h+z1BGTNyLWnfMXAsTAYdKkqy4A7KUNM0LTDsOgEPPfQQHnroofEeDxERjTNfSIY/rCA/Y/AUMNxUgKyqk7JM8Ejl2EyYnW/DuS4f8jPS73kkDAKGOgHAVsJEROmnxxs/j73vVMeQqYCQrMJmMkyqjoFjcVF+Brq8IXiDMjLS7JkkHO1bb72FmTNnYtOmTbj88svTfhckEREBza4AMmKkApp6/Djd7sWdV5UmfL8nGMbcwsnVMXAs9DoJC4qz8HZdFyxG/airCSZDwiCgpaUFr776Kp5//nk899xz+PSnP41NmzahoqJiosZHREQaCoQVuPzhmEv553sFJN4PIATScul7PGWYDbi40I7jrW4UpFGaJOHGQL1ejzVr1mD37t04ePAgysrKcN111+HRRx+dqPEREZGGXP4w4n2Aj6YC7PFTAf6QArvVOGk7Bo5FcY4FBXYTXIH0OTY4ZPIiGAzij3/8I55//nnU1dXhvvvuw0033TQRYyMiIo21uAKwGgdP4E09fpzu8OLOqxOnAnxhGQvy7eM1vLQmSRLKC+2orutCSFZhMqR+EaWEQcDmzZtx5MgR3HDDDfje976HRYsWTdS4iIhIY2FFRZc3hDzb4I5/wzkV0LcvbLJ3DBwLi1GP+dPteL/BiYJMc8rvm0gYBDzzzDPIyMjAyZMn8bOf/Sx6XQgBSZLgcrnGfYBERKQNV291u1gT0/5THZg/3Y4Ce/x8ti+kYFqmeUp0DByLaXYLZuWF0eT0I8+W2vsDEgYBqpr+zRGIiCiizR2EWT94Am/sjqQC7hoiFRCQFZRnT52OgWNROi1ybNAXkmMWZUoVqZ+wICKiMVNUgQ5PEDbz4CBg3zDaBiu9LYenUsfAsTDodVhQkgVfSIGipu7xegYBRERTgDsQjrb+vdD+Ux1YMN2esAJgpGOgZcp1DByLLIsRcwsy0O0LJXsocfG/JhHRFNDhCcKgG/xPfkO3D2c6vLi6PHGZ4LCqoijB0UGKzZFrQ26GCe4UPTbIIICIaJJTVYEWZzBmmd9ogaAEpwJkRYVBp4Pdkrq57VSl00mYV2SHrKoIK6m3zy4pQUBPTw/Wr1+P+fPnY8GCBXjrrbfQ1dWF1atXo7y8HKtXr0Z3d3f09du2bUNZWRnmzZuHV155JRlDJiJKW56QDFlVY5az3XeqAwuKsxKmAjwhdgwcC6tJj/nTs9DjC6Vc+f2kBAFf//rXsWbNGhw/fhzvv/8+FixYgO3bt2PVqlWoqanBqlWrsH37dgDAsWPHUFVVhaNHj2Lv3r3YunUrFCV9ezcTEU20Lk8o5l6Ahm4f6jp9uDrBhkAgsimwgKmAMSnMsqA4x4oef2qlBSY8CHC5XPj73/+OO++8EwBgMpmQk5ODPXv2YPPmzQAiRYpefPFFAMCePXuwceNGmM1mlJaWoqysDIcOHZroYRMRpSUhBFqcgYSpgKvmxu8VEJJVWI16ZLBM8JjNLciEQS8hEE6dD7ITHgScPn0aBQUFuP3223HppZfirrvugtfrRWtrK4qLiwEAxcXFaGtrAwA0NjZi5syZ0fc7HA40NjbGvPfOnTtRWVmJyspKtLe3j/8fhogoxflCCgJhBcYYu/r3nerAwuIs5CdIBbiDYczItaZ85bt0YDLoUFGcDXcgjFTJCkx4ECDLMt555x3cc889ePfdd5GRkRFd+o8lVv4k3l/GLVu2oLq6GtXV1SgoKNBszERE6arbF4rZMKi+NxWQqDZAn0T7BWhksm1GlKbQscEJDwIcDgccDgcuv/xyAMD69evxzjvvoKioCM3NzQCA5uZmFBYWRl9fX18ffX9DQwNKSkomethERGmpxRmIWbFu/6kOSEicCvCHFGRbjbDEaDhEozcrLwP59tTovzDhQcD06dMxc+ZMnDhxAgDwl7/8BQsXLsTatWuxe/duAMDu3buxbt06AMDatWtRVVWFYDCIM2fOoKamBitWrJjoYRMRpZ1AWIEnKMecxPfVdGBhSeJUgC8sY0aOdTyHOCXpdRIWFmehIAVWWJJy6PPRRx/FLbfcglAohDlz5uDJJ5+EqqrYsGEDdu3ahVmzZuGFF14AAFRUVGDDhg1YuHAhDAYDduzYAX2M2tdERDSQM86Sc32XD2e7fNhyzZy47412DIzRcZDGzmYywJaX/LoLkki1Q4saqaysRHV1dbKHQUSUNO/XdyMQVgelA54/dA7PHzqHJ2+7LO5KgCcoI8tqQEVJ9kQMlTQwmnmPFQOJiCahkKyi2xeGNVYq4NTQqYCgrKA4m6mAyY5BABHRJOTqrVV/4Wmqc10+nOtKXCCIHQOnDgYBRESTUJsrALNh8CpA36mARL0CPL0dA2OVGabJhUEAEdEkIysqOjwh2GJU+XuzNxWQlxF/w5/MjoFTBoMAIqJJxh2QoahiUL+As51e1Hf5cE2CVEBYUWHUs2PgVMEggIhokmn3BGGKUSZ4OKkAbyhSG4AdA6cGBgFERJOIqgq0uQLIiNEwaN+pDlSUZCE3USpAEZhmT34RG5oYDAKIiCYRd1CGrKqDNvWd7fSivtuPq8vj91UJygoyzfqYHQdpcmIQQEQ0iXR6gtBLsTsGSgCunBO/V4AnKGNGjm0cR0ephkEAEdEkIYRAa4xUgBAC+091YNGM7ISpAADIy2SZ4KmEQQAR0SThDSkIhBUYL9gUeK7LF0kFJDgV4AvJ7Bg4BTEIICKaJLq9Ieh1g/9Zf/NUB3QSsDJB22BfSGHHwCmIQQAR0STR3ONHhilOKqAkG7lxOgKqQkAnsWPgVMQggIhoEvCHFHhDMkyGgf+sn+30oaHbj6vLE6UCFBTYzYPeS5Mf/4sTEU0CPb7QoGZBALCvtjcVkOBUQCCsYDo7Bk5JDAKIiCaBZlcgZipgX03kVEC8pX5FFTDq2TFwqmIQQESU5oKyAqcvPGhnf12nD409iU8FeIIypmezY+BUxSCAiCjNOX1hxMgEYP+poVMBsqqigB0DpywGAUREaa7VHYT1glUAIQT2nerA4gSpgLCiwmTQIYsdA6csBgFERGksrKjo9AwOAuo6vWjs8eOqBKkAb1DGjGxrzA2FNDUwCCAiSmPugAwAgybyfac6oZMStw1WVHYMnOoYBBARpbE2VwBmfYxUQE07Fs/IjrvrPygryDDrY7YcpqmDQQARUZpSVIF2TxBW0+BUQJMzgKvL4rcNZsdAAhgEEBGlLU9AhqKKQcf73qxJ3CtACAEh2DGQGAQQEaWtdk8AhgsaBvWdCrjEkRM3FeAPK8jLYMdAYhBARJSWhBBocQaRcUEq4EyHF83OQMICQf6wgmKWCSYwCCAiSkueoAxZVWHQD/xnfF9vgaAr4hQIUoUAAORmMBVADAKIiNJStzc0aC/AcFIBvqCCQrsFRj3/+ScGAUREaanJObhh0OlhpAKCioLibJYJpggGAUREacYXkhEIKYM+ze+rSZwKUFQBg05CFjsGUi8GAUREaabHO7hhkBAC+2s7sCRBKoAdA+lCDAKIiNJMk9MP2wWpgNr23lRAefxUADsG0oUYBBARpZFAWIEnKA864x89FVAaOxUQVlSY2TGQLsAggIgojbj84UHXhBDYf6oDS2fmxM33e4MyZuSwYyANxCCAiCiNNDsDg9oG17Z70eIKJGwbLKsC+ZnsGEgDMQggIkoTIVlFty80KAjYd6odep0UNxUQCCuwmw3sGEiDMAggIkoT7kAkFdB/Sb+vQNASR/xUgC8kY0YuywTTYEkLAhRFwaWXXoobb7wRANDV1YXVq1ejvLwcq1evRnd3d/S127ZtQ1lZGebNm4dXXnklWUMmIkqqVlcAZsPAVYBTbR60uoK4uix+x0BVsEwwxZa0IOCRRx7BggULol9v374dq1atQk1NDVatWoXt27cDAI4dO4aqqiocPXoUe/fuxdatW6EoSrKGTUQ04cKKii5vCO2eIGwXNAzaX9sRSQXEKRDEjoGUSFKCgIaGBvzxj3/EXXfdFb22Z88ebN68GQCwefNmvPjii9HrGzduhNlsRmlpKcrKynDo0KFkDJuIaML4QwpanX68X9+D/ac68GFDDywGPXQXpALerImkAuyW+G2DS3JsEzVsSjNJ2SXyjW98Az/60Y/gdruj11pbW1FcXAwAKC4uRltbGwCgsbERV1xxRfR1DocDjY2NEztgIqJxpqoCnpCMHm8ILc4AfCEFkgRYjHrk2Uwxj/bVtHnQ5g5i02WzYt9TCEgSkGNjmWCKbcKDgD/84Q8oLCzE8uXL8cYbbwz5etHb9rK/eOdcd+7ciZ07dwIA2tvbxzROIqLxFlZUuAMyOjxBtLkCkFUBvU6CzWhAfubQ/zzvP9UBQ4JUADsG0lAmPAjYv38/XnrpJbz88ssIBAJwuVz40pe+hKKiIjQ3N6O4uBjNzc0oLCwEEPnkX19fH31/Q0MDSkpKYt57y5Yt2LJlCwCgsrJy/P8wREQjFAgrcPpCaHUH0eUNAQBMeh0yzcYR1fSPngqYmYPMOFUAA4qC6Vl2TcZNk9OEh4fbtm1DQ0MD6urqUFVVhU984hP41a9+hbVr12L37t0AgN27d2PdunUAgLVr16KqqgrBYBBnzpxBTU0NVqxYMdHDJiIaFSEE3IEwznV68faZThys7cRHLW74QwrybCbkZ5hht4wsAADOpwLitQ1WVAEjOwbSEFKmcsQDDzyADRs2YNeuXZg1axZeeOEFAEBFRQU2bNiAhQsXwmAwYMeOHdDrucuViFKX3LvM3+4Jos0VhKyq0EsSbKbhLfMPx76+VECcAkGeoIzibCs7BlJCkoiVdJ8EKisrUV1dnexhENEUEQgrcPnDaHEF0O0LQQjAqNMhw2zQfCIWQuCup6sxK8+G732mIuZrOr1BLLsoF1lxTg3Q5DOaeS9lVgKIiNKJEAKeoIweXxgtTj+8oUj9EqtRj1xr7N38WulLBXxxRexTAX0dA+0sE0xD4N8QIqJh6lvm7/QG0eIMQlZU6HUSrCY98jMmrjnPmzWRVMDlcU4FeIIyZufb2DGQhsQggIgogb5l/shu/mB0mT9zHJb5h0MIgf21kbbBmXE+6atCYJqdHQNpaAwCiIj66b/M3+oKwBOUAUzMMv9wnGz1oN0dxJcuj50KCIQVZJoNsJn4zzsNjX9LiGjKS5Vl/uHYd6odBp2EFXFOBXhDMuYVsTYADQ+DACKakvqW+dvcQXT2W+bPMOlh0KfmjvpIgaBOXDordiqg77AXOwbScDEIIKIpQQgBb0hBtzc0aJk/x2oa0JgnVZ1odaPDE8SXr4idCoh0DDSxYyANG4MAIpq0ZEWFJxipzd/qCiIsq9DpJNhScJl/OPb1ngqIlwrwhxXMLcyc4FFROmMQQESTSiCswB2Q0eoKoMsbgioEDH3L/GlcOEcVAvtrO7FsVm7MVEBfx8BslgmmEWAQQERpL7KbP9KC1x2QIy14DXpkW41pscw/FItRB0UV2PHFS2E16WEx6hAIqwNe4w3KKMpix0AaGQYBRJS2VFWgrtOLs51e6HU6WI16TMtMv2X+RCxGHQKyint+dRgN3X44cq147EvLBwUCIUXF9CxLEkdK6YghIxGlpaCs4MNGJ851+ZCXYUaubXJtiPMEZBw83YluXzgaAABAQ7cf9/zqMEyG8/98K6qAQSexTwCNGFcCiCjtOP1hHGl0AkBabvCLxROUcazJiQ8anPiwyYkz7V4IAP+z5YpoANCnodsPtV/rN3cgjJIcK3TsGEgjxCCAiNKGEAJNPQGcbHUj02xI60/+vpCMo00ufNDgxJFGJ053eKAKwKiXsGB6FjatmIVLHNkoyrLAkWsdEAg4cq3oP98rqkBB1uQIhmhiMQggorQgKypq2jxocQaQazMlpW7/WPhCMo41u3CkMfJpv7Y9MukbdBLmTbfjC5UzsXhGNuZNzxqw1B9WVDz2peWD9gSEZDX6fYtRz46BNCr8W0NEKc8XknGs0QVvSEZ+RvLr9w9HIKzgWLMLHzY48WGjEzVt7uikf3GRHTf3Tvrzp9thNsRf0QiEVViMOjx71+VQBaCTgJCsRjcFeoIySqdlpMUzodTDIICIUlqHO4CjTS6YDXrkpXD+PxBWcLzFjQ8aenCk0YmTbR4oqoBeJ+HiwkysX35+0h9pGiMQVgcdCeyjqCryM1kmmEaHQQARpaS+439nOr3ItZpS7vx7UI5M+n2f9E+2uiGrAjoJKC+043NLZ2CxIxsLpmfBahqfvQuBsAK71ciOgTRq/JtDRCknKCs43uxGty+EaRnmlCj4E5JVnGhx4cNGJz5odOJEy/lJv6wwE+uWlmDxjBwsKLZP2KTsDclYMJ0dA2n0GAQQUUpJleN/YUXFiRY3PmyMfNI/3uJCWIlM+nOmZeIzS0pwyYxsLCzJSsoncSEEBIAcdgykMWAQQEQpIdnH/8K9pw8+bOjBh41OfNTiRkhWIQEoLcjApxcXY/GMbCwsyY5Zu3+i+UIK8jNMCTcVEg0l+X+TiWjKS8bxP1lRcarNgw96P+l/1OxCsPfYXem0DKypmI7FM7KxqCQbmZbU+6cyICsoK2LHQBqb1PubTURTykQd/1NUgVNtnt7l/R4ca3ZFd9zPzrfh+oVFWDwjGxUl2chK8U58qhDQSRJyUnyclPoYBBBR0rS7Azg2huN/FqMOJoMuOin2Pz+vqAK17Z5IcZ5GJ441ueAPKwCAWXk2rJofmfQXzchOu/a7kY6BZhhS7MQEpR8GAUQ04bQ4/tfXXe/2p96OVtLb8cVleLuuE3uPtOJYswu+UGTSn5lrxXXzCnCJIweLSrKQY0vvzXRBWUUROwaSBhgEENGECoQVnGgZ+/E/k0EXDQCASFOde597B//3xoVodgZwbXkBLnFEcvq5k2gHvayoMOp17BhImmAQQEQTRsvjf0FZjdldr6wgE49/afmY7p3KPEEZjjx2DCRtMKFERONOCIHGbj/eOdsN0xg/xfb4QnjkLydxpsMLR651wPccuVYY9JN7clRUgWmZqVs+mdILgwAiGldhRcXxFjdOtrqRazON+vy/ogr88YMm3P3sYbx+oh0fNTnx2C3LooHAhd31JqOQrMJi0qdEnQKaHPg3iYjGjS8k42ijM1rYZrTH/z5qduHxv9XidIcXS2fmYMu1czAz15awu95k5AmGMacgkx0DSTMMAohoXIz1+B8AdPtCeOpAHf56vA3TMk14YM18XDk3PzoJJuquN9kIISAEmAogTTEIICJN9R3/q+v0ImeUx/8UVeCPHzbjuX+cRVBWsX6ZA1+4bOaElxJOppCsIigrCCnng5zpOZZx60hIUxODACLSTP/jf/mjPP53tMmJx/9Wi7pOHy7tXfp35NrGYbSpI6yoCIZVBBUlei3TbECB3YwcmwlWkx5Wo35CyinT1MIggIg0Mdbjf13eEJ48cAZvnGhHgd2Mf71hPlbOyZ90+e9YE77NqEe+3YQcmxFWY2TCZzVAmggMAohoTMba/U9WVPzhw2Y8949zCCsqNlTOxM3LHZNi6V9WVARkFSE5MuELABajHnmZRuRYbbCaDbAa9aNKmRBpgUEAEY1auLcT32i7/33Y6MR//60WZ7t8WH5RLrZcMwclOdah35iCFFUgEI7k8FUhAABmgw45NhNybTZYTQbYTJzwKbUwCCCiUfEGZRxrGt3xv05PEE8eqMPfTraj0G7Gg59agMtL89Jm6V9RBYKygmBYhYrIhG/U65CbYUSOxYQMswFWkx4mAyd8Sm0MAohoxEZ7/E9WVPz+gyY8f6gesqpi42Uz8fllqb30H53w5cgnfAmAXqdDrs0IR64xOuGbDan7ZyCKZ8KDgPr6etx6661oaWmBTqfDli1b8PWvfx1dXV34whe+gLq6OsyePRu//vWvkZubCwDYtm0bdu3aBb1ej5/97Gf4p3/6p4keNhFhbMf/PmjoweN/P436Lh8um52Lr1wzB8XZqbX0rwoR2bQnK71L+hJ0OiDHZkJJjhF2S2TPQyoHLUQjIQnRm7yaIM3NzWhubsayZcvgdruxfPlyvPjii3jqqaeQl5eHBx54ANu3b0d3dzd++MMf4tixY9i0aRMOHTqEpqYmfPKTn8TJkyeh1yf+n7CyshLV1dUT9Kcimvz6H//LtZmGffyv0xPErv1n8GZNB4qyzNhyzRysKM0f59EOTRWi9yy+CllVIQHQ6SRkWYzItRlhtxh7P+Hr0iZNQVPbaOa9CV8JKC4uRnFxMQDAbrdjwYIFaGxsxJ49e/DGG28AADZv3ozrrrsOP/zhD7Fnzx5s3LgRZrMZpaWlKCsrw6FDh7By5cqJHjrRlOX0h/FhYw8kSMM+/hdWVLz0fhOq3j4HVQW+uGIWblo2IynL5kIIBC+c8CUJdqsBBXYrsqyRo3kWIyd8mlqSuiegrq4O7777Li6//HK0trZGg4Pi4mK0tbUBABobG3HFFVdE3+NwONDY2Bjzfjt37sTOnTsBAO3t7eM8eqLJTwiBxh4/Tra6YTcbh70M/l59D/7777Vo6Pbj8tI83HXNHEzPsozzaCOEEAgpkXLCfRO+JEnItBgww25BVu8nfItBz3a8NOUlLQjweDz4/Oc/j5/+9KfIysqK+7pY2Yp4kfqWLVuwZcsWAJFlESIavf7H//Js5mEd/2t3B7Fr32nsr+1EcbYF37txISpn543rOBVVwBeSI+V1Re+EbzagONuC7H7FdzjhEw2WlCAgHA7j85//PG655RbcdNNNAICioiI0NzejuLgYzc3NKCwsBBD55F9fXx99b0NDA0pKSpIxbKIpY6TH/8KKihffbcT/VNdDAPjS5bPwuUsd43pELigr8ARl6CQJ07MtyGV5XaIRm/BDrEII3HnnnViwYAHuv//+6PW1a9di9+7dAIDdu3dj3bp10etVVVUIBoM4c+YMampqsGLFiokeNtGU0e4OoLquC4oK5GWYhwwA3jnXja89/y6ePngWl87Kwc+/uAxfuGzWuAQAQgh4AjI6PEEoqsC8IjtWzs3HxUV2FNjNyDQbGAAQjcCErwTs378fzzzzDBYvXoylS5cCAH7wgx/ggQcewIYNG7Br1y7MmjULL7zwAgCgoqICGzZswMKFC2EwGLBjx44hTwYQ0ciN9PhfmzuAX7x5Bm+d7kRJtgX/9pkKLL8od1zGJisq3EEZiipQmGXGwpwsZFkN3MRHNEYTfkRwovCIINHwBcIKjre40e0NIS8j8fG/sKLid+824tfV9ZAAfKFyJj576YxxKYfrDynwhsMw6nSYlWdDYZaFZ/SJ4kiLI4JElFqcvjA+bIoc/5uWmfj4X/XZLuz8+2k0OwO4cm4+7ry6FIV2bXf9q0LAHZARVlRk24xYXJgzqr4ERDQ0BgFEU9RIjv+1ugJ44s3T+MeZLszIseL7aytw6Sxtl/5Dsgp3MAydBBRnWzE92wK7xajpzyCigRgEEE1Bwz3+F5JV/PadBvzmcAN0OmDzytlYt7REs6V/IQS8IQVBWYHFoMfFhXZMs5vZeIdogjAIIJpivEEZR5uc8A9x/O/QmS488eZptLgCuLpsGu68unTIdMFwKaqAOxCGIgTyM0xYUGxHttXIjX5EE4xBANEUMpzufy3OAHa+WYu367oxM9eKf1+3CEtm5mjy8wNhBd6QDINOwsw8G4qyLLCauNGPKFkYBBBNAX3H/850eJFri338Lygr+O3hBvzmnQbodRJuv3I2PrNk7Ev/5zf6KbBbjVhYnIW8DBMM43CagIhGhkEA0STX//jftEzzoON/Qggcqoss/be6gri2vAB3XDUb+WNc+g8rKtyBMABgerYFxTlWZHGjH1FKYRBANIkNdfyvqcePJ948jeqz3ZiVZ8MPPrsIix05o/55Qgj4QgoCsgKzQYeywkxMs5uT0jmQiIbGIIBokgmEFfT4QmjsCcDpC8FuGXz8LxBW8JvDDfjtOw0w6nW486pS3HhJ8aiX6BVVwB0MQ1EEcjNNmDc9stGPTXuIUhuDAKJJICgrcPrCaOzxo8cXhiQBGSYDCi4o5COEwMHTnfjFvjNocwdx3cUFuP2qUuRlmEb1c/s2+ul1Ehy5VhRlWWAz8Z8VonTB/1uJ0lRIVtHjC6HZ6UeXLwwJgNWoj3vsr7Hbj51v1uKdcz2YnW/Dts8txqIZ2SP+uWpvE5+QosJuNnCjH1EaYxBAlEbCigqnP4ymHj+6vCEIAdhMeuTbBk78FqMOJoMOam9rkN+/14Sf/fUUTAYdvnJNKT69uGTEZXjDigpXIAwBoDjbguJsK7IsbOJDlM4YBBCluLCiwuUPo9kZQIcnCACwGPTIs8X+xG8x6hCQVdz+1Nto6PbDkWvFDz9/CTavvAgfu7gQuSNc+veFZPjDCkwGHeYWZKLAbmYTH6JJgkEAUQqSFRWugIwWZwDtngCEQKTAT5yJv48QAkIA9/zqMBq6/QCAhm4//uW3H+DZuy6Hyy8P6+crqoAnGIasCuRYjSgvsiOHG/2IJh0GAUQpQlEFXP4wWt0BtLmCUIWAWa9HjjVxa99AWMEHDT14u64b1We7seOLl0YDgD4N3X6ow2gaHpQVuAORjX4zciJNfDLM/GeCaLLi/91ESdRXQ7/VFUCbOwhFFTDpdZHjdQkm/hZnANVnu/B2XTc+bOxBWBGwGHVYOjMHFqMejlzrgEDAkWtFvA/xQgh4gjKCsooMsx4Li+3IyzRr1iSIiFIXgwCiCaaqkTK6be4AWl0ByKqAUaeD3WyMu1lPVlQca3b1ftrvik7wJdkW3LCoGJUX5WLRjGwY9TpYjDo89qXl0ZSAI9eKx760HCFZHXRPVyAMVQBFWRbMyLEiy8qNfkRTCYMAogmgqgLuoIwOdxDNTj9kVcCg0yEzwcTf7Q3h8NluvH22C+/V98AXUmDQSVg0IxtrKqbjstl5KMmxDnpfIKzCYtTh2bsuhyoAnRQ5ThgIR4KAvo1+Bp0OpdMyUJhl4UY/oimKQQDROBEiMvF3eoJo7A5AVtWEE78qBE61eVBd14W3z3bjVJsHAJCXYcLVZdNQOTsPSxzZwyrGEwifn/SBvo1+MsKKimybEWWFmcixmUZ8TJCIJhcGAUQa6suvd3pCaHL6EZJV6CUJmWYDDPrBzXO8QRnv1vfg7bouvHO2Gz3+SNGfedPt+NLls1A5Ow9zpmWMeok+JKtwBUPQSRJKsq0ozrEikxv9iKgX/zUgGiMhBLwhBV2eIBp7/AiGVeh0kYnfbjYOem19tx/VdV2oPtuNY80uKKpAhlmP5bNyUTk7D8tm5SLbOrpue2FFhS+kQFYjqwBWox7zi7KQn2mGycCNfkQ0EIMAolHyBmV0eUNo6vHDH1aglyRkmA3IvGDiD8oKPmx0orquG2/XdaHNHSn4Mzvfhs8tnYHK2bmYPz1rVEvzIVmFP6xAUVWovdUDp2ebkWszIcNsYK6fiBJiEEA0Ar5QZOJv7vHDG1Sg00nIMBmQnzHwf6U2dyA66X/Q6ERIVmE26LDEkYP1yx1YflEuCi9o7jMcQVmBP6RAFQIqgEyTHiU5FuTYTMgw69myl4hGhEEA0RD8IQXd3kiO3x0IQ9f7iT8/0xx9jaIKfNTsQvXZblTXdeFslw8AUJRlxvULilA5Ow+LZ2SPaEleCIFg7yf9vh4AdrMRs/JsyLYZYTMZuMRPRGPCIIAohkA4MvE3Ov1w++Voa95pmec/vTv9YRw+Gzm3/865bniDCvQ6CRXFWbjjqtmonJ0HR4512Jv6hBAIhFUEZQWKEJAkwG4xojQ7A1kWI2xmPQv4EJGmGARMUX015kXf7wH0ftiEJAE6SYIETKla8YGwAqcvhEZnAK7eXfo2kwHTej/xi74jfGe7UF3XjZOtbggAOTYjVs7JR+VFeVg6M2fYZXZVIRAIKwjKKlQhIEFCts2A4pwMZFmNyDDp2Z6XiMYVg4ARiDdxCojoBNr/60ETbILvicg3IRApLKMKAVVEXqMK9H59/r2qKqBCQKiA2vu1QO971Mj9Iq85fw9xwT36yzQbkGMzoW/K7/aF4A0qEBDQ9QYFep004Fedrvd67+/1Ol30tQadBKn3V51OgiQBEiToJECSzv86IODo/Trea8dDUFbg9IXR7PSjyxuGJAE2owH5GZGJ3xeS8VZtB94+243Ddd3o8oUAAOWFmdi0YhYqL8rF3MLMhCV++yiqiOT0wwrQ++fNsRkxM9eGTKsBGSYDz+0T0YRiEDAEVRXo9IbgDoQhBNDjC8ETkiOzqCQhOnMDAHq/7tV/or3wn/a4vVxE7+TY+yYJ5ydACecnyPPf67u/dP73va8x6Pu/T+r3PQzqPR+QVdz25KEBZWazrUYEwmo0iAAQzU2rQkCR+wKZ3sCiN5gZ8NreAObCByBFH5UEAREZjxDRZ3jh8+lblYgEHX1BiQ6SDr1BiAS9DtAhEqRIEmDoDUr0uvOBS1/goQqBFmcAnd4QJESO0uVnRDr0NXb78faxLlTXdeFokwuyKmAz6XHprFxcdlEull2Ui1zb0O14FTXyST8gK9EgJzfDhIvybcgwRyb9qbTSQkSph0FAAqoqcKLVja88XR2dHH9+yzJk64zwBJXeCVDq96m+bzIcvEogBKD2+7R/flUh8musSRYC0U/ykWuRm6pAdHKN9XNijyPeSobAx+cVYutz7wxoPXvPrw7jv7+0HO+e64ZBr4NRL8Ggi/xq1Otg6P3VqOv9vUE3rp9iB/xZ+v1eUQQUiH7P7cJn0z/VIZDRb8XDbjbAYtDBE4wc4Tt8NrKbv9kZAADMzLNh7ZISVF6UiwXFWUMuzSuqgD+sINg36esk5GeYUJphQ4bFCJtRz0mfiFIKg4AEOr2haAAARCbHrc++g/9740L88zOHkzw67VwxJz9m61lPUMb/99LRYd9H1/vp26iXBgUOiQIJg67vax2Mun6vveDraMBxwfej74318/R6GPWRlQCrST9oxeMnG5bgJ6+cwKG6bpj0Oix2ZGPdkhIsn52H6VmJj/DJSmTnfkiJFOYx6nXIzzAhLzMDmWYDrEY9m/EQUUpjEJBASFZiTo6OXCv++do50WX1vuX38/ns80v657/f/1rs3+t6cwA6CYAUWdruW/bX9d5Eh8ib+j5Q6uL9nH7pgvNjk2JcA/IzTDFbz2bbjPjh5y9BWFEhKyLyqxr5Nda1eF+HVQFZURFWIr+GFAFvKHz+mnr+e+F+99CSBOC/v7wc3//DsQFB3bd+/T5+cvMSHG1yYfGM7ITFdcKKCn9IQVhVAQGYjTpMs0cK82SaDbAYdZz0iSitMAhIwGSI3Zc902zAjZeUJHFk2hJAzNazeknCwuKs5IxJCMiqGEagMTCQiL6mf2DRG4QUZ1tjBnWFWZaYTXkuLMFrMepRmGVGXgar8RHR5MAgIIH8DBOeuLVy0J6AdncwunM+ruhxu/OfDKULftP/M2OsT5BS9Hv9rw1+c6zXXfj6mN/rvRaUFdhMejxz5+WAiGzSC4RleINK9HUXbiYcb5Ik9aYNACu0mWyzrIaYQZ2u33MIhNXopG8z6lGcbUFuhgk2k56TPhFNOgwCEtDpJMwrsuN/t16FQDjSy91i1EWO2/W+RvSLA6JBwYBr5zfkAec3/vXf5Dbwtf2ODGLw8cPBr+3dkNjvG9HNf9Ex9rtH3337/gy9r/eHVPhCkZr2fUcLVSGg9jtuGC/k6QsQLvyz9P/+gK+k8xsk+78h1qmIkZ6IuDBg6f91SFYHrXj0BXWeoIxMswEzci3IsZpgYwleIpoCGAQMQaeTUGA3D/3CKaJ/zYHh1kno27E/VI2EvqBDEWLI+giqiuh1RYhobQRFiV0boW+LQabZgKduXxHdvW8xSLCaDCzBS0RTEoMAGhGp95z+4MoHqa//EUkALMxDRFMegwCaMvpOTRARUUTarH/u3bsX8+bNQ1lZGbZv357s4RAREaW9tAgCFEXBvffeiz/96U84duwYnn/+eRw7dizZwyIiIkpraREEHDp0CGVlZZgzZw5MJhM2btyIPXv2JHtYREREaS0tgoDGxkbMnDkz+rXD4UBjY2MSR0RERJT+0mJjYP/z531iFa7ZuXMndu7cCQA4fvw4KisrAQBOpxPZ2dkx7x3re8O5duHX7e3tKCgoGOafaHQS/Tm0eu9Qr+Oz5LMczli1ft9on+VIrvNZJv4+n+XIXzvRz7Kurm5Y4x5ApIEDBw6I66+/Pvr1D37wA/GDH/xg2O//yle+MqLvDefahV8vX7582OMZrUR/Dq3eO9Tr+Cz5LIcaw3i8b7TPciTX+SwTf5/PcuSvTYdnmRbpgMsuuww1NTU4c+YMQqEQqqqqsHbt2mG//zOf+cyIvjeca4nuOV7G8jOH+96hXsdnyWc51BjG432jfZYjuc5nmfj7fJYjf206PEtJiBhr7Sno5Zdfxje+8Q0oioI77rgDDz74YLKHNEBlZSWqq6uTPYxJgc9SO3yW2uGz1A6fpXbG+izTYk8AAHzqU5/Cpz71qWQPI64tW7YkewiTBp+ldvgstcNnqR0+S+2M9VmmzUoAERERaSst9gQQERGR9hgEEBERTVEMAoiIiKYoBgHj5PTp07jzzjuxfv36ZA8l7b344ov4yle+gnXr1uHPf/5zsoeT1j766CPcfffdWL9+PR577LFkDyeteb1eLF++HH/4wx+SPZS09sYbb+Caa67B3XffjTfeeCPZw0lrqqriwQcfxNe+9jXs3r17WO9hEDACd9xxBwoLC7Fo0aIB12N1OJwzZw527dqVjGGmhZE8y89+9rN44okn8NRTT+F//ud/kjHclDaSZ7lgwQI8/vjj+PWvf80jWhcYyXMEgB/+8IfYsGHDRA8zLYzkWUqShMzMTAQCATgcjmQMN6WN5Fnu2bMHjY2NMBqNw3+WYyo1NMX87W9/E4cPHxYVFRXRa7Isizlz5oja2loRDAbFJZdcIo4ePRr9/uc///lkDDXljeZZ3n///eLw4cPJGG5KG+mz3LNnj1i5cqV49tlnkzXklDSS5/jqq6+K559/Xjz55JPi97//fRJHnZpG8iwVRRFCCNHS0iK++MUvJmvIKWskz3Lbtm3i8ccfF0IMf+7hSsAIXHvttcjLyxtwjR0OR2ckz1IIgX/5l3/BDTfcgGXLliVpxKlrpH8v165diwMHDuDZZ59NxnBT1kie4+uvv46DBw/iueeewxNPPAFVVZM06tQ0kmep00WmodzcXASDwWQMN6WN5Fk6HA7k5uYCAPR6/bDunzbFglJVrA6H//jHP9DZ2YkHH3wQ7777LrZt24Z//dd/TeIo00O8Z/noo4/itddeg9PpxKlTp3D33XcncZTpId6zfOONN/C73/0OwWAwpYtvpYp4z/G//uu/AABPPfUUpk2bFp3IKL54z/J3v/sdXnnlFfT09OCrX/1qEkeYPuI9y69//ev42te+hjfffBPXXnvtsO7FIGCMRJwOh/n5+Xj88ceTMKL0Fe9Z3nfffbjvvvuSMKL0Fe9ZXnfddbjuuusmfkBpKt5z7HPbbbdN4GjSW7xnedNNN+Gmm25KwojSV7xnabPZRrwXjeHrGDkcDtTX10e/bmhoQElJSRJHlL74LLXDZ6kNPkft8FlqR8tnySBgjMba4ZDO47PUDp+lNvgctcNnqR1Nn6W2+xgnt40bN4rp06cLg8EgZsyYIX7xi18IIYT44x//KMrLy8WcOXPEv//7vyd5lOmBz1I7fJba4HPUDp+ldsb7WbKBEBER0RTFdAAREdEUxSCAiIhoimIQQERENEUxCCAiIpqiGAQQERFNUQwCiIiIpigGAUSkubq6Ojz33HMjfl9PTw9+/vOfR79uamrC+vXrtRwaEfXDIICINJcoCJBlOe77LgwCSkpK8Jvf/Ebz8RFRBIMAoknk6aefxiWXXIIlS5bgy1/+Ms6ePYtVq1bhkksuwapVq3Du3DkAkcY39913H6688krMmTNnwET7ox/9CIsXL8aSJUvwwAMPAABqa2uxZs0aLF++HNdccw2OHz+e8D4PPPAA3nzzTSxduhQPP/wwnnrqKdx88834zGc+g+uvvx4ejwerVq3CsmXLsHjx4mib4wceeAC1tbVYunQpvvOd76Curg6LFi0CAAQCAdx+++1YvHgxLr30Urz++usAIp38brrpJqxZswbl5eX47ne/OzEPm2gy0KawIREl25EjR8TFF18s2tvbhRBCdHZ2ihtvvFE89dRTQgghdu3aJdatWyeEEGLz5s1i/fr1QlEUcfToUTF37lwhhBAvv/yyWLlypfB6vdF7CCHEJz7xCXHy5EkhhBAHDx4UH//4xxPe5/XXXxef/vSno2N78sknxYwZM6L3C4fDwul0CiGEaG9vF3PnzhWqqoozZ86IioqK6Pv6f/3jH/9Y3HbbbUIIIT766CMxc+ZM4ff7xZNPPilKS0tFT0+P8Pv9YtasWeLcuXNaPlqiSYuthIkmib/+9a9Yv349pk2bBgDIy8vDW2+9hd/97ncAgC9/+csDPiV/9rOfhU6nw8KFC9Ha2goAeO2113D77bfDZrNF7+HxeHDgwAHcfPPN0fcGg8GE94ll9erVyMvLAxBphfp//s//wd///nfodDo0NjYmfC8A7Nu3D1/72tcAAPPnz8dFF12EkydPAgBWrVqF7OxsAMDChQtx9uzZAf3WiSg2BgFEk4QQYkCv+1j6f99sNg94b7x7qKqKnJwcvPfeezHvGes+sWRkZER//+yzz6K9vR2HDx+G0WjE7NmzEQgEEo490b37j0Gv1yfcd0BE53FPANEksWrVKvz6179GZ2cnAKCrqwtXXnklqqqqAEQm3quvvjrhPa6//nr88pe/hM/ni94jKysLpaWleOGFFwBEJuP3338/4X3sdjvcbnfc7zudThQWFsJoNOL111/H2bNnh3zftddei2effRYAcPLkSZw7dw7z5s1LOA4iSoxBANEkUVFRgQcffBAf+9jHsGTJEtx///342c9+hieffBKXXHIJnnnmGTzyyCMJ77FmzRqsXbsWlZWVWLp0KX784x8DiAQQu3btwpIlS1BRURHdyBfPJZdcAoPBgCVLluDhhx8e9P1bbrkF1dXVqKysxLPPPov58+cDAPLz83HVVVdh0aJF+M53vjPgPVu3boWiKFi8eDG+8IUv4KmnnhqwAkBEI8dWwkRERFMUVwKIiIimKAYBREREUxSDACIioimKQQAREdEUxSCAiIhoimIQQERENEUxCCAiIpqiGAQQERFNUf8/munTBwDVcvgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(facecolor = 'w',figsize=(8,6))\n", "plt.rcParams['pdf.fonttype'] = 42\n", "sns.lineplot(data=df[(df['time']==180)&(df['concentration']!=500000)], \n", " x='concentration', y='Mean',ci='sd', style='time',markers=True)\n", "plt.xscale('log')\n", "plt.ylim(0,1400)\n", "plt.savefig('dose_titration_scale.pdf')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.7.9" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": false, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": false }, "varInspector": { "cols": { "lenName": 16, "lenType": 16, "lenVar": 40 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": true } }, "nbformat": 4, "nbformat_minor": 4 }