Forecast.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "2df03549",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "\n",
    "weather = pd.read_csv(\"data.csv\", index_col=\"DATE\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8c563f2a",
   "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",
       "    </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",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1970-01-01</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>34.0</td>\n",
       "      <td>24.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-02</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>20.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-03</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>25.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-04</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>18.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-05</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>16.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",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-09</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>100.0</td>\n",
       "      <td>67.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-10</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>103.0</td>\n",
       "      <td>71.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-11</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>84.0</td>\n",
       "      <td>63.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-12</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.01</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>93.0</td>\n",
       "      <td>60.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-13</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>99.0</td>\n",
       "      <td>61.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>38234 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                STATION                                               NAME  \\\n",
       "DATE                                                                         \n",
       "1970-01-01  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-02  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-03  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-04  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-05  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "...                 ...                                                ...   \n",
       "2022-07-09  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-10  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-11  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-12  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-13  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "\n",
       "            PRCP  SNOW  SNWD   TMAX  TMIN  \n",
       "DATE                                       \n",
       "1970-01-01  0.00   0.0   0.0   34.0  24.0  \n",
       "1970-01-02  0.00   0.0   0.0   39.0  20.0  \n",
       "1970-01-03  0.02   0.0   0.0   36.0  25.0  \n",
       "1970-01-04  0.00   0.0   0.0   31.0  18.0  \n",
       "1970-01-05  0.00   0.0   0.0   43.0  16.0  \n",
       "...          ...   ...   ...    ...   ...  \n",
       "2022-07-09  0.00   NaN   NaN  100.0  67.0  \n",
       "2022-07-10  0.00   NaN   NaN  103.0  71.0  \n",
       "2022-07-11  0.00   NaN   NaN   84.0  63.0  \n",
       "2022-07-12  0.01   NaN   NaN   93.0  60.0  \n",
       "2022-07-13  0.00   NaN   NaN   99.0  61.0  \n",
       "\n",
       "[38234 rows x 7 columns]"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "ae3c0687",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "STATION    0.000000\n",
       "NAME       0.000000\n",
       "PRCP       0.016294\n",
       "SNOW       0.246718\n",
       "SNWD       0.246090\n",
       "TMAX       0.001752\n",
       "TMIN       0.001779\n",
       "dtype: float64"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\n",
    "null_pct = weather.apply(pd.isnull).sum()/weather.shape[0]\n",
    "null_pct"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "9ac07d94",
   "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",
       "    </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",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1970-01-01</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>34.0</td>\n",
       "      <td>24.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-02</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>20.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-03</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>25.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-04</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>18.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-05</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>16.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",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-09</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>100.0</td>\n",
       "      <td>67.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-10</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>103.0</td>\n",
       "      <td>71.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-11</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>84.0</td>\n",
       "      <td>63.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-12</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.01</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>93.0</td>\n",
       "      <td>60.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-13</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>99.0</td>\n",
       "      <td>61.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>38234 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                                               name  \\\n",
       "DATE                                                                         \n",
       "1970-01-01  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-02  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-03  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-04  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-05  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "...                 ...                                                ...   \n",
       "2022-07-09  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-10  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-11  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-12  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-13  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "\n",
       "            prcp  snow  snwd   tmax  tmin  \n",
       "DATE                                       \n",
       "1970-01-01  0.00   0.0   0.0   34.0  24.0  \n",
       "1970-01-02  0.00   0.0   0.0   39.0  20.0  \n",
       "1970-01-03  0.02   0.0   0.0   36.0  25.0  \n",
       "1970-01-04  0.00   0.0   0.0   31.0  18.0  \n",
       "1970-01-05  0.00   0.0   0.0   43.0  16.0  \n",
       "...          ...   ...   ...    ...   ...  \n",
       "2022-07-09  0.00   NaN   NaN  100.0  67.0  \n",
       "2022-07-10  0.00   NaN   NaN  103.0  71.0  \n",
       "2022-07-11  0.00   NaN   NaN   84.0  63.0  \n",
       "2022-07-12  0.01   NaN   NaN   93.0  60.0  \n",
       "2022-07-13  0.00   NaN   NaN   99.0  61.0  \n",
       "\n",
       "[38234 rows x 7 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.columns = weather.columns.str.lower()\n",
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "cfdad7bc",
   "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": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather = weather.ffill()\n",
    "weather.apply(pd.isnull).sum()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "d863596d",
   "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": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.apply(lambda x: (x == 9999).sum())\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "03464151",
   "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": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.dtypes\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "ec925c4e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['1970-01-01', '1970-01-02', '1970-01-03', '1970-01-04', '1970-01-05',\n",
       "       '1970-01-06', '1970-01-07', '1970-01-08', '1970-01-09', '1970-01-10',\n",
       "       ...\n",
       "       '2022-07-04', '2022-07-05', '2022-07-06', '2022-07-07', '2022-07-08',\n",
       "       '2022-07-09', '2022-07-10', '2022-07-11', '2022-07-12', '2022-07-13'],\n",
       "      dtype='object', name='DATE', length=38234)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.index\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "8991fb9a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1970    730\n",
       "1971    730\n",
       "1972    732\n",
       "1973    730\n",
       "1974    730\n",
       "1975    730\n",
       "1976    732\n",
       "1977    730\n",
       "1978    730\n",
       "1979    730\n",
       "1980    732\n",
       "1981    730\n",
       "1982    730\n",
       "1983    730\n",
       "1984    732\n",
       "1985    730\n",
       "1986    730\n",
       "1987    730\n",
       "1988    732\n",
       "1989    730\n",
       "1990    730\n",
       "1991    730\n",
       "1992    732\n",
       "1993    730\n",
       "1994    730\n",
       "1995    730\n",
       "1996    608\n",
       "1997    725\n",
       "1998    728\n",
       "1999    730\n",
       "2000    732\n",
       "2001    730\n",
       "2002    730\n",
       "2003    730\n",
       "2004    732\n",
       "2005    730\n",
       "2006    730\n",
       "2007    730\n",
       "2008    732\n",
       "2009    729\n",
       "2010    730\n",
       "2011    729\n",
       "2012    731\n",
       "2013    730\n",
       "2014    730\n",
       "2015    728\n",
       "2016    732\n",
       "2017    730\n",
       "2018    727\n",
       "2019    730\n",
       "2020    731\n",
       "2021    730\n",
       "2022    388\n",
       "Name: DATE, dtype: int64"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.index = pd.to_datetime(weather.index)\n",
    "weather.index.year.value_counts().sort_index()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "a8cd805c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot: xlabel='DATE'>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAGVCAYAAADUsQqzAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABFXElEQVR4nO3deXgT1cI/8G/adGNpWYSWSpGC/C4iiggXKPB6VaoVRUXxXlAUXBAXcAFX7isgvkoRr4h4kU1WZRMRXMAiFtmk7DuyU9aStixtujfL+f1RG5I0W9NJ5iT9fp6nzwPJZHLmJDnznTNnzmiEEAJEREREEglRuwBERERE9hhQiIiISDoMKERERCQdBhQiIiKSDgMKERERSYcBhYiIiKTDgEJERETSYUAhIiIi6WjVLoA3zGYzsrKyUL9+fWg0GrWLQ0RERB4QQqCgoADx8fEICXHdRxKQASUrKwsJCQlqF4OIiIi8cO7cOTRv3tzlMgEZUOrXrw+gYgOjo6NVLg0RERF5Qq/XIyEhwbIfdyUgA0rlaZ3o6GgGFCIiogDjyfAMDpIlIiIi6TCgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYWIiIikw4BCRERE0mFAISIiIukwoBAREZF0GFCIVKDLL8W/VxzAUV2B2kUhIpISAwqRCl5ZvBuLtp3F/VM2qV0UIiIpMaAQqeDPLD0AwGQWKpeEiEhODChEREQkHQYUIiIikg4DChEREUmHAYWIiIikw4BCRERE0mFAISIiIukwoBAREZF0GFCIiIhIOgwoREREJB0GFCIiIpIOAwoRERFJhwGFiIiIpMOAQkRERNJhQCEiIiLpMKAQERGRdBhQiFSg0WjULgIRkdQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYWIiIikU+2AsnHjRjz44IOIj4+HRqPBypUrbZ4XQmDMmDFo1qwZoqKikJycjOPHj9ssc+XKFQwcOBDR0dFo0KABnnvuORQWFtZoQ4iIiCh4VDugFBUVoUOHDpg6darD5ydOnIgpU6Zg+vTp2LZtG+rWrYuUlBSUlpZalhk4cCAOHTqEtWvX4ueff8bGjRsxdOhQ77eCiIiIgoq2ui/o3bs3evfu7fA5IQQmT56M9957Dw8//DAAYMGCBYiNjcXKlSsxYMAAHD58GGlpadixYwc6d+4MAPjiiy9w//334z//+Q/i4+NrsDlEREQUDBQdg5KZmQmdTofk5GTLYzExMejatSsyMjIAABkZGWjQoIElnABAcnIyQkJCsG3bNofrLSsrg16vt/kjIiKi4KVoQNHpdACA2NhYm8djY2Mtz+l0OjRt2tTmea1Wi0aNGlmWsZeamoqYmBjLX0JCgpLFJiIiIskExFU8o0aNQn5+vuXv3LlzaheJiIiIfEjRgBIXFwcAyM7Otnk8Ozvb8lxcXBxycnJsnjcajbhy5YplGXsRERGIjo62+SMiIqLgpWhASUxMRFxcHNLT0y2P6fV6bNu2DUlJSQCApKQk5OXlYdeuXZZl1q1bB7PZjK5duypZHCIiIgpQ1b6Kp7CwECdOnLD8PzMzE3v37kWjRo3QokULvP766/jwww/Rpk0bJCYmYvTo0YiPj0ffvn0BADfddBPuu+8+PP/885g+fToMBgOGDx+OAQMG8AoeIiIiAuBFQNm5cyfuuusuy/9HjhwJABg8eDDmzZuHt99+G0VFRRg6dCjy8vLQs2dPpKWlITIy0vKahQsXYvjw4ejVqxdCQkLQr18/TJkyRYHNISIiomCgEUIItQtRXXq9HjExMcjPz+d4FApI7ceuQWGZEQBwesIDKpeGiMg/qrP/DoireIiCjUbtAhARSY4BhYiIiKTDgEJERETSYUAhIiIi6TCgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYWIiIikw4BCRERE0mFAISIiIukwoBAREZF0GFCIiIhIOgwoREREJB0GFCIiIpIOAwoRERFJhwGFiIiIpMOAQkRERNJhQCEiIiLpMKAQERGRdBhQiIiISDoMKERERCQdBhQiIiKSDgMKERERSYcBhYiIiKTDgEJERETSYUAhIiIi6TCgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYVIDRq1C0BEJDcGFCIiIpIOAwoRERFJhwGFiIiIpMOAQkRERNJhQCEiIiLpMKAQERGRdBhQiIiISDoMKERERCQdBhQiIiKSjuIBxWQyYfTo0UhMTERUVBRat26N//u//4MQwrKMEAJjxoxBs2bNEBUVheTkZBw/flzpohAREVGAUjygfPzxx5g2bRr++9//4vDhw/j4448xceJEfPHFF5ZlJk6ciClTpmD69OnYtm0b6tati5SUFJSWlipdHCIiIgpAWqVXuGXLFjz88MN44IEHAAAtW7bE4sWLsX37dgAVvSeTJ0/Ge++9h4cffhgAsGDBAsTGxmLlypUYMGCA0kUiIiKiAKN4D0r37t2Rnp6OY8eOAQD27duHzZs3o3fv3gCAzMxM6HQ6JCcnW14TExODrl27IiMjw+E6y8rKoNfrbf6IiIgoeCneg/Luu+9Cr9ejbdu2CA0NhclkwkcffYSBAwcCAHQ6HQAgNjbW5nWxsbGW5+ylpqZi3LhxSheViIiIJKV4D8q3336LhQsXYtGiRdi9ezfmz5+P//znP5g/f77X6xw1ahTy8/Mtf+fOnVOwxERERCQbxXtQ3nrrLbz77ruWsSS33HILzpw5g9TUVAwePBhxcXEAgOzsbDRr1szyuuzsbNx2220O1xkREYGIiAili0qkGo3aBSAikpziPSjFxcUICbFdbWhoKMxmMwAgMTERcXFxSE9Ptzyv1+uxbds2JCUlKV0cIiIiCkCK96A8+OCD+Oijj9CiRQvcfPPN2LNnDyZNmoRnn30WAKDRaPD666/jww8/RJs2bZCYmIjRo0cjPj4effv2Vbo4REREFIAUDyhffPEFRo8ejZdffhk5OTmIj4/HCy+8gDFjxliWefvtt1FUVIShQ4ciLy8PPXv2RFpaGiIjI5UuDhEREQUgjbCe4jVA6PV6xMTEID8/H9HR0WoXh6jabn1/DfSlRgDA6QkPqFwaIiL/qM7+m/fiISIiIukwoBAREZF0GFCIiIhIOgwoREREJB0GFCIiIpIOAwoRERFJhwGFiIiIpMOAQkRERNJhQCEiIiLpMKAQERGRdBhQiIiISDoMKERERCQdBhQiIiKSDgMKERERSYcBhYiIiKTDgEJERETSYUAhIiIi6TCgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYWIiIikw4BCRERE0mFAISIiIukwoBAREZF0GFCIiIhIOgwoREREJB0GFCIVaDQatYtARCQ1BhQiIiKSDgMKERERSYcBhYiIiKTDgEJERETSYUAhIiIi6TCgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYWIiIikw4BCRERE0mFAISIiIukwoBAREZF0fBJQLly4gCeffBKNGzdGVFQUbrnlFuzcudPyvBACY8aMQbNmzRAVFYXk5GQcP37cF0UhIiKiAKR4QLl69Sp69OiBsLAw/PLLL/jzzz/x6aefomHDhpZlJk6ciClTpmD69OnYtm0b6tati5SUFJSWlipdHCIiIgpAWqVX+PHHHyMhIQFz5861PJaYmGj5txACkydPxnvvvYeHH34YALBgwQLExsZi5cqVGDBggNJFIpKORqN2CYiI5KZ4D8qPP/6Izp0745///CeaNm2Kjh07YtasWZbnMzMzodPpkJycbHksJiYGXbt2RUZGhsN1lpWVQa/X2/wRERFR8FI8oJw6dQrTpk1DmzZtsGbNGrz00kt49dVXMX/+fACATqcDAMTGxtq8LjY21vKcvdTUVMTExFj+EhISlC42ERERSUTxgGI2m3H77bdj/Pjx6NixI4YOHYrnn38e06dP93qdo0aNQn5+vuXv3LlzCpaYiIiIZKN4QGnWrBnatWtn89hNN92Es2fPAgDi4uIAANnZ2TbLZGdnW56zFxERgejoaJs/IiIiCl6KB5QePXrg6NGjNo8dO3YMN9xwA4CKAbNxcXFIT0+3PK/X67Ft2zYkJSUpXRwiIiIKQIpfxTNixAh0794d48ePx7/+9S9s374dM2fOxMyZMwEAGo0Gr7/+Oj788EO0adMGiYmJGD16NOLj49G3b1+li0NEREQBSPGA8ve//x0rVqzAqFGj8MEHHyAxMRGTJ0/GwIEDLcu8/fbbKCoqwtChQ5GXl4eePXsiLS0NkZGRSheHiIiIApBGCCHULkR16fV6xMTEID8/n+NRKCDd9sGvyCs2AABOT3hA5dIQEflHdfbfvBcPERERSYcBhYiIiKTDgEJERETSYUAhIiIi6TCgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYWIiIikw4BCRERE0mFAISIiIukwoBAREZF0GFCIiIhIOgwoREREJB0GFCIiIpIOAwoRERFJhwGFiIiIpMOAQkRERNJhQCEiIiLpMKAQERGRdBhQiIiISDoMKERERCQdBhQiIiKSDgMKkQo0aheAiEhyDChEREQkHQYUIiIikg4DChEREUmHAYWIiIikw4BCRERE0mFAISIiIukwoBAREZF0GFCIiIhIOgwoREREJB0GFCIiIpIOAwoRERFJhwGFiIiIpMOAQkRERNJhQCEiIiLpMKAQERGRdBhQiFSg0WjULgIRkdQYUIiIiEg6DChEREQkHZ8HlAkTJkCj0eD111+3PFZaWophw4ahcePGqFevHvr164fs7GxfF4WIiIgChE8Dyo4dOzBjxgzceuutNo+PGDECP/30E5YtW4YNGzYgKysLjz76qC+LQkRERAHEZwGlsLAQAwcOxKxZs9CwYUPL4/n5+Zg9ezYmTZqEu+++G506dcLcuXOxZcsWbN261VfFISIiogDis4AybNgwPPDAA0hOTrZ5fNeuXTAYDDaPt23bFi1atEBGRobDdZWVlUGv19v8ESlp1sZTGDxnO8qMJrWLQkRE8FFAWbJkCXbv3o3U1NQqz+l0OoSHh6NBgwY2j8fGxkKn0zlcX2pqKmJiYix/CQkJvig21WIfrT6MDcdy8f3uC2oXhYiI4IOAcu7cObz22mtYuHAhIiMjFVnnqFGjkJ+fb/k7d+6cIuslsldczh4UIiIZKB5Qdu3ahZycHNx+++3QarXQarXYsGEDpkyZAq1Wi9jYWJSXlyMvL8/mddnZ2YiLi3O4zoiICERHR9v8ERERUfDSKr3CXr164cCBAzaPPfPMM2jbti3eeecdJCQkICwsDOnp6ejXrx8A4OjRozh79iySkpKULg4REREFIMUDSv369dG+fXubx+rWrYvGjRtbHn/uuecwcuRINGrUCNHR0XjllVeQlJSEbt26KV0comoRQqhdBCIigg8Ciic+++wzhISEoF+/figrK0NKSgq+/PJLNYpCREREEvJLQFm/fr3N/yMjIzF16lRMnTrVH29PREREAYb34iEiIiLpMKAQERGRdBhQiIiISDoMKERERCQdBhQiIiKSDgMKERERSYcBhcgK52kjIpIDAwoRERFJhwGFiIj84qtNp7Bw2xm1i0EBQpWp7omIqHbJ0Zfiw1WHAQD9OydAG8rjY3KN3xAiIvK5onKT5d8c6kWeYEAhsiLYdBIRSYEBhYiIiKTDgEJERETSYUAhIiK/4nxD5AkGFCIrbDiJfEOjdgEo4DCgEBERkXQYUIiIiEg6DChEREQkHQYUIhXwfDzVZpxviDzBgEJkhc0mkW9omMqpmhhQiIiISDoMKERERCQdBhQiIiKSDgMKBY1dZ67i10O6Gq2DE7UREclBq3YBiJTSb9oWAMD6N+9Ey+vqqlwaIiKqCfagUNDJyitRuwhERFRDDCgUdMw8TUNEFPAYUCjocBIoIqLAx4BCQacmA10Zboh8LxgGo5vMAgaTWe1iBDUGFAo65mBo/YiCjCaIbvAghEDK5I1ISl3HkOJDDCgUdBhPiMjXTuQU4lJhGU5fKlK7KEGLAYWCDxNKQFp94CLm/pGpdjGIqoXNje9wHhQKOjUZR8KzQ+owmQVeXrgbANA1sTHaxUerXCIiz7DN8B32oFDQYYMReKzP4+86e1XFkhC5xzbGPxhQKOhwHpTAFjxDKcmaJkg/WF755zsMKBR0RAAc3gRrY60E1g0FkgBobgIWAwoFHfagBLZguhyVgh8Diu8woFDQefGbXbhaVK52MchL7EEJTmVGk9pFUAwziX8woFBQ4uWqgYv5JDilH85Ruwg+wTEovsOAQorI1pdKNWGRif2uASuEXShBKVh/kWxqfIfzoJAiuo5PBwDsGX0PGtYNV7k0bDSIyD/Y1vgOe1BIUZmX5ehF8bbNCIQrgIhIXWwn/EPxgJKamoq///3vqF+/Ppo2bYq+ffvi6NGjNsuUlpZi2LBhaNy4MerVq4d+/fohOztb6aKQCmT53Toqx+oDF7Fo21n/F4aqh2d4glKwfqw8I+k7igeUDRs2YNiwYdi6dSvWrl0Lg8GAe++9F0VF146sR4wYgZ9++gnLli3Dhg0bkJWVhUcffVTpolAt5mjg2ssLd+PfKw7g3JViFUpErlgHSrb3RAT4YAxKWlqazf/nzZuHpk2bYteuXbjjjjuQn5+P2bNnY9GiRbj77rsBAHPnzsVNN92ErVu3olu3bkoXifxKli4U50/lFRuQ0Mh/RSH3rAOlhoekQcn6Y5Wlp5Xk5vMxKPn5+QCARo0q9gi7du2CwWBAcnKyZZm2bduiRYsWyMjIcLiOsrIy6PV6mz8iV1y1f2a2jlJjPAlOwTQBH1sQ//BpQDGbzXj99dfRo0cPtG/fHgCg0+kQHh6OBg0a2CwbGxsLnU7ncD2pqamIiYmx/CUkJPiy2FQDsuz7XQ1ic1VEWcpPRFTb+TSgDBs2DAcPHsSSJUtqtJ5Ro0YhPz/f8nfu3DmFSuh/+SUG6EsNahcj6DFoEBEFNp/NgzJ8+HD8/PPP2LhxI5o3b255PC4uDuXl5cjLy7PpRcnOzkZcXJzDdUVERCAiIsJXRfWbcqMZHcb9CgA48VFvaEOD7ypvWXKBLOUgz9gMkg2eMwFEVAOK7yGFEBg+fDhWrFiBdevWITEx0eb5Tp06ISwsDOnp6ZbHjh49irNnzyIpKUnp4kjlitX9YYoNwXNfChm56kHh/k9uDCjBL9Cnh2cPrX8o3oMybNgwLFq0CD/88APq169vGVcSExODqKgoxMTE4LnnnsPIkSPRqFEjREdH45VXXkFSUhKv4AkCsvxwvW0AJSk+1RIbjuXCZDbj7raxahfFp4wmMxZvD845iBiofUfxgDJt2jQAwJ133mnz+Ny5c/H0008DAD777DOEhISgX79+KCsrQ0pKCr788kuliyI1WXbkwcplDwobFJJAmdGEwXO2AwD2jb0XMVFhKpfIdxZvP4tTEt2riwKD4gHFkymAIyMjMXXqVEydOlXpt5dabdgxyjoFtKzloqqC6XJUVwyma9/JwjJjUAeU3Wfz1C6CogL9FFWgCL5RmkRwHUj8sQOcszkTwxbthtFk9vl7UWCy/hYGe4C2/8UF+eaSQhhQ1BKkP1BZNkvtcnzw859Ytf8ifjnoeG4fsqX256WG2jyz6hFd8Ey2WVt6/NTAgOJH/Br7j32D7+kOQOkdRXG5UdkVBinrHoTacCoUqGU7NrtNLTOwZ5HcY0AhRQXCkWBt2QEGqtpyL57a3IMS6Jtb2z4vtTCgUFCyH8TG9oRkFuyDLmtVbxEphgGFFBXsDW11+bJhFkIg/XA2zl8t9tl71EZCCKz9MxtZeSU+fZ9a3YMSRNtbSzr8VOGzqe6JZOLpVRKBFLDSD+dgyIKdAIDTEx5QuTQ1Y13rarf3P+2/iFcX7wHg23oNpp20O/Y78UD6nZF62IOikqD9gUqyWbVhorbtp6+oXQSfUPvz2XLikt/fU5Kfjd/UpnBG3mNA8acg2TEGgtrQ/vHrRIGiyjwoqpSCAg0DCilK1oZH1nLVCBNK0Aj2idqCmdo9fsGMAUUlbI98y/XdjDUQQvhnp+DDxitYr4wI1u1yJdibgypjUNgAkgcYUPyoNjS8JrMsDY/dZcbC+hmBf07PQL9pW2A2O19Odjxy8w01vgOB9L1TQqBvbm37vNTCgEKKOpFTqHYR3NKXGLHzzFXsPpuHS0VlahfHayFBFFCsG/zaGbxq2R6vlm0ueYcBhVBSbsKO01eq9CZ4ylV3rdFkxvbMKygzmmz+7WuyHOF4sq/ddeaqV13eQdUjZx1Q1CuFamT5vvpKUH1X7QTztqmNAYUweO52/HN6Br7afMqr17s6+h2/+gj+NSMD73y3Hx+nVfz7zWX7a1Da6pcJsL2s26aMPt4x5Ba676HpN20LftyX5duCBJDa2YMS3DgPCnmDAUUlMv08t2dWzKexaNtZr17valvm/JEJAFi5NwuzNlX8+yeVd8au9n9Kfy6ncos8KsVP+y5We93ckQcPmdoDco8Byz8YUPzIeofy6yGdT98r/XA2pv5+olqnDrz9ydncidbLdfiabS/PtVLK08x4c4oneMjU4KtRlmA/xWMvmLaXBwq+w4Cikne/P+DT9T83fyc+WXMUm457Pium2ctWQzj5t6z8eQ8Un64/iFpG23oKnu3ylEwBzReqXmasTjkosDCgBDmdvtTjZc1mHxbEzzxt8OXZMVR/pxysu/Egyl0eC/4dtu2HGvSbS4pgQPGj6rS77cakoeW7q2p+R9VqtATeTp4kY+MqS5mW7z7v8PFLVQbPenGKJ4h25NZbX1xuVK0capHl+0qe4eflHwwokiour7gUt9enG2q0nur0EHg7x5rNFTLerUI1VRqaAGp5gvXyxsXbz6ldBL+TpyfPP4JpJtng/BXKgQGlmgwmM85eLvbb+5UYfD9nSCVnY1BKyk0ue3JkbGtc5Y5tpy77tSy+4u8elPNXi1Hqh++jwaTuuUbOJKu8qpcZq+v81WJ8s/WMX77P5D0GlGoaPGc77vjkd6Qfzla7KIrLKXA8Z8ed//kd3Ses82iWWI0k5x1cNfj/+fXYteX8UBZ7Sh09+rOmD17IR8+Pf0fvzzf5ZP2BcCWY0oI9lFircjdjlbf93s824r2VBzEl/bi6BSGXGFCqacvJiqPvb7aeUbkknlGiIcjWVwSXdUcchzK1GxtPOOtCD+SuZn9mwZ/3V8zTknnJ2bwuNWP9KcgScil4VZ5C/+OE51c5WgvcViOwMKCo6NXFe9wuU9P7rZiEwIKM0/gzS1+zFXmoJjv8i/klaPnuKvx7heNLsHefvYol28969B6entMvKLUdkOmPhkepTKTEjvzs5WK0fHeV15P0KcW6TnaduapeQVQSwDnZS3JssBKlYJ72HQYUP7Lfofy4LwtHdFWDg/VVDDW9OfDyXecx5odDuH9KzbvmnTWiSg3wS0pdB6BiRtuDF/KrPP/ol1vw7vcH8McJD8aQ2E9176SIo1cerG4xa0yOprnCHZ/8DgBOQyH5R20bJCsLJYLhsWz5b5AaqBhQVHbGwYDbyu5HJRy8oFzPibPfsrNZWmvi/FXnA5FPXVKuQdipwtG6YmNQgujITaYdtDrjklR4Uz+S9WK5cqN3A7Ktf8Mbj+UqVRyyw4DiwI7TV/CvGRk4lFX1KF5pjnZW9lcxbDru/Q+gXMErIpz3oCjPVQNmf1pGSV+sO4H3fzxU5fHiciMGz9mu+qkQtbgLQ0d1Beg/I8PlFVIfrfoT73zn+xtFBiJf7q8v5pfg8ZlbfX57jeqoac+wUo5mF6CwrGbtibczcJN7DCgOvL5kL7ZnXsG/Vyjb/e+ojXeUH1btt71x3FOztytaDm95MtDUH0f1X/5+wu0yVY7YqrH+eVtOo8iu0Zq/5Qw2HMtV5FSIUs2ZTPOgPDtvB7ZlXkH/mVsdPm82C8zalImlO8/hzGUHA21reRvvy8HaY344hIxTlzH0610+e49AFizTDgQjBhQHLvw158e+c3k+f6/8EoNN41RuNONqcXm11qEvNShdLIf82YPiSpGCp8Ccsd+mwrJrdawvNaCg1ACjl71TSs1t448waDYLFJQaYHDTFZ7r5BJ1R+x7CE1mUeOj2EAjREW9Wv7vZDmjyVwlLFfXlaLqtSf+INMpvZqS6UAh2GjVLkBt9+8VB7Dn7FV88s8OAIDuE9JxqdDzBmXNIR1e+HoXXrqzNd65r62viulSodUpF6UOBO1XU93Ju+yPSGt6hGrdCN36/q8AgKRWjbF4aLdqrafcaLa8Xm2frT3mdplBc7ZjsweXYlYnLNl37/ebtgV7/XAwIJN3lx/A0p3XZsw9lKXH7S0aVlnuvs834UROIfaMvgcN64Z79V4y7D7tf37nrtTwFh4K8qZpsL0sXrGikB32oEhg2a5r92upTjgBYBkvMW39SUXLVB2+uCzUvtHQ5Xt+00NA+V4dR41Qhhddwxdqem8lBX3uwSRVnoQTwH0jbXPHa7sPR7Zw4o8hBdbhBADmbM50uFzl5IjefNcqybgDXbHH8T2qAhHn7fEdBhQv/X401+WVJoFijZuBc9Y9D2UOTk0cOJ+PV6zmczntaHyBF6rbBWzfQ7LVrkGv6T7HWSO09k/3Mwpbn4Jz1pQ52ileKSrHD3svOJ2Oe8IvR9y+t7+EVKORrkn3fra+FD/uy6rxdPjHsgucTjyoBkcDLa2vDnFVuzkFpfhh7wWvr0hR0u9HcnBUV+B2OWch8I8Tl/Dq4j0w+XEUrbN3OnghH298u8/tDVtrOleVK+euFOONb/d5VKc18cPeCxj30yGYZRm9/BcGlBro+fHv1VreVRvu7QydNfltbDlxCS+4GTj3y8FrAWbKuqqDUx/872ab/8/943QNSnTNyj1ZNv+/7OY8+pvLbK8OydaXKXru3Vk9P79gJ3acvuLytcMXXQtw1TnYGjAzA68t2YuJaUerPLfFyxkwfaU630NzDfajKZM34tXFe/DVJsc9Dp6697ONeHbeTofjzPxx9Z49+x125qUiDJpzbXC8q6P0h774A68t2Ysv17sfPO5LBy/k45l5O5AyeaPbZZ3tBgd+tQ0/7svyaBJLpTg7/dvni81Yvvu827L4sgPl+QU7sXz3eTw8dbP7hWvgtSV7MfeP01gr2S1cGFAkcdpHU4i7st/BZGj2tme63vn6iv3I+hy961M8y3dX7TK+UuT5wE13XDVC9pPK2Td43s6TUDkB1KoDWVWe++OkXAHFXQ+KdZ3U5LLMvOKK3iilej8cTZR4NNu3R6uO2Pcq2V/p5OooXffXb8NZb56/BnEec1lv1fvMVx246H4hP3H0fbCZ+8mH9Vv53qUG//SO5VXzAg1fY0CxYjCZ8fZ3+1R5b3fd3hfzHXcz1uT8p9ZHfZPOzqdXh8luJ2a/U/Nk0Ou5q9fqrLr7xH9M/N2my9xVI1SddXvTmGXry/B/P/9p89iO0+7H/XydcRrv/3jIL/cbKnBypYkQAmN+OGhzieuby/ZZTls5u1LO3Zgjb3uiSw0mjPx2r8v1eFpdk349iqkeXPLuCfv3tP9dG51s8IHz18JxIE3HoeY9sOwPdtyWxM0CNW1Gi8qMWLrjLC4XVj2gqs6pUyXINp6GAcXKit0X8O1OdQZvufu9jlyqfHDypHck3Ysj1Q/sdqbesA8k9u3z1lPuy/7SN97P+3C5qByLt1+blM1VI+SPWTJnb8606WXz5LMb/cMhzNtyWtV72/xx4jIWZJzBuiM5lseO6Aowf8tpAMDDU/9w+Lo3l7n+vnvbCzNvy2l8v/tCjdeTU1CKKetO4JM1R52OEaqOKgHF7nnr76I161OsTrdErn0OAHXDVJfx6Tb/d1cWd0Wt6U599MqDeGf5AYfzXYX6OTD4OxC5w4Bi5YoC3Vsms7D8VWcgn7sfibPBp8YanND3ZEbW/GLfz7FiNosqR1TWgUQIUWXQnCdzv9h0i3rRIGZZ9Vq5+t1W52iwJr//onLv5sNw1ruhFFcDGq3nj7Hmbt4Ud4Otve1BuWT3vt6ux2C69kIlBnTaf4fsdxSezDPjy14JpQetBlBnj8+tPlhxOuvPi1VPN1p/Dbyddwmo6JH8eusZt59jqGSJQLLiqMtRWr3koNvNGV1+KdqPXYPW/16N1v9ejTb/+4vN0ZWr9sPdD9bRa99beQDZeu/HWYR60Dep9+G08kDFj+6ezzZggN0MpNajyZ+Ytc3mSiH7531lxoZTltM81TlKsj895Q1Hl427yqL2Vy395ubqIqV2ZhknL6P1v1c7fd5ZvX3l5jTg+atuLsf2svwL7W9V4OV6rH86Snze9mtQ8uxrTVf1wU9/ouMHv7o97eY6xNewED7lunAOJxG0HoNSwwp2ddrXOqh+8mvVwfKe6paajtErD+KZeTtcLsceFImFOGgVFm71/N4rszadqjJL6G4Pu9fd7TAcjVH5phplc6RpdESNXq+EI7oCnMwtwja7UxbWXe+O5oCobj7x9tLWPWcrPr/qNL6uBjx7+vv/OK3qJcSuTkfYB7whC3a6XL+jezR5E1peWuj6NJqvuqi9zaf2v09v12O9XUKB8YvuxqB4s45r6/KiQFbm/JEJfakRMzeecrlcdcZXqTkGxZ6jolSnfDXdqbt6ufUuyV39e8LdgH2OQQkwlTuFS4Vl2OnmclJHH61BoSN9699Ljr5UkXEFdcOrTiTs7oe5PfOKoo2Lsx+3u2rbfCIXJX6Y8r6yHC4Hyf4VfkxmgRV7zuNNFzfEq0kD4O14ibQDOjw7b4clbAGOG+X9VgMuPeVoPSdzr91tOkTBFsa610ypG7S5W4+z77r15+hsHTtPX/H4FgD26/CmB8XRKeULeSU2A2mdKTea8ceJSy7H0+w667rNcXRndmd8GU/yiw3IOHkZW09d9vqqlOqc0jrs4NSMtb3n8hz2PunyS7H3XJ7LWGf9PfNHpiv0cY95dTGgWFm47UyVxyq/E0mp6XhsekaV561/0I72PZ6eN3R7isfq313Gp6PftC0erdcVR1Nn/7iv6iWt1v41IwO/H81xuUwlT8bgeHKayZHF28/hZTdH70qoDB+e9KDM33IaI5bu8+oeTp60Pd5m3aU7z2HdkRw88qXr78zDU/9QZKbbXp9u8OrUmDvfWP0+lTrD5249qw84nsjQ+mvrKKBknLyMx6ZnoOv43zwqR45dkLH/XXiyczrloOeux4R1Ht27KvWXwxj41TaMWLrX6TL7zuW5nPPns9/c3zqhkqPtcXTJtzf6/HcTHp+1FQNmbkXyJO/mZKnOabstJy87besOX9Sj79Q/0C01vcpz3VLT0XfqHy4/H3/fo0qJm6EqiQHFyqncqj/wyiMo60Fx1vKsBpE66g1w9rqq71Oz573RPj66ymP2d1J2JP2wZwGl2IOGsSaDsn4/6vn8Il7Xn6UHxb2Vey+4XcbZepxdRmpNydu6O1vXMYVmrKy8wZ2S57Str75RqhfP3XrSnMy0bL1djnZmW/6ap8bbIOXodLMvVU6waD0xoyObjisz/46jU65K3ZzV+j4/nowhdPQVcNc7a19+ZzcA3emHK+gKy4wug4wSV5mpRdWAMnXqVLRs2RKRkZHo2rUrtm+vepmV2tztFJbsuDYOxNE9Rb53MIGYI+66FJVokK3Xse5ItsPbr1eWQwiBXp+ud7gesyW01fzke3Um+KoJ+yNUT1UO1nRVzsoSetIb5Gw11rcRKHByhZK7gcHvrTzg9LX2PKnW5bscf3c9+dwrA5ev9rMXq3lvJmfchWhHv/+TuYXo+H9rLf93VJfe9gxaXu/HsQBlRs93YI7mOBJCYJKbAZyHsmx7RxzVmSdz++QVl6N7ajpe/HqXV22Dp6fcbvtgrfuFrBg9PBBVmtFkRvuxa9B+7Bqnv8skB703gUK1gLJ06VKMHDkSY8eOxe7du9GhQwekpKQgJ8ezo3N/cfcbmPzbccsPxX6gJwD8anU1xZkrzs/RWgcdR9xN9Q64m8mxYk6KSs/OczyIsnLHsvtsHk466FECrl1N4mxuhkqetLHuGnKlOg08mX7bkbeXV4wn8eQUT5iLARfuAmiZ1aRw/1njuLF3dzT+zdazmPDLkSoz2zpel7MRldf++YaTuUgWZFQ9FWqv8vJ3JXe01qvKL1Hm8nd3E60dclCXvT7dYPN/R3UZVsPrNWsacKpjnt3tKVz1HhSWGatMdLb+WK7D22BYO2BXj47C9ndOArG1j1YdRlZ+KdIO6by6gaL91W5Aze4NVcnZdA++mgyzkvVVlledjLe56oepInxFtYAyadIkPP/883jmmWfQrl07TJ8+HXXq1MGcOXPUKpJDnnx1HV0R4YirAVtK3GQvx80lx57c3LByR+pqnpHKBvmCi0tBhRAehQvZLmvzRmUD52qnUtnN6mywrfWstcdzCh0u48kR48EL+R7tvGvSJJ/KdVw+a5VHlLJdFWCvzM0N9jxp3B0Fx5oGDPtq82Yn6mkPg/19wNydErCvk6sODp7c9fY5OyXijvVvw1X744yv2htnByA1DarRkVUvZLBmfQBQk3tcyUojVLjeq7y8HHXq1MF3332Hvn37Wh4fPHgw8vLy8MMPP9gsX1ZWhrKyaztfvV6PhIQE5OfnIzq66jgKb7V8d5XDx8O1IS7vFBoeGgJo4HSZcG3Fl9Td3Ubdvo8H63G3jKfrqHFZ//phOgtvleUwmsxOewbCtSGAcB0A/bG9StRHWKgGGo0GZrNwOt5Eic/XncrXu6pXpcpR03q1lPUv9svaP+8JR+9X3fexf16jqbojsl7GUTndlUMIUWX8mifrsVnGwWfsyTq0IRqbMTDuyurw+b/aRE/X4WiZam+vB8u4ayc8fh8v6tbdeyhRDiXqzPr5AX9PwAcPt6+yfE3o9XrExMR4tP92Hc985NKlSzCZTIiNjbV5PDY2FkeOVJ3/ITU1FePGjfN5uZ7/n0TMcnCXVHc7A3c9KJ7uTNy+jwfr8cc6PHofBepElu1VYh0VOxvXxwJK1Yk/3sMf9arE5+KJmr6PEK6XCbTfv9Es3J5PrOnv36N1BNB3VZZ1KP0engze9yVVAkp1jRo1CiNHjrT8v7IHRWmvJ/8/9OvUHPnFBtSL1KJOuNYmTV4tKkdkWGjF1OtCIEIbWuVoo7jMCINJoNxkRp3wUESFhdo8X2YwQacvhcEk0CAqDNpQDSK0IYj6a04SIQSKykzIKy6HNjQEoSEa1A0PRZ2Iax+VEAJGk0BhmRFCAFHhIdCGhCDM7ghMX2KEWQiYhUB0ZFiV50vKTQgJ0SA8NARXispRP1KLiLBQyzIGoxkGkxn1IrW4XFiOBnXCYDQJm/WUG824XFgGASBCG4IIbSjqhF/bZrNZ4HJROcxCoMxgRsO6YagfGWZT76a/6qug1ICw0BCEa0NQ12p7zWaBMqMZZUYT9CVGXFcvHGGhttsLVAyAM5rMCAnRoH6E1qbOgIorS0oNJoRoNLhcVI74mEjLMpX1Fa7VIK+4ohx1I0Itn0tlfZSbzJZ11AkPRYhGY1MOg9GM4nITNJqK0xxmIdCobrjNd8BkEsgtLENoiAYmsxnX1YuA1u4I3GA0o6DUiDCtBqEaDcK1IVWWKSk3odRggsksEBNVUaeVZTGazCg1mC1dz0azGQ3r2JZDCIH8EgOMJoGo8FBE2L2H2VxR/oJSI8qMJoSGhCAmKqzKEZjBaEap0YSiMhPqRoSiTpgWoaHX3sdoMqOwzAiTWSBEo4HJLNC4Xrjl9I/BWPG8RlNxe4K6EaFVviPW72Uwmat8ttVR9tfn52o9ecXlNr9LeyXlRmg0GoSFhNhsq31ZtX/1nDleR8XnZ/99r1TZ2xaiQZXP3lp+sQFhoRqH21JxG4mKU0Su1lH5GURHhjncHoPRDLMQCA3ROFyP0WRGiEYDo1lYegvtFZYaUW40Iyrcdb2WGc2IiQpzWm+Gv3ag9r9/T5cx/3UbkhKDCXXDtU6XqRg8rEFUeGiV5yuXKSwzom641ul3oHI5jabq6U4hBC4VlqNehBZCCKffxTKDCeUms9PfROVn7KwMlbcKKSg1OpxaonKZglIjoqNs36Ouk233F1UCynXXXYfQ0FBkZ9tOx52dnY24uLgqy0dERCAiwvezntaN0KJtnPMup+sbRCnyPq2a1FNkPW419HzRhEZ1XD7fLMb5tideV7dG61ZKjT+fatRXTbVo7J86cae5H7c5kCj1W/eHQCor2WreUI52QFaqDJINDw9Hp06dkJ5+7fIns9mM9PR0JCUlqVEkIiIikohqp3hGjhyJwYMHo3PnzujSpQsmT56MoqIiPPPMM2oViYiIiCShWkDp378/cnNzMWbMGOh0Otx2221IS0urMnCWiIiIah9VLjOuqepcpkRERERyqM7+m/fiISIiIukwoBAREZF0GFCIiIhIOgwoREREJB0GFCIiIpIOAwoRERFJhwGFiIiIpBMQNwu0Vzl1i16vV7kkRERE5KnK/bYnU7AFZEApKCgAAJ/c0ZiIiIh8q6CgADExMS6XCciZZM1mM7KyslC/fn2nt+T2ll6vR0JCAs6dO8dZahXCOvUN1qtvsF59g/XqG4FWr0IIFBQUID4+HiEhrkeZBGQPSkhICJo3b+7T94iOjg6IDzuQsE59g/XqG6xX32C9+kYg1au7npNKHCRLRERE0mFAISIiIukwoNiJiIjA2LFjERERoXZRggbr1DdYr77BevUN1qtvBHO9BuQgWSIiIgpu7EEhIiIi6TCgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DCimKF4X5ButVGSUlJWoXISgVFBTYfEf5fVVGbf++1orLjMvLyzF79mw0btwYnTt3RqtWrdQuUlAoLy/HlClTEB0djdtuuw1dunRRu0hBgfWqPIPBgFdffRWnT59GkyZN8PLLL6Nr166K38urtjEYDBg+fDgOHjyIxo0bY+DAgejfv7/axQp4/L5WCPoelO+//x5NmzbF3LlzMXLkSNx///2YO3eu2sUKeKtXr0azZs3w3XffYcqUKXjwwQeRmpqqdrECHutVeTqdDl27dsX+/fvx4IMPYv/+/XjxxRfxySefAKi4+ShVX15eHu6++24cPHgQr7zyCgwGA0aPHo2RI0eqXbSAxu+rFRHEzGazSElJEW+99ZYQQohDhw6JsWPHirCwMLF+/XqVSxfYHnvsMfHSSy8JIYTIysoSs2fPFhqNRsydO1eUlZWpXLrAxXpV3nfffSduvvlmcf78eSGEEHl5eeL9998XkZGR4uDBg0KIiraCqmf9+vWiTZs24sCBA0IIIUpLS8XcuXOFRqMRv/zyi8qlC1z8vl4TdAGloKBAFBUVCSGE2Lt3r6hfv77YunWrzTK9e/cWXbp0sXwBqHpOnjwpmjdvLpYsWWLz+NNPPy1uv/32KvVNnjl16hTrVUEmk0kIIcS0adNEfHy8zXMXL14UycnJokePHmoULSgsX75cREVF2TxmNpvFk08+Kdq3by9KSkpUKllg4ve1qqA6xfPWW28hKSkJly5dAgAkJCRAo9EgKysLQMW5fQCYPn06du3ahbS0NNXKGkjWrl2L/fv3W7oWExMTUV5ejqtXrwK4NpDrk08+wcWLF7F69WpLXZNzJ0+etBlMeMMNN7Bea2jmzJlYtGgRTpw4gZCQiuYtNDQUcXFx2LRpk2W5uLg4vPvuu9ixYwfWrl0LgAM7Xdm+fTsA29ML0dHRSEhIwPLlywFU1J9Go8HYsWNx4sQJy+O16pRENX333Xf47bffcPHiRX5fHVE3Hylj2rRpIjo6WjRv3lxoNBrx+++/CyGE0Ol04p///Kd49NFHLcsaDAYhhBBDhw4VHTp0UKG0gWPu3LkiLi5O3HLLLaJ+/fri5ZdftvQ6vfDCCzb1V15eLoQQYsyYMaJFixbCaDSqUeSAMHv2bNGiRQvRqVMn0bVrV/H1119b6sv+e8l69UxaWppo0qSJuO2228QNN9wg2rRpIz799FMhhBD79+8XN910k5gwYYLNaTKdTiceeugh8dRTT6lVbOmtWLFCxMfHi8aNG4vMzEwhxLU29NSpU6JXr17ixRdfFIWFhUKIil4Ag8EgnnnmGXHHHXeoVWzpLViwQDRt2lR06dJFNGnSRPTo0UMsX75cCCHE7t27Rbt27fh9FQHeg7J582a0atUKH3zwAWbMmIENGzagY8eOOH78OAAgNjYW3bp1w9mzZ7FkyRKb1/bv3x8XL17EsWPH1Ci69GbPno2PPvoIn332GdatW4fp06dj/vz5yM3NBQDce++9KCsrw+effw7g2lHS008/jUuXLmHPnj2qlV1mn3/+OVJTUzFx4kR8/vnnuO+++zB48GDMnDkTQgg8+OCDKCkpYb1W01dffYVHHnkEe/bswdq1a/HCCy/gzTffxE8//YRbbrkFd955J5YvX44tW7ZYXhMbG4uwsDDLkSvZWrhwIcaPH4877rgDN910EyZMmAAA0Gq1EEIgMTERd955J3bv3o0VK1YAAEJCQqDVatGwYUNERESgsLBQzU2QjtFotLQB48ePx6ZNm7By5Uq0bt0aX331FUpKStCxY0f07NkT33//fa3/vgb0lv7444+4//77kZmZiQEDBqBVq1a4fPmy5ZQOADzyyCNITEzE9OnTcfnyZWi1WgDA8ePHUa9ePcTExKhVfCkJIWAymbBu3TokJSVhwIABuO666/DEE08gPj7essPs0aMH7r33XkyaNAkXL1603Op7//79uO6661CvXj01N0NKxcXFWLVqleVSzO7du+P9999Hz549MX78ePz666+45557kJKSwnr1gPirizszMxO//fYbHn30UQBAmzZt8MYbb+Dxxx/HG2+8gUuXLuH999+H0WjEzJkzceHCBcs6SkpK0KhRI1XKLyuTyQQAuPHGG9GrVy98/PHHeOihh7B+/XqsX78eQMVlsADw0ksv4frrr8esWbNw9OhRyzpycnIQHx/P76udoqIi5ObmYvDgwXjmmWcQHh6O7t27o127dtDr9ZZTuOPGjYPBYOD3Vd0OnJqpHFQkxLVux8GDB4vk5GSb5VavXi06deok7rnnHpGRkSHOnDkj+vXrJ5544gl2mTvRsWNHMWTIEKHT6YQQQrzyyivib3/7mxg7dqzYsmWLEKJisGz37t3F7bffLpYsWSJOnDgh+vfvL3r37s0Bcg6UlZWJRo0aiUWLFgkhhKWOHnvsMREfHy+efPJJUVBQII4ePSp69OjBenXi2LFjNlcxlJSUiKZNm4qZM2cKIYSlWzwvL0/UqVNHpKamCiGEWLp0qfif//kfccMNN4hPP/1UPPXUU6Jp06Zi06ZN/t8ICdnXqxDX2tWDBw+Khx56SNx///1Vntu0aZPo3bu3aNCggXjzzTfFwIEDRaNGjcTPP/8shKg9V5w4Y1+ve/bssex3KvdhCxcuFLfddpvNKZ1ly5bV+u9rwASUb7/9VgwZMkRMnjxZ7N+/3/K4dUgRomJsxF133SXy8/Mtz5lMJrFv3z7Rvn170bZtW9GkSRPRs2dPcfbsWb9ug4yc1evixYtFQkKCuOeee0Tjxo1F27ZtxQcffCDuuusuceutt4oJEyYIISrOi6akpIh27dqJ+Ph40b17d8u56trMWb0+/vjjom3btpaxPN9884246667xJAhQ8SNN95oWVan04n77ruP9Wpl6dKlomXLluJvf/ub6NKli5g9e7YQQojCwkIxaNAgkZKSYmngK8fujBo1SrRo0cKyjvPnz4uhQ4eKvn37ivvvv18cOXLE/xsiGWf1KoRtuJgzZ45o166dmDNnjhDiWkARouIS4//93/8VgwYNEo8++ijrVVSt16+++srmeet91xNPPCGefvppIYSwCSm1/fsqfUC5dOmSeOyxx0RcXJx48cUXRc+ePcX1118v5s2bZ1nGbDZbPux58+aJ6OhoS0I1m82WH1l+fr44fvy42Llzp/83RDLO6nXu3LmWZXQ6nZg4caK44447hF6vtzz+/PPPi0ceeURkZ2cLISqOYC9evGizI66tnNXr/PnzhRAVR1OtWrUSrVq1EvHx8aJOnTqWwXFarVasWrXKsq7S0lLW619+/fVX0bJlSzF16lSRlpYmRo4cKbRaraXXZN68eaJjx45ixowZQohrO88dO3aIJk2aiB07dtisjz1RFRzVa1hYmJg5c6YoLi4WQlyry/Pnz4vnnntO/P3vfxcFBQVCCFFlbh72SFdwVa+V373KfVNJSYm49dZbxddff+10fbX1+yp9QFm2bFmVOUv69esnWrduLVasWCGEsE2iv/32m0hISBDp6elV1lXbuxqtuarXyh2mwWAQAwYMEB9++KEQ4lpjNHLkSNG6dWvLyH3W6zXO6jUxMdHyfT137pxYs2aNmD9/vuVIPycnR7Rq1UosW7ZMjWJLq/K7NW7cONGpUydLfQkhxMsvvyw6duwo1qxZI/R6vRg4cGCVnqalS5eK+Ph4cerUKX8XXWru6rVz587i+++/r/K6n3/+WXTu3FmMHTtW7Nu3T/Tp04c90Va8qdcLFy6Ili1bimPHjgkhKg5iRowY4b9CS0z6QbKLFi1C8+bNcf3111tGhD/00EM4deoUvvjiC1y6dAkhISGWgV2NGjVCeXm55f/Watt9DFxxVa9Tp05FdnY2tFotLl++jJ07dwIAwsPDkZ2djWPHjmHAgAGoW7cuANarNWf1evr0aXzxxRfIyclB8+bNkZycjEGDBiEsLAwA8PvvvyM8PBw9e/ZUs/jSqfxu/fnnn2jdujXCwsIsAzQ//PBD1K1bF9988w1CQ0MxbNgwhISEYMCAAdiyZQvOnj2L1atXo1OnToiLi1NzM6Tjrl4jIyPxww8/QKfTAbg2cPauu+5Cly5d8MEHH6BTp04wGAxo2rSpOhshoerWKwD89ttvSEhIQLNmzfDaa6+hXbt2OHPmDAwGQ+2Y68QFqQLKxo0bsWbNGhiNRstjbdq0waFDhwDAMiL88OHDuPvuu1FaWmq5vC00NBQA0LFjR5jNZvzxxx9+Lr28vKnXH374AQAwatQorFq1Cj169MDLL7+Mzp07Q6/XY+jQof7fEMnUpF5DQkKQm5uLI0eO4L///S9GjBiBRx99FNddd12tbpTWrl2LV199FZMnT7ZMDgYAvXr1wi+//AKTyWRp9Bs2bIhBgwYhIyMDe/bsQVJSEr766isYjUY8++yz6NKlC/bv349PPvkEUVFRKm6V+ryt18orc0JDQ1FUVISZM2dixowZ+Mc//oHdu3cjLS3NcqVZbeRtvR45cgRAxZVoP//8Mw4ePIiWLVsiPT0dGRkZWL58OcLCwnjwp3IPjhBCiNzcXDFo0CCh0WhEhw4dbLpoT548KZo0aSLuuOMOMXHiRJGUlCQSExNFenq66NChgxg9erQQ4lrXWm5urnjppZfEunXr1NgUqdSkXt977z3LsitWrBDvvPOOeOKJJ8S3336rwpbIRYnvqxBC7Nq1S/Tt21ckJia6PP9cG2RlZYk+ffqIpk2bioEDB4pbbrlFxMTEiG3btgkhhDh69Ki4/vrrLfVnPfYhLi5OTJo0yfL/goICkZmZyVsDiJrX62effWb5/6FDh0TXrl3FggUL/LoNMlKqXouKikSfPn0c3uKCJBiDYjAYxJdffilSUlLE0qVLLZcFlpaWWpbZvHmzGDJkiLj99tvF8OHDRW5urhBCiKeeekr069dPraJLjfXqG0rX6+7du/1afhkVFRWJwYMHi/79+9uMFenSpYvlyga9Xi8+/PBDERUVZRnzUHlQ8o9//EMMGTLE8jqOiaqgdL1SBaXrlRdtOKd6QBFCiK1bt4off/xRCFExuKhJkyZiz549VZazTqHZ2dmiffv2lgGc9pcbE+vVV5SoV+tLNKliiv/KO+BW1s37778vunbtamnYT506JXr06CG6desmTp8+LYQQ4syZM+Kmm26yzLlBtlivvsF69Q8pAor9EU98fLwYOnSo5dJW+0mZysvLxZdffik6duzISzBdYL36ButVedZXO1SG4ieeeEI8//zzNsudP39e3HjjjaJly5aWCe7uvvtuy4SCZIv16husV/+QIqBUqjzi/Pbbb4VWqxW//vqrzfPnz58XX375pejcubPNjJzkGuvVN1ivvtWjRw/LfEcmk8myIzh+/LhYsmSJGDFihM18SOQZ1qtvsF6VpxFCzksGunfvjrp162LhwoVo2rQpcnNz0aRJEyxevBhZWVl444031C5iQGK9+gbrVVmnTp1C9+7dsWrVKnTq1AkAUF5ejvDwcJVLFthYr77BevUNrdoFsGc0GqHVajFr1ix06NABS5YswcmTJ7F582bMnz8fjz/+uNpFDEisV99gvSpLCAGNRoPNmzejXr16lsZ+3Lhx0Ol0GDduHOfd8ALr1TdYr74lbQ8KAHTp0gU7d+5EixYtMGPGDKSkpKhdpKDAevUN1qtyhg8fjrp16yI5ORlDhw5FcXExvv76a9x7771qFy2gsV59g/XqIyqeXnLqxIkTon379qJOnTpVbrBE3mO9+gbrVVklJSXixhtvFBqNRkRERFhuTEk1w3r1Ddar70h3igeomLWwX79+eOedd2r9DJBKYr36ButVWZGRkWjZsiXuueceTJo0CZGRkWoXKSiwXn2D9eo7Up/iIaLayWQyWW5fQcphvfoG69U3GFCIiIhIOlLdLJCIiIgIYEAhIiIiCTGgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYWIFPX0009Do9FAo9EgLCwMsbGxuOeeezBnzhyYzeYqy6ekpCA0NBQ7duwAAJw+fdryemd/8+bNw/r1650+r9Pp/L3ZRKQwBhQiUtx9992Hixcv4vTp0/jll19w11134bXXXkOfPn1gNBoty509exZbtmzB8OHDMWfOHABAQkICLl68aPl74403cPPNN9s81r9/f8s6jh49avPcxYsXeQdZoiAg5b14iCiwRUREIC4uDgBw/fXX4/bbb0e3bt3Qq1cvzJs3D0OGDAEAzJ07F3369MFLL72Ebt26YdKkSYiKirK8FgDq1asHrVZr85i1pk2bokGDBj7fJiLyL/agEJFf3H333ejQoQO+//57AIAQAnPnzsWTTz6Jtm3b4sYbb8R3332ncimJSBYMKETkN23btsXp06cBAL/99huKi4uRkpICAHjyyScxe/bsaq+zefPmqFevnuXv5ptvVrLIRKQSnuIhIr8RQkCj0QAA5syZg/79+0OrrWiGHn/8cbz11ls4efIkWrdu7fE6N23ahPr161v+HxYWpmyhiUgVDChE5DeHDx9GYmIirly5ghUrVsBgMGDatGmW500mE+bMmYOPPvrI43UmJiZyDApREOIpHiLyi3Xr1uHAgQPo168fFi5ciObNm2Pfvn3Yu3ev5e/TTz/FvHnzYDKZ1C4uEamMPShEpLiysjLodDqYTCZkZ2cjLS0Nqamp6NOnDwYNGoROnTrhscceQ/v27W1el5CQgFGjRiEtLQ0PPPCAR++Vk5OD0tJSm8caN27MUz1EAY4BhYgUl5aWhmbNmkGr1aJhw4bo0KEDpkyZgsGDB2PPnj3Yt28fZs2aVeV1MTEx6NWrF2bPnu1xQPnb3/5W5bGMjAx069atxttBROrRCCGE2oUgIiIissYxKERERCQdBhQiIiKSDgMKERERSYcBhYiIiKTDgEJERETSYUAhIiIi6TCgEBERkXQYUIiIiEg6DChEREQkHQYUIiIikg4DChEREUmHAYWIiIik8/8BXXYdeIOXmEcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "weather[\"snwd\"].plot()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "6cc3c220",
   "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>1970-01-01</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>34.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>39.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-02</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>36.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-03</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>25.0</td>\n",
       "      <td>31.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-04</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>18.0</td>\n",
       "      <td>43.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-05</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>40.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>2022-07-09</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>67.0</td>\n",
       "      <td>103.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-10</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>103.0</td>\n",
       "      <td>71.0</td>\n",
       "      <td>84.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-11</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>84.0</td>\n",
       "      <td>63.0</td>\n",
       "      <td>93.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-12</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.01</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>93.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>99.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-13</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>99.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>38234 rows × 8 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                                               name  \\\n",
       "DATE                                                                         \n",
       "1970-01-01  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-02  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-03  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-04  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-05  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "...                 ...                                                ...   \n",
       "2022-07-09  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-10  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-11  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-12  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-13  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "\n",
       "            prcp  snow  snwd   tmax  tmin  target  \n",
       "DATE                                               \n",
       "1970-01-01  0.00   0.0   0.0   34.0  24.0    39.0  \n",
       "1970-01-02  0.00   0.0   0.0   39.0  20.0    36.0  \n",
       "1970-01-03  0.02   0.0   0.0   36.0  25.0    31.0  \n",
       "1970-01-04  0.00   0.0   0.0   31.0  18.0    43.0  \n",
       "1970-01-05  0.00   0.0   0.0   43.0  16.0    40.0  \n",
       "...          ...   ...   ...    ...   ...     ...  \n",
       "2022-07-09  0.00   0.0   0.0  100.0  67.0   103.0  \n",
       "2022-07-10  0.00   0.0   0.0  103.0  71.0    84.0  \n",
       "2022-07-11  0.00   0.0   0.0   84.0  63.0    93.0  \n",
       "2022-07-12  0.01   0.0   0.0   93.0  60.0    99.0  \n",
       "2022-07-13  0.00   0.0   0.0   99.0  61.0     NaN  \n",
       "\n",
       "[38234 rows x 8 columns]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather[\"target\"] = weather.shift(-1)[\"tmax\"]\n",
    "weather\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "93b6dd18",
   "metadata": {},
   "outputs": [],
   "source": [
    "weather = weather.ffill()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "66eea269",
   "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>1970-01-01</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>34.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>39.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-02</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>39.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>36.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-03</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>25.0</td>\n",
       "      <td>31.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-04</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.0</td>\n",
       "      <td>18.0</td>\n",
       "      <td>43.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970-01-05</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>40.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>2022-07-09</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>67.0</td>\n",
       "      <td>103.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-10</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>103.0</td>\n",
       "      <td>71.0</td>\n",
       "      <td>84.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-11</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>84.0</td>\n",
       "      <td>63.0</td>\n",
       "      <td>93.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-12</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.01</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>93.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>99.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-13</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>99.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>99.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>38234 rows × 8 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                                               name  \\\n",
       "DATE                                                                         \n",
       "1970-01-01  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-02  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-03  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-04  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1970-01-05  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "...                 ...                                                ...   \n",
       "2022-07-09  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-10  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-11  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-12  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "2022-07-13  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "\n",
       "            prcp  snow  snwd   tmax  tmin  target  \n",
       "DATE                                               \n",
       "1970-01-01  0.00   0.0   0.0   34.0  24.0    39.0  \n",
       "1970-01-02  0.00   0.0   0.0   39.0  20.0    36.0  \n",
       "1970-01-03  0.02   0.0   0.0   36.0  25.0    31.0  \n",
       "1970-01-04  0.00   0.0   0.0   31.0  18.0    43.0  \n",
       "1970-01-05  0.00   0.0   0.0   43.0  16.0    40.0  \n",
       "...          ...   ...   ...    ...   ...     ...  \n",
       "2022-07-09  0.00   0.0   0.0  100.0  67.0   103.0  \n",
       "2022-07-10  0.00   0.0   0.0  103.0  71.0    84.0  \n",
       "2022-07-11  0.00   0.0   0.0   84.0  63.0    93.0  \n",
       "2022-07-12  0.01   0.0   0.0   93.0  60.0    99.0  \n",
       "2022-07-13  0.00   0.0   0.0   99.0  61.0    99.0  \n",
       "\n",
       "[38234 rows x 8 columns]"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "c5f8bc5f",
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.linear_model import Ridge\n",
    "\n",
    "rr = Ridge(alpha=.1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "4e9b2eba",
   "metadata": {},
   "outputs": [],
   "source": [
    "predictors = weather.columns[~weather.columns.isin([\"target\", \"name\", \"station\"])]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "b7837a22",
   "metadata": {},
   "outputs": [],
   "source": [
    "def backtest(weather, model, predictors, start=3650, step=90):\n",
    "    all_predictions = []\n",
    "    \n",
    "    for i in range(start, weather.shape[0], step):\n",
    "        train = weather.iloc[:i,:]\n",
    "        test = weather.iloc[i:(i+step),:]\n",
    "        \n",
    "        model.fit(train[predictors], train[\"target\"])\n",
    "        \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",
    "        \n",
    "        all_predictions.append(combined)\n",
    "    return pd.concat(all_predictions)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "a0b1f643",
   "metadata": {},
   "outputs": [],
   "source": [
    "predictions = backtest(weather, rr, predictors)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "b24083d4",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "7.111283462839401"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from sklearn.metrics import mean_absolute_error, mean_squared_error\n",
    "\n",
    "mean_absolute_error(predictions[\"actual\"], predictions[\"prediction\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "6a72adc0",
   "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>2000-04-14</th>\n",
       "      <td>27.0</td>\n",
       "      <td>96.375971</td>\n",
       "      <td>69.375971</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-13</th>\n",
       "      <td>27.0</td>\n",
       "      <td>89.923609</td>\n",
       "      <td>62.923609</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1996-02-29</th>\n",
       "      <td>92.0</td>\n",
       "      <td>39.635710</td>\n",
       "      <td>52.364290</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-11-12</th>\n",
       "      <td>73.0</td>\n",
       "      <td>23.967113</td>\n",
       "      <td>49.032887</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-04-13</th>\n",
       "      <td>122.0</td>\n",
       "      <td>74.424496</td>\n",
       "      <td>47.575504</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1986-03-04</th>\n",
       "      <td>48.0</td>\n",
       "      <td>48.002980</td>\n",
       "      <td>0.002980</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-02-12</th>\n",
       "      <td>41.0</td>\n",
       "      <td>41.001795</td>\n",
       "      <td>0.001795</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1997-05-13</th>\n",
       "      <td>73.0</td>\n",
       "      <td>73.000578</td>\n",
       "      <td>0.000578</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2006-08-26</th>\n",
       "      <td>87.0</td>\n",
       "      <td>87.000560</td>\n",
       "      <td>0.000560</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1991-06-11</th>\n",
       "      <td>88.0</td>\n",
       "      <td>87.999882</td>\n",
       "      <td>0.000118</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>34584 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            actual  prediction       diff\n",
       "DATE                                     \n",
       "2000-04-14    27.0   96.375971  69.375971\n",
       "2022-07-13    27.0   89.923609  62.923609\n",
       "1996-02-29    92.0   39.635710  52.364290\n",
       "2000-11-12    73.0   23.967113  49.032887\n",
       "2000-04-13   122.0   74.424496  47.575504\n",
       "...            ...         ...        ...\n",
       "1986-03-04    48.0   48.002980   0.002980\n",
       "1980-02-12    41.0   41.001795   0.001795\n",
       "1997-05-13    73.0   73.000578   0.000578\n",
       "2006-08-26    87.0   87.000560   0.000560\n",
       "1991-06-11    88.0   87.999882   0.000118\n",
       "\n",
       "[34584 rows x 3 columns]"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predictions.sort_values(\"diff\", ascending=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "005f2704",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "prcp   -2.810265\n",
       "snow   -0.496113\n",
       "snwd   -0.002275\n",
       "tmax    0.667245\n",
       "tmin    0.266117\n",
       "dtype: float64"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.Series(rr.coef_, index=predictors)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "6226b172",
   "metadata": {},
   "outputs": [],
   "source": [
    "def pct_diff(old, new):\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[label], weather[col])\n",
    "    return weather\n",
    "    \n",
    "rolling_horizons = [3, 14]\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": 26,
   "id": "00130239",
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "def expand_mean(df):\n",
    "    return df.expanding(1).mean()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "id": "6a99937b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "7.111283462839401"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mean_absolute_error(predictions[\"actual\"], predictions[\"prediction\"])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "a08755da",
   "metadata": {},
   "outputs": [],
   "source": [
    "weather = weather.iloc[14:,:]\n",
    "weather = weather.fillna(0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "6662da5e",
   "metadata": {},
   "outputs": [],
   "source": [
    "predictors = weather.columns[~weather.columns.isin([\"target\", \"name\", \"station\"])]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "304d6dd9",
   "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>2000-04-14</th>\n",
       "      <td>27.0</td>\n",
       "      <td>96.375971</td>\n",
       "      <td>69.375971</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-13</th>\n",
       "      <td>27.0</td>\n",
       "      <td>89.923609</td>\n",
       "      <td>62.923609</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1996-02-29</th>\n",
       "      <td>92.0</td>\n",
       "      <td>39.635710</td>\n",
       "      <td>52.364290</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-11-12</th>\n",
       "      <td>73.0</td>\n",
       "      <td>23.967113</td>\n",
       "      <td>49.032887</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-04-13</th>\n",
       "      <td>122.0</td>\n",
       "      <td>74.424496</td>\n",
       "      <td>47.575504</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1986-03-04</th>\n",
       "      <td>48.0</td>\n",
       "      <td>48.002980</td>\n",
       "      <td>0.002980</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-02-12</th>\n",
       "      <td>41.0</td>\n",
       "      <td>41.001795</td>\n",
       "      <td>0.001795</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1997-05-13</th>\n",
       "      <td>73.0</td>\n",
       "      <td>73.000578</td>\n",
       "      <td>0.000578</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2006-08-26</th>\n",
       "      <td>87.0</td>\n",
       "      <td>87.000560</td>\n",
       "      <td>0.000560</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1991-06-11</th>\n",
       "      <td>88.0</td>\n",
       "      <td>87.999882</td>\n",
       "      <td>0.000118</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>34584 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            actual  prediction       diff\n",
       "DATE                                     \n",
       "2000-04-14    27.0   96.375971  69.375971\n",
       "2022-07-13    27.0   89.923609  62.923609\n",
       "1996-02-29    92.0   39.635710  52.364290\n",
       "2000-11-12    73.0   23.967113  49.032887\n",
       "2000-04-13   122.0   74.424496  47.575504\n",
       "...            ...         ...        ...\n",
       "1986-03-04    48.0   48.002980   0.002980\n",
       "1980-02-12    41.0   41.001795   0.001795\n",
       "1997-05-13    73.0   73.000578   0.000578\n",
       "2006-08-26    87.0   87.000560   0.000560\n",
       "1991-06-11    88.0   87.999882   0.000118\n",
       "\n",
       "[34584 rows x 3 columns]"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predictions.sort_values(\"diff\", ascending=False)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "97d33abc",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "82.5510951730512"
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mean_squared_error(predictions[\"actual\"], predictions[\"prediction\"])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "66f52c9a",
   "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>1990-03-07</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>37.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>48.0</td>\n",
       "      <td>44.666667</td>\n",
       "      <td>-0.171642</td>\n",
       "      <td>24.666667</td>\n",
       "      <td>-0.189189</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>46.357143</td>\n",
       "      <td>-0.201849</td>\n",
       "      <td>26.928571</td>\n",
       "      <td>-0.257294</td>\n",
       "      <td>0.050714</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-08</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>48.0</td>\n",
       "      <td>18.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>43.666667</td>\n",
       "      <td>0.099237</td>\n",
       "      <td>22.666667</td>\n",
       "      <td>-0.205882</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>45.500000</td>\n",
       "      <td>0.054945</td>\n",
       "      <td>25.928571</td>\n",
       "      <td>-0.305785</td>\n",
       "      <td>0.031429</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-09</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>62.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>49.000000</td>\n",
       "      <td>0.265306</td>\n",
       "      <td>23.333333</td>\n",
       "      <td>0.371429</td>\n",
       "      <td>0.006667</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>45.285714</td>\n",
       "      <td>0.369085</td>\n",
       "      <td>25.000000</td>\n",
       "      <td>0.280000</td>\n",
       "      <td>0.012857</td>\n",
       "      <td>0.555556</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-10</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>68.0</td>\n",
       "      <td>57.000000</td>\n",
       "      <td>0.070175</td>\n",
       "      <td>31.000000</td>\n",
       "      <td>0.387097</td>\n",
       "      <td>0.006667</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>46.357143</td>\n",
       "      <td>0.315871</td>\n",
       "      <td>25.857143</td>\n",
       "      <td>0.662983</td>\n",
       "      <td>0.005714</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-11</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>68.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>86.0</td>\n",
       "      <td>63.666667</td>\n",
       "      <td>0.068063</td>\n",
       "      <td>39.333333</td>\n",
       "      <td>0.093220</td>\n",
       "      <td>0.006667</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>49.000000</td>\n",
       "      <td>0.387755</td>\n",
       "      <td>27.714286</td>\n",
       "      <td>0.551546</td>\n",
       "      <td>0.005714</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-12</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>86.0</td>\n",
       "      <td>49.0</td>\n",
       "      <td>85.0</td>\n",
       "      <td>71.666667</td>\n",
       "      <td>0.200000</td>\n",
       "      <td>45.000000</td>\n",
       "      <td>0.088889</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>53.000000</td>\n",
       "      <td>0.622642</td>\n",
       "      <td>30.285714</td>\n",
       "      <td>0.617925</td>\n",
       "      <td>0.005714</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-13</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>85.0</td>\n",
       "      <td>60.0</td>\n",
       "      <td>81.0</td>\n",
       "      <td>79.666667</td>\n",
       "      <td>0.066946</td>\n",
       "      <td>50.666667</td>\n",
       "      <td>0.184211</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>56.357143</td>\n",
       "      <td>0.508238</td>\n",
       "      <td>33.357143</td>\n",
       "      <td>0.798715</td>\n",
       "      <td>0.005714</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-14</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>81.0</td>\n",
       "      <td>54.0</td>\n",
       "      <td>82.0</td>\n",
       "      <td>84.000000</td>\n",
       "      <td>-0.035714</td>\n",
       "      <td>54.333333</td>\n",
       "      <td>-0.006135</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>58.714286</td>\n",
       "      <td>0.379562</td>\n",
       "      <td>35.428571</td>\n",
       "      <td>0.524194</td>\n",
       "      <td>0.005714</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-15</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>82.0</td>\n",
       "      <td>52.0</td>\n",
       "      <td>79.0</td>\n",
       "      <td>82.666667</td>\n",
       "      <td>-0.008065</td>\n",
       "      <td>55.333333</td>\n",
       "      <td>-0.060241</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>61.285714</td>\n",
       "      <td>0.337995</td>\n",
       "      <td>37.142857</td>\n",
       "      <td>0.400000</td>\n",
       "      <td>0.005714</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-16</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>79.0</td>\n",
       "      <td>57.0</td>\n",
       "      <td>72.0</td>\n",
       "      <td>80.666667</td>\n",
       "      <td>-0.020661</td>\n",
       "      <td>54.333333</td>\n",
       "      <td>0.049080</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>63.357143</td>\n",
       "      <td>0.246900</td>\n",
       "      <td>39.214286</td>\n",
       "      <td>0.453552</td>\n",
       "      <td>0.002143</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-17</th>\n",
       "      <td>USW00093721</td>\n",
       "      <td>BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US</td>\n",
       "      <td>0.88</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>72.0</td>\n",
       "      <td>52.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>77.666667</td>\n",
       "      <td>-0.072961</td>\n",
       "      <td>53.666667</td>\n",
       "      <td>-0.031056</td>\n",
       "      <td>0.293333</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>64.214286</td>\n",
       "      <td>0.121246</td>\n",
       "      <td>40.071429</td>\n",
       "      <td>0.297683</td>\n",
       "      <td>0.064286</td>\n",
       "      <td>12.688889</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-07</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.08</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>32.0</td>\n",
       "      <td>23.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>43.000000</td>\n",
       "      <td>-0.255814</td>\n",
       "      <td>26.000000</td>\n",
       "      <td>-0.115385</td>\n",
       "      <td>0.450000</td>\n",
       "      <td>-0.822222</td>\n",
       "      <td>49.214286</td>\n",
       "      <td>-0.349782</td>\n",
       "      <td>26.214286</td>\n",
       "      <td>-0.122616</td>\n",
       "      <td>0.096429</td>\n",
       "      <td>-0.170370</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-08</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.3</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>27.0</td>\n",
       "      <td>56.0</td>\n",
       "      <td>45.333333</td>\n",
       "      <td>0.345588</td>\n",
       "      <td>27.333333</td>\n",
       "      <td>-0.012195</td>\n",
       "      <td>0.456667</td>\n",
       "      <td>-0.956204</td>\n",
       "      <td>50.571429</td>\n",
       "      <td>0.206215</td>\n",
       "      <td>26.357143</td>\n",
       "      <td>0.024390</td>\n",
       "      <td>0.097857</td>\n",
       "      <td>-0.795620</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-09</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>56.0</td>\n",
       "      <td>37.0</td>\n",
       "      <td>63.0</td>\n",
       "      <td>49.666667</td>\n",
       "      <td>0.127517</td>\n",
       "      <td>29.000000</td>\n",
       "      <td>0.275862</td>\n",
       "      <td>0.033333</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>50.785714</td>\n",
       "      <td>0.102672</td>\n",
       "      <td>26.857143</td>\n",
       "      <td>0.377660</td>\n",
       "      <td>0.097857</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-10</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>63.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>60.000000</td>\n",
       "      <td>0.050000</td>\n",
       "      <td>32.333333</td>\n",
       "      <td>0.020619</td>\n",
       "      <td>0.006667</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>51.285714</td>\n",
       "      <td>0.228412</td>\n",
       "      <td>27.285714</td>\n",
       "      <td>0.209424</td>\n",
       "      <td>0.097857</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-11</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>59.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>59.333333</td>\n",
       "      <td>-0.005618</td>\n",
       "      <td>35.333333</td>\n",
       "      <td>0.018868</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>51.928571</td>\n",
       "      <td>0.136176</td>\n",
       "      <td>28.142857</td>\n",
       "      <td>0.279188</td>\n",
       "      <td>0.097857</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-12</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>61.0</td>\n",
       "      <td>33.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>61.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>34.000000</td>\n",
       "      <td>-0.029412</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>52.285714</td>\n",
       "      <td>0.166667</td>\n",
       "      <td>28.500000</td>\n",
       "      <td>0.157895</td>\n",
       "      <td>0.097857</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-13</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.06</td>\n",
       "      <td>0.3</td>\n",
       "      <td>0.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>25.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>52.000000</td>\n",
       "      <td>-0.307692</td>\n",
       "      <td>31.333333</td>\n",
       "      <td>-0.202128</td>\n",
       "      <td>0.020000</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>52.000000</td>\n",
       "      <td>-0.307692</td>\n",
       "      <td>28.428571</td>\n",
       "      <td>-0.120603</td>\n",
       "      <td>0.102143</td>\n",
       "      <td>-0.412587</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-14</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.01</td>\n",
       "      <td>0.1</td>\n",
       "      <td>0.0</td>\n",
       "      <td>40.0</td>\n",
       "      <td>15.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>45.666667</td>\n",
       "      <td>-0.124088</td>\n",
       "      <td>24.333333</td>\n",
       "      <td>-0.383562</td>\n",
       "      <td>0.023333</td>\n",
       "      <td>-0.571429</td>\n",
       "      <td>52.214286</td>\n",
       "      <td>-0.233926</td>\n",
       "      <td>27.785714</td>\n",
       "      <td>-0.460154</td>\n",
       "      <td>0.102857</td>\n",
       "      <td>-0.902778</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-15</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>43.0</td>\n",
       "      <td>23.0</td>\n",
       "      <td>46.0</td>\n",
       "      <td>39.666667</td>\n",
       "      <td>0.084034</td>\n",
       "      <td>21.000000</td>\n",
       "      <td>0.095238</td>\n",
       "      <td>0.023333</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>51.785714</td>\n",
       "      <td>-0.169655</td>\n",
       "      <td>27.642857</td>\n",
       "      <td>-0.167959</td>\n",
       "      <td>0.102857</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-16</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>46.0</td>\n",
       "      <td>26.0</td>\n",
       "      <td>54.0</td>\n",
       "      <td>43.000000</td>\n",
       "      <td>0.069767</td>\n",
       "      <td>21.333333</td>\n",
       "      <td>0.218750</td>\n",
       "      <td>0.003333</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>51.071429</td>\n",
       "      <td>-0.099301</td>\n",
       "      <td>28.000000</td>\n",
       "      <td>-0.071429</td>\n",
       "      <td>0.102857</td>\n",
       "      <td>-1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1990-03-17</th>\n",
       "      <td>USW00024015</td>\n",
       "      <td>AKRON WASHINGTON CO AIRPORT, CO US</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>54.0</td>\n",
       "      <td>25.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>47.666667</td>\n",
       "      <td>0.132867</td>\n",
       "      <td>24.666667</td>\n",
       "      <td>0.013514</td>\n",
       "      <td>0.006667</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>50.500000</td>\n",
       "      <td>0.069307</td>\n",
       "      <td>27.857143</td>\n",
       "      <td>-0.102564</td>\n",
       "      <td>0.104286</td>\n",
       "      <td>-0.808219</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                station                                               name  \\\n",
       "DATE                                                                         \n",
       "1990-03-07  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-08  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-09  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-10  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-11  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-12  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-13  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-14  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-15  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-16  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-17  USW00093721  BALTIMORE WASHINGTON INTERNATIONAL AIRPORT, MD US   \n",
       "1990-03-07  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-08  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-09  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-10  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-11  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-12  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-13  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-14  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-15  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-16  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "1990-03-17  USW00024015                 AKRON WASHINGTON CO AIRPORT, CO US   \n",
       "\n",
       "            prcp  snow  snwd  tmax  tmin  target  rolling_3_tmax  \\\n",
       "DATE                                                               \n",
       "1990-03-07  0.00   0.0   0.0  37.0  20.0    48.0       44.666667   \n",
       "1990-03-08  0.00   0.0   0.0  48.0  18.0    62.0       43.666667   \n",
       "1990-03-09  0.02   0.0   0.0  62.0  32.0    61.0       49.000000   \n",
       "1990-03-10  0.00   0.0   0.0  61.0  43.0    68.0       57.000000   \n",
       "1990-03-11  0.00   0.0   0.0  68.0  43.0    86.0       63.666667   \n",
       "1990-03-12  0.00   0.0   0.0  86.0  49.0    85.0       71.666667   \n",
       "1990-03-13  0.00   0.0   0.0  85.0  60.0    81.0       79.666667   \n",
       "1990-03-14  0.00   0.0   0.0  81.0  54.0    82.0       84.000000   \n",
       "1990-03-15  0.00   0.0   0.0  82.0  52.0    79.0       82.666667   \n",
       "1990-03-16  0.00   0.0   0.0  79.0  57.0    72.0       80.666667   \n",
       "1990-03-17  0.88   0.0   0.0  72.0  52.0    61.0       77.666667   \n",
       "1990-03-07  0.08   0.0   0.0  32.0  23.0    61.0       43.000000   \n",
       "1990-03-08  0.02   0.3   0.0  61.0  27.0    56.0       45.333333   \n",
       "1990-03-09  0.00   0.0   0.0  56.0  37.0    63.0       49.666667   \n",
       "1990-03-10  0.00   0.0   0.0  63.0  33.0    59.0       60.000000   \n",
       "1990-03-11  0.00   0.0   0.0  59.0  36.0    61.0       59.333333   \n",
       "1990-03-12  0.00   0.0   0.0  61.0  33.0    36.0       61.000000   \n",
       "1990-03-13  0.06   0.3   0.0  36.0  25.0    40.0       52.000000   \n",
       "1990-03-14  0.01   0.1   0.0  40.0  15.0    43.0       45.666667   \n",
       "1990-03-15  0.00   0.0   0.0  43.0  23.0    46.0       39.666667   \n",
       "1990-03-16  0.00   0.0   0.0  46.0  26.0    54.0       43.000000   \n",
       "1990-03-17  0.02   0.0   0.0  54.0  25.0    44.0       47.666667   \n",
       "\n",
       "            rolling_3_tmax_pct  rolling_3_tmin  rolling_3_tmin_pct  \\\n",
       "DATE                                                                 \n",
       "1990-03-07           -0.171642       24.666667           -0.189189   \n",
       "1990-03-08            0.099237       22.666667           -0.205882   \n",
       "1990-03-09            0.265306       23.333333            0.371429   \n",
       "1990-03-10            0.070175       31.000000            0.387097   \n",
       "1990-03-11            0.068063       39.333333            0.093220   \n",
       "1990-03-12            0.200000       45.000000            0.088889   \n",
       "1990-03-13            0.066946       50.666667            0.184211   \n",
       "1990-03-14           -0.035714       54.333333           -0.006135   \n",
       "1990-03-15           -0.008065       55.333333           -0.060241   \n",
       "1990-03-16           -0.020661       54.333333            0.049080   \n",
       "1990-03-17           -0.072961       53.666667           -0.031056   \n",
       "1990-03-07           -0.255814       26.000000           -0.115385   \n",
       "1990-03-08            0.345588       27.333333           -0.012195   \n",
       "1990-03-09            0.127517       29.000000            0.275862   \n",
       "1990-03-10            0.050000       32.333333            0.020619   \n",
       "1990-03-11           -0.005618       35.333333            0.018868   \n",
       "1990-03-12            0.000000       34.000000           -0.029412   \n",
       "1990-03-13           -0.307692       31.333333           -0.202128   \n",
       "1990-03-14           -0.124088       24.333333           -0.383562   \n",
       "1990-03-15            0.084034       21.000000            0.095238   \n",
       "1990-03-16            0.069767       21.333333            0.218750   \n",
       "1990-03-17            0.132867       24.666667            0.013514   \n",
       "\n",
       "            rolling_3_prcp  rolling_3_prcp_pct  rolling_14_tmax  \\\n",
       "DATE                                                              \n",
       "1990-03-07        0.000000            0.000000        46.357143   \n",
       "1990-03-08        0.000000            0.000000        45.500000   \n",
       "1990-03-09        0.006667            2.000000        45.285714   \n",
       "1990-03-10        0.006667           -1.000000        46.357143   \n",
       "1990-03-11        0.006667           -1.000000        49.000000   \n",
       "1990-03-12        0.000000            0.000000        53.000000   \n",
       "1990-03-13        0.000000            0.000000        56.357143   \n",
       "1990-03-14        0.000000            0.000000        58.714286   \n",
       "1990-03-15        0.000000            0.000000        61.285714   \n",
       "1990-03-16        0.000000            0.000000        63.357143   \n",
       "1990-03-17        0.293333            2.000000        64.214286   \n",
       "1990-03-07        0.450000           -0.822222        49.214286   \n",
       "1990-03-08        0.456667           -0.956204        50.571429   \n",
       "1990-03-09        0.033333           -1.000000        50.785714   \n",
       "1990-03-10        0.006667           -1.000000        51.285714   \n",
       "1990-03-11        0.000000            0.000000        51.928571   \n",
       "1990-03-12        0.000000            0.000000        52.285714   \n",
       "1990-03-13        0.020000            2.000000        52.000000   \n",
       "1990-03-14        0.023333           -0.571429        52.214286   \n",
       "1990-03-15        0.023333           -1.000000        51.785714   \n",
       "1990-03-16        0.003333           -1.000000        51.071429   \n",
       "1990-03-17        0.006667            2.000000        50.500000   \n",
       "\n",
       "            rolling_14_tmax_pct  rolling_14_tmin  rolling_14_tmin_pct  \\\n",
       "DATE                                                                    \n",
       "1990-03-07            -0.201849        26.928571            -0.257294   \n",
       "1990-03-08             0.054945        25.928571            -0.305785   \n",
       "1990-03-09             0.369085        25.000000             0.280000   \n",
       "1990-03-10             0.315871        25.857143             0.662983   \n",
       "1990-03-11             0.387755        27.714286             0.551546   \n",
       "1990-03-12             0.622642        30.285714             0.617925   \n",
       "1990-03-13             0.508238        33.357143             0.798715   \n",
       "1990-03-14             0.379562        35.428571             0.524194   \n",
       "1990-03-15             0.337995        37.142857             0.400000   \n",
       "1990-03-16             0.246900        39.214286             0.453552   \n",
       "1990-03-17             0.121246        40.071429             0.297683   \n",
       "1990-03-07            -0.349782        26.214286            -0.122616   \n",
       "1990-03-08             0.206215        26.357143             0.024390   \n",
       "1990-03-09             0.102672        26.857143             0.377660   \n",
       "1990-03-10             0.228412        27.285714             0.209424   \n",
       "1990-03-11             0.136176        28.142857             0.279188   \n",
       "1990-03-12             0.166667        28.500000             0.157895   \n",
       "1990-03-13            -0.307692        28.428571            -0.120603   \n",
       "1990-03-14            -0.233926        27.785714            -0.460154   \n",
       "1990-03-15            -0.169655        27.642857            -0.167959   \n",
       "1990-03-16            -0.099301        28.000000            -0.071429   \n",
       "1990-03-17             0.069307        27.857143            -0.102564   \n",
       "\n",
       "            rolling_14_prcp  rolling_14_prcp_pct  \n",
       "DATE                                              \n",
       "1990-03-07         0.050714            -1.000000  \n",
       "1990-03-08         0.031429            -1.000000  \n",
       "1990-03-09         0.012857             0.555556  \n",
       "1990-03-10         0.005714            -1.000000  \n",
       "1990-03-11         0.005714            -1.000000  \n",
       "1990-03-12         0.005714            -1.000000  \n",
       "1990-03-13         0.005714            -1.000000  \n",
       "1990-03-14         0.005714            -1.000000  \n",
       "1990-03-15         0.005714            -1.000000  \n",
       "1990-03-16         0.002143            -1.000000  \n",
       "1990-03-17         0.064286            12.688889  \n",
       "1990-03-07         0.096429            -0.170370  \n",
       "1990-03-08         0.097857            -0.795620  \n",
       "1990-03-09         0.097857            -1.000000  \n",
       "1990-03-10         0.097857            -1.000000  \n",
       "1990-03-11         0.097857            -1.000000  \n",
       "1990-03-12         0.097857            -1.000000  \n",
       "1990-03-13         0.102143            -0.412587  \n",
       "1990-03-14         0.102857            -0.902778  \n",
       "1990-03-15         0.102857            -1.000000  \n",
       "1990-03-16         0.102857            -1.000000  \n",
       "1990-03-17         0.104286            -0.808219  "
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "weather.loc[\"1990-03-07\": \"1990-03-17\"]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "id": "3887add9",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot: >"
      ]
     },
     "execution_count": 35,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGdCAYAAAAxCSikAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA7VUlEQVR4nO3deXxU9b3/8fcsmQkJJEBCEgKBIEZ2iLLEgIJL2qi4YG9byrWVoq21FxUb61X8qbS3S2x78eLCFbWtelstlLZSpRQbEXAhCAQUkV22CGRjSUL2zJzfH8kMCQaSiZOcWV7Px2Me4Mx3Zj7n3JS873e1GIZhCAAAIIBZzS4AAACgPQQWAAAQ8AgsAAAg4BFYAABAwCOwAACAgEdgAQAAAY/AAgAAAh6BBQAABDy72QX4g9vt1rFjx9SrVy9ZLBazywEAAB1gGIYqKyuVnJwsq/XCfSghEViOHTumlJQUs8sAAACdUFhYqIEDB16wTUgEll69eklquuCYmBiTqwEAAB1RUVGhlJQU7+/xCwmJwOIZBoqJiSGwAAAQZDoynYNJtwAAIOARWAAAQMAjsAAAgIBHYAEAAAGPwAIAAAIegQUAAAQ8AgsAAAh4BBYAABDwCCwAACDgEVgAAEDAI7AAAICAR2ABAAABj8DSSaer67Vk/Wf6/FS12aUAABDyCCyd9OcthXrin7s18/mNKq6oNbscAABCGoGlk05U1UuSjp6u0Xdf2qyK2gaTKwIAIHQRWDqptt7l/fuu4xW6+w8Fqmt0XeAdAACgswgsnVTdHFimj+mvaIdNGz47oR8v3y632zC5MgAAQg+BpZNqGpoCy4TUPlrynfGyWy168+Nj+uWqXW22NwxDBYdP6X/X7VfhSSbqAgDgC7vZBQSrmuYeliiHTVem9dNvvjFWP1r2sX77/kElxkTq+1MvkiTtLqrQGx8d0xsfH9Pnp2okSfmfndAf7swwrXYAAIINgaWTPD0skRE2SdKtlw5USUWdcv+5W79YtUuHTlRpy6FT2lNc6X1PlMOm6nqXPthfprIzdYrv6TSldgAAgg1DQp1U7e1hOZv57pp6keZMSZUkvfrhEe0prpTDZtVXRybq2X+/VAWPfkVjB8bKbUj//OS4GWUDABCU6GHppNrmHpYezT0skmSxWPTY9JGyWiz6rPSMrh+dpOtG9VdsVIS3zU1jk7X983K9uf24vpOZ2t1lAwAQlAgsneTpYenhsLV63mq16LEbR573fdPH9tcvVu3S5kMndby8Rv1je3RpnQAAhAKGhDqppo0elo5I7t1DEwb3kWFI/9jOsBAAAB1BYOmklquEfHXTuGRJ0koCCwAAHUJg6QTDMM72sHQisFw/JklWi/RR4Wn2ZAEAoAMILJ1Q73LL1byjbaSPQ0KSlNArUplD4yRJb24/5tfaAAAIRQSWTqitd3v/3pkhIUm6cWzTsNCbHzMsBABAewgsnVDd0ChJslstirB17hZeNypJdqtFu45XaH/JGX+WBwBAyCGwdELNeZY0+6JPtENXpsVLklYyLAQAwAURWDrBuwdLJ+avtORZLfTmx8dkGJzyDADA+RBYOsGzy21n5694fGVkohx2qz4rrdLuosr23wAAQJgisHSCp4elMyuEWuoVGaGrh/WT1NTLAgAA2kZg6YQaP/WwSC2GhbYzLAQAwPkQWDrBH5NuPa4ZnqAoh02FJ2v08eflX/rzAAAIRQSWTjh7jtCXPzsyymHXtSMSJUkrGRYCAKBNBJZOON9JzZ1109j+kprOFnK7GRYCAOBcBJZO8K4S+pKTbj2mDeunXpF2FVXU6t19pX75TAAAQgmBpROq65t2uvVXD4vTbtM3xqdIkl5874BfPhMAgFBCYOmEmuazhPwVWCTpjitSZbNa9MH+E9pxlMm3AAC0RGDphJrms4S+7E63LQ3sE+Wdy/L8u/SyAADQEoGlEzzLmv2xD0tLd00dKkla9clxFZ6s9utnAwAQzAgsneCvnW7PNTI5RlemxcvlNvS79w/69bMBAAhmBJZO8OdOt+e6a+pFkqRlmwt1qqre758PAEAw6lRgWbx4sVJTUxUZGamMjAxt2rTpgu2XL1+u4cOHKzIyUmPGjNGqVatavX7mzBndc889GjhwoHr06KGRI0dqyZIlnSmtW9T46bTmtlxxcbxG9o9RTYNLf9x42O+fDwBAMPI5sCxbtkw5OTlasGCBtm7dqnHjxik7O1slJSVttt+wYYNmzZqlO++8U9u2bdOMGTM0Y8YM7dixw9smJydHq1ev1h//+Eft2rVL999/v+655x698cYbnb+yLuTd6bYLelgsFot+MK2pl+WV/EPePV8AAAhnPgeWJ598Ut///vc1Z84cb09IVFSUfv/737fZ/qmnntJ1112nBx98UCNGjNDPfvYzXXbZZXr22We9bTZs2KDZs2frqquuUmpqqu666y6NGzeu3Z4bs3RlD4sk3TCmvwb07qGyM/X629ajXfIdAAAEE58CS319vQoKCpSVlXX2A6xWZWVlKT8/v8335Ofnt2ovSdnZ2a3aT548WW+88YaOHj0qwzC0du1a7d27V1/96lfb/My6ujpVVFS0enSns3NYvvxZQm2JsFl1xxVDJEm/fe+AXGzXDwAIcz4FlrKyMrlcLiUmJrZ6PjExUUVFRW2+p6ioqN32zzzzjEaOHKmBAwfK4XDouuuu0+LFizV16tQ2PzM3N1exsbHeR0pKii+X8aWdPUuo6+Ysf2tiimIi7TpQVqW8ncVd9j0AAASDgFgl9Mwzz2jjxo164403VFBQoIULF2ru3Ll6++2322w/f/58lZeXex+FhYXdWu/ZOSxd08MiSdFOu759+WBJ0gvvftZl3wMAQDDw6TdufHy8bDabiotb/3/8xcXFSkpKavM9SUlJF2xfU1OjRx55RK+//rqmT58uSRo7dqw++ugj/fd///cXhpMkyel0yul0+lK637jchuobm7fm76I5LB7fnZKq3753UFuPnNaWQyc1IbVvl34fAACByqceFofDofHjx2vNmjXe59xut9asWaPMzMw235OZmdmqvSTl5eV52zc0NKihoUFWa+tSbDab3G63L+V1i5oWq3a6Yh+WlhJ6Reprlw2QJC1Zz3b9AIDw5fOYRk5OjmbPnq0JEyZo0qRJWrRokaqqqjRnzhxJ0u23364BAwYoNzdXkjRv3jxNmzZNCxcu1PTp07V06VJt2bJFL7zwgiQpJiZG06ZN04MPPqgePXpo8ODBWr9+vf7v//5PTz75pB8v1T88JzVbLJLT3vUjat+fepGWbSnU27uKtb/kjC5O6Nnl3wkAQKDxObDMnDlTpaWlevzxx1VUVKT09HStXr3aO7H2yJEjrXpLJk+erNdee02PPvqoHnnkEaWlpWnFihUaPXq0t83SpUs1f/583XbbbTp58qQGDx6sX/ziF7r77rv9cIn+VVt/djjIYrF0+fcN7ddTXxmRqH/tLNaL7x7Qr74+tsu/EwCAQGMxDCPo18xWVFQoNjZW5eXliomJ6dLv2l1UoesWvae4aIcKHvtKl36XR8HhU/q35zbIYbPqvYeuVmJMZLd8LwAAXcmX398BsUoomHg3jevi+SstjR/cRxNT+6je5dZLHxzqtu8FACBQEFh85F3S3MUrhM71g6lDJUmvbjysytqGbv1uAADMRmDxkaeHpatXCJ3rmuEJujihpyrrGvWnTUe69bsBADAbgcVHnh6WyG7uYbFaLbpratOhiL97/6DqGjkUEQAQPggsPqo2qYdFkmakD1BijFPFFXX6+0fHuv37AQAwC4HFR7UN3T/p1sNht+rO5kMRX3j3gNwciggACBMEFh95Dz6M6LpzhC5k1qRB6uW0a3/JGb2zu8SUGgAA6G4EFh/VdMNJzRfSKzJCtzUfivg8hyICAMIEgcVHnkm3UV14UnN75kxJlcNm1eZDp1Rw+KRpdQAA0F0ILD7y9LB09yqhlhJjInXrpU2HIj7PoYgAgDBAYPGRmauEWvp+8xLnvF3FOnyiytRaAADoagQWH9WatNPtuS5O6KmJqX1kGE1nDQEAEMoILD6qrm+UZM6y5nMNS+olSdpfcsbkSgAA6FoEFh+ZdZZQWy7u11MSgQUAEPoILD4y6yyhtlyc0NzDUkpgAQCENgKLjwKqhyWhqYfl8Ilq1Te6Ta4GAICuQ2DxkWeVUGQA9LAkxjjV02mXy23oECuFAAAhjMDio9qGwBkSslgsGprAPBYAQOgjsPjo7FlC5gcWSUojsAAAwgCBxQeGYZydwxIAPSzS2XksBBYAQCgjsPigrtEtw2j6e6D0sLC0GQAQDggsPvAsaZYCKLA097B8VnpGLrdhcjUAAHQNAosPqpuHgxw2q+y2wLh1KX2j5LBbVdfo1tFTNWaXAwBAlwiM37pBwtPDEijzVyTJZrXoovhoSdL+0kqTqwEAoGsQWHxQE2ArhDyYeAsACHUEFh/UBNAeLC0RWAAAoY7A4gPPSc2R9LAAANCtCCw+CKRdblvyBJZ9JWdkGKwUAgCEHgKLD6oDcNKtJA2Jj5bVIlXWNqq0ss7scgAA8DsCiw8C6aTmlpx2mwb1jZLEsBAAIDQRWHwQiMuaPS5O6CVJ2l9KYAEAhB4Ciw88gSXQ5rBITLwFAIQ2AosPPDvdBtoqIanFxNtiAgsAIPQQWHwQFD0sDAkBAEIQgcUHgbrTrSQN7de0PX9pZZ3KaxpMrgYAAP8isPjAu0rIYTe5ki/qFRmhpJhIScxjAQCEHgKLD6oDuIdFktISm4aFPiOwAABCDIHFB4G6063H0H6eHW85tRkAEFoILD4I1LOEPFjaDAAIVQQWH9Q0uCUFbg8LK4UAAKGKwOKDmuYelkDc6VY6G1g+P1XjHb4CACAUEFh8EKhnCXnERTvUOypChiF9Ri8LACCEEFh8EKinNXtYLBalMY8FABCCCCw+CPRVQhITbwEAoYnA0kENLrcaXIakwB0Sks4ubSawAABCCYGlg2paTGIN1CEhiR4WAEBoIrB0kOccIatFctgC97Z5AsuhE1VqdLlNrgYAAP8I3N+8AebsSc12WSwWk6s5v+TYHuoRYVODy9Dhk9VmlwMAgF8QWDrIMyQUqLvcelitFm8vy75ihoUAAKGBwNJB1fWBv0LIwxNY2IsFABAqCCwdVBvgm8a1xMRbAECoIbB0UKBvGteSJ7DsPFZhciUAAPgHgaWDAn1b/pYmDO4ji0XaU1yp4+U1ZpcDAMCXRmDpIM/Bh8EwhyWup1PpKb0lSWt3l5pbDAAAfkBg6SDPsubIIAgsknTNsARJ0ju7S0yuBACAL4/A0kHVnnOEgmBISJKuHt4UWD7YX+adMAwAQLAisHRQbRBNupWkUckxSoxxqqbBpY0HTphdDgAAXwqBpYOCaZWQJFksFl3T3MuylmEhAECQI7B0UDCtEvK42jOPZU+JDMMwuRoAADqPwNJBNUG0063HlIvj5bBbVXiyhk3kAABBjcDSQcHYwxLttOvyi+IksVoIABDcCCwd5JnDEuiHH57rmmH9JBFYAADBjcDSQZ4eliiH3eRKfHPN8ERJ0pbDp1Re02ByNQAAdA6BpYNqvKuEguuWDYqL0sUJPeVyG3p3L7veAgCCU6d++y5evFipqamKjIxURkaGNm3adMH2y5cv1/DhwxUZGakxY8Zo1apVX2iza9cu3XzzzYqNjVV0dLQmTpyoI0eOdKa8LnF2Dktw9bBIYnkzACDo+RxYli1bppycHC1YsEBbt27VuHHjlJ2drZKStn8ZbtiwQbNmzdKdd96pbdu2acaMGZoxY4Z27NjhbfPZZ5/piiuu0PDhw7Vu3Tpt375djz32mCIjIzt/ZX5WE2T7sLTkWd68bm+pXG6WNwMAgo/F8HGDjoyMDE2cOFHPPvusJMntdislJUX33nuvHn744S+0nzlzpqqqqrRy5Urvc5dffrnS09O1ZMkSSdK3vvUtRURE6A9/+EOnLqKiokKxsbEqLy9XTExMpz6jPZf9LE8nq+r1rx9N1SWJvbrkO7pKg8uty36Wp8raRv3tPybrskF9zC4JAACffn/71MNSX1+vgoICZWVlnf0Aq1VZWVnKz89v8z35+fmt2ktSdna2t73b7dY//vEPXXLJJcrOzlZCQoIyMjK0YsUKX0rrctXNpzUH07JmjwibVVMvaV4ttIthIQBA8PEpsJSVlcnlcikxMbHV84mJiSoqKmrzPUVFRRdsX1JSojNnzuiJJ57Qddddp3/961+69dZb9bWvfU3r169v8zPr6upUUVHR6tGV3G5DtQ1uScE5JCRxejMAILiZPoPU7W4KArfccot+9KMfSZLS09O1YcMGLVmyRNOmTfvCe3Jzc/XTn/6022qsbTx72nEw9rBI0lXD+slikXYer1BRea2SYgNnfhAAAO3xqYclPj5eNptNxcXFrZ4vLi5WUlJSm+9JSkq6YPv4+HjZ7XaNHDmyVZsRI0acd5XQ/PnzVV5e7n0UFhb6chk+80y4lYI3sMT1dCo9pbckae0eelkAAMHFp8DicDg0fvx4rVmzxvuc2+3WmjVrlJmZ2eZ7MjMzW7WXpLy8PG97h8OhiRMnas+ePa3a7N27V4MHD27zM51Op2JiYlo9upJnl1un3Sqr1dKl39WVPMNCa5jHAgAIMj4PCeXk5Gj27NmaMGGCJk2apEWLFqmqqkpz5syRJN1+++0aMGCAcnNzJUnz5s3TtGnTtHDhQk2fPl1Lly7Vli1b9MILL3g/88EHH9TMmTM1depUXX311Vq9erXefPNNrVu3zj9X+SXVNgTfwYdtuXp4ghbm7dUH+8tU2+AKumMGAADhy+fAMnPmTJWWlurxxx9XUVGR0tPTtXr1au/E2iNHjshqPdtxM3nyZL322mt69NFH9cgjjygtLU0rVqzQ6NGjvW1uvfVWLVmyRLm5ubrvvvs0bNgw/fWvf9UVV1zhh0v88jw9LME6HOQxKjlGiTFOFVfUadPBk96VQwAABDqf92EJRF29D8vGAyf0rRc2ami/aK154Cq/f353euDPH+uvWz/X3KuH6sHs4WaXAwAIY122D0u4CuZdbs81aUjTpnGbD50yuRIAADqOwNIB3pOag/AcoXNNSO0rSfqo8LTqWizXBgAgkBFYOsAzhyUyBHpYLoqPVly0Q/WNbu04Wm52OQAAdAiBpQPO9rAEf2CxWCyakNo0LLTpIMNCAIDgQGDpgBrPOUIh0MMiSRObh4W2HDppciUAAHQMgaUDauqD+xyhc3kDy+FTcruDfpEYACAMEFg6oLoheE9qbsvI5Bj1iLCpvKZB+0rOmF0OAADtIrB0QG19aOx06xFhs+qywb0lSZsZFgIABAECSwd4VwmFSA+LJE0Y3DQsRGABAAQDAksH1ITIWUItTRrimXjLSiEAQOAjsHRATYicJdRSekpv2awWHT1do6Ona8wuBwCACyKwdICnhyVUVglJUrTTrtHJTec2sLwZABDoCCwdECqnNZ/Ls03/poMEFgBAYCOwdECtdw5L8J8l1NLZDeSYxwIACGwElg7w9rA4Qut2ebbo31NcqfLqBpOrAQDg/ELrN3AX8c5hCYHTmluK7+nURf2iJUlbDjMsBAAIXASWDvCuEgqhSbceE737sTAsBAAIXASWdhiGEZL7sHh4hoXYQA4AEMgILO2od7nlaj4gMJR2uvXwbCC3/fPT3snFAAAEGgJLO2qbT2qWQrOHZVDfKPXr5VSDy9DHhafNLgcAgDYRWNrhGQ6yWy2KsIXe7bJYLJrkWd58mHksAIDAFHq/gf2sur5RUmhOuPXwzGNhAzkAQKAisLTj7JLm0A0sng3kth4+5Z2vAwBAICGwtMOzpDkU5694jOgfo55OuyrrGrW7qMLscgAA+AICSzs8PSyhuELIw2a16LLBTcNCbNMPAAhEBJZ2VIdBD4skTRzMfiwAgMBFYGmHZ2+SUJ50K0kTm/dj2XjgpNzMYwEABBgCSzu8Bx+G2DlC57p0UG9FOWwqO1OnnceZxwIACCwElnaE8jlCLTntNk0eGi9JWr+31ORqAABojcDSDu85QiE86dbjqmH9JEnr9pSYXAkAAK0RWNoRLj0s0tnAUnD4lMqrG0yuBgCAswgs7fDMYQnlZc0eA/tE6eKEnnIb0nv7GRYCAAQOAks7vENCYdDDIklXXeIZFiKwAAACB4GlHTWes4TCoIdFkq4aliCpaeIty5sBAIGCwNKOmjDZh8Vj4pA+inLYVFrJ8mYAQOAgsLTj7D4s4RFYWN4MAAhEBJZ21IbZHBaJ5c0AgMBDYGmHd5VQGAaWrUdOq7yG5c0AAPMRWNoRThvHeXiWN7vcht7fV2Z2OQAAEFjaE04bx7XkWd68lmEhAEAAILC0I9z2YfFgeTMAIJAQWNoRTjvdtsTyZgBAICGwXIDLbai+0S1JinLYTa6mezUtb46TxPJmAID5CCwX4FnSLIXPPiwtTWseFmJ5MwDAbOHVbeCjaKddB355g+oa3YqMCL9s55l461neHNsjwuSKAADhKvx+C/vIarWoh8Mmi8VidindLqUvy5sBAIGBwIILOnt6M8NCAADzEFhwQZ7lzev2lsowWN4MADAHgQUX1HJ586fHWN4MADAHgQUX1PL05rydxSZXAwAIVwQWtOuGMUmSpJXbjzEsBAAwBYEF7frKyEQ57FZ9VlqlXccrzS4HABCGCCxoV6/ICF09rGm10Jvbj5lcDQAgHBFY0CE3jUuWxLAQAMAcBBZ0yDXDExTlsKnwZI0+/rzc7HIAAGGGwIIOiXLYlTUiUZL05scMCwEAuheBBR3mGRb6x/bjcrsZFgIAdB8CCzps6iXx6hVpV1FFrbYcPmV2OQCAMEJgQYc57TZlj2rak4VhIQBAdyKwwCeeYaFVnxxXo8ttcjUAgHBBYIFPJg+NU99oh05U1Sv/wAmzywEAhAkCC3wSYbPqutHNW/V/fNzkagAA4YLAAp/dNLZpWOifO46rvpFhIQBA1yOwwGeThvRVQi+nKmob9d6+UrPLAQCEAQILfGazWjR9bH9JrBYCAHSPTgWWxYsXKzU1VZGRkcrIyNCmTZsu2H758uUaPny4IiMjNWbMGK1ateq8be+++25ZLBYtWrSoM6Whm9zYPCyUt7NYtQ0uk6sBAIQ6nwPLsmXLlJOTowULFmjr1q0aN26csrOzVVJS0mb7DRs2aNasWbrzzju1bds2zZgxQzNmzNCOHTu+0Pb111/Xxo0blZyc7PuVoFtdNqi3BvTuoap6l9bubvv/9gAA+IvPgeXJJ5/U97//fc2ZM0cjR47UkiVLFBUVpd///vdttn/qqad03XXX6cEHH9SIESP0s5/9TJdddpmeffbZVu2OHj2qe++9V6+++qoiIiI6dzXoNhaLRTeOax4W2s6wEACga/kUWOrr61VQUKCsrKyzH2C1KisrS/n5+W2+Jz8/v1V7ScrOzm7V3u126zvf+Y4efPBBjRo1qt066urqVFFR0eqB7udZLfTO7hKGhQAAXcqnwFJWViaXy6XExMRWzycmJqqoqKjN9xQVFbXb/le/+pXsdrvuu+++DtWRm5ur2NhY7yMlJcWXy4CfjEqOUUIvp2ob3CrgbCEAQBcyfZVQQUGBnnrqKb388suyWCwdes/8+fNVXl7ufRQWFnZxlWiLxWLRFWnxkqT39pWZXA0AIJT5FFji4+Nls9lUXFzc6vni4mIlJSW1+Z6kpKQLtn/vvfdUUlKiQYMGyW63y2636/Dhw3rggQeUmpra5mc6nU7FxMS0esAcU9P6SRL7sQAAupRPgcXhcGj8+PFas2aN9zm32601a9YoMzOzzfdkZma2ai9JeXl53vbf+c53tH37dn300UfeR3Jysh588EG99dZbvl4PutmUi5t6WD49VqGyM3UmVwMACFV2X9+Qk5Oj2bNna8KECZo0aZIWLVqkqqoqzZkzR5J0++23a8CAAcrNzZUkzZs3T9OmTdPChQs1ffp0LV26VFu2bNELL7wgSYqLi1NcXFyr74iIiFBSUpKGDRv2Za8PXaxfL6dG9o/RzuMV+mB/mW5JH2B2SQCAEORzYJk5c6ZKS0v1+OOPq6ioSOnp6Vq9erV3Yu2RI0dktZ7tuJk8ebJee+01Pfroo3rkkUeUlpamFStWaPTo0f67CpjqykvitfN4hd7bR2ABAHQNi2EYhtlFfFkVFRWKjY1VeXk581lM8MH+Mt322w+VGOPUxvnXdnjyNAAgvPny+9v0VUIIfuMH95HTblVxRZ32lZwxuxwAQAgisOBLi4ywKeOipnlI7+5ltRAAwP8ILPCLqezHAgDoQgQW+MWVzfuxfHjwhOoa2aYfAOBfBBb4xSWJPc9u03+IbfoBAP5FYIFftNym/12GhQAAfkZggd+wTT8AoKsQWOA3LbfpP8E2/QAAPyKwwG882/RL0vv7GRYCAPgPgQV+deUlLG8GAPgfgQV+1XIeSwic+gAACBAEFvgV2/QDALoCgQV+1XKbfoaFAAD+QmCB353dpp/lzQAA/yCwwO882/RvPMA2/QAA/yCwwO/Yph8A4G8EFvidxWLR1EuaellWfHTU5GoAAKGAwIIuMWvSIEnSim3HVFxRa3I1AIBgR2BBlxg/uI8mpvZRvcutlz44ZHY5AIAgR2BBl/nB1KGSpFc3HlZlbYPJ1QAAghmBBV3mmuEJujihpyrrGvWnTUfMLgcAEMQILOgyVqtFd029SJL0u/cPqr7RbXJFAIBgRWBBl7olPVmJMU4VV9Tp76wYAgB0EoEFXcppt+mOKUMkSc+/e0BuNwciAgB8R2BBl5uVMUi9nHbtLzmjd3aXmF0OACAIEVjQ5WIiI/Tvlzfty/L8u5+ZXA0AIBgRWNAt7pgyRBE2izYfOqWCwyfNLgcAEGQILOgWiTGRuvXSAZKk59cfMLkaAECwIbCg29zVvJFc3q5i7S85Y3I1AIBgQmBBt7k4oae+MjJRhiG9+C69LACAjiOwoFvdPa1pI7nXtx3VmbpGk6sBAAQLAgu61fjBfZUaF6V6l1sf7C8zuxwAQJAgsKDbXTUsQZK0lj1ZAAAdRGBBt7t6eHNg2VMiw2DnWwBA+wgs6HYZQ/oqMsKq4oo67TpeaXY5AIAgQGBBt4uMsGnK0HhJTb0sAAC0h8ACU1zVPCy0jsACAOgAAgtMcfWwfpKkgsOnVF7dYHI1AIBAR2CBKQb2idIliT3lNqT1+0rNLgcAEOAILDDN1c3Lm9exvBkA0A4CC0zj2Y9l3d5Sud0sbwYAnB+BBaaZkNpHPZ12nayq1/aj5WaXAwAIYAQWmCbCZtWVac3LmxkWAgBcAIEFprqa5c0AgA4gsMBUV13StLz548/LVVpZZ3I1AIBARWCBqRJiIjV6QIwkaf1eljcDANpGYIHpPMub2aYfAHA+BBaYzrO8+d29pWp0uU2uBgAQiAgsMF16Sm/1iYpQZW2jth45bXY5AIAARGCB6WxWi6Y1T75lWAgA0BYCCwKCZ3kz+7EAANpCYEFAmJrWTxaLtLuoUsdO15hdDgAgwBBYEBD6RDt0aUpvSQwLAQC+iMCCgHHtiERJ0lufFptcCQAg0BBYEDBuGNNfkvTB/jKdrKo3uRoAQCAhsCBgDImP1qjkGLnchlbvKDK7HABAACGwIKDcODZZkrRy+zGTKwEABBICCwLKjWObhoU2HjjBYYgAAC8CCwJKSt8ojRsYK7chrd5x3OxyAAABgsCCgOMZFnpzO4EFANCEwIKAM715WGjzoZMqrqg1uRoAQCAgsCDgJPfuofGD+8gwpH/QywIAEIEFAWp6854s//iEwAIAILAgQE0f218Wi1Rw+BRnCwEAOhdYFi9erNTUVEVGRiojI0ObNm26YPvly5dr+PDhioyM1JgxY7Rq1Srvaw0NDXrooYc0ZswYRUdHKzk5WbfffruOHWMfjnCWGBOpial9JTEsBADoRGBZtmyZcnJytGDBAm3dulXjxo1Tdna2SkraPrBuw4YNmjVrlu68805t27ZNM2bM0IwZM7Rjxw5JUnV1tbZu3arHHntMW7du1d/+9jft2bNHN99885e7MgS9m5on37KJHADAYhiG4csbMjIyNHHiRD377LOSJLfbrZSUFN177716+OGHv9B+5syZqqqq0sqVK73PXX755UpPT9eSJUva/I7Nmzdr0qRJOnz4sAYNGtRuTRUVFYqNjVV5ebliYmJ8uRwEsNLKOmX88m25Dem9/7xaKX2jzC4JAOBHvvz+9qmHpb6+XgUFBcrKyjr7AVarsrKylJ+f3+Z78vPzW7WXpOzs7PO2l6Ty8nJZLBb17t27zdfr6upUUVHR6oHQ06+XU5dfFCdJWsmwEACENZ8CS1lZmVwulxITE1s9n5iYqKKitg+rKyoq8ql9bW2tHnroIc2aNeu8aSs3N1exsbHeR0pKii+XgSDC2UIAACnAVgk1NDTom9/8pgzD0HPPPXfedvPnz1d5ebn3UVhY2I1VojtdNzpJNqtFnx6r0MGyKrPLAQCYxKfAEh8fL5vNpuLi4lbPFxcXKykpqc33JCUldai9J6wcPnxYeXl5FxzLcjqdiomJafVAaOob7dCUi+MlSf+glwUAwpZPgcXhcGj8+PFas2aN9zm32601a9YoMzOzzfdkZma2ai9JeXl5rdp7wsq+ffv09ttvKy4uzpeyEOJuHONZLcQ8FgAIVz4PCeXk5OjFF1/UK6+8ol27dumHP/yhqqqqNGfOHEnS7bffrvnz53vbz5s3T6tXr9bChQu1e/du/eQnP9GWLVt0zz33SGoKK1//+te1ZcsWvfrqq3K5XCoqKlJRUZHq6+v9dJkIZtmjkuSwWbW7qFJrd7e9fB4AENrsvr5h5syZKi0t1eOPP66ioiKlp6dr9erV3om1R44ckdV6NgdNnjxZr732mh599FE98sgjSktL04oVKzR69GhJ0tGjR/XGG29IktLT01t919q1a3XVVVd18tIQKmKjIvTdKal64d0D+smbnypzaJwiI2xmlwUA6EY+78MSiNiHJfSdqWvUtQvXqbiiTjlfuUT3XZtmdkkAgC+py/ZhAczS02nX/5s+UpK0eO1+FZ6sNrkiAEB3IrAgaNw0tr8yL4pTXaNb/7Vyp9nlAAC6EYEFQcNisei/bhklu9WivJ3FTMAFgDBCYEFQSUvspTuuGCJJ+smbn6q2wWVyRQCA7kBgQdC579o0JcY4dfhEtV5494DZ5QAAugGBBUGHCbgAEH4ILAhKTMAFgPBCYEFQOncC7ju7i9t/EwAgaBFYELRaTsB96K+fqOxMnckVAQC6CoEFQe1HWZcoLaGnSivr9OPlH8vtDvqNmwEAbSCwIKj1cNj07L9fJqfdqnV7SvW79w+aXRIAoAsQWBD0hiX10mM3Nq0a+vVbu7X989PmFgQA8DsCC0LCbRmDdP3oJDW4DN37p22qrG0wuyQAgB8RWBASLBaLnvjaWA3o3UOHT1Tr0RU7FAIHkQMAmhFYEDJioyL09Kx02awW/f2jY/pLwedmlwQA8BMCC0LK+MF9lfOVSyRJj//9U31WesbkigAA/kBgQci5e9pQTR4ap5oGl+59bZvqGjkgEQCCHYEFIcdmteh/Zqarb7RDO49X6Nl39ptdEgDgSyKwICQlxkTqFzNGS5L+d91n+uTzcpMrAgB8GQQWhKzrx/TXjWP7y+U29OPlHzM0BABBjMCCkPbTm0cpLtqhPcWVDA0BQBAjsCCkxfV06ucMDQFA0COwIOSdOzRU3+g2uyQAgI8ILAgLLYeGnnlnn9nlAAB8RGBBWGBoCACCG4EFYeP6Mf01naEhAAhKBBaElf9qMTT00F+3q7q+0eySAAAdQGBBWInr6dQvvzZGFov0+rajuuGp97TtyCmzywIAtIPAgrCTPSpJr96Zof6xkTp0olpfX5KvRW/vVaOLISIACFQEFoSlyRfHa/W8qbp5XLJcbkOL3t6nry/J18GyKrNLAwC0gcCCsBUbFaGnZ12qp76Vrl6Rdn1UeFo3PPWelm0+YnZpAIBzEFgQ9m5JH6DV90/V5Rf1VU2DSw/99RO9s7vY7LIAAC0QWABJA3r30Gvfu1y3ZQySJP1s5S6WPQNAACGwAM2sVosevn644ns6dbCsSv+Xf8jskgAAzQgsQAu9IiP0n9nDJElPvb1PZWfqTK4IACARWIAv+Pr4gRozIFaVdY1a+K89ZpcDABCBBfgCq9Wix28aKUlaurlQnx7j3CEAMBuBBWjDxNS+umlcsgxD+umbO2UYhtklAUBYI7AA5/Hw9cMVGWHVpoMnteqTIrPLAYCwRmABzmNA7x76wdShkqRfrtql2gaXyRUBQPgisAAXcPe0oeofG6mjp2v04rsHzC4HAMIWgQW4gB4Omx6+frgk6X/Xfabj5TUmVwQA4YnAArTj5nHJmjC4j2oaXPrBHwq04yirhgCguxFYgHZYLBb99JZRinbYtP3zct307Pua/7ftOsGmcgDQbQgsQAeMSo7V2w9M0y3pTUud/7SpUFf99zr9/v2DanBx5hAAdDWLEQIbTFRUVCg2Nlbl5eWKiYkxuxyEuM2HTuonb3yqT49VSJLSEnpqwU2jdEVavMmVAUBw8eX3N4EF6ASX29CyzYX673/t0cmqeklS1ogEPXLDCF3Ur6fJ1QFAcCCwAN2kvLpB//P2Xv1h42G53IbsVotuz0zVvGvTFBsVYXZ5ABDQCCxAN9tfcka/XLVL7+wukST1jorQ/dem6bbLByvCxlQxAGgLgQUwyXv7SvXzlbu0p7hSkjS0X7R+cvMoXZnWz+TKACDwEFgAEzW63Fq6uVD/k7dXJ5rnt9yWMUiP3DBC0U67ydUBQODw5fc3fdWAn9ltVn378sFa++BVmp05WJL06odHdMPT72nzoZMmVwcAwYnAAnSRmMgI/fSW0Xr1exlKjo3U4RPV+ubz+crlIEUA8BmBBehiUy6O1+ofTdU3xg+UYUjPv3tANz/7Plv8A4APmMMCdKO8ncWa/7ftKjtTL6tFuvyiON0wpr+uG52k+J5Os8sDgG7FpFsggJ2sqtdjK3boH58c9z5HeAEQjggsQBAoPFmtf3xyXKs+Oa7tn58dHrJapKmX9NN3J6dqalo/Wa0WE6sEgK5DYAGCzPnCy9B+0fru5FR97bKBLIkGEHIILEAQO1hWpf/LP6TlWz7XmbpGSVJMpF0zJ6bo9sxUpfSNMrlCAPAPAgsQAiprG/SXgs/1yoZDOnSiWpJksUjpKb2VNSJR1wxP0PCkXrJYGDICEJwILEAIcbsNrdtbopc+OKT39pW1em1A7x66dkSCrh2RqDEDYhXttMlpt5lUKQD4hsAChKii8lq9s7tEa3YV6/39ZaprdH+hTYTNominXdEOu6KdNvWOcmj6mP76+njmwQAILAQWIAzU1Lv0wf4yrdldrLW7S1VUUXvB9r2cdn1zYoq+O5l5MAACA4EFCEONLreq6l2qrm9UVV2jztS5VFXXqL3FlfpD/mEdKKuS1DQPJmtEouZMSVXmRXHMgQFgmi4PLIsXL9ZvfvMbFRUVady4cXrmmWc0adKk87Zfvny5HnvsMR06dEhpaWn61a9+pRtuuMH7umEYWrBggV588UWdPn1aU6ZM0XPPPae0tLQO1UNgAS7M7Ta0fl+pXvrgkN7dW+p9vk9UhPpEORTTI0KxLR69m5+P6+lQXLSz+U+H+kQ7FGHjRA8A/uHL72+fB7SXLVumnJwcLVmyRBkZGVq0aJGys7O1Z88eJSQkfKH9hg0bNGvWLOXm5urGG2/Ua6+9phkzZmjr1q0aPXq0JOnXv/61nn76ab3yyisaMmSIHnvsMWVnZ2vnzp2KjIz0tUQA57BaLbp6WIKuHpag/SWVennDIf214KhOVTfoVHWDT5/VK9LeKty0fMScE3o8f+8RYVOj21Cjy1Cj2y2X21CDy5DbMJTQy6l+vZz09AC4IJ97WDIyMjRx4kQ9++yzkiS3262UlBTde++9evjhh7/QfubMmaqqqtLKlSu9z11++eVKT0/XkiVLZBiGkpOT9cADD+jHP/6xJKm8vFyJiYl6+eWX9a1vfavdmuhhAXx3pq5RhSerVV7T4H1UNP95urpBJ6vrdeJMnU5W1evEmXqdqq6Xu4sGkCMjrErpE6VBfaOU0rfpz4QYp6znCTFWi0V2q0U2W9OfdqtV9ua/O+02RUZYFRlhk9N+9k87PUNAwOmyHpb6+noVFBRo/vz53uesVquysrKUn5/f5nvy8/OVk5PT6rns7GytWLFCknTw4EEVFRUpKyvL+3psbKwyMjKUn5/focACwHc9nXaN6N/xgO9yGzpdXa9T1a3Djefvp1v897mv1zS4FNEcKmzW5pDRHCBOnKlTbYNb+0rOaF/Jma663OYw0xRgPCHGGWGTw26Vr307hiSX293cY2TI5W7qOWp0GbJIzUHK2hSqrJbmMGX1fmdkiz8jI2yyWS2igwmBzm616P9NH2ne9/vSuKysTC6XS4mJia2eT0xM1O7du9t8T1FRUZvti4qKvK97njtfm3PV1dWprq7O+98VFRW+XAaATrBZLYrr6VScnw9mbHC5dex0jY6crPY+Ck9Wq6yyvs32hgy5DTUPMbmbw0JTaKhvdKuu0a26BpfqGt2qd51d9t3oNtRY71JVvcuv9QPhwmG3Bk9gCRS5ubn66U9/anYZAPwgwmbV4LhoDY6L9vtne0JMbXOAaevP+jb2sukIT6+JZyjK1jw0Jampt6V5zo7LbaihufelrtGluga3alv8WdvglsvduRqA7mSzmjus6lNgiY+Pl81mU3Fxcavni4uLlZSU1OZ7kpKSLtje82dxcbH69+/fqk16enqbnzl//vxWw0wVFRVKSUnx5VIAhAGb1aIeDpt6ONj9Fwh2PsUlh8Oh8ePHa82aNd7n3G631qxZo8zMzDbfk5mZ2aq9JOXl5XnbDxkyRElJSa3aVFRU6MMPPzzvZzqdTsXExLR6AACA0OXzkFBOTo5mz56tCRMmaNKkSVq0aJGqqqo0Z84cSdLtt9+uAQMGKDc3V5I0b948TZs2TQsXLtT06dO1dOlSbdmyRS+88IIkyWKx6P7779fPf/5zpaWleZc1Jycna8aMGf67UgAAELR8DiwzZ85UaWmpHn/8cRUVFSk9PV2rV6/2Tpo9cuSIrC3GuSZPnqzXXntNjz76qB555BGlpaVpxYoV3j1YJOk///M/VVVVpbvuukunT5/WFVdcodWrV7MHCwAAkMTW/AAAwCS+/P5mJyUAABDwCCwAACDgEVgAAEDAI7AAAICAR2ABAAABj8ACAAACHoEFAAAEPAILAAAIeAQWAAAQ8Hzemj8QeTbrraioMLkSAADQUZ7f2x3ZdD8kAktlZaUkKSUlxeRKAACAryorKxUbG3vBNiFxlpDb7daxY8fUq1cvWSwWv352RUWFUlJSVFhYGLbnFHEPmnAfuAcS98CD+8A9kL78PTAMQ5WVlUpOTm51cHJbQqKHxWq1auDAgV36HTExMWH7A+nBPWjCfeAeSNwDD+4D90D6cvegvZ4VDybdAgCAgEdgAQAAAY/A0g6n06kFCxbI6XSaXYppuAdNuA/cA4l74MF94B5I3XsPQmLSLQAACG30sAAAgIBHYAEAAAGPwAIAAAIegQUAAAQ8Aks7Fi9erNTUVEVGRiojI0ObNm0yu6Qu8+677+qmm25ScnKyLBaLVqxY0ep1wzD0+OOPq3///urRo4eysrK0b98+c4rtIrm5uZo4caJ69eqlhIQEzZgxQ3v27GnVpra2VnPnzlVcXJx69uypf/u3f1NxcbFJFfvfc889p7Fjx3o3gsrMzNQ///lP7+uhfv1teeKJJ2SxWHT//fd7nwuH+/CTn/xEFoul1WP48OHe18PhHkjS0aNH9e1vf1txcXHq0aOHxowZoy1btnhfD4d/G1NTU7/ws2CxWDR37lxJ3fOzQGC5gGXLliknJ0cLFizQ1q1bNW7cOGVnZ6ukpMTs0rpEVVWVxo0bp8WLF7f5+q9//Ws9/fTTWrJkiT788ENFR0crOztbtbW13Vxp11m/fr3mzp2rjRs3Ki8vTw0NDfrqV7+qqqoqb5sf/ehHevPNN7V8+XKtX79ex44d09e+9jUTq/avgQMH6oknnlBBQYG2bNmia665Rrfccos+/fRTSaF//efavHmznn/+eY0dO7bV8+FyH0aNGqXjx497H++//773tXC4B6dOndKUKVMUERGhf/7zn9q5c6cWLlyoPn36eNuEw7+NmzdvbvVzkJeXJ0n6xje+IambfhYMnNekSZOMuXPnev/b5XIZycnJRm5urolVdQ9Jxuuvv+79b7fbbSQlJRm/+c1vvM+dPn3acDqdxp/+9CcTKuweJSUlhiRj/fr1hmE0XXNERISxfPlyb5tdu3YZkoz8/Hyzyuxyffr0MX7729+G3fVXVlYaaWlpRl5enjFt2jRj3rx5hmGEz8/BggULjHHjxrX5Wrjcg4ceesi44oorzvt6uP7bOG/ePGPo0KGG2+3utp8FeljOo76+XgUFBcrKyvI+Z7ValZWVpfz8fBMrM8fBgwdVVFTU6n7ExsYqIyMjpO9HeXm5JKlv376SpIKCAjU0NLS6D8OHD9egQYNC8j64XC4tXbpUVVVVyszMDLvrnzt3rqZPn97qeqXw+jnYt2+fkpOTddFFF+m2227TkSNHJIXPPXjjjTc0YcIEfeMb31BCQoIuvfRSvfjii97Xw/Hfxvr6ev3xj3/UHXfcIYvF0m0/CwSW8ygrK5PL5VJiYmKr5xMTE1VUVGRSVebxXHM43Q+32637779fU6ZM0ejRoyU13QeHw6HevXu3ahtq9+GTTz5Rz5495XQ6dffdd+v111/XyJEjw+b6JWnp0qXaunWrcnNzv/BauNyHjIwMvfzyy1q9erWee+45HTx4UFdeeaUqKyvD5h4cOHBAzz33nNLS0vTWW2/phz/8oe677z698sorksLz38YVK1bo9OnT+u53vyup+/73EBKnNQNdYe7cudqxY0erMftwMWzYMH300UcqLy/XX/7yF82ePVvr1683u6xuU1hYqHnz5ikvL0+RkZFml2Oa66+/3vv3sWPHKiMjQ4MHD9af//xn9ejRw8TKuo/b7daECRP0y1/+UpJ06aWXaseOHVqyZIlmz55tcnXm+N3vfqfrr79eycnJ3fq99LCcR3x8vGw22xdmORcXFyspKcmkqszjueZwuR/33HOPVq5cqbVr12rgwIHe55OSklRfX6/Tp0+3ah9q98HhcOjiiy/W+PHjlZubq3Hjxumpp54Km+svKChQSUmJLrvsMtntdtntdq1fv15PP/207Ha7EhMTw+I+nKt379665JJLtH///rD5Wejfv79GjhzZ6rkRI0Z4h8bC7d/Gw4cP6+2339b3vvc973Pd9bNAYDkPh8Oh8ePHa82aNd7n3G631qxZo8zMTBMrM8eQIUOUlJTU6n5UVFToww8/DKn7YRiG7rnnHr3++ut65513NGTIkFavjx8/XhEREa3uw549e3TkyJGQug/ncrvdqqurC5vrv/baa/XJJ5/oo48+8j4mTJig2267zfv3cLgP5zpz5ow+++wz9e/fP2x+FqZMmfKFrQ327t2rwYMHSwqffxs9XnrpJSUkJGj69One57rtZ8Fv03dD0NKlSw2n02m8/PLLxs6dO4277rrL6N27t1FUVGR2aV2isrLS2LZtm7Ft2zZDkvHkk08a27ZtMw4fPmwYhmE88cQTRu/evY2///3vxvbt241bbrnFGDJkiFFTU2Ny5f7zwx/+0IiNjTXWrVtnHD9+3Puorq72trn77ruNQYMGGe+8846xZcsWIzMz08jMzDSxav96+OGHjfXr1xsHDx40tm/fbjz88MOGxWIx/vWvfxmGEfrXfz4tVwkZRnjchwceeMBYt26dcfDgQeODDz4wsrKyjPj4eKOkpMQwjPC4B5s2bTLsdrvxi1/8wti3b5/x6quvGlFRUcYf//hHb5tw+LfRMJpWyg4aNMh46KGHvvBad/wsEFja8cwzzxiDBg0yHA6HMWnSJGPjxo1ml9Rl1q5da0j6wmP27NmGYTQt33vssceMxMREw+l0Gtdee62xZ88ec4v2s7auX5Lx0ksvedvU1NQY//Ef/2H06dPHiIqKMm699Vbj+PHj5hXtZ3fccYcxePBgw+FwGP369TOuvfZab1gxjNC//vM5N7CEw32YOXOm0b9/f8PhcBgDBgwwZs6caezfv9/7ejjcA8MwjDfffNMYPXq04XQ6jeHDhxsvvPBCq9fD4d9GwzCMt956y5DU5rV1x8+CxTAMw3/9NQAAAP7HHBYAABDwCCwAACDgEVgAAEDAI7AAAICAR2ABAAABj8ACAAACHoEFAAAEPAILAAAIeAQWAAAQ8AgsAAAg4BFYAABAwCOwAACAgPf/Ac+vMb2BUPnZAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "(predictions[\"diff\"].round().value_counts().sort_index() / predictions.shape[0]).plot()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "id": "a343aa78",
   "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>1979-12-30</th>\n",
       "      <td>50.0</td>\n",
       "      <td>56.730245</td>\n",
       "      <td>6.730245</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1979-12-31</th>\n",
       "      <td>45.0</td>\n",
       "      <td>50.640905</td>\n",
       "      <td>5.640905</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-01-01</th>\n",
       "      <td>38.0</td>\n",
       "      <td>46.213188</td>\n",
       "      <td>8.213188</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-01-02</th>\n",
       "      <td>44.0</td>\n",
       "      <td>41.945868</td>\n",
       "      <td>2.054132</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1980-01-03</th>\n",
       "      <td>31.0</td>\n",
       "      <td>46.439711</td>\n",
       "      <td>15.439711</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-09</th>\n",
       "      <td>103.0</td>\n",
       "      <td>95.253513</td>\n",
       "      <td>7.746487</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-10</th>\n",
       "      <td>84.0</td>\n",
       "      <td>98.319716</td>\n",
       "      <td>14.319716</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-11</th>\n",
       "      <td>93.0</td>\n",
       "      <td>83.513123</td>\n",
       "      <td>9.486877</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-12</th>\n",
       "      <td>99.0</td>\n",
       "      <td>88.691876</td>\n",
       "      <td>10.308124</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2022-07-13</th>\n",
       "      <td>99.0</td>\n",
       "      <td>92.989566</td>\n",
       "      <td>6.010434</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>34584 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            actual  prediction       diff\n",
       "DATE                                     \n",
       "1979-12-30    50.0   56.730245   6.730245\n",
       "1979-12-31    45.0   50.640905   5.640905\n",
       "1980-01-01    38.0   46.213188   8.213188\n",
       "1980-01-02    44.0   41.945868   2.054132\n",
       "1980-01-03    31.0   46.439711  15.439711\n",
       "...            ...         ...        ...\n",
       "2022-07-09   103.0   95.253513   7.746487\n",
       "2022-07-10    84.0   98.319716  14.319716\n",
       "2022-07-11    93.0   83.513123   9.486877\n",
       "2022-07-12    99.0   88.691876  10.308124\n",
       "2022-07-13    99.0   92.989566   6.010434\n",
       "\n",
       "[34584 rows x 3 columns]"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "predictions\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9413a900",
   "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.10.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}