IBM- Data analysis with Python_ Predict Housing prices.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": 17,
   "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\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Load the csv:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "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": 19,
   "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": 19,
     "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": 20,
   "metadata": {
    "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": 21,
   "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": 21,
     "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": 22,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [],
   "source": [
    "df.drop(['id','Unnamed: 0'],axis=1,inplace=True)\n"
   ]
  },
  {
   "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": 23,
   "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": [
    "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": 24,
   "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": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "mean=df['bathrooms'].mean()\n",
    "df['bathrooms'].replace(np.nan,mean, inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "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\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": 27,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "     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\n"
     ]
    }
   ],
   "source": [
    "sr = df['floors'].value_counts().to_frame()\n",
    "print(sr)"
   ]
  },
  {
   "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": 28,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\sshar\\anaconda3\\lib\\site-packages\\seaborn\\_decorators.py:36: FutureWarning: Pass the following variables as keyword args: x, y. From version 0.12, the only valid positional argument will be `data`, and passing other arguments without an explicit keyword will result in an error or misinterpretation.\n",
      "  warnings.warn(\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='waterfront', ylabel='price'>"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAERCAYAAABxZrw0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAVG0lEQVR4nO3dfZDdVX3H8c9nd0ESHsayrAwkkYhBEEGe1ocWqzwszmorznRUZKqslpq21iS0RUetU8cZ63SmndokVcuOpWxGxQIFBx1YzRaQBwHdhEAIwXGLAbMoXFfkKTGwu9/+ce+GXbLZXELO7/fLue/XzE7u7/7uveebzd1Pzp57fuc4IgQAyE9b2QUAANIg4AEgUwQ8AGSKgAeATBHwAJApAh4AMlW5gLd9ue3Hbd/f5OM/YPsB25tsfyt1fQCwv3DV5sHbfrukZyStiYiT9vDY4yRdJemciHjC9qsi4vEi6gSAqqtcDz4ibpX0m+n32X6t7UHb62zfZvuExqmPSfpKRDzReC7hDgANlQv43eiXtCwizpB0qaSvNu5/naTX2b7D9l22e0urEAAqpqPsAvbE9iGS/kDS1ban7n5F488OScdJOkvSQkm32T4pIn5bcJkAUDmVD3jVf8v4bUScOsu5rZLuiojnJf3c9k9VD/yfFFgfAFRS5YdoIuIp1cP7/ZLkulMap78j6ezG/UeoPmTzUBl1AkDVVC7gbV8p6U5Jx9veavtiSX8q6WLb90raJOm9jYd/X9KY7Qck3SzpkxExVkbdAFA1lZsmCQDYN5L24G3/TeMCpPttX2n7oJTtAQBekKwHb3uBpNslnRgR221fJemGiLhid8854ogjYvHixUnqAYAcrVu37tcR0TXbudSzaDokzbP9vKT5kh6d68GLFy/W8PBw4pIAIB+2H97duWRDNBExKulfJD0i6ZeSnoyIH8xS3FLbw7aHa7VaqnIAoOUkC3jbv6f6bJfXSDpa0sG2P/Tix0VEf0R0R0R3V9esv2UAAPZCyg9ZeyT9PCJqjQuRrlX9ilQAQAFSBvwjkt5qe77rawycK2lzwvYAANOkHIO/W9I1ktZL2thoqz9Ve3jB2NiYli9frrExrvkCWlnSefAR8fmIOCEiToqID0fEjpTtoW5gYEAbN27UmjVryi4FQIkqt1QBXp6xsTENDg4qIjQ4OEgvHmhhBHxmBgYGNDk5KUmamJigFw+0MAI+M0NDQxofH5ckjY+Pa+3atSVXBKAsBHxmenp61NFRv0C5o6ND5513XskVASgLAZ+Zvr4+tbXV/1nb29t10UUXlVwRgLIQ8Jnp7OxUb2+vbKu3t1ednZ1llwSgJPvDln14ifr6+rRlyxZ670CLI+Az1NnZqVWrVpVdBoCSMUQDAJki4AEgUwQ8AGSKgAeATBHwAJApAh4AMkXAA0CmCHgAhWJDmuIQ8AAKxYY0xUkW8LaPt71h2tdTti9J1R6A6mNDmmKl3JP1pxFxakScKukMSdskXZeqPQDVx4Y0xSpqiOZcSf8XEQ8X1F5LY4wTVcWGNMUqKuA/KOnK2U7YXmp72PZwrVYrqJy8McaJqmJDmmIlD3jbB0o6X9LVs52PiP6I6I6I7q6urtTlZI8xTlQZG9IUq4ge/LskrY+Ixwpoq+UxxokqY0OaYhUR8BdqN8Mz2PcY40TV9fX16eSTT6b3XoCkAW97vqTzJF2bsh28gDFOVN3UhjT03tNLGvARsS0iOiPiyZTt4AWMcQKYwpWsmWGME8AU9mTNEJtuA5AI+Cyx6TYAiSEaAMgWAQ8AmSLgASBTBDwAZIqAB4BMEfAAkCkCHgAyRcADQKYIeADIFAEPAJki4AEgUwQ8AGSKgAeATBHwAJCp1Fv2vdL2NbYftL3Z9u+nbA8A8ILU68GvlDQYEe+zfaCk+YnbAwA0JAt424dJerukj0hSRDwn6blU7QEAZko5RHOspJqk/7J9j+2v2z74xQ+yvdT2sO3hWq2WsBwAaC0pA75D0umSvhYRp0l6VtKnX/ygiOiPiO6I6O7q6kpYDgC0lpQBv1XS1oi4u3F8jeqBDwAoQLKAj4hfSfqF7eMbd50r6YFU7QEAZko9i2aZpG82ZtA8JOmjidsDADQkDfiI2CCpO2UbAIDZcSUrAGSKgM/Q2NiYli9frrGxsbJLAVAiAj5DAwMD2rhxo9asWVN2KQBKRMBnZmxsTIODg4oIDQ4O0osHWhgBn5mBgQFNTk5KkiYmJujFAy2MgM/M0NCQxsfHJUnj4+Nau3ZtyRUBKAsBn5menh51dNRnv3Z0dOi8884ruSIAZSHgM9PX16e2tvo/a1tbmy666KKSKwJQFgI+M52dnTr66KMlSUcffbQ6OztLrgiYiWm8xSHgMzM2NqbR0VFJ0qOPPsoPESqHabzFIeAzMzAwoIiQJE1OTvJDhEphGm+xCPjMMIsGVcY03mIR8JlhFg2qjA5IsQj4zEyfRdPe3s4sGlQKHZBiEfCZ6ezsVG9vr2yrt7eXWTSolL6+vp1DNJOTk3RAEku94QdK0NfXpy1btvDDA7Q4evAZ6uzs1KpVq+i9o3IGBgZkW5Jkmw9ZE0sa8La32N5oe4Pt4ZRtAai+oaEhTUxMSKrPouFD1rSK6MGfHRGnRgRb9wEtjg9Zi8UQDYDCMMurWKkDPiT9wPY620tne4DtpbaHbQ/XarXE5QAoE7O8ipV6Fs2ZEfGo7VdJWmv7wYi4dfoDIqJfUr8kdXd3R+J6AJSMWV7FSdqDj4hHG38+Luk6SW9O2R6A6mOWV3GSBbztg20fOnVb0jsl3Z+qPQDATCmHaI6UdF1jzmuHpG9FxGDC9gAA0yTrwUfEQxFxSuPrDRHxj6naArD/YMOP4jBNEkCh2PCjOAQ8gMKw4UexCHgAhWHDj2IR8BlijBNVxYYfxSLgM8QYJ6qqp6dnxmqSrEWTFgGfGcY4UWXnn3/+zk3hI0Lvec97Sq4obwR8ZhjjRJVdf/31M3rw3/3ud0uuKG8EfGYY40SVDQ0NzejB8/5Mi4DPDOtto8p4fxaLgM8M622jynh/FouAzwzrbaPKeH8Wi4DP0Pnnn6/58+czQwGV1NfXp5NPPpneewEI+Axdf/312rZtGzMUUEmsB1+cpgPe9jG2exq3502t9Y5qYR48gClNBbztj0m6RtJljbsWSvpOoprwMjAPHsCUZnvwfy3pTElPSVJE/EzSq1IVhb3HPHgAU5oN+B0R8dzUge0OSWyQXUHMMwYwpdmA/6Htz0qaZ/s8SVdLauoTPNvttu+x/b29LRLN6+vr2zlEMzk5yUwFoIU1G/CfllSTtFHSX0i6QdLnmnzuCkmbX3ppAICXo9mAnyfp8oh4f0S8T9LljfvmZHuhpD+S9PW9LxEvxcDAwIzFnPiQFWhdzQb8/2pmoM+TNNTE8/5N0qckTe7uAbaX2h62PVyr1ZosB7szNDSkiYkJSfVZNHzICrSuZgP+oIh4ZuqgcXv+XE+w/ceSHo+IdXM9LiL6I6I7Irq7urqaLAe709PTs3Otj7a2Nj5kBVpYswH/rO3Tpw5snyFp+x6ec6ak821vkfRtSefY/sZeVYmm8SErgCkdTT7uEklX2360cXyUpAvmekJEfEbSZyTJ9lmSLo2ID+1VlWjaE088scsxl4RDklavXq2RkZGyy9Do6KgkacGCBaXWsWTJEi1btqzUGlJrqgcfET+RdIKkv5L0cUmv39PQC8rxxS9+cc5joGzbt2/X9u17GgDAvjBnD972ORFxk+0/edGp42wrIq5tppGIuEXSLXtXIl6KLVu2zHmM1lWV3uqKFSskSStXriy5kvztaYjmHZJukjTburMhqamAR3EWLlyorVu37jxetGhRidUAKNOcAR8Rn7fdJunGiLiqoJrwMixatGhGwC9cuLDEagCUaY9j8BExKekTBdSCfeDuu++e8xhA62h2muRa25faXmT78KmvpJVhr0ztWL+7YwCto9lpkn+m+pj7x190/7H7thy8XG1tbTuvZJ06BtCamv3pP1HSVyTdK2mDpNWS3pCoJrwMPT09cx4DaB3NBvyApNdLWqV6uL++cR8qZunSpXMeA2gdzQ7RHB8Rp0w7vtn2vSkKAgDsG8324O+x/dapA9tvkXRHmpLwclx22WUzjvv7+0uqBEDZmg34t0j6ke0tjcXD7pT0Dtsbbd+XrDq8ZENDM1dxZrlgoHU1O0TTm7QK7DNTK0nu7hhA62gq4CPi4dSFAAD2LSZJA0CmCPjMHH744XMeA2gdBHxmnnzyyTmPAbQOAj4z05cpmO0YQOsg4AEgU8kC3vZBtn9s+17bm2x/IVVbAIBdNTsPfm/skHRORDxj+wBJt9u+MSLuStgmAKAhWcBHfSHyZxqHBzS+WJwcAAqSdAzedrvtDZIel7Q2InbZXsj2UtvDtodrtVrKcgCgpSQN+IiYiIhTJS2U9GbbJ83ymP6I6I6I7q6urpTlAEBLKWQWTUT8VtItYk0bAChMylk0XbZf2bg9T1KPpAdTtQcAmCnlLJqjJA3Yblf9P5KrIuJ7CdsDAEyTchbNfZJOS/X6AIC5cSUrAGSKgAeATBHwAJApAh4AMkXAA0CmCHgAyBQBDwCZIuABIFMEPABkioAHgEwR8ACQKQIeADJFwANAplIuFwygYfXq1RoZGSm7jEqY+j6sWLGi5EqqYcmSJVq2bFmS1ybggQKMjIzoZ5vu0asPmSi7lNId+Hx94GDHw8MlV1K+R55pT/r6BDxQkFcfMqHPnv5U2WWgQr60/rCkr88YPABkKuWerIts32x7s+1NthlwA4ACpRyiGZf0dxGx3vahktbZXhsRDyRsEwDQkKwHHxG/jIj1jdtPS9osaUGq9gAAMxUyBm97seobcN89y7mltodtD9dqtSLKAYCWkDzgbR8i6X8kXRIRu0whiIj+iOiOiO6urq7U5QBAy0ga8LYPUD3cvxkR16ZsCwAwU8pZNJb0n5I2R8S/pmoHADC7lD34MyV9WNI5tjc0vt6dsD0AwDTJpklGxO2SnOr1AQBz40pWAMgUAQ8AmSLgASBTBDwAZIqAB4BMEfAAkCkCHgAyRcADQKYIeADIFHuyAgUYHR3Vs0+3J9+DE/uXh59u18Gjo8lenx48AGSKHjxQgAULFmjH+C/12dN32RIBLexL6w/TKxak2+iOHjwAZIqAB4BMEfAAkCkCHgAyxYes+9jq1as1MjJSdhkzrFixorS2lyxZomXLlpXWPtDKUu7Jerntx23fn6oNAMDupezBXyHp3yWtSdhG5ZTdWz3rrLN2uW/lypXFFwKgdMl68BFxq6TfpHp9zO6ggw6acTxv3rySKgFQttI/ZLW91Paw7eFarVZ2Ofu9wcHBGcc33nhjSZUAKFvpAR8R/RHRHRHdXV1dZZeTFXrvQGtjFk2GTjnlFEmMvQOtrvQePAAgjZTTJK+UdKek421vtX1xqrYAALtKNkQTERemem1gf/TIM6wHL0mPbav3K4+cP1lyJeV75Jl2HZfw9RmDBwqwZMmSskuojOcaV3q/4hi+J8cp7XuDgAcKUPYFcFUytXQGkwDSyybgq7gGTFmmvg9lrkFTJayHg1aVTcCPjIxow/2bNTH/8LJLKV3bcyFJWvfQYyVXUr72bVxMjdaVTcBL0sT8w7X9hHeXXQYqZN6DN5RdAlAa5sEDQKYIeADIVDZDNKOjo2rf9iS/kmOG9m1jGh0dL7sMoBT04AEgU9n04BcsWKBf7ejgQ1bMMO/BG7RgwZFllwGUgh48AGQqmx68VJ/zzBi81Pa7pyRJkwex7kl9Hjw9eKk6FwNW5UK8VrgALpuAZ62PF4yMPC1JWnIswSYdyXujYtiIpjiOiLJr2Km7uzuGh4fLLmO/x1ofQOuwvS4iumc7xxg8AGSKgAeATBHwAJCppB+y2u6VtFJSu6SvR8Q/pWyvCqowU6EqsxSk1pipAFRVyj1Z2yV9RdK7JJ0o6ULbJ6ZqDy+YN28eMxUAJO3Bv1nSSEQ8JEm2vy3pvZIeSNhm6eitAqiKlGPwCyT9Ytrx1sZ9M9heanvY9nCtVktYDgC0lpQB71nu22XSfUT0R0R3RHR3dXUlLAcAWkvKgN8qadG044WSHk3YHgBgmpQB/xNJx9l+je0DJX1Q0vUJ2wMATJPsQ9aIGLf9CUnfV32a5OURsSlVewCAmZLOg4+IGySxvCMAlIArWQEgUwQ8AGSqUssF265JerjsOjJxhKRfl10EsBu8P/edYyJi1jnmlQp47Du2h3e3RjRQNt6fxWCIBgAyRcADQKYI+Hz1l10AMAfenwVgDB4AMkUPHgAyRcADQKYI+AzZ7rX9U9sjtj9ddj3AFNuX237c9v1l19IKCPjMsFUiKu4KSb1lF9EqCPj87NwqMSKekzS1VSJQuoi4VdJvyq6jVRDw+Wlqq0QA+SPg89PUVokA8kfA54etEgFIIuBzxFaJACQR8NmJiHFJU1slbpZ0FVsloipsXynpTknH295q++Kya8oZSxUAQKbowQNApgh4AMgUAQ8AmSLgASBTBDwAZIqAR0uyfYnt+XvxvBNsb7B9j+3XllUH0AymSaIl2d4iqTsifv0SntMu6ZOS5kXE5190zqr/PE2mrgNoFj147Ndsf8r28sbtL9u+qXH7XNvfsP0128O2N9n+QuPccklHS7rZ9s2N+95p+07b621fbfuQxv1bbP+D7dslXSDpEkl/bvtm24ttb7b9VUnrJS2y/c+277e90fYFjdc4y/Yttq+x/aDtb7pulzqAfYmAx/7uVkl/2LjdLekQ2wdIepuk2yT9fUR0S3qjpHfYfmNErFJ9fZ6zI+Js20dI+pyknog4XdKwpL+d1sbvIuJtEfEtSf8h6csRcXbj3PGS1kTEaY32T5V0iqQeSf9s+6jG405T/T+HEyUdK+nMF9exT78rgAh47P/WSTrD9qGSdqh+GXy36qF/m6QP2F4v6R5Jb1A9YF/srY3777C9QVKfpGOmnf/vOdp/OCLuatx+m6QrI2IiIh6T9ENJb2qc+3FEbG0M4WyQtPgl/j2Bl6yj7AKAlyMinm+MY39U0o8k3SfpbEmvlbRd0qWS3hQRT9i+QtJBs7yMJa2NiAt308yzc5Qw/dxsSzVP2THt9oT42UMB6MEjB7eqHuS3qt5r/0vVe8mHqR7AT9o+UvVtDKc8LenQxu27JJ1pe4kk2Z5v+3V7WccFttttd0l6u6Qf7+E50+sA9ikCHjm4TdJRku5sDI38TtJtEXGv6kMzmyRdLumOac/pl3Sj7ZsjoibpI5KutH2f6oF/wl7UcZ3qv0HcK+kmSZ+KiF/t4Tk769iL9oA5MU0SADJFDx4AMkXAA0CmCHgAyBQBDwCZIuABIFMEPABkioAHgEz9P6beAVZWSuPuAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.boxplot(df['waterfront'],df['price'])"
   ]
  },
  {
   "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": 29,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 8081250.0)"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAESCAYAAAD38s6aAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABGIElEQVR4nO29fXxc53Xf+T33ziswAAgSAEmRokhKlGjJlm1J69iJKiu2s2snqb3bj5pETbpt8yI1m9hxuk7jdBPns8p2N9pkk1hdbyrGSZvUqV1bdRJtGjuJI9Oyu5ITSbZsyaJEiZREiiIBgniZwbzee8/+ce8dzAAzwIDEAIPB+eoDYXDnzr3PXM787nnOOc85oqoYhmEY/Yez2QMwDMMwuoMJvGEYRp9iAm8YhtGnmMAbhmH0KSbwhmEYfYoJvGEYRp/ScwIvIn8gIpMi8kyH+/+QiHxHRJ4Vkf/Y7fEZhmFsFaTX8uBF5A6gAPyRqr5xlX2PAJ8F3qWqMyIyoaqTGzFOwzCMXqfnLHhVfRS41LhNRK4VkS+KyJMi8lURORo99VPAJ1R1JnqtibthGEZEzwl8G44BH1TVW4GPAP9PtP164HoR+a8i8riIvHfTRmgYhtFjJDZ7AKshIjngu4HPiUi8OR39TgBHgDuB/cBXReSNqjq7wcM0DMPoOXpe4AlnGbOq+pYWz50FHlfVGnBaRJ4nFPy/28DxGYZh9CQ976JR1XlC8f6HABLy5ujpPwW+N9o+RuiyObUZ4zQMw+g1ek7gReTTwGPADSJyVkR+AvhR4CdE5GngWeAD0e5/CUyLyHeALwO/oKrTmzFuwzCMXqPn0iQNwzCM9aGrFryI/Hy0AOkZEfm0iGS6eT7DMAxjka5Z8CKyD/gacKOqlkTks8BfqOq/b/easbExPXjwYFfGYxiG0Y88+eSTF1V1vNVz3c6iSQBZEakBA8C5lXY+ePAgTzzxRJeHZBiG0T+IyCvtnuuai0ZVXwN+E3gVeB2YU9W/6tb5DMMwjGa6JvAiMkqY7XIIuAoYFJEfa7HfPSLyhIg8MTU11a3hGIZhbDu6GWR9D3BaVaeihUifJ1yR2oSqHlPV21T1tvHxlm4kwzAM4zLopsC/CrxdRAYkrDHwbuC5Lp7PMAzDaKCbPvivAw8BTwHfjs51rFvnMwzDMJrpahaNqv4q8KvdPIdhGIbRmq1QbMwwDKNnOX5ikgcfPcWZmSJXjw5w7x2HufPoxGYPC+jBWjSGYRhbheMnJvnYw88ymS+zI5tkMl/mYw8/y/ETvdF7yATeMAzjMnnw0VMkXWEglUAk/J10hQcf7Y2itibwhmEYl8mZmSLZpNu0LZt0OTtT3KQRNWMCbxiGcZlcPTpAqeY3bSvVfPaPDmzSiJoxgTcMw7hM7r3jMDVfKVY9VMPfNV+5947Dmz00wATeMAzjsrnz6AT3vf8mJoYyzJVqTAxluO/9N/VMFo2lSRqGYVwBdx6d6BlBX4pZ8IZhGH2KCbxhGEafYgJvGIbRp5jAG4Zh9Ckm8IZhGH2KCbxhGEafYgJvGIbRp5jAG4Zh9Ckm8IZhGH2KrWQ1DGPL0svNNnqBrlnwInKDiHyz4WdeRD7crfMZhrG96PVmG71AN5tuP6+qb1HVtwC3AkXgT7p1PsMwthe93myjF9goH/y7gZdU9ZUNOp9hGH1Orzfb6AU2SuB/BPj0Bp3LMIxtQK832+gFui7wIpIC3g98rs3z94jIEyLyxNTUVLeHYxhGn9DrzTaWcvzEJHcfe5zb73+Eu489viGxgo2w4N8HPKWqF1o9qarHVPU2Vb1tfHx8A4ZjGEY/0OvNNhrZrIDwRqRJ3o25ZwzD6AK93GyjkcaAMMBAKkGx6vHgo6e6Ov6uWvAiMgB8H/D5bp7HMAyjl9msgHBXBV5Vi6q6S1XnunkewzCMXmazAsJWqsAwDKPLbFZA2ATeMAyjy2xWQNhq0RiGYWwAmxEQNgveMAyjTzGBNwzD6FNM4A3DMPoUE3jDMIw+xQTeMAyjTzGBNwzD6FNM4A3DMPoUE3jDMIw+xQTeMAyjTzGBNwzD6FNM4A3DMPoUE3jDMIw+xQTeMAyjTzGBNwzD6FNM4A3DMPqUbvdk3SEiD4nICRF5TkTe0c3zGYZhGIt0u+HHx4EvqupdIpICutuA0DAMw6jTNYEXkWHgDuCfAqhqFah263yGYRhGM9100RwGpoB/JyLfEJFPisjg0p1E5B4ReUJEnpiamuricAzDMLYX3RT4BHAL8Luq+lZgAfjo0p1U9Ziq3qaqt42Pj3dxOIZhGNuLbgr8WeCsqn49+vshQsE3DMMwNoCuCbyqngfOiMgN0aZ3A9/p1vkMwzCMZrqdRfNB4I+jDJpTwD/r8vkMwzCMiK4KvKp+E7itm+cwDMMwWmMrWQ3DMPoUE3jDMIw+xQTeMAyjT+l2kNUwjB7g+IlJHnz0FGdmilw9OsC9dxzmzqMTmz0so8uYBW8Yfc7xE5N87OFnmcyX2ZFNMpkv87GHn+X4icnNHprRZUzgDaPPefDRUyRdYSCVQCT8nXSFBx89tdlDM7qMCbxh9DlnZopkk27TtmzS5exMcZNGZGwUJvCG0edcPTpAqeY3bSvVfPaPWvXufscE3jD6nHvvOEzNV4pVD9Xwd81X7r3j8GYPzegylkVjGH3OnUcnuI/QF392psj+K8yisYycrYMJvGGska0ocHcenViXMcYZOUlXmjJy7ovOYfQW5qIxjDWw3VMOLSNna2ECbxhrYLsLnGXkbC1M4A1jDWx3gbOMnK2FCbxhrIHtLnCWkbO1MIE3jDXQqwJ3/MQkdx97nNvvf4S7jz3etZjAnUcnuO/9NzExlGGuVGNiKMN977/JAqw9iqjqZo+hzm233aZPPPHEZg/DMFYkzqJZj5TD9RpPnNmSTbqUaj41X014twki8qSqtmysZGmShrFG1ivlcL1oDPwCDKQSFKseDz56qqfGaWw8XRV4EXkZyAM+4LW7yxiGcfmcmSmyI5ts2radAr9GezbCgv9eVb24AecxjG3J1aMDTObLdQsetlfg12iPBVkNY4vTq4FfY/PptsAr8Fci8qSI3NNqBxG5R0SeEJEnpqamujwcw+g/LLPFaEdXs2hE5CpVPSciE8BfAx9U1Ufb7W9ZNIZhGGtjpSyarlrwqnou+j0J/Anwtm6ezzAMw1ikawIvIoMiMhQ/Bv5b4Jlunc8wDMNopptZNLuBPxGR+Dz/UVW/2MXzGYZhGA10TeBV9RTw5m4d3zAMw1gZW8lqGG3Yio09DKMRy4M3jBZs98YeRn9gAm8YLdjujT2M/sAE3jBasN0bexj9gQm8YbRguzf2MPoDC7IaRgvuveMwH3v4WYpVr6nGej/Vd2kMIudSLiJCvuJZQLmPMAveMFrQ7/VdGoPIrsCLUwucnCzgChZQ7iPMgjeMNvRaY4/1pDGIfGqqgCsCAhcLVQ6P56xhSJ9gFrxhbEMag8hVP0AERMLHYAHlfsEE3jC2IY1B5JTroAqq4WOwgHK/0LHAi8g1IvKe6HE2LiRmGMbWo7FJyFguha+KHyhjuZQ1DOkjOhJ4Efkp4CHgwWjTfuBPuzQmwzC6TGMQOVC4bnyQIxM5AqXvAsrbmU6DrD9DWMv96wCqejJq4mEYxhaln4PIRkinAl9R1WpU+hcRSRC24zMMYx2wwmZGN+jUB/8VEflXQFZEvg/4HPD/dm9YhrF9sMJmRrfoVOA/CkwB3wbuBf4C+OVuDcowthNW2MzoFp26aLLAH6jq7wGIiBtts0RZw7hCzswU2ZFNNm2zPHRjPejUgv8bQkGPyQJfWv/hGMb2wwqbGd2iU4HPqGoh/iN63NGnT0RcEfmGiPz55QzQMPqdxpx0VV2Wh378xCR3H3uc2+9/hLuPPW6+eaNjOhX4BRG5Jf5DRG4FSh2+9ueA59Y6MMPYLqxU2MwCsMaV0KkP/sPA50TkXPT3XuCHV3uRiOwHfgD418C/uJwBGsZ2oF1OemMAFmAglbBCYEbHdCTwqvp3InIUuAEQ4ISq1jp46e8A/xJoW9ZARO4B7gE4cOBAJ8MxjG2DBWCNK2FFgReRd6nqIyLyD5Y8dUREUNXPr/DaHwQmVfVJEbmz3X6qegw4BnDbbbfZ4inDaODq0QEm82X8QJnKV6j6Aa4Ih8YGN3toxhZgNR/8O6Pff7/Fzw+u8trvAd4vIi8DnwHeJSKfuvyhGsb24947DjNfqnF2pkTNDxDAC5SpQsX88MaqrGjBq+qviogDfEFVP7uWA6vqLwG/BBBZ8B9R1R+7zHEaRsd0Y9n/ZpUSuPPoBLsGU+QrHn6gpFyHsVyaRLQQyvzwxkqsmkWjqgHwsxswFsO4YrqRdbLZmSyFqs914zmO7hnm8HiO4WzS/PBGR3SaJvnXIvIREblaRHbGP52eRFWPq+pqLh1ji9JLedrdWPa/2aUEbCGUcbl0KvA/DvxPwFeAJxp+jG3OZlu3S2lsRRdzpdZuN465FlZbCGUY7eg0D/5GQoG/nbBM8FeBf9utQRlbh27kaTf6u3MpFxEhX/E68n3HWSfxeODyrd14HFP5ChfzFfaMZBjKJK/omK2Ov5pf/86jE9xHeK3PzhTZb+WEjQ7pVOD/EJgHHoj+vjva9kPdGJSxdWjM086Xa0zlK1Q8n7MzJY6fmFyzCMUzgqQruAIvTi0AsG9Hpj47uA/aHvfeOw7zsYefpVj1yCZdSjX/sqzdxnHsGU7z2myZszMl9u1QEq5zxRZ04/EbZz7t3ps15zAuh04F/gZVfXPD318Wkae7MSBja9GYp31utowIOCKIwMcefpa7zs7y2KlLHWefNM4ITk0VcEVA4GKhyuHx3Kqzg/WydpfOTEC4kC9zfr7CLQdGr9iCthWqxkbQqcB/Q0TerqqPA4jIdwH/tXvDMrYKscU8OV8GFFRQYPdQhorn84njL7F/NNuRlQrNM4J4UQ8SPobOfN/rYe0uXUE6nE0ylEkwV6rx6XvefkXHbnV8sBWqxvrTaZD1u4D/T0RejhYuPQa8U0S+LSLf6trojJ4nLpSlhMGZhCtcNZJlOJskX/bwgmBN2SeNGSMp10EVVMPHsHHZI93OXLHMGGMj6FTg3wscIlzZ+s7o8fcTrmb9+90ZmrFVuPPoBLccGOXAzoF6njZAxQtIu80fsdWs1MaMkbFcCl8VP1DGcqkNzR7pduaKZcYYG0GnxcZe6fZAjK3F0gyQdxzeyUNPvdYU3HQdYWSg2Q3RaKW2yyJp9KFfNz6IiFCoeEwMZTZ0BWm3M1cGkg6np8Ob3eGxQX7lB46a/91YV0S1d+p73XbbbfrEE5ZevxmsZSl+YwZIY6bKXbfs47FTl+qCGIv+0v3ue/9NAC2PEddB72faXb/t8N6N9UdEnlTV21o912mQ1ehjWqXsfeShpxnPpVvmn7fLAHns1KVlAcib9+9oaQXffezxbZtFYhk0xkZhAr/FWY8iWEsFx/OV2WKNQtnjuoncsuyXtWSAtMto2cgsks0qFNaObr/3Xnu/vcR2uzadBlmNHmS9ygQsXYp/sVDBEfBVW2a/rEcGSLezSOL6OLf+2l9x76ee5OXpQk+UUoDuvvdeKx3RS2zHa2MCv4VZryJYSwUnzjlPNWTANFqY65EB0s0sksYvcrkWEKgyXahRqHgbXiisFd1875tdGK2X2Y7Xxlw0W5i1TvXbTU+XLu93HcHzlfGhdP21jRZmJxkmq02Fu5ml0vhFjhdLKTCVrzCUWb3Ubjz2Fy7MU/OVVMLhyMTQmsa30vvv5nu3BVTt2Y7XxgS+h1irf3AthbVWq33SKDgHdw4wvVDFdQRVbVnPZaXVop3WWelWfZXGL3LKdfB8RZzFmclK7pB47FXPZ77shftXfV6eLqy6CnfpMVZ6/9167+tZbK3f2I7XxgS+R1hr8SlYW2Gt1TI3lgpOfLO5HAtzPbJEriQY1vhFHsulOTdXgiD0R568kKcWBCQdaVkMLR77dMHDQXAcIVBlvuSxZyTR8j0sHevMQmXZ+5/Kl/nQZ77BcDbZ1eDeehVb60e247Uxge8ROhHFVqJ33/tvWibEAHcfe7xpv7VOT6/EwozPFVeXrPqhoM6Vah29/nJudo00fpGHMgl2eSmmF6rUAiWVEPbvyFILtOUx47HX6+AAEtXCaXW9Wo315ekF9u/I1veZL9WYXqgSqHJg58Ca389asNLC7dmO18YEvkdYTYDbit77b2rKPW+331A6Qanmb8j09OrRAV6eLjBdqCECriPUAiVf9joqIXylM4ClX+RDYzlGB6pU/aDp/bc6Zmz91107slgLp9X1ajXWpONwIV9hOJsCwqwkgEzCrQf3upn3bqWF27Pdrk3XBF5EMsCjQDo6z0Oq+qvdOt9WZzX/YKei124/Va1nbnRzenr8xCRnZ4q8Plepb0s4Yfng0YHkiqL2wJde4JNfO8182cMRGBtMsXsktITbWc8rBTIbz3P7/Y90NIOJrf+hTCK0ugMFheHBJDVfecfhnU2zoxcuzLN3JNt0jN3Dac7OluvXuuIFCDQFrfs9uGf0Bt1Mk6wA74rqyL8FeK+IXHmd1T5ltdS5TtvGtdtvoepz3/tvYmIow1ypxsRQZt2Xxh8/MckvPPQ0r8+Vm7Z7gTKaTTKWS7cVtQe+9AIff+RFSjUfAQKFyUKVC3MlYPlsY605zZ3mnsfVMQ+N5RjJJMgmXUYGkhzcleOuW/bx0FOvNZ2zUPHrFnpMwnW4fiJXv9YDKZexoVS9G1S7c28UvdRD1+guXbPgNSxyU4j+TEY/vVP4psdYzT/YaQbASvt1Oj1dLcDZ7vkHHz1FvuzhRoFJDY1fRGCh6q8oap/82mkcgYTjEAQ+fvRJmSxUcRwhlXCbZhsrzWji5xvH1yrANl+qkXSE2+9/pOl9NF6nxvf67NfmGEy7jGQz9XPuHExyaaHGYDrRNDNqLBwW34x6Ibh3pfENY2vR1WJjIuICTwLXAZ9Q1V9ssc89wD0ABw4cuPWVV6xwZSs6LVB1pYWsGl/v+QEX5ivUgoAj4zk++r43AO2LhP3ynz3D+bkSCdchCKDmB/U7uuvAvh0DLcf74KOneOzUNALRCtrmMaUTDj9z57V86D3X17fd9r/9NcWKFwZOXYfxoTS5dILz8+X6Apal4//+N+2tF0PLpRNMFSqMZJMtC6bF/WCnF6oMR/ucOJ9HgH2joUsmbk8oCNdN5ChUPAbb9JC9kqyky6XVjfjBR08tMwCK1bBS53o0MjE2npWKjW1INUkR2QH8CfBBVX2m3X5WTXJlOhWJte7XSgAaW/CpKo4jTAxlGEy5LYOVE0OhVfuNV2dQwrZ9fqDUotzzoUyCB37krW1vRqcvLhA0fBRFAA1/HxobbBKg4ycmufdTTxKoRrn6YSB0Vy7JQsVnfCjddvzxDebuY48vE7qLhTKXFmrsH82STbq8OFXA85X9o1mGMklOTRXq2TVBNLbGY8cunF6oEtnuRr9QqbF3JItEGUIQvoe5Uo2v/uK7NnSMxvqw6dUkVXVWRI4TNg5pK/DGynTqYulkv3ZT9WLVY89whtMXF+r9VRXwAyXpCqcuLnBkItd0rDgW8GsfeCO/8NDTzBRrqIRq7TrCjoEkv3nXm9vmnA+kEowNppgsVBefjNw7E7n0sljDg4+eYnQgyfRCFQ1CoQ1QLi3UyKVdskm37fjjIG+rrKW5Yg0/0Lro+4HiyOIK2DinvuwFpFxpak+YcIVPfu0040PpnqgS2c6FVfN1w7KpjM2na0FWERmPLHdEJAu8BzjRrfMZa6NdXY6qF1Cq+VT9gNjIi9ME4+Btu2DlnUcn+I273syRiRwigohw3fggv3lX2K99aWCvMSC8eyTLRC5VP6YI7B5KMzGcWSZAZ2aKjOXSXDWSJeEKvoZumqG0y/W7h1ccf3yjaBV0rfgB6cTiVyKuxROvgB3OJtk1mAotd5rbE8aB7E4C4RtBu2B7KuFYJ6ltRDct+L3AH0Z+eAf4rKr+eRfPZ6yBMzNFXKHudki5DmO5FClXqPkauiECDTNaUMZyodAe2jVAsRa0DRi2mj20my3kUm6TNbl7JIvjCDPFRTdJKwGKA8nD2WS9PWDsJoqDqQKUa0H9NYHCifN5BlJu6OJpEXRNOA5DmcWvxFguzWuzJRINJRtSCZeju4eWualKNZ/BJe8n3r4Z1nG7YHtcU2c7LfbZzlhHp23Ke3/7K7w4tYArUl/M46ty3fggH33fG7j/iyd4YbJA0hV2D6WpeAEzxRpDmQSZhMOlhRoVP2Aw5fKTtx9qCoAupZW/u1j1qNR8pgpV/EBJJxySrlCqhe4PESGVcBgbTC0LWsLK3aAe+NIL/Jsvv0htSbTWFZgYTpN03XpXqUaha9WBaq5UYzyXplDxmlYKt+todaU++PWqV77du0Ztp7rvmx5k7RQT+I3jfb/zKCcnC7hOg8AHypGJHF/48B3A4pfk5IV58hWfnYNJUq7Da7Nhnvu+HRkSrrOqcMSLjBoDe/OlKmdnS0wMpZkr1qh4Ab7C6ECCfTsG6uIqUHeBLG3516pEw4OPnuKpV2cQgaF0tFgpCoimXYcju4dWzBq50kD2lWTLrLcob0bmTi+w3W5uJvDGMm6//xFcgYuFapOLJlCWZVM0WuCnpgp40erOhCscHs+1FMxGC2q+VGMw7TKWy9SfPzmZB4Uju4eAZlfR4fEwiPv8+bBcb8KVeiqkG2WsLBXnxi/1K9MLTcHVlOuAhI+P7hleljXSK9Zeu5mOpTCuje12HTc9i8a4PLopPLGPNhZTaE53bKQx46TqB7iO1B/DYiCxbvFP5smXPUYHwswTzw+YzIcZMrsG03WLav+OxXNV/QBHFo85X6pR9Rczcbwo7XHvSOvVsI1B43TCxQsUUfAJF1zBYtC00S/eSwt/tmO98m6wna5jnIbcDhP4HmW9hGe1Jh8XC+XQReIHJByHD7z5qmXHaAzYpVynbsE3CuZgyq2Pt1jx6ouMJvMVskmXobTLQsWn5pWo+mH64YX5CiLCUCZ0/cQWPIQFumKHjiD1VMgL8xXeemB02Rgbg8YVL6inOIqEWTYo7BlO14O2cU2Zp16dQYA9I5k1FwJb7xvwdqxX3g369Tp6fkDVD6jUAipeQMXz8YOVPTAm8D3KetVUb7xJvDxd4N5PPUkuSie89cAIf/HMBfxAySRckq7wieMv8UePv9LUwagx42Qsl6r74PfkFgUz5TrUfJ/pgkexIXtFCJt4V/2ATEIYSGcYiVaZnp0p8fJ0KMyuEy6MGhpMoKr1Al2OE2XzRIuKPG2d0jeUTnByshDeCFRD94yCi3JdlLZZqIQzlMZgaqBhptC52TJX7aCjjk+tru3l3ICX3iDicfVCSYOtTD/UffcDpRqJeMULRd0Lwu/VbLHKyckCJy8UeGEyv+JxTOA3kLVYfI3TzLiuesXzOTtT6qjkLjTfJPLlGtOFGopSrgVM5ss89eoMOweTjOUyzJdqYWMMoFjxOH2xwE/9hycINBRMUcgkXQbSLkcmcqgqC1WflOuQdJTnL+QJNKwc2UigiyJdrCkT0XjmS7WwTk20T9IRcgmHHdlkvUDXYNol5bpcLFTqrqHDOwebgpnxtcyXa+ENYMmMVRU++r43NF2vu489Xr8ucVlgGhY0dWLtXekNuNUN4qGnXquXSrjSwGivxBU2g61W9z0ItMEyDwU9dr1MFyq8cKHAycl8KOgXCkwtKW63EibwG8RaLb54mtm45N6JUho7tRQbbxJT+Up4DKSew+35AefnKlyMUhUdERKuUPYCKtG2GAXKNZ+BtMv73riHx05d4uJCnsl8hdGBRX9nqyljELlIVLW++OZioULScRAnzt4Js1tGB9N88effXr9eCVc4NDZYt8I++r43tLyWr89XltWwAQiAX/6Tb3H1rlxd7J49N4fnB9QCxQFqQeinrxB2fBrKJPiVH7ix42sbk026nLwwv6zZSqt/p3Y3iMdOXbriQGAvxRU2i16t+x7PTmMXS6UWirmqMpWvcHKywAsX8tHvApcWqi2P4wgc3DXIkd05jq1wPhP4DeD4iUk+9JlvsFD1yCRcxofSDGWSK1p88TRzcr5MqI7Ny+I7sRQbfZFVP3R5xIHLkxfydUF0nXBxU+zaCBRokV2lhMHPTxx/if2jWfKlGhVfOT9fadpnKQlHSKdc5soeJ87nSTjhTQRAfEg4i1k0jTOUdlZYowUOoTi6An6LcwOcnatwMSpDMDlfpuorrgNJ18Hzg3oNHIn+10leWSs/78VChXzFX1bCuJWwdjMQuB7uvX5gs2cxqpFl7i1a5zVfCYKA8/Pl0DK/kI8s9ELbjmcJRzg4Nsj1u3McmRji+t05Do8NMpBOkEo4JvCbyfETk3zkoafJlz0UKNZ8zlwqcvXOAXLpxIot8+4D7v3UkyiQdIWxXIbhbBJV7UgIGn2RDlCJFD3pQKXBl6FBaBEE2toCr+9HfINQXp4urhrgAUi7wnA2wWS+ynDGZb7kUfaaZwa1ANQLs2iE5hlKqy9kK3FsdUNqejoAV6T+vv0AEtF7hvC81+waaHvj7cRfPlOssXMw2ZGwdjMQuJ2ySNqxGbOYJp+5F1D1Avwg4NxsqS7mJydDMc9HDd2XknSFa8dzHGkQ84O7BhlMJ0gnHFLxj+uQcFevNGMC32V+/QvPMVtcvDOrhpbm67Ml9u8cWPELfefRCW45MHrZQtBoBb86XUQIy/a6rkMtCtgknDCfvSEuWhf7dgih37ATAoWFis9Q2qXqaUs3CkSFvVyHPSMZXKd5hrJUXJeWOACaxt+KOP2y+fQaLoIinMXEDTmWimGn/vK5Uo1dg+mmM7QT1m4GAvs1i2QtdHsWU6tb5n7df17zA87MFCNfeSjmL04WKFZbzy3TCYdrxwc5snuI63cvink25ZJOhHWD0pGYO0tiW51iAt9lTk8XcQQc16kHTpTQml7pCx2L2rPn5shXPEQhm3IZyiSWNb9YidgKvvXX/opCxa9bFRAu3ReRei78s6/NQQefo4QrYb9SVnZnpBzYPZKhWPXCNnyrHPyqHRmGMoszlOMnJvn1LzzHyakCScdhOJPgyVcu1d1MKVfYO5Kh4gWrulVaPR8QfsmCqFhZzFIx7NRf3mqBTTth7WYgsB+ySK6U9ZzFeLGYR1Z5xQu/Ry9PL3Aycq+8cCHPS5OFuutxKQMpl+smchyZyEWCvijmKbfZMm9c8X2lmMBvEOHiIAcvCNDIamy3dDq2GGu+T6nq1y3qYtXHC5SfufPAmnPhCxWfQMOaL3HdmSBQfF959txc3WIfTDocHh/ihfPzVH1tKYydLn6uBpBLJzg3FwaLZQUZFqFuQTfm1U/my7gSNu1uKidM6C565VKps8EsPR/xTCQgIHQjxQXFlophp2KxVmHtViBwq2WRdIPLncX4gdaDn7FlXqp5nL640JTNcuriAtU2Yj6Yduvulfj3NS3EPNmBi2UlYiMwOX7wTe32MYHvMofHBsP8bFUcB5Li4AfK3pEMDz56il/+s2eWBYBii3G64OE4QkIcAgmX7O8ZyfAX33693nWok+DRg4+eYudgMkqTBHEg8JT6xzPSXUdgoRpw5tICe0YyvNpGPAMNx7K0mFcrzs4UicthrLS3QF1g50o1ilWfiucTBK27PHWK60iYi7/EpZSKygy7jsMH77x2xdTETsWil4R1rTePzQ5Irjed3GyDQJus8ooXUCjXOHVxYTGb5UKB0xcXln1+YoYzCa7fPVT3md+wZ4gDO7Okky5p162LuXuZLpZ2NLoN0aC1Qx+rRdN14kbU+bKHF4SrRdPJ8O7d2C5uvlRj12CKQtVnKl9hz3Cac3PlSKDCcrW+KleNZDg7Wwqnd20KcC39ksbFvgoVj/Nz5bbTSDfKIAl0ceFROwaSDiMDSc7PVVZ200SVISttzhkznEkwkk0yGLXJmynWcJ3QSr+Sj2jKFRyRpvecTjjsHclwfq6MArccGF1WLKzxGkIY+K16PvmyR8ULc/KXthHcqvRrca7GYmv7dmT58e85xHddu6tuoc+Xa7w0WWhKTXx5SWexRkYHknX3yvWRmO8fDcU8dq+kE+vrYmlHozvwL+/70WJ16pXBVvuZwG8AS6v6zSxUqDV0DsqXa5ydKZFwhevGc/VWcQlH6u3vgkDrzS0ai3RBWEMm6QgXF6oUKl6UihW6VxJROuBINkkm6XJutkwtchOtxDU7s8vcH3HVSQhvAG87uJNnXpslX2mXoBjiAKvEQPkX7znCh95zff2De36ujBelbtZWuNGsFgeIcRtmAcNph2It/KOxIuZK5X6/dXaWTxx/CS8ISLvhzS3puvVA68nJPFUvIOkK1+8e3lIWcL8V52qVnjhbrHFyMt+UzXLmUrGtmI/n0hzZneP63Tmum8hx494RrtqRbfKVpxJX5mK5EhortK4k8Oai2QCWTpfjf5yYqXwldEMEioiweyjDa7OlMC9dIYja3w1lkkwVqk1FuoCoIfQ8cYS00fL2AlAN68IknDCoupq4C3Burrxs+9LXnZzMN5UlaMdKewihZfTYqUt8iEV/d9wez5HVJTzZgbvIV6LprJKvhp2bdg9lmhqGrNRyD2D/aLZJBKfyZT5x/CV2DiaZK9ZAoFSD0xcLbVPyetEVstXTKpemJ4arP/NN2SxnZ9rHavYMZ+pifmRiiBv3DrNnJFMX83TCXXcXy5XSym3YChP4TWDpP068CCnO5AhFJ1xANJxNUI16gB4ayzE6UK2n/MWUan4oYM5yXzNAwnFQQks4IaunQbpOuOI10eZ4IcpcsdZxuiQst7YdCd/zfLnGU6/OcPzEZFO3JggXD3lBKNCtZgKOdJ6ymXIdDo0N8tz5PNeN55qm0nHLvQNtWu4pLBPB2O02XwpjJfFMK1/22DOSaJlL34srTLdSWmVjemLFC5icL/N8LOZRAPT1FsZJzL4d2bpVfsOeIW66aoSxXHpdUhI3ksYYw0qYwG8CSwNAroTlcMdyiznU1ag/aCbpNhX+ikViafAoIaHPvJXWVf0Axwlrqu8fzeL5yquXii3tYofQyg80dI8sJc4ZT7pOVKMmXO5fXwG7CrHIxzVrYhdUXIIhdpMUqx5DmUQ9mBv4fn01ali7Pk3F85kvexSrPkLojlppQlGq+Su21lut5d5SEax4obumsYSyRCWPW1nAvbrCtFfTKhvTEyuez+uzJU6cbxbzyXzruiwCXL1zgCMToWUei/nOwVTdKt9MF8uV0hjQR5y2Ot41H7yIXA38EbCH0PA6pqofX+k1W8kHf6VT7Ua/fC6dYKpQqQddpxcqTOarjOdSjOXS9cySbNJlMl/BDwJSblz4KxT/jzz0NBcLretWxFw9mkVESLrCfKnKxSirJp1wGB1IhjXUgbmoFntsx8QiG+eMu1FDbUWb8uFHBxJMFVovt4ZwVWstugml3EgQEQKUHdkk+YqHKoznUlxaqFH2/LCEb9Jh/+hA0zWKRWgg6TBbqjFVqNbjDiuRTjh8/xt387UXp5sC30OZBP/47des6oOP2wsOZRL1lavzJQ8vquUTx0r2jGSW+bBbdbZa2nxks9js7k+N6YnlWlhU78T5+SYxn16hLss1u+Kl/DnesHeYG68aZkc2RTrZ+arPrcpmNfzwgP9ZVZ8SkSHgSRH5a1X9ThfPuSE05qnPFWu8PlfiqVdnOs6qaJepEX/B4pWf+bLH9EIVR4RqtJgnzgqp+gEDLC54aleUqBENAn7tf7iZBx89xVyxSirhRP7LgLlSjZ9+57XcvH8HP/3HT9YjkomoEiQsljfwUdJJF1Wi5dPCbNFbUdxhsVRCypF6G73Yop2N2vP5QeiaCqIlpn6g5Cs+L00VSDhCseI39YL97JNno9z+zgwVVeWR56fqxxeRev2Zm/fv4Ob9O1q2Anzoqdfqvvay5+MVlVuuHuGpM3PUono2rhPORoYyyZYWcC+7QjayOFecnhgKus8rl4p859x8U2ribJu6LK4jHIrFfPcQb9g7xI1XDTOSWRTzreBi2Si6JvCq+jrwevQ4LyLPAfuALS/wDz56Kqp9XotEKsxt/8Txl7h5/44Vvyjt/LB33bIPCIWmWA0XJSUdp14/JZYv1wktEQmUQiWc3s8sVFZ1jyQcmFqo1cf2Cw89XW+KESiUagEPPHKSD73rCLty6Xo7v2LVX2YVB4TClHSEVMphvuyHAUzoKCWyGiyuRK36AaUFH0cExwlnBYkoRTNeEAZhumTVV1Ku8IY9Q5RqPg899Rq5lMtLUwWi5KIVkWh8cyWPlCvcsGe4/lzsKvn0PW9f9u8XFzcbyWbqbQcvFso8dWaO0YEk+bJHqeYTaJg+emgs19IC3kxXyHrNONf6+nr1xFpAyfM4PbXQLOYr1GWB0IBIOGGV03/2PQf5odsORKu513/VZz+yIT54ETkIvBX4eovn7gHuAThw4MBGDOeKOTNTZK5Yq5fwhTANzwt0VX9q7If1fOX03EI9wPrAIycZjxpQx4HNSvRcK+GSyLo9O1NkMl9ZNV3QC8APgvoX9dJCddniIS+A3/6bk+wbTpNIJRjLpXnlUvtMilqgzBQ9XAd8FZKOrOkLV20YQKBadxH5Da6Wpe+p6iunLy4wmHKZbWjr1ynxv9fSrJuVskZaZZnMRfGAdMIlj9e0oOrMTLGeedP4WdishVBXGtzt9PWxmFf9gGLF58XJPN95fT5azh9a5gsd1GX5xisz1PyAkWwSx3FwJDQovn5qhg+/54Z1uirbg64LvIjkgP8MfFhV55c+r6rHIKx4edttt/VOUv4KXD06wOtzpaalxrG7opVINFo/U/kKwxmX2VJYm6XRQn99bnnAqNUF8QOt17WZK4WZLJ2s9nQd4WMPP8v5uVLbfVXhQqHKSFaZbzNNXj6ecKQahK6bTnPTlxK/rn3mTkjFC9oKxWrHrnjh65Sw9HGcrbOSq6SVayVucXhurhTW2NFoZuIrroSpkvd+6kmGMommIPlm1Cm/0uBuu9f/7lde4taDoyxUPU5eKPDMa/ORiyXPi1MFym0i3tlkVJdld4437BnmjfuGuX73UNiAJeHwzt/4Mrty6WVZTlslbbOX6KrAi0iSUNz/WFU/381zbST33nGYp16dwQ/CL7Nq2C90KJNcJhJLrZ+L+QoXCzVcR3DcxUySmNXEsVRrFrb5Faa3y1Cl6vkrWr1O1PJuPJdmepWg7VJi183l0slNoVXZgbUcu/Glr1wqknaF0cHUigXc3nF457JFTgnHIQgCHAl9vrWoxSAQxhCimECx4m16OuRKee4PfOkFPvm10yxU/XpcY2kc6cxMEQflpdkSFS8g4QqZpMvrsyV+6MHHeGlqtbosOa7fHfrKb7pqhOsmcmSTbttVn70cq9hqdE3gJfyX+33gOVX9rW6dZzO48+gEP3PntdGXPs6qCLNQZotVbr//kbqfcqn1s2ckw8vTRbxASbgtluFfrvnbAZ6yor8TwptV0pGwLVgXx3K5rPdwKr4yU6y1LeB2/MQkf/T4K2FRsgCKQUCtUOUH37SHh791njBqQL3scMKRaN1C2K0qXrG81GLeyAVP7QRTg4CPP/IijoQxmlLN5+OPvEgQKD/1zmuZL9f4zrl5ylW/nsGigO8pFS/8HD33+mJP0OFMgiO7hzi6Z4g37B3mTfuGOTyeW3NKYq+mbW5FupkmeTvwVeDbLK5P+Veq+hftXrOV0iShObUsrqEyvCSNr1j12DOcIV/26r1FvSiAmHTD3PQ4QLca66G3CUc6soB7UNu7ggCHxwfbLst/729/hRenFqLUUOqVOK8bH0REOH1xAV81nM05Um+JmE449Xz/w+O5pnTIja790u58FwsVqp4fLoSLunnF6Z7X7BpYtalLNunwQ7ddzY1XDfPGfSMc3DVAOuGuS0riZqdtbiVWSpO0WjTrRLt6HnEZgrmSV28yfblXPMw6ubJ/r4QTBoRXKR+z5YkXZDUGbFfaV4CrRrPLrOkbfvkLaFR1Mqbq+QTAzoEkhYrPzsEkKdfhtdlwBaUroUUjSL3GfWNtl82o/XL8xCT/9isvcWamyMRQhu+5bhef+PJL9Q/kStcoXnHsBRrVzhfGcilA+NpH392V8Rqds1l58NuG4ycmeerVmXrjiMFUuOw9LncbT1+u1Cq+HN/zUvwg7CjV7xa6Nvx0um+ltugvv+vsLI+dukQl8q2rBiSiHq7xPXbvSJaLhQqXFmpRMDW01OPerDsHk+TSoXum0cXQtmH3ZL6jht1LaeXuueP6cap+wFS+wrfOzvLsuXkG0glEhG+emeUbZ2YX3/wSXIGf+HuHuemqYW7eP8JH//O3mcqXGUwvjjm+IRm9jQn8FRJPf2MrMM7wiF2OjaGnKxXUddD3+hg6KdC11WlsMNLJDW2qUKNYDTvc//aXTuI4i68LK1r6xO1kM1GAcHwow2A6scz6XsnF0MonPr1QIV/2OmrY3cjxE5P8yp89Q8IVskmHl6byfPDT3+DgrgEuLlRXrMsylE6Qr3hhx7HI/aTAz737SFOg9affea35xLco5qJZA8dPTHL/F09w6uICAId2DSDRqlI/0LAUrx9csStmo9gKY1wPLvd9xgsiGxtzK6GFe/XOgXoHKlXl/HyZIxNDHVnfrXziZ2dKjA4kGW+wilu5beJSuOdmSzx9Zo77v3CC2VKVmq9tZ3hCWAkzrMcyzJv27eBN+0YYG0rziUdOtsyiadVkfKWmKMbmYT74dSBu3DFTrNW/+J4fdkVyBTJJl1w6wWR+5QYYxtYjtm4dRxhMuQym3fqKVghXtV5aqLF/NFsX7LlSjfFcmnzFW1aO4sxMkVzKRUQoVDz2jw7wwoV59o5kl9WpmVmo8O9//Lv41mtzPPvaHI+9NM2Lk4UVa+QvJZMQfvdHb+V737C7o/37tQFIv2I++HXgwUdPkS+HKxYdETw/qLtfAg394zPFGulEWF1wPdwpRu/gOPC2gztbpvBdWqgxOpCsu1w8X5kt1iiUPa6byDGZL/ORh55GCEtB78gmI9EM+LUPvJE7j05w97HHOT9XxHVcSjWPci1goerh+cr3/fajbcclQDrpkEm4beu3lD3lYw8/y1c7FPherXpprB0T+BVYugK15gf1fF6/YeYTrpBcXOhi2r4xpJywsXdjt6ZuECikXYd3HN7Jg4+eolj1mC/WqEUZOmF5hcVoy8VC1MBFwwYuA6kEr82WQGEwneD0xQXKNR9X4KOf/xbfd+MeXpleaNlkJSbpCknHwXEgl0kwkEiQSTlcWqgyHXXyWomVGl4sZas3ADEWMYFvw7IVqIWwoFctKiuw1EJv9PPmUg5eQNvep8b6UA2itE/Hwe/ytb7l6pF6KWFB6wXTklEN+rCSZok9I9llDVyCIKDmhbO6V6eL9awdX8NVr//h8VeWnU+AG/cO8Q//mwO8ef8Ib9gzzHt++yvLyg2P5dJcLFRIRAXv2t3n1nL/s5Wk/UP/Fkm+QhqnqXEbvbjeSytxh7AAGECxFqzJR2qszErly7wAjkzk6Ga57x3ZBM++nq9/HqYXavVx+boo5NMLVXzfB5Sar4jAi5MFvnM+jx9lqAQsF1tHYCiT4KqRDNeOD3LT3iEOjw8ynE3xT7/7IG89MEom5XL16MCyUhClmo/rOFw3nuOaXe0FOJtc20rSeJGeqi5L8zS2DibwbTgzUyTb0L5tOJskESnNUsGp51xH39xAWXEFoLE2VruSpy8u4HfRgC9WPObLHq/PlsiXa/UbfDyTCzSsghkofOd8gZofPleohB2kVspjSLlhj9xrdg6wKxf2g3Ucp6VLpJ3wHh4bpFTzGcokl7lWIPy8/vQ7r+34/d55dIL73n8TE0MZ5ko1JoYyFmDdopiLpg2tpqk+MJhyOTyeA+CZc3OrNrA2us9aK0uulWp086h4AWdalE9ulZ6YTjQvisqlHErVgMaRurJ483r+/DxIWLpifCiN68gyl0i7csNAPfC7fzSLI3CpuDjL2D+a5eb9O9b0njej6qWx/pjAt6FVtkTc2m2+VONioWLivs1YLZC7dzhNLp1kMl+iVAuYGA7r+9cCpVgLliXk+woa3Ry8AJJu2FQ6zon/lR+4cdk52glvo/DvHs7gONLU3rAXmnsbG48JfANLF3fcdcu+psUdN+8b5s+eft1SILc46YTD4bFBnjufX33nDsmlXSYL1ahgV6jj04Uavmq9wctSYldf0hV8DTuDVf2wHO+uwVTb6patqlA2Cv/dxx4PWzpamuO2xwQ+olXXmoeeeq3ue4wXOpm4b30qXsDBsUHOzZWZ67CpCYDL8iCpCLgSCrgjMFOs1fvnioR7qi4WTWg04kUWA/aZhFN3/cWVJ5fSaWclS3M0YradwMcW0Ldfm6FUU1SVXDrBcNolnUowkFp0wVS8gA995hs88CNvrS90MvqDLzxzftV9lpY4aOXpVwUvajeYdML68vUgfJQ540etCJcWP4vF3Q+UvSPp+vZ2KYmdLkCyNEcjZlsJfGwBzZeqFCqLaRfz5TBLYjhdY3Qww2uzpXpO8XzZiyx3bdu1xth6ZJMO5VqwYoZOp5O1xYJki00/ag056Uvz0x0WxT7hhK9OuIKqrljIq1PLfKs2zNjIJijbhW2VJhlbQO3a3M1XAl69FHZbavxCThWqTC/UMHnvH0pt+oVeDktvBK4rNDYwUuCanWF2iwAIDKRcrtk5wNE9Q+zIJjtKSWyXB98y22aLpTnGxtfSaprHT0xu9tC2NNvKgo8toJX86OZi3z6s57+1SGgt+QpVLyCbdEkEiq+KK8LFQrVejjftLvrbi1WPI7uHO2r0sRbLfKulOVr9m+6wrSz4VhaQYawHKSfMfEk4QirhcGhsMOyAFIT+eS8Is2kgLGExX6queYXoVrTMO2XpwkKwwPB60M2m238A/CAwqapv7NZ5VqPRrzeUTjBXqtV7axrGeiDR/2dKNUYHEsyXw3LBAymXcs1HJKxAiggJDQhUOT9f4ZYDo2v2M281y7xTLDDcHbppwf974L1dPP6qLPXrxUWgsgl31dcaRisE6pZ4TMIRnKjhdr7ic8uBUb76i+/igR95K76GmVqqoTUvAlePZpkYSvPpe97el2J9OVj9m+7QNQteVR8VkYPdOn4nxH49z1dOzy3U3TOWy25cDvWFSY6DH/gkox6tjiNo9F+1pswWq9x+/yMMpRM4ELVGVFIJh6uGsyRcWbWfaTzzfOHCPDU/fO2RiaG+zSxpV4ahH9/rRrLpQVYRuQe4B+DAgQPreuwzM0VcgXNzZTRQE3bjskm6wlUjzeI8mS/jB8pUvkLVD0DDdMeqH+AKnJwsoBreGFxX8IKAc7NFAsKaMw986QUeO3VpWVpgPPOsen4946tU9Xl5utDXJQf61f20mWy6wKvqMeAYhC37Lvc4rXJoh9IJnr+QN2E3rggBrhrJIgKvz5Z4ebrIkfFB5ks1hrNJDkXVHM/OlNg5GHZ2OjVVwHUEVOplpGu+ogJXj2aYKVb5+CMvMjGUYtdgumlVajzznC54OITun0CV+ZLHnpGEZZYYHbPpAn+lxI2wX5gs4KA4jsP5uTJ/e3q6q11+jP5FWKztL4T9dkXg3GwZRckknPpCpqoX1LsleX7ApUKVi4Uqnh9mzTgi+KqkXCcsDQwMZ1OcmirgCMyXPMZymaa0wDidN5wJRCFcgaofWGaJsSa2lMC36vT+0FOvMTlfBlVqCgS2HMm4fFyB8aF0vVPSVKHKcDbB5Hwo7hB2URpIJVioeEzmK+wfzeL5Aa9eKlHxlZQbrmiq+krCod6nt7HLU9UPcCLRjonFO84oSbkOXtQ4RKPGIpZZYqyFrmXRiMingceAG0TkrIj8xJUcr9VKt08cf4mq54eZCuszbGMbsyOb4PrdQxwayzFXqnFoLMfPves6Du7KRcLtcNVIluGoXEC+7OEFYdXGi4UqboMrJhFZ3l6gjOVSuCIEGt4cIBTroKEbFCymBcYZJUOZBAGKFwQEgTKcTVhmibEmuplFc/d6Hq/VSjcvCMiXPVKuw4JvC5iMyyOTcFANG4dM5ct84cN3ND3/IcISvEvztCteQLrBIk+4DhIEeAEgoeXuB2Ev1kNjg0wVKvWaM0OZRH12sLQGTWNGiefPU42yaA7uyllmibEmtoyLplWhpbTrUPEC9o9mKUbNjA1jrYiEgdAgUKptAjetygS4jjAyEH4mY3eK4zgMJITD4zmKVY+JoUy9DEHsYjw7U+TQWI5/9LadTf0GGsXbMkqM9WDLCHyrlW4jA0kuLdRwHSGTgJJV8932xJUdJSrslU26lL1gxR65ioYrmxVSidZey1Z52h9481U89NRrFKseY7kUr82WQWHPcLrlQp1G0bbKicZGsGUEvpUFlXRdfubOA3z27141cTeAUNgPjGZ5bbaMEAZMz8+VVxT4SlRZMukK47l02/1aWdU3799RF/3rxgcREQqV0HJvJ9qdNu4wjCtlywh8KwvqHYd38ntfPUW+Yv53IyQZBS8HUi6D6bAkRaum2DGNpQcChalCheMnJju2tC/HlWKVE42NoucFvlVq5GOn4IUL8zx2anqzh2f0GId2DfDFn39n3Uo+P1fGkbBeTKPQOwKj2SQLNR8/CDNkxnJpEq7UhbZblra11DM2ip4W+KVfsJenCzx2anpZKzXDiJEoPTGe8d37qSfDGuwJh325NMPZZL3nqQJ7d2TrrwGYL1V56tUZbr//EeZLNQbTLiPZsDTBelnaVjnR2Ch6uh5841S2UPF4fa4CmLgb7ZkqVOqP7zw6wS0HRjmwc4DD47l6/nospkv7A+TLtdB3L7Ajm6RY9bmYr5IvLzbAXg9L2yonGhtFTwt83ARgcr7MK9M2fTVWJ1/2mtq8rSSmS587P1cGYPdQBhEhnXBAYCq/eNNYD0u7nxt3GL1FT7torh4d4OXpApP5ilntBhDViNHWszjXEUYHkk0ulNXK0DY+p8C+HZm6pT+WS3NurkTZ81dtiL1WLM/d2Ah6UuAba2FfWqiZuBsMpV3yFb9lJy4hzF/fM5xhKJNY5kJZSUwbn4tXq8YMZ5NUPJ9iNezQZDXKja1Gzwn8A196gX/z5RejJgnGdiEOc7b6Vx9Iurxx3w5OnJ9jruQRaJgFI0AyaoQRU6x6l+1CabXWIpVw+fV/cLOJurEl6SmBz5c9futLJzd7GMYGk3SFD37vdQD8zt+cbKrf7wj883ce5ub9O/jYw8+yczBdF9/5KBOmUZCvxIViXYWMfqOnBP7sTJHxzR6E0XUcCRckeYHytoM7l4noJ792moWqz2DK5SdvP8SH3nM9wDLx/ZUfuBFYX0E237jRT4i2cmpuEum9R3TvP/mdzR6G0WWu2TlQb30XF+IyDOPyEJEnVfW2Vs/1dJqksbUYTLmM51L1uuitmMilSLhied+GsQH0lIvG2FoIMJhyKHnKeC7FWC5NqeaTTSW465Z9fO7Js5ydKaFA0oE9I1mAFQtxGYaxfpjAG024DqQTLqphbfSlfnBormu+v14faHld88bXGIax8ZjA9wkugLCs0XjCgZFMgpKnlKNl+XHYJeEI144P8v1v2tu28UQrWgUiP7SO78UwjPWhqwIvIu8FPk6oP59U1V9faf90wuGanVleuVTq5rC2JOmEgyNQrgX1XPHhTKJuXS+1qtfqAjGBNoz+o2sCLyIu8Ang+4CzwN+JyMOq+p12r3FESLgOrgM7B5JMFWrtdu1bBBjKJHj30XHOz1evyKo2DGN7000L/m3Ai6p6CkBEPgN8AGgr8H6gTAxlSDpCLVAqtYD5LdLMI+kKnq/LVmKmXOFnv/c680cbhrHhdC0PXkTuAt6rqj8Z/f2Pge9S1Z9dst89wD3RnzcAzzvp3LA7PHYAVcVxE+I4mxMrUA008GuAqO9VgoXZ80GlMN+FM40BF7tw3H7DrtPq2DXqjH66Tteoass1ot0UzlbZ0MvuJqp6DDjWxXH0PCLyRLuFCsYidp1Wx65RZ2yX69TNhU5ngasb/t4PnOvi+QzDMIwGuinwfwccEZFDIpICfgR4uIvnMwzDMBromotGVT0R+VngLwnTJP9AVZ/t1vm2ONvaRbUG7Dqtjl2jztgW16mnio0ZhmEY64cVGzMMw+hTTOANwzD6FBP4LiAiV4vIl0XkORF5VkR+Ltq+U0T+WkRORr9HG17zSyLyoog8LyL/XcP2W0Xk29FzD4jICsV4tx4i4orIN0Tkz6O/7RotQUR2iMhDInIi+ky9w65TMyLy89F37RkR+bSIZOwaAapqP+v8A+wFbokeDwEvADcC/yfw0Wj7R4H7o8c3Ak8DaeAQ8BLgRs/9LfAOwnUFXwDet9nvb52v1b8A/iPw59Hfdo2WX6M/BH4yepwCdth1aro++4DTQDb6+7PAP7VrpGbBdwNVfV1Vn4oe54HnCD+EHyD8shL9/u+jxx8APqOqFVU9DbwIvE1E9gLDqvqYhp++P2p4zZZHRPYDPwB8smGzXaMGRGQYuAP4fQBVrarqLHadlpIAsiKSAAYI19xs+2tkAt9lROQg8Fbg68BuVX0dwpsAEFcH2wecaXjZ2Wjbvujx0u39wu8A/xIIGrbZNWrmMDAF/LvIlfVJERnErlMdVX0N+E3gVeB1YE5V/wq7Ribw3UREcsB/Bj6sqivVsGlX1qGjcg9bERH5QWBSVZ/s9CUttvX1NYpIALcAv6uqbwUWCN0N7dh21ynyrX+A0N1yFTAoIj+20ktabOvLa2QC3yVEJEko7n+sqp+PNl+IpoFEvyej7e3KOpyNHi/d3g98D/B+EXkZ+AzwLhH5FHaNlnIWOKuqX4/+fohQ8O06LfIe4LSqTqlqDfg88N3YNTKB7wZR5P33gedU9bcannoY+CfR438C/FnD9h8RkbSIHAKOAH8bTSvzIvL26Jj/Y8NrtjSq+kuqul9VDxKWsXhEVX8Mu0ZNqOp54IyI3BBtejdhyW27Tou8CrxdRAai9/ZuwriXXaPNjvL24w9wO+HU7lvAN6Of7wd2AX8DnIx+72x4zf9CGM1/nobIPXAb8Ez03P9NtPq4n36AO1nMorFrtPz6vAV4Ivo8/Skwatdp2TX6X4ET0fv7D4QZMtv+GlmpAsMwjD7FXDSGYRh9igm8YRhGn2ICbxiG0aeYwBuGYfQpJvCGYRh9igm8YRhGn2ICb2wLokUtXxKRb4rID4vIv7qCYxXWc2yG0S261pPVMHqMtwJJVX0L1EX6f9/UERlGlzEL3tiyiMigiPwXEXk6avTwwyLy3qgxxteihg1/LiITwKeAt0QW/OcIS8t+U0T+eIXj/6mIPBk1krhnyXP/l4g8JSJ/IyLj0ba3iMjjIvItEfkTERkVkTeIyN82vO6giHwrenyriHwlOsdfxnVTDGO9MIE3tjLvBc6p6ptV9Y3AF4HfA/4+8PeAPQCqOgn8JPBVVX2Lqv5DoBQ9/tEVjv/jqnor4fL1D4nIrmj7IPCUqt4CfAX41Wj7HwG/qKo3A98GflVVnwNSInI42ueHgc9Gxej+DXBXdI4/AP71FV8Rw2jABN7YynwbeI+I3C8if4+wXOxpVT2pYQ2OT13h8T8kIk8DjxNWHzwSbQ+A/xQ9/hRwu4iMADtU9SvR9j8kbNQBYYehH4oe/3D02huANwJ/LSLfBH6Z5kqGhnHFmA/e2LKo6gsicithIbf/A/gr1ql+t4jcSViG9h2qWhSR40Cm3VBWOdx/Aj4nIp8Ph60nReRNwLOq+o71GK9htMIseGPLIiJXAUVV/RRhR5/vBg6JyLXRLnev8PJa5CZpxwgwE4n7UeDtDc85wF3R438EfE1V54CZaCYB8I8J3Teo6kuAD/wKi5b/88C4iLwjei9JEblp1TdtGGvALHhjK/Mm4DdEJABqwE8DY8B/EZGLwNcI3SCtOAZ8S0SeauOH/yLwz6OA6POEbpqYBeAmEXkSmCN0u0BYc/zfisgAcAr4Zw2v+U/AbxC6kVDVqojcBTwQuXcShC0Mn13D+zeMFbFywUbfErlZPqKqP7jJQzGMTcFcNIZhGH2KWfDGtiZKffybFk+9W1WnN3o8hrGemMAbhmH0KeaiMQzD6FNM4A3DMPoUE3jDMIw+xQTeMAyjT/n/AaCuELPSf/0DAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.regplot(x='sqft_above', y='price',data=df)\n",
    "plt.ylim(0,)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "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": 30,
   "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": 30,
     "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": 31,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.00046769430149029567"
      ]
     },
     "execution_count": 31,
     "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": 32,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.4928532179037931"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[['sqft_living']]\n",
    "Y = df['price']\n",
    "lm = LinearRegression()\n",
    "lm\n",
    "lm.fit(X,Y)\n",
    "lm.score(X, Y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 7\n",
    "\n",
    "Fit a linear regression model to predict the <code>'price'</code> using the list of features:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "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": 34,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6576926701831082"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X = df[[\"floors\", \"waterfront\",\"lat\" ,\"bedrooms\" ,\"sqft_basement\" ,\"view\" ,\"bathrooms\",\"sqft_living15\",\"sqft_above\",\"grade\",\"sqft_living\"]]\n",
    "Y = df['price']\n",
    "lm2 = LinearRegression()\n",
    "lm2\n",
    "lm2.fit(X,Y)\n",
    "lm2.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": 35,
   "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": 36,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.7513302656738235"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pipe=Pipeline(Input)\n",
    "pipe.fit(X,Y)\n",
    "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": 37,
   "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:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "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",
    "\n",
    "Create and fit a Ridge regression object using the training data, set the regularization parameter to 0.1, and calculate the R^2 using the test data.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.linear_model import Ridge"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.6478759163939107"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Rig = Ridge(alpha=0.1)\n",
    "Rig.fit(x_train,y_train)\n",
    "Rig.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": 43,
   "metadata": {
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.7002744273597521"
      ]
     },
     "execution_count": 43,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pf = PolynomialFeatures(degree=2)\n",
    "x_train_pf = pf.fit_transform(x_train)\n",
    "x_test_pf = pf.fit_transform(x_test)\n",
    "rig = Ridge(alpha=0.1)\n",
    "rig.fit(x_train_pf,y_train)\n",
    "rig.score(x_test_pf,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 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.8"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}