2_Plagiarism_Feature_Engineering(1).ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plagiarism Detection, Feature Engineering\n",
    "\n",
    "In this project, you will be tasked with building a plagiarism detector that examines an answer text file and performs binary classification; labeling that file as either plagiarized or not, depending on how similar that text file is to a provided, source text. \n",
    "\n",
    "Your first task will be to create some features that can then be used to train a classification model. This task will be broken down into a few discrete steps:\n",
    "\n",
    "* Clean and pre-process the data.\n",
    "* Define features for comparing the similarity of an answer text and a source text, and extract similarity features.\n",
    "* Select \"good\" features, by analyzing the correlations between different features.\n",
    "* Create train/test `.csv` files that hold the relevant features and class labels for train/test data points.\n",
    "\n",
    "In the _next_ notebook, Notebook 3, you'll use the features and `.csv` files you create in _this_ notebook to train a binary classification model in a SageMaker notebook instance.\n",
    "\n",
    "You'll be defining a few different similarity features, as outlined in [this paper](https://s3.amazonaws.com/video.udacity-data.com/topher/2019/January/5c412841_developing-a-corpus-of-plagiarised-short-answers/developing-a-corpus-of-plagiarised-short-answers.pdf), which should help you build a robust plagiarism detector!\n",
    "\n",
    "To complete this notebook, you'll have to complete all given exercises and answer all the questions in this notebook.\n",
    "> All your tasks will be clearly labeled **EXERCISE** and questions as **QUESTION**.\n",
    "\n",
    "It will be up to you to decide on the features to include in your final training and test data.\n",
    "\n",
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in the Data\n",
    "\n",
    "The cell below will download the necessary, project data and extract the files into the folder `data/`.\n",
    "\n",
    "This data is a slightly modified version of a dataset created by Paul Clough (Information Studies) and Mark Stevenson (Computer Science), at the University of Sheffield. You can read all about the data collection and corpus, at [their university webpage](https://ir.shef.ac.uk/cloughie/resources/plagiarism_corpus.html). \n",
    "\n",
    "> **Citation for data**: Clough, P. and Stevenson, M. Developing A Corpus of Plagiarised Short Answers, Language Resources and Evaluation: Special Issue on Plagiarism and Authorship Analysis, In Press. [Download]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# NOTE:\n",
    "# you only need to run this cell if you have not yet downloaded the data\n",
    "# otherwise you may skip this cell or comment it out\n",
    "\n",
    "!wget https://s3.amazonaws.com/video.udacity-data.com/topher/2019/January/5c4147f9_data/data.zip\n",
    "!unzip data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# import libraries\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import os"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This plagiarism dataset is made of multiple text files; each of these files has characteristics that are is summarized in a `.csv` file named `file_information.csv`, which we can read in using `pandas`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "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>File</th>\n",
       "      <th>Task</th>\n",
       "      <th>Category</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>g0pA_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>non</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>g0pA_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>cut</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>g0pA_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>light</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>g0pA_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>heavy</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>g0pA_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>non</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             File Task Category\n",
       "0  g0pA_taska.txt    a      non\n",
       "1  g0pA_taskb.txt    b      cut\n",
       "2  g0pA_taskc.txt    c    light\n",
       "3  g0pA_taskd.txt    d    heavy\n",
       "4  g0pA_taske.txt    e      non"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "csv_file = 'data/file_information.csv'\n",
    "plagiarism_df = pd.read_csv(csv_file)\n",
    "\n",
    "# print out the first few rows of data info\n",
    "plagiarism_df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Types of Plagiarism\n",
    "\n",
    "Each text file is associated with one **Task** (task A-E) and one **Category** of plagiarism, which you can see in the above DataFrame.\n",
    "\n",
    "###  Tasks, A-E\n",
    "\n",
    "Each text file contains an answer to one short question; these questions are labeled as tasks A-E. For example, Task A asks the question: \"What is inheritance in object oriented programming?\"\n",
    "\n",
    "### Categories of plagiarism \n",
    "\n",
    "Each text file has an associated plagiarism label/category:\n",
    "\n",
    "**1. Plagiarized categories: `cut`, `light`, and `heavy`.**\n",
    "* These categories represent different levels of plagiarized answer texts. `cut` answers copy directly from a source text, `light` answers are based on the source text but include some light rephrasing, and `heavy` answers are based on the source text, but *heavily* rephrased (and will likely be the most challenging kind of plagiarism to detect).\n",
    "     \n",
    "**2. Non-plagiarized category: `non`.** \n",
    "* `non` indicates that an answer is not plagiarized; the Wikipedia source text is not used to create this answer.\n",
    "    \n",
    "**3. Special, source text category: `orig`.**\n",
    "* This is a specific category for the original, Wikipedia source text. We will use these files only for comparison purposes."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Pre-Process the Data\n",
    "\n",
    "In the next few cells, you'll be tasked with creating a new DataFrame of desired information about all of the files in the `data/` directory. This will prepare the data for feature extraction and for training a binary, plagiarism classifier."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### EXERCISE: Convert categorical to numerical data\n",
    "\n",
    "You'll notice that the `Category` column in the data, contains string or categorical values, and to prepare these for feature extraction, we'll want to convert these into numerical values. Additionally, our goal is to create a binary classifier and so we'll need a binary class label that indicates whether an answer text is plagiarized (1) or not (0). Complete the below function `numerical_dataframe` that reads in a `file_information.csv` file by name, and returns a *new* DataFrame with a numerical `Category` column and a new `Class` column that labels each answer as plagiarized or not. \n",
    "\n",
    "Your function should return a new DataFrame with the following properties:\n",
    "\n",
    "* 4 columns: `File`, `Task`, `Category`, `Class`. The `File` and `Task` columns can remain unchanged from the original `.csv` file.\n",
    "* Convert all `Category` labels to numerical labels according to the following rules (a higher value indicates a higher degree of plagiarism):\n",
    "    * 0 = `non`\n",
    "    * 1 = `heavy`\n",
    "    * 2 = `light`\n",
    "    * 3 = `cut`\n",
    "    * -1 = `orig`, this is a special value that indicates an original file.\n",
    "* For the new `Class` column\n",
    "    * Any answer text that is not plagiarized (`non`) should have the class label `0`. \n",
    "    * Any plagiarized answer texts should have the class label `1`. \n",
    "    * And any `orig` texts will have a special label `-1`. \n",
    "\n",
    "### Expected output\n",
    "\n",
    "After running your function, you should get a DataFrame with rows that looks like the following: \n",
    "```\n",
    "\n",
    "        File\t     Task  Category  Class\n",
    "0\tg0pA_taska.txt\ta\t  0   \t0\n",
    "1\tg0pA_taskb.txt\tb\t  3   \t1\n",
    "2\tg0pA_taskc.txt\tc\t  2   \t1\n",
    "3\tg0pA_taskd.txt\td\t  1   \t1\n",
    "4\tg0pA_taske.txt\te\t  0\t   0\n",
    "...\n",
    "...\n",
    "99   orig_taske.txt    e     -1      -1\n",
    "\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Read in a csv file and return a transformed dataframe\n",
    "def numerical_dataframe(csv_file='data/file_information.csv'):\n",
    "    '''Reads in a csv file which is assumed to have `File`, `Category` and `Task` columns.\n",
    "       This function does two things: \n",
    "       1) converts `Category` column values to numerical values \n",
    "       2) Adds a new, numerical `Class` label column.\n",
    "       The `Class` column will label plagiarized answers as 1 and non-plagiarized as 0.\n",
    "       Source texts have a special label, -1.\n",
    "       :param csv_file: The directory for the file_information.csv file\n",
    "       :return: A dataframe with numerical categories and a new `Class` label column'''\n",
    "    \n",
    "    df=pd.read_csv(csv_file)\n",
    "    df['Category']=[0 if df.loc[i,'Category']=='non'\n",
    "                 else 1 if df.loc[i,'Category']=='heavy'\n",
    "                 else 2 if df.loc[i,'Category']=='light'\n",
    "                 else 3 if df.loc[i,'Category']=='cut'\n",
    "                 else -1\n",
    "                 for i in df.index]\n",
    "    \n",
    "    df['Class']=[-1 if df.loc[i,'Category']==-1\n",
    "                 else 0 if df.loc[i,'Category']==0\n",
    "                 else 1\n",
    "                 for i in df.index]\n",
    "    \n",
    "    \n",
    "    return df\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Test cells\n",
    "\n",
    "Below are a couple of test cells. The first is an informal test where you can check that your code is working as expected by calling your function and printing out the returned result.\n",
    "\n",
    "The **second** cell below is a more rigorous test cell. The goal of a cell like this is to ensure that your code is working as expected, and to form any variables that might be used in _later_ tests/code, in this case, the data frame, `transformed_df`.\n",
    "\n",
    "> The cells in this notebook should be run in chronological order (the order they appear in the notebook). This is especially important for test cells.\n",
    "\n",
    "Often, later cells rely on the functions, imports, or variables defined in earlier cells. For example, some tests rely on previous tests to work.\n",
    "\n",
    "These tests do not test all cases, but they are a great way to check that you are on the right track!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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>File</th>\n",
       "      <th>Task</th>\n",
       "      <th>Category</th>\n",
       "      <th>Class</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>g0pA_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>g0pA_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>g0pA_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>g0pA_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>g0pA_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>g0pB_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>g0pB_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>g0pB_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>g0pB_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>g0pB_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             File Task  Category  Class\n",
       "0  g0pA_taska.txt    a         0      0\n",
       "1  g0pA_taskb.txt    b         3      1\n",
       "2  g0pA_taskc.txt    c         2      1\n",
       "3  g0pA_taskd.txt    d         1      1\n",
       "4  g0pA_taske.txt    e         0      0\n",
       "5  g0pB_taska.txt    a         0      0\n",
       "6  g0pB_taskb.txt    b         0      0\n",
       "7  g0pB_taskc.txt    c         3      1\n",
       "8  g0pB_taskd.txt    d         2      1\n",
       "9  g0pB_taske.txt    e         1      1"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# informal testing, print out the results of a called function\n",
    "# create new `transformed_df`\n",
    "transformed_df = numerical_dataframe(csv_file ='data/file_information.csv')\n",
    "\n",
    "# check work\n",
    "# check that all categories of plagiarism have a class label = 1\n",
    "transformed_df.head(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "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>File</th>\n",
       "      <th>Task</th>\n",
       "      <th>Category</th>\n",
       "      <th>Class</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>90</th>\n",
       "      <td>g4pE_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>91</th>\n",
       "      <td>g4pE_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92</th>\n",
       "      <td>g4pE_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>93</th>\n",
       "      <td>g4pE_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>94</th>\n",
       "      <td>g4pE_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95</th>\n",
       "      <td>orig_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>96</th>\n",
       "      <td>orig_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>97</th>\n",
       "      <td>orig_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98</th>\n",
       "      <td>orig_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99</th>\n",
       "      <td>orig_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              File Task  Category  Class\n",
       "90  g4pE_taska.txt    a         1      1\n",
       "91  g4pE_taskb.txt    b         2      1\n",
       "92  g4pE_taskc.txt    c         3      1\n",
       "93  g4pE_taskd.txt    d         0      0\n",
       "94  g4pE_taske.txt    e         0      0\n",
       "95  orig_taska.txt    a        -1     -1\n",
       "96  orig_taskb.txt    b        -1     -1\n",
       "97  orig_taskc.txt    c        -1     -1\n",
       "98  orig_taskd.txt    d        -1     -1\n",
       "99  orig_taske.txt    e        -1     -1"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "transformed_df.tail(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Tests Passed!\n",
      "\n",
      "Example data: \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>File</th>\n",
       "      <th>Task</th>\n",
       "      <th>Category</th>\n",
       "      <th>Class</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>g0pA_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>g0pA_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>g0pA_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>g0pA_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>g0pA_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             File Task  Category  Class\n",
       "0  g0pA_taska.txt    a         0      0\n",
       "1  g0pA_taskb.txt    b         3      1\n",
       "2  g0pA_taskc.txt    c         2      1\n",
       "3  g0pA_taskd.txt    d         1      1\n",
       "4  g0pA_taske.txt    e         0      0"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# test cell that creates `transformed_df`, if tests are passed\n",
    "\n",
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "\n",
    "# importing tests\n",
    "import problem_unittests as tests\n",
    "\n",
    "# test numerical_dataframe function\n",
    "tests.test_numerical_df(numerical_dataframe)\n",
    "\n",
    "# if above test is passed, create NEW `transformed_df`\n",
    "transformed_df = numerical_dataframe(csv_file ='data/file_information.csv')\n",
    "\n",
    "# check work\n",
    "print('\\nExample data: ')\n",
    "transformed_df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Text Processing & Splitting Data\n",
    "\n",
    "Recall that the goal of this project is to build a plagiarism classifier. At it's heart, this task is a comparison text; one that looks at a given answer and a source text, compares them and predicts whether an answer has plagiarized from the source. To effectively do this comparison, and train a classifier we'll need to do a few more things: pre-process all of our text data and prepare the text files (in this case, the 95 answer files and 5 original source files) to be easily compared, and split our data into a `train` and `test` set that can be used to train a classifier and evaluate it, respectively. \n",
    "\n",
    "To this end, you've been provided code that adds  additional information to your `transformed_df` from above. The next two cells need not be changed; they add two additional columns to the `transformed_df`:\n",
    "\n",
    "1. A `Text` column; this holds all the lowercase text for a `File`, with extraneous punctuation removed.\n",
    "2. A `Datatype` column; this is a string value `train`, `test`, or `orig` that labels a data point as part of our train or test set\n",
    "\n",
    "The details of how these additional columns are created can be found in the `helpers.py` file in the project directory. You're encouraged to read through that file to see exactly how text is processed and how data is split.\n",
    "\n",
    "Run the cells below to get a `complete_df` that has all the information you need to proceed with plagiarism detection and feature engineering."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "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>File</th>\n",
       "      <th>Task</th>\n",
       "      <th>Category</th>\n",
       "      <th>Class</th>\n",
       "      <th>Text</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>g0pA_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>inheritance is a basic concept of object orien...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>g0pA_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>pagerank is a link analysis algorithm used by ...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>g0pA_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>the vector space model also called term vector...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>g0pA_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>bayes theorem was names after rev thomas bayes...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>g0pA_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>dynamic programming is an algorithm design tec...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             File Task  Category  Class  \\\n",
       "0  g0pA_taska.txt    a         0      0   \n",
       "1  g0pA_taskb.txt    b         3      1   \n",
       "2  g0pA_taskc.txt    c         2      1   \n",
       "3  g0pA_taskd.txt    d         1      1   \n",
       "4  g0pA_taske.txt    e         0      0   \n",
       "\n",
       "                                                Text  \n",
       "0  inheritance is a basic concept of object orien...  \n",
       "1  pagerank is a link analysis algorithm used by ...  \n",
       "2  the vector space model also called term vector...  \n",
       "3  bayes theorem was names after rev thomas bayes...  \n",
       "4  dynamic programming is an algorithm design tec...  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "import helpers \n",
    "\n",
    "# create a text column \n",
    "text_df = helpers.create_text_column(transformed_df)\n",
    "text_df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Sample processed text:\n",
      "\n",
      " inheritance is a basic concept of object oriented programming where the basic idea is to create new classes that add extra detail to existing classes this is done by allowing the new classes to reuse the methods and variables of the existing classes and new methods and classes are added to specialise the new class inheritance models the is kind of relationship between entities or objects  for example postgraduates and undergraduates are both kinds of student this kind of relationship can be visualised as a tree structure where student would be the more general root node and both postgraduate and undergraduate would be more specialised extensions of the student node or the child nodes  in this relationship student would be known as the superclass or parent class whereas  postgraduate would be known as the subclass or child class because the postgraduate class extends the student class  inheritance can occur on several layers where if visualised would display a larger tree structure for example we could further extend the postgraduate node by adding two extra extended classes to it called  msc student and phd student as both these types of student are kinds of postgraduate student this would mean that both the msc student and phd student classes would inherit methods and variables from both the postgraduate and student classes  \n"
     ]
    }
   ],
   "source": [
    "# after running the cell above\n",
    "# check out the processed text for a single file, by row index\n",
    "row_idx = 0 # feel free to change this index\n",
    "\n",
    "sample_text = text_df.iloc[0]['Text']\n",
    "\n",
    "print('Sample processed text:\\n\\n', sample_text)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Split data into training and test sets\n",
    "\n",
    "The next cell will add a `Datatype` column to a given DataFrame to indicate if the record is: \n",
    "* `train` - Training data, for model training.\n",
    "* `test` - Testing data, for model evaluation.\n",
    "* `orig` - The task's original answer from wikipedia.\n",
    "\n",
    "### Stratified sampling\n",
    "\n",
    "The given code uses a helper function which you can view in the `helpers.py` file in the main project directory. This implements [stratified random sampling](https://en.wikipedia.org/wiki/Stratified_sampling) to randomly split data by task & plagiarism amount. Stratified sampling ensures that we get training and test data that is fairly evenly distributed across task & plagiarism combinations. Approximately 26% of the data is held out for testing and 74% of the data is used for training.\n",
    "\n",
    "The function **train_test_dataframe** takes in a DataFrame that it assumes has `Task` and `Category` columns, and, returns a modified frame that indicates which `Datatype` (train, test, or orig) a file falls into. This sampling will change slightly based on a passed in *random_seed*. Due to a small sample size, this stratified random sampling will provide more stable results for a binary plagiarism classifier. Stability here is smaller *variance* in the accuracy of classifier, given a random seed."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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>File</th>\n",
       "      <th>Task</th>\n",
       "      <th>Category</th>\n",
       "      <th>Class</th>\n",
       "      <th>Text</th>\n",
       "      <th>Datatype</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>g0pA_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>inheritance is a basic concept of object orien...</td>\n",
       "      <td>train</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>g0pA_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>pagerank is a link analysis algorithm used by ...</td>\n",
       "      <td>test</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>g0pA_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>the vector space model also called term vector...</td>\n",
       "      <td>train</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>g0pA_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>bayes theorem was names after rev thomas bayes...</td>\n",
       "      <td>train</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>g0pA_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>dynamic programming is an algorithm design tec...</td>\n",
       "      <td>train</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>g0pB_taska.txt</td>\n",
       "      <td>a</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>inheritance is a basic concept in object orien...</td>\n",
       "      <td>train</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>g0pB_taskb.txt</td>\n",
       "      <td>b</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>pagerank pr refers to both the concept and the...</td>\n",
       "      <td>train</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>g0pB_taskc.txt</td>\n",
       "      <td>c</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>vector space model is an algebraic model for r...</td>\n",
       "      <td>test</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>g0pB_taskd.txt</td>\n",
       "      <td>d</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>bayes theorem relates the conditional and marg...</td>\n",
       "      <td>train</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>g0pB_taske.txt</td>\n",
       "      <td>e</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>dynamic programming is a method for solving ma...</td>\n",
       "      <td>test</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             File Task  Category  Class  \\\n",
       "0  g0pA_taska.txt    a         0      0   \n",
       "1  g0pA_taskb.txt    b         3      1   \n",
       "2  g0pA_taskc.txt    c         2      1   \n",
       "3  g0pA_taskd.txt    d         1      1   \n",
       "4  g0pA_taske.txt    e         0      0   \n",
       "5  g0pB_taska.txt    a         0      0   \n",
       "6  g0pB_taskb.txt    b         0      0   \n",
       "7  g0pB_taskc.txt    c         3      1   \n",
       "8  g0pB_taskd.txt    d         2      1   \n",
       "9  g0pB_taske.txt    e         1      1   \n",
       "\n",
       "                                                Text Datatype  \n",
       "0  inheritance is a basic concept of object orien...    train  \n",
       "1  pagerank is a link analysis algorithm used by ...     test  \n",
       "2  the vector space model also called term vector...    train  \n",
       "3  bayes theorem was names after rev thomas bayes...    train  \n",
       "4  dynamic programming is an algorithm design tec...    train  \n",
       "5  inheritance is a basic concept in object orien...    train  \n",
       "6  pagerank pr refers to both the concept and the...    train  \n",
       "7  vector space model is an algebraic model for r...     test  \n",
       "8  bayes theorem relates the conditional and marg...    train  \n",
       "9  dynamic programming is a method for solving ma...     test  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "random_seed = 1 # can change; set for reproducibility\n",
    "\n",
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "import helpers\n",
    "\n",
    "# create new df with Datatype (train, test, orig) column\n",
    "# pass in `text_df` from above to create a complete dataframe, with all the information you need\n",
    "complete_df = helpers.train_test_dataframe(text_df, random_seed=random_seed)\n",
    "\n",
    "# check results\n",
    "complete_df.head(10)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Determining Plagiarism\n",
    "\n",
    "Now that you've prepared this data and created a `complete_df` of information, including the text and class associated with each file, you can move on to the task of extracting similarity features that will be useful for plagiarism classification. \n",
    "\n",
    "> Note: The following code exercises, assume that the `complete_df` as it exists now, will **not** have its existing columns modified. \n",
    "\n",
    "The `complete_df` should always include the columns: `['File', 'Task', 'Category', 'Class', 'Text', 'Datatype']`. You can add additional columns, and you can create any new DataFrames you need by copying the parts of the `complete_df` as long as you do not modify the existing values, directly.\n",
    "\n",
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "# Similarity Features \n",
    "\n",
    "One of the ways we might go about detecting plagiarism, is by computing **similarity features** that measure how similar a given answer text is as compared to the original wikipedia source text (for a specific task, a-e). The similarity features you will use are informed by [this paper on plagiarism detection](https://s3.amazonaws.com/video.udacity-data.com/topher/2019/January/5c412841_developing-a-corpus-of-plagiarised-short-answers/developing-a-corpus-of-plagiarised-short-answers.pdf). \n",
    "> In this paper, researchers created features called **containment** and **longest common subsequence**. \n",
    "\n",
    "Using these features as input, you will train a model to distinguish between plagiarized and not-plagiarized text files.\n",
    "\n",
    "## Feature Engineering\n",
    "\n",
    "Let's talk a bit more about the features we want to include in a plagiarism detection model and how to calculate such features. In the following explanations, I'll refer to a submitted text file as a **Student Answer Text (A)** and the original, wikipedia source file (that we want to compare that answer to) as the **Wikipedia Source Text (S)**.\n",
    "\n",
    "### Containment\n",
    "\n",
    "Your first task will be to create **containment features**. To understand containment, let's first revisit a definition of [n-grams](https://en.wikipedia.org/wiki/N-gram). An *n-gram* is a sequential word grouping. For example, in a line like \"bayes rule gives us a way to combine prior knowledge with new information,\" a 1-gram is just one word, like \"bayes.\" A 2-gram might be \"bayes rule\" and a 3-gram might be \"combine prior knowledge.\"\n",
    "\n",
    "> Containment is defined as the **intersection** of the n-gram word count of the Wikipedia Source Text (S) with the n-gram word count of the Student  Answer Text (S) *divided* by the n-gram word count of the Student Answer Text.\n",
    "\n",
    "$$ \\frac{\\sum{count(\\text{ngram}_{A}) \\cap count(\\text{ngram}_{S})}}{\\sum{count(\\text{ngram}_{A})}} $$\n",
    "\n",
    "If the two texts have no n-grams in common, the containment will be 0, but if _all_ their n-grams intersect then the containment will be 1. Intuitively, you can see how having longer n-gram's in common, might be an indication of cut-and-paste plagiarism. In this project, it will be up to you to decide on the appropriate `n` or several `n`'s to use in your final model.\n",
    "\n",
    "### EXERCISE: Create containment features\n",
    "\n",
    "Given the `complete_df` that you've created, you should have all the information you need to compare any Student  Answer Text (A) with its appropriate Wikipedia Source Text (S). An answer for task A should be compared to the source text for task A, just as answers to tasks B, C, D, and E should be compared to the corresponding original source text.\n",
    "\n",
    "In this exercise, you'll complete the function, `calculate_containment` which calculates containment based upon the following parameters:\n",
    "* A given DataFrame, `df` (which is assumed to be the `complete_df` from above)\n",
    "* An `answer_filename`, such as 'g0pB_taskd.txt' \n",
    "* An n-gram length, `n`\n",
    "\n",
    "### Containment calculation\n",
    "\n",
    "The general steps to complete this function are as follows:\n",
    "1. From *all* of the text files in a given `df`, create an array of n-gram counts; it is suggested that you use a [CountVectorizer](https://scikit-learn.org/stable/modules/generated/sklearn.feature_extraction.text.CountVectorizer.html) for this purpose.\n",
    "2. Get the processed answer and source texts for the given `answer_filename`.\n",
    "3. Calculate the containment between an answer and source text according to the following equation.\n",
    "\n",
    "    >$$ \\frac{\\sum{count(\\text{ngram}_{A}) \\cap count(\\text{ngram}_{S})}}{\\sum{count(\\text{ngram}_{A})}} $$\n",
    "    \n",
    "4. Return that containment value.\n",
    "\n",
    "You are encouraged to write any helper functions that you need to complete the function below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Calculate the ngram containment for one answer file/source file pair in a df\n",
    "from sklearn.feature_extraction.text import CountVectorizer\n",
    "def calculate_containment(df, n, answer_filename):\n",
    "    '''Calculates the containment between a given answer text and its associated source text.\n",
    "       This function creates a count of ngrams (of a size, n) for each text file in our data.\n",
    "       Then calculates the containment by finding the ngram count for a given answer text, \n",
    "       and its associated source text, and calculating the normalized intersection of those counts.\n",
    "       :param df: A dataframe with columns,\n",
    "           'File', 'Task', 'Category', 'Class', 'Text', and 'Datatype'\n",
    "       :param n: An integer that defines the ngram size\n",
    "       :param answer_filename: A filename for an answer text in the df, ex. 'g0pB_taskd.txt'\n",
    "       :return: A single containment value that represents the similarity\n",
    "           between an answer text and its source text.\n",
    "    '''\n",
    "    counts = CountVectorizer(analyzer='word', ngram_range=(n,n))\n",
    "    a_text = df.loc[df[\"File\"] == answer_filename][\"Text\"].values[0]\n",
    "    a_task = df.loc[df[\"File\"] == answer_filename][\"Task\"].values[0]\n",
    "    s_text = df.loc[(df[\"Task\"] == a_task) & (df[\"Datatype\"] == 'orig')][\"Text\"].values[0]\n",
    "    ngrams = counts.fit_transform([a_text, s_text])\n",
    "    ngram_array = ngrams.toarray()\n",
    "    result=np.amin(ngram_array,axis=0).sum()/ngram_array[0].sum()\n",
    "    return result"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Test cells\n",
    "\n",
    "After you've implemented the containment function, you can test out its behavior. \n",
    "\n",
    "The cell below iterates through the first few files, and calculates the original category _and_ containment values for a specified n and file.\n",
    "\n",
    ">If you've implemented this correctly, you should see that the non-plagiarized have low or close to 0 containment values and that plagiarized examples have higher containment values, closer to 1.\n",
    "\n",
    "Note what happens when you change the value of n. I recommend applying your code to multiple files and comparing the resultant containment values. You should see that the highest containment values correspond to files with the highest category (`cut`) of plagiarism level."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Original category values: \n",
      " [0, 3, 2, 1, 0]\n",
      "\n",
      "3-gram containment values: \n",
      " [0.009345794392523364, 0.9641025641025641, 0.6136363636363636, 0.15675675675675677, 0.031746031746031744]\n"
     ]
    }
   ],
   "source": [
    "# select a value for n\n",
    "n = 3\n",
    "\n",
    "# indices for first few files\n",
    "test_indices = range(5)\n",
    "\n",
    "# iterate through files and calculate containment\n",
    "category_vals = []\n",
    "containment_vals = []\n",
    "for i in test_indices:\n",
    "    # get level of plagiarism for a given file index\n",
    "    category_vals.append(complete_df.loc[i, 'Category'])\n",
    "    # calculate containment for given file and n\n",
    "    filename = complete_df.loc[i, 'File']\n",
    "    c = calculate_containment(complete_df, n, filename)\n",
    "    containment_vals.append(c)\n",
    "\n",
    "# print out result, does it make sense?\n",
    "print('Original category values: \\n', category_vals)\n",
    "print()\n",
    "print(str(n)+'-gram containment values: \\n', containment_vals)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Tests Passed!\n"
     ]
    }
   ],
   "source": [
    "# run this test cell\n",
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "# test containment calculation\n",
    "# params: complete_df from before, and containment function\n",
    "tests.test_containment(complete_df, calculate_containment)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### QUESTION 1: Why can we calculate containment features across *all* data (training & test), prior to splitting the DataFrame for modeling? That is, what about the containment calculation means that the test and training data do not influence each other?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Answer:**\n",
    "<ol>\n",
    "\n",
    "   <li>Source text are non changing accross train or test data\n",
    "   </li>\n",
    "    \n",
    "   <li> \n",
    "    The test and the train splits are independent from each other and don't affect their containment result.\n",
    "    </li>\n",
    "    \n",
    "</ol>\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Longest Common Subsequence\n",
    "\n",
    "Containment a good way to find overlap in word usage between two documents; it may help identify cases of cut-and-paste as well as paraphrased levels of plagiarism. Since plagiarism is a fairly complex task with varying levels, it's often useful to include other measures of similarity. The paper also discusses a feature called **longest common subsequence**.\n",
    "\n",
    "> The longest common subsequence is the longest string of words (or letters) that are *the same* between the Wikipedia Source Text (S) and the Student Answer Text (A). This value is also normalized by dividing by the total number of words (or letters) in the  Student Answer Text. \n",
    "\n",
    "In this exercise, we'll ask you to calculate the longest common subsequence of words between two texts.\n",
    "\n",
    "### EXERCISE: Calculate the longest common subsequence\n",
    "\n",
    "Complete the function `lcs_norm_word`; this should calculate the *longest common subsequence* of words between a Student Answer Text and corresponding Wikipedia Source Text. \n",
    "\n",
    "It may be helpful to think of this in a concrete example. A Longest Common Subsequence (LCS) problem may look as follows:\n",
    "* Given two texts: text A (answer text) of length n, and string S (original source text) of length m. Our goal is to produce their longest common subsequence of words: the longest sequence of words that appear left-to-right in both texts (though the words don't have to be in continuous order).\n",
    "* Consider:\n",
    "    * A = \"i think pagerank is a link analysis algorithm used by google that uses a system of weights attached to each element of a hyperlinked set of documents\"\n",
    "    * S = \"pagerank is a link analysis algorithm used by the google internet search engine that assigns a numerical weighting to each element of a hyperlinked set of documents\"\n",
    "\n",
    "* In this case, we can see that the start of each sentence of fairly similar, having overlap in the sequence of words, \"pagerank is a link analysis algorithm used by\" before diverging slightly. Then we **continue moving left -to-right along both texts** until we see the next common sequence; in this case it is only one word, \"google\". Next we find \"that\" and \"a\" and finally the same ending \"to each element of a hyperlinked set of documents\".\n",
    "* Below, is a clear visual of how these sequences were found, sequentially, in each text.\n",
    "\n",
    "<img src='notebook_ims/common_subseq_words.png' width=40% />\n",
    "\n",
    "* Now, those words appear in left-to-right order in each document, sequentially, and even though there are some words in between, we count this as the longest common subsequence between the two texts. \n",
    "* If I count up each word that I found in common I get the value 20. **So, LCS has length 20**. \n",
    "* Next, to normalize this value, divide by the total length of the student answer; in this example that length is only 27. **So, the function `lcs_norm_word` should return the value `20/27` or about `0.7408`.**\n",
    "\n",
    "In this way, LCS is a great indicator of cut-and-paste plagiarism or if someone has referenced the same source text multiple times in an answer."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### LCS, dynamic programming\n",
    "\n",
    "If you read through the scenario above, you can see that this algorithm depends on looking at two texts and comparing them word by word. You can solve this problem in multiple ways. First, it may be useful to `.split()` each text into lists of comma separated words to compare. Then, you can iterate through each word in the texts and compare them, adding to your value for LCS as you go. \n",
    "\n",
    "The method I recommend for implementing an efficient LCS algorithm is: using a matrix and dynamic programming. **Dynamic programming** is all about breaking a larger problem into a smaller set of subproblems, and building up a complete result without having to repeat any subproblems. \n",
    "\n",
    "This approach assumes that you can split up a large LCS task into a combination of smaller LCS tasks. Let's look at a simple example that compares letters:\n",
    "\n",
    "* A = \"ABCD\"\n",
    "* S = \"BD\"\n",
    "\n",
    "We can see right away that the longest subsequence of _letters_ here is 2 (B and D are in sequence in both strings). And we can calculate this by looking at relationships between each letter in the two strings, A and S.\n",
    "\n",
    "Here, I have a matrix with the letters of A on top and the letters of S on the left side:\n",
    "\n",
    "<img src='notebook_ims/matrix_1.png' width=40% />\n",
    "\n",
    "This starts out as a matrix that has as many columns and rows as letters in the strings S and O **+1** additional row and column, filled with zeros on the top and left sides. So, in this case, instead of a 2x4 matrix it is a 3x5.\n",
    "\n",
    "Now, we can fill this matrix up by breaking it into smaller LCS problems. For example, let's first look at the shortest substrings: the starting letter of A and S. We'll first ask, what is the Longest Common Subsequence between these two letters \"A\" and \"B\"? \n",
    "\n",
    "**Here, the answer is zero and we fill in the corresponding grid cell with that value.**\n",
    "\n",
    "<img src='notebook_ims/matrix_2.png' width=30% />\n",
    "\n",
    "Then, we ask the next question, what is the LCS between \"AB\" and \"B\"?\n",
    "\n",
    "**Here, we have a match, and can fill in the appropriate value 1**.\n",
    "\n",
    "<img src='notebook_ims/matrix_3_match.png' width=25% />\n",
    "\n",
    "If we continue, we get to a final matrix that looks as follows, with a **2** in the bottom right corner.\n",
    "\n",
    "<img src='notebook_ims/matrix_6_complete.png' width=25% />\n",
    "\n",
    "The final LCS will be that value **2** *normalized* by the number of n-grams in A. So, our normalized value is 2/4 = **0.5**.\n",
    "\n",
    "### The matrix rules\n",
    "\n",
    "One thing to notice here is that, you can efficiently fill up this matrix one cell at a time. Each grid cell only depends on the values in the grid cells that are directly on top and to the left of it, or on the diagonal/top-left. The rules are as follows:\n",
    "* Start with a matrix that has one extra row and column of zeros.\n",
    "* As you traverse your string:\n",
    "    * If there is a match, fill that grid cell with the value to the top-left of that cell *plus* one. So, in our case, when we found a matching B-B, we added +1 to the value in the top-left of the matching cell, 0.\n",
    "    * If there is not a match, take the *maximum* value from either directly to the left or the top cell, and carry that value over to the non-match cell.\n",
    "\n",
    "<img src='notebook_ims/matrix_rules.png' width=50% />\n",
    "\n",
    "After completely filling the matrix, **the bottom-right cell will hold the non-normalized LCS value**.\n",
    "\n",
    "This matrix treatment can be applied to a set of words instead of letters. Your function should apply this to the words in two texts and return the normalized LCS value."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Compute the normalized LCS given an answer text and a source text\n",
    "def lcs_norm_word(answer_text, source_text):\n",
    "    \n",
    "    '''Computes the longest common subsequence of words in two texts; returns a normalized value.\n",
    "       :param answer_text: The pre-processed text for an answer text\n",
    "       :param source_text: The pre-processed text for an answer's associated source text\n",
    "       :return: A normalized LCS value'''\n",
    "    splitted_answer = answer_text.split()\n",
    "    splitted_source = source_text.split()\n",
    "    n_a=len(splitted_answer)+1\n",
    "    n_s=len(splitted_source)+1\n",
    "    matrix =  [[0]*(n_a) for _ in range(n_s)]\n",
    "    for i in range(1,len(matrix)):\n",
    "        for j in range(1,len(matrix[0])):\n",
    "            if splitted_answer[j-1] == splitted_source[i-1]:\n",
    "                matrix[i][j] = 1 + matrix[i-1][j-1]\n",
    "            else:\n",
    "                matrix[i][j] = max(matrix[i-1][j], matrix[i][j-1])\n",
    "    return matrix[-1][-1] /( n_a-1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Test cells\n",
    "\n",
    "Let's start by testing out your code on the example given in the initial description.\n",
    "\n",
    "In the below cell, we have specified strings A (answer text) and S (original source text). We know that these texts have 20 words in common and the submitted answer is 27 words long, so the normalized, longest common subsequence should be 20/27.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "LCS =  0.7407407407407407\n",
      "Test passed!\n"
     ]
    }
   ],
   "source": [
    "# Run the test scenario from above\n",
    "# does your function return the expected value?\n",
    "\n",
    "A = \"i think pagerank is a link analysis algorithm used by google that uses a system of weights attached to each element of a hyperlinked set of documents\"\n",
    "S = \"pagerank is a link analysis algorithm used by the google internet search engine that assigns a numerical weighting to each element of a hyperlinked set of documents\"\n",
    "\n",
    "# calculate LCS\n",
    "lcs = lcs_norm_word(A, S)\n",
    "print('LCS = ', lcs)\n",
    "\n",
    "\n",
    "# expected value test\n",
    "assert lcs==20/27., \"Incorrect LCS value, expected about 0.7408, got \"+str(lcs)\n",
    "\n",
    "print('Test passed!')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This next cell runs a more rigorous test."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Tests Passed!\n"
     ]
    }
   ],
   "source": [
    "# run test cell\n",
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "# test lcs implementation\n",
    "# params: complete_df from before, and lcs_norm_word function\n",
    "tests.test_lcs(complete_df, lcs_norm_word)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Finally, take a look at a few resultant values for `lcs_norm_word`. Just like before, you should see that higher values correspond to higher levels of plagiarism."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Original category values: \n",
      " [0, 3, 2, 1, 0]\n",
      "\n",
      "Normalized LCS values: \n",
      " [0.1917808219178082, 0.8207547169811321, 0.8464912280701754, 0.3160621761658031, 0.24257425742574257]\n"
     ]
    }
   ],
   "source": [
    "# test on your own\n",
    "test_indices = range(5) # look at first few files\n",
    "\n",
    "category_vals = []\n",
    "lcs_norm_vals = []\n",
    "# iterate through first few docs and calculate LCS\n",
    "for i in test_indices:\n",
    "    category_vals.append(complete_df.loc[i, 'Category'])\n",
    "    # get texts to compare\n",
    "    answer_text = complete_df.loc[i, 'Text'] \n",
    "    task = complete_df.loc[i, 'Task']\n",
    "    # we know that source texts have Class = -1\n",
    "    orig_rows = complete_df[(complete_df['Class'] == -1)]\n",
    "    orig_row = orig_rows[(orig_rows['Task'] == task)]\n",
    "    source_text = orig_row['Text'].values[0]\n",
    "    \n",
    "    # calculate lcs\n",
    "    lcs_val = lcs_norm_word(answer_text, source_text)\n",
    "    lcs_norm_vals.append(lcs_val)\n",
    "\n",
    "# print out result, does it make sense?\n",
    "print('Original category values: \\n', category_vals)\n",
    "print()\n",
    "print('Normalized LCS values: \\n', lcs_norm_vals)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "# Create All Features\n",
    "\n",
    "Now that you've completed the feature calculation functions, it's time to actually create multiple features and decide on which ones to use in your final model! In the below cells, you're provided two helper functions to help you create multiple features and store those in a DataFrame, `features_df`.\n",
    "\n",
    "### Creating multiple containment features\n",
    "\n",
    "Your completed `calculate_containment` function will be called in the next cell, which defines the helper function `create_containment_features`. \n",
    "\n",
    "> This function returns a list of containment features, calculated for a given `n` and for *all* files in a df (assumed to the the `complete_df`).\n",
    "\n",
    "For our original files, the containment value is set to a special value, -1.\n",
    "\n",
    "This function gives you the ability to easily create several containment features, of different n-gram lengths, for each of our text files."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "# Function returns a list of containment features, calculated for a given n \n",
    "# Should return a list of length 100 for all files in a complete_df\n",
    "def create_containment_features(df, n, column_name=None):\n",
    "    \n",
    "    containment_values = []\n",
    "    \n",
    "    if(column_name==None):\n",
    "        column_name = 'c_'+str(n) # c_1, c_2, .. c_n\n",
    "    \n",
    "    # iterates through dataframe rows\n",
    "    for i in df.index:\n",
    "        file = df.loc[i, 'File']\n",
    "        # Computes features using calculate_containment function\n",
    "        if df.loc[i,'Category'] > -1:\n",
    "            c = calculate_containment(df, n, file)\n",
    "            containment_values.append(c)\n",
    "        # Sets value to -1 for original tasks \n",
    "        else:\n",
    "            containment_values.append(-1)\n",
    "    \n",
    "    print(str(n)+'-gram containment features created!')\n",
    "    return containment_values\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Creating LCS features\n",
    "\n",
    "Below, your complete `lcs_norm_word` function is used to create a list of LCS features for all the answer files in a given DataFrame (again, this assumes you are passing in the `complete_df`. It assigns a special value for our original, source files, -1.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "# Function creates lcs feature and add it to the dataframe\n",
    "def create_lcs_features(df, column_name='lcs_word'):\n",
    "    \n",
    "    lcs_values = []\n",
    "    \n",
    "    # iterate through files in dataframe\n",
    "    for i in df.index:\n",
    "        # Computes LCS_norm words feature using function above for answer tasks\n",
    "        if df.loc[i,'Category'] > -1:\n",
    "            # get texts to compare\n",
    "            answer_text = df.loc[i, 'Text'] \n",
    "            task = df.loc[i, 'Task']\n",
    "            # we know that source texts have Class = -1\n",
    "            orig_rows = df[(df['Class'] == -1)]\n",
    "            orig_row = orig_rows[(orig_rows['Task'] == task)]\n",
    "            source_text = orig_row['Text'].values[0]\n",
    "\n",
    "            # calculate lcs\n",
    "            lcs = lcs_norm_word(answer_text, source_text)\n",
    "            lcs_values.append(lcs)\n",
    "        # Sets to -1 for original tasks \n",
    "        else:\n",
    "            lcs_values.append(-1)\n",
    "\n",
    "    print('LCS features created!')\n",
    "    return lcs_values\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## EXERCISE: Create a features DataFrame by selecting an `ngram_range`\n",
    "\n",
    "The paper suggests calculating the following features: containment *1-gram to 5-gram* and *longest common subsequence*. \n",
    "> In this exercise, you can choose to create even more features, for example from *1-gram to 7-gram* containment features and *longest common subsequence*. \n",
    "\n",
    "You'll want to create at least 6 features to choose from as you think about which to give to your final, classification model. Defining and comparing at least 6 different features allows you to discard any features that seem redundant, and choose to use the best features for your final model!\n",
    "\n",
    "In the below cell **define an n-gram range**; these will be the n's you use to create n-gram containment features. The rest of the feature creation code is provided."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1-gram containment features created!\n",
      "2-gram containment features created!\n",
      "3-gram containment features created!\n",
      "4-gram containment features created!\n",
      "5-gram containment features created!\n",
      "6-gram containment features created!\n",
      "LCS features created!\n",
      "\n",
      "Features:  ['c_1', 'c_2', 'c_3', 'c_4', 'c_5', 'c_6', 'lcs_word']\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# Define an ngram range\n",
    "ngram_range = range(1,7)\n",
    "\n",
    "\n",
    "# The following code may take a minute to run, depending on your ngram_range\n",
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "features_list = []\n",
    "\n",
    "# Create features in a features_df\n",
    "all_features = np.zeros((len(ngram_range)+1, len(complete_df)))\n",
    "\n",
    "# Calculate features for containment for ngrams in range\n",
    "i=0\n",
    "for n in ngram_range:\n",
    "    column_name = 'c_'+str(n)\n",
    "    features_list.append(column_name)\n",
    "    # create containment features\n",
    "    all_features[i]=np.squeeze(create_containment_features(complete_df, n))\n",
    "    i+=1\n",
    "\n",
    "# Calculate features for LCS_Norm Words \n",
    "features_list.append('lcs_word')\n",
    "all_features[i]= np.squeeze(create_lcs_features(complete_df))\n",
    "\n",
    "# create a features dataframe\n",
    "features_df = pd.DataFrame(np.transpose(all_features), columns=features_list)\n",
    "\n",
    "# Print all features/columns\n",
    "print()\n",
    "print('Features: ', features_list)\n",
    "print()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "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>c_1</th>\n",
       "      <th>c_2</th>\n",
       "      <th>c_3</th>\n",
       "      <th>c_4</th>\n",
       "      <th>c_5</th>\n",
       "      <th>c_6</th>\n",
       "      <th>lcs_word</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0.398148</td>\n",
       "      <td>0.079070</td>\n",
       "      <td>0.009346</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.191781</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.984694</td>\n",
       "      <td>0.964103</td>\n",
       "      <td>0.943299</td>\n",
       "      <td>0.922280</td>\n",
       "      <td>0.901042</td>\n",
       "      <td>0.820755</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0.869369</td>\n",
       "      <td>0.719457</td>\n",
       "      <td>0.613636</td>\n",
       "      <td>0.515982</td>\n",
       "      <td>0.449541</td>\n",
       "      <td>0.382488</td>\n",
       "      <td>0.846491</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0.593583</td>\n",
       "      <td>0.268817</td>\n",
       "      <td>0.156757</td>\n",
       "      <td>0.108696</td>\n",
       "      <td>0.081967</td>\n",
       "      <td>0.060440</td>\n",
       "      <td>0.316062</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0.544503</td>\n",
       "      <td>0.115789</td>\n",
       "      <td>0.031746</td>\n",
       "      <td>0.005319</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.242574</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>0.329502</td>\n",
       "      <td>0.053846</td>\n",
       "      <td>0.007722</td>\n",
       "      <td>0.003876</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.161172</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>0.590308</td>\n",
       "      <td>0.150442</td>\n",
       "      <td>0.035556</td>\n",
       "      <td>0.004464</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.301653</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>0.765306</td>\n",
       "      <td>0.709898</td>\n",
       "      <td>0.664384</td>\n",
       "      <td>0.625430</td>\n",
       "      <td>0.589655</td>\n",
       "      <td>0.553633</td>\n",
       "      <td>0.621711</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>0.759777</td>\n",
       "      <td>0.505618</td>\n",
       "      <td>0.395480</td>\n",
       "      <td>0.306818</td>\n",
       "      <td>0.245714</td>\n",
       "      <td>0.195402</td>\n",
       "      <td>0.484305</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>0.884444</td>\n",
       "      <td>0.526786</td>\n",
       "      <td>0.340807</td>\n",
       "      <td>0.247748</td>\n",
       "      <td>0.180995</td>\n",
       "      <td>0.150000</td>\n",
       "      <td>0.597458</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        c_1       c_2       c_3       c_4       c_5       c_6  lcs_word\n",
       "0  0.398148  0.079070  0.009346  0.000000  0.000000  0.000000  0.191781\n",
       "1  1.000000  0.984694  0.964103  0.943299  0.922280  0.901042  0.820755\n",
       "2  0.869369  0.719457  0.613636  0.515982  0.449541  0.382488  0.846491\n",
       "3  0.593583  0.268817  0.156757  0.108696  0.081967  0.060440  0.316062\n",
       "4  0.544503  0.115789  0.031746  0.005319  0.000000  0.000000  0.242574\n",
       "5  0.329502  0.053846  0.007722  0.003876  0.000000  0.000000  0.161172\n",
       "6  0.590308  0.150442  0.035556  0.004464  0.000000  0.000000  0.301653\n",
       "7  0.765306  0.709898  0.664384  0.625430  0.589655  0.553633  0.621711\n",
       "8  0.759777  0.505618  0.395480  0.306818  0.245714  0.195402  0.484305\n",
       "9  0.884444  0.526786  0.340807  0.247748  0.180995  0.150000  0.597458"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# print some results \n",
    "features_df.head(10)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Correlated Features\n",
    "\n",
    "You should use feature correlation across the *entire* dataset to determine which features are ***too*** **highly-correlated** with each other to include both features in a single model. For this analysis, you can use the *entire* dataset due to the small sample size we have. \n",
    "\n",
    "All of our features try to measure the similarity between two texts. Since our features are designed to measure similarity, it is expected that these features will be highly-correlated. Many classification models, for example a Naive Bayes classifier, rely on the assumption that features are *not* highly correlated; highly-correlated features may over-inflate the importance of a single feature. \n",
    "\n",
    "So, you'll want to choose your features based on which pairings have the lowest correlation. These correlation values range between 0 and 1; from low to high correlation, and are displayed in a [correlation matrix](https://www.displayr.com/what-is-a-correlation-matrix/), below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "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>c_1</th>\n",
       "      <th>c_2</th>\n",
       "      <th>c_3</th>\n",
       "      <th>c_4</th>\n",
       "      <th>c_5</th>\n",
       "      <th>c_6</th>\n",
       "      <th>lcs_word</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>c_1</th>\n",
       "      <td>1.00</td>\n",
       "      <td>0.94</td>\n",
       "      <td>0.90</td>\n",
       "      <td>0.89</td>\n",
       "      <td>0.88</td>\n",
       "      <td>0.87</td>\n",
       "      <td>0.97</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>c_2</th>\n",
       "      <td>0.94</td>\n",
       "      <td>1.00</td>\n",
       "      <td>0.99</td>\n",
       "      <td>0.98</td>\n",
       "      <td>0.97</td>\n",
       "      <td>0.96</td>\n",
       "      <td>0.98</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>c_3</th>\n",
       "      <td>0.90</td>\n",
       "      <td>0.99</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.00</td>\n",
       "      <td>0.99</td>\n",
       "      <td>0.98</td>\n",
       "      <td>0.97</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>c_4</th>\n",
       "      <td>0.89</td>\n",
       "      <td>0.98</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.00</td>\n",
       "      <td>0.99</td>\n",
       "      <td>0.95</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>c_5</th>\n",
       "      <td>0.88</td>\n",
       "      <td>0.97</td>\n",
       "      <td>0.99</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.00</td>\n",
       "      <td>0.95</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>c_6</th>\n",
       "      <td>0.87</td>\n",
       "      <td>0.96</td>\n",
       "      <td>0.98</td>\n",
       "      <td>0.99</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.00</td>\n",
       "      <td>0.94</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lcs_word</th>\n",
       "      <td>0.97</td>\n",
       "      <td>0.98</td>\n",
       "      <td>0.97</td>\n",
       "      <td>0.95</td>\n",
       "      <td>0.95</td>\n",
       "      <td>0.94</td>\n",
       "      <td>1.00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           c_1   c_2   c_3   c_4   c_5   c_6  lcs_word\n",
       "c_1       1.00  0.94  0.90  0.89  0.88  0.87      0.97\n",
       "c_2       0.94  1.00  0.99  0.98  0.97  0.96      0.98\n",
       "c_3       0.90  0.99  1.00  1.00  0.99  0.98      0.97\n",
       "c_4       0.89  0.98  1.00  1.00  1.00  0.99      0.95\n",
       "c_5       0.88  0.97  0.99  1.00  1.00  1.00      0.95\n",
       "c_6       0.87  0.96  0.98  0.99  1.00  1.00      0.94\n",
       "lcs_word  0.97  0.98  0.97  0.95  0.95  0.94      1.00"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "# Create correlation matrix for just Features to determine different models to test\n",
    "corr_matrix = features_df.corr().abs().round(2)\n",
    "\n",
    "# display shows all of a dataframe\n",
    "display(corr_matrix)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## EXERCISE: Create selected train/test data\n",
    "\n",
    "Complete the `train_test_data` function below. This function should take in the following parameters:\n",
    "* `complete_df`: A DataFrame that contains all of our processed text data, file info, datatypes, and class labels\n",
    "* `features_df`: A DataFrame of all calculated features, such as containment for ngrams, n= 1-5, and lcs values for each text file listed in the `complete_df` (this was created in the above cells)\n",
    "* `selected_features`: A list of feature column names,  ex. `['c_1', 'lcs_word']`, which will be used to select the final features in creating train/test sets of data.\n",
    "\n",
    "It should return two tuples:\n",
    "* `(train_x, train_y)`, selected training features and their corresponding class labels (0/1)\n",
    "* `(test_x, test_y)`, selected training features and their corresponding class labels (0/1)\n",
    "\n",
    "** Note: x and y should be arrays of feature values and numerical class labels, respectively; not DataFrames.**\n",
    "\n",
    "Looking at the above correlation matrix, you should decide on a **cutoff** correlation value, less than 1.0, to determine which sets of features are *too* highly-correlated to be included in the final training and test data. If you cannot find features that are less correlated than some cutoff value, it is suggested that you increase the number of features (longer n-grams) to choose from or use *only one or two* features in your final model to avoid introducing highly-correlated features.\n",
    "\n",
    "Recall that the `complete_df` has a `Datatype` column that indicates whether data should be `train` or `test` data; this should help you split the data appropriately."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Takes in dataframes and a list of selected features (column names) \n",
    "# and returns (train_x, train_y), (test_x, test_y)\n",
    "def train_test_data(complete_df, features_df, selected_features):\n",
    "    '''Gets selected training and test features from given dataframes, and \n",
    "       returns tuples for training and test features and their corresponding class labels.\n",
    "       :param complete_df: A dataframe with all of our processed text data, datatypes, and labels\n",
    "       :param features_df: A dataframe of all computed, similarity features\n",
    "       :param selected_features: An array of selected features that correspond to certain columns in `features_df`\n",
    "       :return: training and test features and labels: (train_x, train_y), (test_x, test_y)'''\n",
    "    \n",
    "    # get the training features\n",
    "    train_x = (features_df[selected_features].iloc[complete_df.index[complete_df[\"Datatype\"] == 'train']]).values\n",
    "    # And training class labels (0 or 1)\n",
    "    train_y = ((complete_df[complete_df[\"Datatype\"] == 'train'])['Class']).values\n",
    "    \n",
    "    # get the test features and labels\n",
    "    test_x = (features_df[selected_features].iloc[complete_df.index[complete_df[\"Datatype\"] == 'test']]).values\n",
    "    test_y = ((complete_df[complete_df[\"Datatype\"] == 'test'])['Class']).values\n",
    "    \n",
    "    return (train_x, train_y), (test_x, test_y)\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Test cells\n",
    "\n",
    "Below, test out your implementation and create the final train/test data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Tests Passed!\n"
     ]
    }
   ],
   "source": [
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "test_selection = list(features_df)[:2] # first couple columns as a test\n",
    "# test that the correct train/test data is created\n",
    "(train_x, train_y), (test_x, test_y) = train_test_data(complete_df, features_df, test_selection)\n",
    "\n",
    "# params: generated train/test data\n",
    "tests.test_data_split(train_x, train_y, test_x, test_y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## EXERCISE: Select \"good\" features\n",
    "\n",
    "If you passed the test above, you can create your own train/test data, below. \n",
    "\n",
    "Define a list of features you'd like to include in your final mode, `selected_features`; this is a list of the features names you want to include."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Training size:  70\n",
      "Test size:  25\n",
      "\n",
      "Training df sample: \n",
      " [[0.39814815 0.         0.19178082]\n",
      " [0.86936937 0.44954128 0.84649123]\n",
      " [0.59358289 0.08196721 0.31606218]\n",
      " [0.54450262 0.         0.24257426]\n",
      " [0.32950192 0.         0.16117216]\n",
      " [0.59030837 0.         0.30165289]\n",
      " [0.75977654 0.24571429 0.48430493]\n",
      " [0.51612903 0.         0.27083333]\n",
      " [0.44086022 0.         0.22395833]\n",
      " [0.97945205 0.78873239 0.9       ]]\n"
     ]
    }
   ],
   "source": [
    "# Select your list of features, this should be column names from features_df\n",
    "# ex. ['c_1', 'lcs_word']\n",
    "selected_features = ['c_1', 'c_5', 'lcs_word']\n",
    "\n",
    "\n",
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "\n",
    "(train_x, train_y), (test_x, test_y) = train_test_data(complete_df, features_df, selected_features)\n",
    "\n",
    "# check that division of samples seems correct\n",
    "# these should add up to 95 (100 - 5 original files)\n",
    "print('Training size: ', len(train_x))\n",
    "print('Test size: ', len(test_x))\n",
    "print()\n",
    "print('Training df sample: \\n', train_x[:10])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question 2: How did you decide on which features to include in your final model? "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlcAAAGjCAYAAAAB0DAWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdeXgUVdbA4d/pToBAErKHhLDHgGyyI6BsyjaCijo4Ojrq6CCgM584royK4i4MjqiIuCGOC+q4ooKCoAKyyirKFraQjZCQhGwk3ff7ozohHQIBknTT5LzPU0/SVaeqz610ktP33qoWYwxKKaWUUqpm2LydgFJKKaXUuUSLK6WUUkqpGqTFlVJKKaVUDdLiSimllFKqBmlxpZRSSilVg7S4UkoppZSqQVpcKXUOEpE5IlIr91kRkZtFxIjIwNo4vlJK+TotrpSqhIg0FJG7ROQnEckUkWIRSRORr13FhZ+3c6xNIjJQRB4VkRBv51IZEWnpKvBeOknMHhHZ4qF8zurzpZTyLC2ulKpAROKB9cDzQCHwNDAWmA74A28BT3ktQc8YCEwGKisW3gECgB89mdBZbiAnPl9KqTrmnH73rdTpEpEAYD7QGrjaGPNJhZBnRaQn0LMGnzPIGJN7gm12oL4xJr+mnq+6jDEOwOHtPJRS6mylPVdKubsNaAv8u5LCCgBjzBpjzMzy60TkShFZLiJHXMtyEbmi4r6uoaqlItJVRBaKSDawybWtdC7TpSLysIjswuo5G1Nu/x4i8qmIZIhIkYhsE5F/ncowpYi0E5GZIvKriOSKSL6IrBORv1WIm4PVCwOw25WTEZFHK+Q5sMJ+ESLysojsF5Gjrq8vi0h4hbjS/QeLyD0issvVlu0iclNV7aguETlPRN4RkRRXnntEZKqINKoQV1Pn61HX4/Yi8h/X8+aJyGIRaeuKuUpEfhGRAlc+YyvJ+1oR+UJE9rnOV4aIfCYinSuJLX2ddROR712vyUwReVtEomrmTCqlTkR7rpRyd43r6+xT3UFEJgAvA78DTwAGuBn4TERuN8ZUPFZz4HvgI+B/QGCF7dOwhh9fA3KAba7n+QPwKbAT+DeQCfQBpgBdgD9WkepAoD9Wz9xuoJFrn9kiEmGMedoV9yoQDIwGJgIZrvWbTnIOGgMrgHjgTeAXoCswHhgsIr0q6Z17Cmt48VWgyBU7R0R2GmOWV9GWUg1EJOIE24578ygi3bHO/WHX8x4ALgD+AfQTkQHGmGJX+EBq9ny9DRxxtTsS+CewUEQeBp4DXsE6d7cCr4rIVmPMsnL734n1M58NpAJtsIarl4tIN2PMjgrPFwcsxnqNfQx0A/4K9BCRnmdTb6hS5xxjjC666OJagENAzmnEh2L9w9wJBJdbHwzsAnKBkHLr92AVX7dVcqybXdu2AQ0rbGuA9Q/1R8CvwraJrv0Glls3x/r1dotrVMlz2oClQDbgX279o65jtjxJnuWf70nXugkVYu9wrX+8kv3XA/XKrW+KVWS9fwrnvaXrGFUtWyrstxGrCA6qsH60K/7mWjhfpdu+BKTc+n+41ucCzcutj8TqsXy/wnEqy+d81zmbWWF96evsrhO8Vh7w9u+aLrqcy4sOCyrlLhirt+hUDcHq0ZhhjCnbz/X9i1i9UpdW2CcTa1L8ibxiju9VGAJEu/YLcQ3BRbh6bb52xQw9WaLGmLzS70WkgWu4Lgz4Fqvd7U62fxVGAwc5vsfvVayenNGV7DPTGHO0XH4HgO3AeafxvJ9jnZvKlrTygSLSCegMvAfUr3AOlwF5lDuHtXC+Zhhjyt8e46fSNhhj9pV73oNYBbbbeSjNRyzBrrxLY3tX8nw5WL1h5c10ra/s56GUqiE6LKiUuxwg6DTiW7m+/lrJttLbALSusH6XsSaFn8j2Stad7/r65kn2iz7JNkQkEKsXZQzQrJKQ0JPtX4VWwFpjTEn5lcaYEhHZhjUkVVFiJesOAS1O43mTjDGLKtsgIoUVVpWew8dcS2XKzmEtnK+K7c1yfd1dSWwWFc6DiHQFHscarmxUIb6yYyQaY4rKrzDGFIlIIse/JpVSNUiLK6XcbQH6i0hrY0xl//wrkjN4jqrmulS2vfR57gU2nGC/5CqO+x4wEqt36UesHrQS4A9Yw0We7sk+UYF5Juf0VJQe99/AghPEZJX7vqbP14naW+V5EJHmrhxysAqsbVg9bQb4D8fP28O17aTHVUrVDi2ulHL3P6xJzLcBk04hfpfrawesycPltXd9PZUirSqlk5XzTtRTczJi3dxyJPCOMWZchW0Vhy3hxP+YTyQRaCsifuV7r1xXMSZQM+egukrPoaOqc+iB83W6RmMVUJcbY5ZUyCcca95VRW1EpF75oVcRqY/Vy/h7bSarVF2nc66Ucvc6Vq/APVLJrRTAuuLMdYUgwHdYPQh/F5GgcjFBwN+xJrt/VwN5LQTSgQdEJKySnALKP38lSntH3HotRCQGq5Cs6Ijr63HPdQKfYU3Ernisv7nWf3qKx6lN67F6JseJyHHDYiLiV+7c1vb5Ol0nyudvQJMT7BMMTKiwboJr/Wc1mp1Syo32XClVjjEmX0RGAl9h3UrhW6zi6BBWkTAIGIZ16TzGmMMich/WrRhWue55BNYVcfHA7caY7BrIK09E/oL1T3GbiLyJdYViCNbE6quwejeWnmD/XFdbbhCRAmAN1pye27Hm64RX2GWl6+uzIvIu1tVrW4wxJ/o4meewblPwsoh0wypkumLdVmCba7tXGWOMiNyIdSuGTa5z+CvQEOtndRXwIDDHA+frdH2DNVz8jlgf+ZMF9MMaotxF5X/LdwGTRaQjsA7ojnUrht+BGTWUl1KqElpcKVWBMWana/Lw7cDVwL+whmQygbXATVjzcUrjZ4pICtZ8qNKbSW4ERhtjaqyHwBizUKy7wz8A3IBV7GVh/ROdzknuQ+VyA/AMMMrVhh1YbSumwtWLxpjlInI/MA7rflt+WJPAKy0WjDHZItLPFXM5cAvW1XqzgMnmBHeg9zRjzAbXz/ZBrDzHYd0KYQ/W7SvKD+3W2vk6g7x3icgIrHtkTcLqyVoODABewro1RUVJWJPxpwHXAUeBd4F7yl8JqZSqeeJ+ZbBSSilfJyJ7gD3GmIFeTkWpOknnXCmllFJK1SAtrpRSSil1ThCRN0UkXUQqHZJ33YR3hojsFJFNrjmipdtuEpEdruWmcuu7i8hm1z4zRKTK25locaWUUkqpc8UcYPhJto/A+vSD87A+m/MVANeVwpOxPu2gF9bFIKU3Cn7FFVu638mOD2hxpZRS5xxjTEudb6XqImNM6Q1/T+QKYK6xrMT6OLEYrKvAvzPGZBpjsrCuEh/u2hZsjPnZ9fFVc4Erq8pDiyullFJK1RVNgf3lHie51p1sfVIl609Kb8XgUq/rX+vsZZNTX3rI2yl41djusd5Owav8clK9nYJX2et4+4uTdlUddI7qu7ji7crqnnWPD/foxyFV939t8Ya3bscaois12xhT8QPjT6ay9pozWH9SWlwppZRSyiPEZq/W/q5C6nSKqYqScP8g9jisz2VNwvpQ9PLrl7rWx1USf1I6LKiUUkqpuuIL4C+uqwYvBLKNMSlYHzE2VERCXRPZhwILXdtyReRC11WCfwE+r+pJtOdKKaWUUh5R3Z6rKo8v8j5WD1SEiCRhXQHoD2CMmQV8jfWxUTuxPlLqFte2TBF5HOujrgCmGGNKJ8aPx7oKMQDro6i+qSoPLa6UUkop5RG1XVwZY66rYrsB7jjBtjeBNytZvxboeDp5aHGllFJKKY+o7eLqbKFzrpRSSimlapD2XCmllFLKI8ReN3qutLhSSimllEfY6siwoBZXSimllPKIujLnSosrpZRSSnlEXSmudEK7UkoppVQN0p4rpZRSSnmE2OpGn44WV0oppZTyiLoyLKjFlVJKKaU8QosrpZRSSqkaVFeKq7ox+KmUUkop5SHac6WUUkopj9A7tCuPmD35Fv7Q/wIOZubQ9Y+PeDudGrdv81qWvTcLp3HS/uLhdLtsjNv23Iw0vn/reQpys2nQKIhL/3YvgWGRZduPFuTx/r9up1W3vvS/YYKn069Ry5cv57lnn8XpdDJ69Gj+euutbtuTk5N5dPJksrKyCG7cmKeeeoro6GgvZVszlq1awzMvzMLhdHD1yBHcdsO1bttT0tKZ9ORUco/k4XA4mTjur/Tv04vi4mIem/oCv27bgYjwwP+Np1fXC7zUijP309pNPDX7XZxOJ9cMHcDfxox0256cfogHp88mNy8fh9PJ3TePYUDPCyguKeHhGW+ydedeHA4HV1zSj7FjRnmpFWdu2ZadPPvhQpxOJ1dd1JVbh1/ktj0lM5uH3vqM3IIiHE4nd42+hIs7ncdXqzYz59sVZXHbD6Qx719jadesiaebcMb6xEdwz2XnYxf4bF0Sc37a7ba9SeMGTB7didBG9cguKObhjzeSnlNEj1Zh3D2iXVlcy4hGTPpoI0t/S/d0E2pFXRkW1OLKy+Z+uZyZ8xbz1uO3eTuVGud0Ovjxvy8z6p9PERgWwcdT/o+WXXoT1rRFWcyKD1+nbd9LaNdvCEm/bWDl/+Zw6d/uLdu+6tN3iG3byRvp1yiHw8HTTz3FrFdfJTo6mj9ffz0DBg6kTZs2ZTHTp09n5KhRXH755axetYoZL7zAk0895cWsq8fhcPDE9Jd57fmnaRIZwbV/+zuD+l1Im1bHfv6vvv0ewwb150+jR7Fr917G3/cw3340l4+//AaAT99+lUNZhxl/z7/44LUXsfnQZdwOh5PHX5nLG0/cR3REGGMmPsqgC7sS37xpWcysDz5n+MW9uO6yS9i57wC3T57O4rf+zcJlazhaXMIXM5+koLCIkeMncdmAC2kaHXmSZzy7OJxOnnr/G2bfdQPRocFc9/TrDOzcljaxx9ow+6ufGNqjA9cO6MGu5IPc8dJ7LOj0f1zWuxOX9bZ+77cfSOP/Zs7zqcLKJvDAqPZMmLOGtJxC3hnXhx9+T2f3wbyymInD2/HVhgPM35BMz1Zh3DkkgUf+t5m1uzO5fqZVWAYH+PPZXRezcmeGt5pS4+pKceU7f6lOkYjEiMgXIpIsIkZEWno7p5NZ9st2srLzqg70QemJ22kcFUvjqBjsfv7E9x7A7g0r3WIyk/cRd34XAJq2u4Dd638+tv+eHRTkZNGsQzeP5l0btmzZQrNmzYiLi8Pf359hw4ezdOlSt5jEXbvo3bs3AD179Tpuu6/Z/Ns2mjeNpVlsDP7+/oy4ZCDfL/vZLUZEyMvPByA3L4/IiDAAdu3ZR+/uXQEIDw0hKDCQX3/f7tkGVNOm7Yk0j42mWUwU9fz9+EP/3ny/8he3GBHhSH4hALl5BUSFhVjrgYLCIkocDgqPFuPvZ6dRwwBPN6Fatuw+QPOoUOIiQ/H3szO8RweWbNzmFiMCeQVFABwpKCSycdBxx/lm9RZG9OzokZxrSoe4EPYfyudAVgElDsO3m1MZeL57L3SrqEasTjwEwJrdmQxod3wv9SUdolmxI4PCYqdH8lY155wrrgAnsAC42tuJ1HV5hzPchvgCQyPIyzrkFhPRrDW71i0HIPGXFRQXFlB4JAfjdLJi3mv0GXNu9Oilp6fTpMmxd97RUVGkp6W5xSS0bcuiRYsA+H7xYvLy8jh8+LBH86xJ6QcP0STq2M8/OjKC9Az3d+ATbrmB+d9+zyVX/ZkJ9z7MpLvuAKBtfGuWLPuZkhIHScmpbN2+g9T0gx7Nv7rSD2XRxFUsAkRHhJF2KMst5o4/j+bLJSsY+Je7GDf53zw07gYAhl7Uk4AG9el/w/9xyc0T+etVIwgJCvRo/tWVdjiX6NDGZY+jQ4NJP5zrFjN+1ADmr9rMpfc/z4SX3ufBPw0/7jgL1271ueIqKrg+adkFZY/TsguJDKrvFrMjNZdL2lt/Ewa1jyawgR+NA/zdYoZ1imHhppTaT9iDbDZ7tRZf4RPFlYg0E5FPROSgiBwSkZdOFGuMSTPGzATWeDBFVQljjl8n4v6475jbSN62mQ8fvYPkbZtpFBqO2OxsWTKfFp17EhTmO8MgJ2MqORlS4WTcfffdrFu7lmvHjGHtunVERUVh9+HJn4ZK2ox7m79etJQrRgxh8SfvMnPq4zz4+HPWnLQ/DCM6MoJr/3Ynz774Cl06tve5c1Hpz7xi+39YyehLL2Lp3P8w67F/cv+/Z+N0Otm8PRG7zcYP7/yH7978N299uoD9Kb4/56bCrz/frN7CFX0vYNGzE5l553VMeusznM5j523T7iQa1PPnvKZRnk20miq2Ezjut+H5Bdvo1iqUdyf0pXvLUNKyC3GUa3tEYH3io4P4+RwaEgRrWLA6i6846+dciYgdmA98D9wIOIAeNXTsscBYAHtcX2wRbWvisMolMDSCI5nHehuOZGXQMCTcLaZRaDgj7nwYgOLCAhLXLaN+w0ak7vqNlO2/suX7+RQXFeIoKca/fgP6/PGvHm1DTYmOjiY1NbXscVp6OpFR7v8woqKimP788wDk5+ezeNEigoKOHybxFdGREW69TWkHM4iMcP/5f/LVAmZNexKALh3bc/ToUbKycwgPDeH+f4wri/vz+LtoEdcUXxIdEUZqRmbZ47SMTKLCQ9xiPv72B16bcg8AXc+Pp+hoMVk5R5i/dCUXde+Ev58f4SHBdGt/Hlt27qZZjO8UGdEhQaRlZZc9TsvKITLE/fX86fINvPKP6wG4oE0ziopLyDqST3hwIwAWrPmVET07eC7pGpKWU0R042PDuNGNG5CRW+QWk5FbxL3vbwAgoJ6dwe2bcKSopGz7kI5NWLI1jRJnJe9SfZgvFUjV4Qs9V72AWOBeY0yeMabQGLOsJg5sjJltjOlhjOmhhVXNi2qVQHZaMjkHU3GUFLNz1Q+06nKhW0xBbjbGac0nWPfVPM6/aCgAQ8bez1+mzeXGqW/Td8xttO17qc8WVgAdOnRg3759HEhKori4mIULFjBgwAC3mKysLJyuc/HGG29w5ZVXeiPVGtOxXVv2JR0gKTmV4uJivlm8lEEXuf/8Y6KjWLXO+geza88+io4eJSykMQWFheQXWHORVqxZh5/d7jYR3hd0SmjF3gNpJKUe5GhxCV//uIpBvbu6xcRGhrNyw1YAdu1Lpqi4mLDGQcREhrNq41aMMeQXFrHx9120jovxRjPOWIeWTdmbnklSRhbFJQ4WrP2VgRckuMU0CQtm1e/WVXSJKdZ5CgtqCIDTafh2ne8NCQJsPZBNs/CGxIYE4GcXhnZqwg+/u/c8hjT0L+vJv6V/a774Jclt+7DOMSzcfG4NCYL2XJ1NmgF7jTElVUb6oHeevp3+3dsSERJI4oJpTJn1OXM++8nbadUIm93OxTeM58vpD2GcDtpdNJSwpi1Y/elcIlsm0KrrhSRv28TKj+eACLEJHX3+dgsn4ufnxwMPPsj48eNxOp1cceWVxMfHM/Pll2nfoQMDBw5k7dq1zJgxAwG6d+/Og5MmeTvtavHzszNp4h3c/s9JOJxORl82lPhWLXnp9bfp0C6BQRf14d47xjL5uf8w98NPEBGemHQPIkJm1mFu/+e/EJsQHRHO0w/d5+3mnDY/u52Hxt/IbQ9PtW5FMKQ/57WIY8Y7n9DxvJYMvrAb9912HY/MeJO3P1+IIDw98TZEhOtHXsK/nn+dURMmgYHRQy6mbavm3m7SafGz25j0pxGMf+FdHE7Dlf26EB8bxctfLKF9i1gGXdCWe64ZymP//ZJ3Fq9CgMdvvqJsuHzdjr1EhwYTFxnq3YacAYfT8Nz8rbx0Uw/sNuHzX5JITD/CuMHxbE3O5sffD9LddYWgMbB+TybPzN9atn9MSADRjRuwbk/mSZ5Fnc2ksnkBZxMR6QN8DsSeToElIn5AMdDKGLOnqvh6Xf96dp+IWjT1pYe8nYJXje0e6+0UvMovJ7XqoHOYvY63vzhpl7dT8Jq+i8OrDjrHrXt8eGVTxGpN3J/fqNb/2qR3b/VovmfKF4YFVwMpwDMi0khEGohIv5PtICINgNJLM+q7HiullFLKi3RY8CxhjHGIyChgBrAP66KL94DlJ9mtoNz3v7u++kS1q5RSSp2r9ONvziLGmH3AKc/uNcZoIaWUUkqdZXyp96k6fGFYUCmllFLKZ/hkcSUis0TkSCXLLG/nppRSSqnK6Zyrs5gxZhwwrspApZRSSp01fKlAqg6fLK6UUkop5XtstroxJdonhwWVUkoppc5W2nOllFJKKY+QOtJzpcWVUkoppTyi9OONznU6LKiUUkopj7DZpFpLVURkuIhsE5GdIvJAJdtbiMhiEdkkIktFJM61fpCIbCi3FIrIla5tc0Rkd7ltXarKQ3uulFJKKeURtTksKCJ24GVgCJAErBGRL4wxW8uFTQPmGmPeFpHBwNPAjcaYJUAX13HCgJ3At+X2u9cY8/Gp5qI9V0oppZQ6F/QCdhpjEo0xR4EPgCsqxLQHFru+X1LJdoBrgG+MMflnmogWV0oppZTyCLFJ9RaRsSKyttwyttzhmwL7yz1Ocq0rbyNwtev70UCQiIRXiPkT8H6FdU+6hhKfF5H6VbVTiyullFJKeYRNpFqLMWa2MaZHuWV2ucNXNuZoKjy+BxggIuuBAcABoKTsACIxQCdgYbl9HgTaAT2BMOD+qtqpc66UUkop5RG1fCuGJKBZucdxQHL5AGNMMnAVgIgEAlcbY7LLhYwBPjXGFJfbJ8X1bZGIvIVVoJ2U9lwppZRSyiOqOyxYhTXAeSLSSkTqYQ3vfeH2/CIRIlJa+zwIvFnhGNdRYUjQ1ZuFWPeRuBLYUlUiWlwppZRSyucZY0qAO7GG9H4DPjTG/CoiU0TkclfYQGCbiGwHooEnS/cXkZZYPV8/VDj0uyKyGdgMRABPVJWLDgsqpZRSyiNq+7MFjTFfA19XWPdIue8/Biq9pYIxZg/HT4DHGDP4dPPQ4spl6ksPeTsFr7n3ziqL8HPalJg23k7Bq4Lj2no7Ba8Kiw70dgpeFRrdwtspeE1mWqq3U6hzpI6Ml2lxpZRSSimP0I+/UUoppZRSp017rpRSSinlEbU95+psocWVUkoppTyilu9zddbQ4koppZRSHqHFlVJKKaVUDbLphHallFJKKXW6tOdKKaWUUh6hw4JKKaWUUjVIiyullFJKqRqkt2JQSimllKpBeod2pZRSSil12rTnSimllFIeoR/crJRSSilVg3TOlVJKKaVUDaorVwvWkQ46pZRSSinP0J4rpZRSSnlEXblaUIsrpZRSSnmEzrlSNWbf5rUse28WTuOk/cXD6XbZGLftuRlpfP/W8xTkZtOgURCX/u1eAsMiy7YfLcjj/X/dTqtufel/wwRPp1+rZk++hT/0v4CDmTl0/eMj3k6nxg3u3pqnbx+KzSb8d+EGXvjoZ7ftcVHBvHjXSMIbNyQrt5DxUz8n+VAuAJNvGcSQnvEATPtgGZ/9+JvH86+u/h2imXxtF2w2Yd6y3cxasM1te9Owhjx7Uw/Cg+pxOK+YiW+sJvVwAQAPXN2JQZ2aYBNh2dY0Hpu30RtNqJY+8RHcc9n52AU+W5fEnJ92u21v0rgBk0d3IrRRPbILinn4442k5xTRo1UYd49oVxbXMqIRkz7ayNLf0j3dhDPWs3kIEy5ujU2Eb7am8cEvSW7bo4Lqc8/g8wgJ8Ce3qISnv91GRt5Ra1tgfe4eHE9kYH0AJn35K2m5RR5vQ3XU9df+ieicK1UjnE4HP/73ZS6b+DjXPfEqO1YtJfPAXreYFR++Ttu+l/CnKa/Q4/LrWfm/OW7bV336DrFtO3kwa8+Z++VyRt4x3dtp1AqbTXhuwnDGPPIBfce9ylUDOtC2WYRbzJRbL2Xe4s30v+N1pr3/Ew/fMgiAIT3j6RzfhAF3vs7QiXP4+9V9CAqo541mnDGbwJTru3LzjGUMnbyQy3s2Iz4myC1m0h8788nKvYyYsogZ87dy31UdAejWOpzubcIZ8dh3DHv0Wzq3DKN3QmRlT3PWsgk8MKo9/5i7lmteXMawzjG0imzkFjNxeDu+2nCAP728nNeX7OTOIQkArN2dyfUzV3D9zBWMe2sNhcUOVu7M8EYzzohN4O8D2jDpy1+59b1fGJQQSfPQALeY2/u14rtt6Yz9YD3vrN7HrX1alm27f0gCH64/wK3v/cIdH23gcEGxh1tQPXX9tX8ydptUa/EV51xxJSKXicgyETksIqki8pqIBFW9Z+1IT9xO46hYGkfFYPfzJ773AHZvWOkWk5m8j7jzuwDQtN0F7F5/rHcjfc8OCnKyaNahm0fz9pRlv2wnKzvP22nUim4JsexOzmRv6mGKS5x8+uNWRvRJcItp2zyCHzfsAeCnjXsZcWFC2foVm/fhcBryi4rZkpjG4B5tPN2EarmgVRh704+wPyOPYofhyzX7GXJBrFtMfEwQK1y9MT9vO8ilru0GQ31/O/5+Nur52/GzCxk5hR5vQ3V0iAth/6F8DmQVUOIwfLs5lYHnR7vFtIpqxOrEQwCs2Z3JgHbRxx3nkg7RrNiRQWGx0yN514S20UEkZxeSklNEidOwdMdB+rUOd4tpERrA+v2HAdhwIJu+rcMAaB4agF3gF9e2wmInRSW+03bQ1746B4sroDHwBBALnA/EAVO9lUze4Qy3Ib7A0Ajysg65xUQ0a82udcsBSPxlBcWFBRQeycE4nayY9xp9xtzm0ZxVzYgJD+JARm7Z4+SMHGLC3ev8LbvTGHmRNfwzsm9bghrWJzQogC2JaVzSow0B9f0ICw7gos4taBoR7NH8q6tJSAApmQVlj1MPF9CkQu/Fb/uzGd6tKQDDusYSFOBPSKN6rE/MZOW2g6yeOpJVU0fy09Y0dqXm4kuiguuTln2s/WnZhUQG1XeL2ZGayyXtmwAwqH00gQ38aBzg7xYzrFMMCzel1H7CNSiiUT3Syw3jHTxSRHgj957XxEN5XNzG6sm9qHU4jer5EdzAj7iQAI4cdTB5RDtmXduFsX1b4kMdFoC+9k9Ge67OIiLSTEQ+EZGDInJIRF46Uawx5j1jzAJjTL4xJgt4DejnuWwr5nP8uooXS/QdcxvJ2zbz4aN3kLxtM41Cw345w+EAACAASURBVBGbnS1L5tOic0+Cws6dLuG6pLKLYkyFF8Tk1xfTr2Nzlrx4K307NSc5I4cSh5Ol63ezaM0uvpl2M6/dfyVrfj9AidO33r1X3n73x099vIneCZHMf+gSeidEkpKVj8PppEVkI9rEBNHn/q/oc998+rSNotd5Eccf8CxW2b+Bin8Onl+wjW6tQnl3Ql+6twwlLbsQh/NYVERgfeKjg/jZh4YEofK2V/Tq8j10bhrMrGu70LlpYw4eKcLhNNhtQqeYYGYv382EDzcQ07gBQyvp0Tub1fXX/snUleLqrJ/QLiJ2YD7wPXAj4AB6nMYh+gO/nuDYY4GxAGPufYK+V1xXvWQrERgawZHMg2WPj2Rl0DDEvXu8UWg4I+58GIDiwgIS1y2jfsNGpO76jZTtv7Ll+/kUFxXiKCnGv34D+vzxrzWep6p5yRm5NI041lMVGxFMauYRt5jUzCPc9OT/AGjUwJ9R/dqRm2+9458+bznT51k9mq/edwWJBzI9lHnNSMkqICbs2Lv1JiEBpB0ucItJzy5k/CxrGLxhfTvDuzUlt6CE6y5uzYbETPKLHAAs3ZJK19ZhrN7hO0VGWk4R0Y2PtT+6cQMyKkzKzsgt4t73NwAQUM/O4PZNOFJUUrZ9SMcmLNmaRomzkndpZ7GDeUeJKtdLFxlYn0OuyeqlDuUd5bFvfgeggb+Ni9uEk3fUQcaRo+zMyCMlxzpXyxMPcX50EAt86HqOuv7aPxlfKpCqwxd6rnphDfHda4zJM8YUGmOWncqOIjIEuAmo9DI0Y8xsY0wPY0yP2iisAKJaJZCdlkzOwVQcJcXsXPUDrbpc6BZTkJuNcfVKrPtqHudfNBSAIWPv5y/T5nLj1LfpO+Y22va9VAsrH7J+ezKtY8NoHt0Yfz8bo/u355uV291iwoIDyt7l3jWmH+9+a10VZLMJoUHWH+f2LaPo0DKKJb8kejT/6tq0J4uWUYHEhTfE3y6M6tmMRRvdh7dCA+uVtX/CiHZ8tHwPAAcy8+mVEIHdJvjZhd4JkexM8a2hka0HsmkW3pDYkAD87MLQTk344Xf3q/1CGvqXtf+W/q35osIVdcM6x7Bws28NCQJsS8ulaeMAmgTVx88mDDwvkhW73d8cBDfwK+vhuq57MxZsTbP2Tc8lsL4fjRtY7/27xIWwN8u9MDnb1fXX/sloz9XZoxmw1xhTUmVkOSJyIfAecI0xZntV8bXFZrdz8Q3j+XL6Qxing3YXDSWsaQtWfzqXyJYJtOp6IcnbNrHy4zkgQmxCx3Pudgsn887Tt9O/e1siQgJJXDCNKbM+Z85nP3k7rRrhcBruf2UhHz1xHXabjfe+3ci2fRk8cEN/NuxIYcGqHfTr1IKHbx6EwfDzlv3c9/ICAPztNr6aeiMAuflHGTftC7fhIl/gcBomv7+BuXddjM0mfLR8DztScph4eXs2781i0cYULkyI5N7R1lVSq7dn8Mj76wH4Zl0SfdtFsWDyEIyBH35NZbGPzTtyOA3Pzd/KSzf1wG4TPv8licT0I4wbHM/W5Gx+/P0g3VuFceeQBIyB9XsyeWb+1rL9Y0ICiG7cgHV7fKvHEsBp4MUfd/HMFR2xCSzYmsbezHxu6tWc7elH+HlPJhc0bVx2heCmA9m8+MOusn1fXb6bqVd2QgS2px/h619Tvdia01fXX/sKpOIckLONiPQBPgdiT7XAEpGuwELgVmPMl6eyzwvLE8/uE1GL7r3zCW+n4FVBMb51FV5NC45r6+0UvCosOtDbKXhVaB1u/65NvlW01Ybds6/xaHfQtXNWV+t/7bybe/lE95UvDAuuBlKAZ0SkkYg0EJETTlAXkY7AAuDvp1pYKaWUUqr21ZVhwbO+uDLGOIBRQDywD0gCrj3JLv8EIoE3ROSIa6l0QrtSSimlPKeuFFe+MOcKY8w+4MpTjL0FuKV2M1JKKaWUqpxPFFdKKaWU8n1221k/YFYjfLKVIjKr3JBf+WWWt3NTSimlVOVqe1hQRIaLyDYR2SkiD1SyvYWILBaRTSKyVETiym1ziMgG1/JFufWtRGSViOwQkXkiUuUHvfpkcWWMGWeMCaxkGeft3JRSSilVudosrlw3HX8ZGAG0B64TkfYVwqYBc40xnYEpwNPlthUYY7q4lsvLrX8WeN4Ycx6QBdxaVTt9srhSSimllO+p5Z6rXsBOY0yiMeYo8AFwRYWY9sBi1/dLKtnuRkQEGAx87Fr1NqcwB1yLK6WUUkr5BBEZKyJryy1jy21uCuwv9zjJta68jcDVru9HA0EiUvqZdA1cx1wpIqUFVDhwuNx9Nis75nF0QrtSSimlPMJe2adanwZjzGxg9gk2n8rnpd8DvCQiNwM/AgeA0sKpuTEmWURaA9+LyGYg5xSOeRwtrpRSSinlEbV8r6okrI/MKxUHJJcPMMYkA1cBiEggcLUxJrvcNowxiSKyFOgK/A8IERE/V+/VccesjA4LKqWUUsojannO1RrgPNfVffWAPwFflA8QkQgRKa19HgTedK0PFZH6pTFAP2CrsT4jcAlwjWufm7A+ku+ktLhSSimllM9z9SzdifXZwr8BHxpjfhWRKSJSevXfQGCbiGwHooEnXevPB9aKyEasYuoZY0zpJ6nfD9wtIjux5mC9UVUuOiyolFJKKY/wq+WPsDHGfA18XWHdI+W+/5hjV/6Vj1kBdDrBMROxrkQ8ZVpcKaWUUsojfOnzAatDiyullFJKeYQWV0oppZRSNaiuFFc6oV0ppZRSqgZpz5VSSimlPKKu9FxpceUytnust1PwmikxbbydglflpuzydgpeVdfbn1PHX/+ZcW29nYLXhEUHejuFOkeLK6WUUkqpGqTFlVJKKaVUDaorxZVOaFdKKaWUqkHac6WUUkopj6grPVdaXCmllFLKI7S4UkoppZSqQXWluNI5V0oppZRSNUh7rpRSSinlEXapGz1XWlwppZRSyiNsWlwppZRSStUce92orbS4UkoppZRn2HRCu1JKKaWUOl3ac6WUUkopj9AJ7UoppZRSNUgntCullFJK1SCd0K6UUkopVYN0QrtSSimllDpt2nOllFJKKY/QOVdKKaWUUjVI51ypWrF8+XKee/ZZnE4no0eP5q+33uq2PTk5mUcnTyYrK4vgxo156qmniI6O9lK21Te4e2uevn0oNpvw34UbeOGjn922x0UF8+JdIwlv3JCs3ELGT/2c5EO5AEy+ZRBDesYDMO2DZXz2428ez7+2zZ58C3/ofwEHM3Po+sdHvJ2OR9WFttf113//DtFMvrYLNpswb9luZi3Y5ra9aVhDnr2pB+FB9TicV8zEN1aTergAgAeu7sSgTk2wibBsaxqPzdvojSacsT7xEdxz2fnYBT5bl8Scn3a7bW/SuAGTR3citFE9sguKefjjjaTnFNGjVRh3j2hXFtcyohGTPtrI0t/SPd2EWlFXeq50zpUHORwOnn7qKV6eOZNPPv2UBQsWsGvXLreY6dOnM3LUKD76+GNuHzuWGS+84KVsq89mE56bMJwxj3xA33GvctWADrRtFuEWM+XWS5m3eDP973idae//xMO3DAJgSM94Osc3YcCdrzN04hz+fnUfggLqeaMZtWrul8sZecd0b6fhFed62+v6698mMOX6rtw8YxlDJy/k8p7NiI8JcouZ9MfOfLJyLyOmLGLG/K3cd1VHALq1Dqd7m3BGPPYdwx79ls4tw+idEOmNZpwRm8ADo9rzj7lruebFZQzrHEOryEZuMROHt+OrDQf408vLeX3JTu4ckgDA2t2ZXD9zBdfPXMG4t9ZQWOxg5c4MbzRDVcM5V1yJyCAR2Swih0XkkIh8KiJNvZ0XwJYtW2jWrBlxcXH4+/szbPhwli5d6haTuGsXvXv3BqBnr17Hbfcl3RJi2Z2cyd7UwxSXOPn0x62M6JPgFtO2eQQ/btgDwE8b9zLiwoSy9Ss278PhNOQXFbMlMY3BPdp4ugm1btkv28nKzvN2Gl5xrre9rr/+L2gVxt70I+zPyKPYYfhyzX6GXBDrFhMfE8QKV4/Mz9sOcqlru8FQ39+Ov5+Nev52/OxCRk6hx9twpjrEhbD/UD4HsgoocRi+3ZzKwPPdRyBaRTVideIhANbszmRAu+NHKC7pEM2KHRkUFjs9krcn2G1SrcVXnHPFFbAVGGaMCQFigR3AK95NyZKenk6TJk3KHkdHRZGeluYWk9C2LYsWLQLg+8WLycvL4/Dhwx7Ns6bEhAdxICO37HFyRg4x4e7vXLfsTmPkRVYX+Mi+bQlqWJ/QoAC2JKZxSY82BNT3Iyw4gIs6t6BpRLBH81eqOur6679JSAApmQVlj1MPF9AkNMAt5rf92QzvZr33HdY1lqAAf0Ia1WN9YiYrtx1k9dSRrJo6kp+2prErNRdfERVcn7TsY21Pyy4kMqi+W8yO1FwuaW/9PxjUPprABn40DvB3ixnWKYaFm1JqP2EPsolUa/EVPlFciUgzEflERA66eqNeOlGsMSbNGJNcbpUDiD/BcceKyFoRWfvGG2/UdNqV5VZZDm6P7777btatXcu1Y8awdt06oqKisNvttZ5bbajs96DiOZj8+mL6dWzOkhdvpW+n5iRn5FDicLJ0/W4WrdnFN9Nu5rX7r2TN7wcocZ47797Uua+uv/4rb7/746c+3kTvhEjmP3QJvRMiScnKx+F00iKyEW1iguhz/1f0uW8+fdpG0eu8iOMPeJaqrASo+Nf/+QXb6NYqlHcn9KV7y1DSsgtxOI9FRQTWJz46iJ/PsSFBu1Rv8RVn/YR2EbED84HvgRuxiqUeVezTHNgEBLvi/1ZZnDFmNjAboKCw8PjKp4ZFR0eTmppa9jgtPZ3IqCi3mKioKKY//zwA+fn5LF60iKAg93e7viI5I5emEcdyj40IJjXziFtMauYRbnryfwA0auDPqH7tyM0vAmD6vOVMn7ccgFfvu4LEA5keylyp6qvrr/+UrAJiwo71VDUJCSDtcIFbTHp2IeNnWZP8G9a3M7xbU3ILSrju4tZsSMwkv8gBwNItqXRtHcbqHb5RaKTlFBHd+Fjboxs3ICO3yC0mI7eIe9/fAEBAPTuD2zfhSFFJ2fYhHZuwZGsaJc5a/9fkUb7U+1QdvtBz1QtreO9eY0yeMabQGLPsZDsYY/a5hgUjgIeA3z2QZ5U6dOjAvn37OJCURHFxMQsXLGDAgAFuMVlZWThd71DfeOMNrrzySm+kWiPWb0+mdWwYzaMb4+9nY3T/9nyzcrtbTFhwQNk73LvG9OPdb60rgmw2ITTI+uPUvmUUHVpGseSXRI/mr1R11PXX/6Y9WbSMCiQuvCH+dmFUz2Ys2ug+xBUaWK+s/RNGtOOj5XsAOJCZT6+ECOw2wc8u9E6IZGeK7wwLbj2QTbPwhsSGBOBnF4Z2asIPv7tf7RfS0L+s7bf0b80XvyS5bR/WOYaFm8+tIcG65KzvuQKaAXuNMSVVRlZgjMkUkbeBjSLS9EyOUZP8/Px44MEHGT9+PE6nkyuuvJL4+Hhmvvwy7Tt0YODAgaxdu5YZM2YgQPfu3Xlw0iRvplwtDqfh/lcW8tET12G32Xjv241s25fBAzf0Z8OOFBas2kG/Ti14+OZBGAw/b9nPfS8vAMDfbuOrqTcCkJt/lHHTvnDrMj9XvPP07fTv3paIkEASF0xjyqzPmfPZT95OyyPO9bbX9de/w2mY/P4G5t51MTab8NHyPexIyWHi5e3ZvDeLRRtTuDAhkntHW1cIrt6ewSPvrwfgm3VJ9G0XxYLJQzAGfvg1lcU+NPfI4TQ8N38rL93UA7tN+PyXJBLTjzBucDxbk7P58feDdG8Vxp1DEjAG1u/J5Jn5W8v2jwkJILpxA9bt8a3eylNR25PSRWQ48AJgB143xjxTYXsL4E0gEsgEbjDGJIlIF6z52aUjXk8aY+a59pkDDACyXYe52Riz4aR5VDYP6GwiIn2Az4HYMymORCQO2A+EG2NO+Er1xLDg2Sruqn97OwWvyk3ZVXWQOmcFxfjWVXg1LTiurbdT8Jqw6EBvp+B16x4f7tFxuu93HqzW/9rB8ZEnzNc1jWg7MARIAtYA1xljtpaL+QiYb4x5W0QGA7cYY24UkQTAGGN2iEgssA443xhz2FVczTfGfHyqefrCsOBqIAV4RkQaiUgDEel3omARuUpE2oqITUQigenA+pMVVkoppZSqfXaRai1V6AXsNMYkGmOOAh8AV1SIaQ8sdn2/pHS7MWa7MWaH6/tkIB2rd+uMnPXFlTHGAYzCuuJvH1Y1eu1JdmkKLABygc2AExhdy2kqpZRSqgrVvRVD+av8XcvYcodvijVSVSrJta68jcDVru9HA0EiEl4+QER6AfWA8sMaT4rIJhF5XkTc76tRCV+Yc4UxZh9wSjO7jTEvAi/WbkZKKaWU8rTyV/lX4lTugnEP8JKI3Az8CBwAyqYciUgM8A5wkzGm9P4nDwKpWAXXbOB+YMrJ8vSJ4koppZRSvs9eu+NlSVgXwZWKA8rf97J0yO8qABEJBK42xmS7HgcDXwEPGWNWltun9GqKIhF5C6tAO6mzfliwMiIyS0SOVLLM8nZuSimllKpcLd+hfQ1wnoi0EpF6wJ+AL8oHiEiEiJTWPg9iXTmIK/5TYK4x5qMK+8S4vgrWKNqWqhLxyZ4rY8w4YJy381BKKaXUqTuFSelnzBhTIiJ3AguxbsXwpjHmVxGZAqw1xnwBDASeFhGDNSx4h2v3MUB/INw1ZAjHbrnwrusCOQE2cAr1h08WV0oppZRSFRljvga+rrDukXLffwwcd0sFY8x/gf+e4JiDTzcPLa6UUkop5RF15eNvtLhSSimllEfU8oT2s4YWV0oppZTyCO25UkoppZSqQXWktvLNWzEopZRSSp2ttOdKKaWUUh5hq/Qm6uceLa6UUkop5RF1ZVhQiyullFJKeYRNiyullFJKqZpTV3qudEK7UkoppVQN0p4rpZRSSnmETmivY/xyUr2dgtcEx7X1dgpelZuyy9spKC/Sn3/dddmwUd5Ooc6pK8OCWlwppZRSyiPqyoR2nXOllFJKKVWDtOdKKaWUUh5RRzqutLhSSimllGfoBzcrpZRSStWgOlJbaXGllFJKKc+oKxO960o7lVJKKaU8QnuulFJKKeURUkfGBbW4UkoppZRH1JX7XGlxpZRSSimPqCMdV1pcKaWUUsoz6spE77rSTqWUUkopj9CeK6WUUkp5hE5oV0oppZSqQTqhXSmllFKqBtWR2krnXCmllFJK1STtuVJKKaWUR+iwoFJKKaVUDdIJ7arGLFu1hmdemIXD6eDqkSO47YZr3banpKUz6cmp5B7Jw+FwMnHcX+nfpxfFxcU8NvUFft22AxHhgf8bT6+uF3ipFWemf4doJl/bBZtNmLdsN7MWbHPb3jSsIc/e1IPwoHoczitm4hurST1cAMADV3diUKcm2ERYtjWNx+Zt9EYTatXsybfwh/4XcDAzh65/fMTb6XhUXW47nPvtH9y9NU/fPhSbTfjvwg288NHPbtvjooJ58a6RhDduSFZuIeOnfk7yoVwAJt8yiCE94wGY9sEyPvvxN4/nX13nRwdyzQVNsQms2J3Jd9sPum0PbejPDd3jCKzvR/5RB2+v2c/hgmIAZlzVieTsQgCy8ot59ec9Hs6+9tSVniudc1XLHA4HT0x/mVemPcEX77zG14uWsGv3XreYV99+j2GD+vPxmzOZ9uiDPDH9JQA+/vIbAD59+1Vee/4Zpr00G6fT6fE2nCmbwJTru3LzjGUMnbyQy3s2Iz4myC1m0h8788nKvYyYsogZ87dy31UdAejWOpzubcIZ8dh3DHv0Wzq3DKN3QqQ3mlGr5n65nJF3TPd2Gl5Rl9sO53b7bTbhuQnDGfPIB/Qd9ypXDehA22YRbjFTbr2UeYs30/+O15n2/k88fMsgAIb0jKdzfBMG3Pk6QyfO4e9X9yEooJ43mnHGBBjTpSkzl+/miW+3071ZCE2C6rvFjO4Uw+q9h3l60Q6++S2Nyzs0KdtW7HDyzOIdPLN4xzlVWIF1bqqzVHl8keEisk1EdorIA5VsbyEii0Vkk4gsFZG4cttuEpEdruWmcuu7i8hm1zFnyCl0v53TxZWIvCUiRkTivZXD5t+20bxpLM1iY/D392fEJQP5fpn7OzgRIS8/H4DcvDwiI8IA2LVnH727dwUgPDSEoMBAfv19u2cbUA0XtApjb/oR9mfkUewwfLlmP0MuiHWLiY8JYsVv6QD8vO0gl7q2Gwz1/e34+9mo52/Hzy5k5BR6vA21bdkv28nKzvN2Gl5Rl9sO53b7uyXEsjs5k72phykucfLpj1sZ0SfBLaZt8wh+3LAHgJ827mXEhQll61ds3ofDacgvKmZLYhqDe7TxdBOqpWVYQzLyjnIo7ygOY/gl6TCdY4PdYmKCG7Dt4BEAth/Mo1OF7er0iYgdeBkYAbQHrhOR9hXCpgFzjTGdgSnA0659w4DJQG+gFzBZREJd+7wCjAXOcy3Dq8rlnC2uROQiwOu/kekHD9Ek6liPS3RkBOkZGW4xE265gfnffs8lV/2ZCfc+zKS77gCgbXxrliz7mZISB0nJqWzdvoPUdPeu5bNZk5AAUjILyh6nHi6gSWiAW8xv+7MZ3q0pAMO6xhIU4E9Io3qsT8xk5baDrJ46klVTR/LT1jR2peZ6NH+l1JmJCQ/iQMax39fkjBxiwt17rbfsTmPkRe0AGNm3LUEN6xMaFMCWxDQu6dGGgPp+hAUHcFHnFjSN8K3Co3GAP1n5xWWPswqKaRzg7xZz4HABXZpa7bogNpgAfzuN6tkB8LPZuG9wPP8c2Oa4oszX2USqtVShF7DTGJNojDkKfABcUSGmPbDY9f2SctuHAd8ZYzKNMVnAd8BwEYkBgo0xPxtjDDAXuLKqRHxizpWINANeAC7GKgjfN8bceZJ4P+BF4CbAqxN1DOa4dVKhc/PrRUu5YsQQbv7TNWzYspUHH3+Oz+a+yug/DCNxzz6u/dudxDaJokvH9tjtdk+lXm2V/R6YCqfjqY838dh1XbmmbwtW78ggJSsfh9NJi8hGtIkJos/9XwHwzl396XVeGqt3ZBx/UKXUWaXy3333X/7Jry/m2fHDuO7Szvy8ZR/JGTmUOJwsXb+bbgmxfDPtZg7l5LHm9wOU+NB0CDjB8FWFv32fbk7hj12a0rtFGDsz8sjKP4rDaQU98s1vZBeWEN6oHv+4uDXJ2YVk5B2t9bw9oZbnszcF9pd7nITVE1XeRuBqrJpiNBAkIuEn2Lepa0mqZP1JnfXFlaubbz7wPXAj4AB6VLHbROBHY8ymkw2NishYrK4+Zk59ktv+cn2N5FxedGSEW29T2sEMIiPC3WI++WoBs6Y9CUCXju05evQoWdk5hIeGcP8/xpXF/Xn8XbSIq/JnetZIySogJuxYT1WTkADSDhe4xaRnFzJ+ljVM2rC+neHdmpJbUMJ1F7dmQ2Im+UUOAJZuSaVr6zAtrpTyAckZuTSNONZTFRsRTGrmEbeY1Mwj3PTk/wBo1MCfUf3akZtfBMD0ecuZPm85AK/edwWJBzI9lHnNOFxQTGjDYz1VoQH+ZBcWu8VkF5bw+kpr/m09u40uscEUljjLtgEcyjvKjoNHiAsJOHeKq4rvsE93/3L/t11mG2Nml26uZJeKT3gP8JKI3Az8CBwASk6y76kc8zi+MCzYC4gF7jXG5BljCo0xy04U7Orluh2o8vIbY8xsY0wPY0yP2iisADq2a8u+pAMkJadSXFzMN4uXMuiiC91iYqKjWLVuA2DNsyo6epSwkMYUFBaSX2DNM1qxZh1+djttWrWolTxrw6Y9WbSMCiQuvCH+dmFUz2Ys2pjiFhMaWK/sncyEEe34aPkeAA5k5tMrIQK7TfCzC70TItmZosOCSvmC9duTaR0bRvPoxvj72Rjdvz3frHSfLxoWHFD2u3/XmH68+601yGCzCaFB1puy9i2j6NAyiiW/JHo0/+ram5VPZGA9whv6YxehW1wIm5Jz3GIa1bOX/dce1i6SlXuzAAjwt+PnuqSuUT07rSMakXouzTc1zmot5f9vu5bZ5Y6eBDQr9zgOSHZ7emOSjTFXGWO6Av9yrcs+yb5Jru9PeMzKnPU9V1iN3WuMKTnF+P8AU1wny+v8/OxMmngHt/9zEg6nk9GXDSW+VUteev1tOrRLYNBFfbj3jrFMfu4/zP3wE0SEJybdg4iQmXWY2//5L8QmREeE8/RD93m7OafF4TRMfn8Dc++6GJtN+Gj5Hnak5DDx8vZs3pvFoo0pXJgQyb2jrSsEV2/P4JH31wPwzbok+raLYsHkIRgDP/yayuJNKSd7Op/0ztO30797WyJCAklcMI0psz5nzmc/eTstj6jLbYdzu/0Op+H+Vxby0RPXYbfZeO/bjWzbl8EDN/Rnw44UFqzaQb9OLXj45kEYDD9v2c99Ly8AwN9u46upNwKQm3+UcdO+KBsu8xVOAx9uSOaOi1ojAiv3ZJGaW8Rl7aPZl1XA5pQczosM5PKOTcDAzow8PtxwAIAmQfW5rltTnMa64vq7bemk5hZ5uUU+Yw1wnoi0wuqR+hPg1nMiIhFApjHGCTwIvOnatBD4//buPD6q+urj+OckBAhJgISwhR3ZBFkEFEHFBVG0UmttbWm1tdW61T5qrUutexfbuvR5WrUWrVKtu62CG+4ri4CyirLvCYGQEAIECJnz/DFjSEIkYCZzM5nv+/W6L+be+7uXczKBOfP7/e69f6g0if1U4NfuXmhmJWZ2DPAx8CPC044OyKqPgzc0ZjYSmAzkHEyBZWZbgd3s67ZrDxQAV7r7k191XNmm1Q37B1GP+tw0J+gQArVh9qtBhyASmIyOgV/3E5gJF44POoTA3XfOoJjeeWr3jpI6fdY2S8s4YLxmdgbhTpZk4BF3/72Z3QHMcfcpZvYdwlcI68kOnwAAIABJREFUOuFhwZ+7++7IsT8Fboyc6vfu/mhk+3BgEpAKvAb8wmspnuKh52oWkAf80cxuJTznapi7T/uK9n2oOtyZB4wn4IntIiIiCc/r9+IEd38VeLXatlsqvX4eeP4rjn2EfT1ZlbfPAY44lDga/Jwrdy8nXBz1AtYSHv/83gHab3L3jV8ukc0F7l76VceIiIhIDLjXbYkT8dBzhbuv5SDuK/EVxybIzfZFRESkIYiL4kpEREQagXoeFmwoGvywYE3M7EEz217D8mDQsYmIiEjNzEN1WuJFXPZcufulwKW1NhQREZGGI44KpLqIy+JKRERE4lCCFFdxOSwoIiIi0lCp50pERERiI0F6rlRciYiISGyEVFyJiIiIRE08XfFXF5pzJSIiIhJF6rkSERGR2EiQnisVVyIiIhIbcfR8wLpQcSUiIiKxoZ4rERERkejRhHYREREROWTquRIREZHYSJCeKxVXIiIiEhsqrhJL8raNQYcQmKz26UGHEKhtHQ8LOoRAleStCDoECZDef4kpFVciIiIi0aMJ7SIiIiJyyNRzJSIiIrGhBzeLiIiIRJHu0C4iIiISRZpzJSIiIiKHSj1XIiIiEhOJcrWgiisRERGJDRVXIiIiIlGk4kpEREQkikLlQUcQE5rQLiIiIhJF6rkSERGRmHDdRFREREQkijQsKCIiIhJFofK6LbUws3FmtsTMlpvZDTXs72pm75rZXDNbYGZnRLb/0MzmVVpCZjYksu+9yDm/3NeutjjUcyUiIiIx4eX113NlZsnA/cBYYD0w28ymuPviSs1uAp5197+bWX/gVaC7uz8BPBE5z0BgsrvPq3TcD919zsHGop4rERERaQyOBpa7+0p33wM8DZxVrY0DLSOvWwG5NZxnAvBUXQJRz5WIiIjERh0ntJvZxcDFlTZNdPeJkdedgHWV9q0HRlQ7xW3AG2b2CyANOKWGv+Z77F+UPWpm5cB/gN+5H/gJ1CquREREJDbqOKE9UkhN/IrdVtMh1dYnAJPc/R4zGwk8bmZHuIfvbmpmI4Cd7r6o0jE/dPcNZpZBuLg6H3jsQHGquIqBD+cs4A8TnyAUCvGdU0/gZ+eeWWV/7qYt/PreiZTs2El5KMQvLziXE44aTNnevdz810dYvHwN5eXlnDXmWC4+d3xAWXw9I3tl86tvHE6ywYufrGfSh6uq7O/Qqjm3nj2QzLSmFJeWcfPz89m0bTfDe2Txy9P7VbTrnp3Gjc/N573PN8U6hTo5eVhP7rzkVJKSjH+/Po//e25Glf2d27Xkb1edSZtWLSgq2cVld00md0sJALf+5CTGHtULgLuf/ogXP/g85vHXp4m3/oQzRg9mc+E2jvzuLUGHE3PKv3Hnf3j7dL4zuBNJBtNXFfLm0s1V9me2SOG8YZ1Jb9aEnXvK+dfsdWwtLQPgr98eSG7xLgCKdpbxjxmrYxx9/fH6vVpwPdCl0npn9h/2uxAYB+DuM8ysOZANfPnh8n2qDQm6+4bInyVm9iTh4ccDFleac1XPystD/PbvjzHx9mt46e938soHM1m+dkOVNg8+PZlxxx/Nf//2W+65/nLueCD8nr3+0Wz2lO1lygO/5/n/u51nXnuPDfmba/prGqQkgxvG9+d/HpvDd/72EacN6kiPtmlV2lw9rh+vzNvA9++fxsPvLueKsX0AmLOqkB88MJ0fPDCdSx+dza6ycmYuLwgija8tKcn48+XjOPeWpxl16T/49gkD6Nslu0qbOy48hWfeXsjonz/M3U99yM0/OQmAsUf1YlCvDpxwxcOcevUkfnHOSDJSmwaRRr157KVpnPnze4MOIzDKv/Hmb8C5QzrxwLRV/O6NpQzr0poOGc2qtDl7YEdmrdnKnW8t47XP8/nmgA4V+8rKQ/zx7WX88e1ljaqwioHZQG8z62FmTQkXSlOqtVkLjAEws8OB5sDmyHoS8F3Cc7WIbGtiZtmR1ynAmcAiatHoiiszOzFyCeX2SsuPg4pnwdKVdM1pT5eO7Wia0oQzRo/gnZmfVmljZmzfGf6WUrKjlHZZrcPbgdJdu9lbXs6uPWWkNEkmrUVqrFP42gZ0bs26LTvZUFTK3nLnjYUbOfHw9lXa9GiXxqyVWwCYvaqQE/q13+88Ywa0Z/qyAnaVxdfN54b2yWFVbiFrNm6lbG+IFz5YzOkj+1Rp07drNh/MWw3Ah/PXcPoxfSq2T1+4lvKQs3N3GYtW5nPy8MNinUK9+ujTpRQV7wg6jMAo/8abf/esFhTs2MOWHXsod+fT9VsZlNOySpuOLZuzZPN2AJZu3sHAavsbrVCobssBuPte4ArgdeBzwlcFfmZmd5jZNyPNrgF+ZmbzCfdQXVBp/tRoYL27r6x02mbA62a2AJgHbAAeqi3NRldcReS6e3ql5V9BBbJpSxEdsrMq1ttnZ5G/pahKm5//8Gxeenc6J/7oKi699R5uuvQ8AE497ihSmzdj9HlXMuaCq/npt0+ndUZ6TOOvi3Ytm5FfXFqxnl+8i7bVvr0t21jCmP7hb2wn9W9PevMmtEpNqdLmtIEdeX1BXv0HHGUd22SwoaCkYj23YBsd22RUabNoVT5nHhce/jxzVF8yWjQjMyOVRSvzGTP8MFKbNSGrZSrHDepGp+wE+c9XJM61Sk2haGdZxXpRadl+/69t2FrKkE7hf9ODc1qSmpJMWtNkAJokJXHdyb245sTD9ivK4p2Hyuu01Hp+91fdvY+7H+buv49su8Xdp0ReL3b3Y919sLsPcfc3Kh37nrsfU+18O9x9mLsPcvcB7n6lu9caSFwUV2bWxcz+a2abzWyLmd0XpfNebGZzzGzOxKdfjMYp91PTBQVWbc7dq+/P5OxTjuO9x/6XB2+/huvvmUgoFGLh0pUkJyXx/uP/y5uP3MOjL0xlXV78zDk6mJmFf5m6hKE9Mnni8lEM655JfvEuykP7WmWnN6NX+wxmxNmQIIDV8AOo/vtw68Nvc+wRXXn3bxcyamBXcgu2sbc8xHtzV/HW7BW8dvcFPHT9t5j9xQb2JshjI0TiXU3/91X/z++FhXn0yk7n+jG96dU2naKdeyr+77vltc/58zvLmTR7HecMyiE7rRFNCajnm4g2FA1+QnvkpmAvA+8QnqFfDgyv5bB2ZpYP7AReBG5y9/36nytfdRBaPvOAl1V+Xe2zs9hYUFixnl9QSLs2rau0ef6N93nojl8BcOThvdi9p4yibdt5+b2ZHDdsIClNmtCmdUuG9u/NouWr6NKx1pvDNgj523bTvtW+Ycz2rZpTULK7SpuCkt1c+1T4Pm2pTZM5uX8Htu/eW7F/7BEdeHdxPntD9fL21KvcghI6Ze/rqcrJbsnGwu1V2mws3M6Pf/8fANKapzD+2H6U7Az/jO59Zhr3PjMNgH9cdxYrNxQiIg3f1tIyMlvs66nKTE2heFdZlTbFu/by8Mw1ADRNTmJITkt27Q1V7APYsmMPyzZvp3PrVAp27IlR9PUsQb4kxkPP1dFADnBtpHtul7t/dID2XwBDgI7AycAwILBZkwP79GDNhnzWb9zMnrK9vPrBx5w04sgqbXLatmHmvPANZFeszWV3WRlZrTLo2LYNH89fjLuzc9du5n+xgp6dOwaRxteyeEMxXdq0IKd1Kk2SjVMHduD9L6r2vLVukVLRw/OT0T2Z8un6KvtPG9SR1xfG35AgwNylufTMyaJr+1akNEni7NH9eW3m0iptslqmVuR/1bnH8sQb84HwZPjMjHBh2r97OwZ0b8e7n65ERBq+NUU7aZvelDYtUkg2Y2jn1izI3ValTVrT5IoertP6tWXmmvB0kdSUZJokWUWbntlpbNy2K5bhSxQ0+J4rwpdVrolMVKuVu28ENkZWV5nZdcArwCX1FN8BNUlO5qbLzueim+8iFArx7bGj6d2tM399/L8c0bs7Jx8zlOsumsAtf32Ef01+HcO48+qLMDN+cOYYfvOXhxl/+Y3gcPbY4+nbo2sQaXwt5SHnzy8v5r4fDyc5yZj86XpWbtrOpSf3YnFuMR98sZlhPbK4Ymwf3GHu6kL++PK+pxR0bJ1K+1bN+WR1fPbYlIec6//+Os/9bgLJSUk8+cZ8lqwt4IbzRjNvWR5TP17GsQO7cfMFJ+E4Mxat47r7pwKQkpzEK3edD0DJzj1ceveUKsOljcHjd17C6GF9yW6dzsqpd3PHg5OZ9OKHQYcVM8q/8eYfcnh2Xi4/P64nZjBzdREbS3bzjf7tWVtUysK8bfRum843j+gADssLdvDsvPBV5B0ymjFhaCdCHr7i+s0lm9hYrcc/ntXn428aEqvlJqOBi9zkazKQc7AFVrXjRwCvuXvWgdrV17BgPDjqX1uDDiFQqz/5JOgQAlWStyLoEEQCcfFNVwUdQuDuO2dQjVPE6svudx6r02dts5N/FNN4v654GBacBeQBfzSzNDNrbmbHflXjyK0YulpYF+CPhIszERERCVKCTGhv8MVV5JLH8UAvwjf/Wk/4uT9fZSgwA9gBTCd8s6//qecwRUREpBYeCtVpiRfxMOcKd18LfOsg295LgBPYRUREJLHFRXElIiIijUAcDe3VRYMfFqyJmT1Y7fE2Xy4PBh2biIiIfIUEmXMVlz1X7n4pcGnQcYiIiMjBi6d5U3URlz1XIiIiIg1VXPZciYiISByKo6G9ulBxJSIiIrGh4kpEREQkehLl8TcqrkRERCQ2NKFdRERERA6Veq5EREQkNjTnSkRERCR6XMWViIiISPQkyk1EVVyJiIhITHh5YhRXmtAuIiIiEkXquRIREZGYSJSeKxVXEWXrVwQdQmAy23cLOoRAFXbuG3QIEqCSvMT9t5/oendIDzqEhKM5VyIiIiJRlCg9V5pzJSIiIhJF6rkSERGRmEiUnisVVyIiIhITIT24WURERCR6NKFdREREJIoSZVhQE9pFREREokg9VyIiIhIT6rkSERERiSIPheq01MbMxpnZEjNbbmY31LC/q5m9a2ZzzWyBmZ0R2d7dzErNbF5kebDSMcPMbGHknH81M6stDvVciYiISEyE6rHnysySgfuBscB6YLaZTXH3xZWa3QQ86+5/N7P+wKtA98i+Fe4+pIZT/x24GJgZaT8OeO1AsajnSkRERGLCy0N1WmpxNLDc3Ve6+x7gaeCs6iEALSOvWwG5BzqhmXUEWrr7DHd34DHgW7UFouJKRERE4oKZXWxmcyotF1fa3QlYV2l9fWRbZbcB55nZesK9UL+otK9HZLjwfTM7vtI519dyzv1oWFBERERioq4T2t19IjDxK3bXNBfKq61PACa5+z1mNhJ43MyOAPKAru6+xcyGAS+a2YCDPOd+VFyJiIhITNTzTUTXA10qrXdm/2G/CwnPmcLdZ5hZcyDb3TcBuyPbPzGzFUCfyDk713LO/WhYUERERGKinudczQZ6m1kPM2sKfB+YUq3NWmAMgJkdDjQHNptZ28iEeMysJ9AbWOnueUCJmR0TuUrwR8Dk2gJRz5WIiIjEPXffa2ZXAK8DycAj7v6Zmd0BzHH3KcA1wENmdjXh4b0L3N3NbDRwh5ntBcqBS929MHLqy4BJQCrhqwQPeKUgqLgSERGRGKnvm4i6+6uEJ6pX3nZLpdeLgWNrOO4/wH++4pxzgCMOJQ4VVyIiIhITIT24WaLlo0XL+dOzrxMKhfj2cUdy4bjjquzPKyzmpkdfpKR0N+WhEFedPYbjB/bmlY8XMumN6RXtlm7I55nfXEy/Lh1incLXdlTX1lx+fE+SzHhtcT5Pf7q+yv52Gc341cm9aZ2aQsnuvdz5xhIKduwJ70tvxi9P7kXb9GYA3PjSZ+SX7I55DnUxekB7bv3eEJKSjGc+WsWDU5dU2d8pqwV/+vFw2mQ0ZeuOMq7+5yw2bi0F4IZzBnLSwA4kmfHR4nxuf2Z+ECnUycnDenLnJaeSlGT8+/V5/N9zM6rs79yuJX+76kzatGpBUckuLrtrMrlbSgC49ScnMfaoXgDc/fRHvPjB5zGPvz5NvPUnnDF6MJsLt3Hkd2+p/YBGprHnv3bhHD568kFCHqL/8eMY+o1zq+wvKcjnnUf/QmlJMc3TMjjlZ9eSntW2Yv+e0h089ZtL6DF0FKPPuzzW4debRHn8jYqrelYeCvGHp15j4lXn0T6zJRPufJgTB/XlsJx9/4gmvvIhpw4fwPdOGM6K3M38/L4nmTrwSr4xYiDfGDEQCBdWVz7wTFwVVkkGvzjhMK6fvIjN2/dw/7lDmL5qC2uLSivaXHJsD95csok3v9jEkE6tuHBkd/701lIArh/bhyfmrOPTdVtpnpKE13rxa8OSZHDHD47k/L98yMainUy+cQxvzc9leV5JRZsbvzuI/85cw39nrGFk37Zc9+0j+OUjsxnasw3DDmvD6be/CcBz153EiD5t+Xjp5qDSOWRJScafLx/HOb95ktyCbbz1vz9l6sxlLFlXUNHmjgtP4Zm3F/L02ws5fnA3bv7JSVx29xTGHtWLQb06cMIVD9MspQkv/fl83p69gpLSPQFmFF2PvTSNB555m0d/e1HQoQSiMecfCpXzwb/vZ/w1fyA9K5vn77iS7kNGkNWpW0Wb6c8+TN9RY+h37FjWfz6Pmf+ZxCk/u7Zi/8cvPE5O34FBhF+vvLw86BBiolFeLRiZ9f+kmW01syIzeyKoWBat2kDXdpl0bptJSpNkxg0fwLvzq/ZemMGO0nCPzPbSXbRtlbHfeV6btYjTjzqkId/A9W2fQW7xLvK27WZvyHlv2WaO7dmmSptumanMXbcVgHkbihnVMwuArpmpJBt8Gtm3qyzE7r3x9Y1ncI8s1mzazrqCHZSVOy/NXsfYwTlV2vTqmMH0zzcBMGPJZk6J7HecZinJpDRJomlKMk2SjYJtu2KeQ10M7ZPDqtxC1mzcStneEC98sJjTR/ap0qZv12w+mLcagA/nr+H0Y/pUbJ++cC3lIWfn7jIWrczn5OGHxTqFevXRp0spKt4RdBiBacz5b1q5lFbtcmjVriPJTVLoNeIEVs2bWaVNYe5aOh8eftJKp36DWTV3X6/uptXLKN1WRJcBQ2Mat0RPoyyugP8CG4FuQDvg7qACyd9aQvvMVhXr7TNbsmlrSZU2l40/gZc/Xsgp1/+Fy+97il9/f9x+53l9zuK4K66y05qyqdIw3ubtu2mT1rRKm5VbdnD8YdkAHNezDWlNm9CyeRM6t05l+55ybj29Hw9+bwgXj+pOUq2PymxYOrROJa9wXy/dxq2ldMhMrdLm83XFjBsavtnvaUfmkJGaQuu0psxdWcjMJZuZddeZfHzXmXy4OJ8VG6v+3jR0HdtksKFgX8y5Bdvo2KbqF4dFq/I587h+AJw5qi8ZLZqRmZHKopX5jBl+GKnNmpDVMpXjBnWjU3ZLROLBjq0FVYb40jOz2VG0pUqb7C49WfHJNABWfjqdsl2l7Nq+DQ+FmP7MQ4w8t/H16EH9P7i5oYiL4srMupjZf81ss5ltMbP7DtD2VMI3EbvW3Yvdvczd58Yu2tpVrxFem7WIs0YN5q0/Xc0DV0zgxkdfJBTaNwa2YNV6mjdNoXendrENtI4Ophb6x7TVDOrUkge/N4RBnVqxeftuykNOcpIxsGNLJk5bxeXPzqNjq+ac2q99vcccTTU9N7360OYfnl/AiD5tefmmMYzo05a8op2Uh0J0a5vGYR0zGHn9K4y87mVG9m3H0b2zYxN4lNScf9UfwK0Pv82xR3Tl3b9dyKiBXckt2Mbe8hDvzV3FW7NX8NrdF/DQ9d9i9hcb2BtH/7FKYqtpCkP1fw+jzr2I3CULefa2n5O7ZCFpmW2wpGQWvfsy3QYdRUal4qwxqef7XDUYDX7OVeSmXi8D7wDnE77/xPADHHIMsAT4l5mdDqwEfuXu79dw7osJP+ma+375Uy4af3KUo4f2rTPILyquWM8v2kbb1lW/vb8wbR5//58fADD4sC7sLttL0fadtGmZBsDU2Z9x+lEDoh5bfdu8Yw/tMppVrLdNb8aWHVXnzGzZsYfbX/sCgOYpSRx/WBt27CmnYPselhfsIG9buOdr2sotHN4+g6lxNKc5r6iUjln7eqo6tE4lf2tplTabindx2YPh4YAWzZIZN7QTJaV7mXB8T+atLGTn7vD8hPcWbeTInlnMWlZAvMgtKKFT9r7f9Zzslmws3F6lzcbC7fz49+Grn9OapzD+2H6U7Ay/5/c+M417nwl/s//HdWexckMhIvEgPTOb7YX75kduLyqgReuqUyLSMttw+hU3A1C2q5SVn3xEsxZpbFzxOXlLP2PROy9TtnsX5XvLSGnWnJHf/WlMc6gv8VQg1UU89FwdDeQQ7ona4e673P2jA7TvDJwKvAt0AO4BJpvZfl/73X2iuw939+H1UVgBDOjeiTWbCllfUETZ3nKmzvmMEwdXnXfSIaslH3+xCoCVeZvZU7aXrIwWAIRCzhufxN+QIMCS/BI6tUqlQ0YzmiQZJ/Zuy/RVVT8gWzZvUtHDNWFYF6Yuzg8fu6mE9GZNaNU8XP8P6dyaNUVVC5OGbsHqIrq3S6dzmxakJBvjj+rCW/PzqrTJTG9a8Y328tP78dy01QBsKNzJ0X2ySU4ymiQbI/q0rTIRPh7MXZpLz5wsurZvRUqTJM4e3Z/XZi6t0iarZWpF/ledeyxPvBG+IjIpycjMCBem/bu3Y0D3drz76cqYxi/ydbXr0Yfi/Fy2bd5I+d4yln/8Pj2GHFOlTWlJccUw1yevPMPhx50KwNiLr+dHdz/G+Xf9i1HnXkTfUac0msIKIFQeqtMSLxp8zxXhIb417r73INuXAqvd/Z+R9afN7DeEbxpW6y3ro61JchI3fv90Lvu/JygPOd86dgi9ctpx/5R36d8th5MG9+VX3zmV2//9Eo+//TEG/PaCs7DIJ84ny9bQPrMlndtmxjr0Ogs5/O2DFfzxrCNIMpi6OJ81hTv58dFdWbppOzNWFzI4coUgwIINxfzt/RUVx/5j2iru+tZAzGDppu28+tnGALM5dOUh59an5vHYVceTlGQ8N201y/K2cfU3+7NwTRFvzc/jmD5tufbscOE8a2kBtzwVHsF+7ZP1jOrXjqm3jsUd3v9sI28vyDvQX9fglIec6//+Os/9bgLJSUk8+cZ8lqwt4IbzRjNvWR5TP17GsQO7cfMFJ+E4Mxat47r7pwKQkpzEK3edD0DJzj1cevcUykNxdrloLR6/8xJGD+tLdut0Vk69mzsenMykFz8MOqyYacz5JyUnc/x5l/HSvTfhoXL6HXcqWZ26MeuFx2jbvQ89jjyG3CULmPn8JDAjp88Rjep2CwJWfQ5EQxN5avVkIOdgCiwzuxD4jbv3rLRtIXCTu39lcbX7vSca9g+iHn1jUbfaGzViKxbEV9EWbdvWL6m9USNWkrci6BAkIHfdd1PQIQTuymN7xvRSoS8u+ladPmv7PfxiXFzaFA/DgrOAPOCPZpZmZs3NbL9b11fyApBpZj82s2Qz+w7QCZgWi2BFRESkZokyob3BF1fuXg6MB3oRfpr1euB7B2hfCHwT+BVQDNwAnOXu8TMTWEREpBHycq/TEi/iYc4V7r4W+NYhtP8QaHy3thUREZEGLy6KKxEREYl/8XTFX100+GHBmpjZg2a2vYblwaBjExERkZp5yOu0xIu47Lly90uBS4OOQ0RERA5eKI7mTdVFXBZXIiIiEn/i6Yq/uojLYUERERGRhko9VyIiIhIT8XQ7hbpQcSUiIiIxoTlXIiIiIlGkOVciIiIicsjUcyUiIiIxEYqje1XVhYorERERiQlNaBcRERGJokR5/I2KKxEREYmJROm50oR2ERERkShSz5WIiIjERKL0XKm4EhERkZjQnKsEM+rtNkGHEJjC/I1BhxCorPbpQYcQqG+cNj7oECRAvTsk7u//tVf8LugQAnfl3Edi+ve5bsUgIiIiEj2J8vgbTWgXERERiSIVVyIiIhITXh6q01IbMxtnZkvMbLmZ3VDD/q5m9q6ZzTWzBWZ2RmT7WDP7xMwWRv48udIx70XOOS+ytKstDg0LioiISEzU59WCZpYM3A+MBdYDs81sirsvrtTsJuBZd/+7mfUHXgW6AwXAeHfPNbMjgNeBTpWO+6G7zznYWFRciYiISEzU85yro4Hl7r4SwMyeBs4CKhdXDrSMvG4F5AK4+9xKbT4DmptZM3ff/XUC0bCgiIiIxAUzu9jM5lRaLq60uxOwrtL6eqr2PgHcBpxnZusJ91r9ooa/5hxgbrXC6tHIkODNZma1xameKxEREYkJD9XtPlfuPhGY+BW7ayp6qneVTQAmufs9ZjYSeNzMjnD3EICZDQD+BJxa6ZgfuvsGM8sA/gOcDzx2oDjVcyUiIiIxESr3Oi21WA90qbTemciwXyUXAs8CuPsMoDmQDWBmnYEXgB+5+4ovD3D3DZE/S4AnCQ8/HpCKKxEREYkJL/c6LbWYDfQ2sx5m1hT4PjClWpu1wBgAMzuccHG12cxaA68Av3b3aV82NrMmZvZl8ZUCnAksqi0QDQuKiIhITBzM7RS+9rnd95rZFYSv9EsGHnH3z8zsDmCOu08BrgEeMrOrCQ8ZXuDuHjmuF3Czmd0cOeWpwA7g9UhhlQy8BTxUWywqrkRERKRRcPdXCU9Ur7ztlkqvFwPH1nDc74Cveh7SsEONQ8WViIiIxESiPP5GxZWIiIjERH3eRLQhUXElIiIiMRFyFVciIiIiUVOeIMWVbsUgIiIiEkXquRIREZGYSJApVyquREREJDYSZVhQxVUMjOyVza++cTjJBi9+sp5JH66qsr9Dq+bcevZAMtOaUlxaxs3Pz2fTtt0M75HFL0/vV9Gue3YaNz43n/c+3xTrFL620QPac+v3hpCUZDzz0SoenLqkyv5OWS3404+H0yajKVt3lHH1P2excWspADecM5CTBnYgyYyPFudz+zPzg0ihThL5vQeVoeJ9AAAU3ElEQVQ4vH063xnciSSD6asKeXPp5ir7M1ukcN6wzqQ3a8LOPeX8a/Y6tpaWAfDXbw8kt3gXAEU7y/jHjNUxjr7uEj3/tQvn8NGTDxLyEP2PH8fQb5xbZX9JQT7vPPoXSkuKaZ6WwSk/u5b0rLYV+/eU7uCp31xCj6GjGH3e5bEOv15NvPUnnDF6MJsLt3Hkd2+p/YBGQj1XEWa2GrjI3d+q/3Ciy8xOBP7t7p2DiiHJ4Ibx/bl80mzyt+3i8UtH8v4Xm1i1eUdFm6vH9eOVeRt4eV4uR/XI4oqxfbjlPwuZs6qQHzwwHYCWqSm8eNXxzFxeEFQqhyzJ4I4fHMn5f/mQjUU7mXzjGN6an8vyvJKKNjd+dxD/nbmG/85Yw8i+bbnu20fwy0dmM7RnG4Yd1obTb38TgOeuO4kRfdrycbUPp4Yskd97CD9B9dwhnbjvo1Vs3VnGtSf3YmHeNjaW7HvQ/NkDOzJrzVY+XltEn7ZpfHNABx6bE36ofVl5iD++vSyg6Osu0fMPhcr54N/3M/6aP5Celc3zd1xJ9yEjyOrUraLN9Gcfpu+oMfQ7dizrP5/HzP9M4pSfXVux/+MXHien78Agwq93j700jQeeeZtHf3tR0KFIPdCE9no2oHNr1m3ZyYaiUvaWO28s3MiJh7ev0qZHuzRmrdwCwOxVhZzQr/1+5xkzoD3TlxWwq6z+Hh0QbYN7ZLFm03bWFeygrNx5afY6xg7OqdKmV8cMpkd6Y2Ys2cwpkf2O0ywlmZQmSTRNSaZJslGwbVfMc6iLRH7vAbpntaBgxx627NhDuTufrt/KoJyWVdp0bNmcJZu3A7B08w4GVtsfzxI9/00rl9KqXQ6t2nUkuUkKvUacwKp5M6u0KcxdS+fDhwDQqd9gVs2dse/41cso3VZElwFDYxp3rHz06VKKinfU3rCRKXev0xIvGk1xZWYNcoizXctm5BeXVqznF++ibUazKm2WbSxhTP8OAJzUvz3pzZvQKjWlSpvTBnbk9QV59R9wFHVonUpe4b7cN24tpUNmapU2n68rZtzQTgCcdmQOGakptE5rytyVhcxcsplZd53Jx3edyYeL81mxsYR4ksjvPUCr1BSKdpZVrBeVlu2X24atpQzpFC4oBue0JDUlmbSmyQA0SUriupN7cc2Jh+1XlMSDRM9/x9aCKkN86ZnZ7CjaUqVNdpeerPgk/IzclZ9Op2xXKbu2b8NDIaY/8xAjz1WvTmNT7nVb4sVBF1dmlmxmN5rZCjMrMbNPzKyLhf3FzDaZWbGZLTCzIw5wnh5mttXMkiLrD5vZpkr7/21mV0Ve55jZFDMrNLPlZvazSu1uM7PnI+23AReYWaqZTTKzIjNbDBz1dX4o0WQ1bKv++/GXqUsY2iOTJy4fxbDumeQX76I8tK9VdnozerXPYEa8DQvVkHz1Lx5/eH4BI/q05eWbxjCiT1vyinZSHgrRrW0ah3XMYOT1rzDyupcZ2bcdR/fOjk3gUZLI7z3UnH/1H8ALC/PolZ3O9WN606ttOkU791Tkf8trn/Pnd5YzafY6zhmUQ3Za03qPOZoSPf+aOhmq/58w6tyLyF2ykGdv+zm5SxaSltkGS0pm0bsv023QUWRUKs6kcUiUnqtD6e35JTABOANYCgwCdhJ+avRooA9QDPQDtn7VSdx9VaQYOhL4BDge2G5mh7v755Fz3RNp/hTwGZATOe+bZrbS3d+O7D8L+C7wI6AZcCtwWGRJA147UEJmdjFwMUDXM35B9tAzDvZncdDyt+2mfat9vTXtWzWnoNKcC4CCkt1c+9Q8AFKbJnNy/w5s3723Yv/YIzrw7uJ89obi5xcLIK+olI5Z+3Lv0DqV/K2lVdpsKt7FZQ+GhwJaNEtm3NBOlJTuZcLxPZm3spCdu8sBeG/RRo7smcWsZfFTZCTyew+wtbSMzBb7emoyU1Mo3lVWpU3xrr08PHMNAE2TkxiS05Jde0MV+wC27NjDss3b6dw6lYIde2IUfd0lev7pmdlsL9w3R3J7UQEtWrep0iYtsw2nX3EzAGW7Sln5yUc0a5HGxhWfk7f0Mxa98zJlu3dRvreMlGbNGfndn8Y0B4m+eOp9qotDGRa8CLjJ3Zd42Hx33wKUARmEix9z98/dvbYxjPeBE8ysQ2T9+ch6D6AlMN/MugDHAde7+y53nwc8DJxf6Twz3P1Fdw+5eylwLvB7dy9093XAXw8UhLtPdPfh7j68PgorgMUbiunSpgU5rVNpkmycOrAD739R9Yqv1i1SKr7R/WR0T6Z8ur7K/tMGdeT1hfE3LLRgdRHd26XTuU0LUpKN8Ud14a35VfPITG9akfvlp/fjuWmrAdhQuJOj+2STnGQ0STZG9GlbZSJ8PEjk9x5gTdFO2qY3pU2LFJLNGNq5NQtyt1Vpk9Y0uaKH57R+bZm5pgiA1JRkmiRZRZue2WlsjLM5d4mef7sefSjOz2Xb5o2U7y1j+cfv02PIMVXalJYU46FwMfnJK89w+HGnAjD24uv50d2Pcf5d/2LUuRfRd9QpKqwkrhxKz1UXYEX1je7+jpndB9wPdDWzF4Bfufu26m0reR/4JrAe+AB4j3DRtAv40N1DZpYDFLp75U/UNcDwSuvrqp03p9q2NQeTWH0qDzl/fnkx9/14OMlJxuRP17Ny03YuPbkXi3OL+eCLzQyLXCXmDnNXF/LHlxdXHN+xdSrtWzXnk9WFAWbx9ZSHnFufmsdjVx1PUpLx3LTVLMvbxtXf7M/CNUW8NT+PY/q05dqzw6PIs5YWcMtTcwF47ZP1jOrXjqm3jsUd3v9sI2/H2byjRH7vAUIOz87L5efH9cQMZq4uYmPJbr7Rvz1ri0pZmLeN3m3T+eYRHcBhecEOnp23AYAOGc2YMLQTIQ9fdfnmkk1VrrKLB4mef1JyMsefdxkv3XsTHiqn33GnktWpG7NeeIy23fvQ48hjyF2ygJnPTwIzcvoc0ehut3Agj995CaOH9SW7dTorp97NHQ9OZtKLHwYdVr1LlJ4r81rGML+8FQPh4uk6d598gLbtgGcJF0g3H6Bdb2AW4WG/6cAUYBHwNvCZu98d6blaDbT+ssAysz8AOe5+gZndBvRy9/MqnXcVcJm7T42s/wy49WBuxTDs5qkJ8pbvrzB/e9AhBCqrfXrQIQRq5JCc2htJo9W7Q+L+/l97xe+CDiFwe+Y+UuP0wPpyX6u+dfqsvaJ4SUzj/boOZVjwYeC3ZtY7Mol9kJm1MbOjzGyEmaUAOwj3PpUf6ETuvgwoBc4DPoj0cuUD5xDu1SIyrDcduNPMmpvZIOBC4IkDnPpZ4NdmlmlmnYFfHEJ+IiIiUo90teD+7iVcvLwBbAP+CaQSniP1EFBEeBhuC3D3QZzvfWCLu6+ttG7A3EptJgDdgVzgBcK9UG8e4Jy3R2JYFYnz8YOIQ0RERGJAVwtGuHv3Squ/iyyVrSd85eAhcfcJ1dZ/Bfyq2rb1wJlfcfxtNWzbSfjKwcruOtTYRERERL6uBnnjTREREWl84mlory7qrbgys8+AbjXsusTdDzRvSkRERBqheBraq4t6K67cfUB9nVtERETiT6L0XDWaZwuKiIiINASacyUiIiIxoWFBERERkSgKBR1AjKi4EhERkZhQz5WIiIhIFGlCu4iIiIgcMvVciYiISExoWFBEREQkihJlWFDFlYiIiMREovRcac6ViIiIxES5122pjZmNM7MlZrbczG6oYX9XM3vXzOaa2QIzO6PSvl9HjltiZqcd7DlrouJKRERE4p6ZJQP3A6cD/YEJZta/WrObgGfd/Ujg+8ADkWP7R9YHAOOAB8ws+SDPuR8NC4qIiEhM1POw4NHAcndfCWBmTwNnAYsrtXGgZeR1KyA38vos4Gl33w2sMrPlkfNxEOfcj4orERERiYl6ntDeCVhXaX09MKJam9uAN8zsF0AacEqlY2dWO7ZT5HVt59yPiquIT347zoL8+83sYnefGGQMQVL+iZt/IucOyj/I/K+c+0gQf20Vifb+P+ir6/RZa2YXAxdX2jSx0s+vpnNXL+cmAJPc/R4zGwk8bmZHHODYmqZP1Voias5Vw3Fx7U0aNeWfuBI5d1D+yl8OmrtPdPfhlZbKhel6oEul9c7sG/b70oXAs5FzzQCaA9kHOPZgzrkfFVciIiLSGMwGeptZDzNrSniC+pRqbdYCYwDM7HDCxdXmSLvvm1kzM+sB9AZmHeQ596NhQREREYl77r7XzK4AXgeSgUfc/TMzuwOY4+5TgGuAh8zsasLDexe4uwOfmdmzhCeq7wV+7u7lADWds7ZYzBPkhl4NXaKNu1en/BM3/0TOHZS/8k/s/BsrFVciIiIiUaQ5VyIiIiJRpOJKREREJIpUXImIiIhEkYorkQbAzJLMrHXQcQTFzL5vZs2CjkNiy8yyzGx8ZMkKOh6RaFFx1cBE7rFRHnQc9cnM0s3sATObb2ZPmtlh1fYvDCq2WDCzHDObbGaLzOx8MxtB+KZ0W8xstpl1qe0c8crMeta0AP8HDIq8brTM7PxKr1uY2d/NbLOZFZjZw2aWFmR89c3MHjezwZHXo4EVwO8jy1IzGxVkfCLRoqsFA2BmXQ+wuznwhbs32sLXzP5B+I63DwEnAOcD57j7e5H9Je6eEVyE9cvM/gMsB0LALwk/pf0xIAW4F9jr7j8ILsL6Y2YhwveW+apHYLi7J8cwpJgys23u3jLy+l5gGHBdZPedwAJ3vyqo+OqbmW0B2rl7uZl9DNzr7s9E9n0HuM7djz7gSeKUma3jIB6b4u4H+nyQOKHiKgD6gLFc4HB3L46sjwWeIHwzt1cToLjaDOQQ7jneCbRx962RfW2B+e6eE2CI9cbMHgV6Er5B36JK2/OAIe6eH1hwMVD5d9vMVgAnu/uayHoX4CN37xZkjPXJzIqALu6+3cwKgPaVbtSYBGz9svhsbMzshEqrRwE/Bv4KrAG6AVcAj7n7PQGEJ1Gm4ioAZrYB+DYwp4bdzYFtjby4KgI6uvuuStuOAV4EfgH8s7H+Bwvh/N09s/rryLoRfv8bc3F5LPA34D3glsgHbR4w2N03BRpcPavWc5Xn7h2r7d/u7unBRFf/zOwJwkPg1xHurX3b3SdF9p0PXOvug4KLMDbMbBFwmrtvqLStMzDV3Y8ILjKJlkY79NTAzQGOdPfy6gvh2+7X6anhcWARMLryBnefCZxO+EO3RRBBxdA6M+sQeV39g6Q34edcNVruPo3wN/c84FMzOy/gkGIpzczWmtlaINPM+n25I9JztTW40GLiSmAE4blWbYF/mtkKM1sO/C9wWZDBxVAOsL3atu1ApwBikXqgZwsG4xLC82324+67afxF753AflcGuftcMzuF8Dykxux/gHIAd19Xbd8AwvOuGrXIF4m7zOxp4C+Ee2xr/DfRyJxcbb2w0uueQKMeEnL3AmB0ZCrAMYSHxEqBBcAL7l4UZHwxNAWYYma/A9YTnoP6aw7igcASHzQs2MCZ2QPufnnQcQRF+Sdu/omcOyj/xpy/mTUHbgO+S7gXKxd4Drjd3UsDDE2iRMVVA1d5jkYiUv6Jm38i5w7Kv7Hmb2bJhCezP1l53qk0Lo19+KkxaOzzr2qj/BNXIucOyr9R5h8ZEr9XhVXjpuKq4Uv0rkXln7gSOXdQ/o05/5fMbHzQQUj90YR2ERGR2GoOPG9mM4AqNxd19x8FFpVEjYqrhq9Rdo0fAuWfuBI5d1D+jTn/RZFFGikNCwbIzP5a/VlaZjbKzP630qZ/xzismFH+iZt/IucOyj/R83f3279qCTo2iQ5dLRigyGNQOrn7nkrbmgHr3L1dcJHFhvJP3PwTOXdQ/omeP4CZnUT4uaqdgA3Av939nWCjkmhRz1WwnP3fg+QatjVWyj9x80/k3EH5J3T+ZnYR8AywEfgv4acVPGlmPws0MImahPhFbsA+BH4XeWDplw8uvS2yPREo/8TNP5FzB+Wf6PlfB4x19xvd/R/u/hvg1Mh2aQQ0LBigyIM6XwY6En4MRFfC32DGu/v6IGOLBeWfuPkncu6g/JW/bQE6uHtZpW3NgFx3bxNcZBItKq4CFvnGdjThZ0utA2a5eyI8Yw1Q/omcfyLnDso/kfM3s8nAWuB6d99pZmmEn7naw911/6tGQMWViIhIDJlZR+BpYBThh3dnAdOBCe6eG2RsEh0qrkRERAJgZl0ID43mJsJwaCJRcSUiIhJDZvY/wHvuviDoWKR+6GpBERGR2DqK8PMFC81sspn90syGfXn1pMQ/9VyJiIgEwMy6AydElm8DuHvrAEOSKNGzBUVERGLMzPoSLqpOBI4FlgLvBxmTRI96rkRERGLIzPKBEuB54D1gmruXBBqURJV6rkRERGLrJeA44FtAJpBlZh/oisHGQz1XIiIiATCz9sBowsOD5wEF7t4r2KgkGtRzJSIiEmNmdiThouok4HhgBzAr0KAkatRzJSIiEkNmVgQUAx8QnsT+vrsvDzYqiSYVVyIiIjFkZt3dfXUtbSa4+1MxCkmiTMWViIhIA2Nm29y9ZdBxyNeju8GKiIg0PBZ0APL1qbgSERFpeDSsFMdUXImIiIhEkYorERERkShScSUiItLwrA06APn6dLWgiIhIDJlZf2CLu+ebWTpwLVAO3O3uO4ONTqJBPVciIiKx9STQOvL6bsKPwBkJ/COwiCSq1HMlIiISQ2a21d1bm5kBG4EBQCmwyt3bBRudRIOeLSgiIhJbu80sA+gPrHP3AjNrAjQPOC6JEhVXIiIisfUk8C6QDtwX2TYUWBVYRBJVGhYUERGJMTM7FShz93cj68OBlu7+TrCRSTSouBIREYkhM+sE7HT3okrbsoDm7p4bXGQSLbpaUEREJLZeBDpX29YJeCGAWKQeqOdKREQkhsys2N1bHex2iT/quRIREYmtzWbWq/KGyPqWgOKRKFNxJSIiEluPAP8xszPNrL+ZjQeeBx4OOC6JEg0LioiIxJCZJQHXABcCXQg/R/CfwD2uD+VGQcWViIhIPTOzk79qF1DxQaxbMTQOKq5ERETqmZl91Q1Cv/wQNsDdvWeMQpJ6pOJKREREJIo0oV1EREQkilRciYiIiESRiisRERGRKFJxJSIiIhJFKq5EREREouj/AYE9fxhSTc8/AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 720x432 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "plt.figure(figsize=(10,6))\n",
    "s = sns.heatmap(corr_matrix, annot = True, cmap = 'RdBu', vmin = 0.8, vmax=1)\n",
    "s.set_yticklabels(s.get_yticklabels(), rotation = 0, fontsize =12) \n",
    "s.set_xticklabels(s.get_xticklabels(), rotation = 90, fontsize =12) \n",
    "plt.title('Correlation Heatmap', fontsize =18)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Answer:**\n",
    "<ol>\n",
    "    <li>C1</li>\n",
    "    <li>C5</li>\n",
    "    <li> LCS method</li>\n",
    "</ol>\n",
    "c_3 and c_4,  c_4 and c_5 , c_5 and c_6 and are the highest correlated features.\n",
    "We drop C2 as it has correlation above the threshold (0.97) and select C1 and C5 that have low correlation with each other correlation with each other and lcs_word as it's different measurement of similarity."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Creating Final Data Files\n",
    "\n",
    "Now, you are almost ready to move on to training a model in SageMaker!\n",
    "\n",
    "You'll want to access your train and test data in SageMaker and upload it to S3. In this project, SageMaker will expect the following format for your train/test data:\n",
    "* Training and test data should be saved in one `.csv` file each, ex `train.csv` and `test.csv`\n",
    "* These files should have class  labels in the first column and features in the rest of the columns\n",
    "\n",
    "This format follows the practice, outlined in the [SageMaker documentation](https://docs.aws.amazon.com/sagemaker/latest/dg/cdf-training.html), which reads: \"Amazon SageMaker requires that a CSV file doesn't have a header record and that the target variable [class label] is in the first column.\"\n",
    "\n",
    "## EXERCISE: Create csv files\n",
    "\n",
    "Define a function that takes in x (features) and y (labels) and saves them to one `.csv` file at the path `data_dir/filename`.\n",
    "\n",
    "It may be useful to use pandas to merge your features and labels into one DataFrame and then convert that into a csv file. You can make sure to get rid of any incomplete rows, in a DataFrame, by using `dropna`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "def make_csv(x, y, filename, data_dir):\n",
    "    '''Merges features and labels and converts them into one csv file with labels in the first column.\n",
    "       :param x: Data features\n",
    "       :param y: Data labels\n",
    "       :param file_name: Name of csv file, ex. 'train.csv'\n",
    "       :param data_dir: The directory where files will be saved\n",
    "       '''\n",
    "    # make data dir, if it does not exist\n",
    "    if not os.path.exists(data_dir):\n",
    "        os.makedirs(data_dir)\n",
    "    df = pd.concat([pd.DataFrame(y), pd.DataFrame(x)], axis=1).dropna()\n",
    "    df.to_csv(os.path.join(data_dir, filename), index=False, header=False)\n",
    "    # nothing is returned, but a print statement indicates that the function has run\n",
    "    print('Path created: '+str(data_dir)+'/'+str(filename))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Test cells\n",
    "\n",
    "Test that your code produces the correct format for a `.csv` file, given some text features and labels."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Path created: test_csv/to_delete.csv\n",
      "Tests passed!\n"
     ]
    }
   ],
   "source": [
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "fake_x = [ [0.39814815, 0.0001, 0.19178082], \n",
    "           [0.86936937, 0.44954128, 0.84649123], \n",
    "           [0.44086022, 0., 0.22395833] ]\n",
    "\n",
    "fake_y = [0, 1, 1]\n",
    "\n",
    "make_csv(fake_x, fake_y, filename='to_delete.csv', data_dir='test_csv')\n",
    "\n",
    "# read in and test dimensions\n",
    "fake_df = pd.read_csv('test_csv/to_delete.csv', header=None)\n",
    "\n",
    "# check shape\n",
    "assert fake_df.shape==(3, 4), \\\n",
    "      'The file should have as many rows as data_points and as many columns as features+1 (for indices).'\n",
    "# check that first column = labels\n",
    "assert np.all(fake_df.iloc[:,0].values==fake_y), 'First column is not equal to the labels, fake_y.'\n",
    "print('Tests passed!')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "# delete the test csv file, generated above\n",
    "! rm -rf test_csv"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If you've passed the tests above, run the following cell to create `train.csv` and `test.csv` files in a directory that you specify! This will save the data in a local directory. Remember the name of this directory because you will reference it again when uploading this data to S3."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Path created: plagiarism_data/train.csv\n",
      "Path created: plagiarism_data/test.csv\n"
     ]
    }
   ],
   "source": [
    "# can change directory, if you want\n",
    "data_dir = 'plagiarism_data'\n",
    "\n",
    "\"\"\"\n",
    "DON'T MODIFY ANYTHING IN THIS CELL THAT IS BELOW THIS LINE\n",
    "\"\"\"\n",
    "\n",
    "make_csv(train_x, train_y, filename='train.csv', data_dir=data_dir)\n",
    "make_csv(test_x, test_y, filename='test.csv', data_dir=data_dir)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Up Next\n",
    "\n",
    "Now that you've done some feature engineering and created some training and test data, you are ready to train and deploy a plagiarism classification model. The next notebook will utilize SageMaker resources to train and test a model that you design."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "conda_amazonei_mxnet_p36",
   "language": "python",
   "name": "conda_amazonei_mxnet_p36"
  },
  "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.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}