Codecademy_NBA_Trends_Project.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Codecademy [NBA Trends Project](https://www.codecademy.com/projects/practice/nba-trends)\n",
    "\n",
    "*Analyze National Basketball Association (NBA) data to look at associations between teams, win-rates, playoff appearances, and more.*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this project, you'll analyze data from the NBA (National Basketball Association) and explore possible associations. \n",
    "\n",
    "This data was originally sourced from 538's Analysis of the [Complete History Of The NBA](http://fivethirtyeight.com/interactives/the-complete-history-of-every-nba-team-by-elo) and contains the original, unmodified data from [Basketball Reference](http://www.basketball-reference.com/) as well as several additional variables 538 added to perform their own analysis. \n",
    "\n",
    "You can read more about the data and how it's being used by 538 [here](https://github.com/fivethirtyeight/data/tree/master/nba-elo). For this project we've limited the data to just 5 teams and 10 columns (plus one constructed column, `point_diff`, the difference between `pts` and `opp_pts`).\n",
    "\n",
    "You will create several charts and tables in this project, so you'll need to use `plt.clf()` between plots in your code so that the plots don't layer on top of one another.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "from scipy.stats import pearsonr, chi2_contingency\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "#to make the output look nicer\n",
    "np.set_printoptions(suppress=True, precision = 2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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>game_id</th>\n",
       "      <th>year_id</th>\n",
       "      <th>fran_id</th>\n",
       "      <th>opp_fran</th>\n",
       "      <th>game_location</th>\n",
       "      <th>is_playoffs</th>\n",
       "      <th>pts</th>\n",
       "      <th>opp_pts</th>\n",
       "      <th>game_result</th>\n",
       "      <th>forecast</th>\n",
       "      <th>point_diff</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>194611010TRH</td>\n",
       "      <td>1947</td>\n",
       "      <td>Knicks</td>\n",
       "      <td>Huskies</td>\n",
       "      <td>A</td>\n",
       "      <td>0</td>\n",
       "      <td>68</td>\n",
       "      <td>66</td>\n",
       "      <td>W</td>\n",
       "      <td>0.359935</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>194611020CHS</td>\n",
       "      <td>1947</td>\n",
       "      <td>Knicks</td>\n",
       "      <td>Stags</td>\n",
       "      <td>A</td>\n",
       "      <td>0</td>\n",
       "      <td>47</td>\n",
       "      <td>63</td>\n",
       "      <td>L</td>\n",
       "      <td>0.368899</td>\n",
       "      <td>-16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>194611020PRO</td>\n",
       "      <td>1947</td>\n",
       "      <td>Celtics</td>\n",
       "      <td>Steamrollers</td>\n",
       "      <td>A</td>\n",
       "      <td>0</td>\n",
       "      <td>53</td>\n",
       "      <td>59</td>\n",
       "      <td>L</td>\n",
       "      <td>0.359935</td>\n",
       "      <td>-6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>194611050BOS</td>\n",
       "      <td>1947</td>\n",
       "      <td>Celtics</td>\n",
       "      <td>Stags</td>\n",
       "      <td>H</td>\n",
       "      <td>0</td>\n",
       "      <td>55</td>\n",
       "      <td>57</td>\n",
       "      <td>L</td>\n",
       "      <td>0.620204</td>\n",
       "      <td>-2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>194611070STB</td>\n",
       "      <td>1947</td>\n",
       "      <td>Knicks</td>\n",
       "      <td>Bombers</td>\n",
       "      <td>A</td>\n",
       "      <td>0</td>\n",
       "      <td>68</td>\n",
       "      <td>63</td>\n",
       "      <td>W</td>\n",
       "      <td>0.339290</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        game_id  year_id  fran_id      opp_fran game_location  is_playoffs  \\\n",
       "0  194611010TRH     1947   Knicks       Huskies             A            0   \n",
       "1  194611020CHS     1947   Knicks         Stags             A            0   \n",
       "2  194611020PRO     1947  Celtics  Steamrollers             A            0   \n",
       "3  194611050BOS     1947  Celtics         Stags             H            0   \n",
       "4  194611070STB     1947   Knicks       Bombers             A            0   \n",
       "\n",
       "   pts  opp_pts game_result  forecast  point_diff  \n",
       "0   68       66           W  0.359935           2  \n",
       "1   47       63           L  0.368899         -16  \n",
       "2   53       59           L  0.359935          -6  \n",
       "3   55       57           L  0.620204          -2  \n",
       "4   68       63           W  0.339290           5  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nba = pd.read_csv('./data_files/nba_games.csv')\n",
    "nba.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Subset Data to 2010 Season, 2014 Season\n",
    "nba_2010 = nba[nba.year_id == 2010]\n",
    "# print(nba_2010)\n",
    "nba_2014 = nba[nba.year_id == 2014]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 1\n",
    "The data has been subset for you into two smaller datasets: games from 2010 (named nba_2010) and games from 2014 (named nba_2014). To start, let’s focus on the 2010 data.\n",
    "\n",
    "Suppose you want to compare the knicks to the nets with respect to points earned per game. Using the pts column from the nba_2010 DataFrame, create two series named knicks_pts (fran_id = \"Knicks\") and nets_pts(fran_id = \"Nets\") that represent the points each team has scored in their games.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "21719     93\n",
      "21725    100\n",
      "21728    127\n",
      "21734    117\n",
      "21739     89\n",
      "        ... \n",
      "22102    105\n",
      "22110    103\n",
      "22116     98\n",
      "22118    114\n",
      "22126    113\n",
      "Name: pts, Length: 82, dtype: int64\n"
     ]
    }
   ],
   "source": [
    "knicks_pts = nba_2010.pts[nba_2010.fran_id == 'Knicks']\n",
    "print(knicks_pts)\n",
    "nets_pts = nba_2010.pts[nba_2010.fran_id =='Nets']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 2\n",
    "\n",
    "Calculate the difference between the two teams’ average points scored and save the result as diff_means_2010. Based on this value, do you think fran_id and pts are associated? Why or why not?\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Average Knicks Points in 2010: 102.10975609756098\n",
      "Average Nets Points in 2010: 92.3780487804878\n",
      "Difference between average points scored by Nets and Knicks in 2010: 9.731707317073173\n"
     ]
    }
   ],
   "source": [
    "avg_knicks = knicks_pts.mean()\n",
    "print(f\"Average Knicks Points in 2010: {avg_knicks}\")\n",
    "avg_nets = nets_pts.mean()\n",
    "print(f\"Average Nets Points in 2010: {avg_nets}\")\n",
    "diff_means_2010 = avg_knicks - avg_nets\n",
    "print(f\"Difference between average points scored by Nets and Knicks in 2010: {diff_means_2010}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 3\n",
    "Rather than comparing means, it’s useful look at the full distribution of values to understand whether a difference in means is meaningful. Create a set of overlapping histograms that can be used to compare the points scored for the Knicks compared to the Nets. Use the series you created in the previous step (1) and the code below to create the plot. Do the distributions appear to be the same?\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGdCAYAAADqsoKGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA5z0lEQVR4nO3de1xVdb7/8feW2yYVvOBwKUTQJklJazNjMINaIYZNTg1NaI1pqTOcbAw5nvE2jWQzoebxwfHKQ8PSmtTHGa3pGKXoBOlImgrmFMdsIvECD4JTbC8JiOv3hz/3tNsbZOMFWbyej8d6POK7Pmt9v+ubybvvXnsti2EYhgAAANq5Tm09AAAAgKuBUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEzBu60HcD1duHBBJ0+eVNeuXWWxWNp6OAAAoAUMw9CpU6cUFhamTp2aXo/pUKHm5MmTCg8Pb+thAACAVjh27JhuueWWJvd3qFDTtWtXSRcnJSAgoI1HAwAAWsJutys8PNzxe7wpHSrUXPrIKSAggFADAEA7c7lbR7hRGAAAmAKhBgAAmAKhBgAAmEKHuqcGAND+NTY2qqGhoa2HgavIy8tL3t7eV/y4FUINAKDdOH36tI4fPy7DMNp6KLjKbrrpJoWGhsrX17fV5yDUAADahcbGRh0/flw33XSTevXqxUNUTcIwDNXX1+urr75SWVmZbr311mYfsNccQg0AoF1oaGiQYRjq1auX/P3923o4uIr8/f3l4+Ojo0ePqr6+XlartVXn4UZhAEC7wgqNObV2dcbpHFdhHAAAAG2OUAMAAEyBe2oAAO1aZqa5+7ukT58+Sk9PV3p6+mVrX331VaWnp+ubb7655uO6kbBSAwDANTRhwgQ99NBDTm1/+ctfZLVatXDhwhaf56OPPtKvf/3rqzy6lpkwYYIsFossFot8fHwUFRWl6dOn68yZM5KkgoICWSyWNg9RrNQAAHAdvfzyy5oyZYqWL1+uSZMmtfi4Xr16XcNRXd7999+vV155RQ0NDdq5c6cmTZqkM2fOaOXKlW06ru9ipQYAgOtk4cKFeuaZZ/TGG284BZpLqzmLFi1SaGioevbsqSlTpjg9OblPnz7Kzs52/PzNN9/o17/+tYKDg2W1WjVw4EBt2bLFbb81NTX68Y9/rNGjR+vcuXP6+uuv9fjjjzu+Hn/rrbfqlVdeaXbsfn5+CgkJUXh4uB577DE9/vjjeuutt/Tll1/qnnvukSR1795dFotFEyZMkHRxRSomJkb+/v7q2bOnEhMTHas71wIrNUAHk9lWNwR4oD2MEfDUzJkztXz5cm3ZskWJiYku+99//32Fhobq/fff1+eff67U1FQNHjxYkydPdqm9cOGCkpOTderUKb3++uvq27evPv30U3l5ebnUHj9+XElJSYqNjdWaNWvk7e2t6dOn69NPP9W7776roKAgff755/r22289uh5/f381NDQoPDxcmzZtUkpKig4fPqyAgAD5+/uroqJCY8eO1cKFC/Xwww/r1KlT2rlz5zV9GjShBoBbBQVt13drMw1ZCDeqd999V3/961+1Y8cO3XvvvW5runfvrmXLlsnLy0v9+/fXAw88oB07drgNNdu3b9fevXtVWlqqH/7wh5KkqKgol7rPPvtMI0aM0M9//nP913/9l+MZP+Xl5brzzjsVGxsr6eIqkCf27t2rN954Q/fdd5+8vLzUo0cPSdIPfvADdevWTZL0z3/+U+fPn9cvfvELRURESJJiYmI86sdTfPwEAMA1dscdd6hPnz76wx/+oFOnTrmtGTBggNNKS2hoqKqqqtzWlpSU6JZbbnEEGne+/fZb/fSnP9VDDz2kJUuWOD208N/+7d+0YcMGDR48WL/73e+0e/fuy17Dli1b1KVLF1mtVsXFxWno0KFaunRpk/WDBg3Sfffdp5iYGP3yl7/U6tWr9fXXX1+2nytBqAEA4Bq7+eabVVhYqIqKCt1///1ug42Pj4/TzxaLRRcuXHB7vpa8JsLPz0+JiYl65513dPz4cad9ycnJOnr0qNLT03Xy5Endd999mj59erPnu+eee1RSUqLDhw/r3Llz2rx5s37wgx80We/l5aX8/Hy9++67uv3227V06VLddtttKisru+zYW4tQAwDAddC7d28VFhaqqqpKSUlJstvtrT7XHXfcoePHj+uzzz5rsqZTp0567bXXZLPZdO+99+rkyZNO+3v16qUJEybo9ddfV3Z2tlatWtVsn507d1a/fv0UERHhEsAuvVm7sbHRqd1isegnP/mJnn/+eRUXF8vX11dvvvmmJ5fqEUINAADXyS233KKCggLV1NQoKSlJtbW1rTrPsGHDNHToUKWkpCg/P19lZWV699139d577znVeXl56c9//rMGDRqke++9V5WVlZKkP/zhD/rrX/+qzz//XJ988om2bNmi6OjoVl9XRESELBaLtmzZoq+++kqnT5/Wnj179OKLL2rfvn0qLy/X5s2b9dVXX11RP5fDjcIAgHatvd0gfumjqHvuuUcjRozQtm3bWnWeTZs2afr06Ro7dqzOnDmjfv36af78+S513t7eWr9+vVJTU3XvvfeqoKBAvr6+mjVrlr788kv5+/srISFBGzZsuKJrev755zVz5kw9+eSTeuKJJzRjxgx98MEHys7Olt1uV0REhP7zP/9TycnJre7ncizGtfxu1Q3GbrcrMDBQtbW1CggIaOvhAG2ipV+XbstvPw0fntmq49rbLzd45ty5cyorK1NkZKSsVmtbDwdXWXP/flv6+7tVHz+tWLHC0anNZtPOnTubrS8sLJTNZpPValVUVJRycnKarN2wYYMsFovLI6Vb0y8AAOg4PA41GzduVHp6uubMmaPi4mIlJCQoOTlZ5eXlbuvLyso0atQoJSQkqLi4WLNnz9bUqVO1adMml9qjR49q+vTpSkhIuOJ+AQBAx+JxqFm8eLEmTpyoSZMmKTo6WtnZ2QoPD2/y3Q85OTnq3bu3srOzFR0drUmTJumpp57SokWLnOoaGxv1+OOP6/nnn3f7ACFP+wUAAB2LR6Gmvr5e+/fvV1JSklN7UlJSkw/uKSoqcqkfOXKk9u3b5/ROi3nz5qlXr16aOHHiVelXkurq6mS32502AABgTh6FmurqajU2Nio4ONipPTg42PE1se+rrKx0W3/+/HlVV1dLkv7+978rNzdXq1evvmr9SlJWVpYCAwMdW3h4+GWvEQAAtE+tulH4u49aliTDMFzaLld/qf3UqVP61a9+pdWrVysoKOiq9jtr1izV1tY6tmPHjjV7fgAA0H559JyaoKAgeXl5uayOVFVVuayiXBISEuK23tvbWz179tQnn3yiL7/8Ug8++KBj/6XHQnt7e+vw4cMKDw/3uF/p4iOi/fz8PLlEAADQTnm0UuPr6yubzab8/Hyn9vz8fMXHx7s9Ji4uzqV+27Ztio2NlY+Pj/r3769Dhw6ppKTEsY0ePdrxjonw8PBW9QsAADoWj58onJGRoXHjxik2NlZxcXFatWqVysvLlZaWJuniRz4nTpzQunXrJElpaWlatmyZMjIyNHnyZBUVFSk3N1fr16+XJFmtVg0cONCpj0uvLf9u++X6BQAAHZvHoSY1NVU1NTWaN2+eKioqNHDgQOXl5SkiIkKSVFFR4fTsmMjISOXl5WnatGlavny5wsLCtGTJEqWkpFzVfgEAHVNLn5LdVv1NmDBBa9euVVZWlmbOnOlof+utt/Twww/Lkwf79+nTR+np6UpPT/doDN81fPhwFRYWSrr4CUxERIQmTJigGTNmyMvLS6+++qrS09P1zTfftLqPttKqdz89/fTTevrpp93ue/XVV13ahg0bpgMHDrT4/O7Ocbl+gRvB9f7LFUD7YLVatWDBAv3mN79R9+7d23o4mjx5subNm6dz585py5Ytmjp1qry8vDRjxoy2HtoV4S3dAABcY4mJiQoJCVFWVlazdbt379bQoUPl7++v8PBwTZ06VWfOnJF0cYXl6NGjmjZtmiwWi+Pbv0ePHtWDDz6o7t27q3PnzhowYIDy8vKa7eemm25SSEiI+vTpo2eeeUb33Xef3nrrLRUUFOjJJ59UbW2to49L/7O2YsUK3XrrrbJarQoODtYjjzxy5RNzlRFqAAC4xry8vPTiiy9q6dKlOn78uNuaQ4cOaeTIkfrFL36hjz/+WBs3btSuXbv0zDPPSJI2b96sW265xXEbRkVFhSRpypQpqqur0wcffKBDhw5pwYIF6tKli0fj8/f3V0NDg+Lj45Wdna2AgABHH9OnT9e+ffs0depUzZs3T4cPH9Z7772noUOHXtmkXAOt+vgJAAB45uGHH9bgwYM1d+5c5ebmuux/6aWX9Nhjjznul7n11lu1ZMkSDRs2TCtXrlSPHj3k5eWlrl27KiQkxHFceXm5UlJSFBMTI0luXzXUlAsXLmjbtm3aunWr0tPT5evrq8DAQFksFpc+OnfurJ/97Gfq2rWrIiIidOedd7ZyJq4dVmoAALhOFixYoLVr1+rTTz912bd//369+uqr6tKli2MbOXKkLly4oLKysibPOXXqVP3xj3/UT37yE82dO1cff/zxZcexYsUKdenSRVarVaNHj9avfvUrzZ07t8n6ESNGKCIiQlFRURo3bpz+/Oc/6+zZsy276OuIUAMAwHUydOhQjRw5UrNnz3bZd+HCBf3mN79xem7bwYMHdeTIEfXt27fJc06aNElffPGFxo0bp0OHDik2NlZLly5tdhyPP/64SkpK9M9//lPffvutcnNzddNNNzVZ37VrVx04cEDr169XaGio/vCHP2jQoEE33DekCDUAAFxH8+fP1//8z/+4vJD5rrvu0ieffKJ+/fq5bL6+vpIufgW7sbHR5Zzh4eFKS0vT5s2b9e///u9NvkvxksDAQPXr18/xxP7vaqoPb29vJSYmauHChfr444/15Zdf6m9/+5unl39NEWoAALiOYmJi9Pjjj7uspsyYMUNFRUWaMmWKSkpKdOTIEb399tv67W9/66jp06ePPvjgA504ccLxUuj09HRt3bpVZWVlOnDggP72t78pOjq61ePr06ePTp8+rR07dqi6ulpnz57Vli1btGTJEpWUlOjo0aNat26dLly4oNtuu63V/VwLhBoAAK6zF154weWhe3fccYcKCwt15MgRJSQk6M4779Rzzz2n0NBQR828efP05Zdfqm/fvurVq5ckqbGxUVOmTFF0dLTuv/9+3XbbbVqxYkWrxxYfH6+0tDSlpqaqV69eWrhwobp166bNmzfr3nvvVXR0tHJycrR+/XoNGDCg1f1cCxbDk0cZtnN2u12BgYGqra1VQEBAWw8HJmSmh+8VFLRd38OHZ7bqOBNNP9w4d+6cysrKFBkZKavV2tbDwVXW3L/flv7+ZqUGAACYAqEGAACYAqEGAACYAqEGAACYAqEGAACYAqEGANCudKAv7XYoV+PfK6EGANAuXHrybX19fRuPBNfCpXdJ+fj4tPocvKUbANAueHt766abbtJXX30lHx8fderE/5ebgWEYOnv2rKqqqtStWzeX1zZ4glADAGgXLBaLQkNDVVZWpqNHj7b1cHCVdevWTSEhIVd0DkINAKDd8PX11a233spHUCbj4+NzRSs0lxBqAADtSqdOnXhNAtziA0kAAGAKhBoAAGAKhBoAAGAKhBoAAGAKhBoAAGAKhBoAAGAKhBoAAGAKhBoAAGAKhBoAAGAKhBoAAGAKhBoAAGAKrQo1K1asUGRkpKxWq2w2m3bu3NlsfWFhoWw2m6xWq6KiopSTk+O0f/PmzYqNjVW3bt3UuXNnDR48WK+99ppTTWZmpiwWi9N2pW/zBAAA5uFxqNm4caPS09M1Z84cFRcXKyEhQcnJySovL3dbX1ZWplGjRikhIUHFxcWaPXu2pk6dqk2bNjlqevTooTlz5qioqEgff/yxnnzyST355JPaunWr07kGDBigiooKx3bo0CFPhw8AAEzK47d0L168WBMnTtSkSZMkSdnZ2dq6datWrlyprKwsl/qcnBz17t1b2dnZkqTo6Gjt27dPixYtUkpKiiRp+PDhTsc8++yzWrt2rXbt2qWRI0f+a7De3qzOAAAAtzxaqamvr9f+/fuVlJTk1J6UlKTdu3e7PaaoqMilfuTIkdq3b58aGhpc6g3D0I4dO3T48GENHTrUad+RI0cUFhamyMhIjRkzRl988UWz462rq5PdbnfaAACAOXkUaqqrq9XY2Kjg4GCn9uDgYFVWVro9prKy0m39+fPnVV1d7Wirra1Vly5d5OvrqwceeEBLly7ViBEjHPuHDBmidevWaevWrVq9erUqKysVHx+vmpqaJseblZWlwMBAxxYeHu7J5QIAgHakVTcKWywWp58Nw3Bpu1z999u7du2qkpISffTRR/rTn/6kjIwMFRQUOPYnJycrJSVFMTExSkxM1DvvvCNJWrt2bZP9zpo1S7W1tY7t2LFjLb5GAADQvnh0T01QUJC8vLxcVmWqqqpcVmMuCQkJcVvv7e2tnj17Oto6deqkfv36SZIGDx6s0tJSZWVludxvc0nnzp0VExOjI0eONDlePz8/+fn5teTSAABAO+fRSo2vr69sNpvy8/Od2vPz8xUfH+/2mLi4OJf6bdu2KTY2Vj4+Pk32ZRiG6urqmtxfV1en0tJShYaGenAFAADArDz+9lNGRobGjRun2NhYxcXFadWqVSovL1daWpqkix/5nDhxQuvWrZMkpaWladmyZcrIyNDkyZNVVFSk3NxcrV+/3nHOrKwsxcbGqm/fvqqvr1deXp7WrVunlStXOmqmT5+uBx98UL1791ZVVZX++Mc/ym63a/z48Vc6BwAAwAQ8DjWpqamqqanRvHnzVFFRoYEDByovL08RERGSpIqKCqdn1kRGRiovL0/Tpk3T8uXLFRYWpiVLlji+zi1JZ86c0dNPP63jx4/L399f/fv31+uvv67U1FRHzfHjxzV27FhVV1erV69euvvuu/Xhhx86+gUAAB2bxbh0124HYLfbFRgYqNraWgUEBLT1cGBCmZmZbT2Eq+Y79+lfd8OHZ7bqOBNNP4DvaOnvb979BAAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATIFQAwAATKFVoWbFihWKjIyU1WqVzWbTzp07m60vLCyUzWaT1WpVVFSUcnJynPZv3rxZsbGx6tatmzp37qzBgwfrtddeu+J+AQBAx+FxqNm4caPS09M1Z84cFRcXKyEhQcnJySovL3dbX1ZWplGjRikhIUHFxcWaPXu2pk6dqk2bNjlqevTooTlz5qioqEgff/yxnnzyST355JPaunVrq/sFAAAdi8UwDMOTA4YMGaK77rpLK1eudLRFR0froYceUlZWlkv9jBkz9Pbbb6u0tNTRlpaWpoMHD6qoqKjJfu666y498MADeuGFF1rVrzt2u12BgYGqra1VQEBAi44BPJGZmdnWQ7hqCgraru/hwzNbdZyJph/Ad7T097dHKzX19fXav3+/kpKSnNqTkpK0e/dut8cUFRW51I8cOVL79u1TQ0ODS71hGNqxY4cOHz6soUOHtrpfSaqrq5PdbnfaAACAOXkUaqqrq9XY2Kjg4GCn9uDgYFVWVro9prKy0m39+fPnVV1d7Wirra1Vly5d5OvrqwceeEBLly7ViBEjWt2vJGVlZSkwMNCxhYeHe3K5AACgHWnVjcIWi8XpZ8MwXNouV//99q5du6qkpEQfffSR/vSnPykjI0MF31v/9rTfWbNmqba21rEdO3as2esCAADtl7cnxUFBQfLy8nJZHamqqnJZRbkkJCTEbb23t7d69uzpaOvUqZP69esnSRo8eLBKS0uVlZWl4cOHt6pfSfLz85Ofn58nlwgAANopj1ZqfH19ZbPZlJ+f79Sen5+v+Ph4t8fExcW51G/btk2xsbHy8fFpsi/DMFRXV9fqfgEAQMfi0UqNJGVkZGjcuHGKjY1VXFycVq1apfLycqWlpUm6+JHPiRMntG7dOkkXv+m0bNkyZWRkaPLkySoqKlJubq7Wr1/vOGdWVpZiY2PVt29f1dfXKy8vT+vWrXP6ptPl+gUAAB2bx6EmNTVVNTU1mjdvnioqKjRw4EDl5eUpIiJCklRRUeH07JjIyEjl5eVp2rRpWr58ucLCwrRkyRKlpKQ4as6cOaOnn35ax48fl7+/v/r376/XX39dqampLe4XAAB0bB4/p6Y94zk1uNZ4Ts3VwXNqAHzXNXlODQAAwI2KUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEyBUAMAAEzB49ckAPBMWz6ZFwA6ElZqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKbQq1KxYsUKRkZGyWq2y2WzauXNns/WFhYWy2WyyWq2KiopSTk6O0/7Vq1crISFB3bt3V/fu3ZWYmKi9e/c61WRmZspisThtISEhrRk+AAAwIY9DzcaNG5Wenq45c+aouLhYCQkJSk5OVnl5udv6srIyjRo1SgkJCSouLtbs2bM1depUbdq0yVFTUFCgsWPH6v3331dRUZF69+6tpKQknThxwulcAwYMUEVFhWM7dOiQp8MHAAAm5e3pAYsXL9bEiRM1adIkSVJ2dra2bt2qlStXKisry6U+JydHvXv3VnZ2tiQpOjpa+/bt06JFi5SSkiJJ+vOf/+x0zOrVq/WXv/xFO3bs0BNPPPGvwXp7szoDAADc8milpr6+Xvv371dSUpJTe1JSknbv3u32mKKiIpf6kSNHat++fWpoaHB7zNmzZ9XQ0KAePXo4tR85ckRhYWGKjIzUmDFj9MUXXzQ73rq6OtntdqcNAACYk0ehprq6Wo2NjQoODnZqDw4OVmVlpdtjKisr3dafP39e1dXVbo+ZOXOmbr75ZiUmJjrahgwZonXr1mnr1q1avXq1KisrFR8fr5qamibHm5WVpcDAQMcWHh7e0ksFAADtTKtuFLZYLE4/G4bh0na5enftkrRw4UKtX79emzdvltVqdbQnJycrJSVFMTExSkxM1DvvvCNJWrt2bZP9zpo1S7W1tY7t2LFjl784AADQLnl0T01QUJC8vLxcVmWqqqpcVmMuCQkJcVvv7e2tnj17OrUvWrRIL774orZv36477rij2bF07txZMTExOnLkSJM1fn5+8vPza/Y8aD8yMzPbeggAgBuYRys1vr6+stlsys/Pd2rPz89XfHy822Pi4uJc6rdt26bY2Fj5+Pg42l566SW98MILeu+99xQbG3vZsdTV1am0tFShoaGeXAIAADApjz9+ysjI0Msvv6w1a9aotLRU06ZNU3l5udLS0iRd/Mjnu99YSktL09GjR5WRkaHS0lKtWbNGubm5mj59uqNm4cKF+v3vf681a9aoT58+qqysVGVlpU6fPu2omT59ugoLC1VWVqY9e/bokUcekd1u1/jx46/k+gEAgEl4/JXu1NRU1dTUaN68eaqoqNDAgQOVl5eniIgISVJFRYXTM2siIyOVl5enadOmafny5QoLC9OSJUscX+eWLj7Mr76+Xo888ohTX3PnznV85HD8+HGNHTtW1dXV6tWrl+6++259+OGHjn4BAEDHZjEu3bXbAdjtdgUGBqq2tlYBAQFtPRx4qL3eU1NQ0NYjaH+GD89s1XHt8Y8IYwYur6W/v3n3EwAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAXvth4AAHxfQUFmq467nm+Pbq9vjQfMjJUaAABgCqzUADCNgoLr1xcLNcCNh5UaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCq0KNStWrFBkZKSsVqtsNpt27tzZbH1hYaFsNpusVquioqKUk5PjtH/16tVKSEhQ9+7d1b17dyUmJmrv3r1X3C8AAOg4PA41GzduVHp6uubMmaPi4mIlJCQoOTlZ5eXlbuvLyso0atQoJSQkqLi4WLNnz9bUqVO1adMmR01BQYHGjh2r999/X0VFRerdu7eSkpJ04sSJVvcLAAA6FothGIYnBwwZMkR33XWXVq5c6WiLjo7WQw89pKysLJf6GTNm6O2331ZpaamjLS0tTQcPHlRRUZHbPhobG9W9e3ctW7ZMTzzxRKv6dcdutyswMFC1tbUKCAho0TG4cWRmZrb1EFqloKCtR4BrYfjwzLYeQptpp/8poh1r6e9vj1Zq6uvrtX//fiUlJTm1JyUlaffu3W6PKSoqcqkfOXKk9u3bp4aGBrfHnD17Vg0NDerRo0er+5Wkuro62e12pw0AAJiTR6GmurpajY2NCg4OdmoPDg5WZWWl22MqKyvd1p8/f17V1dVuj5k5c6ZuvvlmJSYmtrpfScrKylJgYKBjCw8Pv+w1AgCA9qlVNwpbLBannw3DcGm7XL27dklauHCh1q9fr82bN8tqtV5Rv7NmzVJtba1jO3bsWJO1AACgffP2pDgoKEheXl4uqyNVVVUuqyiXhISEuK339vZWz549ndoXLVqkF198Udu3b9cdd9xxRf1Kkp+fn/z8/Fp0bQAAoH3zaKXG19dXNptN+fn5Tu35+fmKj493e0xcXJxL/bZt2xQbGysfHx9H20svvaQXXnhB7733nmJjY6+4XwAA0LF4tFIjSRkZGRo3bpxiY2MVFxenVatWqby8XGlpaZIufuRz4sQJrVu3TtLFbzotW7ZMGRkZmjx5soqKipSbm6v169c7zrlw4UI999xzeuONN9SnTx/HikyXLl3UpUuXFvULAAA6No9DTWpqqmpqajRv3jxVVFRo4MCBysvLU0REhCSpoqLC6dkxkZGRysvL07Rp07R8+XKFhYVpyZIlSklJcdSsWLFC9fX1euSRR5z6mjt3ruNrvJfrFwAAdGweP6emPeM5Ne0bz6nBjYTn1ADXzzV5Tg0AAMCNilADAABMgVADAABMgVADAABMweNvPwEAOrb2eKNwexwzPMdKDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMAVCDQAAMIVWhZoVK1YoMjJSVqtVNptNO3fubLa+sLBQNptNVqtVUVFRysnJcdr/ySefKCUlRX369JHFYlF2drbLOTIzM2WxWJy2kJCQ1gwfAACYkMehZuPGjUpPT9ecOXNUXFyshIQEJScnq7y83G19WVmZRo0apYSEBBUXF2v27NmaOnWqNm3a5Kg5e/asoqKiNH/+/GaDyoABA1RRUeHYDh065OnwAQCASXl7esDixYs1ceJETZo0SZKUnZ2trVu3auXKlcrKynKpz8nJUe/evR2rL9HR0dq3b58WLVqklJQUSdKPfvQj/ehHP5IkzZw5s+nBenuzOgMAANzyaKWmvr5e+/fvV1JSklN7UlKSdu/e7faYoqIil/qRI0dq3759amho8GiwR44cUVhYmCIjIzVmzBh98cUXHh0PAADMy6NQU11drcbGRgUHBzu1BwcHq7Ky0u0xlZWVbuvPnz+v6urqFvc9ZMgQrVu3Tlu3btXq1atVWVmp+Ph41dTUNHlMXV2d7Ha70wYAAMypVTcKWywWp58Nw3Bpu1y9u/bmJCcnKyUlRTExMUpMTNQ777wjSVq7dm2Tx2RlZSkwMNCxhYeHt7g/AADQvngUaoKCguTl5eWyKlNVVeWyGnNJSEiI23pvb2/17NnTw+H+S+fOnRUTE6MjR440WTNr1izV1tY6tmPHjrW6PwAAcGPzKNT4+vrKZrMpPz/fqT0/P1/x8fFuj4mLi3Op37Ztm2JjY+Xj4+PhcP+lrq5OpaWlCg0NbbLGz89PAQEBThsAADAnjz9+ysjI0Msvv6w1a9aotLRU06ZNU3l5udLS0iRdXB154oknHPVpaWk6evSoMjIyVFpaqjVr1ig3N1fTp0931NTX16ukpEQlJSWqr6/XiRMnVFJSos8//9xRM336dBUWFqqsrEx79uzRI488IrvdrvHjx1/J9QMAAJPw+Cvdqampqqmp0bx581RRUaGBAwcqLy9PERERkqSKigqnZ9ZERkYqLy9P06ZN0/LlyxUWFqYlS5Y4vs4tSSdPntSdd97p+HnRokVatGiRhg0bpoKCAknS8ePHNXbsWFVXV6tXr166++679eGHHzr6xZXJzMxs6yEA7UpBQWZbD+Gyhg/PbOshANeVxbh0124HYLfbFRgYqNraWj6K+h5CzbXz/3M5cN0Rav6Fv+Lat5b+/ubdTwAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBQINQAAwBS823oAAADAVWZmW4/Ac209ZlZqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKRBqAACAKbQq1KxYsUKRkZGyWq2y2WzauXNns/WFhYWy2WyyWq2KiopSTk6O0/5PPvlEKSkp6tOnjywWi7Kzs69KvwAAoOPwONRs3LhR6enpmjNnjoqLi5WQkKDk5GSVl5e7rS8rK9OoUaOUkJCg4uJizZ49W1OnTtWmTZscNWfPnlVUVJTmz5+vkJCQq9IvAADoWDwONYsXL9bEiRM1adIkRUdHKzs7W+Hh4Vq5cqXb+pycHPXu3VvZ2dmKjo7WpEmT9NRTT2nRokWOmh/96Ed66aWXNGbMGPn5+V2VfgEAQMfiUaipr6/X/v37lZSU5NSelJSk3bt3uz2mqKjIpX7kyJHat2+fGhoarlm/klRXVye73e60AQAAc/Io1FRXV6uxsVHBwcFO7cHBwaqsrHR7TGVlpdv68+fPq7q6+pr1K0lZWVkKDAx0bOHh4S3qDwAAtD+tulHYYrE4/WwYhkvb5erdtV/tfmfNmqXa2lrHduzYMY/6AwAA7YdHL7QMCgqSl5eXy+pIVVWVyyrKJSEhIW7rvb291bNnz2vWryT5+fk1eY8OAAAwF49Wanx9fWWz2ZSfn+/Unp+fr/j4eLfHxMXFudRv27ZNsbGx8vHxuWb9AgCAjsWjlRpJysjI0Lhx4xQbG6u4uDitWrVK5eXlSktLk3TxI58TJ05o3bp1kqS0tDQtW7ZMGRkZmjx5soqKipSbm6v169c7zllfX69PP/3U8c8nTpxQSUmJunTpon79+rWoXwAA0LF5HGpSU1NVU1OjefPmqaKiQgMHDlReXp4iIiIkSRUVFU7PjomMjFReXp6mTZum5cuXKywsTEuWLFFKSoqj5uTJk7rzzjsdPy9atEiLFi3SsGHDVFBQ0KJ+AQBAx2YxLt212wHY7XYFBgaqtrZWAQEBbT2cG0pmZmZbD8G0/n8uB6674cMz23oIN4z2+FccY/6Xlv7+5t1PAADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFAg1AADAFFoValasWKHIyEhZrVbZbDbt3Lmz2frCwkLZbDZZrVZFRUUpJyfHpWbTpk26/fbb5efnp9tvv11vvvmm0/7MzExZLBanLSQkpDXDBwAAJuRxqNm4caPS09M1Z84cFRcXKyEhQcnJySovL3dbX1ZWplGjRikhIUHFxcWaPXu2pk6dqk2bNjlqioqKlJqaqnHjxungwYMaN26cHn30Ue3Zs8fpXAMGDFBFRYVjO3TokKfDBwAAJuVxqFm8eLEmTpyoSZMmKTo6WtnZ2QoPD9fKlSvd1ufk5Kh3797Kzs5WdHS0Jk2apKeeekqLFi1y1GRnZ2vEiBGaNWuW+vfvr1mzZum+++5Tdna207m8vb0VEhLi2Hr16uXp8AEAgEl5FGrq6+u1f/9+JSUlObUnJSVp9+7dbo8pKipyqR85cqT27dunhoaGZmu+f84jR44oLCxMkZGRGjNmjL744otmx1tXVye73e60AQAAc/Io1FRXV6uxsVHBwcFO7cHBwaqsrHR7TGVlpdv68+fPq7q6utma755zyJAhWrdunbZu3arVq1ersrJS8fHxqqmpaXK8WVlZCgwMdGzh4eGeXC4AAGhHWnWjsMVicfrZMAyXtsvVf7/9cudMTk5WSkqKYmJilJiYqHfeeUeStHbt2ib7nTVrlmprax3bsWPHLnNlAACgvfL2pDgoKEheXl4uqzJVVVUuKy2XhISEuK339vZWz549m61p6pyS1LlzZ8XExOjIkSNN1vj5+cnPz6/ZawIAAObg0UqNr6+vbDab8vPzndrz8/MVHx/v9pi4uDiX+m3btik2NlY+Pj7N1jR1Tuni/TKlpaUKDQ315BIAAIBJefzxU0ZGhl5++WWtWbNGpaWlmjZtmsrLy5WWlibp4kc+TzzxhKM+LS1NR48eVUZGhkpLS7VmzRrl5uZq+vTpjppnn31W27Zt04IFC/S///u/WrBggbZv36709HRHzfTp01VYWKiysjLt2bNHjzzyiOx2u8aPH38Flw8AAMzCo4+fJCk1NVU1NTWaN2+eKioqNHDgQOXl5SkiIkKSVFFR4fTMmsjISOXl5WnatGlavny5wsLCtGTJEqWkpDhq4uPjtWHDBv3+97/Xc889p759+2rjxo0aMmSIo+b48eMaO3asqqur1atXL91999368MMPHf0CAICOzWJcumu3A7Db7QoMDFRtba0CAgLaejg3lMzMzLYegmkVFLT1CNBRDR+e2dZDuGG0x7/iGPO/tPT3t8crNfAcgQFAWygoyGzrIVwWwQtXEy+0BAAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApkCoAQAApuDd1gMAAHRcBQWZ16WfzOvTzVVVUNCyuuHDM6/lMNoVQs1V0tx/MC39gwkAuDb4e7hj4OMnAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCoQaAABgCq0KNStWrFBkZKSsVqtsNpt27tzZbH1hYaFsNpusVquioqKUk5PjUrNp0ybdfvvt8vPz0+23364333zzivsFAAAdh8ehZuPGjUpPT9ecOXNUXFyshIQEJScnq7y83G19WVmZRo0apYSEBBUXF2v27NmaOnWqNm3a5KgpKipSamqqxo0bp4MHD2rcuHF69NFHtWfPnlb3CwAAOhaLYRiGJwcMGTJEd911l1auXOloi46O1kMPPaSsrCyX+hkzZujtt99WaWmpoy0tLU0HDx5UUVGRJCk1NVV2u13vvvuuo+b+++9X9+7dtX79+lb1647dbldgYKBqa2sVEBDgyWVfVvNv6W5mJwAAV2D48My2HoJDc78Lr0RLf397e3LS+vp67d+/XzNnznRqT0pK0u7du90eU1RUpKSkJKe2kSNHKjc3Vw0NDfLx8VFRUZGmTZvmUpOdnd3qfiWprq5OdXV1jp9ra2slXZycq+073bg4f76ZnQAAXIG6uqv/O621rsGv1/9/3osnvtw6jEehprq6Wo2NjQoODnZqDw4OVmVlpdtjKisr3dafP39e1dXVCg0NbbLm0jlb068kZWVl6fnnn3dpDw8Pb/oiAQBoR/7+9/ltPQSH+dd4KKdOnVJgYGCT+z0KNZdYLBannw3DcGm7XP3321tyTk/7nTVrljIyMhw/X7hwQf/3f/+nnj17NnvcjcButys8PFzHjh276h+VtVfMiXvMiyvmxBVz4h7z4upGnBPDMHTq1CmFhYU1W+dRqAkKCpKXl5fL6khVVZXLKsolISEhbuu9vb3Vs2fPZmsunbM1/UqSn5+f/Pz8nNq6devW9AXegAICAm6YP1Q3CubEPebFFXPiijlxj3lxdaPNSXMrNJd49O0nX19f2Ww25efnO7Xn5+crPj7e7TFxcXEu9du2bVNsbKx8fHyarbl0ztb0CwAAOhaPP37KyMjQuHHjFBsbq7i4OK1atUrl5eVKS0uTdPEjnxMnTmjdunWSLn7TadmyZcrIyNDkyZNVVFSk3Nxcx7eaJOnZZ5/V0KFDtWDBAv385z/XX//6V23fvl27du1qcb8AAKCDM1ph+fLlRkREhOHr62vcddddRmFhoWPf+PHjjWHDhjnVFxQUGHfeeafh6+tr9OnTx1i5cqXLOf/7v//buO222wwfHx+jf//+xqZNmzzq12zOnTtnzJ071zh37lxbD+WGwZy4x7y4Yk5cMSfuMS+u2vOcePycGgAAgBsR734CAACmQKgBAACmQKgBAACmQKgBAACmQKhpQ3369JHFYnHZpkyZIuniExQzMzMVFhYmf39/DR8+XJ988kkbj/raOn/+vH7/+98rMjJS/v7+ioqK0rx583ThwgVHTUecF+ni48HT09MVEREhf39/xcfH66OPPnLsN/u8fPDBB3rwwQcVFhYmi8Wit956y2l/S66/rq5Ov/3tbxUUFKTOnTtr9OjROn78+HW8iqvvcvOyefNmjRw5UkFBQbJYLCopKXE5h9nmpbk5aWho0IwZMxQTE6POnTsrLCxMTzzxhE6ePOl0jo40J5KUmZmp/v37q3PnzurevbsSExO1Z88ep5r2MCeEmjb00UcfqaKiwrFderjgL3/5S0nSwoULtXjxYi1btkwfffSRQkJCNGLECJ06daoth31NLViwQDk5OVq2bJlKS0u1cOFCvfTSS1q6dKmjpiPOiyRNmjRJ+fn5eu2113To0CElJSUpMTFRJ06ckGT+eTlz5owGDRqkZcuWud3fkutPT0/Xm2++qQ0bNmjXrl06ffq0fvazn6mxsfF6XcZVd7l5OXPmjH7yk59ofjMv5THbvDQ3J2fPntWBAwf03HPP6cCBA9q8ebM+++wzjR492qmuI82JJP3whz/UsmXLdOjQIe3atUt9+vRRUlKSvvrqK0dNu5iTtvw+OZw9++yzRt++fY0LFy4YFy5cMEJCQoz58+c79p87d84IDAw0cnJy2nCU19YDDzxgPPXUU05tv/jFL4xf/epXhmEYHXZezp49a3h5eRlbtmxxah80aJAxZ86cDjcvkow333zT8XNLrv+bb74xfHx8jA0bNjhqTpw4YXTq1Ml47733rtvYr6Xvz8t3lZWVGZKM4uJip3azz0tzc3LJ3r17DUnG0aNHDcNgTgzDMGpraw1Jxvbt2w3DaD9zwkrNDaK+vl6vv/66nnrqKVksFpWVlamyslJJSUmOGj8/Pw0bNky7d+9uw5FeWz/96U+1Y8cOffbZZ5KkgwcPateuXRo1apQkddh5OX/+vBobG2W1Wp3a/f39tWvXrg47L5e05Pr379+vhoYGp5qwsDANHDiwQ8xRU5gXqba2VhaLxfFuwI4+J/X19Vq1apUCAwM1aNAgSe1nTlr1lm5cfW+99Za++eYbTZgwQZIcL+/8/gs7g4ODdfTo0es9vOtmxowZqq2tVf/+/eXl5aXGxkb96U9/0tixYyV13Hnp2rWr4uLi9MILLyg6OlrBwcFav3699uzZo1tvvbXDzsslLbn+yspK+fr6qnv37i41339ZbkfS0efl3Llzmjlzph577DHHyxs76pxs2bJFY8aM0dmzZxUaGqr8/HwFBQVJaj9zwkrNDSI3N1fJyckur1W3WCxOPxuG4dJmJhs3btTrr7+uN954QwcOHNDatWu1aNEirV271qmuo82LJL322msyDEM333yz/Pz8tGTJEj322GPy8vJy1HTEefmu1lx/R5ujluoI89LQ0KAxY8bowoULWrFixWXrzT4n99xzj0pKSrR7927df//9evTRR1VVVdXsMTfanBBqbgBHjx7V9u3bNWnSJEdbSEiIJLkk4KqqKpf/GzWT//iP/9DMmTM1ZswYxcTEaNy4cZo2bZqysrIkddx5kaS+ffuqsLBQp0+f1rFjx7R37141NDQoMjKyQ8+L1LI/FyEhIaqvr9fXX3/dZE1H1FHnpaGhQY8++qjKysqUn5/vWKWROu6cdO7cWf369dPdd9+t3NxceXt7Kzc3V1L7mRNCzQ3glVde0Q9+8AM98MADjrZLv6gufSNKuvg5Z2FhoeLj49timNfF2bNn1amT8x9LLy8vx1e6O+q8fFfnzp0VGhqqr7/+Wlu3btXPf/7zDj8vLbl+m80mHx8fp5qKigr94x//6BBz1JSOOC+XAs2RI0e0fft29ezZ02l/R5wTdwzDUF1dnaR2NCdtdYcyLmpsbDR69+5tzJgxw2Xf/PnzjcDAQGPz5s3GoUOHjLFjxxqhoaGG3W5vg5FeH+PHjzduvvlmY8uWLUZZWZmxefNmIygoyPjd737nqOmI82IYhvHee+8Z7777rvHFF18Y27ZtMwYNGmT8+Mc/Nurr6w3DMP+8nDp1yiguLjaKi4sNScbixYuN4uJixzdWWnL9aWlpxi233GJs377dOHDggHHvvfcagwYNMs6fP99Wl3XFLjcvNTU1RnFxsfHOO+8YkowNGzYYxcXFRkVFheMcZpuX5uakoaHBGD16tHHLLbcYJSUlRkVFhWOrq6tznKMjzcnp06eNWbNmGUVFRcaXX35p7N+/35g4caLh5+dn/OMf/3Ccoz3MCaGmjW3dutWQZBw+fNhl34ULF4y5c+caISEhhp+fnzF06FDj0KFDbTDK68dutxvPPvus0bt3b8NqtRpRUVHGnDlznP6y6YjzYhiGsXHjRiMqKsrw9fU1QkJCjClTphjffPONY7/Z5+X99983JLls48ePNwyjZdf/7bffGs8884zRo0cPw9/f3/jZz35mlJeXt8HVXD2Xm5dXXnnF7f65c+c6zmG2eWluTi59td3d9v777zvO0ZHm5NtvvzUefvhhIywszPD19TVCQ0ON0aNHG3v37nU6R3uYE4thGMY1XQoCAAC4DrinBgAAmAKhBgAAmAKhBgAAmAKhBgAAmAKhBgAAmAKhBgAAmAKhBgAAmAKhBgAAmAKhBgAAmAKhBgAAmAKhBgAAmAKhBgAAmML/A6kDLGrFU3oKAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(knicks_pts, color='blue', label='Knicks Pts', density=True, alpha=0.5)\n",
    "plt.hist(nets_pts, color='black', label='Nets Pts', density=True, alpha=0.5)\n",
    "plt.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 4\n",
    "Now, let’s compare the 2010 games to 2014. Replicate the steps from Tasks 2 and 3 using `nba_2014`. First, calculate the mean difference between the two teams points scored. Save and print the value as `diff_means_2014`. Did the difference in points get larger or smaller in 2014? Then, plot the overlapping histograms. Does the mean difference you calculated make sense?\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 5\n",
    "For the remainder of this project, we’ll focus on data from 2010. Let’s now include all teams in the dataset and investigate the relationship between franchise and points scored per game.\n",
    "\n",
    "Using nba_2010, generate side-by-side boxplots with points scored (pts) on the y-axis and team (fran_id) on the x-axis. Is there any overlap between the boxes? Does this chart suggest that fran_id and pts are associated? Which pairs of teams, if any, earn different average scores per game?\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 6\n",
    "We'd like to know if teams tend to win more games at home compared to away.\n",
    "\n",
    "The variable, `game_result`, indicates whether a team won a particular game ('W' stands for “win” and 'L' stands for “loss”). The variable, `game_location`, indicates whether a team was playing at home or away ('H' stands for “home” and 'A' stands for “away”). \n",
    "\n",
    "Data scientists will often calculate a contingency table of frequencies to help them determine if categorical variables are associated. Calculate a table of frequencies that shows the counts of game_result and game_location.\n",
    "\n",
    "Save your result as `location_result_freq` and print your result. Based on this table, do you think the variables are associated?`\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 7\n",
    "\n",
    "Convert this table of frequencies to a table of proportions and save the result as `location_result_proportions`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 8\n",
    "Using the contingency table created above (Task 6), calculate the expected contingency table (if there were no association) and the Chi-Square statistic.\n",
    "\n",
    "Does the actual contingency table look similar to the expected table — or different? Based on this output, do you think there is an association between these variables?\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*For a 2x2 table, Chi-squared greater than about 4 indicates an association. We're not there*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 9\n",
    "\n",
    "For each game, 538 has calculated the probability that each team will win the game. We want to know if teams with a higher probability of winning (according to 538) also tend to win games by more points. \n",
    "\n",
    "In the data, 538's prediction is saved as `forecast`. The `point_diff` column gives the margin of victory/defeat for each team (positive values mean that the team won; negative values mean that they lost). \n",
    "\n",
    "Using `nba_2010`, calculate the covariance between `forecast` (538's projected win probability) and `point_diff` (the margin of victory/defeat) in the dataset. Save and print your result. Looking at the matrix, what is the covariance between these two variables?\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 10\n",
    "\n",
    "Because 538’s forecast variable is reported as a probability (not a binary), we can calculate the strength of the correlation.\n",
    "\n",
    "Using nba_2010, calculate the correlation between `forecast` and `point_diff`. Call this `point_diff_forecast_corr`. Save and print your result. Does this value suggest an association between the two variables?\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Task 11\n",
    "\n",
    "Generate a scatter plot of `forecast` (on the x-axis) and `point_diff` (on the y-axis). Does the correlation value make sense?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3.9.13 ('myenv')",
   "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.9.13"
  },
  "vscode": {
   "interpreter": {
    "hash": "c3d8a8be75f9e7e90c92f56a26dc284ac664b6f4edfe25f866cebb8c85a664d9"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}