t81_558_class4_training.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# T81-558: Applications of Deep Neural Networks\n",
    "**Class 4: Training a Neural Network**\n",
    "* Instructor: [Jeff Heaton](https://sites.wustl.edu/jeffheaton/), School of Engineering and Applied Science, [Washington University in St. Louis](https://engineering.wustl.edu/Programs/Pages/default.aspx)\n",
    "* For more information visit the [class website](https://sites.wustl.edu/jeffheaton/t81-558/)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Module Video Material\n",
    "\n",
    "Main video lecture:\n",
    "\n",
    "* [Part 4.1: Early Stopping and Feature Vector Encoding](https://www.youtube.com/watch?v=ATuyK_HWZgc&index=12&list=PLjy4p-07OYzulelvJ5KVaT2pDlxivl_BN)\n",
    "* [Part 4.2: Evaluating Classification and Regression Networks](https://www.youtube.com/watch?v=hXkZqGi5mB4&index=13&list=PLjy4p-07OYzulelvJ5KVaT2pDlxivl_BN)\n",
    "* [Part 4.3: Cross-Validation for Neural Networks](https://www.youtube.com/watch?v=SIyMm5DFwQ8&list=PLjy4p-07OYzulelvJ5KVaT2pDlxivl_BN)\n",
    "\n",
    "Weekly video update:\n",
    "\n",
    "* [Weekly Update #4, Feb 4, 2018](https://www.youtube.com/watch?v=rae7VSDZb6I)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Helpful Functions\n",
    "\n",
    "You will see these at the top of every module.  These are simply a set of reusable functions that we will make use of.  Each of them will be explained as the semester progresses.  They are explained in greater detail as the course progresses.  Class 4 contains a complete overview of these functions."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from sklearn import preprocessing\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import shutil\n",
    "import os\n",
    "import requests\n",
    "import base64\n",
    "\n",
    "\n",
    "# Encode text values to dummy variables(i.e. [1,0,0],[0,1,0],[0,0,1] for red,green,blue)\n",
    "def encode_text_dummy(df, name):\n",
    "    dummies = pd.get_dummies(df[name])\n",
    "    for x in dummies.columns:\n",
    "        dummy_name = \"{}-{}\".format(name, x)\n",
    "        df[dummy_name] = dummies[x]\n",
    "    df.drop(name, axis=1, inplace=True)\n",
    "\n",
    "\n",
    "# Encode text values to a single dummy variable.  The new columns (which do not replace the old) will have a 1\n",
    "# at every location where the original column (name) matches each of the target_values.  One column is added for\n",
    "# each target value.\n",
    "def encode_text_single_dummy(df, name, target_values):\n",
    "    for tv in target_values:\n",
    "        l = list(df[name].astype(str))\n",
    "        l = [1 if str(x) == str(tv) else 0 for x in l]\n",
    "        name2 = \"{}-{}\".format(name, tv)\n",
    "        df[name2] = l\n",
    "\n",
    "\n",
    "# Encode text values to indexes(i.e. [1],[2],[3] for red,green,blue).\n",
    "def encode_text_index(df, name):\n",
    "    le = preprocessing.LabelEncoder()\n",
    "    df[name] = le.fit_transform(df[name])\n",
    "    return le.classes_\n",
    "\n",
    "\n",
    "# Encode a numeric column as zscores\n",
    "def encode_numeric_zscore(df, name, mean=None, sd=None):\n",
    "    if mean is None:\n",
    "        mean = df[name].mean()\n",
    "\n",
    "    if sd is None:\n",
    "        sd = df[name].std()\n",
    "\n",
    "    df[name] = (df[name] - mean) / sd\n",
    "\n",
    "\n",
    "# Convert all missing values in the specified column to the median\n",
    "def missing_median(df, name):\n",
    "    med = df[name].median()\n",
    "    df[name] = df[name].fillna(med)\n",
    "\n",
    "\n",
    "# Convert all missing values in the specified column to the default\n",
    "def missing_default(df, name, default_value):\n",
    "    df[name] = df[name].fillna(default_value)\n",
    "\n",
    "\n",
    "# Convert a Pandas dataframe to the x,y inputs that TensorFlow needs\n",
    "def to_xy(df, target):\n",
    "    result = []\n",
    "    for x in df.columns:\n",
    "        if x != target:\n",
    "            result.append(x)\n",
    "    # find out the type of the target column.  Is it really this hard? :(\n",
    "    target_type = df[target].dtypes\n",
    "    target_type = target_type[0] if hasattr(target_type, '__iter__') else target_type\n",
    "    # Encode to int for classification, float otherwise. TensorFlow likes 32 bits.\n",
    "    if target_type in (np.int64, np.int32):\n",
    "        # Classification\n",
    "        dummies = pd.get_dummies(df[target])\n",
    "        return df.as_matrix(result).astype(np.float32), dummies.as_matrix().astype(np.float32)\n",
    "    else:\n",
    "        # Regression\n",
    "        return df.as_matrix(result).astype(np.float32), df.as_matrix([target]).astype(np.float32)\n",
    "\n",
    "# Nicely formatted time string\n",
    "def hms_string(sec_elapsed):\n",
    "    h = int(sec_elapsed / (60 * 60))\n",
    "    m = int((sec_elapsed % (60 * 60)) / 60)\n",
    "    s = sec_elapsed % 60\n",
    "    return \"{}:{:>02}:{:>05.2f}\".format(h, m, s)\n",
    "\n",
    "\n",
    "# Regression chart.\n",
    "def chart_regression(pred,y,sort=True):\n",
    "    t = pd.DataFrame({'pred' : pred, 'y' : y.flatten()})\n",
    "    if sort:\n",
    "        t.sort_values(by=['y'],inplace=True)\n",
    "    a = plt.plot(t['y'].tolist(),label='expected')\n",
    "    b = plt.plot(t['pred'].tolist(),label='prediction')\n",
    "    plt.ylabel('output')\n",
    "    plt.legend()\n",
    "    plt.show()\n",
    "\n",
    "# Remove all rows where the specified column is +/- sd standard deviations\n",
    "def remove_outliers(df, name, sd):\n",
    "    drop_rows = df.index[(np.abs(df[name] - df[name].mean()) >= (sd * df[name].std()))]\n",
    "    df.drop(drop_rows, axis=0, inplace=True)\n",
    "\n",
    "\n",
    "# Encode a column to a range between normalized_low and normalized_high.\n",
    "def encode_numeric_range(df, name, normalized_low=-1, normalized_high=1,\n",
    "                         data_low=None, data_high=None):\n",
    "    if data_low is None:\n",
    "        data_low = min(df[name])\n",
    "        data_high = max(df[name])\n",
    "\n",
    "    df[name] = ((df[name] - data_low) / (data_high - data_low)) \\\n",
    "               * (normalized_high - normalized_low) + normalized_low\n",
    "        \n",
    "# This function submits an assignment.  You can submit an assignment as much as you like, only the final\n",
    "# submission counts.  The paramaters are as follows:\n",
    "# data - Pandas dataframe output.\n",
    "# key - Your student key that was emailed to you.\n",
    "# no - The assignment class number, should be 1 through 1.\n",
    "# source_file - The full path to your Python or IPYNB file.  This must have \"_class1\" as part of its name.  \n",
    "# .             The number must match your assignment number.  For example \"_class2\" for class assignment #2.\n",
    "def submit(data,key,no,source_file=None):\n",
    "    if source_file is None and '__file__' not in globals(): raise Exception('Must specify a filename when a Jupyter notebook.')\n",
    "    if source_file is None: source_file = __file__\n",
    "    suffix = '_class{}'.format(no)\n",
    "    if suffix not in source_file: raise Exception('{} must be part of the filename.'.format(suffix))\n",
    "    with open(source_file, \"rb\") as image_file:\n",
    "        encoded_python = base64.b64encode(image_file.read()).decode('ascii')\n",
    "    ext = os.path.splitext(source_file)[-1].lower()\n",
    "    if ext not in ['.ipynb','.py']: raise Exception(\"Source file is {} must be .py or .ipynb\".format(ext))\n",
    "    r = requests.post(\"https://api.heatonresearch.com/assignment-submit\",\n",
    "        headers={'x-api-key':key}, json={'csv':base64.b64encode(data.to_csv(index=False).encode('ascii')).decode(\"ascii\"),\n",
    "        'assignment': no, 'ext':ext, 'py':encoded_python})\n",
    "    if r.status_code == 200:\n",
    "        print(\"Success: {}\".format(r.text))\n",
    "    else: print(\"Failure: {}\".format(r.text))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Building the Feature Vector\n",
    "\n",
    "Neural networks require their input to be a fixed number of columns.  This is very similar to spreadsheet data.  This input must be completely numeric.  \n",
    "\n",
    "It is important to represent the data in a way that the neural network can train from it.  In class 6, we will see even more ways to preprocess data.  For now, we will look at several of the most basic ways to transform data for a neural network.\n",
    "\n",
    "Before we look at specific ways to preprocess data, it is important to consider four basic types of data, as defined by [Stanley Smith Stevens](https://en.wikipedia.org/wiki/Stanley_Smith_Stevens).  These are commonly referred to as the [levels of measure](https://en.wikipedia.org/wiki/Level_of_measurement):\n",
    "\n",
    "* Character Data (strings)\n",
    "    * **Nominal** - Individual discrete items, no order. For example: color, zip code, shape.\n",
    "    * **Ordinal** - Individual discrete items that can be ordered.  For example: grade level, job title, Starbucks(tm) coffee size (tall, vente, grande) \n",
    "* Numeric Data\n",
    "    * **Interval** - Numeric values, no defined start.  For example, temperature.  You would never say \"yesterday was twice as hot as today\".\n",
    "    * **Ratio** - Numeric values, clearly defined start.  For example, speed.  You would say that \"The first car is going twice as fast as the second.\"\n",
    "\n",
    "The following code contains several useful functions to encode the feature vector for various types of data.  Encoding data:\n",
    "\n",
    "* **encode_text_dummy** - Encode text fields, such as the iris species as a single field for each class.  Three classes would become \"0,0,1\" \"0,1,0\" and \"1,0,0\". Encode non-target predictors this way. Good for nominal.\n",
    "* **encode_text_index** - Encode text fields, such as the iris species as a single numeric field as \"0\" \"1\" and \"2\".  Encode the target field for a classification this way.  Good for nominal.\n",
    "* **encode_numeric_zscore** - Encode numeric values as a z-score.  Neural networks deal well with \"centered\" fields, zscore is usually a good starting point for interval/ratio.\n",
    "\n",
    "*Ordinal values can be encoded as dummy or index.  Later we will see a more advanced means of encoding*\n",
    "\n",
    "Dealing with missing data:\n",
    "\n",
    "* **missing_median** - Fill all missing values with the median value.\n",
    "\n",
    "Creating the final feature vector:\n",
    "\n",
    "* **to_xy** - Once all fields are numeric, this function can provide the x and y matrixes that are used to fit the neural network.\n",
    "\n",
    "Other utility functions:\n",
    "\n",
    "* **hms_string** - Print out an elapsed time string.\n",
    "* **chart_regression** - Display a chart to show how well a regression performs.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Dealing with Outliers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Length before MPG outliers dropped: 398\n",
      "Length after MPG outliers dropped: 388\n",
      "      mpg  cylinders  displacement  horsepower  weight  acceleration  year  \\\n",
      "0    18.0          8         307.0       130.0    3504          12.0    70   \n",
      "1    15.0          8         350.0       165.0    3693          11.5    70   \n",
      "2    18.0          8         318.0       150.0    3436          11.0    70   \n",
      "3    16.0          8         304.0       150.0    3433          12.0    70   \n",
      "4    17.0          8         302.0       140.0    3449          10.5    70   \n",
      "5    15.0          8         429.0       198.0    4341          10.0    70   \n",
      "6    14.0          8         454.0       220.0    4354           9.0    70   \n",
      "7    14.0          8         440.0       215.0    4312           8.5    70   \n",
      "8    14.0          8         455.0       225.0    4425          10.0    70   \n",
      "9    15.0          8         390.0       190.0    3850           8.5    70   \n",
      "10   15.0          8         383.0       170.0    3563          10.0    70   \n",
      "11   14.0          8         340.0       160.0    3609           8.0    70   \n",
      "12   15.0          8         400.0       150.0    3761           9.5    70   \n",
      "13   14.0          8         455.0       225.0    3086          10.0    70   \n",
      "14   24.0          4         113.0        95.0    2372          15.0    70   \n",
      "15   22.0          6         198.0        95.0    2833          15.5    70   \n",
      "16   18.0          6         199.0        97.0    2774          15.5    70   \n",
      "17   21.0          6         200.0        85.0    2587          16.0    70   \n",
      "18   27.0          4          97.0        88.0    2130          14.5    70   \n",
      "19   26.0          4          97.0        46.0    1835          20.5    70   \n",
      "20   25.0          4         110.0        87.0    2672          17.5    70   \n",
      "21   24.0          4         107.0        90.0    2430          14.5    70   \n",
      "22   25.0          4         104.0        95.0    2375          17.5    70   \n",
      "23   26.0          4         121.0       113.0    2234          12.5    70   \n",
      "24   21.0          6         199.0        90.0    2648          15.0    70   \n",
      "25   10.0          8         360.0       215.0    4615          14.0    70   \n",
      "26   10.0          8         307.0       200.0    4376          15.0    70   \n",
      "27   11.0          8         318.0       210.0    4382          13.5    70   \n",
      "28    9.0          8         304.0       193.0    4732          18.5    70   \n",
      "29   27.0          4          97.0        88.0    2130          14.5    71   \n",
      "..    ...        ...           ...         ...     ...           ...   ...   \n",
      "367  28.0          4         112.0        88.0    2605          19.6    82   \n",
      "368  27.0          4         112.0        88.0    2640          18.6    82   \n",
      "369  34.0          4         112.0        88.0    2395          18.0    82   \n",
      "370  31.0          4         112.0        85.0    2575          16.2    82   \n",
      "371  29.0          4         135.0        84.0    2525          16.0    82   \n",
      "372  27.0          4         151.0        90.0    2735          18.0    82   \n",
      "373  24.0          4         140.0        92.0    2865          16.4    82   \n",
      "374  23.0          4         151.0        93.5    3035          20.5    82   \n",
      "375  36.0          4         105.0        74.0    1980          15.3    82   \n",
      "376  37.0          4          91.0        68.0    2025          18.2    82   \n",
      "377  31.0          4          91.0        68.0    1970          17.6    82   \n",
      "378  38.0          4         105.0        63.0    2125          14.7    82   \n",
      "379  36.0          4          98.0        70.0    2125          17.3    82   \n",
      "380  36.0          4         120.0        88.0    2160          14.5    82   \n",
      "381  36.0          4         107.0        75.0    2205          14.5    82   \n",
      "382  34.0          4         108.0        70.0    2245          16.9    82   \n",
      "383  38.0          4          91.0        67.0    1965          15.0    82   \n",
      "384  32.0          4          91.0        67.0    1965          15.7    82   \n",
      "385  38.0          4          91.0        67.0    1995          16.2    82   \n",
      "386  25.0          6         181.0       110.0    2945          16.4    82   \n",
      "387  38.0          6         262.0        85.0    3015          17.0    82   \n",
      "388  26.0          4         156.0        92.0    2585          14.5    82   \n",
      "389  22.0          6         232.0       112.0    2835          14.7    82   \n",
      "390  32.0          4         144.0        96.0    2665          13.9    82   \n",
      "391  36.0          4         135.0        84.0    2370          13.0    82   \n",
      "392  27.0          4         151.0        90.0    2950          17.3    82   \n",
      "393  27.0          4         140.0        86.0    2790          15.6    82   \n",
      "395  32.0          4         135.0        84.0    2295          11.6    82   \n",
      "396  28.0          4         120.0        79.0    2625          18.6    82   \n",
      "397  31.0          4         119.0        82.0    2720          19.4    82   \n",
      "\n",
      "     origin  \n",
      "0         1  \n",
      "1         1  \n",
      "2         1  \n",
      "3         1  \n",
      "4         1  \n",
      "5         1  \n",
      "6         1  \n",
      "7         1  \n",
      "8         1  \n",
      "9         1  \n",
      "10        1  \n",
      "11        1  \n",
      "12        1  \n",
      "13        1  \n",
      "14        3  \n",
      "15        1  \n",
      "16        1  \n",
      "17        1  \n",
      "18        3  \n",
      "19        2  \n",
      "20        2  \n",
      "21        2  \n",
      "22        2  \n",
      "23        2  \n",
      "24        1  \n",
      "25        1  \n",
      "26        1  \n",
      "27        1  \n",
      "28        1  \n",
      "29        3  \n",
      "..      ...  \n",
      "367       1  \n",
      "368       1  \n",
      "369       1  \n",
      "370       1  \n",
      "371       1  \n",
      "372       1  \n",
      "373       1  \n",
      "374       1  \n",
      "375       2  \n",
      "376       3  \n",
      "377       3  \n",
      "378       1  \n",
      "379       1  \n",
      "380       3  \n",
      "381       3  \n",
      "382       3  \n",
      "383       3  \n",
      "384       3  \n",
      "385       3  \n",
      "386       1  \n",
      "387       1  \n",
      "388       1  \n",
      "389       1  \n",
      "390       3  \n",
      "391       1  \n",
      "392       1  \n",
      "393       1  \n",
      "395       1  \n",
      "396       1  \n",
      "397       1  \n",
      "\n",
      "[388 rows x 8 columns]\n"
     ]
    }
   ],
   "source": [
    "import tensorflow.contrib.learn as skflow\n",
    "import pandas as pd\n",
    "import os\n",
    "import numpy as np\n",
    "from sklearn import metrics\n",
    "from scipy.stats import zscore\n",
    "\n",
    "path = \"./data/\"\n",
    "\n",
    "filename_read = os.path.join(path,\"auto-mpg.csv\")\n",
    "df = pd.read_csv(filename_read,na_values=['NA','?'])\n",
    "\n",
    "# create feature vector\n",
    "missing_median(df, 'horsepower')\n",
    "df.drop('name',1,inplace=True)\n",
    "#encode_numeric_binary(df,'mpg',20)\n",
    "#df['origin'] = df['origin'].astype(str)\n",
    "#encode_text_tfidf(df, 'origin')\n",
    "\n",
    "# Drop outliers in horsepower\n",
    "print(\"Length before MPG outliers dropped: {}\".format(len(df)))\n",
    "remove_outliers(df,'mpg',2)\n",
    "print(\"Length after MPG outliers dropped: {}\".format(len(df)))\n",
    "\n",
    "print(df)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Other Examples: Dealing with Addresses\n",
    "\n",
    "Addresses can be difficult to encode into a neural network.  There are many different approaches, and you must consider how you can transform the address into something more meaningful.  Map coordinates can be a good approach.  [Latitude and longitude](https://en.wikipedia.org/wiki/Geographic_coordinate_system) can be a useful encoding.  Thanks to the power of the Internet, it is relatively easy to transform an address into its latitude and longitude values.  The following code determines the coordinates of [Washington University](https://wustl.edu/):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'lat': 38.6471178, 'lng': -90.3026148}\n"
     ]
    }
   ],
   "source": [
    "import requests\n",
    "\n",
    "address = \"1 Brookings Dr, St. Louis, MO 63130\"\n",
    "\n",
    "response = requests.get('https://maps.googleapis.com/maps/api/geocode/json?address='+address)\n",
    "\n",
    "resp_json_payload = response.json()\n",
    "\n",
    "print(resp_json_payload['results'][0]['geometry']['location'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If latitude and longitude are simply fed into the neural network as two features, they might not be overly helpful.  These two values would allow your neural network to cluster locations on a map.  Sometimes cluster locations on a map can be useful.  Consider the percentage of the population that smokes in the USA by state:\n",
    "\n",
    "![Smokers by State](https://raw.githubusercontent.com/jeffheaton/t81_558_deep_learning/master/images/class_6_smokers.png \"Smokers by State\")\n",
    "\n",
    "The above map shows that certian behaviors, like smoking, can be clustered by global region. \n",
    "\n",
    "However, often you will want to transform the coordinates into distances.  It is reasonably easy to estimate the distance between any two points on Earth by using the [great circle distance](https://en.wikipedia.org/wiki/Great-circle_distance) between any two points on a sphere:\n",
    "\n",
    "The following code implements this formula:\n",
    "\n",
    "$\\Delta\\sigma=\\arccos\\bigl(\\sin\\phi_1\\cdot\\sin\\phi_2+\\cos\\phi_1\\cdot\\cos\\phi_2\\cdot\\cos(\\Delta\\lambda)\\bigr)$\n",
    "\n",
    "$d = r \\, \\Delta\\sigma$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Distance, St. Louis, MO to Ft. Lauderdale, FL: 1685.0869618595307 km\n"
     ]
    }
   ],
   "source": [
    "from math import sin, cos, sqrt, atan2, radians\n",
    "\n",
    "# Distance function\n",
    "def distance_lat_lng(lat1,lng1,lat2,lng2):\n",
    "    # approximate radius of earth in km\n",
    "    R = 6373.0\n",
    "\n",
    "    # degrees to radians (lat/lon are in degrees)\n",
    "    lat1 = radians(lat1)\n",
    "    lng1 = radians(lng1)\n",
    "    lat2 = radians(lat2)\n",
    "    lng2 = radians(lng2)\n",
    "\n",
    "    dlng = lng2 - lng1\n",
    "    dlat = lat2 - lat1\n",
    "\n",
    "    a = sin(dlat / 2)**2 + cos(lat1) * cos(lat2) * sin(dlng / 2)**2\n",
    "    c = 2 * atan2(sqrt(a), sqrt(1 - a))\n",
    "\n",
    "    return R * c\n",
    "\n",
    "# Find lat lon for address\n",
    "def lookup_lat_lng(address):\n",
    "    response = requests.get('https://maps.googleapis.com/maps/api/geocode/json?address='+address)\n",
    "    json = response.json()\n",
    "    if len(json['results']) == 0:\n",
    "        print(\"Can't find: {}\".format(address))\n",
    "        return 0,0\n",
    "    map = json['results'][0]['geometry']['location']\n",
    "    return map['lat'],map['lng']\n",
    "\n",
    "\n",
    "# Distance between two locations\n",
    "\n",
    "import requests\n",
    "\n",
    "address1 = \"1 Brookings Dr, St. Louis, MO 63130\" \n",
    "address2 = \"3301 College Ave, Fort Lauderdale, FL 33314\"\n",
    "\n",
    "lat1, lng1 = lookup_lat_lng(address1)\n",
    "lat2, lng2 = lookup_lat_lng(address2)\n",
    "\n",
    "print(\"Distance, St. Louis, MO to Ft. Lauderdale, FL: {} km\".format(\n",
    "        distance_lat_lng(lat1,lng1,lat2,lng2)))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Distances can be useful to encode addresses as.  You must consider what distance might be useful for your dataset.  Consider:\n",
    "\n",
    "* Distance to major metropolitan area\n",
    "* Distance to competitor\n",
    "* Distance to distribution center\n",
    "* Distance to retail outlet\n",
    "\n",
    "The following code calculates the distance between 10 universities and washu:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "School 'Princeton', distance to wustl is: 1354.7554246422899\n",
      "School 'Harvard', distance to wustl is: 1670.4803488515167\n",
      "School 'University of Chicago', distance to wustl is: 418.07074534593164\n",
      "School 'Yale', distance to wustl is: 1508.0574212290953\n",
      "School 'Columbia University', distance to wustl is: 1418.0702619177403\n",
      "School 'Stanford', distance to wustl is: 2781.0147376094565\n",
      "School 'MIT', distance to wustl is: 1672.3056259012085\n",
      "School 'Duke University', distance to wustl is: 1046.5733989340915\n",
      "School 'University of Pennsylvania', distance to wustl is: 1307.0188113752533\n",
      "School 'Johns Hopkins', distance to wustl is: 1184.1983102146019\n"
     ]
    }
   ],
   "source": [
    "# Encoding other universities by their distance to Washington University\n",
    "\n",
    "schools = [\n",
    "    [\"Princeton University, Princeton, NJ 08544\", 'Princeton'],\n",
    "    [\"Massachusetts Hall, Cambridge, MA 02138\", 'Harvard'],\n",
    "    [\"5801 S Ellis Ave, Chicago, IL 60637\", 'University of Chicago'],\n",
    "    [\"Yale, New Haven, CT 06520\", 'Yale'],\n",
    "    [\"116th St & Broadway, New York, NY 10027\", 'Columbia University'],\n",
    "    [\"450 Serra Mall, Stanford, CA 94305\", 'Stanford'],\n",
    "    [\"77 Massachusetts Ave, Cambridge, MA 02139\", 'MIT'],\n",
    "    [\"Duke University, Durham, NC 27708\", 'Duke University'],\n",
    "    [\"University of Pennsylvania, Philadelphia, PA 19104\", 'University of Pennsylvania'],\n",
    "    [\"Johns Hopkins University, Baltimore, MD 21218\", 'Johns Hopkins']\n",
    "]\n",
    "\n",
    "lat1, lng1 = lookup_lat_lng(\"1 Brookings Dr, St. Louis, MO 63130\")\n",
    "\n",
    "for address, name in schools:\n",
    "    lat2,lng2 = lookup_lat_lng(address)\n",
    "    dist = distance_lat_lng(lat1,lng1,lat2,lng2)\n",
    "    print(\"School '{}', distance to wustl is: {}\".format(name,dist))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "# Training with a Validation Set and Early Stopping\n",
    "\n",
    "**Overfitting** occurs when a neural network is trained to the point that it begins to memorize rather than generalize.  \n",
    "\n",
    "![Training vs Validation Error for Overfitting](https://raw.githubusercontent.com/jeffheaton/t81_558_deep_learning/master/images/class_3_training_val.png \"Training vs Validation Error for Overfitting\")\n",
    "\n",
    "It is important to segment the original dataset into several datasets:\n",
    "\n",
    "* **Training Set**\n",
    "* **Validation Set**\n",
    "* **Holdout Set**\n",
    "\n",
    "There are several different ways that these sets can be constructed.  The following programs demonstrate some of these.\n",
    "\n",
    "The first method is a training and validation set.  The training data are used to train the neural network until the validation set no longer improves.  This attempts to stop at a near optimal training point.  This method will only give accurate \"out of sample\" predictions for the validation set, this is usually 20% or so of the data.  The predictions for the training data will be overly optimistic, as these were the data that the neural network was trained on.  \n",
    "\n",
    "![Training with a Validation Set](https://raw.githubusercontent.com/jeffheaton/t81_558_deep_learning/master/images/class_1_train_val.png \"Training with a Validation Set\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Using TensorFlow backend.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Train on 150 samples, validate on 38 samples\n",
      "Epoch 1/1000\n",
      "0s - loss: 1.2092 - val_loss: 1.1779\n",
      "Epoch 2/1000\n",
      "0s - loss: 1.1924 - val_loss: 1.1629\n",
      "Epoch 3/1000\n",
      "0s - loss: 1.1771 - val_loss: 1.1493\n",
      "Epoch 4/1000\n",
      "0s - loss: 1.1629 - val_loss: 1.1371\n",
      "Epoch 5/1000\n",
      "0s - loss: 1.1511 - val_loss: 1.1262\n",
      "Epoch 6/1000\n",
      "0s - loss: 1.1386 - val_loss: 1.1160\n",
      "Epoch 7/1000\n",
      "0s - loss: 1.1279 - val_loss: 1.1055\n",
      "Epoch 8/1000\n",
      "0s - loss: 1.1162 - val_loss: 1.0952\n",
      "Epoch 9/1000\n",
      "0s - loss: 1.1061 - val_loss: 1.0853\n",
      "Epoch 10/1000\n",
      "0s - loss: 1.0962 - val_loss: 1.0751\n",
      "Epoch 11/1000\n",
      "0s - loss: 1.0850 - val_loss: 1.0664\n",
      "Epoch 12/1000\n",
      "0s - loss: 1.0751 - val_loss: 1.0567\n",
      "Epoch 13/1000\n",
      "0s - loss: 1.0643 - val_loss: 1.0463\n",
      "Epoch 14/1000\n",
      "0s - loss: 1.0536 - val_loss: 1.0356\n",
      "Epoch 15/1000\n",
      "0s - loss: 1.0426 - val_loss: 1.0244\n",
      "Epoch 16/1000\n",
      "0s - loss: 1.0322 - val_loss: 1.0126\n",
      "Epoch 17/1000\n",
      "0s - loss: 1.0201 - val_loss: 1.0002\n",
      "Epoch 18/1000\n",
      "0s - loss: 1.0073 - val_loss: 0.9848\n",
      "Epoch 19/1000\n",
      "0s - loss: 0.9924 - val_loss: 0.9692\n",
      "Epoch 20/1000\n",
      "0s - loss: 0.9760 - val_loss: 0.9547\n",
      "Epoch 21/1000\n",
      "0s - loss: 0.9632 - val_loss: 0.9415\n",
      "Epoch 22/1000\n",
      "0s - loss: 0.9506 - val_loss: 0.9297\n",
      "Epoch 23/1000\n",
      "0s - loss: 0.9389 - val_loss: 0.9162\n",
      "Epoch 24/1000\n",
      "0s - loss: 0.9260 - val_loss: 0.9032\n",
      "Epoch 25/1000\n",
      "0s - loss: 0.9136 - val_loss: 0.8903\n",
      "Epoch 26/1000\n",
      "0s - loss: 0.9009 - val_loss: 0.8770\n",
      "Epoch 27/1000\n",
      "0s - loss: 0.8880 - val_loss: 0.8628\n",
      "Epoch 28/1000\n",
      "0s - loss: 0.8744 - val_loss: 0.8486\n",
      "Epoch 29/1000\n",
      "0s - loss: 0.8611 - val_loss: 0.8339\n",
      "Epoch 30/1000\n",
      "0s - loss: 0.8470 - val_loss: 0.8195\n",
      "Epoch 31/1000\n",
      "0s - loss: 0.8335 - val_loss: 0.8045\n",
      "Epoch 32/1000\n",
      "0s - loss: 0.8192 - val_loss: 0.7898\n",
      "Epoch 33/1000\n",
      "0s - loss: 0.8052 - val_loss: 0.7750\n",
      "Epoch 34/1000\n",
      "0s - loss: 0.7908 - val_loss: 0.7605\n",
      "Epoch 35/1000\n",
      "0s - loss: 0.7770 - val_loss: 0.7451\n",
      "Epoch 36/1000\n",
      "0s - loss: 0.7625 - val_loss: 0.7303\n",
      "Epoch 37/1000\n",
      "0s - loss: 0.7484 - val_loss: 0.7155\n",
      "Epoch 38/1000\n",
      "0s - loss: 0.7342 - val_loss: 0.7011\n",
      "Epoch 39/1000\n",
      "0s - loss: 0.7211 - val_loss: 0.6868\n",
      "Epoch 40/1000\n",
      "0s - loss: 0.7065 - val_loss: 0.6724\n",
      "Epoch 41/1000\n",
      "0s - loss: 0.6931 - val_loss: 0.6582\n",
      "Epoch 42/1000\n",
      "0s - loss: 0.6796 - val_loss: 0.6444\n",
      "Epoch 43/1000\n",
      "0s - loss: 0.6665 - val_loss: 0.6311\n",
      "Epoch 44/1000\n",
      "0s - loss: 0.6536 - val_loss: 0.6179\n",
      "Epoch 45/1000\n",
      "0s - loss: 0.6410 - val_loss: 0.6051\n",
      "Epoch 46/1000\n",
      "0s - loss: 0.6287 - val_loss: 0.5927\n",
      "Epoch 47/1000\n",
      "0s - loss: 0.6167 - val_loss: 0.5807\n",
      "Epoch 48/1000\n",
      "0s - loss: 0.6053 - val_loss: 0.5686\n",
      "Epoch 49/1000\n",
      "0s - loss: 0.5941 - val_loss: 0.5572\n",
      "Epoch 50/1000\n",
      "0s - loss: 0.5831 - val_loss: 0.5463\n",
      "Epoch 51/1000\n",
      "0s - loss: 0.5727 - val_loss: 0.5363\n",
      "Epoch 52/1000\n",
      "0s - loss: 0.5626 - val_loss: 0.5259\n",
      "Epoch 53/1000\n",
      "0s - loss: 0.5527 - val_loss: 0.5164\n",
      "Epoch 54/1000\n",
      "0s - loss: 0.5432 - val_loss: 0.5066\n",
      "Epoch 55/1000\n",
      "0s - loss: 0.5340 - val_loss: 0.4976\n",
      "Epoch 56/1000\n",
      "0s - loss: 0.5252 - val_loss: 0.4890\n",
      "Epoch 57/1000\n",
      "0s - loss: 0.5168 - val_loss: 0.4808\n",
      "Epoch 58/1000\n",
      "0s - loss: 0.5086 - val_loss: 0.4731\n",
      "Epoch 59/1000\n",
      "0s - loss: 0.5008 - val_loss: 0.4655\n",
      "Epoch 60/1000\n",
      "0s - loss: 0.4931 - val_loss: 0.4579\n",
      "Epoch 61/1000\n",
      "0s - loss: 0.4857 - val_loss: 0.4507\n",
      "Epoch 62/1000\n",
      "0s - loss: 0.4785 - val_loss: 0.4438\n",
      "Epoch 63/1000\n",
      "0s - loss: 0.4716 - val_loss: 0.4371\n",
      "Epoch 64/1000\n",
      "0s - loss: 0.4649 - val_loss: 0.4307\n",
      "Epoch 65/1000\n",
      "0s - loss: 0.4586 - val_loss: 0.4244\n",
      "Epoch 66/1000\n",
      "0s - loss: 0.4521 - val_loss: 0.4184\n",
      "Epoch 67/1000\n",
      "0s - loss: 0.4459 - val_loss: 0.4126\n",
      "Epoch 68/1000\n",
      "0s - loss: 0.4401 - val_loss: 0.4070\n",
      "Epoch 69/1000\n",
      "0s - loss: 0.4342 - val_loss: 0.4015\n",
      "Epoch 70/1000\n",
      "0s - loss: 0.4284 - val_loss: 0.3958\n",
      "Epoch 71/1000\n",
      "0s - loss: 0.4230 - val_loss: 0.3904\n",
      "Epoch 72/1000\n",
      "0s - loss: 0.4171 - val_loss: 0.3853\n",
      "Epoch 73/1000\n",
      "0s - loss: 0.4117 - val_loss: 0.3801\n",
      "Epoch 74/1000\n",
      "0s - loss: 0.4065 - val_loss: 0.3752\n",
      "Epoch 75/1000\n",
      "0s - loss: 0.4012 - val_loss: 0.3705\n",
      "Epoch 76/1000\n",
      "0s - loss: 0.3962 - val_loss: 0.3658\n",
      "Epoch 77/1000\n",
      "0s - loss: 0.3909 - val_loss: 0.3613\n",
      "Epoch 78/1000\n",
      "0s - loss: 0.3872 - val_loss: 0.3570\n",
      "Epoch 79/1000\n",
      "0s - loss: 0.3813 - val_loss: 0.3520\n",
      "Epoch 80/1000\n",
      "0s - loss: 0.3772 - val_loss: 0.3474\n",
      "Epoch 81/1000\n",
      "0s - loss: 0.3719 - val_loss: 0.3431\n",
      "Epoch 82/1000\n",
      "0s - loss: 0.3670 - val_loss: 0.3388\n",
      "Epoch 83/1000\n",
      "0s - loss: 0.3622 - val_loss: 0.3346\n",
      "Epoch 84/1000\n",
      "0s - loss: 0.3578 - val_loss: 0.3301\n",
      "Epoch 85/1000\n",
      "0s - loss: 0.3534 - val_loss: 0.3259\n",
      "Epoch 86/1000\n",
      "0s - loss: 0.3489 - val_loss: 0.3219\n",
      "Epoch 87/1000\n",
      "0s - loss: 0.3442 - val_loss: 0.3177\n",
      "Epoch 88/1000\n",
      "0s - loss: 0.3405 - val_loss: 0.3138\n",
      "Epoch 89/1000\n",
      "0s - loss: 0.3356 - val_loss: 0.3098\n",
      "Epoch 90/1000\n",
      "0s - loss: 0.3312 - val_loss: 0.3061\n",
      "Epoch 91/1000\n",
      "0s - loss: 0.3278 - val_loss: 0.3021\n",
      "Epoch 92/1000\n",
      "0s - loss: 0.3232 - val_loss: 0.2985\n",
      "Epoch 93/1000\n",
      "0s - loss: 0.3187 - val_loss: 0.2945\n",
      "Epoch 94/1000\n",
      "0s - loss: 0.3146 - val_loss: 0.2908\n",
      "Epoch 95/1000\n",
      "0s - loss: 0.3108 - val_loss: 0.2871\n",
      "Epoch 96/1000\n",
      "0s - loss: 0.3065 - val_loss: 0.2835\n",
      "Epoch 97/1000\n",
      "0s - loss: 0.3026 - val_loss: 0.2800\n",
      "Epoch 98/1000\n",
      "0s - loss: 0.2986 - val_loss: 0.2764\n",
      "Epoch 99/1000\n",
      "0s - loss: 0.2958 - val_loss: 0.2733\n",
      "Epoch 100/1000\n",
      "0s - loss: 0.2913 - val_loss: 0.2695\n",
      "Epoch 101/1000\n",
      "0s - loss: 0.2879 - val_loss: 0.2668\n",
      "Epoch 102/1000\n",
      "0s - loss: 0.2835 - val_loss: 0.2632\n",
      "Epoch 103/1000\n",
      "0s - loss: 0.2809 - val_loss: 0.2597\n",
      "Epoch 104/1000\n",
      "0s - loss: 0.2772 - val_loss: 0.2566\n",
      "Epoch 105/1000\n",
      "0s - loss: 0.2730 - val_loss: 0.2532\n",
      "Epoch 106/1000\n",
      "0s - loss: 0.2690 - val_loss: 0.2500\n",
      "Epoch 107/1000\n",
      "0s - loss: 0.2656 - val_loss: 0.2473\n",
      "Epoch 108/1000\n",
      "0s - loss: 0.2624 - val_loss: 0.2446\n",
      "Epoch 109/1000\n",
      "0s - loss: 0.2589 - val_loss: 0.2411\n",
      "Epoch 110/1000\n",
      "0s - loss: 0.2560 - val_loss: 0.2383\n",
      "Epoch 111/1000\n",
      "0s - loss: 0.2528 - val_loss: 0.2352\n",
      "Epoch 112/1000\n",
      "0s - loss: 0.2494 - val_loss: 0.2323\n",
      "Epoch 113/1000\n",
      "0s - loss: 0.2464 - val_loss: 0.2295\n",
      "Epoch 114/1000\n",
      "0s - loss: 0.2425 - val_loss: 0.2273\n",
      "Epoch 115/1000\n",
      "0s - loss: 0.2399 - val_loss: 0.2253\n",
      "Epoch 116/1000\n",
      "0s - loss: 0.2375 - val_loss: 0.2226\n",
      "Epoch 117/1000\n",
      "0s - loss: 0.2339 - val_loss: 0.2190\n",
      "Epoch 118/1000\n",
      "0s - loss: 0.2310 - val_loss: 0.2161\n",
      "Epoch 119/1000\n",
      "0s - loss: 0.2284 - val_loss: 0.2136\n",
      "Epoch 120/1000\n",
      "0s - loss: 0.2261 - val_loss: 0.2111\n",
      "Epoch 121/1000\n",
      "0s - loss: 0.2219 - val_loss: 0.2098\n",
      "Epoch 122/1000\n",
      "0s - loss: 0.2213 - val_loss: 0.2095\n",
      "Epoch 123/1000\n",
      "0s - loss: 0.2185 - val_loss: 0.2048\n",
      "Epoch 124/1000\n",
      "0s - loss: 0.2152 - val_loss: 0.2022\n",
      "Epoch 125/1000\n",
      "0s - loss: 0.2120 - val_loss: 0.1996\n",
      "Epoch 126/1000\n",
      "0s - loss: 0.2111 - val_loss: 0.1973\n",
      "Epoch 127/1000\n",
      "0s - loss: 0.2079 - val_loss: 0.1952\n",
      "Epoch 128/1000\n",
      "0s - loss: 0.2049 - val_loss: 0.1935\n",
      "Epoch 129/1000\n",
      "0s - loss: 0.2029 - val_loss: 0.1914\n",
      "Epoch 130/1000\n",
      "0s - loss: 0.2006 - val_loss: 0.1894\n",
      "Epoch 131/1000\n",
      "0s - loss: 0.1987 - val_loss: 0.1878\n",
      "Epoch 132/1000\n",
      "0s - loss: 0.1982 - val_loss: 0.1873\n",
      "Epoch 133/1000\n",
      "0s - loss: 0.1941 - val_loss: 0.1832\n",
      "Epoch 134/1000\n",
      "0s - loss: 0.1915 - val_loss: 0.1809\n",
      "Epoch 135/1000\n",
      "0s - loss: 0.1903 - val_loss: 0.1790\n",
      "Epoch 136/1000\n",
      "0s - loss: 0.1886 - val_loss: 0.1774\n",
      "Epoch 137/1000\n",
      "0s - loss: 0.1859 - val_loss: 0.1757\n",
      "Epoch 138/1000\n",
      "0s - loss: 0.1836 - val_loss: 0.1739\n",
      "Epoch 139/1000\n",
      "0s - loss: 0.1818 - val_loss: 0.1723\n",
      "Epoch 140/1000\n",
      "0s - loss: 0.1799 - val_loss: 0.1707\n",
      "Epoch 141/1000\n",
      "0s - loss: 0.1782 - val_loss: 0.1687\n",
      "Epoch 142/1000\n",
      "0s - loss: 0.1763 - val_loss: 0.1674\n",
      "Epoch 143/1000\n",
      "0s - loss: 0.1745 - val_loss: 0.1658\n",
      "Epoch 144/1000\n",
      "0s - loss: 0.1733 - val_loss: 0.1644\n",
      "Epoch 145/1000\n",
      "0s - loss: 0.1710 - val_loss: 0.1619\n",
      "Epoch 146/1000\n",
      "0s - loss: 0.1701 - val_loss: 0.1608\n",
      "Epoch 147/1000\n",
      "0s - loss: 0.1698 - val_loss: 0.1587\n",
      "Epoch 148/1000\n",
      "0s - loss: 0.1667 - val_loss: 0.1584\n",
      "Epoch 149/1000\n",
      "0s - loss: 0.1653 - val_loss: 0.1579\n",
      "Epoch 150/1000\n",
      "0s - loss: 0.1632 - val_loss: 0.1550\n",
      "Epoch 151/1000\n",
      "0s - loss: 0.1621 - val_loss: 0.1529\n",
      "Epoch 152/1000\n",
      "0s - loss: 0.1600 - val_loss: 0.1519\n",
      "Epoch 153/1000\n",
      "0s - loss: 0.1582 - val_loss: 0.1511\n",
      "Epoch 154/1000\n",
      "0s - loss: 0.1568 - val_loss: 0.1509\n",
      "Epoch 155/1000\n",
      "0s - loss: 0.1560 - val_loss: 0.1489\n",
      "Epoch 156/1000\n",
      "0s - loss: 0.1542 - val_loss: 0.1473\n",
      "Epoch 157/1000\n",
      "0s - loss: 0.1529 - val_loss: 0.1458\n",
      "Epoch 158/1000\n",
      "0s - loss: 0.1516 - val_loss: 0.1445\n",
      "Epoch 159/1000\n",
      "0s - loss: 0.1504 - val_loss: 0.1427\n",
      "Epoch 160/1000\n",
      "0s - loss: 0.1489 - val_loss: 0.1416\n",
      "Epoch 161/1000\n",
      "0s - loss: 0.1477 - val_loss: 0.1415\n",
      "Epoch 162/1000\n",
      "0s - loss: 0.1472 - val_loss: 0.1414\n",
      "Epoch 163/1000\n",
      "0s - loss: 0.1454 - val_loss: 0.1383\n",
      "Epoch 164/1000\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0s - loss: 0.1438 - val_loss: 0.1370\n",
      "Epoch 165/1000\n",
      "0s - loss: 0.1428 - val_loss: 0.1361\n",
      "Epoch 166/1000\n",
      "0s - loss: 0.1414 - val_loss: 0.1347\n",
      "Epoch 167/1000\n",
      "0s - loss: 0.1403 - val_loss: 0.1335\n",
      "Epoch 168/1000\n",
      "0s - loss: 0.1393 - val_loss: 0.1327\n",
      "Epoch 169/1000\n",
      "0s - loss: 0.1380 - val_loss: 0.1320\n",
      "Epoch 170/1000\n",
      "0s - loss: 0.1374 - val_loss: 0.1322\n",
      "Epoch 171/1000\n",
      "0s - loss: 0.1363 - val_loss: 0.1298\n",
      "Epoch 172/1000\n",
      "0s - loss: 0.1349 - val_loss: 0.1286\n",
      "Epoch 173/1000\n",
      "0s - loss: 0.1337 - val_loss: 0.1273\n",
      "Epoch 174/1000\n",
      "0s - loss: 0.1328 - val_loss: 0.1263\n",
      "Epoch 175/1000\n",
      "0s - loss: 0.1322 - val_loss: 0.1259\n",
      "Epoch 176/1000\n",
      "0s - loss: 0.1308 - val_loss: 0.1251\n",
      "Epoch 177/1000\n",
      "0s - loss: 0.1302 - val_loss: 0.1233\n",
      "Epoch 178/1000\n",
      "0s - loss: 0.1294 - val_loss: 0.1221\n",
      "Epoch 179/1000\n",
      "0s - loss: 0.1278 - val_loss: 0.1224\n",
      "Epoch 180/1000\n",
      "0s - loss: 0.1274 - val_loss: 0.1229\n",
      "Epoch 181/1000\n",
      "0s - loss: 0.1264 - val_loss: 0.1205\n",
      "Epoch 182/1000\n",
      "0s - loss: 0.1253 - val_loss: 0.1190\n",
      "Epoch 183/1000\n",
      "0s - loss: 0.1244 - val_loss: 0.1182\n",
      "Epoch 184/1000\n",
      "0s - loss: 0.1239 - val_loss: 0.1176\n",
      "Epoch 185/1000\n",
      "0s - loss: 0.1228 - val_loss: 0.1165\n",
      "Epoch 186/1000\n",
      "0s - loss: 0.1226 - val_loss: 0.1149\n",
      "Epoch 187/1000\n",
      "0s - loss: 0.1211 - val_loss: 0.1149\n",
      "Epoch 188/1000\n",
      "0s - loss: 0.1204 - val_loss: 0.1154\n",
      "Epoch 189/1000\n",
      "0s - loss: 0.1198 - val_loss: 0.1153\n",
      "Epoch 190/1000\n",
      "0s - loss: 0.1205 - val_loss: 0.1153\n",
      "Epoch 191/1000\n",
      "0s - loss: 0.1179 - val_loss: 0.1109\n",
      "Epoch 192/1000\n",
      "0s - loss: 0.1181 - val_loss: 0.1097\n",
      "Epoch 193/1000\n",
      "0s - loss: 0.1175 - val_loss: 0.1094\n",
      "Epoch 194/1000\n",
      "0s - loss: 0.1162 - val_loss: 0.1101\n",
      "Epoch 195/1000\n",
      "0s - loss: 0.1149 - val_loss: 0.1100\n",
      "Epoch 196/1000\n",
      "0s - loss: 0.1154 - val_loss: 0.1105\n",
      "Epoch 197/1000\n",
      "0s - loss: 0.1150 - val_loss: 0.1068\n",
      "Epoch 198/1000\n",
      "0s - loss: 0.1130 - val_loss: 0.1061\n",
      "Epoch 199/1000\n",
      "0s - loss: 0.1130 - val_loss: 0.1061\n",
      "Epoch 200/1000\n",
      "0s - loss: 0.1126 - val_loss: 0.1045\n",
      "Epoch 201/1000\n",
      "0s - loss: 0.1108 - val_loss: 0.1054\n",
      "Epoch 202/1000\n",
      "0s - loss: 0.1105 - val_loss: 0.1071\n",
      "Epoch 203/1000\n",
      "0s - loss: 0.1104 - val_loss: 0.1062\n",
      "Epoch 204/1000\n",
      "0s - loss: 0.1090 - val_loss: 0.1033\n",
      "Epoch 205/1000\n",
      "0s - loss: 0.1084 - val_loss: 0.1011\n",
      "Epoch 206/1000\n",
      "0s - loss: 0.1087 - val_loss: 0.1003\n",
      "Epoch 207/1000\n",
      "0s - loss: 0.1083 - val_loss: 0.1020\n",
      "Epoch 208/1000\n",
      "0s - loss: 0.1068 - val_loss: 0.1016\n",
      "Epoch 209/1000\n",
      "0s - loss: 0.1064 - val_loss: 0.1006\n",
      "Epoch 210/1000\n",
      "0s - loss: 0.1058 - val_loss: 0.0997\n",
      "Epoch 211/1000\n",
      "0s - loss: 0.1064 - val_loss: 0.0979\n",
      "Epoch 212/1000\n",
      "0s - loss: 0.1049 - val_loss: 0.0994\n",
      "Epoch 213/1000\n",
      "0s - loss: 0.1044 - val_loss: 0.0982\n",
      "Epoch 214/1000\n",
      "0s - loss: 0.1038 - val_loss: 0.0987\n",
      "Epoch 215/1000\n",
      "0s - loss: 0.1032 - val_loss: 0.0977\n",
      "Epoch 216/1000\n",
      "0s - loss: 0.1024 - val_loss: 0.0961\n",
      "Epoch 217/1000\n",
      "0s - loss: 0.1021 - val_loss: 0.0943\n",
      "Epoch 218/1000\n",
      "0s - loss: 0.1020 - val_loss: 0.0944\n",
      "Epoch 219/1000\n",
      "0s - loss: 0.1011 - val_loss: 0.0938\n",
      "Epoch 220/1000\n",
      "0s - loss: 0.1022 - val_loss: 0.0952\n",
      "Epoch 221/1000\n",
      "0s - loss: 0.1002 - val_loss: 0.0932\n",
      "Epoch 222/1000\n",
      "0s - loss: 0.0999 - val_loss: 0.0928\n",
      "Epoch 223/1000\n",
      "0s - loss: 0.0991 - val_loss: 0.0911\n",
      "Epoch 224/1000\n",
      "0s - loss: 0.0989 - val_loss: 0.0908\n",
      "Epoch 225/1000\n",
      "0s - loss: 0.0986 - val_loss: 0.0905\n",
      "Epoch 226/1000\n",
      "0s - loss: 0.0972 - val_loss: 0.0927\n",
      "Epoch 227/1000\n",
      "0s - loss: 0.0994 - val_loss: 0.0953\n",
      "Epoch 228/1000\n",
      "0s - loss: 0.0975 - val_loss: 0.0916\n",
      "Epoch 229/1000\n",
      "0s - loss: 0.0982 - val_loss: 0.0876\n",
      "Epoch 230/1000\n",
      "0s - loss: 0.0968 - val_loss: 0.0876\n",
      "Epoch 231/1000\n",
      "0s - loss: 0.0957 - val_loss: 0.0877\n",
      "Epoch 232/1000\n",
      "0s - loss: 0.0950 - val_loss: 0.0883\n",
      "Epoch 233/1000\n",
      "0s - loss: 0.0956 - val_loss: 0.0908\n",
      "Epoch 234/1000\n",
      "0s - loss: 0.0961 - val_loss: 0.0912\n",
      "Epoch 235/1000\n",
      "0s - loss: 0.0944 - val_loss: 0.0863\n",
      "Epoch 236/1000\n",
      "0s - loss: 0.0940 - val_loss: 0.0843\n",
      "Epoch 237/1000\n",
      "0s - loss: 0.0940 - val_loss: 0.0847\n",
      "Epoch 238/1000\n",
      "0s - loss: 0.0930 - val_loss: 0.0849\n",
      "Epoch 239/1000\n",
      "0s - loss: 0.0928 - val_loss: 0.0861\n",
      "Epoch 240/1000\n",
      "0s - loss: 0.0926 - val_loss: 0.0844\n",
      "Epoch 241/1000\n",
      "0s - loss: 0.0919 - val_loss: 0.0847\n",
      "Epoch 242/1000\n",
      "0s - loss: 0.0918 - val_loss: 0.0832\n",
      "Epoch 243/1000\n",
      "0s - loss: 0.0911 - val_loss: 0.0833\n",
      "Epoch 244/1000\n",
      "0s - loss: 0.0911 - val_loss: 0.0824\n",
      "Epoch 245/1000\n",
      "0s - loss: 0.0905 - val_loss: 0.0828\n",
      "Epoch 246/1000\n",
      "0s - loss: 0.0903 - val_loss: 0.0842\n",
      "Epoch 247/1000\n",
      "0s - loss: 0.0899 - val_loss: 0.0826\n",
      "Epoch 248/1000\n",
      "0s - loss: 0.0894 - val_loss: 0.0811\n",
      "Epoch 249/1000\n",
      "0s - loss: 0.0893 - val_loss: 0.0806\n",
      "Epoch 250/1000\n",
      "0s - loss: 0.0889 - val_loss: 0.0818\n",
      "Epoch 251/1000\n",
      "0s - loss: 0.0887 - val_loss: 0.0811\n",
      "Epoch 252/1000\n",
      "0s - loss: 0.0883 - val_loss: 0.0822\n",
      "Epoch 253/1000\n",
      "0s - loss: 0.0882 - val_loss: 0.0803\n",
      "Epoch 254/1000\n",
      "0s - loss: 0.0875 - val_loss: 0.0798\n",
      "Epoch 255/1000\n",
      "0s - loss: 0.0873 - val_loss: 0.0788\n",
      "Epoch 256/1000\n",
      "0s - loss: 0.0869 - val_loss: 0.0788\n",
      "Epoch 257/1000\n",
      "0s - loss: 0.0867 - val_loss: 0.0783\n",
      "Epoch 258/1000\n",
      "0s - loss: 0.0863 - val_loss: 0.0776\n",
      "Epoch 259/1000\n",
      "0s - loss: 0.0867 - val_loss: 0.0761\n",
      "Epoch 260/1000\n",
      "0s - loss: 0.0858 - val_loss: 0.0774\n",
      "Epoch 261/1000\n",
      "0s - loss: 0.0860 - val_loss: 0.0782\n",
      "Epoch 262/1000\n",
      "0s - loss: 0.0857 - val_loss: 0.0771\n",
      "Epoch 263/1000\n",
      "0s - loss: 0.0855 - val_loss: 0.0754\n",
      "Epoch 264/1000\n",
      "0s - loss: 0.0849 - val_loss: 0.0754\n",
      "Epoch 265/1000\n",
      "0s - loss: 0.0849 - val_loss: 0.0760\n",
      "Epoch 00264: early stopping\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<keras.callbacks.History at 0x1797ea40c18>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import io\n",
    "import requests\n",
    "import numpy as np\n",
    "import os\n",
    "from sklearn.model_selection import train_test_split\n",
    "from sklearn import metrics\n",
    "from keras.models import Sequential\n",
    "from keras.layers.core import Dense, Activation\n",
    "from keras.callbacks import EarlyStopping\n",
    "\n",
    "path = \"./data/\"\n",
    "    \n",
    "filename = os.path.join(path,\"iris.csv\")    \n",
    "df = pd.read_csv(filename,na_values=['NA','?'])\n",
    "\n",
    "species = encode_text_index(df,\"species\")\n",
    "x,y = to_xy(df,\"species\")\n",
    "\n",
    "# Split into train/test\n",
    "x_train, x_test, y_train, y_test = train_test_split(    \n",
    "    x, y, test_size=0.25, random_state=42)\n",
    "\n",
    "model = Sequential()\n",
    "model.add(Dense(10, input_dim=x.shape[1], activation='relu'))\n",
    "model.add(Dense(5,activation='relu'))\n",
    "model.add(Dense(y.shape[1],activation='softmax'))\n",
    "model.compile(loss='categorical_crossentropy', optimizer='adam')\n",
    "monitor = EarlyStopping(monitor='val_loss', min_delta=1e-3, patience=5, verbose=1, mode='auto')\n",
    "\n",
    "model.fit(x,y,validation_data=(x_test,y_test),callbacks=[monitor],verbose=2,epochs=1000)\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now that the neural network is trained, we can make predictions about the test set.  The following code predicts the type of iris for test set and displays the first five irises. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[  1.13033631e-03   9.03616965e-01   9.52527672e-02]\n",
      " [  9.96876478e-01   3.12355976e-03   1.00982956e-13]\n",
      " [  4.28446478e-09   2.90802354e-03   9.97092009e-01]\n",
      " [  9.09838884e-04   8.96424115e-01   1.02666058e-01]\n",
      " [  8.87830800e-04   9.41267371e-01   5.78448363e-02]]\n"
     ]
    }
   ],
   "source": [
    "from sklearn import metrics\n",
    "import tensorflow as tf\n",
    "\n",
    "pred = model.predict(x_test)\n",
    "print(pred[0:5]) # print first five predictions"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These numbers are in scientific notation.  Each line provides the probability that the iris is one of the 3 types of iris in the data set.  For the first line, the second type of iris has a 91% probability of being the species of iris."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Early Stopping and the Best Weights\n",
    "\n",
    "In the previous section we used early stopping so that the training would halt once the validation set no longer saw score improvements for a number of steps.  This number of steps that early stopping will tolerate no improvement is called *patience*.  If the patience value is large, the neural network's error may continue to worsen while early stopping is patiently waiting.  At some point earlier in the training the optimal set of weights was obtained for the neural network.  However, at the end of training we will have the weights for the neural network that finally exhausted the patience of early stopping.  The weights of this neural network might not be bad, but it would be better to have the most optimal weights during the entire training operation.  \n",
    "\n",
    "The code presented below does this.  An additional monitor is used and saves a copy of the neural network to **best_weights.hdf5** each time the validation score of the neural network improves.  Once training is done, we just reload this file and we have the optimal training weights that were found.\n",
    "\n",
    "This technique is slight overkill for many of the examples for this class.  It can also introduce an occasional issue (as described in the next section).  Because of this, most of the examples in this course will not use this code to obtain the absolute best weights.  However, for the larger, more complex datasets in this course, we will save the absolute best weights as demonstrated here.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Train on 150 samples, validate on 38 samples\n",
      "Epoch 1/1000\n",
      "0s - loss: 2.1069 - val_loss: 1.9295\n",
      "Epoch 2/1000\n",
      "0s - loss: 1.9730 - val_loss: 1.8234\n",
      "Epoch 3/1000\n",
      "0s - loss: 1.8565 - val_loss: 1.7278\n",
      "Epoch 4/1000\n",
      "0s - loss: 1.7548 - val_loss: 1.6441\n",
      "Epoch 5/1000\n",
      "0s - loss: 1.6673 - val_loss: 1.5716\n",
      "Epoch 6/1000\n",
      "0s - loss: 1.5888 - val_loss: 1.5084\n",
      "Epoch 7/1000\n",
      "0s - loss: 1.5213 - val_loss: 1.4517\n",
      "Epoch 8/1000\n",
      "0s - loss: 1.4554 - val_loss: 1.3965\n",
      "Epoch 9/1000\n",
      "0s - loss: 1.3950 - val_loss: 1.3433\n",
      "Epoch 10/1000\n",
      "0s - loss: 1.3374 - val_loss: 1.2913\n",
      "Epoch 11/1000\n",
      "0s - loss: 1.2832 - val_loss: 1.2432\n",
      "Epoch 12/1000\n",
      "0s - loss: 1.2313 - val_loss: 1.1982\n",
      "Epoch 13/1000\n",
      "0s - loss: 1.1849 - val_loss: 1.1570\n",
      "Epoch 14/1000\n",
      "0s - loss: 1.1478 - val_loss: 1.1307\n",
      "Epoch 15/1000\n",
      "0s - loss: 1.1219 - val_loss: 1.1130\n",
      "Epoch 16/1000\n",
      "0s - loss: 1.1072 - val_loss: 1.1070\n",
      "Epoch 17/1000\n",
      "0s - loss: 1.1001 - val_loss: 1.1030\n",
      "Epoch 18/1000\n",
      "0s - loss: 1.0963 - val_loss: 1.0996\n",
      "Epoch 19/1000\n",
      "0s - loss: 1.0923 - val_loss: 1.0969\n",
      "Epoch 20/1000\n",
      "0s - loss: 1.0901 - val_loss: 1.0936\n",
      "Epoch 21/1000\n",
      "0s - loss: 1.0874 - val_loss: 1.0921\n",
      "Epoch 22/1000\n",
      "0s - loss: 1.0858 - val_loss: 1.0910\n",
      "Epoch 23/1000\n",
      "0s - loss: 1.0846 - val_loss: 1.0900\n",
      "Epoch 24/1000\n",
      "0s - loss: 1.0833 - val_loss: 1.0889\n",
      "Epoch 25/1000\n",
      "0s - loss: 1.0823 - val_loss: 1.0881\n",
      "Epoch 26/1000\n",
      "0s - loss: 1.0812 - val_loss: 1.0870\n",
      "Epoch 27/1000\n",
      "0s - loss: 1.0800 - val_loss: 1.0859\n",
      "Epoch 28/1000\n",
      "0s - loss: 1.0792 - val_loss: 1.0851\n",
      "Epoch 29/1000\n",
      "0s - loss: 1.0784 - val_loss: 1.0842\n",
      "Epoch 30/1000\n",
      "0s - loss: 1.0772 - val_loss: 1.0831\n",
      "Epoch 31/1000\n",
      "0s - loss: 1.0763 - val_loss: 1.0821\n",
      "Epoch 32/1000\n",
      "0s - loss: 1.0752 - val_loss: 1.0811\n",
      "Epoch 33/1000\n",
      "0s - loss: 1.0744 - val_loss: 1.0802\n",
      "Epoch 34/1000\n",
      "0s - loss: 1.0735 - val_loss: 1.0794\n",
      "Epoch 35/1000\n",
      "0s - loss: 1.0726 - val_loss: 1.0785\n",
      "Epoch 36/1000\n",
      "0s - loss: 1.0718 - val_loss: 1.0776\n",
      "Epoch 37/1000\n",
      "0s - loss: 1.0710 - val_loss: 1.0767\n",
      "Epoch 38/1000\n",
      "0s - loss: 1.0703 - val_loss: 1.0759\n",
      "Epoch 39/1000\n",
      "0s - loss: 1.0694 - val_loss: 1.0751\n",
      "Epoch 40/1000\n",
      "0s - loss: 1.0687 - val_loss: 1.0743\n",
      "Epoch 41/1000\n",
      "0s - loss: 1.0678 - val_loss: 1.0736\n",
      "Epoch 42/1000\n",
      "0s - loss: 1.0670 - val_loss: 1.0728\n",
      "Epoch 43/1000\n",
      "0s - loss: 1.0661 - val_loss: 1.0720\n",
      "Epoch 44/1000\n",
      "0s - loss: 1.0653 - val_loss: 1.0711\n",
      "Epoch 45/1000\n",
      "0s - loss: 1.0647 - val_loss: 1.0703\n",
      "Epoch 46/1000\n",
      "0s - loss: 1.0636 - val_loss: 1.0694\n",
      "Epoch 47/1000\n",
      "0s - loss: 1.0627 - val_loss: 1.0685\n",
      "Epoch 48/1000\n",
      "0s - loss: 1.0619 - val_loss: 1.0678\n",
      "Epoch 49/1000\n",
      "0s - loss: 1.0611 - val_loss: 1.0668\n",
      "Epoch 50/1000\n",
      "0s - loss: 1.0605 - val_loss: 1.0661\n",
      "Epoch 51/1000\n",
      "0s - loss: 1.0592 - val_loss: 1.0650\n",
      "Epoch 52/1000\n",
      "0s - loss: 1.0584 - val_loss: 1.0641\n",
      "Epoch 53/1000\n",
      "0s - loss: 1.0575 - val_loss: 1.0632\n",
      "Epoch 54/1000\n",
      "0s - loss: 1.0565 - val_loss: 1.0622\n",
      "Epoch 55/1000\n",
      "0s - loss: 1.0554 - val_loss: 1.0612\n",
      "Epoch 56/1000\n",
      "0s - loss: 1.0553 - val_loss: 1.0604\n",
      "Epoch 57/1000\n",
      "0s - loss: 1.0536 - val_loss: 1.0597\n",
      "Epoch 58/1000\n",
      "0s - loss: 1.0527 - val_loss: 1.0589\n",
      "Epoch 59/1000\n",
      "0s - loss: 1.0517 - val_loss: 1.0581\n",
      "Epoch 60/1000\n",
      "0s - loss: 1.0507 - val_loss: 1.0570\n",
      "Epoch 61/1000\n",
      "0s - loss: 1.0497 - val_loss: 1.0556\n",
      "Epoch 62/1000\n",
      "0s - loss: 1.0487 - val_loss: 1.0546\n",
      "Epoch 63/1000\n",
      "0s - loss: 1.0485 - val_loss: 1.0535\n",
      "Epoch 64/1000\n",
      "0s - loss: 1.0467 - val_loss: 1.0524\n",
      "Epoch 65/1000\n",
      "0s - loss: 1.0455 - val_loss: 1.0513\n",
      "Epoch 66/1000\n",
      "0s - loss: 1.0444 - val_loss: 1.0507\n",
      "Epoch 67/1000\n",
      "0s - loss: 1.0436 - val_loss: 1.0500\n",
      "Epoch 68/1000\n",
      "0s - loss: 1.0424 - val_loss: 1.0486\n",
      "Epoch 69/1000\n",
      "0s - loss: 1.0412 - val_loss: 1.0473\n",
      "Epoch 70/1000\n",
      "0s - loss: 1.0401 - val_loss: 1.0461\n",
      "Epoch 71/1000\n",
      "0s - loss: 1.0390 - val_loss: 1.0451\n",
      "Epoch 72/1000\n",
      "0s - loss: 1.0378 - val_loss: 1.0439\n",
      "Epoch 73/1000\n",
      "0s - loss: 1.0367 - val_loss: 1.0426\n",
      "Epoch 74/1000\n",
      "0s - loss: 1.0355 - val_loss: 1.0414\n",
      "Epoch 75/1000\n",
      "0s - loss: 1.0345 - val_loss: 1.0403\n",
      "Epoch 76/1000\n",
      "0s - loss: 1.0333 - val_loss: 1.0394\n",
      "Epoch 77/1000\n",
      "0s - loss: 1.0320 - val_loss: 1.0381\n",
      "Epoch 78/1000\n",
      "0s - loss: 1.0309 - val_loss: 1.0372\n",
      "Epoch 79/1000\n",
      "0s - loss: 1.0294 - val_loss: 1.0357\n",
      "Epoch 80/1000\n",
      "0s - loss: 1.0282 - val_loss: 1.0343\n",
      "Epoch 81/1000\n",
      "0s - loss: 1.0270 - val_loss: 1.0329\n",
      "Epoch 82/1000\n",
      "0s - loss: 1.0257 - val_loss: 1.0322\n",
      "Epoch 83/1000\n",
      "0s - loss: 1.0252 - val_loss: 1.0314\n",
      "Epoch 84/1000\n",
      "0s - loss: 1.0229 - val_loss: 1.0297\n",
      "Epoch 85/1000\n",
      "0s - loss: 1.0225 - val_loss: 1.0276\n",
      "Epoch 86/1000\n",
      "0s - loss: 1.0202 - val_loss: 1.0263\n",
      "Epoch 87/1000\n",
      "0s - loss: 1.0190 - val_loss: 1.0253\n",
      "Epoch 88/1000\n",
      "0s - loss: 1.0174 - val_loss: 1.0240\n",
      "Epoch 89/1000\n",
      "0s - loss: 1.0160 - val_loss: 1.0227\n",
      "Epoch 90/1000\n",
      "0s - loss: 1.0145 - val_loss: 1.0210\n",
      "Epoch 91/1000\n",
      "0s - loss: 1.0141 - val_loss: 1.0192\n",
      "Epoch 92/1000\n",
      "0s - loss: 1.0119 - val_loss: 1.0182\n",
      "Epoch 93/1000\n",
      "0s - loss: 1.0102 - val_loss: 1.0168\n",
      "Epoch 94/1000\n",
      "0s - loss: 1.0086 - val_loss: 1.0155\n",
      "Epoch 95/1000\n",
      "0s - loss: 1.0070 - val_loss: 1.0142\n",
      "Epoch 96/1000\n",
      "0s - loss: 1.0057 - val_loss: 1.0130\n",
      "Epoch 97/1000\n",
      "0s - loss: 1.0041 - val_loss: 1.0112\n",
      "Epoch 98/1000\n",
      "0s - loss: 1.0023 - val_loss: 1.0094\n",
      "Epoch 99/1000\n",
      "0s - loss: 1.0008 - val_loss: 1.0073\n",
      "Epoch 100/1000\n",
      "0s - loss: 0.9995 - val_loss: 1.0055\n",
      "Epoch 101/1000\n",
      "0s - loss: 0.9981 - val_loss: 1.0038\n",
      "Epoch 102/1000\n",
      "0s - loss: 0.9955 - val_loss: 1.0026\n",
      "Epoch 103/1000\n",
      "0s - loss: 0.9942 - val_loss: 1.0017\n",
      "Epoch 104/1000\n",
      "0s - loss: 0.9925 - val_loss: 1.0005\n",
      "Epoch 105/1000\n",
      "0s - loss: 0.9911 - val_loss: 0.9988\n",
      "Epoch 106/1000\n",
      "0s - loss: 0.9891 - val_loss: 0.9964\n",
      "Epoch 107/1000\n",
      "0s - loss: 0.9873 - val_loss: 0.9945\n",
      "Epoch 108/1000\n",
      "0s - loss: 0.9853 - val_loss: 0.9926\n",
      "Epoch 109/1000\n",
      "0s - loss: 0.9842 - val_loss: 0.9904\n",
      "Epoch 110/1000\n",
      "0s - loss: 0.9820 - val_loss: 0.9886\n",
      "Epoch 111/1000\n",
      "0s - loss: 0.9801 - val_loss: 0.9867\n",
      "Epoch 112/1000\n",
      "0s - loss: 0.9782 - val_loss: 0.9851\n",
      "Epoch 113/1000\n",
      "0s - loss: 0.9768 - val_loss: 0.9839\n",
      "Epoch 114/1000\n",
      "0s - loss: 0.9744 - val_loss: 0.9821\n",
      "Epoch 115/1000\n",
      "0s - loss: 0.9729 - val_loss: 0.9798\n",
      "Epoch 116/1000\n",
      "0s - loss: 0.9706 - val_loss: 0.9780\n",
      "Epoch 117/1000\n",
      "0s - loss: 0.9686 - val_loss: 0.9767\n",
      "Epoch 118/1000\n",
      "0s - loss: 0.9669 - val_loss: 0.9745\n",
      "Epoch 119/1000\n",
      "0s - loss: 0.9646 - val_loss: 0.9728\n",
      "Epoch 120/1000\n",
      "0s - loss: 0.9629 - val_loss: 0.9707\n",
      "Epoch 121/1000\n",
      "0s - loss: 0.9605 - val_loss: 0.9691\n",
      "Epoch 122/1000\n",
      "0s - loss: 0.9587 - val_loss: 0.9670\n",
      "Epoch 123/1000\n",
      "0s - loss: 0.9566 - val_loss: 0.9647\n",
      "Epoch 124/1000\n",
      "0s - loss: 0.9546 - val_loss: 0.9626\n",
      "Epoch 125/1000\n",
      "0s - loss: 0.9523 - val_loss: 0.9606\n",
      "Epoch 126/1000\n",
      "0s - loss: 0.9503 - val_loss: 0.9586\n",
      "Epoch 127/1000\n",
      "0s - loss: 0.9483 - val_loss: 0.9567\n",
      "Epoch 128/1000\n",
      "0s - loss: 0.9461 - val_loss: 0.9547\n",
      "Epoch 129/1000\n",
      "0s - loss: 0.9439 - val_loss: 0.9523\n",
      "Epoch 130/1000\n",
      "0s - loss: 0.9423 - val_loss: 0.9505\n",
      "Epoch 131/1000\n",
      "0s - loss: 0.9396 - val_loss: 0.9480\n",
      "Epoch 132/1000\n",
      "0s - loss: 0.9375 - val_loss: 0.9459\n",
      "Epoch 133/1000\n",
      "0s - loss: 0.9354 - val_loss: 0.9439\n",
      "Epoch 134/1000\n",
      "0s - loss: 0.9333 - val_loss: 0.9418\n",
      "Epoch 135/1000\n",
      "0s - loss: 0.9309 - val_loss: 0.9405\n",
      "Epoch 136/1000\n",
      "0s - loss: 0.9289 - val_loss: 0.9384\n",
      "Epoch 137/1000\n",
      "0s - loss: 0.9270 - val_loss: 0.9353\n",
      "Epoch 138/1000\n",
      "0s - loss: 0.9241 - val_loss: 0.9330\n",
      "Epoch 139/1000\n",
      "0s - loss: 0.9220 - val_loss: 0.9309\n",
      "Epoch 140/1000\n",
      "0s - loss: 0.9197 - val_loss: 0.9286\n",
      "Epoch 141/1000\n",
      "0s - loss: 0.9177 - val_loss: 0.9267\n",
      "Epoch 142/1000\n",
      "0s - loss: 0.9152 - val_loss: 0.9245\n",
      "Epoch 143/1000\n",
      "0s - loss: 0.9136 - val_loss: 0.9225\n",
      "Epoch 144/1000\n",
      "0s - loss: 0.9105 - val_loss: 0.9200\n",
      "Epoch 145/1000\n",
      "0s - loss: 0.9093 - val_loss: 0.9172\n",
      "Epoch 146/1000\n",
      "0s - loss: 0.9064 - val_loss: 0.9151\n",
      "Epoch 147/1000\n",
      "0s - loss: 0.9044 - val_loss: 0.9133\n",
      "Epoch 148/1000\n",
      "0s - loss: 0.9018 - val_loss: 0.9115\n",
      "Epoch 149/1000\n",
      "0s - loss: 0.8993 - val_loss: 0.9089\n",
      "Epoch 150/1000\n",
      "0s - loss: 0.8977 - val_loss: 0.9062\n",
      "Epoch 151/1000\n",
      "0s - loss: 0.8950 - val_loss: 0.9041\n",
      "Epoch 152/1000\n",
      "0s - loss: 0.8924 - val_loss: 0.9021\n",
      "Epoch 153/1000\n",
      "0s - loss: 0.8900 - val_loss: 0.9001\n",
      "Epoch 154/1000\n",
      "0s - loss: 0.8879 - val_loss: 0.8979\n",
      "Epoch 155/1000\n",
      "0s - loss: 0.8859 - val_loss: 0.8952\n",
      "Epoch 156/1000\n",
      "0s - loss: 0.8833 - val_loss: 0.8929\n",
      "Epoch 157/1000\n",
      "0s - loss: 0.8816 - val_loss: 0.8905\n",
      "Epoch 158/1000\n",
      "0s - loss: 0.8788 - val_loss: 0.8888\n",
      "Epoch 159/1000\n",
      "0s - loss: 0.8766 - val_loss: 0.8869\n",
      "Epoch 160/1000\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0s - loss: 0.8745 - val_loss: 0.8846\n",
      "Epoch 161/1000\n",
      "0s - loss: 0.8717 - val_loss: 0.8822\n",
      "Epoch 162/1000\n",
      "0s - loss: 0.8694 - val_loss: 0.8794\n",
      "Epoch 163/1000\n",
      "0s - loss: 0.8676 - val_loss: 0.8771\n",
      "Epoch 164/1000\n",
      "0s - loss: 0.8652 - val_loss: 0.8751\n",
      "Epoch 165/1000\n",
      "0s - loss: 0.8626 - val_loss: 0.8729\n",
      "Epoch 166/1000\n",
      "0s - loss: 0.8604 - val_loss: 0.8705\n",
      "Epoch 167/1000\n",
      "0s - loss: 0.8581 - val_loss: 0.8683\n",
      "Epoch 168/1000\n",
      "0s - loss: 0.8562 - val_loss: 0.8666\n",
      "Epoch 169/1000\n",
      "0s - loss: 0.8536 - val_loss: 0.8639\n",
      "Epoch 170/1000\n",
      "0s - loss: 0.8513 - val_loss: 0.8618\n",
      "Epoch 171/1000\n",
      "0s - loss: 0.8493 - val_loss: 0.8597\n",
      "Epoch 172/1000\n",
      "0s - loss: 0.8467 - val_loss: 0.8572\n",
      "Epoch 173/1000\n",
      "0s - loss: 0.8445 - val_loss: 0.8550\n",
      "Epoch 174/1000\n",
      "0s - loss: 0.8423 - val_loss: 0.8527\n",
      "Epoch 175/1000\n",
      "0s - loss: 0.8401 - val_loss: 0.8507\n",
      "Epoch 176/1000\n",
      "0s - loss: 0.8377 - val_loss: 0.8488\n",
      "Epoch 177/1000\n",
      "0s - loss: 0.8356 - val_loss: 0.8467\n",
      "Epoch 178/1000\n",
      "0s - loss: 0.8335 - val_loss: 0.8442\n",
      "Epoch 179/1000\n",
      "0s - loss: 0.8317 - val_loss: 0.8421\n",
      "Epoch 180/1000\n",
      "0s - loss: 0.8289 - val_loss: 0.8399\n",
      "Epoch 181/1000\n",
      "0s - loss: 0.8275 - val_loss: 0.8376\n",
      "Epoch 182/1000\n",
      "0s - loss: 0.8248 - val_loss: 0.8355\n",
      "Epoch 183/1000\n",
      "0s - loss: 0.8222 - val_loss: 0.8337\n",
      "Epoch 184/1000\n",
      "0s - loss: 0.8203 - val_loss: 0.8323\n",
      "Epoch 185/1000\n",
      "0s - loss: 0.8186 - val_loss: 0.8300\n",
      "Epoch 186/1000\n",
      "0s - loss: 0.8164 - val_loss: 0.8273\n",
      "Epoch 187/1000\n",
      "0s - loss: 0.8140 - val_loss: 0.8251\n",
      "Epoch 188/1000\n",
      "0s - loss: 0.8121 - val_loss: 0.8233\n",
      "Epoch 189/1000\n",
      "0s - loss: 0.8097 - val_loss: 0.8212\n",
      "Epoch 190/1000\n",
      "0s - loss: 0.8076 - val_loss: 0.8191\n",
      "Epoch 191/1000\n",
      "0s - loss: 0.8057 - val_loss: 0.8173\n",
      "Epoch 192/1000\n",
      "0s - loss: 0.8037 - val_loss: 0.8154\n",
      "Epoch 193/1000\n",
      "0s - loss: 0.8019 - val_loss: 0.8133\n",
      "Epoch 194/1000\n",
      "0s - loss: 0.7996 - val_loss: 0.8110\n",
      "Epoch 195/1000\n",
      "0s - loss: 0.7979 - val_loss: 0.8090\n",
      "Epoch 196/1000\n",
      "0s - loss: 0.7955 - val_loss: 0.8072\n",
      "Epoch 197/1000\n",
      "0s - loss: 0.7946 - val_loss: 0.8056\n",
      "Epoch 198/1000\n",
      "0s - loss: 0.7914 - val_loss: 0.8032\n",
      "Epoch 199/1000\n",
      "0s - loss: 0.7895 - val_loss: 0.8012\n",
      "Epoch 200/1000\n",
      "0s - loss: 0.7883 - val_loss: 0.7993\n",
      "Epoch 201/1000\n",
      "0s - loss: 0.7861 - val_loss: 0.7973\n",
      "Epoch 202/1000\n",
      "0s - loss: 0.7848 - val_loss: 0.7957\n",
      "Epoch 203/1000\n",
      "0s - loss: 0.7821 - val_loss: 0.7935\n",
      "Epoch 204/1000\n",
      "0s - loss: 0.7803 - val_loss: 0.7915\n",
      "Epoch 205/1000\n",
      "0s - loss: 0.7781 - val_loss: 0.7896\n",
      "Epoch 206/1000\n",
      "0s - loss: 0.7760 - val_loss: 0.7880\n",
      "Epoch 207/1000\n",
      "0s - loss: 0.7743 - val_loss: 0.7861\n",
      "Epoch 208/1000\n",
      "0s - loss: 0.7725 - val_loss: 0.7842\n",
      "Epoch 209/1000\n",
      "0s - loss: 0.7709 - val_loss: 0.7823\n",
      "Epoch 210/1000\n",
      "0s - loss: 0.7696 - val_loss: 0.7803\n",
      "Epoch 211/1000\n",
      "0s - loss: 0.7668 - val_loss: 0.7785\n",
      "Epoch 212/1000\n",
      "0s - loss: 0.7651 - val_loss: 0.7766\n",
      "Epoch 213/1000\n",
      "0s - loss: 0.7633 - val_loss: 0.7748\n",
      "Epoch 214/1000\n",
      "0s - loss: 0.7615 - val_loss: 0.7730\n",
      "Epoch 215/1000\n",
      "0s - loss: 0.7597 - val_loss: 0.7712\n",
      "Epoch 216/1000\n",
      "0s - loss: 0.7580 - val_loss: 0.7694\n",
      "Epoch 217/1000\n",
      "0s - loss: 0.7568 - val_loss: 0.7676\n",
      "Epoch 218/1000\n",
      "0s - loss: 0.7551 - val_loss: 0.7658\n",
      "Epoch 219/1000\n",
      "0s - loss: 0.7528 - val_loss: 0.7640\n",
      "Epoch 220/1000\n",
      "0s - loss: 0.7511 - val_loss: 0.7622\n",
      "Epoch 221/1000\n",
      "0s - loss: 0.7498 - val_loss: 0.7604\n",
      "Epoch 222/1000\n",
      "0s - loss: 0.7487 - val_loss: 0.7589\n",
      "Epoch 223/1000\n",
      "0s - loss: 0.7462 - val_loss: 0.7571\n",
      "Epoch 224/1000\n",
      "0s - loss: 0.7454 - val_loss: 0.7552\n",
      "Epoch 225/1000\n",
      "0s - loss: 0.7428 - val_loss: 0.7535\n",
      "Epoch 226/1000\n",
      "0s - loss: 0.7410 - val_loss: 0.7519\n",
      "Epoch 227/1000\n",
      "0s - loss: 0.7394 - val_loss: 0.7502\n",
      "Epoch 228/1000\n",
      "0s - loss: 0.7379 - val_loss: 0.7485\n",
      "Epoch 229/1000\n",
      "0s - loss: 0.7362 - val_loss: 0.7469\n",
      "Epoch 230/1000\n",
      "0s - loss: 0.7348 - val_loss: 0.7452\n",
      "Epoch 231/1000\n",
      "0s - loss: 0.7329 - val_loss: 0.7435\n",
      "Epoch 232/1000\n",
      "0s - loss: 0.7313 - val_loss: 0.7418\n",
      "Epoch 233/1000\n",
      "0s - loss: 0.7299 - val_loss: 0.7403\n",
      "Epoch 234/1000\n",
      "0s - loss: 0.7282 - val_loss: 0.7386\n",
      "Epoch 235/1000\n",
      "0s - loss: 0.7266 - val_loss: 0.7370\n",
      "Epoch 236/1000\n",
      "0s - loss: 0.7249 - val_loss: 0.7354\n",
      "Epoch 237/1000\n",
      "0s - loss: 0.7239 - val_loss: 0.7340\n",
      "Epoch 238/1000\n",
      "0s - loss: 0.7221 - val_loss: 0.7324\n",
      "Epoch 239/1000\n",
      "0s - loss: 0.7205 - val_loss: 0.7308\n",
      "Epoch 240/1000\n",
      "0s - loss: 0.7193 - val_loss: 0.7292\n",
      "Epoch 241/1000\n",
      "0s - loss: 0.7174 - val_loss: 0.7276\n",
      "Epoch 242/1000\n",
      "0s - loss: 0.7159 - val_loss: 0.7261\n",
      "Epoch 243/1000\n",
      "0s - loss: 0.7147 - val_loss: 0.7246\n",
      "Epoch 244/1000\n",
      "0s - loss: 0.7129 - val_loss: 0.7229\n",
      "Epoch 245/1000\n",
      "0s - loss: 0.7118 - val_loss: 0.7214\n",
      "Epoch 246/1000\n",
      "0s - loss: 0.7103 - val_loss: 0.7199\n",
      "Epoch 247/1000\n",
      "0s - loss: 0.7087 - val_loss: 0.7184\n",
      "Epoch 248/1000\n",
      "0s - loss: 0.7072 - val_loss: 0.7169\n",
      "Epoch 249/1000\n",
      "0s - loss: 0.7060 - val_loss: 0.7155\n",
      "Epoch 250/1000\n",
      "0s - loss: 0.7044 - val_loss: 0.7139\n",
      "Epoch 251/1000\n",
      "0s - loss: 0.7037 - val_loss: 0.7126\n",
      "Epoch 252/1000\n",
      "0s - loss: 0.7015 - val_loss: 0.7110\n",
      "Epoch 253/1000\n",
      "0s - loss: 0.7001 - val_loss: 0.7095\n",
      "Epoch 254/1000\n",
      "0s - loss: 0.6990 - val_loss: 0.7081\n",
      "Epoch 255/1000\n",
      "0s - loss: 0.6976 - val_loss: 0.7067\n",
      "Epoch 256/1000\n",
      "0s - loss: 0.6967 - val_loss: 0.7052\n",
      "Epoch 257/1000\n",
      "0s - loss: 0.6948 - val_loss: 0.7038\n",
      "Epoch 258/1000\n",
      "0s - loss: 0.6940 - val_loss: 0.7028\n",
      "Epoch 259/1000\n",
      "0s - loss: 0.6928 - val_loss: 0.7012\n",
      "Epoch 260/1000\n",
      "0s - loss: 0.6909 - val_loss: 0.6997\n",
      "Epoch 261/1000\n",
      "0s - loss: 0.6904 - val_loss: 0.6984\n",
      "Epoch 262/1000\n",
      "0s - loss: 0.6885 - val_loss: 0.6970\n",
      "Epoch 263/1000\n",
      "0s - loss: 0.6879 - val_loss: 0.6959\n",
      "Epoch 264/1000\n",
      "0s - loss: 0.6858 - val_loss: 0.6945\n",
      "Epoch 265/1000\n",
      "0s - loss: 0.6845 - val_loss: 0.6929\n",
      "Epoch 266/1000\n",
      "0s - loss: 0.6832 - val_loss: 0.6916\n",
      "Epoch 267/1000\n",
      "0s - loss: 0.6822 - val_loss: 0.6903\n",
      "Epoch 268/1000\n",
      "0s - loss: 0.6809 - val_loss: 0.6889\n",
      "Epoch 269/1000\n",
      "0s - loss: 0.6796 - val_loss: 0.6876\n",
      "Epoch 270/1000\n",
      "0s - loss: 0.6790 - val_loss: 0.6863\n",
      "Epoch 271/1000\n",
      "0s - loss: 0.6770 - val_loss: 0.6853\n",
      "Epoch 272/1000\n",
      "0s - loss: 0.6763 - val_loss: 0.6842\n",
      "Epoch 273/1000\n",
      "0s - loss: 0.6752 - val_loss: 0.6830\n",
      "Epoch 274/1000\n",
      "0s - loss: 0.6744 - val_loss: 0.6813\n",
      "Epoch 275/1000\n",
      "0s - loss: 0.6736 - val_loss: 0.6802\n",
      "Epoch 276/1000\n",
      "0s - loss: 0.6715 - val_loss: 0.6790\n",
      "Epoch 277/1000\n",
      "0s - loss: 0.6704 - val_loss: 0.6784\n",
      "Epoch 278/1000\n",
      "0s - loss: 0.6697 - val_loss: 0.6766\n",
      "Epoch 279/1000\n",
      "0s - loss: 0.6679 - val_loss: 0.6753\n",
      "Epoch 280/1000\n",
      "0s - loss: 0.6668 - val_loss: 0.6741\n",
      "Epoch 281/1000\n",
      "0s - loss: 0.6664 - val_loss: 0.6731\n",
      "Epoch 282/1000\n",
      "0s - loss: 0.6647 - val_loss: 0.6717\n",
      "Epoch 283/1000\n",
      "0s - loss: 0.6637 - val_loss: 0.6705\n",
      "Epoch 284/1000\n",
      "0s - loss: 0.6625 - val_loss: 0.6693\n",
      "Epoch 285/1000\n",
      "0s - loss: 0.6614 - val_loss: 0.6682\n",
      "Epoch 286/1000\n",
      "0s - loss: 0.6602 - val_loss: 0.6672\n",
      "Epoch 287/1000\n",
      "0s - loss: 0.6594 - val_loss: 0.6660\n",
      "Epoch 288/1000\n",
      "0s - loss: 0.6582 - val_loss: 0.6647\n",
      "Epoch 289/1000\n",
      "0s - loss: 0.6575 - val_loss: 0.6639\n",
      "Epoch 290/1000\n",
      "0s - loss: 0.6561 - val_loss: 0.6627\n",
      "Epoch 291/1000\n",
      "0s - loss: 0.6550 - val_loss: 0.6615\n",
      "Epoch 292/1000\n",
      "0s - loss: 0.6543 - val_loss: 0.6603\n",
      "Epoch 293/1000\n",
      "0s - loss: 0.6534 - val_loss: 0.6594\n",
      "Epoch 294/1000\n",
      "0s - loss: 0.6519 - val_loss: 0.6583\n",
      "Epoch 295/1000\n",
      "0s - loss: 0.6509 - val_loss: 0.6571\n",
      "Epoch 296/1000\n",
      "0s - loss: 0.6500 - val_loss: 0.6560\n",
      "Epoch 297/1000\n",
      "0s - loss: 0.6489 - val_loss: 0.6551\n",
      "Epoch 298/1000\n",
      "0s - loss: 0.6481 - val_loss: 0.6539\n",
      "Epoch 299/1000\n",
      "0s - loss: 0.6473 - val_loss: 0.6530\n",
      "Epoch 300/1000\n",
      "0s - loss: 0.6459 - val_loss: 0.6518\n",
      "Epoch 301/1000\n",
      "0s - loss: 0.6449 - val_loss: 0.6506\n",
      "Epoch 302/1000\n",
      "0s - loss: 0.6446 - val_loss: 0.6496\n",
      "Epoch 303/1000\n",
      "0s - loss: 0.6437 - val_loss: 0.6485\n",
      "Epoch 304/1000\n",
      "0s - loss: 0.6431 - val_loss: 0.6478\n",
      "Epoch 305/1000\n",
      "0s - loss: 0.6414 - val_loss: 0.6467\n",
      "Epoch 306/1000\n",
      "0s - loss: 0.6404 - val_loss: 0.6458\n",
      "Epoch 307/1000\n",
      "0s - loss: 0.6394 - val_loss: 0.6446\n",
      "Epoch 308/1000\n",
      "0s - loss: 0.6384 - val_loss: 0.6435\n",
      "Epoch 309/1000\n",
      "0s - loss: 0.6379 - val_loss: 0.6426\n",
      "Epoch 310/1000\n",
      "0s - loss: 0.6366 - val_loss: 0.6416\n",
      "Epoch 311/1000\n",
      "0s - loss: 0.6359 - val_loss: 0.6405\n",
      "Epoch 312/1000\n",
      "0s - loss: 0.6351 - val_loss: 0.6395\n",
      "Epoch 313/1000\n",
      "0s - loss: 0.6347 - val_loss: 0.6385\n",
      "Epoch 314/1000\n",
      "0s - loss: 0.6330 - val_loss: 0.6380\n",
      "Epoch 315/1000\n",
      "0s - loss: 0.6324 - val_loss: 0.6376\n",
      "Epoch 316/1000\n",
      "0s - loss: 0.6316 - val_loss: 0.6364\n",
      "Epoch 317/1000\n",
      "0s - loss: 0.6305 - val_loss: 0.6349\n",
      "Epoch 318/1000\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0s - loss: 0.6301 - val_loss: 0.6338\n",
      "Epoch 319/1000\n",
      "0s - loss: 0.6289 - val_loss: 0.6328\n",
      "Epoch 320/1000\n",
      "0s - loss: 0.6281 - val_loss: 0.6323\n",
      "Epoch 321/1000\n",
      "0s - loss: 0.6271 - val_loss: 0.6315\n",
      "Epoch 322/1000\n",
      "0s - loss: 0.6265 - val_loss: 0.6303\n",
      "Epoch 323/1000\n",
      "0s - loss: 0.6255 - val_loss: 0.6295\n",
      "Epoch 324/1000\n",
      "0s - loss: 0.6250 - val_loss: 0.6284\n",
      "Epoch 325/1000\n",
      "0s - loss: 0.6238 - val_loss: 0.6278\n",
      "Epoch 326/1000\n",
      "0s - loss: 0.6231 - val_loss: 0.6269\n",
      "Epoch 327/1000\n",
      "0s - loss: 0.6221 - val_loss: 0.6260\n",
      "Epoch 328/1000\n",
      "0s - loss: 0.6214 - val_loss: 0.6250\n",
      "Epoch 329/1000\n",
      "0s - loss: 0.6206 - val_loss: 0.6241\n",
      "Epoch 330/1000\n",
      "0s - loss: 0.6199 - val_loss: 0.6236\n",
      "Epoch 331/1000\n",
      "0s - loss: 0.6192 - val_loss: 0.6226\n",
      "Epoch 332/1000\n",
      "0s - loss: 0.6182 - val_loss: 0.6215\n",
      "Epoch 333/1000\n",
      "0s - loss: 0.6174 - val_loss: 0.6206\n",
      "Epoch 334/1000\n",
      "0s - loss: 0.6167 - val_loss: 0.6200\n",
      "Epoch 335/1000\n",
      "0s - loss: 0.6158 - val_loss: 0.6190\n",
      "Epoch 336/1000\n",
      "0s - loss: 0.6151 - val_loss: 0.6182\n",
      "Epoch 337/1000\n",
      "0s - loss: 0.6142 - val_loss: 0.6172\n",
      "Epoch 338/1000\n",
      "0s - loss: 0.6136 - val_loss: 0.6162\n",
      "Epoch 339/1000\n",
      "0s - loss: 0.6128 - val_loss: 0.6157\n",
      "Epoch 340/1000\n",
      "0s - loss: 0.6120 - val_loss: 0.6148\n",
      "Epoch 341/1000\n",
      "0s - loss: 0.6115 - val_loss: 0.6141\n",
      "Epoch 342/1000\n",
      "0s - loss: 0.6105 - val_loss: 0.6131\n",
      "Epoch 343/1000\n",
      "0s - loss: 0.6108 - val_loss: 0.6119\n",
      "Epoch 344/1000\n",
      "0s - loss: 0.6092 - val_loss: 0.6111\n",
      "Epoch 345/1000\n",
      "0s - loss: 0.6083 - val_loss: 0.6106\n",
      "Epoch 346/1000\n",
      "0s - loss: 0.6079 - val_loss: 0.6105\n",
      "Epoch 347/1000\n",
      "0s - loss: 0.6072 - val_loss: 0.6096\n",
      "Epoch 348/1000\n",
      "0s - loss: 0.6063 - val_loss: 0.6085\n",
      "Epoch 349/1000\n",
      "0s - loss: 0.6054 - val_loss: 0.6074\n",
      "Epoch 350/1000\n",
      "0s - loss: 0.6046 - val_loss: 0.6063\n",
      "Epoch 351/1000\n",
      "0s - loss: 0.6042 - val_loss: 0.6054\n",
      "Epoch 352/1000\n",
      "0s - loss: 0.6037 - val_loss: 0.6049\n",
      "Epoch 353/1000\n",
      "0s - loss: 0.6028 - val_loss: 0.6043\n",
      "Epoch 354/1000\n",
      "0s - loss: 0.6018 - val_loss: 0.6032\n",
      "Epoch 355/1000\n",
      "0s - loss: 0.6012 - val_loss: 0.6023\n",
      "Epoch 356/1000\n",
      "0s - loss: 0.6006 - val_loss: 0.6015\n",
      "Epoch 357/1000\n",
      "0s - loss: 0.5997 - val_loss: 0.6009\n",
      "Epoch 358/1000\n",
      "0s - loss: 0.5989 - val_loss: 0.6004\n",
      "Epoch 359/1000\n",
      "0s - loss: 0.5985 - val_loss: 0.5999\n",
      "Epoch 360/1000\n",
      "0s - loss: 0.5979 - val_loss: 0.5990\n",
      "Epoch 361/1000\n",
      "0s - loss: 0.5972 - val_loss: 0.5977\n",
      "Epoch 362/1000\n",
      "0s - loss: 0.5964 - val_loss: 0.5969\n",
      "Epoch 363/1000\n",
      "0s - loss: 0.5958 - val_loss: 0.5963\n",
      "Epoch 364/1000\n",
      "0s - loss: 0.5951 - val_loss: 0.5955\n",
      "Epoch 365/1000\n",
      "0s - loss: 0.5944 - val_loss: 0.5948\n",
      "Epoch 366/1000\n",
      "0s - loss: 0.5939 - val_loss: 0.5939\n",
      "Epoch 367/1000\n",
      "0s - loss: 0.5931 - val_loss: 0.5933\n",
      "Epoch 368/1000\n",
      "0s - loss: 0.5925 - val_loss: 0.5927\n",
      "Epoch 369/1000\n",
      "0s - loss: 0.5917 - val_loss: 0.5917\n",
      "Epoch 370/1000\n",
      "0s - loss: 0.5911 - val_loss: 0.5908\n",
      "Epoch 371/1000\n",
      "0s - loss: 0.5904 - val_loss: 0.5902\n",
      "Epoch 372/1000\n",
      "0s - loss: 0.5898 - val_loss: 0.5894\n",
      "Epoch 373/1000\n",
      "0s - loss: 0.5902 - val_loss: 0.5884\n",
      "Epoch 374/1000\n",
      "0s - loss: 0.5883 - val_loss: 0.5878\n",
      "Epoch 375/1000\n",
      "0s - loss: 0.5878 - val_loss: 0.5877\n",
      "Epoch 376/1000\n",
      "0s - loss: 0.5873 - val_loss: 0.5870\n",
      "Epoch 377/1000\n",
      "0s - loss: 0.5866 - val_loss: 0.5861\n",
      "Epoch 378/1000\n",
      "0s - loss: 0.5859 - val_loss: 0.5850\n",
      "Epoch 379/1000\n",
      "0s - loss: 0.5854 - val_loss: 0.5843\n",
      "Epoch 380/1000\n",
      "0s - loss: 0.5847 - val_loss: 0.5835\n",
      "Epoch 381/1000\n",
      "0s - loss: 0.5839 - val_loss: 0.5828\n",
      "Epoch 382/1000\n",
      "0s - loss: 0.5833 - val_loss: 0.5823\n",
      "Epoch 383/1000\n",
      "0s - loss: 0.5831 - val_loss: 0.5823\n",
      "Epoch 384/1000\n",
      "0s - loss: 0.5823 - val_loss: 0.5814\n",
      "Epoch 385/1000\n",
      "0s - loss: 0.5819 - val_loss: 0.5799\n",
      "Epoch 386/1000\n",
      "0s - loss: 0.5809 - val_loss: 0.5791\n",
      "Epoch 387/1000\n",
      "0s - loss: 0.5804 - val_loss: 0.5786\n",
      "Epoch 388/1000\n",
      "0s - loss: 0.5797 - val_loss: 0.5778\n",
      "Epoch 389/1000\n",
      "0s - loss: 0.5789 - val_loss: 0.5773\n",
      "Epoch 390/1000\n",
      "0s - loss: 0.5783 - val_loss: 0.5766\n",
      "Epoch 391/1000\n",
      "0s - loss: 0.5778 - val_loss: 0.5759\n",
      "Epoch 392/1000\n",
      "0s - loss: 0.5771 - val_loss: 0.5754\n",
      "Epoch 393/1000\n",
      "0s - loss: 0.5766 - val_loss: 0.5749\n",
      "Epoch 394/1000\n",
      "0s - loss: 0.5760 - val_loss: 0.5742\n",
      "Epoch 395/1000\n",
      "0s - loss: 0.5754 - val_loss: 0.5730\n",
      "Epoch 396/1000\n",
      "0s - loss: 0.5749 - val_loss: 0.5721\n",
      "Epoch 397/1000\n",
      "0s - loss: 0.5741 - val_loss: 0.5714\n",
      "Epoch 398/1000\n",
      "0s - loss: 0.5735 - val_loss: 0.5709\n",
      "Epoch 399/1000\n",
      "0s - loss: 0.5731 - val_loss: 0.5700\n",
      "Epoch 400/1000\n",
      "0s - loss: 0.5723 - val_loss: 0.5693\n",
      "Epoch 401/1000\n",
      "0s - loss: 0.5717 - val_loss: 0.5687\n",
      "Epoch 402/1000\n",
      "0s - loss: 0.5711 - val_loss: 0.5680\n",
      "Epoch 403/1000\n",
      "0s - loss: 0.5708 - val_loss: 0.5678\n",
      "Epoch 404/1000\n",
      "0s - loss: 0.5700 - val_loss: 0.5667\n",
      "Epoch 405/1000\n",
      "0s - loss: 0.5699 - val_loss: 0.5658\n",
      "Epoch 406/1000\n",
      "0s - loss: 0.5686 - val_loss: 0.5653\n",
      "Epoch 407/1000\n",
      "0s - loss: 0.5681 - val_loss: 0.5650\n",
      "Epoch 408/1000\n",
      "0s - loss: 0.5675 - val_loss: 0.5644\n",
      "Epoch 409/1000\n",
      "0s - loss: 0.5669 - val_loss: 0.5637\n",
      "Epoch 410/1000\n",
      "0s - loss: 0.5662 - val_loss: 0.5628\n",
      "Epoch 411/1000\n",
      "0s - loss: 0.5656 - val_loss: 0.5618\n",
      "Epoch 412/1000\n",
      "0s - loss: 0.5653 - val_loss: 0.5613\n",
      "Epoch 413/1000\n",
      "0s - loss: 0.5652 - val_loss: 0.5603\n",
      "Epoch 414/1000\n",
      "0s - loss: 0.5644 - val_loss: 0.5599\n",
      "Epoch 415/1000\n",
      "0s - loss: 0.5632 - val_loss: 0.5591\n",
      "Epoch 416/1000\n",
      "0s - loss: 0.5631 - val_loss: 0.5588\n",
      "Epoch 417/1000\n",
      "0s - loss: 0.5619 - val_loss: 0.5577\n",
      "Epoch 418/1000\n",
      "0s - loss: 0.5613 - val_loss: 0.5568\n",
      "Epoch 419/1000\n",
      "0s - loss: 0.5608 - val_loss: 0.5560\n",
      "Epoch 420/1000\n",
      "0s - loss: 0.5601 - val_loss: 0.5553\n",
      "Epoch 421/1000\n",
      "0s - loss: 0.5596 - val_loss: 0.5547\n",
      "Epoch 422/1000\n",
      "0s - loss: 0.5589 - val_loss: 0.5542\n",
      "Epoch 423/1000\n",
      "0s - loss: 0.5584 - val_loss: 0.5535\n",
      "Epoch 424/1000\n",
      "0s - loss: 0.5576 - val_loss: 0.5530\n",
      "Epoch 425/1000\n",
      "0s - loss: 0.5570 - val_loss: 0.5521\n",
      "Epoch 426/1000\n",
      "0s - loss: 0.5565 - val_loss: 0.5512\n",
      "Epoch 427/1000\n",
      "0s - loss: 0.5559 - val_loss: 0.5504\n",
      "Epoch 428/1000\n",
      "0s - loss: 0.5551 - val_loss: 0.5500\n",
      "Epoch 429/1000\n",
      "0s - loss: 0.5547 - val_loss: 0.5498\n",
      "Epoch 430/1000\n",
      "0s - loss: 0.5540 - val_loss: 0.5485\n",
      "Epoch 431/1000\n",
      "0s - loss: 0.5537 - val_loss: 0.5474\n",
      "Epoch 432/1000\n",
      "0s - loss: 0.5526 - val_loss: 0.5467\n",
      "Epoch 433/1000\n",
      "0s - loss: 0.5518 - val_loss: 0.5462\n",
      "Epoch 434/1000\n",
      "0s - loss: 0.5517 - val_loss: 0.5459\n",
      "Epoch 435/1000\n",
      "0s - loss: 0.5508 - val_loss: 0.5449\n",
      "Epoch 436/1000\n",
      "0s - loss: 0.5501 - val_loss: 0.5440\n",
      "Epoch 437/1000\n",
      "0s - loss: 0.5494 - val_loss: 0.5432\n",
      "Epoch 438/1000\n",
      "0s - loss: 0.5487 - val_loss: 0.5424\n",
      "Epoch 439/1000\n",
      "0s - loss: 0.5482 - val_loss: 0.5417\n",
      "Epoch 440/1000\n",
      "0s - loss: 0.5474 - val_loss: 0.5410\n",
      "Epoch 441/1000\n",
      "0s - loss: 0.5470 - val_loss: 0.5406\n",
      "Epoch 442/1000\n",
      "0s - loss: 0.5463 - val_loss: 0.5399\n",
      "Epoch 443/1000\n",
      "0s - loss: 0.5459 - val_loss: 0.5388\n",
      "Epoch 444/1000\n",
      "0s - loss: 0.5452 - val_loss: 0.5384\n",
      "Epoch 445/1000\n",
      "0s - loss: 0.5440 - val_loss: 0.5374\n",
      "Epoch 446/1000\n",
      "0s - loss: 0.5434 - val_loss: 0.5365\n",
      "Epoch 447/1000\n",
      "0s - loss: 0.5428 - val_loss: 0.5358\n",
      "Epoch 448/1000\n",
      "0s - loss: 0.5424 - val_loss: 0.5350\n",
      "Epoch 449/1000\n",
      "0s - loss: 0.5415 - val_loss: 0.5344\n",
      "Epoch 450/1000\n",
      "0s - loss: 0.5408 - val_loss: 0.5337\n",
      "Epoch 451/1000\n",
      "0s - loss: 0.5403 - val_loss: 0.5333\n",
      "Epoch 452/1000\n",
      "0s - loss: 0.5395 - val_loss: 0.5325\n",
      "Epoch 453/1000\n",
      "0s - loss: 0.5390 - val_loss: 0.5313\n",
      "Epoch 454/1000\n",
      "0s - loss: 0.5380 - val_loss: 0.5305\n",
      "Epoch 455/1000\n",
      "0s - loss: 0.5372 - val_loss: 0.5298\n",
      "Epoch 456/1000\n",
      "0s - loss: 0.5364 - val_loss: 0.5290\n",
      "Epoch 457/1000\n",
      "0s - loss: 0.5356 - val_loss: 0.5282\n",
      "Epoch 458/1000\n",
      "0s - loss: 0.5351 - val_loss: 0.5275\n",
      "Epoch 459/1000\n",
      "0s - loss: 0.5343 - val_loss: 0.5265\n",
      "Epoch 460/1000\n",
      "0s - loss: 0.5330 - val_loss: 0.5257\n",
      "Epoch 461/1000\n",
      "0s - loss: 0.5321 - val_loss: 0.5249\n",
      "Epoch 462/1000\n",
      "0s - loss: 0.5318 - val_loss: 0.5243\n",
      "Epoch 463/1000\n",
      "0s - loss: 0.5303 - val_loss: 0.5231\n",
      "Epoch 464/1000\n",
      "0s - loss: 0.5296 - val_loss: 0.5223\n",
      "Epoch 465/1000\n",
      "0s - loss: 0.5290 - val_loss: 0.5215\n",
      "Epoch 466/1000\n",
      "0s - loss: 0.5280 - val_loss: 0.5206\n",
      "Epoch 467/1000\n",
      "0s - loss: 0.5269 - val_loss: 0.5199\n",
      "Epoch 468/1000\n",
      "0s - loss: 0.5256 - val_loss: 0.5189\n",
      "Epoch 469/1000\n",
      "0s - loss: 0.5251 - val_loss: 0.5178\n",
      "Epoch 470/1000\n",
      "0s - loss: 0.5235 - val_loss: 0.5167\n",
      "Epoch 471/1000\n",
      "0s - loss: 0.5225 - val_loss: 0.5156\n",
      "Epoch 472/1000\n",
      "0s - loss: 0.5213 - val_loss: 0.5145\n",
      "Epoch 473/1000\n",
      "0s - loss: 0.5204 - val_loss: 0.5134\n",
      "Epoch 474/1000\n",
      "0s - loss: 0.5192 - val_loss: 0.5123\n",
      "Epoch 475/1000\n",
      "0s - loss: 0.5178 - val_loss: 0.5111\n",
      "Epoch 476/1000\n",
      "0s - loss: 0.5167 - val_loss: 0.5100\n",
      "Epoch 477/1000\n",
      "0s - loss: 0.5159 - val_loss: 0.5090\n",
      "Epoch 478/1000\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0s - loss: 0.5140 - val_loss: 0.5077\n",
      "Epoch 479/1000\n",
      "0s - loss: 0.5129 - val_loss: 0.5065\n",
      "Epoch 480/1000\n",
      "0s - loss: 0.5118 - val_loss: 0.5052\n",
      "Epoch 481/1000\n",
      "0s - loss: 0.5112 - val_loss: 0.5038\n",
      "Epoch 482/1000\n",
      "0s - loss: 0.5095 - val_loss: 0.5026\n",
      "Epoch 483/1000\n",
      "0s - loss: 0.5076 - val_loss: 0.5010\n",
      "Epoch 484/1000\n",
      "0s - loss: 0.5063 - val_loss: 0.4997\n",
      "Epoch 485/1000\n",
      "0s - loss: 0.5049 - val_loss: 0.4984\n",
      "Epoch 486/1000\n",
      "0s - loss: 0.5032 - val_loss: 0.4969\n",
      "Epoch 487/1000\n",
      "0s - loss: 0.5019 - val_loss: 0.4955\n",
      "Epoch 488/1000\n",
      "0s - loss: 0.5003 - val_loss: 0.4940\n",
      "Epoch 489/1000\n",
      "0s - loss: 0.4990 - val_loss: 0.4926\n",
      "Epoch 490/1000\n",
      "0s - loss: 0.4972 - val_loss: 0.4911\n",
      "Epoch 491/1000\n",
      "0s - loss: 0.4959 - val_loss: 0.4897\n",
      "Epoch 492/1000\n",
      "0s - loss: 0.4946 - val_loss: 0.4882\n",
      "Epoch 493/1000\n",
      "0s - loss: 0.4932 - val_loss: 0.4868\n",
      "Epoch 494/1000\n",
      "0s - loss: 0.4911 - val_loss: 0.4852\n",
      "Epoch 495/1000\n",
      "0s - loss: 0.4896 - val_loss: 0.4837\n",
      "Epoch 496/1000\n",
      "0s - loss: 0.4884 - val_loss: 0.4823\n",
      "Epoch 497/1000\n",
      "0s - loss: 0.4868 - val_loss: 0.4807\n",
      "Epoch 498/1000\n",
      "0s - loss: 0.4852 - val_loss: 0.4791\n",
      "Epoch 499/1000\n",
      "0s - loss: 0.4834 - val_loss: 0.4777\n",
      "Epoch 500/1000\n",
      "0s - loss: 0.4820 - val_loss: 0.4765\n",
      "Epoch 501/1000\n",
      "0s - loss: 0.4806 - val_loss: 0.4747\n",
      "Epoch 502/1000\n",
      "0s - loss: 0.4787 - val_loss: 0.4731\n",
      "Epoch 503/1000\n",
      "0s - loss: 0.4769 - val_loss: 0.4714\n",
      "Epoch 504/1000\n",
      "0s - loss: 0.4758 - val_loss: 0.4699\n",
      "Epoch 505/1000\n",
      "0s - loss: 0.4741 - val_loss: 0.4684\n",
      "Epoch 506/1000\n",
      "0s - loss: 0.4723 - val_loss: 0.4667\n",
      "Epoch 507/1000\n",
      "0s - loss: 0.4707 - val_loss: 0.4647\n",
      "Epoch 508/1000\n",
      "0s - loss: 0.4691 - val_loss: 0.4630\n",
      "Epoch 509/1000\n",
      "0s - loss: 0.4676 - val_loss: 0.4613\n",
      "Epoch 510/1000\n",
      "0s - loss: 0.4656 - val_loss: 0.4594\n",
      "Epoch 511/1000\n",
      "0s - loss: 0.4640 - val_loss: 0.4576\n",
      "Epoch 512/1000\n",
      "0s - loss: 0.4624 - val_loss: 0.4557\n",
      "Epoch 513/1000\n",
      "0s - loss: 0.4607 - val_loss: 0.4539\n",
      "Epoch 514/1000\n",
      "0s - loss: 0.4587 - val_loss: 0.4521\n",
      "Epoch 515/1000\n",
      "0s - loss: 0.4570 - val_loss: 0.4502\n",
      "Epoch 516/1000\n",
      "0s - loss: 0.4553 - val_loss: 0.4484\n",
      "Epoch 517/1000\n",
      "0s - loss: 0.4536 - val_loss: 0.4466\n",
      "Epoch 518/1000\n",
      "0s - loss: 0.4514 - val_loss: 0.4447\n",
      "Epoch 519/1000\n",
      "0s - loss: 0.4497 - val_loss: 0.4430\n",
      "Epoch 520/1000\n",
      "0s - loss: 0.4478 - val_loss: 0.4412\n",
      "Epoch 521/1000\n",
      "0s - loss: 0.4463 - val_loss: 0.4393\n",
      "Epoch 522/1000\n",
      "0s - loss: 0.4443 - val_loss: 0.4375\n",
      "Epoch 523/1000\n",
      "0s - loss: 0.4426 - val_loss: 0.4357\n",
      "Epoch 524/1000\n",
      "0s - loss: 0.4409 - val_loss: 0.4340\n",
      "Epoch 525/1000\n",
      "0s - loss: 0.4392 - val_loss: 0.4322\n",
      "Epoch 526/1000\n",
      "0s - loss: 0.4376 - val_loss: 0.4305\n",
      "Epoch 527/1000\n",
      "0s - loss: 0.4358 - val_loss: 0.4287\n",
      "Epoch 528/1000\n",
      "0s - loss: 0.4347 - val_loss: 0.4271\n",
      "Epoch 529/1000\n",
      "0s - loss: 0.4328 - val_loss: 0.4253\n",
      "Epoch 530/1000\n",
      "0s - loss: 0.4309 - val_loss: 0.4236\n",
      "Epoch 531/1000\n",
      "0s - loss: 0.4295 - val_loss: 0.4220\n",
      "Epoch 532/1000\n",
      "0s - loss: 0.4281 - val_loss: 0.4203\n",
      "Epoch 533/1000\n",
      "0s - loss: 0.4264 - val_loss: 0.4186\n",
      "Epoch 534/1000\n",
      "0s - loss: 0.4249 - val_loss: 0.4169\n",
      "Epoch 535/1000\n",
      "0s - loss: 0.4234 - val_loss: 0.4153\n",
      "Epoch 536/1000\n",
      "0s - loss: 0.4220 - val_loss: 0.4137\n",
      "Epoch 537/1000\n",
      "0s - loss: 0.4207 - val_loss: 0.4122\n",
      "Epoch 538/1000\n",
      "0s - loss: 0.4192 - val_loss: 0.4105\n",
      "Epoch 539/1000\n",
      "0s - loss: 0.4178 - val_loss: 0.4089\n",
      "Epoch 540/1000\n",
      "0s - loss: 0.4164 - val_loss: 0.4074\n",
      "Epoch 541/1000\n",
      "0s - loss: 0.4152 - val_loss: 0.4059\n",
      "Epoch 542/1000\n",
      "0s - loss: 0.4136 - val_loss: 0.4044\n",
      "Epoch 543/1000\n",
      "0s - loss: 0.4122 - val_loss: 0.4029\n",
      "Epoch 544/1000\n",
      "0s - loss: 0.4109 - val_loss: 0.4014\n",
      "Epoch 545/1000\n",
      "0s - loss: 0.4095 - val_loss: 0.3999\n",
      "Epoch 546/1000\n",
      "0s - loss: 0.4082 - val_loss: 0.3985\n",
      "Epoch 547/1000\n",
      "0s - loss: 0.4068 - val_loss: 0.3971\n",
      "Epoch 548/1000\n",
      "0s - loss: 0.4058 - val_loss: 0.3956\n",
      "Epoch 549/1000\n",
      "0s - loss: 0.4040 - val_loss: 0.3942\n",
      "Epoch 550/1000\n",
      "0s - loss: 0.4028 - val_loss: 0.3928\n",
      "Epoch 551/1000\n",
      "0s - loss: 0.4015 - val_loss: 0.3915\n",
      "Epoch 552/1000\n",
      "0s - loss: 0.4003 - val_loss: 0.3899\n",
      "Epoch 553/1000\n",
      "0s - loss: 0.3987 - val_loss: 0.3886\n",
      "Epoch 554/1000\n",
      "0s - loss: 0.3977 - val_loss: 0.3874\n",
      "Epoch 555/1000\n",
      "0s - loss: 0.3963 - val_loss: 0.3858\n",
      "Epoch 556/1000\n",
      "0s - loss: 0.3950 - val_loss: 0.3844\n",
      "Epoch 557/1000\n",
      "0s - loss: 0.3939 - val_loss: 0.3829\n",
      "Epoch 558/1000\n",
      "0s - loss: 0.3927 - val_loss: 0.3816\n",
      "Epoch 559/1000\n",
      "0s - loss: 0.3914 - val_loss: 0.3803\n",
      "Epoch 560/1000\n",
      "0s - loss: 0.3902 - val_loss: 0.3791\n",
      "Epoch 561/1000\n",
      "0s - loss: 0.3890 - val_loss: 0.3780\n",
      "Epoch 562/1000\n",
      "0s - loss: 0.3885 - val_loss: 0.3771\n",
      "Epoch 563/1000\n",
      "0s - loss: 0.3863 - val_loss: 0.3754\n",
      "Epoch 564/1000\n",
      "0s - loss: 0.3854 - val_loss: 0.3740\n",
      "Epoch 565/1000\n",
      "0s - loss: 0.3845 - val_loss: 0.3730\n",
      "Epoch 566/1000\n",
      "0s - loss: 0.3827 - val_loss: 0.3718\n",
      "Epoch 567/1000\n",
      "0s - loss: 0.3815 - val_loss: 0.3705\n",
      "Epoch 568/1000\n",
      "0s - loss: 0.3806 - val_loss: 0.3695\n",
      "Epoch 569/1000\n",
      "0s - loss: 0.3794 - val_loss: 0.3682\n",
      "Epoch 570/1000\n",
      "0s - loss: 0.3780 - val_loss: 0.3669\n",
      "Epoch 571/1000\n",
      "0s - loss: 0.3774 - val_loss: 0.3655\n",
      "Epoch 572/1000\n",
      "0s - loss: 0.3760 - val_loss: 0.3644\n",
      "Epoch 573/1000\n",
      "0s - loss: 0.3749 - val_loss: 0.3632\n",
      "Epoch 574/1000\n",
      "0s - loss: 0.3740 - val_loss: 0.3621\n",
      "Epoch 575/1000\n",
      "0s - loss: 0.3738 - val_loss: 0.3615\n",
      "Epoch 576/1000\n",
      "0s - loss: 0.3725 - val_loss: 0.3606\n",
      "Epoch 577/1000\n",
      "0s - loss: 0.3711 - val_loss: 0.3588\n",
      "Epoch 578/1000\n",
      "0s - loss: 0.3696 - val_loss: 0.3576\n",
      "Epoch 579/1000\n",
      "0s - loss: 0.3686 - val_loss: 0.3567\n",
      "Epoch 580/1000\n",
      "0s - loss: 0.3677 - val_loss: 0.3555\n",
      "Epoch 581/1000\n",
      "0s - loss: 0.3665 - val_loss: 0.3547\n",
      "Epoch 582/1000\n",
      "0s - loss: 0.3655 - val_loss: 0.3537\n",
      "Epoch 583/1000\n",
      "0s - loss: 0.3647 - val_loss: 0.3529\n",
      "Epoch 584/1000\n",
      "0s - loss: 0.3639 - val_loss: 0.3522\n",
      "Epoch 585/1000\n",
      "0s - loss: 0.3629 - val_loss: 0.3507\n",
      "Epoch 586/1000\n",
      "0s - loss: 0.3618 - val_loss: 0.3496\n",
      "Epoch 587/1000\n",
      "0s - loss: 0.3607 - val_loss: 0.3483\n",
      "Epoch 588/1000\n",
      "0s - loss: 0.3599 - val_loss: 0.3472\n",
      "Epoch 589/1000\n",
      "0s - loss: 0.3594 - val_loss: 0.3460\n",
      "Epoch 590/1000\n",
      "0s - loss: 0.3584 - val_loss: 0.3451\n",
      "Epoch 591/1000\n",
      "0s - loss: 0.3574 - val_loss: 0.3443\n",
      "Epoch 592/1000\n",
      "0s - loss: 0.3563 - val_loss: 0.3435\n",
      "Epoch 593/1000\n",
      "0s - loss: 0.3554 - val_loss: 0.3425\n",
      "Epoch 594/1000\n",
      "0s - loss: 0.3547 - val_loss: 0.3417\n",
      "Epoch 595/1000\n",
      "0s - loss: 0.3537 - val_loss: 0.3406\n",
      "Epoch 596/1000\n",
      "0s - loss: 0.3528 - val_loss: 0.3394\n",
      "Epoch 597/1000\n",
      "0s - loss: 0.3523 - val_loss: 0.3383\n",
      "Epoch 598/1000\n",
      "0s - loss: 0.3515 - val_loss: 0.3373\n",
      "Epoch 599/1000\n",
      "0s - loss: 0.3502 - val_loss: 0.3367\n",
      "Epoch 600/1000\n",
      "0s - loss: 0.3495 - val_loss: 0.3360\n",
      "Epoch 601/1000\n",
      "0s - loss: 0.3487 - val_loss: 0.3352\n",
      "Epoch 602/1000\n",
      "0s - loss: 0.3480 - val_loss: 0.3346\n",
      "Epoch 603/1000\n",
      "0s - loss: 0.3472 - val_loss: 0.3338\n",
      "Epoch 604/1000\n",
      "0s - loss: 0.3463 - val_loss: 0.3327\n",
      "Epoch 605/1000\n",
      "0s - loss: 0.3459 - val_loss: 0.3312\n",
      "Epoch 606/1000\n",
      "0s - loss: 0.3447 - val_loss: 0.3302\n",
      "Epoch 607/1000\n",
      "0s - loss: 0.3442 - val_loss: 0.3294\n",
      "Epoch 608/1000\n",
      "0s - loss: 0.3431 - val_loss: 0.3285\n",
      "Epoch 609/1000\n",
      "0s - loss: 0.3430 - val_loss: 0.3281\n",
      "Epoch 610/1000\n",
      "0s - loss: 0.3415 - val_loss: 0.3274\n",
      "Epoch 611/1000\n",
      "0s - loss: 0.3410 - val_loss: 0.3266\n",
      "Epoch 612/1000\n",
      "0s - loss: 0.3401 - val_loss: 0.3252\n",
      "Epoch 613/1000\n",
      "0s - loss: 0.3392 - val_loss: 0.3242\n",
      "Epoch 614/1000\n",
      "0s - loss: 0.3386 - val_loss: 0.3234\n",
      "Epoch 615/1000\n",
      "0s - loss: 0.3379 - val_loss: 0.3230\n",
      "Epoch 616/1000\n",
      "0s - loss: 0.3370 - val_loss: 0.3222\n",
      "Epoch 617/1000\n",
      "0s - loss: 0.3363 - val_loss: 0.3214\n",
      "Epoch 618/1000\n",
      "0s - loss: 0.3355 - val_loss: 0.3205\n",
      "Epoch 619/1000\n",
      "0s - loss: 0.3348 - val_loss: 0.3194\n",
      "Epoch 620/1000\n",
      "0s - loss: 0.3342 - val_loss: 0.3184\n",
      "Epoch 621/1000\n",
      "0s - loss: 0.3333 - val_loss: 0.3176\n",
      "Epoch 622/1000\n",
      "0s - loss: 0.3326 - val_loss: 0.3169\n",
      "Epoch 623/1000\n",
      "0s - loss: 0.3326 - val_loss: 0.3165\n",
      "Epoch 624/1000\n",
      "0s - loss: 0.3312 - val_loss: 0.3153\n",
      "Epoch 625/1000\n",
      "0s - loss: 0.3306 - val_loss: 0.3144\n",
      "Epoch 626/1000\n",
      "0s - loss: 0.3302 - val_loss: 0.3136\n",
      "Epoch 627/1000\n",
      "0s - loss: 0.3292 - val_loss: 0.3127\n",
      "Epoch 628/1000\n",
      "0s - loss: 0.3286 - val_loss: 0.3120\n",
      "Epoch 629/1000\n",
      "0s - loss: 0.3280 - val_loss: 0.3116\n",
      "Epoch 630/1000\n",
      "0s - loss: 0.3273 - val_loss: 0.3107\n",
      "Epoch 631/1000\n",
      "0s - loss: 0.3265 - val_loss: 0.3101\n",
      "Epoch 632/1000\n",
      "0s - loss: 0.3259 - val_loss: 0.3093\n",
      "Epoch 633/1000\n",
      "0s - loss: 0.3251 - val_loss: 0.3088\n",
      "Epoch 634/1000\n",
      "0s - loss: 0.3245 - val_loss: 0.3082\n",
      "Epoch 635/1000\n",
      "0s - loss: 0.3240 - val_loss: 0.3072\n",
      "Epoch 636/1000\n",
      "0s - loss: 0.3232 - val_loss: 0.3064\n",
      "Epoch 637/1000\n",
      "0s - loss: 0.3228 - val_loss: 0.3058\n",
      "Epoch 638/1000\n",
      "0s - loss: 0.3219 - val_loss: 0.3050\n",
      "Epoch 639/1000\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0s - loss: 0.3215 - val_loss: 0.3041\n",
      "Epoch 640/1000\n",
      "0s - loss: 0.3205 - val_loss: 0.3034\n",
      "Epoch 641/1000\n",
      "0s - loss: 0.3204 - val_loss: 0.3032\n",
      "Epoch 642/1000\n",
      "0s - loss: 0.3195 - val_loss: 0.3022\n",
      "Epoch 643/1000\n",
      "0s - loss: 0.3190 - val_loss: 0.3012\n",
      "Epoch 644/1000\n",
      "0s - loss: 0.3183 - val_loss: 0.3004\n",
      "Epoch 645/1000\n",
      "0s - loss: 0.3175 - val_loss: 0.3000\n",
      "Epoch 646/1000\n",
      "0s - loss: 0.3170 - val_loss: 0.2993\n",
      "Epoch 647/1000\n",
      "0s - loss: 0.3163 - val_loss: 0.2988\n",
      "Epoch 648/1000\n",
      "0s - loss: 0.3156 - val_loss: 0.2981\n",
      "Epoch 649/1000\n",
      "0s - loss: 0.3151 - val_loss: 0.2973\n",
      "Epoch 650/1000\n",
      "0s - loss: 0.3144 - val_loss: 0.2964\n",
      "Epoch 651/1000\n",
      "0s - loss: 0.3139 - val_loss: 0.2956\n",
      "Epoch 652/1000\n",
      "0s - loss: 0.3133 - val_loss: 0.2950\n",
      "Epoch 653/1000\n",
      "0s - loss: 0.3126 - val_loss: 0.2944\n",
      "Epoch 654/1000\n",
      "0s - loss: 0.3120 - val_loss: 0.2941\n",
      "Epoch 655/1000\n",
      "0s - loss: 0.3119 - val_loss: 0.2938\n",
      "Epoch 656/1000\n",
      "0s - loss: 0.3112 - val_loss: 0.2926\n",
      "Epoch 657/1000\n",
      "0s - loss: 0.3104 - val_loss: 0.2918\n",
      "Epoch 658/1000\n",
      "0s - loss: 0.3098 - val_loss: 0.2912\n",
      "Epoch 659/1000\n",
      "0s - loss: 0.3093 - val_loss: 0.2908\n",
      "Epoch 660/1000\n",
      "0s - loss: 0.3086 - val_loss: 0.2902\n",
      "Epoch 661/1000\n",
      "0s - loss: 0.3080 - val_loss: 0.2893\n",
      "Epoch 662/1000\n",
      "0s - loss: 0.3076 - val_loss: 0.2884\n",
      "Epoch 663/1000\n",
      "0s - loss: 0.3070 - val_loss: 0.2878\n",
      "Epoch 664/1000\n",
      "0s - loss: 0.3065 - val_loss: 0.2871\n",
      "Epoch 665/1000\n",
      "0s - loss: 0.3059 - val_loss: 0.2868\n",
      "Epoch 666/1000\n",
      "0s - loss: 0.3053 - val_loss: 0.2866\n",
      "Epoch 667/1000\n",
      "0s - loss: 0.3048 - val_loss: 0.2858\n",
      "Epoch 668/1000\n",
      "0s - loss: 0.3043 - val_loss: 0.2851\n",
      "Epoch 669/1000\n",
      "0s - loss: 0.3037 - val_loss: 0.2845\n",
      "Epoch 670/1000\n",
      "0s - loss: 0.3030 - val_loss: 0.2839\n",
      "Epoch 671/1000\n",
      "0s - loss: 0.3027 - val_loss: 0.2832\n",
      "Epoch 672/1000\n",
      "0s - loss: 0.3021 - val_loss: 0.2826\n",
      "Epoch 673/1000\n",
      "0s - loss: 0.3015 - val_loss: 0.2820\n",
      "Epoch 674/1000\n",
      "0s - loss: 0.3009 - val_loss: 0.2816\n",
      "Epoch 675/1000\n",
      "0s - loss: 0.3007 - val_loss: 0.2815\n",
      "Epoch 676/1000\n",
      "0s - loss: 0.3000 - val_loss: 0.2808\n",
      "Epoch 677/1000\n",
      "0s - loss: 0.2994 - val_loss: 0.2799\n",
      "Epoch 678/1000\n",
      "0s - loss: 0.2987 - val_loss: 0.2789\n",
      "Epoch 679/1000\n",
      "0s - loss: 0.2983 - val_loss: 0.2781\n",
      "Epoch 680/1000\n",
      "0s - loss: 0.2980 - val_loss: 0.2776\n",
      "Epoch 681/1000\n",
      "0s - loss: 0.2973 - val_loss: 0.2772\n",
      "Epoch 682/1000\n",
      "0s - loss: 0.2969 - val_loss: 0.2771\n",
      "Epoch 683/1000\n",
      "0s - loss: 0.2964 - val_loss: 0.2770\n",
      "Epoch 684/1000\n",
      "0s - loss: 0.2959 - val_loss: 0.2762\n",
      "Epoch 685/1000\n",
      "0s - loss: 0.2953 - val_loss: 0.2753\n",
      "Epoch 686/1000\n",
      "0s - loss: 0.2948 - val_loss: 0.2745\n",
      "Epoch 687/1000\n",
      "0s - loss: 0.2942 - val_loss: 0.2737\n",
      "Epoch 688/1000\n",
      "0s - loss: 0.2941 - val_loss: 0.2733\n",
      "Epoch 689/1000\n",
      "0s - loss: 0.2933 - val_loss: 0.2727\n",
      "Epoch 690/1000\n",
      "0s - loss: 0.2931 - val_loss: 0.2720\n",
      "Epoch 691/1000\n",
      "0s - loss: 0.2925 - val_loss: 0.2718\n",
      "Epoch 692/1000\n",
      "0s - loss: 0.2919 - val_loss: 0.2714\n",
      "Epoch 693/1000\n",
      "0s - loss: 0.2916 - val_loss: 0.2711\n",
      "Epoch 694/1000\n",
      "0s - loss: 0.2909 - val_loss: 0.2704\n",
      "Epoch 695/1000\n",
      "0s - loss: 0.2904 - val_loss: 0.2696\n",
      "Epoch 696/1000\n",
      "0s - loss: 0.2900 - val_loss: 0.2692\n",
      "Epoch 697/1000\n",
      "0s - loss: 0.2895 - val_loss: 0.2688\n",
      "Epoch 698/1000\n",
      "0s - loss: 0.2895 - val_loss: 0.2678\n",
      "Epoch 699/1000\n",
      "0s - loss: 0.2885 - val_loss: 0.2677\n",
      "Epoch 700/1000\n",
      "0s - loss: 0.2882 - val_loss: 0.2677\n",
      "Epoch 701/1000\n",
      "0s - loss: 0.2876 - val_loss: 0.2671\n",
      "Epoch 702/1000\n",
      "0s - loss: 0.2873 - val_loss: 0.2662\n",
      "Epoch 703/1000\n",
      "0s - loss: 0.2866 - val_loss: 0.2655\n",
      "Epoch 704/1000\n",
      "0s - loss: 0.2861 - val_loss: 0.2649\n",
      "Epoch 705/1000\n",
      "0s - loss: 0.2856 - val_loss: 0.2644\n",
      "Epoch 706/1000\n",
      "0s - loss: 0.2851 - val_loss: 0.2641\n",
      "Epoch 707/1000\n",
      "0s - loss: 0.2848 - val_loss: 0.2636\n",
      "Epoch 708/1000\n",
      "0s - loss: 0.2845 - val_loss: 0.2634\n",
      "Epoch 709/1000\n",
      "0s - loss: 0.2838 - val_loss: 0.2628\n",
      "Epoch 710/1000\n",
      "0s - loss: 0.2834 - val_loss: 0.2624\n",
      "Epoch 711/1000\n",
      "0s - loss: 0.2830 - val_loss: 0.2616\n",
      "Epoch 712/1000\n",
      "0s - loss: 0.2830 - val_loss: 0.2606\n",
      "Epoch 713/1000\n",
      "0s - loss: 0.2820 - val_loss: 0.2602\n",
      "Epoch 714/1000\n",
      "0s - loss: 0.2815 - val_loss: 0.2599\n",
      "Epoch 715/1000\n",
      "0s - loss: 0.2811 - val_loss: 0.2596\n",
      "Epoch 716/1000\n",
      "0s - loss: 0.2817 - val_loss: 0.2599\n",
      "Epoch 717/1000\n",
      "0s - loss: 0.2803 - val_loss: 0.2591\n",
      "Epoch 718/1000\n",
      "0s - loss: 0.2799 - val_loss: 0.2579\n",
      "Epoch 719/1000\n",
      "0s - loss: 0.2799 - val_loss: 0.2570\n",
      "Epoch 720/1000\n",
      "0s - loss: 0.2792 - val_loss: 0.2570\n",
      "Epoch 721/1000\n",
      "0s - loss: 0.2786 - val_loss: 0.2569\n",
      "Epoch 722/1000\n",
      "0s - loss: 0.2783 - val_loss: 0.2560\n",
      "Epoch 723/1000\n",
      "0s - loss: 0.2776 - val_loss: 0.2556\n",
      "Epoch 724/1000\n",
      "0s - loss: 0.2770 - val_loss: 0.2554\n",
      "Epoch 725/1000\n",
      "0s - loss: 0.2768 - val_loss: 0.2550\n",
      "Epoch 726/1000\n",
      "0s - loss: 0.2763 - val_loss: 0.2550\n",
      "Epoch 727/1000\n",
      "0s - loss: 0.2761 - val_loss: 0.2548\n",
      "Epoch 728/1000\n",
      "0s - loss: 0.2755 - val_loss: 0.2540\n",
      "Epoch 729/1000\n",
      "0s - loss: 0.2751 - val_loss: 0.2531\n",
      "Epoch 730/1000\n",
      "0s - loss: 0.2745 - val_loss: 0.2526\n",
      "Epoch 731/1000\n",
      "0s - loss: 0.2741 - val_loss: 0.2523\n",
      "Epoch 732/1000\n",
      "0s - loss: 0.2739 - val_loss: 0.2518\n",
      "Epoch 733/1000\n",
      "0s - loss: 0.2732 - val_loss: 0.2515\n",
      "Epoch 734/1000\n",
      "0s - loss: 0.2730 - val_loss: 0.2514\n",
      "Epoch 735/1000\n",
      "0s - loss: 0.2725 - val_loss: 0.2507\n",
      "Epoch 736/1000\n",
      "0s - loss: 0.2720 - val_loss: 0.2502\n",
      "Epoch 737/1000\n",
      "0s - loss: 0.2718 - val_loss: 0.2495\n",
      "Epoch 738/1000\n",
      "0s - loss: 0.2718 - val_loss: 0.2488\n",
      "Epoch 739/1000\n",
      "0s - loss: 0.2706 - val_loss: 0.2489\n",
      "Epoch 740/1000\n",
      "0s - loss: 0.2711 - val_loss: 0.2494\n",
      "Epoch 741/1000\n",
      "0s - loss: 0.2704 - val_loss: 0.2488\n",
      "Epoch 742/1000\n",
      "0s - loss: 0.2699 - val_loss: 0.2474\n",
      "Epoch 743/1000\n",
      "0s - loss: 0.2692 - val_loss: 0.2469\n",
      "Epoch 744/1000\n",
      "0s - loss: 0.2688 - val_loss: 0.2464\n",
      "Epoch 745/1000\n",
      "0s - loss: 0.2685 - val_loss: 0.2458\n",
      "Epoch 746/1000\n",
      "0s - loss: 0.2681 - val_loss: 0.2456\n",
      "Epoch 747/1000\n",
      "0s - loss: 0.2675 - val_loss: 0.2452\n",
      "Epoch 748/1000\n",
      "0s - loss: 0.2672 - val_loss: 0.2449\n",
      "Epoch 749/1000\n",
      "0s - loss: 0.2667 - val_loss: 0.2441\n",
      "Epoch 750/1000\n",
      "0s - loss: 0.2667 - val_loss: 0.2439\n",
      "Epoch 751/1000\n",
      "0s - loss: 0.2659 - val_loss: 0.2433\n",
      "Epoch 752/1000\n",
      "0s - loss: 0.2656 - val_loss: 0.2426\n",
      "Epoch 753/1000\n",
      "0s - loss: 0.2654 - val_loss: 0.2423\n",
      "Epoch 754/1000\n",
      "0s - loss: 0.2654 - val_loss: 0.2423\n",
      "Epoch 755/1000\n",
      "0s - loss: 0.2643 - val_loss: 0.2416\n",
      "Epoch 756/1000\n",
      "0s - loss: 0.2640 - val_loss: 0.2410\n",
      "Epoch 757/1000\n",
      "0s - loss: 0.2637 - val_loss: 0.2407\n",
      "Epoch 758/1000\n",
      "0s - loss: 0.2633 - val_loss: 0.2401\n",
      "Epoch 759/1000\n",
      "0s - loss: 0.2630 - val_loss: 0.2399\n",
      "Epoch 760/1000\n",
      "0s - loss: 0.2625 - val_loss: 0.2395\n",
      "Epoch 761/1000\n",
      "0s - loss: 0.2623 - val_loss: 0.2389\n",
      "Epoch 762/1000\n",
      "0s - loss: 0.2619 - val_loss: 0.2389\n",
      "Epoch 763/1000\n",
      "0s - loss: 0.2614 - val_loss: 0.2387\n",
      "Epoch 764/1000\n",
      "0s - loss: 0.2611 - val_loss: 0.2379\n",
      "Epoch 765/1000\n",
      "0s - loss: 0.2609 - val_loss: 0.2376\n",
      "Epoch 766/1000\n",
      "0s - loss: 0.2604 - val_loss: 0.2368\n",
      "Epoch 767/1000\n",
      "0s - loss: 0.2599 - val_loss: 0.2364\n",
      "Epoch 768/1000\n",
      "0s - loss: 0.2597 - val_loss: 0.2360\n",
      "Epoch 769/1000\n",
      "0s - loss: 0.2592 - val_loss: 0.2357\n",
      "Epoch 770/1000\n",
      "0s - loss: 0.2592 - val_loss: 0.2364\n",
      "Epoch 771/1000\n",
      "0s - loss: 0.2585 - val_loss: 0.2362\n",
      "Epoch 772/1000\n",
      "0s - loss: 0.2583 - val_loss: 0.2350\n",
      "Epoch 773/1000\n",
      "0s - loss: 0.2577 - val_loss: 0.2344\n",
      "Epoch 774/1000\n",
      "0s - loss: 0.2572 - val_loss: 0.2339\n",
      "Epoch 775/1000\n",
      "0s - loss: 0.2568 - val_loss: 0.2335\n",
      "Epoch 776/1000\n",
      "0s - loss: 0.2564 - val_loss: 0.2333\n",
      "Epoch 777/1000\n",
      "0s - loss: 0.2561 - val_loss: 0.2332\n",
      "Epoch 778/1000\n",
      "0s - loss: 0.2558 - val_loss: 0.2332\n",
      "Epoch 779/1000\n",
      "0s - loss: 0.2555 - val_loss: 0.2327\n",
      "Epoch 780/1000\n",
      "0s - loss: 0.2551 - val_loss: 0.2321\n",
      "Epoch 781/1000\n",
      "0s - loss: 0.2546 - val_loss: 0.2313\n",
      "Epoch 782/1000\n",
      "0s - loss: 0.2542 - val_loss: 0.2306\n",
      "Epoch 783/1000\n",
      "0s - loss: 0.2539 - val_loss: 0.2302\n",
      "Epoch 784/1000\n",
      "0s - loss: 0.2536 - val_loss: 0.2297\n",
      "Epoch 785/1000\n",
      "0s - loss: 0.2534 - val_loss: 0.2291\n",
      "Epoch 786/1000\n",
      "0s - loss: 0.2530 - val_loss: 0.2288\n",
      "Epoch 787/1000\n",
      "0s - loss: 0.2526 - val_loss: 0.2289\n",
      "Epoch 788/1000\n",
      "0s - loss: 0.2519 - val_loss: 0.2290\n",
      "Epoch 789/1000\n",
      "0s - loss: 0.2524 - val_loss: 0.2285\n",
      "Epoch 790/1000\n",
      "0s - loss: 0.2518 - val_loss: 0.2292\n",
      "Epoch 791/1000\n",
      "0s - loss: 0.2518 - val_loss: 0.2290\n",
      "Epoch 00790: early stopping\n"
     ]
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import io\n",
    "import requests\n",
    "import numpy as np\n",
    "import os\n",
    "from sklearn.model_selection import train_test_split\n",
    "from sklearn import metrics\n",
    "from keras.models import Sequential\n",
    "from keras.layers.core import Dense, Activation\n",
    "from keras.callbacks import EarlyStopping\n",
    "from keras.callbacks import ModelCheckpoint\n",
    "\n",
    "path = \"./data/\"\n",
    "    \n",
    "filename = os.path.join(path,\"iris.csv\")    \n",
    "df = pd.read_csv(filename,na_values=['NA','?'])\n",
    "\n",
    "species = encode_text_index(df,\"species\")\n",
    "x,y = to_xy(df,\"species\")\n",
    "\n",
    "# Split into train/test\n",
    "x_train, x_test, y_train, y_test = train_test_split(    \n",
    "    x, y, test_size=0.25, random_state=42)\n",
    "\n",
    "model = Sequential()\n",
    "model.add(Dense(10, input_dim=x.shape[1], activation='relu'))\n",
    "model.add(Dense(5,activation='relu'))\n",
    "model.add(Dense(y.shape[1],activation='softmax'))\n",
    "model.compile(loss='categorical_crossentropy', optimizer='adam')\n",
    "monitor = EarlyStopping(monitor='val_loss', min_delta=1e-3, patience=5, verbose=1, mode='auto')\n",
    "checkpointer = ModelCheckpoint(filepath=\"best_weights.hdf5\", verbose=0, save_best_only=True) # save best model\n",
    "\n",
    "model.fit(x,y,validation_data=(x_test,y_test),callbacks=[monitor,checkpointer],verbose=2,epochs=1000)\n",
    "model.load_weights('best_weights.hdf5') # load weights from best model"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Potential Keras Issue on Small Networks and Early Stopping\n",
    "\n",
    "You might occasionally see this error:\n",
    "\n",
    "```\n",
    "OSError: Unable to create file (Unable to open file: name = 'best_weights.hdf5', errno = 22, error message = 'invalid argument', flags = 13, o_flags = 302)\n",
    "```\n",
    "\n",
    "Usually you can just run rerun the code and it goes away.  This is an unfortnuate result of saving a file each time the validation score improves (as described in the previous section).  If the errors improve two rapidly, you might try to save the file twice and get an error from these two saves overlapping.  For larger neural networks this will not be a problem because each training step will take longer, allowing for plenty of time for the previous save to complete.  Because of this potential issue, this code is not used with every neural network in this course.  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Calculate Classification Accuracy\n",
    " \n",
    "Accuracy is the number of rows where the neural network correctly predicted the target class.  Accuracy is only used for classification, not regression.\n",
    "\n",
    "$ accuracy = \\frac{\\textit{#} \\ correct}{N} $\n",
    "\n",
    "Where $N$ is the size of the evaluted set (training or validation). Higher accuracy numbers are desired.\n",
    "\n",
    "As we just saw, by default, Keras will return the percent probability for each class. We can change these prediction probabilities into the actual iris predicted with **argmax**."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[1 0 2 1 1 0 1 2 1 1 2 0 0 0 0 1 2 1 1 2 0 2 0 2 2 2 2 2 0 0 0 0 1 0 0 2 1\n",
      " 0]\n"
     ]
    }
   ],
   "source": [
    "pred = np.argmax(pred,axis=1) # raw probabilities to chosen class (highest probability)\n",
    "print(pred)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now that we have the actual iris flower predicted, we can calculate the percent accuracy (how many were correctly classified)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Accuracy score: 1.0\n"
     ]
    }
   ],
   "source": [
    "y_compare = np.argmax(y_test,axis=1) \n",
    "score = metrics.accuracy_score(y_compare, pred)\n",
    "print(\"Accuracy score: {}\".format(score))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Calculate Classification Log Loss\n",
    "\n",
    "Accuracy is like a final exam with no partial credit.  However, neural networks can predict a probability of each of the target classes.  Neural networks will give high probabilities to predictions that are more likely.  Log loss is an error metric that penalizes confidence in wrong answers. Lower log loss values are desired.\n",
    "\n",
    "The following code shows the output of predict_proba:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Numpy array of predictions\n",
      "[  0.0998  67.6963  32.2039]\n",
      "As percent probability\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "array([[ 0.001 ,  0.677 ,  0.322 ],\n",
       "       [ 0.8379,  0.1556,  0.0064],\n",
       "       [ 0.    ,  0.0151,  0.9849],\n",
       "       [ 0.0003,  0.6036,  0.3961],\n",
       "       [ 0.0005,  0.6396,  0.3599]], dtype=float32)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Log loss score: 0.2284613974441431\n"
     ]
    }
   ],
   "source": [
    "from IPython.display import display\n",
    "\n",
    "# Don't display numpy in scientific notation\n",
    "np.set_printoptions(precision=4)\n",
    "np.set_printoptions(suppress=True)\n",
    "\n",
    "# Generate predictions\n",
    "pred = model.predict(x_test)\n",
    "\n",
    "print(\"Numpy array of predictions\")\n",
    "print(pred[0]*100)\n",
    "\n",
    "print(\"As percent probability\")\n",
    "display(pred[0:5])\n",
    "\n",
    "score = metrics.log_loss(y_test, pred)\n",
    "print(\"Log loss score: {}\".format(score))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[Log loss](https://www.kaggle.com/wiki/LogarithmicLoss) is calculated as follows:\n",
    "\n",
    "$ \\text{log loss} = -\\frac{1}{N}\\sum_{i=1}^N {( {y}_i\\log(\\hat{y}_i) \\cdot (1 - {y}_i)\\log(1 - \\hat{y}_i))} $\n",
    "\n",
    "The log function is useful to penalizing wrong answers.  The following code demonstrates the utility of the log function:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\jeffh\\Anaconda3\\lib\\site-packages\\ipykernel_launcher.py:12: RuntimeWarning: divide by zero encountered in log\n",
      "  if sys.path[0] == '':\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtoAAAE0CAYAAAASSJRcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xl8XPV57/Hvo31fLMmSdxtjY4wN\nxMiEJQkyEAJtAs16aRpKuEloSJOUtkmbljRpkiZt6Jb21d423Ca9aZvEoSELNyEESGJICwZsNtsY\n4w0bW7Ily9a+j57+MUdmsGV0RqPRGWk+79dLL82cM8szj49GX//0m98xdxcAAACAqZUTdQEAAADA\nbETQBgAAANKAoA0AAACkAUEbAAAASAOCNgAAAJAGBG0AAAAgDQjaAJDhzOwlM7t6Ch7nGjP7Qcjb\nPmFm56X6nACQzQjaAJA9viTpL0Le9q8kfT6NtQDArEfQBoAsYGbrJVW6++aQd7lX0gYzm5fGsgBg\nViNoA8AMYWaFZvYVM2sOvr5iZoUJ+//AzFqCfR80Mzezs4Pd10l6OOG2l5nZMTNbFFy/wMw6zGyV\nJLn7gKStkq6ZvlcIALMLQRsAZo47JF0i6UJJF0i6WNKnJcnMrpX0e5KulnS2pCtOue9aSbvGrrj7\no5K+KukbZlYs6d8lfdrdX0i4z87geQAAk0DQBoCZ4zckfd7dW929TdLnJN0U7HuPpH919x3u3hfs\nS1QlqfuUbX8qqVLSE5KaJf3jKfu7g/sBACaBoA0AM8d8SQcSrh8Ito3tezlhX+JlSTohqTxxg7sP\nS/p/ktZI+mt391PuUy6pI7WSASB7EbQBYOZolrQk4friYJsktUhamLBv0Sn3fU7SysQNZrZA0mcl\n/aukv06c7x04V9KzKdYMAFmLoA0AM8e3JX3azOrMrFbSZyT9R7Dvbkm3mNm5ZlYS7Et0nxLmbZuZ\nKT6a/TVJH1A8qH8hYX+hpIskPZielwIAsx9BGwBmjj+TtEXx0eltkp4KtsndfyLp7yX9QtIeSY8F\n9xkM9j8lqdPMXh9s/7ikekl/EkwZuUXxoP7GYP/1kja5+9iIOQAgSXb6lDwAwExnZudK2i6p0N1H\ngm3XSPqIu/9aiPs/LukD7r49vZUCwOxF0AaAWcLM3i7px5JKJX1D0miYUA0ASA+mjgDA7PFbktok\n7ZUUk3RbtOUAQHZjRBsAAABIA0a0AQAAgDQgaAMAAABpkBd1Acmora31pUuXRvb8vb29Ki0tjez5\nZzr6N3n0LjX0LzX0b/LoXWroX2roX2q2bt16zN3rUnmMGRW0ly5dqi1btkT2/Js2bVJTU1Nkzz/T\n0b/Jo3epoX+poX+TR+9SQ/9SQ/9SY2YHUn0Mpo4AAAAAaUDQBgAAANKAoA0AAACkAUEbAAAASAOC\nNgAAAJAGBG0AAAAgDQjaAAAAQBoQtAEAAIA0IGgDAAAAaUDQBgAAANKAoA0AAACkAUEbAAAASAOC\nNgAAAJAGBG0AAAAgDQjaAAAAQBoQtAEAAIA0IGgDAAAAaUDQBgAAANKAoA0AAACkAUEbAAAASAOC\nNgAAAJAGBG0AAAAgDQjaAAAAQBoQtAEAAIA0IGgDAAAAaUDQBgAAANKAoA0AAACkAUEbAAAASAOC\nNgAAAJAGBG0AAAAgDQjaAAAAQBoQtAEAAIA0IGgDAAAAaRBp0Daza81sl5ntMbNPRVkLAAAAMJUi\nC9pmlivpHyVdJ2m1pF83s9VR1QMAAABMpShHtC+WtMfd97n7kKSNkm6IsB4AwCmGRkY1POpRlwEA\nM1JehM+9QNLLCdcPSXp9RLUAQFYYHXV19A+rvWdQx3qG1N47qPaeofj13vj39p4htQeXuwZG9L5z\nC/TmqAsHgBnI3KMZqTCzd0t6i7t/MLh+k6SL3f1jp9zuVkm3SlJ9ff1FGzdunPZax/T09KisrCyy\n55/p6N/k0bvUzOb+ubsGYlL3kKtr0NU9HP/eNRT/6g6+x7fFbzfeu75JKiuQKgpMFQWm8uCrosC0\nvHRIa+bNzv6l22w+9qYD/UsN/UvNhg0btrp7YyqPEeWI9iFJixKuL5TUfOqN3P0uSXdJUmNjozc1\nNU1LcePZtGmTonz+mY7+TR69S81M65+7q6t/RG09A2rtHlTb2FfPK6PP8RHnIR3rGdTgyOi4j1NW\nmKeaskLVlBZoQX2hassKVFNaqJqyAtWUFaq2NP69pqxA1SUFys2xcR9npvUvk9C71NC/1NC/6EUZ\ntJ+UtMLMlkk6LOlGSe+NsB4ASKuB4djJwNzWPfjqEB1sPxZcHoqdHp7zc021QTCuKS3U2XPL4tcT\nAvPJy6UFKsrPjeBVAgDGRBa03X3EzD4q6aeSciV93d13RFUPAEzG6KjreN+QWrsG1do98KognRig\n27oH1T0wctr9zaSa0gLVlhWqrrxQy+tKVVdeqLrgel15oeaWF6qurEgVxXkyG3/UGQCQeaIc0Za7\n3yfpvihrAIDxuLu6B0fU2jWgI52DOto1oKPdAzraOaCjXYMnL7d2D2pknFU5ygrzTgbmcxsq9KYV\nhacF6Lry+MhzXi7nDgOA2SjSoA0AURibwnGka0BHuwZ0JAjMp17uG4qddt/yojw1VBSpvqJIlyyv\nOXl5bnmh5lbER55rywtUUsDbKwBkO34TAJhV+odiau7sV0vHgJo7+3Wkc0Bbdw7qX/c9EQ/SXQPq\n6Bs+7X4FeTlBaC7UefMrdOWquaqvKFR9EKQbKoo0t6KQAA0ACI3fGABmjMGRmI50Dqi5Y0Atnf1q\n6Qy+dwyoObg8XoiuKJCW1A1pYXWJGpdWq768SPWVrwTo+opCVRbnM/8ZADClCNoAMoK7q61nUIdO\n9OvQiX41d/SrpaNfzZ3x6Rwtnf061jN02v2qSvI1r7JY8yuLtG5xleZXFWteZVF8W1U8TG/+71+q\nqekNEbwqAEA2I2gDmBajo67W7kEdOtGnwx39QaDu06ET/Tp8ol+HOvo1dMp60OWFeZpXFQ/N582v\n0LzKYs2rKtL84Pu8yiKmcgAAMha/oQBMidio60jXQDw0vypAxy+3dAyctjZ0TWmBFlYXa9W8cl29\nul4Lqoq1sLpYC6tLNK+qSBVF+RG9GgAAUkfQBhBa7+CIDh7vi3+19+nA8V4dPN6vg+29OnSi/7Rl\n7urKC7WwulhrF1TqujXztKA6HqQXVRdrflUxo9EAgFmN33IAThqbJ32wPR6mDwTfxy4f6xl81e0r\nivK0pKZU5y2o1HVr52lRdUkwIh0P0pyZEACQzQjaQJZxd7V1D2rfsV7tT/h6OQjUiWtHm0nzK4u1\naE6xrlo1V4trSrR4TomW1JRoyZxSVZYwtQMAgDMhaAOzVNfAsPa3xUP0K6G6R/vbetWbEKYL8nK0\nJAjPly2v1ZIgTC+uiY9OF+YxKg0AwGQQtIEZbDg2qgPtvdrTmhCkg1CduBRejkkLq0u0rLZUjUvm\n6Ky6Ui2rjX/NqyxWbg7rRwMAMNUI2sAMMBRzPd/cpd2t3drb2qPdwddLx3pf9QHEuvJCLast1dXn\n1p8M0mfVlWrRnBJGpgEAmGYEbSCD9AyOnAzSe1p7tKe1W7tbe3SwvU/+4C8lxUenl9aU6uy5Zbpm\ndb1W1JdpeV2ZltWWqpzl8AAAyBgEbSACw7FR7Wvr1QtHuvTCkW690NKlF4/26HBH/8nb5Oeazqot\n05oFlbqwalhvfv0arZhbrqW1jE4DADATELSBNHKPnw1xZ0s8UO860q2dLV3a29aj4Vh8ykd+rml5\nXZnWL63We+sX6+y5ZVoxt0yL55QoLzdHkrRp0yY1nT8/ypcCAACSRNAGpsjQyKh2t3Zrx+EuPd/S\ndXK0uqNv+ORt5lUW6ZyGcjWdM1fnzivXOQ3lOqu2TAV5ORFWDgAA0oGgDUzCwHBMu450a3tzp7Yf\n7tT2w13adaT75CnGSwpytbK+XNetadCqhgqd01CuVQ3lqiopiLhyAAAwXQjawAR6B0e0s6UrHqib\n4993t/YoFqz2UVmcrzULKnTL5Ut13oJKrZlfoaU1pcphyTwAALIaQRtIEBt17W7t1jMHO/TMy/Gv\nF492a2wFvZrSAq1ZUKmrzp2rtQsqdd78Si2sLpYZoRoAALwaQRtZ7WjXgJ4+GapPaNuhzpNnTaws\nztcFi6p0zep6rV1YpbULKlVfUUioBgAAoRC0kTVGYqN6vqVLT750QlsPHNfTBzvU0jkgKb7yx+p5\nFXrnRQt14aIqXbioSstqSwnVAABg0gjamLV6B0f0zMsdevKl49ry0gk9dfCE+oLR6oXVxVq/dE48\nVC+u0up5FSrKZ21qAAAwdQjamDVO9A7p8f3tevKlE9ry0nFtb+5SbNRlJp3bUKF3X7RQ65fNUeOS\nOWqoLIq6XAAAMMsRtDFjdQ8M64n9x/Xo3nY9urddO1u6JEmFeTl63eIqfaRpuRqXztHrFlepglOT\nAwCAaUbQxozRPxTT1gMn9OjeY3p0b7u2He5UbNRVmJejxqXV+sQ1K3Xp8hqtXVDFCWAAAEDkCNrI\nWO6uXUe79fCuNj38Ypu2vHRCQ7FR5eWYLlgUH7G+dHmN1i2uZn41AADIOARtZJTO/mH9955jJ8P1\nka74qiCrGsp182VLdNnZtVq/dI7KCjl0AQBAZiOtIFLurheP9uihnUe1aVernjrYodioq7wwT29Y\nUaumc+r0ppV1mldZHHWpAAAASSFoY9rFRl1bD5zQAzuO6MGdR3WgvU+SdN78Cn34irN0xcq5et3i\nKuXnMs8aAADMXARtTIv+oZiebh3Rj//zWf3shVYd7x1SQW6OLl1eow+98Sy9eXW96itYcg8AAMwe\nBG2kzcBwTJt2telHzzXrZztb1T8cU3nREV25aq7evLpeV6ysUznL7gEAgFmKoI0pNTQyqv/a06Yf\nPduiB54/qp7BEc0pLdDb1y3QglirPvRrG1h6DwAAZAWCNlLm7npi/3F9/+nD+sn2I+rsH1ZFUZ5+\nZW2D3nr+fF22vEZ5uTnatGkTIRsAAGQNgjYm7eXjffreU4d1z1OHdPB4n0oKcvWW8xr01vPn6Y0r\n6gjVAAAgqxG0kZS+oRH9ZNsRfXfrIT22r12SdNnyGt1+9Qpdu6ZBJQUcUgAAABJBGyHtae3Wf2w+\nqHu2HlL34IgWzynR7715pd7+ugVaNKck6vIAAAAyDkEbZzQ0MqoHnj+i/9h8QJv3HVdBbo6uW9ug\n9168WBcvmyMzi7pEAACAjEXQxmnaewb1b48d0LeeOKi27kEtrC7WH167Su9uXKjassKoywMAAJgR\nCNo4aW9bj772X/t1z9ZDGhwZ1ZWr5uqmS5boTSvrlJvD6DUAAEAyCNrQky8d112P7NNDO48qPzdH\n71y3UB94wzKdPbcs6tIAAABmLIJ2Fntsb7u+8tCLenz/cVWX5OtjV67Qb166hOkhAAAAU4CgnWXc\nXY/ta9dXHtqtJ/Yf19zyQn32bat14/rFKi7Ijbo8AACAWYOgnUW2HjiuL9+/62TA/tO3rdaNFy9W\nUT4BGwAAYKpFErTN7C8lvU3SkKS9km5x944oaskG+4/16s77X9BPth9RXXmhPnf9efpf6xcRsAEA\nANIoqhHtByX9kbuPmNmXJf2RpD+MqJZZq71nUH//s9365uMHVZCXo9+9eqU+9KZlnL0RAABgGkSS\nuNz9gYSrmyW9K4o6ZqvYqOubjx/QX/50l/qGYrpx/SLdfvVK1ZXzIUcAAIDpkglDm/9b0neiLmK2\nePrgCf3JD7dr++EuXX52jT53/Xk6e2551GUBAABkHXP39Dyw2UOSGsbZdYe7/zC4zR2SGiW9w89Q\niJndKulWSaqvr79o48aNaak3jJ6eHpWVZeba0n3DrrtfHNLDL4+ostD066sKdHFDbkadJj2T+5fp\n6F1q6F9q6N/k0bvU0L/U0L/UbNiwYau7N6byGGkL2hM+sdnNkj4s6Sp37wtzn8bGRt+yZUt6C3sN\nmzZtUlNTU2TPfyYPv9imT93znI52DeiWy5fp9qtXqLwoP+qyTpOp/ZsJ6F1q6F9q6N/k0bvU0L/U\n0L/UmFnKQTuqVUeuVfzDj1eEDdk4XffAsL74453a+OTLOntumb73kct14aKqqMsCAACAopuj/Q+S\nCiU9GExt2OzuH46olhlp64ET+vi3n1ZLZ78+fMVy3X71CpbrAwAAyCBRrTpydhTPOxuMjrq++sg+\n/dUDuzS/qkjfve0yrVtcHXVZAAAAOEUmrDqCkNp7BvW7dz+rR15s06+eP09//o61qsjAudgAAAAg\naM8YO5o79aFvbNGx3iF98e1r9N6LF2fUiiIAAAB4NYL2DHDfthb9/t3PqqokX/d8+DKtXVgZdUkA\nAACYAEE7g7m7/u5nu/WVh3Zr3eIq/fNNF2lueVHUZQEAACAEgnaGio26Pv2D7fr2Ewf1znUL9aV3\nrFFhHquKAAAAzBQE7Qw0MBzT7Ruf0f07jui3NyzXJ645h/nYAAAAMwxBO8P0D8X0gW88qUf3tutP\n3rpaH3jDsqhLAgAAwCQQtDPIwHA8ZG/e166/ec8Fese6hVGXBAAAgEkiaGeIgeGYPvRvW/TYvnb9\n9bsJ2QAAADNdTtQFQBqJjeqj33pKv9x9TF9+5/mEbAAAgFmAoB0xd9dn7t2hh3a26gs3nKf3NC6K\nuiQAAABMAYJ2xP7p4b361uMHdVvTct106dKoywEAAMAUIWhH6L5tLbrz/l264cL5+uQ150RdDgAA\nAKYQQTsie1q79Yn/fFbrFlfpznedr5wc1skGAACYTQjaEegeGNat/75VJQV5+qf3XcQZHwEAAGYh\nlvebZu6uT31vmw609+mbH3y96iuKoi4JAAAAacCI9jT7/tOH9ePnWvT716zUJWfVRF0OAAAA0oSg\nPY0Od/Trsz/cofVLq/Vbb1oedTkAAABII4L2NHF3feLuZzXqrr95z4XK5cOPAAAAsxpBe5rc89Rh\nPbavXXf86motmlMSdTkAAABIM4L2NOjoG9KX7tupdYurdON6zvwIAACQDQja0+DL9+9SZ/+wvvj2\ntayXDQAAkCUI2mm2s6VLG588qJsvXapz51VEXQ4AAACmCUE7ze68/wWVF+bpd65aEXUpAAAAmEYE\n7TTavK9dv9jVpo9sOFuVJflRlwMAAIBpRNBOE3fXnfe/oIaKIr3/sqVRlwMAAIBpNmHQNrOPmln1\ndBQzm2zed1xPHezQb29YrqL83KjLAQAAwDQLM6LdIOlJM7vbzK41M5bNCOGfHt6r2rICvbuR5fwA\nAACy0YRB290/LWmFpK9Jer+k3Wb2JTPjHOJnsKO5U4+82KZbLl/GaDYAAECWCjVH291d0pHga0RS\ntaTvmtmdaaxtxvqXX+5XWWGe3nfJkqhLAQAAQETyJrqBmX1c0s2Sjkn6F0mfdPdhM8uRtFvSH6S3\nxJnlRO+QfrytRTeuX6TKYlYaAQAAyFYTBm1JtZLe4e4HEje6+6iZvTU9Zc1c9zx1SEMjo3rv6xdH\nXQoAAAAiNGHQdvfPvMa+nVNbzszm7vrWEwe1bnGVVjVwFkgAAIBsxjraU2jLgRPa19ar976eudkA\nAADZjqA9he59pllF+Tm6bk1D1KUAAAAgYgTtKTISG9V921p01bn1Ki0MM/UdAAAAsxlBe4o8tq9d\n7b1Detv586MuBQAAABmAoD1FfvRsi8oK89R0Tl3UpQAAACADELSnwOio62cvtKrpnDrOBAkAAABJ\nBO0psaO5S8d6BnXlqrlRlwIAAIAMQdCeAr/Y1Soz6YqVTBsBAABAHEF7CvxiV6suWFilmrLCqEsB\nAABAhiBop6izb1jPvNzBaDYAAABeJdKgbWafMDM3s9oo60jFlgPH5S5durwm6lIAAACQQSIL2ma2\nSNKbJR2Mqoap8MT+4yrIzdGFi6qiLgUAAAAZJMoR7b+V9AeSPMIaUvbES8d1/sJKlvUDAADAq0QS\ntM3sekmH3f3ZKJ5/qvQNjWjboU5dvGxO1KUAAAAgw5h7egaUzewhSQ3j7LpD0h9LusbdO83sJUmN\n7n7sDI9zq6RbJam+vv6ijRs3pqXeMHp6elRWVnby+s72mL785IB+96JCXVCXF1ldM8Wp/UN49C41\n9C819G/y6F1q6F9q6F9qNmzYsNXdG1N5jLSlQ3e/erztZrZW0jJJz5qZJC2U9JSZXezuR8Z5nLsk\n3SVJjY2N3tTUlK6SJ7Rp0yYlPv/uR/ZJ2qn3XfdGlvYL4dT+ITx6lxr6lxr6N3n0LjX0LzX0L3rT\nPgzr7tsknTyF4kQj2pls2+FOza8sImQDAADgNKyjnYLtzZ1as6Ay6jIAAACQgSIP2u6+dCaOZvcM\njmj/sV6CNgAAAMYVedCeqZ5v7pK7tJagDQAAgHEQtCdp19FuSdKqeeURVwIAAIBMRNCepL2tPSot\nyFVDRVHUpQAAACADEbQnaW9bj5bPLVOwRCEAAADwKgTtSdrb2qPldSwCDwAAgPERtCehd3BEzZ0D\nWl5XGnUpAAAAyFAE7UnYf6xXkhjRBgAAwBkRtCdhLGgvrWVEGwAAAOMjaE/CoRP9kqSF1cURVwIA\nAIBMRdCehMMdfaoszld5UX7UpQAAACBDEbQn4fCJfi2oYjQbAAAAZ0bQnoTDHf1awLQRAAAAvAaC\ndpLcnRFtAAAATIignaTO/mH1DsX4ICQAAABeE0E7Sc0dA5Kk+YxoAwAA4DUQtJN0rGdQklRXXhhx\nJQAAAMhkBO0ktXXHg3ZtGUEbAAAAZ0bQTtLYiHZtWUHElQAAACCTEbSTdKxnUEX5OSorzIu6FAAA\nAGQwgnaSjvUMqbasUGYWdSkAAADIYATtJB3rGWR+NgAAACZE0E5SWzdBGwAAABMjaCcpPnWED0IC\nAADgtRG0k+Du6uofVmVJftSlAAAAIMMRtJMwPCoNxUZVWUzQBgAAwGsjaCehd9gliaANAACACRG0\nk9A3Ev9eUUTQBgAAwGsjaCehjxFtAAAAhETQTsLY1JEKgjYAAAAmQNBOwtjUEUa0AQAAMBGCdhLG\npo5UFOVFXAkAAAAyHUE7CUwdAQAAQFgE7ST0j7iK83OVn0vbAAAA8NpIjEkYjEklBblRlwEAAIAZ\ngKCdhMGYVEzQBgAAQAgE7SQMxpwRbQAAAIRC0E5CfESbFUcAAAAwMYJ2EoZirpJ8RrQBAAAwMYJ2\nEpijDQAAgLAI2kkYjDlBGwAAAKEQtJMwFBNTRwAAABAKQTsJrDoCAACAsAjaSRgcYdURAAAAhBNZ\n0Dazj5nZLjPbYWZ3RlVHWCOxUY24VMzUEQAAAIQQyfCsmW2QdIOk89190MzmRlFHMvqHY5Kk4gL+\nCAAAAICJRZUab5P0F+4+KEnu3hpRHaENjYxKkgrzGNEGAADAxKIK2islvdHMHjezh81sfUR1hDYc\nc0lSfi4j2gAAAJiYuXt6HtjsIUkN4+y6Q9IXJf1c0u9IWi/pO5LO8nGKMbNbJd0qSfX19Rdt3Lgx\nLfVOpK1vVJ98pF8fXFugNyzIj6SGma6np0dlZWVRlzEj0bvU0L/U0L/Jo3epoX+poX+p2bBhw1Z3\nb0zlMdI2R9vdrz7TPjO7TdL3gmD9hJmNSqqV1DbO49wl6S5Jamxs9KampvQUPIE9rT3SIw9r7Xmr\n1XThgkhqmOk2bdqkqP79Zjp6lxr6lxr6N3n0LjX0LzX0L3pRzYP4gaQrJcnMVkoqkHQsolpCGY7F\n52gXMHUEAAAAIUS1KPTXJX3dzLZLGpJ083jTRjLJyaCdR9AGAADAxCIJ2u4+JOl9UTz3ZI2tOsKH\nIQEAABAGqTGkoRhBGwAAAOGRGkMaW96PqSMAAAAIg9QY0vAIH4YEAABAeKTGkE5OHcmziCsBAADA\nTEDQDonl/QAAAJAMUmNIrDoCAACAZJAaQxpiHW0AAAAkgdQYEh+GBAAAQDJIjSGNLe+Xl8uHIQEA\nADAxgnZII6NB0M6hZQAAAJgYqTGkUY8HbXI2AAAAwiA2hhQLRrRzjakjAAAAmBhBO6STQTuHoA0A\nAICJEbRDGnWXSTJGtAEAABACQTuk2KiLwWwAAACERdAOKeYuBrMBAAAQFkE7pFFGtAEAAJAEgnZI\nsVGaBQAAgPDIjiGNOiPaAAAACI+gHRIfhgQAAEAyCNohxRjRBgAAQBII2iHFPwxJ0gYAAEA4BO2Q\nmDoCAACAZBC0Q4oFZ4YEAAAAwiBoh8Q62gAAAEgGQTukmIugDQAAgNAI2iExog0AAIBkELRD4sOQ\nAAAASAZBO6T4OtokbQAAAIRD0A6JqSMAAABIBkE7pJg7zQIAAEBoZMeQYqMuZo4AAAAgLIJ2SM7y\nfgAAAEgCQTukUc4MCQAAgCQQtENyF1NHAAAAEBpBOySXR10CAAAAZhCCdhIY0AYAAEBYBO2QnAFt\nAAAAJIGgHRI5GwAAAMkgaCeBD0MCAAAgLIJ2WAxpAwAAIAkE7ZBcrKMNAACA8AjaIfFhSAAAACQj\nkqBtZhea2WYze8bMtpjZxVHUkSzmaAMAACCsqEa075T0OXe/UNJngusZjQFtAAAAJCOqoO2SKoLL\nlZKaI6ojNHfmaAMAACC8vIie93ZJPzWzv1I87F8WUR0AAABAWpin6VN+ZvaQpIZxdt0h6SpJD7v7\nPWb2Hkm3uvvVZ3icWyXdGlw9R9KudNQbUq2kYxE+/0xH/yaP3qWG/qWG/k0evUsN/UsN/UvNOe5e\nnsoDpC1ov+aTmnVKqnJ3NzOT1OnuFRPdL2pmtsXdG6OuY6aif5NH71JD/1JD/yaP3qWG/qWG/qVm\nKvoX1RztZklXBJevlLQ7ojoAAACAtIhqjvaHJP2dmeVJGtArU0MAAACAWSGSoO3u/yXpoiieO0V3\nRV3ADEf/Jo/epYb+pYb+TR69Sw39Sw39S03K/YtkjjYAAAAw23EKdgAAACANCNqSzOxaM9tlZnvM\n7FPj7C80s+8E+x83s6UJ+/4o2L7LzN4ynXVnihD9+z0ze97MnjOzn5nZkoR9MTN7Jvi6d3orzwwh\n+vd+M2tL6NMHE/bdbGa7g6+bp7fy6IXo3d8m9O1FM+tI2MexZ/Z1M2s1s+1n2G9m9vdBf58zs3UJ\n+7L92Juod78R9Ow5M3vUzC7AMmE3AAAGUElEQVRI2PeSmW0Ljr0t01d15gjRvyYz60z4Gf1Mwr7X\n/LnPBiH698mE3m0P3u/mBPuy+vgzs0Vm9gsz22lmO8zsd8a5zdS997l7Vn9JypW0V9JZkgokPStp\n9Sm3+Yikfw4u3yjpO8Hl1cHtCyUtCx4nN+rXlIH92yCpJLh821j/gus9Ub+GGdC/90v6h3HuO0fS\nvuB7dXC5OurXlEm9O+X2H5P09YTrWX3sBT14k6R1krafYf+vSPqJJJN0iaTHg+1ZfeyF7N1lYz2R\ndN1Y74LrL0mqjfo1ZHj/miT9aJztSf3cz9avifp3ym3fJunnCdez+viTNE/SuuByuaQXx/m9O2Xv\nfYxoSxdL2uPu+9x9SNJGSTeccpsbJH0juPxdSVeZmQXbN7r7oLvvl7QneLxsMmH/3P0X7t4XXN0s\naeE015jJwhx/Z/IWSQ+6+3F3PyHpQUnXpqnOTJRs735d0renpbIZwt0fkXT8NW5yg6R/87jNkqrM\nbJ449ibsnbs/GvRG4n3vNCGOvTNJ5T1z1kiyf7z3JXD3Fnd/KrjcLWmnpAWn3GzK3vsI2vHmvpxw\n/ZBOb/jJ27j7iKROSTUh7zvbJduDDyj+v8QxRWa2xcw2m9mvpaPADBe2f+8M/nz1XTNblOR9Z6vQ\nrz+YrrRM0s8TNmf7sRfGmXqc7cdesk5933NJD5jZVouf/Rjju9TMnjWzn5jZecE2jr0kmFmJ4kHw\nnoTNHH8Bi08Ffp2kx0/ZNWXvfVGto51JbJxtpy7FcqbbhLnvbBe6B2b2PkmNeuVkRZK02N2bzews\nST83s23uvjcNdWaqMP37/5K+7e6DZvZhxf+6cmXI+85mybz+GyV9191jCduy/dgLg/e+FJnZBsWD\n9hsSNl8eHHtzJT1oZi8EI5R4xVOSlrh7j5n9iqQfSFohjr1kvU3Sf7t74ug3x58kMytT/D8gt7t7\n16m7x7nLpN77GNGO/29kUcL1hYqfuXLc21j8JDuViv/JJsx9Z7tQPTCzqyXdIel6dx8c2+7uzcH3\nfZI2Kf4/y2wyYf/cvT2hZ/9Xr6xBn+3HXzKv/0ad8qdTjr1QztTjbD/2QjGz8yX9i6Qb3L19bHvC\nsdcq6fvKvimHE3L3LnfvCS7fJynfzGrFsZes13rvy9rjz8zyFQ/Z33T3741zkyl77yNoS09KWmFm\ny8ysQPGD8tQVCO6VNPbJ0ncp/qECD7bfaPFVSZYp/r/tJ6ap7kwxYf/M7HWSvqp4yG5N2F5tZoXB\n5VpJl0t6ftoqzwxh+jcv4er1is8nk6SfSrom6GO1pGuCbdkizM+uzOwcxT+08ljCNo69cO6V9JvB\nJ/AvkdTp7i3i2JuQmS2W9D1JN7n7iwnbS82sfOyy4r0bd+WIbGZmDcFnoWRmFyueV9oV8ucekplV\nKv4X5B8mbMv64y84rr4maae7/80ZbjZl731ZP3XE3UfM7KOKNypX8VUJdpjZ5yVtcfd7Ff8H+Xcz\n26P4SPaNwX13mNndiv+CHpH026f8aXrWC9m/v5RUJuk/g/fNg+5+vaRzJX3VzEYVfxP9C3fPqrAT\nsn8fN7PrFT/Gjiu+Conc/biZfUHxXzyS9PlT/jw4q4XsnRT/INDG4D/HY7L+2JMkM/u24qs71JrZ\nIUmflZQvSe7+z5LuU/zT93sk9Um6JdiX1ceeFKp3n1H8szz/J3jfG3H3Rkn1kr4fbMuT9C13v3/a\nX0DEQvTvXZJuM7MRSf2Sbgx+hsf9uY/gJUQqRP8k6e2SHnD33oS7cvzFB1ZukrTNzJ4Jtv2xpMXS\n1L/3cWZIAAAAIA2YOgIAAACkAUEbAAAASAOCNgAAAJAGBG0AAAAgDQjaAAAAQBoQtAEAAIA0IGgD\nAAAAaUDQBoBZyMzWm9lzZlYUnA1uh5mtibouAMgmnLAGAGYpM/szSUWSiiUdcvc/j7gkAMgqBG0A\nmKXMrEDxUwUPSLrM3WMRlwQAWYWpIwAwe82RVCapXPGRbQDANGJEGwBmKTO7V9JGScskzXP3j0Zc\nEgBklbyoCwAATD0z+01JI+7+LTPLlfSomV3p7j+PujYAyBaMaAMAAABpwBxtAAAAIA0I2gAAAEAa\nELQBAACANCBoAwAAAGlA0AYAAADSgKANAAAApAFBGwAAAEgDgjYAAACQBv8DbwF8veD+hP8AAAAA\nSUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1797e29c208>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "from matplotlib.pyplot import figure, show\n",
    "from numpy import arange, sin, pi\n",
    "\n",
    "#t = arange(1e-5, 5.0, 0.00001)\n",
    "#t = arange(1.0, 5.0, 0.00001) # computer scientists\n",
    "t = arange(0.0, 1.0, 0.00001)  # data     scientists\n",
    "\n",
    "fig = figure(1,figsize=(12, 10))\n",
    "\n",
    "ax1 = fig.add_subplot(211)\n",
    "ax1.plot(t, np.log(t))\n",
    "ax1.grid(True)\n",
    "ax1.set_ylim((-8, 1.5))\n",
    "ax1.set_xlim((-0.1, 2))\n",
    "ax1.set_xlabel('x')\n",
    "ax1.set_ylabel('y')\n",
    "ax1.set_title('log(x)')\n",
    "\n",
    "show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Evaluating Regression Results\n",
    "\n",
    "Regression results are evaluated differently than classification.  Consider the following code that trains a neural network for the [MPG dataset](https://github.com/jeffheaton/t81_558_deep_learning/blob/master/datasets_mpg.ipynb). "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Train on 398 samples, validate on 100 samples\n",
      "Epoch 1/1000\n",
      "0s - loss: 150811.4942 - val_loss: 884.9704\n",
      "Epoch 2/1000\n",
      "0s - loss: 12992.9749 - val_loss: 16838.9134\n",
      "Epoch 3/1000\n",
      "0s - loss: 6151.2360 - val_loss: 858.9771\n",
      "Epoch 4/1000\n",
      "0s - loss: 1857.3857 - val_loss: 961.7318\n",
      "Epoch 5/1000\n",
      "0s - loss: 461.6950 - val_loss: 453.1432\n",
      "Epoch 6/1000\n",
      "0s - loss: 387.2971 - val_loss: 215.2402\n",
      "Epoch 7/1000\n",
      "0s - loss: 269.6235 - val_loss: 228.8721\n",
      "Epoch 8/1000\n",
      "0s - loss: 256.1139 - val_loss: 215.1607\n",
      "Epoch 9/1000\n",
      "0s - loss: 254.5415 - val_loss: 210.5670\n",
      "Epoch 10/1000\n",
      "0s - loss: 246.8290 - val_loss: 207.9672\n",
      "Epoch 11/1000\n",
      "0s - loss: 245.7514 - val_loss: 206.1007\n",
      "Epoch 12/1000\n",
      "0s - loss: 245.1531 - val_loss: 204.1335\n",
      "Epoch 13/1000\n",
      "0s - loss: 244.5450 - val_loss: 202.1923\n",
      "Epoch 14/1000\n",
      "0s - loss: 242.1365 - val_loss: 200.1411\n",
      "Epoch 15/1000\n",
      "0s - loss: 236.9681 - val_loss: 197.8704\n",
      "Epoch 16/1000\n",
      "0s - loss: 233.7680 - val_loss: 196.0665\n",
      "Epoch 17/1000\n",
      "0s - loss: 232.9960 - val_loss: 195.0739\n",
      "Epoch 18/1000\n",
      "0s - loss: 230.3290 - val_loss: 191.7835\n",
      "Epoch 19/1000\n",
      "0s - loss: 227.9730 - val_loss: 190.2279\n",
      "Epoch 20/1000\n",
      "0s - loss: 225.6483 - val_loss: 186.8743\n",
      "Epoch 21/1000\n",
      "0s - loss: 222.4072 - val_loss: 185.6266\n",
      "Epoch 22/1000\n",
      "0s - loss: 219.8205 - val_loss: 182.2413\n",
      "Epoch 23/1000\n",
      "0s - loss: 217.5886 - val_loss: 180.0197\n",
      "Epoch 24/1000\n",
      "0s - loss: 214.6021 - val_loss: 177.9827\n",
      "Epoch 25/1000\n",
      "0s - loss: 211.0622 - val_loss: 175.6313\n",
      "Epoch 26/1000\n",
      "0s - loss: 209.4602 - val_loss: 173.5214\n",
      "Epoch 27/1000\n",
      "0s - loss: 208.1123 - val_loss: 171.3492\n",
      "Epoch 28/1000\n",
      "0s - loss: 206.0525 - val_loss: 168.4899\n",
      "Epoch 29/1000\n",
      "0s - loss: 201.1731 - val_loss: 165.9868\n",
      "Epoch 30/1000\n",
      "0s - loss: 199.2211 - val_loss: 163.5638\n",
      "Epoch 31/1000\n",
      "0s - loss: 198.5224 - val_loss: 161.2110\n",
      "Epoch 32/1000\n",
      "0s - loss: 193.5752 - val_loss: 159.2851\n",
      "Epoch 33/1000\n",
      "0s - loss: 192.2848 - val_loss: 156.6271\n",
      "Epoch 34/1000\n",
      "0s - loss: 189.1332 - val_loss: 154.3807\n",
      "Epoch 35/1000\n",
      "0s - loss: 186.0450 - val_loss: 152.1636\n",
      "Epoch 36/1000\n",
      "0s - loss: 183.5400 - val_loss: 149.8624\n",
      "Epoch 37/1000\n",
      "0s - loss: 181.0363 - val_loss: 149.3006\n",
      "Epoch 38/1000\n",
      "0s - loss: 182.1456 - val_loss: 146.0300\n",
      "Epoch 39/1000\n",
      "0s - loss: 178.2740 - val_loss: 145.2686\n",
      "Epoch 40/1000\n",
      "0s - loss: 173.4227 - val_loss: 141.6625\n",
      "Epoch 41/1000\n",
      "0s - loss: 172.7585 - val_loss: 138.9251\n",
      "Epoch 42/1000\n",
      "0s - loss: 171.8932 - val_loss: 139.5851\n",
      "Epoch 43/1000\n",
      "0s - loss: 167.3690 - val_loss: 134.6213\n",
      "Epoch 44/1000\n",
      "0s - loss: 164.6388 - val_loss: 132.6882\n",
      "Epoch 45/1000\n",
      "0s - loss: 163.2543 - val_loss: 130.6598\n",
      "Epoch 46/1000\n",
      "0s - loss: 159.7163 - val_loss: 130.7230\n",
      "Epoch 47/1000\n",
      "0s - loss: 160.2630 - val_loss: 130.3859\n",
      "Epoch 48/1000\n",
      "0s - loss: 155.7280 - val_loss: 129.1296\n",
      "Epoch 49/1000\n",
      "0s - loss: 154.7054 - val_loss: 124.2799\n",
      "Epoch 50/1000\n",
      "0s - loss: 153.4917 - val_loss: 120.9803\n",
      "Epoch 51/1000\n",
      "0s - loss: 148.1566 - val_loss: 120.0594\n",
      "Epoch 52/1000\n",
      "0s - loss: 147.0108 - val_loss: 121.3335\n",
      "Epoch 53/1000\n",
      "0s - loss: 145.7593 - val_loss: 120.7526\n",
      "Epoch 54/1000\n",
      "0s - loss: 144.0165 - val_loss: 114.3476\n",
      "Epoch 55/1000\n",
      "0s - loss: 142.2844 - val_loss: 111.7173\n",
      "Epoch 56/1000\n",
      "0s - loss: 139.9215 - val_loss: 110.6866\n",
      "Epoch 57/1000\n",
      "0s - loss: 137.0591 - val_loss: 110.1224\n",
      "Epoch 58/1000\n",
      "0s - loss: 134.2647 - val_loss: 106.2880\n",
      "Epoch 59/1000\n",
      "0s - loss: 132.1356 - val_loss: 104.8116\n",
      "Epoch 60/1000\n",
      "0s - loss: 130.9806 - val_loss: 110.8378\n",
      "Epoch 61/1000\n",
      "0s - loss: 136.4747 - val_loss: 102.4049\n",
      "Epoch 62/1000\n",
      "0s - loss: 130.5400 - val_loss: 108.2946\n",
      "Epoch 63/1000\n",
      "0s - loss: 131.0337 - val_loss: 98.4562\n",
      "Epoch 64/1000\n",
      "0s - loss: 123.1012 - val_loss: 97.1794\n",
      "Epoch 65/1000\n",
      "0s - loss: 120.5692 - val_loss: 94.7868\n",
      "Epoch 66/1000\n",
      "0s - loss: 119.1059 - val_loss: 93.4718\n",
      "Epoch 67/1000\n",
      "0s - loss: 117.2886 - val_loss: 98.6343\n",
      "Epoch 68/1000\n",
      "0s - loss: 121.2075 - val_loss: 94.3916\n",
      "Epoch 69/1000\n",
      "0s - loss: 116.9649 - val_loss: 90.7541\n",
      "Epoch 70/1000\n",
      "0s - loss: 114.2378 - val_loss: 88.7964\n",
      "Epoch 71/1000\n",
      "0s - loss: 110.7137 - val_loss: 86.2426\n",
      "Epoch 72/1000\n",
      "0s - loss: 110.4869 - val_loss: 92.2123\n",
      "Epoch 73/1000\n",
      "0s - loss: 115.3338 - val_loss: 92.0765\n",
      "Epoch 74/1000\n",
      "0s - loss: 108.4319 - val_loss: 82.5144\n",
      "Epoch 75/1000\n",
      "0s - loss: 104.0291 - val_loss: 83.0405\n",
      "Epoch 76/1000\n",
      "0s - loss: 103.7248 - val_loss: 81.9978\n",
      "Epoch 77/1000\n",
      "0s - loss: 103.3422 - val_loss: 81.6308\n",
      "Epoch 78/1000\n",
      "0s - loss: 100.5476 - val_loss: 78.0558\n",
      "Epoch 79/1000\n",
      "0s - loss: 99.8780 - val_loss: 75.9163\n",
      "Epoch 80/1000\n",
      "0s - loss: 99.3813 - val_loss: 82.9997\n",
      "Epoch 81/1000\n",
      "0s - loss: 99.0231 - val_loss: 75.8720\n",
      "Epoch 82/1000\n",
      "0s - loss: 93.7391 - val_loss: 72.6853\n",
      "Epoch 83/1000\n",
      "0s - loss: 93.0160 - val_loss: 71.5996\n",
      "Epoch 84/1000\n",
      "0s - loss: 92.3786 - val_loss: 70.5565\n",
      "Epoch 85/1000\n",
      "0s - loss: 90.3970 - val_loss: 70.7920\n",
      "Epoch 86/1000\n",
      "0s - loss: 89.6167 - val_loss: 70.5023\n",
      "Epoch 87/1000\n",
      "0s - loss: 87.9159 - val_loss: 67.3543\n",
      "Epoch 88/1000\n",
      "0s - loss: 86.9525 - val_loss: 66.2694\n",
      "Epoch 89/1000\n",
      "0s - loss: 86.2118 - val_loss: 65.6467\n",
      "Epoch 90/1000\n",
      "0s - loss: 85.2408 - val_loss: 64.5314\n",
      "Epoch 91/1000\n",
      "0s - loss: 83.5046 - val_loss: 63.5319\n",
      "Epoch 92/1000\n",
      "0s - loss: 82.6090 - val_loss: 64.9609\n",
      "Epoch 93/1000\n",
      "0s - loss: 81.1910 - val_loss: 62.1564\n",
      "Epoch 94/1000\n",
      "0s - loss: 79.9782 - val_loss: 62.2616\n",
      "Epoch 95/1000\n",
      "0s - loss: 79.5162 - val_loss: 63.0424\n",
      "Epoch 96/1000\n",
      "0s - loss: 82.2127 - val_loss: 59.5661\n",
      "Epoch 97/1000\n",
      "0s - loss: 77.7171 - val_loss: 60.6299\n",
      "Epoch 98/1000\n",
      "0s - loss: 80.9684 - val_loss: 64.2177\n",
      "Epoch 99/1000\n",
      "0s - loss: 82.7562 - val_loss: 59.4830\n",
      "Epoch 100/1000\n",
      "0s - loss: 76.0679 - val_loss: 56.5858\n",
      "Epoch 101/1000\n",
      "0s - loss: 73.6569 - val_loss: 59.7040\n",
      "Epoch 102/1000\n",
      "0s - loss: 73.2199 - val_loss: 55.5563\n",
      "Epoch 103/1000\n",
      "0s - loss: 71.3402 - val_loss: 54.2142\n",
      "Epoch 104/1000\n",
      "0s - loss: 73.6433 - val_loss: 53.4801\n",
      "Epoch 105/1000\n",
      "0s - loss: 70.4874 - val_loss: 54.4136\n",
      "Epoch 106/1000\n",
      "0s - loss: 71.7230 - val_loss: 60.5264\n",
      "Epoch 107/1000\n",
      "0s - loss: 71.9563 - val_loss: 51.8378\n",
      "Epoch 108/1000\n",
      "0s - loss: 67.4697 - val_loss: 51.0925\n",
      "Epoch 109/1000\n",
      "0s - loss: 67.9659 - val_loss: 57.1168\n",
      "Epoch 110/1000\n",
      "0s - loss: 68.3725 - val_loss: 50.1634\n",
      "Epoch 111/1000\n",
      "0s - loss: 67.8489 - val_loss: 49.5179\n",
      "Epoch 112/1000\n",
      "0s - loss: 68.8507 - val_loss: 48.8458\n",
      "Epoch 113/1000\n",
      "0s - loss: 64.2464 - val_loss: 48.3410\n",
      "Epoch 114/1000\n",
      "0s - loss: 63.8923 - val_loss: 48.2327\n",
      "Epoch 115/1000\n",
      "0s - loss: 62.5394 - val_loss: 47.3220\n",
      "Epoch 116/1000\n",
      "0s - loss: 61.8350 - val_loss: 49.5652\n",
      "Epoch 117/1000\n",
      "0s - loss: 61.7264 - val_loss: 50.3741\n",
      "Epoch 118/1000\n",
      "0s - loss: 62.0407 - val_loss: 47.4053\n",
      "Epoch 119/1000\n",
      "0s - loss: 60.3533 - val_loss: 47.2662\n",
      "Epoch 120/1000\n",
      "0s - loss: 62.0854 - val_loss: 50.6114\n",
      "Epoch 121/1000\n",
      "0s - loss: 62.0551 - val_loss: 45.6362\n",
      "Epoch 122/1000\n",
      "0s - loss: 63.6479 - val_loss: 44.3717\n",
      "Epoch 123/1000\n",
      "0s - loss: 60.9967 - val_loss: 46.9744\n",
      "Epoch 124/1000\n",
      "0s - loss: 57.2035 - val_loss: 44.4330\n",
      "Epoch 125/1000\n",
      "0s - loss: 56.9244 - val_loss: 44.5246\n",
      "Epoch 126/1000\n",
      "0s - loss: 56.9935 - val_loss: 44.0118\n",
      "Epoch 127/1000\n",
      "0s - loss: 55.9113 - val_loss: 43.7560\n",
      "Epoch 128/1000\n",
      "0s - loss: 55.2014 - val_loss: 42.3474\n",
      "Epoch 129/1000\n",
      "0s - loss: 55.1030 - val_loss: 42.6942\n",
      "Epoch 130/1000\n",
      "0s - loss: 57.3711 - val_loss: 45.5312\n",
      "Epoch 131/1000\n",
      "0s - loss: 54.8791 - val_loss: 41.4254\n",
      "Epoch 132/1000\n",
      "0s - loss: 56.5617 - val_loss: 42.8268\n",
      "Epoch 133/1000\n",
      "0s - loss: 55.0927 - val_loss: 40.8612\n",
      "Epoch 134/1000\n",
      "0s - loss: 53.8303 - val_loss: 40.6034\n",
      "Epoch 135/1000\n",
      "0s - loss: 53.2116 - val_loss: 40.4556\n",
      "Epoch 136/1000\n",
      "0s - loss: 52.0259 - val_loss: 40.0974\n",
      "Epoch 137/1000\n",
      "0s - loss: 51.5206 - val_loss: 44.9810\n",
      "Epoch 138/1000\n",
      "0s - loss: 51.7043 - val_loss: 39.8432\n",
      "Epoch 139/1000\n",
      "0s - loss: 51.1995 - val_loss: 39.6087\n",
      "Epoch 140/1000\n",
      "0s - loss: 51.4319 - val_loss: 41.9399\n",
      "Epoch 141/1000\n",
      "0s - loss: 51.4475 - val_loss: 39.2047\n",
      "Epoch 142/1000\n",
      "0s - loss: 51.3008 - val_loss: 43.7470\n",
      "Epoch 143/1000\n",
      "0s - loss: 52.6647 - val_loss: 40.3957\n",
      "Epoch 144/1000\n",
      "0s - loss: 50.0021 - val_loss: 38.3514\n",
      "Epoch 145/1000\n",
      "0s - loss: 50.6613 - val_loss: 38.5582\n",
      "Epoch 146/1000\n",
      "0s - loss: 52.3884 - val_loss: 37.9323\n",
      "Epoch 147/1000\n",
      "0s - loss: 51.6568 - val_loss: 37.7600\n",
      "Epoch 148/1000\n",
      "0s - loss: 48.4505 - val_loss: 38.2580\n",
      "Epoch 149/1000\n",
      "0s - loss: 47.9677 - val_loss: 38.1941\n",
      "Epoch 150/1000\n",
      "0s - loss: 47.7968 - val_loss: 37.2814\n",
      "Epoch 151/1000\n",
      "0s - loss: 47.5908 - val_loss: 37.0471\n",
      "Epoch 152/1000\n",
      "0s - loss: 47.4636 - val_loss: 36.8562\n",
      "Epoch 153/1000\n",
      "0s - loss: 46.7298 - val_loss: 36.9440\n",
      "Epoch 154/1000\n",
      "0s - loss: 47.5526 - val_loss: 38.0565\n",
      "Epoch 155/1000\n",
      "0s - loss: 46.5543 - val_loss: 36.3607\n",
      "Epoch 156/1000\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0s - loss: 46.5518 - val_loss: 37.9595\n",
      "Epoch 157/1000\n",
      "0s - loss: 47.1224 - val_loss: 38.1219\n",
      "Epoch 158/1000\n",
      "0s - loss: 47.3709 - val_loss: 37.2413\n",
      "Epoch 159/1000\n",
      "0s - loss: 45.4097 - val_loss: 35.7499\n",
      "Epoch 160/1000\n",
      "0s - loss: 45.2445 - val_loss: 35.7975\n",
      "Epoch 161/1000\n",
      "0s - loss: 44.8730 - val_loss: 35.9459\n",
      "Epoch 162/1000\n",
      "0s - loss: 44.7698 - val_loss: 37.3796\n",
      "Epoch 163/1000\n",
      "0s - loss: 44.7476 - val_loss: 40.9131\n",
      "Epoch 164/1000\n",
      "0s - loss: 45.3167 - val_loss: 40.5059\n",
      "Epoch 165/1000\n",
      "0s - loss: 47.8969 - val_loss: 42.4199\n",
      "Epoch 00164: early stopping\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<keras.callbacks.History at 0x17901796b70>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from sklearn.model_selection import train_test_split\n",
    "import pandas as pd\n",
    "import os\n",
    "import numpy as np\n",
    "from sklearn import metrics\n",
    "from scipy.stats import zscore\n",
    "\n",
    "path = \"./data/\"\n",
    "\n",
    "filename_read = os.path.join(path,\"auto-mpg.csv\")\n",
    "df = pd.read_csv(filename_read,na_values=['NA','?'])\n",
    "\n",
    "cars = df['name']\n",
    "df.drop('name',1,inplace=True)\n",
    "missing_median(df, 'horsepower')\n",
    "x,y = to_xy(df,\"mpg\")\n",
    "# Split into train/test\n",
    "x_train, x_test, y_train, y_test = train_test_split(    \n",
    "    x, y, test_size=0.25, random_state=45)\n",
    "model = Sequential()\n",
    "model.add(Dense(10, input_dim=x.shape[1], activation='relu'))\n",
    "model.add(Dense(10))\n",
    "model.add(Dense(10))\n",
    "model.add(Dense(10))\n",
    "model.add(Dense(1))\n",
    "model.compile(loss='mean_squared_error', optimizer='adam')\n",
    "monitor = EarlyStopping(monitor='val_loss', min_delta=1e-3, patience=5, verbose=1, mode='auto')\n",
    "model.fit(x,y,validation_data=(x_test,y_test),callbacks=[monitor],verbose=2,epochs=1000)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Mean Square Error\n",
    "\n",
    "The mean square error is the sum of the squared differences between the prediction ($\\hat{y}$) and the expected ($y$).  MSE values are not of a particular unit.  If an MSE value has decreased for a model, that is good.  However, beyond this, there is not much more you can determine.  Low MSE values are desired.\n",
    "\n",
    "$ \\text{MSE} = \\frac{1}{n} \\sum_{i=1}^n \\left(\\hat{y}_i - y_i\\right)^2 $\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Final score (MSE): 42.41987609863281\n"
     ]
    }
   ],
   "source": [
    "# Predict\n",
    "pred = model.predict(x_test)\n",
    "\n",
    "# Measure MSE error.  \n",
    "score = metrics.mean_squared_error(pred,y_test)\n",
    "print(\"Final score (MSE): {}\".format(score))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Root Mean Square Error\n",
    "\n",
    "The root mean square (RMSE) is essentially the square root of the MSE.  Because of this, the RMSE error is in the same units as the training data outcome. Low RMSE values are desired.\n",
    "\n",
    "$ \\text{RMSE} = \\sqrt{\\frac{1}{n} \\sum_{i=1}^n \\left(\\hat{y}_i - y_i\\right)^2} $"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Final score (RMSE): 6.513054370880127\n"
     ]
    }
   ],
   "source": [
    "# Measure RMSE error.  RMSE is common for regression.\n",
    "score = np.sqrt(metrics.mean_squared_error(pred,y_test))\n",
    "print(\"Final score (RMSE): {}\".format(score))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Training with Cross-Validation\n",
    "\n",
    "Cross-Validation uses a number of folds, and multiple models, to generate out of sample predictions on the entire dataset.  It is important to note that there will be one model (neural network) for each fold. Each model contributes part of the final out-of-sample prediction.\n",
    "\n",
    "![K-Fold Crossvalidation](https://raw.githubusercontent.com/jeffheaton/t81_558_deep_learning/master/images/class_1_kfold.png \"K-Fold Crossvalidation\")\n",
    "\n",
    "For new data, which is data not present in the training set, predictions from the fold models can be handled in several ways.\n",
    "\n",
    "* Choose the model that had the highest validation score as the final model.\n",
    "* Preset new data to the 5 models and average the result (this is an [enesmble](https://en.wikipedia.org/wiki/Ensemble_learning)).\n",
    "* Retrain a new model (using the same settings as the cross-validation) on the entire dataset.  Train for as many steps, and with the same hidden layer structure.\n",
    "\n",
    "## Regression with Cross-Validation\n",
    "\n",
    "The following code trains the MPG dataset using a 5-fold cross-validation.  The expected performance of a neural network, of the type trained here, would be the score for the generated out-of-sample predictions."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Fold #1\n",
      "Epoch 00255: early stopping\n",
      "Fold score (RMSE): 8.614513397216797\n",
      "Fold #2\n",
      "Epoch 00565: early stopping\n",
      "Fold score (RMSE): 3.9934253692626953\n",
      "Fold #3\n",
      "Epoch 00448: early stopping\n",
      "Fold score (RMSE): 4.960789203643799\n",
      "Fold #4\n",
      "Epoch 00031: early stopping\n",
      "Fold score (RMSE): 15.435755729675293\n",
      "Fold #5\n",
      "Epoch 00406: early stopping\n",
      "Fold score (RMSE): 3.914012908935547\n",
      "Final, out of sample score (RMSE): 8.567543983459473\n"
     ]
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import os\n",
    "import numpy as np\n",
    "from sklearn import metrics\n",
    "from scipy.stats import zscore\n",
    "from sklearn.model_selection import KFold\n",
    "from keras.models import Sequential\n",
    "from keras.layers.core import Dense, Activation\n",
    "\n",
    "path = \"./data/\"\n",
    "\n",
    "filename_read = os.path.join(path,\"auto-mpg.csv\")\n",
    "filename_write = os.path.join(path,\"auto-mpg-out-of-sample.csv\")\n",
    "df = pd.read_csv(filename_read,na_values=['NA','?'])\n",
    "\n",
    "# Shuffle\n",
    "np.random.seed(42)\n",
    "df = df.reindex(np.random.permutation(df.index))\n",
    "df.reset_index(inplace=True, drop=True)\n",
    "\n",
    "# Preprocess\n",
    "cars = df['name']\n",
    "df.drop('name',1,inplace=True)\n",
    "missing_median(df, 'horsepower')\n",
    "\n",
    "# Encode to a 2D matrix for training\n",
    "x,y = to_xy(df,'mpg')\n",
    "\n",
    "# Cross-Validate\n",
    "kf = KFold(5)\n",
    "    \n",
    "oos_y = []\n",
    "oos_pred = []\n",
    "fold = 0\n",
    "for train, test in kf.split(x):\n",
    "    fold+=1\n",
    "    print(\"Fold #{}\".format(fold))\n",
    "        \n",
    "    x_train = x[train]\n",
    "    y_train = y[train]\n",
    "    x_test = x[test]\n",
    "    y_test = y[test]\n",
    "    \n",
    "    model = Sequential()\n",
    "    model.add(Dense(20, input_dim=x.shape[1], activation='relu'))\n",
    "    model.add(Dense(10, activation='relu'))\n",
    "    model.add(Dense(1))\n",
    "    model.compile(loss='mean_squared_error', optimizer='adam')\n",
    "    \n",
    "    monitor = EarlyStopping(monitor='val_loss', min_delta=1e-3, patience=5, verbose=1, mode='auto')\n",
    "    model.fit(x_train,y_train,validation_data=(x_test,y_test),callbacks=[monitor],verbose=0,epochs=1000)\n",
    "    \n",
    "    pred = model.predict(x_test)\n",
    "    \n",
    "    oos_y.append(y_test)\n",
    "    oos_pred.append(pred)        \n",
    "\n",
    "    # Measure this fold's RMSE\n",
    "    score = np.sqrt(metrics.mean_squared_error(pred,y_test))\n",
    "    print(\"Fold score (RMSE): {}\".format(score))\n",
    "\n",
    "\n",
    "# Build the oos prediction list and calculate the error.\n",
    "oos_y = np.concatenate(oos_y)\n",
    "oos_pred = np.concatenate(oos_pred)\n",
    "score = np.sqrt(metrics.mean_squared_error(oos_pred,oos_y))\n",
    "print(\"Final, out of sample score (RMSE): {}\".format(score))    \n",
    "    \n",
    "# Write the cross-validated prediction\n",
    "oos_y = pd.DataFrame(oos_y)\n",
    "oos_pred = pd.DataFrame(oos_pred)\n",
    "oosDF = pd.concat( [df, oos_y, oos_pred],axis=1 )\n",
    "oosDF.to_csv(filename_write,index=False)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Classification with Cross-Validation\n",
    "\n",
    "The following code trains and fits the iris dataset with Cross-Validation.  It also writes out the out of sample (predictions on the test set) results."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Fold #1\n",
      "Epoch 00332: early stopping\n",
      "Fold score (accuracy): 1.0\n",
      "Fold #2\n",
      "Epoch 00201: early stopping\n",
      "Fold score (accuracy): 0.9666666666666667\n",
      "Fold #3\n",
      "Epoch 00411: early stopping\n",
      "Fold score (accuracy): 0.9666666666666667\n",
      "Fold #4\n",
      "Epoch 00164: early stopping\n",
      "Fold score (accuracy): 0.9333333333333333\n",
      "Fold #5\n",
      "Epoch 00277: early stopping\n",
      "Fold score (accuracy): 1.0\n",
      "Final score (accuracy): 0.9733333333333334\n"
     ]
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import os\n",
    "import numpy as np\n",
    "from sklearn import metrics\n",
    "from scipy.stats import zscore\n",
    "from sklearn.model_selection import KFold\n",
    "from keras.models import Sequential\n",
    "from keras.layers.core import Dense, Activation\n",
    "\n",
    "path = \"./data/\"\n",
    "\n",
    "filename_read = os.path.join(path,\"iris.csv\")\n",
    "filename_write = os.path.join(path,\"iris-out-of-sample.csv\")\n",
    "df = pd.read_csv(filename_read,na_values=['NA','?'])\n",
    "\n",
    "# Shuffle\n",
    "np.random.seed(42)\n",
    "df = df.reindex(np.random.permutation(df.index))\n",
    "df.reset_index(inplace=True, drop=True)\n",
    "\n",
    "# Encode to a 2D matrix for training\n",
    "species = encode_text_index(df,\"species\")\n",
    "x,y = to_xy(df,\"species\")\n",
    "\n",
    "# Cross-validate\n",
    "kf = KFold(5)\n",
    "    \n",
    "oos_y = []\n",
    "oos_pred = []\n",
    "fold = 0\n",
    "\n",
    "for train, test in kf.split(x):\n",
    "    fold+=1\n",
    "    print(\"Fold #{}\".format(fold))\n",
    "        \n",
    "    x_train = x[train]\n",
    "    y_train = y[train]\n",
    "    x_test = x[test]\n",
    "    y_test = y[test]\n",
    "    \n",
    "    model = Sequential()\n",
    "    model.add(Dense(50, input_dim=x.shape[1], activation='relu')) # Hidden 1\n",
    "    model.add(Dense(25, activation='relu')) # Hidden 2\n",
    "    model.add(Dense(y.shape[1],activation='softmax')) # Output\n",
    "    model.compile(loss='categorical_crossentropy', optimizer='adam')\n",
    "    monitor = EarlyStopping(monitor='val_loss', min_delta=1e-3, patience=25, verbose=1, mode='auto')\n",
    "\n",
    "    model.fit(x,y,validation_data=(x_test,y_test),callbacks=[monitor],verbose=0,epochs=1000)\n",
    "    \n",
    "    pred = model.predict(x_test)\n",
    "    \n",
    "    oos_y.append(y_test)\n",
    "    pred = np.argmax(pred,axis=1) # raw probabilities to chosen class (highest probability)\n",
    "    oos_pred.append(pred)        \n",
    "\n",
    "    # Measure this fold's accuracy\n",
    "    y_compare = np.argmax(y_test,axis=1) # For accuracy calculation\n",
    "    score = metrics.accuracy_score(y_compare, pred)\n",
    "    print(\"Fold score (accuracy): {}\".format(score))\n",
    "\n",
    "\n",
    "# Build the oos prediction list and calculate the error.\n",
    "oos_y = np.concatenate(oos_y)\n",
    "oos_pred = np.concatenate(oos_pred)\n",
    "oos_y_compare = np.argmax(oos_y,axis=1) # For accuracy calculation\n",
    "\n",
    "score = metrics.accuracy_score(oos_y_compare, oos_pred)\n",
    "print(\"Final score (accuracy): {}\".format(score))    \n",
    "    \n",
    "# Write the cross-validated prediction\n",
    "oos_y = pd.DataFrame(oos_y)\n",
    "oos_pred = pd.DataFrame(oos_pred)\n",
    "oosDF = pd.concat( [df, oos_y, oos_pred],axis=1 )\n",
    "oosDF.to_csv(filename_write,index=False)\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Training with both a Cross-Validation and a Holdout Set\n",
    "\n",
    "If you have a considerable amount of data, it is always valuable to set aside a holdout set before you cross-validate.  This hold out set will be the final evaluation before you make use of your model for its real-world use.\n",
    "\n",
    "![Cross Validation and a Holdout Set](https://raw.githubusercontent.com/jeffheaton/t81_558_deep_learning/master/images/class_3_hold_train_val.png \"Cross Validation and a Holdout Set\")\n",
    "\n",
    "The following program makes use of a holdout set, and then still cross-validates.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Fold #1\n",
      "Epoch 00059: early stopping\n",
      "Fold score (RMSE): 15.772638320922852\n",
      "Fold #2\n",
      "Epoch 00220: early stopping\n",
      "Fold score (RMSE): 3.587620496749878\n",
      "Fold #3\n",
      "Epoch 00402: early stopping\n",
      "Fold score (RMSE): 4.99105167388916\n",
      "Fold #4\n",
      "Epoch 00675: early stopping\n",
      "Fold score (RMSE): 3.76237154006958\n",
      "Fold #5\n",
      "Epoch 00336: early stopping\n",
      "Fold score (RMSE): 5.120087623596191\n",
      "\n",
      "Cross-validated score (RMSE): 8.101740837097168\n",
      "Holdout score (RMSE): 3.4944186210632324\n"
     ]
    }
   ],
   "source": [
    "from sklearn.model_selection import train_test_split\n",
    "import pandas as pd\n",
    "import os\n",
    "import numpy as np\n",
    "from sklearn import metrics\n",
    "from scipy.stats import zscore\n",
    "from sklearn.model_selection import KFold\n",
    "from keras.callbacks import EarlyStopping\n",
    "\n",
    "path = \"./data/\"\n",
    "\n",
    "filename_read = os.path.join(path,\"auto-mpg.csv\")\n",
    "filename_write = os.path.join(path,\"auto-mpg-holdout.csv\")\n",
    "df = pd.read_csv(filename_read,na_values=['NA','?'])\n",
    "\n",
    "# create feature vector\n",
    "missing_median(df, 'horsepower')\n",
    "df.drop('name',1,inplace=True)\n",
    "encode_text_dummy(df, 'origin')\n",
    "\n",
    "# Shuffle\n",
    "np.random.seed(42)\n",
    "df = df.reindex(np.random.permutation(df.index))\n",
    "df.reset_index(inplace=True, drop=True)\n",
    "\n",
    "# Encode to a 2D matrix for training\n",
    "x,y = to_xy(df,'mpg')\n",
    "\n",
    "\n",
    "# Keep a 10% holdout\n",
    "x_main, x_holdout, y_main, y_holdout = train_test_split(    \n",
    "    x, y, test_size=0.10) \n",
    "\n",
    "\n",
    "# Cross-validate\n",
    "kf = KFold(5)\n",
    "    \n",
    "oos_y = []\n",
    "oos_pred = []\n",
    "fold = 0\n",
    "for train, test in kf.split(x_main):        \n",
    "    fold+=1\n",
    "    print(\"Fold #{}\".format(fold))\n",
    "        \n",
    "    x_train = x_main[train]\n",
    "    y_train = y_main[train]\n",
    "    x_test = x_main[test]\n",
    "    y_test = y_main[test]\n",
    "    \n",
    "    model = Sequential()\n",
    "    model.add(Dense(20, input_dim=x.shape[1], activation='relu'))\n",
    "    model.add(Dense(5, activation='relu'))\n",
    "    model.add(Dense(1))\n",
    "    model.compile(loss='mean_squared_error', optimizer='adam')\n",
    "    \n",
    "    monitor = EarlyStopping(monitor='val_loss', min_delta=1e-3, patience=5, verbose=1, mode='auto')\n",
    "    model.fit(x_train,y_train,validation_data=(x_test,y_test),callbacks=[monitor],verbose=0,epochs=1000)\n",
    "    \n",
    "    pred = model.predict(x_test)\n",
    "    \n",
    "    oos_y.append(y_test)\n",
    "    oos_pred.append(pred) \n",
    "\n",
    "    # Measure accuracy\n",
    "    score = np.sqrt(metrics.mean_squared_error(pred,y_test))\n",
    "    print(\"Fold score (RMSE): {}\".format(score))\n",
    "\n",
    "\n",
    "# Build the oos prediction list and calculate the error.\n",
    "oos_y = np.concatenate(oos_y)\n",
    "oos_pred = np.concatenate(oos_pred)\n",
    "score = np.sqrt(metrics.mean_squared_error(oos_pred,oos_y))\n",
    "print()\n",
    "print(\"Cross-validated score (RMSE): {}\".format(score))    \n",
    "    \n",
    "# Write the cross-validated prediction (from the last neural network)\n",
    "holdout_pred = model.predict(x_holdout)\n",
    "\n",
    "score = np.sqrt(metrics.mean_squared_error(holdout_pred,y_holdout))\n",
    "print(\"Holdout score (RMSE): {}\".format(score))    \n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Scikit-Learn Versions: model_selection vs cross_validation\n",
    "\n",
    "Scikit-Learn changed a bit in how cross-validation is handled.  Both versions still work, but you should use the **sklearn.model_selection** import, rather than **sklearn.cross_validation**.  The following shows both the new and old forms of cross-validation.  All examples from this class will use the newer form.\n",
    "\n",
    "The following two sections show both forms:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Fold #1: train=318, test=80\n",
      "Fold #2: train=318, test=80\n",
      "Fold #3: train=318, test=80\n",
      "Fold #4: train=319, test=79\n",
      "Fold #5: train=319, test=79\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\jeffh\\Anaconda3\\lib\\site-packages\\sklearn\\cross_validation.py:41: DeprecationWarning: This module was deprecated in version 0.18 in favor of the model_selection module into which all the refactored classes and functions are moved. Also note that the interface of the new CV iterators are different from that of this module. This module will be removed in 0.20.\n",
      "  \"This module will be removed in 0.20.\", DeprecationWarning)\n"
     ]
    }
   ],
   "source": [
    "# Older scikit-learn syntax for splits/cross-validation\n",
    "# Still valid, but going away.  Do not use.\n",
    "# (Note the red box warning below)\n",
    "\n",
    "from sklearn.cross_validation import train_test_split\n",
    "from sklearn.cross_validation import KFold\n",
    "\n",
    "path = \"./data/\"\n",
    "\n",
    "filename_read = os.path.join(path,\"auto-mpg.csv\")\n",
    "df = pd.read_csv(filename_read,na_values=['NA','?'])\n",
    "\n",
    "kf = KFold(len(df), n_folds=5)\n",
    "\n",
    "fold = 0\n",
    "for train, test in kf:        \n",
    "    fold+=1\n",
    "    print(\"Fold #{}: train={}, test={}\".format(fold,len(train),len(test)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Fold #1: train=318, test=80\n",
      "Fold #2: train=318, test=80\n",
      "Fold #3: train=318, test=80\n",
      "Fold #4: train=319, test=79\n",
      "Fold #5: train=319, test=79\n"
     ]
    }
   ],
   "source": [
    "# Newer scikit-learn syntax for splits/cross-validation\n",
    "# Use this method (as shown above)\n",
    "\n",
    "from sklearn.model_selection import train_test_split\n",
    "from sklearn.model_selection import KFold\n",
    "\n",
    "path = \"./data/\"\n",
    "\n",
    "filename_read = os.path.join(path,\"auto-mpg.csv\")\n",
    "df = pd.read_csv(filename_read,na_values=['NA','?'])\n",
    "\n",
    "kf = KFold(5)\n",
    "\n",
    "fold = 0\n",
    "for train, test in kf.split(df):        \n",
    "    fold+=1\n",
    "    print(\"Fold #{}: train={}, test={}\".format(fold,len(train),len(test)))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# How Kaggle Competitions are Scored\n",
    "\n",
    "[Kaggle](https://www.kaggle.com/) is a platform for competitive data science. Competitions are posted onto Kaggle by companies seeking the best model for their data.  Competing in a Kaggle competition is quite a bit of work, I've [competed in one Kaggle competition](https://www.kaggle.com/jeffheaton).   \n",
    "\n",
    "Kaggle awards \"tiers\", such as:\n",
    "\n",
    "* Kaggle Grandmaster\n",
    "* Kaggle Master\n",
    "* Kaggle Expert\n",
    "\n",
    "Your [tier](https://www.kaggle.com/progression) is based on your performance in past competitions.  \n",
    "\n",
    "To compete in Kaggle you simply provide predictions for a dataset that they post.  You do not need to submit any code.  Your prediction output will place you onto the [leaderboard of a competition](https://www.kaggle.com/c/otto-group-product-classification-challenge/leaderboard/public).\n",
    "\n",
    "![How Kaggle Competitions are Scored](https://raw.githubusercontent.com/jeffheaton/t81_558_deep_learning/master/images/class_3_kaggle.png \"How Kaggle Competitions are Scored\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "An original dataset is sent to Kaggle by the company.  From this dataset, Kaggle posts public data that includes \"train\" and \"test.  For the \"train\" data, the outcomes (y) are provided.  For the test data, no outcomes are provided.  Your submission file contains your predictions for the \"test data\".  When you submit your results, Kaggle will calculate a score on part of your prediction data.  They do not publish want part of the submission data are used for the public and private leaderboard scores (this is a secret to prevent overfitting).  While the competition is still running, Kaggle publishes the public leaderboard ranks.  Once the competition ends, the private leaderboard is revealed to designate the true winners.  Due to overfitting, there is sometimes an upset in positions when the final private leaderboard is revealed.  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Managing Hyperparameters\n",
    "\n",
    "There are many different settings that you can use for a neural network.  These can affect performance.  The following code changes some of these, beyond their default values:\n",
    "\n",
    "* **activation:** relu, sigmoid, tanh\n",
    "* Layers/Neuron Counts\n",
    "* **optimizer:** adam, sgd, rmsprop, and [others](https://keras.io/optimizers/)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Epoch 00065: early stopping\n",
      "Score (RMSE): 3.935321807861328\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "from matplotlib.pyplot import figure, show\n",
    "from sklearn.model_selection import train_test_split\n",
    "import pandas as pd\n",
    "import os\n",
    "import numpy as np\n",
    "from sklearn import metrics\n",
    "from scipy.stats import zscore\n",
    "import tensorflow as tf\n",
    "\n",
    "path = \"./data/\"\n",
    "preprocess = False\n",
    "\n",
    "filename_read = os.path.join(path,\"auto-mpg.csv\")\n",
    "df = pd.read_csv(filename_read,na_values=['NA','?'])\n",
    "\n",
    "# create feature vector\n",
    "missing_median(df, 'horsepower')\n",
    "encode_text_dummy(df, 'origin')\n",
    "df.drop('name',1,inplace=True)\n",
    "if preprocess:\n",
    "    encode_numeric_zscore(df, 'horsepower')\n",
    "    encode_numeric_zscore(df, 'weight')\n",
    "    encode_numeric_zscore(df, 'cylinders')\n",
    "    encode_numeric_zscore(df, 'displacement')\n",
    "    encode_numeric_zscore(df, 'acceleration')\n",
    "\n",
    "# Encode to a 2D matrix for training\n",
    "x,y = to_xy(df,'mpg')\n",
    "\n",
    "# Split into train/test\n",
    "x_train, x_test, y_train, y_test = train_test_split(\n",
    "    x, y, test_size=0.20, random_state=42)\n",
    "\n",
    "model = Sequential()\n",
    "model.add(Dense(100, input_dim=x.shape[1], activation='relu'))\n",
    "model.add(Dense(50, activation='relu'))\n",
    "model.add(Dense(25, activation='relu'))\n",
    "model.add(Dense(1))\n",
    "model.compile(loss='mean_squared_error', optimizer='adam')\n",
    "monitor = EarlyStopping(monitor='val_loss', min_delta=1e-3, patience=5, verbose=1, mode='auto')\n",
    "model.fit(x,y,validation_data=(x_test,y_test),callbacks=[monitor],verbose=0,epochs=1000)\n",
    "\n",
    "# Predict and measure RMSE\n",
    "pred = model.predict(x_test)\n",
    "score = np.sqrt(metrics.mean_squared_error(pred,y_test))\n",
    "print(\"Score (RMSE): {}\".format(score))\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "# Module 4 Assignment\n",
    "\n",
    "You can find the first assignmeht here: [assignment 4](https://github.com/jeffheaton/t81_558_deep_learning/blob/master/assignments/assignment_yourname_class1.ipynb)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}