bigartm_examples.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# BigArtm Experiments"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Данный ноутбук показывает возможности интерфейса для создания экспериментов над тематическими моделями в библиотеке BigArtm."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/opt/anaconda/lib/python2.7/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.\n",
      "  warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')\n"
     ]
    }
   ],
   "source": [
    "from artm_experiments import Experiment, Pool, OptimizationTopicsFilter\n",
    "from artm import *\n",
    "import glob\n",
    "import os\n",
    "%pylab inline\n",
    "%load_ext autoreload\n",
    "%autoreload 2"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Как обычно создаем нужную нам модель, над которой и будем проводить эксперименты."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "batch_vectorizer = BatchVectorizer(data_path='kos', data_format='batches')\n",
    "\n",
    "# эта функция создает новую модель со случайным начальным seed'ом\n",
    "def create_model():\n",
    "    model_artm = ARTM(num_topics=50,\n",
    "                      scores=[PerplexityScore(name='PerplexityScore',\n",
    "                                              use_unigram_document_model=False,\n",
    "                                              dictionary_name='dictionary')],\n",
    "                      regularizers=[SmoothSparseThetaRegularizer(name='SparseTheta', tau=-0.2),\n",
    "                                    DecorrelatorPhiRegularizer(name='DecorrelatorPhi', tau=2.5e+5)])\n",
    "    \n",
    "    if not os.path.exists('kos/dictionary.dict'):\n",
    "        model_artm.gather_dictionary(dictionary_target_name='dictionary', \n",
    "                                     data_path='kos', vocab_file_path='vocab.kos.txt')\n",
    "        model_artm.save_dictionary(dictionary_name='dictionary', dictionary_path='kos/dictionary.dict')\n",
    "    else:\n",
    "        model_artm.load_dictionary(dictionary_name='dictionary', dictionary_path='kos/dictionary.dict')\n",
    "    model_artm.initialize(dictionary_name='dictionary', seed=np.random.randint(1000000))\n",
    "\n",
    "    return model_artm"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Теперь создаем эксперимент, используя оптимизационную выпуклую оболчку. Затем 5 раз создаем тематические модели и передаем созаднные темы в наш класс эксперимента. Для понимания процесса обработки тем можно нарисовать гистограммы расстояний тем до выпуклой оболчки. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEPCAYAAABY9lNGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGHZJREFUeJzt3Xu0JWV95vHvQ3O/NXijiSCNM1HUEZEY4wRH9+AYEY04\nTjSCGnQYdWWJMmaWS2Rpumc0g+gSx1vWJMFRMGK8jFEw4abQEkxUEJBWsTUxON66ExGhlSgqv/mj\n6jTHps7pfc7e++zq3t/PWntRVXtX1e/sLs5z6n2r3kpVIUnS9nabdgGSpH4yICRJnQwISVInA0KS\n1MmAkCR1MiAkSZ1WJCCS7Jbk+iQXtfMHJ7k8yaYklyVZvRJ1SJKGt1JnEGcAX5k3fybwyap6KHAl\n8JoVqkOSNKSJB0SSw4ATgfPmLT4JOL+dPh945qTrkCQtzUqcQbwVeBUw/5btQ6pqC0BVbQYesAJ1\nSJKWYKIBkeRpwJaquhHIIh91vA9J6pndJ7z944BnJDkR2Ac4IMn7gM1JDqmqLUnWAP/UtXISg0OS\nlqGqFvujfCgTPYOoqrOq6kFV9WDgucCVVfUC4GLghe3HTgU+vsg2evVat27d1GvYGWrqa13WZE2z\nUNe4TOs+iDcCT06yCXhSOy9J6pFJNzFtU1WfBj7dTv8A+A8rtW9J0tJ5J/USDQaDaZdwL32sCfpZ\nlzUNx5qG19e6xiHjbK8atyTV5/okqY+SUH3vpJYk7bwMCElSJwNCktTJgJAkdTIgJEmdDAhJ6pE1\na9aSZKTXuHiZqyT1SPMLftTfe17mKkmaIANCktTJgJAkdTIgJEmdDAhJUicDQpLUyYCQJHUyICRJ\nnQwISVInA0KS1GmiAZFkrySfS3JDko1J1rXL1yX5dpLr29cJk6xDkrR0Ex+LKcm+VXVnklXAZ4BX\nAE8FtlbVuTtY17GYJM2UmRqLqarubCf3Anbnnp98fEMOSpLGbuIBkWS3JDcAm4Erqura9q3Tk9yY\n5LwkqyddhyRpaVbiDOLuqno0cBjw2CQPB/4YeHBVHUMTHIs2NUmSVt7uK7WjqrojyQbghO36Hv4M\nuHih9davX79tejAYMBgMJlShJO2sNrSv8ZpoJ3WS+wE/q6rbk+wDXAa8Ebi+qja3n3kl8OtVdUrH\n+nZSS5opfeqknvQZxKHA+Ul2o2nO+mBV/XWSC5IcA9wN3AK8dMJ1SJKWyEeOSlKP9OkMwjupJUmd\nDAhJUicDQpLUyYCQJHUyICRJnQwISVInA0KS1MmAkCR1MiAkSZ0MCElSJwNCktTJgJAkdTIgJEmd\nDAhJUicDQpLUyYCQJHUyICRJnQwISVInA0KS1GmiAZFkrySfS3JDko1J1rXLD05yeZJNSS5LsnqS\ndUiSli5Voz4cewc7SPatqjuTrAI+A7wC+E/ArVX1piSvBg6uqjM71q1J1ydJfZIEGPX3XqiqjFrL\nxJuYqurOdnIvYHean/wk4Px2+fnAMyddhyRpaSYeEEl2S3IDsBm4oqquBQ6pqi0AVbUZeMCk65Ak\nLc3uk95BVd0NPDrJgcBfJnkE9z5/WvB8av369dumB4MBg8FgAlVK0s5sQ/sar4n3QfzSzpLXAXcC\n/wUYVNWWJGuAq6rqYR2ftw9C0kyZmT6IJPebu0IpyT7Ak4GbgYuAF7YfOxX4+CTrkCQt3UTPIJI8\nkqYTerf29cGq+qMk9wE+BBwOfBN4TlX9sGN9zyAkzZQ+nUGsaBPTUhkQkmZNnwLCO6klSZ0MCElS\nJwNCktTJgJAkdTIgJEmdDAhJUicDQpLUyYCQJHUyICRJnQwISVKnJQVE+2yHAydVjCSpP3YYEEku\nTHJgkv2ALwFfSfKqyZcmSZqmYc4gHl5Vd9A8FvQS4EjgBROtSpI0dcMExB5J9qAJiIuq6meMPtSg\nJKnnhgmIPwFuAfYDrk5yBHDHJIuSJE3fkp8HkWaw8lVV9fPJlPRL+/J5EJJmyk71PIgk/zPJQfMW\nHQSsH3XHkqR+G6aJ6anzHwdaVbcBJ06uJElSHwwTEKuS7DU3k2QfYK9FPi9J2gUMExDvBz6V5LQk\npwFXAOcPs/EkhyW5MsmXk2xM8vJ2+bok305yffs6Yfk/giRpEobqpE7yVOBJ7ewVVXXZUBtP1gBr\nqurGJPsDXwBOAn4X2FpV5+5gfTupJc2UPnVS7z7Mh6rqEpqb5JakqjYDm9vpHyW5GXhg+/bIxUuS\nJmfBJqYk17T/3ZrkjnmvrUmWfB9EkrXAMcDn2kWnJ7kxyXlJVi+jdknSBC35Pohl7aRpXtoAvL6q\nPp7k/sD3q6qSvAE4tKpO61iv1q1bt21+MBgwGAwmXq8kTcvympg2tK85/30sTUzD9kEcCzyepupr\nquqGoXeQ7A58Arikqt7W8f4RwMVVdXTHe/ZBSJopfeqDGOZGuT+kuWrpvsD9gPcmee0S9vF/gK/M\nD4e283rOs2hGiZUk9cgOzyCSbAIeVVU/aef3AW6sqofucOPJccDVwEaaSCzgLOAUmv6Iu2nGeXpp\nVW3pWN8zCEkzpU9nEMNcxfRdYG/gJ+38XsB3htl4VX0GWNXx1qVDVSdJmpphAuJ24MtJrqCJtScD\nn0/ydoCqesUE65MkTckwTUynLvZ+VQ11V/Vy2MQkadb0qYlp2KuY9gQe0s5uah8aNHEGhKRZ06eA\n2GETU5IBzVVMt9Dc/Xx4klOr6upRdy5J6q9hmpi+AJxSVZva+YcAH6iqX5t4cZ5BSJoxfTqDGOqZ\n1HPhAFBVXwP2GHXHkqR+G+YqpuuSnAf8eTv/POC6yZUkSeqDYZqY9gJeRjPUBsDfAO+qqrsmXJtN\nTJJmTp+amIYJiDO2H0Opa9kkGBCSZk2fAmKYPoiu+yBeOOqOJUn9tmAfRJKTacZMOjLJRfPeOgD4\nwaQLkyRN12Kd1H8LfI9mBNe3zFu+FbhpkkVJkqZvRR4YtFz2QUiaNTtbH4QkaQYZEJKkTgsGRJJP\ntf89Z+XKkST1xWKd1Icm+U3gGUn+gmagvm2q6vqJViZJmqoFO6mT/A5wGs0d1NsPrVFVdfyEa7OT\nWtLM6VMn9TB3Ur+uql4/6o6Ww4CQNGt2qoAASPIM4Ant7Iaq+sRQG08OAy4ADgHuBv6sqt6e5GDg\ng8ARNM+ZeE5V3d6xvgEhaabsVAGR5GzgscD720UnA9dW1Vk73HiyBlhTVTcm2R/4AnAS8CLg1qp6\nU5JXAwdX1Zkd6xsQkmbKzhYQNwHHVNXd7fwq4IaqOnrJO0s+BryzfT2xqra0IbKhqo7q+LwBIWmm\n9Ckghr0P4qB506uXs6Mka4FjgM8Ch1TVFoCq2gw8YDnblCRNzjAPDDobuCHJVTSXuj4BuFdz0GLa\n5qWPAGdU1Y+SbB+PC8bl+vXrt00PBgMGg8FSdi1JM2BD+xqvYTupDwV+vZ39fPtX/3A7SHYHPgFc\nMvcMiSQ3A4N5TUxXVdXDOta1iUnSTOlTE9PEB+tLcgHw/ar6g3nLzgF+UFXn2EktSfeYmYBIchxw\nNbCR5icu4Czg88CHgMOBb9Jc5vrDjvUNCEkzZWYCYlQGhKRZ06eAWPQqpiSrknx11J1IknY+iwZE\nVf0C2JTkQStUjySpJ4a5zPVg4MtJPg/8eG5hVT1jYlVJkqZumIB43cSrkCT1zrD3QRwB/GpVfTLJ\nvsCqqto68eLspJY0Y3aaTmqAJC+muQv6T9pFDwQ+NuqOJUn9NsxYTC8DjgPuAKiqr+PYSZK0yxsm\nIH5aVXfNzbRDZ9juI0m7uGEC4tNJzgL2SfJk4MPAxZMtS5I0bcM8D2I3mmdT/xbNaK6XAeetRO+x\nndSSZk2fOqmHvYppT+Aomqo3zW9ymiQDQtKs6VNA7PA+iCRPA/438A80ZxBHJnlpVV0y6s4lSf01\nTBPTV4GnV9Xft/P/CvirrkeEjr04zyAkzZg+nUEM00m9dS4cWt8AJn6TnCRpuhZsYkryrHbyuiR/\nTfP8hgKeDVy7ArVJkqZosT6I3543vQV4Yjv9z8A+E6tIktQLPjBIknqkT30Qw1zFdCTwcmDt/M87\n3Lck7dqGGe77Y8C7ae6evnuy5UiS+mKYgPhJVb19ORtP8m7g6cCWqjq6XbYOeDHwT+3HzqqqS5ez\nfUnS5AxzH8QpwK8ClwM/nVteVdfvcOPJ44EfARdsFxBbq+rcIda3D0LSTNmp+iCARwIvAI7nniam\naucXVVXXtA8b2t7IhUuSJmuYgHg28OAxj790epIXANcB/62qbh/jtiVJYzBMQHwJOIh7+gxG9cfA\n/6iqSvIG4Fya0WI7rV+/ftv0YDBgMBiMqQxJ2lVsaF/jNUwfxAbgaJq7p+f3QQx1mWvbxHTxXB/E\nsO+179sHIWmm7Gx9EOtG3EeY1+eQZE1VbW5nn0VzhiJJ6pmJ3kmd5EJgANyXZriOdcC/B46h6fC+\nBXhpVW1ZYH3PICTNlD6dQQzTxLSVe6rdE9gD+HFVHTjqznfEgJA0a/oUEDtsYqqqA7btsqn8JOBx\no+5YktRvy2piSnJDVT16AvVsvx/PICTNlJ3qDGLecyGgecDQY4CfjLpjSVK/DXMV0/znQvycpmP5\npIlUI0nqDZ8HIUk9slM0MSX5w0XWq6p6/ag7lyT112JNTD/uWLYfzbAY9wUMCEnahQ3VxJTkAOAM\nmnD4EPCWqhrX2EyL7dcmJkkzZadoYgJIch/gD4DnAecDx1bVbaPuVJLUf4v1QbyZZqykPwUeWVU/\nWrGqJElTt2ATU5K7aUZv/Tm/fL4Tmk5qh9qQpDHbKZqYqmq3UTcuSdp5GQKSpE4GhCSpkwEhSepk\nQEiSOhkQkqROBoQkqZMBIUnqNNGASPLuJFuS3DRv2cFJLk+yKcllSVZPsgZJ0vJM+gziPcBTtlt2\nJvDJqnoocCXwmgnXIElahokGRFVdA2w/uN9JNAP/0f73mZOsQZK0PNPog3hAVW0BqKrNwAOmUIMk\naQeGeSb1pC06KtX69eu3TQ8GAwaDwYTLkaSdzYb2NV4TfyZ1kiOAi6vq6Hb+ZmBQVVuSrAGuqqqH\nLbCuo7lKmil9Gs11JZqY0r7mXAS8sJ0+Ffj4CtQgSVqiiZ5BJLkQGNA8w3oLsA74GPBh4HDgm8Bz\nquqHC6zvGYSkmdKnM4iJNzGNwoCQNGv6FBDeSS1J6mRASJI6GRCSpE4GhCSpkwEhSepkQEiSOhkQ\nkqROBoQkqZMBIUnqZEBIkjoZEJKkTgaEJKmTASFJ6mRASJI6GRCSpE4GhCSpkwEhSepkQEiSOhkQ\nkqROu09rx0luAW4H7gZ+VlWPnVYtkqR7m1pA0ATDoKpum2INkqQFTLOJKVPevyRpEdP8BV3AFUmu\nTfLiKdYhSeowzSam46rqe0nuTxMUN1fVNdt/aP369dumB4MBg8Fg5SqUpJ3ChvY1XqmqsW90yUUk\n64CtVXXudsurD/VJ0kpJQtPAMtJWqKqMWstUmpiS7Jtk/3Z6P+C3gC9NoxZJUrdpNTEdAvxlkmpr\neH9VXT6lWiRJHXrRxLQQm5gkzZqZb2KSpHFbs2YtSUZ6rVq130jrr1mzdtpfw1h5BiFplzCuv7xH\n20YY9XeWZxCSpN4zICRJnQwISVInA0KS1MmAkCR1MiAkSZ0MCElSJwNCktTJgJAkdTIgJEmdDAhp\nGcYx7s+uNm6Pdj2OxSQtwxjHyxlHOcKxmDrqcCwmSdJkGBCSpE7TeqLc0C699NKR1j/++OPZc889\nx1SNJM2O3vdBrF79lGWv/9OffoW3ve21vOQlLxljVdOxZs1atmz55kjbOOSQI9i8+ZbxFDTj7IPo\nH/sg7lXHyH0QUzuDSHIC8L9omrneXVXndH3u9tuXfwax116nc9dddy17/T5pwmG0g2bLlpGPF0kz\nZCp9EEl2A94JPAV4BHBykqOmUctSbdiwYdoldNgw7QI69fG7sqbhWNNSbJh2ARMzrU7qxwJfr6pv\nVtXPgL8ATppSLUvSz4N0w7QL6NTH78qahmNNS7Fh2gVMzLQC4oHAt+bNf7tdJknqid5fxXTggb+9\n7HXvumsje+zxyDFWI0mzYypXMSV5HLC+qk5o588EavuO6iRe4iFJyzCOq5imFRCrgE3Ak4DvAZ8H\nTq6qm1e8GElSp6k0MVXVL5KcDlzOPZe5Gg6S1CO9vlFOkjQ907oP4oQkX03ytSSv7nj/wCQXJbkx\nycYkL2yXH5bkyiRfbpe/og91zXt/tyTXJ7moDzUlWZ3kw0lubr+z3+hBTa9M8qUkNyV5f5KxjIMy\nRE0HJfloki8m+WyShw+77krX1IPjfMHvqn1/Gsf5Yv9+0zrOF6tpUsf5u5NsSXLTIp95e5Kvt///\nHTPsz9Opqlb0RRNKfw8cAewB3Agctd1nXgOc3U7fD7iVpjlsDXBMu3x/mn6Mo6Zd17z3Xwn8OXBR\nH2oC3gu8qJ3eHThwyv9+vwJ8A9izfe+DwO+tUE1vAl7XTj8U+OSw606hpmkf5511Tfk4X7CmKR7n\nC/37TeQ4b7f1eOAY4KYF3n8q8Fft9G8Anx3lOJ/GGcQwN8kVcEA7fQBwa1X9vKo2V9WNAFX1I+Bm\nxnf/xLLrguavPuBE4Lwx1TNSTUkOBP5dVb0HoP3+7phmTe38KmC/JLsD+wLfXaGaHg5cCVBVm4C1\nSe4/5LorWlMPjvOFvqtpHuedNU35OF/we2IyxzlVdQ1w2yIfOQm4oP3s54DVSQ4Z8ue5l2kExDA3\nyb0TeHiS7wJfBM7YfiNJ1tIk6ed6UtdbgVcx+ihb46rpSOD7Sd7TNgf8aZJ9pllTVX0XeAvw/4Dv\nAD+sqk+uUE1fBJ4FkOSxwIOAw4Zcd6Vr2mZKx/lidU3rOF+opmke5501TfA4H8ZCdS/rOO/r8yCe\nAtxQVb8CPBp4V5L9595spz8CnNH+hTXVupI8DdjS/tWX9jXVmmhOtY8F3lVVxwJ3AmdOs6YkB9H8\n1XIEzWn4/klOWaGa3ggcnOR64GXADcAvVmjfC1m0pike5511Tfk4X+i7muZxvtD3NM3jfHsj/RtN\n4zLX79Ak7ZzD2mXzvQg4G6Cq/iHJPwJHAde1p2wfAd5XVR/vSV3HAc9IciKwD3BAkguq6vemWNO3\ngG9V1XXt5z4CjKMDdpSa1gLfqKofACT5KPCbwIWTrqmqtgL/eW6+rekbNKf/O/p5Vrompnmcd9T1\njbau5zKl43yR72o/pnScL/I9ncBkjvNhfAc4fN78XN17spzjfBwdJ0vsZFnFPZ0le9J0ljxsu8+8\nC1jXTh9C88vuPu38BcC5fatr3meeyPg670b9rj4NPKSdXgecM82aaNpBNwJ70/xl817gZStU02pg\nj3b6xcB7h113pWvqwXG+YF1TPM4X+66mdZwvdExN5Dift9+1wMYF3juRezqpH8c9ndTLOs7HevAt\n4Qc8gebKjK8DZ7bLXgq8pJ0+FLgMuKl9ndwuP47mtPJGmtO564ETpl3XdtsY2/84o9YEPAq4tv2+\nPgqs7kFN62g6XW8Czp/7H2wFanpc+/7NNH9lrl5s3WnW1IPjfMHvaorH+WL/ftM6zheraVLH+YU0\nHd4/penjeNH8mtrPvJMmDL4IHDvKce6NcpKkTn3tpJYkTZkBIUnqZEBIkjoZEJKkTgaEJKmTASFJ\n6mRAaJeTZOuYtnNJktu2H9Y6ydp2eOevJflAe9fzjrZ1VZJj2+lPtIPMLfTZM5LsPfpPII3GgNCu\naFw397wJeH7H8nOAt1TVQ4AfAqctZaNV9fRafMTR/0ozBIg0VQaEZkKSI5J8qn2IyhXtsNUkeXCS\nv2sf+vL6+WcfVXUV0DVI3vHA/22nzwf+Y8f+9m7PLr7cjsWz97z3/jHJfZLs255N3NA+WObZSV5O\nM8DbVUk+Nb5vQFo6A0Kz4h3Ae6rqGJrhCt7RLn8b8NaqehTNEMiLnn0kuS9wW1Xd3S76Ns0v9O39\nPvDjqnoEzbALj5n33tw+TgC+U1WPrqqjgUur6h00g6gNqupJS/0hpXEyIDQr/i3wgXb6fTTjHc0t\n/0g7Pc7RNp9A89Q1qmojzbg4c+aGYN4IPDnJ2UkeX83ooHPvr+RQ2lInA0KzYph+iR3+Uq6qW4GD\nksz9vzPs8OD32nZVfZ3mWQYbgTckee0Q25FWjAGhXVHXL/q/BU5up58P/E07/XfA77TTz11gW9tv\n7yrg2e30qUDX8xquBp4HkOTfAEffa8PJocC/VNWFwJtpwgLgDmDBq5ykleJortrlJPk5zZDIoTlz\nOJemU/m9wH2Bf6Z5yP23k/xrmqagvWmGKD+lqg5vt3M1zcPo9wduBU6rqiuSHEnzTN+DaYbjfn41\nz/mdX8PewHtogmHumdIvq6rr2wfLPKZ9vRm4G7gL+P32/dOB02n6J+yH0NQYEJppSfapqn9pp38X\neG5V3euqJGkWTeORo1Kf/FqSd9KcbdzGvEdISrPOMwhJUic7qSVJnQwISVInA0KS1MmAkCR1MiAk\nSZ0MCElSp/8P2bBpC/baz6MAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fcf0ce51ad0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "50\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEPCAYAAABCyrPIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGJZJREFUeJzt3XmwbWV55/Hvj0EmkcmWq6JcsENQIyIOcaB1F7YtTmiM\nJoqzlJWylZDutO3QMfdWmXJK2c7diS1BbRstozZBjY1o2DGUAyIgKIMjqCjXiChohID36T/WOnA4\nnGHtc87aa597v5+qXay99nrX++xd5+6H9332eleqCkmSdhk6AEnSbDAhSJIAE4IkqWVCkCQBJgRJ\nUsuEIEkCek4ISU5Nsi3JxQv2n5zksiSXJHljnzFIkrrZrefznwa8E/jA3I4kI+ApwAOq6pYkd+05\nBklSB72OEKrqXOC6BbtfCryxqm5pj/lpnzFIkroZooZwBPDoJF9Kck6ShwwQgyRpgb6njJbq84Cq\neniShwIfAQ4fIA5J0jxDJIQfAB8HqKqvJNme5KCqunbhgUlcaEmSVqGqMmmbaUwZpX3MOQM4DiDJ\nEcDuiyWDOVW1YR9btmwZPIadNf6NHLvxD//YsmXL3DfQKh7Df2+tVq8jhCSnAyPgoCTfB7YAfwOc\nluQS4Cbg+X3GIEnqpteEUFUnLvHS8/rsV5I0Oa9U7tFoNBo6hDXZyPFv5NjB+Ie20eNfraxlvqlv\nSWqW45O040rCXE1gwpZrmsdfD0moGS0qS5I2ABOCJAkwIUiSWiYESRJgQpAktUwIkiTAhCBJapkQ\nJEmACUGS1DIhSJIAE4IkqWVCkCQBJgRJUsuEIEkCTAiSpJYJQZIEmBAkSa1eE0KSU5NsS3LxIq/9\naZLtSQ7sMwZJUjd9jxBOAx6/cGeSQ4DHAVf13L8kqaNeE0JVnQtct8hLbwVe0WffkqTJTL2GkOQE\n4AdVdcm0+5YkLW23aXaWZC/gNTTTRbfunmYMkqTFTTUhAPcBNgNfSxLgEOCrSR5WVT9ZrMHWrVtv\n3R6NRoxGo/6jlLSkTZs2s23b6sp/Bx98KNdcc+X6BiTG4zHj8XjN50lVrT2a5TpINgOfqKoHLPLa\n94BjqmqxOgNJqu/4JE2m+X+51f67DBvl3/Tq3+fw7zEJVTXx7EvfPzs9HfgCcESS7yd50YJDCqeM\nJGkm9D5CWAtHCNLscYSwYsvB3+NMjhAkSRuHCUGSBJgQJEktE4IkCTAhSJJaJgRJEmBCkCS1TAiS\nJMCEIElqmRAkSYAJQZLUMiFIkgATgiSpZUKQJAEmBElSy4QgSQJMCJKklglBkgSYECRJrV4TQpJT\nk2xLcvG8fW9OclmSi5J8LMld+oxBktRN3yOE04DHL9j3GeD+VXU08C3g1T3HIEnqoNeEUFXnAtct\n2PfZqtrePv0ScEifMUiSuhm6hvBi4NMDxyBJAnYbquMk/w24uapOX+64rVu33ro9Go0YjUb9BiZJ\nG8x4PGY8Hq/5PKmqtUezXAfJocAnquqoefteCLwEOK6qblqmbfUdn6TJJAFW++8ybJR/06t/n8O/\nxyRUVSZtN40RQtpH8yQ5HngF8OjlkoEkabp6HSEkOR0YAQcB24AtwGuAOwHXtod9qar+4xLtHSFI\nM8YRwootB3+Pqx0h9D5ltBYmBGn2mBBWbDn4e1xtQhj6V0aSpBlhQpAkASYESVLLhCBJAkwIkqSW\nCUGSBJgQJEktE4IkCTAhSJJaJgRJWld7kGTix6ZNm4cOfLKlK5LsAty5qq7vL6Tb9efSFdKMcemK\nFVuuut16fTa9LV2R5PQkd0myD/B14NIkr1hNkJKk2dVlyuh+7YjgaTR3NzsMeF6vUUmSpq5LQtg9\nye40CeHMqrqZ1Y8XJUkzqktC+GvgSmAf4PPtHdCmUkOQJE3PxPdDSFNp2bWqbuknpNv1ZVFZmjEW\nlVdsuep2G6Go/Pok+8/btT+wddKOJEmzrcuU0ROq6udzT6rqOuCJ/YUkSRpCl4Swa5I95p4k2QvY\nY5njJUkbUJeE8H+AzyU5KclJwNnA+7ucPMmpSbYluXjevgOSfCbJFUnOSrLf6kKXJK2nTkXlJE8A\nHts+Pbuqzup08uRY4JfAB6rqqHbfm4Brq+rNSV4JHFBVr1qivUVlacZYVF6x5arbDV1UnvhXRhN3\n0PxM9RPzEsLlwGOqaluSTcC4qo5coq0JQZoxJoQVW6663dAJYckpoyTntv+9Icn18x43JFnLdQh3\nq6ptAFV1DXC3NZxLkrROdlvqhao6tv3vvj3HsGxK3Lp1663bo9GI0WjUcziSdhSbNm1m27arhg6j\nd+PxmPF4vObzdK0hHAMcS/PlfW5VXdi5gztOGV0GjOZNGZ1TVfddoq1TRtKM2UhTRmuN1SmjO574\nz2l+VXQQcFfgfUn+bJLY2secM4EXttsvAP5ugnNJknqy4gghyRXAA6vqxvb5XsBFVfXbK548OR0Y\n0SSTbcAW4Azgb4F7AVcBfzD/wrcF7R0hSDPGEUJ/7YYeISxZQ5jnR8CewI3t8z2Aq7ucvKpOXOKl\nf9+lvSRperokhF8A30hyNk3aexxwXpJ3AFTVH/cYnyRpSrokhP/bPuaM+wlFkjSkrr8yuhNwRPv0\nivYmOb2zhiDNHmsI/bWb+RpCkhHNr4yupHmn90rygqr6/KSdSZJmV5dfGX0VOLGqrmifHwF8qKoe\n3HtwjhCkmeMIob92Q48QOt1TeS4ZAFTVN4HdJ+1IkjTbuhSVz0/yXuCD7fPnAOf3F5IkaQhdpoz2\nAF5Gs3QFwD8B766qf+05NqeMpBnklFF/7YaeMuqSEE6pqrevtK8PJgRp9pgQ+ms3dELoUkN4wSL7\nXjhpR5Kk2bZkDSHJs4ETgcOSnDnvpX2Bn/UdmCRpupYrKn8B+DHNCqdvmbf/BuDiRVtIkjas3m+h\nuRbWEKTZYw2hv3YboYYgSdoJmBAkScAyCSHJ59r/vml64UiShrJcUfnuSR4JnJDkw9z+NphU1QW9\nRiZJmqoli8pJngGcRHOF8sKlKqqqjus5NovK0gyyqNxfu6GLyl2uVH5tVb1u1ZGtgQlBmj0mhP7a\nzXxCaE9+AvDo9um4qj45aUeLnPM/0YxAtgOXAC9auD6SCUGaPSaE/toNnRBW/JVRkjcApwCXto9T\nkrx+8hBvd857ACcDx1TVUTS1jGet5ZySpLXpsvz1k4Cjq2o7QJL3AxcCr1lj37sC+yTZDuwN/GiN\n55MkrUHX6xD2n7e931o7raof0SyH8X3gauDnVfXZtZ5XkrR6XUYIbwAuTHIOzeTYo4FXraXTJPsD\nTwUOBX4BfDTJiVV1+sJjt27deuv2aDRiNBqtpWtJg9qjndefzMEHH8o111y5/uHsIMbjMePxeM3n\n6VpUvjvw0PbpeVV1zZo6bX7S+viqekn7/HnA71bVyxccZ1FZmjFDFWpX811gUXkyXUYIVNWPgTNX\nPLC77wMPT7IncBPwWOAr63h+SdKEBlnLqKrOAz5KU5z+Gk1Kfc8QsUiSGi5/LWkiThn1127oKaNl\nRwhJdk1y+erDkiRtFMsmhKr6DXBFkntPKR5J0kC6FJUPAL6R5DzgV3M7q+qE3qKSJE1dl4Tw2t6j\nkCQNrut1CIcCv1VVn02yN7BrVd3Qe3AWlaWZY1G5v3YzXVRuT/wSmp+I/nW7657AGZN2JEmabV2u\nQ3gZ8CjgeoCq+hZwtz6DkiRNX5eEcNP8+xQk2Y3Vj8EkSTOqS0L4xySvAfZK8jjgb4FP9BuWJGna\nutxCcxeaO5v9B5pqyVnAe6dR7bWorKFs2rSZbduumrjdLrvszfbt/zJxu420mqdF5f7aDV1U7vor\nozsBR9K8yysW3uqyLyYEDWX1XyTDfxn0zYTQX7uhE8KK1yEkeRLwV8B3aN7pYUn+qKo+PXmYkqRZ\n1WXK6HLgyVX17fb5fYBPVdWRvQfnCEEDcYSwNEcI/bUbeoTQpah8w1wyaH0X6P2iNEnSdC05ZZTk\n6e3m+Un+HvgITdp7Jt7MRpJ2OMvVEJ4yb3sb8Jh2+5+BvXqLSJI0CG+QIy3CGsLSrCH0127oGkKX\nXxkdBpwMbJ5/vMtfS9KOpcvy12cAp9Jcnby933AkSUPpkhBurKp3rHfHSfYD3gv8Dk2ieXFVfXm9\n+5EkddPlOoQTgd8CPgPcNLe/qi5YU8fJ+4B/rKrT2gXz9q6q6xccYw1Bg7CGsDRrCP21m/kaAvAA\n4HnAcdw2ZVTt81VJchfg31XVCwGq6hba5bUlScPokhCeCRy+zusXHQb8NMlpwAOB84FTqurX69iH\nJGkCXRLC14H9gZ+sc7/HAC+rqvOTvA14FbBl4YFbt269dXs0GjEajdYxDGnjW+3KrBtphVXYo53+\n0WLG4zHj8XjN5+lSQxgDR9FcnTy/hrDqn50mORj4YlUd3j4/FnhlVT1lwXHWEDSIjVRDWEusO8O8\n/EaKdSPUEO7wf+1rVVXbkvwgyRFV9U3gscCl692PJKm7wa5UTvJAmp+d7k6zYN6LquoXC45xhKBB\nOELoo7+mT0cIS7cbeoTQZcroBm57d3ei+QL/VVXdZeIoJ2RC0FBMCH301/RpQli63dAJYcUpo6ra\nd14nAZ4KPHzSjiRJs21VU0ZJLqyqB/UQz8J+HCFoEI4Q+uiv6dMRwtLtZn6EMO++CNDcUOchwI2T\ndiRJmm1dfmU0/6egtwBX0kwbSZJ2IN4PQVqEU0Z99Nf06ZTR0u1mdsooyZ8v066q6nWTdiZJml3L\nTRn9apF9+wAnAQcBJgRJ2oF0mjJKsi9wCk0y+Ajwlqpaz7WNlurXKSMNwimjPvpr+nTKaOl2Mztl\n1J70QOA/A88B3g8cU1XXrS5ESdIsW66G8JfA04H3AA+oql9OLSpJU+AKorq9JaeMkmynWd30Fm4/\n/glNUdmlK7TD2lmmjHaGaZiNFOvMThlV1S5rC0mStJH4pS9JAkwIkqSWCUGSBJgQJEktE4IkCTAh\nSJJaJgRJEjBwQkiyS5ILkpw5ZBySpOFHCKcAlw4cgySJARNCkkOAJwLvHSoGSdJthhwhvBV4Batf\naESStI663FN53SV5ErCtqi5KMqJZDWpRW7duvXV7NBoxGo36Dk+SNpTxeMx4PF7zeQa5p3KS1wPP\npVlJdS9gX+DjVfX8Bce52qkG4WqnfbQbos+NFevQq50OkhBuF0DyGOBPq+qERV4zIWgQJoQ+2g3R\n58aKdeiEMPSvjCRJM2LwEcJyHCFoKI4Q+mg3RJ8bK1ZHCJKkmWBCkCQBJgRJUsuEIEkCTAiSpJYJ\nQZIEmBAkSS0TgiQJMCFIklomBEkSYEKQJLVMCJIkwIQgSWqZECRJgAlBktQyIUiSABOCJKllQpAk\nASYESVJrkISQ5JAk/5DkG0kuSfLHQ8QhSbpNhriJfZJNwKaquijJnYGvAk+tqssXHFdDxCcNceP6\n1f6tDxHrRrpx/UaKdb2+75JQVZm03SAjhKq6pqouard/CVwG3HOIWCRJjcFrCEk2A0cDXx42Ekna\nue02ZOftdNFHgVPakcIdbN269dbt0WjEaDSaSmyaHZs2bWbbtqsmbrfLLnuzffu/9BBRH/Zop36k\nyY3HY8bj8ZrPM0gNASDJbsAngU9X1duXOMYaggaYI19LW2OdrT43Vqw7ZQ2h9TfApUslA0nSdA31\ns9NHAc8BjktyYZILkhw/RCySpMZgU0ZdOGUkcMpo9vo01r7a7cxTRpKkGWJCkCQBJgRJUsuEIEkC\nTAiSpJYJQZIEmBAkSS0TgiQJGHhxu1n0rnf9L84665xVtX31q0/mkY98xDpHJEnTYUJY4G1v+yu+\n853fAw6fsOWZHHHEGVNPCNNeCfTggw/lmmuunLgdrD5WSdNhQljUE4AHT9jmB8DPeohlec0X7OSX\nu2/fvrrL67dtW/0SzauNtVkKQFLfrCFIkgATgiSpZUKQJAEmBElSy4QgSQJMCJKklglBkgQMmBCS\nHJ/k8iTfTPLKoeKQJDUGSQhJdgHeBTweuD/w7CRHDhFLn8bj8dAhrNF46ADWYDx0AGs0HjqANRoP\nHcAajYcOYBBDjRAeBnyrqq6qqpuBDwNPHSiW3pgQhjQeOoA1Gg8dwBqNhw5gjcZDBzCIoRLCPWnW\nepjzw3afJGkgrmW0wB577M4++/wJu+66/0Ttbrrp2+y+++/3FJUk9S9Vq1lsbI2dJg8HtlbV8e3z\nVwFVVW9acNz0g5OkHUBVTbwq5FAJYVfgCuCxwI+B84BnV9VlUw9GkgQMNGVUVb9J8nLgMzR1jFNN\nBpI0rEFGCJKk2TMTVyqvdJFakhOSfC3JhUnOS/KoIeJcTNcL7JI8NMnNSZ4+zfhW0uGzf0ySnye5\noH382RBxLqXL559k1P7tfD3J6u6P2pMOn/9/aWO/IMklSW5JMtkvHnrUIf67JDkzyUVt/C8cIMwl\ndYh//yQfb79/vpTkfkPEuZgkpybZluTiZY55R5JvtZ//0SuetKoGfdAkpW8DhwK7AxcBRy44Zu95\n2w8ALhs67q6xzzvuc8AngacPHfeEn/1jgDOHjnUN8e8HfAO4Z/v8rkPHPenfz7zjnwx8dui4J/z8\nXw28Ye6zB64Fdhs69gnifzPw2nb7t2fs8z8WOBq4eInXnwB8qt3+XeBLK51zFkYIK16kVlXzb/57\nZ2D7FONbTtcL7E4GPgr8ZJrBddA1/lm9h2WX+E8EPlZVVwNU1U+nHONyJr1A89nAh6YSWTdd4i9g\n33Z7X+DaqrplijEup0v89wP+AaCqrgA2J/k30w1zcVV1LnDdMoc8FfhAe+yXgf2SHLzcOWchIXS6\nSC3J05JcBnwCePGUYlvJirEnuQfwtKr6n8zeF2vXCwQf0Q45PzVLQ2a6xX8EcGCSc5J8Jcnzphbd\nyjpfoJlkL+B44GNTiKurLvG/C7hfkh8BXwNOmVJsXXSJ/2vA0wGSPAy4N3DIVKJbu4Xv72pWuAB4\nFhJCJ1V1RlXdF3ga8BdDxzOBtwHz5yZnLSms5KvAvavqaJp/3GcMHM+kdgOOoRk+Hw+8Nsm/HTak\nVXkKcG5V/XzoQCb0eODCqroH8CDg3UnuPHBMk3gjcECSC4CXARcCvxk2pP7MwpXKV9Nk3TmHtPsW\nVVXnJjk8yYFV9bPeo1tel9gfAnw4SWjmUJ+Q5OaqOnNKMS5nxfir6pfztj+d5H/MyGcP3T7/HwI/\nraobgRuTfB54IM3c8dAm+dt/FrM1XQTd4n8R8AaAqvpOku8BRwLnTyXC5XX5+7+BeTMSbfzfnUp0\na3c1cK95z5f9bgVmoqi8K7cVdu5EU9i574Jj7jNv+xjgB0PH3TX2BcefxmwVlbt89gfP234YcOXQ\ncU8Y/5HA2e2xewOXAPcbOvZJ/n5oCuPXAnsNHfMqPv93A1vm/pZopjAOHDr2CeLfD9i93X4J8L6h\n414Q32bgkiVeeyK3FZUfToei8uAjhFriIrUkf9S8XO8Bfj/J84F/BX4N/MFwEd+mY+y3azL1IJfR\nMf5nJHkpcDPNZ/+Hw0V8e13ir6rLk5wFXEwz1H9PVV06YNi3muDv52nAWVX166FiXUzH+P8CeN+8\nn0b+15qN0WXX+O8LvD/Jdppfq500XMS3l+R0YAQclOT7wBaaxDb3t//3SZ6Y5NvAr2hGa8ufs80e\nkqSd3IYpKkuS+mVCkCQBJgRJUsuEIEkCTAiSpJYJQZIEmBC0A0pywzqd59NJrkty5oL9m9ulkL+Z\n5ENJVryep11L6Zh2+5NJ7rLMsack2XPt70CajAlBO6L1urjmzcBzF9n/JuAtVXUE8HMmvFipqp5c\nVdcvc8if0FxVLU2VCUE7hSSHJvlcu2rr2UkOafcfnuSL7Q1QXjd/dFFV5wC/XOR0x3HbqqPvB35v\nkf72bEcP30jycWDPea99L8mBSfZuRwsXJrk4yTOTnAzcAzgnyefW7xOQVmZC0M7incBp1azaenr7\nHODtwFur6oE0C+EtO7pIchBwXVXN3ZPjhzRf4Au9FPhVVd2fZkmBh8x7ba6P44Grq+pBVXUU8P+q\n6p00C5CNquqxk75JaS1MCNpZPILbVgv938Cj5u3/aLt9+jr292jggwBVdQnNuvpz5pZAvwR4XJI3\nJDm2mpU1517faMukawdgQtDOoktdYcUv4aq6Ftg/ydy/nZWXFF7i3FX1LZrVey8B/iIzdr9q7XxM\nCNoRLfbF/gWaW1BCUyj+p3b7i8Az2u1nLXGuhec7B3hmu/0C4O8Wafd54DkASX4HOOoOJ07uDvy6\nqk4H/pImOQBcDyz5KySpL652qh1OkluAH9F8kRfw32mKwO8DDgL+GXhRVf2wvXvaB2mKvmcBJ1bV\nvdrzfJ7mxup3prkfwUlVdXaSw2juv3sAzR20nlvNPXnnx7Anzf0vjgIuo7l14cuq6oIk36WpKTyE\nJhFsp1na/aXt6y8HXk5TX7COoKkxIWinlmSvufsMJPlD4FlVdYdfDUk7g8FvkCMN7MFJ3kUzmriO\nebdLlHY2jhAkSYBFZUlSy4QgSQJMCJKklglBkgSYECRJLROCJAmA/w/UqPdUjdfxDgAAAABJRU5E\nrkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fcf0d2e9110>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "97\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEPCAYAAABCyrPIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF3VJREFUeJzt3XuQbWV95vHvw8WDgHJLpFWUo5koGsUjUWMGoz1aiWiM\nMIxXSLyMZVmWGmamZkZixZxTpRWCKZ3RqJMYDTCj6CTGGMAgeKGHMIlXQC4CagwY0HOSIERwvIDn\nN3+s1dIc+7L2Xr332t3n+6naxdq797veX6+zez+s912XVBWSJO0zdAGSpNlgIEiSAANBktQyECRJ\ngIEgSWoZCJIkYMKBkOSoJJ9Ocm2Sq5O8rn19e5Kbk1zePk6YZB2SpLVlkuchJJkD5qrqyiQHA18E\nTgReCNxRVW+bWOeSpJHsN8mVV9VOYGe7fGeS64AHtz/OJPuWJI1manMISbYC24DPti+9NsmVSd6b\n5JBp1SFJWt5UAqEdLvowcFpV3Qm8G3h4VW2j2YNw6EiSBjbROQSAJPsBFwAXVtXbl/n50cD5VXXs\nMj/zQkuSNIaqGnlYfhp7CH8CfHlpGLSTzYtOBq5ZqXFVbdjH9u3bB69hb61/I9du/cM/Nnr945ro\npHKS44FTgauTXAEU8AbglCTbgN3AjcCrJlmHJGltkz7K6P8C+y7zo49Psl9J0ug8U3mC5ufnhy6h\nl41c/0auHax/aBu9/nFNfFK5jyQ1y/VJ0ixKQs3opLIkaQMwECRJgIEgSWoZCJIkwECQJLUMBEkS\nYCBIkloGgiQJMBAkSS0DQZIEGAiSpJaBIEkCDARJUstAkCQBBoIkqWUgSJIAA0GS1DIQJEmAgSBJ\nahkIkiTAQJDUw9zcVpKM/Zib2zr0r6AlUlVD17CiJDXL9Ul7uyRAn7/R4N/4+ktCVWXUdu4hSJIA\nA0GS1DIQJEmAgSBJahkIkiTAQJAktQwESRJgIEiSWgaCJAkwECRJLQNBkgQYCJKk1kQDIclRST6d\n5NokVyf5zfb1w5JcnOSGJBclOWSSdUiS1jbRq50mmQPmqurKJAcDXwROBF4O3FpVb0nyeuCwqjp9\nmfZe7VSaYV7tdDbN5NVOq2pnVV3ZLt8JXAccRRMK57RvOwc4aZJ1SJLWNrU5hCRbgW3AZ4Ajq2oX\nNKEBPGBadUiSlrffNDpph4s+DJxWVXcm2XMfccV9xh07dvx4eX5+nvn5+UmUKEkb1sLCAgsLC73X\nM/E7piXZD7gAuLCq3t6+dh0wX1W72nmGS6rqUcu0dQ5BmmHOIcymmZxDaP0J8OXFMGidB7ysXX4p\n8JdTqEOStIpJH2V0PHApcDXN/0YU8Abgc8CfAg8BbgJeUFW3L9PePQRphrmHMJvG3UOY+JBRHwaC\nNNsMhNk0y0NGkqQNwECQJAEGgiSpZSBIkgADQZLUMhAkSYCBIElqGQiSJMBAkCS1DARJEmAgSJJa\nBoIkCTAQJEktA0GSBBgIkqSWgSBJAgwESVLLQJAkAQaCJKllIEiSAANBktQaKRCS7JPk/pMqRpI0\nnDUDIcm5Se6f5CDgGuDLSf7L5EuTNA1zc1tJMtajvy1j9z03t3Ud+tdSqarV35BcWVXbkpwKHAec\nDnyxqo6deHFJrVWfpH6aL/Zx/876tO3bPvj9sLwkVNXIid1lyGj/JPsDJwHnVdVd9PsESJJmUJdA\n+CPgRuAg4NIkRwPfmWRRkqTpW3PI6CcaNPuX+1bV3ZMp6V59OWQkTZhDRpvPxIaMkvxukkOXvHQo\nsGPUjiRJs63LkNGzqur2xSdVdRvw7MmVJEkaQpdA2DfJlsUnSe4LbFnl/ZKkDWi/Du/5APCpJGe1\nz18OnDO5kiRJQ+g0qZzkWcAz2qefqKqLJlrVPf06qSxNmJPKm8+4k8ojH2U0TQaCNHkGwuYzbiCs\nOGSU5LKqekqSO7j3v1iAqiqvaSRJm4h7CNJezj2EzWeSl64gyXFJfjPJ65I8foSi3pdkV5Krlry2\nPcnNSS5vHyeMWrQkaf11OTHtd2iOKjoC+Cng7CS/3XH9ZwHPXOb1t1XVce3j452rlSRNTJfDTk8F\nHldV3wdI8nvAlcCb12pYVZe11z7a03pcN1eStI66DBl9EzhgyfMtwC09+31tkiuTvDfJIT3XJUla\nB10C4V+Aa5Oc3Z6cdg1we5J3JHnHGH2+G3h4VW0DdgJvG2MdkqR11mXI6C/ax6KFPh1W1T8tefrH\nwPmrvX/Hjh0/Xp6fn2d+fr5P99KK5ua2smvXTWO1PfLIo9m588b1LUjqaGFhgYWFhd7r6Xqm8n2A\nR7RPb2hvktOtg2QrcH5VPbZ9PldVO9vl/wg8sapOWaGth51qavoefrlRP6sedrr5rPuJaUtWPE9z\nlNGNNP96D0ny0qq6tEPbc4F54Igk3wC2A/8myTZgd7vOV41atCRp/XW5p/IXgVOq6ob2+SOAD1bV\nz0+8OPcQNEXuIYzVukfbvu037jaftIneU3kxDACq6ivA/qN2JEmabV0mlb+Q5L3A+9vnpwJfmFxJ\nkqQhdBky2gK8BnhK+9JfA++qqh9OuDaHjDRVDhmN1bpH277tN+42n7SJXf46yWlV9fa1XpsEA0HT\nZCCM1bpH277tN+42n7RJziG8dJnXXjZqR5Kk2bba/RBeDJwCPCzJeUt+dD/g25MuTJImpc9JiLB5\nT0RcbVL5b4Bv0Vzh9K1LXr8DuGrZFpK0ATRhMP5w065dm/P6nN4gR2o5hzBW6x5t+7Yff5v3+537\n9T0NE71BjiRp8zMQJEnAKoGQ5FPtf8+cXjmSpKGsNqn8wCT/Gnhukg+xx13OquryiVYmSZqqFSeV\nkzwPeAXNGcp7XqqiqurpE67NSWVNlZPKY7Xu0bZveyeVVzLJM5XfWFVvGruyHgwETZOBMFbrHm37\ntjcQVjKxQGhX/lzgqe3Thaq6YNSOxmEgaJoMhLFa92jbt72BsJKJHXaa5AzgNODL7eO0JL87eomS\npFnWZcjoKmBbVe1un+8LXFFVx068OPcQNEXuIYzVukfbvu3dQ1jJpE9MO3TJ8iGjdiJJmn1dbpBz\nBnBFkkto4vypwOkTrUqSNHVdJ5UfCDyxffq5qto50aru6dchI02NQ0Zjte7Rtm97h4xWMtGjjIZi\nIGiaDISxWvdo27e9gbASL24nSerFQJAkAWsEQpJ9k1w/rWIkaWPYQpKxHnNzW4cufkWrHmVUVT9K\nckOSh1bVN6ZVlCTNth8w7hzELN9trcthp4cB1yb5HPDdxRer6rkTq0qSNHVdAuGNE69CkjS4ruch\nHA38bFV9MsmBwL5VdcfEi/OwU02Rh52O1bpH277thz3sdJY/K5O8uN0rgQ8Df9S+9GDgo6N2JEma\nbV0OO30NcDzwHYCq+irwgEkWJUmavi6B8IOq+uHikyT70W9fS5I0g7oEwv9J8gbgvkl+Gfgz4PzJ\nliVJmrYu90PYh+beyr9CM5NyEfDeacz2OqmsaXJSeazWPdr2be+k8oo9TPgWmvcBjqHZAjcsHUKa\nJANB02QgjNW6R9u+7Q2EFXsYMxDWPA8hya8Cfwj8Hc1WeFiSV1XVhaOXKUmaVV2GjK4HnlNVX2uf\n/wzwsao6ZuLFuYegKXIPYazWPdr2be8ewoo9TPDy13cshkHr60Cnk9KSvC/Jrva+zIuvHZbk4vYa\nSRcl8ZackjQDVgyEJCcnORn4QpK/SvKyJC+lOcLo8x3XfxbwzD1eOx34ZFU9Evg08Ftj1C1JWmcr\nDhklOWu1hlX18k4dNJe9OL+qjm2fXw88rap2JZkDFlYafnLISNPkkNFYrXu07dveIaMVe1jvSeWu\nX/hjeEBV7Wr72JnEs54laQZ0OcroYcDrgK1L37+Ol7/emP9bJUmbTJfLX38UeB/N3MHudehzV5Ij\nlwwZ/eNqb96xY8ePl+fn55mfn1+HEqTNY25uK7t23TR0GRrQwsICCwsLvdfT5bDTz1bVL4zdQbKV\nZg7hse3zM4FvV9WZSV4PHFZVp6/Q1jkETc1GnUMYejx8I47FD73NZnUOoUsgnAL8LHAxzX3jAKiq\nyzsUdS4wDxwB7AK20+xx/BnwEOAm4AVVdfsK7Q0ETY2BMO22/fs2EFboYYKBcAbwGzRnKi8OGVVV\nPX3kKkdkIGiaDIRpt+3ft4GwQg+TunQF8Hzg4dO6fpEkaRhdzlS+Bjh00oVIkobVZQ/hUOD6JJ/n\n3nMI63XYqSRpBnQJhO0Tr0KSNLhO90MYipPKmqZ+E40HsGQHemRHHnk0O3feOFbboSdIh+u73zZ3\nUnmZdh2OMrqDe37z+wD7A9+tqvuPXOWIDARN09DX9NmoR8xszL49ymg5aw4ZVdX9lnQS4ETgyaN2\nJEmabWMNGSW5oqoeP4F69uzHPQRNjXsI0247ZN/uISyny8XtTl7ydB/gCcD3R+1IkjTbuhxl9GtL\nlu8GbqQZNpIkbSIeZSS1HDKadtsh+3bIaDkr7iEk+Z1V2lVVvWnUziRJs2u1IaPvLvPaQcAraK5e\naiBI0ibSacgoyf2A02jC4E+Bt1bVqje2WQ8OGWmaHDKadtsh+3bIaDmrTionORz4T8CpwDnAcVV1\n23glSpJm2WpzCL8PnAy8B3hsVd05taokSVO34pBRkt00Fwq5m3vvG4VmUtlLV2hTccho2m2H7Nsh\no+WsuIdQVV3ulSBJ2iT80pckAQaCJKllIEiSAANBktQyEDRT5ua2kmTsx9zc1qF/BWnD8uJ2minr\ncQjlMIdvbtS6m7436uGbG7XuWT3s1D0ESRJgIEiSWgaCJAkwECRJLQNBkgQYCJKkloEgSQIMBElS\ny0CQJAEGgiSpZSBIkgADQZLUWvEWmpOW5EbgX4DdwF1V9aShapEkDRgINEEwX1W3DViDJKk15JBR\nBu5fkrTEkF/IBXwiyeeTvHLAOiRJDDtkdHxVfSvJT9MEw3VVddmeb9qxY8ePl+fn55mfn59ehdLU\nbGlvdCONbmFhgYWFhd7rmYk7piXZDtxRVW/b43XvmLaX2ZvvmGbfG6Vt/769Y9oSSQ5McnC7fBDw\nK8A1Q9QiSWoMNWR0JPAXSaqt4QNVdfFAtUiSmJEho5U4ZLT3ccjIvme/bf++HTKSJM00A0GSBBgI\nkqSWgSBJAgwESVLLQJAkAQaCJKllIEiSAANBktQyECRJgIEgSWoZCJIkwECQJLUMBEkSYCBIkloG\ngiQJMBAkSS0DQZIEGAiSpJaBIEkCDARJUstAkCQBBsKmNTe3lSRjPebmtg5dfg9bxv69pb1dqmro\nGlaUpGa5vlnWfMGNu+3CUNu9X90A/X7vYdra98Zq27/vSf99JaGqRv6/HPcQJEmAgSBJahkIkiTA\nQJAktQwESRJgIEiSWvsNXcBazjrrA1xwwcVjt3/JS17AiSf+6jpWJEmb08yfh/CYxzyFa675JeCR\nY6zhs5xwwq1ceOH/HruGubmt7Np101ht99nnQHbv/n9Tb3uPcf9tDwB+MHav/Wv3uHb7nuW2/fue\n1fMQZn4PofEs4JfGaHcA8JFePTdhMN4/3u7d439o+rRt9Dnz9ge9+u5Xu2cMS0NxDkGSBBgIkqTW\nYIGQ5IQk1yf5SpLXD1WHJKkxSCAk2Qd4J/BM4OeAFyc5ZohatJqFoQvoYWHoAnpaGLqAnhaGLqCn\nhaELGMRQewhPAr5aVTdV1V3Ah4ATB6pFK1oYuoAeFoYuoKeFoQvoaWHoAnpaGLqAQQwVCA8G/mHJ\n85vb1yRJA5n5w063bNmfgw56Pfvue8TIbe+66xa2bHn0BKqSpM1nkBPTkjwZ2FFVJ7TPTweqqs7c\n432ze9acJM2wcU5MGyoQ9gVuAJ4BfAv4HPDiqrpu6sVIkoCBhoyq6kdJXgtcTDOP8T7DQJKGNdPX\nMpIkTc9MnKm81klqSZ6b5EtJrkjyuSTHD1HncrqeYJfkiUnuSnLyNOtbS4dt/7Qktye5vH389hB1\nrqTL9k8y3352rklyybRrXE2H7f+f29ovT3J1kruTHDpErcvpUP/9k5yX5Mq2/pcNUOaKOtR/aJKP\ntN8/n0kyM0epJHlfkl1JrlrlPe9I8tV2+29bc6VVNeiDJpS+BhwN7A9cCRyzx3sOXLL8WOC6oevu\nWvuS930KuAA4eei6R9z2TwPOG7rWHvUfAlwLPLh9/lND1z3q52fJ+58DfHLoukfc/r8FnLG47YFb\ngf2Grn2E+t8CvLFdfuSMbf+nANuAq1b4+bOAj7XLvwB8Zq11zsIewponqVXV0mspHwzsnmJ9q+l6\ngt3rgA8D/zjN4jroWv+sXoK0S/2nAH9eVbcAVNU/T7nG1Yx6guaLgQ9OpbJuutRfwP3a5fsBt1bV\n3VOscTVd6n808GmAqroB2Jrkp6db5vKq6jLgtlXeciLwP9v3fhY4JMmRq61zFgKh00lqSU5Kch1w\nPvDvp1TbWtasPcmDgJOq6n8we1+sXU8Q/MV2l/Njs7TLTLf6HwEcnuSSJJ9P8htTq25tnU/QTHJf\n4ATgz6dQV1dd6n8n8Ogk3wS+BJw2pdq66FL/l4CTAZI8CXgocNRUqutvz9/vFtY4AXgWAqGTqvpo\nVT0KOAl489D1jOC/A0vHJmctFNbyReChVbWN5o/7owPXM6r9gONodp9PAN6Y5F8NW9JYfg24rKpu\nH7qQET0TuKKqHgQ8HnhXkoMHrmkUvwccluRy4DXAFcCPhi1pcmbhTOVbaFJ30VHta8uqqsuSPDzJ\n4VX17YlXt7outT8B+FCS0IyhPivJXVV13pRqXM2a9VfVnUuWL0zy7hnZ9tBt+98M/HNVfR/4fpJL\ngcfRjB0PbZTP/ouYreEi6Fb/y4EzAKrq75L8PXAM8IWpVLi6Lp//O1gyItHW//WpVNffLcBDljxf\n9bsVmIlJ5X25Z2LnPjQTO4/a4z0/s2T5OOAfhq67a+17vP8sZmtSucu2P3LJ8pOAG4eue8T6jwE+\n0b73QOBq4NFD1z7K54dmYvxW4L5D1zzG9n8XsH3xs0QzhHH40LWPUP8hwP7t8iuBs4eue4/6tgJX\nr/CzZ3PPpPKT6TCpPPgeQq1wklqSVzU/rvcA/y7JS4AfAt8DXjBcxffoWPu9mky9yFV0rP95SV4N\n3EWz7V84XMX31qX+qro+yUXAVTS7+u+pqi8PWPaPjfD5OQm4qKq+N1Sty+lY/5uBs5ccGvlfazb2\nLrvW/yjgnCS7aY5We8VwFd9bknOBeeCIJN8AttME2+Jn/6+SPDvJ14Dv0uytrb7ONj0kSXu5DTOp\nLEmaLANBkgQYCJKkloEgSQIMBElSy0CQJAEGgjahJHes03ouTHJbkvP2eH1reynkryT5YJI1z+dp\nr6V0XLt8QZL7r/Le05Ic0P83kEZjIGgzWq+Ta94C/Poyr58JvLWqHgHczognK1XVc6rqO6u85T/Q\nnFUtTZWBoL1CkqOTfKq9ausnkhzVvv7wJH/b3gDlTUv3LqrqEuDOZVb3dO656ug5wL9dpr8D2r2H\na5N8BDhgyc/+PsnhSQ5s9xauSHJVkucneR3wIOCSJJ9avy0grc1A0N7iD4Czqrlq67ntc4C3A/+t\nqh5HcyG8VfcukhwB3FZVi/fkuJnmC3xPrwa+W1U/R3NJgScs+dliHycAt1TV46vqWODjVfUHNBcg\nm6+qZ4z6S0p9GAjaW/wi91wt9H8Bxy95/cPt8rnr2N9TgfcDVNXVNNfVX7R4CfSrgV9OckaSp1Rz\nZc3Fn2+0y6RrEzAQtLfoMq+w5pdwVd0KHJpk8W9n7UsKr7DuqvoqzdV7rwbenBm7X7X2PgaCNqPl\nvtj/huYWlNBMFP91u/y3wPPa5RetsK4913cJ8Px2+aXAXy7T7lLgVIAkjwGO/YkVJw8EvldV5wK/\nTxMOAN8BVjwKSZoUr3aqTSfJ3cA3ab7IC3gbzSTw2cARwD8BL6+qm9u7p72fZtL3IuCUqnpIu55L\naW6sfjDN/QheUVWfSPIwmvvvHkZzB61fr+aevEtrOIDm/hfHAtfR3LrwNVV1eZKv08wpPIEmCHbT\nXNr91e3PXwu8lmZ+wXkETY2BoL1akvsu3mcgyQuBF1XVTxw1JO0NBr9BjjSwn0/yTpq9idtYcrtE\naW/jHoIkCXBSWZLUMhAkSYCBIElqGQiSJMBAkCS1DARJEgD/Hxj5iVt/1SilAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fcf0d664c90>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "139\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEPCAYAAABCyrPIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF19JREFUeJzt3XuwpHV95/H3h/tF5ZbIMYKMuBo0AUaixCyUnmgZwTXq\nsl4hWWRdyrLEsLu1u6IVM1OltQRTkvUSd+NqkN2IVkISBAwCCr2EjXdALgKiBhLQmSQKCl7B+e4f\n/Rw5jNPnPN3n9PP0OfN+VXXN83T37/l9p8+Z/szz+z2XVBWSJO3SdwGSpNlgIEiSAANBktQwECRJ\ngIEgSWoYCJIkYMqBkOSQJFcluSXJTUne2Dy/KcndSa5rHidMsw5J0vIyzfMQkswBc1V1Q5JHAV8E\nXgK8Eri/qs6dWueSpLHsNs2NV9UWYEuz/ECSW4HHNy9nmn1LksbT2RxCkg3ARuCzzVNnJLkhyQeS\n7NdVHZKkHeskEJrhoguBM6vqAeB9wOFVtZHhHoRDR5LUs6nOIQAk2Q24FLisqt61g9cPAy6pqqN2\n8JoXWpKkCVTV2MPyXewh/Anw5cVh0Ew2LzgJuHlU46qa+cemTZt6r8E6rbHrOpt/oSt40Emda+Xz\nXP2fzfimOqmc5DjgFOCmJNcz/C14C3Byko3ANuBO4HXTrEOStLxpH2X0/4Bdd/DSJ6bZryRpfJ6p\nvArm5+f7LqEV61w9a6FGsM7VtlbqnNTUJ5VXIknNcn3SziwJC3MBE25hRePdGi0JNaOTypKkNcBA\nkCQBBoIkqWEgSJIAA0GS1DAQJEmAgSBJahgIkiTAQJAkNQwESRJgIEiSGgaCJAkwECRJDQNBkgQY\nCJKkhoEgSQIMBElSw0CQJAEGgiSpYSBIkgADQZLUMBAkSYCBIElqGAiSJMBAkCQ1DARJEmAgSJIa\nBoIkCTAQJEkNA0Haic3NbSDJRI+V23PivufmNqxC/9peqqrvGkZKUrNcn7TWDb/YJ/03tpK2K20f\n/G4YLQlVNXZqu4cgSQIMBElSw0CQJAFTDoQkhyS5KsktSW5K8jvN8wckuSLJ7UkuT7LfNOuQJC1v\nqpPKSeaAuaq6IcmjgC8CLwFOA75VVe9I8ibggKo6awftnVSWpshJ5fVpJieVq2pLVd3QLD8A3Aoc\nwjAUzm/edj7w0mnWIUlaXmdzCEk2ABuBzwAHV9VWGIYG8Niu6pAk7dhuXXTSDBddCJxZVQ8k2X5f\nb+S+3+bNm3+6PD8/z/z8/DRKlHozN7eBrVvvmrj9wQcfxpYtd65eQVpzBoMBg8FgxduZ+olpSXYD\nLgUuq6p3Nc/dCsxX1dZmnuHqqnrqDto6h6B1b2Xj+LCS8XTnENanmZxDaPwJ8OWFMGhcDLymWT4V\n+FgHdUiSljDto4yOA64BbmL4X4EC3gJ8Dvgz4FDgLuAVVXXfDtq7h6B1zz2Eydr63TDapHsIXstI\n6pmBMFlbvxtGm+UhI0nSGmAgSJIAA0GS1DAQJEmAgSBJahgIkiTAQJAkNQwESRJgIEiSGgaCJAkw\nECRJDQNBkgQYCJKkhoEgSQIMBElSw0CQJAEGgiSpYSBIkgADQZLUMBAkSYCBIElqjBUISXZJ8php\nFSNJ6s+ygZDkgiSPSbIvcDPw5ST/ZfqlSWpnT5JM9ND45uY2TPx5z81t6Lv8JaWqln5DckNVbUxy\nCnAMcBbwxao6aurFJbVcfdJaN/xiXsnv+Ura99V25X339d2wsp9XN3UnoarGTvw2Q0a7J9kdeClw\ncVU9yMp+CyRJM6hNIPwxcCewL3BNksOA706zKElS95YdMvqZBsP9pV2r6qHplPSIvhwy0rrnkNFk\nbR0yWqKXaQ0ZJflvSfZf9NT+wOZxO5IkzbY2Q0YnVtV9CytVdS/wwumVJEnqQ5tA2DXJngsrSfYG\n9lzi/ZKkNWi3Fu/5MPCpJOc166cB50+vJElSH1pNKic5EXhes3plVV0+1aoe7tdJZa0Jc3Mb2Lr1\nrhVswUnlcds6qbxELxNOKo99lFGXDAStFSv9kjAQxm9rICzRy4SBMHLIKMm1VXV8kvt55N8+QFWV\n1zSSpHXEPQRpFbiH0H3f7iEs0csUL11BkmOS/E6SNyZ5+hhFfTDJ1iQ3LnpuU5K7k1zXPE4Yt2hJ\n0uprc2La7zE8qugg4OeADyX53ZbbPw94wQ6eP7eqjmken2hdrSRpatocdnoKcHRV/RAgye8DNwBv\nX65hVV3bXPtoe153V5JmTJsho28Aey1a3xO4Z4X9npHkhiQfSLLfCrclSVoFbQLhO8AtST7UnJx2\nM3BfkncnefcEfb4POLyqNgJbgHMn2IYkaZW1GTL6q+axYLCSDqvqnxat/i/gkqXev3nz5p8uz8/P\nMz8/v5LuJWndGQwGDAaDFW+n7ZnKewBPaVZvb26S066DZANwSVUd2azPVdWWZvk/As+sqpNHtPWw\nU60JHnbafd8edrpEL6t9YtqiDc8zPMroToY/wUOTnFpV17RoewEwDxyU5O+BTcCvJ9kIbGu2+bpx\ni5Ykrb4291T+InByVd3erD8F+EhV/crUi3MPQWuEewjd9+0ewhK9TPOeygthAFBVXwF2H7cjSdJs\nazOp/IUkHwD+tFk/BfjC9EqSJPWhzZDRnsAbgOObp/4G+KOq+vGUa3PISGuGQ0bd9+2Q0RK9TOvy\n10nOrKp3LffcNBgIWisMhO77NhCW6GWKcwin7uC514zbkSRpti11P4RXAycDT0xy8aKXHg18e9qF\nSZK6tdSk8t8C32R4hdN3Lnr+fuDGHbaQJK1Z3iBHWgXOIXTft3MIS/QyzRvkSJLWPwNBkgQsEQhJ\nPtX8eU535UiS+rLUpPLjkvxL4MVJPsp2dzmrquumWpkkqVMjJ5WTvAx4LcMzlLe/VEVV1XOnXJuT\nyloznFTuvm8nlZfoZYpnKr+1qt42cWUrYCBorTAQuu/bQFiil2kFQrPxFwPPblYHVXXpuB1NwkDQ\nWmEgdN1+L+BHE7XcZZd92Lbt+xP2u2AnDYQkZwPHAh9unno18PmqesvYVY7JQNBaYSCspb7X/57N\nNAPhRmBjVW1r1ncFrq+qoyaqdJziDAStEQbCWurbQBil7XkI+y9a3m/cTiRJs6/NDXLOBq5PcjXD\naHw2cNZUq5Ikda7tpPLjgGc2q5+rqi1Trerhfh0y0prgkNFa6tsho5HtZvkL10DQWmEgrKW+DYRR\nvJaRJAkwECRJjSUDIcmuSW7rqhhJUn+WDISq+glwe5IndFSPJKknbQ47PQC4JcnngO8tPFlVL55a\nVZKkzrUJhLdOvQpJUu/anodwGPDkqvpkkn2AXavq/qkX52GnWiM87HQt9e1hp6Mse5RRktOBC4E/\nbp56PHDRuB1JkmZbm8NO3wAcB3wXoKruAB47zaIkSd1rEwg/qqofL6wk2Y2V7W9JkmZQm0D4v0ne\nAuyd5PnAnwOXTLcsSVLX2twPYReG91b+DYazKZcDH+hittdJZXVpbm4DW7fetYItrM1JTuvutu9Z\nnlRue5TRHsARDD+F2xcPIU2TgaAu9XekkF+sa6ftyvue5UBY9jyEJP8K+J/A1xh+Ek9M8rqqumz8\nMiVJs6rNkNFtwIuq6qvN+pOAj1fVEVMvzj0Edcg9hC7b9tm3ewijtJlUvn8hDBpfB1qdlJbkg0m2\nNvdlXnjugCRXJLk9yeVJvCWnJM2AkYGQ5KQkJwFfSPLXSV6T5FSGRxh9vuX2zwNesN1zZwGfrKpf\nBK4C3jxB3ZKkVTZyyCjJeUs1rKrTWnUwvOzFJVV1VLN+G/CcqtqaZA4YjBp+cshIXXLIqMu2ffbt\nkNEoIyeV237hT+CxVbW16WNLEs96lqQZ0OYooycCbwQ2LH7/Kl7+2l0ASZoBbS5/fRHwQYZzB9tW\noc+tSQ5eNGT0j0u9efPmzT9dnp+fZ35+fhVKkKT1YzAYMBgMVrydNoedfraqfnXiDpINDOcQjmzW\nzwG+XVXnJHkTcEBVnTWirXMI6oxzCF227bNv5xBGtmsRCCcDTwauAH608HxVXdeiqAuAeeAgYCuw\nieEex58DhwJ3Aa+oqvtGtDcQ1BkDocu2ffZtIIxs1yIQzgZ+m+GZygtDRlVVzx27yjEZCOqSgdBl\n2z77NhBGaTOH8HLg8K6uXyRJ6kebM5VvBvafdiGSpH612UPYH7gtyed55BzCah12KkmaAW0CYdPU\nq5Ak9a7V/RD64qSyuuSkcpdt++zbSeVR2pypfD8P/+33AHYHvldVjxm3M0nS7Fo2EKrq0QvLGf4X\n6iXAs6ZZlCSpexMNGSW5vqqePoV6tu/HISN1xiGjLtv22bdDRqO0GTI6adHqLsAzgB+O25Ekaba1\nOcroNxctPwTcyXDYSJK0jniUkdRwyKjLtn327ZDRKCP3EJL83hLtqqreNm5nkqTZtdSQ0fd28Ny+\nwGsZXr3UQJCkdaTVkFGSRwNnMgyDPwPeWVVL3thmNThkpC45ZNRl2z77dsholCUnlZMcCPwn4BTg\nfOCYqrp3shIlSbNsqTmEPwBOAt4PHFlVD3RWlSSpcyOHjJJsY3h104d45P5RGE4qT/3SFQ4ZqUsO\nGXXZts++HTIaZeQeQlW1uVeCJGmd8EtfkgQYCJKkhoEgSQIMBElSw0CQJAEGgiSpYSBIkgADQZLU\nMBAkSYCBIElqGAiSJMBAkCQ1DARJEmAgaB2Zm9tAkokf0s6u1S00++L9EDSOld3PAHbW6/Nbd7d9\nz/L9ENxDkCQBBoIkqWEgSJKAJW6hOW1J7gS+A2wDHqyqY/uqRZLUYyAwDIL5qrq3xxokSY0+h4zS\nc/+SpEX6/EIu4Mokn09yeo91SJLod8jouKr6ZpKfZxgMt1bVtdu/afPmzT9dnp+fZ35+vrsKJWkN\nGAwGDAaDFW9nJk5MS7IJuL+qzt3ueU9MU2uemLaW2vbZtyemjdLLkFGSfZI8qlneF/gN4OY+apEk\nDfU1ZHQw8FdJqqnhw1V1RU+1SJKYkSGjURwy0jgcMlpLbfvs2yGjUTzsU5IEGAiSpIaBIEkCDARJ\nUsNAkCQBBoIkqWEgSJIAA0GS1DAQJEmAgSBJahgIkiTAQJAkNQwESRJgIEiSGgaCJAkwECRJDQNB\nkgQYCJKkhoEgSQIMBElSw0CQJAGwW98FtPGxj13KlVdeNXH7008/laOPPnoVK5Kk9SdV1XcNIyWp\nquLoo5/DjTf+MvCkCbYy4Iwznsx73vPOiWqYm9vA1q13TdQW4OCDD2PLljsnbj+pnbVuWMnvc1bQ\nvq+2ffZt3ZO07eI7NwlVlXHbrYk9hKFXAM+ZoF0B35i41+GX0+Q/wK1bx/6ZrIqds+5+apbWC+cQ\nJEmAgSBJahgIkiTAQJAkNQwESRJgIEiSGgaCJAkwEDTSniSZ6DE3t6Hv4iVNYA2dmKZu/YhJTxDr\n66Q2SSvjHoIkCTAQJEmN3gIhyQlJbkvylSRv6qsOSdJQL4GQZBfgvcALgF8CXp3kiD5qWR2Dvgto\nadB3AS0N+i6ghUHfBbQ06LuAlgZ9F9DSoO8CpqqvPYRjgTuq6q6qehD4KPCSnmpZBYO+C2hp0HcB\nLQ36LqCFQd8FtDTou4CWBn0X0NKg7wKmqq9AeDzwD4vW726ekyT1ZE0cdrrnnruz775vZtddDxq7\n7Y9//DX22OPFU6hKktaXXu6YluRZwOaqOqFZPwuoqjpnu/fN7u3cJGmGTXLHtL4CYVfgduB5wDeB\nzwGvrqpbOy9GkgT0NGRUVT9JcgZwBcN5jA8aBpLUr172ECRJs2cmzlRe7iS1JCcn+VLzuDbJkTNY\n44ub+q5P8rkkx3VdY5s6F73vmUkeTHJSl/Ut6n+5z/M5Se5Lcl3z+N1ZrLN5z3zzc785ydVd19jU\nsNzn+Z+bGq9LclOSh5LsP4N1PibJxUluaOp8Tdc1NnUsV+f+Sf6y+Tf/mSRP66HGDybZmuTGJd7z\n7iR3NJ/nxmU3WlW9PhiG0leBw4DdgRuAI7Z7z7OA/ZrlE4DPzGCN+yxaPhK4dRY/y0Xv+xRwKXDS\nLNYJPAe4uI/fyTHr3A+4BXh8s/5zs1jndu9/EfDJWawTeDNw9sJnCXwL2G0G63wH8NZm+Rd7+jyP\nBzYCN454/UTg483yr7b53pyFPYRlT1Krqs9U1Xea1c/Q/TkLbWr8/qLVRwHbOqxvQdsT/t4IXAj8\nY5fFLdK2zr4vm9qmzpOBv6iqewCq6p87rhHGP9Hz1cBHOqnskdrUWcCjm+VHA9+qqoc6rBHa1fk0\n4CqAqrod2JDk57sssqquBe5d4i0vAf53897PAvslOXipbc5CIIx7ktq/By6bakU/q1WNSV6a5Fbg\nEuDfdVTbYsvWmeQXgJdW1f+gvy/ctj/zX2t2dT/exy457ep8CnBgkquTfD7Jb3dW3cNa/xtKsjfD\nvey/6KCu7bWp873A05J8A/gScGZHtS3Wps4vAScBJDkWeAJwSCfVtbf93+MelvnP9Jo4MW1Bkl8H\nTmO4qzRzquoi4KIkxwNvB57fc0k78t+BxWOiff8vfJQvAk+oqu8nORG4iOGX76zZDTgGeC6wL/Dp\nJJ+uqq/2W9ZIvwlcW1X39V3ICC8Arq+q5yZ5EnBlkqOq6oG+C9vO7wPvSnIdcBNwPfCTfktauVkI\nhHsYpuuCQ5rnHiHJUcD7gROqaqndpGloVeOCqro2yeFJDqyqb0+9uoe1qfMZwEeThOEY7YlJHqyq\nizuqEVrUufgLoKouS/K+Gf087wb+uap+CPwwyTXA0QzHoLsyzu/nq+hnuAja1XkacDZAVX0tyd8B\nRwBf6KTCoTa/n/ezaBSgqfPrnVTX3j3AoYvWl/zeAmZiUnlXHp7A2YPhBM5Tt3vPE4A7gGfNcI1P\nWrR8DPAPs1jndu8/j34mldt8ngcvWj4WuHNG6zwCuLJ57z4M/7f4tFmrs3nffgwnaffu+rMc4/P8\nI2DTwu8AwyGPA2ewzv2A3Zvl04EP9fSZbgBuGvHaC3l4UvlZtJhU7n0PoUacpJbkdcOX6/3AW4ED\ngfc1/7N9sKqOnbEa/02Sfwv8GPgB8Iqu6huzzkc06bpGaF3ny5K8HniQ4ef5ylmss6puS3I5cCPD\nIYP3V9WXZ63O5q0vBS6vqh90Wd+Ydb4d+NCiQyn/a3W7V9i2zqcC5yfZxvAos9d2WSNAkguAeeCg\nJH8PbGIYYAu/m3+d5IVJvgp8j+He19LbbNJDkrSTm4WjjCRJM8BAkCQBBoIkqWEgSJIAA0GS1DAQ\nJEmAgaB1KMn9q7Sdy5Lcm+Ti7Z7f0Fzy+CtJPpJk2fN5mmsdHdMsX5rkMUu898wke638byCNx0DQ\nerRaJ9e8A/itHTx/DvDOqnoKcB9jnpRUVS+qqu8u8Zb/wPCsZ6lTBoJ2CkkOS/Kp5uqpVyY5pHn+\n8CSfbm508rbFexdVdTWwo4uqPZeHrxZ6PvCvd9DfXs3ewy1J/hLYa9Frf5fkwCT7NHsL1ye5McnL\nk7wR+AXg6iSfWr1PQFqegaCdxXuA86pqI3BBsw7wLuAPq+pohheqW3LvIslBwL1VtXC/i7sZfoFv\n7/XA96rqlxheUuAZi15b6OME4J6qenpVHQV8oqrew/ACZPNV9bxx/5LSShgI2ln8Gg9f5fP/AMct\nev7CZvmCVezv2cCfAlTVTQyvn79g4ZLjNwHPT3J2kuNreAXNhddn9bLkWscMBO0s2swrLPslXFXf\nAvZPsvBvZ/lLCo/YdlXdwfDKuDcBb09P942WFhgIWo929MX+twxvHQnDieK/aZY/DbysWX7ViG1t\nv72rgZc3y6cCH9tBu2uAUwCS/DJw1M9sOHkc8IOqugD4A4bhAPBdYORRSNK0eLVTrTtJHgK+wfCL\nvIBzGU4Cfwg4CPgn4LSqujvJv2A4tLMXcDlwclUd2mznGoY3UH8Uw/sIvLaqrkzyRIb32T2A4Z2y\nfquG995dXMNeDO83cRRwK8NbF76hqq5L8nWGcwrPYBgE2xheNv31zetnAGcwnF9wHkGdMRC0U0uy\n98L9AZK8EnhVVf3MUUPSzqD3G+RIPfuVJO9luDdxL4tuiyjtbNxDkCQBTipLkhoGgiQJMBAkSQ0D\nQZIEGAiSpIaBIEkC4P8DShNfvuqAlhMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fcf0ce4ed50>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "179\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEPCAYAAABCyrPIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF8ZJREFUeJzt3X2QZXV95/H3h5lxHFQehECrIKPZKOqqSJSNi6W9GiMm\nRlh8xrhgLMuyfGDX2l3RipnZ1YrRlG406m5cXZ3EoGtMVHwKoHKXkETRADIioMZAAsqoCOFBQXC+\n+8c9zTRjT/e5ffvcc7v7/aq6xbnn3nN+375038+c8zu/30lVIUnSfn0XIEmaDgaCJAkwECRJDQNB\nkgQYCJKkhoEgSQI6DoQkm5N8OcnFSXYm2dasPzjJOUmuTHJ2kgO7rEOStLR0PQ4hyf5V9eMkG4C/\nAV4NPAu4vqremuS1wMFVdUanhUiSFtX5KaOq+nGzuBnYCBRwIrCjWb8DOKnrOiRJi+s8EJLsl+Ri\n4Drg3Kr6CnB4Ve0CqKrrgMO6rkOStLhJHCHsrqrHAEcAxyV5BMOjhLu9res6JEmL2ziphqrqpiQD\n4ARgV5LDq2pXkhng+wttk8SgkKRlqKqMuk3XVxkdOncFUZItwFOBy4GzgNOat50KfHJf+6iqVfvY\ntm1b7zWs1/pXc+3W3/9jtde/XF0fIdwP2JFkP4bh83+r6rNJvgR8NMlvA1cDz+24DknSEjoNhKra\nCRy7wPofAb/aZduSpNE4UrlDs7OzfZcwltVc/2quHay/b6u9/uXqfGDaOJLUNNcnSdMoCTVtncqS\npNXDQJAkAQaCJKlhIEiSAANBktQwECRJgIEgSWoYCJIkwECQJDUMBEkSYCBIkhoGgiQJMBAkSQ0D\nQZIEGAiSpIaBIEkCDARpTZiZ2UqSZT02bLjXsrabmdna94+tFeYd06Q1IAmw3L+V5W4b/PucTt4x\nTZI0FgNBkgQYCJKkhoEgSQIMBElSw0CQJAEGgiSpYSBIkoCOAyHJEUm+mOSyJDuTvKpZvy3JNUku\nah4ndFmHJGlpnY5UTjIDzFTVJUnuDfw9cCLwPODmqnr7Ets7UllqwZHKmm+5I5U3dlHMnKq6Driu\nWb4lyeXAA5qXRy5WktSdifUhJNkKHAN8uVn1yiSXJHlfkgMnVYckaWETCYTmdNHHgNOr6hbgPcCD\nq+oYhkcQi546kiR1r9NTRgBJNjIMgz+tqk8CVNUP5r3lfwOf2tf227dvv2t5dnaW2dnZTuqUpsHM\nzFZ27bq67zK0ygwGAwaDwdj76Xz66yR/Avywql4zb91M079Akv8EPK6qTllgWzuVta4sv3PYTmXt\nsdxO5a6vMjoeOB/YyfA3roDXA6cw7E/YDVwFvKyqdi2wvYGgdcVA0EqYykAYl4Gg9cZA0ErwBjmS\npLEYCJIkwECQJDUMBEkSYCBIkhoGgiQJMBAkSQ0DQZIEGAiSpIaBIEkCDARJUsNAkCQBBoIkqWEg\nSJIAA0GS1DAQJEmAgSBJahgIkiTAQNAaNzOzlSTLeszMbO27fGmivKey1rTl36MY+rhnsPdU1krw\nnsqSpLEYCJIkwECQJDUMBEkSYCBIkhoGgiQJMBAkrRPLHZOynsajOA5Ba5rjELrcdnWNQxjns11N\nPyc4DkGSNKaRAiHJfkkOGOH9RyT5YpLLkuxM8upm/cFJzklyZZKzkxw4auGSpJW1ZCAkOTPJAUnu\nBXwd+EaS/9Jy/3cCr6mqRwCPB16R5GjgDODzVfVQ4IvA65ZXviRppbQ5Qnh4Vd0EnAR8DngQ8KI2\nO6+q66rqkmb5FuBy4AjgRGBH87Ydzb4lST1qEwibkmxi+KV9VlXdwTJ6ZpJsBY4BvgQcXlW7YBga\nwGGj7k+StLI2tnjPHwNXAV8Dzk9yFHDTKI0kuTfwMeD0qrolyd6Bss+A2b59+13Ls7OzzM7OjtK0\nJK15g8GAwWAw9n5Gvuw0w2u3NlTVnS3fvxH4NPC5qnpHs+5yYLaqdiWZAc6rqoctsK2XnWosXnba\n5bar63JMLztdWptO5d9LctC8VQcB20do4/8A35gLg8ZZwGnN8qnAJ0fYnySpA0seISS5uKoes9e6\ni6rq2CV3nhwPnA/sZBjNBbweuBD4KHAkcDXw3Kq6cYHtPULQWDxC6HLb1fUvZ48QltamD2FDks1V\ndXvT0BZgc5udV9XfABv28fKvtitRkjQJbQLhz4AvJPlA8/zF7LlkVJK0RrTqVE7ydOApzdNzq+rs\nTqva066njDQWTxl1ue3qOpXiKaMW203zD2ogaFx9BMLMzFZ27bp6mW2CgdANA6HFdvv6QZNcUFVP\nSHIzd/8UA1RVtZ7TaLkMBI2rj0AYt00DoRsGQovtpvkHNRA0LgOh2zZX09+ngbC0Np3KJDkWeALD\nT/OCqrp41IYkSdOtzcC032V4VdEhwKHAB5P8TteFSZImq83AtCuBR1fVbc3zLcAlzdTV3RbnKSON\nyVNG3ba5mv4+PWW0tDaznX4XuOe855uBa0dtSJI03dr0IfwLcFmScxnG61OBC5O8E6CqXt1hfZKk\nCWkTCB9vHnMG3ZQiSepT25HK9wAe0jy9srlJTufsQ9C47EPots3V9PdpH8LSljxCSDLL8Cqjqxj+\n5hyZ5NSqOn/UxiRJ06vNKaO3Ab9WVVcCJHkI8GHgl7ssTJI0Wa3uqTwXBgBV9U1gU3clSZL60OYI\n4atJ3gd8qHn+QuCr3ZUkSepDm4Fpm4FXMJy6AuCvgXdX1U87rs1OZY3NTuVu21xNf592KrfYrkUg\nnL7X/ZAXXNcFA0HjMhC6bXM1/X0aCEtr04dw6gLrThu1IUnSdNtnH0KSFwCnAA9Kcta8l+4D/Kjr\nwiRJk7VYp/LfAt9jOMPp2+atvxm4tMuiJEmT5w1ytKbZh9Btm6vp79M+hKW16UOQJK0DBoIkCVgk\nEJJ8ofnvWyZXjiSpL4t1Kt8vyb8FnpnkIwxPNN6lqi7qtDJJU25zc15+dIcffhTXXXfVypajse2z\nUznJs4GXMByhvPdUFVVVT+64NjuVNTY7lae3zUn/bdup3GK7FiOV31BVb1x2ZWMwEDQuA2F62zQQ\nutNZIDQ7fybwxObpoKo+3bKo9wPPAHZV1aOadduAlwLfb972+qr6q31sbyBoLAbC9LZpIHSns8tO\nk7wZOB34RvM4Pcnvtdz/B4CnLbD+7VV1bPNYMAwkSZPVZvrr3wCOqardAEl2ABcDr19qw6q6IMlR\nC7y0vJ4oSVJn2o5DOGje8oEr0O4rk1yS5H1JVmJ/kqQxtTlCeDNwcZLzGP7L/onAGWO0+R7gv1dV\nJXkT8HaGVzMtaPv27Xctz87OMjs7O0bT0iiWf1mllrK8z9bLVRc2GAwYDAZj76dtp/L9gMc1Ty+s\nqutaNzA8ZfSpuU7ltq81r9uprLHYwbv22lzud4Kdyktrc4RAVX0POGvJNy4szOszSDIzL1BOBr6+\nzP1KklZQq0BYriRnArPAIUn+CdgG/LskxwC7gauAl3VZgySpHae/1prmKaO116anjJbWyTiEJBuS\nXLH8siRJq8WigVBVPwOuTPLACdUjSepJmz6Eg4HLklwI3Dq3sqqe2VlVkqSJaxMIb+i8CklS79qO\nQzgK+KWq+nyS/YENVXVz58XZqawx2am89tq0U3lpXU5u91LgY8AfN6seAHxi1IYkSdOtzVxGrwCO\nB24CqKpvAYd1WZQkafLaBMLtVfXTuSdJNrL840RJ0pRqEwj/L8nrgS1Jngr8OfCpbsuSJE1am1to\n7sdwNtJfY9gTdDbwvkn09tqprHHZqbz22rRTeWld30LzHsDRDD/NK+efQuqSgaBxGQhrr00DYWmd\nzXaa5DeA/wX8A8P/iw9K8rKq+tzoZUqSplWbU0ZXAM+oqm83z38R+ExVHd15cR4haEweIay9Nj1C\nWFpn4xCAm+fCoPEdoPNBaZKkydrnKaMkJzeLX03yWeCjDOP1OcBXJlCbJGmCFutD+M15y7uAJzXL\nPwC2dFaRJKkX3iBHa5p9CGuvTfsQltblVUYPAl4FbJ3/fqe/lqS1pc30158A3s9wdPLubsuRJPWl\nTSDcVlXv7LwSSVKv2oxDOAX4JeAc4Pa59VV1Ubel2Yeg8dmHsPbatA9haZ31IQCPBF4EPJk9p4yq\neS5JWiPaBMJzgAdPav4iSVI/2oxU/jpwUNeFSJL61eYI4SDgiiRf4e59CF52KklrSJtA2NZ5FZKk\n3jlSWWuaVxmtvTa9ymhpXY5Uvpk9n+I9gE3ArVV1wKiNSZKm15KdylV1n6o6oAmALcCzgPe02XmS\n9yfZleTSeesOTnJOkiuTnJ3kwGVXL0laMW2uMrpLDX0CeFrLTT6wwHvPAD5fVQ8Fvgi8bpQaJEnd\naHPK6OR5T/cDHgvc1mbnVXVBkqP2Wn0ie6bS3gEMGIaEJKlHba4ymn9fhDuBqxh+qS/XYVW1C6Cq\nrkty2Bj7kiStkCUDoape3HENi3bfb9++/a7l2dlZZmdnOy5H0vTa3FwtpPkGgwGDwWDs/ezzstMk\nv7vIdlVVb2zVwPCU0aeq6lHN88uB2aralWQGOK+qHraPbb3sVGPxslPbXIk2V9v30HIvO12sU/nW\nBR4ALwFeO0ptzWPOWcBpzfKpwCdH2JckqSOtBqYluQ9wOsMw+Cjwtqr6fovtzgRmgUMY3pd5G8Mb\n7vw5cCRwNfDcqrpxH9t7hKCxeIRgmyvR5mr7HlruEcKigZDkvsBrgBcyvCLoHVV1w7KrHJGBoHEZ\nCLa5Em2utu+hFR+pnOQPgJOB9wKPrKpbxqhPkjTlFutU3s1wdtM7uXushmGncudTV3iEoHF5hGCb\nK9HmavseWvEjhKoaaRSzJGl180tfkgQYCJKkhoEgSQIMBElSw0CQJAEGgiSpYSBIkgADQavEzMxW\nkoz8kNReq8nt+uJIZc1Z/ojj9TOa1ja7a3O1fQ91Mf21JGkdMRAkSYCBIElqGAiSJMBAkCQ1DARJ\nEmAgaIKWO5bA8QTq1+Zl/97OzGztu/iROA5BE+Pdy2xzPbbZx3eY4xAkSWMxECRJgIEgSWoYCJIk\nwECQJDUMBEkSYCBIkhoGgiQJgI19NZzkKuBfgN3AHVV1XF+1SJJ6DASGQTBbVTf0WIMkqdHnKaP0\n3L4kaZ4+v5ALODfJV5K8tMc6JEn0e8ro+Kr6XpJfYBgMl1fVBT3WI0nrWm+BUFXfa/77gyQfB44D\nfi4Qtm/fftfy7Owss7OzE6pQklaHwWDAYDAYez+9TH+dZH9gv6q6Jcm9gHOA/1ZV5+z1Pqe/XkOc\n/to212Obq2n6676OEA4HPp6kmhr+bO8wkCRNljfI0cR4hGCb67HN1XSE4GWfkiTAQJAkNQwESRJg\nIEiSGgaCJAkwECRJDQNBkgQYCJKkhoEgSQIMBElSw0DQyGZmtpJk5Iek6eZcRhrZ8uckWn3z0Nim\nbTqXkSRp3TEQJEmAgSBJahgIkiTAQJAkNQwESRJgIKyo5V6fv2HDvZa13Tjbzsxs7fvjkjRlHIew\nglbb9fnL/WxX289pm7bZZ5uOQ5AkrToGgiQJMBAkSQ0DQZIEwMa+C+jCbbfdxjXXXLOsbQ899FAO\nOuigFa5IkqbfmgyE008/gx07zmTTpgNG2m737p9y5JGHccUVX+2oMkmaXmsyEK6//kZuv/2t3H77\naSNu+W1uuOGELkqSpKlnH4IkCegxEJKckOSKJN9M8tq+6pAkDfUSCEn2A94FPA14BPCCJEf3UUu3\nBn0XMJbBYNB3CWMY9F3AmAZ9FzCmQd8FjGnQdwG96OsI4TjgW1V1dVXdAXwEOLGnWjo06LuAsRgI\nfRr0XcCYBn0XMKZB3wX0oq9AeADwz/OeX9OskyT1ZE1eZbR58ya2bPlDNm36i5G22737VjZv3tRR\nVZI03XqZ7TTJrwDbq+qE5vkZQFXVW/Z63+qZ6lSSpshyZjvtKxA2AFcCTwG+B1wIvKCqLp94MZIk\noKdTRlX1sySvBM5h2I/xfsNAkvo11TfIkSRNzlSMVF5qkFqSZyb5WpKLk1yY5Pg+6lxI2wF2SR6X\n5I4kJ0+yvqW0+OyflOTGJBc1j9/po859afP5J5ltfne+nuS8Sde4mBaf/39uar8oyc4kdyaZmtkX\nW9R/QJKzklzS1H9aD2XuU4v6D0ryl833z5eSPLyPOheS5P1JdiW5dJH3vDPJt5rP/5gld1pVvT4Y\nhtK3gaOATcAlwNF7vWf/ecuPBC7vu+62tc973xeATwMn9133iJ/9k4Cz+q51jPoPBC4DHtA8P7Tv\nukf9/Zn3/mcAn++77hE//9cBb5777IHrgY191z5C/W8F3tAsP3TKPv8nAMcAl+7j9acDn2mW/w3w\npaX2OQ1HCEsOUquqH897em9g9wTrW0zbAXavAj4GfH+SxbXQtv6Rr1aYkDb1nwL8RVVdC1BVP5xw\njYsZdYDmC4APT6SydtrUX8B9muX7ANdX1Z0TrHExbep/OPBFgKq6Etia5BcmW+bCquoC4IZF3nIi\n8CfNe78MHJjk8MX2OQ2B0GqQWpKTklwOfAr47QnVtpQla09yf+CkqvqfTN8Xa9sBgo9vDjk/M02H\nzLSr/yHAfZOcl+QrSV40seqW1nqAZpItwAnAaINrutWm/ncBD0/yXeBrwOkTqq2NNvV/DTgZIMlx\nwAOBIyZS3fj2/vmuZYkBwNMQCK1U1Seq6mHAScCb+q5nBH8IzD83OW2hsJS/Bx5YVccw/OP+RM/1\njGojcCzDw+cTgDck+Vf9lrQsvwlcUFU39l3IiJ4GXFxV9wceA7w7yb17rmkUvw8cnOQi4BXAxcDP\n+i2pO9MwUvlahqk754hm3YKq6oIkD05y36r6UefVLa5N7Y8FPpIkDM+hPj3JHVV11oRqXMyS9VfV\nLfOWP5fkPVPy2UO7z/8a4IdVdRtwW5LzgUczPHfct1F+95/PdJ0ugnb1vxh4M0BV/UOSfwSOBqbh\nLlRtfv9vZt4Ziab+70ykuvFdCxw57/mi363AVHQqb2BPx849GHbsPGyv9/zivOVjgX/uu+62te/1\n/g8wXZ3KbT77w+ctHwdc1XfdI9Z/NHBu8979gZ3Aw/uufZTfH4Yd49cDW/queRmf/7uBbXO/SwxP\nYdy379pHqP9AYFOz/FLgg33XvVd9W4Gd+3jt19nTqfwrtOhU7v0IofYxSC3Jy4Yv13uBZyX5D8BP\ngZ8Az+2v4j1a1n63TSZe5CJa1v/sJC8H7mD42T+vv4rvrk39VXVFkrOBSxke6r+3qr7RY9l3GeH3\n5yTg7Kr6SV+1LqRl/W8CPjjv0sj/WtNxdNm2/ocBO5LsZni12kv6q/jukpwJzAKHJPknYBvDYJv7\n3f9skl9P8m3gVoZHa4vvs0kPSdI6t2o6lSVJ3TIQJEmAgSBJahgIkiTAQJAkNQwESRJgIGgNSnLz\nCu3nc0luSHLWXuu3NlMhfzPJh5MsOZ6nmUvp2Gb500kOWOS9pye55/g/gTQaA0Fr0UoNrnkr8FsL\nrH8L8LaqeghwIyMOVqqqZ1TVTYu85T8yHFUtTZSBoHUhyVFJvtDM2npukiOa9Q9O8nfNDVDeOP/o\noqrOA25ZYHdPZs+sozuAf79Ae/dsjh4uS/KXwD3nvfaPSe6bZP/maOHiJJcmeU6SVwH3B85L8oWV\n+wSkpRkIWi/+CPhADWdtPbN5DvAO4H9U1aMZToS36NFFkkOAG6pq7p4c1zD8At/by4Fbq+oRDKcU\neOy81+baOAG4tqoeU1WPAv6qqv6I4QRks1X1lFF/SGkcBoLWi8ezZ7bQPwWOn7f+Y83ymSvY3hOB\nDwFU1U6G8+rPmZsCfSfw1CRvTvKEGs6sOff6apsmXWuAgaD1ok2/wpJfwlV1PXBQkrm/naWnFN7H\nvqvqWwxn790JvClTdr9qrT8Ggtaihb7Y/5bhLShh2FH8183y3wHPbpafv4997b2/84DnNMunAp9c\nYLvzgRcCJPnXwKN+bsfJ/YCfVNWZwB8wDAeAm4B9XoUkdcXZTrXmJLkT+C7DL/IC3s6wE/iDwCHA\nD4AXV9U1zd3TPsSw0/ds4JSqOrLZz/kMb6x+b4b3I3hJVZ2b5EEM7797MMM7aP1WDe/JO7+GezK8\n/8WjgMsZ3rrwFVV1UZLvMOxTeCzDINjNcGr3lzevvxJ4JcP+BfsRNDEGgta1JFvm7jOQ5HnA86vq\n564aktaD3m+QI/Xsl5O8i+HRxA3Mu12itN54hCBJAuxUliQ1DARJEmAgSJIaBoIkCTAQJEkNA0GS\nBMD/B44GJa2m0DgQAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fcf0ef54390>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "219\n"
     ]
    }
   ],
   "source": [
    "exp = Experiment(Pool(topics_filter=OptimizationTopicsFilter(eps=0.55, verbose=False), \n",
    "                      save_topics=True))\n",
    "for i in xrange(5):\n",
    "    model_artm = create_model()\n",
    "    model_artm.fit_offline(batch_vectorizer=batch_vectorizer, num_collection_passes=15, num_document_passes=1)\n",
    "    exp.collect_topics(model_artm.get_phi(), model_artm.get_theta())\n",
    "    vals, bins = exp.topics_pool.topics_filter.plot_hist()\n",
    "    print exp.topics_pool.get_basic_topics_count()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "В конце выведем отфильтрованные темы."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "topic249 | [u'theme' u'strongly' u'grassroots' u'counties' u'corps']\n",
      "topic209 | [u'feb' u'research' u'rep' u'makes' u'means']\n",
      "topic64 | [u'industrial' u'settled' u'passion' u'nod' u'soros']\n",
      "topic167 | [u'rnc' u'smith' u'traffic' u'oregon' u'actual']\n",
      "topic182 | [u'murkowski' u'keyes' u'barbara' u'tony' u'lisa']\n",
      "topic3 | [u'employee' u'indictment' u'teachers' u'austin' u'mcgreevey']\n",
      "topic14 | [u'citizen' u'matches' u'egon' u'casts' u'kerr']\n",
      "topic67 | [u'cbs' u'weekend' u'combat' u'duty' u'served']\n",
      "topic7 | [u'operations' u'coalition' u'insurgents' u'command' u'combat']\n",
      "topic199 | [u'language' u'ideas' u'moores' u'substantial' u'providing']\n",
      "topic233 | [u'handling' u'hughes' u'divided' u'karen' u'sees']\n",
      "topic123 | [u'industrial' u'feelings' u'betting' u'salvador' u'banner']\n",
      "topic112 | [u'keever' u'seemann' u'map' u'gerlach' u'rowland']\n",
      "topic169 | [u'threatened' u'berg' u'factor' u'lists' u'letters']\n",
      "topic59 | [u'club' u'murkowski' u'lisa' u'indian' u'slot']\n",
      "topic127 | [u'opinion' u'favorable' u'unfavorable' u'slightly' u'masondixon']\n",
      "topic175 | [u'pdf' u'congressman' u'surveyusa' u'winner' u'thune']\n",
      "topic170 | [u'rep' u'toomey' u'specter' u'advantage' u'dccc']\n",
      "topic150 | [u'womens' u'fraud' u'conclusions' u'poverty' u'pdb']\n",
      "topic228 | [u'delegate' u'dreier' u'snow' u'mcmahon' u'focusing']\n",
      "topic86 | [u'forecast' u'economist' u'economists' u'wall' u'rates']\n",
      "topic142 | [u'conway' u'myers' u'saleh' u'fighters' u'falluja']\n",
      "topic88 | [u'libby' u'probe' u'plame' u'leak' u'thune']\n",
      "topic56 | [u'cheneys' u'smith' u'walmart' u'hometown' u'mercenary']\n",
      "topic13 | [u'tulsa' u'faster' u'sample' u'registration' u'tightening']\n",
      "topic134 | [u'decision' u'raised' u'thought' u'makes' u'rep']\n",
      "topic243 | [u'size' u'chamber' u'margaret' u'qualified' u'commerce']\n",
      "topic222 | [u'hospital' u'hersh' u'clintons' u'powell' u'surgery']\n",
      "topic4 | [u'engage' u'hispanics' u'oreilly' u'commerce' u'rapids']\n",
      "topic166 | [u'electricity' u'buying' u'emissions' u'motion' u'warming']\n",
      "topic193 | [u'chris' u'catholic' u'davis' u'mydd' u'catholics']\n",
      "topic98 | [u'novak' u'falls' u'machines' u'frost' u'columnist']\n",
      "topic25 | [u'redistricting' u'fahrenheit' u'documentary' u'moores' u'luntz']\n",
      "topic145 | [u'georgia' u'north' u'statewide' u'alaska' u'counties']\n",
      "topic181 | [u'secondplace' u'cent' u'property' u'flipflopper' u'figured']\n",
      "topic69 | [u'kerryedwards' u'legal' u'thoughts' u'catholic' u'rightwing']\n",
      "topic102 | [u'hurry' u'deserved' u'cells' u'rooms' u'guild']\n",
      "topic235 | [u'face' u'head' u'system' u'americas' u'moment']\n",
      "topic63 | [u'oregon' u'bounce' u'finish' u'unlike' u'electorate']\n",
      "topic132 | [u'veep' u'duty' u'cleland' u'unit' u'discharge']\n",
      "topic43 | [u'occupation' u'council' u'violence' u'prison' u'wounded']\n",
      "topic60 | [u'soft' u'gov' u'latino' u'material' u'kean']\n",
      "topic143 | [u'chalabi' u'daschle' u'diedrich' u'pelosi' u'speaker']\n",
      "topic17 | [u'daschle' u'dakota' u'cash' u'rnc' u'thune']\n",
      "topic187 | [u'commander' u'scientists' u'husseins' u'debt' u'epa']\n",
      "topic31 | [u'seattle' u'chafee' u'retaliation' u'pleasure' u'plea']\n",
      "topic239 | [u'tribune' u'ballots' u'limbaugh' u'martinez' u'vitter']\n",
      "topic61 | [u'wage' u'cap' u'protesters' u'mcgreevey' u'award']\n",
      "topic219 | [u'taxes' u'science' u'edition' u'debates' u'abc']\n",
      "topic202 | [u'percentage' u'arg' u'conducted' u'smith' u'factor']\n",
      "topic20 | [u'sinclair' u'broadcast' u'channel' u'viewers' u'balanced']\n",
      "topic204 | [u'mouth' u'kerr' u'hunt' u'chandlers' u'catholics']\n",
      "topic164 | [u'couples' u'sharon' u'palestinian' u'adults' u'samesex']\n",
      "topic241 | [u'australia' u'cincinnati' u'hospitals' u'rampant' u'bounds']\n",
      "topic148 | [u'grand' u'inquiry' u'criminal' u'investigating' u'plame']\n",
      "topic5 | [u'invasion' u'postwar' u'uniform' u'taliban' u'safer']\n",
      "topic66 | [u'trial' u'stations' u'claimed' u'declared' u'host']\n",
      "topic54 | [u'sounds' u'highly' u'stand' u'anger' u'providing']\n",
      "topic151 | [u'industry' u'countries' u'nuclear' u'explosives' u'donald']\n",
      "topic1 | [u'english' u'muslims' u'wink' u'wmds' u'bombs']\n",
      "topic213 | [u'moveon' u'kean' u'panel' u'condi' u'rices']\n",
      "topic104 | [u'morrison' u'forget' u'matching' u'cash' u'club']\n",
      "topic237 | [u'guard' u'records' u'duty' u'served' u'active']\n",
      "topic157 | [u'chandler' u'ben' u'bid' u'chicago' u'kerr']\n",
      "topic52 | [u'swift' u'star' u'veterans' u'boat' u'scott']\n",
      "topic6 | [u'pdf' u'surveyusa' u'louisiana' u'tied' u'usa']\n",
      "topic41 | [u'murphy' u'musgrave' u'diedrich' u'matsunaka' u'larry']\n",
      "topic176 | [u'dakota' u'diedrich' u'brad' u'pelosi' u'stephanie']\n",
      "topic124 | [u'missing' u'invasion' u'afraid' u'flight' u'agreed']\n",
      "topic205 | [u'pollution' u'park' u'possibility' u'mydd' u'emissions']\n",
      "topic109 | [u'carry' u'caught' u'targetting' u'firefighters' u'surface']\n",
      "topic224 | [u'deaths' u'occupation' u'wounded' u'chalabi' u'religious']\n",
      "topic106 | [u'boxer' u'dreier' u'schwarzenegger' u'reserves' u'adjust']\n",
      "topic55 | [u'pdf' u'college' u'survey' u'percentage' u'arg']\n",
      "topic153 | [u'geneva' u'global' u'karen' u'liberty' u'burden']\n",
      "topic74 | [u'blogger' u'indicted' u'earle' u'billmon' u'features']\n",
      "topic229 | [u'nrcc' u'murphy' u'usa' u'posting' u'rowland']\n",
      "topic12 | [u'prosecutors' u'rowland' u'fraud' u'subpoena' u'justices']\n",
      "topic38 | [u'signatures' u'jeff' u'pelosi' u'lois' u'filing']\n",
      "topic68 | [u'womens' u'saudi' u'arabia' u'planet' u'epa']\n",
      "topic94 | [u'congressman' u'districts' u'contest' u'congressional' u'helping']\n",
      "topic77 | [u'deal' u'jun' u'reporters' u'efforts' u'chief']\n",
      "topic34 | [u'hints' u'unconventional' u'postconvention' u'playing' u'pleasure']\n",
      "topic126 | [u'memo' u'fbi' u'medicare' u'hoeffel' u'request']\n",
      "topic83 | [u'reform' u'houston' u'debt' u'relief' u'nra']\n",
      "topic217 | [u'nuclear' u'miles' u'delays' u'body' u'pelosi']\n",
      "topic92 | [u'black' u'werent' u'islam' u'factions' u'ideals']\n",
      "topic186 | [u'error' u'blogs' u'arg' u'baseball' u'straw']\n",
      "topic218 | [u'armed' u'accomplished' u'private' u'contractors' u'husseins']\n",
      "topic190 | [u'farmer' u'chedrcheez' u'murphy' u'king' u'stork']\n",
      "topic85 | [u'chalabis' u'berg' u'mosul' u'rewards' u'contractor']\n",
      "topic15 | [u'annoying' u'medals' u'stripped' u'blitzer' u'titan']\n",
      "topic203 | [u'detroit' u'islam' u'richardson' u'vance' u'roosevelt']\n",
      "topic141 | [u'furious' u'treasury' u'rich' u'oneill' u'unemployment']\n",
      "topic105 | [u'josh' u'living' u'threatened' u'scientists' u'species']\n",
      "topic116 | [u'organizations' u'rosenberg' u'simon' u'ndn' u'sees']\n",
      "topic29 | [u'debt' u'meetup' u'born' u'roots' u'damning']\n",
      "topic240 | [u'mars' u'rogers' u'reasonable' u'faster' u'species']\n",
      "topic23 | [u'science' u'cell' u'scientists' u'institute' u'britain']\n",
      "topic76 | [u'rankings' u'cattle' u'pelosi' u'buzz' u'diedrichs']\n",
      "topic144 | [u'clintons' u'ban' u'risk' u'surgery' u'cards']\n",
      "topic238 | [u'register' u'adults' u'uncommitted' u'sproul' u'junkies']\n",
      "topic133 | [u'disaster' u'briefing' u'clarkes' u'strike' u'iran']\n",
      "topic16 | [u'apologies' u'sisters' u'plea' u'plays' u'pledge']\n",
      "topic191 | [u'georgia' u'activists' u'cleland' u'martinez' u'iowans']\n",
      "topic198 | [u'rowland' u'hart' u'nagourney' u'greenwood' u'adam']\n",
      "topic99 | [u'nrsc' u'margins' u'senatorial' u'rosenberg' u'hart']\n",
      "topic27 | [u'cattle' u'vacation' u'strike' u'enron' u'industry']\n",
      "topic221 | [u'veep' u'signatures' u'naders' u'nelson' u'arizona']\n",
      "topic165 | [u'ndn' u'caucus' u'establishment' u'ideological' u'dlc']\n",
      "topic168 | [u'listen' u'body' u'musgrave' u'rogers' u'outraged']\n",
      "topic28 | [u'papers' u'cheneys' u'activities' u'editorial' u'awol']\n",
      "topic223 | [u'guest' u'blogger' u'yellin' u'ilsen' u'frost']\n",
      "topic117 | [u'scientific' u'sales' u'enron' u'emissions' u'manufacturing']\n",
      "topic161 | [u'shares' u'constitutional' u'advertisers' u'patriot' u'partners']\n",
      "topic210 | [u'stork' u'matching' u'stan' u'clerk' u'hardest']\n",
      "topic207 | [u'oregon' u'losses' u'cleveland' u'households' u'manufacturing']\n",
      "topic101 | [u'graham' u'posting' u'face' u'announcement' u'cheneys']\n",
      "topic231 | [u'advice' u'fired' u'abuse' u'fraud' u'paper']\n",
      "topic220 | [u'traffic' u'servers' u'server' u'nominate' u'choose']\n",
      "topic58 | [u'arab' u'armed' u'gen' u'fallujah' u'commander']\n",
      "topic156 | [u'rice' u'adviser' u'powell' u'abc' u'nbc']\n",
      "topic147 | [u'caucus' u'gallup' u'braun' u'caucuses' u'clarks']\n",
      "topic40 | [u'delaware' u'utah' u'bond' u'mccain' u'thurlow']\n",
      "topic32 | [u'testimony' u'clarke' u'wilson' u'journalists' u'plame']\n",
      "topic118 | [u'veterans' u'faith' u'rally' u'hughes' u'pool']\n",
      "topic246 | [u'jones' u'tnr' u'publication' u'sweet' u'cent']\n",
      "topic247 | [u'hoeffel' u'brain' u'guilty' u'specters' u'prosecutors']\n",
      "topic174 | [u'alabama' u'splice' u'fulfill' u'theaters' u'exam']\n",
      "topic137 | [u'tariffs' u'steel' u'korea' u'europeans' u'frost']\n",
      "topic89 | [u'wingnut' u'qpoll' u'requests' u'sterns' u'toomeys']\n",
      "topic115 | [u'geneva' u'supports' u'constitutional' u'opposes' u'banning']\n",
      "topic114 | [u'green' u'nevada' u'districts' u'dccc' u'schrader']\n",
      "topic130 | [u'keyes' u'alan' u'lightbulb' u'weirdness' u'enters']\n",
      "topic215 | [u'enron' u'resolution' u'bell' u'inflation' u'french']\n",
      "topic24 | [u'ballots' u'stupid' u'reform' u'beauprez' u'absentee']\n",
      "topic236 | [u'human' u'rights' u'abortion' u'citizens' u'laws']\n",
      "topic158 | [u'coors' u'schaffer' u'salazar' u'soros' u'birthday']\n",
      "topic9 | [u'divorce' u'chicago' u'prices' u'boxer' u'flu']\n",
      "topic185 | [u'unemployment' u'india' u'crisis' u'estimate' u'economists']\n",
      "topic211 | [u'company' u'fbi' u'cheneys' u'request' u'ties']\n",
      "topic138 | [u'rumsfelds' u'weapon' u'donald' u'guns' u'sharon']\n",
      "topic22 | [u'latino' u'ndn' u'legislation' u'child' u'proposed']\n",
      "topic208 | [u'journalists' u'liberals' u'sanctions' u'blood' u'industry']\n",
      "topic18 | [u'explosives' u'traffic' u'flip' u'examples' u'missing']\n",
      "topic73 | [u'ethics' u'chandler' u'dakota' u'governors' u'ken']\n",
      "topic171 | [u'british' u'combat' u'resistance' u'torture' u'armed']\n",
      "topic42 | [u'tenenbaum' u'demint' u'martinez' u'inez' u'nebraska']\n",
      "topic183 | [u'stem' u'schwarzenegger' u'machines' u'dreier' u'discipline']\n",
      "topic21 | [u'surgery' u'coffee' u'humphreys' u'coburns' u'hospital']\n",
      "topic45 | [u'condis' u'declassified' u'connect' u'apt' u'beers']\n",
      "topic19 | [u'means' u'dkos' u'rep' u'talking' u'advantage']\n",
      "topic65 | [u'tariffs' u'consumer' u'trippi' u'europeans' u'steel']\n",
      "topic173 | [u'afscme' u'trippi' u'seiu' u'aflcio' u'prediction']\n",
      "topic125 | [u'rasmussen' u'traffic' u'means' u'theme' u'error']\n",
      "topic113 | [u'servers' u'unfit' u'creative' u'implement' u'flowing']\n",
      "topic10 | [u'nag' u'funds' u'matching' u'oneill' u'inflation']\n",
      "topic84 | [u'dkos' u'research' u'difference' u'page' u'thousands']\n",
      "topic57 | [u'fiscal' u'destruction' u'nuclear' u'donald' u'choices']\n",
      "topic230 | [u'deal' u'failure' u'osama' u'saudi' u'rnc']\n",
      "topic93 | [u'postdebate' u'mad' u'outreach' u'steve' u'indictments']\n",
      "topic26 | [u'mars' u'limbaugh' u'reserves' u'doctor' u'tobacco']\n",
      "topic62 | [u'matsunaka' u'afscme' u'seiu' u'nag' u'boy']\n",
      "topic44 | [u'income' u'cap' u'definition' u'unfit' u'novak']\n",
      "topic146 | [u'resistance' u'shiite' u'najaf' u'film' u'british']\n",
      "topic35 | [u'bugs' u'calm' u'npr' u'rusty' u'located']\n",
      "topic179 | [u'responding' u'defensive' u'accusations' u'mcgreevey' u'oconnor']\n",
      "topic163 | [u'medical' u'reconstruction' u'missile' u'coffee' u'mars']\n",
      "topic90 | [u'bunning' u'caucuses' u'mongiardo' u'manager' u'volunteers']\n",
      "topic50 | [u'laden' u'film' u'muslim' u'osama' u'juan']\n",
      "topic189 | [u'headtohead' u'conservatism' u'buys' u'dreams' u'compassionate']\n",
      "topic216 | [u'indiana' u'fletcher' u'daniels' u'montana' u'inactive']\n",
      "topic178 | [u'invasion' u'black' u'car' u'child' u'disaster']\n",
      "topic103 | [u'humphreys' u'tulsa' u'burt' u'parity' u'roosevelt']\n",
      "topic80 | [u'upenn' u'drug' u'drugs' u'apologies' u'canadian']\n",
      "topic129 | [u'angry' u'anger' u'moderates' u'ceo' u'thurlow']\n",
      "topic111 | [u'courts' u'divorce' u'houston' u'ruled' u'justices']\n",
      "topic244 | [u'approve' u'braun' u'parenthesis' u'disapprove' u'buzz']\n",
      "topic201 | [u'incidentally' u'indian' u'johnson' u'native' u'thune']\n",
      "topic30 | [u'candidacy' u'momentum' u'bid' u'opportunities' u'flag']\n",
      "topic227 | [u'primaries' u'wilson' u'map' u'dccc' u'schrader']\n",
      "topic226 | [u'reynolds' u'demint' u'contribution' u'firefighters' u'norquist']\n",
      "topic96 | [u'lawyers' u'experts' u'justice' u'crisis' u'medical']\n",
      "topic212 | [u'caucus' u'forces' u'takes' u'west' u'mike']\n",
      "topic232 | [u'lisa' u'contested' u'targeting' u'ideas' u'falls']\n",
      "topic128 | [u'matsunaka' u'nancy' u'stations' u'catholic' u'liberals']\n",
      "topic0 | [u'images' u'chain' u'stealing' u'frost' u'dallas']\n",
      "topic75 | [u'required' u'filed' u'map' u'registrations' u'cuba']\n",
      "topic72 | [u'harris' u'gas' u'climate' u'guantanamo' u'legislation']\n",
      "topic95 | [u'geneva' u'resolution' u'beasley' u'questioning' u'impeachment']\n",
      "topic46 | [u'wall' u'nuclear' u'donald' u'mcclellan' u'missile']\n",
      "topic195 | [u'morrison' u'stations' u'unfit' u'blogger' u'konop']\n",
      "topic37 | [u'spains' u'alexander' u'lindauer' u'wage' u'rodney']\n",
      "topic162 | [u'dollar' u'revenue' u'ann' u'deficit' u'provided']\n",
      "topic108 | [u'richard' u'forces' u'talking' u'coalition' u'hussein']\n",
      "topic184 | [u'con' u'slight' u'idetestthispres' u'reverse' u'hotshotxi']\n",
      "topic159 | [u'shiite' u'najaf' u'council' u'interim' u'bremer']\n",
      "topic192 | [u'surgery' u'leans' u'shots' u'consistently' u'flu']\n",
      "topic121 | [u'income' u'fatalities' u'missile' u'arms' u'cap']\n",
      "topic53 | [u'registrants' u'floridas' u'census' u'shells' u'beach']\n",
      "topic135 | [u'trippi' u'moines' u'des' u'volunteers' u'rankings']\n",
      "topic196 | [u'kerryedwards' u'events' u'bushcheney' u'editors' u'harris']\n",
      "topic51 | [u'limbaugh' u'sinclair' u'fantasy' u'flu' u'vaccine']\n",
      "topic33 | [u'toomey' u'bunning' u'mongiardo' u'hoeffel' u'konop']\n",
      "topic160 | [u'saudi' u'boat' u'swift' u'arabia' u'relations']\n",
      "topic242 | [u'sinclair' u'lawyer' u'lawyers' u'contracts' u'guantanamo']\n",
      "topic194 | [u'pac' u'enron' u'wage' u'indicted' u'corporations']\n",
      "topic110 | [u'approve' u'cbs' u'kerryedwards' u'bushcheney' u'pat']\n",
      "topic119 | [u'flip' u'bond' u'bunnings' u'economists' u'agreement']\n",
      "topic48 | [u'harris' u'cuba' u'schneider' u'enlarge' u'floridas']\n",
      "topic149 | [u'leans' u'lisa' u'pro' u'ditka' u'river']\n",
      "topic71 | [u'alexander' u'bean' u'schrader' u'dccc' u'keever']\n",
      "topic140 | [u'hall' u'gas' u'margins' u'wilson' u'eve']\n",
      "topic136 | [u'explosives' u'takes' u'abuses' u'land' u'womens']\n",
      "topic36 | [u'reagans' u'ronald' u'catholic' u'relations' u'burden']\n",
      "topic248 | [u'corner' u'kids' u'amp' u'husband' u'proudly']\n",
      "topic206 | [u'blair' u'cuba' u'assholes' u'clerks' u'telegraph']\n",
      "topic78 | [u'ryan' u'allegations' u'hoeffel' u'lying' u'appointed']\n",
      "topic79 | [u'clarke' u'richard' u'contributions' u'morrison' u'wilson']\n"
     ]
    }
   ],
   "source": [
    "exp.show_basic_topics()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Также иногда полезно посмотреть на все темы, которые использовались в эксперименте."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "topic0 | [u'images' u'chain' u'stealing' u'frost' u'dallas'] | topic0\n",
      "topic1 | [u'english' u'muslims' u'wink' u'wmds' u'bombs'] | topic1\n",
      "topic10 | [u'nag' u'funds' u'matching' u'oneill' u'inflation'] | topic10\n",
      "topic101 | [u'graham' u'posting' u'face' u'announcement' u'cheneys'] | topic101\n",
      "topic102 | [u'hurry' u'deserved' u'cells' u'rooms' u'guild'] | topic102\n",
      "topic103 | [u'humphreys' u'tulsa' u'burt' u'parity' u'roosevelt'] | topic103\n",
      "topic104 | [u'morrison' u'forget' u'matching' u'cash' u'club'] | topic104\n",
      "topic105 | [u'josh' u'living' u'threatened' u'scientists' u'species'] | topic105\n",
      "topic106 | [u'boxer' u'dreier' u'schwarzenegger' u'reserves' u'adjust'] | topic106\n",
      "topic108 | [u'richard' u'forces' u'talking' u'coalition' u'hussein'] | topic108\n",
      "topic109 | [u'carry' u'caught' u'targetting' u'firefighters' u'surface'] | topic109\n",
      "topic110 | [u'approve' u'cbs' u'kerryedwards' u'bushcheney' u'pat'] | topic110\n",
      "topic111 | [u'courts' u'divorce' u'houston' u'ruled' u'justices'] | topic111\n",
      "topic112 | [u'keever' u'seemann' u'map' u'gerlach' u'rowland'] | topic112\n",
      "topic113 | [u'servers' u'unfit' u'creative' u'implement' u'flowing'] | topic113\n",
      "topic114 | [u'green' u'nevada' u'districts' u'dccc' u'schrader'] | topic114\n",
      "topic115 | [u'geneva' u'supports' u'constitutional' u'opposes' u'banning'] | topic115\n",
      "topic116 | [u'organizations' u'rosenberg' u'simon' u'ndn' u'sees'] | topic116\n",
      "topic117 | [u'scientific' u'sales' u'enron' u'emissions' u'manufacturing'] | topic117\n",
      "topic118 | [u'veterans' u'faith' u'rally' u'hughes' u'pool'] | topic118\n",
      "topic119 | [u'flip' u'bond' u'bunnings' u'economists' u'agreement'] | topic119\n",
      "topic12 | [u'prosecutors' u'rowland' u'fraud' u'subpoena' u'justices'] | topic12\n",
      "topic121 | [u'income' u'fatalities' u'missile' u'arms' u'cap'] | topic121\n",
      "topic123 | [u'industrial' u'feelings' u'betting' u'salvador' u'banner'] | topic123\n",
      "topic124 | [u'missing' u'invasion' u'afraid' u'flight' u'agreed'] | topic124\n",
      "topic125 | [u'rasmussen' u'traffic' u'means' u'theme' u'error'] | topic125\n",
      "topic126 | [u'memo' u'fbi' u'medicare' u'hoeffel' u'request'] | topic126\n",
      "topic127 | [u'opinion' u'favorable' u'unfavorable' u'slightly' u'masondixon'] | topic127\n",
      "topic128 | [u'matsunaka' u'nancy' u'stations' u'catholic' u'liberals'] | topic128\n",
      "topic129 | [u'angry' u'anger' u'moderates' u'ceo' u'thurlow'] | topic129\n",
      "topic13 | [u'tulsa' u'faster' u'sample' u'registration' u'tightening'] | topic13\n",
      "topic155 | [u'totals' u'registrations' u'registration' u'registrants' u'las'] | topic13\n",
      "topic130 | [u'keyes' u'alan' u'lightbulb' u'weirdness' u'enters'] | topic130\n",
      "topic132 | [u'veep' u'duty' u'cleland' u'unit' u'discharge'] | topic132\n",
      "topic133 | [u'disaster' u'briefing' u'clarkes' u'strike' u'iran'] | topic133\n",
      "topic2 | [u'ryan' u'salazar' u'keyes' u'ken' u'jack'] | topic134\n",
      "topic8 | [u'speech' u'air' u'terrorism' u'clinton' u'remember'] | topic134\n",
      "topic39 | [u'christopher' u'mike' u'chemical' u'profile' u'scared'] | topic134\n",
      "topic49 | [u'takes' u'congressional' u'democracy' u'saudi' u'previously'] | topic134\n",
      "topic81 | [u'mike' u'scared' u'releasing' u'mich' u'pictures'] | topic134\n",
      "topic82 | [u'keyes' u'briefing' u'interviews' u'explosives' u'leans'] | topic134\n",
      "topic122 | [u'oregon' u'naders' u'signatures' u'petition' u'registrations'] | topic134\n",
      "topic131 | [u'profile' u'christopher' u'mike' u'scared' u'attempt'] | topic134\n",
      "topic134 | [u'decision' u'raised' u'thought' u'makes' u'rep'] | topic134\n",
      "topic154 | [u'remember' u'reason' u'gore' u'feb' u'line'] | topic134\n",
      "topic172 | [u'chemical' u'crooks' u'falling' u'virginiadem' u'bloomfield'] | topic134\n",
      "topic197 | [u'prison' u'clarke' u'testimony' u'criminal' u'abuse'] | topic134\n",
      "topic245 | [u'salazar' u'matsunaka' u'ken' u'musgrave' u'leans'] | topic134\n",
      "topic135 | [u'trippi' u'moines' u'des' u'volunteers' u'rankings'] | topic135\n",
      "topic136 | [u'explosives' u'takes' u'abuses' u'land' u'womens'] | topic136\n",
      "topic137 | [u'tariffs' u'steel' u'korea' u'europeans' u'frost'] | topic137\n",
      "topic138 | [u'rumsfelds' u'weapon' u'donald' u'guns' u'sharon'] | topic138\n",
      "topic14 | [u'citizen' u'matches' u'egon' u'casts' u'kerr'] | topic14\n",
      "topic140 | [u'hall' u'gas' u'margins' u'wilson' u'eve'] | topic140\n",
      "topic141 | [u'furious' u'treasury' u'rich' u'oneill' u'unemployment'] | topic141\n",
      "topic142 | [u'conway' u'myers' u'saleh' u'fighters' u'falluja'] | topic142\n",
      "topic143 | [u'chalabi' u'daschle' u'diedrich' u'pelosi' u'speaker'] | topic143\n",
      "topic144 | [u'clintons' u'ban' u'risk' u'surgery' u'cards'] | topic144\n",
      "topic145 | [u'georgia' u'north' u'statewide' u'alaska' u'counties'] | topic145\n",
      "topic146 | [u'resistance' u'shiite' u'najaf' u'film' u'british'] | topic146\n",
      "topic11 | [u'delegates' u'caucus' u'braun' u'arg' u'caucuses'] | topic147\n",
      "topic147 | [u'caucus' u'gallup' u'braun' u'caucuses' u'clarks'] | topic147\n",
      "topic148 | [u'grand' u'inquiry' u'criminal' u'investigating' u'plame'] | topic148\n",
      "topic149 | [u'leans' u'lisa' u'pro' u'ditka' u'river'] | topic149\n",
      "topic15 | [u'annoying' u'medals' u'stripped' u'blitzer' u'titan'] | topic15\n",
      "topic150 | [u'womens' u'fraud' u'conclusions' u'poverty' u'pdb'] | topic150\n",
      "topic151 | [u'industry' u'countries' u'nuclear' u'explosives' u'donald'] | topic151\n",
      "topic153 | [u'geneva' u'global' u'karen' u'liberty' u'burden'] | topic153\n",
      "topic91 | [u'powell' u'adviser' u'sept' u'rice' u'abc'] | topic156\n",
      "topic100 | [u'edition' u'nbc' u'baseball' u'lineup' u'sens'] | topic156\n",
      "topic156 | [u'rice' u'adviser' u'powell' u'abc' u'nbc'] | topic156\n",
      "topic157 | [u'chandler' u'ben' u'bid' u'chicago' u'kerr'] | topic157\n",
      "topic158 | [u'coors' u'schaffer' u'salazar' u'soros' u'birthday'] | topic158\n",
      "topic159 | [u'shiite' u'najaf' u'council' u'interim' u'bremer'] | topic159\n",
      "topic16 | [u'apologies' u'sisters' u'plea' u'plays' u'pledge'] | topic16\n",
      "topic160 | [u'saudi' u'boat' u'swift' u'arabia' u'relations'] | topic160\n",
      "topic161 | [u'shares' u'constitutional' u'advertisers' u'patriot' u'partners'] | topic161\n",
      "topic162 | [u'dollar' u'revenue' u'ann' u'deficit' u'provided'] | topic162\n",
      "topic163 | [u'medical' u'reconstruction' u'missile' u'coffee' u'mars'] | topic163\n",
      "topic164 | [u'couples' u'sharon' u'palestinian' u'adults' u'samesex'] | topic164\n",
      "topic165 | [u'ndn' u'caucus' u'establishment' u'ideological' u'dlc'] | topic165\n",
      "topic166 | [u'electricity' u'buying' u'emissions' u'motion' u'warming'] | topic166\n",
      "topic167 | [u'rnc' u'smith' u'traffic' u'oregon' u'actual'] | topic167\n",
      "topic168 | [u'listen' u'body' u'musgrave' u'rogers' u'outraged'] | topic168\n",
      "topic169 | [u'threatened' u'berg' u'factor' u'lists' u'letters'] | topic169\n",
      "topic17 | [u'daschle' u'dakota' u'cash' u'rnc' u'thune'] | topic17\n",
      "topic170 | [u'rep' u'toomey' u'specter' u'advantage' u'dccc'] | topic170\n",
      "topic214 | [u'mcauliffe' u'alexander' u'hastert' u'humphreys' u'rodney'] | topic170\n",
      "topic171 | [u'british' u'combat' u'resistance' u'torture' u'armed'] | topic171\n",
      "topic173 | [u'afscme' u'trippi' u'seiu' u'aflcio' u'prediction'] | topic173\n",
      "topic225 | [u'afscme' u'seiu' u'industrial' u'channel' u'coburns'] | topic173\n",
      "topic174 | [u'alabama' u'splice' u'fulfill' u'theaters' u'exam'] | topic174\n",
      "topic175 | [u'pdf' u'congressman' u'surveyusa' u'winner' u'thune'] | topic175\n",
      "topic176 | [u'dakota' u'diedrich' u'brad' u'pelosi' u'stephanie'] | topic176\n",
      "topic178 | [u'invasion' u'black' u'car' u'child' u'disaster'] | topic178\n",
      "topic179 | [u'responding' u'defensive' u'accusations' u'mcgreevey' u'oconnor'] | topic179\n",
      "topic18 | [u'explosives' u'traffic' u'flip' u'examples' u'missing'] | topic18\n",
      "topic181 | [u'secondplace' u'cent' u'property' u'flipflopper' u'figured'] | topic181\n",
      "topic182 | [u'murkowski' u'keyes' u'barbara' u'tony' u'lisa'] | topic182\n",
      "topic183 | [u'stem' u'schwarzenegger' u'machines' u'dreier' u'discipline'] | topic183\n",
      "topic97 | [u'moment' u'watch' u'slight' u'kids' u'atlanta'] | topic184\n",
      "topic107 | [u'slight' u'reverse' u'kids' u'stocks' u'atlanta'] | topic184\n",
      "topic184 | [u'con' u'slight' u'idetestthispres' u'reverse' u'hotshotxi'] | topic184\n",
      "topic200 | [u'mich' u'falling' u'bloomfield' u'crooks' u'virginiadem'] | topic184\n",
      "topic185 | [u'unemployment' u'india' u'crisis' u'estimate' u'economists'] | topic185\n",
      "topic186 | [u'error' u'blogs' u'arg' u'baseball' u'straw'] | topic186\n",
      "topic187 | [u'commander' u'scientists' u'husseins' u'debt' u'epa'] | topic187\n",
      "topic189 | [u'headtohead' u'conservatism' u'buys' u'dreams' u'compassionate'] | topic189\n",
      "topic19 | [u'means' u'dkos' u'rep' u'talking' u'advantage'] | topic19\n",
      "topic190 | [u'farmer' u'chedrcheez' u'murphy' u'king' u'stork'] | topic190\n",
      "topic191 | [u'georgia' u'activists' u'cleland' u'martinez' u'iowans'] | topic191\n",
      "topic192 | [u'surgery' u'leans' u'shots' u'consistently' u'flu'] | topic192\n",
      "topic193 | [u'chris' u'catholic' u'davis' u'mydd' u'catholics'] | topic193\n",
      "topic194 | [u'pac' u'enron' u'wage' u'indicted' u'corporations'] | topic194\n",
      "topic195 | [u'morrison' u'stations' u'unfit' u'blogger' u'konop'] | topic195\n",
      "topic196 | [u'kerryedwards' u'events' u'bushcheney' u'editors' u'harris'] | topic196\n",
      "topic198 | [u'rowland' u'hart' u'nagourney' u'greenwood' u'adam'] | topic198\n",
      "topic199 | [u'language' u'ideas' u'moores' u'substantial' u'providing'] | topic199\n",
      "topic20 | [u'sinclair' u'broadcast' u'channel' u'viewers' u'balanced'] | topic20\n",
      "topic201 | [u'incidentally' u'indian' u'johnson' u'native' u'thune'] | topic201\n",
      "topic202 | [u'percentage' u'arg' u'conducted' u'smith' u'factor'] | topic202\n",
      "topic203 | [u'detroit' u'islam' u'richardson' u'vance' u'roosevelt'] | topic203\n",
      "topic204 | [u'mouth' u'kerr' u'hunt' u'chandlers' u'catholics'] | topic204\n",
      "topic205 | [u'pollution' u'park' u'possibility' u'mydd' u'emissions'] | topic205\n",
      "topic177 | [u'blair' u'perry' u'land' u'authorization' u'feith'] | topic206\n",
      "topic206 | [u'blair' u'cuba' u'assholes' u'clerks' u'telegraph'] | topic206\n",
      "topic207 | [u'oregon' u'losses' u'cleveland' u'households' u'manufacturing'] | topic207\n",
      "topic208 | [u'journalists' u'liberals' u'sanctions' u'blood' u'industry'] | topic208\n",
      "topic209 | [u'feb' u'research' u'rep' u'makes' u'means'] | topic209\n",
      "topic21 | [u'surgery' u'coffee' u'humphreys' u'coburns' u'hospital'] | topic21\n",
      "topic210 | [u'stork' u'matching' u'stan' u'clerk' u'hardest'] | topic210\n",
      "topic211 | [u'company' u'fbi' u'cheneys' u'request' u'ties'] | topic211\n",
      "topic188 | [u'gov' u'takes' u'stupid' u'undecideds' u'rightwing'] | topic212\n",
      "topic212 | [u'caucus' u'forces' u'takes' u'west' u'mike'] | topic212\n",
      "topic213 | [u'moveon' u'kean' u'panel' u'condi' u'rices'] | topic213\n",
      "topic215 | [u'enron' u'resolution' u'bell' u'inflation' u'french'] | topic215\n",
      "topic216 | [u'indiana' u'fletcher' u'daniels' u'montana' u'inactive'] | topic216\n",
      "topic217 | [u'nuclear' u'miles' u'delays' u'body' u'pelosi'] | topic217\n",
      "topic218 | [u'armed' u'accomplished' u'private' u'contractors' u'husseins'] | topic218\n",
      "topic219 | [u'taxes' u'science' u'edition' u'debates' u'abc'] | topic219\n",
      "topic22 | [u'latino' u'ndn' u'legislation' u'child' u'proposed'] | topic22\n",
      "topic220 | [u'traffic' u'servers' u'server' u'nominate' u'choose'] | topic220\n",
      "topic70 | [u'naders' u'signatures' u'schwarzenegger' u'boxer' u'dreier'] | topic221\n",
      "topic221 | [u'veep' u'signatures' u'naders' u'nelson' u'arizona'] | topic221\n",
      "topic222 | [u'hospital' u'hersh' u'clintons' u'powell' u'surgery'] | topic222\n",
      "topic223 | [u'guest' u'blogger' u'yellin' u'ilsen' u'frost'] | topic223\n",
      "topic224 | [u'deaths' u'occupation' u'wounded' u'chalabi' u'religious'] | topic224\n",
      "topic226 | [u'reynolds' u'demint' u'contribution' u'firefighters' u'norquist'] | topic226\n",
      "topic227 | [u'primaries' u'wilson' u'map' u'dccc' u'schrader'] | topic227\n",
      "topic228 | [u'delegate' u'dreier' u'snow' u'mcmahon' u'focusing'] | topic228\n",
      "topic229 | [u'nrcc' u'murphy' u'usa' u'posting' u'rowland'] | topic229\n",
      "topic23 | [u'science' u'cell' u'scientists' u'institute' u'britain'] | topic23\n",
      "topic230 | [u'deal' u'failure' u'osama' u'saudi' u'rnc'] | topic230\n",
      "topic231 | [u'advice' u'fired' u'abuse' u'fraud' u'paper'] | topic231\n",
      "topic232 | [u'lisa' u'contested' u'targeting' u'ideas' u'falls'] | topic232\n",
      "topic233 | [u'handling' u'hughes' u'divided' u'karen' u'sees'] | topic233\n",
      "topic235 | [u'face' u'head' u'system' u'americas' u'moment'] | topic235\n",
      "topic236 | [u'human' u'rights' u'abortion' u'citizens' u'laws'] | topic236\n",
      "topic152 | [u'service' u'guard' u'records' u'duty' u'active'] | topic237\n",
      "topic237 | [u'guard' u'records' u'duty' u'served' u'active'] | topic237\n",
      "topic238 | [u'register' u'adults' u'uncommitted' u'sproul' u'junkies'] | topic238\n",
      "topic239 | [u'tribune' u'ballots' u'limbaugh' u'martinez' u'vitter'] | topic239\n",
      "topic24 | [u'ballots' u'stupid' u'reform' u'beauprez' u'absentee'] | topic24\n",
      "topic240 | [u'mars' u'rogers' u'reasonable' u'faster' u'species'] | topic240\n",
      "topic241 | [u'australia' u'cincinnati' u'hospitals' u'rampant' u'bounds'] | topic241\n",
      "topic242 | [u'sinclair' u'lawyer' u'lawyers' u'contracts' u'guantanamo'] | topic242\n",
      "topic243 | [u'size' u'chamber' u'margaret' u'qualified' u'commerce'] | topic243\n",
      "topic244 | [u'approve' u'braun' u'parenthesis' u'disapprove' u'buzz'] | topic244\n",
      "topic246 | [u'jones' u'tnr' u'publication' u'sweet' u'cent'] | topic246\n",
      "topic247 | [u'hoeffel' u'brain' u'guilty' u'specters' u'prosecutors'] | topic247\n",
      "topic87 | [u'husband' u'amp' u'corner' u'mom' u'pride'] | topic248\n",
      "topic248 | [u'corner' u'kids' u'amp' u'husband' u'proudly'] | topic248\n",
      "topic249 | [u'theme' u'strongly' u'grassroots' u'counties' u'corps'] | topic249\n",
      "topic25 | [u'redistricting' u'fahrenheit' u'documentary' u'moores' u'luntz'] | topic25\n",
      "topic26 | [u'mars' u'limbaugh' u'reserves' u'doctor' u'tobacco'] | topic26\n",
      "topic27 | [u'cattle' u'vacation' u'strike' u'enron' u'industry'] | topic27\n",
      "topic28 | [u'papers' u'cheneys' u'activities' u'editorial' u'awol'] | topic28\n",
      "topic29 | [u'debt' u'meetup' u'born' u'roots' u'damning'] | topic29\n",
      "topic3 | [u'employee' u'indictment' u'teachers' u'austin' u'mcgreevey'] | topic3\n",
      "topic30 | [u'candidacy' u'momentum' u'bid' u'opportunities' u'flag'] | topic30\n",
      "topic31 | [u'seattle' u'chafee' u'retaliation' u'pleasure' u'plea'] | topic31\n",
      "topic32 | [u'testimony' u'clarke' u'wilson' u'journalists' u'plame'] | topic32\n",
      "topic33 | [u'toomey' u'bunning' u'mongiardo' u'hoeffel' u'konop'] | topic33\n",
      "topic234 | [u'bunning' u'mongiardo' u'kentucky' u'income' u'bunnings'] | topic33\n",
      "topic34 | [u'hints' u'unconventional' u'postconvention' u'playing' u'pleasure'] | topic34\n",
      "topic35 | [u'bugs' u'calm' u'npr' u'rusty' u'located'] | topic35\n",
      "topic36 | [u'reagans' u'ronald' u'catholic' u'relations' u'burden'] | topic36\n",
      "topic37 | [u'spains' u'alexander' u'lindauer' u'wage' u'rodney'] | topic37\n",
      "topic38 | [u'signatures' u'jeff' u'pelosi' u'lois' u'filing'] | topic38\n",
      "topic4 | [u'engage' u'hispanics' u'oreilly' u'commerce' u'rapids'] | topic4\n",
      "topic40 | [u'delaware' u'utah' u'bond' u'mccain' u'thurlow'] | topic40\n",
      "topic41 | [u'murphy' u'musgrave' u'diedrich' u'matsunaka' u'larry'] | topic41\n",
      "topic42 | [u'tenenbaum' u'demint' u'martinez' u'inez' u'nebraska'] | topic42\n",
      "topic43 | [u'occupation' u'council' u'violence' u'prison' u'wounded'] | topic43\n",
      "topic44 | [u'income' u'cap' u'definition' u'unfit' u'novak'] | topic44\n",
      "topic180 | [u'cleveland' u'income' u'engage' u'petitions' u'cspan'] | topic44\n",
      "topic45 | [u'condis' u'declassified' u'connect' u'apt' u'beers'] | topic45\n",
      "topic46 | [u'wall' u'nuclear' u'donald' u'mcclellan' u'missile'] | topic46\n",
      "topic48 | [u'harris' u'cuba' u'schneider' u'enlarge' u'floridas'] | topic48\n",
      "topic5 | [u'invasion' u'postwar' u'uniform' u'taliban' u'safer'] | topic5\n",
      "topic50 | [u'laden' u'film' u'muslim' u'osama' u'juan'] | topic50\n",
      "topic51 | [u'limbaugh' u'sinclair' u'fantasy' u'flu' u'vaccine'] | topic51\n",
      "topic52 | [u'swift' u'star' u'veterans' u'boat' u'scott'] | topic52\n",
      "topic53 | [u'registrants' u'floridas' u'census' u'shells' u'beach'] | topic53\n",
      "topic54 | [u'sounds' u'highly' u'stand' u'anger' u'providing'] | topic54\n",
      "topic55 | [u'pdf' u'college' u'survey' u'percentage' u'arg'] | topic55\n",
      "topic56 | [u'cheneys' u'smith' u'walmart' u'hometown' u'mercenary'] | topic56\n",
      "topic57 | [u'fiscal' u'destruction' u'nuclear' u'donald' u'choices'] | topic57\n",
      "topic58 | [u'arab' u'armed' u'gen' u'fallujah' u'commander'] | topic58\n",
      "topic59 | [u'club' u'murkowski' u'lisa' u'indian' u'slot'] | topic59\n",
      "topic6 | [u'pdf' u'surveyusa' u'louisiana' u'tied' u'usa'] | topic6\n",
      "topic139 | [u'october' u'pdf' u'tied' u'sample' u'surveyusa'] | topic6\n",
      "topic60 | [u'soft' u'gov' u'latino' u'material' u'kean'] | topic60\n",
      "topic61 | [u'wage' u'cap' u'protesters' u'mcgreevey' u'award'] | topic61\n",
      "topic47 | [u'afscme' u'seiu' u'mcentee' u'coup' u'majette'] | topic62\n",
      "topic62 | [u'matsunaka' u'afscme' u'seiu' u'nag' u'boy'] | topic62\n",
      "topic63 | [u'oregon' u'bounce' u'finish' u'unlike' u'electorate'] | topic63\n",
      "topic64 | [u'industrial' u'settled' u'passion' u'nod' u'soros'] | topic64\n",
      "topic65 | [u'tariffs' u'consumer' u'trippi' u'europeans' u'steel'] | topic65\n",
      "topic66 | [u'trial' u'stations' u'claimed' u'declared' u'host'] | topic66\n",
      "topic67 | [u'cbs' u'weekend' u'combat' u'duty' u'served'] | topic67\n",
      "topic68 | [u'womens' u'saudi' u'arabia' u'planet' u'epa'] | topic68\n",
      "topic69 | [u'kerryedwards' u'legal' u'thoughts' u'catholic' u'rightwing'] | topic69\n",
      "topic7 | [u'operations' u'coalition' u'insurgents' u'command' u'combat'] | topic7\n",
      "topic71 | [u'alexander' u'bean' u'schrader' u'dccc' u'keever'] | topic71\n",
      "topic72 | [u'harris' u'gas' u'climate' u'guantanamo' u'legislation'] | topic72\n",
      "topic120 | [u'harris' u'smith' u'rogers' u'methods' u'block'] | topic72\n",
      "topic73 | [u'ethics' u'chandler' u'dakota' u'governors' u'ken'] | topic73\n",
      "topic74 | [u'blogger' u'indicted' u'earle' u'billmon' u'features'] | topic74\n",
      "topic75 | [u'required' u'filed' u'map' u'registrations' u'cuba'] | topic75\n",
      "topic76 | [u'rankings' u'cattle' u'pelosi' u'buzz' u'diedrichs'] | topic76\n",
      "topic77 | [u'deal' u'jun' u'reporters' u'efforts' u'chief'] | topic77\n",
      "topic78 | [u'ryan' u'allegations' u'hoeffel' u'lying' u'appointed'] | topic78\n",
      "topic79 | [u'clarke' u'richard' u'contributions' u'morrison' u'wilson'] | topic79\n",
      "topic80 | [u'upenn' u'drug' u'drugs' u'apologies' u'canadian'] | topic80\n",
      "topic83 | [u'reform' u'houston' u'debt' u'relief' u'nra'] | topic83\n",
      "topic84 | [u'dkos' u'research' u'difference' u'page' u'thousands'] | topic84\n",
      "topic85 | [u'chalabis' u'berg' u'mosul' u'rewards' u'contractor'] | topic85\n",
      "topic86 | [u'forecast' u'economist' u'economists' u'wall' u'rates'] | topic86\n",
      "topic88 | [u'libby' u'probe' u'plame' u'leak' u'thune'] | topic88\n",
      "topic89 | [u'wingnut' u'qpoll' u'requests' u'sterns' u'toomeys'] | topic89\n",
      "topic9 | [u'divorce' u'chicago' u'prices' u'boxer' u'flu'] | topic9\n",
      "topic90 | [u'bunning' u'caucuses' u'mongiardo' u'manager' u'volunteers'] | topic90\n",
      "topic92 | [u'black' u'werent' u'islam' u'factions' u'ideals'] | topic92\n",
      "topic93 | [u'postdebate' u'mad' u'outreach' u'steve' u'indictments'] | topic93\n",
      "topic94 | [u'congressman' u'districts' u'contest' u'congressional' u'helping'] | topic94\n",
      "topic95 | [u'geneva' u'resolution' u'beasley' u'questioning' u'impeachment'] | topic95\n",
      "topic96 | [u'lawyers' u'experts' u'justice' u'crisis' u'medical'] | topic96\n",
      "topic98 | [u'novak' u'falls' u'machines' u'frost' u'columnist'] | topic98\n",
      "topic99 | [u'nrsc' u'margins' u'senatorial' u'rosenberg' u'hart'] | topic99\n"
     ]
    }
   ],
   "source": [
    "exp.show_all_topics(sort_by_closest_topic=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "И оценить расстояния между темами."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.2902280390262604"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "exp.topics_pool.get_dist_between_topics('topic214', 'topic170')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Далее для обработки тем мы воспользуемся навигатором. Чтобы импортировать ранее созданные темы в него, нужно выбрать размер батча, который предстоит оценить асессорам. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "topic249:\n",
      "[u'theme' u'strongly' u'grassroots' u'counties' u'corps']\n",
      "topic209:\n",
      "[u'feb' u'research' u'rep' u'makes' u'means']\n",
      "topic64:\n",
      "[u'industrial' u'settled' u'passion' u'nod' u'soros']\n",
      "topic167:\n",
      "[u'rnc' u'smith' u'traffic' u'oregon' u'actual']\n",
      "topic182:\n",
      "[u'murkowski' u'keyes' u'barbara' u'tony' u'lisa']\n",
      "topic3:\n",
      "[u'employee' u'indictment' u'teachers' u'austin' u'mcgreevey']\n",
      "topic14:\n",
      "[u'citizen' u'matches' u'egon' u'casts' u'kerr']\n",
      "topic67:\n",
      "[u'cbs' u'weekend' u'combat' u'duty' u'served']\n",
      "topic7:\n",
      "[u'operations' u'coalition' u'insurgents' u'command' u'combat']\n",
      "topic199:\n",
      "[u'language' u'ideas' u'moores' u'substantial' u'providing']\n"
     ]
    }
   ],
   "source": [
    "exp.show_next_topics_batch(10)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Выбрав батч, далее мы сохраняем нужные данные в навигатор."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "exp.save_dataset_to_navigator('kos')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "exp.save_next_topics_batch_to_navigator(10, 'kos')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "И после оценок асессоров мы можем загрузить их обратно и свободно работать с ними."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "exp.load_assessments_from_navigator()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "exp.show_assessments()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}