House Sales_in_King_Count_USA.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " <a href=\"https://www.bigdatauniversity.com\"><img src = \"https://ibm.box.com/shared/static/ugcqz6ohbvff804xp84y4kqnvvk3bq1g.png\" width = 300, align = \"center\"></a>\n",
    "\n",
    "<h1 align=center><font size = 5>Data Analysis with Python</font></h1>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# House Sales in King County, USA"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This dataset contains house sale prices for King County, which includes Seattle. It includes homes sold between May 2014 and May 2015."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<b>id</b> :a notation for a house\n",
    "\n",
    "<b> date</b>: Date house was sold\n",
    "\n",
    "\n",
    "<b>price</b>: Price is prediction target\n",
    "\n",
    "\n",
    "<b>bedrooms</b>: Number of Bedrooms/House\n",
    "\n",
    "\n",
    "<b>bathrooms</b>: Number of bathrooms/bedrooms\n",
    "\n",
    "<b>sqft_living</b>: square footage of the home\n",
    "\n",
    "<b>sqft_lot</b>: square footage of the lot\n",
    "\n",
    "\n",
    "<b>floors</b> :Total floors (levels) in house\n",
    "\n",
    "\n",
    "<b>waterfront</b> :House which has a view to a waterfront\n",
    "\n",
    "\n",
    "<b>view</b>: Has been viewed\n",
    "\n",
    "\n",
    "<b>condition</b> :How good the condition is  Overall\n",
    "\n",
    "<b>grade</b>: overall grade given to the housing unit, based on King County grading system\n",
    "\n",
    "\n",
    "<b>sqft_above</b> :square footage of house apart from basement\n",
    "\n",
    "\n",
    "<b>sqft_basement</b>: square footage of the basement\n",
    "\n",
    "<b>yr_built</b> :Built Year\n",
    "\n",
    "\n",
    "<b>yr_renovated</b> :Year when house was renovated\n",
    "\n",
    "<b>zipcode</b>:zip code\n",
    "\n",
    "\n",
    "<b>lat</b>: Latitude coordinate\n",
    "\n",
    "<b>long</b>: Longitude coordinate\n",
    "\n",
    "<b>sqft_living15</b> :Living room area in 2015(implies-- some renovations) This might or might not have affected the lotsize area\n",
    "\n",
    "\n",
    "<b>sqft_lot15</b> :lotSize area in 2015(implies-- some renovations)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You will require the following libraries "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import seaborn as sns\n",
    "from sklearn.pipeline import Pipeline\n",
    "from sklearn.preprocessing import StandardScaler,PolynomialFeatures\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 1.0 Importing the Data "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " Load the csv:  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "file_name='https://s3-api.us-geo.objectstorage.softlayer.net/cf-courses-data/CognitiveClass/DA0101EN/coursera/project/kc_house_data_NaN.csv'\n",
    "df=pd.read_csv(file_name)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "we use the method <code>head</code> to display the first 5 columns of the dataframe."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "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>Unnamed: 0</th>\n",
       "      <th>id</th>\n",
       "      <th>date</th>\n",
       "      <th>price</th>\n",
       "      <th>bedrooms</th>\n",
       "      <th>bathrooms</th>\n",
       "      <th>sqft_living</th>\n",
       "      <th>sqft_lot</th>\n",
       "      <th>floors</th>\n",
       "      <th>waterfront</th>\n",
       "      <th>...</th>\n",
       "      <th>grade</th>\n",
       "      <th>sqft_above</th>\n",
       "      <th>sqft_basement</th>\n",
       "      <th>yr_built</th>\n",
       "      <th>yr_renovated</th>\n",
       "      <th>zipcode</th>\n",
       "      <th>lat</th>\n",
       "      <th>long</th>\n",
       "      <th>sqft_living15</th>\n",
       "      <th>sqft_lot15</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>7129300520</td>\n",
       "      <td>20141013T000000</td>\n",
       "      <td>221900.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1180</td>\n",
       "      <td>5650</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>7</td>\n",
       "      <td>1180</td>\n",
       "      <td>0</td>\n",
       "      <td>1955</td>\n",
       "      <td>0</td>\n",
       "      <td>98178</td>\n",
       "      <td>47.5112</td>\n",
       "      <td>-122.257</td>\n",
       "      <td>1340</td>\n",
       "      <td>5650</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>6414100192</td>\n",
       "      <td>20141209T000000</td>\n",
       "      <td>538000.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>2.25</td>\n",
       "      <td>2570</td>\n",
       "      <td>7242</td>\n",
       "      <td>2.0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>7</td>\n",
       "      <td>2170</td>\n",
       "      <td>400</td>\n",
       "      <td>1951</td>\n",
       "      <td>1991</td>\n",
       "      <td>98125</td>\n",
       "      <td>47.7210</td>\n",
       "      <td>-122.319</td>\n",
       "      <td>1690</td>\n",
       "      <td>7639</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>5631500400</td>\n",
       "      <td>20150225T000000</td>\n",
       "      <td>180000.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>1.00</td>\n",
       "      <td>770</td>\n",
       "      <td>10000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>6</td>\n",
       "      <td>770</td>\n",
       "      <td>0</td>\n",
       "      <td>1933</td>\n",
       "      <td>0</td>\n",
       "      <td>98028</td>\n",
       "      <td>47.7379</td>\n",
       "      <td>-122.233</td>\n",
       "      <td>2720</td>\n",
       "      <td>8062</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>2487200875</td>\n",
       "      <td>20141209T000000</td>\n",
       "      <td>604000.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>3.00</td>\n",
       "      <td>1960</td>\n",
       "      <td>5000</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>7</td>\n",
       "      <td>1050</td>\n",
       "      <td>910</td>\n",
       "      <td>1965</td>\n",
       "      <td>0</td>\n",
       "      <td>98136</td>\n",
       "      <td>47.5208</td>\n",
       "      <td>-122.393</td>\n",
       "      <td>1360</td>\n",
       "      <td>5000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>1954400510</td>\n",
       "      <td>20150218T000000</td>\n",
       "      <td>510000.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>2.00</td>\n",
       "      <td>1680</td>\n",
       "      <td>8080</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>8</td>\n",
       "      <td>1680</td>\n",
       "      <td>0</td>\n",
       "      <td>1987</td>\n",
       "      <td>0</td>\n",
       "      <td>98074</td>\n",
       "      <td>47.6168</td>\n",
       "      <td>-122.045</td>\n",
       "      <td>1800</td>\n",
       "      <td>7503</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 22 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0          id             date     price  bedrooms  bathrooms  \\\n",
       "0           0  7129300520  20141013T000000  221900.0       3.0       1.00   \n",
       "1           1  6414100192  20141209T000000  538000.0       3.0       2.25   \n",
       "2           2  5631500400  20150225T000000  180000.0       2.0       1.00   \n",
       "3           3  2487200875  20141209T000000  604000.0       4.0       3.00   \n",
       "4           4  1954400510  20150218T000000  510000.0       3.0       2.00   \n",
       "\n",
       "   sqft_living  sqft_lot  floors  waterfront  ...  grade  sqft_above  \\\n",
       "0         1180      5650     1.0           0  ...      7        1180   \n",
       "1         2570      7242     2.0           0  ...      7        2170   \n",
       "2          770     10000     1.0           0  ...      6         770   \n",
       "3         1960      5000     1.0           0  ...      7        1050   \n",
       "4         1680      8080     1.0           0  ...      8        1680   \n",
       "\n",
       "   sqft_basement  yr_built  yr_renovated  zipcode      lat     long  \\\n",
       "0              0      1955             0    98178  47.5112 -122.257   \n",
       "1            400      1951          1991    98125  47.7210 -122.319   \n",
       "2              0      1933             0    98028  47.7379 -122.233   \n",
       "3            910      1965             0    98136  47.5208 -122.393   \n",
       "4              0      1987             0    98074  47.6168 -122.045   \n",
       "\n",
       "   sqft_living15  sqft_lot15  \n",
       "0           1340        5650  \n",
       "1           1690        7639  \n",
       "2           2720        8062  \n",
       "3           1360        5000  \n",
       "4           1800        7503  \n",
       "\n",
       "[5 rows x 22 columns]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Question 1 \n",
    "Display the data types of each column using the attribute dtype, then take a screenshot and submit it, include your code in the image. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Unnamed: 0         int64\n",
       "id                 int64\n",
       "date              object\n",
       "price            float64\n",
       "bedrooms         float64\n",
       "bathrooms        float64\n",
       "sqft_living        int64\n",
       "sqft_lot           int64\n",
       "floors           float64\n",
       "waterfront         int64\n",
       "view               int64\n",
       "condition          int64\n",
       "grade              int64\n",
       "sqft_above         int64\n",
       "sqft_basement      int64\n",
       "yr_built           int64\n",
       "yr_renovated       int64\n",
       "zipcode            int64\n",
       "lat              float64\n",
       "long             float64\n",
       "sqft_living15      int64\n",
       "sqft_lot15         int64\n",
       "dtype: object"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.dtypes"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We use the method describe to obtain a statistical summary of the dataframe."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Unnamed: 0</th>\n",
       "      <th>id</th>\n",
       "      <th>price</th>\n",
       "      <th>bedrooms</th>\n",
       "      <th>bathrooms</th>\n",
       "      <th>sqft_living</th>\n",
       "      <th>sqft_lot</th>\n",
       "      <th>floors</th>\n",
       "      <th>waterfront</th>\n",
       "      <th>view</th>\n",
       "      <th>...</th>\n",
       "      <th>grade</th>\n",
       "      <th>sqft_above</th>\n",
       "      <th>sqft_basement</th>\n",
       "      <th>yr_built</th>\n",
       "      <th>yr_renovated</th>\n",
       "      <th>zipcode</th>\n",
       "      <th>lat</th>\n",
       "      <th>long</th>\n",
       "      <th>sqft_living15</th>\n",
       "      <th>sqft_lot15</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>21613.00000</td>\n",
       "      <td>2.161300e+04</td>\n",
       "      <td>2.161300e+04</td>\n",
       "      <td>21600.000000</td>\n",
       "      <td>21603.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>2.161300e+04</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>10806.00000</td>\n",
       "      <td>4.580302e+09</td>\n",
       "      <td>5.400881e+05</td>\n",
       "      <td>3.372870</td>\n",
       "      <td>2.115736</td>\n",
       "      <td>2079.899736</td>\n",
       "      <td>1.510697e+04</td>\n",
       "      <td>1.494309</td>\n",
       "      <td>0.007542</td>\n",
       "      <td>0.234303</td>\n",
       "      <td>...</td>\n",
       "      <td>7.656873</td>\n",
       "      <td>1788.390691</td>\n",
       "      <td>291.509045</td>\n",
       "      <td>1971.005136</td>\n",
       "      <td>84.402258</td>\n",
       "      <td>98077.939805</td>\n",
       "      <td>47.560053</td>\n",
       "      <td>-122.213896</td>\n",
       "      <td>1986.552492</td>\n",
       "      <td>12768.455652</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>6239.28002</td>\n",
       "      <td>2.876566e+09</td>\n",
       "      <td>3.671272e+05</td>\n",
       "      <td>0.926657</td>\n",
       "      <td>0.768996</td>\n",
       "      <td>918.440897</td>\n",
       "      <td>4.142051e+04</td>\n",
       "      <td>0.539989</td>\n",
       "      <td>0.086517</td>\n",
       "      <td>0.766318</td>\n",
       "      <td>...</td>\n",
       "      <td>1.175459</td>\n",
       "      <td>828.090978</td>\n",
       "      <td>442.575043</td>\n",
       "      <td>29.373411</td>\n",
       "      <td>401.679240</td>\n",
       "      <td>53.505026</td>\n",
       "      <td>0.138564</td>\n",
       "      <td>0.140828</td>\n",
       "      <td>685.391304</td>\n",
       "      <td>27304.179631</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>0.00000</td>\n",
       "      <td>1.000102e+06</td>\n",
       "      <td>7.500000e+04</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>290.000000</td>\n",
       "      <td>5.200000e+02</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>290.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1900.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>98001.000000</td>\n",
       "      <td>47.155900</td>\n",
       "      <td>-122.519000</td>\n",
       "      <td>399.000000</td>\n",
       "      <td>651.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>5403.00000</td>\n",
       "      <td>2.123049e+09</td>\n",
       "      <td>3.219500e+05</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>1.750000</td>\n",
       "      <td>1427.000000</td>\n",
       "      <td>5.040000e+03</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>7.000000</td>\n",
       "      <td>1190.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1951.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>98033.000000</td>\n",
       "      <td>47.471000</td>\n",
       "      <td>-122.328000</td>\n",
       "      <td>1490.000000</td>\n",
       "      <td>5100.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>10806.00000</td>\n",
       "      <td>3.904930e+09</td>\n",
       "      <td>4.500000e+05</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>2.250000</td>\n",
       "      <td>1910.000000</td>\n",
       "      <td>7.618000e+03</td>\n",
       "      <td>1.500000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>7.000000</td>\n",
       "      <td>1560.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1975.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>98065.000000</td>\n",
       "      <td>47.571800</td>\n",
       "      <td>-122.230000</td>\n",
       "      <td>1840.000000</td>\n",
       "      <td>7620.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>16209.00000</td>\n",
       "      <td>7.308900e+09</td>\n",
       "      <td>6.450000e+05</td>\n",
       "      <td>4.000000</td>\n",
       "      <td>2.500000</td>\n",
       "      <td>2550.000000</td>\n",
       "      <td>1.068800e+04</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>8.000000</td>\n",
       "      <td>2210.000000</td>\n",
       "      <td>560.000000</td>\n",
       "      <td>1997.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>98118.000000</td>\n",
       "      <td>47.678000</td>\n",
       "      <td>-122.125000</td>\n",
       "      <td>2360.000000</td>\n",
       "      <td>10083.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>21612.00000</td>\n",
       "      <td>9.900000e+09</td>\n",
       "      <td>7.700000e+06</td>\n",
       "      <td>33.000000</td>\n",
       "      <td>8.000000</td>\n",
       "      <td>13540.000000</td>\n",
       "      <td>1.651359e+06</td>\n",
       "      <td>3.500000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>4.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>13.000000</td>\n",
       "      <td>9410.000000</td>\n",
       "      <td>4820.000000</td>\n",
       "      <td>2015.000000</td>\n",
       "      <td>2015.000000</td>\n",
       "      <td>98199.000000</td>\n",
       "      <td>47.777600</td>\n",
       "      <td>-121.315000</td>\n",
       "      <td>6210.000000</td>\n",
       "      <td>871200.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>8 rows × 21 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        Unnamed: 0            id         price      bedrooms     bathrooms  \\\n",
       "count  21613.00000  2.161300e+04  2.161300e+04  21600.000000  21603.000000   \n",
       "mean   10806.00000  4.580302e+09  5.400881e+05      3.372870      2.115736   \n",
       "std     6239.28002  2.876566e+09  3.671272e+05      0.926657      0.768996   \n",
       "min        0.00000  1.000102e+06  7.500000e+04      1.000000      0.500000   \n",
       "25%     5403.00000  2.123049e+09  3.219500e+05      3.000000      1.750000   \n",
       "50%    10806.00000  3.904930e+09  4.500000e+05      3.000000      2.250000   \n",
       "75%    16209.00000  7.308900e+09  6.450000e+05      4.000000      2.500000   \n",
       "max    21612.00000  9.900000e+09  7.700000e+06     33.000000      8.000000   \n",
       "\n",
       "        sqft_living      sqft_lot        floors    waterfront          view  \\\n",
       "count  21613.000000  2.161300e+04  21613.000000  21613.000000  21613.000000   \n",
       "mean    2079.899736  1.510697e+04      1.494309      0.007542      0.234303   \n",
       "std      918.440897  4.142051e+04      0.539989      0.086517      0.766318   \n",
       "min      290.000000  5.200000e+02      1.000000      0.000000      0.000000   \n",
       "25%     1427.000000  5.040000e+03      1.000000      0.000000      0.000000   \n",
       "50%     1910.000000  7.618000e+03      1.500000      0.000000      0.000000   \n",
       "75%     2550.000000  1.068800e+04      2.000000      0.000000      0.000000   \n",
       "max    13540.000000  1.651359e+06      3.500000      1.000000      4.000000   \n",
       "\n",
       "       ...         grade    sqft_above  sqft_basement      yr_built  \\\n",
       "count  ...  21613.000000  21613.000000   21613.000000  21613.000000   \n",
       "mean   ...      7.656873   1788.390691     291.509045   1971.005136   \n",
       "std    ...      1.175459    828.090978     442.575043     29.373411   \n",
       "min    ...      1.000000    290.000000       0.000000   1900.000000   \n",
       "25%    ...      7.000000   1190.000000       0.000000   1951.000000   \n",
       "50%    ...      7.000000   1560.000000       0.000000   1975.000000   \n",
       "75%    ...      8.000000   2210.000000     560.000000   1997.000000   \n",
       "max    ...     13.000000   9410.000000    4820.000000   2015.000000   \n",
       "\n",
       "       yr_renovated       zipcode           lat          long  sqft_living15  \\\n",
       "count  21613.000000  21613.000000  21613.000000  21613.000000   21613.000000   \n",
       "mean      84.402258  98077.939805     47.560053   -122.213896    1986.552492   \n",
       "std      401.679240     53.505026      0.138564      0.140828     685.391304   \n",
       "min        0.000000  98001.000000     47.155900   -122.519000     399.000000   \n",
       "25%        0.000000  98033.000000     47.471000   -122.328000    1490.000000   \n",
       "50%        0.000000  98065.000000     47.571800   -122.230000    1840.000000   \n",
       "75%        0.000000  98118.000000     47.678000   -122.125000    2360.000000   \n",
       "max     2015.000000  98199.000000     47.777600   -121.315000    6210.000000   \n",
       "\n",
       "          sqft_lot15  \n",
       "count   21613.000000  \n",
       "mean    12768.455652  \n",
       "std     27304.179631  \n",
       "min       651.000000  \n",
       "25%      5100.000000  \n",
       "50%      7620.000000  \n",
       "75%     10083.000000  \n",
       "max    871200.000000  \n",
       "\n",
       "[8 rows x 21 columns]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 2.0 Data Wrangling"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Question 2 \n",
    "Drop the columns <code>\"id\"</code>  and <code>\"Unnamed: 0\"</code> from axis 1 using the method <code>drop()</code>, then use the method <code>describe()</code> to obtain a statistical summary of the data. Take a screenshot and submit it, make sure the inplace parameter is set to <code>True</code>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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>price</th>\n",
       "      <th>bedrooms</th>\n",
       "      <th>bathrooms</th>\n",
       "      <th>sqft_living</th>\n",
       "      <th>sqft_lot</th>\n",
       "      <th>floors</th>\n",
       "      <th>waterfront</th>\n",
       "      <th>view</th>\n",
       "      <th>condition</th>\n",
       "      <th>grade</th>\n",
       "      <th>sqft_above</th>\n",
       "      <th>sqft_basement</th>\n",
       "      <th>yr_built</th>\n",
       "      <th>yr_renovated</th>\n",
       "      <th>zipcode</th>\n",
       "      <th>lat</th>\n",
       "      <th>long</th>\n",
       "      <th>sqft_living15</th>\n",
       "      <th>sqft_lot15</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>2.161300e+04</td>\n",
       "      <td>21600.000000</td>\n",
       "      <td>21603.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>2.161300e+04</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "      <td>21613.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>5.400881e+05</td>\n",
       "      <td>3.372870</td>\n",
       "      <td>2.115736</td>\n",
       "      <td>2079.899736</td>\n",
       "      <td>1.510697e+04</td>\n",
       "      <td>1.494309</td>\n",
       "      <td>0.007542</td>\n",
       "      <td>0.234303</td>\n",
       "      <td>3.409430</td>\n",
       "      <td>7.656873</td>\n",
       "      <td>1788.390691</td>\n",
       "      <td>291.509045</td>\n",
       "      <td>1971.005136</td>\n",
       "      <td>84.402258</td>\n",
       "      <td>98077.939805</td>\n",
       "      <td>47.560053</td>\n",
       "      <td>-122.213896</td>\n",
       "      <td>1986.552492</td>\n",
       "      <td>12768.455652</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>3.671272e+05</td>\n",
       "      <td>0.926657</td>\n",
       "      <td>0.768996</td>\n",
       "      <td>918.440897</td>\n",
       "      <td>4.142051e+04</td>\n",
       "      <td>0.539989</td>\n",
       "      <td>0.086517</td>\n",
       "      <td>0.766318</td>\n",
       "      <td>0.650743</td>\n",
       "      <td>1.175459</td>\n",
       "      <td>828.090978</td>\n",
       "      <td>442.575043</td>\n",
       "      <td>29.373411</td>\n",
       "      <td>401.679240</td>\n",
       "      <td>53.505026</td>\n",
       "      <td>0.138564</td>\n",
       "      <td>0.140828</td>\n",
       "      <td>685.391304</td>\n",
       "      <td>27304.179631</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>7.500000e+04</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>290.000000</td>\n",
       "      <td>5.200000e+02</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>290.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1900.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>98001.000000</td>\n",
       "      <td>47.155900</td>\n",
       "      <td>-122.519000</td>\n",
       "      <td>399.000000</td>\n",
       "      <td>651.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>3.219500e+05</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>1.750000</td>\n",
       "      <td>1427.000000</td>\n",
       "      <td>5.040000e+03</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>7.000000</td>\n",
       "      <td>1190.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1951.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>98033.000000</td>\n",
       "      <td>47.471000</td>\n",
       "      <td>-122.328000</td>\n",
       "      <td>1490.000000</td>\n",
       "      <td>5100.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>4.500000e+05</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>2.250000</td>\n",
       "      <td>1910.000000</td>\n",
       "      <td>7.618000e+03</td>\n",
       "      <td>1.500000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>7.000000</td>\n",
       "      <td>1560.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1975.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>98065.000000</td>\n",
       "      <td>47.571800</td>\n",
       "      <td>-122.230000</td>\n",
       "      <td>1840.000000</td>\n",
       "      <td>7620.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>6.450000e+05</td>\n",
       "      <td>4.000000</td>\n",
       "      <td>2.500000</td>\n",
       "      <td>2550.000000</td>\n",
       "      <td>1.068800e+04</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>4.000000</td>\n",
       "      <td>8.000000</td>\n",
       "      <td>2210.000000</td>\n",
       "      <td>560.000000</td>\n",
       "      <td>1997.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>98118.000000</td>\n",
       "      <td>47.678000</td>\n",
       "      <td>-122.125000</td>\n",
       "      <td>2360.000000</td>\n",
       "      <td>10083.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>7.700000e+06</td>\n",
       "      <td>33.000000</td>\n",
       "      <td>8.000000</td>\n",
       "      <td>13540.000000</td>\n",
       "      <td>1.651359e+06</td>\n",
       "      <td>3.500000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>4.000000</td>\n",
       "      <td>5.000000</td>\n",
       "      <td>13.000000</td>\n",
       "      <td>9410.000000</td>\n",
       "      <td>4820.000000</td>\n",
       "      <td>2015.000000</td>\n",
       "      <td>2015.000000</td>\n",
       "      <td>98199.000000</td>\n",
       "      <td>47.777600</td>\n",
       "      <td>-121.315000</td>\n",
       "      <td>6210.000000</td>\n",
       "      <td>871200.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              price      bedrooms     bathrooms   sqft_living      sqft_lot  \\\n",
       "count  2.161300e+04  21600.000000  21603.000000  21613.000000  2.161300e+04   \n",
       "mean   5.400881e+05      3.372870      2.115736   2079.899736  1.510697e+04   \n",
       "std    3.671272e+05      0.926657      0.768996    918.440897  4.142051e+04   \n",
       "min    7.500000e+04      1.000000      0.500000    290.000000  5.200000e+02   \n",
       "25%    3.219500e+05      3.000000      1.750000   1427.000000  5.040000e+03   \n",
       "50%    4.500000e+05      3.000000      2.250000   1910.000000  7.618000e+03   \n",
       "75%    6.450000e+05      4.000000      2.500000   2550.000000  1.068800e+04   \n",
       "max    7.700000e+06     33.000000      8.000000  13540.000000  1.651359e+06   \n",
       "\n",
       "             floors    waterfront          view     condition         grade  \\\n",
       "count  21613.000000  21613.000000  21613.000000  21613.000000  21613.000000   \n",
       "mean       1.494309      0.007542      0.234303      3.409430      7.656873   \n",
       "std        0.539989      0.086517      0.766318      0.650743      1.175459   \n",
       "min        1.000000      0.000000      0.000000      1.000000      1.000000   \n",
       "25%        1.000000      0.000000      0.000000      3.000000      7.000000   \n",
       "50%        1.500000      0.000000      0.000000      3.000000      7.000000   \n",
       "75%        2.000000      0.000000      0.000000      4.000000      8.000000   \n",
       "max        3.500000      1.000000      4.000000      5.000000     13.000000   \n",
       "\n",
       "         sqft_above  sqft_basement      yr_built  yr_renovated       zipcode  \\\n",
       "count  21613.000000   21613.000000  21613.000000  21613.000000  21613.000000   \n",
       "mean    1788.390691     291.509045   1971.005136     84.402258  98077.939805   \n",
       "std      828.090978     442.575043     29.373411    401.679240     53.505026   \n",
       "min      290.000000       0.000000   1900.000000      0.000000  98001.000000   \n",
       "25%     1190.000000       0.000000   1951.000000      0.000000  98033.000000   \n",
       "50%     1560.000000       0.000000   1975.000000      0.000000  98065.000000   \n",
       "75%     2210.000000     560.000000   1997.000000      0.000000  98118.000000   \n",
       "max     9410.000000    4820.000000   2015.000000   2015.000000  98199.000000   \n",
       "\n",
       "                lat          long  sqft_living15     sqft_lot15  \n",
       "count  21613.000000  21613.000000   21613.000000   21613.000000  \n",
       "mean      47.560053   -122.213896    1986.552492   12768.455652  \n",
       "std        0.138564      0.140828     685.391304   27304.179631  \n",
       "min       47.155900   -122.519000     399.000000     651.000000  \n",
       "25%       47.471000   -122.328000    1490.000000    5100.000000  \n",
       "50%       47.571800   -122.230000    1840.000000    7620.000000  \n",
       "75%       47.678000   -122.125000    2360.000000   10083.000000  \n",
       "max       47.777600   -121.315000    6210.000000  871200.000000  "
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.drop([\"id\",\"Unnamed: 0\"], axis = 1,inplace = True) \n",
    "df.describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "we can see we have missing values for the columns <code> bedrooms</code>  and <code> bathrooms </code>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "number of NaN values for the column bedrooms : 13\n",
      "number of NaN values for the column bathrooms : 10\n"
     ]
    }
   ],
   "source": [
    "print(\"number of NaN values for the column bedrooms :\", df['bedrooms'].isnull().sum())\n",
    "print(\"number of NaN values for the column bathrooms :\", df['bathrooms'].isnull().sum())\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "We can replace the missing values of the column <code>'bedrooms'</code> with the mean of the column  <code>'bedrooms' </code> using the method replace. Don't forget to set the <code>inplace</code> parameter top <code>True</code>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean=df['bedrooms'].mean()\n",
    "df['bedrooms'].replace(np.nan,mean, inplace=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "We also replace the missing values of the column <code>'bathrooms'</code> with the mean of the column  <code>'bedrooms' </codse> using the method replace.Don't forget to set the <code> inplace </code>  parameter top <code> Ture </code>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean=df['bathrooms'].mean()\n",
    "df['bathrooms'].replace(np.nan,mean, inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "number of NaN values for the column bedrooms : 0\n",
      "number of NaN values for the column bathrooms : 0\n"
     ]
    }
   ],
   "source": [
    "print(\"number of NaN values for the column bedrooms :\", df['bedrooms'].isnull().sum())\n",
    "print(\"number of NaN values for the column bathrooms :\", df['bathrooms'].isnull().sum())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 3.0 Exploratory data analysis"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Question 3\n",
    "Use the method value_counts to count the number of houses with unique floor values, use the method .to_frame() to convert it to a dataframe.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "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>floors</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1.0</th>\n",
       "      <td>10680</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2.0</th>\n",
       "      <td>8241</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1.5</th>\n",
       "      <td>1910</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3.0</th>\n",
       "      <td>613</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2.5</th>\n",
       "      <td>161</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3.5</th>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     floors\n",
       "1.0   10680\n",
       "2.0    8241\n",
       "1.5    1910\n",
       "3.0     613\n",
       "2.5     161\n",
       "3.5       8"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['floors'].value_counts().to_frame()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 4\n",
    "Use the function <code>boxplot</code> in the seaborn library  to  determine whether houses with a waterfront view or without a waterfront view have more price outliers ."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7f57750ffac8>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAEKCAYAAAC7c+rvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3X10XNV57/HvI1sG8xaDZWzHMjGNRRKaRR3Qsl0SCG8yFjfY3HtDrrldYZLSOCW85brNDenKCi+ha5HV5qV2WmIK3AxdacDQtAhiYYR5bcOLZV5MgAQrRMBgwJKNHYMBS9Zz/5gtMSOk0cxozpyZ0e+z1izN3mefs/fxWPNo77PPPubuiIiIxKEu7gaIiMjEpSAkIiKxURASEZHYKAiJiEhsFIRERCQ2CkIiIhIbBSEREYmNgpCIiMQm0iBkZv/HzJ41s1+b2c/N7EAzO9rMHjOzrWZ2q5lNCWUPCOmusH1exnG+FfJ/a2ZnZuQvDXldZnZ5Rn7BdYiISPlZVCsmmNkc4D+BY939HTNbB6wHzgJ+4e63mNlPgKfd/Toz+xpwnLv/pZmtAP67u/8vMzsW+DmwEPgwcC9wTKjmBaAFSAGbgPPc/blQV9515DqPhoYGnzdvXkn/bUREat3mzZt73X3GWOUmR9yOycBUM+sDDgJeA04D/nfYngSuBK4Dlof3ALcDPzYzC/m3uPt7wO/NrIt0QALocvcXAczsFmC5mT1faB2eIxLPmzePzs7OIk9fRGRiMrOX8ikX2XCcu78K/D3wMungsxvYDOxy9/5QLAXMCe/nAK+EfftD+emZ+cP2GS1/ehF1ZDGzlWbWaWadPT09xZy+iIjkIbIgZGaHk+55HE16GO1goHWEooO9EBtlW6nyc9WRneF+vbs3u3vzjBlj9iZFRKRIUU5MOAP4vbv3uHsf8AvgRGCamQ0OAzYC28L7FDAXIGz/ELAzM3/YPqPl9xZRh4iIxCDKIPQysNjMDgrXdk4HngPuBz4fyiSAO8L7tpAmbL8vXKtpA1aEmW1HA03A46QnIjSFmXBTgBVAW9in0DpERCQGUV4Teoz0xf8ngGdCXdcD3wRWhQkG04Ebwy43AtND/irg8nCcZ4F1pAPY3cBF7r4/XNO5GNgAPA+sC2UptA4Rya23t5dLLrmEHTt2xN0UqTGRTdGuFc3Nza7ZcTLRff/736etrY3ly5ezatWquJsjVcDMNrt781jltGKCiOTU29tLe3s77k57e7t6Q1JSCkIiklMymWRwxGRgYIBkMhlzi6SWKAiJSE4dHR309fUB0NfXxz333BNzi6SWKAiJSE4tLS3U19cDUF9fz5IlS2JukdQSBSERySmRSJC+ywLq6upIJBJj7CGSPwUhEcmpoaGB1tZWzIzW1lamT//ASlciRYt6AVMRqQGJRILu7m71gqTkFIREZEwNDQ2sWbMm7mZIDdJwnIiIxEZBSEREYqMgJCIisVEQEhGR2CgIiYhIbBSEREQkNgpCIiISGwUhkRLRg99ECqcgJFIiyWSSLVu26FEHIgWILAiZ2cfM7KmM1x/M7OtmdoSZdZjZ1vDz8FDezGy1mXWZ2RYzOz7jWIlQfquZJTLyTzCzZ8I+qy2sslhMHSLjoQe/iRQnsiDk7r919wXuvgA4AdgL/DtwObDR3ZuAjSEN0Ao0hddK4DpIBxTgCmARsBC4YjCohDIrM/ZbGvILqkNkvPTgN5HilGs47nTgd+7+ErAcGPwNTQLnhPfLgZs97VFgmpnNBs4EOtx9p7u/CXQAS8O2w9z9EU//9t887FiF1CEyLnrwm0hxyhWEVgA/D+9nuvtrAOHnkSF/DvBKxj6pkJcrPzVCfjF1ZDGzlWbWaWadPT09BZymTFQtLS1Dz9wxMz34TSRPkQchM5sCLANuG6voCHleRH4xdWRnuF/v7s3u3jxjxowxDikCZ5999tBwnLuzbNmymFskUh3K0RNqBZ5w9zdC+o3BIbDwc3vITwFzM/ZrBLaNkd84Qn4xdYiMy5133pnVE2pra4u5RSLVoRxB6DzeH4oDaAMGZ7glgDsy8s8PM9gWA7vDUNoGYImZHR4mJCwBNoRte8xscZgVd/6wYxVSh8i4dHR0ZPWEdE1IJD+RBiEzOwhoAX6RkX0t0GJmW8O2a0P+euBFoAv4Z+BrAO6+E/gusCm8rg55ABcCN4R9fge0F1OHyHi1tLRQX18PQH19va4JieTJBv96k5E1Nzd7Z2dn3M2QCtfb28uKFSvYt28fBxxwALfccgvTp0+Pu1kisTGzze7ePFY5rZggUgINDQ20trZiZrS2tioAieRpctwNEKkViUSC7u5uEonE2IVFBFAQEimZhoYG1qxZE3czRKqKhuNERCQ2CkIiIhIbBSEREYmNgpCIiMRGQUhERGKjICQiIrFREBIRkdgoCImISGwUhEREJDYKQiIiEhsFIRERiY2CkIiIxEZBSEREYqMgJCIisYn68d7TzOx2M/uNmT1vZn9qZkeYWYeZbQ0/Dw9lzcxWm1mXmW0xs+MzjpMI5beaWSIj/wQzeybss9rMLOQXXIeIiJRf1D2hfwDudvePA38CPA9cDmx09yZgY0gDtAJN4bUSuA7SAQW4AlgELASuGAwqoczKjP2WhvyC6hARkXhEFoTM7DDgZOBGAHff5+67gOVAMhRLAueE98uBmz3tUWCamc0GzgQ63H2nu78JdABLw7bD3P0Rd3fg5mHHKqQOERGJQZQ9oT8CeoD/Z2ZPmtkNZnYwMNPdXwMIP48M5ecAr2Tsnwp5ufJTI+RTRB1ZzGylmXWaWWdPT09hZy0iInmLMghNBo4HrnP3TwFv8/6w2EhshDwvIj+XvPZx9+vdvdndm2fMmDHGIUVEpFhRBqEUkHL3x0L6dtJB6Y3BIbDwc3tG+bkZ+zcC28bIbxwhnyLqEBGRGEQWhNz9deAVM/tYyDodeA5oAwZnuCWAO8L7NuD8MINtMbA7DKVtAJaY2eFhQsISYEPYtsfMFodZcecPO1YhdYiISAwmR3z8S4CfmdkU4EXgy6QD3zozuwB4GTg3lF0PnAV0AXtDWdx9p5l9F9gUyl3t7jvD+wuBnwJTgfbwAri2kDpERCQelp5YJqNpbm72zs7OuJshIlJVzGyzuzePVU4rJoiISGwUhERKpLe3l0suuYQdO3bE3RSRqqEgJFIiyWSSLVu2kEwmxy4sIoCCkEhJ9Pb20t7ejrvT3t6u3pBInhSEREogmUwyOMlnYGBAvSGRPCkIiZRAR0cHfX19APT19XHPPffE3CKR6qAgJFICLS0t1NfXA1BfX8+SJUtibpFIdVAQEimBRCJBeJwVZkYikRhjDxEBBSGRkmhoaGDmzJkAzJw5k+nTp8fcIpHqoCAkUgK9vb2kUukni6RSKc2OE8mTgpBICaxdu3Zodpy7s3bt2phbJFIdFIRESmDjxo1Z6XvvvTemlohUFwUhkRIYvhCwFgYWyY+CkEgJnHHGGVnplpaWmFoiUl0UhERK4Ktf/WrOtIiMTEFIpEQy7xMSkfwoCImUQDKZZNKkSQBMmjRJa8eJ5CnSIGRm3Wb2jJk9ZWadIe8IM+sws63h5+Eh38xstZl1mdkWMzs+4ziJUH6rmSUy8k8Ix+8K+1qxdYiMR0dHB/39/QD09/dr7TiRPJWjJ3Squy/IeMzr5cBGd28CNoY0QCvQFF4rgesgHVCAK4BFwELgisGgEsqszNhvaTF1iIyX1o4TKU4cw3HLgcGxiiRwTkb+zZ72KDDNzGYDZwId7r7T3d8EOoClYdth7v6Ip+fD3jzsWIXUITIumWvH1dXVae04kTxFHYQcuMfMNpvZypA3091fAwg/jwz5c4BXMvZNhbxc+akR8oupI4uZrTSzTjPr7OnpKeB0ZaJqaGigtbUVM6O1tVVrx4nkaXLEx/+0u28zsyOBDjP7TY6yI00p8iLyc8lrH3e/HrgeoLm5WXcdSl4SiQTd3d3qBYkUINKekLtvCz+3A/9O+prOG4NDYOHn9lA8BczN2L0R2DZGfuMI+RRRh8i4NTQ0sGbNGvWCRAoQWRAys4PN7NDB98AS4NdAGzD4p2ICuCO8bwPODzPYFgO7w1DaBmCJmR0eJiQsATaEbXvMbHGYFXf+sGMVUoeIiMQgyuG4mcC/h4u1k4F/dfe7zWwTsM7MLgBeBs4N5dcDZwFdwF7gywDuvtPMvgtsCuWudved4f2FwE+BqUB7eAFcW0gdIiISD9NCi7k1Nzd7Z2dnWerq7e3lqquu4sorr9SQjohUNTPbnHFrzqi0YkIFSSaTbNmyRXfbi8iEoSBUIXp7e2lvb8fdaW9v15M5RWRCUBCqEMlkcugZNAMDA+oNiciEoCBUITo6Oujr6wOgr69Pa49VoRdeeIHW1la6urribopI1VAQqhBae6z6XXPNNbz99ttcffXVcTdFpGooCFUIrT1W3V544QW6u7sB6O7uVm9IJE8KQhVCa49Vt2uuuSYrrd6QSH6iXjtOCqC1x6rXYC9otLSIjEw9oQqitceq17x583KmRWRkCkIiJfDtb387K/2d73wnppaIVBcFIZESOOaYY4Z6P/PmzWP+/PnxNkikSigIiZTIxRdfTF1dHZdddlncTRGpGgpCIiXy8MMP4+48+OCDcTdFpGrkHYTM7CNmdkZ4P3XwWUEiorX/RIqVVxAys68AtwNrQ1Yj8B9RNUqk2mjtP5Hi5NsTugj4NPAHAHffChwZVaNEqo3W/hMpTr5B6D133zeYMLPJgJ6GJxK0tLQweXL63u/Jkydr7T+RPOUbhB40s78BpppZC3AbcGc+O5rZJDN70szuCumjzewxM9tqZrea2ZSQf0BId4Xt8zKO8a2Q/1szOzMjf2nI6zKzyzPyC65DZDwSiQQDAwNAejhOq16I5CffIHQ50AM8A3wVWA98O+ce77sMeD4j/T3gh+7eBLwJXBDyLwDedPf5wA9DOczsWGAF8MfAUuCfQmCbBPwj0AocC5wXyhZch4iIxCPfIDQVuMndz3X3zwM3hbyczKwR+G/ADSFtwGmkJzkAJIFzwvvlIU3Yfnoovxy4xd3fc/ffA13AwvDqcvcXw1DhLcDyIusQGZdkMkldXfrXqa6uThMTRPKUbxDaSHbQmQrcm8d+PwL+LzAQ0tOBXe7eH9IpYE54Pwd4BSBs3x3KD+UP22e0/GLqyGJmK82s08w6e3p68jhNmeg6Ojro70//l+vv79fEBJE85RuEDnT3twYT4f1BuXYws88B2919c2b2CEV9jG2lyh+r/vcz3K9392Z3b54xY8YIu0Sjt7eXSy65RPeYVKGWlpastCYmiOQn3yD0tpkdP5gwsxOAd8bY59PAMjPrJj1UdhrpntG0MLsO0vcbbQvvU8DccPzJwIeAnZn5w/YZLb+3iDoqQjKZZMuWLRrKqUJnn312VnrZsmUxtUSkuuQbhL4O3GZmD5vZw8CtwMW5dnD3b7l7o7vPIz2x4D53/zPgfuDzoVgCuCO8bwtpwvb7PH33XxuwIsxsOxpoAh4HNgFNYSbclFBHW9in0Dpipzvuq9ttt92WlV63bl1MLRGpLnkFIXffBHwcuBD4GvCJYcNshfgmsMrMukhfj7kx5N8ITA/5q0jPyMPdnwXWAc8BdwMXufv+cE3nYmAD6dl360LZguuoBLrjvrrde2/2JdKOjo6YWiJSXSxXR8DMTnP3+8zsf4y03d1/EVnLKkRzc7N3dnZGXs/SpUvZu3fvUPqggw7i7rvvjrzecurt7eWqq67iyiuvrLkH95122mlDExMgfcPqfffdF2OLROJlZpvdvXmscmP1hD4bfp49wutz42qhZGlpaaG+vh6A+vr6mrywXcvXvD7zmc9kpU866aSYWiJSXXIGIXe/wszqgHZ3//Kw15+XqY0TQiKRYPCWpbq6upq7477Wr3m99957OdMiMrIxrwm5+wBjTEKQ8WtoaKC1tRUzo7W1teaGq2r9mtcjjzySlf7Vr34VU0tEqku+s+M6zOyvzWyumR0x+Iq0ZRNQIpHguOOOq7leEGiVaREZ2eSxiwDw56Rv6vzasPw/Km1zJraGhgbWrFkTdzMi0dLSwvr16+nr66vZa14iUrh8e0LHkl4s9GngKWAN6QVFRfJS69e8RKQ4+QahJPAJYDXpAPQJ3l8IVEqklpftqfVrXnPnzs2ZFpGR5RuEPubuf+Hu94fXSuBjUTZsIlq7di1PP/00a9euHbtwFarla16XXXZZVnrVqlUxtUSkuuQbhJ40s8WDCTNbBPxXNE2amHp7e4cu1t9zzz012xtas2ZNzfWCAO66666sdFtbW0wtEfmgSh5lyTcILQJ+ZWbdYUHSR4DPmtkzZrYlstZNIGvXrs2awlyrvaFa9cADD2Sl77///ngaEpFK/hKTsVXyjeL5BqGlwNGkV1D4bHh/FulVE87OsZ/kSWuPSSWr5C8xya3SbxTPdwHTl3K9om7kRDAwMJAzLZVt8uTJOdPVrLe3l1/+8pe4O7/85S8r7ktMcqv0G8Xz7QlJxAYfDT1aWipb5uKlI6WrWTKZHDqfvr6+ivsSk9wq/UZxfdNViIULF2alFy1aFFNLpBiNjY0509Vsw4YNWelaW90davuaV6UvjqwgVCFeeil7VLO7uzuehkhR5s+fn5VuamqKqSWlV8tDjYNq+ZpXpd8oriBUIbZt25YzLZWtlhcwfeutt3Kmq12lX7gfr0q/UTyyIGRmB5rZ42b2tJk9a2ZXhfyjzewxM9tqZreGR3MTHt99q5l1he3zMo71rZD/WzM7MyN/acjrMrPLM/ILrkNkPPbv358zXc1qeagRKv/CfSlU8o3iUfaE3gNOc/c/ARYAS8MNr98DfujuTcCbwAWh/AXAm+4+H/hhKIeZHQusIL1W3VLgn8xskplNIr2eXSvpte3OC2UptI5KMHXq1JxpqWy1PDGhlocaofIv3JdCJd8oHlkQ8rTBfnt9eDlwGnB7yE8C54T3y3l/PbrbgdMtPZC5HLjF3d9z998DXcDC8Opy9xfdfR9wC7A87FNoHbF75513cqalsh1wwAE509Xs8ccfz0o/9thjMbUkGi0tLUPXuSZPnlxxF+5rXaTXhEKP5SlgO9AB/A7Y5e6DfyamgDnh/RzgFYCwfTcwPTN/2D6j5U8voo7YHXzwwTnTtaCWZyDV8pNVW1pastK19iWdSCSG7ssbGBioyCGrWhZpEHL3/e6+AGgk3XP5xEjFws+ReiRewvxcdWQxs5Vm1mlmnT09PSPsUnrvvvtuznQtqOUZSLXs7LOzF0VZtmxZTC2RWlSW2XHuvgt4AFgMTDOzwTmejcDgNLAUMBcgbP8QsDMzf9g+o+X3FlHH8PZe7+7N7t48Y8aM4k5astT6DKRadtttt2Wl161bF1NLopFMJoduDq+rq9MfSWUW5ey4GWY2LbyfCpwBPA/cD3w+FEsAd4T3bSFN2H6fp6estAErwsy2o4Em4HFgE9AUZsJNIT15oS3sU2gdsTvppJOy0ieffHJMLYnGRJiBVKs2btyYlR6+zmG16+joGJpI0t/fX5MTEypZlD2h2cD9YZXtTUCHu98FfBNYZWZdpK/H3BjK3whMD/mrgMsB3P1ZYB3wHHA3cFEY5usHLgY2kA5u60JZCq2jEtTyhW2YGDOQatXwv9Mq5O+2kqn0FQVqXZSz47a4+6fc/Th3/6S7Xx3yX3T3he4+393Pdff3Qv67IT0/bH8x41h/6+4fdfePuXt7Rv56dz8mbPvbjPyC64jbww8/nJV+6KGHYmpJNPSLXr3OOOOMrPTwiQrVrtJXFKh1WjGhQtT6NNFa/0UfPtO/Qmb+l8S5556blf7CF74QU0uiUekrCtQ6BaEKkUgkhi6OTpo0qea+pGv9F72Wh6zuvPPOrHQtPjW2klcUqHUKQhWioaGBU089FYBTTz215r6kQb/o1Wr4AxZr8XpeJa8oUOsUhCpIb28vQM1OX9YvenWq9ZmbEi8FoQrR29vL5s2bAdi0aVNNBqJaXjFBRIqjIFQh/u7v/i5nuhasXbuWp59+mrVr18bdFCnA8JmaDz74YEwtkVqkIFQhavl5NJDuBQ1eS9iwYYN6Q1Vk5syZOdMi46EgJGWxdu3aoRlj7q7eUBV5/fXXc6ZFxkNBSMpi+FIvw2dcSeWaNWtWzrTIeCgISVnU8pNHa90bb7yRMy0yHgpCUha1vKJArVuyZMnQ52VmnHnmmTG3qPQ0czM+k8cuIjJ+dXV1Wb2fwdUhJD6rV6+mq6trzHJ9fX1Z1/O2bt3KpZdeOuZ+8+fPz6tcJch81tWqVavibs6Eom8CKYtaXwSzltXX1w+tazh9+vShhWhrhZ51FS/1hKQsWlpa2LBhw1C6Fod0qk0hvZQLL7yQ7u5ubrjhhppb8SKZTA493nv//v3qDZWZekJSFj/60Y+y0t///vdjaokUo76+nqamppoLQKCH2sVNQUjKIpVK5UyLxEVr48VLQUhERGITWRAys7lmdr+ZPW9mz5rZZSH/CDPrMLOt4efhId/MbLWZdZnZFjM7PuNYiVB+q5klMvJPMLNnwj6rLcwjLaaOuH34wx/Oma52kyZNypkWiUutP9W40kXZE+oH/srdPwEsBi4ys2OBy4GN7t4EbAxpgFagKbxWAtdBOqAAVwCLgIXAFYNBJZRZmbHf0pBfUB2VYOfOnTnT1W7hwoVZ6UWLFsXUEpFsw2dq1tpTjStdZEHI3V9z9yfC+z3A88AcYDmQDMWSwDnh/XLgZk97FJhmZrOBM4EOd9/p7m8CHcDSsO0wd3/E0zcx3DzsWIXUEbtDDjkkZ7raDb8G9Morr8TUEpFsZ599dlZ62bJlMbVkYirLNSEzmwd8CngMmOnur0E6UAFHhmJzgMxvplTIy5WfGiGfIuqI3eAD7UZLV7vhQUdBSCrFbbfdlpVet25dTC2ZmCIPQmZ2CPBvwNfd/Q+5io6Q50Xk52xOPvuY2Uoz6zSzzp6enjEOKfkYvNlxtLRIXLS4brwiDUJmVk86AP3M3X8Rst8YHAILP7eH/BQwN2P3RmDbGPmNI+QXU0cWd7/e3ZvdvXnGjBn5n7CMavA+jNHSIjIxRTk7zoAbgefd/QcZm9qAwRluCeCOjPzzwwy2xcDuMJS2AVhiZoeHCQlLgA1h2x4zWxzqOn/YsQqpQyJW69e8pHodeeSRWWk9tK+8ohwT+TTwReAZM3sq5P0NcC2wzswuAF4Gzg3b1gNnAV3AXuDLAO6+08y+C2wK5a5298GpYxcCPwWmAu3hRaF1SPT6+vpypkXisn379qy0HlVRXpEFIXf/T0a+BgNw+gjlHbholGPdBNw0Qn4n8MkR8ncUWodEa/bs2XR3d2elRUS0YoKUxauvvpqV3rbtA5fiRGIxfDhueFqipSAkZTF8IoKG46RSvP766znTEi0FISmLwYeijZYWiYv+b8ZLQUjKQvcJichIFISkLBobG7PSRx11VEwtEZFKoiAkZZE5Mw7gxRdfjKchIlJRFIRERCQ2CkIiIhIbBSEpiylTpuRMi8jEpCAkZbFv376caRGZmBSEREQkNrpZQ2QUq1evpqurq+j9L7300pzb58+fP2YZkVqnICRSQ8YbOEezdetWYOzAWiwF5IlLQUhkFIV8KZ588skfyFu9enUpm5OXrq4ufvPUU8wq8XEHx+13PfVUznLF0EptE5uCkEgJPPTQQ1mB6KGHHoqtLbOAC0Z9ikrluRGt1Ra13t5errrqKq688kqmT58ed3OyKAjJuIxn+EfXTETKI5lMsmXLFpLJJKtWrYq7OVk0O06kRBYsWMCCBQti7QWJDNfb20t7ezvuTnt7Ozt27Ii7SVki6wmZ2U3A54Dt7v7JkHcEcCswD+gGvuDub5qZAf9A+tHbe4EvufsTYZ8E8O1w2GvcPRnyT+D9R3uvBy5zdy+mDilevj2VSrlmIjLRJJPJocdTDAwMVFxvKMqe0E+BpcPyLgc2unsTsDGkAVqBpvBaCVwHQ0HrCmARsBC4wswOD/tcF8oO7re0mDqkPL7yla9kpS+88MKYWiIysXR0dAw9RLKvr4977rkn5hZliywIuftDwM5h2cuBZHifBM7JyL/Z0x4FppnZbOBMoMPdd7r7m0AHsDRsO8zdH/F0iL952LEKqUPK4Itf/GJW+rzzzoupJSITS0tLC/X19QDU19ezZMmSmFuUrdwTE2a6+2sA7v6amQ0+zH0O8EpGuVTIy5WfGiG/mDpeG+9JSX5mzZrF66+/rl6QlIUmzaQlEgna29sBqKurI5FIxNyibJUyMWGk+aReRH4xdXywoNlKM+s0s86enp4xDiv5mjVrFgsWLFAvSKSMGhoaaG1txcxobW2d8FO03zCz2aGHMhvYHvJTwNyMco3AtpB/yrD8B0J+4wjli6njA9z9euB6gObm5nHdxKC/xkTike/vxje/+U0eeeSRofSJJ57ItddeG1WzYpFIJOju7q64XhCUvyfUBgz+KySAOzLyz7e0xcDuMKS2AVhiZoeHCQlLgA1h2x4zWxxmvZ0/7FiF1CEiE9g3vvGNnOla0NDQwJo1ayquFwTRTtH+OeleTIOZpUjPcrsWWGdmFwAvA+eG4utJT53uIj19+ssA7r7TzL4LbArlrnb3wckOF/L+FO328KLQOqKmKcwila2hoYFDDz2UPXv2cOKJJ1bkF3UtiywIuftoA/+nj1DWgYtGOc5NwE0j5HcCnxwhf0ehdYjUilQqxR6qaymc14C3Uqkxy0XpqKOOoru7uyZ7QZVOy/ZUiEpae6yWRbXKNES70rSu/0Wrvr6epqYm9YJioCAkE0pXVxdPPvskTIvg4APpH0+++mRpj7sr/6KNjY3s6u2tugVMpzU2jl1QapKCUAVZsGABEO+1oAnRU5gGA6cMlLwNUal7oFLupJCoFfL7l0qleOeddyJpx9SpU2nM8w+D8fbSFYQkS1dXFy/8+gmOOmR/yY89pS/9Zfpu96YxShbm5bcmlfR4Urmq8aF9hXxJd3V18eunn+bQKWN/Ne/t38/+gWiu/fW/+w4v7dk9Zrk9+/rHXZeCkHzAUYfs59vNb8XdjLxd03lI3E2QMunq6uLZZ55n2kFHjl0qes2NAAAIRUlEQVS4AAP70sOXr/6utCtM79q7fexCwxw6ZTILZx4+dsEK8Pgbb477GApCIlJVph10JKd+fEXczcjL/b+5Je4mVDwFoSLU8pBAKpXi7T2Tqqp38dKeSRyc5xTfVCoFu6vsOssuSHn+U5hfp/RTtAf7B1HMHXudaOaJVKNUKsWeff0l6WGUw559/enfqXFQECpCV1cXTz7zHAMHHVHS49q+9BfH5t+9XtLj1u0dvpi51Kr58+dHctye8AfStKamkh97Gvm3O5VKsXvvnqrpYezaux1PFTZ5oN+9JNdaMu0PzxOaZKWdNdnv4/9jR0GoSAMHHcG7x34u7mbk5cDn7sq7bGNjI+/2v1Z114QOzHMmT2NjIz3WU3Wz4xrn5Hd+Ud1LNHhcreIRrVNOOSXSUZamCP6IGO8fPgpCIlI1Ghsbsfd2VNU1oTmN+Q9iFvJHRJS3U5Tz5mgFoSKkUinq9u4uqIcRp7q9O0ilStu9r2q7IromNNh5LPXltF28/7QsYdfe7SUfjnvr3fQ1mEMOLO2stF17tzMnkitphZk6dWrcTRiVgpB8wMtvRTMx4Y296S/+mQeVdijs5bcmcUyeZaO6ZgIZQx5zSjzkMSfadleTqP4dtm5NXzed89HSBow5TI+szbWyjJOCUBEaGxt5473JVXVNqLFxVl5lo/yy2xe+pA+cV9ov6WPIv91R/uLqukn0dM2r9igIFalu786SD8fZu38AwA88rKTHTc+Oyy8I6Ut64ijkmkIxtw9o0VXJh4JQEaIbEtgDQNNH8wsY+Zul4RwZl0q+piDVTUGoCBoSkFpQ672UKHt66uWVjoKQiEx46unFZ8IFITNbCvwDMAm4wd2vjblJVU1/bUql0v+d6jChgpCZTQL+EWgBUsAmM2tz9+eiqlNf0u+rtr82C70ZsNY/P5EoTKggBCwEutz9RQAzuwVYDkQWhApRbV/SoL82M1Xj5ycSt4kWhOYAr2SkU8CiKCvUl3T10mcnEr0qWs++JEZaQvYDy8Ca2Uoz6zSzzp6enjI0S0RkYppoQSgFzM1INwLbhhdy9+vdvdndm2fMmFG2xomITDQTLQhtAprM7GgzmwKsANpibpOIyIQ1oa4JuXu/mV0MbCA9Rfsmd3825maJiExYEyoIAbj7emB93O0QEZGJNxwnIiIVREFIRERioyAkIiKxMfcP3CYjGcysB3ipjFU2AL1lrK/cdH7Vq5bPDXR+pfYRdx/zHhcFoQpjZp3u3hx3O6Ki86tetXxuoPOLi4bjREQkNgpCIiISGwWhynN93A2ImM6vetXyuYHOLxa6JiQiIrFRT0hERGKjIBQDM7vJzLab2a9H2W5mttrMusxsi5kdX+42jkce53eKme02s6fC6zvlbmOxzGyumd1vZs+b2bNmdtkIZar288vz/Kr58zvQzB43s6fD+V01QpkDzOzW8Pk9Zmbzyt/S4uR5fl8ys56Mz+8v4mjrEHfXq8wv4GTgeODXo2w/C2gn/fyjxcBjcbe5xOd3CnBX3O0s8txmA8eH94cCLwDH1srnl+f5VfPnZ8Ah4X098BiweFiZrwE/Ce9XALfG3e4Sn9+XgB/H3dbBl3pCMXD3h4CdOYosB272tEeBaWY2uzytG788zq9quftr7v5EeL8HeJ70E3szVe3nl+f5Va3wmbwVkvXhNfzC+HIgGd7fDpxuZiM9ELPi5Hl+FUVBqDKN9BjymvkiCP40DBm0m9kfx92YYoRhmk+R/mszU018fjnOD6r48zOzSWb2FLAd6HD3UT8/d+8HdgPTy9vK4uVxfgD/MwwV325mc0fYXjYKQpUpr8eQV7EnSC/p8SfAGuA/Ym5PwczsEODfgK+7+x+Gbx5hl6r6/MY4v6r+/Nx9v7svIP1k5YVm9slhRar688vj/O4E5rn7ccC9vN/ri4WCUGXK6zHk1crd/zA4ZODp5zvVm1lDzM3Km5nVk/6C/pm7/2KEIlX9+Y11ftX++Q1y913AA8DSYZuGPj8zmwx8iCocXh7t/Nx9h7u/F5L/DJxQ5qZlURCqTG3A+WGW1WJgt7u/FnejSsXMZg2OsZvZQtL/D3fE26r8hHbfCDzv7j8YpVjVfn75nF+Vf34zzGxaeD8VOAP4zbBibUAivP88cJ+HK/qVLp/zG3Z9chnp636xmXBPVq0EZvZz0jOMGswsBVxB+gIi7v4T0k9+PQvoAvYCX46npcXJ4/w+D1xoZv3AO8CKavklBz4NfBF4Joy7A/wNcBTUxOeXz/lV8+c3G0ia2STSwXOdu99lZlcDne7eRjoI/4uZdZHuAa2Ir7kFy+f8LjWzZUA/6fP7UmytRSsmiIhIjDQcJyIisVEQEhGR2CgIiYhIbBSEREQkNgpCIiISGwUhkQphZpeG1atfNbMfx90ekXLQfUIileNrQCvwWaB5vAczs8lh7TORiqWekEgFMLOfAH9E+m79wzPyP2JmG8NikxvN7Kgx8n9qZj8ws/uB75nZZzOeG/OkmR0ax/mJjEZBSKQCuPtfkl5f7lTgzYxNPyb9WIjjgJ8Bq8fIBzgGOMPd/wr4a+CisKDlSaRXOBCpGApCIpXtT4F/De//BfjMGPkAt7n7/vD+v4AfmNmlwDQNz0mlURASqS6jrbOVmf/2UKb7tcBfAFOBR83s4xG2TaRgCkIile1XvL+A5p8B/zlGfhYz+6i7P+Pu3wM6AQUhqSiaHSdS2S4FbjKzbwA9vL8i92j5w33dzE4F9gPPAe0Rt1ekIFpFW0REYqPhOBERiY2CkIiIxEZBSEREYqMgJCIisVEQEhGR2CgIiYhIbBSEREQkNgpCIiISm/8PwFzrlFuLGOwAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.boxplot(x=\"floors\", y=\"price\", data=df)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 5\n",
    "Use the function <code> regplot</code>  in the seaborn library  to  determine if the feature <code>sqft_above</code> is negatively or positively correlated with price."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0, 8086221.105223094)"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAELCAYAAABwLzlKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvX18XOd13/k9984bXgYESAJ8pylatCnZsR2JteRGq6iOY8tO1/J25dZMt1ZTZ6lN7SZp1q7tfvzSym0/0jab1EpdV9w4G2nTWFG4SaXtWnYUKzSbrWiblCNbEhGRpmgSfANAgMAM5v3es3/cewcDYAAMQAwGL+eLDz4z88xz73PvYHDPPef5PeeIqmIYhmEYrcBp9QEYhmEY6xczQoZhGEbLMCNkGIZhtAwzQoZhGEbLMCNkGIZhtAwzQoZhGEbLMCNkGIZhtAwzQoZhGEbLaKoREpF/KiKviMjLIvJ1EUmJyE0i8l0ROS0ifyQiibBvMnx9Jnx/T81+Phu2/7WIvK+m/d6w7YyIfKamfcFjGIZhGMuPNCtjgojsAP4SuFVV8yLyFPAN4APAn6jqkyLyH4GXVPWrIvKPgbep6v8iIh8B/gdV/XsicivwdeCdwHbgz4E3hcO8Bvw8MAB8Hzioqq+GYzU8xlznsXnzZt2zZ8+SfjaGYRhrnZMnTw6rau98/WJNPo4Y0CYiZaAduAy8G/jF8P3HgX8BfBW4L3wOcAT49yIiYfuTqloEXheRMwQGCeCMqp4FEJEngftE5NRCx9A5LPGePXs4ceLEIk/fMAxjfSIiP2mkX9PCcap6EfhN4DyB8RkDTgLXVbUSdhsAdoTPdwAXwm0rYf9Nte3TtpmtfdMixjAMwzBaQNOMkIj0EHgeNxGE0TqA99fpGnkhMst7S9U+1xhTEJFDInJCRE4MDQ3V2cQwDMNYCpopTHgP8LqqDqlqGfgT4G8C3SIShQF3ApfC5wPALoDw/Q3ASG37tG1max9exBhTUNXDqnpAVQ/09s4b0jQMwzAWSTON0HngThFpD+d2fg54FfgL4P6wzwPA0+HzZ8LXhO8/H87VPAN8JFS23QTsA75HIETYFyrhEsBHgGfCbRY6hmEYhtECmiZMUNXvisgR4EWgAvwAOAz8v8CTIvKvwravhZt8Dfi/QuHBCIFRQVVfCdVur4b7+biqegAi8gngW4AL/J6qvhLu69MLGcMwDMNoDU2TaK8VDhw4oKaOMwzDWBgiclJVD8zXr9kSbcMwDOMGOdo/yGPHznJhNMeunnYevHsv9+zva/VhLQmWtscwDGMFc7R/kC888wqDmQLdbXEGMwW+8MwrHO0fbPWhLQlmhAzDMFYwjx07S9wV2hMxRILHuCs8duxsqw9tSTAjZBiGsYK5MJqjLe5OaWuLuwyM5lp0REuLGSHDMIwVzK6edvJlb0pbvuyxs6e9RUe0tJgRMgzDWME8ePdeyp6SK1VQDR7LnvLg3XtbfWhLghkhwzCMFcw9+/t46INvoS+dYixfpi+d4qEPvmXNqONMom0YhrHCuWd/35oxOtMxT8gwDMNoGWaEDMMwjJZhRsgwDMNoGWaEDMMwjJZhRsgwDMNoGWaEDMMwjJZhRsgwDMNoGWaEDMMwjJZhRsgwDMNoGZYxwTAMo0HWcnG5VtE0T0hE3iwif1XzOy4ivy4iG0XkORE5HT72hP1FRB4VkTMi8kMRua1mXw+E/U+LyAM17beLyI/CbR4VEQnbFzyGYRjGXKz14nKtomlGSFX/WlXfoarvAG4HcsCfAp8Bvq2q+4Bvh68B3g/sC38PAV+FwKAAXwTuAN4JfDEyKmGfQzXb3Ru2L2gMwzCM+VjrxeVaxXLNCf0c8GNV/QlwH/B42P448KHw+X3AExpwHOgWkW3A+4DnVHVEVUeB54B7w/e6VPUFVVXgiWn7WsgYhmEYc7LWi8u1iuUyQh8Bvh4+36KqlwHCxyigugO4ULPNQNg2V/tAnfbFjGEYhjEna724XKtouhESkQTwQeCP5+tap00X0b6YMaZ2EjkkIidE5MTQ0NA8uzQMYz2w0ovLHe0f5ODh49z1yPMcPHx81cxVLYcn9H7gRVW9Gr6+GoXAwsfokxoAdtVstxO4NE/7zjrtixljCqp6WFUPqOqB3t7eBZyqYRhrlZVcXG41iyaWQ6J9kMlQHMAzwAPAw+Hj0zXtnxCRJwlECGOqellEvgX8mxoxwnuBz6rqiIhkRORO4LvAR4HfWcwYS37GhmGsSVZqcbla0QRAeyJGrlThsWNnV+Tx1tJUIyQi7cDPAw/WND8MPCUiHwPOAx8O278BfAA4Q6Ck+yWA0Nh8Cfh+2O8hVR0Jn/8K8PtAG/Bs+LvgMQzDMFYzF0ZzdLfFp7StFtFEU42QquaATdParhGo5ab3VeDjs+zn94Dfq9N+AnhrnfYFj2EYhrFa2dXTzmCmUPWEYPWIJixtj2EYxipnpYsm5sKMkGEYxipnJYsm5sNyxxmGYawBVqpoYj7MEzIMwzBahhkhwzAMo2WYETIMwzBahhkhwzAMo2WYETIMwzBahhkhwzAMo2WYETIMwzBahhkhwzAMo2WYETIMwzBahhkhwzAMo2WYETIMwzBahhkhwzAMo2WYETIMwzBahhkhwzAMo2WYETIMwzBaRlONkIh0i8gREekXkVMi8i4R2Sgiz4nI6fCxJ+wrIvKoiJwRkR+KyG01+3kg7H9aRB6oab9dRH4UbvOoiEjYvuAxDMMwjOWn2Z7Ql4Fvqup+4O3AKeAzwLdVdR/w7fA1wPuBfeHvIeCrEBgU4IvAHcA7gS9GRiXsc6hmu3vD9gWNYRiGYbSGphkhEekC7ga+BqCqJVW9DtwHPB52exz4UPj8PuAJDTgOdIvINuB9wHOqOqKqo8BzwL3he12q+oKqKvDEtH0tZAzDMAyjBTTTE9oLDAH/p4j8QER+V0Q6gC2qehkgfIzq0e4ALtRsPxC2zdU+UKedRYxhGIZhtIBmGqEYcBvwVVX9aWCCybBYPaROmy6ifS4a2kZEDonICRE5MTQ0NM8uDcMwjMXSTCM0AAyo6nfD10cIjNLVKAQWPg7W9N9Vs/1O4NI87TvrtLOIMaagqodV9YCqHujt7W34hA3DMIyF0TQjpKpXgAsi8uaw6eeAV4FngEjh9gDwdPj8GeCjoYLtTmAsDKV9C3iviPSEgoT3At8K38uIyJ2hKu6j0/a1kDEMwzCMFhBr8v7/CfCfRCQBnAV+icDwPSUiHwPOAx8O+34D+ABwBsiFfVHVERH5EvD9sN9DqjoSPv8V4PeBNuDZ8Bfg4YWMYRiGYbQGCYRlxmwcOHBAT5w40erDMAzDWFWIyElVPTBfP8uYYBiGYbQMM0KGYRhGyzAjZBiGYbSMZgsTDMMwADjaP8hjx85yYTTHrp52Hrx7L/fs75t/Q2NNY56QYRhN52j/IF945hUGMwW62+IMZgp84ZlXONo/OP/GxprGjJBhGE3nsWNnibtCeyKGSPAYd4XHjp1t9aEZLcaMkGEYTefCaI62uDulrS3uMjCaa9ERGSsFM0KGYTSdXT3t5MvelLZ82WNnT3uLjshYKZgRMgyj6Tx4917KnpIrVVANHsue8uDde1t9aEaLMXWcYRhN5579fTxEMDc0MJpj5yLVcaawW3uYETKMNcBquDjfs7/vho4pUtjFXZmisHso3LexOrFwnGGsctaL/NkUdmsTM0KGscpZLxdnU9itTcwIGcYqZ71cnE1htzYxI2QYq5z1cnE2hd3axIyQYaxyVtLF+Wj/IAcPH+euR57n4OHjSzovdc/+Ph764FvoS6cYy5fpS6d46INvMVHCKseK2s2DFbUzVgOROu5G5M9LcQyReq0t7pIve5Q9NUOxTmm0qJ1JtA1jDXCj8ueloFYgAdCeiJErVXjs2NmWH5uxcmlqOE5EzonIj0Tkr0TkRNi2UUSeE5HT4WNP2C4i8qiInBGRH4rIbTX7eSDsf1pEHqhpvz3c/5lwW1nsGIZh3BjrRSBhLC3LMSf0t1T1HTVu2WeAb6vqPuDb4WuA9wP7wt9DwFchMCjAF4E7gHcCX4yMStjnUM129y5mDMMwbpz1IpAwlpZWCBPuAx4Pnz8OfKim/QkNOA50i8g24H3Ac6o6oqqjwHPAveF7Xar6ggYTW09M29dCxjAM4wZZSQIJY/XQbCOkwJ+JyEkRORS2bVHVywDhYxQs3gFcqNl2IGybq32gTvtixpiCiBwSkRMicmJoaGgBp2sY6xdTrxmLodnChJ9R1Usi0gc8JyL9c/SVOm26iPa5aGgbVT0MHIZAHTfPPg3DCFkJAgljddFUT0hVL4WPg8CfEszpXI1CYOFjtJBgANhVs/lO4NI87TvrtLOIMQzDMIwW0DQjJCIdIpKOngPvBV4GngEihdsDwNPh82eAj4YKtjuBsTCU9i3gvSLSEwoS3gt8K3wvIyJ3hqq4j07b10LGMAzDMFpAM8NxW4A/DVXTMeAPVfWbIvJ94CkR+RhwHvhw2P8bwAeAM0AO+CUAVR0RkS8B3w/7PaSqI+HzXwF+H2gDng1/AR5eyBiGYRhGa7CMCfNgGRMMwzAWjmVMMIw1yGooXmcYC8ESmBrGKmG9FK8z1hdmhAxjlbBeitcZ6wszQoaxSrDcbMZaxIyQYawSLDebsRYxYYJhrBIevHsvX3jmFXKlypR6Pas5N1sktHjt6jhlT0nEHPb1pU1wsY4wT8gwVglrLTdbJLR4fTjLeKFCvuwxlitz7lrWBBfrCPOEDGMVsZZys0VCi2vZCg6C4wi+KuP5Cls3xKwY3jrBPCHDMFpCJLQoeT4SphYWgZLnm+BiHWFGyDCMlhAJLRKuQ5S4RRUSrmOCi3VEw0ZIRN4gIu8Jn7dFyUkNwzAWQ1QEL52K4aNUfB/fV7raYqtecGE0TkNzQiLyPxOUw94IvJGgBMJ/BH6ueYdmGMZa5p79fTxEMDdU8cYpheq4PZs6TR23jmhUmPBxglpA3wVQ1dNhoTrDMIxFs5aEFsbiaNQIFVW1FJZlQERizF/F1DCMVYglSTWWk0bnhL4jIv8caBORnwf+GPh/mndYhmG0AkuSaiw3jRqhzwBDwI+ABwmKw32uWQdlGEZrsCSpxnLTaDiuDfg9Vf0/AETEDdtMyG8Ya4gLozm62+JT2mzNjtFMGvWEvk1gdCLagD9f+sMxDKOVWJJUY7lp1AilVDUbvQifN/StFBFXRH4gIv8lfH2TiHxXRE6LyB+JSCJsT4avz4Tv76nZx2fD9r8WkffVtN8btp0Rkc/UtC94DMMwJtfu5EoVVIPH2jU7R/sHOXj4OHc98jwHDx+3uSLjhmnUCE2IyG3RCxG5Hcg3uO2vAadqXj8C/Laq7gNGgY+F7R8DRlX1ZuC3w36IyK3AR4C3APcC/yE0bC7wFeD9wK3AwbDvgscwDCNgriSpJlowmkGjc0K/DvyxiFwKX28D/t58G4nITuAXgH8N/IYEGu93A78Ydnkc+BfAV4H7wucAR4B/H/a/D3hSVYvA6yJyhmDNEsAZVT0bjvUkcJ+InFroGKpqcnPDCJlt7U6taAGgPREjV6pYolHjhmjICKnq90VkP/BmQIB+VS03sOm/A/4ZEKX42QRcV9VK+HoA2BE+3wFcCMeriMhY2H8HcLxmn7XbXJjWfscixxiuPWgROUSQIYLdu3c3cJqGsfYx0YLRDOY0QiLyblV9XkT+zrS39okIqvonc2z7t4FBVT0pIvdEzXW66jzvzdZeL5Q4V//5xp9sUD0MHAY4cOCAeUmGQSBaOHcty3i+QsnzSbgOXW0x9mzqbPWhGauY+TyhnwWeB/77Ou8pMKsRAn4G+KCIfABIAV0EnlG3iMRCT2UnEIX4BoBdwECYkWEDMFLTHlG7Tb324UWMYRjGPLxr70a+d24ER8AJSy4MZkoc/BsbW31oxipmTiOkql8UEQd4VlWfWsiOVfWzwGcBQk/ok6r690Xkj4H7gSeBB4Cnw02eCV+/EL7/vKqqiDwD/KGI/BawHdgHfI/Aq9knIjcBFwnEC78YbvMXCxljIedlGDfCUqbEWe70Oi+cHaG3M0GmMOkJpVMxXjg7wq82bVRjrTOvOk5VfeATSzjmpwlECmcI5mO+FrZ/DdgUtv8GQZYGVPUV4CngVeCbwMdV1Qu9nE8A3yJQ3z0V9l3wGIaxHCyluqwVSrULozk2dybZ29vJ/q1d7O3tZHNn0uaEjBtCGnEEROTzBJLsPwImonZVXfOhrAMHDuiJEydafRjGEtDqxJwHDx9nMFOoqssAcqUKfekUXz90Z8v2tZLHNFYvInJSVQ/M16/RdUL/CPjHwHeAEzW/hrEqWAlrXKJy1rUsVl22lPtqlPkWshrGYmh0ndCtBEboLgJBwn8lKGpnGKuCpV7jEnlVr10dpxwWY9vXl57Tu9rV0z7Dk1hMSpyj/YOM58tcHsuTirn0ppOkU/EbTq8zn6dYW4RuYDTHTivzYCwBjYbjngLGgf8UNh0EulX17zbx2FYEFo5bG9z1yPN0t8WJamJlCmUGxwsUPeWdezYu6GIaeVWlise1iVLQqLA5nSDuutUMA7NtF3eFtrhLvuxR9nTW/ks5dqP7vZFjM4xaGg3HNeoJvVlV317z+i9E5KXFHZphLD+1XkimUObS9QKKkoo5DGYKfPLIS/R2JskUK/POF0Ve1bVsBQfBcQRflfF8ha0bYrN6V0vhSURjb2hLkYy5DGeLFCs+E0WPRz/ytkUbDMuGYLSKRo3QD0TkTlU9DiAidwD/X/MOyzCWlgfv3ssXnnmFXKnC4HhggAA2dyapeMr1XJlsocLNfZ3V+aKHoO4FOMocUPJ83Gq14WDdzHzzMjdazro2a0FXW5yutjiqyli+vGT7jbBsCMZy0Kgw4Q7gv4nIORE5R7DO5mdF5Eci8sOmHZ1hLBG1iTmLnpJwHbZvaKOrLc5wtogj4Kk2VMgtKneQcB2iaLYqJFyn6WUPmlVqwUo4GK2iUSN0L3ATQQaFnw2ffwD429TPpmAYK4579vfx9UN38s49G9m6IUVXeOdf8nwgMCIRc3kBkUosnYrho1R8H99XutpiTVeLNUuhZso3o1U0msD0J80+EMNoFtNVX+/au5EjL14kV6rQFndxHaHiKb3pZHWbyAuYTTEWze1UvHFKoTpuz6bOpqvFmqFQi85xolhuWOlnGEtFQ+q49Yyp41YmjS48nU31df9tO3jh7AgDozk6Ei7XJkp0tcVn9Dny4sU1rRgzVZzRLJZaHWcYK4baC2d3W5zXh7M8+AcnSadiM+7gZ1N9vXB2ZMoq/8io1XoX60Exth7O0VjZmBEyWsKNpNCpvXCO58vV9TK5YmWGsq1R1Vc91drnnn656YqxVqcSaqYqrtXn1mrW+/k3SqPCBMNYMm40hU5typrhbBEHwXWEsq8zlG03ovpqlmLsaP8gBw8f5/Yv/RkP/sFJzl3LtiyVUDPPsdVpklrJej//hWBGyFh2aj2ZRiTR06m9cJY8H5FJiTRMvZO/EdVXMxRjtRenQtnHV+Vatky2WFnw57AUNEsVd6N/49XOej//hWDhOGPZWUgIqF5Io3bhacJ1AkOEVNVttXfy86nJ5gqZNEOJVntxiha7KjCUKZJOxef9HBaSq26+zzEKQTYjH9x6X/y63s9/IZgRMm6Yhca+G03kOV2AUJ3v+eBbeOiDb+GxY2cZy5Wo+MrGjjidyVjdO/nZshTMun+YYoiWMo5fe3FKuA4VTxFncq3SXJ9DqeIxXqgE/Uoe565l58zs0Oh5LvU5wtIla12trPfzXwgWjjNuiMXEvhsNAc0V0ogWnp74/Ht57H+6nT2bOhnLl+lLpxqWF7ciZFIbStzcmcRHKVd8Kp7y6uUxBkbzvGvv1HLZ0XFmCkGuupjj4DjCeL7S0PG24jzX++LX9X7+C8E8IeOGaETiW89TijyZKAT0rr0beezYWT739MvVPjeibGuEaP+ZQpmhTJGS5xN3hLF8ed5tF6t8qg0lplMxOvIu1/MVXIFUzCWdinHkxYu8bWd3dX8LzVU3/dheuzrOtg1t1ffH82WGs0XOXctx8PDxpqi21nvZh/V+/guhaUZIRFLAMSAZjnNEVb8oIjcBTwIbgReBf6CqJRFJAk8AtwPXgL+nqufCfX0W+BjgAb+qqt8K2+8Fvgy4wO+q6sNh+4LHMBbHfIZirpBatE5ntj7pZIx82WtaSGNXTzvnrmW5li0jQlVhlylUONo/OOsFo5Ew3mxMvzj5Cts2JNncmar2mW7Eo9BONXwns+eqq3ds2aLHcLZIbzrFeL7MpbE8AElXFnTsC6UZYb7VxHo//0ZpZjiuCLw7LAHxDuBeEbkTeAT4bVXdB4wSGBfCx1FVvRn47bAfInIr8BHgLQQ57P6DiLgi4gJfAd5PUHTvYNiXhY5hLJ75JL6NhIJm66OqTQtpHO0fZGA0x+WxwAMqhiExgJ72+Kyhqkf//DU+9sQJzo/keH14gsHxwqzhrUiKfdcjz3Pw8PFqiDIKJf7XT7+brrY4mzqSU7ab7t3Ml6vuXXs3Vsf51Sd/QKniTfksN3bEGc2VyZUqDGeLAAhCX1fKVFtGy2maEdKAbPgyHv4q8G7gSNj+OPCh8Pl94WvC939Oggpk9wFPqmpRVV8HzgDvDH/PqOpZVS0ReD73hdssdAxjkcwX+26kDPVsfSZKXjXz9ULne+biaP8gnzryEpfHClPaK77S0xZnc2eyrorp0T9/jS8/fwbPD4yVrzCYLXF1LD/jnBqdK2tknU6UAfymzZ1sSMVoi7tsaI+zZ1NnNbVQNM5EqcK1iRLjNSHFTR1J0qkYfekUhUoQctzenSKdCjzYVqi2ZjPQxvqjqXNCobdyEriZwGv5MXBdVSthlwFgR/h8B3ABQFUrIjIGbArbj9fstnabC9Pa7wi3WegYwzd8suuU+WLfjaiE5urTSEhjtvmZ2dofO3aWTKGCGxajUw3ujkRgouTNGvL73b98HUdAJTBAEYPZEo4j3LS5s9o221zZI9/snzOZar7sMZ4vE3eEux55foacevo5f+XojxFg64YUIkIq5lLyfIazxWqW8HzZY19fmq8fupODh4+3XLV1I+FMY+3RVCOkqh7wDhHpBv4UuKVet/Cxnkeic7TX8+Lm6j/XGFMQkUPAIYDdu3fX2cSoZS5DUTsRX5sgszak1kif2YguaGXPYyxX5vJYnhfPj/KBt27h5Pmxuhe6C6M5Kr5PzHWIOQ7lUB6tCoXKzLGjC/54oYJQ/0s0lC3xi+/cWO3/4vlRPN8nGXPpTSdJp+JUPJ9z13Ls2dSOK/CD86N879w1tqWTJBIJxvJlOpMxFCj7OuW47x+4zgtnR7gwmqOzJuGqr4oAl64X6Gn38PwghFn2PMbzJWKuw1i+TMJ1uOuR50knY1XhxUI/68VQ70bA8tUZtSyLRFtVrwNHgTuBbhGJjN9O4FL4fADYBRC+vwEYqW2fts1s7cOLGGP68R5W1QOqeqC3t3dxJ20AU4vJzRZSa6QP1A/hPHbsLGXP41q2jKcQdx18Vf7zS5cpe17duahdPe3EnKAgnesIcdepGpeORGzK2LVhNUeCOxa/znn2pRO8cHak2l8EHBEqvnLpeoFMoczV8SJxV/B85fJYEQVcEYZzZSZKHl+67610tyfY0Bafctxlz+MrR39cDbmdG8kxmivj+UFxPiHw6AYzwT5jDjgCA9cLFMseQqCkixR2AlUV4FKFOOsxW0jytavj84ZojfVDM9VxvUBZVa+LSBvwHgIhwF8A9xPM4TwAPB1u8kz4+oXw/edVVUXkGeAPReS3gO3APuB7BNeMfaES7iKBeOEXw20WNEazPgMjoJGQ2nx9Zgvh5EoVcsVK9aIP4AqUVRnLlaeozqIL3ZfueyufOvISo7kyKsGf33WE7vY4v3n/26ccR+1d++aOBIPZUvU9ITBKW9JJNnUE80hR/y3pFJfG8gS7V66MFSj7Pju72xjKFKvHq4Dna9VA1lMbjoUGJ/IcPF9xJMiysLkzyaWxPJ4fFixXEBF2dbcRc4WhTJHedHKK1wHQ05Hkm//0TprJbB5P2dOmqh6N1UUzPaFtwF+E5b+/Dzynqv8F+DTwGyJyhmA+5mth/68Bm8L23wA+A6CqrwBPAa8C3wQ+rqpeOOfzCeBbwCngqbAvCx3DWPnMpqArVXyKYf64CNXAEyh6U32W2nmmf3v/29nX14mIICLc3NvBR+98A48dOzvF06oVTWzZ0EZfZ6K6P5HAAPV1par7jvp3tcXZviEwBEpgrPb1dhKL0gzJ5LEmXKdqIOsJFYqeTzI2+a8a5cgreX51HAiMYsyVatnySNzRKq9jNsFJIubYQk6jStM8IVX9IfDTddrPEijbprcXgA/Psq9/DfzrOu3fAL6xFGMYK5sLozlcgbNDWUqeT8J12NyZIOEKRcfB8xU3XD/jo3SlYuTL/qzzTPUm+ut5Wp0Jd8pd+5YNbTiOMJors7Onjba4O+Ui+tixs9WJ/662OF1tcXKlCn3pVHXuS4BCedJAljw4M5Rlz8b2uvNjMcchnZr8V93cmeTi9TwxJ5Cxx1whEXPY2BGf4vnlyx4d044/al8Or2M2wUmU884WchpgGROMVUJnwuXM0ASuCK4E5bgvXi8EHsy7tvGVoz+m4ivJmEPcdciXfRyBc9dyqCqdyRi/fNdNs17oZgsdiQhj+RIXr+cpe4qGarod3SkS4aR/R8Il4Tp87umX55z4v2d/H/cPXOd3/uLMlLEVKFf8al2k6dkk7nv79ikKupgbhA57O5OM5ct1+0Tj/vJdN9Vtb1T0cSP1cOYSnKy1hZxWOwhUlULZpxiKexrFjJCxKqgu56qVp4XzH7/6njfxtp3dPHbsLKevjpMpBh5AplDBEUCE9oQ7Ix1OLbNlfrgylkdE8H2trg9ynWBN0UTJ48O376yWAO9ui5MPhQDRxP/OaZLx3/3L11GlmvVAwnOKuw5dbcEi2a8funPGMUbnFxmmz//CrfP2icadrX0ulkJGvV5S16xXyXnF8ylUfAplj2LFp1TxiabYy149+U59zAgZq4JMscKO7hTD2VI1HLe1K0m2GCwzAgAdAAAgAElEQVQHi+6so3UwV8YKOI7ghAYkU6iwdUNshgw4uoMdyhQZzhTZumFyEWe+7FHylN50glzJQ0Sn7C+dgi8/fwZQUjVSbJg58R9dqCZKFWKOUAlXscUcwQ0Vc/UWvC7k7no272IxXsdSyajXmsdTj/UgOVdVihWfYujpFMo+FT8wNL4q50dynLqc4dTlcV69PM654YmG921GyGgKSx2eiOYX9vZOLgiN5lpqmZLs06mf7DM6ttODGTKFCj3tcbZ2Jbl4vcDAaJ4d3UrMdap1e9ri7oz9FSo+pYkSnq8kY5NS7O3d0JmMzZj4jy5UqZhLxVccUXwFTxVHZUYeuFbfXVs9nMZZi5+V52vVw4keIy9nNFfi1OVxTl3O0H95nP4rGSZK3jx7nB0zQsaScyMX0NmMVzS/MJwtMJYrU/R8Yo7DfW/fPmX7Kck+/ZnJPjsSbvXYcsUKFc/n6nggmY67Dr7vc2W8yE2b2km4DlfHC2QL2aD4XE0YTcMXjgAqOI7gowxliriOzJj4Pz2YIVesUPK0KrGGIPOC7ytdHfEpeeBePD86JRNCo3fXS2X8rR5O46yFz6pYqTE4Zb8aTitVfM4MZjl1ZZxXLwUGZ3q6q4iuVIxbtnVxy7Y0b92xgQ83mJnTjJCx5Cw2PFHPeH3qyEts6kiQLXng+wxnS6gGZQ/irvCVoz/mieM/qSquImPV1RZjOFPCD9cBpVPBRT7hOtUFo7kahVpVRy1CwoVc2SfuCl0pl6HsZB42V8CJVq0qbO5IMJovBytYRSlU/LoZFzKFChXfr6YIiuZtXUeqeeCiFD5xV6ZkQtjezZxVV+f6/BZr/OulEzIZdX1uJONHK/B9rRqcQiUwOkH6KuXSWIH+y+O8GobWzgxmqfgzRQYxR3hjXye3bE1z6/Yubtnaxfbu4IYp5ji0Jdw6I9fHjJDREAu5w15snZ7pxsvzldFcmUyxws29nZwZyqIKO3vaUKVakiBXrPD6cJaPPf59/PAi7wik4i4xR0i4QV63d+3dGKjoPJ/p/1eRDQLIlZU+N1DgXc9XiElgNKJsCe0xQeIuHUmXzZ0p2hIxhrNFihWtZlwAOHj4OBdGc2EuOChWmMGvvftmfvU9b6r2j84/KtuATJb+nu/ueimN/5EXL3L/bTt44ezIDYkK1oNqbKULMMqePyW0VqoEN1/ZYoX+MKx26krwONv/6NauFLdsS3PLti5u3dbFzX2dJGIOriMkY8Har2T4G3MXtvzUjJAxLwu9w15snZ7psfWhTBFHAmMkIpTDUNZPruWqGQdirgQKnUxximHxNSiB/aF3bOPKeInTgxlePD+Kqs4wQBFlz6+uvWmLu7w+NoGD4MQcYgRj797YPmXNT1ScLuYGxxcZoNrP68pYoe7dJMBXv/Njnjo5wK6edk4PZuhIuJwdypIve9XjLAKnr2ZIp2J8/hdurbufep8fBAqmF8+PzkiGWstsxuuFsyPVmk+LodXzWsvJShFg1AoIChWPQjnIJ+j5yuvDE7x6eZxTl8fpv5zh/EhuZuJMoD3h8uataW4NQ2v7t3axsSOBIzJpbOJudTnEjWJGyJiXR77Zz+B4AU+DcFZvOjmlzPZ0Hrx7Lw/+wUkUxUHQMOoV1emZ7Z91emy95PlVo/Hq5fGqRBrCuRQNUtXMZlQU+M8vXWbv5g4y+TLFBtYupFMxxgoV+q9kqoZDvEnv6vxIjoHRwAObvp4nusDXejQQhC5mM0L5ss/V8QLD2SKlis/oBMRcB3dapm6kTqbdaUz//MbzZS5eLxBz5jYCzZpYXw+qsVpa4fVVwlpYQWhtUiY9lCmG4oEgtHb6aoZCZaZs2hHYs7mDW7Z2ceu2NPu3dbF7Yzsx15ni3UTeTjMwI2TMydH+QfqvZIDgIljxPS6M5NjZ0zbrReqe/X10Jl0KZb8mu0GKdGqmaqyW6bF1Aco+QTis5oocGQRfp7bXQzVYsDpfvyDppzBWqNCVchnPT8bOoi19BUcDr2x6ddhapl/U50tP6EpgqKNDjOnknBFAKuawry89b9n06fM4VzPBBPJ84oZmTayvRdXYbCyX1xfJo4theK3s+eTLHq9dzVTVaq9eHme4JsdhLRs7EkFYbeukl9OWCAxMIjQ2keFZrlJrZoSMOXn42VNT7sBVgxrrl67nObBn06zbvWlL14wLW65UmfPCNj22HnMdKp6H6zpUKn41Yej0/w1nutcwDX8eAwRQ8aE94bAxFaM3naK/OI43zXMSwPNhZ08K15n0BKcbg+mpfuZbtlecdoda9oM5K4FARBEasfnKpk+fx4kyO0Rrl6bvI6JZE+trQTXWKM3w+jxfJ41OKCCo+D4XatbknLqc4exwtu73PxFzeFNfZ6hYC4xOXzpJompoXJJxh4TrBGKbFmFGyJiT16/liDnBxTciqHdD3YtUdEF+5dIYmWIFUWhLuKRTMRIxd94LWxRbP9o/yIN/cBJxZHL1tUBcBAT2b+3ilYtj9Yv71CEyYLOxuSNOKhHDFei/Mj5r2hHXEdKpOKrK6avj3Pvb3+H0UJa449CVinHyJyOUwm0TrrBtQ6paNnw2pr8bd5zquiQRIR5eIGYrmw7153Gihbu11DMCzZpYX22qsRthKby+UmVyIWih7FH2fMZy5VA0EITV+q+MM1GsvyZnZ09bYGxCxdrezR2k4jGS8akhNbeFBqceZoSMeXFEcFypSowhSEtTr95PVGAuX/KqHkqu5FHxlY/fs7uhC1u1Jg+BJFrEwQslpOXwlu9HF8cA6Ig77O1N89qV8bpzPo3+u6XiLteyRfLluf2WKJv1cLZIpuiRLeVwJRBeDE4LgZQ85Scj+Tn3V9c4SuDdRfNIW7uSdcumz3fRW4gRaMbE+kpXjS0lC/X6IgFBoTzp6RTKHj8eyvLqpdDLuTLOpeuzr8nZHxqcW7Z1sX9rmp6OxIyQ2kKVahHLOb9lRsiYk72bOzg9mMV1AvdeNZgMjcecGYqr6O78WraC4wgxcfAlyPK8dUOKb/zocrU66Fxf7Gg/Wzekgn9CAfWV6TbGEZgo+VwYmSAVdyl6dTTQUFWuzcXF6wXi7twmS4HNnQlypQrXJkoIUKwsbTkqRwIhgxfOhe3r6yRbnMzCvZCy6SvBCDRq3Fa7lHs+gz89z1qx7HF5LM+pyxlevTxO/+VxTg9m635PXUd4Y28oHtgehNV29rSTChVqS6lUg+VXNYrVdJubAwcO6IkTJ1p9GC3jaP8gnzryUnWxZVQttC+dZHNnknzZYzxfZlNHgtNDEyRdoegFRdqEQO7sqbJ9Q4qB63n2bOqY8k8azWHUXnw+9/TLdLfFEREGxwtczRTrHpsbVTqd5yvcHg+yJ5TmMUSuI3MKGJKu0NeVoiPh8tpgtqp6a2DKaQrVuS2C+a1oezfM2uA6QqHi0ZGI8ehHfrruvFMkQihVPDKFCsVKEL77+D1vrK47Wi3UXvRqvxvNqvjaLKp/o5EJdnS389F3vYG/cdNGCmWfsXyJv76SCSXSgadzfZY1OX3pZFUefcu2Lt60JU26LU7CdaaE1ppFFMadPp/bl04tSLIvIidV9cC8/cwIzc16N0Iw+c81MJpjLF+mPeHSG+ZsyxTKDIzmiblBiYWS51PxtFo22/cDT8gLUwXs25Ku7ncoE0iTozxtflQdFIg70NOe4Hq+QrkmDNgoUXodCIyLI+B5OqdIYL55o97OBP/2/rfz2LGz/ODCaKhomwwRLhZXpqrhIo9sU0eCRMzl/tt2VDMp1F6kb9+9gW+8fJWK75N0HTa0B1khejuTDIWS77grvGlL14r2LJbqotcqpudZy5U8Xh/OVoUDpy6P85Nr9dfkpOIO+8OQWqRY27qhLTQ2y69UA7jrkeerN4ERqspYvsx//fS7G95Po0bIwnHGvNSGVKIvaETtgtItG4KS1k5Y6kAkuOSnU3GGsiV2dk9NNjo6USIQhukMD6Tsw2C2VC1Ut1Bqt/F9RaMSEHPsa65htqSTdKZiVW9kSzrJpbFCWFL8xoyQp4HhKXtBqp5A0p6sFsT73b98fUaJ7lypwrf7h9jZ0zZlXdBQJs94vhycv0C+DK8PZ2cNp6yEMNhqk3JPz7MWhdWidTl/fTUzpWhhhBCtyUlX1Wo396VpT7grRqkGy69qNCNkLIh6C0qjC2dXeCEZzhbJl4Oy0lHKnJ72oARDLaXwouvXsTK1oa75JNgOQT632RaEBvV/FnO2ASIwkiuRKZQZy5fZ15dmMFNg+4Y2hrNFKjp3GG+6p1OPhOug6iEiUzKFRyW6d9cpkz29fTgb3BBE2b8bKWOxEjIarGQpd22etWLF53ouCKtFaW5OXR5ncJZwcU97vGpsbtnaxVu2b6CnIxHO4wSezkpTqsHyqxqbZoREZBfwBLCVYKnEYVX9sohsBP4I2AOcA/6uqo5K4Pt9GfgAkAP+oaq+GO7rAeBz4a7/lao+HrbfDvw+0EZQ5vvXVFUXM4bRGNO/oG44R7S5MwkEF2zPV+Kuw1u2b5hS0G36F7t20el0JgvICTu622ZNMRL18WZxlxLhHIunHtqAQYMpNfOC/YtMST1Um2j0ps0d5Mse50dyRGm2I0+mWPEYzBTxdGp4sB75skfCdWbI+eYq0T29PbohgMm1VNPLWNSyUjIarCQpd22etVypwtmhiSkJPc8OT9S94Yi7wr6+NLduT1cFBG/Y1B5IpG9QqbbcLLegpWlzQiKyDdimqi+KSBo4CXwI+IfAiKo+LCKfAXpU9dMi8gHgnxAYiDuAL6vqHaFBOQEcILgunARuD43K94BfA44TGKFHVfVZEfnfFjLGXOexVueEbiQMUztH1JmMMZQtsqEtTsXzuRhKSnd0p4iF5a/b4m5wMfZ9Eq5Le9JlX1+aVy6NMV6or2iLiLvC7o3tjOdLDGfLKIFMuqc9TsVXRidKhHk+gamBsVRNLN3zfUqeEnOEHd2pOaXTSTcwOL5OigVUoaKKg+I4Dr2dCYYyRYpeUC5iU3scRNjQFp9yIfU9jyuZQLo9m6dWe64xASWQw8cch3Qqxj+48w1zzgl5YVnzUsVDkSD/HVQ9oUidOH2OZali/0tB7XdquVR80/OsDY4XefnS9aqH038lqDdVjx3dbdWMA2/dESwG7UzGSMaDUu/NSnGzmmj5nJCqXgYuh88zInIK2AHcB9wTdnscOAp8Omx/QgOreFxEukNDdg/wnKqOAIjIc8C9InIU6FLVF8L2JwiM3LMLHSM81nVDbRjGFfjB+VE+9sT32dfbyWfef8uC69VAcNcU1MAJLtLRfEmp4qMECzedULjQjsu79m7ku69fm/dYu5IufekUY7kSiZgTLujzGcuX+ZWffSNPnbjAwPVCVW3mMJmhoFDxcYBk3EVEcCUQPpyfZ+1O7XqjmBusUXJEUF8RR6h4PpeuF6rjqMLwRJmEK+SKHkXPpyPh8st33cRTJwfoaY8xnJ07gzgExjVX8nAlSA0U5Yt7287uGSW6I29sY0ecsVyZQsULEkyGYclICeg60JNK1vUsVlIYbDkSgEYCgkLZI1uscOryOK9cGq8KCC5er/+96EzGJlPdbO/ip3Z0saWrrapUS7jLKxxYayzLnJCI7AF+GvgusCW66KvqZRGJvnk7gAs1mw2EbXO1D9RpZxFjrCsjFIVhPF+5PBYUdHNFODeSm3NOYLY5hPtvCz72sheo2Fz1iblOUI0x3NZ1gjtD8ZVsscJXv/PjecNiMQeyJZ8H797Lp468VC0G52uQ+PPR50/TFnd5w8Y2hrOlKZmnI3ygUPYCT0lgSzrBeL5CyfPnXTsEVKXPXihPFxEUneJxRaG24MLvc8vWNPmyx5EXL9KZcLk6NrfRCw+NiWKwwDcRc6vzQlF47OuH7pzyN4mSpG5oS7G5MxB8DGcLDGdLuI7ganCn7/uByvDT9+6f8Tdd7jDYUnjfC9l2MvtAhfPX8vzo4vWqRPq1q5lZ//5B6Cyoi3Pov7uJ//HATtpqwmpmcJaWphshEekE/m/g11V1fI4/YL03dBHtcx5OI9uIyCHgEMDu3bvn2eXqI1IjvT48US2JoERzObNnx46MV8VTXh+bqM5BPPr8aXrTyaoBqChUZlECRHNGkRJsrj9YxQdBeezYWUbCsNv09zNFD9cps7e3k5cvjdXdj0bjKIxOlCn7StwRyg2q2qI5AF81UP1NO+7aiHbJU84MZqvZHRai3o7WPNUKOGZTidVTlI3lAlXclnQqyMod/n3OXZvgc0+/zK5jUy/eyxn7v9GCe/Nt6/taLdB2baLIDwfGePXSeLV0wWhu9jU5+7elOXM1S8X36W5PBJJ+hHy5wl+eucbH371viT8No5amGiERiRMYoP+kqn8SNl+NQmBhuG0wbB8AdtVsvhO4FLbfM639aNi+s07/xYwxBVU9DByGYE6o4RNeJURhmCg/GUyWwK530YvuQr93boSYBEYmJsG6oMjbuTxWXyFUS2B8Ji+wjYibfVVO/mRkTnXZ9XyFVKbQkJQ7SmffgGJ7Vubbpl7K/LmoPY7oZmA8X6arbfZCdvVCaVHJ80tjeZwwllfyFTylWK7wgwujfOyJE7ypr7PqGS1XHZwbEUHU23aiWOYrR8+wf1sXr10d56ULY2Gqmwznhifqr8mJObw5lEe/ZXsX79jVw86eIKz287/1HfrSqSleTnti7qzvxtLQTHWcAF8DTqnqb9W89QzwAPBw+Ph0TfsnRORJAtHAWGhEvgX8GxHpCfu9F/isqo6ISEZE7iQI830U+J3FjLHU577SicIwbjhxLYCPsrkzNeOiV3sXmnSFfDkwOm44x9PoRTxfnpl0sZFLta/MyGZdixNMnZAr1U/qOBvziQSWk+lH4iv8ZCRHd1uMrrZE3fBYtUpszULVmOPg+z6OBGtNyjWZx4ez5UApKPD68MSyS7FnWwt0+uo47/93xzg7PAHATZvaZ8xLvnZ1nFyxQtkPhCXtyRiqMDCa52/95tG63y0Bdm9qD6XRXbx91wb2b+uiIxGrq1TbvbFjxcyPrTea6Qn9DPAPgB+JyF+Fbf+cwDA8JSIfA84DHw7f+waBau0MgXz6lwBCY/Ml4Pthv4cikQLwK0xKtJ8Nf1noGOuNKAzzyDf7eW0wS9yF7ekUxYrHlfFgLczBw8d58O69U+5C+7pSnLsW3BlWfB+R5iuA5jMVGi703NAWn1dpt9oYy1f4Rz9zU90FpkdevEhb3GG84JMr++THi9y5p4fv/eQ6wmSxv1rVoOPUD7sux4LVep7bcLbIWKHC9XyWaLnMmaEJ/tc//iv+1Yd+irft3MDXv3ue0TDMqAQ3D4VKTWjNCwxQd1uwJufW7V381M4NvGNnN5s6Ew0XY1tJMvH1hqXtmYe1KtGOiC5Ap6+Okyl6bOyIs6kjWf0nzJUqdCRchrPBYlOvZp4jWqeyEpyKRhaErkbetXfTjNQ1Bw8f5/XhLNcmSjhIMM8WKvj60kmGMkU8DbJQuE6QicERSMZcfA28iZs2dzCWL/Ol+966LHnb6uWHGxjNVwv+RR6bF6ZuSsZcyp4/q8cqwN984ybec+sW3rGrm729HaRCefRihQOtkImvZSx33BKx1o1QxGz5u85fm6A8Wc4HuNEkNQsnmodaSzSyaDaiM+nwUzt6qhfFux55Pkh55Gk1xYsSzLft6+0kV/ZDAUmwbqviK64EcnNV2N4dFOXrC/P/LVfetqP9g/zH7/yYC6M5ejuTnBnKBrVx5lnIG3eDtU8VT/FVSbhCe8LlxOffa0q1FUzL1wkZq4t6MfuK508xQK2yA2vNAEHjBgggW/Q5dXmMTx55KUhOmgmSk8YccAjS9lRCefzpoQn29XYg4RqtfX2dZPIlLmeKCLB1Q7LqHdVmLK+lLe5yejDDwcPHFxSiqxfW+5s3byZbLPPyxXF+eHGMrrY4cl34q4EaJeO0z0IkUK39s/ft54kXzjFeKNORiFUNTmQkzQCtDcwIGQB0JlzODGWDEI4EJRhqF2220g5ML3mwHrmeryDAtWwJRyar24IXzpUEC1OjbA9lz+dL9721ajhqQ021tYl2HZs5V3NtokimUGEwU2hYTn20f5DPP/1ycAwxh9ODGT7+hy/SGyZ6LdVRDAqTHmEkMPE1yLn2yN95G/fs72NTR4IvPPNKNRehzdWsPcwIrTOO9g/yyDf7p6iRPvBT27gWhndAKaywyRUlTMtGa41hq6mVcUefRdkPLuCuA6449HWlqvLnR77ZP6/goN6E/MhEmZ72+Lxy6ornMzJR4gfnR/mX/+VVRidKlLypGdEnrk1KnLtSMd6xq7sqHHj77m5+eOE6//ufvVb9Pt4cquOAqieWTsaq6YRsrmbtYUZoHREVqBvNlatqpNeuZvnrq6dxwhxp8xV+ayUr98iWl+ga7xAaJJFqwb10KgitVTyfc9dy7NnUHi5MzvLgH5wknYqxry89Jd3SRLFczby9ry/N9VypmpA2IhVzOH8ty4vnR3nxJ6P82StX+OHFsbolC2Yj4coMtd97bt3Ke27dOqXf9MWpkfdT69kZawczQuuIx46dJVOohEXegonr6BLi6+SK/bWqNFtrKEHJ9dt298wIqV0dL1al9eP5MtcmgiSquWIQZvvkkZcQoKstzrYNbVPCXI8dO8vV8Txxx2Gi5JEre0wUK5Q9n7/zH/7brMeTjDm0J9xZsxOM5suLXpzaiuzexvJgRmiNUztZPJQpUvb86rqJ2vIHCtW4vRmg1YESVHu9nisF67dUqzcTAF3J4O88nC0GUm4Hyr7SnogFyToVOpIxXh+eoFAO8tZ98o9foqstztnh2TMFBFVzhc5kjPZELDA8E0VGcuVZs04DeD4NZSBYbUXujBvDjNAaZnpYYzhbxFcoh6l2pk/0m+1ZfYznyyTjLp1xh7Hi1MwB40WfK2P5ID2TSDU1k+/7lCs+vsL5sOy0Etx8DE+UGA69pggBbt3exS/esZvbd/fwj37/e2zsSE5Rp/WmU1ybKBFznTmzUTSSgWAlZfc2mo8VvVjD1IY1RIQt6VRQilvXt9JsreAQJExtT8TIhGmLIsVZIkxLc22iRCwsbFcOFxv3X8ngReFXZt58OBKUM79pUzu3bkuzt7eDrlScv3/HG9i/rYs3bOqckSonX/ZwHYebezur843TEWhI1fbg3XurC6VVg0dTxK1dzAitYS6M5mirKf/c1RYnFlXcbNExGUuHT5As9exQtnpTEXm4vmpV8pyvaNgW9J833KrQ15WiMxXHdWYmtZ3NSOwNq8z2ThM1RHzoHdsamtO5Z38fD33wLUEdqXyZvnRqyTM4GCsHC8etYeqFNTyCdDt7ezs5O5RlYoGJP42VR7Ey8284PSQ2fa1VZ8IhX/Kp3TISpPjAK5fGaIu79KaDxa21obDZSkAAfOGZV+hMBZLqoWypKidPp2JcGS9xtH+wYUNkRmd9YEZoDVNvDUhUMno8XzYDtEaYq3LEzu4gNc/F6wUSrhMUugNyZX/GwqtaD0k1KFI4MJqnpz3O53/h1in7nc1IRMapVPHZvzVdLf0elHhvvIaQsX4wI7QGqVXEdSaC0tbRQr+37eji6Zcu25zQCiUqGrgUf57OpMvl8SIOQdZsJShdnnBl1lxtAkGetjCRaMwVNnUkZq20W28xbNT34OHjQTl3k1obc2BGaI1Rf6FfkMIF4FNHXjIDtIIpe8rO7hQD1wuL3kdUrr3iKahSVgDFkWCeqOxPGrkZVWIJKu22Jxz29nZWMxVMp5Fqpya1NhrBjNAqJboL/dHFUfJlRVXpTMboSrokE7EpJbhdR3j42VP0dCTnXMdhrAxuxABBEEqrhKXIp7eH9qhqfOrdj3iqbO0MwnizSaMbWVBqUmujEcwIrUKiu9DxfIlscXJCYLxQYbxQoS1WpKKCg1Srp54eytKdLdZNJGmsPWo9nHoGp/a5M+19QUmnYnNKoxvxclZDobjlKOhnzI1JtFch0V3obJVE85UgrFP0fAoVP0wqCdcmyg2V1DZWP7VGRoD4tP/0N2xsq2auRqA94fKGje1Bu+PMK43e1dNed63QDBXdCpZaRzdz07OFH+0fbPWhrSvME1qFRHehNrdjzEe1YmlULkEg7jgMZ0vV8FzSDeZ/IKjVc9vunnkL2jXq5axkqbXlqFsZNM0TEpHfE5FBEXm5pm2jiDwnIqfDx56wXUTkURE5IyI/FJHbarZ5IOx/WkQeqGm/XUR+FG7zqIQ5RBYzxmqj3l2oYUwn7gSGJxaW1HYkKERU8YP5orgbLFkuVHzG86UFZSZY6V5OI0xfzA0mnGgFzfSEfh/498ATNW2fAb6tqg+LyGfC158G3g/sC3/vAL4K3CEiG4EvAgcIbtpOisgzqjoa9jkEHAe+AdwLPLvQMZp29k0gil+fHsyQKVTWfX0dY3aC9D0OiiKqFCo+HQmXQtlDJMiijggx9fFVuTJe5LbdPQuaE1nJXk4jmHBiZdA0T0hVjwEj05rvAx4Pnz8OfKim/QkNOA50i8g24H3Ac6o6Ehqe54B7w/e6VPUFVVUCQ/ehRY6xKqiNX2/tSrGxI265d4wqkVcTEXMExwmEKY7r8M49G3n0Iz+Np4GSUlXxfUUEdvW00ZdO8vVDd65qo7JQLEfdymC554S2qOplAFW9LCLRN34HcKGm30DYNlf7QJ32xYxx+UZPajmI4tcVTzk9kqG0RAsajbVB3HHwfI+461DxfBxH0PCnVFau50p88shL+H5QEFwIitht72oj5gp96dSs+z7aP8jDz57i9bBK6t7NHXz63v2r3mDNln5otZ/XamOlCBPq3dPrItoXM8bMjiKHCEJ97N69e57dLg8XRnO4AhdG8yZIMKrEw4wGN20OhAWDmQKerwxlipQ8H1SJucL1fImxXBnHCer6iEDF97l0PYePUCx7vP/fHSNTrEyRKh/tH+STR17iek013igPqVoAAA/wSURBVNODWT515CX+7f1vX/UX7NUeUlwLLLdE+2oUAgsfIy3kALCrpt9O4NI87TvrtC9mjBmo6mFVPaCqB3p7exd0gkvB0f5BDh4+zl2PPM/Bw8c52j9IOhkzA2RMQYCetjjj+TLfOzfC6ESR8XwZ1xFu2tzB7o3tOI7D5s4E4/kKjiMkXJdYaE08P1DNbWyPcWW8yOnBLK4wRar82LGzZIsVXBFcxwl/hUwhUJEZxo2y3J7QM8ADwMPh49M17Z8QkScJxAJjYSjtW8C/iRRuwHuBz6rqiIhkRORO4LvAR4HfWcwYTTzXBXO0f5BHvtnPa4NZHBTHcbgyVuB7r1+zaqdGVYgiBNVNY45wPV9BUVIxp5qOp1QJEo96vuL5yuB4UMww7oZiBEcoVsIcckCuFBS9Q2A4W2Jvb2dVqnxhNIfna/B+dBwCFc83FZmxJDTNCInI14F7gM0iMkCgcnsYeEpEPgacBz4cdv8G8AHgDJADfgkgNDZfAr4f9ntIVSOxw68QKPDaCFRxz4btCxqjldSu1k4nYwxli2QLlcl8X74tLTUCBNjSlSQZc7gYpvVRnZwX3NyZpD0RY6JYYTBTpKc9zrWwQqoXpukpeUoCnRKXTrhOtfIqYfE7mJQq7+ppZzhbRMMQXjAuxBzHVGTGkiA6WzpdA4ADBw7oiRMnlny/tQkg2+IuZwazVHzFD/8eFnYzItrjDm/Y1EF3e4KB0RwdYWb01wazpGIOmzuTdIUpdM4OZSlUPFIxl4oXfJ/K075McVdAFUXY0d3GcLYYJDuVQFUXeUJ96RQP3r13xpyQr9DTHl8Tc0JG8xCRk6p6YL5+K0WYsO6YvlrbU8URqhUwDQOCbAZlXxnKFHj21++e8t7Bw8dnrHMpVnySNd5NzHEAj4oPjiP4qmxIxehNpxjKFom5wubOROBdKWztSk6RKt+zv4/fvP/tU9Rx+3rXhjrOWBmYEWoR0xNAJlyHsucjIoia/NoIcJwgAW2pzqRgvdQ5riNsaI8znq9Q8YJ1QI7j0B4Ttm5I0ZdOVVPyROHggdEcN/d2ICJki5MeUGRkTEFmNBNLYNoipqfe6U0nq5PHKbs1WLPItMf58MMEb4nYzH/VeqlzPn7PG4m7LulUDB+l4vv4vtLVFpuxEPOe/X18/dCdfOm+twZlPooVWytjLDt2uWsR0+9iXUfoaY+TdIWBsWKrD89oEjEHEAmKyzWw4LhU8Ym7Qm9nsu779byUt+3s5rFjZ6l444EYIeawZ1NnXePSSHE6w2gmZoRaxPTV2p3JGKrKwFip1YdmNBHHCTwaATa0xbhep2pphCuBFNtXGMoWOdo/WF1AOlcNnIWEzyyTtNFqzAgtM9MvIO/au5Fv/KjIqSuZVh+asQzctKmdz7z/liAP4HiBmAPT6wwKkx5T3A3UbzFXqotDl9JzsRLcRqsxI7SMTA99nLuW5YWz11p9WMYyIiJVL/jBPziJiNCRmJRZqyqnrmR489Y0UrNAdDxf4sXzo8E2wNYNKUTkhj0XyyRttBoTJiwjtaGPbLHCZZv7WXdki0E13Hv293Hb7h52b2xnb29ndZ1PvuzRkXCniFYyhTIXrxcQCYQKviqXrhfIFIJQ3o14LpZJ2mg1ZoSWkaiI1uB4gXPXLNyxHulMTnocsxmAX77rpintV8aCDAlb0ikSroMgiMBQJriJuRHPZS0UpzNWNxaOW0Z29bRz7lqWwYx5QGuZ2YoNBokKJt+Zq5RApHAbGM2hwI7uVNVbujSWRzRYmLoUnoutAzJaiRmhZSASI7x2dZyRibItRF2jOIBPfQOUjDls7UoyUZpaln02A1DbHmVGAKqG6GqmgKjMWFhq/P/t3XtwVOUZx/HvbzcX7gRERbloUNTxggHUqsV6o4pWpZ3RgmOttcM4Y8dbHe3o6B/6R8fadjrq1KoUbWu1SqWOZfAu4oy2FhUBUbmDQ1AUHAgYVJKQp3+87yabmMSEbHKSnOczs7PnnH335Lxv3uTZ856z7+N6Gw9CXey+V9Zw/2vrqd1b79Px9FEZQBlxwOAStlfXNEk4WJoVR4wcAtAwH1tHNf9OWS4JnQ+bub7Arwl1oddWbeWehWvZU+cBqLcrzoobp47nRxWNGeFFGGJTRlw0YSTF2Syjhw/gmIOHcNDQ0vAF5IElnb7g79dtXF/mZ0Jd6PZnVnjw6cVKizLUm5GReOgnkxv+6ZePWMOcNzayuybcyTZrSjnXTT2iyVxsh+43iEtPHM6bG7YXJHW0X7dxfZWncvgWnUnlUH7rs3jz9l7HjRraMISWm/TTOdc+nsohQblPxB6Aeq7Sogw1dfWt3iRywKAS/86Mc93Ag1CB5c+K4HqO3G8jmxGjyvpRlM1Qu9e4eNIonlqymc07vsKA4gyMHNofwO88c64beBAqEDOjcvtX/Pq5lez8qoaaOj8N6g5ZAYo3CWQyDXOz5QJH/nWa3CSxu2v2Ngkw1009ItE6OJdmqQtCkqYB9wJZYI6Z/WZf9rPr61req9zJssodvLupimWVVWzf7TNgF0JRBob2K6K6pp49ebN7ZhRmHDj7qP35dFdNuy74+wV953q2VAUhSVngfuD7wGbgbUnzzezDtt5Xt7eeNZ9Vs7RyB8s2VbG0sor1W6tbvJ6QUUhANri0mAElWfqXZNmwrbrFzJhpUJwVRRnxdW3j9Zch/YqYNaW8yawAnkzNuXRKVRACTgLWmdkGAElPAtOBVoPQhm27Oe6Ol5pMKJkjYPyBg5g0dhgVY8qYOHYYH2//kjsWfEhxVg0pl/easf+gYrZVt547pi8oysDIwaUM7l9CdTuzdHrQcS7d0haERgGVeeubge+09YbdNXUNAWjEoBIqxgxj0iFlVIwpY8LosiYTUgIcOXIwmYyafMIvzojaemNPbT279nwzmPUErc13VpIV15x5uF83cc51iVR9T0jSJcC5ZjYrrl8OnGRm1zYrdxVwVVw9Eljdwu5GAJ+35+dmSgcNyQ4ZMRYzI5MtUibTfcHfAKzerH6v1dV+Xb+76tP6PdW7CvgT2t0OKeBt0cjbolFa2+IQM9v/2wql7UxoMzAmb3008EnzQmY2G5jd1o4kvdOeL2L1dd4OjbwtGnlbNPK2aFva5o57GxgvqVxSCTATmJ/wMTnnXGql6kzIzOokXQO8SLhF+xEz+yDhw3LOudRKVRACMLPngOcKsKs2h+tSxNuhkbdFI2+LRt4WbUjVjQnOOed6lrRdE3LOOdeDeBDqIEnTJK2WtE7SLUkfT1eQNEbSIkkrJX0g6fq4fbiklyWtjc/D4nZJui+2yXuSJuXt64pYfq2kK5KqU2dIykpaKmlBXC+XtDjWaW68yQVJpXF9XXz90Lx93Bq3r5Z0bjI16RxJZZLmSVoV+8YpKe4Tv4x/G+9LekJSv7T2i04zM3+080G4mWE9MA4oAZYDRyd9XF1Qz4OASXF5MLAGOBr4LXBL3H4LcHdcPh94nvCd15OBxXH7cGBDfB4Wl4clXb99aI8bgX8AC+L6P4GZcflB4Oq4/Avgwbg8E5gbl4+OfaUUKI99KJt0vfahHf4GzIrLJUBZGvsE4UvvG4H+ef3hZ2ntF519+JlQxzRM+2NmNUBu2p8+xcy2mNm7cfkLYCXhD2864R8R8fmHcXk68KgF/wPKJB0EnAu8bGbbzWwH8DIwrRur0mmSRgM/AObEdQFnAfNikebtkGufecDZsfx04Ekz22NmG4F1hL7Ua0gaAnwPeBjAzGrMrIoU9omoCOgvqQgYAGwhhf2iEDwIdUxL0/6MSuhYukUcOpgILAYONLMtEAIVkJv4rbV26QvtdQ/wKyA3nfd+QJWZ1cX1/Do11De+vjOW7wvtMA7YBvwlDk3OkTSQFPYJM/sY+D2wiRB8dgJLSGe/6DQPQh3TUqa6Pnt7oaRBwL+AG8ysral+WmuXXt1eki4AtprZkvzNLRS1b3mtV7dDVARMAh4ws4nAbsLwW2v6bFvE617TCUNoBwMDgfNaKJqGftFpHoQ6pl3T/vQFkooJAehxM3s6bv4sDqkQn7fG7a21S29vr+8CF0n6iDD0ehbhzKgsDsNA0zo11De+PhTYTu9vBwh12Gxmi+P6PEJQSlufAJgKbDSzbWZWCzwNnEo6+0WneRDqmFRM+xPHqx8GVprZH/Jemg/k7ma6Avh33vafxjuiTgZ2xqGZF4FzJA2Lnx7Pidt6BTO71cxGm9mhhN/1q2Z2GbAIuDgWa94Oufa5OJa3uH1mvEuqHBgPvNVN1SgIM/sUqJR0ZNx0NiEFSqr6RLQJOFnSgPi3kmuL1PWLgkj6zoje9iDc9bOGcCfLbUkfTxfVcQphWOA9YFl8nE8Yx14IrI3Pw2N5EZIFrgdWACfk7evnhAuu64Ark65bJ9rkDBrvjhtH+GexDngKKI3b+8X1dfH1cXnvvy22z2rgvKTrs49tUAG8E/vFM4S721LZJ4A7gVXA+8DfCXe4pbJfdPbhMyY455xLjA/HOeecS4wHIeecc4nxIOSccy4xHoScc84lxoOQc865xHgQcs45lxgPQs4lKH5R8RVJyyTNkHSDpAH7uK+PJI0o9DE615VSl97buR5mIlBsZhUQAgnwGPBlkgflXHfxMyHnCkzSQEnPSloek57NUEiGuErSGzHZ2wJJBxACTkU8E7qeMCHmIkmL2tj/A5LeiUnV7mz28s2S3oqPw2P5QyQtjMnlFkoaK2loPHPKxDIDJFVKKpZ0mKQXJC2R9Lqko7qoqZzzIORcF5gGfGJmx5vZscALwJ+BC4HTgJEAZrYVmAW8bmYVZnYvYQLLM83szDb2f5uZnQBMAE6XNCHvtV1mdhLwR8Jkq8TlR81sAvA4cJ+Z7SQkVDs9lrkQeNHChJyzgWvNbDJwE/CnzjSGc23xIORc4a0Apkq6W9JphCn/N5rZWgvzZD3Wyf3/WNK7wFLgGEKGzpwn8p5PicunEDLDQpjnbEpcngvMiMszgbkxfcepwFOSlgEPETLtOtcl/JqQcwVmZmskTSZM+noX8BIFyhMTZ1u+CTjRzHZI+ithgsyGH9/KMi1snw/cJWk4MBl4lZAbpyp3jcq5ruZnQs4VmKSDgS/N7DFCBs5TgXJJh8Uil7bx9i+AwW28PoSQUG6npAP5ZjK1GXnPb8bl/xLOdAAuA94AMLNqwqzO9xJmCN9rIXnhRkmXxLpI0vFt1de5zvAzIecK7zjgd5LqgVrgamAE8KykzwlB4NhW3jsbeF7SlpauC5nZcklLgQ+ADcB/mhUplbSY8AEzF+yuAx6RdDMhRfeVeeXnEtIMnJG37TLgAUm3A8WEhH7L21Nx5zrKUzk4180knQHcZGYXJH0sziXNh+Occ84lxs+EnOuh4rBaabPNl5vZiiSOx7mu4EHIOedcYnw4zjnnXGI8CDnnnEuMByHnnHOJ8SDknHMuMR6EnHPOJeb/jn0cK85Iaa8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.regplot(x=\"sqft_above\", y=\"price\", data=df)\n",
    "plt.ylim(0,)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "We can use the Pandas method <code>corr()</code>  to find the feature other than price that is most correlated with price."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "zipcode         -0.053203\n",
       "long             0.021626\n",
       "condition        0.036362\n",
       "yr_built         0.054012\n",
       "sqft_lot15       0.082447\n",
       "sqft_lot         0.089661\n",
       "yr_renovated     0.126434\n",
       "floors           0.256794\n",
       "waterfront       0.266369\n",
       "lat              0.307003\n",
       "bedrooms         0.308797\n",
       "sqft_basement    0.323816\n",
       "view             0.397293\n",
       "bathrooms        0.525738\n",
       "sqft_living15    0.585379\n",
       "sqft_above       0.605567\n",
       "grade            0.667434\n",
       "sqft_living      0.702035\n",
       "price            1.000000\n",
       "Name: price, dtype: float64"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.corr()['price'].sort_values()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module 4: Model Development"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Import libraries "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "from sklearn.linear_model import LinearRegression\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "We can Fit a linear regression model using the  longitude feature <code> 'long'</code> and  caculate the R^2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.00046769430149007363"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[['long']]\n",
    "Y = df['price']\n",
    "lm = LinearRegression()\n",
    "lm\n",
    "lm.fit(X,Y)\n",
    "lm.score(X, Y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question  6\n",
    "Fit a linear regression model to predict the <code>'price'</code> using the feature 'sqft_living' then calculate the R^2. Take a screenshot of your code and the value of the R^2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.49285321790379316"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[['sqft_living']]\n",
    "Y = df['price']\n",
    "lm = LinearRegression()\n",
    "lm\n",
    "lm.fit(X,Y)\n",
    "lm.score(X, Y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 7\n",
    "Fit a linear regression model to predict the 'price' using the list of features:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "features =[\"floors\", \"waterfront\",\"lat\" ,\"bedrooms\" ,\"sqft_basement\" ,\"view\" ,\"bathrooms\",\"sqft_living15\",\"sqft_above\",\"grade\",\"sqft_living\"]     "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "the calculate the R^2. Take a screenshot of your code"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6576951666037494"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Multi_X = df[features]\n",
    "lm1 = LinearRegression()\n",
    "lm1\n",
    "lm1.fit(Multi_X,Y)\n",
    "lm1.score(Multi_X, Y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### this will help with Question 8\n",
    "\n",
    "Create a list of tuples, the first element in the tuple contains the name of the estimator:\n",
    "\n",
    "<code>'scale'</code>\n",
    "\n",
    "<code>'polynomial'</code>\n",
    "\n",
    "<code>'model'</code>\n",
    "\n",
    "The second element in the tuple  contains the model constructor \n",
    "\n",
    "<code>StandardScaler()</code>\n",
    "\n",
    "<code>PolynomialFeatures(include_bias=False)</code>\n",
    "\n",
    "<code>LinearRegression()</code>\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "Input=[('scale',StandardScaler()),('polynomial', PolynomialFeatures(include_bias=False)),('model',LinearRegression())]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 8\n",
    "Use the list to create a pipeline object,  predict the 'price', fit the object using the features in the list <code> features </code>, then fit the model and calculate the R^2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Pipeline(memory=None,\n",
       "     steps=[('scale', StandardScaler(copy=True, with_mean=True, with_std=True)), ('polynomial', PolynomialFeatures(degree=2, include_bias=False, interaction_only=False)), ('model', LinearRegression(copy_X=True, fit_intercept=True, n_jobs=None,\n",
       "         normalize=False))])"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pipe=Pipeline(Input)\n",
    "pipe"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/opt/conda/envs/Python36/lib/python3.6/site-packages/sklearn/preprocessing/data.py:645: DataConversionWarning: Data with input dtype int64, float64 were all converted to float64 by StandardScaler.\n",
      "  return self.partial_fit(X, y)\n",
      "/opt/conda/envs/Python36/lib/python3.6/site-packages/sklearn/base.py:467: DataConversionWarning: Data with input dtype int64, float64 were all converted to float64 by StandardScaler.\n",
      "  return self.fit(X, y, **fit_params).transform(X)\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "Pipeline(memory=None,\n",
       "     steps=[('scale', StandardScaler(copy=True, with_mean=True, with_std=True)), ('polynomial', PolynomialFeatures(degree=2, include_bias=False, interaction_only=False)), ('model', LinearRegression(copy_X=True, fit_intercept=True, n_jobs=None,\n",
       "         normalize=False))])"
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pipe.fit(Multi_X,Y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/opt/conda/envs/Python36/lib/python3.6/site-packages/sklearn/pipeline.py:511: DataConversionWarning: Data with input dtype int64, float64 were all converted to float64 by StandardScaler.\n",
      "  Xt = transform.transform(Xt)\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "0.7513404627519746"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pipe.score(Multi_X,Y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module 5: MODEL EVALUATION AND REFINEMENT"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "import the necessary modules  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "done\n"
     ]
    }
   ],
   "source": [
    "from sklearn.model_selection import cross_val_score\n",
    "from sklearn.model_selection import train_test_split\n",
    "print(\"done\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "we will split the data into training and testing set"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "number of test samples : 3242\n",
      "number of training samples: 18371\n"
     ]
    }
   ],
   "source": [
    "features =[\"floors\", \"waterfront\",\"lat\" ,\"bedrooms\" ,\"sqft_basement\" ,\"view\" ,\"bathrooms\",\"sqft_living15\",\"sqft_above\",\"grade\",\"sqft_living\"]    \n",
    "X = df[features ]\n",
    "Y = df['price']\n",
    "\n",
    "x_train, x_test, y_train, y_test = train_test_split(X, Y, test_size=0.15, random_state=1)\n",
    "\n",
    "\n",
    "print(\"number of test samples :\", x_test.shape[0])\n",
    "print(\"number of training samples:\",x_train.shape[0])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 9\n",
    "Create and fit a Ridge regression object using the training data, setting the regularization parameter to 0.1 and calculate the R^2 using the test data. \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.linear_model import Ridge"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6478759163939115"
      ]
     },
     "execution_count": 43,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Ridge_obj = Ridge(alpha=0.1)\n",
    "Ridge_obj.fit(x_train,y_train)\n",
    "Ridge_obj.score(x_test, y_test)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 10\n",
    "Perform a second order polynomial transform on both the training data and testing data. Create and fit a Ridge regression object using the training data, setting the regularisation parameter to 0.1.  Calculate the R^2 utilising the test data provided. Take a screenshot of your code and the R^2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.7002744265869922"
      ]
     },
     "execution_count": 46,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Poly_obj=PolynomialFeatures(degree=2)\n",
    "x_train_pr=Poly_obj.fit_transform(x_train)\n",
    "x_test_pr=Poly_obj.fit_transform(x_test)\n",
    "Ridge_obj = Ridge(alpha=0.1)\n",
    "Ridge_obj.fit(x_train_pr,y_train)\n",
    "Ridge_obj.score(x_test_pr, y_test)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<p>Once you complete your notebook you will have to share it. Select the icon on the top right a marked in red in the image below, a dialogue box should open, select the option all&nbsp;content excluding sensitive code cells.</p>\n",
    "        <p><img width=\"600\" src=\"https://s3-api.us-geo.objectstorage.softlayer.net/cf-courses-data/CognitiveClass/DA0101EN/coursera/project/save_notebook.png\" alt=\"share notebook\"  style=\"display: block; margin-left: auto; margin-right: auto;\"/></p>\n",
    "        <p></p>\n",
    "        <p>You can then share the notebook&nbsp; via a&nbsp; URL by scrolling down as shown in the following image:</p>\n",
    "        <p style=\"text-align: center;\"><img width=\"600\"  src=\"https://s3-api.us-geo.objectstorage.softlayer.net/cf-courses-data/CognitiveClass/DA0101EN/coursera/project/url_notebook.png\" alt=\"HTML\" style=\"display: block; margin-left: auto; margin-right: auto;\" /></p>\n",
    "        <p>&nbsp;</p>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<h2>About the Authors:</h2> \n",
    "\n",
    "<a href=\"https://www.linkedin.com/in/joseph-s-50398b136/\">Joseph Santarcangelo</a> has a PhD in Electrical Engineering, his research focused on using machine learning, signal processing, and computer vision to determine how videos impact human cognition. Joseph has been working for IBM since he completed his PhD."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Other contributors: <a href=\"https://www.linkedin.com/in/michelleccarey/\">Michelle Carey</a>, <a href=\"www.linkedin.com/in/jiahui-mavis-zhou-a4537814a\">Mavis Zhou</a> "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3.6",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.9"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}