{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "# HIDDEN\n", "\n", "from datascience import *\n", "path_data = '../../../data/'\n", "import matplotlib\n", "matplotlib.use('Agg', warn=False)\n", "%matplotlib inline\n", "import matplotlib.pyplot as plots\n", "plots.style.use('fivethirtyeight')\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Sampling from a Population ###\n", "\n", "The law of averages also holds when the random sample is drawn from individuals in a large population.\n", "\n", "As an example, we will study a population of flight delay times. The table `united` contains data for United Airlines domestic flights departing from San Francisco in the summer of 2015. The data are made publicly available by the [Bureau of Transportation Statistics](http://www.transtats.bts.gov/Fields.asp?Table_ID=293) in the United States Department of Transportation.\n", "\n", "There are 13,825 rows, each corresponding to a flight. The columns are the date of the flight, the flight number, the destination airport code, and the departure delay time in minutes. Some delay times are negative; those flights left early." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "
Date | Flight Number | Destination | Delay | \n", "
---|---|---|---|
6/1/15 | 73 | HNL | 257 | \n", "
6/1/15 | 217 | EWR | 28 | \n", "
6/1/15 | 237 | STL | -3 | \n", "
6/1/15 | 250 | SAN | 0 | \n", "
6/1/15 | 267 | PHL | 64 | \n", "
6/1/15 | 273 | SEA | -6 | \n", "
6/1/15 | 278 | SEA | -8 | \n", "
6/1/15 | 292 | EWR | 12 | \n", "
6/1/15 | 300 | HNL | 20 | \n", "
6/1/15 | 317 | IND | -10 | \n", "
... (13815 rows omitted)
" ], "text/plain": [ "Date | Flight Number | Destination | Delay\n", "6/1/15 | 73 | HNL | 257\n", "6/1/15 | 217 | EWR | 28\n", "6/1/15 | 237 | STL | -3\n", "6/1/15 | 250 | SAN | 0\n", "6/1/15 | 267 | PHL | 64\n", "6/1/15 | 273 | SEA | -6\n", "6/1/15 | 278 | SEA | -8\n", "6/1/15 | 292 | EWR | 12\n", "6/1/15 | 300 | HNL | 20\n", "6/1/15 | 317 | IND | -10\n", "... (13815 rows omitted)" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "united = Table.read_table(path_data + 'united_summer2015.csv')\n", "united" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "One flight departed 16 minutes early, and one was 580 minutes late. The other delay times were almost all between -10 minutes and 200 minutes, as the histogram below shows." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "-16" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "united.column('Delay').min()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "580" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "united.column('Delay').max()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/choldgraf/anaconda/envs/textbook/lib/python3.6/site-packages/matplotlib/axes/_axes.py:6462: UserWarning: The 'normed' kwarg is deprecated, and has been replaced by the 'density' kwarg.\n", " warnings.warn(\"The 'normed' kwarg is deprecated, and has been \"\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbUAAAEcCAYAAABAuSr7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XtcTfn+P/DXLhKFna4uJZJLqTFFHJfcM8YxNePRRMZtkHDGkVzCzMiMU1Lul46R5JAvJqGR4TBuhTC+Y4Q5cyLlNl1cNtomjdq/P/q1v/ZU29rau91evZ6PR4+HvT5r7/V+78a8rLU+ay2JTCZTgIiISASM9F0AERGRtjDUiIhINBhqREQkGgw1IiISDYYaERGJBkONiIhEg6FGRESiwVAjIiLRYKhpUVZWlr5L0Bn2ZpjYm+ESc3+67I2hRkREosFQIyIi0WCoERGRaDDUiIhINBhqREQkGgw1IiISDYYaERGJBkONiIhEo4G+CzBUAX+LqLRs2d/99VAJERFV4J4aERGJBkONiIhEQ2+htmXLFvTu3Rv29vawt7fH0KFDcfTo0WrXz83NhVQqrfRz/PjxWqyaiIjqMr2dU2vVqhWWLl0KJycnlJWV4X/+538wduxYnDp1Cl27dq32ffv27VMZt7CwqI1yiYjIAOgt1EaMGKHy+osvvsDWrVtx6dIltaHWokUL2Nra6ro8IiIyQHXinFppaSn27dsHuVwOLy8vteuOGzcOHTp0wLBhw3Dw4MFaqpCIiAyBRCaTKfS18evXr8PHxwfFxcUwMzPDli1bMGzYsCrXffToEXbt2oVevXqhQYMGOHz4MFauXInY2FgEBASo3Y4unt0TEpFQadnqRRO1vh0iIlLl7Oxc7ZheQ62kpAT37t3Ds2fPcPDgQWzfvh2HDh2Ci4uLoPeHhobi/PnzOHfunI4rray669TUfdmGLCsri70ZIPZmuMTcny570+vhRxMTE7Rv3x7dunXDkiVL4Obmhk2bNgl+v6enJ7Kzs3VYIRERGZI6cU6tQllZGUpKSgSvn5mZyUkjRESkpLfZj+Hh4fDx8UHr1q1RVFSEpKQkpKenY+/evQCApUuX4vLly0hJSQEA7Nq1Cw0bNoS7uzuMjIxw5MgRxMXFITw8XF8tEBFRHaO3UMvPz0dQUBAKCgrQrFkzuLq6IikpCYMHDwYA5OXl4fbt2yrviYmJwd27d2FsbAwnJyds2LDhjZNEiIio/tBbqMXGxmo0HhgYiMDAQF2WREREBq5OnVMjIiKqCYYaERGJBkONiIhEg6FGRESiwVAjIiLRYKgREZFoMNSIiEg0GGpERCQaDDUiIhINhhoREYkGQ42IiESDoUZERKLBUCMiItHQ+C79t27dQnp6OgoLC+Hv74+2bduipKQE+fn5sLW1hYmJiS7qJCIieiPBoVZWVoaQkBDs2LEDCoUCEokEPXr0UIZanz59MG/ePHz22We6rJeIiKhagg8/rly5Ejt37sTixYtx7NgxKBQK5Zi5uTlGjhyJQ4cO6aRIIiIiIQSHWmJiIj755BOEhoaiffv2lcZdXFxw69YtrRZHRESkCcGh9uDBA3h6elY73rhxYxQVFWmlKCIiorchONRsbGxw586dasevXLkCe3t7rRRFRET0NgSH2gcffID4+HiVQ4wSiQQAcOzYMezevRt+fn7ar/BPtmzZgt69e8Pe3h729vYYOnQojh49qvPtEhFR3Sc41MLCwtCmTRt4e3tj6tSpkEgkWLVqFYYMGYKAgAB07doVc+bM0WWtAIBWrVph6dKlOH36NE6ePAlvb2+MHTsW165d0/m2iYiobhMcas2aNcO///1vzJkzBwUFBTA1NUVGRgbkcjnCwsJw+PBhNG7cWJe1AgBGjBiBoUOHon379ujQoQO++OILmJub49KlSzrfNhER1W0aXXxtamqK0NBQhIaG6qoejZSWluLAgQOQy+Xw8vLSdzlERKRngvfURo4cidOnT1c7fubMGYwcOVIrRb3J9evX0bp1a9jY2CAkJAQ7d+6Eq6trrWybiIjqLsF7aunp6Rg/fny14w8fPsTZs2e1UtSbODs7Iy0tDc+ePcPBgwcxffp0HDp0CC4uLlWun5WVpfUa5HJ5rW2rrmBvhom9GS4x91eT3pydnasd0/jej9W5f/8+zMzMtPVxapmYmCgvAO/WrRv+93//F5s2bcKGDRuqXF/dF/C2qutVF9uqC7KystibAWJvhkvM/emyN7WhlpqaisOHDytfJyQk4NSpU5XWk8lkOH36tNqLs3WprKwMJSUletk2ERHVHWpD7ddff8XBgwcBlF+TdvnyZfz8888q60gkEjRp0gR9+vRBZGSk7ir9/8LDw+Hj44PWrVujqKgISUlJSE9Px969e3W+bSIiqtvUhtqcOXOU155ZWFhg/fr18Pf3r5XCqpOfn4+goCAUFBSgWbNmcHV1RVJSEgYPHqzXugAgJCKhysOSezYs0kM1RET1j+Bzak+ePNFlHYLFxsbquwQiIqqj+ORrIiISDcF7ahYWFsp7Parz+PHjGhVERET0tgSH2vz58yuFWmlpKe7cuYPDhw+jQ4cOGDZsmNYLJCIiEkpwqC1cuLDasby8PAwZMgQdOnTQSlFERERvQyvn1Ozs7PDpp58iOjpaGx9HRET0VrQ2UaRJkybIzc3V1scRERFpTCuhduPGDWzevBlOTk7a+DgiIqK3Ivicmru7e5WzH58+fYpnz56hSZMmSExM1GpxREREmhAcan369KkUahKJBFKpFO3atcOoUaNgYWGh9QKJiIiEEhxqvJMHERHVdbyjCBERiYZGz1OTyWRISkpCTk4OZDIZFAqFyrhEIqn2mWZERES6JjjUfvjhB0yYMAFyuRxNmzaFVCqttI6Q22gRERHpiuBQ+/zzz2FjY4MdO3bA1dVVlzURERG9FcHn1LKzszFt2jQGGhER1VmCQ83JyQlFRUW6rIWIiKhGBIfa4sWLER8fj5ycHB2WQ0RE9PYEn1M7ceIELCws0LNnT3h7e6N169YwNjZWWUcikSAmJkbrRRIREQkhONTi4+OVfz5+/HiV6zDUiIhInwSH2pMnT3RZBxERUY3p7Y4iq1atwsCBA2Fvbw8nJycEBATgxo0bat+Tm5sLqVRa6ae6PUciIqpfNLqjiDalp6dj8uTJ8PDwgEKhQEREBPz8/HDhwoU33hh537596Nq1q/I1b6RMRESAmlBzd3eHkZERLl26hIYNG1b76JnXSSQSXLlyRdCGk5OTVV5v3rwZDg4OyMjIwPDhw9W+t0WLFrC1tRW0HSIiqj+qDbWKR80YGRmpvNaVoqIilJWVVXn7rT8bN24ciouL4eTkhBkzZsDX11dndRERkeGoNtT+/KgZXT96JiwsDG5ubvDy8qp2HXNzc3z99dfo1asXGjRogMOHD2PSpEmIjY1FQECATusjIqK6TyKTyRRvXk23Fi1ahOTkZBw5cgSOjo4avTc0NBTnz5/HuXPnql0nKyurhhVWFhKRIHjd1Ysman37RET1lbOzc7VjGk8UycnJqfbRMwDw4YcfavR5CxcuRHJyMr777juNAw0APD09kZiYqHYddV/A2zIzM6u0TC6XV7lcF9uvbVlZWaLooyrszTCJuTdA3P3psjfBoZaTk4OZM2fi/PnzAFBloEkkEo1CbcGCBdi/fz++++47dOzYUfD7XpeZmclJI0REBECDUJs5cyZ+/PFHhISEwNPTE82aNavRhufOnYs9e/Zg586dkEqlyM/PB1C+B2Rubg4AWLp0KS5fvoyUlBQAwK5du5QzMY2MjHDkyBHExcUhPDy8RrUQEZE4CA61y5cvY/bs2Vi4cKFWNhwXFwcAlWYuLliwQLmNvLw83L59W2U8JiYGd+/ehbGxMZycnLBhwwZOEiEiIgAahFrr1q3RvHlzrW1YJpO9cZ0/z7gMDAxEYGCg1mogIiJxEXybrFmzZmHHjh18phoREdVZgvfUJkyYgLKyMnh4eOCDDz5Aq1atqnz0zKxZs7ReJBERkRCCQy0zMxMrV65EYWEhtm7dWuU6DDUiItInwaEWEhKC58+fY82aNVqZ/UhERKRtgkPt+vXrWLhwISZMmKDLeoiIiN6a4Ikijo6OVV5wTUREVFcIDrWwsDBs2bIFubm5uqyHiIjorQk+/HjmzBk0b94cXl5e8Pb2RuvWrauc/RgTE6P1IomIiIQQHGrx8fHKPx8/frzKdRhqRESkT4JD7cmTJ7qsg4iIqMYEn1MjIiKq6xhqREQkGgw1IiISDYYaERGJBkONiIhEg6FGRESiISjUXrx4gRYtWvAaNCIiqtMEhVqTJk1gZWXFO/MTEVGdJvjwo5+fH/bv34+ysjJd1kNERPTWBN9R5K9//SvS0tLw3nvvYfz48XB0dETjxo0rrefp6anVAomIiIQSHGq+vr7KP1+6dAkSiURlXKFQQCKR4PHjx9qrrgqrVq3Cd999h5s3b8LExATdu3fHkiVL4OLiotPtEhFR3Sc41DZu3KjLOgRLT0/H5MmT4eHhAYVCgYiICPj5+eHChQuwsLDQd3lERKRHgkMtMDBQl3UIlpycrPJ68+bNcHBwQEZGBoYPH66nqoiIqC54q+vUbt26hYyMDDx9+lTb9WisqKgIZWVlkEql+i6FiIj0TCKTyRRCV/7222+xdOlSPHjwAACwf/9+9O/fH48ePYKPjw8+//xzfPjhhzortioTJ07ErVu3cOrUqUoPLa2QlZWl9e2GRCQIXnf1oola3z4RUX3l7Oxc7Zjgw48HDx5EUFAQBg4ciODgYHzxxRfKMUtLS3Ts2BG7d++u1VBbtGgRMjIycOTIkWoDDVD/BbwtMzOzSsvkcnmVy3Wx/dqWlZUlij6qwt4Mk5h7A8Tdny57E3z4ceXKlRgwYACSk5OrPL/WvXt3XLt2TavFqbNw4ULs27cPKSkpcHR0rLXtEhFR3SU41P773//ir3/9a7Xj1tbWePjwoVaKepMFCxYoA61jx461sk0iIqr7BB9+bNKkCeRyebXjt2/fhqWlpVaKUmfu3LnYs2cPdu7cCalUivz8fADlhwPNzc11vn0iIqq7BO+peXt7Y9euXSgpKak09ttvv2H79u0YNGiQVourSlxcHJ4/fw5fX1906tRJ+bN+/Xqdb5uIiOo2wXtqX3zxBQYPHowBAwbAz88PEokEx44dw8mTJ7F9+3YYGxtjwYIFuqwVACCTyXS+DSIiMkyC99ScnJxw9OhR2NraYvny5VAoFNi4cSPWrl0LNzc3HDlyBPb29rqslYiISC3Be2oA0KlTJ+zfvx8ymQzZ2dkoKyuDo6MjrKysdFUfERGRYBqFWgWpVAoPDw9t10JERFQjGoWaTCbDxo0bcfToUdy5cwcA4ODggGHDhmHmzJm8VRUREemV4HNq2dnZ6Nu3L2JiYvDq1Sv069cP/fr1w6tXrxATE4M+ffrg1q1buqyViIhILcF7avPmzcOzZ89w8OBBeHt7q4ydPn0a48aNw4IFC5CUlKT1IomIiIQQvKd2/vx5BAcHVwo0AOjfvz+mTZuGc+fOabU4IiIiTQgOtebNm6s9ZyaVStG8eXOtFEVERPQ2BIfauHHjsHPnTjx//rzS2NOnT7Fz506MHz9eq8URERFpQvA5NWdnZ0gkEnTv3h1jxoxB+/btAZQ/MHT37t2wtraGs7Mz9u/fr/K+2n6+GhER1V+CQy0oKEj557Vr11YaLygoQFBQEBSK/3vmqEQiYagREVGtERxq3333nS7rICIiqjHBoda3b19d1kFERFRjgieKEBER1XUMNSIiEg2GGhERiQZDjYiIRIOhRkREoiE41KKionDjxo1qx3/55RdERUVppSgiIqK3ITjUli9fjuvXr1c7zlAjIiJ909rhx6KiIjRs2FCj95w9exajR49Gly5dIJVKkZiYqHb93NxcSKXSSj/Hjx+vSelERCQSai++vnbtGjIzM5Wvz58/j1evXlVaTyaTIT4+Hs7OzhptXC6Xw8XFBWPGjEFwcLDg9+3btw9du3ZVvrawsNBou0REJE5qQ+3QoUPKQ4oSiQTbtm3Dtm3bqlxXKpXim2++0WjjPj4+8PHxAQDMmDFD8PtatGgBW1tbjbZFRETipzbUJk6ciPfeew8KhQKDBg3CokWLMHTo0ErrmZmZoV27dmjQQPBdt2pk3LhxKC4uhpOTE2bMmAFfX99a2S4REdVtalPIzs4OdnZ2AMpvaNypUydYW1vXSmFVMTc3x9dff41evXqhQYMGOHz4MCZNmoTY2FgEBARU+76srCyt1yKXywUv18X29UEsfVSFvRkmMfcGiLu/mvSm7lSXQd3Q2NLSEp999pny9bvvvovHjx9j7dq1akNN03N9QpiZmVVaJpfLq1yui+3XtqysLFH0URX2ZpjE3Bsg7v502ZtGxwt/+OEH7NixAzk5OZDJZCrPTgPKz7tduXJFqwW+iaen5xtnTRIRUf0gONTWrVuH8PBw2NjYwMPDAy4uLrqsS7DMzExOGiEiIgAahNo///lPeHt749tvv9X4erTqFBUVITs7GwBQVlaGe/fu4erVq7CwsIC9vT2WLl2Ky5cvIyUlBQCwa9cuNGzYEO7u7jAyMsKRI0cQFxeH8PBwrdRDRESGTXCoyWQy+Pr6ai3QAOCnn37CyJEjla8jIyMRGRmJMWPGIDY2Fnl5ebh9+7bKe2JiYnD37l0YGxvDyckJGzZsUHs+jYiI6g/Boebp6an1mTj9+vWDTCardjw2NlbldWBgIAIDA7VaAxERiYfg22TFxMTg0KFD2Lt3ry7rISIiemuC99TGjx+PkpISBAcHIyQkBC1btoSxsbHKOhKJBBkZGVovkoiISAjBoWZlZQVra2t06NBBl/UQERG9NcGhlpqaqss6iIiIaoxPviYiItHQKNQeP36MZcuWYdiwYfDw8MDFixeVy6OiovDrr7/qpEgiIiIhBB9+zM3NxfDhw/H48WO4uLggJycHv//+O4DyR8EkJyfj4cOHiI6O1lmxRERE6ggOtSVLlkChUCAjIwNNmzatNGHk/fff53k3IiLSK8GHH0+dOoWpU6fC0dEREomk0njbtm3x4MEDrRZHRESkCcGh9vLlS0il0mrHnz59CiMjzjshIiL9EZxCXbp0wdmzZ6sdT01Nhbu7u1aKIiIiehuCQ2369OnYv38/YmJi8OTJEwDld9b/73//iylTpuDHH3/EzJkzdVYoERHRmwieKOLv74979+4hIiICERERAIBRo0YBAIyMjLB06VIMHz5cN1USEREJoNGTr0NCQuDv74+UlBRkZ2ejrKwM7dq1w8iRI+Ho6KijEomIiITRKNQAoE2bNpgxY4YuaiEiIqoRwefUMjIysGrVqmrHV69erbzDCBERkT4I3lOLiopSO6X/2rVrSE9Px759+7RSGBERkaYE76ldvXoVXl5e1Y736NEDP//8s1aKIiIiehuCQ+3FixdV3knkdUVFRTUuiIiI6G0JDrUOHTrgxIkT1Y4fP34c7du310pR6pw9exajR49Gly5dIJVKkZiYqPNtEhGRYRAcauPHj8exY8cwf/585cXXQPljZ+bNm4cTJ05g3LhxOinydXK5HC4uLli+fDkaN26s8+0REZHhEDxRZOrUqcjMzMSWLVsQFxcHGxsbAEBBQQEUCgUCAwMxffp0nRVawcfHBz4+PgBgMJcWBPwtotKyPRsW6aESIiJx0+g6tXXr1ikvvs7JyQEAODo6wtfXF3379tVFfURERIIJCrWSkhJcunQJdnZ26NevH/r166fruoiIiDQmKNQaNGgAPz8/REREwMnJSdc1aV1WVpbWP1Mul2u0/M90UZOuGWLNQrE3wyTm3gBx91eT3pydnasdExRqRkZGcHBwMNgp++q+gLdlZmZWaZlcLq9yeVV0UZMuZWVlGVzNQrE3wyTm3gBx96fL3gTPfgwODkZCQgIKCwt1UggREVFNCZ4o8uLFCzRp0gQeHh4YMWIEHB0dK02pl0gkmDVrltaLfF1RURGys7MBlD/P7d69e7h69SosLCxgb2+v020TEVHdJjjUwsPDlX/es2dPlevURqj99NNPGDlypPJ1ZGQkIiMjMWbMGMTGxup020REVLcJDrW6cl/Hfv36QSaT6bsMIiKqgwSHmoODgy7rICIiqjGNHxJ669YtpKeno7CwEP7+/mjbti1KSkqQn58PW1tbmJiY6KJOIiKiNxIcamVlZQgJCcGOHTugUCggkUjQo0cPZaj16dMH8+bNw2effabLeomIiKoleEr/ypUrsXPnTixevBjHjh2DQqFQjpmbm2PkyJE4dOiQTookIiISQnCoJSYm4pNPPkFoaGiVj5hxcXHBrVu3tFocERGRJgSH2oMHD+Dp6VnteOPGjQ32jiNERCQOgkPNxsYGd+7cqXb8ypUrvPiZiIj0SnCoffDBB4iPj1c5xCiRSAAAx44dw+7du+Hn56f9ComIiAQSHGphYWFo06YNvL29MXXqVEgkEqxatQpDhgxBQEAAunbtijlz5uiyViIiIrUEh1qzZs3w73//G3PmzEFBQQFMTU2RkZEBuVyOsLAwHD58uNK9IImIiGqTRhdfm5qaIjQ0FKGhobqqh4iI6K29MdSKi4tx+PBh5ObmokWLFhg2bBjs7OxqozYiIiKNqA213377De+//z5yc3OVF1s3adIEu3fvRr9+/WqlQCIiIqHUnlNbtmwZ7ty5gxkzZmDPnj2IjIyEqakpFixYUFv1ERERCaZ2T+3UqVMYM2YMli1bplxmY2ODKVOm4P79+2jdurXOCyQiIhJK7Z5afn4+evbsqbKsV69eUCgUuHfvnk4LIyIi0pTaUCstLYWpqanKsorXxcXFuquKiIjoLbxx9mNOTg4uX76sfP3s2TMAQFZWFszNzSutr+7+kERERLr0xlCLjIxEZGRkpeXz589XeV3xjLXHjx9rrzoiIiINqA21jRs31lYdRERENaY21AIDA3VeQFxcHNatW4f8/Hx07twZkZGR6N27d5XrpqWlYeTIkZWWX7x4ER07dtR1qUREVMdpdJssbUtOTkZYWBhWrlyJXr16IS4uDv7+/sjIyFD7GJuMjAxYWFgoX1tZWdVGuVoV8LeIKpfv2bColishIhIPwTc01oWNGzciMDAQEyZMQKdOnRAdHQ1bW1vEx8erfZ+1tTVsbW2VP8bGxrVUMRER1WV6C7WSkhJcuXIFgwYNUlk+aNAgXLhwQe17BwwYgE6dOuGDDz7AmTNndFkmEREZEL0dfnz06BFKS0thbW2tstza2hoFBQVVvsfOzg6rVq2Ch4cHSkpKsGfPHvj6+iI1NbXa83BERFR/6PWcmqacnZ3h7OysfO3l5YU7d+5g3bp1akMtKytL67XI5XKNlguli1q1pS7XVlPszTCJuTdA3P3VpLfXc+DP9BZqlpaWMDY2RmFhocrywsJC2NjYCP4cT09PJCcnq11H3RfwtszMzCotk8vlVS7XhC5q1YasrKw6W1tNsTfDJObeAHH3p8ve9HZOzcTEBN26dcPJkydVlp88ebLS/SbVyczMhK2trbbLIyIiA6TXw48zZ87EtGnT4OnpiZ49eyI+Ph55eXmYNGkSAGDatGkAgM2bNwMANm3aBAcHB3Tp0gUlJSXYu3cvUlNT8a9//UtvPRARUd2h11D76KOP8PjxY0RHRyM/Px9dunTB3r174eDgAACVngTwxx9/4Msvv8SDBw9gamqqXN/Hx0cf5etEVdev8do1IiJh9D5RZMqUKZgyZUqVY6mpqSqv//73v+Pvf/97bZRFREQGSK8XXxMREWkTQ42IiESDoUZERKLBUCMiItFgqBERkWgw1IiISDQYakREJBoMNSIiEg29X3xNb8anZBMRCcM9NSIiEg2GGhERiQYPPxow3vyYiEgV99SIiEg0GGpERCQaDDUiIhINnlMTGU7/J6L6jKFWT3BSCRHVBzz8SEREosE9tXqsukOVVVn2d38dVkJEpB3cUyMiItEw2D21uLg4rFu3Dvn5+ejcuTMiIyPRu3dvfZclWiERCTAzM6u0nOfliKguMchQS05ORlhYGFauXIlevXohLi4O/v7+yMjIgL29vb7Lq1c0OYTJACQiXTPIUNu4cSMCAwMxYcIEAEB0dDR++OEHxMfHY8mSJVrdlib/0yb1OAOTiHRNIpPJFPouQhMlJSVo2bIltm7dCj8/P+XyuXPn4saNGzh8+LAeqyMiIn0yuIkijx49QmlpKaytrVWWW1tbo6CgQE9VERFRXWBwoUZERFQdgws1S0tLGBsbo7CwUGV5YWEhbGxs9FQVERHVBQYXaiYmJujWrRtOnjypsvzkyZPo2bOnnqoiIqK6wCBnP86cORPTpk2Dp6cnevbsifj4eOTl5WHSpEn6Lo2IiPTI4PbUAOCjjz5CZGQkoqOj0a9fP2RkZGDv3r1wcHDQW01xcXFwd3eHra0t+vfvj3PnzumtFqHOnj2L0aNHo0uXLpBKpUhMTFQZVygUiIyMROfOnWFnZ4cRI0bgl19+UVlHJpMhKCgIDg4OcHBwQFBQEGQyWW22UcmqVaswcOBA2Nvbw8nJCQEBAbhx44bKOoba25YtW9C7d2/Y29vD3t4eQ4cOxdGjR5XjhtpXVVatWgWpVIp58+Yplxlqf5GRkZBKpSo/HTt2VI4bal+vy8vLQ3BwMJycnGBra4uePXsiPT1dOV5bPRpkqAHAlClTkJmZiYKCApw+fRp9+vTRWy0VF4OHhobizJkz8PLygr+/P+7evau3moSQy+VwcXHB8uXL0bhx40rja9euxcaNGxEVFYUTJ07A2toaH374IZ4/f65cZ8qUKbh69SqSkpKQlJSEq1evYtq0abXZRiXp6emYPHkyjh49ipSUFDRo0AB+fn548uSJch1D7a1Vq1ZYunQpTp8+jZMnT8Lb2xtjx47FtWvXABhuX3926dIlJCQkwNXVVWW5Iffn7OyMX3/9Vfnz+j98DbkvoDyMhg0bBoVCgb179+LChQtYsWKFyiz12urR4K5Tq4sGDx4MV1dXrFu3TrnMw8MDvr6+Wr8YXFdat26NFStWYOzYsQDK/1XVuXNnTJ06FXNeb8jCAAANdElEQVTnzgUA/P7773B2dsbXX3+NSZMm4ddff0XPnj1x5MgR9OrVCwBw/vx5DB8+HJcuXYKzs7Pe+nldUVERHBwckJiYiOHDh4uqNwBwdHTEkiVLMHHiRFH09fTpU/Tv3x/r1q1DVFQUXFxcEB0dbdC/t8jISKSkpOD8+fOVxgy5rwpfffUVzp49q3LU4HW12aPB7qnVFSUlJbhy5QoGDRqksnzQoEG4cOGCnqqqudzcXOTn56v01bhxY/Tu3VvZ18WLF2Fubq4yQadXr14wMzOrU70XFRWhrKwMUqkUgHh6Ky0txb59+yCXy+Hl5SWavmbPng1fX194e3urLDf0/nJyctC5c2e4u7vj008/RU5ODgDD7wsAUlNT4enpiUmTJqFDhw7o27cvvvnmGygU5ftMtdmjQU4UqUvEejF4fn4+AFTZ12+//QYAKCgogKWlJSQSiXJcIpHAysqqTvUeFhYGNzc3eHl5ATD83q5fvw4fHx8UFxfDzMwMO3fuhKurq/IvvqH2BQDbt29HdnY2vvnmm0pjhvx76969OzZt2gRnZ2c8fPgQ0dHR8PHxQUZGhkH3VSEnJwdbt27FjBkzMHv2bGRmZmLBggUAgKCgoFrtkaFGorZo0SJkZGTgyJEjMDY21nc5WuHs7Iy0tDQ8e/YMBw8exPTp03Ho0CF9l1VjWVlZ+Oqrr3DkyBE0bNhQ3+Vo1dChQ1Ved+/eHd26dcOuXbvQo0cPPVWlPWVlZXj33XeVp1veeecdZGdnIy4uDkFBQbVaCw8/1pBYLwa3tbUFALV92djY4NGjR8pDDED5sfOHDx/Wid4XLlyIffv2ISUlBY6Ojsrlht6biYkJ2rdvj27dumHJkiVwc3PDpk2bDL6vixcv4tGjR+jVqxcsLS1haWmJs2fPIi4uDpaWlmjRogUAw+3vdebm5ujcuTOys7MN/vcGlP+d6tSpk8qyjh074t69e8pxoHZ6ZKjVkFgvBm/bti1sbW1V+iouLsb58+eVfXl5eaGoqAgXL15UrnPx4kXI5XK9975gwQJloL0+dRow/N7+rKysDCUlJQbf14gRI3Du3DmkpaUpf959912MGjUKaWlp6NChg0H397ri4mJkZWXB1tbW4H9vQPm5r5s3b6osu3nzpvJRYLXZo3FYWFh4DXohAE2bNkVkZCTs7OxgamqK6OhonDt3Dhs2bEDz5s31XV61ioqK8J///Af5+fnYsWMHXFxc0KxZM5SUlKB58+YoLS3FmjVr4OTkhNLSUixevBj5+flYs2YNGjVqBCsrK/z4449ISkqCm5sb7t+/j5CQEHh4eOh1qvHcuXOxe/duJCQkoE2bNpDL5ZDL5QDK/xEikUgMtrfw8HCYmJigrKwM9+/fR2xsLPbu3Yvw8HBlL4bYFwCYmprC2tpa5efbb7+Fg4MDxo4da9C/t88//1z5e7t58ybmzZuH7OxsrF69GlKp1GD7qtCmTRtERUXByMgIdnZ2OH36NJYtW4aQkBB4enrW6u+OU/q1JC4uDmvXrkV+fj66dOmCiIgIvV47J0RaWhpGjhxZafmYMWMQGxsLhUKB5cuXIyEhATKZDJ6enoiJiYGLi4tyXZlMhvnz5+P7778HAAwfPhwrVqxQzjTUh+q2vWDBAixcuBAADLa36dOnIy0tDQUFBWjWrBlcXV0xa9YsDB48GIDh9lWdESNGKKf0A4bb36effopz587h0aNHsLKyQvfu3bF48WJ07twZgOH29bqjR4/iq6++ws2bN9GmTRtMnToV06ZNU078qK0eGWpERCQaPKdGRESiwVAjIiLRYKgREZFoMNSIiEg0GGpERCQaDDUiIhINhhqJWmJiosqDGVu1agU3NzeMHTsW+/fvV7kljybS0tIglUqRlpam5YrfLCUlBc7Oznjx4oVWPq/iO8rNzdXK51UlLS0NkZGRKCsre6v3//zzz2jZsmWdf0Yh6R9DjeqF7du349ixY9i7dy8WL16MRo0aYfLkyfDz88Pvv/+u7/IEe/XqFb766ivMmjULTZo00cpnDhs2DMeOHYOdnZ1WPq8q6enpiIqKeutQe+eddzBgwAD84x//0HJlJDYMNaoX3Nzc0KNHD/Tt2xejR49GfHw8tm3bhjNnzuDLL7/Ud3mCpaam4s6dO/jkk0+09plWVlbo0aMHGjVqpLXP1IVJkyYhKSlJ+agSoqow1Kje8vX1xfvvv49//etfKofyXrx4gSVLlsDd3R3W1tZwd3dHTEzMG/cyTpw4AX9/f3Tq1AktW7bEX/7yF6xfvx6lpaXKdQICAtCvX79K783JyYGFhQXi4+PVbmPHjh0YPHgwLCwsVJZLpVIsW7YM69evR9euXdGyZUt8/PHHKCwsRGFhISZOnAgHBwe4urpizZo1Ku+t6vCjm5sbgoKCsG/fPnh5eaFVq1YYMGBApSc3jxgxAiNGjKhUp5ubG6ZPnw6g/KnPUVFRAMoDtOJQcAWh3/egQYPQtGlT7Nq1S+13RPUbn6dG9ZqPjw9SU1Px008/oU+fPnj16hVGjRqF//znP5g3bx5cXV1x6dIlREdH48mTJ2oPf+Xk5MDb2xtBQUFo1KgRrly5gqioKDx69Ajh4eEAgMmTJ+Pjjz/G5cuX4enpqXzv9u3bYWZmBn9//2o//+XLl0hPT8fixYurHN+9ezdcXFywcuVKFBQUYNGiRQgODkZRURGGDBmCiRMn4sCBAwgPD4eLiwt8fHzUfjfnzp1DVlaW8nDtP/7xDwQEBODq1asa3Ytv/PjxePDgAXbs2FHpuXaafN8NGjSAl5cXjh8/jtDQUMHbp/qFoUb1Wps2bQD831OVk5KScP78eaSmpipvSN2/f38AQFRUFGbPnl3p6b0VPv30U+WfFQoFevfujZKSEqxfvx5ffvkljIyMMGTIEDg6OmLbtm3KUPvjjz+QmJgIf39/NG3atNpaMzMzUVxcjK5du1Y53qhRI+zatQsNGpT/tf7ll1+wadMmLF68GPPmzQMA9O3bF4cOHcKBAwfeGGrPnz9Henq6MsBsbW0xcOBAHDt2TG34/lnr1q3RqlUrAOUPx6yoD9D8+3Zzc8P69etRVlYGIyMeaKLK+F8F1WsVsx8r7iT+ww8/wN7eHj179sSrV6+UP4MGDcIff/yBS5cuVftZeXl5mD17Nrp27Qpra2tYWVlh2bJlePr0qfLhiEZGRpg0aRKSk5Px9OlTAOXnyQoKCjBx4kS1tVacS7K0tKxyfODAgSqBUfEcuYo7+APlezvt27fH/fv31W4LKH++1et7ZBV3U6948KM2aPp9W1lZ4eXLl3jy5InWaiBx4Z4a1WsV/3N//cm8d+/ehZWVVZXrP378uMrlZWVlGDNmDPLy8hAWFgZnZ2c0btwYqampiImJQXFxsXLdcePGISIiAnv27EFQUBDi4+Ph6emJd955R22tL1++BIBqJ3T8+ZBgw4YNq11e8Vnq/Pm8XcV2X++lpjT9vhs3bgwABjVjlWoXQ43qtaNHj8LU1BTdunUDALRo0QJt27ZFQkJCles7ODhUufz27dv46aefsHnzZgQEBCiXVzwX6nUtWrSAn58fEhISMHjwYKSlpWHdunVvrLVFixYAyp85VVeYmpri+fPnlZYLrVHT77tiD626vVUihhrVWwcPHsT333+P4OBg5TVfgwcPRkpKCszMzJSH74SomD1ZsXcElJ8r+/bbb6tcf8qUKRg6dCg+++wzNGvWDKNGjXrjNpydnQGUT0jR5PH2umRvb4+UlBSUlJTAxMQEAHD27NlKQVexl/f777+rnDfU9PvOzc1FmzZtlHtsRH/GUKN6ITMzE48ePUJJSQnu3buHo0eP4sCBAxg4cCCWLFmiXO/jjz9GYmIifH19MXPmTLi5uaGkpAS3b9/G999/j8TExCoveu7UqRPs7e3x9ddfw9jYGA0aNMCmTZuqradHjx5wd3fHuXPnEBQUJOhCant7e9jb2+Py5csqe4P69NFHHyEhIQF/+9vfEBgYiNzcXGzcuBHNmjVTWa9Tp04AgA0bNmDo0KEwNjbGu+++q/H3/eOPP6J379612iMZFoYa1QsTJkwAUH64zMrKCu+88w7i4+Ph6+urnCQClO9pJScnY/Xq1di+fTtyc3PRpEkTtGvXDj4+Psq9kT8zMTFBYmIi5s+fj+DgYFhYWGDs2LGwt7fHrFmzqnyPn58frl69ikmTJgnu46OPPsKBAwewYsUKDbrXHW9vb6xevRrr169HSkoK3N3d8c0332DcuHEq67333nuYMmUKtm7dihUrVkChUEAmk2n0fd+7dw/Xrl2r9pIGIgCQyGSyt7v5HRHVyLBhw2BkZFTlebfq3L59G927d8ehQ4fwl7/8RYfV1T1r1qzB1q1bceXKFZVr3Yhexz01olr08uVL/Pzzzzh16hQuXLig8d0x2rVrh7Fjx2LNmjX1KtSKi4vxz3/+E0uWLGGgkVoMNaJalJeXBx8fHzRv3hyhoaF4//33Nf6MxYsXY9u2bXjx4oXWbmpc1925cwfBwcEYPXq0vkuhOo6HH4mISDR4RxEiIhINhhoREYkGQ42IiESDoUZERKLBUCMiItFgqBERkWj8PzwDCXEnQtk9AAAAAElFTkSuQmCC\n", "text/plain": [ "