Skip to content

Instantly share code, notes, and snippets.

@jph00
Last active November 1, 2021 05:04
Show Gist options
  • Save jph00/0a7273c3b76f3ceceda70eb0efaf4f43 to your computer and use it in GitHub Desktop.
Save jph00/0a7273c3b76f3ceceda70eb0efaf4f43 to your computer and use it in GitHub Desktop.
nbs/uk-young-covid.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"id": "03a380d9",
"cell_type": "code",
"source": "import requests, pandas as pd, matplotlib as mpl, matplotlib.pyplot as plt",
"execution_count": 2,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"id": "b2a48a1e",
"cell_type": "code",
"source": "endpoint = ('https://api.coronavirus.data.gov.uk/v1/data?'\n 'filters=areaType=nation;areaName=england&'\n 'structure={\"date\":\"date\",\"newCases\":\"newCasesBySpecimenDateAgeDemographics\"}')\ndata = requests.get(endpoint, timeout=20).json()['data']",
"execution_count": 3,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"id": "9bcf2046",
"cell_type": "code",
"source": "rows = [(o['date'],p['age'],p['cases']) for o in data for p in o['newCases']]\ndf = pd.DataFrame(rows, columns=['date','age','cases'])\ndf['date']= pd.to_datetime(df['date'])",
"execution_count": 4,
"outputs": []
},
{
"metadata": {
"scrolled": false,
"trusted": true
},
"id": "8c29237c",
"cell_type": "code",
"source": "df = df.query(\"age not in ['00_59','60+'] & '2020-03-01'<=date<'2021-10-23'\")\nyoung = df.query(\"age in ['00_04','05_09','10_14','15_19']\"\n ).groupby(pd.Grouper(key='date', freq='W-FRI')).sum()\n\nyoung.plot(y='cases', legend=False, title='England Weekly Cases, Age<20', figsize=(10,6))\nplt.gca().yaxis.set_major_formatter(mpl.ticker.StrMethodFormatter('{x:,.0f}'))",
"execution_count": 5,
"outputs": [
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 720x432 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAmwAAAGPCAYAAADhrYyUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABUx0lEQVR4nO3deXxcdb3/8ddnMtmXJmnT0ibdN9pCC6W07FaQTVRAQXEDFcEFV7hXwevv4kO9CnpdAfUiICAqICIgIIil7FBoge6le7N0y76vM9/fH+ekTNMkTbPNZOb9fDzmkcz3nO85n5mk00++qznnEBEREZHYFYh2ACIiIiLSOyVsIiIiIjFOCZuIiIhIjFPCJiIiIhLjlLCJiIiIxDglbCIiIiIxTgmbiBzCzKaYmTOz4BBc+zNm9tJgX/cIY1hqZqVHekxEJFqUsImMYGa208yazawh4nFrtOPqLzN72sy+FfG80E8cuys7KjpR9s7MPmFmK/2fxR4z+6eZnRbtuHpiZlPNLGxmv4nCvcea2V/MbLeZ1ZrZy2a2pMs5nzCzXWbWaGaPmFn+cMcpEguUsImMfB90zmVFPL4S7YAG4AXgPRHPzwA2dVO2xTm3dzgD6wszuxb4JfAjYBwwCfgNcGEUwzqcy4Fq4DIzSx2OG5rZOP/bLOAN4AQgH7gHeMLMsvzz5gH/B3wa7/1swns/RRKOEjaRONXZ9Whm/2tm1Wa2w8zOjzg+1cxeMLN6M/u3md1mZvf1cK3PmtlG/9ztZvaFiGNLzazUzK4zs/1+q9JnI46PNrPHzKzOzF4HpvcS9gvAqWbW+dl0Ol4CtKhL2Qv+tY82s2fMrMrM3jGzj0bcN9V/7cVmts/Mfmdm6T28vq+Z2QYzK+pS/p9m9rcuZbeY2S+7ucYo4PvANc65h51zjc65dufcP5xz/+mfs9jMXjWzGv99utXMUvxjZma/8N/DWjNbY2bHHO61mNkYM3vcv2aVmb0Y8V71xeXAd4F24INdXtM5/vtaa2a/MbPnzezzEcc/5/9eVPuto5N7uomZHeW/nxuA7wE457Y7537unNvjnAs5524HUoDZfrVPAv9wzr3gnGsA/h/wYTPLPoLXJxIXlLCJxLclwDvAGOAnwJ1mZv6xPwOvA6Px/gP9dC/X2Q98AMgBPgv8wswWRhw/ChgFFAJXAreZWZ5/7DagBRgPfM5/9OR1IBVY4D8/A3gG2Nql7AUzy/SP/RkYC3wc+I3fKgNwMzALOA6Y4cf2311vaGb/D/gM8B7nXNexa/cB55lZrn9uEPgY8MduYj8ZSAP+3svrCwHfxPt5nAycBXzZP3aO/9pmAbn+fSr78FquA0qBArxWqO8Afdpz0MxOB4qA+4EH8ZK3zmNjgIeAG/B+R94BTok4fpF/rw/7934R+EuX6yeb2YfN7B9+/fnA14BreojnOLyEbatfNA9Y3XncObcNaPPfC5HE4pzTQw89RugD2Ak0ADURj6v8Y58Btkacm4H3H/lReF11HUBGxPH7gPv876f45wZ7uO8jwNf975cCzZHn4iV4JwFJeC03R0cc+xHwUi+v6Tng63hdZKV+2U0RZWFgMl5C82KXuv8H3AgY0AhMjzh2MrAjIuYy4OfAS8CoiPOWdt7Xf/7PiPf0A8CGHuL+JLD3CH9+3wD+7n9/JrDZf98CEecc7rV8H3gUmNGP3587gEcirtkOjPWfXw682iWOEuDzEe/LlRHHA3hdlpP95z/wfw9ewEvSsw8TSw6wFrghomwZ8MUu55UBS6P9b08PPYb7oRY2kZHvIudcbsTj9xHHDozzcs41+d9mAROAqogy8P4z7paZnW9mr/ldbjXA+/FaiTpVOuc6Ip43+fcpAIJdrr3rMK/nBbyWptPxkin8r51lJc65XXhJ2xK/K7DGj+uTeAlpAV6Cuiri2FN+eadc4Grgx8652l7iuQf4lP/9p+i+dQ281rAx1svMWjOb5Xdf7jWzOrzkdQyAc+5Z4Fa8Fsl9Zna7meX04bX8FK9F6l9+d/X1vbyWyFjSgUuBP/n3fxUoBj7hnzKBiJ+bc87hteR1mgz8KiKmKrykrtA/PhtIBt4G1jjn6g8Tyz+A15xzP4441ICXyEXKAXq8lki8UsImkpj2APlmlhFRNrG7E80biP434H+Bcc65XOBJvP+cD6ccryUv8tqTDlPnBbzE7Ay8bjaAl4FT/bIX/LIS4PkuyWqWc+5LQAVeq9+8iGOjnHNZEfepxmsx+4OZndpLPI8A8/3xZB/AT3C68Spe1+9FvVzrt3iTKGY653LwuhQPvI/OuV87507A6wqcBfzn4V6Lc67eOXedc24a3hi0a83srF5i6HQxXvLzGz+B3IuXbHV2i+7B6y4FvDF2kc/x3v8vdHn/051zr/hxfRSvC7cCeMDM1pvZ9d2ME0zFe4/LgC9wsPW82xWOmU3D6zLf3IfXJxJXlLCJJCC/hWol8D0zSzGzk+ky4DxCCt5/kuVAh3kTF87p431CwMP+fTLMbC5wxWGqvYLX+vUp/ITNOVft3/9TvJuwPQ7MMrNP+2Olks3sRDOb45wLA7/HG2s3Fg4sB3Jul/iew2uV+7t1WU4i4pwWvLFcfwZed84V93BeLd64stvM7CL/9Sb7rZM/8U/LBuqABjM7GvhSZ30/9iVmlozXBdoChA73WszsA2Y2w0+o6vDGyYX8Y3eb2d09vM9XAHcBx+IlVsfhJcXHmdmxwBPAsf5rCeKNO4tcSuV3wA2dYwbNbJSZXdrlPdnlnPs+3ri7L+O1uq03s+/5dZL997YZuNx/rZH+BHzQzE73xyx+H3i4t9Y6kXilhE1k5PuHHbwOW2+D3iN9Em/cUiXwQ+ABoLXrSf5/jl/DG5Rejddl9tgRxPcVvO7RvcDdwB96O9nvpl2FlySuizj0It7kghci4joHuAzY7V//Zr8ewLfxugpf87sf/827sw8j7/cM3kSKx8zshB7CugcvsempO7TzWj8HrsWbdVmO1wr1FbwWJID/wHv/6vGSsAciquf4ZdV43caVeK2ah3stM/3nDXitfL/xE1HwWjZf7hqnmRXiTXj4pXNub8RjFV536xXOuQq8LtOf+LHMxUvyW/3X+ne89/t+P6Z1wPld7+Wf65xzzzvnPovX1dr5fpyC12p5DlAT8Tt8ul9vPfBFvMRtP17C++Wu1xdJBOYNSxCRRGdmDwCbnHM3RjuWWGNmk/C6Mo9yztVFO56+MG+5kNXAfOdc+yBcL4A3hu2TzrnlA72eiBwZtbCJJCi/C266mQXM7Dy8xV0fiXJYMcdPVK4F7h8pyRqAc67NOTdnIMmamZ1rZrn+OLPO8XavDVqQItJng75PoIiMGEfhjS8bjddy8iXn3FvRDSm2+OOm9uF1UZ4X5XCi4WS8sXspwAa8GcnN0Q1JJDGpS1REREQkxqlLVERERCTGKWETERERiXFxP4ZtzJgxbsqUKdEOQ0REROSwVq1aVeGcK+haHvcJ25QpU1i5cmW0wxARERE5LDPrdvs+dYmKiIiIxDglbCIiIiIxTgmbiIiISIxTwiYiIiIS45SwiYiIiMQ4JWwiIiIiMU4Jm4iIiEiMO2zCZmZ3mdl+M1sXUfZTM9tkZmvM7O9mlhtx7AYz22pm75jZuT1cc6qZrTCzLWb2gJml+OVmZr/2668xs4URdc7zr7nVzK4f0KsWERERGUH60sJ2N3Bel7JngGOcc/OBzcANAGY2F7gMmOfX+Y2ZJXVzzZuBXzjnZgLVwJV++fnATP9xNfBb/7pJwG3+8bnAx/17iYiIiMS9wyZszrkXgKouZf9yznX4T18DivzvLwTud861Oud2AFuBxZF1zcyAM4GH/KJ7gIsi6t/rPK8BuWY23r/GVufcdudcG3C/f66IiIhI3BuMMWyfA/7pf18IlEQcK/XLIo0GaiISvshzeqrfl+uKiIiIxKUBJWxm9l9AB/CnzqJuTnNdq/VyTk/H+nLdyLiuNrOVZrayvLy8p9NEREREYkZzW6jHY/1O2MzsCuADwCedc53JUykwMeK0ImB3l6oVeF2dwW7O6al+X657gHPudufcIufcooKCQza8FxEREYkpL2+t4PSfLO/xeL8SNjM7D/g28CHnXFPEoceAy8ws1cym4k0eeN2vc6+ZLfaTu+XAJX6dK4BHI+pf7s8WPQmodc7tAd4AZvqzS1PwJjY81p/YRURERGKFc47fPreNT9+5gryM5B7P68uyHn8BXgVmm1mpmV0J3ApkA8+Y2dtm9jv/puuBB4ENwFPANc65zva9+cAe//tvA9ea2Va8MW13+uVPAtvxJiv8Hviyf90O4CvA08BG4EH/XiIiIiIjUkNrB1/+05vc/NQmzj9mPI9cc2qP59q7vZlDx8xygDudc5cO+c26WLRokVu5cuVw31ZERESkR1v3N/CFP65kZ2UT1593NJ8/fSpmhpmtcs4t6np+sLuLDDbnXB0w7MmaiIiISKx5at0erntwNWnJSdx35RJOnj76sHWGJWETEREREbj75R187x8bWDAxl999aiHjR6X3qZ4SNhEREZFh8uDKUhZMzOXBL5xEarC7zaC6p83fRURERIZJaXUTC4pGHVGyBkrYRERERIZFbXM7dS0dFOX1rRs0khI2ERERkWFQWu0tXTsxL+OI6yphExERERkGJVXNABQpYRMRERGJTQda2PLVJSoiIiISk0qrm8lKDTIqvectqHqihE1ERERkGJRWN1GUl46ZHXFdJWwiIiIiw6C0urlf49dACZuIiIjIkHPOUVLV1K8lPUAJm4iIiMiQq2lqp7EtxMR8tbCJiIiIxKQSf4aoWthEREREYlRpdecabErYRERERGJS6YEWNnWJioiIiMSkkqpmctL6twYbKGETERERGXKl1U39nnAASthEREREhlxJdXO/x6+BEjYRERGRIeWc83c5UAubiIiISEyqbGyjpT3MRLWwiYiIiMSmkqqBzRAFJWwiIiIiQ+rAGmz5amETERERiUklA1yDDZSwiYiIiAyp0upm8jKSyUoN9vsaSthEREREhlBpdfOA1mADJWwiIiIiQ6q0qmlAa7CBEjYRERGRIRMOO0prmgc0fg2UsImIiIgMmfKGVto6BrYGGyhhExERERkypYMwQxSUsImIiIgMmc412CYOYA02UMImIiIiMmQ6dzkozFULm4iIiEhMKq1uZkxWCukpSQO6jhI2ERERkSFSUt004PFroIRNREREZMiUVjcPeA02UMImIiIiMiRCYcfumoHvcgBK2ERERESGxL66FtpDTi1sIiIiIrGqc0kPjWETERERiVGdS3oMdJcDUMImIiIiMiQ6W9gm5CphExEREYlJpdVNjM1OJS15YGuwQR8SNjO7y8z2m9m6iLJ8M3vGzLb4X/Mijt1gZlvN7B0zO7eHa041sxV+/QfMLMUvNzP7tV9/jZktjKhznn/NrWZ2/cBetoiIiMjQKqluGpQZotC3Fra7gfO6lF0PLHPOzQSW+c8xs7nAZcA8v85vzKy7tPJm4Bd+/WrgSr/8fGCm/7ga+K1/3STgNv/4XODj/r1EREREYtJgrcEGfUjYnHMvAFVdii8E7vG/vwe4KKL8fudcq3NuB7AVWBxZ0cwMOBN4qIf69zrPa0CumY33r7HVObfdOdcG3O+fKyIiIhJzOkJh9tS2MHEQZohC/8ewjXPO7QHwv471ywuBkojzSv2ySKOBGudcRzfn9FS/L9cVERERiQl7alsIhQdnDTYY/EkH1k2ZO4JzejrWl+u+ewOzq81spZmtLC8v7+k0ERERkSFRUu0t6TEYa7BB/xO2fX5XJf7X/X55KTAx4rwiYHeXuhV4XZ3Bbs7pqX5frnuAc+5259wi59yigoKCPr8oERERkcHQuaTHxPzotrA9Blzhf38F8GhE+WVmlmpmU/EmD7wOYGb3mtli55wDlgOX9FD/cn+26ElArd/l+gYw059dmoI3seGxfsYuIiIiMqRKq5sxg/GjBidhCx7uBDP7C7AUGGNmpcCNwE3Ag2Z2JVAMXArgnFtvZg8CG4AO4BrnXMi/1Hxgj//9t4H7zeyHwFvAnX75k8D78SYrNAGf9a/bYWZfAZ4GkoC7nHPrB/C6RURERIZMaVUT43PSSAkOzuizwyZszrmP93DorB7O/x/gfyLLzCwH2OKcK/HP2U6X2aN+uQOu6eG6T+IldCIiIiIxzVvSY3DGr8Ew7XTgnKtzzl06HPcSERERibaS6qZBmyEK2ppKREREZFC1dYTZW9dC0SDtcgBK2EREREQG1Z7aZpxDLWwiIiIisWrzvgZACZuIiIhITGppD/Hjf25kwqg0jpuYO2jXPewsURERERHpm/99+h22lzdy35VLyEgZvDRLLWwiIiIig+D1HVXc+fIOPnXSJE6bOWZQr62ETURERGSAmto6+M+HVjMxL4Mbzp8z6NdXwiYiIiLSg5e3VnDmz55j2cZ9vZ530z83UVzVxE8vmU9m6uCPOFPCJiIiItKD13dUsb28kSvvWckPHt9AW0f4kHNe3lrBva/u4rOnTGXJtNFDEocSNhEREZEeVDW2kZMW5DOnTOHOl3bwkd++ws6KxgPH61va+dZDa5g2JpNvnTd7yOJQwiYiIiLSg8rGVgqyU/neh+Zx+6dPoLiqiQ/c8hKPvl0GwA8f38ie2mb+96MLSEtOGrI4tKyHiIiISA8qG9oYnZkKwDnzjmJe4Si+/pe3+Pr9b/Pwm2U8v7mcL75nOgsn5Q1pHGphExEREelBZWMbo7NSDjwvzE3n/qtP4qtnzuCFLeXMGpfFN8+eOeRxqIVNREREpAdVjW3kZ6YcVBZMCnDdObM5/5jxjMlOITU4dF2hB+455HcQERERGYFCYUd1Uxujs1K7PT53Qs6wxaIuUREREZFuVDe14RyM7tLCFg1K2ERERES6UdnQBnDQGLZoUcImIiIi0o3KxlaAQ8awRYMSNhEREZFudLawjelhDNtwUsImIiIi0o2qRi9hUwubiIiISIyqbGjFDPIylLCJiIiIxKTKxjbyMlJICli0Q1HCJiIiItIdb1uq6LeugRI2ERERkW51t8tBtChhExEREelGRWNrTMwQBSVsIiIiIt1SC5uIiIhIDGsPhalpao+JXQ5ACZuIiIjIIaqb/G2p1MImIiIiEpve3UdUY9hEREREYlIs7XIASthEREREDlHR4G38PkZj2ERERERi07stbOoSFREREYlJlQ1tBAxy05OjHQqghE1ERETkEJX+GmyBGNhHFJSwiYiIiByisqGV0THSHQpK2EREREQOEUu7HIASNhEREZFDVDa2xcwuB6CETUREROQQXpeoEjYRERGJYc457nxpB/vqWqIdyrBr6whT19IRM7scgBI2ERER6UZpdTM/eHwDD75REu1Qhl3nPqJxM4bNzL5pZuvNbJ2Z/cXM0sws38yeMbMt/te8Hupe4Z+zxcyuiCifamYr/PIHzCzFLzcz+7WZbTWzNWa2cCCxi4iISM9KqpsA2F7RGOVIhl+s7XIAA0jYzKwQ+BqwyDl3DJAEXAZcDyxzzs0ElvnPu9bNB24ElgCLgRsjErubgV/49auBK/3y84GZ/uNq4Lf9jV1ERER6V1bdDMD28oYoRzL8Ym2XAxh4l2gQSDezIJAB7AYuBO7xj98DXNRNvXOBZ5xzVc65auAZ4DwzM+BM4KFu6l8I3Os8rwG5ZjZ+gPGLiIhIN0r9hG1beSPOuShHM7wqG7yELS5miTrnyoD/BYqBPUCtc+5fwDjn3B7/nD3A2G6qFwKRneKlftlooMY519GlvLc6IiIiMsjKaryEraG1g/31rVGOZnhV+i1scTFL1O/CvBCYCkwAMs3sU32t3k2Z66W8tzrdxXa1ma00s5Xl5eV9DElEREQ6lVY3EfS3ZdqWYN2ilQ2tBANGTlps7CMKA+sSfR+wwzlX7pxrBx4GTgH2dXZV+l/3d1O3FJgY8bwIrzu1Aq+rM9ilvLc6h3DO3e6cW+ScW1RQUNCvFyciIpLIymqaWTjZG16+rTyxJh5UNbaRF0P7iMLAErZi4CQzy/DHnp0FbAQeAzpnfV4BPAreJAUzW+aXPw2cY2Z5fkvdOcDTzuskXw5c0rW+f93L/dmiJ+F1we4ZQPwiIiLSjVDYsaemhUWT88hISWLb/sRqYatoaIup7lDwJg30i3NuhZk9BLwJdABvAbcDWcCDZnYlXlJ3qV9lvH8ezrkqM/sB8IZ/7PvOuSr/+28D95vZD/1r3umXPwm8H9gKNAGf7W/sIiIi0rN9dS10hB1FeRlMK8hMuKU9qhpbY2rCAQwgYQNwzt2ItzxHpFa81rauTgJui6h7F3BXN9fcjrfUR9dyB1wzkHhFRETk8DpniBblpTO9IIuVO6ujHNHwqmxsY0FebrTDOMiw7XTgnLvVOffYcN1PRERE+qesxls0t9BP2MpqmmluC0U5quFT1dAWU7scgLamEhERkS5Kq7wWtsLcdKYVZAKwvSIxxrG1doSob+2IqV0OQAmbiIiIdFFW08yYrFTSkpOYXpAFwPYEmSkai7scgBI2ERER6aKsppmivHQApo7JxCxx1mKLxV0OQAmbiIiIdFFa3Uyhn7ClJSdRlJeeMGuxxeIuB6CETURERCKEw+6gFjaAaWOyEmYT+MoGbxuu0VnqEhUREZEYVdHQSltHmKLcdxO26QVZbC9vJByO/03g3x3DphY2ERERiVGl/qbvhREtbNPHZtLcHmJPXUu0who2FQ1tJCcZOWkDWqp20ClhExERkQPeXTQ340BZ50zRRNiiqqqxlfzMFLxdN2OHEjYRERE5oKz63TXYOh1Yiy0BxrFVNrQxOsaW9AAlbCIiIhKhtLqJvIxkMlPf7RIsyEolOy2YEDNFKxvbYm5JD1DCJiIiIhHKapoPGr8GYGZML8hKiLXYKhtbY25JD1DCJiIiIhFKq5spys04pHxaQWZC7Hbg7SOqLlERERGJUc45yqoPbWEDb+LB3roWGlo7ohDZ8GhpD9HYFlKXqIiIiMSuqsY2mttDBy2a2+ndPUXjt1s0Vnc5ACVsIiIi4iurOXSGaKcZY72ZovE8ji1WdzkAJWwiIiLiK+tmDbZOk/IzSQpYXI9jq4zRXQ5ACZuIiIj4OhfN7W4MW0owwKT8jDhvYfMStjEawyYiIiKxqqymmezUIKPSk7s9Pr0gk23747eFrarR6xJVC5uIiIjErNLqpm5b1zpNK8hiR2UjoTjdBL6yoY2UYICs1NjaRxSUsImIiIivtLq52xminaYXZNLWET4w1i3eVDa2MToG9xEFJWwiIiLiK6tu7nbCQacDm8DH6Ti2yobWmFyDDZSwiYiICFDb3E59a0e3S3p0iveEraoxNnc5ACVsIiIigjd+Dei1SzQvM4W8jOS43QS+oqEtJhfNBSVsIiIiwrtrsPU26QCI603gqxqVsImIiEgMK+1l0dxI0wuy4nJ7qqa2DprbQ+RrDJuIiIjEqrKaZtKTk8jL6H4Ntk7TCjKpaGijtql9mCIbHgcWzdUYNhEREYlVpdVNFOWlH3ZJiwMTDyriq5UtlrelAiVsIiIigtfCdrjxawDTx/oJ2/74Stg6dznQsh4iIiISsw63aG6niXnpJCcZ2yvia6Zohd8lOlpdoiIiIhKLGls7qGlqpzC39wkHAMGkAJNHZ8ZhC5ufsKmFTURERGJRWU3flvToNL0gM+6W9qhqbCM1GCAjJSnaoXRLCZuIiEiC68uiuZGmFWSxq7KJjlB4KMMaVhUNrYzJSo3JfURBCZuIiEjC61w0t6iXbakiTRmdQUfYsae2ZSjDGlbetlSx2R0KSthEREQSXml1MynBAGOy+jbgflJ+JgC7KpuGMqxhVdnQFrPj10AJm4iISMIrrWmmMDedQKBv3YGTRnuTE4qr4idhUwubiIiIxLS+LunR6aicNFKSAuyqio+lPVraQ5T7Y9hilRI2ERGRBFdW7bWw9VVSwCjKT6c4DrpEnXPc8PBa2kNhls4qiHY4PVLCJiIiksBa2kNUNLQeUQsbwOT8jLgYw3b3Kzv5+1tlfOOsWZwyY0y0w+mREjYREZEEdqRrsHWalJ9BcVUTzrmhCGtYrNheyQ+f2Mj75ozjq2fOiHY4vRpQwmZmuWb2kJltMrONZnaymeWb2TNmtsX/mtdD3Sv8c7aY2RUR5VPNbIVf/oCZpfjlZma/NrOtZrbGzBYOJHYRERHxxq8BFOUdfpeDSJNGZ9LQ2kF1U/tQhDXk9tQ2c82f32RyfgY//9iCPk+4iJaBtrD9CnjKOXc0sADYCFwPLHPOzQSW+c8PYmb5wI3AEmAxcGNEYncz8Au/fjVwpV9+PjDTf1wN/HaAsYuIiCS8zjXYjmQMG3hdogC7KkfexIPWjhBfvO9NmttC3H75CeSkJUc7pMPqd8JmZjnAGcCdAM65NudcDXAhcI9/2j3ARd1UPxd4xjlX5ZyrBp4BzjNveeEzgYe6qX8hcK/zvAbkmtn4/sYvIiIi3i4HwYAxLiftiOpNHsFLe9z46HpWl9Tws48uYMbY7GiH0ycDaWGbBpQDfzCzt8zsDjPLBMY55/YA+F/HdlO3ECiJeF7ql40GapxzHV3Ke6sjIiIi/VRc1cT43DSSjrBLcOKBFraRlbD9eUUx979RwjXvnc55x4ycdp+BJGxBYCHwW+fc8UAj3XR/9qC73wrXS3lvdQ69uNnVZrbSzFaWl5f3MSQREZHEs3FPHUcflXPE9dKSkxiXkzoiWtjCYceWffX88dWd3PjYOs6YVcC1Z8+OdlhHJDiAuqVAqXNuhf/8IbyEbZ+ZjXfO7fG7LPf3UHdpxPMi4DmgAq+rM+i3shUBuyPqTOxSZzfdcM7dDtwOsGjRopE7fUVERGQINbZ2sL2ikQ8t6F+H1eT8zJhci622uZ23iqt5q7iGN4urebukhvoWr/Nu1rgsfn3ZcUfcohht/U7YnHN7zazEzGY7594BzgI2+I8rgJv8r48CmFkh3hi0s4CngR9FTDQ4B7jBOefMbDlwCXB/ZH3gMeArZnY/3mSF2s6uVxERETlyG/fU4RzMm3DkLWzgbVH14pbY6skqr2/ljJ8sp7k9RMBg1rhsPrhgAsdPzGXh5DymjcnEGzI/sgykhQ3gq8Cf/KU3tgOfxetmfdDMrgSKgUv9c8cDHQDOuSoz+wHwhn/s+865Kv/7bwP3m9kPgbfwJzUATwLvB7YCTf69REREpJ/W764D4JjCUf2qPyk/g311rbS0h0hLThrM0PrtreJqmttD3PyRY7lg/gSyUgea6sSGAb0K59zbwKJuDp3VTdlJwG0Rde8C7urmmtvxlvroWu6Aa/obq4iIiBxsXVktY7JSGJfTvz00O2eKllQ1MXNcbMy2XFtWS1LA+NCCQtJTYiOJHAzDlnY6524drnuJiIjI4a3fXcfcCaP63UU4KWKmaKwkbGtKa5k5NiuukjXQ1lQiIiIJqbUjxOZ99RzTz/FrAJNHZwKwK0ZmijrnWFtWy/yi/nXxxjIlbCIiIglo894GOsKu3+PXAPIykslKDVIcI7sdlNU0U9XYxrFFudEOZdApYRMREUlA63fXAv2fIQpgZgc2gY8Fa0u91zR/AElorFLCJiIikoDW7a4lOy14YBxaf00enREzXaJrympJTjKOHh8b4+kGkxI2ERGRBLSurI6543MGvCbZpNEZlFY1EwpHf536taW1zBqXTWowviYcgBI2ERGRhNMRCrNpb92Axq91mpSfQVsozN66lkGIrP/iecIBKGETERFJONsrGmlpD3NMYf/Hr3WanO/NFI32FlUlVc3UNrdzbGFuVOMYKkrYREREEsy6ss4JBwNvjepcPLe4KrozRdeU1QCohU1ERETiw/rddaQlB5g2JnPA1xo/Ko1gwNgV5Ra2taW1pCQFmBUjC/gONiVsIiIiCWZdWS1zxucQTBp4GhBMClCUlx71maJrSmuZMz6blGB8pjbx+apERESkW+GwY8PuugGtv9bVxPwMSqKYsIXDjnVltRwbp92hoIRNREQkoZRUN1Hf2sExgzB+rdPk0RlD0iXqnOOVrRU41/uSITsrG6lv7WB+nE44ACVsIiIiCWVdWR3AoCzp0Wlyfia1ze3UNrUP2jUBnttczifuWMHja/b0et5afxKFWthEREQkLqzbXUswYMwclzVo15zkzxTdNcgzRd/aVQ3AX14v7vW8NaW1pAYDzBw7eK8p1ihhExERSSDrd9cN+m4AndtbDfaeom/7e4O+sq2SXb1sML+2rJZ5EwZnEkWsit9XJiIiIgdxzrG+rHZQFsyN1JmwDeY4Nucca0preM+sAgIGD7xR0u15obD3mo6Nww3fIylhExERSRD76lqpbGwblAVzI2WmBhmTlTqoux2UVDVT09TOOfPG8d7ZY/nrqlLaQ+FDzttR0UBjW4hji3IH7d6xSAmbiIhIgujc4WCwW9jAnyk6iGPY3i6tAWBBUS6XLZ5EeX0rz27af8h5a/xu03jd4aCTEjYREZEEsW53LWZw9FGDn7BNys8Y1Ba2NSU1pAQDzD4qm/fOLmBsdmq33aJrSmtJT05iekH8TjgAJWwiIiIJY/3uOqaNySQzNTjo156Un8GeuhZaO0KDcr01pd5EguSkAMGkAJcuKuK5d/azp7b5oPPW+mPykgI2KPeNVUrYREREEoQ34WBoug4nj87AOSitbj78yYfREQqztqyWBRHj0j62aBJhBw++UXrQeet313JsHC+Y20kJm4iISAKoamxjd23LoG5JFWmyvxbbYHSLbi1voLk9dNC4tEmjMzhtxhgeXFlCKOwOnNfSHo778WughE1ERCQhrN/tTzgY5BminSYeWNpj4BMP1pR4sS6YmHtQ+WWLJ1JW08yLW8oBWFsa/zscdFLCJiIikgA6t6Qa7CU9OhVkpZKRkkRx1cC7RFeX1pCdGmTq6MyDys+eO468jOQDkw/WltWS1c158UgJm4jIECiubOK7j6ylpX1wBmCLDNS63bUU5aUzKiN5SK5vZt5M0UFY2mN1aQ3HFo0i0GUiQWowiY8sLOKZDfsor29lTak34aDrefFICZuIyCALhR3ffPBt7nutmNd3VEU7HBEANuyuG7Lu0E6T8jMGvNtBS3uITXvqD+kO7XTZ4ol0hB0PvFHMhj11cb/DQSclbCIig+zuV3ayyt+0ev3uuihHIwL1Le3sqGgckgVzI3ktbE2E/UkB/bFxTx0dYceCHsalzRibzYlT8vjd89tp6wjH/Q4HnZSwiYgMoh0Vjfz06U2cdfRYCnPTDwz0FommzfvqAZgzfmgTtsmjM2jtCFPe0Nrva6wuqQFgfi+J2MdOnERDa4d3nlrYRETkSITDjm89tJqUpAA/+vCxzJuQwwa1sEkM6OymnDJmaAfnT/IH/w+kW3RNaS1jslIZPyqtx3MuOHY82WlBstOCB5YTiXdK2EREBsk9r+7kjZ3V/PcH5zEuJ415E0axvaLxQEuASLTsqmzCDIry0of0PpMHYWmP1aU1HDdxFGY9TyRIT0niurNnceVpU3s9L54M/t4UIiIJaFdlIz956h2Wzi7gIwsLAQ4sULpxTx0nTsmPZniS4Eqqmhifk0ZqMGlI7zMhN52kgLFlf0O/6te1tLOtvJELjys87LmfOXVqv+4xUqmFTURkgLyu0DUEA8aPP3zsgb/4O7cAWl+mcWwSXbuqmpg0DF2HKcEAS2cV8MAbJdS1tB9x/XX+QriJsHPBkVLCJiIyQPet2MWKHVX8vw/MZfyod7ucxuWkMjozRTNFJep2VTYxOX94Fpf95tmzqG1u5w8v7Tziuqv9hG1Bgsz8PBJK2EREBqCkqomb/rmJM2YVcOmiooOOmRlzJ+SwTgmbRFFTWwcVDa3D0sIGXsvyOXPHccdL26ltOrJWtjWlNUzKzyAvM2WIohu5lLCJiAzAdx9ZR8AO7gqNNG/CKLbsq6e1QzseSHQUV3kzNiflD99sym+ePYv6lg7ueGn7EdVbXVKj7tAeKGETEemnupZ2XthSzmdPnUJhbvez7+ZNyKEj7Niyr3+DsEUGqrhy+BO2OeNzuODY8dz10g6qGtv6VKe8vpXdtS3qDu2BEjYRkX5atbMa5+Dk6aN7POfAxAMtoCtR0tnCNtzrlX39fTNpag9x+wt9a2VbU1oD0OOWVIlOCZuISD+t2FFFcpJx/MS8Hs+ZnJ9BVmpQEw8kaoqrmshJC5KbMbzjwmaNy+ZDCyZwzys7qejDzgerS2oIGEO+fdZIpYRNRKSfXt9RyfyiXNJTel7bKhAw5ozPZp2W9pAo2VU5PEt6dOdrZ82ktSPE757bdthzV5fWMnNsNhkpWiK2O0rYRET6obktxNqy2j4tiDtvwig27qknNIANsUX6q7hq+Jb06Gp6QRYXHV/IH1/bxf66lh7Pc86xprSGBRM14aAnA07YzCzJzN4ys8f95/lm9oyZbfG/dttXYGZX+OdsMbMrIsqnmtkKv/wBM0vxy83Mfm1mW81sjZktHGjsIiL99VZJNe0hx5KpfUnYcmhuD7Gjov/b9Yj0RyjsKK1uYuIwTjjo6utnzaQj7PhNL61sJVXNVDe197rhe6IbjBa2rwMbI55fDyxzzs0ElvnPD2Jm+cCNwBJgMXBjRGJ3M/ALv341cKVffj4w039cDfx2EGIXEemX13dUYQYnTOl5/FqneRM08UCiY09tM+0hF9UN0iePzuSShUX8eUUxe2qbuz1ndeeEAyVsPRpQwmZmRcAFwB0RxRcC9/jf3wNc1E3Vc4FnnHNVzrlq4BngPPMWMToTeKib+hcC9zrPa0CumY0fSPwiIv31+o4q5o7PISct+bDnzhyXRUpSgA2aeCDDLBprsHXnK2fOwOG45dmthLsZGrCmtIaUpACzj8qOQnQjw0BH9v0S+BYQ+Q6Pc87tAXDO7TGzsd3UKwRKIp6X+mWjgRrnXEeX8t7q7BngaxAROSJtHWHeLK7m44sn9en85KQAs47KYp1a2GSYRWMNtu5MzM/go4sm8qcVxfx5RTEpSQFSgwFSkwOkBpOobmpj7oQcUoIaWt+TfidsZvYBYL9zbpWZLT3S6t2UuV7Ke6vTXWxX43WbMmlS3z5QRUT6at3uWlrawyzuw4SDTvPGj+LpDXtxznW7I4LIUNhV1UQwYEzoYWHn4fSd989hekEWdS3ttHaEaW0P09oR8r7vCPOhBROiHWJMG0gL26nAh8zs/UAakGNm9wH7zGy837o2HtjfTd1SYGnE8yLgOaACr6sz6LeyFQG7I+pM7FJnN91wzt0O3A6waNEiTcsSkUH1+o4qAE7sw4SDTscU5vDAyhJ217b0uCuCyGArrmqiKC+dpED0/0jITA3yudOmRjuMEavfbY/OuRucc0XOuSnAZcCzzrlPAY8BnbM+rwAeBTCzQjNb5pc/DZxjZnn+ZINzgKedcw5YDlzStb5/3cv92aInAbWdXa8iIsPp9R1VTC/IZExWap/rzO2ceKD12GQYFVc2MWl0dJb0kME1FJ3FNwFnm9kW4Gz/OcB4oAPAOVcF/AB4w3983y8D+DZwrZltxRvTdqdf/iSwHdgK/B748hDELiLSq1DY8cbOKhZP7Xk7qu7MGZ+NGazTxAMZRsVVTUzKV4tuPBiU5YSdc8/hdWninKsEzurmtJOA2yLq3AXc1c21tuMt9dG13AHXDEa8IiL9tWlvHfUtHX1afy1SRkqQaWMy2aCJBzJMapvaqW1uj9qiuTK4hm3/B+fcrcN1LxGRodKf8WudjikcdaC+yFDrXNIjmovmyuDR/FkRkSPwxs4qCnPT+zVxYN6EHPbUtlDV2DYEkYkcbFeVt7NGNBfNlcGjhE1EpI+cc7y+o+qIu0M7accDGU67YmQNNhkcSthERPpoe0UjFQ1tLO53wpYDwLoyTTyQoVdS1cSYrBQyU4dt9JMMISVsIiJ91Dn+rL8JW25GCoW56Wphk2Gxq7JJrWtxRAmbiEgfvb6jijFZqUwd0/9Zd/Mm5GhPURkW3pIeStjihRI2EZE+en1HFYun5g1oa6l5E0axo7KRhtaOw58s0k9tHWH21DZr0dw4ooRNRMTX2hHCW/LxUKXVTZTVNB/R/qHdmTchB+dg4x61ssnQKa1uIuxgslrY4oYSNhERoLktxCk/fpbLbn+NEn/9qkhv7Owcv3ZkOxx0Na/Qm3igLapkKHWuwTZJS3rEDSVsIiLAmtIaKhvbeGNnFef/6kUeXFlyUGvb6zuqyEkLMvuo7AHd56icNEZnprBe49hkCHUmbGphix9K2EREgFXF1QA8es1pzJuQw7ceWsNV966ivL4VgBU7qjhxSj5Jgf6PXwMwM06YnMdzm8tpD4UHHLdId4orm0hLDlCQnRrtUGSQKGETEQHe3FXNtIJMji0axV+uOonvXjCHF7aUc+4vX+DPK4rZXt7Yr+2ouvOxEydSXt/KMxv2Dcr1RLra5c8QHcgEGYktSthEJOE551i1q5oTJuUBEAgYnz99Gk989TQKc9P5zt/XAv1ff62rpbPHUpibzn2v7RqU64l0VaIlPeKOEjYRSXg7KhqpbmrnhMl5B5XPHJfNw18+hW+8byZLZxdwbOGoQblfUsD4xJJJvLKtkm3lDYNyTZFOzjl/DTYt6RFPlLCJSMJbtcsbv9Y1YQNITgrwjffN4u7PLiY5afA+Mi9dVEQwYPx5RfGgXVMEoLyhlaa2kDZ9jzNK2EQk4b1ZXE1OWpDpBVnDds+x2Wmce8xRPLSqlJb20LDdV+Jf57I06hKNL0rYRCThrdpVzcLJeQQGOAP0SH1qyWRqm9t5fM2eYb2vxLddlVqDLR4pYRORhFbb3M7mfQ0HJhwMp5Om5TO9IFOTD2RQFVc1YQZFeenRDkUGkRI2EUlobxX3PH5tqJkZn1wymbdLalinnQ9kkBRXNjE+J43UYFK0Q5FBpIRNRBLam7uqCRgsmJgblft/ZGERackB/qTJBzJIiqua1B0ah5SwiUhCW1VczZzxOWSmBqNy/1EZyXxw/gQefbuM+pb2qMQg8WWX1mCLS0rYRCRhdYTCvF1cE5Xu0EifOmkyTW0hHnmrLKpxyMjX1NZBeX0rk0drDbZ4o4RNRBLWO/vqaWwLRT1hm180imMKc7jvteKDNpwXOVIlVc0ATFQLW9xRwiYiCetNf8HchVGYIRrJzPjUksm8s6/+wCK+Iv2xq7IRgMlK2OKOEjYRSVirdlUzNjs1JpY/+NBxE8hODWqJDxmQYi2aG7eUsIlIwlpVXM0Jk/MwG94Fc7uTkRLkwwsLeXLtXiobWqMdjoxQxVVNZKcFyc1IjnYoMsiUsIlIQtpf10JJVXPUx69F+sSSybSFwtr5QPptV2UTk0dnxMQfITK4lLCJSEJ6018wd2EMJWyzj8pm9rhsnlDCJv0QDjvWltUya1x2tEORIaCETUQS0qpd1aQEA8ybkBPtUA5ywfzxvLGrir21LdEORUaYd/bVU9XYxinTx0Q7FBkCSthEJCGt2lXNgqJRMbd9z/uPHY9z8M91amWTI/PKtkoATp4+OsqRyFBQwiYiCaelPcS6srqY6g7tNGNsFkcfpW5ROXKvbqtgyugMCnOjP+tZBp8SNhFJOOt319IWCnNClNdf68kH5o9n5a5q9tQ2RzsUGSE6QmFWbK/iZHWHxi0lbCKScDoXp43FFjbwukUBnly7N8qRyEixbncd9a0dnKLu0LilhE1EEs6qXdVMGZ3BmKzUaIfSrWkFWcwdn8MTa3ZHOxQZIV7ZVgHASdOUsMUrJWwiEhN2VDQe2FZnKDnnWLWrJmZb1zpdMH88bxbXUFajblE5vFe3VTJ7XDYF2bH5R4gMnBI2EYm6/fUtXHjrS7znp8/x6TtX8NS6vXSEwkNyr5KqZioaWmNqwdzuXOB3i/5zrSYfSO9aO0K8sbNKs0PjnBI2EYm6Hz6+kZb2MF9aOp1t+xv44n2rOO3m5fzy35sHfT2yVcVVADGfsE0Zk8kxhTna9UAO6+3iGlrawxq/FueUsIkIpdVNrCurjcq9X9hczmOrd/OlpdP59nlH88K33svvL1/ErKOy+eW/t3Dqzc9yzZ/fpL6lfVDut3xTOdlpQWaOjf3V4D8wfwJvl9RQ4m/oLdKdV7ZVEjBYovFrcU0Jm0iCC4cdn7v7DT78m1d4u6RmWO/d0h7i/z26jqljMvnS0ukABJMCnD13HPd+bjHP/+dSrjxtKk+s2cNty7cN+H7byht4fM1uPr54EkmB2N9r8UC3qBbRlV68uq2SYwpHMSpdG77HMyVsIgnuXxv2sXlfA0kB44t/XEV5feuw3fu25VvZVdnE/1x0DGnJh+44MHl0Jt95/xwuPr6Qu17eMeAB+Lc+u5XUYBJXnzFtQNcZLhPzM1hQNEqL6EqPmttCvFVSrfFrCUAJm0gCc85x2/KtTMrP4MEvnExNcxvX/OlN2jqGZsB/pK376/nd89u4+PhCTpnR+2Kf150zC4Cf/eudft9vW3kDj75dxuUnT47Z5Ty6c8H88awuraW4Ut2icqiVu6poDzntH5oA+p2wmdlEM1tuZhvNbL2Zfd0vzzezZ8xsi/+125G9ZnaFf84WM7sionyqma3wyx8wsxS/3Mzs12a21czWmNnC/sYuIp4XtlSwtqyWLy2dzrFFo7j5I/N5fWcVP3xiw5De1znHf/19HenJSXzn/XMOe35RXgafPWUKf3+rjPW7+zfW7pZlW0gNJnHVCGld69S5iO4Tmi0q3XhlWyXBgLEoxifRyMANpIWtA7jOOTcHOAm4xszmAtcDy5xzM4Fl/vODmFk+cCOwBFgM3BiR2N0M/MKvXw1c6ZefD8z0H1cDvx1A7CIC3PbsVo7KSePDCwsBuPC4Qq46fSr3vrqLB1eW9Fq3uS3E/rr+zeB8aFUpK3ZUcf35c/q8btSX3zuDUenJ3PTPTUd8v637G3hs9W4uP2Vkta6Bl6weNzGXJ9ZqEV051CvbKjluYi6ZqcFohyJDrN8Jm3Nuj3PuTf/7emAjUAhcCNzjn3YPcFE31c8FnnHOVTnnqoFngPPMzIAzgYe6qX8hcK/zvAbkmtn4/sYvkuhe31HF6zur+MJ7ppEafHf82LfPO5pTZ4zmu39f1+0khIqGVn7+r3c4+aZlnP6T5Ty/ufyI7lvV2MaPntzICZPzuOzEiX2uNyo9ma+eOZMXt1Qc8T1veXYLaclJXH36yGpd6/SB+eNZV1bHzoqhX1hYRo66lnbWltZoOY8EMShj2MxsCnA8sAIY55zbA15SB4ztpkohEPnne6lfNhqocc51dCnvrY6I9MOty7cyOjOFy06cdFB5MCnArR9fyNicVL74x1Xsr/da0baXN3DDw2s55aZnuWX5Vk6cks/0giyuumclyzbu6/N9b/rnRupbOvifi48hcIQzNT910iQm5qfz4yc3Egq7PtXZur/ea107eQqjR1jrWqfz1S0q3Xh9exVhhzZ8TxADTtjMLAv4G/AN51xdX6t1U+Z6Ke+tTncxXW1mK81sZXn5kf0lLpII1pTW8MLmcq48fSrpKYfOzszLTOH2Ty+iprmNL/xxFVfdu5Kzfv48f3uzlI8sLOLf176H31++iD9ftYSjx2fzxftW8fT63jcqD4cd979ezIMrS7ny9KkcfVTOEcedGkziW+cezaa99Tz8Zmmf6vx62VbSk5O46vSpR3y/WFGYm87CSbk8/GbpkO0AISPPK9sqSQ0GOH5SbrRDkWEwoITNzJLxkrU/Oece9ov3dXZV+l/3d1O1FIjsCykCdgMVeF2dwS7lvdU5hHPudufcIufcooKCgiN/YSJx7rblW8lJC/Lpkyb3eM7cCTn85JIFvFVcwxs7q/jqe2fwyvVn8uMPH8v0giwAcjNSuO/zSzimcBTX/OnNHpef2LC7jo/+36tc//BaFk3O4+tnzex37B+YP54FRaP42b8209wW6vXcLfvq+ceakd261unqM6azrbyRO1/aEe1QJEa8sq2CRVPyul0SR+LPQGaJGnAnsNE59/OIQ48BnbM+rwAe9c8vNLNlfvnTwDlmludPNjgHeNo554DlwCVd6/vXvdyfLXoSUNvZ9Soifbd5Xz1Pr9/HZ06ZQnZa7wttfmjBBJ76xum8cv2ZXHvO7G4H7OekJXPv5xZz/KRcvvqXN3n07bIDx+pa2vneY+v5wC0vsr2ikZ98ZD4PfuFkMlL6P0DazPjO++ewt66Fu17uPXn51bItpCePnHXXenPuvHGcPXccv/j3Zi3xIVQ2tLJpb72W80ggA2lhOxX4NHCmmb3tP94P3AScbWZbgLP95wDj8WaW4pyrAn4AvOE/vu+XAXwbuNbMtuKNabvTL38S2A5sBX4PfHkAsYskrN8s30pGShKfPbVvXYRHH5Vz2AQrOy2Zuz+7mMVT8/nmA2/z0KpSHn6zlDP/93nueXUnn1gyiWevew8fPXHiEY9b686SaaN535xx/Pa5bVQ2dL/Q7+Z99Tyxdg9XnDKF/MyUAd8z2syM7184j2AgwH89shbv71tJVK9t9/7L1IK5iaPff+Y6516i+3FlAGd1U3YScFtE/buAu7q57na8pT66ljvgmn4FKyIA7Kps5LHVu7nytKnkDXISk5ka5A+fWcxV967kP/66GoAFE3O56zOLmF+UO6j3Arj+/KM595cvcNW9K1k4KY8x2akUZKUyJjuVMVkp3LJsKxnJSVw1QmeGdmf8qHS+dd5s/vvR9TzydhkXH18U7ZAkSl7ZVkFWapD5haOiHYoMk2FbuMU5d+tw3UtEuve757cRTAoMWRKTnpLEHVcs4uanNjF7XDYfXTQ4LWrdmTE2i2+fN5s/vraL+1bsoqX90MH4X146PS5a1yJ9aslkHnmrjB88vpH3zBobd69P+ubVbZUsnppPMEkbFiUKrbQnkiD21bXw0KpSPnbiRMbmpA3ZfdKSk7jxg/OG7PqRrj5jOlefMR3nHI1tISrqWylvaKWivpX6lg4+sCD+lmoMBIwff3g+F/z6RX74xAZ+/tHjoh2SDLM9tc1sr2jkE0smHf5kiRtK2EQSxLOb9tMeclxx8pRohzLozIys1CBZqUGmjMmMdjhDbvZR2XzxPdO5dflWLj6+kNNnajZ8Inlmg7fuocavJRa1pYokiJe2VnBUThozxmZFOxQZBF85cwZTx2TyX39fd9jlTSR+hMKOu17awYKiUcwdf+RrGcrIpYRNJAGEw45XtlZw6owxeCvyyEiXlpzEjy4+luKqJn65bHO0w5Fh8syGfeysbOLqM6br33KCUcImkgA27Kmjuqmd02dqzaZ4cvL00Xxs0UTueHEH63fXRjscGQa3v7CNifnpnDtvXLRDkWGmhE0kAby0tQKAU2ZozEu8ueH9R5ObnsyNj67X2mxxbtWuKt4sruHKU6dqdmgC0k9cJAG8vLWC2eOyGZs9dLNDJTpyM1L41nmzWbmrmsdWd7tbn8SJ/3t+O6PSk/noiRMPf7LEHSVsInGupT3E6zuqOHWGukPj1aUnTGR+0Sh+9ORGGls7oh2ODIHt5Q08s3Efnz5p8oC2dpORSwmbSJx7c1c1rR1hTpup7tB4FQgYN35wHvvqWrl1+dZohyND4I6XdpAcCHDFKVOiHYpEiRI2kTj34tYKggFj8VQlbPHshMl5fHhhIXe+uIMdFY3RDkcGUUVDK39bVcqHFxZSkJ0a7XAkSpSwicS5l7dWcPykXLJS1Y0S764/72hSggF++PiGaIcig+iPr+6itSPM5+NoX1w5ckrYROJYTVMba8tqOW2GVsJPBGNz0vjaWTNYtmk/yzftj3Y4Mgia20Lc++pO3jdnrBa9TnBK2ETi2KvbKnEOjV9LIJ85ZSrTxmTy/cc30NYRjnY4MkAPvVlKdVM7V58xPdqhSJQpYROJYy9trSArNcj8otxohyLDJCUY4L8/OJcdFY3c9fKOaIcjAxAKO+54cTsLJuZy4pS8aIcjUaaETSSOvby1gpOm5ZOsRTYTytLZY3nfnLHcsmwL++paoh2O9NMzG/ayq7KJL5wxTdtQiRI2kXhVUtXEzsomrb+WoL57wVzaQ44fPblROyCMQKGw47fPbWNSfgbnzjsq2uFIDFDCJhKnXva3ozpNCVtCmjImky++ZxqPvr2b6x5cTUt7KNohyRH43fPbWF1ayzfeN5OkgFrXBDTPXyROvbS1gnE5qZpZlsC+8b5ZBJMC/PyZzWzZ38DvPn0Chbnp0Q7rgPZQWN313VhdUsMvntnMBceO5+LjC6MdjsQI/UsRiUPhsOOVbZWcOmOMxr4ksEDA+NpZM7nj8kXsrGjkQ7e8xGvbK6MdFo2tHfz4nxuZ999Pc88rO6MdTkxpbO3g6/e/xdjsVH508bH69ysHKGETiUMb99ZR1dim7lAB4H1zx/HIV05lVEYyn7pjBfe8sjMq49qcc/xj9W7O+tnz/N/z2xmbk8r/PLGR9btrhz2WWPW9x9azq6qJn3/sOEZlJEc7HIkhSthE4tBLW7zxa5pwIJ2mF2TxyDWnsnR2ATc+tp7/fGjNsI5r27Kvnk/esYKv/uUtRmel8LcvncJjXzmN3IxkvvaXt2hu0xi7J9bs4a+rSvny0umcNE1rJ8rBlLCJxKGXtlYwc2wW43LSoh2KxJCctGRu//Qivn7WTB5aVcrZv3ieR98uIxweuta2htYOfvTkRs7/1YusK6vlBxcdw2NfOY0TJueRn5nCLz52HNsrGvnBE4m9nVZZTTM3PLyGBRNz+cb7ZkU7HIlBSthE4kxLe4g3dlapdU26FQgY3zx7FvdduYSs1GS+fv/bXHjbywdmFQ+W5rYQv39hO2f8ZDm3v7CdjywsYvl/LOXTJ00+aNbjqTPGcPUZ0/jzimKeWrd3UGMYKUJhx7UPvE0o7PjVx47TRAzplmaJisSZN4uraWkPa/ya9Oq0mWN44qun8cjbZfzsX5v55B0rOGNWAdefdzRzJ+T0+7qtHSH+sqKY257bRnl9K6fPHMN158zmuIm5Pda57uzZvLK1kusfXsNxE3M5alRitQz/7vltrNhRxU8vmc+UMZnRDkdilNJ4kTjz8tYKkgLGkmn50Q5FYlwgYHx4YRHLrnsP371gDmtKa7jglhe59oG3qWxoPaJrtYfC/HlFMUt/+hzf+8cGpo7J5IGrT+KPVy7pNVkDbzutX112HG0dYb7ptzQlgnDYm4Txi2c2c8H88VxyQlG0Q5IYphY2kThS39LOP1bv4biJuWSnaYaZ9E1achKfP30aly6ayO+e38adL+7ghS3l/OjiYznnMKvsO+f4x5o9/PTpTZRUNXP8pFx+eskCTp0x+oiWpJhWkMX3PjSPbz20httf2M6XlsbvZuehsOPxNbu5bflWNu9rYNa4LH50kZbwkN4pYROJE+Gw47oHV1NW08xNHzk22uHICDQqPZlvn3c0Fx43gWsfWM3Vf1zFRxYWceOH5pLTzR8Aq3ZV88MnNvBWcQ1zxufwh88cw9LZBf1OPC49oYjnN5fzs3+9wynTR7PgMC1zI017KMwjb5Xxm+e2saOikZljs/jVZcdxwbHjCWrcmhyGxfsec4sWLXIrV66MdhgiQ+6WZVv42TOb+e4Fc/j86dOiHY6McG0dYW55dgu/eW4b47JT+emlCw5MZCmpauKmpzbxxJo9jM1O5T/Onc1HFhYNyhZKtU3tvP/XL2IGf7nqJCbmZwz4mtHWHgrz15Wl/Oa5rZRWNzN3fA5fPXMG5847ioC2nZIuzGyVc27RIeVK2ERGvmc37ePKe1Zy4YIJ/OJjx6lrRQbN2yU1XPvg22wvb+TykyeTnpzEH17eSVLAuPqMaXzhPdPISBnczprVJTVcftfrpAYD3HvlYo4+qv+TIKLJOcfT6/fyk6feYXtFI8dNzOVrZ83gvbPH6t+o9EgJm0ic2lHRyIdufYmJeRn87UunkJ6SFO2QJM60tIe4+alN/OHlnZjBJQuLuO6c2UM6m3Pzvnouv/N1mto6uPMzJ3LilJE1iWblzip+9ORG3iyuYcbYLK4/72jOmqNETQ5PCZtIHGpo7eDi216moqGVx75yWlx0H0nserukhrTkwLC1eJVWN3H5na9TVtPMbZ9YyPvmjhuW+w7E1v0N/OSpTfxrwz7GZqdy7dmzuOSEIo1Rkz7rKWHTpAOREco5x3/+dTXbyhv445VLlKzJkDvc8hyDrSgvg79+8WQ+e/cbfOG+Vdz8kfkxufRFeX0ry9/Zz7837GPZpv2kJydx3dmzuPL0qYPeXSyJS79JIiPUb57bxj/X7eW/3j9HuxpI3BqdlcqfrzqJL/xxJf/x19VUNbZy9RnRXfLDOcemvfUs27iPf2/cz+rSGpyD8aPS+OwpU/ji0umMyUqNaowSf9QlKhLjwmHHnroWtu1vYFu5/9jfyGs7Kvng/An86jJNMpD419oR4toHVvPE2j3MLxrFkqn5LJ46mhOn5JGbkXLI+eGwo6S6iY176tm0t462jjCj0pMPPHL8r5mpQRpbO6hv6aC+pZ2GiO/rWztoaOk4UNbQ0kF9awcVDa2U13sLCy8oGsVZc8Zx1pyxzB2fo3+LMmAawyYygjjneG5zObc/v53VpTU0tYUOHMtOCzJjbBYLinL59nlHa5KBJIxQ2HHHi9tZtnE/b5fU0BYKAzB7XDaLp+YzrSCTrfsb2LS3nk176mj0/92YQcDsiHdQSAkGyEkLkpUaJKvza6qX6C2emsd7jx7L2OzE2kZLhp4SNpERIBx2/HvjPm5dvpU1pbVMGJXGucccxYyxWUwv8B5jslL0V7wkvJb2EKtLanhjZxUrdlSxalc1TW0hslODzBmfw9Hjs5kzPoc543OYPS6btOQAjW0hapvbqWtup9Z/NLV1kJESJDstSHZqsvc1zUvQUoP6Y0iGnxI2kRgWDjv+uW4vtzy7hU1765mUn8E1753OxccXkRLU7DKRw+kIhalsbGNsdqr+oJERTbNERWKMc45t5Y0s37SfB1aWsHV/A9MKMvn5RxfwoQUTtAyAyBEIJgUYl6PuSYlfSthEhlFrR4gV26t4dtN+nt20n+KqJgCOKczhlo8fz/uPHT8o2/uIiEh8UcImMgRa2kOUVjdTUtVEsf/YXt7Aih1VNLWFSA0GOHXGGK4+YxrvPXoshbnp0Q5ZRERi2IhL2MzsPOBXQBJwh3PupiiHFDUdoTCtHZ0PbzZUwAwzMIyAPzMKoD0cpiPk6Ai5d78Ph8lICZKfkUJ2WjCuNyF2ztEecrSFwrR1RDxCIdo6HO2hMB3hMG0d3vvSHgrTHnI0tHRQ2dhKZWMblQ1tVDW2UdnYRk1TG2HnvPebiPfdjMbWDvbWtRA5PDQtOcCk/AwuPr6Qs+aM5eRpYzS7U0RE+mxEJWxmlgTcBpwNlAJvmNljzrkNwxVDRyhMc3uIlvYwLe0hmttDNLeF6PCToFDY0RF2hJwjFPK+B0fYgXMQdo6w/z95a3uY+tYOGv1Hg/9obA3R2hGitd1LxFr8r53JWUu79/2RTlHvTcAgLyOF3Ixk8jNTyElLJilgJAWMQMAImJFkHPje4EBi6CUq3nXaQ+5AMtQeCh9IkByQkhQgOclITgqQHAwc9Dwl6D1SI74PmNERdrR3hGkP+0mVn0h5CdWh37d1hGlq66CpzfvZNLV5P5+mtg4G8nalBAOMzkwhPzOF0VmpTBmdQZKZ//MEh/ezdc6RnhxkUn4Gk0anMyk/g4n5GRRkaSC0iIj034hK2IDFwFbn3HYAM7sfuBDoMWGrbmzjzyuKCTlHOOwlVGEXkVgd+BqmI+y1QDW2dlDX0k5ds/e1cxp4Q2sH7aGhmVWbGgyQlRokMzVIRkoS6SlJpAYD5GWmkBoMkBr0nqcmB0gLJpGa/G5ZWnISKcEABoTdu4mDw5t96PAG5CYHzPuaZAQDAZICRlNbB1WNbdQ0tVPV5LUcVTW2sae25cD7FHIO57w1kELhd6/tHDic/9XTmYSlBANeYuYnYABNbR3dJljtflLX6n/tSWfcwSQjNRggGAiQHPQTvqTO+xkZKUFGZ6WSkZLkvZfJ3nuaGnw3GUzxE8aUYIDUiFiDSXbgWsEkIys1SH5mClmpQSVcIiISNSMtYSsESiKelwJLup5kZlcDVwOkHDWD7/x97WEvHPRbk5IC3n/SOenJ5KR5/1lPGZ1JTrq3YKKXACSRlpJEWjBAekoSaUEvYYq8RlLASy4CgXe7y95tmfKeRyZpyZoRCBzcdRkKuwMtcMGAKWESEZGENdIStu7+xz6kycs5dztwO8CC4xe6p284i0AAkuzdLr4kM4J+i03AH3sk0WdmpARNa4+JiIhEGGkJWykwMeJ5EbC7twrJSQGOGqW1eURERGTkGmnNGG8AM81sqpmlAJcBj0U5JhEREZEhNaJa2JxzHWb2FeBpvGU97nLOrY9yWCIiIiJDakQlbADOuSeBJ6Mdh4iIiMhwGWldoiIiIiIJRwmbiIiISIxTwiYiIiIS45SwiYiIiMQ4JWwiIiIiMU4Jm4iIiEiMU8ImIiIiEuOUsImIiIjEOCVsIiIiIjHOnHPRjmFImVk98M4w3GoMUDEM9xkF1A7DfWIhhuF6T3szXK91OH+u0X5fh/O16nd15N5H7+nQiPb7qn//hzfbOZfdtXDEbU3VD+845xYN9U3MbOUw3ed259zVQ32fWIhhuN7Tw8QwXK912H6u0X5fh/m16nd15N5H7+nQ3Ev//gf/PoP6nprZyu7K1SU68vwj2gEQGzEMl+F6rXpPR/69ok2/q4NP7+ng07//flLCNsI456L+CxgLMQyX4Xqtek9H/r2iTb+rg0/v6eDTv//+S4SE7fY4u08i0Xs6NPS+Dj69p4NP7+nQ0Ps6+Ab7Pe32enE/6UBERERkpEuEFjYRERGREU0J2xEws4vNzJnZ0dGOJR6ZWcNhjj9nZlGdNTYSmFmRmT1qZlvMbJuZ/crMUno5/xtmljGcMY5Uh/sdlSOnz9Who8/UwRMLn6tK2I7Mx4GXgMuOpJKZJQ1NOCIHMzMDHgYecc7NBGYBWcD/9FLtG4ASNokWfa5KTIuVz1UlbH1kZlnAqcCV+B8sZrbUzF4ws7+b2QYz+52ZBfxjDWb2fTNbAZwcvchHFv89fTzi+a1m9pkohjTSnAm0OOf+AOCcCwHfBD5nZplm9r9mttbM1pjZV83sa8AEYLmZLY9i3COGmWWZ2TIze9N/Ly/0y6eY2UYz+72ZrTezf5lZerTjjWX6XB16+kwdFDHxuZoIC+cOlouAp5xzm82syswW+uWLgbnALuAp4MPAQ0AmsM4599/RCFYS1jxgVWSBc67OzIqBzwNTgeOdcx1mlu+cqzKza4H3Oueivar8SNECXOy/r2OA18zsMf/YTODjzrmrzOxB4CPAfdEKdAS4CH2uSuyLic9VtbD13ceB+/3v7/efA7zunNvuZ9x/AU7zy0PA34Y3RBEM6G7qtwFnAL9zznUAOOeqhjOwOGLAj8xsDfBvoBAY5x/b4Zx72/9+FTBl2KMbWfS5KiNBTHyuqoWtD8xsNF6T6DFm5oAkvB/ekxz6Q+x83uJ/2MiR6eDgPyTSohXICLUer1XnADPLASYC2+n+Q0eOzCeBAuAE51y7me3k3d/T1ojzQoC6RHugz9Vho8/UgYuJz1W1sPXNJcC9zrnJzrkpzrmJwA68v/oWm9lUf4zFx/AGz0r/7QLmmlmqmY0Czop2QCPMMiDDzC6HAwOzfwbcDfwL+KKZBf1j+X6deuCQjYalR6OA/X6y9l5gcrQDGqH0uTo89Jk6cDHxuaqErW8+Dvy9S9nfgE8ArwI3AevwPmy6nid94P+ytzrnSoAHgTXAn4C3ohrYCOO8lbAvBi41sy3AZrwxV98B7gCKgTVmthrv9xe8VbX/qUkHvev8HcX7vVxk3gbNnwQ2RTWwkUufq0NIn6mDJ1Y+V7XTwQCY2VLgP5xzH4hyKCOemS0Afu+cWxztWES6o9/R4aHP1cGh39f4oxY2iToz+yLewOLvRjsWke7od1RGEv2+xie1sImIiIjEOLWwiYiIiMQ4JWwiIl2Y2UQzW+7vXLDezL7ul+eb2TPm7Sf4jJnl+eVnm9kqf7XzVWZ2ZsS1/sfMSkz7kIrIAKhLVESkCzMbD4x3zr1pZtl4i+BeBHwGqHLO3WRm1wN5zrlvm9nxwD7n3G4zOwZ42jlX6F/rJLylFbY457Ki8XpEZORTwiYichhm9ihwq/9Y6pzb4yd1zznnZnc514AKYIJzrjWivEEJm4j0l7pERUR6YWZTgOOBFcA459weAP/r2G6qfAR4KzJZExEZKG1NJSLSAzPLwlvM9Rv+Zs+HO38ecDNwzjCEJyIJRC1sIiLdMLNkvGTtT865h/3ifX5XaOc4t/0R5xfhrch/uXNu23DHKyLxTQmbiEgX/ji0O4GNzrmfRxx6DLjC//4K4FH//FzgCeAG59zLwxiqiCQITToQEenCzE4DXgTWAmG/+Dt449geBCbh7R94qXOuysy+C9wAbIm4zDnOuf1m9hO8/QUnALuBO5xz3xuWFyIicUMJm4iIiEiMU5eoiIiISIxTwiYiIiIS45SwiYiIiMQ4JWwiIiIiMU4Jm4iIiEiMU8ImItIDM/uemf1HL8cvMrO5wxmTiCQmJWwiIv13EaCETUSGnNZhExGJYGb/BVwOlADlwCqgFrgaSAG2Ap8GjgMe94/V4m36DnAbUAA0AVc55zYNY/giEqeUsImI+MzsBOBuYAkQBN4Efgf8wTlX6Z/zQ2Cfc+4WM7sbeNw595B/bBnwRefcFjNbAvzYOXfm8L8SEYk3wWgHICISQ04H/u6cawIws8f88mP8RC0XyAKe7lrRzLKAU4C/eluRApA61AGLSGJQwiYicrDuuh3uBi5yzq02s88AS7s5JwDUOOeOG7LIRCRhadKBiMi7XgAuNrN0M8sGPuiXZwN7zCwZ+GTE+fX+MZxzdcAOM7sUwDwLhi90EYlnGsMmIhIhYtLBLqAU2AA0At/yy9YC2c65z5jZqcDvgVbgEiAM/BYYDyQD9zvnvj/sL0JE4o4SNhEREZEYpy5RERERkRinhE1EREQkxilhExEREYlxSthEREREYpwSNhEREZEYp4RNREREJMYpYRMRERGJcUrYRERERGLc/wfFSO1JInHgjQAAAABJRU5ErkJggg==\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"metadata": {
"trusted": false
},
"id": "92807341",
"cell_type": "code",
"source": "young[young.index<'2020-10-01'].cases.sum()",
"execution_count": 15,
"outputs": [
{
"data": {
"text/plain": "15619"
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
]
},
{
"metadata": {
"trusted": false
},
"id": "4940ef2c",
"cell_type": "code",
"source": "young[young.index>'2020-10-01'].cases.sum()",
"execution_count": 16,
"outputs": [
{
"data": {
"text/plain": "1152837"
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
}
]
},
{
"metadata": {
"trusted": false
},
"id": "5637c9eb",
"cell_type": "code",
"source": "",
"execution_count": null,
"outputs": []
}
],
"metadata": {
"_draft": {
"nbviewer_url": "https://gist.github.com/0a7273c3b76f3ceceda70eb0efaf4f43"
},
"gist": {
"id": "0a7273c3b76f3ceceda70eb0efaf4f43",
"data": {
"description": "nbs/uk-young-covid.ipynb",
"public": true
}
},
"kernelspec": {
"name": "python3",
"display_name": "Python 3 (ipykernel)",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.9.5",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": false,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment