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": 48,
   "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": 30,
   "metadata": {
    "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": 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": 19,
   "metadata": {
    "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": 19,
     "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": 20,
   "metadata": {
    "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": 20,
     "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": 31,
   "metadata": {
    "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": 31,
     "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": 37,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "13"
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.bedrooms.isnull().sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "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": 38,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean=df['bedrooms'].mean()\n",
    "df['bedrooms'].replace(np.nan,mean, inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 39,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.bedrooms.isnull().sum()"
   ]
  },
  {
   "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": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean=df['bathrooms'].mean()\n",
    "df['bathrooms'].replace(np.nan,mean, inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {
    "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": 45,
   "metadata": {
    "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": 45,
     "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": 121,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x2cb2a5086d0>"
      ]
     },
     "execution_count": 121,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAERCAYAAABxZrw0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAVG0lEQVR4nO3dfZDdVX3H8c9nd8EkPIxlWRnyIBGDIII8rQ8tVnlYnNXWONNRkamyWmraWpPQFh21Th1nrNOZdmqTVC07lrIZFQsUHHRgNVtAHgR0EwIhBMctBsxGYV2Rp8TA7n77x70bdslmcwk5v98v575fMzu5v/u7955vNnc/OXvu+Z3jiBAAID8tZRcAAEiDgAeATBHwAJApAh4AMkXAA0CmCHgAyFTlAt72FbYft/1Ag4//gO0HbW+2/a3U9QHAwcJVmwdv++2SnpG0NiJO2cdjT5B0taTzIuIJ26+KiMeLqBMAqq5yPfiIuE3Sb6beZ/u1tvttr7d9u+2T6qc+JukrEfFE/bmEOwDUVS7g96JX0vKIOEvSZZK+Wr//dZJeZ/tO23fb7i6tQgComLayC9gX24dL+gNJ19ievPsV9T/bJJ0g6RxJCyXdbvuUiPht0XUCQNVUPuBV+y3jtxFx+gzntkm6OyKel/Rz2z9VLfB/UmSBAFBFlR+iiYinVAvv90uSa06rn/6OpHPr9x+t2pDNw6UUCgAVU7mAt32VpLsknWh7m+1LJP2ppEts3ydps6T31h/+fUmjth+UdIukT0bEaBl1A0DVVG6aJADgwEjag7f9N/ULkB6wfZXtOSnbAwC8IFkP3vYCSXdIOjkidtq+WtKNEXHl3p5z9NFHx+LFi5PUAwA5Wr9+/a8jomOmc6ln0bRJmmv7eUnzJG2f7cGLFy/W4OBg4pIAIB+2H9nbuWRDNBExLOlfJD0q6ZeSnoyIH8xQ3DLbg7YHR0ZGUpUDAE0nWcDb/j3VZru8RtJ8SYfZ/tCLHxcRvRHRGRGdHR0z/pYBANgPKT9k7ZL084gYqV+IdJ1qV6QCAAqQMuAflfRW2/NcW2PgfElbErYHAJgi5Rj8PZKulbRB0qZ6W72p2sMLRkdHtWLFCo2Ocs0X0MySzoOPiM9HxEkRcUpEfDgidqVsDzV9fX3atGmT1q5dW3YpAEpUuaUK8PKMjo6qv79fEaH+/n568UATI+Az09fXp4mJCUnS+Pg4vXigiRHwmRkYGNDY2JgkaWxsTOvWrSu5IgBlIeAz09XVpba22gXKbW1tuuCCC0quCEBZCPjM9PT0qKWl9s/a2tqqiy++uOSKAJSFgM9Me3u7uru7ZVvd3d1qb28vuyQAJTkYtuzDS9TT06OtW7fSeweaHAGfofb2dq1evbrsMgCUjCEaAMgUAQ8AmSLgASBTBDwAZIqAB4BMEfAAkCkCHgAyRcADKBQb0hSHgAdQKDakKU6ygLd9ou2NU76esn1pqvYAVB8b0hQr5Z6sP42I0yPidElnSdoh6fpU7QGoPjakKVZRQzTnS/q/iHikoPaaGmOcqCo2pClWUQH/QUlXzXTC9jLbg7YHR0ZGCionb4xxoqrYkKZYyQPe9qGSlkq6ZqbzEdEbEZ0R0dnR0ZG6nOwxxokqY0OaYhXRg3+XpA0R8VgBbTU9xjhRZWxIU6wiAv4i7WV4BgceY5youp6eHp166qn03guQNOBtz5N0gaTrUraDFzDGiaqb3JCG3nt6SQM+InZERHtEPJmyHbyAMU4Ak7iSNTOMcQKYxJ6sGWLTbQASAZ8lNt0GIDFEAwDZIuABIFMEPABkioAHgEwR8ACQKQIeADJFwANApgh4AMgUAQ8AmSLgASBTBDwAZIqAB4BMEfAAkCkCHgAylXrLvlfavtb2Q7a32P79lO0BAF6Qej34VZL6I+J9tg+VNC9xewCAumQBb/tISW+X9BFJiojnJD2Xqj0AwHQph2iOlzQi6b9s32v767YPe/GDbC+zPWh7cGRkJGE5ANBcUgZ8m6QzJX0tIs6Q9KykT7/4QRHRGxGdEdHZ0dGRsBwAaC4pA36bpG0RcU/9+FrVAh8AUIBkAR8Rv5L0C9sn1u86X9KDqdoDAEyXehbNcknfrM+geVjSRxO3BwCoSxrwEbFRUmfKNgAAM+NKVgDIFAGfodHRUa1YsUKjo6NllwKgRAR8hvr6+rRp0yatXbu27FIAlIiAz8zo6Kj6+/sVEerv76cXDzQxAj4zfX19mpiYkCSNj4/TiweaGAGfmYGBAY2NjUmSxsbGtG7dupIrAlAWAj4zXV1damurzX5ta2vTBRdcUHJFAMpCwGemp6dHLS21f9aWlhZdfPHFJVcEoCwEfGba29s1f/58SdL8+fPV3t5eckXAdEzjLQ4Bn5nR0VENDw9LkrZv384PESqHabzFIeAz09fXp4iQJE1MTPBDhEphGm+xCPjMMIsGVcY03mIR8JlhFg2qjA5IsQj4zEydRdPa2sosGlQKHZBiEfCZaW9vV3d3t2yru7ubWTSolJ6ent1DNBMTE3RAEku94QdK0NPTo61bt/LDAzQ5evAZam9v1+rVq+m9o3L6+vpkW5Jkmw9ZE0sa8La32t5ke6PtwZRtAai+gYEBjY+PS6rNouFD1rSK6MGfGxGnRwRb9wFNjg9Zi8UQDYDCMMurWKkDPiT9wPZ628tmeoDtZbYHbQ+OjIwkLgdAmZjlVazUs2jOjojttl8laZ3thyLitqkPiIheSb2S1NnZGYnrAVAyZnkVJ2kPPiK21/98XNL1kt6csj0A1ccsr+IkC3jbh9k+YvK2pHdKeiBVewCA6VIO0Rwj6fr6nNc2Sd+KiP6E7QEApkjWg4+IhyPitPrXGyLiH1O1BeDgwYYfxWGaJIBCseFHcQh4AIVhw49iEfAACsOGH8Ui4DPEGCeqig0/ikXAZ4gxTlRVV1fXtNUkWYsmLQI+M4xxosqWLl26e1P4iNB73vOekivKGwGfGcY4UWU33HDDtB78d7/73ZIryhsBnxnGOFFlAwMD03rwvD/TIuAzw3rbqDLen8Ui4DPDetuoMt6fxSLgM8N626gy3p/FIuAztHTpUs2bN48ZCqiknp4enXrqqfTeC0DAZ+iGG27Qjh07mKGASmI9+OI0HPC2j7PdVb89d3Ktd1QL8+ABTGoo4G1/TNK1ki6v37VQ0ndSFYX9xzx4AJMa7cH/taSzJT0lSRHxM0mvSlUU9h/z4AFMajTgd0XEc5MHttsksUF2BTHPGMCkRgP+h7Y/K2mu7QskXSOpoU/wbLfavtf29/a3SDSup6dn9xDNxMQEMxWAJtZowH9a0oikTZL+QtKNkj7X4HNXStry0ksDALwcjQb8XElXRMT7I+J9kq6o3zcr2wsl/ZGkr+9/iXgp+vr6pi3mxIesQPNqNOD/V9MDfa6kgQae92+SPiVpYm8PsL3M9qDtwZGRkQbLwd4MDAxofHxcUm0WDR+yAs2r0YCfExHPTB7Ub8+b7Qm2/1jS4xGxfrbHRURvRHRGRGdHR0eD5WBvurq6dq/10dLSwoesQBNrNOCftX3m5IHtsyTt3Mdzzpa01PZWSd+WdJ7tb+xXlWgYH7ICmNTW4OMulXSN7e3142MlXTjbEyLiM5I+I0m2z5F0WUR8aD/rRIOeeOKJPY65JByStGbNGg0NDZVdhoaHhyVJCxYsKLWOJUuWaPny5aXWkFpDPfiI+ImkkyT9laSPS3r9voZeUI4vfvGLsx4DZdu5c6d27tzXAAAOhFl78LbPi4ibbf/Ji06dYFsRcV0jjUTErZJu3b8S8VJs3bp11mM0r6r0VleuXClJWrVqVcmV5G9fQzTvkHSzpJnWnQ1JDQU8irNw4UJt27Zt9/GiRYtKrAZAmWYN+Ij4vO0WSTdFxNUF1YSXYdGiRdMCfuHChSVWA6BM+xyDj4gJSZ8ooBYcAPfcc8+sxwCaR6PTJNfZvsz2IttHTX4lrQz7ZXLH+r0dA2gejU6T/DPVxtw//qL7jz+w5eDlamlp2X0l6+QxgObU6E//yZK+Iuk+SRslrZH0hlRFYf91dXXNegygeTQa8H2SXi9ptWrh/vr6faiYZcuWzXoMoHk0OkRzYkScNuX4Ftv3pSgIAHBgNNqDv9f2WycPbL9F0p1pSsLLcfnll0877u3tLakSAGVrNODfIulHtrfWFw+7S9I7bG+yfX+y6vCSDQxMX8WZ5YKB5tXoEE130ipwwEyuJLm3YwDNo6GAj4hHUhcCADiwmCQNAJki4DNz1FFHzXoMoHkQ8Jl58sknZz0G0DwI+MxMXaZgpmMAzYOAB4BMJQt423Ns/9j2fbY32/5CqrYAAHtqdB78/tgl6byIeMb2IZLusH1TRNydsE0AQF2ygI/aQuTP1A8PqX+xODkAFCTpGLztVtsbJT0uaV1E7LG9kO1ltgdtD46MjKQsBwCaStKAj4jxiDhd0kJJb7Z9ygyP6Y2Izojo7OjoSFkOADSVQmbRRMRvJd0q1rQBgMKknEXTYfuV9dtzJXVJeihVewCA6VLOojlWUp/tVtX+I7k6Ir6XsD0AwBQpZ9HcL+mMVK8PAJgdV7ICQKYIeADIFAEPAJki4AEgUwQ8AGSKgAeATBHwAJApAh4AMkXAA0CmCHgAyBQBDwCZIuABIFMEPABkKuVywQDq1qxZo6GhobLLqITJ78PKlStLrqQalixZouXLlyd5bQIeKMDQ0JB+tvlevfrw8bJLKd2hz9cGDnY9MlhyJeV79JnWpK9PwAMFefXh4/rsmU+VXQYq5Esbjkz6+ozBA0CmUu7Jusj2Lba32N5smwE3AChQyiGaMUl/FxEbbB8hab3tdRHxYMI2AQB1yXrwEfHLiNhQv/20pC2SFqRqDwAwXSFj8LYXq7YB9z0znFtme9D24MjISBHlAEBTSB7wtg+X9D+SLo2IPaYQRERvRHRGRGdHR0fqcgCgaSQNeNuHqBbu34yI61K2BQCYLuUsGkv6T0lbIuJfU7UDAJhZyh782ZI+LOk82xvrX+9O2B4AYIpk0yQj4g5JTvX6AIDZcSUrAGSKgAeATBHwAJApAh4AMkXAA0CmCHgAyBQBDwCZIuABIFMEPABkij1ZgQIMDw/r2adbk+/BiYPLI0+36rDh4WSvTw8eADJFDx4owIIFC7Rr7Jf67Jl7bImAJvalDUfqFQvSbXRHDx4AMkXAA0CmCHgAyBQBDwCZ4kPWA2zNmjUaGhoqu4xpVq5cWVrbS5Ys0fLly0trH2hmKfdkvcL247YfSNUGAGDvUvbgr5T075LWJmyjcsrurZ5zzjl73Ldq1ariCwFQumQ9+Ii4TdJvUr0+ZjZnzpxpx3Pnzi2pEgBlK/1DVtvLbA/aHhwZGSm7nINef3//tOObbrqppEoAlK30gI+I3ojojIjOjo6OssvJCr13oLkxiyZDp512miTG3oFmV3oPHgCQRsppkldJukvSiba32b4kVVsAgD0lG6KJiItSvTZwMHr0GdaDl6THdtT6lcfMmyi5kvI9+kyrTkj4+ozBAwVYsmRJ2SVUxnP1K71fcRzfkxOU9r1BwAMFKPsCuCqZXDqDSQDpZRPwVVwDpiyT34cy16CpEtbDQbPKJuCHhoa08YEtGp93VNmllK7luZAkrX/4sZIrKV/rDi6mRvPKJuAlaXzeUdp50rvLLgMVMvehG8suASgN8+ABIFMEPABkKpshmuHhYbXueJJfyTFN645RDQ+PlV0GUAp68ACQqWx68AsWLNCvdrXxISummfvQjVqw4JiyywBKQQ8eADKVTQ9eqs15ZgxeavndU5KkiTmse1KbB08PXqrOxYBVuRCvGS6AyybgWevjBUNDT0uSlhxPsEnH8N6oGDaiKY4jouwaduvs7IzBwcGyyzjosdYH0Dxsr4+IzpnOMQYPAJki4AEgUwQ8AGQq6YestrslrZLUKunrEfFPKdurgirMVKjKLAWpOWYqAFWVck/WVklfkfQuSSdLusj2yanawwvmzp3LTAUASXvwb5Y0FBEPS5Ltb0t6r6QHE7ZZOnqrAKoi5Rj8Akm/mHK8rX7fNLaX2R60PTgyMpKwHABoLikD3jPct8ek+4jojYjOiOjs6OhIWA4ANJeUAb9N0qIpxwslbU/YHgBgipQB/xNJJ9h+je1DJX1Q0g0J2wMATJHsQ9aIGLP9CUnfV22a5BURsTlVewCA6ZLOg4+IGyWxvCMAlIArWQEgUwQ8AGSqUssF2x6R9EjZdWTiaEm/LrsIYC94fx44x0XEjHPMKxXwOHBsD+5tjWigbLw/i8EQDQBkioAHgEwR8PnqLbsAYBa8PwvAGDwAZIoePABkioAHgEwR8Bmy3W37p7aHbH+67HqASbavsP247QfKrqUZEPCZYatEVNyVkrrLLqJZEPD52b1VYkQ8J2lyq0SgdBFxm6TflF1HsyDg89PQVokA8kfA56ehrRIB5I+Azw9bJQKQRMDniK0SAUgi4LMTEWOSJrdK3CLparZKRFXYvkrSXZJOtL3N9iVl15QzlioAgEzRgweATBHwAJApAh4AMkXAA0CmCHgAyBQBj6Zk+1Lb8/bjeSfZ3mj7XtuvLasOoBFMk0RTsr1VUmdE/PolPKdV0iclzY2Iz7/onFX7eZpIXQfQKHrwOKjZ/pTtFfXbX7Z9c/32+ba/Yftrtgdtb7b9hfq5FZLmS7rF9i31+95p+y7bG2xfY/vw+v1bbf+D7TskXSjpUkl/bvsW24ttb7H9VUkbJC2y/c+2H7C9yfaF9dc4x/attq+1/ZDtb7pmjzqAA4mAx8HuNkl/WL/dKelw24dIepuk2yX9fUR0SnqjpHfYfmNErFZtfZ5zI+Jc20dL+pykrog4U9KgpL+d0sbvIuJtEfEtSf8h6csRcW793ImS1kbEGfX2T5d0mqQuSf9s+9j6485Q7T+HkyUdL+nsF9dxgL8vAAGPg956SWfZPkLSLtUug+9ULfRvl/QB2xsk3SvpDaoF7Iu9tX7/nbY3SuqRdNyU8/89S/uPRMTd9dtvk3RVRIxHxGOSfijpTfVzP46IbfUhnI2SFr/kvynwErWVXQDwckTE8/Vx7I9K+pGk+yWdK+m1knZKukzSmyLiCdtXSpozw8tY0rqIuGgvzTw7SwlTz820VPOkXVNuj4ufPRSAHjxycJtqQX6bar32v1Stl3ykagH8pO1jVNvGcNLTko6o375b0tm2l0iS7Xm2X7efdVxou9V2h6S3S/rxPp4ztQ7ggCLgkYPbJR0r6a760MjvJN0eEfepNjSzWdIVku6c8pxeSTfZviUiRiR9RNJVtu9XLfBP2o86rlftN4j7JN0s6VMR8at9PGd3HfvRHjArpkkCQKbowQNApgh4AMgUAQ8AmSLgASBTBDwAZIqAB4BMEfAAkKn/B6beAVavjmuTAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.boxplot(x=df[\"waterfront\"], y=df[\"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": 120,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 8081250.0)"
      ]
     },
     "execution_count": 120,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAESCAYAAAD38s6aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO29fXgc93Xf+zkz+wosAIIEQFJ8EUmJEmXZsi3zcexGlRnb6bWT1L63j/KiNr1N00RqmlhJep3G6U3s5yq390Y3aROr103FOGnjOrVrq06iprGTODItu1dKItGWLVmUKJGSSFEkQBAvu9jXmTn3j5lZ7AK7wILE4mVxPnogLGZnZ3473P3O+Z1zfueIqmIYhmH0Hs56D8AwDMPoDibwhmEYPYoJvGEYRo9iAm8YhtGjmMAbhmH0KCbwhmEYPcqGE3gR+T0RGReRZzrc/4dE5Dsi8qyI/Oduj88wDGOzIBstD15E7gQKwKdU9Y3L7HsY+BzwblWdEpExVR1fi3EahmFsdDacBa+qjwFXGreJyA0i8iUReUpEviYiR6KnfhL4hKpORa81cTcMw4jYcALfhuPAh1T1bcCHgX8Xbb8JuElE/oeIPCEi71u3ERqGYWwwEus9gOUQkRzwt4DPi0i8OR39TgCHgWPAXuBrIvJGVZ1e63EahmFsNDa8wBPOMqZV9S0tnjsPPKGqNeCsiDxPKPh/s5YDNAzD2IhseBeNqs4SivcPAkjIm6On/wj4nmj7CKHL5sy6DNQwDGODseEEXkQ+AzwO3Cwi50XknwD/APgnIvI08CzwwWj3PwMmReQ7wFeAX1DVyfUYt2EYxkZjw6VJGoZhGKtDVy14Efn5aAHSMyLyGRHJdPN8hmEYxjxds+BFZA/wdeANqloSkc8Bf6qq/7Hda0ZGRvTAgQNdGY9hGEYv8tRTT11W1dFWz3U7iyYBZEWkBvQBF5ba+cCBAzz55JNdHpJhGEbvICKvtHuuay4aVX0N+A3gVeB1YEZV/7xb5zMMwzCa6ZrAi8gwYbbLQeA6oF9EfrTFfveIyJMi8uTExES3hmMYhrHl6GaQ9b3AWVWdiBYifYFwRWoTqnpcVY+q6tHR0ZZuJMMwDOMq6KbAvwq8Q0T6JKwx8B7guS6ezzAMw2igmz74vwIeBk4C347Odbxb5zMMwzCa6WoWjap+DPhYN89hGIZhtGYzFBszDMPYsJw4Nc5Dj53h3FSRfcN93HvnIY4dGVvvYQEbsBaNYRjGZuHEqXE++sizjOfLbMsmGc+X+egjz3Li1MboPWQCbxiGcZU89NgZkq7Ql0ogEv5OusJDj22MorYm8IZhGFfJuaki2aTbtC2bdDk/VVynETVjAm8YhnGV7Bvuo1Tzm7aVaj57h/vWaUTNmMAbhmFcJffeeYiarxSrHqrh75qv3HvnofUeGmACbxiGcdUcOzLG/R+4lbGBDDOlGmMDGe7/wK0bJovG0iQNwzCugWNHxjaMoC/ELHjDMIwexQTeMAyjRzGBNwzD6FFM4A3DMHoUE3jDMIwexQTeMAyjRzGBNwzD6FFM4A3DMHoUE3jDMIwexVayGoaxadnIzTY2Al2z4EXkZhH5ZsPPrIj8XLfOZxjG1mKjN9vYCHSz6fbzqvoWVX0L8DagCPxht85nGMbWYqM329gIrJUP/j3AS6r6yhqdzzCMHmejN9vYCKyVwP8I8Jk1OpdhGFuAjd5sYyPQdYEXkRTwAeDzbZ6/R0SeFJEnJyYmuj0cwzB6hI3ebGMhJ06Nc/fxJ7jjgUe5+/gTaxIrWAsL/v3ASVW91OpJVT2uqkdV9ejo6OgaDMcwjF5gozfbaGS9AsJrkSZ5N+aeMQyjC2zkZhuNNAaEAfpSCYpVj4ceO9PV8XfVgheRPuB7gS908zyGYRgbmfUKCHdV4FW1qKo7VHWmm+cxDMPYyKxXQNhKFRiGYXSZ9QoIm8AbhmF0mfUKCFstGsMwjDVgPQLCZsEbhmH0KCbwhmEYPYoJvGEYRo9iAm8YhtGjmMAbhmH0KCbwhmEYPYoJvGEYRo9iAm8YhtGjmMAbhmH0KCbwhmEYPYoJvGEYRo9iAm8YhtGjmMAbhmH0KCbwhmEYPYoJvGEYRo/S7Z6s20TkYRE5JSLPicg7u3k+wzAMY55uN/z4OPAlVb1LRFJAdxsQGoZhGHW6JvAiMgjcCfwYgKpWgWq3zmcYhmE0000XzSFgAvgPIvINEfmkiPQv3ElE7hGRJ0XkyYmJiS4OxzAMY2vRTYFPALcDv62qbwXmgI8s3ElVj6vqUVU9Ojo62sXhGIZhbC26KfDngfOq+lfR3w8TCr5hGIaxBnRN4FX1InBORG6ONr0H+E63zmcYhmE00+0smg8BfxBl0JwB/nGXz2cYhmFEdFXgVfWbwNFunsMwDMNoja1kNQzD6FFM4A3DMHoUE3jDMIwepdtBVsMwNgAnTo3z0GNnODdVZN9wH/feeYhjR8bWe1hGlzEL3jB6nBOnxvnoI88yni+zLZtkPF/mo488y4lT4+s9NKPLmMAbRo/z0GNnSLpCXyqBSPg76QoPPXZmvYdmdBkTeMPocc5NFckm3aZt2aTL+aniOo3IWCtM4A2jx9k33Eep5jdtK9V89g5b9e5exwTeMHqce+88RM1XilUP1fB3zVfuvfPQeg/N6DKWRWMYPc6xI2PcT+iLPz9VZO81ZtFYRs7mwQTeMFbIZhS4Y0fGVmWMcUZO0pWmjJz7o3MYGwtz0RjGCtjqKYeWkbO5MIE3jBWw1QXOMnI2FybwhrECtrrAWUbO5sIE3jBWwFYXOMvI2VyYwBvGCtioAnfi1Dh3H3+COx54lLuPP9G1mMCxI2Pc/4FbGRvIMFOqMTaQ4f4P3GoB1g2KqOp6j6HO0aNH9cknn1zvYRjGksRZNKuRcrha44kzW7JJl1LNp+arCe8WQUSeUtWWjZUsTdIwVshqpRyuFo2BX4C+VIJi1eOhx85sqHEaa09XBV5EXgbygA947e4yhmFcPeemimzLJpu2baXAr9GetbDgv0dVL6/BeQxjS7JvuI/xfLluwcPWCvwa7bEgq2FscjZq4NdYf7ot8Ar8uYg8JSL3tNpBRO4RkSdF5MmJiYkuD8cweg/LbDHa0dUsGhG5TlUviMgY8BfAh1T1sXb7WxaNYRjGylgqi6arFryqXoh+jwN/CLy9m+czDMMw5umawItIv4gMxI+BvwM8063zGYZhGM10M4tmJ/CHIhKf5z+r6pe6eD7DMAyjga4JvKqeAd7creMbhmEYS2MrWQ2jDZuxsYdhNGJ58IbRgq3e2MPoDUzgDaMFW72xh9EbmMAbRgu2emMPozcwgTeMFmz1xh5Gb2BBVsNowb13HuKjjzxLseo11VjvpfoujUHkXMpFRMhXPAso9xBmwRtGC3q9vktjENkVeHFijtPjBVzBAso9hFnwhtGGjdbYYzVpDCKfmSjgioDA5UKVQ6M5axjSI5gFbxhbkMYgctUPEAGR8DFYQLlXMIE3jC1IYxA55Tqogmr4GCyg3Ct0LPAicr2IvDd6nI0LiRmGsflobBIykkvhq+IHykguZQ1DeoiOBF5EfhJ4GHgo2rQX+KNuDcowjO7SGEQOFG4c7efwWI5A6bmA8lam0yDrTxPWcv8rAFU9HTXxMAxjk9LLQWQjpFOBr6hqNSr9i4gkCNvxGYaxClhhM6MbdOqD/6qI/EsgKyLfC3we+G/dG5ZhbB2ssJnRLToV+I8AE8C3gXuBPwV+uVuDMoythBU2M7pFpy6aLPB7qvo7ACLiRtssUdYwrpFzU0W2ZZNN2ywP3VgNOrXg/5JQ0GOywJdXfziGsfWwwmZGt+hU4DOqWoj/iB539OkTEVdEviEif3I1AzSMXqcxJ11VF+Whnzg1zt3Hn+COBx7l7uNPmG/e6JhOBX5ORG6P/xCRtwGlDl/7s8BzKx2YYWwVlipsZgFY41ro1Af/c8DnReRC9Pdu4IeXe5GI7AW+H/hXwD+/qhEaxhagXU56YwAWoC+VsEJgRsd0JPCq+jcicgS4GRDglKrWOnjpbwH/Amhb1kBE7gHuAdi/f38nwzGMLYMFYI1rYUmBF5F3q+qjIvL3Fjx1WERQ1S8s8dofAMZV9SkROdZuP1U9DhwHOHr0qC2eMowG9g33MZ4v4wfKRL5C1Q9wRTg40r/eQzM2Acv54N8V/f67LX5+YJnXfjfwARF5Gfgs8G4R+fTVD9Uwth733nmI2VKN81Mlan6AAF6gTBQq5oc3lmVJC15VPyYiDvBFVf3cSg6sqr8E/BJAZMF/WFV/9GoHahid0o1l/+tVSuDYkTF29KfIVzz8QEm5DiO5NIloIZT54Y2lWDaLRlUD4GfWYCyGcc10I+tkvTNZClWfG0dzHNk1yKHRHIPZpPnhjY7oNE3yL0TkwyKyT0S2xz+dnkRVT6jqci4dY5OykfK0u7Hsf71LCdhCKONq6VTgfxz4Z8BXgScbfowtznpbtwtpbEUXc63WbjeOuRKWWwhlGO3oNA/+DYQCfwdhmeCvAf++W4MyNg/dyNNu9HfnUi4iQr7ideT7jrNO4vHA1Vu78Tgm8hUu5yvsGsowkEle0zFbHX85v/6xI2PcT3itz08V2WvlhI0O6VTgfx+YBR6M/r472vZD3RiUsXlozNPOl2tM5CtUPJ/zUyVOnBpfsQjFM4KkK7gCL07MAbBnW6Y+O7gf2h733jsP8dFHnqVY9cgmXUo1/6qs3cZx7BpM89p0mfNTJfZsUxKuc80WdOPxG2c+7d6bNecwroZOBf5mVX1zw99fEZGnuzEgY3PRmKd9YbqMCDgiiMBHH3mWu85P8/iZKx1nnzTOCM5MFHBFQOByocqh0dyys4PVsnYXzkxAuJQvc3G2wu37h6/ZgrYVqsZa0KnAf0NE3qGqTwCIyHcB/6N7wzI2C7HFPD5bBhRUUGDnQIaK5/OJEy+xdzjbkZUKzTOCeFEPEj6Gznzfq2HtLlxBOphNMpBJMFOq8Zl73nFNx251fLAVqsbq02mQ9buA/09EXo4WLj0OvEtEvi0i3+ra6IwNT1woSwmDMwlXuG4oy2A2Sb7s4QXBirJPGjNGUq6DKqiGj2Htske6nblimTHGWtCpwL8POEi4svVd0ePvI1zN+ne7MzRjs3DsyBi37x9m//a+ep42QMULSLvNH7HlrNTGjJGRXApfFT9QRnKpNc0e6XbmimXGGGtBp8XGXun2QIzNxcIMkHce2s7DJ19rCm66jjDU1+yGaLRS22WRNPrQbxztR0QoVDzGBjJruoK025krfUmHs5Phze7QSD+/8v1HzP9urCqiunHqex09elSffNLS69eDlSzFb8wAacxUuev2PTx+5kpdEGPRX7jf/R+4FaDlMeI66L1Mu+u3Fd67sfqIyFOqerTVc50GWY0eplXK3ocffprRXLpl/nm7DJDHz1xZFIC8be+2llbw3cef2LJZJJZBY6wVJvCbnNUogrVQcDxfmS7WKJQ9bhzLLcp+WUkGSLuMlrXMIlmvQmHt6PZ732jvdyOx1a5Np0FWYwOyWmUCFi7Fv1yo4Aj4qi2zX1YjA6TbWSRxfZy3/eqfc++nn+LlycKGKKUA3X3vG610xEZiK14bE/hNzGoVwVooOHHOeaohA6bRwlyNDJBuZpE0fpHLtYBAlclCjULFW/NCYa3o5ntf78JoG5mteG3MRbOJWelUv930dOHyftcRPF8ZHUjXX9toYXaSYbLcVLibWSqNX+R4sZQCE/kKA5nlS+3GY3/h0iw1X0klHA6PDaxofEu9/26+d1tA1Z6teG1M4DcQK/UPrqSw1nK1TxoF58D2PibnqriOoKot67kstVq00zor3aqv0vhFTrkOnq+IMz8zWcodEo+96vnMlr1w/6rPy5OFZVfhLjzGUu+/W+99NYut9Rpb8dqYwG8QVlp8ClZWWGu5zI2FghPfbK7GwlyNLJFrCYY1fpFHcmkuzJQgCP2Rpy/lqQUBSUdaFkOLxz5Z8HAQHEcIVJkteewaSrR8DwvHOjVXWfT+J/Jl7vvsNxjMJrsa3FutYmu9yFa8NibwG4RORLGV6N3/gVsXCTHA3cefaNpvpdPTa7Ew43PF1SWrfiioM6VaR6+/mptdI41f5IFMgh1eism5KrVASSWEvduy1AJtecx47PU6OIBEtXBaXa9WY315co6927L1fWZLNSbnqgSq7N/et+L3sxKstHB7tuK1MYHfICwnwG1F7wO3NuWet9tvIJ2gVPPXZHq6b7iPlycLTBZqiIDrCLVAyZe9jkoIX+sMYOEX+eBIjuG+KlU/aHr/rY4ZW/91147M18Jpdb1ajTXpOFzKVxjMpoAwKwkgk3Drwb1u5r1baeH2bLVr0zWBF5EM8BiQjs7zsKp+rFvn2+ws5x/sVPTa7aeq9cyNbk5PT5wa5/xUkddnKvVtCScsHzzcl1xS1B788gt88utnmS17OAIj/Sl2DoWWcDvrealAZuN57njg0Y5mMLH1P5BJhFZ3oKAw2J+k5ivvPLS9aXb0wqVZdg9lm46xczDN+ely/VpXvACBpqB1rwf3jI1BN9MkK8C7ozrybwHeJyLXXme1R1kuda7TtnHt9pur+tz/gVsZG8gwU6oxNpBZ9aXxJ06N8wsPP83rM+Wm7V6gDGeTjOTSbUXtwS+/wMcffZFSzUeAQGG8UOXSTAlYPNtYaU5zp7nncXXMgyM5hjIJskmXob4kB3bkuOv2PTx88rWmcxYqft1Cj0m4DjeN5erXui/lMjKQqneDanfutWIj9dA1ukvXLHgNi9wUoj+T0c/GKXyzwVjOP9hpBsBS+3U6PV0uwNnu+YceO0O+7OFGgUkNjV9EYK7qLylqn/z6WRyBhOMQBD5+9EkZL1RxHCGVcJtmG0vNaOLnG8fXKsA2W6qRdIQ7Hni06X00XqfG9/rs12foT7sMZTP1c27vT3JlrkZ/OtE0M2osHBbfjDZCcO9a4xvG5qKrxcZExAWeAm4EPqGqv9hin3uAewD279//tldescKVrei0QNW1FrJqfL3nB1yarVALAg6P5vjI+28B2hcJ++U/foaLMyUSrkMQQM0P6nd014E92/pajvehx87w+JlJBKIVtM1jSiccfvrYDdz33pvq247+n39BseKFgVPXYXQgTS6d4OJsub6AZeH4v+9Nu+vF0HLpBBOFCkPZZMuCaXE/2Mm5KoPRPqcu5hFgz3DokonbEwrCjWM5ChWP/jY9ZK8lK+lqaXUjfuixM4sMgGI1rNS5Go1MjLVnqWJja1JNUkS2AX8IfEhVn2m3n1WTXJpORWKl+7USgMYWfKqK4whjAxn6U27LYOXYQGjVfuPVKZSwbZ8fKLUo93wgk+DBH3lr25vR2ctzBA0fRRFAw98HR/qbBOjEqXHu/fRTBKpRrn4YCN2RSzJX8RkdSLcdf3yDufv4E4uE7nKhzJW5GnuHs2STLi9OFPB8Ze9wloFMkjMThXp2TRCNrfHYsQtnI1SJbHejn6vU2D2URaIMIQjfw0ypxtd+8d1rOkZjdVj3apKqOi0iJwgbh7QVeGNpOnWxdLJfu6l6seqxazDD2ctz9f6qCviBknSFM5fnODyWazpWHAv41Q++kV94+GmmijVUQrV2HWFbX5LfuOvNbXPO+1IJRvpTjBeq809G7p2xXHpRrOGhx84w3Jdkcq6KBqHQBihX5mrk0i7ZpNt2/HGQt1XW0kyxhh9oXfT9QHFkfgVsnFNf9gJSrjS1J0y4wie/fpbRgfSGqBLZzoVV83XNsqmM9adrQVYRGY0sd0QkC7wXONWt8xkro11djqoXUKr5VP2A2MiL0wTj4G27YOWxI2P8+l1v5vBYDhFBRLhxtJ/fuCvs174wsNcYEN45lGUsl6ofUwR2DqQZG8wsEqBzU0VGcmmuG8qScAVfQzfNQNrlpp2DS44/vlG0CrpW/IB0Yv4rEdfiiVfADmaT7OhPhZY7ze0J40B2J4HwtaBdsD2VcKyT1Baimxb8buD3Iz+8A3xOVf+ki+czVsC5qSKuUHc7pFyHkVyKlCvUfA3dEIGGGS0oI7lQaA/u6KNYC9oGDFvNHtrNFnIpt8ma3DmUxXGEqeK8m6SVAMWB5MFsst4eMHYTxcFUAcq1oP6aQOHUxTx9KTd08bQIuiYch4HM/FdiJJfmtekSiYaSDamEy5GdA4vcVKWaT/+C9xNvXw/ruF2wPa6ps5UW+2xlrKPTFuV9v/lVXpyYwxWpL+bxVblxtJ+PvP8WHvjSKV4YL5B0hZ0DaSpewFSxxkAmQSbhcGWuRsUP6E+5/MQdB5sCoAtp5e8uVj0qNZ+JQhU/UNIJh6QrlGqh+0NESCUcRvpTi4KWsHQ3qAe//AL/9isvUlsQrXUFxgbTJF233lWqUehadaCaKdUYzaUpVLymlcLtOlpdqw9+teqVb/WuUVup7vu6B1k7xQR+7Xj/bz3G6fECrtMg8IFyeCzHF3/uTmD+S3L60iz5is/2/iQp1+G16TDPfc+2DAnXWVY44kVGjYG92VKV89MlxgbSzBRrVLwAX2G4L8GebX11cRWou0AWtvxrVaLhocfOcPLVKURgIB0tVooComnX4fDOgSWzRq41kH0t2TKrLcrrkbmzEdhqNzcTeGMRdzzwKK7A5UK1yUUTKIuyKRot8DMTBbxodWfCFQ6N5loKZqMFNVuq0Z92Gcll6s+fHs+DwuGdA0Czq+jQaBjEff5iWK434Uo9FdKNMlYWinPjl/qVybmm4GrKdUDCx0d2DS7KGtko1l67mY6lMK6MrXYd1z2Lxrg6uik8sY82FlNoTndspDHjpOoHuI7UH8N8ILFu8Y/nyZc9hvvCzBPPDxjPhxkyO/rTdYtq77b5c1X9AEfmjzlbqlH15zNxvCjtcfdQ69WwjUHjdMLFCxRR8AkXXMF80LTRL76RFv5sxXrl3WArXceqFyz5vAn8BmW1hGe5Jh+XC+XQReIHJByHD775ukXHaAzYpVynbsE3CmZ/yq2Pt1jx6ouMxvMVskmXgbTLXMWn5pWo+mH64aXZCiLCQCZ0/cQWPIQFumKHjiD1VMhLsxXeun940Rgbg8YVL6inOIqEWTYo7BpM14O2cU2Zk69OIcCuocyKC4Gt9g14K9Yr7wa9eh1rfkDFC6jUfCpeQNULu5UthQn8BmW1aqo33iRenixw76efIhelE75t/xB/+swl/EDJJFySrvCJEy/xqSdeaepg1JhxMpJL1X3wu3LzgplyHWq+z2TBo9iQvSKETbyrfkAmIfSlMwxFq0zPT5V4eTIUZtcJF0YN9CdQ1XqBLseJsnmiRUWetk7pG0gnOD1eCG8EqqF7RsFFuTFK2yxUwhlKYzA10DBT6MJ0meu20VHHp1bX9mpuwAtvEPG4NkJJg81ML9R992Ix9wIqnk81MloApotVTo8XOH2pwAvj+SWPYwK/hqzE4mucZsZ11Suez/mpUkcld6H5JpEv15gs1FCUci1gPF/m5KtTbO9PMpLLMFuqhY0xgGLF4+zlAj/5n54k0FAwRSGTdOlLuxwey6GqzFV9Uq5D0lGev5Qn0LByZCOBzot0saaMReOZLdXCOjXRPklHyCUctmWT9QJd/WmXlOtyuVCpu4YObe9vCmbG1zJfroU3gAUzVlX4yPtvabpedx9/on5d4rLANCxo6sTau9YbcKsbxMMnX6uXSrjWwOhGiSusB5ut7rsfKBXPp1IL6pa5F4Qf5MlChRcuFTg9ng8F/VKBiQXF7ZbCBH6NWKnFF08zG5fcO1FKY6eWYuNNYiJfCY+B1HO4PT/g4kyFy1GqoiNCwhXKXkAl2hajQLnm05d2ef8bd/H4mStcnssznq8w3Dfv72x8TUwQuUhUtb745nKhQtJxECfO3gmzW4b703zp599Rv14JVzg40l+3wj7y/ltaXsvXZyuLatgABMAv/+G32LcjVxe7Zy/M4PkBtUBxgFoQ+ukrhB2fBjIJfuX739DxtY3JJl1OX5pd1Gyl1b9TuxvE42euXHMgcCPFFdaLjVr3PQi0LuIVL3S11PwAVWUiX+H0eIEXLuWj3wWuzFVbHscROLCjn8M7cxxf4nwm8GvAiVPj3PfZbzBX9cgkXEYH0gxkkktafPE0c3y2TKiOzcviO7EUG32RVT90ecSBy9OX8nVBdJ1wcVPs2ggUaOHbU8Lg5ydOvMTe4Sz5Uo2Kr1ycrTTts5CEI6RTLjNlj1MX8ySc8CYCID4knPksmsYZSjsrrNECh1AcXQG/xbkBzs9UuByVIRifLVP1FdeBpOvg+UG9Bo5E/+skr6yVn/dyoUK+4i8qYdxKWLsZCFwN914vsN6zmNjVGLtZKrV5Mb84Ww4t80v5yEIvtO14lnCEAyP93LQzx+GxAW7amePQSD+ZaPW5Cfw6cuLUOB9++GnyZQ8FijWfc1eK7NveRy6dWLJl3v3AvZ9+CgWSrjCSyzCYTaKqHQlBoy/SASqRoicdqDT4MjQILYJAW1vg9f2IbxDKy5PFJfeNSbvCYDbBeL7KYMZltuRR9ppnBrUA1AuzaITmGUqrL2QrcWx1Q2p6OgBXpP6+/QAS0XuG8LzX7+hre+PtxF8+VayxvT/ZkbB2MxC4lbJI2rHWsxjVMM4UBkFDQa/5ih8EXJgu1cX89Hgo5vmooftCkq5ww2iOww1ifmBHP5mozET8k044pFynaW1JK0zgu8yvffE5povzd2bV0NJ8fbrE3u19S36hjx0Z4/b9w1ctBI1W8KuTRYSwbK/rOtQiH1/CCfPZG+KidbFvhxBONTshUJir+AykXaqetnSjQFTYy3XYNZTBdZpnKAvFdWGJA6Bp/K2I0y+bT6/hIijCWUzckGOhGHbqL58p1djRn246Qzth7WYgsFezSFZCt2cxjS6Wut/cDzg3VYx85aGYvzheoFhtPbdMJxxuGO3n8M4BbtoZivn12/vIJBOLxDzpXl3ZMBP4LnN2sogj4LhOvXSuElrTS32hY1F79sIM+YqHKGRTLgOZxKLmF0sRW8Fv+9U/p1CJo/HhONwojTDOhX/2tRlY2iAAwhuC54fZJ0vJfMqBnUMZilUvbMO3zMGv25ZhIDM/Qzlxapxf++JznJ4okHQcBjMJnnrlSo8v/iAAACAASURBVN3NlHKF3UMZKl6wrFul1fMB4ZcsiIqVxSwUw0795a0W2LQT1m4GAnshi+RaWc1ZTKv0xJof8PLkHKcj98oLl/K8NF6oux4X0pdyuXEsx+GxXCToOfYN95FJuqQTDunEvIXuOh18CTvEBH6NCP/RHLwgQCOrsd3S6dhirPk+papft6iLVR8vUH762P4V58IXKj6BhjVf4rozQaD4vvLshZm6xd6fdDg0OsALF2ep+tpSGDtd/FwNIJdOcGEmDBbLEjIsQt2CbsyrH8+XcSVs2t1UTpjQXfTKlVJng1l4PuKZSEBA6EaKC4otFMNOxWKlwtqtQOBmyyLpBlc7i2lMT4yt9HLN5+zluaZsljOX59ouMupPu3X3Svx77/ZQzFOuQzoZulfSieVdLEsRG4HJ0QNvarePCXyXOTTSH+Znq+I4kBQHP1B2D2V46LEz/PIfP7MoABRbjJMFD8cREuIQSLhkf9dQhj/99uv1rkOdBI8eeuwM2/uTUZokiAOBp9Q/npHuOgJz1YBzV+bYNZTh1TbiGWg4loXFvFpxfqpIXA5jqb0F6gI7U6pRrPpUPJ8gaN3lqVNcR8Jc/AUupVRUZth1HD507IYlUxM7FYuNJKwrvXmsd0BytenkZhunJ1a9ed/5XKXGmctz89kslwqcvTy36PMTM5hJcNPOgSaf+d7hLOmkO2+Vu6Flvpo0ug3RoLVDH6tF03XiRtT5socXhKtF08nQp9bYLm62VGNHf4pC1WciX2HXYJoLM+VIoMJytb4q1w1lOD9d4sCO/rYFuBZ+SeNiX4WKx8WZcttppBtlkAQ6v/CoHX1Jh6G+JBdnKku7aaLKkJU254wZzCQYyibpj9rkTRVruE5opV/LRzTlCo5I03tOJxx2D2W4OFNGgdv3Dy8qFtZ4DSEM/FY9n3zZo+KFOfkL2whuVnq1OFdjsbU927L8+Hcf5LsO7aj7zmfLNV4aLzSlJr68oLNYI8N9ybp75aaxUNT3bMs2B0Bdh8RV+stXQqM78M/u/wfF6sQr/a32M4FfAxZW9Zuaq1Br6ByUL9c4P1Ui4Qo3jubqreISjtTb3wWB1ptbNBbpgrCGTNIRLs9VKVTClaVBELpXElE64FA2SSbpcmG6TC1yEy3F9duzi9wfcdVJCG8Abz+wnWdemyZfaZegGOIAy8RA+efvPcx9772p/sG9OFPGi1I3a0vcaJaLA8S4DbOAwbRDsRb+0VgRc6lyv986P80nTryEFwSk3fDmlnTdeqD19HieqheQdIWbdg5uKgu414pztUpPnC5WeXGi0JTNcu5Ksa2Yj+bSHN6Z46adOW4cy3HzrkGu25atW+NxFouziv7yldBYoXUpgTcXzRqwcLoc/+PETOQroRsiUESEnQMZXpsuhXnpCkHU/m4gk2SiUG0q0gVEDaFniSOkjZa3F4BqWBcm4YRB1eXEXYALM+VF2xe+7vR4vqksQTuW2kMILaPHz1zhPub93XF7PEeWl/BkB+4iX4mms0q+GnZu2jmQaWoYslTLPYC9w9kmEZzIl/nEiZfY3p9kplgDgVINzl4utE3J24iukM2cVtkqPfHKXLXuXokt8/NT7WM1uwYzdTE/PDbAkV0D7BzK1MU81WFK4lrSym3YChP4dWDhP068CCnO5AhFJ1xANJhNUI16gB4cyTHcV62n/MWUan4oYM5iXzNAwnFQQks4IcunQbpOuOI10eZ4IcpMsdZxuiQstrYdCd/zbLnGyVenOHFqvKlbE4SLh7wgFOhWMwFHOk/ZTLkOB0f6ee5inhtHc01f2Ljl3v42LfcUFolg7HabLYWxknimlS977BpKtMyl34grTDdTWuXC9MTx2fK8mEcB0NdbGCcxe7Zl61b5TTsHeMPuQUYG0k2W+Vq4WK6VxhjDUpjArwMLA0CuhOVwR3LzOdTVqD9oJuk2Ff6KRWJh8Cghoc+8ldZV/QDHCWuq7x3O4vnKq1eKLe1ih9DKD1RbVqqLc8aTrhPVqAmX+9dXwC5DLPJxzZrYBRWXYIjdJMWqx0AmUQ/mBr5fX40a1q5PU/F8ZssexaqPELqjlppQlGr+kq31lmu5t1AEK17ormksoSxRyeNWFvBGXWG6UdMqG9MTq364YOj5i81iPp5vXZdFgH3b+zg8FlrmN+0c4JbrBtnenyLtunUxXy8Xy7XSGNBHnLY63jUfvIjsAz4F7CI0vI6r6seXes1m8sFf61S70S+fSyeYKFTqQdfJuQrj+SqjuRQjuXQ9sySbdBnPV/CDgJQbF/4Kxf/DDz/N5ULruhUx+4aziAhJV5gtVbkcZdWkEw7DfcmwhjowE9Vijz/6scjGOeNu1FBb0aZ8+OG+BBOF1sutIVzVWotuQik3EkSEAGVbNkm+4qEKo7kUV+ZqlD0/LOGbdNg73Nd0jWIR6ks6TJdqTBSq9bjDUqQTDt/3xp18/cXJpsD3QCbBP3zH9cv64OP2ggOZRH3l6mzJw4tq+cSxkl1DmUU+7FadrRY2H1kv1rv7U1P1xJrP+akipxaI+eQSdVmu3xEv5Q/95bfsHmBbX2pFqz43K+vV8MMD/jdVPSkiA8BTIvIXqvqdLp5zTWjMU58p1nh9psTJV6c6zqpol6kRf8HilZ/5ssfkXBVHhGq0mCfOCqn6AX3ML3hqV5SoEQ0CfvV/uY2HHjvDTLFKKuHUU8RmSjV+6l03cNvebfzUHzxVj0gmokqQMF/ewEdJJ11UQ8FMJYTporekuMN8qYSUI/U2erFFOx215/OD0DUVREtM/UDJV3xemiiQcIRixW/qBfu5p85Huf2dGSqqyqPPT9SPLyL1+jO37d3GbXu3tWwF+PDJ1+q+9rLn4xWV2/cNcfLcDLWono3rhLORgUyypQW8kV0ha1mcqzE9sVzzeWWyyHOvzzalJk63qcviOsLBWMx35jiya5Bbdg8ymE1Gn8WrX/XZi3RN4FX1deD16HFeRJ4D9gCbXuAfeuxMVPu8FolUmNv+iRMvcdvebUt+Udr5Ye+6fQ8QCk2xGi5KSjpOvX5KLF+uE354JVAKlXB6PzVXWdY9knBgYq5WH9svPPx0vSlGoFCqBTz46Gnue/dhduTS9XZ+xaq/yCoOCIUp6QiplMNs2Q8DmNBRSmQ1mF+JWvUDSnM+jgiOE84KElGKZrwgDMJ0yaqvpFzhll0DlGo+D598jVzK5aWJAlFy0ZJINL6ZkkfKFW7eNVh/LnaVfOaedyz694uLmw1lM/W2g5cLZU6em2G4L0m+7FGq+QQapo8eHMm1tIDX0xWyWjPOlb4+CKIgaC2gVPM4e3mO71yYbUpNbFeXBUIDIuGEVU5/7J0H+MGj+8IGMV1Y9dmLrIkPXkQOAG8F/qrFc/cA9wDs379/LYZzzZybKjJTrNVL+EKYhucFuqw/NfbDer5ydmauHmB98NHTjEYNqOPAZiV6rpVwSWTdnp8qMp6vLJsu6AXgB0H9i3plrrpo8ZAXwG/+5Wn2DKZJpBKM5NK8cqV9JkUtUKaKHq4DvgpJR1Y0Da42DCBQrbuI/AZXy8L3VPWVs5fn6E+5TDe09euU+N9rYdbNUlkjrbJMZqJ4QDrhksdrWlB1bqpYz7xp/Cys10Koaw3udvr6xvTEYtXjxfFCaJk3rACd66Auy8lXpqj5AUPZJK7jIAKlqsdfvzzFz/+dm1fpqmwNui7wIpID/ivwc6o6u/B5VT0OYcXLo0ePbpyk/CXYN9zH6zOlpqlg7K5oJRKN1s9EvsJgxmW6FNZmabTQX59ZHDBqdUH8QOt1bWZKYSZLJ6s9XUf46CPPcnGm1HZfVbhUqDKUVWbbTJMXjyccqQah66bT3PSFxK9rn7kTUvGCtkKx3LErXvg6JSx9HGfrLOUqaeVaiVscXpgphTV2NJqZ+IorYarkvZ9+ioFMoilIvh51yq81uNvq9XOVGr/91Ze4/cAwcxWP05fyPHshH7lY8rw4UaDcJuKdTUZ1WXbmuHnnAG/aM8SNO3P1Bizv/tcnGMmlm4yFvlT7yqtGe7oq8CKSJBT3P1DVL3TzXGvJvXce4uSrU/hB+GVWDfuFDmSSi0RiofVzOV/hcqGG6wiOO59JErOcOJZqzcI2u8T0dhGqVD1/SavXiVrejebSTC4TtF1I7Lq5Wjq5KbQqO7CSYze+9JUrRdKuMNyfWrKA2zsPbV+0yCnhOARBgCNhJkYtajEIhDGEKCZQrHjrng65VJ77g19+gU9+/SxzVb8e11gYRzo3VcRBeWm6RMULSLhCJulycabMDz/0OC9NLFeXJcwvf8PuQd64Z5BDozmyKbftqs+NHKvYbHRN4CW8/f4u8Jyq/ptunWc9OHZkjJ8+dkP0pY+zKsIslOlilTseeLTup1xo/ewayvDyZBEvUBJui2X4V2v+doCnLOnvhPBmlXQkbAvWxbFcLas9nIqvTBVrbQu4nTg1zqeeeCUsShZAMQioFar8wJt28ci3LhJGDaiXHU44Eq1bCLtVxSuWF1rMa7ngqZ1gahDw8UdfxJEwRlOq+Xz80RfxA+Un33UD+XKN71yYpVz16xksCvieUvHCz9Fzr8/3BB3MJDi8c4Cbo5TEN143yMGR/nqRrU5TEjdq2uZmpJtpkncAXwO+zfz6lH+pqn/a7jWbKU0SmlPL4hoqgwvS+IpVj12DGfJlr95b1IsCiEk3zE2PA3TLsRp6m3CkIwt4A2p7VxDg0Gh/22X57/vNr/LixFyUGkq9EueNo/2ICGcvz+GrhrM5R+otEdMJp57vf2g015QOuda1X9qd73KhQtXzw4VwUTevON3z+h19yzZ1ySYdfvDoPm7ZHYr59Tv6SEdifq0pieudtrmZWCpN0mrRrBLt6nnEZQhmSl69yfTVXvEw6+Ta/r0SThgQXqZ8zKYnXpDVGLBdal8BrhvOLrKmb/7lL6JR1cmYqucTANv7khQqPtv7k6Rch9emwxWUroQWjSD1GveNtV3Wo/bLiVPj/PuvvsS5qSJjAxm++8YdfOIrL9U/kEtdo3jFsRdoVDtfGMmlAOHrH3lPV8ZrdM565cFvGU6cGufkq1P1xhH9qXDZe1zuNp6+XKtVfDW+54X4QdhRqtctdG346XTfSm3eX37X+WkeP3OFSuRbVw1IRD1c43vs7qEslwsVrszVomBqaKnHvVm39yfJpUP3TKOLoW3D7vF8Rw27F9LK3fOum0epeAET+QrfOj/Nsxdm6UsnEBG+eW6ab5ybnn/zC3AFfvyOg9x63SC37dnGR77wbS4XyvSn58cc35CMjY0J/DUST39jKzDO8IjLPzeGnq5VUFdB3+tj6KRA12anscFIJze0iUKNYjVsivybXz6N48y/Lqxo6RO3k81EzRpGBzL0pxOLrO+lXAytfOKTcxXyZa+jht2NnDg1zq/88TMkXCGbdHhposCHPvMNDuzo4/Jcdcm6LAPpBPmKF3Yci9xPCtz37hv52e+dT0f8Z8duMJ/4JsVcNCvgxKlxHvjSKc5cngPg4I4+JFpV6gcaluL1g2t2xawVm2GMq8HVvs84JtjYmFsJLdx92/vqHahUlYuzZQ6PDXRkfbfyiZ+fKjHcl2S0wSpu57apegGvTRd5+vwMD3zxFNPFKjVf287whLAS5s27Brhl9yBv2jPEm/YMMTKQ5re/8mLLLJpWTcaXaopirB/mg18F4sYdU8Va/Yvv+WFXJFcgk3TJpROM55dugGFsPmLr1nGE/pRLf9qtr2iFcFXrlbkae4ezdcGeKdUYzaXJV7xF5SjOTRXJpVxEhELFY+9wHy9cmmX3UHZRnZqpuQqf+ol38PS5aZ55bYbHX5rkxfHCkjXyF5JJCP/u77+Nd79hZ0f792oDkF7FfPCrwEOPnSFfDlcsOiJ4flB3vwQa+senijXSibC64Gq4U4yNg+PA2w9sb5nCd2WuxnBfsu5y8XxlulijUPa4cSzHeL7Mhx9+GiEsBb0tm4xEM+BXP/hGjh0Z4+7jT3BxpojruJRqHuVawFzVw/OV9/zrr7YdlwDppEMm4bat31L2lI/9t2c7FviNWvXSWDkm8EuwcAVqzQ/qvRX9hplPuEJyfqGLafvakHLCxt6N3Zq6QaCQdh3eeWg7Dz12hmLVY7ZYoxZl6ITlFeajLZcLUQMXDRu49KUSvDZdAoX+dIKzl+co13xcgV/6wrf43lt38crkXMsmKzFJV0g6Do4TNjLPJhNkUg5Tc1Umo05eS7FUw4uFbOYGIEYzJvBtWLQCtRAW9KpFZQUWWuiNft5cysELaNv71FgdqkGU9uk4+F2+1rfvG6qXEha0XjAtGdWgDytpltg1lF3UwCUIAmpeOKt7dbJYz9rxFV6frfCpx19ZdD4Bbtk9wA8d3cdte7dxy+5Bvvc3v7qo3PBILs3lQoVEVPCu3X1uJfc/W0naO1hdzTY0TlPjNnpxvZdW4g5hATCAYi1YkY/UWJqllsx4ARwey9HNCrHbsgmefT1f/zxMztXq4/J1Xsgn56r4vg8oNV8RgRfHC3znYh4/ylAJWCy2jsBAJsF1QxluGO3n1t0DHBrtZyib4se++yC3Xz9MNuWyb7hvUSmIUs3HdRxuHM1x/Y72ApxNdn6B7r3zUH2RnqouSvM0Ng8m8G04N1Uk29C+bTCbJBEpzULBqedcR9/cQFlyBaCxMpa7kmcvz+F30YAvVjxmyx6vT5fIl2v1G3w8kws0rIIZKHznYoGaHz5XqIQdpJbKY0i5YY/c67f3sSMX9oN1HKelS6Sd8B4a6adU8xnIJBe5ViD8vP7Uu27o+P0eOzLG/R+4lbGBDDOlGmMDGQuwblLMRdOGVtNUH+hPuRwazQHwzIWZZRtYG91npZUlV0o1unlUvIBXW5RPbpWemE40L4rKpRxK1YDGkboyf/N6/uIsSFi6YnQgjevIIpdIu3LDQD3wu3c4iyNwpTg/y9g7nOW2vdtW9J7Xo+qlsfqYwLehVbZE3NpttlTjcqFi4r7FWC6Qu3swTS6dZDxfolQLGBsM6/vXAqVYCxYl5PsKGt0cvACSbtiHNM6J/5Xvf8Oic7QT3kbh3zmYwXGkqb3hRmjubaw9JvANLFzccdfte5oWd9y2Z5A/fvp1S4Hc5KQTDodG+nnuYn75nTskl3YZL1Sjgl2hjk8Waviq9QYvC4ldfUlX8DXsDFb1w3K8O/pTbatbtqpC2Sj8dx9/ImzpaGmOWx4T+IhWXWsePvla3fcYL3Qycd/8VLyAAyP9XJgpM9NhUxMAl8VBUhFwJRRwR2CqWKv3zxUJ91SdL5rQaMSLzAfsMwmn7vqLK08upNPOSpbmaMRsOYGPLaBvvzZFqaaoKrl0gsG0SzqVoC8174KpeAH3ffYbPPgjb60vdDJ6gy8+c3HZfRaWOGjl6VcFL2o3mHTC+vL1IHyUOeNHrQgXFj+Lxd0PlN1D6fr2dimJnS5AsjRHI2ZLCXxsAc2WqhQq82kXs+UwS2IwXWO4P8Nr06V6TvFs2Yssd23btcbYfGSTDuVasGSGTqeTtfmCZPNNP2oNOekL89Md5sU+4YSvTriCqi5ZyKtTy3yzNsxYyyYoW4UtlSYZW0Dt2tzNVsIsCW/BF3KiUGVyrobJe+9QatMv9GpYeCNwXalXE42fv357mN0iAAJ9KZfrt/dxZNcA27LJjlIS2+XBt8y22WRpjrHxtbCa5olT4+s9tE3NlrLgYwtoKT+6udi3Dqv5by0SWku+htUes0mXRKD4qrgiXC5U6+V40+68v71Y9Ti8c7CjRh8rscw3W5qj1b/pDlvKgm9lARnGapBywsyXhCOkEg4HR/rDDkhB6J/3gjCbBsISFrOl6opXiG5Gy7xTFi4sBAsMrwbdbLr9e8APAOOq+sZunWc5Gv16A+kEM6VavbemYawGEv1/qlRjuC/BbDksF9yXcinXfETCCqSIkNCAQJWLsxVu3z+8Yj/zZrPMO8UCw92hmxb8fwTe18XjL8tCv15cBCqbcJd9rWG0QqBuicckHMGJGm7nKz637x/ma7/4bh78kbfia5ippRpa8yKwbzjL2ECaz9zzjp4U66vB6t90h65Z8Kr6mIgc6NbxOyH263m+cnZmru6esVx242qoL0xyHPzAJxn1aHUcQaP/qjVluljljgceZSCdwIGoNaKSSjhcN5gl4cqy/UzjmecLl2ap+eFrD48N9GxmSbsyDL34XteSdQ+yisg9wD0A+/fvX9Vjn5sq4gpcmCmjgZqwG1dN0hWuG2oW5/F8GT9QJvIVqn4AGqY7Vv0AV+D0eAHV8MbguoIXBFyYLhIQ1px58Msv8PiZK4vSAuOZZ9Xz6xlfparPy5OFni450Kvup/Vk3QVeVY8DxyFs2Xe1x2mVQzuQTvD8pbwJu3FNCHDdUBYReH26xMuTRQ6P9jNbqjGYTXIwquZ4fqrE9v6ws9OZiQKuI6BSLyNd8xUV2DecYapY5eOPvsjYQIod/emmVanxzHOy4OEQun8CVWZLHruGEpZZYnTMugv8tRI3wn5hvICD4jgOF2fK/PXZya52+TF6F2G+tr8Q9tsVgQvTZRQlk3DqC5mqXlDvluT5AVcKVS4Xqnh+mDXjiOCrknKdsDQwMJhNcWaigCMwW/IYyWWa0gLjdN5wJhCFcAWqfmCZJcaK2FQC36rT+8MnX2N8tgyq1BQIbDmScfW4AqMD6XqnpIlClcFsgvHZUNwh7KLUl0owV/EYz1fYO5zF8wNevVKi4ispN1zRVPWVhEO9T29jl6eqH+BEoh0Ti3ecUZJyHbyocYhGjUUss8RYCV3LohGRzwCPAzeLyHkR+SfXcrxWK90+ceIlqp4fZiqszrCNLcy2bIKbdg5wcCTHTKnGwZEcP/vuGzmwIxcJt8N1Q1kGo3IB+bKHF4RVGy8XqrgNrphEZHl7gTKSS+GKEGh4c4BQrIOGblAwnxYYZ5QMZBIEKF4QEATKYDZhmSXGiuhmFs3dq3m8VivdvCAgX/ZIuQ5zvi1gMq6OTMJBNWwcMpEv88Wfu7Pp+fsIS/AuzNOueAHpBos84TpIEOAFgISWux+EvVgPjvQzUajUa84MZBL12cHCGjSNGSWeP0s1yqI5sCNnmSXGitg0LppWhZbSrkPFC9g7nKUYNTM2jJUiEgZCg0CptgnctCoT4DrCUF/4mYzdKY7j0JcQDo3mKFY9xgYy9TIEsYvx/FSRgyM5/v7btzf1G2gUb8soMVaDTSPwrVa6DfUluTJXw3WETAJKVs13yxNXdpSosFc26VL2giV75CoarmxWSCVaey1b5Wl/8M3X8fDJ1yhWPUZyKV6bLoPCrsF0y4U6jaJtlRONtWDTCHwrCyrpuvz0sf187m9eNXE3gFDY9w9neW26jBAGTC/OlJcU+EpUWTLpCqO5dNv9WlnVt+3dVhf9G0f7EREKldBybyfanTbuMIxrZdMIfCsL6p2HtvM7XztDvmL+dyMkGQUv+1Iu/emwJEWrptgxjaUHAoWJQoUTp8Y7trSvxpVilRONtWLDC3yr1MjHz8ALl2Z5/Mzkeg/P2GAc3NHHl37+XXUr+eJMGUfCejGNQu8IDGeTzNV8/CDMkBnJpUm4Uhfablna1lLPWCs2tMAv/IK9PFng8TOTi1qpGUaMROmJ8Yzv3k8/FdZgTzjsyaUZzCbrPU8V2L0tW38NwGypyslXp7jjgUeZLdXoT7sMZcPSBKtlaVvlRGOt2ND14BunsoWKx+szFcDE3WjPRKFSf3zsyBi37x9m//Y+Do3m6vnrsZgu7A+QL9dC373AtmySYtXncr5KvjzfAHs1LG2rnGisFRta4OMmAOOzZV6ZtOmrsTz5stfU5m0pMV343MWZMgA7BzKICOmEAwIT+fmbxmpY2r3cuMPYWGxoF82+4T5eniwwnq+Y1W4AUY0YbT2Lcx1huC/Z5EJZrgxt43MK7NmWqVv6I7k0F2ZKlD1/2YbYK8Xy3I21YEMKfGMt7CtzNRN3g4G0S77it+zEJYT567sGMwxkEotcKEuJaeNz8WrVmMFskornU6yGHZqsRrmx2dhwAv/gl1/g337lxahJgrFViMOcrf7V+5Iub9yzjVMXZ5gpeQQaZsEIkIwaYcQUq95Vu1BarbVIJVx+7e/dZqJubEo2lMDnyx7/5sun13sYxhqTdIUPfc+NAPzWX55uqt/vCPzTdx3itr3b+Ogjz7K9P10X39koE6ZRkK/FhWJdhYxeY0MJ/PmpIqPrPQij6zgSLkjyAuXtB7YvEtFPfv0sc1Wf/pTLT9xxkPveexPAIvH9le9/A7C6gmy+caOXEG3l1Fwn0rsP6+5/9FvrPQyjy1y/va/e+i4uxGUYxtUhIk+p6tFWz23oNEljc9GfchnNpep10VsxlkuRcMXyvg1jDdhQLhpjcyFAf8qh5CmjuRQjuTSlmk82leCu2/fw+afOc36qhAJJB3YNZQGWLMRlGMbqYQJvNOE6kE64qIa10Rf6waG5rvneen2gxXXNG19jGMbaYwLfI7gAwqJG4wkHhjIJSp5Sjpblx2GXhCPcMNrP971pd9vGE61oFYi8bxXfi2EYq0NXBV5E3gd8nFB/Pqmqv7bU/umEw/Xbs7xypdTNYW1K0gkHR6BcC+q54oOZRN26XmhVr9QFYgJtGL1H1wReRFzgE8D3AueBvxGRR1T1O+1e44iQcB1cB7b3JZko1Nrt2rMIMJBJ8J4jo1ycrV6TVW0Yxtammxb824EXVfUMgIh8Fvgg0Fbg/UAZG8iQdIRaoFRqAbObpJlH0hU8XxetxEy5ws98z43mjzYMY83pWh68iNwFvE9VfyL6+x8C36WqP7Ngv3uAe6I/bwaed9K5QXdwZD+qiuMmxHHWJ1agGmjg1wBR36sEc9MXg0phtgtnGgEud+G4vYZdp+Wxa9QZvXSdrlfVlmtEuymcrbKhF91NVPU4cLyL49jwiMiT7RYqGPPYdVoeu0adsVWuUzcXOp0H9jX8UktQkQAABi1JREFUvRe40MXzGYZhGA10U+D/BjgsIgdFJAX8CPBIF89nGIZhNNA1F42qeiLyM8CfEaZJ/p6qPtut821ytrSLagXYdVoeu0adsSWu04YqNmYYhmGsHlZszDAMo0cxgTcMw+hRTOC7gIjsE5GviMhzIvKsiPxstH27iPyFiJyOfg83vOaXRORFEXleRP6nhu1vE5FvR889KCJLFOPdfIiIKyLfEJE/if62a7QAEdkmIg+LyKnoM/VOu07NiMjPR9+1Z0TkMyKSsWsEqKr9rPIPsBu4PXo8ALwAvAH4f4CPRNs/AjwQPX4D8DSQBg4CLwFu9NxfA+8kXFfwReD96/3+Vvla/XPgPwN/Ev1t12jxNfp94Ceixylgm12npuuzBzgLZKO/Pwf8mF0jNQu+G6jq66p6MnqcB54j/BB+kPDLSvT7f44efxD4rKpWVPUs8CLwdhHZDQyq6uMafvo+1fCaTY+I7AW+H/hkw2a7Rg2IyCBwJ/C7AKpaVdVp7DotJAFkRSQB9BGuudny18gEvsuIyAHgrcBfATtV9XUIbwJAXB1sD3Cu4WXno217oscLt/cKvwX8CyBo2GbXqJlDwATwHyJX1idFpB+7TnVU9TXgN4BXgdeBGVX9c+wamcB3ExHJAf8V+DlVXaqGTbuyDh2Ve9iMiMgPAOOq+lSnL2mxraevUUQCuB34bVV9KzBH6G5ox5a7TpFv/YOE7pbrgH4R+dGlXtJiW09eIxP4LiEiSUJx/wNV/UK0+VI0DST6PR5tb1fW4Xz0eOH2XuC7gQ+IyMvAZ4F3i8insWu0kPPAeVX9q+jvhwkF367TPO8FzqrqhKrWgC8Afwu7Ribw3SCKvP8u8Jyq/puGpx4B/lH0+B8Bf9yw/UdEJC0iB4HDwF9H08q8iLwjOub/2vCaTY2q/pKq7lXVA4RlLB5V1R/FrlETqnoROCciN0eb3kNYctuu0zyvAu8Qkb7ovb2HMO5l12i9o7y9+APcQTi1+xbwzejn+4AdwF8Cp6Pf2xte878TRvOfpyFyDxwFnome+3+JVh/30g9wjPksGrtGi6/PW4Ano8/THwHDdp0WXaP/AzgVvb//RJghs+WvkZUqMAzD6FHMRWMYhtGjmMAbhmH0KCbwhmEYPYoJvGEYRo9iAm8YhtGjmMAbhmH0KCbwxpYgWtTyZRH5poj8sIj8y2s4VmE1x2YY3aJrPVkNY4PxViCpqm+Bukj/X+s7JMPoLmbBG5sWEekXkf8uIk9HjR5+WETeFzXG+HrUsOFPRGQM+DTwlsiC/zxhadlvisgfLHH8PxKRp6JGEvcseO5fi8hJEflLERmNtr1FRJ4QkW+JyB+KyLCI3CIif93wugMi8q3o8dtE5KvROf4srptiGKuFCbyxmXkfcEFV36yqbwS+BPwO8HeBvw3sAlDVceAngK+p6ltU9QeBUvT4Hyxx/B9X1bcRLl+/T0R2RNv7gZOqejvwVeBj0fZPAb+oqrcB3wY+pqrPASkRORTt88PA56JidP8WuCs6x+8B/+qar4hhNGACb2xmvg28V0QeEJG/TVgu9qyqntawBsenr/H494nI08AThNUHD0fbA+C/RI8/DdwhIkPANlX9arT99wkbdUDYYeiHosc/HL32ZuCNwF+IyDeBX6a5kqFhXDPmgzc2Lar6goi8jbCQ2/8N/DmrVL9bRI4RlqF9p6oWReQEkGk3lGUO91+Az4vIF8Jh62kReRPwrKq+czXGaxitMAve2LSIyHVAUVU/TdjR528BB0XkhmiXu5d4eS1yk7RjCJiKxP0I8I6G5xzgrujx3we+rqozwFQ0kwD4h4TuG1T1JcAHfoV5y/95YFRE3hm9l6SI3LrsmzaMFWAWvLGZeRPw6yISADXgp4AR4L+LyGXg64RukFYcB74lIifb+OG/BPzTKCD6PKGbJmYOuFVEngJmCN0uENYc//ci0gecAf5xw2v+C/DrhG4kVLUqIncBD0bunQRhC8NnV3IBDGMprFyw0bNEbpYPq+oPrPdYDGM9MBeNYRhGj2IWvLGliVIf/7LFU+9R1cm1Ho9hrCYm8IZhGD2KuWgMwzB6FBN4wzCMHsUE3jAMo0cxgTcMw+hR/n/v9xM3Lt6WdAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.regplot(x=df[\"sqft_above\"],y=df[\"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": 57,
   "metadata": {
    "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": 57,
     "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": 58,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.00046769430149007363"
      ]
     },
     "execution_count": 58,
     "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": 60,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.49285321790379316"
      ]
     },
     "execution_count": 60,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x = df[['sqft_living']]\n",
    "y = df.price\n",
    "lr = LinearRegression()\n",
    "lr.fit(x, y)\n",
    "lr.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": 62,
   "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": 81,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.65765288394285"
      ]
     },
     "execution_count": 81,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[features]\n",
    "y = df.price\n",
    "lr.fit(X, y)\n",
    "lr.score(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": 82,
   "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": 89,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.7473257069312564"
      ]
     },
     "execution_count": 89,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[features]\n",
    "y = df.price\n",
    "pipe = Pipeline(Input)\n",
    "pipe.fit(X, y)\n",
    "pipe.score(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": 90,
   "metadata": {
    "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": 107,
   "metadata": {
    "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": 108,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.linear_model import Ridge"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 110,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6478759163939112"
      ]
     },
     "execution_count": 110,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "rm = Ridge(alpha=0.1)\n",
    "rm.fit(x_train, y_train)\n",
    "rm.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": 113,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.7002744255607272"
      ]
     },
     "execution_count": 113,
     "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",
    "\n",
    "rr = Ridge(alpha = 0.1)\n",
    "rr.fit(X_train_pr, y_train)\n",
    "rr.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 3",
   "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.8.2"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}