House Sales in King County, USA.ipynb

{
    "cells": [
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": " <a href=\"https://www.bigdatauniversity.com\"><img src = \"https://ibm.box.com/shared/static/ugcqz6ohbvff804xp84y4kqnvvk3bq1g.png\" width = 300, align = \"center\"></a>\n\n<h1 align=center><font size = 5>Data Analysis with Python</font></h1>"
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "# House Sales in King County, USA"
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "This dataset contains house sale prices for King County, which includes Seattle. It includes homes sold between May 2014 and May 2015."
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "<b>id</b> : A notation for a house\n\n<b> date</b>: Date house was sold\n\n\n<b>price</b>: Price is prediction target\n\n\n<b>bedrooms</b>: Number of bedrooms\n\n\n<b>bathrooms</b>: Number of bathrooms\n\n<b>sqft_living</b>: Square footage of the home\n\n<b>sqft_lot</b>: Square footage of the lot\n\n\n<b>floors</b> :Total floors (levels) in house\n\n\n<b>waterfront</b> :House which has a view to a waterfront\n\n\n<b>view</b>: Has been viewed\n\n\n<b>condition</b> :How good the condition is overall\n\n<b>grade</b>: overall grade given to the housing unit, based on King County grading system\n\n\n<b>sqft_above</b> : Square footage of house apart from basement\n\n\n<b>sqft_basement</b>: Square footage of the basement\n\n<b>yr_built</b> : Built Year\n\n\n<b>yr_renovated</b> : Year when house was renovated\n\n<b>zipcode</b>: Zip code\n\n\n<b>lat</b>: Latitude coordinate\n\n<b>long</b>: Longitude coordinate\n\n<b>sqft_living15</b> : Living room area in 2015(implies-- some renovations) This might or might not have affected the lotsize area\n\n\n<b>sqft_lot15</b> : LotSize area in 2015(implies-- some renovations)"
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "You will require the following libraries: "
        },
        {
            "metadata": {},
            "cell_type": "code",
            "source": "import pandas as pd\nimport matplotlib.pyplot as plt\nimport numpy as np\nimport seaborn as sns\nfrom sklearn.pipeline import Pipeline\nfrom sklearn.preprocessing import StandardScaler,PolynomialFeatures\nfrom sklearn.linear_model import LinearRegression\n%matplotlib inline",
            "execution_count": 1,
            "outputs": []
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "# Module 1: Importing Data Sets "
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": " Load the csv:  "
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "file_name='https://s3-api.us-geo.objectstorage.softlayer.net/cf-courses-data/CognitiveClass/DA0101EN/coursera/project/kc_house_data_NaN.csv'\ndf=pd.read_csv(file_name)",
            "execution_count": 2,
            "outputs": []
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "\nWe use the method <code>head</code> to display the first 5 columns of the dataframe."
        },
        {
            "metadata": {},
            "cell_type": "code",
            "source": "df.head()",
            "execution_count": 3,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 3,
                    "data": {
                        "text/plain": "   Unnamed: 0          id             date     price  bedrooms  bathrooms  \\\n0           0  7129300520  20141013T000000  221900.0       3.0       1.00   \n1           1  6414100192  20141209T000000  538000.0       3.0       2.25   \n2           2  5631500400  20150225T000000  180000.0       2.0       1.00   \n3           3  2487200875  20141209T000000  604000.0       4.0       3.00   \n4           4  1954400510  20150218T000000  510000.0       3.0       2.00   \n\n   sqft_living  sqft_lot  floors  waterfront  ...  grade  sqft_above  \\\n0         1180      5650     1.0           0  ...      7        1180   \n1         2570      7242     2.0           0  ...      7        2170   \n2          770     10000     1.0           0  ...      6         770   \n3         1960      5000     1.0           0  ...      7        1050   \n4         1680      8080     1.0           0  ...      8        1680   \n\n   sqft_basement  yr_built  yr_renovated  zipcode      lat     long  \\\n0              0      1955             0    98178  47.5112 -122.257   \n1            400      1951          1991    98125  47.7210 -122.319   \n2              0      1933             0    98028  47.7379 -122.233   \n3            910      1965             0    98136  47.5208 -122.393   \n4              0      1987             0    98074  47.6168 -122.045   \n\n   sqft_living15  sqft_lot15  \n0           1340        5650  \n1           1690        7639  \n2           2720        8062  \n3           1360        5000  \n4           1800        7503  \n\n[5 rows x 22 columns]",
                        "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 \u00d7 22 columns</p>\n</div>"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 1 \nDisplay the data types of each column using the attribute dtype, then take a screenshot and submit it, include your code in the image. "
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "df.dtypes",
            "execution_count": 6,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 6,
                    "data": {
                        "text/plain": "Unnamed: 0         int64\nid                 int64\ndate              object\nprice            float64\nbedrooms         float64\nbathrooms        float64\nsqft_living        int64\nsqft_lot           int64\nfloors           float64\nwaterfront         int64\nview               int64\ncondition          int64\ngrade              int64\nsqft_above         int64\nsqft_basement      int64\nyr_built           int64\nyr_renovated       int64\nzipcode            int64\nlat              float64\nlong             float64\nsqft_living15      int64\nsqft_lot15         int64\ndtype: object"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "We use the method describe to obtain a statistical summary of the dataframe."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "df.describe()",
            "execution_count": 7,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 7,
                    "data": {
                        "text/plain": "        Unnamed: 0            id         price      bedrooms     bathrooms  \\\ncount  21613.00000  2.161300e+04  2.161300e+04  21600.000000  21603.000000   \nmean   10806.00000  4.580302e+09  5.400881e+05      3.372870      2.115736   \nstd     6239.28002  2.876566e+09  3.671272e+05      0.926657      0.768996   \nmin        0.00000  1.000102e+06  7.500000e+04      1.000000      0.500000   \n25%     5403.00000  2.123049e+09  3.219500e+05      3.000000      1.750000   \n50%    10806.00000  3.904930e+09  4.500000e+05      3.000000      2.250000   \n75%    16209.00000  7.308900e+09  6.450000e+05      4.000000      2.500000   \nmax    21612.00000  9.900000e+09  7.700000e+06     33.000000      8.000000   \n\n        sqft_living      sqft_lot        floors    waterfront          view  \\\ncount  21613.000000  2.161300e+04  21613.000000  21613.000000  21613.000000   \nmean    2079.899736  1.510697e+04      1.494309      0.007542      0.234303   \nstd      918.440897  4.142051e+04      0.539989      0.086517      0.766318   \nmin      290.000000  5.200000e+02      1.000000      0.000000      0.000000   \n25%     1427.000000  5.040000e+03      1.000000      0.000000      0.000000   \n50%     1910.000000  7.618000e+03      1.500000      0.000000      0.000000   \n75%     2550.000000  1.068800e+04      2.000000      0.000000      0.000000   \nmax    13540.000000  1.651359e+06      3.500000      1.000000      4.000000   \n\n       ...         grade    sqft_above  sqft_basement      yr_built  \\\ncount  ...  21613.000000  21613.000000   21613.000000  21613.000000   \nmean   ...      7.656873   1788.390691     291.509045   1971.005136   \nstd    ...      1.175459    828.090978     442.575043     29.373411   \nmin    ...      1.000000    290.000000       0.000000   1900.000000   \n25%    ...      7.000000   1190.000000       0.000000   1951.000000   \n50%    ...      7.000000   1560.000000       0.000000   1975.000000   \n75%    ...      8.000000   2210.000000     560.000000   1997.000000   \nmax    ...     13.000000   9410.000000    4820.000000   2015.000000   \n\n       yr_renovated       zipcode           lat          long  sqft_living15  \\\ncount  21613.000000  21613.000000  21613.000000  21613.000000   21613.000000   \nmean      84.402258  98077.939805     47.560053   -122.213896    1986.552492   \nstd      401.679240     53.505026      0.138564      0.140828     685.391304   \nmin        0.000000  98001.000000     47.155900   -122.519000     399.000000   \n25%        0.000000  98033.000000     47.471000   -122.328000    1490.000000   \n50%        0.000000  98065.000000     47.571800   -122.230000    1840.000000   \n75%        0.000000  98118.000000     47.678000   -122.125000    2360.000000   \nmax     2015.000000  98199.000000     47.777600   -121.315000    6210.000000   \n\n          sqft_lot15  \ncount   21613.000000  \nmean    12768.455652  \nstd     27304.179631  \nmin       651.000000  \n25%      5100.000000  \n50%      7620.000000  \n75%     10083.000000  \nmax    871200.000000  \n\n[8 rows x 21 columns]",
                        "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 \u00d7 21 columns</p>\n</div>"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "# Module 2: Data Wrangling"
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 2 \nDrop 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>"
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "df.drop([\"id\",\"Unnamed: 0\"], axis=1, inplace = True)\ndf.describe()",
            "execution_count": 9,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 9,
                    "data": {
                        "text/plain": "              price      bedrooms     bathrooms   sqft_living      sqft_lot  \\\ncount  2.161300e+04  21600.000000  21603.000000  21613.000000  2.161300e+04   \nmean   5.400881e+05      3.372870      2.115736   2079.899736  1.510697e+04   \nstd    3.671272e+05      0.926657      0.768996    918.440897  4.142051e+04   \nmin    7.500000e+04      1.000000      0.500000    290.000000  5.200000e+02   \n25%    3.219500e+05      3.000000      1.750000   1427.000000  5.040000e+03   \n50%    4.500000e+05      3.000000      2.250000   1910.000000  7.618000e+03   \n75%    6.450000e+05      4.000000      2.500000   2550.000000  1.068800e+04   \nmax    7.700000e+06     33.000000      8.000000  13540.000000  1.651359e+06   \n\n             floors    waterfront          view     condition         grade  \\\ncount  21613.000000  21613.000000  21613.000000  21613.000000  21613.000000   \nmean       1.494309      0.007542      0.234303      3.409430      7.656873   \nstd        0.539989      0.086517      0.766318      0.650743      1.175459   \nmin        1.000000      0.000000      0.000000      1.000000      1.000000   \n25%        1.000000      0.000000      0.000000      3.000000      7.000000   \n50%        1.500000      0.000000      0.000000      3.000000      7.000000   \n75%        2.000000      0.000000      0.000000      4.000000      8.000000   \nmax        3.500000      1.000000      4.000000      5.000000     13.000000   \n\n         sqft_above  sqft_basement      yr_built  yr_renovated       zipcode  \\\ncount  21613.000000   21613.000000  21613.000000  21613.000000  21613.000000   \nmean    1788.390691     291.509045   1971.005136     84.402258  98077.939805   \nstd      828.090978     442.575043     29.373411    401.679240     53.505026   \nmin      290.000000       0.000000   1900.000000      0.000000  98001.000000   \n25%     1190.000000       0.000000   1951.000000      0.000000  98033.000000   \n50%     1560.000000       0.000000   1975.000000      0.000000  98065.000000   \n75%     2210.000000     560.000000   1997.000000      0.000000  98118.000000   \nmax     9410.000000    4820.000000   2015.000000   2015.000000  98199.000000   \n\n                lat          long  sqft_living15     sqft_lot15  \ncount  21613.000000  21613.000000   21613.000000   21613.000000  \nmean      47.560053   -122.213896    1986.552492   12768.455652  \nstd        0.138564      0.140828     685.391304   27304.179631  \nmin       47.155900   -122.519000     399.000000     651.000000  \n25%       47.471000   -122.328000    1490.000000    5100.000000  \n50%       47.571800   -122.230000    1840.000000    7620.000000  \n75%       47.678000   -122.125000    2360.000000   10083.000000  \nmax       47.777600   -121.315000    6210.000000  871200.000000  ",
                        "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>"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "We can see we have missing values for the columns <code> bedrooms</code>  and <code> bathrooms </code>"
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "print(\"number of NaN values for the column bedrooms :\", df['bedrooms'].isnull().sum())\nprint(\"number of NaN values for the column bathrooms :\", df['bathrooms'].isnull().sum())\n",
            "execution_count": 10,
            "outputs": [
                {
                    "output_type": "stream",
                    "text": "number of NaN values for the column bedrooms : 13\nnumber of NaN values for the column bathrooms : 10\n",
                    "name": "stdout"
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "\nWe 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>"
        },
        {
            "metadata": {},
            "cell_type": "code",
            "source": "mean=df['bedrooms'].mean()\ndf['bedrooms'].replace(np.nan,mean, inplace=True)",
            "execution_count": 11,
            "outputs": []
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "\nWe 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>"
        },
        {
            "metadata": {},
            "cell_type": "code",
            "source": "mean=df['bathrooms'].mean()\ndf['bathrooms'].replace(np.nan,mean, inplace=True)",
            "execution_count": 12,
            "outputs": []
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "print(\"number of NaN values for the column bedrooms :\", df['bedrooms'].isnull().sum())\nprint(\"number of NaN values for the column bathrooms :\", df['bathrooms'].isnull().sum())",
            "execution_count": 13,
            "outputs": [
                {
                    "output_type": "stream",
                    "text": "number of NaN values for the column bedrooms : 0\nnumber of NaN values for the column bathrooms : 0\n",
                    "name": "stdout"
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "# Module 3: Exploratory Data Analysis"
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 3\nUse 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"
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "df2=df[\"floors\"].value_counts().to_frame()\ndf2",
            "execution_count": 19,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 19,
                    "data": {
                        "text/plain": "     floors\n1.0   10680\n2.0    8241\n1.5    1910\n3.0     613\n2.5     161\n3.5       8",
                        "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>"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 4\nUse 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."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "sns.boxplot(x=\"waterfront\", y =\"price\", data=df)\n#Houses with no waterfront have more outliers",
            "execution_count": 21,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 21,
                    "data": {
                        "text/plain": "<matplotlib.axes._subplots.AxesSubplot at 0x7f8f3eb28c50>"
                    },
                    "metadata": {}
                },
                {
                    "output_type": "display_data",
                    "data": {
                        "text/plain": "<Figure size 432x288 with 1 Axes>",
                        "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAEKCAYAAAC7c+rvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAHnRJREFUeJzt3X2UXVWZ5/HvL4lAoiJQFCyoShvspFVaBOEKmbanGyGEwmkTZpa0pGdNbjtM1zQiRB27Bcc1GV8Xrp4lQ5iWNiMZKjMCRkaHwpWXqfDSvgGmEl5iiE5KDFAJDWUlRjAIJHnmj7sr3Cpu6s3cs6+5v89ad91znrPP2buyKnmy99lnH0UEZmZmOUzJ3QAzM2teTkJmZpaNk5CZmWXjJGRmZtk4CZmZWTZOQmZmlo2TkJmZZeMkZGZm2dQ1CUn6mKQtkn4s6XZJx0g6TdJDkrZJ+oako1LZo9N+Xzo+q+o616X4TyVdXBXvSLE+SddWxSdch5mZFU/1WjFBUhvwfeD0iHhR0ipgNfA+4FsRcYekfwAejYibJX0YeGdE/LWky4F/GREflHQ6cDtwLnAqsB74g1TN/wMuAvqBDcCiiHg81TXuOkb7OU488cSYNWvWYf2zMTM70m3cuPEXEdE6VrlpdW7HNGC6pFeAGcAzwAXAX6TjXcB/Bm4GFqZtgDuB/yZJKX5HRLwE/FxSH5WEBNAXEU8ASLoDWChp60TriFEy8axZs+jt7Z3kj29m1pwkPTmecnUbjouIHcB/AZ6iknz2ABuBX0bEvlSsH2hL223A0+ncfal8S3V8xDmHirdMoo5hJHVK6pXUOzAwMJkf38zMxqFuSUjS8VR6HqdRGUZ7PXBJjaJDvRAd4tjhio9Wx/BAxPKIKEVEqbV1zN6kmZlNUj0nJswDfh4RAxHxCvAt4I+A4yQNDQO2AzvTdj8wEyAdfxOwqzo+4pxDxX8xiTrMzCyDeiahp4C5kmakezsXAo8D9wEfSGXKwF1puzvtk47fm+7VdAOXp5ltpwFzgB9RmYgwJ82EOwq4HOhO50y0DjMzy6Ce94QeonLzfxOwOdW1HPgk8PE0waAFuCWdcgvQkuIfB65N19kCrKKSwNYCV0XE/nRP5yPAOmArsCqVZaJ1WDEGBwe55pprGBwczN0UM2sQdZuifaQolUrh2XGHx5e//GXuvvtuFixYwMc+9rHczTGzOpK0MSJKY5XziglWiMHBQdauXUtEsHbtWveGzAxwErKCdHV1ceDAAQD279/PypUrM7fIzBqBk5AVYv369ezbV3l0a9++ffT09GRukZk1AichK8S8efOYNq0ya37atGlcdNFFmVtkZo3AScgKUS6XmTKl8us2depUFi9enLlFZtYInISsEC0tLXR0dCCJjo4OWlpes1qSmTWhei9ganZQuVxm+/bt7gWZ2UFOQlaYlpYWli1blrsZZtZAPBxnZmbZOAmZmVk2TkJmZpaNk5CZmWXjJGRmZtk4CZmZWTZOQmZmlo2TkJk1Pb9wMR8nITNrel1dXWzevNmvGMmgbklI0lslPVL1+ZWkj0o6QVKPpG3p+/hUXpKWSeqT9Jiks6uuVU7lt0kqV8XPkbQ5nbNMklJ8wnWYWXPyCxfzqlsSioifRsRZEXEWcA6wF/g2cC1wT0TMAe5J+wCXAHPSpxO4GSoJBVgKnAecCywdSiqpTGfVeR0pPqE6zKx5+YWLeRU1HHch8LOIeBJYCHSleBdwadpeCKyMigeB4ySdAlwM9ETErojYDfQAHenYsRHxQEQEsHLEtSZShxXA4+7WiPzCxbyKSkKXA7en7ZMj4hmA9H1SircBT1ed059io8X7a8QnU8cwkjol9UrqHRgYmMCPaaPxuLs1Ir9wMa+6JyFJRwELgG+OVbRGLCYRn0wdwwMRyyOiFBGl1tbWMS5p4+Fxd2tUfuFiXkX0hC4BNkXEs2n/2aEhsPT9XIr3AzOrzmsHdo4Rb68Rn0wdVmced7dG5Rcu5lVEElrEq0NxAN3A0Ay3MnBXVXxxmsE2F9iThtLWAfMlHZ8mJMwH1qVjz0uam2bFLR5xrYnUYXXmcXdrZOVymTPOOMO9oAzqmoQkzQAuAr5VFb4euEjStnTs+hRfDTwB9AH/HfgwQETsAj4HbEifz6YYwJXA19I5PwPWTKYOqz+Pu1sjG3rhontBxVNlYpkdSqlUit7e3tzN+J03ODjIokWLePnllzn66KO57bbb/Bfe7AgmaWNElMYq5xUTrBAedzezWqblboA1j3K5zPbt2z3ubmYHOQlZYYbG3c3Mhng4zszMsnESMjOzbJyEzMwsGychMzPLxknIzMyycRIyM7NsnITMzCwbJyEzM8vGScjMzLJxEjIzs2ychMzMLBsnITMzy8ZJyMzMsnESMjOzbOr9eu/jJN0p6SeStkr6Z5JOkNQjaVv6Pj6VlaRlkvokPSbp7KrrlFP5bZLKVfFzJG1O5yyTpBSfcB1mZla8eveEbgTWRsTbgDOBrcC1wD0RMQe4J+0DXALMSZ9O4GaoJBRgKXAecC6wdCippDKdVed1pPiE6jAzszzqloQkHQv8CXALQES8HBG/BBYCXalYF3Bp2l4IrIyKB4HjJJ0CXAz0RMSuiNgN9AAd6dixEfFARASwcsS1JlKHmZllUM+e0FuAAeB/SHpY0tckvR44OSKeAUjfJ6XybcDTVef3p9ho8f4acSZRxzCSOiX1SuodGBiY2E9tZmbjVs8kNA04G7g5It4F/JpXh8VqUY1YTCI+mnGdExHLI6IUEaXW1tYxLmlmZpNVzyTUD/RHxENp/04qSenZoSGw9P1cVfmZVee3AzvHiLfXiDOJOszMLIO6JaGI+CfgaUlvTaELgceBbmBohlsZuCttdwOL0wy2ucCeNJS2Dpgv6fg0IWE+sC4de17S3DQrbvGIa02kDjMzy2Bana9/NfB1SUcBTwAfopL4Vkm6AngKuCyVXQ28D+gD9qayRMQuSZ8DNqRyn42IXWn7SuBWYDqwJn0Arp9IHWZmlocqE8vsUEqlUvT29uZuhpnZ7xRJGyOiNFY5r5hgZmbZOAlZYQYHB7nmmmsYHBzM3RQzaxBOQlaYrq4uNm/ezMqVK3M3xcwahJOQFWJwcJC1a9cSEaxdu9a9ITMDnISsIF1dXRw4cACA/fv3uzdkZoCTkBVk/fr17Nu3D4B9+/bR09OTuUVm1gichKwQ8+bNY9q0ymNp06ZN46KLLsrcIjNrBE5CVohyucyUKZVftylTprB48eLMLTKzRuAkZIVoaWnh1FNPBeDUU0+lpaUlc4vMXuXHB/JxErJCDA4OsmPHDgB27tzpv+zWUPz4QD5OQlaIrq4uhpaIOnDggP+yW8Pw4wN5OQlZITw7zhqVHx/Iy0nICuHZcdao/B+kvJyErBDVs+OmTp3q2XHWMPwfpLychKwQLS0tdHR0IImOjg7PjrOGUS6XDw7HHThwwP9BKli9X2pndlC5XGb79u3+S25mB7knZIVpaWlh2bJl7gVZQ+nq6kISAJI8MaFgdU1CkrZL2izpEUm9KXaCpB5J29L38SkuScsk9Ul6TNLZVdcpp/LbJJWr4uek6/elczXZOsysOa1fv579+/cDldlxnphQrCJ6Qu+NiLOqXvN6LXBPRMwB7kn7AJcAc9KnE7gZKgkFWAqcB5wLLB1KKqlMZ9V5HZOpw8yalycm5JVjOG4h0JW2u4BLq+Iro+JB4DhJpwAXAz0RsSsidgM9QEc6dmxEPBCVpyBXjrjWROowsyblmZt51TsJBfB/JW2U1JliJ0fEMwDp+6QUbwOerjq3P8VGi/fXiE+mjmEkdUrqldQ7MDAwgR/XzH7XeOZmXvWeHfeeiNgp6SSgR9JPRimrGrGYRHw04zonIpYDywFKpdJY1zSz33GeuZlPXXtCEbEzfT8HfJvKPZ1nh4bA0vdzqXg/MLPq9HZg5xjx9hpxJlGHmTUxz9zMp25JSNLrJb1xaBuYD/wY6AaGZriVgbvSdjewOM1gmwvsSUNp64D5ko5PExLmA+vSseclzU2z4haPuNZE6jAzswzqORx3MvDtNGt6GnBbRKyVtAFYJekK4CngslR+NfA+oA/YC3wIICJ2SfocsCGV+2xE7ErbVwK3AtOBNekDcP1E6jAzszw0tLy+1VYqlaK3tzd3M8ysjgYHB/nMZz7D0qVLPSR3mEjaWPVoziF5xQQza3p+qV0+TkJm1tT8Uru8nITMrKn5pXZ5OQlZYQYHB7nmmmv8P01rKH6pXV5OQlYYj7tbI5o3b96wVbS9dlyxnISsEB53t0a1YMEChmYJRwTvf//7M7eouTgJWSE87m6Nqru7e1hP6O67787coubiJGSF8Li7Nar169cP6wn5d7NYTkJWCL+zxRqVfzfzchKyQvidLdao/LuZl5OQFcLvbLFG5d/NvJyErDALFixgxowZnn1kDadcLnPGGWe4F5SBk5AVpru7m71793r2kTUcv08on3EnIUlvljQvbU8feleQ2Xj4OSEzq2VcSUjSXwF3Al9NoXbg/9SrUXbk8XNCZlbLeHtCVwHvAX4FEBHbgJPq1Sg78vg5ITOrZbxJ6KWIeHloR9I0wG/Ds3HzsxhmVst4k9A/SvoUMF3SRcA3gXHdXZY0VdLDkr6T9k+T9JCkbZK+IemoFD867fel47OqrnFdiv9U0sVV8Y4U65N0bVV8wnVYfZXL5YPDcQcOHPAsJDMDxp+ErgUGgM3AvwdWA58e57lLgK1V+18CboiIOcBu4IoUvwLYHRGzgRtSOSSdDlwO/CHQAXwlJbapwN8DlwCnA4tS2QnXYWZmeYw3CU0HVkTEZRHxAWBFio1KUjvwL4CvpX0BF1CZ5ADQBVyathemfdLxC1P5hcAdEfFSRPwc6APOTZ++iHgiDRXeASycZB1WZ11dXcMWifTEBDOD8SehexiedKYD68dx3n8F/hY4kPZbgF9GxL603w+0pe024GmAdHxPKn8wPuKcQ8UnU8cwkjol9UrqHRgYGMePaWNZv349+/fvByqz4zwxwcxg/EnomIh4YWgnbc8Y7QRJfwY8FxEbq8M1isYYxw5XfKz6Xw1ELI+IUkSUWltba5xiEzVv3ryD63NNmTLFExPMDBh/Evq1pLOHdiSdA7w4xjnvARZI2k5lqOwCKj2j49LsOqg8b7QzbfcDM9P1pwFvAnZVx0ecc6j4LyZRh9WZJyaYWS3Txi4CwEeBb0oa+sf8FOCDo50QEdcB1wFIOh/4RET8a0nfBD5AJTGVgbvSKd1p/4F0/N6ICEndwG2SvgycCswBfkSlVzNH0mnADiqTF/4inXPfROoY55+B/RZ27979mn0vkWIAN910E319fVnbsGPHDgDa2trGKFl/s2fP5uqrr87djMKMqycUERuAtwFXAh8G3j5imG0iPgl8XFIflfsxt6T4LUBLin+cyow8ImILsAp4HFgLXBUR+9M9nY8A66jMvluVyk64Dqu/z3/+86Pum+X04osv8uKLYw3uWD1otI6ApAsi4l5J/6rW8Yj4Vt1a1iBKpVL09vbmbsbvvPPPP/81sfvvv7/wdpjVsmTJEgBuvPHGzC05ckjaGBGlscqNNRz3p8C9QK219wM44pOQHR7t7e309/cf3J85c+Yopc2sWYyahCJiqaQpwJqIWFVQm+wINHPmzGFJqL29PWNrzKxRjHlPKCIOULn3YjZpDz300Kj7ZtacxjtFu0fSJyTNlHTC0KeuLbMjysh7j56UaGYw/ina/5bKPaAPj4i/5fA2x45UU6ZMObhiwtC+mdl4/yU4ncpioY8CjwA3UVlQ1Gxc5s2bN+q+mTWn8SahLuDtwDIqCejtvLoQqNmYOjs7R903s+Y03uG4t0bEmVX790l6tB4NMjOz5jHentDDkuYO7Ug6D/hBfZpkR6KvfvWrw/aXL1+eqSVm1kjGm4TOA34oaXtakPQB4E8lbZb0WN1aZ0eM9euHv/nDr3IwMxj/cFxHXVthR7yhFbQPtW9mzWlcSSginqx3Q8zMrPn4YQ0zM8vGScgKccIJJ4y6b2bNyUnICrFnz55R982sOTkJWSGql+yptW9mzaluSUjSMZJ+JOlRSVskfSbFT5P0kKRtkr4h6agUPzrt96Xjs6qudV2K/1TSxVXxjhTrk3RtVXzCdZiZWfHq2RN6CbggrbRwFtCRHnj9EnBDRMwBdgNXpPJXALsjYjZwQyqHpNOBy6msVdcBfEXSVElTqaxndwmVte0WpbJMtA4zM8ujbkkoKl5Iu69LnwAuAO5M8S7g0rS9kFfXo7sTuFCSUvyOiHgpIn4O9AHnpk9fRDwRES8DdwAL0zkTrcPMzDKo6z2h1GN5BHgO6AF+BvwyIvalIv1AW9puA54GSMf3AC3V8RHnHCreMok6zMwsg7omoYjYHxFnAe1Uei5vr1UsfdfqkcRhjI9WxzCSOiX1SuodGBiocYqZmR0OhcyOi4hfAvcDc4HjJA2t1NAO7Ezb/cBMgHT8TcCu6viIcw4V/8Uk6hjZ3uURUYqIUmtr6+R+aDMzG1M9Z8e1SjoubU8H5gFbgfuAD6RiZeCutN2d9knH743KO6C7gcvTzLbTgDnAj4ANwJw0E+4oKpMXutM5E63DzMwyGO8CppNxCtCVZrFNAVZFxHckPQ7cIenzwMPALan8LcD/lNRHpXdyOUBEbJG0Cngc2AdcFRH7ASR9BFgHTAVWRMSWdK1PTqQOMzPLo25JKCIeA95VI/4ElftDI+O/AS47xLW+AHyhRnw1sPpw1GFmZsXziglmZpaNk5CZmWXjJGRmZtk4CZmZWTZOQmZmlo2TkJmZZeMkZGZm2TgJmZlZNk5CZmaWjZOQmZll4yRkZmbZOAmZmVk2TkJmZpaNk5CZmWVTz/cJmVkDu+mmm+jr68vdjIYw9OewZMmSzC1pDLNnz+bqq68upC4nIbMm1dfXx7YtD/N7b9ifuynZHfVKZVDopSd7M7ckv6demFpofU5CZk3s996wn0+d/avczbAG8sVNxxZaX93uCUmaKek+SVslbZG0JMVPkNQjaVv6Pj7FJWmZpD5Jj0k6u+pa5VR+m6RyVfwcSZvTOcskabJ1mJlZ8eo5MWEf8B8i4u3AXOAqSacD1wL3RMQc4J60D3AJMCd9OoGboZJQgKXAeVRe2b10KKmkMp1V53Wk+ITqMDOzPOqWhCLimYjYlLafB7YCbcBCoCsV6wIuTdsLgZVR8SBwnKRTgIuBnojYFRG7gR6gIx07NiIeiIgAVo641kTqMDOzDAqZoi1pFvAu4CHg5Ih4BiqJCjgpFWsDnq46rT/FRov314gziTrMzCyDuichSW8A/jfw0YgY7Q6oasRiEvFRmzOecyR1SuqV1DswMDDGJc3MbLLqmoQkvY5KAvp6RHwrhZ8dGgJL38+leD8ws+r0dmDnGPH2GvHJ1DFMRCyPiFJElFpbW8f/A5uZ2YTUc3acgFuArRHx5apD3cDQDLcycFdVfHGawTYX2JOG0tYB8yUdnyYkzAfWpWPPS5qb6lo84loTqcPMzDKo53NC7wH+DbBZ0iMp9ingemCVpCuAp4DL0rHVwPuAPmAv8CGAiNgl6XPAhlTusxGxK21fCdwKTAfWpA8TrcPMzPKoWxKKiO9T+x4MwIU1ygdw1SGutQJYUSPeC7yjRnxwonWYmVnxvICpmZll4yRkZmbZOAmZmVk2TkJmZpaNk5CZmWXjJGRmZtk4CZmZWTZOQmZmlo2TkJmZZeMkZGZm2dRz7Tgza2A7duzg189P5Yubjs3dFGsgTz4/ldfv2FFYfe4JmZlZNu4JmTWptrY2Xtr3DJ86e7R3TVqz+eKmYzm6rbgXTrsnZGZm2TgJmZlZNk5CZmaWjZOQmZllU7eJCZJWAH8GPBcR70ixE4BvALOA7cCfR8RuSQJupPLq7b3AX0bEpnROGfh0uuznI6Irxc/h1Vd7rwaWRERMpo4j3U033URfX1/uZrzGkiVLstQ7e/Zsrr766ix1m9lw9ewJ3Qp0jIhdC9wTEXOAe9I+wCXAnPTpBG6Gg0lrKXAecC6wVNLx6ZybU9mh8zomU4eZmeVTt55QRHxX0qwR4YXA+Wm7C7gf+GSKr4yIAB6UdJykU1LZnojYBSCpB+iQdD9wbEQ8kOIrgUuBNROtIyKeOZw/dyNqhP/1n3/++a+J3XjjjcU3xMwaStH3hE4e+kc/fZ+U4m3A01Xl+lNstHh/jfhk6rACHHPMMcP2p0+fnqklZtZIGmVigmrEYhLxydTx2oJSp6ReSb0DAwNjXNbGY+3atcP216xZk6klZtZIik5Cz6ZhNtL3cyneD8ysKtcO7Bwj3l4jPpk6XiMilkdEKSJKra2tE/oBbWzuBZnZkKKTUDdQTttl4K6q+GJVzAX2pKG0dcB8ScenCQnzgXXp2POS5qZZb4tHXGsidVhBzjzzTM4880z3gszsoHpO0b6dygSBEyX1U5nldj2wStIVwFPAZan4aipTp/uoTJ/+EEBE7JL0OWBDKvfZoUkKwJW8OkV7Tfow0TrMzCyfes6OW3SIQxfWKBvAVYe4zgpgRY14L/COGvHBidZhZmZ5eBVtsyb21At+nxDAs3srdyZOnnEgc0vye+qFqcwpsD4nIbMmNXv27NxNaBgvpxVFjn6z/0zmUOzvhpOQWZNqhIeYG8XQElJ+gLp4TkJ11qjrtuUw9OeQa824RuM17MychOqur6+PR368lf0zTsjdlOymvFx5NnjjE89mbkl+U/fuGruQWRNwEirA/hkn8OLb3pe7GdZApv9kde4mmDWERlm2x8zMmpCTkJmZZePhuDrbsWMHU/fu8fCLDTN17yA7duzL3Qyz7NwTMjOzbNwTqrO2tjb+6aVpnphgw0z/yWra2k7O3Qyz7NwTMjOzbNwTKsDUvbt8TwiY8ptfAXDgGK9VVnlOyD0haIwHuhvpQepme4jZSajOvD7Xq/r6ngdg9lv8jy+c7N+NBuIXLeajyhsO7FBKpVL09vbmbsYRwetzmTUPSRsjojRWOd8TMjOzbJyEzMwsm6ZLQpI6JP1UUp+ka3O3x8ysmTXVxARJU4G/By4C+oENkroj4vG8LauvRph9BI0zA6nZZh+ZNbJm6wmdC/RFxBMR8TJwB7Awc5uaxvTp0z0LycyGaaqeENAGPF213w+cl6kthfH/+s2sUTVbT0g1Yq+Zoy6pU1KvpN6BgYECmmVm1pyaLQn1AzOr9tuBnSMLRcTyiChFRKm1tbWwxpmZNZtmS0IbgDmSTpN0FHA50J25TWZmTaup7glFxD5JHwHWAVOBFRGxJXOzzMyaVlMlIYCIWA14NVEzswbQbMNxZmbWQJyEzMwsGychMzPLxq9yGIOkAeDJ3O04gpwI/CJ3I8xq8O/m4fXmiBjzGRcnISuUpN7xvGPErGj+3czDw3FmZpaNk5CZmWXjJGRFW567AWaH4N/NDHxPyMzMsnFPyMzMsnESskL4terWqCStkPScpB/nbkszchKyuqt6rfolwOnAIkmn522V2UG3Ah25G9GsnISsCH6tujWsiPgusCt3O5qVk5AVodZr1dsytcXMGoiTkBVhXK9VN7Pm4yRkRRjXa9XNrPk4CVkR/Fp1M6vJScjqLiL2AUOvVd8KrPJr1a1RSLodeAB4q6R+SVfkblMz8YoJZmaWjXtCZmaWjZOQmZll4yRkZmbZOAmZmVk2TkJmZpaNk5BZA5H0UUkzJnHe2yQ9IulhSb9/GNpxqReZtSI4CZk1lo8CE0pCaZXyS4G7IuJdEfGzqmOSNJm/55dSWfHcrK6chMzqQNLfSrombd8g6d60faGk/yXpZkm9krZI+kw6dg1wKnCfpPtSbL6kByRtkvRNSW9I8e2S/pOk7wMfpJK8/p2k+yTNkrRV0leATcBMSYskbZb0Y0lfqmrnC5K+IOlRSQ9KOlnSHwELgL9LvavfumdldihOQmb18V3gn6ftEvAGSa8D/hj4HvAfI6IEvBP4U0nvjIhlVNbUe29EvFfSicCngXkRcTbQC3y8qo7fRMQfR8RtwD8AN0TEe9OxtwIrI+JdwCvAl4ALgLOAd0u6NJV7PfBgRJyZ2vxXEfFDKssq/U1EnFXdszI73JyEzOpjI3COpDcCL1FZFqZEJTF9D/hzSZuAh4E/pPbQ19wU/4GkR4Ay8Oaq498Ypf4nI+LBtP1u4P6IGEhLKH0d+JN07GXgO1VtnjWRH9LstzUtdwPMjkQR8Yqk7cCHgB8CjwHvBX4feBH4BPDuiNgt6VbgmBqXEdATEYsOUc2vR2lC9bFar9IY8kq8unbXfvxvghXMPSGz+vkulWTzXSq9n78GHgGOpZIk9kg6mcprz4c8D7wxbT8IvEfSbABJMyT9wSTa8RCVIb8T0ySGRcA/jnFOdTvM6sZJyKx+vgecAjwQEc8CvwG+FxGPUhmG2wKsAH5Qdc5yYI2k+yJiAPhL4HZJj1FJSm+baCMi4hngOuA+4FFgU0TcNcZpdwB/c7imfJsdilfRNjOzbNwTMjOzbJyEzMwsGychMzPLxknIzMyycRIyM7NsnITMzCwbJyEzM8vGScjMzLL5/7PzJ25ACb4KAAAAAElFTkSuQmCC\n"
                    },
                    "metadata": {
                        "needs_background": "light"
                    }
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 5\nUse 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."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "sns.regplot (x=\"sqft_above\", y =\"price\", data=df)\nplt.ylim(0,)",
            "execution_count": 23,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 23,
                    "data": {
                        "text/plain": "(0, 8086228.609231744)"
                    },
                    "metadata": {}
                },
                {
                    "output_type": "display_data",
                    "data": {
                        "text/plain": "<Figure size 432x288 with 1 Axes>",
                        "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAELCAYAAABwLzlKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvXucHNd13/k9Vf2aRw9mAMzgDYGQIIGkniQiUTE/NCw7EilnRW2WSgRnI1qRF1xHiu14pZWUjywlVOwPGXvjiI6tEGs5JtexaBprh9ysKIUWDcHeEJIAyjRFAiJAEAIGr5kB5tE9/a46+0dV9fTM9Mz0DKan53G++Mynu2/fqlvV3ahT59zfPUdUFcMwDMNoBU6rD8AwDMNYu5gRMgzDMFqGGSHDMAyjZZgRMgzDMFqGGSHDMAyjZZgRMgzDMFqGGSHDMAyjZZgRMgzDMFpGU42QiPwLEXlZRH4oIl8XkZSI3CQi3xWR0yLyJyKSCPsmw9dnwvd31ezn82H7j0TkAzXtd4dtZ0TkczXt8x7DMAzDWHqkWRkTRGQb8NfALaqaF5EngW8AHwT+TFWfEJH/CLyoql8VkX8GvF1V/1cR+SjwP6rqPxKRW4CvA+8GtgJ/Abw5HOZV4O8B/cD3gQOq+ko4VsNjzHYeGzdu1F27di3qZ2MYhrHaOXHixJCq9s7VL9bk44gBbSJSBtqBy8D7gJ8L338M+FfAV4F7w+cAh4H/ICIStj+hqkXgdRE5Q2CQAM6o6lkAEXkCuFdETs53DJ3FEu/atYvjx48v8PQNwzDWJiLy40b6NS0cp6oXgd8CzhMYn1HgBDCiqpWwWz+wLXy+DbgQblsJ+2+obZ+yzUztGxYwhmEYhtECmmaERKSHwPO4iSCM1gHcU6dr5IXIDO8tVvtsY0xCRA6KyHEROT44OFhnE8MwDGMxaKYw4WeA11V1UFXLwJ8BfxfoFpEoDLgduBQ+7wd2AITvrwOu17ZP2Wam9qEFjDEJVT2kqvtUdV9v75whTcMwDGOBNNMInQfuEJH2cG7np4FXgL8E7gv73A88FT5/OnxN+P5z4VzN08BHQ2XbTcAe4HsEQoQ9oRIuAXwUeDrcZr5jGIZhGC2gacIEVf2uiBwGXgAqwA+AQ8D/CzwhIv8mbPtauMnXgP8rFB5cJzAqqOrLodrtlXA/n1RVD0BEPgV8C3CBP1DVl8N9fXY+YxiGYRitoWkS7dXCvn371NRxhmEY80NETqjqvrn6NVuibRiGYdwgR04N8OjRs1wYzrGjp50H7trN/r19rT6sRcHS9hiGYSxjjpwa4ItPv8xApkB3W5yBTIEvPv0yR04NtPrQFgUzQoZhGMuYR4+eJe4K7YkYIsFj3BUePXq21Ye2KJgRMgzDWMZcGM7RFncntbXFXfqHcy06osXFjJBhGMYyZkdPO/myN6ktX/bY3tPeoiNaXMwIGYZhLGMeuGs3ZU/JlSqoBo9lT3ngrt2tPrRFwYyQYRjGMmb/3j4e/NCt9KVTjObL9KVTPPihW1eNOs4k2oZhGMuc/Xv7Vo3RmYp5QoZhGEbLMCNkGIZhtAwzQoZhGEbLMCNkGIZhtAwzQoZhGEbLMCNkGIZhtAwzQoZhGEbLMCNkGIZhtAwzQoZhGEbLsIwJhmEYDbKai8u1iqZ5QiLyFhH5m5q/MRH5FRFZLyLPisjp8LEn7C8i8oiInBGRvxWR22r2dX/Y/7SI3F/TfruIvBRu84iISNg+7zEMwzBmY7UXl2sVTTNCqvojVX2nqr4TuB3IAX8OfA74tqruAb4dvga4B9gT/h0EvgqBQQG+BLwHeDfwpciohH0O1mx3d9g+rzEMwzDmYrUXl2sVSzUn9NPAa6r6Y+Be4LGw/THgw+Hze4HHNeAY0C0iW4APAM+q6nVVHQaeBe4O3+tS1edVVYHHp+xrPmMYhmHMymovLtcqlsoIfRT4evh8k6peBggfo4DqNuBCzTb9Ydts7f112hcyhmEYxqys9uJyraLpRkhEEsCHgD+dq2udNl1A+0LGmNxJ5KCIHBeR44ODg3Ps0jCMtcByLy535NQABw4d486Hn+PAoWMrZq5qKTyhe4AXVPVq+PpqFAILH6NPqh/YUbPdduDSHO3b67QvZIxJqOohVd2nqvt6e3vncaqGYaxWlnNxuZUsmlgKifYBJkJxAE8D9wMPhY9P1bR/SkSeIBAhjKrqZRH5FvAbNWKE9wOfV9XrIpIRkTuA7wIfA35nIWMs+hkbhrEqWa7F5WpFEwDtiRi5UoVHj55dlsdbS1ONkIi0A38PeKCm+SHgSRH5BHAe+EjY/g3gg8AZAiXdxwFCY/Nl4PthvwdV9Xr4/BeBPwTagGfCv3mPYRiGsZK5MJyjuy0+qW2liCaaaoRUNQdsmNJ2jUAtN7WvAp+cYT9/APxBnfbjwFvrtM97DMMwjJXKjp52BjKFqicEK0c0YWl7DMMwVjjLXTQxG2aEDMMwVjjLWTQxF5Y7zjAMYxWwXEUTc2GekGEYhtEyzAgZhmEYLcOMkGEYhtEyzAgZhmEYLcOMkGEYhtEyzAgZhmEYLcOMkGEYhtEyzAgZhmEYLcOMkGEYhtEyzAgZhmEYLcOMkGEYhtEyzAgZhmEYLcOMkGEYhtEyzAgZhmEYLcOMkGEYhtEymmqERKRbRA6LyCkROSki7xWR9SLyrIicDh97wr4iIo+IyBkR+VsRua1mP/eH/U+LyP017beLyEvhNo+IiITt8x7DMAzDWHqa7Ql9Bfimqu4F3gGcBD4HfFtV9wDfDl8D3APsCf8OAl+FwKAAXwLeA7wb+FJkVMI+B2u2uztsn9cYhmEYRmtomhESkS7gLuBrAKpaUtUR4F7gsbDbY8CHw+f3Ao9rwDGgW0S2AB8AnlXV66o6DDwL3B2+16Wqz6uqAo9P2dd8xjAMwzBaQDM9od3AIPCfROQHIvL7ItIBbFLVywDhY1SPdhtwoWb7/rBttvb+Ou0sYAzDMAyjBTTTCMWA24Cvquq7gHEmwmL1kDptuoD22WhoGxE5KCLHReT44ODgHLs0DMMwFkozjVA/0K+q3w1fHyYwSlejEFj4OFDTf0fN9tuBS3O0b6/TzgLGmISqHlLVfaq6r7e3t+ETNgzDMOZH04yQql4BLojIW8KmnwZeAZ4GIoXb/cBT4fOngY+FCrY7gNEwlPYt4P0i0hMKEt4PfCt8LyMid4SquI9N2dd8xjAMwzBaQKzJ+//nwH8WkQRwFvg4geF7UkQ+AZwHPhL2/QbwQeAMkAv7oqrXReTLwPfDfg+q6vXw+S8Cfwi0Ac+EfwAPzWcMwzAMozVIICwzZmLfvn16/PjxVh+GYRjGikJETqjqvrn6WcYEwzAMo2WYETIMwzBahhkhwzAMo2U0W5hgGIYBwJFTAzx69CwXhnPs6Gnngbt2s39v39wbGqsa84QMw2g6R04N8MWnX2YgU6C7Lc5ApsAXn36ZI6cG5t7YWNWYETIMo+k8evQscVdoT8QQCR7jrvDo0bOtPjSjxZgRMgyj6VwYztEWdye1tcVd+odzLToiY7lgRsgwjKazo6edfNmb1JYve2zvaW/RERnLBTNChmE0nQfu2k3ZU3KlCqrBY9lTHrhrd6sPzWgxpo4zDKPp7N/bx4MEc0P9wzm2L1AdZwq71YcZIcNYBayEi/P+vX03dEyRwi7uyiSF3YPhvo2ViYXjDGOFs1bkz6awW52YETKMFc5auTibwm51YkbIMFY4a+XibAq71YkZIcNY4ayVi7Mp7FYnZoQMY4WznC7OR04NcODQMe58+DkOHDq2qPNS+/f28eCHbqUvnWI0X6YvneLBD91qooQVjhW1mwMramesBCJ13I3InxfjGCL1WlvcJV/2KHtqhmKN0mhRO5NoG8Yq4Eblz4tBrUACoD0RI1eq8OjRsy0/NmP50tRwnIicE5GXRORvROR42LZeRJ4VkdPhY0/YLiLyiIicEZG/FZHbavZzf9j/tIjcX9N+e7j/M+G2stAxDMO4MdaKQMJYXJZiTuinVPWdNW7Z54Bvq+oe4Nvha4B7gD3h30HgqxAYFOBLwHuAdwNfioxK2OdgzXZ3L2QMwzBunLUikDAWl1YIE+4FHgufPwZ8uKb9cQ04BnSLyBbgA8CzqnpdVYeBZ4G7w/e6VPV5DSa2Hp+yr/mMYRjGDbKcBBLGyqHZRkiB/yYiJ0TkYNi2SVUvA4SPUbB4G3ChZtv+sG229v467QsZYxIiclBEjovI8cHBwXmcrmGsXUy9ZiyEZgsTfkJVL4lIH/CsiJyapa/UadMFtM9GQ9uo6iHgEATquDn2aRhGyHIQSBgri6Z6Qqp6KXwcAP6cYE7nahQCCx+jhQT9wI6azbcDl+Zo316nnQWMYRiGYbSAphkhEekQkXT0HHg/8EPgaSBSuN0PPBU+fxr4WKhguwMYDUNp3wLeLyI9oSDh/cC3wvcyInJHqIr72JR9zWcMwzAMowU0Mxy3CfjzUDUdA/5YVb8pIt8HnhSRTwDngY+E/b8BfBA4A+SAjwOo6nUR+TLw/bDfg6p6PXz+i8AfAm3AM+EfwEPzGcMwDMNoDZYxYQ4sY4JhGMb8sYwJhrEKWQnF6wxjPlgCU8NYIayV4nXG2sKMkGGsENZK8TpjbWFGyDBWCJabzViNmBEyjBWC5WYzViMmTDCMFcIDd+3mi0+/TK5UmVSvZyXnZouEFq9eHaPsKYmYw56+tAku1hDmCRnGCmG15WaLhBavD2UZK1TIlz1Gc2XOXcua4GINYZ6QYawgVlNutkhocS1bwUFwHMFXZSxfYfO6mBXDWyOYJ2QYRkuIhBYlz0fC1MIiUPJ8E1ysIcwIGYbREiKhRcJ1iBK3qELCdUxwsYZo2AiJyBtE5GfC521RclLDMIyFEBXBS6di+CgV38f3la622IoXXBiN09CckIj8LwTlsNcDbyQogfAfgZ9u3qEZhrGa2b+3jwcJ5oYq3hilUB23a0OnqePWEI0KEz5JUAvouwCqejosVGcYhrFgVpPQwlgYjRqhoqqWwrIMiEiMuauYGoaxArEkqcZS0uic0HdE5F8CbSLy94A/Bf6f5h2WYRitwJKkGktNo0boc8Ag8BLwAEFxuC8066AMw2gNliTVWGoaDce1AX+gqv8ngIi4YZsJ+Q1jFXFhOEd3W3xSm63ZMZpJo57QtwmMTkQb8BeLfziGYbQSS5JqLDWNGqGUqmajF+Hzhn6VIuKKyA9E5L+Gr28Ske+KyGkR+RMRSYTtyfD1mfD9XTX7+HzY/iMR+UBN+91h2xkR+VxN+7zHMAxjYu1OrlRBNXisXbNz5NQABw4d486Hn+PAoWM2V2TcMI0aoXERuS16ISK3A/kGt/1l4GTN64eB31bVPcAw8Imw/RPAsKq+CfjtsB8icgvwUeBW4G7g90LD5gK/C9wD3AIcCPvOewzDMAJmS5JqogWjGTQ6J/QrwJ+KyKXw9RbgH821kYhsB34W+HXgVyXQeL8P+Lmwy2PAvwK+CtwbPgc4DPyHsP+9wBOqWgReF5EzBGuWAM6o6tlwrCeAe0Xk5HzHUFWTmxtGyExrd2pFCwDtiRi5UsUSjRo3RENGSFW/LyJ7gbcAApxS1XIDm/574H8HohQ/G4ARVa2Er/uBbeHzbcCFcLyKiIyG/bcBx2r2WbvNhSnt71ngGEO1By0iBwkyRLBz584GTtMwVj8mWjCawaxGSETep6rPicg/mPLWHhFBVf9slm3/PjCgqidEZH/UXKerzvHeTO31Qomz9Z9r/IkG1UPAIYB9+/aZl2QYBKKFc9eyjOUrlDyfhOvQ1RZj14bOVh+asYKZyxP6SeA54H+o854CMxoh4CeAD4nIB4EU0EXgGXWLSCz0VLYDUYivH9gB9IcZGdYB12vaI2q3qdc+tIAxDMOYg/fuXs/3zl3HEXDCkgsDmRIH/s76Vh+asYKZ1Qip6pdExAGeUdUn57NjVf088HmA0BP6tKr+YxH5U+A+4AngfuCpcJOnw9fPh+8/p6oqIk8Dfywi/w7YCuwBvkfg1ewRkZuAiwTihZ8Lt/nL+Ywxn/MyjBthMVPiLHV6nefPXqe3M0GmMOEJpVMxnj97nV9q2qjGamdOdZyq+sCnFnHMzxKIFM4QzMd8LWz/GrAhbP9VgiwNqOrLwJPAK8A3gU+qqhd6OZ8CvkWgvnsy7DvvMQxjKVhMdVkrlGoXhnNs7Eyyu7eTvZu72N3bycbOpM0JGTeENOIIiMivEUiy/wQYj9pVddWHsvbt26fHjx9v9WEYi0CrE3MeOHSMgUyhqi4DyJUq9KVTfP3gHS3b13Ie01i5iMgJVd03V79G1wn9U+CfAd8Bjtf8GcaKYDmscYnKWdeyUHXZYu6rUeZayGoYC6HRdUK3EBihOwkECX9FUNTOMFYEi73GJfKqXr06RjksxranLz2rd7Wjp32aJ7GQlDhHTg0wli9zeTRPKubSm06STsVvOL3OXJ5ibRG6/uEc263Mg7EINBqOexIYA/5z2HQA6FbVf9jEY1sWWDhudXDnw8/R3RYnqomVKZQZGCtQ9JR371o/r4tp5FWVKh7XxktBo8LGdIK461YzDMy0XdwV2uIu+bJH2dMZ+y/m2I3u90aOzTBqaTQc16gn9BZVfUfN678UkRcXdmiGsfTUeiGZQplLIwUUJRVzGMgU+PThF+ntTJIpVuacL4q8qmvZCg6C4wi+KmP5CpvXxWb0rhbDk4jGXteWIhlzGcoWKVZ8xosej3z07Qs2GJYNwWgVjRqhH4jIHap6DEBE3gP8f807LMNYXB64azdffPplcqUKA2OBAQLY2Jmk4ikjuTLZQoU39XVW54sehLoX4ChzQMnzcavVhoN1M3PNy9xoOevarAVdbXG62uKoKqP58qLtN8KyIRhLQaPChPcA/11EzonIOYJ1Nj8pIi+JyN827egMY5GoTcxZ9JSE67B1XRtdbXGGskUcAU+1oUJuUbmDhOsQRbNVIeE6TS970KxSC1bCwWgVjRqhu4GbCDIo/GT4/IPA36d+NgXDWHbs39vH1w/ewbt3rWfzuhRd4Z1/yfOBwIhEzOYFRCqxdCqGj1LxfXxf6WqLNV0t1iyFminfjFbRaALTHzf7QAyjWUxVfb1393oOv3CRXKlCW9zFdYSKp/Smk9VtIi9gJsVYNLdT8cYoheq4XRs6m64Wa4ZCLTrH8WK5YaWfYSwWDanj1jKmjlueNLrwdCbV1323beP5s9fpH87RkXC5Nl6iqy0+rc/hFy6uasWYqeKMZrHY6jjDWDbUXji72+K8PpTlgT86QToVm3YHP5Pq6/mz1yet8o+MWq13sRYUY2vhHI3ljRkhoyXcSAqd2gvnWL5cXS+TK1amKdsaVX3VU6194akfNl0x1upUQs1UxbX63FrNWj//RmlUmGAYi8aNptCpTVkzlC3iILiOUPZ1mrLtRlRfzVKMHTk1wIFDx7j9y/+NB/7oBOeuZVuWSqiZ59jqNEmtZK2f/3wwI2QsObWeTCOS6KnUXjhLno/IhEQaJt/J34jqqxmKsdqLU6Hs46tyLVsmW6zM+3NYDJqlirvR73ils9bPfz5YOM5YcuYTAqoX0qhdeJpwncAQIVV1W+2d/FxqstlCJs1QotVenKLFrgoMZoqkU/E5P4f55Kqb63OMQpDNyAe31he/rvXznw9mhIwbZr6x70YTeU4VIFTnez50Kw9+6FYePXqW0VyJiq+s74jTmYzVvZOfKUvBjPuHSYZoMeP4tRenhOtQ8RRxJtYqzfY5lCoeY4VK0K/kce5adtbMDo2e52KfIyxestaVylo///lg4TjjhlhI7LvRENBsIY1o4enxX3s/j/7Pt7NrQyej+TJ96VTD8uJWhExqQ4kbO5P4KOWKT8VTXrk8Sv9wnvfunlwuOzrOTCHIVRdzHBxHGMtXGjreVpznWl/8utbPfz6YJ2TcEI1IfOt5SpEnE4WA3rt7PY8ePcsXnvphtc+NKNsaIdp/plBmMFOk5PnEHWE0X55z24Uqn2pDielUjI68y0i+giuQirmkUzEOv3CRt2/vru5vvrnqph7bq1fH2LKurfr+WL7MULbIuWs5Dhw61hTV1lov+7DWz38+NM0IiUgKOAokw3EOq+qXROQm4AlgPfAC8E9UtSQiSeBx4HbgGvCPVPVcuK/PA58APOCXVPVbYfvdwFcAF/h9VX0obJ/3GMbCmMtQzBZSi9bpzNQnnYyRL3tNC2ns6Gnn3LUs17JlRKgq7DKFCkdODcx4wWgkjDcTUy9OvsKWdUk2dqaqfaYa8Si0Uw3fycy56uodW7boMZQt0ptOMZYvc2k0D0DSlXkd+3xpRphvJbHWz79RmhmOKwLvC0tAvBO4W0TuAB4GfltV9wDDBMaF8HFYVd8E/HbYDxG5BfgocCtBDrvfExFXRFzgd4F7CIruHQj7Mt8xjIUzl8S3kVDQTH1UtWkhjSOnBugfznF5NPCAimFIDKCnPT5jqOqRv3iVTzx+nPPXc7w+NM7AWGHG8FYkxb7z4ec4cOhYNUQZhRL/6rPvo6stzoaO5KTtpno3c+Wqe+/u9dVxfumJH1CqeJM+y/UdcYZzZXKlCkPZIgCC0NeVMtWW0XKaZoQ0IBu+jId/CrwPOBy2PwZ8OHx+b/ia8P2flqAC2b3AE6paVNXXgTPAu8O/M6p6VlVLBJ7PveE28x3DWCBzxb4bKUM9U5/xklfNfD3f+Z7ZOHJqgM8cfpHLo4VJ7RVf6WmLs7EzWVfF9MhfvMpXnjuD5wfGylcYyJa4Opqfdk6NzpU1sk4nygB+08ZO1qVitMVd1rXH2bWhs5paKBpnvFTh2niJsZqQ4oaOJOlUjL50ikIlCDlu7U6RTgUebCtUWzMZaGPt0dQ5odBbOQG8icBreQ0YUdVK2KUf2BY+3wZcAFDVioiMAhvC9mM1u63d5sKU9veE28x3jKEbPtk1ylyx70ZUQrP1aSSkMdP8zEztjx49S6ZQwQ2L0akGd0ciMF7yZgz5/f5fv44joBIYoIiBbAnHEW7a2Fltm2mu7OFvnpo1mWq+7DGWLxN3hDsffm6anHrqOf/ukdcQYPO6FCJCKuZS8nyGssVqlvB82WNPX5qvH7yDA4eOtVy1dSPhTGP10VQjpKoe8E4R6Qb+HLi5XrfwsZ5HorO01/PiZus/2xiTEJGDwEGAnTt31tnEqGU2Q1E7EV+bILM2pNZIn5mILmhlz2M0V+byaJ4Xzg/zwbdu4sT50boXugvDOSq+T8x1iDkO5VAerQqFyvSxowv+WKGCUP9HNJgt8XPvXl/t/8L5YTzfJxlz6U0nSafiVDyfc9dy7NrQjivwg/PDfO/cNbakkyQSCUbzZTqTMRQo+zrpuO/rH+H5s9e5MJyjsybhqq+KAJdGCvS0e3h+EMIsex5j+RIx12E0XybhOtz58HOkk7Gq8GK+n/VCqHcjYPnqjFqWRKKtqiPAEeAOoFtEIuO3HbgUPu8HdgCE768Drte2T9lmpvahBYwx9XgPqeo+Vd3X29u7sJM2gMnF5GYKqTXSB+qHcB49epay53EtW8ZTiLsOvir/5cXLlD2v7lzUjp52Yk5QkM51hLjrVI1LRyI2aezasJojwR2LX+c8+9IJnj97vdpfBBwRKr5yaaRAplDm6liRuCt4vnJ5tIgCrghDuTLjJY8v3/tWutsTrGuLTzrusufxu0deq4bczl3PMZwr4/lBcT4h8OgGMsE+Yw44Av0jBYplDyFQ0kUKO4GqCnCxQpz1mCkk+erVsTlDtMbaoZnquF6grKojItIG/AyBEOAvgfsI5nDuB54KN3k6fP18+P5zqqoi8jTwxyLy74CtwB7gewTXjD2hEu4igXjh58Jt5jVGsz4DI6CRkNpcfWYK4eRKFXLFSvWiD+AKlFUZzZUnqc6iC92X730rnzn8IsO5MirB1+86Qnd7nN+67x2TjqP2rn1jR4KBbKn6nhAYpU3pJBs6gnmkqP+mdIpLo3mC3StXRguUfZ/t3W0MZorV41XA87VqIOupDUdDgxN5Dp6vOBJkWdjYmeTSaB7PDwuWK4gIO7rbiLnCYKZIbzo5yesA6OlI8s1/cQfNZCaPp+xpU1WPxsqimZ7QFuAvw/Lf3weeVdX/CnwW+FUROUMwH/O1sP/XgA1h+68CnwNQ1ZeBJ4FXgG8Cn1RVL5zz+RTwLeAk8GTYl/mOYSx/ZlLQlSo+xTB/XIRq4AkUvck+S+0802/e9w729HUiIogIb+rt4GN3vIFHj56d5GnViiY2rWujrzNR3Z9IYID6ulLVfUf9u9ribF0XGAIlMFZ7ejuJRWmGZOJYE65TNZD1hApFzycZm/ivGuXIK3l+dRwIjGLMlWrZ8kjc0SqvYybBSSLm2EJOo0rTPCFV/VvgXXXazxIo26a2F4CPzLCvXwd+vU77N4BvLMYYxvLmwnAOV+DsYJaS55NwHTZ2Jki4QtFx8HzFDdfP+ChdqRj5sj/jPFO9if56nlZnwp10175pXRuOIwznymzvaaMt7k66iD569Gx14r+rLU5XW5xcqUJfOlWd+xKgUJ4wkCUPzgxm2bW+ve78WMxxSKcm/qtu7ExycSRPzAlk7DFXSMQc1nfEJ3l++bJHx5Tjj9qXwuuYSXAS5byzhZwGWMYEY4XQmXA5MziOK4IrQTnuiyOFwIN57xZ+98hrVHwlGXOIuw75so8jcO5aDlWlMxnjF+68acYL3UyhIxFhNF/i4kiesqdoqKbb1p0iEU76dyRcEq7DF5764awT//v39nFf/wi/85dnJo2tQLniV+siTc0mce87tk5S0MXcIHTY25lkNF+u2yca9xfuvKlue6OijxuphzOb4GS1LeS02kEBpYpPoeJV19w1ghkhY0VQXc5VK08L5z9+6WfezNu3d/Po0bOcvjpGphh4AJlCBUcAEdoT7rR0OLXMlPnhymgeEcH3tbo+yHWCNUXjJY+P3L69WgK8uy1OPhQCRBP/26dIxn//r19HlWrWAwnPKe46dLUFi2S/fvCOaccYnV9kmH7tZ2+Zs0807kzts7EYMuq1krpmrUrOfV8pVnyKFY9COXiM/o+UvXrynfqYETJWBJlihW3dKYaypWo4bnNXkmwxWA7bvPp0AAAgAElEQVQW3VlH62CujBZwHMEJDUimUGHzutg0GXB0BzuYKTKUKbJ53cQiznzZo+QpvekEuZKHiE7aXzoFX3nuDKCkaqTYMH3iP7pQjZcqxByhEq5iizmCGyrm6i14nc/d9UzexUK8jsWSUa82j6cea0VyXg4zixTKHoXQq1VVfFXOX89x8tIYJ69keOXyGOeGxhverxkhoyksdngiml/Y3TuxIDSaa6llUrJPp36yz+jYTg9kyBQq9LTH2dyV5OJIgf7hPNu6lZjrVOv2tMXdafsrVHxK4yU8X0nGJqTYW7uhMxmbNvEfXahSMZeKrzii+AqeKo7KtDxwrb67tno4jbMaPyvV0MspT3g6FT/wboZzJU5eHuPk5QynLo9x6kqG8ZI3xx5nxoyQsejcyAV0JuMVzS8MZQuM5soUPZ+Y43DvO7ZO2n5Ssk9/erLPjoRbPbZcsULF87k6Fkim466D7/tcGSty04Z2Eq7D1bEC2UI2KD5XE0bT8IUjgAqOI/gog5kiriPTJv5PD2TIFSuUPK1KrCHIvOD7SldHfFIeuBfOD0/KhNDo3fViGX+rh9M4q+Gz8nylUPaqnk6x4qOqlCo+rw1meSU0Oicvj01LdxXRlYpx85Yubt6S5q3b1vGRBjNzmhEyFp2FhifqGa/PHH6RDR0JsiUPfJ+hbAnVoOxB3BV+98hrPH7sx1XFVWSsutpiDGVK+OE6oHQquMgnXKe6YDRXo1Cr6qhFSLiQK/vEXaEr5TKYncjD5go40apVhY0dCYbz5WAFqyiFil8340KmUKHi+9UUQdG8retINQ9clMIn7sqkTAhbu5m16upsn99CjX+9dEImo67PjWT8aBXFSo3BKfuUvcDoXBotcOryGK+EBue1wSzlOiKDmCO8sa+TW0Kjc/PmLrZ2p3Adh2TcmSbNnw0zQkZDzOcOe6F1eqYaL89XhnNlMsUKb+rt5MxgFlXY3tOGKtWSBLlihdeHsnzise/jhxd5RyAVd4k5QsIN8rq9d/f6QEXn+ZNyv8GEDQLIlZU+N1DgjeQrxCQwGlG2hPaYIHGXjqTLxs4UbYkYQ9kixYpWMy4AHDh0jAvDuTAXHBQrTOOX3/cmfuln3lztH51/VLYBmSj9Pdfd9WIa/8MvXOS+27bx/NnrNyQqWAuqseUuwIgEBIWyR6ESGB1flWyxwqnLwTxOFF6b6f/opq5kaHACo7OnL00yXPOVijkkYg7JWPB6vpgRMuZkvnfYC63TMzW2Ppgp4khgjESEchjK+vG1XDXjQMwVChWfQqY4ybD4GpTA/vA7t3BlrMTpgQwvnB8OJ1Lrn2fZ86trb9riLq+PjuMgODGHGMHYO9e3T1rzExWni7nB8UUGqPbzujJaoDLDoF/9zms8eaKfHT3tnB7I0JFwOTuYJV/2qsdZBE5fzZBOxfi1n72l7n7qfX4AFc/nhfPD05Kh1jKT8Xr+7PVqzaeF0Op5raVkOQkwyl4kHgjmc0oVH89XXh8aD8NqY5y6nOH89dz0xJlAe8Jl7+Y0N2/pqj5u6EwSd4VkzCUZd0jGnCBl1CIUITAjZMzJw988xcBYAU+DcFZvOjmpzPZUHrhrNw/80QkUxUHQMOoV1emZ6T/r1Nh6KQwR+AqvXB6ryj8hnEvRIFXNTEZFgf/y4mV2b+wgky9TbGDtQjoVY7RQ4dSVTNVwiDfhXZ2/nqN/OPDApq7niS7wtR4NBKGLmYxQvuxzdazAULZIqeIzPA4x18GdkqkbqZNpdwpTP7+xfJmLIwVizuxGoFkT62tFNRbRCq8vEhDUzud4fjA3eTI0OK9cznD6aoZCZbps2hHYtbGDmzd3ccuWNHu3dLFzfTupeGhs3Amj06yqN2aEjFk5cmqAU1cyQHARrPgeF67n2N7TNuNFav/ePjqTLoWyX5PdIEU6NV01VsvU2LoAZZ8gHFZzRY4Mgq+T2+uhGixYnatfkPRTGC1U6Eq5jOUnYmfRlr6Co4FXNrU6bC1TL+pzpSd0JTDU0SHGdGLOCCAVc9jTl56zbPrUeZyrmWACeS5xQ7Mm1lejamwmlsrrq3g+hYpPsewFCs1KkBXk1auZqlrtlctjDNXkOKxlfUeiOodz85Y0b9mcpiuVqBqaZMwlGXOCec8lwoyQMSsPPXNy0h24alBj/dJInn27Nsy43Zs3dU27sOVKlVkvbFNj6zHXoeJ5uK5DpeJXE4ZOvSFzpnoNU/DnMEAAFR/aEw7rUzF60ylOFcfwpnhOAng+bO9J4ToTnuBUYzA11c9cy/aKU+5Qy34wZyUQiChCIzZX2fSp8zhRZodo7dLUfUQ0a2J9NajGGqUZXl89mXTJC24CI6XaySsZzg5m6/7+EzGHN/d1hvM4gdHZsi5FKh4LDE48CKnF3CUppjAjZoSMWXn9Wo6YE1x8I4J6N9S9SEUX5JcvjZIpVhCFtoRLOhUjEXPnvLBFsfUjpwZ44I9OII5MrL4WiIuAwN7NXbx8cbR+cZ86RAZsJjZ2xEklYrgCp66M1VUEQTC/lU7FUVVOXx3j7t/+DqcHs8Qdh65UjBM/vk4p3DbhClvWpeZMYTL13bjjVNcliQjx8K50prLpUH8eJ1q4W0s9I9CsifWVqBpbKIvh9Xm+Vo1NFF4bGS9x8spYVSJ96soY48X6a3K297RNqNW2dPGm3k7ak5HBcasprZYbZoSMOXFEcFypSowhSEtTr95PVGAuX/KqHkqu5FHxlU/u39nQha1ak4dAEi3i4GmwOrsc3vK9dHEUgI64w+7eNK9eGas759NoUCEVd7mWLZIvz+63RNmsh7JFMkWPbCmHK4HwYmBKCKTkKT++np91f3WNowTeXTSPtLkrWbds+lwXvfkYgWZMrC931dhishCvL8qzFsmkc6UKZwayNV7OGJdG6q/JSadi3ByKBqK/jelkGFJbuFItYinnt8wIGbOye2MHpweyuE7g3qsGcel4zJmmuIruzq9lKziOEBMHX4Isz5vXpfjGS5er1UFn+2FH+9m8LhX8JxRQX5lqYxyB8ZLPhevjpOIuRa+OBhqqyrXZuDhSIO7ObrIU2NiZIFeqcG28hADFyuKWo3IkEDJ44VzYnr5OssWJLNzzKZu+HIxAo8ZtpUu55zL4tTLpYsUnX6pwaTTPK5cynLwyxqnLY5weqL8mx3WEN/Z2TBiczWlu2thBKuFW53AWUziw1KpGsZpus7Nv3z49fvx4qw+jZRw5NcBnDr9YXWwZVQvtSyfZ2JkkX/YYy5fZ0JHg9OA4SVcoekGRNiGQO3uqbF2Xon8kz64NHZP+k0ZzGLUXny889UO62+KICANjBa5minWPzY0qnc7xE26PB9kTSnMYIteRWQUMSVfo60rRkXB5dSBbVb01MOU0iercFsH8VrS9G2ZtcB2hUPHoSMR45KPvqjvvFIkQShWPTKFCsRKE7z65/43VdUcrhdqLXu1vo1kVX5tF9B31D+fY2t3Gz793F39n93oKZY/h8RI/CvOqRZ7OyCxrciLhwM1burhlSxddbfGqPDrhNlc4EIVxp87n9qVT85Lsi8gJVd03Zz8zQrOz1o0QTP7PNZov055w6Q1ztmUKZfqH88TcoMRCyfOpeFotm+37gSfkhakC9mxKV/c7mAmkyVGeNj+qDgrEHehpTzCSr1CuCQM2SpReBwLj4gh4ns4qEphr3qi3M8Fv3vcOHj16lh9cGA4VbRMhwoXiymQ1XOSRbehIkIi53HfbtmomhdqL9O071/GNH16l4vskXYd17UFWiN7OJIOh5DvuCm/e1LWsPYvFuui1iloBQaHiMV6s8NpgtroA9JXLY5y/Vn9NTlvc5S2b09y8Jc0tW7q4dWsXW7vbq8KBZMyt5ixcKu58+LnqTWCEqjKaL/NXn31fw/tp1AhZOM6Yk9qQSvQDjahdULppXVDS2glLHYgEl/x0Ks5gtsT27snJRofHSwTCMJ3mgZR9GMiWqoXq5kvtNr6vaFQCYpZ9zTbMpnSSzlSs6o1sSie5NFoIS4rfmBHyNDA8ZS9I1RNI2pPVgni//9evTyvRnStV+PapQbb3tE1aFzSYyTOWLwfnL5Avw+tD2RnDKcshDLbSpNxTZdKXRvK8cmmsui7nR1czk4oWRgjhmpxQIn3r1i7evDlNeyJWDam1WqkGS69qNCNkzIt6C0qjC2dXeCEZyhbJl4Oy0lHKnJ72oARDLaXwouvXsTK1oa65JNgOQT63mRaEBvV/FnK2ASJwPVciUygzmi+zpy/NQKbA1nVtDGWLVHT2MN5UT6ceCddB1UNEJmUKj0p076xTJntq+1A2uCGIsn83UsZiOWQ0WM5SblWl5PlB9oGyx2i+zCuXQrXalSDzwMAM4eKe9nhVGn3L1nW8bds61nckFkU40EyWWtXYNCMkIjuAx4HNBEslDqnqV0RkPfAnwC7gHPAPVXVYAt/vK8AHgRzw86r6Qriv+4EvhLv+N6r6WNh+O/CHQBtBme9fVlVdyBhGY0z9gbrhHNHGziQQXLA9X4m7DrduXTepoNvUH3btotOpTBSQE7Z1t82YYiTq483gLiXCORZPPbQBgwaTauYF+xeZlHqoNtHoTRs7yJc9zl/PEaXZjjyZYsVjIFPE08nhwXrkyx4J15km55utRPfU9uiGACbWUk0tY1HLcslosJyk3LUy6Xy5wmsD47x8abQ6j/PaDGty4q6wpy/NLVsDL+ftO9bxhvWReGDxUtwsBUstaGnanJCIbAG2qOoLIpIGTgAfBn4euK6qD4nI54AeVf2siHwQ+OcEBuI9wFdU9T2hQTkO7CO4LpwAbg+NyveAXwaOERihR1T1GRH5t/MZY7bzWK1zQjcShqmdI+pMxhjMFlnXFqfi+VwMJaXbulPEwvLXbXE3uBj7PgnXpT3psqcvzcuXRhkr1Fe0RcRdYef6dsbyJYayZZRAJt3THqfiK8PjJcI8n8DkwFiqRjHk+T4lT4k5wrbu1KzS6aQbGBxfJ8QCqlBRxUFxHIfezgSDmSJFLygXsaE9DiKsa4tPupD6nseVTCDdnslTqz3XmIASyOFjjkM6FeOf3PGGWeeEvLCseanioUiQ/w6qnlCkTpw6x7JYsf/FoPY3tZQqvkgmXQzTJ710cYSTkWLtSlBvqh7butuqwoG3bVvHrVu76EzGm57iZiXR8jkhVb0MXA6fZ0TkJLANuBfYH3Z7DDgCfDZsf1wDq3hMRLpDQ7YfeFZVrwOIyLPA3SJyBOhS1efD9scJjNwz8x0jPNY1Q20YxhX4wflhPvH499nT28nn7rl53vVqILhrCmrgBBfpaL6kVPFRgoWbTihcaMflvbvX893Xr815rF1Jl750itFciUTMoVQJqjuO5sv84k++kSePX6B/pFBVmzlMZCgoVHwcIBl3ERFcCYQP5+dYu1O73ijmBmuUHBHUV8QRKp7PpZFCdRxVGBovk3CFXNGj6Pl0JFx+4c6bePJEPz3tMYays2cQh8C45koergSpgaJ8cW/f3j2tRHfkja3viDOaK1OoeDgiJMKwZKQEdB3oSSXrehbLKQy2FAlAa2XS2WKFVy6P8sOLY1UBwcWR+r+LzmSMvZtD4cC2dbxzxzo2daVakuJmNbIkc0Iisgt4F/BdYFN00VfVyyIS/fK2ARdqNusP22Zr76/TzgLGWFNGKArDeL5yeTQo6OaKcO56btY5gZnmEO67LfjYg5ok4KpPzHWCwljhtq4TxL/FD1LIf/U7r80ZFos5kC35PHDXbj5z+MVqMThfg8Sfjzx3mra4yxvWtzGULU3KPB3hA4WyF3hKApvSCcbyFUqeP+faIaAqffZCebqIoOgkjysKtQUXfp+bN6fJlz0Ov3CRzoTL1dHZjV54aIwXgwW+iZhbnReKwmNfP3jHpO8kSpK6ri3Fxs5A8DGULTCULeE6gquBV+P7gcrws3fvnfadLnUYbDG87/lsG5Wjzpcq/PhajhcvjHDySmBwXr2amfH7DwQCQsxx+Gc/uZv/ad8O2sKMA8tBOLDaaLoREpFO4P8GfkVVx2ZxU+u9oQton/VwGtlGRA4CBwF27tw5xy5XHpEa6fWh8WpJBCWay5k5O3ZkvCqe8vroeHUO4pHnTtObTlYNQEWhMoMSIJozipRgs31hFR8E5dGjZ7keht2mvp8perhOmd29nfzw0mjd/Wg0jsLweJmyr8Qdodygqi2an/JVA9XflOOujWiXPOXMQLaa3WE+6u1ozVOtgGMmlVg9RdloLlDFbUqngqzc4fdz7to4X3jqh+w4OvnivZSx/xstuDfXtrUy6WvjRV68MBrO5QQiguFcfU+0L51k75Y0Z65m8Xyf7vYEruPgSOAVHnl1iAf2v2mRPw2jlqYaIRGJExig/6yqfxY2X41CYGG4bSBs7wd21Gy+HbgUtu+f0n4kbN9ep/9CxpiEqh4CDkEwJ9TwCa8QojBMlJ8MJkpg17voRXeh3zt3nZgERiYmwbqgyNu5PFpfIVRLYHwmLrCNiJt9VU78+Pqs6rKRfIVUptCQlDtKZ9+AYntG5tqmXsr82ag9juhmYCxfpqtt5kJ29UJpUcnzS6N5nDCWV/IVPKVYrvCDC8N84vHjvLmvs+oZLVUdnBsRQdTbdrxY5veOvMYtW7v40ZUML/aPVBeCnhsar/sdpWIOb96c5tatXbxtWzfv3NHN9vVtJGMOP/1/fIf1HalJcznLWSa+mmimOk6ArwEnVfXf1bz1NHA/8FD4+FRN+6dE5AkC0cBoaES+BfyGiPSE/d4PfF5Vr4tIRkTuIAjzfQz4nYWMsdjnvtyJwjBuOHEtgI+ysTM17aJXexeadIV8OTA6bjjH0+hFPF+ennSxkUu1r0zLZl2LE0ydkCvVT+o4E3OJBJaSqUfiK/z4eo7uthhdbYm64bFqldiahaoxx8H3fRwJ5inKNZnHh7LlQCko8PrQ+JJLsWdaC3T66hj3/PujnB0aB+CmDe3T5iVPD2QYLwQebMwR2pMxVKF/OM9P/uaRur8tAXZuaOfmzYFw4B0713Hzli46wjU5UyMyy2l+bK3RTE/oJ4B/ArwkIn8Ttv1LAsPwpIh8AjgPfCR87xsEqrUzBPLpjwOExubLwPfDfg9GIgXgF5mQaD8T/jHfMdYaURjm4W+e4tWBLHEXtqZTFCseV8aCtTAHDh3jgbt2T7oL7etKce5acGdY8X1Emh8fn8tUaLjQc11bfE6l3UpjNF/hn/7ETXUXmB5+4SJtcYexgk+u7JMfK3LHrh6+9+MRhIlif7WqQcepH3ZdigWr9S7yQ9kio4UKI/ks0dz+mcFx/rc//Rt+48Nv4+07u/njY+cZHi9Vy7ZXfKVQqQmteYEB6m6Ls3dLOjA42wMvZ2M62XCKm+UkE19rWNqeOVitEu2I6AJ0+uoYmaLH+o44GzqS1f+EuVKFjoTLUDZYbOrVzHNE61SWg1PRyILQlch7d2+YlrrmwKFjvD6U5dp4CQcJ5tlCBV9fOslgpoinQRYK1wkyMTgCyZiLr4E3cdPGDkbzZb5871uXJG9bvfxw/cP5asG/yGPzwtRNyZhb/b3VQ4CfeOMGPvDWzbxrRw839baTisduKMVNq2TiqxXLHbdIrHYjFDFT/q7z18YpT5TzAW40Sc38ieahVhONLJqN6Ew6vG1bT/WieOfDzwUpjzyt3uUrwXzbnt5OcmU/FJAE67YqvuJKIDdXha3dQVG+vjD/31LlbTtyaoD/+J3XuDCco7czyZnBbFAbZ46FvHE3WPtU8RRflYQrtCdcTnzxA4t6fMbi0vJ1QsbKol7MvuL5kwxQq+zAajNA0LgBAsgWfU5eHuXTh18MkpNmguSkMQccgrQ9lVAef3pwnD29HUi4RmtPXyeZfInLmSICbF6XrHpHtRnLa2mLu5weyHDg0LF5hejqhfV+Ys9GssUKL18c5aVLo3S1xZER4W/6a5SMUz4LkSBX32fv2ctj//0cY/kyHcmJY4yMpLE6MCNkANCZcDkzmA1COBKUYKhdtNlKOzC15MFaZCRfQYBr2RKOTFS3BS+cKwkWpkbZHsqez5fvfWvVcNSGmmprE+04On2u5tp4kUyhwkCm0LCc+sipAX7tqR/iOoEK7fRAhk/+8Qv0ppNcHi1MK2EO4eLi8HuNBCa+BjnXHvoHb2f/3j562hI2V7PKMSO0xjhyaoCHv3lqkhrpg2/bwrUwvANKYZlNrihhWjZaawxbTa2MO/osyn5wAXcdcMWhrytVlT8//M1TcwoO6k3IXx8v09Men1NO7fnKtWyRH5wf5l//11cYHi9R8iZnRB+/NiFx7krFeNfOHt62bR237ezmnTu7efH8CP/2Wz+q/h7fFKrjgKonlk7GqumEbK5m9WFGaA0RFagbzpWraqRXr2b50dXTOGGOtLkKv7WS5XtkS0t0jXcIDZJIteBeOhWErSqez7lrOXZtaA8XJmd54I9OkE7F2NOXnpRuabxYrmbe3tOXZiRXqiakjWiLu5y/luWF88Oc+PEwf/HyVV68OFK3ZMFMJFzh43931yQD8lM3b+Knbt40qd/UxamR91Pr2RmrBzNCa4hHj54lU6iERd6CievoEhJJYKOkncvYFhkhSlBy/badPdNCalfHilVp/Vi+zLXxIIlqrhiE2T59+EUE6GqLs2Vd26Qw16NHz3J1LE/ccRgveeTDXGtlz+cf/N5/n/F4UjGHtoQ7Y3aC4Xx5wYtTW5Hd21gazAitcmoniwczRcqeX61jUlv+QAkyCgftrThSY74oQbXXkVwpWL+lWr2ZAOhKBt/zULYYSLkdKPtKeyIWJOtU6EjGeH1onEI5yFv36T99ka62OGeHZs4U4DpCwhU6kjHaEzE6Ei7DuRLXxkszZp0G8HwaykCw0orcGTeGGaFVzNSwxlC2iK9QDlPtTJ3oN9uz8hjLl0nGXTrjDqPFyZkDxoo+V0bzQXomkWpqJt/3KVd8fKVadloJbj6GxksMhV5ThAC3bu3iH9/xBm5/Qzcf/0/fp6c9MSnrwMbOZLVU+2zZKBrJQGDZC9YWlhJ2FVMb1hARNqVTQSluXdtKs9WCQ5AwtT0RIxOmLYoUZ4kw2/O18RKxsLBdOVz8efJKBi8KvzL95sMR2NSV5KaN7dyyJc3u3g7SqTgH3r2TN2/qYuf6jmmpcvJlD9dxeFNvJzOtFxVoSNX2wF27qwulVYNHU8StXswIrWIuDOdoqyn/3NUWJxZV3GzRMRmLh0+QLPVsTbXPyMP1VauS53xFw7ag/5w3IAp96RSdyTiuMz2p7UxGYndYZbZ3iqgh4sPv3NLQnM7+vX08+KFbgzpS+TJ96dSiZ3Awlg8WjlvF1AtreATpdnb3dnJ2MMv4PBN/GsuPYmX6dzg1JDZ1rVVnwiFf8qndMhKk+MDLl0Zpi7v0poPFrbWhsJlKQAB88emX6UwFkurBbKkqJ0+nYlwZK3Hk1EDDhsiMztrAjNAqpt4akKhk9Fi+bAZolTBb5Yjt3SkEuDhaIOE6QaE7IFf2py28qhWkqAZF4fqH8/S0x/m1n71l0n5nMhKRcSpVfPZuTldLvwcl3huvIWSsHcwIrUJqFXGdiaC0dbTQ7+3bunjqxcs2J7RMiYoGLsbX05l0uTxWxCFIWqoEpcsTrsyYq00gyNMWJhKNucKGjsSMlXbrLYaN+h44dCwo525Sa2MWzAitMuov9AtSuAB85vCLZoCWMWVP2d6don+ksOB9ROXaK56CKmUFUBwJ5onK/oSRm1YllqDSbnvCYXdvZzVTwVQaqXZqUmujEcwIrVCiu9CXLg6TLyuqSmcyRlfSJZmITSrB7TrCQ8+cpKcjOes6DmN5cCMGCIJQWiUsRT61PbRHVeNT737EU2VzZ5AgdCZpdCMLSk1qbTSCGaEVSHQXOpYvkS1OTAiMFSqMFSq0xYpUVHCQavXU04NZurPF6oJUY3VT6+HUMzi1z50p7wtKOhWbVRrdiJezEgrFLUVBP2N2TKK9AonuQmeqJJqvBGGdoudTqPhhUkm4Nl5uqKS2sfKpNTICxKf8T3/D+rZq5moE2hMub1jfHrQ7zpzS6B097XXXCk1T0S1jqXV0Mzc1W/iRUwOtPrQ1hXlCK5DoLtTmdoy5qFYsjcolCMQdh6FsqRqeS7rB/A8EtXpu29kzZ0G7Rr2c5Sy1thx1y4OmeUIi8gciMiAiP6xpWy8iz4rI6fCxJ2wXEXlERM6IyN+KyG0129wf9j8tIvfXtN8uIi+F2zwiYQ6RhYyx0qh3F2oYU4k7geGJhSW1HQkKEVX8YL4o7gZLlgsVn7F8aV6ZCZa7l9MIUxdzgwknWkEzPaE/BP4D8HhN2+eAb6vqQyLyufD1Z4F7gD3h33uArwLvEZH1wJeAfQQ3bSdE5GlVHQ77HASOAd8A7gaeme8YTTv7JhDFr08PZMgUKmu+vo4xM0H6HgdFEVUKFZ+OhEuh7CESZFFHhJj6+KpcGSty286eec2JLGcvpxFMOLE8aJonpKpHgetTmu8FHgufPwZ8uKb9cQ04BnSLyBbgA8Czqno9NDzPAneH73Wp6vOqqgSG7sMLHGNFUBu/3tyVYn1H3HLvGFUiryYi5giOEwhTHNfh3bvW88hH34WngZJSVfF9RQR29LTRl07y9YN3rGijMl8sR93yYKnnhDap6mUAVb0sItEvfhtwoaZff9g2W3t/nfaFjHH5Rk9qKYji1xVPOX09Q2mRFjQaq4O44+D5HnHXoeL5OI6g4b9SWRnJlfj04Rfx/aAguBAUsdva1UbMFfrSqRn3feTUAA89c5LXwyqpuzd28Nm79654gzVT+qGVfl4rjeUiTKh3T68LaF/IGNM7ihwkCPWxc+fOOXa7NFwYzuEKXBjOmyDBqBIPMxrctDEQFgxkCni+MpgpUvJ8UCXmCiP5EqO5Mo4T1PURgYrvc2kkh49QLHvc8++PkilWJkmVj5wa4NOHX2Skphrv6Xskp+4AAA/rSURBVIEsnzn8Ir953ztW/AV7pYcUVwNLLdG+GoXAwsdIC9kP7Kjptx24NEf79jrtCxljGqp6SFX3qeq+3t7eeZ3gYnDk1AAHDh3jzoef48ChYxw5NUA6GTMDZExCgJ62OGP5Mt87d53h8SJj+TKuI9y0sYOd69txHIeNnQnG8hUcR0i4LrHQmnh+oJpb3x7jyliR0wNZXGGSVPnRo2fJFiu4IriOE/4JmUKgIjOMG2WpPaGngfuBh8LHp2raPyUiTxCIBUbDUNq3gN+IFG7A+4HPq+p1EcmIyB3Ad4GPAb+zkDGaeK7z5sipAR7+5ileHcjioDiOw5XRAt97/ZpVOzWqQhQhqG4ac4SRfAVFScWcajqeUiVIPOr5iucrA2NBMcO4G4oRHKFYCXPIAblSUPQOgaFsid29nVWp8oXhHJ6vwfvRcQhUPN9UZMai0DQjJCJfB/YDG0Wkn0Dl9hDwpIh8AjgPfCTs/g3gg8AZIAd8HCA0Nl8Gvh/2e1BVI7HDLxIo8NoIVHHPhO3zGqOV1K7WTidjDGaLZAuViXxfvi0tNQKEoNBcMuZwMUzrozoxL7ixM0l7IsZ4scJApkhPe5xrYYVUL0zTU/KUBDopLp1wnWrlVcLidzAhVd7R085QtoiGIbxgXIg5jqnIjEVBdKZ0ugYA+/bt0+PHjy/6fmsTQLbFXc4MZKn4ih9+HxZ2MyLa4w5v2NBBd3uC/uEcHWFm9FcHsqRiDhs7k3SFKXTODmYpVDxSMZeKF/yeylN+THFXQBVF2NbdxlC2GCQ7lUBVF3lCfekUD9y1e9qckK/Q0x5fFXNCRvMQkROqum+ufstFmLDmmLpa21PFEaoVMA0DgmwGZV8ZzBR45lfumvTegUPHpq1zKVZ8kjXeTcxxAI+KD44j+KqsS8XoTacYzBaJucLGzkTgXSls7kpOkirv39vHb933jknquD29q0MdZywPzAi1iKkJIBOuQ9nzERFETX5tBDhOkIC2VGdSsF7qHNcR1rXHGctXqHjBOiDHcWiPCZvXpehLp6opeaJwcP9wjjf1diAiZIsTHlBkZExBZjQTS2DaIqam3ulNJ6uTxym7NVi1yJTHufDDBG+J2PT/qvVS53xy/xuJuy7pVAwfpeL7+L7S1RabthBz/94+vn7wDr5871uDMh/Fiq2VMZYcu9y1iKl3sa4j9LTHSbpC/2ix1YdnNImYA4gExeUaWHBcqvjEXaG3M1n3/Xpeytu3d/Po0bNUvLFAjBBz2LWhs65xaaQ4nWE0EzNCLWLqau3OZAxVpX+01OpDM5qI4wQejQDr2mKM1KlaGuFKIMX2FQazRY6cGqguIJ2tBs58wmeWSdpoNWaElpipF5D37l7PN14qcvJKptWHZiwBN21o53P33BzkARwrEHNgap1BYcJjiruB+i3mSnVx6GJ6LlaC22g1ZoSWkKmhj3PXsjx/9lqrD8tYQkSk6gU/8EcnEBE6EhMya1Xl5JUMb9mcRmoWiI7lS7xwfjjYBti8LoWI3LDnYpmkjVZjwoQlpDb0kS1WuGxzP2uObDGohrt/bx+37exh5/p2dvd2Vtf55MseHQl3kmglUyhzcaSASCBU8FW5NFIgUwhCeTfiuVgmaaPVmBFaQqIiWgNjBc5ds3DHWqQzOeFxzGQAfuHOmya1XxkNMiRsSqdIuA6CIAKDmeAm5kY8l9VQnM5Y2Vg4bgnZ0dPOuWtZBjLmAa1mZio2GCQqmHhntlICkcKtfziHAtu6U1Vv6dJoHtFgYepieC62DshoJWaEloBIjPDq1TGuj5dtIeoqxQF86hug/7+9Ow+uqjzjOP793ZuFNQRERNkMruOCLGrVYt2oolVpZ7TgWGvtMM7YcaujHR39Q//oWNtOR51alaJtrVap1LEM7iLOaGtRERCVVXAARcGBsEQhCXn6x/ve5CYmMSE3OUnO85m5c892T8775k2ee95zzvuUFmUYXlZKVXXjtOwtBYD85bmREYD6QPTFrj3I9I0HS53raTwIdbL7X13NA69/TM2+Oh+Op5fKAMqIYQNL2La7ulHCwdKsOHJ4GUD9eGzt1fSZslwSOu82c72BXxPqRK+v3MK9C9awt9YDUE9XnBU3TTmCH41vyAgvQhebMuLiccMpzmYZOaQfxx5SxsGDSsMDyP1LOnzB36/buN7Mz4Q60R3PLvfg04OVFmWoMyMj8fBPJtX/068YuprZb66nqjrcyTZzcgXXTzmy0Vhshx4wgMtOGsJb67YVJHW0X7dxvZWncvgWHUnlUHHbc3j19lzHjxhU34WWG/TTOdc2nsohQblvxB6Auq/SogzVtXUt3iQybECJPzPjXBfwIFRg+aMiuO4j99vIZsSI8j4UZTPU7DMumTiCpxdvYtP2rzGgOAPDB/UF8DvPnOsCHoQKxMzYuO1rfv38CnZ8XU11rZ8GdYWsAMWbBDKZ+rHZcoEj/zpNbpDYqup9jQLM9VOOTLQMzqVZ6oKQpKnAfUAWmG1mv9mf/ezcU8OyjZUs3VDJko2VLN1YybYqHwG7EIoyMKhPEbur69ibN7pnRmHEgXOOPpDPd1a36YK/X9B3rntLVRCSlAUeAL4PbALekTTPzD5q7XO1++pY9cUulsag896G7azbWtXs9YSMQgKygaXF9CvJ0rcky7qtu5vNjJkGxVlRlBF7ahquv5T1KWLm5IpGowJ4MjXn0ilVQQg4GVhrZusAJD0FTANaDELrtlZx/J0vNxpQMicjOHzYACaOHsz4UeVMGD2YT7d9xZ3zP6I4q/qUy/vMOHBAMVt3t5w7pjcoysDwgaUM7FvC7jZm6fSg41y6pS0IjQA25s1vAr7T2geqqmvrA9ABA0qYMGowE8eUM35UOeNGljcakBLgqOEDyWTU6Bt+cUbU1Bl7a+rYufebwaw7aGm8s5KsuPasw/26iXOuU6TqOSFJlwLnmdnMOH8FcLKZXddku6uBq+PsUcCqZnY3FPiyLT83UzqgLFs2dDRmRiZbpEym64K/AVidWd0+q63ZU1dV+Xnd3t07C/gT2lwPKeB10cDrokFa62KMmR34bRul7UxoEzAqb34k8FnTjcxsFjCrtR1JerctD2L1dl4PDbwuGnhdNPC6aF3axo57BzhCUoWkEmAGMC/hY3LOudRK1ZmQmdVKuhZ4iXCL9qNm9mHCh+Wcc6mVqiAEYGbPA88XYFetdteliNdDA6+LBl4XDbwuWpGqGxOcc851L2m7JuScc64b8SDUTpKmSlolaa2kW5M+ns4gaZSkhZJWSPpQ0g1x+RBJr0haE98Hx+WSdH+sk/clTczb15Vx+zWSrkyqTB0hKStpiaT5cb5C0qJYpjnxJhcklcb5tXH9oXn7uC0uXyXpvGRK0jGSyiXNlbQyto1TU9wmfhn/Nj6Q9KSkPmltFx1mZv5q44twM8PHwFigBFgGHJP0cXVCOQ8GJsbpgcBq4Bjgt8CtcfmtwD1x+gLgBcIzr6cAi+LyIcC6+D44Tg9Ounz7UR83Af8A5sf5fwIz4vRDwDVx+hfAQ3F6BjAnTh8T20opUBHbUDbpcu1HPfwNmBmnS4DyNLYJwkPv64G+ee3hZ2ltFx19+ZlQ+9QP+2Nm1UBu2J9excw2m9l7cXoXsILwhzeN8I+I+P7DOD0NeMyC/wHlkg4GzgNeMbNtZrYdeAWY2oVF6TBJI4EfALPjvICzgblxk6b1kKufucA5cftpwFNmttfM1gNrCW2px5BUBnwPeATAzKrNrJIUtomoCOgrqQjoB2wmhe2iEDwItU9zw/6MSOhYukTsOpgALAIOMrPNEAIVkBv4raV66Q31dS/wKyA3nPcBQKWZ1cb5/DLVlzeu3xG37w31MBbYCvwldk3OltSfFLYJM/sU+D2wgRB8dgCLSWe76DAPQu3TXKa6Xnt7oaQBwL+AG82staF+WqqXHl1fki4EtpjZ4vzFzWxq37KuR9dDVARMBB40swlAFaH7rSW9ti7ida9phC60Q4D+wPnNbJqGdtFhHoTap03D/vQGkooJAegJM3smLv4idqkQ37fE5S3VS0+vr+8CF0v6hND1ejbhzKg8dsNA4zLVlzeuHwRso+fXA4QybDKzRXF+LiEopa1NAEwB1pvZVjOrAZ4BTiOd7aLDPAi1TyqG/Yn91Y8AK8zsD3mr5gG5u5muBP6dt/yn8Y6oU4AdsWvmJeBcSYPjt8dz47IewcxuM7ORZnYo4Xf9mpldDiwELombNa2HXP1cEre3uHxGvEuqAjgCeLuLilEQZvY5sFHSUXHROYQUKKlqE9EG4BRJ/eLfSq4uUtcuCiLpOyN62otw189qwp0styd9PJ1UxsmEboH3gaXxdQGhH3sBsCa+D4nbi5As8GNgOXBi3r5+Trjguha4KumydaBOzqTh7rixhH8Wa4GngdK4vE+cXxvXj837/O2xflYB5yddnv2sg/HAu7FdPEu4uy2VbQK4C1gJfAD8nXCHWyrbRUdfPmKCc865xHh3nHPOucR4EHLOOZcYD0LOOecS40HIOedcYjwIOeecS4wHIeecc4nxIORcguKDiq9KWippuqQbJfXbz319ImlooY/Ruc6UuvTeznUzE4BiMxsPIZAAjwNfJXlQznUVPxNyrsAk9Zf0nKRlMenZdIVkiCslvRmTvc2XNIwQcMbHM6EbCANiLpS0sJX9Pyjp3ZhU7a4mq2+R9HZ8HR63HyNpQUwut0DSaEmD4plTJm7TT9JGScWSDpP0oqTFkt6QdHQnVZVzHoSc6wRTgc/M7AQzOw54EfgzcBFwOjAcwMy2ADOBN8xsvJndRxjA8iwzO6uV/d9uZicC44AzJI3LW7fTzE4G/kgYbJU4/ZiZjQOeAO43sx2EhGpnxG0uAl6yMCDnLOA6M5sE3Az8qSOV4VxrPAg5V3jLgSmS7pF0OmHI//VmtsbCOFmPd3D/P5b0HrAEOJaQoTPnybz3U+P0qYTMsBDGOZscp+cA0+P0DGBOTN9xGvC0pKXAw4RMu851Cr8m5FyBmdlqSZMIg77eDbxMgfLExNGWbwZOMrPtkv5KGCCz/se3ME0zy+cBd0saAkwCXiPkxqnMXaNyrrP5mZBzBSbpEOArM3uckIHzNKBC0mFxk8ta+fguYGAr68sICeV2SDqIbyZTm573/lac/i/hTAfgcuBNADPbTRjV+T7CCOH7LCQvXC/p0lgWSTqhtfI61xF+JuRc4R0P/E5SHVADXAMMBZ6T9CUhCBzXwmdnAS9I2tzcdSEzWyZpCfAhsA74T5NNSiUtInzBzAW764FHJd1CSNF9Vd72cwhpBs7MW3Y58KCkO4BiQkK/ZW0puHPt5akcnOtiks4EbjazC5M+FueS5t1xzjnnEuNnQs51U7FbrbTJ4ivMbHkSx+NcZ/Ag5JxzLjHeHeeccy4xHoScc84lxoOQc865xHgQcs45lxgPQs455xLzf8WPEgRAo4XkAAAAAElFTkSuQmCC\n"
                    },
                    "metadata": {
                        "needs_background": "light"
                    }
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "\nWe can use the Pandas method <code>corr()</code>  to find the feature other than price that is most correlated with price."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "df.corr()['price'].sort_values()",
            "execution_count": 24,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 24,
                    "data": {
                        "text/plain": "zipcode         -0.053203\nlong             0.021626\ncondition        0.036362\nyr_built         0.054012\nsqft_lot15       0.082447\nsqft_lot         0.089661\nyr_renovated     0.126434\nfloors           0.256794\nwaterfront       0.266369\nlat              0.307003\nbedrooms         0.308797\nsqft_basement    0.323816\nview             0.397293\nbathrooms        0.525738\nsqft_living15    0.585379\nsqft_above       0.605567\ngrade            0.667434\nsqft_living      0.702035\nprice            1.000000\nName: price, dtype: float64"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "# Module 4: Model Development"
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "\nWe can Fit a linear regression model using the  longitude feature <code>'long'</code> and  caculate the R^2."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "X = df[['long']]\nY = df['price']\nlm = LinearRegression()\nlm.fit(X,Y)\nlm.score(X, Y)",
            "execution_count": 25,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 25,
                    "data": {
                        "text/plain": "0.00046769430149007363"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question  6\nFit 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."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "X = df[['sqft_living']]\nY = df['price']\nlm = LinearRegression()\nlm.fit(X,Y)\nlm.score(X,Y)",
            "execution_count": 27,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 27,
                    "data": {
                        "text/plain": "0.49285321790379316"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 7\nFit a linear regression model to predict the <code>'price'</code> using the list of features:"
        },
        {
            "metadata": {},
            "cell_type": "code",
            "source": "features = [\"floors\", \"waterfront\",\"lat\" ,\"bedrooms\" ,\"sqft_basement\" ,\"view\" ,\"bathrooms\",\"sqft_living15\",\"sqft_above\",\"grade\",\"sqft_living\"]\nZ=df[features]\nY=df[\"price\"]\nlm = LinearRegression()\nlm.fit(Z,Y)\nlm.predict(Z)",
            "execution_count": 41,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 41,
                    "data": {
                        "text/plain": "array([283850.64176653, 662015.89176653, 307084.89176653, ...,\n       303822.64176653, 428176.14176653, 303694.64176653])"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "Then calculate the R^2. Take a screenshot of your code."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "lm.score(Z,Y)",
            "execution_count": 42,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 42,
                    "data": {
                        "text/plain": "0.657679183672129"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### This will help with Question 8\n\nCreate 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\nThe 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"
        },
        {
            "metadata": {},
            "cell_type": "code",
            "source": "Input=[('scale',StandardScaler()),('polynomial', PolynomialFeatures(include_bias=False)),('model',LinearRegression())]",
            "execution_count": 46,
            "outputs": []
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 8\nUse 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."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "pipe = Pipeline(Input)\nY=df['price']\nZ=df[features]\npipe.fit(Z, Y)\nypipe = pipe.predict(Z)\nprint(\"The predicted values for the model is:  \", ypipe[0:5])\nprint(\"The Rsqaured value for the model is: \", pipe.score(Z,Y))",
            "execution_count": 55,
            "outputs": [
                {
                    "output_type": "stream",
                    "text": "The predicted values for the model is:   [349649.75 559166.25 449506.75 393246.75 521711.75]\nThe Rsqaured value for the model is:  0.7513408553309376\n",
                    "name": "stdout"
                },
                {
                    "output_type": "stream",
                    "text": "/opt/conda/envs/Python36/lib/python3.6/site-packages/sklearn/preprocessing/data.py:645: DataConversionWarning: Data with input dtype int64, float64 were all converted to float64 by StandardScaler.\n  return self.partial_fit(X, y)\n/opt/conda/envs/Python36/lib/python3.6/site-packages/sklearn/base.py:467: DataConversionWarning: Data with input dtype int64, float64 were all converted to float64 by StandardScaler.\n  return self.fit(X, y, **fit_params).transform(X)\n/opt/conda/envs/Python36/lib/python3.6/site-packages/sklearn/pipeline.py:331: DataConversionWarning: Data with input dtype int64, float64 were all converted to float64 by StandardScaler.\n  Xt = transform.transform(Xt)\n/opt/conda/envs/Python36/lib/python3.6/site-packages/sklearn/pipeline.py:511: DataConversionWarning: Data with input dtype int64, float64 were all converted to float64 by StandardScaler.\n  Xt = transform.transform(Xt)\n",
                    "name": "stderr"
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "# Module 5: Model Evaluation and Refinement"
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "Import the necessary modules:"
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "from sklearn.model_selection import cross_val_score\nfrom sklearn.model_selection import train_test_split\nprint(\"done\")",
            "execution_count": 56,
            "outputs": [
                {
                    "output_type": "stream",
                    "text": "done\n",
                    "name": "stdout"
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "We will split the data into training and testing sets:"
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "features =[\"floors\", \"waterfront\",\"lat\" ,\"bedrooms\" ,\"sqft_basement\" ,\"view\" ,\"bathrooms\",\"sqft_living15\",\"sqft_above\",\"grade\",\"sqft_living\"]    \nX = df[features]\nY = df['price']\n\nx_train, x_test, y_train, y_test = train_test_split(X, Y, test_size=0.15, random_state=1)\n\n\nprint(\"number of test samples:\", x_test.shape[0])\nprint(\"number of training samples:\",x_train.shape[0])",
            "execution_count": 57,
            "outputs": [
                {
                    "output_type": "stream",
                    "text": "number of test samples: 3242\nnumber of training samples: 18371\n",
                    "name": "stdout"
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 9\nCreate 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"
        },
        {
            "metadata": {},
            "cell_type": "code",
            "source": "from sklearn.linear_model import Ridge",
            "execution_count": 59,
            "outputs": []
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "RidgeModel = Ridge(alpha=0.1)\nRidgeModel.fit(x_train,y_train)\nRidgeModel.score(x_train,y_train)",
            "execution_count": 65,
            "outputs": [
                {
                    "output_type": "execute_result",
                    "execution_count": 65,
                    "data": {
                        "text/plain": "0.6594378534950245"
                    },
                    "metadata": {}
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "### Question 10\nPerform 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."
        },
        {
            "metadata": {
                "jupyter": {
                    "outputs_hidden": false
                }
            },
            "cell_type": "code",
            "source": "pr = PolynomialFeatures(degree=2)\nx_train_pr=pr.fit_transform(x_train)\nx_test_pr=pr.fit_transform(x_test)\nRidgeModel2=Ridge(alpha=0.1)\nRidgeModel2.fit(x_train_pr, y_train)\nRsqu_test = RidgeModel2.score(x_test_pr, y_test)\nprint(\"The Rsquared value of the test data is:  \", Rsqu_test)",
            "execution_count": 69,
            "outputs": [
                {
                    "output_type": "stream",
                    "text": "The Rsquared value of the test data is:   0.7002744279699229\n",
                    "name": "stdout"
                }
            ]
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "<p>Once you complete your notebook you will have to share it. Select the icon on the top right a marked in red in the image below, a dialogue box should open, and select the option all&nbsp;content excluding sensitive code cells.</p>\n        <p><img width=\"600\" src=\"https://s3-api.us-geo.objectstorage.softlayer.net/cf-courses-data/CognitiveClass/DA0101EN/coursera/project/save_notebook.png\" alt=\"share notebook\"  style=\"display: block; margin-left: auto; margin-right: auto;\"/></p>\n        <p></p>\n        <p>You can then share the notebook&nbsp; via a&nbsp; URL by scrolling down as shown in the following image:</p>\n        <p style=\"text-align: center;\"><img width=\"600\"  src=\"https://s3-api.us-geo.objectstorage.softlayer.net/cf-courses-data/CognitiveClass/DA0101EN/coursera/project/url_notebook.png\" alt=\"HTML\" style=\"display: block; margin-left: auto; margin-right: auto;\" /></p>\n        <p>&nbsp;</p>"
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "<h2>About the Authors:</h2> \n\n<a href=\"https://www.linkedin.com/in/joseph-s-50398b136/\">Joseph Santarcangelo</a> has a PhD in Electrical Engineering, his research focused on using machine learning, signal processing, and computer vision to determine how videos impact human cognition. Joseph has been working for IBM since he completed his PhD."
        },
        {
            "metadata": {},
            "cell_type": "markdown",
            "source": "Other contributors: <a href=\"https://www.linkedin.com/in/michelleccarey/\">Michelle Carey</a>, <a href=\"www.linkedin.com/in/jiahui-mavis-zhou-a4537814a\">Mavis Zhou</a> "
        },
        {
            "metadata": {},
            "cell_type": "code",
            "source": "",
            "execution_count": null,
            "outputs": []
        }
    ],
    "metadata": {
        "kernelspec": {
            "name": "python3",
            "display_name": "Python 3.6",
            "language": "python"
        },
        "language_info": {
            "name": "python",
            "version": "3.6.9",
            "mimetype": "text/x-python",
            "codemirror_mode": {
                "name": "ipython",
                "version": 3
            },
            "pygments_lexer": "ipython3",
            "nbconvert_exporter": "python",
            "file_extension": ".py"
        },
        "widgets": {
            "state": {},
            "version": "1.1.2"
        }
    },
    "nbformat": 4,
    "nbformat_minor": 4
}