ml3lab.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Bach Chorales Harmony\n",
    "\n",
    "Учитывая музыкальный поток, задача анализа гармонии музыки состоит в том, чтобы связать метку аккорда с каждой точкой времени. Такие ярлыки показывают основную гармонию, указывая фундаментальную ноту (корень) и режим, используя имена аккордов, такие как C minor.\n",
    "\n",
    "Набор данных состоит из 60 хоралов (5665 событий) J.S. Бах (1675-1750). Каждое событие каждого хорала помечено, используя 1 из 101 аккордовых меток.\n",
    "Количество атрибутов: 17 (имя последовательности, номер события, наличие метки (x12), аккорд).\n",
    "\n",
    "Информация о высотных классах была извлечена из загруженных источников MIDI от (JSB Chorales) [http://www.jsbchorales.net]. Информация о метре были рассчитаны через программу Meter, которая является частью Melisma музыкальный анализатор (Melisma) [http://www.link.cs.cmu.edu/music-analysis/]. Рукописные ярлыки вручную аннотируются экспертом-человеком.\n",
    "\n",
    "## Описание данных\n",
    "\n",
    "   * ***Choral ID***: соответствует именам файлов из (Bach Central) [http://www.bachcentral.com];\n",
    "   * ***Event number***:  индекс (начиная с 1) события внутри хорала;\n",
    "   * ***Pitch classes(C..B) (x12)***: YES / NO в зависимости от того, присутствует ли данный высотный класс [https://ru.wikipedia.org/wiki/Высотный_класс];\n",
    "   * ***Pitch class Bass***: высотный класс басовой ноты;\n",
    "   * ***Meter***: целые числа от 1 до 5, меньшие числа обозначают менее акцентированные события, более высокие числа обозначают более акцентированные события;\n",
    "   * ***Chord label***: аккорд, резонирующий во время данного события."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Подключаем используемые в работе библиотеки"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy  as np\n",
    "import matplotlib.pyplot  as plt\n",
    "import pandas as pd\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Считываем данные из датасета"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "data = pd.read_csv('source/jsbach_chorals_harmony.data', sep = ',', engine = 'python')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Контроль ввода"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Choral ID</th>\n",
       "      <th>Event number</th>\n",
       "      <th>Pitch class C</th>\n",
       "      <th>Pitch class C#</th>\n",
       "      <th>Pitch class D</th>\n",
       "      <th>Pitch class D#</th>\n",
       "      <th>Pitch class E</th>\n",
       "      <th>Pitch class F</th>\n",
       "      <th>Pitch class F#</th>\n",
       "      <th>Pitch class G</th>\n",
       "      <th>Pitch class G#</th>\n",
       "      <th>Pitch class A</th>\n",
       "      <th>Pitch class A#</th>\n",
       "      <th>Pitch class B</th>\n",
       "      <th>Pitch class Bass</th>\n",
       "      <th>Meter</th>\n",
       "      <th>Chord label</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>1</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>2</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>E</td>\n",
       "      <td>5</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>3</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>E</td>\n",
       "      <td>2</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>4</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>5</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>2</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>6</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>4</td>\n",
       "      <td>D_m</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>7</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>2</td>\n",
       "      <td>D_m</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>8</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>A</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>9</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>A</td>\n",
       "      <td>2</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>10</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>Bb</td>\n",
       "      <td>5</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>11</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>Bb</td>\n",
       "      <td>1</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>12</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>2</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>13</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>1</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>14</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>3</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>15</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>E</td>\n",
       "      <td>2</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>16</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>4</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>17</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>18</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>Bb</td>\n",
       "      <td>5</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>19</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>Bb</td>\n",
       "      <td>2</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>20</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>Bb</td>\n",
       "      <td>1</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>21</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>A</td>\n",
       "      <td>3</td>\n",
       "      <td>C_M7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>22</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>G</td>\n",
       "      <td>2</td>\n",
       "      <td>C_M7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>23</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>G</td>\n",
       "      <td>1</td>\n",
       "      <td>C_M7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>24</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>A</td>\n",
       "      <td>4</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>25</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>B</td>\n",
       "      <td>2</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>26</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>3</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>27</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>2</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>28</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>5</td>\n",
       "      <td>D_m7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>29</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>F</td>\n",
       "      <td>2</td>\n",
       "      <td>D_m7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>000106b_</td>\n",
       "      <td>30</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>1</td>\n",
       "      <td>D_m7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5635</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>80</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>5</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5636</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>81</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5637</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>82</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>4</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5638</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>83</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>G</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5639</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>84</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>A</td>\n",
       "      <td>5</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5640</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>85</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>A</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5641</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>86</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>Bb</td>\n",
       "      <td>4</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5642</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>87</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>3</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5643</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>88</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>5</td>\n",
       "      <td>D_m</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5644</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>89</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>3</td>\n",
       "      <td>D_m</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5645</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>90</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>E</td>\n",
       "      <td>4</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5646</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>91</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F#</td>\n",
       "      <td>3</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5647</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>92</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>G</td>\n",
       "      <td>5</td>\n",
       "      <td>G_m</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5648</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>93</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>G</td>\n",
       "      <td>3</td>\n",
       "      <td>G_m</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5649</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>94</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>4</td>\n",
       "      <td>D_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5650</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>95</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>A</td>\n",
       "      <td>5</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5651</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>96</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>Bb</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5652</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>97</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>4</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5653</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>98</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>3</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5654</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>99</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>5</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5655</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>100</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>E</td>\n",
       "      <td>3</td>\n",
       "      <td>BbM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5656</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>101</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>4</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5657</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>102</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>3</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5658</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>103</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>E</td>\n",
       "      <td>5</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5659</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>104</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>3</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5660</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>105</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>G</td>\n",
       "      <td>4</td>\n",
       "      <td>G_m</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5661</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>106</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>G</td>\n",
       "      <td>3</td>\n",
       "      <td>G_m</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5662</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>107</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>5</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5663</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>108</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>C</td>\n",
       "      <td>3</td>\n",
       "      <td>C_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5664</th>\n",
       "      <td>015505b_</td>\n",
       "      <td>109</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>4</td>\n",
       "      <td>F_M</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5665 rows × 17 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     Choral ID  Event number Pitch class C Pitch class C# Pitch class D  \\\n",
       "0     000106b_             1           YES             NO            NO   \n",
       "1     000106b_             2           YES             NO            NO   \n",
       "2     000106b_             3           YES             NO            NO   \n",
       "3     000106b_             4           YES             NO            NO   \n",
       "4     000106b_             5           YES             NO            NO   \n",
       "5     000106b_             6            NO             NO           YES   \n",
       "6     000106b_             7            NO             NO           YES   \n",
       "7     000106b_             8           YES             NO            NO   \n",
       "8     000106b_             9           YES             NO            NO   \n",
       "9     000106b_            10            NO             NO           YES   \n",
       "10    000106b_            11            NO             NO           YES   \n",
       "11    000106b_            12           YES             NO           YES   \n",
       "12    000106b_            13           YES             NO           YES   \n",
       "13    000106b_            14            NO             NO           YES   \n",
       "14    000106b_            15            NO             NO           YES   \n",
       "15    000106b_            16           YES             NO            NO   \n",
       "16    000106b_            17           YES             NO            NO   \n",
       "17    000106b_            18            NO             NO           YES   \n",
       "18    000106b_            19            NO             NO           YES   \n",
       "19    000106b_            20            NO             NO           YES   \n",
       "20    000106b_            21            NO             NO            NO   \n",
       "21    000106b_            22           YES             NO            NO   \n",
       "22    000106b_            23            NO             NO            NO   \n",
       "23    000106b_            24           YES             NO            NO   \n",
       "24    000106b_            25            NO             NO           YES   \n",
       "25    000106b_            26           YES             NO            NO   \n",
       "26    000106b_            27           YES             NO            NO   \n",
       "27    000106b_            28           YES             NO           YES   \n",
       "28    000106b_            29            NO             NO           YES   \n",
       "29    000106b_            30            NO             NO           YES   \n",
       "...        ...           ...           ...            ...           ...   \n",
       "5635  015505b_            80           YES             NO            NO   \n",
       "5636  015505b_            81           YES             NO            NO   \n",
       "5637  015505b_            82           YES             NO            NO   \n",
       "5638  015505b_            83           YES             NO            NO   \n",
       "5639  015505b_            84           YES             NO            NO   \n",
       "5640  015505b_            85           YES             NO            NO   \n",
       "5641  015505b_            86            NO             NO           YES   \n",
       "5642  015505b_            87           YES             NO           YES   \n",
       "5643  015505b_            88            NO             NO           YES   \n",
       "5644  015505b_            89            NO             NO           YES   \n",
       "5645  015505b_            90           YES             NO            NO   \n",
       "5646  015505b_            91           YES             NO            NO   \n",
       "5647  015505b_            92            NO             NO           YES   \n",
       "5648  015505b_            93           YES             NO            NO   \n",
       "5649  015505b_            94            NO             NO           YES   \n",
       "5650  015505b_            95           YES             NO            NO   \n",
       "5651  015505b_            96           YES             NO            NO   \n",
       "5652  015505b_            97           YES             NO            NO   \n",
       "5653  015505b_            98           YES             NO            NO   \n",
       "5654  015505b_            99            NO             NO           YES   \n",
       "5655  015505b_           100            NO             NO           YES   \n",
       "5656  015505b_           101           YES             NO            NO   \n",
       "5657  015505b_           102            NO             NO           YES   \n",
       "5658  015505b_           103           YES             NO            NO   \n",
       "5659  015505b_           104           YES             NO            NO   \n",
       "5660  015505b_           105            NO             NO           YES   \n",
       "5661  015505b_           106            NO             NO           YES   \n",
       "5662  015505b_           107           YES             NO            NO   \n",
       "5663  015505b_           108           YES             NO            NO   \n",
       "5664  015505b_           109            NO             NO            NO   \n",
       "\n",
       "     Pitch class D# Pitch class E Pitch class F Pitch class F# Pitch class G  \\\n",
       "0                NO            NO           YES             NO            NO   \n",
       "1                NO           YES            NO             NO           YES   \n",
       "2                NO           YES            NO             NO           YES   \n",
       "3                NO            NO           YES             NO            NO   \n",
       "4                NO            NO           YES             NO            NO   \n",
       "5                NO            NO           YES             NO            NO   \n",
       "6                NO            NO           YES             NO            NO   \n",
       "7                NO            NO           YES             NO            NO   \n",
       "8                NO            NO           YES             NO            NO   \n",
       "9                NO            NO           YES             NO            NO   \n",
       "10               NO            NO           YES             NO            NO   \n",
       "11               NO            NO           YES             NO            NO   \n",
       "12               NO            NO           YES             NO           YES   \n",
       "13               NO            NO           YES             NO            NO   \n",
       "14               NO           YES           YES             NO           YES   \n",
       "15               NO            NO           YES             NO            NO   \n",
       "16               NO            NO           YES             NO            NO   \n",
       "17               NO            NO           YES             NO            NO   \n",
       "18               NO            NO           YES             NO           YES   \n",
       "19               NO            NO           YES             NO            NO   \n",
       "20               NO           YES            NO             NO           YES   \n",
       "21               NO           YES            NO             NO           YES   \n",
       "22               NO           YES            NO             NO           YES   \n",
       "23               NO            NO           YES             NO            NO   \n",
       "24               NO            NO           YES             NO           YES   \n",
       "25               NO           YES            NO             NO           YES   \n",
       "26               NO           YES            NO             NO           YES   \n",
       "27               NO            NO           YES             NO            NO   \n",
       "28               NO           YES           YES             NO            NO   \n",
       "29               NO           YES           YES             NO            NO   \n",
       "...             ...           ...           ...            ...           ...   \n",
       "5635             NO            NO           YES             NO            NO   \n",
       "5636             NO            NO           YES             NO           YES   \n",
       "5637             NO            NO           YES             NO            NO   \n",
       "5638             NO            NO           YES             NO           YES   \n",
       "5639             NO            NO           YES             NO            NO   \n",
       "5640             NO           YES            NO             NO           YES   \n",
       "5641             NO            NO           YES             NO            NO   \n",
       "5642             NO            NO           YES             NO            NO   \n",
       "5643             NO            NO           YES             NO            NO   \n",
       "5644             NO            NO           YES             NO            NO   \n",
       "5645             NO           YES            NO             NO           YES   \n",
       "5646             NO           YES            NO            YES            NO   \n",
       "5647             NO            NO            NO             NO           YES   \n",
       "5648             NO           YES            NO             NO           YES   \n",
       "5649             NO            NO            NO            YES            NO   \n",
       "5650             NO            NO           YES             NO            NO   \n",
       "5651             NO            NO           YES             NO            NO   \n",
       "5652             NO           YES            NO             NO           YES   \n",
       "5653             NO           YES            NO             NO            NO   \n",
       "5654             NO            NO           YES             NO            NO   \n",
       "5655             NO           YES            NO             NO           YES   \n",
       "5656             NO            NO           YES             NO            NO   \n",
       "5657             NO            NO           YES             NO            NO   \n",
       "5658             NO           YES            NO             NO           YES   \n",
       "5659             NO           YES           YES             NO           YES   \n",
       "5660             NO            NO            NO             NO           YES   \n",
       "5661             NO            NO            NO             NO           YES   \n",
       "5662             NO           YES            NO             NO           YES   \n",
       "5663             NO           YES            NO             NO           YES   \n",
       "5664             NO            NO           YES             NO            NO   \n",
       "\n",
       "     Pitch class G# Pitch class A Pitch class A# Pitch class B  \\\n",
       "0                NO           YES             NO            NO   \n",
       "1                NO            NO             NO            NO   \n",
       "2                NO            NO             NO            NO   \n",
       "3                NO           YES             NO            NO   \n",
       "4                NO           YES             NO            NO   \n",
       "5                NO           YES             NO            NO   \n",
       "6                NO           YES             NO            NO   \n",
       "7                NO           YES             NO            NO   \n",
       "8                NO           YES             NO            NO   \n",
       "9                NO            NO            YES            NO   \n",
       "10               NO            NO            YES            NO   \n",
       "11               NO           YES            YES            NO   \n",
       "12               NO            NO            YES            NO   \n",
       "13               NO            NO            YES            NO   \n",
       "14               NO            NO            YES            NO   \n",
       "15               NO           YES             NO            NO   \n",
       "16               NO           YES             NO            NO   \n",
       "17               NO            NO            YES            NO   \n",
       "18               NO            NO            YES            NO   \n",
       "19               NO           YES            YES            NO   \n",
       "20               NO           YES            YES            NO   \n",
       "21               NO            NO            YES            NO   \n",
       "22               NO            NO            YES            NO   \n",
       "23               NO           YES             NO            NO   \n",
       "24               NO            NO             NO           YES   \n",
       "25               NO            NO             NO            NO   \n",
       "26               NO            NO             NO            NO   \n",
       "27               NO           YES             NO            NO   \n",
       "28               NO           YES             NO           YES   \n",
       "29               NO           YES             NO            NO   \n",
       "...             ...           ...            ...           ...   \n",
       "5635             NO           YES             NO            NO   \n",
       "5636             NO            NO            YES            NO   \n",
       "5637             NO           YES             NO            NO   \n",
       "5638             NO            NO            YES            NO   \n",
       "5639             NO           YES             NO            NO   \n",
       "5640             NO           YES             NO            NO   \n",
       "5641             NO            NO            YES            NO   \n",
       "5642             NO            NO             NO            NO   \n",
       "5643             NO           YES             NO            NO   \n",
       "5644             NO            NO            YES            NO   \n",
       "5645             NO            NO             NO            NO   \n",
       "5646             NO           YES             NO            NO   \n",
       "5647             NO            NO            YES            NO   \n",
       "5648             NO            NO            YES            NO   \n",
       "5649             NO           YES             NO            NO   \n",
       "5650             NO           YES             NO            NO   \n",
       "5651             NO           YES            YES            NO   \n",
       "5652             NO            NO             NO            NO   \n",
       "5653             NO           YES             NO            NO   \n",
       "5654             NO            NO            YES            NO   \n",
       "5655             NO            NO            YES            NO   \n",
       "5656             NO           YES             NO            NO   \n",
       "5657             NO           YES            YES            NO   \n",
       "5658             NO            NO             NO            NO   \n",
       "5659             NO            NO             NO            NO   \n",
       "5660             NO            NO            YES            NO   \n",
       "5661             NO           YES             NO            NO   \n",
       "5662             NO            NO             NO            NO   \n",
       "5663             NO            NO            YES            NO   \n",
       "5664             NO           YES             NO            NO   \n",
       "\n",
       "     Pitch class Bass  Meter Chord label  \n",
       "0                   F      3         F_M  \n",
       "1                   E      5         C_M  \n",
       "2                   E      2         C_M  \n",
       "3                   F      3         F_M  \n",
       "4                   F      2         F_M  \n",
       "5                   D      4         D_m  \n",
       "6                   D      2         D_m  \n",
       "7                   A      3         F_M  \n",
       "8                   A      2         F_M  \n",
       "9                  Bb      5         BbM  \n",
       "10                 Bb      1         BbM  \n",
       "11                  C      2         BbM  \n",
       "12                  C      1         BbM  \n",
       "13                  D      3         BbM  \n",
       "14                  E      2         BbM  \n",
       "15                  F      4         F_M  \n",
       "16                  F      3         F_M  \n",
       "17                 Bb      5         BbM  \n",
       "18                 Bb      2         BbM  \n",
       "19                 Bb      1         BbM  \n",
       "20                  A      3        C_M7  \n",
       "21                  G      2        C_M7  \n",
       "22                  G      1        C_M7  \n",
       "23                  A      4         F_M  \n",
       "24                  B      2         F_M  \n",
       "25                  C      3         C_M  \n",
       "26                  C      2         C_M  \n",
       "27                  F      5        D_m7  \n",
       "28                  F      2        D_m7  \n",
       "29                  F      1        D_m7  \n",
       "...               ...    ...         ...  \n",
       "5635                F      5         F_M  \n",
       "5636                F      3         F_M  \n",
       "5637                F      4         F_M  \n",
       "5638                G      3         F_M  \n",
       "5639                A      5         F_M  \n",
       "5640                A      3         F_M  \n",
       "5641               Bb      4         BbM  \n",
       "5642                C      3         BbM  \n",
       "5643                D      5         D_m  \n",
       "5644                D      3         D_m  \n",
       "5645                E      4         C_M  \n",
       "5646               F#      3         C_M  \n",
       "5647                G      5         G_m  \n",
       "5648                G      3         G_m  \n",
       "5649                D      4         D_M  \n",
       "5650                A      5         F_M  \n",
       "5651               Bb      3         F_M  \n",
       "5652                C      4         C_M  \n",
       "5653                C      3         C_M  \n",
       "5654                D      5         BbM  \n",
       "5655                E      3         BbM  \n",
       "5656                F      4         F_M  \n",
       "5657                F      3         F_M  \n",
       "5658                E      5         C_M  \n",
       "5659                F      3         C_M  \n",
       "5660                G      4         G_m  \n",
       "5661                G      3         G_m  \n",
       "5662                C      5         C_M  \n",
       "5663                C      3         C_M  \n",
       "5664                F      4         F_M  \n",
       "\n",
       "[5665 rows x 17 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Контроль размеров таблицы и наименования признаков(столбцов)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(5665, 17)\n",
      "Index(['Choral ID', 'Event number', 'Pitch class C', 'Pitch class C#',\n",
      "       'Pitch class D', 'Pitch class D#', 'Pitch class E', 'Pitch class F',\n",
      "       'Pitch class F#', 'Pitch class G', 'Pitch class G#', 'Pitch class A',\n",
      "       'Pitch class A#', 'Pitch class B', 'Pitch class Bass', 'Meter',\n",
      "       'Chord label'],\n",
      "      dtype='object')\n"
     ]
    }
   ],
   "source": [
    "print(data.shape)\n",
    "print(data.columns)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Основные сведения по данным, из таблицы видно, что пропусков в данных нет, уже хорошо, большинство признаков бинарные."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Choral ID</th>\n",
       "      <th>Pitch class C</th>\n",
       "      <th>Pitch class C#</th>\n",
       "      <th>Pitch class D</th>\n",
       "      <th>Pitch class D#</th>\n",
       "      <th>Pitch class E</th>\n",
       "      <th>Pitch class F</th>\n",
       "      <th>Pitch class F#</th>\n",
       "      <th>Pitch class G</th>\n",
       "      <th>Pitch class G#</th>\n",
       "      <th>Pitch class A</th>\n",
       "      <th>Pitch class A#</th>\n",
       "      <th>Pitch class B</th>\n",
       "      <th>Pitch class Bass</th>\n",
       "      <th>Chord label</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>unique</th>\n",
       "      <td>62</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>16</td>\n",
       "      <td>102</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>top</th>\n",
       "      <td>002908ch</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>D_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>freq</th>\n",
       "      <td>206</td>\n",
       "      <td>3875</td>\n",
       "      <td>4711</td>\n",
       "      <td>3300</td>\n",
       "      <td>4956</td>\n",
       "      <td>3540</td>\n",
       "      <td>4381</td>\n",
       "      <td>4253</td>\n",
       "      <td>3523</td>\n",
       "      <td>5006</td>\n",
       "      <td>3290</td>\n",
       "      <td>4644</td>\n",
       "      <td>3874</td>\n",
       "      <td>689</td>\n",
       "      <td>503</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Choral ID Pitch class C Pitch class C# Pitch class D Pitch class D#  \\\n",
       "count       5665          5665           5665          5665           5665   \n",
       "unique        62             2              2             2              2   \n",
       "top     002908ch            NO             NO            NO             NO   \n",
       "freq         206          3875           4711          3300           4956   \n",
       "\n",
       "       Pitch class E Pitch class F Pitch class F# Pitch class G  \\\n",
       "count           5665          5665           5665          5665   \n",
       "unique             2             2              2             2   \n",
       "top               NO            NO             NO            NO   \n",
       "freq            3540          4381           4253          3523   \n",
       "\n",
       "       Pitch class G# Pitch class A Pitch class A# Pitch class B  \\\n",
       "count            5665          5665           5665          5665   \n",
       "unique              2             2              2             2   \n",
       "top                NO            NO             NO            NO   \n",
       "freq             5006          3290           4644          3874   \n",
       "\n",
       "       Pitch class Bass Chord label  \n",
       "count              5665        5665  \n",
       "unique               16         102  \n",
       "top                   D         D_M  \n",
       "freq                689         503  "
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.describe(include = ['object'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Слияние двух таблиц выше, максимальная полнота сведений о данных"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Choral ID</th>\n",
       "      <th>Event number</th>\n",
       "      <th>Pitch class C</th>\n",
       "      <th>Pitch class C#</th>\n",
       "      <th>Pitch class D</th>\n",
       "      <th>Pitch class D#</th>\n",
       "      <th>Pitch class E</th>\n",
       "      <th>Pitch class F</th>\n",
       "      <th>Pitch class F#</th>\n",
       "      <th>Pitch class G</th>\n",
       "      <th>Pitch class G#</th>\n",
       "      <th>Pitch class A</th>\n",
       "      <th>Pitch class A#</th>\n",
       "      <th>Pitch class B</th>\n",
       "      <th>Pitch class Bass</th>\n",
       "      <th>Meter</th>\n",
       "      <th>Chord label</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>5665</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>unique</th>\n",
       "      <td>62</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>16</td>\n",
       "      <td>NaN</td>\n",
       "      <td>102</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>top</th>\n",
       "      <td>002908ch</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>D</td>\n",
       "      <td>NaN</td>\n",
       "      <td>D_M</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>freq</th>\n",
       "      <td>206</td>\n",
       "      <td>NaN</td>\n",
       "      <td>3875</td>\n",
       "      <td>4711</td>\n",
       "      <td>3300</td>\n",
       "      <td>4956</td>\n",
       "      <td>3540</td>\n",
       "      <td>4381</td>\n",
       "      <td>4253</td>\n",
       "      <td>3523</td>\n",
       "      <td>5006</td>\n",
       "      <td>3290</td>\n",
       "      <td>4644</td>\n",
       "      <td>3874</td>\n",
       "      <td>689</td>\n",
       "      <td>NaN</td>\n",
       "      <td>503</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>NaN</td>\n",
       "      <td>53.374404</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>3.134863</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>NaN</td>\n",
       "      <td>37.268208</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.109710</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>NaN</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>NaN</td>\n",
       "      <td>24.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>NaN</td>\n",
       "      <td>48.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>NaN</td>\n",
       "      <td>75.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>4.000000</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>NaN</td>\n",
       "      <td>207.000000</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>5.000000</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Choral ID  Event number Pitch class C Pitch class C# Pitch class D  \\\n",
       "count       5665   5665.000000          5665           5665          5665   \n",
       "unique        62           NaN             2              2             2   \n",
       "top     002908ch           NaN            NO             NO            NO   \n",
       "freq         206           NaN          3875           4711          3300   \n",
       "mean         NaN     53.374404           NaN            NaN           NaN   \n",
       "std          NaN     37.268208           NaN            NaN           NaN   \n",
       "min          NaN      1.000000           NaN            NaN           NaN   \n",
       "25%          NaN     24.000000           NaN            NaN           NaN   \n",
       "50%          NaN     48.000000           NaN            NaN           NaN   \n",
       "75%          NaN     75.000000           NaN            NaN           NaN   \n",
       "max          NaN    207.000000           NaN            NaN           NaN   \n",
       "\n",
       "       Pitch class D# Pitch class E Pitch class F Pitch class F#  \\\n",
       "count            5665          5665          5665           5665   \n",
       "unique              2             2             2              2   \n",
       "top                NO            NO            NO             NO   \n",
       "freq             4956          3540          4381           4253   \n",
       "mean              NaN           NaN           NaN            NaN   \n",
       "std               NaN           NaN           NaN            NaN   \n",
       "min               NaN           NaN           NaN            NaN   \n",
       "25%               NaN           NaN           NaN            NaN   \n",
       "50%               NaN           NaN           NaN            NaN   \n",
       "75%               NaN           NaN           NaN            NaN   \n",
       "max               NaN           NaN           NaN            NaN   \n",
       "\n",
       "       Pitch class G Pitch class G# Pitch class A Pitch class A#  \\\n",
       "count           5665           5665          5665           5665   \n",
       "unique             2              2             2              2   \n",
       "top               NO             NO            NO             NO   \n",
       "freq            3523           5006          3290           4644   \n",
       "mean             NaN            NaN           NaN            NaN   \n",
       "std              NaN            NaN           NaN            NaN   \n",
       "min              NaN            NaN           NaN            NaN   \n",
       "25%              NaN            NaN           NaN            NaN   \n",
       "50%              NaN            NaN           NaN            NaN   \n",
       "75%              NaN            NaN           NaN            NaN   \n",
       "max              NaN            NaN           NaN            NaN   \n",
       "\n",
       "       Pitch class B Pitch class Bass        Meter Chord label  \n",
       "count           5665             5665  5665.000000        5665  \n",
       "unique             2               16          NaN         102  \n",
       "top               NO                D          NaN         D_M  \n",
       "freq            3874              689          NaN         503  \n",
       "mean             NaN              NaN     3.134863         NaN  \n",
       "std              NaN              NaN     1.109710         NaN  \n",
       "min              NaN              NaN     1.000000         NaN  \n",
       "25%              NaN              NaN     2.000000         NaN  \n",
       "50%              NaN              NaN     3.000000         NaN  \n",
       "75%              NaN              NaN     4.000000         NaN  \n",
       "max              NaN              NaN     5.000000         NaN  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.describe(include = 'all')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Подготовка данных\n",
    "\n",
    "   Выбросим признак Choral ID, поскольку от названия произведения мало что зависит, также мало что зависит от порядкового номера события, значит Event Number выбрасываем туда-же. \n",
    "   Задача состоит в предсказании резонирующего в данный момент аккорда по встречающимся или не встречающимся высотным классам основного исполнения и баса. Большинство признаков бинарных, остаётся только придумать, что делать с категориальным признаком Pitch Class Bass, например векторизовать его."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "data = data.drop('Choral ID', axis=1)\n",
    "data = data.drop('Event number', axis=1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Поскольку столбец Chord label содержит ответ, уберём из таблицы и его"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "y = data['Chord label']\n",
    "data = data.drop('Chord label', axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Pitch class C</th>\n",
       "      <th>Pitch class C#</th>\n",
       "      <th>Pitch class D</th>\n",
       "      <th>Pitch class D#</th>\n",
       "      <th>Pitch class E</th>\n",
       "      <th>Pitch class F</th>\n",
       "      <th>Pitch class F#</th>\n",
       "      <th>Pitch class G</th>\n",
       "      <th>Pitch class G#</th>\n",
       "      <th>Pitch class A</th>\n",
       "      <th>Pitch class A#</th>\n",
       "      <th>Pitch class B</th>\n",
       "      <th>Pitch class Bass</th>\n",
       "      <th>Meter</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>E</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>E</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>YES</td>\n",
       "      <td>NO</td>\n",
       "      <td>NO</td>\n",
       "      <td>F</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  Pitch class C Pitch class C# Pitch class D Pitch class D# Pitch class E  \\\n",
       "0           YES             NO            NO             NO            NO   \n",
       "1           YES             NO            NO             NO           YES   \n",
       "2           YES             NO            NO             NO           YES   \n",
       "3           YES             NO            NO             NO            NO   \n",
       "4           YES             NO            NO             NO            NO   \n",
       "\n",
       "  Pitch class F Pitch class F# Pitch class G Pitch class G# Pitch class A  \\\n",
       "0           YES             NO            NO             NO           YES   \n",
       "1            NO             NO           YES             NO            NO   \n",
       "2            NO             NO           YES             NO            NO   \n",
       "3           YES             NO            NO             NO           YES   \n",
       "4           YES             NO            NO             NO           YES   \n",
       "\n",
       "  Pitch class A# Pitch class B Pitch class Bass  Meter  \n",
       "0             NO            NO                F      3  \n",
       "1             NO            NO                E      5  \n",
       "2             NO            NO                E      2  \n",
       "3             NO            NO                F      3  \n",
       "4             NO            NO                F      2  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['F', 'E', 'D', 'A', 'Bb', 'C', 'G', 'B', 'F#', 'C#', 'Eb', 'D#',\n",
       "       'G#', 'A#', 'Ab', 'Db'], dtype=object)"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data['Pitch class Bass'].unique()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Имеются повторяющиеся с точки зрения музыки, но разные по названию ноты (классы), приведём их к стандарту как в заголовке таблицы для основного звучания(C, C# = Db, D, D# = Eb, E, F, F#, G, G# = Ab, A, A# = Bb, B)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "d = {'Bb' : 'A#', 'Db' : 'C#', 'Eb': 'D#', 'Ab' : 'G#'}\n",
    "data = data.replace({'Pitch class Bass': d})\n",
    "#data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Разобьём таблицу на части по типу признака"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['Pitch class C', 'Pitch class C#', 'Pitch class D', 'Pitch class D#', 'Pitch class E', 'Pitch class F', 'Pitch class F#', 'Pitch class G', 'Pitch class G#', 'Pitch class A', 'Pitch class A#', 'Pitch class B', 'Pitch class Bass']\n",
      "['Meter']\n"
     ]
    }
   ],
   "source": [
    "categorical_columns = [c for c in data.columns if data[c].dtype.name == 'object']\n",
    "numerical_columns   = [c for c in data.columns if data[c].dtype.name != 'object']\n",
    "print(categorical_columns)\n",
    "print(numerical_columns)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['Pitch class C', 'Pitch class C#', 'Pitch class D', 'Pitch class D#', 'Pitch class E', 'Pitch class F', 'Pitch class F#', 'Pitch class G', 'Pitch class G#', 'Pitch class A', 'Pitch class A#', 'Pitch class B'] ['Pitch class Bass']\n"
     ]
    }
   ],
   "source": [
    "data_describe = data.describe(include = [object])\n",
    "binary_columns    = [c for c in categorical_columns if data_describe[c]['unique'] == 2]\n",
    "nonbinary_columns = [c for c in categorical_columns if data_describe[c]['unique'] > 2]\n",
    "print(binary_columns, nonbinary_columns)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Бинарный признак делаем с циферками"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "data.at[data['Pitch class C'] != 'YES', 'Pitch class C'] = 0\n",
    "data.at[data['Pitch class C'] == 'YES', 'Pitch class C'] = 1\n",
    "data.at[data['Pitch class C#'] != 'YES', 'Pitch class C#'] = 0\n",
    "data.at[data['Pitch class C#'] == 'YES', 'Pitch class C#'] = 1\n",
    "data.at[data['Pitch class D'] != 'YES', 'Pitch class D'] = 0\n",
    "data.at[data['Pitch class D'] == 'YES', 'Pitch class D'] = 1\n",
    "data.at[data['Pitch class D#'] != 'YES', 'Pitch class D#'] = 0\n",
    "data.at[data['Pitch class D#'] == 'YES', 'Pitch class D#'] = 1\n",
    "data.at[data['Pitch class E'] != 'YES', 'Pitch class E'] = 0\n",
    "data.at[data['Pitch class E'] == 'YES', 'Pitch class E'] = 1\n",
    "data.at[data['Pitch class F'] != 'YES', 'Pitch class F'] = 0\n",
    "data.at[data['Pitch class F'] == 'YES', 'Pitch class F'] = 1\n",
    "data.at[data['Pitch class F#'] != 'YES', 'Pitch class F#'] = 0\n",
    "data.at[data['Pitch class F#'] == 'YES', 'Pitch class F#'] = 1\n",
    "data.at[data['Pitch class G'] != 'YES', 'Pitch class G'] = 0\n",
    "data.at[data['Pitch class G'] == 'YES', 'Pitch class G'] = 1\n",
    "data.at[data['Pitch class G#'] != 'YES', 'Pitch class G#'] = 0\n",
    "data.at[data['Pitch class G#'] == 'YES', 'Pitch class G#'] = 1\n",
    "data.at[data['Pitch class A'] != 'YES', 'Pitch class A'] = 0\n",
    "data.at[data['Pitch class A'] == 'YES', 'Pitch class A'] = 1\n",
    "data.at[data['Pitch class A#'] != 'YES', 'Pitch class A#'] = 0\n",
    "data.at[data['Pitch class A#'] == 'YES', 'Pitch class A#'] = 1\n",
    "data.at[data['Pitch class B'] != 'YES', 'Pitch class B'] = 0\n",
    "data.at[data['Pitch class B'] == 'YES', 'Pitch class B'] = 1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Векторизуем единственный оставшийся категориальный признак с помощью метода get_dummies из pandas:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Index(['Pitch class Bass_A', 'Pitch class Bass_A#', 'Pitch class Bass_B',\n",
      "       'Pitch class Bass_C', 'Pitch class Bass_C#', 'Pitch class Bass_D',\n",
      "       'Pitch class Bass_D#', 'Pitch class Bass_E', 'Pitch class Bass_F',\n",
      "       'Pitch class Bass_F#', 'Pitch class Bass_G', 'Pitch class Bass_G#'],\n",
      "      dtype='object')\n"
     ]
    }
   ],
   "source": [
    "data_nonbinary = pd.get_dummies(data[nonbinary_columns])\n",
    "print(data_nonbinary.columns)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Нормализация количественных признаков\n",
    "\n",
    "Многие алгоритмы машинного обучения чувствительны к масштабированию данных. К таким алгоритмам, например, относится метод ближайших соседей, машина опорных векторов и др.\n",
    "\n",
    "В этом случае количественные признаки полезно нормализовать. Это можно делать разными способами. Например, каждый количественный признак приведем к нулевому среднему и единичному среднеквадратичному отклонению:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Meter</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>5665.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>3.134863</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>1.109710</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>2.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>3.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>4.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>5.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             Meter\n",
       "count  5665.000000\n",
       "mean      3.134863\n",
       "std       1.109710\n",
       "min       1.000000\n",
       "25%       2.000000\n",
       "50%       3.000000\n",
       "75%       4.000000\n",
       "max       5.000000"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data_numerical = data[numerical_columns]\n",
    "data_numerical.describe()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_numerical = (data_numerical - data_numerical.mean(axis = 0))/data_numerical.std(axis = 0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Meter</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>5.665000e+03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>-1.053605e-15</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>1.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>-1.923803e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>-1.022667e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>-1.215302e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>7.796064e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>1.680743e+00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              Meter\n",
       "count  5.665000e+03\n",
       "mean  -1.053605e-15\n",
       "std    1.000000e+00\n",
       "min   -1.923803e+00\n",
       "25%   -1.022667e+00\n",
       "50%   -1.215302e-01\n",
       "75%    7.796064e-01\n",
       "max    1.680743e+00"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data_numerical.describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Соединяем всё в одну таблицу"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(5665, 25)\n"
     ]
    }
   ],
   "source": [
    "data = pd.concat((data_numerical, data_nonbinary, data[binary_columns]), axis = 1)\n",
    "print(data.shape)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Meter</th>\n",
       "      <th>Pitch class Bass_A</th>\n",
       "      <th>Pitch class Bass_A#</th>\n",
       "      <th>Pitch class Bass_B</th>\n",
       "      <th>Pitch class Bass_C</th>\n",
       "      <th>Pitch class Bass_C#</th>\n",
       "      <th>Pitch class Bass_D</th>\n",
       "      <th>Pitch class Bass_D#</th>\n",
       "      <th>Pitch class Bass_E</th>\n",
       "      <th>Pitch class Bass_F</th>\n",
       "      <th>...</th>\n",
       "      <th>Pitch class D</th>\n",
       "      <th>Pitch class D#</th>\n",
       "      <th>Pitch class E</th>\n",
       "      <th>Pitch class F</th>\n",
       "      <th>Pitch class F#</th>\n",
       "      <th>Pitch class G</th>\n",
       "      <th>Pitch class G#</th>\n",
       "      <th>Pitch class A</th>\n",
       "      <th>Pitch class A#</th>\n",
       "      <th>Pitch class B</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>5.665000e+03</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "      <td>5665.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>-1.053605e-15</td>\n",
       "      <td>0.121447</td>\n",
       "      <td>0.055075</td>\n",
       "      <td>0.079259</td>\n",
       "      <td>0.105914</td>\n",
       "      <td>0.038482</td>\n",
       "      <td>0.121624</td>\n",
       "      <td>0.044307</td>\n",
       "      <td>0.119682</td>\n",
       "      <td>0.081553</td>\n",
       "      <td>...</td>\n",
       "      <td>0.417476</td>\n",
       "      <td>0.125154</td>\n",
       "      <td>0.375110</td>\n",
       "      <td>0.226655</td>\n",
       "      <td>0.249250</td>\n",
       "      <td>0.378111</td>\n",
       "      <td>0.116328</td>\n",
       "      <td>0.419241</td>\n",
       "      <td>0.180229</td>\n",
       "      <td>0.316152</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>1.000000e+00</td>\n",
       "      <td>0.326675</td>\n",
       "      <td>0.228147</td>\n",
       "      <td>0.270166</td>\n",
       "      <td>0.307754</td>\n",
       "      <td>0.192374</td>\n",
       "      <td>0.326880</td>\n",
       "      <td>0.205795</td>\n",
       "      <td>0.324618</td>\n",
       "      <td>0.273707</td>\n",
       "      <td>...</td>\n",
       "      <td>0.493186</td>\n",
       "      <td>0.330923</td>\n",
       "      <td>0.484194</td>\n",
       "      <td>0.418704</td>\n",
       "      <td>0.432617</td>\n",
       "      <td>0.484958</td>\n",
       "      <td>0.320647</td>\n",
       "      <td>0.493478</td>\n",
       "      <td>0.384412</td>\n",
       "      <td>0.465014</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>-1.923803e+00</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>-1.022667e+00</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>-1.215302e-01</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>7.796064e-01</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>1.680743e+00</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>8 rows × 25 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "              Meter  Pitch class Bass_A  Pitch class Bass_A#  \\\n",
       "count  5.665000e+03         5665.000000          5665.000000   \n",
       "mean  -1.053605e-15            0.121447             0.055075   \n",
       "std    1.000000e+00            0.326675             0.228147   \n",
       "min   -1.923803e+00            0.000000             0.000000   \n",
       "25%   -1.022667e+00            0.000000             0.000000   \n",
       "50%   -1.215302e-01            0.000000             0.000000   \n",
       "75%    7.796064e-01            0.000000             0.000000   \n",
       "max    1.680743e+00            1.000000             1.000000   \n",
       "\n",
       "       Pitch class Bass_B  Pitch class Bass_C  Pitch class Bass_C#  \\\n",
       "count         5665.000000         5665.000000          5665.000000   \n",
       "mean             0.079259            0.105914             0.038482   \n",
       "std              0.270166            0.307754             0.192374   \n",
       "min              0.000000            0.000000             0.000000   \n",
       "25%              0.000000            0.000000             0.000000   \n",
       "50%              0.000000            0.000000             0.000000   \n",
       "75%              0.000000            0.000000             0.000000   \n",
       "max              1.000000            1.000000             1.000000   \n",
       "\n",
       "       Pitch class Bass_D  Pitch class Bass_D#  Pitch class Bass_E  \\\n",
       "count         5665.000000          5665.000000         5665.000000   \n",
       "mean             0.121624             0.044307            0.119682   \n",
       "std              0.326880             0.205795            0.324618   \n",
       "min              0.000000             0.000000            0.000000   \n",
       "25%              0.000000             0.000000            0.000000   \n",
       "50%              0.000000             0.000000            0.000000   \n",
       "75%              0.000000             0.000000            0.000000   \n",
       "max              1.000000             1.000000            1.000000   \n",
       "\n",
       "       Pitch class Bass_F      ...        Pitch class D  Pitch class D#  \\\n",
       "count         5665.000000      ...          5665.000000     5665.000000   \n",
       "mean             0.081553      ...             0.417476        0.125154   \n",
       "std              0.273707      ...             0.493186        0.330923   \n",
       "min              0.000000      ...             0.000000        0.000000   \n",
       "25%              0.000000      ...             0.000000        0.000000   \n",
       "50%              0.000000      ...             0.000000        0.000000   \n",
       "75%              0.000000      ...             1.000000        0.000000   \n",
       "max              1.000000      ...             1.000000        1.000000   \n",
       "\n",
       "       Pitch class E  Pitch class F  Pitch class F#  Pitch class G  \\\n",
       "count    5665.000000    5665.000000     5665.000000    5665.000000   \n",
       "mean        0.375110       0.226655        0.249250       0.378111   \n",
       "std         0.484194       0.418704        0.432617       0.484958   \n",
       "min         0.000000       0.000000        0.000000       0.000000   \n",
       "25%         0.000000       0.000000        0.000000       0.000000   \n",
       "50%         0.000000       0.000000        0.000000       0.000000   \n",
       "75%         1.000000       0.000000        0.000000       1.000000   \n",
       "max         1.000000       1.000000        1.000000       1.000000   \n",
       "\n",
       "       Pitch class G#  Pitch class A  Pitch class A#  Pitch class B  \n",
       "count     5665.000000    5665.000000     5665.000000    5665.000000  \n",
       "mean         0.116328       0.419241        0.180229       0.316152  \n",
       "std          0.320647       0.493478        0.384412       0.465014  \n",
       "min          0.000000       0.000000        0.000000       0.000000  \n",
       "25%          0.000000       0.000000        0.000000       0.000000  \n",
       "50%          0.000000       0.000000        0.000000       0.000000  \n",
       "75%          0.000000       1.000000        0.000000       1.000000  \n",
       "max          1.000000       1.000000        1.000000       1.000000  \n",
       "\n",
       "[8 rows x 25 columns]"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.describe()\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Всё в порядке, 25 и ожидалось"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Х и У"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "X = data\n",
    "feature_names = X.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "count     5665\n",
       "unique     102\n",
       "top        D_M\n",
       "freq       503\n",
       "Name: Chord label, dtype: object"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "y.describe()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "pandas.core.indexes.base.Index"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(feature_names)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['Meter', 'Pitch class Bass_A', 'Pitch class Bass_A#',\n",
       "       'Pitch class Bass_B', 'Pitch class Bass_C', 'Pitch class Bass_C#',\n",
       "       'Pitch class Bass_D', 'Pitch class Bass_D#', 'Pitch class Bass_E',\n",
       "       'Pitch class Bass_F', 'Pitch class Bass_F#', 'Pitch class Bass_G',\n",
       "       'Pitch class Bass_G#', 'Pitch class C', 'Pitch class C#',\n",
       "       'Pitch class D', 'Pitch class D#', 'Pitch class E', 'Pitch class F',\n",
       "       'Pitch class F#', 'Pitch class G', 'Pitch class G#', 'Pitch class A',\n",
       "       'Pitch class A#', 'Pitch class B'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "feature_names"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(5665, 25)\n",
      "(5665,)\n"
     ]
    }
   ],
   "source": [
    "print(X.shape)\n",
    "print(y.shape)\n",
    "N, d = X.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "pandas.core.frame.DataFrame"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(X)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "pandas.core.series.Series"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3965 1700\n"
     ]
    }
   ],
   "source": [
    "from sklearn.model_selection import train_test_split\n",
    "X_train, X_test, y_train, y_test = train_test_split(X, y, test_size = 0.3, random_state = 42)\n",
    "\n",
    "N_train, _ = X_train.shape \n",
    "N_test,  _ = X_test.shape \n",
    "\n",
    "print(N_train, N_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "5665"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "3965+1700"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Задание №3. Построить  1) random forests, 2) extremely random forests, 3) adaboost, 4) gradient boosted trees.\n",
    "\n",
    "### 1) random forests\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.11778058007566204 0.2723529411764706\n"
     ]
    }
   ],
   "source": [
    "from sklearn import ensemble\n",
    "rf = ensemble.RandomForestClassifier(n_estimators = 1000)\n",
    "rf.fit(X_train, y_train)\n",
    "\n",
    "err_train = np.mean(y_train != rf.predict(X_train))\n",
    "err_test  = np.mean(y_test  != rf.predict(X_test))\n",
    "\n",
    "print(err_train, err_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 100  200  300  400  500  600  700  800  900 1000 1100]\n"
     ]
    }
   ],
   "source": [
    "tree_arr = np.linspace(100, 1100, 11)\n",
    "tree_arr = tree_arr.astype('int')\n",
    "print(tree_arr)\n",
    "test_err = []\n",
    "train_err = []\n",
    "for tree in tree_arr:\n",
    "    rf = ensemble.RandomForestClassifier(n_estimators = tree)\n",
    "    rf.fit(X_train, y_train)\n",
    "\n",
    "    y_train_pred = rf.predict(X_train)\n",
    "    y_test_pred = rf.predict(X_test)\n",
    "    train_err.append(np.mean(y_train != y_train_pred))\n",
    "    test_err.append(np.mean(y_test != y_test_pred))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEaCAYAAADQVmpMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3X+UVeV97/H3Z/hZgiIC5iojMyQxCagU5UA0NhoTf0CsoDEaFZJ4YyRNam/bLF3RS5JbaWmtNIlJNY2Y2jQBtWhiLia6QL0QkyxpGNQYfmhECzKSCqJEEUSB7/1j79Ezw5mZM/ucPWdm+LzWOmvOfvbzPPvZ82N/5tn7nH0UEZiZmWVRV+sBmJlZ7+UQMTOzzBwiZmaWmUPEzMwyc4iYmVlmDhEzM8vMIWJmZpk5RKzPkLRR0m5JO4seN9V6XJWS1CgpJPWv9VjM2vIvpfU150bEg51VktQ/IvZ2VtbVPmqlJ43FDi6eidhBQdJlkn4l6ZuSXgL+pp2yOklfkbRJ0lZJP5A0LO2jZUZwuaTngP9XYjvrJf1p0XJ/SS9KOlHSYEkLJW2XtEPSKknvLGP4D6dfd6Szq5NLjT3d3mfTMbwsaamkhqKxvF/SA5JekvSUpIsyf0PNUg4RO5h8AHgWOAKY107ZZenjdOBdwFCg7Smx04BxwNkltnEHcEnR8tnAixHxKPAZYBhwNDAC+DNgdxnjPjX9elhEDI2IR0qNXdJ5wP8GPg6MAn6RjgdJ7wAeAG5P618CfEfSsWVs36xdDhHra36S/pff8riiaN2WiPjniNgbEbvbKZsJfCMino2IncC1wMVtrkf8TUS8VtRHsduB6ZKGpMuXpmUAb5KEx3siYl9ErI6IVyrY17Zj/zzwDxGxPj219ffAxHQ28qfAxoj4t7T+o8CPgE9UsH0zXxOxPue8Dq6JbC6j7ChgU9HyJpK/k+LTTqX6ASAiNkhaD5wr6V5gOnBCuvqHJLOQOyUdBiwE5kTEm+3114m242gAviXp60VlAkan6z4gaUfRuv7pmMwyc4jYwaTULavblm0hOeC2GAPsBV4A6jvop1jLKa06YF1EbABIw+I64DpJjcB9wFPAv2YYd6nyzcC8iFjUtmI6G/l5RJzZybbMusSns8xauwP4a0ljJQ0lOSX0H1185dOdwFnAF3j7VBaSTpd0vKR+wCskp7f2ldHfNmA/yTWajnwXuLblOoekYZIuTNf9FHivpE9JGpA+Jksa14X9MjuAQ8T6mnvbvE/kni62v43kFM/DwH8BrwN/0ZUOIuL3wCPAB4H/KFr1P4C7SQJkPfBzklNaSPqupO+2098ukov+v0qv85zUTr17gH8kOV32CrAGmJaue5Uk2C4mmW39d1p3UFf2zawt+UOpzMwsK89EzMwsM4eImZll5hAxM7PMHCJmZpaZQ8TMzDLL9c2GkqYC3wL6Ad+LiOvbrP8S8DmSN3NtAz4bEZvSdTcA55AE3QPAX0YHLyUbOXJkNDY25rEbZmZ91urVq1+MiFFZ2+cWIukbqm4GzgSagVWSlkTEuqJqjwGFiNgl6QvADcAnJX0QOAWYkNb7JclN71a0t73GxkaampqqvyNmZn2YpE2d12pfnqezpgAb0hvZvUHyLt4ZxRUiYnn6RiqAlbS+rcRgYCDJm6EGkNx2wszMepA8Q2Q0rW8Q15yWtedy4H6A9FbXy4Hfp4+lEbE+p3GamVlGeYaISpSVvKYhaRZQAOany+8h+byGepLg+YikU0u0my2pSVLTtm3bqjZwMzMrT54h0kxy2+sW9ST37GlF0hnAHGB6ROxJi88HVkbEzvQzHe4HDrhfUEQsiIhCRBRGjcp8XcjMzDLKM0RWAcekd0MdSHLjtyXFFSSdANxCEiBbi1Y9B5yWfrToAJKL6j6dZWZWLYsWQWMjk2BSJd3kFiLprbOvBJaSBMDiiFgraa6k6Wm1+SQfP3qXpMcltYTM3cAzwG+B3wC/iYh7qzrA9BtIXV3yddEBH8FgZtY3LVoEs2fDpopemAX0obv4FgqFKPslvi3fwF273i4bMgQWLICZM/MZoJlZniJg3z54440DH2++2Xr5/PPhheQFrwWgKaLUNeyy9J0QkaKpoQHmzes4CHbuhGOOgf/+7wPXNTTAxo25jfEtixbBnDnw3HMwZkznYzbra3rL38C+fQcegNs7MFdjXaX9ZjieVxoifevjcTdtSmYYADNmwPr1sHYtrFuXfF27tuPp26ZN8NvfwnHHgTJ/TzvWdhZUPOae+EdkVm2l/gY+9zl4/nk4++x8Dq5ZD9r79+fzPairg4EDSz8GDGi9PHgwHHpo6XUdteto3ac/DVu3dj7OMvStmUjLQr9+yX8QLQYOhPe/H449Nnl861vQ0UuCGxrg3HOTx2mnwaAKPvxt//5kdvPEE8nj+uth9+4D6x11VPJHZNZbRMCrr8LLL8OOHcnXcp4/9VQ+B+dBg7IfaKt1cC53W/36VX//u6IoyD0TKWXfPvi7v0sCY/x4eNe7oH/RrjY2lr4mcsMNyQ/53nvhX/8VbroJhg5N/js691z42Mdg2bL2p+E7diQzmSeeaP11585kvdT+dHPLFnjve5NtnX02fPjDybZrobtONfSWUxq10F3fm7174Q9/ePtAX24gtCx3FAZ1dXDYYclj+PDkcdRRyRmCUiT48Y+zHaD79cvv7EFf1PK7NGdOxRfX++ZMpJxrG539ke7eDQ89lATKT3+aHOQh+cMo/sPp3z85/fXSS0lfLYYPhz/+Yzj+eJgwIXm0zIRK/dCGD4eTT4YVK5JwGzAA/uRP3g6VCROSbeet1IsOBg9Ovj/Tp7ffrquWLEm+/6+//naZX9yQ6OoLP15/vesh0PL81Vc7HsvAgW8HwPDhrQOhs+eHHFL6d7axsfTfQHddk7RWJK2OiELm9n0uRPI4EEXAo4/CRz+a/NfW1oABcOGFb4fFhAnJf1yl/jPq7ACxZw/88pewdGnyeOKJpM473wlnnZUEyplnwhFHvN1fuf+xRsArrySBuGVLcvqs+OuWLbBqVetTgd2tf3+YOBFGjIDDD2/9KFU2fHjrWWa5ajULikjOwe/enfwO7N7d+vmuXcn56lKnW4cMSWaobQNhz54D6xYbOvTAA325gTB4cPX/w/erI3sUh0iq7FdnVaKurvTpKKlr53i7cgD7/e+TU2hLl8IDD8CLLyblJ54I9fVJefFBZNAguOIKePe7S4fEa68duI1hw2D06CT4Hnyw/XH/8Ifl72NnPvWp9tdNnZrM7FoeL7/c8atOhg3rOGzaLi9fDldddeBB7JvfhHPO6fgAX42ySkJ60qSOD/ptl4cNS/7J6Wl8KrPHcIikuvQ+kaxqPQ3fvz+ZEbXMUn7xi47rDx6cBENLQJR6fuSR8I53vN2mu/axK9vZvz/5z7s4WF56CbZv73j55Zfze3UNJP88/NEfJY8hQ0p/zVJ2wQXJPw/lfG/MKlRpiPTNC+t5mTev9DR83rzu2X5dHRQKyWPOnI5nRtu3J/+RdvVURHftY1e2U1f39iyiK/bvT04/tg2fSy9tv80tt5R/sB80KJ+LufPn1/b3zKwrIqJPPCZNmhTdYuHCiIaGCCn5unBh92y3lIaGiCRGWj8aGirrt7v2sVbfy7y+b9XUk37PrE8DmqKCY69PZ/VmvkCZjb9vZm+p9HRWN7xm1HIzc2Zy4GtoSE6rNDT4QFgOf9/MqsYzETOzg5hnImZmVjMOETMzy8whYmZmmeUaIpKmSnpK0gZJ15RY/yVJ6yQ9IekhSQ1F68ZIWiZpfVqnMc+xmplZ1+UWIpL6ATcD04DxwCWSxrep9hhQiIgJJB+Je0PRuh8A8yNiHDAFqM7N783MrGrynIlMATZExLMR8QZwJzCjuEJELI+IlhfrrwTqAdKw6R8RD6T1dhbVMzOzHiLPEBkNbC5abk7L2nM5cH/6/L3ADkk/lvSYpPnpzMbMzHqQPEOk1E2FSr4pRdIsko/6nZ8W9Qc+BFwFTAbeBVxWot1sSU2SmrZ19EmFZmaWizxDpBk4umi5HtjStpKkM4A5wPSI2FPU9rH0VNhe4CfAiW3bRsSCiChERGHUqFFV3wEzM+tYniGyCjhG0lhJA4GLgSXFFSSdANxCEiBb27QdLqklGT4CrMtxrGZmlkFuIZLOIK4ElgLrgcURsVbSXEktn7M6HxgK3CXpcUlL0rb7SE5lPSTptySnxm7Na6xmZpaN751lZnYQ872zzMysZhwiZmaWmUPEzMwyc4iYmVlmDhEzM8vMIWJmZpk5RMzMLDOHiJmZZeYQMTOzzBwiZmaWmUPEzMwyc4iYmVlmDhEzM8vMIWJmZpk5RMzMLDOHiJmZZZZriEiaKukpSRskXVNi/ZckrZP0hKSHJDW0WX+opOcl3ZTnOM3MLJvcQkRSP+BmYBowHrhE0vg21R4DChExAbgbuKHN+r8Ffp7XGM3MrDJ5zkSmABsi4tmIeAO4E5hRXCEilkfErnRxJVDfsk7SJOCdwLIcx2hmZhXIM0RGA5uLlpvTsvZcDtwPIKkO+DpwdW6jMzOzivXPsW+VKIuSFaVZQAE4LS36InBfRGyWSnXzVrvZwGyAMWPGVDRYMzPrujxDpBk4umi5HtjStpKkM4A5wGkRsSctPhn4kKQvAkOBgZJ2RkSri/MRsQBYAFAoFEoGlJmZ5SfPEFkFHCNpLPA8cDFwaXEFSScAtwBTI2JrS3lEzCyqcxnJxfcDXt1lZma1lds1kYjYC1wJLAXWA4sjYq2kuZKmp9Xmk8w07pL0uKQleY3HzMyqTxF94yxQoVCIpqamWg/DzKxXkbQ6IgpZ2/sd62ZmlplDxMzMMnOImJlZZg4RMzPLzCFiZmaZOUTMzCwzh4iZmWXmEDEzs8wcImZmlplDxMzMMnOImJlZZg4RMzPLzCFiZmaZOUTMzCwzh4iZmWXmEDEzs8xyDRFJUyU9JWmDpAM+3lbSlyStk/SEpIckNaTlEyU9Imltuu6TeY7TzMyyyS1EJPUDbgamAeOBSySNb1PtMZLPT58A3A3ckJbvAj4dEccCU4EbJR2W11jNzCybPGciU4ANEfFsRLwB3AnMKK4QEcsjYle6uBKoT8t/FxFPp8+3AFuBUTmO1czMMsgzREYDm4uWm9Oy9lwO3N+2UNIUYCDwTFVHZ2ZmFeufY98qURYlK0qzgAJwWpvyI4EfAp+JiP0l2s0GZgOMGTOm0vGamVkX5TkTaQaOLlquB7a0rSTpDGAOMD0i9hSVHwr8DPhKRKwstYGIWBARhYgojBrls11mZt0tzxBZBRwjaaykgcDFwJLiCpJOAG4hCZCtReUDgXuAH0TEXTmO0czMKpBbiETEXuBKYCmwHlgcEWslzZU0Pa02HxgK3CXpcUktIXMRcCpwWVr+uKSJeY3VzMyyUUTJyxS9TqFQiKamploPw8ysV5G0OiIKWdv7HetmZpaZQ8TMzDJziJiZWWYOETMzy8whYmZmmTlEzMwsM4eImZll5hAxM7PM8rwBo5lZj/fmm2/S3NzM66+/Xuuh5Grw4MHU19czYMCAqvbrEDGzg1pzczOHHHIIjY2NSKVuPt77RQTbt2+nubmZsWPHVrVvn84ys4Pa66+/zogRI/psgABIYsSIEbnMthwiZnbQ68sB0iKvfXSImJnV0I4dO/jOd77T5XYf+9jH2LFjRw4j6hqHiJlZFyxaBI2NUFeXfF20qLL+2guRffv2ddjuvvvu47DDDqts41XQaYhI6ifpr7tjMGZmPdmiRTB7NmzaBBHJ19mzKwuSa665hmeeeYaJEycyefJkTj/9dC699FKOP/54AM477zwmTZrEsccey4IFC95q19jYyIsvvsjGjRsZN24cV1xxBcceeyxnnXUWu3fvrnRXy9bpq7MiYp+kGcA3u2E8ZmY181d/BY8/3v76lSthz57WZbt2weWXw623lm4zcSLceGP7fV5//fWsWbOGxx9/nBUrVnDOOeewZs2at15Fddttt3H44Yeze/duJk+ezAUXXMCIESNa9fH0009zxx13cOutt3LRRRfxox/9iFmzZpWzyxUr9yW+v5J0E/AfwGsthRHxaC6jMjPrgdoGSGflWUyZMqXVy3C//e1vc8899wCwefNmnn766QNCZOzYsUycmHz466RJk9i4cWP1BtSJckPkg+nXuUVlAXyko0aSpgLfAvoB34uI69us/xLwOWAvsA34bERsStd9BvhKWvXvIuLfyxyrmVkmHc0YILkGsmnTgeUNDbBiRXXG8I53vOOt5ytWrODBBx/kkUceYciQIXz4wx8u+TLdQYMGvfW8X79+Pet0FkBEnN7VjiX1A24GzgSagVWSlkTEuqJqjwGFiNgl6QvADcAnJR0O/B+gQBJWq9O2L3d1HGZm1TJvXnINZNeut8uGDEnKszrkkEN49dVXS677wx/+wPDhwxkyZAhPPvkkK1euzL6hnJT16ixJwyR9Q1JT+vi6pGGdNJsCbIiIZyPiDeBOYEZxhYhYHhEtP46VQH36/GzggYh4KQ2OB4Cp5e6UmVkeZs6EBQuSmYeUfF2wICnPasSIEZxyyikcd9xxXH311a3WTZ06lb179zJhwgS++tWvctJJJ1W4B9VX7ums24A1wEXp8qeAfwM+3kGb0cDmouVm4AMd1L8cuL+DtqPbNpA0G5gNMGbMmA66NjOrjpkzKwuNUm6//faS5YMGDeL+++8vua7lusfIkSNZs2bNW+VXXXVVdQfXiXJD5N0RcUHR8nWSOngNAwCl3h4ZJStKs0hOXZ3WlbYRsQBYAFAoFEr2bWZm+Sn3zYa7Jf1Jy4KkU4DOrtw0A0cXLdcDW9pWknQGMAeYHhF7utLWzMxqq9yZyJ8BPyi6DvIy8JlO2qwCjpE0FngeuBi4tLiCpBOAW4CpEbG1aNVS4O8lDU+XzwKuLXOsZmbWTToNEUl1wPsi4o8lHQoQEa901i4i9kq6kiQQ+gG3RcRaSXOBpohYAswHhgJ3pTcHey4ipkfES5L+liSIAOZGxEtZdtDMzPJTzjvW96dhsLic8GjT9j7gvjZlXyt6fkYHbW8juaBvZmY9VLnXRB6QdJWkoyUd3vLIdWRmZtbjlRsinwX+HHgYWJ0+mvIalJnZwSLrreABbrzxRnYVv/OxBsq5i28dMCsixrZ5vKsbxmdm1rNU+V7wvT1Eyr0m8k/Ayd0wHjOznqvlXvAtB+6We8FD5ncgFt8K/swzz+SII45g8eLF7Nmzh/PPP5/rrruO1157jYsuuojm5mb27dvHV7/6VV544QW2bNnC6aefzsiRI1m+fHmVdrJryn2J7zJJFwA/jgi/qc/M+qYa3Au++Fbwy5Yt4+677+bXv/41EcH06dN5+OGH2bZtG0cddRQ/+9nPgOSeWsOGDeMb3/gGy5cvZ+TIkV3d06op95rIl4DFwB5Jr0h6VVKXXqllZtbr5Xwv+GXLlrFs2TJOOOEETjzxRJ588kmefvppjj/+eB588EG+/OUv84tf/IJhwzq7dWH3KXcmMgyYCYyNiLmSxgBH5jcsM7MaqPG94COCa6+9ls9//vMHrFu9ejX33Xcf1157LWeddRZf+9rXSvTQ/cqdidwMnARcki6/CtyUy4jMzHqqefOSe78Xq/Be8MW3gj/77LO57bbb2LlzJwDPP/88W7duZcuWLQwZMoRZs2Zx1VVX8eijjx7QtlbKnYl8ICJOlPQYQES8LGlgjuMyM+t5Wi6ez5kDzz0HY8YkAVLBbX2LbwU/bdo0Lr30Uk4+OXkd09ChQ1m4cCEbNmzg6quvpq6ujgEDBvAv//IvAMyePZtp06Zx5JFH1uzCusq5Ti7pP0k+3XBVGiajgGURcULeAyxXoVCIpia/dcXMumb9+vWMGzeu1sPoFqX2VdLqiChk7bPc01nfBu4BjpA0D/gl8PdZN2pmZn1DuR+Pu0jSauCjJJ/1cV5ErM91ZGZm1uOVe02EiHgSeDLHsZiZWS9T7uksM7M+62B4D3Ve++gQMbOD2uDBg9m+fXufDpKIYPv27QwePLjqfZd9OsvMrC+qr6+nubmZbdu21XoouRo8eDD19fVV7zfXEJE0FfgWyScbfi8irm+z/lTgRmACcHFE3F207gbgHJLZ0gPAX/q+XWZWbQMGDGDs2LG1HkavldvpLEn9SN7pPg0YD1wiaXybas8BlwG3t2n7QeAUknA5DpgMnJbXWM3MLJs8ZyJTgA0R8SyApDuBGcC6lgoRsTFdt79N2wAGAwNJXlI8AHghx7GamVkGeV5YHw1sLlpuTss6FRGPAMuB36ePpaXelyJptqQmSU19/XymmVlPlGeIqERZWdc0JL0HGAfUkwTPR9LrJ607i1gQEYWIKIwaNaqiwZqZWdflGSLNwNFFy/XAljLbng+sjIidEbETuJ/kLsJmZtaD5Bkiq4BjJI1N7/h7MbCkzLbPAadJ6i9pAMlFdd9mxcysh8ktRCJiL3AlsJQkABZHxFpJcyVNB5A0WVIzcCFwi6S1afO7gWeA3wK/AX4TEffmNVYzM8umrFvB9wa+FbyZWdd1163gzczMDuAQMTOzzBwiZmaWmUPEzMwyc4iYmVlmDhEzM8vMIWJmZpk5RMzMLDOHiJmZZeYQMTOzzBwiZmaWmUPEzMwyc4iYmVlmDhEzM8vMIWJmZpk5RMzMLLNcQ0TSVElPSdog6ZoS60+V9KikvZI+0WbdGEnLJK2XtE5SY55jNTOzrsstRCT1A24GpgHjgUskjW9T7TngMuD2El38AJgfEeOAKcDWvMZqZmbZ9M+x7ynAhoh4FkDSncAMYF1LhYjYmK7bX9wwDZv+EfFAWm9njuM0M7OM8jydNRrYXLTcnJaV473ADkk/lvSYpPnpzKYVSbMlNUlq2rZtWxWGbGZmXZFniKhEWZTZtj/wIeAqYDLwLpLTXq07i1gQEYWIKIwaNSrrOM3MLKM8Q6QZOLpouR7Y0oW2j0XEsxGxF/gJcGKVx2dmZhXKM0RWAcdIGitpIHAxsKQLbYdLaplefISiaylmZtYz5BYi6QziSmApsB5YHBFrJc2VNB1A0mRJzcCFwC2S1qZt95GcynpI0m9JTo3dmtdYzcwsG0WUe5miZysUCtHU1FTrYZiZ9SqSVkdEIWt7v2PdzMwyc4iYmVlmDhEzM8vMIWJmZpk5RMzMLDOHiJmZZeYQMTOzzBwiZmaWmUPEzMwyc4iYmVlmDhEzM8vMIWJmZpk5RMzMLDOHiJmZZeYQMTOzzBwiZmaWWa4hImmqpKckbZB0TYn1p0p6VNJeSZ8osf5QSc9LuinPcZqZWTa5hYikfsDNwDRgPHCJpPFtqj0HXAbc3k43fwv8PK8xmplZZfKciUwBNkTEsxHxBnAnMKO4QkRsjIgngP1tG0uaBLwTWJbjGM3MrAJ5hshoYHPRcnNa1ilJdcDXgas7qTdbUpOkpm3btmUeqJmZZZNniKhEWZTZ9ovAfRGxuaNKEbEgIgoRURg1alSXB2hmZpXpn2PfzcDRRcv1wJYy254MfEjSF4GhwEBJOyPigIvzZmZWO3mGyCrgGEljgeeBi4FLy2kYETNbnku6DCg4QMzMep7cTmdFxF7gSmApsB5YHBFrJc2VNB1A0mRJzcCFwC2S1uY1HjMzqz5FlHuZomcrFArR1NRU62GYmfUqklZHRCFre79j3czMMnOImJlZZg4RMzPLzCFiZmaZOUTMzCwzh4iZmWXmEDEzs8wcImZmlplDxMzMMnOImJlZZg4RMzPLzCFiZmaZOUTMzCwzh4iZmWXmEDEzs8wcImZmllmuISJpqqSnJG2QdMDH20o6VdKjkvZK+kRR+URJj0haK+kJSZ/Mc5xmZpZNbiEiqR9wMzANGA9cIml8m2rPAZcBt7cp3wV8OiKOBaYCN0o6LK+xmplZNv1z7HsKsCEingWQdCcwA1jXUiEiNqbr9hc3jIjfFT3fImkrMArYkeN4zcysi/I8nTUa2Fy03JyWdYmkKcBA4JkS62ZLapLUtHo1NDbCokXl9btoUVK/rq572lW7jzz6yrPPWm6np227HD19fJXq6/vXG7T8DGDSpIo6iohcHsCFwPeKlj8F/HM7db8PfKJE+ZHAU8BJnW9vUkDEkCERCxdGhxYuTOrB248821W7jzz6yrPPWm6np227HD19fJXq6/vXG7T+GUyKqOBYr0gO1lUn6WTgbyLi7HT52jS0/qFE3e8DP42Iu4vKDgVWAP8QEXd1vr1CQBMAgwbBSSe1X3flStiz58DyvNpVu488+sqzz1pup6dtuxw9fXyV6uv71xu0/hkUiGhS1r7yPJ21CjhG0lhJA4GLgSXlNEzr3wP8oJwAaavUL2g56/NqV+0+8ugrzz5ruZ2etu1y9PTxVaqv719vUNXvdSXTmM4ewMeA35Fcz5iTls0FpqfPJ5NcK3kN2A6sTctnAW8Cjxc9Jna8rUlvTY0bGjqeyjU0tJ5K592u2n3k0VeefdZyOz1t2+Xo6eOrVF/fv96g9c+gstNZuYZIdz58TaSyvvLss5bb6WnbLkdPH1+l+vr+9QbVvCZS84N/tR4wKRoayv9FXLgwSWMpuqVdtfvIo688+6zldnratsvR08dXqb6+f71By8+gx15Y726FQiGamppqPQwzs15F0uqIKGRt73tnmZlZZg4RMzPLzCFiZmaZOUTMzCwzh4iZmWXWZ16dJWkbsClD05HAixnaDQP+kKFdtfvIo688+6zldnratsvR08dXqb6+f73B+yLikKyN87wVfLeKiFFZ2klqyvLyNkkLImJ2lm1Ws488+sqzz1pup6dtuxw9fXyV6uv71xtIqui9ET6dld29PaSPPPrKs89abqenbbscPX18lerr+9fn9ZnTWVllnYmYmfUFlR4DPROBBbUegJlZDVV0DDzoZyJmZpadZyJmZpaZQ8TMzDJziJiZWWYOkSKSxkn6rqS7JX2h1uMxM+tuks6TdKuk/yvprM7q9/kQkXSbpK2S1rQpnyrpKUkbJF0DEBHrI+LPgIsAv+zXzPqELh4HfxIRVwCXAZ/srO8+HyLA94GpxQWS+gE3A9OA8cAlksan66YDvwQe6t5hmpnl5vt04TiY+kpIiOIFAAACW0lEQVS6vkN9PkQi4mHgpTbFU4ANEfFsRLwB3AnMSOsviYgPAjO7d6RmZvnoynFQiX8E7o+IRzvru8/cO6uLRgObi5abgQ9I+jDwcWAQcF8NxmVm1l1KHgeBvwDOAIZJek9EfLejTg7WEFGJsoiIFcCK7h2KmVlNtHcc/Dbw7XI76fOns9rRDBxdtFwPbKnRWMzMaqEqx8GDNURWAcdIGitpIHAxsKTGYzIz605VOQ72+RCRdAfwCPA+Sc2SLo+IvcCVwFJgPbA4ItbWcpxmZnnJ8zjoGzCamVlmfX4mYmZm+XGImJlZZg4RMzPLzCFiZmaZOUTMzCwzh4iZmWXmEDGrgKTDJH2x1uMwqxWHiFllDgMOCJH0Nttmfd7BegNGs2q5Hni3pMeBN4GdwO+BicB4SbOA/wUMBP4T+GJE7Es/Me46kjtGPwP8z4jYWYsdMKuE37FuVgFJjcBPI+K49KMEfgYcFxH/JWkccAPw8Yh4U9J3gJUkHzPwY2BaRLwm6cvAoIiYW5OdMKuAZyJm1fXriPiv9PlHgUnAKkkAfwRsBU4i+SS5X6XlA0nua2TW6zhEzKrrtaLnAv49Iq4triDpXOCBiLikW0dmlgNfWDerzKvAIe2sewj4hKQjACQdLqmB5JTWKZLek5YPkfTebhmtWZV5JmJWgYjYLulXktYAu4EXitatk/QVYJmkOpIL738eESslXQbcIWlQWv0rwO+6efhmFfOFdTMzy8yns8zMLDOHiJmZZeYQMTOzzBwiZmaWmUPEzMwyc4iYmVlmDhEzM8vMIWJmZpn9f3tXcjnqsyF2AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.semilogx(tree_arr, train_err, 'b-o', label = 'train')\n",
    "plt.semilogx(tree_arr, test_err, 'r-o', label = 'test')\n",
    "plt.xlim([np.max(tree_arr), np.min(tree_arr)])\n",
    "plt.title('Error vs. tree')\n",
    "plt.xlabel('tree')\n",
    "plt.ylabel('error')\n",
    "plt.legend()\n",
    "pass"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 2) extremely random forests"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.11778058007566204 0.2776470588235294\n"
     ]
    }
   ],
   "source": [
    "ert = ensemble.ExtraTreesClassifier(n_estimators = 100).fit(X_train, y_train)\n",
    "\n",
    "err_train = np.mean(y_train != ert.predict(X_train))\n",
    "err_test  = np.mean(y_test  != ert.predict(X_test))\n",
    "\n",
    "print(err_train, err_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 100  200  300  400  500  600  700  800  900 1000 1100]\n"
     ]
    }
   ],
   "source": [
    "tree_arr_ert = np.linspace(100, 1000, 10)\n",
    "tree_arr_ert = tree_arr_ert.astype('int')\n",
    "print(tree_arr_ert)\n",
    "test_err = []\n",
    "train_err = []\n",
    "for tree in tree_arr_ert:\n",
    "    ert = ensemble.ExtraTreesClassifier(n_estimators = tree)\n",
    "    ert.fit(X_train, y_train)\n",
    "\n",
    "    y_train_pred = rf.predict(X_train)\n",
    "    y_test_pred = rf.predict(X_test)\n",
    "    train_err.append(np.mean(y_train != y_train_pred))\n",
    "    test_err.append(np.mean(y_test != y_test_pred))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 77,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEaCAYAAADQVmpMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAHrhJREFUeJzt3X+8VXWd7/HXm9+hhoroNY4CTWbgj4uyIc2bZuMPqAmtzERpYvJG5Tj3MePDHunVmqRxauTWOI6W4lyncUANLRsqfIg6kNVDJg5Kxg/NI6Ec6coRI0UUBT/3j72ObQ6bc/b57r3O3mfzfj4e+8FZ3/X9ftdncXC/XWvttbYiAjMzsxQD6l2AmZn1Xw4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RaxqSNkh6VdK2kteN9a6rWpLGSgpJg+pdi1lX/kdpzeYjEfFgT50kDYqInT219XaOemmkWmzf4iMR2ydImiXpF5L+UdKLwFf30jZA0tWSnpG0WdLtkkZkc3QeEVws6VngP8tsZ52kPytZHiTpBUknShomab6kLZK2Sloh6bAKyn84+3NrdnR1crnas+19Jqvh95LulzSmpJb3SHpA0ouSnpR0fvJfqFnGIWL7kvcC64FDgWv30jYre50OvBPYH+h6Suw0YDxwdplt3AnMKFk+G3ghIh4FPg2MAI4ARgKfB16toO5Tsz8PjIj9I+KRcrVLOhf438DHgFHAz7J6kLQf8ABwR9Z/BvBtScdUsH2zvXKIWLP5YfZ/+Z2vz5as2xQR/xwROyPi1b20XQR8KyLWR8Q24Erggi7XI74aEa+UzFHqDmC6pOHZ8oVZG8AbFMPjXRGxKyJWRsRLVexr19o/B3w9ItZlp7b+HpiYHY38GbAhIv416/8o8H3gvCq2b+ZrItZ0zu3mmsjGCtreATxTsvwMxf9OSk87lZsHgIhok7QO+IikHwHTgROy1f9O8SjkLkkHAvOBqyLijb3N14OudYwB/knSN0vaBIzO1r1X0taSdYOymsySOURsX1LukdVd2zZRfMPtdCSwE3geaOlmnlKdp7QGAGsjog0gC4trgGskjQUWA08C/zeh7nLtG4FrI2JB147Z0chPI+LMHrZl1is+nWW2uzuBv5E0TtL+FE8Jfa+Xn3y6CzgL+AJ/PJWFpNMlHSdpIPASxdNbuyqYrwN4k+I1mu7cDFzZeZ1D0ghJn8jW/Rh4t6RPSRqcvSZLGt+L/TLbg0PEms2Putwncm8vx99G8RTPw8BvgdeAv+rNBBHxO+AR4H3A90pW/TfgHooBsg74KcVTWki6WdLNe5lvO8WL/r/IrvOctJd+9wL/QPF02UvAamBatu5lisF2AcWjrf+X9R3am30z60r+UiozM0vlIxEzM0vmEDEzs2S5hoikqdmdsW2Sriiz/jJJayU9LumhLnfXXidpTXb37Q2SlGetZmbWe7mFSPYJlJsoXtibAMyQNKFLt8eAQkQcT/GC43XZ2PcBpwDHA8cCkyneJWxmZg0kzyORKUBbdufv6xQ/9nhOaYeIWJp98gRgObt/Dn8YMITip0cGU/ycvpmZNZA8bzYcze531LZTfNbP3lwM3AcQEY9IWgr8juIdtzdGxLquAyTNBmYD7LfffpPe85731Kh0M7N9w8qVK1+IiFGp4/MMkXLXMMp+nljSTKBAdspK0rsoPuCu88jkAUmnRsTDpeMiYh4wD6BQKERra2uNSjcz2zdIeqbnXnuX5+msdorPCerUQvEmp91IOgO4CpgeETuy5o8CyyNiW/YQvPuAsjdYmZlZ/eQZIiuAo7LHRwyheKfsotIOkk4AbqEYIJtLVj0LnJZ9F8Ngikcoe5zOMjOz+sotRLJnDV0K3E8xABZGxBpJcyRNz7rNpfh9DXdLWiWpM2TuAZ4Gfg38CvhVRPwor1rNzCxN0zz2xNdEzMx6T9LKiCikjm+eO9ZXroSxY2HBHk/BLm/BgmL/AQP6Zlyt58hjrjznrOd2Gm3blWj0+qrV7PvXH2S/g0kwqap5IqIpXpMgAiKGD4+YPz+6NX9+sV/nmLzH1XqOPObKc856bqfRtl2JRq+vWs2+f/1Bye9gEkRU8d7bPKezpHjrZNbQoXBSNx/mWr4cduzYsz2vcbWeI4+58pyznttptG1XotHrq1az719/UPI7KACtEcmPlWqe01mlyv0DrWR9XuNqPUcec+U5Zz2302jbrkSj11etZt+//qCWf9fVHMY00mtS6aHxmDHdH8qNGbP7oXTe42o9Rx5z5TlnPbfTaNuuRKPXV61m37/+oOR3UO3prOY7Ehk+HK69tvs+115b7NdX42o9Rx5z5TlnPbfTaNuuRKPXV61m37/+oNzvIFU1CdRIr0md/ydT6cW5+fOL/aW+GVfrOfKYK88567mdRtt2JRq9vmo1+/71B9nvwBfWM75PxMys93yfiJmZ1Y1DxMzMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMwsWa4hImmqpCcltUm6osz6yyStlfS4pIckjSlZd6SkJZLWZX3G5lmrmZn1Xm4hImkgcBMwDZgAzJA0oUu3x4BCRBxP8XvVrytZdzswNyLGA1OAzXnVamZmafI8EpkCtEXE+oh4HbgLOKe0Q0QsjYjt2eJyoAUgC5tBEfFA1m9bST8zM2sQeYbIaGBjyXJ71rY3FwP3ZT+/G9gq6QeSHpM0Nzuy2Y2k2ZJaJbV2dHTUrHAzM6tMniFS7usWyz4yWNJMit/SODdrGgS8H7gcmAy8E5i1x2QR8yKiEBGFUaNG1aJmMzPrhTxDpB04omS5BdjUtZOkM4CrgOkRsaNk7GPZqbCdwA+BE3Os1czMEuQZIiuAoySNkzQEuABYVNpB0gnALRQDZHOXsQdJ6jy8+CCwNsdazcwsQW4hkh1BXArcD6wDFkbEGklzJE3Pus0F9gfulrRK0qJs7C6Kp7IekvRriqfGbs2rVjMzS+NvNjQz24f5mw3NzKxuHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlizXEJE0VdKTktokXVFm/WWS1kp6XNJDksZ0Wf92Sc9JujHPOs3MLE1uISJpIHATMA2YAMyQNKFLt8eAQkQcD9wDXNdl/deAn+ZVo5mZVSfPI5EpQFtErI+I14G7gHNKO0TE0ojYni0uB1o610maBBwGLMmxRjMzq0KeITIa2Fiy3J617c3FwH0AkgYA3wS+2N0GJM2W1CqptaOjo8pyzcyst/IMEZVpi7IdpZlAAZibNV0CLI6IjeX6vzVZxLyIKEREYdSoUVUVa2ZmvTcox7nbgSNKlluATV07SToDuAo4LSJ2ZM0nA++XdAmwPzBE0raI2OPivJmZ1U+eIbICOErSOOA54ALgwtIOkk4AbgGmRsTmzvaIuKikzyyKF98dIGZmDSa301kRsRO4FLgfWAcsjIg1kuZImp51m0vxSONuSaskLcqrHjMzqz1FlL1M0e8UCoVobW2tdxlmZv2KpJURUUgd7zvWzcwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLFmuISJpqqQnJbVJ2uM70iVdJmmtpMclPSRpTNY+UdIjktZk6z6ZZ51mZpYmtxCRNBC4CZgGTABmSJrQpdtjQCEijgfuAa7L2rcDfx4RxwBTgeslHZhXrWZmlibPI5EpQFtErI+I14G7gHNKO0TE0ojYni0uB1qy9t9ExFPZz5uAzcCoHGs1M7MEeYbIaGBjyXJ71rY3FwP3dW2UNAUYAjxdZt1sSa2SWjs6Oqos18zMeivPEFGZtijbUZoJFIC5XdoPB/4d+IuIeHOPySLmRUQhIgqjRvlAxcysrw3Kce524IiS5RZgU9dOks4ArgJOi4gdJe1vB34CXB0Ry3Os08zMEuV5JLICOErSOElDgAuARaUdJJ0A3AJMj4jNJe1DgHuB2yPi7hxrNDOzKuQWIhGxE7gUuB9YByyMiDWS5kiannWbC+wP3C1plaTOkDkfOBWYlbWvkjQxr1rNzCyNIspepuh3CoVCtLa21rsMM7N+RdLKiCikjvcd62ZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWLM8HMJqZNbw33niD9vZ2XnvttXqXkqthw4bR0tLC4MGDazqvQ8TM9mnt7e0ccMABjB07FqncN1j0fxHBli1baG9vZ9y4cTWd26ezzGyf9tprrzFy5MimDRAASYwcOTKXoy2HiJnt85o5QDrltY8OETOzOtq6dSvf/va3ez3uQx/6EFu3bs2hot5xiJiZ9cKCBTB2LAwYUPxzwYLq5ttbiOzatavbcYsXL+bAAw+sbuM10GOISBoo6W/6ohgzs0a2YAHMng3PPAMRxT9nz64uSK644gqefvppJk6cyOTJkzn99NO58MILOe644wA499xzmTRpEscccwzz5s17a9zYsWN54YUX2LBhA+PHj+ezn/0sxxxzDGeddRavvvpqtbtasR4/nRURuySdA/xjbyeXNBX4J2Ag8C8R8Y0u6y8D/iewE+gAPhMRz2TrPg1cnXX9u4j4t95u38ysN/76r2HVqr2vX74cduzYvW37drj4Yrj11vJjJk6E66/f+5zf+MY3WL16NatWrWLZsmV8+MMfZvXq1W99iuq2227j4IMP5tVXX2Xy5Ml8/OMfZ+TIkbvN8dRTT3HnnXdy6623cv755/P973+fmTNnVrLLVav0I76/kHQj8D3glc7GiHh0bwMkDQRuAs4E2oEVkhZFxNqSbo8BhYjYLukLwHXAJyUdDPwtUAACWJmN/X0v9s3MrKa6BkhP7SmmTJmy28dwb7jhBu69914ANm7cyFNPPbVHiIwbN46JE4vfID5p0iQ2bNhQu4J6UGmIvC/7c05JWwAf7GbMFKAtItYDSLoLOAd4K0QiYmlJ/+VAZ3SeDTwQES9mYx8ApgJ3VlivmVmvdXfEAMVrIM88s2f7mDGwbFltathvv/3e+nnZsmU8+OCDPPLIIwwfPpwPfOADZT+mO3To0Ld+HjhwYGOdzgKIiNMT5h4NbCxZbgfe203/i4H7uhk7uusASbOB2QBHHnlkQolmZpW79triNZDt2//YNnx4sT3VAQccwMsvv1x23R/+8AcOOugghg8fzhNPPMHy5cvTN5STikJE0giKp5dOzZp+CsyJiD90N6xMW+xl/pkUT12d1puxETEPmAdQKBTKzm1mVisXXVT886qr4Nln4cgjiwHS2Z5i5MiRnHLKKRx77LG87W1v47DDDntr3dSpU7n55ps5/vjjOfrooznppJOq3IPaq/R01m3AauD8bPlTwL8CH+tmTDtwRMlyC7CpaydJZwBXAadFxI6SsR/oMnZZhbWameXmoouqC41y7rjjjrLtQ4cO5b777iu7rvO6xyGHHMLq1avfar/88strW1wPKr1P5E8i4m8jYn32ugZ4Zw9jVgBHSRonaQhwAbCotIOkE4BbgOkRsblk1f3AWZIOknQQcFbWZmZmDaTSEHlV0v/oXJB0CtDtlZuI2AlcSvHNfx2wMCLWSJojaXrWbS6wP3C3pFWSFmVjXwS+RjGIVlA8dfZiL/bLzMz6QKWnsz4P3J5dGwH4PfDpngZFxGJgcZe2r5T8fEY3Y2+jeBrNzMwaVI8hImkAcHRE/HdJbweIiJdyr8zMzBpej6ezIuJNiqeliIiXHCBmZtap0msiD0i6XNIRkg7ufOVamZmZNbxKQ+QzwF8CDwMrs1drXkWZme0rUh8FD3D99dezvfTOxzqo5Cm+A4CZETGuy6unj/iamTWfGj8Lvr+HSCVP8X1T0v8BTu6DeszMGlfns+A737g7nwUPyXcglj4K/swzz+TQQw9l4cKF7Nixg49+9KNcc801vPLKK5x//vm0t7eza9cuvvzlL/P888+zadMmTj/9dA455BCWLl3a88ZyUOlHfJdI+jjwg4jw40XMrDnV4VnwpY+CX7JkCffccw+//OUviQimT5/Oww8/TEdHB+94xzv4yU9+AhSfqTVixAi+9a1vsXTpUg455JDe7mnNVHpN5DJgIbBD0kuSXpbkT2mZ2b4l52fBL1myhCVLlnDCCSdw4okn8sQTT/DUU09x3HHH8eCDD/KlL32Jn/3sZ4wYMaLnyfpIpUciI4CLgHERMUfSkcDh+ZVlZlYHdX4WfERw5ZVX8rnPfW6PdStXrmTx4sVceeWVnHXWWXzlK18pM0Pfq/RI5CbgJGBGtvwycGMuFZmZNaprry0++71Ulc+CL30U/Nlnn81tt93Gtm3bAHjuuefYvHkzmzZtYvjw4cycOZPLL7+cRx99dI+x9VLpkch7I+JESY8BRMTvs4cqmpntO3J4Fnzpo+CnTZvGhRdeyMknFz/HtP/++zN//nza2tr44he/yIABAxg8eDDf+c53AJg9ezbTpk3j8MMPr9uFdVVynVzSf1H8dsMVWZiMApZExAl5F1ipQqEQra2+dcXMemfdunWMHz++3mX0iXL7KmllRBRS56z0dNYNwL3AoZKuBX4O/H3qRs3MrDlU+vW4CyStBP6U4rcOnhsR63KtzMzMGl6l10SIiCeAJ3KsxczM+plKT2eZmTWtfeEe6rz2MdcQkTRV0pOS2iRdUWb9qZIelbRT0nld1l0naY2kdZJukKQ8azWzfdOwYcPYsmVLUwdJRLBlyxaGDRtW87krPp3VW5IGUry/5EygHVghaVFErC3p9iwwC7i8y9j3AacAx2dNPwdOA5blVa+Z7ZtaWlpob2+no6Oj3qXkatiwYbS0tNR83txCBJgCtEXEegBJdwHnAG+FSERsyNa92WVsAMOAIRQv5A8Gns+xVjPbRw0ePJhx48bVu4x+K8/TWaOBjSXL7VlbjyLiEWAp8Lvsdb8/DWZm1njyDJFy1zAqOuko6V3AeKCFYvB8UNKpZfrNltQqqbXZD0XNzBpRniHSDhxRstwCbKpw7EeB5RGxLSK2AfdRfHbXbiJiXkQUIqIwatSoqgs2M7PeyTNEVgBHSRqXPWfrAmBRhWOfBU6TNEjSYIoX1X06y8ysweQWIhGxE7gUuJ9iACyMiDWS5kiaDiBpsqR24BPALZLWZMPvAZ4Gfg38CvhVRPwor1rNzCxNRQ9g7A/8AEYzs97rqwcwmpmZ7cEhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpYs1xCRNFXSk5LaJF1RZv2pkh6VtFPSeV3WHSlpiaR1ktZKGptnrWZm1nu5hYikgcBNwDRgAjBD0oQu3Z4FZgF3lJnidmBuRIwHpgCb86rVzMzSDMpx7ilAW0SsB5B0F3AOsLazQ0RsyNa9WTowC5tBEfFA1m9bjnWamVmiPE9njQY2liy3Z22VeDewVdIPJD0maW52ZLMbSbMltUpq7ejoqEHJZmbWG3mGiMq0RYVjBwHvBy4HJgPvpHjaa/fJIuZFRCEiCqNGjUqt08zMEuUZIu3AESXLLcCmXox9LCLWR8RO4IfAiTWuz8zMqpRniKwAjpI0TtIQ4AJgUS/GHiSp8/Dig5RcSzEzs8aQW4hkRxCXAvcD64CFEbFG0hxJ0wEkTZbUDnwCuEXSmmzsLoqnsh6S9GuKp8ZuzatWMzNLo4hKL1M0tkKhEK2trfUuw8ysX5G0MiIKqeN9x7qZmSVziJiZWTKHiJmZJXOImJlZMoeImZklc4iYmVkyh4iZmSVziJiZWTKHiJmZJXOImJlZMoeImZklc4iYmVkyh4iZmSVziJiZWTKHiJmZJXOImJlZslxDRNJUSU9KapN0RZn1p0p6VNJOSeeVWf92Sc9JujHPOs3MLE1uISJpIHATMA2YAMyQNKFLt2eBWcAde5nma8BP86rRzMyqk+eRyBSgLSLWR8TrwF3AOaUdImJDRDwOvNl1sKRJwGHAkhxrNDOzKuQZIqOBjSXL7VlbjyQNAL4JfDGHuszMrEbyDBGVaYsKx14CLI6Ijd11kjRbUquk1o6Ojl4XaGZm1RmU49ztwBElyy3ApgrHngy8X9IlwP7AEEnbImK3i/MRMQ+YB1AoFCoNKDMzq5E8Q2QFcJSkccBzwAXAhZUMjIiLOn+WNAsodA0QMzOrv9xOZ0XETuBS4H5gHbAwItZImiNpOoCkyZLagU8At0hak1c9ZmZWe4pojrNAhUIhWltb612GmVm/ImllRBRSx/uOdTMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0uWa4hImirpSUltkvb4jnRJp0p6VNJOSeeVtE+U9IikNZIel/TJPOs0M7M0uYWIpIHATcA0YAIwQ9KELt2eBWYBd3Rp3w78eUQcA0wFrpd0YF61mplZmkE5zj0FaIuI9QCS7gLOAdZ2doiIDdm6N0sHRsRvSn7eJGkzMArYmmO9ZmbWS3mezhoNbCxZbs/aekXSFGAI8HSZdbMltUpqXbkSxo6FBQsqm3fBgmL/AQP6Zlyt58hjrjznrOd2Gm3blWj0+qrV7PvXH3T+DmDSpKomiohcXsAngH8pWf4U8M976ftd4Lwy7YcDTwIn9by9SQERw4dHzJ8f3Zo/v9gP/vjKc1yt58hjrjznrOd2Gm3blWj0+qrV7PvXH+z+O5gUUcV7vaL4Zl1zkk4GvhoRZ2fLV2ah9fUyfb8L/Dgi7ilpezuwDPh6RNzd8/YKAa0ADB0KJ520977Ll8OOHXu25zWu1nPkMVeec9ZzO4227Uo0en3Vavb96w92/x0UiGhV6lx5ns5aARwlaZykIcAFwKJKBmb97wVuryRAuir3D7SS9XmNq/UcecyV55z13E6jbbsSjV5ftZp9//qDmv5dV3MY09ML+BDwG4rXM67K2uYA07OfJ1O8VvIKsAVYk7XPBN4AVpW8Jna/rUlvHRqPGdP9odyYMbsfSuc9rtZz5DFXnnPWczuNtu1KNHp91Wr2/esPdv8dVHc6K9cQ6cuXr4lUN1eec9ZzO4227Uo0en3Vavb96w9qeU2k7m/+tXrBpBgzpvJ/iPPnF9NYij4ZV+s58pgrzznruZ1G23YlGr2+ajX7/vUHnb+Dhr2w3tcKhUK0trbWuwwzs35F0sqIKKSO97OzzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLFnTfDpLUgfwTMLQQ4AXEsaNAP6QMK7Wc+QxV55z1nM7jbbtSjR6fdVq9v3rD46OiANSB+f5KPg+FRGjUsZJak35eJukeRExO2WbtZwjj7nynLOe22m0bVei0eurVrPvX38gqap7I3w6K92PGmSOPObKc856bqfRtl2JRq+vWs2+f02vaU5npUo9EjEzawbVvgf6SATm1bsAM7M6quo9cJ8/EjEzs3Q+EjEzs2QOETMzS+YQMTOzZA6REpLGS7pZ0j2SvlDveszM+pqkcyXdKuk/JJ3VU/+mDxFJt0naLGl1l/apkp6U1CbpCoCIWBcRnwfOB/yxXzNrCr18H/xhRHwWmAV8sqe5mz5EgO8CU0sbJA0EbgKmAROAGZImZOumAz8HHurbMs3McvNdevE+mLk6W9+tpg+RiHgYeLFL8xSgLSLWR8TrwF3AOVn/RRHxPuCivq3UzCwfvXkfVNE/APdFxKM9zd00z87qpdHAxpLlduC9kj4AfAwYCiyuQ11mZn2l7Psg8FfAGcAISe+KiJu7m2RfDRGVaYuIWAYs69tSzMzqYm/vgzcAN1Q6SdOfztqLduCIkuUWYFOdajEzq4eavA/uqyGyAjhK0jhJQ4ALgEV1rsnMrC/V5H2w6UNE0p3AI8DRktolXRwRO4FLgfuBdcDCiFhTzzrNzPKS5/ugH8BoZmbJmv5IxMzM8uMQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOEbMqSDpQ0iX1rsOsXhwiZtU5ENgjRLLHbJs1vX31AYxmtfIN4E8krQLeALYBvwMmAhMkzQT+FzAE+C/gkojYlX1j3DUUnxj9NPAXEbGtHjtgVg3fsW5WBUljgR9HxLHZVwn8BDg2In4raTxwHfCxiHhD0reB5RS/ZuAHwLSIeEXSl4ChETGnLjthVgUfiZjV1i8j4rfZz38KTAJWSAJ4G7AZOIniN8n9ImsfQvG5Rmb9jkPErLZeKflZwL9FxJWlHSR9BHggImb0aWVmOfCFdbPqvAwcsJd1DwHnSToUQNLBksZQPKV1iqR3Ze3DJb27T6o1qzEfiZhVISK2SPqFpNXAq8DzJevWSroaWCJpAMUL738ZEcslzQLulDQ063418Js+Lt+sar6wbmZmyXw6y8zMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMws2f8HPMNm0txu95EAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.semilogx(tree_arr_ert, train_err, 'b-o', label = 'train')\n",
    "plt.semilogx(tree_arr_ert, test_err, 'r-o', label = 'test')\n",
    "plt.xlim([np.max(tree_arr_ert), np.min(tree_arr_ert)])\n",
    "plt.title('Error vs. tree')\n",
    "plt.xlabel('tree')\n",
    "plt.ylabel('error')\n",
    "plt.legend()\n",
    "pass"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 3) adaboost"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.8007566204287516 0.8241176470588235\n"
     ]
    }
   ],
   "source": [
    "ada = ensemble.AdaBoostClassifier(n_estimators = 100)\n",
    "ada.fit(X_train, y_train)\n",
    "\n",
    "err_train = np.mean(y_train != ada.predict(X_train))\n",
    "err_test = np.mean(y_test != ada.predict(X_test))\n",
    "\n",
    "print(err_train, err_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 100  200  300  400  500  600  700  800  900 1000 1100]\n"
     ]
    }
   ],
   "source": [
    "tree_arr_ada = np.linspace(100, 1000, 10)\n",
    "tree_arr_ada = tree_arr_ada.astype('int')\n",
    "print(tree_arr_ada)\n",
    "test_err = []\n",
    "train_err = []\n",
    "for tree in tree_arr_ada:\n",
    "    ert = ensemble.AdaBoostClassifier(n_estimators = tree)\n",
    "    ert.fit(X_train, y_train)\n",
    "\n",
    "    y_train_pred = rf.predict(X_train)\n",
    "    y_test_pred = rf.predict(X_test)\n",
    "    train_err.append(np.mean(y_train != y_train_pred))\n",
    "    test_err.append(np.mean(y_test != y_test_pred))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEaCAYAAADQVmpMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAHrhJREFUeJzt3X+8VXWd7/HXm9+hhoroNY4CTWbgj4uyIc2bZuMPqAmtzERpYvJG5Tj3MePDHunVmqRxauTWOI6W4lyncUANLRsqfIg6kNVDJg5Kxg/NI6Ec6coRI0UUBT/3j72ObQ6bc/b57r3O3mfzfj4e+8FZ3/X9ftdncXC/XWvttbYiAjMzsxQD6l2AmZn1Xw4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RaxqSNkh6VdK2kteN9a6rWpLGSgpJg+pdi1lX/kdpzeYjEfFgT50kDYqInT219XaOemmkWmzf4iMR2ydImiXpF5L+UdKLwFf30jZA0tWSnpG0WdLtkkZkc3QeEVws6VngP8tsZ52kPytZHiTpBUknShomab6kLZK2Sloh6bAKyn84+3NrdnR1crnas+19Jqvh95LulzSmpJb3SHpA0ouSnpR0fvJfqFnGIWL7kvcC64FDgWv30jYre50OvBPYH+h6Suw0YDxwdplt3AnMKFk+G3ghIh4FPg2MAI4ARgKfB16toO5Tsz8PjIj9I+KRcrVLOhf438DHgFHAz7J6kLQf8ABwR9Z/BvBtScdUsH2zvXKIWLP5YfZ/+Z2vz5as2xQR/xwROyPi1b20XQR8KyLWR8Q24Erggi7XI74aEa+UzFHqDmC6pOHZ8oVZG8AbFMPjXRGxKyJWRsRLVexr19o/B3w9ItZlp7b+HpiYHY38GbAhIv416/8o8H3gvCq2b+ZrItZ0zu3mmsjGCtreATxTsvwMxf9OSk87lZsHgIhok7QO+IikHwHTgROy1f9O8SjkLkkHAvOBqyLijb3N14OudYwB/knSN0vaBIzO1r1X0taSdYOymsySOURsX1LukdVd2zZRfMPtdCSwE3geaOlmnlKdp7QGAGsjog0gC4trgGskjQUWA08C/zeh7nLtG4FrI2JB147Z0chPI+LMHrZl1is+nWW2uzuBv5E0TtL+FE8Jfa+Xn3y6CzgL+AJ/PJWFpNMlHSdpIPASxdNbuyqYrwN4k+I1mu7cDFzZeZ1D0ghJn8jW/Rh4t6RPSRqcvSZLGt+L/TLbg0PEms2Putwncm8vx99G8RTPw8BvgdeAv+rNBBHxO+AR4H3A90pW/TfgHooBsg74KcVTWki6WdLNe5lvO8WL/r/IrvOctJd+9wL/QPF02UvAamBatu5lisF2AcWjrf+X9R3am30z60r+UiozM0vlIxEzM0vmEDEzs2S5hoikqdmdsW2Sriiz/jJJayU9LumhLnfXXidpTXb37Q2SlGetZmbWe7mFSPYJlJsoXtibAMyQNKFLt8eAQkQcT/GC43XZ2PcBpwDHA8cCkyneJWxmZg0kzyORKUBbdufv6xQ/9nhOaYeIWJp98gRgObt/Dn8YMITip0cGU/ycvpmZNZA8bzYcze531LZTfNbP3lwM3AcQEY9IWgr8juIdtzdGxLquAyTNBmYD7LfffpPe85731Kh0M7N9w8qVK1+IiFGp4/MMkXLXMMp+nljSTKBAdspK0rsoPuCu88jkAUmnRsTDpeMiYh4wD6BQKERra2uNSjcz2zdIeqbnXnuX5+msdorPCerUQvEmp91IOgO4CpgeETuy5o8CyyNiW/YQvPuAsjdYmZlZ/eQZIiuAo7LHRwyheKfsotIOkk4AbqEYIJtLVj0LnJZ9F8Ngikcoe5zOMjOz+sotRLJnDV0K3E8xABZGxBpJcyRNz7rNpfh9DXdLWiWpM2TuAZ4Gfg38CvhVRPwor1rNzCxN0zz2xNdEzMx6T9LKiCikjm+eO9ZXroSxY2HBHk/BLm/BgmL/AQP6Zlyt58hjrjznrOd2Gm3blWj0+qrV7PvXH2S/g0kwqap5IqIpXpMgAiKGD4+YPz+6NX9+sV/nmLzH1XqOPObKc856bqfRtl2JRq+vWs2+f/1Bye9gEkRU8d7bPKezpHjrZNbQoXBSNx/mWr4cduzYsz2vcbWeI4+58pyznttptG1XotHrq1az719/UPI7KACtEcmPlWqe01mlyv0DrWR9XuNqPUcec+U5Zz2302jbrkSj11etZt+//qCWf9fVHMY00mtS6aHxmDHdH8qNGbP7oXTe42o9Rx5z5TlnPbfTaNuuRKPXV61m37/+oOR3UO3prOY7Ehk+HK69tvs+115b7NdX42o9Rx5z5TlnPbfTaNuuRKPXV61m37/+oNzvIFU1CdRIr0md/ydT6cW5+fOL/aW+GVfrOfKYK88567mdRtt2JRq9vmo1+/71B9nvwBfWM75PxMys93yfiJmZ1Y1DxMzMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMwsWa4hImmqpCcltUm6osz6yyStlfS4pIckjSlZd6SkJZLWZX3G5lmrmZn1Xm4hImkgcBMwDZgAzJA0oUu3x4BCRBxP8XvVrytZdzswNyLGA1OAzXnVamZmafI8EpkCtEXE+oh4HbgLOKe0Q0QsjYjt2eJyoAUgC5tBEfFA1m9bST8zM2sQeYbIaGBjyXJ71rY3FwP3ZT+/G9gq6QeSHpM0Nzuy2Y2k2ZJaJbV2dHTUrHAzM6tMniFS7usWyz4yWNJMit/SODdrGgS8H7gcmAy8E5i1x2QR8yKiEBGFUaNG1aJmMzPrhTxDpB04omS5BdjUtZOkM4CrgOkRsaNk7GPZqbCdwA+BE3Os1czMEuQZIiuAoySNkzQEuABYVNpB0gnALRQDZHOXsQdJ6jy8+CCwNsdazcwsQW4hkh1BXArcD6wDFkbEGklzJE3Pus0F9gfulrRK0qJs7C6Kp7IekvRriqfGbs2rVjMzS+NvNjQz24f5mw3NzKxuHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlizXEJE0VdKTktokXVFm/WWS1kp6XNJDksZ0Wf92Sc9JujHPOs3MLE1uISJpIHATMA2YAMyQNKFLt8eAQkQcD9wDXNdl/deAn+ZVo5mZVSfPI5EpQFtErI+I14G7gHNKO0TE0ojYni0uB1o610maBBwGLMmxRjMzq0KeITIa2Fiy3J617c3FwH0AkgYA3wS+2N0GJM2W1CqptaOjo8pyzcyst/IMEZVpi7IdpZlAAZibNV0CLI6IjeX6vzVZxLyIKEREYdSoUVUVa2ZmvTcox7nbgSNKlluATV07SToDuAo4LSJ2ZM0nA++XdAmwPzBE0raI2OPivJmZ1U+eIbICOErSOOA54ALgwtIOkk4AbgGmRsTmzvaIuKikzyyKF98dIGZmDSa301kRsRO4FLgfWAcsjIg1kuZImp51m0vxSONuSaskLcqrHjMzqz1FlL1M0e8UCoVobW2tdxlmZv2KpJURUUgd7zvWzcwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLFmuISJpqqQnJbVJ2uM70iVdJmmtpMclPSRpTNY+UdIjktZk6z6ZZ51mZpYmtxCRNBC4CZgGTABmSJrQpdtjQCEijgfuAa7L2rcDfx4RxwBTgeslHZhXrWZmlibPI5EpQFtErI+I14G7gHNKO0TE0ojYni0uB1qy9t9ExFPZz5uAzcCoHGs1M7MEeYbIaGBjyXJ71rY3FwP3dW2UNAUYAjxdZt1sSa2SWjs6Oqos18zMeivPEFGZtijbUZoJFIC5XdoPB/4d+IuIeHOPySLmRUQhIgqjRvlAxcysrw3Kce524IiS5RZgU9dOks4ArgJOi4gdJe1vB34CXB0Ry3Os08zMEuV5JLICOErSOElDgAuARaUdJJ0A3AJMj4jNJe1DgHuB2yPi7hxrNDOzKuQWIhGxE7gUuB9YByyMiDWS5kiannWbC+wP3C1plaTOkDkfOBWYlbWvkjQxr1rNzCyNIspepuh3CoVCtLa21rsMM7N+RdLKiCikjvcd62ZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWLM8HMJqZNbw33niD9vZ2XnvttXqXkqthw4bR0tLC4MGDazqvQ8TM9mnt7e0ccMABjB07FqncN1j0fxHBli1baG9vZ9y4cTWd26ezzGyf9tprrzFy5MimDRAASYwcOTKXoy2HiJnt85o5QDrltY8OETOzOtq6dSvf/va3ez3uQx/6EFu3bs2hot5xiJiZ9cKCBTB2LAwYUPxzwYLq5ttbiOzatavbcYsXL+bAAw+sbuM10GOISBoo6W/6ohgzs0a2YAHMng3PPAMRxT9nz64uSK644gqefvppJk6cyOTJkzn99NO58MILOe644wA499xzmTRpEscccwzz5s17a9zYsWN54YUX2LBhA+PHj+ezn/0sxxxzDGeddRavvvpqtbtasR4/nRURuySdA/xjbyeXNBX4J2Ag8C8R8Y0u6y8D/iewE+gAPhMRz2TrPg1cnXX9u4j4t95u38ysN/76r2HVqr2vX74cduzYvW37drj4Yrj11vJjJk6E66/f+5zf+MY3WL16NatWrWLZsmV8+MMfZvXq1W99iuq2227j4IMP5tVXX2Xy5Ml8/OMfZ+TIkbvN8dRTT3HnnXdy6623cv755/P973+fmTNnVrLLVav0I76/kHQj8D3glc7GiHh0bwMkDQRuAs4E2oEVkhZFxNqSbo8BhYjYLukLwHXAJyUdDPwtUAACWJmN/X0v9s3MrKa6BkhP7SmmTJmy28dwb7jhBu69914ANm7cyFNPPbVHiIwbN46JE4vfID5p0iQ2bNhQu4J6UGmIvC/7c05JWwAf7GbMFKAtItYDSLoLOAd4K0QiYmlJ/+VAZ3SeDTwQES9mYx8ApgJ3VlivmVmvdXfEAMVrIM88s2f7mDGwbFltathvv/3e+nnZsmU8+OCDPPLIIwwfPpwPfOADZT+mO3To0Ld+HjhwYGOdzgKIiNMT5h4NbCxZbgfe203/i4H7uhk7uusASbOB2QBHHnlkQolmZpW79triNZDt2//YNnx4sT3VAQccwMsvv1x23R/+8AcOOugghg8fzhNPPMHy5cvTN5STikJE0giKp5dOzZp+CsyJiD90N6xMW+xl/pkUT12d1puxETEPmAdQKBTKzm1mVisXXVT886qr4Nln4cgjiwHS2Z5i5MiRnHLKKRx77LG87W1v47DDDntr3dSpU7n55ps5/vjjOfrooznppJOq3IPaq/R01m3AauD8bPlTwL8CH+tmTDtwRMlyC7CpaydJZwBXAadFxI6SsR/oMnZZhbWameXmoouqC41y7rjjjrLtQ4cO5b777iu7rvO6xyGHHMLq1avfar/88strW1wPKr1P5E8i4m8jYn32ugZ4Zw9jVgBHSRonaQhwAbCotIOkE4BbgOkRsblk1f3AWZIOknQQcFbWZmZmDaTSEHlV0v/oXJB0CtDtlZuI2AlcSvHNfx2wMCLWSJojaXrWbS6wP3C3pFWSFmVjXwS+RjGIVlA8dfZiL/bLzMz6QKWnsz4P3J5dGwH4PfDpngZFxGJgcZe2r5T8fEY3Y2+jeBrNzMwaVI8hImkAcHRE/HdJbweIiJdyr8zMzBpej6ezIuJNiqeliIiXHCBmZtap0msiD0i6XNIRkg7ufOVamZmZNbxKQ+QzwF8CDwMrs1drXkWZme0rUh8FD3D99dezvfTOxzqo5Cm+A4CZETGuy6unj/iamTWfGj8Lvr+HSCVP8X1T0v8BTu6DeszMGlfns+A737g7nwUPyXcglj4K/swzz+TQQw9l4cKF7Nixg49+9KNcc801vPLKK5x//vm0t7eza9cuvvzlL/P888+zadMmTj/9dA455BCWLl3a88ZyUOlHfJdI+jjwg4jw40XMrDnV4VnwpY+CX7JkCffccw+//OUviQimT5/Oww8/TEdHB+94xzv4yU9+AhSfqTVixAi+9a1vsXTpUg455JDe7mnNVHpN5DJgIbBD0kuSXpbkT2mZ2b4l52fBL1myhCVLlnDCCSdw4okn8sQTT/DUU09x3HHH8eCDD/KlL32Jn/3sZ4wYMaLnyfpIpUciI4CLgHERMUfSkcDh+ZVlZlYHdX4WfERw5ZVX8rnPfW6PdStXrmTx4sVceeWVnHXWWXzlK18pM0Pfq/RI5CbgJGBGtvwycGMuFZmZNaprry0++71Ulc+CL30U/Nlnn81tt93Gtm3bAHjuuefYvHkzmzZtYvjw4cycOZPLL7+cRx99dI+x9VLpkch7I+JESY8BRMTvs4cqmpntO3J4Fnzpo+CnTZvGhRdeyMknFz/HtP/++zN//nza2tr44he/yIABAxg8eDDf+c53AJg9ezbTpk3j8MMPr9uFdVVynVzSf1H8dsMVWZiMApZExAl5F1ipQqEQra2+dcXMemfdunWMHz++3mX0iXL7KmllRBRS56z0dNYNwL3AoZKuBX4O/H3qRs3MrDlU+vW4CyStBP6U4rcOnhsR63KtzMzMGl6l10SIiCeAJ3KsxczM+plKT2eZmTWtfeEe6rz2MdcQkTRV0pOS2iRdUWb9qZIelbRT0nld1l0naY2kdZJukKQ8azWzfdOwYcPYsmVLUwdJRLBlyxaGDRtW87krPp3VW5IGUry/5EygHVghaVFErC3p9iwwC7i8y9j3AacAx2dNPwdOA5blVa+Z7ZtaWlpob2+no6Oj3qXkatiwYbS0tNR83txCBJgCtEXEegBJdwHnAG+FSERsyNa92WVsAMOAIRQv5A8Gns+xVjPbRw0ePJhx48bVu4x+K8/TWaOBjSXL7VlbjyLiEWAp8Lvsdb8/DWZm1njyDJFy1zAqOuko6V3AeKCFYvB8UNKpZfrNltQqqbXZD0XNzBpRniHSDhxRstwCbKpw7EeB5RGxLSK2AfdRfHbXbiJiXkQUIqIwatSoqgs2M7PeyTNEVgBHSRqXPWfrAmBRhWOfBU6TNEjSYIoX1X06y8ysweQWIhGxE7gUuJ9iACyMiDWS5kiaDiBpsqR24BPALZLWZMPvAZ4Gfg38CvhVRPwor1rNzCxNRQ9g7A/8AEYzs97rqwcwmpmZ7cEhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpYs1xCRNFXSk5LaJF1RZv2pkh6VtFPSeV3WHSlpiaR1ktZKGptnrWZm1nu5hYikgcBNwDRgAjBD0oQu3Z4FZgF3lJnidmBuRIwHpgCb86rVzMzSDMpx7ilAW0SsB5B0F3AOsLazQ0RsyNa9WTowC5tBEfFA1m9bjnWamVmiPE9njQY2liy3Z22VeDewVdIPJD0maW52ZLMbSbMltUpq7ejoqEHJZmbWG3mGiMq0RYVjBwHvBy4HJgPvpHjaa/fJIuZFRCEiCqNGjUqt08zMEuUZIu3AESXLLcCmXox9LCLWR8RO4IfAiTWuz8zMqpRniKwAjpI0TtIQ4AJgUS/GHiSp8/Dig5RcSzEzs8aQW4hkRxCXAvcD64CFEbFG0hxJ0wEkTZbUDnwCuEXSmmzsLoqnsh6S9GuKp8ZuzatWMzNLo4hKL1M0tkKhEK2trfUuw8ysX5G0MiIKqeN9x7qZmSVziJiZWTKHiJmZJXOImJlZMoeImZklc4iYmVkyh4iZmSVziJiZWTKHiJmZJXOImJlZMoeImZklc4iYmVkyh4iZmSVziJiZWTKHiJmZJXOImJlZslxDRNJUSU9KapN0RZn1p0p6VNJOSeeVWf92Sc9JujHPOs3MLE1uISJpIHATMA2YAMyQNKFLt2eBWcAde5nma8BP86rRzMyqk+eRyBSgLSLWR8TrwF3AOaUdImJDRDwOvNl1sKRJwGHAkhxrNDOzKuQZIqOBjSXL7VlbjyQNAL4JfDGHuszMrEbyDBGVaYsKx14CLI6Ijd11kjRbUquk1o6Ojl4XaGZm1RmU49ztwBElyy3ApgrHngy8X9IlwP7AEEnbImK3i/MRMQ+YB1AoFCoNKDMzq5E8Q2QFcJSkccBzwAXAhZUMjIiLOn+WNAsodA0QMzOrv9xOZ0XETuBS4H5gHbAwItZImiNpOoCkyZLagU8At0hak1c9ZmZWe4pojrNAhUIhWltb612GmVm/ImllRBRSx/uOdTMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0uWa4hImirpSUltkvb4jnRJp0p6VNJOSeeVtE+U9IikNZIel/TJPOs0M7M0uYWIpIHATcA0YAIwQ9KELt2eBWYBd3Rp3w78eUQcA0wFrpd0YF61mplZmkE5zj0FaIuI9QCS7gLOAdZ2doiIDdm6N0sHRsRvSn7eJGkzMArYmmO9ZmbWS3mezhoNbCxZbs/aekXSFGAI8HSZdbMltUpqXbkSxo6FBQsqm3fBgmL/AQP6Zlyt58hjrjznrOd2Gm3blWj0+qrV7PvXH3T+DmDSpKomiohcXsAngH8pWf4U8M976ftd4Lwy7YcDTwIn9by9SQERw4dHzJ8f3Zo/v9gP/vjKc1yt58hjrjznrOd2Gm3blWj0+qrV7PvXH+z+O5gUUcV7vaL4Zl1zkk4GvhoRZ2fLV2ah9fUyfb8L/Dgi7ilpezuwDPh6RNzd8/YKAa0ADB0KJ520977Ll8OOHXu25zWu1nPkMVeec9ZzO4227Uo0en3Vavb96w92/x0UiGhV6lx5ns5aARwlaZykIcAFwKJKBmb97wVuryRAuir3D7SS9XmNq/UcecyV55z13E6jbbsSjV5ftZp9//qDmv5dV3MY09ML+BDwG4rXM67K2uYA07OfJ1O8VvIKsAVYk7XPBN4AVpW8Jna/rUlvHRqPGdP9odyYMbsfSuc9rtZz5DFXnnPWczuNtu1KNHp91Wr2/esPdv8dVHc6K9cQ6cuXr4lUN1eec9ZzO4227Uo0en3Vavb96w9qeU2k7m/+tXrBpBgzpvJ/iPPnF9NYij4ZV+s58pgrzznruZ1G23YlGr2+ajX7/vUHnb+Dhr2w3tcKhUK0trbWuwwzs35F0sqIKKSO97OzzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLFnTfDpLUgfwTMLQQ4AXEsaNAP6QMK7Wc+QxV55z1nM7jbbtSjR6fdVq9v3rD46OiANSB+f5KPg+FRGjUsZJak35eJukeRExO2WbtZwjj7nynLOe22m0bVei0eurVrPvX38gqap7I3w6K92PGmSOPObKc856bqfRtl2JRq+vWs2+f02vaU5npUo9EjEzawbVvgf6SATm1bsAM7M6quo9cJ8/EjEzs3Q+EjEzs2QOETMzS+YQMTOzZA6REpLGS7pZ0j2SvlDveszM+pqkcyXdKuk/JJ3VU/+mDxFJt0naLGl1l/apkp6U1CbpCoCIWBcRnwfOB/yxXzNrCr18H/xhRHwWmAV8sqe5mz5EgO8CU0sbJA0EbgKmAROAGZImZOumAz8HHurbMs3McvNdevE+mLk6W9+tpg+RiHgYeLFL8xSgLSLWR8TrwF3AOVn/RRHxPuCivq3UzCwfvXkfVNE/APdFxKM9zd00z87qpdHAxpLlduC9kj4AfAwYCiyuQ11mZn2l7Psg8FfAGcAISe+KiJu7m2RfDRGVaYuIWAYs69tSzMzqYm/vgzcAN1Q6SdOfztqLduCIkuUWYFOdajEzq4eavA/uqyGyAjhK0jhJQ4ALgEV1rsnMrC/V5H2w6UNE0p3AI8DRktolXRwRO4FLgfuBdcDCiFhTzzrNzPKS5/ugH8BoZmbJmv5IxMzM8uMQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOEbMqSDpQ0iX1rsOsXhwiZtU5ENgjRLLHbJs1vX31AYxmtfIN4E8krQLeALYBvwMmAhMkzQT+FzAE+C/gkojYlX1j3DUUnxj9NPAXEbGtHjtgVg3fsW5WBUljgR9HxLHZVwn8BDg2In4raTxwHfCxiHhD0reB5RS/ZuAHwLSIeEXSl4ChETGnLjthVgUfiZjV1i8j4rfZz38KTAJWSAJ4G7AZOIniN8n9ImsfQvG5Rmb9jkPErLZeKflZwL9FxJWlHSR9BHggImb0aWVmOfCFdbPqvAwcsJd1DwHnSToUQNLBksZQPKV1iqR3Ze3DJb27T6o1qzEfiZhVISK2SPqFpNXAq8DzJevWSroaWCJpAMUL738ZEcslzQLulDQ063418Js+Lt+sar6wbmZmyXw6y8zMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMws2f8HPMNm0txu95EAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.semilogx(tree_arr_ada, train_err, 'b-o', label = 'train')\n",
    "plt.semilogx(tree_arr_ada, test_err, 'r-o', label = 'test')\n",
    "plt.xlim([np.max(tree_arr_ada), np.min(tree_arr_ada)])\n",
    "plt.title('Error vs. tree')\n",
    "plt.xlabel('tree')\n",
    "plt.ylabel('error')\n",
    "plt.legend()\n",
    "pass"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 4) gradient boosted trees\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.296343001261034 0.3947058823529412\n"
     ]
    }
   ],
   "source": [
    "gbt = ensemble.GradientBoostingClassifier(n_estimators = 100)\n",
    "gbt.fit(X_train, y_train)\n",
    "\n",
    "err_train = np.mean(y_train != gbt.predict(X_train))\n",
    "err_test = np.mean(y_test != gbt.predict(X_test))\n",
    "\n",
    "print(err_train, err_test)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[ 100  200  300  400  500  600  700  800  900 1000 1100]\n"
     ]
    }
   ],
   "source": [
    "tree_arr_gbt = np.linspace(100, 1000, 10)\n",
    "tree_arr_gbt = tree_arr_gbt.astype('int')\n",
    "print(tree_arr_gbt)\n",
    "test_err = []\n",
    "train_err = []\n",
    "for tree in tree_arr_gbt:\n",
    "    ert = ensemble.GradientBoostingClassifier(n_estimators = tree)\n",
    "    ert.fit(X_train, y_train)\n",
    "\n",
    "    y_train_pred = rf.predict(X_train)\n",
    "    y_test_pred = rf.predict(X_test)\n",
    "    train_err.append(np.mean(y_train != y_train_pred))\n",
    "    test_err.append(np.mean(y_test != y_test_pred))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEaCAYAAADQVmpMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAHrBJREFUeJzt3X+8VXWd7/HXm9+hhoroNY78aDIDf1yUDWneNBt/QE1oZSZKE5M3Kse5jxkf9kiv1ow0To3cGsfJUpzrNA6ooWWXCh+gXsjqIRMHJOOH5pFQjnQFMVJEUfBz/9jr2OawD2ef797r7H027+fjsR+c9V3f73d9Fgf327XWXmsrIjAzM0vRr94FmJlZ3+UQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQsaYhaaOkVyXtKHl9q951VUvSGEkhaUC9azHrzP8ordl8JCIe6q6TpAERsbu7tp7OUS+NVIsdWHwkYgcESTMl/ULSP0l6Efi7Ltr6SbpO0jOStki6U9KwbI6OI4LLJD0L/N8y21kv6c9KlgdIekHSKZKGSJonaZuk7ZJWSDqqgvIfyf7cnh1dnVau9mx7n8lq+L2kxZJGl9TyHkkPSnpR0pOSLkr+CzXLOETsQPJeYANwJHBDF20zs9dZwDuBg4HOp8TOBMYB55XZxt3A9JLl84AXImIV8GlgGHAMMBz4PPBqBXWfkf15aEQcHBGPlqtd0gXA/wQ+BowAfpbVg6SDgAeBu7L+04FvSzq+gu2bdckhYs3mh9n/5Xe8PluybnNE/EtE7I6IV7touxT4ZkRsiIgdwDXAxZ2uR/xdRLxSMkepu4BpkoZmy5dkbQBvUAyPd0XEnohYGREvVbGvnWv/HPC1iFifndr6B2BCdjTyZ8DGiPi3rP8q4PvAhVVs38zXRKzpXLCfayKbKmh7B/BMyfIzFP87KT3tVG4eACKiTdJ64COSfgRMA07OVv8HxaOQeyQdCswDro2IN7qarxud6xgN/LOkb5S0CRiZrXuvpO0l6wZkNZklc4jYgaTcI6s7t22m+IbbYRSwG3geaNnPPKU6Tmn1A9ZFRBtAFhbXA9dLGgMsAp4E/ndC3eXaNwE3RMT8zh2zo5GfRsQ53WzLrEd8Ostsb3cDfyNprKSDKZ4S+l4PP/l0D3Au8AX+eCoLSWdJOlFSf+Aliqe39lQw31bgTYrXaPbnVuCajusckoZJ+kS27sfAuyV9StLA7DVJ0rge7JfZPhwi1mx+1Ok+kft7OP4Oiqd4HgF+C7wG/FVPJoiI3wGPAu8Dvley6r8A91EMkPXATyme0kLSrZJu7WK+nRQv+v8iu85zahf97gf+keLpspeANcDUbN3LFIPtYopHW/8v6zu4J/tm1pn8pVRmZpbKRyJmZpYs1xCRNCW7qalN0tVl1l8paZ2kxyU93OnGqBslrc1unLpZkvKs1czMei63EMkuHt5C8ZzseGC6pPGduj0GFCLiJIrnim/Mxr4POB04CTgBmETxBi8zM2sgeR6JTAbaspu2Xqf4iZXzSztExNLsoiHAcvb+COUQYBDFC38DKX7E0szMGkieITKSvW+Gas/aunIZ8ABA9liHpcDvstfiiFifU51mZpYoz5sNy13DKPtRMEkzgALZKStJ76L4bKKOI5MHJZ0REY90GjcLmAVw0EEHTXzPe95To9LNzA4MK1eufCEiRqSOzzNE2ik+4qFDC8XPp+9F0tnAtcCZEbEra/4osDx7dhGSHgBO5Y9PMwUgIuYCcwEKhUK0trbWeh/MzJqapGe679W1PE9nrQCOze78HUTxJqeFpR0knQzcBkyLiC0lq54Fzsweoz2Q4hGKT2eZmTWY3EIke0zEFcBiigGwICLWSpotaVrWbQ7FR23fK2m1pI6QuQ94Gvg18CvgVxHxo7xqNTOzNE1zx7pPZ5mZ9ZyklRFRSB3vO9bNzCxZ84TIypUwZgzM3+cp2OXNn1/s369f74yr9Rx5zJXnnPXcTqNtuxKNXl+1mn3/+oLsdzARJlY1T0Q0xWsiREDE0KER8+bFfs2bV+zXMSbvcbWeI4+58pyznttptG1XotHrq1az719fUPI7mAgRVbz3Ns81ESneuiIyeDCcWvZp2UXLl8OuXfu25zWu1nPkMVeec9ZzO4227Uo0en3Vavb96wtKfgcFoDUi+dmEzXM6q1S5f6CVrM9rXK3nyGOuPOes53YabduVaPT6qtXs+9cX1PLvuprDmEZ6TSw9NB49ev+HcqNH730onfe4Ws+Rx1x5zlnP7TTativR6PVVq9n3ry8o+R1Uezqr+Y5Ehg6FG27Yf58bbij2661xtZ4jj7nynLOe22m0bVei0eurVrPvX19Q7neQqpoEaqTXxI7/k6n04ty8ecX+Uu+Mq/UcecyV55z13E6jbbsSjV5ftZp9//qC7HfgC+sZ32xoZtZzvtnQzMzqxiFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJcg0RSVMkPSmpTdLVZdZfKWmdpMclPSxpdMm6UZKWSFqf9RmTZ61mZtZzuYWIpP7ALcBUYDwwXdL4Tt0eAwoRcRLF71W/sWTdncCciBgHTAa25FWrmZmlyfNIZDLQFhEbIuJ14B7g/NIOEbE0InZmi8uBFoAsbAZExINZvx0l/czMrEHkGSIjgU0ly+1ZW1cuAx7Ifn43sF3SDyQ9JmlOdmSzF0mzJLVKat26dWvNCjczs8rkGSLlvimr7NMeJc2g+AVbc7KmAcD7gauAScA7gZn7TBYxNyIKEVEYMWJELWo2M7MeyDNE2oFjSpZbgM2dO0k6G7gWmBYRu0rGPpadCtsN/BA4JcdazcwsQZ4hsgI4VtJYSYOAi4GFpR0knQzcRjFAtnQae5ikjsOLDwLrcqzVzMwS5BYi2RHEFcBiYD2wICLWSpotaVrWbQ5wMHCvpNWSFmZj91A8lfWwpF9TPDV2e161mplZGn8plZnZAcxfSmVmZnXjEDEzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZLmGiKQpkp6U1Cbp6jLrr5S0TtLjkh6WNLrT+rdLek7St/Ks08zM0uQWIpL6A7cAU4HxwHRJ4zt1ewwoRMRJwH3AjZ3WfxX4aV41mplZdfI8EpkMtEXEhoh4HbgHOL+0Q0QsjYid2eJyoKVjnaSJwFHAkhxrNDOzKuQZIiOBTSXL7VlbVy4DHgCQ1A/4BvDF3KozM7OqDchxbpVpi7IdpRlAATgza7ocWBQRm6Ry07w1bhYwC2DUqFFVFWtmZj2XZ4i0A8eULLcAmzt3knQ2cC1wZkTsyppPA94v6XLgYGCQpB0RsdfF+YiYC8wFKBQKZQPKzMzyk2eIrACOlTQWeA64GLiktIOkk4HbgCkRsaWjPSIuLekzk+LF930+3WVmZvWV2zWRiNgNXAEsBtYDCyJiraTZkqZl3eZQPNK4V9JqSQvzqsfMzGpPEc1xFqhQKERra2u9yzAz61MkrYyIQup437FuZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyXINEUlTJD0pqU3SPt+RLulKSeskPS7pYUmjs/YJkh6VtDZb98k86zQzszS5hYik/sAtwFRgPDBd0vhO3R4DChFxEnAfcGPWvhP484g4HpgC3CTp0LxqNTOzNHkeiUwG2iJiQ0S8DtwDnF/aISKWRsTObHE50JK1/yYinsp+3gxsAUbkWKuZmSXIM0RGAptKltuztq5cBjzQuVHSZGAQ8HSZdbMktUpq3bp1a5XlmplZT+UZIirTFmU7SjOAAjCnU/vRwH8AfxERb+4zWcTciChERGHECB+omJn1tgE5zt0OHFOy3AJs7txJ0tnAtcCZEbGrpP3twE+A6yJieY51mplZojyPRFYAx0oaK2kQcDGwsLSDpJOB24BpEbGlpH0QcD9wZ0Tcm2ONZmZWhdxCJCJ2A1cAi4H1wIKIWCtptqRpWbc5wMHAvZJWS+oImYuAM4CZWftqSRPyqtXMzNIoouxlij6nUChEa2trvcswM+tTJK2MiELqeN+xbmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpYsz6f4mpk1vDfeeIP29nZee+21epeSqyFDhtDS0sLAgQNrOq9DxMwOaO3t7RxyyCGMGTMGqdzXIPV9EcG2bdtob29n7NixNZ3bp7PM7ID22muvMXz48KYNEABJDB8+PJejLYeImR3wmjlAOuS1jw4RM7M62r59O9/+9rd7PO5DH/oQ27dvz6GinnGImJn1wPz5MGYM9OtX/HP+/Orm6ypE9uzZs99xixYt4tBDD61u4zXQbYhI6i/pb3qjGDOzRjZ/PsyaBc88AxHFP2fNqi5Irr76ap5++mkmTJjApEmTOOuss7jkkks48cQTAbjggguYOHEixx9/PHPnzn1r3JgxY3jhhRfYuHEj48aN47Of/SzHH3885557Lq+++mq1u1qxbj+dFRF7JJ0P/FNPJ5c0BfhnoD/wrxHx9U7rrwT+O7Ab2Ap8JiKeydZ9Grgu6/r3EfHvPd2+mVlP/PVfw+rVXa9fvhx27dq7bedOuOwyuP328mMmTICbbup6zq9//eusWbOG1atXs2zZMj784Q+zZs2atz5Fdccdd3D44Yfz6quvMmnSJD7+8Y8zfPjwveZ46qmnuPvuu7n99tu56KKL+P73v8+MGTMq2eWqVfoR319I+hbwPeCVjsaIWNXVAEn9gVuAc4B2YIWkhRGxrqTbY0AhInZK+gJwI/BJSYcDfwsUgABWZmN/34N9MzOrqc4B0l17ismTJ+/1Mdybb76Z+++/H4BNmzbx1FNP7RMiY8eOZcKECQBMnDiRjRs31q6gblQaIu/L/pxd0hbAB/czZjLQFhEbACTdA5wPvBUiEbG0pP9yoCM6zwMejIgXs7EPAlOAuyus18ysx/Z3xADFayDPPLNv++jRsGxZbWo46KCD3vp52bJlPPTQQzz66KMMHTqUD3zgA2U/pjt48OC3fu7fv39jnc4CiIizEuYeCWwqWW4H3ruf/pcBD+xn7MjOAyTNAmYBjBo1KqFEM7PK3XBD8RrIzp1/bBs6tNie6pBDDuHll18uu+4Pf/gDhx12GEOHDuWJJ55g+fLl6RvKSUUhImkYxdNLZ2RNPwVmR8Qf9jesTFt0Mf8MiqeuzuzJ2IiYC8wFKBQKZec2M6uVSy8t/nnttfDsszBqVDFAOtpTDB8+nNNPP50TTjiBt73tbRx11FFvrZsyZQq33norJ510EscddxynnnpqlXtQe5WezroDWANclC1/Cvg34GP7GdMOHFOy3AJs7txJ0tnAtcCZEbGrZOwHOo1dVmGtZma5ufTS6kKjnLvuuqts++DBg3nggQfKruu47nHEEUewZs2at9qvuuqq2hbXjUrvE/mTiPjbiNiQva4H3tnNmBXAsZLGShoEXAwsLO0g6WTgNmBaRGwpWbUYOFfSYZIOA87N2szMrIFUGiKvSvpvHQuSTgf2e+UmInYDV1B8818PLIiItZJmS5qWdZsDHAzcK2m1pIXZ2BeBr1IMohUUT5292IP9MjOzXlDp6azPA3dm10YAfg98urtBEbEIWNSp7SslP5+9n7F3UDyNZmZmDarbEJHUDzguIv6rpLcDRMRLuVdmZmYNr9vTWRHxJsXTUkTESw4QMzPrUOk1kQclXSXpGEmHd7xyrczMzBpepSHyGeAvgUeAldmrNa+izMwOFKmPgge46aab2Fl652MdVPIU337AjIgY2+nV3Ud8zcyaT42fBd/XQ6SSp/i+Kel/Aaf1Qj1mZo2r41nwHW/cHc+Ch+Q7EEsfBX/OOedw5JFHsmDBAnbt2sVHP/pRrr/+el555RUuuugi2tvb2bNnD1/+8pd5/vnn2bx5M2eddRZHHHEES5cu7X5jOaj0I75LJH0c+EFE+PEiZtac6vAs+NJHwS9ZsoT77ruPX/7yl0QE06ZN45FHHmHr1q284x3v4Cc/+QlQfKbWsGHD+OY3v8nSpUs54ogjerqnNVPpNZErgQXALkkvSXpZkj+lZWYHlpyfBb9kyRKWLFnCySefzCmnnMITTzzBU089xYknnshDDz3El770JX72s58xbNiw7ifrJZUeiQwDLgXGRsRsSaOAo/Mry8ysDur8LPiI4JprruFzn/vcPutWrlzJokWLuOaaazj33HP5yle+UmaG3lfpkcgtwKnA9Gz5ZeBbuVRkZtaobrih+Oz3UlU+C770UfDnnXced9xxBzt27ADgueeeY8uWLWzevJmhQ4cyY8YMrrrqKlatWrXP2Hqp9EjkvRFxiqTHACLi99lDFc3MDhw5PAu+9FHwU6dO5ZJLLuG004qfYzr44IOZN28ebW1tfPGLX6Rfv34MHDiQ73znOwDMmjWLqVOncvTRR9ftwroquU4u6T8pfrvhiixMRgBLIuLkvAusVKFQiNZW37piZj2zfv16xo0bV+8yekW5fZW0MiIKqXNWejrrZuB+4EhJNwA/B/4hdaNmZtYcKv163PmSVgJ/SvFbBy+IiPW5VmZmZg2v0msiRMQTwBM51mJmZn1MpaezzMya1oFwD3Ve+5hriEiaIulJSW2Sri6z/gxJqyTtlnRhp3U3Slorab2kmyUpz1rN7MA0ZMgQtm3b1tRBEhFs27aNIUOG1Hzuik9n9ZSk/hTvLzkHaAdWSFoYEetKuj0LzASu6jT2fcDpwElZ08+BM4FledVrZgemlpYW2tvb2bp1a71LydWQIUNoaWmp+by5hQgwGWiLiA0Aku4BzgfeCpGI2Jite7PT2ACGAIMoXsgfCDyfY61mdoAaOHAgY8eOrXcZfVaep7NGAptKltuztm5FxKPAUuB32WuxPw1mZtZ48gyRctcwKjrpKOldwDighWLwfFDSGWX6zZLUKqm12Q9FzcwaUZ4h0g4cU7LcAmyucOxHgeURsSMidgAPUHx2114iYm5EFCKiMGLEiKoLNjOznskzRFYAx0oamz1n62JgYYVjnwXOlDRA0kCKF9V9OsvMrMHkFiIRsRu4AlhMMQAWRMRaSbMlTQOQNElSO/AJ4DZJa7Ph9wFPA78GfgX8KiJ+lFetZmaWpqIHMPYFfgCjmVnP9dYDGM3MzPbhEDEzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkDhEzM0vmEDEzs2QOETMzS+YQMTOzZA4RMzNLlmuISJoi6UlJbZKuLrP+DEmrJO2WdGGndaMkLZG0XtI6SWPyrNXMzHoutxCR1B+4BZgKjAemSxrfqduzwEzgrjJT3AnMiYhxwGRgS161mplZmgE5zj0ZaIuIDQCS7gHOB9Z1dIiIjdm6N0sHZmEzICIezPrtyLFOMzNLlOfprJHAppLl9qytEu8Gtkv6gaTHJM3Jjmz2ImmWpFZJrVu3bq1ByWZm1hN5hojKtEWFYwcA7weuAiYB76R42mvvySLmRkQhIgojRoxIrdPMzBLlGSLtwDElyy3A5h6MfSwiNkTEbuCHwCk1rs/MzKqUZ4isAI6VNFbSIOBiYGEPxh4mqePw4oOUXEsxM7PGkFuIZEcQVwCLgfXAgohYK2m2pGkAkiZJagc+AdwmaW02dg/FU1kPS/o1xVNjt+dVq5mZpVFEpZcpGluhUIjW1tZ6l2Fm1qdIWhkRhdTxvmPdzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLJlDxMzMkjlEzMwsmUPEzMySOUTMzCyZQ8TMzJI5RMzMLFmuISJpiqQnJbVJurrM+jMkrZK0W9KFZda/XdJzkr6VZ51mZpYmtxCR1B+4BZgKjAemSxrfqduzwEzgri6m+Srw07xqNDOz6uR5JDIZaIuIDRHxOnAPcH5ph4jYGBGPA292HixpInAUsCTHGs3MrAp5hshIYFPJcnvW1i1J/YBvAF/MoS4zM6uRPENEZdqiwrGXA4siYtP+OkmaJalVUuvWrVt7XKCZmVVnQI5ztwPHlCy3AJsrHHsa8H5JlwMHA4Mk7YiIvS7OR8RcYC5AoVCoNKDMzKxG8gyRFcCxksYCzwEXA5dUMjAiLu34WdJMoNA5QMzMrP5yO50VEbuBK4DFwHpgQUSslTRb0jQASZMktQOfAG6TtDaveszMrPYU0RxngQqFQrS2tta7DDOzPkXSyogopI73HetmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWzCFiZmbJHCJmZpbMIWJmZskcImZmlswhYmZmyRwiZmaWLNcQkTRF0pOS2iTt8x3pks6QtErSbkkXlrRPkPSopLWSHpf0yTzrNDOzNLmFiKT+wC3AVGA8MF3S+E7dngVmAnd1at8J/HlEHA9MAW6SdGhetZqZWZoBOc49GWiLiA0Aku4BzgfWdXSIiI3ZujdLB0bEb0p+3ixpCzAC2J5jvWZm1kN5ns4aCWwqWW7P2npE0mRgEPB0mXWzJLVKal25EsaMgfnzK5t3/vxi/379emdcrefIY64856zndhpt25Vo9Pqq1ez71xd0/A5g4sSqJoqIXF7AJ4B/LVn+FPAvXfT9LnBhmfajgSeBU7vf3sSAiKFDI+bNi/2aN6/YD/74ynNcrefIY64856zndhpt25Vo9Pqq1ez71xfs/TuYGFHFe72i+GZdc5JOA/4uIs7Llq/JQutrZfp+F/hxRNxX0vZ2YBnwtYi4t/vtFQJaARg8GE49teu+y5fDrl37tuc1rtZz5DFXnnPWczuNtu1KNHp91Wr2/esL9v4dFIhoVepceZ7OWgEcK2mspEHAxcDCSgZm/e8H7qwkQDor9w+0kvV5jav1HHnMleec9dxOo227Eo1eX7Waff/6gpr+XVdzGNPdC/gQ8BuK1zOuzdpmA9OynydRvFbyCrANWJu1zwDeAFaXvCbsf1sT3zo0Hj16/4dyo0fvfSid97haz5HHXHnOWc/tNNq2K9Ho9VWr2fevL9j7d1Dd6axcQ6Q3X74mUt1cec5Zz+002rYr0ej1VavZ968vqOU1kbq/+dfqBRNj9OjK/yHOm1dMYyl6ZVyt58hjrjznrOd2Gm3blWj0+qrV7PvXF3T8Dhr2wnpvKxQK0draWu8yzMz6FEkrI6KQOt7PzjIzs2QOETMzS+YQMTOzZA4RMzNL5hAxM7NkTfPpLElbgWcShh4BvJAwbhjwh4RxtZ4jj7nynLOe22m0bVei0eurVrPvX19wXEQckjo4z0fB96qIGJEyTlJrysfbJM2NiFkp26zlHHnMleec9dxOo227Eo1eX7Waff/6AklV3Rvh01npftQgc+QxV55z1nM7jbbtSjR6fdVq9v1rek1zOitV6pGImVkzqPY90EciMLfeBZiZ1VFV74EH/JGImZml85GImZklc4iYmVkyh4iZmSVziJSQNE7SrZLuk/SFetdjZtbbJF0g6XZJ/0fSud31b/oQkXSHpC2S1nRqnyLpSUltkq4GiIj1EfF54CLAH/s1s6bQw/fBH0bEZ4GZwCe7m7vpQwT4LjCltEFSf+AWYCowHpguaXy2bhrwc+Dh3i3TzCw336UH74OZ67L1+9X0IRIRjwAvdmqeDLRFxIaIeB24Bzg/678wIt4HXNq7lZqZ5aMn74Mq+kfggYhY1d3cTfPsrB4aCWwqWW4H3ivpA8DHgMHAojrUZWbWW8q+DwJ/BZwNDJP0roi4dX+THKghojJtERHLgGW9W4qZWV109T54M3BzpZM0/emsLrQDx5QstwCb61SLmVk91OR98EANkRXAsZLGShoEXAwsrHNNZma9qSbvg00fIpLuBh4FjpPULumyiNgNXAEsBtYDCyJibT3rNDPLS57vg34Ao5mZJWv6IxEzM8uPQ8TMzJI5RMzMLJlDxMzMkjlEzMwsmUPEzMySOUTMqiDpUEmX17sOs3pxiJhV51BgnxDJHrNt1vQO1AcwmtXK14E/kbQaeAPYAfwOmACMlzQD+B/AIOA/gcsjYk/2jXHXU3xi9NPAX0TEjnrsgFk1fMe6WRUkjQF+HBEnZF8l8BPghIj4raRxwI3AxyLiDUnfBpZT/JqBHwBTI+IVSV8CBkfE7LrshFkVfCRiVlu/jIjfZj//KTARWCEJ4G3AFuBUit8k94usfRDF5xqZ9TkOEbPaeqXkZwH/HhHXlHaQ9BHgwYiY3quVmeXAF9bNqvMycEgX6x4GLpR0JICkwyWNpnhK63RJ78rah0p6d69Ua1ZjPhIxq0JEbJP0C0lrgFeB50vWrZN0HbBEUj+KF97/MiKWS5oJ3C1pcNb9OuA3vVy+WdV8Yd3MzJL5dJaZmSVziJiZWTKHiJmZJXOImJlZMoeImZklc4iYmVkyh4iZmSVziJiZWbL/D0amXDGAdIwWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.semilogx(tree_arr_gbt, train_err, 'b-o', label = 'train')\n",
    "plt.semilogx(tree_arr_gbt, test_err, 'r-o', label = 'test')\n",
    "plt.xlim([np.max(tree_arr_gbt), np.min(tree_arr_gbt)])\n",
    "plt.title('Error vs. tree')\n",
    "plt.xlabel('tree')\n",
    "plt.ylabel('error')\n",
    "plt.legend()\n",
    "pass"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Наблюдается ярковыраженное переобучение во всех 4х примерах."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}