week-6-local-regression-assignment.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Predicting house prices using k-nearest neighbors regression\n",
    "In this notebook, you will implement k-nearest neighbors regression. You will:\n",
    "  * Find the k-nearest neighbors of a given query input\n",
    "  * Predict the output for the query input using the k-nearest neighbors\n",
    "  * Choose the best value of k using a validation set"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Fire up GraphLab Create"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#import graphlab\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import sklearn\n",
    "from sklearn.linear_model import LinearRegression, Lasso, Ridge \n",
    "from sklearn.neighbors import NearestNeighbors\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Load in house sales data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For this notebook, we use a subset of the King County housing dataset created by randomly selecting 40% of the houses in the full dataset."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "dtype_dict = {'bathrooms':float, 'waterfront':int, 'sqft_above':int, 'sqft_living15':float, 'grade':int, \n",
    "              'yr_renovated':int, 'price':float, 'bedrooms':float, 'zipcode':str, 'long':float, 'sqft_lot15':float, \n",
    "              'sqft_living':float, 'floors':float, 'condition':int, 'lat':float, 'date':str, 'sqft_basement':int, \n",
    "              'yr_built':int, 'id':str, 'sqft_lot':int, 'view':int}\n",
    "\n",
    "sales     =  pd.read_csv('kc_house_data_small.csv', dtype= dtype_dict)\n",
    "train      = pd.read_csv('kc_house_data_small_train.csv', dtype= dtype_dict)\n",
    "test       = pd.read_csv('kc_house_data_small_test.csv', dtype= dtype_dict)\n",
    "validation = pd.read_csv('kc_house_data_validation.csv', dtype= dtype_dict)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Import useful functions from previous notebooks"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To efficiently compute pairwise distances among data points, we will convert the SFrame into a 2D Numpy array. First import the numpy library and then copy and paste `get_numpy_data()` from the second notebook of Week 2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def get_numpy_data(df, features, output):\n",
    "    df['constant'] = 1 \n",
    "    features = ['constant'] + features # this is how you combine two lists\n",
    "    feature_matrix = df[features]\n",
    "    output_array = df[output]\n",
    "    return(feature_matrix, output_array)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We will also need the `normalize_features()` function from Week 5 that normalizes all feature columns to unit norm. Paste this function below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def normalize_features(feature_matrix):\n",
    "    norms = np.linalg.norm( feature_matrix, axis= 0 )\n",
    "    normalized =   np.array(feature_matrix/norms)\n",
    "    return ( normalized, norms)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Split data into training, test, and validation sets"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "((5527, 21), (1741, 21), (1435, 21))\n"
     ]
    }
   ],
   "source": [
    "#(train_and_validation, test) = sales.random_split(.8, seed=1) # initial train/test split\n",
    "#(train, validation) = train_and_validation.random_split(.8, seed=1) # split training set into training and validation sets\n",
    "\n",
    "print(train.shape, test.shape , validation.shape)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Extract features and normalize"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Using all of the numerical inputs listed in `feature_list`, transform the training, test, and validation SFrames into Numpy arrays:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "feature_list = ['bedrooms',  \n",
    "                'bathrooms',  \n",
    "                'sqft_living',  \n",
    "                'sqft_lot',  \n",
    "                'floors',\n",
    "                'waterfront',  \n",
    "                'view',  \n",
    "                'condition',  \n",
    "                'grade',  \n",
    "                'sqft_above',  \n",
    "                'sqft_basement',\n",
    "                'yr_built',  \n",
    "                'yr_renovated',  \n",
    "                'lat',  \n",
    "                'long',  \n",
    "                'sqft_living15',  \n",
    "                'sqft_lot15']\n",
    "features_train, output_train = get_numpy_data(train, feature_list, 'price')\n",
    "features_test,   output_test = get_numpy_data(test, feature_list, 'price')\n",
    "features_valid, output_valid = get_numpy_data(validation, feature_list, 'price')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In computing distances, it is crucial to normalize features. Otherwise, for example, the `sqft_living` feature (typically on the order of thousands) would exert a much larger influence on distance than the `bedrooms` feature (typically on the order of ones). We divide each column of the training feature matrix by its 2-norm, so that the transformed column has unit norm.\n",
    "\n",
    "IMPORTANT: Make sure to store the norms of the features in the training set. The features in the test and validation sets must be divided by these same norms, so that the training, test, and validation sets are normalized consistently."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "features_train, norms = normalize_features(features_train) # normalize training set features (columns)\n",
    "features_test  = np.array(features_test / norms) # normalize test set by training set norms\n",
    "features_valid = np.array(features_valid / norms) # normalize validation set by training set norms"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 1. Compute a single distance"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To start, let's just explore computing the \"distance\" between two given houses.  We will take our **query house** to be the first house of the test set and look at the distance between this house and the 10th house of the training set.\n",
    "\n",
    "To see the features associated with the query house, print the first row (index 0) of the **test** feature matrix. You should get an 18-dimensional vector whose components are between 0 and 1."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "feature0 = features_test[0]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now print the 10th row (index 9) of the **training** feature matrix. Again, you get an 18-dimensional vector with components between 0 and 1."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "train_feature9 = features_train[9]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***QUIZ QUESTION ***\n",
    "\n",
    "What is the Euclidean distance between the query house and the 10th house of the training set? \n",
    "\n",
    "Note: Do not use the `np.linalg.norm` function; use `np.sqrt`, `np.sum`, and the power operator (`**`) instead. The latter approach is more easily adapted to computing multiple distances at once."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.059723593713980776"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "distance0_9 = np.sqrt( np.dot( feature0 - train_feature9, feature0- train_feature9 ))\n",
    "distance0_9"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 2. Compute multiple distances"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Of course, to do nearest neighbor regression, we need to compute the distance between our query house and *all* houses in the training set.  \n",
    "\n",
    "To visualize this nearest-neighbor search, let's first compute the distance from our query house (`features_test[0]`) to the first 10 houses of the training set (`features_train[0:10]`) and then search for the nearest neighbor within this small set of houses.  Through restricting ourselves to a small set of houses to begin with, we can visually scan the list of 10 distances to verify that our code for finding the nearest neighbor is working.\n",
    "\n",
    "Write a loop to compute the Euclidean distance from the query house to each of the first 10 houses in the training set."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 0, 10, 20, 30, 40, 50, 60, 70, 80, 90])"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x= np.array(range(10))\n",
    "(x*10).reshape(10,)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0.060274709162955915, 0.085468811476437465, 0.061499464352793153, 0.053402739792943625, 0.05844484060170442, 0.059879215098128345, 0.0546314049677546, 0.055431083236146067, 0.052383627840220298, 0.059723593713980776]\n"
     ]
    }
   ],
   "source": [
    "dist_list =[]\n",
    "for i in range(10):\n",
    "    diff = features_test[0] - features_train[i]\n",
    "    dist_list.append( np.sqrt(np.dot(diff, diff )))\n",
    "\n",
    "print(dist_list)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x10ba9f410>]"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAD8CAYAAAB3u9PLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8XHW9//HXJ/vS7E23tM2kdKML3ZIubC6IFPVSUPaC\niFUEQRD0Il7v9SqKv4sbolYW2WW3gBdl6VUBUShtkzR0L6RLpklLm6aTpUnTbJ/fH3NS0pA2k2Qm\nZybzeT4efXTmzHfOfGag855zvuf7/YqqYowxxsS4XYAxxpjwYIFgjDEGsEAwxhjjsEAwxhgDWCAY\nY4xxWCAYY4wBLBCMMcY4LBCMMcYAFgjGGGMccW4X0BfDhw9Xj8fjdhnGGBNRSkpKDqhqbm/tIioQ\nPB4PxcXFbpdhjDERRUQqAmlnp4yMMcYAFgjGGGMcFgjGGGMACwRjjDEOCwRjjDGABYIxxhiHBYIx\nxhggwEAQkcUisk1EykXkth4eTxSRZ5zHV4uIx9keLyKPisgGEdkiIt/t8pxdzvYyEYmKwQUvrd/L\nvvpmt8swxpge9RoIIhILLAfOBaYBl4nItG7NlgE+VZ0I3AXc6Wy/CEhU1ZnAPOBrnWHh+ISqzlbV\nwgG9iwiwv76Z658s5Z43trtdijHG9CiQI4T5QLmq7lDVFuBpYEm3NkuAR53bK4CzREQABVJFJA5I\nBlqA+qBUHmFKvT4AiisOulyJMcb0LJBAyAN2d7lf6WzrsY2qtgF1QA7+cGgE9gJe4Oeq2vmNqMD/\niUiJiFzT73cQIUq9tQBs3lNPQ3Ory9UYY8xHhbpTeT7QDowBCoBvicgE57HTVXUu/lNR14vImT3t\nQESuEZFiESmurq4OcbmhU1LhIyk+hg6FdU44GGNMOAkkEKqAcV3uj3W29djGOT2UAdQAlwOvqmqr\nqu4H3gIKAVS1yvl7P/AC/vD4CFW9X1ULVbUwN7fXyfrC0pG2djZU1XHBnLHExghrd9lpI2NM+Akk\nENYCk0SkQEQSgEuBF7u1eRG4yrl9IfCaqir+00SfBBCRVGAhsFVEUkUkrcv2TwMbB/pmwtWmPfW0\ntHVw5qThTBudboFgjAlLvQaC0ydwA7AS2AI8q6qbROR2ETnPafYgkCMi5cAtQOelqcuBYSKyCX+w\nPKyq64GRwL9E5F1gDfCSqr4azDcWTkor/B3Kc/OzKPRkUba7lpa2DperMsaYYwW0HoKqvgy83G3b\n97vcbsZ/iWn35x06zvYdwKy+FhupSr0+8jKTGZmexHxPNg+/tYuNe+qYOz7L7dKMMeYoG6kcYqpK\nSYWPefn+L/9CTzYAxXbayBgTZiwQQmxPXTP76o8wd3wmALlpiRQMT2XNTp/LlRljzLEsEEKss/9g\nXn720W2F+VmUVByko0PdKssYYz7CAiHEOscfTB2ddnRbUUE2vqZWtlcfcrEyY4w5lgVCiK3z+pg1\nNpP42A8/6iKnH2HtLjttZIwJHxYIIdTc2s6mPfXMzT/2aiJPTgrDhyXaeARjTFixQAih9ZV1tHUo\n87pdXioiFHmyLBCMMWHFAiGESroMSOuuyJNNpe8we+sOD3ZZxhjTIwuEECr1+igYnkp2asJHHrN+\nBGNMuLFACBFVpbTCd9zRyCePTiM1IdYGqBljwoYFQoh4DzZR09jC3PzMHh+Pi41hbn4Wa3ZaIBhj\nwoMFQoiUHB2Qdvz5igrzs9m2r4G6w7ZgjjHGfRYIIVLq9TEsMY5JI9KO26aoIAvVD0czG2OMmywQ\nQqSkopY54zOJjZHjtpkzLos4WzDHGBMmLBBC4NCRNrZ9UM+cXqa3Tk6IZUZehgWCMSYsWCCEwLu7\na+nQE/cfdCryZPHu7jqaW9sHoTJjjDk+C4QQ6OwTmD2u5yuMuiryZNPS3sGGqrpQl2WMMSdkgRAC\nJV4fk0cOIyM5vte2hUcHqNlpI2OMuywQgqyjQ1nnrQ14eczs1AQmjhhGsY1YNsa4zAIhyHYcOETd\n4dYe5y86niJPFsW7bMEcY4y7LBCCrLSiFiDgIwTwD1Crb27jvf0NoSrLGGN6FVAgiMhiEdkmIuUi\nclsPjyeKyDPO46tFxONsjxeRR0Vkg4hsEZHvBrrPSFVS4SMzJZ4Jw1MDfs78AqcfwaaxMMa4qNdA\nEJFYYDlwLjANuExEpnVrtgzwqepE4C7gTmf7RUCiqs4E5gFfExFPgPuMSKVeH3PGZRJzggFp3Y3N\nSmZkeqLNfGqMcVUgRwjzgXJV3aGqLcDTwJJubZYAjzq3VwBniYgACqSKSByQDLQA9QHuM+LUNbXy\n/v5DAY0/6Mq/YE42a3cdRNX6EYwx7ggkEPKA3V3uVzrbemyjqm1AHZCDPxwagb2AF/i5qh4McJ8R\nZ91uZ0GcPvQfdCryZLO3rpmqWlswxxjjjlB3Ks8H2oExQAHwLRGZ0JcdiMg1IlIsIsXV1dWhqDFo\nSit8xAjMCmBAWndFNh7BGOOyQAKhChjX5f5YZ1uPbZzTQxlADXA58KqqtqrqfuAtoDDAfQKgqver\naqGqFubm5gZQrntKvbVMHZVOamJcn587ZVQaaYlx1o9gjHFNIIGwFpgkIgUikgBcCrzYrc2LwFXO\n7QuB19R/MtwLfBJARFKBhcDWAPcZUdo7lHVeX5/7DzrFxgjzPFl2pZExxjW9BoLTJ3ADsBLYAjyr\nqptE5HYROc9p9iCQIyLlwC1A52Wky4FhIrIJfwg8rKrrj7fPYL6xwfbevgYaW9qPu0JaIIo82by/\n/xC+xpYgVmaMMYEJ6NyGqr4MvNxt2/e73G7Gf4lp9+cd6mn78fYZyY6ukDY+u9/76OxHKKnw8alp\nI4NSlzHGBMpGKgdJqdfH8GEJjMtO7vc+ThmbQUJsjHUsG2NcYYEQJKUVPuaOz8I//KJ/kuJjmTnW\nFswxxrjDAiEIag4dYVdNU58mtDueIk82G6pswRxjzOCzQAiCUq9/Qrv+XmHUVZEni9Z2pWx37YD3\nZYwxfWGBEASlXh9xMcLMvIwB76sw3ya6M8a4wwIhCEoqfEzPyyApPnbA+8pIiWfKyDTWVtgANWPM\n4LJAGKDW9g7WV9Yyrx/zFx1PUUEWpRU+2m3BHGPMILJAGKAte+tpbu0Y0IC07oo82Rw60saWvfVB\n26cxxvTGAmGAjg5IC0KHcieb6M4Y4wYLhAEq9dYyOiOJ0Rn9H5DW3ZjMZPIykym2ie6MMYPIAmGA\nSit8QRl/0F2RJ8sWzDHGDCoLhAH4wFnQpj8L4vSm0JPN/oYjeA82BX3fxhjTEwuEASj1Br//oNOH\n/Qh22sgYMzgsEAagtMJHYlwM00anB33fk0YMIyM53gaoGWMGjQXCAJR4ff4ZSuOC/zHGxAiF+Vms\nrbBAMMYMDguEfmpubWdTVX1I+g86FRVks6O6kQOHjoTsNYwxppMFQj9t2lNHS3tHSK4w6lTk8e/b\nLj81xgwGC4R+Kq3wz0YayiOEGXkZJMbZgjnGmMFhgdBPJRU+xmenkJuWGLLXSIyLZda4TIotEIwx\ng8ACoR9UlRKvj7njgzd/0fHM92SzcU89jUfaQv5axpjoZoHQD5W+w1Q3HAnJ+IPuCj1ZtHfYgjnG\nmNALKBBEZLGIbBORchG5rYfHE0XkGefx1SLicbYvFZGyLn86RGS289gbzj47HxsRzDcWSp0D0uaE\nsP+g07z8LGLEJrozxoRer4EgIrHAcuBcYBpwmYhM69ZsGeBT1YnAXcCdAKr6hKrOVtXZwJXATlUt\n6/K8pZ2Pq+r+ILyfQVFa4SMlIZapo9JC/lppSfFMHZVugWCMCblAjhDmA+WqukNVW4CngSXd2iwB\nHnVurwDOEhHp1uYy57kRr9Rby6yxmcTFDs4ZtyJPFuu8tbS2dwzK6xljolMg32h5wO4u9yudbT22\nUdU2oA7I6dbmEuCpbtsedk4X/VcPARKWmlra2Ly3flD6DzoVFWTT1NLO5j22YI4xJnQG5SeuiCwA\nmlR1Y5fNS1V1JnCG8+fK4zz3GhEpFpHi6urqQaj2xNZX1tHeoUFdIa03tmCOMWYwBBIIVcC4LvfH\nOtt6bCMicUAGUNPl8UvpdnSgqlXO3w3Ak/hPTX2Eqt6vqoWqWpibmxtAuaHVuULanHGDd4QwMj2J\n8dkpFgjGmJAKJBDWApNEpEBEEvB/ub/Yrc2LwFXO7QuB19RZ2UVEYoCL6dJ/ICJxIjLcuR0PfA7Y\nSARY5/UxITeVrNSEQX3dQk8Wxbt8tmCOMSZkeg0Ep0/gBmAlsAV4VlU3icjtInKe0+xBIEdEyoFb\ngK6Xpp4J7FbVHV22JQIrRWQ9UIb/COP3A343IaaqlHprmTcIl5t2N9+TTU1jCzsONA76axtjokNc\nII1U9WXg5W7bvt/ldjNw0XGe+wawsNu2RmBeH2t13a6aJg42toR0QrvjKXT6EYp3HeSk3GGD/vrG\nmKHPRir3QWf/wWBeYdTppNxUslMTbAU1Y0zIWCD0QanXR1pSHBNd+IUu4iyYYx3LxpgQsUDog9IK\nH3PGZxET486QifkF2VTUNLG/vtmV1zfGDG0WCAFqaG5l276GQZnh9HgKj45HsNNGxpjgs0AIUNnu\nWlTd6T/oNH1MOsnxsXbayBgTEhYIASqtqEUEZo9z7wghPjaGOeMzLRCMMSFhgRCgEq+PKSPTSEuK\nd7WOQk82W/bW09Dc6modxpihxwIhAB0dyjqvz5XxB93N92TTof4ZV40xJpgsEAJQXn2IhuY25row\nQrm7OeMziY0RW2fZGBN0FggBKHVxQFp3qYlxTB+TzpqdFgjGmOCyQAhASYWP7NQEPDkpbpcCQGF+\nNmW7a2lpswVzjDHBY4EQgFKvj7njMwmXNXzmF2RxpK2DjXvq3C7FGDOEWCD0oraphe3VjcwJg/6D\nTvPynQFqdtrIGBNEFgi9WOdczRMO/QedctMSmTA81UYsG2OCygKhFyUVPmJjhFPGZrhdyjEKPVkU\nVxyko8MWzDHGBIcFQi9KvT6mjU4nJSGgpSMGTaEnm9qmVrZXH3K7FGPMEGGBcAJt7R2U7a51dUK7\n45nvTHS3xsYjGGOCxALhBLbta6CppT0sRih3l5+TwvBhiRRbP4IxJkgsEE6gc0BaOIxQ7k5EmF+Q\nZQPUjDFBY4FwAqXeWkakJTI2K9ntUnpUmJ9NVe1h9tQedrsUY8wQYIFwAiUVPuaOzwqbAWndzS/o\nXDDHjhKMMQMXUCCIyGIR2SYi5SJyWw+PJ4rIM87jq0XE42xfKiJlXf50iMhs57F5IrLBec6vJcy+\ndasbjuA92BRW4w+6mzoqjdSEWOtHMMYERa+BICKxwHLgXGAacJmITOvWbBngU9WJwF3AnQCq+oSq\nzlbV2cCVwE5VLXOecw/wVWCS82dxEN5P0JR6nf6D/PC7wqhTXGwMc/Oz7AjBGBMUgRwhzAfKVXWH\nqrYATwNLurVZAjzq3F4BnNXDL/7LnOciIqOBdFV9R1UVeAw4v5/vISRKK3wkxMYwfUx4DUjrrsiT\nzbZ9DdQdtgVzjDEDE0gg5AG7u9yvdLb12EZV24A6IKdbm0uAp7q0r+xln64q9fqYnpdOUnys26Wc\nUJEnG9UPr4gyxpj+GpROZRFZADSp6sZ+PPcaESkWkeLq6uoQVPdRLW0dvFtZx7wwvNy0u9njMomL\nERugZowZsEACoQoY1+X+WGdbj21EJA7IAGq6PH4pHx4ddLYf28s+AVDV+1W1UFULc3NzAyh34Dbv\nraelrSMsB6R1l5wQy4y8DFtBzRgzYIEEwlpgkogUiEgC/i/3F7u1eRG4yrl9IfCa0zeAiMQAF+P0\nHwCo6l6gXkQWOn0NXwT+d0DvJIhKwmiFtEDML8jm3d11NLe2u12KMSaC9RoITp/ADcBKYAvwrKpu\nEpHbReQ8p9mDQI6IlAO3AF0vTT0T2K2qO7rt+uvAA0A5sB14ZUDvJIhKvT7yMpMZmZ7kdikBKczP\noqW9gw1VtmCOMab/AprCU1VfBl7utu37XW43Axcd57lvAAt72F4MzOhDrYOmtMJHoTN5XCTorHXN\nzoMURVDdxpjwYiOVu9lTe5i9dc1hOcPp8WSnJjBxxDDrRzDGDIgFQjedA9Iipf+gU5Enm+IKH+22\nYI4xpp8sELopraglKT6Gk0enu11KnxR5smhobuO9fQ1ul2KMiVAWCN2UeH2cMjaT+NjI+mg6+w7s\ntJExpr8i61svxJpb29m8py4s1z/ozdisZEalJ7HGJrozxvSTBUIXG6rqaG3XiOs/AP+COUUF2azd\neRBnCIgxxvSJBUIXH66QFjlXGHVV5Mnig/pmKn22YI4xpu8sELooqfDhyUkhZ1ii26X0S2G+049Q\nYf0Ixpi+s0BwqCql3tqImL/oeKaMSiMtKY41O60fwRjTdxYIjt0HD3Pg0JGI7FDuFBsjzMvPsiuN\njDH9YoHgiNQBad0VebJ5f/8hfI0tbpdijIkwFgiOkgofwxLjmDwyze1SBuToeARbMMcY00cWCI5S\nr4/Z4zKJjem+8mdkOWVsBgmxMbbOsjGmzywQgMYjbWzZWx+xl5t2lRQfyyljMywQjDF9ZoEAvFtZ\nS4cS0VcYdVVUkM3GqjoOt9iCOcaYwFkg8OGAtDnjhkggeLJobVfKdte6XYoxJoJYIACl3lomjRhG\nRkq826UExbzx2YjYRHfGmL6J+kDwD0jzRfT4g+4yUuKZMjKNNRYIxpg+iPpA2HGgkdqm1ogff9Bd\noSeL0gofbe0dbpdijIkQUR8IJZ0T2uVH/hVGXRV5smlsaWfrB7ZgjjEmMFEfCOu8PjKS45kwfJjb\npQRV5wA1u/zUGBOogAJBRBaLyDYRKReR23p4PFFEnnEeXy0ini6PnSIiq0Rkk4hsEJEkZ/sbzj7L\nnD8jgvWm+qKkwsec8ZnERPiAtO7GZCaTl5lsgWCMCVivgSAiscBy4FxgGnCZiEzr1mwZ4FPVicBd\nwJ3Oc+OAx4FrVXU68HGgtcvzlqrqbOfP/oG+mb6qO9zK+/sPMW8IdSh3VeTJYu0uny2YY4wJSCBH\nCPOBclXdoaotwNPAkm5tlgCPOrdXAGeJiACfBtar6rsAqlqjqmEzWqpsdy06hAakdVdUkE11wxG8\nB5vcLsUYEwECCYQ8YHeX+5XOth7bqGobUAfkAJMBFZGVIlIqIrd2e97Dzumi/3ICZFCVVviIEZg1\nbmh1KHfq7EdYs9NOGxljehfqTuU44HRgqfP3BSJylvPYUlWdCZzh/Lmypx2IyDUiUiwixdXV1UEt\nrtTrY8qodIYlxgV1v+FiYu4wMlPiKd5lM58aY3oXSCBUAeO63B/rbOuxjdNvkAHU4D+aeFNVD6hq\nE/AyMBdAVaucvxuAJ/GfmvoIVb1fVQtVtTA3NzfQ99Wr9g6lzFvLvCF2uWlXMTFCYX6WdSwbYwIS\nSCCsBSaJSIGIJACXAi92a/MicJVz+0LgNfX3ZK4EZopIihMUHwM2i0iciAwHEJF44HPAxoG/ncC9\nv7+BhiNtQ2qEck+KPNnsONDIgUNH3C7FGBPmeg0Ep0/gBvxf7luAZ1V1k4jcLiLnOc0eBHJEpBy4\nBbjNea4P+CX+UCkDSlX1JSARWCki653tVcDvg/rOelFa4Z/4baiNUO6usHPBHDtKMMb0IqCT56r6\nMv7TPV23fb/L7WbgouM893H8l5523dYIzOtrscFUUuEjJzWB8dkpbpYRcjPzMkiMi2HtLh+LZ4x2\nuxxjTBiL2pHK67w+5uZn4cLFTYMqIS6G2eMyrR/BGNOrqAyEg40t7DjQOOT7DzoVebLZtKeexiNt\nbpdijAljURkI67z+yzCHev9Bp6KCbNo7lHVeWzDHGHN8URkIJRU+4mKEU8ZmuF3KoJg7PpMYsYnu\njDEnFpWBUOr1MX1MOknxsW6XMijSkuI5eXQ6xRUWCMaY44u6QGhr7+Dd3XVDdv6i4ynyZFNaUUur\nLZhjjDmOqAuErR80cLi1PWo6lDsVebI53NrO5j31bpdijAlTURcInSukRUuHcqcij//9Wj+CMeZ4\nojIQRqUnMSYz2e1SBtWI9CTyc1IsEIwxxxV1gVDq9UXd0UGnwvxsim3BHGPMcURVIOyvb6bSd5g5\n44fuDKcnUuTJosYZlGeMMd1FVSCURtmAtO6KCvwT3a21BXOMMT2IqkAoqfCREBfD9DHRMSCtuwnD\nU8lJTWCtLZhjjOlBVAVCqbeWU/IySIiLqrd9lIhQ6LEFc4wxPYuab8Yjbe1sqIy+AWndFXmy8R5s\nYl99s9ulGGPCTNQEwqY99bS0d0TdgLTuFk7IAeCJdypcrsQYE26iJhBKnQFpc4fwGsqBmJGXwefn\n5LH8je1HZ301xhiIpkDw+hiXncyItCS3S3HdD5ZMZ1R6Ejc/U2ZrJBhjjoqKQFBVSip8UX+6qFN6\nUjy/vHgWFQeb+PFLW9wuxxgTJqIiEPbUNbOv/kjUjj/oyYIJOVxz5gSeWuPlb5v3uV2OMSYMREUg\ndE5oZ0cIx7rl7MmcPDqd255fz4FDR9wuxxjjsoACQUQWi8g2ESkXkdt6eDxRRJ5xHl8tIp4uj50i\nIqtEZJOIbBCRJGf7POd+uYj8WkK42n1phY/k+FimjkoL1UtEpMS4WO6+dDb1zW3c9tx6m+PImCjX\nayCISCywHDgXmAZcJiLTujVbBvhUdSJwF3Cn89w44HHgWlWdDnwcaHWecw/wVWCS82fxQN/M8ZR6\nfcwal0FcbFQcEPXJ5JFpfGfxVP62ZT9Pr93tdjnGGBcF8g05HyhX1R2q2gI8DSzp1mYJ8KhzewVw\nlvOL/9PAelV9F0BVa1S1XURGA+mq+o76f5Y+BpwfhPfTo8/PyeOy+eNDtfuId/WpHk6bmMOP/rKZ\nXTbxnTFRK5BAyAO6/nSsdLb12EZV24A6IAeYDKiIrBSRUhG5tUv7yl72CYCIXCMixSJSXF1dHUC5\nH/Wl0wpYMrvH3RsgJkb4+UWziIsRvvlMGW22zOagaW5t563yA/zy/7bxv2VVbpdjolzcIOz/dKAI\naAL+LiIl+AMjIKp6P3A/QGFhoZ3kDpHRGcncccFMvvHUOpa/vp2bPjXJ7ZKGpObWdtZ5a3lnRw2r\ndtRQ5q2lpUsAl+8/xC1nTyaEXWomwqzaXsOD/9rBLy6aTUZKfEhfK5BAqALGdbk/1tnWU5tKp98g\nA6jB/8v/TVU9ACAiLwNz8fcrjO1ln2aQ/dusMfx9yz5+/dr7fGxKLrPHRfeo7mBoaevg3cpaVm2v\nYdX2Gkq9Po60dRAjMH1MBl86zcOiCTnMHZ/FHS9v5jevlVPb1MoPz5tOTIyFQrTbV9/MN54qJSM5\nntjY0P//EEggrAUmiUgB/i/tS4HLu7V5EbgKWAVcCLymqioiK4FbRSQFaAE+BtylqntFpF5EFgKr\ngS8CvwnKOzID8sMlM1iz8yA3P1PGSzeeTkpCqA8ih5bW9g7WV9bxzo4a3tlRQ/EuH4db2xGBk0el\nc8XCfBZNyKGoIJuM5GN/7d35hVPITEng/jd3UHe4lV9cPIt4uxAiarW2d3DDk6U0tbTz1FcXMiwx\n9P8We30FVW0TkRuAlUAs8JCqbhKR24FiVX0ReBD4g4iUAwfxhwaq6hORX+IPFQVeVtWXnF1/HXgE\nSAZecf4Yl2Ukx/OLi2dz+QPvcMdLW7jjgplulxTW2to72LSnnlU7/EcAxbsO0tjSDsCUkWlcUjSO\nhRNyWDghm8yUhBPuS0T4j8+cTFZKAne+upWG5lZ+t3QeyQmxg/FWTJj56atbWbvLx92XzmbSyMG5\nZF4i6drzwsJCLS4udruMqPCTl7dw/5s7eOhLhXxy6ki3ywkb7R3Klr31/j6A7TWs2XmQBmc+qIkj\nhrFoQg6LTsphQUE2OcMS+/06T6728r0/baAwP4sHrir6yNGEGdpe2bCX654o5apF+fxwyYwB709E\nSlS1sNd2FgimJ0fa2lny27c4cKiFld88Y0BfbpGso0PZtq/B3wewwx8AdYf9Q2kmDE9lgRMACydk\nB33ixL+s38PNz5QxcUQaj315Prlp0fnfINrsqD7Eeb99i4kjhvHs1xYFZUEvCwQzYFs/qOe837zF\nx6bkcv+V86LiyhdV5f39h44eAbyzowZfkz8AxmenHD0CWDghh1EZoZ859x/vVXPtH0oYmZ7IH5Yt\nYFx2Sshf07inqaWNC5a/zf6GZl668QzGZCYHZb+BBoL1GJrjmjoqnVsXT+HHL23h2eLdXFI09Ab3\nqSo7DjQePQJYvaOGA4daAMjLTOask0eyaEIOC0/KIS9I/zj74mOTc3n8Kwu4+uE1XHjv2/xh2QIm\nD9L5ZDO4VJXvvbCR9/Y38NiX5wctDPrCAsGc0JdPK+C1rfv54Z83s3BCDvk5qW6XFDTNre3c/EwZ\nr2z8AIBR6UmcMSn36FFAuPwan5efxbPXLuLKB9dw8X2rePhLRcyxiRqHnCdWe3lhXRU3f2oyZ0zK\ndaUGO2VkerWn9jCLf/Xm0XOaQ2FOqNqmFr7yaDElXh83nTWJJbPz8OSkhPVpMW9NE1c8uJoDh45w\n/5WFnD5puNslmSBZX1nLhfes4tSJOTx0VVHQx6AEesoo8v9lm5Abk5nMj86fQam3lnve2O52OQNW\n6WviC/e8zfrKOn572Vy++anJFAxPDeswABifk8KKaxcxPjuFLz+yllc27HW7JBMEvsYWrnu8lNy0\nRO66eLarAxItEExAlszO47xZY7j77++zvrLW7XL6bfOeej7/u7epbjjCY8vm89lTRrtdUp+MSE/i\nmWsWMSMvneufLOWZtV63SzID0NGh3PxsGdUNR/jd0rlkpZ54rEqoWSCYgP1oyQxy0xL55jNlHHYG\nX0WSt8oPcPF9q4iNEVZcdyoLJ+S4XVK/ZKTE8/hXFnD6pFy+89wG7vtH5B+1Ravfvl7OG9uq+f6/\nTWNWGEwVY4FgApaREs8vLprFjupGfvJyZK3F/L9lVXzp4TXkZSbz/NdPjfgrdVIS4njgi4V87pTR\n/L9XtnLnq1ttgaMI8+Z71dz1t/e4YE4eSxeExxV8FgimT06dOJxlpxfwh3cqeH3bfrfL6ZWqcv+b\n27np6TL869K2AAAPhElEQVTmjvdfrTM6Y/Av5wuFhLgY7r50DpcvGM89b2znP17YSHuHhUIkqKo9\nzE1Pr2PyiDTuuGBG2PRfWSCYPvv3c6YwZWQat65Yz8HGFrfLOa6ODuX2v2zmJy9v5bOnjOaxZfOH\n3BQQsTHCHefP4PpPnMRTa7zc+NQ6WtpsPYtw1tLWwfVPlNLartxzxdywmkDSAsH0WVJ8LHddMpu6\npla++3x4rsXc3NrON55ax8Nv7WLZ6QX85tI5JMYNzUniRIR/P2cq3/vMyby0YS/LHl1LU0ub22WZ\n47jjpc2U7a7lpxeewoTcYW6XcwwLBNMv08ak861PT2blpn38saSy9ycMorqmVr740Bpe2rCX//zs\nyfzX56ZFxdoCXz1zAj+98BTeKj/AFQ+sprYpfI/eotX/llXx6KoKvnJ6AZ+ZGX5XuFkgmH77yhkT\nWFCQzQ9f3MTug01ulwP4B9FddN/brPP6+PVlc/jKGRPcLmlQXVw4jt8tncfGqnouue8d9tU3u12S\ncby/r4HbnttAkSeL75w71e1yemSBYPotNkb4xcWziBHh5mfKXO/Q3PqBf4zB3tpmHr16PufNGuNq\nPW5ZPGMUD19dRKWviQvvfZuKmka3S4p6h460ce3jJaQmxvLby+eG7cJH4VmViRhjs1K4/fzpFFf4\nuNfF6+FXba/hontXoSjPXruIUydG97QOp00czpNfXcih5jYuvHcVW/bWu11S1FJVvvPcenYeaOQ3\nl81lZHroZ8ntLwsEM2Dnz87js6eM5q6/vsfGqrpBf/2/rN/DVQ+tYWR6Es9//TROHp0+6DWEo1nj\nMnn2a4uIFeGS+1ZRUnHQ7ZKi0iNv7+Kl9Xv59jlTWHRSeA+GtEAwAybiv/Rx+DD/KObm1sEbxfzg\nv3byjafWMWtcBiuuXeTKFNXhbNLINFZct4icYYksfWA1b0TA2JGhpKTCxx0vbeFTJ4/k2jNPcruc\nXlkgmKDITEng5xfNonz/If7nla0hf72ODuWOlzbzo79s5pxpo/jDsgW9rlkcrcZmpfDs1xYxYfgw\nvvpYMX9+d4/bJUWFA4eOcP0TpYzJTPb3tUXAlW4WCCZoTp80nKtP8/DI27t4873qkL3OkbZ2bnqm\njN//cydXLcpn+dK5JMUPzTEGwZKblsjTX1vInHFZ3Pj0Op5YXeF2SUNae4dy09Pr8DW1cM8VcyNm\nQKQFggmq7yyeyqQRw/j2H9/FF4JRzPXNrXzpobX8+d093HbuVH5w3nRiI+CXVzhIT4rnsWXz+eSU\nEXzvhY0sf708LAcVAq5fsTZQd/31Pd4qr+FHS2YwfUyG2+UELKAx0yKyGLgbiAUeUNX/6fZ4IvAY\nMA+oAS5R1V0i4gG2ANucpu+o6rXOc94ARgOHncc+rap2gjPCJcXH8qtLZ3P+8rf43p82sPzyuUGb\np2VffTNXPbSG8v2HuOuSWVwwZ2xQ9htNkuJjuffKefz7H9/lZyu3UdvUwn985mRX5tJpPNJGRU0T\nFTWN7KppYteBRnbVNFJR08T+hmbOOnkkt54zhUkRNhHha1v38dvXy7m4cCwXF41zu5w+6TUQRCQW\nWA6cDVQCa0XkRVXd3KXZMsCnqhNF5FLgTuAS57Htqjr7OLtfqqq2BNoQM31MBrecPYU7X93K86VV\nfGHewL+439/XwJceXkttUwsPX13k2hKDQ0F8bAy/vHg2mSkJ/P6fO6k73MpPLpgZkpXwGppbqahp\nOvpF3/mlv6umieqGI8e0HT4sEU9OCqdNHE5aUhzPlVRyzq/e5KJ54/jm2ZMiYlLC3Qeb+ObTZUwb\nnc7tS2a4XU6fBXKEMB8oV9UdACLyNLAE6BoIS4AfOLdXAL+VcJm+z7jimjMn8PrW/fz3i5uYX5A9\noPWJ1+46yLJH1pIYH8szX1vEjLzIOQQPVzExwn//2zQykuO5++/vU3e4lbsvndOvvpi6w61U1DSy\n80DjMV/+FTWNHDh07GnDEWmJeHJS+cSUXPJzUvHkpJKfk0J+TgppSceeZ7/xrEksf72cP6yq4E9l\nVVx9WgHXfewkMlLC83x8c2s71z1RggL3XjEvIvu1el1TWUQuBBar6lec+1cCC1T1hi5tNjptKp37\n24EFwDBgE/AeUA/8p6r+02nzBpADtAPPAT/WHooRkWuAawDGjx8/r6LCOsMixe6DTZx79z+ZNjqd\np65Z2K9z/a9u3MuNT5cxNiuZR6+eHzYL3w8lD7+1kx/+eTOnnpTD/V8sZFjisb8TVZXaptajX/T+\nL37/r/yKmkZ8Ta3HtB+dkUR+TgqenFQ8w1Px5KSQ73zx92dmz90Hm7jrr+/xQlkV6UnxXP+Jk/ji\nIk/YfeF+9/kNPLXGy++/WMjZ00a6Xc4xAl1TOdSB0AAMU9UaEZkH/AmYrqr1IpKnqlUikoY/EB5X\n1cdOVEthYaEWF9sZpkjyXEkl3/rju3xn8VSu+3jfrsN+9O1d/ODPm5gzLpMHrypyfXnBoeyFdZV8\n+4/rmTEmnSsXefB2ntevaWTXgUbqmz+cPVUExmQk4xme4vzK9/9dMDyV8dkpIfui3rynnp+u3Mob\n26oZnZHEzWdP5gtzx4bFRQUrSir59h/f5bqPn8R3FoffPEWBBkIgcV0FdO0ZGets66lNpYjEARlA\njfOL/wiAqpY4QTEZKFbVKmd7g4g8if/U1AkDwUSez8/N4+9b9/HLv27jjEnDAzrdo6r8dOU27nlj\nO2dPG8mvL51DckJ4/Rocai6YM5a0xHiuf7KUb//xXWIE8rKS8eSksmR2Xpdf/CmMzQrdl/6JTBuT\nziNXz+ft7Qe485Wt3LpiPQ/8cwe3njOVs04e4doiM5v31PO9FzawaEIO3zp7sis1BEsgRwhx+E/5\nnIX/i38tcLmqburS5npgpqpe63Qqf15VLxaRXOCgqraLyATgn8BM/KePMlX1gIjEA08Bf1PVe09U\nix0hRCZfYwvn/OpNMpLj+fM3Tj/hl0lLWwe3Pbee59dVsXTBeG5fMiMsfgFGi/31zRw60sbYrBQS\n4sL3qnRV5ZWNH/DzldvYcaCRIk8Wt507lXn52YNaR31zK+f95l80tbTz0o1nkJuWOKivH6hAjxB6\n/S+uqm3ADcBK/JeQPquqm0TkdhE5z2n2IJAjIuXALcBtzvYzgfUiUoa/s/laVT0IJAIrRWQ9UIY/\naH7fp3doIkZWqn8U8/v7D3Hnq8cfxXzoSBvLHl3L8+uq+PanJ/Pj8y0MBtuI9CQm5A4L6zAA/3Qp\nn5k5mpU3n8kdF8xgV00TX7hnFV99rJj39zUMSg2qyreffZfdvsMsXzo3bMOgL3o9QggndoQQ2X7w\n4iYeeXsXf1g2/yOXje6vb+bqR9ay9YMG/ufzM7moMLKu3zbuampp46F/7eTef+ygqaVtUC5Vve8f\n2/l/r2zlPz97ctivuxG0TuVwYoEQ2Zpb2/nsr//JoSNtrPzmmUfnHtpefYirHlrDwcYWfrd0Lh+f\nMsLlSk2kOtjYcvRSVRFCdqnq6h01XP7Aas6ZPjKogy9DJWinjIwJlqT4WO6+dA41h1r43p82oqqU\nVPj4wj1v09zaztPXLLQwMAOSnZrAf31uGn//1sf47MzR3Pfmds782evc/+b2oM3Cu7++mRueWkd+\ndgp3fuGUsA+DvrAjBDPolr9ezs9WbuOKheP5Y3ElozOSeOzLCxifY2MMTHAF+1LVtvYOLn9gNRsq\n6/jT9acxZVRkTKthRwgmbF37sZMozM/i8Xe8TB2dznPXnWphYEKi81LVJ7+6gBFpidy6Yj3n3v0m\nf9u8r18T+/1s5TbW7DzITz4/I2LCoC/sCMG44oO6Zl5YV8VVp+b3a/SqMX010EtVX934Adc+XsIV\nC8fz4/Nnhrja4LJOZWOM6UFrewfPFu/mV397n+qGI5w9rfdZVXceaOS83/yLCbmpPHvtIhLjImug\npJ0yMsaYHsTHxrB0QT7/+PeP8+1PT2bV9hrO+dWbfGfFevbWHf5I+8Mt7Vz3eAmxscLypXMjLgz6\nwgLBGBOVUhLiuOGTk3jz1k/wpVMLeGFdFR//2Rv8zytbqXMm7FNV/vNPG9m2r4FfXTKbsVlDu6/L\nTt4aY6JadmoC3/+3aVx9moe7/voe9725nafWePn6x08iMS6G50orufGsSVFxSbT1IRhjTBddL1UF\nOGPScB65en5ET6MSzNlOjTEmanReqrpqew2vbtzLTZ+aHNFh0BcWCMYY04NFJ+Ww6KQct8sYVNap\nbIwxBrBAMMYY47BAMMYYA1ggGGOMcVggGGOMASwQjDHGOCwQjDHGABYIxhhjHBE1dYWIVAMV/Xz6\ncOBAEMuJdPZ5fMg+i2PZ5/GhofJZ5Ktqbm+NIioQBkJEigOZyyNa2OfxIfssjmWfx4ei7bOwU0bG\nGGMACwRjjDGOaAqE+90uIMzY5/Eh+yyOZZ/Hh6Lqs4iaPgRjjDEnFk1HCMYYY05gyAeCiCwWkW0i\nUi4it7ldj5tEZJyIvC4im0Vkk4jc5HZN4UBEYkVknYj8xe1a3CQimSKyQkS2isgWEVnkdk1uEpGb\nnX8nG0XkKRFJcrumUBvSgSAiscBy4FxgGnCZiExztypXtQHfUtVpwELg+ij/PDrdBGxxu4gwcDfw\nqqpOBWYRxZ+JiOQBNwKFqjoDiAUudbeq0BvSgQDMB8pVdYeqtgBPA0tcrsk1qrpXVUud2w34/8Hn\nuVuVu0RkLPBZ4AG3a3GTiGQAZwIPAqhqi6rWuluV6+KAZBGJA1KAPS7XE3JDPRDygN1d7lcS5V+A\nnUTEA8wBVrtbiet+BdwKdLhdiMsKgGrgYef02QMikup2UW5R1Srg54AX2AvUqer/uVtV6A31QDA9\nEJFhwHPAN1W13u163CIinwP2q2qJ27WEgThgLnCPqs4BGoGo7XMTkSz8ZxMKgDFAqohc4W5VoTfU\nA6EKGNfl/lhnW9QSkXj8YfCEqj7vdj0uOw04T0R24T+d+EkRedzdklxTCVSqaucR4wr8ARGtPgXs\nVNVqVW0FngdOdbmmkBvqgbAWmCQiBSKSgL9T6EWXa3KNiAj+c8RbVPWXbtfjNlX9rqqOVVUP/v83\nXlPVIf8rsCeq+gGwW0SmOJvOAja7WJLbvMBCEUlx/t2cRRR0sse5XUAoqWqbiNwArMR/lcBDqrrJ\n5bLcdBpwJbBBRMqcbf+hqi+7WJMJH98AnnB+PO0Arna5Hteo6moRWQGU4r86bx1RMGrZRiobY4wB\nhv4pI2OMMQGyQDDGGANYIBhjjHFYIBhjjAEsEIwxxjgsEIwxxgAWCMYYYxwWCMYYYwD4/x9p9F95\nqvOwAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10b4a4890>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(range(10), dist_list)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*** QUIZ QUESTION ***\n",
    "\n",
    "Among the first 10 training houses, which house is the closest to the query house?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(8, 0.052383627840220298)\n"
     ]
    }
   ],
   "source": [
    "print(np.argmin(dist_list), min(dist_list) )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It is computationally inefficient to loop over computing distances to all houses in our training dataset. Fortunately, many of the Numpy functions can be **vectorized**, applying the same operation over multiple values or vectors.  We now walk through this process."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Consider the following loop that computes the element-wise difference between the features of the query house (`features_test[0]`) and the first 3 training houses (`features_train[0:3]`):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[  0.00000000e+00  -3.87821276e-03  -1.20498190e-02  -1.05552733e-02\n",
      "   2.08673616e-04  -8.52950206e-03   0.00000000e+00  -5.10236549e-02\n",
      "   0.00000000e+00  -3.47633726e-03  -5.50336860e-03  -2.48168183e-02\n",
      "  -1.63756198e-04   0.00000000e+00  -1.70254220e-05   1.29876855e-05\n",
      "  -5.14364795e-03   6.69281453e-04]\n",
      "[  0.00000000e+00  -3.87821276e-03  -4.51868214e-03  -2.26610387e-03\n",
      "   7.19763456e-04   0.00000000e+00   0.00000000e+00  -5.10236549e-02\n",
      "   0.00000000e+00  -3.47633726e-03   1.30705004e-03  -1.45830788e-02\n",
      "  -1.91048898e-04   6.65082271e-02   4.23090220e-05   6.16364736e-06\n",
      "  -2.89330197e-03   1.47606982e-03]\n",
      "[  0.00000000e+00  -7.75642553e-03  -1.20498190e-02  -1.30002801e-02\n",
      "   1.60518166e-03  -8.52950206e-03   0.00000000e+00  -5.10236549e-02\n",
      "   0.00000000e+00  -5.21450589e-03  -8.32384500e-03  -2.48168183e-02\n",
      "  -3.13866046e-04   0.00000000e+00   4.70885840e-05   1.56292487e-05\n",
      "   3.72914476e-03   1.64764925e-03]\n"
     ]
    }
   ],
   "source": [
    "for i in xrange(3):\n",
    "    print features_train[i]-features_test[0]\n",
    "    # should print 3 vectors of length 18"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The subtraction operator (`-`) in Numpy is vectorized as follows:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[  0.00000000e+00  -3.87821276e-03  -1.20498190e-02  -1.05552733e-02\n",
      "    2.08673616e-04  -8.52950206e-03   0.00000000e+00  -5.10236549e-02\n",
      "    0.00000000e+00  -3.47633726e-03  -5.50336860e-03  -2.48168183e-02\n",
      "   -1.63756198e-04   0.00000000e+00  -1.70254220e-05   1.29876855e-05\n",
      "   -5.14364795e-03   6.69281453e-04]\n",
      " [  0.00000000e+00  -3.87821276e-03  -4.51868214e-03  -2.26610387e-03\n",
      "    7.19763456e-04   0.00000000e+00   0.00000000e+00  -5.10236549e-02\n",
      "    0.00000000e+00  -3.47633726e-03   1.30705004e-03  -1.45830788e-02\n",
      "   -1.91048898e-04   6.65082271e-02   4.23090220e-05   6.16364736e-06\n",
      "   -2.89330197e-03   1.47606982e-03]\n",
      " [  0.00000000e+00  -7.75642553e-03  -1.20498190e-02  -1.30002801e-02\n",
      "    1.60518166e-03  -8.52950206e-03   0.00000000e+00  -5.10236549e-02\n",
      "    0.00000000e+00  -5.21450589e-03  -8.32384500e-03  -2.48168183e-02\n",
      "   -3.13866046e-04   0.00000000e+00   4.70885840e-05   1.56292487e-05\n",
      "    3.72914476e-03   1.64764925e-03]]\n"
     ]
    }
   ],
   "source": [
    "print features_train[0:3] - features_test[0]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that the output of this vectorized operation is identical to that of the loop above, which can be verified below:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.]\n",
      "[ 0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.]\n",
      "[ 0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.]\n"
     ]
    }
   ],
   "source": [
    "# verify that vectorization works\n",
    "results = features_train[0:3] - features_test[0]\n",
    "print results[0] - (features_train[0]-features_test[0])\n",
    "# should print all 0's if results[0] == (features_train[0]-features_test[0])\n",
    "print results[1] - (features_train[1]-features_test[0])\n",
    "# should print all 0's if results[1] == (features_train[1]-features_test[0])\n",
    "print results[2] - (features_train[2]-features_test[0])\n",
    "# should print all 0's if results[2] == (features_train[2]-features_test[0])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Aside: it is a good idea to write tests like this cell whenever you are vectorizing a complicated operation."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Perform 1-nearest neighbor regression\n",
    "\n",
    "Now that we have the element-wise differences, it is not too hard to compute the Euclidean distances between our query house and all of the training houses. \n",
    "First, write a single-line expression to define a variable `diff` such that `diff[i]` gives \n",
    "* **the element-wise difference between the features of the query house and the `i`-th training house.**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "diff = features_train - features_test[0]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To test the code above, run the following cell, which should output a value -0.0934339605842:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "-0.0934339987465\n"
     ]
    }
   ],
   "source": [
    "print diff[-1].sum() # sum of the feature differences between the query and last training house\n",
    "# should print -0.0934339605842"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The next step in computing the Euclidean distances is to take these feature-by-feature differences in `diff`, square each, and take the sum over feature indices.  That is, \n",
    "* **compute the sum of square feature differences for each training house (row in `diff`).**\n",
    "\n",
    "By default, `np.sum` sums up everything in the matrix and returns a single number. To instead sum only over a row or column, we need to specifiy the `axis` parameter described in the `np.sum` [documentation](http://docs.scipy.org/doc/numpy-1.10.1/reference/generated/numpy.sum.html). In particular, `axis=1` computes the sum across each row.\n",
    "\n",
    "Below, we compute this sum of square feature differences for all training houses and verify that the output for the 16th house in the training set is equivalent to having examined only the 16th row of `diff` and computing the sum of squares on that row alone."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.00330705902846\n",
      "0.00330705902846\n"
     ]
    }
   ],
   "source": [
    "print np.sum(diff**2, axis=1)[15] # take sum of squares across each row, and print the 16th sum\n",
    "print np.sum(diff[15]**2) # print the sum of squares for the 16th row -- should be same as above"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "With this result in mind, write a single-line expression to \n",
    "* **compute the Euclidean distances between the query house and all houses in the training set**. Assign the result to a variable `distances`.\n",
    "\n",
    "**Hint**: Do not forget to take the square root of the sum of squares."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "distances = np.sqrt(np.sum( (features_train - features_test[0])**2, axis =1))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To test the code above, run the following cell, which should output a value 0.0237082324496:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.0237082324167\n"
     ]
    }
   ],
   "source": [
    "print distances[100] # Euclidean distance between the query house and the 101th training house\n",
    "# should print 0.0237082324496"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now you are ready to write a function that computes the distances from a query house to all training houses. The function should take two parameters: (i) the matrix of training features and (ii) the single feature vector associated with the query."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [],
   "source": [
    "def cal_distance(features_train, query_house ):\n",
    "    distances = np.sqrt(np.sum( (features_train - query_house)**2, axis =1))\n",
    "    return distances"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.023708232416678195"
      ]
     },
     "execution_count": 40,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cal_distance(features_train, query_house =features_test[0])[100]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*** QUIZ QUESTIONS ***\n",
    "\n",
    "1.  Take the query house to be third house of the test set (`features_test[2]`).  What is the index of the house in the training set that is closest to this query house?\n",
    "2.  What is the predicted value of the query house based on 1-nearest neighbor regression?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "382\n"
     ]
    }
   ],
   "source": [
    "OneNN = np.argmin(cal_distance(features_train, query_house = features_test[2] ) )\n",
    "print(OneNN)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "249000"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "int(output_train[OneNN])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Perform k-nearest neighbor regression"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For k-nearest neighbors, we need to find a *set* of k houses in the training set closest to a given query house. We then make predictions based on these k nearest neighbors."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Fetch k-nearest neighbors\n",
    "\n",
    "Using the functions above, implement a function that takes in\n",
    " * the value of k;\n",
    " * the feature matrix for the training houses; and\n",
    " * the feature vector of the query house\n",
    " \n",
    "and returns the indices of the k closest training houses. For instance, with 2-nearest neighbor, a return value of [5, 10] would indicate that the 6th and 11th training houses are closest to the query house.\n",
    "\n",
    "**Hint**: Look at the [documentation for `np.argsort`](http://docs.scipy.org/doc/numpy/reference/generated/numpy.argsort.html)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def kNN( k, train, query_house):\n",
    "    distances = np.sqrt(np.sum( (train - query_house)**2, axis =1))\n",
    "    indices = np.argsort( distances)[:k]\n",
    "    return indices "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*** QUIZ QUESTION ***\n",
    "\n",
    "Take the query house to be third house of the test set (`features_test[2]`).  What are the indices of the 4 training houses closest to the query house?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([ 382, 1149, 4087, 3142]),\n",
       " array([ 0.0028605 ,  0.00322584,  0.00350216,  0.00359315]))"
      ]
     },
     "execution_count": 48,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "kNN( k =4 , train =  features_train, query_house = features_test[2] )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Make a single prediction by averaging k nearest neighbor outputs"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now that we know how to find the k-nearest neighbors, write a function that predicts the value of a given query house. **For simplicity, take the average of the prices of the k nearest neighbors in the training set**. The function should have the following parameters:\n",
    " * the value of k;\n",
    " * the feature matrix for the training houses;\n",
    " * the output values (prices) of the training houses; and\n",
    " * the feature vector of the query house, whose price we are predicting.\n",
    " \n",
    "The function should return a predicted value of the query house.\n",
    "\n",
    "**Hint**: You can extract multiple items from a Numpy array using a list of indices. For instance, `output_train[[6, 10]]` returns the prices of the 7th and 11th training houses."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def kNN_predict( k, query_house, features_train = features_train,  output_train = output_train):\n",
    "    distances = np.sqrt(np.sum( (features_train - query_house)**2, axis =1))\n",
    "    indices = np.argsort( distances)[:k]\n",
    "    \n",
    "    return output_train[indices].mean()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*** QUIZ QUESTION ***\n",
    "\n",
    "Again taking the query house to be third house of the test set (`features_test[2]`), predict the value of the query house using k-nearest neighbors with `k=4` and the simple averaging method described and implemented above."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "413987"
      ]
     },
     "execution_count": 53,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "int(kNN_predict( k= 4, query_house= features_test[2] ))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "249000"
      ]
     },
     "execution_count": 54,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "int(kNN_predict( k= 1, query_house= features_test[2] ))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Compare this predicted value using 4-nearest neighbors to the predicted value using 1-nearest neighbor computed earlier."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Make multiple predictions"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Write a function to predict the value of *each and every* house in a query set. (The query set can be any subset of the dataset, be it the test set or validation set.) The idea is to have a loop where we take each house in the query set as the query house and make a prediction for that specific house. The new function should take the following parameters:\n",
    " * the value of k;\n",
    " * the feature matrix for the training houses;\n",
    " * the output values (prices) of the training houses; and\n",
    " * the feature matrix for the query set.\n",
    " \n",
    "The function should return a set of predicted values, one for each house in the query set.\n",
    "\n",
    "**Hint**: To get the number of houses in the query set, use the `.shape` field of the query features matrix. See [the documentation](http://docs.scipy.org/doc/numpy-1.10.1/reference/generated/numpy.ndarray.shape.html)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def kNN_predict_set( k, query_set, features_train = features_train,  output_train = output_train):\n",
    "    pred =[]\n",
    "    for query_house in query_set :\n",
    "        distances = np.sqrt(np.sum( (features_train - query_house)**2, axis =1))\n",
    "        indices = np.argsort( distances)[:k]\n",
    "        pred.append(output_train[indices].mean())\n",
    "    return pred"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*** QUIZ QUESTION ***\n",
    "\n",
    "Make predictions for the first 10 houses in the test set using k-nearest neighbors with `k=10`. \n",
    "\n",
    "1. What is the index of the house in this query set that has the lowest predicted value? \n",
    "2. What is the predicted value of this house?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[881300.0, 431860.0, 460595.0, 430200.0, 766750.0, 667420.0, 350032.0, 512800.70000000001, 484000.0, 457235.0]\n",
      "(6, 350032.0)\n"
     ]
    }
   ],
   "source": [
    "pred10 = kNN_predict_set( k=10,  query_set= features_test[:10])\n",
    "print(pred10)\n",
    "print( np.argmin(pred10), min(pred10))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Choosing the best value of k using a validation set"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There remains a question of choosing the value of k to use in making predictions. Here, we use a validation set to choose this value. Write a loop that does the following:\n",
    "\n",
    "* For `k` in [1, 2, ..., 15]:\n",
    "    * Makes predictions for each house in the VALIDATION set using the k-nearest neighbors from the TRAINING set.\n",
    "    * Computes the RSS for these predictions on the VALIDATION set\n",
    "    * Stores the RSS computed above in `rss_all`\n",
    "* Report which `k` produced the lowest RSS on VALIDATION set."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "(Depending on your computing environment, this computation may take 10-15 minutes.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[105453830251561.0, 83445073504025.5, 72692096019202.562, 71946721652091.688, 69846517419718.594, 68899544353180.82, 68341973450051.094, 67361678735491.5, 68372727958976.102, 69335048668556.734, 69523855215598.836, 69049969587246.18, 70011254508263.68, 70908698869034.344, 71106928385945.156]\n"
     ]
    }
   ],
   "source": [
    "kvals = range(1,16)\n",
    "rss_all = []\n",
    "for k in k_list:\n",
    "    pred  = kNN_predict_set( k= k ,  query_set= features_valid , \n",
    "                            features_train = features_train,  output_train = output_train)\n",
    "    RSS = np.dot( pred - output_valid, pred - output_valid) \n",
    "    rss_all.append(RSS)\n",
    "\n",
    "print(rss_all)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To visualize the performance as a function of `k`, plot the RSS on the VALIDATION set for each considered `k` value:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 88,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(67361678735491.5, 8)\n"
     ]
    }
   ],
   "source": [
    "print(min(rss_all), np.argmin(rss_all)+1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x10c507ed0>]"
      ]
     },
     "execution_count": 86,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEDCAYAAADZUdTgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuUXGWd7vHvkzsJt2AaDLmDEWgUuqFFMIB4vBCCgsPo\nLEIYQRKzHMXjcXTNwDBrcJjV6qwZz3g84miWICIBvI0zOJMDskRuA2g6kAQDBkIgIeHWDnINEpL8\nzh/vrulKp7q7ulPdu6r281lrr6rae1fVL1ldT+1697vfVxGBmZkVw6i8CzAzs5Hj0DczKxCHvplZ\ngTj0zcwKxKFvZlYgDn0zswKp29CXdLWk5yT9pop9T5V0v6Qdkj5SYfv+krZI+sbwVGtm1hjqNvSB\na4D5Ve67GbgQuL6P7X8H3Ln3JZmZNba6Df2IuBN4vnydpMMl3SxplaS7JB2Z7ftERKwFdvV+HUnH\nA4cAPx+Jus3M6lndhn4flgGfiYjjgS8A3+xvZ0mjgK9m+5qZFd6YvAuolqR9gXcBP5JUWj1+gKd9\nClgREVvKnmNmVlgNE/qkXyUvRETbIJ5zEnCKpE8B+wLjJL0SEZcMS4VmZnWuYZp3IuIl4HFJHwVQ\ncuwAz1kUETMjYjapiedaB76ZFVndhr6kG4B7gSOy7paLgUXAYklrgHXA2dm+75C0Bfgo8G1J6/Kq\n28ysnslDK5uZFUfdHumbmVnt1d2J3ClTpsTs2bPzLsPMrKGsWrXqdxHRMtB+dRf6s2fPpqurK+8y\nzMwaiqRN1ezn5h0zswJx6JuZFYhD38ysQBz6ZmYF4tA3MyuQpgn95cth9mwYNSrdLl+ed0VmZvWn\n7rpsDsXy5bB0KWzblh5v2pQeAyxalF9dZmb1ZsAj/YGmLcwGPvu6pA2S1ko6rmzbTkmrs+WmWhZe\n7rLLegK/ZNu2tN7MzHpU07xzDf1PW3gGMDdblgL/XLbttYhoy5azhlzlADZvHtx6M7OiGjD0K01b\n2MvZpCGLIyLuAw6UNLVWBVZj5szBrTczK6panMidBjxZ9nhLtg5ggqQuSfdJ+nAN3quizk6YOHH3\ndRMnpvVmZtZjuE/kzoqIrZIOA26T9GBEPNZ7J0lLSU1DzBzC4XnpZO3FF8MLL8CMGfDlL/skrplZ\nb7U40t8KzCh7PD1bR0SUbjcCtwPtlV4gIpZFREdEdLS0DDhIXEWLFsH116f7113nwDczq6QWoX8T\n8LGsF8+JwIsR8bSkyZLGA0iaAswDHqrB+/WpLZs994EHhvNdzMwa14DNO9m0hacBU7IpCS8HxgJE\nxLeAFcACYAOwDfh49tSjSFMX7iJ9uXwlIoY19KdOhUMOgdWrh/NdzMwa14ChHxELB9gewKcrrL8H\nePvQSxua9nYf6ZuZ9aVphmEoaWuDdevg9dfzrsTMrP40Xei3t8OOHSn4zcxsd00Z+uB2fTOzSpou\n9A8/HPbd1+36ZmaVNF3ojxoFxx7r0Dczq6TpQh9SE8+aNbBrV96VmJnVl6YN/Vdegcf2GPDBzKzY\nmjb0wU08Zma9NWXot7bCmDEOfTOz3poy9MePh6OPdrdNM7PemjL0wcMxmJlV0rSh39YGzz4LTz+d\ndyVmZvWjaUPfJ3PNzPbUtKFfGlvf7fpmZj2aNvT33z8NyeAjfTOzHk0b+pCO9h36ZmY9mjr029vT\nVbkvvph3JWZm9WHA0Jd0taTnJP2mj+2S9HVJGyStlXRc2bYLJD2aLRfUsvBqlE7mrl070u9sZlaf\nqjnSvwaY38/2M4C52bIU+GcASQeR5tN9J3ACcLmkyXtT7GC5B4+Z2e4GDP2IuBN4vp9dzgaujeQ+\n4EBJU4HTgVsj4vmI+D1wK/1/edTcm98MBx/s0DczK6lFm/404Mmyx1uydX2tHzFSOtp3t00zs6Qu\nTuRKWiqpS1JXd3d3TV+7vT3Nl7t9e01f1sysIdUi9LcCM8oeT8/W9bV+DxGxLCI6IqKjpaWlBiX1\naGuDN97wROlmZlCb0L8J+FjWi+dE4MWIeBq4BfiApMnZCdwPZOtGlE/mmpn1GDPQDpJuAE4Dpkja\nQuqRMxYgIr4FrAAWABuAbcDHs23PS/o7YGX2UldERH8nhIfFW96SJkp3u76ZWRWhHxELB9gewKf7\n2HY1cPXQSqsNT5RuZtajLk7kDre2tnSk74nSzazoChH6nijdzCwpTOiD2/XNzAoR+kcf7YnSzcyg\nIKE/fjy0tjr0zcwKEfrgidLNzKBgof/ss/DMM3lXYmaWn8KEfmnOXB/tm1mROfTNzAqkMKF/wAFw\n2GHutmlmxVaY0AefzDUzK1Tot7XBhg3w0kt5V2Jmlo9ChX7pytw1a/Ktw8wsL4UMfbfrm1lRFSr0\np071ROlmVmyFCn0ptes79M2sqAoV+uCJ0s2s2KoKfUnzJa2XtEHSJRW2z5L0C0lrJd0uaXrZtp2S\nVmfLTbUsfija29NE6Q89lHclZmYjb8DQlzQauBI4A2gFFkpq7bXbPwLXRsQxwBXAl8u2vRYRbdly\nVo3qHjJfmWtmRVbNkf4JwIaI2BgR24EbgbN77dMK3Jbd/2WF7XVj7lyYNMmhb2bFVE3oTwOeLHu8\nJVtXbg1wTnb/j4D9JL0pezxBUpek+yR9uNIbSFqa7dPV3d09iPIHrzRRurttmlkR1epE7heAd0t6\nAHg3sBXYmW2bFREdwHnA1yQd3vvJEbEsIjoioqOlpaVGJfWtvd0TpZtZMVUT+luBGWWPp2fr/ltE\nPBUR50REO3BZtu6F7HZrdrsRuB1o3/uy905bG7z8MmzcmHclZmYjq5rQXwnMlTRH0jjgXGC3XjiS\npkgqvdalwNXZ+smSxpf2AeYBufebKV2Z63Z9MyuaAUM/InYAFwO3AA8DP4yIdZKukFTqjXMasF7S\nI8AhQGe2/iigS9Ia0gner0RE7qH/trelidLdrm9mRaOIyLuG3XR0dERXV9ewv8+xx8K0abBixbC/\nlZnZsJO0Kjt/2q/CXZFb4uEYzKyIChv67e1pknRPlG5mRVLo0Ae365tZsRQ29I89Nt26icfMiqSw\noX/ggTBnjkPfzIqlsKEPPVfmmpkVReFD/9FH09W5ZmZFUOjQLw2z7InSzawoCh36Ho7BzIqm0KF/\n6KHQ0uJ2fTMrjkKHvpSO9n2kb2ZFUejQh9Su/5vfeKJ0MyuGwoe+J0o3syJx6Hs4BjMrkMKH/lve\nAhMnul3fzIqh8KE/enQah8ehb2ZFUPjQB0+UbmbFUVXoS5ovab2kDZIuqbB9lqRfSFor6XZJ08u2\nXSDp0Wy5oJbF10p7exqK4fHH867EzGx4DRj6kkYDVwJnAK3AQkmtvXb7R+DaiDgGuAL4cvbcg4DL\ngXcCJwCXS5pcu/JrozQcg5t4zKzZVXOkfwKwISI2RsR24Ebg7F77tAK3Zfd/Wbb9dODWiHg+In4P\n3ArM3/uya+ttb0tt+w59M2t21YT+NODJssdbsnXl1gDnZPf/CNhP0puqfC6SlkrqktTV3d1dbe01\nM2ECtLa626aZNb9ancj9AvBuSQ8A7wa2AjurfXJELIuIjojoaGlpqVFJg+PhGMysCKoJ/a3AjLLH\n07N1/y0inoqIcyKiHbgsW/dCNc+tF21t8PTT8OyzeVdiZjZ8qgn9lcBcSXMkjQPOBW4q30HSFEml\n17oUuDq7fwvwAUmTsxO4H8jW1R0Ps2xmRTBg6EfEDuBiUlg/DPwwItZJukLSWdlupwHrJT0CHAJ0\nZs99Hvg70hfHSuCKbF3dKfXgcbu+mTUzRUTeNeymo6Mjurq6cnnvOXPghBPgBz/I5e3NzIZM0qqI\n6BhoP1+RW8Ync82s2Tn0y3iidDNrdg79MqWTuWvX5luHmdlwceiX8XAMZtbsHPplpk2DKVMc+mbW\nvBz6ZUoTpbvbppk1K4d+L+3taaL0N97IuxIzs9pz6PfS1gbbt3uidDNrTg79Xjwcg5k1M4d+L3Pn\nponS3a5vZs3Iod/L6NFwzDE+0jez5uTQr8ATpZtZs3LoV9DeDi+95InSzaz5OPQrKJ3Mdbu+mTUb\nh34FnijdzJqVQ7+CCRPgqKMc+mbWfBz6ffBwDGbWjKoKfUnzJa2XtEHSJRW2z5T0S0kPSForaUG2\nfrak1yStzpZv1fofMFza2+Gpp+C55/KuxMysdgYMfUmjgSuBM4BWYKGk1l67/TVp7tx20sTp3yzb\n9lhEtGXLJ2tU97DzMMtm1oyqOdI/AdgQERsjYjtwI3B2r30C2D+7fwDwVO1KzIdD38yaUTWhPw14\nsuzxlmxduS8C50vaAqwAPlO2bU7W7HOHpFMqvYGkpZK6JHV1d3dXX/0wmjwZZs92u76ZNZdanchd\nCFwTEdOBBcD3JY0CngZmZs0+fw5cL2n/3k+OiGUR0RERHS0tLTUqae+1tflI38yaSzWhvxWYUfZ4\nerau3GLghwARcS8wAZgSEa9HxH9l61cBjwFv3duiR0ppovRXXsm7EjOz2qgm9FcCcyXNkTSOdKL2\npl77bAbeCyDpKFLod0tqyU4EI+kwYC6wsVbFD7f2doiANWvyrsTMrDYGDP2I2AFcDNwCPEzqpbNO\n0hWSzsp2+zzwCUlrgBuACyMigFOBtZJWAz8GPhkRzw/HP2Q4eDgGM2s2StlcPzo6OqKrqyvvMgBY\nvhw+9rE02uasWdDZCYsW5V2VmdmeJK2KiI6B9vMVuX1YvhyWLu0ZXnnTpvR4+fJ86zIz2xsO/T5c\ndhls27b7um3b0nozs0bl0O/D5s2DW29m1ggc+n2YOXNw683MGoFDvw+dnWmC9HL77JPWm5k1Kod+\nHxYtgmXLUq8dKa1buNC9d8yssTn0+7FoETzxBOzcCW9/Ozz4YN4VmZntHYd+FSRYvBhWroS1a/Ou\nxsxs6Bz6VTr/fBg3Dq66Ku9KzMyGzqFfpTe9Cc45B77/ffjDH/KuxsxsaBz6g7BkCfz+9/DTn+Zd\niZnZ0Dj0B+E974E5c+A738m7EjOzoXHoD8KoUXDRRXDbbbCxYQaINjPr4dAfpAsvTOF/9dV5V2Jm\nNngO/UGaPh3OOAO++13YsSPvaszMBsehPwSLF8NTT8HNN+ddiZnZ4Dj0h+CDH4SDD3affTNrPFWF\nvqT5ktZL2iDpkgrbZ0r6paQHJK2VtKBs26XZ89ZLOr2Wxedl7NjUtv+zn8Ezz+RdjZlZ9QYM/Wxi\n8yuBM4BWYKGk1l67/TVp7tx20sTp38ye25o9PhqYD3yzNFF6o7voojQmz/e+l3clZmbVq+ZI/wRg\nQ0RsjIjtwI3A2b32CWD/7P4BwFPZ/bOBGyPi9Yh4HNiQvV7DO+IIOOWU1MRTZ9MMm5n1qZrQnwY8\nWfZ4S7au3BeB8yVtAVYAnxnEc5G0VFKXpK7u7u4qS8/fkiXw6KNw1115V2JmVp1anchdCFwTEdOB\nBcD3JVX92hGxLCI6IqKjpaWlRiUNv498BPbf31fomlnjqCaYtwIzyh5Pz9aVWwz8ECAi7gUmAFOq\nfG7DmjgRzjsPfvQjeOGFvKsxMxtYNaG/EpgraY6kcaQTszf12mcz8F4ASUeRQr872+9cSeMlzQHm\nAr+uVfH1YMmSNOrm9dfnXYmZ2cAGDP2I2AFcDNwCPEzqpbNO0hWSzsp2+zzwCUlrgBuACyNZR/oF\n8BBwM/DpiNg5HP+QvBx3HLS1uc++mTUGRZ11Peno6Iiurq68yxiUK6+Eiy+GVavSl4CZ2UiTtCoi\nOgbaz1fk1sB558H48T7aN7P659CvgcmTU0+e5cvhtdfyrsbMrG8O/RpZsgRefBF+8pO8KzEz65tD\nv0be/W44/HD32Tez+ubQrxEpDbl8xx3pKl0zs3rk0K+hCy6A0aM9q5aZ1S+Hfg0deiiceSZccw28\n8Ube1ZiZ7cmhX2OLF6cx9lesyLsSM7M9OfRrbMECmDrVJ3TNrD459GtszJg0q9aKFbC1aYaWM7Nm\n4dAfBhddBLt2eVYtM6s/Dv1h8Ja3wGmnpWEZdu3Kuxozsx4O/WGyZAls3Ai33553JWZmPRz6w+Sc\nc+DAAz0Im5nVF4f+MNlnHzj//DQWz/PP512NmVni0B9GixfD66+n0TfNzOqBQ38YtbXB8cenPvt1\nNleNmRVUVaEvab6k9ZI2SLqkwvZ/krQ6Wx6R9ELZtp1l23rPrdv0liyBtWvTrFpmZnkbMPQljQau\nBM4AWoGFklrL94mIz0VEW0S0Af8X+Jeyza+VtkXEWRTMwoWpfd9X6JpZPajmSP8EYENEbIyI7cCN\nwNn97L+QNDm6AQccAB/9KFx/Pbz6at7VmFnRVRP604Anyx5vydbtQdIsYA5wW9nqCZK6JN0n6cN9\nPG9ptk9Xd3d3laU3jiVL4OWX4Uc/yrsSMyu6Wp/IPRf4cUTsLFs3K5uh/Tzga5IO7/2kiFgWER0R\n0dHS0lLjkvJ38snw1re6z76Z5a+a0N8KzCh7PD1bV8m59GraiYit2e1G4HagfdBVNrjSrFp33w2/\n/W3e1ZhZkVUT+iuBuZLmSBpHCvY9euFIOhKYDNxbtm6ypPHZ/SnAPOChWhTeaD72sTQCp4/2zSxP\nA4Z+ROwALgZuAR4GfhgR6yRdIam8N865wI0Ru/VIPwrokrQG+CXwlYgoZOi/+c3woQ/BtdfC9u15\nV2NmRaWos6uGOjo6oqurK+8yhsWKFWk6xZ/8JI3NY2ZWK5JWZedP++UrckfQ6afDtGnus29m+XHo\nj6DRo+HjH4ebb4Ynnxx4fzOzWnPoj7CLLkrj8FxzTd6VmFkROfRH2Jw50NoKf/u3MGoUzJ7tUTjN\nbOSMybuAolm+HDZsgJ3Z5WubNsHSpen+okX51WVmxeAj/RF22WV7dtncti2tNzMbbj7SH2GbN1de\nv2kTvO998Pa3wzHHpNvWVpg4cWTrM7Pm5iP9ETZzZuX1kybBSy/Bt7+dTva+4x2w335wxBFplM4r\nroB//Vd47DHYtavyayxfns4R+FyBmfXFR/ojrLMzteFv29azbuLEFPaLFqW2/o0b08QrDz6Ybh94\nIF3QVbqObtKk9Eug/FfB+vXwuc/1vK7PFZhZJb4iNwfLl6c2/M2b05F/Z+fAwfzKK7BuXc8XQelL\nYaBJ12fNgieeqFnpZlanqr0i16HfwCLgqadS+J9xRuV9pL6bg8yseXgYhgKQ0rAO8+enI/pKRo2C\nb30LXn99ZGszs+qN5Pk4h36T6Ozcs6fP+PHpD+jP/gzmzk3nDTzCp1l9Wb48nX/btCn9ei+djxuu\n4HfoN4lFi2DZsnTEL6Xbq66CRx+FW25Jvwg++UmHv9neqNUR+Y4d8OKLsHUr/MVf7N6xA4b32h23\n6RdEBPz853D55fCrX6UTyJddBhdeCOPG5V2dWf0rHZGXB/S4cfCnfwpHH506W5SWl1/e/XHvbX/4\nw8DvN9jzcT6RaxVFpCP/L34xhf+sWSn8L7jA4W/W27ZtcP/96bPyN3+z5xF5bxMmwL779r3st9+e\n6y67DH73uz1fa7A976oNfffTLxgpnfg9/fQU/pdfno5eOjsd/lZsu3alOax/9au0/PrXqWt0aZys\nvkip6/SkSTB27ODfd9KkytfudHYO/rWqUVWbvqT5ktZL2iDpkgrb/0nS6mx5RNILZdsukPRotlxQ\ny+Jt6Erhf999aUavQw5Jf3hHHJEmeXnjjbwrNNs7A7W/P/MM/Nu/pYOd970PJk9OzTQXXQQ33ABv\nehP85V+mfZ55pu8ecjNnwoEHDi3wofL5uGXLhvGiyojodwFGA48BhwHjgDVAaz/7fwa4Ort/ELAx\nu52c3Z/c3/sdf/zxYSNv166IFSsi3vGOCIiYPTviO9+J2L494rrrImbNipDS7XXX5V2tWf+uuy5i\n4sT0t1xaxo+PWLgw4qMfjZg5s2f9mDERxx0X8clPRnz3uxEPPRSxc2d1rzlxYv18HoCuGCDPI5U9\nYOifBNxS9vhS4NJ+9r8HeH92fyHw7bJt3wYW9vd+Dv187doV8R//0RP+U6ZEjBtXv3/oZr3t2hVx\n6KG7/82WL7NmRfzJn0R89asRd98dsW1b9a9dzwdA1YZ+Nc0704Dyyf22ZOv2IGkWMAe4bTDPlbRU\nUpekru7u7ipKsuEiwYIFqU3z3/89DQLnoaAN6ndAv+3bUzPlV78Kf/zHMHVqulK9EimdHP3BD+DP\n/xzmzYN99qn+vRYtSs/ftSvdNuK4VrXup38u8OOIGODUx+4iYllEdERER0tLS41LsqGQ4Mwz+27b\n37QJzjsPvvENWLUq9Tu2+jAc4TzSFxD157/+C372M7jkEjj1VDjgADjpJPjCF2D1anj/++Gggyo/\nt69Rboukmt47W4EZZY+nZ+sqORf4dK/nntbrubdXX57lbebM9AHvbeJEuOOOdMKr9PiEE+Bd70rL\niSemE2E2snr3Jd+0CZYsSRfpnXxy6if+6qu731Za1/v2mWd6Rnkt2bYNPvEJWLkSZszYfZk6FUaP\nrr7mvgYgjIBHHoH//M+eZf36tG3sWDjuuHTF+bx56e9u6tTK/w8wvD1iGsmA/fQljQEeAd5LCvGV\nwHkRsa7XfkcCNwNzsvYlJB0ErAKOy3a7Hzg+IvocG9L99OtLXx+eZcvSkf6TT8I99/Qsq1f3dHE7\n8sj0QTzppHR75JHp6LP0uoMdadT6tn176k++YAH8/veDe+748am/+KRJfd9edVXfz580KX05lBs9\nGg49dM8vg+nTe+4ffHA6aOj99zV+PHz4w2ndPfekI3tIR++lg4p589KcE/01zRTtb6ymF2dJWgB8\njdST5+qI6JR0BenEwU3ZPl8EJkTEJb2eexHwV9nDzoj4bn/v5dCvP4P58Lz6KnR17f5FUBr++cAD\n0xfApEnp53n5IHClL5Jm/lDW0vPPp//b0tHvypX9X+UppV9mvcN80iQYU8Xv/dmzK//imzULHn8c\nXnghHQD0tWzZsuegf+PGpQOEvvrBv/WtPQE/b17qTjzKA8f0yVfkWl0o/Ty/5x649950u25d5X0P\nOigd+R12WPpy8UViSURqnilv4vjtb9O2MWNSE0cpGD/72TSeS297O69Cf7/4qvmijoDu7j2/DP7h\nHyrv7yHBB6/a0B+we89IL+6y2fykvrvTlZZRo1KXuPe8J2Lx4ojOzogbb4z49a8jfve71C2vknru\nUldJpXpfey11Jfz7v48466zUbbb0/zJ5csSZZ0Z86UsRd9yxZ3fD4exLPhz/t7Nm9d2t0gaHKrts\n+kjfRlxfTQXTpsH116fpInsvzz67+777759+EZQvTzwBX//67s0c9dxsVOnoedSodJRbavKYO7fn\nKL7aJo5Gasve218Q1sPNO1a3hvJBf/XV1HZc6Qth48b+J4k59NDUpizV9t8xVJs3w513pl4nr7yy\n5/b994drr03nPw4+eOTrG2mN9CVVzxz6Vtdq+UHftSt1KZw+fc9uhSUtLelI+eST09LePjLnDErn\nNO68E+66K91W+pVTzu3ZNhQOfSucvpqNDjoIPvQhuPtueOyxtG6ffeCd7+z5EjjxxHSRz97auTPN\nWVwe8s89l7Ydcgicckq6oOjUU+Gss9KXXm+ezN6GwkMrW+F0dlZuNvr613t+RTz9dOr9cvfdafny\nl1NQS3DMMT1fAiefnH45lPT1y2T79nRFcink7747zYgEKbxPP70n5OfO3b2J6Utf8gVENvJ8pG9N\nZbDNRq+8ksYZKn0J3Htvz4VGs2al8B83LnUlLT9BPHZsCvHHH4fXXkvrjjyyJ+BPOaW6S/7dnm21\n4uYdsyHYsQPWrOn5Erjrrj17DpWMHQuf+lQK+ZNPLsZJV6tfDn2zGohIQwpU+pj4hKvVk2pD3xc1\nm/VD6ruZxiM2WiNy6JsNoLMznWAt5xOu1qgc+mYDGPE5TM2GkbtsmlVh0SKHvDUHH+mbmRWIQ9/M\nrEAc+mZmBeLQNzMrEIe+mVmB1N0VuZK6gQEGnx1xU4Df5V3EIDRSvY1UKzRWvY1UKzRWvfVY66yI\naBlop7oL/Xokqauay5vrRSPV20i1QmPV20i1QmPV20i19ubmHTOzAnHom5kViEO/OsvyLmCQGqne\nRqoVGqveRqoVGqveRqp1N27TNzMrEB/pm5kViEPfzKxAHPr9kDRD0i8lPSRpnaTP5l3TQCSNlvSA\npH/Pu5aBSDpQ0o8l/VbSw5JOyrumvkj6XPY38BtJN0iakHdN5SRdLek5Sb8pW3eQpFslPZrdTs6z\nxnJ91PsP2d/CWkk/lXRgnjWWVKq1bNvnJYWkKXnUNhQO/f7tAD4fEa3AicCnJbXmXNNAPgs8nHcR\nVfo/wM0RcSRwLHVat6RpwP8EOiLibcBo4Nx8q9rDNcD8XusuAX4REXOBX2SP68U17FnvrcDbIuIY\n4BHg0pEuqg/XsGetSJoBfADYPNIF7Q2Hfj8i4umIuD+7/zIplKblW1XfJE0HzgS+k3ctA5F0AHAq\ncBVARGyPiBfyrapfY4B9JI0BJgJP5VzPbiLiTuD5XqvPBr6X3f8e8OERLaofleqNiJ9HxI7s4X3A\n9BEvrII+/m8B/gn4C6ChesM49KskaTbQDvwq30r69TXSH2EjTNc9B+gGvps1R31H0qS8i6okIrYC\n/0g6onsaeDEifp5vVVU5JCKezu4/AxySZzGDdBHw//Iuoi+Szga2RsSavGsZLId+FSTtC/wE+F8R\n8VLe9VQi6YPAcxGxKu9aqjQGOA7454hoB16lvpof/lvWFn426YvqUGCSpPPzrWpwIvXNbogjUkmX\nkZpWl+ddSyWSJgJ/BfxN3rUMhUN/AJLGkgJ/eUT8S9719GMecJakJ4Abgf8h6bp8S+rXFmBLRJR+\nOf2Y9CVQj94HPB4R3RHxBvAvwLtyrqkaz0qaCpDdPpdzPQOSdCHwQWBR1O9FRIeTDgDWZJ+36cD9\nkt6ca1VVcuj3Q5JIbc4PR8T/zrue/kTEpRExPSJmk04y3hYRdXs0GhHPAE9KOiJb9V7goRxL6s9m\n4ERJE7O/ifdSpyede7kJuCC7fwHwbznWMiBJ80nNk2dFxLa86+lLRDwYEQdHxOzs87YFOC77m657\nDv3+zQNHEEINAAAAh0lEQVT+lHTUvDpbFuRdVBP5DLBc0lqgDfhSzvVUlP0a+TFwP/Ag6XNTV5fh\nS7oBuBc4QtIWSYuBrwDvl/Qo6dfKV/KssVwf9X4D2A+4NfusfSvXIjN91NqwPAyDmVmB+EjfzKxA\nHPpmZgXi0DczKxCHvplZgTj0zcwKxKFvZlYgDn0zswL5/1t6GGhuXsDrAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10c31a090>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "\n",
    "kvals = range(1, 16)\n",
    "plt.plot(kvals, rss_all,'bo-')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***QUIZ QUESTION ***\n",
    "\n",
    "What is the RSS on the TEST data using the value of k found above?  To be clear, sum over all houses in the TEST set."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "133118823551516.81"
      ]
     },
     "execution_count": 90,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pred_test  = kNN_predict_set( k= 8 ,  query_set= features_test , \n",
    "                            features_train = features_train,  output_train = output_train)\n",
    "RSS_test = np.dot( pred_test - output_test, pred_test - output_test) \n",
    "RSS_test"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 92,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "14.124239470850629"
      ]
     },
     "execution_count": 92,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.log(RSS_test)/np.log(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}