textbook/Choosing_a_Sample_Size.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# HIDDEN\n",
    "from datascience import *\n",
    "import numpy as np\n",
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plots\n",
    "plots.style.use('fivethirtyeight')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Choosing a Sample Size ###\n",
    "Candidate A is contesting an election. A polling organization wants to estimate the proportion of voters who will vote for her. Let's suppose that they plan to take a simple random sample of voters, though in reality their method of sampling would be more complex. How can they decide how large their sample should be, to get a desired level of accuracy?\n",
    "\n",
    "We are now in a position to answer this question, after making a few assumptions:\n",
    "- The population of voters is very large and that therefore we can just as well assume that the random sample will be drawn with replacement.\n",
    "- The polling organization will make its estimate by constructing an approximate 95% confidence interval for the percent of voters who will vote for Candidate A.\n",
    "- The desired level of accuracy is that the width of the interval should be no more than 1%. That's pretty accurate! For example, the confidence interval (33.2%, 34%) would be fine but (33.2%, 35%) would not.\n",
    "\n",
    "We will work with the sample proportion of voters for Candidate A. Recall that a proportion is a mean, when the values in the population are only 0 (the type of individual you are not counting) or 1 (the type of individual you are counting)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Width of Confidence Interval ###\n",
    "If we had a random sample, we could go about using the bootstrap to construct a confidence interval for the percent of voters for Candidate A. But we don't have a sample yet – we are trying to find out how big the sample has to be so that our confidence interval is as narrow as we want it to be.\n",
    "\n",
    "In situations like this, it helps to see what theory predicts.\n",
    "\n",
    "The Central Limit Theorem says that the probabilities for the sample proportion are roughly normally distributed, centered at the population proportion of 1's, with an SD equal to the SD of the population of 0's and 1's divided by the square root of the sample size.\n",
    "\n",
    "So the confidence interval will still be the \"middle 95%\" of a normal distribution, even though we can't pick off the ends as the 2.5th and 97.5th percentiles of bootstrapped proportions. \n",
    "\n",
    "Is there another way to find how wide the interval would be? Yes, because we know that for normally distributed variables, the interval \"center $\\pm$ 2 SDs\" contains 95% of the data.\n",
    "\n",
    "The confidence interval will stretch for 2 SDs of the sample proportion, on either side of the center. So the width of the interval will be 4 SDs of the sample proportion.\n",
    "\n",
    "We are willing to tolerate a width of 1% = 0.01. So, using the formula developed in the last section, \n",
    "\n",
    "$$\n",
    "4 \\times \\frac{\\mbox{SD of the 0-1 population}}{\\sqrt{\\mbox{sample size}}}\n",
    "~ \\le ~ 0.01\n",
    "$$\n",
    "\n",
    "So\n",
    "\n",
    "$$\n",
    "\\sqrt{\\mbox{sample size}} ~ \\ge ~ 4 \\times \\frac{\\mbox{SD of the 0-1 population}}{0.01}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### The SD of a collection of 0's and 1's ###\n",
    "If we knew the SD of the population, we'd be done. We could calculate the square root of the sample size, and then take the square to get the sample size. But we don't know the SD of the population. The population consists of 1 for each voter for Candidate A, and 0 for all other voters, and *we don't know what proportion of each kind there are.* That's what we're trying to estimate.\n",
    "\n",
    "So are we stuck? No, because we can *bound* the SD of the population. Here are histograms of two such distributions, one for an equal proportion of 1's and 0's, and one with 90% 1's and 10% 0's. Which one has the bigger SD? "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsIAAAEWCAYAAABsVfF9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlUFFf+NvCngCgCIq00YFhEBREVIRA3oqiYBBcQiRqj\nI2owuEBMdFxwiXFjBFyQTIK4IOokbnEP6rhMgkYDQnSi4rgRDeIKKDYqCirU+4c/+rVtlga7aaCf\nzzmeI7duV38vZcjD7Vu3BJlMJoKIiIiISMfoabsAIiIiIiJtYBAmIiIiIp3EIExEREREOolBmIiI\niIh0EoMwEREREekkBmEiIiIi0kkMwkRERESkk7QehJOTkzF8+HC0a9cOEokEW7ZsUeoTEREBZ2dn\nNG/eHL6+vrh06ZLC8WfPnmH69Olo3bo1rK2tMXz4cNy+fbumhkBEREREdZDWg3BBQQHat2+PyMhI\nGBkZKR2PiYlBXFwcli5diqSkJEilUgQEBKCgoEDeZ+bMmdi/fz8SEhLw73//G48ePcKwYcMginxW\nCBERERGVTahNT5azsbHB0qVLMXz4cHlb27ZtMX78eEyZMgUAUFhYCEdHR4SHh2P06NF4+PAhHBwc\nEBcXh8GDBwMAbt26BRcXF+zcuRO9e/fWyliIiIiIqHbT+oxwRTIzM5Gdna0QZg0NDeHp6YnU1FQA\nwB9//IEXL14o9LG2toaTk5O8DxERERHR62p1EM7JyYEgCJBKpQrtUqkUOTk5AIDc3Fzo6+ujadOm\n5fYhIiIiInpdrQ7CRERERESaUquDsIWFBURRRG5urkJ7bm4uLCws5H2Ki4uRl5dXbh8iIiIiotfV\n6iBsb28PS0tLJCUlydsKCwuRkpKCrl27AgDc3NxgYGCg0OfWrVu4fPmyvA8RERER0eu0HoQLCgqQ\nnp6Oc+fOoaSkBDdv3kR6ejpu3rwJAJg4cSJiYmKQmJiICxcuICQkBCYmJvIdIkxNTREYGIh58+bh\n2LFjOHv2LCZMmAAXFxf07NlTm0OrdTIyMrRdgtZw7LpHV8cN6PbYiYiqwkDbBfzxxx/w8/ODIAgA\nXj48IyIiAsOHD0dsbCy+/PJLFBYWYsaMGZDJZPDw8MCuXbtgbGwsP0dkZCQMDAwQFBSEwsJC9OzZ\nE6tXr5afk4iIiIjodbVqH2HSrIyMDDg6Omq7DK3g2HVv7Lo6bkC3x05EVBVaXxpBRERERKQNDMJE\nREREpJMYhImIiIhIJzEIExEREZFOYhAmIiIiIp3EIExEREREOolBmIiIiIh0EoMwEREREekkBmEi\nIiIi0kkMwkRERESkkxiEiYiIiAAMGDAAvr6+2i6jQqIoYuHChejQoQOaNm2KkSNHarukOo1BmIiI\niBRs3rwZEolE/sfc3Bzt27dHaGgo7ty5o+3y3khaWhoiIyPx8OFDpWOCIEBPr3ZHox9//BErVqzA\ngAEDsGrVKoSEhJTb988//8Ts2bPRr18/WFtbQyKR4PTp02X2lUgk2LJli6bKBgDk5+dj8uTJcHBw\ngLW1NXx9fXHmzBmVXhsZGanwb7L0T/Pmzd+oJoM3ejUREREpyc59gHsPlINWTTKXmMJSKqn26wVB\nwKxZs2Bvb4/CwkKkpqZiy5YtSE5ORkpKCgwNDdVYbc1JTU3FkiVL8Le//Q2mpqYKx/bs2aOlqlT3\n22+/QSKRICoqqtK+aWlpWLNmDRwdHdGuXbtyQ3BNEEURH3/8MS5cuIBJkybB3Nwc8fHx8PPzQ1JS\nEhwcHCo9hyAIWLZsGRo3bixv09fXf6O6GISJiIjU7N6Dh1gYs1mrNXw9ecQbBWEA8Pb2hoeHBwAg\nMDAQZmZmWLlyJQ4cOICPPvqozNc8efIERkZGb/S+mqBKXQYGtT8W5ebmKgTBivTv3x8DBw6EiYkJ\nNm/erNUgvGfPHqSlpWHDhg3w9/cHAPj7+8PDwwOLFy9GQkKCSufx8/ODVCpVW121e/6fiIiIag0v\nLy+Ioojr168DADZt2gSJRIJff/0VM2bMQJs2bWBjYyPvn5WVhU8//RStWrVC8+bN4e3tjf379yuc\n88SJE5BIJNi+fTsWL14MZ2dnvP322xg8eDCuXbumVMOJEyfQv39/WFtbw87ODsOGDcPFixcV+kRE\nREAikeDixYsYN24cWrZsCU9PT0RGRuLrr78GAHTs2BESiQRNmzbFb7/9BuDlGmE/Pz+Fcz19+hRz\n585Fhw4dYGlpCQ8PD8TExEAURYV+EokEU6dOxf79++Hp6QlLS0t069YNP//8s0rf28reJysrCxKJ\nBIcOHZL//dXay2JmZgYTExOV3v91BQUF+Oqrr+Dq6gpLS0s4ODjA19cXKSkpCjVnZGQgLy+v0vP9\n9NNPMDc3l4dgAGjWrBkCAgJw8OBBFBUVqVRXSUkJHj16pPT9r67a/6sPERER1QqlwbRp06YAXn5U\nDQBhYWEwMzPDtGnT5Gtv7927hw8//BBPnjzB+PHj0axZM/z4448IDAxEfHy80oxyaeibNGkSZDIZ\nVq1ahYEDB+LEiRMwMzMDAPz6668YPHgw7O3tMXPmTBQVFSE+Ph4+Pj44evQoWrVqpVBXUFAQWrRo\ngblz5+LZs2fw8vLC1atXsXPnTkRGRsrH4eTkpPC6V40YMQLHjh1DYGAgXF1dcezYMSxYsAA3btzA\n8uXLFfqmpqbi4MGDCAoKgomJCVavXo3Ro0fj/Pnz8jGUp7L3MTc3x5o1a/Dtt9/izp07iIiIgCiK\n8trf1OtjnzJlCn766ScEBwfDyckJMpkMp0+fxvnz59GtWzcAwOnTp+Hn54eZM2ciLCyswvOfO3cO\nHTt2VGp3d3fHxo0bceXKFbi4uFR4DlEU8e677+Lx48cwMjJC//79ER4eDktLyyqO9v9jECYiIqIy\nPXz4EHl5eSgsLMTJkyexdOlSGBkZwcfHR6Ff48aNsX//foUbzaKjo5GTk4N9+/bB09MTADBmzBj0\n6tULc+bMgb+/v8L6znv37uHUqVPyj/27d+8Of39/xMbGYs6cOQCAr776Ck2aNMGRI0fkwfKjjz5C\n165dsXDhQmzYsEGhLicnJ6W2jh07YufOnejfvz9sbW0rHP+BAwdw9OhRzJ49G9OnTwfwMlyHhoZi\n/fr1CA4ORtu2beX9MzIykJqaCnt7e/kYunfvjh07duCzzz574/cZOnQoduzYgfz8fAwZMqTC2qvq\n9Vndw4cPY9SoUVi0aFGFrxMEocxfIF6XnZ2NLl26KLVbWVkBAO7evVthEDYzM8O4cePQuXNnNGjQ\nACkpKVi7di1Onz6NY8eOqbxc5HUMwkT1WXE+rJvKIBRd0nYlNc666SOdHDdggAZv8Uc7vTlRFBVm\nbQVBQNu2bREVFSUPL6VGjRqltNvCkSNH4OrqKg/BAGBoaIixY8ciLCwMZ8+ehbu7u/zY8OHDFcKM\nl5cXnJ2dcejQIcyZMwfZ2dlIT0/HpEmTFGZXW7VqhX79+uHnn3+GKIryUCYIAj799NM3+h4cOXIE\n+vr6GD9+vEL7559/js2bN+Pw4cMKQdjLy0seggGgffv2aNy4MTIzM9X6PjXB1NQUp0+fxp07d8rd\nmaF79+4qLYsAXi6jaNiwoVJ7w4YNIYoinj59WuHrJ0yYoPC1n58f3N3dERwcjDVr1mDq1Kkq1fE6\n/rQkqs/EIjTKXYBGT4y1XUmNKyko0MlxFzdoAwODYG2XQfWAIAhYsmQJHB0dYWhoCBsbG1hbW5fZ\n79XwV+rGjRsYOHCgUnubNm0giiKysrIUgnDpsoZXtW7dGsePH5efD0CZuwu0adMGiYmJuH//PszN\nzeXtLVu2rHygFbhx4wYsLCyUdpdwdHSEnp4esrKyFNrL+v6YmZlBJpOp9X1qwqJFixASEoIOHTqg\nY8eO6NOnDz755BOVdncoS6NGjcpcB1xUVARBENCoUaMqn3PIkCH46quvcPTo0WoHYd4sR0RERGV6\n55130LNnT3Tp0qXMkFeqOiGmJtR0XeVt5aWuG7tqkr+/P86cOYPo6GjY2dlhzZo18PT0xM6dO6t1\nPktLS2RnZyu13717FwCUPmVQlbW1NR48eFCt1wIMwkRERKQBtra2yMjIUGq/cuUKBEGAnZ2dQvvV\nq1eV+l69elXer3Q9b3nnNDY2RrNmzSqtS5X1rKVsbW2Rk5ODR48eKbRnZGSgpKREaQzVVVPvU1VS\nqRSjR4/Gxo0bcf78edjb2yMiIqJa53JxccG5c+eU2k+dOgVDQ0O0adOmWufNyspS+BSgqhiEiYiI\nSO18fHxw9uxZnDx5Ut5WVFSEhIQEWFpaws3NTaH/1q1bkZ+fL//62LFjuHjxovzGPEtLS7i6umLr\n1q0KSw3++usvHDx4EB988IFKIdfY+OWSqcqWK5SOobi4GKtXr1Zoj42NhSAI+PDDDys9hypq6n1U\nVVJSovTkPVNTU7Ro0ULhGlVl+zR/f3/cv39f4aEl9+/fx969e+Hj46OwfvjmzZtKv/Dcv39f6Zzx\n8fG4d+8e3n//fZXH9jquESYiIiIlqn6cX16/yZMnY+fOnfj4448xbtw4mJubY9u2bbhy5Qri4+OV\nbq6TSqXo27cvRo4cCZlMhtWrV+Ptt99GaGiovM+iRYswePBgfPDBBxg1ahQKCwuxbt06NGrUCF99\n9ZVK9b7zzjsQRRELFizAkCFD0KBBA/Ts2bPM2eS+ffuid+/eiIiIQFZWlnxbs3379iEoKEhtN7Bp\n6n0ePnyI1atXQxAEnDt3DqIoYtOmTUhKSkKTJk0QHFz2/QSPHj1Cu3bt4Ofnhw4dOsDU1BQpKSn4\n+eefMW7cOHm/qmyf5u/vj7i4OEyaNAmXL1+WP1muuLgYs2fPVug7fvx4JCcnKyx5cHFxQUBAANq1\nawdDQ0OkpKRg165dcHFxwdixY6v1/QEYhImIiKgMqi4hKK+fubk5Dh06hPnz5yMhIQFPnz6Fs7Mz\nvv/+e/Tv31/pHJMnT0ZGRga+++475Ofnw9PTE1FRUQo7RPTo0QO7d+/G4sWLERkZCX19fXh6emLe\nvHll3mxXFjc3N8yfPx/x8fH4/PPPUVJSgsTERLz33ntljmfTpk2IiIjArl27sG3bNtjY2GDevHn4\n4osvlMZQ1vdC1e3FVH2fsmosj0wmw+LFixV20ijdTs7W1rbcIGxkZITg4GAkJSXh4MGDePHiBezs\n7BAeHq60e4Oq49PT08P27dsxb948rF27Fk+fPoW7uzvi4uLg6OiodM7Xf1H6+OOPkZaWhsTERBQV\nFcHW1hZffvklpk6d+kZrwQWZTFb3VnBTtWRkZCj9Y9MVOjv2FzkQ/wySfxSoSwoKCnRy3MUN2uDa\n82BY22hnTSG9lJ37APcePKy8owaZS0zf+BHLNeHEiRPw8/NDQkICAgICtF0O6RjOCBMREamZpVRS\nJ0Ioka7jzXJEREREpJMYhImIiEirqrKlGZE6cWkEERERaU1VHtNLpG6cESYiIiIincQgTEREREQ6\niUGYiIiIiHQSgzARERER6SQGYSIiIiLSSQzCRERERKSTan0QLikpQXh4OFxdXWFlZQVXV1eEh4ej\npKREoV9ERAScnZ3RvHlz+Pr64tKlS1qqmIiIiIjqglofhFesWIGEhAQsXboUv//+O6KiorBu3TpE\nR0fL+8TExCAuLg5Lly5FUlISpFIpAgICUFBQoMXKiYiIiKg2q/VBOC0tDX379sWHH34IW1tb9O3b\nF3379sWpU6fkfVatWoUpU6bA19cXbdu2RVxcHB4/fowdO3ZosXIiIiKqSwYMGABfX19tl1EhURSx\ncOFCdOjQAU2bNsXIkSO1XVKdVuuDcLdu3XD8+HFkZGQAAC5duoTjx4/Dx8cHAJCZmYns7Gz07t1b\n/hpDQ0N4enoiNTVVKzUTERHVZZs3b4ZEIpH/MTc3R/v27REaGoo7d+5ou7w3kpaWhsjISDx8+FDp\nmCAI0NOr3dHoxx9/xIoVKzBgwACsWrUKISEh5fb9888/MXv2bPTr1w/W1taQSCQ4ffp0mX0lEgm2\nbNmiqbIBAPn5+Zg8eTIcHBxgbW0NX19fnDlzRuXX7927F97e3mjevDlatmyJsWPH4tatW29UU61/\nxPLkyZPx+PFjdOnSBfr6+iguLsbUqVPx6aefAgBycnIgCAKkUqnC66RSKe7evauNkomISMcJz7Mh\nFN/Tag2ivjnEtyyr/XpBEDBr1izY29ujsLAQqamp2LJlC5KTk5GSkgJDQ0M1VltzUlNTsWTJEvzt\nb3+DqampwrE9e/ZoqSrV/fbbb5BIJIiKiqq0b1paGtasWQNHR0e0a9eu3BBcE0RRxMcff4wLFy5g\n0qRJMDc3R3x8PPz8/JCUlAQHB4cKX79hwwZMmTIF3bt3x6JFi5CXl4e4uDicOnUKx44dg5mZWbXq\nqvVBeOfOndi6dSsSEhLg5OSE9PR0hIWFoUWLFm/8cUDpLLMu0cUxl9LFsVs2BRoDOrteXifHXfIU\n0NfNf++Ojo7aLkFOKL4Hw9yFWq2hUPr1GwVhAPD29oaHhwcAIDAwEGZmZli5ciUOHDiAjz76qMzX\nPHnyBEZGRm/0vpqgSl0GBrU+FiE3NxeNGzdWqW///v0xcOBAmJiYYPPmzVoNwnv27EFaWho2bNgA\nf39/AIC/vz88PDywePFiJCQklPva58+fY8GCBejatSsSExPl7T4+PvD29saKFSuwYMGCatVV66/4\nvHnz8MUXX2DQoEEAAGdnZ2RlZWHFihUYOXIkLCwsIIoicnNzYW1tLX9dbm4uLCwsKjx3bfqhWRMy\nMjJ0bsyldHbsL3Ig5gLGxsbarqTGFRQU6OS4ixs0Ap7r3s83qhleXl6IjY3F9evXAQCbNm3C559/\njr1792Lfvn3Ys2cP7t27h7y8PABAVlYW5s2bh2PHjuHp06dwdnbG1KlTMWDAAPk5T5w4AT8/P6xZ\nswYZGRn4/vvvkZ+fj27dumHp0qVo1aqVQg0nTpzA4sWLcfbsWejr66Nbt26YP38+nJ2d5X0iIiKw\nZMkSJCcnY8WKFThy5AiaNGmCTz75BFFRURAEAR07dgTwcuY7MTER7733HgYMGAA9PT2FsPX06VMs\nXrwYu3fvRm5uLmxsbBAYGIgvv/wSgiDI+0kkEgQFBcHb2xv/+Mc/cPXqVbRq1Qrh4eHo06dPpd/b\nyt4nKysLrq6uEAQBoihCIpEo1F6W6s6SAi9/hkZERCAxMRF3795F48aN0bZtW8yZMwfdunWT13zz\n5k00a9YMTZs2rfB8P/30E8zNzeUhGACaNWuGgIAAbNu2DUVFRWjYsGGZr7148SJkMpk8C5ZydXWF\nk5MTdu7cWX+D8JMnT5TW6+jp6cm3T7O3t4elpSWSkpLg5uYGACgsLERKSgrCw8NrvF4iIqL66tq1\nawAgDz2lQTAsLAxmZmaYNm2afO3tvXv38OGHH+LJkycYP348mjVrhh9//BGBgYGIj49XmlGOiYmB\nKIqYNGkSZDIZVq1ahYEDB+LEiRPyQPfrr79i8ODBsLe3x8yZM1FUVIT4+Hj4+Pjg6NGj8tBcWldQ\nUBBatGiBuXPn4tmzZ/Dy8sLVq1exc+dOREZGysfh5OSk8LpXjRgxAseOHUNgYCBcXV1x7NgxLFiw\nADdu3MDy5csV+qampuLgwYMICgqCiYkJVq9ejdGjR+P8+fOVhtLK3sfc3Bxr1qzBt99+izt37iAi\nIgKiKMprf1Ovj33KlCn46aefEBwcDCcnJ8hkMpw+fRrnz5+XB+HTp0/Dz88PM2fORFhYWIXnP3fu\nnPyXj1e5u7tj48aNuHLlClxcXMp87bNnzwAAjRo1UjrWqFEjXLp0Cbm5uUrLZFVR64Nw3759ERMT\nAzs7O7Rt2xZnz57FypUrMWLECHmfiRMnIjo6Gg4ODmjdujWWLVsGExMTDB48WIuVExER1W0PHz5E\nXl4eCgsLcfLkSSxduhRGRkbyG9ZLNW7cGPv371eYuIqOjkZOTg727dsHT09PAMCYMWPQq1cvzJkz\nB/7+/tDX15f3v3fvHk6dOiX/2L979+7w9/dHbGws5syZAwD46quv0KRJExw5ckQeLD/66CN07doV\nCxcuxIYNGxTqcnJyUmrr2LEjdu7cif79+8PW1rbC8R84cABHjx7F7NmzMX36dAAvw3VoaCjWr1+P\n4OBgtG3bVt4/IyMDqampsLe3l4+he/fu2LFjBz777LM3fp+hQ4dix44dyM/Px5AhQyqsvapKZ/FL\nHT58GKNGjcKiRYsqfJ0gCGX+AvG67OxsdOnSRandysoKAHD37t1yg3Dr1q2hp6eHlJQUjBo1St5+\n//59XL58GQBw+/btagXh2n1rJIClS5di4MCBmDZtGrp27Yqvv/4aY8aMwVdffSXv8+WXXyIkJAQz\nZsxAnz59kJOTg127dunkx6JERETqIIoiPvroI7Ru3Rrt27fHZ599BktLS2zdulUeXkqNGjVK6dPb\nI0eOwNXVVR6CgZe7Oo0dOxbZ2dk4e/asQv/hw4crrH318vKCs7MzDh06BOBlkEpPT8eIESMUZldb\ntWqFfv364eeff4YoivJ2QRDkN9ZX15EjR6Cvr4/x48crtH/++ecQRRGHDx9WaPfy8pKHYABo3749\nGjdujMzMTLW+T00wNTXF6dOnK9wlpHv37sjLy8OMGTMqPd/Tp0/LXPrQsGFDiKKIp0+flvtaiUSC\nwYMHY/v27VixYgUyMzNx+vRpjB49Gs+fPwfwcjVAddT6GWFjY2MsXrwYixcvrrBfWFhYpdPyRERE\npBpBELBkyRI4OjrC0NAQNjY2CvfivNrv1fBX6saNGxg4cKBSe5s2bSCKIrKysuDu7i5vf30tMPBy\nJvD48ePy8wEoc3eBNm3aIDExEffv34e5ubm8vWXLlpUPtAI3btyAhYWF0u4Sjo6O0NPTQ1ZWlkJ7\nWd8fMzMzyGQytb5PTVi0aBFCQkLQoUMHdOzYEX369MEnn3xS6e4O5WnUqBGKioqU2ouKiiAIQpnL\nHl4VHR2NZ8+eYdGiRVi4cCEEQcD777+Pv/3tb9iwYQNMTEyqVVetnxEmIiIi7XjnnXfQs2dPdOnS\npcyQV6qyEKMtNV3Xq0s9XvXqTHVd4e/vjzNnziA6Ohp2dnZYs2YNPD09sXPnzmqdz9LSEtnZ2Urt\npVvdvv4pw+tMTEywYcMGXLhwAQcOHMDp06fx448/Ij8/H3p6etX+pYdBmIiIiNTO1ta2zG38rly5\nAkEQYGdnp9B+9epVpb5Xr16V9ytdz1veOY2NjdGsWbNK61JlPWspW1tb5OTk4NGjRwrtGRkZKCkp\nURpDddXU+1SVVCrF6NGjsXHjRpw/fx729vaIiIio1rlcXFxw7tw5pfZTp07B0NAQbdq0Uek8VlZW\n6NatG1q2bIni4mKcOHECnTt3rvaWfSoH4aioKFy4cKHc4xcvXlRpc2ciIiKq/3x8fHD27FmcPHlS\n3lZUVISEhARYWlrKd3oqtXXrVuTn58u/PnbsGC5evCi/Mc/S0hKurq7YunWrwlKDv/76CwcPHsQH\nH3ygUsgtvX+osuUKpWMoLi7G6tWrFdpjY2MhCAI+/PDDSs+hipp6H1WVlJQoPXnP1NQULVq0ULhG\nT58+RUZGhtKNdmXx9/fH/fv3FR5acv/+fezduxc+Pj4K64dv3ryp0l7oK1asQG5uLiZNmqTKsMqk\n8hrhyMhItGrVCu3atSvzeGkQ5jpdIiKiuk/Vj/PL6zd58mTs3LkTH3/8McaNGwdzc3Ns27YNV65c\nQXx8vNLNdVKpFH379sXIkSMhk8mwevVqvP322wgNDZX3WbRoEQYPHowPPvgAo0aNQmFhIdatW4dG\njRop3ERfkXfeeQeiKGLBggUYMmQIGjRogJ49e5Y5m9y3b1/07t0bERER8n18jx07hn379iEoKEhh\nx4g3oan3efjwIVavXg1BEHDu3DmIoohNmzYhKSkJTZo0QXBwcJmve/ToEdq1awc/Pz906NABpqam\nSElJwc8//4xx48bJ+1Vl+zR/f3/ExcVh0qRJuHz5svzJcsXFxZg9e7ZC3/HjxyM5ORkPHjyQt8XE\nxOB///sfPDw8YGhoiP/85z84cOAAxo4di379+lXr+wOo8Wa5x48f46233lLX6YiIiOosUd8chdKv\ntV7Dm1B1CUF5/czNzXHo0CHMnz8fCQkJ8gdqfP/99+jfv7/SOSZPnoyMjAx89913yM/Ph6enJ6Ki\nohR2iOjRowd2796NxYsXIzIyEvr6+vD09MS8efPKvNmuLG5ubpg/fz7i4+Px+eefo6SkROGhFK+P\nZ9OmTYiIiMCuXbuwbds22NjYyB/29foYyvpeqLq9mKrvU1aN5ZHJZFi8eLG8vyAI8u3kbG1tyw3C\nRkZGCA4ORlJSEg4ePIgXL17Azs4O4eHhmDBhQrXGp6enh+3bt2PevHlYu3Ytnj59Cnd3d8TFxSk9\nAEgQBKVflNq3b48DBw7gyJEjeP78Odq2bYvY2FgMHz5cpe9FeQSZTFbur3znz59Heno6ACAkJASf\nfvopOnXqpNRPJpMhISEBDRs2xIkTJ96oINIcnX26GnR47C9yIP4ZpJNbCeruk+Xa4NrzYFjbaGdN\nIVFVlT5ZLiEhAQEBAdouh3RMhTPC+/btk6/7FQQB69evx/r168vsa2ZmhjVr1qi/QiIiIiIiDagw\nCI8ZMwZ9+/aFKIrw9vbG7Nmz8cEHHyj1MzY2RsuWLWFgUOu3JSYiIiIiAlBJELayspLv65aYmAgn\nJ6dqPb6OiIiIqDxV2dKMSJ1UnsLt3r27JusgIiIiHVT6mF4ibSg3CIeGhkIQBHzzzTfQ19dX2L6k\nPIIg4LvvvlNrgUREREREmlBuEP7111+hp6eHkpIS6Ovr49dff630owt+tEFEREREdUW5Qbh027Ty\nviYiIiIiqstUfsQyEREREVF9Uq39zh4/fgyZTFbmYxVtbW3fuCgiIiIiIk1TOQgXFhYiKioK33//\nfYV3d/LOTyIiIiKqC1QOwlOnTsWWLVswYMAAdOvWTeHZ30REREREdY3KQTgxMRGjRo1CTEyMJush\nIiIiIqosspcAAAAfxklEQVQRKt8sJwgCXF1dNVkLEREREVGNUTkI9+/fH0ePHtVgKURERERENUfl\nIDx16lT89ddf+OKLL3Dq1CncvXsXubm5Sn+IiIiIiOoCldcId+rUCcDLB2v88MMP5fbjrhFERERE\nVBeoHIRnzJjBRygTERERUb2hchCeNWuWJusgIiIiIqpRfMQyEREREekklWeEo6KiKu0jCAJmzJjx\nRgUREREREdUElYNwZGRkuccEQYAoigzCRERERFRnqByEHzx4oNRWUlKCrKwsxMfHIzk5GTt27FBr\ncUREREREmvJGa4T19PRgb2+P8PBwtG7dmrPBRERERFRnqO1mOU9PTxw+fFhdpyMiIiIi0ii1BeE/\n/vgDenrchIKIiIiI6gaV1whv2bKlzPb8/HwkJycjMTERo0aNUlthRERERESapHIQDgkJKfdYs2bN\nMGXKFK4RJiIiIqI6Q+UgfPbsWaU2QRBgZmaGxo0bq7UoIiIiIiJNU3lRr52dndIfW1vbGgnB2dnZ\nmDhxIhwcHGBlZYVu3bohOTlZoU9ERAScnZ3RvHlz+Pr64tKlSxqvi4iIiIjqrlp/d1t+fj58fHwg\nCAJ27NiBtLQ0REVFQSqVyvvExMQgLi4OS5cuRVJSEqRSKQICAlBQUKDFyomIiIioNlN5aYS2fPPN\nN2jevDlWrlwpb7Ozs1Pos2rVKkyZMgW+vr4AgLi4ODg6OmLHjh0YPXp0jdZLRERERHVDrZ8RPnDg\nADw8PBAUFARHR0f06NEDa9eulR/PzMxEdnY2evfuLW8zNDSEp6cnUlNTtVEyEREREdUBtT4IZ2Zm\nYt26dWjZsiV27dqFiRMnYsGCBYiPjwcA5OTkQBAEhaUSACCVSpGTk6ONkomIiIioDqj1SyNKSkrg\n4eGBuXPnAgBcXFxw9epVxMfH47PPPtNydURERERUV6kUhJ88eQIbGxvMnj0b06ZN03RNCiwtLdGm\nTRuFtjZt2mD16tUAAAsLC4iiiNzcXFhbW8v75ObmwsLCosJzZ2RkqL/gWk4Xx1xKF8du2RRoDOjs\njaM6Oe6Sp4C+bv57d3R01HYJRFTHqBSEjYyMYG5uDlNTU03Xo6Rr165KP9AzMjJga2sLALC3t4el\npSWSkpLg5uYGACgsLERKSgrCw8MrPLeu/dDMyMjQuTGX0tmxv8iBmAsYGxtru5IaV1BQoJPjLm7Q\nCHiuez/fiIiqQ+U1woMGDcLu3btRUlKiyXqUhISE4NSpU1i+fDn++usv7NmzB2vWrEFwcLC8z8SJ\nExETE4PExERcuHABISEhMDExweDBg2u0ViIiIiKqO1ReI+zr64vjx4+jb9++GDVqFOzt7dGoUSOl\nfh4eHmot8J133sGmTZuwYMECLFu2DDY2Npg7dy6CgoLkfb788ksUFhZixowZkMlk8PDwwK5du3Ry\nNoiIiIiIVCPIZDJRlY4SiUTxhYKg8LUoihAEAXl5eeqrjtRKZ5cHQIfH/iIH4p9BOvlLoe4ujWiD\na8+DYW1jV3lnIiIdp/KMcGxsrCbrICIiIiKqUSoH4REjRmiyDiIiIiKiGlWtB2pcvXoVJ0+eRH5+\nvrrrISIiIiKqEVUKwtu3b0eHDh3QqVMn9O/fH2fOnAEA3L9/Hx4eHti9e7dGiiQiIiIiUjeVg/De\nvXsxbtw4tGnTBgsXLoQo/v977Jo1a4Y2bdpg69atGimSiIiIiEjdVA7Cy5cvR69evbBr164y1wu/\n++67OH/+vFqLIyIiIiLSFJWD8JUrV+Dr61vucalUinv37qmlKCIiIiIiTVM5CBsZGaGgoKDc43/9\n9ReaNWumlqKIiIiIiDRN5SDs5eWFzZs349mzZ0rH7ty5g40bN8Lb21utxRERERERaYrK+wjPnTsX\nffr0Qa9evTBo0CAIgoAjR44gKSkJGzduhL6+PsLCwjRZKxERERGR2qg8I9y6dWscOnQIlpaWiIyM\nhCiKiI2NxTfffAMXFxccPHgQtra2mqyViIiIiEhtVJ4RBgAnJyfs3r0bMpkM165dQ0lJCezt7WFu\nbq6p+oiIiIiINKJKQbiUmZkZ3N3d1V0LEREREVGNqVIQlslkiI2NxaFDh5CVlQUAsLOzg4+PD0JD\nQ2FmZqaRIomIiIiI1E3lNcLXrl1D9+7dsWzZMrx48QI9evRAjx498OLFCyxbtgzvvfcerl69qsla\niYiIiIjURuUZ4enTp+Phw4fYu3cvvLy8FI4dO3YMgYGBCAsLw44dO9ReJBERERGRuqk8I5ySkoIJ\nEyYohWAA6NmzJ8aPH4/k5GS1FkdEREREpCkqzwg3adKkwjXAZmZmaNKkiVqKIs0ogQGuZt7Rdhla\nUVBYopNjt2hShEYlwJOnRdoupcYV6+i4gefQ01d5joOISKepHIQDAwPxww8/IDAwEI0bN1Y4lp+f\njx9++AGjRo1Se4GkPnn5BYhJ2KrtMrSioKAAxsbG2i6jxgUPeQdWBTehb6Cv7VJqXPGLYp0cd1Mr\nS+i/rXvjJiKqDpWDsKOjIwRBwLvvvovhw4ejVatWAICrV69i69atkEqlcHR0xO7duxVeFxAQoN6K\niYiIiIjUQOUgPG7cOPnfv/nmG6XjOTk5GDduHERRlLcJgsAgTERERES1kspBODExUZN1EBERERHV\nKJWDcPfu3TVZBxERERFRjeKtxURERESkkxiEiYiIiEgnMQgTERERkU5iECYiIiIincQgTEREREQ6\nSeUgHBUVhQsXLpR7/OLFi4iKilJLUUREREREmqZyEI6MjMT//ve/co8zCBMRERFRXaK2pRGPHz/G\nW2+9pa7TERERERFpVIUP1Dh//jzS09PlX6ekpODFixdK/WQyGRISEuDo6Kj+ComIiIiINKDCILxv\n3z75cgdBELB+/XqsX7++zL5mZmZYs2aN+iskIiIiItKACoPwmDFj0LdvX4iiCG9vb8yePRsffPCB\nUj9jY2O0bNkSBgYqP7GZiIiIiEirKkyuVlZWsLKyAgAkJibCyckJUqm0RgojIiIiItIklW+W6969\ne60IwdHR0ZBIJJgxY4ZCe0REBJydndG8eXP4+vri0qVLWqqQiIiIiOqCKq1l+Pnnn/H9998jMzMT\nMpkMoigqHBcEAWfOnFFrga/6/fffsXHjRnTo0EGhPSYmBnFxcVi5ciUcHBwQFRWFgIAAnDp1CsbG\nxhqrh4iIiIjqLpWD8D//+U/Mnz8fFhYWcHd3R7t27TRZl5L8/HyMGzcOsbGxiIyMVDi2atUqTJky\nBb6+vgCAuLg4ODo6YseOHRg9enSN1klEREREdYPKQXjVqlXw8vLC9u3btbJf8OTJkxEQEIDu3bsr\ntGdmZiI7Oxu9e/eWtxkaGsLT0xOpqakMwkRERERUJpWDsEwmg7+/v1ZC8MaNG5GZmYl169YpHcvJ\nyYEgCErrl6VSKe7evVtTJRIRERFRHaNyEPbw8EBGRoYmaynTn3/+iUWLFuHQoUPQ01Pbg/AAQCvj\n0baCggJtl6A1ujj25//3AJziF8VarkQ7dHHcJcUi9KGbP9/4UCciqiqVg/CyZcswdOhQuLm54eOP\nP9ZkTQrS0tKQl5eHLl26yNuKi4uRnJyM9evXIyUlBaIoIjc3F9bW1vI+ubm5sLCwqPDcuvZD896p\n8zp782BBQYFOjv0tAwOgCNA30Nd2KTWu+EWxTo5bT18AoHs/34iIqkPlIDxq1Cg8e/YMEyZMwJQp\nU9C8eXPo6yv+T0YQBJw8eVKtBfr6+sLd3V2hLSQkBA4ODpg6dSocHBxgaWmJpKQkuLm5AQAKCwuR\nkpKC8PBwtdZCRERERPWHykHY3NwcUqkUDg4OmqxHiampKUxNTRXajIyMYGZmBicnJwDAxIkTER0d\nDQcHB7Ru3RrLli2DiYkJBg8eXKO1EhEREVHdoXIQ3r9/vybrqBJBEBS+/vLLL1FYWIgZM2ZAJpPB\nw8MDu3bt0smPwomIiIhINVV6oEZtkZiYqNQWFhaGsLAwLVRDRERERHVRlbZhyMvLQ3h4OHx8fODu\n7o60tDR5e1RUFC5fvqyRIomIiIiI1E3lGeHr16+jX79+yMvLQ7t27ZCZmYmnT58CAJo2bYpdu3bh\n3r17WLp0qcaKJSIiIiJSF5WD8Lx58yCKIk6ePInGjRsr3TTXv3//WrWOmIiIiIioIiovjTh69CiC\ng4Nhb2+vdLMaALRo0QK3b99Wa3FERERERJqichAuKiqCmZlZucfz8/PV/uQ3IiIiIiJNUTm5Ojs7\n47fffiv3+P79+9GxY0e1FEVEREREpGkqB+GJEydi9+7dWLZsGR48eAAAKCkpwZUrV/DZZ5/h1KlT\nCA0N1VihRERERETqpPLNckOHDsXNmzexePFiLF68GADkT27T09PDggUL0K9fP81USURERESkZlV6\noMaUKVMwdOhQ/PTTT7h27RpKSkrQsmVL+Pn5wd7eXkMlEhERERGpX5WfLGdjY4OQkBBN1EJERERE\nVGNUXiN88uRJREdHl3t8xYoV8ifNERERERHVdirPCEdFRVW4fdr58+dx4sQJ7Ny5Uy2FERERERFp\nksozwufOnUPnzp3LPd6pUyecPXtWLUUREREREWmaykH4yZMnZT5R7lWPHz9+44KIiIiIiGqCykHY\nwcEBv/zyS7nH//Of/6BVq1ZqKYqIiIiISNNUDsKjRo3CkSNHMGPGDPkDNQAgLy8P06dPxy+//ILA\nwECNFElEREREpG4q3ywXHByM9PR0rF27FvHx8bCwsAAA5OTkQBRFjBgxAhMnTtRYoURERERE6lSl\nfYT/+c9/yh+okZmZCQCwt7eHv78/unfvron6iIiIiIg0QqUg/OzZM/z++++wsrJCjx490KNHD03X\nRURERESkUSqtETYwMMCgQYMqvFmOiIiIiKguUSkI6+npwc7OjtujEREREVG9ofKuERMmTMCGDRuQ\nm5uryXqIiIiIiGqEyjfLPXnyBEZGRnB3d8eAAQNgb2+PRo0aKfQRBAFffPGF2oskIiIiIlI3lYPw\n/Pnz5X/ftm1bmX0YhImIiIiorlA5CJ89e1aTdRARERER1SiVg7CdnZ0m6yAiIiIiqlFVeqAGAFy9\nehUnTpxAbm4uhg4dihYtWuDZs2fIzs6GpaUlGjRooIk6iYiIiIjUSuUgXFJSgilTpuD777+HKIoQ\nBAGdOnWSB+H33nsP06dPx6RJkzRZLxERERGRWqi8fdry5cvxww8/YM6cOThy5AhEUZQfMzExgZ+f\nH/bt26eRIomIiIiI1E3lILxp0yaMHDkSU6dORatWrZSOt2vXDlevXlVrcUREREREmqJyEL59+zY8\nPDzKPd6oUSM+eY6IiIiI6gyVg7CFhQWysrLKPX7mzBnY2tqqpSgiIiIiIk1TOQgPHDgQCQkJCssf\nBEEAABw5cgRbt27FoEGD1F8hEREREZEGqByEZ86cCRsbG3h5eSE4OBiCICA6Ohrvv/8+hg0bhg4d\nOuDvf/+7JmslIiIiIlIblYOwqakpDh8+jL///e/IycmBoaEhTp48iYKCAsycORMHDhxAo0aN1F5g\ndHQ0vL29YWdnBwcHB3zyySe4ePGiUr+IiAg4OzujefPm8PX1xaVLl9ReCxERERHVHyoHYQAwNDTE\n1KlTcfz4cdy+fRt3795FSkoKZsyYAUNDQ40UmJycjODgYBw+fBiJiYkwMDDAoEGDIJPJ5H1iYmIQ\nFxeHpUuXIikpCVKpFAEBASgoKNBITURERERU91X6QI3CwkIcOHAA169fR9OmTeHj4wMrK6uaqA0A\nsGPHDoWvV69eDTs7O6SmpsLHxwcAsGrVKkyZMgW+vr4AgLi4ODg6OmLHjh0YPXp0jdVKRERERHVH\nhUH4zp076N+/P65fvy5/gIaRkRG2bt2KHj161EiBr3v06BFKSkpgZmYGAMjMzER2djZ69+4t72No\naAhPT0+kpqYyCBMRERFRmSpcGhEeHo6srCyEhIRg27ZtiIiIgKGhIcLCwmqqPiUzZ86Eq6srOnfu\nDADIycmBIAiQSqUK/aRSKXJycrRRIhERERHVARXOCB89ehTDhw9HeHi4vM3CwgKfffYZbt26BWtr\na40X+KrZs2cjLS0NBw8elG/dRkRERERUHRUG4ezsbHTp0kWhrWvXrhBFETdv3qzRIDxr1izs2bMH\n+/btg52dnbzdwsICoigiNzdXoZ7c3FxYWFhUeM6MjAyN1Vtb6fINhLo49ucvXgAAil8Ua7kS7dDF\ncZcUi9CHbv58c3R01HYJRFTHVBiEi4uLlXaDKP26sLBQc1W9JiwsDHv37sW+ffvQunVrhWP29vaw\ntLREUlIS3Nzc5LWlpKQozGSXRdd+aN47dR7GxsbaLkMrCgoKdHLsbxkYAEWAvoG+tkupccUvinVy\n3Hr6Lz8t07Wfb0RE1VHprhGZmZk4ffq0/OuHDx8CeDnbYGJiotTfw8NDjeUB06ZNw48//ohNmzbB\n1NRUvu7X2NhYHmwmTpyI6OhoODg4oHXr1li2bBlMTEwwePBgtdZCRERERPWHIJPJxPIOSiSSMtfi\niqKo1F7alpeXp9YCy6shLCxM4aa9qKgobNiwATKZDB4eHli2bBnatm2r1lrqupRT5xGTkKjtMrRC\nV2eEg4e8A6uCmTo5M6qrM8JNrd6F8PZ8WFg213YpRES1XoUzwrGxsTVVR7kePHigUr/XgzERERER\nUUUqDMIjRoyoqTqIiIiIiGpUlR6xTERERERUXzAIExEREZFOYhAmIiIiIp3EIExEREREOolBmIiI\niIh0EoMwEREREekkBmEiIiIi0kkMwkRERESkkxiEiYiIiEgnMQgTERERkU5iECYiIiIincQgTERE\nREQ6iUGYiIiIiHQSgzARERER6SQGYSIiIiLSSQzCRERERKSTDLRdANUcG6keQj7pqO0ytOLFi+cw\nMHhL22XUONvmEjy5pO0qiIiIaicGYR1i2vABLAq+1nYZWlH8ohj6BvraLqPGFf2p7QqIiIhqLy6N\nICIiIiKdxCBMRERERDqJQZiIiIiIdBKDMBERERHpJAZhIiIiItJJDMJEREREpJMYhImIiIhIJzEI\nExEREZFOYhAmIiIiIp3EIExEREREOolBmIiIiIh0EoMwEREREekkBmEiIiIi0kkMwkRERESkkxiE\niYiIiEgnMQgTERERkU6qV0E4Pj4erq6usLKyQq9evZCSkqLtkoiIiIiolqo3QXjXrl2YNWsWpk2b\nhuPHj6Nz584YOnQobt26pe3SiIiIiKgWqjdBeOXKlRg5ciQCAwPh6OiIJUuWwNLSEgkJCdoujYiI\niIhqoXoRhJ8/f44zZ86gV69eCu3e3t5ITU3VTlFEREREVKvViyB8//59FBcXw8LCQqFdKpUiJydH\nS1URERERUW1moO0CtCkjI0PbJdSoZk2kaO6xVttlEJEG6eu/BVnBM537+QYAjo6O2i6BiOqYehGE\nmzVrBn19faXZ39zcXKVZ4lfp2g/NjIwMODq6a7sMrXg5dt263qUiIo5i1qxe2i7jjb1VsA7Pjceq\n3F+Xr/ndPN0dOxFRVdSLpRFvvfUW3NzccPToUYX2pKQkdO3aVTtFEdUCx4/fxsqV1/Dbb7e1Xcob\n0S86jkaPFkK/6Ddtl0JERPVIvQjCABAaGorNmzfjX//6F65cuYKwsDBkZ2djzJgx2i6NSGuiov6L\nR4+KERHxX22X8kYaPo6CIOaj4eMIbZdCRET1SL1YGgEAAQEBePDgAZYvX47s7Gw4Oztj+/btsLGx\n0XZpRFpx/PhtpKffAwCkp9/Hb7/dxnvvva3lqqpOv+g49J+nv/z783ToF/2G4obvabkqIiKqD+pN\nEAaAoKAgBAUFabsMolohKuq/yM9/DgDIz3+GiIj/Yt++uheEGz6Ogp6YDwDQ+79Z4ScN92m5KiIi\nqg/qzdIIIvr/Xp0NLlU6K1yXvDobLG/7v1lhIiKiN8UgTFQPvTobXKp0VrgueXU2uJQe1woTEZGa\nMAgT1TNlzQaXqkuzwmXNBsuPcVaYiIjUoF6tESYi4MoVGTw9m0MQBABAQcFjGBubAABEUcSlS7I6\ncdOc3osrKG7giWIIZRwVoffiEm+aIyKiN8IgTFTPjB3bDmPHtpN/XVcfLPHceGyVHqBBRERUVYJM\nJhO1XQQRERERUU3jGmEiIiIi0kkMwkRERESkkxiEiYiIiEgnMQgTERERkU5iECYiIiIinVTvg/DG\njRvh5+eHFi1aQCKR4MaNGyq9bu/evejatSssLS3RrVs37Nu3T8OVqt+zZ88wffp0tG7dGtbW1hg+\nfDhu3674YQqbN2+GRCJB06ZNIZFI5H9/9uxZDVVdPfHx8XB1dYWVlRV69eqFlJSUCvtfuHABAwYM\nQPPmzdG+fXssWbKkhipVr6qMOysrS35NX722v/zySw1WrB7JyckYPnw42rVrB4lEgi1btlT6mvpw\nzas67vp0zaOjo+Ht7Q07Ozs4ODjgk08+wcWLFyt9XX247kSkOfU+CD958gR9+vTBrFmz5A8YqExa\nWhrGjh2LYcOG4cSJExgyZAjGjBmD//63bj2edubMmdi/fz8SEhLw73//G48ePcKwYcMgihXvmGds\nbIwrV67I/1y+fBkNGjSooaqrbteuXZg1axamTZuG48ePo3Pnzhg6dChu3bpVZv9Hjx4hICAAVlZW\nOHr0KCIiIvDtt98iNja2hit/M1UdNwAIgoDdu3crXFsvL68arFo9CgoK0L59e0RGRsLIyKjS/vXl\nmld13ED9uebJyckIDg7G4cOHkZiYCAMDAwwaNAgymazc19SX605EmqMz+wifOXMG3t7eOHv2LGxt\nbSvsGxQUBJlMhl27dsnbBg0aBKlUirVr12q6VLV4+PAhHBwcEBcXh8GDBwMAbt26BRcXF+zcuRO9\ne/cu83WbN29GWFiYyjPntcH7778PFxcXrFixQt7m4eGBQYMGYe7cuUr9161bhwULFuDPP/+UB/xl\ny5Zh/fr1+N///ldjdb+pqo47KysLrq6uSEpKgpubW02WqlE2NjZYunQphg8fXm6f+nLNX6XKuOvr\nNQde/lJgZ2eHzZs3w8fHp8w+9fG6E5F61fsZ4epIS0tTCop9+vRBamqqliqqujNnzuDFixcK47C2\ntoaTk1Ol43j69ClcXFzQvn17DBs2DOfOndN0udX2/PlznDlzBr169VJo9/b2Lnecv//+O7p166Yw\ny92nTx/cuXMHWVlZmixXbaoz7lKBgYFwdHRE3759sXfvXg1WWXvUh2v+JurjNX/06BFKSkpgZmZW\nbh9dv+5EVDkG4TLk5OTAwsJCoU0qlSInJ0dLFVVdTk4O9PX10bRpU4X2ysbh6OiI7777Dlu2bMG6\ndetgaGiIvn374q+//tJ0ydVy//59FBcXV+l6lXd9RVGsM9e4OuM2MTFBeHg4NmzYgO3bt8PLywtB\nQUHYvn17TZSsVfXhmldHfb7mM2fOhKurKzp37lxuH1297kSkOgNtF1Ad4eHhWL58ebnHBUFAYmIi\n3nvvvRqsqmaoOvbq6tSpEzp16iT/unPnzujRowdWr16NyMjIap+XtK9p06YIDQ2Vf+3m5oYHDx7g\nm2++wdChQ7VYGWlKfb3ms2fPRlpaGg4ePKjyvR9ERGWpk0E4NDQUn3zySYV9bGxsqn1+CwsLpdmC\n3NxcpZkFbVB17GlpaSguLkZeXp7CrHBubi48PT1Vfj89PT24ubnh2rVr1a5Zk5o1awZ9ff0qXa/y\nrq8gCLXiGquiOuMui7u7OzZt2qTu8mqd+nDN1aWuX/NZs2Zhz5492LdvH+zs7Crsy+tORJWpk0sj\nJBIJHBwcKvxjaGhY7fN37twZR48eVWhLSkpCly5d3rDyN6fq2N3c3GBgYICkpCT5a2/duoXLly+j\na9euVXrP8+fPw9LSUt1DUYu33noLbm5uZV6v8sbZuXNnpKSkKGwJ98svv6B58+aV/o+1tqjOuMty\n7ty5Wntt1ak+XHN1qcvXPCwsDLt370ZiYiJat25daX9edyKqTJ0MwlWRk5OD9PR0ZGRkQBRFXLp0\nCenp6Qpb7gwcOBCLFi2Sfz1hwgT8+uuviImJQUZGBqKjo3HixAmEhIRoYwjVYmpqisDAQMybNw/H\njh3D2bNnMWHCBLi4uKBnz57yfq+PPSoqCr/88gsyMzORnp6O0NBQXLx4EWPHjtXGMFQSGhqKzZs3\n41//+heuXLmCsLAwZGdn49NPPwUALFiwAP7+/vL+Q4YMgZGREUJCQnDx4kX89NNP+OabbxQ+Qq4L\nqjruLVu2YMeOHbhy5Qr+/PNPfPvtt0hISMD48eO1NYRqKygoQHp6Os6dO4eSkhLcvHkT6enpuHnz\nJoD6e82rOu76dM2nTZuGLVu2YO3atTA1NUVOTg5ycnJQUFAg71NfrzsRaU6dXBpRFQkJCYiKioIg\nCBAEAcOGDQMAxMbGyrcdun79usLsQOfOnbFu3Tr84x//QEREBFq2bIn169fjnXfe0coYqisyMhIG\nBgYICgpCYWEhevbsidWrVyusqXt97Pn5+Zg8eTJycnJgamqKjh074t///net3nopICAADx48wPLl\ny5GdnQ1nZ2ds374d1tbWAIDs7Gxcv35d3t/U1BS7d+/GtGnT4O3tDTMzM0yaNKlO/aIDVH3cwMut\no27evAk9PT04ODggNjYWQ4YM0Ub5b+SPP/6An5+f/N9yREQEIiIiMHz4cMTGxtbba17VcQP155qv\nW7cOgiAoBF3g5SxxWFgYgPr73zoRaY7O7CNMRERERPSqer80goiIiIioLAzCRERERKSTGISJiIiI\nSCcxCBMRERGRTmIQJiIiIiKdxCBMRERERDqJQZiIiIiIdBKDMBERERHpJAZhIiIiItJJ/w9Vb0bn\n1t+1HwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1146e7588>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# HIDDEN \n",
    "pop_50 = make_array(1, 1, 1, 1, 1, 0, 0, 0, 0, 0)\n",
    "pop_90 = make_array(1, 1, 1, 1, 1, 1, 1, 1, 1, 0)\n",
    "\n",
    "coins = Table().with_columns(\n",
    "   \"Proportion of 1's: 0.5\", pop_50,\n",
    "   \"Proportion of 1's: 0.9\", pop_90,\n",
    ")\n",
    "coins.hist(bins=np.arange(-0.5, 1.6, 1))\n",
    "plots.scatter(0.5, -0.02, marker='^', color='darkblue', s=60)\n",
    "plots.scatter(0.9, -0.02, marker='^', color='gold', s=60)\n",
    "plots.ylim(-0.05, 1);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Remember that the possible values in the population are only 0 and 1.\n",
    "\n",
    "The blue histogram (50% 1's and 50% 0's) has more spread than the gold. The mean is 0.5. Half the deviations from mean are equal to 0.5 and the other half equal to -0.5, so the SD is 0.5.\n",
    "\n",
    "In the gold histogram, all of the area is being squished up around 1, leading to less spread. 90% of the deviations are small: 0.1. The other 10% are -0.9 which is large, but overall the spread is smaller than in the blue histogram.\n",
    "\n",
    "The same observation would hold if we varied the proportion of 1's or let the proportion of 0's be larger than the proportion of 1's. Let's check this by calculating the SDs of populations of 10 elements that only consist of 0's and 1's, in varying proportions. The function `np.ones` is useful for this. It takes a positive integer as its argument and returns an array consisting of that many 1's."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table border=\"1\" class=\"dataframe\">\n",
       "    <thead>\n",
       "        <tr>\n",
       "            <th>Population Proportion of 1's</th> <th>Population SD</th>\n",
       "        </tr>\n",
       "    </thead>\n",
       "    <tbody>\n",
       "        <tr>\n",
       "            <td>0.1                         </td> <td>0.3          </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "        <tr>\n",
       "            <td>0.2                         </td> <td>0.4          </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "        <tr>\n",
       "            <td>0.3                         </td> <td>0.458258     </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "        <tr>\n",
       "            <td>0.4                         </td> <td>0.489898     </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "        <tr>\n",
       "            <td>0.5                         </td> <td>0.5          </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "        <tr>\n",
       "            <td>0.6                         </td> <td>0.489898     </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "        <tr>\n",
       "            <td>0.7                         </td> <td>0.458258     </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "        <tr>\n",
       "            <td>0.8                         </td> <td>0.4          </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "        <tr>\n",
       "            <td>0.9                         </td> <td>0.3          </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "</table>"
      ],
      "text/plain": [
       "Population Proportion of 1's | Population SD\n",
       "0.1                          | 0.3\n",
       "0.2                          | 0.4\n",
       "0.3                          | 0.458258\n",
       "0.4                          | 0.489898\n",
       "0.5                          | 0.5\n",
       "0.6                          | 0.489898\n",
       "0.7                          | 0.458258\n",
       "0.8                          | 0.4\n",
       "0.9                          | 0.3"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sd = make_array()\n",
    "for i in np.arange(1, 10, 1):\n",
    "    # Create an array of i 1's and (10-i) 0's\n",
    "    population = np.append(np.ones(i), 1-np.ones(10-i))\n",
    "    sd = np.append(sd, np.std(population))\n",
    "    \n",
    "zero_one_sds = Table().with_columns(\n",
    "    \"Population Proportion of 1's\", np.arange(0.1, 1, 0.1),\n",
    "    \"Population SD\", sd\n",
    ")\n",
    "\n",
    "zero_one_sds"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Not surprisingly, the SD of a population with 10% 1's and 90% 0's is the same as that of a population with 90% 1's and 10% 0's. That's because you switch the bars of one histogram to get the other; there is no change in spread.\n",
    "\n",
    "More importantly for our purposes, the SD increases as the proportion of 1's increases, until the proportion of 1's is 0.5; then it starts to decrease symmetrically."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAFmCAYAAACYz2bqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlYVHX/PvB7WNw1dgYV3JgENyCMxh2Xx+xRQ1Jw1y/6\ne9RETQ01MDNckUe0xC1BUFyBFNFMLbcaCklzTVNxAxcciKUHs2IY5veHF1PTcHTAYQaY+3VdXJd8\n5nPOeb8BuTnnzDlHVFRUpAIREVEFzIxdABER1VwMCSIiEsSQICIiQQwJIiISxJAgIiJBDAkiIhLE\nkCAiIkFGC4nY2Fh4eHhALBbD19cX6enpgnOzs7NhbW2t8WFjY4OTJ0+q56SlpVU459atW4Zoh4io\nTrIwxkb379+P0NBQrFmzBlKpFDExMQgICEBGRgZatGhR4TIikQj79+9Hx44d1WPW1tZaczIyMmBl\nZaUes7Ozq54miIhMgFH2JDZu3Ihx48Zh/PjxkEgkiIyMhKOjI+Li4gSXUalUsLKygr29vfrDwkI7\n4+zs7DTmiESi6myFiKhOM3hIKBQKXLx4Eb6+vhrj/fr1Q0ZGxnOXLQ+VQYMGITU1Vet1lUoFX19f\nuLm5wc/PDzKZTJ+lExGZHIOHRH5+PpRKJRwcHDTG7e3tkZubW+EyTZo0wbJly7Bt2zYkJyejd+/e\nmDRpEpKTk9VzxGIx1q5di4SEBOzcuRMSiQR+fn44c+ZMtfZDRFSXGeWcRGXZ2NggODhY/bmnpycK\nCwvx6aefIiAgAADg6uoKV1dX9ZyuXbsiOzsb69atg1QqNXjNRER1gcH3JGxtbWFubq6115CXl6e1\nd/E8r732Gu7cufPcOd7e3i+cQ0REwgweEpaWlvD09MTp06c1xk+dOlWpv/gvX74MR0fHl55DRETC\njHK4KTg4GNOmTYOXlxekUim2bt0KuVyOoKAgAEB4eDjOnz+vPjm9Z88eWFpaokuXLjAzM8ORI0cQ\nFxeH8PBw9To3bdoEFxcXuLu7o6SkBImJiThy5Ah27NhhjBZrtMzMTEgkEmOXYRSm2rup9g2Ydu/6\nYJSQ8Pf3R2FhIaKioiCXy+Hu7o7k5GT1NRJyuRxZWVkay6xevRoPHjyAmZkZXF1dsWHDBowYMUL9\nukKhwOLFi/Ho0SM0aNAAbm5uSE5ORv/+/Q3aGxFRXSLik+lMjyn/ZWWqvZtq34Bp964PvHcTEREJ\nYkgQEZEghgQREQliSBARkSCGBBERCWJIEBGRIIYEEREJYkgQEZEghgQREQliSBARkSCGBBERCWJI\nEBGRIIYEEREJYkgQEZEghgQREQliSBARkSCGBBERCWJIEBGRIIYEEREJYkgQEZEghgQREQliSBAR\nkSCGBBERCWJIEBGRIIYEEREJYkgQEZEghgQREQliSBARkSCGBBERCWJIEBGRIIYEEREJYkgQEZEg\nhgQREQliSBARkSCGBBERCWJIEBGRIIYEEREJYkgQEZEghgQREQkyWkjExsbCw8MDYrEYvr6+SE9P\nF5ybnZ0Na2trjQ8bGxucPHlSY15aWhp8fX0hFovh5eWF+Pj46m6DiKhOM0pI7N+/H6GhoQgJCYFM\nJoOPjw8CAgLw8OFDwWVEIhFSUlJw8+ZN3Lx5Ezdu3EDv3r3Vr2dlZWHkyJGQSqWQyWSYM2cO5s+f\nj0OHDhmiJSKiOskoIbFx40aMGzcO48ePh0QiQWRkJBwdHREXFye4jEqlgpWVFezt7dUfFhYW6tfj\n4uLg5OSEiIgISCQSTJgwAaNHj8b69esN0RIRUZ1k8JBQKBS4ePEifH19Ncb79euHjIyM5y5bHiqD\nBg1Camqqxmtnz55F3759Ncb69++PCxcuQKlU6qV2IiJTY/CQyM/Ph1KphIODg8a4vb09cnNzK1ym\nSZMmWLZsGbZt24bk5GT07t0bkyZNQnJysnpObm5uhessLS1Ffn6+/hsh0pE8rxDvfbQZ49/7L977\naDNyfykydklEOrN48RTjs7GxQXBwsPpzT09PFBYW4tNPP0VAQIARKyN6sRXRiXiUWwAzMxEe5RZg\nRfRefBI+zdhlEenE4CFha2sLc3Nzrb2GvLw8rT2B53nttdewa9cu9ecODg4VrtPCwgK2traC68nM\nzNR5m3WJqfYNGL737Ic5UJT+dcgz68GfRvn683tuWiQSiV7WY/CQsLS0hKenJ06fPg0/Pz/1+KlT\npzBs2DCd13P58mU4OjqqP/fx8cHhw4c15pw8eRJeXl4wNzcXXI++vpC1SWZmpkn2DRind5cWTuo9\nibIyFZo72Bi8Bn7PTbN3fTDKu5uCg4Oxe/duJCQk4ObNm1iwYAHkcjmCgoIAAOHh4RoBsmfPHnz+\n+ee4efMmbt26hejoaMTFxWHq1KnqOUFBQcjJyUFoaChu3ryJhIQE7N27FzNnzjR4f0R/t3DWKLRw\ntEE9Sws0d7DBwlmjjF0Skc6Mck7C398fhYWFiIqKglwuh7u7O5KTk9GiRQsAgFwuR1ZWlsYyq1ev\nxoMHD2BmZgZXV1ds2LABI0aMUL/eqlUrJCUlISwsDPHx8RCLxYiMjMSQIUMM2hvRPznYWfEcBNVa\noqKiIpWxiyDDMuXdb1Pt3VT7Bky7d33gvZuIiEgQQ4KIiAQxJIiISFCtuJiO6GXJ8wqxIjoR2Q9z\n4NLCCQtnjYKDnZWxy6p2pto36Q/3JMgklF/1rChVqq96NgWm2jfpD0OCTELBr8UwMxMBAMzMRMgv\nKjZyRYZhqn2T/jAkyCTYvNIUZWXP3u1dVqaCzStNjVyRYZhq36Q/DAkyCeVXPVtamJvUVc+m2jfp\nDy+mM0GmfHGRqfZuqn0Dpt27PnBPgoiIBDEkiIhIEEOCiIgEMSSIiEgQQ4KIiAQxJIiISBBDgoiI\nBDEkiIhIEEOCiIgEMSSIiEgQQ4KIiAQxJIiISBBDgoiIBDEkiIhIEEOCiIgEMSSIiEgQQ4KIiAQx\nJIiISBBDgoiIBDEkiIhIEEOCiIgEMSSIiEgQQ4KIiAQxJIiISBBDgoiIBFkYuwAyDfK8QqyITkTB\nr8WweaUpFs4aBQc7K2OXRdWM3/faj3sSZBArohPxKLcAJYpSPMotwIrovcYuiQyA3/fajyFBBlHw\nazHMzEQAADMzEfKLio1cERkCv++1H0OCDMLmlaYoK1MBAMrKVLB5pamRKyJD4Pe99mNIkEEsnDUK\nLRxtUM/SAs0dbLBw1ihjl0QGwO977We0E9exsbGIjo6GXC6Hm5sbVq5ciW7dur1wudu3b6NPnz4Q\niUS4f/++ejwtLQ1Dhw7VmCsSifDDDz/A1dVV7/VT5TjYWeGT8GnGLoMMjN/32k+nkPj222+RkZEB\nuVwOkUgER0dHSKVS9OzZs0ob3b9/P0JDQ7FmzRpIpVLExMQgICAAGRkZaNGiheByCoUCkydPRo8e\nPfD9999rvS4SiZCRkQErq7/ePWFnZ1elGomI6AUhcefOHUycOBFXr16FSqXSeE0kEsHDwwPx8fFo\n3bp1pTa6ceNGjBs3DuPHjwcAREZG4sSJE4iLi8OiRYsEl/voo4/QqVMndO/evcKQAJ6FgrW1daXq\nISKiigmek3jy5An8/f2Rk5ODFStW4Mcff4RcLsfjx4/x448/YtmyZcjOzsY777yD3377TecNKhQK\nXLx4Eb6+vhrj/fr1Q0ZGhuByx44dw9dff43IyEjBOSqVCr6+vnBzc4Ofnx9kMpnOdRERkTbBkNix\nYwcKCwvx9ddfY9q0aWjbti3q1auH+vXro23btpg+fTq++uor/PLLL9i5c6fOG8zPz4dSqYSDg4PG\nuL29PXJzcytcJicnB7Nnz0ZMTAwaNWpU4RyxWIy1a9ciISEBO3fuhEQigZ+fH86cOaNzbUREpEnw\ncNPRo0cxbtw4tGnTRnDhdu3aYdy4cfjyyy8xderUaikQAKZOnYrJkyfDy8sLALQOfQGAq6urxgnq\nrl27Ijs7G+vWrYNUKq222oiI6jLBkLh+/Tr+85//vHAFPXr0wL59+3TeoK2tLczNzbX2GvLy8rT2\nLsrJZDKkp6cjIiICwLOQKCsrg729PaKiojBhwoQKl/P29kZKSspz68nMzNS59rrEVPsGTLd3U+0b\nMM3eJRKJXtYjGBJFRUWwt7d/4Qrs7OxQVFSk8wYtLS3h6emJ06dPw8/PTz1+6tQpDBs2rMJl0tPT\nNT4/fPgw1qxZg5MnT0IsFgtu6/Lly3B0dHxuPfr6QtYmmZmZJtk3YLq9m2rfgGn3rg+CIVFSUgIL\nixe/Q9bc3BwKhaJSGw0ODsa0adPg5eUFqVSKrVu3Qi6XIygoCAAQHh6O8+fPIzU1FQDg5uamsfz5\n8+dhZmaG9u3bq8c2bdoEFxcXuLu7o6SkBImJiThy5Ah27NhRqdqIiOgvz02BI0eO4Nq1a89dQVZW\nVqU36u/vj8LCQkRFRUEul8Pd3R3JycnqayTkcnml16tQKLB48WI8evQIDRo0gJubG5KTk9G/f/9K\n10dERM+IioqKtM8CA5W61kAkEqGgoEBvRVH1MuXdb1Pt3VT7Bky7d30Q3JO4dOmSIesgIqIaSDAk\nXFxcDFkHERHVQFW6wd/169dx48YN9T2ciIiobhK84jo5ORlTpkzRGp89eza6d++OoKAg/Pvf/8Zb\nb71VqdtyEBFR7SEYEomJiRCJRBpjBw4cwPbt29G7d2/s2rULixcvxvnz57F27dpqL5SIiAxP8HDT\nzz//rHVH1uTkZDRp0gQJCQlo1qwZ3nrrLRQVFeHw4cP48MMPq71YIiIyLME9ifz8fDg7O2uMyWQy\n9OjRA82aNVOP9ezZs0rXShARUc0nGBLW1tYoLCxUf3716lUUFxeja9euGvMaNGigdViKiIjqBsGQ\n6NixIxITE9Wfp6SkQCQSYcCAARrzbt269cL7IxERUe0keE5i9uzZePvtt9G3b1/Y29vjxIkT6Nat\nGzw9PTXmffnll1pjRERUNwjuSfTs2RMJCQlo2LAhHjx4gNGjR2P79u0ac+RyOR4+fIi333672gsl\nIiLDe+7FdEOGDMGQIUMEX3d0dERaWpreiyIioppBcE+CiIiIIUFERIIYEkREJIghQUREghgSREQk\niCFBRESCdH6exL1795CSkoIHDx7gjz/+0HhNJBJh/fr1ei+OiIiMS6eQ+OKLLxAUFISysjLY29uj\nXr16Gq/z3k1ERHWTTiGxYsUK9OzZEzExMbCzs6vumoiIqIbQ6ZzEvXv3MHPmTAYEEZGJ0SkkJBIJ\nCgoKqrsWIiKqYXQKiSVLlmDNmjW4d+9eNZdDREQ1iU7nJCIiIlBQUAAfHx+0a9cOVlZWGq+LRCJ8\n+eWX1VIgEREZj04hYWZmBldX1+quhYiIahidQuLw4cPVXQcREdVAvOKaiIgE6XzF9ePHj7F+/Xp8\n9913KCwshLW1NXr16oXg4GA+45qIqI7SaU/i1q1b6NWrFz777DM0btwY3t7eaNy4MTZv3oxevXrh\n9u3b1V0nEREZgU57EosXL0bTpk1x/PhxtGrVSj2enZ2Nd955B4sXL8bOnTurrUgiIjIOnfYkZDIZ\nFi5cqBEQAODi4oIPPvgAMpmsWoojIiLj0ikkFAoFmjRpUuFrTZo0gUKh0GtRRERUM+gUEp07d8aW\nLVtQVlamMa5SqbB161Z07ty5WoojIiLj0umcxPz58zFy5Ej4+PjA398fYrEYcrkcqampuH37NpKS\nkqq7TiIiMgKdQmLAgAFITEzEsmXLEBUVBZVKBZFIBE9PTyQmJqJfv37VXScRERmBztdJDBgwAAMG\nDMDTp09RVFQEKysrNGrUqDprIyIiI9M5JMo1atSI4UBEZCIEQ2LVqlWYMGECnJycsGrVqueuRCQS\nYf78+XovjoiIjEswJCIiIjBgwAA4OTkhIiLiuSupSkjExsYiOjoacrkcbm5uWLlyJbp16/bC5W7f\nvo0+ffpAJBLh/v37Gq+lpaXhww8/xPXr1+Hk5IRZs2YhKCioUnUREdFfBN8CW1hYCG9vb/W/n/dR\n2afW7d+/H6GhoQgJCYFMJoOPjw8CAgLw8OHD5y6nUCgwefJk9OjRQ+u1rKwsjBw5ElKpFDKZDHPm\nzMH8+fNx6NChStVGRER/0ek6ifv37wteMFdaWqr1F/2LbNy4EePGjcP48eMhkUgQGRkJR0dHxMXF\nPXe5jz76CJ06dYKfn5/Wa3Fxceq9HolEggkTJmD06NFYv359pWojIqK/6BQSHh4euHz5coWv/fTT\nT/Dw8NB5gwqFAhcvXoSvr6/GeL9+/ZCRkSG43LFjx/D1118jMjKywtfPnj2Lvn37aoz1798fFy5c\ngFKp1Lm+ukyeV4j3PtqM+RE78d5Hm5H7S5GxSyKqNvx51w+dQkKlUgm+plAoYGam+2Mp8vPzoVQq\n4eDgoDFub2+P3NzcCpfJycnB7NmzERMTI/jOqtzc3ArXWVpaivz8fJ3rq8tWRCfiUW4BFKVKPMot\nwIrovcYuiaja8OddPwRPXBcVFaGo6K/kffToEWxtbTXm/P7779izZ0+1P09i6tSpmDx5Mry8vAA8\nP7QqKzMzU2/rqumyH+ZAUfpsr+r3358i68GfJtV/OVPsGTC9vk39510ikehlPYIhsXnzZqxatQoi\nkQgikQgTJ06scJ5KpUJoaKjOG7S1tYW5ubnWXkNeXp7WnkA5mUyG9PR09busVCoVysrKYG9vj6io\nKEyYMAEODg4VrtPCwkIr3P5OX1/I2sClhRMe5Rbg99+fomHDRmjuYGNS/QPPflGaWs+AafbNn3f9\nEAyJwYMHw8XFBSqVCjNmzEBISAjatGmjMad+/fpo3749OnXqpPMGLS0t4enpidOnT2ucgD516hSG\nDRtW4TLp6ekanx8+fBhr1qzByZMnIRaLAQA+Pj5az+I+efIkvLy8YG5urnN9ddnCWaOwInovsh78\nieYONlg4a5SxSyKqNvx51w/BkOjcubP67q4ikQiDBg2CjY2NXjYaHByMadOmwcvLC1KpFFu3boVc\nLldf0xAeHo7z588jNTUVAODm5qax/Pnz52FmZob27durx4KCghAbG4vQ0FAEBQXhzJkz2Lt3L7Zu\n3aqXmusCBzsrfBI+zST/qiTTw593/dDpthxjxozR60b9/f1RWFiIqKgoyOVyuLu7Izk5GS1atAAA\nyOVyZGVlVWqdrVq1QlJSEsLCwhAfHw+xWIzIyEgMGTJEr7UTEZkSUVFRkU5ngX/++WckJCTg1q1b\n+OOPPzRXIhLh4MGD1VIg6Z8p/2Vlqr2bat+AafeuDzrtSZw7d059juL27dvo2LEjioqK8ODBA7Ro\n0ULrXAUREdUNOl3gsGTJEgwdOhRnzpyBSqVCdHQ0rly5ggMHDkCpVCIkJKS66yQiIiPQKSSuXr2K\nwMBAiEQiAFA/xrRPnz4ICQnBkiVLqq9CIiIyGp1CQqFQoFGjRjAzM4O1tTUeP36sfs3V1RU///xz\ntRVIRETGo1NItGnTBg8ePAAAdOzYETt37kRZWRnKysqwa9cuwYvgiIiodtMpJAYNGoTvv/8eAPD+\n++/j+PHjcHZ2RuvWrfH5558jODi4WoskIiLj0OndTX+/7Yavry++/vprHDp0CE+fPsWAAQPQr1+/\naiuQiIiMp9LPuAae3Tq8MrcHJyKi2kn3e3wTEZHJEdyT6NKli/otry8iEolw8eJFvRVFREQ1g2BI\n9OjRQ+eQICKiukkwJDZt2mTIOoiIqAbiOQkiIhKk07ub9uzZ88I5o0ePfuliiIioZtEpJKZPn17h\n+N/PWTAkiIjqHp1C4tKlS1pjBQUFOHbsGD7//HNs2bJF74UREZHx6RQSLi4uFY55enpCpVJhw4YN\niI2N1XtxRERkXC994rpbt2746quv9FELERHVMC8dEufOnUPjxo31UQsREdUwOh1uWrVqldaYQqHA\ntWvX8NVXX+E///mP3gsjIiLj0ykkIiIitMbq168PZ2dnvP/++5g7d67eCyMiIuPTKSQKCwuruw4i\nIqqBeMU1EREJ0vl5EkqlEnv27MHZs2eRk5MDJycn+Pj4YNSoUTA3N6/OGomIyEh02pPIzs6GVCrF\nzJkzceLECeTl5eHEiROYMWMGunXrhuzs7Oquk4iIjECnkJg3bx6Ki4tx9OhR/PTTTzh16hR++ukn\nHDlyBP/73/8wf/786q6TiIiMQKeQkMlk+Oijj/DGG29ojEulUixatAgymaxaiiMiIuPSKSQaN24M\ne3v7Cl+zt7dHw4YN9VoUERHVDDqFRGBgIOLi4ip8LT4+HqNGjdJrUUREVDPo9O6mtm3bIjU1Fd27\nd8fQoUPh4OCA3NxcHDx4EMXFxfjXv/6FHTt2qOePHz++2gomIiLDERUVFaleNMna2lr3FYpEKCgo\neKmiqHplZmZCIpEYuwyjMNXeTbVvwLR714cqP0+CiIjqvio/T4KIiOo+na+4BoBr167hu+++Q2Fh\nIaytrdGzZ0+4u7tXV21ERGRkOoVEaWkppk+fjs8//xwq1V+nMEQiEUaMGIFNmzbx1hxERHWQTm+B\njYiIQEpKCsLCwnDp0iU8fvwYly5dQlhYGFJSUip83gQREdV+Ou1JJCUlISQkBCEhIeoxFxcXhISE\nQKlUYteuXQgLC6u2IomIyDh02pN4/Pix1i05yr3xxht4/PixXosiIqKaQaeQEIvFOHPmTIWvZWRk\nQCwWV3rDsbGx8PDwgFgshq+vL9LT0wXn3rhxA0OHDsWrr74KsVgMT09PLF26FAqFQj0nLS0N1tbW\nGh82Nja4detWpWsjIqJndDrcFBgYiKioKJiZmSEwMBCOjo6Qy+XYv38/oqKi8N5771Vqo/v370do\naCjWrFkDqVSKmJgYBAQEICMjAy1atNCaX69ePYwZMwZdunTBK6+8gp9++gmzZs2CUqnExx9/rJ4n\nEomQkZEBKysr9ZidnV2laiMior/oFBIffPAB7t27h5UrV2o871qlUmHEiBFYsGBBpTa6ceNGjBs3\nTn37jsjISJw4cQJxcXFYtGiR1vw2bdqgTZs26s9btmyJgICACvc+7OzsKnWFOBERCdMpJCwsLBAb\nG4u5c+ciPT1dfZ1E9+7dK32dhEKhwMWLFzFz5kyN8X79+iEjI0Onddy5cwcnTpzA4MGDNcZVKhV8\nfX3x559/on379ggJCUGvXr0qVR8REf1Fp5DIz89H48aN0aFDB3To0OGlNpifnw+lUgkHBweNcXt7\ne3zzzTfPXfbNN9/EpUuXUFJSgokTJ2rsdYjFYqxduxZeXl5QKBTYu3cv/Pz88OWXX0Iqlb5UzURE\npkowJJRKJSIjI7F582YUFxfD3NwcgwYNQnR0tMYxf0OKj4/HkydP8NNPP2HRokVYu3Yt5syZAwBw\ndXWFq6urem7Xrl2RnZ2NdevWMSSIiKpIMCTi4uIQGRmJXr16wcvLC3fv3sXhw4fRtGlTbNy4scob\ntLW1hbm5OXJzczXG8/LytPYu/ql58+YAgFdffRWlpaWYNWsW3nvvPZiZVfwmLW9vb6SkpDx3nZmZ\nmZWovu4w1b4B0+3dVPsGTLN3fd35VjAktm/fjokTJ+KTTz5Rj8XHx2PevHn45JNPUK9evSpt0NLS\nEp6enjh9+jT8/PzU46dOncKwYcN0Xo9SqVR/CIXE5cuX4ejo+Nz1mOIthE351smm2rup9g2Ydu/6\nIBgSWVlZWLZsmcaYv78/5s6di/v376Ndu3ZV3mhwcDCmTZsGLy8vSKVSbN26FXK5HEFBQQCA8PBw\nnD9/HqmpqQCAxMRENGjQAB06dIClpSUuXLiApUuXYtiwYbC0tAQAbNq0CS4uLnB3d0dJSQkSExNx\n5MgRjYchERFR5QiGxJMnT9C0aVONsfLPi4uLX2qj/v7+KCwsRFRUFORyOdzd3ZGcnKy+RkIulyMr\nK+uvIi0ssGbNGty9excqlQrOzs6YMmUK3n33XfUchUKBxYsX49GjR2jQoAHc3NyQnJyM/v37v1St\nRESmTPDJdNbW1khISEDnzp3VY0qlEl27dsXu3bu13vraunXrai2U9MeUd79NtXdT7Rsw7d714blv\ngZ04cWKF42PHjtUa4yNLiYjqHsGQ2LBhgyHrICKiGkgwJMaMGWPIOoiIqAbS6S6wRERkmhgSREQk\niCFBRESCGBJERCSIIUFERIIYEkREJIghQUREghgSREQkiCFBRESCGBJERCSIIUFERIIYEkREJIgh\nQUREghgSREQkiCFBRESCGBJERCSIIUFERIIYEkREJIghQUREghgSREQkiCFBRESCGBJERCSIIUFE\nRIIYEkREJIghQUREghgSREQkiCFBRESCGBJERCSIIUFERIIYEkREJIghQUREghgSREQkiCFBRESC\nGBJERCSIIUFERIIYEkREJMhoIREbGwsPDw+IxWL4+voiPT1dcO6NGzcwdOhQvPrqqxCLxfD09MTS\npUuhUCg05qWlpcHX1xdisRheXl6Ij4+v7jaIiOo0o4TE/v37ERoaipCQEMhkMvj4+CAgIAAPHz6s\ncH69evUwZswYpKSk4Ny5c4iIiMCOHTuwfPly9ZysrCyMHDkSUqkUMpkMc+bMwfz583Ho0CFDtUVE\nVOdYGGOjGzduxLhx4zB+/HgAQGRkJE6cOIG4uDgsWrRIa36bNm3Qpk0b9ectW7ZEQECAxt5HXFwc\nnJycEBERAQCQSCQ4d+4c1q9fj6FDh1ZzR0REdZPB9yQUCgUuXrwIX19fjfF+/fohIyNDp3XcuXMH\nJ06cQM+ePdVjZ8+eRd++fTXm9e/fHxcuXIBSqXzpuomITJHB9yTy8/OhVCrh4OCgMW5vb49vvvnm\nucu++eabuHTpEkpKSjBx4kSNvY7c3FytkLC3t0dpaSny8/O1tmcs8rxCrIhORMGvxbB5pSkWzhoF\nBzsrY5dFRHpUl/6f16p3N8XHx+Pbb79FbGwsvvrqK6xdu9bYJVXaiuhEPMotQImiFI9yC7Aieq+x\nSyIiPatL/88Nvidha2sLc3Nz5Obmaozn5eW98K/95s2bAwBeffVVlJaWYtasWXjvvfdgZmYGBweH\nCtdpYWFKcCoAAAAYaUlEQVQBW1tbwXVmZmZWsZOqyX6YA0XpX4e/sh78afAaAMP3XZOYau+m2jdg\nmv/PJRKJXtZj8JCwtLSEp6cnTp8+DT8/P/X4qVOnMGzYMJ3Xo1Qq1R9mZmbw8fHB4cOHNeacPHkS\nXl5eMDc3F1yPvr6QunJp4YRHuQUwMxOhrEyF5g42Bq8hMzPT4NusKUy1d1PtGzBO7zXh/7m+GOVw\nU3BwMHbv3o2EhATcvHkTCxYsgFwuR1BQEAAgPDxcI0ASExORmpqKzMxM3Lt3DykpKVi6dCmGDRsG\nS0tLAEBQUBBycnIQGhqKmzdvIiEhAXv37sXMmTON0aKghbNGoYWjDepZWqC5gw0Wzhpl7JKISM/q\n0v9zo7wF1t/fH4WFhYiKioJcLoe7uzuSk5PRokULAIBcLkdWVtZfRVpYYM2aNbh79y5UKhWcnZ0x\nZcoUvPvuu+o5rVq1QlJSEsLCwhAfHw+xWIzIyEgMGTLE4P09j4OdFT4Jn2bsMoioGtWl/+eioqIi\nlbGLIMPioQfT691U+wZMu3d9qFXvbiIiIsNiSBARkSCGBBERCWJIEBGRIIYEEREJYkgQEZEghgQR\nEQliSBARkSCGBBERCWJIEBGRIIYEEREJYkgQEZEghgQREQliSBARkSCGBBERCWJIEBGRIIYEEREJ\nYkgQEZEghgQREQliSBARkSCGBBERCWJIEBGRIIYEEREJYkgQEZEghgQREQliSBARkSCGBBERCWJI\nEBGRIIYEEREJYkgQEZEghgQREQliSBARkSCGBBERCWJIEBGRIIYEEREJYkgQEZEghgQREQliSBAR\nkSCGBBERCTJaSMTGxsLDwwNisRi+vr5IT08XnJuWloYxY8bAzc0NzZs3R48ePbBz506tOdbW1hof\nNjY2uHXrVnW3QkRUZ1kYY6P79+9HaGgo1qxZA6lUipiYGAQEBCAjIwMtWrTQmv/DDz+gY8eOmD17\nNhwdHXHixAnMnj0bDRs2xPDhw9XzRCIRMjIyYGVlpR6zs7MzSE9ERHWRUUJi48aNGDduHMaPHw8A\niIyMxIkTJxAXF4dFixZpzZ87d67G55MmTYJMJsPBgwc1QgJ4FgrW1tbVVzwRkQkx+OEmhUKBixcv\nwtfXV2O8X79+yMjI0Hk9xcXFGnsMAKBSqeDr6ws3Nzf4+flBJpPpo2QiIpNl8JDIz8+HUqmEg4OD\nxri9vT1yc3N1WsfRo0fx7bffIigoSD0mFouxdu1aJCQkYOfOnZBIJPDz88OZM2f0Wj8RkSkxyuGm\nl3HmzBlMmTIFkZGR8PT0VI+7urrC1dVV/XnXrl2RnZ2NdevWQSqVGqPUGksikRi7BKMx1d5NtW/A\ntHvXB4PvSdja2sLc3FxrryEvL09r7+Kf0tPTERgYiIULF+L//u//Xrgtb29v3Llz52XKJSIyaQYP\nCUtLS3h6euL06dMa46dOnXruX/zfffcdAgMDERoaiqlTp+q0rcuXL8PR0fFlyiUiMmlGOdwUHByM\nadOmwcvLC1KpFFu3boVcLlefYwgPD8f58+eRmpoKAJDJZBg1ahT+3//7fxg+fLh6L8Tc3By2trYA\ngE2bNsHFxQXu7u4oKSlBYmIijhw5gh07dhijRSKiOsEoIeHv74/CwkJERUVBLpfD3d0dycnJ6msk\n5HI5srKy1PP37NmD33//HdHR0YiOjlaPOzs749KlSwCevWtq8eLFePToERo0aAA3NzckJyejf//+\nhm2OiKgOERUVFamMXQQREdVMdfbeTZW57QcAXLt2DYMHD4aTkxM6duyIyMhIA1Wqf/q+5UltUdnv\nebnbt2+jZcuWcHZ2ruYKq09Vet+4cSN8fHzg6OgId3d3LFmyxACV6l9lez9x4gQGDhwIZ2dntGvX\nDmPGjMHt27cNVK1+fP/99xg9ejQ6dOgAa2tr7Nmz54XLVPV3XJ0MifLbfoSEhEAmk8HHxwcBAQF4\n+PBhhfOLi4vh7+8PsViM06dPY+XKlYiOjsaGDRsMXPnLq2zv5bc8SUhIQHp6OiZPnozZs2dj3759\nBq785VS273IKhQKTJ09Gjx49DFSp/lWl97CwMMTHx2PJkiX44YcfkJSUhO7duxuwav2obO9ZWVkY\nO3YsevToAZlMhtTUVPz5558IDAw0cOUv57fffkPHjh0RERGBRo0avXD+y/yOq5OHmwYMGIDOnTtj\n7dq16jFvb28MGzaswtt+bN26FeHh4bh16xbq1asHAFi9ejXi4+Nx9epVg9WtD5XtvSJBQUEoKyvD\n9u3bq6tMvatq36GhoSguLkb37t2xYMEC3L9/3xDl6lVle8/MzET37t2Rnp6ucW1RbVTZ3lNTUzF5\n8mTk5eVBJBIBePbGGD8/P9y+fbtW3tKnZcuW+O9//4vRo0cLznmZ33F1bk+iKrf9OHv2LLp166b+\n4gFA//79kZOTg+zs7OosV6+q85YnNVlV+z527Bi+/vrrWn1osSq9HzlyBG3atMFXX30FT09PdOnS\nBe+++y5++eUXA1SsP1Xp/bXXXoOlpSUSEhJQVlaG4uJi7N69G97e3rUyIHT1Mr/j6lxIVOW2H7m5\nuRXOV6lUOt8qpCaorlue1HRV6TsnJwezZ89GTEyMTrvrNVVVer937x6ys7ORkpKCzZs3Y8uWLcjM\nzHzuX6I1UVV6d3Z2xv79+7FixQo4ODigVatWuH79Ovbu3WuIko3mZX7H1bmQoKoTuuVJXTR16lRM\nnjwZXl5eAJ7dHNJUlJWVoaSkBFu2bIFUKoVUKsVnn32Gc+fO4fz588Yur1rl5uZi5syZGD16NE6d\nOoXDhw+jSZMmmDhxorFLq7HqXEhU5bYfDg4OFc4XiUQvvFVITWLIW57UJFXpWyaTYdWqVbCzs4Od\nnR1mzZqFJ0+ewN7eHgkJCYYoWy+q0rujoyMsLCzQpk0b9Vi7du1gbm5eq87JVKX3mJgYNG7cGB9/\n/DE6d+6Mbt264bPPPsN3331XqUOytc3L/I6rcyFRldt++Pj4ID09HSUlJeqxkydPwsnJCS4uLtVZ\nrl4Z8pYnNUlV+k5PT4dMJkNaWhrS0tIQFhaGRo0aIS0tDX5+fgaoWj+q0rtUKkVpaSnu3bunHrt7\n9y6USmWd/3n//fffYW5urjFmZvbs12BZWVm11FkTvMzvuDoXEsCz237s3r0bCQkJuHnzJhYsWKB1\n24+//yIYMWIEGjVqhOnTp+Pnn3/GwYMH8emnnyI4ONhYLVRZZXuXyWQIDAzEpEmT1Lc8yc3NRX5+\nvrFaqJLK9u3m5qbx4eTkBDMzM7Rv3x6vvPKKsdqoksr27uvrCw8PD8yYMQOXL1/GpUuXMGPGDPj4\n+KgPv9UWle194MCBuHTpEiIjI3Hnzh1cvHgRwcHBaNmyZa06xPrbb7/hypUruHz5MsrKyvDgwQNc\nuXIFDx48AKDf33G17lbhuqjsbT+aNWuGlJQUhISEoF+/frCyssLMmTMxffp0Y7VQZdVxy5PaoLJ9\n1yWV7V0kEiExMRELFizAkCFD0KBBA/Tt2xfLly83VgtVVtnee/fujdjYWHz66aeIjo5Gw4YN0bVr\nV+zbtw8NGzY0VhuVduHCBQwdOlT9Nt6VK1di5cqVGD16NDZs2KDX33F18joJIiLSjzp5uImIiPSD\nIUFERIIYEkREJIghQUREghgSREQkiCFBRESCGBJERCSIIVEL7d69G9bW1uoPZ2dn9OzZEzExMVAq\nlcYuD++++y66dOlS6eWuXLmCiIgIFBUVab1mbW2NVatW6aM8naWlpWl8nZ2cnCCVShEZGYk//vjD\noLXoy+HDhyt80Ex5r999950Rqnqxp0+fYsqUKZBIJLC2tkZYWJjg3PXr12PUqFFwc3N77s+Nrk90\nM3V18oprUyASibB9+3Y0b94cxcXFOHDgAObPn49ffvkFoaGhRq+t/ErQyrhy5QpWrVqFUaNGaT3P\n4vjx42jevLm+StSZSCRCZGQkvLy88PTpU5w8eRKrVq3C3bt3sWnTJoPX87IOHz6Mb7/9Vut2DJ6e\nnjh+/Djat29vpMqeLzY2FikpKdiwYQPatWsHR0dHwbk7duxAs2bNMGTIEMTFxRmwyrqJIVGLde7c\nGa1btwbw7H48d+7cwebNm40eElWlUqkEw8Xb29vA1TyjUqkgkUjU2+/Vqxfy8vKwe/durFy5UvDh\nTCUlJRoPeDG2F9XTpEkTo32NdXH9+nWIxWKdHjNafjdXpVKJrVu3VndpdR4PN9UhXl5eKC4uVt+c\nr7S0FMuWLUOXLl3g4OCALl26YNmyZSgtLVUvk52dDWtra2zduhULFy6ERCJB8+bNMXLkSK0nVlW0\n616+/It221euXIk+ffrAxcUF7dq1w9tvv41z586pX9+9ezdmzJih7sPa2ho2NjbqW1dXtO3jx49j\n4MCB6jtZjh07Frdu3dKYM3jwYLz11lv45ptv0KdPHzRv3hzdu3fHF198ocuXtELlN8G7e/euxjaO\nHj2K3r17QywWq385FRcXY968eXB3d4ejoyNef/11bNy4UWN95Yd6Dh48iOnTp6N169ZwcXHBlClT\nUFhYqDG3Mus7dOgQ3nvvPbi6uuLVV1/F9OnTsWfPHjx69Eh9CM3DwwPAsxs9VnS4acOGDXj99dfh\n4OAANzc3zJs3D8XFxRpzrK2tsXz5cnz22Wfw8PCAs7MzBg8ejOvXr+v09UxMTETPnj0hFovRrl07\nTJ06FXK5XGP9e/fuxYMHD9Q/F/o4LPbPP0jOnz8Pf39/tG3bFk5OTvD09MS8efNeeju1Hfck6pC7\nd+/C3NwcjRs3BgBMmzYNqampeP/99yGVSpGRkYHVq1cjOzsbW7Zs0Vh27dq16Ny5MzZu3Ii8vDws\nWbIEw4cPx5kzZ7RurVwVOTk5ePfdd9GyZUs8ffoUSUlJGDx4ME6fPg13d3cMGjQIISEhiIqKQkJC\ngvrQklgsrnB9x48fx8iRI+Hr64tt27bhyZMnWL58Od566y3IZDL1ciKRCHfv3kVoaCjmzp0LGxsb\nREdHIygoCGfPnlXviVVGeTiU3y1WJBLh1q1b+OCDDzBv3jy0bt0a1tbWUKlUCAwMxJUrVxAWFoYO\nHTrg2LFjWLhwIQoKCvDhhx9qrDcsLAy+vr6Ii4vD7du3sWTJEjx+/BgHDx4EgEqv74MPPsCAAQOw\nZcsW/PHHH+jQoQPy8/Nx4cIF7N27FyqVSr13UdEhwiVLlmDt2rWYMmUKBg0ahOvXr2P58uW4evUq\nvvzyS425SUlJkEgkWLVqFUpKSrBo0SKMHTsWZ8+eVd+KuyLbtm3DnDlzMGLECHz88cfIycnBkiVL\ncP78eXzzzTdo1KgRjh8/jpUrV+Lq1avYtWsXVCqVXg6LFRQUqP/922+/Yfjw4Xj99dexefNmNG7c\nGNnZ2fjhhx9eeju1HUOiFlMqlVAqlSguLkZKSgoOHz6Mf//732jQoAGuXbuGffv2ITQ0FPPnzwfw\n7JCUubk5VqxYgdmzZ6NDhw7qdTVr1kxjb6Bdu3YYNGgQ9uzZg3Hjxr10revWrVP/u6ysDP3798el\nS5eQkJCAlStXwsbGRv0QnL8fRhOybNkytGnTBsnJyepfQl27dkXXrl2xfv16LFu2TD23oKAAR48e\nVa+zS5cuaN++PVJSUjBnzpwX1q5SqaBUKtXnJLZt24YuXbqgbdu2Gts4cOAAOnbsqB47evQozpw5\ng02bNmHUqFEAnn0PfvvtN6xfvx7BwcEaz1Xu0KED1q9fDwDqO3VOmTIF3377LXr37o1jx45Van3e\n3t749NNPNXqxtbVFvXr18Nprrz2356KiImzYsAFjx45V78H17dsXtra2mDp1Ko4ePYpBgwap51ta\nWiIxMVH9B4VKpUJQUBB+/PFHvP766xVuo6ysDCtWrEDv3r0RExOjHpdIJHjrrbewc+dOTJkyBd7e\n3rCxsdGp7qrKzMzEr7/+io8//lj9/6JHjx617pGu1YGHm2oplUqFrl27ws7ODm3atMG8efMwcuRI\n9e2+09PTIRKJtI7hBgYGQqVSae2uDx06VOPzN954Ay1atMDZs2f1Uu/p06cxdOhQtG3bFra2trCz\ns8Pt27e1Dg/p4unTp7h8+TL8/f01/kpt1aoV3njjDa3e2rVrpxE6dnZ2sLe3V997/0Xeeecd2NnZ\nwcXFBUFBQejduzd27typMcfFxUUjIIBn3wNzc3OMGDFCYzwwMBAlJSVaf6X+82FHw4YNg5mZmfp7\n8P3331dqfYMHD9apv4qcPXsWCoUCAQEBGuPDhw+HhYWF1te4b9++GnucHTt2hEqleu7XODMzE3l5\neVr9SKVSODs7G/SdVm3btsUrr7yC2bNnIykpCQ8fPjTYtms67knUUiKRCLt27YKTkxOaNm0KZ2dn\njROT5cey//kukPLP/3msu6JHGNrb2yMnJ+ela7106RICAwMxYMAArF+/HmKxGGZmZpg5c2aV3kpa\nVFQElUpV4aEoR0dH/PjjjxpjFZ1crlevns7bjoqKgpeXFxo2bAgXF5cKnztQUS2FhYWwtraGhYXm\nfzNHR0eoVKoXfg8sLS1hZWWl/h4UFRVVan1Ch+p0Ub6uf67D3NwcNjY2Wtv659e4/GfxeV9joW0A\nz3r65zaqU7NmzXDo0CH897//RUhICIqLi+Hu7o4PPvgAb7/9tsHqqIkYErWYu7u74GGZ8sMOcrlc\nY075CcG/H5YAoPX8W+DZM3D/fr1D/fr1oVAoNOb8/biukEOHDsHS0hI7d+7U+Mu/qKhI8N1Bz2Nl\nZQWRSKRxcrOcXC7X6u1ltW3btkpPLbO2tkZhYSFKS0s1frHr+j1QKBQoKiqCk5NTldb3MsrPqcjl\nco3j/0qlEgUFBXrZ1t9/Rv9JLpcb/Cl5nTp1wvbt21FWVoYLFy5gzZo1mDRpEtLS0uDm5mbQWmoS\nHm6qo7p37w6VSoX9+/drjCclJUEkEqFnz54a4+UnR8udOXMGDx8+hI+Pj3rM2dkZ165d05h37Nix\nF14T8fTpU62T3998843WoYj69esDePYc4udp1KgRPD09kZqaCpXqr2dmlZ9o7NWr13OXr4yqXO9R\nrkePHlAqlThw4IDGeFJSEurXr6/xtQWgNS8lJQUqlUo9rzLrE6q7fv36L/z6AsDrr7+OevXqaf38\n7Nu3D0qlUuvnpyokEgkcHBy0tpGRkYH79+/r9ftYGWZmZvD29kZYWBiUSiVu3LhhlDpqCu5J1FHu\n7u4YMWIEIiIioFAo8MYbb6jf3TRixAi4u7trzC8uLsbo0aMRFBSEvLw8LF26FBKJRH2CFHh2bD4q\nKgpRUVHo2rUr0tPTsW/fvhfWMmDAAGzevBnTpk1Tv0119erV6kdMlmvfvj1UKhViYmIwevRoWFpa\nolOnTlqHVwBg4cKFGDlyJAIDAzF58mQ8efIEERERsLKy0uuzyf8eQpX1r3/9C1KpFHPnzkVeXh7c\n3d1x7Ngx7Ny5E3PnztX6a/znn39GcHAwhg8fjszMTCxfvhy9evVS/7KszPqE6m7fvj0KCwsRFxcH\nLy8v1K9fX32i9u/LWFlZYcaMGVi7di0aNmyIgQMH4vr161ixYgW6deuGN998s8pfl3JmZmYICwvD\nnDlzMGXKFIwcORIPHz7E8uXLIZFIMHbs2Cqt9+LFi8jKykJZWRkA4MaNG0hNTQUAvPnmm2jQoIHW\nMseOHcO2bdswePBgtGrVCr/99hs+++wzNG3aVCvMTQ1Dog7btGkTWrdujV27diEqKgpisRhz5szB\nggULtObOnTsXd+7cwfTp0/H777+jV69eiIyM1NgDmDt3Lv73v/+pnxE8cOBAbNmyBf3799da39//\nku3Xrx9WrVqFDRs24IsvvoC7uzs2b96M1atXa8zr1KkTQkNDsX37diQkJKCsrAyXLl2Cs7Oz1ls0\n+/fvj6SkJKxatQqTJk2CpaUlevXqhfDwcK3zMBX9Va3rVeG67kkIbSM5ORlLlizBunXrUFBQABcX\nF6xYsQLTpk3TmhsREYEjR45g0qRJKCsrw6BBgzSuDans+ioyYcIE/Pjjj1i6dCl+/fVXjWeZ/3OZ\nRYsWwc7ODvHx8YiLi4ONjQ1Gjx6Njz76SGtbQv2/yMSJE9GoUSOsW7cOY8eORePGjTFw4ECEh4dr\nnfvR9XuxZcsW7N27V73MgQMH1Htf5T9P/9SuXTs0atQIq1evhlwuR5MmTfDaa6/hwIED6sN9porP\nuDZx2dnZ8PDwwLp16zB+/Hhjl2OS0tLS8PbbbyMlJQV9+vQxdjlEGnhOgqgGeJnDWkTViSFBL3Vy\nlvSD3wOqqXi4iYiIBHFPgoiIBDEkiIhIEEOCiIgEMSSIiEgQQ4KIiAQxJIiISND/B6MzOyV/1aMk\nAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x11499f6d8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "zero_one_sds.scatter(\"Population Proportion of 1's\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Summary:** The SD of a population of 1's and 0's is at most 0.5. That's the value of the SD when 50% of the population is coded 1 and the other 50% are coded 0."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### The Sample Size ###\n",
    "We know that \n",
    "$$\n",
    "\\sqrt{\\mbox{sample size}} ~ \\ge ~ 4 \\times \\frac{\\mbox{SD of the 0-1 population}}{0.01}\n",
    "$$\n",
    "and that the SD of the 0-1 population is at most 0.5, regardless of the proportion of 1's in the population. So it is safe to take\n",
    "\n",
    "$$\n",
    "\\sqrt{\\mbox{sample size}} ~ \\ge ~ 4 \\times \\frac{0.5}{0.01} ~=~ 200\n",
    "$$\n",
    "\n",
    "So the sample size should be at least $200^2 = 40,000$. That's an enormous sample! But that's what you need if you want to guarantee great accuracy with high confidence no matter what the population looks like."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [Root]",
   "language": "python",
   "name": "Python [Root]"
  },
  "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
}