{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*This post was entirely written using the IPython notebook. Its content is BSD-licensed. You can see a static view or download this notebook with the help of nbviewer at [20161015_SimpleStatisticalTestERAttendance.ipynb](http://nbviewer.ipython.org/urls/raw.github.com/flothesof/posts/master/20161015_SimpleStatisticalTestERAttendance.ipynb).*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this notebook, I am trying to answer a question arising from the lecture of the following article: J. S. Furyk et al., “Emergency department presentations during the State of Origin rugby league series: a retrospective statewide analysis,” The Medical journal of Australia, vol. 197, no. 11, pp. 663–666, 2012. You can find the PDF [here](https://www.mja.com.au/system/files/issues/197_11_101212/fur11557_fm.pdf).\n",
    "\n",
    "The authors of the study conclude they have a statistically significantly lower emergency room attendance during game days. The question is: can we replicate their conclusion from their data?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# What is this study about? "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This article is about determining whether there were fewer visits to emergency rooms during rugby event days. The data used in this article is as follows: \n",
    "\n",
    "\n",
    "> Main outcome measures: Number of patients presenting to Queensland EDs on 24 game days and 80 control days.\n",
    "\n",
    "> In total, there were 49 702 ED presentations on game days and 172 351 presentations on control days.\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's transform this data in a neat little table:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>number of days</th>\n",
       "      <th>attendees</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>rugby game</th>\n",
       "      <td>24</td>\n",
       "      <td>49702</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>control</th>\n",
       "      <td>80</td>\n",
       "      <td>172351</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "            number of days  attendees\n",
       "rugby game              24      49702\n",
       "control                 80     172351"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "df = pd.DataFrame(data=[[24, 49702], [80, 172351]], \n",
    "                  index=['rugby game', 'control'], \n",
    "                  columns=['number of days', 'attendees'])\n",
    "df"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, let's do some simple math. What's the average number of people in the ER for these two categories of days?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>number of days</th>\n",
       "      <th>attendees</th>\n",
       "      <th>attendees per day</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>rugby game</th>\n",
       "      <td>24</td>\n",
       "      <td>49702</td>\n",
       "      <td>2070.916667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>control</th>\n",
       "      <td>80</td>\n",
       "      <td>172351</td>\n",
       "      <td>2154.387500</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "            number of days  attendees  attendees per day\n",
       "rugby game              24      49702        2070.916667\n",
       "control                 80     172351        2154.387500"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['attendees per day'] = df['attendees'] / df['number of days']\n",
    "df"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We see that on average, there were less people in the ER on game days than on control days. 83 to be precise:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "83.470833333333303"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['attendees per day'].loc['control'] - df['attendees per day'].loc['rugby game']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Is 83 people (more or less), a lot, compared to roughly 2100 per day?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.039523809523809524"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "83 / 2100"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It's 4 percent. The article says that this 4 percent difference is statistically significant. How can we prove that?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This is a tough question. We need a statiscal model to answer it.\n",
    "\n",
    "Why? Because intuitively you know that there are many reasons to observe variations in numbers and doing it over 24 / 80 days averages the result, but maybe not enough. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Single day model "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's build our model for the number of patients for a single day. It is composed of:\n",
    "\n",
    "- a mean attendance\n",
    "- a standard deviation\n",
    "\n",
    "We'll use a gaussian distribution based on these two values. Let's write a function that simulates a single day attendance:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "np.random.seed(42)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def single_day(mean, std):\n",
    "    \"Simulates the number of ER attendees for a single day.\"\n",
    "    return np.random.normal(loc=mean, scale=std)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If we say that on average the standard deviation is a 100 people per day, we can simulate a distribution for a single day:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2149.6714153011235"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "single_day(2100, 100)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can do this for 24 days in a row:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[2086.1735698828816,\n",
       " 2164.7688538100692,\n",
       " 2252.3029856408025,\n",
       " 2076.5846625276663,\n",
       " 2076.5863043050817,\n",
       " 2257.921281550739,\n",
       " 2176.743472915291,\n",
       " 2053.052561406505,\n",
       " 2154.2560043585963,\n",
       " 2053.6582307187537,\n",
       " 2053.427024642974,\n",
       " 2124.1962271566035,\n",
       " 1908.67197553422,\n",
       " 1927.5082167486967,\n",
       " 2043.7712470759027,\n",
       " 1998.7168879665576,\n",
       " 2131.4247332595273,\n",
       " 2009.1975924478788,\n",
       " 1958.769629866471,\n",
       " 2246.5648768921556,\n",
       " 2077.4223699513464,\n",
       " 2106.7528204687924,\n",
       " 1957.5251813786542,\n",
       " 2045.5617275474817]"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dist = [single_day(2100, 100) for _ in range(24)]\n",
    "dist"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's ask pandas about this distribution's mean and standard deviation:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "count      24.000000\n",
       "mean     2080.898268\n",
       "std        96.697184\n",
       "min      1908.671976\n",
       "25%      2035.127833\n",
       "50%      2076.585483\n",
       "75%      2137.132551\n",
       "max      2257.921282\n",
       "dtype: float64"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.Series(dist).describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One thing that appears in the above data is that the mean of the data does not necessarily reflect the mean we supplied to the generating function (2080 instead of 2100). The same goes for the standard deviation. This is a sample size effect which disappears if we work with large samples. We'll try 10x, 100x and 1000x larger samples below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "count     240.000000\n",
       "mean     2100.234578\n",
       "std        99.303645\n",
       "min      1775.873266\n",
       "25%      2029.066005\n",
       "50%      2106.225511\n",
       "75%      2162.773027\n",
       "max      2485.273149\n",
       "dtype: float64"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.Series([single_day(2100, 100) for _ in range(240)]).describe()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "count    2400.000000\n",
       "mean     2103.999784\n",
       "std        98.471270\n",
       "min      1798.048784\n",
       "25%      2037.499539\n",
       "50%      2102.853160\n",
       "75%      2168.237130\n",
       "max      2419.310757\n",
       "dtype: float64"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.Series([single_day(2100, 100) for _ in range(2400)]).describe()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "count    24000.000000\n",
       "mean      2099.943087\n",
       "std         99.887928\n",
       "min       1707.759975\n",
       "25%       2032.403398\n",
       "50%       2100.169540\n",
       "75%       2167.442840\n",
       "max       2547.908425\n",
       "dtype: float64"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.Series([single_day(2100, 100) for _ in range(24000)]).describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We see a progression: 2080, 2100, 2103, 2099. This is a wiggly thing. That's the reason for which we need to do statistical tests to see if it's a significant thing to have a difference in 4% on average."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Running statistical tests "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "How do statistical tests work? I'll quote an illustration from Allen Downey found in his wonderful article: [\"There is still only one test\"](http://allendowney.blogspot.fr/2016/06/there-is-still-only-one-test.html)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "![onlyonetest](https://lh4.googleusercontent.com/Bud31guq0w0FvylY57VMR0zHkYqxIpYAfOqgZietyvv1n2ToNEHwHKZWYix8pwct8kDKsZKiwvOWm6PIFEL3gBIQmbakQYHwVT02nn9_H8Fht_zaSBlrRNcqwZa950Vb5nt-5B84)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2135.125"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['attendees'].sum() / (df['number of days'].sum())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Our **data** is the number of ER people on game days (49 702) and on control days (172 351).\n",
    "\n",
    "Our **model under the H0 hypothesis** (\"there is no effect of rugby games\") is that on average 2135.125 people visit the ER. Using our previously discussed random model, we can simulate a lot of samples of 24 / 80 days under this assumption.\n",
    "\n",
    "**Test statistic**: using the data from the 80 control days and 24 game days, we'll simply substract the number of attendees in the control case from the number of attendees in the game case, normalizing by the number of days in control and game cases.\n",
    "\n",
    "Let's use the [HypothesisTest class found in Allen Downey's book Think Stats2](http://greenteapress.com/thinkstats2/html/thinkstats2010.html) for this."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "class HypothesisTest():\n",
    "    def __init__(self, data):\n",
    "        self.data = data\n",
    "        self.MakeModel()\n",
    "        self.actual = self.TestStatistic(data)\n",
    "\n",
    "    def PValue(self, iters=1000):\n",
    "        \"Returns p-value of actual data based on simulated data.\"\n",
    "        self.test_stats = [self.TestStatistic(self.RunModel()) \n",
    "                           for _ in range(iters)]\n",
    "\n",
    "        count = sum(1 for x in self.test_stats if x >= self.actual)\n",
    "        return count / iters\n",
    "\n",
    "    def TestStatistic(self, data):\n",
    "        \"Test statistic for the current test.\"\n",
    "        raise UnimplementedMethodException()\n",
    "\n",
    "    def MakeModel(self):\n",
    "        pass\n",
    "\n",
    "    def RunModel(self):\n",
    "        \"Returns a simulated data sample.\"\n",
    "        raise UnimplementedMethodException()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Our only needed simulation parameter is going to be the standard deviation of the attendance or \"by how many people can attendance vary per day\". We'll use the arbitrarily value of 100 (for now)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "variance_param = 100 "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can generate our model easily with the function already defined previously:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "class ERTest(HypothesisTest):\n",
    "    \"Model for replicating the Furyk et al. (2012) study.\"\n",
    "    def __init__(self, data, variance_param):\n",
    "        super().__init__(data)\n",
    "        self.variance_param = variance_param\n",
    "        \n",
    "    def TestStatistic(self, simulated_data):\n",
    "        \"The statistic is the relative difference in control vs game day attendees.\"\n",
    "        sum_game, sum_control = simulated_data\n",
    "        test_stat = (sum_control / 80. - sum_game / 24.)\n",
    "        return test_stat\n",
    "\n",
    "    def RunModel(self):\n",
    "        \"Returns a simulated data set of 104 observations under the H0 hypothesis.\"\n",
    "        simulated_data = (sum(single_day(2135.125, self.variance_param) for _ in range(24)),\n",
    "                         sum(single_day(2135.125, self.variance_param) for _ in range(80)))\n",
    "        return simulated_data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Good, now we're ready to run our simulation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "data = (49702, 172351)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "ert = ERTest(data, variance_param)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "83.4708333333333"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ert.actual"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0001"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ert.PValue(iters=10000)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "What this means is that using the above hypotheses, the likelihood to observe our actual result, with a difference in means of 83 people per day and the given sample sizes (24 / 80) is very unlikely to be due to chance. \n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To get a better feeling for the distribution of the test statistic, let's plot it:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "plt.style.use('bmh')\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def plot_test_stat(test, title=\"\"):\n",
    "    \"Plots the test statistic distribution and observed value.\"\n",
    "    plt.hist(test.test_stats, bins=30, cumulative=False, normed=True)\n",
    "    ylim = plt.ylim()\n",
    "    plt.vlines(test.actual, *ylim, label='observed test stat')\n",
    "    plt.legend(loc='upper left')\n",
    "    plt.xlabel('test statistic: control - game days (in people per day)')\n",
    "    plt.title(title)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEdCAYAAADn46tbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmYFsW1/z9fEFQQIRgHFR1EIa6IwSUm0Zg4BtEkYtw1\nMW5ZDBpNbqKi0YvmZlOvcU00Rq9x4xJ/aBI1KhrxmowiS3AMiaiAyCIICsoisjic3x9VLzTvvO+8\nPUPPTA+cz/PMM13Vp5aurrdP1znVVTIzHMdxHCdJh7augOM4jpM/XDk4juM4DXDl4DiO4zTAlYPj\nOI7TAFcOjuM4TgNcOTiO4zgNcOWwGSJphKTX27oeTjbE+zmtreuRBZLOkrSmrevhuHLYXLkOOKSt\nK9FcJO0g6X5J/5K0RtJTjcg9KGlJ/PtfSdsXyWwh6VpJ8yStkPR3SYNa50oyZVP5YMnYdK6lXePK\nYTNCgQ5mtsLMFrd1fTaCLYFFwPXAX0sJSBLwF6APUAN8EfgE8Mci0f8Gzga+BRwIvAH8VVJVi9R8\nM0RSp7aug9N0XDnkEEnflPS+pC2L4odLmpUI3yFpenzjnSHpZ5I6J86PkDRN0smSpgKrgP7FZghJ\nu0p6SNJbkj6Q9E9JXy8q+1lJv5N0haT5khZJukdSlyK5UyRNkvShpHcl/UVS98T570maGs+/July\nSR2b0j5mNsvMLjKzu4G3y4h9Edgf+JqZTTKzCcAZwGckfS7WZRvgO8BwM/uLmb1CUBSrgPPKlS/p\ncElrJX1Z0vh4LVMkfaFIbndJoyW9J2mxpDGS9i2SOSa210pJCyT9Otmmku6W9LSkH0iaG+/Pg5I+\n1lgbSfqipNrYN+ZK+h9JPSukWSvpwljn5THdhUUyXSXdlKjLPyR9NXG+T8zn9HjvlwM/aaTM/4rX\nvVTSSOBjRecb7ZsKZqj3JG1VlG6EpNcau16nAmbmfzn7A7YFPgBOKor/F/CTeCzgvwhvu9XAl4G3\ngBEJ+RExn2eBg4B+wDYx/vWE3L7AsPi/L3A+sBo4PCHzLLCY8Lb+CeBIwtv71QmZs2O6y4E9Y37f\nA3rG81cBM4FjCW/0Q4A3i/K4CljbhLa6G3iqRPxVwPQS8bOBy+PxF4B6YOcimXtL5Zk4fziwFngN\nOBrYA7gTWA70ijJVwHzgVmBvoD9wE/AOsF2U2Q9YQxi97AEcBcwC7im6viXAn2I+nwNeBx4qus/J\n+3lEvO/DgN2AA4BngP+r0JZrgXdjun7x3q0BvlLUD8YCnwZ2Bb4JrAS+EM/3ifnMBk6P4T5lyrsI\nWAZ8PZb3I+A9YHXavglsFfvhGYk0iv3sR239W27Pf21eAf8rc2Pgf4FHE+ED44OsfyNpvg+8lgiP\nAD4CehfJbfAwKZPXn4DfJsLPAnVFMr8Bnk+EZwE3lclv6/jAGlwUfwbwXiJ8PvDvJrRTOeXwW6C2\nRPwE4JZ4fFps0y2KZK4FpjRSZkE5nJWI60hC0RGU0wtF6QRMBy6M4fuAF4tkjo112iVxfUuBbRIy\nX4zl71bqfsZ79fOifKtjmv0aua61wO+L4h4AnovHnwdWAN2KZO4CHo7HBeVweYp7N4f4spOI+38k\nlEPKvnkT8LdE+CiCwto+bT/yv4Z/W+DklXuBP0n6uJm9C3wDGG9mSXPQt4BzCW9wXYEtCA+gJAvM\n7K3GCpK0NeEB82VgR6Bz/Hu2SLSuKDwPGBzz2B7YBXi6TDH7EBTEQ8EdsI6OQGdJ25nZIjP7NfDr\nxuqbAWkcnpVkDHhxXcCsXtIEwnVCUOYHSlpWlG4rwlsyhJHAM0XnnyPcw70JD0+AV8xseULm+UT6\nN0rU7SDgU5K+V6LO/YF/NnJdLxaFn2e9WehAgr9nXtE97EQYzSSZ2EgZSOoG9AbGFZ2qBYYm5NL0\nzd8CUyTtaWavEkYzj5jZO43VwWkcVw75ZQxhuHy6pF8DpwD/WTgp6SSCyeIS4G+Et8uTgZ8W5fNB\nirL+G/gK8APCj/wD4FcE81aS1UVho6HfqtxDtSB3IlBq2mXWDvL5BEd0Mb1Y76eYH//vAMwtI9MU\nxPrr70Bwlp9PQ4W9JHFcrr02RoF1AK4hjEyKaep1FV/T+wQlUXxNxX2jUr8rpK90nRX7ppm9Iul5\n4JuSriGMvo6pkK9TAVcOOcXM1kYH3RkE++m2wB8SIocBk83spkKEpL7NLO4w4AEzeyjmI4JfIfWD\nxMzekTSXMKT/SwmRfxOG+rub2Zhm1rMpPA9cKWl3M5sBIGkvwuimNsr8g/BQO4pgGilc+5HA7RXy\nF2E68KsxXUfCG/s98fwk4ExgnpmtKpPHvwkmqiSfJ5hlXknE7SVpm8To4bOEh+rUMvlOAvYxs1Kj\nikocwobX/plEOZOAHsDWFpz3zcbMlkp6i3AtTyZOHVokmrZv/ha4keCzmGtmxSMyp6m0tV3L/8r/\nAQMID4rJwP8rOnc+wQF6LMHpeBHB2VmfkCnpWyiOJ9h5XyE83PYGfkd4QxybkHkWuKMonx8DbyTC\n5xJm+lxBcEjvE+tZcEhfEfMdRviB700YEf0ykccFwNQUbTOQMBvpz8ALMTwwcV6Eh9m4eF2fIpg6\naovyuQFYAHwp1uf3hBFbr0bKLvgcXiU4pPeMbbYc2CHKVBFGI08QHnh94v+fAock7u9qgpN/D4KD\nfhYJuz/B5/A+8HBsz88RHOEPN3I/Px/vw/WxXXaLed8JbNnIdRUc0uez3iG9Gjg2ITMmXvdQgoN4\nULxn58bzBZ/DZ1Lcw+8TRrwFh/QPCSPIpEO6Yt+MclsS+v+HhNlnbf77be9/bV4B/6twg4JiqAe+\nXBS/BXBb/DG/D9xPeOg2RznsHB9iy4gznuKPMKkcxlJBOcS404CX4o/0HeBRYNvE+bPjNa0gPITH\nAd8pqttHKdplbWyXwt/a5LVHmV6E0daS2EYjgY8XyXQEfknwn6wA/g58skLZhxfuCUEBfUiYSXZE\nkdwuBNPOgigzk+BL6pOQGUJQWh9GuVsJb+aF83cDTwH/Eeu4HHgQ+Fhj95nwRv5UvPZlhFHKr4AO\nFdr0QsK3IB/EvnBRkcyWwM+BGYSR4DzgceDz8Xyf2DZplIMIynJhrOODhJecpHKo2DcTsr8iKMWy\nit3/0v8pNmqjSBpCGLJ1AO4ys2uKzneOnf4AwsPqFDObHedVjyZo/bvN7MJEmtOAy2KHnAd83dr3\nh1nOZoKkwwnKchczm9fCZd1NmG02uCXLiWWtJfwOR7Z0WS2BpD8QZp6d0NZ12RSo+BGcpA6Et5mj\nCMPa0yTtWSR2LrDYzPoTlMi1MX4lwZTww6I8O0a5w81sf2AKYWjqOO2FYoes00ZI6iHpKOA4wujB\nyYA0X0gfDEyz8FXqGmAUialmkaGsd8SNJs4SsbBMwwuEoV6Swg+rW3QwbUsYPThOeyHNbKL2Rnu9\nppcIJqlrzOz5SsJOOtLMVurN+vnWEJxsB5eTsTDf+31JPcuZiczsI0nDCCOG5YSpjcOaWnnHaQvM\n7DmCr6I1yjq7NcqJZbXKNWWNmTV3lp7TCGlGDqWGz8VvGMUyKiGz/qS0BfBdwuyS3gQlcXmKujiO\n4zitQJqRw1zCp/cFdqahCWgOYWbGvOhP2NbM3mskz/0BM7M3Y/hB4NJSgscee6ytXLmSHXbYAYCu\nXbvSr18/9t9/fwDq6sJHux6uHC4c56U+7T3s7entmadwXV0dY8aET4h22GEHunbtym233dZs31jF\n2UrxYf8awY8wn7A2zWlmNjUhMwzY18yGSToVOM7MTk2cPxM40My+F8M7EqYA7mdmiyT9hDB97+Li\n8r/xjW/YTTfdVBztNINf/vKXDB8+vK2rscng7Zktm3J79uwZFsRdvLj1JmRedNFF3Hvvvc1WDhVH\nDtGHcAFhznRhKutUSVcDE83sMcLXpfcpLAO9CEgqhplAN8L6OUMJC6+9GtP/XdJqwoc/Z5Uq/+23\nm7OKgVOK2bNnt3UVNim8PbPF2zNfpFo+w8yeJHzBmYwbkTheRVjXp1Taks4iM7sDuCN1TR3HcZxW\nI/eb/Rx11FFtXYVNhtNPP72tq7BJ4e2ZLd6e2TJw4MCNSp/qC+m25JlnnrFBg9rjlr6O4ziBtvA5\nTJ48mZqampbzObQ1dXV1lFIOZsbChQupr69vg1q1T5YsWUL37t0rCzp07NiRqqoqivYt2IDa2loO\nPbR4EVGnuXh75ovcK4dyLFy4kG7dutGlS5fKwg4AO+20U1tXod2wYsUKFi5cSK9evdq6Ko7TJuTe\n51CYz1tMfX29KwanxejSpUvFUam/5WaLt2e+yL1ycBzHcVqf3JuVyvkcHKelWb66npfnFW8BDdt0\n7sjy1fXUTRjH/gd/utE8qrbpzI7bbtlSVdykcJ9Dvsi9ctjUOP/88+nduzeXX95+lpK65pprmDlz\nJrffXmnnzE2LNxZ/yK8mv9sgfsSRfbn6rzNZOuMttn13eqN5XHdMP1cOTrsk92alcj4Hp3UpN2vn\n+eefZ999982kjGOPPZb777+/VdJec801fPe7321WWQW23d37Zpb4qCFf5F45OE2nNaf3mlmj0z0d\nx2mf5F45JFdqbC+8/vrrHHvssfTt25fPfvazPPnkkxucX7RoEccffzzV1dUce+yxzJ07d925yy+/\nnD322INdd92Vz33uc7z66qsArF69miuvvJL99tuPvfbaix/96EesWhX2UCq8vd98883stddefO97\n3+OQQw7h6aefXpdvfX09/fv3Z8qUKQBMnDiRIUOG0LdvXw4//HCef379HimzZ8/mK1/5Cn369OGE\nE04o++HOihUrOOWUU3j77beprq6murqaBQsWYGbceOONHHDAAfTv359zzz2XJUuWALBq1SrOO+88\n+vXrR9++fTnyyCN59913+dnPfsa4ceO49NJLqa6uLrkAW1PTXnbZZQwYMIA+ffpQU1PDiy++CMAz\nzzzDDTfcwB//+Eeqq6s5/PDDm3aDI0tntL++mWdqa2vbugpOgk3S51D4GjELmvpF40cffcTpp5/O\nGWecwcMPP8y4ceP42te+xrPPPsvuu+8OwOjRo/nDH/7AAQccwH/+53/y7W9/m8cff5yxY8cyfvx4\nJk2aRLdu3Zg2bdq6j9ZGjBjBnDlzqK2tpWPHjnz729/muuuu44orrgDCdx9Llizhn//8J2vXruWW\nW25h9OjRfPGLXwTCA3G77bZjwIABzJs3j9NOO43f/va31NTU8Nxzz3HmmWcyYcIEevbsybe+9S0+\n9alP8fDDDzNp0iROOeUUvvSlLzW41i5duvDggw9y3nnnrVM6ALfddhtPPPEEf/nLX9huu+0YPnw4\nP/rRj/jd737H//7v/7Js2TL+/e9/07lzZ6ZMmcJWW23Fj3/8Y8aPH8/JJ5/M17/+9ZJt29S0Bxxw\nAMOHD6dbt27cfvvtnH322bz88svU1NTwgx/8gDfffJPbbrutSffXcTYXcj9yaG8+h0mTJrFixQou\nuugitthiCw477DCOOuooHnrooXUygwcP5pBDDqFTp05cccUVTJo0iXnz5tGpUyeWL1/Oa6+9hpnR\nv39/qqqqALj//vv52c9+xrbbbkvXrl256KKLNsizY8eODB8+nE6dOrHllltywgkn8MQTT7By5UoA\nHnroIU488UQgKKfBgwdTU1MDwOGHH87+++/P008/zdy5c6mrq+Oyyy6jU6dOfPrTn2bIkCFNaoN7\n7rmHK664gh122IFOnTpx8cUX88gjj7B27Vo6derE4sWLmTFjBpLYb7/92GabbVLl29S0J554It27\nd6dDhw4MGzaMVatWMX164w7kpuA+h2xxn0O+2CRHDq25fkkx8+fPb/Al8i677ML8+fPXhXv37r3u\nuGvXrvTo0YO3336bww47jG9+85tccsklvPXWW3zpS1/iJz/5CStXrmTFihV84QtfWJdu7dq1JNfF\n2m677ejUqdO6cN++fdljjz148sknOeqoo3jiiSfWzZCaM2cOf/rTn9aZu8yM+vp6Pve5z/H222/T\no0cPtt566w3qP29e+i2+586dyxlnnEGHDh3W5d+pUycWLlzIKaecwrx58zj33HNZunQpJ510Elde\neSUdO1beofLUU09tUtpbb72V+++/nwULFgCwfPlyFi1alPo6HGdzJvcjh/bmc9hxxx0bPEjnzp3L\njjvuuC781ltvrTtevnw577333rqd7r71rW8xduxYXnjhBaZPn84tt9zCdtttR5cuXXjhhRd44403\neOONN3jzzTeZNWvWunxKOYWPP/54HnroIR5//HH23HNP+vTpAwTldMopp6zLa+bMmcyePZsLL7yQ\nHXbYgffff58PP/xwg/qXo1S5vXv35sEHH9wg/7lz57LDDjuwxRZbcPHFFzNu3DjGjBnDmDFjGDVq\nVNm8knTs2DF12hdffJFbbrmF3//+98ycOZOZM2fSrVu3dQo1Cye6+xyyxX0O+SL3yqG9ccABB9Cl\nSxduvvlmPvroI2praxkzZgwnnHDCOpmnn36a8ePHs3r1an7+859z4IEHstNOO/HSSy/xj3/8g48+\n+oitttqKLbfcko4dOyKJM844g8svv5x33w3z7ufNm8fYsWMbrcvxxx/Ps88+y913373OpARw0kkn\nMWbMGMaOHcvatWtZuXIlzz//PPPnz2fnnXdm//3355e//CVr1qzhxRdfbOBQT7L99tvz3nvvsXTp\n0nVxZ511Fj/96U/XKZV3332XJ554AggPgFdeeYW1a9fStWtXOnXqtO7Nf/vtt99A4RXTlLTLli1j\niy22oGfPnqxevZprr72W5cuXrztfVVXF7NmzyfuqxI7TVqRSDpKGSHpV0uuSGuz1LKmzpFGSpkka\nJ6k6xveUNFbSMkk3F6XpJOm3kl6T9Iqkr5Yqu735HDp16sQDDzzA008/Tb9+/bjkkku4/fbb1zmj\nJXHiiSdyzTXX0K9fP6ZMmcIdd4Q9j5YtW8b3v/99dtttNz75yU+y3XbbccEFFwBw1VVXsdtuuzF4\n8GB23XVXTjjhBGbMmNFoXXr16sVBBx3EpEmT+OpX1zdv7969uf/++7nhhhvo378/AwcO5NZbb2Xt\n2rUA3HHHHUyaNIndd9+d6667jtNOO61sGf379+f4449n0KBB7LbbbixYsIDzzjuPo48+mhNOOIE+\nffowZMgQJk+eDMCCBQs4++yz2XXXXfnMZz7DoYceysknh32ivvOd7/DnP/+Z3Xffncsuu6xBWU1J\ne+SRR1JTU8NBBx3EJz/5SbbeeusNzHlDhw7FzNh999054ogjGr+pZXCfQ7a4zyFfpNlDugPwOmEP\n6XnAROBUM3s1IfNdYEDcQ/oU4KtmdqqkLsD+wL6EPaYvTKS5CuhgZv8Zwz3NrIGzoNx+DvPmzfNV\nRp0W5cnJr/OryR80iC98IZ2G647px8CdumVdNaed0R73c0gzcjgYmGZms8xsDTAKGFokMxS4Jx6P\nJigSzGyFmb0ArCqR7znALwqBUooB2p/Pwdl8cJ9DtrjPIV+kUQ69gTmJ8NwYV1LGzOqB9yWV/dhA\nUmHHmZ9K+oekP0jaPn21HcdxnJYkzVTWUsOSYltUsYxKyBSXuzPwdzP7oaQfANcD3ygWnD59OsOG\nDaO6uhqA7t27M2DAAHbbbbcUVXecjaMwOij4F5bOqGNKjwVAL7bdff+S55PhugnjWPbxLuvs6YW3\nYw83DB966KG5qk+W4QItWV5tbS0jR44EoLq6mqqqqnXfMjWHND6HQ4CrzGxIDA8HzMyuScg8EWXG\nS+oIzDezqsT5M4EDinwOy8ysWzzeGXjCzAYUl+8+B6etcJ+DkxWbqs9hItBPUh9JnYFTgUeKZB4F\nzozHJwGl5lgWV/JRSYWvuo4EXilVuPscnLZibYUXJ/c5ZIv7HPJFRbOSmdVLugB4iqBM7jKzqZKu\nBiaa2WPAXcB9kqYBiwgKBABJM4FuQGdJQ4HBcabT8JjmBuAd4OymVLxjx46sWLHCtwp1WoQlyz5g\n9vtr2roajtNmpFo+w8yeBPYoihuROF4FnFwmbd8y8bOBisthlvvOoaqqioULF/L+++9XysLZjFm+\nup43Fn9YWRDo02MrZr2/krVmzH5/DX+cuZrSLreAf+eQLf6dQ75ot2srSaJXr15tXQ0n57w8b1nJ\n3dxKMeLIKn41Obn2ku9T4Wy+5F45+B7S2eF79GbL0hl1FUcPHTtQch/qYnyvae+feSP3ysFx2jNL\nVtanmtnke007eSP3C++1t7WV8oy/lWWL+xyyxftnvsi9cnAcx3Fan9wrB//OITt8Hnm2+HcO2eL9\nM1/kXjk4juM4rU/ulYP7HLLDbbrZ4j6HbPH+mS9yrxwcx3Gc1if3ysF9DtnhNt1scZ9Dtnj/zBe5\nVw6O4zhO65N75eA+h+xwm262uM8hW7x/5ovcKwfHcRyn9cm9cnCfQ3a4TTdb3OeQLd4/84WvreQ4\nOSDtAn3gi/Q5rUPulYP7HLLDbbrZkqXPIe0CfbDpLtLn/TNfpDIrSRoi6VVJr0u6tMT5zpJGSZom\naZyk6hjfU9JYScsk3Vwm70ck/XPjLsNxHMfJkorKQVIH4FbgKGAf4DRJexaJnQssNrP+wI3AtTF+\nJXAF8MMyeX8VWNpY+e5zyA636WaL+xyyxftnvkgzcjgYmGZms8xsDTAKGFokMxS4Jx6PBmoAzGyF\nmb0ArCrOVFJX4AfAT5tZd8dxHKeFSKMcegNzEuG5Ma6kjJnVA+9L6lkh3/8C/htodINf9zlkh9t0\ns8W/c8gW75/5Io1DutRGulZBRiVk1p+UBgL9zOw/JO1apgwARo8ezZ133kl1dTUA3bt3Z8CAAes6\nUmEo6mEPlwrXTRjH0hlvrXuQF0xBGxOe0mMB0CuV/JRJL7J0xvyK+XNk39Tl1014h4HHDc5F+3o4\nXbhAS5ZXW1vLyJEjAaiurqaqqoqamhqai8zKPsODgHQIcJWZDYnh4YCZ2TUJmSeizHhJHYH5ZlaV\nOH8mcICZXRjD5xF8EauBTkAV8LyZHVFc/vXXX2/nnHNOsy/QWc/muEfvy/OWcfHj01PJjjiyb6oZ\nQwW5NHtINzXPNFx3TD8G7tQtlWx7YlPunz17BkPK4sWLW63MyZMnU1NTU/bFuxJpzEoTgX6S+kjq\nDJwKPFIk8yhwZjw+CRhbIp91lTSz281sZzPbDTgUeK2UYnAcx3HahopmJTOrl3QB8BRBmdxlZlMl\nXQ1MNLPHgLuA+yRNAxYRFAgAkmYC3YDOkoYCg83s1bQVdJ9Ddmyqb2VthfscssX7Z75I9RGcmT0J\n7FEUNyJxvAo4uUzavhXyngXsl6YejuM4TuvgayttRvg88mzx7xyyxftnvsi9cnAcx3Fan9wrB/c5\nZIfbdLPFfQ7Z4v0zX+ReOTiO4zitT+6Vg/scssNtutniPods8f6ZL3KvHBzHcZzWJ/fKwX0O2eE2\n3Wxxn0O2eP/MF7lXDo7jOE7rk3vl4D6H7HCbbra4zyFbvH/mi9wrB8dxHKf1yb1ycJ9DdrhNN1vc\n55At3j/zRe6Vg+M4jtP65F45uM8hO9ymmy3uc8gW75/5ItWqrI6TN+YvXcXC5asryq2uX9sKtXGc\nTY/cKwf3OWTHpmTTXbh8daod3kYc2eiK8RuF+xyyZVPqn5sCuTcrOY7jOK1PKuUgaYikVyW9LunS\nEuc7SxolaZqkcZKqY3xPSWMlLZN0c0J+a0mPSZoqaYqkn5cr230O2eE23Wxxn0O2eP/MFxWVg6QO\nwK3AUcA+wGmS9iwSOxdYbGb9gRuBa2P8SuAK4Iclsr7OzPYCPgkcKumo5l2C4ziOkzVpRg4HA9PM\nbJaZrQFGAUOLZIYC98Tj0UANgJmtMLMXgFVJYTP70Myei8cfAZOBnUsV7j6H7HCbbra4zyFbvH/m\nizTKoTcwJxGeG+NKyphZPfC+pJ5pKiCpB/AV4Jk08o7jOE7Lk2a2kkrEWQUZlZBpmLHUERgJ3Ghm\nb5aSuemmm+jatSvV1dUAdO/enQEDBqx7yyjYKT1cOZy06eahPhsT7rbbQGC93b/wFl8cnjLpRZbO\nmF/2fHPCU3osAHpt4HPY2PKJs6rSlF834R0GHje4RdvX+2e24QIt3X4jR44EoLq6mqqqKmpqamgu\nMmv8GS7pEOAqMxsSw8MBM7NrEjJPRJnx8YE/38yqEufPBA4wswuL8r4LWGpmPyhX/vXXX2/nnHNO\nMy7NKaa2tnaTGbq/PG9Z6qmsV/91Zqo808oW5JbOqKtoWmpqnmn41Zf7kebzjaptOrPjtlumyjMP\nbEr9s5iePYMhZfHixa1W5uTJk6mpqSn1cp+KNCOHiUA/SX2A+cCpwGlFMo8CZwLjgZOAsSXy2aCS\nkn4KbGtm5zZWuPscsmNT/eG1FW3lc1iysj6VIrnumH7tSjl4/8wXFZWDmdVLugB4iuCjuMvMpkq6\nGphoZo8BdwH3SZoGLCIoEAAkzQS6AZ0lDQUGA8uAy4Gpkl4imKBuNbP/yfbyHMdxnOaQ6gtpM3sS\n2KMobkTieBVwcpm05T5RTfWNRV1dHYMGDUoj6lRgUx62twVpzEpOerx/5gv/QtpxHMdpQO6Vg/sc\nssPfyrLFRw3Z4v0zX+ReOTiO4zitT+6Vg6+tlB2+dk22+NpK2eL9M1/kXjk4juM4rU/ulYP7HLLD\nbbrZ4j6HbPH+mS9yrxwcx3Gc1if3ysF9DtnhNt1scZ9Dtnj/zBe5Vw6O4zhO65N75eA+h+xwm262\nuM8hW7x/5ovcKwfHcRyn9cm9cnCfQ3a4TTdb3OeQLd4/80XulYPjOI7T+uReObjPITvcppst7nPI\nFu+f+SL3ysFxHMdpfVIpB0lDJL0q6XVJl5Y431nSKEnTJI2TVB3je0oaK2mZpJuL0gyS9M+Y543l\nynafQ3a4TTdb3OeQLd4/80XFzX4kdQBuBWqAecBESX82s1cTYucCi82sv6RTgGsJu8GtBK4A9o1/\nSW4DvmlmEyQ9LukoMxuz8ZfktFfmL13FwuWrU8muTrOJsuM4zSbNTnAHA9PMbBaApFHAUCCpHIYC\nhZ3hRhPkLHKvAAAgAElEQVSUCWa2AnhBUv9khpJ2ALqZ2YQYdS9wHNBAObjPITvybtNduHw1Fz8+\nPZXsiCPLbTDYerjPIVvy3j83N9KYlXoDcxLhuTGupIyZ1QPvS+pZIc+5FfJ0HMdx2og0ykEl4qyC\njErINDVPwH0OWeI23Wxxn0O2eP/MF2nMSnOB6kR4Z4LvIckcYBdgnqSOwLZm9l6FPHepkCcAzz33\nHJMmTaK6OlShe/fuDBgwYN0QtNChPLxphAsP3ILJplyYaFaqJD9l0ossnTG/Yn5NCU/psQDolWn5\naa+nqeW39f308IYKryXLq62tZeTIkQBUV1dTVVVFTU0NzUVmjb3gQ3zYv0ZwSM8HJgCnmdnUhMww\nYF8zGybpVOA4Mzs1cf5M4EAz+14ibjzwPWAi8BfgZjN7srj8Z555xgYNGtTsC3TaDy/PW9Ykn8PV\nf52ZmdymmOd1x/Rj4E7dUuXptCw9ewYr++LFi1utzMmTJ1NTU1PKSpOKiiMHM6uXdAHwFMEMdZeZ\nTZV0NTDRzB4D7gLukzQNWESYqQSApJlAN6CzpKHA4DjTaRjwe2Ar4PFSisFxHMdpG9KYlYgP7j2K\n4kYkjlcBJ5dJW3JaiZn9AxhQqey6ujp85JANtbW1PiMkQ5bOqPMZSxni/TNf+BfSjuM4TgNyrxz8\nO4fs8LeybPFRQ7Z4/8wXuVcOjuM4TuuTe+Xg3zlkh88jzxb/ziFbvH/mi9wrB8dxHKf1yb1ycJ9D\ndrhNN1vc55At3j/zRe6Vg+M4jtP65F45uM8hO9ymmy3uc8gW75/5ItVHcI7jtD86dghLkqShapvO\n7Ljtli1cI6c9kXvl4D6H7HCbbrbk3eewZGV96vWarjumX5srB++f+SL3ZiXHcRyn9cm9cnCfQ3a4\nTTdb3OeQLd4/80XulYPjOI7T+uReObjPITvcppstefc5tDe8f+aL3CsHx3Ecp/XJvXJwn0N2uE03\nW9znkC3eP/NFKuUgaYikVyW9LunSEuc7SxolaZqkcZKqE+cui/FTJQ1OxP9A0r8k/VPSA5I6Z3NJ\njuM4zsZSUTlI6gDcChwF7AOcJmnPIrFzgcVm1h+4Ebg2pt2bsEPcXsDRwG8U2Imwf/QgM9uP8L3F\nqZTAfQ7Z4TbdbHGfQ7Z4/8wXaUYOBwPTzGyWma0BRgFDi2SGAvfE49HAEfH4WGCUmX1kZm8C02J+\nAB2BrpK2ALoA85p9FY7jOE6mpFEOvYE5ifDcGFdSxszqgSWSepZI+xbQ28zmAdcDs2Pc+2b211KF\nu88hO9ymmy3uc8gW75/5Io1yUIk4SylTMl5SD8Joow+wE7CNpNNT1MVxHMdpBdKsrTQXqE6Ed6ah\nCWgOsAswT1JHoLuZvSdpbowvTnsk8IaZLQaQ9DDwGWBkceHTp09n2LBhVFeHKnTv3p0BAwass08W\n3jY8XDl86KGH5qo+pcKFt/GCPb9cmCP7ppKfMulFls6YXzG/poSn9FgA9GLb3ffPrPy015MsP6vr\nCfQDvH+2VLhAS5ZXW1vLyJHhEVpdXU1VVRU1NTU0F5kVDwKKBMLD/jWgBpgPTABOM7OpCZlhwL5m\nNkzSqcBxZnZqdEg/AHyKYGJ6GugPHATcFf+vAu4GJprZr4vLf+aZZ2zQoEHNvkCn/fDyvGVc/Pj0\nVLIjjuybalG5tHKbe57XHdOPgTt1SyXrNJ2ePXsCsHjx4lYrc/LkydTU1JSy3qSiolkp+hAuAJ4C\n/k1wME+VdLWkL0exu4CPS5oGfB8YHtO+AjwIvAI8DgyzwASC4/ol4GWC+emOUuW7zyE73KabLe5z\nyBbvn/ki1ZLdZvYksEdR3IjE8SrClNVSaX8B/KJE/NXA1U2prOM4jtM65P4Laf/OITt8Hnm2+HcO\n2eL9M1/kXjk4juM4rU/ulYP7HLLDbbrZ4j6HbPH+mS9yv02o0/6Zv3QVC5evrii3un5tK9TGcZw0\n5F45uM8hO9rKprtw+epUU1RHxLn+7QX3OWSL+xzyRe7NSo7jOE7rk3vl4D6H7HCbbra4zyFbvH/m\ni9wrB8dxHKf1yb1ycJ9DdrhNN1vc55At3j/zRe6Vg+M4jtP65F45uM8hO9ymmy3uc8gW75/5IvfK\nwXEcx2l9cq8c3OeQHW7TzRb3OWSL9898kXvl4DiO47Q+uVcO7nPIDrfpZov7HLLF+2e+yL1ycBzH\ncVqfVMpB0hBJr0p6XdKlJc53ljRK0jRJ4yRVJ85dFuOnShqciO8u6f/F+H9L+lSpst3nkB1u080W\n9zlki/fPfFFROUjqANwKHAXsA5wmac8isXOBxWbWH7gRuDam3ZuwQ9xewNHAbyQV9jS9CXjczPYC\nBgJTcRzHcXJBmpHDwcA0M5tlZmuAUcDQIpmhwD3xeDRwRDw+lrDn9Edm9iYwDThYUjfgMDO7GyCe\nX1qqcPc5ZIfbdLPFfQ7Z4v0zX6RRDr2BOYnw3BhXUsbM6oElknqWSPtWjNsNeFfS3ZImS7pD0tbN\nvAbHcRwnY9Ls56AScZZSplz8FsAg4HwzmyTpRmA4MKJYePr06QwbNozq6uDG6N69OwMGDFhnnyy8\nbXi4cvjQQw9tk/JnvLsC2B5Y/7ZdsNc3N0zc+6GS/JRJL7J0xvyNLi8ZntJjAdCLbXffP7Py015P\nsvws27Njh368PG8ZdRPGAbD/wZ8GaBB+c8oktuvaaZPqn60RLtCS5dXW1jJy5EgAqqurqaqqoqam\nhuYis+LnfJGAdAhwlZkNieHhgJnZNQmZJ6LMeEkdgflmVlUsK+lJggJ4ExhnZrvF+EOBS83sK8Xl\nP/PMMzZo0KBmX6DT9rw8b1nqzX6u/uvMVHmmlfU8s83zumP6MXCnbqnydNbTs2dPABYvXtxqZU6e\nPJmamppSL+ipSGNWmgj0k9RHUmfgVOCRIplHgTPj8UnA2Hj8CHBqnM3UF+gHTDCzBcAcSZ+IcjXA\nK6UKd59DdrhNN1vc55At3j/zRUWzkpnVS7oAeIqgTO4ys6mSrgYmmtljwF3AfZKmAYsICgQze0XS\ng4QH/xpgmK0fqlwIPCCpE/AGcHbG1+Y4juM0k1R7SJvZk8AeRXEjEserCFNWS6X9BfCLEvEvAwdV\nKtu/c8gOn0eeLf6dQ7Z4/8wX/oW04ziO04DcKwf3OWSH23SzxX0O2eL9M1/kXjk4juM4rU/ulYP7\nHLLDbbrZ4j6HbPH+mS9yrxwcx3Gc1if3ysF9DtnhNt1scZ9Dtnj/zBe5Vw6O4zhO65PqO4e2xH0O\n2ZGlTXf+0lUsXL46lezq+rWZlZsn3OeQLe5zyBe5Vw5OPlm4fHWq9ZIgrNvjOE77IvdmJfc5ZIfb\ndLPFfQ7Z4v0zX+ReOTiO4zitT+6Vg/scssNtutniPods8f6ZL3KvHBzHcZzWJ/fKwX0O2eE23Wxx\nn0O2eP/MF7lXDo7jOE7rk3vl4D6H7HCbbra4zyFbvH/mi1TKQdIQSa9Kel3SpSXOd5Y0StI0SeMk\nVSfOXRbjp0oaXJSug6TJkoq3HXUcx3HakIrKQVIH4FbgKGAf4DRJexaJnQssNrP+wI3AtTHt3oQd\n4vYCjgZ+Iym54fVFlNk7uoD7HLLDbbrZ4j6HbPH+mS/SjBwOBqaZ2SwzWwOMAoYWyQwF7onHo4Ej\n4vGxwCgz+8jM3gSmxfyQtDNwDHDnRl2B4ziOkzlplENvYE4iPDfGlZQxs3pgiaSeJdK+lUh7A3Ax\nYI0V7j6H7HCbbra4zyFbvH/mizRrK6lEXPEDvZxMyXhJXwIWmlmdpM+XkQNg9OjR3HnnnVRXBzdG\n9+7dGTBgwLqOVBiKerh1w912GwisN60UHpSlwlN6LAB6pZZPEyau11RJfsqkF1k6Y/5Gl9fc60lb\nftrraev2rJswjmUf79Lm/a+9hQu0ZHm1tbWMHDkSgOrqaqqqqqipqaG5yKzRF3ckHQJcZWZDYng4\nYGZ2TULmiSgzXlJHYL6ZVRXLSnoSGEEwQ30d+AjYGugGPGxm3ygu//rrr7dzzjmn2RforKe2tjaz\nt7OX5y1r0sJ7V/91ZmZyeclz6Yy6iqOHPNQzyzyvO6YfA3fqlirPppJl/8wbPXv2BGDx4sWtVubk\nyZOpqakp++JdiTQjh4lAP0l9gPnAqcBpRTKPAmcC44GTgLEx/hHgAUk3EMxJ/YAJZjYeuBxA0uHA\nD0spBsdx8kXHDuHFIA1V23Rmx223bOEaOS1FReVgZvWSLgCeIvgo7jKzqZKuBiaa2WPAXcB9kqYB\niwgKBDN7RdKDhBlJa4BhVmmoUoT7HLJjU30rays2R5/DkpX1qUcj1x3Tr0nKwftnvki1n4OZPQns\nURQ3InG8ijBltVTaXwC/aCTv54Dn0tTDcRzHaR1y/4W0f+eQHT6PPFv8O4ds8f6ZL3KvHBzHcZzW\nJ/fKwX0O2eE23WzZHH0OLYn3z3yRe+XgOI7jtD65Vw7uc8gOt+lmi/scssX7Z75INVvJ2XyYv3QV\nC5evrii3un5tK9TGcZy2IvfKwX0O2ZHGprtw+epUXz6PiMstbM64zyFb3OeQL3JvVnIcx3Fan9wr\nB/c5ZIfbdLPFfQ7Z4v0zX+ReOTiO4zitT+6Vg/scssNtutniPods8f6ZL3KvHBzHcZzWJ/fKwX0O\n2eE23Wxxn0O2eP/MF7mfyuo4Tvsk7d4Pvu9DPsm9cnCfQ3a4TTdb3OfQOGn3fijs++D9M1/k3qzk\nOI7jtD6plIOkIZJelfS6pEtLnO8saZSkaZLGSapOnLssxk+VNDjG7SxprKRXJE2RdGG5st3nkB1u\n080W9zlki/fPfFFROUjqANwKHAXsA5wmac8isXOBxWbWH7gRuDam3ZuwQ9xewNHAbyQJ+Aj4DzPb\nG/g0cH6JPB3HcZw2Is3I4WBgmpnNMrM1wChgaJHMUOCeeDwaOCIeHwuMMrOPzOxNYBpwsJm9bWZ1\nAGa2HJgK9C5VuPscssNtutniPods8f6ZL9Ioh97AnER4Lg0f5OtkzKweWCKpZ4m0bxWnlbQrsD8w\nvgn1dhzHcVqQNLOVVCLOUso0mlbSNoSRxkVxBNGAm266ia5du1JdHdwY3bt3Z8CAAeveMgp2Sg9X\nDidtuuXk6yaMY+mMt9a9FRfs6sVh4qqs5c4nw1N6LAB6pZZPE05b/pRJL7J0xvyNLq/U9SR9Dhtb\n/ubcnnUT3mHgcYNT9c/2Gi7Q0r/vkSNHAlBdXU1VVRU1NTU0F5kVP+eLBKRDgKvMbEgMDwfMzK5J\nyDwRZcZL6gjMN7OqYllJTwIjotwWwGPAE2Z2U7nyr7/+ejvnnHOafYHOempraysO3V+etyz1kt1p\npik2Rba95bl0Rl1F01Ie6pn3PK87ph8Dd+qWqn+2V3r27AnA4sWLW63MyZMnU1NTU+oFPRVpzEoT\ngX6S+kjqDJwKPFIk8yhwZjw+CRgbjx8BTo2zmfoC/YAJ8dz/AK80phjAfQ5Zsqn+8NoK9zlki/fP\nfFHRrGRm9ZIuAJ4iKJO7zGyqpKuBiWb2GHAXcJ+kacAiggLBzF6R9CDwCrAGGGZmJumzwNeAKZJe\nIpiaLjezJ1vgGjd70u7uBr7Dm+M4gVRfSMeH9h5FcSMSx6sIU1ZLpf0F8IuiuOeBjmnKrqurY9Cg\nQWlEnTIUdndLawZx0pGmPZ30bMpmpfaIfyHtOI7jNCD3ysF9Dtnhb7nZ4u2ZLT5qyBe5Vw6O4zhO\n65P7VVnd55AdbiPPFm/PbCgs7V03YRz7H/zpRmV9ee/WI/fKwXGcTZvC0t5LZ7zFtu82/o1NYXlv\np+XJvXJwn0N50k5RLUxP9bfcbPH2zBZvz3yRe+XglKcwRbUSPj3VcZymknuHtO/nkB2+/0C2eHtm\ni7dnvsi9cnAcx3Fan9ybldznkB1u080Wb89sSdOehZlNlfBZTRtP7pWD4zhOgcLMpkr4rKaNJ/dm\nJfc5ZIfbdLPF2zNbvD3zRe6Vg+M4jtP65F45uM8hO9xGni3entni7Zkv3OeQM3zvBcfZeNxxvfGk\nUg6ShgA3sn6zn2uKzncG7gUOAN4FTjGz2fHcZcA5wEeEvaKfSpNngc1tbaW0H7ZB0z9u87WAssXb\nM1uybE93XG88FZWDpA7ArUANMA+YKOnPZvZqQuxcYLGZ9Zd0CnAtYXvQvQmbAO0F7Az8VVJ/QCny\nBGD69HQPyrzT1KUuWoIV86b7wyxDvD2zxdszW+rq6qipqWl2+jQjh4OBaWY2C0DSKGAokHyQDwUK\nO8ONBm6Jx8cCo8zsI+DNuI3owQTlUClPAD744IOmXlMuycNSF/UfbhptmRe8PbPF2zNbXn755Y1K\nn0Y59AbmJMJzCQ/4kjJxz+klknrG+HEJubdinFLkmXvcP+A47Zu0vgnY/PwTaZSDSsRZSply8aVm\nSRXnCcDbb79dtmJmxlorXUgxq+vXslWnVNtWN8kE9OMxb6TKMw+L3616r3xbOk3H2zNb2qI90/om\nAH715X6pngvbdO7I8tX1Jc8VK6I8KxyZlXwmrxeQDgGuMrMhMTwcsKQDWdITUWa8pI7AfDOrKpaV\n9CTB/KRKeRb47ne/a0nT0sCBA316azOpq6vztssQb89s8fbcOOrq6jYwJXXt2pXbbrstzbtzSdIo\nh47AawTn8XxgAnCamU1NyAwD9jWzYZJOBY4zs4JD+gHgUwRz0tNAf8LIodE8HcdxnLajolkp+hAu\nAJ5i/bTTqZKuBiaa2WPAXcB90eG8CDg1pn1F0oPAK8AaYJgFbVQyzxa4PsdxHKcZVBw5OI7jOJsf\nuVk+Q9KJkv4lqV7SoKJzl0maJmmqpMGJ+CGSXpX0uqRLW7/W7QNJIyTNlTQ5/g1JnCvZtk7jeN/b\nOCS9KellSS9JmhDjPibpKUmvSRojqXtb1zOvSLpL0gJJ/0zElW0/STfH33mdpFSOndwoB2AK8FXg\nuWSkpL1Y/yHd0cBvFCh8nHcUsA9wmqQ9W7fK7Ypfmdmg+PcklG/btqxke8D7XiasBT5vZp80s8I0\n9uHAX81sD2AscFmb1S7/3E3of0lKtp+ko4Hdzaw/8B3g9jQF5EY5mNlrZjaNhjNThxI/pDOzN4HC\nh3TrPs4zszVA4UM6pzSlHvrl2tZpHO97G49o+PwZCtwTj+8BjmvVGrUjzKwWeK8ourj9hibi743p\nxgPdJfWqVEZulEMjFH+EV/iQrtTHeb1bsV7tjfPjkPLOxHCzXNs6jeN9b+MxYIykiZK+GeN6mdkC\nADN7G9i+zWrXPqkqar+qGN+s33mrrsoq6WkgqbFE6CQ/NrNHyyUrEdekD+k2BxprW+A3wE/MzCT9\nFLge+CbpPnB0GuLttvF8xszelrQ98JSk1/A2bCma1V9bVTmY2RebkWwusEsivDNhsT4B1SXiN0ua\n0La/AwqKuFzbOo0zF+97G0V8s8XM3pH0J4KpboGkXma2QNIOwMI2rWT7o1z7Net3nlezUlLTPUJY\n4bWzpL5AP8JHcxOBfpL6xCXDT42yThGxoxQ4HvhXPC7Xtk7jeN/bCCR1kbRNPO4KDCZMSHkEOCuK\nnQn8uU0q2H4QDZ+VZ8Xjs1jffo8A34B1K168XzA/NUZuNvuRdBxhNdePA49JqjOzo/1Duky4Nk5f\nWwu8SZix0NhHik4jlPswtI2r1Z7oBfxRkhGeQQ+Y2VOSJgEPSjoHmA2c1JaVzDOSRgKfB7aTNJuw\nLNEvgf9X3H5m9rikYyRNBz4Azk5Vhj8LHMdxnGLyalZyHMdx2hBXDo7jOE4DXDk4juM4DXDl4DiO\n4zTAlYPjOI7TAFcOjuM4TgNcOeQASafHLVSbk/Zfkj6XdZ3ygqRn47ztTRpJMyUd0cplDpb0cCK8\nTNKurVmHjUHS4ZLmVJZse+Ky+fdlkM/3JP0iizpVol0rh6x+UJLOlPT31kgbv6pdG5d9BsDMRhb2\n066Q9m5JP0nGmdm+Zva3ptW6ddhcHuztmJ8B6x40ZtYtrs7bnmhPH2plUdc7gK9L+ngGeTVKu1YO\nGVJYpK410hbkN/t9ExT2J3faAEkHAtua2cS2rkt7J/mi19KY2SrgceJyGC1dWLv8I6xPXk/4HHwp\n8KMYfwjwPGGt85eAwxNpzgJmRPkZwGnAnsCHhOUjlgGLy5SXOi1wDDAZWALMAkYk8pkV670s5vUp\nwjoyf0/I3AAsAN4H6oC9gW8Bq4GVMd2fo+xM4Ih43AG4HJgey54I9E7Znocm2m0W8I0Yv21s64Wx\nrB8n0pwJ/B24Dlgc2+WoeO6nwEfAiljfm2P8WmAY8DowI8Z9hrCm03vAeODTiTKeBc7ZiH5yCWGR\nsbnAubH83VLcpz5R9izCUgSLCMuOHAi8HK/3lqKyziEsRbIIeAKobqReZxCWMnkn3rPkfTwIeCG2\nx1uEZWW2iOduBf67KK9HgAvj8aXxWpcCU4EvlCn/SuCOorhk29wdy3os5jUO6Fsmr0JbfSvW9y3g\nPxLnRdiIZnq83lFAj8T5YwnrfS0mbFKzZ+LczJj237Fd7wI6x3OHA7MTsjsCowl9dQbwvUba/27g\nNsISKEtjP6tOnN8znlsU2/GkorS/Af5C+B0fUSL/XYH/i31rTLyH9ybOPwjMj/f4/4C9Y/yBwNtA\nh4TsCcBLifDpwDNZPEcb/e20dAEtWvnQcb6QCO8EvMv6B1RNDG8HdIk3ql881wvYKx6fCfytkXKa\nlBb4HLBPPN43doJjEz+keuLSJcV5EBYhmwh0i+E9COvcFzrlT0q0QeGhcjHhwVWo5wDgY/H4UeCS\nMte3S/yBnAx0BD4G7BfP3Qv8MbZBH+A14OxEvVcRHooCzgPeSuTb4MFOeIiMAboDW8ayFscO34Gw\niN3iRL2brRyAIQTFsCewFetfKAoPwEr3aS3hIdAZOJLwIvBw7E87ERT4YVH+OILC+wTrlfTzZeq1\nN+Gh8lmgE2EJ9dWJ+ziIsEppYeXhf7P+4X8QMDeR13bAcsKaZJ8gKLJCf6mm/AP9QeCHRXHJtrmb\n8Ns5IF7P/cDIMnkV2uqB2M77Eh7Qhev5PkHZ7Riv97ZCXrHOy4EjCH3vYsKmUwVlOBP4Z2zvHkAt\n8TdAQjnEtppEWKK+I+HhPB34Ypk63034TRfuwY3EFzRCX59NeDsXsD9Bqe2VSPsecEgMdy6R/wuE\nl6ZOwGGE31dSOZwVy+kE/IoNH/7/Ij7DYvhh4PuJ8CeBd7N+nja4hpYuoEUrn3gwxvAlwD1FMk8S\n3tK6EB46XwW2KpJJoxyalTbK3ABcn/gh1bPhm0FSOXwBeJUwolBRPpWUw6vAl5vRjsOBh0rEdyCM\nVPZIxH0bGJuo9+uJc1sTHhJVMVxOORyeCH8deLFI5gXWj1w2RjncBfwsEd6dxAMw5X3aIXH+XTZ8\ngxzN+of240SlmWi7D4BdSpRzJYkHbexfqyjxBhrPX5S8PwRlUROPzwceS1zf24SXoi0qtM1TwLdL\n3Jukcrgjce5o4JUyeRWUQ/9E3DXA7+LxK2z4ErcjQRl2AK4g7EZYOCfCyOdzif79raJ6TIvHSeXw\nKeDNEv36rjJ1vrvoHnQlWAB6E16SniuSvx24MpH294207S7x+rZOxD1AQjkUyfeI7Vd4IbwEuD8e\n94z9qFdCvh+wpjm/iab8bWo+hz7AyZIWx7/3CG8GO5rZCuAU4LvAfEmPStojTaZNTSvpYEljJS2U\n9D7BHJHKgWRmzxKG878G3pZ0e2F54xTsAryRUrY43YwS8R8nvNnMTsTNYsNdpN4uHJjZh/GwUn3n\nJo53inkmKS6jJHF2zdL4t3MJkZ3YcAesOSR8PZI+leI+JfcU+JAwWkiGC9faB7ip0PcI5ggrcx0b\n1Cv2r0WJevWPfWx+rNfPiup1L0GpEv/fF/OZQXhLv4qwtv9ISTuWKB/Cm2+3MucKvJ04XkHj99XY\n8L7OIlwnhLb5Y6JtCqsA96Lo/lt4+s1hw3Yrl2+SaqB30W//MtbvhlaK5D34gNAmO8X6HlKU1+ls\nuJlWY7OkdgLeS/weCvUGgo9C0i8lTY/3dyah/Qr3+H7gy5K6EBTV32zDJba7EUY9LUp7Vw5WFJ5D\n0M4949/HLMzAuBbAzJ42s8HADgTzyB1l8mlYUNPSjgT+RLD39wB+y/qHUpqybjWzAwmb1+9BGGqn\nSTuH8PbYVOYQ3kaKeZfwI+6TiOtDsCmnoVx9k/HzCCaAJNVpyoj3dtv4N7eEyHzCxibJfJM8QPn7\n1FTmAN8p6nvbmNmLZeq1bvOV+BDYLnH+NoKde/dYrx8X1et+YKik/Qgmsz8VTpjZKDM7jPX37Jdl\n6vtPgkknK8SGG8pUs35DmdnA0UVt09XM5keZPkV57cKGCiGZbx9Kb1QzB3ijqIzuZvaVRuqcvAfb\nEEyc82Je/1eU17ZmdkEibWO/xfnAxyRtnYhL9r2vAV8hjBR7EPr/ur0ZzGwewcdzPAnln2Avgvm4\nRWnvyuFtYLdE+H7gK3H+dgdJW8W50DtJqpL0lfhDXEOwc9bHdAuAnSV1KlVIM9JuQ3hzWCPpYMJb\nR4F3CEPIkg9xSQfGkccWhDfTlUVl7VYqXeRO4L8k9Yt5DZD0sUbkCzwA1Eg6UVJHST0lDTSztQTb\n9M8kbSOpD/ADGnbWclSqLwRzTH9Jp8ayTyF0/nLbxjaFB4GzJe0Z792VbPijbuw+QdMUxe3A5ZL2\nBpDUXdKJZWRHE94MPxP7zU+KyuoGLDWzFZL2JIxY12FmbxHs6/cRzE2rYpmfkPQFhQ2IVhP6Tz2l\neZywH0CWXClpa0n7EPYMGBXjfwv8XFJ1rOf2ko6N5x4EvhTrvYWkHxH6/LhEvudL6i2pJ2E0MIqG\nTACWSrok/u47StonzsoqxzHxHnQG/gsYH9v2MeATkr4e69Qp/i7TWhpmE+7P1THtoQRlUGAbghnx\nPYXNjn5BQ2VzH8G8tC/B55fkcMKEh5alpe1WLflHmOUwi+AP+I8YdxDB+7+I8HB6lPD2uEOMf4+i\nWTqms9QAAAJVSURBVBEE08mjMc3CEuU0KS1hdsGbhKHfI8DNbOiMuopgrlhMcDwmfQ5HEN4KlkaZ\n+4Autt7W+FJM93CMe4OGs5XeiGWPB3aK5x4HhjfSlp8FXmT9zJ0zbL099L5Yl1k0nK1U7IxPOjUP\nIYyyFgE3Fp9PpPkM4cf0HsEZn5ytNJaNm610KeFNbi7BbFRPnMFFeDMreZ8o7RuaTbSFx/C9wOWJ\n8NcIb+Tvx7a6s5F6nRFl3iE88JL38TDCyGEp8FzsL8Xt/LVYv2R9BsR7voQw6nuEhM+kRB3GAweV\nuXcb+LcomhlUlE/B5/BNwohvHglnN0HxfZ/gE1tCcDj/NHF+KMGP8h7Bx7RX4tzMeA//Tej3/0P0\n+xXXifA7HRnv9yKC76qcH+duwmSDpwiTA/4P6JM435+gJBbGe/RX1k/S2KBtyuS/K/C3eA/HFPWt\nroTR3tJ4fV+n6HdB8N8tAf6nKN+tCCOb7Zv7m0j755v9OJsN8S18CrClhVFRu0XSYcB9ZrbrRuTx\nReC7Znb8RtalD0G5dcq6XSXNBM41s7EZ53s3MMfM/jPLfLNEYee2byevXWEHwp3NbHhLl5+bbUId\npyVQ2H72L4Sh/DXAI5uAYuhEmMH0u43Jx8yeBp7OpFL+UWemSDoBWFusFM3s1taqQ3v3OThOJb5D\nMAtMI+6T3bbV2Tji6Oc9wsyZm9q4OklaygTR3vLdaCQ9S5it2KZ91c1KjuM4TgN85OA4juM0wJWD\n4ziO0wBXDo7jOE4DXDk4juM4DXDl4DiO4zTAlYPjOI7TgP8PiOZOcFmZt44AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10d3727f0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_test_stat(ert, \"variance: {} people per day\".format(variance_param))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "What this means is that if we had randomly taken 24 days and 80 days with a standard deviation of 100 people, there's almost no chance of observing a difference in number of people of 83."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "What if the standard deviation was larger? Luckily we can repeat our experiment by changing only this value."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.2327"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "variance_param = 500\n",
    "ert1 = ERTest(data, variance_param)\n",
    "ert1.PValue(iters=10000)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEdCAYAAADwwTuSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmYFcXVuN/DsKiAKOiAogMoRFFRgxtZlOgoolHI5waa\nGDXGxKiJJp8LGA1qNIr+4hJNNCZ+xI0QPjVREzYFQz6UbYKDE0EFRLZhEUEWkW04vz+q7tC09965\nM3fmzlRz3ue5z+2qPlVdp6u7T1edqmpRVQzDMAwjF5o1dgEMwzCMcDCjYRiGYeSMGQ3DMAwjZ8xo\nGIZhGDljRsMwDMPIGTMahmEYRs6Y0TB2QUSGicgHjV0Oo37w9TmvsctRH4jI5SKyrbHLsbtjRsOI\n8wDQp7ELUVdEpK+I7Ij9qkTkezG5TiIyWkTW+d+fRWT/mExzEblfRCpFZJOI/J+I9C6sRvVCUiZj\nKcnRJVjMaBgAiKOZqm5S1TWNXZ48UeBYoJP/HQA8n9opIgL8A+gClAJnAF8C/hrL5/8BVwBXAccD\nHwKvi0hxA5d/t0FEWjR2GYzaYUYjMETk+yLyqYi0isUPEZFFkfCTIjLfvyEvEJF7RKRlZP8wEZkn\nIheJyFxgC9Aj3p0hIl1F5EURWSYin4nIOyLyndix3xCRP4jIbSKyXEQ+EZGnRWSvmNwgESkTkc9F\nZLWI/ENE2kX2/1hE5vr974vIrSJSVMdTtVpVV0V+WyL7zsAZlW+rapmqzgAuBb4qIqf4srQBfggM\nUdV/qOocnAHZAlyd6aCRls45IjLd61IhIqfG5A4VkRdEZK2IrBGR8SJyVEzmbH++NovIShH5bfSc\nisgIEXlNRH4qIkt9/YwWkX2znRgROUNEpvhrY6mI/I+ItK8hzQ4R+Ykv80af7icxmdYi8kikLP8W\nkf+K7O/i87nE1/1G4K4sx/yl13u9iIwE9o3tz3ptiuvOWisie8TSDROR97Ppa2RBVe0X0A/YG/gM\nuDAW/x/gLr8twC9xb8clwDnAMmBYRH6Yz+cN4ASgO9DGx38QkTsKuMb/dwOuBbYCfSMybwBrgF/j\n3thPBz4B7ozIXOHT3Qoc7vP7MdDe778DWAgMwLUA+gMfxfK4A9hRw/npC+zAtQpWAm8C343J3AHM\nT5N2MXCr3z4VqAIOisk8A0zI4fjvA2cBhwF/BDYCHb1MMbAceAw4AugBPAJ8DHTwMkcD23CtncOA\nM4FFwNORY40A1gF/8/mcAnwAvBir52h9nubr/RrgEOA4YCLwzxrO6w5gtU/X3dfdNuDc2HUwCfgK\n0BX4PrAZONXv7+LzWQxc4sNdMhzvemAD8B1/vBuBtcDWXK9NYA/cdXhpJI346+zGxr6XQ/01egHs\nV4dKgz8Dr0bCx/sHXI8saW4A3o+EhwHbgc4xuV0eMhny+hvw+0j4DaA8JvM74M1IeBHwSIb89vQP\nsn6x+EuBtZHwtcC7NZTtS7iWwPFAb+Dn/sEVNT6/B6akSTsDeNRvX+zPafOYzP1ARZbjp4zG5ZG4\nIiIGEGe03oqlE2A+8BMffhaYFpMZ4Mt0sA+PANYDbSIyZ/jjH5KuPn1d/SqWb4lPc3QWvXYAf4rF\nPQ9M9tvfADYBbWMyTwEv+e2U0bg1h2t8Cf4lKBL3v0SMRo7X5iPAvyLhM/31sH9d7j37Kc0xQuQZ\n4G8isp+qrga+C0xX1Wi30lXAlbg3vtZAc9yDKcpKVV2W7UAisifuwXMOzjfQ0v/eiImWx8KVQD+f\nx/7AwcBrGQ5zJM5wvOjcDdUUAS1FpIOqfqKqvwV+m628qvoB7m07xSzfb/4zEblLVauypSc3R2tN\nMgpMi5SpSkRm4PQEZ9COF5ENsXR74N6qwbUcJsb2T8bV4RG4hyrAHFXdGJF5M5L+wzRlOwE4SUR+\nnKbMPYB3sug1LRZ+k53dS8cDrYDKWB22YNf6AJiZ5RiISFugMzA1tmsKMDAil8u1+XugQkQOV9X3\ncK2fV1T142xlMDJjRiNMxuOa3ZeIyG+BQcAvUjtF5EJc18fNwL9wb6MXAXfH8vksh2P9P+Bc4Ke4\nm/8z4EFcN1mUrbGw8kWfWaaHbUruAiDd8NB8HfNvAbcD+wMrcF1DpWnkOvr9eBlwjvSlGWRqg7BT\n/2bA67iWU9yQr4tsZzpf+Ri2ZsBwXEsmTm31iuv0Kc54xHWKXxs1XXep9DXpWeO1qapzRORN4Psi\nMhzXWju7hnyNLJjRCBBV3eEdg5fi+mf3Bv4SETkZmKWqj6QiRKRbHQ93MvC8qr7o8xFcF1DODxhV\n/VhEluK6Bv6RRuRdXJfBoao6vo7lzMZxwOe4Pnlwb8i3i8ihqroAQER64lpDU7zMv3EPuzNxXSwp\n3U8HnqjheIIbtvyeT1eEe8N/2u8vAy4DKnVXB32Ud3FdXVG+gevemROJ6ykibSKtja/hHrZzM+Rb\nBhypqulaITXRh111/2rkOGXAPsCe6gYN1BlVXS8iy3C6jIvs+npMNNdr8/fAwzifyFJVjbfgjNrQ\n2P1j9qvbD+iFe4DMAv43tu9anON1AM7ZeT3OyVoVkUnru4jH4/qR5+AeekcAf8C9UU6KyLwBPBnL\n5+fAh5HwlbiRR7fhHOFH+nKmHOG3+Xyvwd34R+BaUPdF8rgOmFvDebkBOA/nPP4S8BPcG2g0H8E9\n5KZ6vU7CdZlMieX1EM6Z/k1fnj/hWngdsxw/5dN4D+cIP9yfs41AJy9TjGu9jMU9CLv4/7uBPpH6\n3YobXHAYbmDAIiJ+BZxP41PgJX8+T8E54F/KUp/f8PXwa+AYf330xznrW2XRK+UIv5adjvCtwICI\nzHiv90CcY7q3r7Mr/f6UT+OrOVzfN+BayClH+H/jWpxRR3iN16aXa4W7/j/HjYZr9Ps35F+jF8B+\neVSeMxhVwDmx+ObA4/4m/xR4DvcwrovROMg/3DbgR2D5mzNqNCZRg9HwcRcDb/ub92PgVWDvyP4r\nvE6bcA/nqcAPY2XbXsM5uRH39rsR92Y5E/heGrmOuNbZOn+ORgL7xWSKgPtw/plNwP8BX67h+H1T\ndYIzTJ/jRradFpM7GNdFtNLLLMT5qrpEZPr78n/u5R7Dvcmn9o8AJgA/82XcCIwG9s1Wz7g3+Ale\n9w24Vs2DQLMseu3AGeC/4ozwMuD6mEwr4FfAAlzLsRIYA3zD7+/iz00uRkNwRnSVL+No3MtP1GjU\neG1GZB/EGcuMBt9+uf3En9CsiEh/XPOuGfCUqg6P7W/pL/jjcA+qQaq62O8bCnwPN1LnelWdIG6O\nwb9wTqvmwAuqeqeXH4G78dbhmtmXq2o255xhNBlEpC/OiB6sqpUNfKwRuNFv/RryOP5YO4DvqOrI\nhj5WQyAif8GNhDu/scsSOjX6NESkGe4NpxT35jBTRF5WNxIhxZXAGlXtISKDcMMSB4vIETgHbE/c\nW8HrItJDVbeIyKmqusn3974pImPVTbIC+G9Vjc/ONYxQiDuCjUZCRPbBdT9+CzdHxciTXGaEnwjM\nU9VFqroNGEVk2JtnIDudfC+ws3IGAKNUdbuqfoQbGXMigKpu8jKtcMYr2uSxmepGyOQyuik0QtXp\nbVzX1nBVfbMmYaNmchk91ZmdY8LBOfBOzCSjbkz6Or8sQXys9TIfl2rB/Bs4FPitqkbHbt8tIrfj\nxqkP8cbKMJo8qjoZ5wspxLGuKMRx/LEKolN9o6p1HTVoZCCXN/p0Te34W0cmmYxpVXWHqn4Z1211\nku/KAmckeuJGRHQAbsmhjIZhGEYByKWlsRS3zECKg3C+jShLcKNBKr2Pop2qrvVj8w/OllbdmOx/\n4kaKzFHVlT5+m3f0/Xe6Qg0YMEA3b95Mp06dAGjdujXdu3fn2GOPBaC83E1QDjX8wgsvJEqf3Um/\n1HZTKY/pt3vrV15ezvjxbvpTp06daN26NY8//nid/W41jp7yRuB9nCN8OW59notVdW5E5hrgKFW9\nRkQGA99S1ZQj/HmcI6ozbhmJHrgWxDZVXeeXAhiPG0c/RkQ6qeoKP1HnQeBzVb01Xq7vfve7+sgj\nj8SjE8N9993HkCFDGrsYDUaS9WtI3dq3d4vRrlnTeKvXJ7nuIPn6XX/99TzzzDN1Nho1tjS8j+I6\n3Lju1JDbuSJyJzBTVf+OmzH7rLgltT8BBvu0c0RkNG4CzjbgGlVVETkAeNr7NZoBf1HVMf6Qz4vI\nfriurXIyLEO9YkVdVnIIh8WLFzd2ERqUJOuXZN3A9NvdyWkZEVUdh5uVGo0bFtneghtamy7tvcC9\nsbgK3GzRdPLp1gQyDMMwmgDBDm0988wzG7sIDcoll1zS2EVoUJKsX5J1A9MvdI455pi80uc0I7wp\nMnHiRO3dO8TPNRtG3WkKPg0jbGbNmkVpaWnD+TSaKuXl5aQzGqrKqlWrqKqq6bMJTZt169bRrl27\nmgUDJVT9ioqKKC4uJvbNiF2YMmUKX/96fEHW5GD67d4EazQysWrVKtq2bctee+1Vs3AT5sADD2zs\nIjQooer32WefMWfhUrbvsQ8AxW1acsDerbKmWb5+C6s2xj8psSu55GMYTYFgjUZqPHKcqqqq4A2G\n0XRp3bo1i9Yu5eG33ac5Hji7+xce9vG31FUbt3LTmPlZ802XT1Ml6W/hSdcvX4J1hBtGY9EsS9eU\nYSSdYI1GdNamYTQlpkyZUrNQwJh+uzfBGo2kce211/KrX/2qsYtRK4YPH87VV6ede2kYRkIJ1mhk\n8mkYhSXTKKI333yTo446ql6OMWDAAJ577rmCpB0+fDg/+tGPcpYvagazKzfs8mt7yDG7hLdW7ahL\n0ZssSe/zT7p++RKsI9yoPVVVVRQVFWaFa1XNOiw1KazbXMWdry/MKjPsdFud20gOwbY0QvRpfPDB\nBwwYMIBu3brxta99jXHjxu2y/5NPPuG8886jpKSEAQMGsHTp0up9t956K4cddhhdu3bllFNO4b33\n3IcTt27dyu23387RRx9Nz549ufHGG9myZQuw823/N7/5DT179uTHP/4xffr04bXXXqvOt6qqih49\nelBRUQHAzJkz6d+/P926daNv3768+ebO79YsXryYc889ly5dunD++ednnGC2adMmBg0axIoVKygp\nKaGkpISVK1eiqjz88MMcd9xx9OjRgyuvvJJ169YBsGXLFq6++mq6d+9Ot27dOP3001m9ejX33HMP\nU6dO5ZZbbqGkpCTtQnK1TTt06FB69epFly5dKC0tZdq0aQBMnDiRhx56iL/+9a+UlJTQt2/f2lWw\nZ/2C8K7N2pD0Pv+k65cvu1VLIzWbtj6o7Yzc7du3c8kll3DppZfy0ksvMXXqVL797W/zxhtvcOih\nhwJuufC//OUvHHfccfziF7/gBz/4AWPGjGHSpElMnz6dsrIy2rZty7x586onxg0bNowlS5YwZcoU\nioqK+MEPfsADDzzAbbfdBrh5K+vWreOdd95hx44dPProo7zwwgucccYZgHtQdujQgV69elFZWcnF\nF1/M73//e0pLS5k8eTKXXXYZM2bMoH379lx11VWcdNJJvPTSS5SVlTFo0CC++c1vfkHXvfbai9Gj\nR3P11VdXGyOAxx9/nLFjx/KPf/yDDh06MGTIEG688Ub+8Ic/8Oc//5kNGzbw7rvv0rJlSyoqKthj\njz34+c9/zvTp07nooov4zne+k/bc1jbtcccdx5AhQ2jbti1PPPEEV1xxBbNnz6a0tJSf/vSnfPTR\nRzz++OO1ql/D2F0ItqURmk+jrKyMTZs2cf3119O8eXNOPvlkzjzzTF588cVqmX79+tGnTx9atGjB\nbbfdRllZGZWVlbRo0YKNGzfy/vvvo6r06NGD4uJiAJ577jnuuece9t57b1q3bs3111+/S55FRUUM\nGTKEFi1a0KpVK84//3zGjh3L5s2bAXjxxRe54IILAGe0+vXrR2mpWzOyb9++HHvssbz22mssXbqU\n8vJyhg4dSosWLfjKV75C//79a3UOnn76aW677TY6depEixYtuOmmm3jllVfYsWMHLVq0YM2aNSxY\nsAAR4eijj6ZNmzY55VvbtBdccAHt2rWjWbNmXHPNNWzZsoX587PPo6gNex8a1rVZW5Le5590/fJl\nt2ppNOZ6PcuXL//CLOiDDz6Y5cuXV4c7d+5cvd26dWv22WcfVqxYwcknn8z3v/99br75ZpYtW8Y3\nv/lN7rrrLjZv3symTZs49dRTq9Pt2LGD6HpiHTp0oEWLFtXhbt26cdhhhzFu3DjOPPNMxo4dy623\nus+VLFmyhL/97W/V3WaqSlVVFaeccgorVqxgn332Yc8999yl/JWV8e9xZWbp0qVceumlNGvWrDr/\nFi1asGrVKgYNGkRlZSVXXnkl69ev58ILL+T222/PyQczePDgWqV97LHHeO6551i5ciUAGzdu5JNP\nPslZD8PYnQm2pRGaT+OAAw74wgN26dKlHHDAAdXhZcuWVW9v3LiRtWvXVn+Z8KqrrmLSpEm89dZb\nzJ8/n0cffZQOHTqw11578dZbb/Hhhx/y4Ycf8tFHH7Fo0aLqfNI5o8877zxefPFFxowZw+GHH06X\nLl0AZ7QGDRpUndfChQtZvHgxP/nJT+jUqROffvopn3/++S7lz0S643bu3JnRo0fvkv/SpUvp1KkT\nzZs356abbmLq1KmMHz+e8ePHM2rUqIx5RSkqKso57bRp03j00Uf505/+xMKFC1m4cCFt27atNrT1\n4bw3n0bYJF2/fAnWaITGcccdx1577cVvfvMbtm/fzpQpUxg/fjznn39+tcxrr73G9OnT2bp1K7/6\n1a84/vjjOfDAA3n77bf597//zfbt29ljjz1o1aoVRUVFiAiXXnopt956K6tXu2UtKisrmTRpUtay\nnHfeebzxxhuMGDGiumsK4MILL2T8+PFMmjSJHTt2sHnzZt58802WL1/OQQcdxLHHHst9993Htm3b\nmDZt2hcc+VH2339/1q5dy/r166vjLr/8cu6+++5qY7N69WrGjh0LuBt1zpw57Nixg9atW9OiRYvq\nlsL++++/iyGMU5u0GzZsoHnz5rRv356tW7dy//33s3Hjxur9xcXFLF68mFBXfzaMhiZYoxGaT6NF\nixY8//zzvPbaa3Tv3p2bb76ZJ554otoJLiJccMEFDB8+nO7du1NRUcGTTz4JuAfdDTfcwCGHHMKX\nv/xlOnTowHXXXQfAHXfcwSGHHEK/fv3o2rUr559/PgsWLMhalo4dO3LCCSdQVlbGf/3Xf1XHd+7c\nmeeee46HHnqIHj16cMwxx/DYY4+xY4ebZ/Dkk09SVlbGoYceygMPPMDFF1+c8Rg9evTgvPPOo3fv\n3hxyyCGsXLmSq6++mrPOOovzzz+fLl260L9/f2bNmgXAypUrueKKK+jatStf/epX+frXv85FF7nv\nev3whz/k5Zdf5tBDD2Xo0KFfOFZt0p5++umUlpZywgkn8OUvf5k999xzl27BgQMHoqoceuihnHba\nadkrNQPm0wibpOuXL4n7nkZlZWWwK6gaYTBu1gc8OOszwM3ByGWeRk0yD5zdnWMObFvjse17Gka+\n5Ps9jWBbGqH5NIzdB/NphE3S9cuXYI2GYRiGUXiCNRqh+TSM3QfzaYRN0vXLl2CNhmEYhlF4gjUa\n5tMwGosdNQweMZ9G2CRdv3zJyWiISH8ReU9EPhCRW9Lsbykio0RknohMFZGSyL6hPn6uiPTzca1E\nZLqIvC0iFSIyLCLfVUSmicj7IvJnEanVrPWioiI2bdpUmySGkTPrNnzG4k+3NXYxDKPRqPGBLCLN\ngMeAUqASmCkiL6vqexGxK4E1qtpDRAYB9wODReQI4CKgJ3AQ8LqI9FDVLSJyqqpuEpEi4E0RGauq\nM4DhwK9V9X9F5HGf9+/j5crk0yguLmbVqlV8+umnuZ8Fw0jDxq1VfLhm5wz4Haos/nQbf124Fcg8\nYtF8GmGTdP3yJZe3+BOBeaq6CEBERgEDgajRGAikWgsvAI/67QHAKFXdDnwkIvN8ftNVNdUcaOXL\nkWrznwakZo09DdxBGqORCRGhY8eOuYobRkZmV27gwVmr0+xJ/ndCDCMTuXRPdQaWRMJLfVxaGVWt\nAtaJSPs0aZel0opIMxF5G1gBvKaqM0WkA7BWVVOfOlsKpJ2pl3SfRtL7VZOsn/k0wibp+uVLLi2N\ndK9VcU9gJpmMab1h+LKI7A38zXdlrUqTJq3XcfLkyZSVlVFS4twn7dq1o1evXtVNy1TFhxpOfYei\nqZRnd9RvwepNwP7ATkOQ6nqqNgz+q3yZ9ucaLp8xlQ377VVj+VI0hfNj4TDCU6ZMYeTIkQCUlJRQ\nXFxc/fmDulDjMiIi0ge4Q1X7+/AQQFV1eERmrJeZ7n0Uy1W1OC4rIuOAYao6PXaMXwAbVfVBEfkY\n6KiqO/yxh6nqWfFyZVpGxDDqi9mVG7hpTPbvbNgyIkZoFGIZkZlAdxHpIiItgcHAKzGZV4HL/PaF\nQGqZ1VdwDvGWItIN6A7MEJH9RKQdgIjsCZwOzPVpJvk88Hm+XHu1DMMwjIagRqPhfRTXAROAd3GO\n7bkicqeInOPFngL2847uG4AhPu0cYDQwBxgDXKOuaXMA8IaIlAPTgfGqOtbnNQT4mYh8ALT3eX8B\n82mETZL1M59G2CRdv3zJaQ6Eqo4DDovFDYtsb8ENrU2X9l7g3lhcBZC2b0lVFwIn5VIuwzAMo7AE\n+7nXpK89lfSx4k1Bv+Xrt7Bq49aM+7dW7ci4Lxs2TyNskq5fvgRrNAwjX1Zt3JrV0T3Mj4wyDGMn\ntvZUEyXp/apJ1s98GmGTdP3yxVoahtEEKGrmhvhmo7hNywKVxjAyE6zRMJ9G2CRZv7r4NNZtrspp\nLkdTIMl1B8nXL1+C7Z4yDMMwCk+wRsN8GmGTZP3MpxE2SdcvX4I1GoZhGEbhCdZomE8jbJKsn83T\nCJuk65cvwRoNwzAMo/AEazTMpxE2SdbPfBphk3T98iVYo2EYhmEUnmCNhvk0wibJ+plPI2ySrl++\nBGs0DMMwjMITrNEwn0bYJFk/82mETdL1y5dgjYZhGIZReII1GubTCJsk62c+jbBJun75EqzRMAzD\nMApPsEbDfBphk2T9zKcRNknXL1+CNRqGYRhG4QnWaJhPI2ySrJ/5NMIm6frlS7BGwzAMwyg8wRoN\n82mETZL1M59G2CRdv3zJyWiISH8ReU9EPhCRW9Lsbykio0RknohMFZGSyL6hPn6uiPTzcQeJyCQR\nmSMiFSLyk4j8MBFZKiKz/K9/fShqGIZh5E+N3wgXkWbAY0ApUAnMFJGXVfW9iNiVwBpV7SEig4D7\ngcEicgRwEdATOAh4XUR6ANuBn6lquYi0Af4tIhMieT6oqg9mK5f5NMImyfqZTyNskq5fvuTS0jgR\nmKeqi1R1GzAKGBiTGQg87bdfAE7z2wOAUaq6XVU/AuYBJ6rqClUtB1DVjcBcoHMkP6mLMoZhGEbD\nkovR6AwsiYSXsusDfhcZVa0C1olI+zRpl8XTikhX4FhgeiT6WhEpF5E/iki7dIUyn0bYJFk/82mE\nTdL1y5cau6dI/9avOcpkTeu7pl4ArvctDoDfAXepqorI3cCDuO6vXZg8eTJlZWWUlDj3Sbt27ejV\nq1d10zJV8aGGKyoqmlR5kqjfgtWbgP2BnQ/6VNfS+gXlVOyzEuiYcT8Ap3fLuj/XcEXZNNYvWJ5V\nvnzGx6Ro7PqzcDjhKVOmMHLkSABKSkooLi6mtLSUuiKq8ed/TECkD3CHqvb34SGAqurwiMxYLzNd\nRIqA5apaHJcVkXHAMC/XHPg7MFZVH8lw7C7Aq6p6dHzfxIkTtXfv3nXR2TAAmF25gZvGzM+4f9jp\n3bjz9YVZ8yikzANnd+fUo7oAsGbNmqyyhpGJWbNmUVpaWmcXQC7dUzOB7iLSRURaAoOBV2IyrwKX\n+e0LgUl++xWcQ7yliHQDugMz/L7/AebEDYaIdIoEzwP+k6syhmEYRsNSo9HwPorrgAnAuzjH9lwR\nuVNEzvFiTwH7icg84AZgiE87BxgNzAHGANf4bqevAd8GThORt2NDa+8XkXdEpBzoC/w0XbnMpxE2\nSdbPfBphk3T98iUXnwaqOg44LBY3LLK9BTe0Nl3ae4F7Y3FvAkUZ5L+bS5kMwzCMwhPsjHCbpxE2\nSdbP5mmETdL1y5dgjYZhGIZReII1GubTCJsk62c+jbBJun75kpNPwzBCY/n6LazauDWrzNaqHQUq\njWEkh2CNhvk0wqah9Vu1cWvWORjg5kY0BObTCJuk65cvwRoNw9jdKIp0Js+u3JBWprhNSw7Yu1WB\nSmTsjphPo4mS9H7VJOvXUD6NdZurqrdvGjM/7a+mLrn6IMl1B8nXL1+CNRqGYRhG4QnWaJhPI2yS\nrJ/5NMIm6frlS7BGwzAMwyg8wRoN82mETZL1s3kaYZN0/fIlWKNhGIZhFJ5gjYb5NMImyfqZTyNs\nkq5fvgRrNAzDMIzCE6zRMJ9G2CRZP/NphE3S9cuXYI2GYRiGUXiCNRrm0wibJOtnPo2wSbp++RKs\n0TAMwzAKT7BGw3waYZNk/cynETZJ1y9fgjUahmEYRuEJ1miYTyNskqyf+TTCJun65UuwRsMwDMMo\nPMEaDfNphE2S9TOfRtgkXb98ycloiEh/EXlPRD4QkVvS7G8pIqNEZJ6ITBWRksi+oT5+roj083EH\nicgkEZkjIhUi8pOI/L4iMkFE3heR8SLSrj4UNQzDMPKnRqMhIs2Ax4AzgSOBi0Xk8JjYlcAaVe0B\nPAzc79MeAVwE9ATOAn4nIgJsB36mqkcAXwGujeQ5BHhdVQ8DJgFD05XLfBphk2T9zKcRNknXL19y\naWmcCMxT1UWqug0YBQyMyQwEnvbbLwCn+e0BwChV3a6qHwHzgBNVdYWqlgOo6kZgLtA5TV5PA9+q\ntVaGYRhGg5CL0egMLImEl7LzAf8FGVWtAtaJSPs0aZfF04pIV+BYYJqPKlbVlT6vFcD+6QplPo2w\nSbJ+jenTKGoGsys3ZP0tX78lr2Mkue4g+frlS/McZCRNnOYokzWtiLTBtUyuV9XPcihLNZMnT6as\nrIySEuc5qVftAAAgAElEQVQ+adeuHb169apuWqYqPtRwRUVFkypPaPqVz5jK+gXLqruKUg/yaLhi\nn5VAxzrvB+D0bln35xquKJvG+gXLcygvWfNbt7kbd76+MOvxHji7OwvemVmv59vCTTc8ZcoURo4c\nCUBJSQnFxcWUlpZSV0Q1/vyPCYj0Ae5Q1f4+PARQVR0ekRnrZaaLSBGwXFWL47IiMg4Y5uWaA38H\nxqrqI5G85gLfUNWVItIJeENVe8bLNXHiRO3du3edFTeSzezKDdw0Zn5WmWGnuwdsXfc3hsy5vQ8B\n4Pj7J9Y5nwfO7s4xB7bNKmMkl1mzZlFaWpruhT4ncumemgl0F5EuItISGAy8EpN5FbjMb1+Ic2Dj\n5Qb70VXdgO7ADL/vf4A5UYMRSXO5374MeDlHXQzDMIwGpkaj4X0U1wETgHdxju25InKniJzjxZ4C\n9hORecANuBFQqOocYDQwBxgDXKOqKiJfA74NnCYib4vILBHp7/MaDpwhIu8DpwP3pSuX+TTCJsn6\n2TyNsEm6fvmSi08DVR0HHBaLGxbZ3oIbWpsu7b3AvbG4N4GiDPJrcMbCMAzDaGIEOyPc5mmETZL1\ns3kaYZN0/fIlWKNhGIZhFJ5gjYb5NMImyfqZTyNskq5fvgRrNAzDMIzCE6zRMJ9G2CRZP/NphE3S\n9cuXYI2GYRiGUXiCNRrm0wibfPRbvn5Ljesrba3aUY+lrR3m0wibpOuXLznN0zCMpsSqjVtzWiLE\nMIz6J9iWhvk0wibJ+plPI2ySrl++BGs0DMMwjMITrNEwn0bYJFk/82mETdL1y5dgjYZhGIZReII1\nGubTCJsk62c+jbBJun75EqzRMAzDMApPsEbDfBphk2T9zKcRNknXL1+CNRqGYRhG4QnWaJhPI2yS\nrJ/5NMIm6frlS7BGwzAMwyg8wRoN82mETZL1M59G2CRdv3wJ1mgYhmEYhSdYo2E+jbBJsn7m0wib\npOuXL8EaDcMwDKPwBGs0zKcRNknWz3waYZN0/fIlJ6MhIv1F5D0R+UBEbkmzv6WIjBKReSIyVURK\nIvuG+vi5ItIvEv+UiKwUkXdieQ0TkaUiMsv/+uejoGEYhlF/1Gg0RKQZ8BhwJnAkcLGIHB4TuxJY\no6o9gIeB+33aI4CLgJ7AWcDvRER8mhE+z3Q8qKq9/W9cOgHzaYRNkvUzn0bYJF2/fMmlpXEiME9V\nF6nqNmAUMDAmMxB42m+/AJzmtwcAo1R1u6p+BMzz+aGqU4C1GY4pGeINwzCMRiQXo9EZWBIJL/Vx\naWVUtQpYJyLt06RdliZtOq4VkXIR+aOItEsnYD6NsEmyfubTCJuk65cvuXwjPN1bv+Yok0vaOL8D\n7lJVFZG7gQdx3V+7MHnyZMrKyigpce6Tdu3a0atXr+qmZariQw1XVFQ0qfI0Nf1SD+ZUV1A8XFE2\njfULlmfcv35BORX7rAQ61nk/AP5b5DWVp/7KS72Up7Hr38KFC0+ZMoWRI0cCUFJSQnFxMaWlpdQV\nUc3+DBeRPsAdqtrfh4cAqqrDIzJjvcx0ESkClqtqcVxWRMYBw1R1ug93AV5V1aMzHDvj/okTJ2rv\n3r3roLIROrMrN3DTmPlZZYad3o07X1+Yl0x95FHfMuf2PgSA4++fWOd8Hji7O8cc2DarjJFcZs2a\nRWlpaZ1dALl0T80EuotIFxFpCQwGXonJvApc5rcvBCb57VeAwX50VTegOzAjkk6ItUZEpFMkeB7w\nn1wUMQzDMBqeGo2G91FcB0wA3sU5tueKyJ0ico4XewrYT0TmATcAQ3zaOcBoYA4wBrhGfdNGREYC\nbwFfEpHFInKFz+t+EXlHRMqBvsBP05XLfBphk2T9zKcRNknXL19y8Wngh70eFosbFtneghtamy7t\nvcC9aeIvySD/3VzKZBiGYRSeYGeE2zyNsEmyfjZPI2ySrl++5NTSMIxCsXz9FlZt3JpVZmvVjgKV\nxjCMOMEajfLycpI8emrKlCmJfuPJpN+qjVtzGhnVlFm/oDzRrY3d9do0HMEaDcMw6kZRMzdsORvF\nbVpywN6tClQiIySCNRrm0wibJOvX1FsZ6zZX5TSXI5PRSHLdQfL1y5dgHeGGYRhG4QnWaNg8jbBJ\nsn42TyNskq5fvgRrNAzDMIzCE6zRMJ9G2CRZv6bu08iXJNcdJF+/fAnWaBiGYRiFJ1ijYT6NsEmy\nfubTCJuk65cvwRoNwzAMo/AEazTMpxE2SdbPfBphk3T98iVYo2EYhmEUnmCNhvk0wibJ+plPI2yS\nrl++BGs0DMMwjMITrNEwn0bYJFk/82mETdL1y5dgjYZhGIZReII1GubTCJsk62c+jbBJun75EqzR\nMAzDMApPsEbDfBphk2T9zKcRNknXL1+CNRqGYRhG4QnWaJhPI2ySrJ/5NMIm6frlS05GQ0T6i8h7\nIvKBiNySZn9LERklIvNEZKqIlET2DfXxc0WkXyT+KRFZKSLvxPLaV0QmiMj7IjJeRNrlo6BhGIZR\nf9RoNESkGfAYcCZwJHCxiBweE7sSWKOqPYCHgft92iOAi4CewFnA70REfJoRPs84Q4DXVfUwYBIw\nNF25zKcRNknWz3waYZN0/fIll5bGicA8VV2kqtuAUcDAmMxA4Gm//QJwmt8eAIxS1e2q+hEwz+eH\nqk4B1qY5XjSvp4Fv5aaKYRiG0dDkYjQ6A0si4aU+Lq2MqlYB60SkfZq0y9KkjVOsqit9XiuA/dMJ\nmU8jbJKsn/k0wibp+uVL8xxkJE2c5iiTS9o6MXnyZMrKyigpce6Tdu3a0atXr+qmZariQw1XVFQ0\nqfIUSr+2hxwD7Hzwprp6ahuuKJvG+gXLs8pX7LMS6Fjn/QCc3q3A5aUg5SmfMZUN++3VZK4XC9c9\nPGXKFEaOHAlASUkJxcXFlJaWUldENfszXET6AHeoan8fHgKoqg6PyIz1MtNFpAhYrqrFcVkRGQcM\nU9XpPtwFeFVVj47kNRf4hqquFJFOwBuq2jNerokTJ2rv3r3rrLjRNJlduYGbxszPKjPs9G7c+frC\nBpcp1HFqI3Nu70MAOP7+iQ16rAfP6U7Vjsz7i9u05IC9W2XNw2iazJo1i9LS0nQv9DmRS0tjJtDd\nP+CXA4OBi2MyrwKXAdOBC3EObIBXgOdF5CFct1R3YEYknfDF1sgrwOXAcJ/nyznqYhhGPbFuc1VW\nw/LA2d3NaOym1OjT8D6K64AJwLs4x/ZcEblTRM7xYk8B+4nIPOAG3AgoVHUOMBqYA4wBrlHftBGR\nkcBbwJdEZLGIXOHzGg6cISLvA6cD96Url/k0wibJ+iXdp5F0/ZJ8bdYHubQ0UNVxwGGxuGGR7S24\nobXp0t4L3Jsm/pIM8mtwxsIwDMNoYgQ7I9zmaYRNkvVL+jyNpOuX5GuzPgjWaBiGYRiFJ1ijYT6N\nsEmyfknv80+6fkm+NuuDYI2GYRiGUXiCNRrm0wibJOuX9D7/pOuX5GuzPgjWaBiGYRiFJ1ijYT6N\n8Fi+fguzKzcwu3IDT/9tQvV29Lc12zTkQEh6n3/S9UvivVef5DRPwzDqg1Ubt1YvEbJ+wTL2Xv3F\n5UKG+bWTDMNomgTb0jCfRtgkuV88ybpB8vVL+r2XL8EaDcMwDKPwBGs0zKcRNknuF0+ybpB8/ZJ+\n7+VLsEbDMAzDKDzBGg3zaYRNkvvFk6wbJF+/pN97+RKs0TAMwzAKT7BGw3waYZPkfvEk6wbJ1y/p\n916+BGs0DMMwjMITrNEwn0bYJLlfPMm6QfL1S/q9ly/BGg3DMAyj8ARrNMynETZJ7hdPsm6QfP2S\nfu/lS7BGwzAMwyg8wRoN82mETZL7xZOsGyRfv6Tfe/kSrNEwDMMwCk+wRsN8Gk2L6LcyMv2i38pI\ncr94knWD5OsX2r1XaHL6noaI9AcexhmZp1R1eGx/S+AZ4DhgNTBIVRf7fUOB7wHbgetVdUK2PEVk\nBNAXWAcocLmqvpOnnkYDE/1WRibsWxmGET41Gg0RaQY8BpQClcBMEXlZVd+LiF0JrFHVHiIyCLgf\nGCwiRwAXAT2Bg4DXRaQHIDXk+d+q+tds5TKfRtgkuV88ybpB8vVL+r2XL7m0NE4E5qnqIgARGQUM\nBKJGYyAwzG+/ADzqtwcAo1R1O/CRiMzz+UkNeQbbbWYYuwNFzWB25YasMsVtWnLA3q0KVCKjUORi\nNDoDSyLhpbgHf1oZVa0SkXUi0t7HT43ILfNxUkOed4vI7cBEYIiqbosXqry8nN69e+dQ/DCZMmVK\not941i8oT+wba5J1A6ffus3duPP1hVnlHji7e5BGI+n3Xr7kYjQkTZzmKJMpPl1LIpXnEFVdKSIt\ngD8AtwB3x4UnT55MWVkZJSUlALRr145evXpVV3bKmRVquKKiokmVp6Zw+Yyp7rvf/mGZcpZGwxX7\nrAQ6ArCpcv4X9gPg/R7p0tcmXFE2jfULludcnrrsb5zy0iTKk0t5y2d8zDHf6gc0/vW5O4enTJnC\nyJEjASgpKaG4uJjS0lLqiqjGn/8xAZE+wB2q2t+HhwAadYaLyFgvM11EioDlqloclxWRcbhuLKkp\nTx/fF+ffGBAv18SJEzXJLY3QmF25ISdHeE1vp01JpimVJSVzbu9DADj+/omNWp5c8njg7O4cc2Db\nrDJG4Zk1axalpaXpXuhzIhffwUygu4h08aOkBgOvxGReBS7z2xcCk/z2KziHeEsR6QZ0B2Zky1NE\nOvl/Ab4F/KeuyhmGYRj1S41GQ1WrgOuACcC7OMf2XBG5U0TO8WJPAft5R/cNwBCfdg4wGpgDjAGu\nUUfaPH1ez4vIbGA20IE0XVNg8zRCJ8lj/ZOsGyRfv6Tfe/mS0zwNVR0HHBaLGxbZ3oIbWpsu7b3A\nvbnk6ePr3tlmGIZhNCjBDm21eRphk+TRRUnWDZKvX9LvvXwJ1mgYhmEYhSdYo2E+jbBJcr94knWD\n5OuX9HsvX4I1GoZhGEbhCdZomE8jbJLcL55k3SD5+iX93suXYI2GYRiGUXiCNRrm0wibJPeLJ1k3\nSL5+Sb/38iWneRrG7s3y9VtYtXFrVpnoB5YMw0guwRoN82kUjob4wFKS+8WTrBskX7+mdO81RYI1\nGoZhNG3smxvJJFijYd/TCJskf3MiybqB92nk0LJct7kqyG9uJP3ey5dgHeGGYRhG4QnWaJhPI2yS\n/CaeZN0g+fol/d7Ll2CNhmEYhlF4zKfRRClkv2pNQ2obYjhtkvv9k6wb5O7TyIWm6Cw3n0Z2gjUa\nRv1R05Da2g6nNYxcCdVZvjsTbPeU+TTCJslv4knWDZKvX9LvvXwJ1mgYhmEYhSdYo2FrT4VNktcv\nSrJukHz9kn7v5UuwRsMwDMMoPMEaDfNphE2S+8WTrBskX7+k33v5YqOnEo6tUGuETlMclrs7k5PR\nEJH+wMO4lslTqjo8tr8l8AxwHLAaGKSqi/2+ocD3gO3A9ao6IVueItIVGAXsC8wCLlXV7fEy2TyN\n3GiIFWrrgyTPZUiyblC/8zRyodDDcm2eRnZq7J4SkWbAY8CZwJHAxSJyeEzsSmCNqvbAGYL7fdoj\ngIuAnsBZwO/EkS3P4cCvVfUw4FOf9xeYPz/7gzB0KioqGrsIDcqmyuTWX5J1g+Trl/R7L99BRLm0\nNE4E5qnqIgARGQUMBN6LyAwEhvntF4BH/fYAYJRvKXwkIvN8fpIlz9OAi336p4E7gN/HC/XZZ5/l\npmGgrFu3rkaZkLueqj5Pbv0lWTdIvn653HshM3v27LzS52I0OgNLIuGluAd/WhlVrRKRdSLS3sdP\njcgt83GSLk8R6QCsVdUdkfgDc9Rlt6Opdj0ZRqExv0fhyMVoSJo4zVEmU3y6brGUfDxN/FgArFix\nIl10g7J5WxUti7L36G3bobRqnl0mlxbCu/M+rPEmaKqtiFzYsrbw9VcokqwbNE39cvF7PHhO9xrv\nu+I2LVm8eHF9Fi1xiGraZ/JOAZE+wB2q2t+HhwAadYaLyFgvM11EioDlqloclxWRcbhuLMmUp4h8\nDHRU1R3+2MNU9ax4uX70ox9ptIvqmGOOSdQw3PLy8kTpEyfJ+iVZNzD9QqO8vHyXLqnWrVvz+OOP\np3uhz4lcjEYR8D5QCiwHZgAXq+rciMw1wFGqeo2IDAa+paqDvSP8eeAkXLfUa0APXEsjnudgVX1P\nRP4CvKSqfxGRx4HZqvpEXRU0DMMw6o8au6e8j+I6YAI7h8fOFZE7gZmq+nfgKeBZ7+j+BBjs084R\nkdHAHGAbcI06K5Uuz5RjfQgwSkR+Cbzt8zYMwzCaADW2NAzDMAwjRRDLiIjIj0XkPRGpEJH7IvFD\nRWSeiMwVkX6R+P5e/gMRuaVxSl07RORGEdnhR52l4n7j9SsXkWMj8Zd53d4Xke82TolzQ0Tu9/VT\nLiIvisjekX2Jqb8UIZcdQEQOEpFJIjLH328/8fH7isgEf82NF5F2kTRpr9OmjIg0E5FZIvKKD3cV\nkWlevz+LSHMf31JERnn9popISeOWvGZEpJ2I/K+/r94VkZPqtf5UtUn/gG/gurGa+/B+/r8nrvuq\nOdAVmI9zsDfz212AFkA5cHhj61GDjgcB44CFQHsfdxbwD799EjDNb+8LLADaAfukthtbhyy6nQ40\n89v3Aff67SOSUn8RXYMte0SHTsCxfrsNzvd4OG7S7c0+/hbgPr+d9jpt6j/gp8BzwCs+/BfgQr/9\nOPBDv/0j4Hd+exBu3lmjl78G3f4EXOG3m/tnRb3VXwgtjR/hFNwOoKqrffxA/MRBVf0ISE0crJ6M\nqKrbcEuSDCx8sWvFQ8BNsbiBuKVZUNXpQDsR6YibRT9BVdep6qc4g9q/kIWtDar6uu6cdzMNZyAh\nMvEzAfWXIuSyA6CqK1S13G9vBObi6mwgbrIt/j+lV6brtMkiIgcBZwN/jESfBrzot58GvuW3o3q/\ngBu802QRkbbAyao6AsDfX+uox/oLwWh8CTjFNx3fEJHjfHx80mFq4mC6yYidC1LSOiAi5wJLVDW+\ndkEmPTLpHQLfA8b47UTUX4yQy/4F/Dpwx+KMfUdVXQnOsADFXizE6zH1kqYAGSYVp3TYZeIy8Gm0\nC7kJcgiwWkRG+O63J0VkL+qx/prEKrci8hoQtW6Cq9DbcGXcR1X7iMgJwP/iTkxtJw42GjXodytw\nRrpkacLZJkw2Gln0+7mqvuplfg5sU9U/R2TiNMn6qwVNrm7qioi0wb1ZX6+qG0Ukkx5B6Swi3wRW\nqmq5iHwjFU3mScWZ7sOmSnOgN3CtqpaJyEO4Ean1Vn9NwmioarqHJgAicjXwkpebKSJV/s1gKRB1\nSh0EVOJOQrr4RiOTfiJyFK4/f7aICK6ss0TkRJx+B0fEU3osxfl5ovFv1H+pcydb/YFz3OO6A06L\nRGfSr8nVXy3IdE0GhXcCvwA8q6ov++iVItJRVVeKSCdglY/PVI9Nla8BA0TkbGBPoC1ukdV2ItLM\ntzaiOqT0qxQ3Z21vVV3bCOXOlaW4nosyH34RZzTqr/4a22mTg1PnB8CdfvtLwCK/nXKktgS6sdOR\nWsROZ2RLnDOyZ2PrkaOuC4F9/fbZ7HRQ9SG9Izy1vU9jlz2LTv2Bd4EOsfgk1l+wZY/p8QzwYCxu\nOHCL3x7CTkdq2us0hB/Ql10d4YP89uPA1X77GnY6wgcThiN8MvAlvz3M11291V+jK5jDCWgBPAtU\nAGVA38i+of4mnQv0i8T3x436mAcMaWwdaqHrh/jRUz78mNdvNtA7En+51+0D4LuNXe4adJoHLMJ9\nG2VW6gZMYv2FXnZf/q8BVd7gve3rrD/QHnjd6/YakReVTNdpU//FjEY3YLq/p/4CtPDxrYDRvj6n\nAV0bu9w56HUMMNPX4Uu4F8x6qz+b3GcYhmHkTAijpwzDMIwmghkNwzAMI2fMaBiGYRg5Y0bDMAzD\nyBkzGoZhGEbOmNEwDMMwcsaMRhNGRC4R94ncuqT9j4icUt9lair4dci+19jlaGhEZKGInFazZL0e\ns5+IvBQJb/DrUAWBiPQVkSU1SzY+IjJMRJ6th3x+LCL31keZaiKRRqO+bjT/3Yr/K0RaEeki7nsa\n1XWiqiPVf0e9hrQjROSuaJyqHqWq/6pdqQvD7vLAD5h7gOoHkKq2VbcScUiENAGtPsr6JPAdEdmv\nHvLKSiKNRj2Sz+JktU2bbUHC3Qq/xo/RCIjI8bj1lWY2dllCJ/oC2NCo6hbcCtIN/1G2xp7y3gBT\n6J/BLYPwGbAeuNHH9wHeBNbilkfoG0lzOW4Np/X+/2Lch2c+x33bfAOwJsPxck6LW+dlFrAOt7TG\nsEg+i3y5N/i8TgIuA/4vIvMQsBL4FLdEwBHAVcBWYLNP97KXXQic5reb4VbTne+PPRPonOP5/Hrk\nvC3CL1sC7O3P9Sp/rJ9H0lwG/B/wALDGn5cz/b67ge3AJl/e3/j4Hbh1fj4AFvi4rwIz/LGnA1+J\nHOMN4Ht5XCc3s3MByCv98Q/JoZ66eNnLgcXAJ8APgeNxyzCsAR6NHet7wBwvOxYoyVKuS4GPgI99\nnUXr8QTgLX8+lgGPsvPjZI8B/y+W1yvAT/z2LV7X9bhlW07NcPzbgSdjcdFzM8If6+8+r6lAtwx5\npc7VVb68y4CfRfYLbh2k+V7fUey6vMUA4D/+nE4i8kErf16G4NY1+wR4Cmjp9/UFFkdkD8AtwLjK\nX4s/znL+R+DWnprg9XsjWl+4e3uCP+Zc/IebIml/B/wDdx+flib/rsA//bU13tfhM5H9o4Hlvo7/\nCRzh448HVuA/aObjzgfejoQvASbWx3M0673T0AdojJ+/oE6NhA8EVrPzwVXqwx2AvXwFdvf7OuIX\nmcM9/P6V5Ti1SgucAhzpt4/yF8eAyA1Whf9uezwPoB/uYd/Whw/DrZGfuljvSnMOUg+bm3APtFQ5\ne7FzYcRX8V/0SqPfwf7GuQi3GN++wNF+3zPAX/056IJb0+aKSLm34B6WAlwNLIvk+4UHPu7hMh63\nTk4rf6w1/kZohlssbk2k3HU2Gri1lCpxD4A92PmikXow1lRPO3APh5a4LxN+jlvjp4O/1lbiPoQD\n7mM+H+AW20wZ7zczlOsI3MPma7g1136NeyFI1WNv3IeeUisBv8tOo3ACsDSSVwdgI7CfP/biyPVS\nQuYH/Wjgv2Nx0XMzAnfvHOf1eQ4YmSGv1Ll63p/no3AP7pQ+N+CM4AFe38dTefkyb8StjFyEu4bn\nsdNILgTe8ed7H2AK/h4gYjT8uSoDfu7z6YozUmdkKPMI3D2dqoOH8S9uuGt9Me5tXnDfGvmYnff8\nCNzDvo8Pt0yT/1u4l6kWwMm4+ytqNC73x2kBPMiuRuE/+GeYD78E3BAJfxlYXd/P0y/o0NAHaIwf\nkQemD98MPB2TGYd7q9sL9zD6L2CPmEwuRqNOab3MQ8CvIzdYFbu+SUSNxqnAe7gWiMTyqclovAec\nU4fzOAR4MU18M1zL5rBI3A+ASZFyfxDZtyfu4VHsw5mMRt9I+DvEVtz0N9x3M+VRC72eAu6JhA8l\n8mDMsZ46RfavZtc3zhfY+TAfgzemkXP3GXBwmuPcTuQB7K+vLaR5Y/X7r4/WD86IlPrta4G/R/Rb\ngXtZal7DuZkA/CBN3USNxpORfWcBczLklTIaPSJxw4E/+O057PpydwDOSDbDfWtmVGSf4FpKp0Su\n76ti5Zjnt6NG4yTgozTX9VMZyjwiVgetcT0GnXEvT5Nj8k8At0fS/inLuT3Y67dnJO55IkYjJr+P\nP3+pF8Wbgef8dnt/HXWMyHfHfbMm72dott/u4tPoAlwkImv8by3uTeIAVd2E+/bvj4DlIvKqiByW\nS6a1TSsiJ4rIJBFZJSKf4ro1cnJcqeobuG6B3wIrROQJ/6GcXDgYt4JubTkY15yPsx/uTWhxJG4R\nu37xa0VqQ1U/95s1lXdpZPtAn2eU+DHS4kf7rPe/g9KIHMiuXytbQsSXJCIn5VBPqyLbn+NaF9Fw\nStcuwCOpaw/XraEZ9NilXP76+iRSrh7+Glvuy3VPrFzP4Iwt/v9Zn88C3Fv9HbjvKowUkQPSHB/c\nm3LbDPtSrIhsbyJ7vSq71usinJ7gzs1fI+dmDu4B3ZFY/at7Ki5h1/OWKd8oJUDn2L0/lJ1frktH\ntA4+w52TA315+8TyuoRdP0CWbdTWgbgvBH4eiavWUUSaich9IjLf1+9C3PlL1fFzwDn+S3wX4V4o\no9ddW1wrqUFJqtHQWHgJzpq397991Y0IuR9AVV9T1X5AJ1w3y5MZ8vnigWqXdiTwN5w/YR/g9+x8\nWOVyrMdU9XjgSFz3VOq74jWlXYJ726wtS3BvL3FW427uLpG4Lrg+61zIVN5ofCWuKyFKSS7H8HW7\nt/8tTSOynJ3fKk/lG+V5MtdTbVkC/DB27bVR1WkZylX9QRz/cOgQ2f84rh/9UF+un8fK9RwwUESO\nxnW9/S21Q1VHqerJ7Kyz+zKU9x1c11B9Iez6kZ8Sdn7kZzFwVuzctFbV5V6mSyyvg9nVUETz7UL6\njwctAT6MHaOdqp6bpczROmiD6yqt9Hn9M5bX3qp6XSRttntxObCviOwZiYtee98GzsW1LPfBXf/V\nXxVU1UqcD+k8Ii8FEXriuqEblKQajRW4T8KmeA44148/byYie/ix3AeKSLGInOtv0G24ftQqn24l\ncJCItEh3kDqkbYN709jmv853SWTfx7imaNqHu4gc71sqzXFvsptjxzokXTrPH4Ffikh3n1cvEdk3\ni3yK54FSEblARIpEpL2IHKPu62ajgXtEpI2IdAF+yhcv4kzUVF5w3To9RGSwP/Yg3E3xao7HyMZo\n4AoROdzX3e3serNnqyeonQF5ArhVRI4AEJF2InJBBtkXcG+SX/XXzV2xY7UF1qvqJhE5HNfCrUZV\nl+H675/FdVtt8cf8koicKiItcd0jn7Pz2okzhl2/DFkf3C4ie4rIkcAVOIc3OGP8KxEp8eXcX0QG\n+PGEBmIAAALoSURBVH2jgW/6cjcXkRtx1/zUSL7Xikhncd/sHhrJN8oMYL2I3Ozv+yIROdKPEsvE\n2b4OWgK/BKb7c/t34Esi8h1fphb+vsy1Z2Ixrn7u9Gm/jjMSKdrguiPXikhr3LDnuBF6FtdNdRTO\npxilL26gRcPS0P1fjfHDjbpYhPM3/MzHnYAbjfAJ7qH1Ku5ts5OPX0tslAauC+ZVn2ZVmuPUKi1u\ntMNHuCbkK8Bv2NUJdgeu22MNzuEZ9WmchnuLWO9lngX20p19mW/7dC/5uA/54uipD/2xpwMH+n1j\nyPKxIFw33jR2jiS6VHf2tz7ry7KIL46eig8CiDpT++BaZZ8AD8f3R9J8FXeTrcUNAoiOnppEfqOn\nbsG9+S3FdT9V4UeU4d7k0tYT6X1Pi/F97T78DHBrJPxt3Bv8p/5c/TFLuS71Mh/jHoTRejwZ19JY\nj/s62x1pzvO3ffmi5enl63wdrpX4ChGfTJoyTAdOyFB3u/jPiI1UiuWT8ml8H9dCrCTiZMcZxBtw\nPrd1OEf33ZH9A3F+mrU4H1bPyL6Fvg7fxV33/4P3K8bLhLtPR/r6/gTnG8vkJxqBG+QwATco4Z9A\nl8j+HjjjscrX0evsHByyy7nJkH9X4F++DsfHrq3WuNbheq/fd4jdFzj/4Drgf2L57oFrCe1f13si\n1599hMnY7fFv7RVAK3WtqGARkZNx3/bumkceZwA/UtXz8ixLF5zRa1Hf51VEFgJXquqkes53BO4b\n27+oz3zrExGZjxusMCkSdx1wkKoOaejjN2/oAxhGU0REvoUbT98GN6LnlQQYjBa4EVV/yCcfVX0N\n90nQ+mC3n6xan4jI+cCOuLFU1ccKVYak+jQMoyZ+iOtemIfzR13TuMXJD99aWosbyfNIIxcnSkN1\nZYSWb96IyBu40ZONeq1a95RhGIaRM9bSMAzDMHLGjIZhGIaRM2Y0DMMwjJwxo2EYhmHkjBkNwzAM\nI2fMaBiGYRg58/8BkQk26ec/F34AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10d38ab38>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_test_stat(ert1, \"variance: {} people per day\".format(variance_param))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here, we see that the observed difference would be likely to be observed by chance. Having a difference of 83 people per day is likely to be normal over 24/80 days if numbers fluctuate by 500 people every day. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The point of all of this is: if we don't know by how much numbers fluctuate every day and how (here I used a Gaussian assumption but really, I have no idea), we can't say if this is statistically significant. I wonder why the authors have not mentioned this basic statistic of their sample in the article."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Conclusions and last plot "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "To conclude, let's do a parametric plot of P-values, depending on the unknown parameter, the standard deviation of the number of people per day at the ER in Queensland, Australia."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "pvalues = []\n",
    "variances = [10, 50, 100, 150, 200, 250, 300, 350, 400, 500, 750, 1000]\n",
    "for variance_param in variances:\n",
    "    ert2 = ERTest(data, variance_param)\n",
    "    pvalues.append(ert2.PValue(iters=1000))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>variance (people per day)</th>\n",
       "      <th>p-value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>10.0</td>\n",
       "      <td>0.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>50.0</td>\n",
       "      <td>0.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>100.0</td>\n",
       "      <td>0.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>150.0</td>\n",
       "      <td>0.006</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>200.0</td>\n",
       "      <td>0.040</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>250.0</td>\n",
       "      <td>0.088</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>300.0</td>\n",
       "      <td>0.115</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>350.0</td>\n",
       "      <td>0.140</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>400.0</td>\n",
       "      <td>0.178</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>500.0</td>\n",
       "      <td>0.245</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>750.0</td>\n",
       "      <td>0.300</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>1000.0</td>\n",
       "      <td>0.391</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    variance (people per day)  p-value\n",
       "0                        10.0    0.000\n",
       "1                        50.0    0.000\n",
       "2                       100.0    0.000\n",
       "3                       150.0    0.006\n",
       "4                       200.0    0.040\n",
       "5                       250.0    0.088\n",
       "6                       300.0    0.115\n",
       "7                       350.0    0.140\n",
       "8                       400.0    0.178\n",
       "9                       500.0    0.245\n",
       "10                      750.0    0.300\n",
       "11                     1000.0    0.391"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df2 = pd.DataFrame(data= [variances, pvalues], index=['variance (people per day)', 'p-value']).transpose()\n",
    "df2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x111051b38>"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAESCAYAAAD9gqKNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VPXV+PHPSSAhEAghEJZAwuqGKKKyuKEGFR9tqWvR\np1W0Kta1Vfu4Pa2lfeHSSrXWn0qLSq1StNGn2kVk0aJRUBDQiCAkQFYWSUIgELLMnN8f9yaZCQmZ\nQCYzmTnv1ysv567zncN1ztzzvfd+RVUxxhhjYkLdAGOMMeHBEoIxxhjAEoIxxhiXJQRjjDGAJQRj\njDEuSwjGGGOAABOCiEwVkY0isklE7j/MeleKiFdExvnMe1BENovIBhG5sD0abYwxpv11aW0FEYkB\nngUygRJglYi8raobm6yXCNwJrPSZdzxwNXA8MBhYKiKj1G5+MMaYsBPIGcJ4YLOq5qtqLbAQmNbM\ner8GngCqfeZNAxaqap2qbgM2u/szxhgTZgJJCGlAoc90kTuvgYiMBQar6r9b2ba46bbGGGPCQ6sl\nI0CamddQ8hERAZ4Crm/rtsYYY8JHIAmhCEj3mR6M05dQrycwGviPmxwGAO+IyHcD2BaAH//4x5qX\nl8eAAQMA6NGjByNHjmTs2LEArFu3DiCip3Nzc7nyyivDpj2hnM7Kyoq6f3+Lh8UjkOl169bx3nvv\nATBgwAB69OjB888/39wP7yMirfXvikgs8A1Op/J24DPgGlXd0ML6HwD3qOpaETkBeA2YgFMqWgIc\n0qm8bNkyHTdu3KE7iyKPP/44DzzwQKibERYsFv4sHv4sHo3uvvtuXnnllXZLCK2eIaiqR0TuABbj\n9Dm8qKobRGQWsEpV/9l0E9xSkap+LSJvAF8DtcBtdoVR8woKCkLdhLBhsfBn8fBn8YD8ggLmzJ3f\nbE3+aARSMkJVFwHHNpn3SAvrnt9k+jHgsSNtoDHGmEb5BQXc+PCTMGk6J37zervu2+5UDhPXXntt\nqJsQNiwW/iwe/qI9Hr99/mWYNJ3YuAROPvnkdt23JYQwcdZZZ4W6CWHDYuHP4uEvmuOxsqCCj7ft\nITYuAWjseG4vAZWMQqmyspKKigqcC5giV0VFBUlJSR32fqpKUlISiYmJHfaegcrOzo7q/+mbsnj4\ni8Z4bN9XzfMrilhZsJdaL3hqqhqSQnsK64RQWloKwKBBgyI+IQwaNKhD309VKSsro7q6mpSUlA59\nb2NMYKrrvLzx5U5e/2InNR6le9cY7rrphyx85WU8k64BurXr+4V1Qqiuru7wL8poISKkpKRQUnLI\nbSEhF22//lpj8fAXLfH4tKCC51YUsX1fDQCZI5O5aXwaKd27ctGxfZkzdz5wabu+Z1gnBGOMiTbb\n91XzwopiVhRUADA0uRt3nDGEkwY2lncz0tN5ZvYvWLNmTbu+t3Uqm7CTnZ0d6iaEFYuHv0iNR02d\nl1fXbOfmrA2sKKige9cYZk5I47nLjvNLBsFkZwhh7uOPP2bmzJl89dVXoW6KMSZIPi2o4PmVRZTs\ndcpD549I5uYJTnmoI1lC6AQivUO9qWipEQfK4uEvkuLRtDyUkdyNO88YzEkDe4akPZ02IdTfur27\nspq+ifHcO3MGGenprW/Yzvswxpi2qqnz8kbOLhau29Fw9dAPxw1k2uh+dIkJ3Q/ATtmHUH/r9vq0\nKewecyXr06Zw48NPkt+GZ5y0xz7Gjh3L008/zaRJkxgxYgR33nknNTU1h6z3+9//nhkzZvjNe+CB\nB3jwwQcBWLBgARMnTiQ9PZ1TTz2V+fPnt/ieKSkpbNu2rWH69ttv59FHH22Yfu+995g8eTLDhg3j\n4osv5uuvvw7484SLSK0RHymLh7/OHo/PCiu45a0NvPL5dmo8ynkjknnxyhO4YkxqSJMBdNIzhDlz\n5zfcug0QG5eAZ9J0LvnZUwy64LqA9lGy5BX6Tz50H3PmzueZ2b8IuC1ZWVm89dZbdO/enenTp/Pk\nk0/y0EMP+a1zxRVX8OSTT1JZWUliYiJer5d33nmHV199FYB+/frxxhtvkJ6ezooVK7jqqqs49dRT\nGTNmzCHvd7jy0RdffMFdd93FwoULGTt2LG+88QbXXnstq1atomvXjq1FGmP87dhXzfMri1mR75aH\nenfjjjMGc/Kg0JSHmtMpzxB2V1YfcpdebFwCqt6A96HqbXYfpZXVLWzRvJtvvpmBAweSlJTEPffc\nw1tvvXXIOoMHD+akk07i3/92BpRbvnw53bt3p/6R3xdccAHpbqlq0qRJnHfeeaxYsaKFdrf8sNi/\n/OUvzJgxg1NOOQUR4fvf/z7x8fGsXr26TZ8p1CKpRtweLB7+Ols8auq8vLZ2BzdlbWBFfgUJXWO4\nZUIaz19+XFglA+ikZwh9E+PZ2eTWbU9NFZmjUnjmplMC2sddW//B+mb2kZIY36a2+N44N2TIELZv\n3873v//9hi/0p556iiuuuIIrrriCN998k6uvvpo333yTK664omG7JUuW8Nvf/pa8vDy8Xi8HDx7k\nhBNOaFM7AAoLC3n99df505/+BDjJo66uju3bt7d5X8aYo/dZoXNzWf3VQ+eNSOaW8Wmk9AjPM/ZO\neYZw78wZsGIhnpoqwPkiZ8VCZ34H7gOguLi44XVhYSEDBw7k9ddfp6CggIKCgoYv/mnTpvHxxx9T\nUlLCv/71r4bR0Wpqarjhhhu466672Lx5M1u3bmXKlCktngl0796dAwcONEzv2rWr4XVaWhr33HMP\nW7ZsYcuWLWzdupXCwkIuv/zyNn2mUOvsNeL2ZvHw1xnisWNfNb9csoX/fW8LJXtryOjdjd/+10ge\nPG9o2CYD6KQJISM9nZdm38fo4qX0y8lidPFSXpp9X5uuEGqPfQC8+OKLlJSUUF5eztNPP81ll13W\n7HopKSmcccYZ3HHHHQwdOpRRo0YBTkKoqakhJSWFmJgYlixZwgcffNDi+40ZM4Y333wTr9fL0qVL\n+eSTTxqWXXfddbz88st8/vnnAOzfv58lS5awf//+Nn0mY8yRqS8P3Zy1gU/qy0PjB4Vleag5nbJk\nBI23bod6H1deeSVXXHEFO3fu5L/+67+49957D7vubbfdxqxZsxrmJSYm8vjjj3PDDTdQU1PD1KlT\nufjii1vcx6OPPsptt93GvHnzuOSSS7jkkksaltVf9XT//fezZcsWEhISmDBhAmecccZRfcaO1tlq\nxMFm8fAXrvFYVbiX/7eiiJK9Tj/keSOSuXn8IPr2iAtxywLX6pjKHaGlMZVLSkrC+uF2Y8eO5Zln\nnuGcc84JdVOOWLjH2Jhwt3NfDS+sLOJj9+qhdPfqobEdcEawZs0aMjMz2+1a1U5ZMjKRrTPUiDuS\nxcNfuMSjxuNlwdod3JT1NR+75aGbxw/ihcuP65BkEAwBlYxEZCrwNE4CeVFVn2iyfCZwO+AB9gG3\nqOpGEckANgAb3VVXqupt7dX4UIu2R0oYYxyri/by/z4potgtD507vDe3TEjrVOWh5rSaEEQkBngW\nyARKgFUi8raqbvRZ7TVVneuu/x3gKaC+EJ6rqofWgyLA2rVrQ92EiBSuNeJQsXj4C2U8dlU65aHs\nbY3lodvPGMwpnfSMoKlAzhDGA5tVNR9ARBYC02j81Y+qVvqsnwj43iFmP6ONMZ1ajcfLmzm7WLB2\nB9UepVuXGH4wbgCXje5H19jIqbwH8knSgEKf6SJ3nh8RuU1EcoHHgbt8Fg0Vkc9F5AMRsZ86plXh\nUiMOFxYPfx0dj9VFe5n55kZeXr2dao8yeXhvXrzqeK4+qX9EJQMI7AyhuV/4h1yapKrPAc+JyHTg\n58AMYDuQrqrlIjIO+LuInNDkjKJF8fHxlJaW0qdPH6vXt7P6MZXj49t2Z7Yx0cIpDxWTvW0PAEOS\n4rnjjCGckhYZ5aHmBJIQigDfu7UG4/QltOR14AUAVa0BatzXa0QkDzgG8Bv3LSsri3nz5jU8zycp\nKYkxY8Zw1llnUVlZyfr164mNjSUpKQmAioqKhvVs+simPR4PQ4cOJTExseEXV31tNtTT9fPCpT2h\nnrZ4dGw8PvjwQz7csofVZFBd5+Xgti+5YFQfHrj8ErrGxoT082dnZ7NgwQIA0tPTSU1NJTMzk/bS\n6n0IIhILfIPTqbwd+Ay4RlU3+KwzUlVz3dffAX6uquNFpC9QpqpeERkOLAfGqOoe3/do6T4EY4zp\nSE2vHpo8rDe3TEyjX5hePdTh9yGoqge4A1gMrAcWquoGEZklIpe6q90hIl+JyBrgJ8D17vxzgC9F\nZC3wBjCzaTIwDqsTN7JY+LN4+AtGPHZV1vCrpVt5aFEexXurGZIUzxMXj+ThzGFhmwyCIaD7EFR1\nEXBsk3mP+Lz+SQvbvQUc+jxoY4wJA7UeL29+tYvX1u6kus7rXD10ygAuOzGyrh4KVFg/usIYY4Jl\nddFenltRRFGFUx46Z5hzc1lqYuc5I2jvklGnfbidMcYciV2VNcz9tJiPtjrV68FJ8dxxxmDGpfUK\ncctCL/rOicKU1YkbWSz8WTz8HWk8aj1eFn6xgx9lbeCjrXuI7xLDj04fxNzLj7Nk4LIzBGNMxFtT\nvJdnP+nc5aGOYH0IxpiItauyhj9+WsyHPuWh2ycN5tTBkXFGYH0IxhjTilqPl7e++pZX1+6gus5L\nfJcY/vuU/lx+YipxUXj1UKAsMmHC6sSNLBb+LB7+WovHmuK9zHxrIy+uKqG6zsvZw3rz4pXHM/3k\nAZYMWmFnCMaYiPDt/hr+uLKY5T7lodsmDea0CCkPdQTrQzDGdGq1Hi//55aHDkZZecj6EIwxxrW2\neB/PflJIoXv10FlDe3PrRLt66EhFdvrsRKxO3Mhi4c/i4S87O5tv99cwe9lW7n83l8KKagYnxfPo\n1BH8YsowSwZHwc4QjDGdRq3Hywd55fwmd4NTHooVrj1lAFeMifzyUEewhBAmbNzcRhYLfxYPR0N5\nqHYw4OWsoUncOnGwnRG0I0sIxpiwtnu/8+yh5Vucq4fSesVz+xl29VAw2DlWmLA6cSOLhb9ojUet\nx8sbX+7kxr9tYPmWPcTHCjecNpDrBuy2ZBAkdoZgjAk7a0v28f8+KaJgz0EAzhqaxMwJg+nfM47s\n7NwQty5y2X0IxpiwsXu/8+yh/7jloUG9nGcPnT7EzgiaY/chGGMiTp1X+b+vdvHq2h1U1TpXD10z\ndgBXnmRXD3Uki3SYiNY6cXMsFv4iPR7rSvbx47c28qfPSqiq9XJmRhLzrjyBa09p/tlDkR6PUAro\nDEFEpgJP4ySQF1X1iSbLZwK3Ax5gH3CLqm50lz0I3AjUAXer6uL2a74xprMq3V/LHz8r5oO8csAp\nD902KY3xQ5JC3LLo1WofgojEAJuATKAEWAVMr//Cd9dJVNVK9/V3gNtU9WIROQF4DTgdGAwsBUZp\nkze1PgRjIlt+QQFz5s5nd2U1KT3iGT35Et7d0cW/PDQmlbguVrRoi1D0IYwHNqtqPoCILASmAQ0J\noT4ZuBIBr/v6u8BCVa0DtonIZnd/n7ZD240xnUB+QQE3PvwkTJpObFwCO2uqWPr8H0k98zLOG3sM\nt05MY0DP+FA30xBYH0IaUOgzXeTO8yMit4lILvA4cFcL2xY3t62xuqgvi4W/zh6POS/Mb0gGALFx\nCQy68HrSij/mlxcMb3My6OzxCGeBnCE0dzpySJ1JVZ8DnhOR6cDPgRmBbpuVlcW8efNIT08HICkp\niTFjxjTcsl9/AETydE5OTli1J5TTOTk5YdWeUE931ngMG3MaS3PLWfL5BpJ0CL1GjAVgb946AES9\nURWP9pjOzs5mwYIFAKSnp5OamkpmZibtJZA+hInAL1V1qjv9AKBNO5Z91hegXFV7N11XRBYBj6iq\nX8nI+hCMiQwVB+tYvqWcpZvL2PjtAQBKlrxC/8lXNZwhAHhqqhhdvJRnZv8iVE2NCKHoQ1gFjBSR\nDGA7MB24xncFERmpqvW3D16K0wkN8A7wmog8hVMqGgl81h4NN8aEh5o6LysLK1i2uZzPCivwuL8x\nE7rGcNbQ3pxw0k3M+cMLeCZdQ2xcAp6aKlixkHtn3xfahptDtJoQVNUjIncAi2m87HSDiMwCVqnq\nP4E7RGQKUAOUA9e7234tIm8AXwO1OFcfhf7W6DCUnZ1tT7V0WSz8hWM8vKqs37mfpZvL+GjrHipr\nPADECJw2uCdTRvZhUkYSCV1jAThx9s+YM3c+pZXVpCTGc+/s+8hwS8RtFY7xiBQB3YegqouAY5vM\ne8Tn9U8Os+1jwGNH2kBjTPgoqjjI0s1lLMstZ2dlTcP8kSkJZI7sw3kjkunTvesh22Wkp1t5qBOw\nZxkZYw5rT1Uty7fsYWluGd+4/QIAfXt0JXNEMpmj+jA0OeEwezDBYs8yMsYEXXWdl5UFFSzdXMbq\nor1+/QJnD+1N5qg+nDwwkRhpt+8iEwYsIYQJq4s2slj466h4eFX5akclSzeX8+HWcg7UOpeFxgic\nPrgXU0YlMymjN91CfDexHR/BYwnBmChXsOcgyzaX8X6ef7/AqL4JTBnZh3OHJ5PcTL+AiTzWh2BM\nFCqvquU/eeUsyy1n0+7GfoF+PbqSObIPU0b2IT25WwhbaAJhfQjGmCNSXeflk/wK3s8tY1XRXrzu\nb8HuXWM4e1hvpozswxjrF4hqlhDChNVFG1ks/B1NPLyqfLm9kmW5zv0Cvv0CE4b0ItO9XyC+Ez1l\n1I6P4LGEYEwEyi+vYmluOe/nlvHt/tqG+cf2607myD5MHt6b5ATrFzD+rA/BmAhRfqCWD9znCOWW\nVjXM758Yx/kjk8kc2Yf03tYvEEmsD8EY0+BgnZcV+XtYurmcz4sb+wV6xMVyzrDeZI7sw4kDeli/\ngAmIJYQwYXXRRhYLf03j4fE6/QJLc8vI3raHKrdfIFZgYnovpozsw8T0pIgdfcyOj+CxhGBMJ7G1\nrIpluWW8n1vO7gON/QLH+fQL9LZ+AXMULCGECfvF08hi0aj0QC07ko7hx/+3kTyffoEBPePIHNmH\nzJHJDE6Krn4BOz6CxxKCMWGmqtbDJ/kVLMstY03xvoZ+gcS4WM4Z7twvMLp/D8T6BUw7s4QQJqwu\n2igaY+HxKutK9rEsr5zsrXs4WOf0C3SJEdL2bWLGtAuZMKRXxPYLtEU0Hh8dxRKCMSG0pbSKpbll\nfJBXTqlPv8DxqU6/wLnDk/lydSVnDesdwlaaaGH3IRjTwUr31/J+XhnLcsvYUnawYf5An36BtCjr\nFzBHxu5DMKYTqqr18PG2CpbmlrGupLFfoGd8LJOHJZM5KpkTUq1fwISWJYQwYXXRRpESC49XWVuy\nj2W5ZWRvq6Dap19gUnovMkf1YfyQXsTFHr5fIFLi0V4sHsETUEIQkanA00AM8KKqPtFk+U+Bm4Ba\n4FvgRlUtdJd5gC8AAfJV9Xvt13xjwk9e6QGW5Zbzfl4ZZQfqGuafkNqDKaP6cM6w3vTqZr/FTPhp\ntQ9BRGKATUAmUAKsAqar6kafdSYDn6rqQRG5FThXVae7y/aqaq/DvYf1IZjObvf+Gt7PLWdZbhlb\nyxv7BQb1qu8X6MOgXvEhbKGJRKHoQxgPbFbVfAARWQhMAxoSgqou91l/JfDfPtNWFDUR6UCNh+xt\ne1iWW8a6kkrqf1r1jI9l8vBkpozsw/Gp3a1fwHQagSSENKDQZ7oIJ0m05EfAuz7T8SLyGVAHPKGq\nb7e5lVHA6qKNwjkWHq+ypngfS3PL+GTbHqrd0ee7xggT0pPIHJnM+CG96NpKv0BbhHM8QsHiETyB\nJITmft40W2cSkR8ApwKTfWanq+oOERkGvC8iX6rqVt/tsrKymDdvHunp6QAkJSUxZsyYhn/07Oxs\ngIiezsnJCav2hHI6JycnpO//5ptv8vo/FtElKZW+ifGcM240dd2S2JNyHB/klZP/1WoAeo0Yy4n9\nezBo32ZOGpjIheeNjch4hNt0NMcjOzubBQsWAJCenk5qaiqZmZm0l0D6ECYCv1TVqe70A4A207E8\nBfg9cI6qlrawr5eBf6jqW77zrQ/BhIv8ggJufPhJmDSd2LgEPDVV7Fr6Cn0mfY/45P4ApPWKJ3NU\nHzJHJDPQ+gVMCIWiD2EVMFJEMoDtwHTgGt8VROQU4AXgIt9kICK9gQOqWiMifYEzAL9EYkw4mTN3\nfkMyAIiNSyB1ynWUfpTFTXf+lMyRfTiun/ULmMjUaqFTVT3AHcBiYD2wUFU3iMgsEbnUXe03QA/g\nbyKyVkT+7s4/HlgtImuBZcBjvlcnmUb1p4UmtLHYXVndkAzqxcYlcHxqAnecMYTjQ3DzmB0b/iwe\nwRPQxdCqugg4tsm8R3xeX9DCdiuAk46mgcZ0pMpaL56aKr+k4Kmpol+ilYZM5LNHJ4YJu2qiUahi\n8d6mUiqGTaZk8Z/x1DhjD3hqqmDFQu6dOSMkbQI7NpqyeASP3S5pDLCyoIKnPiogPrk/99xxKzkf\n/IPSympSEuO5d/Z9ZLhXwBkTySwhhAm7trpRR8fi6537mb1sK16Fa07uzw2nD4Ipp3TY+7fGjg1/\nFo/gsZKRiWr55VX8fHEe1R7lomP6MOO0gaFukjEhYwkhTNgvnkYdFYtv99fw4KI89lV7mJjei5+c\nlR6Wl5PaseHP4hE8lhBMVNp7sI6H3s1j9/5aTkjtwUPnDyM2JvySgTEdyRJCmLBrqxsFOxbVdV4e\nWbKF/D0HyejdjV9dOJxuYTxWsR0b/iwewRO+/xcYEwQer/Lo+9tYv3M/fXt0ZfbUETY2gTEuSwhh\nwuqijYIVC1Xl99mFrCiooGd8LI9NHUFqYlxQ3qs92bHhz+IRPJYQTNSY//l2Fm0qJT5W+PWFI8hI\nTmh9I2OiiCWEMGF10UbBiMXb67/lr+t2EiPwcOYwTujfo93fI1js2PBn8QgeSwgm4i3fUs5zK4oA\n+OnZ6UxMTwpxi4wJT5YQwoTVRRu1ZyzWluzjN//JR4EbThvIRcektNu+O4odG/4sHsFjCcFErNzd\nB5i1ZAu1XuV7o/sx/eT+oW6SMWHNEkKYsLpoo/aIxfa91Tz8Xh4Har1MHt6bWyemheVdyIGwY8Of\nxSN4LCGYiFN+oJYHF+VSXlXHKYMS+dnkDGI6aTIwpiNZQggTVhdtdDSxOFDj4eH38ijZW8PIlAR+\nMWU4cbGd+zC3Y8OfxSN4Ovf/Kcb4qPV4mbV0K7mlVQzqFcfsi0bQIy421M0yptOwhBAmrC7a6Ehi\n4VXlyQ8LWFuyj97duvDo1JEkd+8ahNZ1PDs2/Fk8gieghCAiU0Vko4hsEpH7m1n+UxFZLyLrRGSJ\niAzxWXa9u903InJdezbeGHAeSTF3ZTEf5JWT0DWG2VNHMKiXjYFsTFuJqh5+BZEYYBOQCZQAq4Dp\nqrrRZ53JwKeqelBEbgXOVdXpIpIMrAbGAQJ8DoxT1Qrf91i2bJmOGzeuHT+WiSavf7GTF1eV0CVG\nmH3RCE5J6xnqJhnTIdasWUNmZma7XTERyBnCeGCzquarai2wEJjmu4KqLlfVg+7kSiDNfX0RsFhV\nK1R1D7AYmNo+TTcGFm8q5cVVJQjwP5MzLBkYcxQCSQhpQKHPdBGNX/jN+RHwbgvbFreybdSyumij\nQGPxaUEFv/uoAIAfTxrMuSOSg9mskLFjw5/FI3gCeRB8c6cjzdaZROQHwKnA5LZsm5WVxbx580hP\nTwcgKSmJMWPGNFxeVn8ARPJ0Tk5OWLUnlNM5OTmtrp9ffpCFu/viVTidfPqW7wf6hUX7QxGPaJqO\n5nhkZ2ezYMECANLT00lNTSUzM5P2EkgfwkTgl6o61Z1+AFBVfaLJelOA3wPnqGqpO286Tn/Cre70\nC8AHqvq677bWh2DaoqD8ID/95yb2VXu46Jg+3HN2eI6FbEywhaIPYRUwUkQyRCQOmA6847uCiJwC\nvAB8tz4ZuN4DLhCRJLeD+QJ3njFH5Nv9NTy4KJd91R4mDOnFT86yZGBMe2k1IaiqB7gDp0N4PbBQ\nVTeIyCwRudRd7TdAD+BvIrJWRP7ublsO/BrnSqNPgVlu57JpwuqijVqKxb7qOh5alMe3+2s5IbUH\nD2cOIzYm8pOBHRv+LB7BE9Bgsqq6CDi2ybxHfF5fcJht5wPzj6x5xjiq67w8sngL+eUHSe/djV9d\nOJxuXey+SmPak/0fFSbs+SyNmsbC41UefX8bX+3cT98eXXl06gh6dQvot0xEsGPDn8UjeCwhmLCm\nqjzzcSErCiroGR/LY1NHkJoYF+pmGRORLCGECauLNvKNxStrdvDuN6XExwq/unA4GckJIWxZaNix\n4c/iETyWEEzYeufrb3lt7Q5iBB46fxij+yeGuknGRLRW70PoCHYfggHILyhgztz57K6sxqvC9iFn\nEZ/cn3vOTmfqsZ1vLGRjgq2970OInp45E9byCwq48eEnYdJ0YuMS8NRUsWvxn7n91lssGRjTQaxk\nFCaivS46Z+78hmSwN28dsXEJDLrwejZ9+K9QNy3kov3YaMriETyWEExY2F1ZTWycf4dxbFwCpfur\nQ9QiY6KPJYQwEe3XVqf0iMdTUwVArxFjAfDUVJGSaAPdRPux0ZTFI3gsIZiwMHDCRZQs/nNDUvDU\nVMGKhdw7c0ZoG2ZMFLGEECaiuS66LLeM/3wbx4AzLyMtdxExS37P6OKlvDT7PjLcR6JHs2g+Nppj\n8Qgeu8rIhFTu7gM87Q5yc88lp/KdE6aSnZ1tZQFjQsDOEMJENH4BVhysY9bSrVR7lIuO6cOlx/cF\nojMWh2Px8GfxCB5LCCYknAfWbWVnZQ3H9uvOnWcMsXENjAkxSwhhItrqoi+uKmFtSSW9u3XhF1OG\nEefzKOtoi0VrLB7+LB7BYwnBdLgP8srJytlFrMD/Zg6jXw97eqkx4cASQpiIlrpoXukBfvdhPgC3\nThzMSQMPfWBdtMQiUBYPfxaP4LGEYDrMXp9O5AtG9eG7J/QNdZOMMT4sIYSJSK+LerzKYx9sY8e+\nGkb1TeCuM1vuRI70WLSVxcOfxSN4AkoIIjJVRDaKyCYRub+Z5WeLyOciUisilzdZ5hGRNSKyVkT+\n3l4NN53JwTIkAAAcm0lEQVTL/M+383nxPpK6deGRKcOJt/GQjQk7rd6YJiIxwLNAJlACrBKRt1V1\no89q+cD1wH3N7GK/qtpgB62I5Lroh1vKef2LncQIPHz+0FaHwIzkWBwJi4c/i0fwBHKn8nhgs6rm\nA4jIQmAa0JAQVLXAXdbcaDt2cXkU21pWxZMfOnci3zIhjbGDeoa4RcaYlgRy3p4GFPpMF7nzAhUv\nIp+JyCciMq1NrYsikVgX3Vddx6ylWzhY5yVzZDKXje4X0HaRGIujYfHwZ/EInkDOEJr7hd+WcTfT\nVXWHiAwD3heRL1V1q+8KWVlZzJs3j3T3QWZJSUmMGTOm4dSw/gCI5OmcnJywas/RTnu9ynsHBlGy\nt4ae327g9FGDERka0PY5OTkhb384TVs8LB7109nZ2SxYsACA9PR0UlNTyczMpL20OqayiEwEfqmq\nU93pBwBV1SeaWfdl4B+q+lYL+2p2uY2pHHleXl3CX9ftpFd8LM9+71gG9LRxDYxpb+09pnIgJaNV\nwEgRyRCROGA68M5h1m9onIj0drdBRPoCZwBfH0V7TSeQvXUPf11X34k8zJKBMZ1EqwlBVT3AHcBi\nYD2wUFU3iMgsEbkUQEROE5FC4ErgBRHJcTc/HlgtImuBZcBjTa5OMq7OXhfNLyjgrod/xeW3P8jt\nD86iunwnN50+iFPS2t6J3Nlj0d4sHv4sHsET0HgIqroIOLbJvEd8Xq8GhjSz3QrgpKNsowlz+QUF\n3PjwkzBpOrHDEuhbU0X5+3/htKuPD3XTjDFtYHcHhYnOfG31nLnznWQQlwBAbFwCyef/kN/98c9H\ntL/OHItgsHj4s3gEjyUEc9Tyyw40JIN6sXEJlFZWh6hFxpgjYQkhTHTWuujiTaVsKavGU1PlN99T\nU0VK4pF1JnfWWASLxcOfxSN4LCGYI+LxKnNXFvHkhwX0Oe0iKv7zakNS8NRUwYqF3DtzRmgbaYxp\nk1bvQ+gIdh9C57K/xsPs97eyumgfsQJ3njmE0d0PMGfufEorq0lJjOfemTPIcG80NMYER3vfhxDQ\nVUbG1CuuOMjPF2+hqKKapG5d+HnmsIZBbp6Z/YsQt84YczSsZBQmOkNd9POivdz59iaKKqoZltyN\nZ6Yd0+yIZ0erM8SiI1k8/Fk8gsfOEEyrVJW3v97NCyuL8CpMykji/skZdI+LDXXTjDHtyPoQzGHV\nerw8+0kR735TCsA1J/fn+tMGEtPCaGfGmI5jfQimw+ypquXXy7aRs6OSuFjh3nPSOW9En1A3yxgT\nJNaHECbCrS66tayKO9/eRM6OSlK6d2XOpaM6LBmEWyxCzeLhz+IRPHaGYA7xSf4envhPPlW1Xo7t\n151fThlOSo+uoW6WMSbIrA/BAM4D6ua8MJ+vd+5jR2UdKaddxNTTjuOes9OJ72InksaEI+tDMO0u\nv6CAGx7+LTLpGmKHJ9C/porK5a9yzdXHWzIwJorY/+1hIpR10cefe8lJBj5PK02c/IMjflrp0bIa\nsT+Lhz+LR/BYQohyO/fV8Gl+hT2t1BhjCSFchOIZ73mlB7j7H99Q46Vdn1Z6tOx59/4sHv4sHsFj\nCSFKrS3ex73/3EzZgTrOuvh7eFf81Z5WakyUs4QQJjqyLrost4yH38vjQK2XycN788wPz2b+7J8x\nungp/XKyGF28lJdm3xeyp5VajdifxcOfxSN4ArrKSESmAk/jJJAXVfWJJsvPdpefBHxfVd/yWXY9\n8DCgwGxVfaWd2m7aSFX5W84u5n1WAsAVJ/bj5glpxIiQkZ5uTys1Jsq1eh+CiMQAm4BMoARYBUxX\n1Y0+66QDvYD7gHfqE4KIJAOrgXGAAJ8D41S1wvc97D6E4PN4lbmfFvP39d8CcMuENK4ckxriVhlj\njkYo7kMYD2xW1XwAEVkITAMaEoKqFrjLmmaXi4DF9QlARBYDU4HXj77pJlA1dV6eWJ7PR1v30DVG\nuG9yBueNSA51s4wxYSaQPoQ0oNBnusidF4im2xa3YduoEqy66L7qOh5clMdHW/fQvWsMs6eOCPtk\nYDVifxYPfxaP4AnkDKG505FAn3cR0LZZWVnMmzePdLcTMykpiTFjxjRcXlZ/AETydE5OTrvv/5ix\n43n4vTxyVq8kKb4LT911JcP6JITF5z3cdE5OTli1J9TTFg+LR/10dnY2CxYsACA9PZ3U1FQyMzNp\nL4H0IUwEfqmqU93pBwBt2rHsLnsZ+IdPH8J04FxVvdWdfgH4QFX9SkbWh9D+tpZV8fCiPHYfqCWj\ndzdmTx1BamJcqJtljGlH7d2HEEjJaBUwUkQyRCQOmA68c5j1fRv3HnCBiCS5HcwXuPNMEH1Rso97\n/rmZ3QdqOXFAD373nVGWDIwxrWo1IaiqB7gDWAysBxaq6gYRmSUilwKIyGkiUghcCbwgIjnutuXA\nr3GuNPoUmKWqe4LzUTq39qqLLt9SzkOL8thf4+HsYb15fOpIesZ3rmcYWo3Yn8XDn8UjeAL6plDV\nRcCxTeY94vN6NTCkhW3nA/OPuIWmVfkFBcyZO58NOyvZvq+WlNMu4uozT+TWiWnExthQl8aYwHSu\nn44R7Eifz+L36Oph9Y+ufo1Lrj6+0yYDe1aNP4uHP4tH8NijKzq5OXPnN/Po6v8O2aOrjTGdlyWE\nMHGkddGSiqqIe3S11Yj9WTz8WTyCxxJCJ6aq7KqsC6tHVxtjOi9LCGHiSOqiy7fsQUZnsn3xnyPq\n0dVWI/Zn8fBn8Qge61TupPZV1/HciiLik/vzwN0/ZvXSdyitrCYlMZ57Q/joamNM52UJIUxkZ2e3\n6ZfPvM9K2HOwjhMH9OAH547iuvPGBrF1HautsYh0Fg9/Fo/gsZJRJ/Tl9kre/aaUrjHCT85MJ0Y6\n5+WlxpjwYgkhTAT6i6fG4+Xp7AIAvn9yf9KTuwWzWSFhv/78WTz8WTyCxxJCJ/P6FzspqqhmcFI8\n08f2D3VzjDERxBJCmAjk2uqC8oMsXLcTgJ+clU5cbGT+89l15v4sHv4sHsETmd8oEcirytMfF1Dr\nVS4+NoWTBiaGuknGmAhjCSFMtFYXfe+bUr7asZ/khC7cNH5QB7UqNKxG7M/i4c/iETyWEDqBsgO1\n/OmzEgBunTi40z3O2hjTOVhCCBOHq4s+v7KIyhoPpw/uxbnDe3dgq0LDasT+LB7+LB7BYwkhzH1W\nWMHyLXuI7xLDnWcORuyeA2NMkFhCCBPN1UWraj384eMiAK4fN4ABPaPjgXVWI/Zn8fBn8QgeSwhh\n7C9rdrCzsoaRKQlcdmJqqJtjjIlwASUEEZkqIhtFZJOI3N/M8jgRWSgim0VkhYiku/MzROSAiKxx\n/55r7w8QKZrWRTfvPsBbX+0iRuAnZ6d32tHPjoTViP1ZPPxZPIKn1ctVRCQGeBbIBEqAVSLytqpu\n9FntR0CZqo4Ske8DvwGmu8tyVXVcO7c7onm8ylMfFeBVuPzEfhzTt3uom2SMiQKBnCGMBzarar6q\n1gILgWlN1pkG1I/ZmIWTPOpFz0/bo+BbF/37+m/JLa0iNbEr1586MIStCg2rEfuzePizeARPIAkh\nDSj0mS5y5zW7jqp6gD0i0sddNlREPheRD0TE/iVbsWNfNfM/3w7AnWcMIaFrbIhbZIyJFoEkhOZ+\n4Wsr64i7znYgXVVPBe4FFoiIPXOhGdnZ2agqz35SRHWdl3OG9WZCelKomxUSViP2Z/HwZ/EInkBu\neS0CfIffGozTl+CrEBgClIhILNBLVcvdZTUAqrpGRPKAY4A1vhtnZWUxb9480t1RvpKSkhgzZkzD\nqWH9ARDJ0zk5OXgGjeazwr3U5n/JuKEZwLCwaV9HTufk5IRVe0I9bfGweNRPZ2dns2DBAgDS09NJ\nTU0lM9O3Qn90RLXpj/0mKzhf8N/g9AtsBz4DrlHVDT7r3AacqKq3ich04HuqOl1E+uJ0NntFZDiw\nHBijqnt832PZsmU6blx09zvvq67jpqwNlFfVcfdZQ7jkuL6hbpIxJsytWbOGzMzMduunbfUMQVU9\nInIHsBinxPSiqm4QkVnAKlX9J/Ai8BcR2QyU0niF0TnAr0SkFvAAM5smg2iXX1DAnLnzWVe8l7Iq\nDxMvmsbFx0bOcJjGmM4joPsQVHWRqh6rqqNU9XF33iNuMkBVq1X1anf5RFXd5s5/S1VPVNVTVPU0\nVf130D5JJ5RfUMCNDz/J+rQp1AwcTf/JV7FpWRaFhYWtbxzBrEbsz+Lhz+IRPHancgjNmTsfJk0n\nNi4BgNi4BLqeea0z3xhjOpglhBDaXVndkAx6jXDKRLFxCZRWVoeyWSFn15n7s3j4s3gEjyWEECqv\n8uCpqfKb56mpIiUxOh5iZ4wJL5YQQuTNnF0cGHUuJYv/jKemir1565zksGIh986cEermhZTViP1Z\nPPxZPILHEkII/CevnLmfFhOf3J9f3HM7o4uX0nvLckYXL+Wl2feRkZ7e+k6MMaadtXofQkeIpvsQ\n1pXs4+FFedR6lZvGD+Lqk/qHuknGmE6qve9DsDOEDrS1rIpfLtlCrVf53uh+XDXGxjgwxoQPSwgd\nZFdlDQ8tyuNArZezh/Vm5oQ0v+EwrS7ayGLhz+Lhz+IRPJYQOsC+6joeXpRH6YFaThzQg/snZ0TV\ngDfGmM7BEkKQ1dR5eWTJFvL3HCSjdzdmXTCcuC6Hht2urW5ksfBn8fBn8QgeSwhB5PEqTyzP56sd\n++nbvSuzp46gZ3wgD5g1xpiOZwkhSFSVF1YW89HWPfSIi2X21BGkJsa1uL7VRRtZLPxZPPxZPILH\nEkKQ/O3LXbz99bd0jRF+OWUYw/okhLpJxhhzWJYQgmDp5jLmrXLGEPqfczM4eVDPVrexumgji4U/\ni4c/i0fwWEJoZ58X7WXOh/kA3DoxjcnDk0PcImOMCYwlhHaUu/sAv1q2FY/ClWNSufzEwG88s7po\nI4uFP4uHP4tH8FhCaCc79lXzv+/lUVXr5dzhvblp/KBQN8kYY9rEroE8CvXDX+7ce5At5dV0P/kC\nxp8wgvsmZxAjbbvxzOqijSwW/iwe/iwewWMJ4QjVD3/JpOnEDkugV00V3y59hRsvu5+4WDvxMsZ0\nPgF9c4nIVBHZKCKbROT+ZpbHichCEdksIitEJN1n2YPu/A0icmF7Nj5UVJXZz750yPCX/aZcx/Mv\n/+WI9ml10UYWC38WD38Wj+BpNSGISAzwLHARMBq4RkSOa7Laj4AyVR0FPA38xt32BOBq4HjgYuA5\nkTbWUkJMVdmxr5qPtu7hpVUlPPhuLle/9hUfbytvSAb1jmb4y5ycnPZobkSwWPizePizeDRat25d\nu+4vkJLReGCzquYDiMhCYBqw0WedacAj7uss4A/u6+8CC1W1DtgmIpvd/X3a9E3uevhX3DtzRkCD\nw9TX7ndXVtM3MT7g7VrbR/qQIWzfV8Pm3QfI3X2ATburyC09wL5qzyHbd42NxVNT5ZcUjmb4y4qK\niiPaLhJZLPxZPPxZPBp98cUX7bq/QBJCGlDoM12E86Xe7Dqq6hGRChHp485f4bNesTvvEOvTpvCD\nB3/Dr352J2mDh7TYmOKiQn7x2z8Qd+a1xMYlsLOmKqDtWtvHd+95lEFnXoan56GXiiZ168KovgmM\nSunOqL7dGdk3gYMX3M2P/ncOHrds1DD85ez7AmqDMcaEm0ASQnMlnqbDrLW0TiDbAk65Je7Ma7nz\nsRcYdMF1LTamZMkr9J98rV/tPpDtWttH3/Ovo3D53xj93RsZ1df94k9JYFTf7vTr0ZVDKl09M3hp\n9n3MmTuf0spqUhLjufcohr8sKCg4ou0ikcXCn8XDn8UjeAJJCEWA77fcYKCkyTqFwBCgRERigSRV\nLReRInf+4bZl3bp1HJ/jnPqcOABOjlnL2LFjm2/NuB+6L3zzSjeY+ENayDVt3EctUOH8lUFRmROA\nlsy44tKG16W7d1O6e3dgbWjitNNOY82aNUe0baSxWPizePiL5nisW7fOr0zUo0ePdt1/q2Mqu1/w\n3wCZwHbgM+AaVd3gs85twImqepuITAe+p6rT3U7l14AJOKWiJcAoDYeBnI0xxvhp9QzB7RO4A1iM\nc1XSi6q6QURmAatU9Z/Ai8Bf3E7jUmC6u+3XIvIG8DXOT+/bLBkYY0x4avUMwRhjTHSwW2o7gIgM\nFpH3ReRrEckRkbvc+ckislhEvhGR90QkyWebZ9wb+taJSAsdKp2XiMSIyBoRecedHioiK91Y/FVE\nurjzW7zpMVKISJKI/M29eXO9iEyI8mPjpyLylYh8KSKvucdA1BwfIvKiiOwUkS995rX5eBCR692b\nib8RkYCuuLGE0DHqgHtU9QRgEnC7e3PfA8BSVT0WeB94EEBELgZGuDf6zQReCE2zg+punFJivSeA\nOW4s9uDc7Agt3PQYYX4P/FtVjwdOxrnHJyqPDREZBNwJjFPVk3DK2tcQXcfHyzg3Avtq0/EgIsnA\nL4DTcfpwH/FNIi1SVfvr4D/g78AUnP/x+7vzBgAb3NcvAN/3WX9D/XqR8IdztdkS4FzgHXfet0CM\n+3oi8K77ehEwwX0dC3wb6va3cyx6AnnNzI/WY2MQkA8k4ySDd4ALgF3RdHwAGcCXR3o84PTjPu8z\n/3nf9Vr6szOEDiYiQ4GxwEqcf+CdAKq6A6i/K67pzYAt3tDXST0F/Az3ul8RSQHKVdXrLi+i8fP6\n3fQI7HFveowUw4HdIvKyW0L7o4h0J0qPDVUtAeYABTifrQJYA+yJ0uOjXmqAx0N9bI7oOLGE0IFE\nJBHn0R53q2olLd84EfANfZ2NiFwC7FTVdTR+TuHQz6w+y/x2QYTEwtUFGAf8P1UdB+zHKQ9E3bEB\nICK9cR6Fk4FzttAD5zloTUXL8dGalj7/ER0nlhA6iNsJlgX8RVXfdmfvFJH+7vIBOKfF4GT5Vm/o\n66TOBL4rIluAvwLn49R+k9wHKYL/522IhXtPTC9VLe/YJgdVEVCoqqvd6TdxEkQ0HhvglFK3qGqZ\n+4v//4AzgN5RenzUa+vxEMgNxYewhNBxXgK+VtXf+8x7B5jhvp4BvO0z/zoAEZmIc7q8s2OaGVyq\n+pCqpqvqcJw65/uq+gPgA+Aqd7Xr8Y/F9e7rq3A61CKG++9aKCLHuLMygfVE4bHhKgAmikg398nI\n9fGItuOj6VlzW4+H94AL3CvYknH6Yd5r9V1D3XkSDX84v4o9wDpgLU5NdCrQB1iKcyf4EqC3zzbP\nArnAFzhXXIT8cwQhLpNp7FQehvMU3E3A60BXd3488AawGaffZWio2x2EOJwMrHKPj7eApGg+NnCe\nnLwB+BL4M9A1mo4PYAHOr/lqnAR5A04ne5uOB5zEsdmN2XWBvLfdmGaMMQawkpExxhiXJQRjjDGA\nJQRjjDEuSwjGGGMASwjGGGNclhCMMcYAlhDMURCRs0RkQ+trdiwRWSAi3w1xGx4Rkb+Esg2BEpEP\nROTGdtjPmyJyYXu0yYSGJQRzxFQ1W51HNocNERkDnKSq74S6LUTXM3UAHgceDXUjzJGzhGCOiPvc\nmHA0E2cc76jX0f9GqroK6Cki4zryfU37sYQQRUTkfhH5W5N5vxeRp93XM9xR3faKSK6I3OKz3mQR\nKRSR/xGR7cBL9fOa7D/X3f4rEfmez7LrReQjEfmtiJSJSJ6ITPVZniwiL4lIsYiUishbPssuFZG1\nIlIuItnuWUBLLgaWN3nfbHdUqT3u5zvfZ3kvEZknIiXu5/u1+wwdxPG/IrJNRHaIyHwR6eUuyxAR\nr4jc7La5WETuOUzsJ4rIx+5nWCsikw+z7lYReUCc0dNKxRlBKy6QeLjb/o+IfAFU+jwQznf/F4gz\nOlu5iPwBn2fmiMhwEVkmIrtFZJeIvOrzme8Tkawm+/qDiPzOZ9Zy4JKWPpsJc6F+bof9ddwfztMP\nK4FEdzoG55kpp7vTF+M+CwY4G+dRzGPd6clALU5JoCvOM2QmAwU++7+CxkE8rnLfq376epxns9yI\n8wV0K1Dss+2/cJ5+2gtnoJOz3fnjgJ3Aae52PwS24j7Lpsnn6w54gRSfede77b7L3e/VOCNu9XaX\n/x14DugG9MV5Hs7N7rIbcZ4Dk+Hu+03gFXdZhvter7nbnojzBMrz3eWP+KybBuwGLnKnM93plBb+\nnbbiPMdnENAbyAZ+FUg83Ndr3G3jm9l3Cs4YA5e58fiJG58b3eUj3PZ1cdf9D/A7d9kAYB/OE0Vx\nt99Zf4y4834KZIX6WLe/I/yOCHUD7K+D/8HhQ+AH7usLgM2HWff/gDvd15OBg75fxDRJCM1svxb4\njvv6emCTz7IE9ws11f2iqav/ommyj+eAWU3mbcRNGE3mD8J5iGCcz7zrgaIm630K/Lf73gd9vzhx\nnsC6zH29FLjVZ9kxQA1OIq1PCKN8lj8B/Ml97ZsQ/gf4c5M2LAJ+2ELctuImJXf64vp/p9bi4W57\n/WH+TX4IfNJkXmF9Qmhm/WnA5z7T/wJ+5L6+FPiqyfo34Qz1GPJj3f7a/mclo+jzV5wxanH/u6B+\ngYhcLM5A5aUiUo7zRdTXZ9tvVbW2pR2LyHU+pYxyYHST7XfUv1DVKvdlIs7z3MtUdW8zu80A7nXL\nTGXufgfjfPk3tcf9b88m84ubTOe722fgnO1s99n3C0A/d7364Rx9t+uCM0QhOJ3GRc3st7nPcHWT\nz3AmMLCZdeu1tN9A4uG7bVOD8B9JC99pEeknziD2RSKyB3gV/3/DV4AfuK//G2h6JVVPGv8dTCdj\nCSH6/A04V0TScMoGCwDcGnUWziDl/VQ1GXgX/2eyt3jVjIikA38EblPVZHf79U22b0kh0Ke+Vt3M\nstmq2sf9S1bVRFV9vemKqnoAyMP5Je+r6dCB6TilskKcM4QUn333Vmdwd9x1Mny2y8Apr9SPPyD4\nD05Sv9/mPsMrTT5DT1U93IDwvvvN8NlvIPE43NVN2/EfOKXpez2Gc+Zzoqr2xvny9/03/DtwkoiM\nxjlDaNqBfzzOY5hNJ2QJIcqo6m6cjr+XcUam+sZdFOf+7VZVr4hcDLTlmvIeOF8ku0UkRkRuwKmr\nB9KmHTjJ5zkR6S0iXUTkbHfxn4BbRWQ8gIj0EJH/EpEeLezu3zilLF+pInKnu9+rgOOAf7vvuxh4\nSkR6up3Iw0XkHHe7vwI/FZGh4gx/OhtYqI1j+wL8XEQS3C/IG4CFzbTpVeA7InKhG5tu4nTIN3c2\nUe92EUkTZ3zgB33229Z4NPUv4AQR+Z6IxIrI3Tglu3o9cfp+9ro/Gn7mu7GqVuP0pSwAPlXVpmcj\nk3H+LU0nZAkhOi3A6Ths+HWnzhjPdwF/E5EynFr6281vfihV3YAzOPpKnNLQaJzO0MNu5vP6hzj9\nCBtxfoHf7e73c+Bm4Fm3XZtoHCGrOX+isaRR71NgFE5H7q+BK7RxmMXrcBLh10AZzhlU/RfkSzgl\nkQ9xzjwO4MTI13KcwUmWAL9R1WWHfEjnS3Ma8BDwLU4J6D4O///fApxklev+zXb31Vo8Dnvvg6qW\n4nT4P4ETjxH4/zvNAk7FKfv8A+fLv6k/A2NwykcNROR0oFIbhwM1nYwNkGMijoi8Cryhqu+IyPU4\nnaDntLZdG98jA9iC08nubW39Nu57K06bw3I4SBEZgjOi2QD3h0T9/CycTvXWh2o0YalLqBtgTHtT\nZ4zmjhBI/0hEce9ruBendFbpu0xVrwxNq0x7sYRgzJEL1ul1WJ62i0h3nHLeVpwr0EyEsZKRMcYY\nwDqVjTHGuCwhGGOMASwhGGOMcVlCMMYYA1hCMMYY47KEYIwxBoD/DyTX6l5zBzNYAAAAAElFTkSu\nQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1110637f0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "df2.plot.line(x='variance (people per day)', y='p-value', style='-o')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The point that this plot makes is: depending on the attendance at the ER, the numbers used in the study can or cannot be statistically significant. If the fluctuation in attendance for a 2100 people ER is above 200 people per day (roughly 10%), the numbers from the study would be judged non significant. On the other hand, if the numbers are around 50 people per day (low fluctuation) then there is a real effect. It would be interesting to contact the author about this point."
   ]
  }
 ],
 "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.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}