Skip to content

Instantly share code, notes, and snippets.

@jasonreich
Last active June 7, 2021 16:41
Show Gist options
  • Save jasonreich/f132aa1703ecded98967aadbed26efa4 to your computer and use it in GitHub Desktop.
Save jasonreich/f132aa1703ecded98967aadbed26efa4 to your computer and use it in GitHub Desktop.
England and York
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "import pandas as pd\nimport numpy as np\nfrom uk_covid19 import Cov19API\n%matplotlib inline",
"execution_count": 1,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "england_only = [\n \"areaType=nation\",\n \"areaName=England\"\n]\nyork_only = [\n \"areaType=utla\",\n \"areaName=York\"\n]\ncases = {\n \"date\": \"date\",\n \"newCases\": \"newCasesBySpecimenDate\"\n}\napi_e = Cov19API(filters=england_only, structure=cases)\ndf_e = api_e.get_dataframe()\napi_y = Cov19API(filters=york_only, structure=cases)\ndf_y = api_y.get_dataframe()",
"execution_count": 2,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "df_pop = pd.read_excel(\n 'ukmidyearestimates20192020ladcodes.xls',\n sheet_name='MYE2 - Persons',\n header=4\n)\ndef get_population(name):\n return df_pop[df_pop['Name'].str.lower() == name.lower()]['All ages'].values[0]",
"execution_count": 3,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "def index_and_window(df, name):\n population = get_population(name)\n df['date'] = pd.to_datetime(df['date'])\n df = df.set_index('date')\n df = df[::-1]\n df[name] = df['newCases'].rolling(7).sum() * 100000 / population\n df = df[-93:-2].drop('newCases', axis=1)\n # df = df.rename(columns = {\"newCases\": name})\n return df\n\ndf_e = index_and_window(df_e, 'England')\ndf_y = index_and_window(df_y, 'York')\n\ndf = df_e.join(df_y)",
"execution_count": 4,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "df.plot.line(logy=True, ylabel=\"7-day cases per 100,000\", style='.-')",
"execution_count": 5,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 5,
"data": {
"text/plain": "<AxesSubplot:xlabel='date', ylabel='7-day cases per 100,000'>"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 432x288 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEpCAYAAABbU781AAAASXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMisyNTQ1LmdlZjY3M2M0NmUsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvfRvX6wAAAAlwSFlzAAALEwAACxMBAJqcGAAAXfZJREFUeJztnXd4W+XVwH9HtmNn7x3Hzg4ZkEUIIwTCSikEaMveswMK3YXSr4wuOmhLB7TMAA1QNoS9SQKZTkL2dOzYGc5ydhwPne+P98qWFVmWbcnyOL/n0XN1h+49ule6575niqpiGIZhGLHAl2gBDMMwjKaDKRXDMAwjZphSMQzDMGKGKRXDMAwjZphSMQzDMGKGKRXDMAwjZiQnWoBE06VLF83MzEy0GIZhGI2KrKysnaraNXR5s1cqmZmZLFy4MNFiGIZhNCpEJDfccjN/GYZhGDHDlIphGIYRM0ypGIZhGDGj2ftUDMNoXpSUlJCfn09RUVGiRWkUpKWl0adPH1JSUqLavkkqFRHpD9wNtFfVbyVaHsMwGg75+fm0bduWzMxMRCTR4jRoVJVdu3aRn59Pv379ovpMXM1fItJBRF4WkdUiskpETqzlfp4Uke0isjzMuikiskZE1ovInQCqmq2qN0az7+37j5CVW1gbsQzDaIQUFRXRuXNnUyhRICJ07ty5RqO6ePtUHgLeU9WhwHHAquCVItJNRNqGLBsYZj/TgCmhC0UkCfgX8DVgGHC5iAyriYAF+4q48vG5plgMoxlhCiV6anqu4qZURKQdcCrwBICqFqvqnpDNJgFviEia95mbgb+H7ktVZwK7wxxmPLDeG5kUAy8AF0Qp3/ki8ihAUYmfRz5bz6drtvOvT9ebgjEMI64kJSUxatSo8tcDDzxQ6321adMmJjLl5OQwYsSIOu8nnj6V/sAO4CkROQ7IAu5Q1YOBDVT1JRHpB7wgIi8BNwBn1eAYvYG8oPl84AQR6Qz8FhgtInep6u9DP6iqM4AZqT0H3SzAR6u289Gq7QAk+4Rvn9qPli2SOXFAF8ZmdKzJ9zYMw4hIy5YtWbJkSaLFiAvxNH8lA2OAR1R1NHAQuDN0I1X9I1AEPAJMVdUDNThGuHGZquouVf2Oqg4Ip1CC6d4ujZe/cyLXnphZvrNSv/Kvz7L58wdrufQ/c3hn6VaycgttFGMYzZT6+v9nZmZyzz33MGbMGEaOHMnq1asB2LFjB2eddRZjxozh29/+NhkZGezcubPSZw8cOMAZZ5xR/tk33ngDcCOQY445hptvvpnhw4dz9tlnc/jwYfe9srI47rjjOPHEE/nXv/4Vk+8Qz5FKPpCvqvO8+ZcJo1REZCIwAngNuAe4rYbHSA+a7wNsqYmQ3dqmMjazE4jwv4WbKCn1gwh+v6I4BfO95xbhE1CF5CTh7nOPoU1aMht2HOT0IV1J8vmYm72LCf07MzajI1k5u5mTvavKUU5WbmGl7Q3DSAz3zVjByi37Im6zv6iE1dv241fwCQzt0Za2aVWH1w7r1Y57zh8ecZ+HDx9m1KhR5fN33XUXl156KQBdunRh0aJFPPzww/z5z3/m8ccf57777mPy5MncddddvPfeezz66KNH7TMtLY3XXnuNdu3asXPnTiZMmMDUqVMBWLduHc8//zyPPfYYl1xyCa+88gpXXXUV119/Pf/4xz+YNGkSP/3pTyPKHC1xUyqquk1E8kRkiKquAc4AVgZvIyKjgceArwMbgf+KyG9U9ZdRHmYBMMgzoW0GLgOuqI28YzM6Mv2mCczN3kXHVi24/60VlJT6SU7yMbh7G5Ztdj+8kjLl3hkVX+ORzzZU2k9qso8jpX4AfLKW84/rxZAebcnffZh+XVtTVFzG3z9ZR2mZkpLs4+fnDKGotIwJ/c3MZhgNkX1FpfjVvferm4+kVKIhkvnrG9/4BgBjx47l1VdfBWD27Nm89tprAEyZMoWOHY++V6gqv/jFL5g5cyY+n4/NmzdTUFAAQL9+/cqV2NixY8nJyWHv3r3s2bOHSZMmAXD11Vfz7rvv1ul7QfzzVL4PTBeRFkA2cH3I+lbAxaq6AUBErgWuC92JiDwPnAZ0EZF84B5VfUJVS0XkNuB9IAl4UlVX1FbYsRkdy2/sQ3q0LR9NAFz5+NxyJTNxUBc+XrUdDbOPzq1bsHVvEYr7Ab65ZEvY7QCKS/38+m0XEJfkW8dDl42iZ/uWMRvF2IjIMCJT3YgC3P8o8P9PSfbx0GWj4/p/Sk1NBZwzv7S0FHAKozqmT5/Ojh07yMrKIiUlhczMzPJQ4MA+A/s9fPgwqhqXKLi4KhVVXQKMi7D+i5D5EtzIJXS7yyPs4x3gndpLGZ5gBQOUj2ICSmb2+p2UlPpJ8gmIUFbmfnC3TR5UPspJSfZxwahevLQwv3zo/LURPfho1XZKytxoJvAEVOZXbntucbmZLSXJx8NXjaFjyxTmbtxdYVqLUlHMWreD659aQJlfaZHs457zhlF4uMQUjGHUkGArRqL+P6eccgovvvgiP//5z/nggw8oLDzat7N37166detGSkoKn376Kbm5YYsIl9OhQwfat2/P7NmzOeWUU5g+fXpMZG2SGfXxoDolE/yDCx3lvLFkS7mSueGU/txwSv+wZrYh3duydPNeAIrL/Nz0dEVJfp/AqYO68MWGXeWKYvpNE8L+wPN2H+KH/1tCqaexjpT6+cXrLm+0RZKP6TefgE+kksyfr9nO/JzdTBrclfH9OtsoxzCCCP3/15VQn8qUKVMihhXfc889XH755fzvf/9j0qRJ9OzZk7ZtK6X4ceWVV3L++eczbtw4Ro0axdChQ6uV46mnnuKGG26gVatWnHPOObX+PsFINMOqpsy4ceM03v1UIt2gg9eBM7MVl/pJ9gnH9ulAVm5hWPOZAOcd15OhPdpVGsW8uiift5ZuoaRMKS1TSv2VR0QAack+SvyK36/4ROjYKoWdB4vL16cm+Sj2RlItkn08d3N45WUYjZFVq1ZxzDHHJFqMGnHkyBGSkpJITk5mzpw5fPe7363XkORw50xEslT1KEuUjVTqgUhPOdWNgIJtub889xjuf3sVxaV+FJjx1VZmfLUVn7j9LNq0hzK/IsBDl42id8dWR42IfD6hc5sWbN7jbK1lqvhxn1GcsurctgVbvPVHSv3c/doyTh/SjTOHdTflYhgJYNOmTVxyySX4/X5atGjBY48d5SVoMJhSaWBEUjJjMzpyTK/2zM3exYYdB3ht0ebygIAlnkIBZyrLKzzM1FG9qww8uOKxuZSU+WmR7OMnZw+t5Ae67XTnFyou9aMKq7ftZ/W2/fxn5gZ+ce4xHJfegflBfh7DMOLLoEGDWLx4caLFiApTKg2cUCUTmM/KLeSdZVvLFcGvzhteSTEElEdV+3nu5glV+oGC57fsOczz8zfhV6e8fvP2qvIkUTONGYYRivlU6sGnEi9CfTXxcK6HhlOOy+jE7PUVmbwDurbh0uP7UFLmt1wbo1HQGH0qicZ8Ks2EqkYxsT5GqJ9nYe5uikv9+ETYtvcwv3tnNQKkpqyvMiLNMIzmgSkVo1oi+Xm+3LCTBz9YiwIlpX7mZu8ypWIYzRjrUW/UmLEZHbn19IGMzejISQO6kJrs/YxEjvLlGIZRGVXllFNOqVQS5cUXX2TKlKNaRh3FtGnTuO22mpRHrH9MqRh1YmxGR567eQIjerVDBHq2T0u0SIbRoBER/v3vf/OjH/2IoqIiDh48yN13311tleBAyZaGjikVo86MzejIf64Zhw/hoY/WJVocw4g9efNh1oNuGgNGjBjB+eefzx/+8Afuu+8+rrrqKn784x9z7LHHMmHCBJYuXQrAvffeyy233MLZZ5/NNddcU2kfb7/9NieeeOJRJfATjflUjJjQu0NLrpzQl6e/zOGWSf0Z0DU23egMI668eydsWxZ5myP7oGA5qB/EB91HQGq7qrfvMRK+Vn0nx0DflBYtWnDKKacwevRoXn/9dT755BOuueaa8oz5rKwsZs+eTcuWLZk2bRoAr732Gn/5y1945513wlYsTiTNdqQSaCe8d+/eRIvSZLj19IGkJPn43n8XWTMzo+lQtNcpFHDTotjcM1q3bs2ll17K1Vdfzdy5c7n66qsBmDx5Mrt27SJwb5o6dSotW7Ys/9ynn37KH/7wB95+++0Gp1CgGY9UAu2Ex40bd3OiZWkq5O46RJlfWVOwn8sfncvzt1h4sdHAiWJEQd58eHoqlBVDUgv45uOQPj4mh/f5fPh8vrCl7QNl6Vu3bl1pef/+/cnOzmbt2rWMG1dlEfiE0WxHKkbsmZu9C7/35ygu8zNz7fYES2QYMSB9PFz7Jky+201jpFCCOfXUU8tLz3/22Wd06dKFdu3Cm9gyMjJ49dVXueaaa1ixotbto+JGsx2pGLFnQv/OtEj2UVzqx6+wLN9Mi0YTIX18XJRJgHvvvZfrr7+eY489llatWvH0009H3H7IkCFMnz6diy++mBkzZjBgwIC4yVZTrExLIy7T0hAJlIrJ3nGAVxZt5t9XjWXKiB6JFsswyrEyLTUn9Jxl5RZywohBm0v37+wTuq2NVIyYEsi+Ly71s6ZgPz99+SuW5u/hjGOsbL5hNAWycgu59D9zSGrdsXe49eZTMeJCi2QfN5/Sn/1FpTz82QaufHyuRYQZRhPgxYV5rqtsFe3tTakYcSN/z+Hy312xVxfMMIzGi6ryVd6eipkwmFIx4saE/p0r6oJ584bREGjuvuSaEHyuPlu7g9Xb9nPLqf0pO7hnS7jtTakYcWNsRkem3zyB8ZkdUYVubVMTLZJhkJaWxq5du0yxRIGqsmvXLtLS0lBVHvxgDemdWvKTs4dQdmDXtnCfMUe9EVfGZnTk75ePYeIfP+GJ2Ru5d+rwRItkNHP69OlDfn4+O3bsSLQojYK0tDT69OnDe8u3sXzzPv588XG0SK56PGJKxYg7PdqnceGo3rywYBO3nzGITq1bJFokoxmTkpJCv379Ei1Go6LMrzz44VoGdG3NRaPDBn2VY+avYGJcidSo4JZT+1NU4ufZObmJFsUwjBry0EdrWb/9ABeO6k2Sr4qwLw8bqQTI+QKeuQC0DJJS41aOobkyqHtbzhjajcdnZwPKKYO6Wt6KYTQCFubs5h+frAfgX5+t56SBXSL+d22kAlB8EN68DfwlrgppaRFsnJloqZocpw3pyv6iUv720TrLWzGMRsJz8zcRCGkoiSI1wEYq/jJ49iLYvRF8yeAvBRQWPgWFOZB5iuuPsG0ZbFkMvUa7+YKVsGcT9JtoI5oo2VdUAmD97A2jEbGuYD8ASQIpyb5qUwNMqexYBZs3wSVPQ9uesHEW7MmBRc/A4mfdKxLJLc1UFiUT+nehRdJ6isv8iPWzN4wGz+pt+1i2eR9XnZBBzw5pTOjfudoHQTN/lZW4bm5tezrFcOqPoWMmSJJbH+j0FjhV5fOes6q0CHJmJUDwxsfYjI48f8sE+nVpTWqyj4HdrDukYTRknpy9kbQUHz85ZzC3nj4wKsuCKRVwJrBgxZA50TXjkSTntD/+ZkhODZlP8zZWaNkpIWI3RsZmdOSfV4zmYHEZj87ckGhxDMOogp0HjvD6ki18c0wfOrSKPg3AzF+IUyCZEysWBZry5Mxyy9PHQ/dhR8+v+9CZxz6+3/lXhnzNzGBRMLxXey4Y1YsnZ+dw7YmZdGuXVv2HDMOoV56bt4niUj/Xn5xZo89F1U9FRDoBqqqNIlxHRPoDdwPtVfVbkbYdNzRdF374Su2VwZIX4PVvu/fJaXDtDFMsUZC76yBnPPg5Zx7TnZF92kdlqzUMo36Ym72TG6ctZFD3trx+68lhtxGRLFU9qp9xleYvEekrIi+IyA5gHrBARLZ7yzKjFU5EkkRksYi8Fe1nwuzjSe/Yy8OsmyIia0RkvYjcCaCq2ap6Y1Q7b9O9bkpg/2Yq/CtHzL8SJRmdW3PmMd15b8U2HvxgjYUYG0YDISu3kKufmM/B4jJWbtlX4/9lJJ/K/4DXgB6qOkhVBwI9gdeBF2pwjDuAVeFWiEg3EWkbsmxgmE2nAVPCfD4J+BfwNWAYcLmIDKuBbHUnc6LnXxFAnc/FiIr+XVsD4Nfo4t8Nw4g/c7N3UlLmLFhl/pr/LyMplS6q+j9VLQssUNUyVX0BiCoWVET6AF8HHq9ik0nAGyKS5m1/M/D30I1UdSawO8znxwPrvZFJMU7ZXRClbOeLyKN799axj3rA/3LaXdCxvyvzsndz3fbZTDjjmO4keyUfkpOqj383DCP+tE9zTnkhuryUUCIplSwReVhEThCRXt7rBBF5GFgc5f7/BvwM8IdbqaovAe8BL4jIlcANwCXRi09vIC9oPh/oLSKdReTfwGgRuauKY89Q1Vvat29fg8NVQfp4OO3ncOWLzgT23CXw+Z+shlg1jM3oyJPXHk9qso+RvdubT8UwGgCz1++kdYskfnDmIKbfNKHG/8tISuUaYBlwH/A+8IH3fjlwdXU7FpHzgO2qmhVpO1X9I1AEPAJMVdUD0YnuDhN+l7pLVb+jqgNU9fc12F/d6DIITvo+FCyHT38LT081xVINpw7pyu1nDGJhbiFZueEGo4Zh1Bd5uw/xwcptXHNSJnecObhWD3pVKhVVLVbVR1R1iqqOVNUR3vuHVfVIFPs+GZgqIjk4s9RkEflv6EYiMhEYgfPf3FND+fOB9KD5PkDYbmT1RnIgnluhrNgc91Fw/cmZdGmTyh/fW2ONkwwjgTwzJwcR4eoJGbXeR8TkRxE5R0QeEZE3ReQN7/1RDvNwqOpdqtpHVTOBy4BPVPWqkP2PBh7D+UGuBzqJyG9qIP8CYJCI9BORFt5x3qzB52NP5kSXKAkgUjn/xQhLqxbJfH/yQOZt3M3PXl5qUWCGkQAOHinlhQV5fG1ED3p1aFnr/UQKKf4bLnLrc+CPwJ+897eLyEO1PmJlWgEXq+oGVfUD1wJHNdwQkeeBOcAQEckXkRsBVLUUuA1nnlsFvKiqK2IkW+1IHw/XvgUdMqB1F+h9VBi3EYahPdoiwEtZ+RZebBgJ4JVF+ewvKuX6k+vWwCxSRv25qjo4dKGI/A9Yi1M4UaGqnwGfhVn+Rch8CW7kErrd5RH2/Q7wTrSy1Avp42Hy/8GrN8HGz2HA6YmWqMGzMEiJHCmxCsaGUZ/4/cojn22gR/s0qKMJOpL5q0hEwmUFHo9zrBuROOZ8VxMsa1qiJWkUTOjfmdQUXyDbhzK/+VYMo754fHY2W/cWUbC3iCufmFcnS0EkpXId8A8RWSkiH3ivVcA/vHVGJFLSYNQVsPotOLA90dI0eMZmdGT6TRP44VmDGdqjLQ9/tp6VW/YlWizDaBb8d+4moHKvo9oSKfprkaqeAEwG7gJ+AZyuqidUFyZseIy51jX9WjI90ZI0CsZmdOT2MwbxzI3j6dCyBdc+NY8/v7/G/CuGEUdWbtnHpt2HSPZJ1I24IlFd9JcAGd4rHcjwlhnR0HUwZJwMWU+DP2z+pxGGbm3T+MFZg9ixv5h/frreHPeGEUeemL2RlilJPHnt8fzo7CG1SngMJlL019nAOuBe4FxcuZX7gHXeOiMaxl4HhRthxu2WCFkDdh0oLs9stbpghhEftu8r4s2vNnPJuD6cOqRr1I24IhEp+ush4ExVzQleKCL9cNFWx9TpyM2Ftj3ddPGzsOxlaz0cJRP6dyYl2UdxqR+fz1oPG0Y8eHZuLqV+rXMYcTCRzF/JuIz1UDYDKTGToKmTPx8rjV9zxmZ05LmbTqB1iyTGZXS08GLDiDGHi8v479xczjymO5ldWsdsv5FGKk/ieqi8QEXRxnRc1voTMZOgqRMojV96GPB7/e2NaBiX2YkLRvfm9cWbKSopIy0lKdEiGUaT4W8fraXwUAkTB3WJ6X4jRX/9HrgC95h9InCS9/7Kei3S2NgJlMaf8D3wJcOS5xItUaPinOE9OFRcxux1OxMtimE0GV5btJlHZ2YD8Lt3VsU0ECZij3pVXUUVDbaMGpA+3r1adYZPfg2z/gJoRb97o0pO7N+ZtmnJvL9iG2cO655ocQyj0ZKVW8isdTvYuPMgM77aQiC9OBAIEysTc0SlEg4ReRo4BPxLVY9q72tE4OQ7XM7Kx/eB+FyXSHPcR6RFso/JQ7vx0aoCSsv8JCdFjII3DCMMWbmFXP7oXIrLXGrD2IwOLN+8j9Iyf53zUkKpzT/0n8BHRNFTxQghKQUGTHbv1W+l8aNkyvAeFB4qYUGO5aoYRm14YcGmcoXiE5g8tDvP3TwhJnkpodR4pKKqC3Al51+JmRTNiWMvhQVP4PrZt7DS+FEwaUhXUpN9vL9iGycOsNBiw6gJc7N38cbizQhOoQRGJmPjFFVZpVIRkWTgRuAioBeuLMwW4A3gCa+isFFT0sdD/9NdqPFVr5rpKwpatUhm4qCufLBiG/ecPwwr6mAY1ZOVW8iri/J5JSufvp1b84tzh7Jq6/5yhRIvIo1UngX24DLqA/kqfXA9T/4LXBo3qZo6x3wdsj+BNt0SLUmj4Zzh3floVQG/emMFF47ubXkrhhGBrNxCrnhsLkdK/Qhw55ShTB7anclD4x/sEsmnMkZVv6uqc1U133vNVdXvAqPjLllTpu+JbrppbmLlaER0beu6aT47N9dqgRlGNcxcu50jpc6HIgJrCvbX27EjKZVCEblYRMq3ERGfiFwK2D+6LnQ9BlLbw6Y5iZak0bBiy77yWmDFVgvMMKpEVVnoBbX4xEVQ1meZo0jmr8uAPwAPi0ghLvGxPfCpt86oLT4f9D3BRio1YEL/zqQm+ygq9aMK48z8ZRhhmT5vE19s2MUVJ/Sld4eWcfehhFKlUvEKSV4KICKdAVFVS2uOFX0nwLoP4OAuaG0RTdUxNqMj02+ewPR5uby6aDNZmwo5wYpMGkYlXlqYxz1vrmBUegd+c8EIfL76D2qJGFIsIu2BKUBvQEVkC/C+qu6pB9maNgG/St48GHpu/I+XNx/WfwwDz2i0EWeBEMhDR8p46KN1nDeyF307t0q0WIbRIPhy/U5+9vJSFFi1dR+L8/YkJKAlUj+Va4BFwGlAK6A1cDqQ5a1rsIhIfxF5QkReTrQsVdJrjMtTqQ+/St58eOpr8PkD8PR50fV1yZsPsx5skD1g7pk6DAGufWo+WTm7Ey2OYTQIHvp4XXnpldKyxPkdI41U7gbGho5KRKQjMA94JtKORSQNmAmkesd5WVXvqY2QIvIkcB6wXVVHhKybguv9kgQ8rqoPqGo2cGODViopadBrdP34VRY949oaA5R6WfyRRisbZ8GzF7ms/6QWDa6UzJY9RZT6lY07D3Lxf+bwyFVjOWd4j0SLZRgJI2fnQbJyC0kSATTmpVdqQiSl4qQ7Gj/lDUIicgSYrKoHRCQFmC0i76pq+V1URLoBh1V1f9Cygaq6PmRf03DlYSopMhFJAv4FnIXLpVkgIm+q6soo5Es8fSfAnIeh5DCktIzPMfZvg1UzqHQ5+54c+TOzHgS/l9taFoUSqmfmZu/Cr+67+BW+/9xirj05k/ZpyZw4oIvlsBjNjvvfWklaShJ/u3QUawrin+AYiUhK5bfAIhH5gIp+Kn1xN/BfV7djVVXggDeb4r1CldQk4Lsicq6qFonIzbgM/kpOBlWdKSKZYQ4zHljvjUzwer9cADQSpXIifPEQbF4EmdXc6GtDWSm8cpNTDN94FNa8CyteBX9x1Z85sAM2zauYT0pucKVkJvTvTItkHyWlrsBknw4tecwr452WvJ7pN8e2lpFhNGQ+WV3AJ6u384tzh3LmsO4Jr+YdqZ/K08A44HPcqKMY+AwYp6rTotm5iCSJyBJgO/Chqs4LXq+qLwHvAS+IyJXADcAlNZC/NxUKD9xopbeIdBaRfwOjReSuKmQ7X0Qe3bt3bw0OF2PST3DTYL/Kuo/go/ti48uY8X03yjjxVjj2ErjwYUhtB1+9UPVnPn8Ayo7AuX8GXwr0Gd+gRingRYLd5IrhPXfzBL4xtnf50Lmo1M+nq7cnVD7DqC/mZu/kRy9+Ra8OaVx3UuxaAlfLnIfp3VbC2pyr66dSCES4A0VGVcuAUSLSAXhNREaElstX1T96I4xHgAGqeiDMrqoinBlOVXUX8J1qZJsBzBg3btzNNThebGnVCTpkwJLnoVN/18N+zdtu3RcPwdS/Q5fBTjEEeq/kza+Yh8rrgnn3zoqGYF/+Ewad7bYZdgGseA3O/RO0CGkhumMtLHwKxt0A42+GQ7vgs9+7YwbvP1iGBCmc0GJ4qSnrKS7141d4Y8lmRqV3SLgZwDDiSVZuIVc9Pp9Sv3LoSBnLNu+tn9/6u3fCvEfo0UZ6h1td4yrFACKyTFVHRru9qu4Rkc9w4cmVlIqITARGAK8B9wC31UCUfFyL4wB9cEUvGwd582FvPmgZvHy967ESQMvgjVvdMlXwJUHmqZAzE/xlFduqH5JT4doZ7gavCp/+FuY9UrGvYL/IcZfD4mdh9dtu9BLMR/dCSis47U43f+JtsPBJ+OCXcMP7rt5Dzmx49kIoK/OOm3gnfmDkMjd7F61Tk/j9O6u5+ZmFiJdNHOvS3obREPh4VQGlfudRKPPHttFWWMLdW8IQqUrxN6paBVQbaiMiXYEST6G0BM7EZegHbzMaeAz4OrAR+K+I/EZVf1nd/j0WAINEpB+wGZfpf0WUn008ObPchQJAXFn8Fa87JZCUAt1HwuYFbrW/FLI/pdwtpWUV+yktgg/+z5m5vvgbbM6CwV+D7M+8fQWV2O97InToC189X1mpLHjCjZLGXg+tvZ7VqW3g9Lthxu3w4jXQuhssfQHKGp4TP3jksq7gANPnbUI1UNJlZ6U/W1ZuIXOzd9koxmjUZO84CFQuZx83cufA+3fDliwYPAWyP0c5EC6QK+JI5X/AdMJHgKVFIUZP4GkvQssHvKiqb4Vs0wq4WFU3AIjItcB1oTsSkedx+TJdRCQfuEdVn1DVUhG5DXgfF1L8pKquiEK2hkHmRPe0H7jxj7vBvYLNW09PrVg/5QF4704370sCxIvSEpdEmecF1vmS4ZQfwsQfHW2m8vng2Mtg1p9h3xZo2xPm/8cNacEpm1FXVGzfZbDb/6o33XznQVC40Sk5X1KDc+IDfGNMH17OyueIZw6buXYHx/Rsx8KcQvILD/HW0q341f0ZLxrdmyE92rLzQDEnD+jMSQO7sDR/L3M27GRE7/ZkdG7Nyi17ydl1yJSQ0WDYsOMAH64qYMrwHozs0z6+v81Fz8CMO5xVxJcMp/wIJv6Ygj9OCGsVEtWwygYRyQKuDdcyWETyVDU9zMcaHePGjdOFCxcmToDq/BOh66vyqax6C778O6AgSTD5bpj44/DH3LUB/jEGug93o46dayvWhX521oPwyW/cD0p8MPmX0PckeP5S6JAJ35kZy7MRM7JyC5mzYScF+47w/PxN5WaC2iJAaoqZ0oyGwXeezWLWuh189tPTyyt4x5RN85yZfNsy2LqkYnnQ/UFEslR1XOhHI41UfgDsq2LdRXUQ1wgmvZroqtD14eYDzH/0aHNXOA7tcgqiYAUgMPISNxIpKzn6s5kTISm18n7Tx7uR0Ef3wvZV0O2Ymn7ruBNsDktOEp76IgdwyiHJJ6i6BLELjuvFS1n5+NWt69upFZt2H0K9bQd3b8OaggMoFdWRTakYiSQrdzfvrdjGj84aHB+Fkjcfpp1bkTB9zFRXpzDc/SEMkQpKVtk8XVUT+GhvhCV9vHOaRxOVlRN0acUH3Ya6aK9wn61qv6OvgU9/73wxX/9zfL5TjDjv2F48P38TJaV+UpJ9/Oq84RQeKi63Qb/x1Zbydd+eNID731pRPn/tSf24f8YKijxTWpc2LRL8bYzmTFbObm5/YTEdWqVw08Q4hRBnf16hUCQJeo2Ck74fdcRnraK/jAZKdaOeAFWNPqr6bLh1rTvD8ItczsuZ90Bq27rLHyeCo8PC2Z5D1w3p0fao+U9WF/DSwjx+985qsncc5OzhPWzEYtQrWbmFXPbYXErKlJQkYdXW/fH5Dfo8tSC+6O4PIVTpU2kuJNynkihikWuSvxAePwO+/iAcf1Ns5WuAvLlkM7e/sASAtGSfZe4b9cpDH63lrx+tAyBJ4EdnD+HW0wfG/kD//SZsWQInfAf6T6ry/lAbnwpe18dvqeqLMRHWaDjU4MmjSnqPhZ7HwRd/h8N7oV/ikiHrg7zCw/jE1RsrKvUzZ8NOUypGvbF5z2GgDiHE0TxIFua6FhmTfgaTflorOavLqPd7IbumVIyjEYEBZ8Dsv7gIsQaSDBkvAjXHjpT4UWDd9poUfzCM2pO3+xBvLNnCSQM6c/LALjUPIc6bD9O+7pztwcnSoSx6xv2vx9S+u0k0PpUPReQnuLyVg4GFqmqNLAz3AwXA36CSIeNBsG9maf4e3liyhb6dWpGWkmQ5LEZcuW/GSpJ8woOXHEfP9rWoaL7hU/f/BJcs/fkf4KQ7XHJ1YORSVuLCiAedDe371FrWaJTKDd701qBlCvSv9VGNpsOAyfDZA4BGFW7Y2AmEKpeU+fnGw1/wj0/W47NyMEacyMot5Ll5uXy0qoCfTxlaO4UCcCSQHSJuJLL+I/dCKkYu+7fBgQJXVaMOVKtUVLUeS18ajY708ZBxEuxYDZe/0GRHKaGkJPk4bUg3lm3eh1+hxHJYjBiTlVvIlY/NpajUjwBj+nao3Y7KSmDlG9D1GBh5sfN9LnkOsp4C1I1cZv/VVdhIbQdp7eskd5Wl7wOISCsR+aWIPOrNDxKR8+p0VKNp0fdEOLzHOe2bEacN6UaLZPcXUuCEfp0SK5DRpPhyw06KSv2AG1wszC2s3Y6W/g/25sFZ98GpP3YPfqOugOSWLmxYfLDmHZc5f2S/6/pah9Yb1SoV4ClcL5WTvPl84De1PqLR9Og61BW43LUh0ZLUK2MzOvL8zRM4e1h3/ArzNpqb0YgNfr+SleOUiODMq7UqGOkvg1l/gR4jna8kQCCpefIv4fr34LgrvRVa4RutJdH4VAao6qUicjmAqh4WkWjaCRvNha6D3XTnGug+LLGy1DNjMzryn6vHcscLS3jwgzW0S0tmX1GpOe6NWpOVs5vfvrOKRZv2cNnx6aR3alX739PK12H3Brh4mhvuBBOcViDiusJGU+apGqJRKsVe6Xp1x5YBuE6QhuHoPAgQ2LEm0ZIkBBHhd98YyYKc3fzfGyvMcW/Umqyc3Vzy6FzK/EqST7h4bB/GZtbSrLppHrx3F7Tr4+p3RaImZZ6qIRrz1z24lr/pIjId+Bj4Wa2PaDQ9WrSCjhnOWd9MaZOazNnDXJuhYMe9YUTL3kMl/PzVZZQFKmqrMre2JtW8+fD0eS6a60CB67FUHenjXXXyOgbbRBP99aGILAIm4Mx7d6jqzjod1Wh6dB3q2hE3Y6aO6sXz83MpLlP8wKj0DokWyWgEZOUW8vrizby3Yiu7DxSTHFRFu9aNt3JmVTTTU3+95o9FW1ByEnAKzgSWgmv9axgVdBkMGz6BslJIap51SsdmdOT5W07k2Tk5vPnVFn7/zirOGtadUwZ1NTOYEZYFObu5/NG5lPoVAX73jZEM7t627p1Jg30i9Zw/Vu2/X0QeBgYCz3uLvi0iZ6rqrRE+ZjQ3ug51Tr7CHOgShyJ3jYRAcmRG51Y89PF6lm/ZxyOfb4i5f8VaIjd+1hbs5/bnF5c3kPMJ7D5YXKkXUK3pNABQ6H86nP6Les0fi+aRchIwQr1yxiLyNLAsrlIZjY+uQ91055pmrVQCtEhOQnBD+yMldU+MDFYi+4tKuOXZLErL/BYQ0AiZl72Lf3yynjnZO2mVkkRKkuD319HcFcqmL930tDvrPSE5GqWyBugL5Hrz6cDSuElkNE66DHLTHath6NcTK0sDYEL/zqSmVBSfLCopO2qbaEcb7y3fym3PLQ7bEvlIiZ//Lci1UUsj4cnZ2fz6rVUobmTy98tH065li9hfv5wvXHJjrzGx2V8NiEapdAZWiUggxfJ4YI6IvAmgqtXEqhnNgrR20K53sw0rDiVQfPLLDTt5f/k2/jMzm7OGdefYPh0AeGlhHj9/ZSl+hWSf8PtvjCSjcyvmb9zNyN7tOXikjPdWbGPT7kMsydtTad+j0zuwfMteSsoUBV5cuBnBlY6Zdv3xpKYkmZJpYLy3fCsPfrC2UmVrAVZu3c+tpw+M/XXKnQ3px0Ny/XcqjUap/CruUhhNg65DTKkEEbCNXzG+L1P/+QXXP7WAM47pRs7Og8zPqSi5UepXfvpy1YP/yUO78sX6XZSWuRbHvzzPJZjOzd7F8s17eXf5NhQoLvNzxePzCKS4paZUYRqLRYM2IypUlT+9v5qHP8sGXHOtJJ+PMr8/tuauYA4XwrblcNpdsd93FEQTUvx5fQhiNAG6DIFFT4PfD75oUqCaB53bpPKDMwfx05eX8uLCfADOGNqNL9bvpKTMT3KSj3GZnfhy/U5CDVxJAmMzOnHr6YOOGn2MzehIVm4hn67ZTkmpH59PGNK9Lcu3uIq0RSV+PlpVUFmp5M6Bp893YaZJLZp0/5tEc+BIKXe/tow3lmyptPxb4/rQu0PL+I0kc+cACpknx37fUdA8Yz+N+NB1CJQccsXrOmYkWpoGxfb9R8q7RiYJjMnoyPdOH1iuKACycndTUuonyefKk5eVVTzNVhURFNzjJbCfKx+fS3GpH7/Cc/M20at9WkXpmC//Dn6Xv6BlxUgT6H/TECPhXlqYx+/eWcWeQyVcMT6dVxdvpqTUXc9vjukTXzlzv4CkVOh9VKffesGUihE7yiPA1ppSCSHQNTJwYwmnKEKVQ7Q3yqr207tDS/743mr+740VALTzFfFl6kza4JTbEU0iN+04hsb+69YbH68q4OZnFqIawdxXjxQeLOYXry3j3eXbAGiR5OObY9P55tj0+lN8ObOhzzhISYvvcaqgRkpFRDoC6apq0V/G0XQd4qY7VsOgsxIrSwMjdERR1agjeHltbz7B+9m48wAPfbwegJt8b9JGD7BD21KiKdxWcjuDcrvS98D6BvWUHy17D5VwpxfsAM7c9985iYmEm79xF//+PJt52bs4VFwR6Vfmd+HkcXHGh6NoL2xbChN/Ev9jVUE0yY+fAVO9bZcAO0Tkc1X9UXxFMxodrTpB667NugZYJGKS1FZDTh3cjf/MzKZT6U5uTnqHnB5T+GRLElfIB3ylA1i0ML+8tPpzN08Aoh8hJZJDxaVcP20+hYdKaJHko6TMhW6/tsRFwtXnqGXWuh1c++R8/OrChG+fPIj/zNpQaVQaN0KDLjbNc/6yBPlTILqRSntV3SciNwFPqeo9ImIjFSM8VgOsQREYIbV89w5StyuZlzzA6Us+I23m21zRv4hns1u7BM1SP9c9OY9DxX78qg1ayczN3smdrywjd9chHrlqDF3bpjE3exert+1jxldbvbwgP+8t3xp3mXfsP8KPX/yqfLQkQAtPocX9vOXNhyfP8YIuUuG6t1wosS8F+iTOTxaNUkkWkZ7AJcDdcZbHaOyktYd1H7onpr4nWPhqvKjBeR27533Y9gYM/yZ06ke/kcUwE24YsI+X8tpSXOrHJ0K7likMKl7KhKRVzC07hhum+ThwpAz1lEyi/RXgnPJXPDYPv0JKktC1bVr5CDArt5APVxaUJ5w+9WUOB4pK6dkhjZMHxrb+WlZuIe+v2MabSzazxxstBYcJ18uodO37TqEAlB2B9++GQzuhbQ8oWJ6w/1s0SuV+4H3gC1VdICL9gXXxFctolOTNdz90fwk8NQU69ofCbFCF5DQLX40VuXNg2tcBdU+okc5r3nx4/bvu/Zq33XzvsZDSin7F65l+0zfLn6hbb8+i/1u/JpkyjtCC7yfdy0f+TMA9+f/unVWcNKAzpw3pVn4Tr+9RzLQvN5aPCvx+rVT+JthvNaR7W56cvZHnF+QB8NDH67l/6nCG9GjLvI27w8oc7fdxim0uR7xWv7+7aCRDesSgCGRNCTjixQcI5Afy0wWenpqw/1s0eSovAS8FzWcD34ynUEYjJWeWaysM7gnq4I6gJ6niei2/3aRZ+kLFea7uvK58I+galFRs22MkbP2KsV8LeqLOXYJKGQKkUcp9xxUye15/jpS6J/+s3EKycgv55yfrGd23A0vz95abyupjFLPzwBE+W7MDEdcIKpy/IniEsKZgP3M27kIVyvzK3a8vL98u2Sf85qIRZHRqxax1O/Gr8sTsjZSWabX+mNcXby5XKD6BwkMxKgJZU/ZthaQ0OPUn0H+S60W/4HEqtQRuiEpFRAYDjwDdVXWEiBwLTFVV61NvVCZzontyDrQkPet+eO/nUFrk2pXWY/ntJk3xoYr3vqTI5/Ww1+RJkiqXQO95HCx5rnKiauuu5dn4IkLvUWcz/ViXdLllz2Gen78Jv7oimcs27y2vRVZcWveCmdWhqvzf68s5UuLnb5eMIn/P4WpHBRP6dyY1EMad5GNMRke+3OAap5X6lTtfCV8Xt6jEz8y1O8Lue272Ll7OcsENInXoHR8LcmZBv4kw6acVyxZPj0lL4LoQjfnrMeCnwH8AVHWpiDwHmFIxKhOuJWn3Yc78cniPM7sYdWfrV9DjWNi3GVq0qfq8+ssg+3OXBDf03Mr+l57HwfxHYXd2RVXpPZsAH/Q6DrYshtR2jO1W4a94ZVF+eUTTr84bzn0zVrhRjLp6ZPHkraVbeXf5Nn42ZQgXjO4d1WfCJYYuenwuJaWuisH4zE7M9qoYCODzuWrBCryclUf/Lq3Llddxfdrz6MxsHvxwLT3bp3HPecNZu31/4gIY9m9z+WCjr6pYFsOWwHUhGqXSSlXni0jwstI4yRMTPL/P3bjItW8lWp5mRfr4yj/m9PFwxq/gxWtgw6cw6MzEydYU2LXBtReY8gC06QYv3wArXoORYX7m2Z86xXPOb2H4RZXX9TzOTbcuqVAqa96FjBPh0v/C346Fz34HlzwDhM+zGdKjLf9bsIkXF+bzxpItnDSwS8y+Zrl/o18ntuw9zE9fXsrAbm24ZWL/Gu2nugTTBV4Vg4CiLDxUTJvUZP7y4Rru+N8SwFM4AmWeL2fH/iN0atOCW4cnsMVDzmw3DR2NhP7/EkA0SmWniHgdX0BEvgVsre5DIpIOPAP0APzAo6r6UG2EFJEngfOA7ao6ImTdFOAhIAl4XFUf8Pw+N4rIy7U5nhFjBn8NWnWGxc+aUqkra99z08FToEMGdHsQPv0dDLvw6I6bi6dDy44w5Nyj99N1qDORbP3KKaQ9m2D7Cjj7Ny7f6MTvwed/cOs9BRQuOXNsRke6tk3lX59u4KSBnblgVHSjiFACSuSEfp1YW7CfX72x4qhS/3m7D/FV/t46jQwiKZng5dv2FvHI5xsAd+Pr2T6NzXuKUKC0LP7mvmrZ+Dmktq94OGhARKNUbgUeBYaKyGZgI3BV5I8AbjTzY1VdJCJtgSwR+VBVVwY2EJFuwGFV3R+0bKCqrg/Z1zTgnzglRdC2ScC/gLOAfGCBiLwZfAyjAZDcAo69FOY/Bgd3QesE2aCbAmveha7HQKd+bn7y3fDCFfDy9XDS9yueUg/thtVvwbgbIDn16P0kpUD34U5pAKwJKKuvuemE78G8/8A7P4XB50Q0p/zwzMHMy97Nna8sY/nmvUwZ0bPayKqs3ELmbNjJ8N7t2X+4hJ+8tLQ8gbEq4nEzr8rBfuaw7jz15cbyUcz3Th/E/W+tqJ+ExmjYOAsyTnI+tQZGNNFf2cCZItIa8AUrgGo+txVvRKOq+0VkFdAbCL7hTwK+KyLnqmqRiNwMXAScG7KvmSKSGeYw44H1noyIyAvABSHHMBoCo6+CuQ/DspdgwncSLU3j5HAh5H4JJ99esaxVF0Bg1Zuw7gO4doa7+S972TlsR0d4/ut5HKx43YV8r30XOg+sMIW17AAjLoKFT0H+goihy8lJPm4+tT/ffjaLx2Zt5Jk5uTx3c0X0VFZuIZc9OoeSMteHvU1qEvuPHN20LMBJ/TuxaNMeSsrCF9esD6oy9zWIRNC9+VC4EcbfkjgZIhBN9NcdwFPAfuAxERkD3KmqH0R7EE8hjAbmBS9X1ZdEpB/wgoi8BNyAG3VES28gL2g+HzhBRDoDvwVGi8hdqvr7MDKdD5w/cKC1vq0Xug+HnqOcCeyEb7vQGaNmrP/YhRIHm7NyZ0OgcXFpkVMC6ePdee55nAsdroqex0HWNChY4Z58Q5V96+5uqv5qQ1TXbz9QXoX5SKmfj1ZWlNz/56frKfEcEgq0Tk0uVyoCnDSgMwtyCsuTB398jitxWZvimrGkKnNfwtk4y037NcxoymjMXzeo6kMicg7QDbgep2SiUioi0gZ4BfiBqu4LXa+qf/RGGI8AA1T1QOg2kXYfZpmq6i4g4uOwqs4AZowbN+7mGhzPqAujr4J3fuJMKsdeknCHYqNjzbtuZBIc7ZU50Zm3yo64EcdXz0FhjisqeOJtkfcXsMfP/qtLWA2YvgIMPANm/Rn8peBLjhiiGqjCHCi5/9rifEb17cCjM7PJyi3EJ+7PmpLs4/YzBlcyJf3obFeINFy/mAAN4mbeUMiZBS07QbfhiZYkLNEolcCN+1xc7a+vRKJ7zBSRFJxCma6qr1axzURgBPAacA9QzT+hEvlAetB8H2BLFdsaiSbgB1jwOCz+r2XY14SyElj/IQw9r7IdPTiMNP0EyHoalr3o1i14HIZdUPU57jbc5a8sfwXSOrjPB5M+Hq58CZ67DPocH/FaBZuLOrZqwa/fWsm3n80CIMkn3Dd1OHsPl0Q0JZniiAJV2DjTFYxsoI3wolEqWSLyAdAPuMtzuvur+5CneJ4AVqnqX6rYZjQuD+bruACA/4rIb1T1l1HKvwAY5JnQNgOXAVdE+Vmjvgk4hROc8dsoWfiUK2veOYy5NjiMNG+eK9uh/soZ9OFISYMOfZ19vvfYo6PHAAZMdgEAsx50raID7Q3CEGweWrFlL9PnbXIrVNl7uIRbTx8YdlujBqx8wzXBO2ZqoiWpkmhU3Y3AncDxqnoIaIEzgVXHycDVwGQRWeK9QmMbWwEXq+oGVfUD1wK5oTsSkeeBOcAQEckXkRsBVLUUN7J5H1gFvKiqK6KQzUgEmRO9OkUkNOO30ZE3H973+o1//gc3XxWBqgahGfRV7XePd+PPmVX1fid8D1Jawqywz4Zh+caYPqSl+EiS8OVUjFqQNx9eucm9X/h45N9BAokm+ssvIhuBwSISdSsxVQ14ECNt80XIfAlu5BK63eUR9vEO8E60chkJJH089J8MmxfAlS/bKCVasj9zfg2ofvRRk6zqnFkVdcH8ZVXvt3VnGHs9zPs3nHZnhRkzAtE0JTNqgGqF7wugrLTBjvSjif66CbgD569YAkzAjRomx1Uyo2nSfZiLWOpzfKIlaTyUHHZT8UU3wos2qzpzoqseHU2tqJO+D/P/Ay/fCF97IKr9m4krBgQqf+fMhry5FRWJG/BIPxqfyh3A8cBcVT1dRIYC98VXLKPJ0qa7C309ss/1XjEi4/fD6rehU38YdZULI43V02lNRjV789zT8pYsePr8inwYI37kzXfnurTIzZ/wXVc5YdMXDbo/UTRKpchLTEREUlV1tYhU7a0zjEi08XIfDmw3pRIN6z5wtb6+8ZgLw4410Y5qcmZBIN+99EiDNb00KXJmuXMNboTSpitkTHCvBkw0jvp8EekAvA58KCJvYGG7Rm1p081NDxQkVo7Gwpd/h3Z9ji4IWd8EAgAAUFf52IgvmRMrwoaTUhusuSuUapWKql6kqntU9V7g/3BhwhfGWS6jqVI+UjGlUi35WZD7hSvumJSSWFkCprKxXuDnLmv+GnfSx0PXYdC+T6PK6apWqYjIBC83BVX9HPgUV3LFMGpO+Uhle2LlaOjkzYc3b4OU1jDmmkRL40gfD+f91eW0zHnY+XuM+KHqQr4Hnd1oFApEZ/56BAgunXLQW2YYNSetA/hSTKlEIm8+PH0ebF/pyq9sX5VoiSoQgRNvhd0bYN37iZamabM3H47she4jqt+2ARGNUhFVLa9I7SUpRuPgN4yj8fncaMWUStUsf6XCQavqOckbEMdcAO3TYc6/Ei1J06bAy+NugkolW0RuF5EU73UHkB1vwYwmTJtu5lMJR958ePE6WPCEm482L6W+SUp2laZzZrnioA00s7vRU7DcTbsdk1g5akg0I47vAH8HfomLKfwYaJiF/I3GQZvusM8CCCuROxemnetK24sPzrrfZdE31HyE7l6F3PmPwqJnG5UjudFQsMJ190xrl2hJakQ0ZVq24wo1GkZsaNMNtixJtBQNi6XPO4UCgDiFMvHHCRUpIlsWU97HpczyVuJCwYpGZ/qC6MxfhhFb2nSHgztcvSnDkdrWTaMpBNkQCPRxAef36XtSYuUJkDffVVQOmORC5xsLJUUubLt7w+yZEglzuBv1T5vu7qn80G6XJWw4x3xKK5j4k9iWYokX6eNdqZa5j8CKV539P+PExMpUXtbkiOs5M+gcF6Hm9zsF2JhMdDtWu2KfTVGpiEiSqtojpRE7WnuK5ECBKZUAO1a7G8ipDdjkFUr6eFcY9HAhfHy/ayDWrmfi5Ckva6LOfLjm7Yp1ja1/TyON/ILozF/rReRPIjIs7tIYzYOqsuobq6kiFlTTAKvBIgJff9DdzF+9pX6uX1W/k74nUV6fLDkNzn3QTcEtzzglvnLFkoIVkNwyqjYDDY1ozF/H4hz1j4uID3gSeCFcv3nDiIpwWfXlpovixmeqqCuHdjsF23VooiWpHZ0HwKjLIWuaK9Eez+uXNx+mnef6iiSFHOfwbjc99lI4/ia3vOexMPNBWPce7N8ae3niRcFyF0oc3Dq6kRBN7a/9qvqYqp4E/AzXR36riDwtImF6mxpGNQRGKgeDlEq56cJfYapoLuxY46aNVakAtOvtvYnj9VOFz37vos00zHGypkHbnnDBwxWKJn08XDYdeoyE938BRw6E3XWDQtUplUboT4Hoan8lichUEXkNeAh4EOgPzMA6Lhq1IbWNq2kVPFLJnOhMKeCezhp69FMs2bHaTRuj+StA/9Pi1yo6bz58/keYfjFs+KTiOCIVx9mTB+s+hNFXueTMYJKS4et/gX2b4a0fNnwT64ECOLTLKcJGSDTmr3W4IpJ/UtUvg5a/LCKnxkcso8kTmlWfPh5ad4MD29yNobmYvsCNVFJauxL3jZX08XDsZbD0Bbj61dhdv9BGVWOuhVFXwjs/djXRWndxyxf/101HX121fIPPgWUvepUKGrCJNZBJ31RHKsCxqnpjiEIBQFVvj4NMRnOgTffKSqVon1MoAEV7EyNTotix2o1SfI08baz/JGeWatUldvvMmVWhUPBBxwzoewJc8ZJTDO/f7fq1L34WBp7h1ldFd+/JP5zprCERiPzq1jhjo6L5FZeKyK0i8rCIPBl4xV0yo2kTWlRy+0o3bdEGti1PjEyJYseaxu1PCdBlkJvGstdKizbeG3EBAAFzV7ueMOmnsOYdeP8uZ9oae13kfQ0+J34mulhSsML5qFp1SrQktSIapfIs0AM4B/gc6APsj6dQRjMg1Py1bZmbDrvQ3ZRKDidErHqnaC/s39K4/SkBOntKZefa2OzPXwaLnoHW3eG0u442V034HnTq7+qPpbSGVp0j7y99vDOdAVz+fP2avnLnRO/LyVsAKS0btt8nAtEolYGq+n/AQVV9Gvg60Dg9SEbDoU13lzQXKPFesNz1Whl8tjNPNKQeIvFkh3cDbgojlbR20KYH7Fwfm/0tftb9Ls79A5z286OVQHIqjL3BvS85BM9+o/ob8bAL3VTq0dS46i14agp8/Gt4empkGb/8FxRmw64N1W/bQInmzJZ40z0iMgJoD2TGTSKjeRDIVTm4w023LXfRLoEM4oJmYgLb4SnPpjBSAWcCi8VIpWgffPIb6HtihSIIh78YdxvT6Pwkvce4af6CussYLSte894olB6G7M/Cb7d4Onxwd8W2DdnvE4FolMqjItIR15/+TWAl8Me4SmU0fYKz6v1lzqfSfQR07OdMGc3Fr7Jjjcuc7tA30ZLEhoBSqejr56hptYS3f+weOEZdWRFqHo5AYctoC3G26gSdB8LmrOjkiAXJLb033vf46n/w6W8rzsWmefDMhfDG96DXKFcFoLEUFg1DNKXvH/fefo7LTzGMulOeVb8Ddm905oseI1wEVPdhzWikstrdiBth5nRYOg+Coj0uzyIQ7hvIgi874m6Y186I7M9Y9ZYL/QXXBKzrkKq3Tx/vfC05s6LvPdPneFj/kVN8kRRWrDiyF9r2hvE3wq5sWPJfl3fz+Z+gfR/XNhh1iuTs30JSSs2+TwOjSqUiIj+K9EFV/UvsxTGaDcEjlVLPKR8wfXUfActfrb8/fSLZsQYyGkjZ+FjQZbCb7lxboVRyZjlTDrgyPNUVdgwoFIiuEGT6+JrdfHuPha+ehz2bIocgx4rtq6D3aNcfZ9aDOAORH1AoOUh5vTKAvLluu0aoTAJEMn+19V7jgO8Cvb3Xd4DGGUBtNBzKKxVvd6YuSapwVvcY4Z7u9uYlTr54EGoCOrLffcem4k8B6OJVbtoZFFYcbMIRqjfpBO6x8TIB9TneTTcvjO1+w1FSBLuzK3JOgs11yS1h8j1u2ojNXaFUOVJR1fsAROQDYIyq7vfm7wVeqhfpjKZLciq07OhGKnvz3BNuildRNpCkVrCi6fgaApnhZcUV2dwBk1dTiPwK0D7dmbiCnfUdMgB131v91SvR/Vtc0Mbwi+JjAuo+3MmYvxBGfDO2+w5l1zrXOyjQZz6cua77sEZt7golGkd9X6A4aL4Yi/4yYkFrL1dl23I3OgnQ3Xuqa0rO+kBmuPpdGHXOrKZRSDIUXxJ0GgC7gsKKAxFM5/zWVRde9VbVny8rcTlL/SbFzwyUlAK9RjulEm8CofHB2fHp4yt/t9D5Rk60yY/zReReb5QyD3g6rlIZzYM23dyNdV9+5WZEqW1dFFjBssTJFmsqtdv1Q/FB56RPSvWe5JsQoWHFG2dCansYez10zIRlEQwdO1Y75dtrdHxl7D0Wtn7lfDzxZPtK8KW49gDNhGhK3/8WuB4o9F7Xq+rv4y2Y0Qxo0x12ek/rPUI63PUY0bRGKgHT3tDzoNcY51vJegZadoAtixIqWszpMggKcysSWzfOhMxTXLXgEd+CjZ9XLtETzJbFbhpvpdLneBeNVpMHl9o0kdu+yp2PpJSay9hIiSqtVFUXqepDwGFVXRxnmYzmQiACDI5um9p9pHNwFh+sX5niRc5sNz33z3DDe27kUlTozH+NNHO6SroMdn6E3RtdhFXhRujnFTQfebEzAZYnBIawZYkb1XSMc8fDPuPcND/KfJWNs+DJc1xCZk2u1/aVFf6UZkJNaxV8Jy5SGM2TQK5Kqy6VFQx4IxeFgpX1LlZcyPnC+Rra9XRBCgMmU54M10gzp6uksxcBtmuduxlDhVLpNtQ9QFRlAtuy2HVrjHfF5na9XUOvaDPr5/zTKcOaVDg+csApVVMqEWniSQNGvRJQJD1GHJ2PUl6upQn4VfxlkPulMwEF6D+p0WdOV0mXoMKSG2e6h4bgG+vIb7mb+e6NlT9XWuySXuNt+gL3e+s9Nrqw4pKiyson2usVCMRopCXsa0tNlcr5cZHCaJ4ERiqlRUebEzr0hZRWrh5SLExDGz51vcrruq/cOfDp72u2n23LXN5N8I0oEFo6+e6G2yyqtqS2daOAneucUuk3sfJDQyCM960fVj6P21e6UUB9KBVwfpXd2a4kzLKX3euT3xx9bRc94yoEdOzvFMrVr0d3vQLtHGyk4hCR20UkPXiZqubHXySj2XB4j5tumne0nTp/gVM2mxfW3eew9gN49sKa28NDyZsPT58Hnz/gck6i3U/uF26aeXLl5U0slLQSXQa5Uij7t1SYvgLs3+aqBGd/Wvl61JeTPkBaOzdd8Di8cqN7zfwTTPt6hUwlRTD7r84HdsodTum17Rbd/rev8uq6ZcZF/IZKpJHKr4F5IjJLRL4nIl3rSyijmbAnxytBHqYia86siqKEdfU5rH3Pe1PHjn85s8BfWnOZcma7vh/tetXuuI2RzoMqKlD3m1R5XfC1DeTsgFMqaR1c2HF9cKiQiltg0EiqrLiiPfHiZ51iPO3n0ONYt2zr0uj2v32l8yE19o6eNSTSt83GNeT6NTAWWCki74nItSLStl6kM5o2mRNdnkY4v0LmxIowTF9y3XwO5Vn5Ujf/RfoJQTMS3X78ZW6kEuxPaQ4EaoCltoODOyuvy5zo/EkAKGR452brElelt77qvfULrnCcUvFbRFwl4XmPwkf3QbfhTjF2G+bWb4vSz7d9VbPzp0BkpaKq6lfVD1T1RqAX8DAwBadwDKNuRPIrpI+HS6e796OvqpuJKK29m3bMqJv/IqmFm6a0cm1ue4+t/jMFK1x3x4xmplS0zE2P7INnLqhsKgxc92OmAuqqGpcUuUi/+jJ9Bcsx+W647m247i33/vLnoXVnePenULzfRbHlL3C5Rl2HRKdUDu2GA9uanT8FIpe+r/S4oKoluH4qb4pIy/AfMYwaEqnC7OCzXXjqvi11O0bADIPUTTltmuumZ9wD7/0c8uZVX2E4kJ8S6k9p6gSPTsJVGk4fD996Ev4+Bmb+2UWI+UvqV6kE5AiVC9y1nf1X995fViF/j5EVYdKRCJRn6dr8lEqkkcqlVa1Q1WbSQNxIOOkT3B88tOlTTThQ4KZ7890NorbkzXP2/lFXuFHL6rer/0zObJfI175P7Y/bGBnytepDppNS4OTbIX8+fPl3t6y+lUpVDDk3fPXgHiOdjyXUpBdKM438gghKRVVj0BPUMOpI3xPg8O7KpdRrSqAkiL/ERR7VBlU3Uul7oosa6jcJVs2IrOxy58KGj6Hr4NodszGTPt4146ouZHr01S5faeXr7ia+b2u9ilklVZlme3gVtKszgW383Plo9ja/gNnmFZZgND7SJ7hp3tza7+PAds8BS+17tOzOhkM7K5z1x5wHe3KdzyQcefPhmakuLHrDJ02rDEu0RBMynZLm+VZwzdpC/S+JJJz8gQiwSEolb74bxZYdaVjfp54wpWI0bLoMgpadXC5LbTm43fXQAFc2ozYE/Cl9PSU3+GuAVG0CC+526Pc3rTIssaZV54r3Db1kTatO0K4PbIsQVpz9mSvnAg3/+8QBUypGw0bEjQ7qOlIJRGrtya3dPvLmuiiyLl6Dqbbd3RPs6hnht8+cWNGEq6mVYYk1A89oXN0Pe4yMPFIJRBuKr3F8nxhjSsVo+PSd4Jo+HdhR/bahFB+E4gMunLh117qNVNJPqJzINvQ8d3P54P+ONnGkj3f+l1adml4ZlljT2ErW9Bjp6pqVVBGvtGeTy6069WeN4/vEGFMqRsMnYHLKq4UJLOCkb93NJUHuqYVP5dBudxOplPxIRUTXl/8IX/6laC/0Gtvsbiq1ojGVrOkx0pm3AhFeoaz70I1OTr+rcXyfGGNKxWj49BzlzAi1MYEFclTaBJRKLUYqAWUWUG4BCgNVdsOUmVGFwhzoFOe+IEb9EykCrDDHNZ4bdHa9itSQMKViNHxS0lz+Qm2c9YEclYBS2ZvnHOc1YdNc1xK215jKyzMnVkSVhdrOD+122eTxbjZl1D8dM135mXBKZd2Hbjr4nHoVqSFhSsVoHKSf4GpDlRTV7HPB5q/26W5EcbCKVrZVkTcPeh4HLVqFyDQeRl3p3l/xYmVTR2GOm9ZXcUSj/hCp2lm/7gNXPLQZ9aQPxZSK0TjoO8EphPd+XrO4/wPbAYHWXaBDhltWExNY6RHIX+j6q4c77sDJbtqyfeXlAdOYmb+aJq27wuZFrr9OgOJDrn9MMzZ9gSkVo7Hg8yoWZz1ds54oB7e7CKyklIpqxTVRKrP/6jLxN80Pf9xAiPGOkAIUga6GAUVmNB0CyY3+EpfgGvhN5Mx2ya6mVAyjEVCwDFfjVCv34KiOA9sr2hZ38HrORatUtn7lih0CVfZi6TzA5SPsXFN5eWEOtOlxtMnMaPzkzKqc3LjyDfd+3QeugnVGMyseGoIpFaNxUKkHhx9aRNnS58B2Z6oAaNHaZW9Ho1T2b4PnL3dNoyIVRkxOdc74HaFKZaOZvpoqmRPdbyEQpLHiddfFdN370P80F1jSjIlU+t4wGg6BBLl1H8JXz8PH97vs+GEXRM4FOFBQOb8kmrDijbPgtVvg4G646UNn0siZ5W4m4Y7VZfDRBS93b3Q3GKPpEfgt5syC1Pbw3p3w2GT3uxp+UaKlSzimVIzGQ6D3Rfdh8NJ1MOefrr/4tTPC3+xVPfNXUE/x9umwY3XVx9g0zxUB1DL3NFpaFLnnC7gqxOs/grJS59AvOezKo9tIpekS/JvYtxlm/8W9n/cfV2mhGSY9BjDzl9H42J1N+U83UsG+4gOu8m2wUgmMVKoqWf/5Hyq6FgaaM1VHlyHOaRsIIw6MhCycuHmQ2obynoZlJc2ugGQoplSMxkfmREj2WvuKr+qCfcE5KgE6ZLjRx8EwdcTWf+zK1EtSzYobdvUiwALO+kDklyU+Ng8C/r7GUhAzzpj5y2h8BBpATb/Y9beoytQQUCqhIxVwNcCCl+/OhpdvgG7DYMrvYfPCqn0ooXQZ5KY71sDQr1uOSnMj2McS7W+mCWNKxWicpI+HzFMqeoGHI7hES4DysOJc6OOVw8/+HF65CfylcNl0pwz6T4pelrT20LanKzoJzgzWom3lPiFG06Y6v1szwsxfRuOlx7FuhHHkQPj15cUku1csax+Sq7JpHjx7kUuSLCsObxaLhi6DKsKKd290/hSR2u3LMBoxplSMxkuPkYBWXYL8QIHzuQSPGNLaQcuOFW2FZ/+15o75cHQZ4sKKVb0clcza7ccwGjlNSqmISH8ReUJEXk60LEY90GOEm1bV2vXAdmjVpaIDY4D26W6kkr/QZUGLr+5O1q5DoHi/Cy8tzDUnvdFsafBKRUSeFJHtIrI8ZPkUEVkjIutF5E4AVc1W1RsTI6lR77RPd/6Mqlq7huaoBOjQ1/liXrwW2veGK16qe9fBLoPdNPtzKDtiTnqj2dIYHPXTgH8CzwQWiEgS8C/gLCAfWCAib6pqFXYQo0ki4vwq25aHX3+wCqWSlOrMX75kuPFD6D0GBp1ZN1kCYcXr3ndTy1ExmikNfqSiqjOB3SGLxwPrvZFJMfACcEG0+xSRW0RkoYgs3LGjlo5Zo2HQfQQUrHD+kFAObK+cowKuouyqN70ZcRFfsaBNd1eyY/0nbt7MX0YzpcErlSroDQQ3G88HeotIZxH5NzBaRO6q6sOq+qiqjlPVcV27do23rEY86THSZc3v2lB5ebgSLVC5wqz6Y5f9LOIiwIr3uxFQIMrMMJoZjcH8FY5wsZqqqruA79S3MEYCCfQLL1jmanAFKNrrfBuhSiVQYbasOPbZz12HuKTJ9umuBphhNEMa6y8/Hwh+FOwDbEmQLEYi6TrUNfDatgxGfLNiebgcFYhv9nPAWW/+FKMZ01iVygJgkIj0AzYDlwFXJFYkIyEkt3AjhNAIsEA2fesw5s14ZT8HnPVH9jnfjWVYG82QBu9TEZHngTnAEBHJF5EbVbUUuA14H1gFvKiqKxIpp5FAeowMo1QCdb+6H719vCgtctPNi2rW8tgwmhANfqSiqpdXsfwd4J16FsdoiPQY6Rp3BTvmwxWTjDe7sylveRwoyW+jFaOZ0eBHKoZRLd0DmfVBo5WD212WfMtO9SeHlUA3jIY/UjGMaglEgG1bBgPPcO8PFDh/iq8en5usBLphmFIxmgCtOkGrrrDsJcg4yd3Md3g94+vbYW4l0I1mjpm/jMZP3nw4vAsKlsOTU+CxMyB/HhzYZg5zw6hnmq1SEZHzReTRvXv3JloUo67kzKroOa9lTrkEiNTD3jCMmNNslYqqzlDVW9q3b59oUYy6EuwgT24JU/7gpuYwN4x6x3wqRuMnnIO8+zBzmBtGAjClYjQNQh3k5jA3jITQbM1fhmEYRuwxpWIYhmHEDFMqhmEYRswwpWIYhmHEDFMqhmEYRswwpWIYhmHEDNFAJnIzRUT2A2uCFrUH9oZ5X918TbZN1Gcb6nG6ADsb8Pepy2eb+nGCr11DldGOE/1na/JfHKKqbQlFVZv1C1gYMv9ouPfVzddk20R9tgEfZ2E9HacxnIvGdpwq/z8NSEY7TvTzNfkvVrr2gZeZv45mRhXvq5uvybaJ+mxDPU59yJSozzb144TSEGW040Q/X9ttyzHzl8hCVR2XaDmaM3YNGi927ZoWNbmeVW1rIxV4NNECGHYNGjF27ZoWNbmeYbdt9iMVwzAMI3bYSCUIEZkiImtEZL2I3OktGyUic0VkiYgsFJE6VSkUkSdFZLuILA9adq+IbPaOsUREzo3Bd0kXkU9FZJWIrBCRO7zlfxKR1SKyVEReE5EOcTrOcSIyR0SWicgMEWkXg+8U7vrE9NxVcX1ies4iHCce5yzs9fHWfd87nytE5I8xOFa46/Nr77wtEZEPRKRXHY8R7rx1EpEPRWSdN+0Yg+8S7jj/C/qd5YjIkhgcp8rr463/iYioiHSp43HSRGS+iHzlHec+b/nF3rxfRGJjxgznvW+OLyAJ2AD0B1oAXwHDgA+Ar3nbnAt8VsfjnAqMAZYHLbsX+EmMv09PYIz3vi2w1vs+ZwPJ3vI/AH+I03EWAJO85TcAv47T9Ynpuavi+sT0nEU4TkzPWTXX53TgIyDVW9ctTtenXdA2twP/jsN5+yNwp/f+znhdn5D1DwK/itf18ebTgfeBXKBLHY8jQBvvfQowD5gAHAMMAT4DxtX1+6ha9Fcw44H1qpqtqsXAC8AFgAKBJ8b2wJa6HERVZwK767KPKI+zVVUXee/3A6uA3qr6gaqWepvNBfrE4zi4H+pMb7MPgW/W5ThUfX1iSrjrE+tzVtVxiP05i3R9vgs8oKpHvHXb63iosNdHVfcFbdMa93+qNVWctwuAp733TwMX1uUYEY4DgIgIcAnwfAyOU9X1Afgr8DPqeM68fauqHvBmU7yXquoqVV0T4aM1xpRKBb2BvKD5fG/ZD4A/iUge8Gfgrjgd/zbPTPBkLIbvwYhIJjAa93QSzA3Au3E6znJgqrfqYtxTV12o6vpAHM9dGGJ6zkKI9TmrRMj1GQxMFJF5IvK5iBxfx91XeX1E5Lfe/+dK4Fd1PE44uqvqVnA3aaBbHI4RzESgQFXXxXKnwddHRKYCm1X1qxjuP8kz2W0HPlTV0PtBTDClUoGEWaa4J7ofqmo68EPgiTgc+xFgADAK2IobWscEEWkDvAL8IPipUUTuBkqB6XE6zg3ArSKShRvWF9f1EGGWKXE8d0cJEONzFoZYn7NywlyfZKAjzgTyU+BF7wm81ocIs0wBVPVu7/8zHbitDsdoKFxODEYpwQRfH9xv7G5irIBVtUxVR+FG2uNFZEQs9x/AlEoF+VR+MuyDM3VdC7zqLXsJN8yPKapa4F1wP/BYrI4hIim4H+p0VX01aPm1wHnAleoZWWN9HFVdrapnq+pY3B9wQx0PE/b6xOvchRLrcxaOOJwzoMrfQT7wqmcWmQ/4cSU6aktV/59gniMGJr0wFIhITwBvWldTXpWISDLwDeB/Mdxn6PUZAPQDvhKRHNy5XCQiPWJxPFXdg/OhTInF/kIxpVLBAmCQiPQTkRbAZcCbuD/GJG+byUBMh7xQ/kcIcBHODFLXfQpuVLVKVf8StHwK8HNgqqoeiuNxunlTH/BL4N91PFTY6xOPcxdKrM9ZhOPE+pxVeX2A13G/Z0RkMM65vvOoHURPVddnUNA2U4HVdThGVbyJe/jDm74Rh2MEOBNYrar5sdhZuOujqstUtZuqZqpqJk5hj1HVbXU4TlfxohZFpCXe96ir/GGJhbe/qbxw0V1rcU+Id3vLTgGycNEs84CxdTzG8zgzTQnux3Ij8CywDFiK+4P0jMF3OQVnflgKLPFe5wLrcbbvwLK6RuNUdZw7vHO5FngALycqDtcnpueuiusT03MW4TjxOGdVXZ8WwH9xSngRMDlO1+cV7xhLcWU+esfhvHUGPsY98H0MdIrH9fGWTwO+U9f9V3d9QrbJoe7RX8cCi73jLMeLXMM9iOUDR4AC4P26fidLfjQMwzBihpm/DMMwjJhhSsUwDMOIGaZUDMMwjJhhSsUwDMOIGaZUDMMwjJhhSsUwDMOIGaZUDMMwjJhhSsUwDMOIGaZUDMMwjJhhSsUwDMOIGaZUDMMwjJhhSsUwDMOIGaZUDMMwjJhhSsUwDMOIGaZUDMMwjJhhSsUwDMOIGaZUjHpBRC4SERWRoYmWxYgO73o9GzSfLCI7ROStRMpl1A0RORDP/ZtSMeqLy4HZuN7lUSMiSfERx4iCg8AIr6c5wFnA5gTKYzQCTKkYcUdE2gAn4/qJX+YtO01EZorIayKyUkT+LSI+b90BEblfROYBJyZOcgN4F/i69/5yXO92AERkvIh8KSKLvekQb/ksERkVtN0XInJsfQptRMb7/70VNP9PEbnOe58jIveJyCIRWVZT64IpFaM+uBB4T1XXArtFZIy3fDzwY2AkMAD4hre8NbBcVU9Q1dn1LaxRiReAy0QkDTgWmBe0bjVwqqqOBn4F/M5b/jhwHYCIDAZSVXVpvUlsxIKdqjoGeAT4SU0+aErFqA8ux92c8KaXe+/nq2q2qpbhnoBP8ZaXAa/Ur4hGODxlkIm7Zu+ErG4PvCQiy4G/AsO95S8B54lICnADMK1ehDViyaveNAt3/aMmOeaiGEYQItIZmIyzzSuQBCjuBqUhmwfmizxFYzQM3gT+DJwGdA5a/mvgU1W9SEQygc8AVPWQiHwIXABcAoyrT2GNqCil8qAiLWT9EW9aRg31hI1UjHjzLeAZVc1Q1UxVTQc24kYl40Wkn+dLuRTnyDcaHk8C96vqspDl7alw3F8Xsu5x4O/AAlXdHV/xjFqQCwwTkVQRaQ+cEasdm1Ix4s3lwGshy14BrgDmAA8Ay3GKJnQ7owGgqvmq+lCYVX8Efi8iX+BGoMGfyQL2AU/Vg4hGlIhIMnBEVfOAF4GlwHRgccyOoRpqgTCM+CMipwE/UdXzEiyKEQdEpBfOHDZUVf0JFsfwEJHjgMdUdXy8jmEjFcMwYoqIXIOLErvbFErDQUS+gwuI+WVcj2MjFcMwDCNW2EjFiDkiki4in4rIKhFZISJ3eMs7iciHIrLOm3b0lp8lIlleolWWiEwO2tdvRSQv3qUlDMOIDTZSMWKOiPQEeqrqIhFpi4t1vxAXIbRbVR8QkTuBjqr6cxEZDRSo6hYRGQG8r6q9vX1NwEWqrFPVNon4PoZhRI8pFSPuiMgbwD+912mqutVTPJ+p6pCQbQXYCfRS1SNByw+YUjGMho+Zv4y44iXFjcY5brur6lYAb9otzEe+CSwOViiGYTQeLKPeiBteIclXgB+o6j43CIm4/XDgD8DZ9SCeYRhxwEYqRlzw6j69AkxX1UAdoQLP7BXwu2wP2r4PLvnxGlXdUN/yGoYRG0ypGDHH84s8AaxS1b8ErXoTuNZ7fy3whrd9B+Bt4C5V/aIeRTUMI8aYo96IOSJyCjALWAYEkt9+gfOrvAj0BTYBF6vqbhH5JXAXsC5oN2er6nYR+SOupEsvYAvwuKreWy9fxDCMGmNKxTAMw4gZZv4yDMMwYoYpFcMwDCNmmFIxDMMwYoYpFcMwDCNmmFIxDMMwYoYpFcNIMCJyr4j8JML6C0VkWH3KZBi1xZSKYTR8LgRMqRiNAstTMYwEICJ3A9cAecAOXHuAvcAtQAtgPXA1MAp4y1u3F1dwE+BfQFfgEHCzqq6uR/ENo0pMqRhGPSMiY4FpwAm4oq6LgH8DT6nqLm+b3+B6zPxDRKYBb6nqy966j4HvqOo6ETkB+L2qTj76SIZR/1iVYsOofyYCr6nqIQARedNbPsJTJh2ANsD7oR/0Kj+fBLwUVPU5Nd4CG0a0mFIxjMQQzkQwDbhQVb8SkeuA08Js4wP2qOqouElmGHXAHPWGUf/MBC4SkZZeu+XzveVtga1e24Arg7bf761DVfcBG0XkYnAVoUXkuPoT3TAiYz4Vw0gAQY76XCAfWAkcBH7mLVsGtFXV60TkZOAx4AjwLVzl50eAnkAK8IKq3l/vX8IwwmBKxTAMw4gZZv4yDMMwYoYpFcMwDCNmmFIxDMMwYoYpFcMwDCNmmFIxDMMwYoYpFcMwDCNmmFIxDMMwYoYpFcMwDCNm/D+hHlcAVK36SwAAAABJRU5ErkJggg==\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "df.plot.line(logy=False, ylabel=\"7-day cases per 100,000\", style='.-')",
"execution_count": 6,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 6,
"data": {
"text/plain": "<AxesSubplot:xlabel='date', ylabel='7-day cases per 100,000'>"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 432x288 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEpCAYAAACUUUmzAAAASXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMisyNTQ1LmdlZjY3M2M0NmUsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvfRvX6wAAAAlwSFlzAAALEwAACxMBAJqcGAAAVgJJREFUeJztnXd8VFX2wL9nUukldAiEDlIEg4hKEVTsfcW6tlXXXbdXXd3V9efuquu6zbboqqhYFxFQUQRBQIGQ0HsgBBJKIKG3lJn7++O+SSZhMpkkM5lJcr6fz3zee/eVe+a9mfPuPffcc8QYg6IoitJ4cEVaAEVRFKVuUcWvKIrSyFDFryiK0shQxa8oitLIUMWvKIrSyFDFryiK0siIjbQAwdCuXTuTkpISaTEURVHqFRkZGfnGmPYVy+uF4k9JSSE9PT3SYiiKotQrRGSHv3I19SiKojQyVPEriqI0MlTxK4qiNDLCauMXkdbAq8BgwAD3AJuB94EUIBuYZIw5GE45FEWpXxQXF5Obm8upU6ciLUq9IDExkW7duhEXFxfU8eEe3P0n8Lkx5jsiEg80BX4HzDPGPCUiDwEPAb8NsxyKotQjcnNzadGiBSkpKYhIpMWJaowxFBQUkJubS8+ePYM6J2ymHhFpCYwF/usIV2SMOQRcA0xxDpsCXFvVtfYdLSRjh3YKFKWxcOrUKZKSklTpB4GIkJSUVK3eUTht/L2A/cDrIrJSRF4VkWZAR2PMHgBn2aGqC+UdOcVtry5V5a8ojQhV+sFT3XsVTsUfC5wFvGSMGQ4cx5p1gkJE7heRdBFJBygq8bA0qyA8kiqKolQgJiaGYcOGlX6eeuqpGl+refPmIZEpOzubwYMH1/o64bTx5wK5xphlzvb/sIo/T0Q6G2P2iEhnYJ+/k40xk4HJAAmd+xqPgR0Fx0nPPsCy7QcY1SuJ1B5twii+oiiNmSZNmrBq1apIixEWwtbiN8bsBXJEpL9TdCGwAZgJ3OmU3QnMqOpaHVskMK5fOz5Iz2XSf5bwtzmb1fSjKEo5MnYc5IX5W8OuF1JSUnjsscc466yzGDJkCJs2bQJg//79XHzxxZx11ll8//vfp0ePHuTn55c799ixY1x44YWl586YYdVfdnY2AwcO5L777mPQoEFMnDiRkydP2u+VkcGZZ57JueeeywsvvBCS7xBur54fA1Mdj54s4G7sy+YDEfkesBO4saqLdGiZyBt3j+S+N9OZu9F2EAqLPSzZlq+tfkVp4Pxx1no27D4S8Jijp4rZtPcoHgMugQGdWtAisXLXxjO6tOSxqwYFvObJkycZNmxY6fbDDz/MTTfdBEC7du1YsWIFL774Is8++yyvvvoqf/zjH5kwYQIPP/wwn3/+OZMnTz7tmomJiUyfPp2WLVuSn5/PqFGjuPrqqwHIzMzk3Xff5ZVXXmHSpElMmzaN22+/nbvvvpt///vfjBs3jl//+tcBZQ6WsCp+Y8wqYISfXRdW91oiwg8u6MOizHwKSzwY4P3lOew/WsjVw7rqC0BRGjFHTpXgcdKHe4zdDqT4gyGQqef6668HIDU1lY8++giAxYsXM336dAAuvfRS2rQ5XScZY/jd737HwoULcblc7Nq1i7y8PAB69uxZ+qJJTU0lOzubw4cPc+jQIcaNGwfAd7/7XWbPnl2r7wX1JEibl9QebXjnvlEszcpne/5x/pexiylLdvBO2k7eu/9cVf6K0gCpqmUO1sxz26tLKS7xEBfr4p83Dw+rPkhISADsAHBJSQlglXpVTJ06lf3795ORkUFcXBwpKSmlbpjea3qve/LkSYwxYfFuqnchG1J7tOHB8X3p2a45Lud+FLsNf/1iU1A3XlGUhkdqjzZMvXcUv5jYn6n3jopII3D06NF88MEHAMyZM4eDB08fazh8+DAdOnQgLi6O+fPns2OH3+CZpbRu3ZpWrVqxePFiwL44QkG9avH7MqpXEvGxLopLPCDC0qwD3PFaGmentOH8Pu1J7dGGjB0HWZpVoB5AitIISO3RJqT/84o2/ksvvTSgS+djjz3GLbfcwvvvv8+4cePo3LkzLVq0KHfMbbfdxlVXXcWIESMYNmwYAwYMqFKO119/nXvuuYemTZtyySWX1Pj7+CL1oZU8YsQI4y8ef6li79mWaSt38c6ynYAd3JkwoANfb9mP22OIj3Ux9d5RAKUvAt91fSkoSnSxceNGBg4cGGkxqkVhYSExMTHExsayZMkSfvCDH9SpO6i/eyYiGcaY08ZZ622LH8q/4ZduP4BgI8F5DKXePwCnij088HY6B44X4/aUf9HFuoTX7z6bMX3LJ6nR3oKiKNVh586dTJo0CY/HQ3x8PK+88kqkRaqUeq34fRnVK4mEOFfp4M5Dlw7gz7M3UVziweUSjOE0pQ9Q4jF8b0o6Vw7pjMsldGyZwMETxby/PAePxxAX6+IPV57B4ZPF+hJQFKVS+vbty8qVKyMtRlA0GMXvHdzxbaUP6da6nGnHO+of4xIQwe32EBPjon+H5ny0cpff6xaVeHj043UAJMS6eOe+yAwcKYqihIoGo/jh9MGditu+LwYos/EvzSpg/Z4jpZM/bkztxozVuykq8QCU+gcXlnj4eOUuVfyKotRrGpTirwp/LwYvXg+huFgXk87uzqSzu7M0q4A2TeN54pP1FJV4MMZOGuvWpgklHhMS04+OJSiKUtc0KsVfGf7MRN5ygP6dWrA0q4D+nVrw50838pfZmxDsy+Kd+8p7C1VHeS/Zls9try7DY9SMpChK3aGK3yGQD7DvvnW5h/nHvEwM1vTzvTfSOFroxhj/bqOBFPmL87eVMyM9MWs9EwZ0YHTf9voCUJQIYoxhzJgxPPLII1x22WUAfPDBB7z22mt8/vnnAc994403SE9P5/nnn68LUWuEKv5qMqZfe15euI2iEg8uERLiYjh00k7ZPlXs4VcfriL34Mly8wf8KfG9h0+Rln2gdPaxMbA69zCrcw/zwvxtvHv/6eepWUhR6gYR4eWXX+bGG29k/PjxuN1uHnnkkSqVvjd8Q7Sjir+aVDQLAdz2ylKK3B4EIffgSYrdthlf7CSP8aekn/l8EwZ4/paz2F5wnN2HTvJu2k48BorcHp6bs4U3vzfSeiABHyzfyUMfrcVjIMYl3DO6J11bJ3LsVAnn9m6nLwJFyUmD7EWQMgaSR9b6coMHD+aqq67i6aef5vjx49x+++388pe/JCsri6ZNmzJ58mSGDh3K448/zu7du8nOzqZdu3ZMnDix9BqffvopTz75JLNmzaJdu3a1lilUqOKvAad5C91X9iJwezzc+soySjwGDzCwc4vTzl+58yAfrdzFDy/ozeVDOwO2NT9tRa4NQYHwzbZ8rn1hMX06tGDbvmOs2XW49Hy3x/DKwqzS7cS4rRGLT6IoYWf2Q7B3beBjCo9A3jowHhAXdBwMCS0rP77TELis6oxa3rj78fHxjB49muHDh/Pxxx/z1Vdfcccdd5TOzM3IyGDx4sU0adKEN954A4Dp06fz3HPP8dlnn/mN1BlJVPGHgIovgve/fy5vL8nm07V7+N1Ha7nqzAIuHdyZ1B5tSM8+wI/fXUnrpnH8cHyfctco7Un0bMu8TXm8uCCLtbtsHPJLB3ViweZ9FLut59Flgzrz8apdGKyJaXHmflX8SuPl1GGr9MEuTx0OrPiDpFmzZtx00000b96cd999l2nTpgEwYcIECgoKOHzYNsiuvvpqmjRpUnre/PnzSU9PZ86cObRsWXs5Qo0q/jDgfRGM7JXEwx+t5ZVF23l10XYGd23Jht1HcRtDXIywee/RSucdLN1u7f8eAzECQ7q14r6xvcqZmGav30Nhsc1NMMN5CYzRgWGloRFEy5ycNJhyNbiLICYebng1JOYeAJfLhcvl8hv91xsyuVmzZuXKe/XqRVZWFlu2bGHECH8pSSJLvQvLXJ84cLyobPAWyNx3DLfz4/F4TMDk8d7oozECcbGu0gHdB8f3KX1BTL13FL+6pD+TRnQlK/8E/5ibya2vaEpKpRGSPBLunAkTHrHLECl9X8aOHVsaFnnBggW0a9eu0tZ8jx49+Oijj7jjjjtYv359yGWpLdriDyO+oaNtzJ9BPPHJ+tJtb8vdH5XNLah4TGqPNrwwf2tpgLrCEg8zdHax0hhJHhkWhe/l8ccf5+6772bo0KE0bdqUKVOmBDy+f//+TJ06lRtvvJFZs2bRu3fvsMlWXep1WOb6QEUXzHC4ZHqzDxWVePAYiI0RHhjXmyZxMer6qdRL6mNY5kjTaMIy1weqih8Uqjq8vYMzOrfguS+38PxXWwFIjKt8LoGiKI0TVfwNhHKzi3cfYd2uI6UeP/M25qniVxSlFFX8DZDzerfjhbitpaafd9J20rV1Ew5pTgFFUVDF3yDxNf10adWEpz/fyCMfr0OABDX9KPUEY0ypu6QSmOqO1arib6D4mn625x/jX19txRA4jISiRAuJiYkUFBSQlJSkyr8KjDEUFBSQmJgY9Dmq+BsB4/p34IUF23B7DHExgd1IFSUa6NatG7m5uezfvz/SotQLEhMT6datW9DHq+JvBKT2aMNDlw7gT59t5DeX9tfWvhL1xMXF0bNnz0iL0WDRmbuNhNtGdSc+xsWew6ciLYqiKBFGFX8joWl8LGf3bMPCLfmRFkVRlAijir8RMbZvezbnHWWvtvoVpVGjir8RMbZfewAWZuqAmaI0ZsKq+EUkW0TWisgqEUl3ytqKyJcikuksdaSxjhjQqQUdWiSwcIsqfkVpzNRFi3+8MWaYT6Cgh4B5xpi+wDxnW6kDRIQxfduzeGs+bk/0B+dTFCU8RMLUcw3gjWc6Bbg2AjI0Wsb2a8ehE8Ws9UnlqChK4yLcit8Ac0QkQ0Tud8o6GmP2ADjLDmGWQfFhTN/2iKDmHkVpxAQ1gUtE2gLGGFPd1E7nG2N2i0gH4EsR2RTsic6L4n6A7t27V7NapTLaNotnSNdWfLZ2DzEu0aBtitIIqbTFLyLdReQ9EdkPLAOWi8g+pywlmIsbY3Y7y33AdGAkkCcinZ06OgP7Kjl3sjFmhDFmRPv27av1pZTA9O3QnE17j/K3OZu57VVN1agojY1App73scq6kzGmrzGmD9AZ+Bh4r6oLi0gzEWnhXQcmAuuAmcCdzmF3AjNqLL1SIxJiYwCbyN0btE1RlMZDIMXfzhjzvjHG7S0wxriNMe8BwUT56ggsFpHVQBrwqTHmc+Ap4GIRyQQudraVOuSGs7rijXdYVe5fRVEaHoFs/Bki8iLW8ybHKUvGttJXVnVhY0wWcKaf8gLgwuqLqoSK1JS23HJOd95ZtpO/3Xim2vgVpZERqMV/B7AW+CPwBTDHWV8HfDf8oinh5KcX9sUlsH73kUiLoihKHVNpi98YUwS85HyUBkbHlomM79+B/2Xk8ouL+xEbo9E7FKWxEPDfLiKXiMhLIjJTRGY465fWlXBKeLnp7GT2HS1kwWb16VeUxkSlLX4R+QfQD3gTyHWKuwE/EZHLjDE/Db94SjgZP6AD7Zon8N7yHC46o2OkxVEUpY4INLh7uTGmX8VCEXkf2AKo4q/nxMW4uCG1K68u2s6+I6fo0DL4nJ2KokQ3adsLiGme1MnfvkCmnlMiMtJP+dmABnRvINw0Ihm3x/CrD1frRC5FaSAszz7ApP8sJaZZm67+9gdq8d8FvORMwvKaepKBI84+pQFw8EQxLoGFmfmkZR9g6r2j1L1TUeo501fssivif38gr54VwDki0gno6lwi1xizN9RCKpFjaVYBxonQXOTM4lXFryj1m+NFJXbFGL/x16vy6hGgh/NJBno4ZUoDYVSvJOJj7c/AJaKzeBWlAbA29zBndmuF+/ih3f72BwrSNhHIBB4HLgeuwE7gynT2KQ2A1B5teOe+UbRuGseZya21ta8o9Zxt+4+RlX+cG1K74T5W4NdCE8jG/0/gImNMtm+hiPQEPgMGhkxSJaKk9mjDVUO78L+MXApL3KVB3BRFqX/M3ZAHwIUDK3fRDmTqiaVsUNeXXUBcbQRToo8xfdtxstjNih2HIi2Koii14MsNeQzq0pKurZtUekygFv9r2Bj871E+SNvNwH9DJqUSFZzbO4lYl7Aocz/n9lY7v6LURwqOFZKx8yA/mdA34HGVtviNMX8BbsV685wLnOes3+bsUxoQLRLjOKt7GxZl5kdaFEVRasi8TfswBi6uYiZ+wNSLxpiNwMZQCqZEL2P6tuO5uVsoOFZIUvOESIujKEo1mbshj86tEhnUpWXA46odklFEpjjB2gbXWDolKhnTrz3GwDfbNCOXotQ3ThW7WZSZz0UDO1KV131NYvE+D8xFY/I3OIZ0bUXrpnEs2qLROhWlvjHl22xOFrtJade0ymOrrfiNMcuNMdOMMb+tkXS1JScNFv3NLpWQEuMSzu/TjoWZ+yub8KcoShSSseMgz3yxGYC/frG5yrhbgSZwxYrI90XkcxFZIyKrRWS2iDwgIpFx59yxFF6/HL56EqZcrco/DIzt2468I4Vk7jsWaVEURQmSpVn5uD22sVbshF4JRKAW/1vAME6fuXsm8HbtRa0Bi58DTzEYD7iLIHtRRMRoyIzu2x6Av3y2UaN1Kko9oWtra94RIC7WVWXolUBePWcZY/pXKMsFlorIltoIWSPyt0LWfBCXVfzGDa171rkYDZ29h08hwPzN+1mSVaDROhWlHpB/rBCA+8f1YuIZnar8zwZq8R8UkRtFpPQYEXGJyE1A3TcFZ/wQ4prATVPhnAcgoRV89kv4/HdlJp+K9n8dD6g2vl3EYLqMiqJEnq+37KdPh+Y8fNnAoBpqgVr8NwNPAy+KyEFsL6IVMN/ZV3cc2wc5mXDdZBhwuf10TYWP7oOlL8DSF6FpEpwoAIwVtXQbiE2EO2dCsr+8Moov3midhSUeRKN1KkrUc7LIzbLtB7j9nB5BnxNo5m62MeYmY0x77Mzdc40xHZyy7SGQN3iO7obu58HQSWVlh3Os2cdKC64Yuzxt2+h4QDVI7dGGd+49h7bN4unXsbmaeRQlylm2vYCiEg9j+7UL+pyq4vG3ckw7dwLfFZGbRKR17cSsAcbA7hWQu7ysLGUMxCSAxEBsE7jgd3bpux2T4P0i9nglKFJT2nLHuT3YuPcoeUc0y6aiRDMLt+QTH+vinJ7B984DuXPeAawALgCaAs2A8UCGs69ucReXb7Unj7TmmwmP2OWIu07fvusTSOprTT3tB9S5yPWZK4d2wRj4dM2eSIuiKEoAFmbu55yebWkSH3w49UA2/keAVGPMId9CEWkDLAPerImQNUMgJv70VnvyyPJ2e3/b1/8HXpkAaf+Bsb+uG3EbAH06NGdApxZ8smY394xW7ylFiUZ2HTrJ1n3HuPns5GqdF8jUI5QZzX3xUGkK3zDRsnPNB2e7pkK/S+Hb5+HU4dDL1oC56swurNh5iNyDJyItiqIofljohFcZ2699tc4LpPj/BKxwArL9zvm8jDX//KmmgtaI5h1r55FzwUNw6hDM+b26d1aDq4Z2AdTcoyjRysIt++ncKpG+HZpX67xKTT3GmCkiMhO4BOiKbeUvAB42xtSvKZ1dhkP3c2HFFDv4GxOv7p1B0D2pKUO7teL95TmUeAyjeiWpl4+iRAlp2wuYt2kfY/q0qzIaZ0UCevUYYw4aY94zxvzNGPOss14tpS8iMSKyUkQ+cbbbisiXIpLpLOtGk3QaYpfGre6d1WBYcmuy8o/ztzmbue3VpRrGQVGigIwdB7n91TSKSjwszNxf7f9lTcIyIyJrq3H4TymfzOUhYJ4xpi8wz9kOP0NuLPP7j4lT984g8XoKeIzO5FWUaGFpVgHFbg8AHo+p9v8ykDvn9ZV8bgA6BXNxEemGDe72qk/xNcAUZ30KcG21JK4pySPh6hfs+uAb1cwTJBPP6ITL6UUGE/xJUZTw4/s/rMn/MlCL/33gauCqCp8rgcQgr/8P4DdYTyAvHY0xewCcZYdqSVwbht8KA66EjTPgpJosgiG1RxvuH9sLgCeuHqQ2fkWJAjq2TMAAF/RrX6NAioEU/xrgWWPM3RU/wKGqLiwiVwL7jDEZ1ZKo7Pz7RSRdRNL37w9hRqgLHobCI7DkxdBds4Hzg3F9iI9xsXHv0UiLoigK8NWmfQD8/qozatQYC6T4fwYcqWTfdUFc+3zgahHJBt4DJojI20CeiHQGcJb7/J1sjJlsjBlhjBnRvn31fFQD0mkwDLwalr4EJw6E7roNmFZN47hwYAdmrd5NidtT9QmKooSVuRv30bNdM3q3r54bp5dAQdoWGWN2VrIvvaoLG2MeNsZ0M8akYKN5fmWMuR2YiY39g7OcUW2pa8sFD0HRUfjgDvXpD5Jrh3cl/1gRi7bmR1oURWnUHCssYem2Ai4cUHMreY28emrJU8DFIpIJXOxs1y1Fx62HT/YimHKVKv8guKB/e1o1iePjlbsiLYqiNGoWZ+6nyO3hwoEda3yNQLF6QoYxZgF28hfGmALgwrqot1KyF5UFo/D69KuXT0ASYmO4YmhnPlqRy7HCEpon1MlPR1GUCszduI+WibGMSKm5o0VVYZldIjIp0DH1kpQxEBtv1zVkc9BcP7wrp4o9fLFub6RFUZRGidtjmL9pHxf070BcTM0NNlXN3PUAP6rx1aOV5JFw5ywbsrlZe+h2dqQlqhek9mhD+xbx/PurTJ3BqygRYFXOIQqOF3HhwNp5wQfzyvhSRH4lIslOuIW2ItK2VrVGA8kj4Zzvw9G9ULAt0tLUC1bsPMTB48VkF5zg1lc0fIOi1DXzNuYR4xIu6Bd+xX8P8CCwEMhwPlV69dQL+l5sl1u/jKwc9YSlWQV4jB0cKdTwDYpS53yyeg9dWiWydf+xWl2nSsVvjOnp59OrVrVGC21SoF0/yJwTaUnqBd5E7N44gM2qkfFHUZTaMXvtHnYePEHuwZO1DphYpeIXkaYi8qiITHa2+zqzchsGfSdC9jfWxVMJSGqPNky9dxQ/u6gvHVrE825ajk7oUpQ64r3lOYB1SKxtwMRgTD2vA0XAec52LvBkjWuMNvpcBO5C2K5hmoMhtUcbfnpRP564Zgib844ydZnfOX6KooSYXQdPIAIxUvuAicE4Y/c2xtwkIrcAGGNOSnWj/kczPc6DuGbW3NP/0khLU2+4ZFBHzuudxDOfb2L/0ULGD+igAdwUJUzkHDjB1v3Huf2cHnRunVjrpEjBtPiLRKQJzpQnEekNFNa4xmgjNgF6jbMDvMZfimHFHyLCjanJHC9y8/z8rZqkRVHCyOfO3Jn7x/biwfF9at3ICkbxPwZ8DiSLyFRs8pTf1KrWaKPvxXBoJ+RvibQk9Yrdh0+WrmuSFkUJH5+u3cPgri3pntQ0JNcLxqvnS+B64C7gXWCEE4Kh4dDHcev84lGN21MNRvVKIi7GWv1iYjRJi6KEg12HTrIq5xCXD+kcsmsGO+d3HDa+znig4cU3OLoHENg6B6Zcrco/SFJ7tOGl21MBuGlEstr4FSUMzF67B4Ar6lLxi8iLwAPAWmAd8H0ReSFkEkQDvonX3YWaiL0aXDSwI73bNyP34IlIi6IoDZJP1+5hUJeW9EhqFrJrBuPVMw4YbIzxDu5Owb4EGg4pY+wgb8kpu91jdGTlqWec0yuJWatskpbYWgSOUhSlPLsPnWTlzkP8+pL+Ib1uMP/SzUB3n+1kbFrGhoM3aNuAK8B4bJIWJWhG9UriaGEJG/ZUlrBNUZSaMHmhjSPWI0SDul6CUfxJwEYRWSAiC4ANQHsRmSkiM0MqTSRJHgnfeQNaJcP8v6hrZzUY1dPG7FuWpaksFSVULNy8nynf7gDgVx+uDqm7dDCmnj+ErLZoJzYexvwSPvkZbJ1bFsRNCUiHlon0ateMpVkF3De2YYRxUpRIsn73YX783srSfFFed+lQOVBUqfiNMV+HpKb6wrDbYNFz8MUjsGcN9Byj2bmC4JxeSXyyejdujyHG1XAmditKXZKx4yCvLspi7sY8WiTEER/rwu321DpEQ0U0f15FYuNh8A3wzd/hqyftoO+dM1X5V8GoXm15N20nG/ccYXDXVpEWR1HqHRk7DnLTf5ZQ4jG4BJ6dNJRWTeJZmlVQ6xANFVEXDH/Ee92mPGU5eZWAnNPTtkZ09q6i1Iy3lmRT4rHGHQE27jlKao82IQnRUJFqKX4RaSMiQ0MqQTTSaxyIc2ti4jUnbxB0apVISlJTluoAr6JUm6z9x5izPi9k0TerokpTj+PJc7Vz7Cpgv4h8bYz5RdikijTJI+GM62DjDPjudDXzBMmoXkl8tnaP2vkVpRqcKCrhgbczSIhz8febh7F137GQm3YqEkyLv5Ux5gg2Xs/rxphU4KKwSRQt9L8UPCWQ2DLSktQbRvVK4sipEv44a71G6lSUIMjIPsANL35LZt4x/n3LWVwyqFNYTDsVCUbxx4pIZ2AS8ElYpYkmugy3y90rIytHPaKpk4rxrSU7NEyzolRBxo6D3PzKUjbuPUqMS2hSh6lMg1H8TwBfANuMMctFpBeQGV6xooC2vSG+BexeFWlJ6g2Z+2wC6FCkhlOUhs68jXkUu+1grjGmTv8vwfjxfwh86LOdBdwQTqGiApcLOp+pLf5q4A3TXOw2uFyiYZoVJQCrcw8D4KqDwdyKBBOds5+IzBORdc72UBF5NPyiRQFdhkHeOnAXR1qSekFqjza8d98oOrdKpGViHIO66PiIovhj/qZ9fLM1n5vPTuaXE/sz9d5RdRrWPBhTzyvAw0AxgDFmDXBzOIWKGroMtxE792+KtCT1htSUtjw3aRgFx4t4a8mOSIujKFHHN1vz+el7K+nWpglPXDO4TgZzKxKM4m9qjKmYmaQkHMJEHZ2H2WVd2flz0mDR3+p9Iphzeycxtl97XliwlSOntLekKF4ydhzkjtfSOHKqhH1HC1m763BE5AhG8ec7Cda98fi/A+wJq1TRQttekNCybuz8OWkw5SqY9392Wc+V/28u6c+hE8U8MWs9L8zfqh4+igJ8usbGswJwuyPnABFMrJ4HgcnAABHZBWwHbq/qJBFJBBYCCU49/zPGPCYibYH3gRQgG5hkjIlOreAd4N2zKvx1bV9UlgimxAkTUdXEsZw0e1xK9AWSG9y1Fef1TuJ/GbtwCcTHuurcjqko0cb63baFXxezcwMRjFdPFnCRiDQDXMaYYLOUFAITjDHHRCQOWCwis7ETweYZY54SkYeAh4Df1lD+8NNlGCybbAd4Y+LCV0+xT+pCoeowEdvmw9vX2/WY6Awkd0aXlny7rQCPgaIQh5VVlPrG6pxDLNt+kEkjutEjqVnYZ+cGIhivnp+KSEvgBPB3EVkhIhOrOs9Yjjmbcc7HANcAU5zyKcC1NRG8zug8zObh3bcxfHUc3gVpr0DHIdDjfHuX2lYR1/6bf9hsYSZ6A8ldNrgzCbH2J+YxsGDTPv45d4uafZRGyV+/2EzbZvH84apBERnQ9SUYG/89TsiGiUAH4G7gqWAuLiIxIrIK2Ad8aYxZBnQ0xuwBcJYdKjn3fhFJF5H0/fv3B1NdeAj3DF5jYNZPwFMMN70Jlz0NeGDjrMrPOXkIcpeXbbtiojKQXGqPNrxz3yh+ObEfY/u2Y/mOg/x9bia3vqKzepXGxbdb81m8NZ8fXtCb5gmRj4YfjOL3Rtu6HBurZ7VPWUCMMW5jzDCgGzBSRAYHK5gxZrIxZoQxZkT79u2DPS30tO0FCa3K2/l3LIWv/xqaAdh5T9hsX6l327o6DoakPrB+euXnLHsZio7Dlf+A+ObQvj90O7v2soSB1B5t+PGEvpzTKwlv3LbCEg/zN+2LrGCKUkcYY3jmi810bpXI7aN61F3Faz+kawvp5G9XMIo/Q0TmYBX/FyLSAvBUp35jzCFgAXApkOfE/sFZRrcGEIG2PWHzbNixxHrdvH4pzH8SXr8Mtnx5uhum73YgF83012Hxc876a/YYERh0nTXdHPPT0zl5CJa8CAOuhBF3w4Tfw961p5t6osw1dFSvJOJjXaXK/7O1uzl8Ul09lYbPy19vY1XOIa4d1oXEuDqKx5P+Oky7l07Npau/3WKqSCouIi5gGJBljDkkIklAV2ciV6Dz2gPFzjlNgDnA08A4oMBncLetMeY3ga41YsQIk56eHlDOsJGTZhW8pwTb0alwv1xOt83jgZhYOPs+WP4KuEusVxACxn36AOzOpdZt011ktyUGJjxic/7mrYeXzoMrnoOzv1e+vvl/hq+fhgcWQ6chUHwK/jUM2vSEuz+zLw6va6i7KKoGfjN2HGRpVgFxMcJfv9hMz3bNuHxwZ8b0a6+DvkqDZHn2ASa9vAQDJMbVkWebj24ZMfkY6bvdp1logvHq8YjIdqCf46IZLJ2BKSISg+1ZfGCM+URElgAfiMj3gJ3AjdW4Zt2TvcgqdQAM9L4IdnzjKNVYaNEFDm63u91FsPSFsnM9Ph2jkpOw8i1r01/9Dqx+H5q1hxMFjseQT8KXDmdAu37W3OOr+LfOhcX/sAPAnYbYsrhEGP0LmP1rmPlj6H85LP57mWuoO0jX0DogtUeb0h99cYmHv87Zwpa8TF76ehvv3Ff+D+F9SUTS80FRasuUb7PDljD9NHLSYNVUWPVeqW4xHPPbsg8mEcu9wE+xdvpVwChgCTAh0HlOj2C4n/IC4MIqv0S0kDLG5t11F1nlfIHjeer1nweYcqVV3q5YGHk/pE22PQRxWvyeEsDAijdh5dvWEweBy5+FZu1O98X3mnsW/hWO5kGLjrBzGbxzk71Wbrp9yN7jOwy0y5Vv2Q9i6zaeqB34RaS0/1RY4uH1b+zLc2lWAc0TYnjy042UuA2xMcK/bxlO+xaJ5V4EFV8M+qJQoo0St4eMHQcR6iAQW04avHGFY0Eo0y15z4za7e/wYIaXfwqcDSw1xowXkQHAH0MocnSTPNKaSioqZ98W9J2flN9/xjXlXwzZi6DbSFjyAmyZbcvEBfs3woBf+m+ND7rOmnRm/RR6nAtLX3JeINilbys+N40yM5TAuT+C/pfB1O9Yr6QoaO1XZFSvJBLiXBSVeDDAJ2v28NnaPXgqtE+K3YYH3l5Ruh0jwtk925CefRC3x+ASoW/H5mzJs9NLdKKYEi1MW5HLnsOn+O2l/fEYQt8o8U7gTGgFS/7tYzYu0y27jpq9/k4NRvGfMsacEhFEJMEYs0lE+odO+npA8sjAyrPifn/bYHsOWQvKeg+BWuKFRwGxL4ots6FpW3DF2VZ8xXNTxkBsYtl1z7ja1jnyPvj2eTiyB1p2rsk3DxupPdow9d5RLM0q4OyUNvxjbibfbiubvh4j9kUWE+NiaNdWpDvun25jSM8+WJqU2m0MOQdOlL4wwt6dVpQgKCxx8695WzkzuTUPjOuNSIhTkXrH8bwm3fgWlesHPwSj+HNFpDXwMfCliBwE/HYflCqorPfgj3JeOi4454c2Cby/cyu7bupd8M0/rflnXMDx84jga/P/5UQXGa8spcTtIS7WxR+uHMTBE0WlXePbXl1KcUnZvic+WV+6/cgVZ/DHWespLPHgAYYnt47cl1IU4L20HHYdOslTNwwJvdIH+1/3Kn0Ezvsx9B4fdAiXKr16yh0sMg5oBXxujCmqsdDVJKJePZEiJw2mXF3Wiq+pZ86b10J+JvxsjbX3RzGB7PRV2fQzdhzk7aU7mL5yF+f3TuK8PkmM6tVOW/5KnfPt1nzuezOd7klN+ewnY8Kj+LfOhbdvAMT29ivRDyKSYYwZcVp5EO6co4D13hg9jh//Gc4s3DqhUSp+CE0Qtg0z4IM74NYPoN8loZUvCnl0+jreXrYDARLqyn1OURwydhzkpv8socRjiI9x8e79Yfr9LX0JPn/IOpMMubFS/VCZ4g9mAtdLwDGf7eNOmRJukkdav/7aDM72vxyad7QTOhoBnVsnAHaYu0jz/ip1zJwNe8vGnzxh+v0ZY//PXVPh8r/WSD8EFbLB+HQLjDEeghsbUKKBmDgYfjts+Ry+fCxqZvKGi1G92pHoExiuQ/OECEukNCZW7bROCDV23wxmxv3OJZC/2YZ5qSHBKPAsEfkJZa38HwJZNa5RqXs6DwOMHehd9p+omckbDlJ7tGHqfaOYtzGPD9JzePbLzYzt356OLasz91BRqs/8zftYtv0gt53TnS6tm1TffdM3SoArDr47HXr68c5Jf90miBp8fY1lDUbxPwD8C3gU24OeB9xf4xqVuqcg01kxUTWTN1x4vYWuOrML33npW26ZvJRrhnVhdF8NDaGEh1PFbh6fuZ5e7Zvx2FWDiI8NxphSgU2f+szVKYZ3b7a99SZtrcdO8kg4ccCO2511B8Q3q7G8VUpnjNlnjLnZGNPBGNPRGHOrMSa6A6sp5UkZY2MBQVA+vg2FgZ1b8pML+5KVf5y/z83ktlc1HLQSHh6bsZ4dBSe4Y1SPmil9gAOOIUVi7P+0eScbiXfBn+H1yyFrISx4yuYH6XbaeG21qKGESr0ieSSk3mnXb57aoFv7FSnxmNIY4oXFOtirhJ4Plu/k/fQcAJ76fFPNGheHc+04XP/LbbDGuz6F4bc6YV+wPYCpN0Laf+z2rJ/VarxOFX9jof/ldhnbuGzd3tAQYO2UTesqLK7SKFiTe4jfz1hful1cU0+yRc9Zb53Lni7z5EsZY6PrSoxdNvUxU9Yy614wqRf1n9IQ6DjILvPWRVaOOsYbGuJnF/YluU0T/vVVJjkHTlR9oqJUweqcQ9z26jJaNYklIdZV8wTqh3JsAMfht0Pr7mXl3hn5Ex6Buz6BSW/ahpvXFFQLk20wg7tbReR/2OxbG2pckxJZWnS2g0SNTPFD2WDvNcO7cvXzi7nztWVcM6yrDvYqNea95Tv5w4z1tG4Sy/QHR7P38KmaRYfNSbMTsYzHtvQrUjHu152zaj+pk+AU/1DgZuBVJynLa8B7Th5epb4gAp0Gw97Gp/i99GzXjB+N78NfZm/i73MzeWnBNqbepzN7leAxxvDn2Zt4ZaEdiD18soS9h0+VizsVNL6B1iQGju6B1smBz6kqYGSQBOPVc9QY84ox5jzgN8BjwB4RmSIifWotgVJ3dBwC+zaCxx1pSSKG72DvqRIPi7b4SW+pKBXI2HGQp2dv4roXvy1V+mBj7tfYYSB7EZQUlt+uI4JJxBIDXAHcDaQAfwOmAmOAz4B+YZRPCSWdBttMYAXboH3jfGzewd7CEg/G2JjpHmMY17+DtvwVv6RnH+CWV5ZS7LYBDK4e1oU56/eWRoetcXKVlDG2J25MnbtZB2PqyQTmA381xnzrU/4/ERkbHrGUsNBxsF3mrW20it83D8DuQyeYuiyHf321lf8szDot/WNt0axg9Rfvs2vbLJ5/zsssVfougf4dW3DnuSm1f7btBwBiU6le9HidulkHZeM3xhzzt8MY85MQy6OEk/b9bXrIvetg8A2RliZieO2xL8zfWtrgKizxsChzf60UdMaOg3yzdT892jZj1+GTPDdnCx5jNCtYPSNjx0FufWUphSU2Z3aLhBjiYgSPx5S28Gtk069I1gIwbhj/SJ3PrQlG8ZeIyIPAIKDUCdwYc0/YpFLCQ2yCTeLeCD17/DGqVxIJsTb9o8fAgs37+eEFfcrNvJy/eR/Ltx/gwoEdT8sDADZHcJ8Ozflmaz5vLd2Bvyjnp4o9vPT1Nu4+L4VVOYe0BxClZOw4yLyNeczdmFeq9AX43uhejOnXPvS9t8w5Nm1iBCZUBqP43wI2AZcATwC3ARvDKZQSRjoOhh3fRFqKqMDX7HOisIQXFmzj7jeWc26vtrRrkcCna/awKDMfgBcXbKNfx+Zs238cj8fgza1RMUcwWHPAxDM6MX/zPopKPCAwd0MeczfkARDjEq4d3oW2TeO5dHBnfQlEAYu35nPna2m4PWUmHbA5nMf0ax+aFr4vxthkKr3H2wi6dUwwir+PMeZGEbnGGDNFRN4Bvgi3YEqY6DQY1n5ggz01bRtpaSKO7x/6yKkS3lq6g2+2WmUf65LSFPYAuQdPlioG35a9AFed2Zk5G/JKB/zuG9uL+8b2sq3Enm35aOUupi7bCYDbY5iWsQuA177J5u3vjeTc3u3KC7ZhFuzfZNNtNqIQG5Fg7oY8fv7+qnJK/+aR3elakwibwZK3zrpv9r049NcOgmAUf7GzPCQig4G9WO8epT5SOsC7Dnrq2LwvnVollip6Aa4/qyszV+8uVeaPXnFGaa7fGJeACG4nR/Cd5/XkzvN6nmYOKFUaIkxbkUtxiQcRwe0xGOxL4MF3VvDrSwZw4LjNMZy67yPMp78ABLPwWVx3zar3yj8aB7rnbsjjmS82sSXvGN3bNKXQ7Sl9njec1S28cmbOscs+F4WvjgAEo/gni0gb4PfATKA58IewSqWEj05D7HKvKv6KeF09vYr+prO7c9PZ3csprP6dWpxm4/er6Cvga1Zq0zTe5wXiwu02PPzRWgCaSiFLEx+jJQAGd0kRe1bNoWs9VvxzN+Zx/5vpGAMJsa6ITZrzvnx6t2/GjFW7mb1uL2B7dn+9cSixMa66ezllzoXOZ0KLTuGtpxKqVPzGmFed1a+BXuEVRwk7zTtAsw46wOsHX+VcmTKvaOutjoLwPdf3BbJwyz7+NW8rBvhVzPu0NEfxGEEwFBPLG7u60Xp+Zr1MHr9131F+8f6q0rGQUyUenvl8I2P6tufc3nX3fSp66vjmPzfGkL7jIA+O7xMeeSrmzj55EHKWweifh76uIKlU8YvILwKdaIx5LvTiKHVCp8Gwd22kpYhKQj6IF2Q9/1mYxXD3eu6J/ZzNPW5hTlYRP46ZxqNF9zAtuz1kbyHGlcm9o3vSq30z9h4+VRprKBrNKBk7DvK/jFxmrtpFbIwQH+uixG29p5ZtP8iy7QeJi8nkvftGkZoS/rGmt5Zkl/PUuebMLnweiklYVZGTBm9cAe5iG2Dtzpk2BLNxQ9+J4akzCAK1+Fs4y/7A2VgzD8BVwMJwCqWEmY6DYOlL8PVf7eAhhCTwk1KBii29Skjt0YYZl3voMe+fFMV3pP9tf6Nw6074YBrju8fwUTal4wH/WZjFWbKFUa6N/OmrgSSkjCJ9x0HcnuiZL5Cx4yA3T15CsduGx3jxtrPo0DLRTpo7eJJ30nZigGK34QdTV3DP6J64PZ6w9GiWZuXzzOebWbHzECJW6cfHuvjuuSl8NxSTsKpi+0IbQhnsrPlv/w1Fx+1LwHjCU2cQVKr4jTF/BBCROcBZxpijzvbjwId1Ip0SHmKb2hRv85+E+X8qmzbubZGo8q89pflTPXb+RKD7mpNG/y+/a2MomULIW8/QM0ZC296MjttMQtyY0pbp93sW8OCOJ4mlhELiuT37dxR77CzsU8Ue/jUvkxtSu5Jz4ASjellPobruDUzLyC030zUr/ziXDelc2juZttIOcrtcQonbw1OzNwEQ48rkN5f054wuLVmTe7hUZn9zJ4L5Pp+u2c2P3l1pIyK4hMeuGsjRU+6gxmRCRhPv9R3b0saZZfveui5i/7dgBne7A0U+20WoV0/9xlPss2HKfBPdhQ0+H2+dseaDsvypVeU53jCzLHCep6Ts2JTRtF7/MVPv+Q9Ls+3Er06rnyduRwkiEG+KeXTQAW7aaAekEfh6y36+Lg08twWX805PiKub3kB2/nFmrd6FYJV+RTNKxXGUb7fl89ycLaU9mr84LwGc8ycM6MDXW/ZT4jb2u2DnTiTEugKG2Ji1eje//GB1mdutMRw95ebB8XUcV3L/ZhuHZ/QvoPcESP8vrHnf7otg/utgJ3Clich0Z/taYErYJFLCT//LrKnHXQSuGECs0jcGeoyOtHQNg8O55bcDBeDK32KXFRNs9BwLK6aQmpBD6vjhtmxXW3AUoEsMw0eO493RqaWKdM76vUxemFU698A7qGrTTuaHVfEfPlHMPW8sJzbGxfO3DiW74ITflnnF8Y0X4rbaHk2Mi/P7tuOrjftKFfzcjWXpvd0+cycKSzz8Y+4WHhzfm4wdZbOhv92az1OzN7Fm12H6dWzOjoITlLjDbMevDI/bJkbvdwmMf9iWidgXvbsoovmvg/Hq+ZOIzMZG4wS42xizsqrzRCQZeBPoBHiAycaYf4pIW+B9bK8hG5hkjNEM2HWJN7OP1/4M8O2/YOMsOHUooqI1CE4dsbbdvpfYP3jWfDiy2/+xB7LsDM4zroPOQ8qPB6Q4L+HsxdDFUfw7v4WEVsjQGyH9dVjzHqk3XFROkU5Zkl021wCh2O3BAFvyjmGMQXxdWmqBrwnm6Kli/vDxenYfPsk7941iZM/gBmwr9gAAvtmaX2raeujSAfx59iZK3B5inbkTJW5rG1+Umc9iZ2a1S4QeSU3IyrfZ1WJdwpPXDibGVYcumhXZuRSO7YVB15WVVfzvRah3HUyLH2PMCmCFiNwfjNJ3KAF+aYxZISItgAwR+RK4C5hnjHlKRB4CHgJ+WwPZldpQMaHDd16H50fAgr9Yb4MQKYdGydoPofg4jPsNdB4Gr14In/7CRmFs3r78sQuftVP2L3vqdJ/uFp0gqS9sXwTn/RiO7IHNs+G8H8HFT0DzjnaM5oxrYOBVgH9FujQrn3W7jjBj1W4SYl30SGrmVxFW9A6qbLtXu2as23WYl7/Owl0hOFFcjDgvnOCp2AOo6FI7pFtrvzb+d9N28r8M27NyG8P+o2UWaWMMy7PD6KIZDOun23GzvpeULw9RMpXaEJTi9+EBYHIwBxpj9gB7nPWjIrIR6ApcA1zgHDYFWIAq/sgTEwdjfw0zHoQtn1tzkFJ9jIGM123Sm66p9gV67UsweRx8eJe18/Z0WnoF22D1e3DOA5VP5EkZDeumgbsEVr5l3QDPutPuG/1z2PQJzPgR7F5lTQrJI/3ONfB4DPdOWc4H6bkIp9v8M3Yc5DbHzz02RrhsSGc+W7MHt8fgcgnDk1uzYudBv7GJfPF4DEuzCmqlbP3JX9nciU/WlM2sfvjygaUT4yJi2vHFa+bpOxESmkdOjkqoMgNXBWrUDBSRFGA4sAzo6LwUvC+HDpWcc7+IpItI+v79miWpThh6E7RJgTmPwsK/Wc8UpXrsWmHnSIy4q6zX1PEMm0h7x2L46gnr153+hn0RuGLg/J9Wfr2U0VB4BPasgowp0OsCSOpt98XEwagfWvPcomdhytWVPjOXSzjLUZgG6wE0b2Ne6f4P0ndyqsRT6mY5c9VuSnzCSqzOPVSq9AW4YkgnEuNsgvF4x0+/xsnGa4i3d/OLif2Zeu8obj2ne7ntiLq17vgWju8rb+aJIqrb4r+quhWISHNgGvAzY8yRYO2LxpjJOL2LESNGVNHOUEJCTBwM/o5VIl89WbUbonI6Ga9BXDMYMql8ecuu4I0E5C6CT35m112xcGgHtOjo/3reMZi5j8ORXLjkT+X3H9lVdt2SwF5Z5/ZuR2LcVgqLrYJ/e+kOmsbHMGdDHmtyDyPYd1VcjItfXNyP577cUjow+ocrB5VrTd8zuhf3jO5VIzfLUFJV7yBibPgYYpvYXlgUEmjm7k+A6caYHG+ZMSa3suMruUYcVulPNcZ85BTniUhnY8weEekM7Kv8CkqdE+tNueCJqLtZvWTbV7D6fRt4K7Fl+X09x9p76y7CutA6k3eMCXyPW3SEVsn2mMQ2MOCK8vtTxtjrlpwEPGVB+Pzga//v1roJz325hWfnWI+iGJfw+ysGcryozM99RErbSuMUVRbOQgF2LIFV70K3syG+WaSl8UsgU8//ActEZJGI/FBE2gc49jTENu3/C2ysEN5hJuAYKbkTmFGd6yphptc461YIEXU3q3fkpMHUSXaOxLavTje5eL05JjwCV/zdKuuK7puVXdfrEVR0FHZX8K3wXvecH9jew8q38ZsNxiG1RxseHN+Ha4Z35fqzupbZbo3heJG73GCo99jKthU/5KTBm1fbwf2cpVFrLg1k6skCUoGLgJuAP4pIBvAu8JF3Jm8Azge+C6wVkVVO2e+Ap4APROR7wE7gxpqLr4Sc5JEw7FY7kHjbh9raD5bMOWUT43wnYfni683R8YzgXPqyF5Up8sp6B97rtuhoTUJfP23NdlVce3Tf9rz09bboGAyt7+SkwbqPrGOEN0SDxx21PeZAit8YYzzAHGCOY7a5DLgFeBYI2AMwxiym8sHgC2sgq1JX9B5vFX8TbdkFza4Vdimu4HpKwbr0pYyxYy3BTPg598ew6j3rkisuiAk8RlNZNFKlmmxfBG9eYz2uoF70mAMp/nJK2xhTjDXTzBSRJmGVSoksSX3tMn+LjeSpBGbjJ7BtHgz7LiT1DO3EnOpM+ImJhT4TIH+THUMIYowmagZD6zMZr5dX+ql32HGZKA56GEjx31TZDmPMyTDIokQLSU48k4KtkZWjPnDiAHzyc5vg5srnIDY+9HVUZ8LPoOsgbbI1N7lcUdvibFAUOlZv75jNmbdGrcL3Eig655a6FESJIuKbQqvuZTFkFP/kpFmlf6IAvvtReJR+dUkeCXd+AtPutf79rZIjLVHDxuOxk+d6jrOOEVHcyveluhO4lMZCuz6q+AOx6VN47VKbyUwEiqOoE9zjXLhjhm31f/KzgF4+Si3Zu8ZO1DrzFhjzy3qh9EEVv1IZ7fpB/lZVGhXZ9Cm8ehG8f3uZXdfrbRNNtOsDFz5mvUw+vCtq3QrrPZlf2mWEkqbXFFX8in/a9bW+yJVFlWyMbPkC3rsVcpfbuAcx8cH54keKLmdZ754NH8MbV6ryDweZc2zk1IrB96Kc6oZsUBoLXs+egkxo1TWyskQLK98qWxeB4bdFt/fGzm/K1t2FNgFINMiZkwbrP4YOA20a0Lz1cHA79Ls0OuQLlhMHYFe6DW5Yz1DFr/innU3nR36mDQzmJcg8sg2SE07aiPrivZEyxvryu4use+eqd6HTmXBif/ifX2W/k5w0eP3yClngHL79N9z1aXTfU1+2fWXva5+LIy1JtVHFr/inRSeIb1F+gDcnDaZc5UwmamQB3NzFdiCv32WQfHb9ePH5zgFI6gezfwWzfhzU5K5akZMGU66096xiPRtnlSl9cdnYQnvXURobasWb0X9fvWR+CU3aQtezIi1JtVHFr/hHxPHsySwry15kI0BiGl9+3pxlNjzy8Nth4JWRliZ4fOcA7EqHb/4R9OSuauFt4TfvCN/8y/mdUP534vFA1te23NtrGvE9+Pyhsl7JmvftpMGi49H9cvV4YOuXdlDXFRNpaaqNKn6lctr1g2wfO3HKGPtCMMa21qJxQDNcZM4BV5z11a6vDLgCljxv3Ty9sXxCgbcnWHLKbscm2oBxnhKrzONb2PL0/8Le1XD+zyGxRZli98Yt6jgEZv8GZv8WcEV3WPDdK+38jb71z8wDqviVQLTra1tgRcdteNkWne0fWVw2RHDXEZGWsO7InGv94xNaRFqSmpM8Eq76N8z4AZz7o9Ap1A0zypQ+Auf/xNq9t3wOa/8HX/4Bju6xNvyuI+Cix8qn9vTtleQss/kgIhEWvDrjVxmv22Vi/Qx3oe6cSuWUevY4oRs2OBG0x/7aDhDuSo+MXHXN4VzYt75eDuKdxrBboHmn0IXj8KaPBNsgiE209yl5JFz4B7h3no0auvg5a9vPW2fdYSuj3yW2ZwW211BXvcodS+D1y2De/wXMZAbAhlk2/DXAB9+tl26yqviVyvH17AGbPLrTUDj3QWufXT89crLVJVvn2mXfiZGVIxSIQN+LYNt8m8fXl5w0WFSNlJvrPoL/jLP2/GtfhgmPnm6aad7eZnXz4i4OPNkteSTc9j/7++o5tu5a+6umWtMUxia1+epJWPhs2b3ISYM5f4B3b4UP77DHQVmvpJ6hph6lctr2sq24/C1waKdt4V/4GCS2soNaG2bAxD/ZYGANmcwvbeyi9v0jLUlo6DvRtlhz06DHebbM64lTUmhb7XfOCqx0M+fC/+626zEJNg/wsFv8H9vvEljyQnChpQF6X2AH0Ve9AycP1k148BMFzooL8MD2r+0HsSbOo3soVfY9zoddGY7XUpRO3quCBv6PVWpFXCK07m5b/F4zz6BrneV1Nt9roG57faRiq7ekCLIW2FZykPmio55eF1gzSuacsrJSjy2svX57Fa3YZS+XrXsTz1SGb/axYAdrU++2cnjNSOHEXQI7l0Kv8XDho7bu0qj0pixdJlhvpD4X2hdjdb5PlKEtfiUw7fpZxX9oB3QeZnsBYGdZxiRYc0/3cyIqYsjwnVzkiisLdFZ0rGGYebwktoLkUbbVftHjtqxZh/LHFB2r/HxjYN8G2xtEQpt4xkvnodA1FdJfh3MeCO9LN3shnDwAZ99rXXVz0uwLx9tDmfD7MpdT73et7veJMlTxK4Fp18/OUPSUwEV/LCtPbGld2TZ8DJf8uWGYe8qlTyyGd2+G1j1sKy+uaWRlCzV9L4a5j9lYTC27wKZPrNvlqB9Yb5zlr8LZ34NW3U4/N3uR7e2N/Y3tFYbL3z71bpj5I9i5pMwkVRXZ39jWe89qyLR+OsQ3Lwu05i/5TbCpMusJDeDfqoSVpD7OoBfQtmf5fYOus7bPT38eGs+G6g4uhvo6h3fZpTd9YpM2kLfWRuF856Z66b1RKV7/861zrb16y+cw+mfWfDHpTZsv9oM7YaGf+5j+OiS2hjG/CG8o4sHXQ0JL+PqvZc8z0LNd8Sa8cTl8FYRnjhd3sZ1N3P8y+xLzkjyy/HeruF3P0Ra/EhiPu2z9o+/bgS7vj79pO7vMmAKr36+dvbN0mn+JVbo1vVZOmv3TuwurF5bgxAHYOBNSxtrBxZQx1s49/8nwzHSNNB3OgJZdbS9n4ywbeuCc79t9bXvCiLvtZK9dGc5gr3Mfj+23x599L8SFOQNrfDM7HrFxph1ncTmmJeM+/dnmZsCnv3JONME/r+1f2wHkQdeF73tEIar4lcAc21u2XvHPtDvD2VGNP1plZM7xmeZfi2tlL3ImE1UzrMSS5+1EtcufsVEjvSwMMtF5fUPEmja8bowj7ik/Oa3Uk6bCfVw11ZrBRtxdN3K27OKseGyYBC8lJ23Pwxg7yXD1u9C0DRwvcMZogpxZvv5ja+LqfWE4pI9a1NSjBKbvRNvi8xd33hv9EawiqY1iTGhVtl4bJdtxMKUeGMYEN7v4eAEs+49t9fkq/Zp4o9Qn2vQsM+Oteqe8aaTnWPvcwfZ4YptYxZvxBnQ/r+5cWwffYEM3iMsOuLvisGpLYPU78NolNhRE8Um4/G827WRia2g/sOrn5TXzDLi8vJmnEaAtfiUwySOt65q/ga3kkXDXJzD9ATi+35oPasqpQ2Xrk96quZLdvhAQ+2fe9Cls/ixwfJ2cNJj7uPViGffb0/fXc++NgJSGWaBsYpWvTfvOWTb5zLqP7D3as8rGzR96U93J6M0h7P39gSPnKFj6oh2UBme+yWYYcyWc/1OY90fYvwXa96v82ssm299dx0Fh/xrRhrb4laoJNLCVPBKun2wjVy5/peZ15KbZgTyAUwdrdo2jebD8v1Yx3fwOjPy+9TfPXuz/eG9wsR3f2B5N4ZGa1Vtf6XOhbclXlkUseSRc+Hu4b55NxrPmfVv+zT/qdqDb9/fnXU853yp4f/IPv93OU8h4o/Jr5qTBl7+36/P/3LAG7oNAFb9Se7qNsPFZvvkXFB6t/vkeN+xaAUNutAPGvhOLqsM3/7T2+HG/sdsXPWbNGdO+B/P/cvqf23fSElGYNzfcBGvKatrWmly8VBV2oa6oTP7mHWDAldYUVHzK/7lrPijLmRwt36cOUcWvhIYLHraTYNImV//cfRusqaX7KOtmuHVeeW+iYDi619p6h95kwweA9Qo578d239dPn+7i5w0zDXasoiEN3gZLsG6KfScG7h1EisrkH3G39dbxzjj3paSwLP5StH2fOkJt/Epo6JYKfS+BRX+3f6w+FwVvG/cq425nW1vt6ndtDyD57ODrn/2QrbdfhRm2pWMHfjyPWna1A5e9xsP43zVcW34o8DepKZpJGWtnmX/7LziSW17mr5+xYxUX/5/1AKoP3yfEqOJXQsfAqyDzC9u6/uYfdlAumD9UTho0aw9tUmw4AXFZc08wir/4FMx7AjY4kUKn/8AqdG+9pXlnC09PHuNtDV7+rM02pgSmPg10u1zQe4KdgZy33knqMsva/hf/3eZLPv8nkZYyYqipRwkdx/dRGtyqpBA+/QWsnVb1LNrcNEg+x5pdmraFbiNtWrtAbF8EUyfBcwNh6Qtl5RXD5Ho9j1p2sZPPuvm8TNZPh05DVOk3VJomOSvGejDN/i28MwmatIZL/xJJySKOKn4ldKSMKfP5d8XCvk0w7R6Y92TlU+iP58OBrPIKue9FNrXdsX2nH1980v6Bp1xpexenDsGoB6v2Trngd3B4p437Aja5Sm5ao5ux2ajoc5Hzu3D8/nevsG7HhUdtqPFGTNgUv4i8JiL7RGSdT1lbEflSRDKdZf3MW6b4x9fL4u7ZPl1pT+UJK7xhnX1NCN5ImN4BOICshfD2jfC3geVDAiPQLKlq7xRv3Jd0J2We18xzxrU1+KJKvaD09/ioHewVR9153I3Oi6ci4WzxvwFcWqHsIWCeMaYvMM/ZVhoSvl4W/S71SaMX499zIifN9g66DC8r6zTUxo759nnI+ho++zW8eRVsnQOFh+HcH5/ewq/KOyW+mfX42TDDxuXxZhPzegApDRPv7+LMW+xYTyP14qlI2AZ3jTELRSSlQvE1wAXO+hRgAeBnuqTSIEgeaWPav3sLtOhU3pzjJSfN2tl9A37lLodTh6176JtXVzhBbEyWmniYjLjbTjL7+hlbx4WP1firKfWM+uaVFGbq2sbf0RizB8BZdqjsQBG5X0TSRSR9//79dSagEmJSzoeJT9jp9JkVBmzdJdbumlwhkUv2Ihtnx0v/K6rfwvdHx0F24HjZS3a7YphppWHTwEIr14aoHdw1xkw2xowwxoxo3759pMVRasOw22wKxwV/Ka/Q89ZB8YnTewIpY5zAXDFW4Y/+WeiCpfnG7Zn+g0Y3VV9RoO79+PNEpLMxZo+IdAb8uG0oDY6YOBj7a5j5Y/j4BzYEcPJIWPuhsz++/PGVdctD0VLzjjlAw4uxryhBUteKfyZwJ/CUs/Qzn1ppkCT1wYbSfdfGSeky3Cb5APjofjsGUDHyZzgUcu/xdgJPQ4yxryhBEjbFLyLvYgdy24lILvAYVuF/ICLfA3YCN4arfiXK2LkEO7nL2OBYeesojZtfly1vHeRTlLB69dxSya7GlepGsXjt9t6W9qVPwecPRablXZ9CDyhKGNBYPUrd4K+l3fEMbXkrSgRQxa/UHRVb2tryVpSIELXunIqiKEp4UMWvKIrSyFDFryiK0shQxa8oitLIUMWvKIrSyFDFryiK0sgQ4xs0K0oRkaPAZp+iVsBhP+vV3Q7XsfWhnuqe2w7Ir2MZG1o9kZLR99k1hO/T0GSsbj3V+S/2N8a0oCLGmKj/AOkVtif7W6/udriOrQ/11ODc9LqWsaHVE0EZK/3/1NPv06BkrEE91fkvlnv23k99NfXMqmS9utvhOrY+1FPdc2t63erWo/c89PVUpL5/n4YmY7ieZaXUF1NPujFmRKTlaMzoM6i/6LNrWFTneVZ2bH1p8U+OtACKPoN6jD67hkV1nqffY+tFi19RFEUJHfWlxV+KiFwqIptFZKuIPOSUDRORpSKyysnTW6vIXyLymojsE5F1PmWPi8gup45VInJ5CL5LsojMF5GNIrJeRH7qlP9VRDaJyBoRmS4ircNUz5kiskRE1orILBFpGYLv5O/5hPTeVfJ8QnrPAtQTjnvm9/k4+37s3M/1IvJMCOry93z+z7lvq0Rkjoh0CUE9/u5dWxH5UkQynWWbMNXzvs9vLVtEVtWyjkqfj7P/VyJiRKRdLetJFJE0EVnt1PNHp/xGZ9sjIqEx2fkb8Y3WDxADbAN6AfHAauAMYA5wmXPM5cCCWtYzFjgLWOdT9jjwqxB/n87AWc56C2CL830mArFO+dPA02GqZzkwzim/B/i/MD2fkN67Sp5PSO9ZgHpCes+qeD7jgblAgrOvQ5ieT0ufY34CvByme/cM8JCz/lC4nlGF/X8D/hCO5+NsJwNfADuAdrWsR4DmznocsAwYBQwE+gMLgBG1vWfG1D+vnpHAVmNMljGmCHgPuAabysnb8moF7K5NJcaYhcCB2lwjyHr2GGNWOOtHgY1AV2PMHGNMiXPYUqBbOOrB/pgWOod9CdxQm3qo/PmEFH/PJ9T3rLJ6CP09C/R8fgA8ZYwpdPbVNke13+djjDnic0wzSlOj1ZxK7t01wBRnfQpwbZjqAUBEBJgEvFvLOip7PgB/B35DaO6ZMcYcczbjnI8xxmw0xmwOcGq1qW+KvyuQ47Od65T9DPiriOQAzwIPh6n+Hzld4tdC0U31RURSgOHYt7wv9wCzw1TPOuBqZ9eN2NZLbajs+UAY750fQnrPKhDqe1aOCs+nHzBGRJaJyNcicnYtL1/p8xGRPzn/n9uAP9SynsroaIzZA1aZAh3CVI+XMUCeMSYzVBf0fT4icjWwyxizOoTXj3FMU/uAL40xFfVBSKhvil/8lBlsy+jnxphk4OfAf8NQ90tAb2AYsAfbhQwJItIcmAb8zLf1JSKPACXA1DDVcw/woIhkYLuwRbWtwk+ZIYz37jQBQnzP/BDqe1aKn+cTC7TBdvd/jc1X7e8eB12FnzIDYIx5xPn/TAV+VIs6oolbqGVr3xff54P9jT1CiF+Sxhi3MWYYtsc6UkQGh/L6Xuqb4s+lfAurG9ascyfwkVP2IbZLG1KMMXnOQ/EAr4SqDhGJw/6YphpjPvIpvxO4ErjNOEa/UNdjjNlkjJlojEnF/kG21bIav88nXPeuIqG+Z/4Iwz0DKv0d5AIfOSaANMCDna5fUyr7//jyDiEwX1VCnoh0BnCWtTVdVYqIxALXA++H6HoVn09voCewWkSysfdyhYh0CkV9xphDWJv+paG4XkXqm+JfDvQVkZ4iEg/cDMzE/njHOcdMAELWtfPi/cE6XIft8tf2moLtnWw0xjznU34p8FvgamPMiTDW08FZuoBHgZdrWZXf5xOOe1eRUN+zAPWE+p5V+nyAj7G/Z0SkH3ZANv+0CwRPZc+nr88xVwObalFHIGZiG2k4yxlhqgfgImCTMSa3thfy93yMMWuNMR2MMSnGmBTsS/UsY8zeWtTTXhxvNBFp4v0OtZXfL6EYIa7LD9ZrZwu2pfWIUzYayMB6KSwDUmtZx7tYk0Qx9oF+D3gLWAuswf6AO4fgu4zGdrXXAKucz+XAVqwt1ltWKy+LAPX81LmXW4CncOZ1hOH5hPTeVfJ8QnrPAtQTjntW2fOJB97GvihXABPC9HymOXWswU757xqme5cEzMM2zOYBbcNRj1P+BvBAba8f6PlUOCab2nv1DAVWOvWsw/FGwjaWcoFCIA/4orbfSSdwKYqiNDLqm6lHURRFqSWq+BVFURoZqvgVRVEaGar4FUVRGhmq+BVFURoZqvgVRVEaGar4FUVRGhmq+BVFURoZqvgVRVEaGar4FUVRGhmq+BVFURoZqvgVRVEaGar4FUVRGhmq+BVFURoZqvgVRVEaGar4FUVRGhmq+JVSROQ6ETEiMiDSsijB4Tyvt3y2Y0Vkv4h8Ekm5lNohIsfCeX1V/IovtwCLsblYg0ZEYsIjjhIEx4HBTo5WgIuBXRGUR6kHqOJXABCR5sD52NyoNztlF4jIQhGZLiIbRORlJ8k4InJMRJ4QkWXAuZGTXAFmA1c467dg89ACICIjReRbEVnpLPs75YtEZJjPcd+IyNC6FFoJjPP/+8Rn+3kRuctZzxaRP4rIChFZW91euip+xcu1wOfGmC3AARE5yykfCfwSGAL0Bq53ypsB64wx5xhjFte1sEo53gNuFpFEbMLuZT77NgFjjTHDgT8Af3bKXwXuAhCRfkCCMWZNnUmshIJ8Y8xZwEvAr6pzoip+xcstWAWCs7zFWU8zxmQZY9zYluRop9wNTKtbERV/OAo7BfvMPquwuxXwoYisA/4ODHLKPwSuFJE44B7gjToRVgklHznLDOzzD5rYkIui1DtEJAmYgLUVGyAGMFglYioc7t0+5bwMlOhgJvAscAGQ5FP+f8B8Y8x1IpICLAAwxpwQkS+Ba4BJwIi6FFYJihLKN84TK+wvdJZuqqnLtcWvAHwHeNMY08MYk2KMSQa2Y1v3I0Wkp2Pbvwk7+KtEH68BTxhj1lYob0XZYO9dFfa9CvwLWG6MORBe8ZQasAM4Q0QSRKQVcGGoLqyKXwFrIpheoWwacCuwBHgKWId9GVQ8TokCjDG5xph/+tn1DPAXEfkG25PzPScDOAK8XgciKkEiIrFAoTEmB/gAWANMBVaGrA5jKvbkFcUiIhcAvzLGXBlhUZQwICJdsKafAcYYT4TFURxE5EzgFWPMyHDVoS1+RWmEiMgdWO+fR1TpRw8i8gDWieLRsNajLX5FUZTGhbb4Gykikiwi80Vko4isF5GfOuVtReRLEcl0lm2c8otFJMOZLJIhIhN8rvUnEckJ9zRzRVFCg7b4Gyki0hnobIxZISItsL7A12I9Pw4YY54SkYeANsaY34rIcCDPGLNbRAYDXxhjujrXGoX1QMg0xjSPxPdRFCV4VPErAIjIDOB553OBMWaP83JYYIzpX+FYAfKBLsaYQp/yY6r4FSX6UVOPgjOxZzh2sK+jMWYPgLPs4OeUG4CVvkpfUZT6g87cbeQ4wdmmAT8zxhyxjfmAxw8CngYm1oF4iqKEAW3xN2KcOC3TgKnGGG/cjzzHxOMdB9jnc3w37ASuO4wx2+paXkVRQoMq/kaKY6f/L7DRGPOcz66ZwJ3O+p3ADOf41sCnwMPGmG/qUFRFUUKMDu42UkRkNLAIWAt4J/D8Dmvn/wDoDuwEbjTGHBCRR4GHgUyfy0w0xuwTkWew4R26ALuBV40xj9fJF1EUpdqo4lcURWlkqKlHURSlkaGKX1EUpZGhil9RFKWRoYpfURSlkaGKX1EUpZGhil9RgkBEHheRXwXYf62InFGXMilKTVHFryih4VpAFb9SL1A/fkWpBBF5BLgDyAH2Y0NXHwbuB+KBrcB3gWHAJ86+w9ggdgAvAO2BE8B9xphNdSi+olSKKn5F8YOIpAJvAOdggxmuAF4GXjfGFDjHPInNUfBvEXkD+MQY8z9n3zzgAWNMpoicA/zFGDPh9JoUpe7R6JyK4p8xwHRjzAkAEZnplA92FH5roDnwRcUTnYin5wEf+kQ7TQi3wIoSLKr4FaVy/HWH3wCuNcasFpG7gAv8HOMCDhljhoVNMkWpBTq4qyj+WQhcJyJNnNSUVznlLYA9Tkjr23yOP+rswxhzBNguIjeCjYQqImfWneiKEhi18StKJfgM7u4AcoENwHHgN07ZWqCFMeYuETkfeAUoBL6DjXj6EtAZiAPeM8Y8UedfQlH8oIpfURSlkaGmHkVRlEaGKn5FUZRGhip+RVGURoYqfkVRlEaGKn5FUZRGhip+RVGURoYqfkVRlEaGKn5FUZRGxv8Dgh+tV4i31nsAAAAASUVORK5CYII=\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "df['Days'] = range(0, len(df))\ndf",
"execution_count": 7,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 7,
"data": {
"text/plain": " England York Days\ndate \n2021-03-06 59.804259 42.256597 0\n2021-03-07 59.049200 43.680977 1\n2021-03-08 58.470025 42.731390 2\n2021-03-09 57.636794 42.731390 3\n2021-03-10 57.478676 42.731390 4\n... ... ... ...\n2021-05-31 34.496444 26.113628 86\n2021-06-01 37.378106 25.638834 87\n2021-06-02 40.554685 28.012800 88\n2021-06-03 43.987097 28.487594 89\n2021-06-04 47.110378 29.437180 90\n\n[91 rows x 3 columns]",
"text/html": "<div>\n<style scoped>\n .dataframe tbody tr th:only-of-type {\n vertical-align: middle;\n }\n\n .dataframe tbody tr th {\n vertical-align: top;\n }\n\n .dataframe thead th {\n text-align: right;\n }\n</style>\n<table border=\"1\" class=\"dataframe\">\n <thead>\n <tr style=\"text-align: right;\">\n <th></th>\n <th>England</th>\n <th>York</th>\n <th>Days</th>\n </tr>\n <tr>\n <th>date</th>\n <th></th>\n <th></th>\n <th></th>\n </tr>\n </thead>\n <tbody>\n <tr>\n <th>2021-03-06</th>\n <td>59.804259</td>\n <td>42.256597</td>\n <td>0</td>\n </tr>\n <tr>\n <th>2021-03-07</th>\n <td>59.049200</td>\n <td>43.680977</td>\n <td>1</td>\n </tr>\n <tr>\n <th>2021-03-08</th>\n <td>58.470025</td>\n <td>42.731390</td>\n <td>2</td>\n </tr>\n <tr>\n <th>2021-03-09</th>\n <td>57.636794</td>\n <td>42.731390</td>\n <td>3</td>\n </tr>\n <tr>\n <th>2021-03-10</th>\n <td>57.478676</td>\n <td>42.731390</td>\n <td>4</td>\n </tr>\n <tr>\n <th>...</th>\n <td>...</td>\n <td>...</td>\n <td>...</td>\n </tr>\n <tr>\n <th>2021-05-31</th>\n <td>34.496444</td>\n <td>26.113628</td>\n <td>86</td>\n </tr>\n <tr>\n <th>2021-06-01</th>\n <td>37.378106</td>\n <td>25.638834</td>\n <td>87</td>\n </tr>\n <tr>\n <th>2021-06-02</th>\n <td>40.554685</td>\n <td>28.012800</td>\n <td>88</td>\n </tr>\n <tr>\n <th>2021-06-03</th>\n <td>43.987097</td>\n <td>28.487594</td>\n <td>89</td>\n </tr>\n <tr>\n <th>2021-06-04</th>\n <td>47.110378</td>\n <td>29.437180</td>\n <td>90</td>\n </tr>\n </tbody>\n</table>\n<p>91 rows × 3 columns</p>\n</div>"
},
"metadata": {}
}
]
},
{
"metadata": {
"scrolled": true,
"trusted": true
},
"cell_type": "code",
"source": "df_window = df[-7:]\n# Model as: cases = a * b ^ day\n# : ln cases = ln a + day * ln b\nlog_params = np.polyfit(df_window['Days'], np.log(df_window[['England', 'York']]), 1)\n# Model as: 2 * a * b ^ 0 = a * b ^ d\n# : 2 = b ^ d\n# : ln 2 = d * ln b\n#. : ln 2 / ln b = d\ndoubling_rates = (np.log(2) / log_params[0]).round()\nprint(f'England doubles every {doubling_rates[0]} days. York doubles every {doubling_rates[1]} days.')\n# Using CDC estimate of serial interval\nr_numbers = np.exp(log_params[0] * 3.96)\nprint(f'R Numbers: {r_numbers}')",
"execution_count": 8,
"outputs": [
{
"output_type": "stream",
"text": "England doubles every 11.0 days. York doubles every 26.0 days.\nR Numbers: [1.29044784 1.10968025]\n",
"name": "stdout"
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "def doubling_time(values):\n days = np.arange(stop=len(values))\n log_params = np.polyfit(days, np.log(values), 1)\n return (np.log(2) / log_params[0])\n\ndef r_value(values):\n days = np.arange(stop=len(values))\n log_params = np.polyfit(days, np.log(values), 1)\n return np.exp(3.96 * log_params[0])\n \ndf[['England', 'York']].rolling(7).apply(r_value).dropna().plot.line()",
"execution_count": 9,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 9,
"data": {
"text/plain": "<AxesSubplot:xlabel='date'>"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 432x288 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEpCAYAAACDc9l6AAAASXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMisyNTQ1LmdlZjY3M2M0NmUsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvfRvX6wAAAAlwSFlzAAALEwAACxMBAJqcGAAAV9JJREFUeJztnXd4VFXawH9n0hsJpNACJLSE0HuVpiIiInZx7V3X1XV1v7Wsdd1Vd9e1K/beBVSqKKJIk15SIAkpJBBSIRBC+vn+ODMQQspMcqfm/J5nnpu599xz3rl38s6573mLkFKi0Wg0GvfH5GwBNBqNRmMMWqFrNBqNh6AVukaj0XgIWqFrNBqNh6AVukaj0XgIWqFrNBqNh+DtrIEjIiJkTEyMs4bXaDQat2Tr1q1FUsrIxo45TaHHxMSwZcsWZw2v0Wg0bokQIrupY9rkotFoNB6CVugajUbjIWiFrtFoNB6C02zojVFdXU1ubi4VFRXOFsUt8Pf3Jzo6Gh8fH2eLotFoXACXUui5ubmEhIQQExODEMLZ4rg0UkqKi4vJzc0lNjbW2eJoNBoXwKVMLhUVFYSHh2tlbgVCCMLDw/XTjEajOYlLKXRAK3Mb0NeqlRxKhLpaZ0uh0RiOyyl0Z+Pl5cWwYcNOvp599tlW9xUcHGyITFlZWQwaNMiQvto9h7Nh/iTY9aWzJdFoDMelbOiuQEBAADt27HC2GBp7cWgXIOHgDhh2tbOl0WgMRc/QrSQmJobHH3+cESNGMHjwYPbs2QNAYWEh5557LiNGjOD222+nV69eFBUVnXZuWVkZZ5999slzv/vuO0DNvAcMGMCtt97KwIEDmTFjBidOnABg69atDB06lPHjx/Paa6859sN6MgV7zNtk58qh0dgBl52hP7k4ieSDRw3tM6FbBx6/cGCzbU6cOMGwYcNOvn/ooYe48sorAYiIiGDbtm28/vrr/Pe//+Wdd97hySefZPr06Tz00EOsWLGCt95664w+/f39WbRoER06dKCoqIhx48YxZ84cANLS0vj88895++23ueKKK1iwYAHXXHMNN954I6+88gpTpkzhr3/9q3EXob1jUeQFKc6VQ6OxAy6r0J1FcyaXSy65BICRI0eycOFCANauXcuiRYsAmDlzJh07djzjPCklDz/8MGvWrMFkMnHgwAHy8/MBiI2NPfkDMnLkSLKysigtLeXIkSNMmTIFgGuvvZbly5cb+THbLxZFXl4EZQUQHOVceTQaA3FZhd7STNoZ+Pn5AWrhtKamBlDKuiU+/fRTCgsL2bp1Kz4+PsTExJx0N7T0aen3xIkTSCm1B4s9qKmC4jToPhIObFWzda3QNR6EtqG3kUmTJvHVV18BsHLlSg4fPnxGm9LSUqKiovDx8WH16tVkZzeZLA2AsLAwQkNDWbt2LaB+EDQGUJwOdTUw6DL1XptdNB6GVugNsNjQLa8HH3yw2faPP/44K1euZMSIESxfvpyuXbsSEhJyWps//OEPbNmyhVGjRvHpp58SHx/fohzvv/8+f/zjHxk/fjwBAQFt+kwaMxb7eexkCIzQC6Maj0NYYzKwB6NGjZIN86GnpKQwYMAAp8jTWiorK/Hy8sLb25sNGzZw5513OtTt0R2vmdNY9Q9Y+wI8kgefXAo1FXDLT86WSqOxCSHEVinlqMaOuawN3V3Yv38/V1xxBXV1dfj6+vL22287WyRNUxSkQHhf8PaDqATY8SnU1YFJP6hqPAOt0NtIv3792L59u7PF0FhDQTJ0Har+jhoAVWVQmgMdezlXLo3GIPTURNM+qDoOh7PUzBxObfXCqMaD0Apd0z4o3AtINTOHU1u9MKrxILRC17QPCs0h/5aZuX8HCO2hZ+gaj0IrdE37oCAZvPygU71iIFED9Axd41FohV4PKSWTJk06Lcz+q6++YubMmS2e+8EHH3D33XfbUzxNWyhIgcj+YPI6tS9qABSlQm218+TSaAxEK/R6CCGYP38+f/nLX6ioqOD48eM88sgjLWY7tKQB0LgwBSmnzC0WohKgtgpKMpwjk0ZjMNptsQGDBg3iwgsv5LnnnuP48eNcc8013H///WRkZBAYGMhbb73FkCFDeOKJJzh48CBZWVlEREQwY8aMk30sXbqUp59+msWLFxMREeHET6MB4MQROHrg1EKohfoLo5FxDhdLozEa11Xoyx+EQ7uN7bPLYDi/5QpElrznvr6+TJo0ieHDh/Ptt9/y888/c911152MBN26dStr164lICCADz74AIBFixbxv//9j2XLljWaeVHjBBouiFqIiANhgvxkGHix4+XSaAzGdRW6EwkKCuLKK68kODiYzz//nAULFgAwffp0iouLKS0tBWDOnDmn5VlZvXo1W7ZsYeXKlXTo0MEpsmsawbLw2XCG7uMPnfrohVGNx+C6Ct2KmbQ9MZlMmEymRtPjWlLbBgUFnba/d+/eZGRkkJqayqhRjaZa0DiDghTwDVZuig2JGgD5SY6XSaOxA3pRtAUmT558Mn3tL7/8QkRERJOz7169erFw4UKuu+46kpK0knAZClKU4m4sx3xUgloUrT7heLk0GoPRCr0FnnjiCbZs2cKQIUN48MEH+fDDD5ttHxcXx6effsrll1/Ovn37HCSlpkmkVDPwhuYWC1EDAGmOJNVo3BudPtfN0desBY4XwX/6wHnPwPi7zjxemAqvjYaLXoPh1zhePo3GRppLn9viDF0I8Z4QokAIkdhCu9FCiFohxGWtFVSjMZyjB9Q2rBH7OUB4H/AJgrxdjpNJo7ET1phcPgCaDZUUQngBzwE/GCCTRmMcx1QxboI7N37c5KXcWfN2OEwkjQdTVgC5W6Cu1inDt+jlIqVcI4SIaaHZn4AFwGgjhNJoDKOsBYUO0G0YbPtI/RPWTw2g0djKwtsgYzUEdIS+50C/86D/DPAPdcjwbV4UFUJ0By4G5lvR9jYhxBYhxJbCwsJG2zjLpu+O6GtlBWWH1LY5hd51KFSXQ1GaY2TSeCbH8iHzVxgwB/qfD/tWw8JbVLlDB2GEl8uLwN+klC0+Y0gp35JSjpJSjoqMjDzjuL+/P8XFxVpRWYGUkuLiYvz9/Z0timtzLF/NjnyauU5dh6lt3k6HiKTxUJK/BVkH0x6Bi9+AB9Jg9C3qe1XrmHxPRgQWjQK+MAfbRACzhBA1Uspvbe0oOjqa3Nxcmpq9a07H39+f6OhoZ4vh2pTlQ3CX5ttE9AfvAGVHH3qlQ8TSeCCJCyBqIETFq/cmE0SPgc3vQMk+h+QLarNCl1KeTDAthPgAWNIaZQ7g4+NDbGxsyw01Gmspy4fgqObbeHlDl0F6hq5pPUf2Q87vMP3R0/c7OAGcNW6LnwMbgDghRK4Q4mYhxB1CiDvsLp1G01bK8iGkhRk6KLNL3i6oq7O7SBoPJHGh2g5qYC+P6A/CSyWAcwDWeLnMs7YzKeUNbZJGozESKZUNvbkFUQtdh8Lmt9WjcUQ/+8um8SwSF0D3kadXxAK1dhPuuARwOvRf47lUHoWaE9Yp9G7D1FabXTS2UpQGh3adOTu34MBSh1qhazyXsgK1tcbkEhmvao4e3G5fmTSeR+ICQDSdUz9qIJRkQtVxu4uiFbrGczlm8UFvYVEUwMsHOg/UM3SNbUipFHqvidChW+NtHJgATit0jedyMkrUihk6KLNL3i71T6rRWMOh3arQ+KBLmm7TeaDaOsDsohW6xnOxKPQQK2zooBZGK0vhcKb9ZNJ4FokLlBdLwtym23SMUXEOBSl2F0crdI3ncuyQsov7h1nX3hIxenCHnQTSeBRSQsr3EDsZgsKbbmfyUj7oDqiMpRW6xnMpK1AeLo1VKmqMqAFg8tGZFzXWUZCsql0lzGm5bVSCnqFrNG2i7JB1C6IWvP2gc4JeGNVYR8oSQEDcBS237Zygvo/lJXYVSSt0jedSVmCdy2J9ug5TJhe9MKppiZTF0GOsdWs09VMA2BGt0DWey7FD1gUV1afrUKg4onJzaDRNUZIJ+bthwIXWtY9KUFs7pwDQCl3jmdRUwYkS2xW6JWJUBxhpmmPPErUdMNu69iFd1eK8nqFrNK3guCVK1EaF3nkQePnCgS0tt9V4Hutego8uatnklrJYlS7sGGNdv0KYF0a1QtdobKelWqJN4e2n7Oi5WqG3O6SELe9Dxi+wf2PT7Y4dgpxNqjKRLXQ2e7rYcX1GK3SNZ2JNLdGmiB6tTC611cbKpHFtitNPBZVt/aDpdnuWAhLirTS3WIgaoBLGlea2VsIW0Qpd45lYaona6uUCED0KaipUWLem/bB3udr2O0+VkztxuPF2e5ZApz6nPFesJcqSAsB+/uhaoWs8k2P5gICgM2vXtkiPMWqrzS7ti9QfoPNgmPaw+kHf9fWZbU4chsw1yrvF2oA1C5bSdAX2ixjVCl3jmZTlQ2C4yqJoKx26K6+E3M3Gy6VxTU4chv0boP95ytOp6zBldmlo7079AepqrHdXrE9AR/Xd0jN0jcZGyqysVNQYQiizS+4mY2XSuC7pq0DWQv+Z6v3I69VM+sDWU22OF8Oa/0KHaOg2onXjRCW0KVfQzpwjzR7XCl3jmZTl2+6yWJ/oMXA4C8oKDRNJ48KkroDACOhuVtSDLgOfoFOLo5XH4NPLoDQHLnkLTK1Unb2nQNFeOJJj86l1dZK/ftN8Wgqt0DWeybF86/OgN0b0aLXVZhfPp7YG0n5U5haTl9rn30HlOE9coGbmX16jcvxc9j7ETGz9WP1mqG36jzafumpPAan5Zc220Qpd43lIaTa52JCYqyHdhoHJWyv09kDuJpXuof95p+8feQNUl8PbU5Vv+pyXIX5W28aK6A9hPdUPiA1IKXn9l3SiOwY0204rdI3nceIw1FW3zmXRgk+AigTUCt3z2btcpU3uPe30/d1HqsjhI/vh3Kdg+DVtH0sI6HsuZPwKNZVWn7Yxo4Tt+49w+5Q+zbbTCl3jedhSS7Q5okfDgW3qkVzjuaT+oMwo/h1O3y8EzHkF5r4BE+81brx+M6D6OGSvt/qU139JJyLYj8tHRjfbTit0jedhay3Rpogerf7xCu1fmEDjJEoy1CKlxbulId1HwLCrjR0z9ixVSctKs8vu3FJ+Syvi5kmx+Pt4NdtWK3SN53GylqgBCh202cWTSV2ptg3t5/bENwhiJlm9MPr6L+mE+HtzzbieLbbVCl3jeZycobfR5NIxRrmy5WiF7rHsX6/uc6fejh2337lQlKryqjdDekEZK5IOcd34XoT4txwkpxW6xvM4lq98iP1C2taPEGqWrmfonkt+klr4dDQn3Rd/arbZ67+k4+dt4saJsVZ1qxW6xvMoO9S2oKL69BgNxWl2rwWpcQJV5cqG7gyFHt5HPRU0Y0ffmXOERdsPcN34GCKC/azqVit0jedRVtD6sP+GWOzo9UPANZ5B4R6QdSpPuTPoN0Ml+qo+ccahujrJY98nERHsx5+m97W6S63QNZ5Ha2qJNkWXIWqbn2hMfxrXwVI9yE4z9MJjlTy3Yg//+zEV2VhRi77nQs0JyFp3xqFvtuayM+cID50fb5Xt3IJ3WwTWaFySsgLoM92YvgLC7J4hT+Mk8pPBO8D6MnLWdnu0gjd/zeCzTdlUVNcBUF5ZwyMXDEDUT7kbM1GNn7YS+p1zcnfpiWqeW7GHkb06cvHw7jaNrRW6xrOoPgGVpcbZ0EEVMrBzLUiNE8hPVDnKTc37dtvCxxuy+MfSFGrrJBcN68Yfp/Xl4w3ZvLM2kxB/H+49p9+pxj4BKlnX3uVw/nMn86u/+FMqJeVVfDhnzOk/AFagFbrGs7Bksetg28ymWaISIPM3FTHqpf9lPIaCZEP9z9MLjvHUkmTG9Q7n6bmD6BUeBMBjsxM4VlHDCz+lEuLvzU2T6nmsJMxVmR5zt0CP0ew9dIyPNmRz9ZieDOoearMM+tup8SyK09Q2vF/z7WwhKgFqK5VHRGR/4/rVOI+yAjheeKosXBupq5M8vDCRQF9vXrhy2GleKSaT4LlLB3O8soanliRTJyXXjY/B19ukkn15+UHSQlYf78Wj3yUS4u/NAzPiWiWHXhTVeBZFqWobYb1nQItYakdqs4vnkG8uA9fZGIX+9dYcNmWV8PCs+EZdDL29TLw0bxhT4yJ5emkKk/+9mrfW7OMYgZyImcaRzV9y8we/4+dt4p3rRtExyLdVcrSo0IUQ7wkhCoQQjS7zCyH+IITYZX6tF0IMbZUkGo0RFKUrDxd/2x9XmyQyDoRJL4x6EgYq9KKySv61bA9jYjpx+cgeTbbz8/bi/RtG8/6No4mNCOJfy/Yw4ZmfeTitP2G1xfxv3AmW3zuZUTGdWi2LNTP0D4AmMtcAkAlMkVIOAf4BvNVqaTSatlKcZqy5BdTiVafedi3uq3EwBcnqhz8oos1dPb0kmfKqGv51ySBMpuYXMYUQTIuL4vPbxrH47kmck9AZ/4QLqPMOYK73RmWGaQMtni2lXAM0GSYnpVwvpTxsfrsRaD6/o0ZjT4pSIcJghQ5mTxc9Q/cY8hPV2kgb+S2tkG93HOTOKX3oG2VbqonB0aG8cOUwnrlqHKa48yH52zanajbahn4zsLypg0KI24QQW4QQWwoLda1GjcEcL1bFLeyi0BPUomgjUX0aN6O2Bgr3ttncIqXk2eV76BUeyF3T2rhmM+hSKC+GzF/b1I1hCl0IMQ2l0P/WVBsp5VtSylFSylGRkZFGDa3RKOzh4WIhaoAKEy/ca3zfGsdSkgE1FW1W6OvSi0k6eJS7pvZpMU95i/Q9B/w6QNLCNnVjiEIXQgwB3gEuklIWG9GnRmMzRWaFbqSHiwWLe5s2u7g/BcYsiL65Zh+RIX7MtTGas1F8/CH+AkhZbFNpuoa0WaELIXoCC4FrpZSpbe1Po2k1Rang5QthvYzvu1Nv1bd2XYS6OmdL0Dbyk0B4QUTrfL0BEg+oKkI3TYzFz9ugSNNBl0JFKez7udVdWOO2+DmwAYgTQuQKIW4WQtwhhLjD3OQxIBx4XQixQwixpdXSaDRtoTgdOvUxNJT7JF7eSgG09xl68vfwr66w9kWoq3W2NK0jPxnC+6pZcSt5a00GwX7eXD225SpCVtN7KgR0hM3vtnpx1Bovl3lSyq5SSh8pZbSU8l0p5Xwp5Xzz8VuklB2llMPMr1GtkkSjaStFafYxt1jonKBn6HuWKJPAT4/DezOV37+7kZ/YppS5OSXlLN2dx9VjexIaYH0mxBbx8lHFqNN/hE8ublUOfh0pqvEMaqvhcKZ9FkQtRA2AowfgxBH7jeHKSKlSvSbMgUveUSau+ZNg8zvOlsx6Ko/Bkew22c/fXZuJScCNE2OMk8vCpPtg7huw/3d4ayocsi1ts1boGs/gcBbU1UCEHXOtWPyWC/fYbwxX5kg2HM2FXpNgyOVw10boORaW3g9H85wtnXVYTGatzOFScryKLzbvZ87Q7nQNDTBQsHoMuxpuXA61VfDuubB3hdWnaoWu8QxOerjYeYYOp8LG2xvZ69U2ZqLaduiqZpRwKoeOq3My5L91JpePNmRRUV3HbZPtXFQ6eiTc9otajF96v3o6sgKt0DWewUkfdDva0EN7gG9I+10YzVqnFu0iB5zaZzFxFbuJLf3gdnUPQ21fzCw5XsW7v2VybkJn4rq0sQC5NYR0gfF3q6eiA9usOkUrdI1nUJQKQZGqwpC9EKJ9pwDIXgs9J4CpntoI6Qo+gVC8z3lyWUtdrSom0Xf66Z/BSl5bnc7xqhr+el7r3R1tJm4mmLwh5TurmmuFrvEMitLtaz+3YKleZOUjsMdQekCtU1jMLRZMJuUqanlCcmVyN8PxAhgwx/ZTD5fz8YZsLhsZTf/ODpidWwjoqNwZk7+z6junFbrGsVQdhx8fV77ARlKcZl9zi4WoBDhRAmX59h/LlbDYz3tNPPNYeB/3MLmkLFbBYf1m2HzqCz+mgYA/n+OEAicD5qgf00O7WmyqFbrGcZSXwEdzYd2Lxrq6lZeoxEb2XBC1YFlMe2MiPD8A/pcAL4/wnBwvJZmw5C9QVX76/uy1KtdIl8FnnhPRDw5nQ02VY2RsDVKqoKje08C/g02n7jl0lIXbc7lxQgzdwuzk2dIc8bNVZGvy9y021Qpd4xiOHYIPLoC8HRDSDQ7tNq7vIjsm5WpIj3FqoSrufOh7NsScBSX7IHON/cd2BHuWwpZ3Yd1Lp+/PXg89xzUehRveF2Stcmt0VfJ2Qul+GHChzaf+Z8VeQvy8uXNqHzsIZgVB4RAzySqzi1boGvtTkgHvzoAj++EPX6vAlPwk40LHix3gsmjB2xfO+ydc9Kp6XTxfzVw9xTf9cJbarnvxVMHtsgK16NyYuQVOmbqKXNiOnrJYzXLjZtl02qbMElbtKeDOqX0JC2xdWThDSJijvuctfM+0QtfYl9oa+HCOitC7/nu1wNNlMFQfV4/3RlCUBiYf+yTlagkhVIk6TzG5HM6CDtGAgB8fU/uy16ltzKTGzwk3z1xd2Y6eslgt6AaFW33KsYpqHlywiy4d/LlhQoz9ZLOG+AsBoWbpzaAVusa+5O2A0hyY9R/oPlLts9hhrVjksYqiNHM2RG9j+rOVCA9S6EeyofsImPRnlZs7e716+QRB1ybKBQd0hMAI11XohXuhaK9N3i1SSh74eifZJeW8eNUwAnztkPDNFkI6Q8/xWqFrnEzGL2rbe+qpfZHxyrfWKDt6cZpjzC1NERmn3OFakUzJpairU4ubHXvBhHvUTH35/0Hmb9BjjEoe1RThfV1XoaeYFxPjL7D6lDd+3ccPSfk8dH4843pbP6u3KwkXtZgcTit0jX3J/BU6Dz69GK+3n5rV5tuWeKhRqiuU6cbZCh3cJ/y9KcoOQW0ldIwB30A490n1o1uYcqb/eUMiXFmhL4boMdChm1XNf0sr5L8/7GX2kK7cPCnWzsLZgBULulqha+xHdYXKGhc7+cxjXQa3fYZeWwMLboa66tOfAByNRaG7u9nlsNlLpWOM2g66VHn1gErI1RzhfZVvfsVRu4nXGmRJJuTtJKXjFF74MZU7Pt7KrR9t4estOZSWV5/Wtq5OkniglHs+307fqGCeu3QIQggnSd4Iod0henSzTZxkdNS0C3J+VzO+3lPOPNZlMOz6AsoKIbgV9WXr6mDxPSo/98znnKvQQ3uCd4AHKPQstQ2LUVshYM4ryo0xuoUyBxZPl5J90G24vSS0meXfvMMs4PYt3cghjZjwIKpq6vgxOZ+HvXZzVr9I+kUFs/tAKbsPlHKsooYQP2/evHYUQX4uqB5nvwAMafKwC0qs8Rgyf1WuYr0mnHnMsjCavxuCp9vWr5Sw8hHY8SlMfQjG3dHyOfbEZFImnyJPUOgCwnqc2hfZH85/ruVzLTEARekuo9BX7y3AO2cDhQHRvH7rpfSNCsbfxwspJTtzS1m2O4+lu/L4La2Q+C4dmDO0G0Ojw5jUL8I5AUTW0FhgVz20QtfYj8w1yrPFr5HcFyc9XXZDHxsV+pr/wsbXYeydMOVvbZfTCCLjYP9GZ0vRNo5kQ4fuao3DVjrFAsJl7Oil5crlcIl3Np36TSGye+jJY0IIhvUIY1iPMB46P57aOom3l2dYnz3jU2hcj4qjKuVnY+YWgMBOyovCVjt68T5Y/TQMuRLO+5cyC7gCkXHKPbOyzNmStJ7DWcrDpTV4+0FYT5dR6E8uSUKWFRApi/DqPqLJdkIIj1HmoBW6xl5kr1Ph4I0tiFroMsjmElskf6u20x9tVQpUuxEZr7bu7OlyOOvUgmhrCO/rElkXf0zOZ+G2Azw0tELtcBETkCNwof8IjUeRuQa8/ZW7WFN0GawUYPUJ6/tN+lat9Ne387oCEW7u6VJdAcfymlTohccq+WpLDrd/vIXL568nvaCRJ5GIfuoJyomphQ8fr+LhRbuJ7xLChZEFgICuTS8iehpaoWvsQ8avKpmTj3/TbboMVrN4awtGFO9T0aUDLzZGRiPpFKvSD7jrwuiR/WrbIH3CzpwjzH1tHaP/+RP/980uduaUsq/wOFe9tYG9h46d3kd4X6gqc2pq4X8sSebw8Sqev2Io3od2qB+ZxtZwPBSt0DXGU1YIBUnNm1vg9IVRa7CYWxIuarVodsPLR+U0cdcZ+pEGPujA3kPHuO69TRQcreD+c/uz9J5JbHhoOl/dPh6TEMx7eyPJB+v5nVtyujgpSddvaYUs3H6AO6b0YWC3UJV2oh2ZW0ArdI09yPxVbWOnNt8uLEbVd7RWoSctUiac0Og2CGdH3DlJl8UH3azQ9xeXc+27v+PvY+LL28fzp7P7MbBbKEII+kYF8+Xt4/HzNjHv7Y3szi1V5zqxvuiJqloeWZRI74gg7p7eV6VrPpYHXYc5XBZnohW6xngy14BfaNPJnCyYTOaFUSsUevE+1W7gXENEtAsRcXA4U9mj3Y3DWSo4KjiKgqMVXPPu71TV1vHxzWPp0SnwjOaxEUF8dft4Qvy9+cM7GzlUWmF2efR3ikJ/cVUq+0vK+dclg/H38YKDO9QBPUPXaNpI5q8q94c12Q+7DFY5Xerqmm+XtEhtXdHcYiEyDmSdipZ0N8wui6Unarj23U0Ul1Xy4Y1jmq2f2aNTIJ/cPJaKmjqeWZ5Sr76oYz9/0sFS3vktkytH9TiVSOvgdkC0GIjjaWiFrjGWkgylHKwNxe8yWC2kHclqvl3St9BjrOuaW6BeThc3LHZxOBs6xvDflXvZV1jG29ePYmiPsBZPi4kI4vbJvflux0E2ZZaYk3Q5zoZeWyd5aOFuOgb68NCs+FMHDm5X98Mv2GGyuAJaoWuMJX2V2vY527r2nQepbeJC9WNQW31mm6J0lSIgYa4hItqN8L4gTO5nR5cSDmdR6t+NzzftZ96YnkzoE9HyeWbumtqXbqH+PP59EnWd+qof9Mbuox34cH0Wu3JLefzCgacqCknZLhdEQSt0jdHs+1lFDFo8HloiKkGVcPv5H/DycHi6M7w0FBbfC1nrlCkm2Q3MLQA+Acrtz90U+onDUHWMH/MC8PU28aez+9p0eoCvF49ckEBK3lE2HOkIdTWnFlntyKHSCp5fuZcp/SOZPaTrqQPH8pTrZDtbEAWdy0VjJDVVakF08OXWh+T7+MO9O1Xi/sNZ6lW4F3Z9DVs/gNAeUFtlNrd0t6PwBhEZ734K/bAqBbjigD+3TI0lKqSZ2IEmmDW4C+N7h/NaYjYTQZmd7Jyj/h9Lk6mpkzx10cDT09we3K627XCGrhW6xjhyNyl7eF8rzS0WAjupepX1a1ZWHYc9y2D3V8qMM/3vxspqLyLjIP0nlavdWSXxbESWZCGAo/7duHVy71b1IYTgiTkDufzlA+CL+lGzoiBDa1mTWsjSXXn85dz+9AoPOv3gwR3K9NXOFkRBK3SNkaSvUulyWwoosgbfIBhyuXrVVquSde5AZJwquHHYyVWUbCAzPZnewIVTxhPi30yZuRaI6xLCJePiObA1nOCcJEJbPqVVVFTX8th3icRGBHH7lEZ+gA5uV09Kvme6W3o62obeGFI6NR+F27Jvlao96W/wv7KXj+tkVWwJi+/9b8+7xXeork6yN2U3JYRyxcT4lk9ogXvP7keWiOZwtkH1YhvhzV8zyCou5x8XDcLPu0Hx5na8IApaoZ9JXR18eQ28eRaUHnC2NO5DWSHk7bTeu8VT6TwQpj0COz+Hnx53tjQt8v3Og4ScyIWOvc5Ujq2gY5Avwd0TiKrcz6aMIgMkPJ2souO89ks6Fw7txqR+jXjiHD0Axwu1QteY2fi6KmtWsAfenWHdApee0UPGarXta2OxCk9k8l9h9C2w7iVY/6qzpWmSyppa/vPDXvr4FNOxe3/D+h0wZDSBopJ3l65BGvx/8dSSZHy9TPz9ggGNN7AsiLZDDxewQqELId4TQhQIIRpNXC0ULwsh0oUQu4QQTWeTd3UO7oCfnoD42XDrKuVd8d55kLul6XMqy+DDC+Gdc1T+iPZK+ioIDIeu7jszklKybHceF76ylk82ZrdeGQkB5/9buVmufAR2fQV1tVBeoqIonZS8qiEfb8jm0JEyushCRFvyoDfAt4tSthV5e1iZbFzmxdV7Cvh5TwH3nt2Pzh2a8MTJ26nWcboMMmxcd8KaGfoHwMxmjp8P9DO/bgPeaLtYTqCyTFWQD4pUhXG7DoWbfwD/MKWw9y4/45TqijIqProMmb2euvwkeOdcNbNvb9TVKf/z3tNcq+iEDeSUlHPzh1u469Nt5Bwu5+/fJnLnJ9s4Ul7Vug5NXnDJ2xBzFiy8FZ7qBP+OhVdGwKujIOMXQ+W3ldLyal75OZ2LekuErG19paLGMOeGHxtSyH9+2EtNbQtpHaygqqaOp5Yk0zsyiOsnxDTdsCBFxUD4uGhNUDvT4n+flHINUNJMk4uAj6RiIxAmhOjaTHvXZMXf1OzpkjeVGx1Ap95w80oVAfj5VVR8eTPfrd3B7R9vYfqzK/j9X+fhk7uReyrv5NITj1J+ohz53gzIWuvcz+Jo8hPheIHt7oougJSSd37LYMYLa9iYUcyjsxPY/Mg5PDwrnp9S8pn10m9szmru698M3n5w1Wcw7e+qmPXMZ2HuG+AbrFIZOJHXf03naEU19ww3R1caOEMnKBwCI7iw2zHSC8pYsC23zV2+vy6TzKLjPDY7AV/vZtRW4V6IMM585G4YMZ3qDuTUe59r3uc+JC2C7Z/AWX+B2MmUnqhmbVoRC7fl8saWYzzd5SW+DLwKU/Iipvw4i9jsr3nF6wUmmpLYOvxp5t10H1HxYzn36KMcqAlFfjRXBca0F/ZZwv3dz36+PPEQTy9NYWLfcH76yxRunhSLj5eJ2yb3YcGdE/DxNnHlmxt4f11m6wbw7wBT/gpTH4Rxd8Kwq9V12ru85YRkduLAkRO8vy6Li4d3J8ZkNokYqdABIuOIrslheM8wXvgxjRNVta3uquBoBS+vSuOcAVFMjYtqumFNlUofEdl2bx13xQjn3sb8yRo1PgohbkOZZejZs6fVA9TWSWrq6gxZhW+UDa9B5ACY+hArkw7x8KLdFJWdetQO9vOmT9R1VA24lEsO/pcH8+ZDDXDhS4wZeQMA4/uE89WWKC5bHMzLpv8xZuEtKupx8gPu43LXWtJXqZwsIV2cLYlNVNbU8uzyPcR3CeHNa0fhZTr9Pg3tEcaSP03i/q928uTiZLKLy3l0dsIZ7Wwm/gJI+V4t4EWPbFtfreD5lWqh//4ZcbB5IXj5qtS3RhLRH5G0iIevjOfyNzfy5pp9/Pmc1s2cn1uxl+payd8vSGi+Yck+VQHLkiStHWKEQs8F6hd4jAYONtZQSvkW8BbAqFGjml1xyis9wZrUQtakFbEuvYjyylrG9u7E9PgopsdHnRkd1loqSuHAVirH38cjC5P5ZmsuCV078PwVw+jZKZCoED+C/OpdJnku7P5GRQHWK4UmhODK0T0ZGxvO/Z+Hc03B81y8+mkVYDL7RfD2NUZeV6OmCnJ+hzG3OVsSm/lofTb7S8r56KYxTSrpEH8f3rhmJP9alsK7azPJPXyCl+cNI9C3Df86/Waohbu9Sx2u0BMPlLJo+wFum9yb7mEBakbbMVbZ/I0kMg4qjjA6spYLhnRl/q/7uGJUD7qF2Wbb3rb/MAu25XLn1D7ERLTwP2/xSGvHJhcjFPr3wN1CiC+AsUCplDKvNR2dqKplya6DfLZpP9v3HwGgcwc/zhnQmdAAH37ZW8CTi5N5cnEysRFBTOobwaR+EYzvE06H1ka4Za0DWcf9m8NYVpbL3dP6cs/Z/Zq20wmhohebICYiiA9vO4tr3vYhJ78z9+z4VNVrvPJjCOjYOhldmcIU5Q3U3fEzzbZQcryKl39OY0r/SCb3j2y2rZdJ8OjsBHqFB/LE90lc+eZG3r1+FFFNeVq0RGAn6DVBpTY4+7HW9dEKqmvr+NuCXYQH+XHXFHMCruIMtVZkNJGnimY/dP5IfkrO59nle3h5nvVeUOVVNTzw1U66dPDnj9OsSBhWlAoIrdCbQwjxOTAViBBC5AKPAz4AUsr5wDJgFpAOlAM32irEnkNH+WJTDgu35XK0ooY+kUE8eH480+Ki6N85+GTinUdnJ5BdfJyf9xSwJrWQBdty+XhjNiYBQ6LDGNu7E+NiwxkV09HqEOaMzcvoKn3Z6xPPN3eOZkTPtivdYD9vPrxpLFe+VUduSRee3f8mpi+vhesXe575JW+n2rZUncjFeHlVGscra3ikKX/mRrhufAzRHQO4+7PtXPz6et6/cXSzBSCaJf4CWPGgmiHbQ6E2wuur95F08ChvXjuS0EAfFTtRkmF97npbiDiVGz469ixum9ybV35O5/oJvRjZq5NVXTy1OJnM4uN8estYgv2smHsW7oGwHu0y5N9Ci1dJSjmvheMS+KOtA9fWST5cn8XXW3NIPHAUXy8TMwd14Q9jezImttPp2dPq0Ss8iBsnxnLjxFiqaurYvv8wa9OL2LCvmPfWZvLmrxmYBAyODmNK/0imxkUyNDqs0Ufqd37L4Ky01ez1H8Q3d09TX3KDCA304eObx3Llm3U8e6ych7PeQ6avQvQ7x7AxXIK8nSr9bcfYFpseq6gmp+QEhWWVjOrV8XRTlgPZV1jGJxuzuWpMT5sV8vT4znx1+3hu/GAzl76xnjevGcmEvtbnDj9J3Cyl0Pcsgwl3236+jSQfPMorP6cxZ2g3zhtoXus4dghqTkC4HX5QOnRT9WKLUgG4Y0ofvtqSw5OLk/n2romYWliHWL47jy8253Dn1D7W52YvTD31Q9JOcVrGo5S8ozz+fRIDu3Xg8QsTuGhYdzoF2WZn9vU2MbZ3OGPNZadOVNWyff9hNmaWsDatkFd/TuPlVWl0DPRhfJ9wBnUPZXD3UAZ2C+XVn9NZvG4bt/jnUj3pFnwMVOYWIkP8+OSWsVw9v5Jryxdz/LO/8sWwj5kxqAtjYjrh7eWePtunkbcTugxp0v/8WEU1f/5iB1uyD1N64lTRgxE9w/js1nGq/qODeWbZHvx9vLivlYt0g7qH8u0fJ3Lj+5u47r1NPHvpEC4baWMlpY691ELyXvsr9OraOh74eidhgb48OWfgqQOWUnn2eEIQAiL7n7RrB/l587eZ8fzlq50s3H6g2et18MgJHly4m6HRofzlXCvvUV2tqpTUe4oR0rstTlPo4cF+LLvnLBK6dTCszwBfLyb0jWBC3wj+cm5/Dh+v4rf0In7ZW8DmrBKW7T49kvN/8QWQBT59pxomQ0O6hQXw3b3T2bPiz4zd+QjFW77h6g2jmRoXyfs3jG7yScQtqK2BQ4kw6qbGD9dJ7vl8O2vSirhydA96dQokumMgJeVVPPZdIvd9uYPXrh7R4mzNSHbmHOGnlHz+el4ckSF+re6ne1gAX98xgTs/2coDX++k5Hglt022sqiHhbhZ8Nt/4Xix8t22E6+tTic5T5laOtafNJVkqK29TD4RcadSQgBzh3Xnww3ZPLt8D2NiOtEz/EzTSG2d5L4vd1BdW8dLVw3Hx9pJz5FsqKlo1x4u4ESF3jXU31Bl3hgdg3yZM7Qbc4Z2A+BIeRWJB46y+0Ap3TsGMCdzmYoE7TLErnKEBvgw9qI74cCHvCSX0ifuSl78OYPvdx7komHu5bJ/GsVp6pG9Cfv5P5emsHpvIU/PHcQ1406PRKysruXppSk8u2IPD8+y3o7dVt7+LYMQP2+uG9/2yMjQAB8+uHEM9325g38t20N1rbRu8c5C/CxY829I+0H5p9uB7fsP8+rP6Vw0rJ6pxUJJBph8oIOd6rRG9oednylPMv9QTCbBs5cMZt7bG7n8zfV8estY+kadMnkdr6zh0e8S+T2zhP9cNqRlr5b6FCrTTns3uXjAM7/1hAX6MqlfBHdO7cOcIV1V+HXsWca7bDWGyQumPYKpOJV7IrczuHsoTy9N4WiFY2ov2oVmFkQ/2ZjNe+syuWli7BnKHODmSbFcN74Xb63J4OMNWXYWVJF7uJzliYeYN7Znm/J+18fX28RLVw3jomHd+M8Pe3nxp1Trc8B0Hab8v/csBVSe7+93HuQfS5LJKSlvs2z7Csu4+cMtdAn154kLB57ZoHifCiiyVyGOkwujqSd3DejagS9vG09tHVzx5kYSD5QCsDu3lNmvrGXR9gPcM72v7SasIrPLYmT79XCBdqbQT6MkA47m2meFvykGzIGuQzH9+gz/vLA/RWWVvPBjasvnuSp5O8E74IxCDr+lFfL490lMj49q0otECMFjsxM4Oz6Kx79PYm2a8alWG/L+OlWZ54bmcoG0Am8vE/+7YhiXjYzmxZ/S+O/KvdTVWaHUhaC673nUpa/iH4u2MvZfq7jn8+28uzaTi19fx/b9h1stU/7RCq57dxMC+PjmsaebWiyUZNrXw8Zi/ig6PWNpXJcQvr5jPAE+Xsx7eyNPL0nmkjfWUVFdy2e3jOMvM+JsN0UWpkJwZ890DbYB5yn0sgKnDQ2cSo4UO9VxY5pMMP0xOLKfIQXfcfWYnny4Poukg6WOk8FI8naqMl/1nnCKyyq5+7Pt9IsK5uV5w5uNqvT2MvHK1cOJ7hjIv3/YY3iq1focrajmy805XDCkq83BLdbgZRL8+9IhzBvTg9dW72PsM6v4+7e7WZtWRHVtHTkl5Xy/8yBPfJ/E5fPXM/nfqxnw6Apu3hiJqeYEGdt+ZHL/SD65eSw/3jeZAF8vrnprIysSbQ/pKC2v5rp3N3GkvIoPbhxDbGOmC4vLoj0VescY8PJT7oQNiI0I4qs7xhMR7Mc7azM5Z0Bnlt97FuP7tHItoXBPu/Y/t+C8ul5HD6pq4876Rc38VT3uWlud3ij6ng3dRsCOT/m/a25gReIhHv02kW/umODQxcE2U1cHebtg6FWn7X5pVRpllTW8Mm+4Vb7Dgb7e3Dm1Dw8t3M3a9CLO6td8kE9r+XJTDmWVNdwyyX4KzGQS/HPuYCb2jWDZ7jwWbD3AJxv34+MlqK5VP1YBPl4M6t6B4T3DiAj2o5t/F1j7HK9PhYCzTwXdLLprIrd+tIU7P93GQ+fHc+tZva2atR6vrOGWjzaTWXSc928czeDoJqpHleVD9XH7fv9NXiqxXWHjT6HdwwJYcOcEkg8eZWLf8NY7CEip3COHXNEGYT0DJxZqlCoHyODLHD90XZ2qTt//fMcH+ggBfabB2hcJ9a7ioVkDeODrnXy9NYcrR1uf38bpHM6EqmOn2c/T8o/x6e/7uXpMT/rZ4N99yYjuvPhTKq+tTreLQq+ureP9dZmMje3UtIIzCJNJMHtIN2YP6caJqlp+TS1kU2YJsZFBDO8RRnyXkDPdVVP6ElCw67RdEcF+fH7rOP7ylVpwPXikgsdmJzT7o799/2Hu+3IH2SXlvDJvOBOb848/6eHScvxAm4jsf6roRCN0CvJtvPKQLRw7BJVH23VSLgvOM7mYvBvNMe4QDu1STwfO8lntMVYlETq4nUtHdGdEzzBe+imNagPyRjuMvB1qW0+h/2tZCoG+Xvz5HNuKI/t5e3HrWb3ZmFHC1uzW242bYtnuPA6WVnDrWY6JyLQQ4OvFzEFdeOzCBK4d14tB3UMbjz3oNgIObjtjt7+PF6/OG8FNE2P5YH0W9365g6qaM78jNbV1vPhTKpfN30B1reTzW8cxe0i35oWzt8uihch4OJwNFUftN4bFpKNNLk5U6P6hkP6jqujuaDJ/VdtYJyn06NFqm/M7Qgj+OK0vB0srWLa7VSlwnEPeTpWlzzwrWpNayOq9hfxpel/Cg2337543pidhgT688Uu6oWJKKXl3bSa9I4OYHt9M6lVn0n0EHMuDo2fef5NJ8OjsAfzfzDgW7zzIzR9upqyyBikl2cXH+WpLDpfN38CLP6ko0OV/Potxva2wQxfvU5OqUDs/FcZORj2N/2S/MczRqO3dBx2caXLx7wAVhSpTX8wkx46dulIpog5OqsMR2EnNJnI2ATAtLoo+kUG8tSaDOUO7uUewUd5OiEoAb19qauv459IUenYKbL6aTDME+Xlz44RYXvgplZS8owzoakyMwpq0InbllvLPiwe57hqFpaDxwe2NfieFENw1tS8RwX48tHA3F7z8G5XVdRw6WgGoiORX5g3nwqEtzMrrU5IBYb3s57JoocdYCIxQdXoHXWKfMQr3gl+o8nJp5zhvhu7XQc3wHG12KT0A2etOS33rFHqMUT9mUmIyCW6b3Jukg0dZv6/YuXJZg5RKoZvNLV9tyWVv/jEeOj++TTnrr5/QiyBfL974ZZ9BYkqeX7mX7mEBXD6yR8snOIsuQ1Q63WZszQBXjOrBW9eOpIO/D6NjO/H03EH8eN9kfn/obNuUOTguKZjJSwVQpa6Emkr7jFGUqmbn7jARsjPOU+jCpGbmqSscO27SQkDCICcsxtanxzhlxzcXDL5oWHcigv14a02Gc+WyhtIcJXvXodTU1vHCT6mMienEzEFtK3ARFujLNeN6sWTXQbKLj7dZzJXJ+ezKLeXec5pJh+wK+AZC1IBG7egNOXtAZxb/aRKvzBvONeN60a9ziO1PHhaXRUd5eMXPVgvomb/Zp//CPe0+oMiCc7/l/c+H4nQoMtZu2iy7v1GPuBE2hGjbgx5j1Tbnd0AtgN0woRe/phay99AxJwpmBScjRIexKbOEwmOV3DQpxhBTkaUE3DPL2lZsu7ZO8r+VqfSOCOKS4W6QXqHbMDiwTSlbe3O8EKrKHJa2l9gpqo7qnsXG911eoj5POw/5t+BchR43U20dNUsvSlPeGYObLlDhMML7Kh98s0IH+MPYXgT4ePH2by4+S8/bqUwEnRNYkXQIfx8TU/obs+AY1cGfe87ux4qkQ/yYnN/qfpbsOsje/GP8+dz+7pHVstsIOFGiiqHYG0d5uFjw8Yd+56pUwXWtry3aKCcXRLXLIjhboYf1hKiBjlPou78BBAy00+KMLZhMED3m5MIoqGRiV47uwXc7DpBvXvBySfJ2QmQ8dV7+rEzKZ0r/SAJ8jcuHc9vk3sR1DuGx7xIpq6yx+XzlxpdGfJcQZg920sK3rXQfobZWmF3aTLEd0+Y2RfxsOF4AuZuN7dfisqhNLoCzFTpA//Mge72yydoTKWH31yoZl7O8WxrSY4zKc1FecnLXTRNjqa2TvL8uy3lytUTeTug6hJ25Rzh0tOLMLH5txMfLxDOXDubQ0YqTBY1tYcG2XDKLjnP/jDjX9WxpSNRA5SRwwAEKvSRDPWGFOTCQrd+5KrNjisFml/xk8Am0v/ulm+B8hR53vgqySV9l33EOblcJ/V3B3GKh5zi1rTdr6RkeyHkDu/DF5v1UVBv8eGoEZYUqbLzLYH5IysfbJDg73nh3sRE9O3LN2F58sD6LnTlHrD6vorqWl1elM7RHGOcMcFG/88bw9lUFL1rwdDGEkgylzL2ML+rSJP6hKpBvz1Jj1wn2b4DoUU0WWGlvOP8qdB+p/FTtbXbZ/Y2aAQ240L7j2EK3EWqmVM+ODqp25ZHyahbvPOgkwZqhIAkAGTWQFYl5jO8Tbmjpvvr8dWYcUSF+PLhwt9VRtC+vSuPAkRP8bWYrMvY5m+4j1NNPnZ0jhkv2OT6HESizy+FMKEg2pr+Ko5CfCD3HG9OfB+B8hW7yUmaX1JVQWWafMepqlbtivxmulV7TNxC6DjnNjg4wrncn+ncO5qMN2XbNQNgq8pVC32fqRVZxueHmlvp08PfhyTkDSck7apU7Z/LBo7y5JoPLR0ZbX4fSleg2XOUkKbaj15eU9k+b2xRxswABKUuM6S9nE8g6rdDr4XyFDjDyBqgshR2f2qf/7HUqtNoZicBaosdYOLD1tBQIQgiuHdeL3QdK2WGDucEh5CdDUBTLMmoQAmYk2Dc677yBXZg1uAsvmiNIm6Kmto4HF+6iY6BPkznYXZ5uloVRO5pdyovVj4YzFHpIZ7VuZJT74v4N6gnXkkpD4yIKvccYpdg2vKbqVBrN7m+UH2z/mcb33VZ6jIHqcvXoWI+LR0QT7OfNxxuynSRYE+QnKnfFxEOM6NmRqA7+dh1OCMHTcwcTGuDLfV/uoLKm8XWFD9ZnsSu3lCfmDCQs0LZi4y5DZJxa4LOnp8tJDxcnmFwAEubCod2wf2Pb+9q/QT3h+gW3vS8PwTUUOsCEP6lCr0YHH0gJqT+oVXYf4wsbtJmTAUanm12C/by5dER3luzKo6jMTiHTtlJXC4V7OBYaR3LeUWba0dxSn05Bvjx7yWD2HDrGSz+lnXF8f3E5/125l3MGRHGBu7gpNobJS6VTsKeni8Wc44wZOsDI6yG4C6x8tG2LozWVkLsFek4wTjYPwHUUetws9SVb97Kxq+B5O6HsEPQ7z7g+jSQ0WhXpbWTGcu34XlTV1vHl5hwnCNYIJRlQU8G2CpU3xJ7284ack9CZK0ZFM//XfWzNPuXmmVFYxv8t2Im3ycQ/5g5yv4XQhnQbodI72+NJFdSCpLe//fOgN4VvEEx7CHI3qYRdreXgDqithF7afl4f11HoJi8Y/0f1uLl/g3H9pq0EhJqhuyo9x6nP3OCHrG9UCBP6hPPZ7/upcYVc6eYF0aUFnRjQtQM9wwMdOvyjsxPoGhrA/V/t5N8r9nDO/35l+vO/sjGjhEdnD6BrqAs+gdlKt+FQU6GUuj3IT1R5YxxRGL0phl2jQvV/eqL16bP3r1dbvSB6Gq6j0AGGXg0BnWD9K8b1mbpCuUYGubDXQ6/xatH2cOYZh64b34sDR06wao+Ta7AC5CchhYnvDoQwq42JuFpDiL8P/718KNkl5by5JoPOHfx4cs5A1j843b2qPTVH37PVDHrbR/bpPz8JOg+0T9/W4uUN5z6pzD/bPmxdH9kbILyfa/9fOwHXUui+gTDmVpVSt+hMW6nNlBUqe6QrLobWp9dEtc0+88nknAGd6Rrqz3trz1T2DqcgmWK/HtSa/LhitHPS0Y7vE84Pf57M1r+fw6e3jOP6CTF2KfrsNAI7qeC3nV8YHz1dVqASWXUeZGy/raH/TPW9/+VZqLQxGV1dHeRs1OaWRnAthQ4w+lYVALThtbb3lf4jIKH/jLb3ZU8i4tSTSfb6Mw55e5m45aze/J5ZwuaskkZOdhx1hxLZVtGN8wZ2obOdvVuao3/nEPf1ZLGGsbdDzQnY9rGx/Vo8qZw9QweVu/zcp9QPzPpXbTu3MAUqSvWCaCO4nkIPjoQhl8OuL6Gmqm19pf4AIV1VAQFXxmRStsD9Zyp0gKvH9CQ8yJdXfnZgmuGGVJZhOpLFrupo/jDOQ8wbrkqXwdBrEmx629jshOY1EKJcQKGDCtlPmAu//RcW3XlKvpawTHz0DP0MXE+hg3ocqy5XATetpbYa9v2sFkPdwfOh1wTlRdJIXckAXy9uOas3a1ILbcprYigFKQCUhvRjvDU1KzVtY+ztULrf2Ipe+UlqghPkQvdv9gsw6mZI/hbemAAfXwxZa5s/Z/8G9TnCejlERHfCNRV6r4mAgMw1re9j/wYVEeeq7ooNscw2mpilXzOuJ6EBPry62jmz9Ny9WwAYOmqi+7sGugNxsyC0B2x607g+8xNdw9xSn8BOMOvfcF8SnP2Y+tH5YLZK4tUYUqq1pp7j3WOi5mBcU6EHdlIRYG1R6Kk/KFt876mGiWVXugwFn6BGF0ZBeXjcODGGH5Pzmw2BtxfZyZs4Lv2ZMVGHWTsEL28YfYv6H8g3IJlVbbUqpuxqCt1CYCc46364Z7ty3VxwS+MpEI7sh2MH1ROt5gxcU6EDxE5WwQdV5a07P22lqlnqLmHBXt7Qc2yjC6MWbpgQQ7CfN685eJZeWl6NT3EKRUF96RDg59Cx2zUjrgPvAGNm6cXpUFvlGh4uzeEbBPO+UBlYP7sSjjQIqrOYoLT/eaO4sEKfor6ADVLLWkVJpipN5S7mFgs9J6j0tOWNe7OEBfpy7fheLN2dR3qBnTJTNsKCrTn0I4eQXkMdNqYGNWsdcgXs/LLJ74TVWBYcXXWGXp+QzvCHr6D6hFLqJ46owhjvzoAVf1PlG6PcNAGbnXFdhd5zHJi8W2d2SVuptq7urtgQy2NkMz9iN0+Kxc/bxMurDPDTt4KqmjqWrd9GR1FGp9jhDhlTU49RNykXxraEyYOyn5t8VDCOOxA1AK74SFX0ej4evrwGjh2C8/8Dt69xbqSrC+O6Ct0vREV42qrQpVQuj5HxzktA1Fq6j1R2/+x1TTaJCPbjpomxfL/zINv227lsH/Dq6nSCSs2FeN1hdudpdB2qqgu1NYd4fpLK5ujtRv77fabBRa8rU+TlHyr7+tjblFlG0yhWKXQhxEwhxF4hRLoQ4sFGjocKIRYLIXYKIZKEEDcaIl3sZJXbpaLU+nNyflfujqNvMUQEh+Ljr5R6M3Z0gLum9SUqxI8nFydTV2e/AhjJB4/y+up0Los+onbox1zHIwTEXwgZq1WFntbiCiH/rWHolXDddzBwrp6VW0GLCl0I4QW8BpwPJADzhBAJDZr9EUiWUg4FpgLPCyHaPhWInawqkjTh+dEoG14F/zAYdnWbh3cKvSaoDJHNVG8K9vPm/2bGszPnCN/uOGAXMapr6/jrNzsJC/TlvIgSlRHSlao9tScGXKjWk9J/bN355SVw9IB7KnSNTVgzQx8DpEspM6SUVcAXwEUN2kggRCgH5WCgBGh7/s/oMeDlZ73ZpSRDPZqOvtl9H8t6TYC6mtMKRzfGJcO7MzQ6lOdW7OF4pfGpVt9ak0HSwaM8PXcgvsUp0Lnhb7jGYfQYA0GRrTe7WGp4aoXu8Vij0LsD9X2Hcs376vMqMAA4COwG7pVSnpHvVQhxmxBiixBiS2FhYcsj+/gr+5m1Cv33N9VC6uhbrWvvikSPAWFq0exiMgkeuzCB/KOVvPHLPkNFSMtXhSQuGNKVmfHhru2/3B4weUH8BWqxv7rC9vNPeri4uMuips1Yo9AbC8dqaLg9D9gBdAOGAa8KITqccZKUb0kpR0kpR0VGRlonYexkyN8Nx4ubb3fiiEpmNPgy6ODGVWv8O6iFMCt+xEb26sRFw7rx1m8Z5JS00l+/ASeqarn/650E+3vz1JyByv2zrlorA2cTfyFUlUHmr7afm58IgeEQbN/6rxrnY41CzwXq50qNRs3E63MjsFAq0oFMIN4QCWOnqG3Wb8232/YhVB+HcXcZMqxT6XuOCqqyIn3q32bGYxJwy4dbSDxgw+JxI1TV1HHHJ1tJPFDKs5cMJjzYr16GPq3QnUrsZPDrACnf236uZUFUh8p7PNYo9M1APyFErHmh8yqg4bdqP3A2gBCiMxAHZBgiYbfhqsBzczPW2mplbomdrFIGuDv9ZqjF4H0/t9i0W1gAb1wzkpLyKi56bR3PrdhDRbXtGfpq6yT3fbWDX1MLeeaSwcywlJc7tFutY4T3tblPjYF4+0L/81SkpC3l6epqVWI1/YPcLmhRoUspa4C7gR+AFOArKWWSEOIOIcQd5mb/ACYIIXYDq4C/SSmLDJHQy0ctFKauaDpaLmmRWsUff7chQzqd7iOVR0madV4N0+Ki+Om+KVwyvDtv/LKPWS//dlrdzZaQUvL3bxNZuiuPh2fFn179Jz9JuSt6edv6KTRGEz8byotVcQdrOZylMpfqNZB2gVV+6FLKZVLK/lLKPlLKf5r3zZdSzjf/fVBKOUNKOVhKOUhK+YmhUk64RyXC/+TSM33S9/0Mi/+sZiB9XbhuqC2YvJTZJe1HVZ3FCkIDffjP5UP56KYxVFbXcdn8DTy9JLnF2bqUkmeW7+HzTfv547Q+3Da5z+kN8hOhi57duQR9z1FPSymLrWsvpXpyBT1Dbye4bqRofWLPgis+VoVzP738VMmqxIXw6RUqIvSahapQhKfQbwaUF0FeIxnnmmFy/0h+uG8yV4/pyTtrM5n1UtOz9RNVtdz92XbeWpPBteN68cCMuNMbHMs3lywb3NpPoTESv2BVc3TP0jMKip9BbQ18e5dK7DX6VrXQrvF43EcDxs2Ey96H3C3w2VWw4XX45iZV9eSGJSqhjyfR52xAWG12qU+wnzf/vHgwn94ylsoaNVv/y5c72JRZgjQrgoNHTnDZ/PUsS8zjofPjeeqigWfmOc/frbb6cd11GHAhlOacylfUGFXl8OUfYOdnMPVhmPUfvSDaThCypV96OzFq1Ci5ZcsW20/c/Q0svFUtGvY3K3nfQOMFdAXeOUd9zltbXhxtirLKGv63MpWvtuRQVllD74ggZg/pymebcqioruXlecOYHt/Ej+HaF+Gnx+FvWTpK1FWoroC3p6n1pLs2qIyM9TlxRGUozPkdLnheBdlpPAohxFYp5ajGjrnPDN3C4Mvg8g/grAfgyk88V5mDMrsc2AZlVgRhNUGwnzePXZjApkfO5j+XDaFTkC8v/5xOkJ8Xi+6a0LQyB2U/1yH/roWPP1z8plocXXLf6aaXiqNqnenAVvU/opV5u8M9XRcSLlIvT6ffubD6n7BvFQy9qk1dBfp6c/moHlw+qgd5pSfoGOiLv08LyY7yk/SCqCvSdQhMewhWPQW7v1Y50yvL4LMrIG+HSjsbf4GzpdQ4AfebobcnugyFoKjm7aWtoGtoQMvKvKZSRYlq7wjXZOKfocdYWPoAFKXD51cpM8ul72hl3o7RCt2VMZnULD19lW3BJEZQuEclCdMLoq6JyQsunq/u0fyJkLUW5s6HgRc7WzKNE9EK3dXpdy5UHIEDrVhAbguHzCH/XbTLosvSqTec/6yKlJ7zisodrmnXuKcNvT3RexoIL0j9QZXlcxT5iapAsbtVfWpvjLgOBl3qvumiNYaiZ+iuTkCYymfTmmLZbeHQbpUDXVeJcX20MteY0QrdHeg8UCVYclTMgJTuW7JMo2nHaIXuDkQlwIkSKMt3zHjH8tR4OuRfo3ErtEJ3Byzl3yylxOzNyQVR7bKo0bgTWqG7A1FmhZ7vIIWuc7hoNG6JVujuQFCECjAqSHHMeIcSIawn+Ic6ZjyNRmMIWqG7C1EDHGdyyU/UEaIajRuiFbq70HmgOXrTuoIXraZ4nwr5d6TPu0ajMQSt0N2FqAGqlNiRLPuOs/trQMCgy+w7jkajMRyt0N0FRyyMSgm7vlQVokK7228cjUZjF7RCdxci49XWngujB7ZCSQYM0TlBNBp3RCt0d8EvGMJ6QUGS/cbY9SV4+6syZxqNxu3QCt2dsKQAsAe11ZC4AOLO1+6KGo2bohW6OxE1AIrTVfEJo9n3syprps0tGo3bohW6OxGVoAoaFKUZ3/fOLyCgE/Q9x/i+NRqNQ9AK3Z2weLoYbXapOAp7l6m82l4+xvat0Wgchlbo7kR4XzB5274wWl6ilHZTpCyGmgptbtFo3Bxdscid8PaFiP7Wz9AL98K6l5T3ijBB76nKgyVuFnj5qn7yE2HT29AxFqJH2VV8jUZjX7RCdzeiBkDu5ubbHNgGvz0Pe5aoMnKjblamlJTFkLYSuAeoVyzDrwNc8DwIYU/JNRqNndEK3d2ISlDuhZXHwC/k9GNlhfDTE7DjE/APg8n/B2PvgKBwdXzG06q0XOoKNWPvPEi5QoZGa2Wu0XgAWqG7GycXRvdAj9Hq79oa2Pw2rH5G5XuZcA9M/iv4dzj9XCGg6xD10mg0HodW6O7GyepFSeDlDUnfQtIiOJINfabDzOcgsr9TRdRoNM5BK3R3I7Qn+ATBsr9CbZXyeomdAjOfUYud2nSi0bRbtEJ3N0wmGHkDFKfBgDkQfwEEdnK2VBqNxgXQCt0dmfkvZ0ug0WhcEB1YpNFoNB6CVQpdCDFTCLFXCJEuhHiwiTZThRA7hBBJQohfjRVTo9FoNC3RoslFCOEFvAacC+QCm4UQ30spk+u1CQNeB2ZKKfcLIaLsJK9Go9FomsCaGfoYIF1KmSGlrAK+AC5q0OZqYKGUcj+AlLLAWDE1Go1G0xLWKPTuQE6997nmffXpD3QUQvwihNgqhLiusY6EELcJIbYIIbYUFha2TmKNRqPRNIo1Cr0xx2bZ4L03MBK4ADgPeFQIcUZ0i5TyLSnlKCnlqMjISJuF1Wg0Gk3TWOO2mAv0qPc+GjjYSJsiKeVx4LgQYg0wFEg1REqNRqPRtIg1M/TNQD8hRKwQwhe4Cvi+QZvvgLOEEN5CiEBgLGDH8vQajUajaUiLM3QpZY0Q4m7gB8ALeE9KmSSEuMN8fL6UMkUIsQLYBdQB70gpE5vrd+vWrWVCiL31doUCpfq91e8jgCIXkke/N+69vrft672t9zuOppBSOuUFbGnw/i393qb3+vp56Ht9b9vde1vv92nt679cKVJ0sX5v0/uGOFse/d7Y964ki35v3/cNsbX9SYRZ4zscIcQWKaWuedZK9PXzXPS9bV/Yer+ba+/MGfpbThzbE9DXz3PR97Z9Yev9brK902boGo1GozEWV7Khn0QI8Z4QokAIkVhv3xNCiAPmBGA7hBCzDBqrhxBitRAixZxY7F7z/v8IIfYIIXYJIRaZ89XYc7yhQogNQojdQojFQogOLfVlw5hnJFez1/U0993Y/bPX9WxsLHtey0bvn/nYn8zXOUkI8W8Dx2zs/v3DfC13CCFWCiG6GTRWY9ezkxDiRyFEmnnb0Yixmhnvy3rfyywhxA4Dx2vy/pmPPyCEkEKICAPG8hdCbBJC7DSP9aR5/+Xm93VCCGNNa02tljrzBUwGRgCJ9fY9ATxgh7G6AiPMf4eggqESgBmAt3n/c8Bzdh5vMzDFvP8m4B8GjecF7AN6A77ATvN4drmezdw/e13Pxsayy7Vs4f5NA34C/MzHoux8/zrUa3MPMN+O1/PfwIPmvx806t41NV6D488Dj9n7/pnf90C5Z2cDEQaMJYBg898+wO/AOGAAyvXwF2CUUZ9NStfycjmJlHINUOKgsfKklNvMfx9DBUR1l1KulFLWmJttREXI2m081A1eY272I3CpEeNhXXI1Q2ns/tnxejb2XbHXtWzu/t0JPCulrDQfMypBXaP3T0p5tF6bIM5Mx9EqmrieFwEfmv/+EJhrxFjNjAeAEEIAVwCfGzheU/cP4AXg/zDuWkopZZn5rY/5JaWUKVLKvc2c2mpcUqE3w93mx8z3jHzssyCEiAGGo35J63MTsNzO4yUCc8yHLuf0dAttobnkana9ns1gl+tZD3tdy9NocP/6o6KlfxdC/CqEGG3QME3ePyHEP4UQOcAfgMcMGq8xOksp80ApRMBR6bHPAvKllGn26Lz+/RNCzAEOSCl3GjyGl9lkVAD8KKVsqFsMxZ0U+htAH2AYkId6FDMMIUQwsAD4c/3ZjxDiEaAG+NTO490E/FEIsRX1KFhl1FCN7JPY+Xo2KYydrmcD7HUtT9LI/fMGOqIeqf8KfGWeYbZ5qEb2SQAp5SNSyh6oa3m3AWO5GvMwcHZen/r3D/V9fAQ7/ChKKWullMNQT6RjhBCDjB6jPm6j0KWU+eaLUwe8jXoUNQQhhA/q5n4qpVxYb//1wGzgD9JsCLPXeFLKPVLKGVLKkagv8T6Dhms0uZo9r2dT2Ot6NsSO1xJo8vuSi6oJIKWUm1ApMNq8sIZ1yfE+w0CzUiPkCyG6Api3dq93IITwBi4BvrRD3w3vXx8gFtgphMhCXeNtQoguRo0ppTyCspnPNKrPxnAbhW75Qpm5GPVYbUS/AngXSJFS/q/e/pnA34A5UspyI8ZqYbwo89YE/B2Yb9CQjSZXs9f1bAp7Xc8mxrLXtWzy/gHfAtPNbfqjFjCLzujAdpq6f/3qtZkD7DFgrKb4Hrje/Pf1qGR89uYcYI+UMtfIThu7f1LK3VLKKClljJQyBvUjOkJKeaiNY0UKszeXECIA82dqS58tYuQKq1Ev1KwqD6hGXdybgY+B3agEYN8DXQ0aaxLqEXYXsMP8mgWko2yXln1GeRE0Nd69qBX3VOBZzDECBo05y9zvPuAR8z67XM9m7p+9rmdjY9nzWjZ1/3yBT1A/jNuA6Xa+fwvMY+1ChYZ3t+P1DAdWAWnmbSd7flfM+z8A7jBqnJbuX4M2WRjj5TIE2G4eKxGztw5qApULVAL5wA9GfT4dWKTRaDQegtuYXDQajUbTPFqhazQajYegFbpGo9F4CFqhazQajYegFbpGo9F4CFqhazQajYegFbpGo9F4CFqhazQajYegFbpGo9F4CFqhazQajYegFbpGo9F4CFqhazQajYegFbpGo9F4CFqhazQajYegFbpGo9F4CFqhazQajYegFboLI4S4WAghhRDxzpZFYwzm+/lxvffeQohCIcQSZ8qlsT9CiDJ7j6EVumszD1iLqiNpNUIIL/uIozGA48Agc41JgHOBA06UR+NBaIXuogghgoGJqJqOV5n3TRVCrBFCLBJCJAsh5psLISOEKBNCPCWE+B0Y7zzJNVawHLjA/Pc8VF1NAIQQY4QQ64UQ283bOPP+34QQw+q1WyeEGOJIoTVtx/w/vKTe+1eFEDeY/84SQjwphNgmhNjdmidzrdBdl7nACillKlAihBhh3j8GuB8YDPQBLjHvDwISpZRjpZRrHS2sxia+AK4SQvijCgn/Xu/YHmCylHI48BjwL/P+d4AbAIQQ/QE/KeUuh0mscRRFUsoRwBvAA7aerBW66zIP9Y+PeTvP/PcmKWWGlLIWNbObZN5fi6oEr3FxzIo4BnVPlzU4HAp8LYRIBF4ABpr3fw3MFkL4ADcBHzhEWI2jWWjebkV9R2zC21BRNIYghAgHpqNsrRLwAiTqn182aG55X2FW8hr34Hvgv8BUILze/n8Aq6WUFwshYoBfAKSU5UKIH4GLgCuAUY4UVmMYNZw+kfZvcLzSvK2lFfpZz9Bdk8uAj6SUvaSUMVLKHkAmajY+RggRa7adX4laNNW4H+8BT0kpdzfYH8qpRdIbGhx7B3gZ2CylLLGveBo7kQ0kCCH8hBChwNlGdq4VumsyD1jUYN8C4GpgA/AskIhS8g3badwAKWWulPKlRg79G3hGCLEO9WRW/5ytwFHgfQeIqDEQIYQ3UCmlzAG+AnYBnwLbDR1HyoZP8BpXRQgxFXhASjnbyaJonIAQohvKBBMvpaxzsjgaGxBCDAXellKOsec4eoau0bgBQojrUN4wj2hl7l4IIe5AOTD83e5j6Rm6RqPReAZ6hu4iCCF6CCFWCyFShBBJQoh7zfs7CSF+FEKkmbcdzfvPFUJsNQcgbBVCTK/X1z+FEDmOCDXWaDSug56huwhCiK5AVynlNiFECMoPdS7K06FESvmsEOJBoKOU8m9CiOFAvpTyoBBiEPCDlLK7ua9xqNX0NCllsDM+j0ajcTxaobsoQojvgFfNr6lSyjyz0v9FShnXoK0AioBuUsrKevvLtELXaNoP2uTigpgDSoajFsE6SynzAMzbqEZOuRTYXl+ZazSa9oeOFHUxzEm5FgB/llIeVZPvZtsPBJ4DZjhAPI1G48LoGboLYc7TsQD4VEppyemQbza1WOzsBfXaR6MCi66TUu5ztLwajca10ArdRTDbwd8FUqSU/6t36HvgevPf1wPfmduHAUuBh6SU6xwoqkajcVH0oqiLIISYBPwG7AYsgSMPo+zoXwE9gf3A5VLKEiHE34GHgLR63cyQUhYIIf6NShPQDTgIvCOlfMIhH0Sj0TgNrdA1Go3GQ9AmF41Go/EQtELXaDQaD0ErdI1Go/EQtELXaDQaD0ErdI1Go/EQtELXtFuEEE8IIZqsrC6EmCuESHCkTBpNW9AKXaNpmrmAVugat0H7oWvaFUKIR4DrgBygEJWmuBS4DfAF0oFrgWHAEvOxUlQCNIDXgEigHLhVSrnHgeJrNM2iFbqm3SCEGAl8AIxFJabbBswH3pdSFpvbPI3KM/+KEOIDYImU8hvzsVXAHVLKNCHEWOAZKeX0M0fSaJyDzraoaU+cBSySUpYDCCG+N+8fZFbkYUAw8EPDE81ZMCcAX9fLgOlnb4E1GlvQCl3T3mjskfQDYK6UcqcQ4gZgaiNtTMARKeUwu0mm0bQRvSiqaU+sAS4WQgSYy/xdaN4fAuSZ0xf/oV77Y+ZjSCmPAplCiMtBZccUQgx1nOgaTctoG7qmXVFvUTQbyAWSgePA/5n37QZCpJQ3CCEmAm8DlcBlqCyYbwBdAR/gCynlUw7/EBpNE2iFrtFoNB6CNrloNBqNh6AVukaj0XgIWqFrNBqNh6AVukaj0XgIWqFrNBqNh6AVukaj0XgIWqFrNBqNh6AVukaj0XgI/w/50WYgH0bZwQAAAABJRU5ErkJggg==\n"
},
"metadata": {
"needs_background": "light"
}
}
]
}
],
"metadata": {
"_draft": {
"nbviewer_url": "https://gist.github.com/f132aa1703ecded98967aadbed26efa4"
},
"gist": {
"id": "f132aa1703ecded98967aadbed26efa4",
"data": {
"description": "England and York",
"public": true
}
},
"kernelspec": {
"name": "python3",
"display_name": "Python 3",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.9.2+",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment