{
"cells": [
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"##this script will take a dataframe and generate a KDE plot with medians\n",
"## Venice Servellita 2022-01-09"
]
},
{
"cell_type": "code",
"execution_count": 144,
"metadata": {},
"outputs": [],
"source": [
"import pandas as pd\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import seaborn as sns\n",
"import matplotlib.font_manager\n",
"%matplotlib inline\n",
"import seaborn as sns\n",
"from datetime import datetime\n",
"import matplotlib.dates as mdates\n",
"plt.rcParams['svg.fonttype'] = 'none'\n",
"import warnings\n",
"warnings.filterwarnings('ignore')\n",
"sns.set()"
]
},
{
"cell_type": "code",
"execution_count": 145,
"metadata": {},
"outputs": [],
"source": [
"data1 = pd.read_csv(\"data1.csv\")\n",
"df_1 = pd.DataFrame(data=data1)"
]
},
{
"cell_type": "code",
"execution_count": 146,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" sample_id | \n",
" days | \n",
" variant | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" 162 | \n",
" 27 | \n",
" delta | \n",
"
\n",
" \n",
" 1 | \n",
" 183 | \n",
" 18 | \n",
" delta | \n",
"
\n",
" \n",
" 2 | \n",
" B1 | \n",
" 23 | \n",
" delta | \n",
"
\n",
" \n",
" 3 | \n",
" B11 | \n",
" 10 | \n",
" delta | \n",
"
\n",
" \n",
" 4 | \n",
" B13 | \n",
" 13 | \n",
" delta | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" sample_id days variant\n",
"0 162 27 delta\n",
"1 183 18 delta\n",
"2 B1 23 delta\n",
"3 B11 10 delta\n",
"4 B13 13 delta"
]
},
"execution_count": 146,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df_1.head()"
]
},
{
"cell_type": "code",
"execution_count": 147,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" variant | \n",
" days | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" delta_Med | \n",
" 15 | \n",
"
\n",
" \n",
" 1 | \n",
" omicron_Med | \n",
" 11 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" variant days\n",
"0 delta_Med 15\n",
"1 omicron_Med 11"
]
},
"execution_count": 147,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"medians = df_1.groupby(['variant'])['days'].median()\n",
"df_med = medians.add_suffix('_Med').reset_index()\n",
"df_med"
]
},
{
"cell_type": "code",
"execution_count": 148,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"15"
]
},
"execution_count": 148,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"delta_med = df_med.iloc[0][1]\n",
"delta_med"
]
},
{
"cell_type": "code",
"execution_count": 149,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"11"
]
},
"execution_count": 149,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"om_med = df_med.iloc[1][1]\n",
"om_med"
]
},