{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Analysis of 5K race results from races run by both Jesse and Hugh\n", "This notebook analyzes 5K races run by Jesse Bloom and Hugh Haddox.\n", "\n", "### import modules" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%matplotlib inline\n", "import scipy.stats\n", "import pandas\n", "import matplotlib\n", "import matplotlib.pyplot as plt\n", "import seaborn\n", "seaborn.set_context('notebook', font_scale=2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Read in race results\n", "First, we read in a table of race results manually compiled from the internet." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "year race Hugh Jesse\n", "2013 Shore Run 22.22 19.42\n", "2015 Dawg Dash 20.17 19.27\n", "2017 Shore Run 20.37 19.93\n" ] } ], "source": [ "results = pandas.read_csv('race_results.csv')\n", "print(results.to_string(index=False))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot race results over time\n", "Next, we look at how the race results have changed over time for each runner, drawing a dashed line at the critical threshold of 20 minutes." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhIAAAFrCAYAAACE+GArAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XdclWX/wPHPgcOUoQxBQURzoOLAGaKl5Shz5cqVaI+j\nHNkwy6f015P11O/3lGWYIy1XWjhQ05ZpinugKBJOFAeyBNnrrN8f5DEehgc4cBjf9+vFi9d939d9\n39/DCw7fc93X9b0UOp1OhxBCCCFEOZiZOgAhhBBC1FySSAghhBCi3CSREEIIIUS5SSIhhBBCiHKT\nREIIIYQQ5SaJhBBCCCHKTWnqAGqipKQMU4cghBBCVClXV/ti90uPhBBCCCHKTRIJIYQQQpSbJBJC\nCCGEKDdJJIQQQghRbpJICCGEEKLcJJEQQgghRLlJIiGEEEKIcpNEQgghhBDlJomEEEIIIcpNEglR\naVRqLVqdztRhCCGEqERSIlsYVVauit9P3+ZwRBz3M/JQmivo1NKVZ7p70byxg6nDE0IIYWQKnU4+\nMpaVrLVRvNTMPP5301kS7ucUOWamUDB1cBseb+dugsiEEEJUlKy1ISrdtz9fLDaJANDqdHzz00Xu\npRZ/XAghRM0kiYQwirjkLCKvp5TaRqPVceBcbBVFJIQQoipIIiGM4srtVIPaXb5lWDshhBA1gyQS\nwii0WsOG2mg0MiRHCCFqE0kkhFE0dTdsRoZ3o+IH6wghhKiZJJEQRtGskT1N3R6dJPTp1LgKohFC\nCFFVJJEQRqFQKHjpuTbYWpVemiTyRukDMoUQQtQskkgIo2nS0I53J3WhSytXzBQK/X5L5cNfs5DQ\n60RE3zNFeEIIISqBFKQqBylI9WgZ2fkkp+dibanExtKcxRvCSEnPA8DGSsnCwK64O9maOEohhBCG\nKqkglSQS5SCJRNndjM/g39+dQaXWAtDI2Zb3JnXF5hGPQoQQYvbs6Wg0Gnr3fpING74FwNnZhZs3\nY1i7dhMtW7bWtz106CD//Oc8/vnP/2HQoCEA9OrVlUmTXsLdvRFbtmwmNvYOzs4uDB48jEmTXsLM\nrKDX9KOP3ics7BSffRbEV18tJSLiHGZmCrp168Hs2a/j7t5Ifx+VSsWmTev57befiY+Pw97egZ49\nezF16iu4uLjo240aNYT27Tvi5OTErl0hWFtb89FHn9KxY6eq+NEZVUmJhLyLiyrR1N2eKc/68PXu\nKADikrNZvTuK2SPbF3oMIoQQxbl27Qrx8XH84x8vk5KSTGpqKjdvxhh8/t69v5CTk82IEWNwdnZm\nz54fWbNmJfXr12f48FH6dllZWcyZM50ePXoye/ZrXL16mV27QkhMTOTrr9cBoNVqefvtNzhz5hTP\nPjuYF14Yz927sezYsY1Tp06wevV6nJ0fJhNHjhyiYcOGzJr1GrGxd/DxaWOsH0u1IImEqDKPt3Pn\nVkImv566BcC5a/fYdfgGzz/R3MSRCSGqu5ycHN5//98EBPQG4JtvVpXp/KSkRDZu3IKXV1MA+vV7\nhqFDB/Dbbz8XSiSys7OYPHkqU6e+rN+XkZHO/v2/c/v2LZo08eK3337m1KnjLFz4AQMHDtK369dv\nINOmBbJmzUrefvu9v8WezeLFn9C8eYtyvfbqTgZbiio1qs9jtGvmpN/efSyGM5cTTRiREKKm6Ny5\na7nPbdXKR59EANjZ2dG4sQfJyclF2vbrN7DIuQApKQVtDx7cj6WlFd27P05qaqr+q2FDd1q2bM3h\nw6GFzq9fv0GtTSJAeiREFTMzUzBjaDsWrz9NUmouAGv2XMTNyRZPVzsTRyeEqK5sbGyxsbEp9/lO\nTk5F9imVFuTm5j6yrYWFBQAajQaA2Ng75OfnMWTIgBLvl5eXi5WVdYn3rk0kkRBVzs7GgjkjO/DR\nhjPkqTTkqTQEbY9gYWA37GwsTB2eEKIaMjc3rANdo1EXu//BgEpDKBSlt9VqtTg7u/Dee/8qsY25\n+cN/r2Zm5gbfuyaSRxvCJDxd7Zg6+OGAo6TUXFbtikSj1ZowKiFETWFuXvDPOT9fVWh/SkrlF71z\nd29MWloqHTp0olu3HoW+tFotlpaWKJV153N6tUkkkpKSWLRoEU8++SS+vr4EBAQwb948bt++Xahd\nZmYm//d//0f//v3x9fWlR48ezJw5k4sXLxp8r5ycHL744gv69+9Phw4dGDRoEJs2bUJmwlatLq0b\nMqSnt377z5j7bD943XQBCSFqDCcnZwCuXr2k36fVajlwYF+l3/uJJ/qgVqv57rt1hfZfvXqZ+fNf\nY/36bys9huqkWqRMSUlJjB49mri4OAICAhg0aBA3btxgz549HD58mODgYLy9vcnJyWHChAlcunQJ\nPz8/+vXrR3x8PHv37uXIkSOsXbuWLl26lHovjUbD3LlzCQ0N5cknn2TgwIEcOnSIDz74gDt37vD2\n229X0asWAMN6N+N2YibnrhVUu/z11C2auNnh387dxJEJIaqzPn2eIihoCcuXB5GamoqjY31+//1X\nUlPvV/q9Bw8exr59v7F27Wpu3IimS5fupKQkExKyBRsbG2bOfLXSY6hOqkUiERQURFxcHO+88w5T\npkzR79+1axfz58/nk08+YeXKlWzcuJFLly7x4osv8t57D6fWnDp1ismTJ/P++++ze/fuUu/1888/\nExoayksvvaRPGubOncvUqVNZu3Ytw4cPp3Xr1qVeQxiPmULBtCFt+XBDGHHJ2QCs++USjZxt8TZw\nRVEhRN3j4ODIZ58FsWrVV2zYsBZbW1ueeqofgwcP46WXJlbqvZVKJZ99FsR3361j377fOHr0MA4O\njvj5deWll6bTvPljlXr/6qZaVLbs2bMnOp2Oo0ePFhkQ079/f+Li4oiIiGDMmDFERkYSFhaGnV3h\nEf6TJ0/m+PHjHDp0CDc3txLvNXbsWC5cuMCxY8dwdHTU7w8PD2fs2LFMmjSJd999t9R4pbKl8cWn\nZLN4fRg5eQUDpZwcrFgU2A2HepYmjkwIIQRU48qWGo2GGTNmoFQqix1Va2lpiUqlQq1WM3bsWJKT\nk4skEQ/aQUFVspLk5+dz4cIFfHx8CiURAB06dMDGxobTp09X8BWJ8nB3smXG0LYs3RqBDkhJz2P5\njgvMG+eH0sDR2kIIIaqeyRMJc3NzAgMDiz0WHR3N9evX8fLywtLSklGjRhXbLiUlhbCwMGxtbfH0\n9CzxXrGxsajVary8vIqNw93dnZiYmHK9DlFxHR5zYcSTzdkeWjDg8sqdNL7ff5UXB8ijJiGEqK6q\n7Uc9rVbL4sWL0Wq1jBkzptS2//nPf8jKymLYsGH6nonipKamAmBvX3z3jL29PTk5OajVxc9DFpVv\n0ONN6ebTUL994Gwsh87fNWFEQgghSmPyHoni6HQ6Fi1axPHjx/H19S2xxwJg+fLlhISE4OHhweuv\nv17qdR8kCCUlGw/25+XllToHuEEDW5TK2l1gxJTmT+rG/GWHuXE3HYDv9l6m7WOutGlWu6vDCSFE\nTVTtEgm1Ws3ChQsJCQmhSZMmLF++vMR//EuXLmX58uXUr1+fVatWFRn38N+srKyAguVfi5Ofn49C\noXhkGdb797MNeCWiIl4Z2o4P1oeRmaNCrdHx0dqTLJrcjQb2VqYOTQgh6qSSBltWq0cbOTk5zJw5\nk5CQELy9vdmwYUOxMzA0Gg3vvvsuy5cvx9nZmfXr19OyZctHXv9BopGZmVns8YyMDGxtbctUSlVU\nDpf6NrwyrJ1+ifG0rHyWhVxApdaYODIhhBB/V23+Y6alpREYGEhoaCht27Zl8+bNNG7cuEi7/Px8\nZs2axbZt2/Dw8GDz5s34+PgYdA8PDw8sLCy4c+dOkWMajYb4+HiaNWtW4dcijKONtxMvPPVwxbwb\ncels+O2yVCAVQohqpFokEnl5ecyYMYPz58/TvXt3Nm7ciLOzc5F2Op2ON998kwMHDtCyZUu+//57\nvL29Db6PUqmkY8eOREVFFemViIiIICcnBz8/v4q+HGFE/bp6EuD7sMrl0Qvx7D9TNBEUQghhGtUi\nkViyZAnh4eH4+fmxevXqYutEAGzcuJG9e/fStGnTEh97PMrw4cPJz88nKChIv0+lUrF06VIARo8e\nXb4XISqFQqFg0jOtadbo4bO5H/Zf4+LNyi+DK4QQ4tFMXtkyKSmJvn37olKpGDlyJI0aNSq23eTJ\nk+nfvz/3799nwIABtGrVqth2Y8eOxdXVFYB169aRkZFBYGAgDg4F5ZY1Gg0TJkwgPDycnj170q5d\nOw4fPsylS5cKlc0uPWapbFnVUtJz+WB9GOlZ+UDBUuSLArviUr/0gbFCCCGMo6TBliZPJPbt28es\nWbMe2W7nzp0MHz7coHZt2hQsT/3UU08RGxvL/v37CxWqyszMJCgoiF9++YXU1FS8vLwYN24c48aN\nM2igpSQSpnHtThr/u/ksGm3Br6xXQzsWvNgFKwuZiiuEEJWt2iYSNZEkEqYTei6W9b9e1m93b9OQ\nGUPbofhrdocQQojKUSOmfwrxKE928qCPn4d++9TFRH49ecuEEQkhRN0miYSoccb3a0lLz4fFx7Yd\njObC9WQTRiSEqG1+/nk3vXp15ZtvVpXarlevrowaNaTS4xk1agjPPNOn0u9THpJIiBpHaW7GzOfb\n66tc6oCVu/4kIUUqjgohRFWTRELUSI71LJk9oj0WyoJf4Zw8NV9ujyAnTxZcE0KIqiSJhKixmjVy\nYPIzD6uaxiVns2ZPFFoZPyyEEFVGEglRo/n7ujOgWxP9dvjVe/x45IYJIxJCGCI1M4+9p24R/MdV\nfjoeQ2ItWAzx7NkwevXqytKlnxU59tFH79OrV1euXr1caP+BA/uYNm0S/fv3ZvjwZ1mxIojTp0/S\nq1dXfv55d5Hr3LwZw4IFbzJw4JMMGPAkb7wxu8g1q1q1W/1TiLIa3fcx7iRlEhVTUO3yx6MxeLnZ\n07mVq4kjE0L8N61OR0jodX47dUtfEwZge+h1evq6E/hMayyUdaM2zJYt3/Pll5/h7OzCwIHPoVar\n2L49mAMH9hXbvmA5iSk0aeLFsGEjuHXrJkeOHCIy8gKbNm3F1bVhFb+CApJIiBrP3MyMl4f58sG6\n09xLywVg9Z4o3nuxCx6uxZdbF0KYRkjodX4+cbPYY8ci48lXaZj5fPsqjqpk4eFnHjlzozwSExNY\nuXIZHh6erFjxDU5OBetLjRgxmunTJxd7jkql4tlnBzN//rv6fUuXfsbWrd/zxx+/88ILE4wepyEk\nkRC1gp2NBa+O7MCHG8PIV2nJy9cQFHKBhYFdqWdtYerwhBAUPM747VTpdV/CLidxIy6dZo0cqiiq\n0p07d5Zz584a/bp//PE7+fl5vPjiFH0SAdCqlQ/PPjuE3bt3FHteYOA/Cm0HBPRm69bvuXs31ugx\nGkrGSIhaw7OhHVOfa6vfTryfw6pdf6LVyuBLIaqDU1EJhR5nlOTohbgqiMYwU6ZM48iRsBK/yuvi\nxSgA2rZtV+RYhw4diz3H0tISNzf3QvscHQtq6uTk5JQ7loqSRELUKl19GjK4Z1P9duSNFLaHRpsw\nIiHEA/cz8wxrl2FYu5osLS0VACcnlyLHXFyK7oOCRKIkplztQhIJUesM792cjo897Cr85eQtTkTF\nmzAiIQQUPII0Zrvq5sGaPzqdtsix3NzcQtv16tUDICsrs0jbrKysSoiu8kgiIWodM4WCaUPa4e5k\nq9+37udL3IyXxdaEMKVuPobNKujR1q2SI6kcFhYFCVBxjxliY+8U2m7dumCV6osX/yzSNioqshKi\nqzySSIhaydZayZyR7bGxKphGlq/WsiwkgvTsfBNHJkTd1bCBLT193Utt09LTkTZNG1RRRMbl4dEE\nc3Nzzp4NK9QDcezYEa5cuVSo7YABz2JhYcGGDd+Smpqq33/9ejS7doVUWczGILM2RK3VyLke04e0\n48ttEeiA5PQ8VuyI5M2xnVCaSw4thCkEPtOafJWGsMtJRY619HRkzsgO+kcENU2DBg3o3ftJDh78\ng2nTJuHv34u7d2M5ciSUDh06ERFxTt/W3b0R//jHDFauXMbkyePo1etJ8vJyOXhwP1ZWBesImZnV\njPcpSSRErdaxhQvDn2jOjkPXAbh8O5Xg/deYMKCViSMTom6yUJrzynBfbsRlcDQyjtSMPOxsLOje\n1o02TRtgVkOTiAcWLFiEi4srBw7sY9u2YJo3f4wPP/w/7t69UyiRAJg4cTINGjixZcv3/Pzzjzg4\nODJmzHjq12/A0qWfYm1tbaJXUTYKnSmHetZQSUnyrL0m0el0rNgZWegT0JRnfejdsbEJoxJC1GVp\naaloNJpCNSQe+OabVaxdu5qvv15H27a+JoiueK6u9sXurxn9JkJUgEKh4KXn2uDpWk+/b+Pey0TH\nppkwKiFEXXb2bBhDhw7k22+/LrQ/NTWVn3/ejb29Ay1a1IyeU+mRKAfpkaiZElNzWLzuNFm5BUuN\nO9pZsiiwGw3srUwcmRCirsnOziYwcCwJCfE8/nhPmjdvQUZGOocOHSQtLZV3332fgQMHmTrMQkrq\nkZBEohwkkai5/oxJYUnwOR781j/W2IH54ztjoZTOOSFE1bp37x6bN6/n2LEjJCYmYmtrQ+vWbRk/\n/kW6dOlm6vCKkETCiCSRqNn2nrrFD39c02/37tCIyc/61NiR4kIIURVkjIQQf+nfrQn+7R4WvDkc\nEccfZ0234I0QQtRkkkiIOkehUBD4jA/e7g+z6x/2X+XyrfsmjEoIIWomSSREnWRpYc7sEe1xsC0o\naavR6li+M5LktNxHnCmEEOLvJJEQdZaTgzUzn2+PuVnB2IiMbBVBIRHkqTQmjkwIIWoOSSREndaq\nSX3G9384V/tWQibrf7lk0iV5hRCiJpFEQtR5ff08eLLTwyqXJ6IS+O3UbRNGJIQQNYckEkIAE/q3\nooWno35768FrRF5PNmFEQghRM0giIQSgNDdj1nBffZVLnQ5W7vqThPvZJo5MCCGqt2qTSCQlJbFo\n0SKefPJJfH19CQgIYN68edy+XXIXc3Z2Nk899RQfffRRme711ltv0bp162K/Pv3004q+FFFDOdpZ\nMXtEe/0S49l5apZtv0BOntrEkQkhRPVVLZYRT0pKYvTo0cTFxREQEMCgQYO4ceMGe/bs4fDhwwQH\nB+Pt7V3oHLVazbx584iNLXshoUuXLuHi4sLYsWOLHOvSpUt5X4aoBZo1ciDwmdZ889NFAGLvZfHN\nTxeZ+bxvjV/eWAhRMQ9W5fznP/+HQYOGmDqcaqNaJBJBQUHExcXxzjvvMGXKFP3+Xbt2MX/+fD75\n5BNWrlyp35+amsobb7zB0aNHy3wvlUrFjRs36NOnD3PmzDFK/KJ2CWjfiJsJGewLuwPA2StJ7DkW\nw9CAZiaOTAghqp9q8Whj3759ODk5ERgYWGj/sGHD8PLy4siRI2i1WgD27NnDoEGDOHr0KAEBAWW+\nV3R0NCqVitatWxsldlE7vfBUC9o0baDf3nn4BuFXk0wYkRBCVE8mTyQ0Gg0zZsxg9uzZmJkVDcfS\n0hKVSoVaXfCcOjg4GGtra1auXMmMGTPKfL/Lly8DSCIhSmVuZsbLw9rh4mit37d6dxR372WZMCoh\nhKh+TP5ow9zcvEhPxAPR0dFcv34dLy8vLC0tAZg5cyadO3fGysqKkydPlvl+DxKJGzduMHbsWC5f\nvoy1tTV9+vThtddew83N7RFXEHWFva0ls0e059/fnSFfpSU3X0PQ9ggWBnbF1trC1OEJUaOl5aVz\nJuEcqXnp1LOwpXPDjrjaOps6rDJLTr7H2rVrOHr0EKmp93FxcaVv335MnvwPbG3r6dup1Wo2bPiW\n0NA/iI29g4WFJW3atGX8+El07dq9zO0AsrIy2bBhLQcO7CMpKRFHx/r06vUE//jHDBo0cKqyn4HJ\neyRKotVqWbx4MVqtljFjxuj3+/v7Y2VlVe7rPkgkli9fjqenJy+88ALe3t6EhIQwevRo4uPjKxy7\nqD283Ox5aVAb/XbC/RxW/RiFViuVL4UoD61Oy67oX3jv2L/Zfm0P+28f4sfrv/L+if9lQ1QwKo3K\n1CEaLD4+nmnTAtm1azutW/swZsx4mjRpyubNG5g9ezo5OTn6tl988R++/fZr7O0dGDFiDE891Y+o\nqEjefHMOZ8+GlbldZmYmr7zyDzZtWk+jRh6MHj0WX9/2/PjjDqZNC+TevXtV9nMoU4+ERqMhPDyc\nU6dOER8fT3JyMhqNBkdHR7y9vfHz86NLly5YWFTs05pOp2PRokUcP34cX1/fEnssysPa2hpvb2+W\nLVtGy5Yt9ftXrFjBF198wYcffsiyZctKvUaDBrYoleZGi0lUb8+52pOcmc/W/VcBuHA9mV/D7hD4\nXFsTRyZEzbM5Yid7bx4o9tjJ+DMolDreCJhWxVEZpl69gg+x9vbWuLra8+67b5KUlMjKlSvp06eP\nvt2GDRv46KOP+OGHdcyfP5/MzEx+/HEH3bp147vvvtO3u3BhPKNGjeKnn3YwcGBfg9sBLF++hOvX\no1m0aBETJkzQt92/fz8zZ85k5covWLp0aSX/RAoYlEjcvHmTtWvXsnv3brKzs0tch0ChUODg4MCw\nYcP4xz/+Ua7HBGq1moULFxISEkKTJk1Yvny5/rGGMXz11VfF7p8xYwbbtm3jwIEDZGVlUa9evWLb\nAdyXIkV1zsAunly6kcKFv6pdbvvjKi72lnRvI4/ChDBUWl46P176vdQ2J+6cJSw6iqYOTaooKsNl\nZeUBkJGRy8WLNzh06BD+/gG0a9eFpKQMfbuBA4exevUatm8PYcqUV8jIyECn03Hnzl0uXbqBs7ML\nAO7u3gQH78TNzZ2kpAyD26nVanbs2EmzZs0ZMGBooXt36NCd9u078vvvvxMTE0e9enZGe/2urvbF\n7i81kcjIyODTTz9l+/btqNVqGjZsSJ8+fWjdujWenp7Y29uj1Wq5f/8+iYmJhIeHc/bsWTZs2EBw\ncDCjR49m7ty52NsXf/P/lpOTw9y5cwkNDcXb25u1a9dW2ZgFMzMzfHx8uHPnDvHx8Tz22GNVcl9R\nM5iZKZgxtC2LN5whIaUgkfz254u4O9ni5WbY77cQdd2ZhHNoddpHtjsRd6ZaJhJ/d+VKweJ+aWlp\nfPPNqiLHLSwsSExMICkpEVfXhjz1VH/279/LyJGDad++I48/3pOePXvTrFlz/Tn29vYGtbt16yY5\nOdlotdpi752fn49GoyE6+hodOnSqnB/A35SYSBw+fJgFCxaQlZXF6NGjGTVqFO3atTPoosePHyc4\nOJjg4GD27dvHxx9/jL+/f6nnpKWlMW3aNM6fP0/btm1Zs2YNzs7GHXiTm5vLlStXsLS0xMfHp9jj\nQIXGYIjay9bagjkj2vPhhjBy8zXkq7QEbb/Aosldsbc1Xq+ZELVVal66Qe3S8tIqOZKKy8ws6AX4\n888L/PnnhRLbpaen4+rakIULP8DHpy0///wj4eFnCA8/w4oVQfj4tOXtt9+lZcuCmYSGtHtw75s3\nY1i7dnWp964KJSYSr7zyCiNHjuTVV18t8z90f39//P39SUhI4PPPP2fatGlERkaW2D4vL48ZM2Zw\n/vx5unfvzooVK7CzM153zAOJiYmMHj2aVq1asXv37kLHcnJyiIqKwsnJCQ8PD6PfW9QOjV3qMW1I\nW4K2F7xxJKfnsmJnJG+80ElfWlsIUbx6FrZGbVdZYmPvcP58OL6+HfDyaqrf/+CxvpWVlf6R++TJ\nU5k69eVHXlOpVDJu3ETGjZtIfHw8YWEn+OOPfZw6dYL5819n69YfUSqVBrWzsSn4+QwcOIiFCz+o\nhJ9A2ZT4zrdt2zb+9a9/VahXwM3NjU8++YStW7eW2m7JkiWEh4fj5+fH6tWrKyWJAPDy8sLHx4cr\nV67w448/6vfrdDo+++wzUlJSGDduHAophSxK4dfSleG9H1a5vHQrlS0HrpkwIiFqhs4NOxrUrotb\n5XfHl+bw4YP8+9//4sKF84X2Z2QUfMJ3cHDksccKButfuhRV7DXWrVvDpk3rUalU3LoVw/LlX3L0\n6GEA3N3dGTx4OEuWLKNLl24kJSVy926swe28vJpiaWnJ5cuXih2zuH37Ftav/4a0tFSj/UxKU2KP\nRHFd/+XVpk2bEo8lJSWxadMmAJo3b87q1cV300yfPr3MjxzWrVtHRkYGgYGBODg4APCvf/2LwMBA\n5s+fz969e/Hw8CAsLIzIyEi6devGyy8/OrMUYnBPb24nZHLmSkG1y31hd2jqZk9A+0YmjkyI6svV\n1pke7l04GX+mxDaPOXrTukGLKoyqqFatCv7/nThxjOeeGwoULK8QFnYKc3NzfHzaYm9vT6dOnTlx\n4hgHDuyjb99++vP37v2FNWtW0qlTZyZMCESptGDz5g2cPn2Cbt166HszVCoVycn3sLS0xNnZmbS0\nNIPaWVlZ8dRT/fn115/44YdNjBs3UX/v8PAzfPnlZzRs6M6LLz5ccqIyGaUgVUREBHfv3qVdu3Y0\naVK2ATLnz59HpSqYN7x9+/YS2wUGBpY5kdiwYQOxsbE8//zz+kSiU6dObN26lWXLlnHy5EkOHjyI\np6cnc+fOZerUqUadISJqLzOFgpeea0N8Sjaxf1W7XP/rZRo516N5YwcTRydE9TWu9QjyNfmEJxUd\nV/CYozfTOwSavFe4U6fOtGnTlgMH9vHGG7Np0aIVp06d4ObNGEaNGqufQDB//j+ZOXMaCxe+w+OP\n96R588e4desmx44dwcHBkTfffAeAxo09GDZsBLt2hTBp0gv4+/fCzEzByZPHiYm5wZQp06hXz456\n9ewMagcwa9ZrREZG8NVXX3DkSCht2/qSmJjAoUMHMDdXsmDBwmKrRVcGha6kuZzFOH78OCtWrGDq\n1Kk88cQTALz++uv8+uuvQMHMh6lTp/L6669XTrTVxN+n2oi6LfF+NovXh5GVW1DCvb6dJf8zuRuO\ndjJgV4iS6HQ6bmbc5mTcGdL+qmzZxa0TrRo8hpmieow1SktLZfnyLzl27AiZmRm4uzdi8OBhjBv3\nYqF/0AkVX2FSAAAgAElEQVQJ8axb9w0nThwlNfU+zs4u+Pl1YfLkqXh4eOrbqdVqdu7cxk8//Uhc\n3F00Gg3Nmj3GiBGjeeaZ58rcDiA9PY0NG9Zy6NABkpISqV+/Ae3a+RIYOJWWLVsZ/WdS0vRPgxOJ\nc+fOMXHiRDQaDQsWLGDSpEn8/vvvzJkzBxsbG5544gnCwsJISUnhiy++YODAgUZ9AdWJJBLi7yJv\nJPP5lvM8+Etq4eHIW+P8sFBWjzdEIYQwhpISCYPf6b799lt9sagHVbR27NiBQqHgf/7nf1i6dCnb\nt2/H2tqaH374wThRC1ED+DZzZnSfh890r8Wmsen3KyUWbhNCiNrE4ETi7NmzdOjQgfHjx2Nubk5+\nfj7Hjx9HqVQyYMAAoGCEaefOnfnzzz8rLWAhqqOB3ZvweLuHxdMOnb/LwXN3TRiREEJUDYMTibS0\nNBo3bqzfDgsLIycnhw4dOmBr+3DOr62tLdnZUkJa1C0KhYLJz/jQ9G9VLjf/foUrt6tm+pUQQpiK\nwYmEm5tboZUxDxw4gEKhICAgQL9Pp9Nx6dIlXF1djRulEDWApYU5s0e0x962YNE6jVbH8h0XSEnP\nNXFkQghReQxOJNq3b8/58+fZvn07R48eZefOnQD6xxr5+fl8+umn3Llzhx49elROtEJUc86O1swc\n7ou5WcH0tfRsFUEhF8hXaUwcmRBCVA6DZ21ER0czceJEUlMLump1Oh2DBg1iyZIlAPTu3Zt79+7h\n6OhIcHAwTZs2Le1yNZrM2hCP8sfZO3y394p+27+dG1MHtzX5/HghhCivCk//BLh9+zYrVqzg3r17\ndO3alSlTpmBhUdCNO3HiRJycnJg/fz6enp6PuFLNJomEeBSdTse6Xy5xOCJOv++Fp1owsLuXCaMS\nQojyM0oiURqdTldnPm1JIiEMoVJr+b/vzxIdW1CfX6GAN17oRDtvJxNHJoQQZVfhOhJ/p1KpiIiI\n4LfffuPMmYKa6X8fiCmEAAulGbOeb4+jXUHZdZ0OVu6MJDE1x8SRCSGE8ZQpkVCpVHz++ef07NmT\nF154gddee43g4GAA5s+fz4gRI7h582alBCpETVTfzorZI9qjNC/orcvKVRO0PYLcfLWJIxNCCOMw\nOJFQqVRMmzaNr7/+GpVKhZ+fX6HKfWq1mqioKCZMmEBiYmKlBCtETfRYY0deHNhavx2blMU3P12U\nypdCiFrB4ERi48aNnDhxgr59+/LHH3+wefPmIsfHjRvHvXv3WLNmjdEDFaIm692hMU93eTgI+czl\nJPYcl947IUTNZ3AisXPnTpydnVmyZAlOTkUHiymVSt577z0aNWrEkSNHjBqkELXBC0+1wMervn57\n56HrnLt2z4QRCSFExRmcSNy8eZMuXbpgbW1dYhtzc3N8fX25e1fWGBDivynNzXh5uC/ODgV/Qzpg\n9e4/iUvOMm1gQghRAQYnElZWVqSkpDyyXVJSElZWVhUKSojaysHWkjkj22P51xLjOXkagrZfIDtX\nBl8KIWomgxMJX19fIiIiiI6OLrHNlStXiIyMxNfX1yjBCVEbebnZM2VQG/12fEo2X+/+E61WBl8K\nIWoegxOJyZMnk5+fz/Tp0zl48GChFT51Oh3Hjh1j1qxZaDQaxo8fXynBClFb9GjrxrM9Hla5jIhO\nZueR6yaMSAghyqdMlS2/+uorgoKCClWwtLa2RqPRoFKp0Ol0TJ48mXfeeadSgq0upLKlMAatVscX\n284Tef3hI8OZw33p6tPQhFEJIUTxjFYi+9ixY3zzzTecOXOG3NyC5ZGVSiWdOnUiMDCQ/v37Vzza\nak4SCWEsWbkqFq8PI/F+QbVLSwsz3n2xK00a2pk4MiGEKMzoa21otVpSU1PRaDTUr19fv3hXXSCJ\nhDCm2KRMPtx4hrz8gqXGXRytWTS5G3Y2dedvSghR/VV4rY0FCxawdevWhyeameHk5ISrq2uhJOLr\nr79m8uTJ5Y9UiDrGw9WO6YPb6rfvpeWyYmckGq3WhFEJIYRhDE4kduzYwenTpx/Z7uTJk5w9e7ZC\nQQlR1/i1cmVYr2b67Ys377P1QMkzpIQQorpQlnRg3rx5RdbMOHbsGJMmTSrxYpmZmVy8eJHGjRsb\nL0Ih6oghAd7cSsgg/GpBtcu9p2/TpKEdAe0bmTgyIYQoWYmJRN++fXnzzTf12wqFgnv37nHvXukl\nfZVKJXPmzDFehELUEWYKBVMHt+WjjWe4e6+g2uX6Xy/T2KUezRo5mDg6IYQoXqmDLcPDw9Fqteh0\nOiZOnEivXr145ZVXir+QQoGVlRWenp44OjpWWsDVgQy2FJUp4X42i9eFkZ1XUO2ygb0ViyZ3w7Ge\npYkjE0LUZRWetbFgwQI6d+7M6NGjjRpYTSSJhKhsF64n88XW8zz462zh6cj8cX4ozQ0e1iSEEEZl\n9OmfdZkkEqIq/HLiJlsPPhxw2cfPg0kDW5swIiFEXVZSIlHiGIn/tnPnzjLdcPjw4WVqL4Qo7Jke\nXtxMyODUxYJBzwfDY/Fys6NPJw8TRyaEEA8Z3CPh4+NTqDR2SXQ6HQqFgosXL5YpkKSkJIKCgggN\nDSU5ORlHR0f8/f2ZO3cuTZo0Kfac7OxsBg8ezNNPP827775r8L1ycnJYtWoVP/30EwkJCXh6ejJh\nwgTGjx9v0GuUHglRVfJUGj7eeIZbiZkAmJspmD/ej5ae9U0cmRCirqlwj8TIkSOL/Ser0WjIyMjg\nwoULJCQkMGjQILp3716m4JKSkhg9ejRxcXEEBAQwaNAgbty4wZ49ezh8+DDBwcF4e3sXOketVjNv\n3jxiY2PLdC+NRsPcuXMJDQ3lySefZODAgRw6dIgPPviAO3fu8Pbbb5fpekJUJisLc2aPbM8H68LI\nzFGh0er4akckiwK74uRgberwhBDCeGMktFotn332GRs2bOD7778v01LiixYtIjg4mHfeeYcpU6bo\n9+/atYv58+fTt29fVq5cqd+fmprKG2+8wdGjRwGYNGmSwT0Su3fvZt68ebz00kv6pEGlUjF16lRO\nnjzJrl27aN269OfQ0iMhqtqlm/f59IdzaP/6c23WyJ53JnTGQmlu4siEEHVFhUtkP4qZmRnz5s3D\nzc2NpUuXluncffv24eTkRGBgYKH9w4YNw8vLiyNHjqD9q1zwnj17GDRoEEePHiUgIKDMcW7atAml\nUsnLL7+s32dhYcFrr72GTqdj27ZtZb6mEJXNp2kDxj7dQr99Iy6D9b9eRsZKCyFMzeBHG4ZQKBS0\nbduWY8eOGXyORqNhxowZKJVKzMyK5jWWlpaoVCrUajWWlpYEBwdjbW3NypUrsbW11fdKGCI/P58L\nFy7g4+NTpNZFhw4dsLGxMagMuBCm8HQXT24lZHLkQhwAxyLjaepmT/9uxY8hEkKIqmDURAIgOjq6\nTJ+SzM3Ni/RE/P1a169fx8vLC0vLgmI8M2fOpHPnzlhZWXHy5MkyxRYbG4tarcbLy6vYONzd3YmJ\niSnTNYWoKgqFghcHtuJuchbX76YDEPzHNTxc69HW28nE0Qkh6iqDH20kJCSU+HX37l0iIiJ46623\nuH79Op06dapwYFqtlsWLF6PVahkzZox+v7+/P1ZWVuW6ZmpqKgD29sU/57G3tycnJwe1Wl2u6wtR\n2SyU5sx6vr2+yqVWp2PFzkiSUnNMHJkQoq4yuEfiySeffOTUSJ1Oh1KpZPbs2RUKSqfTsWjRIo4f\nP46vr2+JPRZl9SBBeNC78d8e7M/Ly0OpLPlH06CBLUoZ5CZMxNXVnvde6sGC5UdRa7Rk5apZsetP\n/jOnN9ZWRu9kFEKIUhn8rlPaip5mZmbY2trSqlUrXnzxRTp27FjugNRqNQsXLiQkJIQmTZqwfPny\nEv/xl9WDngyVSlXs8fz8fBQKBTY2NqVe5/79bKPEI0R5OdezYOKAVqz75RIAMXHp/O+G07wyrJ1B\ntVCEEKKsKlxH4o8//jBaMCXJycnR13jw9vZm7dq1uLm5Ge36DwZYZmZmFns8IyMDW1vbYgd9ClHd\nPNGxMTcTMjhwtqCWStilRH52s+M5f2/TBiaEqFOqzX/MtLQ0AgMDCQ0NpW3btmzevLnUXpDy8PDw\nwMLCgjt37hQ5ptFoiI+Pp1mzZka9pxCVadzTLWnV5GGVy5DQ60RE3zNhREKIuqbMD1Szs7OJiYkh\nOzu71NkZ3bp1M/iaeXl5zJgxg/Pnz9O9e3dWrFiBnZ1dWUN7JKVSSceOHYmIiCAzM7PQPSIiIsjJ\nycHPz8/o9xWisijNzZg53JcP1p8mJT0PHbDqxygWBnbF3cnW1OEJIeoAgxMJjUbDJ598wg8//PDI\nWQ0KhYKoqCiDg1iyZAnh4eH4+fmxevVqrK0rr/Tv8OHDCQsLIygoiAULFgAFYyYeFNGSZdJFTeNQ\nz5LZI9rz8XdnUam15OSpCdoewXuTumIjgy+FEJXM4HeZtWvXsnHjRgCaNWuGi4tLqTMbDJWUlMSm\nTZsAaN68OatXry623fTp08s87XPdunVkZGQQGBiIg4MDACNGjGD79u2sW7eOK1eu0K5dOw4fPsyl\nS5d46aWXHlkeW4jqyNvdgSnP+vD17oIEPi45m9W7o5g9sj1mMvhSCFGJDM4EQkJCUCqVrF69Gn9/\nf6MFcP78ef0siu3bt5fYLjAwsMyJxIYNG4iNjeX555/XJxLm5uasWbOGoKAgfvnlF86cOYOXlxeL\nFi1i3Lhx5X8hQpjY4+3cuZWQya+nbgFw7to9dh2+wfNPNDdxZEKI2szgRbs6dOhAt27d+Oabbyo7\npmpPFu0S1ZVGq+WLLef5M+a+ft+s533p0rqhCaMSQtQGFV60q379+jI/XYhqztzMjBnDfHGt/3Cc\n0Zo9F7mTVPyUZyGEqCiDE4mnn36a8PBw7t2TqWVCVGd2NhbMGdkBK4uC6qt5Kg1B2yPIzCm+EJsQ\nQlSEwYnE66+/jru7OzNmzODEiRNkZ0t1RyGqK09XO6YObqPfTkrNZdWuSDRarQmjEkLURgaPkRg1\nahRZWVmFVsc0MzMr8XFHZGSkUQKsjmSMhKgpdhy6zu5jMfrtZ7p7MeapFqYLSAhRY1W4RHZxiYFG\noyl/REKISjesdzNuJ2Zy7lrBI8lfT92iiZsd/u3cTRyZEKK2MLhHQjwkPRKiJsnJU/PhhjDikgse\nR1oozVgwsTPe7g4mjkwIUZNUeNaGEKJmsrFSMmdkB32VS5Vay7KQC6Rn5Zs4MiFEbSCJhBB1gLuT\nLTOGtuXBiKaU9DyW77iAWiODL4UQFVPiow1fX18A9uzZg7e3t37bUDLYUojq56fjMWwPva7f7tvZ\ngxcHSFl4IcSjlXmw5YOFuR7kGY9aqEsIUf0NerwptxIyOX0pEYADZ2Np6mbPEx0bmzgyIURNJYMt\ny0F6JERNlpev4d/fneF2YkG1S3MzBW+P70wLT0cTRyaEqM5ksKUQAgArS3PmjGiPnY0FABqtjq92\nXOB+Rp6JIxNC1ETl6pFITU0lL6/0Nx03N7dyB1XdSY+EqA0uxqTwWfB5tH+9BTRr5MA7E/ywUJqb\nODIhRHVUUo+EwYmERqPh888/Z+vWraSnp5faVqFQEBUVVfYoawhJJERt8fvp23y//6p+O6C9Oy8N\naiML9AkhiqhwZcs1a9awZs0aAOzs7HBwkGI2QtR0/bp6cishg6OR8QAcvRBPUzd7+nVtYuLIhBA1\nhcGJxI4dO7C0tGTVqlX4+/tXZkxCiCqiUCiY9Exr7iZncSOuoKfth/3X8HC1o03TBiaOTghRExg8\n2DIuLo6ePXtKEiFELWOhNGfW8+1xqGcJgFanY8XOSO6l5pg4MiFETWBwIuHl5UVWVlZlxiKEMBEn\nB2tmPe+LuVnB2IjMHBXLQi6Qp5KF+YQQpTM4kXjxxRc5c+YMZ8+ercx4hBAm0tKzPhMGtNJv30rM\nZO3PF5FSM0KI0hg8RmLMmDHExMQQGBjI0KFD8fHxoX79+iW2HzJkiFECFEJUnT6dPLiVkMnB8FgA\nTl1MpKmbPc8+3tTEkQkhqiuDE4mUlBTCwsJQqVSEhIQ8sr0kEkLUTOP7tSQ2KZOrd9IA2HYwGs+G\ndrRv7mziyIQQ1ZHBicTHH39MREQENjY2dO/enQYNGshccyFqIaW5GTOfb88H605zPyMPHbBy158s\nCuyKm5OtqcMTQlQzBhekevzxx1EqlezYsQNXV9fKjqtak4JUoi64EZfOJ5vOolIXLDXeyNmW9yZ1\nxcbK4M8fQohapMJrbeTl5eHn51fnkwgh6opmjRwIfObhEuNxydms2ROlL6kthBBQhkTC19eXmJiY\nSgxFCFHd9PRtxIBuD6tchl+9x49HbpgwIiFEdWNwIjFnzhyio6NZtmwZWq22MmMSQlQjo/s+Rlvv\nh1Uufzwaw9krSSaMSAhRnRg8RiI4OJh9+/Zx5MgRXFxc8PX1pX79+iiVRZ+XKhQKPvjgA6MHW13I\nGAlR12TmqPhg3WnupeUCBUuRv/diFzxc7UwcmRCiqlR49U8fHx8UCoVBxWkUCgUXL14sW4Q1iCQS\noi66nZjJRxvDyFcV9Eg2bGDDwsCu1LO2MHFkQoiqUOFEYseOHWW64fPPP1+m9jWJJBKirgq7lMjy\nnZH6bd9mTrw2uiNmZjIVXIjarsKJRGVLSkoiKCiI0NBQkpOTcXR0xN/fn7lz59KkSeEljXfu3Mm6\ndeuIiYnBwcGBZ599lldffZV69eoZdK+33nqLH3/8sdhj06ZNY968eY+IVRIJUXeFHIpmz7Gb+u1n\ne3gxum8LE0YkhKgKZU4kUlJScHJyMsrNk5OTcXYuuSpeUlISo0ePJi4ujoCAAFq3bs2NGzc4ePAg\njo6OBAcH4+3tDcCqVatYsmQJrVu35oknnuDKlSuEhobi5+fHhg0bsLS0fGQ8Q4YMISUlhbFjxxY5\n1qVLF3r27Fnq+ZJIiLpMq9Px5bYIIqKT9fumD23L423dTRiVEKKylTmR6NatG7NmzWLixInFDqg0\nRH5+Pt9++y2rV6/mzJkzJbZbtGgRwcHBvPPOO0yZMkW/f9euXcyfP5++ffuycuVKYmNjGTBgAO3b\nt2fjxo1YWBQ8m126dCnLly9n4cKFTJw4sdSYVCoVfn5+9OnTh2XLlpXrdUkiIeq67Fw1H24IIz4l\nGwBLpRkLJnahqXvxbzRCiJqvzAWp5s6dy+eff86AAQNYv349aWlpBt8sLi6OL7/8kn79+rFq1apH\nPirYt28fTk5OBAYGFto/bNgwvLy8OHLkCFqtli1btqBWq5kxY4Y+iQB4+eWXsbOzY+vWrY+MLTo6\nGpVKRevWrR/ZVghRPFtrJXNGtsfGyhyAfLWWZSERpGfnmzgyIURVK7GrYeLEiQQEBPCvf/2Ljz/+\nmP/85z907dqVrl274uPjg4eHB3Z2dmi1Wu7fv09iYiLh4eGEhYURFRWFRqMhICCARYsW0bRpySsH\najQaZsyYgVKpxMysaF5jaWmJSqVCrVZz+vRpALp3716ojZWVFZ06deLIkSNkZGRgb1/yp6LLly8D\nSCIhRAU1cq7HtCHtCNoWgQ5ITs9jxY5I3hzbCaW5wSVqhBA1XKnPLJo1a8a6des4cuQIq1ev5sSJ\nE5w4caLExbp0Oh1mZmYEBAQwc+ZM/Pz8HhmAubl5kZ6IB6Kjo7l+/TpeXl5YWlpy69YtXFxcih1U\n6eHhAcCNGzfo0KFDifd7kEjcuHGDsWPHcvnyZaytrenTpw+vvfYabm5uj4xZCFGgUwsXhj/RnB2H\nrgNw+XYqwfuvMWFAKxNHJoSoKgYNfujVqxe9evUiLi6Oo0ePcvr0ae7evcv9+/dRqVTUr18fb29v\n/dgDd/eKD7rSarUsXrwYrVbLmDFjAEhNTcXT07PY9g96ITIzM0u97oNEYvny5fTv359OnTpx/vx5\nQkJCOHr0KFu2bDFK/ELUFYP9m3I7IYOwywXVLvefvYOXmx29OzY2cWRCiKpQplGUjRo1YtSoUYwa\nNaqy4gEKejYWLVrE8ePH8fX11fdYqNXqEmdlPNifl5dX6rWtra3x9vZm2bJltGzZUr9/xYoVfPHF\nF3z44YePHITZoIEtSqV5WV6SELXa/MDuzA86TExcOgAb916hbUtXfJoaZ+aXEKL6qnbrAavVahYu\nXEhISAhNmjRh+fLl+iTB2toalUpV7Hn5+QWDvGxsbEq9/ldffVXs/hkzZrBt2zYOHDhAVlZWqTUp\n7t/PNuSlCFGnvDysHYvXnSYrV41ao+XDb0+yKLAbDeytTB2aEMIIKryMeFXIyclh5syZhISE4O3t\nzYYNGwqNWXBwcCAjo/iplw/2lzbQsjRmZmb4+PigVquJj48v1zWEqMsa1rfh5eG+PBhClZaZz/Id\nF1CpZZE/IWqzapNIpKWlERgYSGhoKG3btmXz5s00blz4Gau3tzfJycnk5uYWOT82NhYzM7NSZ4jk\n5uYSERHBpUuXSjwOBbNAhBBl187biRf+VuUy+m463+29bNAaPUKImqlaJBJ5eXnMmDGD8+fP0717\ndzZu3FhsJcwuXbqg1WoJCwsrcv65c+do0aIFdnYlr0aYmJjI6NGjeeutt4ocy8nJISoqCicnJ/0M\nECFE2fXv1gT/dg97Eg9HxPHH2VgTRiSEqEzVIpFYsmQJ4eHh+Pn5sXr16hKTgcGDB2Nubs6yZcv0\nYyIAVq5cSWZmJi+88EKp9/Hy8sLHx4crV64UWmtDp9Px2WefkZKSwrhx40qc3iqEeDSFQkHgMz54\n/63K5Q/7r3L51n0TRiWEqCwmX7QrKSmJvn37olKpGDlyJI0aNSq23fTp07GysuLTTz9l9erVPPbY\nY/Tt25dr165x8OBBOnfuzPr16wvN6li3bh0ZGRkEBgbi4OAAwLlz5wgMDCQvL49+/frh4eFBWFgY\nkZGRdOvWjW+//faR63VIiWwhHi0lPZcP1p0mPbtggLS9rQWLArvh7Ght4shEXafWaDl7JYnwq/fI\ny9fQsIENvTo0wtO15B5tUY1X/9y3bx+zZs16ZLvTp0/j4OCATqdj8+bNbN68mVu3buHq6kr//v2Z\nPXt2kYGWTz31FLGxsezfv79Q/YkrV66wbNkyTp48SVZWFp6engwdOpSpU6catOiXJBJCGObK7VT+\n8304Gm3B24yXmx0LJnbBykKmTwvTSLifzRdbzpNwP6fIsae7eDKuX0vMpFe6WEZLJNRqNXv37uXU\nqVPEx8fTuXNnpk+fzrZt2/D19cXHx8coAVdnkkgIYbgD4bFs/O2yfvvxtm5MG9JWHiGKKpebr2bR\nN6e4l1Z0wP4DQwO8Gd67eRVGVXOUlEiUqY5EZGQkr732GrGxseh0OhQKhf6RwQ8//EBUVBRvv/12\niSWvhRB1T18/D24lZBB67i4AJ6IS8HKz55keXiaOTNQ1xyLjS00iAH47fZuB3b2wsap2ZZaqLYMH\nW965c4eXXnpJv5T34sWLC03p6tevH5aWlnzyySecPHmyUoIVQtRME/q3ooWHo35768FrRF5PNmFE\norbTanWkpOdy9U4qJ/6M56fjMew+GvPI8/LyNVyQ380yMTjl+uqrr8jIyOCTTz5h2LBhACxcuFB/\n/OWXX6Zjx45MmTKFtWvX0qNHD+NHK4SokZTmZsx63pcP1odxPyMPnQ5W7vqThZO74tbA1tThiRpI\npdaSkpFLctpfX+kPv99Ly+V+Rp5+bE5ZZeUUX0FZFM/gROLIkSO0adNGn0QUx9/fn06dOnHx4kWj\nBCeEqD0c7ayYPaI9H393FrVGS3aemmXbL/DPF7tIN7IoIidPTXJaLvfSiyYKyWm5pGXlP/oi5dTA\nXmYWlYXBf72pqakGLQvu4uJCVFRUhYISQtROzRo5EPhMa775qeDDRuy9LL756SIzn/eVkfJ1iE6n\nIyNbpU8K7hWTKGTnqY1yL2tLc5wdrXF2sEat0RIVU3o9E4d6lvg2l8XmysLgRMLV1ZVr1649st3V\nq1dxcXGpUFBCiNoroH0jbiZksC/sDgBnrySx51gMQwOamTgyYSwarZb7GXnF9iTcS88jJT3XaGuw\nONha4ORgrU8WnB2tcXmw7WiNrZVSP0NIrdHyv5vOEn03vcTrjXyiOUrzalGrscYwOJHo3bs3W7Zs\nYfPmzYwfP77YNps2beLmzZuMGTPGaAEKIWqfF55qQWxSFhdvFnw63Hn4Bk0a2uHX0tXEkQlD5Ks0\nhROE9MJjFe5n5KM1QokiM4WCBvaW+gTh78mCs0PBl2UZapIozc14bUxHVu+OIiK68IBKK0tzxvRt\nQe+OjUs4W5TE4DoSCQkJDBs2jLS0NPr370+PHj1YvHgxAQEBDB8+nNDQUH766Sfs7OzYsWNHoQJQ\ntY3UkRCi4jKy81m8Pkw/Hc/a0pz3JnWlsUs9E0dWt+l0OrL/Gp9Q0hiFjGzjDEa0UJrh5GCNi4NV\n0STB0ZoG9laYm1VO78CdxEzCryaRq9Lg1sCWbj4NZazOIxilIFVUVBSvvvoqd+7cQaFQ6GtJQMEv\nn4uLC1988QVdu3Y1TtTVlCQSQhjHrYQM/v3dGfJVBd3cbg1sWBjYFVtrCxNHVntpdTrSMvOL71H4\n63tuvsYo97K1UhZJEFz+1rNgb2shhclqEKNVtlSpVOzdu5eTJ08SHx+PRqOhYcOGdO3aleeeew5r\n69o/2lUSCSGM59TFBFbu+lO/3b65M3NHdcDMTP7BlIdaoyXlwfiEYsYopGTkotYYZ2UERztL/SOG\n4sYoyCf82qXarrVRE0kiIYRxbTsYzc8nbuq3Bz3elFF9HjNhRNVXbr76vwYvPkwUUtLzSM3Iwxhv\n6uZmChrYWxX0IPxXouDsaI2TvTUWShmUWJcYpUT2A6mpqeTl5ZXaxs3NrTyXFkLUQSOeaM7txEx9\nRcGfT9zEy82O7m3q1vuITqcjM0dVaPDif49RyMo1zrRISwuzorMc/va9vp2V9AoJgxjcI6HRaPj8\n88FN5aIAACAASURBVM/ZunUr6eklT50BUCgUtbqWhPRICGF82bkqFq8P06/KaGlhxj8ndsHLrfhP\nQTWRVqsjNTOv2LoJD8YqPBgvUlF2NhZFexL+NkahnrVSxieIMqnwo41Vq1bx+eefA2BnZ6dfrKsk\nf/zxRxlDrDkkkRCicty9l8WHG8L0g/2cHaxZNLkr9raWJo7MMCq1hpT0vIe9CP+VKFSkbPPfKYD6\n9lbFJgoF362wtpTxCcK4KpxIPPPMM9y9e5dVq1bh7+9v1OBqGkkkhKg84VeTCNp+Qb/t41WfN8d2\nqrRpgGWRnasupsDSw+10I5VtVporCooslfDooYG9lRRNElWuwolEx44d8ff3Z+XKlUYNrCaSREKI\nyvXj0RvsPHxDv+3X0gUbKyXpWfk41rPkcV932jRtYNSy2jqdjvSsfH1ikJL+sDLjg0cROZVQtrm4\nRMGhnqWUDBfVToUHW3p5eZGVlWW0gIQQoiSDe3pzKyGTs1eSAAi/eq/Q8aOR8bRp2oBZz7fH1tqw\ntzG1RktqRl6hxKDwo4c81BrjlW0u/pFD0bLNQtR0BvdIbNmyhffff5/vvvuOzp07V3Zc1Zr0SAhR\n+bJzVbz51VHyShl82L65M6+P6QhAnkpT7ODFv49PMMZkd2OXbRaipqhwj8SYMWOIiYkhMDCQoUOH\n4uPjQ/369UtsP2TIkLJHKYQQf4lLzi41iQC4cD2Zd1efICNbRWaO8co2l1RgycnBqlLLNgtRExmc\nSKSkpBAWFoZKpSIkJOSR7SWREEJUxOlLiQa1i0vOLtN1pWyzEMZlcCLx8ccfExERgY2NDd27d6dB\ngwbyxyaEqDRZ5exhcLSzLLbA0oPvUrZZCOMy+C/q8OHDuLi4sGPHDlxdZalfIUTlcnIwbN2eTi1c\n6NfVU8o2C2EiBv/F5eXl4efnJ0mEEKJK9PR1f2QbBTCuX0vaejvh1sBWkgghTMDgvzpfX19iYmIq\nMRQhhHjIzcmWvp09Sm3zdBdPXOvbVFFEQojiGJxIzJkzh+joaJYtW4ZWa5y51kIIUZrx/VoyoFsT\nzP9r8ShzMwUDuzdh7NMtTRSZEOIBg+tIBAcHs2/fPo4cOYKLiwu+vr7Ur18fpbLoMAuFQsEHH3xg\n9GCrC6kjIUTVSsvM49SlRDKy83GwtaRbGzcc69WM9TeEqC0qXCLbx8cHhUKBIc0VCgUXL14sW4Q1\niCQSQggh6poKF6T6+OOPjRaMEEIIIWoHg3skxEPSIyGEEKKuqXCPRGVLSkoiKCiI0NBQkpOTcXR0\nxN/fn7lz59KkSZNCbXfu3Mm6deuIiYnBwcGBZ599lldffZV69eoZdK+cnBxWrVrFTz/9REJCAp6e\nnkyYMIHx48dLkS0hhBCiDErskVi4cCEKhYK5c+fi7OzMwoULDb9oGQdbJiUlMXr0aOLi4v6/vXsP\ni6ra/wf+nmEYhrsKiAyICDaQYorcVOiYllpqppaXUiP5lVhq1qNmpubx5DG/p5OPqZmleQzT8p56\nrJOlqYE64g28oYmYOiBXkQEGZobZvz+QSZwBZASGy/v1PD7KXmuvvXbu4O1nr70HUVFRCAwMRHp6\nOg4dOgRXV1ds2bIFfn5+AIAvv/wSy5YtQ2BgIP72t7/hypUrOHz4MEJCQhAfHw+ptOYFWOXl5Xjz\nzTdx+PBh9OvXDwqFAkeOHMHly5cRGxuLOXPmPMR8WZEgIqLWpc6LLSsXV/7444/o3LkzgoKCHvpg\ndV1s+eGHH2LLli14//33MWnSJOP23bt347333kP//v2xZs0aqFQqDBo0CN27d8fGjRtha2sLAPjs\ns8+wevVqLFiwABMmTKjxWHv37sWsWbOqhAadTofXX38dSqUSu3fvRmBgYI1jMEgQEVFrU+cgsWvX\nLgDAwIED4eTkZPz6YY0cOfKh+/bt2xeCICAxMRHiBz5Vb+DAgcjMzERKSgo+++wzrFmzBmvWrEH/\n/v2NfcrKytC3b1/4+Phg9+7dNR5r3LhxOHfuHI4ePQpXV1fj9jNnzmDcuHF49dVXMW/evBrHYJAg\nIqLWps5rJB4MAnUJBnVRXl6OuLg4SCQSkxABAFKpFDqdDnq9HklJSQCAiIiIKn3s7OzQs2dPJCQk\nQK1Ww9nZ/MlqtVqcO3cOQUFBVUIEADzxxBOwt7c3HoOIiIhq99Bvtpw7dy62bdtWa7+vvvqqyu2J\n2tjY2CAmJgbjx483aUtLS8O1a9fg6+sLqVSKGzduwN3d3eyiSm/vilfppqenV3sslUoFvV4PX19f\ns/Po0KEDXwNORERUBw8dJHbt2vVQ/1pXKpU4derUI00KAAwGAz766CMYDAaMGTMGAFBQUFBttaFy\ne1FRUbVjFhQUVOlrbgyNRgO9Xv8oUyciImo9hGrMnDlTmDhxovFXYGCgEBUVVWXbg79GjhwpBAUF\nCQMGDKhu2IdiMBiEefPmCQqFQhg1apRQVlYmCIIgBAYGCs8//7zZfZYvXy4oFArh4MGD1Y574sQJ\nQaFQCB999JHZ9ldeeUVQKBRCUVFRjfNbsGCBAMDkl0qlEgRBEBYuXMh2trOd7Wxne4trN6faxZb7\n9u3DzJkzjV8/7OuxJRIJFi9ejBEjRtTa1xy9Xo8FCxZg586d6NixIzZt2gRPT08AQM+ePeHl5YWf\nfvrJZL9PPvkE69atwzfffIPevXubHTslJQWjR4/GuHHjsGjRIpP20aNH49y5c7h48aLZ9RqVuNiS\niKj5Ss45j1/+PIT0whsAAFuxLcI9QzCk8zNoK2tj5dk1XXVebDl06FDI5XIYDAYIgoAJEyYgOjoa\nb775ptn+IpEIdnZ28PHxMVnI+LA0Gg1mzJiBw4cPw8/PD//5z3+MIQIAXFxcoFab/yFeub262xYA\njPOq7vaHWq2Gg4NDjSGCiIiar19vHMauq/uqbNMZdDiaeQLn8y7h3V5vor2Du5Vm1zzV+GbLkJAQ\n459HjhyJXr16ITQ0tEEmcvfuXbzxxhtITk5G165dsW7dOri5uVXp4+fnh6SkJJSWlkImk1VpU6lU\nEIvF6NSpU7XH8Pb2hq2tLW7dumXSVl5ejtu3byMgIKB+ToiIiJoUVVGmSYi4X6FWjU2p2/BuL/P/\nYCbzHvqf3h9//DFGjx7dIJMoKytDXFwckpOTERERgY0bN5qECAAIDQ2FwWDAyZMnTfY/e/YsunTp\nAicnp2qPI5FI0KNHD1y8eNGkKpGSkgKNRlMlPBERUctxRHWs1j5XC9KhKspshNm0HE2ihr9s2TKc\nOXMGISEhWLt2bbVhYNiwYbCxscGqVaug1WqN29esWYOioiKMHTu21mONGDECWq0WK1euNG7T6XT4\n7LPPAKDBwhIREVnX9YI/H67f3RsNPJOWxeof2pWTk4NNmzYBAPz9/bF27Vqz/SZPnoyAgADExsZi\n7dq1GDFiBPr374+rV6/i0KFD6NWrl/Ex0UobNmyAWq1GTEwMXFxcAACjRo3Cjh07sGHDBly5cgXd\nunXD77//jtTUVMTGxtb6emwiImpeNHoNkm6fRZYm5+F24Gc31onVP0b8119/xdSpU2vtl5SUBBcX\nFwiCgM2bN2Pz5s24ceMGPDw8MHDgQEybNs1koeWAAQOgUqlw4MAB+Pj4GLcXFRVh5cqV+Omnn1BQ\nUABfX1+8/PLLePnllx9qoSWf2iAiatoEQcCf6ptIVClxMusstAbdQ++7IHImOjh61t6xlanzZ21Q\n9RgkiIiaporqwxkkZCgtWusQ2LYL3g6Z3AAza/7q/PgnERFRcyAIAq4X3kBChhKnspKhM1N9kNnI\nENEhBHY2dvjlxiGz47S1a4MJj3OdXF0xSBARUbNUotPgRNZpJKqUyCi+bbaPn4svouSRCPXsATsb\nKYCKqsOvNw4j9c4fAAB7iQy9O4RhkF9/uEirfxcRmcdbGxbgrQ0iIusQBAHphX8iQaXE6eyUGqoP\nvRAlj4CPs7zasUr1pSgr18HJ1gE2YpuGnHaLwFsbRETUbJXoSnDi9hkkZlRffejs0glR3pHo1f4J\nY/WhJjKJDDKJrNZ+VDMGCSIiapIEQcC1u38iMUOJ09nJ0BlMP5nZXiJDRIdQRMkj4O3kZYVZEoME\nERE1KcW6Epy4fRoJGUrcLs4y28fftROi5b0R0r47pA9RfaCGwyBBRERWJwgC0u5eR4JKiTM5KdCb\nrT7YI7JDL0TJIyF36mCFWZI5DBJERGQ1RbpinMg8hYSME8gqyTbbJ8DVD1HySIS0fwJSG9tGniHV\nhkGCiIgalSAIuFqQjsQMJc7knDNbfXCQ2CPSKxRR8kh48S2TTRqDBBERNYoibTGUt08hMUOJrBLz\nn3sR4NoZ0d6RCPHoDltWH5oFBgkiImowgiDgj4JrSMxQ4mz2OeiFcpM+jhKHe9WHCH7GRTPEIEFE\nRPVOrS0yVh+yS3LN9nmsjT+i5JHo6RHM6kMzxiBBRET1QhAEXLmTVlF9yDmPcnPVB1sH9O4Qhr7y\nCHRwbG+FWVJ9Y5AgIqJHotYW4XjmSSRmKJGjyTPb57E2/oj27o0eHsGwFfNHT0vCv00iIqozg2Aw\nVh+Scy6YrT442Toan7zwdPCwwiypMTBIEBHRQyvUqu9VH04gt5rqg6JtF0TLI/AEqw+tAv+GiYio\nRpXVh4QMJVJqqD708QpHX3k42rP60KowSBARkVl3y9RQ3lv7kFuab7ZPUNvHEOUdiSfcu0LC6kOr\nxL91IiIyMggGXM6/WlF9yL0Ag2Aw6eNs64TeXmGIkkfCw8HNCrOkpoRBgoiIcLesEMcyT+Joxgnk\nsfpAdcArgYiolTIIBlzK/wOJGUqcy71ovvogdUIfr3BEySPgbs/qA5likCAiamUKyu4an7zIL71j\n0i6CCEHtHkO0PBLd3bvCRmxjhVlSc8EgQUTUClRUH64gUaXEubxLZqsPLlJn9PUKRx95BNzt21lh\nltQcMUgQEbVgBWV3cSwjCYkZJ3CnrMCkXQQRHm+nQLR3JILdHmf1geqMQYKIqIUxCAZczLuMhAwl\nzudeggDBpI+r1Bl95BHo6xUON1Yf6BEwSBARtRB3SgtwNDMJxzKSqq0+dHULRJQ8EsFuQaw+UL1g\nkCAiasbKDeW4mH8ZiRlKnM9NNVt9aGPnij5e4ejjFQ43+7ZWmCW1ZAwSRETN0J3SAhzNOIGjmUko\nKLtr0i6CCN3cghDtHYmu7QJZfaAGwyBBRNRMlBvKcSEvFYkZSlzIu1xt9aHiyYtwtJOx+kANr8kF\niaysLAwZMgTTp0/Ha6+9VqWtpKQEX3zxBX788UdkZ2dDLpdj+PDhiI2Nhb29/UONP3v2bOzZs8ds\n2xtvvIFZs2Y96ikQEdWrPM0dHMs8gaMZSbirLTRpF0GEYPcgRMlZfaDG16SCRHFxMaZPn46ioiKT\nNo1Gg4kTJ+L8+fMICAjAuHHjoFKpsHLlShw5cgRff/01nJycaj1Gamoq3N3dMW7cOJO20NDQejkP\nIqJHVW4ox/m8VCRkHMelvCtmqw9t7dqgr7xi7UNbWRsrzJKoCQUJlUqF6dOn48KFC2bb161bh/Pn\nz+Ppp5/G8uXLIZVKAQA//PAD5syZg88//xxz5syp8Rg6nQ7p6el46qmnMH369Ho/ByKiR5Wnyb/3\n5MUJ3NWqTdorqg+PI1oeia5ugRCLxFaYJdFfmkSQ2LBhA1asWIHS0lL07t0bx48fN+mzb98+iEQi\nLFy40BgiAGDEiBGIj4/Hd999hxkzZkAmk1V7nLS0NOh0OgQGBjbIeRARWaLcUI5zeZeQqFLiUn71\n1YcoeQT6yMPRxs7VCrMkMq9JBIn4+Hh4e3tj0aJFuH79utkgcevWLXh5ecHT09OkLSgoCBcuXMD5\n8+cRFhZW7XEuX74MAAwSRNQk5GrycTTjBI5lJqHQTPVBLBKju9vjiPKOxOPtFKw+UJPUJILEokWL\n0LdvX9jY2OD69etm+0ilUuh0OrNtlWsqVCrVQwWJ9PR0jBs3DpcvX4ZMJsNTTz2Fd955x2xIISKq\nT+WGcqTkXkRihhKp+X+YrT60k7VFlDwCvb3CWH2gJq9JBIknn3yy1j7BwcFQKpU4c+YMQkJCjNuL\nioqQlJRk/HNNKoPE6tWrMXDgQPTs2RPJycnYuXMnEhMTsXXrVnTo0OERzoSIyLxcTR4S71Uf1FrT\n71VikRhPuHdFlDwSQe0eY/WBmo0mESQexqRJk6BUKvHuu+/i73//O8LDw6FSqbBkyRKUlpYCAATB\nNNnfTyaTwc/PD6tWrcJjjz1m3P7FF19g+fLlWLx4MVatWlXrXNq2dYBEwseriKhm+nI9Tmak4Ne0\nBKRkXTLbx8PRDU/7R6F/575oa8/qAzU/zSZI9O/fH3PmzMGnn36KuLg44/awsDBMmTIFy5Ytq3Gh\nJQB8/vnnZrfHxcVh+/bt+O2331BcXAxHR8cax7lzp6TuJ9DKZJXkIFeTD5mNHfxcOvK5dmpVskty\ncTTjBI5nnoRaV131oRui5ZEIbNcFYpEY+iIgp8h0nQRRU+Hh4Wx2e7MJEgAQGxuLgQMH4siRIygp\nKUG3bt3Qp08fY0Dw8PCwaFyxWIygoCDcunULt2/fRkBAQH1Ou1X5s/AmdvyxF2l3rxu3uUid8bTv\n3/B0x79BJBJZb3JEDUhv0CMl9yISVMdx+c5Vs33cZe0QJY9EpFcYXO3Mf1Mmam6aVZAAgI4dO2L8\n+PFVtp0/fx4A0KVLl2r3Ky0txZUrVyCVShEUFGS2HQDs7Ozqcbaty7W717HizFroDFUXxRZq1dh1\ndR/yS+9gjGKElWZH1DCyS3JwNCMJxzKTUKQrNmkXi8To4RGMaHkkFG0DuPaBWpxmEyQ++eQTbN26\nFT/99BPc3d2N2wsLC3H8+HEEBATA29u72v2zs7MxevRoKBQK7N27t0qbRqPBxYsX0a5duxrHoOoJ\ngoBNl7abhIj7Hb51FGGeIfB37dSIMyOqfzqDHik555GQcQJXqqs+2LshWh6JSK9QuEhZfaCWq9kE\niYCAABQWFuL777/HtGnTAAAGgwEff/wxNBoNJk+eXOP+vr6+CAoKQmpqKvbs2YPhw4cDqPgB+Omn\nnyI/Px9Tp05l6d1CaXev43ZJdq399qT9hOEBz8LZ1hnOUifIJKwAUfORVZKDxAwllJmnzFYfbEQ2\n6OHRDVGsPlAr0myCxPDhw7F582asWrUKqamp6NSpE44dO4YLFy7g+eefNwaDShs2bIBarUZMTAxc\nXFwAVLyvIiYmBu+99x72798Pb29vnDx5EufPn0d4eDimTJlijVNrEVRFmQ/V74+Ca/j01Grj11Kx\nLZylznCROsNF6gRnqdO9r53gInWGs7QicLhInWBnY8egR41OZ9AjOfscEjKU+KPgmtk+HvZuiJJH\nordXGJyltX/mD1FL0myChEQiwfr16/HZZ5/ht99+Q0JCAjp16oSFCxdi7NixEIurJv/4+HioVCqM\nHDnSGCR69uyJbdu2YdWqVVAqlTh06BB8fHwwY8YMvP7661VevU11Yyu27FLSGnTIK81HXmn+QxzD\n9r7A8VfAeDBwOEudIWPooEeUVZyNhAwllLdPoVhn+qSWjcgGPT2CEe0dicfaBPB6o1ZLJNT28gUy\nkZPDR7QelKfJx8Jj/2f2LX33c5A4wCAYUFpe2qDzqQgdNQWOv35n6KBKunIdzuacR2IN1Yf2Du4V\nT150CGX1gVqVFvH4JzVdbvbt0NMjGGdyzlXbRyq2xbzId9HGzhXach3U2iIUatVQa9X3/lwEtU6N\nwjK18c9qbRE0+rqHDp1Bh7zSO8grvVNrX1uxpErAqBpAnO8LHU6Q2cgYOlqg28VZSMw4AWXmKRTr\nTasPEpENerbvjmh5JLq08ec1QHQfBgmqN68EvYi80nzcUKtM2mzFEvy/4AnGzw2Q2tjCzb4t3Ozb\n1jqurlxXJVgUav/6vVBbVCWIaPSaOs9bZ9Ajv/QO8h8idEjEEjjbOsHF7l7gsP0reLjYOVe03fva\nXsLQ0ZRpy3U4m3MOCSol0u6mm+3j6eBhrD44SWt+UR1Ra8VbGxbgrY3qact1OJp5AscykpCryYOd\njR16eHTDUx2j4elg2QvD6kJXroNaV1Rr4CjUqi0KHXVhDB0PrN8w/doJ9hJ7ho5GklmchURVxdqH\nEjPXgEQsQYhHd0TJI9GlTWf+vRDdU92tDQYJCzBItAw6gx5FVQLHg7da/vrd3A+c+iQR2ZgsGDUX\nOFykzgwdFtCW63AmOwUJGUpcu++tq/fzdGiPaO9IRHToBSdbVh+IHsQgUY8YJFofvUEPtbaomgpH\n1d/N3WOvTxKRDZzMBg7ThaUOrTx0ZBTdrnjvw+3TZitQFdWHJxDtHYkAV79W/d+KqDYMEvWIQYJq\nojfoUaQrrggcZeq/AkjlGo8yNQp1FUHE3GOF9clGZFOlqvHgEyv3b3eUOLSIH6Taci1OZ6cgMUOJ\na3f/NNung6MnouUV1QdHW4dGniFR88QgUY8YJKi+lBvKodZVvb2irmZth7k3KdYnsUhsDBnOUie4\n3Hv7qIudM1xsqwYRB1v7Rntro96gx+nsFCRlnYFaWwQXqTMiPEPQs313SO57f4mqKBOJGUqcuH3a\n7JM+tmIJerXvgSh5JPxdO7WI0ETUmBgk6hGDBFlDZegwHzjUVW69FOtKan2nx6MQi8RVnlCp6dHZ\nRwkdd8vU+Dx5ndk3p3Z0kuP17hPxx51rSMxQIr3whtkxvBw9ES3vjYgOIXBg9YHIYgwS9YhBgpq6\nckP5vdsrNQeOxgsdjmbfzfHg1462DsbQIQgCPj31ebUBAQDEEMFgZu62YluEtu+BKO9IdHbxZfWB\nqB4wSNQjBglqSSpCR4nZwFEZRArvu73S0KHDydYRzlInSEQ2+FN9q077yx07IMo7EhGeveBga99A\nsyRqnRgk6hGDBLVWBsGAIl3xX4tGteqqazzK/vq6SNuwoaOSWCRGRIdeiJZHwo/VB6IGwyBRjxgk\niGpnEAwo1pU8UOEw/7taW2Rx6IjwDEFMt5frefZE9CB+1gYRNSqxSGxcB+ENrxr73h86KgNG0u3T\nuJh/pdbjdHD0rK8pE5EFGCSIyOruDx2VurTpjA+PLq2xUiEWiRHpFdoYUySiajTOg+BERHXUTtYW\ngzv1r7HP4E4DjB8ER0TWwYoEETVZw/wHQ2ojxf4/f0NpeZlxu8xGhsF+/THQ9ynrTY6IAHCxpUW4\n2JKocZXqy3A+9yLUumI4S50Q7PY4ZBI7a0+LqFXhUxv1iEGCiIham+qCBNdIEBERkcUYJIiIiMhi\nDBJERERkMQYJIiIishiDBBEREVmMQYKIiIgsxiBBREREFmOQICIiIosxSBAREZHFGCSIiIjIYgwS\nREREZDF+1gYRERFZjBUJIiIishiDBBEREVmMQYKIiIgsxiBBREREFmOQICIiIosxSBAREZHFJNae\nADUtOTk5WLlyJQ4fPoy8vDy4urqiT58+mDFjBjp27Fil7w8//IANGzbg+vXrcHFxwXPPPYe3334b\njo6O1Y5vMBgwduxYeHh4YPXq1SbtgiBg9+7diI+PR3p6Ouzt7REdHY23334bPj4+9X6+1PCsfU3d\nuHEDAwcOrHb/lJQU2NnZWX6C1OiseU29//772LVrV43zGzlyJJYuXWr5CTYzDBJklJOTg9GjRyMz\nMxNRUVEYMmQI0tPT8d///he///47tmzZAj8/PwDAl19+iWXLliEwMBATJkzAlStXsGHDBiQnJyM+\nPh5SqdTsMRYvXoyUlBQ8/fTTZtuXL1+ONWvWwN/fH2PGjEF+fj5+/PFH/Pbbb9i6dSs6d+7cUKdP\nDaApXFOpqakAgCFDhsDf39+k3cbGpn5OlhqFta+pZ555Bt7e3mb327JlC3JychAeHl5v59ssCET3\nLFiwQFAoFML69eurbP/hhx8EhUIhxMXFCYIgCLdu3RK6du0qjB07VtBqtcZ+y5cvFxQKhbBx40aT\nsTUajTBr1ixBoVAICoVCePPNN036pKWlCQqFQhg9enSVcY8cOSIoFAph+vTp9XWq1EisfU0JgiCs\nWLFCUCgUwqVLl+rxzMhamsI1Zc7+/fsFhUIhzJo1y8Iza764RoKMfv31V7Rr1w4xMTFVtr/wwgvw\n9fVFQkICDAYDtm7dCr1ej7i4ONja2hr7TZkyBU5OTti2bVuV/Y8ePYqhQ4diz549iI6Orvb4qamp\n8PLyQmxsbJVxn3zySbi6uuLs2bP1dKbUWKx9TQHA5cuXYWtri4CAgPo7MbKapnBNPaiwsBALFiyA\nm5sb5s+fb/nJNVO8tUEAgPLycsTFxUEikUAsNs2XUqkUOp0Oer0eSUlJAICIiIgqfezs7NCzZ08k\nJCRArVbD2dkZALBnzx4UFxdjyZIliIyMrLYEPWTIEAwZMsRke25uLgoLC7lGoplpCtcUUBEkOnfu\nXOWHCTVPTeWaetCaNWtw584dLF68GK6uro9whs0TgwQBqLhP/GDCr5SWloZr167B19cXUqkUN27c\ngLu7u9nFSpX3DtPT0/HEE08AAF566SXMnz8fTk5OuHXr1kPPSaPRICUlxbhoafLkyXU9LbKipnBN\nlZSU4ObNm4iIiMCiRYtw+PBh5ObmIiAgAJMmTcLw4cPr4UypsTSFa+pBWVlZ+Pbbb+Hv749Ro0ZZ\ncFbNH4ME1chgMOCjjz6CwWDAmDFjAAAFBQXVVgcq031RUZFxW1hYWJ2P++BK+7lz5+LZZ5+t8zjU\n9DTmNXXlyhUIggClUomCggIMHjwYd+7cwcGDBzF79mxcv34db7/99iOeEVmbtb5PAcDmzZtRVlaG\n2NjYVrtwl0GCqiUIAj788EMcO3YMwcHBxn8J6PX6alc7V24vKyt7pGPr9Xq89tpr0Gg0OHjwIJYu\nXYri4mJMnTr1kcYl62rsa0qtVqNz586IiorCvHnzjOXwrKwsvPzyy1i9ejUGDRqEoKAgC8+IFT/o\n7QAACRlJREFUrM2a36e0Wi22bt0KDw8PvPDCC480VnPGxZZkll6vxwcffIBt27ahY8eOWL16tfF/\nPplMBp1OZ3Y/rVYLALC3t3+k4/v7+2Pu3Ln4xz/+gb1798LPzw8rVqxASkrKI41L1mONa+rJJ5/E\n//73PyxYsKDKPXVPT09MnToVgiBg3759FpwNNQXW/j515MgR5OfnY+TIkdWGltaAQYJMaDQavPXW\nW9i5cyf8/PwQHx8PT09PY7uLiwvUarXZfSu3V5YO60Pbtm3x1ltvAQAOHDhQb+NS42lq1xQAdOvW\nDQDqdD+cmo6mcE0dPHgQAFr9bVcGCari7t27iImJweHDh9G1a1ds3rwZcrm8Sh8/Pz/k5eWhtLTU\nZH+VSgWxWIxOnTrV+dh//PEH9u7da7bcWDmHO3fu1Hlcsi5rXlM3b97E8ePHodFoTNoqj8W3WjY/\n1rym7nfkyBF07NjRGEpbKwYJMiorK0NcXBySk5MRERGBjRs3ws3NzaRfaGgoDAYDTp48abL/2bNn\n0aVLFzg5OdX5+N988w1mzZqFxMREk7bLly8DAHx9fes8LlmPta+pFStWICYmBr///rtJ26lTpwAA\nwcHBdR6XrMfa11QllUqFnJwc9OrVy+IxWgoGCTJatmwZzpw5g5CQEKxdu7ba/8mGDRsGGxsbrFq1\nynivEah4lrqoqAhjx4616PjPPfccgIpv/vf/K+LmzZtYvXo17OzsMGzYMIvGJuuw9jU1ePBgAMDn\nn3+OkpIS4/Zr167hq6++gqurK6+pZsba11SlS5cuAQC6du36SOO0BHxqgwBUvL9+06ZNACoWOq5d\nu9Zsv8mTJyMgIACxsbFYu3YtRowYgf79++Pq1as4dOgQevXqZXz8qq6ioqIwatQo7Ny5E0OHDsWA\nAQOgVquxf/9+lJaWYunSpejQoYPF50iNqylcU8888wwGDRqE/fv3Y9iwYRgwYAAKCwvxyy+/QKvV\nYuXKlWjTpo3F50iNqylcU5Vu3LgBAGjfvv0jjdMSMEgQACA5Odm4wnnHjh3V9ouJiYGdnR1mzpwJ\nLy8vbN68GfHx8fDw8MBrr72GadOmPdLq5SVLliA4OBjff/89vvvuO8hkMoSFhSEuLg6hoaEWj0uN\nr6lcU8uXL0d8fDx27NiB77//Hvb29oiMjMRbb71lfBkRNQ9N5ZoCKt5TAdT/IuDmSCQIgmDtSRAR\nEVHzxDUSREREZDEGCSIiIrIYgwQRERFZjEGCiIiILMYgQURERBZjkCAiIiKLMUgQERGRxRgkiIiI\nyGIMEkRERGQxBgkiIiKyGIMEERERWYxBgogazPr16xEYGIi5c+eabf/1119N2vfu3YtXXnkFvXr1\nQs+ePfHSSy9h27ZtMPexQBqNBuvWrcOYMWMQHh6O4OBgREVFYdq0aUhJSanSd+fOnQgMDMTGjRvx\nr3/9C6GhoQgNDcXChQvr96SJWhkGCSJqMMOHD4dEIsH+/ftRVlZm0v7DDz8AAEaOHAkAmDdvHmbN\nmoVLly6he/fu6N27N9LT0zF//nzMnj27SpgoLS3F+PHj8cknnyArKwthYWGIioqCSCTCL7/8glde\neQXnzp0zOea3336Lb775BmFhYfDz80Pnzp0b6OyJWgd+jDgRNRh3d3f069cPBw4cwIEDBzBkyBBj\nW0FBAQ4dOoSOHTsiPDwc27Ztw/bt2/H444/jiy++gJeXFwAgPz8fU6ZMwd69exEeHo6xY8cCqAgE\nFy5cwODBg7Fs2TJIJBXfzrRaLWbNmoWff/4ZW7ZsQffu3avM6fr161i9ejWefvppAIDBYGiM/xRE\nLRYrEkTUoF588UUAwJ49e6ps37dvH3Q6HUaMGAGRSISvv/4aALB06VJjiACAdu3a4Z///CeAilsl\nlWQyGfr164eZM2caQwQASKVSjBo1CgCQmZlpMh9vb29jiAAAsZjfBokeBSsSRNSg+vXrB3d3d/z+\n++/Iz89Hu3btAFTc1hCJRBgxYgSys7ORnp6ONm3aICgoyGSMxx57DJ6enrh+/TpycnLg4eGBCRMm\nYMKECVX6FRYW4sqVKzhy5AiAiurEg8yNT0SWY5AgogYlkUgwfPhwrF+/Hvv27cPEiRORlpaGlJQU\nREREwMfHx7gwsqCgAIGBgTWOl5mZCQ8PDwBAdnY2Nm3aBKVSifT0dBQUFAAARCJRtfu7urrW05kR\nEcAgQUSN4KWXXsL69euxd+9eTJw40Xibo/IWRHl5OQDAzc0Nffv2rXEsR0dHAMDx48cxZcoUaDQa\nyOVyREREwN/fH8HBwZBIJJgyZYrZ/Xkrg6h+MUgQUYMLCAhAjx49kJKSgtu3b2P//v1wcHDA4MGD\nAcBYYXB0dMS///3vWscTBAHz58+HRqPB4sWLMXr06CrtBw8erP+TICKzGM2JqFG8+OKLEAQB69at\nw7Vr1/Dss8/CwcEBAODj4wMvLy/cunULaWlpJvsWFBRg+PDhiI2NRXFxMfLz83Hz5k14eHiYhAgA\nSExMBMAnMogaA4MEETWKoUOHQiaTYdOmTQD+endEpZiYGBgMBsyePRsqlcq4vbS0FB988AEuX74M\ne3t7ODo6ok2bNpDJZMjNzUVycrKxryAI2LFjB7777jsAMPvuCiKqX7y1QUSNwsnJCYMGDcKePXuM\n7464X0xMDM6cOYOff/4ZQ4cORffu3eHk5ISzZ88iPz8ffn5+WLRoEQDAxsYGEydOxNq1azF+/HhE\nRETAwcEBly5dwq1bt/Dqq69i48aNyM3NtcapErUqDBJE1GhCQkKwZ88e47sj7icWi7F8+XLs2rUL\n27dvx4ULFyAIAnx8fDB27FhMmjSpyhMX77zzDjw8PLB9+3acPn0aUqkUISEhWLJkCSIjI3H06FFc\nvXoVqampfOSTqAGJBHMvsCciagCvvvoqkpKScODAAcjlcmtPh4jqAddIEFGDKi0tBVDxoVlKpRJP\nPfUUQwRRC8JbG0TUoF577TVcvHgRZWVlsLW1xYwZM6w9JSKqR6xIEFGDCg4OhiAI8Pf3x6pVq7he\ngaiF4RoJIiIishgrEkRERGQxBgkiIiKyGIMEERERWYxBgoiIiCzGIEFEREQWY5AgIiIii/1/avLc\nL9k/0MgAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "results_melt = pandas.melt(results, id_vars=['year'], value_vars=['Hugh', 'Jesse'], \n", " var_name='runner', value_name='time')\n", "\n", "seaborn.pointplot(x='year', y='time', data=results_melt, hue='runner', fit_reg=False)\n", "plt.ylabel('time (minutes)')\n", "plt.axhline(y=20, color='black', linestyle='--', linewidth=1)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Examine if there is a significant difference between runners\n", "We perform statistical tests to determine if the runners have significantly different times.\n", "\n", "First, we plot the distributions for the two runners." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhIAAAFrCAYAAACE+GArAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlcVPX+P/DXDDAsIiCCqCyOkUJKKqAYkgtu93u9mXYT\nxVxwSUnLqJu5/LrqvZllqZSBiElKkBYuqFnXrmkK4oIboobKFTFjBAQUZRlgtt8fPZgblxkYjoMs\n83o+Hj0e8jmfOed9eGi8+JzP+XxEGo1GAyIiIiIBxC1dABEREbVdDBJEREQkGIMEERERCcYgQURE\nRIIxSBAREZFgDBJEREQkmHlLF9AWFRWVtXQJRERET5Szc0ed7RyRICIiIsEYJIiIiEgwBgkiIiIS\njEGCiIiIBGOQICIiIsEYJIiIiEgwBgkiIiISjEGCiIiIBGOQIJMir1aiRqFq6TKIiNoNrmxJ7Z5K\nrcbPF2U4niFDfkklAKC3mz3GDPKAv5dzC1dHRNS2iTQajaali2hruER226FSqxG19wou55ToPD5+\niBQvDXvqCVdFRNT2cIlsMklHL8j0hggAOHjqNm7mPXyCFRERtS8MEtRuaTQaHLuY12i/nw3oQ0RE\nujFIULtVLleg8IG80X45dzkiQUQkFIMEtVtmYpFB/cQiw/oREVF9DBLUbtlYWUDaVffkoD/q09Px\nCVRDRNQ+MUhQuzZ6oFuDx83EIoz0a7gPERHpxyBB7doQn256w4SZWIRZf/aGq1OHJ1wVEVH7wXUk\nBOA6Em3PL7n38fPFPNwuKINYJILPU44Y5ecGty62LV0aEVGboG8dCQYJARgkiIjI1HBBKiIiIjI6\nBgkiIiISjEGCiIiIBGOQICIiIsEYJIiIiEgwBgkiIiISjEGCiIiIBDNv6QJqFRUVISoqCikpKSgp\nKYG9vT0CAwMREREBd3d3bb/y8nLExMTgp59+Qn5+Pjp06AB/f38sWrQIzzzzjEHXksvl2LJlC374\n4QcUFhbCzc0N06ZNwyuvvAIRN3AiIiIyWKtYkKqoqAghISHIz89HUFAQvLy8kJubi+PHj8Pe3h5J\nSUmQSqWQy+UIDQ3F9evX4evrC19fXxQUFODw4cMwMzPD9u3b4e/v3+C1VCoVFixYgJSUFAwfPhy9\ne/dGamoqbty4gTlz5mDp0qUG1MsFqYiIyLS06pUtV65ciaSkJCxbtgyzZ8/Wth84cABLlixBcHAw\nYmNj8cUXX2DDhg2YMWMG/v73v2v7nT17FrNmzYKnpycOHjzY4LUOHjyIxYsX1wkNCoUCr776KtLT\n03HgwAF4eXk1eA4GCSIiMjWtemXLI0eOwNHREWFhYXXaJ0yYAA8PD6SlpUGtVuPw4cMQiUR46623\n6vQLCAhAQEAAsrOzUVhY2OC1duzYAXNzc7z22mvaNgsLC7z11lvQaDTYs2eP8W6MiIionWvxORIq\nlQrh4eEwNzeHWFw/10gkEigUCiiVSoSGhqKkpAS2tvU3WpJIJACAiooKvdeqqanBlStX4O3tDXt7\n+zrH+vXrB2tra5w7d+4x74iIiMh0tHiQMDMzqzcSUSsnJwe3bt2Ch4cHJBIJJk2apLPf/fv3cf78\nedjY2MDNTfeW0QAgk8mgVCrh4eGhs46uXbvi9u3bgu6DiIjIFLWKRxu6qNVqrF69Gmq1GpMnT26w\n77p161BRUYEJEyZoRyZ0KS0tBQB07Kj7OU/Hjh0hl8uhVCqFF05ERGRCWnxEQheNRoOVK1fi9OnT\n8PHx0TtiAQAxMTFITk6Gq6sr3n777QbPWxsQ9IWN2vbq6mqYm+v/1nTqZANzc7PGboOIiKjda3VB\nQqlUYsWKFUhOToa7uztiYmL0/uDfuHEjYmJi4ODggC1bttSb9/C/LC0tAfz+loYuNTU1EIlEsLa2\nbvA8Dx5UGnAnRERE7Ye+tzZaVZCQy+WIiIhASkoKpFIptm/fDhcXl3r9VCoVVq5ciT179qBz587Y\ntm0bevXq1ej5a4NGeXm5zuNlZWWwsbHROemTiIiI6ms1QeLhw4eYN28eMjMz0adPH8TFxaFz5871\n+tXU1ODNN9/EsWPH4Orqim3btkEqlRp0DVdXV1hYWCAvL6/eMZVKhYKCAnh6ej7urRAREZmMVvGr\nd3V1NcLDw5GZmYmAgAAkJibqDBEajQbvvPMOjh07hl69euGbb74xOEQAgLm5Ofr374+srKx6oxKX\nL1+GXC6Hr6/v494OERGRyWgVQSIyMhIZGRnw9fXF1q1bda4TAQCJiYk4fPgwevTogYSEBJ2PPRoz\nceJE1NTUICoqStumUCiwceNGAEBISIiwmyAiIjJBLb5EdlFREYKDg6FQKPDyyy+jW7duOvvNmjUL\nY8aMwYMHDzB27Fj07t1bZ7/Q0FA4OzsDAOLj41FWVoawsDDY2dkB+P0RxrRp05CRkYEhQ4agb9++\nOHHiBK5fv869NoiIiPRotXttHDlyBK+//nqj/fbv34+JEyca1K92F9CRI0dCJpPh6NGjdRaqKi8v\nR1RUFA4dOoTS0lJ4eHhg6tSpmDp1qkETLRkkiIjI1LTaINEWMUgQEZGpadWbdhEREVHbxCBBRERE\ngjFIEBERkWAMEkRERCQYgwQREREJxiBBREREgjFIEBERkWAMEkRERCQYgwQREREJxiBBREREgjFI\nEBERkWAMEkRERCQYgwQREREJxiBBREREgjFIEBERkWAMEkRERCQYgwQREREJxiBBREREgjFIEBER\nkWAMEkRERCQYgwQREREJxiBBREREgjFIEBERkWAMEkRERCQYgwQREREJxiBBREREgjFIEBERkWAM\nEkRERCQYgwQREREJZt7SBdQqKipCVFQUUlJSUFJSAnt7ewQGBiIiIgLu7u46P1NZWYkXXngBo0aN\nwnvvvWfwtd5991189913Oo/NmzcPixcvFnQPREREpqZVBImioiKEhIQgPz8fQUFBGDduHHJzc/H9\n99/jxIkTSEpKglQqrfMZpVKJxYsXQyaTNfl6169fh5OTE0JDQ+sd8/f3F3obREREJqdVBImoqCjk\n5+dj2bJlmD17trb9wIEDWLJkCdauXYvY2Fhte2lpKf72t7/h5MmTTb6WQqFAbm4uRowYgUWLFhml\nfiIiIlPVKuZIHDlyBI6OjggLC6vTPmHCBHh4eCAtLQ1qtRoA8P3332PcuHE4efIkgoKCmnytnJwc\nKBQKeHl5GaV2IiIiU9biIxIqlQrh4eEwNzeHWFw/10gkEigUCiiVSkgkEiQlJcHKygqxsbGwsbFp\n8qjEjRs3AIBBgoiIyAhaPEiYmZnVG4molZOTg1u3bsHDwwMSiQQAsHDhQvj5+cHS0hLp6elNvl5t\nkMjNzUVoaChu3LgBKysrjBgxAm+99RZcXFyE3wwREZGJaRWPNnRRq9VYvXo11Go1Jk+erG0PDAyE\npaWl4PPWBomYmBi4ublhypQpkEqlSE5ORkhICAoKCh67diIiIlPR4iMSumg0GqxcuRKnT5+Gj4+P\n3hELIaysrCCVShEdHY1evXpp2zdv3ozPPvsMH3zwAaKjoxs8R6dONjA3NzNaTURERG1Vk4KESqVC\nRkYGzp49i4KCApSUlEClUsHe3h5SqRS+vr7w9/eHhYWF4IKUSiVWrFiB5ORkuLu7IyYmRvtYwxg2\nbdqksz08PBx79uzBsWPHUFFRgQ4dOug9x4MHlUarh4iIqC1wdu6os92gIPHrr79i+/btOHjwICor\nK6HRaHT2E4lEsLOzw4QJEzB37twmzzeQy+WIiIhASkoKpFIptm/f/sTmLIjFYnh7eyMvLw8FBQXw\n9PR8ItclIiJqyxoMEmVlZVi/fj327t0LpVKJLl26YMSIEfDy8oKbmxs6duwItVqNBw8e4N69e8jI\nyMDFixeRkJCApKQkhISEICIiAh076k4xf/Tw4UPMmzcPmZmZ6NOnD+Li4tC5c2ej3SgAVFVVITs7\nGxKJBN7e3jqPA3isORhERESmRG+QOHHiBJYvX46KigqEhIRg0qRJ6Nu3r0EnPX36NJKSkpCUlIQj\nR47go48+QmBgoN7+1dXVCA8PR2ZmJgICArB582bY2to2/W4ace/ePYSEhKB37944ePBgnWNyuRxZ\nWVlwdHSEq6ur0a9NRETUHul9a2PBggUYNWoUjhw5glWrVhkcIoDf36z47LPPcOTIETz33HOYN29e\ng/0jIyORkZEBX19fbN26tVlCBAB4eHjA29sb2dnZdfba0Gg02LBhA+7fv4+pU6dCJBI1y/WJiIja\nG5FGz4SH69ev6xz+F+LatWt45plndB4rKipCcHAwFAoFXn75ZXTr1k1nv/nz59d75JCeno6ZM2di\n5syZOjftio+PR1lZGcLCwmBnZwcAuHTpEsLCwlBdXY3Ro0fD1dUV58+fx9WrVzFo0CBs27at0cmd\nRUVlhtw2ERFRu9HkyZbGChEA9IYIAMjMzIRCoQAA7N27V2+/sLCwJs9dSEhIgEwmw0svvaQNEgMG\nDMDu3bsRHR2N9PR0HD9+HG5uboiIiMCrr75q1DdEiIiI2ju9IxJNcfnyZdy9exd9+/bVu+V3e8IR\nCSIiMjX6RiSatLLl6dOnMXPmTKSmpmrb3n77bUyZMgVvv/02/u///g+ffvrp41VKREREbYbBQeLS\npUuYN28ezp07h9u3bwMAfvrpJxw6dAhWVlYYO3YsHBwc8MUXX+Df//53c9VLRERErYjBQWLbtm3a\nVSenTZsGANi3bx9EIhFWrVqFjRs3Yu/evbCyssK3337bbAUTERFR62HwEtkXL15Ev3798MorrwAA\nampqcPr0aZibm2Ps2LEAgK5du8LPzw9XrlxpnmqJiIioVTF4ROLhw4fo3r279uvz589DLpejX79+\nsLGx0bbb2NigspJ7URAREZkCg4OEi4tLnS22jx07BpFIhKCgIG2bRqPB9evX4ezsbNwqiYiIqFUy\nOEg8++yzyMzMxN69e3Hy5Ens378fALSPNWpqarB+/Xrk5eVh8ODBzVMtERERtSoGryORk5OD6dOn\no7S0FMDvow/jxo1DZGQkAGDo0KEoLi6Gvb09kpKS0KNHj+aruoVxHQkiIjI1+taRaNKCVL/99hs2\nb96M4uJiDBw4ELNnz4aFhQUAYPr06XB0dMSSJUvg5uZmnKpbKQYJIiIyNUYJEg3RaDQms9kVgwQR\nEZmaJu+10RCFQoFr164hPz8fTk5O8Pf3R0FBgd4Nt4iIiKh9alKQUCgUiI6Oxs6dO1FeXg4AGD9+\nPPz9/bFkyRJUVFTg008/bdfzI4iIiOi/DH5rQ6FQYN68efjiiy+gUCjg6+uLPz4VUSqVyMrKwrRp\n03Dv3r1mKZaIiIhaF4ODRGJiIs6cOYPg4GD8/PPP2LlzZ73jU6dORXFxMeLi4oxeKBEREbU+BgeJ\n/fv3o3PnzoiMjISjo2O94+bm5vj73/+Obt26IS0tzahFEhERUetkcJD49ddf4e/vDysrK719zMzM\n4OPjg7t37xqlOCIiImrdDA4SlpaWuH//fqP9ioqKYGlp+VhFERERUdtgcJDw8fHB5cuXkZOTo7dP\ndnY2rl69Ch8fH6MUR0RERK2bwUFi1qxZqKmpwfz583H8+PE6O3xqNBqcOnUKr7/+OlQqlXarcSIi\nImrfmrSy5aZNmxAVFVVnBUsrKyuoVCooFApoNBrMmjULy5Yta5ZiWwuubElERKbGaEtknzp1Cl9+\n+SUuXLiAqqoqAL+/sTFgwACEhYVhzJgxj19tK8cgQUREpsboe22o1WqUlpZCpVLBwcFBu3mXKWCQ\nICIiU6MvSBg8R2L58uXYvXv3fz8oFsPR0RHOzs51QsQXX3yBWbNmCa+UiIiI2gyDg8S+fftw7ty5\nRvulp6fj4sWLj1UUERERtQ16N+1avHhxvT0zTp06hZkzZ+o9WXl5Oa5du4bu3bsbr0IiIiJqtfQG\nieDgYLzzzjvar0UiEYqLi1FcXNzwCc3NsWjRIuNVSERERK1Wg5MtMzIyoFarodFoMH36dDz//PNY\nsGCB7hOJRLC0tISbmxvs7e2breDWgJMtiYjI1OibbKl3RAIAfH19tX9+6aWX4OfnB39/f+NWRkRE\nRG2W4Nc/TRlHJIiIyNQIGpH4o/379zfpghMnTmxS/6KiIkRFRSElJQUlJSWwt7dHYGAgIiIi4O7u\nrvMzlZWVeOGFFzBq1Ci89957Bl9LLpdjy5Yt+OGHH1BYWAg3NzdMmzYNr7zySp1VO4mIiKhhBgeJ\nZcuWGfRDVqPRQCQSNSlIFBUVISQkBPn5+QgKCsK4ceOQm5uL77//HidOnEBSUhKkUmmdzyiVSixe\nvBgymczg6wCASqVCREQEUlJSMHz4cPzpT39Camoq3n//feTl5WHp0qVNOh8REZEpMzhIvPzyyzqD\nhEqlQllZGa5cuYLCwkKMGzcOAQEBTSoiKioK+fn5WLZsGWbPnq1tP3DgAJYsWYK1a9ciNjZW215a\nWoq//e1vOHnyZJOuAwD/+te/kJKSgjlz5mhDQ0REBF599VVs374dEydOhJeXV5PPS0REZIqMNkdC\nrVZjw4YNSEhIwDfffNOkrcSHDBkCjUaDkydPQiyuu0bWmDFjkJ+fj8uXL0MsFuP777/Hhx9+iJKS\nEgQFBeHkyZOYOXOmwY82QkNDceXKFZw6darO2yUZGRkIDQ016FycI0FERKbmsZfIboxYLMbixYvh\n4uKCjRs3Gvw5lUqF8PBwvPHGG/VCBABIJBIoFAoolUoAQFJSEqysrBAbG4vw8PAm1VhTU4MrV67A\n29u73iuq/fr1g7W1tUGrdxIREdHvDH60YQiRSIQ+ffrg1KlTBn/GzMwMYWFhOo/l5OTg1q1b8PDw\ngEQiAQAsXLgQfn5+sLS0RHp6epPqk8lkUCqV8PDw0FlH165dcfv27Sadk4iIyJQZbUSiVk5ODozx\ntEStVmP16tVQq9WYPHmytj0wMBCWlpaCzllaWgoA6NhR9/BMx44dIZfLtaMfRERE1DCDRyQKCwv1\nHlOpVCguLkZiYiJu3bqFIUOGPFZRGo0GK1euxOnTp+Hj46N3xKKpagNC7ejG/6ptr66uhrm5/m9N\np042MDc3M0pNREREbZnBQWL48OGNvv6p0Whgbm6ON954Q3BBSqUSK1asQHJyMtzd3RETE6P3B39T\n1Y5kKBQKncdramogEolgbW3d4HkePKg0Sj1ERERtxWMvSNXQjp5isRg2Njbo3bs3ZsyYgf79+ze9\nQvy+UFTtGg9SqRTbt2+Hi4uLoHPpUjvBsry8XOfxsrIy2NjY6Jz0SURERPUZHCR+/vnn5qwDDx8+\nxLx585CZmYk+ffogLi4OnTt3Nuo1XF1dYWFhgby8vHrHVCoVCgoK4OnpadRrEhERtWet4lfv6upq\nhIeHIzMzEwEBAUhMTDR6iAB+3+K8f//+yMrKqjcqcfnyZcjl8joblREREVHDmvz6Z2VlJW7fvo3K\nysoG384YNGiQweeMjIxERkYGfH19sXXrVlhZWTW1LINNnDgR58+fR1RUFJYvXw7g9zkTtWtfhISE\nNNu1iYiI2huDg4RKpcLatWvx7bffNvp6pEgkQlZWlkHnLSoqwo4dOwAATz31FLZu3aqz3/z585v8\n2md8fDzKysoQFhYGOzs7AMBf//pX7N27F/Hx8cjOzkbfvn1x4sQJXL9+HXPmzOHy2ERERE1gcJDY\nvn07EhMTAQA9e/aEk5NTg69IGiozM1P7FsXevXv19gsLC2tykEhISIBMJsNLL72kDRJmZmaIi4tD\nVFQUDh06hAsXLsDDwwMrV67E1KlThd8IERGRCTJ4r41x48bhzp072Lp1KwIDA5u7rlaNe20QEZGp\neey9NvLy8jB48GCTDxFERET0XwYHCQcHh0YXpCIiIiLTYnCQGDVqFDIyMlBcXNyc9RAREVEbYnCQ\nePvtt9G1a1eEh4fjzJkzqKzkMtFERESmzuDJlpMmTUJFRUWdbbbFYrHexx1Xr141SoGtESdbEhGR\nqXnsvTZ0BQOVSiW8IiIiImrzDB6RoP/iiAQREZmax379k4iIiOh/MUgQERGRYHrnSPj4+AAAvv/+\ne0ilUu3XhmrPky2JiIjod3qDRO3GXLVTKBrbqIuIiIhMDydbCsDJlkREZGo42ZKIiIiMTtA+4KWl\npaiurm6wj4uLi6CCiIiIqO0wOEioVCp8+umn2L17Nx49etRgX5FIhKysrMcujoiIiFo3g4NEXFwc\n4uLiAAC2traws7NrtqKIiIiobTA4SOzbtw8SiQRbtmxBYGBgc9ZEREREbYTBky3z8/MxZMgQhggi\nIiLSMjhIeHh4oKKiojlrISIiojbG4CAxY8YMXLhwARcvXmzOeoiIiKgNMXiOxOTJk3H79m2EhYXh\nxRdfhLe3NxwcHPT2Hz9+vFEKJCIiotbL4CBx//59nD9/HgqFAsnJyY32Z5AgIiJq/wwOEh999BEu\nX74Ma2trBAQEoFOnThCJRM1ZGxEREbVyBgeJEydOwMnJCfv27YOzs3Nz1kRERERthMGTLaurq+Hr\n68sQQURERFoGBwkfHx/cvn27GUshIiKitsbgILFo0SLk5OQgOjoaarW6OWsiIiKiNkKk0Wg0hnRM\nSkrCkSNHkJaWBicnJ/j4+MDBwQHm5vWnWYhEIrz//vtGL7a1KCoqa+kSiIiInihn54462w0OEt7e\n3hCJRDCku0gkwrVr15pWYRvCIEFERKZGX5Bo0uufzamoqAhRUVFISUlBSUkJ7O3tERgYiIiICLi7\nu9fpu3//fsTHx+P27duws7PDn//8Z7z55pvo0KGDQdd699138d133+k8Nm/ePCxevPix74eIiMgU\nGDwi0ZyKiooQEhKC/Px8BAUFwcvLC7m5uTh+/Djs7e2RlJQEqVQKANiyZQsiIyPh5eWFYcOGITs7\nGykpKfD19UVCQgIkEkmj1xs/fjzu37+P0NDQesf8/f0xZMiQRurliAQREZmWJo9I3L9/H46Ojka5\neElJCTp37qz3eFRUFPLz87Fs2TLMnj1b237gwAEsWbIEa9euRWxsLGQyGT7//HP4+voiMTERFhYW\nAICNGzciJiYGu3btwvTp0xusRaFQIDc3FyNGjMCiRYuMcn9ERESmSu9bG3/6058QHx8PpVIp+OQ1\nNTWIjY3F2LFjG+x35MgRODo6IiwsrE77hAkT4OHhgbS0NKjVauzatQtKpRLh4eHaEAEAr732Gmxt\nbbF79+5Ga8rJyYFCoYCXl5ewmyIiIiItvSMSERERWLduHRISEhAWFoaJEyfC3t7eoJPm5+dj9+7d\n2LNnD8rKyrBkyRK9fVUqFcLDw2Fubg6xuH6ukUgkUCgUUCqVOHfuHAAgICCgTh9LS0sMGDAAaWlp\nKCsrQ8eOuodfAODGjRsAwCBBRERkBHqDxPTp0xEUFIR//vOf+Oijj7Bu3ToMHDgQAwcOhLe3N1xd\nXWFrawu1Wo0HDx7g3r17yMjIwPnz55GVlQWVSoWgoCCsXLkSPXr00FuAmZlZvZGIWjk5Obh16xY8\nPDwgkUhw584dODk56ZxU6erqCgDIzc1Fv3799F6vNkjk5uYiNDQUN27cgJWVFUaMGIG33noLLi4u\nej9LREREdTX41kbPnj0RHx+PtLQ0bN26FWfOnMGZM2f0btal0WggFosRFBSEhQsXwtfXV3BharUa\nq1evhlqtxuTJkwEApaWlcHNz09m/dhSivLy8wfPWBomYmBiMGTMGAwYMQGZmJpKTk3Hy5Ens2rUL\nXbt2FVw3ERGRKTHo9c/nn38ezz//PPLz83Hy5EmcO3cOd+/exYMHD6BQKODg4ACpVApfX1+MGDHi\nsX8QazQarFy5EqdPn4aPj492xEKpVOp9K6O2vbq6usFzW1lZQSqVIjo6Gr169dK2b968GZ999hk+\n+OADREdHN3iOTp1sYG5u1pRbIiIiapcMXkcCALp164ZJkyZh0qRJzVUPlEolVqxYgeTkZLi7uyMm\nJkYbEqysrKBQKHR+rqamBgBgbW3d4Pk3bdqksz08PBx79uzBsWPHUFFR0eCaFA8eVBpyK0RERO2G\nvtc/Dd5r40mQy+VYuHAhkpOTIZVKkZCQUGfOgp2dHcrKdK/hUNve0ETLhojFYnh7e0OpVKKgoEDQ\nOYiIiExNqwkSDx8+RFhYGFJSUtCnTx/s3LkT3bt3r9NHKpWipKQEVVVV9T4vk8kgFosbnNhZVVWF\ny5cv4/r163qPA7+/BUJERESNaxVBorq6GuHh4cjMzERAQAASExN1LmDl7+8PtVqN8+fP1/v8pUuX\n8PTTT8PW1lbvde7du4eQkBC8++679Y7J5XJkZWXB0dFR+wYIERERNaxVBInIyEhkZGTA19cXW7du\n1RsGXnjhBZiZmSE6Olo7JwIAYmNjUV5ejilTpjR4HQ8PD3h7eyM7O7vOXhsajQYbNmzA/fv3MXXq\nVL1vpRAREVFdLb7XRlFREYKDg6FQKPDyyy+jW7duOvvNnz8flpaWWL9+PbZu3QpPT08EBwfj5s2b\nOH78OPz8/PDVV1/VeasjPj4eZWVlCAsLg52dHQDg0qVLCAsLQ3V1NUaPHg1XV1ecP38eV69exaBB\ng7Bt27ZG9+vgXhtERGRqHnsb8eZy5MgRvP766432O3fuHOzs7KDRaLBz507s3LkTd+7cgbOzM8aM\nGYM33nij3kTLkSNHQiaT4ejRo3XWn8jOzkZ0dDTS09NRUVEBNzc3vPjii3j11VcN2vSLQaLtKXlY\nhZRMGXLzy2AmFqFvT0cE+XSDjVWTXlwiIjJZrTZItEUMEm3L8UsyfP3vbKj/5696BytzLHq5H3q7\nO7RQZUREbYfRgoRSqcThw4dx9uxZFBQUwM/PD/Pnz8eePXvg4+MDb29voxTcmjFItB2/3L6PyG8v\nQd9fcmtLc3zw6mB06sg3dYiIGtLkbcR1uXr1Kt566y3IZDJoNBqIRCLt3INvv/0WWVlZWLp0qd69\nM4ietB/P/Ko3RACAvFqJYxky/HXYU0+sJiKi9sTgtzby8vIwZ84cyGQyjB07FqtXr8YfBzNGjx4N\niUSCtWtCHrhgAAAfuElEQVTXIj09vVmKJWqKqholsm4/aLRfRnbRE6iGiKh9MjhIbNq0CWVlZVi7\ndi02btyIkJCQOsdfe+01bN68GRqNBtu3bzd6oURNVaNQNzgaUataoWr2WoiI2iuDg0RaWhqeeeYZ\nTJgwQW+fwMBADBgwANeuXTNKcUSPw9baAnYdGn8Lp7uT/n1ViIioYQYHiYa28P4jJycnPHjQ+HAy\nUXMTi0UY2k/3uiR/NHxA90b7EBGRbgYHCWdnZ9y8ebPRfv/5z3/g5OT0WEURGcufB/eAm7P+EYeB\n3l0w4Gn+fSUiEsrgIDF06FDk5uZi586devvs2LEDv/76K55//nmjFEf0uGyszLF0mh+GD+gOSwsz\nbbt9BwkmDu2J117syyXRiYgeg8HrSBQWFmLChAl4+PAhxowZg8GDB2P16tUICgrCxIkTkZKSgh9+\n+AG2trbYt2+fQY9B2iquI9E2yauVuFtcAbFYBPcutjA3axVbzRARtQlGWZAqKysLb775JvLy8iAS\nibRrSQC/b3zl5OSEzz77DAMHDjRO1a0UgwQREZkao61sqVAocPjwYaSnp6OgoAAqlQpdunTBwIED\n8Ze//AVWVlZGKbg1Y5AgIiJTw702jIhBgoiITI1RlsiuVVpaiurq6gb7uLi4CDk1ERERtSEGBwmV\nSoVPP/0Uu3fvxqNHjxrsKxKJkJWV9djFERERUetmcJCIi4tDXFwcAMDW1la7WRcRERGZLoODxL59\n+yCRSLBlyxYEBgY2Z01ERETURhgcJPLz8zFkyBCGCGqTKquUOHU1H7n5ZTATi+DzlCP8ejtzLQki\nosdkcJDw8PBARUVFc9ZC1Cwu3LiHuB+uobrmv7t8pl3JR2c7S7w5qT/cu9i2YHVERG2bwb+OzZgx\nAxcuXMDFixebsx4io7ope4jYA7/UCRG1Sh5VY0PSJZRV1rRAZURE7YPBIxKTJ0/G7du3ERYWhhdf\nfBHe3t5wcHDQ23/8+PFGKZDocRw68ytUav1LpTyqqEFq5l38JVD65IoiImpHDA4S9+/fx/nz56FQ\nKJCcnNxofwYJamnVChUyb5Y02u/stXsMEkREAhkcJD766CNcvnwZ1tbWCAgIQKdOnbhrIrVq1TUq\nqA1YuLWySvkEqiEiap8MDhInTpyAk5MT9u3bB2dn5+asicgoOlibo4OVOSoaCQpdOlk/oYqIiNof\ngydbVldXw9fXlyGC2gwzsRhBz3ZrtN/wAd2fQDVERO2TwUHCx8cHt2/fbsZSiIzvz8/1gJO9/h1p\n+0g7wd+L4ZiISCiDg8SiRYuQk5OD6OhoqNXq5qyJyGjsO0iwbJoffHs5QfyHOT2WFmYY6eeKiEn9\nYCbmolREREIZvI14UlISjhw5grS0NDg5OcHHxwcODg4wN68/zUIkEuH99983erGtBbcRb5vuP6rC\nrwVlEItF6OXmABsrQZvfEhGZJH3biBscJLy9vSESiWBId5FIhGvXrjWtwjaEQYKIiEyNviDRpNc/\niYiIiP7I4BEJ+i+OSBARkal57BGJ5lZUVISoqCikpKSgpKQE9vb2CAwMREREBNzd3ev03b9/P+Lj\n43H79m3Y2dnhz3/+M95880106NDBoGvJ5XJs2bIFP/zwAwoLC+Hm5oZp06bhlVde4SJbRERETaB3\nRGLFihUQiUSIiIhA586dsWLFCsNP2sTJlkVFRQgJCUF+fj6CgoLg5eWF3NxcHD9+HPb29khKSoJU\nKgUAbNmyBZGRkfDy8sKwYcOQnZ2NlJQU+Pr6IiEhARKJpMFrqVQqLFiwACkpKRg+fDh69+6N1NRU\n3LhxA3PmzMHSpUsNqJcjEkREZFqaPNmydnLlv/71L/Ts2RPe3t4GX6ypky1XrlyJpKQkLFu2DLNn\nz9a2HzhwAEuWLEFwcDBiY2Mhk8kwduxYPPvss0hMTISFhQUAYOPGjYiJicGKFSswffr0Bq918OBB\nLF68uE5oUCgUePXVV5Geno4DBw7Ay8urwXMwSBARkalpcpDYt28fAGDMmDGwtbXVfm2ol156yeC+\nQ4YMgUajwcmTJyH+n3f6x4wZg/z8fFy+fBkbN25EbGwsYmNjERwcrO1TXV2NIUOGwM3NDQcOHGjw\nWqGhobhy5QpOnToFe3t7bXtGRgZCQ0Mxc+ZMvPfeew2eg0GCiIhMTZPnSPxvEGhKMGgKlUqF8PBw\nmJub1wsRACCRSKBQKKBUKnHu3DkAQEBAQJ0+lpaWGDBgANLS0lBWVoaOHXXfbE1NDa5cuQJvb+86\nIQIA+vXrB2tra+01iIiIqHEGL+m3fPly7N69u9F+X3zxRZ3HE40xMzNDWFgYpk2bVu9YTk4Obt26\nBQ8PD0gkEty5cwdOTk46J1W6uroCAHJzc/VeSyaTQalUwsPDQ2cdXbt25TLgRERETWBwkNi3b59B\nv62np6fjwoULj1UUAKjVaqxevRpqtRqTJ08GAJSWluodbahtLy8v13vO0tLSOn11nUMul0Op5LbS\nREREhtD7aGPx4sW4d+9enbZTp05h5syZek9WXl6Oa9euoXv3x9tNUaPRYOXKlTh9+jR8fHwQFhYG\nAFAqlXrfyqhtr66u1nve2oBgyDl0Lf1dq1MnG5ibmzV+I0RERO2c3p+WwcHBeOedd7Rfi0QiFBcX\no7i4uOETmptj0aJFggtSKpVYsWIFkpOT4e7ujpiYGO0PeCsrKygUCp2fq6mpAQBYW1vrPbelpSUA\nNHgOkUjU4DkA4MGDykbvg4iIgF8f/YZjv53E9fvZUGvU6GHvjuGuQ+Dj9ExLl0ZN1OTJln/5y1/Q\nvXt3qNVqaDQaTJ8+Hc8//zwWLFigs79IJIKlpSXc3NzqTWQ0lFwuR0REBFJSUiCVSrF9+3a4uLho\nj9vZ2aGsTPcbE7Xt+h5bANDWpe/xR1lZGWxsbHRO+iQioqY5ITuNpBv7ocF/Xw7MKrmBrJIbGOY6\nBFO8JrZgdWQsDa5s6evrq/3zSy+9BD8/P/j7+zdLIQ8fPsS8efOQmZmJPn36IC4uDp07d67TRyqV\n4ty5c6iqqoKVlVWdYzKZDGKxGD169NB7DVdXV1hYWCAvL6/eMZVKhYKCAnh6ehrnhoiITNhvZXfr\nhYg/SpWdgtTOHYO7Nc/PFHpyDP7V+6OPPkJISEizFFFdXY3w8HBkZmYiICAAiYmJ9UIEAPj7+0Ot\nVuP8+fP1Pn/p0iU8/fTTsLW11Xsdc3Nz9O/fH1lZWfVGJS5fvgy5XF4nPBERkTCpeSf1hohax/NO\nPqFqqDm1ijH8yMhIZGRkwNfXF1u3btUbBl544QWYmZkhOjpaOycCAGJjY1FeXo4pU6Y0eq2JEyei\npqYGUVFR2jaFQoGNGzcCQLOFJSIiU3LjQU6jfe6U5UGurHoC1VBzavFNu4qKirBjxw4AwFNPPYWt\nW7fq7Dd//nx4enpizpw52Lp1KyZOnIjg4GDcvHkTx48fh5+fn/Y10Vrx8fEoKytDWFgY7OzsAAB/\n/etfsXfvXsTHxyM7Oxt9+/bFiRMncP36dcyZM6fR5bGJiKhxjY1G/LEntW0tvo34kSNH8Prrrzfa\n79y5c7Czs4NGo8HOnTuxc+dO3LlzB87OzhgzZgzeeOONehMtR44cCZlMhqNHj8LNzU3bXl5ejqio\nKBw6dAilpaXw8PDA1KlTMXXqVIMmWnKJbCKihm3/ZSfOF15qsE+3Di74++B3GuxDrUeT99og/Rgk\niIgaduvhbWy4ENNgnym9X8Iwt8AnVBE9Ln1BolXMkSAiovblKXspXug5Vu9x3y798Lzr4CdYETUX\njkgIwBEJIiLDXC2+hmO/peHGg5vQQAN32+4Y5jYEz3UbCLGIv8u2JXy0YUQMEkRETaPW/L64oZmY\n2wu0VU1e2ZKIiMhYxCIxIGrpKqg5cFyJiIiIBGOQICIiIsEYJIiIiEgwBgkiIiISjEGCiIiIBGOQ\nICIiIsEYJIiIiEgwBgkiIiISjEGCiIiIBGOQICIiIsEYJIiIiEgwBgkiIiISjEGCiIiIBGOQICIi\nIsEYJIiIiEgwBgkiIiISjEGCiIiIBGOQICIiIsEYJIiIiEgwBgkiIiISjEGCiIiIBGOQICIiIsEY\nJIiIiEgwBgkiIiISjEGCiIiIBGt1QaKwsBD+/v6Ij4+vd6yyshIbNmzAqFGj8Oyzz+JPf/oTNm3a\nBLlcbvD53333XXh5een8b/369Ua8EyIiovbPvKUL+KOKigosWrQI5eXl9Y7J5XLMmDEDV69ehaen\nJ0JDQyGTyRAVFYXU1FR8+eWXsLW1bfQa169fh5OTE0JDQ+sd8/f3N8p9EBERmYpWEyRkMhkWLVqE\nX375RefxuLg4XL16FaNGjcJnn30GiUQCANi/fz+WLl2KTZs2YenSpQ1eQ6FQIDc3FyNGjMCiRYuM\nfg9ERESmplU82oiPj8f48eNx/fp1PPfcczr7/PDDDxCJRFi1apU2RADAxIkT0bdvX3zzzTeoqqpq\n8Do5OTlQKBTw8vIyav1ERESmqlUEiYSEBLi6uuLrr7/GhAkTdPbJy8tDt27d4OLiUu+Yt7c35HI5\nrl692uB1bty4AQAMEkREREbSKoLEP//5T+zfvx9+fn56+0gkEigUCp3HaudUyGSyBq9TGyRyc3MR\nGhoKX19fBAYGYvny5SgsLBRYPRERkelqFUFi6NChMDMza7CPj48PioqKkJGRUae9vLwc586d0/65\nIbVBIiYmBm5ubpgyZQqkUimSk5MREhKCgoKCx7gLIiIi09NqJls2Zvbs2UhPT8fbb7+Nf/zjHxg0\naBBkMhk+/PBD7dwIjUbT4DmsrKwglUoRHR2NXr16ads3b96Mzz77DB988AGio6MbraVTJxuYmzcc\nfIiIiExBmwkSwcHBWLp0KTZs2IDw8HBt+8CBA/Haa68hMjISVlZWDZ5j06ZNOtvDw8OxZ88eHDt2\nDBUVFejQoUOD53nwoLLpN0BERNSGOTt31NneZoIEAMyZMwdjxoxBamoqKisr0bdvXwQGBmoDgrOz\ns6DzisVieHt7Iy8vDwUFBfD09DRm2UREJqu0+iHSZGdw7f5/oNaoILXzwFDXQHS37drSpZGRtKkg\nAQDu7u6YNm1anbbatzWefvppvZ+rqqpCdnY2JBIJvL29dR4HAEtLSyNWS0Rkuq4UZ+HLqzugUP93\novydMhlOyM5g4tPjMNpjeAtWR8bSKiZbGmLdunUYNGgQiouL67Q/evQIZ86cgaenJ1xdXfV+/t69\newgJCcG7775b75hcLkdWVhYcHR0bPAcRERnmXmUxvrz6dZ0QUUsDDfbd/AFXi6+1QGVkbG0mSHh6\neuLRo0f49ttvtW1qtRofffQR5HI55s+f3+DnPTw84O3tjezsbHz33Xfado1Ggw0bNuD+/fuYOnUq\nRCJRs90DEZGpSJWdgkKtbLDP0TupT6gaak5t5tHGiy++iJ07dyI6OhrXr19Hjx49cPr0afzyyy8Y\nP348XnzxxTr94+PjUVZWhrCwMNjZ2QH4fb2KsLAwLFmyBIcPH4arqyvOnz+Pq1evYtCgQXjttdda\n4taIiNqdX4qvN9onuzQH1aoaWJpJGu1LrVebGZEwNzfHtm3bMG3aNGRlZWHHjh1QqVRYtWoVPv74\nY4jFdW8lISEB0dHRePTokbZtwIAB2L17N8aOHYtz585hx44dqKioQEREBLZt21Zn6W0iIhKusdGI\nWioD+1HrJdI0tvgC1VNUVNbSJRARtWqxl7fjSiNzIBytOuH9wGV8pNxG6Hv9s82MSBARUdsx1DWw\n0T7Pdx/MENEOMEgQEZHR9e3sjSHdBuk9/pS9FMHuQ59gRdRc+GhDAD7aICJqnEajQYrsFFJ+O4l7\n8t9f3e9oYYug7gH4k3QUJGYWLVwhNYW+RxsMEgIwSBARGU6j0aCk6j7UGjUcrTrBXNxmXhikP2gX\nS2QTEVHbIxKJ4GTduaXLoGbCORJEREQkGEckiIioWSnVSvxWdhdqjRrdbV1gbW7d0iWRETFIEBFR\ns1CpVTh0+yhOyE6jXFEBAJCILTCoqx8mev4ZNhY2LVwhGQMnWwrAyZZERA1Ta9SIu/o1Mouu6jzu\natsNb/u9xtGJNoQLUhER0RNzuegXvSECAGTl+Th658QTrIiaC4MEEREZXdrd9Eb7nLqbDrVG/QSq\noebEIEFEREaXX1HYaJ+HNWWQK6ueQDXUnBgkiIjI6CzNLBvtI4IIFmKubtnWMUgQEZHR9Xfu22if\nPp29uEx2O8AgQURERjfMNRBWZlZ6j4sgwmiP4U+wImouDBJERGR0nawcsKD/bNjoeL3TTGSGad6T\n0LuTZwtURsbGdSQE4DoSRESGkSurkF5wAddKsqGGGtKO7ghyHQwHS/uWLo2aiLt/GhGDBBERmRou\nSEVERERGxyBBREREgjFIEBERkWAMEkRERCQYgwQREREJxiBBREREgjFIEBERkWAMEkRERCQYgwQR\nEREJxiBBREREgjFIEBERkWDca4OIiIgE44gEERERCcYgQURERIIxSBAREZFgDBJEREQkGIMEERER\nCcYgQURERIIxSFCrlZycDC8vL0RFRTXYz8vLCyNHjmz2ekaOHImBAwc2+3WIWpOoqCh4eXkhOTm5\npUuhVopBgoiIiARjkCAiIiLBGCSIiIhIMAYJalfS09Ph5eWFNWvW1Du2bNkyeHl54dq1a3Xaf/zx\nR0yaNAm+vr4YOnQo1q9fj1OnTul9LpyTk4OFCxfC398ffn5+mDt3br1zErVnRUVF+Mc//oFhw4bB\nx8cHI0eOxLp161BeXl6nn1KpRHR0NMaPH48BAwYgICAAc+fOxenTpwX1A4Dy8nKsX78eo0ePho+P\nD4YOHYpVq1ahpKSkWe+Z9DNv6QKIWtJXX32FDz/8EM7OzpgwYQIUCgW+/vpr/Pjjjzr7V1VVITQ0\nFFKpFFOmTEFubi5+/vlnZGRk4NChQ3BxcXnCd0D0ZN29exdTp05FYWEhgoOD4enpiWvXriEuLg6n\nTp3Cjh07YGNjAwBYvXo1vv32WwQEBGDYsGEoKyvDv/71L8ydOxfbt2/H4MGDm9SvrKwMr7zyCrKz\nsxEYGIixY8ciLy8Pu3btwokTJ/Dtt9+iS5cuLfa9MVUMEtTqnT17ttE3N4QoKCjAhg0b4OHhgW++\n+QZOTk4AgGnTpmHy5Mk6P6NQKPDXv/4V77//vrZtzZo1SEhIwKFDhzBr1iyj10nUmvzjH/9AYWEh\nYmNjMWLECG17QkIC1qxZg+joaCxZsgTl5eXYtWsXBg0ahMTERG2/kJAQTJo0CTt27MDgwYMN7gcA\nkZGRyM7OxsqVKzFt2jRt36NHj2LhwoVYs2YNNm7c2PzfBKqDQYJavbNnz+Ls2bNGP++hQ4dQXV2N\n8PBwbYgAgD59+uCll17Crl27dH5uwYIFdb4eOXIkEhIS8Ntvvxm9RqLW5N69e0hNTcXw4cPrhAgA\nmD59OrZt24Z9+/ZhyZIlUKvV0Gg0yM/PR1FREZydnQEAzz77LI4cOYKuXbsCgMH9lEol9u/fj169\netUJEQAwatQo+Pn54aeffkJ5eTlsbW2b+TtBf8QgQa3eG2+8gUWLFuk97uXlJei8V65cAQD069ev\n3jE/Pz+dQUIikaBbt2512hwcHAAAlZWVguogaiuysrKg0WhQWlqqc5TQwsIC+fn5KCwshIuLC8aN\nG4cffvgBwcHB8PX1xbBhwxAcHIynn35a+xk7OzuD+uXm5qKyshIqlUrntaurq6FSqXDjxg34+/s3\nzzeAdGKQIJP14MEDAKgzGlFL33NWS0tLvefTaDTGKYyolXr06BEA4NKlS7h06ZLefqWlpXBxccHH\nH38MHx8fJCcna0cW169fDx8fH3zwwQd45plnAMCgfrXXvnXrFqKjo/Ve++HDh0a8YzIEgwS1KyKR\nCIDuH+pyubzO17XDnxUVFXB0dKxz7H9nnxO1d3fu3MH58+fh6+uLnj17atvVajUAwMrKChKJBACw\ncOFCRERENHpOCwsLzJkzB3PmzMHdu3dx8uRJ/Pjjj0hLS0N4eDiOHj0KCwsLg/p16NABADBhwgR8\n8sknzfAdIKH4+ie1KxYWFgB0P2b43zkMffv2BQBcvny5Xt/MzMxmqI6o9Tp69CiWL1+Oixcv1mmv\nHQlwcHDQPka8evWqznPExMRg69atqKmpwa1bt7Bu3TocO3YMANC9e3eEhITgyy+/xHPPPYfCwkLk\n5eUZ3K9nz56QSCT45ZdfdP6isGPHDmzevFk70khPDoMEtSs9evSAmZkZzpw5U2cE4vjx4/jll1/q\n9B0/fjwsLCwQGxuL+/fva9v/85//ICkp6YnVTNQa1D5mSE1N1bbV1NTg1KlTMDMzg4+PD9zd3TFo\n0CCkpqbWe0X6u+++w8aNG5GamgqJRAILCwvExcVh48aNqKmpqXPOoqIiSCQSODs7G9zP0tIS48aN\nw82bN7F9+/Y61z579izWrFmDPXv2wN7evjm+PdQAPtqgdsXR0RGjR4/Gv//9b4SEhGD48OH47bff\n8PPPP8Pf3x8XLlzQ9nV1dcWbb76JDRs2YMKECRg1ahSqqqrw73//WzsXQixm1ibTEBAQgGeffRY/\n/vgj5s6dC29vb6SlpeHWrVuYMWMG7OzsAADvv/8+pk2bhoiICAwbNgy9evVCbm4ujh8/DgcHB6xa\ntQoA4O7ujilTpiApKQkvvPAChg8fDrFYjBMnTiAnJwevv/46bG1tYWtra1A/AFi6dCkyMjLw8ccf\n4+jRo+jXrx8KCwtx+PBhmJub48MPP+S/2RbAIEHtzocffoguXbrgxx9/RGJiInr16oXPP/8cd+7c\nqRMkAGD+/Pno3LkzvvrqK+zduxcODg4ICwuDo6Mj1qxZA2tr6xa6C6InSywWY+vWrVi3bh2OHz+O\n9PR0uLq6YvHixZg7d66231NPPYXk5GTExMQgJSUFp0+fRpcuXTBhwgS8/vrrcHd31/ZdsWIFevXq\nhb1792Lfvn1QqVTo1asXPv74Y0ycOLHJ/RwdHbFr1y5s2bIFP/30ExITE+Ho6IiRI0di4cKF8Pb2\nfjLfLKpDpOFUczJRDx48gEql0vnWxueff45NmzZh9+7dOl8PJSKi33EMiExWeno6goKC6r1Kdv/+\nfezbtw/29vb8DYeIqBF8tEEma+jQoXB1dcWmTZtw5coV9O7dGw8fPsSRI0fw4MEDrF27Vvu6GxER\n6cZHG2TS7t27h7i4OBw/fhwFBQWwsbGBj48P5s6di8DAwJYuj4io1WOQICIiIsE4R4KIiIgEY5Ag\nIiIiwRgkiIiISDAGCSIiIhKMQYKIiIgEY5AgIiIiwRgkiIiISDAGCSJqsmXLlsHLywvnzp3Da6+9\nhn79+mHIkCHYu3cvvLy8MGvWLJ2f69OnD7y8vOq0jRw5EoMHD0ZVVRUiIyMxcuRI+Pj4YOTIkYiM\njERlZaXOa2dnZ2PPnj2YOHEi+vXrh+eeew7vvvsu8vLydF47NTUVc+bMwaBBg9CvXz+MHz8eX375\nZZ2tqwEgOTkZXl5eSExMxCeffAJ/f3/4+/trd7UkoroYJIhIsL///e/IyMjAsGHDYG1tjbKyMkHn\nUavVmDdvHrZv3w5XV1cEBQWhuLgYW7Zswd/+9jedn9m4cSPee+89iMViDBs2DGZmZvjuu+/wyiuv\noLy8vE7fTZs2Yd68eTh79ix69eqFYcOGobi4GJ988gleffXVemECAL7++mt89dVXGDhwIKRSKXr2\n7Cno3ojaO+61QUSClZSU4ODBg+jWrRvUajXOnTsn6DyPHj2CTCbDgQMH8NRTTwEAsrOzMWnSJBw7\ndgw5OTnw9PSs85mUlBRs3rwZI0eOBACUl5djypQpuHnzJg4ePIipU6cCAE6dOoXPP/8c3bt3xxdf\nfIFevXoBACorK/HOO+/g559/RnR0dL3Acvv2bcTExGDUqFEAfg87RFQfRySISLAxY8agW7duAACx\n+PH+d7JgwQJtiACA3r17Y9CgQQCAmzdv1uv/f//3f9oQAQC2traYMGFCvf5ffvklgN9HT2pDBADY\n2NhgzZo1sLKywo4dO+qNSri6umpDBPD490fUXvFfBhEJ9r/zHR5Hv3796rU5OzsDAORyeb1j/fv3\nr9fm5OQEANp5FSqVCufPnwcADB48uF5/R0dH9OnTB+Xl5cjKyqpzjFvIExmGjzaISDAHBwejncvO\nzq5em5mZGQDdjxU6duyot3/tXoSlpaWoqqoCAPj7+zd4/fz8fAwYMED7tb29vYGVE5k2BgkiEkwk\nEhncV6PRNDjPoCnnMrS/SqUCAFhbW2P06NEN9q0d/ajFRxlEhmGQICKjqf3hW/sD/I/Ky8u1IwVP\nioODAywsLKBUKvHxxx9rRyyIyHgYuYnIaGxsbAAARUVF9Y5dunTpSZcDiUSC/v37Q6FQ4PTp0/WO\nKxQKTJ06FTNmzIBMJnvi9RG1BwwSRGQ0PXv2hIWFBXJzc5GWlqZtv3fvHj755JMWqSksLAwAsGrV\nKty4cUPbrlQq8eGHH+LixYsoKyuDq6tri9RH1Nbx0QYRGY2NjQ1CQ0ORmJiI+fPn47nnnoOFhQXS\n09MhlUrRu3dvZGdnP9Gaxo4di7CwMHz11Vd4+eWX4ePjAycnJ1y9ehX5+flwdHREZGTkE62JqD3h\niAQRGdXy5cuxdOlSSKVSnD17FteuXUNISAi+/vprWFtbt0hN/+///T9s2rQJgwYNQk5ODlJTU2Fl\nZYUZM2Zg//79ddavIKKmEWme9OwnIiIiajc4IkFERESCMUgQERGRYAwSREREJBiDBBEREQnGIEFE\nRESCMUgQERGRYAwSREREJBiDBBEREQnGIEFERESCMUgQERGRYP8f04KzZc07HzEAAAAASUVORK5C\nYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "seaborn.stripplot(data=results_melt, x='runner', y='time', s=10)\n", "plt.ylabel('time (minutes)')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Then we test if they are significantly different using the non-parametric [Mann-Whitney test](https://en.wikipedia.org/wiki/Mann%E2%80%93Whitney_U_test)." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "MannwhitneyuResult(statistic=0.0, pvalue=0.04042779918502612)" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "scipy.stats.mannwhitneyu(results['Hugh'], results['Jesse'])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Therefore, we can conclude with *P < 0.05* that Jesse's times are significantly faster than Hugh's." ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.4.3" } }, "nbformat": 4, "nbformat_minor": 2 }