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\n",
    "\n",
    "\n",
    "<b>bathrooms</b>: Number of bathrooms\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",
    "from sklearn.linear_model import LinearRegression\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module 1: Importing Data Sets "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " Load the csv:  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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": 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>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": 12,
     "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": 13,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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": 13,
     "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": 14,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module 2: 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 <code>inplace</code> parameter is set to <code>True</code>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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": 15,
     "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": 16,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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 <code>replace()</code>. Don't forget to set the <code>inplace</code> parameter to <code>True</code>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "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>'bathrooms' </code> using the method <code>replace()</code>. Don't forget to set the <code> inplace </code>  parameter top <code> True </code>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean=df['bathrooms'].mean()\n",
    "df['bathrooms'].replace(np.nan,mean, inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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": [
    "# Module 3: Exploratory Data Analysis"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 3\n",
    "Use the method <code>value_counts</code> to count the number of houses with unique floor values, use the method <code>.to_frame()</code> to convert it to a dataframe.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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": 20,
     "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": 22,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='waterfront', ylabel='price'>"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAERCAYAAABxZrw0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAVG0lEQVR4nO3dfZDdVX3H8c9nd0ESHsayrAwkkYhBEEGe1ocWqzwszmorznRUZKqslpq21iS0RUetU8cZ63SmndokVcuOpWxGxQIFBx1YzRaQBwHdhEAIwXGLAbMoXFfkKTGwu9/+ce+GXbLZXELO7/fLue/XzE7u7/7uveebzd1Pzp57fuc4IgQAyE9b2QUAANIg4AEgUwQ8AGSKgAeATBHwAJApAh4AMlW5gLd9ue3Hbd/f5OM/YPsB25tsfyt1fQCwv3DV5sHbfrukZyStiYiT9vDY4yRdJemciHjC9qsi4vEi6gSAqqtcDz4ibpX0m+n32X6t7UHb62zfZvuExqmPSfpKRDzReC7hDgANlQv43eiXtCwizpB0qaSvNu5/naTX2b7D9l22e0urEAAqpqPsAvbE9iGS/kDS1ban7n5F488OScdJOkvSQkm32T4pIn5bcJkAUDmVD3jVf8v4bUScOsu5rZLuiojnJf3c9k9VD/yfFFgfAFRS5YdoIuIp1cP7/ZLkulMap78j6ezG/UeoPmTzUBl1AkDVVC7gbV8p6U5Jx9veavtiSX8q6WLb90raJOm9jYd/X9KY7Qck3SzpkxExVkbdAFA1lZsmCQDYN5L24G3/TeMCpPttX2n7oJTtAQBekKwHb3uBpNslnRgR221fJemGiLhid8854ogjYvHixUnqAYAcrVu37tcR0TXbudSzaDokzbP9vKT5kh6d68GLFy/W8PBw4pIAIB+2H97duWRDNBExKulfJD0i6ZeSnoyIH8xS3FLbw7aHa7VaqnIAoOUkC3jbv6f6bJfXSDpa0sG2P/Tix0VEf0R0R0R3V9esv2UAAPZCyg9ZeyT9PCJqjQuRrlX9ilQAQAFSBvwjkt5qe77rawycK2lzwvYAANOkHIO/W9I1ktZL2thoqz9Ve3jB2NiYli9frrExrvkCWlnSefAR8fmIOCEiToqID0fEjpTtoW5gYEAbN27UmjVryi4FQIkqt1QBXp6xsTENDg4qIjQ4OEgvHmhhBHxmBgYGNDk5KUmamJigFw+0MAI+M0NDQxofH5ckjY+Pa+3atSVXBKAsBHxmenp61NFRv0C5o6ND5513XskVASgLAZ+Zvr4+tbXV/1nb29t10UUXlVwRgLIQ8Jnp7OxUb2+vbKu3t1ednZ1llwSgJPvDln14ifr6+rRlyxZ670CLI+Az1NnZqVWrVpVdBoCSMUQDAJki4AEgUwQ8AGSKgAeATBHwAJApAh4AMkXAA0CmCHgAhWJDmuIQ8AAKxYY0xUkW8LaPt71h2tdTti9J1R6A6mNDmmKl3JP1pxFxakScKukMSdskXZeqPQDVx4Y0xSpqiOZcSf8XEQ8X1F5LY4wTVcWGNMUqKuA/KOnK2U7YXmp72PZwrVYrqJy8McaJqmJDmmIlD3jbB0o6X9LVs52PiP6I6I6I7q6urtTlZI8xTlQZG9IUq4ge/LskrY+Ixwpoq+UxxokqY0OaYhUR8BdqN8Mz2PcY40TV9fX16eSTT6b3XoCkAW97vqTzJF2bsh28gDFOVN3UhjT03tNLGvARsS0iOiPiyZTt4AWMcQKYwpWsmWGME8AU9mTNEJtuA5AI+Cyx6TYAiSEaAMgWAQ8AmSLgASBTBDwAZIqAB4BMEfAAkCkCHgAyRcADQKYIeADIFAEPAJki4AEgUwQ8AGSKgAeATBHwAJCp1Fv2vdL2NbYftL3Z9u+nbA8A8ILU68GvlDQYEe+zfaCk+YnbAwA0JAt424dJerukj0hSRDwn6blU7QEAZko5RHOspJqk/7J9j+2v2z74xQ+yvdT2sO3hWq2WsBwAaC0pA75D0umSvhYRp0l6VtKnX/ygiOiPiO6I6O7q6kpYDgC0lpQBv1XS1oi4u3F8jeqBDwAoQLKAj4hfSfqF7eMbd50r6YFU7QEAZko9i2aZpG82ZtA8JOmjidsDADQkDfiI2CCpO2UbAIDZcSUrAGSKgM/Q2NiYli9frrGxsbJLAVAiAj5DAwMD2rhxo9asWVN2KQBKRMBnZmxsTIODg4oIDQ4O0osHWhgBn5mBgQFNTk5KkiYmJujFAy2MgM/M0NCQxsfHJUnj4+Nau3ZtyRUBKAsBn5menh51dNRnv3Z0dOi8884ruSIAZSHgM9PX16e2tvo/a1tbmy666KKSKwJQFgI+M52dnTr66KMlSUcffbQ6OztLrgiYiWm8xSHgMzM2NqbR0VFJ0qOPPsoPESqHabzFIeAzMzAwoIiQJE1OTvJDhEphGm+xCPjMMIsGVcY03mIR8JlhFg2qjA5IsQj4zEyfRdPe3s4sGlQKHZBiEfCZ6ezsVG9vr2yrt7eXWTSolL6+vp1DNJOTk3RAEku94QdK0NfXpy1btvDDA7Q4evAZ6uzs1KpVq+i9o3IGBgZkW5Jkmw9ZE0sa8La32N5oe4Pt4ZRtAai+oaEhTUxMSKrPouFD1rSK6MGfHRGnRgRb9wEtjg9Zi8UQDYDCMMurWKkDPiT9wPY620tne4DtpbaHbQ/XarXE5QAoE7O8ipV6Fs2ZEfGo7VdJWmv7wYi4dfoDIqJfUr8kdXd3R+J6AJSMWV7FSdqDj4hHG38+Luk6SW9O2R6A6mOWV3GSBbztg20fOnVb0jsl3Z+qPQDATCmHaI6UdF1jzmuHpG9FxGDC9gAA0yTrwUfEQxFxSuPrDRHxj6naArD/YMOP4jBNEkCh2PCjOAQ8gMKw4UexCHgAhWHDj2IR8BlijBNVxYYfxSLgM8QYJ6qqp6dnxmqSrEWTFgGfGcY4UWXnn3/+zk3hI0Lvec97Sq4obwR8ZhjjRJVdf/31M3rw3/3ud0uuKG8EfGYY40SVDQ0NzejB8/5Mi4DPDOtto8p4fxaLgM8M622jynh/FouAzwzrbaPKeH8Wi4DP0Pnnn6/58+czQwGV1NfXp5NPPpneewEI+Axdf/312rZtGzMUUEmsB1+cpgPe9jG2exq3502t9Y5qYR48gClNBbztj0m6RtJljbsWSvpOoprwMjAPHsCUZnvwfy3pTElPSVJE/EzSq1IVhb3HPHgAU5oN+B0R8dzUge0OSWyQXUHMMwYwpdmA/6Htz0qaZ/s8SVdLauoTPNvttu+x/b29LRLN6+vr2zlEMzk5yUwFoIU1G/CfllSTtFHSX0i6QdLnmnzuCkmbX3ppAICXo9mAnyfp8oh4f0S8T9LljfvmZHuhpD+S9PW9LxEvxcDAwIzFnPiQFWhdzQb8/2pmoM+TNNTE8/5N0qckTe7uAbaX2h62PVyr1ZosB7szNDSkiYkJSfVZNHzICrSuZgP+oIh4ZuqgcXv+XE+w/ceSHo+IdXM9LiL6I6I7Irq7urqaLAe709PTs3Otj7a2Nj5kBVpYswH/rO3Tpw5snyFp+x6ec6ak821vkfRtSefY/sZeVYmm8SErgCkdTT7uEklX2360cXyUpAvmekJEfEbSZyTJ9lmSLo2ID+1VlWjaE088scsxl4RDklavXq2RkZGyy9Do6KgkacGCBaXWsWTJEi1btqzUGlJrqgcfET+RdIKkv5L0cUmv39PQC8rxxS9+cc5joGzbt2/X9u17GgDAvjBnD972ORFxk+0/edGp42wrIq5tppGIuEXSLXtXIl6KLVu2zHmM1lWV3uqKFSskSStXriy5kvztaYjmHZJukjTburMhqamAR3EWLlyorVu37jxetGhRidUAKNOcAR8Rn7fdJunGiLiqoJrwMixatGhGwC9cuLDEagCUaY9j8BExKekTBdSCfeDuu++e8xhA62h2muRa25faXmT78KmvpJVhr0ztWL+7YwCto9lpkn+m+pj7x190/7H7thy8XG1tbTuvZJ06BtCamv3pP1HSVyTdK2mDpNWS3pCoJrwMPT09cx4DaB3NBvyApNdLWqV6uL++cR8qZunSpXMeA2gdzQ7RHB8Rp0w7vtn2vSkKAgDsG8324O+x/dapA9tvkXRHmpLwclx22WUzjvv7+0uqBEDZmg34t0j6ke0tjcXD7pT0Dtsbbd+XrDq8ZENDM1dxZrlgoHU1O0TTm7QK7DNTK0nu7hhA62gq4CPi4dSFAAD2LSZJA0CmCPjMHH744XMeA2gdBHxmnnzyyTmPAbQOAj4z05cpmO0YQOsg4AEgU8kC3vZBtn9s+17bm2x/IVVbAIBdNTsPfm/skHRORDxj+wBJt9u+MSLuStgmAKAhWcBHfSHyZxqHBzS+WJwcAAqSdAzedrvtDZIel7Q2InbZXsj2UtvDtodrtVrKcgCgpSQN+IiYiIhTJS2U9GbbJ83ymP6I6I6I7q6urpTlAEBLKWQWTUT8VtItYk0bAChMylk0XbZf2bg9T1KPpAdTtQcAmCnlLJqjJA3Yblf9P5KrIuJ7CdsDAEyTchbNfZJOS/X6AIC5cSUrAGSKgAeATBHwAJApAh4AMkXAA0CmCHgAyBQBDwCZIuABIFMEPABkioAHgEwR8ACQKQIeADJFwANAplIuFwygYfXq1RoZGSm7jEqY+j6sWLGi5EqqYcmSJVq2bFmS1ybggQKMjIzoZ5vu0asPmSi7lNId+Hx94GDHw8MlV1K+R55pT/r6BDxQkFcfMqHPnv5U2WWgQr60/rCkr88YPABkKuWerIts32x7s+1NthlwA4ACpRyiGZf0dxGx3vahktbZXhsRDyRsEwDQkKwHHxG/jIj1jdtPS9osaUGq9gAAMxUyBm97seobcN89y7mltodtD9dqtSLKAYCWkDzgbR8i6X8kXRIRu0whiIj+iOiOiO6urq7U5QBAy0ga8LYPUD3cvxkR16ZsCwAwU8pZNJb0n5I2R8S/pmoHADC7lD34MyV9WNI5tjc0vt6dsD0AwDTJpklGxO2SnOr1AQBz40pWAMgUAQ8AmSLgASBTBDwAZIqAB4BMEfAAkCkCHgAyRcADQKYIeADIFHuyAgUYHR3Vs0+3J9+DE/uXh59u18Gjo8lenx48AGSKHjxQgAULFmjH+C/12dN32RIBLexL6w/TKxak2+iOHjwAZIqAB4BMEfAAkCkCHgAyxYes+9jq1as1MjJSdhkzrFixorS2lyxZomXLlpXWPtDKUu7Jerntx23fn6oNAMDupezBXyHp3yWtSdhG5ZTdWz3rrLN2uW/lypXFFwKgdMl68BFxq6TfpHp9zO6ggw6acTxv3rySKgFQttI/ZLW91Paw7eFarVZ2Ofu9wcHBGcc33nhjSZUAKFvpAR8R/RHRHRHdXV1dZZeTFXrvQGtjFk2GTjnlFEmMvQOtrvQePAAgjZTTJK+UdKek421vtX1xqrYAALtKNkQTERemem1gf/TIM6wHL0mPbav3K4+cP1lyJeV75Jl2HZfw9RmDBwqwZMmSskuojOcaV3q/4hi+J8cp7XuDgAcKUPYFcFUytXQGkwDSyybgq7gGTFmmvg9lrkFTJayHg1aVTcCPjIxow/2bNTH/8LJLKV3bcyFJWvfQYyVXUr72bVxMjdaVTcBL0sT8w7X9hHeXXQYqZN6DN5RdAlAa5sEDQKYIeADIVDZDNKOjo2rf9iS/kmOG9m1jGh0dL7sMoBT04AEgU9n04BcsWKBf7ejgQ1bMMO/BG7RgwZFllwGUgh48AGQqmx68VJ/zzBi81Pa7pyRJkwex7kl9Hjw9eKk6FwNW5UK8VrgALpuAZ62PF4yMPC1JWnIswSYdyXujYtiIpjiOiLJr2Km7uzuGh4fLLmO/x1ofQOuwvS4iumc7xxg8AGSKgAeATBHwAJCppB+y2u6VtFJSu6SvR8Q/pWyvCqowU6EqsxSk1pipAFRVyj1Z2yV9RdK7JJ0o6ULbJ6ZqDy+YN28eMxUAJO3Bv1nSSEQ8JEm2vy3pvZIeSNhm6eitAqiKlGPwCyT9Ytrx1sZ9M9heanvY9nCtVktYDgC0lpQB71nu22XSfUT0R0R3RHR3dXUlLAcAWkvKgN8qadG044WSHk3YHgBgmpQB/xNJx9l+je0DJX1Q0vUJ2wMATJPsQ9aIGLf9CUnfV32a5OURsSlVewCAmZLOg4+IGySxvCMAlIArWQEgUwQ8AGSqUssF265JerjsOjJxhKRfl10EsBu8P/edYyJi1jnmlQp47Du2h3e3RjRQNt6fxWCIBgAyRcADQKYI+Hz1l10AMAfenwVgDB4AMkUPHgAyRcADQKYI+AzZ7rX9U9sjtj9ddj3AFNuX237c9v1l19IKCPjMsFUiKu4KSb1lF9EqCPj87NwqMSKekzS1VSJQuoi4VdJvyq6jVRDw+Wlqq0QA+SPg89PUVokA8kfA54etEgFIIuBzxFaJACQR8NmJiHFJU1slbpZ0FVsloipsXynpTknH295q++Kya8oZSxUAQKbowQNApgh4AMgUAQ8AmSLgASBTBDwAZIqAR0uyfYnt+XvxvBNsb7B9j+3XllUH0AymSaIl2d4iqTsifv0SntMu6ZOS5kXE5190zqr/PE2mrgNoFj147Ndsf8r28sbtL9u+qXH7XNvfsP0128O2N9n+QuPccklHS7rZ9s2N+95p+07b621fbfuQxv1bbP+D7dslXSDpEkl/bvtm24ttb7b9VUnrJS2y/c+277e90fYFjdc4y/Yttq+x/aDtb7pulzqAfYmAx/7uVkl/2LjdLekQ2wdIepuk2yT9fUR0S3qjpHfYfmNErFJ9fZ6zI+Js20dI+pyknog4XdKwpL+d1sbvIuJtEfEtSf8h6csRcXbj3PGS1kTEaY32T5V0iqQeSf9s+6jG405T/T+HEyUdK+nMF9exT78rgAh47P/WSTrD9qGSdqh+GXy36qF/m6QP2F4v6R5Jb1A9YF/srY3777C9QVKfpGOmnf/vOdp/OCLuatx+m6QrI2IiIh6T9ENJb2qc+3FEbG0M4WyQtPgl/j2Bl6yj7AKAlyMinm+MY39U0o8k3SfpbEmvlbRd0qWS3hQRT9i+QtJBs7yMJa2NiAt308yzc5Qw/dxsSzVP2THt9oT42UMB6MEjB7eqHuS3qt5r/0vVe8mHqR7AT9o+UvVtDKc8LenQxu27JJ1pe4kk2Z5v+3V7WccFttttd0l6u6Qf7+E50+sA9ikCHjm4TdJRku5sDI38TtJtEXGv6kMzmyRdLumOac/pl3Sj7ZsjoibpI5KutH2f6oF/wl7UcZ3qv0HcK+kmSZ+KiF/t4Tk769iL9oA5MU0SADJFDx4AMkXAA0CmCHgAyBQBDwCZIuABIFMEPABkioAHgEz9P6beAVZWSuPuAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.boxplot(x='waterfront', 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": 24,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='sqft_above', ylabel='price'>"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAESCAYAAAD38s6aAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABCC0lEQVR4nO29e3hkd3nn+XnPqauk0qVbUrf7hruN7TY2xuCeBILjdMDZAZKYffYhCZ6BnZmQ2Jth4pAsmZBMArPOZDfesJOBLDNxD8lMCIQEPJB4CBAuTuMwayfYxrfGbbfdNm71TeputVSlup/z7h/nnFJJqpJKl5JKpffzPGpJp6rO+el01ff8zvt+f+8rqophGIbRfTgbPQDDMAyjPZjAG4ZhdCkm8IZhGF2KCbxhGEaXYgJvGIbRpZjAG4ZhdCkdJ/Ai8sciMi4iz7T4/J8Wke+JyDER+bN2j88wDGOzIJ3mgxeRW4Ec8ClVvWGJ514NfA54i6pOisioqo6vxzgNwzA6nY6bwavqQ8Cl+m0icpWIfFVEHhORvxORg+FDPw98QlUnw9eauBuGYYR0nMA34Qjwi6p6M/BB4D+G268BrhGR/yEij4jI2zZshIZhGB1GbKMHsBQi0gf8EPB5EYk2J8PvMeBq4DCwB/g7EblBVS+v8zANwzA6jo4XeIK7jMuqelODx8aAR1S1ArwkIs8RCP531nF8hmEYHUnHh2hUdZpAvH8KQAJeFz78l8CPhtuHCUI2JzdinIZhGJ1Gxwm8iHwWeBi4VkTGROR9wD8F3iciTwLHgHeGT/8b4KKIfA/4W+BXVfXiRozbMAyj0+g4m6RhGIaxNrR1Bi8ivxwuQHpGRD4rIql2Hs8wDMOYpW0zeBHZDXwbeI2qFkTkc8CXVfW/NnvN8PCwXnnllW0Zj2EYRjfy2GOPXVDVkUaPtdtFEwPSIlIBeoAziz35yiuv5NFHH23zkAzDMLoHEfl+s8faFqJR1dPAR4FXgLPAlKp+rV3HMwzDMObSNoEXkSECt8t+YBfQKyLvafC8O0XkURF5dGJiol3DMQzD2HK0M8l6G/CSqk6EC5G+QLAidQ6qekRVD6nqoZGRhmEkwzAMYwW0U+BfAd4oIj0S1Bh4K/BsG49nGIZh1NHOGPzfA/cDjwNPh8c60q7jGYZhGHNpq4tGVT8CfKSdxzAMw1grjh4f576HTnJqMs/eoR7uuvUAhw+ObvSwVkzHlSowDMPYCI4eH+fDDxxjPFtkMB1nPFvkww8c4+jxzdtmwgTeMAwDuO+hk8RdoScRQyT4HneF+x7avPULTeANwzCAU5N50nF3zrZ03GVsMr9BI1o9JvCGYRjA3qEeChVvzrZCxWPPUM8GjWj1mMAbhmEAd916gIqn5MtVVIPvFU+569YDGz20FWMCbxiGARw+OMo9t1/PaCbFVKHCaCbFPbdfv6ldNJuhZZ9hGMa6cPjg6KYW9PnYDN4wDKNLMYE3DMPoUkzgDcMwuhQTeMMwjC7FBN4wDKNLMYE3DMPoUkzgDcMwuhQTeMMwjC7FBN4wDKNLsZWshmF0HN3WeGOjaNsMXkSuFZEn6r6mReQD7TqeYRjdQTc23tgo2tmT9TlVvUlVbwJuBvLAF9t1PMMwuoNubLyxUaxXDP6twIuq+v11Op5hGJuUbmy8sVGsVwz+3cBnGz0gIncCdwLs27dvnYZjGEansneoh/FskZ7ErDx1WuONzZIjaPsMXkQSwO3A5xs9rqpHVPWQqh4aGRlp93AMw+hwOr3xxmbKEaxHiObtwOOqen4djmUYxian0xtvbKYcwXqEaO6gSXjGMAyjEZ3ceOPUZJ7BdHzOtk7NEbR1Bi8iPcCPAV9o53EMwzDWi83UnLutAq+qeVXdrqpT7TyOYRjGetHpOYJ6rFSBYRjGMuj0HEE9VqrAMAxjmXRyjqAem8EbhmF0KSbwhmEYXYoJvGEYRpdiAm8YhtGlmMAbhmF0KSbwhmEYXYoJvGEYRpdiAm8YhtGlmMAbhmF0KSbwhmEYXYoJvGEYRpdiAm8YhtGlmMAbhmF0KSbwhmEYXYoJvGEYRpfS7pZ9gyJyv4gcF5FnReRN7TyeYRiGMUu7G358DPiqqr5LRBJA5zUtNAzD6FLaJvAi0g/cCvxzAFUtA+V2Hc8wDMOYSztDNAeACeC/iMh3ReSTItI7/0kicqeIPCoij05MTLRxOIZhGFuLdgp8DHgD8J9U9fXADPCh+U9S1SOqekhVD42MjLRxOIZhGFuLdgr8GDCmqn8f/n4/geAbhmEY60DbBF5VzwGnROTacNNbge+163iGYRjGXNrtovlF4DOhg+Yk8C/afDzDMAwjpK0Cr6pPAIfaeQzDMAyjMbaS1TAMo0sxgTcMw+hS2h2DNwxjnTl6fJz7HjrJqck8e4d6uOvWAxw+OLrRwzI2AJvBG0YXcfT4OB9+4Bjj2SKD6Tjj2SIffuAYR4+Pb/TQjA3ABN4wuoj7HjpJ3BV6EjFEgu9xV7jvoZMbPTRjAzCBN4wu4tRknnTcnbMtHXcZm8xv0IiMjcQE3jC6iL1DPRQq3pxthYrHniEr5LoVMYE3jC7irlsPUPGUfLmKavC94il33Xpgo4dmbADmojGMLuLwwVHuIYjFj03m2bNCF405cboDE3jDWCadLn6HD46uajyREyfuyhwnzj3hvo3Ng4VoDGMZbAUbojlxugcTeMNYBltB/MyJ0z2YwBvGMtgK4mdOnO7BBN4wlsFWED9z4nQPJvCGsQw6RfyOHh/njiOPcMu9D3LHkUfWNAdw+OAo99x+PaOZFFOFCqOZFPfcfr0lWDchoqobPYYahw4d0kcffXSjh2EYixK5aFZjQ1zt8SOXSzruUqh4VDw1Ed6iiMhjqtqw70ZbbZIi8jKQBTyg2mwQhrGZWK0NcbXUJ3oBehIx8uUq9z100gTemMN6+OB/VFUvrMNxDGNLcGoyz2A6PmdbtyV6jbXBYvCGscnYColeY21ot8Ar8DUReUxE7mzzsQxjS9ApiV6j82l3iObNqnpGREaBr4vIcVV9qP4JofDfCbBv3742D8cwNj9rVW/G6H7WzUUjIv8WyKnqR5s9x1w0hmEYy2MxF03bQjQi0isimehn4H8CnmnX8QzDMIy5tDNEswP4oohEx/kzVf1qG49nGIZh1NE2gVfVk8Dr2rV/wzAMY3GsHrxhNKHT674bxlKYD94wGrAV6r4b3Y8JvGE0YCvUfTe6HxN4w2jAVqj7bnQ/FoM3jAbsHephPFusFfSCzV0OwPIJWxObwRtGA7qpHMD8fMJLF3Lc9enHOPTvvr7mteSNzsIE3jAa0E1NL+rzCdlilYszZXxV8qWqJY+7HAvRGEYTNrru+1pRX174Qq6EgyAOVHy1WvJdjs3gDaPLqS8vXPZ8REAVEm7w8bfkcfdiAm8YXU59PiHhOniqqMJIJgls7uSxsTgm8IbR5dTnE9JxB0eE7X1x+pKxTZ08NpbGYvCGsQWozydsdNNwY/0wgTeMDWKjvOndkjw2lsZCNIaxAVitG2M9aFngReRVInJb+HM6auZhGMbysVo3xnrQksCLyM8D9wP3hZv2AH/ZpjEZRtdjtW6M9aDVGfz7gTcD0wCqegKwIJ5hrJB6b3qE2RWNtaZVgS+pajn6RURiQEvdukXEFZHvisiXVjJAw+hGlqp1c/T4OHcceYRb7n3Q6sUYK6ZVgf+WiPwGkBaRHwM+D/z3Fl/7S8CzKxmcYXQri9W6sQSssVa0apP8EPA+4GngLuDLwCeXepGI7AF+HPgd4FdWOEbD6Eqa2RXrE7CA1YsxVkyrAp8G/lhV/zMEYZdw21IZof8A/GugqeNGRO4E7gTYt29fi8MxjO6lvjhYhCVgjZXQaojmmwSCHpEGvrHYC0TkJ4BxVX1sseep6hFVPaSqh0ZGRlocjmF0L1ECNluscHIix/Fz07wwnqMvaesSjeXR6jsmpaq56BdVzYnIUun+NwO3i8g7gBTQLyKfVtX3rHCshtEya71KdD1Xnd516wF+9f4nmcxXcAQEqPrKRK7E0ePjFqYxWqbVGfyMiLwh+kVEbgYKi71AVX9dVfeo6pXAu4EHTdyN9WCtk5TrnfQ8fHCU7b0JYq6gQNx12D2YZiAdt4VQxrJodQb/AeDzInIm/P0K4GfaMiLDWCVrnaTciKRnruzx6pE+RKS2TVUtDm8si5YEXlW/IyIHgWsJ7hiPq2ql1YOo6lHg6EoGaGxONrLJ81onKTci6dltTb+NjWFRgReRt6jqgyLyv8x76GoRQVW/0MaxGZuUKKQRd2VOSOMeWJHIRxeL589PU/GURMzh6tFM04vGWoljdNyJbIkL2RI7B1JkUvEV72853HXrAT78wDHy5SrpuEuh4lnddmPZLBWD/5Hw+082+PqJNo7L2MTUhzRypSrnpoqcvpzn7j//7rLj1tHF4qULOaaLVQoVj6l8hZcv5prGwZdaJbqc445ni+zsT1L1lbHJAtOF8qqbZLSySrWbmn4bG4eoLl5xQEQc4F2q+rl2D+bQoUP66KOPtvswRpu55d4HGUzHyZWqnLlcJAgjK54Pgz1xRvqSZEvVlkI3dxx5hPFskXNTRaqe4jiCr0rMEXYOpBjNpPjsnW9c8LrVNrWIjhvdBUwXKpzPFlGFN+wbWnHIqf7upn5mbuJtrBQReUxVDzV6bMkYvKr6IvKvgLYLvNEdRCGSiWwJEXBE8H1wRbmcr5ArVnn1aF9LoZso/l32fNww4SgSNI9eLA6+2qYW8+Pu/ek4mVSMqUKl4QWlVWyVqrGetGqT/LqIfFBE9orItuirrSMzNi1RiKRY9QDF9xUfBREcAU+15Rro0aKfhOsQ3WyqQsJ12hoHb1e1RysTbKwnrQr8zwL/EvgW8Gjdl2EsIIof9yZieD7EXGHXQBo/VOiEO/u2W0rcootFJhXDR6n6Pr6v9KdjbU06rkUcvxFWJthYT1oV+NcAnwCeBJ4A/gC4vk1jMrqAwwdH+fi7X8+uwXToPonhOoKvMJJJ1p63lLhFF4v9w30MpGKk4y4DPXGu3N7X1rh1u5KcbzqwjbHJAt87O8WL41ku5IrmjjHaRqsLnf6EoNnHx8Pf7wi3/XQ7BmV0Dsv1s89//rvesJuHT15ibDLPldt6uDhTxnUEVV1g/Wt2rI1qEr3Wxz16fJz7Hz/NUE+cbLFKqepzaabC+w/vs/i70RaWdNEAiMiTqvq6pbatFnPRdBaNHB9ThUpTF0wrDpFm7pat4C6Z78wByJerTZ1AhtEKq3LRhHxXRN6oqo+EO/xB4H+s1QCN9rGaFaXzHR9Vb3EXTCsOkY2qgb6RK2sj2rUithP+to3GzkFjWo3B/yDw/4nIyyLyMvAw8CMi8rSIPNW20RmrYrVFsuY7Pi7kSou6YFbjEGmHuyRaUHTzb3+Nuz79GC9fzG1oh6R2JFit+5Odg8VoVeDfBuwnWNn6I+HP7yBYzfqT7RmasVrqZ8Wt2hLrmS9IZc8HmrtgViNgay1+9R/6YsXHV+VirkKuVF32eVgr2uHMWe3/cTdg56A5rRYb+367B2KsPcsJCTS6xZ1fD8V1hKqnTV0wS9VPWew2eq1rr9R/6KNFUgpMZEtkUvFF7w6WW/umlXMZhabuCce20hW287HuT3YOFsNaxGwylhNrbLXoVtPiYLdfzz23X18TpKVcMIsJ2FIFyNZa/Oo/9AnXoeop4szehTS7O4jGWa56TBerwXPLXq32zVIF01r5O9cyNmxVJ+0cLEZLLpr1wlw0i7Ncp0mrz1+Ou2OlNV5W4yBZSQKt/njThQpnpoL+NA4gIlR8n6tH+vjQ269reC5arX0zf2yTMyUqvq55DZvFzk23u4+WYqufg7Vw0RgdQCtOk8V86JEgQyBk0XNOjGfZ2Z+ac6xmt7grnYGemszjCpycyFH2fBKuw3BfYsnb6JWWHq4P+WRSMbZXE1ycKVPxlURM2DOYpuLrgn0tp/ZNo7G9fHGGPYNB++LowiKAr7rqssmNaEfYZ7Nh56A5JvCbiKVijY0E5/7HTy/woc9/TrZYJe6WGO6bFfm1vsXNJGOcGM/hOhLE8n3l9OUiV4/2Lfq6ldon53/o9w/3MdRTpuz5C+4i6vcV3e7XwjrSvPZNo7HFHYfz2RL96UTgOiJoqpp0nbYVFtuohWCdhJ2DxrRN4EUkBTwEJMPj3K+qH2nX8bYCS8UaWxHDRs8Z6olzaaZCTyLWluYSH//G8zx3PouvQfNoh6DPKARt6BZ73SMvXUQVHIHh3gQ7BtIN7y4WS2xGRGWM65m/r2jmn0nFuDhTxvcVFPp741Q85U0HttXufiayJXb2J+fsb0d/krHLRfLlKmXPJ6iULLXE9Hon/8wfvrVp5wy+BLxFVXMiEge+LSJfiRZLGctnKadJK26CRs8Z7ktS9XxGM6lFb3EXE4tmj338G8/zsQdfwK/TcR/w1Wf3YJqZ8lxrZET0ukj/fYXxXBmATDo+ZybdLIzzrrHLPHzyUm1MfYngnNVfIC/OlJgpedxy74O1cUeJ5ao3Tb7iUwlLCoz0wp8+8n3603EG03EuZEucvlwEhP7wnMZch2tG+xjsSTA2GYRn1rMTVD1r3VnL2Hy0TeA1mJrlwl/j4VfnZHQ3IUvFGltxEzR7ztU7+hdNdtaLhSvw3Vcmed+nvsPVI32847VXcP/jpxsKySe//RKOgIahDsLvngZiOJpJzTlGdJE4O1UEVVwJnhsxniuTSsTm3F00uiu5kCvyiaMvsmcoXRvTdKFSewOm4y4XZ0qMZ8uM9CXmXhjesBsI3qyerwxnEmzvTfLCRI6qp/QmA7/1zoEUY5MFzmeLZFIxLuRKTOYrZFIxBnsSvP/wVdz/+OmmrqN2Y7XnjbbG4EXEBR4DXg18QlX/vsFz7gTuBNi3b187h9MVLBZrbMVLvlK/eSQWnq+cnQoaebgivHwpzyeOvsi23jgD6UCs64VkpuwRc0Ach4rn1y7xvrLAI18/2xybLDQdy7vesHvOncPjr0zi+T7JmMtIJkkmFWcqX8EL3SzZYoWJbIli1SPpuoz0JZkqVJgpeYz0JRjJzI67/sJQv0AqGXPxfMWRwEsPwXdflVJF+f7FGcqesq03zvbeZC3/0SjJ3Q5xbXQHZf5wo60Cr6oecJOIDAJfFJEbVPWZec85AhyBwCbZzvF0O624CVp1HMwXjMhp89KFmVqXpmiGqwpT+cqcJG0kJL1hWCTmCOBQ9X1UwXVkTvJ3/mzTEWphndDMUrs4PHzyEncze1GIxlP1lTOXi+wahJLnk4q5ZIuVWtvAmCOUPZ+Zssdvv/MGfvOvnlkggPUXhvkLpBJucJEqVf05+3QcqYl7dA6ii9zDJy+1vZBYs1BMo5CU+cO3FuviolHVyyJylKDkwTNLPN1YBa24CZZ6zmJOm7Ln4zqB4kbuEiDs3jRLJCQ/dfM2PvbgC1R9H0cCYReF22/cyX0PneQ3/+oZ9g718Pz5aa4YSNdeP9ybqMXcI2FXYLTOWhldFHZkUoEdUYNnnZsqEnMcMqnYgraByZjUlrE3CldFFwZYuEBq10Ca05cLwUzeBTQQ/x2ZFKcvF5pe5NpNs1CMiFDx/DVbHWxsPlqtRbNsRGQknLkjImngNuB4u45nrB33PXSSctXj3FSR585nOTdVJB13uDRTwRUJWvCFbfiG+5JkUjFijtOwxsrdt13DL73l1aTjLlU/EL3bb9zJY69MzSkOlSt5XMiVamPYMZAOZ/2BsIvAjkxyToI1KlDWn46zayBNzA0EV4H3H76KRMylWPEoV30KFY+S55MvB3/XifPTDWvDRBcGCJLPPornK3FHiLnCYE+8dpyoU1V/Ok4y5lAKV8lGrNdsuVmhtlyp2pamJcbmoZ0z+CuAPwnj8A7wOVX9UhuPZ6wRz5+fZrpYxUFwJag/U/Z8UjFhz1Avz4/niLuwK5MKxc7l/Yf3NY01333bNdx92zW1/d9x5JEFM85tvYFVszc5a9Uc7IkjBA2v03GXC7kSY5MFpgoV7jjyyJwQRH86Tn86XlsdGx3v9795Yk5mXwnuNqph/Ke+FMOeoR7e+bpd3P/46TkLpCbzFXqSMUYzKX7rx1/DfQ+dXDDzz6RiVPO6otnyaq2MiyXXzR/e3VTmTSrm004XzVPA69u1f6N9VELbiuPMruT0fUXE4SsfuLUmSCfOT1MOi3F95ZlzqGpLNqlGyb/tvUlyxSoT2RIzZY+k67CtN06h4jGRLaGqC5KY810xkahGXvXHX5kkHsbHAaJQvu/Dtkyc+x46yWfvfOMCAbxxz+CcBVK/20Bw5yeqE7HFL3LNWAsr41oXautktrKv3/eVUtWnVPUoVX2KFQ/PX/wTZytZjQUkYg6Fsoevsys50WA7zArPhx84xoArVD2fE+OBI3b3YKqpSEUfzolsiQvZ0hx/+IVcibKn7BlKMeT5nL5c5Nx0id2DKWKuwyuX8ghwIVdmulBlJJOkPx0n4Tqh5zwQ1Tcd2FazbHq+j1PL0AZE8fjtvckFZQdaFY7FEtV3L/Ncr4WVcass1d9qvv5IyEuVQMyj2XrF8zk5McPxc1meP59ddB8m8Jucdsxorh7N8PLFHNOFaq1uTH9v0Og6ol6YTk4EJQjQQIAPjPTVRCp67onxLNlilaGeODv7k5y+XGRsssDuQSXmOkzmK2zrjTfc33BfsnZXkYo7NbfMFQNJylWfr3zg1tq46sM/yZhL1VeE6LVurXBYfXx8JcKxVqGPtbIyboVQTDf7+quhOyuamZergUXX85VXLuV57lyW585lOX4+y8mJXO3zsBQm8JuY1c5oml0c7rr1AL96/5OhBTJ4k1Xn3fLXC1O9syYqx5uOu5w4P10bX75UDTzlM2V2DaTZPZjmfLbI6csFehIxKp7PVD7wm8/f34VcCSFMtiJByAjl/HSJ1+8bmvM31Rc1K1V9PF+DYl9A1fcblh14/JXJ2orTqGHEegmHlbptnW7x9atqbWZeqnoUK35oH1bOTBVrYv7c+WCGXqw0jrMP9yW4dkeGxZp1mMBvYlYzo5l/cXj5Yo67Pv0YfclgIVCp4oEEpXV9VSbzZT54/5O1xhf1wpRwg1k1dbbJQsWj7Clxz+Nirkq+4tdE+pVLeXoSLpmky6W8z0gmiecrxYrH9y8GH9aqr7giJGMOZc+v+eKDXEDwIanqwjhzVNRMJKjgGLlqHAnEIBFzuHJ735xQjq/BRSDy0C/VEKT+HK7k7qn+dZlkjKlCBaDr4+erZbNeDMvz4uYVL5g4TWRLNSGPwi3ZsAfBfPpTMa7dmQm+dgTfh/uSuI7wZ3c2P7YJfAfSqnDUz2ii1ZqlqsfYZIGjx8cXFZv6i0O2WOFiroKiFCs+L1/KU/WUPUNpVAk85hLMwl+6kOPn//RRqt5sQjXmBHF6EWFnX5ILuSKXZiqUqz7TxSoxR2riTvi96ikTuQpxJ7gw9STcuXVpFKqqDCUcvFKwmGqkN85M2avN8A9s6wXmlj7OFiuB+M+f9Ch8/N2vr52T+lBO5HdHZjs+LSUcK717mv+6QsVDgLgjTBUqK46fb4Xk42ZIJnu+BmJe8WsJUc9XpgqV2TBLKOYXZ8oN95GOu1yzo49rdmQ4GIr6FQMp4m4wQUnGnNr3mLu4090afnQYy2leEDWn8MKYdDSzdRxZ0vMcVVYUEU5O5MJYddBQG4L9+Eqt2FfcDVaKikjDzL0rMNQTZ6pQpeIrcQeqfhRWCah/VTLmUKr6JGMO1+zIBCGVihcWIoNUzAmSpI7D/u1BJ6nILhmdk3e9YXdtFh5tf/livqmzIB132N6XZO9QD8fOTFH1fCphdctKuCJXCY6dScX4vXe9run5a9TAZCJbJF/26E/Hm4rsahqfNGMrNbxYacOZdlALtUQz9IpfW1h24nyO43WhlrNTxYb7iLvCVSN9XLtzVsz3DvWQjM2KeTLukHCbi7k1/NgkHD0+zt1//l1mylVSdXVVmoVdohnN+HQR0DkrK2Phas1W2vlFZW1rdsK6MgER0WM0mBAIwQrVbMmrlRWoDxs2klsJv0pVnxPnsxTDKXcwmw3252mwsw+9/TpgoUukUYjKFWhcnxIKFZ/z00XGp4uUPcV1grLFVc+fLYsQ/rPUtGd+PHi6UAnKC6uyb1tP0xl9O+LI3Zx8nE99aYtTk/laIn89/s6oTEWp4lGs+kHopeLxwniuNit/7lyWVy7lG75/HIH9w71zwiz7h3tJx2Mk4/Uzc7eWg1otJvAdQjQLyxWrqMBM2SN/Mc9IX4LR/lTT7kr3AHd9+jGUYDYw3JeiPx1HVRcVjfrbXQcohQIed+aKs9NA7OejRBeA1u4Gg7g3OA54PjVxj/ZV8UGrQdxdoNYfdv4Mt1EtmUYXoHpcEUrh8TwfYvU1b4BXbe9peFGdHwKZX+clWoWbirmLJmrbEUfuluRjK6yXVbLec16sRDP0IEd0vC5u/tKFmaZ3jHuH0nPi5q8e7SOTis8R8kTMWTMxb4QJfIcQlQdQZjVKgYlcGccR9g837nx0+OAob9g3tGzRqPdOv3IxH87CwXUdKn44m27T+67qKTE3jNvT+LIQ1Htx2DmQwnVm70bqhXa6UKHq+bVqkDD34tSIYoPgvK+zdyGRL3+pTlnzF1mVqsFdUNTYY/4+ItoRR96syceV0K67lfokaORuGbtU4Lnzs3HzFyZylBe8fwJ29Cdrs/Jrd2a4ZkeGoZ5EEGIJhbzdYt4IE/gO4dRknmyxGoQm6mYECkzmK/xuAwGIxO7YmSmypSqikE64ZFIxEjF3SdGIvNM3//bXyJUC760XirsrwUrWgzv7OXZmqhafXopmgj3/cQfwQ4dOM3YNBguhVJUT56d52+9/ixMTOeJOcDs7XawyXaxybrrEaF8Cx5EW7yFm8QnyAb5qzQEES3fKAuYssupJuPQm3doFYv4+ItqxKGkzJB/XirW4W5nvOS9VPM5OF+ckQU+czzZtRjPUE58TZrl2Z4bRTGqOmCdjTm0l+GK0OzluAt8h7B3qCasgSlj61q/NLDOpWMNyvh9+4BgVz6NQ9mqhlHw5qLPy/sP7Wrbs5UrBqtVkzAmbcQTFxDxPOXZmqhbCGEi6QezRa1ySQOq+NxNaJZgpl5ZYqCHCnFWu2ZJHrpzHlaDk7/yZeK3y5DIRwPd9fIT+dKxhY45mojJVqNQWWUX/H62I7FovSlrORWOzu22We7dS7zkvhonQ8excr/lz57JcDq2q8+lNunOE/NodGXYPpkklXJKuW0uAtiLm81mPcJMJfIdw160HgsYVqriOEJdAbPuSLuWqzy33PkgmGQhQruwxXajQm3SZLlRxHCEmDr4oMVfIpGJ88tsv8bnHxuZ8iBt9uO976CTbeuOhTRLEAb+q1OSzToenSs3Sl3NZahbtok0ToRECNbG9OFMOE7Jr7/iKksd7BsNGIQ2siq2IykaXC2jlotENS/2Xulspz4ubT86U57hZnjuXZTxbarjvZMzh6tG5jpYrh3tJxV2SMTeIna9QzBuxHslxs0l2EB//xvN84uiLVH2fpOuQirtMFauM9AWxvKD/Z1Dv5fTlYm2FZtwVhKAtXCXsOuSrct3O/jmWwk898n1ypSqeH1xE+sILxhUDaXKlKuemig1i1Msj6UrTGf5y6E/FGEjH6U24PD+eI+YENs2lEr6t4kCt5G8ijPVH1kJgzoUwWhRVrnpki1VK1cCH//7DV82pktnptMOiuRHUJiqXZtg12MN737iPQ/u3cXmmwvPng+X80Qz99OXGncFijnBgZNbRcvCKfl492ktvIj5rT4w5SLsSUcy1KkeoBp75v/u1t7S8n8VskibwHUa9z3eqUKEn4TKSSdW86mggShAs4/d8Je46YUMLpeIH4pNwHQ6MBInZfLnKmcsFKl6wOtQPLwQRcQeGehJcLlSphB2XVsP8PMJK+JXbrubu267hjiOP8N1Tk6jPgnGvlqgBSVTTPV+uknAdZsreAk/5zfsG+PIz52sX34GeoNzBSF+SiVyJctUn7grX7OhfdsPy9WKtBGUjmL+8P1eqBmJ+blbMX74403ACIMC+7T3BrDwU8+uu6KcvGauFWNot5o1Yqwuu+eA3EfW32tEHEhbWe9k1kA5WmBJYunyps97oQjdHoeKTCBtVzBfJih/EsF1Z0mXYEqsVdyFo1n3jnkFOTebZkUlyZqoYVoZsvu9olt/K/pXgvAz3BVUpCX8/MZ5jz1B6wW3zN49PzNk+XagwkS0EbppwNVehAi9dyC0Ie3RKaGQzuW0iz3mx4pEve7xwPsuz57I830LBrV2Dqbq4eT/X7+pnW29idhWou/5i3oj1SI6bwHcwi9V76U/HmS6WmS56eKoIQk/cQeLS0M0RvZ2rfvMQjK9L+94brUqtJ+7IqmbZMUdwnSBZ/OEHjpFJxmoXtAu5Es3SAG64incpBGo+5JFMcoHYAQ27I82UPfbVbb+QK+FIUDs/EZu9g8oWq+wciM2Jo3bKQqROddvUe84LZY+XLszwvbPTs46W8aULbl27M5iZ37h7gOFMsvZ/3Cli3oj1yNu0TeBFZC/wKWAnQaj4iKp+rF3H60bqP5DDfYlaDD6q9zJd9BjpSzDcl1ywfH/+h3jPUJozlwuLircjsGeoh9OXC01nwvWlB+aTcJ1wBarX0sUC5jpuYo4QDxtb+6qcuVwg7ga5gv50nP3DvVycKXF+uoQrgEiQkHaEwZ4EpyaDO5qY09wPH9zB+PzLW65qeJ72b+9p2Ki6d97Cpmj1L8yuFxAJts+37XXKQqSNTgRH1HvOxy4VeOb01LILbl23s5/X7hlg92B6zqKhzUa7yzy3cwZfBf53VX1cRDLAYyLydVX9XhuP2ZGsNP46/wN59WgfqspM2WOmFIh7tMinJxFjIlvkD8MFU55CTIRXj/bxWz9+kL96YowvTjZOOEX4BPH93YMpzk4Fy/njjnDVSC/veO0VfOzBF4LyuwJxJ/gwReWBndoMOrB5+hrUdllK5OsfEpSK51P1lZgEC6/KVZ9JrxyWAQiel3CE0f4kIjKnbEHF85nIlZcM0zgID5+8xM37Bvjm8QlmyoGA/9wt+7lxz2DDWe5bD47w5WfO4/mBnRRVfKRmLY0aoyRcZ0HYo5NCI+tdN75aF2o5P13k6dNTHD/besGtSMxv2DPA/u09pMK72c0o5htBO1v2nQXOhj9nReRZYDewpQS+3q8+la9wdqrA469MtuTAWOzCcMu9D9bqngfldIVy1UeBhCs1YZ3IlXhq7DIPPHl2ybH2JRxGMylOnJ9GRALB9ZXvXwpmmoPpWK1Xa9X3qQ+B+hqV5/XmCPtyYvrBrDt4geM4eJFrZt4iK0+V05eL7BlK1+qRPH9+Gkdai/8rynPnpvmHlwOH0r5twaz9/sdPc+OewQV9WiMXzbbeOFP5CsWqhyNCIoz5R1ZL14GhVHJB2GMjQiMbkdT1/aB3b7HicWmmzNNjUxwLQy3PnctybrpxwS0B9gyluflVQ1x3RT837hngqpE+epIm5qtlXVw0InIl8BBwg6pON3teN7po7jjyCC9fzHExV0EkmOl5vuKIcN97bm76oYsuDPXWPBHYkUlSrPpcCme0QUgiEJr6pfMQ1k53gucUlljDL8C1O/r40Nuv484/fXS2uFjd47sHU1R9ZTJfqdVzWYy4KzgIpSUaAy+GG9aracZIX5xs0VuWvTPK1ToSdHmqdxs1cjA0cjtcyBW5kCvPKdIGwT6vGQ3OY6N2hesRGllNdcnlXBgiz/lUocIzp6d45vR0SwW3RjMppgoVUnGH/lRwTj2Ff/uTr+G21+xc7Z+/5dhQF42I9AH/DfhAI3EXkTuBOwH27dvX7uGsO6cm80zlK7VeoBAkBKu+Lkiyza+z4oRFxxyCeuoVTxm7PHcWNN8VIwQXkNpKWK+1EgMKXJwp84uffXyBuEePn5kqsq030VIyk3C8gi5ZvmAxlro2TOQar0BsRrCAKtx3uPJ3ulCplSJuFBdvFEOfyldq6wnq3TuqcG66yF2ffoxMKlZrkLKeoZGVJnUXc/v88DUjlMLJxvFz0zw9NsWzZ4PFQ60U3Lruin5eu3uA63cP8Mt//gQXckV6k7PnNF+u8kffftkEfo1pq8CLSJxA3D+jql9o9BxVPQIcgWAG387xbAR7h3o4O1UgXlfnRDVwcsxv+lz/4To7VcD3A5+24wrlpapoRfuGWnPe6PdWmSkFnZeaEd0xNLoALDaeTmL+eBQYu1xgD0H+oVFcPJOM8cJ4Di+sVzOSSdbuShyEivqzZYYVpgpVYmGbwo2wRDZL6p44Pz2nOcr82XmQx/C4kK1Q9nzijpBOxPjIA8e4+VVDSxbcGs0kOVgn5jfuCRwt82uZn5kqdETSeSvQTheNAH8EPKuq/75dx+l0aiUIfK35zH2CxUlThQq33Psge4d6uJwvz5l1pWKBNS9YyLQ8oVypqC4m7hCENvrTcS7MlNfEL98peL5y+nKBXYPpBXHxo8fHGZvM10JAFc+jeCmPEy7mipKrQO3ER/V2Kr4umD2vR2y8UVL3Qq7E5UKF754K3osXciU+eP+TfPRdr+PNVw9TKFf53pmpwMES1sOvVpVCtcKlfKWWh4mICm5Fbpab9g6ycyDVUi3zTko6dzvtnMG/GXgv8LSIPBFu+w1V/XIbj9lxHD44yvsPXxWWIAgcGHHXIVvyGM0karfCL1/Ms2dwtuztSCZJ/mI+XLfUIWqq8P2LM10l7hFVP7CYzhfb3/3KsxQqflA3niiZDINJl1zZD9YgzFsg5tQ5amB2drpeC54aJXUv5ErB5MILSll4VZ9LlTIf+Ivvcs3Ofp4/l2Uqsic2+P+9ed8gB8OZ+U17B9m7LU0qHltR+dtO9eN3I1aqYJ2YX4KgN+ky3Dcr6CfGs3iRG6YuBLLe/ztJl6aLiTYjy4n/xx3YMZCeM7O+9je/gqriOrMhBs/3EZHZC7cXOIpcB3w/uNOJOU6t3HGUvAXWrRbM0ePj/OG3XmRsMs9IJsmTY1O1x1pZmyAyu95hIOXynd/6x2tay7yTWu9tdqwWTYfRqCbIuan8shOGa030wW619vtmoxWxH+mL059OMFWoMNKX5Nlz2ZpTKYojl6tB79jdg2n6EkEHp6gejapP2YNtvXG29ybnOFiiDlT1/+/ThTLnpkuMZJLLCtk0CvXccvUwl/Jlnjo1xVNjl3nmzPSiBbeic3L7Tbt47e4Bvvj4aSayBXIlj4of5BsyqRj7h/s2VTGyrYbVoukgjh4fZ7pQ4dxUkWTMoTfhMl2srrqK41qgBOGFhOt0xHjWmlYuWhO5CtPFoPnJxbDG/KxTKei45Wmw1mAwHQ8F3OejdQ26m81O9z40N/acLVY4fblIbJkhm6PHx/mtv3om6MDlCM+dm+YXPvM4A+E+llqtHE3EfR+u2dHHx979egBePdIXlIdIJyx00iXYDL6NfPwbz/PJb79UWyn51oMjPPbKVOhUKOOjS9oAN4LV2Bq7jfnnIuoTu3dbT11DkiIzJY/+dHzRWfh8f/oL40GF0N2D6VrBs0Yhm6PHx/mPR1/gpYsz9MRjTObLFEPxbfb/tGswxcGd/bx2dz+v2ztIrlDlni8dI1fyqPo+MSeYnb/3ja/i4ZOX5pRGfvjkJQudbCJsBr8BfPwbz/OxB1/AkaA2Sq5U5YtPnEWAnoTLtt4E55s0HthoTNxnmX8u4q7Dzv5kTdyzxQoXsmUU2Leth5cu5BZ44GG2vnwU0pkKe7ruHkzVxB2ChOwrF3M8e3aaJ165zFePnePhFy/UifnSnatSMeF33nkDP3rdjjnb+9Pxhit065O+9z9+uqUFUcbmwGbwa0h9XPTsVBFUScRcqp6/pnXMjfVHgETMqTU4r3qB1TAf9u1MuMLOsIQzBFU1rxhMM1WoBDF8V5jKVyh5wez5/Yev4stPn+XkhRxeuN6h9l7xFK/J59IJLYyisNjN396h9JI13rulAchWx2bw68B8C9xYWNjL85t/WI3NgwIjfQku58u8ODGz4PGqwrmpQrDq2KHmgT99uYDvB0XYRIJVzBXP5/e/cQKnrgyD5yllb24VxXTcper7ZJIx0okYPQmXhCuMZ4tLJuTHligsB51T5dJoHybwa8T85eFRsa2K5y+4zbcY9+bDEShWfar5xiESUaXsQzImNQ981fMpV/2wdPLCdoPz8y9CUBL3V992kNfvG+TaHRne+0f/sGCWnS15xJ0g2dvsxrCV95ctOOp+rEzbGnFqMj+nUcRwbwJo/EEzcd98+Aq5YpXL+WDmHPnEa41UQmtpORT1UtXn2XPZmgA3ayV3RX+KA8O9vOaKDNfv6qc/Hec9b3wV1+8aIOY63HXrASqeki9XUQ2+Vzxl12CafduaC3E6vvRHu9m+zTXTPZjArxF7h3pqHYEgWDDjdmYjGWOF+Ko1+2i0VmC+bodNt5asSS8EdwXTxQqnJvN8/2KeizOlBbPnwwdHuef262sVGEczKa4Z7SPmOmRS8QUhlmjfv/AjVy359zTatyVYuwsL0awRjZZfIzDSGydf9pkpd9Hy0C1KtYUia62G3yI/faESNPguez7j2TJ3/KNtC547vxJllO/Jl6vsGUrjCFyK7iwIaqvfuGewpb9pvRuAGOuLzeDXiEazocFUjIlcxcS9S1jMtRJz4FXbe+hJuMRdqa1+bbS6P5q9R3d4fhizH+lL8PDJSwuef/T4OHcceYRb7n2QO448AjDnvbajP8VwX4KrRnq5flc/iZjDhx84xtHj46v+m43Njc3gV0Gj5eKRvezo8XF+9k++s8EjNJrhSmB7XKoRylIIEJSpCcoD58tebQbfqEa6AMm4Q7Hik4wJvsLBnf1A0PJwvoOlaYGy26+vvdfuOPIIZc/f8KbeRudhAt8CjYQcmNNx6dxUcU4rvvseOrlkUSdj4wjCI6tfRhw0WAFXlAu5uQ4bZWGDcgVKFb/mskrU1Ulv5GBppXmH2R2NZpjAL0F967zJfIWxyQIPn7wIBPEtxxFcCbr6eKp84uiL3LhnkOfPN+1MaHQRPqHIhxfzeFgHPqKpi0qDGX6mN4aqNq370op4m93RaIbF4JfgvodOUq56TOTKC/qQ+hA2XQ5scZ6vVD2fe796nFw31dztYtayn7MAriuL7lMIVq0m4w67+pPsH+5b1MEy350FC8V7M9gd5+cRLD+wPtgMfglOTebJFqtNe07C7CwtWrD67Llsw+Sa0XmsRdHM6P8/6oebjLnEwhXMCdeh7Pm4jqB+ULLgwEgfqspUobJkSYBWmmMcPjjKPdCx9dXXq9GJsRAT+CXYO9TDmRaWfc/H4u9bi96EQ7GqbO9NMNyX5HhYR74vGWMiVwobkENVgyt/qyGUVsW7k+2OK20CbqyedvZk/WPgJ4BxVb2hXcdpB/VJ1UwyZitPjUWR8N/+lEO+7DFVqNCTcIk5wmS+giOCalAN0veVC7kicddtOYTSyeLdCpYE3jjaGYP/r8Db2rj/thDdTo5niwym45Q7sWC70TE4QDLmUPWUXMkj4Qp/92tv4ePvfj0zZQ9FiTmC60ittMFMydtSK0ZbySMY7aFtM3hVfUhErmzX/ttFdDtZ9ZQTl7KUF2mqYGxN4qFjCiAZC+oPiUCl6jNZqHLo3329Vo/GEfAFUjGH4b4kmVSMqUKlaUOQ+x46yfPnp6l4SiLm1GrKb+aLgTXZ3jg2PAYvIncCdwLs27dvg0cT3E66AqcmCxZHN2oIQbngTDrOaCbFifEsU/kKvioiUPWCqpGiylS+AjLXA1/1fM5MFXCzwnBPnDuOPNJ0XcV0MSgbXCh7vHwxt+kTkp2eBO5m2trwI5zBf6nVGPxGNPyYv4hpbDIf1PA2cTdChKBbkwiMTxcpeUpvGGMvVX3KXmCRjRqWx10HR4RKuD1a8BTdGSKwoz85pyl3T9yh4ivnpoJmIo4j+BqEd3YOpKwJh9EUa/jRgKPHx7n3q8cDt4OAiAQfLlN2ow5XIBV3EYEzl4soSirm0JNwmciVGUgFC5Uqntb8kpWwa1PMkdr7SYmEP7BPTheqDPelao6Sly7muXq0L7BUSpi2laD8sCUkjZWyJQR+/iw96kV59nI+KPla+8cwgiJhuwdTlKo+E7ky/ekY49OBuAMM9yXpT8cpVjwmw5Z89USNXqLEatwV4q7DgZE+jp+bxgmFOyLqI1CoeGGjkNm7gYTrWELSWDFtc9GIyGeBh4FrRWRMRN7XrmMtxnxXzHi2yCeOvki56q3JIheju+hNBLNzX2H/cB+/9JZXc+X2PkpesGhp10C61iS74imOBAIdm6fySrDoKeYEgj/clwQCwW5Ug+bAcC8VT8mkYvgoVd/H95X+dMwSksaKaaeL5o527Xs5NFpkUfV9ssXqEq80thKpmIMqlD1lIC5zGlbfTeMG1aWqTzJcqRpzHSRcvRq06Atm4K8ezTCRKxFzAy98JhWr3RXU16D5rR8/CATv16o3TTl00Vy5vc8SksaK6foQTaNFFknXCT6cMafWocfY2ogETbF9Xyk3aOzRyOrnOsJAT5zpQpWqp8RcB6dBYjQKEY5N5tk/3Mc/+YFtPHzyUkNHiQm5sZZ0vcA3qrQ30BPnQq5s4r4FiHqnLpU7VzTMxQR14ufTyOr3ztft4v7HT5NJwcWZMn7Yr6+/Nz4nrFK/ErVR6WkTdaNddL3AN5p55cvBbbHR3SQc8AlCI0u10iuGNdoH03GuHs3Uti8lyDfuGWw5rGJFt4z1pusFfv7Mqy8ZY6pg8fetgIfQm3ApVjwcaHpRFyARdleaKlZ504GgL2orgrycOjFWdMtYb7pW4OfPvHb2Jzh2ZopTK6gMaWxOrh7p5UNvv473feo7uCILmnFA4JqBoOl1wnXIpGJ85ZlzPHzyEo+/MokAOwdSiMiqBdmKbhnrTVcK/PyZ13Pnpnn4ZGWjh2WsMyLC4YOjXD3Sx8uX8nh+sAJ1uC9JzBXGJgvsH+5DZNbjOF0o8/x4jiu39wRlCAgWOO0ahEwqvipBts5LxnrTlR2d6m+Fc6Uql/Im7luRiVwJgA+9/TpGMyn2beth/3AvMVeoeMqB4d4FVQ7PT5dq752E6yAE7pqJbLCv1QjyZui8ZHQXXSnwpybzpOMu2WLQQ9XoTpZqmlUOXVKHD45yz+3XM5pJzWmP92tvO7hQcH2fHZlgUdJwXxIfRVUpVb1VC3KzcVj83WgXXReiOXp8nOlChbNTBXx/ceeEsflp5o5xHUi4s5eAZsnQ+dbH+jh9tGL1fLaIqDCaSa3a1rjZm3cYm4uuEvgo9u46YH06upuYI7VKjfUEjhiHoZ44+4f7ltzPfMGN3kORrTbmis20jU1L1wj80ePj/MJnHqNQMWXvBmJO84bYcVfY3psgW6xSqvpI2FBjtD+16oYSVrvc6Ca6QuCPHh/ng/c/aeK+yXEFrgyLbt1z+/U8NXaZP/jbF2r+9YQrDKTjvPeNr1qw1B/WTpQtjGJ0C10h8Pc9dJJLufJGD8NYBQlXEJkb5z58cLS2UnS+cN/dYB8myoYxl00v8EePj/P4K5PY3L1zcWTxWjBDPTEyqUTDOLfNpg1j5WxqgY8SYkvZ5YyNIxkmPKu+ki97FMrenMRoJulycOeAxbkNow1sSoH3fOXEeJbf/uvvcTlfxrduTOtKzAlWiVY8JRlzqPrB93RMyJWDe6nRTJLehMtM2WtLnNwwjKXpeIFXVV6+mOepscs8NTbFk6cuc+zM9IIViMbaIsDugSQ//Y/21RKaABdyJcqe0hN3+cUf3c/dt12zrP2aoBvG+tFWgReRtwEfA1zgk6r6u0u95uxUgSdPTfHU2GWeHLvM02NTTDfpvhRzhLgrtRohr1zK28KmOhzgA7ddzd23XTOn6cRyZ8+NEpqGYXQ+om0Kb4iICzwP/BgwBnwHuENVv9fsNb27r9GR9/5+48eSLjfsGuCmvYPcuGeQG/cM8ML5LB/5798j7grpuMuz56ZxwtDBVmC4N85IJkWuVLWQh2FsUUTkMVU91Oixds7gfwB4QVVPhoP4c+CdQFOBr/qR39nhNbv6QzEf4MY9gxwY7sVx5qZT927rQURqM9PeRIyehEsy5vL9S5uvBOtwb5yP/tRNJtKGYawJ7RT43cCput/HgB+c/yQRuRO4E2Bkz36+9Iu3cO3ODHG3tTpo89uhffiBY8RcYaQvzkRufatIJhxhx0CKQsWjXPVJuMLVO/ptZm0YxobQToFv5F5cEDtR1SPAEYBDhw7pDbsHVnzA+mXm5arPwR1JRISXLmSJwviOQE88aLbthT04FwvoDPclSMddC4EYhrHpaKfAjwF7637fA5xp4/EAWxhjGIYR0c568N8BrhaR/SKSAN4NPNDG4xmGYRh1tG0Gr6pVEflXwN8Q2CT/WFWPtet4hmEYxlza6oNX1S8DX27nMQzDMIzGdGXLPsMwDMME3jAMo2sxgTcMw+hSTOANwzC6lLbVolkJIjIBfH/e5mHgwgYMpxOxczGLnYsAOw+zbNVz8SpVHWn0QEcJfCNE5NFmhXS2GnYuZrFzEWDnYRY7FwuxEI1hGEaXYgJvGIbRpWwGgT+y0QPoIOxczGLnIsDOwyx2LubR8TF4wzAMY2Vshhm8YRiGsQJM4A3DMLqUjhZ4EXmbiDwnIi+IyIc2ejxrjYjsFZG/FZFnReSYiPxSuH2biHxdRE6E34fqXvPr4fl4TkT+cd32m0Xk6fCxj4tIo4YrHY2IuCLyXRH5Uvj7Vj0PgyJyv4gcD98bb9rC5+KXw8/GMyLyWRFJbdVzsSJUtSO/CEoMvwgcABLAk8BrNnpca/w3XgG8Ifw5Q9Ck/DXA/w18KNz+IeDe8OfXhOchCewPz48bPvYPwJsIOml9BXj7Rv99KzgfvwL8GfCl8Peteh7+BPi58OcEMLgVzwVB28+XgHT4++eAf74Vz8VKvzp5Bl9r2q2qZSBq2t01qOpZVX08/DkLPEvwpn4nwYec8Pv/HP78TuDPVbWkqi8BLwA/ICJXAP2q+rAG7+ZP1b1mUyAie4AfBz5Zt3krnod+4FbgjwBUtayql9mC5yIkBqRFJAb0EHSF26rnYtl0ssA3atq9e4PG0nZE5Erg9cDfAztU9SwEFwEg6kHY7JzsDn+ev30z8R+Afw34ddu24nk4AEwA/yUMV31SRHrZgudCVU8DHwVeAc4CU6r6NbbguVgpnSzwLTXt7gZEpA/4b8AHVHV6sac22KaLbN8UiMhPAOOq+lirL2mwbdOfh5AY8AbgP6nq64EZgjBEM7r2XISx9XcShFt2Ab0i8p7FXtJgW1eci5XSyQK/IU271xsRiROI+2dU9Qvh5vPhbSXh9/Fwe7NzMhb+PH/7ZuHNwO0i8jJBKO4tIvJptt55gOBvGFPVvw9/v59A8LfiubgNeElVJ1S1AnwB+CG25rlYEZ0s8F3ftDvM5P8R8Kyq/vu6hx4A/ln48z8D/qpu+7tFJCki+4GrgX8Ib1OzIvLGcJ//a91rOh5V/XVV3aOqVxL8Pz+oqu9hi50HAFU9B5wSkWvDTW8FvscWPBcEoZk3ikhP+De8lSBPtRXPxcrY6CzvYl/AOwicJS8C/2ajx9OGv+8WglvFp4Anwq93ANuBbwInwu/b6l7zb8Lz8Rx1TgDgEPBM+Nj/S7hKebN9AYeZddFsyfMA3AQ8Gr4v/hIY2sLn4v8Ajod/x58SOGS25LlYyZeVKjAMw+hSOjlEYxiGYawCE3jDMIwuxQTeMAyjSzGBNwzD6FJM4A3DMLoUE3jDMIwuxQTe2BKEi1++ISJPiMjPiMhvrGJfubUcm2G0i9hGD8Aw1onXA3FVvQlqIv1/buiIDKPN2Aze2LSISK+I/LWIPBk2hPgZCZrEHBeRb4eNHb4kIqPAp4Gbwhn85wlK0D4hIp9ZZP9/KSKPhQ0n7pz32P8jIo+LyDdFZCTcdpOIPCIiT4nIF0VkSESuE5F/qHvdlSLyVPjzzSLyrfAYfxPVVzGMtcIE3tjMvA04o6qvU9UbgK8C/xn4SeCHgZ0AqjoO/Bzwd6p6k6r+FFAIf/6ni+z/Z1X1ZoJl7neLyPZwey/wuKq+AfgW8JFw+6eAX1PVG4GngY+o6rNAQkQOhM/5GeBzYZG5PwDeFR7jj4HfWfUZMYw6TOCNzczTwG0icq+I/DBBWdmXVPWEBjU4Pr3K/d8tIk8CjxBUKbw63O4DfxH+/GngFhEZAAZV9Vvh9j8haNwBQSeinw5//pnwtdcCNwBfF5EngN9kbsVDw1g1FoM3Ni2q+ryI3ExQoO3/Ar7GGtX5FpHDBOVq36SqeRE5CqSaDWWJ3f0F8HkR+UIwbD0hIq8Fjqnqm9ZivIbRCJvBG5sWEdkF5FX10wSdf34I2C8iV4VPuWORl1fCMEkzBoDJUNwPAm+se8wB3hX+/E+Ab6vqFDAZ3kkAvJcgfIOqvgh4wG8xO/N/DhgRkTeFf0tcRK5f8o82jGVgM3hjM/Na4PdExAcqwC8Aw8Bfi8gF4NsEYZBGHAGeEpHHm8Thvwr8b2FC9DmCME3EDHC9iDwGTBGEXSCoTf6HItIDnAT+Rd1r/gL4PYIwEqpaFpF3AR8PwzsxgraFx5bx9xvGoli5YKNrCcMsH1TVn9jgoRjGhmAhGsMwjC7FZvDGlia0Pn6zwUNvVdWL6z0ew1hLTOANwzC6FAvRGIZhdCkm8IZhGF2KCbxhGEaXYgJvGIbRpfz/iQx77nEb/7UAAAAASUVORK5CYII=\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)"
   ]
  },
  {
   "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": 25,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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": 25,
     "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": [
    "\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": 26,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.00046769430149007363"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[['long']]\n",
    "Y = df['price']\n",
    "lm = LinearRegression()\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 <code>'sqft_living'</code> then calculate the R^2. Take a screenshot of your code and the value of the R^2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.49285321790379316"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[['sqft_living']]\n",
    "Y = df['price']\n",
    "lm = LinearRegression()\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 <code>'price'</code> using the list of features:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "features =[\"floors\", \"waterfront\",\"lat\" ,\"bedrooms\" ,\"sqft_basement\" ,\"view\" ,\"bathrooms\",\"sqft_living15\",\"sqft_above\",\"grade\",\"sqft_living\"]     "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Then calculate the R^2. Take a screenshot of your code."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6576485656702705"
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lm2=LinearRegression()\n",
    "lm2.fit(df[features], df['price'])\n",
    "lm2.score(df[features],df['price'])"
   ]
  },
  {
   "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": 32,
   "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 to predict the 'price', fit the object using the features in the list <code>features</code>, and calculate the R^2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jupyterlab/conda/envs/python/lib/python3.6/site-packages/sklearn/preprocessing/data.py:625: DataConversionWarning: Data with input dtype int64, float64 were all converted to float64 by StandardScaler.\n",
      "  return self.partial_fit(X, y)\n",
      "/home/jupyterlab/conda/envs/python/lib/python3.6/site-packages/sklearn/base.py:465: 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",
      "/home/jupyterlab/conda/envs/python/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.7513408481278352"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pipe= Pipeline(Input)\n",
    "pipe.fit(df[features], df['price'])\n",
    "pipe.score(df[features], df['price'])\n"
   ]
  },
  {
   "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": 37,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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 sets:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "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, set the regularization parameter to 0.1, and calculate the R^2 using the test data. \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.linear_model import Ridge"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6478759163939115"
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "RidgeModel= Ridge(alpha=0.1)\n",
    "RidgeModel.fit(x_train,y_train)\n",
    "RidgeModel.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, set the regularisation parameter to 0.1, and calculate the R^2 utilising the test data provided. Take a screenshot of your code and the R^2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.7002744283956152"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pr = PolynomialFeatures(degree=2)\n",
    "x_train_pr = pr.fit_transform(x_train)\n",
    "x_test_pr = pr.fit_transform(x_test)\n",
    "RidgeModelPr=Ridge(alpha=0.1)\n",
    "RidgeModelPr.fit(x_train_pr, y_train)\n",
    "RidgeModelPr.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, and 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",
   "language": "python",
   "name": "conda-env-python-py"
  },
  "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.13"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}