
{
"cell_type": "code",
"execution_count": 150,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAbIAAAFcCAYAAACtPTvuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABhFUlEQVR4nO3deVyU5f7/8dfs7CDKoriLgru5a4llpaXiQtiiadvRTp2yPB1bzcr0eCpLKy3Lb786laa2CFkqmmbmrliJG+auiKyCbMNs9+8PjpOoICrDzQyf5+Ph48HMvb3vAflw3fd1X5dGURQFIYQQwk1p1Q4ghBBCXA8pZEIIIdyaFDIhhBBuTQqZEEIItyaFTAghhFuTQiaEEMKteVwhS05OVjtCrbZ37161I1yTZTcfYdnNR65p20vO+eaby/55OHf9XgtxtTyukInKmc1mtSPUuLp4zlB3z1vUPVLIhBBCuDUpZEIIIdyaFDIhhBBuTQqZEEIItyaFTAghhFuTQiaEEMKtSSETQgjh1qSQCSGEcGtSyIQQQrg1KWRCCCHcmksL2fLlyxk8eDADBw5k4cKFlyzfv38/cXFxDBo0iJdeegmbzQbAzp07iYuLIzY2lr///e/k5+cDsH37dnr16sXw4cMZPnw4L7zwgivjCyGEcAMuK2QZGRnMnj2bRYsWkZCQwJIlSzh06FC5dSZPnszUqVNJSkpCURSWLl0KwAsvvMCbb77J8uXLiYyM5JNPPgFgz549PPzwwyQmJpKYmMjMmTNdFV8IIYSbcFkh27x5M7179yYoKAgfHx8GDRrEqlWrnMvT0tIwm8106dIFgLi4OOfyFStWEBkZidVqJSMjg4CAAABSUlLYuHGjs6WWnp7uqvjCA5h/P0DWC3PwfWk+mU/MoGjNFhSHQ+1YQohq5rJClpmZSUhIiPN1aGgoGRkZFS4PCQlxLjcYDKSmptK/f3+2bdvGkCFDAPD392fs2LEsX76c/v37M2nSJFfFF27MUVRCxmOvc2bMcyjFZmztW6Dx9yVnyvuk3/0MitWmdkQhRDXSu2rHDocDjUbjfK0oSrnXV1oeFRXF5s2bWbx4MZMmTWLx4sVMmzbNufy+++7j7bffpqCgAH9//3LHljnJKueOn09BgQGoQvaiEvyem4cS5E/pk/HkG8u2Ow3QcgiG1dsJ+iMVe6sIDrrh53C13PF77W66deumdoQ6z2WFLDw8nJ07dzpfZ2VlERoaWm55VlaW83V2djahoaGUlpby66+/cttttwEwbNgw3njjDRwOBx999BETJkxAp9M5t7vw6/PkB6tiycnJbvn5nPAvm1SzW7eWFa6jWKykDX8SfXQrAh6Jc/5hdDA1lTZRUWUrtWuHbvP36E5m0rVTZzQGl/0XUJ27fq+FuFouu7TYt29ftmzZQm5uLiUlJaxevZqYmBjn8oiICEwmk/MvxsTERGJiYtDr9bz22mvs2bMHgJUrV9K1a1e0Wi1r1qwhKSkJgISEBDp37oyPj4+rTkG4meyX30ej1xHw8MhyrfuL6UKDQQO5MxfUYDohhKu47M/RsLAwJk2axLhx47BarcTHx9OpUyfGjx/PxIkT6dixI7NmzWLKlCkUFhbSvn17xo0bh06nY/bs2UydOhW73U5YWBgzZswA4I033uDll19m3rx5BAcH8+abb7oqvnAzRWu2ULTiV0Jm/QuN9kp/n2nQNwrj3Jc/4HfX7ZjaR9ZIRiGEa2gURVHUDlGd5HJK5dz181l2c9mlxZHrL7206Cgo4kTf+wl67B5MnaMuWV7u0iJgfGI8AGeHjMO8az8RP35QaQvOXbnr91qIqyUjewi3l/v2fzG1j7xsEauMz8AbsWflUrJuu4uSCSFqghQy4dasx05T8OVy/McMueptNTot/vEDyZ25AA+7MCFEnSKFTLi13P/8Hz539kMXHHhN23v17YI9r4CSX6WbuhDuSgqZcFuWQycoXrsVv9ibr3kfGq0W38H9yP9wSfUFE0LUKClkwm2dffszfIf0R+vrfV378b65B+ade7EeOVVNyYQQNUkKmXBLtrQMildvxvfOm657X1qTEe9bepL/xffVkEwIUdOkkAm3lDd/Kd639ELrVz0PxPsM6EXhklUoNhmHUQh3I4VMuB1HUQkFX63Ad3C/atunoUk4upBgitduq7Z9CiFqhhQy4XYKv12DsV0r9GH1q3W/3v27U7Bk1ZVXFELUKlLIhFtRFIW8Bd/gM+jGat+3V5/OFP+8DUdhcbXvWwjhOlLIhFsp/W0/SkERpk5tqn3fugA/jG1bUZS0qdr3LYRwHSlkwq2c+2I53jf3rMLAwNfGu28XCr9Z45J9CyFcQwqZcB92B0WJP+MzoKfLDuHVowMlW37HUVTismMIIaqXFDLhNuy5+RiimqOrH+SyY2j9fDBGNaf4ZxlIWAh3IYVMuA17Zi7e/Vw/LYlX9w4ULV/v8uMIIaqHFDLhFhSLDUdhMV69O7n8WF69OlK8dqs8HC2Em5BCJtyCPScPrb8vWpPR5cfS1Q9CF1IP8859Lj+WEOL6SSETbsGelYs2wK/Gjmfq0pbiNZtr7HhCiGsnhUzUera0DJSSUrR+1zfK/dXw6t5OnicTwk1IIRO1XkHiz2gDfEGjqbFjGlo3x56ejS09q8aOKYS4NlLIRK1XlLiuRi8rAmh0Wkyd21D8y84aPa4Q4upJIRO1mu1MNpY/T1z35JnXwtixDcVrt9b4cYUQV0cKmajVilb8ilf39jV6WfE8U+coSjbsRHE4avzYQoiqk0ImarXC5evx6tFBlWPrw+qj9fXBsueQKscXQlSNFDJRa9nPFVK6ax+mG9qqlsHYsTUlm3apdnwhxJVJIRO1Vsm67RjbR6L1NqmWwdQhkuL1O1Q7vhDiyqSQiVqraOWveHVrp2oGY4fWmLfvQbHKcFVC1FZSyEStpNjtFK/bhle39qrm0AX4oQ+vT+kfqarmEEJUTAqZqJVKd+13jnmoNmO7SIo3yn0yIWorKWSiVipaswVTl2i1YwBgbNcSsxQyIWotKWSiVipevQlTV3Xvj51nbNcKc/I+mdZFiFpKCpmodWyZuVhPnMEY1VztKEDZfTJdSD1KU/5UO4oQ4jKkkIlap+SXHZg6tUaj16kdxcnYrhXmzb+rHUMIcRkuLWTLly9n8ODBDBw4kIULF16yfP/+/cTFxTFo0CBeeuklbP+7dLNz507i4uKIjY3l73//O/n5+QCcO3eOCRMmcOeddzJmzBiysmRkck9UvHYbpk5Rascox9S2JSWbflM7hhDiMlxWyDIyMpg9ezaLFi0iISGBJUuWcOhQ+aF+Jk+ezNSpU0lKSkJRFJYuXQrACy+8wJtvvsny5cuJjIzkk08+AWDOnDl0796dlStXMmrUKGbMmOGq+EIlisNB8fodtaajx3nG6BaYd+xBURS1owghLuKyQrZ582Z69+5NUFAQPj4+DBo0iFWrVjmXp6WlYTab6dKlCwBxcXHO5StWrCAyMhKr1UpGRgYBAQEArF+/ntjYWACGDh3Khg0bsFqtrjoFoQLLviNovU3ow+qrHaUcXYN6aLxNWP88rnYUIcRF9K7acWZmJiEhIc7XoaGh7N69u8LlISEhZGRkAGAwGEhNTeWhhx5Cr9fzz3/+85Jt9Ho9fn5+5ObmEhYWVu7YycnJrjotj1CbPx/Ttz+jaxLKwdTyDyAXl7QA4GDq0Wva74X7a1VSDMDh1Kt7yNnUqAGpS3/Acmefa8qghtr8vfYU3bp1UztCneeyQuZwONBcMPWGoijlXl9peVRUFJs3b2bx4sVMmjSJxYsXX3IMRVHQai9tVMoPVsWSk5Nr9edz+o1FePXrgXdU+Xtkh7zLOn60ibr6e2cHU1PLbWf09rmmfRX1zMZ+KpewWvz5Xai2f6+FqC4uu7QYHh5erjNGVlYWoaGhFS7Pzs4mNDSU0tJSfvrpJ+f7w4YNI/V/fzmHhoaSnZ0NgM1mo6ioiKCgIFedgqhhisWKeXsKpo6t1Y5yWca2LTFvT1E7hhDiIi4rZH379mXLli3k5uZSUlLC6tWriYmJcS6PiIjAZDI5L30kJiYSExODXq/ntddeY8+ePQCsXLmSrl27AtC/f38SEhKAsvto3bt3x2AwuOoURA0z79qPPiIMrb+v2lEuS984DHveOWyZuWpHEUJcwGWXFsPCwpg0aRLjxo3DarUSHx9Pp06dGD9+PBMnTqRjx47MmjWLKVOmUFhYSPv27Rk3bhw6nY7Zs2czdepU7HY7YWFhzt6JTz31FM8//zxDhgzB39+fWbNmuSq+UEHJxl0Y27dSO0aFNFotprYtMe/Yg9+QmCtvIISoES4rZACxsbHOXobnLViwwPl1dHQ033zzzSXbde/ene++++6S94OCgpg/f371BxW1QskvO/G5rbfaMSpliGyGedtuKWRC1CIysoeoFRzmUkr/OICxXUu1o1TKGN0C89bdV15RCFFjpJCJWqE0eR/6ZhFofbzVjlIpQ+tmWPYfQSm1qB1FCPE/UshErVC8cRemWt4aA8oe1m4cJgMIC1GLSCETtYJ5028Y29Xejh4XMrRphnnnHrVjCCH+RwqZUJ1isVL6eyrG6BZqR6kSY+tmcp9MiFpECplQXekfqegjQtH6+agdpUqM0S0w79yrdgwhxP9IIROqK9n8O8a2tf/+2Hm68AYopRZspzPVjiKEQAqZqAVKNv2Gsa17XFYE0Gg0Za2y5H1qRxFCIIVMqExxOCjduRdjW/fo6HGeoWUTzDukw4cQtYEUMqEqS+oxNP6+6OoFqB3lqhijmkshE6KWkEImVGXenuJWlxXPM0Q2xbL3MIrVpnYUIeo8KWRCVeYtv2Ns01ztGFdN6+uNLiwYy77DakcRos6TQiZUZd6+x616LF7I2Ka5dPgQohaQQiZUY8vMxZ53Dn3jMLWjXBNDqyaYd8hEm0KoTQqZUE3pzj0Yo1qg0brnj6G0yISoHdzzN4jwCCXb92CMbKp2jGumbxqOPSMHe16B2lGEqNOkkAnVmLfvxhDVXO0Y10yj02GIbErpb/vVjiJEnSaFTKhCsdqw7Dnslj0WL2Ro1VQuLwqhMilkQhWWfYfRhQWj9a3dE2leibF1U3kwWgiVSSETqjAn78PYupnaMa6bsU0zSn/fj6IoakcRos6SQiZUYd6RgsGNO3qcp60fBFotthPpakcRos6SQiZU4SktMo1Gg7FNc0p3SYcPIdQihUzUOHteAfYz2eibNVQ7SrUwtGqCOVkm2hRCLVLIRI0r/f0AhshmaHQ6taNUC2ObZtLhQwgVSSETNc78234MrRqrHaPaGFo1wbL/qIyEL4RKpJCJGle6c69bj+hxMa2PN7rw+jISvhAqkUImalzpH6ke0WPxQsbIpphlhA8hVCGFTNQo25lsFHMpurD6akepVoZIeTBaCLVIIRM1qvT3AxhaN0Oj0agdpVoZWzejVIaqEkIVUshEjTL/th9DS8/p6HGevllDbKczsZ8rVDuKEHWOFDJRo0p/O4ChZRO1Y1Q750j4vx9QO4oQdY4UMlFjFEWh9I9UjK08r5ABGFo2kZHwhVCBSwvZ8uXLGTx4MAMHDmThwoWXLN+/fz9xcXEMGjSIl156CZut7Dmc5ORk4uPjGT58OA888ABpaWkAbN++nV69ejF8+HCGDx/OCy+84Mr4oprZ07NAUdA2CFI7iksYWzeldLt0+BCiprmskGVkZDB79mwWLVpEQkICS5Ys4dChQ+XWmTx5MlOnTiUpKQlFUVi6dKnz/enTp5OYmEhsbCzTp08HYM+ePTz88MMkJiaSmJjIzJkzXRVfuEDpH6kYWzf1uI4e5xnaNMf8m4yEL0RNc1kh27x5M7179yYoKAgfHx8GDRrEqlWrnMvT0tIwm8106dIFgLi4OFatWoXFYuGpp54iOjoagKioKNLTy0YWT0lJYePGjcTGxvL3v//d+b5wD+bfDqBvHqF2DJfRhdQDRcF2KkPtKELUKS4rZJmZmYSEhDhfh4aGkpGRUeHykJAQMjIyMBqNDB8+HACHw8HcuXO57bbbAPD392fs2LEsX76c/v37M2nSJFfFFy5Q+vsBj+yxeJ5Go8EY1Vy64QtRw/Su2rHD4Sh3CUlRlHKvr7TcYrHw/PPPY7PZePTRRwGYNm2ac/l9993H22+/TUFBAf7+/uWOnZycXO3n40nU+nwCfttHdq9olNTUq962uKQFAAdTj17TsQ9ecMxWJcUAHL6GHFdiCPQmb8U6zE0Cq33f10L+L7het27d1I5Q57mskIWHh7Nz507n66ysLEJDQ8stz8rKcr7Ozs52Li8qKuKxxx4jKCiIDz/8EIPBgMPh4KOPPmLChAnoLhg1XXeZEdTlB6tiycnJqnw+tjPZnFQgsle3a7pHdsi77PvcJirqqrc9mJpabjujt8817+tKSs1Q+P3PtK8FP4Nqfa+FqGkuu7TYt29ftmzZQm5uLiUlJaxevZqYmBjn8oiICEwmk/MvxsTEROfyyZMn06xZM+bMmYPRaCwLqtWyZs0akpKSAEhISKBz5874+Pi46hRENSpN+RNDpOd29DjP0Lopln1HZCR8IWqQy1pkYWFhTJo0iXHjxmG1WomPj6dTp06MHz+eiRMn0rFjR2bNmsWUKVMoLCykffv2jBs3jn379rF27VoiIyMZOXIkUHZ/bcGCBbzxxhu8/PLLzJs3j+DgYN58801XxRfVrHR3KoZmjdSO4XIXjoRv6lz9LT4hxKVcVsgAYmNjiY2NLffeggULnF9HR0fzzTfflFverl07Uiu4d9G6dWsWL15c/UGFy5X+dgBTx9Zqx6gRxtbNMCfvk0ImRA2RkT1EjSjd8yeGFp7bY/FChtZNMW/frXYMIeoMKWTC5ez5BTjOnkPXsIHaUWqEsU1zzDulC74QNUUKmXA5y55DGFo0RqOtGz9u+sbh2LPPYs/NVzuKEHVC3fjNIlRVmvIn+uae39HjPI1O+79W2V61owhRJ0ghEy5X+scBDHWokMH/ZoyWQiZEjZBCJlyuNKXudPQ4z9imGebtKWrHEKJOkEImXEoptWA7loahSbjaUWqUIao5pX+kotjtakcRwuNJIRMuZTl4HF3DUDQmo9pRapQuwA9dkD+W1GNqRxHC40khEy5VuvdQnbs/dp6hTXNK5T6ZEC4nhUy4VGnKQfRN69ZlxfOMbZpRsk0ejBbC1aSQCZey7P4TgwdPplkZY3QLzDv2qB1DCI8nhUy4jKIoWPYfrrOFTN+kIfbMXOw5eWpHEcKjSSETLmNPzwKtFm29ALWjqEKj02KMao45We6TCeFKUsiEy5TuPYyhRYTHz0FWGUPrZpi3y+VFIVxJCplwGcu+w+ibNFQ7hqqM0S0o2fqH2jGE8GhSyITLlKYcxNCsjheyNs2xpPwpM0YL4UJSyITLWPYerrPPkJ2n9fVG3zCE0pSDakcRwmNJIRMuoZRasJ1MR9+4bj5DdiFDVHMZd1EIF5JCJlzCcvA4ukahaIwGtaOozhjVnJLNcp9MCFeRQiZcwrLvMIamdfv+2HnGtq0wb09BURS1owjhkaSQCZco3XtILiv+jy6kHhqdFtvRNLWjCOGRpJAJl7DsPVTneyyep9FoMLZtKeMuCuEiUsiES1gOHEXfrG73WLyQMao55s2/qx1DCI8khUxUO/vZczgKS9CFBqsdpdYwtmtFyVZpkQnhClLIRLWz7DuMoUWjOj001cX0TRpizzmLLTNX7ShCeBwpZKLaWfYfqfNDU11Mo9NibNcKswxXJUS1k0Imql3p3sPoG4epHaPWMUa1oGTzb2rHEMLjSCET1c6y/zAG6ehxCVOHSEo2/a52DCE8jhQyUa0URcFy8Bj6pvIM2cUMLZtgO5GOPa9A7ShCeBQpZKJa2U9nojEY0AX6qx2l1tHodRijW8h9MiGqmRQyUa1K9x+t8yPeV8bYtiUlG+U+mRDVqUqF7Mknn2Tz5s2uziI8gOXAERmaqhLG9pGUbExWO4YQHqVKhez222/ngw8+YNCgQXzyySfk5eW5OJZwV5a9h9A3kR6LFTG2bor1yCns+XKfTIjqUqVCNmzYML788ks++OADcnJyiI+PZ/LkyezeXflIBcuXL2fw4MEMHDiQhQsXXrJ8//79xMXFMWjQIF566SVstrJZdJOTk4mPj2f48OE88MADpKWVDbZ67tw5JkyYwJ133smYMWPIysq62vMVLmbZf0R6LFZCYzBgjG6JeYvcJxOiulT5HpnD4eD48eMcO3YMu91O/fr1efXVV3nvvfcuu35GRgazZ89m0aJFJCQksGTJEg4dOlRuncmTJzN16lSSkpJQFIWlS5c6358+fTqJiYnExsYyffp0AObMmUP37t1ZuXIlo0aNYsaMGdd63sIFFLsd6+FT6JvIpcXKlN0n26V2DCE8RpUK2ezZs+nfvz//93//x+DBg1m9ejXPP/88X3755WVbWgCbN2+md+/eBAUF4ePjw6BBg1i1apVzeVpaGmazmS5dugAQFxfHqlWrsFgsPPXUU0RHRwMQFRVFeno6AOvXryc2NhaAoUOHsmHDBqxW6zWfvKhetuPpaIP80Xp7qR2lVjN1ak3JBrlPJkR10VdlpdzcXBYsWOAsLuf5+Pjw9ttvX3abzMxMQkJCnK9DQ0PLXYq8eHlISAgZGRkYjUaGDx8OlLUC586dy2233XbJNnq9Hj8/P3JzcwkLK39PJjlZfklUxlWfj2Hzbkz1/TmYmlrt+y4uaQHAwdSj17T9hZlalRQDcNgFOavEYcf36Cl2/bwBJcDXpYeS/wuu161bN7Uj1HlVKmR2u/2SIjZx4kTee+89brrppstu43A4yg0aqyhKuddXWm6xWHj++eex2Ww8+uijlz2GoihotZc2KuUHq2LJycku+3xyf9mDNaoVAVFR1b7vQ946ANpcw74PpqaW287o7XPN+6ouOR3bEF6k4HeL635WXfm9FqI2qbSQvfLKK2RkZJCcnExu7l+jdttsNk6ePFnpjsPDw9m5c6fzdVZWFqGhoeWWX9hZIzs727m8qKiIxx57jKCgID788EMMBgNQ1qrLzs4mPDwcm81GUVERQUFBVT9b4VKWvYcxtIhQO4ZbMLVrRcn6HfgN7a92FCHcXqX3yOLj4xk4cCB+fn4MGjTI+W/YsGEsWLCg0h337duXLVu2kJubS0lJCatXryYmJsa5PCIiApPJ5Lz0kZiY6Fw+efJkmjVrxpw5czAajc5t+vfvT0JCAgArVqyge/fuziIn1Gc5cBR9Uxn1viqMndpQ8qtc9hOiOlTaIuvYsSMdO3bkxhtvvOQ+1JWEhYUxadIkxo0bh9VqJT4+nk6dOjF+/HgmTpxIx44dmTVrFlOmTKGwsJD27dszbtw49u3bx9q1a4mMjGTkyJFAWUtswYIFPPXUUzz//PMMGTIEf39/Zs2ade1nLqqVYrVhO56GQUa9rxJDiwjs2XnYTmeibxR65Q2EEBWqtJA99dRTvPvuu/ztb3+77PLly5dXuvPY2FhnL8PzLmzJRUdH880335Rb3q5dO1IruAkfFBTE/PnzKz2mUIf16Cl0IcFoTMYrryzQaLUYO7eheEMyAffeqXYcIdxapYVs/PjxALz88ss1Eka4L8uBo/L82FUytY+k5OftUsiEuE6V3iPr0KEDAD179qRhw4b07NmT4uJiduzYQdu2bWskoHAPlgNH0UfIZcWrYeocRcmvySiKonYUIQB49913nf0QrsX69et59913qy9QFVXpgeipU6eyYMECDh8+zJQpUzh16hQvvviiq7MJN2LZd1haZFdJF94A9HosB67t2TghqttTTz3FiBEjrnn7lJQU8vPzqy9QFVXpObI9e/bwzTff8PHHHzNy5EieeeYZ4uLiXJ1NuBHLgaP43NZH7RhuRaPRYOrchpJfdmBq21LtOMIDPPPMM7Rv356HH34YgEWLFrF161ZCQ0P5448/KCoqQlEUpk+fTrdu3Xj++efJy8vj5MmT3HzzzeTk5NC6dWseeeQRvvnmG5YsWYLVaiU/P5/x48czevRovvvuO9asWYNWq+X48eN4eXnxxhtvUFhYyOLFi7Hb7fj7+zNp0qQaO+8qtcjOP3i8adMmevfuDYDZbHZpMOE+FIsV26kz6COk993VMnWKovinrWrHEB5i1KhRLFu2zPl62bJlREVFkZmZyZIlS1ixYgUjR44s1+nObDbz448/MnnyZOd7RUVFfP3113z88cckJCQwe/Zs3nrrLefyHTt28PLLL/PDDz/QuXNnPv74Yzp37sy9997L4MGDa7SIQRVbZE2bNmX8+PGcOnWKnj178swzzxCl4qgIonaxHjmFLrQ+GqM803e1TJ3akDfvKxzmUrReJrXjCDfXq1cvSktLSUlJwdvbm9zcXB5//HGOHj3K4sWLOXnyJNu2bcPX96+h0S43+ouvry/z58/nl19+4dixYxw4cIDi4mLn8vbt2xMeXnYroV27dqxZs8b1J1eJKrXIZs6cydChQ/niiy8wGAx0796df//7367OJtyEPAh97bR+PhiaNcK8LUXtKMIDaDQa4uPjSUxM5NtvvyU+Pp5ffvnFOczfrbfeyn333VduGx8fn0v2c+bMGUaMGEFaWhrdunXj6aefLrfcy+uvgcE1Go3qHZaqVMh8fHzo3r07586dY+/evXTq1IkjR464OptwE6UHjsplxetg6tyG4rVyeVFUj5EjR7Ju3TqSkpKIi4tj06ZN3HLLLYwePZoOHTrw008/YbfbK93Hnj17CA4O5vHHH+emm27i559/BrjidjqdzjmvZE2qUiF79913GTx4ME888QRPPvkkTz75JBMnTnR1NuEmLPsPo28sPRavlalLtBQyUW1CQkJo164dbdq0ISwsjHvvvZft27cTGxvLyJEjadKkCadOncLhcFS4j/OjOd1xxx3ceeedpKenExwczPHjxys9du/evdm4cSOvv/56dZ9WpTRKFdqEAwYM4KuvvrrqYarUICN+V84Vn8+J3qMJenI0huauGzB4xRNlo98Pnlv5X4SXc8no90+UPehvmVv5eKE1RbE7yHhoCk1+/W+1Dlcl/xdEXVGlFlnDhg3dooiJmqdYrFhPnpHxAq+DRqfFdEM0xWu3qR1FCLdUpV6Lffr04c033+TWW28td5Ovffv2Lgsm3IPl8En0YdJj8XqZOkdRtHozAWNjr7yyEKKcKhWy7777DoBVq1Y539NoNKxdu9Y1qYTbsKYekx6L1cDUtS35ny5DsdrQGKr031II8T9V+h+zbt06V+cQbqr0wFG5rFgNdEEB6BuGYt6egveNN6gdRwi3UqV7ZEVFRUybNo0HHniAvLw8pk6dSlFRkauzCTdg2S9jLFYXU5doilZvVjuGEG6nSoVs+vTp+Pv7k5OTg8lkorCwkKlTp7o6m3AD1tRjGKSQVQuv7u0oXr1J7RhCuJ0qXVrcv38/M2fO5JdffsHb25tZs2YxdOhQV2cTtZyzx6I8DF0tDJFNsefkYz1+GkOzRmrHETXsWJe7sKdlVvt+dRGhNP/926va5vnnn6dnz54VDg4fFRVFamoqu3fvJikpqdw4jWqoUiHTass33Ox2+yXvibrHeuQU+tBg6bFYTTRaLaaubSles4XAv92ldhxRw+xpmTRKeK/a93t6hOsGrzh06BA5OTku239VVamQ9ejRg7feeguz2cyvv/7Kl19+Sa9evVydTdRyltRjcn+smnl1a0fRyl+lkIkapSgK//nPf1i/fj2hoaHY7XZ69uxJQkIC//3vf3E4HLRv355XXnkFk6lscOtz587x3nvvUVxczIcffsjYsWN58cUXycjIIDMzkz59+jBjxgw0Go3L81epWfWvf/0LHx8f/P39mTNnDtHR0Tz77LOuziZqudIDR2RW6Gpm6hKNeedeHIXFV15ZiGqSlJTEvn37+OGHH3j33Xc5ceIEJSUlLF26lMWLF5OYmEj9+vX55JNPnNsEBAQwceJEBgwYwGOPPcb69etp27YtS5YsISkpiR07drB3794ayX/FFtmaNWv45JNPSE1NxcvLi6ioKLp27eqsyqLusu4/gqFNc7VjeBStjzfG6BYUr9+B39D+ascRdcT27dsZOHAgBoOB4OBgYmJiUBSF48ePc/fddwNgtVpp165dhfsYOnQou3fv5rPPPuPIkSPk5eWVm/rFlSotZCtXrmT27NlMnDiR6OhoNBoNKSkpzJgxg9LSUgYOHFgjIUXtZEk9hs+tvdWO4XFMXcsuL0ohEzXl4qlY9Ho9drudO++8kylTpgBlj2FVNvr9F198QVJSEnfffTd9+/bl4MGDNTa9S6WXFj///HM+++wzhg4dSmRkJK1atWLEiBF89NFHfPrppzUSUNROitWG7WS69Fh0Aa8eHSheswXlClNmCFFd+vTpw8qVK7FYLOTn5/Prr78CZVfkcnJyUBSFV199lf/+97/ltrtw2pZNmzZxzz33MGzYMEpLSzlw4EClI+xXp0pbZEVFRTRqdGk34BYtWlBaWuqyUKL2sx49hS4kGI3JqHYUj6MPq48uOBDzjr149+6kdhxRB9x2222kpKQwdOhQGjRoQKtWrfD39+eJJ57ggQcewOFw0LZtWyZMmFBuu06dOjF37lxmzZrFAw88wKuvvsrHH3+Mn58fN9xwA6dOnaqR/JUWMp1OV+EytWcEFeqyyBiLLmXq1p6iFRukkNUhuohQl3SV11XxqsmkSZOYNGnSJe+PGjXqkvdSU1OBskbNmjVrnO8nJSVdY8rrI6OTimtiOXAUfaMQtWN4LK+eHch7fxH1X/tHjXRfFuq72oeWxV8qLWSpqal07dr1kvcVRcFisbgslKj9LPuPYGjVRO0YHsvQqglKiRnrwWMYo1qoHUeIWq3SQnZhk1GIC1kOHMXnlh5qx/BYGo0Gr54dKVzxK8FSyISoVKWFLCLCdVPXC/elWG3YTpyWh6FdzKtHBwq+XUPwpHFqRxGiVpMBE8VVsx5Lkx6LNcDYPhLbsdPY0jLUjiJErSaFTFw1y4GjMsZiDdDodZh6tKdoxa9qRxGiVpNCJq6aJfWYSx6EVhwKpXl2SnJtKA55vAMou0+2fL3aMYSo1Vza/X758uV8+OGH2Gw2HnjgAcaMGVNu+f79+3nppZcoKiqie/fuvPbaa+j1f0WaM2cOOp2OJ598EigbD+zJJ58kPLysNdCuXTtmzpzpylMQl2HZd7jaeiw67Aqnfyni+MpCcvea0eoAjQa72UG9diaaDPSnyW2+QMXPNHoyry7R5L23EHv2WXQN6qkdR7jQf5umUnjSWu379Wti4IETUdW+3/PGjx/P9OnTCQtT7565ywpZRkYGs2fP5rvvvsNoNHLvvffSq1cvIiMjnetMnjyZ6dOn06VLF1588UWWLl3K6NGjKSgoYObMmfz444/87W9/c66/Z88eHn74YR599FFXxRZVUF09FrNTSvj9rWx0Ri3hfbxpOcIPvW9ZwbKbHeQdLOX4D+dI/ewsFq9mGAPq3gUEjcmI6Ya2FK3aRMD9MpmtJys8aWXEz82rfb8Jtxyr9n1eaMGCBS7df1W4rJBt3ryZ3r17ExQUBMCgQYNYtWoVTzzxBABpaWmYzWa6dOkCQFxcHO+99x6jR49m7dq1NG/enIceeqjcPlNSUsjOzuaHH34gIiKCV155hYYNZXSJmlQdPRYVReHAZ2c5mlhAyzh/6nf0gose+tV5aanfyZv6nbzJP1jK1vlWjAFa7BbQGetWQfPu1ZHCxHVSyIRLzZ8/n++//x6dTseNN97I6NGjefLJJ2nZsiWHDh2iXbt23HDDDSxbtoz8/HzmzZtHq1atGDBgAJ9//jkhISG89tprJCcnYzAYePzxxxk8eDADBgygU6dO7N+/n0WLFrF+/Xo+/fRTNBoN7du35+WXX8bX15ebbrqJQYMGkZycjE6nY86cOTRpUrUrPy77jZCZmUlIyF8jP4SGhpKRkVHh8pCQEOfyESNGMGHChEuGyPL392fs2LEsX76c/v37X3Y4FeFa1qOn0DWod809Fh12heR/Z3F6QzFd/lmf+p28LyliFwtsY8KvsQGbWWHzs2ewFtetwXRN3dtj3rEHe36B2lGEh/rll19Yt24d3377LcuWLeP48eP8+uuvpKamMn78eBITE9m1axdpaWksWbKEoUOHsmTJknL7+OKLLyguLmblypV8+umnzJs3zzlwRkxMDElJSWRnZzN//ny++OILli9fjre3N3PnzgUgKyuLPn36kJCQQI8ePVi4cGGV87usReZwOMoNraMoSrnXV1p+OdOmTXN+fd999/H2229TUFCAv79/ufWSk5OvN75Hu57Px7DxD4z1/Dn4v7HWrobigFOf6LDlaGgQa+dMXjHkVW3bUmtzdPWgVClm3T+O0fSfNnRXUUsvzNuqpGyOpMPXcA5qMbVoyN75C7HednWXdOX/gut169ZN7QjXbevWrQwZMgRvb28A7rrrLhISEmjQoIFzDrLw8HD69OkDQKNGjS4ZEHjHjh3cfffdaLVaQkJC+PHHH53LOnfu7FznlltuoV69svu999xzDy+88IJzvX79+gHQunVrdu7cWeX8Litk4eHh5YJkZWURGhpabnlWVpbzdXZ2drnlF3M4HHz00UeXtNQuN7CxJ/xguUpycvJ1fT6563ZjjW5JQNTV3zze/X422nOldHksCO1VXh48+b+JXDs/0IhDS8+R86lCnzfD0emvPA7hwdRU2lyQ1+jtA1Duvdqu+NZ8LL8fpuFzf6/yNtf7vRZ1x+WmW7HZbBiN5f9arGwgeb1eX64xcvz4ceetn/MTMV98HEVRnNPAXLjexfOjXYnLLi327duXLVu2kJubS0lJCatXryYmJsa5PCIiApPJ5PyLMTExsdzyS4JqtaxZs8Y5unJCQgKdO3fGx8fHVacgLsOy7/A1PUN2bHk+Z7YUE/1g4FUXsXK0GlqNCsBhcbD73ew6MwuDV8+OlGz5HUdBkdpRhAfq3bs3P/74I2azGZvNxrfffkvv3lc3aW6PHj1YsWIFiqKQk5PD/ffff8mYvD179mTdunXk5eUBsHTpUnr16nXd+V3WIgsLC2PSpEmMGzcOq9VKfHw8nTp1Yvz48UycOJGOHTsya9YspkyZQmFhIe3bt2fcuMqH4nnjjTd4+eWXmTdvHsHBwbz55puuii8qYDlwFJ8BPa9qm7OppexdkEfHfwSj97n+bvQanYY2Y4LY/X4Ox38soPnQgOveZ22n9fXG2CGSoqRN+MfLzOyeyK+JwSU9DP2aGK64zi233ML+/fu56667sNls3HTTTdxyyy18/vnnVT7O6NGjmT59OsOGDQPg5Zdfxs/Pr9w60dHRPProo4wdOxar1Ur79u157bXXru6ELkOjeNiftHI5pXLX8/koVhtHWwwk/POZVe7sYStxsO5vp2hymx8NbvC+puMCbJ5T1vLu+3Sx872SM1ZSPjjLTXMaEtiy4jyXXFp8YjwAlrnqdxu+GsXrtmE5cJSGC9+o0vryf0HUFXWrH7O4LtYjJ696jMWUD3LwizBeVxGriHe4gWZD/dkxLRO7pWamVFeTV8+OlGz6TS4vCnERKWSiyiz7j17VrNBZv5dwZlMJLUf6X3nlaxTWwwtTkJb9n5x12TFqC62fD6YOrSlK2qR2FCFqFSlkospKDxyt8oPQdouD397KouVIf3TeLvwx02hoFR/AiaRCzqaWuu44tYRX704UfivzBApxISlkosqupsfin0vy8W6gJ7ijl4tTgcFPR/Oh/uz6TxZ2m0fd8r1EWe/FP+ThaCEuIIVMVJnlwBEMVbi0WJxp5fDSfJoPc90lxYuFdPNC56Xh6Hf5NXZMNWh9vTF1aiNTuwhxASlkokqUUgu2tEz0ESFXXHfv/LOE3+iLV32XTq5QnkZDyxEBpH6RR0m27crruzGvPl3k8qIQF5BCJqrEcugE+oYN0BgqfyblbGopWbtKiLil5h9U9w7TE9rLh70f5db4sWuSV48OmHfuxZ6Tp3YUIWoFKWSiSiwHjmJo2uiK6+39MJcmA/3QmdT50Wpymy+ZO0rI3W9W5fg1QettwqtbO4p+3KB2FCFqBSlkokos+w6jb1z5rNBZu0ooOm0lrFf1PzNWVTovLU0H+ZHyfo5HD1/l1fcGCr5OUjuGELWCFDJRJZa9h9FX0iJTFIW9C8paYxrdlQfydaWwnt5Y8h2kbyq+8spuyqtrWyz7DmNLz7ryykJ4OClkokosB49haFpx1/vMnSVY8u2E3OD67vZXpNPQbIgfez/MxWH3zFaZxmjAq2dHChLWqR1FCNVJIRNX5CgqwZ6Ziy68wWWXK4rCgU/P0uRWP1C5NXZevbYm9H5aTqzy3OetvG7qSuE3q9WOIYTqpJCJK7IcPIa+STiaCuYiyv7NjDnXToMutaA1dp5GQ7M7/dj/aR4Oy5VXd0emjq2xncrAcvik2lGEUJUUMnFFln1HKh1jMfWLPCJu8a01rbHz/Jsb8Y3Qk7vh+qeOqY00Oh3eN94gz5SJOk8Kmbgiy/7D6BtffozFs6mlFBy3EtpNvZ6KlWk60I+cVVqsJXa1o7iEd79uFCxN8ugemkJciRQycUWlew5haHb5HosHF+bRqL8PGn3tao2d5xthwNTYwZFvz6kdxSUMrZuCzUbp7wfUjiKEaqSQiSsqexj60kuLhaetZP9WoupzY1UR2MvB4a/PYS3yvFaZRqPBq183CpasUjuKEKqRQiYqZc8+i1JqQdsg6JJlh7/OJ6y3Dzqv2v1jZAiGwDYmDntoq8wnphuFy9ai2Dx7jEkhKlK7fwMJ1VkOHMXQIgKNpvylQ0uhnZNrCml4Y82PqXgtmt7uy+FvPLNVpm8Uij68AcXrd6odRQhVSCETlSrdd+SyHT2O/1hAvbZeGIPco0egV6ieetEmDn/nma0y75tuoGDxCrVjCKEKKWSiUpa9f15yf8xhVzjy7Tka9avd98Yu1vhWXw5/fQ7F4Xk9/Lxu6krx2m04CorUjiJEjZNCJipVuucQ+uYR5d47s7kYg78Wv6ZGlVJdG+8wPUFRJsxZnnd5URfgh6lTGwq/X692FCFqnBQyUSHFbsf654lLWmSHvs6nYV/3uDd2sSa3+VCSZfPIVpl3TDcKFv2odgwhapwUMlEh67HTaAP90Pr+dQkx/6iFwhNW6neuRcNRXQXvMAN6by3mbM9rlXl1a4/l4DGsx0+rHUWIGiWFTFTIsu8whosuKx5Zlk9479r7AHRVmOppKcm0edxoHxqDHu+busozZaLOkUImKlS67zD6Jn9N3WItspO2roiwPu7VyeNiWqMWvbeWo8s8b2R871t6UvDVChSHQ+0oQtQYKWSiQpbdBzE0++v+2MnVhQRFmTAGukeX+8qY6mk5tDTf41plhpaN0RgNmLf8oXYUIWqMFDJRIcv+I84ei4qicCShgPDe7t0aO09r1BLYyuhxrTKNRoP3zT0498VytaMIUWOkkInLchQUYc86i75hCAA5e8zYSx0ERrpXl/vKNL7d1yNbZd79u1OUtAmKStSOIkSNkEImLqt032H0zRuh0ZX9iBxdVkB4Hx/QuG8nj4v5hBsIjDRy1MNG+9AF+mPqEoVx/S61owhRI6SQicuy7D3snLrFkm/nzNZiQru5Z5f7yjS5zZdDSz1vDEafAb0wrtisdgwhaoQUMnFZpbsPOh+EPr6qgPodvND7uX8nj4t5hxsIijJx+Jt8taNUK1PnaLS5BZSm/Kl2FCFcTgqZuKzSPX+ibx6Boigc+76AMA/p5HE5TW735fA3BVjOeU6rTKPTYu0ezbnPv1c7ihAu59JCtnz5cgYPHszAgQNZuHDhJcv3799PXFwcgwYN4qWXXsJ20XxKc+bM4f3333e+PnfuHBMmTODOO+9kzJgxZGVluTJ+naXY7VgPHsfQvBE5u80oikJAC4PasVzGK0RP/Y4mDn7lWa0yW7doCr/7CYd0+hAezmWFLCMjg9mzZ7No0SISEhJYsmQJhw4dKrfO5MmTmTp1KklJSSiKwtKlSwEoKCjgxRdf5NNPPy23/pw5c+jevTsrV65k1KhRzJgxw1Xx6zTrkVNo6wWg9fXmaGIBYb08q5PH5TQZ6MexH85RkuM5k1Mq9fwxRLegMGGd2lGEcCmXFbLNmzfTu3dvgoKC8PHxYdCgQaxa9dfQOWlpaZjNZrp06QJAXFycc/natWtp3rw5Dz30ULl9rl+/ntjYWACGDh3Khg0bsFqtrjqFOqt0z58YWjbGcs5zO3lczBikI6yHNwf+e1btKNXK97be5H/6ndoxhHApvat2nJmZSUhIiPN1aGgou3fvrnB5SEgIGRkZAIwYMQKg3GXFi7fR6/X4+fmRm5tLWFj5iR+Tk5Or9Vw8zZU+H6/VG9D4GkhdeBSv5pCeexpyayZbRcylzQE4eeLkNW1/4XaBpebL7ysKTv5Xj77bWUzhnjE6/gl/PT4nz/D74gTsrZuoHccjdevWTe0IdZ7LCpnD4UBzweUoRVHKvb7S8qpQFAWt9tJGpfxgVSw5OfmKn8/pGZ/j1a8bh+d50WKoH4FNTTWUrmInTWUZmjS9+l/GJ0+cLLedyeRV4b4MtxVStNpOxxnhlyxzNwdTU2nTti0Fd8YQuGU/ofeOUDuSEC7hskuL4eHh5TpjZGVlERoaWuHy7OzscssvJzQ0lOzsbABsNhtFRUUEBQVVb/A6TlEUSvceotgegr3EQWArzxnJoyoa3uTD2QMWsnd7TgcJn9t7U7j8F+x5njUclxDnuayQ9e3bly1btpCbm0tJSQmrV68mJibGuTwiIgKTyeS8zJWYmFhu+eX079+fhIQEAFasWEH37t0xGDy3N50a7GeywaFwfIOW0J7eoPXsTh4X0xq1NB/sR8p7uR4z+aYuKACvbu04J5NuCg/lskIWFhbGpEmTGDduHCNGjGDo0KF06tSJ8ePHk5KSAsCsWbOYOXMmd9xxB8XFxYwbN67SfT711FP8/vvvDBkyhEWLFjF16lRXxa+zSncfRN+iCac3FBPWw3OfHatMg65eKIrC8ZWe04LxveMmzn3yLYrdc56VE+I8l90jA4iNjXX2MjxvwYIFzq+jo6P55ptvKtz+ySefLPc6KCiI+fPnV29IUU7pH6mUKCEERhoxeMB0LddEo6HlcH/2LThLoxhfjP7u/zkYopqj8fai+Ket+A66Ue04QlQrGdlDlFP6234yjgURWkdbY+f5NjUS3MHE3gUqd9esJhqNBt87+5H34RK1owhR7aSQiXJKfvuTvLwG1Gurfk9FtTUb7Ef6hmJy95nVjlItvG+8AUvqUUr3HVY7ihDVSgqZcLJnn8WeX0xQj3A0urrVyeNy9D46WsT6s+uNbOwWh9pxrpvGoMd30E3kfbBY7ShCVCspZMKpeOcBCq2hhPX0UTtKrdGgqxfGQC2pn+epHaVa+A66kaIVv2LLyFE7ihDVRgqZcDqzZDdWv3C8QlzaB8i9aDRExgdwbHmBR1xi1Ab44h3TjfwFFXeyEsLdSCETToW/pOAV3VTtGLWOIUBHyzh/dkzL9IgJOP1ib+bcfxNxFBarHUWIaiGFTABw7pgF07lj+PdurnaUWql+Z28CWxnZ9WY2iuLeD0rrwxtg6hRF/mcJakcRolpIIRMAHHj/CHq9FW2D+mpHqbVajPCn4KiFQ0vdf94yvxEDyP9gMQ5zqdpRhLhuUsgEDrtC+qI/0DRqfNUDN9clWoOW6Afr8edX+ZzZ4t6X5QwtG6NvHkHB4pVqRxHiukkhE5xIKqSe6Rj6Fo3VjlLrmYJ1RD8YRPLMLLfv/OF/122cnf0FitVzJhMVdZMUMsHe+bk08D+GJkLmq6oK/+ZGWt8TwNYXM8j7030vzRmjW6IPq0/B0iS1owhxXaSQ1XFF6VZO/1KI6dxRNI2kkFVVvfZetBzpz+bJZ8hOcd8pX/zib+fsrE+lVSbcmhSyOm7/p2dp1rMQTCY0/v5qx3Er9Tt7E3lPINteyuT4KvccKd/UPhJdSD25VybcmhSyOkxxKOz7+CxNWp1CGyHPj12Lem1NdHisHgc+O8vOf2diKXC/58z877mT3Lc+RSm1qB1FiGsihawOO7WuCK1Jg1fBIWjYSO04bsunoYHOk+pjL3Hw0wOnOPp9Pnar+4zNaIxugb5JOPn/TVQ7ihDXRMYiqsP2fJBDszv9UDbuQ3vrILXjuDWdSUvLuwIJ7WHhxKpCDnyWT7MhfkT09yWgpRFNFWbathbZKUq3UXzaSvEZOyVZNkrz7FgLHTjsoNWBwV+LT6ge/2ZG6kUb8W1sqJZHJgLuG0zO9I8IGD0ErZ+MtSncixSyOqo4w8aptUXc9v8a4Pj6FJqGEWpH8gh+TY20mxBM0WkrWTtK2PpiBrYSB/4tjPg1MmAM1KI1aHBYFSyFDiy5dkqy7BRn2rCXOvBuoMergQ5jkA5TgBavEB2+EXq0eg2KQ8FWrGDOtnN2fwF7PrSi1UGj/r40HeJPYAvjNec2tIjA1LE1eR8sJvjZh6vxExHC9aSQ1VH7/18uDWN80WccxhrWEI3BoHYkj+LbyIDvcAPNh4Ml305xuo3SXDvWYgeKDTR6MHhr8I4yEtpLh1ewDoO/Fq6mdaUoFJ+xkbXLzKZJ6QS2NtHub/WoF3Vtc8n53zeY7OfeIeCB4ejDZIQX4T6kkNVBikNhz/yzdHuxAfZ9++X5MRczBuowBuqqf8caDT4NDTQbYqDJQD8yt5ew5bkzhPf2ocNjwVe9O31Yfbxv7kHuf/6P0NnPVX9eIVxEOnvUQSdWF2Lw01IvyoRyYC+axtJj0d1pDRrCb/Sh63MNsJY4WPtwGoX7rv7emX/8QIp+3CCzSAu3IoWsDkqZ+79OHoqCI3UfmibN1I4kqonOW0vLkQG0vjuAtE/1HPhv7lWN1q/188E/fiDZL77r9qP8i7pDClkdU3oG0jcV0/hWX5TTp8BoQhMYpHYsUc0Co0yE32cjbX0xO1/PxG6p+uMAPnfciC0tk6IVv7owoRDVRwpZHZPxrY7Gt/qh99biOLAXbVNpjXkqnS90eKwe5lw7W54v6z1ZFRqdjsCHRpAz5T0cJe47lqSoO6SQ1SF2i4PMBB3Nh/oB4NiXAo2lkHkyrVFL1NggtEYNm589U+ViZuochaFFBGff+9LFCYW4flLI6pDD35zD1FAhoHnZ80aO/XvRNm2ubijhchqdhtZ3B2Dw07LlhTNVvswYMG445xZ8i/XIKRcnFOL6SCGrQ/54N4egfmVjASo5WVBcBCGhKqcSNUKrITI+AK1Ow45XM3HYr9yRQxdSD78RA8h69m3p+CFqNSlkdUTWrhIKT1jx71T2C8mxLwVNs5ZotPIjUGfoNLS+LxBzrp3d7+VUqTj5xt6M9UQ6RQnraiCgENdGfovVEb/Pzqb5MH80/3su157yBxq5rFjnaPQaoh8IImtnCYe/PVeF9XUEPXo32S+9iz3PPaeqEZ5PClkdUJxh49j3BTQf7Od8T9m3G23zliqmEmrReWtp+0g9Dn6ZR8b24iuub4xugalHB3JenVcD6YS4elLI6oA983No1N/HOUySkpuDUlAAoWEqJxNqMQXriBoXRPKMLArTrjwPWcCYoRSv2ULJpt9qIJ0QV0cKmYezmR2kzMul5cgA53uOlN/RtIiU+2N1XEBLI41v92Pri5lX7Jav9fUm8G93kfnUTBzF5hpKKETVyG8yD/fnV/kEtjIScMEUH/Y/ktHIZUUBNLzRG58wPb+/k33Fzh9ePTtiaB5B7r8/rqF0QlSNSwvZ8uXLGTx4MAMHDmThwoWXLN+/fz9xcXEMGjSIl156CZvNBsDp06cZM2YMd9xxB4899hhFRUUAbN++nV69ejF8+HCGDx/OCy+84Mr4bk9RFH57K7tcawxFwZHyO9qWkeoFE7WHRkPLu/zJ3Wfm+Iord+YIeCSOgm/WYN65twbCCVE1LitkGRkZzJ49m0WLFpGQkMCSJUs4dOhQuXUmT57M1KlTSUpKQlEUli5dCsBrr73G6NGjWbVqFR06dOCDDz4AYM+ePTz88MMkJiaSmJjIzJkzXRXfI5xIKsRhVwjt4eV8T5+dWfZFsMw3JcroTGWjf+z96Cz5Ryu/X6YL8CPwkZFkPD5dhq8StYbLCtnmzZvp3bs3QUFB+Pj4MGjQIFatWuVcnpaWhtlspkuXLgDExcWxatUqrFYrO3bsYNCgQeXeB0hJSWHjxo3Exsby97//nfT0dFfF9wi73sii1V0BaC6YrNH70AG0LVuXe08In3ADzYb6s+OVTGzmyu+Xefe9AUNEKLkzF9RQOiEq57JClpmZSUhIiPN1aGgoGRkZFS4PCQkhIyODs2fP4ufnh16vL/c+gL+/P2PHjmX58uX079+fSZMmuSq+28tMLiHvgIXGA3zLve+duhdNq9YqpRK1WVgPL7zD9KTMzbniugHj76JgaRLm7Sk1kEyIyrlshmiHw1Hur35FUcq9rmj5xesBztfTpk1zvnfffffx9ttvU1BQgL+/f7n1k5OTq/Vc3FHqv/T4xyj8efjPv960Wmly4iine/VDOXFSvXDXwFzaHICT15j7wu0CS83XtS93crXnaOoFaQv12JvkEdCl8paZbnAfTjzyMgXznwUv0/XEdGvdunVTO0Kd57JCFh4ezs6dO52vs7KyCA0NLbc8KyvL+To7O5vQ0FCCg4MpKCjAbrej0+mc2zkcDj766CMmTJiATvfXtPEXfn1eXf/BOptaym+7j3DjlxHovf9qdNv/2EVRcAMat3a/FtlJU9kvyiZNm1z9tidOltvOZPK65n25k4vPu6qCx1k48N882g+IwDukkl8RUVGcPZVDve+3EvLGP68jqRDXx2WXFvv27cuWLVvIzc2lpKSE1atXExMT41weERGByWRytp4SExOJiYnBYDDQvXt3VqxYAUBCQgIxMTFotVrWrFlDUlKS8/3OnTvj4+PjqlNwWzunZ9JyhH+5IgbgSN5GaaOmKqUS7sK/hZHwPj4kz8xCcVTeJT/wkTiKfviF4l92VrqeEK7kskIWFhbGpEmTGDduHCNGjGDo0KF06tSJ8ePHk5JSdl191qxZzJw5kzvuuIPi4mLGjRsHwCuvvMLSpUsZPHgwO3fu5OmnnwbgjTfe4PPPP2fIkCF8++23TJ8+3VXx3Vb+EQvHfyygxcjyl1sVRcG+cyulTZqrE0y4lca3+WLJt19xPEatnw+Bj91D5hMzsOfLWIxCHRrFw+ZnSE5OrtOXFtc+cgrFBm0fCir3vuPkcSzTXuDMyNE0aeZ+rbLNc8pa3n2fvvLYgBe7+BJb6JwnAMh8em71hKulrvXS4nnmLBu738/lpncbEnjBA/WXk7fgG7QmI2EfvXLNxxPiWsnIHh4k/4iFo8vO0eou/0uW2XduQ9umLUi3e1FFXiF6mg31Z+e0zCtOxhkwNhbz9hQKv/+5htIJ8RcpZB5kx7RMWgwLwBhwaQcYx7aNaNq0VSGVcGdhPbwwBmrZt+BspetpvUwETRxD1rPvYDuTXUPphCgjhcxDnE0t5djyAlrGX9oaU3KyUNJPo2nRSoVkwq1pNLSKD+Dk2kKydpVUuqqxTXN8bu9D5j+mozgqb8EJUZ2kkHmIrS9mEDkqAKPfpa0x+/bNaNu0RXOZRxWEuBKDn47IUYHs/HcWlnP2Stf1jx+ILTOX/E++q6F0Qkgh8wiZO0tI31hEixGXtsYA7Js3oIluV8OphCep19ZE/Q4mfnu78lHyNXod9SaO4eyb/w9L6tEaTCjqMilkbk5RFDZOSidqXBB6r0u/nUpOFsqJY2hatVEhnfAkzYb4ce6QhROrKu9mr28Uiv+YIZwZ/ypK6ZUn7RTiekkhc3PHfiig6LSVpnf4XXa5bdMGtNEd0BgMNZxMeBqtQUub+wPZ8+FZCk9WXqB8buuDrl4AOdM/qqF0oi6TQubG7BYHG59Kp/2j9dDqLt+t3vHrOjTtO9VwMuGpfBoaaDLQjx2vZWG3VtyhQ6PREPT3eyj8dg3F63fUYEJRF0khc2N/vJuDTyMDYT0vP0yX4+RxlNwc6a0oqlXDG73R+WrY99EVuuQH+BL0xGgy/zEde3bl6wpxPaSQuanCNCu7ZmbR4e/1KlzHvm412s5d0Wjl2yyqkUZD61EBnFpXxJktRZWuauochXe/bmQ8Ll3yhevIbzg3tfHpdJrHBuDX+PL3vhSbDfuGdWg7193huoTr6P10tLk/kF1vZFOcYa10Xf/7BmNPzyZv/pIaSifqGilkbuj4ygLObCmm9X0BFa7jSN6GJjgYTYOQCtcR4noEtDTSKMaH7VMzK79fptcRNGkseXO+xPzb/hpMKOoKKWRuxlJoZ/2jaXR+uv5lu9ufZ1v5PdquPWswmaiLIm7xReulIWVubqXr6cPqEzg+noyHX5ZR8kW1k0LmZjb96wzBHb0I7e5d4TqO06dQjh9F065jDSYTdZJGQ5t7A8nYWszxKzxf5t23C6Yu0WT+Y0alD1ULcbWkkLmRkz8VciyxgA6PVdzBA8D2wzK0XXug0btsAnAhnHTeWqIfDGLPB7nkHjBXum7AA8OxHj9N3txFNZRO1AVSyNyEOcfGTw+coss/6192PMXzlHP5ODauR9ujbw2mE3WdT0MDkaMC2DYlk5JsW4XraQx66j3zAHlzv6Jk464aTCg8mRQyN6AoCmsfSqPhTT6E9qj4kiKAbeVytG07oPG//LiLQrhKcEcvwnt5s/X5DGwlFXf+0IcEU++p+8kY/yrWUxk1mFB4KilkbuCPOTnkH7bQ7pHKLykqRUXYV32P9sb+NZRMiPIa3+aLqb6OHdMycdgrvg9m6hyFb2x/ztz/HI7iyi9HCnElUshqudMbi0iekUX3KQ3QGSuf3dm2IhFtZBSa+g1qKJ0QF9FoiBwVQGmend/fqXykfN9ht6ALq0/m46/Lw9Liukghq8UKTlpYFX+CG56tj2/Dygf9VfLzsP+4DG2/W2oonRCXp9FriB4XRO6e0kpnlj4/HqP1aBq5Mz6uwYTC00ghq6Us5+wsv+M4reICKxxL8UK2bxah7dBZWmOiVtB5a2k3Poi09YWkflFJMTMaqPfcwxR8s4b8/ybWYELhSaSQ1UL2UgcrRp4gqI2RVqOu3GnDceIY9l9/RttvQA2kE6JqDH46OjwazLEfCji4sOJipgv0p/7Lj5I78/8oWrGhBhMKTyGFrJZx2BRW33cSxaHQ8YlgNJrK74spDgfWj99H1/82NH6Xn5NMCLUYAnV0eDyYo8sL2PdJboX3zPQNQwh+4W9kPv0GxRt21nBK4e6kkNUiDpvCmtEnKUq30e2FBhXOMXYh++ofobgITfdeNZBQiKtnDNTR8fFgTv9SxO45ORX2ZjRGNqXevx4k42+vULL595oNKdyaFLJawmZ2sDLuBAUnrfR8LQSd8crfGkd6GrbFn6MbFi9TtYhazRCgo8NjweQeKGXblAysJfbLrmdqH0m9SeM48+BL0jITVSa//WoBc66N728/hrXIQa9poVUqYorFgnXW9LJLiiGhNZBSiOuj89bSfnzZs5AbHk+n8PTlp38xdWrjbJkV/ij3zMSVSSFTWe4+M1/3PIxvYz3dX2qA1nDly4mKomD9cA6awCA0PXrXQEohqodGryHy7gAadPViw+OnOb3x8hNzmtpHEvzSBLL+9RZ5C76p4ZTC3UghU1Hql3l81+8oreID6fBoMBrtlYsYgG3JFyjHDqONveuKnUGEqHU0Ghr18yX6wSB2v5vDrreysBZfeqnRGNmUBjOeIv+jr8n855sopRYVwgp3IIVMBeYcG0l3n2Dbyxn0fSOMZndUvbehNeFr7L+sRXfPODRGowtTCuFa/s2NdPlnfcw5dtY9lMaZbcWXrKMPq0+DmU9hPXySU4Mfw3rstApJRW0nhawGKQ6F/Z+dZWHbP1E00P/DhgRGVq0YKYqCddFn2Ff/gH7sIzIosPAIOm8tre8NpMWIAP54O5utL56h8FT5lpfWx5t6kx/Cq0cHTt0+nvwvlst8ZqIcmbCqBiiKwsk1hWx5LgOHTaH366EERZmqvn1xEdZ576Ckp6F/8FE0vvK8mPAs9dqaCIysz+kNxfzy2GkaxvgSdX+Qc2g2jUaDX+zNmDq2Ie+DxRQuTSLkrWcwRrdQObmoDaSQuZDDpnAk4Ry73sjCnGMnemwQjW72uar7Wva9u7HOfRtti5boxo1HY6h8zEUh3JXWoKXxrX6E9/Im7Zdifp6QRlg3b1qNCqReOxMajQZD80Y0mPkUxas2kTbsCXzv7Ee9fz2IoUm42vGFiqSQuUDuPjOpX+Zx4LM8fML0tBjhT6ObfNBU4QHn8xzpp7Et+hRH6j50d8SijWrnwsRC1B56Px3NhvjTeIAvGdtK2DEtE51JQ9M7/IkY4ItvuAHfITF49+9O4bK1nLr5IXwG9iHw0bvx6hKtdnyhApcWsuXLl/Phhx9is9l44IEHGDNmTLnl+/fv56WXXqKoqIju3bvz2muvodfrOX36NJMnTyYnJ4cWLVowa9YsfH19OXfuHP/61784efIkwcHBzJkzh5CQEFeeQpXYShykby7mxKoCjn1fQGm+nYhbfOk9I5SAFlXvkKHY7Tj27ca+6gcce3ej69kX/WOTpFOHqJN03loa3exLoxgfzh2xkvVbCX9+lY9XiI7wPj6Edvcm+O4h+I28laLVmzgz7gW0Af743z0Iv6H9MbRsrPYpiBriskKWkZHB7Nmz+e677zAajdx777306tWLyMhI5zqTJ09m+vTpdOnShRdffJGlS5cyevRoXnvtNUaPHs2QIUOYN28eH3zwAZMnT2bOnDl0796djz/+mISEBGbMmMGcOXNcdQqXZTln52xqKTkppWTvKuHM1mLO7i8lMNJIgy7edH66PkFRxip3pVfO5uI4sBf778k4kreh8fNH07kr+onPojFV/T6aEB5LqyEg0khApJGWcQoFx62c3V/K7vdzKEqz4t/USL22XQi6swe+pFG8OYW8Dxaj9fPBO6Yb3r07YeraDkOLCDQ6ndpnI1zAZYVs8+bN9O7dm6CgIAAGDRrEqlWreOKJJwBIS0vDbDbTpUsXAOLi4njvvfcYNWoUO3bsYN68ec7377//fiZPnsz69etZuHAhAEOHDmXatGlYrVYM13nfyFbioOi0Fcs5B5Z8O6V5Dsw5Nkoy7RSdtlKYZqXwhJWC41asRQ78mxkIaG7Ev7mBqPuDCIoyoveuuAOoUlSIcuIYSt5ZlJxslIx0HGknUU4eB0spmqbN0TZtgX7s39A0UL+FKURtpdFpCGhpJKBl2VUKu8VB0UkbhacspG82U5wRSHFGLzTantSz5hKYlIb3qkSMJfPRmgsgtCHaiAj0TSMwNAnFJzoUY8NgtPUC0AX5ow0KQOstf0C6G5cVsszMzHKX/UJDQ9m9e3eFy0NCQsjIyODs2bP4+fmh1+vLvX/xNnq9Hj8/P3JzcwkLCyt37OTk5GsP7v+/f03KXvr+71/5I/w1NbsNyK7KPsMNQOj//ql3v8urdyRZqh392rVeXPbAbBZX/0vm4nPOWrygmlLVbu76vb4WBiDwf//+0gy4gVKgtNKt7VB8tuzfNT6m1q1bt2vbUFQLlxUyh8NRrneeoijlXle0/OL1gAp7+SmKgvaiwXLlB0oIIeoWlz0QHR4eTlbWX38PZmVlERoaWuHy7OxsQkNDCQ4OpqCgALvdfsl2oaGhZGeXtX9sNhtFRUXOS5dCCCHqJpcVsr59+7JlyxZyc3MpKSlh9erVxMTEOJdHRERgMpmclwETExOJiYnBYDDQvXt3VqxYAUBCQoJzu/79+5OQkADAihUr6N69+3XfHxNCCOHeNIoLx3pZvnw5H330EVarlfj4eMaPH8/48eOZOHEiHTt25MCBA0yZMoXCwkLat2/PzJkzMRqNpKWl8fzzz5OTk0PDhg155513CAwMJC8vj+eff56TJ0/i7+/PrFmzaNxYuthWVWFhIffeey/z58/36M+tovP88ssvSUpK4osvvlAxnWvMnTuXlStXAmV/8PXq1Yt33nnHuTwjI4POnTvz0UcfqRVRCNdRRJ3w+++/K0OHDlXat2+vnDx5Uu04LlPRef75559Kv379lPvvv1/FdK6xadMm5Z577lFKS0sVi8WijBs3Tlm9erVzeWZmpnLrrbcqR48eVS+kEC4kgwbXEUuXLuWVV14pd5/SE13uPC0WC1OnTmXixIkqJnOdkJAQnn/+eYxGIwaDgVatWnH69F/d7958803uvfdemjdvrl5IIVxIhqiqI2bMmKF2hBpxufN8++23ueuuuzz2cmrr1q2dXx87doyVK1fy1VdfOV9v3769znz/Rd0kLTLh0TZt2kR6ejp33XWX2lFc7s8//+Thhx/m2Wefdba+lixZwujRozHKMGfCg0khEx7thx9+4M8//2T48OFMmTKFPXv28PTTT6sdq9olJyfz4IMP8swzzzBy5Ejn+2vXrmXw4MEqJhPC9eTSovBoM2fOdH69bds25s6dW+Pjc7paeno6//jHP5g9ezZ9+vRxvp+bm4vZbKZJkyYqphPC9aSQCeHmPvnkE0pLS/nPf/7jfO/ee++lffv2hIfLPF3C87n0OTIhhBDC1eQemRBCCLcmhUwIIYRbk0ImhBDCrUkhE0II4dakkAkhhHBrUsiq2apVqxg7dqxqx//uu++4+eabeeSRR9i9ezdTp05VLcu1OnnyJE8++aTaMa7b+Qewq8u2bdvo1KkTw4cPZ8SIEQwfPpy4uDjWrVvnXOfcuXNMnz6d2NhY53pff/21c/nzzz9Pv379GD58OMOHDyc2NpZbb72VBQtcP2v2V199xccffwzA119/zcKFCy95vyLjx4/n0KFDADz88MPk5ua6NqxwK/IcmYdJSEhg0qRJDB8+nO+++46MjAy1I12106dPc/ToUbVjXLfNmzdzzz33VOs+mzZtSmJiovP1gQMHuO+++1i7di2+vr7cf//9xMbGsmzZMvR6PWlpaTz44IMAjBo1CoAHH3yQRx55xLmP06dPM3jwYAYMGECrVq2qNe+F7rvvPufXycnJzjEiL3y/IhcW2k2bNlV/OOHWpJBVg3fffZfly5cTFBREs2bNnO8fPXqUadOmUVRURFZWFtHR0cyZM4ekpCQWLVrE4sWLgbJfJHfffTfr1q1j/vz5rFmzBoPBQL169Zg5c+YlI9ZXtN+3336blJQUTp06xalTp/j6668pKCjghRdeYObMmaxbt44PP/wQq9WKl5cXzz33HDfccAPvv/8+v//+O5mZmURFRTFr1iznsWw2G6+//jq7du3CYDDQuHFjZs6cyeeff86hQ4d4++23Adi5cyfTp09n7ty5PPDAA9x4443s2bMHu93OxIkTWbJkCUeOHKFDhw688847nD59mrFjx9KvXz/++OMPFEVh6tSp3HDDDUyZMoWMjAweeeQRPvnkE3766Sfmzp2Lw+HA19eXF154gU6dOvH+++9z4sQJMjIyyMrKon379vTq1YuEhAROnTrF5MmTGTp06CXfr8r2l5aWRlZWFmlpaYSFhfHWW28RGhrq/H4ZDAZMJhPTpk0jMjKSjIwMpk2bRnp6OlarlSFDhvD3v/+d2bNnk5mZyb/+9S/efPNNOnfuXC7DvHnz+PHHH9HpdLRo0YKXX36ZkJAQxo4dS2BgIEeOHOG+++67Yus+OjoaLy8v0tLSOHToED4+PowfP965PCIigjlz5mC1Wivcx5kzZ1AUBT8/v0uWjR07lnbt2pGcnMzZs2cZPny4cxaBij7Hw4cP89JLL2GxWFAUhfj4eMaMGcP777/P2bNn6dOnD+vWrWPTpk14eXmRm5vL2bNnGTBgAG+88QbLly8HylqXt956Kz/99BMjR47k3XffZdGiRQA88MADvPzyyzz77LOsW7cOrVZLSUkJAwYM4McffyQ4OLjSz014GFUnkfEAa9asUQYPHqwUFBQoVqtVmTBhgnPOq//85z9KQkKCoiiKYrFYlKFDhyqrVq1SSktLlT59+igHDx5UFEVR5syZo8yaNUs5ffq00rVrV6W0tFRRFEX55JNPlDVr1lxyzIr2qyiKcv/99ysrV65UFEVRvv32W2XChAmKoijK0aNHlaFDhyq5ubmKoijKwYMHlRtvvFEpKipS3nvvPWXQoEGK1Wq95Fg7duxQ7rjjDsXhcCiKoihvvvmmkpycrGRnZytdu3ZVzp49qyiKokyePFn56quvlJMnTypt2rRRfvrpJ0VRFGXq1KnKLbfcohQUFChms1m58cYbleTkZOd633//vaIoirJ+/XrlxhtvVCwWi7J161ZlyJAhiqIoyqFDh5S+ffsqJ06cUBRFUTZv3qzceOONSkFBgfLee+8pt9xyi3Lu3DmlpKRE6dGjhzJz5kzn92XgwIGXnM+V9nfrrbcqBQUFiqIoyqOPPqq8++67is1mU9q3b69kZGQoiqIoy5YtUxYvXqwoiqKMHTtWWbt2raIoimI2m5WxY8cqP/74o6IoinLLLbcou3fvviTDN998o9xzzz1KUVGRoiiK8t577ykPP/yw8/v3wgsvXLKNoijlPpfzkpKSlL59+yrFxcXKtGnTlDfeeOOy25733HPPKTfddJMybNgwZcCAAUrPnj2Vxx57TNmyZctl17///vuV8ePHKxaLRcnPz1cGDRqkrFu3rtLP8YUXXlA++ugjRVHK5kJ7+umnFbvdrrz33nvKa6+95szxf//3f87zf+211xSHw1HuM1u4cKHyzDPPXPJZtmnTRsnJyVEURVGGDRumrF+/XlEURfn666+VSZMmVXr+wjNJi+w6bdmyhdtvv9351+xdd93lnIF48uTJbNq0iQULFnDs2DEyMzMpLi7GaDQyatQovv76a5577jmWLVvGF198QVhYGNHR0YwcOZKYmBhiYmLKjZ13XkX7rcymTZvIzMx0XmYC0Gg0nDhxAoAuXbqg11/649CmTRt0Oh2jRo3ipptuYtCgQXTq1AmAm2++mcTEREaMGMHGjRt55ZVXOHv2LAaDgQEDBgBll8JuuOEG5+cTGhpKfn4+oaGhBAYGEhsbC5TNaqzT6UhNTS13/K1bt9K7d2/neIF9+vQhODjYee+pb9+++Pv7O/fdr18/53Hz8vIuOZ8r7a9nz57OrO3atSM/Px+dTscdd9zBvffey80338xNN91E//79KS4uZseOHeTn5/Puu+8CUFxczIEDByodqHfDhg3ExcXh4+MDwLhx45g/fz4WiwWA7t27V7jtiRMnGD58OFDWWg4PD+eDDz7A29sbjUaDUoWBes5fWiwuLmbSpEkYjUZ69epV4fr33HMPBoMBg8HAHXfcwcaNG2nZsmWFn+Ptt9/Oc889x+7du+nTpw9TpkxBq73y7XiNRsNdd93FsmXL6NixI9999x3PPvtspduMGTOGpUuX0r9/f5YsWXLF9YVnkkJWDS785aHT6Zxf//Of/8Rut3PnnXdy8803k56e7lz33nvvJT4+np49e9K6dWvnL4Qvv/ySlJQUtmzZwr///W/69et3yX/OyvZbEYfDQZ8+fcoNmJuenk5oaChr1qxx/lK9WEBAAImJiezatYutW7fy9NNP88gjjzBmzBjGjBnDq6++il6vZ+DAgfj6+joLmUajce7DYDBcdt8XflbnM17uvQv3BWWft81mA7hkepLLFeOr2Z+Xl5fz/QsLw6xZszh48CCbN2/m448/JjExkRkzZqAoCosXL8bb2xsoG6jXZDJdVQaHw+E8PlDh9wIuvUd2oS5dujg7UFxo7dq17Ny5k+eee67c+z4+Prz55psMHjyYzz77jIceeuiy+73wM1UUBa1WW+nneMstt5CUlMTmzZvZsmUL8+bN47vvvqvwnC4UHx/PyJEjGTVqFAUFBfTs2bPS9WNjY3nnnXfYunUrxcXF9OjRo0rHEZ5Fei1ep5iYGFatWsW5c+dwOBzlfsls3LiRf/zjH86/zv/44w/sdjsADRs2pEuXLvz73/923uw+cOAAQ4cOpVWrVjz66KM8+OCDpKSkXHLMyvZ7IZ1O5/wF2adPHzZt2sThw4cB+OWXXxg2bBhms7nS8/v555958MEHueGGG3jyyScZMWKEs/XStWtXtFotn3zyCffee+9VfW5Q9kt/w4YNAKxbtw6DweBsAZ6/p9OnTx82btzIyZMngbIWcHp6+iX3nKrqWvaXm5tL//79CQoK4sEHH+Tpp58mJSUFPz8/unTpwqeffgqU3dM53/ECyn/+F+rXrx/ffvutsxX9xRdf0KNHj+ueM2zgwIEUFhayYMEC58/DyZMn+c9//lNhJ47AwECee+453nvvvQo7Bn3//fc4HA7y8/NZuXIlAwYMqPRzfOaZZ1ixYgVDhgzhlVdewc/Pz9nyP6+izyYsLIxOnToxdepU4uPjL5vnwm29vb0ZNmwYL7744jX9DArPIC2y69S/f39SU1O56667CAgIIDo6mrNnzwIwadIk/vGPf+Dj44Ofnx89evQo9x86Li6O119/nf79+wNlN+7vvPNO7rrrLnx8fPDy8mLKlCmXHPNK+z2vS5cuzJs3jyeeeIK5c+cybdo0/vnPf6IoCnq9ng8//BBfX99Kzy8mJoYNGzYwdOhQfHx8CAwM5PXXXy93DitWrCA6OvqqPzuTyURiYiKzZs3Cy8uLefPmodPpiIyMxGQyER8fz9dff80rr7zCE088gd1ux8vLi/nz5zsvJ16tyMjIq95fcHAwjz32GA8++CBeXl7odDqmT58OlLXUXn/9dWJjY7FYLAwdOpRhw4YBcPvttzN58mReffVVbrrpJuf+4uPjSU9PZ9SoUTgcDpo1a1aug821MhqNfPrpp7z11lvExsai0+nQ6XQ89thjxMXFVbjdsGHD+Prrr3njjTd45513LlluNpuJj4+nqKiI0aNHOy93V/Q5Pv7447z00kssWbIEnU7HbbfdRo8ePdi2bZtznzExMeVG67/QqFGjeOqpp/jwww8vu/yOO+5g7NixvP/++7Rp04a4uDiWLl3KiBEjruLTEp5ERr9XicPhYNq0aTRq1IgJEyaoHeea2Gw2nnjiCYYNG3bVkzeeOnWK2NhYfvvtNxelE9Vh7NixjBkzhjvuuEPtKJelKAoLFiwgLS2N1157Te04QiVyaVEFhYWF9OrVi/T0dMaNG6d2nGty6NAh+vTpQ7169WrtLznh+W699VbWrVvHU089pXYUoSJpkQkhhHBr0iITQgjh1qSQCSGEcGtSyIQQQrg1KWRCCCHcmhQyIYQQbk0KmRBCCLf2/wEZTZ8bXc3o9AAAAABJRU5ErkJggg==\n",
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sns.set(rc = {'figure.figsize':(8,8)})\n",
"plt.style.use('seaborn-whitegrid')\n",
"sns.displot(df_1, x=\"days\", hue=\"variant\",kind=\"kde\",fill=True, palette=[\"#DC143C\",\"darkviolet\"])\n",
"plt.axvline(delta_med, c='red')\n",
"plt.axvline(om_med, c='darkviolet')\n",
"plt.xticks(np.arange(1, 32, step=13))\n",
"plt.xlabel('days after symptom onset or PCR positivity')\n",
"plt.savefig('KDE_days.pdf', bbox_inches='tight',dpi=1200)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"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"
}
},
"nbformat": 4,
"nbformat_minor": 4
}