notebooks/sklearn.ipynb

{
 "metadata": {
  "name": "",
  "signature": "sha256:29ddb8a14db89c6d65232c988ad3f9bde13b707aa10b0b4d1f7bc5af9ebc6a77"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%load_ext load_style\n",
      "%load_style talk.css\n",
      "from IPython.display import Image\n",
      "from talktools import website"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "html": [
        "<style>\n",
        "\n",
        ".rendered_html {\n",
        "    font-family: \"proxima-nova\", helvetica;\n",
        "    font-size: 130%;\n",
        "    line-height: 1.5;\n",
        "}\n",
        "\n",
        ".rendered_html h1 {\n",
        "    margin: 0.25em 0em 0.5em;\n",
        "    color: #015C9C;\n",
        "    text-align: center;\n",
        "    line-height: 1.2; \n",
        "    page-break-before: always;\n",
        "}\n",
        "\n",
        ".rendered_html h2 {\n",
        "    margin: 1.1em 0em 0.5em;\n",
        "    color: #26465D;\n",
        "    line-height: 1.2;\n",
        "}\n",
        "\n",
        ".rendered_html h3 {\n",
        "    margin: 1.1em 0em 0.5em;\n",
        "    color: #002845;\n",
        "    line-height: 1.2;\n",
        "}\n",
        "\n",
        ".rendered_html li {\n",
        "    line-height: 1.5; \n",
        "}\n",
        "\n",
        "/*.prompt {\n",
        "    font-size: 120%; \n",
        "}*/\n",
        "\n",
        ".CodeMirror-lines {\n",
        "    font-size: 110%; \n",
        "}\n",
        "\n",
        "/*.output_area {\n",
        "    font-size: 120%; \n",
        "}*/\n",
        "\n",
        "/*#notebook {\n",
        "    background-image: url('files/images/witewall_3.png');\n",
        "}*/\n",
        "\n",
        "h1.bigtitle {\n",
        "    margin: 4cm 1cm 4cm 1cm;\n",
        "    font-size: 300%;\n",
        "}\n",
        "\n",
        "h3.point {\n",
        "    font-size: 200%;\n",
        "    text-align: center;\n",
        "    margin: 2em 0em 2em 0em;\n",
        "    #26465D\n",
        "}\n",
        "\n",
        ".logo {\n",
        "    margin: 20px 0 20px 0;\n",
        "}\n",
        "\n",
        "a.anchor-link {\n",
        "    display: none;\n",
        "}\n",
        "\n",
        "h1.title { \n",
        "    font-size: 250%;\n",
        "}\n",
        "\n",
        "</style>"
       ],
       "metadata": {},
       "output_type": "display_data",
       "text": [
        "<IPython.core.display.HTML at 0x1044eb950>"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Machine Learning with scikit-learn"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Broadly speaking, Machine Learning is a field, related to Artificial Intelligence (AI), concerned about developping algorithms that 'learn from data', i.e. automatically adjust their performance from exposure to information encoded in data. This learning is achieved via **tunable parameters** that are automatically adjusted according to performance criteria.\n",
      "\n",
      "There are two major classes of ML (actually 3 with [**reinforcement learning**](http://en.wikipedia.org/wiki/Reinforcement_learning)): \n",
      "\n",
      "[**Supervised learning**](http://en.wikipedia.org/wiki/Supervised_learning)\n",
      ": Algorithms which learn from a training set of *labeled* examples to generalize to the set of all possible inputs. \n",
      "\n",
      "There are two classes of supervised learning algorithms: \n",
      "\n",
      "1. [**classification**](http://en.wikipedia.org/wiki/Classification_(machine_learning): when the label is encoded into a discrete, categorical variable (a *class*, a *label*). Example of a classification algorithm is the [Support Vector Machine](http://en.wikipedia.org/wiki/Support_vector_machine)\n",
      "\n",
      "2. [**regression**](http://en.wikipedia.org/wiki/Regression_analysis): when the label is encoded into a continuous variable. Example of a regression is the simple [linear regression](http://en.wikipedia.org/wiki/Linear_regression_model). \n",
      "\n",
      "[**Unsupervised learning**](http://en.wikipedia.org/wiki/Unsupervised_learning)\n",
      ": Algorithms which learn from a training set of *unlabeled* examples, using the features of the inputs to categorize inputs together according to some statistical criteria. \n",
      "\n",
      "One can also divide Unsupervised learning algorithms into: \n",
      "\n",
      "1. [**Dimensionality reduction**](http://en.wikipedia.org/wiki/Dimensionality_reduction): learning a more compact  representation (i.e. reducing the dimensions) of the data. One example is [Principal Component Analysis](http://en.wikipedia.org/wiki/Principal_component_analysis) (or Empirical Orthogonal Functions decomposition) \n",
      "2. [**Clustering**](http://en.wikipedia.org/wiki/Cluster_analysis): separating the data into clusters. One example of widely used algorithm is the [k-means](http://en.wikipedia.org/wiki/K-means) clustering method."
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "### Some vocabulary \n",
      "\n",
      "+ ** Instance **: sample, observation\n",
      "+ ** Feature **: explanatory variable, independent variable, predictor\n",
      "+ ** target **: dependent variable, predictand, can be categorical (*label*, *class*) or continuous."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import numpy as np\n",
      "import pandas as pd\n",
      "from matplotlib import pyplot as plt\n",
      "from IPython.display import Image, HTML\n",
      "%matplotlib inline"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 2
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "Image(url='http://scikit-learn.org/stable/_static/ml_map.png', width=900)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "html": [
        "<img src=\"http://scikit-learn.org/stable/_static/ml_map.png\" width=\"900\"/>"
       ],
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 3,
       "text": [
        "<IPython.core.display.Image at 0x1044ebbd0>"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "The [scikit-learn](http://scikit-learn.org/stable/) package is an open-source library that provides a robust set of machine learning algorithms for Python. It is built upon the core Python scientific stack (*i.e.* NumPy, SciPy, Cython), and has a simple, consistent API, making it useful for a wide range of statistical learning applications."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#website('http://scikit-learn.org/stable/', width=1000)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 4
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "### Representing Data in `scikit-learn`\n",
      "\n",
      "Most machine learning algorithms implemented in scikit-learn expect data to be stored in a\n",
      "**two-dimensional array or matrix**.  The arrays can be\n",
      "either ``numpy`` arrays, or in some cases ``scipy.sparse`` matrices.\n",
      "The size of the array is expected to be `[n_samples, n_features]`\n",
      "\n",
      "- **n_samples:**   The number of samples: each sample is an item to process (e.g. classify).\n",
      "  A sample can be a document, a picture, a sound, a video, an astronomical object,\n",
      "  a row in database or CSV file,\n",
      "  or whatever you can describe with a fixed set of quantitative traits.\n",
      "- **n_features:**  The number of features or distinct traits that can be used to describe each\n",
      "  item in a quantitative manner.  Features are generally real-valued, but may be boolean or\n",
      "  discrete-valued in some cases.\n",
      "\n",
      "The number of features must be fixed in advance. However it can be very high dimensional\n",
      "(e.g. millions of features) with most of them being zeros for a given sample. This is a case\n",
      "where `scipy.sparse` matrices can be useful, in that they are\n",
      "much more memory-efficient than numpy arrays."
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "$$\n",
      "{\\rm feature~matrix:~~~} {\\bf X}~=~\\left[\n",
      "\\begin{matrix}\n",
      "x_{11} & x_{12} & \\cdots & x_{1D}\\\\\n",
      "x_{21} & x_{22} & \\cdots & x_{2D}\\\\\n",
      "x_{31} & x_{32} & \\cdots & x_{3D}\\\\\n",
      "\\vdots & \\vdots & \\ddots & \\vdots\\\\\n",
      "\\vdots & \\vdots & \\ddots & \\vdots\\\\\n",
      "x_{N1} & x_{N2} & \\cdots & x_{ND}\\\\\n",
      "\\end{matrix}\n",
      "\\right]\n",
      "$$\n",
      "\n",
      "$$\n",
      "{\\rm label~vector:~~~} {\\bf y}~=~ [y_1, y_2, y_3, \\cdots y_N]\n",
      "$$\n",
      "\n",
      "Here there are $N$ samples and $D$ features.\n",
      "Several example datasets are available in the ``sklearn.datasets`` module"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "### The `scikit-learn` interface\n",
      "\n",
      "A great feature about scikit-learn is the consistant API (Application Programmer Interface) which means that once you have learned how one particular algorithm is implemented in scikit-learn, using another algorithm will look very familiar. \n",
      "\n",
      "All objects within scikit-learn share a uniform common basic API consisting of three complementary interfaces: an **estimator** interface for building and \ufb01tting models, a **predictor** interface for making predictions and a **transformer** interface for converting data.\n",
      "\n",
      "The **estimator** interface is at the core of the library. It de\ufb01nes instantiation mechanisms of objects and exposes a fit method for learning a model from training data. All supervised and unsupervised learning algorithms (*e.g.*, for classi\ufb01cation, regression or clustering) are o\ufb00ered as objects implementing this interface. Machine learning tasks like feature extraction, feature selection or dimensionality reduction are also provided as estimators.\n",
      "\n",
      "Scikit-learn strives to have a uniform interface across all methods; given a scikit-learn *estimator*\n",
      "object named `model`, the following methods are available:\n",
      "\n",
      "- Available in **all Estimators**\n",
      "  + `model.fit()` : fit training data. For supervised learning applications,\n",
      "    this accepts two arguments: the data `X` and the labels `y` (e.g. `model.fit(X, y)`).\n",
      "    For unsupervised learning applications, this accepts only a single argument,\n",
      "    the data `X` (e.g. `model.fit(X)`).\n",
      "- Available in **supervised estimators**\n",
      "  + `model.predict()` : given a trained model, predict the label of a new set of data.\n",
      "    This method accepts one argument, the new data `X_new` (e.g. `model.predict(X_new)`),\n",
      "    and returns the learned label for each object in the array.\n",
      "  + `model.predict_proba()` : For classification problems, some estimators also provide\n",
      "    this method, which returns the probability that a new observation has each categorical label.\n",
      "    In this case, the label with the highest probability is returned by `model.predict()`.\n",
      "  + `model.score()` : for classification or regression problems, most (all?) estimators implement\n",
      "    a score method.  Scores are between 0 and 1, with a larger score indicating a better fit.\n",
      "- Available in **unsupervised estimators**\n",
      "  + `model.transform()` : given an unsupervised model, transform new data into the new basis.\n",
      "    This also accepts one argument `X_new`, and returns the new representation of the data based\n",
      "    on the unsupervised model.\n",
      "  + `model.fit_transform()` : some estimators implement this method,\n",
      "    which more efficiently performs a fit and a transform on the same input data.\n",
      "    \n",
      "The **predictor** interface extends the notion of an estimator by adding a `predict` method that takes an array `X_test` and produces predictions based on the learned parameters of the estimator. In the case of supervised learning estimators, this method typically returns the predicted labels (for *classification*) or values (*regression) computed by the model. Some unsupervised learning estimators may also implement the predict interface, such as k-means, where the predicted values are the cluster labels.\n",
      "\n",
      "Since it is common to modify or \ufb01lter data before feeding it to a learning algorithm, some estimators in the library implement a **transformer** interface which de\ufb01nes a transform method. It takes as input some new data `X_test` and yields as output a transformed version. Preprocessing, feature selection, feature extraction and dimensionality reduction algorithms are all provided as transformers within the library."
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "What we're going to do during this session is give an example of **supervised** learning, and more specifically we're going to see how to solve a **classification** problem in scikit-learn, with a focus on how one evaluates the performance of a model. \n",
      "\n",
      "We're going to use a dataset that comes with scikit-learn, which consists in representation of hand-written digits (8 x 8 pixels normalized images) with the associated label (the correct digit)\n",
      "\n",
      "This example is treated in a more comprehensive manner by [Olivier Grisel](http://ogrisel.com/) (see his notebooks [here](https://github.com/ogrisel/parallel_ml_tutorial))"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn.datasets import load_digits\n",
      "digits = load_digits()\n",
      "print(digits.DESCR)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " Optical Recognition of Handwritten Digits Data Set\n",
        "\n",
        "Notes\n",
        "-----\n",
        "Data Set Characteristics:\n",
        "    :Number of Instances: 5620\n",
        "    :Number of Attributes: 64\n",
        "    :Attribute Information: 8x8 image of integer pixels in the range 0..16.\n",
        "    :Missing Attribute Values: None\n",
        "    :Creator: E. Alpaydin (alpaydin '@' boun.edu.tr)\n",
        "    :Date: July; 1998\n",
        "\n",
        "This is a copy of the test set of the UCI ML hand-written digits datasets\n",
        "http://archive.ics.uci.edu/ml/datasets/Optical+Recognition+of+Handwritten+Digits\n",
        "\n",
        "The data set contains images of hand-written digits: 10 classes where\n",
        "each class refers to a digit.\n",
        "\n",
        "Preprocessing programs made available by NIST were used to extract\n",
        "normalized bitmaps of handwritten digits from a preprinted form. From a\n",
        "total of 43 people, 30 contributed to the training set and different 13\n",
        "to the test set. 32x32 bitmaps are divided into nonoverlapping blocks of\n",
        "4x4 and the number of on pixels are counted in each block. This generates\n",
        "an input matrix of 8x8 where each element is an integer in the range\n",
        "0..16. This reduces dimensionality and gives invariance to small\n",
        "distortions.\n",
        "\n",
        "For info on NIST preprocessing routines, see M. D. Garris, J. L. Blue, G.\n",
        "T. Candela, D. L. Dimmick, J. Geist, P. J. Grother, S. A. Janet, and C.\n",
        "L. Wilson, NIST Form-Based Handprint Recognition System, NISTIR 5469,\n",
        "1994.\n",
        "\n",
        "References\n",
        "----------\n",
        "  - C. Kaynak (1995) Methods of Combining Multiple Classifiers and Their\n",
        "    Applications to Handwritten Digit Recognition, MSc Thesis, Institute of\n",
        "    Graduate Studies in Science and Engineering, Bogazici University.\n",
        "  - E. Alpaydin, C. Kaynak (1998) Cascading Classifiers, Kybernetika.\n",
        "  - Ken Tang and Ponnuthurai N. Suganthan and Xi Yao and A. Kai Qin.\n",
        "    Linear dimensionalityreduction using relevance weighted LDA. School of\n",
        "    Electrical and Electronic Engineering Nanyang Technological University.\n",
        "    2005.\n",
        "  - Claudio Gentile. A New Approximate Maximal Margin Classification\n",
        "    Algorithm. NIPS. 2000.\n",
        "\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "X, y = digits.data, digits.target\n",
      "\n",
      "print(\"data shape: %r, target shape: %r\" % (X.shape, y.shape))\n",
      "print(\"labels: %r\" % list(np.unique(y)))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "data shape: (1797, 64), target shape: (1797,)\n",
        "labels: [0, 1, 2, 3, 4, 5, 6, 7, 8, 9]\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "def plot_gallery(data, labels, shape, interpolation='nearest'):\n",
      "    f,ax = plt.subplots(1,5,figsize=(16,5))\n",
      "    for i in range(data.shape[0]):\n",
      "        ax[i].imshow(data[i].reshape(shape), interpolation=interpolation, cmap=plt.cm.gray_r)\n",
      "        ax[i].set_title(labels[i])\n",
      "        ax[i].set_xticks(()), ax[i].set_yticks(())"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 7
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "subsample = np.random.permutation(X.shape[0])[:5]\n",
      "images = X[subsample]\n",
      "labels = ['True label: %d' % l for l in y[subsample]]\n",
      "plot_gallery(images, labels, shape=(8, 8))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAA4sAAAC3CAYAAAC7Q8soAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADilJREFUeJzt3W1sZXldB/Dvb3cCK09bIkRdlOmoUQLoNBFNjBoaMbww\n6HZXjRpjtpjIGzU7GEWNCdslGH1hSOeNxoeYboyC4rqd+Eaj0Q7EhOiS7RBIgBC2QyAYQejysBBZ\n9u+L3kJn/tuxM/Of9vb280km83DP/Z5z5/567vmee3pbrbUAAADAXrcd9QYAAAAwfZRFAAAAOsoi\nAAAAHWURAACAjrIIAABAR1kEAACgoyzegKparqp3H3DZlar6yxtczw3fF26GGWfWmXFmnRlnlpnv\nw3Osy2JVfaGqPj/59XRVPbnn7z9/1Ns3cTM/yPK67ltVr6mqD1bVF6vqX6vqpTexbqaAGb+SGZ89\nZvxKVfWcqvqjqvpUVW1X1cWbWDdTwIx/XVW9vKoerarPTOb736vqh29i3Rwx8/11szrfx7osttae\n11p7fmvt+UkuJ3nd7t9ba2/fXa6qTh3dVqYO475V9aIkDyf53SQvTPJokr+5iXUzBcz4ngXN+Ewy\n450/TTKX5GXZmfNzN7FupoAZv8InkvxMkm/Mzny/I8nf3cS6OWLm+wozOd/Huizup6oWq+rjVfWm\nqvpkkr+oqvuufrt6cgbk2yd/fnZV/WFVXa6q/6qqP66qOw64vvNV9bGqemJyRmHvWYSW5I6qekdV\nfa6q3ltV37vnvndV1cNV9d9V9dGq+rUbfNj3Jnl/a+3h1tr/JllJcraqvusG85hiZtyMz7qTOONV\n9bIkP5HkDa21/2k7HruRLKbfSZzx1toTrbXHW2stye1Jnk7yyRvJYrqZ79mZ75ksixPflJ1W/9Ik\nb8j/f2bgD5J8Z5Kzk99fkuTNB1zXf0zu98Ikf53knVX1rMltleTuJH+75/b1qrq9qm5L8g9JHkty\nV5LXJDlXVa99ppVU1aWq+rl9tuEVSS7t/qW19mSSjyR55QEfA8ePGTfjs+6kzfgPZOfM/Ftq5zLU\n91XVvQfcfo6nkzbju8tsJ/lSkjcl+ekDbj/Hj/megfme5bL4dJIHWmtfaa19+VoLVlUl+eUkv95a\n226tfSHJ7ye55jDsaq39VWvts621p1trb0vy7CTfvWeRR1trf99a+2qStyW5I8kPJvn+JC9qrb21\ntfZUa+3xJH++33pba2dba+/YZzOem+RzV/3b55I87yCPgWPJjJvxWXfSZvxbs3PyYzvJtyT51SQP\nTd5xZDadtBnfXWYuyZ3ZuUzvnZPHxuwx3zMw30d5/fCt9qnJpWoH8eIkz0ny3j3PZ+WAZbqqfiPJ\nL2XnjERL8oIkL9qzyMd3/9Baa1X18T3L3lVVn92z7O1J3nXA7d7rC5P17nVnks/fQBbHgxk347Pu\npM34l5J8JclbW2tPJ3lXVf1bktcm+eAN5DH9TtqMf01r7cmq+u0kv5Lke5K872bymErmewbme5bL\n4tWfXvTF7AxhkqSqvnnPbZ/Ozov0y1tr13VtcVX9SJLfTPKjrbUPTP7tM7nyrfZv27P8bdk5e/yJ\nJF9N8nhrbb/vubqeT2/6QJL79qznuUm+Y/LvzCYzbsZn3Umb8d2DiavPQt/MJ/kx3U7ajF/t9uyU\ngSdvIoPpZb5nYL5n+TLUq11K8oqqOjv5ZtmV3RsmZ3D/LMlqVb04SarqJftdr3yV5yd5Ksmnq+pZ\nVfXm9O9+fF9V3VM7nwR1LsmXk7wnyX8m+XztfPPvN0yunX5lVb1qcr/redv6kSSvrKp7J4/vgSSb\nrbUPX0cGx5sZZ9bN+oxfTPKxJL9TVaeq6oeSLCb5p+vI4Hib6Rmvqh+rqoVJxguyczngh1prHzlo\nBsea+T6GZrksXnEmYHJA+ZYk/5LkQ0nefdUyv5WdD8t4T1U9keSfk1zrLMPuff9x8uvDSbayc1bk\nY1ctu57kZ5N8JskvJLm3tfbVyXXTr0uykOSjST6VnY9Nf8Ge+35tG6vq/bXPz6xprX06yU8l+b3J\nel6VA17nzbFlxs34rDtpM/5Udj6E4cez832Lf5LkF50QmWknasaz82Nh3p6d+f5Qdi49/Ml9luX4\nM98zMN/VmqtbAAAAuNIsv7MIAADADVIWAQAA6CiLAAAAdJRFAAAAOtf8OYtV5dNvOBSttev5ePlh\nzDiHxYwz645ixs03h8U+nFm334xfsyxO7jh+awbY2NgYlrWysjIsa2tra1hWkqyurg7LWlpaGpY1\nUtWR7H+/ZlpnfORz/8Y3vnFY1p133jksK0nW19eHZS0uLg7LGsmM33pra2vDskZ+7SVjt21hYWFY\n1khHOePTOt+bm5vDska+fs/NzQ3LSsY+zmk1K/vw7e3tITm7Ru4rRx4/j+wIydhji5GvByNda8Zd\nhgoAAEBHWQQAAKCjLAIAANBRFgEAAOgoiwAAAHSURQAAADrKIgAAAB1lEQAAgI6yCAAAQEdZBAAA\noKMsAgAA0FEWAQAA6CiLAAAAdJRFAAAAOsoiAAAAHWURAACATrXW9r+xql3r9qO0tLQ0LOvChQvD\nsqbZY489NjRvYWFhSE5VpbVWQ8Kuf91TO+Pz8/PDsi5fvjws6/777x+WlSRra2vDsjY3N4dlJeOe\nAzP+zLa3t4dljdofjc5KkvX19aF50+ioZnzkfG9tbQ3J2XXmzJlhWSP3u+fPnx+WlSSPP/74sKyR\nr3sjzco+fHFxcUjOrpHP1+htG+n1r3/9sKyRXy/J4RyneGcRAACAjrIIAABAR1kEAACgoywCAADQ\nURYBAADoKIsAAAB0lEUAAAA6yiIAAAAdZREAAICOsggAAEBHWQQAAKCjLAIAANBRFgEAAOgoiwAA\nAHSURQAAADrKIgAAAB1lEQAAgI6yCAAAQEdZBAAAoHPqsFa0ubk5NO/ChQvDsk6fPj0sa2NjY1hW\nkpw5c2ZY1uhtW1hYGJp33G1vb09t3qtf/ephWSsrK8OyRltbWxuaN82PdRaMfL5Gfr2srq4Oy+L4\n2NraGpo3cr+7vLw8LOv8+fPDspJkfn5+aB63zuh927QeB45+nCN7wtzc3LCsw+KdRQAAADrKIgAA\nAB1lEQAAgI6yCAAAQEdZBAAAoKMsAgAA0FEWAQAA6CiLAAAAdJRFAAAAOsoiAAAAHWURAACAjrII\nAABAR1kEAACgoywCAADQURYBAADoKIsAAAB0lEUAAAA6yiIAAACdU4e1ovX19cNa1XU7d+7csKz5\n+flhWaPNzc0d9SbMtM3NzaF5TzzxxLCsjY2NYVmjjfyaGf0ccGuNfL4WFhaGZU3zfpxbZ/TzPnIm\nFxcXh2WNtry8PCxrdXV1WFbiuOdqI2dymm1tbQ3NG/n/dhxn0juLAAAAdJRFAAAAOsoiAAAAHWUR\nAACAjrIIAABAR1kEAACgoywCAADQURYBAADoKIsAAAB0lEUAAAA6yiIAAAAdZREAAICOsggAAEBH\nWQQAAKCjLAIAANBRFgEAAOgoiwAAAHSURQAAADrKIgAAAJ1Th7WixcXFoXkPPvjgsKyFhYVhWZub\nm8OyRpufnz/qTZhpo5/7s2fPDs2bViP/37a3t4dlceuNfO63traGZVXVsKxk7NfyysrKsKwkWVpa\nGpp3nI1+jVxdXR2WNXK+Rx7zJGOPx0YfKy4vLw/N43g4d+7c0LwzZ84Myxp9nDI3Nzc075l4ZxEA\nAICOsggAAEBHWQQAAKCjLAIAANBRFgEAAOgoiwAAAHSURQAAADrKIgAAAB1lEQAAgI6yCAAAQEdZ\nBAAAoKMsAgAA0FEWAQAA6CiLAAAAdJRFAAAAOsoiAAAAHWURAACAjrIIAABA59RRb8A02N7eHpa1\ntrY2LCtJTp8+PSxrcXFxWBa33tzc3FFvwjPa2toamre+vj4sa3V1dVgWt97CwsKwrEuXLg3LeuSR\nR4ZlJWNfF86dOzcsK0mWlpaG5vF1I5/3kfvJ0TY2NoZljX594WSan58fmvfAAw8Myxq9zx359bcf\n7ywCAADQURYBAADoKIsAAAB0lEUAAAA6yiIAAAAdZREAAICOsggAAEBHWQQAAKCjLAIAANBRFgEA\nAOgoiwAAAHSURQAAADrKIgAAAB1lEQAAgI6yCAAAQEdZBAAAoKMsAgAA0FEWAQAA6Jw6rBUtLi4O\nzTt9+vSwrHvuuWdY1mj33XffUW8CBzR6xldWVoZljdy2ra2tYVlJMjc3NyxraWlpWBa33sLCwrCs\n9fX1YVmjjZxxjo+R+/CNjY1hWaP3kxcvXhyWZR9+vIycy2k+Ttnc3ByWNfJ177B4ZxEAAICOsggA\nAEBHWQQAAKCjLAIAANBRFgEAAOgoiwAAAHSURQAAADrKIgAAAB1lEQAAgI6yCAAAQEdZBAAAoKMs\nAgAA0FEWAQAA6CiLAAAAdJRFAAAAOsoiAAAAHWURAACAjrIIAABAR1kEAACgU621/W+sate6/Sht\nbGwMy1peXh6WNdrIxzk/Pz8sa6SqSmutjmjdUzvjI+fyoYceGpZ1+vTpYVlJsra2NixrcXFxWNZI\nZvzW29zcHJa1tLQ0LCtJLl++PCzr/vvvH5aVJKurq0NyjmrGp3m+R87kyNeDS5cuDctKkrvvvntY\n1sjXgySZm5sbkmMf/sxGvuZevHhxWNbo45SRx88rKyvDspJxz8G1Ztw7iwAAAHSURQAAADrKIgAA\nAB1lEQAAgI6yCAAAQEdZBAAAoKMsAgAA0FEWAQAA6CiLAAAAdJRFAAAAOsoiAAAAHWURAACAjrII\nAABAR1kEAACgoywCAADQURYBAADoKIsAAAB0lEUAAAA61Vrb/8aq/W+EgVprdRTrNeMcFjPOrDuK\nGTffHBb7cGbdfjN+zbIIAADAyeQyVAAAADrKIgAAAB1lEQAAgI6yCAAAQEdZBAAAoPN/pMqvpQK1\nPHwAAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0x109214750>"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "example of hand-written digit classification with Support Vector Machines (SVM)"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "We are importing the svm.SVC (Support Vector **Classifier** class) from scikit-learn"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn.svm import SVC"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 4,
     "metadata": {},
     "source": [
      "instanciation"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc = SVC()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 10
    },
    {
     "cell_type": "heading",
     "level": 4,
     "metadata": {},
     "source": [
      "fitting"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc.fit(X, y)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 11,
       "text": [
        "SVC(C=1.0, cache_size=200, class_weight=None, coef0=0.0, degree=3, gamma=0.0,\n",
        "  kernel='rbf', max_iter=-1, probability=False, random_state=None,\n",
        "  shrinking=True, tol=0.001, verbose=False)"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "heading",
     "level": 4,
     "metadata": {},
     "source": [
      "scoring"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc.score(X,y)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 12,
       "text": [
        "1.0"
       ]
      }
     ],
     "prompt_number": 12
    },
    {
     "cell_type": "heading",
     "level": 4,
     "metadata": {},
     "source": [
      "prediction"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "y_hat = svc.predict(X)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 13
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "np.alltrue(y_hat == y)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 14,
       "text": [
        "True"
       ]
      }
     ],
     "prompt_number": 14
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "#### Have we got a perfect model ???"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Here we are making an important methodological mistake: we are using **all the instances available** to **train** the model, and using the **same** instances to **evaluate** the model in terms of accuracy. It tell us (almost) nothing about the actual performance in production of the model, just how well it can reproduce the data it's been exposed too ...\n",
      "\n",
      "A way to work around that is to train the model over a **subset** of the available instances (the **training set**), calculate the train score, and test the model (i.e. calculate the test score) over the remaining of the instances (the **test set**).\n",
      "\n",
      "**Cross-validation** consists into repeating this operation several times using successive splits of the original dataset into training and test sets, and calculating a summary statistic of the train and test scores over the iterations (usually average).\n",
      "\n",
      "\n",
      "Several splits can be used: \n",
      "\n",
      "+ **Random split**: a given percentage of the data is selected at random (with replacement) \t\n",
      "+ **K-folds**: the dataset is divided into K exhaustive splits, each split is used as the test set, while the K-1 splits are using as the training set\n",
      "+ **Stratified K-folds**:  for classification mainly. The folds are constructed so that the class distribution is approximately the same in each fold (e.g. the relative frequency of each class is preserved)\n",
      "+ **Leave One Out**: like K-fold with K = 1. One instance is left out, the model is built on the N-1 remaining instances, this procedure is repeated until all the instances have been used.\n"
     ]
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "cross-validation in scikit-learn"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn.cross_validation import train_test_split\n",
      "\n",
      "X_train, X_test, y_train, y_test = train_test_split(X, y, \\\n",
      "                                test_size=0.25, random_state=1)\n",
      "\n",
      "print(\"train data shape: %r, train target shape: %r\"\n",
      "      % (X_train.shape, y_train.shape))\n",
      "print(\"test data shape: %r, test target shape: %r\"\n",
      "      % (X_test.shape, y_test.shape))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "train data shape: (1347, 64), train target shape: (1347,)\n",
        "test data shape: (450, 64), test target shape: (450,)\n"
       ]
      }
     ],
     "prompt_number": 15
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc = SVC().fit(X_train, y_train)\n",
      "train_score = svc.score(X_train, y_train) \n",
      "train_score"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 16,
       "text": [
        "1.0"
       ]
      }
     ],
     "prompt_number": 16
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "test_score = svc.score(X_test, y_test)\n",
      "test_score"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 17,
       "text": [
        "0.40888888888888891"
       ]
      }
     ],
     "prompt_number": 17
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Ok that seems more like a 'normal' result ...\n",
      "\n",
      "- if the **test** data score is **not as good as** the **train** score the model is **overfitting**\n",
      "\n",
      "- if the **train score is not close to 100%** accuracy the model is **underfitting**\n",
      "\n",
      "Ideally **we want to neither overfit nor underfit**: `test_score ~= train_score ~= 1.0`. \n",
      "\n",
      "When setting up a Support Vector Machine classifier, one needs to set up 2 parameters (hyper-parameters) which are NOT tuned at the fitting stage (they are NOT learned). These are **C** and **$\\gamma$** (see the [relevant section](http://en.wikipedia.org/wiki/Support_vector_machine#Parameter_selection) in the wikipedia article). What we did before is to instanciate the SVC class without specifying these parameters, which means that the default are used. Let's try something else. "
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc_2 = SVC(C=100, gamma=0.001).fit(X_train, y_train)\n",
      "svc_2"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 18,
       "text": [
        "SVC(C=100, cache_size=200, class_weight=None, coef0=0.0, degree=3,\n",
        "  gamma=0.001, kernel='rbf', max_iter=-1, probability=False,\n",
        "  random_state=None, shrinking=True, tol=0.001, verbose=False)"
       ]
      }
     ],
     "prompt_number": 18
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc_2.score(X_train, y_train)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 19,
       "text": [
        "1.0"
       ]
      }
     ],
     "prompt_number": 19
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc_2.score(X_test, y_test)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 20,
       "text": [
        "0.99111111111111116"
       ]
      }
     ],
     "prompt_number": 20
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "sum(svc_2.predict(X_test) == y_test) / float(len(y_test))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 21,
       "text": [
        "0.99111111111111116"
       ]
      }
     ],
     "prompt_number": 21
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Could be luck (we only used one train / test split here): Now we're going to use **cross validation** to  repeat the train / test split several times to as to get a more accurate estimate of the real test score by averaging the values found of the individual runs\n",
      "\n",
      "scikit-learn provides a very convenient interface to do that: ```sklearn.cross_validation```"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn import cross_validation"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 22
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "cross_validation."
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "ename": "SyntaxError",
       "evalue": "invalid syntax (<ipython-input-23-3af9688d9f43>, line 1)",
       "output_type": "pyerr",
       "traceback": [
        "\u001b[0;36m  File \u001b[0;32m\"<ipython-input-23-3af9688d9f43>\"\u001b[0;36m, line \u001b[0;32m1\u001b[0m\n\u001b[0;31m    cross_validation.\u001b[0m\n\u001b[0m                     ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
       ]
      }
     ],
     "prompt_number": 23
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "cross_validation.ShuffleSplit?"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 24
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "cv = cross_validation.ShuffleSplit(len(X), n_iter=3, test_size=0.2,\n",
      "    random_state=0)\n",
      "\n",
      "for cv_index, (train, test) in enumerate(cv):\n",
      "    print(\"# Cross Validation Iteration #%d\" % cv_index)\n",
      "    print(\"train indices: {0}...\".format(train[:10]))\n",
      "    print(\"test indices: {0}...\".format(test[:10]))\n",
      "    \n",
      "    svc = SVC(C=100, gamma=0.001).fit(X[train], y[train])\n",
      "    print(\"train score: {0:.3f}, test score: {1:.3f}\\n\".format(\n",
      "        svc.score(X[train], y[train]), svc.score(X[test], y[test])))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "# Cross Validation Iteration #0\n",
        "train indices: [1109  940  192  260 1148  966 1720  554  308  512]...\n",
        "test indices: [1081 1707  927  713  262  182  303  895  933 1266]...\n",
        "train score: 1.000, test score: 0.992\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "# Cross Validation Iteration #1\n",
        "train indices: [1642  586  142   15  701  472  380 1405 1551  450]...\n",
        "test indices: [1014  755 1633  117  181  501  948 1076   45  659]...\n",
        "train score: 1.000, test score: 0.997\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "# Cross Validation Iteration #2\n",
        "train indices: [ 396 1083 1184 1569  560 1502 1722 1162 1316 1685]...\n",
        "test indices: [ 795  697  655  573  412  743  635  851 1466 1383]...\n",
        "train score: 1.000, test score: 0.989\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n"
       ]
      }
     ],
     "prompt_number": 26
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "There's a wrapper for estimating cross validated scores directly, you just have to pass the cross validation method instanciated before"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn.cross_validation import cross_val_score\n",
      "\n",
      "svc = SVC(C=100, gamma=0.001)\n",
      "\n",
      "cv = cross_validation.ShuffleSplit(len(X), n_iter=10, test_size=0.2,\n",
      "    random_state=0)\n",
      "\n",
      "test_scores = cross_val_score(svc, X, y, cv=cv, n_jobs=4) # n_jobs = 4 if you have a quad-core machine ...\n",
      "test_scores"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 27,
       "text": [
        "array([ 0.99166667,  0.99722222,  0.98888889,  0.99722222,  0.98611111,\n",
        "        0.99722222,  0.98611111,  0.99444444,  0.98611111,  0.99166667])"
       ]
      }
     ],
     "prompt_number": 27
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Cross validation can be used to estimate the best hyperparameters for a model\n",
      "\n",
      "Let's see what happens when we fix C but vary $\\gamma$"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "n_iter = 5 # the number of iterations should be more than that ... \n",
      "\n",
      "gammas = np.logspace(-7, -1, 10) # should be more fine grained ... \n",
      "\n",
      "cv = cross_validation.ShuffleSplit(len(X), n_iter=n_iter, test_size=0.2)\n",
      "\n",
      "train_scores = np.zeros((len(gammas), n_iter))\n",
      "test_scores = np.zeros((len(gammas), n_iter))\n",
      "\n",
      "for i, gamma in enumerate(gammas):\n",
      "    for j, (train, test) in enumerate(cv):\n",
      "        C = 1\n",
      "        clf = SVC(C=C, gamma=gamma).fit(X[train], y[train])\n",
      "        train_scores[i, j] = clf.score(X[train], y[train])\n",
      "        test_scores[i, j] = clf.score(X[test], y[test])"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 28
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "f, ax = plt.subplots(figsize=(12,8))\n",
      "#for i in range(n_iter):\n",
      "#    ax.semilogx(gammas, train_scores[:, i], alpha=0.2, lw=2, c='b')\n",
      "#    ax.semilogx(gammas, test_scores[:, i], alpha=0.2, lw=2, c='g')\n",
      "ax.semilogx(gammas, test_scores.mean(1), lw=4, c='g', label='test score')\n",
      "ax.semilogx(gammas, train_scores.mean(1), lw=4, c='b', label='train score')\n",
      "\n",
      "ax.grid()\n",
      "ax.fill_between(gammas, train_scores.min(1), train_scores.max(1), color = 'b', alpha=0.2)\n",
      "ax.fill_between(gammas, test_scores.min(1), test_scores.max(1), color = 'g', alpha=0.2)\n",
      "\n",
      "ax.set_ylabel(\"score for SVC(C=%4.2f, $\\gamma=\\gamma$)\" % ( C ),fontsize=16)\n",
      "ax.set_xlabel(r\"$\\gamma$\",fontsize=16)\n",
      "best_gamma = gammas[np.argmax(test_scores.mean(1))]\n",
      "best_score = test_scores.mean(1).max()\n",
      "ax.text(best_gamma, best_score+0.05, \"$\\gamma$ = %6.4f | score=%6.4f\" % (best_gamma, best_score),\\\n",
      "        fontsize=15, bbox=dict(facecolor='w',alpha=0.5))\n",
      "[x.set_fontsize(16) for x in ax.xaxis.get_ticklabels()]\n",
      "[x.set_fontsize(16) for x in ax.yaxis.get_ticklabels()]\n",
      "ax.legend(fontsize=16,  loc=0)\n",
      "ax.set_ylim(0, 1.1)\n",
      "ax.grid()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAu8AAAH5CAYAAAAx5NBKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl8U1XeP/DPyd4k3Re6pNBF2VFQtoJAC6iIouO4O/KA\n+zIw6jCjoyLL87jgvo6i4zwyij9HHnVYRPZ9q6BC2UGW0rK0pXTf0izn90fa0DRpm7Yp6fJ5v155\nJTn33Jtvepvkm5PvPVdIKUFERERERO2fwt8BEBERERGRd5i8ExERERF1EEzeiYiIiIg6CCbvRERE\nREQdBJN3IiIiIqIOgsk7EREREVEHofJ3AJeSEILzYhIRERFRm5NSirbYbpcbeZdS+vUye/bsdrG9\n5qzXVN+WLm9Ou6//bu1h/7WHfdeV9l9nfO011ae5y9rrvuuI+4+vvbaLozO+9trz/utorz1v+l6K\n115bUs6ZM6dNH6A9mTt37pz28HwTEhLaxfaas15TfVu63Nv2jRs3IjU1tdHHuFR8uf/aw77zpk9n\n2X+d8bXXVJ/mLmuv+w7oePuPrz1XXe29syXL2uv+62ivPW/6tvVrb+7cuZgzZ87cRh+khURbfzto\nT4QQsis9385kzpw5aA9fvKhluP86Lu67jo37r2Pj/uu4hBCQLJuhrqw9jDxQy3H/dVzcdx0b91/H\nxv1HnnDknYiIiIjIhzjyTkRERERETN6JiIiIiDqKLjXPe1OEaJNfN6gLYDkWERERXQpM3uthEkbN\nxS99REREdKkweSciog6D0+YRtRxfP50Dk3ciIupQmIAQNR9fN50HD1glIiIiIuogmLwTEREREXUQ\nTN6JiIia6eDBgxg3bhwMBgPi4uIwe/Zs2O12n63r7fab6rdo0SLceOONiI2NRWBgIAYPHox///vf\nXsWZmZkJhUKBrKwsr/p3Ni3dx4sXL8YVV1wBnU6HpKQkvPPOO259vvrqKwwaNAiBgYEwmUyYMmUK\nzp0759JnwYIFUCgUbpdPP/3UbXtWqxXz5s3D5ZdfDp1Oh/j4ePz5z39u+ZOndo01753Y4sWLcfLk\nSTz99NNtsv3MzEwsWLAAU6ZMQWJiYps8BhFRe1NYWIjx48ejf//+WLp0KY4dO4YZM2bAbrfjf/7n\nf1q9rrfb96bfu+++i6SkJLz//vuIiIjA8uXLce+99yI/Px/Tpk1ruz9SB9fSfbxt2zb8/ve/x0MP\nPYS3334b6enpePbZZ6FQKPDkk08CAL7//ntMnjwZ06ZNw9tvv42zZ89i5syZuPHGG/HLL7+4zWC2\nYcMGBAQEOO97+rydOnUqNmzYgDlz5qB3797IysrCoUOHfPTXoHZHStllLo6n27Cmlnc0U6ZMkSaT\nqc22v2HDBimEkOvWrWuzx+gIOtv/DVF7Nnv2bH+HIF955RUZFhYmS0tLnW2vv/661Ov1sqSkpMXr\n1rZ5u31v+l24cMEthnvvvVcmJiY2+TxPnjwphRDy1KlTTfa91KxWq6yurm6z7bd0H1933XVy9OjR\nLm0zZsyQYWFh0mKxSCmlvOOOO+TgwYNd+ixdulQKIeThw4edbZ9//rkUQsjy8vJGY12xYoVUq9Xy\n0KFDjfZrD6+drqQmN2iTfJZlM9RqsgPMjW+z2WCz2fwdBhF1AitWrMD1118Po9HobLvrrrtQWVmJ\nTZs2tXpdb7fvTb+wsDC3GAYOHIizZ88281l7duDAAUyYMAHh4eEwGo3o27cvPvroI5c+mzdvRlpa\nGgIDAxESEoK0tDTs2bPHuXzRokUYMGAAdDodunfvjpkzZ7q9X0+dOhVDhgzB4sWL0a9fPwQEBGDn\nzp0AgC1btmDMmDEwGAyIiIjAI488grKyslY9r5bu44yMDFx77bUubddeey0KCwuxY8cOZ1tQUJBL\nn+DgYACeP0+b+oz93//9X4wbNw69e/dutB91Hkzem0nMFZfk0lpTp07FF198gTNnzjjr5JKSkpzL\nz58/j8ceewwmkwk6nQ59+vTBP/7xD5dt5OTkYMqUKYiLi4NOp0NsbCwmTZqE8+fPY+PGjRg7diwA\nxxtT7WNs3ry5wZhWrVqFESNGICQkBIGBgejdu7fbz48ZGRm49dZbERERAb1ej969e2PevHkufd55\n5x306tULWq0WsbGxmD59OkpLS136KBQKzJw5E/PmzUNiYiK0Wi32798PANi0aRPGjRuHoKAgGI1G\nTJgwAQcOHGj+H5mIOiSr1drkpTFHjhxxS5S6d+8OvV6PI0eOtHjdw4cPAwAOHz7s1fZbGseOHTvQ\nq1evRuP01qRJk6BWq/HVV19h2bJlmD59ukvivHHjRowbNw5arRZffPEFFi1ahFGjRjm/PKxevRp3\n3303Bg8ejKVLl2L69Ol488033Up6hBDIzMzEs88+ixdeeAErV65EQkICtm3bhvHjxyM2Nhbfffcd\n3n33Xfz444+4//77XdZv7j5v6d+2qqoKGo3Gpa32fm0ZyyOPPIJt27bhyy+/RElJCY4ePYqZM2c2\nmIAnJydDrVajd+/eHuvdd+7cicsvvxzTpk1DcHAwDAYDbrvtNrcaeuo8WPPeSc2aNQv5+fnYtWsX\nli1bBgDQarUAgJKSElxzzTUwm82YO3cuEhMTsXLlSjz++OMwm83ON83JkycjOzsbb775JuLj45GT\nk4P169ejoqICV199Nf7+97/jj3/8Iz744AMMGTIEANCnTx+P8Zw4cQI333wz7rzzTsyZMwcajQZH\njx7FyZMnnX127tyJ1NRU9OzZE++++y5MJhOOHj2Kffv2Ofs8//zzmDdvHqZNm4ZJkybhwIEDePHF\nF5GRkYFNmza51AouWLAAycnJePvtt2EwGBATE4Ply5fjlltuwaRJk/DVV19BSonXXnsNo0aNwt69\ne2EymXy7I4ioTW3fvh0PPfQQzp07h5deegl//OMfAQBffvklJk2ahJCQEJf+CxYswAMPPNDkdhs7\nMLGwsNBtuwAQGhqKwsLCRrfrzbpFRUVebb8lcaxbtw5LlizB559/3mic3sjPz0dmZiaWLVuGfv36\nAQDS0tJc+jz33HMYNGgQVq5c6Wy77rrrnLdnzZqFtLQ0Zzy1y5577jnMnDkTcXFxAByjzxcuXMC6\ndetwxRVXONe/++67cc011+Drr792tsXFxWHcuHE4ePAg+vbt26J93tJ9fNlll2HXrl0ubbW/EBQU\nFAAAxo8fj88++wwPPvggpkyZAgAYMWIEvv32W5f1YmNj8dJLL2Ho0KGw2Wz4+uuv8dhjj6GiogJP\nPfWUs9+5c+ewYMECDBw4EN988w1KSkrwzDPP4NZbb0V6enqTz5s6HibvnVRSUhIiIiKg0WgwdOhQ\nl2XvvfcesrKysH//fiQnJwMAxo4di6KiIsydOxdPPPEEFAoF0tPT8eqrr+Kee+5xrnv77bc7b9cm\n6n369HF7jPp+/fVXWCwWfPzxx86fIVNTU136/OUvf0FkZCTS09Oh0+nc+hQUFOCtt97C1KlT8f77\n7wNwjPpHRkZi8uTJ+OGHHzBp0iSXba5evdr5pQUAnnzySaSlpeE///mPsy0tLQ1JSUl46623PM4K\nQETt0/79+/Hyyy9j3rx5UKlUmDt3LsaMGYP+/fujtLTUY/J188034+eff/ZDtP6XmZmJe++9F7/7\n3e/wX//1X63eXlhYGOLj4/Hoo4/iT3/6E1JTUxEVFeVcXl5ejp07dzrfr+uz2WzYvXs33nvvPZf2\nO++8E88++yzS09Nx2223OdtNJpNL4l5RUYH09HR88MEHLiPnI0eOhFqtxs8//4y+ffte0n3+2GOP\n4bHHHsNnn32G2267DTt37nR+rigUjmKH5cuX4+GHH8af//xn3HDDDcjJycGcOXNw6623Yu3atc5+\n1113ncsXneuvvx5VVVV4+eWX8eSTTzoHq2rLapYsWYLQ0FAAQExMDMaMGYMNGza4faGijo/Jexe0\ncuVKDB8+HAkJCS5veNdddx0+++wzHDx4EP3798eQIUPw+uuvw263Iy0tDf3793c7Ct5bgwYNglqt\nxl133YUHHngAo0aNcnmTr6iowPbt2/HMM884E/f60tPTYbFYcN9997m033XXXbj//vuxefNml+R9\nwoQJLon7b7/9hhMnTuCFF15wed4BAQEYPnx4oyU/RNT+LFu2DEuWLIFK5fgou+qqq/DJJ5+guLgY\n48aN87hOWFiYW71xc4WGhqK4uNitvbCw0Jk8tWZdb7ffnDgKCgpwww03IDExEV999VWjMXpLoVBg\n9erVeOGFF/DAAw+gsrISI0eOxPvvv4+BAweisLAQUkrExMR4XD8/Px8WiwXdunVzaa+9XztSXb+9\nVmFhIWw2G5544gk88cQTLsuEEDh9+jSAlu3zlu7jBx54ABkZGXj88cfxyCOPwGAwYN68eZg+fTqi\no6MBAH/7299w++2349VXX3WuN3DgQPTu3RtLlizBrbfe2uD2b7vtNixatAhZWVno0aOH8/klJye7\nxDVy5EhoNBocPHiQyXsnxOS9meTs9n9wZlPy8vJw/PhxqNVqt2VCCFy4cAEA8M0332Du3Ll4/fXX\n8dRTTyEmJgaPPfYYZs6c2ewkPjk5GatWrcJrr72GyZMnw2w2Y+jQoXjttdcwevRoFBYWwm63N1q2\nUvtGXv+DQKVSITw83O2Nvn6/vLw8AMCDDz6IBx980G37tW+ERNQxPPfccy73o6OjUVlZid9++w0j\nR470uI4vymZ69+7tNg1fdnY2Kioqmjxo0Jt1vd2+t/0qKipw0003wWq14ocffmhwgKQlevXqhW+/\n/RY2mw2bN2/Gs88+ixtvvBFnzpxBaGgoFApFgwfHRkREQK1WO9+ba+Xm5gLwfLBtXSEhIRBCYO7c\nuZg4caLb8tjYWAAt2+ct3ccKhQIffPABXnrpJZw+fRqJiYk4ePAgAGD48OEAHGWkteUytXr27ImA\ngACcOHGi0Rg9ffb26dMHVVVVbu1SyhYPuFH7xuS9C4qIiEB0dLTbT5W1evbsCQCIjIzEhx9+iA8/\n/BC//fYbFixYgNmzZyMyMhKPPfZYsx83NTUVqampsFgs2Lp1K2bNmoUbb7wRp06dcr7J146UeFL7\nRn7u3DmX2nqr1YoLFy64vdHXf9MKDw8HAMybNw/jx4932379g4yIqOMpLS2FwWBocLkvSihuuOEG\nvPHGGygrK3OWAX7zzTfQ6/UYM2ZMq9f1dvve9LNarbjjjjtw/PhxbN++HREREa167g1RKpVIS0vD\n008/jT/84Q/Ouv1hw4bhiy++8DinvFKpxNVXX41Fixbh0UcfdbYvWrQICoUCKSkpLv3rv6cbDAYM\nHz4chw8fxsyZMxuMrSX7vDX7GHDMHlM7g8xHH32EkSNHOj9bExIS8Ouvv7r0P3ToECorK5GQkNDo\ndr/99ltERES4DDbddNNNmD17Ni5cuOD8nNu8eTMsFgsGDhzo9XOmDqSt5qBsjxd0sXneH3nkERke\nHu7WPmfOHBkeHi7z8vKavc2wsDA5bdo0KaWU27dvl0II+cMPP7QoviVLlkghhPz555+llFKOGTNG\nxsfHy8rKSo/9CwoKpFarlQ8//LBL+8KFC93iEELIF1980W0biYmJcvLkyS2KtyGd7f+GqD1raq7q\nhx56qMl5sVursLBQxsTEyGuvvVauXbtWfvLJJ9JoNLq95/zrX/+SSqVSZmVlNWtdb7fvTb+HH35Y\nCiHk+++/L3fs2OFyMZvNjT7PpuZ5z8jIkNdee6385z//KdevXy+/++47eeWVV8pBgwY5+2zevFlq\nNBo5YcIE+f3338uVK1fK2bNnO9+vV69eLYUQ8v7775crV66Ub7zxhtTpdPLxxx93eawpU6a4zY0u\npZRbt26VWq1WTp48WS5evFiuW7dOfv755/KOO+6QR48ebfT5Ncabv62n/Zueni7feOMNuWbNGvnd\nd9/J22+/XQYHB8t9+/Y5+/z973+XCoVCzpgxQ65Zs0YuXLhQ9uzZUyYlJcmKigpnv9tuu02++eab\ncsWKFXLZsmXyvvvuk0II+eGHH7rEWlJSIrt37y5TUlLksmXL5FdffSVNJpO87rrrXPpxnvdLC204\nz7vfE+pLeelqyft7770nhRDy448/ljt37pR79+6VUkpZXFws+/TpI3v16iXnz58v169fL5ctWybf\neOMNecstt0gppSwqKpKDBw+W7777rlyxYoVcu3atnD59uhRCyKVLl0oppczPz5dqtVreeuutcuvW\nrXLXrl0uJ7So6+OPP5b33nuvXLhwody4caP89ttv5dVXXy1NJpOsqqqSUkq5a9cuqdfr5cCBA+WX\nX34p169fLz/77DM5ffp053aef/55KYSQTz31lFy1apV89913ZWBgoNtJMRpK3n/88UepVqvlXXfd\nJb/99lu5ceNG+c0338gnn3xSvv322y36O3e2/xui9qypBGTGjBmXJI6DBw/KsWPHyoCAABkbGytn\nzZol7Xa7S58FCxZIhULhlvx6s643fRrqZ7PZnMsTEhKkQqGQQgiXi6e46msqec/Ly5OTJ0+WSUlJ\nUqfTyejoaHnvvffK7Oxsl36bNm2So0ePlnq9XoaEhMixY8fKjIwM5/JvvvlGDhgwQGo0GhkfHy9n\nzpzp8hyklHLq1KlyyJAhHuP46aef5IQJE2RQUJA0GAyyb9++csaMGbK4uLjR59eUpvaBp/37yy+/\nyCFDhkij0SiDgoLkTTfdJPfv3++27U8++UReeeWV0mg0yri4OHn33XfLkydPuvR5/vnnZa9evaRe\nr5cBAQFy8ODBcuHChR5jPXbsmJw4caI0GAwyNDRU3n///bKoqMilD5P3S4vJO5P3FikvL5f33HOP\nDA0NlUIIlzPqFRYWyqefflomJiZKjUYjo6Ki5OjRo+V7770npZTSbDbLRx99VPbr18/5JjR06FD5\n9ddfuzzGJ598IpOSkqRKpZIKhUJu2rTJYyw7duyQt9xyi4yPj5darVbGxMTIO++8021kZPfu3XLS\npEkyJCREBgQEyD59+sjXX3/dpc8777wje/XqJTUajYyNjZXTpk1z+9LQUPJeG8tNN90kQ0NDpU6n\nkwkJCfKee+6R6enp3v1h6+ls/zdE7VljCci5c+fkU089demC6eTa8xlWqfmYvF9abZm8C8f2uwYh\nhGzs+Qoh0JX+HuQb/L8hunTmzJmDOXPmeFy2cOFC6PV6/P73v7+0QXVSmZmZSEpKQmZmJrp37+7v\ncKiVGnvtkO/V5AZtcsQwz7BKRESdwrlz5zwejE5E1JlwthkiIuoU/vrXv/o7hE6HUw0StT9M3omI\niMhNQkICbDabv8MgonpYNkNERERE1EEweSciIiIi6iD8nrwLIUxCiA+EEDuEEBVCCLsQwqvD2oUQ\nOiHEG0KIczXrbhdCjGrrmImIiIiI/KE91LxfBuAOAD8D2Azgumas+08AEwH8BcAJANMArBJCpEgp\nM3wdKBER+R+nuyOirszv87yLOpOvCyEeAvApgAQpZVYT610JYDeA+6WU/6ppUwI4AOCIlPIWD+tw\nnnfyOf7fEBERUV2dep73RrPpxt0MwALgmzrbsgH4N4DrhRBqH4RHRERERNRu+D15b4V+AE5IKavq\ntR8EoIGjHIeIiIiIqNNoDzXvLRUGoNBDe0Gd5V3a4sWLcfLkSTz99NM+3/bUqVOxadMmnDx50ufb\nJiKi5rHZgNJSoKTEcSkudlxXVQFSAnb7xUtj95vTtzPdZ+UjdSQdOXmnJixevBjr1q1rk+R91qxZ\nKC0t9fl2iYgaY7aaYbaZIaV0nv1TQEAIAYGa+zW3PS1vbB1/sFpdk+66iXf9i6f22rbycr+ET0R+\n0JGT90IAnqaUrB1xL/CwzGWWgtTUVKSmpvo6rg7JbDZDq9V63T8pKakNo7m0as8gqFQq/RwJEdUy\nW82oslah0lqJEnMJiiqLUGIugVVaAcCZdAO+PWhcIRQQEFAoFC4JvkIoHMtq7ttsApWlalRWqFFR\nqkJ5mRIVZSpUlKlQVqJERbkK5aVKlJcpUV6qQFmpEmWlCpSXKlFWpkBZqeNSUd6Rq1eJ6KKNNZe2\n5/fZZupq5mwzswC8ACC4bt27EGIOgL8BCJRSWuqt0+rZZi7V4Exrd8vUqVPxxRdfuLQlJCTgxIkT\n2LhxI8aOHYvvvvsOP/74IxYvXgyr1YrCwkIcO3YMc+fOxbZt25CTk4OYmBhcf/31eOWVVxASEuKy\n/bplM5mZmUhKSsL8+fNx+vRpfPbZZ6isrMSoUaPw8ccfIy4urtF4V61ahblz5+LgwYOw2WyIi4vD\nH/7wB7z44ovOPhkZGZgzZw62bNmCiooKdO/eHVOnTsXf/vY3Z5933nkH8+fPR2ZmJsLDw3Hbbbfh\nlVdeQWBgoLOPQqHA888/D6PRiE8++QTZ2dn45ZdfcOWVV2LTpk347//+b+zatQt2ux3XXHMN3nrr\nLfTr16/B2DnbDFHLNZak1yboKoUKOpUOGqUGCtG6ZNdqgSOhLlOirKROcl2mdCTWJY7kutyZaNf0\ncfZ1XFdV8ss+ETWm7Wab6cgj70sBzAFwJ4AvAEAIoQJwF4BV9RP3rmbWrFnIz8/Hrl27sGzZMgBw\nG1mfPn06Jk6ciK+++gpVVY7vP+fOnYPJZMLbb7+N8PBwnDhxAq+88gomTpyI7du3u6zv6WfmV199\nFSNHjsTnn3+O3NxczJgxA/fddx82bNjQYKwnTpzAzTffjDvvvBNz5syBRqPB0aNHXerpd+7cidTU\nVPTs2RPvvvsuTCYTjh49in379jn7PP/885g3bx6mTZuGSZMm4cCBA3jxxReRkZGBTZs2ucS7YMEC\nJCcn4+2334bBYEBMTAyWL1+OW265BZMmTcJXX30FKSVee+01jBo1Cnv37oXJZGrGHiCiupqTpAfr\nghtM0kuLlcg5o3GMZJcoUV7mGNWuTazLSmsTbIUz0a5tKytVwlzVeUe6DUYbDEYbjEGOa0OgDRqt\nhFIpHQNPQkIIu+NaYYdEzW0hoVBIQAAKBZz3hQCUSgVUSgXUSgWUCiXUKiXUSse1UqGASiGcfZQK\nBVRKx33Hdhzbq7348n5bbJvIlwYObLttt4vkXQhxe83Nq2uuJwoh8gHkSSk3CyF6ADgOYK6U8n8A\nQEq5RwjxDYB3a6aFzATwOIAeAO65pE+gHUpKSkJERAQ0Gg2GDh3qsc+wYcPw6aefurSNGjUKo0Zd\nPEltSkoKkpOTMXr0aOzZswcD6/w3ehptTkxMxMKFC533z58/j7/+9a/IyclBdHS0xzh+/fVXWCwW\nfPzxxzAajQDgVs70l7/8BZGRkUhPT4dOp3PrU1BQgLfeegtTp07F+++/DwC49tprERkZicmTJ+OH\nH37ApEmTXLa5evVqly80Tz75JNLS0vCf//zH2ZaWloakpCS89dZbeOeddzzGT0QXtTZJN1cJZGVq\ncCZLi7NZWpzNrrmdrcGZU1qUlrSLjy2fEUJCb7TDWJNsGwMd147k2+5MwmvbjQ206w12KNrge4nN\nboNd2mGTNljtVtjsNtikow0AJKRLCRMAKIUSWpUWWqUWGpXG+auJTqWDSqHyeGntLypEXUl7eRdc\nVOe2BPBRze2NAMYCEHBMa1n/u/H9AF4G8BKAEAB7AEyQUu5py2A7i1tvvdWtrbq6Gm+++Sa++OIL\nZGVlOUfkAeDo0aMuybsnEydOdLnfv39/AEBWVlaDyfugQYOgVqtx11134YEHHsCoUaMQFRXlXF5R\nUYHt27fjmWeecSbu9aWnp8NiseC+++5zab/rrrtw//33Y/PmzS7J+4QJE1wS999++w0nTpzACy+8\nAKvV6mwPCAjA8OHDsXnz5kafN1FX09Ik3WoF8s5pcCxLgzPZrgn6mSwNLuRp/Pm0vCfsgLYE0JZC\n1F7rii/e111sh64YQlsKaEsgdLXXJY711eWwKSSKARQ39ZhVNZd878PUKDUI1gYjRBeCYG0wgnUX\nb4foQhCsC3berm1XKy+eJkWpUEIJJdTw/tQpdmmHzW6DxW5BVVUVCmWhI+mHzfEJX6Nu4i+EgFap\ndSb93iT8SgVLl6hrahfJu5Sy0a/cUspMeJiTvqbWfUbN5ZLoTKXNMTExbm3PPfccPvzwQ8yePRsj\nRoxAYGAgsrOz8fvf/94lkW9IWJjrDJ21CXJj6yYnJ2PVqlV47bXXMHnyZJjNZgwdOhSvvfYaRo8e\njcLCQtjt9kbLVgoKCjw+J5VKhfDwcOfyWvX75eXlAQAefPBBPPjgg27b79GjR4OPTdSZNTdJF1Dg\nwnlVvVFzLU6f0uBsthZ55zSwWf1YoyBsNUl3/Uux53adp/ZiQF0BKBwfCLUfCy36eLDXXNpIlbUK\nJeYSZJdke72OQW1wSe69Sfx1qosDKwqhgEKpaFbCL6V0ju6XWktRVFUEu7TDarc22F8IAY1SA61K\nC51S55L4q5VqKBVKj0k/UUfH/+IuzFPN+r///W9MmTIFzz//vLOtpKSkzWOpnfnHYrFg69atmDVr\nFm688UacOnUKoaGhUCgUOH36dIPr135pOHfuHPr06eNst1qtuHDhgtuXivrPPTw8HAAwb948jB8/\n3m37Gk0HGQ0kaqGmknQhBJRCCZ1KB4U5DBdO6xwJepYGp7O0OHPKkaznnNG2bV250gyEZAL6800k\n2A0k5Jpy999wyUW5pRzllnKcLTvr9TpapfZisq+9mPS7JP5a1y8DerX+4tSdQkAlmpdc1yb8NrsN\n5ZZylJhLnF8AGqNRapxJf+1of7A2GN2M3bx+bCJ/YvLeiWm1WlRWVjZrncrKSqhUrv8Wn3/+uce+\nbTEvslqtRlpaGv7617/id7/7HU6ePImrr74a11xzDRYuXIhZs2Z5LJ1JSUmBRqPBv//9b4wdO9bZ\n/s0338BqtTY5JWjv3r2RkJCA/fv345lnnvH10yJqN7xJ0m3VahSeC8b5M91xLluH7FMaZ3J+LluL\nstI2/OgQNiDoNBByEgg96X5tPOcc8ab2w2wzI7c8F7nluV6vo1KoGhzJr1/KU7vcqDE6S7BakvAD\ncNbtV1kCTtOcAAAgAElEQVSqUF5dDqvdipP2kxibONalZIiovWLy3on169cP//jHPzB//nxcffXV\n0Ol0GDBgQKPrTJgwAf/6178wYMAAJCcn4/vvv8eOHTs89vXV9Ijz58/Hli1bMHHiRJhMJuTn5+PV\nV19FXFycs2b+zTffxJgxY5CSkoIZM2YgLi4OJ06cQEZGBt5//32EhoZixowZePXVV2EwGHDDDTfg\n0KFDePHFFzFq1CjceOONTcbx97//Hbfccguqq6txxx13ICIiArm5udi+fTt69OjRJie7ImorjSXp\ndqsC+bk6nD9tRN6ZUORkmWrKXLTIOa1BQX4b/9JkyG04OQ/KBlQWGNQGxATGIMYYg9jAWMQYhyAm\n8GbEGmMRbYx2KdPobGrfWyUkpJSQdQpyau/X7VO3vf5ti93iGEmvLkdRVRGKqopQbC52XFcVX7xd\nc11iLnEejNrWrHYr8ivykV/hfRG/UigRpA3yOJJfv6yn9naQNsilPr62jh91SubzK/JRWl2KsIAu\nf3J26gCYvHdiDz30ENLT0/H888+jqKjIOc870PCo+QcffAApJV544QUAwI033oivv/7abcYaIYTX\nI+9N9Rs4cCBWrlyJ5557Dnl5eQgLC8OoUaPw9ddfO2vmBw8ejG3btmHWrFmYPn06zGYzEhIScP/9\n9zu38/LLLyMyMhLz58/HRx99hIiICEyZMgWvvvqqV3HecMMN2Lx5M15++WU8/PDDqKysRHR0NFJS\nUnDPPV1+AiNqpzwl6cVVJTh/XoHc03qcPaVHTrYBead7IPd0AHJOa5Gfq4bN1oalLdrihpPzkExA\nUwGD2oC4wDhngh4TeAVijRMQE+hI1gM1gR7fO2x2R1lEU0lr3fsNtXtapz0QQkABxwmhakeZnSeJ\nqjmBVG1b7Qmk6q5T97bFZkFpdSmqrFUwBTmOG5KQUAolNEoN1Ao1NEqN829tl3aUVZc5k3uXZL/+\ndZ3lTZWq+IpN2lBYVYjCqkKv1xEQCNQGOpN9T/X8epUepiATk3fqENrVSZrami9O0kRUH/9v6FKo\nn6Rn55TiyDErsk9qcC47AOey9cg747jk52hRbW7DmThq687rJ+ahJxy3AwoRpAtCtDEascZYZ0Lu\nSNJjEGuMRaA20G2zVrsVVrsVFpvFMS2hdJz9uO6sJGqlGjqVDkqhdElaIVyTWW8S3dqzqNY9c2rt\nNYAm2wC4LPfU1tzttBW7tMNsNcNsM8NsNaPcUo5ScynKqstQYalwGW0XQjiTerVC3eSsLlJKVFgq\n3BJ8t+S/qhhF5otJv9lmbrPn2xIxxhgcnnYYQdogf4dCnUBNbtAmL2om767LmYRRs/H/hnypNkkv\nLKvEvkOVOHikGsdPSJzN1iI3OwB5Z/U4fzYAleVtWJvrRd15cEDQxVHz2tKWmsQ8JjAGRo3Rubna\nWUNqk3KL3fM59HQqHQJUATBoDAhQB0Cv1l9MIpVqrxJJaplqW7Uzua+yVKG0utSR3FvKYLFd3F8C\nAkqF0nnQZ2tmb6myVjVcxlPl+UtAuaXcF0+3Qa+Nfw3PjORxT9R6TN59hMk7tQX+31BrVFnM2LSj\nDN8vrcbeDCD3rA755wJQWtjGNd1N1J2HGA0uI+V1r6ON0c7kvHa0vO6lLiklVAqVMxkPUAfAoDZA\nq9K6JeZtOfJMLWe1W11G7UurS1FmLnOW49SlFErn/qxbjuMrFpvFmeA3mvibL7aVmL2fMW3CZROw\n4g8rfBozdU1M3n2EyTu1Bf7fUHPYpR0nT5fj28VVWLHSjl93BLdNoq4p8ZyYh54AQjIRGqR1Scgv\nHhjqSNgDVAHOUfLaUfOaCnIICEhIKIQCaoUaBvXFkXK9Wg+1Uu1ST83R8s6rLctxfMVqt6LUXOqS\n9Ncm/mdLz+K7Q985+wZrg1HwbAHP+EqtxuTdR5i8U1vg/w01pazSjFXrK7BkmQWbN2px6mgQ0Nr3\n9Ibqzmuuw8IEYgKjnQl53bKWSH0kVErVxdFyaXWc9KbOBOgKoUCAqs5oucbgON19zUh5bQLG0XJq\njD/KcZrDZrdh7BdjXcpx9j++H/2i+l2Sx6fOqy2Td842Q0TkY1JK7D1chu+XVGPlKom9u4JRVR7a\nvI00UXceFlmNuOCLZS3RxmjEGEchynAbIvQRUCqUsNqtHqf9s8EGrUILo8YIg9oAvUbvTMZrE3Oe\niZJ8oTYZD4TnA5SbW47j6y+NSoUSA6MHYlv2Nmfb+pPrmbxTu8Z3ZyIiHygsqcaS5RVYttyK7Vt0\nyMlyT1YapC4DEtcDyWuAiENA6ElERFchNiTCOVNLN0M3RBmuQpR+AsL0YVAqlC5zfdeOmtetKw9Q\nBSBAHeBSg6xWqlkSQO2CSqGCSqOCAQYAQAxinMs8leOUVZeh1FyKwqpCn5bjDIoe5JK8bzy1EdOH\nTW/lsyNqO0zeiYhawGaT2P5zORYvrcaaNQKHM4JgqQ7xfgPd9gCXrQQuW4WQyw9iVNIwDI0dClPQ\n1QjXh0OjdD1RUu0UibV15XqV3nHQZ53RSJ4dkjoLhVAgQO348umJL8txBkYPdLm/PWu7o4yMJWHU\nTjF5r4cvViJqyNkcC75fVokffrRi5/YAFOYZm16pVkC+Y2T9spVA8mrExSoxqvsoDImbhJ7hM2BQ\nGxChj4BBY+AUiURN8GU5Tvfg7tAoNai2VQMAcspzcKr4FBJCEi7FUyFqNibvdfCgQyKqq7paYt3m\nSvxnaTU2blDg+CEj7DYvT+AibIApvSZZXwXE/oI+Ub0wwjQCQ+LmoHtwdwRpgxAXGIdwfTiMGiMH\nD4h8oDnlOBWWCpRWl6JnWE/sP7/f2W/tibV46KqHLnnsRN5g8k5EVENK4NgJC75dUoWVK+34dWcA\nyor1APTebSAo+2KynrQOSn0pro65GsPihmNI3HREGaIQoY9AXFAcgrXBDZYEEFHbaKgcZ0jsEJfk\nfVPmJibv1G4xeSeiLq2sDPhxTSWWLrNiy2YFso4bAHhZO66sAhI2OZL1y1YCkYcQoA7AcNNwDI17\nFFfHXI3QgFDnAafBumDO4kLUDqUmpOLzjM+d9+sewErU3vBThIi6FJsNyNhnxXeLzViz1o69vwbA\nXNmMEfCIQxeT9R6bAU0lwgLCMMI0AkPj7sGV0VciRBsCU7AJ4QHhCNIGsRyGqJ0bkzAGCqFwzmJz\nsugk8sryEGWM8nNkRO6YvBNRp3f+PLB0eRV++NGK7dtUyDurg9dvf9piIHGdc2YYhGQBAOKD4jEi\n/mYMiR2C3hG9EWWIQmxgLEIDQqFXe1lmQ0TtQkxgDJJDk/FbwW/OtrUn1+LeAff6MSoiz5i8E1Gn\nYzYD6T/Z8N0SM9ZvkDiyPwBWi877DcT8fDFZN6UDSisAoE9EHww33Y9hccOQGJKImMAYRBujEawL\ndpvakYg6Do1Sg4HdBrok7xtObmDyTu0Sk3ci6vDsdiA7G/jPUjOWr7Di551qFF3QwOsDTQ25QPJq\nR8KetAYwngfgmEbuqpirMCxuGIaahiI+KB5xgXGINEQiSBvEkx0RdSLXdL8G/3fo/5z3WfdO7RWT\ndyLqkEpKgE2bbPh+aTW2bgWOH9FBSi0AbdMrKyxA/Paa0fWVQLcMQOGYKlav1mNYXBqGxg3F4JjB\niAuKgynIhNCAUBg1zZjXnYg6lHFJ41zuH7lwBMVVxQjWBfspIiLPmLwTUYdQXQ0cPQp8v7gaq9fa\nsPsXNSrKVAC8PNg05OTFaRwT1wO6UueisIAwpJhSMCxuGAZFD4IpyISYwBiE6EKgVXnxZYCIOrzE\n0ETEB8UjuyQbgGNO+PWZ63Fr71v9HBmRKybvRNQu2e1Afj6weo0dS5dXI32HQHamFoCXteWqCiBx\nw8WZYcJ/A+pM+mIKMiHFlILhpuG4stuVMAWZEGmIRLA2mGczJeqCAlQBGBA1wJm8A8DGkxuZvFO7\nw+SdiNqNigogIwP47j8WbNxkx769KlRXKQF4ebBp1L6LyXr3rYDa7LK4b0RfDI0bipT4FPSN6AtT\nkIlnNyUiAIAQAsNNw/HjsR+dbVuzt/oxIiLPhJTS3zFcMkII2ZWeL1F7V10NnD0LLF9ux4+rLfh5\npwJ5OV6eIAkAdIWOA0wvW+k44DT4jMtilUKFQdGDMCxuGFJMKegV0YtnNyWiBqVnpyPlf1Oc99UK\nNYqeLYJew+lfqXmEEJBStsmoEEfeieiSsdsdB5pu3w4s/cGKrdvsOHxQBZtVAa8ONIUdiNt5cRrH\n2F2A0ubSQ6/WY0jsEAw3DUeKKQWXh1/Os5sSkVd6R/RGhD4C+RX5AACL3YItWVtw/WXX+zkyoov4\nSUZEbaq6Gjh2DFi8xI6166zY/asCRYUqeP32Yzx7MVlPWgvoC9y6hAWEOZP14XHDkRSWxLObElGz\nBWoDMSBqADZkbnC2rT+5nsk7tStM3omozWzdCjz+hA0H9isgpQJeHWyqNDvq1Wtnhum2z+VA01rd\ng7tjeNxwpMSnIMWU4pzOkWc3JaKWUiqUGBw72CV535rFundqX1jzTkRtoqQESEwECtwHyt2F/XYx\nWU/YCGjLPXZznOF0OEbGj8TQuKE8uykR+dyK31Zg4v+b6LxvUBtw4ZkLnDaWmoU170TU4axd20ji\nrikDEtfVzAyzCgg74bGbWqHGwOiBGG4ajjE9xuCKblfw7KZE1KYGxw5GoCYQpdWOc0GUW8rx0+mf\nMDphtJ8jI3Jg8k5EbWLNmnoNoceBvv/nGGGP3w6oLB7XM6gNGBI7BCmmFIxNHIteEb14dlMiumSC\ntEHoF9UP6afTnW3rMtcxead2g8k7EbWJ1aslXIrVJ04DLl/psW94QLizHGZ84nh0D+nOs5sSkV9o\nVVoM7DbQJXnfcmqLHyMicsXknYh87vhx4MSJOom70gz02OzSp3twd6SYUjC6x2ikJaShm7Ebz25K\nRO3CmIQxmP/LfOf9PTl7YLaaOaBA7QKTdyLyObeSmfhtEJpK9Inoi5T4FIxLHIcUUwrPbkpE7dI1\n8ddAq9TCbHOcpbmwqhB7c/diSNwQP0dGxOSdiNqAW/KevBp/u+Zv+NOwP/HspkTU7oUGhKJPZB/s\nydnjbFt7Yi2Td2oXOF0DEfmU1QqsXVtvStbkNXj4qocRbYxm4k5E7Z5erceAqAEubVuyWPdO7QOT\ndyLyqZ9/BkpK6pTBBOSjR+9CJIYm+i8oIqJmEEJghGmES9uenD2oslb5KSKii5i8E5FPuZXMJK3F\nuKQ0v8RCRNRSYxPHQikuHkB/ruwcjl045seIiByYvBORT61aVa8heQ3GJo71SyxERC0VZYzC5eGX\nu7StPrHaT9EQXcTknYh8pqQE+OmnevXuSWuQlsiRdyLqWIwaI/pH9ndp43zv1B4weScin9m4EbBa\n69S7hx9G78sMiA2M9VtMREQtoVKo3GaX2Zu3F5WWSj9FROTA5J2IfGZ1/V+Uk1djbAJLZoioY0pL\ncP3V8GThSWQXZ/spGiIHJu9E5DOsdyeiziQxNBEJwQnO+xISa0+s9V9ARGDyTkQ+kpUFHKs7EYPC\nAiRsRGpCqr9CIiJqFaPGiH5R/VzaNmdthpSygTWI2h6TdyLyCbcpIk07MLDHZQjXh/slHiKi1tKp\ndBgUPcilbV/uPlRYKvwUERGTdyLykR9/rNeQvIb17kTU4Y3pMcbl/tGCo8grz/NTNERM3onIB2w2\nYP36eo3Jq1nvTkQdXv9u/RFtjHbet9qtWHdinR8joq6OyTsRtdru3UBRUZ0GXSEUcbsxqscov8VE\nROQLnuZ735q9lXXv5DdM3omo1VasqNeQuA5D469GkDbIL/EQEfmKXq1H/yjX5H1v7l6UVZf5KSLq\n6pi8E1GrrVxZr4H17kTUSSiEAiO7j3RpO5x/GAWVBX6KiLo6Ju9E1Crl5cDOnfV+Pma9OxF1IkNi\nhyBEF+K8X2mtxJZTW/wYEXVlTN6JqFU2bACsVnGxIfQYNBFnMSJ+hP+CIiLyoSBtEPpFus73vjV7\nK+zS7qeIqCtj8k5ErbJ8eb2G5DUYET8CAeoAv8RDRORrRo3Rre59X+4+lJpL/RQRdWVM3omoVdbW\nP1N48mrWuxNRp6JWqjE0dqhL24HzB1BUVdTAGkRth8k7EbXY6dPAsWN1GoQNSNjAenci6nSGm4ZD\nr9Y77xebi5F+Ot2PEVFXxeSdiFrshx/qNcT9BEOQFUPihvglHiKithJpiESfiD4ubduzt8Nqt/op\nIuqqmLwTUYt5miJyVI9R0Cg1fomHiKiteKp735+3n3XvdMkxeSeiFrHbgc2b6zWy3p2IOqkAdQAG\nRg90aTtw/gAKqwr9FBF1VUzeiahFdu0CCut+ZmlKgLidrHcnok5rhGkE1Aq1835ueS725OzxY0TU\nFTF5J6IWcat3T9yAEIPRbWSKiKizMAWZ0DO8p0tb+ul0WGwWP0VEXRGTdyJqkTVr6jUkr0ZqQiqU\nCqVf4iEiamtGrXvd+4G8AyitZt07XTpM3omo2UpKgF9/la6NSWtY705EnZpRY8SAqAEubQfOH8CF\nigt+ioi6IibvRNRsq1cDFou42BCcCYT/xnp3IurUFEKBFFMKFOJi+nSq+BSOXjjqx6ioq2HyTkTN\ntnx5vYbkNYg0RKJvZF+/xENEdKkkhiYiKTTJpW3nmZ0wW81+ioi6GibvRNQsDU0ROS5pHIQQHtch\nIuosgnXB6BfZz6Vtf95+1r3TJcPknYia5dgx4MSJui12IHE9692JqEswaoxuyfuB8weQX5Hvp4io\nq2HyTkTNsmxZvYbYXwB9AevdiahL0Cg1GBI7xKXtWMExnCw86aeIqKth8k5EzeJpisj4oHi3GlAi\nos6qV0QvmIJMzvs2acOenD2otFT6MSrqKpqVvAshxgsh3hBCbBRCHBZCHBJCbKhpu7YlAQgh4oUQ\n3wohioQQxUKI74QQ8V6umyCE+EIIkSWEqBBCHBFC/I8QQt+SWIiocWVlwK5d7lNEst6diLqScH24\nx9IZ1r3TpdBk8i6EUAghHhVCHAOwGsB/AbACyACwF4AdwBQAq4QQx4UQjwshvPpSUJNkrwfQs2a7\nkwFcDmBDUwm4EMIIYB2AkQBeAHADgM8AzADwv948PhE1T3o6UFBQJ0lXlwHxO1jvTkRdilFjRL8o\n9+T9fPl5P0VEXYnKiz4ZAIIBfA7g/0kpj3jqJIToBeBeAM8BeALAAE/96nkYQCKAnlLKEzXb2Qvg\nNwCPAninkXVH1qx7vZSy9of8TUKIMAB/EULopJRVXsRARF5yq3dP2ASoqpGWmOaXeIiI/CFAFYCB\n3Qa6tB06fwinS067JfVEvubNCPk/AFwmpZzdUOIOAFLKI1LK2QAuA/Cpl49/M4AdtYl7zXYyAWwD\ncEsT69aeg724XnsxAFFzISIfMZuB7dvrNSavxuVhl7vUfhIRdXZCCPSL7IeIgAhnm9lmxsHzB1Fe\nXe7HyKgraDJ5l1K+L6WsrtsmhHhSCBHTQP9qKeUHXj5+PwD7PbQfBNDU2V7W1Kz7uhCijxDCKIQY\nC+BPAOZLKXnUCJEP5eUB+/Z5qHdPHOefgIiI/KibsZt76Uwe696p7bV0tpmPAIwQQjxSUy7TUqEA\nCj20F9Qsa5CU0gJgHAAdgAMASgCsBbAMwPRWxEREHqxcCZjNdX7QCjwDRB7iFJFE1CUZNUb0j+rv\n0nYw/yByy3L9FBF1Fd7UvLupSZy/AwAhxN1CiJlwHDy6UkqZ48P4GiSEMABYAcAI4D4AWQCGAZgF\nwAZH3T0R+YDNBqxbV68xeTUggNSEVH+ERETkV0aNEf0jXZP3A3kHkFueCyklZ+CiNtOi5F0IoQLw\nLIB4OGacWQSgCsBdQogRcMxGs0RKuaiJTRXC8wh7GByj7415CMBVcNTj19bMbxVCFAP4VAgxX0q5\nt/5Kc+bMcd5OTU1FampqEw9DRMXFwC+/SLgcSpK0Bv2j+iPSEOm3uIiI/EWpUOKK6CsQqAl0lsqU\nVJfgeMFxjIgfAaPG6OcI6VLauHEjNm7ceEkeq0XJO4B/ArBKKR+s174GwHs10zh6U05zAEB/D+19\n4ah7b0xfAIV1D3atsavmujccXyxc1E3eicg7R44Ax4/Xa0xai/GJf/BLPERE7UE3Qzf0i+yH9DPp\nzrb9eftRXFXM5L2LqT8gPHfu3DZ7rJbWvN8GYGFDC6WUZVLKX7zYzlIAw4UQibUNQogEACNqljXm\nNIBQIURyvfZhNddnvHh8ImqClMCKFYCUdUbdo3cDxvOsdyeiLi1EF4K+ka7zaxzOP4zccta9U9tp\nafK+H0C0Dx7/HwAyASwRQtwshLgZwBI46tc/qe0khOghhLAKIV6ss+4COA5S/VEI8V9CiDQhxF8B\nvAHgZynlNh/ER9TllZUBO3fWa0xaA4VQYHSP0X6JiYioPfB00Oq+vH04X34edmn3U1TU2bU0eX8C\nwB+FEIGteXApZQWAsQCOAvgSjtH84wDG1iyrJWpiFXXWzQYwFMBuAC8BWA7gQTiS/mtbExcRXZSf\nD/y6u94UkcmrcVXMVQjWBfsnKCKidkCr0mJA1ABolVpn2/mK8zhbehZl1WV+jIw6s5bONvOrEOI+\nANMAvNqaAGqS8Nub6JMJD180pJRHAdzdmscnosalpwPn8+qUzKgqge5bMT7xaf8FRUTUTsQGxaJP\nZB/sydnjbDtw/gDGV41HkDbIj5FRZ9XSkXdIKTOllK1K3ImofausBLZsqdfYYzOgNrPenYgIQHhA\nuFvd+6H8Q8gpuyQzZ1MX1OLknYg6v6Ii4Ndf6zUmrYFaocbI7iP9EhMRUXviab73/bn7UVhZCJvd\n5qeoqDNrcfIuhLje020i6jyysoC9e93r3YfGDYVerfdPUERE7YherUe/yH5QCqWzLaskC/mV+c75\n34l8qTUj73c3cJuIOgGLBdi6FaisrFPvbsgBuu3D+KTx/guMiKgdEULAFGTC5WGXu7QfOn8IRZVF\nfoqKOjOWzRCRR46zqtZrTFoLCGBc4ji/xERE1B5FGiLRL6qfS9vh/MPIKWfdO/kek3ci8ig3F9iz\nx71kJkAVgGGmYZ5XIiLqggK1gW517/vy9qG4qhgWm8VPUVFnxeSdiNzY7cDRo8CRI/UWJK3FiPgR\n0Cg1fomLiKg9MmqMbjPOHL1wFOXV5ax7J59j8k5EbkpLHbPM2O116t2j9gFB51jvTkRUj0qhginY\nhISQBGebhMTh/MMoqCjwX2DUKTF5JyI3+fnA7t31GpPWAGC9OxGRJ5H6SLfSmcMXDiO3PNdPEVFn\nxeSdiNycPg3s3m13bUxejSBtEAbFDPJPUERE7VhoQCj6RPRxaduXuw+l5lJU26r9FBV1RkzeichF\neTlw/Dhw9mydtwelGeixGaO6j4JKofJfcERE7ZRRY8SAbgNc2g6cP4BqezVKzCV+ioo6IybvROSi\nqAjYs6deY/w2QFPJenciogboVDqYgkyINkY72yx2C04UnMCFigt+jIw6m9Yk75YGbhNRB3b6NJCR\n4T5FJMB6dyKixkToIzAgynX0/XD+YeSUcb538p3WJO9P1Ln9x9YGQkT+V13tOFj1l1/rLUheg/CA\ncLeTkBAR0UWRhki3uve9eXtRaa1ElbXKT1FRZ9Pi5F1Kaa1zmyPvRJ1AUZFjfvfysjpTRAbkA9G7\nMSZhDBSClXZERA0xaozoH+U648ze3L2w2W0oNXO+d/INfhITkdO5c8C+ffUak9YCColrk671S0xE\nRB2FXq1H9+DuCNGFONvKLeXILs7G+YrzfoyMOpMWJe9CiCeFEDG+DoaI/MdmA3JzgV9+dZ8iEmC9\nOxFRUxRCgXB9uFvd+9ELR5FbxvneyTdaOvL+EYARQohHhBC9fBkQEflHSYnjzKr79wvXBclrEGuM\nxWVhl/knMCKiDiTKEIW+kX1d2jLyMmC2mVFhqfBTVNSZtCh5l1JapJTfSSk/BTBICPGlEGKqECK6\nyZWJqF3KywMOHgTstjrJe/hhIPg0UhNSIYRoeGUiIgIABGoC3ZL33ed2AxKseyefaNHZVoQQKgDP\nAogHsBfAIgBVAO4SQowAYAWwREq5yFeBElHbkRI4exbYu1cCqJOk15TMXJvMenciIm8YNUYkhyZD\nr9Y7R9oLqwqRV56HvPI8dDN283OE1NG19FSJ/wRglVI+WK99DYD3hBBGACynIeogysoAsxn4aWf9\n5H0NAGBs4lj/BEZE1MGolWqE6ELQL7Ifdp3d5Ww/cuEI4oLiIKXkL5nUKi2teb8NwMKGFkopy6SU\nv7Rw20R0iRUUAOfPA9lZdd4SFBYgYSMSQxLRPbi7/4IjIupgIvQR6Bfpel6MjNwMWO1WlFvK/RQV\ndRYtTd73A2B9O1Encfo0cOBAvUbTDkBbhtSEVH+ERETUYYUFhKFPpOvJmvbk7IEQAiVVJX6KijqL\nlibvTwD4oxAi0JfBENGlV1XlmGXm5188TxF5XfJ1foiKiKjjMmqM6BXeCyrFxerkM6VnUFxVjNxy\nThlJrdPS2WZ+BXAfgGm+DYeILrWiIsBuB3burLegpt49LSHt0gdFRNSBBagDEKgNRN8I11lnjlw4\ngvyKfNilvYE1iZrW4jOsSikzpZSv+jIYIrr0zp51lM2UltR5O9AVArE/o3dEb86MQETUAuEB4egX\nVa/uPScDNmlDWXWZn6KizqDFyTsRdXxWq+NA1YyMegsS1wEKO0fdiYhaKFIfid4RvV3adufshkIo\nUFxV7KeoqDNg8k7UhRUXO+Z435Fuc11QUzLDenciopYxao3oF9kPos70u8cLj8Nis7DunVqFyTtR\nF5abC9hswL699d4KkldDIRScaYaIqIWMGiMCNYG4LOwyl/Yj+UdQUFEAm93WwJpEjWPyTtRF1Z5V\n9ehRwGqtc8KQ0GNAaCYGRA1AiC7EfwESEXVgCqFAiC4EA6IGuLTvyd0DO+yse6cWY/JO1EWVlDhq\n3oOSxFQAACAASURBVHfu8jxFZFoi692JiFoj0hCJvpGuM87U1r0XVhb6KSrq6HyWvAshZgshHhFC\n6Hy1TSJqOxcuACoVsCNdui6oqXe/Pul6P0RFRNR5BOuC3U7WdOj8ISigYN07tZgvR95nA5gP4JQQ\n4m8+3C4RtYHTp4GKCuDUSeXFRmEDEjZAJVQY1WOU/4IjIuoEjBojwnRhMAWZnG02acOxgmMoqiqC\n1W71Y3TUUfkyeU8CMADAiwD6+3C7RORjFRWOy+7d9RbE/QQEFOOqmKtg0Bj8EhsRUWehUWqgV+tx\nZbcrXdp35+yGhESpudRPkVFH5rPkveakTQeklJ9KKe/z1XaJyPcKCwGFAti2w/MUkWMTx/ohKiKi\nzifSEOl2ptU9OXugFEoUVrHunZqPB6wSdUFnzgA6HbBrp3BdUHOw6vXJrHcnIvKFcH2428ma9ubt\nhUapwbnSc36Kijoyr5J3IYReCPG0EGKjECJPCGGpueQKITbULNO3dbBE1HrV1Y6R99OngaLCOm8B\nmhIgbie0Si1S4lP8FyARUSdi1BgRExiD8IBwZ1uVtQqZRZkoqy5Dta3aj9FRR9Rk8i6EiAewF8Dr\nNU3fAnit5vIdAFFze68QonsbxUlEPlJcc1bun36qN8tM4npAacXQuKHQqrSXPjAiok4oQBUAjVKD\ngdEDXdpZ904tpfKiz7sAKgBcLqXM9NRBCJEAYElN39/7KDYiagPnzgFabW29e523ANa7ExH5nBAC\nYQFh6BfZD+tOrnO2787ZjUk9J+FC5QWE68Mb2QKRK2/KZsYDmNlQ4g44DlaFY5aZ8b4Ji4jags0G\n5OQASiWwN0PpujDJkbxPSJ7gh8iIiDqvKEOUW917Rk4GdCodckpz/BQVdVTeJO+y6S4t6ktEl1hJ\nCWC3A3v3ApbqOgerBmcC4b/BoDZgcNxgv8VHRNQZGTVG9AjuAaPG6GwrNhfjTMkZVFgrUGWt8mN0\n1NF4k7yvBfCSECKpoQ5CiEQALwFY46vAiMj3zp8H1GpgR7qHKSIFkGJKgUrhTTUdERF5y6gxQqlQ\nepzvXUCw7p2axZvk/WkAAQCOCiG2CCE+FkK8UnP5WAixBcDRmj5Pt2WwRNQ6Z84ARiOQXv9g1Zop\nIscljfNDVEREnZtSoUSwLhj9o1zPYbk7ZzfUCjXyK/L9FBl1RE0OsUkps4UQVwJ4GMDNAG4FEFqz\nuBDAAQB/AfAPKWVFWwVKRK1TVgaYzYDFAhz/re5L3+6YaQasdyciaiuR+kj0Dnete9+Tswd6tR45\nZTnoE9nHT5FRR+PV7+M1Sfl7NRci6oAKChxnVd25s96C2J8BfQFCtCG4IvoKv8RGRNTZhehCkBSW\nBI1S45zbPbc89/+zd+fhcd/Xfe/fZzbMDAYAsQwAUiTFRaQkiiRIepPt2JEtb2lSJ03zpG6b2O5t\n0iZNnbVN2pvk1nbT2yxNYqe+vcntzdLErdPGaWzfxvFuWRtJiSJIkRJFShR3YifWGWAwy/f+MYMB\nfgMQmAGBGSyf1/Pg4fy+v+8Mjh5xwIPfnN85DCQHCPqDTKYniQQjNY5S1gNNWBXZJG7ehPp6eOZ4\n2nui0CLy7Tvfjs/0I0FEZDXEQjFC/tCCpTM4GEuN1SgyWW/0L7XIJjA1le80U1cHz50sedsXWkQ+\nvlv17iIiq6UuUEckEJl30+qZ3jOEA2EGEgM1ikzWm7KSdzOLmtnPmdkTZtZvZunCV5+ZfbtwLrra\nwYrI8oyOghm8/jrcGZrT3z04ATuOA/A9D3xPjaITEdkc2qJt8/q9n+45TTQYpS/Rh3PquC1LW7Lm\n3cx2AN8G7geeAT4P3CmcbgEOAL8B/JSZvds5d32VYhWRZbp1C8Jh+JuvZIE5yfuu70Bgmng0zoNt\nD9YsPhGRzaAt2sYDLQ/gMx85lwPg2ug1RlOjZF2WZDpJfai+xlHKWlfODaufApLAvrtNWTWzXcAX\nC3t/cIViE5EVkMnk+7u3tsLxkyXJe6FF5Dt2vgMzW/gFRERkRcRCMeqD9TzY+iAXBi8U18/0nuFw\nx2HGp8eVvMuSyimbeQ/wK3dL3AEK5361sFdE1pDRUXAu3yLyzOmS39cL9e7v2aO3rojIaosGo5gZ\nXZ3zhzWFA2F6x3trFJmsJ+Uk75UUYKlYS2SN6e+HUAhefBGmU3Pe8g03IZ6/8vOBB9TfXURktZkZ\nbZE2Hml7xLM+0+99cHJQde+ypHKS928Av2Zme+62wcx2A78GfH2lAhORe+cc3L6dbxH59LMLtIg0\nuK/hPnY3765NgCIim0y8Ps7+tv2etYtDF0mmk2RyGSamJ2oUmawX5dS8/xzwLeCSmR0HzpOfrAr5\nSasHgUeBq4W9IrJGjI/ny2UCATh5suRkoWTmnfe/s/qBiYhsUg11DTTWNbJryy6ujlwFIOdynOs7\nx/7W/Yylxmioa6htkLKmLXnl3Tl3A+gCfgGYBv5O4fEvFB5PA/8C6CrsFZE1YnAQ/H4YGYHXLpXW\nu38DUL27iEg1xUIxcHC086hnfabuvW+ir0aRyXpRVp9351zSOfdp59zjzrlO51xd4auzsPZp51xy\ntYMVkcrcugWxGJw8mcO5Od1kOrshlh8Ionp3EZHqCfgCNIYbORj3TlqdW/c+00ZSZCGasCqyQSWT\nkEjkb1Z95kTGe7LQInJP8x62NWyrQXQiIptXPBqfN6zp/MB50rk0zjnGU+M1ikzWgxVL3s3sR83s\nwyv1eiJyb0ZGwOfL37T63MmSt7rq3UVEaqY50kxLpIWO+o7i2nR2mpcHXsbMGE2N1jA6WetW8sr7\nHxe+RGQNuHULIhG4dg0G++fUuwcmYefTALx393trFJ2IyOYVC8XAFq57jwai6vcui1rJ5P3xwpeI\n1Fg6DUNDEI3CM8dLWkTe/yQEUxjG+x54X20CFBHZxMKBMHX+Og53HPasd/d2EwlGGJ4aJpPL3OXZ\nstmV0yqyLM6576zUa4nIvRkZmX18/ETJjU+FkpkH2x6kLdpWxahERGRGW7RtXt372d6zZHNZHPm6\n9+ZIc42ik7VMN6yKbEB9fRAOQyYDZ06X/I5euFn1u+//7hpEJiIikB/W1BnrpKmuqbiWSCd4bfg1\nAr4Aw1PDizxbNrOVvGH1nWb2rZV6PRFZnlwOenryU1VffNExNemfPVnfCx3nAHjvHtW7i4jUSiwU\nw2c+jnQe8ax393QTCUTonVDduyxsJa+8twOPreDricgyjI3lr7j7fPD0iWnvyb1fBwO/+TWcSUSk\nhqLBKGbGkQ5v8n6m9wzhQJix1BjpbPouz5bNbMmadzPbWeZrqXhWZA0YGIBgMP/45ImSk4V694Pt\nB2kKNyEiIrXhMx8tkRYejj/sWe/u7cY5B8D49DgtkZZahCdrWDk3rF6t4PXcMuMQkRVy82Z+qurY\nGFx6JeQ9uTefvD+267HqByYiIh7xaJydTTsJB8JMZaYAGJoc4sbYDRpCDQwlh5S8yzzlJO9TwJPA\n55fY9wbgn95zRCKybBMTkEpBYyN8+zsZXG7OW7z9HDTkayhVMiMiUnuNdY2YGYfbD/Pc7eeK6929\n3XzPA99DX6KPfa37ahihrEXlJO9ngYxz7v9dbJOZjaDkXaSmhofzte4Az57I4HmLF0pmgr4g7979\n7uoHJyIiHrFQDIAjnUc8yfuZ3jN8/4Pfz0BygFQmRV2grlYhyhpUzg2rp4A3rnYgInLvbt7Md5kB\neP65krd3oUXkkc4jRIPRKkcmIiKlgv4gDaEGDnUc8qx393bnH7h83bvIXOUk778BfMjMbLFNzrnP\nO+cq7l5jZjvM7PNmNmJmo2b2l2a2o4LnP2xmf2FmA2aWNLNXzOynK41DZL1LpWB0NN/f/eZN6Ls9\np97dn8pPVkX17iIia0lbtI3dW3YT8M1+Unpz7CYDiQGC/iCDycEaRidr0ZLJtnPupnPuCTdz6/MK\nMrMo8C1gP/Bh4EeBfcC3C+eWev4bgZNAEPjHwPcAv42GT8kmNHeq6lPPprwndz4NoUkA3rf3fVWM\nSkREFtMSacHv8/Nw2/yuM/XBevV7l3nKqXlfTT8O7Ab2O+deBzCzF4FXydfP/+7dnmhmPuBPga87\n5/7unFPfWb1wRdaunh6IRPKPj5/IeU8W6t0jgQjv2PmOKkcmIiJ3M7fu/Vz/ueJ6d28379v7PkZT\no0ymJ4kEI7UKUdaY5ZS5+M3sW2a2b+7jZX7/DwLHZxJ3AOfcVeAZ4PuXeO5jwEPA7yzze4tsGNks\n9Pfn690zGTjzQtC7oVDvfmzrMd34JCKyhkSCEUL+EF0dXZ71M71nio9V9y5zLae8xMgnzg0lj5fj\nEeD8AusvAweWeO53Ff6MmNkJM5s2sz4z+7SZhZcZj8i6NDqaT+DN4PxLWZKJOR+qRQegM/+PwLt2\nvatGEYqIyN20RlrZ37rfs/bandcYT40T8ocYSAzUKDJZi2pdG94MDC+wfqdwbjHbCn/+d+ArwHuA\n3wR+DPhvKxWgyHrQ1wehwv2pTx+f9p7c/U3w5W9ZUb27iMjaE4/GCflDPNDyQHHN4Tjbd5ZoMEpf\noq+G0claU+vk/V7MxP5nzrmPO+eedM79NvAJ4AfM7KEaxiZSNc7B7dv5qaoAJ0+WbCiUzMRCMd66\n463VDU5ERJYUq8v/AD/aedSzfrrnNAFfgOnsNInpRC1CkzWo1jesDrPwFfYW8lffFzNU+PPrJetf\nB34d6AJeKX3Sxz/+8eLjxx57jMcee6y8SEXWqPFxSKchEMhPWL34UknV2N78W+RN297kaUUmIiJr\nQywUwzCOdBzhL17+i+L6TN27YYxPj1Mfqq9ViLKEJ554gieeeKIq36vW/5K/BBxcYP0A+br3xSxU\nK7+kucm7yEZw5w74/fnHJ55PkcvNuSG17QI03QTQVFURkTXKZz62hLfMaxf58uDLTGWmCAfC9E30\n0RnrrFGEspTSC8Kf+MQnVu171bps5kvAo2a2e2bBzHYBbyucW8zfACngAyXrM8fPr0yIImvb3Kmq\nzxzPek/umf1g6v1731/FqEREpBLx+jgN4Qbua7ivuJbJZTjff55IIMJAcoBVGLkj61Ctk/f/DFwF\nvmhmHzSzDwJfBK4DfzCzyczuN7OMmf3qzJpz7g7w74GfMLN/Z2bvMbN/Bfwq8Cdz20+KbFSTk/my\nmbrCxfZTz5W8pQv17i3hFo5tPVbl6EREpFxN4SYyuQxHOo941s/0nsHv85PJZUikVfcuNU7enXNJ\n4N3AJeDPgM8Cl4F3F87NMPKxWsnzPwn8IvDDwF+TH+z0m+SHP4lseMPD4Cu8i2/dztFzc069uy8N\nu/Izy95835vx+/w1iFBERMoRC8XAzb9ptbu3G8jXvY9OjdYiNFljal3zjnPuBvBDS+y5yl1+0XDO\n/S6LTGIV2ch6eiAazT9+6vgUEJ09uf041E0AqncXEVnrQv4Q0WCUQ+2HPOsv9r1IJpchEojQl+jj\nvsb77vIKsllUfOXdOZehcLV87uOVDkxEFpdOw+AgRAoTs0+cKKmFLJTMALz/AdW7i4isdfH6OK2R\nVloiLcW1ycwkF4cuEglGGEwOknO5GkYoa8Gyymacc0845yZKH4tI9YyO5nu8m+Wnq3afCnk3FFpE\ndtR3cLB9oaZOIiKylrRGW0nn0vPq3rt7uvGZj5zLMTGtlGuzWzJ5NzOfmf1tMzu0yJ5DhT12tz0i\nsrJ6e2dvVD1/YZrEeHD2ZHgYtp0C4C33vQWf1fredBERWUosFAObX/de7PduqnuX8q68/0Pgz4Hx\nRfZMAJ8D/v5KBCUii8vl8sn7TIvIp59Nezfs/ib48h+tPr7n8SpHJyIiyxEJRAj6ghxuP+xZP9N3\nBucc0UCU3oneGkUna0U5yfuPAn9cuGl0Qc65K8AfAh9eobhEZBFjY5DJzA5nOvlcyYY59e7v2/O+\n6gUmIiLLZma0RFrY3rid+uDsNNWRqRGujlwlEogwPDlMNpdd5FVkoysneT8GfLWMfd8E3nRv4YhI\nOQYHIVDoFZVIOC6ej3g3FOrddzTuYF/rvipHJyIiy9Ve304ml6Gro8uzfrr3NGZGjhzj04sVQ8hG\nV07y3gAMl7FvuLBXRFbZrVsQi+Ufn3hhkmxmzlu5+TVovgrAW7a/Rf3dRUTWkVgohsMtOKwJwG9+\nRiZHahGarBHlJO+DwP1l7NtR2CsiqyiRgGQSgoX7U589XvLx6ZySmcd3q95dRGQ9iYXyV2buNqwp\nGozSM9FT9bhk7SgneX8G+EgZ+z4KPH1P0YjIku7cma11B3j+uZIr64WSGYD371V/dxGR9cTv89MU\nbmJP8x6CvtkuYr0TvfRO9BIOhBlLjZHOphd5FdnIyknefxd43Mw+ZWah0pNmFjKzTwGPo0mnIqvu\n9u3ZLjO3e9PcvjZnqqplYNe3AdjbvJedTTtrEKGIiNyLeDROzuXmzeiYufoOqO59E1syeXfOHQd+\nAfhp4KaZfdbM/l3h67PATeBjwM8X9orIKkmlYHgYwuH88VPHU94N9z0HkXwPYNW7i4isT1vCW8i4\nzILDmiBf934neacWockaUNbkFufcp4DHgOeBHwT+deHr7xTWHnPOfXqVYhSRgtGS2RzHT9y93v09\ne95ThYhERGSllVP33pfoq3pcsjYEyt3onHsSeNLM/EBbYXnIOZdZlchEZJ7btyFS6AqZy8HZUwu3\niPSZj/fsVvIuIrIe1QXqiAQiPBJ/BJ/5yLn80L0rI1cYmRphS3gLA4kBprPThPzzKpplg6t4Zrpz\nLuuc6yt8KXEXqZJsFvr7IVoocT//yiTjo3N+aIfG8mUzwIOtD7KtYVsNohQRkZXQFm3DzNjX4p3V\nMdMyEoOx1FgNIpNaqzh5F5HaGB3NX233Fd61zxyf9m7Y/S3w53+ffst9qncXEVnP2qJtpLKpu5bO\nBH1BhpJDtQhNamzJ5N3MvmBmB5faN2f/YTP7q3sLS0RKDQzM9nYHOHHSvBvmtIh87973VikqERFZ\nDbFQDMM4unXh5L0+WE/vRG8tQpMaK+fK+y3gtJl9x8z+iZkdmHvS8g6a2T8zs6eAU8Dt1QhWZLNy\nLl/vPjNVNTGZ5eK5eu+mws2qAV+Ad+96d5UjFBGRlRQNRjEzDrcf9qxfHLxIMp0k6A8ymZlkKjNV\nowilVsppFflTQBfwKvA7wHkzy5jZgJkNAmngReC3CnuOFJ4jIitkYiLfJjJQuMX8xKkkmfScspim\nq9DyGgAH4geI18erH6SIiKwYM6Mt0kY0GPXM7Mi6LC/2vVg8Hk+p3/tmU1a3GefcBeDHzOxngXcA\nbwI6Cqf7yF9tf8o5p79BIqtgaGi21h3g2ZMl94rv/RoUqmjeuv2tqncXEdkA4vVxBicHOdp5lOuj\n14vrZ3rP8Oj2Rwn5QgwkB3TBZpMpu1UkgHNuAvibwpeIVMmtW7NTVQFeOBn0blC9u4jIhtNQ10DO\n5TjaeZQvXvxicb1Y9x6qp2+ijwPxA3d7CdmAyrlh1Wdmf9vMDi2y51Bhj91tj4gsz+QkjI/PTlXt\n6Z/i5pXYnB25fKcZIBwI886d76x+kCIisuJioRg45k1aPd9/nnQ2TcAXIJVNkUwnaxSh1EI5N6z+\nQ+DPgcVKYiaAzwF/fyWCEpFZIyMw99fip0+kvBu2nYJofkz2wfhBWiItVYxORERWS8AXoDHcSFuk\njfb69uJ6KpviwuAFAAxT3fsmU07y/qPAHzvnrt5tg3PuCvCHwIdXKC4RKZg7VRXg2RNZ74Y5JTOP\nbn9U9e4iIhtIPBpnKjs17+r76Z7TANT56+hP9NciNKmRcpL3Y8BXy9j3TfI3sorICslk8jerzkxV\nzeZyvHhq4RaRAB944ANVjE5ERFZbc6SZdDY9b1jTzKTVaDBKf6If51wtwpMaKCd5bwCGy9g3XNgr\nIitkZCQ/VXWmbOalSwlG79TNbghOwPbjQL428tHtj9YgShERWS2xUAyMecn72b6z5FwOv89PJpch\nkU7UKEKptnKS90Hg/jL27SjsFZEV0t8PodDs8TMn0t4Nu56AQH7tcPthtoS3VC84ERFZdeFAmDp/\nHTubdtJY11hcH58e5/Kdy0C+J/zY1FitQpQqKyd5fwb4SBn7Pgo8fU/RiEhRLuedqgpw8mRJQ6c5\n9e7q7y4isjG1RdtIZVJ0dXR51mdaRtb56+hL9NUiNKmBcpL33wUeN7NPmVmo9KSZhczsU8Djhb0i\nsgLGx/M17/5CPp6Ymubii43eTXtmk/f3P/D+KkYnIiLVEq+Pk8qm5pXOzCTv0WCUweQgOZerRXhS\nZUsOaXLOHTezXwB+B/gHZvY14Frh9P3A+4BW4Oedc8dXLVKRTWZwEAJz3qHPnU6STs0pi2m4CfF8\nq7At4S28cdsbqxyhiIhUQyyU/wh2oeTdOYfPfGRdlonpCU9pjWxMZU1Ydc59ysxOA78E/CBQGBfD\nJPAE8OvOuadWJUKRTap0qurTxzPeDXu/DoUqmq6OLv3AFhHZoKLBKGbGg60PUuevI5XNz/sYTA5y\na/wW2xu34zMfo1Oj+rdgEyinbAYA59yTzrnvJd9RZmvhq8k5971K3EVWVjKZ/5q5WdU5R/epkqq1\nOS0i37b9bap3FxHZoHzmoyXSQsZlONTuHXg/UzoTCURU975JLJm8m9kfmdl3zxw757LOub7CV2ax\n54rI8gwPe6eq9gwluPFaSSfW3d8sPnz/XtW7i4hsZPFonKnMFEe3lpTO9Mwm73eSd8jmsgs9XTaQ\ncq68/z3g22Z2xcw+aWYPrHZQIpvdzZslJTMnpnBuTjbfeRpiAwC0R9s53Hm4yhGKiEg1NdY1knXZ\neZNWZ4Y1mRk5ckxMT9QiPKmicpL3TuAfA1eBXwYumdkzZvZPzKxpNYMT2Yymp/PDmSKR2bUTJ0om\n581pEdnV2UVDneajiYhsZDM3rR5uP4zfZsskr49dZzCZH7PjMx/Dk+XM1ZT1bMnk3Tk37pz7Y+fc\nu4DdwK8ALcDvAz1m9t/N7G+ZWdn18yJydyMjMHfKdTqb4dwLMe+mOS0i37b9bQR8Zd17LiIi61TQ\nH6Qh1IDf5+ehtoc852auvtcH61X3vglUlHA756475/5P59zDwKPAHwHvBv4XcMvMfnsVYhTZVHp7\nIRyePb74+gTDA3MuwwcmYefsPDTVu4uIbA5t0TaS6eRd+72HA2GGJ4fJ5HRL4ka27KvlzrnnnHP/\nHNhGfjhTO/CzKxWYyGaUzeaT97n17k8dT3k33f8kBPNr9zXcx8PtD1cxQhERqZWWSAvpXHpe3ftM\n8g6AwXhqvMqRSTUt+7N2M9sHfBj4EfLDmsaA/7FCcYlsSmNjkMuBb86v1c8/V9ICck6LyCOdR4p1\nkCIisrHN/LwvTd5fHXqViekJYqEYAV+A4alhmiPNtQhRqqCiK+9m1mJm/8zMTgAXgX8NvAL8A6DT\nOfdPViFGkU2jvx+CwdnjiVSSi2e3eDfNrXffoXp3EZHNIhKMEPKHaAg1sGfLnuK6w3G272x+TyBC\nz3hPrUKUKiinz3vIzH7QzL4A9ACfAWLkp63ucM59j3Puz51zU6scq8iG5hzcvg2xORfSn+9Okpqc\nk5zX90LHueLh+/eo3l1EZDNpjbQymZmc1+995qbVcCDMxPQE09npWoQnVVDOlfde4PPA24E/AN7k\nnDvonPst55x+tRNZIYkEpFIQmJOrP3087d209+tQaPe+e8tu9rburV6AIiJSc8VhTSU3rZ7uOV18\n7HCqe9/Ayvm8/UngT4C/ds6ll9grIst05453qmrWZTl7KuLdNKdkRvXuIiKbT6xu4br3lwdeJpVJ\nUReoI+gLMjQ5RGu0tRYhyiorp8/7DzjnvrBQ4m5mW8zsjWa2fXXCE9k8btzwlsz03Zng2qslc9Dm\nDGd66463qt5dRGSTiYViGEZHfQdbY1uL6+lcmpcGXgIgGozSO95bqxBllZVT8/5+M/t1s7nXBMHM\nfhnoB54DrpvZ58xMmYTIMkxNwfi4t7/70ycncbk5b7v2c9CQ/2HsM5/q3UVENiGf+dgS3sJUZuqu\nLSND/hDJTJKpjG5H3IjKqXn/CeBB52ZnPprZe4F/C1wg39v994EfRn3eRZZlZMRbMgNw4oTzLsxp\nEbmvZR87t+ysQmQiIrLWxOvjTGXn173P3LQKYJjq3jeocq6UHwV+rWTtHwEp4AMzN60WLsz/feA/\nrGSAIpvB7dveq+6p7BQvnS4pmZlT797V0aV6dxGRTaop3EQml+HY1mOe9bN9Z8nkMgR8AYK+IIPJ\nQeL18RpFKaulnCvv7cBrJWvvBZ4u6TbzZeDBlQpMZLPIZGBgwDtV9dWrCYZ6o7ML/lR+smrB23e+\nXfXuIiKbVCwUAwf3N91Pc3h2GFMyneTVoVcBqA/V0zuhuveNqJzkfRwophWFyaqtwImSfWNAyShI\nEVnK6Gi+x/vcspknn530btr5NITyawFfgPfufm8VIxQRkbUk5A8RDUbJ5DJ3rXsP+AKksikm05ML\nvYSsY+Uk7xeBH5hz/P2FP79Wsm8X0LcCMYlsKn19EArNHudcjheeD3k3zSmZebjtYTobOqsUnYiI\nrEXx+jjJdPKuyTsADsZSY1WOTFZbOZ+7/w7wP82shXxy/hHgHPBMyb6/BZxd2fBENjbnoKcHGhpm\n18amJrh0ttm7cc7Nql0dXTTUNSAiIptXa7SVayPXFrxp1TmHmREOhBlIDNAR66hRlLIayunz/gXy\nXWTeBPwo+XKZH3LO5Wb2mNlW8nXwX16lOEU2pPFxSKfBP6fg7PS5BJOJYPHYooPQOdtBQPXuIiIS\nC8XAYH/rfqLB2XukhqeGuTZ6DYBIMEJfoo85DQNlAyinbAbn3O8553Y65xqcc487514tOd/jnGt1\nzv3B6oQpsjENDnoTd4BnT2Q8x27P18GX/8Fb56/jXbveVa3wRERkjYoEIgR9QQzjcPthz7m5PG8j\nMQAAIABJREFUde/pXJpkOlmLEGWVlJW8i8jquHXLO1U1nZvmxVMlLSDn1Ls/0v4I7fXtVYpORETW\nKjOjJdLCZGZyft17z2zdu2GMT6vf+0ai5F2kRpJJSCS8N6v2j45z9ZWS/u57Z5P3Ix1HVO8uIiIA\ntNe3k8qkOLr17sOawoEwveNqGbmRKHkXqZGREfCVvAOfPpEkl5tdtLaL0HSzePy2HW9TvbuIiAD5\nuneH45H4I55/G25P3KZvIt8AMBqMMjg5qLr3DUTJu0iN3LoFkcjssXOO554zzx6356vFx/XBet5x\n/zuqFZ6IiKxxM5O2w4EwB+IHPOdm6t595iOTyzAxPVH1+GR1KHkXqYHpaRgaguicIaqT2QQXTt+9\nReShjkPEoxpzLSIieX6fn6ZwE1OZKY51HvOcm1s648Onfu8bSNnJu5nVmdkXzOydqxmQyGYwOjp/\n7fKNCfpv1c8u+DKw6zvFQ9W7i4hIqXg0zmR6/k2rp3tPFx+HA2F6J1T3vlGUnbw751LA45U8R0QW\n1tMD4bB37clnvSOsbccJqJv9mPPtO9TfXUREvLaEt5BxGbo6ujBmSy9fH36dkakRIF/3PjQ5RDaX\nrVWYsoIqTcSfBR5djUBENotcDvr6oH7ORfZMLkP38xHPPrfnK8XHW+q28OgOvfVERMRrpu69oa6B\nfS37POfO9uUH35sZzjnVvW8QlSbvPw/8mJl9zMy2m5nfzHxzv1YjSJGNZGwMsllvp5mx1BiXzrZ4\nN85pEXm44zCtkdYqRSgiIutFXaCOSCBCJpeZVzozt+7dzBhNLVCzKetOpcn2OWAP8GngOpAGMnO+\n0isancgGNDAAgZLql+6XxkmOzzZ890XGYNup4nFXZ1fx6oqIiMhcbdE2kukkRzu9/d5nOs4ARANR\n9XvfICotoP3kEufVRFRkCTdvQkPJfafHT+Q8x273N8A3u/b2HW8n6A9WIzwREVln2qJt3Bi7MW9Y\n04WBC0ymJ4kEI0SCEQaTg2RyGd0/tc5V9H/POffxVYpDZFOYmIBUChobZ9emsknOn2707Jtb794W\nbeON295YrRBFRGSdiYViGEZbtI0djTu4MXYDgKzLcq7/HG++780AOBzjqXGaI82LvZyscRXXqJvZ\nNjP7bTN73swuF/78LTPrXI0ARTaS4eH5U1UHRse58nJpf/fZeveuji5aIiX18CIiIgXRYBQzI+dy\n8+re55bOBHwBhqeGqx2erLCKkncz2w+cAT4GTADPAwngZ4CzZrZvkaeLbHo3b3q7zAA8e2qCbGb2\nrRhsuw7NV4vHRzqPqN5dRETuysxoi7QxmZ6cV/c+96bVSCCifu8bQKVX3n8DGAX2O+fe5Zz7kHPu\nMWBfYf03Vzg+kQ0jlcoPZ5rb3z3rsjx/0lu9ltn9Zc/x23a8TfXuIiKyqHh9nKns1Lzk/Vz/OTK5\nDJAf1jSWGiOdVX+R9azS5P1dwP/hnLs6d9E5dw34N4XzIrKAkZH5a8nMOK90e0ti3J6vFh9va9g2\n7yNQERGRUg11DeRcju2N2z2thacyU7wy+Ipn7/j0eLXDkxVUafIeAu72f3yicL4iZrbDzD5vZiNm\nNmpmf2lmO5bxOv/KzHJm9lSlzxWphtu3IeKdw8Tlm6P0Xp9tPWO+LOz6dvG4q6OLLeEt1QpRRETW\nqVgoBi5fQnNs6zHPudM9p4uPAxZgKDlU7fBkBVWavJ8FPlY6jKlw/JPk6+HLZmZR4FvAfuDDwI+S\nL8H5duFcua+zB/gVoB+1q5Q1KJPJ93cvrXd/5sS05zi0/UWIzA7R6OrooiFU0ldSRESkRMAXoDHc\nSCqTWvSm1fpQPX2JvmqHJyuo0kafnwD+GrhgZv8d6AE6gR8mn3R/b4Wv9+PAbvI19K8DmNmLwKvA\nPwV+t8zX+b+BPwMeovL/JpFVNzYGuRyYza6lslO8+II3m0/t/pLn+O071d9dRETKE4/GuT56fV7d\n+9m+s+RcDp/5CPlDDCQHSGVS1AXqahSp3IuKrrw7575CPkEfB34Z+L/IX/EeB77XOffVRZ6+kA8C\nx2cS98L3uAo8A3x/OS9gZv8AOAL8a8DQlXdZg/r6IFiSg49Pj3Gxu9W7uPdrxYe7tuzi4baHqxCd\niIhsBM2RZtLZNHub93q6lI2lxnh9+PXZjU517+tZxX3enXNfcc69EWgEdgKNzrk3LyNxB3gEOL/A\n+svAgaWebGbN5K/O/6JzboHbAUVqz7l8vXuspNvjmQujTIzOXvXwh5Nw33PF466OLg3SEBGRssVC\nMTDw+/x0dXR5zs1tGRn0BxlMDlY7PFkhSybvZjZsZscKj//IzHYDOOcSzrmbzrnEPXz/ZmChaQF3\nCueW8lvAK865/3IPMYisqvFxSKchMKegyznHiRPefcE9z4A/UzxWvbuIiFQiHAhT568jk8vMK53x\n1L0H69XvfR0r58p7FJjpTP1RIL5q0VTAzN5B/gbXn6x1LCKLuXMH/H7vWiIzzoXT3t9Pp3Z9ofjY\nZz7euv2tqncXEZGKtEULw5q2zk/enctXFgf9QVKZFJPpyVqEKPeonJs7rwM/bmYzn+8fM7Pw3TY7\n556s4PsPs/AV9hbyV98X8wfAHwK3zGyml14A8JlZEzDpnJsufdLHP/7x4uPHHnuMxx57rIJwRSq3\n0FTVwYlRLr+0zbs4p979geYH2NeqgcUiIlKZeH2c2+O3OdB2gDp/HalsCoD+RD+3x29zX+N9xb3j\n0+NEgpG7vZRU4IknnuCJJ56oyvcqJ3n/98D/A3ykcPyfFtnrAP8i50u9BBxcYP0A+br3xTxU+PqJ\nBc4NAz8L/F7pibnJu8hqm5yEiQmIl3xedfzUBJn07FulrqWfVMtrxePDnYdV7y4iIhWbuVE16A/y\nSPsjnh7v3b3dxeQ95A8xkBigvb69JnFuNKUXhD/xiU+s2vdaMnl3zv2RmX2FQv914KeBVxZ/Vtm+\nBPwHM9vtnLsCYGa7gLcBv7TEc9+Ft7OMAZ8iXwr0MeDyCsUosmzDw972kADp3DTdp7wfXtneb+T/\nBhcc6TiiencREalYNBjFzHDOcbTzqCd5P9N7hu/b/33FfX2JPh7hkVqFKstUVk9059xt4LaZ/Snw\n5bmtHe/Rfwb+OfBFM/uVwtq/JV+q8wczm8zsfvLJ+Cecc/+2ENN3Sl/MzEYBf4WlOyKr5vZtiJaM\nG0tkxnjltLdF5NSuvyo+9puft9z3FtW7i4hIxXzmoyXSkq97X+Sm1YAvwHR2msR0gvpQfenLyBpW\naZ/3j65g4o5zLgm8G7hEfsjSZ8kn6e8unJthhVht3ouUvCTq8y5rRDoNQ0Pzk/fXe+5w+0rj7II5\n2P3N4uHDbQ+za8uu6gQpIiIbTjwaZyozxaH2Q/hsNtW7NnqNO5OztxQapn7v61DNp5E6524AP7TE\nnquU8YuGc+5dKxSWyD0bHc33eC/17ImM57h+x6skorMdUw93HqYl2rLa4YmIyAbVWNdI1mWpD9Xz\nUOtDvDw4exthd283j+9+HMi3luyb6KMz1lmrUGUZKh7SJCLl6e2FupLJ05OZBC+fbvKsZXd/2XN8\npOOIZzKeiIhIJeb+G3Kk84jn3NxhTZFAhIHkQLGFpKwPSt5FVkEuBz0981tEjk2P8Ep3Sb377tl6\n95A/xBu3vZGQP1SNMEVEZAMK+oM0hBqYzk4vWvfu9/nJ5DIk0vcyb1OqTcm7yCoYG4Nsdv5wprMX\nRxi7M9tpJlCXgu3Hi8cH4wfZ0bijWmGKiMgG1RZtI5lOzrvyfmnoEhPTE8VjwxidGq12eHIPlLyL\nrILBQQiWNIvJ5DI8/5w3m6/f2w2BdPH4cIfq3UVE5N61RFpI59I0R5o9TRByLse5vnPF40ggQl+i\nrwYRynJVnLyb2TEz+yszGzKzrJkdK6z/ezP7wMqHKLL+3Lo1v2RmoRaRk7u+6Dk+0ql6dxERuXdz\n/y1ZrHQmEowwmBwk53JVi03uTUXJu5l9F/As8CDw3/C2bsyx8LRTkU0lkchPVi298j44cYfXznuv\nqk/v/kLxcTQY5UjnEdW7i4jIPYsEI4T8IbK57KLJu8985FzOU0oja1ulV95/HfgqcBD4uZJzp4E3\nrERQIuvZnTvgW+CddbI7QTo1WzYTaR6BttlhxYfbD3Nfw33VCFFERDaB1kgrk5n5w5peGniJ6ex0\n8dhMde/rSaXJ+zHg951b8LOVQSB+7yGJrG8LlcxMZSd58VSDZy38wHHPZ1eHOw7THGmuQoQiIrIZ\nzAxr2tqwlY76juL6dHaalwdme79HA1F6J3prEaIsQ6XJ+xQQucu5TkC/tsmmlkrByAiEw971ifQo\nF0taRE7c/z88x10dXTTUeRN8ERGR5YrVlVn3HogwPDlMNpetWmyyfJUm708DP2tmnsmsZmbAPwa+\ntVKBiaxHo6NgNn/9av8ANy83etbSu/66+LiprolDHYdU7y4iIismFophGM65eS0j5ybvZkaOHOPT\n49UOUZah0uT9V8nXtZ8FfqWw9mHg28BbgU+sXGgi68/t2/OvuudcjpPP5XBuNqvfsvMGxAaKx10d\nXRpPLSIiK8pnPraEtzCVmeLY1mOec2d7z3qutPvNz8jkSLVDlGWoKHl3zp0F3gH0Ar9cWP7ngAPe\n6Zx75W7PFdnoslno759f757MjHOhpEWkf893PMeHOg7RElF/dxERWVnx+jhT2Sl2b9lNU11TcT2R\nTvDqnVeLx9FglJ6JnlqEKBWquM+7c+60c+5xoBHYATQ5597lnOte4qkiG9roKORy88tmRqbvzKt3\nH73/v3qOj3QcUb27iIisuKZwE5lcBjObVzpzpvdM8XE4EGYsNUY6my59CVljyk7ezazOzL5gZu8E\ncM5NOuduOecSqxeeyPrR3z+/tzvA+UvDDA/M3ucdCGXIbP928bgt2saDbQ+q3l1ERFZcLBTL10ew\n+E2rM1T3vvaVnbw751LA45U8R2SzcA56eiBWMhx1Opui+wVvEXzzAxchmCoeH+k4onp3ERFZFSF/\niGgwSjqbXvDKu3OueOw3P3eSd6odolSo0kT8WeDR1QhEZD2bmMi3iQwEvOuJzBivnG7zrOX2fM1z\nrHp3ERFZTfH6OMl0kofaHiIcmL2gNDQ5xI2xG8XjaDBKX6KvFiFKBSpN3n8e+DEz+5iZbTczv5n5\n5n6tRpAia93Q0MJTVfsSfbx2zpuYD+/4rOdY/d1FRGQ1tUZbmc5OE/AFONxx2HPudM/p4uO6QB3j\nqXHP9FVZeypNts8Be4BPA9eBNJCZ86W7HGRTunVrfsmMc47TL6ZITc5ejq9vSpJrn/1Bua1hG3tb\n9qreXUREVk0sFCtO9C6te5970yoABmOpsSpFJssRWHqLxyeXOO+WOC+y4UxOwvg4xOPe9WR2gpdO\nNXnWmh46Q2JON5quji7a69urEKWIiGxWkUCEoC9INpdd8qbVoC/IUHKItqi35FPWjoqSd+fcx1cp\nDpF1a2Rk4amqY9PDXOz2/vCb3vX/eY4PdxymNeJtIykiIrKSzIyWSAtjqTEOth8k4AuQyWUAuDV+\ni4HEAPH6/BWo+mA9vRO9PNj2YC1DlkVUXKNuZtvM7LfN7Hkzu1z487fMTO0yZFO6fRsikfnr1wf7\nufaq98r74PY/9Rwfbj+sencREVl17fXtpDIpwoEwD7c97Dk39+p70B9kMjPJVGaq2iFKmSpK3s1s\nP3AG+BgwATwPJICfAc6a2b4Vj1BkDUunYXAQolHveiaX5oVTAVxu9pJ8645BaLhdPN69ZTc7t+xU\nvbuIiKy6WCiGK1Q3l7aMXLDfe0r93teqSq+8/wYwCuwvTFX9kHPuMWBfYf03Vzg+kTVtdDTf4720\nbGYiM8Yrp73lMPX7nvMcH+44rHp3ERGpilhotqvCsc5jnnOlyXvIF2IgOVCVuKRylSbv7wL+D+fc\n1bmLzrlrwL8pnBfZNPr6oK5u/vpwqp+LZ7z17on7/6fn+FDHIdW7i4hIVfh9fprCTUxlpujq7MKY\nvep0+c5lT4eZ+lA9fRPq975WVZq8h4C7fY4yUTgvsinkcvmpqvX188+9cmWEod7ZQvhAIMfQ1v9W\nPDaMrvYuz5UQERGR1RSPxplMT9JY18jelr3FdYfjbO/Z4nHAFyCVTZFMJ2sRpiyh0uT9LPCx0mFM\nheOfJF8PL7IpjI9DJgN+v3d9MpPgxRcaPWsdD16H0GTxeH/rfjpiHdQFFrhsLyIisgq2hLeQcfku\nM0u1jDRMde9rVKXJ+yeA9wAXzOyTZvaTZvYJ4CXgfYXzIpvC4CAEFmi2OpEendciMrj3ac/x4Y7D\ndMQ6VjM8ERERj7mf9pbetFo6rKnOX0dfQqUza1Glfd6/YmbfC/wa8Mvk53U54AXge51zX135EEXW\nnlwObt6cP1UVoDfZw6UzXZ610R2f8xwfbj9MS6RlNUMUERHxqAvUEQlEyOQy8668vzz4MlOZKcKB\nMADRYJSBxADOOWyhYSZSMxX3eXfOfcU590agEdgJNDrn3qzEXTaLoSF4+mlIpSAY9J7LuiznzmeY\nSs7+XlzfOM1w698Uj/3m50D7ARpC6u8uIiLV1RZtI5lO0l7fzn0N9xXXM7kM5/vPF4/9Pj+ZXIZE\nOlGLMGURFSfvM5xzCefcTeec/q/KppBMQnc3nDyZr3NvXaBRTCI9xoXT3ivqHQcugc8Vjw/EDxCP\nxlXvLiIiVdcWbSOVTQFl1L2bMTY1hqwtlQ5p+ldm9h/vcu73zOxfrkxYImtHJgOXL8OTT8LwMHR0\nQDi88N7R9BAXu+Pexd3f8hyqv7uIiNRKLBQrtolcaliT6t7XpkqvvH8UOHeXc2eBf3RP0YisIc5B\nfz889RS89hq0tEBT0+LPuT7Ux9VXvJuGtv8Xz/Hh9sO0RtXfXUREqi8ajGJm5Fxu3pX3c33nyOQy\nnr2DyUFyLlftMGURlSbvO4FLdzn3OrDrnqIRWSPGx+HUKXjhhfwQpra2+S0hS01lJznXHSGXm72x\np337OKP1p4vHIX+Ih+IPqd5dRERqwsxoi7QxmZ5kZ9NOT/OEycwkF4cuFo995iPrskxMT9QiVLmL\nSpP3JLD9LufuA1L3Fo5IbU1Pwyuv5G9ITSSgvX3hCaoLSWTG5rWIbHnovOf4UPshmsPNqncXEZGa\nidfHmcpOYWbzS2d6vKUzPvMxOjVazfBkCZUm708B/8LMPBW/heNfKJwXWXecg9u383Xt169DPA4N\nFV4cH5zqmZe8p3f9jedY9e4iIlJrDXUNxVKYY53HPOdK694jgYjq3teYivq8Ax8HjgMXzey/AjfJ\nX4n/EaAV1bzLOjQyAi+/DKOj0Nw8v/1jOXIux+Ub4/TfihbX/P4cfVv/1LPvUPyQ6t1FRKSmYqFY\nfkoPCw9ryrkcPstf340EItxJ3iGby+L3LVE/KlVR0ZV359xZ4DHgGvCLwGcKf74OfLdz7szdny2y\ntkxNwfnzcPx4vqNMe/vyEneAZGacl083e9Z2PDTEhP9a8TgajLKvbZ/q3UVEpKYCvgCN4UZSmRT7\nWvZRH6wvnhtNjXJ15Grx2MzIkVPd+xqynCFNzznn3kl+SNMO8kOaHnPOnVrx6ERWQS6XL4158kno\n7c2XyNTXL/28xYylh+eVzDTse8FzfKTjCI2hRtW7i4hIzcWjcSYzk/h9fro6vFPBS0tnfOZjeHK4\nmuHJIirt8+43syCAcy7pnLsFvNPMfsHMji7xdJGam5mO+vLL+baPzc2wElOf+yZucemMtxxm8v4v\neY4PdRyiI9Zx799MRETkHjVHmkln08DCpTNz1QfrVfe+hlRa8/45YAr4MICZ/QTwnwrn0mb2fc65\nr69gfCIrIpmES5fyN6U2NeWvtq+U6WyKVy4ZifHZmpv6hgy3Wj5XrCkEOBg/qHp3ERFZE2KhGIVZ\nTRzb6r1p9XTPac9xOBBmIDFAJpch4Ks0dZSVVmnZzFuAue0z/iXwh8AW4H8C//sKxSWyIuZOR71z\nZ/HpqMu1UIvI+w/eYtKNFI8b6xrZ07JH9e4iIrImhANh6vx1ZHIZDsQPEPKHiuf6En30jPd4n2Aw\nnhqvcpSykEqT93byHWYws33AbuAzzrkx4E+Awysancg9qHQ66nLdSfVz8bT3Un7d3hOe42Odx4gF\nY6p3FxGRNaMtmh/WFPKHeCT+iOdcad17wBdgeEp172tBpcn7GDBzifG7gaFCBxqALLDC1zRFKjc+\nDs8/X9l01OVyznFrdIDXL3h/Mxjf8Ree40Mdh2iPqb+7iIisHfH6OKlsfr7mvGFNC/R7n3c1Xmqi\n0uT9WeCXzOz7gJ8Dvjzn3F4KV+VFauFepqMuVzI7wcUzDWQzs2+l9m1TXI982bPvYPwgbdG20qeL\niIjUTCwUKz4+2untO1J602o4EGZieoLp7HRVYpO7qzR5/yXyw5i+BNSRH9o040PkBziJVNVKTEdd\nrvH0CBfPeJPy7QevMJ2bLB63RlrZ2bTT80NSRESk1qLBKGaGc47DHYeLg5kAroxcmdce0uG4NnKt\n9GWkyiod0nTJObcPaHfOPeCcuzLn9M+Qv4FVpGpGRvJDls6ehVgMWltXpvVjuYameubdrOrb8x3P\n8Ru2voFIMEI4oKoyERFZO3zmoyXSwmRmklgoxv7W/Z7zpVffWyOtvDr0Kq/dea2aYUqJioc0ATjn\nBhdYe9E5N3DvIYksbSWnoy5XJpfmak+SnuuzE558PsedbX/u2XeoXf3dRURkbYpH40xlpoD5pTML\nDWuK18e5NHhJCXwNLSt5F6mV1ZiOulwTmTEudXv7tu95aJzrzls9dqjjEK0R9XcXEZG1p7GukazL\nAksPawIl8GuBOu3LujE0BC+9lB+41NwMgRr/7R2dHuRit7dFZPzhS7zmZm/m2RbbRkesg4Y69XcX\nEZG1Z+79WEc6vMn7xaGLJNNJosGoZ31uAg/wQMsDqx+oFOnKu6x5ySScOQPPPZdv+RiP1z5xz+TS\n9CVvc7Hkyntul3fA8LGtx4gEVO8uIiJrU9AfpCHUwHR2mtZovsHCjKzL8mLfiws+T1fga0fJu6xZ\npdNR29tXfjpqpdK5aW4lr/DC0BNcuxxmfHS20D4SzdLX+nnPftW7i4jIWtcWbSOZTgJLt4ycSwl8\nbSh5lzVpZjrq5curOx21XNPZFDcSl3lh8Dtcn3iV+kAjV1/c4dmzv2uY69NnPWsH2w+q3l1ERNa0\nlkgL6VwamJ+8n+49vehzlcBXn2reZU0ZH88PWhoczCfsjY21jSeVnaJv8jo3k1cwMxqDW/BZflzr\n+Re8wTU/eI4c2eLxri27aIm2qN5dRETWtMWGNb3U/xLT2WlC/tBdn68a+OpS8i5rwvQ0vP46XLkC\nkUi+RKaWprJJeievczt5Fb8FaAq2eIZXTKeMS+e8Q5emdvy15/hYp+rdRURk7YsEI4T8IbK5LNsa\nttFe305/oh+AVDbFhcELdHV0LfoaSuCrR2UzUlOl01Hb2qo3HXUhk5kEl8de4vTQk/RN3qQp2Epj\nsNmTuANcPBcjnZ5da+tMcSP6vzx7DrUfor2+xr+FiIiIlKE10spkZhIzm9cysrun+y7P8lIJTXXo\nyrvUzMgIvPwyjI7mWz9We8jSXMnMBLeSV+ifvEnIH2JLsA1bZFTruZKSmQePDPJs6pJn7UD7Adqi\n3umrIiIia1E8Gqd3opdYKMbRzqN87fLXiucWu2m1lK7Arz4l71J1U1Pw2mtw40Z+wFItS2Qm0mPc\nTr7OYKqHoK+O5lB80aR9Rmm9e+SBUzhc8Xh/634a6xpV7y4iIutCrO7ude9n+s6QzWXx+/xlvZYS\n+NWl5F2qJpeDmzfzN6T6fPl+7WXkyatiPD3CzcRl7qT6qfOHaQ6V/xvEyJ0AN16fHVhh5kjc9wVI\nze55w9Y3qN5dRETWjVgohmE459jTvIfGukbGUmMATExPcHn4Mvtb95f9ekrgV4+Sd6mKO3fg/Pna\nT0cdmx7mRvI1RlKD1PkjtNRVdtl/dDjAn35mp2dt9/4kr/Mdz9rB+EHVu4uIyLrhMx9bwluYykwR\nCUbo6ujiqetPFc9393ZXlLzPvKYS+JWnG1ZlVc1MRz15snbTUZ1zjEwPcW74BOdGTpDKJGmpa6c+\nUH5JSy4L3/hinF/86CM8/2Sz59zeI730pa4Uj/3m58H4g6p3FxGRdSVeH2cqOwVUNqxpMbqJdeXp\nyrusikwm3z3m0iUIhWpT1+6cYzQ9xLXxiyQyY4T99bRUUB4z48qlKH/y6Z28frF+3rlAMEdD19dh\nenbtQPwA9cF61buLiMi60hRuIpPLAPOT9+7ebpxzZd0XVkpX4FeWkndZcQMD+RKZdDo/HdVf3v0t\nKybncoxMD3Jt4hKT2XEi/gaaKyyPAUhM+Pn8H2/jm1+K49z8H1Y7dif5337uOl91X/Ek78e2HiMc\nCKveXURE1pVYKMZM74WH2h6izl9HKpu/oWswOcit8Vtsb9y+rNdWAr9ylLzLiqn1dNScy3En1c/1\nxEUmM0nqA40V3Yg6wzl49lstfO73tzM6PL9/ZTiS5Qc/cpv3/kA/gQB85uQpz/mD8YN01Hcs+79D\nRESkFkL+EPWhetLZNEF/kEPthzjVM/tv3Ome08tO3kEJ/EpR8i73rNbTUbMuy52pPq4lLjKdnaI+\n2FTxjagzbl0L81/+4w4unFn4N483vXOYH/nJG7TE0wD0Jm4xlL5VPB/0BdnXuk/17iIisi61Rdvo\nGe+hyd/E0a1HPcl7d283H3zwg/f0+krg752Sd1k256C3F156Kd8Gsq0t3wKyWjK5DEOpXq4lLpHJ\nTRMLNFEfWN7l/tSU8cX/upUv/0UH2cz8/4j2rSk+/LHrdL15zLN+dtB71f1wx2HCgXD+o0cREZF1\npjXayrWRa8DK3bRaSgn8vVHyLssyOppP2msxHTWTSzMw1cONxKtkXJqGwBYCgaZlv94CNkEBAAAg\nAElEQVTp40382Wd2MNhXN+9cIJjj+z7Uy9/+UC+hOjfv/Ctj3uT9DVvfQMgfIhKMLDseERGRWomF\nYsWbUg+1H8JvfrIuC8CNsRsMJgdX5NNlJfDLp+RdKjI1BZcv5zvJVHs6ajo3zcDUbW4kXiWHI+Zv\nIuBb/l/hwb4Qf/aZHZw+vmXB8wffMMZHPnadzu2pBc8757iUeN6z9kj8ETpjncuOSUREpJYigQgB\nX4BsLkskGOGhtod4aeCl4vkzvWd4z573rMj3UgK/PErepSy5HNy6BRcuVH866nQ2Rd/UTW4lXsfh\naAg24bfl/9XNpI2/+Xw7X/jsNqZT80tktrRO8yM/eZM3f/fwov+N18evMZYZLB5HAhH2NO+hNdq6\n7NhERERqycxoibQwlhojFopxtPOoJ3nv7u1eseQdlMAvx5pI3s1sB/C7wHsAA74B/Kxz7sYSz3sT\n8BPAO4D7gEHgKeBXnHNXVzPmzeTOnXyJTCJR3emoqewUfZPXuZm8gpnRGNyCz+6t7+SFMzH+5Pd2\ncvv6/LIW8zne9wP9/N2P3CZSn1vytc4Nea+6H+08SsAfoCGk/u4iIrJ+tde3M5AcIEaMI51H+Oy5\nzxbPdfd0r/j3UwJfmZon72YWBb4FTAIfLiz/GvBtMzvsnEsu8vQfBh4GPg2cI5/A/ypwysyOOOdu\nrl7kG18ymR+y1NMDDQ35q+3VMJVN0jt5nZ7kNXzmpynYgs/u7U7Y0eEAn/uD7TzzjYWvij9wYIKP\n/vR17n9gsuzXfGXcm7wf23qMOn+d6t1FRGRdi4ViOJe/z+tI5xHPuVfvvMp4anzFBxEqgS9fzZN3\n4MeB3cB+59zrAGb2IvAq8E/JX5G/m990zg3MXTCzZ4Arhdf9N6sS8QZXq+mok5kEt5NX6Zu6gd8C\nNK5A0p7Lwrf+V5y/+KNtJBPz/7rHGjL8vR+/yTs/MFRRp5x0JseriRc8awfbD6reXURE1r25HdO2\nhLewp3kPrw+/DoDDcbbvLN+187tW/PsqgS/PWkjePwgcn0ncAZxzVwtJ+PezSPJemrgX1q6b2QCw\nbTWC3cicy09HffnlfO/2ak1HTWYmuJW8Qv/kTUL+EFuCbcsav1zqyqUof/Lpnbx+sX7B8+/8wCAf\n+vGbNDRly35NB4yNwc3JV0lmR4vrDaEGdjTuUL27iIise36fn6ZwE1OZKcKBMEc7jxaTd4BTt0+t\nSvIOSuDLsRaS90eAv1pg/WXghyp9MTN7GGgHLtxjXJvK0FB+OurYWH46akMVyrYn0mPcTr7OYKqH\noK+O5lB8RZL2xISfz//xNr75pTjOzX+9HbuTfPRnrrP/YKKy101Aahru2waXRrwlM2/Y+gZ8Pp/q\n3UVEZEOIR+NcGb5STN7/8sJfFs999txnybosP/WmnyIcCK/491YCv7i1kLw3A8MLrN8pnCubmQWA\n3wf6gT+899A2vpGRfHnM0FA+Ya9GicxEepQbide4k+qnzh+mObQy39Q5ePZbLXzu97czOjy/8Xw4\nkuUHP3Kb9/5Af0U33aamZ2/WPfAIRCPwqa94+7ur3l1ERDaSLeEtZFwGKFygMh85N9vM4XPnP8eJ\nmyf45GOf5OH4wyv+/ZXA391aSN5X0meAR4H/v707j4+7qvc//jqZmexbk+5tmpaWorQWKD8QuUCL\nXhZFERQFFaFcaK+4gd7f1etVSlsu/sR73dEHi8gmIIIoyI/FH0tAoChYKVvtnqXNvmcymcks5/fH\nN2mTzjTrdzKZ9P18PPJIO99lzvQ0yTtnPuec86y1HcOdfCTr6oKdO6GhAXJzJya0d/a2URPYRXuo\nmSxPDiVZ7j3p/qps7v5ZGdveSLzD6klntHHp1TWUzAiP+J7RKHR2QXY2LFsGxX37QEViEbbUbxl0\n7rKZy5iVP2vM7RcREZlMBta9z8ibwZdO+hI3//VmLAc3LNzbvpc1j65h7cq1rDl+zbj2XklEAT6x\nyRDe20g8wl6CM/o+IsaY7+FMUr3MWvvM4c7bsGHDgT+vXr2a1atXj/QppoRAAPbsgZoaJ5QmO7Rb\na+kIt1LTvZPOcBs5GbmuhvZQ0PDofXN44qFZRCPxM05nzglx+VeqWXFy54jvGbPQ1QkYWLzYWWUn\nY0D1zav7XiUQPrgIUklOCXPz57qy45yIiMhkkOXNIsebQyQWwZvh5fLjLmf5zOVsqNhAnb/uwHlR\nG+WWv93CSzUvsXHVRsqLy11tR7oE+IqKCioqKibkuUz/UkCpYox5Fsi01p5+yOMVgLXWnjmCe3wb\nuAH4srX2F0OcZ1P9elMlFILKSti7F3w+p649mZssOaG9hWr/DvzhDrI9eeR4E08cHastm4u49+Yy\nmhuy4o55fTE+ekk9H7uknsyskfe53w+9YZg3z6lt9w2ovtnVuovbt9zOs3ufHXTN2YvP5munfI1V\n5atUNiMiIlPGO43v0NjdSGHWwXe1/b1+frD5B/xxxx/jzs/yZHHtKddy0XsvcmUO20AxG6Opu4ml\n05dO2gA/kDEGm2jinQsmw8j7Y8D/GGMWWWv3AhhjFgKnAt8c7mJjzFdxgvt/DhXcj1ThsDPKvmuX\nE9ZLSxnVkoijFbMx2nubqfLvoCfaRY6ngGkujrQDNDdkcu/NZWzZXJzw+PITO7n8K9XMnh8a8T1D\nIfB3w/RSKC+HnAEZfG/bXm7bchvP7Hlm0NuF/U6ee7Lq3UVEZMqZnjudms7B+2XmZ+Zz/arrWVW+\nihv/fCNtwYPTFkPREDe9fBMvVr3IdWdcx8w8937+p8sI/ESYDCPvucBWnE2avtP38A1AHnBgkyZj\nTDmwG9horb2h77FLgPuBp4CNOLuz9uuw1g5aceZIGnmPRqG2FrZvh1gMiouTu+xjzMZoDTVS3b2d\nnkiAPG8hWR53Z6BHwoYnH57JH349l95Q/G8gxaW9XHr1Pk5e1TbidxUiEaeuPS8XFh0FRQNK5ivb\nK/nlll/y9O6nE4Z2gPfPez8bV2+krKiMY2ccO5aXJSIiMil193bzUvVLhy0Lbe1p5cY/38gLVS/E\nHSvMKuQ//uk/OHvx2a62KV1G4JM58p7y8A5gjCnDWc/9LJwA/gxwrbW2esA5C4E9wAZr7aa+x+7E\n2ZU10T9OhbX2g4c8z5QP77GYMwn1H/9w1movLmZUK6uMVtRGaQ02UNW9nd5okDxfEZkZ8WUs47Xt\njXzu+ukCaqvjR7dNhuXsCxr55OW15OTFElwdr7+u3WTAokUwffrBuvaajhpu33I7T+1+atDM+oGO\nKT2GL5z4BU5bcBrNgWZWzl3p6giDiIhIqllreWbPMxRnFx9200RrLX/c8Ud+sPkHdIfjl2A+Z/E5\nfOPUb1CUXeRau9IhwE/58D5RpnJ4t9ZZ7nHbNmdZw6IiZ3fUZInEIrSE6qnq3kEk1ku+twhfhvtP\n2NHm5YFb5/PyM4k3P1pyrJ81X62mfEnPiO/p90M4AvPnw9w5B3+52de5jzv+fgdP7HyCqE28cdPS\nkqWsO3Edq8pXHajnawo0qd5dRESmpC21W/D3+snLHHreWm1XLRsqNsStxgbOmvHXr7qeU+af4lq7\nJnuAV3h3yVQN721tzkh7ezsUFjqryCRLJBamKVhHTfdOIjZMgbcYb0b8murjFYvCc4/P4KFfzSXQ\nHf/WQX5BhIvX7uOMc1tGXMMfDEKgx6n7L19wsK69tquWO/5+B4/vePywoX3xtMWsO3EdZy48c9Do\nQyQWIRAOcOaiYedVi4iIpJ2ajhq2NW+jNGf4HcSjsSj3v30/v3jtF4Rj8Uszf+rYT3HN+69xbWOn\nyRzgFd5dMtXCe2ens8FSUxPk5TkfyRKO9dIUrKWmeycxLPmeItfXc+23d0cud/1kAXu2J35BZ5zb\nzCVr91FQlDhoH+pAXXueUyLTX9de76/njr/fwWPbHztsaD+q+CjWrlzLh476UMK3DDtDnczMm8my\nmctG9uJERETSSHuwnVf3vcqM3BkjvmZX6y7WP7+eHa074o4tKFrAptWbWD5zuSvtm6wBXuHdJVMl\nvHd3w+7dsH+/M3pcUJC85+qNhmgI7mN/9x4slgJfER6TnNDe7ffw8J1zefaxGST6/162KMCaa6pZ\nujy+pi6R/rr2DA8sWgilfXXtDf4G7tp6F7//x++JxCIJry0vKmftyrWcddRZeDIOP9O3qbtJ9e4i\nIjJlRWIRnt3z7Kj3MglHw9y25Tbu3np33Pwxj/FwxfFXcNXKq1wZCJyMAV7h3SXpHt6DQWeDpaoq\nyMpy6tqTJRQN0tBTw77AHowxFHqLyTDJWa7GWnjluRIeuGU+HW3xJTjZOVE+cXktZ13QOOLJt4nq\n2pu6m7hr6108su2RhG/nASwoXMBVK6/inMXnDBna+zUHmjm9/HRyfbkja5iIiEiaeaXmFbDOxk2j\n9Ub9G1xfcT37u/bHHXvv9PeyafUmFk1bNO42TrYAr/DuknQN7729TmDfs8dZ7rG4OHkbLAWjPdT3\nVFEXqCLDeMj3Fh12hrkb9ldlc/fPytj2RmHC4yed0calV9dQMiNx2D5Uf1379OlOXXt2thOw7956\nN49se4RQNPHa7/MK5rF25VrOXXLuiEcBVO8uIiJHgp0tO6nuqKY4O/H+KsMJhAP86NUf8ft//D7u\nWJYniy+f/GUuXnbxuPPGZArwCu8uSbfwHok4pTE7djij08lcq70n0k1toJKGYA0e4016aA8FDY/e\nN4cnHppFNBL/PDPnBrn8yzWsOLlzRPcLR6CrC/Lznbr2wgJn/dl7tt7DQ+8+dNjQPjd/LleuvJLz\njj5v1G/dqd5dRESOBM2BZl7f/zoz8kZe957IS9UvccOLN9DS0xJ37OS5J7N+1Xpm588e13NMlgCv\n8O6SdAnvsRjU1zsryITDyV2rPRDxsz+wl8aefWR6MsnzFLm+pfGhtmwu4t6by2huiH/7zeuL8dFL\n6vnYJfVkZg3fV9EY+LucuvajFjkryXSE2rln6z389t3fEowEE143O382V55wJR89+qP4PGNbLac5\n0MwJc05QvbuIiExpwUiQisqKUU1aPZz2YDvffem7PLf3ubhj+Zn5fOPUb/DhJR8eVxaZDAFe4d0l\nkz28W+usHLNtG/T0OKHd5/4qjAD4w53UBvbQHKrDl5FFnqcw6aG9uSGTe28uY8vmxG+7LT+xk8u/\nUs3s+YlHyQ/V5XfenVhQBrNngz/Szn1v3ceD7zxIIBxIeM2svFlccfwVnH/M+WR6xrcufVOgiTPK\nz1C9u4iITHnP732eXF+uKxNMrbU8uetJvv/K9/H3+uOOf2jRh/jWad8ac5kOpD7AK7y7ZDKH95YW\nZ6S9szO5a7X7wx3UdO+iNdRIliebPG/iWnM3RcKGJx+eyR9+PZfeUHyJTHFpL5devY+TV7WNqJa/\nJ+j8cjNzBpSVQa/p5L637uM3b/8m4e5u4GwQseb4NVxwzAVjmnAzUDASpCvURVF2ER8o+8C47iUi\nIpIO3mp4i5ZACwVZ7i1xV++vZ+MLG3mt9rW4Y6U5pVx3xnWctuC0Md8/lQFe4d0lkzG8t7c7Ne0t\nLc6SjzlJ2qSzs7eNmsAu2kPNZHlyyPMmcX3JAba9kc9dP11AbXX8CzMZlrMvaOSTl9eSkxdLcPVg\nA+vaj1oEZHZx/9v3c/9b9x82tJfmlLLm+DVc+J4Lx70pRHdvN93hbgqyClgybQkz8maMaEUaERGR\ndFfvr2dr/dZRLxk5nJiN8eA7D3LzX29OOD/twvdcyNdO+dqY3+VOVYBXeHfJZArvfj/s2gV1dZCb\n6wRSt1lr6Qi3UtO9k85wGzkZueR4k/BECXS0eXng1vm8/EziHdmWHOtnzVerKV/SM+y9ojFnvXav\nFxYdBVn5fh585zfc99Z9dPV2JbymJKeEy1ZcxkXHXjSu0G6tpTPUSSgaojS3lMXTFlOSU5L0EiMR\nEZHJxN/r5+Xql10P7/32tu1lfcV6tjVvizs2v3A+G1dv5LhZx43p3qkI8ArvLpkM4b2nx1nysbo6\neWu1O6G9hWr/DrrCHeR48sjxJnH71QFiUXju8Rk89Ku5BLrj6+LyCyJcvHYfZ5zbQsYIFrPp8kM0\nAmULoKC0m99te5D73rqPjlBHwvOLs4u5bMVlfOrYT5HjG/vbGDEboz3YTtRGmVcwj/Licgqzkl9i\nJCIiMhnFbIxn9jxDSXbyBrAisQi/3PJL7nzjzridzzNMBpcfdznrVq4b00ITEx3gFd5dksrwHgod\nXKvd53NCu9v/92M2RntvM1X+HfREu8jxFJDtSVIdTgJ7d+Ry108WsGd74l8Uzji3mUvW7qOgKJrw\n+ECBHmfN9lkzoWR2gMd2/ZZ737z3sKG9KKuIz6/4PJ9e9ulxTSCNxCIHnmNh0ULmF83XhFQRERHg\n9drX6Qx2UpSdxF0igbcb32Z9xXqqO6rjji0tWcqmMzeNKYBPZIBXeHdJKsJ7OAw1NU6JjDHOCjIj\nGXEejZiN0RpqpLp7Oz2RAHneQrI8SZrxmkC338PDd87l2cdmkOj/admiAGuuqWbp8sR16QP1hp2S\nosICmFMW5Knqh7l76920BdsSnl+QWcClKy7l4mUXk5859pKg/kmomZ5MFpcsZk7BnHGvRiMiIjKV\ndIY6ebfxXdpCbWR7spP6jnQwEuQnf/kJD737UNwxX4aPL570RT67/LOjnns2UQFe4d0lExneo1Go\nrYXt251125OxwVLURmkNNlDVvZ3eaJA8XxGZGeNbSWU0rIVXnivhgVvm09EW/xZWdk6UT1xey1kX\nNA67Tn1/XbvPB3MWBKloeIR7tt6dcCMHcNaC/dz7Psdnln9mXKFdk1BFRERGpyPYwd72vdT56/Bl\n+CjKSt7Gjq/ue5WNL2ykKdAUd2zl7JVsWL2BuQVzR3XPiQjwCu8umYjwHotBY6OzVnsoBNOmub/B\nUiQWoSVUT1X3DiKxXvK9RfgyJnaUeH9VNnf/rIxtbyT+rfukM9q49OoaSmaEh7yPxRlpj0Vh1rwQ\nm9t/zz1v3kVzoDnh+Xm+PD77vs/y2eWfHfNyVZqEKiIiMn7dvd3UdNZQ1V6FMYairCJX1oE/VEew\ng5teuYk/7f5T3LE8Xx7/9oF/42NLPzaqn+PJDvAK7y5JZni39uBa7V1dTk17lsuD4JFYmKZgHTXd\nO4nYMAXeYrwZSdrF6TBCQcOj983hiYdmEY3E/5Y9c26Qy79cw4qTO4e9V39de8mMXt7oeZR7376T\nxu7GhOfm+nK5ZNklfO59nxtzrV00FqUj1KFJqCIiIi4KRULUddWxu2034ViYwszCce+pksjTu5/m\nppdvojMUnzFWla/i26d/m5KckhHfL5kBXuHdJckK721tTnlMW1tyNlgKx3ppCtZS072TGJZ8T3J+\nsx3Ols1F3HtzGc0N8V+QXl+Mj11Sz0cvqScza+h/4/669rz8MO/GHuO+d39FQ3dDwnOzvdlcsuwS\nLl1x6Zh3WtMkVBERkeSLxCI0+hvZ2bqTQDhAfma+6z9vG7sb2fTCJl7d/2rcsWnZ0/j26d9m9cLV\nI75fsgK8wrtL3A7vnZ3OBktNTZCX53y4qTcaoiG4j/3de7BYCnxFeMzEh/bmhkzuvbmMLZsTh+fl\nJ3Zy+VeqmT0/fnOFgaJR512JDG+EXZ7H+c32O6jz1yU8N8uTxaeXfZrPr/j8qH6LHkiTUEVERCae\ntZaWnhZ2texKyuRWay0Pb3uYH7/644QbO52/9Hy+/oGvj3hOXDICvMK7S9wK793dsHs37N/v7Iha\n4PJmpaFokIaeavYF9mKModBbTIaZ+EmUkbDhyYdn8odfz6U3FF8iU1zay6VX7+PkVW1DLntpcUJ7\nJBpht+cJHt5zB7Vd+xOem+XJ4pPv/SSXHXfZmDeC0CRUERGRySGZk1ur2qtYX7Ged5reiTs2N38u\nG1ZvYOWclSO6l9sBXuHdJeMN78Ggs057VVVyNlgKRgPU91RTF6jCmAwKvMVJm709nG1v5HPXTxdQ\nWx2/TrzJsJx9QSOfvLyWnLzYkPcJ9EB3T4TdPMWjNXewr6sm4XmZnkw+8Z5PsOb4NWMK7ZqEKiIi\nMnkla3JrJBbhrjfu4vYtt8dt7GQwXLriUr5w4hdGVIPvZoBXeHfJWMN7b+/BtdozMpwVZNzMhD2R\nbmoDlTQEa/AYL/ne5C25NJyONi8P3Dqfl58pTXh8ybF+1ny1mvIlPUPeJ9QLXf4o28J/4sn626np\nit9oAZy1Wi94zwVccfwVzMybOer2ahKqiIhI+kjW5NZtTdu4ruI6Ktsr444tnraYTWdu4pjSY4a9\nj1sBXuHdJaMN75GIUxqzY4ezmozba7UHIn72B/bS2LOPTE8meZ6ipI4Ux6LQ3uajrdlHa1Om87nZ\n+dzWnElrs4+Wxkwi4fhfHPILIly8dh9nnNsy5CZT0Si0d8Z4u/sZnm65jZquyoTneTO8fPyYj3PF\n8VcwO3/2qF+LJqGKiIikr2RMbg1Ggvz8tZ/zwNsPxB3zZnj51xP/lctWXDZsKa0bAV7h3SUjDe+x\nGNTXO8s+hsNOaHdzrXZ/uJPawB6aQ3X4MrLI8xSOO7T3howTwJt8tLZkJgzo7a0+bGz0z3PGuc1c\nsnYfBUXRw55jgY7OGG+0P8f/a7uN6q49Cc/zGA/nH3M+/3L8vzCnYM6o26JJqCIiIlNHMia3vrb/\nNTa8sCHhSnYrZq1g0+pNzC+cP+Q9xhvgFd5dMlx4t9ZZOWbbNujpcUK7z8Vl1LvC7ezr3k1rqJEs\nTzZ53uH/c1oL3V2eASPkiUfN/V3ur0JTtijAmmuqWbq8e8jz/N2W15oreLbtNqr8OxOe4zEezjv6\nPK484UrmFc4bdVv8vX4C4YAmoYqIiExRbk5u7Qp18d+v/DdP7Hoi7liON4evnfI1LnzPhUMOno4n\nwCu8u2So8N6/wVJnp/trtXf2tlET2EV7qJksTw55Xmd5mmgUOlp9tA4oW+kfPW9rzqStxQnq4d6J\nrX+fVtrLuRc1cPaFjUO+4xAMWf5S/yJ/ar2V6u4dCc/JMBl8eMmHueqEqygrKhtVOzQJVURE5Mjj\n5uTWZ/c+y3f//N0DpbYDnVZ2Gt854ztDLpQx1gCv8O6SROG9o8PZYKmlBfLzIdelsumeHsvu2na2\n1dSxvz5MoLWArpZ8J5A3Z9LW5KO9bWxlLOORXxihZEYv00rDzufpYaZN76Wk//OMMLl50SEn5EYi\nllfrXubJplupDmxLeI7BcO6Sc7nqhKsoLy4fVRs1CVVERETcmtzaHGjmhhdv4OWal+OOFWUV8Z+n\n/ycfWvShw14/lgCv8O6SgeHd73dWj6mrcwJ7/sjW8cdaJ/A3NUFDg/O5sdH5cB6zNDZauromdrQ8\nI8MyrdQJ4NNmhCmZPjigl0zvpXh6mMzMsfd3zFr+UruZ/1t/K1U98WuqghPaz1p8FmtPWMuiaYtG\ndX9NQhUREZFDRWNRGvwN45rcaq3lD9v/wA83/5CeSPyKeR9Z8hH+/dR/pyAr8eY9ow3wCu8uMcbY\nQMCyZw9UV8ev1R6JOCPwg8N4fEAPDb2RqOsys6MHR8YHjJD3/71kei+FxRGSVQJureX1+r/wWO2t\nVPa8ddjz/nnRP7N25VoWlywe1f01CVVERESG48bk1n2d+1hfsZ43G96MOzYrbxYbVm3gpHknJbx2\nNAFe4d0lxhj7ve9Z2tudHT8PDeUtLc5KMxOpoCh8YGR82vTwgZKWA3+ePnwZS7JYa9na8jqPVN3C\n3p6thz3vzIVnsm7lOo4uPXpU99ckVBERERmL8Uxujcai3Pvmvdzyt1uIxCJxxz+z/DN86aQvke2N\nnwA50gCv8O4SY4x1FjVMPo/HUlwSZtqAspX+zyUzwkwrHX8ZSzK90/o3Hq68lV2BLYc9Z1X5Ktad\nuG5Emx700yRUERERcct4Jrdub9nO+ufXs7ttd9yxRcWL2LR6E++d8d64YyMJ8ArvLnErvGfnxCgq\n7aGoNMj0GVFKZ0TjRs8LiyNDbmY00SxOvT4DPg/8p7AWeqLd7PG/w2M1d7LD/9ph73X6gtNZt3Jd\nwv/Qh6NJqCIiIpIsY53cGoqEuOVvt/DrN3+NPSQjeoyHtSvXsub4NXG/EAwX4BXeXTKS8F5cDKWl\nMH2689H/5+KSCN7CJkL52/FmBynwFeMxI18E3pgBoXmCz+8f1DYZYInRHm6gIVRJQ7CKhl7nc32w\nkvZw05DPcWrZqaxbuY7lM5ePuF2ahCoiIiITZayTW7fUbeH6iuup89fFHVs+czkbV22MWz1vqACv\n8O4SY4w97jjL7NkwaxbMns2BP8+Z43zOzubAiHlGBvRGe2nqqaWyaycWe/DtGDMgFPd/ZsDfDz1m\nDjnO4HNGeo/+c4a7RyAcoKazmsr2Sqo6qpzP7VVUdVQRio5uxu0p805h3YnrWDFrxYiv0SRUERER\nSZWxTG719/r54eYf8tiOx+KOZXmyuPaUa7novRcNKvU9XIBXeHfJcDusDhSKhNjXuY89bXsGh/ZJ\nxFpLQ3fDoHBe2eF8TrQl8GidPPdk1p24juNnHz/iazQJVURERCaTjmAHle2V1Pnr8GZ4h53cWlFZ\nwY1/vpG2YFvcsVPmncL6VeuZmTfzwGOJArzCu0tGEt6DkSDV7dXsbd+LMYbirOKUh89gJEh1x+BR\n9Mr2Sqo7qhOuVTpWHuNhXuE8jik9houOvYgT55w4ous0CVVEREQmu9FMbm3taeXGP9/IC1UvxB0r\nzCrkm//0Tc5ZfM6Bxw4N8ArvLhkqvAfCAao7qqlqryLDZFCcXTziJYfcYK2lOdAcF9CrOqoS1l+N\nR2FWIQuLFlJeXE55UTkLixeysHgh8wrm4fOMvI5fk1BFREQk3Yx0cqu1lsd3Ps7/vPI/dIe7446f\nvfhsvnnqNynKdjYNGhjgjy49WuHdDYnCe3dvN5XtldR01ozorZTxCkVC1HTWDFJIS3cAAA3oSURB\nVArn/bXoif5jjFWGyXACdV84Ly8uPxDYp2VPG9eouCahioiISLobNLk1EiDfl3hya21XLRsqNrCl\nPn757Bm5M1h/xno+UPYB4GCAP++Y8xTe3TAwvHeFuqhsr2Rf1z4yPZkUZRa5VubRP0ni0BH0qvYq\nartq45YiGo/8zPyDAX3AKPr8wvmuTxDVJFQRERGZakYyuTVmY9z/1v38/LWfE46F4+7xqWM/xVdP\n/io5vhy6e7tZtWiVwrsbjDG2I9jBnrY91HXVkeXNojCzcMyhvTfay77OfYMCev9nf6/fvXZjmFsw\nd9DoeX9YL80pTXptuSahioiIyJFguMmtu1p3sf759exo3RF37YLCBWw6cxOLihcpvLvFGGOf2PEE\n2d7hlwvqZ62lLdg2aCWX/oC+v2s/MRtzrX15vjzKi8rjatHLCstGtNGAmzQJVURERI5UQ01uDUfD\n3LblNu7eendcDvQYD5973+e45xP3KLy7wRhjX9//esJjkVgkbhS9P7B3hjpdbcec/DmDylz6P0/P\nnZ7ScGytJRwL4+/1axKqiIiIHPGGmty6tWEr659fz/6u/fEXbkDh3Q3GGPvM7mfiSlwq2yvZ37mf\nqI269lw53pxBI+j9nxcULSDbm+3a84xUJBYZ/GEjcZvNZpgMcn25zMidoUmoIiIiIn36J7fuat1F\nd6T7wOTWQDjAj1/9MY/845HBF2xQeHeFMcaywd17zsqbdaAWfeCqLjPzZk7IKLq19kAgD8fCRGKR\nA2/hWCymbwvWTE8mOb4c8nx55HhzyMvMw+fxkenJxJfhw+fxTbpNqEREREQmk8NNbn2p+iVuePEG\nWnpanBM3KLy7YqzhPcuTlXAUvbyonBxfjuvt7BeNRQeF8v7R8oGh3BjjhHFfHjm+HHJ9uWR7swcF\n80xPpmrVRURERFx06ORWay03vXwTz+59VuHdLcOF95l5MxMG9Fn5s1xf+z0SixCOhg+E80QTX30e\nHzleJ5D3j5pnejKdUN4XzjVaLiIiIpI6Aye3AlRUVvC9s76n8O4GY4zNvCGTBUULBi25uLDIqUXP\ny8wb93PEbOxAKO8P5gNZazHGkO3NJi8z78CoebYve1AJS6Ync0J3eBURERGRseuf3LqzdSdnLzlb\n4d0Nxhj7l5q/jHmN8oETPsPRcMIJrh7jITcz90Bdea43lyxvFj6P70AJi8/jG+9LEREREZFJKBqL\n4vV4kxbej7iai0TBvX+JxIGlLDFiB+rK+2V7s8nx5lCUXeSUsnhzBpWw+DJ82rxIRERE5AiW7Cx4\nxIX35p7mwy6RmOfLIzcv9+CkzwEj5b4MnyZ9ioiIiEhKHXFlM7WdtVoiUURERESSxhijmnc3GGPs\nkfR6RURERGTiJTO8azkTEREREZE0ofAuIiIiIpImFN5FRERERNKEwruIiIiISJpQeBcRERERSRMK\n7yIiIiIiaULhXUREREQkTSi8i4iIiIikCYV3EREREZE0ofAuIiIiIpImFN5FRERERNKEwruIiIiI\nSJpQeBcRERERSRMK7yIiIiIiaULhXUREREQkTSi8i4iIiIikiZSHd2NMmTHmYWNMuzGmwxjzO2NM\n2QivzTbG/Lcxps4YEzDGvGKMOT3ZbRYRERERSYWUhndjTC7wHLAUuAz4PHA08HzfseHcAVwFfAc4\nD6gDnjbGHJecFouIiIiIpI6x1qbuyY25BvgBsNRau6fvsYXATuAb1tofDXHtccDfgSustXf3PeYB\n3gG2W2s/nuAam8rXKyIiIiJTnzEGa61Jxr1TXTZzPrC5P7gDWGsrgZeBuPCd4Now8OCAa6PAb4Bz\njDE+11srKVNRUZHqJsg4qP/Sl/ouvan/0pv6TxJJdXhfBryd4PF3gWNHcO0ea20wwbWZwJLxN08m\nC30DS2/qv/Slvktv6r/0pv6TRFId3qcBbQkeb+07NpSSIa7tPz7puP2FONb7jea64c4d6/HRPj4Z\nuNm2ydB3IzlnqvTfVPzaG+6c0R6brH0H6dd/+tob7Ej73jmWY5O1/9Lta28k56b7116qw/sRR18E\nY398MjjSfgANdTzd+m8qfu0Nd47Cu/v3U3gfmyPte6fCu/v3U245KNUTVuuB31trrz7k8V8An7TW\nzhri2geB46y17znk8U/j1L0vs9ZuO+SYZquKiIiISNIla8KqNxk3HYV3gOUJHj8Wp3Z9uGsvMMZk\nH1L3fizQC+w69IJk/SOKiIiIiEyEVJfNPAacYoxZ1P9A31KRp/YdG+5aH/DpAdd6gYuBp621Ybcb\nKyIiIiKSSqkum8kFtgI9OBstAdwA5AErrLWBvvPKgd3ARmvtDQOufwA4B/h3oBK4GvgIcKq19o0J\nehkiIiIiIhMipSPvfeH8g8AO4F7g1zgh/YP9wb2PwWnroWUvVwB3Av8FPA7MA85VcBcRERGRqSil\nI++TTd8Ora8PeMiHU0O/wlqbaD16mWSMMZnA93HegQkCexPttiuTkzGmAlgAdPQ99Dtr7X+lrkUy\nWsaYK4A7gAustcOVP8ok0bcIxHuAGBAFrrPWPpnaVslwjDHFOAOfR+NUMTQAX7TW7k5pw2TEjDHX\nAZfi7E/0CWvto8Ndk+oJq5NK3w6tJ/T/3RhzCfBNBfe08l3Aa61dCmCMmZni9sjoWOBahb701Ddn\n6Spgc2pbImOwzlrbAWCMOR540RhTZDXCN9lZ4IfW2ucAjDFfAX4JnJnSVslo/AnnF7Bf4fTnsFI9\nYXWyW4vzRSBpoG8OxVrgP/ofs9Y2pq5FMkZaFSoNGWMygNuBr+Cs+CVppD+49ykGWhTcJz9rbUd/\ncO+zGViYoubIGFhr/2Kt3TuaaxTeD8MYsxg4Bee3IUkPS3B23f2WMeavxpiXjDHnpbpRMmrfN8a8\naYx5yBizNNWNkRH7OvCStXZLqhsiY2OM+aExZjfwB+CyVLdHxuRanP6TKSztwrsxZr4x5mfGmM3G\nmIAxJmaMWXCYc8uMMQ8bY9qNMR3GmN8ZY8pG+FRXAg8fMhoh45Tk/vPi1EvvtNaejNOHvxq4FKmM\nzwR8/V1mrT3GWrsCeAL4U9+IroxTMvvOGLMc+ARw48CHXX4JR7SJ+Nlnrf26tXYxzpLLDxpjCtx+\nHUeiicotxpjrcUbdv+Vi8494E5g7RywdfyguAT4FtAAvHu6kvhKK54ClOCMIn8eZ0PF83zGMMZ83\nxvy97+PqAdd6+q5RyYz7ktl/VTj1Yr8GsNZux1mK9IQETyFjk9SvP2ttTf89rLV3Avk4v5DJ+CWr\n774InIYTGnYaY/bivGt5W98xcUfSf/b1s9Y+DXT3XSfjNxG55TvAucCHD9m4UsZvwr72Rsxam1Yf\n9K2Q0/fnq3Bmxi9IcN41QAQ4asBjC4Ew8LVhnuN84B+pfq1T8SPZ/YczWvuxvj/PAeqAo1P9uqfK\nRzL7D8gCpg/4+0eAesCT6tc9FT4m4nvngPOfB85P9WueSh9J/trLBhYN+PsHcFYtyUv1654KHxPw\nc+964FWgMNWvdSp+TNT3TqAC+PhI2pR2I++27xWOwPnAZmvtngHXVgIvA8MtHXglGnVPignov6uB\na4wxbwJPAV+31u4cY3PlEEnuv0LgKePUu78B/G/go9ZZBUrGaYK+d0qSJLn/coD7jDFvGWP+Dvwf\n4EJrbfc4mix9ktl3xphlOOG9BHihb0T3r+NrsQyU7O+dxpgNxpga4P3AL40x1caYuUM90VReKnIZ\n8PsEj78LXDTUhVbrgk8GY+o/a20V8M/JapSM2Kj7z1rbBPyvZDZKRmTM3zv7WWu1TF3qjOVrrw04\nNZmNkhEZS9+9Q3qWQE9FY80tG4ANo3miqdzh03BWHjlUa98xmdzUf+lN/Ze+1HfpTf2XvtR36W3C\n+m8qh3cRERERkSllKof3NhL/plOC81uQTG7qv/Sm/ktf6rv0pv5LX+q79DZh/TeVw/s7wPIEjx+L\nU38kk5v6L72p/9KX+i69qf/Sl/ouvU1Y/03l8P4YcMrADXqMMQtxJuU8lqI2ycip/9Kb+i99qe/S\nm/ovfanv0tuE9Z8Z+Qo4k4cxpn/W7oeAfwW+CDQDjdbaF/vOycXZoKcH+E7f+TcAecAKa21gQhst\nB6j/0pv6L32p79Kb+i99qe/S22Trv3QN77EBf7Uc3Ia7wlr7wQHnlQE/As7qO+cZ4FprbfVEtVXi\nqf/Sm/ovfanv0pv6L32p79LbZOu/tAzvIiIiIiJHoqlc8y4iIiIiMqUovIuIiIiIpAmFdxERERGR\nNKHwLiIiIiKSJhTeRURERETShMK7iIiIiEiaUHgXEREREUkTCu8iIiIiImlC4V1EREREJE0ovIuI\niIiIpAmFdxERERGRNKHwLiIigxhjVhpjXjXG+I0xdxtjvH2P+4wx16S6fSIiRzKFdxEROcAYsxD4\nKfAT4GNABPhm3+Hzgd+mpGEiIgKAN9UNEBGRSeUzwIettV19f3/eGHNL359LrbV1KWqXiIgAxlqb\n6jaIiMgkZoz5ItAK/N1auz3V7REROZKpbEZERIZTA5yn4C4iknoK7yIiMpwg8GyqGyEiIgrvIiIy\nvGOBp1LdCBERUXgXEZHhnWqtrU91I0REROFdRESGYIzJx1kuUkREJgGFdxERGcoxwKOpboSIiDi0\nVKSIiIiISJrQyLuIiIiISJpQeBcRERERSRMK7yIiIiIiaULhXUREREQkTSi8i4iIiIikCYV3ERER\nEZE0ofAuIiIiIpImFN5FRERERNLE/wcp+50sc52aJQAAAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0x109605f10>"
       ]
      }
     ],
     "prompt_number": 31
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Grid-search"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "You can search the (hyper) parameter space and find the best hyperparameters using grid search in scikit-learn"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn.grid_search import GridSearchCV"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc_params = {\n",
      "    'C': np.logspace(-1, 2, 4),\n",
      "    'gamma': np.logspace(-4, 0, 5),\n",
      "}"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gs_svc = GridSearchCV(SVC(), svc_params, cv=3, n_jobs=4)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gs_svc.fit(X, y)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gs_svc.best_params_, gs_svc.best_score_"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Exercise: predicting the quality of a wine given a set of physicochemical measurements"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Two datasets were created, using red and white wine samples.\n",
      "The inputs include objective tests (e.g. PH values) and the output is based on sensory data\n",
      "(median of at least 3 evaluations made by wine experts). Each expert graded the wine quality \n",
      "between 0 (very bad) and 10 (very excellent).\n",
      "\n",
      "P. Cortez, A. Cerdeira, F. Almeida, T. Matos and J. Reis. \n",
      "Modeling wine preferences by data mining from physicochemical properties.\n",
      "In Decision Support Systems, Elsevier, 47(4):547-553. ISSN: 0167-9236.\n",
      "\n",
      "This dataset is available from the UC Irvine Machine Learning Repo [http://archive.ics.uci.edu/ml/datasets/Wine+Quality](http://archive.ics.uci.edu/ml/datasets/Wine+Quality)"
     ]
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "You can try several **classification** approaches for the quality (10 discrete classes for `quality`) or you can try (using either statsmodels or sklearn)\n",
      "**regressions** approaches: e.g. predicting the alcohol content given the other (or subset thereof) measurements."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "wine  = pd.read_csv('./data/winequality-red.csv', sep=';')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 51
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "wine.head()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "html": [
        "<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
        "<table border=\"1\" class=\"dataframe\">\n",
        "  <thead>\n",
        "    <tr style=\"text-align: right;\">\n",
        "      <th></th>\n",
        "      <th>fixed acidity</th>\n",
        "      <th>volatile acidity</th>\n",
        "      <th>citric acid</th>\n",
        "      <th>residual sugar</th>\n",
        "      <th>chlorides</th>\n",
        "      <th>free sulfur dioxide</th>\n",
        "      <th>total sulfur dioxide</th>\n",
        "      <th>density</th>\n",
        "      <th>pH</th>\n",
        "      <th>sulphates</th>\n",
        "      <th>alcohol</th>\n",
        "      <th>quality</th>\n",
        "    </tr>\n",
        "  </thead>\n",
        "  <tbody>\n",
        "    <tr>\n",
        "      <th>0</th>\n",
        "      <td>  7.4</td>\n",
        "      <td> 0.70</td>\n",
        "      <td> 0.00</td>\n",
        "      <td> 1.9</td>\n",
        "      <td> 0.076</td>\n",
        "      <td> 11</td>\n",
        "      <td> 34</td>\n",
        "      <td> 0.9978</td>\n",
        "      <td> 3.51</td>\n",
        "      <td> 0.56</td>\n",
        "      <td> 9.4</td>\n",
        "      <td> 5</td>\n",
        "    </tr>\n",
        "    <tr>\n",
        "      <th>1</th>\n",
        "      <td>  7.8</td>\n",
        "      <td> 0.88</td>\n",
        "      <td> 0.00</td>\n",
        "      <td> 2.6</td>\n",
        "      <td> 0.098</td>\n",
        "      <td> 25</td>\n",
        "      <td> 67</td>\n",
        "      <td> 0.9968</td>\n",
        "      <td> 3.20</td>\n",
        "      <td> 0.68</td>\n",
        "      <td> 9.8</td>\n",
        "      <td> 5</td>\n",
        "    </tr>\n",
        "    <tr>\n",
        "      <th>2</th>\n",
        "      <td>  7.8</td>\n",
        "      <td> 0.76</td>\n",
        "      <td> 0.04</td>\n",
        "      <td> 2.3</td>\n",
        "      <td> 0.092</td>\n",
        "      <td> 15</td>\n",
        "      <td> 54</td>\n",
        "      <td> 0.9970</td>\n",
        "      <td> 3.26</td>\n",
        "      <td> 0.65</td>\n",
        "      <td> 9.8</td>\n",
        "      <td> 5</td>\n",
        "    </tr>\n",
        "    <tr>\n",
        "      <th>3</th>\n",
        "      <td> 11.2</td>\n",
        "      <td> 0.28</td>\n",
        "      <td> 0.56</td>\n",
        "      <td> 1.9</td>\n",
        "      <td> 0.075</td>\n",
        "      <td> 17</td>\n",
        "      <td> 60</td>\n",
        "      <td> 0.9980</td>\n",
        "      <td> 3.16</td>\n",
        "      <td> 0.58</td>\n",
        "      <td> 9.8</td>\n",
        "      <td> 6</td>\n",
        "    </tr>\n",
        "    <tr>\n",
        "      <th>4</th>\n",
        "      <td>  7.4</td>\n",
        "      <td> 0.70</td>\n",
        "      <td> 0.00</td>\n",
        "      <td> 1.9</td>\n",
        "      <td> 0.076</td>\n",
        "      <td> 11</td>\n",
        "      <td> 34</td>\n",
        "      <td> 0.9978</td>\n",
        "      <td> 3.51</td>\n",
        "      <td> 0.56</td>\n",
        "      <td> 9.4</td>\n",
        "      <td> 5</td>\n",
        "    </tr>\n",
        "  </tbody>\n",
        "</table>\n",
        "</div>"
       ],
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 52,
       "text": [
        "   fixed acidity  volatile acidity  citric acid  residual sugar  chlorides  \\\n",
        "0            7.4              0.70         0.00             1.9      0.076   \n",
        "1            7.8              0.88         0.00             2.6      0.098   \n",
        "2            7.8              0.76         0.04             2.3      0.092   \n",
        "3           11.2              0.28         0.56             1.9      0.075   \n",
        "4            7.4              0.70         0.00             1.9      0.076   \n",
        "\n",
        "   free sulfur dioxide  total sulfur dioxide  density    pH  sulphates  \\\n",
        "0                   11                    34   0.9978  3.51       0.56   \n",
        "1                   25                    67   0.9968  3.20       0.68   \n",
        "2                   15                    54   0.9970  3.26       0.65   \n",
        "3                   17                    60   0.9980  3.16       0.58   \n",
        "4                   11                    34   0.9978  3.51       0.56   \n",
        "\n",
        "   alcohol  quality  \n",
        "0      9.4        5  \n",
        "1      9.8        5  \n",
        "2      9.8        5  \n",
        "3      9.8        6  \n",
        "4      9.4        5  "
       ]
      }
     ],
     "prompt_number": 52
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Below an example of classification (using the same SVC classifier)  \n",
      "\n",
      "you need to add the cross-validation step"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "quality = wine.pop('quality')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 53
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "y = quality.values"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 54
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "X = wine.values"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 55
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn.preprocessing import StandardScaler as scaler"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 56
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "scaler = scaler()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 57
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "scaler.fit(X)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 58,
       "text": [
        "StandardScaler(copy=True, with_mean=True, with_std=True)"
       ]
      }
     ],
     "prompt_number": 58
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "Xscaled = scaler.transform(X)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 59
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn.svm import SVC"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 60
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc = SVC()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 61
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc.fit(Xscaled, y)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 62,
       "text": [
        "SVC(C=1.0, cache_size=200, class_weight=None, coef0=0.0, degree=3, gamma=0.0,\n",
        "  kernel='rbf', max_iter=-1, probability=False, random_state=None,\n",
        "  shrinking=True, tol=0.001, verbose=False)"
       ]
      }
     ],
     "prompt_number": 62
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "y_hat = svc.predict(Xscaled)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 63
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "y_hat "
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 64,
       "text": [
        "array([5, 5, 5, ..., 6, 6, 6])"
       ]
      }
     ],
     "prompt_number": 64
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "y"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 65,
       "text": [
        "array([5, 5, 5, ..., 6, 5, 6])"
       ]
      }
     ],
     "prompt_number": 65
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "svc.score(X, y)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 66,
       "text": [
        "0.42589118198874298"
       ]
      }
     ],
     "prompt_number": 66
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from sklearn.metrics import confusion_matrix"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 67
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "confusion_matrix(y, y_hat)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 68,
       "text": [
        "array([[  0,   0,   7,   3,   0,   0],\n",
        "       [  0,   1,  38,  13,   1,   0],\n",
        "       [  0,   0, 549, 129,   3,   0],\n",
        "       [  0,   0, 169, 450,  19,   0],\n",
        "       [  0,   0,  10, 115,  74,   0],\n",
        "       [  0,   0,   0,  12,   6,   0]])"
       ]
      }
     ],
     "prompt_number": 68
    }
   ],
   "metadata": {}
  }
 ]
}