House_Sales_in_King_Count_USA (1).ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<center>\n",
    "    <img src=\"https://cf-courses-data.s3.us.cloud-object-storage.appdomain.cloud/IBMDeveloperSkillsNetwork-DA0101EN-SkillsNetwork/labs/FinalModule_Coursera/images/IDSNlogo.png\" width=\"300\" alt=\"cognitiveclass.ai logo\"  />\n",
    "</center>\n",
    "\n",
    "<h1 align=center><font size = 5>Data Analysis with Python</font></h1>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# House Sales in King County, USA\n"
   ]
  },
  {
   "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.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "| Variable      | Description                                                                                                 |\n",
    "| ------------- | ----------------------------------------------------------------------------------------------------------- |\n",
    "| id            | A notation for a house                                                                                      |\n",
    "| date          | Date house was sold                                                                                         |\n",
    "| price         | Price is prediction target                                                                                  |\n",
    "| bedrooms      | Number of bedrooms                                                                                          |\n",
    "| bathrooms     | Number of bathrooms                                                                                         |\n",
    "| sqft_living   | Square footage of the home                                                                                  |\n",
    "| sqft_lot      | Square footage of the lot                                                                                   |\n",
    "| floors        | Total floors (levels) in house                                                                              |\n",
    "| waterfront    | House which has a view to a waterfront                                                                      |\n",
    "| view          | Has been viewed                                                                                             |\n",
    "| condition     | How good the condition is overall                                                                           |\n",
    "| grade         | overall grade given to the housing unit, based on King County grading system                                |\n",
    "| sqft_above    | Square footage of house apart from basement                                                                 |\n",
    "| sqft_basement | Square footage of the basement                                                                              |\n",
    "| yr_built      | Built Year                                                                                                  |\n",
    "| yr_renovated  | Year when house was renovated                                                                               |\n",
    "| zipcode       | Zip code                                                                                                    |\n",
    "| lat           | Latitude coordinate                                                                                         |\n",
    "| long          | Longitude coordinate                                                                                        |\n",
    "| sqft_living15 | Living room area in 2015(implies-- some renovations) This might or might not have affected the lotsize area |\n",
    "| sqft_lot15    | LotSize area in 2015(implies-- some renovations)                                                            |\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You will require the following libraries:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "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\n",
    "\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module 1: Importing Data Sets\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Load the csv:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [],
   "source": [
    "file_name='https://cf-courses-data.s3.us.cloud-object-storage.appdomain.cloud/IBMDeveloperSkillsNetwork-DA0101EN-SkillsNetwork/labs/FinalModule_Coursera/data/kc_house_data_NaN.csv'\n",
    "df=pd.read_csv(file_name)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We use the method <code>head</code> to display the first 5 columns of the dataframe.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "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": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 1\n",
    "\n",
    "Display the data types of each column using the function dtypes, then take a screenshot and submit it, include your code in the image.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "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\n"
     ]
    }
   ],
   "source": [
    "print(df.dtypes)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We use the method describe to obtain a statistical summary of the dataframe.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "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": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module 2: Data Wrangling\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 2\n",
    "\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>\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "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": 39,
     "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>\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "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": [
    "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>\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean=df['bedrooms'].mean()\n",
    "df['bedrooms'].replace(np.nan,mean, inplace=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "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>\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean=df['bathrooms'].mean()\n",
    "df['bathrooms'].replace(np.nan,mean, inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "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\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 3\n",
    "\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": 44,
   "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": 44,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df['floors'].value_counts().to_frame()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 4\n",
    "\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.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='waterfront', ylabel='price'>"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAERCAYAAABxZrw0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAVG0lEQVR4nO3dfZDdVX3H8c9nd0ESHsayrAwkkYhBEEGe1ocWqzwszmorznRUZKqslpq21iS0RUetU8cZ63SmndokVcuOpWxGxQIFBx1YzRaQBwHdhEAIwXGLAbMoXFfkKTGwu9/+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",
    "\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.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='sqft_above', ylabel='price'>"
      ]
     },
     "execution_count": 46,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAESCAYAAAD38s6aAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAABBW0lEQVR4nO29fXhkd3Xn+Tn31quk0kt3S91ud3faDbYbDDbYHmInHqcHnI0DCZ7dhyR4k3kLxM4uAyGzZCAzhOw6Ozt4mZ0MTNiJPSSTsCQmxIHAECABPE1D1ibYBhsbt19Q2+53vbRaUqne7z3zx723VFWqkkpqlVQqnc/ztFW6davur65V5557zvecI6qKYRiG0Xs4m70AwzAMozOYgTcMw+hRzMAbhmH0KGbgDcMwehQz8IZhGD2KGXjDMIwepesMvIj8oYhMiMhTbe7/8yLyAxF5WkT+tNPrMwzD2CpIt+ngReRWIAt8UlVfs8K+VwKfAd6oqjMiMqaqExuxTsMwjG6n6zx4VT0GXKjdJiKvEJGviMhjIvJNETkcPvUrwMdVdSZ8rRl3wzCMkK4z8C24H3i3qt4AvA/4f8PtVwFXicjfisgjInL7pq3QMAyjy4ht9gJWQkQGgB8D/lxEos3J8GcMuBI4AuwDvikir1HVixu8TMMwjK6j6w08wV3GRVV9XZPnTgGPqGoZOCEizxIY/O9s4PoMwzC6kq4P0ajqHIHx/jkACbgufPovgX8Qbt9FELIZ34x1GoZhdBtdZ+BF5AHgYeBqETklIu8AfhF4h4g8ATwN3BHu/tfAtIj8APhvwG+o6vRmrNswDKPb6DqZpGEYhrE+dNSDF5FfDwuQnhKRB0Qk1cnjGYZhGIt0zIMXkcuBbwGvVtW8iHwG+JKq/lGr1+zatUsPHjzYkfUYhmH0Io899tiUqo42e67TKpoYkBaRMtAHnFlu54MHD/Loo492eEmGYRi9g4i81Oq5joVoVPU08O+Al4GzwKyq/k2njmcYhmHU0zEDLyIjBGqXK4C9QL+I/FKT/e4SkUdF5NHJyclOLccwDGPb0ckk623ACVWdDAuRPktQkVqHqt6vqjeq6o2jo03DSIZhGMYa6KSBfxm4SUT6JOgx8CbgmQ4ezzAMw6ihkzH4bwMPAo8D3w+PdX+njmcYhmHU01EVjar+NvDbnTyGYRjGenH0+AT3HRvn5EyO/SN93H3rIY4cHtvsZa2ZrmtVYBiGsRkcPT7Bh77wNBPzBYbTcSbmC3zoC09z9PjWHTNhBt4wDAO479g4cVfoS8QQCX7GXeG+Y1u3f6EZeMMwDODkTI503K3blo67nJrJbdKKLh0z8IZhGMD+kT7yZa9uW77ssW+kb5NWdOmYgTcMwwDuvvUQZU/JlSqoBj/LnnL3rYc2e2lrxgy8YRgGcOTwGPe89RrGMilm82XGMinuees1W1pFsxVG9hmGYWwIRw6PbWmD3oh58IZhGD2KGXjDMIwexQy8YRhGj2IG3jAMo0cxA28YhtGjmIE3DMPoUczAG4Zh9Chm4A3DMHoUM/CGYRg9ilWyGobRdfTa4I3NomMevIhcLSLfq/k3JyLv7dTxDMPoDXpx8MZm0cmZrM+q6utU9XXADUAO+FynjmcYRm/Qi4M3NouNisG/Cfihqr60QcczDGOL0ouDNzaLjYrBvx14oNkTInIXcBfAgQMHNmg5hmF0K/tH+piYL9CXWDRP3TZ4Y6vkCDruwYtIAngr8OfNnlfV+1X1RlW9cXR0tNPLMQyjy+n2wRtbKUewESGanwYeV9XzG3AswzC2ON0+eGMr5Qg2IkRzJy3CM4ZhGM3o5sEbJ2dyDKfjddu6NUfQUQ9eRPqAnwQ+28njGIZhbBRbaTh3Rw28quZUdaeqznbyOIZhGBtFt+cIarFWBYZhGKug23MEtVirAsMwjFXSzTmCWsyDNwzD6FHMwBuGYfQoZuANwzB6FDPwhmEYPYoZeMMwjB7FDLxhGEaPYgbeMAyjRzEDbxiG0aOYgTcMw+hRzMAbhmH0KGbgDcMwehQz8IZhGD2KGXjDMIwexQy8YRhGj2IG3jAMo0fp9Mi+YRF5UESOi8gzInJzJ49nGIZhLNLpgR8fBb6iqm8TkQTQfUMLDcMwepSOGXgRGQRuBf4pgKqWgFKnjmcYhmHU08kQzSFgEvgvIvJdEfmEiPQ37iQid4nIoyLy6OTkZAeXYxiGsb3opIGPAdcD/0lVXw8sAB9o3ElV71fVG1X1xtHR0Q4uxzAMY3vRSQN/Cjilqt8Of3+QwOAbhmEYG0DHDLyqngNOisjV4aY3AT/o1PEMwzCMejqtonk38CehgmYc+GcdPp5hGIYR0lEDr6rfA27s5DEMwzCM5lglq2EYRo9iBt4wDKNH6XQM3jCMDebo8QnuOzbOyZkc+0f6uPvWQxw5PLbZyzI2AfPgDaOHOHp8gg994Wkm5gsMp+NMzBf40Bee5ujxic1emrEJmIE3jB7ivmPjxF2hLxFDJPgZd4X7jo1v9tKMTcAMvGH0ECdncqTjbt22dNzl1Exuk1ZkbCZm4A2jh9g/0ke+7NVty5c99o1YI9ftiBl4w+gh7r71EGVPyZUqqAY/y55y962HNntpxiZgKhrD6CGOHB7jHoJY/KmZHPvWqKIxJU5vYAbeMFZJtxu/I4fHLmk9kRIn7kqdEuee8L2NrYOFaAxjFWwHGaIpcXoHM/CGsQq2g/EzJU7vYAbeMFbBdjB+psTpHczAG8Yq2A7Gz5Q4vYMZeMNYBd1i/I4en+DO+x/hlnsf4s77H1nXHMCRw2Pc89ZrGMukmM2XGcukuOet11iCdQsiqrrZa6hy44036qOPPrrZyzCMZYlUNJciQ7zU40cql3TcJV/2KHtqRnibIiKPqWrTuRsmkzSMVXKpMsRLpTbRC9CXiJErVbjv2LgZeKOOjhp4EXkRmAc8oNLqKmMYRvucnMkxnI7Xbeu1RK+xPmyEB/8PVHVqA45jGNuC/SN9TMwXqh489F6i11gfLMlqGFuMbkn0Gt1Ppw28An8jIo+JyF3NdhCRu0TkURF5dHJyssPLMYytj6lcjHbpqIpGRPaq6hkRGQO+CrxbVY+12t9UNIZhGKtjORVNRz14VT0T/pwAPge8oZPHMwzDMBbpmIEXkX4RyUSPgf8BeKpTxzMMwzDq6aSKZjfwORGJjvOnqvqVDh7PMAzDqKFjBl5Vx4HrOvX+htFpur3vu2GshMkkDaMJ26Hvu9H7mIE3jCZsh77vRu9jBt4wmrAd+r4bvY8ZeMNownbo+270PtZN0jCacPeth/jQF54mV6rUteTdqu0AahPGAwkXEWG+WLHkcY9jHrxhNKGX2gHUJoxdgRcmF3h+IosrWPK4xzEP3jBasNl939eL2oTx+GQWVwQEprIlDo0OWC/5HsY8eMPocWoTxiXPRwREgsdgyeNexgy8YfQ4tQnjhOugCqrBY7DkcS9jBt4wepza/vG7BhJ4qni+smsgYb3kexyLwRtGj3Pk8Bj3QHVQ+CtH+xERssUKY5mUqWh6GDPwhrFJbGSvm15JGBurw0I0hrEJWK8bYyNo28CLyI+IyG3h43TU690wjNVjvW6MjaAtAy8ivwI8CNwXbtoH/GWH1mQYPY/1ujE2gnY9+HcBPw7MAajq84AF9AxjjVivG2MjaNfAF1W1FP0iIjGgrWndIuKKyHdF5ItrWaBh9CK10kVVNbmi0RHaVdF8Q0T+FZAWkZ8E/lfgv7b52l8DngEG17A+w+hJGqWL+xpUNDZNylgPRHVlR1xEHOAdBIOzBfhr4BO6wotFZB/wx8C/Af6Fqv7McvvfeOON+uijj7a5dMPoTSKFTdyVuk6WW7XZmdFZROQxVb2x2XPtevBp4A9V9T+Hb+iG21bKCP0H4F8CprgxjDapVdgA9CVi1hDMWBPtxuC/TmDQI9LA15Z7gYj8DDChqo+tsN9dIvKoiDw6OTnZ5nIMo3eJFDbzhTLjk1mOn5vj7MU8z0/Mb/bSjC1GuwY+parZ6Jfw8Urp/h8H3ioiLwKfBt4oIp9q3ElV71fVG1X1xtHR0TaXYxi9y/6RPqYXipy5WKDiK64jlH1lvlCxQihjVbRr4BdE5ProFxG5Acgv9wJV/U1V3aeqB4G3Aw+p6i+teaWGsQqOHp/gzvsf4ZZ7H+LO+x+5ZMO43u+3HHffeogLC2UURQANuvoy0he3QihjVbQbg38v8Ocicib8/TLgFzqyIsO4RGqTlLVtAO6BNcWw1/v9VuLI4TEGki6Fsk/J80m4DrsGUmRSMSuEMlZFWwZeVb8jIoeBqwlUNMdVtdzuQVT1KHB0LQs0tiabKfNb7yTlZiQ9r9o9yMR8oXpMgFypYoVQxqpYNkQjIm8Mf/5PwM8CVwFXAj8bbjOMJWx2I631bgOwGW0FrBDKWA9W8uB/AniIwLg3osBn131FxpZnvT3e6G7gufNzlD0lEXO4cizT8q5g/0jfEu93LW0AouNOzheZmi+yZyhFJhVf8/s1vu9ydzcrFUIZRjssa+BV9bfDIqcvq+pnNmhNxhbn5EyO4XRgCOcLZSbnixQrHqdm8hw9PrEqIxXdDZQqHnOFCgD5kseL09mWcfC7bz3Eh77wNLlSpa5QaDXeb23cfc9gktMXC5yayXP5sBJznTV706uJ51sPd+NSWTEGr6q+iPxzwAy80RaRB+35ypmLBUTAEUEE3vfgE4wOJJkvVtqKzUd3A9PZCg6C4wi+KnP5CnuGYk3vCtbD+228CwHh/HyBc3NFrj8wsmZv2oqYjI2kXRXNV0XkfcCfAQvRRlW90JFVGVuayIOemCsACiookEnGmMmVyRYqvHJsoC01SnQ3UPJ8XBEARKDk+cvGwS/V+629CwEYTMfJpGLM5ss8cNdN6/a+YG2Cjc7Rrg7+lwkajH0DeLTmn2Es4cjhMe556zUoQaIm5gp7h9LkSh6OgKfa9pCLqK1uwnWIOh+pQsJ1Otpet1PtfK1NsLGRtGvgXw18HHgC+B7wH4FrOrQmowc4cniM6w+McGBHH4dGBxgMvXAIjHPESt5rpCbJpGL4KBXfx/eVwXSso6qSTqlYbj60g1MzeX5wdpYfTswzlS2YOsboGO2GaP6YYNjHx8Lf7wy3/XwnFmV0D6vVs9fun0kGIQ0IDLnrCBVPGc0kq/vXeq+tjhXF0yveHKVQRXNw50DHh1Svt4rl6PEJHnz8NCN9ceYLFYoVnwsLZd515IDF342O0G674CdU9bqVtl0q1i64u2jWtnY2X26ZJG22/1y+zM7+BAslj/6Ey/RCicF0fEkbXKDnW+Teef8jTYuXxjKpS4rrG9ub5doFtxui+a6IVP8CReRHgb9dj8UZ3UvjYOiKp1zMlTkxtdC0gKnZIOnBdJyR/iTffP8b+cqv/wQfedt1jGVSzObLjGVSVQO+HYZQ2xxWY6NpN0Tzo8A/FpGXw98PAM+IyPcBVdVrO7I645K5lJYBjYqPqWxxSZK0VuLXjkKklbqlE+qS2s8+kHARkbblmZ1gvQqwGrHpT3YOWtGuB387cAVBZetPhI/fDPwMzatcjS7gUlsGNCo+VkqSXopCZL3VJbWf3RV4YXKB5yeyuMKGt06I6ETidrPbQnQDdg5a026zsZc6vRBj/VlNUU0zD6ixInSlJOlKFaTLeVnrUX3a6rOPT2YDDb3AVLbEodGBZYuLVtsaodXrGz9nJxK3Vjhl52A52g3RGF3Cam5F2w17tCyff+s13PPWa6oG6eCOPqYXSriOoKpLjPByBmylEv31Nn61n71aJCWLdyGtwj9raY3Q1rms+ZzraXSscMrOwXKYgd9CrLYvebsx3+U8oAfuuqnuvaMLTCsj3MqAteNlrafxq/3sCdeh4gVqMVeE8ckshYpHfyK2pDfOWlojrPZzriediutvJewctMYM/BaiHePRmFicq9GhRx73zYd2cOf9j1TvAp6fmGfPYKruWK08oLUa4UvxstaSQKsN+ewaSHD6YgH1AyPvqeCI0Jdwl1wgV9saoXFtz52f47KhxfHFc/kyU9kiL07nuPP+R9Y9+bfeoa2tiJ2D1rSbZDW6gJVkdo3JprKvKEFSNJIlvu36y3nw8dN1Can5QoXphWLd+663BxTNGY2GSI9PZpleKK54jLUm0KJ2CWOZFL7CK0f7ScYDJU3Cddg7lGY0k1oixVxNa4Rma8sWPaaywbmcy5c5M5un5PkkXelI8q/2czZKT7cLdg5a0zEPXkRSwDEgGR7nQVX97U4dbzuw0q1oMw8fYLgvwZffeysQFNs07jPSF+fCQpm+RKwjHtDHvvYcj790gWIYJnGAEj4T8yXu/Hs7ln3dRx96Ac9XHIFd/Ql2D6WbhjyWS2xG3HLvQwyn40jomcPSu4jIG8ykYkwvlPB9BYXB/viSu5+5fJm+hMtQOlU9lzv6g3PZn4xVDb0gjA2mOhausbbCdg5a0UkPvgi8Max2fR1we22xlLF6VpLZtVNI02yfXQNJMkm3Ix5QZKQj4w7gA6gyOpDg4fHmDUlrjTuArzCRLXF+Nr/kM7Xr5bcjxYy8wSt2DTCUCi54Q31xDu4cWHL3s1CqML1QqobBAHb2J8mkYoxlUhQqPnFH2Du8OChko5N/Gzks3Og+OubBa9ADIRv+Gg//rdwXwWjJSkqTdpJNrfa5cvfgiuXyy8XCWz33iW+dwJHQqANIEPLwNLiwtDJ20eu8hu0T2RKZdLzuM7XKTdz7leN1a7r50A4efPx0Xax2Ll8m7gi33PtQU8+/9nM9/a1Z+pOLHnsq5lLyfKayRQbTi5OerhzL8MBdNzVtTbCRyb+NHhZudB8dTbKKiAs8BrwS+LiqfrvJPncBdwEcOHCgk8vpCZa7FW0n2bTWhFStsXAFvvvyDO/45He4cnSAN7/2Mh58/HRTQ7JQ8og5QaKytu2Rr0uNXa0xnStUWt5ezubL/NZbXl19zeMvz+D5PsmYy2gmSSYVp+L5vDid4+DOvup6/+7FaS7LJEkkEszmywwkYyhQ9rVu3W87dZGHxy/w3Pk5skWPHf1xdvYnOTdbIF/ySMaCOyDPV8qeUvY85vIlihWfmVyZ2XyZO+9/pOkFpVPJv2YXWNOHG201G7vkg4gMA58D3q2qT7Xaz5qNXToryRhXu09kMC7mSpQ8v25Kk6riOIIq7OiPs2tgUYkTNdF6+sws+bKHIJRDDboCjsC+kb5qKKixUdkPzs4RRmcQ6m/9XrUnw5ffe2tdaMb3g/YJqrB3OMW52QII7BlMLVlvFIK679j4Eg97KlvgwkKZfSNpzs0WKHk+QhBmmZwvVtU1vgYXLc/38TW4eMVcYddAgp39yaoxf9v1l/Pw+IWOzlVt1uSt7CkLxTKXDaXrcg6qymy+zDff/8Z1XYOxeSzXbGxDZJKqelFEjhK0PGhp4I1Lp51k00r7NLu1f3E6x77hFFPZUnUEnxJ4saowmyvXGfgo1vzOW67gow+9gCNKzIFyGKvZO1Qf52/0Nnf1J5jIlqrvFxn5sYEE2WKl7jW7MynOzOYRBVDOzRYo+z77htNMzheXrDdSzjSTbs7myni+0peIVY25ApPzRXYNJDkzm6dQ8Um4AiqICPuH05yfL4BSPQeRt/zw+IWOd4ps5amXvaAYzfTh25dOqmhGgXJo3NPAbcC9nTqesX40MxhxVzg/V8RTxXUCjzCSDwIUKs2Tl++57SogiKkvlDwGUy7vvOUKrt03zH3Hxvng559qqh/fPZRmJl+m7GnV4x8bSDIQJjBhUbMuiWA9U9kipXD/K0cHKPsaGOmG9UYXn2b5iKLnkwpDMFGBlDiBDn4wHadY8Tg/X0SBuCvsGkgxmI5zZjZP493wRiVUW9UYJGJONSlv+vDtSSc9+MuAPw7j8A7wGVX9YgePZ6wTJ2dyuALjk1lKnk/CdcgkXS7kysQdJwiHAD7KroEUxYpHJdfakLzntquqhh6a3yFE+vHRzOJdwM7+BDO5IFySjrtMZYucmslXY9wDCbfqoQ6m4wym49XQUJRrEKBQrqZ4KXnwwmSWgzv6muYjYo5DJhXeRYQeOz7EHSFXqpCIuRzenaHk+XUXBtcJPPpa2vWWL7UTYsvEedg/Zz173xhbi06qaJ4EXt+p9zc6x0DC5YXJBVwR3LAP/EyuwmWZJJl0gucmssRdGE7GOD8fjJwbScc4e7FA0fPpTwReeitD0uwOYUd/nMn5IhfzQYhEAF+VdNxlcr6IqlLytJrwnJgvMJcvV2PzjZW69x0bZ2ahSNmr96oVKFd8pheC8E9tr519I33ccd3eamI0k4qxsxJcZPqSseqFA1hyYRhIxhC4pOT1WpUuyyXOe00fbm2BV4e1KjCWUE3KSfgPQCGTTlSTmx/+8jM8P5kl7jjs7ItzIRdowfcNp4i5Dg8+fppr9w03LUb6uxcvkHSD4p9IH55wnSCpquD7iqfgCuzoTxBzHU7N5OsSudHFIeE6DPclqgY6Uq7EXcHzlZgrVSMv4X/irsNgOl7Xayda22ceO0UmGasmI6/YNcCHmxiRRrlqpOpZrbe8HkqXTnSp7EZM9lnPQrHCU6dnl93HDLyxhPlihcvDhGoUotkzmKwmN6MJTAfDZOT4ZDYMUbRux1v75UzFHEqez5mLBfYOQyYV5/x8kYTrcOXuDOOTWSphBWn0fmXP59xskalsiYTrMJpJMhDOfI2qdKG+UrfsB/mCshdUwiZjLori+dq0xUNkOCIP+HfueE1Lw9HKM16toVmvToi95qk3YzvLPosVj+Nn53ny1EWeODXLk6cu8sJEtqo0a4UZ+C3Opdyytnrt/pE+XpzO1u1b8nwO7hyo/r6kHW+YyGxsxxsd4/GXZxCB3ZlUNbatKGdmcjiuQ6Ec9GuZL5SXvN9cvowG9j7oSR9KNXcOxOvWFK0ryh9UPA0uFFD9IjTrK7OZhsM6IbbPdmkL7PnKCxNZnjh1kSdPXeTJU7M8c3ZuSbixHczAb2Eu5Za18bUvTme5+1OPMZB0Sccczs8VcR3BESiUPRZKHtlCpdoRcUk73tDjjlQ1+XIwZPs3HnyC+UKFYiUw/C9dyOGEYRLP8ykrpBxIxhwqns9L08GXteIrrgjJmMNUNpA6oqB+oD/3US4slPm3/2N9jDuTjPH8RJYoylQrbKn4ftO+Ms1CRu0YjvWIB1snxPbpxYuhqvLSdC405oFn/tTpuSUtNSIuG0px7b4hrt03zHX7hnnt5UMML6NNNAPfhbRrOGo9z/lCmcn5IoWKx3s+/V0+9vbXt923fL5QZjpbRlEKZZ+Fkgca9E4v+4qvgUzR87V6EbnhwBCPvzxDxfeJSeBViwh7BpLVYqGK5+MpxBypK1bSMNZeUYg7cOXuDOdm80zW6N6D55WRhMNMroIrwkh/nIWSVw0bpePBxaS29fF8oYyqUvEbPzHEHYe+pMvBnQN1sfpmIaOVDMelXlxr//9uRDFUL9ALF8PzcwWeOBkY88ioz9b0MqplpC8eGvLAoF+7f6gqEW6XDalkbRerZG1dldis+VfUHTFbrFSrNUHxfNg7nF62YVhtZ8Uo5i0EA7Uh8CyiCk0INN8K7B1Kc/pivhr6EAKP2nWEoVSM2XyFsq/Ea4qaavK0VZIxh2LFJxlzuCqMuxfLHj5BOCUVc/B8H8dxSMYc+pPukkrZhOuwUPLqztWL07lqg7JG0nGHnQNJ9o/0cWomx+R8oOt3oHohg+DYmVSMj7ztupbnr1mfmcn5ArmSx2A63vLCvJr/v6thu6hL2qnC7hYu5kpVrzyKm5+fKzbddyAZ4zWXD3LdvuHAmO8bYt9IfRVyKza9ktVon9XEg6Nb1tpqTd+HZEyq1Zqt/vhrb3ejmLf6QYjF85VCQ7yv5AVG+9TFfJ0BVQID3p9wmS96uI4QTMcLlPKN+0WvjAQ6xYrP+GR2sQc7gcd/aHSgqmT5nTte09Rzizu65Fy5TRqUReTLPsPpOCemspybKxJzIOYGoSFfa9YnK3fFa4wHz+XLQXthVQ7s6Gvp0Xci3r+d1CW1lc8nZ3LVXv6b/TkjRcv3T89WjXkUbmwk4Tq8au8gr4tCLfuHOLRrAMdZ2ZivFjPwXUSrxlmt4sHRLWuh4hFzAuMeFR+tFEOuvd2NO0KxEhi5sk9dDxiRxd8Dj3yp6VOF2Xyl4Rld8ij6GXeC93Qc8P2gCtZXKFR8BEilAmVOseLTlwiqShv16nffeogPfv6pJUk3Vrgjffb8fPUC5fnBnUmtw9+fcFuqgGo95NoiK6Da+z0VC4aKtDLcnUgUbid1STdczFajaHEErtqdqYubX70nQyK2MbOWzMB3CdEfru/7eH7QhTE3nWN0ILGkPW5EpH9+z6e/S67kkYwtls7nSpVlY8i12unJ+SK5sh8mP4ViJfhLjYVa8pVYbZAvDALhiDCQcqvDraPnZgsVXAcchP6kWx0A/sBdN1UN7Qc//xRz+SDOX1v92sp7hzBUpDWKGqAcXlyi30czSaD5pKxao9JYZFUML07R6xvfI6ITicLtoi6Bjb+Yeb7y/MR8NdSykqLlil39XLtviNdePsR1+4e5Zu9g3f/rjcYMfJdw37FxShUPH0GrJjDof+5DtZCmlsjYQRA7z5eCcv9ixSMRc1dMPkXa6dt/9xsUyguUPaVUiYZTB7eSh/YM8Ny5uWqPl+Vo7PrYap+KD/0Jh76Ey1S21PQ1qnD5jjSZVHCx+vCXn6krrhpMxSiUPeYKFaazRfYOp4m5Dl6T5Gr1PQli7bXUfk+D2PtiX/eVJmXVFln1JVz6k2719Y3vEdGJRGEvqkta0cmLmary8oVc4JWHidCnzsySK7WvaBnqizfdd7MwA98lnJzJMV+oEBPBdYWK71e136MDyZbJurLnkS95wVANhVzJo+Ir7zpyoC2P5ujxCZ6fzFYlicEwDsX3lWLFR1WrSpiRvnhdh8dGojh79LgZAiTCJOZkC+MOgYGPjGXQ2z1P3HWqyp7adZR9ODmTZ8cavlxxV4KEMoFKR1WXGN1WRqW2yCr6/7GS4e5E1WkvqEvaZT0vZo2Klu+fnuVirnOKlmZ0OjluBr5L2D/Sx7nZAjFHcERwHRdfAxnL+NTCkolDkVc5na3gOEJMHHwJSvMzqRif+NYJPvPYqRX/aO47Nk7cCZKbIkGC1Kv4KIEW/akzc8Cid550pW78Xi0xRwIVzjJuvE8QlzzXQk3QjLOzhSBW77d2z32FqYXmX87l8FVxHYd3H3lFS6liO0ZlNYZ7vatOV3Psra62WevFbDWKlv6Ey2vCEMtrLx/iun3D7N/RnqJlNWxEPsFkkl3C0eMT3P2px0KDEwyuCHqtK/GYwytHB5heKHJhocxA0iVb9NgzmOTMbCFQrhB4ouVwQLWvyqv2DNbJ8IAlX+4Pfv4pXIGzs0V8Xaz8XCtu+B1Yrugu4UBpmVAKQMIVrtqdYSpbrCpe/Jr4+XrghGEoBa4/MFI1do1GMNLMlypetWjLdYR3HXlFXZfMbqdTEs2NZiWpZKRoqdWav3xheUXLddVQyxCHRgeqldSdpJnUNuqGupoZAsvJJM3AdxEf+9pzfPzoD6n4PsmwOtRXuHw4jQicuVhA0UDKqErFU2KOhP3SBd9Xyn5gfBKuw6HRoIw/UspMLZTIFoNBEL6vdaqWvkSMuUKjEqY9amPv8dCLv1RDvH8k6A0/my9TCCdC+eEFbD2Ivr6uI1weNkiLJjBFBVC1RvCGA0N86anz1f83Q31BNezoQJLJbJFSxSfuClftHuxar3i9DEo3Uar4HD83xxMnu1/R0khtLUrEWiZumQ5+i/Ce266qDsI4NZNjYr7I5UNJBtNxxieD8nsHoeT51YIjX4OJSr7USEF0qZrj+Lk5IrPWqIwp+6FyRZb3vFtR+5JyWDB1KezoC1RDD9x1E7fc+xAjfXHOzBZwpJ00bnsowYVt73C6OjA7V6rwiW+dYDSTXKLS+PrxSfaNpKvb5/JlJufz1T45COTLcGIq2/Q2uxtCI1tdbeP5yg8ns9W4+ZOnLvLM2flq/6NGrtjVz2svH+Lafd2haGlkI5Lj3fNpDaA+Pht5XMCSYqTAKCnn5ooMpmOUwjFyJU+bqjm8sC1AsxBMLGzgFbUkWMlJXsnMXooJTsUc5gplHn95hqPHJ6pfgr1D6UAhtIwOsp21R/upwtV7Buu8p3TcZaHkcSDu1u3fbPtUtogjUPaURMyp3kHNFyrsGYq17KS5mYVIW0lt04uKlkY2IjluBr6LaSxGisITUcl+yQtK/VNxtzq9B5YOoyh7SkwCTXszAxh59EG4Is3LF3LLqmCWI+E6lD1/1UZewuNH4SaR4HPUhkyu2NXPCxNZyl6g2a/4ixcT12FZiWTdsSQoDW82r7S/oYCp1fZgGPfi+0U/S56/xCvulkKkblbbrEbRMtygaLlu3xBjg5euaNloNqKPfydnsu4HPgnsIRBP3K+qH+3U8bqZtd6e1/4BzObLzBcqjPTFyaRiTGULTMyXGB1IVL3C9z34BOm4y9nZfCBtFOGVYwP81lsO874Hn2CqhcSx2j5AguKmXQNxprLl5vp0Whv5hOuEyangTqIdamWVTnhr4CmIr5y5mOf3vzHOjr4YL10s4fmKE7ZC2DmQZNdAknzZYzZfJlf0KKjXlgfv+TCUjjM5V6DkadAwLRzV985brqhOdKo1gm86PMqXnjqP5yvJmAOhtDKSloo0b0UM3RMa6ZbBIGtRtNR6551QtGwWne7j30kPvgL8b6r6uIhkgMdE5Kuq+oMOHrPrqNWrz+bKnJ3N8/jLM20pMGovDFeOZbj50I6qlG+h6JFJuswXKkwvlHBEKIXyxoQr1XDFZLbIk6cuMr2Mfj1i90CCsUyK2VyJRMypvl867nD7Nbv5ytPnyYcdxOJOkPSsteMVz1+Ud7ZBo2Y+0NuD+griIPjkyh6FWa/aH8bzFdeBifki09kSjiMc2tXPhBYYzSQ4OZNnBaUmjgRKi2zJw5VwglX4/teGzZ5qjWCkotnRH2c2V6ZQ8XBESIShrehi5jowkkou8Yq7KTSy0YNBcqUKT52eqzPmK/Vo2QxFS6+yYSoaEfk88Huq+tVW+/SiiubO+x/hxeks09ly0IQrbLvriHDfL93Q8svWTNI2my8zOpBkvljh3GwBX5W44yASNO2qLZ2HYPSdOIEhzZeXj18kYw5jmSS/c8dr+I0Hn+DCQqnOeGeSLsm4S39YfZoreU2NaPRVjHq+q0Kx3dhJzXvEXYeyvxiDaVZEFUkp82WPUzP56gVpxfcP7xQcgVTcrVMbNVOUNFOfTGULTGVL1YlRGia7D+/J8P7bDy9JsG60PHGtd42XkgyuKlpq4ubPT8y3vKtKuEJfMkbMEWKOw/95xzXcds2e1XxMgy5Q0YjIQYIB3N9u8txdwF0ABw4c2IjlbCgnZ3LM5srVbo8QaMUrvi4bg43ithVPOTG7QKHi4/vKbK7E3uF0NVm6nPGMLibtTIIpVoJh2fcdG2/aPmC+GAz9iDlBaOSlFrriWm/cq/isRYCmLE6Gqusj37BfyVMm5goslLzqQJGVEBb7kXlh5e9cvsxgunVTt2YhltlcoJ7ZnUkxlS1SqPhA0LekscPhRodG1prUXc3rVqtoObizr9oG93PfDUJgA8nFc5orVfiDv33RDPw603EDLyIDwF8A71XVucbnVfV+4H4IPPhOr2ej2T/Sx9nZoMw+QjXwmBuNSa33NDlfZDDlcjFfwUGquvWyDy9dyC97TM8P4sqRKqZdceGZi3nGpxZa7huEfEordmyse03bezZnpSMt1zqh3fc7Mxucz5grTcMmmWSMFyayeKrVebBFL4jbn5nNV4vSIDhHz5yd5e5PPUYmFasmvzcyNLLWpG6r1/3+N37IFaP9bSta9gymqtLEqPHWcF+i+vwf/f8vdkVOYjvQUQMvInEC4/4nqvrZTh6rW7n71kNhC2DFDRNxPkrcdZjNl6stCGonDA2n40zNF5nKloOxea6gy7VJbKDWi1pNXdB8sXnYpZaxTIKzs+23Geh2lEDqeHY2z56h9BJFydHjE0xmi1TCCuGy53NqJh/KIn0ccSirXxc+ms1XiLlCrljZFElkq6TuU6dnuPZ//2sWSoEq6J23XFGXB4rm2b4wEXjiDoLrCCemFviJjxxteqy1KFq6KSfR63RSRSPAHwDPqOq/79Rxup0jh8d415FXhBWqgQIj7jrMFz3GMosKmI8f/SE7+uMMpYMvx56hFC9O56j4SszdmBubdo4yOd87xr2WUljF2miE7zs2TswRXKGaTI27wthAgrPzpaD3Z8OJUwLJZ9nXJd7zRhQ8NTOgp2ZyzBd94q4ScwKD+tGHXiBf9rj5Fbt48tRFLubKZIuLrZs96iuH16tHSzfLNXuNjiVZReQW4JvA91m8U/9XqvqlVq/pxSRrRG3/jNl8eckIumfOzpFwhdEwplvyfCp+YDzWr37TWI59Q0m+9Zu31W274Xf+hrlCECaTmjuwoVSM0UyKE1ML1fh/3A1mu0YJ5laTqTqdbG2W1B2fWkCAmONUq5+9Fb77kdneM5jkk+/40XVVtGyl0XvdjvWi6TKa9aAYn8yyUPKqxnyzjHpMgmHYvUg753Qg6fDay0eqktRHxqeBoNo3FuZRihUPEWE4HSNbDMId84UKCFQ8DeoJHIe9w6lqP/uotWwzNc5Ccfk5rs1Y6U7ga0+f46MPvcCpmRwxx2Ey2/rOyxG4cizDyZkcqZgTtpwOportGkjgK6vqjWJsLJuuojEWOXp8grl8mXOzhWCYdCIog18IE1abadwhMO69esfQzmfKFn0ee+kCD49PkwgHjUM0KCTIUfgKLsplYfuEmVyZdDhaUGM+JQ929McZSAbhmSj80DhicL5QZmo+UCwtN8e1kUa1y/m5PB/47Pd5y2svo+T5KypahEDRpSipuMujH7yNvkSsKge9fGRpMzJja2IGvoO0ajvbn3TJlzwKFa/qtdey2cZ1s4+/2USx9pKndf1tyn7g7brO4uzV0UwKEeq88NqCtNrww/5j9bHxyfkiCCRdZ9k5rrWoKh/9+vMUyh65knJqJl8NEf3B355Ysn+kaFkoVnh4fDpcv4Stl4WfevUY7/ijRzk5kyOTjDGbD9oDWGy8NzAD3yGOHp/gfQ8+QbZYwfOVqWyRR1+6wK6BBLsGUiRjbrVH9XY3qN1MbQQz8nxrZ682euEnprI8/vJMnUQSgmKp587PkS167OiPs7M/Wa2I3TXQeo7r+bkCf/LIS/zF46eZyhZDCezyfzHJmPCrt76CX7zpR+oULR/72nN84lsnqiqaNx0e5bGXZ6t3Avly4GzEHWE2X7bYeA9gMfgOcfvvfoMXJhdww8ZZFc+vxrb7Ey6jmSQvX8it6wALo/MMpmJ1CfLxySwlzyfhOuwaSFY19XFHuGw4zWy+jEC1kCoK6WRSQQfQvoRbHRpe8X0u5krEHIfDlw0u26NlJV61J1MdJ9iKXuwPvx2xGPwGURuSOX0xjyvgOE6dcQeCmHuLfhxGdyLA1bsH+MBPv4oPfeFpJucLzBcq1RDbcDow3g6COFQlkqcv5kGhPxnjxNRC0PZZhJ19Ca7bP8xffPd01YjXXutPztQXs6XjLumES1/cZTZfYn65vsnACxPZFT9TtzRBMzqHGfh1ojHxdWomT0XBEV1RjmZsDd782su479g4k/OFJb19LubLeH5Q6Rp1lYSgwtXzglg5BDH0iirPTmR5tokRFoKpQ9ftDwqHPvb159k1kMB1FiuhZ3IlhKAVRas7wEobf3NWcNT7mIFfJxrLvFMxh0LFX1NvdKP7SLhSnctaqDHu0RSsqFVBxfcBYSAZeO+VsMWE38ISJ1yHdMIlHXdxBPaN9PGZX725+vxfPXk2NMKLBr5Y8UnHHcYGg2K4ZrhtyNWt4Kj32ZxhhD3IyZlctYsjwO7BVDA5aBPXZKwfxbCx2cVc0Cc/KmGonYLlazCEpOIr0wslLiyUlm3bI8DVezIc2NHHaCbJjv4EZ2frQzN333qIsqfkShVUg5+uIwz1xcmk4sQbLLkQfKlfGXbIXI4jh8e4563XBC2i82XGMqktN4DbWB7z4NeJxtvdwXQc9+Li85ZM3doI4Pk+kfMeGW5l+d5rTrRPzbbI61fguXNz+ERjGGMc3FlvmJt1orzjur3VoSR7h1KcmslXWzsrQdO0N7/2srY+10b3hzc2FjPw60Sz210f2NkfJ1fyq4VMxtYk6uS5HDv64pQ8P+hLr4qEc1obK4O9mmKykqckYsEg9Yn5Enf+vR1L3reZEa4dSrJ3OM3Z2TxKoKkf6ovz4OOnuXbfsBnvbY4Z+EugsZDpbddfXlfgcnGhyGS2+VxJY3NxBWKu03Yf+ZUQIBFzqoNXyj4k3ObtkqMEqcNiaCfhBiMDHx6/wHsa9m/VlqB2OHsi5iyRO270zFej+7AY/BqJVDMT84VqR8gHHz/N3bce4pvvfyN333qIC7nKym9kbApRSGMw1b6P0yxvWTvBaiZXJlfyqgNWgklPTV4TxuxdR+hLuBzeM8ih0QF2DSSbzgho/Dv70Bee5ujxieo+jfkfMLmjEWAefBscPT7Bh7/8DC9MLuD5SswVYo4wlI4juJyYDfXNjvDhLz/DkcNj3Hds3BKsXU6x4q/Kg2/2/7NxgpU0PNfsohDlYzxfuWxosYq1mUSxneEdJnc0WmEGfgWilgMzNTNKy14wBi9fLtYNeqh4ynMTWY4en+C580uGVxldRqRyWq8yhaAdr9T1UG986yj2HnOCRzFXUNWWEsV2ipG2gtxxI/rgG0uxEM0K3HdsPOgn08II1Cokok6Dv/LJR5lesNh7t+M3ylvWAdcVYg3fKmHRk0+4wo/s6OPwngzD6fiKEsX9I33ky/UJ+kbvvNvlju2EmYzOYB78CpycyVFaZSKubJrILcN6/Z+qqmIqPum4S8wPKphdCRQyvkbJVWEwHfSIv3L34Io9X9r1zrtZ7rjWGbHGpdMxD15E/lBEJkTkqU4dYyPYb3FMYwWEYIi6Q6CkuWJXPxVfA4mkr9UOlAoUKj5T2ULbIZRu987bwZLAm0cnPfg/An4P+GQHj9ERauOFmWRs3WK0Ru8RGXdVECfoGjqbL9OXcCmUg8lPjgiOKOVwAMdC0eNjb7+2bSPdzd55O1gSePPomIFX1WMicrBT798papuGuQInphdMDWPUEVWiQpAslbAltOMHM1C/+f43cvT4BO/45HdwWVTTxFxh71AKX5tPbIoci+fOz1H2lETMqfaU38oGfiskgXuVjvaDDw38F1X1Ne3s3w394KMe2XP5MlPZkhl3YwkCOI6gvhKPOdVh3OWKj+sG8tlSxScbzmkVEVKxoF98zBXijjDSn6xTlAB86AtPU6p4TC+UggMp7MokiLvulgvLNGJDtjtHV/eDF5G7gLsADhw4sMmrCeKFrmDG3agj7gp7h9JkUrHqtKMXp7PM5SuUPB+HoIBJgNlcOTDsDqgPoFQ8P+gNj9KfjFH2tU5R0hd3iLvCdLaCg+A4gq/KXL7CnqHYlk9IbvUw01Zl0w28qt4P3A+BB79Rx408iqdOXyRXDnqHDCRjJFzhQtgx0DAg0MvvHUojEgzSUIKpXBVP2TOUIh13eWEyi2g4zFqCn76v+G5QM+F5Sl/coeIHMfiRPq2bw3piOseVYwPVgSAQXDBKnm8JSWPNbLqB3whaDb+eL9RPxpkrWGsBIyDuCv0Jl/mix1gmASinZgoAXD6cohz2eS9V/LrB14FqBmKOgyNQqSgJV1DgFWMZjp+bw5Fg4HYmFRQwRQqTfNkj4TpUPK2GfRKuYwlJY810Uib5APAwcLWInBKRd3TqWMvRrMji40d/SKniMZs3g24sJekKQ6kYr7psiF974ys5uHOAc3NFYq5w+XCawXSCvkSMmCNMzBcZ6YvXtSTwFcqeXzccO5rwFP0seYu1Ffmyx6Fd/ZQ9JZOK4aNUfB/fVwbTMUtIGmumkyqaOzv13quhWZFFxfeZL1SsR7tRhxCEY67Y1c9Xfv0nqtvfA9xy70MMp+OILJry+UKl+rfkCnVtgZWg10zMCQz+roGg58yugSSnL+aJOfUtCn7rLYeB4O+14s2FbYQdDu4csISksWZ6PkTTrJdHMmwTG03iMYwoTOL5PieajMFrpuUuVnySrkPJ84m5DhJWr0ZTnlThlWMZJrPFas+ZmCsM98UZHUhWk7W1BtwMubGe9LyBb/bFHOqLc2GhjKvKCsPpDQNoruWORufN5StUPCXmOjiqxBxhz1CKsUyKB+66aYlE8Lfe8moz5MaG0PPNxprNtIy7Lm9+zW4z7tuAaEZp4+zSRhTFD73vQ7v6lzzfrGXAu468grjrrhg3P3J4jAfuuqk6J+C+Y+Pccu9D3Hn/I9Zwy+goPe/BN5tpuWcwwee/d3azl2Z0mGjK0lAqxoXc8t09C+Wgl/tgOsb7bz9c3b5Sm9todF47cfPaKulaDfw9WGjG6Aw9a+CbSSMfHofnJ+Z5eLy02cszNoCrdw8w0p/kuy/P4DpBT5ha9YoAqZigIni+4jpC3F28qW3HIK+mgMe6KhobTU8a+MYv5ovTWR4en6529DO2ByLC3bceCnrChE2/kq6Dj7J3KMW5uSJ7R9JLZpne+5Xj3HdsnMdfnkGAPUOpuqKktRrkdoZ3GMZ60pMx+FpPKVuscHa2CJhx325kixWOHB7jytEBHEfwQhXL3qE0sdBTb2xjW/F8npvIMjFfCGPyypmLBeYL5er+azXI7QzvMIz1pCcNfNR/er5Q5tRMfrOXY2wSA8nAM//AT7+KsUyKAzv6uGJXPzFXKHvKoV39Swzu+bli1TlIuA5C0Clycj5wEi7FIDdL+FsRk9FJetLAR57SudlCXTWh0VtEzb2a4QpEnVJbDc14/+2Hlxpc32d3ZrEoyUdRVYoV75INci8M7zC2Fj0Xgz96fIKZhSIvTi/grW7SnrHFaNbpOlLO7BlMslBa9M5bJUMbFVbxmqHZg2G8/Px8AVFhLJO65KpS66pobCQ9Y+CPHp/g3q8c59nz86havL2XcQR29SeYWihVK5FdgVfvHarukytVGMukVnyvRoMbJeijgqaYK+ZpG1uWnjDw0Zfy5QumRugFBPiHr7sMgM+F9QpRnxgFBhIumXScscEU+bLHXD5o77weE4Oa1U1YLxhjq9ITBv6+Y+NMZwubvQzjEkjGHCq+8oaDO+oM6hW7nuMT3zrBQsmjP+HyzluuqBYX1Zb+w/oZZQujGL1CR0f2rZa1juy75d6HTC2zhUm6wr4dfdXeLYZhtM9yI/u2vIrm6PEJ5vLLl6Ebm4cj0J9o/WcmwI6BhMkFDaMDbOkQzdHjE3zw80+t2EjK2FiEwKiXfcgkXa7cPcjdtx7iyVMX+f1vjJMrewhBA7DBdNx6nhtGh9hSBr5U8Xnu/DxPnLrIkydn+a9PniFXspaQG4Ur4GlgwDOpGO+85Qrec9tV1ecb2+I2Gu0jh8fq9jcMo7N01MCLyO3ARwEX+ISqfrjd1/q+Mj6V5XsnZ3ny1EWeODXLM2fnKFWai9tjjlhRUwvcUH3SKB9NxYTXH9jBzYd28KXvn+XEdA7fV+KuQ1/C4crdg2GTtgttJS8tOWkY3UXHkqwi4gLPAT8JnAK+A9ypqj9o9Zorr7lOf/kjn+aJUxd56vQc2WLzmaljmSTX7hvm2XNzVHyf4XSCmOvwg7NzeNvMyAvBxe2VYwMc3jPAl546Xx0AHXfg3W+8kvfcdtWK3rVhGFuT5ZKsnfTg3wC8oKrj4SI+DdwBtDTwL1/Icd+x8bptQ+k41+4bCv8Nc92+YfYMBQUskf695PlBq1dHgtiuA/lKbxp6R+CO6y7jd99+fdPnf7fF68y7NoztRycN/OXAyZrfTwE/2riTiNwF3AWQ3PNK3nBwR2DM9w9z3b4hDuzoqxt0XEtjUcoVu/qZzBYZSsc5M7NAvvkNwKYRd+DGgzur3vPHvrZU420xasMw1otOhmh+DvgpVX1n+Ps/At6gqu9u9ZobbrxRH1uDDr6W2lBEf8JFRHj23DyX2pYmSjACJFxh92BwFxEdI1usWOjDMIwNZ7NCNKeA/TW/7wPOLPeC9RA7WijCMAwjoJOFTt8BrhSRK0QkAbwd+EIHj2cYhmHU0DEPXlUrIvLPgb8mkEn+oao+3anjGYZhGPV0VAevql8CvtTJYxiGYRjN2fK9aAzDMIzmmIE3DMPoUczAG4Zh9Chm4A3DMHqUrhr4ISKTwEsNm3cBU5uwnG7EzsUidi4C7Dwssl3PxY+o6mizJ7rKwDdDRB5tVaW13bBzsYidiwA7D4vYuViKhWgMwzB6FDPwhmEYPcpWMPD3b/YCugg7F4vYuQiw87CInYsGuj4GbxiGYayNreDBG4ZhGGvADLxhGEaP0tUGXkRuF5FnReQFEfnAZq9nvRGR/SLy30TkGRF5WkR+Ldy+Q0S+KiLPhz9Hal7zm+H5eFZEfqpm+w0i8v3wuY9JqzFYXYyIuCLyXRH5Yvj7dj0PwyLyoIgcD/82bt7G5+LXw+/GUyLygIiktuu5WBOq2pX/CFoM/xA4BCSAJ4BXb/a61vkzXgZcHz7OEAwpfzXwfwMfCLd/ALg3fPzq8DwkgSvC8+OGz/0dcDPB3JQvAz+92Z9vDefjXwB/Cnwx/H27noc/Bt4ZPk4Aw9vxXBCM/TwBpMPfPwP80+14Ltb6r5s9+OrQblUtAdHQ7p5BVc+q6uPh43ngGYI/6jsIvuSEP/9h+PgO4NOqWlTVE8ALwBtE5DJgUFUf1uCv+ZM1r9kSiMg+4C3AJ2o2b8fzMAjcCvwBgKqWVPUi2/BchMSAtIjEgD6CqXDb9Vysmm428M2Gdl++SWvpOCJyEHg98G1gt6qeheAiAEQzCFudk8vDx43btxL/AfiXUDc+dzueh0PAJPBfwnDVJ0Skn214LlT1NPDvgJeBs8Csqv4N2/BcrJVuNvDNYmQ9qekUkQHgL4D3qurccrs22abLbN8SiMjPABOq+li7L2mybcufh5AYcD3wn1T19cACQRiiFT17LsLY+h0E4Za9QL+I/NJyL2myrSfOxVrpZgO/6qHdWxERiRMY9z9R1c+Gm8+Ht5WEPyfC7a3OyanwceP2rcKPA28VkRcJQnFvFJFPsf3OAwSf4ZSqfjv8/UECg78dz8VtwAlVnVTVMvBZ4MfYnudiTXSzge/5od1hJv8PgGdU9d/XPPUF4J+Ej/8J8Pma7W8XkaSIXAFcCfxdeJs6LyI3he/5j2te0/Wo6m+q6j5VPUjw//khVf0lttl5AFDVc8BJEbk63PQm4Adsw3NBEJq5SUT6ws/wJoI81XY8F2tjs7O8y/0D3kygLPkh8K83ez0d+Hy3ENwqPgl8L/z3ZmAn8HXg+fDnjprX/OvwfDxLjRIAuBF4Knzu9wirlLfaP+AIiyqabXkegNcBj4Z/F38JjGzjc/F/AMfDz/H/EShktuW5WMs/a1VgGIbRo3RziMYwDMO4BMzAG4Zh9Chm4A3DMHoUM/CGYRg9ihl4wzCMHsUMvGEYRo9iBt7YFoTFL18Tke+JyC+IyL+6hPfKrufaDKNTxDZ7AYaxQbweiKvq66BqpP+vTV2RYXQY8+CNLYuI9IvIX4nIE+FAiF+QYEjMcRH5VjjY4YsiMgZ8Cnhd6MH/OUEL2u+JyJ8s8/5/KSKPhQMn7mp47v8RkcdF5OsiMhpue52IPCIiT4rI50RkREReJSJ/V/O6gyLyZPj4BhH5RniMv476qxjGemEG3tjK3A6cUdXrVPU1wFeA/wz8LPD3gT0AqjoBvBP4pqq+TlV/DsiHj39xmff/ZVW9gaDM/T0isjPc3g88rqrXA98Afjvc/kng/ap6LfB94LdV9RkgISKHwn1+AfhM2GTuPwJvC4/xh8C/ueQzYhg1mIE3tjLfB24TkXtF5O8TtJU9oarPa9CD41OX+P7vEZEngEcIuhReGW73gT8LH38KuEVEhoBhVf1GuP2PCQZ3QDCJ6OfDx78QvvZq4DXAV0Xke8AHqe94aBiXjMXgjS2Lqj4nIjcQNGj7t8DfsE59vkXkCEG72ptVNSciR4FUq6Ws8HZ/Bvy5iHw2WLY+LyKvBZ5W1ZvXY72G0Qzz4I0ti4jsBXKq+imCyT8/BlwhIq8Id7lzmZeXwzBJK4aAmdC4HwZuqnnOAd4WPv6fgW+p6iwwE95JAPwjgvANqvpDwAN+i0XP/1lgVERuDj9LXESuWfFDG8YqMA/e2Mq8FviIiPhAGfhfgF3AX4nIFPAtgjBIM+4HnhSRx1vE4b8C/GqYEH2WIEwTsQBcIyKPAbMEYRcIepP/voj0AePAP6t5zZ8BHyEII6GqJRF5G/CxMLwTIxhb+PQqPr9hLIu1CzZ6ljDM8j5V/ZlNXophbAoWojEMw+hRzIM3tjWh9PHrTZ56k6pOb/R6DGM9MQNvGIbRo1iIxjAMo0cxA28YhtGjmIE3DMPoUczAG4Zh9Cj/HdO4e6wy+Ov4AAAAAElFTkSuQmCC\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, ci = None)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can use the Pandas method <code>corr()</code>  to find the feature other than price that is most correlated with price.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "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": 47,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.corr()['price'].sort_values()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module 4: Model Development\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can Fit a linear regression model using the  longitude feature <code>'long'</code> and  caculate the R^2.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.00046769430149029567"
      ]
     },
     "execution_count": 48,
     "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",
    "\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.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.49285321790379316"
      ]
     },
     "execution_count": 49,
     "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",
    "\n",
    "Fit a linear regression model to predict the <code>'price'</code> using the list of features:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "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.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6576951666037495"
      ]
     },
     "execution_count": 51,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[features]\n",
    "Y= df['price']\n",
    "lm = LinearRegression()\n",
    "lm.fit(X, Y)\n",
    "lm.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": 52,
   "metadata": {},
   "outputs": [],
   "source": [
    "Input=[('scale',StandardScaler()),('polynomial', PolynomialFeatures(include_bias=False)),('model',LinearRegression())]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 8\n",
    "\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.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Pipeline(memory=None,\n",
       "     steps=[('scale', StandardScaler(copy=True, with_mean=True, with_std=True)), ('polynomial', PolynomialFeatures(degree=2, include_bias=False, interaction_only=False)), ('model', LinearRegression(copy_X=True, fit_intercept=True, n_jobs=None,\n",
       "         normalize=False))])"
      ]
     },
     "execution_count": 53,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pipe=Pipeline(Input)\n",
    "pipe"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Pipeline(memory=None,\n",
       "     steps=[('scale', StandardScaler(copy=True, with_mean=True, with_std=True)), ('polynomial', PolynomialFeatures(degree=2, include_bias=False, interaction_only=False)), ('model', LinearRegression(copy_X=True, fit_intercept=True, n_jobs=None,\n",
       "         normalize=False))])"
      ]
     },
     "execution_count": 54,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pipe.fit(X,Y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.751340143833189"
      ]
     },
     "execution_count": 55,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pipe.score(X,Y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module 5: Model Evaluation and Refinement\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Import the necessary modules:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "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:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "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",
    "\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": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.linear_model import Ridge"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6478759163939116"
      ]
     },
     "execution_count": 59,
     "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",
    "\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.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.7002744266408322"
      ]
     },
     "execution_count": 60,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from sklearn.preprocessing import PolynomialFeatures\n",
    "from sklearn.linear_model import Ridge\n",
    "pr = PolynomialFeatures(degree=2)\n",
    "x_train_pr = pr.fit_transform(x_train)\n",
    "x_test_pr = pr.fit_transform(x_test)\n",
    "poly = Ridge(alpha=0.1)\n",
    "poly.fit(x_train_pr, y_train)\n",
    "poly.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://cf-courses-data.s3.us.cloud-object-storage.appdomain.cloud/IBMDeveloperSkillsNetwork-DA0101EN-SkillsNetwork/labs/FinalModule_Coursera/images/share_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://cf-courses-data.s3.us.cloud-object-storage.appdomain.cloud/IBMDeveloperSkillsNetwork-DA0101EN-SkillsNetwork/labs/FinalModule_Coursera/images/url_notebook.png\" alt=\"HTML\" style=\"display: block; margin-left: auto; margin-right: auto;\" /></p>\n",
    "        <p>&nbsp;</p>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<h2>About the Authors:</h2> \n",
    "\n",
    "<a href=\"https://www.linkedin.com/in/joseph-s-50398b136/?utm_medium=Exinfluencer&utm_source=Exinfluencer&utm_content=000026UJ&utm_term=10006555&utm_id=NA-SkillsNetwork-Channel-SkillsNetworkCoursesIBMDeveloperSkillsNetworkDA0101ENSkillsNetwork20235326-2021-01-01\">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.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Other contributors: <a href=\"https://www.linkedin.com/in/michelleccarey/?utm_medium=Exinfluencer&utm_source=Exinfluencer&utm_content=000026UJ&utm_term=10006555&utm_id=NA-SkillsNetwork-Channel-SkillsNetworkCoursesIBMDeveloperSkillsNetworkDA0101ENSkillsNetwork20235326-2021-01-01\">Michelle Carey</a>, <a href=\"https://www.linkedin.com/in/jiahui-mavis-zhou-a4537814a?utm_medium=Exinfluencer&utm_source=Exinfluencer&utm_content=000026UJ&utm_term=10006555&utm_id=NA-SkillsNetwork-Channel-SkillsNetworkCoursesIBMDeveloperSkillsNetworkDA0101ENSkillsNetwork20235326-2021-01-01\">Mavis Zhou</a>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Change Log\n",
    "\n",
    "| Date (YYYY-MM-DD) | Version | Changed By      | Change Description                           |\n",
    "| ----------------- | ------- | --------------- | -------------------------------------------- |\n",
    "| 2020-12-01        | 2.2     | Aije Egwaikhide | Coverted Data describtion from text to table |\n",
    "| 2020-10-06        | 2.1     | Lakshmi Holla   | Changed markdown instruction of Question1    |\n",
    "| 2020-08-27        | 2.0     | Malika Singla   | Added lab to GitLab                          |\n",
    "\n",
    "<hr>\n",
    "\n",
    "## <h3 align=\"center\"> © IBM Corporation 2020. All rights reserved. <h3/>\n",
    "\n",
    "<p>\n"
   ]
  },
  {
   "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.7.12"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}