Recommendations_with_IBM.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Recommendations with IBM\n",
    "\n",
    "In this notebook, you will be putting your recommendation skills to use on real data from the IBM Watson Studio platform. \n",
    "\n",
    "\n",
    "You may either submit your notebook through the workspace here, or you may work from your local machine and submit through the next page.  Either way assure that your code passes the project [RUBRIC](https://review.udacity.com/#!/rubrics/2322/view).  **Please save regularly.**\n",
    "\n",
    "By following the table of contents, you will build out a number of different methods for making recommendations that can be used for different situations. \n",
    "\n",
    "\n",
    "## Table of Contents\n",
    "\n",
    "I. [Exploratory Data Analysis](#Exploratory-Data-Analysis)<br>\n",
    "II. [Rank Based Recommendations](#Rank)<br>\n",
    "III. [User-User Based Collaborative Filtering](#User-User)<br>\n",
    "IV. [Content Based Recommendations (EXTRA - NOT REQUIRED)](#Content-Recs)<br>\n",
    "V. [Matrix Factorization](#Matrix-Fact)<br>\n",
    "VI. [Extras & Concluding](#conclusions)\n",
    "\n",
    "At the end of the notebook, you will find directions for how to submit your work.  Let's get started by importing the necessary libraries and reading in the data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "[nltk_data] Downloading package wordnet to\n",
      "[nltk_data]     C:\\Users\\ahmed\\AppData\\Roaming\\nltk_data...\n",
      "[nltk_data]   Package wordnet is already up-to-date!\n",
      "[nltk_data] Downloading package stopwords to\n",
      "[nltk_data]     C:\\Users\\ahmed\\AppData\\Roaming\\nltk_data...\n",
      "[nltk_data]   Package stopwords is already up-to-date!\n",
      "[nltk_data] Downloading package abc to\n",
      "[nltk_data]     C:\\Users\\ahmed\\AppData\\Roaming\\nltk_data...\n",
      "[nltk_data]   Package abc is already up-to-date!\n",
      "[nltk_data] Downloading package punkt to\n",
      "[nltk_data]     C:\\Users\\ahmed\\AppData\\Roaming\\nltk_data...\n",
      "[nltk_data]   Unzipping tokenizers\\punkt.zip.\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>article_id</th>\n",
       "      <th>title</th>\n",
       "      <th>email</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1430.0</td>\n",
       "      <td>using pixiedust for fast, flexible, and easier...</td>\n",
       "      <td>ef5f11f77ba020cd36e1105a00ab868bbdbf7fe7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1314.0</td>\n",
       "      <td>healthcare python streaming application demo</td>\n",
       "      <td>083cbdfa93c8444beaa4c5f5e0f5f9198e4f9e0b</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1429.0</td>\n",
       "      <td>use deep learning for image classification</td>\n",
       "      <td>b96a4f2e92d8572034b1e9b28f9ac673765cd074</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1338.0</td>\n",
       "      <td>ml optimization using cognitive assistant</td>\n",
       "      <td>06485706b34a5c9bf2a0ecdac41daf7e7654ceb7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1276.0</td>\n",
       "      <td>deploy your python model as a restful api</td>\n",
       "      <td>f01220c46fc92c6e6b161b1849de11faacd7ccb2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   article_id                                              title  \\\n",
       "0      1430.0  using pixiedust for fast, flexible, and easier...   \n",
       "1      1314.0       healthcare python streaming application demo   \n",
       "2      1429.0         use deep learning for image classification   \n",
       "3      1338.0          ml optimization using cognitive assistant   \n",
       "4      1276.0          deploy your python model as a restful api   \n",
       "\n",
       "                                      email  \n",
       "0  ef5f11f77ba020cd36e1105a00ab868bbdbf7fe7  \n",
       "1  083cbdfa93c8444beaa4c5f5e0f5f9198e4f9e0b  \n",
       "2  b96a4f2e92d8572034b1e9b28f9ac673765cd074  \n",
       "3  06485706b34a5c9bf2a0ecdac41daf7e7654ceb7  \n",
       "4  f01220c46fc92c6e6b161b1849de11faacd7ccb2  "
      ]
     },
     "execution_count": 55,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import project_tests as t\n",
    "import pickle\n",
    "import nltk\n",
    "from nltk.tokenize import word_tokenize\n",
    "from nltk.stem import WordNetLemmatizer\n",
    "from nltk.corpus import stopwords\n",
    "%matplotlib inline\n",
    "nltk.download('wordnet')\n",
    "nltk.download('stopwords')\n",
    "nltk.download('abc')\n",
    "nltk.download('punkt')\n",
    "\n",
    "df = pd.read_csv('data/user-item-interactions.csv')\n",
    "df_content = pd.read_csv('data/articles_community.csv')\n",
    "del df['Unnamed: 0']\n",
    "del df_content['Unnamed: 0']\n",
    "\n",
    "# Show df to get an idea of the data\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>doc_body</th>\n",
       "      <th>doc_description</th>\n",
       "      <th>doc_full_name</th>\n",
       "      <th>doc_status</th>\n",
       "      <th>article_id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Skip navigation Sign in SearchLoading...\\r\\n\\r...</td>\n",
       "      <td>Detect bad readings in real time using Python ...</td>\n",
       "      <td>Detect Malfunctioning IoT Sensors with Streami...</td>\n",
       "      <td>Live</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>No Free Hunch Navigation * kaggle.com\\r\\n\\r\\n ...</td>\n",
       "      <td>See the forest, see the trees. Here lies the c...</td>\n",
       "      <td>Communicating data science: A guide to present...</td>\n",
       "      <td>Live</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>☰ * Login\\r\\n * Sign Up\\r\\n\\r\\n * Learning Pat...</td>\n",
       "      <td>Here’s this week’s news in Data Science and Bi...</td>\n",
       "      <td>This Week in Data Science (April 18, 2017)</td>\n",
       "      <td>Live</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>DATALAYER: HIGH THROUGHPUT, LOW LATENCY AT SCA...</td>\n",
       "      <td>Learn how distributed DBs solve the problem of...</td>\n",
       "      <td>DataLayer Conference: Boost the performance of...</td>\n",
       "      <td>Live</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Skip navigation Sign in SearchLoading...\\r\\n\\r...</td>\n",
       "      <td>This video demonstrates the power of IBM DataS...</td>\n",
       "      <td>Analyze NY Restaurant data using Spark in DSX</td>\n",
       "      <td>Live</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            doc_body  \\\n",
       "0  Skip navigation Sign in SearchLoading...\\r\\n\\r...   \n",
       "1  No Free Hunch Navigation * kaggle.com\\r\\n\\r\\n ...   \n",
       "2  ☰ * Login\\r\\n * Sign Up\\r\\n\\r\\n * Learning Pat...   \n",
       "3  DATALAYER: HIGH THROUGHPUT, LOW LATENCY AT SCA...   \n",
       "4  Skip navigation Sign in SearchLoading...\\r\\n\\r...   \n",
       "\n",
       "                                     doc_description  \\\n",
       "0  Detect bad readings in real time using Python ...   \n",
       "1  See the forest, see the trees. Here lies the c...   \n",
       "2  Here’s this week’s news in Data Science and Bi...   \n",
       "3  Learn how distributed DBs solve the problem of...   \n",
       "4  This video demonstrates the power of IBM DataS...   \n",
       "\n",
       "                                       doc_full_name doc_status  article_id  \n",
       "0  Detect Malfunctioning IoT Sensors with Streami...       Live           0  \n",
       "1  Communicating data science: A guide to present...       Live           1  \n",
       "2         This Week in Data Science (April 18, 2017)       Live           2  \n",
       "3  DataLayer Conference: Boost the performance of...       Live           3  \n",
       "4      Analyze NY Restaurant data using Spark in DSX       Live           4  "
      ]
     },
     "execution_count": 72,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Show df_content to get an idea of the data\n",
    "df_content.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### <a class=\"anchor\" id=\"Exploratory-Data-Analysis\">Part I : Exploratory Data Analysis</a>\n",
    "\n",
    "Use the dictionary and cells below to provide some insight into the descriptive statistics of the data.\n",
    "\n",
    "`1.` What is the distribution of how many articles a user interacts with in the dataset?  Provide a visual and descriptive statistics to assist with giving a look at the number of times each user interacts with an article.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>email</th>\n",
       "      <th>article_counts</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>910</th>\n",
       "      <td>2b6c0f514c2f2b04ad3c4583407dccd0810469ee</td>\n",
       "      <td>364</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2426</th>\n",
       "      <td>77959baaa9895a7e2bdc9297f8b27c1b6f2cb52a</td>\n",
       "      <td>363</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>985</th>\n",
       "      <td>2f5c7feae533ce046f2cb16fb3a29fe00528ed66</td>\n",
       "      <td>170</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3312</th>\n",
       "      <td>a37adec71b667b297ed2440a9ff7dad427c7ac85</td>\n",
       "      <td>169</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2680</th>\n",
       "      <td>8510a5010a5d4c89f5b07baac6de80cd12cfaf93</td>\n",
       "      <td>160</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                         email  article_counts\n",
       "910   2b6c0f514c2f2b04ad3c4583407dccd0810469ee             364\n",
       "2426  77959baaa9895a7e2bdc9297f8b27c1b6f2cb52a             363\n",
       "985   2f5c7feae533ce046f2cb16fb3a29fe00528ed66             170\n",
       "3312  a37adec71b667b297ed2440a9ff7dad427c7ac85             169\n",
       "2680  8510a5010a5d4c89f5b07baac6de80cd12cfaf93             160"
      ]
     },
     "execution_count": 73,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_n = df.groupby(['email']).size().reset_index(name='article_counts')\n",
    "df_n = df_n.sort_values(by=['article_counts'], ascending=False)\n",
    "df_n.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3X28VWWd9/HPV0QpNYE8eiMPgkqZNkl2fJjbpjFtFM0Gu0dLx4K8ncjSO216EJqZQss7qzEbm9QwSS0NHXuQMUtJ0aaZkSdDBMk8oyRHSFAERUcK/M0f13Vocdxnn73g7Icj3/frtV97rd+61tq/vfTsH+vpuhQRmJmZ1WqnZidgZmb9iwuHmZmV4sJhZmaluHCYmVkpLhxmZlaKC4eZmZXiwmFWgaRRkjZIGtBLu2MkdTYqr1Yh6WpJ/9DsPKw5XDisz0kKSQd2i02T9L0G5rBb/uG/o8b2yyW9q2s+Ip6IiN0jYnMdc7xO0hdrbNvQ/dftsz8k6ZfFWEScExFfaEY+1nwuHNZvSdq5yuJTgY3A8ZKGbeM2XjV2lO9pjeHCYQ0naS9Jt0taJ2mtpH+TtFNetq+kH0haI+lxSR8vrDdN0q2SvifpOeBDVT5mEnA1sBg4s9vnL5d0oaTFwAuSvg+MAv41H6V8RtLofOS0c15nqKTvSFop6VlJP+7hu/WYfy/7pOvzJkl6QtLTkv4uLxsPfBZ4f87vwRzfU9K1klZJelLSF7tOreWjhH+XdLmktcA0SQdIukfSM3n7N0oaXMhhpKQf5tyfkfTPkt6U9+Of5s9el9tudbQk6cOSOvJ/z1mS9i0sC0nnSHo077tvSlJedqCk+yStzzndXMv+suZy4bBm+CTQCbQB+5B+FCMXj38FHgSGA8cBF0g6obDuBOBWYDBwY6WNSxoFHJOX3whMrNDsDODdwOCIOAN4AnhPPj31lQrtvwu8FjgE2Bu4vMLn1pJ/b94OvDGv+zlJb4qInwH/H7g553dobns9sAk4EHgrcDzwN4VtHQk8lvO9BBDwJWBf4E3ASGBazn0AcDvwW2B0zn9mRCwDzgH+M3/2YLqRdGze7vuAYXkbM7s1Oxk4HDg0t+vaJ18A7gKGACOAb9S8p6xpXDisGf5A+oHZLyL+EBH/FqnTtMOBtoi4OCJ+HxGPAdcApxfW/c+I+HFEvBwR/93D9icCiyPiYeD7wCGS3tqtzRURsaLKNrbIp7pOBM6JiGdzzvdVaFpL/r25KCL+OyIeJBWgQys1krRPzumCiHghIlaTilnxs1ZGxDciYlPeZkdEzI6IjRGxBvga8Oe57RGkgvLpvL2XImKr6xpVnAnMiIgHImIjMJV0hDK60ObSiFgXEU8Ac4BxOf4HYD9g35KfaU3kwmH1sBkY2C02kPQjAfBVoAO4S9Jjkqbk+H7AvvkU1rp8WuSzpKOSLiuKG82nT7peo3J4IvloJCJWAveRTl0VraB2I4G1EfFsL+1qyb83vytMvwjsXuWzBgKrCp/1LdLRRZfu+2pvSTPzaa3ngO8Be+XFI4HfRsSmErl22Zd0lAFARGwAniEdtXTp6Xt9hnQkNE/SUkn/dxs+3xrMF8ysHp4gne5YVoiNAX4DEBHPk05XfVLSIcAcSfNJP3SPR8TYKtveqjvniNjqh1XS/wbGAlMlfTKH9yAddXyq8MPYvVvoat1ErwCGShocEet6addb/tuqe34rSBf/96ryY999nS/l2Fsi4hlJpwD/XNjeKEk7V9heb11oryQVMiDd0Qa8Hniyl/WIiN8BH87rvR34uaRfRERHb+ta8/iIw+rhZuDvJY2QtJPSba7vIV2bQNLJ+aKogOdIRyibgXnAc/nC9WskDZD0ZkmHl/jsScBs4GDS6ZBxwJtJ1ydOrLLeU8D+lRZExCrgp8CVkoZIGijpHRWa9kX+1fIbna+jdOV0F3CZpNfl/XyApD+vso09gA3AOknDgU93y30VcKnSrcyDJB1d+OwRknbpYbs3AWdJGidpV9L1mLkRsby3LyXpNEkj8uyzpCJVt1ugrW+4cFg9XAz8B/BL0o/BV4AzI2JJXj4W+DnpR+w/gSsj4t78zMR7SD/2jwNPA98G9qzlQyUNIl14/UZE/K7wepx0cbv76aqiL5GK3TpJn6qw/IOkU22/BlYDF3RvsL359+Jf8vszkh7I0xOBXYCHSfv5VtK1o55cBBwGrAd+AvywQu4Hko4YO4H358X3AEuB30l6uvtGI+Ju4B+AH5CKzwHUfl3ncGCupA3ALOD8/N/LWpg8kJOZmZXhIw4zMyul7oUjn+f9laTb8/x1Sg9GLcqvcTkuSVfkh4gWSzqssI1J+eGhRyVVO91gZmZ11oi7qs4n3V3zukLs0xFxa7d2J5LOfY8lPbh0FXCkpKHA54F20oWzhZJm1XBrpJmZ1UFdjzjy3RLvJl0g7M0E4IZI7gcG5wevTgBmR0TXffSzgfF1S9rMzKqq9xHH10kP+OzRLX6JpM8BdwNT8tOmw9n6gaXOHOspvhVJk4HJALvtttvbDjrooL76DmZmO4SFCxc+HRFtvbWrW+GQdDKwOiIWSjqmsGgq6SnSXYDpwIWk2zdVYTNRJb51IGJ63h7t7e2xYMGC7crfzGxHI+m3vbeq76mqo4G/lLSc1OHZsZK+FxGr8umojcB3SH3kQDqSGFlYfwTpidSe4mZm1gR1KxwRMTUiRkTEaNLDQPdExAfydQvyU8OnAF0Phc0CJua7q44C1uenY+8kjakwRNIQUg+gd9YrbzMzq64ZfVXdKKmNdApqEanLZoA7gJNInd+9CJwFEBFrJX0BmJ/bXRwRaxubspmZdXlVPjnuaxxmZuVJWhgR7b2185PjZmZWiguHmZmV4sJhZmaluHCYmVkpLhxmZlaKh46tYPSUn9TUbvml765zJmZmrcdHHGZmVooLh5mZleLCYWZmpbhwmJlZKS4cZmZWiguHmZmV4sJhZmaluHCYmVkpLhxmZlZK3QuHpAGSfiXp9jw/RtJcSY9KulnSLjm+a57vyMtHF7YxNccfkXRCvXM2M7OeNeKI43xgWWH+y8DlETEWeBY4O8fPBp6NiAOBy3M7JB1MGnr2EGA8cKWkAQ3I28zMKqhr4ZA0Ang38O08L+BY4Nbc5HrSuOMAE/I8eflxuf0EYGZEbIyIx0lDyx5Rz7zNzKxn9T7i+DrwGeDlPP96YF1EbMrzncDwPD0cWAGQl6/P7bfEK6yzhaTJkhZIWrBmzZq+/h5mZpbVrXBIOhlYHRELi+EKTaOXZdXW+WMgYnpEtEdEe1tbW+l8zcysNvXsVv1o4C8lnQQMAl5HOgIZLGnnfFQxAliZ23cCI4FOSTsDewJrC/EuxXXMzKzB6nbEERFTI2JERIwmXdy+JyLOBOYAp+Zmk4Db8vSsPE9efk9ERI6fnu+6GgOMBebVK28zM6uuGQM5XQjMlPRF4FfAtTl+LfBdSR2kI43TASJiqaRbgIeBTcC5EbG58WmbmRk0qHBExL3AvXn6MSrcFRURLwGn9bD+JcAl9cvQzMxq5SfHzcysFBcOMzMrxYXDzMxKceEwM7NSXDjMzKwUFw4zMyvFhcPMzEpx4TAzs1JcOMzMrBQXDjMzK8WFw8zMSnHhMDOzUlw4zMysFBcOMzMrxYXDzMxKqeeY44MkzZP0oKSlki7K8eskPS5pUX6Ny3FJukJSh6TFkg4rbGuSpEfza1JPn2lmZvVXz4GcNgLHRsQGSQOBX0r6aV726Yi4tVv7E0nDwo4FjgSuAo6UNBT4PNAOBLBQ0qyIeLaOuZuZWQ/qOeZ4RMSGPDswv6LKKhOAG/J69wODJQ0DTgBmR8TaXCxmA+PrlbeZmVVX12sckgZIWgSsJv34z82LLsmnoy6XtGuODQdWFFbvzLGe4mZm1gR1LRwRsTkixgEjgCMkvRmYChwEHA4MBS7MzVVpE1XiW5E0WdICSQvWrFnTJ/mbmdkrNeSuqohYB9wLjI+IVfl01EbgO8ARuVknMLKw2ghgZZV498+YHhHtEdHe1tZWh29hZmZQ37uq2iQNztOvAd4F/Dpft0CSgFOAJXmVWcDEfHfVUcD6iFgF3AkcL2mIpCHA8TlmZmZNUM+7qoYB10saQCpQt0TE7ZLukdRGOgW1CDgnt78DOAnoAF4EzgKIiLWSvgDMz+0ujoi1dczbzMyqqFvhiIjFwFsrxI/toX0A5/awbAYwo08TNDOzbeInx83MrBQXDjMzK8WFw8zMSnHhMDOzUlw4zMysFBcOMzMrxYXDzMxKceEwM7NSXDjMzKwUFw4zMyvFhcPMzEpx4TAzs1JcOMzMrBQXDjMzK8WFw8zMSnHhMDOzUuo5dOwgSfMkPShpqaSLcnyMpLmSHpV0s6RdcnzXPN+Rl48ubGtqjj8i6YR65WxmZr2r5xHHRuDYiDgUGAeMz2OJfxm4PCLGAs8CZ+f2ZwPPRsSBwOW5HZIOBk4HDgHGA1fm4WjNzKwJ6lY4ItmQZwfmVwDHArfm+PXAKXl6Qp4nLz9OknJ8ZkRsjIjHSWOSH1GvvM3MrLq6XuOQNEDSImA1MBv4L2BdRGzKTTqB4Xl6OLACIC9fD7y+GK+wTvGzJktaIGnBmjVr6vF1zMyMOheOiNgcEeOAEaSjhDdVapbf1cOynuLdP2t6RLRHRHtbW9u2pmxmZr1oyF1VEbEOuBc4Chgsaee8aASwMk93AiMB8vI9gbXFeIV1zMyswXotHJKOlrRbnv6ApK9J2q+G9dokDc7TrwHeBSwD5gCn5maTgNvy9Kw8T15+T0REjp+e77oaA4wF5tX6Bc3MrG/VcsRxFfCipEOBzwC/BW6oYb1hwBxJi4H5wOyIuB24EPhbSR2kaxjX5vbXAq/P8b8FpgBExFLgFuBh4GfAuRGxucbvZ2ZmfWzn3puwKSJC0gTgnyLiWkmTelspIhYDb60Qf4wKd0VFxEvAaT1s6xLgkhpyNTOzOqulcDwvaSrwAeAd+RmKgfVNy8zMWlUtp6reT3qY7+yI+B3pVtiv1jUrMzNrWbUccXwiIi7smomIJyQdUseczMyshdVyxPEXFWIn9nUiZmbWP/R4xCHpo8DHgP3znVFd9gD+o96JmZlZa6p2quom4KfAl8i3xmbPR8TaumZlZmYtq8fCERHrSf1FnZHvpNont99d0u4R8USDcjQzsxbS68VxSecB04CngJdzOIC31C8tMzNrVbXcVXUB8MaIeKbeyZiZWeur5a6qFaRTVmZmZjUdcTwG3CvpJ6QHAQGIiK/VLSszM2tZtRSOJ/Jrl/wyM7MdWK+FIyIuakQiZmbWP9RyV9UcKo+4d2xdMjIzs5ZWy6mqTxWmBwF/BWzqoa2Zmb3K9XpXVUQsLLz+PSL+Fjiyt/UkjZQ0R9IySUslnZ/j0yQ9KWlRfp1UWGeqpA5Jj0g6oRAfn2MdkqZU+jwzM2uMWk5VDS3M7gS8DfhfNWx7E/DJiHhA0h7AQkmz87LLI+Ifu33OwcDpwCHAvsDPJb0hL/4mqbPFTmC+pFkR8XANOZiZWR+r5VTVQtI1DpGKwePA2b2tFBGrgFV5+nlJy0hjefRkAjAzIjYCj+chZLtGCuzIIwciaWZu68JhZtYEtdxVNWZ7P0TSaNIwsnOBo4HzJE0EFpCOSp4lFZX7C6t18sdCs6Jb/BWnyiRNBiYDjBo1antTNjOzHvR6jUPSQEkfl3Rrfp0nqeahYyXtDvwAuCAingOuAg4AxpGOSC7ralph9agS3zoQMT0i2iOiva2trdb0zMyspFpOVV1FGmP8yjz/wRz7m95WzAXmB8CNEfFDgIh4qrD8GuD2PNsJjCysPgJYmad7ipuZWYPVUjgOj4hDC/P3SHqwt5UkCbgWWFbsnkTSsHz9A+C9wJI8PQu4SdLXSBfHxwLzSEccYyWNAZ4kXUD/6xryNjOzOqilcGyWdEBE/BeApP2BzTWsdzTp6OQhSYty7LOk8T3GkU43LQc+AhARSyXdQrrovQk4NyI25888D7gTGADMiIilNX4/MzPrY7UUjk8DcyQ9RvrX/37AWb2tFBG/pPL1iTuqrHMJcEmF+B3V1jMzs8ap5a6quyWNBd5IKgS/zrfMmpnZDqjHwiHpA4Ai4ru5UCzO8Q9LeiEibmpUkmZm1jqq3Y77SeDHFeI352VmZrYDqlY4BkTE892D+VmMmp/jMDOzV5dqhWOgpN26B3O/Ux7QycxsB1WtcFwL3Jq7CwG2dB0yMy8zM7MdUI8XxyPiHyVtAO7L3YYE8AJwaURc1agEzcystVS9HTcirgauzoVDla55mJnZjqWWBwCJiA31TsTMzPqHXnvHNTMzK3LhMDOzUmoZj+O1kv4hd4GOpLGSTq5/amZm1opqOeL4DrAR+NM83wl8sW4ZmZlZS6ulcBwQEV8B/gAQEf9N5V5vzcxsB1BL4fi9pNeQh2uVdADpCMTMzHZAtRSOzwM/A0ZKuhG4G/hMbytJGilpjqRlkpZKOj/Hh0qaLenR/D4kxyXpCkkdkhZLOqywrUm5/aOSJm3TNzUzsz5Ry3gcsyU9ABxFOkV1fkQ8XcO2NwGfjIgHcv9WCyXNBj4E3B0Rl0qaAkwBLgROJA0XOxY4kjSu+ZGShpKKVzvpqGehpFkR8WzJ72pmZn2g2ngch3ULdY0TPkrSqIh4oNqG87jiq/L085KWAcOBCcAxudn1wL2kwjEBuCEiArhf0mBJw3Lb2RGxNuc1GxgPfL/G72hmZn2o2hHHZVWWBXBsrR+SO0d8KzAX2CcXFSJilaS9c7PhwIrCap051lO8+2dMBiYDjBo1qtbUzMyspGqdHL6zLz4g93P1A+CCiHhO6vGGrEoLokp860DEdGA6QHt7+yuWm5lZ36jlAcBzJQ0uzA+R9LFaNi5pIKlo3BgRP8zhp/IpKPL76hzvBEYWVh8BrKwSNzOzJqjlrqoPR8S6rpl8UfrDva2kdGhxLbAsIr5WWDQL6LozahJwWyE+Md9ddRSwPp/SuhM4PhesIcDxOWZmZk1QS++4O0lSvmiNpAHUNgLg0cAHgYckLcqxzwKXArdIOht4AjgtL7sDOAnoAF4EzgKIiLWSvgDMz+0u7rpQbmZmjVdL4biT9EN/Nenawjmk5zqqiohf0vMT5sdVaB/AuT1sawYwo4ZczcyszmopHBcCHwE+SioEdwHfrmdSZmbWump5APBl0sN4Hi7WzMyqPgB4S0S8T9JDVL799S11zczMzFpStSOO8/O7x94wM7Mterwdt+vpbuBjEfHb4guo6TkOMzN79anlOY6/qBA7sa8TMTOz/qHaNY6Pko4sDpC0uLBoD+Df652YmZm1pmrXOG4Cfgp8idT1eZfn/QCemdmOq1onh+slPQ/8Sb6uYWZmVv0aR36G40FJ7qfczMyA2p4cHwYslTQPeKErGBF/WbeszMysZdVSOC6qexZmZtZv1NLlyH3FeUlHA38N3Fd5DTMzezWr5YgDSeNIxeJ9wOOkwZnMzGwHVO05jjcApwNnAM8ANwPqqyFlzcysf6p2V9WvSeNmvCci3h4R3wA217phSTMkrZa0pBCbJulJSYvy66TCsqmSOiQ9IumEQnx8jnVImtL9c8zMrLGqFY6/An4HzJF0jaTj6HlgpkquA8ZXiF8eEePy6w4ASQeTjm4OyetcKWlAHm3wm6QuTg4GzshtzcysSap1cvijiHg/cBBwL/AJYB9JV0k6vrcNR8QvgFqfMJ8AzIyIjRHxOGn42CPyqyMiHouI3wMzc1szM2uSXjs5jIgXIuLGiDgZGAEsYusuSMo6T9LifCprSI4NB1YU2nTmWE/xV5A0WdICSQvWrFmzHemZmVk1tfSOu0VErI2Ib0XEsdv4eVcBBwDjgFXAZTle6RRYVIlXym16RLRHRHtbW9s2pmdmZr2p6XbcvhIRT3VNS7oGuD3PdgIjC01HACvzdE9xMzNrglJHHNtL0rDC7HuBrjuuZgGnS9pV0hhgLDAPmA+MlTRG0i6kC+izGpmzmZltrW5HHJK+DxwD7CWpE/g8cEx+mDCA5cBHACJiqaRbgIeBTcC5EbE5b+c84E5gADAjIpbWK2czM+td3QpHRJxRIXxtlfaXAJdUiN8B3NGHqZmZ2XZo6KkqMzPr/1w4zMysFBcOMzMrxYXDzMxKceEwM7NSXDjMzKwUFw4zMyvFhcPMzEpx4TAzs1JcOMzMrBQXDjMzK8WFw8zMSnHhMDOzUlw4zMysFBcOMzMrpW6FQ9IMSaslLSnEhkqaLenR/D4kxyXpCkkdkhZLOqywzqTc/lFJk+qVr5mZ1aaeRxzXAeO7xaYAd0fEWODuPA9wImm42LHAZOAqSIWGNHLgkcARwOe7io2ZmTVH3QpHRPwCWNstPAG4Pk9fD5xSiN8Qyf3A4Dw++QnA7IhYGxHPArN5ZTEyM7MGavQ1jn0iYhVAft87x4cDKwrtOnOsp/grSJosaYGkBWvWrOnzxM3MLGmVi+OqEIsq8VcGI6ZHRHtEtLe1tfVpcmZm9keNLhxP5VNQ5PfVOd4JjCy0GwGsrBI3M7MmaXThmAV03Rk1CbitEJ+Y7646ClifT2XdCRwvaUi+KH58jpmZWZPsXK8NS/o+cAywl6RO0t1RlwK3SDobeAI4LTe/AzgJ6ABeBM4CiIi1kr4AzM/tLo6I7hfczcysgepWOCLijB4WHVehbQDn9rCdGcCMPkzNzMy2Q6tcHDczs37ChcPMzEpx4TAzs1JcOMzMrBQXDjMzK8WFw8zMSnHhMDOzUlw4zMysFBcOMzMrxYXDzMxKceEwM7NSXDjMzKwUFw4zMyvFhcPMzEpx4TAzs1KaUjgkLZf0kKRFkhbk2FBJsyU9mt+H5LgkXSGpQ9JiSYc1I2czM0uaecTxzogYFxHteX4KcHdEjAXuzvMAJwJj82sycFXDMzUzsy1a6VTVBOD6PH09cEohfkMk9wODJQ1rRoJmZta8whHAXZIWSpqcY/tExCqA/L53jg8HVhTW7cwxMzNrgrqNOd6LoyNipaS9gdmSfl2lrSrE4hWNUgGaDDBq1Ki+ydLMzF6hKUccEbEyv68GfgQcATzVdQoqv6/OzTuBkYXVRwArK2xzekS0R0R7W1tbPdM3M9uhNbxwSNpN0h5d08DxwBJgFjApN5sE3JanZwET891VRwHru05pmZlZ4zXjVNU+wI8kdX3+TRHxM0nzgVsknQ08AZyW298BnAR0AC8CZzU+5cpGT/lJTe2WX/ruOmdiZtY4DS8cEfEYcGiF+DPAcRXiAZzbgNTMzKwGrXQ7rpmZ9QMuHGZmVooLh5mZleLCYWZmpbhwmJlZKS4cZmZWiguHmZmV4sJhZmaluHCYmVkpzeodd4firknM7NXERxxmZlaKC4eZmZXiwmFmZqW4cJiZWSm+ON5CfBHdzPoDH3GYmVkp/eaIQ9J44J+AAcC3I+LSJqfUND4yMbNm6heFQ9IA4JvAXwCdwHxJsyLi4eZm9urhYmRmteoXhQM4AujIw84iaSYwAXDhqKLWYtDMbfZ1IXIBNKu//lI4hgMrCvOdwJHFBpImA5Pz7AZJj2zjZ+0FPL2N6zbSqyJPfbmBmVT/3FfF/mwR/SFHcJ6V7FdLo/5SOFQhFlvNREwHpm/3B0kLIqJ9e7dTb86zbznPvtMfcgTnuT36y11VncDIwvwIYGWTcjEz26H1l8IxHxgraYykXYDTgVlNzsnMbIfUL05VRcQmSecBd5Jux50REUvr9HHbfbqrQZxn33Kefac/5AjOc5spInpvZWZmlvWXU1VmZtYiXDjMzKwUF45M0nhJj0jqkDSl2fkUSVou6SFJiyQtyLGhkmZLejS/D2lCXjMkrZa0pBCrmJeSK/L+XSzpsCbnOU3Sk3mfLpJ0UmHZ1JznI5JOaGCeIyXNkbRM0lJJ5+d4S+3TKnm21D6VNEjSPEkP5jwvyvExkubm/XlzvuEGSbvm+Y68fHST87xO0uOF/Tkux5v2t7RFROzwL9IF9/8C9gd2AR4EDm52XoX8lgN7dYt9BZiSp6cAX25CXu8ADgOW9JYXcBLwU9IzOUcBc5uc5zTgUxXaHpz/++8KjMn/XwxoUJ7DgMPy9B7Ab3I+LbVPq+TZUvs075fd8/RAYG7eT7cAp+f41cBH8/THgKvz9OnAzQ3anz3leR1waoX2Tftb6nr5iCPZ0qVJRPwe6OrSpJVNAK7P09cDpzQ6gYj4BbC2W7invCYAN0RyPzBY0rAm5tmTCcDMiNgYEY8DHaT/P+ouIlZFxAN5+nlgGanXhJbap1Xy7ElT9mneLxvy7MD8CuBY4NYc774/u/bzrcBxkio9fNyoPHvStL+lLi4cSaUuTar9ITRaAHdJWpi7VgHYJyJWQfpDBvZuWnZb6ymvVtzH5+VD/RmFU30tkWc+TfJW0r8+W3afdssTWmyfShogaRGwGphNOtpZFxGbKuSyJc+8fD3w+mbkGRFd+/OSvD8vl7Rr9zyzhv93d+FIeu3SpMmOjojDgBOBcyW9o9kJbYNW28dXAQcA44BVwGU53vQ8Je0O/AC4ICKeq9a0QqxhuVbIs+X2aURsjohxpN4mjgDeVCWXlslT0puBqcBBwOHAUODCZufZxYUjaekuTSJiZX5fDfyI9AfwVNfhaX5f3bwMt9JTXi21jyPiqfzH+jJwDX88ddLUPCUNJP0Y3xgRP8zhltunlfJs1X2ac1sH3Eu6JjBYUtfDz8VctuSZl+9J7ac4+zrP8fmUYETERuA7tND+dOFIWrZLE0m7Sdqjaxo4HlhCym9SbjYJuK05Gb5CT3nNAibmO0KOAtZ3nX5phm7nhN9L2qeQ8jw932EzBhgLzGtQTgKuBZZFxNcKi1pqn/aUZ6vtU0ltkgbn6dcA7yJdj5kDnJqbdd+fXfv5VOCeyFejm5Dnrwv/WBDpOkxxfzb3b6nRV+Nb9UW6U+E3pHOgf9fsfAp57U+6I+VBYGlXbqRzr3cDj+b3oU3I7fukUxJ/IP0r6Oye8iIdXn8z79+HgPYm5/ndnMdi0h/isEL7v8t5PgKc2MA830465bAYWJRfJ7XaPq2SZ0vtU+AtwK9yPkuAz+X4/qTC1QH8C7Brjg/K8x15+f5NzvOevD+XAN/jj3c58Wl2AAAEnklEQVReNe1vqevlLkfMzKwUn6oyM7NSXDjMzKwUFw4zMyvFhcPMzEpx4TAzs1JcOKylSApJlxXmPyVpWh9t+zpJp/becrs/5zSlnmPndIuPVqGH3irrf7Z+2W2Vy18X5tslXVHvz7VXBxcOazUbgf8jaa9mJ1IkaUCJ5mcDH4uId27jx5UuHCXzAxgNbCkcEbEgIj5e9nNtx+TCYa1mE2mM5U90X9D9iEHShvx+jKT7JN0i6TeSLpV0Zh7j4CFJBxQ28y5J/5bbnZzXHyDpq5Lm5w7lPlLY7hxJN5EetOqezxl5+0skfTnHPkd6QO5qSV/t6UtK+pCkH0r6mdK4EF/J8UuB1yiNv3Bjjn0gf5dFkr7VVSQkbZB0saS5wJ9K+lz+DkskTc9PHCPpQEk/Vxrv4YG8Py4F/ixv8xP5u96e2w+V9OO8L+6X9JYcn6bUeeG9kh6T9PEc303ST/L2l0h6f+//ma1fa/QTh375Ve0FbABeRxqDZE/gU8C0vOw6CuMTABvy+zHAOtI4EbsCTwIX5WXnA18vrP8z0j+YxpKeIh8ETAb+PrfZFVhAGjfiGOAFYEyFPPcFngDagJ1JT/mekpfdS4WneUn/yl+Spz8EPJa/4yDgt8DI4vfK028C/hUYmOevBCbm6QDeV2g7tDD9XeA9eXou8N48PQh4bf5utxfab5kHvgF8Pk8fCyzK09OA/8j7aC/gGVIX4H8FXFPY1p7N/v/Ir/q+fMRhLSdST6s3AGVOncyP1CncRlJXDHfl+EOkH+wut0TEyxHxKOmH+yBS/18Tlbq1nkvq4mNsbj8v0hgS3R0O3BsRayJ1wX0jacCoMu6OiPUR8RLwMLBfhTbHAW8D5uf8jiN1mQGwmdTRYJd3Ko1c9xDpB/8QpX7OhkfEjwAi4qWIeLGXvN5OKjxExD3A6yXtmZf9JNK4Gk+TOlvch7SP3yXpy5L+LCLWl9oL1u/s3HsTs6b4OvAAqVfQLpvIp1fzaZhdCss2FqZfLsy/zNb/n3fvYydIff/8v4i4s7hA0jGkI45K+mKAn2LOm6n89yjg+oiYWmHZSxGxGdLwo6SjkfaIWJFvKBi0jXlW67b7FTlHxG8kvY3UX9WXJN0VERdvw+daP+EjDmtJEbGWNMTn2YXwctK/viGNgjZwGzZ9mqSd8nn+/Umd7t0JfFSpq3AkvUGpJ+Jq5gJ/LmmvfM3hDOC+bcinkj905ULq1PBUSXvn3IZKqnRkMii/P600TsapsOXorVPSKXn9XSW9FnieNOxrJb8AzsztjwGejirjgkjaF3gxIr4H/CNpmF57FfMRh7Wyy4DzCvPXALdJmkf6Qe3paKCaR0g/8PsA50TES5K+TTqd9UA+kllDL0PxRsQqSVNJXXQLuCMi+qpr++nAYkkPRMSZkv6eNALkTqQefs8lXRMp5rNO0jWk00bLSUMFdPkg8C1JF+f1TyP1xLpJ0oOkaz+/KrSfBnxH0mLgRf7Y1XhP/gT4qqSX8/Y/WvobW7/i3nHNzKwUn6oyM7NSXDjMzKwUFw4zMyvFhcPMzEpx4TAzs1JcOMzMrBQXDjMzK+V/ACQnNlw0DcumAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "x = df_n['article_counts']\n",
    "plt.hist(x, bins=30);\n",
    "plt.xlabel('Number of Interactions');\n",
    "plt.ylabel('Article Counts')\n",
    "plt.title('User-Article Interactions');\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "article_counts    3.0\n",
       "dtype: float64"
      ]
     },
     "execution_count": 75,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_n.median()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "50% of individuals interact with 3 number of articles or fewer.\n",
      "The maximum number of user-article interactions by any 1 user is 364.\n"
     ]
    }
   ],
   "source": [
    "# Fill in the median and maximum number of user_article interactios below\n",
    "\n",
    "median_val = 3  # 50% of individuals interact with ____ number of articles or fewer.\n",
    "print(\"50% of individuals interact with {} number of articles or fewer.\".format(median_val))\n",
    "max_views_by_user = df_n['article_counts'].max()  # The maximum number of user-article interactions by any 1 user is ______.\n",
    "print(\"The maximum number of user-article interactions by any 1 user is {}.\".format(max_views_by_user))\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`2.` Explore and remove duplicate articles from the **df_content** dataframe.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 77,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<class 'pandas.core.frame.DataFrame'>\n",
      "RangeIndex: 1056 entries, 0 to 1055\n",
      "Data columns (total 5 columns):\n",
      "doc_body           1042 non-null object\n",
      "doc_description    1053 non-null object\n",
      "doc_full_name      1056 non-null object\n",
      "doc_status         1056 non-null object\n",
      "article_id         1056 non-null int64\n",
      "dtypes: int64(1), object(4)\n",
      "memory usage: 41.3+ KB\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>doc_body</th>\n",
       "      <th>doc_description</th>\n",
       "      <th>doc_full_name</th>\n",
       "      <th>doc_status</th>\n",
       "      <th>article_id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Skip navigation Sign in SearchLoading...\\r\\n\\r...</td>\n",
       "      <td>Detect bad readings in real time using Python ...</td>\n",
       "      <td>Detect Malfunctioning IoT Sensors with Streami...</td>\n",
       "      <td>Live</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>No Free Hunch Navigation * kaggle.com\\r\\n\\r\\n ...</td>\n",
       "      <td>See the forest, see the trees. Here lies the c...</td>\n",
       "      <td>Communicating data science: A guide to present...</td>\n",
       "      <td>Live</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>☰ * Login\\r\\n * Sign Up\\r\\n\\r\\n * Learning Pat...</td>\n",
       "      <td>Here’s this week’s news in Data Science and Bi...</td>\n",
       "      <td>This Week in Data Science (April 18, 2017)</td>\n",
       "      <td>Live</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>DATALAYER: HIGH THROUGHPUT, LOW LATENCY AT SCA...</td>\n",
       "      <td>Learn how distributed DBs solve the problem of...</td>\n",
       "      <td>DataLayer Conference: Boost the performance of...</td>\n",
       "      <td>Live</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Skip navigation Sign in SearchLoading...\\r\\n\\r...</td>\n",
       "      <td>This video demonstrates the power of IBM DataS...</td>\n",
       "      <td>Analyze NY Restaurant data using Spark in DSX</td>\n",
       "      <td>Live</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            doc_body  \\\n",
       "0  Skip navigation Sign in SearchLoading...\\r\\n\\r...   \n",
       "1  No Free Hunch Navigation * kaggle.com\\r\\n\\r\\n ...   \n",
       "2  ☰ * Login\\r\\n * Sign Up\\r\\n\\r\\n * Learning Pat...   \n",
       "3  DATALAYER: HIGH THROUGHPUT, LOW LATENCY AT SCA...   \n",
       "4  Skip navigation Sign in SearchLoading...\\r\\n\\r...   \n",
       "\n",
       "                                     doc_description  \\\n",
       "0  Detect bad readings in real time using Python ...   \n",
       "1  See the forest, see the trees. Here lies the c...   \n",
       "2  Here’s this week’s news in Data Science and Bi...   \n",
       "3  Learn how distributed DBs solve the problem of...   \n",
       "4  This video demonstrates the power of IBM DataS...   \n",
       "\n",
       "                                       doc_full_name doc_status  article_id  \n",
       "0  Detect Malfunctioning IoT Sensors with Streami...       Live           0  \n",
       "1  Communicating data science: A guide to present...       Live           1  \n",
       "2         This Week in Data Science (April 18, 2017)       Live           2  \n",
       "3  DataLayer Conference: Boost the performance of...       Live           3  \n",
       "4      Analyze NY Restaurant data using Spark in DSX       Live           4  "
      ]
     },
     "execution_count": 77,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_content.info()\n",
    "df_content.head()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                              doc_body  \\\n",
      "366  GO SERVERLESS WITH APEX AND COMPOSE'S MONGODB\\...   \n",
      "\n",
      "                                       doc_description  \\\n",
      "366  While AWS's Lambda service treats each functio...   \n",
      "\n",
      "                                     doc_full_name doc_status  article_id  \n",
      "366  Go Serverless with Apex and Compose's MongoDB       Live         365  \n",
      "                                              doc_body  \\\n",
      "694  Skip to content * Unix\\r\\n * R\\r\\n\\r\\n15 PAGE ...   \n",
      "\n",
      "                                       doc_description  \\\n",
      "694  For Beginners in R, here is a 15 page example ...   \n",
      "\n",
      "              doc_full_name doc_status  article_id  \n",
      "694  15 Page Tutorial for R       Live         692  \n",
      "                                              doc_body  \\\n",
      "764  Elvis Dohmatob Home Blog Publications Photos C...   \n",
      "\n",
      "                                       doc_description  \\\n",
      "764  In this post, I’ll demo variational auto-encod...   \n",
      "\n",
      "                                         doc_full_name doc_status  article_id  \n",
      "764  Variational auto-encoder for \"Frey faces\" usin...       Live         761  \n",
      "                                              doc_body  \\\n",
      "975  Compose The Compose logo Articles Sign in Free...   \n",
      "\n",
      "                                       doc_description  \\\n",
      "975  In this third article in our the series, we'll...   \n",
      "\n",
      "                           doc_full_name doc_status  article_id  \n",
      "975  Launching RESTHeart into Production       Live         970  \n",
      "                                              doc_body  \\\n",
      "976  Cloudant allows custom Javascript to be run se...   \n",
      "\n",
      "                                       doc_description  \\\n",
      "976  Cloudant allows custom Javascript to be run se...   \n",
      "\n",
      "                               doc_full_name doc_status  article_id  \n",
      "976  Defensive coding in Map/Index functions       Live         971  \n"
     ]
    }
   ],
   "source": [
    "# Find and explore duplicate articles\n",
    "df_content.duplicated(['article_id']).sum()\n",
    "# Checking if there are any article ids which are duplicates\n",
    "df_cont_dupl = np.where(df_content.duplicated(['article_id']) == True)\n",
    "df_cont_dupl[0]\n",
    "for id in df_cont_dupl[0]:\n",
    "    print(df_content[df_content['article_id'] == id])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1051, 5)"
      ]
     },
     "execution_count": 79,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Remove any rows that have the same article_id and keeyp keep the first only\n",
    "df_cont_dropDupl = df_content.drop_duplicates(subset='article_id', keep='first')\n",
    "df_cont_dropDupl.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`3.` Use the cells below to find:\n",
    "\n",
    "**a.** The number of unique articles that have an interaction with a user.  \n",
    "**b.** The number of unique articles in the dataset (whether they have any interactions or not).<br>\n",
    "**c.** The number of unique users in the dataset. (excluding null values) <br>\n",
    "**d.** The number of user-article interactions in the dataset."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>article_id</th>\n",
       "      <th>title</th>\n",
       "      <th>email</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1430.0</td>\n",
       "      <td>using pixiedust for fast, flexible, and easier...</td>\n",
       "      <td>ef5f11f77ba020cd36e1105a00ab868bbdbf7fe7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1314.0</td>\n",
       "      <td>healthcare python streaming application demo</td>\n",
       "      <td>083cbdfa93c8444beaa4c5f5e0f5f9198e4f9e0b</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1429.0</td>\n",
       "      <td>use deep learning for image classification</td>\n",
       "      <td>b96a4f2e92d8572034b1e9b28f9ac673765cd074</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1338.0</td>\n",
       "      <td>ml optimization using cognitive assistant</td>\n",
       "      <td>06485706b34a5c9bf2a0ecdac41daf7e7654ceb7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1276.0</td>\n",
       "      <td>deploy your python model as a restful api</td>\n",
       "      <td>f01220c46fc92c6e6b161b1849de11faacd7ccb2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   article_id                                              title  \\\n",
       "0      1430.0  using pixiedust for fast, flexible, and easier...   \n",
       "1      1314.0       healthcare python streaming application demo   \n",
       "2      1429.0         use deep learning for image classification   \n",
       "3      1338.0          ml optimization using cognitive assistant   \n",
       "4      1276.0          deploy your python model as a restful api   \n",
       "\n",
       "                                      email  \n",
       "0  ef5f11f77ba020cd36e1105a00ab868bbdbf7fe7  \n",
       "1  083cbdfa93c8444beaa4c5f5e0f5f9198e4f9e0b  \n",
       "2  b96a4f2e92d8572034b1e9b28f9ac673765cd074  \n",
       "3  06485706b34a5c9bf2a0ecdac41daf7e7654ceb7  \n",
       "4  f01220c46fc92c6e6b161b1849de11faacd7ccb2  "
      ]
     },
     "execution_count": 80,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 81,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "article_id     714\n",
       "title          714\n",
       "email         5148\n",
       "dtype: int64"
      ]
     },
     "execution_count": 81,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.nunique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 82,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>title</th>\n",
       "      <th>email</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>article_id</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1429.0</th>\n",
       "      <td>937</td>\n",
       "      <td>937</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1330.0</th>\n",
       "      <td>927</td>\n",
       "      <td>927</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1431.0</th>\n",
       "      <td>671</td>\n",
       "      <td>671</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1427.0</th>\n",
       "      <td>643</td>\n",
       "      <td>643</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1364.0</th>\n",
       "      <td>627</td>\n",
       "      <td>627</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "            title  email\n",
       "article_id              \n",
       "1429.0        937    937\n",
       "1330.0        927    927\n",
       "1431.0        671    671\n",
       "1427.0        643    643\n",
       "1364.0        627    627"
      ]
     },
     "execution_count": 82,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(df.groupby(['article_id']).count()).sort_values('title', ascending=False).head()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>article_id</th>\n",
       "      <th>title</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>email</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0000b6387a0366322d7fbfc6434af145adf7fed1</th>\n",
       "      <td>13</td>\n",
       "      <td>13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>001055fc0bb67f71e8fa17002342b256a30254cd</th>\n",
       "      <td>4</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>00148e4911c7e04eeff8def7bbbdaf1c59c2c621</th>\n",
       "      <td>3</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>001a852ecbd6cc12ab77a785efa137b2646505fe</th>\n",
       "      <td>6</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>001fc95b90da5c3cb12c501d201a915e4f093290</th>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                          article_id  title\n",
       "email                                                      \n",
       "0000b6387a0366322d7fbfc6434af145adf7fed1          13     13\n",
       "001055fc0bb67f71e8fa17002342b256a30254cd           4      4\n",
       "00148e4911c7e04eeff8def7bbbdaf1c59c2c621           3      3\n",
       "001a852ecbd6cc12ab77a785efa137b2646505fe           6      6\n",
       "001fc95b90da5c3cb12c501d201a915e4f093290           2      2"
      ]
     },
     "execution_count": 83,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.groupby(['email']).count().head()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 84,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(45993, 3)"
      ]
     },
     "execution_count": 84,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.shape\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 85,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1051, 5)"
      ]
     },
     "execution_count": 85,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_cont_dropDupl.shape\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The number of unique articles that have at least one interaction is 714.\n",
      "The number of unique articles on the IBM platform is 1051.\n",
      "The number of unique users is 5148.\n",
      "The number of user-article interactions is 45993.\n"
     ]
    }
   ],
   "source": [
    "unique_articles = df.nunique()['article_id'] # The number of unique articles that have at least one interaction\n",
    "total_articles = df_cont_dropDupl.shape[0] # The number of unique articles on the IBM platform\n",
    "unique_users = df.nunique()['email'] # The number of unique users\n",
    "user_article_interactions = df.shape[0] # The number of user-article interactions\n",
    "print(\"The number of unique articles that have at least one interaction is {}.\".format(unique_articles))\n",
    "print(\"The number of unique articles on the IBM platform is {}.\".format(total_articles))\n",
    "print(\"The number of unique users is {}.\".format(unique_users))\n",
    "print(\"The number of user-article interactions is {}.\".format(user_article_interactions))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`4.` Use the cells below to find the most viewed **article_id**, as well as how often it was viewed.  After talking to the company leaders, the `email_mapper` function was deemed a reasonable way to map users to ids.  There were a small number of null values, and it was found that all of these null values likely belonged to a single user (which is how they are stored using the function below)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 87,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>title</th>\n",
       "      <th>email</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>article_id</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1429.0</th>\n",
       "      <td>937</td>\n",
       "      <td>937</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "            title  email\n",
       "article_id              \n",
       "1429.0        937    937"
      ]
     },
     "execution_count": 87,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_view_most = df.groupby(['article_id']).count()\n",
    "df_view_most.loc[df_view_most['email'] == df_view_most['email'].max()]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 88,
   "metadata": {},
   "outputs": [],
   "source": [
    "most_viewed_article_id = \"1429.0\" # The most viewed article in the dataset as a string with one value following the decimal \n",
    "max_views = 937 # The most viewed article in the dataset was viewed how many times?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 89,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>article_id</th>\n",
       "      <th>title</th>\n",
       "      <th>user_id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1430.0</td>\n",
       "      <td>using pixiedust for fast, flexible, and easier...</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1314.0</td>\n",
       "      <td>healthcare python streaming application demo</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1429.0</td>\n",
       "      <td>use deep learning for image classification</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1338.0</td>\n",
       "      <td>ml optimization using cognitive assistant</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1276.0</td>\n",
       "      <td>deploy your python model as a restful api</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   article_id                                              title  user_id\n",
       "0      1430.0  using pixiedust for fast, flexible, and easier...        1\n",
       "1      1314.0       healthcare python streaming application demo        2\n",
       "2      1429.0         use deep learning for image classification        3\n",
       "3      1338.0          ml optimization using cognitive assistant        4\n",
       "4      1276.0          deploy your python model as a restful api        5"
      ]
     },
     "execution_count": 89,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "## No need to change the code here - this will be helpful for later parts of the notebook\n",
    "# Run this cell to map the user email to a user_id column and remove the email column\n",
    "\n",
    "def email_mapper():\n",
    "    coded_dict = dict()\n",
    "    cter = 1\n",
    "    email_encoded = []\n",
    "    \n",
    "    for val in df['email']:\n",
    "        if val not in coded_dict:\n",
    "            coded_dict[val] = cter\n",
    "            cter+=1\n",
    "        \n",
    "        email_encoded.append(coded_dict[val])\n",
    "    return email_encoded\n",
    "\n",
    "email_encoded = email_mapper()\n",
    "del df['email']\n",
    "df['user_id'] = email_encoded\n",
    "\n",
    "# show header\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "It looks like you have everything right here! Nice job!\n"
     ]
    }
   ],
   "source": [
    "## If you stored all your results in the variable names above, \n",
    "## you shouldn't need to change anything in this cell\n",
    "\n",
    "sol_1_dict = {\n",
    "    '`50% of individuals have _____ or fewer interactions.`': median_val,\n",
    "    '`The total number of user-article interactions in the dataset is ______.`': user_article_interactions,\n",
    "    '`The maximum number of user-article interactions by any 1 user is ______.`': max_views_by_user,\n",
    "    '`The most viewed article in the dataset was viewed _____ times.`': max_views,\n",
    "    '`The article_id of the most viewed article is ______.`': most_viewed_article_id,\n",
    "    '`The number of unique articles that have at least 1 rating ______.`': unique_articles,\n",
    "    '`The number of unique users in the dataset is ______`': unique_users,\n",
    "    '`The number of unique articles on the IBM platform`': total_articles\n",
    "}\n",
    "\n",
    "# Test your dictionary against the solution\n",
    "t.sol_1_test(sol_1_dict)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### <a class=\"anchor\" id=\"Rank\">Part II: Rank-Based Recommendations</a>\n",
    "\n",
    "Unlike in the earlier lessons, we don't actually have ratings for whether a user liked an article or not.  We only know that a user has interacted with an article.  In these cases, the popularity of an article can really only be based on how often an article was interacted with.\n",
    "\n",
    "`1.` Fill in the function below to return the **n** top articles ordered with most interactions as the top. Test your function using the tests below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 91,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>article_id</th>\n",
       "      <th>title</th>\n",
       "      <th>user_id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1430.0</td>\n",
       "      <td>using pixiedust for fast, flexible, and easier...</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1314.0</td>\n",
       "      <td>healthcare python streaming application demo</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1429.0</td>\n",
       "      <td>use deep learning for image classification</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1338.0</td>\n",
       "      <td>ml optimization using cognitive assistant</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1276.0</td>\n",
       "      <td>deploy your python model as a restful api</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   article_id                                              title  user_id\n",
       "0      1430.0  using pixiedust for fast, flexible, and easier...        1\n",
       "1      1314.0       healthcare python streaming application demo        2\n",
       "2      1429.0         use deep learning for image classification        3\n",
       "3      1338.0          ml optimization using cognitive assistant        4\n",
       "4      1276.0          deploy your python model as a restful api        5"
      ]
     },
     "execution_count": 91,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 92,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_top_articles(n, df=df):\n",
    "    '''\n",
    "    INPUT:\n",
    "    n - (int) the number of top articles to return\n",
    "    df - (pandas dataframe) df as defined at the top of the notebook \n",
    "    \n",
    "    OUTPUT:\n",
    "    top_articles - (list) A list of the top 'n' article titles \n",
    "    \n",
    "    '''\n",
    "    # Your code here\n",
    "    top_articles = df.groupby(['article_id', 'title']).size()\\\n",
    "                        .reset_index(name='counts').sort_values('counts', ascending=False)[:n].title.tolist()\n",
    "    return top_articles # Return the top article titles from df (not df_content)\n",
    "\n",
    "def get_top_article_ids(n, df=df):\n",
    "    '''\n",
    "    INPUT:\n",
    "    n - (int) the number of top articles to return\n",
    "    df - (pandas dataframe) df as defined at the top of the notebook \n",
    "    \n",
    "    OUTPUT:\n",
    "    top_articles - (list) A list of the top 'n' article titles \n",
    "    \n",
    "    '''\n",
    "    # Your code here\n",
    "    top_articles = df.groupby(\"article_id\").count()[\"title\"].sort_values(ascending=False).index[:n].astype('str')\n",
    "    return top_articles # Return the top article ids"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 93,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['use deep learning for image classification', 'insights from new york car accident reports', 'visualize car data with brunel', 'use xgboost, scikit-learn & ibm watson machine learning apis', 'predicting churn with the spss random tree algorithm', 'healthcare python streaming application demo', 'finding optimal locations of new store using decision optimization', 'apache spark lab, part 1: basic concepts', 'analyze energy consumption in buildings', 'gosales transactions for logistic regression model']\n",
      "Index(['1429.0', '1330.0', '1431.0', '1427.0', '1364.0', '1314.0', '1293.0',\n",
      "       '1170.0', '1162.0', '1304.0'],\n",
      "      dtype='object', name='article_id')\n"
     ]
    }
   ],
   "source": [
    "print(get_top_articles(10))\n",
    "print(get_top_article_ids(10))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 94,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Your top_5 looks like the solution list! Nice job.\n",
      "Your top_10 looks like the solution list! Nice job.\n",
      "Your top_20 looks like the solution list! Nice job.\n"
     ]
    }
   ],
   "source": [
    "# Test your function by returning the top 5, 10, and 20 articles\n",
    "top_5 = get_top_articles(5)\n",
    "top_10 = get_top_articles(10)\n",
    "top_20 = get_top_articles(20)\n",
    "\n",
    "# Test each of your three lists from above\n",
    "t.sol_2_test(get_top_articles)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### <a class=\"anchor\" id=\"User-User\">Part III: User-User Based Collaborative Filtering</a>\n",
    "\n",
    "\n",
    "`1.` Use the function below to reformat the **df** dataframe to be shaped with users as the rows and articles as the columns.  \n",
    "\n",
    "* Each **user** should only appear in each **row** once.\n",
    "\n",
    "\n",
    "* Each **article** should only show up in one **column**.  \n",
    "\n",
    "\n",
    "* **If a user has interacted with an article, then place a 1 where the user-row meets for that article-column**.  It does not matter how many times a user has interacted with the article, all entries where a user has interacted with an article should be a 1.  \n",
    "\n",
    "\n",
    "* **If a user has not interacted with an item, then place a zero where the user-row meets for that article-column**. \n",
    "\n",
    "Use the tests to make sure the basic structure of your matrix matches what is expected by the solution."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 95,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create the user-article matrix with 1's and 0's\n",
    "def create_user_item_matrix(df):\n",
    "    '''\n",
    "    INPUT:\n",
    "    df - pandas dataframe with article_id, title, user_id columns\n",
    "    \n",
    "    OUTPUT:\n",
    "    user_item - user item matrix \n",
    "    \n",
    "    Description:\n",
    "    Return a matrix with user ids as rows and article ids on the columns with 1 values where a user interacted with \n",
    "    an article and a 0 otherwise\n",
    "    '''\n",
    "    # Fill in the function here\n",
    "    # unstack the user-item interaction dataframe\n",
    "    user_item = df.drop_duplicates().groupby(['user_id', 'article_id']).size().unstack()\n",
    "    \n",
    "    # fill missing values with 0\n",
    "    user_item = user_item.fillna(0)\n",
    "    \n",
    "    # convert int\n",
    "    user_item = user_item.astype('int')\n",
    "    \n",
    "    return user_item # return the user_item matrix \n",
    "\n",
    "user_item = create_user_item_matrix(df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 96,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "You have passed our quick tests!  Please proceed!\n"
     ]
    }
   ],
   "source": [
    "## Tests: You should just need to run this cell.  Don't change the code.\n",
    "assert user_item.shape[0] == 5149, \"Oops!  The number of users in the user-article matrix doesn't look right.\"\n",
    "assert user_item.shape[1] == 714, \"Oops!  The number of articles in the user-article matrix doesn't look right.\"\n",
    "assert user_item.sum(axis=1)[1] == 36, \"Oops!  The number of articles seen by user 1 doesn't look right.\"\n",
    "print(\"You have passed our quick tests!  Please proceed!\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`2.` Complete the function below which should take a user_id and provide an ordered list of the most similar users to that user (from most similar to least similar).  The returned result should not contain the provided user_id, as we know that each user is similar to him/herself. Because the results for each user here are binary, it (perhaps) makes sense to compute similarity as the dot product of two users. \n",
    "\n",
    "Use the tests to test your function."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 97,
   "metadata": {},
   "outputs": [],
   "source": [
    "def find_similar_users(user_id, user_item=user_item):\n",
    "    '''\n",
    "    INPUT:\n",
    "    user_id - (int) a user_id\n",
    "    user_item - (pandas dataframe) matrix of users by articles: \n",
    "                1's when a user has interacted with an article, 0 otherwise\n",
    "    \n",
    "    OUTPUT:\n",
    "    similar_users - (list) an ordered list where the closest users (largest dot product users)\n",
    "                    are listed first\n",
    "    \n",
    "    Description:\n",
    "    Computes the similarity of every pair of users based on the dot product\n",
    "    Returns an ordered\n",
    "    \n",
    "    '''\n",
    "    # Compute similarity of each user to the provided user\n",
    "    user_similr = user_item.loc[user_id,:].dot(user_item.T)\n",
    "\n",
    "    # Sort by similarity\n",
    "    user_similr  = user_similr.sort_values(ascending = False)\n",
    "\n",
    "    # Create list of just the ids\n",
    "    # Remove the own user's id\n",
    "    most_similar_users = user_similr.loc[~(user_similr.index ==  user_id)].index.values.tolist()\n",
    "       \n",
    "    return most_similar_users # return a list of the users in order from most to least similar\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 98,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The 10 most similar users to user 1 are: [3933, 23, 3782, 203, 4459, 131, 3870, 46, 4201, 5041]\n",
      "The 5 most similar users to user 3933 are: [1, 23, 3782, 4459, 203]\n",
      "The 3 most similar users to user 46 are: [4201, 23, 3782]\n"
     ]
    }
   ],
   "source": [
    "# Do a spot check of your function\n",
    "print(\"The 10 most similar users to user 1 are: {}\".format(find_similar_users(1)[:10]))\n",
    "print(\"The 5 most similar users to user 3933 are: {}\".format(find_similar_users(3933)[:5]))\n",
    "print(\"The 3 most similar users to user 46 are: {}\".format(find_similar_users(46)[:3]))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`3.` Now that you have a function that provides the most similar users to each user, you will want to use these users to find articles you can recommend.  Complete the functions below to return the articles you would recommend to each user. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 99,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_article_names(article_ids, df=df):\n",
    "    '''\n",
    "    INPUT:\n",
    "    article_ids - (list) a list of article ids\n",
    "    df - (pandas dataframe) df as defined at the top of the notebook\n",
    "    \n",
    "    OUTPUT:\n",
    "    article_names - (list) a list of article names associated with the list of article ids \n",
    "                    (this is identified by the title column)\n",
    "    '''\n",
    "    # Your code here\n",
    "    article_names = []\n",
    "    # Select articles with the same article_id and drop duplicates \n",
    "    article_names = df[df['article_id'].isin(article_ids)]['title'].drop_duplicates().values.tolist()\n",
    "\n",
    "    return article_names # Return the article names associated with list of article ids\n",
    "\n",
    "\n",
    "def get_user_articles(user_id, user_item=user_item):\n",
    "    '''\n",
    "    INPUT:\n",
    "    user_id - (int) a user id\n",
    "    user_item - (pandas dataframe) matrix of users by articles: \n",
    "                1's when a user has interacted with an article, 0 otherwise\n",
    "    \n",
    "    OUTPUT:\n",
    "    article_ids - (list) a list of the article ids seen by the user\n",
    "    article_names - (list) a list of article names associated with the list of article ids \n",
    "                    (this is identified by the doc_full_name column in df_content)\n",
    "    \n",
    "    Description:\n",
    "    Provides a list of the article_ids and article titles that have been seen by a user\n",
    "    '''\n",
    "    # Your code here\n",
    "    # Get all articles for this user id\n",
    "    user_idx = user_item.loc[user_id,:]\n",
    "    # Get articles user interacted with\n",
    "    article_ids = user_idx[user_idx == 1].index.values.astype('str').tolist()\n",
    "    # Get article names \n",
    "    article_names = get_article_names(article_ids)\n",
    "    \n",
    "    return article_ids, article_names # return the ids and names\n",
    "\n",
    "\n",
    "def user_user_recs(user_id, m=10):\n",
    "    '''\n",
    "    INPUT:\n",
    "    user_id - (int) a user id\n",
    "    m - (int) the number of recommendations you want for the user\n",
    "    \n",
    "    OUTPUT:\n",
    "    recs - (list) a list of recommendations for the user\n",
    "    \n",
    "    Description:\n",
    "    Loops through the users based on closeness to the input user_id\n",
    "    For each user - finds articles the user hasn't seen before and provides them as recs\n",
    "    Does this until m recommendations are found\n",
    "    \n",
    "    Notes:\n",
    "    Users who are the same closeness are chosen arbitrarily as the 'next' user\n",
    "    \n",
    "    For the user where the number of recommended articles starts below m \n",
    "    and ends exceeding m, the last items are chosen arbitrarily\n",
    "    \n",
    "    '''\n",
    "    # Your code here\n",
    "    # Get most similar users \n",
    "    most_similar_users = find_similar_users(user_id)\n",
    "    # Get article ids \n",
    "    user_article_ids = set(get_user_articles(user_id)[0])\n",
    "    recs = []\n",
    "    # Create recommendations for this user\n",
    "    for user_neighb in most_similar_users:\n",
    "        neighb_article_ids = set(get_user_articles(user_neighb)[0])\n",
    "        recs += list(set(neighb_article_ids) - set(user_article_ids))\n",
    "        if len(recs) > m:\n",
    "            break\n",
    "    recs = recs[:m]\n",
    "    \n",
    "    return recs # return your recommendations for this user_id    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 100,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['visualize data with the matplotlib library',\n",
       " 'predicting churn with the spss random tree algorithm',\n",
       " 'apache spark lab, part 3: machine learning',\n",
       " 'using github for project control in dsx',\n",
       " 'analyzing data by using the sparkling.data library features',\n",
       " '502    forgetting the past to learn the future: long ...\\nName: title, dtype: object',\n",
       " 'spark 2.1 and job monitoring available in dsx',\n",
       " 'higher-order logistic regression for large datasets',\n",
       " 'recommender systems: approaches & algorithms',\n",
       " 'using bigdl in dsx for deep learning on spark']"
      ]
     },
     "execution_count": 100,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Check Results\n",
    "get_article_names(user_user_recs(1, 10)) # Return 10 recommendations for user 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 101,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "If this is all you see, you passed all of our tests!  Nice job!\n"
     ]
    }
   ],
   "source": [
    "# Test your functions here - No need to change this code - just run this cell\n",
    "assert set(get_article_names(['1024.0', '1176.0', '1305.0', '1314.0', '1422.0', '1427.0'])) == set(['using deep learning to reconstruct high-resolution audio', 'build a python app on the streaming analytics service', 'gosales transactions for naive bayes model', 'healthcare python streaming application demo', 'use r dataframes & ibm watson natural language understanding', 'use xgboost, scikit-learn & ibm watson machine learning apis']), \"Oops! Your the get_article_names function doesn't work quite how we expect.\"\n",
    "assert set(get_article_names(['1320.0', '232.0', '844.0'])) == set(['housing (2015): united states demographic measures','self-service data preparation with ibm data refinery','use the cloudant-spark connector in python notebook']), \"Oops! Your the get_article_names function doesn't work quite how we expect.\"\n",
    "assert set(get_user_articles(20)[0]) == set(['1320.0', '232.0', '844.0'])\n",
    "assert set(get_user_articles(20)[1]) == set(['housing (2015): united states demographic measures', 'self-service data preparation with ibm data refinery','use the cloudant-spark connector in python notebook'])\n",
    "assert set(get_user_articles(2)[0]) == set(['1024.0', '1176.0', '1305.0', '1314.0', '1422.0', '1427.0'])\n",
    "assert set(get_user_articles(2)[1]) == set(['using deep learning to reconstruct high-resolution audio', 'build a python app on the streaming analytics service', 'gosales transactions for naive bayes model', 'healthcare python streaming application demo', 'use r dataframes & ibm watson natural language understanding', 'use xgboost, scikit-learn & ibm watson machine learning apis'])\n",
    "print(\"If this is all you see, you passed all of our tests!  Nice job!\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`4.` Now we are going to improve the consistency of the **user_user_recs** function from above.  \n",
    "\n",
    "* Instead of arbitrarily choosing when we obtain users who are all the same closeness to a given user - choose the users that have the most total article interactions before choosing those with fewer article interactions.\n",
    "\n",
    "\n",
    "* Instead of arbitrarily choosing articles from the user where the number of recommended articles starts below m and ends exceeding m, choose articles with the articles with the most total interactions before choosing those with fewer total interactions. This ranking should be  what would be obtained from the **top_articles** function you wrote earlier."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 102,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_top_sorted_users(user_id, df=df, user_item=user_item):\n",
    "    '''\n",
    "    INPUT:\n",
    "    user_id - (int)\n",
    "    df - (pandas dataframe) df as defined at the top of the notebook \n",
    "    user_item - (pandas dataframe) matrix of users by articles: \n",
    "            1's when a user has interacted with an article, 0 otherwise\n",
    "    \n",
    "            \n",
    "    OUTPUT:\n",
    "    neighbors_df - (pandas dataframe) a dataframe with:\n",
    "                    neighbor_id - is a neighbor user_id\n",
    "                    similarity - measure of the similarity of each user to the provided user_id\n",
    "                    num_interactions - the number of articles viewed by the user - if a u\n",
    "                    \n",
    "    Other Details - sort the neighbors_df by the similarity and then by number of interactions where \n",
    "                    highest of each is higher in the dataframe\n",
    "     \n",
    "    '''\n",
    "    # Your code here\n",
    "    # Column names\n",
    "    colName = ['neighbor_id', 'similarity', 'num_interactions'] \n",
    "    # Create dataframe to hold top users\n",
    "    neighbors_df  = pd.DataFrame(columns= colName) \n",
    "    \n",
    "    # Populate the dataframe\n",
    "    for id in user_item.index.values:\n",
    "        if id != user_id:\n",
    "            neighbor_id = id\n",
    "            \n",
    "            # Get user to user similarity\n",
    "            similarity = user_item[user_item.index == user_id].dot(user_item.loc[id].T).values[0]\n",
    "            \n",
    "            # Get number of interactions for user ot article\n",
    "            num_interactions = user_item.loc[id].values.sum()\n",
    "            neighbors_df.loc[neighbor_id] = [neighbor_id, similarity, num_interactions]\n",
    "    \n",
    "    neighbors_df['similarity'] = neighbors_df['similarity'].astype('int')\n",
    "    neighbors_df['neighbor_id'] = neighbors_df['neighbor_id'].astype('int')\n",
    "    neighbors_df = neighbors_df.sort_values(by = ['similarity', 'neighbor_id'], ascending = [False, True])\n",
    "    \n",
    "    return neighbors_df # Return the dataframe specified in the doc_string\n",
    "\n",
    "\n",
    "def user_user_recs_part2(user_id, m=10):\n",
    "    '''\n",
    "    INPUT:\n",
    "    user_id - (int) a user id\n",
    "    m - (int) the number of recommendations you want for the user\n",
    "    \n",
    "    OUTPUT:\n",
    "    recs - (list) a list of recommendations for the user by article id\n",
    "    rec_names - (list) a list of recommendations for the user by article title\n",
    "    \n",
    "    Description:\n",
    "    Loops through the users based on closeness to the input user_id\n",
    "    For each user - finds articles the user hasn't seen before and provides them as recs\n",
    "    Does this until m recommendations are found\n",
    "    \n",
    "    Notes:\n",
    "    * Choose the users that have the most total article interactions \n",
    "    before choosing those with fewer article interactions.\n",
    "\n",
    "    * Choose articles with the articles with the most total interactions \n",
    "    before choosing those with fewer total interactions. \n",
    "   \n",
    "    '''\n",
    "    # Your code here\n",
    "    # Get similar users\n",
    "    neighbours = get_top_sorted_users(user_id)\n",
    "    top_similar_users = neighbours['neighbor_id'].values.tolist()\n",
    "    recs = []  # recommended article Id's\n",
    "    \n",
    "    # Get articles read by the user\n",
    "    user_article_ids = list(set(get_user_articles(user_id)[0]))\n",
    "    \n",
    "    for neighbour_id in top_similar_users:\n",
    "        recs += df[df['user_id'] == neighbour_id]['article_id'].values.tolist()\n",
    "    \n",
    "    recs = list(set(recs))\n",
    "    \n",
    "    # Selecting articles not seen by User_id\n",
    "    recs = [ x for x in recs if x not in user_article_ids]\n",
    "    \n",
    "    recs_df = df[df.article_id.isin(recs)][['article_id', 'title']].drop_duplicates().head(m) \n",
    "    recs = recs_df['article_id'].values.tolist() # Get ids\n",
    "    rec_names = recs_df['title'].values.tolist() # Get title\n",
    "    \n",
    "    return recs, rec_names"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 103,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The top 10 recommendations for user 20 are the following article ids:\n",
      "[1430.0, 1314.0, 1429.0, 1338.0, 1276.0, 1432.0, 593.0, 1185.0, 993.0, 14.0]\n",
      "\n",
      "The top 10 recommendations for user 20 are the following article names:\n",
      "['using pixiedust for fast, flexible, and easier data analysis and experimentation', 'healthcare python streaming application demo', 'use deep learning for image classification', 'ml optimization using cognitive assistant', 'deploy your python model as a restful api', 'visualize data with the matplotlib library', 'upload files to ibm data science experience using the command line', 'classify tumors with machine learning', 'configuring the apache spark sql context', 'got zip code data? prep it for analytics. – ibm watson data lab – medium']\n"
     ]
    }
   ],
   "source": [
    "# Quick spot check - don't change this code - just use it to test your functions\n",
    "rec_ids, rec_names = user_user_recs_part2(20, 10)\n",
    "print(\"The top 10 recommendations for user 20 are the following article ids:\")\n",
    "print(rec_ids)\n",
    "print()\n",
    "print(\"The top 10 recommendations for user 20 are the following article names:\")\n",
    "print(rec_names)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`5.` Use your functions from above to correctly fill in the solutions to the dictionary below.  Then test your dictionary against the solution.  Provide the code you need to answer each following the comments below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 104,
   "metadata": {},
   "outputs": [],
   "source": [
    "### Tests with a dictionary of results\n",
    "\n",
    "user1_most_sim = get_top_sorted_users(1).iloc[0].neighbor_id # Find the user that is most similar to user 1 \n",
    "user131_10th_sim = get_top_sorted_users(131).iloc[9].neighbor_id # Find the 10th most similar user to user 131"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 105,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This all looks good!  Nice job!\n"
     ]
    }
   ],
   "source": [
    "## Dictionary Test Here\n",
    "sol_5_dict = {\n",
    "    'The user that is most similar to user 1.': user1_most_sim, \n",
    "    'The user that is the 10th most similar to user 131': user131_10th_sim,\n",
    "}\n",
    "\n",
    "t.sol_5_test(sol_5_dict)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`6.` If we were given a new user, which of the above functions would you be able to use to make recommendations?  Explain.  Can you think of a better way we might make recommendations?  Use the cell below to explain a better method for new users."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "I will go for  user_user_recs_part2. It's a good start recommending articles from the most active users and make sure that these articles are the most interacted articles as well. For new users, we can ask them about their preferences, then recommend top articles that are matching their preferences. Once we have more data on them, we can move to matrix factorization.\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`7.` Using your existing functions, provide the top 10 recommended articles you would provide for the a new user below.  You can test your function against our thoughts to make sure we are all on the same page with how we might make a recommendation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 106,
   "metadata": {},
   "outputs": [],
   "source": [
    "new_user = '0.0'\n",
    "\n",
    "# What would your recommendations be for this new user '0.0'?  As a new user, they have no observed articles.\n",
    "# Provide a list of the top 10 article ids you would give to \n",
    "new_user_recs = new_user_recs = get_top_article_ids(10) # Your recommendations here\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 107,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "That's right!  Nice job!\n"
     ]
    }
   ],
   "source": [
    "assert set(new_user_recs) == set(['1314.0','1429.0','1293.0','1427.0','1162.0','1364.0','1304.0','1170.0','1431.0','1330.0']), \"Oops!  It makes sense that in this case we would want to recommend the most popular articles, because we don't know anything about these users.\"\n",
    "\n",
    "print(\"That's right!  Nice job!\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### <a class=\"anchor\" id=\"Content-Recs\">Part IV: Content Based Recommendations (EXTRA - NOT REQUIRED)</a>\n",
    "\n",
    "Another method we might use to make recommendations is to perform a ranking of the highest ranked articles associated with some term.  You might consider content to be the **doc_body**, **doc_description**, or **doc_full_name**.  There isn't one way to create a content based recommendation, especially considering that each of these columns hold content related information.  \n",
    "\n",
    "`1.` Use the function body below to create a content based recommender.  Since there isn't one right answer for this recommendation tactic, no test functions are provided.  Feel free to change the function inputs if you decide you want to try a method that requires more input values.  The input values are currently set with one idea in mind that you may use to make content based recommendations.  One additional idea is that you might want to choose the most popular recommendations that meet your 'content criteria', but again, there is a lot of flexibility in how you might make these recommendations.\n",
    "\n",
    "### This part is NOT REQUIRED to pass this project.  However, you may choose to take this on as an extra way to show off your skills."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 108,
   "metadata": {},
   "outputs": [],
   "source": [
    "def tokenize(x):\n",
    "    '''\n",
    "    Tokenize a string into words. \n",
    "    \n",
    "    Args:\n",
    "        x(string): string to tokezine.\n",
    "            \n",
    "    Returns:\n",
    "        (list): list of lemmatized words\n",
    "    '''\n",
    "    # get stop words\n",
    "    stop_words = (set(stopwords.words('english')) | set(nltk.corpus.abc.words()))\n",
    "    # split each article title into individual words\n",
    "    tokens = word_tokenize(x) \n",
    "    lemmatizer = WordNetLemmatizer()\n",
    "    clean_tokens=[]\n",
    "    for token in tokens:\n",
    "        #clean each token from whitespace and punctuation, and conver to root word\n",
    "        clean_token = lemmatizer.lemmatize(token).lower().strip()\n",
    "        clean_tokens.append(clean_token)\n",
    "        \n",
    "    filtered = [word for word in clean_tokens if word not in stop_words and word.isalpha()]\n",
    "    return filtered\n",
    "\n",
    "def make_content_recs(data_id, user_id=True, m=10, df=df):\n",
    "    '''\n",
    "    This recommender goes through each article title and nltk library to finds the most common words\n",
    "    (related to content) throughout all the articles.\n",
    "    \n",
    "    The recommender will look at the sums of words in the title of each article \n",
    "    and based on the number of matches and popularity of an article.\n",
    "       \n",
    "    Args:\n",
    "        data_id (str) - id of either user or article\n",
    "        user_id (bool) - if true, make recs based on user\n",
    "        m (int) - number of recommendations to give based on term\n",
    "    Returns:\n",
    "        recs (list) - list of article ids that are recommended\n",
    "        rec_names (list) - list of article names that are recommended   \n",
    "    '''\n",
    "    \n",
    "    if(user_id):\n",
    "        user_id = data_id\n",
    "        try:\n",
    "            # get past articles read by the user\n",
    "            article_ids, _ = get_user_articles(user_id)\n",
    "        except KeyError: # user does not exist\n",
    "            print('User Doesn\\'t Exist, Recommending Top Articles')\n",
    "            recs = get_top_article_ids(m)\n",
    "            return recs, get_article_names(recs)\n",
    "    \n",
    "    else:\n",
    "        article_ids = data_id\n",
    "        \n",
    "    title_data = df.drop_duplicates(subset='article_id') #drop duplicates \n",
    "    titles = title_data[title_data.article_id.isin(list(map(float, article_ids)))].title # get article titles\n",
    "    \n",
    "    # Tokenize the words in each article title\n",
    "    title_words=[]\n",
    "    tokenized = tokenize(titles.str.cat(sep=' '))\n",
    "    title_words.extend(tokenized)\n",
    "    \n",
    "    #Find the highest occuring words\n",
    "    common_words = pd.value_counts(title_words).sort_values(ascending=False)[:10].index\n",
    "\n",
    "    top_matches={}\n",
    "    # measure of similarity: count number of occurences of each common word in other article titles\n",
    "    for word in common_words:\n",
    "        word_count = pd.Series(title_data.title.str.count(word).fillna(0)) #gets occurences of each word in title\n",
    "        top_matches[word] = word_count\n",
    "    \n",
    "    # Most common words\n",
    "    top_matches = pd.DataFrame(top_matches) \n",
    "    top_matches['top_matches'] = top_matches.sum(axis=1)\n",
    "    top_matches['article_id'] = title_data.article_id.astype(float)\n",
    "    \n",
    "    # Get most interacted with articles\n",
    "    article_occurences = pd.DataFrame({'occurences':df.article_id.value_counts()})\n",
    "\n",
    "    # Sort matches by most popular articles\n",
    "    top_matches = top_matches.merge(article_occurences, left_on='article_id', right_index=True)\n",
    "    top_matches.sort_values(['top_matches', 'occurences'], ascending=False, inplace=True)    \n",
    "    \n",
    "    # Drop already read articles\n",
    "    recs_df = top_matches[~top_matches.article_id.isin(list(map(float, article_ids)))]\n",
    "    \n",
    "    # Get rec id and names\n",
    "    recs = recs_df.article_id[:m].values.astype(str)\n",
    "    rec_names = get_article_names(recs)\n",
    "    \n",
    "    return recs, rec_names"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`2.` Now that you have put together your content-based recommendation system, use the cell below to write a summary explaining how your content based recommender works.  Do you see any possible improvements that could be made to your function?  Is there anything novel about your content based recommender?\n",
    "\n",
    "### This part is NOT REQUIRED to pass this project.  However, you may choose to take this on as an extra way to show off your skills."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This content based recommender scans through previously interacted articles. The nltk library finds the most common words in the title of each article. Based on these most common words, the recommender looks at the sums of words relevant words in the title of each article, and based on the number of matches in the titles as well as the general popularity of the article it gives back the best recommendations.\n",
    "\n",
    "If the user has not read any articles yet, then we can't really give any content based recommendations, and just return back some of the most popular articles.\n",
    "\n",
    "There is a lot of potential improvement and optimization for this recommender. For example one could construct a custom NLTK corpus which would filter out article words. Currently I used a combination of a couple standard NLTK corpora. Furthermore, If df_content had information for all articles we could expand this recommender to look through not only the title but also the body of the articles."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`3.` Use your content-recommendation system to make recommendations for the below scenarios based on the comments.  Again no tests are provided here, because there isn't one right answer that could be used to find these content based recommendations.\n",
    "\n",
    "### This part is NOT REQUIRED to pass this project.  However, you may choose to take this on as an extra way to show off your skills.\n",
    "\n",
    "I used the NLTK library to search for articles with similar keywords. If the user has no history yet, then no content-based recommendation is given, and we will return some of the most popular articles.\n",
    "\n",
    "We can improve this further by looking at similar keywords but semantically not just the exact keyword. Also, if the user doesn't like a specific article because it has deep learning contents, it doesn't mean that he or she will disklike every article with deep learning content. It will be interesting to augment content-based recommendation with some ML algorithm that can handle such situations."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 109,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "User Doesn't Exist, Recommending Top Articles\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(Index(['1429.0', '1330.0', '1431.0', '1427.0', '1364.0', '1314.0', '1293.0',\n",
       "        '1170.0', '1162.0', '1304.0'],\n",
       "       dtype='object', name='article_id'),\n",
       " ['healthcare python streaming application demo',\n",
       "  'use deep learning for image classification',\n",
       "  'apache spark lab, part 1: basic concepts',\n",
       "  'predicting churn with the spss random tree algorithm',\n",
       "  'analyze energy consumption in buildings',\n",
       "  'visualize car data with brunel',\n",
       "  'use xgboost, scikit-learn & ibm watson machine learning apis',\n",
       "  'gosales transactions for logistic regression model',\n",
       "  'insights from new york car accident reports',\n",
       "  'finding optimal locations of new store using decision optimization'])"
      ]
     },
     "execution_count": 109,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# make recommendations for a brand new user\n",
    "make_content_recs('0.0', user_id=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 110,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array(['1159.0', '390.0', '1332.0', '1324.0', '241.0', '1422.0', '969.0',\n",
       "        '729.0', '14.0', '254.0'], dtype='<U32'),\n",
       " ['got zip code data? prep it for analytics. – ibm watson data lab – medium',\n",
       "  'ibm watson facebook posts for 2015',\n",
       "  'intents & examples for ibm watson conversation',\n",
       "  'use r dataframes & ibm watson natural language understanding',\n",
       "  'pixiedust 1.0 is here! – ibm watson data lab',\n",
       "  'i am not a data scientist – ibm watson data lab',\n",
       "  'flightpredict ii: the sequel  – ibm watson data lab',\n",
       "  'introducing ibm watson studio ',\n",
       "  'apple, ibm add machine learning to partnership with watson-core ml coupling',\n",
       "  'analyze facebook data using ibm watson and watson studio'])"
      ]
     },
     "execution_count": 110,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# make a recommendations for a user who only has interacted with article id '1427.0'\n",
    "make_content_recs(['1427.0'], user_id=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### <a class=\"anchor\" id=\"Matrix-Fact\">Part V: Matrix Factorization</a>\n",
    "\n",
    "In this part of the notebook, you will build use matrix factorization to make article recommendations to the users on the IBM Watson Studio platform.\n",
    "\n",
    "`1.` You should have already created a **user_item** matrix above in **question 1** of **Part III** above.  This first question here will just require that you run the cells to get things set up for the rest of **Part V** of the notebook. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 111,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the matrix here\n",
    "user_item_matrix = pd.read_pickle('user_item_matrix.p')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 112,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th>article_id</th>\n",
       "      <th>0.0</th>\n",
       "      <th>100.0</th>\n",
       "      <th>1000.0</th>\n",
       "      <th>1004.0</th>\n",
       "      <th>1006.0</th>\n",
       "      <th>1008.0</th>\n",
       "      <th>101.0</th>\n",
       "      <th>1014.0</th>\n",
       "      <th>1015.0</th>\n",
       "      <th>1016.0</th>\n",
       "      <th>...</th>\n",
       "      <th>977.0</th>\n",
       "      <th>98.0</th>\n",
       "      <th>981.0</th>\n",
       "      <th>984.0</th>\n",
       "      <th>985.0</th>\n",
       "      <th>986.0</th>\n",
       "      <th>990.0</th>\n",
       "      <th>993.0</th>\n",
       "      <th>996.0</th>\n",
       "      <th>997.0</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>user_id</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 714 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "article_id  0.0  100.0  1000.0  1004.0  1006.0  1008.0  101.0  1014.0  1015.0  \\\n",
       "user_id                                                                         \n",
       "1           0.0    0.0     0.0     0.0     0.0     0.0    0.0     0.0     0.0   \n",
       "2           0.0    0.0     0.0     0.0     0.0     0.0    0.0     0.0     0.0   \n",
       "3           0.0    0.0     0.0     0.0     0.0     0.0    0.0     0.0     0.0   \n",
       "4           0.0    0.0     0.0     0.0     0.0     0.0    0.0     0.0     0.0   \n",
       "5           0.0    0.0     0.0     0.0     0.0     0.0    0.0     0.0     0.0   \n",
       "\n",
       "article_id  1016.0  ...  977.0  98.0  981.0  984.0  985.0  986.0  990.0  \\\n",
       "user_id             ...                                                   \n",
       "1              0.0  ...    0.0   0.0    1.0    0.0    0.0    0.0    0.0   \n",
       "2              0.0  ...    0.0   0.0    0.0    0.0    0.0    0.0    0.0   \n",
       "3              0.0  ...    1.0   0.0    0.0    0.0    0.0    0.0    0.0   \n",
       "4              0.0  ...    0.0   0.0    0.0    0.0    0.0    0.0    0.0   \n",
       "5              0.0  ...    0.0   0.0    0.0    0.0    0.0    0.0    0.0   \n",
       "\n",
       "article_id  993.0  996.0  997.0  \n",
       "user_id                          \n",
       "1             0.0    0.0    0.0  \n",
       "2             0.0    0.0    0.0  \n",
       "3             0.0    0.0    0.0  \n",
       "4             0.0    0.0    0.0  \n",
       "5             0.0    0.0    0.0  \n",
       "\n",
       "[5 rows x 714 columns]"
      ]
     },
     "execution_count": 112,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# quick look at the matrix\n",
    "user_item_matrix.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`2.` In this situation, you can use Singular Value Decomposition from [numpy](https://docs.scipy.org/doc/numpy-1.14.0/reference/generated/numpy.linalg.svd.html) on the user-item matrix.  Use the cell to perform SVD, and explain why this is different than in the lesson."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 113,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of Nans in the users to item interactions matrix is: 0\n",
      "Number of Nans in the users to latent features matrix is: 0\n",
      "Number of Nans in the segma matrix is: 0\n",
      "Number of Nans in the items to latent features matrix is: 0\n"
     ]
    }
   ],
   "source": [
    "# Perform SVD on the User-Item Matrix Here\n",
    "\n",
    "u, s, vt =  np.linalg.svd(user_item_matrix)# use the built in to get the three matrices\n",
    "print(\"Number of Nans in the users to item interactions matrix is: {}\".format(np.isnan(user_item_matrix).sum().sum()))\n",
    "print(\"Number of Nans in the users to latent features matrix is: {}\".format(np.isnan(u).sum().sum()))\n",
    "print(\"Number of Nans in the segma matrix is: {}\".format(np.isnan(s).sum().sum()))\n",
    "print(\"Number of Nans in the items to latent features matrix is: {}\".format(np.isnan(vt).sum().sum()))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can use Singular Value Decomposition because **there are no missing values (NANs) in our data**.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`3.` Now for the tricky part, how do we choose the number of latent features to use?  Running the below cell, you can see that as the number of latent features increases, we obtain a lower error rate on making predictions for the 1 and 0 values in the user-item matrix.  Run the cell below to get an idea of how the accuracy improves as we increase the number of latent features."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 114,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmQAAAFNCAYAAACuWnPfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3Xl8VdW5//HPkzlhCkNAhkAYVZxQEUVlsGprbaveVnulzhO1t9a2t9rWttf22v7ubW8H2169ba1VrFqHarW0tVWqgrMMyqyBAIEElASSACFAhvP8/tg7cAghnEBOdobv+/Xar7P32vvs/Zx1TjgPa62ztrk7IiIiIhKdlKgDEBEREenulJCJiIiIREwJmYiIiEjElJCJiIiIREwJmYiIiEjElJCJiIiIREwJmYh0SWY2y8x+ENG1zcweNLNKM5sfRQwi0rkoIRNpJTObG37RZkYdS2diZsVmttnMesSV3WhmcyMMK1nOBs4Hhrn7pKY7zexaM3uttSc1s+lmVtoWAYbn+56ZPXKIY4rNbJeZVcctQ47wum36OkS6AiVkIq1gZgXAFMCBi9r52mnteb0kSQO+HHUQrWVmqa18ygig2N13JiOeCHzK3XvGLZuiDKaL/C2I7EcJmUjrXA28BcwCronfYWbZZvZTM1tvZtvM7DUzyw73nW1mb5hZlZmVmNm1YflcM7sx7hz7tZyYmZvZF81sNbA6LPtFeI7tZrbIzKbEHZ9qZt8yszVmtiPcn29m95rZT5vE+xcz+0rTF2hmvzaznzQp+7OZ/Xu4/g0z2xiev9DMzm1F/f0YuM3Mcpu5bkH4etPiyvbWT1g3r5vZ3WE9rjWzM8PyEjMrM7Nrmpx2gJnNCWOdZ2Yj4s59TLivInwdn43bN8vMfmVmz5nZTuCcZuIdYmazw+cXmdlNYfkNwP3A5LA16T9bUT+Y2XVm9l4Y81oz+3xY3gP4OzAkvqXKzFLM7Jvhe77VzJ40s35N6vQaM9tgZlvM7NvhvguAbwH/Gp5rSWviDM9xRtzneomZTT+C17FfF7M1aUWzoKXuG2a2FNhpZmnh8542s3IzW2dmt8YdP8nMFoZ/J5vN7GetfX0i7crdtWjRkuACFAH/BpwK1AGD4vbdC8wFhgKpwJlAJjAc2AHMANKB/sCE8DlzgRvjznEt8FrctgNzgH5Adlh2ZXiONOBrwIdAVrjvdmAZcDRgwEnhsZOATUBKeNwAoCY+/rhrTgVKAAu3+wK7gCHheUuAIeG+AmB0gnVXDJwH/An4QVh2IzA37lwOpMU9Z2/9hHVTD1wX1u8PgA1hvWcCHw3ruWd4/Kxwe2q4/xeNdQv0CF/HdWE9ngJsAY6Le+424CyC/7hmNfN65gH/B2QBE4By4Nzm3sdmnnvQ/cAngNHh+zctfJ9OCfdNB0qbHP8Vgv8kDAtf52+Ax5rU6W+B7PDzsAc4Ntz/PeCRRN63ZsqHAluBC8M6Oj/czjvM1zGr8XPR3DFhHIuB/PC1pACLgDuBDGAUsBb4WHj8m8BV4XpP4Iyo//3QoqWlRS1kIgkys7MJuqKedPdFwBrgc+G+FOB64MvuvtHdG9z9DXffA1wB/NPdH3P3Onff6u6LW3Hp/3b3CnffBeDuj4TnqHf3nxJ8CR8dHnsj8B13L/TAkvDY+QQJRmNr1uUEidDmZq73KsGXeGPL26XAmx50UzWE1xtvZunuXuzua1rxWiD4Av2SmeW18nkA69z9QXdvAJ4g+HK+y933uPsLQC0wJu74v7n7K+H78G2CVqt84JMEXYoPhvX4DvB0+Fob/dndX3f3mLvvjg8iPMfZwDfcfXf4ft4PXHUYr2k/7v43d18Tvn/zgBfY91405/PAt929NHyd3wMutf279f7T3Xe5+xJgCUFi1hrPhq1gVWb2bFh2JfCcuz8X1tEcYCFBgnY4ryMRv3T3kvBv4TSC5O8ud69197UEiefl4bF1wBgzG+Du1e7+1hFeWySplJCJJO4a4AV33xJu/4F93ZYDCFpKmktO8g9SnqiS+A0z+1rYFbTNzKqAPuH1D3Wthwi+RAkfH27uIHd34HGCFj0Iks5Hw31FBC0y3wPKzOxxa+UAb3dfDvwV+GZrnheKTyAbE9SmZT3jtvfWnbtXAxUELX0jgNPjkowqgsT5qOae24whQIW774grW0/QanREzOzjZvZW2BVaRZDgDGjhKSOAZ+Jex3sEifOguGM+jFuvYf86SsQl7p4bLpfEXfeyJnV4NjD4MF9HIuLfkxEE3Z7x1/8W+173DcA44H0zW2BmnzzCa4sklQZGiiTAgrFgnwVSzazxyy0TyDWzkwi6CXcTdNE0HYtTQtBl2JydQE7c9lHNHONxcUwBvkHQ0rXC3WNmVknQLdR4rdHA8mbO8wiwPIz3WODZZo5p9Bjwgpn9EDgd+Je9wbj/AfiDmfUm6B77Ea1vGfou8A4QP66tcQB8DrA9XG+uPlojv3HFzHoSdP1uIqinee5+fgvP9Rb2bQL6mVmvuKRsOLDxSIK14Je7TxOMVfyzu9eFLVKN729zMZUA17v7682cr+AQl2zpNR5KCfCwu9/UzHUP53W06m8hvP46dx/bXHDuvhqYEbZefxp4ysz6e9f5oYV0MWohE0nMJQStDuMJxgtNIEhqXgWudvcY8ADws3CgcaqZTQ6/mB4FzjOzz4YDkfub2YTwvIuBT5tZjpmNIfhffUt6EYyjKgfSzOxOoHfc/vuB75vZWAucaGb9Ady9FFhA0DL2dGMXaHPc/d3wGvcDz7t7FYCZHW1mHwlf126CFqmGQ1ffAecvIuhyvDWurJwgobkyrL/rCZLLI3GhBT+oyAC+D7zt7iUELXTjzOwqM0sPl9PM7NgE4y8B3gD+28yyzOxEgvfu0VbEZuFz9y4EY6EyCeq+3sw+TjA2rtFmoL+Z9Ykr+zXw/yz8wYKZ5ZnZxQnGsBkoCJOW1noE+JSZfSx8v7LCgfjDDvN1LCZ4v/qZ2VEELbEtmQ9sDwf6Z4cxHG9mpwGY2ZVmlhf+bVaFz2n1Z1WkvSghE0nMNcCD7r7B3T9sXIB7gCvC8Tq3EbSULSDoGvsRwSD6DQTdNV8LyxezbwzP3QTjnjYTdCke6gv9eYJfqK0i6CLbzf7dOD8DniQYr7Md+B3BAOhGDwEncJDuyiYeIxiE/4e4skzghwQD4D8EBhJ0E2FmV5jZigTO2+gugsH18W4i+GHCVuA4gqTnSPyBoDWuguCHGFcAhK1aHyUYb7SJ4LX8iOD1JWoGwaD5TcAzwHfDcVSJOpMgoW263ErwHlYSdBfPbnyCu79P8L6sDbvphhD8WGE2QYvmDoIB/qcnGMMfw8etZvZOK2JvTEovJnj/ywk+h7cTfOZ3HMbreJigdbmY4PP7xCGu3wB8iuA/R+sIPpP3E3ThA1wArDCzaoI6urzpWECRjqTxV1Qi0g2Y2VSClo2CsOVAREQ6ALWQiXQTZpZOMCnr/UrGREQ6FiVkIt1AODaqiuAXcD+POBwREWlCXZYiIiIiEVMLmYiIiEjElJCJiIiIRKzTTQw7YMAALygoiDoMERERkUNatGjRFnc/5K3iOl1CVlBQwMKFC6MOQ0REROSQzGx9Isepy1JEREQkYkrIRERERCKmhExEREQkYkrIRERERCKmhExEREQkYkrIRERERCKmhExEREQkYklLyMzsATMrM7PlB9lvZvZLMysys6VmdkqyYhERERHpyJLZQjYLuKCF/R8HxobLTOBXSYxFREREpMNK2kz97v6KmRW0cMjFwO/d3YG3zCzXzAa7+wfJiklERKQricWc+pgTc6ch5jS4E4vFr0N9LEYsBg3hMXuPjVsPHiHmjntwbmfvSvzDAfvdg33uvvcxFp4r5sGRMd+37e57r9NYvvfk8ddtcr19+9vOScNyGT+kdxue8fBFeeukoUBJ3HZpWHZAQmZmMwla0Rg+fHi7BCciIl1XQ8ypa4hRH3PqG2LUNTj1sRj1DfvK6xqC7dqGGLX1wbKnvoE94Xp8efz2nnCpa4jtu05DkDjVx/aVBY++Xyzx++ITpvpYmGiFyU5jciVH5psfP0YJGWDNlDX76XL3+4D7ACZOnKhPoIhIF+QeJD81exqoqWugZk89O2v3Pe6qa2BXbT27aoP9u2sbqNlbHj7WBWW7w8dd4fq+5Mupi8UOaHVpC2aQkZpCRloKmWkpZKSmkJaaQlqKkZZqpKakkJ5qwXZKCumpKWSlG+mpKaSmGOmNx6QYqeGSkmKk2r7t1BQjxYzUFMLyFFJT2O+4FNv/uWmN6yns3Zdq+5+7cT3FAAMLv6It/Ka2va/xYOXBVooF10gxwywob9xOseD5+8qCRwji368um6nb/fc3l0K0Xs+sjnNL7ygjKQXy47aHAZsiikVERA5TfUOMHbvr2bG7nu276/Z73NG4vSt43FlbT01tAzv31LOrLnhs3K6pbaC+Fa0+ZpCdnkpORipZ6an7rffvkcGwvsF6VnpqkBylGGmpjUlRCmmptnc9PdX2Jk/pqcG+xvKMMLnKSNuXbGWmpR5QnpZiexMWkdaKMiGbDdxiZo8DpwPbNH5MRKR9NcSc6j37EqbGBGrHnrp9SdauOrbH72uSdNXUNhzyOtnpqfTKSqNnZho5mankZKTRv0cG+f1yyElPpUdmGjkZcY8ZwXE9MoLtnIw0sjOChCs7PZXsjFQy01KUAEmXkbSEzMweA6YDA8ysFPgukA7g7r8GngMuBIqAGuC6ZMUiItIduDvbdtWxdWctFTtr2Vq9h607a9laHWxvqd5DRbhvW5iAVe+pP+R5M9JS6J2VRu+sdHplpdErK51BvbPotbessTyN3tnpceVpe/elp2raS5GWJPNXljMOsd+BLybr+iIind3uugaqauqo2FlLZU247KylYmcdlTVh0rVzD1ura9m6M9h3sC6/XllpDOiZSb+wVeqE7CaJVLjemFDF78tMS23nVy7S/XSc0WwiIt3A7roGNlbtYmPlLkord1G+Y8/eZKtiZ+1+CVhLXYG9MtPo1zMjHCuVw4T8XPr1yKB/z0z698gI1zPo3yNIwjLS1EIl0pEpIRMRaUM1tfV7k63Sql2UVtbs267cxZbqPQc8p1dWGv16ZJCbk8GAnhmMHdSTvjlBUtU3J4O+Oen0DZOs3Jx0crOVYIl0NUrIRERaYXddA6WVNZRU7KKksoaSihpKK3exsSpIuCp21u53fEZqCkNysxjWN4dzjxnIsL7ZDOuXzdDcHIb2zWZgr0yNrxIRJWQiIvHqG2J8sG03JRU1YcK1L/EqCbsY42WmpTC0bzbD+uZw/NA+DM3NDpKusCyvZ+YBcyyJiDSlhExEuh13Z/P2PRSVVbO6bAery6pZV76TksoaPti2e78Z0FMMhuRmk983h3OOziO/bw75/XLI7xeUDVDCJSJtQAmZiHRZsZizadsuVpdVU7R5X/JVtLmaHXHTPfTJTmdUXg9OGd53b6KV3y+H/L45DM7NUpeiiCSdEjIR6fQaYk5JRQ2rwxavIPmqpqisml11+36pOKBnJmMH9uSSk4cydlBPxgzsydiBvRjQM0MTjIpIpJSQiUinUdcQY/3WGorKdrA6TLpWl1Wzprya2vrY3uMG98lizMCeXD4pn7EDewXJV15P+vbIiDB6EZGDU0ImIh3OnvoGirfUBF2Mm6v3jvVat2UndQ37xncN65vN2IE9mTJ2AGMG9ty79M5KjzB6EZHWU0ImIpHaXdfAik3bWVxSxeKSKlZs2sb6rTV7B9abwYh+OYwZ2Itzjx3E2LCbcfTAHuRk6J8wEeka9K+ZiLSbWMxZt3UnizdU7U3A3vtg+97b/Qzuk8UJQ/vwiRMG7x3fNSqvB1npunWPiHRtSshEJGm2Vu/Zm3gtLqliSUkV23cHv27skZHKicNyuWnqKCbk5zIhP5dBvbMijlhEJBpKyESkTTTEnPc+2M7C4goWbahicUklJRW7gGAur3GDevGJEweHyVdfxgzsSarm7xIRAZSQichhqqmtZ/GGKhYUV7JwfQXvbqiiOpzb66jeWZw8PJcrTx/BhPxcThjWR+O9RERaoH8hRSQhZTt2s6i4koXrK1lYXMHyTdtpiDlmcPSgXlxy8hBOK+jHxIJ+DM3NjjpcEZFORQmZiBzA3Vm7ZScLiyuCFrDiCoq31gDBvRtPys/l5mmjmFjQj1OG96VPtqaZEBE5EkrIRAQIbqo9v7iCOSs3M2flZkorg/FffXPSmVjQj8+dPpyJBf04fkgfMtJ0KyERkbakhEykG6veU88rq8qZs3IzL71fxrZddWSkpXD2mAF8YfpozhjVn1EDeui2QiIiSaaETKSbKdu+mznvBa1gbxRtpbYhRm5OOuceO5CPjh/ElLF59MjUPw0iIu1J/+qKdHHuzuqyauas3MwLKzezpKQKgOH9crhq8gjOHz+IiSP6kpaqbkgRkagoIRPpgmIx592SKv6+7APmvLeZ9eGA/JOG9eG2j47j/PFHMW5QT3VFioh0EErIRLqIxiTsuWUf8NyyD/hg224yUlOYPLo/N00ZxXnHDuKoPpoJX0SkI1JCJtKJHSwJmzpuALd/7GjOGz+I3lmakkJEpKNTQibSySgJExHpepSQiXQC8UnY35d9wCYlYSIiXYoSMpEObG15NX94ewPPNUnCblMSJiLSpSghE+lg3J1XV2/hwdfX8XJhOempxrRxeUrCRES6MCVkIh3ErtoG/vRuKbNeL2Z1WTUDembylfPGcsXpI8jrlRl1eCIikkRJTcjM7ALgF0AqcL+7/7DJ/hHAA0AeUAFc6e6lyYxJpKPZVLWL37+5nsfmb2DbrjqOH9qbn152Ep88aTCZaalRhyciIu0gaQmZmaUC9wLnA6XAAjOb7e4r4w77CfB7d3/IzD4C/DdwVbJiEuko3J1F6yt58PVi/rHiQ9ydjx13FNefPZKJI/pqwlYRkW4mmS1kk4Aid18LYGaPAxcD8QnZeOCr4frLwLNJjEckcrX1Mf62bBMPvl7M0tJt9M5K44azR3L15BEM65sTdXgiIhKRZCZkQ4GSuO1S4PQmxywBPkPQrfkvQC8z6+/uW5MYl0i721K9hz+8vYGH31pP+Y49jMrrwfcvOZ7PnDKUnAwN5RQR6e6S+U3QXJ+LN9m+DbjHzK4FXgE2AvUHnMhsJjATYPjw4W0bpUgSfbBtF7+eu4bHFpRQWx9j2rg8rru0gKlj80hJUbekiIgEkpmQlQL5cdvDgE3xB7j7JuDTAGbWE/iMu29reiJ3vw+4D2DixIlNkzqRDqekooZfzVvDUwtLibnz6VOGMnPqaMYM7Bl1aCIi0gElMyFbAIw1s5EELV+XA5+LP8DMBgAV7h4D7iD4xaVIp7V+607ufbmIP72zETP47MR8bp42mvx+Gh8mIiIHl7SEzN3rzewW4HmCaS8ecPcVZnYXsNDdZwPTgf82MyfosvxisuIRSaY15dXc+1IRf16yidQU48ozRvD5aaMY3Cc76tBERKQTMPfO1QM4ceJEX7hwYdRhiACwavMO7nmpiL8s3URmWgpXnj6CmVNHMbB3VtShiYhIB2Bmi9x94qGO08+7RA7Dyk3buefl1Ty37ENyMlKZOXUUN00ZxYCemlFfRERaTwmZSCssK93GL19azZyVm+mVmcYt54zh+rNH0q9HRtShiYhIJ6aETCQByzdu46cvFPJyYTm9s9L4ynljue7MkfTJ0Y2+RUTkyCkhE2nB6s07+NmcVfx9+Yf0yU7n9o8dzdWTR9ArS4mYiIi0HSVkIs0oqajh7n+u4tl3N5Kdnsqt547lxikj6a1ETEREkkAJmUiczdt3878vreaJBSWkmHHjlFHcPG20xoiJiEhSKSETASp21vLreWt46I1iGmLO5ZPyueWcsRzVR9NXiIhI8ikhk25t++467n91HQ+8to6a2nouOXkoXzl3HMP7a2Z9ERFpP0rIpFvaVdvAQ28W8+t5a6iqqePjxx/Fv58/jrGDekUdmoiIdENKyKRbqa2P8fiCDfzvS0WU79jDtHF53PbRozlhWJ+oQxMRkW5MCZl0G2+u2co3nl7KhooaJhX0497PncKkkf2iDktEREQJmXR9O/fU86N/vM/v31zPiP45zLruNKaNy8PMog5NREQEUEImXdxba7dy+1NLKKnYxbVnFvD1C44mJ0MfexER6Vj0zSRdUk1tPT/6+/s89OZ6hvfL4YmZZ3D6qP5RhyUiItIsJWTS5by9diu3PxWMFVOrmIiIdAb6lpIuo6a2nv/5RyGz3ihmeL8cHp95BmeoVUxERDoBJWTSJcxfV8HtTy1h/Va1iomISOejbyzp1HbVNvA/z7/PrDeKGdY3m8duOoPJo9UqJiIinYsSMum0FhRXcPsfl1C8tYZrJo/g6xccQ49MfaRFRKTz0beXdDq7ahv48fOFPPjGOrWKiYhIl6CETDqVTVW7uOaB+awuq+bqySP4hlrFRESkC9A3mXQaRWU7uOp386neXc/DN0xiyti8qEMSERFpE0rIpFNYXFLFdQ/OJzUlhSc+P5nxQ3pHHZKIiEibUUImHd4rq8q5+ZFFDOiZycM3TGJE/x5RhyQiItKmlJBJhzZ7ySa+9uRixgzsxUPXn8bAXllRhyQiItLmlJBJh/X7N4v57uwVnFbQj99ePZE+2elRhyQiIpIUSsikw3F37v7nan754mrOHz+I/51xMlnpqVGHJSIikjRKyKRDaYg53529nEfe2sBnJw7jv/7lBNJSU6IOS0REJKmS+k1nZheYWaGZFZnZN5vZP9zMXjazd81sqZldmMx4pGPbU9/ArY+/yyNvbeDmaaP50WdOVDImIiLdQtJayMwsFbgXOB8oBRaY2Wx3Xxl32HeAJ939V2Y2HngOKEhWTNJxVe+p5+aHF/Fa0Ra+feGx3DR1VNQhiYiItJtkdllOAorcfS2AmT0OXAzEJ2QONE4o1QfYlMR4pIPaWr2H62YtYMWm7fzkspO49NRhUYckIiLSrpKZkA0FSuK2S4HTmxzzPeAFM/sS0AM4L4nxSAdUWlnD1Q/MZ2PlLu676lTOPXZQ1CGJiIi0u2QO0LFmyrzJ9gxglrsPAy4EHjazA2Iys5lmttDMFpaXlychVInCqs07uPRXb7Jlxx4eufF0JWMiItJtJTMhKwXy47aHcWCX5A3AkwDu/iaQBQxoeiJ3v8/dJ7r7xLw83b+wK3hnQyWX/fpNYu488fnJnFbQL+qQREREIpPMhGwBMNbMRppZBnA5MLvJMRuAcwHM7FiChExNYF3cuxsquer+t+mbk87TXziTYwfrvpQiItK9JW0MmbvXm9ktwPNAKvCAu68ws7uAhe4+G/ga8Fsz+ypBd+a17t60W1O6kOUbt3HNA/MZ0CuTJz4/mUG9dSskERGRpE4M6+7PEUxlEV92Z9z6SuCsZMYgHUfhhzu46ndv0ysrnUdvPF3JmIiISEizbkq7WFtezRX3v01GWgp/uOl0hvXNiTokERGRDkMJmSRdSUUNn/vt27g7j954BiP694g6JBERkQ5F97KUpNpUtYsZv32L3fUNPHbTGYwZ2DPqkERERDoctZBJ0pRt380V97/Ntpo6Hr7+dP2aUkRE5CDUQiZJsbV6D1fc/zabt+/m4RsmccKwPlGHJCIi0mGphUza3LaaOq763Xw2VNTwu2tO49QRmvRVRESkJUrIpE3t2F3H1Q/Op6ismvuunsjk0f2jDklERKTDU0Imbaamtp7rZy1gxcZt3HvFKUwbp9tciYiIJEIJmbSJ3XUN3PT7hSxaX8kvLj+Z88frRuEiIiKJ0qB+OWK19TG+8Mgi3lizlZ9edhKfOHFw1CGJiIh0KmohkyNS3xDjS4+9w8uF5fzXv5zAp08ZFnVIIiIinY4SMjlsDTHn359cwvMrNvPdT41nxqThUYckIiLSKSkhk8Pi7nzn2WXMXrKJb1xwDNedNTLqkERERDotJWRyWP5v7hoem1/CF88ZzRemj446HBERkU5NCZm02uwlm/jx84VcMmEIt3306KjDERER6fSUkEmrLCyu4LY/LmHSyH786NITMbOoQxIREen0lJBJwoq37OSm3y9kWG429111KplpqVGHJCIi0iUoIZOEVO6s5bpZCzAzHrzuNHJzMqIOSUREpMvQxLBySLvrGpj58EI2Vu3isZtOZ0T/HlGHJCIi0qWohUxaFIs5X39qKQuKK7n7sxM4dUS/qEMSERHpcpSQSYt+NmfV3rnGdEskERGR5FBCJgf15IIS7nm5iBmT8rl52qiowxEREemylJBJs15bvYVvPbOMKWMHcNfFx2t6CxERkSQ6ZEJmZreYWd/2CEY6hsIPd/CFRxYxZmBP/u+KU0hPVd4uIiKSTIl80x4FLDCzJ83sAlNTSZdWtn03189aQHZGKg9cexq9stKjDklERKTLO2RC5u7fAcYCvwOuBVab2X+ZmW5g2MXU1NZzw0MLqayp5YFrT2NIbnbUIYmIiHQLCfVFubsDH4ZLPdAXeMrM/ieJsUk7aog5tz62mBWbtvG/M07m+KF9og5JRESk2zjkxLBmditwDbAFuB+43d3rzCwFWA18PbkhSnv4wd9W8s/3NvOfFx3HuccOijocERGRbiWRmfoHAJ929/Xxhe4eM7NPtvREM7sA+AWQCtzv7j9ssv9u4JxwMwcY6O65iQYvbWPW6+t48PVirj9rJNecWRB1OCIiIt1OIgnZc0BF44aZ9QLGu/vb7v7ewZ5kZqnAvcD5QCnBDwNmu/vKxmPc/atxx38JOLn1L0GOxEvvb+auv67k/PGD+PYnjo06HBERkW4pkTFkvwKq47Z3hmWHMgkocve17l4LPA5c3MLxM4DHEjivtJGqmlq+/tRSjjmqN7+4fAKpKfoBrYiISBQSScgsHNQPBF2VJNayNhQoidsuDcsOvIDZCGAk8FIC55U28v/+9h6VNXX85LKTyMnQfeZFRESikkhCttbMbjWz9HD5MrA2gec119zizZQBXA485e4NzZ7IbKaZLTSzheXl5QlcWg7ljaIt/HFRKTdNGcX4Ib2jDkdERKRbSyQhuxk4E9hI0Mp1OjAzgeeVAvlx28OATQc59nJa6K509/u74WWFAAAa0klEQVTcfaK7T8zLy0vg0tKS3XUN3PHMMkb0z+Er542NOhwREZFu75D9VO5eRpAwtdYCYKyZjSRI5i4HPtf0IDM7mmBeszcP4xpyGH7x4mrWb63h0RtPJys9NepwREREur1E5iHLAm4AjgOyGsvd/fqWnufu9WZ2C/A8wbQXD7j7CjO7C1jo7rPDQ2cAj8ePU5PkWblpO/e9spbLTh3GWWMGRB2OiIiIkNjg/IeB94GPAXcBVwAHne4inrs/RzBtRnzZnU22v5fIueTINcScO/60lL456ZriQkREpANJZAzZGHf/D2Cnuz8EfAI4IblhSTLMeqOYJaXbuPNTx5GbkxF1OCIiIhJKJCGrCx+rzOx4oA9QkLSIJClKKmr46QuFnHN0Hp86cXDU4YiIiEicRLos7zOzvsB3gNlAT+A/khqVtCl35zvPLgfgB/9yAmaaAFZERKQjaTEhC28gvt3dK4FXgFHtEpW0qdlLNjFvVTl3fnI8Q3Ozow5HREREmmixyzKclf+WdopFkqByZy13/WUlJ+Xn6sbhIiIiHVQiXZZzzOw24AmC+1gC4O4VB3+KdBQ/+Nt7bNtVxyOfPkH3qhQREemgEknIGucb+2JcmaPuyw7vtdVbePqdUv5t+miOHazbI4mIiHRUiczUP7I9ApG2tau2gW89s4yRA3pw67m6PZKIiEhHlshM/Vc3V+7uv2/7cKSt/PzFVWyoqOGxm87Q7ZFEREQ6uES6LE+LW88CzgXeAZSQdVDLN27j/lfX8a8T85k8un/U4YiIiMghJNJl+aX4bTPrQ3A7JemA6hti3PGnZfTNyeBbF+r2SCIiIp1BIi1kTdUAGpTUQc16o5hlG7dxz+dOpk9OetThiIiISAISGUP2F4JfVUIwb9l44MlkBiWHJ7g90irOPWYgnzhBt0cSERHpLBJpIftJ3Ho9sN7dS5MUjxwmd+fbzy4nxeD7lxyv2yOJiIh0IokkZBuAD9x9N4CZZZtZgbsXJzUyaZU/L97EK6vK+d6nxjNEt0cSERHpVFq8dVLoj0AsbrshLJMOomJnLXf9dSUT8nO5anJB1OGIiIhIKyWSkKW5e23jRriekbyQpLV+8kIh23fV8cPP6PZIIiIinVEiCVm5mV3UuGFmFwNbkheStEZJRQ1PLihhxqThHHOUbo8kIiLSGSUyhuxm4FEzuyfcLgWanb1f2t89LxWRkmL82zmjow5FREREDlMiE8OuAc4ws56AufuO5IcliVi/dSdPvVPKVWeMYHAfDeQXERHprA7ZZWlm/2Vmue5e7e47zKyvmf2gPYKTlv3yxSLSUowvTFfrmIiISGeWyBiyj7t7VeOGu1cCFyYvJEnEui07eebdUq48YwSDemdFHY6IiIgcgUQSslQzy2zcMLNsILOF46Ud/PLF1WSkpXDzNLWOiYiIdHaJDOp/BHjRzB4Mt68DHkpeSHIoRWXV/HnxRm6cMoq8XsqNRUREOrtEBvX/j5ktBc4DDPgHMCLZgcnB/fLF1WSlp/L5qaOiDkVERETaQCJdlgAfEszW/xngXOC9pEUkLVq1eQd/WbqJqycX0L+nWsdERES6goO2kJnZOOByYAawFXiCYNqLc9opNmnGL15cTU56KjPVOiYiItJltNRl+T7wKvApdy8CMLOvtktU0qz3P9zO35Z+wBfPGU2/Hrp7lYiISFfRUpflZwi6Kl82s9+a2bkEY8gSZmYXmFmhmRWZ2TcPcsxnzWylma0wsz+05vzdzS/+uZqemWncNEWtYyIiIl3JQVvI3P0Z4Bkz6wFcAnwVGGRmvwKecfcXWjqxmaUC9wLnE9xuaYGZzXb3lXHHjAXuAM5y90ozG3jEr6iLWrFpG39f/iG3njuW3By1jomIiHQlhxzU7+473f1Rd/8kMAxYDDTb2tXEJKDI3de6ey3wOHBxk2NuAu4NJ5vF3ctaFX038vN/rqZXVho3nD0y6lBERESkjSX6K0sA3L3C3X/j7h9J4PChQEncdmlYFm8cMM7MXjezt8zsgtbE010sK93GnJWbufHsUfTJTo86HBEREWljiUwMe7iaG2/mzVx/LDCdoPXtVTM7Pv5WTQBmNhOYCTB8+PC2j7SD+/k/V9EnO53rzi6IOhQRERFJgla1kLVSKZAftz0M2NTMMX929zp3XwcUEiRo+3H3+9x9ortPzMvLS1rAHdGSkipefL+Mm6aMpHeWWsdERES6omQmZAuAsWY20swyCOY0m93kmGeBcwDMbABBF+baJMbU6dz9z1Xk5qRz7VkaOyYiItJVJS0hc/d64BbgeYKZ/Z909xVmdpeZXRQe9jyw1cxWAi8Dt7v71mTF1NksWl/J3MJyPj91ND0zk9m7LCIiIlFK6re8uz8HPNek7M64dQf+PVykiZ//cxX9emRw9WTdOlRERKQrS2aXpRyBBcUVvLp6CzdPG0UPtY6JiIh0aUrIOqi756xiQM8MrjxDrWMiIiJdnRKyDuittVt5Y81Wbp42mpwMtY6JiIh0dUrIOqC756wir1emWsdERES6CSVkHcwba7bw9roK/m36aLLSU6MOR0RERNqBErIOxN25e84qBvXOZMak7ndHAhERke5KCVkH8lrRFhYUV3LLOWPUOiYiItKNKCHrIBpbx4b0yeKzp+Uf+gkiIiLSZSgh6yBeWb2FdzZU8cWPjCEzTa1jIiIi3YkSsg7i3peLGNIni8tOVeuYiIhId6OErAN4Z0Ml89dVcMOUUWSk6S0RERHpbvTt3wH8Zt4a+mSnc7nGjomIiHRLSsgitqa8mhdWbuaqM0bonpUiIiLdlBKyiN3/6lrSU1O49qyCqEMRERGRiCghi1DZ9t08vWgjl506jAE9M6MOR0RERCKihCxCD75RTH0sxk1TRkUdioiIiERICVlEduyu45G31vPx4wdTMKBH1OGIiIhIhJSQReSx+RvYsbuemVPVOiYiItLdKSGLQG19jN+9to7Jo/pzUn5u1OGIiIhIxJSQReDZxRvZvH0PN08fHXUoIiIi0gEoIWtnsZhz3ytrOXZwb6aOHRB1OCIiItIBKCFrZy+9X0ZRWTU3TxuFmUUdjoiIiHQASsja2a/nrWFobjYXnjA46lBERESkg1BC1o4WFlewcH0lN04ZSXqqql5EREQCygra0W9eWUtuTjr/qpuIi4iISBwlZO2kqGwHc1Zu5urJBeRk6CbiIiIiso8SsnZy3ytryUpP4ZrJI6IORURERDoYJWTtYPP23Tzz7kY+OzGf/rqJuIiIiDSR1ITMzC4ws0IzKzKzbzaz/1ozKzezxeFyYzLjicoDr62jIebceLZukyQiIiIHStpgJjNLBe4FzgdKgQVmNtvdVzY59Al3vyVZcURt++46Hn17AxeeMJjh/XOiDkdEREQ6oGS2kE0Citx9rbvXAo8DFyfxeh3So29toHpPPTdP022SREREpHnJTMiGAiVx26VhWVOfMbOlZvaUmXWp+SD21DfwwOvrOHvMAI4f2ifqcERERKSDSmZC1tx9gbzJ9l+AAnc/Efgn8FCzJzKbaWYLzWxheXl5G4eZPM++u5HyHXvUOiYiIiItSmZCVgrEt3gNAzbFH+DuW919T7j5W+DU5k7k7ve5+0R3n5iXl5eUYNtaLOb85pW1HDekN2eN6R91OCIiItKBJTMhWwCMNbORZpYBXA7Mjj/AzOJv6HgR8F4S42lXc97bzNrynXx+2mjdRFxERERalLRfWbp7vZndAjwPpAIPuPsKM7sLWOjus4FbzewioB6oAK5NVjztyd359bw15PfL5sLjj4o6HBEREengknoPH3d/DniuSdmdcet3AHckM4YoLCiu5N0NVdx18XGk6SbiIiIicgjKFpLgN/PW0K9HBped2qV+NCoiIiJJooSsja3avIMX3y/jmskFZGekRh2OiIiIdAJKyNrYb+atJTs9lat1E3ERERFJkBKyNvTBtl38efFG/vW0fPr2yIg6HBEREekklJC1oT+9s5H6mHPD2SOjDkVEREQ6ESVkbWheYTnHDelNfj/dRFxEREQSp4SsjWzbVceiDZWcc/TAqEMRERGRTkYJWRt5bfUWGmLO9KM7x62dREREpONQQtZG5haW0TsrjQn5uVGHIiIiIp2MErI24O7MXVXOlHF5mplfREREWk3ZQxtYsWk75Tv2aPyYiIiIHBYlZG1g3qpyAKaN0/gxERERaT0lZG1gbmEZxw/tTV6vzKhDERERkU5ICdkR2lZTx6L1mu5CREREDp8SsiP0alE5MUfTXYiIiMhhU0J2hOYWltMnO50J+X2jDkVEREQ6KSVkRyAWc+YWljN1XB6pKRZ1OCIiItJJKSE7Ais/2M6W6j1M168rRURE5AgoITsCcwvLAJiqhExERESOgBKyIzC3sJwThvbRdBciIiJyRJSQHaaqmlre2VDJOfp1pYiIiBwhJWSH6dXVW4g5TNP8YyIiInKElJAdprmF5eTmpDMhPzfqUERERKSTU0J2GGIxZ96qMqaO1XQXIiIicuSUkB2GFZu2s6W6VrPzi4iISJtQQnYYNN2FiIiItCUlZIdh7qpyThzWhwE9Nd2FiIiIHDklZK1UVVPLuxsqma5fV4qIiEgbSWpCZmYXmFmhmRWZ2TdbOO5SM3Mzm5jMeNrCK+F0Fxo/JiIiIm0laQmZmaUC9wIfB8YDM8xsfDPH9QJuBd5OVixtaW5hGX1z0jlpmKa7EBERkbaRzBaySUCRu69191rgceDiZo77PvA/wO4kxtImYjFnXmE5U8dpugsRERFpO8lMyIYCJXHbpWHZXmZ2MpDv7n9NYhxtZvmmbWzdqekuREREpG0lMyFrrgnJ9+40SwHuBr52yBOZzTSzhWa2sLy8vA1DbJ25heWYwdSxSshERESk7SQzISsF8uO2hwGb4rZ7AccDc82sGDgDmN3cwH53v8/dJ7r7xLy86JKhuYVlnDgsl/6a7kJERETaUDITsgXAWDMbaWYZwOXA7Mad7r7N3Qe4e4G7FwBvARe5+8IkxnTYKnfW8m5JFdM1GayIiIi0saQlZO5eD9wCPA+8Bzzp7ivM7C4zuyhZ102WV1aX45ruQkRERJIgLZknd/fngOealN15kGOnJzOWIzWvsJy+OemcqOkuREREpI1ppv4ExGLOvFXlTNN0FyIiIpIESsgSsGxj43QXul2SiIiItD0lZAnYO92FBvSLiIhIEighS8DcVWWcNCyXfj0yog5FREREuiAlZIdQsbOWxSVV+nWliIiIJI0SskN4de90Fxo/JiIiIsmhhOwQ5haW069HBicO7RN1KCIiItJFKSFrQfx0Fyma7kJERESSRAlZC5Zu3EbFzlqNHxMREZGkUkLWgrmFZZjBlLFKyERERCR5lJC1YG5hORPyNd2FiIiIJJcSsoPYWr2HJaVVTB+nX1eKiIhIcikhO4hXV28Jp7tQd6WIiIgklxKyg5hbWEb/HhmcoOkuREREJMmUkDWjQdNdiIiISDtSQtaMpaVVVNbUMU3dlSIiItIOlJA1Y25hOSkGUzXdhYiIiLQDJWTNmLsqmO6ir6a7EBERkXaghKyJrdV7WFpapZuJi4iISLtRQtaEprsQERGR9pYWdQAdzSdOHMyQ3GyOH6LpLkRERKR9KCFrIj01hUkj+0UdhoiIiHQj6rIUERERiZgSMhEREZGIKSETERERiZgSMhEREZGIKSETERERiZgSMhEREZGIKSETERERiZgSMhEREZGIKSETERERiZgSMhEREZGImbtHHUOrmFk5sL4NTjUA2NIG5+kOVFeJUT0lTnWVGNVT4lRXiVNdJaat6mmEu+cd6qBOl5C1FTNb6O4To46jM1BdJUb1lDjVVWJUT4lTXSVOdZWY9q4ndVmKiIiIREwJmYiIiEjEunNCdl/UAXQiqqvEqJ4Sp7pKjOopcaqrxKmuEtOu9dRtx5CJiIiIdBTduYVMREREpEPodgmZmV1gZoVmVmRm34w6nqiZ2QNmVmZmy+PK+pnZHDNbHT72DcvNzH4Z1t1SMzslusjbl5nlm9nLZvaema0wsy+H5aqrJswsy8zmm9mSsK7+MywfaWZvh3X1hJllhOWZ4XZRuL8gyvjbm5mlmtm7ZvbXcFv11AwzKzazZWa22MwWhmX6+2uGmeWa2VNm9n74b9Zk1dWBzOzo8PPUuGw3s69EVVfdKiEzs1TgXuDjwHhghpmNjzaqyM0CLmhS9k3gRXcfC7wYbkNQb2PDZSbwq3aKsSOoB77m7scCZwBfDD87qqsD7QE+4u4nAROAC8zsDOBHwN1hXVUCN4TH3wBUuvsY4O7wuO7ky8B7cduqp4M7x90nxE1FoL+/5v0C+Ie7HwOcRPD5Ul014e6F4edpAnAqUAM8Q1R15e7dZgEmA8/Hbd8B3BF1XFEvQAGwPG67EBgcrg8GCsP13wAzmjuuuy3An4HzVVeHrKcc4B3gdIIJFtPC8r1/i8DzwORwPS08zqKOvZ3qZxjBP/gfAf4KmOrpoHVVDAxoUqa/vwPrqTewrulnQ3V1yHr7KPB6lHXVrVrIgKFASdx2aVgm+xvk7h8AhI8Dw3LVHxB2FZ0MvI3qqllhN9xioAyYA6wBqty9Pjwkvj721lW4fxvQv30jjszPga8DsXC7P6qng3HgBTNbZGYzwzL9/R1oFFAOPBh2hd9vZj1QXR3K5cBj4XokddXdEjJrpkw/M01ct68/M+sJPA18xd23t3RoM2Xdpq7cvcGDboBhwCTg2OYOCx+7ZV2Z2SeBMndfFF/czKHdup7inOXupxB0G33RzKa2cGx3rqs04BTgV+5+MrCTfV1uzenOdQVAOE7zIuCPhzq0mbI2q6vulpCVAvlx28OATRHF0pFtNrPBAOFjWVjerevPzNIJkrFH3f1PYbHqqgXuXgXMJRh3l2tmaeGu+PrYW1fh/j5ARftGGomzgIvMrBh4nKDb8ueonprl7pvCxzKCcT6T0N9fc0qBUnd/O9x+iiBBU10d3MeBd9x9c7gdSV11t4RsATA2/BVTBkET5eyIY+qIZgPXhOvXEIyXaiy/OvylyRnAtsZm3a7OzAz4HfCeu/8sbpfqqgkzyzOz3HA9GziPYFDxy8Cl4WFN66qxDi8FXvJwgEZX5u53uPswdy8g+LfoJXe/AtXTAcysh5n1alwnGO+zHP39HcDdPwRKzOzosOhcYCWqq5bMYF93JURVV1EPpItg4N6FwCqCMS3fjjqeqJfwQ/gBUEeQ/d9AMC7lRWB1+NgvPNYIfqW6BlgGTIw6/nasp7MJmqaXAovD5ULVVbN1dSLwblhXy4E7w/JRwHygiKBrIDMszwq3i8L9o6J+DRHU2XTgr6qng9bPKGBJuKxo/Ldbf38Hra8JwMLwb/BZoK/q6qB1lQNsBfrElUVSV5qpX0RERCRi3a3LUkRERKTDUUImIiIiEjElZCIiIiIRU0ImIiIiEjElZCIiIiIRU0Im0s2ZmZvZT+O2bzOz77XRuWeZ2aWHPvKIr3OZmb1nZi83KS8ws+WtOM8l4U3jDzeOAjP7XAv7dpnZ4rgloy2vISKdlxIyEdkDfNrMBkQdSDwzS23F4TcA/+bu5xzhZS8BDjshAwqAlpKlNe4+IW6pTcI1mtXK+hSRdqaETETqgfuArzbd0bSFy8yqw8fpZjbPzJ40s1Vm9kMzu8LM5pvZMjMbHXea88zs1fC4T4bPTzWzH5vZAjNbamafjzvvy2b2B4KJF5vGMyM8/3Iz+1FYdifBxL2/NrMfJ/KCzeym8NpLzOxpM8sxszMJ7mf347D1anS4/CO8ofWrZnZMXL380szeMLO1cXX0Q2BK+PwD6vMgsfQwswfCeN41s4vD8oLwmu+Ey5nNXcPMrjWze+LO91czmx6uV5vZXWb2NjDZzE4N37dFZva87bs9zK1mtjJ8Lx5PJG4RaWNRz5KrRYuWaBegGugNFBPcH/E24HvhvlnApfHHho/TgSpgMJAJbAT+M9z3ZeDncc//B8F//sYS3A0iC5gJfCc8JpNgVvGR4Xl3AiObiXMIsAHII7iB8kvAJeG+uTQzazZBa9LyZsr7x63/APjSQV7vi8DYcP10gtsVNR73x/B1jQeK4urlrwep5wJgF/vu9HBvWP5fwJXhei7BnUR6EMwgnhWWjwUWNncN4FrgnrjtvwLTw3UHPhuupwNvAHnh9r8CD4Trm9h3R4DcqD+TWrR0x6XxBrYi0o25+3Yz+z1wK0HSkIgFHt7HzczWAC+E5cuA+K7DJ909Bqw2s7XAMQT3IjwxrmWpD0HSUQvMd/d1zVzvNGCuu5eH13wUmEpwa5jWOt7MfkCQAPUEnm96gJn1BM4E/mhmjcWZcYc8G76ulWY2KMHrrnH3CU3KPkpwk/Hbwu0sYDhBknSPmU0AGoBxCV4jXgPwdLh+NHA8MCd8PakEt02D4BY7j5rZsxxefYrIEVJCJiKNfg68AzwYV1ZPOLTBgm/x+EHoe+LWY3HbMfb/t6Xp/dmc4J5wX3L3/RKhsKtt50His4OUH45ZBK1rS8zsWoJWp6ZSgKpmEqhG8a//SGIz4DPuXrhfYfDDis3ASWEsuw/y/L3vUSgrbn23uzfEXWeFu09u5hyfIEhuLwL+w8yOc/f61r4QETl8GkMmIgC4ewXwJMEA+UbFwKnh+sUE3V6tdZmZpYTjykYBhQQtUl8ws3QAMxtnZj0OcZ63gWlmNiAcoD4DmHcY8QD0Aj4Ir39FXPmOcB/uvh1YZ2aXhTGamZ10iPPufX4rPA98KUx4MbOTw/I+wAdhK9xVBC1azV2jGJgQ1nE+MOkg1ykE8sxscniddDM7zsxSgHx3fxn4OvtaDUWkHSkhE5F4PwXif235W4IkaD7BGKqDtV61pJAgcfo7cLO77wbuB1YC71gwLcVvOESLfdg9egfwMrAEeMfd/5zA9Y82s9K45TLgPwgSvDnA+3HHPg7cHg6uH02QrN1gZkuAFQRJaUuWAvXhjwUSGtQPfJ8g0V0a1sX3w/L/A64xs7cIuisb677pNV4H1hF0Ff+EoJXzAB78ovNS4Efh61lM0CWbCjxiZsuAd4G73b0qwdhFpI2Ye9PeBBERERFpT2ohExEREYmYEjIRERGRiCkhExEREYmYEjIRERGRiCkhExEREYmYEjIRERGRiCkhExEREYmYEjIRERGRiP1/NX6KKuQfljoAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 720x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "num_latent_feats = np.arange(10,700+10,20)\n",
    "sum_errs = []\n",
    "\n",
    "for k in num_latent_feats:\n",
    "    # restructure with k latent features\n",
    "    s_new, u_new, vt_new = np.diag(s[:k]), u[:, :k], vt[:k, :]\n",
    "    \n",
    "    # take dot product\n",
    "    user_item_est = np.around(np.dot(np.dot(u_new, s_new), vt_new))\n",
    "    \n",
    "    # compute error for each prediction to actual value\n",
    "    diffs = np.subtract(user_item_matrix, user_item_est)\n",
    "    \n",
    "    # total errors and keep track of them\n",
    "    err = np.sum(np.sum(np.abs(diffs)))\n",
    "    sum_errs.append(err)\n",
    "    \n",
    "plt.figure(figsize=(10,5))\n",
    "plt.plot(num_latent_feats, 1 - np.array(sum_errs)/df.shape[0]);\n",
    "plt.xlabel('Number of Latent Features');\n",
    "plt.ylabel('Accuracy');\n",
    "plt.title('Accuracy vs. Number of Latent Features');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`4.` From the above, we can't really be sure how many features to use, because simply having a better way to predict the 1's and 0's of the matrix doesn't exactly give us an indication of if we are able to make good recommendations.  Instead, we might split our dataset into a training and test set of data, as shown in the cell below.  \n",
    "\n",
    "Use the code from question 3 to understand the impact on accuracy of the training and test sets of data with different numbers of latent features. Using the split below: \n",
    "\n",
    "* How many users can we make predictions for in the test set?  \n",
    "* How many users are we not able to make predictions for because of the cold start problem?\n",
    "* How many articles can we make predictions for in the test set?  \n",
    "* How many articles are we not able to make predictions for because of the cold start problem?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 115,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\ahmed\\Anaconda3\\lib\\site-packages\\ipykernel_launcher.py:41: DeprecationWarning: \n",
      ".ix is deprecated. Please use\n",
      ".loc for label based indexing or\n",
      ".iloc for positional indexing\n",
      "\n",
      "See the documentation here:\n",
      "http://pandas.pydata.org/pandas-docs/stable/indexing.html#ix-indexer-is-deprecated\n"
     ]
    }
   ],
   "source": [
    "df_train = df.head(40000)\n",
    "df_test = df.tail(5993)\n",
    "\n",
    "def create_test_and_train_user_item(df_train, df_test):\n",
    "    '''\n",
    "    INPUT:\n",
    "    df_train - training dataframe\n",
    "    df_test - test dataframe\n",
    "    \n",
    "    OUTPUT:\n",
    "    user_item_train - a user-item matrix of the training dataframe \n",
    "                      (unique users for each row and unique articles for each column)\n",
    "    user_item_test - a user-item matrix of the testing dataframe \n",
    "                    (unique users for each row and unique articles for each column)\n",
    "    test_idx - all of the test user ids\n",
    "    test_arts - all of the test article ids\n",
    "    \n",
    "    '''\n",
    "    # Your code here\n",
    "    # create user item matrix for the train dataset\n",
    "    user_item_train = create_user_item_matrix(df_train)\n",
    "    \n",
    "    # create the test dataset\n",
    "    user_item_test = create_user_item_matrix(df_test)\n",
    "    \n",
    "    # get the ids of the train dataset and test dataset\n",
    "    train_idx = set(user_item_train.index)\n",
    "    test_idx = set(user_item_test.index)\n",
    "    \n",
    "    # get shared rows\n",
    "    shared_rows = train_idx.intersection(test_idx)\n",
    "    \n",
    "    # get columns in train and test datasets\n",
    "    train_arts = set(user_item_train.columns)\n",
    "    test_arts = set(user_item_test.columns)\n",
    "    \n",
    "    # get shared columns\n",
    "    shared_cols = train_arts.intersection(test_arts)\n",
    "    \n",
    "    # Creating new user-item matrix for tets with common values\n",
    "    user_item_test = user_item_test.ix[shared_rows, shared_cols]\n",
    "    \n",
    "    return user_item_train, user_item_test, test_idx, test_arts\n",
    "\n",
    "user_item_train, user_item_test, test_idx, test_arts = create_test_and_train_user_item(df_train, df_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 117,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Awesome job!  That's right!  All of the test movies are in the training data, but there are only 20 test users that were also in the training set.  All of the other users that are in the test set we have no data on.  Therefore, we cannot make predictions for these users using SVD.\n"
     ]
    }
   ],
   "source": [
    "# Replace the values in the dictionary below\n",
    "a = 662 \n",
    "b = 574 \n",
    "c = 20 \n",
    "d = 0 \n",
    "\n",
    "\n",
    "sol_4_dict = {\n",
    "    'How many users can we make predictions for in the test set?': c, \n",
    "    'How many users in the test set are we not able to make predictions for because of the cold start problem?': a, \n",
    "    'How many movies can we make predictions for in the test set?': b,\n",
    "    'How many movies in the test set are we not able to make predictions for because of the cold start problem?': d\n",
    "}\n",
    "\n",
    "# This should be article not movies. it was bugging me and wasted some time on it\n",
    "t.sol_4_test(sol_4_dict)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`5.` Now use the **user_item_train** dataset from above to find U, S, and V transpose using SVD. Then find the subset of rows in the **user_item_test** dataset that you can predict using this matrix decomposition with different numbers of latent features to see how many features makes sense to keep based on the accuracy on the test data. This will require combining what was done in questions `2` - `4`.\n",
    "\n",
    "Use the cells below to explore how well SVD works towards making predictions for recommendations on the test data.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 118,
   "metadata": {},
   "outputs": [],
   "source": [
    "# fit SVD on the user_item_train matrix\n",
    "u_train, s_train, vt_train = np.linalg.svd(user_item_train) # fit svd similar to above then use the cells below"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 119,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Use these cells to see how well you can use the training \n",
    "# decomposition to predict on test data\n",
    "def svd_algorithm(u_train, s_train, vt_train):\n",
    "    \"\"\" Return the results of the svd algorithm.\n",
    "    \n",
    "        Args:\n",
    "            u_train (np.array): user item interaction matrix  \n",
    "            s_train (np.array): sigma matrix\n",
    "            vt_train (np.array): v transpose matrix\n",
    "            \n",
    "        Returns:\n",
    "            Dataframe: dataframe of user article interaction counts sorted in descending order\n",
    "        \n",
    "        \"\"\"\n",
    "\n",
    "    num_latent_feats = np.arange(10,700+10,20)\n",
    "    sum_errs_train = []\n",
    "    sum_errs_test = []\n",
    "    all_errs = []\n",
    "\n",
    "    for k in num_latent_feats:\n",
    "        \n",
    "        row_idxs = user_item_train.index.isin(test_idx)\n",
    "        col_idxs = user_item_train.columns.isin(test_arts)\n",
    "        u_test = u_train[row_idxs, :]\n",
    "        vt_test = vt_train[:, col_idxs]\n",
    "        \n",
    "        # Split data \n",
    "        s_train_lat, u_train_lat, vt_train_lat = np.diag(s_train[:k]), u_train[:, :k], vt_train[:k, :]\n",
    "        u_test_lat, vt_test_lat = u_test[:, :k], vt_test[:k, :]\n",
    "\n",
    "        # Dot product:\n",
    "        user_item_train_preds = np.around(np.dot(np.dot(u_train_lat, s_train_lat), vt_train_lat))\n",
    "        user_item_test_preds = np.around(np.dot(np.dot(u_test_lat, s_train_lat), vt_test_lat))\n",
    "        all_errs.append(1 - ((np.sum(user_item_test_preds)+np.sum(np.sum(user_item_test))) \\\n",
    "                             /(user_item_test.shape[0]*user_item_test.shape[1])))\n",
    "\n",
    "        # Calculate the error of each prediction\n",
    "        diffs_train = np.subtract(user_item_train, user_item_train_preds)\n",
    "        diffs_test = np.subtract(user_item_test, user_item_test_preds)\n",
    "\n",
    "        # Get total Error\n",
    "        err_train = np.sum(np.sum(np.abs(diffs_train)))\n",
    "        err_test = np.sum(np.sum(np.abs(diffs_test)))\n",
    "\n",
    "        sum_errs_train.append(err_train)\n",
    "        sum_errs_test.append(err_test)\n",
    "        \n",
    "    # Plot accuracy for train and test vs number of latent features\n",
    "    plt.figure(figsize=(15,10))\n",
    "\n",
    "    # Latent features and training\n",
    "    plt.plot(num_latent_feats, 1 - np.array(sum_errs_train)/(user_item_train.shape[0]*user_item_test.shape[1]), label='Train', color='darkred')\n",
    "\n",
    "    # Latent features and testing\n",
    "    plt.plot(num_latent_feats, 1 - np.array(sum_errs_test)/(user_item_test.shape[0]*user_item_test.shape[1]), label='Test', color='darkblue')\n",
    "\n",
    "    plt.plot(num_latent_feats, all_errs, label='Total Error', color = \"orange\")\n",
    "    plt.xlabel('Number of Latent Features')\n",
    "    plt.ylabel('Accuracy')\n",
    "    plt.legend();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 120,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4gAAAJQCAYAAADMqko4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3XeYVdW9//H3mkYdigOMICgoJEAEEcfeSxDs9VqjsaFJjC16LYkNxfYziRqNigaNscWr14heFUs0UbENCKggiooRQem9TFu/P85hnBkGGJUze8r79Tz7OXuvtfY+302KfFzr7B1ijEiSJEmSlJV0AZIkSZKkhsGAKEmSJEkCDIiSJEmSpDQDoiRJkiQJMCBKkiRJktIMiJIkSZIkwIAoSZIkSUozIEqSJEmSAAOiJEmSJCktJ+kC6kOnTp1iz549ky5DkiRJkhIxfvz4eTHGzhsa1ywCYs+ePSkuLk66DEmSJElKRAjhi7qMc4mpJEmSJAkwIEqSJEmS0gyIkiRJkiSgmfwGsTalpaXMnDmTVatWJV1Kk9OyZUu6d+9Obm5u0qVIkiRJ+g6abUCcOXMm+fn59OzZkxBC0uU0GTFG5s+fz8yZM+nVq1fS5UiSJEn6DprtEtNVq1ZRUFBgONzIQggUFBQ4MytJkiQ1Qs02IAKGwwzxz1WSJElqnJp1QJQkSZIkfcuAmJD58+czaNAgBg0axKabbspmm21WeVxSUlKna5xyyilMmzYtw5VKkiRJai6a7UNqklZQUMDEiRMBuOqqq2jbti0XXnhhtTExRmKMZGXVnuPvu+++jNcpSZIkqflwBrGBmT59OltvvTVnnXUWgwcPZvbs2QwfPpyioiJ+8pOfMGLEiMqxu+22GxMnTqSsrIwOHTpwySWXsM0227DzzjszZ86cBO9CkiRJUmPkDCLwz/POY056Nm9j6TJoEPvccsv3OnfKlCncd9993HXXXQDccMMNbLLJJpSVlbH33ntz1FFH0b9//2rnLF68mD333JMbbriBCy64gNGjR3PJJZf84PuQJEmS1Hw4g9gAbbXVVmy//faVx4888giDBw9m8ODBTJ06lSlTpqx1TqtWrRg2bBgA2223HTNmzKivciVJkiQ1Ec4gwvee6cuUNm3aVO5/8skn3Hrrrbzzzjt06NCBE088sdZ3DObl5VXuZ2dnU1ZWVi+1SpIkSWo6nEFs4JYsWUJ+fj7t2rVj9uzZjB07NumSJEmSJDVRziA2cIMHD6Z///5svfXWbLnlluy6665JlyRJkiSpiQoxxqRryLiioqJYXFxcrW3q1Kn069cvoYqaPv98JUmSpIYjhDA+xli0oXEuMZUkSZIkAQZESZIkSVJaRgNiCGF0CGFOCOGDdfSHEMJtIYTpIYTJIYTBVfpODiF8kt5OrtK+XQjh/fQ5t4UQQibvQZIkSZKai0zPIN4PDF1P/zCgT3obDtwJEELYBLgS2BHYAbgyhNAxfc6d6bFrzlvf9SVJkiRJdZTRgBhj/DewYD1DDgUeiClvAR1CCF2B/YEXY4wLYowLgReBoem+djHGN2Pq6ToPAIdl8h4kSZIkqblI+jUXmwFfVjmemW5bX/vMWtolSZJUT2KMxIoKKsrKiOXlxPJyKsrLqx+XlVFRZb+29qp9FeXlxIqKVHtFRbX9ivJySH9WG1Nl7DrHVNlI111zf839rK+fGmMrx8QI6S3W8bNm21rXqvFnXeWg1v11jllHX2VbzeOq569j7HrPX3MPNfvq2L/Wfi33s857W7tz3X0ZNuD009n+N79J7Pt/qKQDYm2/H4zfo33tC4cwnNRSVDbffPPvW1/GzJ8/n3333ReAr7/+muzsbDp37gzAO++8Q15eXp2uM3r0aA444AA23XTTjNUqSVJzFmOkorSU8tJSKkpLUwFnzf6a9hpta47L19dWVkYsK/s2MK3ZXxOYatlfX1/N/arhqmqgqhagqo6p0bbOsWVlqUDTmIVAyMpKben9qm01+6kxds0YQkj11+FzXX2V31dzXI16v90N36l9XX2VbbV8Vu0LtY2tpa3acS1t69qvHFfL/lrXqc16+pJ6VEnbrl0T+d6NJemAOBPoUeW4OzAr3b5XjfZX0+3daxm/lhjjKGAUpN6DuLEK3lgKCgqYOHEiAFdddRVt27blwgsv/M7XGT16NIMHDzYgSpIatBgj5SUlVJSUUF5lqygtXat9rTElJalQVVt7zTE1g1st+3U9jmVlqc/y8nr/8wrZ2WTl5FR+1txfX1+1c/PyyMnKSu1nZ6fCTXZ2asvKSrVVba+lLSs7G9a0Z2VVXrva91TZr6ylxpgNjataY9Xvq62+WvvXMYYQKj8rA5+kdUo6II4Bzg4hPErqgTSLY4yzQwhjgeuqPJhmCHBpjHFBCGFpCGEn4G3gJOBPiVSeQX/961+54447KCkpYZddduH222+noqKCU045hYkTJxJjZPjw4RQWFjJx4kSOOeYYWrVq9Z1mHiVJzc+akFa6bBklS5dSsmwZpenPkqVLU+3p/bIVK74NY6tXfxvGVq+u/lmlv6y2cWv6S0szc1MhkJ2XR1Zu7refublkVdlqHue2aUOLdfRn5+YScnLWf42cnHW2VWuv0bbWddf0Z2cT1gS/NWHJECMpIRkNiCGER0jNBHYKIcwk9WTSXIAY413As8ABwHRgBXBKum9BCOEa4N30pUbEGNc87OYXpJ6O2gp4Lr39IOed908mTpzzQy9TzaBBXbjlln2+83kffPABTz75JOPGjSMnJ4fhw4fz6KOPstVWWzFv3jzef/99ABYtWkSHDh3405/+xO23386gQYM2av2SpGTFigpKV6ygdPny1LZs2bf7y5enwl26vVrAq/JZsmzZWm0VZWV1riErN5fsFi3IzstLbWv2059Z6f28/Pz19lfrS29Z69iverwm9K13THZ2Bv9TkKTmJ6MBMcZ43Ab6I/CrdfSNBkbX0l4MbL1RCmyAXnrpJd59912KiooAWLlyJT169GD//fdn2rRpnHvuuRxwwAEMGTIk4UolSVVn5KqFtqqBrkpbrf219S1fTtmKFd+plpxWrcht25a8tm3Jy88nt21bWnbsSH6PHuTl55PXtm2qf81++nPN2Gr7+fnktGyZWuonSWpWkl5i2iB8n5m+TIkxcuqpp3LNNdes1Td58mSee+45brvtNp544glGjRqVQIWS1HiVl5amZtaWLKFkyRJWpz8r95cuXSvAldQMejUC33f5fVrIzia3TZtUQGvTpnJr0b49bTfbrFpbbo0xtZ2X26ZNKtS1aUNWjv9IlyT9cP7TpIHZb7/9OOqoozj33HPp1KkT8+fPZ/ny5bRq1YqWLVty9NFH06tXL8466ywA8vPzWbp0acJVS1JmVZSXU7JkCasWLmT1okVrB7s67petXFmn78tu0YK8tm3JqRrK2ral7WabVTteK7jV+KwW6Nq2JTsvz9+WSZIaNANiAzNgwACuvPJK9ttvPyoqKsjNzeWuu+4iOzub0047jRgjIQRuvPFGAE455RROP/10H1IjqcErW72a1QsXVoa8VTX2Kz8XLmTVokXVPlcvWbLBd1pl5eTQon178tq1Sy2jbNeO1oWFdOzTJ9XWrh0t0p9V9yvb8vOdjZMkNXthvS+YbCKKiopicXFxtbapU6fSr1+/hCpq+vzzlZquGCOly5axYu5cVs6dy4p581g5dy4r581Lta05nj//2xm/hQspW7VqvdfNad2alh070qJDB1p27Fhtv/KzY0daduhAXvv2a4W97BYtnJ2TJGkdQgjjY4xFGxrnvyKVpGauoqyMlfPnVw93NcNejePykpJar5Wdl0erTp1o1bkzrQoKKOjff52hb83+mtCX7QoISZISZ0CUpCaqbNUqls2ezfLZs1k2a9Y6P1ctWLDOa7Ro3z4V9jp1ot3mm1M4eHDlcevOnWnVuTOt1wTCTp3Iy893Fk+SpEbMgChJjUzpypUsryX4LZs9m+VVPlctXLjWuVk5ObTZdFPadOtGh622YrPddqN1YWEq7K0JfVVmAJ3VkySpeTEgSlIDEisqWP711yz+/HMWz5jB4s8/Z8mMGSz54otUCJw1i9WLFq11XlZuLm26dqVt165s8qMf0WPPPWnbrVuqrcpnq4IC320nSZLWyYAoSfUoxsjKefOqB8Cq+198Qfnq1dXOabPpprTr2ZNN+vZl8332WSv0tenalVabbGLwkyRJP5gBUZI2slULF1ab/Vv8+eeVgXDJjBmULl9ebXyrggLa9epF54ED6X3oobTv1Yt2PXumPrfYgtxWrRK6E0mS1NwYEBMyf/589t13XwC+/vprsrOz6dy5M0Ct7zNcsGABjz32GGedddZ6r1tWVkanTp1YVGMJWllZGS1atGDAgAGVbSeccAIXXXTRxrgdqVmJMbLim29YOH06i6ZPZ+Enn7Aovb/o009ZvXhxtfF57drRvlcvOvbuTc+f/rRaAGzfsyd5+fkJ3YkkSVJ1BsSEFBQUMHHiRACuuuoq2rZty4UXXrjO8QsWLOCuu+7aYEBcn/z8/MrvXJeysjJyqrwguuZxXc+TGrsYI8u//rpaAFw4fTqLPvmEhdOnU7psWeXYkJ1N+5496dCnD1132ikV/NJbu549admxo0/2lCRJjYJ/o2+AbrrpJh544AEAzjzzTH79619zySWXMG3aNAYNGsTQoUO57LLLOOyww1i0aBFlZWVcd911HHTQQd/r+7p3786ZZ57J888/z3nnncett97KnnvuyWuvvcYRRxzBoYceyqmnnsr8+fMpLCzkvvvuo3v37px44okUFhYyYcIEtt9+e2666aaN+ccgZVyMkeWzZ1cLflVnA6suBc3KyaF9r1506N2bzXbfnY59+tChd2869ulDuy22IDs3N8E7kSRJ2jgMiADjz4OF659Z+846DoLtbvnOp73zzjs89NBDvPPOO5SXl7PDDjuw5557csMNNzB9+vTKGcDS0lKeeuop8vPzmTNnDrvuuusGA+LSpUsZNGhQ5fHvfvc7jjrqKADatGnDG2+8AcCtt97KkiVL+Pe//w3AsGHDOP300znhhBMYNWoU5513Ho8//jgAn376KS+//DJZPhxDDVjJ8uUsnDaN+VOnsmDqVOZ/9FFlICxbsaJyXFZODu233JIOvXvTY6+9KgNgx969yd98c0OgJElq8gyIDcxrr73GkUceSevWrQE47LDDeP311xkyZEi1cTFGLr74Yl5//XWysrL48ssvmTdvHh06dFjntde3xPSYY46pdnzsscdW7r/99ts888wzAJx00klcfvnllX1HH3204VANxop581IBsEoQXDB1Kku++KJyTMjOpsOWW9Lxxz9m8332qQyBHXr3pt3mm5PlUmlJktSM+Tch+F4zfZkSY6zTuAceeIDFixczYcIEcnJy6N69O6tWrfre39umTZv1Htf1PCnTYkUFS2fO/DYEVvlcOW9e5bicVq3YpG9fuu2yCwNOO41N+vWjoF8/OvTuTU6LFgnegSRJUsNlQGxg9thjD84880wuuugiysvLeeqpp/j73/9Ofn4+S5curRy3ePFiunTpQk5ODi+++CJfffVVxmraaaedeOyxxzjuuON48MEH2WOPPTL2XdIaFWVlLPzkk7WC4IJp06r9NrDlJptQ0K8fvQ87jIJ+/SqDYLvNN/e9gJIkSd+RAbGB2WGHHTjuuOPYfvvtAfjFL35R+WqKoqIiBgwYwIEHHsgFF1zAwQcfTFFREYMHD6ZPnz4bvHbN3yAeeOCBjBw5coPn3X777Zx22mlcf/31lQ+pkTamstWrmffBB8yZMIFvJkxgznvvMXfSJMqqzIrn9+hBQb9+DNh998oguEnfvrTu3NknhEqSJG0koa5LGhuzoqKiWFxcXK1t6tSp9OvXL6GKmj7/fLUuJcuWMXfSpMog+M2ECcz/8EMqysoAaNG+PV223ZYugwfTZZttKOjfn0369iWvbduEK5ckSWq8QgjjY4xFGxrnDKKkjFm5YAFz3nuvMgjOmTCBBR9/DOl/MdWqc2cKt9uOLQ84gC6DB1M4eDDte/VyRlCSJCkhBkRJG8Wy2bOrBcFv3nuPJTNmVPbn9+hB4eDB9D3+eArTM4Rtu3UzDEqSJDUgzTogxhj9y2kGNIdly81dyfLlfP3uu8x+802+GjeOb4qLWf7115X9Hfv0oeuOO7LNWWdROHgwXbbdltadOiVYsSRJkuqi2QbEli1bMn/+fAoKCgyJG1GMkfnz59OyZcukS9FGEmNkyRdfMOvNN5k1bhyz3nyTORMnEsvLAdikb1+2+OlPKdxuOwoHD6bzNtvQol27hKuWJEnS99FsA2L37t2ZOXMmc+fOTbqUJqdly5Z079496TL0PZWtXs2cCRP4Kh0GZ40bx/LZswHIbdOGrjvuyI6XXEK3XXah60470WqTTRKuWJIkSRtLsw2Iubm59OrVK+kypMQtmz3729nBceP4Zvx4yktKAGjfqxeb77MP3XbZhW4770znAQPIymm2/7chSZLU5Pk3PakZqSgrY+7kyXw1blzl7wfXPEgmu0ULNi0qYttzzmGzdCBss+mmyRYsSZKkemVAlJqwGCNzJ03i8+eeY8aLLzL77bcpW7ECgLbdutFtl10YfM45dNt5Z7psuy05LVokXLEkSZKSZECUmphVCxcy48UXU6Hw+ecrny7aZdAgBp5+Ot123pluu+xCfo8ePqBJkiRJ1RgQpUYuVlTwzXvv8flzz/H5c88x+623iBUVtOzYkS2GDKHX0KH03H9/2nbtmnSpkiRJauAMiFIjtHL+fGa88EJqlnDsWFbMmQNAYVERO152Gb2GDaPrDjv4QBlJkiR9J/7tUWoEKsrL+Wb8+G9nCd95B2KkVUFBapZw2DB67r8/bbp0SbpUSZIkNWIGRKmBWj5nDl9UmSVcOX8+hEDXHXZg5yuuYMthwygsKiIrOzvpUiVJktREGBClBiLGyJz33uOTf/yDz597jm/Gj0/NEnbuTK8DDqDX0KFsMWQIrTt1SrpUSZIkNVEGRClhiz7/nI8efpgpDz3EgqlTCVlZdN1pJ3YdMYJeQ4dSOHgwISsr6TIlSZLUDBgQpQSsnD+faf/zP0x98EG+euMNALrvvjvb3XUXPzrqKFoVFCRcoSRJkpojA6JUT0pXruSzZ55hyoMP8vlzz1FRWkpB//7sft119D3uONr37Jl0iZIkSWrmDIhSBlWUl/Plq68y5cEH+eSJJyhZupS23box+Jxz6H/iiXTeZhtfVi9JkqQGw4AobWQxRuZMnMjUhx7io0ceYdmsWeTl5/Ojo46i3wkn0GOvvXzyqCRJkhokA6K0kSyeMYOpDz/M1IceYv6UKWTl5tJr2DD2PvFEtjzoIHJbtUq6REmSJGm9DIjSD7BywQI+/p//YcqDD/LV668DsNluu7HfnXfy46OP9mEzkiRJalQMiNJ3VF5ayqdjxjDlb3/js2efpaK0lE369WO3kSPpd/zxPmxGkiRJjZYBUaqjRZ9/zvv33MP7o0ez4ptvaNO1K9v++tf0P/FEugwa5MNmJEmS1OgZEKX1qCgr49NnnmHS3XczY+xYQghsedBBbHPmmfTcf38fNiNJkqQmxYAo1WLJl1/y/r338v5f/sKyr76ibbdu7HzFFQw47TTa9eiRdHmSJElSRhgQpbSK8nJmjB3LpLvu4rP/+z9ijPTcf3/2u+MOtjzwQLJy/J+LJEmSmjb/xqtmb9ns2bz/l78w+Z57WPqf/9C6sJAdLr6YAWecQYdevZIuT5IkSao3BkQ1S7Gigi9efplJd93Fp2PGUFFWxub77steN99M70MPJTsvL+kSJUmSpHpnQFSzsnzOHD68/34m3X03iz/7jFYFBQw+7zy2GT6cjn36JF2eJEmSlCgDopq8GCNf/utfTLrrLj753/+lorSU7nvswW7XXEOfI44gp2XLpEuUJEmSGgQDopqskqVLmXzvvUy++24WTJtGiw4dGPTLXzJw+HA69e+fdHmSJElSg2NAVJNTtno1k+++mzevvZaVc+fSbeedGXr//fz46KPJbd066fIkSZKkBsuAqCajorycqQ8/zBtXXMGSGTPosffe7HH99XTdccekS5MkSZIaBQOiGr0YI589+yyvXXop895/ny7bbsuQu+9mi5/+lBBC0uVJkiRJjUZWJi8eQhgaQpgWQpgeQriklv4tQggvhxAmhxBeDSF0r9J3Ywjhg/R2TJX2+0MIn4cQJqa3QZm8BzVsX40bx9/33JMnDzqIshUrOOiRR/hZcTE9hwwxHEqSJEnfUcZmEEMI2cAdwE+BmcC7IYQxMcYpVYbdDDwQY/xrCGEf4HrgZyGEA4HBwCCgBfCvEMJzMcYl6fMuijE+nqna1fDN+/BDXrvsMj4dM4bWhYXs9+c/M+D008nOzU26NEmSJKnRyuQS0x2A6THGzwBCCI8ChwJVA2J/4Pz0/ivAP6q0/yvGWAaUhRAmAUOBxzJYrxqBxV98wbgrr+TDBx4gLz+f3a69lsHnnUdemzZJlyZJkiQ1eplcYroZ8GWV45nptqomAUem9w8H8kMIBen2YSGE1iGETsDeQI8q541ML0v9YwihRW1fHkIYHkIoDiEUz507d2PcjxK0Yt48XrngAkb/6Ed89OijFF1wAWd89hk7/fa3hkNJkiRpI8nkDGJtPwCLNY4vBG4PIfwc+DfwFVAWY3whhLA9MA6YC7wJlKXPuRT4GsgDRgEXAyPW+qIYR6X7KSoqqvm9aiRKli1j/C238O7/+3+ULlvGT04+mV2uuop2m2+edGmSJElSk5PJgDiT6rN+3YFZVQfEGGcBRwCEENoCR8YYF6f7RgIj030PA5+k22enT18dQriPVMhUE1NeUsLke+7hzWuuYcU339D7sMPYbeRIX3AvSZIkZVAmA+K7QJ8QQi9SM4PHAsdXHZBeProgxlhBamZwdLo9G+gQY5wfQhgIDAReSPd1jTHODqlHVB4GfJDBe1A9ixUVfPToo7x++eUs/uwzuu+xB4c9+STddt456dIkSZKkJi9jATHGWBZCOBsYC2QDo2OMH4YQRgDFMcYxwF7A9SGESGqJ6a/Sp+cCr6VfU7AEODH9wBqAh0IInUktYZ0InJWpe1D9iTEyY+xYXrv0UuZMnEjnbbbhiGefpdfQob6uQpIkSaonIcam//O8oqKiWFxcnHQZWoevx4/n1d/8hpn/+hfte/Vit2uvpe+xxxKyMvqaTkmSJKnZCCGMjzEWbWhcJpeYSutVtno1b159Ne/ceCOtOnVinz/9iW2GDyc7Ly/p0iRJkqRmyYCoRHwzYQLPnXwy8z74gK1PPZW9//AHWrRvn3RZkiRJUrNmQFS9Ki8p4a2RI3lr5Ehad+nC4c88w1YHHph0WZIkSZIwIKoezZ08medOPpk5EyfS/8QT2ee222jZsWPSZUmSJElKMyAq4yrKynjnxhsZd/XVtOzYkUOffJI+hx2WdFmSJEmSajAgKqPmTZnCcyefzDfFxfz4mGPY9/bbad2pU9JlSZIkSaqFAVEZUVFeTvEf/sAbl19OXtu2HPzYY/z46KOTLkuSJEnSehgQtdEt+Phjnv/5z5n15pv0Ofxw9rvzTtoUFiZdliRJkqQNMCBqo4kVFUy47TZeu/RSclq14sCHHqLvcccRQki6NEmSJEl1YEDURrHo0095/pRTmPnaa2x50EEMGTWKtl27Jl2WJEmSpO/AgKgfJFZUMPGuu/jXRReRnZvL0Pvv5ycnneSsoSRJktQIGRD1vS2eMYOxp53Gf/75T3oOGcL+f/kL+d27J12WJEmSpO/JgKjvLMbI+/feyysXXADAkFGjGHD66c4aSpIkSY2cAVHfydKZMxl7+unMGDuWHnvvzdDRo2nfs2fSZUmSJEnaCAyIqpMYIx8+8ACvnHsu5aWl7Hv77Qz6xS8IWVlJlyZJkiRpIzEgaoNKli/npbPOYsqDD9J9990Zet99dNhqq6TLkiRJkrSRGRC1XvOmTOHpo49m/tSp7DpiBDtedhlZ2dlJlyVJkiQpAwyIWqcpDz3EC8OHk9e2LUe/+CJb7Ltv0iVJkiRJyiADotZStmoVr5x3HpPuvpvuu+/OQY8+Sttu3ZIuS5IkSVKGGRBVzaLPPuPpo4/mmwkT2P6//5vdR44kK8f/mkiSJEnNgX/zV6VP/vEPnv/5zyEEDnvqKXofckjSJUmSJEmqR76jQJSXlvLqhRfy1OGH06FPH342YYLhUJIkSWqGnEFs5pbOnMkzxx7LV2+8waBf/pK9/vAHclq0SLosSZIkSQkwIDZjM158kf87/njKVq7koEceoe+xxyZdkiRJkqQEucS0GaooL+eNq67i8f33p3VhIScWFxsOJUmSJDmD2NwsnzOHZ084gS9eeon+J53Efn/+M3lt2iRdliRJkqQGwIDYjMx8/XWeOeYYVi1YwP733svWp55KCCHpsiRJkiQ1EC4xbQZijLx78838fa+9yGndmuPfeosBp51mOJQkSZJUjTOITdyqRYt4/uc/Z/pTT9HnyCMZ+pe/0KJ9+6TLkiRJktQAGRCbsK/Hj+fpo49m6ZdfsvcttzD4nHOcNZQkSZK0TgbEJijGyKS77uKV886jdWEhx772Gt122inpsiRJkiQ1cAbEJqZk2TJePPNMpj78ML2GDWPYAw/QulOnpMuSJEmS1AgYEJuQirIyxhx5JF+89BK7XXstO156KSHL5xBJkiRJqhsDYhPy6m9+w4wXXmD/e+9lwGmnJV2OJEmSpEbG6aUmYtKoUUy47Ta2O/98w6EkSZKk78WA2AT859VXeflXv6Ln0KHsedNNSZcjSZIkqZEyIDZyiz77jDFHHkmH3r05+NFHycpx1bAkSZKk78eA2IitXrKEJw8+GGLk8KefpkX79kmXJEmSJKkRc7qpkaooL+eZ445j4ccfc9TYsXTs3TvpkiRJkiQ1cgbERurfF1/M588+y3533snm++yTdDmSJEmSmgCXmDZC7993H8W//z2DfvUrBp11VtLlSJIkSWoiDIiNzMzXX+fFM89ki/32Y59bbkm6HEmSJElNiAGxEVk8YwZPHXEE7XuMrS7dAAAgAElEQVT25ODHHvOJpZIkSZI2KgNiI1GydClPHnIIFSUlHP7007Ts2DHpkiRJkiQ1MU5BNQKxooL/O/FE5k+ZwpHPPssmP/5x0iVJkiRJaoIMiI3Aa7/9LZ+OGcM+t91GzyFDki5HkiRJUhPlEtMGbsqDD/LODTcwcPhwtj377KTLkSRJktSEGRAbsFlvvcXY00+nx157se/ttxNCSLokSZIkSU2YAbGBWvLll/zjsMNou9lmHPL442Tn5iZdkiRJkqQmzt8gNkAly5fzj0MPpWzFCv7rn/+kVUFB0iVJkiRJagYMiA1MrKjguZNPZu6kSRz+9NN06t8/6ZIkSZIkNRMGxAZm3NVX88kTT7DX73/PlgcckHQ5kiRJkpqRjP4GMYQwNIQwLYQwPYRwSS39W4QQXg4hTA4hvBpC6F6l78YQwgfp7Zgq7b1CCG+HED4JIfw9hJCXyXuoTx/9/e+8OWIEW59yCtudf37S5UiSJElqZjIWEEMI2cAdwDCgP3BcCKHmesmbgQdijAOBEcD16XMPBAYDg4AdgYtCCO3S59wI/DHG2AdYCJyWqXuoT18XF/P8z3/OZrvtxn533ukTSyVJkiTVu0zOIO4ATI8xfhZjLAEeBQ6tMaY/8HJ6/5Uq/f2Bf8UYy2KMy4FJwNCQSk37AI+nx/0VOCyD91Avln71Ff849FBaFxZyyBNPkNOiRdIlSZIkSWqGMhkQNwO+rHI8M91W1STgyPT+4UB+CKEg3T4shNA6hNAJ2BvoARQAi2KMZeu5ZqNSunIl/zjsMFYvXszhY8bQpkuXpEuSJEmS1ExlMiDWtkYy1ji+ENgzhPAesCfwFVAWY3wBeBYYBzwCvAmU1fGaqS8PYXgIoTiEUDx37tzveQuZFWNk7Kmn8s348Rz48MN0Hjgw6ZIkSZIkNWOZDIgzSc36rdEdmFV1QIxxVozxiBjjtsBv022L058jY4yDYow/JRUMPwHmAR1CCDnrumaVa4+KMRbFGIs6d+68Me9ro3lr5Eg+evRRdr/uOnofckjS5UiSJElq5jIZEN8F+qSfOpoHHAuMqToghNAphLCmhkuB0en27PRSU0IIA4GBwAsxxkjqt4pHpc85GXgqg/eQMR8/8QRvXH45/U88kR0uvjjpciRJkiQpcwEx/TvBs4GxwFTgsRjjhyGEESGENdNlewHTQggfA4XAyHR7LvBaCGEKMAo4scrvDi8GLgghTCf1m8S/ZOoeMqV0xQpe+tWv6LrTTgy55x6fWCpJkiSpQQipSbmmraioKBYXFyddRjVz33+f1p0702bTTZMuRZIkSVITF0IYH2Ms2tC4nA0NUGZ0HjAg6RIkSZIkqZpM/gZRkiRJktSIGBAlSZIkSYABUZIkSZKUZkCUJEmSJAEGREmSJElSmgFRkiRJkgQYECVJkiRJaQZESZIkSRJgQJQkSZIkpRkQJUmSJEmAAVGSJEmSlGZAlCRJkiQBBkRJkiRJUpoBUZIkSZIEGBAlSZIkSWkGREmSJEkSYECUJEmSJKUZECVJkiRJgAFRkiRJkpRmQJQkSZIkAQZESZIkSVKaAVGSJEmSBBgQJUmSJElpBkRJkiRJEmBAlCRJkiSlGRAlSZIkSYABUZIkSZKUZkCUJEmSJAEGREmSJElSmgFRkiRJkgQYECVJkiRJaQZESZIkSRJgQJQkSZIkpRkQJUmSJEmAAVGSJEmSlGZAlCRJkiQBBkRJkiRJUpoBUZIkSZIEGBAlSZIkSWkGREmSJEkSYECUJEmSJKUZECVJkiRJgAFRkiRJkpRmQJQkSZIkAQZESZIkSVKaAVGSJEmSBBgQJUmSJElpBkRJkiRJEmBAlCRJkiSlGRAlSZIkSYABUZIkSZKUZkCUJEmSJAEZDoghhKEhhGkhhOkhhEtq6d8ihPByCGFyCOHVEEL3Kn03hRA+DCFMDSHcFkII6fZX09ecmN66ZPIeJEmSJKm5yFhADCFkA3cAw4D+wHEhhP41ht0MPBBjHAiMAK5Pn7sLsCswENga2B7Ys8p5J8QYB6W3OZm6B0mSJElqTjI5g7gDMD3G+FmMsQR4FDi0xpj+wMvp/Veq9EegJZAHtABygW8yWKskSZIkNXuZDIibAV9WOZ6ZbqtqEnBkev9wID+EUBBjfJNUYJyd3sbGGKdWOe++9PLSy9csPa0phDA8hFAcQiieO3fuxrgfSZIkSWrSMhkQawtuscbxhcCeIYT3SC0h/QooCyH0BvoB3UmFyn1CCHukzzkhxjgA2D29/ay2L48xjooxFsUYizp37vzD70aSJEmSmrhMBsSZQI8qx92BWVUHxBhnxRiPiDFuC/w23baY1GziWzHGZTHGZcBzwE7p/q/Sn0uBh0ktZZUkSZIk/UCZDIjvAn1CCL1CCHnAscCYqgNCCJ1CCGtquBQYnd7/D6mZxZwQQi6p2cWp6eNO6XNzgYOADzJ4D5IkSZLUbGQsIMYYy4CzgbHAVOCxGOOHIYQRIYRD0sP2AqaFED4GCoGR6fbHgU+B90n9TnFSjPFpUg+sGRtCmAxMJLUk9Z5M3YMkSZIkNSchxpo/C2x6ioqKYnFxcdJlSJIkSVIiQgjjY4xFGxqXySWmkiRJkqRGxIAoSZIkSQIMiJIkSZKkNAOiJEmSJAkwIEqSJEmS0gyIkiRJkiTAgChJkiRJSjMgSpIkSZIAA6IkSZIkKc2AKEmSJEkCDIiSJEmSpDQDoiRJkiQJMCBKkiRJktIMiJIkSZIkwIAoSZIkSUozIEqSJEmSAAOiJEmSJCnNgChJkiRJAgyIkiRJkqQ0A6IkSZIkCTAgSpIkSZLSDIiSJEmSJMCAKEmSJElKMyBKkiRJkgADoiRJkiQpzYAoSZIkSQIMiJIkSZKkNAOiJEmSJAkwIEqSJEmS0gyIkiRJkiTAgChJkiRJSjMgSpIkSZIAA6IkSZIkKc2AKEmSJEkCDIiSJEmSpDQDoiRJkiQJqENADCGcHULoWB/FSJIkSZKSU5cZxE2Bd0MIj4UQhoYQQqaLkiRJkiTVvw0GxBjj74A+wF+AnwOfhBCuCyFsleHaJEmSJEn1qE6/QYwxRuDr9FYGdAQeDyHclMHaJEmSJEn1KGdDA0II5wAnA/OAe4GLYoylIYQs4BPgvzNboiRJkiSpPmwwIAKdgCNijF9UbYwxVoQQDspMWZIkSZKk+laXJabPAgvWHIQQ8kMIOwLEGKdmqjBJkiRJUv2qS0C8E1hW5Xh5uk2SJEmS1ITUJSCG9ENqgNTSUuq2NFWSJEmS1IjUJSB+FkI4J4SQm97OBT7LdGGSJEmSpPpVl4B4FrAL8BUwE9gRGJ7JoiRJkiRJ9W+DS0VjjHOAY+uhFkmSJElSguryHsSWwGnAT4CWa9pjjKdmsC5JkiRJUj2ryxLTvwGbAvsD/wK6A0szWZQkSZIkqf7VJSD2jjFeDiyPMf4VOBAYkNmyJEmSJEn1rS4BsTT9uSiEsDXQHuiZsYokSZIkSYmoS0AcFULoCPwOGANMAW6sy8VDCENDCNNCCNNDCJfU0r9FCOHlEMLkEMKrIYTuVfpuCiF8GEKYGkK4LYQQ0u3bhRDeT1+zsl2SJEmS9MOsNyCGELKAJTHGhTHGf8cYt4wxdokx3r2hC4cQsoE7gGFAf+C4EEL/GsNuBh6IMQ4ERgDXp8/dBdgVGAhsDWwP7Jk+505Sr9nok96G1ulOJUmSJEnrtd6AGGOsAM7+ntfeAZgeY/wsxlgCPAocWmNMf+Dl9P4rVfojqSem5gEtgFzgmxBCV6BdjPHNGGMEHgAO+571SZIkSZKqqMsS0xdDCBeGEHqEEDZZs9XhvM2AL6scz0y3VTUJODK9fziQH0IoiDG+SSowzk5vY2OMU9Pnz9zANSVJkiRJ38MG34MIrHnf4a+qtEVgyw2cV9tvA2ON4wuB20MIPwf+DXwFlIUQegP9SL1SA1IhdQ9gZR2umfryEIaTWorK5ptvvoFSJUmSJEkbDIgxxl7f89ozgR5VjrsDs2pcexZwBEAIoS1wZIxxcTrcvRVjXJbuew7YidQ7Gbuv75pVrj0KGAVQVFRUa4iUJEmSJH1rgwExhHBSbe0xxgc2cOq7QJ8QQi9SM4PHAsfXuHYnYEH6t46XAqPTXf8BzgghXE9qJnJP4JYY4+wQwtIQwk7A28BJwJ82dA+SJEmSpA2ryxLT7avstwT2BSaQekDMOsUYy0IIZwNjgWxgdIzxwxDCCKA4xjgG2Au4PoQQSS0xXbOM9XFgH+B9UktIn48xPp3u+wVwP9AKeC69SZIkSZJ+oJB6GOh3OCGE9sDfYoyHZKakja+oqCgWFxcnXYYkSZIkJSKEMD7GWLShcXV5imlNK0i9f1CSJEmS1ITU5TeIT/Ptk0KzSL278LFMFiVJkiRJqn91+Q3izVX2y4AvYowz1zVYkiRJktQ41SUg/geYHWNcBRBCaBVC6BljnJHRyiRJkiRJ9aouv0H8H6CiynF5uk2SJEmS1ITUJSDmxBhL1hyk9/MyV5IkSZIkKQl1CYhzQwiVr7QIIRwKzMtcSZIkSZKkJNTlN4hnAQ+FEG5PH88ETspcSZIkSZKkJGwwIMYYPwV2CiG0BUKMcWnmy5IkSZIk1bcNLjENIVwXQugQY1wWY1waQugYQri2PoqTJEmSJNWfuvwGcViMcdGagxjjQuCAzJUkSZIkSUpCXQJidgihxZqDEEIroMV6xkuSJEmSGqG6PKTmQeDlEMJ96eNTgL9mriRJkiRJUhLq8pCam0IIk4H9gAA8D2yR6cIkSZIkSfWrLktMAb4GKoAjgX2BqRmrSJIkSZKUiHXOIIYQfgQcCxwHzAf+Tuo1F3vXU22SJEmSpHq0viWmHwGvAQfHGKcDhBDOr5eqJEmSJEn1bn1LTI8ktbT0lRDCPSGEfUn9BlGSJEmS1AStMyDGGJ+MMR4D9AVeBc4HCkMId4YQhtRTfZIkSZKkerLBh9TEGJfHGB+KMR4EdAcmApdkvDJJkiRJUr2q61NMAYgxLogx3h1j3CdTBUmSJEmSkvGdAqIkSZIkqekyIEqSJEmSAAOiJEmSJCnNgChJkiRJAgyIiYoxJl2CJEmSJFUyICZg1aoyfvazZ7n99veSLkWSJEmSKhkQE9CiRTbz5q3gsste48svlyRdjiRJkiQBBsREhBD485/3o6Ii8stfvuRSU0mSJEkNggExIb16dWDEiF155pnPeOKJj5MuR5IkSZIMiEk699ztGDy4kF//+p8sXLgq6XIkSZIkNXMGxATl5GRxzz1DmDNnBZdc8u+ky5EkSZLUzBkQEzZ4cCHnn78do0ZN5rXXZiZdjiRJkqRmzIDYAFx99S707NmO4cNfYPXqsqTLkSRJktRMGRCTsvxLWP4FAG3a5HHnnT/lo48WcP31bydcmCRJkqTmyoCYhPLVMLYIxp9f2TR0aC+OP74f1133NlOnzk+wOEmSJEnNlQExCdktoM/ZMPNJmPdWZfMf/7gX+fl5DB/+AhUVvhtRkiRJUv0yICal7/nQshAmXgwxFQa7dGnD73+/F6+//hX33DM54QIlSZIkNTcGxKTktoWtr4A5/4ZZz1U2n3zyT9hnn8357//+F7NmLUuwQEmSJEnNjQExSb3PgLZbwaRLIVYAEELg7rt/SklJBeee+8+EC5QkSZLUnBgQk5SVCwOvhUWTYcbDlc29e3fkiit25vHHP2bMmOkJFihJkiSpOTEgJm2L/4KOg2Hy71JPN0278MIitt66E7/85UssWbJ6PReQJEmSpI3DgJi0kAWDbki9E/GTuyqbc3Ozuffe/Zk1axm/+93rCRYoSZIkqbkwIDYEXX8KhfvCh9dC6ZLK5h137MrZZ2/L7be/x1tvzUqwQEmSJEnNgQGxoRh0A6yeB1NvrtY8cuTubLZZPmec8QKlpeUJFSdJkiSpOTAgNhQFRbD5f8FHf4CV31Q25+fncccd+/LBB/O4+ebiBAuUJEmS1NQZEBuSgdemHlTzwTXVmg85pDdHHfUjrr56HJ98sjCh4iRJkiQ1dQbEhqRdn9S7EaffDUurv97ittv2oWXLHM488wVijAkVKEmSJKkpMyA2NFtfAVl5MPnyas1du7blxhv34JVXvuSvf/0woeIkSZIkNWUGxIam1abQ93z44lFYMKFa1xlnDGS33TbjN795lTlzlidUoCRJkqSmyoDYEPW7CFoUwMRLqjVnZQVGjRrC0qUlnH/+q8nUJkmSJKnJMiA2RHnt4Se/ha9fhK9frtbVr18Bl122Iw8/PJXnn/88oQIlSZIkNUUZDYghhKEhhGkhhOkhhEtq6d8ihPByCGFyCOHVEEL3dPveIYSJVbZVIYTD0n33hxA+r9I3KJP3kJg+v4DWm8PEiyFWVOu69NId6dt3E37xixdZvrwkoQIlSZIkNTUZC4ghhGzgDmAY0B84LoTQv8awm4EHYowDgRHA9QAxxldijINijIOAfYAVwAtVzrtoTX+McWKm7iFR2S1h4DWwYDz85/FqXS1a5DBq1BBmzFjClVeOS6hASZIkSU1NJmcQdwCmxxg/izGWAI8Ch9YY0x9Ys4bylVr6AY4CnosxrshYpQ1VzxOg/dYw6bdQUVqta/fduzN8+ED++MfxTJjwTUIFSpIkSWpKMhkQNwO+rHI8M91W1STgyPT+4UB+CKGgxphjgUdqtI1ML0v9YwihxcYquMHJyoZB18Oy6fDpX9bqvvHGPejSpTVnnPECZWUVtVxAkiRJkuoukwEx1NJW8w3vFwJ7hhDeA/YEvgLKKi8QQldgADC2yjmXAn2B7YFNgItr/fIQhocQikMIxXPnzv3eN5G4bgdC593h/auhrPqrLTp0aMmf/rQPEyZ8w623jk+oQEmSJElNRSYD4kygR5Xj7sCsqgNijLNijEfEGLcFfptuW1xlyH8BT8YYS6ucMzumrAbuI7WUdS0xxlExxqIYY1Hnzp03zh0lIQQYdCOs+ho+umWt7iOP/BEHH7wVV1zxBp9/viiBAiVJkiQ1FZkMiO8CfUIIvUIIeaSWio6pOiCE0CmEsKaGS4HRNa5xHDWWl6ZnFQkhBOAw4IMM1N6wdN4Zuh8KU2+CVfOqdYUQuOOOfcnKCvzyly8RY81JWkmSJEmqm4wFxBhjGXA2qeWhU4HHYowfhhBGhBAOSQ/bC5gWQvgYKARGrjk/hNCT1Azkv2pc+qEQwvvA+0An4NpM3UODss11ULYMPrxura4ePdpx3XW78/zzM3jkkY8SKE6SJElSUxCaw4xTUVFRLC4uTrqMH+6t02DGg3Dwx9Bmi2pd5eUV7LrrI3z22SKmTj2VgoJWCRUpSZIkqaEJIYyPMRZtaFwml5hqYxtwFRBg8pVrdWVnZzFq1BAWLlzN+ee/4lJTSZIkSd+ZAbExadMDfnwOfP4ALHp/re6BAztzySU78Le/TeHoo8ewZMnqBIqUJEmS1FgZEBub/pdAbnuYeFmt3SNG7MpNN+3BP/4xnaKiB5k8uRG/4kOSJElSvTIgNjYtNoH+F8OsZ2DOa2t1hxC46KIdeOWVY1i2rIQdd3yI++9v+g96lSRJkvTDGRAbox+fA626wcSLYR2/Ndx99+68995J7LJLN0455XlOO+15Vq4srXWsJEmSJIEBsXHKaZ16YM28N+GrMescVljYhhdeOIrf/W4nRo/+gJ13fpjp0xfWX52SJEmSGhUDYmO15SnQ7scw6TKoKFvnsOzsLK65Zjf+7/+O4Msvl7Lddn/jf//343osVJIkSVJjYUBsrLJyYJvrYPGU1FNNN+CAA7ZkwoSf0bfvJhx55Bh+85tXKC0tr4dCJUmSJDUWBsTGrPvhULAjvH8llK3c4PAttmjPa68dx69/vS1/+MN49trr78ycubQeCpUkSZLUGBgQG7MQYNANsGImfHJHnU7Jy8vmttv25dFHD2Ly5Llsu+0DvPjijMzWKUmSJKlRMCA2doV7Qddh8OF1ULKozqcdc0xf3n33RAoLW7P//o8zYsQ4KipqfyKqJEmSpObBgNgUDLo+FQ6n3PidTuvbt4C33z6BE0/sz5VXjuOAA55g3rwVGSpSkiRJUkNnQGwKOm4DPY+HabfCiq++06lt2uTx178OY9SoIbz66pdsu+3fePPNWRkqVJIkSVJDZkBsKgZeA7EM3r/6O58aQuCMMwYybtzx5OZm/f/27jw+qur+//jrk4QECPsOIoR93xQBERDZKmKFKm6gItFqrW2/9lfrWq1a11pr9Vu/tS4BXIDibhGQXRQFAVkSQPZ9DTshhGzn98e9gSEsssxkZpL38/GYx8zce+fcz+Ro8nhzzj2XHj3G8sorC3BOU05FREREREoSBcTiolwDaHwPrE2BAyvOqYmLLqrJggW3ctVVDbjvvhnccMN/OXDgSJALFRERERGRSKWAWJy0fhRiy8Cih+AcR/8qVy7Np58O4q9/7cEnn6yiY8f3WLIkPciFioiIiIhIJFJALE5K14BWj8DmT2FuMuRln1MzZsYf/9iJGTNuJCMjm86d32fkyLQgFysiIiIiIpFGAbG4afkQtP4zrB0JM/uf1a0vCuvevS4LF95G1651GD58Enfe+SWHD+cEr1YREREREYkoCojFjRm0fQIufQfSv4bJXSFj3Tk3V7NmIpMnD+ZPf+rC22+n0rv3B6Sn61YYIiIiIiLFkQJicdXgVrhiCmRthy87w66559xUbGwMf/lLNz766BoWLtxJly7v8+OPu4NYrIiIiIiIRAIFxOKs5uXQ7zsoVR6m9YSNH51Xc9de25SZM28kIyOHSy8dzcyZG4NTp4iIiIiIRAQFxOKuQjPoNwcqd4BvBsOyF895hVOAzp1rM3fuUGrXTqRfvw8ZNUqL14iIiIiIFBcKiCVB6erQezrUuwEWPQDzfgX5577YTFJSRb79dgjdu9fl9tsn8ec/z8adR+gUEREREZHIoIBYUsSWhsvGQMuHYfUbMPNqyDlwzs1VqlSaiROvIzm5NU899R233DKBI0dyg1iwiIiIiIgUNQXEksRioP2z0Plt2DEdpnSDQ+d+HWF8fCxvvfUznn22O6NHL6dv3w/ZvftwEAsWEREREZGipIBYEjVKhismwqEN3gqnexacc1NmxsMPd2bs2Kv5/vttXHrpaFat2hvEYkVEREREpKgoIJZUtfpA328hNgGm9IDNn59Xczfe2Jxp025gz54sunR5n2++2RykQkVEREREpKgoIJZklVpBv7lQsRXMGgQ/vnJeK5xedtkFzJkzhGrVytC79weMGbM8iMWKiIiIiEioKSCWdGVqQp+ZUHcQ/HAfLPgd5J/7YjONG1fmu++G0KVLbYYM+YKnn/5OK5yKiIiIiEQJBUSBuLLQ/UNocT+s/Kc3mpiTcc7NValShsmTB3PLLS157LHZJCdPIjs7L4gFi4iIiIhIKCggisdioMOLcMm/YNskmNodMrecc3MJCXG8805/nniiKyNHLuXKKz9k796sIBYsIiIiIiLBpoAox2vyK7h8PBxc7a1wunfxOTdlZvz5z115553+fPPNFrp2Hc3atfuCWKyIiIiIiASTAqKcqM6V0PcbMPPulbhlwnk1d+utrZgy5Xp27MikS5f3mTNna5AKFRERERGRYFJAlJOr3M5b4bR8E5j1c1j1r/Nq7vLLL+S774ZQvnw8V1wxjg8+WBGkQkVEREREJFgUEOXUytaBPrOgzgCY92v44Q+Qf+6LzTRrVoU5c4Zy0UU1uOGG//LXv36vFU5FRERERCKIAqKcXqly0P0TaPo7+PHv8P1d4PLPubnq1csybdoN3HhjMx58cBZ33z2FnBytcCoiIiIiEgniwl2ARIGYWOj4CsRXgrSnIK4cXPwP7xrFc1C6dByjR19No0aVePbZuaxfv58PPriGihUTgly4iIiIiIicDY0gyplr8wQ0/3+w8lVY8qfzaiomxnjmme68/fbPmDFjE927j2Hr1nO/96KIiIiIiJw/BUQ5c2bQ4W/Q+C5Y+iwsfe68m0xObsOECdeybt1+unYdzcqVe4JQqIiIiIiInAsFRDk7ZtDx/yBpKCx+BFb873k32bdvEjNm3EhmZg6XXTaGefO2BaFQERERERE5WwqIcvZiYqHLSKg7CBb8DtaMOO8mO3asxezZQyhXrhRXXDGOyZPXn3ebIiIiIiJydhQQ5dzExMFlY6FWP/j+Ttgw7rybbNKkMt9+O4RGjSoxYMDHjB69PAiFioiIiIjImVJAlHMXmwA9PoFql8G3Q2HL+PNusnbtcsyadROXXVaHoUO/4JVXFgShUBERERERORMKiHJ+4spCz/FQuT18PRi2TzvvJitWTGDSpMFce20T7rtvBg8/PAvnXBCKFRERERGR01FAlPNXqgJcMQnKN4FZAyH92/NusnTpOMaN+zl33dWW55//njvu+JLc3PwgFCsiIiIiIqeigCjBkVAVek2BMnVg5lWw54fzbjI2NobXX+/L449fyogRaVx77WdkZuYEoVgRERERETkZBUQJnjK1oNdUKFURZvwM9i877ybNjCefvIzXXuvN+PFr6NfvQ/buzQpCsSIiIiIiUpgCogRXYj3oPQ0sDqb3gYNrgtLsr3/dgf/85+fMm7ed7t3HsGXLwaC0KyIiIiIixyggSvCVb+xNN83Phum94dCmoDR7/fXNmDjxOjZuPEjXrqP58cfdQWlXREREREQ8CogSGpVawxVfQvZebyTx8I6gNNurVz1mzryRrKw8unUby9y524LSroiIiIiIKCBKKFW5GC7/AjI3w4x+cGRPUJq96KKazJ59MxUrxtOr13+YNGldUNoVERERESnpFBAltGp0g8s/gwM/wsz+kBOcawcbN67M7NlDaNq0Cj//+Se89975L4gjIiIiIlLSKSBK6NXqA93GwZ4F8NXVkJsZnGZrJfLVVzfSvfsF3HrrBP7+9/lBaVdEREREpKQKafxSqg0AACAASURBVEA0syvNbIWZrTazh06yv76ZTTOzJWY208zq+tuvMLNFAY8sMxvk72tgZnPNbJWZ/cfM4kP5HSRI6g6ES9+FnV/D19dB3pGgNFuhQgITJlzH4MFN+cMfZvLAA1/hnAtK2yIiIiIiJU3IAqKZxQKvAf2BlsDNZtay0GF/A95xzrUFngKeA3DOzXDOtXfOtQd6AZnAZP8zLwAvO+eaAHuBO0L1HSTIkm6GTm/Atknw7RDIzw1Ks6VLxzF27NXcc087XnxxHsOHTyInJy8obYuIiIiIlCShHEHsBKx2zq11zmUDY4GBhY5pCUzzX884yX6AwcBE51ymmRleYPzQ3zcKGBT0yiV0Gt8JF70Mmz6GOcng8oPSbGxsDK+91ocnn+zKqFFLGTToUw4dyg5K2yIiIiIiJUUoA+IFQOAN8Db72wItBq7zX/8CKG9mVQsdcxMwxn9dFdjnnCsYejpZmwCY2V1mNt/M5qenp5/jV5CQaH4ftP0LrH8X5t0LQZoSamY8/nhXXn+9L5MmradPnw/YvftwUNoWERERESkJQhkQ7STbCieB+4HLzWwhcDmwBTg679DMagNtgC/Pok1vo3NvOOc6Ouc6Vq9e/Wxrl1Br9Si0fBBWvw6LHghaSAS4++52fPDBz/nhh5107z6GTZsOBK1tEREREZHiLJQBcTNwYcD7usDWwAOcc1udc9c65zoAj/rb9gcccgPwiXMux3+/C6hkZnGnalOihBm0ew6a3AvL/wZLHg9qSLz22qZ8+eV1bNmSQdeuY1i+fHfQ2hYRERERKa5CGRDnAU38VUfj8aaKfh54gJlVM7OCGh4GUgq1cTPHppfivOUpZ+BdlwgwDPgsBLVLUTCDjq9Cw2RY+jTM/y3kB29xmZ496zFr1k3k5ubTrdsYvvtO/5YgIiIiInI6IQuI/nWCv8GbHrocGOecW2pmT5nZNf5hPYEVZrYSqAk8U/B5M0vCG4H8qlDTDwL/z8xW412T+HaovoMUAYuBzm9Ci/th1Wvw7c1BuwUGQLt2NZg9+2aqVClN797j+OKLNUFrW0RERESkuLGScM+4jh07uvnzdRP1iLf8b7Dwj1CzF/T4BEpVCFrTO3ce4qqrPmbRop2kpFzJbbe1ClrbIiIiIiKRzswWOOc6/tRxoZxiKnJ2WtwPl74DO2fB1J5weEfQmq5RI5EZM26kZ88LGTZsIi+++H3Q2hYRERERKS4UECWyNLgVLv8cDqyAKZfBweBNCS1fPp4vvriWG25oxgMPzOL++2eSn1/8R9BFRERERM6UAqJEnjr9ofc0yN7rhcQ9C4PWdEJCHGPGXM1vftOBl16az7BhE8nJCd7COCIiIiIi0UwBUSJTtS7Q9xuIiYepl8OOGUFrOibGePXVXjz9dDfee28ZAwd+yqFD2UFrX0REREQkWikgSuSq2AL6fQuJF8KMK2Hjh0Fr2sx49NEuvPFGP778cj29e3/A7t2Hg9a+iIiIiEg0UkCUyFa2LvT5Gqp0hG9ugFWvB7X5X/6yLR99dA2LFu2kW7cxbNx4IKjti4iIiIhEEwVEiXwJVaDXFKgzAObdA0uegCDenmXQoCZMnjyYbdsO0bXraJYu3RW0tkVEREREookCokSHuLLevREbDoe0J2HeryE/eIvL9OhxIbNm3UR+vqNbtzHMnr0laG2LiIiIiEQLBUSJHjFx0PltaPkQrH4dZt8AeVlBa75t2+rMnn0z1auXpU+fDxg/Pni32BARERERiQYKiBJdzKD9c3DRy7DpY5jRH7L3B635Bg0qMXv2zbRuXZVBgz5l5Mi0oLUtIiIiIhLpFBAlOjW/D7q+D+nfeLfBOLwtaE1Xr16W6dNvpFevegwfPokXXpiLC+I1jyIiIiIikUoBUaJX0hC4fDxkrIbJl8HB1UFrunz5eMaPv5abb27OQw99zR/+MJP8fIVEERERESneFBAlutX5GfSaDrkHYHJX2LMgaE3Hx8fy3nsD+N3vLuLllxdw660TyM4O3sI4IiIiIiKRRgFRol+1TtB3NsSWgak9Yfu0oDUdE2P84x9X8Nxz3Rk9ejnXXPMJGRnZQWtfRERERCSSKCBK8VChGfT7FhKTYOZVsGFc0Jo2Mx56qDNvv/0zpkzZQK9e40hPzwxa+yIiIiIikUIBUYqPshdA31lQtRPMvglWvhbU5pOT2/DJJwNJTd1Ft25jWL16b1DbFxEREREJNwVEKV7iK8MVk+GCn8P838DixyCIK5Bec01jpkwZzM6dmbRpM4pnn52j6xJFREREpNhQQJTiJ64MdP8IGt0BS5+GRQ8ENSR261aX1NTbueqqBjz66Dd06PAO33yzOWjti4iIiIiEiwKiFE8xcdDpTWhyLyz/Gyx6MKghsW7d8nz00UD++99fkJGRTffuY7nzzi/Zvftw0M4hIiIiIlLUFBCl+DKDjv/rh8QXYdFDQQ2JAFdf3Yhly4bzxz9ewsiRaTRvnsI77yzFBfk8IiIiIiJFQQFRirejIfHXsPyvsPjhoIfExMR4/vrXy/nhh9to3LgSw4ZNpHfvcaxYsSeo5xERERERCTUFRCn+zKDjP6HJPbDsBVj8SNBDIkDbttWZPXsIr7/el4ULd9K27Sj+/OfZZGXlBv1cIiIiIiKhoIAoJUNBSGz8K1j2PCx+NCQhMSbGuPvudvz4YzKDBzflqae+o23bUUybtiHo5xIRERERCTYFRCk5LAYueQ0a3w3LnoMlfwpJSASoWTOR998fwOTJg3HO0afPB9xyyxfs3HkoJOcTEREREQkGBUQpWSwGLvk/aHwXLH0WlgT3PomF9e2bRGrq7Tz2WBfGjVtBs2YpvPHGYvLztYiNiIiIiEQeBUQpeSwGLvkXNPolLH0m5CGxdOk4nnqqG0uWDKNdu+rcffcUuncfQ2pqesjOKSIiIiJyLhQQpWSyGOj0OjS60w+Jj4c0JAI0b16VGTNuZOTIK1mxYi8XXfQuDz00i8zMnJCeV0RERETkTCkgSsllMdDp335IfBpS/xzykGhmDBvWmhUrkrnttpa88ML3tGo1ggkT1ob0vCIiIiIiZ0IBUUq2oyHxDkj7C6Q+USSnrVq1DG+/fSVffXUjZcrEMWDAx1x//eds3ZpRJOcXERERETkZBUQRi4FOb0DDZEh7CpY8UWSn7tHjQhYtGsYzz3Rj/Pi1NG+ewj/+sYAjR3TvRBEREREpegqIIuCFxM5vQsPhkPYkpD5ZZKeOj4/lkUe6kJZ2O1271uH3v59B8+YpjBqVRl5efpHVISIiIiKigChSwGKg81vQ8HZvqmnqU0V6+kaNKjFx4nVMnjyYatXKcPvtk2jTZiQff7wSF+JrI0VEREREQAFR5HgWA53eggbDvEVrUv9StKc3o2/fJL7//hY++ugaAK677nM6dXqPKVPWKyiKiIiISEgpIIoUFhMLnd+GBrdB6uOQ9nSRl2BmXHttU1JTb2fEiCvZuTOTfv0+pHfvccyZs7XI6xERERGRkkEBUeRkYmKhc4oXEpc8BmnPhKWM2NgYbr+9NStX3sGrr/Zi6dLdXHrpaAYO/ITU1PSw1CQiIiIixZcCosipFITEpFthyZ9g6bNhKyUhIY7f/vYi1qy5k2ee6cZXX22mXbtR3HLLF6xZsy9sdYmIiIhI8aKAKHI6MbHQZYQXEhc/CkufC2s55crF88gjXVi79k4efLATH3+8iubNU7jnnim6h6KIiIiInDcFRJGfcjQk3gKLH4Glz4e7IqpUKcNzz/VgzZo7ueuutrz1ViqNG7/Fgw9+xe7dh8NdnoiIiIhEKQVEkTMREwtdRkLSUFj8MCx7IdwVAVC7djlee60PK1YkM3hwU158cR4NG77J009/R0ZGdrjLExEREZEoo4AocqZiYqHLKKg/BBY9BPN/CzkHw10VAA0bVuKdd64iNfV2evWqx2OPzaZhwzd55ZUFZGXlhrs8EREREYkSCogiZyMmFi4dBU1/Bytfg/EtYNPHECH3J2zVqhqffDKIuXOH0rZtde67bwZNm75NSkoqubn54S5PRERERCKcAqLI2YqJg46vQL/vIKEafH0dfHUNHNoQ7sqO6tSpNlOn3sDUqddTu3Yid9zxJa1bj+Tdd5cqKIqIiIjIKSkgipyrap3hyvnQ4W+wYzqMbwnLXoT8nHBXdlTv3vWZM2con3wykPj4GG67bSJNm77N668v0tRTERERETmBAqLI+YiJgxZ/gKuXQ60+sOgBmNQRds0Jd2VHmRmDBjVh0aJhfP75L6hRoyz33DOVhg3f5KWX5mkxGxERERE5SgFRJBgS68Hln0H3TyB7D0zuCt/fA9mRcxP7mBjj5z9vxHffDWHatBto2bIq99//FfXrv8GTT37Lnj26PYaIiIhISaeAKBJMFw6CAcug2X2w5g0Y3xzWj4mYRWzAG1Hs1aseU6fewJw5Q+nW7QKeeOJb6td/gwce+Irt2w+Fu0QRERERCRMFRJFgK1UeLv47/Gw+lK0H3w6BGT+Dg6vDXdkJOneuzWef/YIlS4ZxzTWNeeml+SQlvcG9905l/fr94S5PRERERIqYAqJIqFTp4K102vGfsHsufNEa0p6GvCPhruwEbdpU5/33B7BiRTK33daKN99cQuPGbzFs2AR+/HF3uMsTERERkSJiLoKmvoVKx44d3fz588NdhpRkmVvhh9/DxnFQoTlc8jrUvDzcVZ3S5s0Heeml+fz734vJysrl2mub8MgjXbjooprhLk1EREREzoGZLXDOdfyp4zSCKFIUytaBbv+BnhO9EcRpPWHOcMjaFe7KTqpu3fK8/PIVbNhwF48+2oWpUzdy8cXv0r//h3z99eZwlyciIiIiIaKAKFKU6lwJA9Kg5cOw7j34ojmsGRFRi9gEql69LH/5Szc2bLiL557rzoIFO+jRYyzdu49h0qR1lIQZCCIiIiIliQKiSFGLKwvtn4X+i6BCC5ib7I0o7l8e7spOqWLFBB56qDPr19/Fq6/2YsOGA/Tv/xEXX/wuH364gvx8BUURERGR4iCkAdHMrjSzFWa22sweOsn++mY2zcyWmNlMM6sbsK+emU02s+VmtszMkvztI81snZkt8h/tQ/kdREKmUivo8xV0fgv2pcHEdrD4T5AbufcjLFu2FL/97UWsXn0nKSk/IyMjh+uv/y+NGr3JU099y4YNWvlUREREJJqFbJEaM4sFVgJ9gc3APOBm59yygGM+AMY750aZWS9guHPuVn/fTOAZ59wUMysH5DvnMs1spP+ZD8+0Fi1SIxEvKx0W3g/r3oHEJGjzJCQNhZjYcFd2Wnl5+Xz88SreeGMJ06ZtAKB37/okJ7dm0KDGlClTKswVioiIiAhExiI1nYDVzrm1zrlsYCwwsNAxLYFp/usZBfvNrCUQ55ybAuCcy3DOZYawVpHwKl0dLh0FvadDfBWYM8wbUdz0acRenwgQGxvD9dc3Y8qU61m37pc88URXVq/ey5AhX1Cnzuvce+9UFizYrmsVRURERKJEKAPiBcCmgPeb/W2BFgPX+a9/AZQ3s6pAU2CfmX1sZgvN7EV/RLLAM/601JfNLCFUX0CkyNW8Aq6cD90+AJcLX/8CJneB7dN++rNhVr9+RR5/vCtr1vySadNuYMCAhqSkpNGx43u0b/8Or7yygF279O88IiIiIpEslAHRTrKt8DDC/cDlZrYQuBzYAuQCcUB3f/8lQEPgdv8zDwPN/e1VgAdPenKzu8xsvpnNT09PP79vIlKUzKDeYLgqDTqnwOHtML0PTOsDu+aGu7qfFBNj9OpVj/feG8C2bb/iX//qQ0JCLPfdN4M6dV5n8ODPmDBhLbm5+eEuVUREREQKCeU1iJcCTzjnfua/fxjAOffcKY4vB/zonKtrZl2A551zPf19twJdnHP3FvpMT+B+59zVp6tF1yBKVMs7Aqv/DWlPw5F0qDsI2v4FKrUOd2VnJTU1nREj0nj33WXs2nWYOnXKMWxYK4YPb02TJpXDXZ6IiIhIsRYJ1yDOA5qYWQMziwduAj4PPMDMqplZQQ0PAykBn61sZtX9972AZf5navvPBgwC0kL4HUTCLzYBmv0OrlkLbZ+GHdNhQlv49jbIWBvu6s5YmzbV+fvfr2DLll/x0UfX0KFDDV544XuaNn2bHj3GMnJkGhkZ2eEuU0RERKREC9kIIoCZXQX8A4gFUpxzz5jZU8B859znZjYYeA5v6uks4F7n3BH/s32Bl/Cmqi4A7nLOZZvZdKC6v30R8CvnXMbp6tAIohQrR/bAshdg5avg8qDRL6H1n6BM7XBXdta2bs3g3XeXkpKSxsqVeylXrhQ33tic4cNb07VrHbx/BxIRERGR83WmI4ghDYiRQgFRiqXMrbD0aVj9JsSU8kYZWzwACVXCXdlZc87x7bdbSUlJ5T//WcGhQzk0bVqZ5OTW3HZbK2rXLhfuEkVERESimgJiAAVEKdYOroHUJ2D9+1CqArT4IzT7HygVnaEqIyObDz5YQUpKGt98s4XYWKN//wYkJ7dhwICGxMdH9r0hRURERCKRAmIABUQpEfalwpLHYPNnULoGtHoUGt/tXcMYpVau3MOIEWmMGrWUbdsOUb16GW69tSXDh7emdevqP92AiIiIiAAKiMdRQJQSZdccWPwI7JgBifWhzROQdAvExIW7snOWm5vP5MnrSUlJ5fPP15CTk0+nTrVITm7DTTc1p2LF6A3BIiIiIkVBATGAAqKUOM7Bjmmw6BHYMw8qNPdWQL3wF2ChXLw49NLTM3n//eWkpKSSmrqL0qXjuO66JiQnt6FnzwuJidHCNiIiIiKFKSAGUECUEss52PwpLPkT7F8GFVtDq4eh3g1RPaII3sI2CxbsICUlldGjf2T//iMkJVVg+PDWDBvWivr1K4a7RBEREZGIoYAYQAFRSrz8PNgwFpY9B/uXQrlG0PJBaHBbVF+jWODw4Rw+/XQ1KSlpTJ26ATPo06c+ycltGDSoMaVLR3cYFhERETlfCogBFBBFfC4ftvwX0p7xpp6WqQMt7ofGd0FcYrirC4r16/czatRSRoxIY8OGA1SqlMCQIS1ITm7NRRfV1L0VRUREpERSQAyggChSSME1imnPwM6ZkFAVmt0HTX8D8ZXCXV1Q5Oc7Zs7cREpKKh99tIqsrFzatq1OcnJrhg5tQbVqZcNdooiIiEiRUUAMoIAochrp38LS52DreIgrD03v9cJimZrhrixo9u3LYuzYH0lJSWPevO2UKhVD+/Y16NCh4FGTNm2qUbZsqXCXKiIiIhISCogBFBBFzsDexV5Q3DjOuy6x0S+96aeJ9cJdWVClpqbz/vvL+f77bSxcuJN9+44AEBNjNGtWmQ4dagYExxpUqVImzBWLiIiInD8FxAAKiCJn4cBKWPYCrHvHe9/gVmj5EFRoGt66QsA5x8aNB1i4cOdxj82bDx49pl698ieExrp1y+taRhEREYkqCogBFBBFzsGhjbD8b7DmTcg7AvWuh1aPQOV24a4s5NLTM1m06PjQuHLlHgp+XVatWua4wNihQw2aNKlMbGx032NSREREii8FxAAKiCLnIWsn/PgPWPUa5ByAOgO8oFi9a7grK1IZGdksWZJ+XHBMTd1FdnYeAGXLxtG2bXUuvbQO/fs3oEePuiQk6PYaIiIiEhkUEAMoIIoEQfY+WPkarHgZjuyGGj29oFirD5TQ6ZY5OXksX76HhQt3HA2Nc+du48iRPMqWjaN37/r079+A/v0bkJRUMdzlioiISAmmgBhAAVEkiHIPweo3vemnh7dAlUug+f+DCwdBbOlwVxd2mZk5zJixkYkT1zFhwjrWrdsPQIsWVejfvwFXXdWQbt0u0OiiiIiIFCkFxAAKiCIhkHcE1r0Ly56HjDUQXxnqD4FGyVC5Q4kdVQzknGPlyr1+WFzLV19tJjs7j8TEUvTpc2x0sV69CuEuVURERIo5BcQACogiIeTyYcd0WDMCNn0E+UegUjsvKCYNhYSq4a4wYhw6lM306ZuYOHEtEyasY8OGAwC0alX16OjiZZddQHx8bJgrFRERkeJGATGAAqJIEcneCxvGwpoU2DMfYuKh7kBomAy1+kKMgk8B5xw//rjn6OjirFmbycnJp3z5+ONGF+vWLR/uUkVERKQYUEAMoIAoEgZ7l8DaEbD+XW9RmzIXQMNh0HA4lG8c7uoizsGD2Uyf7l27OHHiWjZu9O7F2KZNNfr3b8CVVzbg4otrUqFCQpgrFRERkWikgBhAAVEkjPKyYct/vbC4baI3JbVGD29Usd5giEsMd4URxznHsmW7j44ufv31FnJz8wFo3LiSf+/FmkfvwVizpn6GIiIicnoKiAEUEEUiROZWWPcOrE2Bg6sgrhzUv9ELi9Uu1cI2p3DgwBG++WaLfysN75Yaa9fuP7q/du1EOnSoQfv2NY6Gx4YNK2L6eYqIiIhPATGAAqJIhHEO0md7QXHjOO/WGRWaeUGxwa1Qpna4K4x4+/ZlsXhx+nGhcdmy3eTleb/TK1SIDwiMXmhs0aIKpUrpOlAREZGSSAExgAKiSATLyYCNH3hhMf0bsFioc5UXFi8YADGlwl1h1MjKyiUtbddxoXHJknQyM3MBSEiIpXXrascFx3btqpOYGB/mykVERCTUFBADKCCKRIkDK2DtSFg3Cg5vg4Tq3ohiw2So1Crc1UWlvLx8Vq7cy6JFO/3g6D127z4MeLN669evQJMmlWnSpDKNG1c6+rpBg4q65YaIiEgxoYAYQAFRJMrk58K2L2HN294CNy4XqnbygmL9myC+YrgrjGrOOTZvPng0LK5YsYdVq/ayatU+9u8/cvS4mBgLCI+VaNzYey4Ij5quKiIiEj0UEAMoIIpEsax0WP+ed2/F/WkQWxouHOzdLqNmT7CYcFdYbDjn2L37MKtW7WPVqr2sXr336OtVq/Zy4ED20WNjY+24kcfA0cekpAoKjyIiIhFGATGAAqJIMeAc7FngXau4fjTk7IfEJC8oNrwdEuuFu8JizTnHrl2Hj4bF1av3HR11XLVqLwcPHh8ek5Iq0rx5leOud2zQQCurioiIhIsCYgAFRJFiJvcwbP7EG1XcMQ0wqNXHm4J64SBvlFGKjHOO9PTMgJFH73nZst0sX35sZdWKFRNo3776cfdwbNGiKnFxGgUWEREJNQXEAAqIIsVYxnpvYZu1IyBzI5SqBElDoVEyVO6geyuG2eHDOQErq+48urLq4cPHVlZt06bacaGxbdvqlC2r1WtFRESCSQExgAKiSAng8mHHdG9UcdPHkH8EKrXzpqAmDYXS1cJdofhyc/NZuXLPcaFx4cKd7N2bBXiL4zRrVvm40NihQw2qVCkT5spFRESilwJiAAVEkRImey9sGOuFxT3zvXspXjDQG1Ws1Q9itIBKpHHOsXHjARYtSj96D8eFC3eyadPBo8fUq1ee9u29EcbExNCMMJYrV4ratctRq1YitWolUrt2okYzRUSkWFBADKCAKFKC7V3iTT9d/y4c2Q1lLoAGt3nXLFZuDwlVwl2hnMauXZknhMYVK/ZQlH+6ypePp3btY4Ex8Nl7XY5atcpSrVpZYmI0pVlERCKTAmIABUQRIS/bu6fi2hTYNsmbkgqQWN8LipU7HHuUratrFyNYTk7e0YVvgsk5x8GD2Wzbdojt271HwevCz4GrthaIjTVq1iwcHr3natXKhCw8VqlS+ui5KlcurZViRUTkpBQQAyggishxjuz2bpmxdyHsWeg9H1wF+L8PE6qdGBrLN9HUVDnq0KFstm/P9ENjxklD5Pbth9ixI5P8/KL7OxsfH0utWmUDRjZPHPWsXTuRmjUTiY/Xf88iIiWJAmIABUQR+Uk5GbBvsRcY9y3ynvenQb4/UhRbFiq3OxYYq3SAiq0hNiG8dUtEy8vLZ9euw+zefTgk02Lz8x179mQFhNMTw+quXYdP+tmqVctQq1bZUwbJihUTNJAuInIOatQoywUXlA93GSdQQAyggCgi5yQvGw4sP36kce8iyPUXTrE4qNjSD43+iGOVi6FUufDWLRIgJyePnTszTzJdNoPt2zOPC5VHjuSFu1wRkaj34IOdeP75HuEu4wQKiAEUEEUkaFw+ZKwtFBoXQtYOb39MKajeHer0hzpXQYUWup5RooJzjv37jxwNkAcOnHidpYiI/LQmTSrRsmXk3V5LATGAAqKIhNzhbV5g3DkTtk70pqcClK3nBcU6/aFmL40uioiISFgoIAZQQBSRIndoo7da6tYJsH0a5GZATDzU6AG1C0YXm2l0UURERIqEAmIABUQRCau8bEj/xguL2ybC/mXe9sSkgNHFKyAuMaxlioiISPGlgBhAAVFEIsqhDd401K0TYftUyMuEmASocfmxaxfLN9HoooiIiASNAmIABUQRiVh5RyD9az8wToADP3rbyzX0gmLt/lCzJ8SVDWuZIiIiEt0UEAMoIIpI1MhYd2x0ccd0b3QxtjRU7wFVLzl2D8bEBhphFBERkTOmgBhAAVFEolJeFuyc5YfFad61i86/T12pisfuvVgQGis0926zISIiIlLImQbEuKIoRkREzkFsaajdz3sA5B72bp8ReA/G1f+GvMPe/pgEqNTmWGCs3AEqtdX0VBERETljCogiItEirow3zbTqJce25efBwRWwd9Gx4LjpI1jzprffYqB80+NHGit3gISq4fkOIiIiEtEUEEVEollMLFRs6T2ShnjbnIPMTcePNKZ/AxvGHPtc2QuPhcaKLb1VU8s3hlLlw/M9REREJCIoIIqIFDdmkFjPe9QdeGx71i7Yt+hYaNy7ELb8Fwi4Fr10TT8s+oGx4HW5xlCqXJF/FRERESlaCogiIiVF6WpQq4/3KJB7CA6uhoOrAp5XwbZJsHZboc/XOjE4lm+s8CgiIlKMKCCKiJRkcYlQuZ33KCwnAzJWHx8cD672VlXNGnH8sWVqe0HxS6rOcwAAD5lJREFUaHBsAhWa+Sur6k+NiIhItNBfbREROblS5fxbabQ/cV/OQchYExAcC8LjF5C149hxWllVREQkqiggiojI2StV/jTh8QAcXOPdt7HgWsdNHxZaWbXZ8aGxcgdIqFK030FEREROoIAoIiLBVaqCF/yqdIAGQ71tzkHmxuMXyEmfBRtGH/tc2XrHB8bKHaBsXW/RHRERESkSIQ2IZnYl8AoQC7zlnHu+0P76QApQHdgD3OKc2+zvqwe8BVyIt8TeVc659WbWABgLVAF+AG51zmWH8nuIiMh5MoPE+t7jwkHHtmelH7uHY8Fj8+ccXVk1oerxgbFye+++jjGxYfkaIiIixZ055376qHNp2CwWWAn0BTYD84CbnXPLAo75ABjvnBtlZr2A4c65W/19M4FnnHNTzKwckO+cyzSzccDHzrmxZvY6sNg596/T1dKxY0c3f/78UHxNEREJtpwM2LfkWGDcsxD2p0G+/2+BMQkQWzo0544r5y24U7qW91ym1vHvS9fytoXq/CIiIiFiZguccx1/6rhQjiB2AlY759b6BY0FBgLLAo5pCfzefz0D+NQ/tiUQ55ybAuCcy/C3G9AL8O8GzSjgCeC0AVFERKJIqXJQvav3KJCXDQeWe4Fx/1LIzwn+eZ2D3INweBsc3gx75kHWTo67T+TRGiudIjz6obK0/xxfRVNkRUQkqoQyIF4AbAp4vxnoXOiYxcB1eNNQfwGUN7OqQFNgn5l9DDQApgIPAZWBfc653IA2LzjZyc3sLuAugHr16gXj+4iISLjExp/6dhyhlJ8LR9Lh8HYvOGYFPm+HrG2w+3tvW17miZ+PKeUFx4Rq3uI8oRBf+fQhtVRFhVQRETljoQyIJ/trVPifYe8H/mlmtwOzgC1Arl9Xd6ADsBH4D3A78PkZtOltdO4N4A3wppiedfUiIiIxcX7gqo33J+kUnIPcjONDZEGAPLwdjuzmFH+uzo/Lh+w9cGCld978k1ySH1v6xAB5sim0pWvqnpUiIhLSgLgZb4GZAnWBrYEHOOe2AtcC+NcZXuec229mm4GFAdNTPwW64C1oU8nM4vxRxBPaFBERKXJm3q0/SpWHCk3DU4NzkLPvWDg92YjnwRWw8ysvVJ7AvJHOgsAYX5GT/1uviIicVt2BkDTkp4+LUKEMiPOAJv6qo1uAmzh27SAAZlYN2OOcywcexguABZ+tbGbVnXPpeNcdznfOOTObAQzGW8l0GPBZCL+DiIhIdDDzppvGV4aKLU9/bN4RyNpxfHgsHCgzNxZN3SIixU3VTuGu4LyELCA653LN7DfAl3i3uUhxzi01s6fwwt7nQE/gOTNzeFNM7/U/m2dm9wPT/IVpFgD+HZZ5EBhrZk8DC4G3Q/UdREREiqXYBEis5z1EREQChOw2F5FEt7kQEREREZGS7ExvcxGiJdVEREREREQk2iggioiIiIiICKCAKCIiIiIiIj4FRBEREREREQEUEEVERERERMSngCgiIiIiIiKAAqKIiIiIiIj4FBBFREREREQEUEAUERERERERnwKiiIiIiIiIAAqIIiIiIiIi4lNAFBEREREREUABUURERERERHwKiCIiIiIiIgIoIIqIiIiIiIhPAVFEREREREQABUQRERERERHxKSCKiIiIiIgIoIAoIiIiIiIiPgVEERERERERARQQRURERERExKeAKCIiIiIiIgCYcy7cNYScmaUDG4rodNWAXUV0Lgk+9V/0Ut9FN/VfdFP/RTf1X3RT/0W3ouy/+s656j91UIkIiEXJzOY75zqGuw45N+q/6KW+i27qv+im/otu6r/opv6LbpHYf5piKiIiIiIiIoACooiIiIiIiPgUEIPvjXAXIOdF/Re91HfRTf0X3dR/0U39F93Uf9Et4vpP1yCKiIiIiIgIoBFEERERERER8SkgBomZXWlmK8xstZk9FO565ERmlmJmO80sLWBbFTObYmar/OfK/nYzs1f9/lxiZheFr3IBMLMLzWyGmS03s6Vm9j/+dvVhFDCz0mb2vZkt9vvvSX97AzOb6/fff8ws3t+e4L9f7e9PCmf9AmYWa2YLzWy8/159F0XMbL2ZpZrZIjOb72/T788oYGaVzOxDM/vR/xt4qfouOphZM///uYLHATO7L9L7TwExCMwsFngN6A+0BG42s5bhrUpOYiRwZaFtDwHTnHNNgGn+e/D6son/uAv4VxHVKKeWC/zBOdcC6ALc6/9/pj6MDkeAXs65dkB74Eoz6wK8ALzs999e4A7/+DuAvc65xsDL/nESXv8DLA94r76LPlc459oHLKmv35/R4RVgknOuOdAO7/9D9V0UcM6t8P+faw9cDGQCnxDh/aeAGBydgNXOubXOuWxgLDAwzDVJIc65WcCeQpsHAqP816OAQQHb33GeOUAlM6tdNJXKyTjntjnnfvBfH8T7A3kB6sOo4PdDhv+2lP9wQC/gQ3974f4r6NcPgd5mZkVUrhRiZnWBAcBb/ntDfVcc6PdnhDOzCkAP4G0A51y2c24f6rto1BtY45zbQIT3nwJicFwAbAp4v9nfJpGvpnNuG3gBBKjhb1efRjB/yloHYC7qw6jhT1FcBOwEpgBrgH3OuVz/kMA+Otp//v79QNWirVgC/AN4AMj331dFfRdtHDDZzBaY2V3+Nv3+jHwNgXRghD/F+y0zS0R9F41uAsb4ryO6/xQQg+Nk/zKq5WGjm/o0QplZOeAj4D7n3IHTHXqSberDMHLO5fnTbOrizbxocbLD/Gf1X4Qws6uBnc65BYGbT3Ko+i6yXeacuwhvCtu9ZtbjNMeqDyNHHHAR8C/nXAfgEMemI56M+i4C+ddoXwN88FOHnmRbkfefAmJwbAYuDHhfF9gaplrk7OwoGLr3n3f629WnEcjMSuGFw/edcx/7m9WHUcafHjUT71rSSmYW5+8K7KOj/efvr8iJU8SlaFwGXGNm6/EuoeiFN6Kovosizrmt/vNOvGugOqHfn9FgM7DZOTfXf/8hXmBU30WX/sAPzrkd/vuI7j8FxOCYBzTxV3SLxxtC/jzMNcmZ+RwY5r8eBnwWsP02fzWpLsD+gqkAEh7+NUxvA8udc38P2KU+jAJmVt3MKvmvywB98K4jnQEM9g8r3H8F/ToYmO50496wcM497Jyr65xLwvv7Nt05NxT1XdQws0QzK1/wGugHpKHfnxHPObcd2GRmzfxNvYFlqO+izc0cm14KEd5/pt/ZwWFmV+H9i2oskOKceybMJUkhZjYG6AlUA3YAfwY+BcYB9YCNwPXOuT1+GPkn3qqnmcBw59z8cNQtHjPrBnwNpHLsOqhH8K5DVB9GODNri3chfizeP06Oc849ZWYN8UalqgALgVucc0fMrDTwLt61pnuAm5xza8NTvRQws57A/c65q9V30cPvq0/8t3HAaOfcM2ZWFf3+jHhm1h5vgah4YC0wHP/3KOq7iGdmZfGuK2zonNvvb4vo//cUEEVERERERATQFFMRERERERHxKSCKiIiIiIgIoIAoIiIiIiIiPgVEERERERERARQQRURERERExKeAKCIiIWNmzsxeCnh/v5k9EaS2R5rZ4J8+8rzPc72ZLTezGYW2J5lZ2lm0M8jMWp5HHUlmNuQ0+w6b2aKAR3wwzyEiIiWDAqKIiITSEeBaM6sW7kICmVnsWRx+B/Br59wV53naQcA5B0QgCThdeFvjnGsf8MgOwTlO6ix/niIiEsEUEEVEJJRygTeA3xfeUXgE0Mwy/OeeZvaVmY0zs5Vm9ryZDTWz780s1cwaBTTTx8y+9o+72v98rJm9aGbzzGyJmd0d0O4MMxsNpJ6knpv99tPM7AV/2+NAN+B1M3vxTL6wmf3SP/diM/vIzMqaWVfgGuBFf3Svkf+YZGYL/O/QPODn8qqZfWtmawN+Rs8D3f3Pn/DzPEUtiWaW4tez0MwG+tuT/HP+4D+6nuwcZna7mf0zoL3xZtbTf51hZk+Z2VzgUjO72O+3BWb2pZnV9o/7nZkt8/ti7JnULSIi4RMX7gJERKTYew1YYmZ/PYvPtANaAHuAtcBbzrlOZvY/wG+B+/zjkoDLgUbADDNrDNwG7HfOXWJmCcBsM5vsH98JaO2cWxd4MjOrA7wAXAzsBSab2SDn3FNm1gu43zk3/wxr/9g596bf7tPAHc65/zWzz4HxzrkP/X3TgF8551aZWWfg/4Befhu18YJpc+Bz4EPgIb+Oq09x3kZmtsh/Pds5dy/wKDDdOZdsZpWA781sKrAT6OucyzKzJsAYoGPhc5jZ7af5nolAmnPucTMrBXwFDHTOpZvZjcAzQLLfZgPn3BG/BhERiWAKiCIiElLOuQNm9g7wO+DwGX5snnNuG4CZrQEKAl4qEDjVc5xzLh9YZWZr8QJVP6BtwMhbRaAJkA18Xzgc+i4BZjrn0v1zvg/0AD49w3oDtfaDYSWgHPBl4QPMrBzQFfjAzAo2JwQc8qn/vZaZWc0zPO8a51z7Qtv6AdeY2f3++9JAPWAr8E8zaw/kAU3P8ByB8oCP/NfNgNbAFP/7xALb/H1LgPfN7FPO7ecpIiJFSAFRRESKwj+AH4ARAdty8S91MC9VBC6qciTgdX7A+3yO/9vlCp3HAQb81jl3XDDzp0YeOkV9dort52IkMMg5t9gfget5kmNigH0nCXQFAr//+dRmwHXOuRXHbfQWCtqBN1IbA2Sd4vNH+8hXOuB1lnMuL+A8S51zl56kjQF4Yfsa4DEza+Wcyz3bLyIiIkVD1yCKiEjIOef2AOPwFnwpsB5vSifAQKDUOTR9vZnF+NclNgRW4I3Y3eNPe8TMmppZ4k+0Mxe43Myq+Quu3Iw3ZfJclAe2+ecfGrD9oL8P59wBYJ2ZXe/XaGbW7ifaPfr5s/Al8Fs/gGNmHfztFYFt/ijlrXgjfic7x3qgvf8zvhBviu7JrACqm9ml/nlKmVkrM4sBLnTOzQAe4NioqoiIRCgFRBERKSovAYGrmb6JF8q+Bzpz6tG901mBF+Qm4l3PlwW8BSwDfjDvNhT/5idmzPjTWR8GZgCLgR+cc5+dwfmbmdnmgMf1wGN4gXMK8GPAsWOBP/qLxTTCC493mNliYCleSD6dJUCuv/jNGS1SA/wFL3gv8X8Wf/G3/x8wzMzm4E0vLfjZFz7HbGAd3tTev+GNAp/AXzF1MPCC/30W4U2hjQXeM7NUYCHwsnNu3xnWLiIiYWDOFZ6dIyIiIiIiIiWRRhBFREREREQEUEAUERERERERnwKiiIiIiIiIAAqIIiIiIiIi4lNAFBEREREREUABUURERERERHwKiCIiIiIiIgIoIIqIiIiIiIjv/wMk6t5sNwrSZAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x720 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Call the svd algorithm\n",
    "svd_algorithm(u_train, s_train, vt_train)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "`6.` Use the cell below to comment on the results you found in the previous question. Given the circumstances of your results, discuss what you might do to determine if the recommendations you make with any of the above recommendation systems are an improvement to how users currently find articles? "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* The Test accuracy decreases as the number of latent features increases for the testing dataset.\n",
    "* In this project, only 20 users had records of old interactions.\n",
    "* To solve the cold-start problem, we can deploy a rank based recommendation method or content based recommendation."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "<a id='conclusions'></a>\n",
    "### Extras\n",
    "Using your workbook, you could now save your recommendations for each user, develop a class to make new predictions and update your results, and make a flask app to deploy your results.  These tasks are beyond what is required for this project.  However, from what you learned in the lessons, you certainly capable of taking these tasks on to improve upon your work here!\n",
    "\n",
    "\n",
    "## Conclusion\n",
    "\n",
    "> Congratulations!  You have reached the end of the Recommendations with IBM project! \n",
    "\n",
    "> **Tip**: Once you are satisfied with your work here, check over your report to make sure that it is satisfies all the areas of the [rubric](https://review.udacity.com/#!/rubrics/2322/view). You should also probably remove all of the \"Tips\" like this one so that the presentation is as polished as possible.\n",
    "\n",
    "\n",
    "## Directions to Submit\n",
    "\n",
    "> Before you submit your project, you need to create a .html or .pdf version of this notebook in the workspace here. To do that, run the code cell below. If it worked correctly, you should get a return code of 0, and you should see the generated .html file in the workspace directory (click on the orange Jupyter icon in the upper left).\n",
    "\n",
    "> Alternatively, you can download this report as .html via the **File** > **Download as** submenu, and then manually upload it into the workspace directory by clicking on the orange Jupyter icon in the upper left, then using the Upload button.\n",
    "\n",
    "> Once you've done this, you can submit your project by clicking on the \"Submit Project\" button in the lower right here. This will create and submit a zip file with this .ipynb doc and the .html or .pdf version you created. Congratulations! "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 121,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 121,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from subprocess import call\n",
    "call(['python', '-m', 'nbconvert', 'Recommendations_with_IBM.ipynb'])"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}