{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "tags": [
     "remove_input"
    ]
   },
   "outputs": [],
   "source": [
    "path_data = '../../../../data/'\n",
    "\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "\n",
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "plt.style.use('fivethirtyeight')\n",
    "\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Iteration ###\n",
    "It is often the case in programming – especially when dealing with randomness – that we want to repeat a process multiple times. For example, recall the game of betting on one roll of a die with the following rules:\n",
    "\n",
    "- If the die shows 1 or 2 spots, my net gain is -1 dollar.\n",
    "- If the die shows 3 or 4 spots, my net gain is 0 dollars.\n",
    "- If the die shows 5 or 6 spots, my net gain is 1 dollar.\n",
    "\n",
    "The function `bet_on_one_roll` takes no argument. Each time it is called, it simulates one roll of a fair die and returns the net gain in dollars."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "def bet_on_one_roll():\n",
    "    \"\"\"Returns my net gain on one bet\"\"\"\n",
    "    x = np.random.choice(np.arange(1, 7))  # roll a die once and record the number of spots\n",
    "    if x <= 2:\n",
    "        return -1\n",
    "    elif x <= 4:\n",
    "        return 0\n",
    "    elif x <= 6:\n",
    "        return 1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Playing this game once is easy:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "-1"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bet_on_one_roll()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To get a sense of how variable the results are, we have to play the game over and over again. We could run the cell repeatedly, but that's tedious, and if we wanted to do it a thousand times or a million times, forget it.  \n",
    "\n",
    "A more automated solution is to use a `for` statement to loop over the contents of a sequence. This is called *iteration*. A `for` statement begins with the word `for`, followed by a name we want to give each item in the sequence, followed by the word `in`, and ending with an expression that evaluates to a sequence. The indented body of the `for` statement is executed once *for each item in that sequence*."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "cat\n",
      "dog\n",
      "rabbit\n"
     ]
    }
   ],
   "source": [
    "for animal in np.array(['cat', 'dog', 'rabbit']):\n",
    "    print(animal)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It is helpful to write code that exactly replicates a `for` statement, without using the `for` statement.  This is called ***unrolling* the loop**.  \n",
    "\n",
    "A `for` statement simple replicates the code inside it, but before each iteration, it assigns a new value from the given sequence to the name we chose.  For example, here is an unrolled version of the loop above."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "cat\n",
      "dog\n",
      "rabbit\n"
     ]
    }
   ],
   "source": [
    "animal = np.array(['cat', 'dog', 'rabbit']).item(0)\n",
    "print(animal)\n",
    "animal = np.array(['cat', 'dog', 'rabbit']).item(1)\n",
    "print(animal)\n",
    "animal = np.array(['cat', 'dog', 'rabbit']).item(2)\n",
    "print(animal)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Notice that the name `animal` is arbitrary, just like any name we assign with `=`.\n",
    "\n",
    "Here we use a `for` statement in a more realistic way: we print the results of betting five times on the die as described earlier. This is called *simulating* the results of five bets. We use the word *simulating* to remind ourselves that we are not physically rolling dice and exchanging money but using Python to mimic the process.\n",
    "\n",
    "To repeat a process `n` times, it is common to use the sequence `np.arange(n)` in the `for` statement. It is also common to use a very short name for each item. In our code we will use the name `i` to remind ourselves that it refers to an item."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1\n",
      "-1\n",
      "0\n",
      "1\n",
      "0\n"
     ]
    }
   ],
   "source": [
    "for i in np.arange(5):\n",
    "    print(bet_on_one_roll())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this case, we simply perform exactly the same (random) action several times, so the code in the body of our `for` statement does not actually refer to `i`."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Augmenting Arrays ###\n",
    "\n",
    "While the `for` statement above does simulate the results of five bets, the results are simply printed and are not in a form that we can use for computation. Any array of results would be more useful. Thus a typical use of a `for` statement is to create an array of results, by augmenting the array each time.\n",
    "\n",
    "The `append` method in `NumPy` helps us do this. The call `np.append(array_name, value)` evaluates to a new array that is `array_name` augmented by `value`. When you use `append`, keep in mind that all the entries of an array must have the same type."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['Cat', 'Dog', 'Another Pet'], dtype='<U11')"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pets = np.array(['Cat', 'Dog'])\n",
    "np.append(pets, 'Another Pet')\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This keeps the array `pets` unchanged:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['Cat', 'Dog'], dtype='<U3')"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pets"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "But often while using `for` loops it will be convenient to mutate an array – that is, change it – when augmenting it. This is done by assigning the augmented array to the same name as the original."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['Cat', 'Dog', 'Another Pet'], dtype='<U11')"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pets = np.append(pets, 'Another Pet')\n",
    "pets"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Example: Betting on 5 Rolls ###\n",
    "We can now simulate five bets on the die and collect the results in an array that we will call the *collection array*. We will start out by creating an empty array for this, and then append the outcome of each bet. Notice that the body of the `for` loop contains two statements. Both statements are executed for each item in the given sequence."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 0.,  0., -1.,  1., -1.])"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "outcomes = np.array([])\n",
    "\n",
    "for i in np.arange(5):\n",
    "    outcome_of_bet = bet_on_one_roll()\n",
    "    outcomes = np.append(outcomes, outcome_of_bet)\n",
    "    \n",
    "outcomes"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let us rewrite the cell with the `for` statement unrolled:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([-1.,  1.,  1.,  1.,  0.])"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "outcomes = np.array([])\n",
    "\n",
    "i = np.arange(5).item(0)\n",
    "outcome_of_bet = bet_on_one_roll()\n",
    "outcomes = np.append(outcomes, outcome_of_bet)\n",
    "\n",
    "i = np.arange(5).item(1)\n",
    "outcome_of_bet = bet_on_one_roll()\n",
    "outcomes = np.append(outcomes, outcome_of_bet)\n",
    "\n",
    "i = np.arange(5).item(2)\n",
    "outcome_of_bet = bet_on_one_roll()\n",
    "outcomes = np.append(outcomes, outcome_of_bet)\n",
    "\n",
    "i = np.arange(5).item(3)\n",
    "outcome_of_bet = bet_on_one_roll()\n",
    "outcomes = np.append(outcomes, outcome_of_bet)\n",
    "\n",
    "i = np.arange(5).item(4)\n",
    "outcome_of_bet = bet_on_one_roll()\n",
    "outcomes = np.append(outcomes, outcome_of_bet)\n",
    "\n",
    "outcomes"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The contents of the array are likely to be different from the array that we got by running the previous cell, but that is because of randomness in rolling the die. The process for creating the array is exactly the same.\n",
    "\n",
    "By capturing the results in an array we have given ourselves the ability to use array methods to do computations. For example, we can use `np.count_nonzero` to count the number of times money changed hands."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.count_nonzero(outcomes)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Example: Betting on 300 Rolls ###\n",
    "Iteration is a powerful technique. For example, we can see the variation in the results of 300 bets by running exactly the same code for 300 bets instead of five."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "outcomes = np.array([])\n",
    "\n",
    "for i in np.arange(300):\n",
    "    outcome_of_bet = bet_on_one_roll()\n",
    "    outcomes = np.append(outcomes, outcome_of_bet)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "The array `outcomes` contains the results of all 300 bets."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "300"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(outcomes)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To see how often the three different possible results appeared, we can use the array `outcomes` and `DataFrame` methods."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAikAAAEfCAYAAACEWFYzAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAiC0lEQVR4nO3deXSM9+LH8c+InUjQLLbYEksspWqtWkLVUlJLbb1XT6qE8Ku1lpYWV0lstbS0BK2l11ZKa6teLkFLF26caplYao1UmCB2md8fTud2bkQyWcyTmffrHOd0nu93Hp/v5JSPZxuTxWKxCgAAwGDyODsAAADAo1BSAACAIVFSAACAIVFSAACAIVFSAACAIVFSAACAIVFSAACAIVFSAACAIVFS3IjZbHZ2BKdw13VLrN1dufPa4VooKQAAwJAoKQAAwJAoKQAAwJAoKQAAwJDyOjsAAABGl5ycrPv37zs7hsvJmzevihQpkvb4E8wCAECuc+fOHUmSl5eXk5O4nuTkZN25c0cFChR45DinewAAeIzbt2+rcOHCzo7hkgoXLqzbt2+nOU5JAQAgHSaTydkRXFJ6nyslBQAAGBIlBQAAGBIlBQAAFzRw4ED16NEjS/to3Lixpk6dmk2JHMfdPQAAOMh76fkn+vtZwso4/J7IyEhZrdYcSPPkUFIAAHBBrnDLNKd7AABwQX893dOhQweNGDFCkyZNUqVKlRQYGKhx48YpJSXFNv+PP/5Qr1695O/vr5o1a2r58uWp9pmUlKQhQ4YoMDBQZcuWVfv27XXo0CFJD2/Vbty4sQYNGmSbf/HiRVWqVEnz5s3L1BooKQAAuIG1a9fKw8ND33zzjaZPn64FCxZo/fr1tvGIiAidOnVKX375pVauXKlVq1bpzJkztnGr1aoePXro4sWLWr16tfbs2aMmTZqoU6dOio+PV8GCBbVo0SKtW7dOX375paxWqwYMGKCaNWtq8ODBmcrM6R4AANxA1apV9c4770iSAgMD9dlnn2n37t3q1q2b4uLitGPHDm3btk2NGjWSJC1YsEB16tSxvX/Pnj06cuSI4uLiVKhQIUnSuHHjtG3bNq1evVpDhgxRzZo19d5772nIkCE6ePCgYmNjtW/fvkw/Z8ZtS8qTvujJGApLe1m3e2Ht7il3rT0zF4XCcTVq1LB77e/vrz/++EOSdOzYMeXJk0f16tWzjQcEBKhUqVK21//5z3908+ZNBQYG2u3n9u3bOnXqlO31wIEDtXXrVs2fP1+ffvqpSpcunenMbltSAABwJ/ny5bN7bTKZbHf/ZOQuoJSUFPn6+mrr1q2pxjw9PW3/nZiYqOPHj8vDw0MnT57MUmZKCgAAbq5q1apKSUnRzz//rIYNG0qSzp49q4sXL9rmPP3000pISFCePHlUoUKFNPf1f//3f6pYsaKmTZumN954QyEhIXanjRzBhbMAALi5oKAgtW7dWsOGDbNdSxIREWG79kSSWrRooUaNGql3797asWOHTp8+rYMHD2rKlCnav3+/JGnJkiXat2+fPv74Y4WGhqp379564403dPPmzUzloqQAAADNnz9fAQEB6tSpk3r16qVXXnlFAQEBtnGTyaQ1a9bo+eef15AhQ1S/fn2FhYUpLi5OpUqVktls1rhx4xQVFWU70jJ16lSZTCa9/fbbmcpkslgsuftxdJnknhfOAoDxGP3C2aSkJJd4MJpRPe7z5UgKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJEoKAADpyMgTWeG49D5XSgoAAI9RsGDBTD+MDI938+ZNFSxYMM1xHosPAMBjFChQQPfv31dSUpKzo7icvHnzqkCBAmmPP8EsAADkSkWKFHF2BLfE6R4AAGBIhisp+/btU8+ePVW9enV5e3tr5cqV6b7nl19+Ufv27eXv76/q1asrKiqKi5wAAMjlDFdSkpOTFRwcrMjISLtvX0zLtWvX1LlzZ/n6+mrnzp2KjIzUvHnz9OGHHz6BtAAAIKcY7pqUNm3aqE2bNpKkiIiIdOevXbtWt27d0oIFC1SoUCEFBwfr+PHjmj9/vgYPHiyTyZTTkQEAQA4w3JEURx08eFCNGze2O+rSqlUrXbx4Ub///rsTkwEAgKww3JEURyUkJKh06dJ223x8fGxjFSpUcEIqAEBGmc3mLO8jKCgoG5LAaHJ9SZGU6pTOnxfNcqoHAIyPgoG05PrTPb6+vkpISLDbdvnyZUn/PaICAAByn1xfUho0aKDvvvtOt2/ftm3btWuXSpUqpfLlyzsxGQAAyArDlZQbN24oNjZWsbGxSklJ0blz5xQbG6uzZ89KkiZOnKhOnTrZ5nfr1k2FChVSRESEjh49qk2bNmn27NmKiIjgdA8AALmY4UrKoUOH1KxZMzVr1ky3bt3S1KlT1axZM02ZMkWSFB8fr1OnTtnme3l5acOGDbp48aJatmypt956S4MGDdLgwYOdtQQAAJANTBaLxS0fzeq99LyzIwAAJFnCyjg7AgzKcEdSAAAAJEoKAAAwKEoKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJEoKAAAwJLf9gkF3ZDabFRQU5OwYT5y7rlti7awdyN04kgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAyJkgIAAAzJ4ZJy7do1zZgxQ506dVKTJk30448/SpKuXLmiOXPm6MSJE9keEgAAuJ+8jky+cOGC2rdvr/Pnz6ty5co6fvy4kpOTJUklSpTQsmXLdOHCBUVFReVIWAAA4D4cKikTJkzQtWvXtHv3bvn5+SkwMNBuvEOHDvrmm2+yNSAAAHBPDp3u+fbbbxUeHq7g4GCZTKZU4xUqVNCFCxeyLRwAAHBfDpWUmzdvys/P77HjKSkpWQ4FAADgUEmpXLmyfvrppzTHv/32WwUHB2c5FAAAgEMl5bXXXtOqVau0atUq2xETk8mk5ORkjR8/Xnv27FHfvn1zJCgAAHAvDl04279/f/36668aOHCgPD09JUmvv/66LBaLHjx4oPDwcPXo0SNHggIAAPfiUEmRpA8++EA9e/bUhg0bdPLkSaWkpKhixYrq0qWLGjdunBMZAQCAG3K4pEhSw4YN1bBhw+zOAgAAYMNj8QEAgCE5fCRly5YtWr58uU6fPi2LxSKr1Wo3bjKZ9Ouvv2ZbQAAA4J4cKilRUVGKioqSl5eXatasqUqVKuVULgAA4OYcKimLFi1S8+bNtWrVKhUoUCCnMgEAADh2Tcq9e/fUqVMnCgoAAMhxDpWUkJAQHTp0KKeyAAAA2JgsFos1/WkPXb58WZ07d1aHDh306quvqmzZso/8osHcwHvpeWdHAOAElrAyzo6Q48xms4KCgpwdA8gyh65Jeeqpp9S1a1dNmjRJ06ZNe+Qck8mkxMTEbAkHAADcl0MlZcKECZo7d67Kly+vevXqqVixYjmVCwAAuDmHSsqyZcvUvn17rVixIqfyAAAASHLwwtmUlBS1atUqp7IAAADYOFRS2rVrp7179+ZUFgAAABuHSsqIESNkNps1ZMgQ/fjjj4qPj9cff/yR6hcAAEBWOXRNSv369SVJR44c0fLly9Ocd+XKlaylAgAAbs+hkjJq1Khc+1wUAACQuzhUUsaOHZtTOQAAAOw4dE3K/7p+/bquX7+eXVkAAABsHC4pZ86cUXh4uCpVqqTy5curfPnyqlSpkgYMGKAzZ87kREYAAOCGHDrdYzab9eKLLyopKUktWrRQ1apVZbVaZTabtXbtWu3YsUPbt29XYGBgTuUFAABuwqGSMnHiRFmtVu3atUu1a9e2Gzty5IhCQ0M1ceLEx975AwAAkBEOne7Zu3evwsPDUxUUSapVq5b69eunmJiYbAsHAADcl0Ml5e7du4/9UkEvLy/dvXs3y6Gio6NVu3Zt+fn5qXnz5tq/f/9j5//yyy9q3769/P39Vb16dUVFRclqtWY5BwAAcB6HSkpwcLBWr16tW7dupRq7c+eOVq9ereDg4CwFWr9+vcaMGaMRI0Zoz549atCggV555RWdPXv2kfOvXbumzp07y9fXVzt37lRkZKTmzZunDz/8MEs5AACAc5ksFkuGDzls2bJFf/vb31SlShX17dtXQUFBkqTjx49ryZIlMpvNWrFihdq1a5fpQK1atVKNGjU0d+5c27ZnnnlGoaGheu+991LNX7x4sSZMmKDjx4+rUKFCkqTp06dryZIlOnr0aJoPn/Neej7TGQHkXpawMs6OkOPMZrPtz2cgN3PoSEr79u21cOFCWSwWjRo1Sl26dFGXLl00evRoWSwWLVy4MEsF5e7duzp8+LBCQkLstoeEhOjAgQOPfM/BgwfVuHFjW0GRHhadixcv6vfff890FgAA4FwO3d0jSd26ddPLL7+sw4cP256LEhAQoDp16ihvXod3ZycxMVEPHjyQj4+P3XYfHx8lJCQ88j0JCQkqXbp0qvl/jlWoUCFLmQAAgHNkqlXkzZtXzz77rJ599tnsziNJqU7RWK3Wx35n0KPmP2o7AJjNZmdHeCLcZZ1/4vSWa3KopCxbtkw7duxI8zkoffr0Udu2bdW7d+9MhSlZsqQ8PDxSHTW5fPlyqqMrf/L19X3kfElpvgeA+3KHv8y4JgWuwqFrUpYsWSI/P780x/39/RUdHZ3pMPnz51edOnW0a9cuu+27du1Sw4YNH/meBg0a6LvvvtPt27ft5pcqVUrly5fPdBYAAOBcDpWUEydOqEaNGmmOV69eXXFxcVkKNGjQIH3++edatmyZjh07ptGjRys+Pl5hYWGSHj71tlOnTrb53bp1U6FChRQREaGjR49q06ZNmj17tiIiIjjdAwBALubQ6R6TyaTExMQ0x69cuaKUlJQsBerSpYuuXLmi6dOn69KlS6pevbrWrFmjgIAASVJ8fLxOnTplm+/l5aUNGzZo5MiRatmypby9vTVo0CANHjw4SzkAAIBzOfSclI4dOyohIUG7d+9WwYIF7cZu3bqlli1bqkSJEtqyZUu2B81uPCcFcE88JwXIPRw63TN8+HDbNyFv3LhRZrNZcXFx2rhxo9q1ayez2azhw4fnVFYAAOBGHDrd07JlS82fP1+jRo2yXSMiPbzl19PTU/PmzVPr1q2zPSQAAHA/Dj8npWfPnurQoYN27typ06dPy2q1qmLFigoJCZGnp2dOZAQAAG7IoZKyb98+Va1aVU899ZRCQ0NTjScmJuq3337Tc889l20BAQCAe3LompSOHTumeobJX+3evVsdO3bMcigAAACHSsqfj5tPy927d5Unj0O7BAAAeKR0T/dcu3ZNSUlJttdXrlzR2bNnU82zWCz64osvVKpUqexNCAAA3FK6JWX+/PmaNm2apIcPcxs7dqzGjh37yLlWq1Xjx4/P3oQAAMAtpVtSWrRooYIFC8pqtWrSpEnq0qWLatWqZTfHZDKpcOHCqlu3bo59MzIAAHAv6ZaURo0aqVGjRpKkO3fuqGPHjo/9/h4AAIDs4NAtyGPGjMmpHAAAAHYcKilRUVHpzjGZTBo1alSmAwEAAEgOfsFg8eLF096RySSr1SqTyaQrV65kS7icxBcMAu6JLxgEcg+HjqRcvXo11baUlBSdOXNGn3zyiQ4cOKB169ZlW7ic5A5/UP0vd/2Dy13XLbF2d1074Cqy/OS1PHnyqEKFCpo6darKly/PdSsAACBbZOvjYZ9//nlt3749O3cJAADcVLaWFLPZnO6j8wEAADLC4W9BfpSkpCTFxMRo0aJFevnll7MjFwAAcHMOlZSXXnpJJpMp1Xar1SoPDw917do1Q7cpAwAApMehkvLVV1+l2mYymeTt7a2AgAB5enpmWzAAAODeHCopTZs2zakcAAAAdhwqKZJ0+vRpHTt2TNevX1fRokVVrVo1VahQIQeiAQAAd5bhkrJx40ZFRkbq2LFjqcaqVaum0aNHKzQ0NFvDAQAA95WhkjJ58mTNmjVLnp6e6tGjh2rVqqWiRYvqxo0bOnLkiLZu3aqwsDANHz5c48aNy+nMAADADaRbUv71r39p5syZ6tSpk+bOnSsvL69Uc65du6YhQ4Zo1qxZatKkiUJCQnIkLAAAcB/pPsztk08+UY0aNbR06dJHFhRJKlasmBYvXqzg4GAtWLAg20MCAAD3k25J+emnn9S9e3flyfP4qXny5FH37t31888/Z1s4AADgvtItKcnJySpRokSGdla8eHElJydnORQAAEC6JaVUqVI6evRohnZ29OhR+fv7ZzkUAABAuiWldevWWr58uU6ePPnYeSdPntSKFSv0wgsvZFs4AADgvtItKcOHD5eHh4fatm2rVatW6e7du3bjd+/e1apVq9S+fXt5eHho2LBhORYWAAC4j3RvQS5VqpTWrl2rPn36KCIiQsOGDVNgYKA8PT11/fp1xcXF6c6dO/L19dWaNWtUunTpJ5EbAAC4uAw9zK1+/fr6/vvvtXTpUm3btk2//fabbty4oaJFi6p27dpq27atwsLC5O3tncNxAQCAu8jwY/G9vLw0dOhQDR06NAfjAAAAPJTuNSkAAADOQEkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGREkBAACGZLJYLFZnh3AG76XnnR0BAAzBElbG2RGAR+JICgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRDlZRNmzapS5cuqly5sry9vRUTE5Oh9+3du1fNmzeXn5+fnn76aS1ZsiSHkwIAgJxmqJJy8+ZNNWjQQO+//36G33P69Gl1795dDRo00J49ezR8+HCNGjVKGzduzMGkAAAgp+V1doC/6tmzpyQpMTExw+9ZunSp/P39NX36dElS1apV9eOPP+rDDz9UaGhojuQEAAA5z1BHUjLj4MGDCgkJsdvWqlUrHTp0SPfu3XNSKgAAkFWGOpKSGQkJCWrRooXdNh8fH92/f1+JiYny9/d3TjAAyCXMZrOzI2RZUFCQsyMgBzitpKxZs0bDhg2zvV67dq2aNGmSqX2ZTCa711ar9ZHbAQCp8Rc8jMppJaVdu3Z69tlnba9LlSqVqf34+voqISHBbtvly5eVN29elShRIksZAQCA8zitpHh6esrT0zPL+2nQoIE2b95st23Xrl2qW7eu8uXLl+X9AwAA5zDUhbNXr15VbGysjh49Kkk6deqUYmNjdenSJduc8PBwhYeH216HhYXpwoULGjNmjI4dO6Zly5bp888/1+DBg594fgAAkH0MVVK2bNmiZs2aqWPHjpKkN998U82aNbN7ONu5c+d07tw52+sKFSpozZo12r9/v55//nnNmDFDUVFR3H4MAEAuZ7JYLFZnh3AG76XnnR0BAAzBElbG2RGARzLUkRQAAIA/UVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhUVIAAIAhue23ILsjs9msoKAgZ8d44tx13RJrZ+1A7saRFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEgmi8VidXYIAACA/8WRFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEiUFAAAYEhuU1Kio6NVu3Zt+fn5qXnz5tq/f7+zI2W7WbNmqWXLlipXrpwqV66sHj166OjRo3ZzrFarpk6dqmrVqsnf318dOnTQr7/+6qTEOWPmzJny9vbWW2+9ZdvmyuuOj4/XgAEDVLlyZfn5+alhw4bau3evbdxV1/7gwQNNnjzZ9v917dq1NXnyZN2/f982x1XWvm/fPvXs2VPVq1eXt7e3Vq5caTeekXXeuXNHb731lipVqqTSpUurZ8+eOn/+/JNcBuAwtygp69ev15gxYzRixAjt2bNHDRo00CuvvKKzZ886O1q22rt3r/r27avt27dr06ZNyps3r15++WVdvXrVNmfOnDn66KOPFBUVpZ07d8rHx0edO3fW9evXnZg8+/zwww/67LPPVKNGDbvtrrpui8WiF198UVarVWvWrNGBAwc0bdo0+fj42Oa46tpnz56t6OhoRUVF6eDBg4qMjNSiRYs0a9Ys2xxXWXtycrKCg4MVGRmpQoUKpRrPyDrHjh2rr776SosXL9aWLVt0/fp19ejRQw8ePHiSSwEc4hbPSWnVqpVq1KihuXPn2rY988wzCg0N1XvvvefEZDnrxo0bCggI0MqVK9WuXTtZrVZVq1ZN/fr108iRIyVJt27dUlBQkP7xj38oLCzMyYmzJikpSc2bN9ecOXM0bdo0BQcHa/r06S697kmTJmnfvn3avn37I8ddee09evRQ8eLF9fHHH9u2DRgwQFevXtXq1atddu1lypTRtGnT9Oqrr0rK2M84KSlJgYGB+uijj9S9e3dJ0rlz51SrVi2tW7dOrVq1ctp6gMdx+SMpd+/e1eHDhxUSEmK3PSQkRAcOHHBSqifjxo0bSklJkbe3tyTp999/16VLl+w+i0KFCqlJkyYu8VkMHTpUoaGhat68ud12V1735s2bVa9ePYWFhSkwMFBNmzbVwoULZbU+/LeHK6+9UaNG2rt3r44fPy5J+u233xQTE6MXXnhBkmuv/a8yss7Dhw/r3r17dnPKli2rqlWrutRnAdeT19kBclpiYqIePHhgd/hbknx8fJSQkOCkVE/GmDFjVKtWLTVo0ECSdOnSJUl65Gdx8eLFJ54vO3322Wc6efKkPvnkk1Rjrrzu06dPa/HixYqIiNDQoUN15MgRjR49WpLUv39/l1770KFDdePGDTVs2FAeHh66f/++Ro4cqTfeeEOSa//c/yoj60xISJCHh4dKliyZao6r/zmI3M3lS8qfTCaT3Wur1Zpqmyt5++239f3332vbtm3y8PCwG3O1z8JsNmvSpEnaunWr8ufPn+Y8V1u3JKWkpKhu3bq205ZPP/20Tp48qejoaPXv3982zxXXvn79eq1atUrR0dGqVq2ajhw5ojFjxiggIEB9+vSxzXPFtT9KZtbpqp8FXIfLn+4pWbKkPDw8Uv1r4fLly6n+5eEqxo4dqy+++EKbNm1ShQoVbNv9/PwkyeU+i4MHDyoxMVGNGzdWyZIlVbJkSe3bt0/R0dEqWbKkSpQoIcn11i09/JlWrVrVbluVKlV07tw527jkmmt/9913NXjwYHXt2lU1atRQz549NWjQIH3wwQeSXHvtf5WRdfr6+urBgwdKTExMcw5gRC5fUvLnz686depo165ddtt37dqlhg0bOilVzhk9erTWrVunTZs2qUqVKnZj5cuXl5+fn91ncfv2bX333Xe5+rPo0KGD9u/fr5iYGNuvunXrqmvXroqJiVFgYKBLrlt6eF1GXFyc3ba4uDiVK1dOkuv+zCXp5s2bqY4Senh4KCUlRZJrr/2vMrLOOnXqKF++fHZzzp8/r2PHjrnUZwHX4xanewYNGqTw8HDVq1dPDRs21JIlSxQfH59rr+5Py8iRI7V69WqtWLFC3t7etnPVRYoUUdGiRWUymTRw4EDNnDlTQUFBCgwM1IwZM1SkSBF169bNyekzz9vb23Zx8J8KFy6s4sWLKzg4WJJcct2SFBERoTZt2mjGjBnq0qWLYmNjtXDhQo0fP16SXPZnLklt27bV7NmzVb58eVWrVk2xsbH66KOP1LNnT0mutfYbN27o5MmTkh6e4jt37pxiY2NVvHhxlStXLt11enl56e9//7veffdd+fj4qHjx4nrnnXdUo0YNtWjRwokrAx7PLW5Blh4+zG3OnDm6dOmSqlevrilTpui5555zdqxs9b9/Uf9p9OjRGjt2rKSH56AjIyP16aefymKxqF69epoxY4btL3NX0aFDB9styJJrr3v79u2aNGmS4uLiVLZsWfXr10/h4eG2aw1cde3Xr1/X+++/r6+//lqXL1+Wn5+funbtqlGjRqlgwYKSXGftMTEx6tixY6rtvXr10oIFCzK0ztu3b2v8+PFat26dbt++rWbNmmnmzJkqW7bsk1wK4BC3KSkAACB3cflrUgAAQO5ESQEAAIZESQEAAIZESQEAAIZESQEAAIZESQEAAIZESQEAAIZESQGQroULF2rlypXOjgHAzfAwNwDpql+/vnx9fbV582ZnRwHgRjiSAgAADImSAhhEfHy8hg4dquDgYPn6+qpWrVp68803df36dUnS2bNn1a9fP1WqVEl+fn5q2rSp/vnPf9rtIyYmRt7e3oqJiUm1f29vb02dOtX2eurUqfL29taJEyc0bNgwVaxYUWXKlNFrr72mK1eu2ObVqlVLZrNZ+/bts32ZY61atXLoUwCA/3KLb0EGjO7SpUtq1aqVLl++rD59+ig4OFjx8fH6+uuvdeXKFd29e1dt27bV1atX1b9/f/n7+2v9+vUaOHCgLBaLBg4cmOnfu2/fvvLz89M777yjEydOaOHChcqXL5+io6MlPSwzI0eOVLFixTRixAhJD79ZGwByGiUFMIAJEybowoUL2rx5s5o0aWLbPnbsWFmtVo0fP17nz5/Xxo0b1bx5c0nS66+/rnbt2mny5Mnq3bu3vLy8MvV7V6lSRQsXLrS9tlqtWrRokWbOnCkvLy+99NJLmjhxonx8fNSjR4+sLRQAHMDpHsDJUlJStHnzZrVu3dquoPzJZDJp+/btql27tq2gSFL+/Pk1cOBAJScna+/evZn+/fv27Wv3+rnnntODBw907ty5TO8TALIDJQVwssuXL+vatWsKDg5Oc86ZM2dUpUqVVNurVq1qG8+scuXK2b329vaWJF29ejXT+wSA7EBJAZzMan34FACTyZTp9/4prX08ePAgzX14eHhkaN8A8KRRUgAn8/HxUbFixXT06NE05wQEBOj48eOptpvNZtu49N+jIElJSXbzsnKkRcpcgQKArKKkAE6WJ08edejQQTt27NCBAwdSjVutVr344ouKjY3Vnj17bNvv3bunjz/+WIULF1bTpk0lPSwrHh4eqW5B/uuFsZlRuHBhWSyWLO0DABzF3T2AAbz33nv697//rdDQUL322muqXr26EhIS9NVXX2nFihUaNmyY1q9fr169eik8PFx+fn7asGGDfvjhB02ZMsV2Z0+xYsXUtWtXRUdHy2QyKSgoSDExMTp9+nSW8tWtW1effvqpIiMjFRgYqCJFiqhdu3bZsHIASBslBTAAf39/ffvtt3r//fe1fv16JSUlyd/fXyEhISpZsqSKFi2q7du3a+LEiVq6dKlu3rypwMBALViwQL169bLbV1RUlO7fv68VK1YoT548atOmjdatW6fAwMBM5xszZowuXryo+fPn69q1aypXrhwlBUCO47t7AACAIXFNCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMCRKCgAAMKT/B6n/BXVeVYbAAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "outcome_table = pd.DataFrame({'Outcome':outcomes})\n",
    "\n",
    "outcome_table = outcome_table.reset_index()\n",
    "\n",
    "outcome_table = outcome_table.groupby(['Outcome']).count().reset_index()\n",
    "\n",
    "outcome_table.plot.barh('Outcome')\n",
    "\n",
    "plt.legend(bbox_to_anchor=(1.04,1), loc=\"upper left\")\n",
    "\n",
    "plt.xlabel('count')\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Not surprisingly, each of the three outcomes -1, 0, and 1 appeared about about 100 of the 300 times, give or take. We will examine the \"give or take\" amounts more closely in later chapters."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "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.6.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}