main.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 193,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "\n",
    "weather = pd.read_csv(\"weather.csv\",index_col=\"DATE\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 194,
   "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>STATION</th>\n",
       "      <th>NAME</th>\n",
       "      <th>ACMH</th>\n",
       "      <th>ACSH</th>\n",
       "      <th>AWND</th>\n",
       "      <th>FMTM</th>\n",
       "      <th>PGTM</th>\n",
       "      <th>PRCP</th>\n",
       "      <th>SNOW</th>\n",
       "      <th>SNWD</th>\n",
       "      <th>...</th>\n",
       "      <th>WT13</th>\n",
       "      <th>WT14</th>\n",
       "      <th>WT15</th>\n",
       "      <th>WT16</th>\n",
       "      <th>WT17</th>\n",
       "      <th>WT18</th>\n",
       "      <th>WT19</th>\n",
       "      <th>WT21</th>\n",
       "      <th>WT22</th>\n",
       "      <th>WV01</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1980-04-01</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>20.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>24.0</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-02</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>70.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2242.0</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-03</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>30.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1654.0</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-04</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>80.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>624.0</td>\n",
       "      <td>0.57</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-05</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>40.0</td>\n",
       "      <td>50.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1548.0</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-20</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>12.08</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.67</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-21</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>13.65</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.12</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-22</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>19.69</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-23</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>12.08</td>\n",
       "      <td>NaN</td>\n",
       "      <td>39.0</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-24</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>15912 rows × 45 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                STATION                              NAME  ACMH   ACSH   AWND  \\\n",
       "DATE                                                                            \n",
       "1980-04-01  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  20.0   10.0    NaN   \n",
       "1980-04-02  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  70.0  100.0    NaN   \n",
       "1980-04-03  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  30.0   20.0    NaN   \n",
       "1980-04-04  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  80.0  100.0    NaN   \n",
       "1980-04-05  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  40.0   50.0    NaN   \n",
       "...                 ...                               ...   ...    ...    ...   \n",
       "2023-10-20  USW00094789  JFK INTERNATIONAL AIRPORT, NY US   NaN    NaN  12.08   \n",
       "2023-10-21  USW00094789  JFK INTERNATIONAL AIRPORT, NY US   NaN    NaN  13.65   \n",
       "2023-10-22  USW00094789  JFK INTERNATIONAL AIRPORT, NY US   NaN    NaN  19.69   \n",
       "2023-10-23  USW00094789  JFK INTERNATIONAL AIRPORT, NY US   NaN    NaN  12.08   \n",
       "2023-10-24  USW00094789  JFK INTERNATIONAL AIRPORT, NY US   NaN    NaN    NaN   \n",
       "\n",
       "            FMTM    PGTM  PRCP  SNOW  SNWD  ...  WT13  WT14  WT15  WT16  WT17  \\\n",
       "DATE                                        ...                                 \n",
       "1980-04-01   NaN    24.0  0.04   0.0   0.0  ...   NaN   NaN   NaN   1.0   NaN   \n",
       "1980-04-02   NaN  2242.0  0.00   0.0   0.0  ...   NaN   NaN   NaN   NaN   NaN   \n",
       "1980-04-03   NaN  1654.0  0.00   0.0   0.0  ...   NaN   NaN   NaN   NaN   NaN   \n",
       "1980-04-04   NaN   624.0  0.57   0.0   0.0  ...   NaN   NaN   NaN   1.0   NaN   \n",
       "1980-04-05   NaN  1548.0  0.00   0.0   0.0  ...   NaN   NaN   NaN   NaN   NaN   \n",
       "...          ...     ...   ...   ...   ...  ...   ...   ...   ...   ...   ...   \n",
       "2023-10-20   NaN     NaN  0.67   0.0   0.0  ...   NaN   NaN   NaN   NaN   NaN   \n",
       "2023-10-21   NaN     NaN  0.12   0.0   0.0  ...   NaN   NaN   NaN   NaN   NaN   \n",
       "2023-10-22   NaN     NaN  0.00   0.0   0.0  ...   NaN   NaN   NaN   NaN   NaN   \n",
       "2023-10-23   NaN    39.0  0.00   0.0   0.0  ...   NaN   NaN   NaN   NaN   NaN   \n",
       "2023-10-24   NaN     NaN   NaN   NaN   NaN  ...   NaN   NaN   NaN   NaN   NaN   \n",
       "\n",
       "            WT18  WT19  WT21  WT22  WV01  \n",
       "DATE                                      \n",
       "1980-04-01   1.0   NaN   NaN   NaN   NaN  \n",
       "1980-04-02   NaN   NaN   NaN   NaN   NaN  \n",
       "1980-04-03   NaN   NaN   NaN   NaN   NaN  \n",
       "1980-04-04   NaN   NaN   NaN   NaN   NaN  \n",
       "1980-04-05   NaN   NaN   NaN   NaN   NaN  \n",
       "...          ...   ...   ...   ...   ...  \n",
       "2023-10-20   NaN   NaN   NaN   NaN   NaN  \n",
       "2023-10-21   NaN   NaN   NaN   NaN   NaN  \n",
       "2023-10-22   NaN   NaN   NaN   NaN   NaN  \n",
       "2023-10-23   NaN   NaN   NaN   NaN   NaN  \n",
       "2023-10-24   NaN   NaN   NaN   NaN   NaN  \n",
       "\n",
       "[15912 rows x 45 columns]"
      ]
     },
     "execution_count": 194,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 195,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "STATION    0.000000\n",
       "NAME       0.000000\n",
       "ACMH       0.630970\n",
       "ACSH       0.630907\n",
       "AWND       0.086287\n",
       "FMTM       0.363751\n",
       "PGTM       0.333710\n",
       "PRCP       0.000063\n",
       "SNOW       0.000063\n",
       "SNWD       0.000126\n",
       "TAVG       0.589492\n",
       "TMAX       0.000063\n",
       "TMIN       0.000063\n",
       "TSUN       0.998177\n",
       "WDF1       0.631222\n",
       "WDF2       0.369218\n",
       "WDF5       0.374434\n",
       "WDFG       0.678167\n",
       "WDFM       0.999937\n",
       "WESD       0.624183\n",
       "WSF1       0.631033\n",
       "WSF2       0.369218\n",
       "WSF5       0.374497\n",
       "WSFG       0.635747\n",
       "WSFM       0.999937\n",
       "WT01       0.617584\n",
       "WT02       0.941868\n",
       "WT03       0.929927\n",
       "WT04       0.982843\n",
       "WT05       0.977438\n",
       "WT06       0.991139\n",
       "WT07       0.993213\n",
       "WT08       0.860671\n",
       "WT09       0.992836\n",
       "WT11       0.999120\n",
       "WT13       0.862682\n",
       "WT14       0.944256\n",
       "WT15       0.997360\n",
       "WT16       0.683446\n",
       "WT17       0.996229\n",
       "WT18       0.943942\n",
       "WT19       0.999937\n",
       "WT21       0.999686\n",
       "WT22       0.996921\n",
       "WV01       0.999937\n",
       "dtype: float64"
      ]
     },
     "execution_count": 195,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "null_pct = weather.apply(pd.isnull).sum()/weather.shape[0]\n",
    "null_pct"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 196,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['STATION', 'NAME', 'PRCP', 'SNOW', 'SNWD', 'TMAX', 'TMIN'], dtype='object')"
      ]
     },
     "execution_count": 196,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "vld_clms = weather.columns[null_pct < 0.05]\n",
    "vld_clms"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 197,
   "metadata": {},
   "outputs": [],
   "source": [
    "weather = weather[vld_clms].copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 198,
   "metadata": {},
   "outputs": [],
   "source": [
    "weather.columns = weather.columns.str.lower()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 199,
   "metadata": {},
   "outputs": [],
   "source": [
    "weather = weather.ffill()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 200,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "station    0\n",
       "name       0\n",
       "prcp       0\n",
       "snow       0\n",
       "snwd       0\n",
       "tmax       0\n",
       "tmin       0\n",
       "dtype: int64"
      ]
     },
     "execution_count": 200,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.apply(pd.isnull).sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 201,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "station     object\n",
       "name        object\n",
       "prcp       float64\n",
       "snow       float64\n",
       "snwd       float64\n",
       "tmax       float64\n",
       "tmin       float64\n",
       "dtype: object"
      ]
     },
     "execution_count": 201,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.dtypes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 202,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['1980-04-01', '1980-04-02', '1980-04-03', '1980-04-04', '1980-04-05',\n",
       "       '1980-04-06', '1980-04-07', '1980-04-08', '1980-04-09', '1980-04-10',\n",
       "       ...\n",
       "       '2023-10-15', '2023-10-16', '2023-10-17', '2023-10-18', '2023-10-19',\n",
       "       '2023-10-20', '2023-10-21', '2023-10-22', '2023-10-23', '2023-10-24'],\n",
       "      dtype='object', name='DATE', length=15912)"
      ]
     },
     "execution_count": 202,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.index"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 203,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "DatetimeIndex(['1980-04-01', '1980-04-02', '1980-04-03', '1980-04-04',\n",
       "               '1980-04-05', '1980-04-06', '1980-04-07', '1980-04-08',\n",
       "               '1980-04-09', '1980-04-10',\n",
       "               ...\n",
       "               '2023-10-15', '2023-10-16', '2023-10-17', '2023-10-18',\n",
       "               '2023-10-19', '2023-10-20', '2023-10-21', '2023-10-22',\n",
       "               '2023-10-23', '2023-10-24'],\n",
       "              dtype='datetime64[ns]', name='DATE', length=15912, freq=None)"
      ]
     },
     "execution_count": 203,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.index = pd.to_datetime(weather.index)\n",
    "weather.index"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 204,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "DATE\n",
       "1980    275\n",
       "1981    365\n",
       "1982    365\n",
       "1983    365\n",
       "1984    366\n",
       "1985    365\n",
       "1986    365\n",
       "1987    365\n",
       "1988    366\n",
       "1989    365\n",
       "1990    365\n",
       "1991    365\n",
       "1992    366\n",
       "1993    365\n",
       "1994    365\n",
       "1995    365\n",
       "1996    366\n",
       "1997    365\n",
       "1998    365\n",
       "1999    365\n",
       "2000    366\n",
       "2001    365\n",
       "2002    365\n",
       "2003    365\n",
       "2004    366\n",
       "2005    365\n",
       "2006    365\n",
       "2007    365\n",
       "2008    366\n",
       "2009    365\n",
       "2010    365\n",
       "2011    365\n",
       "2012    366\n",
       "2013    365\n",
       "2014    365\n",
       "2015    365\n",
       "2016    366\n",
       "2017    365\n",
       "2018    365\n",
       "2019    365\n",
       "2020    366\n",
       "2021    365\n",
       "2022    365\n",
       "2023    297\n",
       "Name: count, dtype: int64"
      ]
     },
     "execution_count": 204,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.index.year.value_counts().sort_index()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 205,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<Axes: xlabel='DATE'>"
      ]
     },
     "execution_count": 205,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGwCAYAAAAJ/wd3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABJtElEQVR4nO3de3wU1f0//teGQAAhwQAhUKJFbYWPKFZqkdb6wUq5VK0on4+WT1uxH1utBVtLq5ZvBdTWxktV9COi7U9Aa613sN6i3FEBFSRAECNEkFsS5JIEArntnt8fmGU3md2dM3vOzJnZ1/PxyOMBu7Mz58ycOfOeM2fOCQkhBIiIiIhckuV1AoiIiCizMPggIiIiVzH4ICIiIlcx+CAiIiJXMfggIiIiVzH4ICIiIlcx+CAiIiJXZXudgLYikQj27NmD7t27IxQKeZ0cIiIiskEIgUOHDqFfv37IykretmFc8LFnzx4UFRV5nQwiIiJyYOfOnejfv3/SZYwLPrp37w7gWOJzc3M9Tg0RERHZUVdXh6Kiouh1PBnjgo/WRy25ubkMPoiIiHzGTpcJdjglIiIiVzH4ICIiIlcx+CAiIiJXMfggIiIiVzH4ICIiIlcx+CAiIiJXMfggIiIiVzH4ICIiIlcx+CAiIiJXMfggIiIiVzH4ICIiIlcx+CAiIiJXMfggIqM1NIcRiQivk0EBEIkINDSHvU4GgcEHERnscGMLBs94Cz+c9a7XSaEAGP/YSgyaXoLao81eJyXjMfggImOtrtiPlohA2e46r5NCAbBuRw2EAN7dss/rpGQ8Bh9ERETkKgYfRERE5CoGH0REROQqBh9EZKxQyOsUEJEODD6IiIjIVQw+iIgoowhw3BivMfggIiIiVzH4ICJjsc8HUTAx+CAiIiJXMfggIiIiVzH4ICIiIlcx+CAiY4XATh+knuDLLp5j8EFERESuYvBBRERErmLwQUTm4lMXokBi8EFERESuYvBBRERErmLwQURERK5i8EFERBmFb9p6j8EHERERuYrBBxERZRS+ROU9Bh9ERETkKgYfRGQs3qESBRODDyIiInIVgw8iIsoofNvFeww+iIiIyFUMPojIWKEQe30QBRGDDyIiInIVgw8iIiJyFYMPIjIWH7oQBRODDyIiInKVVPBRXFyMc889F927d0dBQQHGjRuH8vLyuGVGjBiBUCgU9/fLX/5SaaKJiIicEoIv23pNKvhYvnw5Jk2ahNWrV2PhwoVobm7GqFGjUF9fH7fcL37xC1RWVkb/7r33XqWJJiIiIv/Kllm4pKQk7v/z5s1DQUEB1q5diwsuuCD6edeuXVFYWKgmhUSUsfimLVEwpdXno7a2FgCQn58f9/k///lP9OrVC4MHD8bUqVNx5MiRhOtobGxEXV1d3B8REREFl1TLR6xIJIKbbroJ3/nOdzB48ODo5//zP/+Dk08+Gf369cOGDRtw6623ory8HC+//LLleoqLi3HHHXc4TQYRERH5jOPgY9KkSSgrK8O7774b9/l1110X/feZZ56Jvn374qKLLkJFRQVOPfXUduuZOnUqpkyZEv1/XV0dioqKnCaLiIiIDOco+Jg8eTJee+01rFixAv3790+67LBhwwAAW7dutQw+cnJykJOT4yQZRERE5ENSwYcQAjfeeCPmz5+PZcuWYcCAASl/U1paCgDo27evowQSERFRsEgFH5MmTcIzzzyDV155Bd27d0dVVRUAIC8vD126dEFFRQWeeeYZ/OAHP0DPnj2xYcMG/Pa3v8UFF1yAs846S0sGiIiIyF+kgo/Zs2cDODaQWKy5c+fimmuuQadOnbBo0SLMnDkT9fX1KCoqwvjx43HbbbcpSzARERH5m/Rjl2SKioqwfPnytBJERNQqxNldiAKJc7sQERGRqxh8EBERkasYfBCRsTi8OlEwMfggIiIiVzH4ICIiIlcx+CAiIiJXMfggImOxywdRMDH4ICIiIlcx+CAiooySYrxMcgGDDyIiInIVgw8iMhc7fRAFEoMPIiIichWDDyIiInIVgw8iIiJyFYMPIiIichWDDyIiyigCfNfWaww+iIiIyFUMPojIWCG+a0sUSAw+iIiIyFUMPoiIiMhVDD6IyBcEJ+QgRfg4z3sMPojIWCFeI0gDvu3iPQYfRERE5CoGH0REROQqBh9E5Avs8kEUHAw+iMhY7PJBFEwMPoiIiMhVDD4sNIcjaA5HvE4GBUBjSxjhCJ8XqMC9SBQcDD7aaAlHMOwvi3H+PUsQ4UWD0tDQHMY37lyI7z+w3Ouk+FaI79qSBuw/5L1srxNgmr2HGnGgvgkAcKQ5jG453EXkzObKOhxpCuOzffVeJ4WIyChs+SAiIiJXMfhIgsM5E5mD5yNRcDD4aIOPmEkV9lcgIrLG4CMJ3mcRERGpx+CjDc52SKqwJBGZiU/wvMfgg4iMFfvkitcLouBg8JEEo2NKB7t8EBFZY/DRBi8YREREejH4SIYtH5QG9h9Siy2RRMHB4KMNXi6IzMHzkSiYGHwQERGRqxh8JCH43IXSwP5DRETWGHy0xQsGacChwZ2Jf9WW+5AoKBh8JMHrBRERkXoMPtrgGwqkAwNZIqLjGHwQkS8wgCMKDgYfSbCuI1VYlpxiSyRREDH4aINvKJAqLEtEZuLNgPcYfCTBNxSIiIjUkwo+iouLce6556J79+4oKCjAuHHjUF5eHrdMQ0MDJk2ahJ49e6Jbt24YP348qqurlSZaJ96skg4MZJ1h6xFRMEkFH8uXL8ekSZOwevVqLFy4EM3NzRg1ahTq6+ujy/z2t7/Fq6++ihdeeAHLly/Hnj17cMUVVyhPOJHp+OZU+hizEQVTtszCJSUlcf+fN28eCgoKsHbtWlxwwQWora3FE088gWeeeQbf+973AABz587FoEGDsHr1apx33nnqUu4C1nukCssSEdFxafX5qK2tBQDk5+cDANauXYvm5maMHDkyuszAgQNx0kknYdWqVZbraGxsRF1dXdyfl0Js5yVFWJSIiKw5Dj4ikQhuuukmfOc738HgwYMBAFVVVejUqRN69OgRt2yfPn1QVVVluZ7i4mLk5eVF/4qKipwmSTk2+ZIqLEvOxA2vzn1IFBiOg49JkyahrKwMzz77bFoJmDp1Kmpra6N/O3fuTGt96eLNKqnClg8iM7EDuPek+ny0mjx5Ml577TWsWLEC/fv3j35eWFiIpqYm1NTUxLV+VFdXo7Cw0HJdOTk5yMnJcZIMIt/gpGhERMdJtXwIITB58mTMnz8fS5YswYABA+K+Hzp0KDp27IjFixdHPysvL8eOHTswfPhwNSl2ES8YlA6+7aIWz0ei4JBq+Zg0aRKeeeYZvPLKK+jevXu0H0deXh66dOmCvLw8XHvttZgyZQry8/ORm5uLG2+8EcOHD/fNmy5sKicyB09HomCSCj5mz54NABgxYkTc53PnzsU111wDAHjwwQeRlZWF8ePHo7GxEaNHj8ajjz6qJLGu440WKcJHzEREx0kFH3Y66XTu3BmzZs3CrFmzHCfKS2wqJ1XYikZEZI1zuxCRsWLH3WHrEanCouQ9Bh9ERETkKgYfRC7gXTuROfhE1HsMPpLg9YLSwQpOLZ6PRMHB4KMtXjFIA45R4QxPR6JgYvBBpAnfdiEissbgIwk+pyciIlKPwQeRCxjIOhM/qy13IqnBkuQ9Bh9J8Dk9pYfPXYiIrDD4IHIBw1hn2NhBFEwMPog0YYdTIiJrDD6S4F0XqcL+Cs7E9fnwLhlEpBiDDyJN2PBBRGSNwUcSvNMiVViWiAzCE9JzDD6INAmx00faQuCstkRBxOCDCMf6ZDQ0hzWuX9uqiYh8h8FHEuwkmDlu/Nc6DJxWgu376r1OChFR4DH4IALw2oZKAMBTqz73OCVERMHH4CMJNnyQMixLjsR1m+E+JAoMBh9EMVQOqc/upkRE1hh8ELmA8wQRmYPno/cYfBBpwjdtiYisMfggiqGrnw/7D6WPd6tEwcHgIwleMCgdIfb6IDISz03vMfggiqHrUQnjWCKi4xh8EMVQ2drFPh/pi5vVlhEcUWAw+EiCz5iJiIKHdbv3GHwQuYBD9ftfY0sY4QiPI5EKDD6S4PWCiACgvrEFZ854G5f837teJ4UoEBh8ELmAcWz6vNyHaz4/iKZwBJsr6zxMBVFwMPgg0oQdTomIrDH4INIk9rEdH+ERER3H4IOIjBU7GBQ77RIFB4OPJFjVkSp8tc/f+AQtWBjHeo/BR1sslBmNd9dERPox+EiCdzukDGMaIqIoBh9EZKy44dW9SwYRKcbgg4goBb42TaQWgw8iF/CunYjoOAYfSfCCQUSZ4mhT2OskuIZ1u/cYfBDF0FUp8SWa9HEf6nPnqx9j0PQSrP38oNdJoQzB4IOIKMPNeW8bAOD+t8s9TgllCgYfRC7gIGPkB2xdIrcw+CAiX2AARxQcDD6S4GiXpAqLkr+FMmTIwUwJ8DLjaJqNwQdRDAYJRET6MfggcgFjGmcYDLorU/Z3hmTTaAw+iMgfeMUgCgwGH0QxMuWZN8nJlOHVWfrJLdLBx4oVK3DppZeiX79+CIVCWLBgQdz311xzDUKhUNzfmDFjVKXXVTwRSRV2XiYiOk46+Kivr8eQIUMwa9ashMuMGTMGlZWV0b9//etfaSWSiIhcwBiZXJIt+4OxY8di7NixSZfJyclBYWGh40QRBQ0bPtLHXUgUHFr6fCxbtgwFBQU4/fTTccMNN2D//v0Jl21sbERdXV3cH5FXGCSQlQzp8pExfZ54nntPefAxZswYPPXUU1i8eDHuueceLF++HGPHjkU4bD1jYnFxMfLy8qJ/RUVFqpNEREREBpF+7JLKj370o+i/zzzzTJx11lk49dRTsWzZMlx00UXtlp86dSqmTJkS/X9dXZ0xAQijYyJzeHk+ZkpVwDqP3KL9VdtTTjkFvXr1wtatWy2/z8nJQW5ubtwfkVd01b2s1ImIjtMefOzatQv79+9H3759dW+KiEiLzOnzkRkypW+LyaQfuxw+fDiuFWPbtm0oLS1Ffn4+8vPzcccdd2D8+PEoLCxERUUFbrnlFpx22mkYPXq00oQTERGRP0kHH2vWrMGFF14Y/X9rf42JEydi9uzZ2LBhA5588knU1NSgX79+GDVqFP70pz8hJydHXaqJfIZ3WunjPtSPg+GRW6SDjxEjRiQtoG+99VZaCTILT0QiQuY8d8kQjLG8x7ldiGLoqpRY2ZEfsJiSWxh8EJEvMIAjCg4GH0Rx9FzheN0kP8iUAC9Dsmk0Bh9EmmRKRZ4JQj7s9NHQHGYHUjIWg48keN6SKrwIkJs+31+PgdNKMPmZdV4nhcgSgw8iMlbs67UM3+x7cuXnAIDXN1Z6nBIiaww+iFzACyf5AcspuYXBRxscyCiz8ekIEZF+DD6SCPmvjxkRkfO6K1Oi70zJp8EYfCTB8pl5OMiYudhplyg4GHwQEQWM44YPpakgSozBB1EMfY/aWK37GR/BEqnF4IMoBlv2zcVjY5/TYClT9nGGZNNoDD6IXJAplToRkR0MPpLg9SLz8FVrsuK3py4hh00fLP/kFgYfRC5gla5XcziC5nDE62QYg28GJRf03XO0Kex1ElJi8EEUI+iVUhCFIwLDi5fg23cvQTjCA5gOln//u6fkEwyaXoJVFfu9TkpSDD6IXMBKXZ8D9U3Yd7gRXxxqRN3RZq+TYwSnj13I/2YvqwAA3PXGxx6nJDkGH0QxWGf7m64YL1Mu5pkSJPOxlPcYfCTB8pl5VB5zdt5TK9GxiY0LdF1U/HaxyoxQiZIxvcgy+CByAQMRfXihVYellNzC4IMoBitff+Px+xIjsqRYTrzH4IPIBaY3gfpZbH8MXfs5c/p8sKAGhemHksEHEfkCH13ZF2LTBxmOwUcSrOxIFdPvQkzF/UbkjOmnDoMPIvK12Ht8N24Y+GiCKH0MPohi6LqusBVNHze6Y2RIlw+2NJFrGHxkoIbmMCIchpp8xtaFUaJYNzSHORx7G0ebzZ8ThIKBwUeGOVjfhIHTSvCjv6/2OilGYguF/8R2rrR79OobW3Dm7W/h0v97V3p7fmgdcNpSs+PAEZTtrlWbGAP54RgGHYOPJIJYQN/+uAoA8MG2Ax6nJLMEsSz52QfbD6A5LPBxZZ2t5TPkqQsA4JElW71OAilget8kBh9E5Asqq9KgBxPp5I+tf8HwafUhr5OQFIMPIvK3uLld7P0kK1N6kDpg+A2zEhmQRZjenYnBR4bh4EMUNE7iiHRiD8PrdACZ83YO+ReDD6JYul619cMVKwDsPjJgEJ4Yiyq5gcFHErxgUDpYftRK1IEubpAx249d0k9PUGVCuTW9M2YmYPCRaVjpeoKd+AwT8POALTtkOgYfRDEYIviPk2OWzsU5+HfNQc8fmYDBB5ELAn+9MoTd3Sz72CWTOnCyrJIbGHwQkS8oHecjk6IJIgMx+EiCz+mJUmswaD4QNx6J+KFWCPqrxOR/DD7aCHqTI+/3ktN18QpqsVr+6RcYOK0EMxd96nVSADh728XeJIuZc+YEv08LmYDBRxLsMU6U3LQFZQCAmYu2aN+Wymti0J+6BDx7aWN85T0GHxmGz7q9wbtJfdLdtTwy8bg/yA0MPjIML4LJce/I8W8sezzhsucETyGi9DH4IHJBUK9Xfo09ZIMm3wVZvkuwu/gygfcYfCQRxALKxy7kX951Bs6k1o5Myit5h8EHkQuCWqGbFsza3c9O5oOJLh/Am5JYwc4dmYLBBxE5ZlboYZ9s0GRYjEXkeww+Mgzr0OT0tVDwfjJddo6Nk1aJoLdkyGKndHKDdPCxYsUKXHrppejXrx9CoRAWLFgQ970QAtOnT0ffvn3RpUsXjBw5Elu26B8DgIg8YEI0m+6rtrzWErlOOvior6/HkCFDMGvWLMvv7733Xjz88MN47LHH8P777+OEE07A6NGj0dDQkHZi3RbESonNx8mpPOQBLD6Bkc5p4Id6gad5cn44hkGXLfuDsWPHYuzYsZbfCSEwc+ZM3HbbbbjssssAAE899RT69OmDBQsW4Ec/+lF6qXWooTmMnOws4zrHeYEnnTe4392haz/7ueaoqm1AQfccZNmcypdlldygtM/Htm3bUFVVhZEjR0Y/y8vLw7Bhw7Bq1SrL3zQ2NqKuri7uT6Ut1YcwcFoJfvfCeqXrJSJ3L8pqZ7WNWW/AL7bnFS/GKf/vDdvLsw8MuUFp8FFVVQUA6NOnT9znffr0iX7XVnFxMfLy8qJ/RUVFKpOEx1d8BgB4+aPdStfrV2z88UZQq3PTWhPt7uegz9tk2GExTlDPRz/x/G2XqVOnora2Nvq3c+dOr5NERBmEd/rxTGgJWrBuN/702scp37wRQuAvb2zGS2t3uZSyYDlQ34RbX9yAtZ8fcH3b0n0+kiksLAQAVFdXo2/fvtHPq6urcfbZZ1v+JicnBzk5OSqTQeSYrtcMTajQdTDhBttJ8CD72MW0Fh6dTCirNz1XCgAYfkpPjPyPPgmXW1mxH3/7snV7/ND+biQtUO54dRNeKd2D59bsxPa7L3Z120pbPgYMGIDCwkIsXrw4+lldXR3ef/99DB8+XOWmbMucKsOeDKpDyQVulicTLop+EZTHSgePNCX9/kB98u8puW376j3btnTLx+HDh7F169bo/7dt24bS0lLk5+fjpJNOwk033YQ///nP+NrXvoYBAwZg2rRp6NevH8aNG6cy3URa6Lq+ceAmdzjZz7K/8MOhTOdRkkmPoVK1ODkNfv1wDN3gZYgqHXysWbMGF154YfT/U6ZMAQBMnDgR8+bNwy233IL6+npcd911qKmpwfnnn4+SkhJ07txZXaqJiBRhYOhfQWnhyUTSwceIESOSnqyhUAh33nkn7rzzzrQSpgofM5AMFhc5blb+du7Ibb/tEvADzYsy2eLhieD52y66pXNTwxuizKPtsYum9VL6eGzimVTvpbo0On7swqMOwNubrcAHHxSPd0Skku4bJ+np7m0uH3seSG8j4Bcuk3KXqnyxNvOvwAcfQW9eJcU01bwm3U0GjZN9K1sv+K0aCUq957d8CCGwYVcN6htbvE6KLV7u38AHHxTPbyczUSttARwDw3gG7Y9ULbVOk6qrLL26oRI/fOQ9XP7oe3o2ECAMPohcENSmevMG37K3n+MGGZM8NkFvxSqvPuR1EqJSFa8uHTu4kxCb5n90bKTVT6sPe5wSe9jnQ6N0+jgE9YJBlOlk6wXjYqwU0kluYS6HRSD9Ah98EHkl7pV0xrGucNIqEfSWDFm+uunyWVBoGi9bLhl8EMXwVcVrADerLpVBQvxjF4plUjBm3mO9YOFjF6KAM6g+V8qEa4Pb+zaox7KVSfkzoHiRJgw+iGKYdNdH7rAzvHomjY/jp+HmM+eo6MFXbQ3lo3OQyBOuzmqrcnh1yd9k0uM4k3JqQsuaDL89JvIyqA588OGzsqCd306OoAhqIOvXFoGgnwZp5c+gsurX8kWpBT74IKLM4cbbLn56LOGEn3LHmyn/YvDRhp9OPPKPoDbbm1b3O9nPdn7DO3Bv6CpfQQ8gbWOfD33SKbymVaykH+skc/HY2BeUFgGZXDCgkMdXbTVqWx6PNoUd/1YFme1T+hqaw4hEvK+UeNz1ib3oNDRHHKxA6+K+Y9JFXGZWWxOSHYyQzx2BDz5iTX15AwZNL8GmPbWebP+dLV9g0PQSPLDwU0+2D2TWyVFzpAkDp5XgysdX2f6NyscjxW9+Ev33df9YiwP1TcrWbQrTytO4We+h+I3NKZeTHXw2IA0JthhwDdfOhEDFBHzVVqPYnfuvD3YCAB5dVuFJWqYtKAMAPLx4iyfbzzSLNu8FAKz5/KAn21/4cXWb/1d5kg6tDLwqP77iM6+T4GteX5jjW17sl6+w4oTf8uJ6fPUPr+N/532YsvV00cfVGDNzBT6pMmdSPtMFPvggMoXXlboO5oUe8oJ4XNLhdefo2OOR8rFLzPcbdtUoS0M4IvD8mmMz1C75ZC/KUrSW//ypNfik6hB21xxVloagY/DhIhM6gRmQhIwVxGtcEMqT7MU26MGK1/mLxCRApni1hO0nXDaLYQP6jenAQca0cr5zg1nciNTRXXV5fReeibwOPpxuXuXNXds1BbUUss+H2zwqSQG4SQw8nRWv15W6Dm625iXaf+nu1iAeFz+Lf+zCWlMnBh9EGSCId/FBuDRIH5XgHUajyDx2iUj0D4klP6qt3PJ+wccumSIINTVRDN6Ymimd42LSOB+p+CmtJmLLR4YwoZ7mMNHJ6azKglhPBqE88QIWz+u9IfO2i3tp9XqvBE9GBh92m79ZKZFKgSxNLsYeuh5b8TSP5/X+iD3OKYMPh2/GkPcCH3yY1CzMzlMUNH4t0elcX/3QdyedFimv8/fBtgPRf6fKR2ygtOPAEdvb8DqPumzfV+91EmzzdfDR2BI24v1rIQQamv0xd4eJ8U/GzHvi9S2lBiaUJzd2qwn5jKXznPGymK749AtcM/dD28tndzh+CZvy/HocaWrRkSxfnLqHGpox4q/LpH7j5Q2xb4OPo01hnHX72/jBQ+9I/1b1c+qr53yAgdNK8MWhxhTbpbb+vX4PBk0vwbz3tnmdFAD+qGRM4mafD13HprHFX8Hvyq37MGh6Ce5765PUC/vMm2VtpiBIUbxyO2fH/X//4eDNn2RXZW2D9G84q60D63YeRGNLBOXVycfSd2PnvrNlHwDg9Q17XNhasPz6X+sAALe/+rHHKSEnTGsRcMKAxlMpd3x5rsxaqmeOKpN2h1fFy6R9EFS+DT7ssipEtjucqk0KZbgglqcgBB+y2DqmT9vylOqxQNtDkcnHxm+nYuCDDzelKvcmVNQGJCFjBbFiDMartqmXMSmfuuuRIJbTdAV1l3CcD43MqTLMYEIARMHhZnnSdQEI4psP6R0X7/ZH22TLZsP+MAqSKw4oL/dD4IMPk5h090SJ6DsbOW6MmaSH2taTDGP4qZi6lVY/7RO/yMjgQ1dBSrVeM1odjEhERgpi/WXC2DXptlzYeuzifTajdO9zP5fTTA4SnBQLPnZxwG4rQnpzHDj/ralMqkTJ/9wsTrpajgL52CWN33rZQte2ftL3qI285tvgg9LHxwDt6dwlQdzdWT4NZmOPhZ3jEsRjl4i/supOaoNaV3Kcj4DwW/EM6PlkrCDubhMeu7jN6wuR7j3uZfbatmh7va/9xV/nYkYGH17VlyZU1N6ngIIkCOXJf30+vE5BADCo8VxGBh9elTuddYaTeW68Pv0yZk6XL/npLi4cEZ4PO360KRz/eETTdprCEU1r9idd5dTJ+Z4qJW2T6p8zjDIy+LDPH0W5vrEFZ854G5f837sply3bUxf9t5cXw1UV+zFoeoln20/EH0dcvzEzV+DsOxamnDAx9i68bHetsu2/s+ULDJpeggcWfpp64TQP2vjZK9EsEYAEvYzoyN+/PtiBQdNL8MKanRrWrl/Qj7kXAh98ZGcFPotY8/lBNIUj2FxZl3LZhxdvcSFFqf3pNc7lYrItew/jaHMYG1MGFMejj0eWbFW2/RmvbAIALP/0C2XrTKbKwaRcXrHz2CWtRzMarrRTX94IALj5xQ1Jl2uXbo+u+n5qpXQqEhFYWu7O+WUl8FfmbjnZqRfKUME/vUi32ItFEF9ZzUR+Oort53bRM29XEGOR9btqPN1+4IMPqzsAfYOMJV+xro5i7H+mTibc8ajk5qu2fA36mEwaKZkBrX2y56JsH0HVgh98eJ0AF8gENXF3qjyvXRXE/R17IVSaP4NPXD8cx3QCFJMCcFOGvvdDEGTC25QyfBt82N7PaRwQ1eegCWUjNgl+OKHIW6mKrAll2om2ZZ/nwnFe7gnpieQcJtag+EoZ2X3n9bnr2+BDF51Rv67mUpn1mhIdG5IMVwXxAufucQze/nNC9z436cJsUlpIrcAHH+mcp365QEo9don5N09sUsmr4qRiu346F+yc7n6pu9oy5eaoHR+UD/ld5+2+Dn7wYVBZNiEtJqTBZDrrGD9d4OzS1efD5GIaxBasWCblL/UgY87SalIeVfFbR2Tlwcftt9+OUCgU9zdw4EDVm7HNbwfECZkcZsL+MJVfqjupCj1DipNRgaPmOwij8irJz2l3m9c3oloGwTjjjDOwaNGi4xvJ9m6sjXR2sPIOp2pX52zFfNslKe6TeKl2R5aLNRiPjTtM2s26+uAF8SbM62BClpaoIDs7G4WFhTpWDeDYyGyNLfaGQ9Z9PKTmK5AsHS3hCCIC6JRtv4EqEhHIsvnCt52mRyGO7evOHTvYTkMmONoURpdOcvvELxdPpw0fqYZiT9dRzeu3TcNxbGwJIzsrCx3cHDglEZ+UU8AqqXYHGfPXPFh2+C340NLnY8uWLejXrx9OOeUU/PjHP8aOHTsSLtvY2Ii6urq4v1Qum/UeJs75wFZaLAcZU1SU1u04GDc/icqLixACF9y7FOfetUhq3olU87vIls9rn1yDgdNKUFl7VPKXwfXXt8oxaHoJVm7d53VStJApxrHn17tb9ymd36Wtnz7xAd7eVKVl3V5eXBqawzjnzoUY+cByW8vrvsaY1B/ClJT44cbB2M66CSgPPoYNG4Z58+ahpKQEs2fPxrZt2/Dd734Xhw4dsly+uLgYeXl50b+ioqKU20g934Q77n/bxqRXDjW2RLCntgG1R5uxpyb5hT+2CfHjFPO7yA4ytuSTvQCAlz/anXphw+g6FR9ZemwOkzsDOj+NTFN323386DI187skqkhvean93CB+uDAkU151CPVNYWzbV29r+Ux61dYr7YZtNyYMkpPsXPY6VFEefIwdOxb//d//jbPOOgujR4/GG2+8gZqaGjz//POWy0+dOhW1tbXRv5071c566LdoULe4txM8TIcbgp4/XeRaPoJxfpk0qqfXTNoTqQ6LW4fND8XDb2ei9p6gPXr0wNe//nVs3Wp9R5STk4OcnBxt29dZN7Zdd6ro2GlSUhV8p8Ore8mUdLSl9VVbo6r1xNJ52cUPlTSQXjpVZ1H2XLA3zoehJ5g0PQXKL+VUhlWWhEhcvrwuI9rH+Th8+DAqKirQt29f3ZuylNYcBwrTAZgxsVz8IGMBPANjBKX6dZvJQVIQi6xpb154WS+kW0dqmzRUz2q1MzndyoOP3//+91i+fDm2b9+OlStX4vLLL0eHDh0wYcIE1Zuyxc1Zbf0gNtoN+m5wkj+Zild60iuf7HCplo82J5iqPJp1OTaLnTvWdPafx5Odxkn52MWlWsyvN2om9/lQ/thl165dmDBhAvbv34/evXvj/PPPx+rVq9G7d2/Vm7LFzR2c8vGI0/Wm+F6m+czrAkfBEpTWfZlLi9cTTjp/fCs8b2pPJd1WILuHxp+hhDyT86k8+Hj22WdVrzItVuea4ecfAHf6cfg0mLfNyW6RqZxNfjyhimxAnQn7hNxjSmkyJR3JWLVymFzHB39uF4Pu9aUubG4UGoMLpgq6H7uYtG6V5B676EuHXUEKeOyUEd5sxHArTz7ddyafG4EPPqzYPQn9Mry63y4WZDaZCsuk4D4d5vQtULMe68EV/UdXvyrp9fpy75kdcAY++PDrxVbbmzFxHU4NLpkKmHboTa4IYpkQzJp83qo+b+T7fBi8c1zmWsOHD85d2TR6fY4FPvgwSYsB3chlRzi1SwihfW6PVnbn07GbPbfSbYJIJPVxcjq8ul1S8yE55Ndjamvf29jnVot48dhP9ljLjp3kFj8EH1ZMTnfggw+dvbvf2RI/t0ey4xyOCJTurHG0neq6Bke/s5JOT/lkfi45B4zTu7ey3bUYNL0E/2/+Rke/b6v2SLPj+XlMPrETGf/YSgyaXoLaI80Jl0nnIvXWpuqk3y//9AsMml6CmYucTU1gN20Dp5Vg76HE541Jxy5u1GGdfY60rdna4s3VGDS9BLOWqhlyX6UgTixnJVk+vW5BC37w4XUCvrS/vlFq+diC8UqpunlVdI3zsfjLOWBeWrtL4Vrbm7loCwDgmfcTT1bYys6xX/bpXqMuRLqt21EDIYAVW75IuIzM7ujUQa4KmbagDMDx46jTa+srJZb2z8Bafhnh9A8vH7tBuO+tcsfr0DW8+sqt+6XW+4un1khN8GkKk+u24AcfaXS8UnkXkpVGZaCyAGXSCKdOcidzmGTX75e9LVMsOmTpqUJk7srcLsY6t6d1eH+/FMAY8ueYvV84mZx0ZUXygMVvvI5Pgx98uLitZCd3B8kjbcqzzkyj9cLil0NqQIdTt3l5bGT7YdnZ59Y3XeYXwHSLk87j2NTiw5YPrxOQROCDD1MY0/IRW9GpW61rpAZfs7W+gFw9Jalq/ZPdexm6u5Wy0ypkwm5WkYZU5a9tQOWbAN8lJrduBz74yMoy4TQEQmnsaV13LAaXSyUcDTKmMSTzw50nIDnOhxmnl6u0PhrRWf58UPzSnlhO6/4zewdaJc/kFAc++AhC3Wh4mc9YpldGTslkK50WPTom/m0XG8tzlycU0FPSMZP3R+CDj3TOVFOOW6p0yF0EOciYV0yuCGJJJdO0nexDbgUTfil/MoKYJ6UM3j+BDz6CUDdqO8ECPqaFnSQHoXyoprNFJ+0OhUpS4d56dbAVrASkecSkesegpNiWdJwPvu2il5s7ONmBTuckCmLzvtcFX4XAvmorsayugYp0vvKcLtXnY/zr73aWd7bPg97SqZvp1bDV8TU5zcEPPgJwb8u7veNkjqaTG0STT1aVkl1AY79KtTukB8gyJOqUfUtCZ7nQsUush1dXv52kaXCQr7blI1XA1PZbvQOBWaelxeE2m1oi+GjHQYQ1TruRbM1eXxuDH3xYve+e4IiEIwJNhoxip6uicDq3S0tE+O49d0dvu2gd58Peyt2Y+ySZ2Ao/VZLNCCX0cjpPjBBC+limutgebQoHotVQlzte/dj1bf7pNWfbnPJ8Ka54dCUeXOhsqgE7ZFrpdtfYmxpDlcAHHzIufvgdDC9eEv2/yumcF29OPudFWxt2HR+BL5JqrAWJ2uiLQ3LDvLd6aPEWfPPPC30XgKQy/ZVNcf9f9Zm+UQz/b0nqOS6mv1KGQdNLsN7hPEBKxBS3jyvrki4alLddPqmyzueug0cwcFoJfvn02uhnduuFP7y0EYOmlyRct5Vk65656FMMml7Sbk4pK+mM7GwS2TrY6fxZ6Xhy1eeOfvfahmPD/z++okJJOmRftW1bRjY5GPU1HYEPPmSqxk+qDmlLx60vbZBa/q9vH58PIdUJmH9CJydJkn4GXNfQgp0HjzjalhfsHPsD9U3ON6ChNn/qy4rM6cRrqqWaFExX7FF3NPHEdzo8vvwzy8//9cGxOYR2HZS/K3xuzc6k6z7O3nxL6c6HE8S+Y7ry5Gb/CVNGVXa7dAQ/+HCzw6lHJ7fTPDpJrp/qLx8l1Six+031CKd2FeR21rRma4lyqeK5eMp9GPcoVE2p9fp5vlNtU+2n+sapDMiipeAHHy6ehLomgEv12MVdJqXFW3pHuvSOTIdT2cjX7tIyAxOruGB7OZV97KOroJYpp0xKs7aO/1rLHl+19Y6rr9qq48YcLCad2HapntuF2pPpcKpLwn5MLjd9u1FBx71qq7FLletvu7hwBrqZJW2PXXSuJ8nKvW4dC3zwkc7uNea9eE2vAQbxGXCsYOdOn7iWD9WPXQyNCL2cUyX+RoOlNlbQ6yhAc5+PJN+x5UMzN8cVSFaI0hpkzKAKSdWJYug1iNCmz0eKZXW97eL2fJBeXuNk53axtU7LgT7UrFsrVgxKJStPXu/qwAcfbvLb7LNBv6nw+uRKh5fHRkh0+tAV2ydsEta0vUT9qtwoQ7H7UFX/LuvYw38nfMoUu5gl0/efVSuRTJ8Pt+sc3wUfss1w1ndQ7gcJ6VTSJgUJBiXFc0FtEpYa4VRXIhKtWNEub3voZFYrPf6PpmWJUkleVtnnQyt353ZJ8p1krRLXFMsqiTwi85qoSm5Xi16ObOt01GG761S9bp3atXgZlGY/7L+2TE6y74IPVwuA9Ain6hIX3xSrbLVx/Hgyed1DO5YPd1+U3WOvOo92j57UxHIKEpnw3HW5V57Wjq/a1mzNjV3n5o2ZH893k1tnfRd8mMzg42zJyYlrN49ez08CAEcdzsmhk1eVgczxiH/bRUNibHA7yEw8yJjVst69GWOXG/tPpkw5HgjRl5d8b1jtKR1dARpbwkomw/Nd8KGiKOqqUB9JMhR1SxoHS/aCVVJW6XhbqaScZVIcG5J70PQS6fls7Fj+6Re2ltu69xBuW1DmaBtrPz9oa7nP9yceav7D7QcsP7/2yTUJf/P/vXN8GG6VRXRZ+V4Mml6Ch2wOzx17jLM7JK+hdJ1LWQlqJl39J1TPLPo/f1/tKB0mTGwIAIcbW5J+//JHuzBoegmeXm1vXhNz2ivtKX7jk3afydbDemfYTV/7Y5I6fw3NYQy5422Mnrki7e37LvhQwQ+xdDrPgX/3/Hpby+mq6O5769i8NH94eaPyddttzfjbilTzaSR231vtKx4rye4c7n7Teh1LPtmb8Dd/fn2zre3K+uP8Y0HYgzbni4ktFzd+7zQdSUrJ9ZYPxU9dVlY4m6BQ2Z1+mhPLfbAtefqnfFnH2A3w7Q55IPsGhq46bN7K7Wmvw8l8QG5yMgxF2e5aNDRHsHXv4bS377vgQzb69NujEBXsZlnHrjFld7tx3P1SthK1IiQSm638E3KUpsXdcXfsHyCtj8MkLqC6+nf5lZ93h9uPWC1ntZUY54Ov2jogH5DYW96Ugp86HS52ulJ1J+L18HpJqDgJTcqdbCtC7Pnh1bxCiYqHTE5kkh6WClQkEgG51gyt83x4WKGZdD5kEqlxPmysT2W17bvgQ8X547e7C10XAB0VHTuIqaHy2MiOFhq7ZdVlxP7bLtZLyvX5sL90ojrBlMc/Stbt8rkZu+d0zbztZkDll5bOWDJpZsuHC/xWhqTvtDS+RmlSZZCMIckwgtUQ6MkuRGa87ZI+kyveRII6fYFJr8i7xYQipStYV8F3wYfWd/oNEj/ImH/4Ka06GfVUSTotsY9dlKbENhX7TybtMq2LsrtEKgjSeQZ5eXLaPJ7t+iEYVKOYlBZr7dOX7Frn7DKormLzXfChgg9ijzZvu+h67OKXlbrL/zmIZ9XyYfdO1LM+HwrWIffYxXpZt4NIdSOcOnmR0n+CmCeVTHnN24rvgg+rCkVXBzC/XkftV7ryGUw5zof0GvVoN3eHXw+mAlZ9Puw/dvHmeCfs86GpY2gkQTOJ249/vAr2TOXn3WFC2uVa3dzlv+BDQf8HEwqFDJPS65c+H14z6Rm39NsuCf6tgt2WBNlOslZkAhVTOqGrSoblyKw+eNul7SCC+w43Jl3+8/31DlMEvLPF3oCFrVTsv3U7DuKV0t1Yk2AQwrb++f7n+NncD2wNWGadPrVvWmX02y4q+O3uQtezRh27wfznojYEIAuxZCuMuHEnPLsqWyfa8nFCgiSqaHGwnqDN/NdhLdPt9tsuMYmwUwbrG1uwaU9d3GczU4zKKzMw356a+EG/fvrEB7Z/60z7/X35oyvxm2dL8V+PrbK1hj/OL8PS8i/w1y8HbpROgcF1WeCDj4aW9iNimnxAWjXEjORp0tsuqrbtRKPFsUyk7XH3wzGPdbixxfaQ383hSNI7I8u3XZKsOvYipTr2cLNFKFnSjzS5Vz6SXfSPNoWlHnO5panFOh0Nzc7m9bDzk7bHRLUvDiVvRUnFy0Pz1qYqR79jnw+FUkXCsSprj0aHlo5ltY9lT/pUzYHpeHNjJT7cfjD6f6dDNafriXe3Sf8mdi/WHm1OuJzsJeho07E5BeyoPdKM1zfom9/GljSvset21OAHD72TcrlwRODbdy/B8OIlCS8K6bR8pDorHlhob8h2WYkeu6jo87H/cCOufDz+zlNmEsJ/r9+T9Pu9dQ1x/080B9D2ffUYNL0EF9y3NPqZssculi029n//y6fXoqklPqCtb2zBWbe/jR8+8q50emLXtbJin+Uyujv3prt+L8NCp9tW/aaVykPku+DjseUV7T5LtMteWLPL+gvbLQOJF9R5cbvpuVKp5Z1GrKl+96fXPpb+Tez3bSuvdGz8ck4BO5aWJ54/xQ5THh2VVx9Kucz++kZ8cagR+w43ouZIk+Uylm+7JKlFdN6J270AqLgQJTqOb5a1v4vsnpOdIB3tE5LqBmj+ut1x/++QIJKymj/Ew1He29nepj/FB9sPoCkcafdoRNafX9Mzh5FpVB5LO+tKd3Ns+VAoUf1lysUlEbcidCf7IdUw1Lr2rcw+sX7eHXyJ3hCxfNvF9mMX1cGHzQnGFNxjyTwdkBlePZW2WUw4eqpFFk3qj9Y2earueo0aA8dgKqY5kBvaX3LdaZbVYAcfEoXcoHM+/YpXY15SngQG7UcveVG/JtpmluSrIzpHOO1gs8ZRcoGSSHyLws4tbc/fxK/xyvXFSZfsxaLtMVA1KaCKeXscbTfNLXjZH0dH67aTVaqcGDLgwUeisQLSX7fOgpju8bUbJTvJgmfjPmharxWTAtGUYtKaqNxYPXaxS3WH0w4Op1ZvJdWBTuJzJ50oE2mb9paI9eNCnW+kqAlsQkn+51yiIED3jMeuDxiX7u8VvGigc5yPdOvJQAcfiVh2OFWwDlPYbUJ2UnhSvW6uKyiTe+xiVfGafMScs5Mr2ccusfzw2CXRBTth2i0+V9VZ1zod9pcNaDG1xfSnMZ52OHXlsYvcNtLdH4EIPnw1/LgN6Z6EOtOd6g4xkyvPWF481050wZZt+dDZ4TRR58u2VDTNm1IWEwY2Fp/pHOcj3XWoKtO+7fPhcnkSCf4ttQ7FLR8Z/baLjIRNtxZHRDrq01gQdTc/tnLSxOvVYxeZYu/Xus2JuMMh8dgl6dsuMUdRdTm32/0k4SNTiW0lWlamJcJJP4G2aZdpVdHZGV72WLbvcKqoz4eStbhP9ti4/bZL+huRXJwdTuWpOI6G3FSlxUnZSVVxm3K32ZaKi5aJ7FSIWRZnedKOaDHfKR9kzMXb3kSPXVT277BDZnteDXZmAuPH+XC75SNmg3aOXbrznkkHV1JLt6ct+Jg1axa++tWvonPnzhg2bBg++ED3ULbtyT03NofJdwaqXrWVrQjSfdXWbSbN7SL92CXm36r7fNjtcKpibpdERVH366xtk57onLEKxIx61bZN+pSdVyacoA64PdNA3GMXp2+7KH7V1vi5XZ577jlMmTIFM2bMwEcffYQhQ4Zg9OjR2Ls3vcGfZCV+7JL+urV2YDT43ExVOeraLTK7xKQLv25x+zvBvpfv8xF7x6WW7T4fCraV+LGL/Vw5qWzbjfMh0+dDfnMJ0pD+2y7pjvMh229Hd9yV9qu20i0D6jLkNPDR++p2er+3HtYvTQ888AB+8Ytf4Gc/+xkA4LHHHsPrr7+OOXPm4A9/+IOtdew+eAR1YXvJ211z1LJS25dgLP+6hmbsOhg/5HFLuP2e3He4sd1yx7+zHk0y0fIyyx1qaJFavrqufT7tpKO6rgE9unZMuVysvXUNSdd9sM0om4mWrT3Sfuh1VXnce6ih3We7Dx5Fdgd7lU/t0fblI5FEy9UkGVre7rrtLFsdM5T37pqjONTYfruHG9uXp72HEpftvTHnTe2RJqn0AqnKdrPN5dqnGTg2/0fb31XVtj/eALD/sHXarcpSovS0Lc/Jlm3VduqFgwn24f769utOdX7ZTcdBi3VX1jbYPgeslq+2OK+SpSHR9AqJzq9Ey6sqf9V17dO/88AR2wHmwXq5c6GytgHdEoyc25bVemOvSYnKUCyrcl2dpDy1PW+SXe+sfrO75ig6tilPh+rs75+QUHwL39TUhK5du+LFF1/EuHHjop9PnDgRNTU1eOWVV+KWb2xsRGPj8Z1WV1eHoqIiFN30PLJyuqpMGhEREWkSaTyCnTOvRG1tLXJzc5Muq7zlY9++fQiHw+jTp0/c53369MEnn3zSbvni4mLccccd7T7vlJ2FDtntnwo1WswXkmOxHAA0hSMJm4asftN23YnWm2zdiX4js26ZPMquP3Y52XUCx46L1Y1C67KhUHxznN39kSo9Au3nikm0fEQINLdpyZJJR6q02NmHidYr+5tUxyj2dzLrzc4KJX0Ekmqdydad7r4DjnXSTDTqqJ1zFzj2iCfbIo/N4YhlM7bVeq3KXaJlk63fbppV1guq6xzZfeHGOSCzvEwdIrPeRPm0qitl6xuZ+tpq/TLnb6cOWSlbgWL3odW6w2H7PTm0PHaRMXXqVEyZMiX6/9aWj4+mfT9l5ERERERmqKurQ95f7S2rPPjo1asXOnTogOrq6rjPq6urUVhY2G75nJwc5OTkqE4GERERGUr52y6dOnXC0KFDsXjx4uhnkUgEixcvxvDhw1VvjoiIiHxGy2OXKVOmYOLEifjmN7+Jb33rW5g5cybq6+ujb78QERFR5tISfFx11VX44osvMH36dFRVVeHss89GSUlJu06oRERElHmUv2qbrrq6OuTl5dl6VYeIiIjMIHP9zsi5XYiIiMg7DD6IiIjIVQw+iIiIyFUMPoiIiMhVDD6IiIjIVQw+iIiIyFUMPoiIiMhVDD6IiIjIVQw+iIiIyFVahldPR+uAq3V1dR6nhIiIiOxqvW7bGTjduOBj//79AICioiKPU0JERESy9u/fj7y8vKTLGBd85OfnAwB27NiRNPHnnnsuPvzwQ1vrlFlW17rr6upQVFSEnTt32p6zhnl0vqyudZuUR137g3lMPy1+zKMJ55fs8syje+mws2xtbS1OOumk6HU8GeOCj6ysY91Q8vLyku7oDh062K40ZJbVve7c3Fwt62Ye3V23CXnUuT8A5jGddfsxjyadX8xje37KY+t1POkytrZooEmTJmlZVve6TUkH85jeuk1Ih879IYt5dGdZJ8ubkA7mMf3lTUiHyjSHhJ2eIS6SmZLXT4Kar1jMYzAwj8HAPAaDn/Iok1bjWj5ycnIwY8YM5OTkeJ0UpYKar1jMYzAwj8HAPAaDn/Iok1bjWj6IiIgo2Ixr+SAiIqJgY/BBRERErmLwQURERK5i8EFERESuYvAhYcWKFbj00kvRr18/hEIhLFiwIO776upqXHPNNejXrx+6du2KMWPGYMuWLXHLVFVV4ac//SkKCwtxwgkn4JxzzsFLL70Ut8xXv/pVhEKhuL+7775bd/YAqMljRUUFLr/8cvTu3Ru5ubm48sorUV1dHbfMgQMH8OMf/xi5ubno0aMHrr32Whw+fFh39gC4l0cvj2NxcTHOPfdcdO/eHQUFBRg3bhzKy8vjlmloaMCkSZPQs2dPdOvWDePHj2+Xhx07duDiiy9G165dUVBQgJtvvhktLS1xyyxbtgznnHMOcnJycNppp2HevHm6swfAvTwuW7as3XEMhUKoqqryTR5//etfY+jQocjJycHZZ59tua0NGzbgu9/9Ljp37oyioiLce++9urIV5Vb+tm/fbnkMV69erTN7ANTkcf369ZgwYQKKiorQpUsXDBo0CA899FC7bXl1LjrB4ENCfX09hgwZglmzZrX7TgiBcePG4bPPPsMrr7yCdevW4eSTT8bIkSNRX18fXe7qq69GeXk5/v3vf2Pjxo244oorcOWVV2LdunVx67vzzjtRWVkZ/bvxxhu15w9IP4/19fUYNWoUQqEQlixZgvfeew9NTU249NJLEYlEouv68Y9/jE2bNmHhwoV47bXXsGLFClx33XWByiPg3XFcvnw5Jk2ahNWrV2PhwoVobm7GqFGj4srib3/7W7z66qt44YUXsHz5cuzZswdXXHFF9PtwOIyLL74YTU1NWLlyJZ588knMmzcP06dPjy6zbds2XHzxxbjwwgtRWlqKm266CT//+c/x1ltvBSaPrcrLy+OOZUFBgS/y2Op///d/cdVVV1lup66uDqNGjcLJJ5+MtWvX4r777sPtt9+Ov/3tb9ryBriXv1aLFi2KO4ZDhw5Vnqe2VORx7dq1KCgowNNPP41Nmzbhj3/8I6ZOnYpHHnkkuoyX56IjghwBIObPnx/9f3l5uQAgysrKop+Fw2HRu3dv8fe//z362QknnCCeeuqpuHXl5+fHLXPyySeLBx98UFva7XKSx7feektkZWWJ2tra6DI1NTUiFAqJhQsXCiGE+PjjjwUA8eGHH0aXefPNN0UoFBK7d+/WnKt4uvIohDnHUQgh9u7dKwCI5cuXCyGOpbdjx47ihRdeiC6zefNmAUCsWrVKCCHEG2+8IbKyskRVVVV0mdmzZ4vc3FzR2NgohBDilltuEWeccUbctq666ioxevRo3VlqR1cely5dKgCIgwcPupeZBJzkMdaMGTPEkCFD2n3+6KOPihNPPDGaZyGEuPXWW8Xpp5+uPhNJ6Mrftm3bBACxbt06XUm3Ld08tvrVr34lLrzwwuj/TToX7WDLhyKNjY0AgM6dO0c/y8rKQk5ODt59993oZ9/+9rfx3HPP4cCBA4hEInj22WfR0NCAESNGxK3v7rvvRs+ePfGNb3wD9913X7umbi/YyWNjYyNCoVDcIDOdO3dGVlZWdJlVq1ahR48e+OY3vxldZuTIkcjKysL777/vRlYSUpXHVqYcx9raWgDHJ25cu3YtmpubMXLkyOgyAwcOxEknnYRVq1YBOHaczjzzTPTp0ye6zOjRo1FXV4dNmzZFl4ldR+syretwk648tjr77LPRt29ffP/738d7772nOzuWnOTRjlWrVuGCCy5Ap06dop+NHj0a5eXlOHjwoKLUp6Yrf61++MMfoqCgAOeffz7+/e9/q0m0JFV5rK2tjZvAzaRz0Q4GH4q0FpapU6fi4MGDaGpqwj333INdu3ahsrIyutzzzz+P5uZm9OzZEzk5Obj++usxf/58nHbaadFlfv3rX+PZZ5/F0qVLcf311+Mvf/kLbrnlFi+yFcdOHs877zyccMIJuPXWW3HkyBHU19fj97//PcLhcHSZqqqqdk3W2dnZyM/Pd+U5ejKq8giYcxwjkQhuuukmfOc738HgwYMBHDsGnTp1Qo8ePeKW7dOnT/QYVFVVxV2UW79v/S7ZMnV1dTh69KiO7FjSmce+ffvisccew0svvYSXXnoJRUVFGDFiBD766CPNuYrnNI922NkPuunMX7du3XD//ffjhRdewOuvv47zzz8f48aNcz0AUZXHlStX4rnnnot7VG3KuWiXcbPa+lXHjh3x8ssv49prr0V+fj46dOiAkSNHYuzYsRAxg8hOmzYNNTU1WLRoEXr16oUFCxbgyiuvxDvvvIMzzzwTADBlypTo8meddRY6deqE66+/HsXFxZ4OsWsnj71798YLL7yAG264AQ8//DCysrIwYcIEnHPOObZmOvSayjyachwnTZqEsrKydq0yQaIzj6effjpOP/306P+//e1vo6KiAg8++CD+8Y9/KN9eIkE/jjrz16tXr7jz8dxzz8WePXtw33334Yc//KHy7SWiIo9lZWW47LLLMGPGDIwaNUph6tzF4EOhoUOHorS0FLW1tWhqakLv3r0xbNiw6OOFiooKPPLIIygrK8MZZ5wBABgyZAjeeecdzJo1C4899pjleocNG4aWlhZs3749rhL0Qqo8AsCoUaNQUVGBffv2ITs7Gz169EBhYSFOOeUUAEBhYSH27t0bt96WlhYcOHAAhYWFrubHioo8WvHiOE6ePDnaobd///7RzwsLC9HU1ISampq4O67q6uroMSgsLMQHH3wQt77WHvixy7R986C6uhq5ubno0qWLjiy1ozuPVr71rW+5GgSkk0c7Eh3H1u90050/K8OGDcPChQvTWocMFXn8+OOPcdFFF+G6667DbbfdFvedCeeiDPNvRX0oLy8PvXv3xpYtW7BmzRpcdtllAIAjR44AQLsWgA4dOrR7SyJWaWkpsrKyXOldb1eiPMbq1asXevTogSVLlmDv3r3RO4zhw4ejpqYGa9eujS67ZMkSRCIRDBs2zLU8pJJOHq24eRyFEJg8eTLmz5+PJUuWYMCAAXHfDx06FB07dsTixYujn5WXl2PHjh0YPnw4gGPHaePGjXGB4sKFC5Gbm4v/+I//iC4Tu47WZVrXoZNbebRSWlqKvn37Ks5ReyryaMfw4cOxYsUKNDc3Rz9buHAhTj/9dJx44onpZyQBt/JnxW/HcNOmTbjwwgsxceJE3HXXXe224+W56Ih3fV3959ChQ2LdunVi3bp1AoB44IEHxLp168Tnn38uhBDi+eefF0uXLhUVFRViwYIF4uSTTxZXXHFF9PdNTU3itNNOE9/97nfF+++/L7Zu3Sr++te/ilAoJF5//XUhhBArV64UDz74oCgtLRUVFRXi6aefFr179xZXX321L/IohBBz5swRq1atElu3bhX/+Mc/RH5+vpgyZUrcMmPGjBHf+MY3xPvvvy/effdd8bWvfU1MmDAhMHn0+jjecMMNIi8vTyxbtkxUVlZG/44cORJd5pe//KU46aSTxJIlS8SaNWvE8OHDxfDhw6Pft7S0iMGDB4tRo0aJ0tJSUVJSInr37i2mTp0aXeazzz4TXbt2FTfffLPYvHmzmDVrlujQoYMoKSkJTB4ffPBBsWDBArFlyxaxceNG8Zvf/EZkZWWJRYsW+SKPQgixZcsWsW7dOnH99deLr3/969Hy3/p2S01NjejTp4/46U9/KsrKysSzzz4runbtKh5//PFA5G/evHnimWeeEZs3bxabN28Wd911l8jKyhJz5szRmj9Vedy4caPo3bu3+MlPfhK3jr1790aX8fJcdILBh4TWV+7a/k2cOFEIIcRDDz0k+vfvLzp27ChOOukkcdttt8W9uiaEEJ9++qm44oorREFBgejatas466yz4l69Xbt2rRg2bJjIy8sTnTt3FoMGDRJ/+ctfRENDg2/yeOutt4o+ffqIjh07iq997Wvi/vvvF5FIJG6Z/fv3iwkTJohu3bqJ3Nxc8bOf/UwcOnQoMHn0+jha5Q+AmDt3bnSZo0ePil/96lfixBNPFF27dhWXX365qKysjFvP9u3bxdixY0WXLl1Er169xO9+9zvR3Nwct8zSpUvF2WefLTp16iROOeWUuG3o5FYe77nnHnHqqaeKzp07i/z8fDFixAixZMkSX+XxP//zPy3Xs23btugy69evF+eff77IyckRX/nKV8Tdd98dmPzNmzdPDBo0SHTt2lXk5uaKb33rW3GvtpqexxkzZliu4+STT47bllfnohMhIWJ6QxIRERFpxj4fRERE5CoGH0REROQqBh9ERETkKgYfRERE5CoGH0REROQqBh9ERETkKgYfRERE5CoGH0REROQqBh9ERETkKgYfRCTlmmuuQSgUQigUQseOHdGnTx98//vfx5w5cywnSBw9ejQ6dOiADz/8EACwffv26O8T/c2bNw/Lli1L+H1VVZXb2SYihRh8EJG0MWPGoLKyEtu3b8ebb76JCy+8EL/5zW9wySWXoKWlJbrcjh07sHLlSkyePBlz5swBABQVFaGysjL697vf/Q5nnHFG3GdXXXVVdB3l5eVx31VWVho1wzMRycv2OgFE5D85OTkoLCwEAHzlK1/BOeecg/POOw8XXXQR5s2bh5///OcAgLlz5+KSSy7BDTfcgPPOOw8PPPAAunTpEv0tAHTr1g3Z2dlxn8UqKChAjx49tOeJiNzDlg8iUuJ73/sehgwZgpdffhkAIITA3Llz8ZOf/AQDBw7EaaedhhdffNHjVBKRCRh8EJEyAwcOxPbt2wEAixYtwpEjRzB69GgAwE9+8hM88cQT0uvs378/unXrFv0744wzVCaZiDzAxy5EpIwQAqFQCAAwZ84cXHXVVcjOPlbNTJgwATfffDMqKipw6qmn2l7nO++8g+7du0f/37FjR7WJJiLXMfggImU2b96MAQMG4MCBA5g/fz6am5sxe/bs6PfhcBhz5szBXXfdZXudAwYMYJ8PooDhYxciUmLJkiXYuHEjxo8fj3/+85/o378/1q9fj9LS0ujf/fffj3nz5iEcDnudXCLyEFs+iEhaY2MjqqqqEA6HUV1djZKSEhQXF+OSSy7B1VdfjaFDh+K//uu/MHjw4LjfFRUVYerUqSgpKcHFF19sa1t79+5FQ0ND3Gc9e/bk4xciH2PwQUTSSkpK0LdvX2RnZ+PEE0/EkCFD8PDDD2PixIlYt24d1q9fj7///e/tfpeXl4eLLroITzzxhO3g4/TTT2/32apVq3DeeeelnQ8i8kZICCG8TgQRERFlDvb5ICIiIlcx+CAiIiJXMfggIiIiVzH4ICIiIlcx+CAiIiJXMfggIiIiVzH4ICIiIlcx+CAiIiJXMfggIiIiVzH4ICIiIlcx+CAiIiJX/f9Dr2NOj5BzKQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "weather[\"snwd\"].plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 206,
   "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>station</th>\n",
       "      <th>name</th>\n",
       "      <th>prcp</th>\n",
       "      <th>snow</th>\n",
       "      <th>snwd</th>\n",
       "      <th>tmax</th>\n",
       "      <th>tmin</th>\n",
       "      <th>target</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1980-04-01</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>52.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>49.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-02</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>37.0</td>\n",
       "      <td>58.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-03</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>56.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-04</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.57</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>56.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>58.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-05</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>61.0</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-20</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.67</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>66.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-21</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.12</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>61.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-22</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>47.0</td>\n",
       "      <td>62.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-23</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-24</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>15912 rows × 8 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                              name  prcp  snow  snwd  \\\n",
       "DATE                                                                          \n",
       "1980-04-01  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.04   0.0   0.0   \n",
       "1980-04-02  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-03  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-04  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.57   0.0   0.0   \n",
       "1980-04-05  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "...                 ...                               ...   ...   ...   ...   \n",
       "2023-10-20  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.67   0.0   0.0   \n",
       "2023-10-21  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.12   0.0   0.0   \n",
       "2023-10-22  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-23  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-24  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "\n",
       "            tmax  tmin  target  \n",
       "DATE                            \n",
       "1980-04-01  52.0  33.0    49.0  \n",
       "1980-04-02  49.0  37.0    58.0  \n",
       "1980-04-03  58.0  42.0    56.0  \n",
       "1980-04-04  56.0  45.0    58.0  \n",
       "1980-04-05  58.0  44.0    61.0  \n",
       "...          ...   ...     ...  \n",
       "2023-10-20  64.0  60.0    66.0  \n",
       "2023-10-21  66.0  57.0    61.0  \n",
       "2023-10-22  61.0  47.0    62.0  \n",
       "2023-10-23  62.0  44.0    62.0  \n",
       "2023-10-24  62.0  44.0     NaN  \n",
       "\n",
       "[15912 rows x 8 columns]"
      ]
     },
     "execution_count": 206,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather[\"target\"] = weather.shift(-1)[\"tmax\"]\n",
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 207,
   "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>station</th>\n",
       "      <th>name</th>\n",
       "      <th>prcp</th>\n",
       "      <th>snow</th>\n",
       "      <th>snwd</th>\n",
       "      <th>tmax</th>\n",
       "      <th>tmin</th>\n",
       "      <th>target</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1980-04-01</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>52.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>49.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-02</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>37.0</td>\n",
       "      <td>58.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-03</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>56.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-04</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.57</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>56.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>58.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-05</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>61.0</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-20</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.67</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>66.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-21</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.12</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>61.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-22</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>47.0</td>\n",
       "      <td>62.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-23</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-24</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>15912 rows × 8 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                              name  prcp  snow  snwd  \\\n",
       "DATE                                                                          \n",
       "1980-04-01  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.04   0.0   0.0   \n",
       "1980-04-02  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-03  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-04  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.57   0.0   0.0   \n",
       "1980-04-05  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "...                 ...                               ...   ...   ...   ...   \n",
       "2023-10-20  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.67   0.0   0.0   \n",
       "2023-10-21  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.12   0.0   0.0   \n",
       "2023-10-22  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-23  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-24  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "\n",
       "            tmax  tmin  target  \n",
       "DATE                            \n",
       "1980-04-01  52.0  33.0    49.0  \n",
       "1980-04-02  49.0  37.0    58.0  \n",
       "1980-04-03  58.0  42.0    56.0  \n",
       "1980-04-04  56.0  45.0    58.0  \n",
       "1980-04-05  58.0  44.0    61.0  \n",
       "...          ...   ...     ...  \n",
       "2023-10-20  64.0  60.0    66.0  \n",
       "2023-10-21  66.0  57.0    61.0  \n",
       "2023-10-22  61.0  47.0    62.0  \n",
       "2023-10-23  62.0  44.0    62.0  \n",
       "2023-10-24  62.0  44.0    62.0  \n",
       "\n",
       "[15912 rows x 8 columns]"
      ]
     },
     "execution_count": 207,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather = weather.ffill()\n",
    "weather"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 208,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.linear_model import Ridge\n",
    "\n",
    "rr = Ridge(alpha=.1)\n",
    "\n",
    "# import xgboost as xgb\n",
    "# rr = xgb.XGBRegressor(objective=\"reg:linear\", random_state=42)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 209,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['prcp', 'snow', 'snwd', 'tmax', 'tmin'], dtype='object')"
      ]
     },
     "execution_count": 209,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predictors = weather.columns[~weather.columns.isin([\"target\",\"name\",\"station\"])]\n",
    "predictors"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 210,
   "metadata": {},
   "outputs": [],
   "source": [
    "def backtest(weather, model, predictors, start=3650, step=90):\n",
    "    all_predictions = []\n",
    "    for i in range(start, weather.shape[0], step):\n",
    "        # train on first 10 years\n",
    "        train = weather.iloc[:i,:]\n",
    "        test = weather.iloc[i:(i+step),:]\n",
    "        model.fit(train[predictors], train[\"target\"])\n",
    "        preds = model.predict(test[predictors])\n",
    "        preds = pd.Series(preds, index=test.index)\n",
    "        combined = pd.concat([test[\"target\"], preds], axis=1)\n",
    "        combined.columns = [\"actual\",\"prediction\"]\n",
    "        combined[\"diff\"] = (combined[\"prediction\"] - combined[\"actual\"]).abs()\n",
    "        all_predictions.append(combined)\n",
    "    return pd.concat(all_predictions)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 211,
   "metadata": {},
   "outputs": [],
   "source": [
    "predictions = backtest(weather, rr, predictors)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 212,
   "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>actual</th>\n",
       "      <th>prediction</th>\n",
       "      <th>diff</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1990-03-30</th>\n",
       "      <td>45.0</td>\n",
       "      <td>45.015728</td>\n",
       "      <td>0.015728</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-31</th>\n",
       "      <td>44.0</td>\n",
       "      <td>49.328511</td>\n",
       "      <td>5.328511</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-04-01</th>\n",
       "      <td>48.0</td>\n",
       "      <td>50.468676</td>\n",
       "      <td>2.468676</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-04-02</th>\n",
       "      <td>46.0</td>\n",
       "      <td>52.745045</td>\n",
       "      <td>6.745045</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-04-03</th>\n",
       "      <td>48.0</td>\n",
       "      <td>50.790173</td>\n",
       "      <td>2.790173</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-20</th>\n",
       "      <td>66.0</td>\n",
       "      <td>68.625909</td>\n",
       "      <td>2.625909</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-21</th>\n",
       "      <td>61.0</td>\n",
       "      <td>68.594439</td>\n",
       "      <td>7.594439</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-22</th>\n",
       "      <td>62.0</td>\n",
       "      <td>61.379808</td>\n",
       "      <td>0.620192</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-23</th>\n",
       "      <td>62.0</td>\n",
       "      <td>60.312942</td>\n",
       "      <td>1.687058</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-24</th>\n",
       "      <td>62.0</td>\n",
       "      <td>60.312942</td>\n",
       "      <td>1.687058</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>12262 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            actual  prediction      diff\n",
       "DATE                                    \n",
       "1990-03-30    45.0   45.015728  0.015728\n",
       "1990-03-31    44.0   49.328511  5.328511\n",
       "1990-04-01    48.0   50.468676  2.468676\n",
       "1990-04-02    46.0   52.745045  6.745045\n",
       "1990-04-03    48.0   50.790173  2.790173\n",
       "...            ...         ...       ...\n",
       "2023-10-20    66.0   68.625909  2.625909\n",
       "2023-10-21    61.0   68.594439  7.594439\n",
       "2023-10-22    62.0   61.379808  0.620192\n",
       "2023-10-23    62.0   60.312942  1.687058\n",
       "2023-10-24    62.0   60.312942  1.687058\n",
       "\n",
       "[12262 rows x 3 columns]"
      ]
     },
     "execution_count": 212,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predictions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 213,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "5.1088323691185"
      ]
     },
     "execution_count": 213,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from sklearn.metrics import mean_absolute_error\n",
    "mean_absolute_error(predictions[\"actual\"],predictions[\"prediction\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 214,
   "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>station</th>\n",
       "      <th>name</th>\n",
       "      <th>prcp</th>\n",
       "      <th>snow</th>\n",
       "      <th>snwd</th>\n",
       "      <th>tmax</th>\n",
       "      <th>tmin</th>\n",
       "      <th>target</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1980-04-01</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>52.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>49.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-02</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>37.0</td>\n",
       "      <td>58.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-03</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>56.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-04</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.57</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>56.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>58.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-05</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>61.0</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-20</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.67</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>66.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-21</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.12</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>61.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-22</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>47.0</td>\n",
       "      <td>62.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-23</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-24</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>15912 rows × 8 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                              name  prcp  snow  snwd  \\\n",
       "DATE                                                                          \n",
       "1980-04-01  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.04   0.0   0.0   \n",
       "1980-04-02  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-03  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-04  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.57   0.0   0.0   \n",
       "1980-04-05  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "...                 ...                               ...   ...   ...   ...   \n",
       "2023-10-20  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.67   0.0   0.0   \n",
       "2023-10-21  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.12   0.0   0.0   \n",
       "2023-10-22  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-23  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-24  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "\n",
       "            tmax  tmin  target  \n",
       "DATE                            \n",
       "1980-04-01  52.0  33.0    49.0  \n",
       "1980-04-02  49.0  37.0    58.0  \n",
       "1980-04-03  58.0  42.0    56.0  \n",
       "1980-04-04  56.0  45.0    58.0  \n",
       "1980-04-05  58.0  44.0    61.0  \n",
       "...          ...   ...     ...  \n",
       "2023-10-20  64.0  60.0    66.0  \n",
       "2023-10-21  66.0  57.0    61.0  \n",
       "2023-10-22  61.0  47.0    62.0  \n",
       "2023-10-23  62.0  44.0    62.0  \n",
       "2023-10-24  62.0  44.0    62.0  \n",
       "\n",
       "[15912 rows x 8 columns]"
      ]
     },
     "execution_count": 214,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 215,
   "metadata": {},
   "outputs": [],
   "source": [
    "def pct_diff(new, old):\n",
    "    return (new - old)/old\n",
    "\n",
    "def compute_rolling(weather, horizon, col):\n",
    "    label = f\"rolling_{horizon}_{col}\"\n",
    "    weather[label] = weather[col].rolling(horizon).mean()\n",
    "    weather[f\"{label}_pct\"] = pct_diff(weather[col], weather[label])\n",
    "    return weather\n",
    "# days for rolling mean\n",
    "rolling_horizons = [3,14]\n",
    "\n",
    "for horizon in rolling_horizons:\n",
    "    for col in [\"tmax\",\"tmin\",\"prcp\"]:\n",
    "        weather = compute_rolling(weather, horizon, col)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 216,
   "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>station</th>\n",
       "      <th>name</th>\n",
       "      <th>prcp</th>\n",
       "      <th>snow</th>\n",
       "      <th>snwd</th>\n",
       "      <th>tmax</th>\n",
       "      <th>tmin</th>\n",
       "      <th>target</th>\n",
       "      <th>rolling_3_tmax</th>\n",
       "      <th>rolling_3_tmax_pct</th>\n",
       "      <th>rolling_3_tmin</th>\n",
       "      <th>rolling_3_tmin_pct</th>\n",
       "      <th>rolling_3_prcp</th>\n",
       "      <th>rolling_3_prcp_pct</th>\n",
       "      <th>rolling_14_tmax</th>\n",
       "      <th>rolling_14_tmax_pct</th>\n",
       "      <th>rolling_14_tmin</th>\n",
       "      <th>rolling_14_tmin_pct</th>\n",
       "      <th>rolling_14_prcp</th>\n",
       "      <th>rolling_14_prcp_pct</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1980-04-01</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>52.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-02</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>37.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-03</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>56.0</td>\n",
       "      <td>53.000000</td>\n",
       "      <td>0.094340</td>\n",
       "      <td>37.333333</td>\n",
       "      <td>0.125000</td>\n",
       "      <td>0.013333</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-04</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.57</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>56.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>54.333333</td>\n",
       "      <td>0.030675</td>\n",
       "      <td>41.333333</td>\n",
       "      <td>0.088710</td>\n",
       "      <td>0.190000</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-05</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>57.333333</td>\n",
       "      <td>0.011628</td>\n",
       "      <td>43.666667</td>\n",
       "      <td>0.007634</td>\n",
       "      <td>0.190000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-20</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.67</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>64.333333</td>\n",
       "      <td>-0.005181</td>\n",
       "      <td>53.333333</td>\n",
       "      <td>0.125000</td>\n",
       "      <td>0.223333</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>64.071429</td>\n",
       "      <td>-0.001115</td>\n",
       "      <td>50.857143</td>\n",
       "      <td>0.179775</td>\n",
       "      <td>0.240000</td>\n",
       "      <td>1.791667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-21</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.12</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>64.666667</td>\n",
       "      <td>0.020619</td>\n",
       "      <td>55.666667</td>\n",
       "      <td>0.023952</td>\n",
       "      <td>0.263333</td>\n",
       "      <td>-0.544304</td>\n",
       "      <td>64.000000</td>\n",
       "      <td>0.031250</td>\n",
       "      <td>51.071429</td>\n",
       "      <td>0.116084</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>0.050000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-22</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>47.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>63.666667</td>\n",
       "      <td>-0.041885</td>\n",
       "      <td>54.666667</td>\n",
       "      <td>-0.140244</td>\n",
       "      <td>0.263333</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>64.000000</td>\n",
       "      <td>-0.046875</td>\n",
       "      <td>50.857143</td>\n",
       "      <td>-0.075843</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-23</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>63.000000</td>\n",
       "      <td>-0.015873</td>\n",
       "      <td>49.333333</td>\n",
       "      <td>-0.108108</td>\n",
       "      <td>0.040000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>64.000000</td>\n",
       "      <td>-0.031250</td>\n",
       "      <td>50.500000</td>\n",
       "      <td>-0.128713</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-24</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>61.666667</td>\n",
       "      <td>0.005405</td>\n",
       "      <td>45.000000</td>\n",
       "      <td>-0.022222</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>63.857143</td>\n",
       "      <td>-0.029083</td>\n",
       "      <td>50.000000</td>\n",
       "      <td>-0.120000</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>15912 rows × 20 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                              name  prcp  snow  snwd  \\\n",
       "DATE                                                                          \n",
       "1980-04-01  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.04   0.0   0.0   \n",
       "1980-04-02  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-03  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-04  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.57   0.0   0.0   \n",
       "1980-04-05  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "...                 ...                               ...   ...   ...   ...   \n",
       "2023-10-20  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.67   0.0   0.0   \n",
       "2023-10-21  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.12   0.0   0.0   \n",
       "2023-10-22  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-23  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-24  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "\n",
       "            tmax  tmin  target  rolling_3_tmax  rolling_3_tmax_pct  \\\n",
       "DATE                                                                 \n",
       "1980-04-01  52.0  33.0    49.0             NaN                 NaN   \n",
       "1980-04-02  49.0  37.0    58.0             NaN                 NaN   \n",
       "1980-04-03  58.0  42.0    56.0       53.000000            0.094340   \n",
       "1980-04-04  56.0  45.0    58.0       54.333333            0.030675   \n",
       "1980-04-05  58.0  44.0    61.0       57.333333            0.011628   \n",
       "...          ...   ...     ...             ...                 ...   \n",
       "2023-10-20  64.0  60.0    66.0       64.333333           -0.005181   \n",
       "2023-10-21  66.0  57.0    61.0       64.666667            0.020619   \n",
       "2023-10-22  61.0  47.0    62.0       63.666667           -0.041885   \n",
       "2023-10-23  62.0  44.0    62.0       63.000000           -0.015873   \n",
       "2023-10-24  62.0  44.0    62.0       61.666667            0.005405   \n",
       "\n",
       "            rolling_3_tmin  rolling_3_tmin_pct  rolling_3_prcp  \\\n",
       "DATE                                                             \n",
       "1980-04-01             NaN                 NaN             NaN   \n",
       "1980-04-02             NaN                 NaN             NaN   \n",
       "1980-04-03       37.333333            0.125000        0.013333   \n",
       "1980-04-04       41.333333            0.088710        0.190000   \n",
       "1980-04-05       43.666667            0.007634        0.190000   \n",
       "...                    ...                 ...             ...   \n",
       "2023-10-20       53.333333            0.125000        0.223333   \n",
       "2023-10-21       55.666667            0.023952        0.263333   \n",
       "2023-10-22       54.666667           -0.140244        0.263333   \n",
       "2023-10-23       49.333333           -0.108108        0.040000   \n",
       "2023-10-24       45.000000           -0.022222        0.000000   \n",
       "\n",
       "            rolling_3_prcp_pct  rolling_14_tmax  rolling_14_tmax_pct  \\\n",
       "DATE                                                                   \n",
       "1980-04-01                 NaN              NaN                  NaN   \n",
       "1980-04-02                 NaN              NaN                  NaN   \n",
       "1980-04-03           -1.000000              NaN                  NaN   \n",
       "1980-04-04            2.000000              NaN                  NaN   \n",
       "1980-04-05           -1.000000              NaN                  NaN   \n",
       "...                        ...              ...                  ...   \n",
       "2023-10-20            2.000000        64.071429            -0.001115   \n",
       "2023-10-21           -0.544304        64.000000             0.031250   \n",
       "2023-10-22           -1.000000        64.000000            -0.046875   \n",
       "2023-10-23           -1.000000        64.000000            -0.031250   \n",
       "2023-10-24                 NaN        63.857143            -0.029083   \n",
       "\n",
       "            rolling_14_tmin  rolling_14_tmin_pct  rolling_14_prcp  \\\n",
       "DATE                                                                \n",
       "1980-04-01              NaN                  NaN              NaN   \n",
       "1980-04-02              NaN                  NaN              NaN   \n",
       "1980-04-03              NaN                  NaN              NaN   \n",
       "1980-04-04              NaN                  NaN              NaN   \n",
       "1980-04-05              NaN                  NaN              NaN   \n",
       "...                     ...                  ...              ...   \n",
       "2023-10-20        50.857143             0.179775         0.240000   \n",
       "2023-10-21        51.071429             0.116084         0.114286   \n",
       "2023-10-22        50.857143            -0.075843         0.114286   \n",
       "2023-10-23        50.500000            -0.128713         0.114286   \n",
       "2023-10-24        50.000000            -0.120000         0.114286   \n",
       "\n",
       "            rolling_14_prcp_pct  \n",
       "DATE                             \n",
       "1980-04-01                  NaN  \n",
       "1980-04-02                  NaN  \n",
       "1980-04-03                  NaN  \n",
       "1980-04-04                  NaN  \n",
       "1980-04-05                  NaN  \n",
       "...                         ...  \n",
       "2023-10-20             1.791667  \n",
       "2023-10-21             0.050000  \n",
       "2023-10-22            -1.000000  \n",
       "2023-10-23            -1.000000  \n",
       "2023-10-24            -1.000000  \n",
       "\n",
       "[15912 rows x 20 columns]"
      ]
     },
     "execution_count": 216,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 217,
   "metadata": {},
   "outputs": [],
   "source": [
    "weather = weather.iloc[14:,:]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 218,
   "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>station</th>\n",
       "      <th>name</th>\n",
       "      <th>prcp</th>\n",
       "      <th>snow</th>\n",
       "      <th>snwd</th>\n",
       "      <th>tmax</th>\n",
       "      <th>tmin</th>\n",
       "      <th>target</th>\n",
       "      <th>rolling_3_tmax</th>\n",
       "      <th>rolling_3_tmax_pct</th>\n",
       "      <th>rolling_3_tmin</th>\n",
       "      <th>rolling_3_tmin_pct</th>\n",
       "      <th>rolling_3_prcp</th>\n",
       "      <th>rolling_3_prcp_pct</th>\n",
       "      <th>rolling_14_tmax</th>\n",
       "      <th>rolling_14_tmax_pct</th>\n",
       "      <th>rolling_14_tmin</th>\n",
       "      <th>rolling_14_tmin_pct</th>\n",
       "      <th>rolling_14_prcp</th>\n",
       "      <th>rolling_14_prcp_pct</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1980-04-15</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>63.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>61.333333</td>\n",
       "      <td>0.027174</td>\n",
       "      <td>48.666667</td>\n",
       "      <td>0.006849</td>\n",
       "      <td>0.186667</td>\n",
       "      <td>-0.892857</td>\n",
       "      <td>58.285714</td>\n",
       "      <td>0.080882</td>\n",
       "      <td>46.428571</td>\n",
       "      <td>0.055385</td>\n",
       "      <td>0.319286</td>\n",
       "      <td>-0.937360</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-16</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>35.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>55.666667</td>\n",
       "      <td>-0.119760</td>\n",
       "      <td>43.000000</td>\n",
       "      <td>-0.186047</td>\n",
       "      <td>0.186667</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>58.285714</td>\n",
       "      <td>-0.159314</td>\n",
       "      <td>46.285714</td>\n",
       "      <td>-0.243827</td>\n",
       "      <td>0.319286</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-17</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>55.000000</td>\n",
       "      <td>-0.036364</td>\n",
       "      <td>39.000000</td>\n",
       "      <td>-0.153846</td>\n",
       "      <td>0.006667</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>57.928571</td>\n",
       "      <td>-0.085080</td>\n",
       "      <td>45.642857</td>\n",
       "      <td>-0.276995</td>\n",
       "      <td>0.319286</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-18</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>53.666667</td>\n",
       "      <td>0.099379</td>\n",
       "      <td>37.666667</td>\n",
       "      <td>0.194690</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>58.142857</td>\n",
       "      <td>0.014742</td>\n",
       "      <td>45.642857</td>\n",
       "      <td>-0.014085</td>\n",
       "      <td>0.278571</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-19</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>56.666667</td>\n",
       "      <td>0.023529</td>\n",
       "      <td>41.000000</td>\n",
       "      <td>0.097561</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>58.142857</td>\n",
       "      <td>-0.002457</td>\n",
       "      <td>45.714286</td>\n",
       "      <td>-0.015625</td>\n",
       "      <td>0.278571</td>\n",
       "      <td>-1.000000</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-20</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.67</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>64.333333</td>\n",
       "      <td>-0.005181</td>\n",
       "      <td>53.333333</td>\n",
       "      <td>0.125000</td>\n",
       "      <td>0.223333</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>64.071429</td>\n",
       "      <td>-0.001115</td>\n",
       "      <td>50.857143</td>\n",
       "      <td>0.179775</td>\n",
       "      <td>0.240000</td>\n",
       "      <td>1.791667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-21</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.12</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>64.666667</td>\n",
       "      <td>0.020619</td>\n",
       "      <td>55.666667</td>\n",
       "      <td>0.023952</td>\n",
       "      <td>0.263333</td>\n",
       "      <td>-0.544304</td>\n",
       "      <td>64.000000</td>\n",
       "      <td>0.031250</td>\n",
       "      <td>51.071429</td>\n",
       "      <td>0.116084</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>0.050000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-22</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>47.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>63.666667</td>\n",
       "      <td>-0.041885</td>\n",
       "      <td>54.666667</td>\n",
       "      <td>-0.140244</td>\n",
       "      <td>0.263333</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>64.000000</td>\n",
       "      <td>-0.046875</td>\n",
       "      <td>50.857143</td>\n",
       "      <td>-0.075843</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-23</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>63.000000</td>\n",
       "      <td>-0.015873</td>\n",
       "      <td>49.333333</td>\n",
       "      <td>-0.108108</td>\n",
       "      <td>0.040000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>64.000000</td>\n",
       "      <td>-0.031250</td>\n",
       "      <td>50.500000</td>\n",
       "      <td>-0.128713</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-24</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>61.666667</td>\n",
       "      <td>0.005405</td>\n",
       "      <td>45.000000</td>\n",
       "      <td>-0.022222</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>63.857143</td>\n",
       "      <td>-0.029083</td>\n",
       "      <td>50.000000</td>\n",
       "      <td>-0.120000</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>15898 rows × 20 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                              name  prcp  snow  snwd  \\\n",
       "DATE                                                                          \n",
       "1980-04-15  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.02   0.0   0.0   \n",
       "1980-04-16  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-17  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-18  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-19  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "...                 ...                               ...   ...   ...   ...   \n",
       "2023-10-20  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.67   0.0   0.0   \n",
       "2023-10-21  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.12   0.0   0.0   \n",
       "2023-10-22  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-23  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-24  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "\n",
       "            tmax  tmin  target  rolling_3_tmax  rolling_3_tmax_pct  \\\n",
       "DATE                                                                 \n",
       "1980-04-15  63.0  49.0    49.0       61.333333            0.027174   \n",
       "1980-04-16  49.0  35.0    53.0       55.666667           -0.119760   \n",
       "1980-04-17  53.0  33.0    59.0       55.000000           -0.036364   \n",
       "1980-04-18  59.0  45.0    58.0       53.666667            0.099379   \n",
       "1980-04-19  58.0  45.0    60.0       56.666667            0.023529   \n",
       "...          ...   ...     ...             ...                 ...   \n",
       "2023-10-20  64.0  60.0    66.0       64.333333           -0.005181   \n",
       "2023-10-21  66.0  57.0    61.0       64.666667            0.020619   \n",
       "2023-10-22  61.0  47.0    62.0       63.666667           -0.041885   \n",
       "2023-10-23  62.0  44.0    62.0       63.000000           -0.015873   \n",
       "2023-10-24  62.0  44.0    62.0       61.666667            0.005405   \n",
       "\n",
       "            rolling_3_tmin  rolling_3_tmin_pct  rolling_3_prcp  \\\n",
       "DATE                                                             \n",
       "1980-04-15       48.666667            0.006849        0.186667   \n",
       "1980-04-16       43.000000           -0.186047        0.186667   \n",
       "1980-04-17       39.000000           -0.153846        0.006667   \n",
       "1980-04-18       37.666667            0.194690        0.000000   \n",
       "1980-04-19       41.000000            0.097561        0.000000   \n",
       "...                    ...                 ...             ...   \n",
       "2023-10-20       53.333333            0.125000        0.223333   \n",
       "2023-10-21       55.666667            0.023952        0.263333   \n",
       "2023-10-22       54.666667           -0.140244        0.263333   \n",
       "2023-10-23       49.333333           -0.108108        0.040000   \n",
       "2023-10-24       45.000000           -0.022222        0.000000   \n",
       "\n",
       "            rolling_3_prcp_pct  rolling_14_tmax  rolling_14_tmax_pct  \\\n",
       "DATE                                                                   \n",
       "1980-04-15           -0.892857        58.285714             0.080882   \n",
       "1980-04-16           -1.000000        58.285714            -0.159314   \n",
       "1980-04-17           -1.000000        57.928571            -0.085080   \n",
       "1980-04-18                 NaN        58.142857             0.014742   \n",
       "1980-04-19                 NaN        58.142857            -0.002457   \n",
       "...                        ...              ...                  ...   \n",
       "2023-10-20            2.000000        64.071429            -0.001115   \n",
       "2023-10-21           -0.544304        64.000000             0.031250   \n",
       "2023-10-22           -1.000000        64.000000            -0.046875   \n",
       "2023-10-23           -1.000000        64.000000            -0.031250   \n",
       "2023-10-24                 NaN        63.857143            -0.029083   \n",
       "\n",
       "            rolling_14_tmin  rolling_14_tmin_pct  rolling_14_prcp  \\\n",
       "DATE                                                                \n",
       "1980-04-15        46.428571             0.055385         0.319286   \n",
       "1980-04-16        46.285714            -0.243827         0.319286   \n",
       "1980-04-17        45.642857            -0.276995         0.319286   \n",
       "1980-04-18        45.642857            -0.014085         0.278571   \n",
       "1980-04-19        45.714286            -0.015625         0.278571   \n",
       "...                     ...                  ...              ...   \n",
       "2023-10-20        50.857143             0.179775         0.240000   \n",
       "2023-10-21        51.071429             0.116084         0.114286   \n",
       "2023-10-22        50.857143            -0.075843         0.114286   \n",
       "2023-10-23        50.500000            -0.128713         0.114286   \n",
       "2023-10-24        50.000000            -0.120000         0.114286   \n",
       "\n",
       "            rolling_14_prcp_pct  \n",
       "DATE                             \n",
       "1980-04-15            -0.937360  \n",
       "1980-04-16            -1.000000  \n",
       "1980-04-17            -1.000000  \n",
       "1980-04-18            -1.000000  \n",
       "1980-04-19            -1.000000  \n",
       "...                         ...  \n",
       "2023-10-20             1.791667  \n",
       "2023-10-21             0.050000  \n",
       "2023-10-22            -1.000000  \n",
       "2023-10-23            -1.000000  \n",
       "2023-10-24            -1.000000  \n",
       "\n",
       "[15898 rows x 20 columns]"
      ]
     },
     "execution_count": 218,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 219,
   "metadata": {},
   "outputs": [],
   "source": [
    "weather = weather.fillna(0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 220,
   "metadata": {},
   "outputs": [],
   "source": [
    "def expandMean(dframe):\n",
    "    return dframe.expanding(1).mean()\n",
    "for col in [\"tmax\",\"tmin\",\"prcp\"]:\n",
    "    weather[f\"month_avg_{col}\"] = weather[col].groupby(weather.index.month, group_keys=False).apply(expandMean)\n",
    "    weather[f\"day_avg_{col}\"] = weather[col].groupby(weather.index.day_of_year, group_keys=False).apply(expandMean)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 221,
   "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>station</th>\n",
       "      <th>name</th>\n",
       "      <th>prcp</th>\n",
       "      <th>snow</th>\n",
       "      <th>snwd</th>\n",
       "      <th>tmax</th>\n",
       "      <th>tmin</th>\n",
       "      <th>target</th>\n",
       "      <th>rolling_3_tmax</th>\n",
       "      <th>rolling_3_tmax_pct</th>\n",
       "      <th>...</th>\n",
       "      <th>rolling_14_tmin</th>\n",
       "      <th>rolling_14_tmin_pct</th>\n",
       "      <th>rolling_14_prcp</th>\n",
       "      <th>rolling_14_prcp_pct</th>\n",
       "      <th>month_avg_tmax</th>\n",
       "      <th>day_avg_tmax</th>\n",
       "      <th>month_avg_tmin</th>\n",
       "      <th>day_avg_tmin</th>\n",
       "      <th>month_avg_prcp</th>\n",
       "      <th>day_avg_prcp</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1980-04-15</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>63.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>61.333333</td>\n",
       "      <td>0.027174</td>\n",
       "      <td>...</td>\n",
       "      <td>46.428571</td>\n",
       "      <td>0.055385</td>\n",
       "      <td>0.319286</td>\n",
       "      <td>-0.937360</td>\n",
       "      <td>63.000000</td>\n",
       "      <td>63.000000</td>\n",
       "      <td>49.000000</td>\n",
       "      <td>49.000000</td>\n",
       "      <td>0.020000</td>\n",
       "      <td>0.020000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-16</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>35.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>55.666667</td>\n",
       "      <td>-0.119760</td>\n",
       "      <td>...</td>\n",
       "      <td>46.285714</td>\n",
       "      <td>-0.243827</td>\n",
       "      <td>0.319286</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>56.000000</td>\n",
       "      <td>49.000000</td>\n",
       "      <td>42.000000</td>\n",
       "      <td>35.000000</td>\n",
       "      <td>0.010000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-17</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>55.000000</td>\n",
       "      <td>-0.036364</td>\n",
       "      <td>...</td>\n",
       "      <td>45.642857</td>\n",
       "      <td>-0.276995</td>\n",
       "      <td>0.319286</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>55.000000</td>\n",
       "      <td>53.000000</td>\n",
       "      <td>39.000000</td>\n",
       "      <td>33.000000</td>\n",
       "      <td>0.006667</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-18</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>53.666667</td>\n",
       "      <td>0.099379</td>\n",
       "      <td>...</td>\n",
       "      <td>45.642857</td>\n",
       "      <td>-0.014085</td>\n",
       "      <td>0.278571</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>56.000000</td>\n",
       "      <td>59.000000</td>\n",
       "      <td>40.500000</td>\n",
       "      <td>45.000000</td>\n",
       "      <td>0.005000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-04-19</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>58.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>56.666667</td>\n",
       "      <td>0.023529</td>\n",
       "      <td>...</td>\n",
       "      <td>45.714286</td>\n",
       "      <td>-0.015625</td>\n",
       "      <td>0.278571</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>56.400000</td>\n",
       "      <td>58.000000</td>\n",
       "      <td>41.400000</td>\n",
       "      <td>45.000000</td>\n",
       "      <td>0.004000</td>\n",
       "      <td>0.000000</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-20</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.67</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>64.333333</td>\n",
       "      <td>-0.005181</td>\n",
       "      <td>...</td>\n",
       "      <td>50.857143</td>\n",
       "      <td>0.179775</td>\n",
       "      <td>0.240000</td>\n",
       "      <td>1.791667</td>\n",
       "      <td>65.047302</td>\n",
       "      <td>64.090909</td>\n",
       "      <td>50.463415</td>\n",
       "      <td>49.545455</td>\n",
       "      <td>0.118330</td>\n",
       "      <td>0.215227</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-21</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.12</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>64.666667</td>\n",
       "      <td>0.020619</td>\n",
       "      <td>...</td>\n",
       "      <td>51.071429</td>\n",
       "      <td>0.116084</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>0.050000</td>\n",
       "      <td>65.048006</td>\n",
       "      <td>63.750000</td>\n",
       "      <td>50.468242</td>\n",
       "      <td>48.340909</td>\n",
       "      <td>0.118331</td>\n",
       "      <td>0.056818</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-22</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>47.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>63.666667</td>\n",
       "      <td>-0.041885</td>\n",
       "      <td>...</td>\n",
       "      <td>50.857143</td>\n",
       "      <td>-0.075843</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>65.045018</td>\n",
       "      <td>62.954545</td>\n",
       "      <td>50.465683</td>\n",
       "      <td>48.431818</td>\n",
       "      <td>0.118244</td>\n",
       "      <td>0.116364</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-23</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>63.000000</td>\n",
       "      <td>-0.015873</td>\n",
       "      <td>...</td>\n",
       "      <td>50.500000</td>\n",
       "      <td>-0.128713</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>65.042773</td>\n",
       "      <td>62.068182</td>\n",
       "      <td>50.460914</td>\n",
       "      <td>48.272727</td>\n",
       "      <td>0.118156</td>\n",
       "      <td>0.146818</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2023-10-24</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>61.666667</td>\n",
       "      <td>0.005405</td>\n",
       "      <td>...</td>\n",
       "      <td>50.000000</td>\n",
       "      <td>-0.120000</td>\n",
       "      <td>0.114286</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>65.040531</td>\n",
       "      <td>61.954545</td>\n",
       "      <td>50.456153</td>\n",
       "      <td>47.386364</td>\n",
       "      <td>0.118069</td>\n",
       "      <td>0.099318</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>15898 rows × 26 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                              name  prcp  snow  snwd  \\\n",
       "DATE                                                                          \n",
       "1980-04-15  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.02   0.0   0.0   \n",
       "1980-04-16  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-17  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-18  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "1980-04-19  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "...                 ...                               ...   ...   ...   ...   \n",
       "2023-10-20  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.67   0.0   0.0   \n",
       "2023-10-21  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.12   0.0   0.0   \n",
       "2023-10-22  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-23  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2023-10-24  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "\n",
       "            tmax  tmin  target  rolling_3_tmax  rolling_3_tmax_pct  ...  \\\n",
       "DATE                                                                ...   \n",
       "1980-04-15  63.0  49.0    49.0       61.333333            0.027174  ...   \n",
       "1980-04-16  49.0  35.0    53.0       55.666667           -0.119760  ...   \n",
       "1980-04-17  53.0  33.0    59.0       55.000000           -0.036364  ...   \n",
       "1980-04-18  59.0  45.0    58.0       53.666667            0.099379  ...   \n",
       "1980-04-19  58.0  45.0    60.0       56.666667            0.023529  ...   \n",
       "...          ...   ...     ...             ...                 ...  ...   \n",
       "2023-10-20  64.0  60.0    66.0       64.333333           -0.005181  ...   \n",
       "2023-10-21  66.0  57.0    61.0       64.666667            0.020619  ...   \n",
       "2023-10-22  61.0  47.0    62.0       63.666667           -0.041885  ...   \n",
       "2023-10-23  62.0  44.0    62.0       63.000000           -0.015873  ...   \n",
       "2023-10-24  62.0  44.0    62.0       61.666667            0.005405  ...   \n",
       "\n",
       "            rolling_14_tmin  rolling_14_tmin_pct  rolling_14_prcp  \\\n",
       "DATE                                                                \n",
       "1980-04-15        46.428571             0.055385         0.319286   \n",
       "1980-04-16        46.285714            -0.243827         0.319286   \n",
       "1980-04-17        45.642857            -0.276995         0.319286   \n",
       "1980-04-18        45.642857            -0.014085         0.278571   \n",
       "1980-04-19        45.714286            -0.015625         0.278571   \n",
       "...                     ...                  ...              ...   \n",
       "2023-10-20        50.857143             0.179775         0.240000   \n",
       "2023-10-21        51.071429             0.116084         0.114286   \n",
       "2023-10-22        50.857143            -0.075843         0.114286   \n",
       "2023-10-23        50.500000            -0.128713         0.114286   \n",
       "2023-10-24        50.000000            -0.120000         0.114286   \n",
       "\n",
       "            rolling_14_prcp_pct  month_avg_tmax  day_avg_tmax  month_avg_tmin  \\\n",
       "DATE                                                                            \n",
       "1980-04-15            -0.937360       63.000000     63.000000       49.000000   \n",
       "1980-04-16            -1.000000       56.000000     49.000000       42.000000   \n",
       "1980-04-17            -1.000000       55.000000     53.000000       39.000000   \n",
       "1980-04-18            -1.000000       56.000000     59.000000       40.500000   \n",
       "1980-04-19            -1.000000       56.400000     58.000000       41.400000   \n",
       "...                         ...             ...           ...             ...   \n",
       "2023-10-20             1.791667       65.047302     64.090909       50.463415   \n",
       "2023-10-21             0.050000       65.048006     63.750000       50.468242   \n",
       "2023-10-22            -1.000000       65.045018     62.954545       50.465683   \n",
       "2023-10-23            -1.000000       65.042773     62.068182       50.460914   \n",
       "2023-10-24            -1.000000       65.040531     61.954545       50.456153   \n",
       "\n",
       "            day_avg_tmin  month_avg_prcp  day_avg_prcp  \n",
       "DATE                                                    \n",
       "1980-04-15     49.000000        0.020000      0.020000  \n",
       "1980-04-16     35.000000        0.010000      0.000000  \n",
       "1980-04-17     33.000000        0.006667      0.000000  \n",
       "1980-04-18     45.000000        0.005000      0.000000  \n",
       "1980-04-19     45.000000        0.004000      0.000000  \n",
       "...                  ...             ...           ...  \n",
       "2023-10-20     49.545455        0.118330      0.215227  \n",
       "2023-10-21     48.340909        0.118331      0.056818  \n",
       "2023-10-22     48.431818        0.118244      0.116364  \n",
       "2023-10-23     48.272727        0.118156      0.146818  \n",
       "2023-10-24     47.386364        0.118069      0.099318  \n",
       "\n",
       "[15898 rows x 26 columns]"
      ]
     },
     "execution_count": 221,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 222,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['prcp', 'snow', 'snwd', 'tmax', 'tmin', 'rolling_3_tmax',\n",
       "       'rolling_3_tmax_pct', 'rolling_3_tmin', 'rolling_3_tmin_pct',\n",
       "       'rolling_3_prcp', 'rolling_3_prcp_pct', 'rolling_14_tmax',\n",
       "       'rolling_14_tmax_pct', 'rolling_14_tmin', 'rolling_14_tmin_pct',\n",
       "       'rolling_14_prcp', 'rolling_14_prcp_pct', 'month_avg_tmax',\n",
       "       'day_avg_tmax', 'month_avg_tmin', 'day_avg_tmin', 'month_avg_prcp',\n",
       "       'day_avg_prcp'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 222,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predictors = weather.columns[~weather.columns.isin([\"target\",\"name\",\"station\"])]\n",
    "predictors"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 223,
   "metadata": {},
   "outputs": [],
   "source": [
    "predictions = backtest(weather, rr, predictors)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 224,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4.770979960120657"
      ]
     },
     "execution_count": 224,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mean_absolute_error(predictions[\"actual\"],predictions[\"prediction\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 225,
   "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>actual</th>\n",
       "      <th>prediction</th>\n",
       "      <th>diff</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2007-03-26</th>\n",
       "      <td>78.0</td>\n",
       "      <td>50.094404</td>\n",
       "      <td>27.905596</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1998-03-26</th>\n",
       "      <td>80.0</td>\n",
       "      <td>52.300436</td>\n",
       "      <td>27.699564</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2003-04-15</th>\n",
       "      <td>86.0</td>\n",
       "      <td>59.468078</td>\n",
       "      <td>26.531922</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-12-23</th>\n",
       "      <td>16.0</td>\n",
       "      <td>41.459917</td>\n",
       "      <td>25.459917</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1997-02-26</th>\n",
       "      <td>71.0</td>\n",
       "      <td>45.759537</td>\n",
       "      <td>25.240463</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2015-06-07</th>\n",
       "      <td>73.0</td>\n",
       "      <td>73.004960</td>\n",
       "      <td>0.004960</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-02-20</th>\n",
       "      <td>43.0</td>\n",
       "      <td>42.995145</td>\n",
       "      <td>0.004855</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2003-03-12</th>\n",
       "      <td>48.0</td>\n",
       "      <td>47.996195</td>\n",
       "      <td>0.003805</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2020-12-21</th>\n",
       "      <td>45.0</td>\n",
       "      <td>45.000988</td>\n",
       "      <td>0.000988</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-12-02</th>\n",
       "      <td>53.0</td>\n",
       "      <td>52.999734</td>\n",
       "      <td>0.000266</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>12248 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            actual  prediction       diff\n",
       "DATE                                     \n",
       "2007-03-26    78.0   50.094404  27.905596\n",
       "1998-03-26    80.0   52.300436  27.699564\n",
       "2003-04-15    86.0   59.468078  26.531922\n",
       "2022-12-23    16.0   41.459917  25.459917\n",
       "1997-02-26    71.0   45.759537  25.240463\n",
       "...            ...         ...        ...\n",
       "2015-06-07    73.0   73.004960   0.004960\n",
       "2000-02-20    43.0   42.995145   0.004855\n",
       "2003-03-12    48.0   47.996195   0.003805\n",
       "2020-12-21    45.0   45.000988   0.000988\n",
       "1990-12-02    53.0   52.999734   0.000266\n",
       "\n",
       "[12248 rows x 3 columns]"
      ]
     },
     "execution_count": 225,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predictions.sort_values(\"diff\",ascending=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 226,
   "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>station</th>\n",
       "      <th>name</th>\n",
       "      <th>prcp</th>\n",
       "      <th>snow</th>\n",
       "      <th>snwd</th>\n",
       "      <th>tmax</th>\n",
       "      <th>tmin</th>\n",
       "      <th>target</th>\n",
       "      <th>rolling_3_tmax</th>\n",
       "      <th>rolling_3_tmax_pct</th>\n",
       "      <th>...</th>\n",
       "      <th>rolling_14_tmin</th>\n",
       "      <th>rolling_14_tmin_pct</th>\n",
       "      <th>rolling_14_prcp</th>\n",
       "      <th>rolling_14_prcp_pct</th>\n",
       "      <th>month_avg_tmax</th>\n",
       "      <th>day_avg_tmax</th>\n",
       "      <th>month_avg_tmin</th>\n",
       "      <th>day_avg_tmin</th>\n",
       "      <th>month_avg_prcp</th>\n",
       "      <th>day_avg_prcp</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DATE</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2007-03-20</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>52.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>43.666667</td>\n",
       "      <td>0.190840</td>\n",
       "      <td>...</td>\n",
       "      <td>27.928571</td>\n",
       "      <td>0.181586</td>\n",
       "      <td>0.120000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>48.782082</td>\n",
       "      <td>47.740741</td>\n",
       "      <td>34.041162</td>\n",
       "      <td>34.481481</td>\n",
       "      <td>0.119891</td>\n",
       "      <td>0.155185</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-21</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>26.0</td>\n",
       "      <td>54.0</td>\n",
       "      <td>44.000000</td>\n",
       "      <td>-0.090909</td>\n",
       "      <td>...</td>\n",
       "      <td>28.857143</td>\n",
       "      <td>-0.099010</td>\n",
       "      <td>0.116429</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>48.771463</td>\n",
       "      <td>46.666667</td>\n",
       "      <td>34.031439</td>\n",
       "      <td>34.037037</td>\n",
       "      <td>0.119746</td>\n",
       "      <td>0.144815</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-22</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>54.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>48.666667</td>\n",
       "      <td>0.109589</td>\n",
       "      <td>...</td>\n",
       "      <td>30.500000</td>\n",
       "      <td>0.311475</td>\n",
       "      <td>0.119286</td>\n",
       "      <td>-0.664671</td>\n",
       "      <td>48.777778</td>\n",
       "      <td>48.777778</td>\n",
       "      <td>34.038647</td>\n",
       "      <td>34.407407</td>\n",
       "      <td>0.119650</td>\n",
       "      <td>0.088148</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-23</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.14</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>46.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>51.666667</td>\n",
       "      <td>0.180645</td>\n",
       "      <td>...</td>\n",
       "      <td>32.785714</td>\n",
       "      <td>0.403050</td>\n",
       "      <td>0.129286</td>\n",
       "      <td>0.082873</td>\n",
       "      <td>48.792521</td>\n",
       "      <td>49.777778</td>\n",
       "      <td>34.053076</td>\n",
       "      <td>34.037037</td>\n",
       "      <td>0.119674</td>\n",
       "      <td>0.080000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-24</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>41.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>56.000000</td>\n",
       "      <td>-0.053571</td>\n",
       "      <td>...</td>\n",
       "      <td>33.857143</td>\n",
       "      <td>0.210970</td>\n",
       "      <td>0.131429</td>\n",
       "      <td>-0.695652</td>\n",
       "      <td>48.797590</td>\n",
       "      <td>49.259259</td>\n",
       "      <td>34.061446</td>\n",
       "      <td>34.888889</td>\n",
       "      <td>0.119578</td>\n",
       "      <td>0.112593</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-25</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>37.0</td>\n",
       "      <td>48.0</td>\n",
       "      <td>55.666667</td>\n",
       "      <td>-0.047904</td>\n",
       "      <td>...</td>\n",
       "      <td>33.785714</td>\n",
       "      <td>0.095137</td>\n",
       "      <td>0.115000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>48.802647</td>\n",
       "      <td>51.888889</td>\n",
       "      <td>34.064982</td>\n",
       "      <td>35.074074</td>\n",
       "      <td>0.119434</td>\n",
       "      <td>0.001481</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-26</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>48.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>78.0</td>\n",
       "      <td>51.333333</td>\n",
       "      <td>-0.064935</td>\n",
       "      <td>...</td>\n",
       "      <td>34.214286</td>\n",
       "      <td>0.052192</td>\n",
       "      <td>0.115000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>48.801683</td>\n",
       "      <td>52.925926</td>\n",
       "      <td>34.067308</td>\n",
       "      <td>37.481481</td>\n",
       "      <td>0.119291</td>\n",
       "      <td>0.043704</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-27</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.03</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>78.0</td>\n",
       "      <td>45.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>59.666667</td>\n",
       "      <td>0.307263</td>\n",
       "      <td>...</td>\n",
       "      <td>34.928571</td>\n",
       "      <td>0.288344</td>\n",
       "      <td>0.117143</td>\n",
       "      <td>-0.743902</td>\n",
       "      <td>48.836735</td>\n",
       "      <td>56.666667</td>\n",
       "      <td>34.080432</td>\n",
       "      <td>39.222222</td>\n",
       "      <td>0.119184</td>\n",
       "      <td>0.166667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-28</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>46.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>63.333333</td>\n",
       "      <td>0.010526</td>\n",
       "      <td>...</td>\n",
       "      <td>35.357143</td>\n",
       "      <td>0.301010</td>\n",
       "      <td>0.117143</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>48.854916</td>\n",
       "      <td>57.592593</td>\n",
       "      <td>34.094724</td>\n",
       "      <td>40.444444</td>\n",
       "      <td>0.119041</td>\n",
       "      <td>0.178148</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-29</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>65.0</td>\n",
       "      <td>66.333333</td>\n",
       "      <td>-0.140704</td>\n",
       "      <td>...</td>\n",
       "      <td>35.571429</td>\n",
       "      <td>0.096386</td>\n",
       "      <td>0.095714</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>48.864671</td>\n",
       "      <td>57.111111</td>\n",
       "      <td>34.100599</td>\n",
       "      <td>41.074074</td>\n",
       "      <td>0.118898</td>\n",
       "      <td>0.131111</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-30</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>65.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>62.000000</td>\n",
       "      <td>0.048387</td>\n",
       "      <td>...</td>\n",
       "      <td>36.428571</td>\n",
       "      <td>0.070588</td>\n",
       "      <td>0.018571</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>48.883971</td>\n",
       "      <td>53.962963</td>\n",
       "      <td>34.106459</td>\n",
       "      <td>40.370370</td>\n",
       "      <td>0.118756</td>\n",
       "      <td>0.234444</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-03-31</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>58.333333</td>\n",
       "      <td>-0.091429</td>\n",
       "      <td>...</td>\n",
       "      <td>37.428571</td>\n",
       "      <td>0.068702</td>\n",
       "      <td>0.017857</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>48.888889</td>\n",
       "      <td>54.296296</td>\n",
       "      <td>34.113501</td>\n",
       "      <td>39.703704</td>\n",
       "      <td>0.118614</td>\n",
       "      <td>0.164444</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-04-01</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.11</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>51.0</td>\n",
       "      <td>55.666667</td>\n",
       "      <td>-0.119760</td>\n",
       "      <td>...</td>\n",
       "      <td>38.500000</td>\n",
       "      <td>0.090909</td>\n",
       "      <td>0.025714</td>\n",
       "      <td>3.277778</td>\n",
       "      <td>58.819322</td>\n",
       "      <td>54.407407</td>\n",
       "      <td>43.495609</td>\n",
       "      <td>40.296296</td>\n",
       "      <td>0.127880</td>\n",
       "      <td>0.176296</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-04-02</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.16</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>51.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>51.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>39.500000</td>\n",
       "      <td>0.088608</td>\n",
       "      <td>0.037143</td>\n",
       "      <td>3.307692</td>\n",
       "      <td>58.809524</td>\n",
       "      <td>54.000000</td>\n",
       "      <td>43.494987</td>\n",
       "      <td>40.407407</td>\n",
       "      <td>0.127920</td>\n",
       "      <td>0.104815</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2007-04-03</th>\n",
       "      <td>USW00094789</td>\n",
       "      <td>JFK INTERNATIONAL AIRPORT, NY US</td>\n",
       "      <td>0.01</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>53.0</td>\n",
       "      <td>41.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>51.000000</td>\n",
       "      <td>0.039216</td>\n",
       "      <td>...</td>\n",
       "      <td>40.071429</td>\n",
       "      <td>0.023173</td>\n",
       "      <td>0.037857</td>\n",
       "      <td>-0.735849</td>\n",
       "      <td>58.802253</td>\n",
       "      <td>53.666667</td>\n",
       "      <td>43.491865</td>\n",
       "      <td>40.555556</td>\n",
       "      <td>0.127772</td>\n",
       "      <td>0.173704</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>15 rows × 26 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                              name  prcp  snow  snwd  \\\n",
       "DATE                                                                          \n",
       "2007-03-20  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   1.0   \n",
       "2007-03-21  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2007-03-22  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.04   0.0   0.0   \n",
       "2007-03-23  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.14   0.0   0.0   \n",
       "2007-03-24  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.04   0.0   0.0   \n",
       "2007-03-25  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2007-03-26  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2007-03-27  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.03   0.0   0.0   \n",
       "2007-03-28  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2007-03-29  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2007-03-30  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2007-03-31  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.00   0.0   0.0   \n",
       "2007-04-01  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.11   0.0   0.0   \n",
       "2007-04-02  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.16   0.0   0.0   \n",
       "2007-04-03  USW00094789  JFK INTERNATIONAL AIRPORT, NY US  0.01   0.0   0.0   \n",
       "\n",
       "            tmax  tmin  target  rolling_3_tmax  rolling_3_tmax_pct  ...  \\\n",
       "DATE                                                                ...   \n",
       "2007-03-20  52.0  33.0    40.0       43.666667            0.190840  ...   \n",
       "2007-03-21  40.0  26.0    54.0       44.000000           -0.090909  ...   \n",
       "2007-03-22  54.0  40.0    61.0       48.666667            0.109589  ...   \n",
       "2007-03-23  61.0  46.0    53.0       51.666667            0.180645  ...   \n",
       "2007-03-24  53.0  41.0    53.0       56.000000           -0.053571  ...   \n",
       "2007-03-25  53.0  37.0    48.0       55.666667           -0.047904  ...   \n",
       "2007-03-26  48.0  36.0    78.0       51.333333           -0.064935  ...   \n",
       "2007-03-27  78.0  45.0    64.0       59.666667            0.307263  ...   \n",
       "2007-03-28  64.0  46.0    57.0       63.333333            0.010526  ...   \n",
       "2007-03-29  57.0  39.0    65.0       66.333333           -0.140704  ...   \n",
       "2007-03-30  65.0  39.0    53.0       62.000000            0.048387  ...   \n",
       "2007-03-31  53.0  40.0    49.0       58.333333           -0.091429  ...   \n",
       "2007-04-01  49.0  42.0    51.0       55.666667           -0.119760  ...   \n",
       "2007-04-02  51.0  43.0    53.0       51.000000            0.000000  ...   \n",
       "2007-04-03  53.0  41.0    44.0       51.000000            0.039216  ...   \n",
       "\n",
       "            rolling_14_tmin  rolling_14_tmin_pct  rolling_14_prcp  \\\n",
       "DATE                                                                \n",
       "2007-03-20        27.928571             0.181586         0.120000   \n",
       "2007-03-21        28.857143            -0.099010         0.116429   \n",
       "2007-03-22        30.500000             0.311475         0.119286   \n",
       "2007-03-23        32.785714             0.403050         0.129286   \n",
       "2007-03-24        33.857143             0.210970         0.131429   \n",
       "2007-03-25        33.785714             0.095137         0.115000   \n",
       "2007-03-26        34.214286             0.052192         0.115000   \n",
       "2007-03-27        34.928571             0.288344         0.117143   \n",
       "2007-03-28        35.357143             0.301010         0.117143   \n",
       "2007-03-29        35.571429             0.096386         0.095714   \n",
       "2007-03-30        36.428571             0.070588         0.018571   \n",
       "2007-03-31        37.428571             0.068702         0.017857   \n",
       "2007-04-01        38.500000             0.090909         0.025714   \n",
       "2007-04-02        39.500000             0.088608         0.037143   \n",
       "2007-04-03        40.071429             0.023173         0.037857   \n",
       "\n",
       "            rolling_14_prcp_pct  month_avg_tmax  day_avg_tmax  month_avg_tmin  \\\n",
       "DATE                                                                            \n",
       "2007-03-20            -1.000000       48.782082     47.740741       34.041162   \n",
       "2007-03-21            -1.000000       48.771463     46.666667       34.031439   \n",
       "2007-03-22            -0.664671       48.777778     48.777778       34.038647   \n",
       "2007-03-23             0.082873       48.792521     49.777778       34.053076   \n",
       "2007-03-24            -0.695652       48.797590     49.259259       34.061446   \n",
       "2007-03-25            -1.000000       48.802647     51.888889       34.064982   \n",
       "2007-03-26            -1.000000       48.801683     52.925926       34.067308   \n",
       "2007-03-27            -0.743902       48.836735     56.666667       34.080432   \n",
       "2007-03-28            -1.000000       48.854916     57.592593       34.094724   \n",
       "2007-03-29            -1.000000       48.864671     57.111111       34.100599   \n",
       "2007-03-30            -1.000000       48.883971     53.962963       34.106459   \n",
       "2007-03-31            -1.000000       48.888889     54.296296       34.113501   \n",
       "2007-04-01             3.277778       58.819322     54.407407       43.495609   \n",
       "2007-04-02             3.307692       58.809524     54.000000       43.494987   \n",
       "2007-04-03            -0.735849       58.802253     53.666667       43.491865   \n",
       "\n",
       "            day_avg_tmin  month_avg_prcp  day_avg_prcp  \n",
       "DATE                                                    \n",
       "2007-03-20     34.481481        0.119891      0.155185  \n",
       "2007-03-21     34.037037        0.119746      0.144815  \n",
       "2007-03-22     34.407407        0.119650      0.088148  \n",
       "2007-03-23     34.037037        0.119674      0.080000  \n",
       "2007-03-24     34.888889        0.119578      0.112593  \n",
       "2007-03-25     35.074074        0.119434      0.001481  \n",
       "2007-03-26     37.481481        0.119291      0.043704  \n",
       "2007-03-27     39.222222        0.119184      0.166667  \n",
       "2007-03-28     40.444444        0.119041      0.178148  \n",
       "2007-03-29     41.074074        0.118898      0.131111  \n",
       "2007-03-30     40.370370        0.118756      0.234444  \n",
       "2007-03-31     39.703704        0.118614      0.164444  \n",
       "2007-04-01     40.296296        0.127880      0.176296  \n",
       "2007-04-02     40.407407        0.127920      0.104815  \n",
       "2007-04-03     40.555556        0.127772      0.173704  \n",
       "\n",
       "[15 rows x 26 columns]"
      ]
     },
     "execution_count": 226,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.loc[\"2007-03-20\":\"2007-04-03\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 227,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<Axes: xlabel='diff'>"
      ]
     },
     "execution_count": 227,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGwCAYAAAC3qV8qAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABK1klEQVR4nO3deXhTVf4G8PcmTdI93WjT0lLK1gqUZQrUKtQFpBQEUUZZBRTBBVCpInZGEZcRB2YUF2b4Oa4jKOCMoKAyVLaClMVC2alQCi3QBQpNutBmu78/2gYCBVqa5Cbp+3me+7S59yT55prHvpx77jmCKIoiiIiIiFyITOoCiIiIiJqLAYaIiIhcDgMMERERuRwGGCIiInI5DDBERETkchhgiIiIyOUwwBAREZHL8ZC6AHsxm804e/Ys/Pz8IAiC1OUQERFRE4iiiIqKCkREREAmu34/i9sGmLNnzyIqKkrqMoiIiOgWFBYWIjIy8rrH3TbA+Pn5Aag7Af7+/hJXQ0RERE2h0+kQFRVl+Tt+PW4bYBouG/n7+zPAEBERuZibDf9o9iDezMxMDB8+HBERERAEAatXr77mDRvbFi5caGnTvn37a46/8847Vq+zf/9+DBgwAJ6enoiKisKCBQuaWyoRERG5qWYHmKqqKvTs2ROLFy9u9HhRUZHV9tlnn0EQBIwaNcqq3RtvvGHVbubMmZZjOp0OgwcPRnR0NLKzs7Fw4ULMmzcPH3/8cXPLJSIiIjfU7EtIqampSE1Nve5xjUZj9fj777/HPffcgw4dOljt9/Pzu6Ztg2XLlkGv1+Ozzz6DUqlEt27dkJOTg3fffRfTpk1r9Dm1tbWora21PNbpdE39SERERORi7DoPTElJCX788UdMmTLlmmPvvPMOgoOD0bt3byxcuBBGo9FyLCsrC8nJyVAqlZZ9KSkpyM3NxcWLFxt9r/nz50OtVls23oFERETkvuwaYL788kv4+fnhoYcestr/7LPPYvny5di0aROefPJJvP3223jppZcsx4uLixEWFmb1nIbHxcXFjb5Xeno6tFqtZSssLLTxpyEiIiJnYde7kD777DOMHz8enp6eVvvT0tIsv/fo0QNKpRJPPvkk5s+fD5VKdUvvpVKpbvm5RERE5Frs1gOzdetW5Obm4oknnrhp28TERBiNRpw8eRJA3TiakpISqzYNj683boaIiIhaD7sFmE8//RQJCQno2bPnTdvm5ORAJpMhNDQUAJCUlITMzEwYDAZLm4yMDMTGxiIwMNBeJRMREZGLaHaAqaysRE5ODnJycgAA+fn5yMnJQUFBgaWNTqfDt99+22jvS1ZWFhYtWoR9+/bhxIkTWLZsGWbNmoUJEyZYwsm4ceOgVCoxZcoUHDp0CCtWrMD7779vdemJiIiIWq9mj4H57bffcM8991geN4SKSZMm4YsvvgAALF++HKIoYuzYsdc8X6VSYfny5Zg3bx5qa2sRExODWbNmWYUTtVqN9evXY/r06UhISEBISAjmzp173VuoiYiIqHURRFEUpS7CHnQ6HdRqNbRaLZcSICIichFN/ftt19uoiYiIiOyBAUYCNQaT1CUQERG5NAYYB/s+5wy6vfY/fLotX+pSiIiIXBYDjINtP14Gk1nE/J+O4MBprdTlEBERuSQGGAcrqagBABjNIp5fsReX9LycRERE1FwMMA5WoqtbMVsmAHnnqvDOz0ckroiIiMj1MMA4WKmurgdmzpA4AMCXWaewObdUypKIiIhcDgOMA+mNZpRV6QEAD/eJwqSkaADA7P/sx4X6/URERHRzDDAOdL6y7vKRQi4g0FuBl1NvQ6dQX5yrqMWfvjsAN51TkIiIyOYYYByopP7yUaifJwRBgJdSjkWje8FDJmDdoWJ8m31a4gqJiIhcAwOMAzUM4A31V1n2dW+rxqz7ugAAXv/hEArKqiWpjYiIyJUwwDhQaf0t1GF+nlb7n7qrI/q2D0SV3oRZK3NgNJmlKI+IiMhlMMA4UMMlpLAremAAQC4T8O4jveCr8kD2qYtYsiVPivKIiIhcBgOMA5VaLiF5XnMsKsgbr4/oBgBY9Msx7D9d7sjSiIiIXAoDjAOVVNQFmLBGAgwAPPSHthgar6mfpTeHs/QSERFdBwOMA5Va7kJSNXpcEAT8ZWQ8Qv1UOHGuCm//xFl6iYiIGsMA40CXx8A03gMDAIE+Svzt4Z4AgK92nMKmo5yll4iI6GoMMA5SazThYrUBwLWDeK+W3KUNJt/RHkDdLL1l9RPgERERUR0GGAdpGMCr9JBB7aW4afuXU+PQOdQX5ytrkc5ZeomIiKwwwDhIqWUArwqCINy0vadCjkVjekEhF7D+cAm+/Y2z9BIRETVggHGQ0iuWEWiqbhFqpN0XCwCYt+YQTpVV2aU2IiIiV8MA4yDXm8TuZqYld0C/mCBU602YtYKz9BIREQEMMA7TMAdMc3pggIZZenvCT+WBPQXl+OdmztJLRETEAOMgTbmF+noiA73xxsj6WXo3HMO+wnJblkZERORyGGAc5NwVg3hvxchebTGsRzhMZhGzVuSgWm+0ZXlEREQuhQHGQVrSAwM0zNLbHRp/T5w4z1l6iYiodWOAcZCShoUcr7OMQFMEeF+epXfpjgJk/n7OJrURERG5GgYYB6gxmKC9VDcLb2MrUTdH/84hmJQUDQB46T/7oa2f3ZeIiKg1YYBxgIZZeD0VMvh7erT49V5OvQ0xIT4o1tXg9TWHWvx6REREroYBxgFKKi6Pf2nKLLw346WU428P94RMAL7bewbrDha3+DWJiIhcCQOMAzT0wIQ1cw6YG0mIDsSTd3UEAPx51QGc54KPRETUijDAOEDDHUhtbvEW6ut5flBnxGn8UFalx5+44CMREbUiDDAOYLmEZMMeGABQecjx7iOXF3z8bs8Zm74+ERGRs2KAcQDLJSQb98AAQNcIfzw/qAsAYN4Ph3C2/JLN34OIiMjZMMA4QEsnsbuZJ5M7oHe7AFTUGvHSf/bDbOalJCIicm8MMA7QEGBC7dADAwAechn+/nBPeCpk2Hb8PJbuPGWX9yEiInIWDDAOUGpZB8k+PTAA0KGNL14eEgcAePunI8g/X2W39yIiIpIaA4ydVeuNqKipW3ixJcsINMXEpPZI6hCMGoMZL6zMgYmXkoiIyE01O8BkZmZi+PDhiIiIgCAIWL16tdXxyZMnQxAEq23IkCFWbS5cuIDx48fD398fAQEBmDJlCiorK63a7N+/HwMGDICnpyeioqKwYMGC5n86J9AwgNdbKYevquWz8N6ITCZg4cM94KvywJ6CcvxfZp5d34+IiEgqzQ4wVVVV6NmzJxYvXnzdNkOGDEFRUZFl++abb6yOjx8/HocOHUJGRgbWrl2LzMxMTJs2zXJcp9Nh8ODBiI6ORnZ2NhYuXIh58+bh448/bm65krtyAK8tZuG9mchAb8wd3hUA8F7G7zharLP7exIRETlas7sEUlNTkZqaesM2KpUKGo2m0WNHjhzBunXrsHv3bvTp0wcA8OGHH2Lo0KH429/+hoiICCxbtgx6vR6fffYZlEolunXrhpycHLz77rtWQccVlFS0fBXq5no4IRLrDxXjlyOlmLViH76ffieUHrxaSERE7sMuf9U2b96M0NBQxMbG4umnn0ZZWZnlWFZWFgICAizhBQAGDRoEmUyGnTt3WtokJydDqVRa2qSkpCA3NxcXL15s9D1ra2uh0+msNmdQaudbqBsjCALefigegd4KHCnS4YMNxxz23kRERI5g8wAzZMgQ/Pvf/8aGDRvw17/+FVu2bEFqaipMJhMAoLi4GKGhoVbP8fDwQFBQEIqLiy1twsLCrNo0PG5oc7X58+dDrVZbtqioKFt/tFtSKkEPTN37eeIvD8YDAP6x+Tj2FjQe/IiIiFyRzQPMmDFjMGLECMTHx2PkyJFYu3Ytdu/ejc2bN9v6raykp6dDq9VatsLCQru+X1PZexK7GxkaH44HekXALAIvrNyHS3qTw2sgIiKyB7sPjOjQoQNCQkJw/PhxAIBGo0FpaalVG6PRiAsXLljGzWg0GpSUlFi1aXh8vbE1KpUK/v7+VpszsPckdjfzxojuCPNX4cT5Kvx13VFJaiAiIrI1uweY06dPo6ysDOHh4QCApKQklJeXIzs729Jm48aNMJvNSExMtLTJzMyEwWCwtMnIyEBsbCwCAwPtXbJNXV4HyfE9MACg9lbgr6N6AAC+2H4S24+fl6QOIiIiW2p2gKmsrEROTg5ycnIAAPn5+cjJyUFBQQEqKysxe/Zs7NixAydPnsSGDRvwwAMPoFOnTkhJSQEA3HbbbRgyZAimTp2KXbt24ddff8WMGTMwZswYREREAADGjRsHpVKJKVOm4NChQ1ixYgXef/99pKWl2e6TO4iUl5Aa3B0birH92gEAZv9nP3Q1hps8g4iIyLk1O8D89ttv6N27N3r37g0ASEtLQ+/evTF37lzI5XLs378fI0aMQJcuXTBlyhQkJCRg69atUKkuX0JZtmwZ4uLiMHDgQAwdOhT9+/e3muNFrVZj/fr1yM/PR0JCAl544QXMnTvX5W6hrqw1oqp+3ImjB/Fe7c/DbkNUkBfOlF/Cm2sOS1oLERFRSwmiKLrlfPM6nQ5qtRparVay8TAnzlXi3r9vga/KAwdfT5Gkhivtyr+A0R9nQRSBTyb2waCuYTd/EhERkQM19e83Zzezo5L68S9SDeC9Wr+YIDzRPwYA8Nd1R+Gm2ZWIiFoBBhg7Kq2oH//iJ934l6vNHNgZngoZjpVWYg/nhiEiIhfFAGNHlwfwOkcPDAD4eyowLL5usPTyXc4xVw4REVFzMcDYkdS3UF/P2H51sxSv3V+ECt6RRERELogBxo4aFnJsI/EdSFdLiA5Ep1BfXDKY8MO+s1KXQ0RE1GwMMHbkDHPANEYQBIzpW9cLw8tIRETkihhg7EiKlaib6sHebaGQCzhwRouDZ7RSl0NERNQsDDB2Ioqi5TZqZxrE2yDYV4XB3erWlVqxm70wRETkWhhg7KSi1ohLhoZZeJ2vBwaA5TLS6pwzXKmaiIhcCgOMnTTcgeTn6QEvpVziahp3Z8cQRAZ6oaLGiJ8OFEldDhERUZMxwNiJM49/aSCTCRjdp64XhpeRiIjIlTDA2ElJhfNNYteYh/tEQSYAu05ewPHSSqnLISIiahIGGDuxDOB10vEvDTRqT9wbFwoAWPkbe2GIiMg1MMDYScMcMKFOfAmpwei+7QAA/80+Db3RLHE1REREN8cAYyelTnwL9dXuiW2DUD8Vyqr0+OVIidTlEBER3RQDjJ00rETtrLdQX8lDLsPDfSIBAN/sKpC4GiIioptjgLETZ57ErjGj+9RdRtp2/DwKL1RLXA0REdGNMcDYQd0svM5/G/WV2gV7485OwRBF4FsO5iUiIifHAGMHuktG1NYPhnW2lahvZEz9YN6Vv52GySxKXA0REdH1McDYQcMcMAHeCngqnHMW3sYM7haGQG8FinU12PJ7qdTlEBERXRcDjB003IEU6kK9LwCg8pDjoT80DOblZSQiInJeDDB24GrjX67UsMDjxqOlluUQiIiInA0DjB2UuNAt1FfrHOaHhOhAmMwivs0+LXU5REREjWKAsQNXmsSuMaPre2FW/lYIMwfzEhGRE2KAsQNXvoQEAPf3CIevygOnyqqx40SZ1OUQERFdgwHGDi4HGNfsgfFWemBErwgAwPLdHMxLRETOhwHGDkor6i4htXHBMTANxtbPCbPuYDEuVuklroaIiMgaA4yNiaLo8mNgACA+Uo1uEf7Qm8xYtfeM1OUQERFZYYCxsfJqA/Qm15uFtzENt1Qv310AUeRgXiIich4MMDbWcAt1kI8SKg/XmYW3MSN6tYWnQobfSyqxt7Bc6nKIiIgsGGBsrMRFZ+FtjNpLgaHx4QCA5bsKJK6GiIjoMgYYG2uYvTbURW+hvtrYfnWDedfsK0JFjUHiaoiIiOowwNhYwx1IYW7QAwMAfaID0bGNDy4ZTFizr0jqcoiIiAAwwNicq09idzVBEDCm/pbqFbt5GYmIiJwDA4yNufokdo156A9toZAL2Hdai8NndVKXQ0RExABja5ZBvG7SAwMAwb4qDO6qAcBeGCIicg4MMDZW6maXkBo0LPC4au8Z1BhMEldDREStHQOMDZnNIs5Vus9t1Ffq3ykEbQO8oKsx4ueDHMxLRETSYoCxoYvVehhMdTPWuvosvFeTyQRLL8w3u7jAIxERSavZASYzMxPDhw9HREQEBEHA6tWrLccMBgPmzJmD+Ph4+Pj4ICIiAhMnTsTZs2etXqN9+/YQBMFqe+edd6za7N+/HwMGDICnpyeioqKwYMGCW/uEDtQw/iXEVwmF3P2y4cN9IiETgF35F3DiXKXU5RARUSvW7L+yVVVV6NmzJxYvXnzNserqauzZswevvvoq9uzZg++++w65ubkYMWLENW3feOMNFBUVWbaZM2dajul0OgwePBjR0dHIzs7GwoULMW/ePHz88cfNLdehGpYRCHXhVahvJFzthbtjQwEAK3azF4aIiKTj0dwnpKamIjU1tdFjarUaGRkZVvs++ugj9OvXDwUFBWjXrp1lv5+fHzQaTaOvs2zZMuj1enz22WdQKpXo1q0bcnJy8O6772LatGmNPqe2tha1tbWWxzqd42/3LXXDW6ivNqZvFDYeLcWK3woxc2Bn+Kqa/RUiIiJqMbtf59BqtRAEAQEBAVb733nnHQQHB6N3795YuHAhjEaj5VhWVhaSk5OhVCot+1JSUpCbm4uLFy82+j7z58+HWq22bFFRUXb5PDdyeR0k9+yBAYCBt4WhQ4gPyqsN+CrrlNTlEBFRK2XXAFNTU4M5c+Zg7Nix8Pf3t+x/9tlnsXz5cmzatAlPPvkk3n77bbz00kuW48XFxQgLC7N6rYbHxcXFjb5Xeno6tFqtZSssdPwljtIK9++BkcsEPHNPJwDAJ1tPoFpvvMkziIiIbM9u/f8GgwGPPPIIRFHEP//5T6tjaWlplt979OgBpVKJJ598EvPnz4dKdWt//FUq1S0/11bccRK7xjzQKwIfbDiGggvV+HpnAZ4Y0EHqkoiIqJWxSw9MQ3g5deoUMjIyrHpfGpOYmAij0YiTJ08CADQaDUpKSqzaNDy+3rgZZ+Cuk9hdTSGX4Zm7OwIA/i/zBCe2IyIih7N5gGkIL8eOHcMvv/yC4ODgmz4nJycHMpkMoaF1d7gkJSUhMzMTBoPB0iYjIwOxsbEIDAy0dck209AD486XkBo89IdItA3wwrmKWt6RREREDtfsAFNZWYmcnBzk5OQAAPLz85GTk4OCggIYDAb88Y9/xG+//YZly5bBZDKhuLgYxcXF0Ov1AOoG6C5atAj79u3DiRMnsGzZMsyaNQsTJkywhJNx48ZBqVRiypQpOHToEFasWIH333/f6tKTszFZzcLr3j0wAKD0kOGp+l6YJVvyUGtkLwwRETmOIIqi2JwnbN68Gffcc881+ydNmoR58+YhJiam0edt2rQJd999N/bs2YNnnnkGR48eRW1tLWJiYvDoo48iLS3NagzL/v37MX36dOzevRshISGYOXMm5syZ0+Q6dTod1Go1tFrtTS9h2cK5ilr0/csvEATg2Fup8HDDieyuVmMw4a6Fm1Ciq8VfHuyO8YnRUpdEREQurql/v5sdYFyFowPMwTNa3P/hNrTxU2H3nwfZ/f2cxWfb8vHG2sOIDPTCphfvdssZiImIyHGa+vebf21spDXcQt2Ysf3aIcRXhdMXL2HV3jNSl0NERK0EA4yNWAbwtoLxL1fyUsoxLbnusuHiTcdhNJklroiIiFoDBhgbKam/hdrd54BpzPjEaAT5KHGqrBpr9p+9+ROIiIhaiAHGRi4vI9C6LiEBgI/KA1P61/XCfLjxOExmtxxWRUREToQBxkbOVbSOSeyuZ2JSNNReCpw4V4WfDhRJXQ4REbk5BhgbaU2T2DXGz1OBx++s64X5aONxmNkLQ0REdsQAYyMlrWQZgRuZfGd7+Kk8kFtSgfWHG190k4iIyBYYYGzAaDLjfMMsvK20BwYA1F4KTLqjPYC6sTBuOsUQERE5AQYYGyir0sMsAjIBCPZpvQEGAKb0j4G3Uo5DZ3XYeLRU6nKIiMhNMcDYQGn9+Jc2firIZYLE1Ugr0EeJR5PqlhT4YMMx9sIQEZFdMMDYAMe/WJs6oAM8FTLsO61F5rHzUpdDRERuiAHGBkrqb6FuDatQN0WIr8qysCN7YYiIyB4YYGygtd9C3ZgnkztA6SFD9qmLyMork7ocIiJyMwwwNlDKS0jXCPX3xJi+UQCADzYek7gaIiJyNwwwNmBZB6kVLiNwI0/d1REKuYAdJy5gV/4FqcshIiI3wgBjA6UVDZeQ2ANzpYgAL/wxoa4X5kP2whARkQ0xwNiAZSFHjoG5xjN3d4RcJmDrsfPYU3BR6nKIiMhNMMC0kMFkRlkVe2CuJyrIGw/1bgsA+HADe2GIiMg2GGBa6HxlLUQR8JAJCPJWSl2OU5p+TyfIBGBT7jkcOK2VuhwiInIDDDAtVHLFLLyyVj4L7/W0D/HBiJ4RADgWhoiIbIMBpoUsdyDx8tENzbi3EwQBWH+4BEeKdFKXQ0RELo4BpoUsdyDxFuob6hTqh6Hx4QCAjzYel7gaIiJydQwwLcRJ7Jpu5r2dAAA/HSzCsZIKiashIiJXxgDTQpcXcmQPzM3EafyR0i0Mogh8tIm9MEREdOsYYFro8hww7IFpipn3dgYArNl3FoUXqiWuhoiIXBUDTAtxGYHm6d5WjQGdQ2AWgc9+zZe6HCIiclEMMC10jssINNsTAzoAAFbsLoS22iBxNURE5IoYYFpAbzSjrEoPgAGmOZI7hyA2zA/VehO+3lUgdTlEROSCGGBa4FxlXe+LQi4g0FshcTWuQxAETE2u64X5/Nd86I1miSsiIiJXwwDTApfHv3hCEDgLb3OM6BmBUD8VSitq8cO+s1KXQ0RELoYBpgVKLbPwcgBvcyk9ZJh8Z3sAwCdbT0AURWkLIiIil8IA0wINt1CH+XH8y60Y3y8a3ko5jhZXYOux81KXQ0RELoQBpgVKKziJXUuovRV4pE8UAOBfW09IXA0REbkSBpgW4CR2LTelfwxkArD12Hku8khERE3GANMCJVwHqcWigryR2r1ukUf2whARUVMxwLRAaUMPDGfhbZEnBsQAqFteoFhbI3E1RETkChhgWqCkgj0wttC7XSD6tQ+CwSTii+0npS6HiIhcAAPMLaoxmFBePw0+B/G2XEMvzNc7T6Gy1ihxNURE5OyaHWAyMzMxfPhwREREQBAErF692uq4KIqYO3cuwsPD4eXlhUGDBuHYsWNWbS5cuIDx48fD398fAQEBmDJlCiorK63a7N+/HwMGDICnpyeioqKwYMGC5n86O2pYA0npIYPai7PwttSg28IQE+IDXY0RK3cXSl0OERE5uWYHmKqqKvTs2ROLFy9u9PiCBQvwwQcfYMmSJdi5cyd8fHyQkpKCmprLYxvGjx+PQ4cOISMjA2vXrkVmZiamTZtmOa7T6TB48GBER0cjOzsbCxcuxLx58/Dxxx/fwke0jytvoeYsvC0nkwmY0r+uF+bTbfkwmri8ABER3YDYAgDEVatWWR6bzWZRo9GICxcutOwrLy8XVSqV+M0334iiKIqHDx8WAYi7d++2tPn5559FQRDEM2fOiKIoiv/4xz/EwMBAsba21tJmzpw5YmxsbJNr02q1IgBRq9Xe6se7oR/3nxWj56wVR/3jV7u8fmtUXWsUe7+xXoyes1Zcs++M1OUQEZEEmvr326ZjYPLz81FcXIxBgwZZ9qnVaiQmJiIrKwsAkJWVhYCAAPTp08fSZtCgQZDJZNi5c6elTXJyMpRKpaVNSkoKcnNzcfHixUbfu7a2FjqdzmqzpxIuI2BzXko5Hr09GgDwr0wuL0BERNdn0wBTXFwMAAgLC7PaHxYWZjlWXFyM0NBQq+MeHh4ICgqyatPYa1z5HlebP38+1Gq1ZYuKimr5B7oByyR2XEbAph5NiobSQ4Z9p7XYfbLxsEpEROQ2dyGlp6dDq9VatsJC+w4ELeUt1HYR4qvCqD9EAgA+zuTEdkRE1DibBhiNRgMAKCkpsdpfUlJiOabRaFBaWmp13Gg04sKFC1ZtGnuNK9/jaiqVCv7+/labPTVMYsdbqG2vYTDvL0dKkHeu8iatiYioNbJpgImJiYFGo8GGDRss+3Q6HXbu3ImkpCQAQFJSEsrLy5GdnW1ps3HjRpjNZiQmJlraZGZmwmAwWNpkZGQgNjYWgYGBtiz5lnEZAfvpFOqLQbfVXWb8dFu+xNUQEZEzanaAqaysRE5ODnJycgDUDdzNyclBQUEBBEHA888/j7feegs//PADDhw4gIkTJyIiIgIjR44EANx2220YMmQIpk6dil27duHXX3/FjBkzMGbMGERERAAAxo0bB6VSiSlTpuDQoUNYsWIF3n//faSlpdnsg7eUZRAvlxGwi6kDOgAA/pt9GmWVtRJXQ0REzqbZAea3335D79690bt3bwBAWloaevfujblz5wIAXnrpJcycORPTpk1D3759UVlZiXXr1sHT83JPxbJlyxAXF4eBAwdi6NCh6N+/v9UcL2q1GuvXr0d+fj4SEhLwwgsvYO7cuVZzxUjpkt4EXU3dbLFcido++sUEoUekGrVGM77acUrqcoiIyMkIopveq6rT6aBWq6HVam0+HuZUWRXuWrgZngoZjrwxhBPZ2cmafWcx85u9CPJRYvvL98JTIZe6JCIisrOm/v12m7uQHKm0omEAryfDix2ldtegbYAXLlTp8d89p6Uuh4iInAgDzC2wDODlHDB25SGX4fGG5QW25sNsdsvOQiIiugUMMLfAMokdb6G2u9F9o+Dn6YET56uw4WjpzZ9AREStAgPMLSi13IHEHhh781V5YFxiOwDAv7ZyYjsiIqrDAHMLLs8Bwx4YR3jsjhh4yATsyr+AnMJyqcshIiInwABzC64cxEv2p1F7YkSvujmC2AtDREQAA8wt4UrUjvdE/7qJ7X4+UITCC9USV0NERFJjgLkFl9dBYg+Mo3SN8MeAziEwi8Bnv3J5ASKi1o4Bppmqao2oqK2fhZfLCDjUE/XLC6zYXQhtteEmrYmIyJ0xwDRTw/gXb6UcvioPiatpXZI7hyA2zA/VehO+3lUgdTlERCQhBphmunIVas7C61iCIGBqcl0vzOe/5qPGYJK4IiIikgoDTDM19MDw8pE0RvSMQLjaE6UVtVjKRR6JiFotBphmKr2iB4YcT+khw7MDOwMA/rE5D5X145GIiKh1YYBpJk5iJ72HEyIRE+KDC1V6fMJ5YYiIWiUGmGayrIPEZQQk4yGX4YXBXQAAn2zNx4UqvcQVERGRozHANFOf9oEY1iMct4X7S11Kqza0ezi6RfijstaIf2w6LnU5RETkYAwwzTQxqT0Wj/sD+ncOkbqUVk0mEzA7JRYA8O8dp3C2/JLEFRERkSMxwJDLuqtLG/SLCYLeaMYHG45JXQ4RETkQAwy5LEEQMGdIXS/Mt9mnceJcpcQVERGRozDAkEtLiA7CwLhQmMwi/p7xu9TlEBGRgzDAkMt7MSUWggD8uL8IB89opS6HiIgcgAGGXN5t4f54oGcEAGDB/3IlroaIiByBAYbcwqz7usBDJiDz93PYcaJM6nKIiMjOGGDILUQH+2BMvygAwIJ1RyGKosQVERGRPTHAkNt49t7O8FTIsKegHBuOlEpdDhER2REDDLmNUH9PTL4jBgDwt/W5MJvZC0NE5K4YYMitPH1XR/h5euBocQV+2HdW6nKIiMhOGGDIrai9FXjqro4AgHczfofeaJa4IiIisgcGGHI7j93ZHiG+KhRcqMaK3QVSl0NERHbAAENux1vpgWcHdgIAfLDxOKr1RokrIiIiW2OAIbc0pm87RAV54VxFLb7YflLqcoiIyMYYYMgtKT1kmDWoCwBgyeY8aKsNEldERES2xABDbuuBXm3RJcwXuhoj/i8zT+pyiIjIhhhgyG3JZQJeHBwLAPj815MoraiRuCIiIrIVBhhya/d1DUPvdgG4ZDDho43HpS6HiIhshAGG3JogCJidUtcL882uAhReqJa4IiIisgUGGHJ7d3QMwYDOITCYRLyX8bvU5RARkQ0wwFCr8FJKHABgVc4ZHC3WSVwNERG1FAMMtQrxkWoMjddAFIG//Y+9MERErs7mAaZ9+/YQBOGabfr06QCAu++++5pjTz31lNVrFBQUYNiwYfD29kZoaChmz54No5GzqVLLpN0XC5kA/HKkBNmnLkpdDhERtYCHrV9w9+7dMJlMlscHDx7Efffdh4cfftiyb+rUqXjjjTcsj729vS2/m0wmDBs2DBqNBtu3b0dRUREmTpwIhUKBt99+29blUivSKdQXf0yIxMrfTmPh/47im6m3QxAEqcsiIqJbYPMemDZt2kCj0Vi2tWvXomPHjrjrrrssbby9va3a+Pv7W46tX78ehw8fxtKlS9GrVy+kpqbizTffxOLFi6HX621dLrUyzw3qAqVchh0nLuC7PWekLoeIiG6RXcfA6PV6LF26FI8//rjVv3SXLVuGkJAQdO/eHenp6aiuvnxra1ZWFuLj4xEWFmbZl5KSAp1Oh0OHDl33vWpra6HT6aw2oqu1DfDCc4M6AwDmfn8QJ89XSVwRERHdCrsGmNWrV6O8vByTJ0+27Bs3bhyWLl2KTZs2IT09HV999RUmTJhgOV5cXGwVXgBYHhcXF1/3vebPnw+1Wm3ZoqKibPthyG08dVdHJMYEoUpvwrPL90JvNEtdEhERNZNdA8ynn36K1NRUREREWPZNmzYNKSkpiI+Px/jx4/Hvf/8bq1atQl5ey9aqSU9Ph1artWyFhYUtLZ/clFwm4L3RvaD2UmD/aS3e+4V3JRERuRq7BZhTp07hl19+wRNPPHHDdomJiQCA48frpnnXaDQoKSmxatPwWKPRXPd1VCoV/P39rTai64kI8MI7D8UDAJZsycP24+clroiIiJrDbgHm888/R2hoKIYNG3bDdjk5OQCA8PBwAEBSUhIOHDiA0tJSS5uMjAz4+/uja9eu9iqXWqHU+HCM7RcFUQRmrczBxSoOEicichV2CTBmsxmff/45Jk2aBA+Py3dq5+Xl4c0330R2djZOnjyJH374ARMnTkRycjJ69OgBABg8eDC6du2KRx99FPv27cP//vc/vPLKK5g+fTpUKpU9yqVW7NX7u6JDGx+U6Gox57/7IYqi1CUREVET2CXA/PLLLygoKMDjjz9utV+pVOKXX37B4MGDERcXhxdeeAGjRo3CmjVrLG3kcjnWrl0LuVyOpKQkTJgwARMnTrSaN4bIVryVHvhgTG8o5TKsP1yCZTsLpC6JiIiaQBDd9J+cOp0OarUaWq2W42Hopj7ZegJv/XgEKg8Z1s7sj85hflKXRETUKjX17zfXQiIC8PidMUju0ga1RjNmfrMXNQbTzZ9ERESSYYAhAiCTCfjbwz0Q7KPE0eIK/HXdUalLIiKiG2CAIaoX6ueJvz3cEwDw+a8nselo6U2eQUREUmGAIbrCPXGhmHxHewDA7P/sw7mKWmkLIiKiRjHAEF3l5dQ4xGn8cL5Sjxe+3Qez2S3HuRMRuTQGGKKreCrk+HBsb6g8ZMj8/Rw++zVf6pKIiOgqDDBEjegc5odX76+b+fmv647i4BmtxBUREdGVGGCIrmN8Yjvc1zUMBpOI55bvRbXeKHVJRERUjwGG6DoEQcBfR/VAmL8Keeeq8ObaI1KXRERE9RhgiG4gyEeJdx/pBUEAvtlVgHUHi6QuiYiIwABDdFN3dgrBk8kdAQBz/nsAZ8svSVwRERExwBA1Qdp9XdAjUg3tJQNmrciBibdWExFJigGGqAmUHjJ8MKY3vJVy7My/gCVb8qQuiYioVWOAIWqi9iE+eOOB7gCAdzN+563VREQSYoAhaoZRf2iLofEamMwiXll9kLP0EhFJhAGGqBkEQcBrw7vBV+WBnMJyrPitUOqSiIhaJQYYomYK8/fErPu6AKibpfdClV7iioiIWh8GGKJbMCkpGnEaP5RXG7Bg3VGpyyEianUYYIhugYdchrdG1g3oXb67ENmnLkpcERFR68IAQ3SL+rQPwh8TIgEAr64+CKPJLHFFREStBwMMUQukp8ZB7aXA4SIdlu44JXU5REStBgMMUQsE+6owOyUWAPD39b+jVFcjcUVERK0DAwxRC43t1w49ItWoqDXi7Z+4YjURkSMwwBC1kFwm4K2R3SEIwOqcs8jKK5O6JCIit8cAQ2QDPSIDMD6xHQBg7vcHoTdyQC8RkT0xwBDZyOzBcQj2UeJYaSU++zVf6nKIiNwaAwyRjai9FXg5NQ4A8P4vx3C2/JLEFRERuS8GGCIbGvWHSPRtH4hLBhPeXHtY6nKIiNwWAwyRDclkAt4c2R1ymYCfDxZjc26p1CUREbklBhgiG4vT+GPyHe0BAK/9cAg1BpO0BRERuSEGGCI7eH5QZ4T5q3CqrBr/t+WE1OUQEbkdBhgiO/DzVOCVYV0BAIs3H8epsiqJKyIici8MMER2cn+PcPTvFAK90Yx5PxyCKIpSl0RE5DYYYIjsRBAEvP5ANyjkAjblnsP6wyVSl0RE5DYYYIjsqGMbX0xL7gAAeGPNYVTrjRJXRETkHhhgiOxsxj2d0TbAC2fKL+HDjcelLoeIyC0wwBDZmZdSjteG1w3o/WTrCRwvrZS4IiIi18cAQ+QA93UNw8C4UBhMIuZ+f5ADeomIWsjmAWbevHkQBMFqi4uLsxyvqanB9OnTERwcDF9fX4waNQolJdaDGwsKCjBs2DB4e3sjNDQUs2fPhtHIsQPkugRBwLwR3aDykGF7Xhl+2HdW6pKIiFyaXXpgunXrhqKiIsu2bds2y7FZs2ZhzZo1+Pbbb7FlyxacPXsWDz30kOW4yWTCsGHDoNfrsX37dnz55Zf44osvMHfuXHuUSuQwUUHemH5PJwDAX348Al2NQeKKiIhcl10CjIeHBzQajWULCQkBAGi1Wnz66ad49913ce+99yIhIQGff/45tm/fjh07dgAA1q9fj8OHD2Pp0qXo1asXUlNT8eabb2Lx4sXQ6/X2KJfIYaYld0D7YG+UVtQi/b8HeCmJiOgW2SXAHDt2DBEREejQoQPGjx+PgoICAEB2djYMBgMGDRpkaRsXF4d27dohKysLAJCVlYX4+HiEhYVZ2qSkpECn0+HQoUPXfc/a2lrodDqrjcjZeCrk+PsjvaCQC/jxQBE+3ZYvdUlERC7J5gEmMTERX3zxBdatW4d//vOfyM/Px4ABA1BRUYHi4mIolUoEBARYPScsLAzFxcUAgOLiYqvw0nC84dj1zJ8/H2q12rJFRUXZ9oMR2UhCdCBevb/urqT5Px/FjhNlEldEROR6bB5gUlNT8fDDD6NHjx5ISUnBTz/9hPLycqxcudLWb2UlPT0dWq3WshUWFtr1/Yha4tHbo/Fg77YwmUXM+HoPirU1UpdERORS7H4bdUBAALp06YLjx49Do9FAr9ejvLzcqk1JSQk0Gg0AQKPRXHNXUsPjhjaNUalU8Pf3t9qInJUgCHj7wXjEafxwvlKPZ5ZlQ280S10WEZHLsHuAqaysRF5eHsLDw5GQkACFQoENGzZYjufm5qKgoABJSUkAgKSkJBw4cAClpaWWNhkZGfD390fXrl3tXS6Rw3gp5VgyIQF+nh7YU1COv/x4WOqSiIhchs0DzIsvvogtW7bg5MmT2L59Ox588EHI5XKMHTsWarUaU6ZMQVpaGjZt2oTs7Gw89thjSEpKwu233w4AGDx4MLp27YpHH30U+/btw//+9z+88sormD59OlQqla3LJZJU+xAfLBrdCwDwZdYprNp7WtqCiIhchM0DzOnTpzF27FjExsbikUceQXBwMHbs2IE2bdoAAN577z3cf//9GDVqFJKTk6HRaPDdd99Zni+Xy7F27VrI5XIkJSVhwoQJmDhxIt544w1bl0rkFAbeFoZn762bHyb9uwM4UsQ76IiIbkYQ3XQiCp1OB7VaDa1Wy/Ew5PRMZhGPfbEbmb+fQ3SwN36Y0R9qL4XUZREROVxT/35zLSQiJyCXCXh/dC+0DfDCqbJqvLAyB2azW/7bgojIJhhgiJxEoI8SSyYkQOkhwy9HSvGPzcelLomIyGkxwBA5kfhINd56oDsA4O8ZvyPz93MSV0RE5JwYYIiczCN9ozC2XxREEXh2+V4UXqiWuiQiIqfDAEPkhF4b3g09ItUorzbgmWV7UGMwSV0SEZFTYYAhckKeCjn+Mf4PCPRW4MAZLV77/voLmRIRtUYMMEROKjLQGx+M7Q1BAFb8VojluwqkLomIyGkwwBA5sQGd2+DFwbEAgLnfH8K+wnJpCyIichIMMERO7um7OmLQbWHQm8x4ZtkeXKjSS10SEZHkGGCInJxMJuDvj/RE+2BvnCm/hOeW74WJk9wRUSvHAEPkAtReCix5NAGeChm2HjuP9zJ+l7okIiJJMcAQuYg4jT/eeagHAOCjTcex4UiJxBUREUmHAYbIhYzs3RaTkqIBAC9+uw8luhqJKyIikgYDDJGL+dOw29A13B8Xqw14fnkOx8MQUavEAEPkYlQecnw4rje8FHJknSjDki15UpdERORwDDBELqhjG1+8/kA3AMC7Gb9jT8FFiSsiInIsBhgiF/VwQiTu7xEOk1nEc8v3QldjkLokIiKHYYAhclGCIOAvD8YjMtALhRcu4c+rDkIUOR6GiFoHBhgiF6b2UuD9Mb0hlwlYs+8s/pN9WuqSiIgcggGGyMUlRAci7b4uAIDXfjiEvHOVEldERGR/DDBEbuCpuzoiqUMwqvUmPPvNXtQaTVKXRERkVwwwRG5ALhPw3uheCPRW4NBZHRauy5W6JCIiu2KAIXITGrUnFv6xJwDgk2352JRbKnFFRET2wwBD5EYGdQ27vNTAyn0oreBSA0TknhhgiNxM+tDbEKfxQ1mVHi+s3AczlxogIjfEAEPkZjwVcnw4tjc8FTJsPXYen2w7IXVJREQ2xwBD5IY6h/nhteF1Sw0sWJeLfYXl0hZERGRjDDBEbmpM3ygMjdfAaBbx7PK9qKw1Sl0SEZHNMMAQuSlBEDD/wR5oG+CFU2XVmLv6oNQlERHZDAMMkRtTeyuwaEwvyATgu71nsGovlxogIvfAAEPk5vq2D8JzA+uWGnhl1UGcKquSuCIiopZjgCFqBWbc2wn9YoJQVb/UgN5olrokIqIWYYAhagXkMgGLRveC2kuBfae1+HsGlxogItfGAEPUSkQEeOGvo3oAAP5vywn8sO+sxBUREd06BhiiVmRIdw3GJ7YDADz7zV78edUB1Bi4cjURuR4GGKJWZt6Ibnjqro4AgGU7CzDio23ILa6QuCoiouZhgCFqZRRyGV5OjcNXU/qhjZ8Kv5dUYsRH2/DVjlMQRa6bRESugQGGqJUa0LkNfn5uAO6JbYNaoxmvrj6IJ7/KxsUqvdSlERHdFAMMUSsW4qvCZ5P74tX7u0IhF7D+cAlS39+KHSfKpC6NiOiGbB5g5s+fj759+8LPzw+hoaEYOXIkcnOtb9m8++67IQiC1fbUU09ZtSkoKMCwYcPg7e2N0NBQzJ49G0Yj13IhsjVBEDClfwxWPXMnOoT4oFhXg7H/2oF31+fCaOJ8MUTknGweYLZs2YLp06djx44dyMjIgMFgwODBg1FVZT3759SpU1FUVGTZFixYYDlmMpkwbNgw6PV6bN++HV9++SW++OILzJ0719blElG97m3VWDOzPx7pEwlRBD7YeByjP96B0xerpS6NiOgagmjnUXvnzp1DaGgotmzZguTkZAB1PTC9evXCokWLGn3Ozz//jPvvvx9nz55FWFgYAGDJkiWYM2cOzp07B6VSedP31el0UKvV0Gq18Pf3t9nnIWoN1uw7iz99dwAVtUb4eXrgnYd6YFiPcKnLIqJWoKl/v+0+Bkar1QIAgoKCrPYvW7YMISEh6N69O9LT01FdfflfeVlZWYiPj7eEFwBISUmBTqfDoUOHGn2f2tpa6HQ6q42Ibs3wnhH46bkB6N0uABU1Rkz/eg9e/u9+VOt5GZeInINdA4zZbMbzzz+PO++8E927d7fsHzduHJYuXYpNmzYhPT0dX331FSZMmGA5XlxcbBVeAFgeFxcXN/pe8+fPh1qttmxRUVF2+ERErUdUkDdWPpmEGfd0giAAy3cXYviH23D4LP9xQETS87Dni0+fPh0HDx7Etm3brPZPmzbN8nt8fDzCw8MxcOBA5OXloWPHjrf0Xunp6UhLS7M81ul0DDFELaSQy/BiSizu6BSMWStykHeuCiMX/4o/DY3DpDvaQxAEqUskolbKbj0wM2bMwNq1a7Fp0yZERkbesG1iYiIA4Pjx4wAAjUaDkpISqzYNjzUaTaOvoVKp4O/vb7URkW3c0TEEPz+XjEG3hUJvMmPemsN4dnkOV7UmIsnYPMCIoogZM2Zg1apV2LhxI2JiYm76nJycHABAeHjdIMGkpCQcOHAApaWlljYZGRnw9/dH165dbV0yETVBkI8S/5rYB6+P6AaFXMCafWfx5Fe/4ZKeaykRkePZ/C6kZ555Bl9//TW+//57xMbGWvar1Wp4eXkhLy8PX3/9NYYOHYrg4GDs378fs2bNQmRkJLZs2QKg7jbqXr16ISIiAgsWLEBxcTEeffRRPPHEE3j77bebVAfvQiKyn825pXhqaTZqDGb0ax+ETyb3gb+nQuqyiMgNNPXvt80DzPWuiX/++eeYPHkyCgsLMWHCBBw8eBBVVVWIiorCgw8+iFdeecWq0FOnTuHpp5/G5s2b4ePjg0mTJuGdd96Bh0fThu0wwBDZ1+6TF/D457tRUWtE97b++PKxfgj2VUldFhG5OMkCjLNggCGyv4NntJj02S6UVenRsY0PvpqSiIgAL6nLIiIX5jTzwBCR++reVo2VTyUhQu2JvHNVeHhJFvLPV938iURELcQAQ0Qt0rGNL759+g50CPHBmfJLeHjJds4VQ0R2xwBDRC3WNsALK59KQtdwf5yv1GPMx1nIPnVB6rKIyI0xwBCRTYT4qvDNtNvRJzoQuhojJnyyC5m/n5O6LCJyUwwwRGQzai8FvpqSiLu6tMElgwlTvtyNnw4USV0WEbkhBhgisikvpRz/mtgHw+LDYTCJmPH1HqzcXSh1WUTkZhhgiMjmlB4yfDC2N8b0jYJZBF767358svWE1GURkRthgCEiu5DLBMx/KB7TkjsAAN768QjeXZ8LN516iogcjAGGiOxGEASkp8ZhdkrdsiIfbDyO19cchtnMEENELcMAQ0R2JQgCpt/TCW8+0A0A8MX2k3jx230wmriSNRHdOgYYInKIR5PaY9HoXpDLBHy39wzG/msHDp7RSl0WEbkoBhgicpiRvdvi/yYkwFMhw+6TFzH8o2148dt9KNHVSF0aEbkYBhgicqhBXcOw8YW7MbJXBEQR+E/2ady9cDPe/+UYLulNUpdHRC6Cq1ETkWT2FlzEWz8eQfapiwCAcLUnXhoSiwd6toVMJkhcHRFJoal/vxlgiEhSoihi7f4ivPPzUZwpvwQA6Bmpxqv3d0Wf9kESV0dEjsYAwwBD5FJqDCZ8ui0f/9h0HFX1l5KG9QjHy0PiEBXkLXF1ROQoDDAMMEQu6VxFLd7NyMWK3YUwi3Wz+j5+Zwym39MRfp4KqcsjIjtjgGGAIXJpR4p0eOvHw/j1eBkAIMRXibT7YjG6bxTkHB9D5LYYYBhgiFyeKIrYcKQUb/90BCfOVwEA4jR+eGVYV/TvHCJxdURkDwwwDDBEbkNvNGPpjlN4f8MxaC8ZAAADOofgyeSOuLNTMASBPTJE7oIBhgGGyO2UV+vx/oZj+CrrFIz16yl1DffHtOQOGNYjHAo5p7YicnUMMAwwRG6r8EI1Pt2WjxW7C3HJUHfHUoTaE4/3j8GYfu3gq/KQuEIiulUMMAwwRG7vYpUey3aewhfbT+F8ZS0AwM/TA+MS2+GxO2KgUXtKXCERNRcDDAMMUatRYzBh9d4z+HjrCZw4VzfYVyEXMKJnW0xL7oBYjZ/EFRJRUzHAMMAQtTpms4iNR0vx8dYT2JV/wbL/ri5t8GRyByR15IBfImfHAMMAQ9Sq7S24iH9tPYF1B4tRP94X3dv6Y+qADhgWHw4PDvglckoMMAwwRATgVFkVPt2Wj5W/FaLGYAYAtA3wwth+UUju0gbdItScGI/IiTDAMMAQ0RUuVunx1Y5T+HL7SZRV6S37A7wVuLNTCAZ0CkH/ziGIDOS6S0RSYoBhgCGiRtQYTPgh5ywyjpRgR14ZKmqNVsdjQnwwoHMI+ncKQVLHYK6/RORgDDAMMER0E0aTGftOlyPz9/PYdvw8cgrLYTJf/l+iXCagV1QABnQOwYDOIegZGcCxM0R2xgDDAENEzaSrMSArrwzbjtUFmvz69Zca+Kk8cHvHYCR3DkHPqAAEeisR6KOEj1LOu5uIbIQBhgGGiFqo8EI1th0/bwk0DeswXU0hFxDgrUSgt8LyM9BbiQBvJQK8FVfst27D3hyiazHAMMAQkQ2ZzCIOntFi2/Hz2HrsHPLPV+FitQF6o/mWXs9TIcMf2gUiMSYYiR2C0CsqAJ4KuY2rJnI9DDAMMERkZ6Io4pLBhIvVBlys0kN7yYCL1XpcrDagvKr+Z7X+8r76n7oaA67+P6/SQ4ZeUQG4PSYIiR2C8Yd2gfBSMtBQ68MAwwBDRE7KZBZx4lwlduRfwM4TZdiZfwHnKmqt2ijkAnpEBiCxPtD0iQ6EDxeppFaAAYYBhohchCiKyD9fhZ1XBJoibY1VG7lMQPe26voemiD0jgpEgLeCg4fJ7TDAMMAQkYsSRRGFFy5hR34Zdp64gJ35ZTh98dI17eQyAQFeivqBwsrLA4h9lJZ9gd4KqL2UCPRpaKOAyoOXpsh5McAwwBCRGzlTfqmud6Y+0Jwsq77l1/JWyhHorUSIrxLBvioE+9T9rHusRLCPCsG+SoT4qhDko4SCd0uRA7lFgFm8eDEWLlyI4uJi9OzZEx9++CH69evXpOcywBCRO6sxmFBe3TBoWA9ttaFuMHG13jJYuOFnw/HySwarifqaSu2lqAs09cEm2FeJIG8lPOQyyGVC3SYIkMkEyIW6niGZTICHTIBMECxtrvzdQybAR+UBf08F/L084O+lgK/SAzKuS9XqNfXvt9OOCFuxYgXS0tKwZMkSJCYmYtGiRUhJSUFubi5CQ0OlLo+ISFKeCjk0ajk0as8mP8dsFlFRa0R5tR4Xquq2sko9zlXWoqxSj7Kqup/nK2tRVn/cZBahvWSA9pIBJ85V3fxNWkAQAF9LqFHAz/OKgOOpgL9nXdDx96w75uvpAW+lB3xVHvBRyeGj9ICPygNKD/YYtQZO2wOTmJiIvn374qOPPgIAmM1mREVFYebMmXj55Zdv+nz2wBARtYy5PryUVdXifKXeEnLOV+pxsUoPo1mEyWyGyQyYRREmswiTKMJsrv+9/rHJLFqOm82A0WyG0SyistYI3SUjKmoMqL3F+XQao5DX9e7UBRq59e/1IcdbJYdMECy1Gq+sseGzmFH3+URc/kz1n89oFiGXCVDIBSg95HU/5TIoPWRQyOs2pYcMSrlg+V0hl1m1qeuVAgSh4Wfd7wIAmSBAEC7/FHDjtpZ2DW1lVz6n4XnWrw3A+r+ZKMJouvo8wOq/5dXnoWdUACICvGz23w5w8R4YvV6P7OxspKenW/bJZDIMGjQIWVlZjT6ntrYWtbWXb0PU6XR2r5OIyJ3JZAICfeqWS+hk547vGoMJFTV1YUZXY4TuUt18ORVX/N4QdhqOV9YaUa03oarWiMpaoyUEGUwiyqsNKK9ufOZksp0Px/a2eYBpKqcMMOfPn4fJZEJYWJjV/rCwMBw9erTR58yfPx+vv/66I8ojIiIb81TI4amQo42f6pZfw2gyo0pvQrXeiKpaI6pqL4ebar2p/qcRlfX7RRGQywC5TFb30zKOp/7n1WN75DLIBQFyWV0vhlkUoTeJMBjNMJjM0Df8NIlXPb78u6H+WK3RDJNZhCiKMIuAiLq7z0SxrjfLLF5+LKLusbnuAcyiCBH1P0VYniNecUy84jWsX6vh9esuvlz5Wa8cs2Q9dgnwkMnq21mfp0BvpS3+898SpwwwtyI9PR1paWmWxzqdDlFRURJWREREjuQhl0HtJYPaSyF1KeQAThlgQkJCIJfLUVJSYrW/pKQEGo2m0eeoVCqoVLee3ImIiMh1OOVQbaVSiYSEBGzYsMGyz2w2Y8OGDUhKSpKwMiIiInIGTtkDAwBpaWmYNGkS+vTpg379+mHRokWoqqrCY489JnVpREREJDGnDTCjR4/GuXPnMHfuXBQXF6NXr15Yt27dNQN7iYiIqPVx2nlgWorzwBAREbmepv79dsoxMEREREQ3wgBDRERELocBhoiIiFwOAwwRERG5HAYYIiIicjkMMERERORyGGCIiIjI5TDAEBERkcthgCEiIiKX47RLCbRUwwTDOp1O4kqIiIioqRr+bt9soQC3DTAVFRUAgKioKIkrISIiouaqqKiAWq2+7nG3XQvJbDbj7Nmz8PPzgyAINntdnU6HqKgoFBYWco2lW8Rz2DI8fy3Hc9gyPH8tx3N4faIooqKiAhEREZDJrj/SxW17YGQyGSIjI+32+v7+/vzStRDPYcvw/LUcz2HL8Py1HM9h427U89KAg3iJiIjI5TDAEBERkcthgGkmlUqF1157DSqVSupSXBbPYcvw/LUcz2HL8Py1HM9hy7ntIF4iIiJyX+yBISIiIpfDAENEREQuhwGGiIiIXA4DDBEREbkcBphmWrx4Mdq3bw9PT08kJiZi165dUpfkMubNmwdBEKy2uLg4qctyWpmZmRg+fDgiIiIgCAJWr15tdVwURcydOxfh4eHw8vLCoEGDcOzYMWmKdVI3O4eTJ0++5js5ZMgQaYp1QvPnz0ffvn3h5+eH0NBQjBw5Erm5uVZtampqMH36dAQHB8PX1xejRo1CSUmJRBU7l6acv7vvvvua7+BTTz0lUcWuhQGmGVasWIG0tDS89tpr2LNnD3r27ImUlBSUlpZKXZrL6NatG4qKiizbtm3bpC7JaVVVVaFnz55YvHhxo8cXLFiADz74AEuWLMHOnTvh4+ODlJQU1NTUOLhS53WzcwgAQ4YMsfpOfvPNNw6s0Llt2bIF06dPx44dO5CRkQGDwYDBgwejqqrK0mbWrFlYs2YNvv32W2zZsgVnz57FQw89JGHVzqMp5w8Apk6davUdXLBggUQVuxiRmqxfv37i9OnTLY9NJpMYEREhzp8/X8KqXMdrr70m9uzZU+oyXBIAcdWqVZbHZrNZ1Gg04sKFCy37ysvLRZVKJX7zzTcSVOj8rj6HoiiKkyZNEh944AFJ6nFFpaWlIgBxy5YtoijWfecUCoX47bffWtocOXJEBCBmZWVJVabTuvr8iaIo3nXXXeJzzz0nXVEujD0wTaTX65GdnY1BgwZZ9slkMgwaNAhZWVkSVuZajh07hoiICHTo0AHjx49HQUGB1CW5pPz8fBQXF1t9H9VqNRITE/l9bKbNmzcjNDQUsbGxePrpp1FWViZ1SU5Lq9UCAIKCggAA2dnZMBgMVt/DuLg4tGvXjt/DRlx9/hosW7YMISEh6N69O9LT01FdXS1FeS7HbRdztLXz58/DZDIhLCzMan9YWBiOHj0qUVWuJTExEV988QViY2NRVFSE119/HQMGDMDBgwfh5+cndXkupbi4GAAa/T42HKObGzJkCB566CHExMQgLy8Pf/rTn5CamoqsrCzI5XKpy3MqZrMZzz//PO688050794dQN33UKlUIiAgwKotv4fXauz8AcC4ceMQHR2NiIgI7N+/H3PmzEFubi6+++47Cat1DQww5DCpqamW33v06IHExERER0dj5cqVmDJlioSVUWs1ZswYy+/x8fHo0aMHOnbsiM2bN2PgwIESVuZ8pk+fjoMHD3Lc2i263vmbNm2a5ff4+HiEh4dj4MCByMvLQ8eOHR1dpkvhJaQmCgkJgVwuv2Z0fUlJCTQajURVubaAgAB06dIFx48fl7oUl9PwneP30bY6dOiAkJAQfievMmPGDKxduxabNm1CZGSkZb9Go4Fer0d5eblVe34PrV3v/DUmMTERAPgdbAIGmCZSKpVISEjAhg0bLPvMZjM2bNiApKQkCStzXZWVlcjLy0N4eLjUpbicmJgYaDQaq++jTqfDzp07+X1sgdOnT6OsrIzfyXqiKGLGjBlYtWoVNm7ciJiYGKvjCQkJUCgUVt/D3NxcFBQU8HuIm5+/xuTk5AAAv4NNwEtIzZCWloZJkyahT58+6NevHxYtWoSqqio89thjUpfmEl588UUMHz4c0dHROHv2LF577TXI5XKMHTtW6tKcUmVlpdW/wvLz85GTk4OgoCC0a9cOzz//PN566y107twZMTExePXVVxEREYGRI0dKV7STudE5DAoKwuuvv45Ro0ZBo9EgLy8PL730Ejp16oSUlBQJq3Ye06dPx9dff43vv/8efn5+lnEtarUaXl5eUKvVmDJlCtLS0hAUFAR/f3/MnDkTSUlJuP322yWuXno3O395eXn4+uuvMXToUAQHB2P//v2YNWsWkpOT0aNHD4mrdwFS3wblaj788EOxXbt2olKpFPv16yfu2LFD6pJcxujRo8Xw8HBRqVSKbdu2FUePHi0eP35c6rKc1qZNm0QA12yTJk0SRbHuVupXX31VDAsLE1UqlThw4EAxNzdX2qKdzI3OYXV1tTh48GCxTZs2okKhEKOjo8WpU6eKxcXFUpftNBo7dwDEzz//3NLm0qVL4jPPPCMGBgaK3t7e4oMPPigWFRVJV7QTudn5KygoEJOTk8WgoCBRpVKJnTp1EmfPni1qtVppC3cRgiiKoiMDExEREVFLcQwMERERuRwGGCIiInI5DDBERETkchhgiIiIyOUwwBAREZHLYYAhIiIil8MAQ0RERC6HAYaIiIhcDgMMETm9u+++G88//zwAoH379li0aJHlWHFxMe677z74+PggICDguvuIyL1wLSQicim7d++Gj4+P5fF7772HoqIi5OTkQK1WX3cfEbkXBhgicilt2rSxepyXl4eEhAR07tz5hvuIyL3wEhIROZWqqipMnDgRvr6+CA8Px9///ner41deQmrfvj3++9//4t///jcEQcDkyZMb3UdE7oc9METkVGbPno0tW7bg+++/R2hoKP70pz9hz5496NWr1zVtd+/ejYkTJ8Lf3x/vv/8+vLy8oNfrr9lHRO6HAYaInEZlZSU+/fRTLF26FAMHDgQAfPnll4iMjGy0fZs2baBSqeDl5QWNRmPZ39g+InIvvIRERE4jLy8Per0eiYmJln1BQUGIjY2VsCoickYMMERERORyGGCIyGl07NgRCoUCO3futOy7ePEifv/9dwmrIiJnxDEwROQ0fH19MWXKFMyePRvBwcEIDQ3Fn//8Z8hk/LcWEVljgCEip7Jw4UJUVlZi+PDh8PPzwwsvvACtVit1WUTkZARRFEWpiyAiIiJqDvbLEhERkcthgCEiIiKXwwBDRERELocBhoiIiFwOAwwRERG5HAYYIiIicjkMMERERORyGGCIiIjI5TDAEBERkcthgCEiIiKXwwBDRERELuf/AdKyd2LWHvkNAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "predictions[\"diff\"].round().value_counts().sort_index().plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 228,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4.770979960120657"
      ]
     },
     "execution_count": 228,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mean_absolute_error(predictions[\"actual\"],predictions[\"prediction\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "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.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}