ch3_probability.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Probability"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Hot Hands"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Basketball players who make several baskets in succession are described as having a hot hand. Fans and players have long believed in the hot hand phenomenon, which refutes the assumption that each shot is independent of the next. However, a 1985 [paper](https://www.cs.colorado.edu/~mozer/Teaching/syllabi/7782/readings/gilovich%20vallone%20tversky.pdf) by Gilovich, Vallone, and Tversky collected evidence that contradicted this belief and showed that successive shots are independent events. This paper started a great controversy that continues to this day, as you can see by Googling hot hand basketball.\n",
    "\n",
    "We do not expect to resolve this controversy today. However, in this lab we'll apply one approach to answering questions like this. The goals for this lab are to (1) think about the effects of independent and dependent events, (2) learn how to simulate shooting streaks in Python, and (3) to compare a simulation to actual data in order to determine if the hot hand phenomenon appears to be real."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Getting started"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Our investigation will focus on the performance of one player: Kobe Bryant of the Los Angeles Lakers. His performance against the Orlando Magic in the 2009 NBA finals earned him the title Most Valuable Player and many spectators commented on how he appeared to show a hot hand. Let's load some data from those games and look at the first several rows."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>vs</th>\n",
       "      <th>game</th>\n",
       "      <th>quarter</th>\n",
       "      <th>time</th>\n",
       "      <th>description</th>\n",
       "      <th>basket</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>9:47</td>\n",
       "      <td>Kobe Bryant makes 4-foot two point shot</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>9:07</td>\n",
       "      <td>Kobe Bryant misses jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>8:11</td>\n",
       "      <td>Kobe Bryant misses 7-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>7:41</td>\n",
       "      <td>Kobe Bryant makes 16-foot jumper (Derek Fisher...</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>7:03</td>\n",
       "      <td>Kobe Bryant makes driving layup</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    vs  game quarter  time                                        description  \\\n",
       "0  ORL     1       1  9:47            Kobe Bryant makes 4-foot two point shot   \n",
       "1  ORL     1       1  9:07                          Kobe Bryant misses jumper   \n",
       "2  ORL     1       1  8:11                   Kobe Bryant misses 7-foot jumper   \n",
       "3  ORL     1       1  7:41  Kobe Bryant makes 16-foot jumper (Derek Fisher...   \n",
       "4  ORL     1       1  7:03                    Kobe Bryant makes driving layup   \n",
       "\n",
       "  basket  \n",
       "0      H  \n",
       "1      M  \n",
       "2      M  \n",
       "3      H  \n",
       "4      H  "
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import warnings\n",
    "warnings.filterwarnings(\"ignore\")\n",
    "\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import io\n",
    "import requests\n",
    "\n",
    "df_url = 'https://raw.githubusercontent.com/akmand/datasets/master/openintro/kobe.csv'\n",
    "url_content = requests.get(df_url, verify=False).content\n",
    "kobe = pd.read_csv(io.StringIO(url_content.decode('utf-8')))\n",
    "\n",
    "kobe.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's check the number of games played."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "5\n"
     ]
    }
   ],
   "source": [
    "print(kobe['game'].nunique())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Five games were played. NBA games are played in four quarters. However, if the score is equal at the end of the 4th quarter, overtimes are played to determine the winning team. Let's check if any game went to overtime."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['1' '2' '3' '4' '1OT']\n"
     ]
    }
   ],
   "source": [
    "print(kobe['quarter'].unique())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>vs</th>\n",
       "      <th>game</th>\n",
       "      <th>quarter</th>\n",
       "      <th>time</th>\n",
       "      <th>description</th>\n",
       "      <th>basket</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td>ORL</td>\n",
       "      <td>2</td>\n",
       "      <td>1OT</td>\n",
       "      <td>4:13</td>\n",
       "      <td>Kobe Bryant misses 22-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>54</th>\n",
       "      <td>ORL</td>\n",
       "      <td>2</td>\n",
       "      <td>1OT</td>\n",
       "      <td>2:17</td>\n",
       "      <td>Kobe Bryant makes 11-foot two point shot</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>80</th>\n",
       "      <td>ORL</td>\n",
       "      <td>4</td>\n",
       "      <td>1OT</td>\n",
       "      <td>4:13</td>\n",
       "      <td>Kobe Bryant makes 11-foot jumper</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>81</th>\n",
       "      <td>ORL</td>\n",
       "      <td>4</td>\n",
       "      <td>1OT</td>\n",
       "      <td>3:32</td>\n",
       "      <td>Kobe Bryant makes 19-foot jumper</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82</th>\n",
       "      <td>ORL</td>\n",
       "      <td>4</td>\n",
       "      <td>1OT</td>\n",
       "      <td>2:49</td>\n",
       "      <td>Kobe Bryant misses 10-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>83</th>\n",
       "      <td>ORL</td>\n",
       "      <td>4</td>\n",
       "      <td>1OT</td>\n",
       "      <td>1:58</td>\n",
       "      <td>Kobe Bryant misses 18-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>84</th>\n",
       "      <td>ORL</td>\n",
       "      <td>4</td>\n",
       "      <td>1OT</td>\n",
       "      <td>0:47</td>\n",
       "      <td>Kobe Bryant misses 15-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     vs  game quarter  time                               description basket\n",
       "53  ORL     2     1OT  4:13         Kobe Bryant misses 22-foot jumper      M\n",
       "54  ORL     2     1OT  2:17  Kobe Bryant makes 11-foot two point shot      H\n",
       "80  ORL     4     1OT  4:13          Kobe Bryant makes 11-foot jumper      H\n",
       "81  ORL     4     1OT  3:32          Kobe Bryant makes 19-foot jumper      H\n",
       "82  ORL     4     1OT  2:49         Kobe Bryant misses 10-foot jumper      M\n",
       "83  ORL     4     1OT  1:58         Kobe Bryant misses 18-foot jumper      M\n",
       "84  ORL     4     1OT  0:47         Kobe Bryant misses 15-foot jumper      M"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "kobe[kobe['quarter'] == '1OT']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can see that Game 2 and Game 4 went to overtime (encoded as \"1OT\")."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Before we begin, we need to make sure the data is properly sorted by `game` and `quarter`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>vs</th>\n",
       "      <th>game</th>\n",
       "      <th>quarter</th>\n",
       "      <th>time</th>\n",
       "      <th>description</th>\n",
       "      <th>basket</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>9:47</td>\n",
       "      <td>Kobe Bryant makes 4-foot two point shot</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>9:07</td>\n",
       "      <td>Kobe Bryant misses jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>8:11</td>\n",
       "      <td>Kobe Bryant misses 7-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>7:41</td>\n",
       "      <td>Kobe Bryant makes 16-foot jumper (Derek Fisher...</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>7:03</td>\n",
       "      <td>Kobe Bryant makes driving layup</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>128</th>\n",
       "      <td>ORL</td>\n",
       "      <td>3</td>\n",
       "      <td>4</td>\n",
       "      <td>3:33</td>\n",
       "      <td>Bryant Layup Shot: Missed</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>129</th>\n",
       "      <td>ORL</td>\n",
       "      <td>3</td>\n",
       "      <td>4</td>\n",
       "      <td>2:02</td>\n",
       "      <td>Bryant 3pt Shot: Missed</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>130</th>\n",
       "      <td>ORL</td>\n",
       "      <td>3</td>\n",
       "      <td>4</td>\n",
       "      <td>00:23.9</td>\n",
       "      <td>Bryant 3pt Shot: Missed</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>131</th>\n",
       "      <td>ORL</td>\n",
       "      <td>3</td>\n",
       "      <td>4</td>\n",
       "      <td>00:06.9</td>\n",
       "      <td>Bryant 3pt Shot: Missed</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>132</th>\n",
       "      <td>ORL</td>\n",
       "      <td>3</td>\n",
       "      <td>4</td>\n",
       "      <td>00:00.5</td>\n",
       "      <td>Bryant Layup Shot: Made (31 PTS)</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>133 rows × 6 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      vs  game quarter     time  \\\n",
       "0    ORL     1       1     9:47   \n",
       "1    ORL     1       1     9:07   \n",
       "2    ORL     1       1     8:11   \n",
       "3    ORL     1       1     7:41   \n",
       "4    ORL     1       1     7:03   \n",
       "..   ...   ...     ...      ...   \n",
       "128  ORL     3       4     3:33   \n",
       "129  ORL     3       4     2:02   \n",
       "130  ORL     3       4  00:23.9   \n",
       "131  ORL     3       4  00:06.9   \n",
       "132  ORL     3       4  00:00.5   \n",
       "\n",
       "                                           description basket  \n",
       "0              Kobe Bryant makes 4-foot two point shot      H  \n",
       "1                            Kobe Bryant misses jumper      M  \n",
       "2                     Kobe Bryant misses 7-foot jumper      M  \n",
       "3    Kobe Bryant makes 16-foot jumper (Derek Fisher...      H  \n",
       "4                      Kobe Bryant makes driving layup      H  \n",
       "..                                                 ...    ...  \n",
       "128                        Bryant Layup Shot: Missed        M  \n",
       "129                          Bryant 3pt Shot: Missed        M  \n",
       "130                          Bryant 3pt Shot: Missed        M  \n",
       "131                          Bryant 3pt Shot: Missed        M  \n",
       "132                 Bryant Layup Shot: Made (31 PTS)        H  \n",
       "\n",
       "[133 rows x 6 columns]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "kobe"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We know that five games were played in total. Yet the last rows shows the records of Game 3. Therefore, we should sort the data by `game` and `quarter` in ascending orders.\n",
    "\n",
    "First, let's replace \"1OT\" with \"5\" for the `quarter` column to make sorting easier."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "kobe['quarter'] = kobe['quarter'].replace('1OT', '5') \n",
    "kobe['quarter'] = kobe['quarter'].astype(int)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Then, we can sort the data by `game` and `quarter` in ascending orders, and finally begin talking about the concept of probability."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>vs</th>\n",
       "      <th>game</th>\n",
       "      <th>quarter</th>\n",
       "      <th>time</th>\n",
       "      <th>description</th>\n",
       "      <th>basket</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>9:47</td>\n",
       "      <td>Kobe Bryant makes 4-foot two point shot</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>9:07</td>\n",
       "      <td>Kobe Bryant misses jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>8:11</td>\n",
       "      <td>Kobe Bryant misses 7-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>7:41</td>\n",
       "      <td>Kobe Bryant makes 16-foot jumper (Derek Fisher...</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>ORL</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>7:03</td>\n",
       "      <td>Kobe Bryant makes driving layup</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>128</th>\n",
       "      <td>ORL</td>\n",
       "      <td>5</td>\n",
       "      <td>4</td>\n",
       "      <td>6:22</td>\n",
       "      <td>Kobe Bryant misses 18-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>129</th>\n",
       "      <td>ORL</td>\n",
       "      <td>5</td>\n",
       "      <td>4</td>\n",
       "      <td>4:26</td>\n",
       "      <td>Kobe Bryant misses 27-foot three point jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>130</th>\n",
       "      <td>ORL</td>\n",
       "      <td>5</td>\n",
       "      <td>4</td>\n",
       "      <td>3:12</td>\n",
       "      <td>Kobe Bryant misses 27-foot three point jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>131</th>\n",
       "      <td>ORL</td>\n",
       "      <td>5</td>\n",
       "      <td>4</td>\n",
       "      <td>2:38</td>\n",
       "      <td>Kobe Bryant makes 9-foot two point shot</td>\n",
       "      <td>H</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>132</th>\n",
       "      <td>ORL</td>\n",
       "      <td>5</td>\n",
       "      <td>4</td>\n",
       "      <td>2:06</td>\n",
       "      <td>Kobe Bryant misses 13-foot jumper</td>\n",
       "      <td>M</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>133 rows × 6 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      vs  game  quarter  time  \\\n",
       "0    ORL     1        1  9:47   \n",
       "1    ORL     1        1  9:07   \n",
       "2    ORL     1        1  8:11   \n",
       "3    ORL     1        1  7:41   \n",
       "4    ORL     1        1  7:03   \n",
       "..   ...   ...      ...   ...   \n",
       "128  ORL     5        4  6:22   \n",
       "129  ORL     5        4  4:26   \n",
       "130  ORL     5        4  3:12   \n",
       "131  ORL     5        4  2:38   \n",
       "132  ORL     5        4  2:06   \n",
       "\n",
       "                                           description basket  \n",
       "0              Kobe Bryant makes 4-foot two point shot      H  \n",
       "1                            Kobe Bryant misses jumper      M  \n",
       "2                     Kobe Bryant misses 7-foot jumper      M  \n",
       "3    Kobe Bryant makes 16-foot jumper (Derek Fisher...      H  \n",
       "4                      Kobe Bryant makes driving layup      H  \n",
       "..                                                 ...    ...  \n",
       "128                  Kobe Bryant misses 18-foot jumper      M  \n",
       "129      Kobe Bryant misses 27-foot three point jumper      M  \n",
       "130      Kobe Bryant misses 27-foot three point jumper      M  \n",
       "131            Kobe Bryant makes 9-foot two point shot      H  \n",
       "132                  Kobe Bryant misses 13-foot jumper      M  \n",
       "\n",
       "[133 rows x 6 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "kobe.sort_values(by = ['game', 'quarter'], ascending = [True, True], ignore_index = True, inplace = True)  \n",
    "# Note: time is already in descending order for each game and quarter, so we do not need to sort that one.\n",
    "kobe"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this data frame, every row records a shot taken by Kobe Bryant. If he hit the shot (made a basket), a hit, `H`, is recorded in the column named basket, otherwise a miss, `M`, is recorded.\n",
    "\n",
    "Just looking at the string of hits and misses, it can be difficult to gauge whether or not it seems like Kobe was shooting with a hot hand. One way we can approach this is by considering the belief that hot hand shooters tend to go on shooting streaks. For this lab, we define the length of a shooting streak to be the ***number of consecutive baskets made until a miss occurs***.\n",
    "\n",
    "For example, in Game 1 Kobe had the following sequence of hits and misses from his nine shot attempts in the first quarter:\n",
    "\n",
    "##### `H M | M | H H M | M | M | M`\n",
    "\n",
    "To verify this use the following command:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0    H\n",
       "1    M\n",
       "2    M\n",
       "3    H\n",
       "4    H\n",
       "5    M\n",
       "6    M\n",
       "7    M\n",
       "8    M\n",
       "Name: basket, dtype: object"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "kobe['basket'][0:9]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Within the nine shot attempts, there are six streaks, which are separated by a \"|\" above. Their lengths are one, zero, two, zero, zero, zero (in order of occurrence)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class = 'exercise'>\n",
    "<h4>Exercise 1</h4>\n",
    "What does a streak length of 1 mean, i.e. how many hits and misses are in a streak of 1? What about a streak length of 0?\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The custom function `calc_streak()` can be used to calculate the lengths of all shooting streaks."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[1, 0, 2, 0, 0, 0, 3, 2, 0, 3, 0, 1, 3, 0, 0, 0, 0, 0, 1, 1, 0, 4, 1, 0, 1, 0, 1, 0, 1, 4, 3, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 2, 0, 1, 2, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 2, 0, 0, 0, 3, 0, 1, 0, 1, 2, 1, 0, 1, 0, 0, 1]\n"
     ]
    }
   ],
   "source": [
    "def calc_streak(x):\n",
    "    all_shoots = list(x)\n",
    "    streak, count = [], 0\n",
    "    \n",
    "    for shoot in all_shoots:\n",
    "        if shoot == 'H':\n",
    "            count += 1\n",
    "        elif shoot == 'M':\n",
    "            streak.append(count)\n",
    "            count = 0\n",
    "    \n",
    "    # If the last shoot is a hit, make sure to append the last count to the streak list as well.\n",
    "    if all_shoots[-1] == 'H':  \n",
    "        streak.append(count)\n",
    "        \n",
    "    return streak  \n",
    "    \n",
    "kobe_streak = calc_streak(kobe['basket'])\n",
    "print(kobe_streak)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can look at the distribution of the lengths of all shooting streaks  with a barplot."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABmwAAAN7CAYAAABI6LUwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/TGe4hAAAACXBIWXMAAB7CAAAewgFu0HU+AAB1nklEQVR4nOzdeZSW5Z3n/08tAopCQbEJxijQiEhMVMQ1ikTtjFvsxKZts3TsaZ0QkpN0Tk5PMmlbTTSZRHNGJ9HqmE6b7ixjjJm4J+koQXtUXInTCkIoMRGQRQoFEShq+f3hjxoIO9ZTV1Xxev11U89138+3zvG5jvrmvp+q9vb29gAAAAAAAFBMdekBAAAAAAAA9nWCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGG1pQfoadrb29PW1lZ6DHqxmpqaJElra2vhSQA6h30N6G3sa0BvY18Dehv7Gl2luro6VVVVnXY9wWYPtbW1Zfny5aXHoJeqrq7OiBEjkiQrV64UB4Eez74G9Db2NaC3sa8BvY19ja40fPjwjkDYGTwSDQAAAAAAoDDBBgAAAAAAoDDBBgAAAAAAoDDBBgAAAAAAoDDBBgAAAAAAoLDaUm/8wx/+MHfffXfHn6+88socddRROz1nzpw5eeCBB9LY2Jg1a9ZkwIABGTNmTM4888wcc8wxlR4ZAAAAAACgIooEm5deein33Xffbq9va2vLLbfckpkzZ27186ampjQ1NeXJJ5/M1KlTc/nll6e62k1DAAAAAABAz9LlwaatrS3f+c530tramoEDB+b111/f5Tm33XZbR6w5/PDDc8EFF2T48OFZvnx57r777ixatCgzZ87MgAEDcskll1T6VwAAAAAAAOhUXX47yi9+8Ys0NjZm1KhROeOMM3a5funSpbnnnnuSJGPGjMlXvvKVnHLKKRk7dmxOOeWUfPnLX86YMWOSJPfcc0+WLVtW0fkBAAAAAAA6W5cGm1dffTU/+clPkiSXXXZZamt3fYPP/fffn9bW1iTJpZdemj59+mz1et++fXPppZcmSVpbW3Pvvfd28tQAAAAAAACV1aXB5p/+6Z+yYcOGnH766ZkwYcIu17e3t+fJJ59MkowaNSrjxo3b7rpx48Zl5MiRSZKnnnoq7e3tnTc0AAAAAABAhXVZsHn00UfzzDPP5MADD8xHP/rR3TpnxYoVWb16dZLkyCOP3OnazQGoqakpK1eufHvDAgAAAAAAdKEuCTbr1q3L97///STJhz/84QwYMGC3zlu8eHHH8ahRo3a6dvMdNn98HgAAAAAAQHe36y+R6QQ//OEP89prr+WII47I1KlTd/u8VatWdRzX19fvdO2QIUO2e96evs+O1NXVpaamJklSXd2lT5JjH7LlP1v+OQN6A/sa0NvY14Dexr4G9Db2NXqyigebefPmZebMmampqclll12Wqqqq3T53w4YNHcf9+vXb6dq+fftu97zdMX369F2uaWhoSH19fWpqajJixIg9uj7sjWHDhpUeAaBT2deA3sa+BvQ29jWgt7Gv0dNUNDG2tLTklltuSXt7e84999wceuihe3R+c3Nzx3Ft7c7b0n777bfd8wAAAAAAALq7it5h87//9//OkiVLMmTIkFx00UV7fH6fPn06jltaWna6dtOmTds9b3c0NDTsck1dXV2SpLW1NStXrtyj68Puqq6u7ij/K1asSFtbW+GJAN4e+xrQ29jXgN7Gvgb0NvY1utLQoUM7vkqlM1Qs2CxZsiR33nlnkuSv//qvd/lIs+3Z8pxdPeZs48aN2z1vd+zq+3H+mA85XaGtrc0/a0CvYl8Dehv7GtDb2NeA3sa+Rk9TsWBz3333paWlJcOHD8/GjRvzyCOPbLPm5Zdf7jh+7rnn8tprryVJjjvuuPTr12+rkLJq1aqdvt+rr77acbynAQYAAAAAAKCkigWbzY8oW758eW688cZdrv/Zz37Wcfztb387/fr1yyGHHNLxsyVLluz0/KVLl3Ycb3keAAAAAABAd1ddeoCdGTZsWAYNGpQkmTdv3k7Xbn598ODBGTp0aMVnAwAAAAAA6CwVu8NmxowZmTFjxk7X3H777bnjjjuSJFdeeWWOOuqorV6vqqrK8ccfn3/7t3/LkiVLsmDBgowbN26b6yxYsKDjDpxJkyalqqqqk34LAAAAAACAyuvWd9gkyTnnnJPq6rfGvPXWW9Pc3LzV683Nzbn11luTJDU1NTn33HO7fEYAAAAAAIC3o9sHm5EjR+aCCy5IkjQ2NuaKK67Io48+msbGxjz66KO54oor0tjYmCQ5//zzc/DBB5ccFwAAAAAAYI9V7JFoneniiy/O66+/nt/85jdZtGhRbrjhhm3WTJ06NRdffHHXDwcAAAAAAPA29YhgU11dnenTp+eEE07IAw88kMbGxqxduzYHHXRQxowZk7POOivHHHNM6TEBAAAAAAD2StFgM23atEybNm231x977LE59thjKzgRAAAAAABA1+v232EDAAAAAADQ2wk2AAAAAAAAhQk2AAAAAAAAhQk2AAAAAAAAhQk2AAAAAAAAhQk2AAAAAAAAhdWWHoDebf27/6L0CD1OY+kBeqj9n/1J6REAAAAAAPaaO2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKE2wAAAAAAAAKq63Uhd98883MmTMnjY2NaWxsTFNTU9asWZPm5ub0798/hxxySI455phMnTo1Bx100Hav8fzzz+fqq6/erfe76KKLMm3atM78FQAAAAAAALpExYLNwoULc+ONN273tTVr1mTu3LmZO3du7r777nz605/Oe97znkqNAgAAAAAA0K1VLNgkSX19fY466qiMHj06Q4YMSV1dXdrb27Nq1arMnj07TzzxRNauXZtvfOMb+epXv5rDDjtsh9eaPn16xowZs8PXBw4cWIHfAAAAAAAAoPIqFmwmTpyYhoaGHb5+8skn54knnsj111+flpaW3HHHHfn85z+/w/XDhg3LoYceWolRAQAAAAAAiqqu2IWrd33pyZMnZ+TIkUmSefPmVWoUAAAAAACAbq1iwWZ37b///kmSTZs2FZ4EAAAAAACgjKLBZunSpXnppZeSJKNGjSo5CgAAAAAAQDEV+w6bHdm4cWOampry9NNP56677kpra2uS5JxzztnpebfddltWrVqV1157LX379s3QoUMzYcKEnH322R2PVQMAAAAAAOiJuiTYzJo1KzfffPMOX7/wwgtz6qmn7vQa8+fP7zhuaWnJunXr8tJLL+UXv/hFPvShD+XP//zPU1VVtVfzrVq1apdr6urqUlNTk2T3vp8H6Fo+l9A9bfnZ9DkFegP7GtDb2NeA3sa+Rk/W5XfYbOmwww7L5ZdfnrFjx+5wzaBBgzJ58uSMHz8+w4cPT3V1dV599dU888wzeeihh9La2po77rgjLS0tueSSS/ZqjunTp+9yTUNDQ+rr61NTU5MRI0bs1fvsixpLD8A+w+cSur9hw4aVHgGgU9nXgN7Gvgb0NvY1epqq9vb29kq/ybp16zruYmlubs7y5cvz2GOP5Yknnsjw4cPz8Y9/PMcdd9w2523YsCG1tbWprd1+V1q4cGGuueaavPnmm6mqqsrXv/71HHbYYXs837Rp03a5ZnOwYc80Dn1v6RHYR4xZ+e+lRwAAAAAA2GtdEmx25OGHH85NN92U5K27XKZMmbJX1/j2t7+dJJk6dWo+8YlP7PE19uSRaK2trVm5cuUev8e+at27/rz0COwj+v/HT0uPAGxHdXV1x99oWrFiRdra2gpPBPD22NeA3sa+BvQ29jW60tChQzu+SqUzFH0k2mmnnZann346jz32WL73ve9l0qRJOfDAA/foGqecckq+973vZf369Zk3b95ezbGnd874kEP343MJ3V9bW5vPKtCr2NeA3sa+BvQ29jV6muLfunT88ccnSTZu3Jjf/va3e3x+TU1NRo4cmSRpamrqzNEAAAAAAAC6RPFgM2DAgI5jjxoDAAAAAAD2RcWDzZZ3xfTr12+Pz29tbc0rr7ySJBk0aFCnzQUAAAAAANBVigebxx57rOP40EMP3ePzH3300bz55ptJkgkTJnTaXAAAAAAAAF2lYsFm1qxZaW5u3umae++9N3PmzEmSDBs2LEceeWTHa2+88Uaef/75nZ6/cOHC/PM//3OSpKqqKmefffbbnBoAAAAAAKDr1Vbqwj/96U/zr//6rznhhBMyfvz4DB8+PP369cuGDRvyhz/8If/+7/+e+fPnvzVEbW0uv/zyVFf/v3705ptv5uqrr8473/nOHH/88Rk9enTq6upSXV2dV199Nc8880wefvjhtLS0JEnOP//8jB49ulK/DgAAAAAAQMVULNgkb90l8+CDD+bBBx/c4Zr6+vpMnz49Rx999HZf//3vf5/f//73Ozy/uro6H/rQh3LRRRe97XkBAAAAAABKqFiw+dKXvpRnnnkmL7zwQpYvX57XXnstb7zxRvr06ZMBAwbksMMOy3HHHZeTTjopffv23eb8wYMH53Of+1wWLFiQhQsXpqmpKWvXrk1zc3MOOOCAjBw5MkcddVSmTp2aYcOGVerXAAAAAAAAqLiKBZuRI0dm5MiROe+88/bq/Nra2px44ok58cQTO3kyAAAAAACA7qV610sAAAAAAACoJMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgMMEGAAAAAACgsNpKXfjNN9/MnDlz0tjYmMbGxjQ1NWXNmjVpbm5O//79c8ghh+SYY47J1KlTc9BBB+3yevPnz8+vfvWrvPDCC3n99ddzwAEH5LDDDsvpp5+eU089tVK/BgAAAAAAQMVVLNgsXLgwN95443ZfW7NmTebOnZu5c+fm7rvvzqc//em85z3v2eG1br/99vzsZz9Le3t7x89ef/31PPvss3n22Wfzf/7P/8nnPve59OnTp7N/DQAAAAAAgIqrWLBJkvr6+hx11FEZPXp0hgwZkrq6urS3t2fVqlWZPXt2nnjiiaxduzbf+MY38tWvfjWHHXbYNtf49a9/nTvuuCNJMnz48PzZn/1ZDj300KxevTr3339/nn/++TzzzDNpaGjIZz7zmUr+OgAAAAAAABVRsWAzceLENDQ07PD1k08+OU888USuv/76tLS05I477sjnP//5rda88cYb+dGPfpQkGTJkSK699toMGDCg4/Xjjjsu1113XZ5++uk88sgjOfPMM3PUUUdV5hcCAAAAAACokOqKXbh615eePHlyRo4cmSSZN2/eNq8/+OCDefPNN5MkH/7wh7eKNZvf42/+5m863uvuu+9+u2MDAAAAAAB0uYoFm921//77J0k2bdq0zWtPPvlkx5oTTjhhu+fX19fnXe96V5Lkueeey/r16ys0KQAAAAAAQGUUDTZLly7NSy+9lCQZNWrUVq+1tLRk4cKFSZJx48altnbHT2+bMGFCkreiT2NjY2WGBQAAAAAAqJCKfYfNjmzcuDFNTU15+umnc9ddd6W1tTVJcs4552y1bunSpWlra0uybcz5Y1u+vmTJkkycOLGTpwYAAAAAAKicLgk2s2bNys0337zD1y+88MKceuqpW/2sqamp43jw4ME7vX59fX3H8apVq/Z4vt05p66uLjU1NUl27/t5gK7lcwnd05afTZ9ToDewrwG9jX0N6G3sa/RkXX6HzZYOO+ywXH755Rk7duw2r235XTT9+vXb6XX69u3bcbxhw4Y9nmP69Om7XNPQ0JD6+vrU1NRkxIgRe/we+yoPqKOr+FxC9zds2LDSIwB0Kvsa0NvY14Dexr5GT9Mlweb444/P9ddfnyRpbm7O8uXL89hjj+WJJ57IjTfemI9//OM57rjjtjpn06ZN/2/InXx/TZLst99+HcfNzc2dODkAAAAAAEDldUmw6d+/f/r379/x57Fjx+aUU07Jww8/nJtuuinf+MY3Mn369EyZMqVjzZYRpqWlZafX3zLu9OnTZ4/na2ho2OWaurq6JElra2tWrly5x+8BVNayZctKjwBsR3V1dcffaFqxYkXH99MB9FT2NaC3sa8BvY19ja40dOjQjq9S6QxFH4l22mmn5emnn85jjz2W733ve5k0aVIOPPDAJMn+++/fsW5XjznbuHFjx/GuHp+2PVt+B87u8CGH7sfnErq/trY2n1WgV7GvAb2NfQ3obexr9DTFv3Xp+OOPT/JWdPntb3/b8fPBgwd3HDc1Ne30GqtWreo43tP4AgAAAAAAUFrxYDNgwICO4y0fNTZy5MhUV7813pIlS3Z6jS1fHzVqVCdPCAAAAAAAUFnFg82Wd89s+Tiz2trajB07NkmyYMGCnX6Pzdy5c5O89b03Y8aMqdCkAAAAAAAAlVE82Dz22GMdx4ceeuhWr21+XNr69evz+OOPb/f8VatW5T/+4z+SJBMnTtzqu28AAAAAAAB6gooFm1mzZqW5uXmna+69997MmTMnSTJs2LAceeSRW73+vve9LwcccECS5Mc//nHWrl271ettbW35p3/6p44vjrrgggs6a3wAAAAAAIAuU1upC//0pz/Nv/7rv+aEE07I+PHjM3z48PTr1y8bNmzIH/7wh/z7v/975s+f/9YQtbW5/PLLO76zZrMDDzwwH/7wh/Pd7343K1euzH/7b/8tH/zgB3PooYdm9erVue+++/L8888nSU455ZQcddRRlfp1AAAAAAAAKqZiwSZJ3njjjTz44IN58MEHd7imvr4+06dPz9FHH73d188666ysXr06P/vZz7J8+fI0NDRss+aYY47J9OnTO21uAAAAAACArlSxYPOlL30pzzzzTF544YUsX748r732Wt5444306dMnAwYMyGGHHZbjjjsuJ510Uvr27bvTa02bNi3vfve786tf/Srz5s3L66+/nv79++ed73xnpkyZklNPPbVSvwYAAAAAAEDFVSzYjBw5MiNHjsx5553XKdc74ogjcsQRR3TKtQAAAAAAALqT6l0vAQAAAAAAoJIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMIEGwAAAAAAgMJqK3nxxsbGzJkzJy+88EIWL16cNWvWpKamJoMHD84RRxyRqVOnZvz48Tu9xqxZs3LzzTfv1vt98pOfzJQpUzphcgAAAAAAgK5TsWBz5ZVXZt68edv8vKWlJa+88kpeeeWVzJo1K6eddlo+8YlPpLa2ou0IAAAAAACg26pYJWlqakqSDBo0KCeddFLGjx+fIUOGpK2tLQsWLMi9996bpqamPPzww2ltbc1nPvOZXV7zS1/6UgYNGrTD1+vr6zttfgAAAAAAgK5SsWAzatSo/OVf/mVOPPHEVFdv/VU548aNy2mnnZYrrrgir7zySh555JGcddZZmTBhwk6vefDBB2fYsGGVGhkAAAAAAKCI6l0v2Ttf+MIXcvLJJ28TazYbMGBAPvaxj3X8efbs2ZUaBQAAAAAAoFurWLDZHUcddVTH8fLlywtOAgAAAAAAUE7RYNPS0tJxvKM7cQAAAAAAAHq7in2Hze6YO3dux/GoUaN2ub6hoSFLly7NmjVrcsABB2TEiBF517velbPPPjuDBw+u5KgAAAAAAAAVUyzYtLW15c477+z488knn7zLc55//vmO47Vr12bt2rX53e9+l3vuuScf//jHc9ZZZ+3VLKtWrdrlmrq6utTU1CRxNxB0Rz6X0D1t+dn0OQV6A/sa0NvY14Dexr5GT1Ys2Nx3331ZuHBhkmTy5MkZPXr0DtcOHz48kydPzrhx41JfX58kWbFiRWbPnp3HH388mzZtyne/+91UVVXlzDPP3ONZpk+fvss1DQ0Nqa+vT01NTUaMGLHH77Gvaiw9APsMn0vo/oYNG1Z6BIBOZV8Dehv7GtDb2NfoaYoEm7lz5+bHP/5xkmTgwIG57LLLdrh28uTJOf3001NVVbXVz8eOHZuTTz45Tz/9dK6//vq0trbmX/7lXzJp0qTU1dVVcnwAAAAAAIBOVdXe3t7elW/48ssv5x/+4R+ybt267LfffvnSl76UCRMmvK1r/uxnP8tPfvKTJMnFF1+cD37wg3t0/p48Eq21tTUrV67cqzn3Reve9eelR2Af0f8/flp6BGA7qqurO/5G04oVK9LW1lZ4IoC3x74G9Db2NaC3sa/RlYYOHdrxVSqdoUvvsFmxYkWuueaarFu3LtXV1fnsZz/7tmNNkpx55pm5/fbb097enrlz5+5xsNn8mLXd5UMO3Y/PJXR/bW1tPqtAr2JfA3ob+xrQ29jX6Gm67FuXmpqa8pWvfCWrV69OVVVVpk+fnuOPP75Trj1w4MAceOCBHe8DAAAAAADQk3RJsFmzZk2uueaaLF++PEly6aWX5vTTT+/U9/jj77gBAAAAAADoKSoebN58881ce+21Wbx4cZLkkksuyfvf//5OfY81a9Zk7dq1SZJBgwZ16rUBAAAAAAAqraLBZuPGjfna176WRYsWJUk++MEP5sILL+z093nggQfS3t6eJJ3ynTgAAAAAAABdqWLBpqWlJddff33mz5+fJDnnnHNy8cUX79E1VqxY0RF7duTpp5/OHXfckSTp06dPzjjjjL0bGAAAAAAAoJDaSl34hhtuyLPPPpskmThxYqZOnZo//OEPOx6ktjYjR47c6mcrV67M1VdfnXHjxuW4447LO9/5zgwcODBJsnz58syePTuPP/54x901H/3oRzN48OAK/UYAAAAAAACVUbFg88QTT3QcP/fcc/n85z+/0/VDhw7NTTfdtN3XFixYkAULFuzw3L59++av/uqvcuaZZ+7dsAAAAAAAAAVVLNh0htGjR+fTn/50FixYkBdffDGrV6/O2rVr09ramv79++cd73hHJk6cmPe9730dd94AAAAAAAD0NBULNrfffvvbvsb++++f9773vXnve9/bCRMBAAAAAAB0T9WlBwAAAAAAANjXCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACFCTYAAAAAAACF1ZYeAAB6kvXv/ovSI/Q4jaUH6KH2f/YnpUcAAAAAupA7bAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAAqrreTFGxsbM2fOnLzwwgtZvHhx1qxZk5qamgwePDhHHHFEpk6dmvHjx+/29ebMmZMHHnggjY2NWbNmTQYMGJAxY8bkzDPPzDHHHFPB3wQAAAAAAKByKhZsrrzyysybN2+bn7e0tOSVV17JK6+8klmzZuW0007LJz7xidTW7niUtra23HLLLZk5c+ZWP29qakpTU1OefPLJTJ06NZdffnmqq900BAAAAAAA9CwVCzZNTU1JkkGDBuWkk07K+PHjM2TIkLS1tWXBggW5995709TUlIcffjitra35zGc+s8Nr3XbbbR2x5vDDD88FF1yQ4cOHZ/ny5bn77ruzaNGizJw5MwMGDMgll1xSqV8JAAAAAACgIioWbEaNGpW//Mu/zIknnrjNXS/jxo3LaaedliuuuCKvvPJKHnnkkZx11lmZMGHCNtdZunRp7rnnniTJmDFjcvXVV6dPnz5JkrFjx2bSpEm56qqr0tjYmHvuuSdTp07NiBEjKvVrAQAAAAAAdLqKPT/sC1/4Qk4++eQdPqJswIAB+djHPtbx59mzZ2933f3335/W1tYkyaWXXtoRazbr27dvLr300iRJa2tr7r333s4YHwAAAAAAoMsU/cKXo446quN4+fLl27ze3t6eJ598Mslbd+yMGzduu9cZN25cRo4cmSR56qmn0t7eXoFpAQAAAAAAKqNosGlpaek43t6dOCtWrMjq1auTJEceeeROr7X5cWpNTU1ZuXJlJ04JAAAAAABQWUWDzdy5czuOR40atc3rixcv3unrW9p8h80fnwcAAAAAANDd1ZZ647a2ttx5550dfz755JO3WbNq1aqO4/r6+p1eb8iQIds9b3fszvq6urrU1NQk2f7dQEBZPpdAb2Nfg+5py8+mzynQG9jXgN7GvkZPVizY3HfffVm4cGGSZPLkyRk9evQ2azZs2NBx3K9fv51er2/fvts9b3dMnz59l2saGhpSX1+fmpqajBgxYo+uvy9rLD0A+wyfS7qKfY2uYl+D7m/YsGGlRwDoVPY1oLexr9HTFEmMc+fOzY9//OMkycCBA3PZZZdtd11zc3PHcW3tztvSfvvtt93zAAAAAAAAursuv8Pm5ZdfznXXXZfW1tbst99++du//dsMHDhwu2v79OnTcdzS0rLT627atGm75+2OhoaGXa6pq6tLkrS2tmblypV7dH2g8pYtW1Z6BIBOZV+D7qm6urrjb2quWLEibW1thScCeHvsa0BvY1+jKw0dOrTjq1Q6Q5cGmxUrVuSaa67JunXrUl1dnc9+9rOZMGHCDtdv+Ri0XT3mbOPGjds9b3fs6vtx/pgPOXQ/PpdAb2Nfg+6vra3NZxXoVexrQG9jX6On6bJHojU1NeUrX/lKVq9enaqqqkyfPj3HH3/8Ts/ZMqSsWrVqp2tfffXV7Z4HAAAAAADQ3XVJsFmzZk2uueaaLF++PEly6aWX5vTTT9/leYccckjH8ZIlS3a6dunSpds9DwAAAAAAoLureLB58803c+2112bx4sVJkksuuSTvf//7d+vcYcOGZdCgQUmSefPm7XTt5tcHDx6coUOHvo2JAQAAAAAAulZFg83GjRvzta99LYsWLUqSfPCDH8yFF1642+dXVVV1PDZtyZIlWbBgwXbXLViwoOMOnEmTJqWqqurtDQ4AAAAAANCFKhZsWlpacv3112f+/PlJknPOOScXX3zxHl/nnHPOSXX1W2PeeuutaW5u3ur15ubm3HrrrUmSmpqanHvuuW9zcgAAAAAAgK5VW6kL33DDDXn22WeTJBMnTszUqVPzhz/8YceD1NZm5MiR2/x85MiRueCCC3LnnXemsbExV1xxRT7wgQ9k+PDhWb58ee66666OO3jOP//8HHzwwZX5hQAAAAAAACqkYsHmiSee6Dh+7rnn8vnPf36n64cOHZqbbrppu69dfPHFef311/Ob3/wmixYtyg033LDNmqlTp+7VHTwAAAAAAAClVSzYdKbq6upMnz49J5xwQh544IE0NjZm7dq1OeiggzJmzJicddZZOeaYY0qPCQAAAAAAsFcqFmxuv/32Tr/msccem2OPPbbTrwsAAAAAAFBSdekBAAAAAAAA9nWCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGG1lbz466+/noULF2bhwoVpbGxMY2Nj1q5dmyQ5/fTTM2PGjF1eY9asWbn55pt36/0++clPZsqUKW9nZAAAAAAAgC5X0WBz2WWXVfLyAAAAAAAAvUJFg82WhgwZklGjRuXZZ5/d62t86UtfyqBBg3b4en19/V5fGwAAAAAAoJSKBpuLLrooY8aMyZgxY1JXV5cVK1bkU5/61F5f7+CDD86wYcM6cUIAAAAAAIDyKhpspk2bVsnLAwAAAAAA9ArVpQcAAAAAAADY1wk2AAAAAAAAhVX0kWidraGhIUuXLs2aNWtywAEHZMSIEXnXu96Vs88+O4MHDy49HgAAAAAAwF7pUcHm+eef7zheu3Zt1q5dm9/97ne555578vGPfzxnnXXWXl131apVu1xTV1eXmpqaJEl1tRuToLvxuQR6G/sadE9bfjZ9ToHewL4G9Db2NXqyHhFshg8fnsmTJ2fcuHGpr69PkqxYsSKzZ8/O448/nk2bNuW73/1uqqqqcuaZZ+7x9adPn77LNQ0NDamvr09NTU1GjBixx++xr2osPQD7DJ9Luop9ja5iX4Pub9iwYaVHAOhU9jWgt7Gv0dN0+2AzefLknH766amqqtrq52PHjs3JJ5+cp59+Otdff31aW1vzL//yL5k0aVLq6urKDAsAAAAAALAXun2wOeCAA3b6+nHHHZeLLrooP/nJT7Jx48bMnDkzH/zgB/foPRoaGna5ZnMEam1tzcqVK/fo+kDlLVu2rPQIAJ3KvgbdU3V1dcff1FyxYkXa2toKTwTw9tjXgN7GvkZXGjp0aMdXqXSGbh9sdseZZ56Z22+/Pe3t7Zk7d+4eB5vNj1nbXT7k0P34XAK9jX0Nur+2tjafVaBXsa8BvY19jZ6mV3zr0sCBA3PggQcmSZqamgpPAwAAAAAAsGd6RbBJss133AAAAAAAAPQUvSLYrFmzJmvXrk2SDBo0qPA0AAAAAAAAe6ZXBJsHHngg7e3tSZIJEyYUngYAAAAAAGDPdOtgs2LFiixatGina55++unccccdSZI+ffrkjDPO6IrRAAAAAAAAOk1tJS/+wgsvZNmyZR1/XrNmTcfxsmXLMmvWrK3WT5kyZas/r1y5MldffXXGjRuX4447Lu985zszcODAJMny5csze/bsPP744x1313z0ox/N4MGDK/PLAAAAAAAAVEhFg82DDz6Yhx56aLuvzZ8/P/Pnz9/qZ38cbDZbsGBBFixYsMP36du3b/7qr/4qZ5555l7PCgAAAAAAUEpFg83bNXr06Hz605/OggUL8uKLL2b16tVZu3ZtWltb079//7zjHe/IxIkT8773va/jzhsAAAAAAICepqLBZsaMGZkxY8Zen7///vvnve99b9773vd24lQAAAAAAADdS3XpAQAAAAAAAPZ1gg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhgg0AAAAAAEBhtZW8+Ouvv56FCxdm4cKFaWxsTGNjY9auXZskOf300zNjxow9ut6cOXPywAMPpLGxMWvWrMmAAQMyZsyYnHnmmTnmmGMq8SsAAAAAAABUXEWDzWWXXdYp12lra8stt9ySmTNnbvXzpqamNDU15cknn8zUqVNz+eWXp7raTUMAAAAAAEDPUtFgs6UhQ4Zk1KhRefbZZ/f43Ntuu60j1hx++OG54IILMnz48Cxfvjx33313Fi1alJkzZ2bAgAG55JJLOnt0AAAAAACAiqposLnooosyZsyYjBkzJnV1dVmxYkU+9alP7dE1li5dmnvuuSdJMmbMmFx99dXp06dPkmTs2LGZNGlSrrrqqjQ2Nuaee+7J1KlTM2LEiE7/XQAAAAAAACqlos8PmzZtWo477rjU1dXt9TXuv//+tLa2JkkuvfTSjlizWd++fXPppZcmSVpbW3Pvvffu9XsBAAAAAACU0K2/8KW9vT1PPvlkkmTUqFEZN27cdteNGzcuI0eOTJI89dRTaW9v77IZAQAAAAAA3q5uHWxWrFiR1atXJ0mOPPLIna6dMGFCkqSpqSkrV66s+GwAAAAAAACdpVsHm8WLF3ccjxo1aqdrN99h88fnAQAAAAAAdHe1pQfYmVWrVnUc19fX73TtkCFDtnvenr7PjtTV1aWmpiZJUl3drTsX7JN8LoHexr4G3dOWn02fU6A3sK8BvY19jZ6sWwebDRs2dBz369dvp2v79u273fN2x/Tp03e5pqGhIfX19ampqcmIESP26Pr7ssbSA7DP8Lmkq9jX6Cr2Nej+hg0bVnoEgE5lXwN6G/saPU23TozNzc0dx7W1O29L++2333bPAwAAAAAA6O669R02ffr06ThuaWnZ6dpNmzZt97zd0dDQsMs1dXV1SZLW1tasXLlyj64PVN6yZctKjwDQqexr0D1VV1d3/E3NFStWpK2trfBEAG+PfQ3obexrdKWhQ4d2fJVKZ+jWwWbLx6Dt6jFnGzdu3O55u2NX34/zx3zIofvxuQR6G/sadH9tbW0+q0CvYl8Dehv7Gj1Nt34k2pYhZdWqVTtd++qrr273PAAAAAAAgO6uWwebQw45pON4yZIlO127dOnS7Z4HAAAAAADQ3XXrYDNs2LAMGjQoSTJv3rydrt38+uDBgzN06NCKzwYAAAAAANBZunWwqaqqyvHHH5/krTtsFixYsN11CxYs6LgDZ9KkSamqquqyGQEAAAAAAN6ubh1skuScc85JdfVbY956661pbm7e6vXm5ubceuutSZKampqce+65XT4jAAAAAADA21FbyYu/8MILWbZsWcef16xZ03G8bNmyzJo1a6v1U6ZM2eYaI0eOzAUXXJA777wzjY2NueKKK/KBD3wgw4cPz/Lly3PXXXdl0aJFSZLzzz8/Bx98cEV+FwAAAAAAgEqpaLB58MEH89BDD233tfnz52f+/Plb/Wx7wSZJLr744rz++uv5zW9+k0WLFuWGG27YZs3UqVNz8cUXv92RAQAAAAAAulxFg01nqa6uzvTp03PCCSfkgQceSGNjY9auXZuDDjooY8aMyVlnnZVjjjmm9JgAAAAAAAB7paLBZsaMGZkxY0anXe/YY4/Nscce22nXAwAAAAAA6A6qSw8AAAAAAACwrxNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAAChNsAAAAAAAACqstPQAAAFDO+nf/RekRepzG0gP0UPs/+5PSIwAAQLfmDhsAAAAAAIDCBBsAAAAAAIDCBBsAAAAAAIDCBBsAAAAAAIDCBBsAAAAAAIDCBBsAAAAAAIDCBBsAAAAAAIDCBBsAAAAAAIDCaksPsDumTZu2W+smTJiQq666qrLDAAAAAAAAdDJ32AAAAAAAABTWI+6w2ezss8/O2WefvcPX+/Xr14XTAAAAAAAAdI4eFWwGDBiQQw89tPQYAAAAAAAAncoj0QAAAAAAAAoTbAAAAAAAAAoTbAAAAAAAAArrUd9hM3v27Dz22GNZuXJlqqurU1dXl3HjxmXKlCmZOHFi6fEAAAAAAAD2So8KNosXL97qz8uWLcuyZcvy8MMP5/jjj8+MGTNywAEH7PF1V61atcs1dXV1qampSZJUV7sxCbobn0ugt7GvAb2NfQ26py0/mz6nQG9gX6Mnq2pvb28vPcSufPSjH81xxx2Xd73rXRk1alT69euXNWvWZO7cufn1r3+dtWvXJkkmTJiQv//7v09t7Z51qGnTpu1yTUNDQ+rr6/dq/n1Z49D3lh6BfcSYlf9eegT2EfY1uop9ja5iX6Or2NcAAGDnesQdNv/4j/+Y/v37b/Pzo48+Ou9///vzta99LYsWLcrcuXPzb//2bznnnHMKTAkAAACUJkTTVYRoADpbjwg224s1m9XV1eVzn/tcPvvZz6a1tTW//OUv9zjYNDQ07HJNXV1dkqS1tTUrV67co+sDlbds2bLSIwB0Kvsa0NvY14Dexr4G3VN1dXWGDRuWJFmxYkXa2toKT0RvNnTo0I6vUukMPSLY7Mrw4cNz9NFHZ86cOVm2bFmampoyePDg3T5/Tx915kMO3Y/PJdDb2NeA3sa+BvQ29jXo/tra2nxW6VF6zbcuHXLIIR3HTU1NBScBAAAAAADYM70m2FRVVZUeAQAAAAAAYK/0mmCzePHijuM9eRwaAAAAAABAab0i2KxYsSL/9//+3yRvfZ+NYAMAAAAAAPQk3T7YPPXUU2ltbd3h66+99lq++c1vpqWlJUnyp3/6p101GgAAAAAAQKeoLT3Artx666357ne/mxNOOCHjxo3LsGHD0qdPn6xZsyZz587Nr3/966xduzZJMn78eMEGAAAAAADocbp9sEmS1atX55e//GV++ctf7nDNCSeckE984hPZb7/9unAyAAAAAACAt6/bB5sZM2Zk7ty5WbBgQZYvX561a9dm/fr16devX+rr6zNu3LhMmTIl48aNKz0qAAAAAADAXun2wWbChAmZMGFC6TEAAAAAAAAqprr0AAAAAAAAAPs6wQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKAwwQYAAAAAAKCw2tIDAAAAAACwfevf/RelR+hxGksP0EPt/+xPSo+wz3OHDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGG1pQfYUytXrswvfvGLPPPMM1m1alVqa2szYsSInHTSSfnTP/3T9O3bt/SIAAAAAAAAe6RHBZunnnoq3/rWt7J+/fqOn23cuDGNjY1pbGzMgw8+mC9+8YsZMWJEwSkBAAAAAAD2TI8JNosWLcoNN9yQ5ubm9OvXLxdeeGEmTpyY5ubmPPLII3nwwQfzyiuv5Gtf+1r++3//79l///1LjwwAAAAAALBbekyw+f73v5/m5ubU1NTk7//+7zNu3LiO1yZOnJiDDz44P/zhD/PKK6/knnvuybRp0wpOCwAAAAAAsPuqSw+wOxYuXJh58+YlSc4444ytYs1m5513XkaNGpUk+cUvfpGWlpYunREAAAAAAGBv9Yhg88QTT3Qcn3HGGdtdU11dndNPPz1Jsm7dujz//PNdMhsAAAAAAMDb1SOCzfz585Mkffv2zejRo3e4bsKECducAwAAAAAA0N31iGCzePHiJMmIESNSU1Ozw3UjR47c5hwAAAAAAIDurrb0ALvS3NyctWvXJknq6+t3uvbAAw9M3759s3HjxqxatWq332N31tbV1XXEourqHtG5YJ/icwn0NvY1oLexrwG9jX0N6G3sa+V1+2CzYcOGjuN+/frtcn2/fv2ycePGrc7blenTp+9yzY9//OMkSU1NTUaMGLHb197XtTz7s9IjsI+oHTGs9AjsI+xrdBX7Gl3FvkZXsa/RVexrdBX7Gl3FvkZXsa+V1+2DTXNzc8dxbe2ux928ZsvzOkNVVVWnXm9fUTvShxzoXexrQG9jXwN6G/sa0NvY12Df0e2DTZ8+fTqOW1padrl+85otz9uVhoaGXa5xOxhdobW1Na+99lqSrR/DB9BT2deA3sa+BvQ29jWgt7Gv0ZN1+2Cz5WPQducxZ5vX7M7j0zbb1XfjQFd57bXXOh7R19DQ4J9NoMezrwG9jX0N6G3sa0BvY1+jJ+v2t4306dMnBx10UJJk1apVO137xhtvZOPGjUlEGAAAAAAAoOfo9sEmSQ455JAkybJly9La2rrDdUuXLt3mHAAAAAAAgO6uRwSbI444IkmycePGvPjiiztcN3fu3G3OAQAAAAAA6O56RLCZPHlyx/FvfvOb7a5pa2vLQw89lCTp379/jjrqqC6ZDQAAAAAA4O3qEcFm7NixOfLII5O8FWwWLFiwzZp77703S5YsSZL8p//0n1JbW9ulMwIAAAAAAOytHhFskuTjH/94+vTpk9bW1lxzzTX5+c9/ngULFuS5557LLbfckh/+8IdJkoMPPjjnn39+4WkBAAAAAAB2X4+5DeXwww/PZz/72XzrW9/K+vXr87/+1//aZs3BBx+cL37xi9l///0LTAgAAAAAALB3qtrb29tLD7EnVq5cmfvvvz/PPPNMmpqaUltbmxEjRuTEE0/M+9///vTt27f0iAAAAAAAAHukxwUbAAAAAACA3qbHfIcNAAAAAABAbyXYAAAAAAAAFCbYAAAAAAAAFCbYAAAAAAAAFCbYAAAAAAAAFCbYAAAAAAAAFCbYAAAAAAAAFCbYAAAAAAAAFCbYAAAAAAAAFCbYAAAAAAAAFFZbegAAAAAAAOhMzc3NWbNmTZJkyJAhhaeB3SPYQAFtbW157rnnMn/+/CxZsiSrVq3Khg0b0tzcnD59+qRfv36pr6/PqFGjcsQRR2TixImprnZDHAAAAADsjqeeeio33nhjqqqqctttt5UeB3aLYANdaMOGDbn33ntz//33Z926dbt9Xv/+/XPOOefkvPPOS79+/So4IQAAAPuidevW5be//W1WrlyZfv365dBDD82RRx6Zqqqq0qMBvC3t7e2lR4DdJthAF1m2bFm+/vWvZ+nSpdu81rdv3/Tt2zf77bdfNm3alI0bN2bjxo0dr69bty4//elP88gjj+SLX/xihg0b1pWjAwDsU1544YXMnj07y5YtS3V1dQ455JCcdNJJOfzww3d57iuvvJJrr702VVVV+da3vtUF0wLs3B/+8If89Kc/zQsvvJDm5uaMGjUq5557bk455ZSONQ899FD++Z//ORs2bNjq3OHDh+fyyy/PxIkTu3psgK3MnTt3j89ZvHhxx/G8efO2CTcTJkx423NBZ6tqlxih4pqbm/N3f/d3eeWVV5Iko0aNymmnnZYJEyZk1KhR6d+//zbnrFu3LkuWLMncuXPz8MMPZ8mSJUmSkSNH5utf/3r69OnTpb8DQKXMmjUrDQ0NblMHimtpacnNN9+cRx55ZLuvn3DCCfmbv/mbDBgwYIfXePnll/P5z38+SfKTn/ykInMC7K5nn3023/jGN9LS0rLNa3/2Z3+Wiy++OM8880yuu+66tLW1bfcatbW1+eIXvyjaAEX9xV/8Radez39/0l25wwa6wC9/+cuOWPORj3wk55133i5vK+/fv3/GjRuXcePG5QMf+EDuueee/OhHP8rSpUvzq1/9Kueff35XjA7QZfwdEqC073znOzuMNUny+OOPZ/78+fnbv/3bjB8/vgsnA9hzb7zxRr71rW91xJqRI0emX79+eemll9LW1paf//znmTx5cr7//e+nra0tJ510UiZNmpQBAwakqakps2fPzpw5c9LS0pKGhob8j//xP/zFQaDX8N+fdFeCDXSBxx57LEnyvve9b69CS1VVVS644IIsW7YsDz74YB577DHBBgCgE73wwgt5+OGHkyQHH3xwPvaxj+Woo47Kpk2bMnfu3Nx1111ZuHBhXnvttVx77bX57Gc/m+OOO67w1AA7NnPmzKxduzbV1dX53Oc+l+OPPz5JsnTp0lx11VV5/fXX853vfCfLly/Ppz71qbz3ve/d6vwpU6bk3nvvzQ9+8IO8+uqreeyxx3L66aeX+FUAOgwcODBTp05NTU3NLte+/PLLefzxx5MkF110UaVHg04h2EAXWL58eZJs9YzgvXHKKafkwQcfzLJlyzpjLAAA/n8zZ85MkgwePDjXXHNNDjzwwCRvfdfg5MmTc/zxx+fee+/Nj3/84zQ3N+f666/PJz/5yW3+BydAdzFnzpwkbz3OcXOsSd660+aCCy7ID37wg7z00kuZPHnyDvey8847L0888UTmz5+fp59+WrABijn//PNz33335fXXX88TTzyRyy67LEceeeROz3n00Uc7gs2f//mfd8WY8LYJNtAFNm3alCRv+/bxzedv7/nDAF3t5ptv7pTrbI7aACXNnz8/yVv/M2BzrNlSVVVVzj///IwZMybXX3991q1bl5tuuinr16/P2Wef3dXjAuzS5i/bPvHEE7d5bdKkSfnBD36QJDnppJN2ep1TTz018+fPz+9///vOHxJgN33kIx/Jqaeemu985zt58cUXc9VVV2XKlCn56Ec/ut1/d4OeSrCBLjBkyJAsXbo0zz//fP7kT/5kr6/z/PPPJ0nq6+s7azSAvfbQQw+VHgGg06xevTpJMm7cuJ2umzBhQq6++upce+21Wb16db73ve9lw4YNueCCC7piTIDdtm7duiRv3Tn4x7b82YgRI3Z6nUMPPTRJ8vrrr3fidAB77rDDDstXv/rV3H///bn99tsza9asPPXUU/nIRz6SM844o/R40CmqSw8A+4JjjjkmSfLzn/88Cxcu3Ktr/O53v8vPf/7zra4HAEDnaG1tTZJUV+/6P5He8Y535Mtf/nKGDRuWJPnRj36U22+/vaLzAeypzU9o2N6+tuXTH3b1JIj9998/SdLc3NyJ0wHsnaqqqpx77rn55je/mWOOOSZvvPFG/vEf/zFXXnllx52F0JO5wwa6wDnnnJOZM2dm/fr1+Yd/+IdMnTo1U6ZMyejRo3f6PwXa2try4osvZtasWZk5c2ZaW1tzwAEH5JxzzunC6QG278ADD8wbb7yRd7/73bnsssv2+jqzZ8/OD3/4w06cDGDPDRgwIE1NTXn11VczevToXa4fNmxYrr766lxzzTVZsmRJfvazn2XDhg3+difQbQwYMCDr16/vuINwb23YsCFJctBBB3XGWACdYsiQIfnCF76QRx99NN///vfzwgsv5O/+7u9y3nnn5aKLLnrbX0sApQg20AWGDBmST33qU7nxxhvT3NycX//61/n1r3+dfv36ZcSIERk8eHD69euX2tratLS0ZMOGDWlqasqyZcs6/uU4eetvPn3605/OkCFDCv42AG8ZO3Zsfvvb32bJkiUZOnToXl/Hf/wD3cEhhxySpqamzJs3L5MnT96tcwYPHtwRbV566aXcd999vuMB6DYGDx6c5cuXZ9WqVdt9ffr06R3rdmbz9w0OHDiwcwcE6AQnn3xy3vOe9+QHP/hBZs6cmbvuuiuPPvpo/vqv/7r0aLBXPBINusikSZPy5S9/OePHj+/42YYNG/LSSy/lmWeeyaOPPpqHH344jz76aJ555pm89NJLW8Wa8ePH5ytf+UqOPfbYEuMDbGPMmDFJkldffTVr1qwpPA3A23PEEUckeeuuv/b29t0+76CDDsqVV17Z8d03zz33XEXmA9hTm+8WfPHFF7f7+pQpUzJlypQccMABO73OCy+8kGTX33UDUMoBBxyQ//Jf/kuuuuqqjBw5MitXrszXv/713HbbbaVHgz0m2EAXOvzwwzu+pPbCCy/M+PHjd/g3yw888MCMHz8+F154Ya699tpcffXVOeyww7p2YICdGDt2bMfx3n4/F0B3sfkvxTQ1NeXxxx/fo3MPOOCAXHHFFTn66KMrMRrAXtn872rPP//8Xl9j06ZNHXvi5rAN0F0deeSRue666/KhD30otbW1HXcIQk/ikWhQwNixY7f6H52bNm3K+vXrs2nTpuy3337Zf//9s99++xWcEGDXttzHGhsb9/oOwBEjRuT000/vrLEA9sro0aMzfvz4rF69Og899FBOPPHEPTq/T58++a//9b/mhhtuyJNPPlmhKQF23zHHHJO///u/T5K0t7enqqpqj6/xyCOPpF+/funXr5+nPQA9Qm1tbaZNm5ZTTz019957b1paWkqPBHukqn1P7vcHAAAAAACg03kkGgAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGGCDQAAAAAAQGG1pQcAAICe5Pnnn8/VV1+dJJkwYUKuuuqqsgN1orVr1+aXv/xl5syZk1deeSVvvvlm2tvbkySf/OQnM2XKlLIDJpk1a1ZuvvnmJMnpp5+eGTNmFJ6o+1qxYkU+9alPJUmGDh2am266qfBE9BYzZszIypUrkyTf/va3M2zYsMITAQD0DoINAEAPdtVVV2Xu3LlJkosuuijTpk0rPBE91fLly3PllVemqampItd/9dVX8/DDD2fevHlZvHhx3njjjbS0tKRv376pq6vL8OHDc9hhh+VP/uRPMnHixPTr168ic7Dvef755/PYY49l4cKFWblyZd58881UV1dn//33z5AhQ3LwwQdnzJgxOfLIIzN69OhUVVWVHhkAgH2UYAMAwD7LHQj/zy233NIRa/r06ZN3vetdGTx4cKqr33qK8qhRo/bqus3Nzbntttty//33p62tbZvX33zzzbz55ptZunRp5syZkySpqanJ1VdfnXHjxu3lb9P7+Gd1zy1evDgNDQ353e9+t81rra2t2bRpU9asWZMXX3wxjzzySJLkHe94R775zW9u93o33XRTHnrooSTd544zAAB6F8EGAAD2catXr85//Md/JEn222+/XHfddTn44IPf9nVbWlpy3XXX5dlnn+34WW1tbcaMGZPhw4enT58+Wb9+fVauXJmXXnopzc3NSd76n+kbNmx42+/PvmvRokX58pe/nHXr1nX8bODAgRkzZkzq6uqSJG+88UZefvnlLFu2rOPRf1uuBwCAribYAADAPm7RokUdx+PHj++UWJMkd955Z0esqaqqygUXXJALL7ww/fv332ZtS0tL5s6dm0cffbTjbgfenmHDhuX2228vPUaXa2lpyf/8n/+zI74MGjQo//k//+dMmjSp446xLa1ZsyZPPvlkHn744axYsaKrxwUAgA6CDQAA7OO2vKtg0KBBnXLNlpaW3HfffR1/njZtWj70oQ/tcH1tbW2OPvroHH300f9fe3ceFVXd/wH8zb4KiGyhsijE6r6BopaSlllqUpqPlqm/tINlHj2eRAtzyad6sjS3Hh9yybJcskxNTBRNUVEzUNZhdUEBZZlhnxnm9wdn7plxFgZZRvP9Oqdz7r1z73e+986dwb6f+/18MGPGDK3p04gMcenSJdy+fRtAU3q/2NhYeHp66tzfwcEBo0ePxujRo3H37t2O6iYRERERkQYGbIiIiIiInnByuVxYbquC6zk5OUIgyMzMDOPGjTP4WG0zcIgMpZqCb+DAgXqDNQ/y8PBojy4RERERERmEARsiIiIiUnPv3j0kJiYiNTUVxcXFkEgksLKygqurK0JDQxEZGdnsAKi24tz19fVITEzE2bNncffuXVRXV8PR0REBAQF4/vnnERgYaHAfb9++jfj4eKSkpKCsrAzm5uZwcXHBgAEDEBkZCRcXF6SlpeHjjz8GAAQHB2PFihXC8YmJidi8ebNam6WlpXjttde0vp8haaWuX7+OEydOICcnB+Xl5bCyskL37t0RHh6OyMhImJu3/T+96+rqcPLkSVy9ehU3b96ERCKBpaUlnJ2dERwcjBEjRsDf31/rsarXR9Xp06eFz05p5MiRiI6OblHfysrKhGV7e3vY2Ni06PiWSk5OxsmTJ1FYWIjKykrY2trCx8cHI0aMwPDhww0ORCkUCly4cAEXL15Ebm4uKioqAABOTk7w8/PDkCFDMGTIkBYFtkpLS3Hy5EmkpqaipKQEVVVVsLe3h5ubG/r06YNRo0bBxcVF67GtuVdLSkowf/58AICrqys2bdqk9RjVtpTHFxUV4fjx40hJScG9e/dgamoKNzc39OvXD+PHj4eDg4PB53/+/HmcPn0a+fn5kEgkcHBwQNeuXTF8+HAMHz4cZmZmWn8zHlZ5ebmw7Orq+tDtKEVHR6O0tFRt2+bNmzU+FwCIiopSu54rVqxAeno6ACA2NhYhISEoLy/HqVOncOnSJdy7dw9isRg2NjbYsWOHRnsymQxJSUm4fPky8vLyIBaLoVAo4ODgAH9/fwwdOhSDBg0y6H7My8tDSkoKsrKycOvWLVRWVkImk8He3h7u7u4ICQnBc889p/NefFi//PILfvjhBwCAqakp5s2bp/H5Ks8zOTkZBQUFqKyshFQqhaWlJRwdHeHh4QE/Pz8MGDAAfn5+bdo/IiIiokcJAzZEREREBABobGzEvn37cOjQIUilUrXXZDIZqqurUVBQgN9//x0vv/wypk6davCg9a1bt/DFF18IaYqU7t+/j6SkJCQlJWkMdOpy5MgRfP/995DJZMK2+vp6VFdXo7CwEMeOHcP8+fPbPUCgJJPJEBcXh4SEBLXtUqkUGRkZyMjIwKlTp7Bs2bIWDXI358qVK/jmm2+EgILq+1ZXV+PmzZuIj49HREQE5s6dCysrqzZ7b0Oo1goRi8Woq6uDtbV1m79PTU0NNm7ciMuXL6ttF4vFSE1NRWpqKs6ePYvFixfD0tJSb1t37tzBV199pVbTR6m4uBjFxcU4d+4cevTogYULF8Ld3b3Z/v388884cOCAxneqsrISlZWVEIlE+PXXX/Hqq69i4sSJzZ9wBzh+/Dh27typ0efCwkIUFhYiISEBMTEx6Nmzp952ampqsG7dOqSmpqptLysrQ1lZGa5du4aEhAQsWrSoTfuv+rv0qNWkuXTpEjZv3qyWhlCXtLQ0bN26FcXFxRqvlZaWorS0FElJSfD398eiRYvg7Oyss62lS5ciNzdX62vKezE7OxuHDh3C1KlTMWHCBMNPSgeFQoFdu3YJqREtLCywcOFCDBw4UG2/oqIifP755xp/H4CmoHRdXR2Ki4uRkpKCAwcOYMOGDZwJRURERP9YDNgQERERERobG/Hll1/i4sWLwjZnZ2f4+fnBwcEBdXV1EIlEKC4uhlwux8GDByEWizF37txm2y4rK8OqVatQXl4OOzs7BAYGwsnJCRKJBNevX0dNTQ0AYP/+/ejWrRuGDh2qs62jR49i586dwrqFhQWCg4PRpUsXVFdXIz09HRKJBOvWrcO0adN0ttO1a1eMHTsWtbW1OHPmDADAxsYGI0aMaPZ8HvTNN9/g9OnTMDExgb+/Pzw9PaFQKCASiVBUVAQAyM/Px6ZNm7B06dIWt69NUlISNmzYINR5MTU1RWBgIDw8PFBXV4eMjAxhlsHZs2dRUlKCjz76SC1g4ezsjLFjxwJoGjC9du0agKZrExoaqvZ+umbp6KMazFAoFDh06JBBAbmWkMvl+OKLL3Dt2jWYm5sjICAA7u7uaGhoQGZmJu7duwcA+Pvvv7Fz50783//9n862bt26hRUrVkAsFgvbvLy84OPjAwAoKCjAjRs3ADTNVFi+fDk+/vhjvbPN4uLiEB8fL6xbW1sjJCQETk5OqKioQFpaGurq6iCVSvHDDz+goqICM2fOVGujLe9VQyQmJuJ///sfAMDT0xM9evSApaUlioqKkJWVBYVCAYlEgs8++wxffvklbG1ttbYjlUqxZs0aiEQiYVvnzp0RFBQEa2trFBcXIzMzE1lZWfjPf/4DNze3NjsH1XvvypUruHXrFrp16/bQ7Y0cOVL4vVIGFXr16qX1s9c3+yMrKwv79u2DXC5Hp06dEBQUhE6dOqGyshIFBQVq+54/fx4bNmwQ0hVaWlrC398frq6uMDU1xZ07d5CdnQ25XA6RSIRly5Zh7dq1cHJy0vreyu+ChYUFunXrBg8PD9ja2kKhUKCiogIikQgSiQRyuRzff/89ALQqaCOXy7FlyxbhnrW1tcWSJUsQHBystl9tbS1WrVqF+/fvA2gKtvn6+qJr166wtrZGfX09ysvLUVBQAIlE8tD9ISIiInpcMGBDRERERNi3b58QrHFycsLs2bMxePBgjRk058+fxzfffIOamhokJCSgV69eegMsAITZBRMmTEBUVJTaTI+qqiqsW7cO169fBwDs2bMH4eHhWmfu3Lp1C7t37xbWe/fujejoaHTu3FnYJpPJsHfvXvzyyy/Ys2ePzj75+/vD398fJSUlwoCivb09Zs+erfdcHiQSiZCeno6ePXti/vz56Nq1q/CaQqHA77//LqQ5unr1KtLT0zUGLFvq7t272Lp1qxCs8fPzw3vvvaf2xHljYyOOHj2K7777DgqFAtnZ2di9ezdmzZol7PPUU08J55uYmCgEbPz8/Fp8HbTx8fGBq6urkErqwIEDKC4uxtixY+Hv798mtXIuXrwIqVSKfv36Ye7cuWozDORyOX744Qf89ttvAIATJ05gwoQJWgMDMpkM69evF4I1jo6OePfdd9G7d2+1/VJSUrBhwwZIJBJUVlZi/fr1WLNmjdZ0d0lJSWrBmmeeeQYzZ85UC3DU1NQgLi4Of/75J4CmgGRQUBCGDBki7NNW96qhtm3bBgcHB8yfPx99+/ZVey09PR2ffvopamtrUV5ejqNHjyIqKkprOwcOHBCCNSYmJpg+fTpefPFFtZlXxcXFWL9+PbKzs7XOanpYgwYNwu+//w4AaGhoQGxsLCZMmICIiAi9s1B0UQYaN23aJARshg8f3uK0bfv27UNjYyOmTJmCCRMmqN03qrOZbt68iU2bNkEul8PExATjx4/HK6+8olHbqbi4GJs2bUJmZibu37+PLVu26AwKDxkyBP3790doaKjWmWaNjY04c+YM4uLiUF9fjx9//BHh4eEPFUirr6/HunXrcPXqVQBNf1NiYmKE4KeqU6dOCcGabt26YfHixVoDYQqFArm5uUhMTGyX9JJEREREjwrT5nchIiIion+ykpISHDx4EEDTQPDKlSt11ugIDw/H4sWLhfV9+/ZBoVDobV8qlWLixIn417/+pZGWy97eHu+9956wvbi4GDk5OVrb2b9/v5AGzdvbG0uWLFEL1gCAubk5pk2bhhdeeEEjnVN7kEqleOqppxAbG6sWrAGaBqnHjRuHsLAwYdu5c+da/Z779+9HXV0dgKYC6cuXL9dID2Rqaorx48djxowZwrb4+PgOTQ9lYmKi9v4KhQJ//vknli9fjjlz5mDt2rXYu3cvrly5gqqqqod6D6lUiqCgICxZskRjIN7MzAzTp08X0nYpFAokJSVpbefPP/9EYWGhcFxMTIxGsAYA+vTpg6VLl8LMzAxA08wpbZ9pY2OjULMDAMLCwvDOO+9ozEaxtbXF/Pnz1VJE7d69WwjGGcuHH36oEawBmmpBvf7668K6rvu5qqpKCJQBwOuvv46XXnpJLVgDNM2EWbZsGVxdXdv0+xoaGooBAwYI6xKJBLt378Y777yD999/Hxs3bsTRo0eRk5MjzGDpCHK5HFOmTMHkyZM1gg4WFhbC8vbt29HQ0AAAmDFjBmbMmKERrAGarl9MTIwwe+jq1atqM5pUzZkzB/3799eZFtDU1BTPPPMM3nnnHaGvx48fb/E5VlVVYfXq1UKwxt3dHStXrtQarAGAjIwMYXnmzJk6Z6yZmJjAz88Pc+bMafMaO0RERESPEgZsiIiIiJ5wR48eFQaIo6Kimq0NEBoaij59+gAAbt++3eyT8Q4ODjqfwgeanr7u37+/sK4tYFNVVYXk5GRhfcaMGXrrkUydOrXDathMmzZNb22WZ599VljWFYwyVHV1Nc6fPy+sT58+XWdKKgAYN24cunfvDqApYHHixIlWvX9LhYWFYd68eWqD0UDTAPrVq1exf/9+fPrpp5g9ezZiYmJw7NgxYaDaUG+++aYQQHmQiYmJQddf9bqMGTMGvr6+Ot/Pz88Po0ePFta1DWqnpqYKwTFzc3PMmjVL54wiExMTzJkzRziH4uJijZovHWn06NHw9vbW+frIkSOFvhYVFQkpDVWdPXtWCMC4urripZde0tmenZ1dm6fKA4AFCxZg8ODBatsUCgWKiopw5swZ7NixAzExMZg5cya++uorYZZfe+rcuXOzacYKCgqEvvj6+uLFF1/Uu7+1tTUmT54srCtnaz2ssLAw4fdMOevOUGVlZYiNjUVWVhaApsD6qlWr9P5Nqa2tFZbbssYXERER0eOKARsiIiKiJ5zySWgAiIiIMOgY1RonmZmZevcdMGBAs8XeVQfIlSm0VGVnZwuza5ycnDRqrDzIxsYGgwYN0rtPW7CwsFB7kl+b5s6tJbKysoSB8E6dOjX73qampmoBi7S0tFa9/8MYNWoU1q9fjzFjxmidJQA0DaTn5OTg22+/xYIFCwzup7u7O3r06KF3H9Un+7Vd/9raWrVi7KrXS5dRo0YJy7m5ucKMJyXVwf9+/frprCui5OzsrDajxRifk1J4eLje121sbIQaMQqFQqiNoio9PV2tPV0BNaWwsDCNoF5rWVtbY/Hixfjggw/Qu3dvnQGz+vp6JCUlYeXKlfjss88eeraXIcLCwpq9Fqq/x8OGDTModaDq76EyWKJPYWEhEhMTsW/fPuzatQtxcXHCf9u3bxfe88aNGwbP9ioqKsKHH36ImzdvAgCCgoLw8ccfN3vvd+nSRVj+448/DHovIiIion8yJn8lIiIieoJJJBLcuXMHQNNMgH379hl03K1bt4RlZf0BXby8vJptz97eXljW9sS+akHunj17aqRW0sbPz0+o+dFePD09m62noHpuqk+TPwzV6+Dn59fs4C8ABAQECMv5+flQKBRtUj+mJVxcXDBnzhzMnDkTIpEIGRkZyM3NRV5ensb9c//+faxevRoffPCBMJNLF+XsIX06deokLGu7/oWFhcKgtLW1td7ZJUo+Pj6wsrJCfX09GhsbUVhYqHGdlVS36xMQEIArV65oHN/RDPm+ql7T5r6v/v7+zbZnZWUFLy8vtcBZW+nfvz/69+8PsViMtLQ0ZGVlIT8/H/n5+RqBtsuXLyM2NharV69ulxl6zQUXgabgtFJaWppBQV7VtJTaAmhKiYmJOHjwoPCb3xy5XI6amhq13zBt8vLysG3bNkgkEgBNQfqFCxc2G6gHmgJ6p06dAtAUsMnLy8PIkSPRt2/fZmd7EhEREf0TMWBDRERE9AQrLy8XlmUymVqRdENVV1frfV1fyi4l1aCHtpoSymLwgPoT2foYul9rtMW5tYTqdXB1dTXoGNX9ZDIZamtrDep3ezA3N0dQUBCCgoKEbSUlJbhw4QKOHj2KsrIyAE3XaePGjdi4caNG3SNVhpyHalDLkHvLkGCWqakpXFxchAL0yoFqbW0aWm9Dtbj7g+11pPa4poZwdnZul4CNkoODA8LDw4UZRHK5HCKRCKdOncKZM2eE87h58yb27NmDWbNmtUsfmqP6m6w628ZQ2n6PFQoFtmzZgsTExBa3V1dX12zAZv369cL1i4iIQHR0tEHBZADo27cvnn/+eRw7dgxA04w15X3g6OiIwMBAhISEYNCgQR3ym05ERERkbAzYEBERET3BtD0d31LNBSHaYjaH6pPw+gbwVemrK9NWOnqmSltch7q6OqMFbLRxc3PDyy+/jMjISHz++edCOrDKykokJSXpTVHW1vdWS+4Z1ev/4Mydh2lTX3sd6Z/8fVVlZmaGwMBABAYGYtSoUVizZo3Q74SEBEyfPt2gGSItYUh7rf1N1pbCLCEhQS1Y07dvXwwbNgy+vr7o0qULrKys1ALL0dHRwsweQ1KimZmZCX8HSktL0dDQ0KIZSrNmzUJoaCh++eUXtTpTlZWVuHjxIi5evIjt27dj8ODBePPNNw0OghIRERE9jhiwISIiInqCqQ6S2tjYYOfOnUbsjW6q/ayvrzfoGEP3e5w8zHV4MO1TRw+MG8rW1hbz589HdHS0MEickZFhUE2Z1lC9Hg9eK31Ur/+Dg9MP06a+9h431tbWQuDhcfi+BgQEYNKkSdizZw8AQCqVIicnB8HBwR3eF9UA1+LFizF48OBWt/nbb78Jy6+99hqioqL07t/SgOH777+PLVu2QCKRICsrC5988gliYmJadB8PHjwYgwcPxr1794TUdZmZmUL6TYVCgYsXLyI9PR2rVq2Cp6dni/pIRERE9LhoPvk3EREREf1jOTo6Csu1tbWPbJBDtWaGMm1Wc5qrrfM4Uk2ppK9WhSrVGhjm5uaPdDCgS5cuanVpKioq2v09Va9pWVmZWj0QXRobG9Wuv+r9+WCbhn5OJSUlOtt73DyO39e+ffuqrXfEvaeNk5NTm/bh3r17Qs0aOzs7TJw4Ue/+NTU1zaa5fJCXlxdiY2OFz10ZtHmYmWIuLi4YOXIk3n77baxbtw6bN2/Ga6+9JgSyJBIJdu3a1eJ2iYiIiB4XDNgQERERPcE6d+6sVhcgKyvLiL3RzcfHR1jOzc01aFBdNbWOLh2d0qy1VK9DTk6OQemKVD9TX1/fR/6cLSwshGXVNE3txdvbG6amTf9bVFtbixs3bjR7TGFhoRDcNDU1hbe3t9rrvr6+wrJqEXl9VPdTPV7pUf/cVKnepyKRqNn9GxoacPPmzXbsUfNU7ztA+73XEZ+Bn5+fsJyZmdnq9lRr4nh6ejb7ncrMzDTo9/VBXl5e+Oijj9okaKPKxcUFUVFRePvtt4VtKSkpkEqlrWqXiIiI6FHFgA0RERHRE27AgAHC8vHjx43YE90CAgKEgcby8nJcv35d7/51dXW4dOlSs+2qDtLKZLLWdbIDBAQECH0Wi8X466+/9O7f2NioVrsiJCSkPbvXalKpFLdv3xbWO6JWhY2NDXr27CmsG1KY/eTJk8Kyn5+fRpq50NBQYfnq1auorKzU215ZWZlagXnV45Uep3tV9T47f/58s3WuLly4gIaGhvbull6FhYVq69ruvY74DFR/j5OTk1s9y0Y1yGTINW7N3wBvb2+tQZuWpBrUZeDAgcKyXC5HVVVVq9skIiIiehQxYENERET0hBs/frwwwyA5OdmgAWuljkobZG9vrzZgt3v3br1PWP/0008GFe+2s7MTBjTFYvEjPxBuZ2eH8PBwYf27777T+wT7sWPHhBkjJiYmiIyMbPc+KolEIhw6dKhFafZ+/fVXtfN5ME1Ve1G9LvHx8RqD96ry8vJw4sQJYf25557T2Kd3795wc3MD0BSE2rFjh872FAoFvv32WyGo4e7ujl69emns9zjdqxEREUJwo6SkBEeOHNG5b01NDX766ac2ff/Dhw8jNTXV4P3r6+tx8OBBYd3R0VFtlpDSw6R6ayk/Pz8h4NXQ0ICNGzca/FnLZDKNQIabm5tw39y4cQPFxcU6j09KSmo2CNwcbUGbNWvW6AzaiMVig9pVTZlnYmLy2KcNJCIiItKFARsiIiKiJ5yHhwdeeeUVYX3Lli3YtWuXzoE0uVyOlJQUfP3111iyZElHdROvvvqqMMsmPz8fn332mUbASCaT4ccff8SRI0c0UhxpY2FhgaeeegpA03kZMivH2KKiooQZHXfu3MGaNWs0BmEbGxtx9OhRtVoPY8eOFYIIHaG6uhq7d+9GdHQ0du7ciby8PJ2plsRiMXbs2IG9e/cK23x9fTssYDN8+HAhrZlMJsOaNWu0zuJKTU3FJ598IgRXfH19MWzYMI39TE1NMW3aNGH93Llz2Lp1q8agdW1tLTZv3ozk5GRh2/Tp04UAqqrH6V61t7fH+PHjhfXvv/8eR44c0UjhV1JSgjVr1qC0tNSg76uhcnJysHr1aixduhTx8fF6A8sikQgrVqxQS4U3YcIErZ+Ban2ly5cvt1vQ7K233hK+46mpqYiNjdWbWq6oqAj79+9HdHS0RlpLBwcH+Pv7A2gKDq5btw5FRUVq+zQ2NuLYsWP4+uuvYWpq2urPoiVBm+XLl2P9+vW4evWqzutZVFSETZs2Ceu9evXqkHSJRERERMbAf+UQERER/UP88ccfLRrEnTJlijBr5dVXX0VpaSlOnz4NhUKBw4cP49ixY+jZsyfc3d1haWmJ2tpalJaWqtXv6MinnLt3745p06YJQYiUlBRER0cjODgYLi4uqK6uRnp6OsRiMczNzfH6668L++qrPTFkyBDh6foNGzYgMTERHh4eMDMzE/Z544032vHMWsbDwwPz5s3Dhg0b0NjYiOzsbLz//vsIDAyEu7s76urqkJmZqTYDwN/fH9OnTzdKf8ViMY4cOYIjR47A1tYWPXr0gJOTE2xsbFBXV4e7d+8iLy9PLW2Wo6Mj3n33Xa2D5u3B3NwcCxYswIoVKyAWi1FRUYGVK1fC29tbmGlRUFCgNvPG0dERCxYs0DlwPHToUGRkZCA+Ph5AUxq18+fPIyQkBI6OjqisrMT169fVBrHHjRuHIUOG6Ozn43SvRkVF4dq1a8jJyYFCocDOnTvx22+/ITAwENbW1igpKUFGRgbkcjmefvppuLm54ezZswDarlZMbm4ucnNzERcXB3d3d3Tv3h2dOnWCmZkZxGIxCgoKUFJSonbM4MGD8cILL2htr1+/frC0tERDQwMKCgqwcOFCBAcHw87OTtinT58+6NOnT6v67eXlhQULFuCrr75CfX09RCIRli1bBnd3d/j6+sLe3h5SqRRisRiFhYXNzvaZMmUKVq9eDYVCgfz8fCxatAiBgYFwc3MTfi+UtW6mTp2KhIQElJaWtuocvL298eGHH2LVqlWQSCRCerSYmBi1FIJyuRznzp3DuXPnYGlpCW9vb7i5ucHGxgbV1dUoKSlBbm6usL+lpSVmzJjRqr4RERERPcoYsCEiIiL6h6isrGy2VoYq1dQ5JiYmiI6ORo8ePbB3715UV1dDJpMhKytL44lt1WMCAgJa3e+WGD9+PBobG7Fnzx7I5XJIpVKkpKSo7WNra4v58+fD0tJS2GZjY6OzzQkTJiA5ORm3b9+GXC5XqyWi9KgMgisNHToUVlZW2Lp1KyorKyGXy5GWloa0tDSNfYcNG4Z58+apXY+O4ObmhuDgYGRmZgozK2pqapqtP9SvXz/Mnj27Q2cDAUC3bt2wcuVKrF+/Hvn5+QCa6ppoS4/m6+uLhQsXwsPDQ2+bs2fPhpOTE37++WdIpVLU1tbi8uXLGvtZWFggKioKkyZN0tve43SvWlhYYNmyZfjiiy+Ez7ysrAxJSUlq+wUEBGDRokVqs8H0fV8NERoaipycHLVgTHFxsd50YJaWlpg4cSImTZqkFgBTZWtrizfeeANxcXFQKBRa27S2tm51wAZoqmWzatUqbN26FXl5eQadg6urK7p06aKxvVevXpg9eza2b98OuVyu9ffCxMQEkydPxqRJk5CQkNDq/gOAj4+PWtAmMzNTI2ijGrxpaGiASCTSOZvIzc0N7777rjAbjoiIiOifiAEbIiIiIhK88MILeOaZZ3DmzBmkpqaisLAQYrEYUqkU1tbW6NKlC7p164aQkBD069evQ4rCP+jll19G//79ER8fj9TUVNy/fx8WFhZwcXHBgAEDEBkZCRcXF7WBYdUn4B9ka2uLtWvXIj4+Hn/99Rdu376N6urqZgulG9uAAQPw9ddf4+TJk/jrr79w8+ZNSCQSWFpaonPnzggJCcHIkSOFdEgdzdPTU5ixkp6ejszMTNy4cQN3796FRCKBVCqFlZUV7Ozs0LVrV/Ts2RNDhw6Fl5eXUfqr7PPatWtx4cIFXLx4ETk5OUJqQGVqqbCwMAwZMsTgWSCTJ0/GiBEjkJCQgJSUFJSUlKCmpga2trZwd3dHnz59MHr0aIO+S4/bvWpnZ4ePPvoISUlJOHPmDPLy8lBVVYVOnTqhW7duGD58OCIiImBubq4WQNb3fTVEZGQkIiMjcePGDaSnp0MkEqGoqAilpaWoqamBQqGAjY0NnJyc4OXlhdDQUISFhcHe3r7ZtseMGQMvLy+cOHECIpEIZWVlaGho0JnurzV8fHzw73//GykpKbh06RKysrJQVlaGmpoamJubw8HBAZ6envD390efPn3w9NNP67wvx4wZg8DAQBw+fBhpaWkoLy+HpaUlnJ2dERoaimeffRa+vr7tcg76gjaff/45srOzkZaWhpycHBQVFaG8vBz19fWwsrKCk5MTfHx8MHDgQISHh7dp6jwiIiKiR5GJoj3+ZUlEREREZGR79uwR0kdNmzYNEydONG6HiEinuXPnCmm5/vvf/8LJycm4HSIiIiIiMoKOSQpNRERERNSBFAoFLly4IKz7+fkZsTdEpI9qDZUuXbowWENERERETywGbIiIiIjoH+fIkSO4c+cOAMDZ2RnBwcFG7hERaSOTybBz505hPSIiwoi9ISIiIiIyLtawISIiIqLHxoULF5CdnY3IyEh4enpqvF5TU4NDhw4JqdAAYPz48TA15XNKRB1t27Zt6NmzJ8LDw2FjY6Px+o0bN7Bt2zbk5uYCaCpAP3bs2I7uJhERERHRI4MBGyIiIiJ6bNTV1eHw4cM4fPgwPDw84OXlBQcHB8hkMty7dw8ikQj19fXC/qGhoRg3bpwRe0z05Lp9+zb++OMPxMXFwcfHBx4eHrC2tkZtbS0KCwtx69YtKEuqmpiY4K233oKLi4uRe01EREREZDwM2BARERHRY+nu3bu4e/eu1tdMTEwQERGBuXPncnYNkZFJpVKIRCKIRCKtr9vZ2WH27NlMh0ZERERETzwThfKRJiIiIiKiR5xMJkNqair+/vtv5Ofno6KiAhKJBPX19bC1tYWLiwuCgoIwcuRI+Pr6Gru7RE+0iooKJCcnIz09HUVFRZBIJJBIJACATp06oXv37ujduzeeffZZ2NnZGbm3RERERETGx4ANERERERERERERERGRkTE/BBERERERERERERERkZExYENERERERERERERERGRkDNgQEREREREREREREREZGQM2RERERERERERERERERsaADRERERERERERERERkZExYENERERERERERERERGRkDNgQEREREREREREREREZGQM2RERERERERERERERERsaADRERERERERERERERkZExYENERERERERERERERGRkDNgQEREREREREREREREZGQM2RERERERERERERERERsaADRERERERERERERERkZExYENERERERERERERERGRkDNgQEREREREREREREREZGQM2RERERERERERERERERsaADRERERERERERERERkZExYENERERERERERERERGRkDNgQEREREREREREREREZ2f8D9S0QvYdfluIAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 1000x500 with 1 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 445,
       "width": 822
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline \n",
    "%config InlineBackend.figure_format = 'retina'\n",
    "plt.style.use('ggplot')\n",
    "plt.rcParams['figure.figsize'] = (10,5)\n",
    "\n",
    "kobe_streak_counts = pd.Series(kobe_streak).value_counts().sort_index()\n",
    "\n",
    "kobe_streak_counts.plot(kind = 'bar', color = 'crimson')\n",
    "plt.xlabel('Length of Shooting Streaks')\n",
    "plt.show(); "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that instead of making a histogram, we chose to make a bar plot. A bar plot is preferable here since our variable is discrete – counts – instead of continuous."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class = 'exercise'>\n",
    "<h4>Exercise 2</h4>\n",
    "Describe the distribution of Kobe's streak lengths from the 2009 NBA finals. What was his typical streak length? How long was his longest streak of baskets?\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Compared to What?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We've shown that Kobe had some long shooting streaks, but are they long enough to support the belief that he had hot hands? What can we compare them to?\n",
    "\n",
    "To answer these questions, let's return to the idea of *independence*. Two processes are independent if the outcome of one process doesn't effect the outcome of the second. If each shot that a player takes is an independent process, having made or missed your first shot will not affect the probability that you will make or miss your second shot.\n",
    "\n",
    "A shooter with a hot hand will have shots that are *not* independent of one another. Specifically, if the shooter makes his first shot, the hot hand model says he will have a *higher* probability of making his second shot.\n",
    "\n",
    "Let's suppose for a moment that the hot hand model is valid for Kobe. During his career, the percentage of time Kobe makes a basket (i.e. his shooting percentage) is about 45%, or in probability notation,\n",
    "#### `P(shot 1 = H) = 0.45`\n",
    "If he makes the first shot and has a hot hand (not independent shots), then the probability that he makes his second shot would go up to, let's say, 60%,\n",
    "#### `P(shot 2 = H|shot 1 = H) = 0.60`\n",
    "As a result of these increased probabilites, you'd expect Kobe to have longer streaks. Compare this to the skeptical perspective where Kobe does *not* have a hot hand (where each shot is independent of the next). If he hit his first shot, the probability that he makes the second is still 0.45.\n",
    "#### `P(shot 2 = H|shot 1 = H) = 0.45`\n",
    "In other words, making the first shot did nothing to effect the probability that he'd make his second shot. If Kobe's shots are independent, then he'd have the same probability of hitting every shot regardless of his past shots: 45%.\n",
    "\n",
    "Now that we've phrased the situation in terms of independent shots, let's return to the question: how do we tell if Kobe's shooting streaks are long enough to indicate that he has hot hands? We can compare his streak lengths to someone without hot hands: an independent shooter."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Simulations in Python"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "While we don't have any data from a shooter we know to have independent shots, that sort of data is very easy to simulate in Python. In a simulation, you set the ground rules of a random process and then the computer uses random numbers to generate an outcome that adheres to those rules. As a simple example, you can simulate flipping a fair coin with the following:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tails\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "\n",
    "outcomes = np.array(['heads', 'tails'])\n",
    "print(np.random.choice(outcomes))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The vector `outcomes` can be thought of as a hat with two slips of paper in it: one slip says `heads` and the other says `tails`. The function `random.choice()` draws one slip from the hat and tells us if it was a head or a tail.\n",
    "\n",
    "Run the second command listed above several times. Just like when flipping a coin, sometimes you'll get a heads, sometimes you'll get a tails, but in the long run, you'd expect to get roughly equal numbers of each.\n",
    "\n",
    "If you wanted to simulate flipping a fair coin 100 times, you could either run the function 100 times or, more simply, adjust the `size` argument, which governs how many samples to draw (the `replace = True` argument indicates we put the slip of paper back in the hat before drawing again). Save the resulting array of heads and tails in a new object called sim_fair_coin."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "sim_fair_coin = np.random.choice(outcomes, size = 100, replace = True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To view the results of this simulation, use `unique()` to count up the number of heads and tails:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[['heads' '48']\n",
      " ['tails' '52']]\n"
     ]
    }
   ],
   "source": [
    "unique, counts = np.unique(sim_fair_coin, return_counts = True)\n",
    "print(np.asarray((unique, counts)).T)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Since there are only two elements in `outcomes`, the probability that we \"flip\" a coin and it lands heads is 0.5. Say we're trying to simulate an unfair coin that we know only lands heads 20% of the time. We can adjust for this by adding an argument called `p`, which provides the probabilities associated with each entry in the array."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[['heads' '22']\n",
      " ['tails' '78']]\n"
     ]
    }
   ],
   "source": [
    "sim_unfair_coin = np.random.choice(outcomes, size = 100, replace = True, p = [0.2, 0.8])\n",
    "unique, counts = np.unique(sim_unfair_coin, return_counts = True)\n",
    "print(np.asarray((unique, counts)).T)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`p = [0.2, 0.8]` indicates that for the two elements in the outcomes array, we want to select the first one, `heads` with probability 0.2 and the second one, `tails` with probability 0.8. Another way of thinking about this is to think of the outcome space as a bag of 10 chips, where 2 chips are labeled \"head\" and 8 chips \"tail\". Therefore at each draw, the probability of drawing a chip that says \"head\" is 20%, and \"tail\" is 80%."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class = 'exercise'>\n",
    "<h4>Exercise 3</h4>\n",
    "In your simulation of flipping the unfair coin 100 times, how many flips came up heads?\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In a sense, we've shrunken the size of the slip of paper that says \"heads\", making it less likely to be drawn and we've increased the size of the slip of paper saying \"tails\", making it more likely to be drawn. When we simulated the fair coin, both slips of paper were the same size. This happens by default if you don't provide a `p` argument; all elements in the `outcomes` array have an equal probability of being drawn."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Simulating the Independent Shooter"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Simulating a basketball player who has independent shots uses the same mechanism that we use to simulate a coin flip. To make a valid comparison between Kobe and our simulated independent shooter, we need to align both their shooting percentage and the number of attempted shots."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "outcomes = np.array(['H', 'M'])\n",
    "sim_basket = np.random.choice(outcomes, size = 133, replace = True, p = [0.45, 0.55])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "With the results of the simulation saved as `sim_basket`, we have the data necessary to compare Kobe to our independent shooter. We can look at Kobe's data alongside our simulated data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['H' 'M' 'M' 'H' 'H' 'M' 'M' 'M' 'M' 'H' 'H' 'H' 'M' 'H' 'H' 'M' 'M' 'H'\n",
      " 'H' 'H' 'M' 'M' 'H' 'M' 'H' 'H' 'H' 'M' 'M' 'M' 'M' 'M' 'M' 'H' 'M' 'H'\n",
      " 'M' 'M' 'H' 'H' 'H' 'H' 'M' 'H' 'M' 'M' 'H' 'M' 'M' 'H' 'M' 'M' 'H' 'M'\n",
      " 'H' 'H' 'H' 'H' 'M' 'H' 'H' 'H' 'M' 'H' 'M' 'H' 'M' 'M' 'M' 'M' 'M' 'M'\n",
      " 'H' 'M' 'H' 'M' 'M' 'M' 'M' 'H' 'H' 'M' 'M' 'H' 'M' 'H' 'H' 'M' 'H' 'M'\n",
      " 'M' 'M' 'H' 'M' 'M' 'M' 'M' 'H' 'M' 'H' 'M' 'M' 'H' 'M' 'M' 'H' 'H' 'M'\n",
      " 'M' 'M' 'M' 'H' 'H' 'H' 'M' 'M' 'H' 'M' 'M' 'H' 'M' 'H' 'H' 'M' 'H' 'M'\n",
      " 'M' 'H' 'M' 'M' 'M' 'H' 'M']\n"
     ]
    }
   ],
   "source": [
    "print(np.array(kobe['basket']))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['H' 'H' 'M' 'H' 'M' 'M' 'H' 'M' 'H' 'M' 'M' 'H' 'H' 'M' 'H' 'H' 'M' 'M'\n",
      " 'H' 'M' 'H' 'H' 'M' 'H' 'H' 'H' 'H' 'H' 'H' 'M' 'H' 'M' 'M' 'M' 'M' 'H'\n",
      " 'M' 'H' 'H' 'M' 'H' 'M' 'M' 'M' 'M' 'H' 'M' 'M' 'H' 'M' 'M' 'M' 'M' 'H'\n",
      " 'H' 'H' 'M' 'M' 'M' 'M' 'M' 'H' 'M' 'M' 'H' 'H' 'M' 'M' 'M' 'H' 'H' 'H'\n",
      " 'H' 'H' 'M' 'H' 'M' 'H' 'M' 'M' 'M' 'M' 'H' 'M' 'M' 'H' 'H' 'H' 'M' 'M'\n",
      " 'M' 'H' 'M' 'H' 'H' 'H' 'H' 'M' 'H' 'M' 'H' 'H' 'H' 'M' 'H' 'M' 'H' 'M'\n",
      " 'M' 'M' 'H' 'H' 'M' 'M' 'H' 'M' 'H' 'H' 'M' 'H' 'M' 'M' 'M' 'M' 'H' 'H'\n",
      " 'M' 'M' 'M' 'M' 'H' 'M' 'H']\n"
     ]
    }
   ],
   "source": [
    "print(sim_basket)  # Notice that every time you rerun the simulation above, the results of sim_basket changes."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Both data sets represent the results of 133 shot attempts, each with the same shooting percentage of 45%. We know that our simulated data is from a shooter that has independent shots. That is, we know the simulated shooter does not have a hot hand."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## On Your Own\n",
    "### Comparing Kobe Bryant to the Independent Shooter\n",
    "\n",
    "Using `calc_streak()`, compute the streak lengths of `sim_basket`.\n",
    "\n",
    "<ol>\n",
    "  <li>Describe the distribution of streak lengths. What is the typical streak length for this simulated independent shooter with a 45% shooting percentage? How long is the player's longest streak of baskets in 133 shots?</li><br>\n",
    "  <li>If you were to run the simulation of the independent shooter a second time, how would you expect its streak distribution to compare to the distribution from the question above? Exactly the same? Somewhat similar? Totally different? Explain your reasoning.</li><br>\n",
    "      <li>How does Kobe Bryant's distribution of streak lengths compare to the distribution of streak lengths for the simulated shooter? Using this comparison, do you have evidence that the hot hand model fits Kobe's shooting patterns? Explain.</li>\n",
    "</ol>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<div class = \"license\">\n",
    "This lab was adapted by David Akman and Imran Ture from OpenIntro by Andrew Bray and Mine Çetinkaya-Rundel.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.11.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}