RF_draft1.ipynb

{
 "cells": [
  {
   "cell_type": "code",
<<<<<<< HEAD
   "execution_count": 5,
=======
   "execution_count": 76,
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import pandas as pd\n",
    "import sys\n",
    "\n",
    "import csv\n",
    "from scipy.stats import skew, kurtosis\n",
    "from numpy import nan\n",
    "import math\n"
   ]
  },
  {
   "cell_type": "code",
<<<<<<< HEAD
   "execution_count": 6,
=======
   "execution_count": 77,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "import pyeeg as p"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Requirement already satisfied: mne in c:\\users\\aai\\anaconda3\\lib\\site-packages (0.21.0)\n",
      "Requirement already satisfied: numpy>=1.11.3 in c:\\users\\aai\\anaconda3\\lib\\site-packages (from mne) (1.18.5)\n",
      "Requirement already satisfied: scipy>=0.17.1 in c:\\users\\aai\\anaconda3\\lib\\site-packages (from mne) (1.5.0)\n",
      "Note: you may need to restart the kernel to use updated packages.\n"
     ]
    }
   ],
   "source": [
    "pip install mne\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {},
   "outputs": [],
   "source": [
    "import mne"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
   "metadata": {},
   "outputs": [],
   "source": [
    "def all_psd(data):\n",
    "    fs = 256                                \n",
    "    N = data.shape[1] # total num of points \n",
    "\n",
    "    # Get only in postive frequencies\n",
    "    fft_vals = np.absolute(np.fft.rfft(data))\n",
    "\n",
    "    n_rows = fft_vals.shape[0]\n",
    "    n_cols = fft_vals.shape[1]\n",
    "    psd_vals = np.zeros(shape=(n_rows, n_cols))\n",
    "\n",
    "    for i in range(n_rows):\n",
    "        for j in range(n_cols):\n",
    "            psd_vals[i][j] = (N/fs) * fft_vals[i][j] * fft_vals[i][j];\n",
    "\n",
    "\n",
    "    # Get frequencies for amplitudes in Hz\n",
    "    fft_freq = np.fft.rfftfreq(data.shape[1], 1.0/fs)\n",
    "\n",
    "    # Define EEG bands\n",
    "    eeg_bands = {'Delta': (0, 4),\n",
    "                 'Theta': (4, 8),\n",
    "                 'Alpha': (8, 12),\n",
    "                 'Beta': (12, 30),\n",
    "                 'Gamma': (30, 45)}\n",
    "\n",
    "    # Take the mean of the fft amplitude for each EEG band\n",
    "    eeg_band_fft = dict()\n",
    "    psd_vals_list = []\n",
    "    for band in eeg_bands:  \n",
    "        freq_ix = np.where((fft_freq >= eeg_bands[band][0]) & \n",
    "                           (fft_freq <= eeg_bands[band][1]))[0]\n",
    "        eeg_band_fft[band] = psd_vals[:,freq_ix]\n",
    "        for i in range(len(eeg_band_fft[band])): \n",
    "            psd_vals_list.append(np.mean(eeg_band_fft[band][i]) * 1000000)\n",
    "    return psd_vals_list"
   ]
  },
  {
   "cell_type": "code",
<<<<<<< HEAD
   "execution_count": 13,
=======
   "execution_count": 99,
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
   "metadata": {},
   "outputs": [],
   "source": [
    "seizure_var=0\n",
    "def add_row(ip_psd, start, end):\n",
    "    row_to_add = []\n",
    "    \n",
    "    ip_psd = ip_psd[:,start:end]\n",
    "    ip_psd = ip_psd[:][0]\n",
    "    \n",
    "    \n",
    "    p = all_psd(ip_psd)\n",
    "    p.append(seizure_var)\n",
    "    \n",
    "    return p"
   ]
  },
  {
   "cell_type": "code",
<<<<<<< HEAD
   "execution_count": 8,
=======
   "execution_count": 100,
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Extracting EDF parameters from C:\\Users\\AAI\\Desktop\\btp\\chb01_03.edf...\n",
      "EDF file detected\n",
      "Setting channel info structure...\n",
      "Creating raw.info structure...\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
<<<<<<< HEAD
      "<ipython-input-8-e0b81f09e547>:1: RuntimeWarning: Channel names are not unique, found duplicates for: {'T8-P8'}. Applying running numbers for duplicates.\n",
=======
      "<ipython-input-100-e0b81f09e547>:1: RuntimeWarning: Channel names are not unique, found duplicates for: {'T8-P8'}. Applying running numbers for duplicates.\n",
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
      "  data = mne.io.read_raw_edf('chb01_03.edf')\n"
     ]
    }
   ],
   "source": [
    "data = mne.io.read_raw_edf('chb01_03.edf')\n",
    "header = ','.join(data.ch_names)\n",
    "df = pd.DataFrame(data[:][0])\n",
    "df = df.transpose()"
   ]
  },
  {
   "cell_type": "code",
<<<<<<< HEAD
   "execution_count": 14,
=======
   "execution_count": 102,
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(256000, 23)"
      ]
     },
<<<<<<< HEAD
     "execution_count": 14,
=======
     "execution_count": 102,
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "start = temp = 1996*256 \n",
    "duration = 10*256\n",
    "end = 2996*256\n",
    "\n",
    "df1 = df.iloc[start:end, :]\n",
    "df1 = pd.DataFrame(df1)\n",
    "df1.shape"
   ]
  },
  {
   "cell_type": "code",
<<<<<<< HEAD
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
=======
   "execution_count": 103,
   "metadata": {},
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
   "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>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>7</th>\n",
       "      <th>8</th>\n",
       "      <th>9</th>\n",
       "      <th>...</th>\n",
       "      <th>106</th>\n",
       "      <th>107</th>\n",
       "      <th>108</th>\n",
       "      <th>109</th>\n",
       "      <th>110</th>\n",
       "      <th>111</th>\n",
       "      <th>112</th>\n",
       "      <th>113</th>\n",
       "      <th>114</th>\n",
       "      <th>115</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>178.786422</td>\n",
       "      <td>187.299138</td>\n",
       "      <td>373.456785</td>\n",
       "      <td>147.897919</td>\n",
       "      <td>433.100880</td>\n",
       "      <td>439.940230</td>\n",
       "      <td>333.774992</td>\n",
       "      <td>416.115508</td>\n",
       "      <td>394.691618</td>\n",
       "      <td>368.760198</td>\n",
       "      <td>...</td>\n",
       "      <td>0.410709</td>\n",
       "      <td>3.580521</td>\n",
       "      <td>0.574664</td>\n",
       "      <td>0.676765</td>\n",
       "      <td>0.389804</td>\n",
       "      <td>3.855179</td>\n",
       "      <td>8.851754</td>\n",
       "      <td>0.729672</td>\n",
       "      <td>0.410709</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>246.545195</td>\n",
       "      <td>146.974897</td>\n",
       "      <td>407.334715</td>\n",
       "      <td>258.835636</td>\n",
       "      <td>733.011971</td>\n",
       "      <td>259.251316</td>\n",
       "      <td>281.086734</td>\n",
       "      <td>447.975677</td>\n",
       "      <td>726.113512</td>\n",
       "      <td>372.659878</td>\n",
       "      <td>...</td>\n",
       "      <td>0.408117</td>\n",
       "      <td>3.408203</td>\n",
       "      <td>0.704777</td>\n",
       "      <td>0.795886</td>\n",
       "      <td>0.492413</td>\n",
       "      <td>4.095119</td>\n",
       "      <td>8.313062</td>\n",
       "      <td>0.824171</td>\n",
       "      <td>0.408117</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>154.673070</td>\n",
       "      <td>152.911484</td>\n",
       "      <td>346.467292</td>\n",
       "      <td>296.324842</td>\n",
       "      <td>448.180389</td>\n",
       "      <td>407.671629</td>\n",
       "      <td>272.390824</td>\n",
       "      <td>534.975166</td>\n",
       "      <td>525.624321</td>\n",
       "      <td>372.344138</td>\n",
       "      <td>...</td>\n",
       "      <td>0.321299</td>\n",
       "      <td>3.311277</td>\n",
       "      <td>0.696812</td>\n",
       "      <td>0.748675</td>\n",
       "      <td>0.366602</td>\n",
       "      <td>3.904908</td>\n",
       "      <td>8.709865</td>\n",
       "      <td>0.767086</td>\n",
       "      <td>0.321299</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>142.616877</td>\n",
       "      <td>123.756485</td>\n",
       "      <td>274.717083</td>\n",
       "      <td>262.639245</td>\n",
       "      <td>786.032023</td>\n",
       "      <td>453.893509</td>\n",
       "      <td>325.054145</td>\n",
       "      <td>585.962765</td>\n",
       "      <td>691.101114</td>\n",
       "      <td>508.397870</td>\n",
       "      <td>...</td>\n",
       "      <td>0.395649</td>\n",
       "      <td>3.612864</td>\n",
       "      <td>0.795568</td>\n",
       "      <td>0.716539</td>\n",
       "      <td>0.344873</td>\n",
       "      <td>3.920525</td>\n",
       "      <td>8.164512</td>\n",
       "      <td>0.807174</td>\n",
       "      <td>0.395649</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>260.544638</td>\n",
       "      <td>196.312737</td>\n",
       "      <td>576.709094</td>\n",
       "      <td>318.495857</td>\n",
       "      <td>711.977651</td>\n",
       "      <td>452.092105</td>\n",
       "      <td>368.535053</td>\n",
       "      <td>774.207210</td>\n",
       "      <td>544.819770</td>\n",
       "      <td>521.200692</td>\n",
       "      <td>...</td>\n",
       "      <td>0.350373</td>\n",
       "      <td>3.773968</td>\n",
       "      <td>0.753706</td>\n",
       "      <td>0.902003</td>\n",
       "      <td>0.380239</td>\n",
       "      <td>4.014634</td>\n",
       "      <td>7.777630</td>\n",
       "      <td>0.833583</td>\n",
       "      <td>0.350373</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>234.320907</td>\n",
       "      <td>261.550282</td>\n",
       "      <td>470.721632</td>\n",
       "      <td>286.447461</td>\n",
       "      <td>751.604378</td>\n",
       "      <td>469.626384</td>\n",
       "      <td>357.002351</td>\n",
       "      <td>577.517453</td>\n",
       "      <td>785.321363</td>\n",
       "      <td>489.468379</td>\n",
       "      <td>...</td>\n",
       "      <td>0.346077</td>\n",
       "      <td>3.542306</td>\n",
       "      <td>0.545660</td>\n",
       "      <td>0.555900</td>\n",
       "      <td>0.331486</td>\n",
       "      <td>3.823273</td>\n",
       "      <td>8.410936</td>\n",
       "      <td>0.756244</td>\n",
       "      <td>0.346077</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>273.379424</td>\n",
       "      <td>259.283579</td>\n",
       "      <td>521.113884</td>\n",
       "      <td>279.343219</td>\n",
       "      <td>571.597933</td>\n",
       "      <td>461.339238</td>\n",
       "      <td>522.689588</td>\n",
       "      <td>474.163002</td>\n",
       "      <td>568.110176</td>\n",
       "      <td>406.117779</td>\n",
       "      <td>...</td>\n",
       "      <td>0.352900</td>\n",
       "      <td>3.593283</td>\n",
       "      <td>0.635650</td>\n",
       "      <td>0.667032</td>\n",
       "      <td>0.392782</td>\n",
       "      <td>3.953978</td>\n",
       "      <td>8.084922</td>\n",
       "      <td>0.754650</td>\n",
       "      <td>0.352900</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>312.157337</td>\n",
       "      <td>327.162636</td>\n",
       "      <td>605.796649</td>\n",
       "      <td>393.135495</td>\n",
       "      <td>1041.770233</td>\n",
       "      <td>647.134773</td>\n",
       "      <td>511.948216</td>\n",
       "      <td>713.981772</td>\n",
       "      <td>1216.144852</td>\n",
       "      <td>716.313123</td>\n",
       "      <td>...</td>\n",
       "      <td>0.339203</td>\n",
       "      <td>3.759761</td>\n",
       "      <td>0.659787</td>\n",
       "      <td>0.565934</td>\n",
       "      <td>0.356684</td>\n",
       "      <td>3.929324</td>\n",
       "      <td>8.775190</td>\n",
       "      <td>0.805310</td>\n",
       "      <td>0.339203</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>90.867944</td>\n",
       "      <td>183.301948</td>\n",
       "      <td>179.605055</td>\n",
       "      <td>206.120312</td>\n",
       "      <td>340.883022</td>\n",
       "      <td>351.363763</td>\n",
       "      <td>279.098579</td>\n",
       "      <td>504.991956</td>\n",
       "      <td>267.143485</td>\n",
       "      <td>271.414462</td>\n",
       "      <td>...</td>\n",
       "      <td>0.363857</td>\n",
       "      <td>3.533879</td>\n",
       "      <td>0.541001</td>\n",
       "      <td>0.655530</td>\n",
       "      <td>0.331713</td>\n",
       "      <td>3.923096</td>\n",
       "      <td>8.709387</td>\n",
       "      <td>0.794712</td>\n",
       "      <td>0.363857</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>121.823402</td>\n",
       "      <td>143.578508</td>\n",
       "      <td>188.337742</td>\n",
       "      <td>315.325243</td>\n",
       "      <td>290.538303</td>\n",
       "      <td>463.355738</td>\n",
       "      <td>185.048609</td>\n",
       "      <td>490.961345</td>\n",
       "      <td>216.813588</td>\n",
       "      <td>398.905451</td>\n",
       "      <td>...</td>\n",
       "      <td>0.355697</td>\n",
       "      <td>3.714707</td>\n",
       "      <td>0.488094</td>\n",
       "      <td>0.660672</td>\n",
       "      <td>0.367932</td>\n",
       "      <td>4.054801</td>\n",
       "      <td>8.720214</td>\n",
       "      <td>0.790561</td>\n",
       "      <td>0.355697</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>201.869968</td>\n",
       "      <td>207.343892</td>\n",
       "      <td>358.151280</td>\n",
       "      <td>313.966200</td>\n",
       "      <td>575.218551</td>\n",
       "      <td>700.470319</td>\n",
       "      <td>288.476645</td>\n",
       "      <td>446.116792</td>\n",
       "      <td>485.959500</td>\n",
       "      <td>569.900688</td>\n",
       "      <td>...</td>\n",
       "      <td>0.352160</td>\n",
       "      <td>3.532321</td>\n",
       "      <td>0.597552</td>\n",
       "      <td>0.723606</td>\n",
       "      <td>0.449655</td>\n",
       "      <td>4.010461</td>\n",
       "      <td>8.305585</td>\n",
       "      <td>0.858595</td>\n",
       "      <td>0.352160</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>127.500698</td>\n",
       "      <td>140.006731</td>\n",
       "      <td>311.099870</td>\n",
       "      <td>335.134552</td>\n",
       "      <td>385.807141</td>\n",
       "      <td>420.605587</td>\n",
       "      <td>234.936613</td>\n",
       "      <td>432.448664</td>\n",
       "      <td>297.271949</td>\n",
       "      <td>323.821490</td>\n",
       "      <td>...</td>\n",
       "      <td>0.397548</td>\n",
       "      <td>3.479982</td>\n",
       "      <td>0.665714</td>\n",
       "      <td>0.629633</td>\n",
       "      <td>0.393463</td>\n",
       "      <td>3.792875</td>\n",
       "      <td>7.799580</td>\n",
       "      <td>0.731764</td>\n",
       "      <td>0.397548</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>117.045834</td>\n",
       "      <td>115.357530</td>\n",
       "      <td>134.405483</td>\n",
       "      <td>234.935440</td>\n",
       "      <td>235.938726</td>\n",
       "      <td>237.162789</td>\n",
       "      <td>148.821140</td>\n",
       "      <td>320.215899</td>\n",
       "      <td>181.596172</td>\n",
       "      <td>181.022392</td>\n",
       "      <td>...</td>\n",
       "      <td>0.366877</td>\n",
       "      <td>3.452338</td>\n",
       "      <td>0.547187</td>\n",
       "      <td>0.605098</td>\n",
       "      <td>0.378353</td>\n",
       "      <td>3.587219</td>\n",
       "      <td>7.997030</td>\n",
       "      <td>0.695318</td>\n",
       "      <td>0.366877</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>195.575947</td>\n",
       "      <td>199.071527</td>\n",
       "      <td>274.309503</td>\n",
       "      <td>282.684451</td>\n",
       "      <td>524.998364</td>\n",
       "      <td>451.207193</td>\n",
       "      <td>293.616764</td>\n",
       "      <td>448.114314</td>\n",
       "      <td>343.557154</td>\n",
       "      <td>355.621789</td>\n",
       "      <td>...</td>\n",
       "      <td>0.312612</td>\n",
       "      <td>3.511218</td>\n",
       "      <td>0.571026</td>\n",
       "      <td>0.620062</td>\n",
       "      <td>0.364408</td>\n",
       "      <td>3.688348</td>\n",
       "      <td>8.361690</td>\n",
       "      <td>0.753523</td>\n",
       "      <td>0.312612</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>195.393033</td>\n",
       "      <td>270.996408</td>\n",
       "      <td>296.393823</td>\n",
       "      <td>357.205989</td>\n",
       "      <td>349.385642</td>\n",
       "      <td>385.753793</td>\n",
       "      <td>249.576777</td>\n",
       "      <td>490.626714</td>\n",
       "      <td>349.490917</td>\n",
       "      <td>266.248300</td>\n",
       "      <td>...</td>\n",
       "      <td>0.405448</td>\n",
       "      <td>3.814354</td>\n",
       "      <td>0.570337</td>\n",
       "      <td>0.706242</td>\n",
       "      <td>0.392385</td>\n",
       "      <td>3.601962</td>\n",
       "      <td>8.723980</td>\n",
       "      <td>0.745340</td>\n",
       "      <td>0.405448</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>176.691208</td>\n",
       "      <td>177.027086</td>\n",
       "      <td>398.185835</td>\n",
       "      <td>223.864522</td>\n",
       "      <td>517.424334</td>\n",
       "      <td>362.294943</td>\n",
       "      <td>281.813067</td>\n",
       "      <td>438.809774</td>\n",
       "      <td>493.636955</td>\n",
       "      <td>436.563591</td>\n",
       "      <td>...</td>\n",
       "      <td>0.345647</td>\n",
       "      <td>3.607638</td>\n",
       "      <td>0.752208</td>\n",
       "      <td>0.576956</td>\n",
       "      <td>0.341811</td>\n",
       "      <td>3.690010</td>\n",
       "      <td>8.485560</td>\n",
       "      <td>0.796332</td>\n",
       "      <td>0.345647</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>280.072969</td>\n",
       "      <td>163.121439</td>\n",
       "      <td>317.006730</td>\n",
       "      <td>188.162560</td>\n",
       "      <td>912.977836</td>\n",
       "      <td>375.896532</td>\n",
       "      <td>256.447720</td>\n",
       "      <td>443.113740</td>\n",
       "      <td>524.124246</td>\n",
       "      <td>264.714965</td>\n",
       "      <td>...</td>\n",
       "      <td>0.357802</td>\n",
       "      <td>3.766675</td>\n",
       "      <td>0.687498</td>\n",
       "      <td>0.711523</td>\n",
       "      <td>0.374456</td>\n",
       "      <td>3.939251</td>\n",
       "      <td>7.935518</td>\n",
       "      <td>0.798783</td>\n",
       "      <td>0.357802</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>173.215078</td>\n",
       "      <td>149.073981</td>\n",
       "      <td>150.037932</td>\n",
       "      <td>152.980695</td>\n",
       "      <td>568.715017</td>\n",
       "      <td>366.064307</td>\n",
       "      <td>247.497097</td>\n",
       "      <td>331.551166</td>\n",
       "      <td>472.490192</td>\n",
       "      <td>244.326979</td>\n",
       "      <td>...</td>\n",
       "      <td>0.422197</td>\n",
       "      <td>3.184529</td>\n",
       "      <td>0.604299</td>\n",
       "      <td>0.719281</td>\n",
       "      <td>0.390349</td>\n",
       "      <td>3.709922</td>\n",
       "      <td>7.721513</td>\n",
       "      <td>0.792122</td>\n",
       "      <td>0.422197</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>172.190633</td>\n",
       "      <td>202.396157</td>\n",
       "      <td>299.976177</td>\n",
       "      <td>160.355398</td>\n",
       "      <td>689.230772</td>\n",
       "      <td>564.925225</td>\n",
       "      <td>340.747540</td>\n",
       "      <td>437.707263</td>\n",
       "      <td>682.712240</td>\n",
       "      <td>824.878818</td>\n",
       "      <td>...</td>\n",
       "      <td>0.364102</td>\n",
       "      <td>3.716139</td>\n",
       "      <td>0.672627</td>\n",
       "      <td>0.819458</td>\n",
       "      <td>0.404210</td>\n",
       "      <td>3.872358</td>\n",
       "      <td>8.345090</td>\n",
       "      <td>0.848888</td>\n",
       "      <td>0.364102</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>253.430027</td>\n",
       "      <td>268.879732</td>\n",
       "      <td>390.696573</td>\n",
       "      <td>340.222539</td>\n",
       "      <td>777.303251</td>\n",
       "      <td>626.193163</td>\n",
       "      <td>366.015566</td>\n",
       "      <td>605.850189</td>\n",
       "      <td>588.125692</td>\n",
       "      <td>589.722847</td>\n",
       "      <td>...</td>\n",
       "      <td>0.452971</td>\n",
       "      <td>3.702314</td>\n",
       "      <td>0.557211</td>\n",
       "      <td>0.682076</td>\n",
       "      <td>0.457097</td>\n",
       "      <td>3.775675</td>\n",
       "      <td>8.739288</td>\n",
       "      <td>0.772543</td>\n",
       "      <td>0.452971</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>299.869998</td>\n",
       "      <td>227.688413</td>\n",
       "      <td>358.779289</td>\n",
       "      <td>151.046274</td>\n",
       "      <td>992.183611</td>\n",
       "      <td>517.627875</td>\n",
       "      <td>364.925738</td>\n",
       "      <td>474.377177</td>\n",
       "      <td>1047.339851</td>\n",
       "      <td>607.868023</td>\n",
       "      <td>...</td>\n",
       "      <td>0.419587</td>\n",
       "      <td>3.258614</td>\n",
       "      <td>0.618107</td>\n",
       "      <td>0.667598</td>\n",
       "      <td>0.411717</td>\n",
       "      <td>3.785812</td>\n",
       "      <td>8.555000</td>\n",
       "      <td>0.867401</td>\n",
       "      <td>0.419587</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>176.532272</td>\n",
       "      <td>106.073423</td>\n",
       "      <td>295.705094</td>\n",
       "      <td>246.425916</td>\n",
       "      <td>381.403013</td>\n",
       "      <td>274.749131</td>\n",
       "      <td>276.640917</td>\n",
       "      <td>336.243612</td>\n",
       "      <td>305.400304</td>\n",
       "      <td>152.363553</td>\n",
       "      <td>...</td>\n",
       "      <td>0.389708</td>\n",
       "      <td>3.440791</td>\n",
       "      <td>0.593434</td>\n",
       "      <td>0.684052</td>\n",
       "      <td>0.388469</td>\n",
       "      <td>3.629259</td>\n",
       "      <td>8.198102</td>\n",
       "      <td>0.838601</td>\n",
       "      <td>0.389708</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>228.204807</td>\n",
       "      <td>316.282246</td>\n",
       "      <td>341.836231</td>\n",
       "      <td>280.807598</td>\n",
       "      <td>711.831743</td>\n",
       "      <td>758.850904</td>\n",
       "      <td>410.210644</td>\n",
       "      <td>612.924933</td>\n",
       "      <td>398.635739</td>\n",
       "      <td>585.706312</td>\n",
       "      <td>...</td>\n",
       "      <td>0.347397</td>\n",
       "      <td>3.557350</td>\n",
       "      <td>0.892613</td>\n",
       "      <td>0.688156</td>\n",
       "      <td>0.383629</td>\n",
       "      <td>4.118524</td>\n",
       "      <td>8.289979</td>\n",
       "      <td>0.772936</td>\n",
       "      <td>0.347397</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>166.710607</td>\n",
       "      <td>202.597143</td>\n",
       "      <td>392.102092</td>\n",
       "      <td>274.927670</td>\n",
       "      <td>743.108185</td>\n",
       "      <td>633.448967</td>\n",
       "      <td>392.219767</td>\n",
       "      <td>812.956043</td>\n",
       "      <td>591.870036</td>\n",
       "      <td>472.259010</td>\n",
       "      <td>...</td>\n",
       "      <td>0.442770</td>\n",
       "      <td>3.740394</td>\n",
       "      <td>0.653449</td>\n",
       "      <td>0.603795</td>\n",
       "      <td>0.363151</td>\n",
       "      <td>3.721382</td>\n",
       "      <td>8.873737</td>\n",
       "      <td>0.761885</td>\n",
       "      <td>0.442770</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>313.305186</td>\n",
       "      <td>188.154894</td>\n",
       "      <td>390.240316</td>\n",
       "      <td>262.611060</td>\n",
       "      <td>610.824073</td>\n",
       "      <td>672.614760</td>\n",
       "      <td>459.503370</td>\n",
       "      <td>559.197594</td>\n",
       "      <td>624.338175</td>\n",
       "      <td>533.952562</td>\n",
       "      <td>...</td>\n",
       "      <td>0.466055</td>\n",
       "      <td>3.497552</td>\n",
       "      <td>0.594804</td>\n",
       "      <td>0.740771</td>\n",
       "      <td>0.461666</td>\n",
       "      <td>3.822586</td>\n",
       "      <td>8.592550</td>\n",
       "      <td>0.775473</td>\n",
       "      <td>0.466055</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>192.454194</td>\n",
       "      <td>162.745680</td>\n",
       "      <td>200.520404</td>\n",
       "      <td>291.521940</td>\n",
       "      <td>585.399265</td>\n",
       "      <td>354.628730</td>\n",
       "      <td>312.605028</td>\n",
       "      <td>531.314156</td>\n",
       "      <td>524.096684</td>\n",
       "      <td>272.855330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.415696</td>\n",
       "      <td>3.346305</td>\n",
       "      <td>0.557133</td>\n",
       "      <td>0.573472</td>\n",
       "      <td>0.416878</td>\n",
       "      <td>3.943003</td>\n",
       "      <td>8.188199</td>\n",
       "      <td>0.833220</td>\n",
       "      <td>0.415696</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>289.190744</td>\n",
       "      <td>299.188237</td>\n",
       "      <td>301.265454</td>\n",
       "      <td>262.134823</td>\n",
       "      <td>692.545951</td>\n",
       "      <td>669.218277</td>\n",
       "      <td>438.199987</td>\n",
       "      <td>714.247490</td>\n",
       "      <td>667.919096</td>\n",
       "      <td>548.131484</td>\n",
       "      <td>...</td>\n",
       "      <td>0.340344</td>\n",
       "      <td>3.520040</td>\n",
       "      <td>0.690583</td>\n",
       "      <td>0.593332</td>\n",
       "      <td>0.421714</td>\n",
       "      <td>3.918895</td>\n",
       "      <td>8.064142</td>\n",
       "      <td>0.741595</td>\n",
       "      <td>0.340344</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>262.605588</td>\n",
       "      <td>244.991274</td>\n",
       "      <td>667.564831</td>\n",
       "      <td>408.848913</td>\n",
       "      <td>473.418976</td>\n",
       "      <td>532.245267</td>\n",
       "      <td>331.669166</td>\n",
       "      <td>479.071785</td>\n",
       "      <td>445.218419</td>\n",
       "      <td>397.654977</td>\n",
       "      <td>...</td>\n",
       "      <td>0.289427</td>\n",
       "      <td>3.392331</td>\n",
       "      <td>0.601149</td>\n",
       "      <td>0.613482</td>\n",
       "      <td>0.391231</td>\n",
       "      <td>3.736227</td>\n",
       "      <td>7.861983</td>\n",
       "      <td>0.797338</td>\n",
       "      <td>0.289427</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>154.738055</td>\n",
       "      <td>222.056288</td>\n",
       "      <td>389.591496</td>\n",
       "      <td>290.912424</td>\n",
       "      <td>472.113157</td>\n",
       "      <td>501.991482</td>\n",
       "      <td>424.592671</td>\n",
       "      <td>657.931235</td>\n",
       "      <td>453.669150</td>\n",
       "      <td>444.013337</td>\n",
       "      <td>...</td>\n",
       "      <td>0.409222</td>\n",
       "      <td>3.558121</td>\n",
       "      <td>0.691285</td>\n",
       "      <td>0.729957</td>\n",
       "      <td>0.362660</td>\n",
       "      <td>3.640856</td>\n",
       "      <td>8.568836</td>\n",
       "      <td>0.791814</td>\n",
       "      <td>0.409222</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>147.971812</td>\n",
       "      <td>155.431087</td>\n",
       "      <td>334.884672</td>\n",
       "      <td>400.147118</td>\n",
       "      <td>203.585776</td>\n",
       "      <td>277.576426</td>\n",
       "      <td>315.079381</td>\n",
       "      <td>417.870902</td>\n",
       "      <td>209.603915</td>\n",
       "      <td>273.144925</td>\n",
       "      <td>...</td>\n",
       "      <td>0.395049</td>\n",
       "      <td>3.560522</td>\n",
       "      <td>0.789697</td>\n",
       "      <td>0.962958</td>\n",
       "      <td>0.445118</td>\n",
       "      <td>3.925553</td>\n",
       "      <td>8.646597</td>\n",
       "      <td>0.771469</td>\n",
       "      <td>0.395049</td>\n",
       "      <td>0.0</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70</th>\n",
       "      <td>261.134719</td>\n",
       "      <td>258.302304</td>\n",
       "      <td>507.139181</td>\n",
       "      <td>350.101826</td>\n",
       "      <td>507.061498</td>\n",
       "      <td>665.581335</td>\n",
       "      <td>251.476918</td>\n",
       "      <td>555.539360</td>\n",
       "      <td>430.354020</td>\n",
       "      <td>432.256490</td>\n",
       "      <td>...</td>\n",
       "      <td>0.425885</td>\n",
       "      <td>3.683649</td>\n",
       "      <td>0.957510</td>\n",
       "      <td>0.740335</td>\n",
       "      <td>0.372962</td>\n",
       "      <td>4.348339</td>\n",
       "      <td>8.573490</td>\n",
       "      <td>0.876375</td>\n",
       "      <td>0.425885</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>71</th>\n",
       "      <td>230.820184</td>\n",
       "      <td>227.389183</td>\n",
       "      <td>388.022951</td>\n",
       "      <td>117.647582</td>\n",
       "      <td>291.595835</td>\n",
       "      <td>223.590412</td>\n",
       "      <td>135.671843</td>\n",
       "      <td>211.902017</td>\n",
       "      <td>245.123641</td>\n",
       "      <td>269.888973</td>\n",
       "      <td>...</td>\n",
       "      <td>0.372161</td>\n",
       "      <td>3.775336</td>\n",
       "      <td>0.728939</td>\n",
       "      <td>0.784113</td>\n",
       "      <td>0.471682</td>\n",
       "      <td>3.614408</td>\n",
       "      <td>7.997206</td>\n",
       "      <td>0.667619</td>\n",
       "      <td>0.372161</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72</th>\n",
       "      <td>199.441131</td>\n",
       "      <td>216.596367</td>\n",
       "      <td>500.309393</td>\n",
       "      <td>128.812704</td>\n",
       "      <td>337.809754</td>\n",
       "      <td>215.469931</td>\n",
       "      <td>163.186875</td>\n",
       "      <td>249.513322</td>\n",
       "      <td>221.490544</td>\n",
       "      <td>330.834000</td>\n",
       "      <td>...</td>\n",
       "      <td>0.377273</td>\n",
       "      <td>3.630585</td>\n",
       "      <td>0.533531</td>\n",
       "      <td>0.768592</td>\n",
       "      <td>0.326154</td>\n",
       "      <td>3.480983</td>\n",
       "      <td>8.328907</td>\n",
       "      <td>0.758970</td>\n",
       "      <td>0.377273</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73</th>\n",
       "      <td>138.730279</td>\n",
       "      <td>191.174348</td>\n",
       "      <td>212.830176</td>\n",
       "      <td>166.411185</td>\n",
       "      <td>397.164134</td>\n",
       "      <td>298.816490</td>\n",
       "      <td>261.386305</td>\n",
       "      <td>284.328117</td>\n",
       "      <td>351.088574</td>\n",
       "      <td>288.612735</td>\n",
       "      <td>...</td>\n",
       "      <td>0.417536</td>\n",
       "      <td>3.335675</td>\n",
       "      <td>0.627466</td>\n",
       "      <td>0.738017</td>\n",
       "      <td>0.383787</td>\n",
       "      <td>3.360009</td>\n",
       "      <td>8.394821</td>\n",
       "      <td>0.674719</td>\n",
       "      <td>0.417536</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74</th>\n",
       "      <td>124.189176</td>\n",
       "      <td>166.689290</td>\n",
       "      <td>311.292913</td>\n",
       "      <td>179.981363</td>\n",
       "      <td>314.544481</td>\n",
       "      <td>449.685500</td>\n",
       "      <td>259.772225</td>\n",
       "      <td>359.374297</td>\n",
       "      <td>248.261614</td>\n",
       "      <td>543.192428</td>\n",
       "      <td>...</td>\n",
       "      <td>0.411413</td>\n",
       "      <td>3.793024</td>\n",
       "      <td>0.489882</td>\n",
       "      <td>0.649907</td>\n",
       "      <td>0.437163</td>\n",
       "      <td>3.789730</td>\n",
       "      <td>8.660178</td>\n",
       "      <td>0.737800</td>\n",
       "      <td>0.411413</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75</th>\n",
       "      <td>177.371730</td>\n",
       "      <td>216.114620</td>\n",
       "      <td>465.725745</td>\n",
       "      <td>294.850166</td>\n",
       "      <td>271.478315</td>\n",
       "      <td>244.936852</td>\n",
       "      <td>146.350015</td>\n",
       "      <td>419.355925</td>\n",
       "      <td>115.308624</td>\n",
       "      <td>283.393541</td>\n",
       "      <td>...</td>\n",
       "      <td>0.420551</td>\n",
       "      <td>3.995588</td>\n",
       "      <td>0.574573</td>\n",
       "      <td>0.761602</td>\n",
       "      <td>0.482047</td>\n",
       "      <td>4.297029</td>\n",
       "      <td>8.685668</td>\n",
       "      <td>0.904984</td>\n",
       "      <td>0.420551</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76</th>\n",
       "      <td>212.262949</td>\n",
       "      <td>256.926485</td>\n",
       "      <td>290.546175</td>\n",
       "      <td>246.348148</td>\n",
       "      <td>481.660190</td>\n",
       "      <td>569.281736</td>\n",
       "      <td>239.155450</td>\n",
       "      <td>348.768704</td>\n",
       "      <td>346.084324</td>\n",
       "      <td>428.659085</td>\n",
       "      <td>...</td>\n",
       "      <td>0.409450</td>\n",
       "      <td>3.577424</td>\n",
       "      <td>0.671235</td>\n",
       "      <td>0.693807</td>\n",
       "      <td>0.402342</td>\n",
       "      <td>3.590282</td>\n",
       "      <td>8.865859</td>\n",
       "      <td>0.678016</td>\n",
       "      <td>0.409450</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>77</th>\n",
       "      <td>92.208480</td>\n",
       "      <td>116.958007</td>\n",
       "      <td>125.609129</td>\n",
       "      <td>137.531339</td>\n",
       "      <td>257.522973</td>\n",
       "      <td>336.428016</td>\n",
       "      <td>145.278474</td>\n",
       "      <td>144.172997</td>\n",
       "      <td>201.550187</td>\n",
       "      <td>189.609357</td>\n",
       "      <td>...</td>\n",
       "      <td>0.339736</td>\n",
       "      <td>3.042537</td>\n",
       "      <td>0.669043</td>\n",
       "      <td>0.702060</td>\n",
       "      <td>0.346723</td>\n",
       "      <td>3.644536</td>\n",
       "      <td>7.998100</td>\n",
       "      <td>0.850882</td>\n",
       "      <td>0.339736</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>78</th>\n",
       "      <td>133.263194</td>\n",
       "      <td>151.053398</td>\n",
       "      <td>318.424137</td>\n",
       "      <td>348.731678</td>\n",
       "      <td>447.613387</td>\n",
       "      <td>354.625726</td>\n",
       "      <td>286.582090</td>\n",
       "      <td>516.643654</td>\n",
       "      <td>576.053540</td>\n",
       "      <td>446.460535</td>\n",
       "      <td>...</td>\n",
       "      <td>0.366773</td>\n",
       "      <td>2.765730</td>\n",
       "      <td>0.729007</td>\n",
       "      <td>0.750096</td>\n",
       "      <td>0.443597</td>\n",
       "      <td>3.813667</td>\n",
       "      <td>7.997914</td>\n",
       "      <td>0.785152</td>\n",
       "      <td>0.366773</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79</th>\n",
       "      <td>249.617943</td>\n",
       "      <td>216.492815</td>\n",
       "      <td>611.193861</td>\n",
       "      <td>353.540734</td>\n",
       "      <td>411.746527</td>\n",
       "      <td>581.867881</td>\n",
       "      <td>268.298119</td>\n",
       "      <td>545.626529</td>\n",
       "      <td>262.406775</td>\n",
       "      <td>423.683028</td>\n",
       "      <td>...</td>\n",
       "      <td>0.353116</td>\n",
       "      <td>3.197137</td>\n",
       "      <td>0.951716</td>\n",
       "      <td>0.671189</td>\n",
       "      <td>0.388950</td>\n",
       "      <td>3.456160</td>\n",
       "      <td>8.127271</td>\n",
       "      <td>0.813018</td>\n",
       "      <td>0.353116</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>80</th>\n",
       "      <td>287.285656</td>\n",
       "      <td>114.657563</td>\n",
       "      <td>238.444302</td>\n",
       "      <td>218.676056</td>\n",
       "      <td>528.903014</td>\n",
       "      <td>499.132091</td>\n",
       "      <td>218.285677</td>\n",
       "      <td>383.158155</td>\n",
       "      <td>433.604277</td>\n",
       "      <td>822.036696</td>\n",
       "      <td>...</td>\n",
       "      <td>0.342787</td>\n",
       "      <td>3.912154</td>\n",
       "      <td>0.955616</td>\n",
       "      <td>0.729682</td>\n",
       "      <td>0.402445</td>\n",
       "      <td>3.712143</td>\n",
       "      <td>7.867940</td>\n",
       "      <td>0.743499</td>\n",
       "      <td>0.342787</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>81</th>\n",
       "      <td>235.658934</td>\n",
       "      <td>228.465556</td>\n",
       "      <td>165.119692</td>\n",
       "      <td>164.523954</td>\n",
       "      <td>451.406382</td>\n",
       "      <td>721.150827</td>\n",
       "      <td>170.684493</td>\n",
       "      <td>271.820094</td>\n",
       "      <td>497.363975</td>\n",
       "      <td>1026.909076</td>\n",
       "      <td>...</td>\n",
       "      <td>0.424034</td>\n",
       "      <td>3.421094</td>\n",
       "      <td>0.589036</td>\n",
       "      <td>0.778009</td>\n",
       "      <td>0.367237</td>\n",
       "      <td>3.714159</td>\n",
       "      <td>7.961158</td>\n",
       "      <td>0.768835</td>\n",
       "      <td>0.424034</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82</th>\n",
       "      <td>1361.341189</td>\n",
       "      <td>560.704333</td>\n",
       "      <td>377.013651</td>\n",
       "      <td>302.907704</td>\n",
       "      <td>965.564365</td>\n",
       "      <td>1069.805076</td>\n",
       "      <td>574.720248</td>\n",
       "      <td>628.673307</td>\n",
       "      <td>642.301453</td>\n",
       "      <td>1005.881923</td>\n",
       "      <td>...</td>\n",
       "      <td>0.377603</td>\n",
       "      <td>3.103968</td>\n",
       "      <td>0.577075</td>\n",
       "      <td>0.789814</td>\n",
       "      <td>0.446599</td>\n",
       "      <td>3.761996</td>\n",
       "      <td>7.830576</td>\n",
       "      <td>0.748145</td>\n",
       "      <td>0.377603</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>83</th>\n",
       "      <td>129.011223</td>\n",
       "      <td>146.677036</td>\n",
       "      <td>173.006724</td>\n",
       "      <td>165.247269</td>\n",
       "      <td>226.673642</td>\n",
       "      <td>313.343127</td>\n",
       "      <td>138.280991</td>\n",
       "      <td>243.445029</td>\n",
       "      <td>163.412541</td>\n",
       "      <td>350.368765</td>\n",
       "      <td>...</td>\n",
       "      <td>0.340245</td>\n",
       "      <td>3.619466</td>\n",
       "      <td>0.588140</td>\n",
       "      <td>0.572551</td>\n",
       "      <td>0.380310</td>\n",
       "      <td>3.682308</td>\n",
       "      <td>7.828496</td>\n",
       "      <td>0.680661</td>\n",
       "      <td>0.340245</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>84</th>\n",
       "      <td>304.801163</td>\n",
       "      <td>220.782539</td>\n",
       "      <td>237.866398</td>\n",
       "      <td>220.630839</td>\n",
       "      <td>521.646422</td>\n",
       "      <td>406.008103</td>\n",
       "      <td>166.187554</td>\n",
       "      <td>269.904440</td>\n",
       "      <td>485.759566</td>\n",
       "      <td>388.666356</td>\n",
       "      <td>...</td>\n",
       "      <td>0.415399</td>\n",
       "      <td>3.355364</td>\n",
       "      <td>0.604846</td>\n",
       "      <td>0.650073</td>\n",
       "      <td>0.373266</td>\n",
       "      <td>4.147951</td>\n",
       "      <td>8.147062</td>\n",
       "      <td>0.781019</td>\n",
       "      <td>0.415399</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>85</th>\n",
       "      <td>711.107350</td>\n",
       "      <td>223.594562</td>\n",
       "      <td>276.553212</td>\n",
       "      <td>219.691153</td>\n",
       "      <td>656.107442</td>\n",
       "      <td>393.870792</td>\n",
       "      <td>337.285762</td>\n",
       "      <td>382.956131</td>\n",
       "      <td>650.180586</td>\n",
       "      <td>487.718554</td>\n",
       "      <td>...</td>\n",
       "      <td>0.372288</td>\n",
       "      <td>3.301264</td>\n",
       "      <td>0.789874</td>\n",
       "      <td>0.829736</td>\n",
       "      <td>0.383227</td>\n",
       "      <td>3.529338</td>\n",
       "      <td>8.029581</td>\n",
       "      <td>0.725724</td>\n",
       "      <td>0.372288</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86</th>\n",
       "      <td>135.675406</td>\n",
       "      <td>109.834773</td>\n",
       "      <td>164.740253</td>\n",
       "      <td>206.855546</td>\n",
       "      <td>223.931121</td>\n",
       "      <td>377.184648</td>\n",
       "      <td>160.657789</td>\n",
       "      <td>262.759996</td>\n",
       "      <td>225.824638</td>\n",
       "      <td>272.525883</td>\n",
       "      <td>...</td>\n",
       "      <td>0.326184</td>\n",
       "      <td>3.751195</td>\n",
       "      <td>0.660056</td>\n",
       "      <td>0.752874</td>\n",
       "      <td>0.487828</td>\n",
       "      <td>4.082743</td>\n",
       "      <td>8.174097</td>\n",
       "      <td>0.702796</td>\n",
       "      <td>0.326184</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>87</th>\n",
       "      <td>161.980780</td>\n",
       "      <td>155.727929</td>\n",
       "      <td>151.253113</td>\n",
       "      <td>147.473257</td>\n",
       "      <td>308.917323</td>\n",
       "      <td>588.494566</td>\n",
       "      <td>146.279572</td>\n",
       "      <td>264.623908</td>\n",
       "      <td>469.286730</td>\n",
       "      <td>653.984137</td>\n",
       "      <td>...</td>\n",
       "      <td>0.472736</td>\n",
       "      <td>3.661577</td>\n",
       "      <td>0.681873</td>\n",
       "      <td>1.012933</td>\n",
       "      <td>0.367742</td>\n",
       "      <td>3.858784</td>\n",
       "      <td>8.166584</td>\n",
       "      <td>0.809441</td>\n",
       "      <td>0.472736</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>88</th>\n",
       "      <td>277.188918</td>\n",
       "      <td>273.849546</td>\n",
       "      <td>245.642152</td>\n",
       "      <td>214.986933</td>\n",
       "      <td>506.417361</td>\n",
       "      <td>998.192086</td>\n",
       "      <td>217.223270</td>\n",
       "      <td>304.823340</td>\n",
       "      <td>531.644168</td>\n",
       "      <td>737.109612</td>\n",
       "      <td>...</td>\n",
       "      <td>0.492119</td>\n",
       "      <td>2.933393</td>\n",
       "      <td>0.932485</td>\n",
       "      <td>0.958366</td>\n",
       "      <td>0.409567</td>\n",
       "      <td>3.819302</td>\n",
       "      <td>8.139295</td>\n",
       "      <td>0.822220</td>\n",
       "      <td>0.492119</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>289.601539</td>\n",
       "      <td>282.423360</td>\n",
       "      <td>158.986420</td>\n",
       "      <td>133.064344</td>\n",
       "      <td>458.467570</td>\n",
       "      <td>811.732916</td>\n",
       "      <td>256.585002</td>\n",
       "      <td>204.001035</td>\n",
       "      <td>708.910683</td>\n",
       "      <td>718.531933</td>\n",
       "      <td>...</td>\n",
       "      <td>0.382146</td>\n",
       "      <td>3.817297</td>\n",
       "      <td>0.713851</td>\n",
       "      <td>0.580734</td>\n",
       "      <td>0.385617</td>\n",
       "      <td>3.551648</td>\n",
       "      <td>8.646808</td>\n",
       "      <td>0.680684</td>\n",
       "      <td>0.382146</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>90</th>\n",
       "      <td>234.369949</td>\n",
       "      <td>173.429167</td>\n",
       "      <td>184.428196</td>\n",
       "      <td>177.303277</td>\n",
       "      <td>468.564661</td>\n",
       "      <td>496.952122</td>\n",
       "      <td>234.807989</td>\n",
       "      <td>265.488176</td>\n",
       "      <td>253.531772</td>\n",
       "      <td>289.464300</td>\n",
       "      <td>...</td>\n",
       "      <td>0.409852</td>\n",
       "      <td>3.895180</td>\n",
       "      <td>0.515565</td>\n",
       "      <td>0.638427</td>\n",
       "      <td>0.350565</td>\n",
       "      <td>4.311694</td>\n",
       "      <td>8.528738</td>\n",
       "      <td>0.811338</td>\n",
       "      <td>0.409852</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>91</th>\n",
       "      <td>178.131305</td>\n",
       "      <td>176.227246</td>\n",
       "      <td>202.702128</td>\n",
       "      <td>133.476574</td>\n",
       "      <td>352.496692</td>\n",
       "      <td>560.309340</td>\n",
       "      <td>175.676208</td>\n",
       "      <td>196.443775</td>\n",
       "      <td>287.358088</td>\n",
       "      <td>419.611402</td>\n",
       "      <td>...</td>\n",
       "      <td>0.452356</td>\n",
       "      <td>3.177358</td>\n",
       "      <td>0.676688</td>\n",
       "      <td>0.684944</td>\n",
       "      <td>0.415825</td>\n",
       "      <td>3.693091</td>\n",
       "      <td>8.008673</td>\n",
       "      <td>0.731318</td>\n",
       "      <td>0.452356</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92</th>\n",
       "      <td>304.607781</td>\n",
       "      <td>371.436987</td>\n",
       "      <td>266.015597</td>\n",
       "      <td>160.723638</td>\n",
       "      <td>457.374406</td>\n",
       "      <td>632.608562</td>\n",
       "      <td>273.757381</td>\n",
       "      <td>268.855794</td>\n",
       "      <td>374.833377</td>\n",
       "      <td>610.540396</td>\n",
       "      <td>...</td>\n",
       "      <td>0.415381</td>\n",
       "      <td>3.622058</td>\n",
       "      <td>0.772391</td>\n",
       "      <td>0.779936</td>\n",
       "      <td>0.372088</td>\n",
       "      <td>3.766363</td>\n",
       "      <td>8.377726</td>\n",
       "      <td>0.737688</td>\n",
       "      <td>0.415381</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>93</th>\n",
       "      <td>324.336720</td>\n",
       "      <td>193.037482</td>\n",
       "      <td>149.427305</td>\n",
       "      <td>79.716942</td>\n",
       "      <td>671.158529</td>\n",
       "      <td>436.239843</td>\n",
       "      <td>152.367038</td>\n",
       "      <td>185.765710</td>\n",
       "      <td>413.048123</td>\n",
       "      <td>415.002397</td>\n",
       "      <td>...</td>\n",
       "      <td>0.315938</td>\n",
       "      <td>3.460561</td>\n",
       "      <td>0.923417</td>\n",
       "      <td>0.813748</td>\n",
       "      <td>0.369319</td>\n",
       "      <td>3.980435</td>\n",
       "      <td>7.931937</td>\n",
       "      <td>0.736686</td>\n",
       "      <td>0.315938</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>94</th>\n",
       "      <td>319.877060</td>\n",
       "      <td>248.020385</td>\n",
       "      <td>176.858312</td>\n",
       "      <td>112.162292</td>\n",
       "      <td>420.255576</td>\n",
       "      <td>462.475224</td>\n",
       "      <td>182.493990</td>\n",
       "      <td>213.649302</td>\n",
       "      <td>395.688687</td>\n",
       "      <td>444.165837</td>\n",
       "      <td>...</td>\n",
       "      <td>0.350740</td>\n",
       "      <td>3.708630</td>\n",
       "      <td>0.682530</td>\n",
       "      <td>0.845809</td>\n",
       "      <td>0.408128</td>\n",
       "      <td>3.987005</td>\n",
       "      <td>7.685725</td>\n",
       "      <td>0.757725</td>\n",
       "      <td>0.350740</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95</th>\n",
       "      <td>370.749741</td>\n",
       "      <td>272.409768</td>\n",
       "      <td>166.026281</td>\n",
       "      <td>128.237646</td>\n",
       "      <td>459.416067</td>\n",
       "      <td>783.189814</td>\n",
       "      <td>197.098012</td>\n",
       "      <td>354.068130</td>\n",
       "      <td>446.461677</td>\n",
       "      <td>648.490928</td>\n",
       "      <td>...</td>\n",
       "      <td>0.340319</td>\n",
       "      <td>3.261493</td>\n",
       "      <td>0.650505</td>\n",
       "      <td>0.611201</td>\n",
       "      <td>0.395495</td>\n",
       "      <td>3.773324</td>\n",
       "      <td>7.898320</td>\n",
       "      <td>0.705151</td>\n",
       "      <td>0.340319</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>96</th>\n",
       "      <td>623.755544</td>\n",
       "      <td>326.300221</td>\n",
       "      <td>153.184294</td>\n",
       "      <td>154.258128</td>\n",
       "      <td>534.045163</td>\n",
       "      <td>1327.872210</td>\n",
       "      <td>183.077307</td>\n",
       "      <td>233.344511</td>\n",
       "      <td>534.873844</td>\n",
       "      <td>679.320297</td>\n",
       "      <td>...</td>\n",
       "      <td>0.349802</td>\n",
       "      <td>3.498607</td>\n",
       "      <td>0.555803</td>\n",
       "      <td>0.602244</td>\n",
       "      <td>0.369587</td>\n",
       "      <td>3.519048</td>\n",
       "      <td>8.399984</td>\n",
       "      <td>0.708840</td>\n",
       "      <td>0.349802</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>97</th>\n",
       "      <td>252.039309</td>\n",
       "      <td>236.559220</td>\n",
       "      <td>197.874672</td>\n",
       "      <td>155.375508</td>\n",
       "      <td>560.624437</td>\n",
       "      <td>497.415952</td>\n",
       "      <td>275.978856</td>\n",
       "      <td>311.177271</td>\n",
       "      <td>487.525446</td>\n",
       "      <td>527.299013</td>\n",
       "      <td>...</td>\n",
       "      <td>0.365616</td>\n",
       "      <td>4.161234</td>\n",
       "      <td>0.666866</td>\n",
       "      <td>0.646731</td>\n",
       "      <td>0.396200</td>\n",
       "      <td>4.125439</td>\n",
       "      <td>7.826471</td>\n",
       "      <td>0.772273</td>\n",
       "      <td>0.365616</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98</th>\n",
       "      <td>534.123192</td>\n",
       "      <td>277.817120</td>\n",
       "      <td>265.240364</td>\n",
       "      <td>176.544917</td>\n",
       "      <td>802.175504</td>\n",
       "      <td>1186.208958</td>\n",
       "      <td>290.126728</td>\n",
       "      <td>427.346321</td>\n",
       "      <td>675.797607</td>\n",
       "      <td>1055.705640</td>\n",
       "      <td>...</td>\n",
       "      <td>0.458701</td>\n",
       "      <td>3.763322</td>\n",
       "      <td>0.642405</td>\n",
       "      <td>0.730825</td>\n",
       "      <td>0.359031</td>\n",
       "      <td>3.895707</td>\n",
       "      <td>8.423729</td>\n",
       "      <td>0.791314</td>\n",
       "      <td>0.458701</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99</th>\n",
       "      <td>496.699587</td>\n",
       "      <td>290.759369</td>\n",
       "      <td>258.886416</td>\n",
       "      <td>186.935119</td>\n",
       "      <td>434.017592</td>\n",
       "      <td>909.756813</td>\n",
       "      <td>315.683917</td>\n",
       "      <td>454.139586</td>\n",
       "      <td>546.748199</td>\n",
       "      <td>724.887901</td>\n",
       "      <td>...</td>\n",
       "      <td>0.435276</td>\n",
       "      <td>3.306024</td>\n",
       "      <td>0.670142</td>\n",
       "      <td>0.658668</td>\n",
       "      <td>0.439468</td>\n",
       "      <td>3.543439</td>\n",
       "      <td>8.295603</td>\n",
       "      <td>0.845966</td>\n",
       "      <td>0.435276</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>100 rows × 116 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "            0           1           2           3            4            5    \\\n",
       "0    178.786422  187.299138  373.456785  147.897919   433.100880   439.940230   \n",
       "1    246.545195  146.974897  407.334715  258.835636   733.011971   259.251316   \n",
       "2    154.673070  152.911484  346.467292  296.324842   448.180389   407.671629   \n",
       "3    142.616877  123.756485  274.717083  262.639245   786.032023   453.893509   \n",
       "4    260.544638  196.312737  576.709094  318.495857   711.977651   452.092105   \n",
       "5    234.320907  261.550282  470.721632  286.447461   751.604378   469.626384   \n",
       "6    273.379424  259.283579  521.113884  279.343219   571.597933   461.339238   \n",
       "7    312.157337  327.162636  605.796649  393.135495  1041.770233   647.134773   \n",
       "8     90.867944  183.301948  179.605055  206.120312   340.883022   351.363763   \n",
       "9    121.823402  143.578508  188.337742  315.325243   290.538303   463.355738   \n",
       "10   201.869968  207.343892  358.151280  313.966200   575.218551   700.470319   \n",
       "11   127.500698  140.006731  311.099870  335.134552   385.807141   420.605587   \n",
       "12   117.045834  115.357530  134.405483  234.935440   235.938726   237.162789   \n",
       "13   195.575947  199.071527  274.309503  282.684451   524.998364   451.207193   \n",
       "14   195.393033  270.996408  296.393823  357.205989   349.385642   385.753793   \n",
       "15   176.691208  177.027086  398.185835  223.864522   517.424334   362.294943   \n",
       "16   280.072969  163.121439  317.006730  188.162560   912.977836   375.896532   \n",
       "17   173.215078  149.073981  150.037932  152.980695   568.715017   366.064307   \n",
       "18   172.190633  202.396157  299.976177  160.355398   689.230772   564.925225   \n",
       "19   253.430027  268.879732  390.696573  340.222539   777.303251   626.193163   \n",
       "20   299.869998  227.688413  358.779289  151.046274   992.183611   517.627875   \n",
       "21   176.532272  106.073423  295.705094  246.425916   381.403013   274.749131   \n",
       "22   228.204807  316.282246  341.836231  280.807598   711.831743   758.850904   \n",
       "23   166.710607  202.597143  392.102092  274.927670   743.108185   633.448967   \n",
       "24   313.305186  188.154894  390.240316  262.611060   610.824073   672.614760   \n",
       "25   192.454194  162.745680  200.520404  291.521940   585.399265   354.628730   \n",
       "26   289.190744  299.188237  301.265454  262.134823   692.545951   669.218277   \n",
       "27   262.605588  244.991274  667.564831  408.848913   473.418976   532.245267   \n",
       "28   154.738055  222.056288  389.591496  290.912424   472.113157   501.991482   \n",
       "29   147.971812  155.431087  334.884672  400.147118   203.585776   277.576426   \n",
       "..          ...         ...         ...         ...          ...          ...   \n",
       "70   261.134719  258.302304  507.139181  350.101826   507.061498   665.581335   \n",
       "71   230.820184  227.389183  388.022951  117.647582   291.595835   223.590412   \n",
       "72   199.441131  216.596367  500.309393  128.812704   337.809754   215.469931   \n",
       "73   138.730279  191.174348  212.830176  166.411185   397.164134   298.816490   \n",
       "74   124.189176  166.689290  311.292913  179.981363   314.544481   449.685500   \n",
       "75   177.371730  216.114620  465.725745  294.850166   271.478315   244.936852   \n",
       "76   212.262949  256.926485  290.546175  246.348148   481.660190   569.281736   \n",
       "77    92.208480  116.958007  125.609129  137.531339   257.522973   336.428016   \n",
       "78   133.263194  151.053398  318.424137  348.731678   447.613387   354.625726   \n",
       "79   249.617943  216.492815  611.193861  353.540734   411.746527   581.867881   \n",
       "80   287.285656  114.657563  238.444302  218.676056   528.903014   499.132091   \n",
       "81   235.658934  228.465556  165.119692  164.523954   451.406382   721.150827   \n",
       "82  1361.341189  560.704333  377.013651  302.907704   965.564365  1069.805076   \n",
       "83   129.011223  146.677036  173.006724  165.247269   226.673642   313.343127   \n",
       "84   304.801163  220.782539  237.866398  220.630839   521.646422   406.008103   \n",
       "85   711.107350  223.594562  276.553212  219.691153   656.107442   393.870792   \n",
       "86   135.675406  109.834773  164.740253  206.855546   223.931121   377.184648   \n",
       "87   161.980780  155.727929  151.253113  147.473257   308.917323   588.494566   \n",
       "88   277.188918  273.849546  245.642152  214.986933   506.417361   998.192086   \n",
       "89   289.601539  282.423360  158.986420  133.064344   458.467570   811.732916   \n",
       "90   234.369949  173.429167  184.428196  177.303277   468.564661   496.952122   \n",
       "91   178.131305  176.227246  202.702128  133.476574   352.496692   560.309340   \n",
       "92   304.607781  371.436987  266.015597  160.723638   457.374406   632.608562   \n",
       "93   324.336720  193.037482  149.427305   79.716942   671.158529   436.239843   \n",
       "94   319.877060  248.020385  176.858312  112.162292   420.255576   462.475224   \n",
       "95   370.749741  272.409768  166.026281  128.237646   459.416067   783.189814   \n",
       "96   623.755544  326.300221  153.184294  154.258128   534.045163  1327.872210   \n",
       "97   252.039309  236.559220  197.874672  155.375508   560.624437   497.415952   \n",
       "98   534.123192  277.817120  265.240364  176.544917   802.175504  1186.208958   \n",
       "99   496.699587  290.759369  258.886416  186.935119   434.017592   909.756813   \n",
       "\n",
       "           6           7            8            9    ...       106       107  \\\n",
       "0   333.774992  416.115508   394.691618   368.760198  ...  0.410709  3.580521   \n",
       "1   281.086734  447.975677   726.113512   372.659878  ...  0.408117  3.408203   \n",
       "2   272.390824  534.975166   525.624321   372.344138  ...  0.321299  3.311277   \n",
       "3   325.054145  585.962765   691.101114   508.397870  ...  0.395649  3.612864   \n",
       "4   368.535053  774.207210   544.819770   521.200692  ...  0.350373  3.773968   \n",
       "5   357.002351  577.517453   785.321363   489.468379  ...  0.346077  3.542306   \n",
       "6   522.689588  474.163002   568.110176   406.117779  ...  0.352900  3.593283   \n",
       "7   511.948216  713.981772  1216.144852   716.313123  ...  0.339203  3.759761   \n",
       "8   279.098579  504.991956   267.143485   271.414462  ...  0.363857  3.533879   \n",
       "9   185.048609  490.961345   216.813588   398.905451  ...  0.355697  3.714707   \n",
       "10  288.476645  446.116792   485.959500   569.900688  ...  0.352160  3.532321   \n",
       "11  234.936613  432.448664   297.271949   323.821490  ...  0.397548  3.479982   \n",
       "12  148.821140  320.215899   181.596172   181.022392  ...  0.366877  3.452338   \n",
       "13  293.616764  448.114314   343.557154   355.621789  ...  0.312612  3.511218   \n",
       "14  249.576777  490.626714   349.490917   266.248300  ...  0.405448  3.814354   \n",
       "15  281.813067  438.809774   493.636955   436.563591  ...  0.345647  3.607638   \n",
       "16  256.447720  443.113740   524.124246   264.714965  ...  0.357802  3.766675   \n",
       "17  247.497097  331.551166   472.490192   244.326979  ...  0.422197  3.184529   \n",
       "18  340.747540  437.707263   682.712240   824.878818  ...  0.364102  3.716139   \n",
       "19  366.015566  605.850189   588.125692   589.722847  ...  0.452971  3.702314   \n",
       "20  364.925738  474.377177  1047.339851   607.868023  ...  0.419587  3.258614   \n",
       "21  276.640917  336.243612   305.400304   152.363553  ...  0.389708  3.440791   \n",
       "22  410.210644  612.924933   398.635739   585.706312  ...  0.347397  3.557350   \n",
       "23  392.219767  812.956043   591.870036   472.259010  ...  0.442770  3.740394   \n",
       "24  459.503370  559.197594   624.338175   533.952562  ...  0.466055  3.497552   \n",
       "25  312.605028  531.314156   524.096684   272.855330  ...  0.415696  3.346305   \n",
       "26  438.199987  714.247490   667.919096   548.131484  ...  0.340344  3.520040   \n",
       "27  331.669166  479.071785   445.218419   397.654977  ...  0.289427  3.392331   \n",
       "28  424.592671  657.931235   453.669150   444.013337  ...  0.409222  3.558121   \n",
       "29  315.079381  417.870902   209.603915   273.144925  ...  0.395049  3.560522   \n",
       "..         ...         ...          ...          ...  ...       ...       ...   \n",
       "70  251.476918  555.539360   430.354020   432.256490  ...  0.425885  3.683649   \n",
       "71  135.671843  211.902017   245.123641   269.888973  ...  0.372161  3.775336   \n",
       "72  163.186875  249.513322   221.490544   330.834000  ...  0.377273  3.630585   \n",
       "73  261.386305  284.328117   351.088574   288.612735  ...  0.417536  3.335675   \n",
       "74  259.772225  359.374297   248.261614   543.192428  ...  0.411413  3.793024   \n",
       "75  146.350015  419.355925   115.308624   283.393541  ...  0.420551  3.995588   \n",
       "76  239.155450  348.768704   346.084324   428.659085  ...  0.409450  3.577424   \n",
       "77  145.278474  144.172997   201.550187   189.609357  ...  0.339736  3.042537   \n",
       "78  286.582090  516.643654   576.053540   446.460535  ...  0.366773  2.765730   \n",
       "79  268.298119  545.626529   262.406775   423.683028  ...  0.353116  3.197137   \n",
       "80  218.285677  383.158155   433.604277   822.036696  ...  0.342787  3.912154   \n",
       "81  170.684493  271.820094   497.363975  1026.909076  ...  0.424034  3.421094   \n",
       "82  574.720248  628.673307   642.301453  1005.881923  ...  0.377603  3.103968   \n",
       "83  138.280991  243.445029   163.412541   350.368765  ...  0.340245  3.619466   \n",
       "84  166.187554  269.904440   485.759566   388.666356  ...  0.415399  3.355364   \n",
       "85  337.285762  382.956131   650.180586   487.718554  ...  0.372288  3.301264   \n",
       "86  160.657789  262.759996   225.824638   272.525883  ...  0.326184  3.751195   \n",
       "87  146.279572  264.623908   469.286730   653.984137  ...  0.472736  3.661577   \n",
       "88  217.223270  304.823340   531.644168   737.109612  ...  0.492119  2.933393   \n",
       "89  256.585002  204.001035   708.910683   718.531933  ...  0.382146  3.817297   \n",
       "90  234.807989  265.488176   253.531772   289.464300  ...  0.409852  3.895180   \n",
       "91  175.676208  196.443775   287.358088   419.611402  ...  0.452356  3.177358   \n",
       "92  273.757381  268.855794   374.833377   610.540396  ...  0.415381  3.622058   \n",
       "93  152.367038  185.765710   413.048123   415.002397  ...  0.315938  3.460561   \n",
       "94  182.493990  213.649302   395.688687   444.165837  ...  0.350740  3.708630   \n",
       "95  197.098012  354.068130   446.461677   648.490928  ...  0.340319  3.261493   \n",
       "96  183.077307  233.344511   534.873844   679.320297  ...  0.349802  3.498607   \n",
       "97  275.978856  311.177271   487.525446   527.299013  ...  0.365616  4.161234   \n",
       "98  290.126728  427.346321   675.797607  1055.705640  ...  0.458701  3.763322   \n",
       "99  315.683917  454.139586   546.748199   724.887901  ...  0.435276  3.306024   \n",
       "\n",
       "         108       109       110       111       112       113       114  115  \n",
       "0   0.574664  0.676765  0.389804  3.855179  8.851754  0.729672  0.410709  0.0  \n",
       "1   0.704777  0.795886  0.492413  4.095119  8.313062  0.824171  0.408117  0.0  \n",
       "2   0.696812  0.748675  0.366602  3.904908  8.709865  0.767086  0.321299  0.0  \n",
       "3   0.795568  0.716539  0.344873  3.920525  8.164512  0.807174  0.395649  0.0  \n",
       "4   0.753706  0.902003  0.380239  4.014634  7.777630  0.833583  0.350373  0.0  \n",
       "5   0.545660  0.555900  0.331486  3.823273  8.410936  0.756244  0.346077  0.0  \n",
       "6   0.635650  0.667032  0.392782  3.953978  8.084922  0.754650  0.352900  0.0  \n",
       "7   0.659787  0.565934  0.356684  3.929324  8.775190  0.805310  0.339203  0.0  \n",
       "8   0.541001  0.655530  0.331713  3.923096  8.709387  0.794712  0.363857  0.0  \n",
       "9   0.488094  0.660672  0.367932  4.054801  8.720214  0.790561  0.355697  0.0  \n",
       "10  0.597552  0.723606  0.449655  4.010461  8.305585  0.858595  0.352160  0.0  \n",
       "11  0.665714  0.629633  0.393463  3.792875  7.799580  0.731764  0.397548  0.0  \n",
       "12  0.547187  0.605098  0.378353  3.587219  7.997030  0.695318  0.366877  0.0  \n",
       "13  0.571026  0.620062  0.364408  3.688348  8.361690  0.753523  0.312612  0.0  \n",
       "14  0.570337  0.706242  0.392385  3.601962  8.723980  0.745340  0.405448  0.0  \n",
       "15  0.752208  0.576956  0.341811  3.690010  8.485560  0.796332  0.345647  0.0  \n",
       "16  0.687498  0.711523  0.374456  3.939251  7.935518  0.798783  0.357802  0.0  \n",
       "17  0.604299  0.719281  0.390349  3.709922  7.721513  0.792122  0.422197  0.0  \n",
       "18  0.672627  0.819458  0.404210  3.872358  8.345090  0.848888  0.364102  0.0  \n",
       "19  0.557211  0.682076  0.457097  3.775675  8.739288  0.772543  0.452971  0.0  \n",
       "20  0.618107  0.667598  0.411717  3.785812  8.555000  0.867401  0.419587  0.0  \n",
       "21  0.593434  0.684052  0.388469  3.629259  8.198102  0.838601  0.389708  0.0  \n",
       "22  0.892613  0.688156  0.383629  4.118524  8.289979  0.772936  0.347397  0.0  \n",
       "23  0.653449  0.603795  0.363151  3.721382  8.873737  0.761885  0.442770  0.0  \n",
       "24  0.594804  0.740771  0.461666  3.822586  8.592550  0.775473  0.466055  0.0  \n",
       "25  0.557133  0.573472  0.416878  3.943003  8.188199  0.833220  0.415696  0.0  \n",
       "26  0.690583  0.593332  0.421714  3.918895  8.064142  0.741595  0.340344  0.0  \n",
       "27  0.601149  0.613482  0.391231  3.736227  7.861983  0.797338  0.289427  0.0  \n",
       "28  0.691285  0.729957  0.362660  3.640856  8.568836  0.791814  0.409222  0.0  \n",
       "29  0.789697  0.962958  0.445118  3.925553  8.646597  0.771469  0.395049  0.0  \n",
       "..       ...       ...       ...       ...       ...       ...       ...  ...  \n",
       "70  0.957510  0.740335  0.372962  4.348339  8.573490  0.876375  0.425885  0.0  \n",
       "71  0.728939  0.784113  0.471682  3.614408  7.997206  0.667619  0.372161  0.0  \n",
       "72  0.533531  0.768592  0.326154  3.480983  8.328907  0.758970  0.377273  0.0  \n",
       "73  0.627466  0.738017  0.383787  3.360009  8.394821  0.674719  0.417536  0.0  \n",
       "74  0.489882  0.649907  0.437163  3.789730  8.660178  0.737800  0.411413  0.0  \n",
       "75  0.574573  0.761602  0.482047  4.297029  8.685668  0.904984  0.420551  0.0  \n",
       "76  0.671235  0.693807  0.402342  3.590282  8.865859  0.678016  0.409450  0.0  \n",
       "77  0.669043  0.702060  0.346723  3.644536  7.998100  0.850882  0.339736  0.0  \n",
       "78  0.729007  0.750096  0.443597  3.813667  7.997914  0.785152  0.366773  0.0  \n",
       "79  0.951716  0.671189  0.388950  3.456160  8.127271  0.813018  0.353116  0.0  \n",
       "80  0.955616  0.729682  0.402445  3.712143  7.867940  0.743499  0.342787  0.0  \n",
       "81  0.589036  0.778009  0.367237  3.714159  7.961158  0.768835  0.424034  0.0  \n",
       "82  0.577075  0.789814  0.446599  3.761996  7.830576  0.748145  0.377603  0.0  \n",
       "83  0.588140  0.572551  0.380310  3.682308  7.828496  0.680661  0.340245  0.0  \n",
       "84  0.604846  0.650073  0.373266  4.147951  8.147062  0.781019  0.415399  0.0  \n",
       "85  0.789874  0.829736  0.383227  3.529338  8.029581  0.725724  0.372288  0.0  \n",
       "86  0.660056  0.752874  0.487828  4.082743  8.174097  0.702796  0.326184  0.0  \n",
       "87  0.681873  1.012933  0.367742  3.858784  8.166584  0.809441  0.472736  0.0  \n",
       "88  0.932485  0.958366  0.409567  3.819302  8.139295  0.822220  0.492119  0.0  \n",
       "89  0.713851  0.580734  0.385617  3.551648  8.646808  0.680684  0.382146  0.0  \n",
       "90  0.515565  0.638427  0.350565  4.311694  8.528738  0.811338  0.409852  0.0  \n",
       "91  0.676688  0.684944  0.415825  3.693091  8.008673  0.731318  0.452356  0.0  \n",
       "92  0.772391  0.779936  0.372088  3.766363  8.377726  0.737688  0.415381  0.0  \n",
       "93  0.923417  0.813748  0.369319  3.980435  7.931937  0.736686  0.315938  0.0  \n",
       "94  0.682530  0.845809  0.408128  3.987005  7.685725  0.757725  0.350740  0.0  \n",
       "95  0.650505  0.611201  0.395495  3.773324  7.898320  0.705151  0.340319  0.0  \n",
       "96  0.555803  0.602244  0.369587  3.519048  8.399984  0.708840  0.349802  0.0  \n",
       "97  0.666866  0.646731  0.396200  4.125439  7.826471  0.772273  0.365616  0.0  \n",
       "98  0.642405  0.730825  0.359031  3.895707  8.423729  0.791314  0.458701  0.0  \n",
       "99  0.670142  0.658668  0.439468  3.543439  8.295603  0.845966  0.435276  0.0  \n",
       "\n",
       "[100 rows x 116 columns]"
      ]
     },
<<<<<<< HEAD
     "execution_count": 15,
=======
     "execution_count": 103,
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\n",
    "res = pd.DataFrame()\n",
    "\n",
    "#first iteration run in 'w' mode, all subsequent iteration run in 'a' mode\n",
    "#with open('test.csv', 'w') as file:\n",
    "#    writer = csv.writer(file)\n",
    "while temp < end:    \n",
    "    row = add_row(data, temp, temp + duration)\n",
    "    res=res.append(pd.Series(row),ignore_index=True)\n",
    "#        writer.writerow(row)\n",
    "    temp += duration\n",
    "\n",
    "#res.columns = ['Variance', 'Kurtosis', 'Skewness', 'Petrosian Fractal Dimension', 'Hjorth Mobility', 'Hjorth Complexity', 'Spectral Entropy', 'Label']\n",
    "res"
   ]
  },
  {
   "cell_type": "code",
<<<<<<< HEAD
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_final = pd.DataFrame()"
=======
   "execution_count": 85,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_final = pd.DataFrame()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 104,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_final = df_final.append(res)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 105,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(104, 116)"
      ]
     },
     "execution_count": 105,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_final.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 114,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.\n",
      " 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.\n",
      " 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.\n",
      " 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.\n",
      " 0. 0. 0. 0. 1. 1. 1. 1.]\n"
     ]
    }
   ],
   "source": [
    "X = df_final.iloc[:,0:115].values\n",
    "#print(X);\n",
    "y = df_final.iloc[:,-1].values\n",
    "print(y);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 115,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.ensemble import RandomForestClassifier\n",
    "from sklearn.model_selection import GridSearchCV"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 116,
   "metadata": {},
   "outputs": [],
   "source": [
    "rf_cs = RandomForestClassifier(n_jobs=-1,random_state=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 117,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\AAI\\anaconda3\\lib\\site-packages\\sklearn\\model_selection\\_split.py:670: UserWarning: The least populated class in y has only 4 members, which is less than n_splits=10.\n",
      "  warnings.warn((\"The least populated class in y has only %d\"\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'n_estimators': 50}\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "RandomForestClassifier(n_estimators=50, n_jobs=-1, random_state=0)"
      ]
     },
     "execution_count": 117,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "param_grid = { 'n_estimators': [10,50,100,200,250,500,700],}\n",
    "    \n",
    "CV_rfc = GridSearchCV(estimator=rf_cs, param_grid=param_grid, cv= 10)\n",
    "CV_rfc.fit(X, y)\n",
    "print (CV_rfc.best_params_)\n",
    "op_RF = CV_rfc.best_estimator_\n",
    "op_RF.fit(X,y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 118,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Requirement already satisfied: boruta in c:\\users\\aai\\anaconda3\\lib\\site-packages (0.3)\n",
      "Requirement already satisfied: scipy>=0.17.0 in c:\\users\\aai\\anaconda3\\lib\\site-packages (from boruta) (1.5.0)\n",
      "Requirement already satisfied: scikit-learn>=0.17.1 in c:\\users\\aai\\anaconda3\\lib\\site-packages (from boruta) (0.23.1)Note: you may need to restart the kernel to use updated packages.\n",
      "\n",
      "Requirement already satisfied: numpy>=1.10.4 in c:\\users\\aai\\anaconda3\\lib\\site-packages (from boruta) (1.18.5)\n",
      "Requirement already satisfied: joblib>=0.11 in c:\\users\\aai\\anaconda3\\lib\\site-packages (from scikit-learn>=0.17.1->boruta) (0.16.0)\n",
      "Requirement already satisfied: threadpoolctl>=2.0.0 in c:\\users\\aai\\anaconda3\\lib\\site-packages (from scikit-learn>=0.17.1->boruta) (2.1.0)\n"
     ]
    }
   ],
   "source": [
    "pip install boruta"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 119,
   "metadata": {},
   "outputs": [],
   "source": [
    "from boruta import BorutaPy\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 120,
   "metadata": {},
   "outputs": [],
   "source": [
    "# define Boruta feature selection method\n",
    "feat_selector = BorutaPy(op_RF, n_estimators='auto', verbose=2, random_state=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 121,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Iteration: \t1 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t115\n",
      "Rejected: \t0\n",
      "Iteration: \t2 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t115\n",
      "Rejected: \t0\n",
      "Iteration: \t3 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t115\n",
      "Rejected: \t0\n",
      "Iteration: \t4 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t115\n",
      "Rejected: \t0\n",
      "Iteration: \t5 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t115\n",
      "Rejected: \t0\n",
      "Iteration: \t6 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t115\n",
      "Rejected: \t0\n",
      "Iteration: \t7 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t115\n",
      "Rejected: \t0\n",
      "Iteration: \t8 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t90\n",
      "Rejected: \t25\n",
      "Iteration: \t9 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t90\n",
      "Rejected: \t25\n",
      "Iteration: \t10 / 100\n",
      "Confirmed: \t0\n",
      "Tentative: \t90\n",
      "Rejected: \t25\n",
      "Iteration: \t11 / 100\n",
      "Confirmed: \t1\n",
      "Tentative: \t89\n",
      "Rejected: \t25\n",
      "Iteration: \t12 / 100\n",
      "Confirmed: \t1\n",
      "Tentative: \t89\n",
      "Rejected: \t25\n",
      "Iteration: \t13 / 100\n",
      "Confirmed: \t1\n",
      "Tentative: \t89\n",
      "Rejected: \t25\n",
      "Iteration: \t14 / 100\n",
      "Confirmed: \t2\n",
      "Tentative: \t86\n",
      "Rejected: \t27\n",
      "Iteration: \t15 / 100\n",
      "Confirmed: \t2\n",
      "Tentative: \t86\n",
      "Rejected: \t27\n",
      "Iteration: \t16 / 100\n",
      "Confirmed: \t2\n",
      "Tentative: \t86\n",
      "Rejected: \t27\n",
      "Iteration: \t17 / 100\n",
      "Confirmed: \t2\n",
      "Tentative: \t83\n",
      "Rejected: \t30\n",
      "Iteration: \t18 / 100\n",
      "Confirmed: \t2\n",
      "Tentative: \t83\n",
      "Rejected: \t30\n",
      "Iteration: \t19 / 100\n",
      "Confirmed: \t7\n",
      "Tentative: \t78\n",
      "Rejected: \t30\n",
      "Iteration: \t20 / 100\n",
      "Confirmed: \t7\n",
      "Tentative: \t78\n",
      "Rejected: \t30\n",
      "Iteration: \t21 / 100\n",
      "Confirmed: \t7\n",
      "Tentative: \t78\n",
      "Rejected: \t30\n",
      "Iteration: \t22 / 100\n",
      "Confirmed: \t11\n",
      "Tentative: \t74\n",
      "Rejected: \t30\n",
      "Iteration: \t23 / 100\n",
      "Confirmed: \t11\n",
      "Tentative: \t74\n",
      "Rejected: \t30\n",
      "Iteration: \t24 / 100\n",
      "Confirmed: \t11\n",
      "Tentative: \t74\n",
      "Rejected: \t30\n",
      "Iteration: \t25 / 100\n",
      "Confirmed: \t11\n",
      "Tentative: \t74\n",
      "Rejected: \t30\n",
      "Iteration: \t26 / 100\n",
      "Confirmed: \t13\n",
      "Tentative: \t72\n",
      "Rejected: \t30\n",
      "Iteration: \t27 / 100\n",
      "Confirmed: \t13\n",
      "Tentative: \t72\n",
      "Rejected: \t30\n",
      "Iteration: \t28 / 100\n",
      "Confirmed: \t13\n",
      "Tentative: \t71\n",
      "Rejected: \t31\n",
      "Iteration: \t29 / 100\n",
      "Confirmed: \t16\n",
      "Tentative: \t66\n",
      "Rejected: \t33\n",
      "Iteration: \t30 / 100\n",
      "Confirmed: \t16\n",
      "Tentative: \t66\n",
      "Rejected: \t33\n",
      "Iteration: \t31 / 100\n",
      "Confirmed: \t16\n",
      "Tentative: \t66\n",
      "Rejected: \t33\n",
      "Iteration: \t32 / 100\n",
      "Confirmed: \t17\n",
      "Tentative: \t63\n",
      "Rejected: \t35\n",
      "Iteration: \t33 / 100\n",
      "Confirmed: \t17\n",
      "Tentative: \t63\n",
      "Rejected: \t35\n",
      "Iteration: \t34 / 100\n",
      "Confirmed: \t20\n",
      "Tentative: \t60\n",
      "Rejected: \t35\n",
      "Iteration: \t35 / 100\n",
      "Confirmed: \t20\n",
      "Tentative: \t60\n",
      "Rejected: \t35\n",
      "Iteration: \t36 / 100\n",
      "Confirmed: \t20\n",
      "Tentative: \t60\n",
      "Rejected: \t35\n",
      "Iteration: \t37 / 100\n",
      "Confirmed: \t22\n",
      "Tentative: \t58\n",
      "Rejected: \t35\n",
      "Iteration: \t38 / 100\n",
      "Confirmed: \t22\n",
      "Tentative: \t58\n",
      "Rejected: \t35\n",
      "Iteration: \t39 / 100\n",
      "Confirmed: \t22\n",
      "Tentative: \t58\n",
      "Rejected: \t35\n",
      "Iteration: \t40 / 100\n",
      "Confirmed: \t26\n",
      "Tentative: \t54\n",
      "Rejected: \t35\n",
      "Iteration: \t41 / 100\n",
      "Confirmed: \t26\n",
      "Tentative: \t54\n",
      "Rejected: \t35\n",
      "Iteration: \t42 / 100\n",
      "Confirmed: \t26\n",
      "Tentative: \t54\n",
      "Rejected: \t35\n",
      "Iteration: \t43 / 100\n",
      "Confirmed: \t27\n",
      "Tentative: \t53\n",
      "Rejected: \t35\n",
      "Iteration: \t44 / 100\n",
      "Confirmed: \t27\n",
      "Tentative: \t53\n",
      "Rejected: \t35\n",
      "Iteration: \t45 / 100\n",
      "Confirmed: \t27\n",
      "Tentative: \t53\n",
      "Rejected: \t35\n",
      "Iteration: \t46 / 100\n",
      "Confirmed: \t28\n",
      "Tentative: \t52\n",
      "Rejected: \t35\n",
      "Iteration: \t47 / 100\n",
      "Confirmed: \t28\n",
      "Tentative: \t52\n",
      "Rejected: \t35\n",
      "Iteration: \t48 / 100\n",
      "Confirmed: \t28\n",
      "Tentative: \t52\n",
      "Rejected: \t35\n",
      "Iteration: \t49 / 100\n",
      "Confirmed: \t30\n",
      "Tentative: \t50\n",
      "Rejected: \t35\n",
      "Iteration: \t50 / 100\n",
      "Confirmed: \t30\n",
      "Tentative: \t50\n",
      "Rejected: \t35\n",
      "Iteration: \t51 / 100\n",
      "Confirmed: \t31\n",
      "Tentative: \t49\n",
      "Rejected: \t35\n",
      "Iteration: \t52 / 100\n",
      "Confirmed: \t31\n",
      "Tentative: \t49\n",
      "Rejected: \t35\n",
      "Iteration: \t53 / 100\n",
      "Confirmed: \t31\n",
      "Tentative: \t49\n",
      "Rejected: \t35\n",
      "Iteration: \t54 / 100\n",
      "Confirmed: \t31\n",
      "Tentative: \t49\n",
      "Rejected: \t35\n",
      "Iteration: \t55 / 100\n",
      "Confirmed: \t31\n",
      "Tentative: \t49\n",
      "Rejected: \t35\n",
      "Iteration: \t56 / 100\n",
      "Confirmed: \t31\n",
      "Tentative: \t49\n",
      "Rejected: \t35\n",
      "Iteration: \t57 / 100\n",
      "Confirmed: \t31\n",
      "Tentative: \t49\n",
      "Rejected: \t35\n",
      "Iteration: \t58 / 100\n",
      "Confirmed: \t31\n",
      "Tentative: \t49\n",
      "Rejected: \t35\n",
      "Iteration: \t59 / 100\n",
      "Confirmed: \t32\n",
      "Tentative: \t48\n",
      "Rejected: \t35\n",
      "Iteration: \t60 / 100\n",
      "Confirmed: \t32\n",
      "Tentative: \t48\n",
      "Rejected: \t35\n",
      "Iteration: \t61 / 100\n",
      "Confirmed: \t32\n",
      "Tentative: \t48\n",
      "Rejected: \t35\n",
      "Iteration: \t62 / 100\n",
      "Confirmed: \t32\n",
      "Tentative: \t48\n",
      "Rejected: \t35\n",
      "Iteration: \t63 / 100\n",
      "Confirmed: \t32\n",
      "Tentative: \t48\n",
      "Rejected: \t35\n",
      "Iteration: \t64 / 100\n",
      "Confirmed: \t32\n",
      "Tentative: \t48\n",
      "Rejected: \t35\n",
      "Iteration: \t65 / 100\n",
      "Confirmed: \t33\n",
      "Tentative: \t47\n",
      "Rejected: \t35\n",
      "Iteration: \t66 / 100\n",
      "Confirmed: \t33\n",
      "Tentative: \t47\n",
      "Rejected: \t35\n",
      "Iteration: \t67 / 100\n",
      "Confirmed: \t37\n",
      "Tentative: \t43\n",
      "Rejected: \t35\n",
      "Iteration: \t68 / 100\n",
      "Confirmed: \t37\n",
      "Tentative: \t43\n",
      "Rejected: \t35\n",
      "Iteration: \t69 / 100\n",
      "Confirmed: \t37\n",
      "Tentative: \t43\n",
      "Rejected: \t35\n",
      "Iteration: \t70 / 100\n",
      "Confirmed: \t38\n",
      "Tentative: \t42\n",
      "Rejected: \t35\n",
      "Iteration: \t71 / 100\n",
      "Confirmed: \t38\n",
      "Tentative: \t42\n",
      "Rejected: \t35\n",
      "Iteration: \t72 / 100\n",
      "Confirmed: \t38\n",
      "Tentative: \t41\n",
      "Rejected: \t36\n",
      "Iteration: \t73 / 100\n",
      "Confirmed: \t38\n",
      "Tentative: \t41\n",
      "Rejected: \t36\n",
      "Iteration: \t74 / 100\n",
      "Confirmed: \t38\n",
      "Tentative: \t41\n",
      "Rejected: \t36\n",
      "Iteration: \t75 / 100\n",
      "Confirmed: \t39\n",
      "Tentative: \t40\n",
      "Rejected: \t36\n",
      "Iteration: \t76 / 100\n",
      "Confirmed: \t39\n",
      "Tentative: \t40\n",
      "Rejected: \t36\n",
      "Iteration: \t77 / 100\n",
      "Confirmed: \t39\n",
      "Tentative: \t39\n",
      "Rejected: \t37\n",
      "Iteration: \t78 / 100\n",
      "Confirmed: \t39\n",
      "Tentative: \t39\n",
      "Rejected: \t37\n",
      "Iteration: \t79 / 100\n",
      "Confirmed: \t39\n",
      "Tentative: \t39\n",
      "Rejected: \t37\n",
      "Iteration: \t80 / 100\n",
      "Confirmed: \t40\n",
      "Tentative: \t38\n",
      "Rejected: \t37\n",
      "Iteration: \t81 / 100\n",
      "Confirmed: \t40\n",
      "Tentative: \t38\n",
      "Rejected: \t37\n",
      "Iteration: \t82 / 100\n",
      "Confirmed: \t40\n",
      "Tentative: \t38\n",
      "Rejected: \t37\n",
      "Iteration: \t83 / 100\n",
      "Confirmed: \t41\n",
      "Tentative: \t37\n",
      "Rejected: \t37\n",
      "Iteration: \t84 / 100\n",
      "Confirmed: \t41\n",
      "Tentative: \t37\n",
      "Rejected: \t37\n",
      "Iteration: \t85 / 100\n",
      "Confirmed: \t41\n",
      "Tentative: \t37\n",
      "Rejected: \t37\n",
      "Iteration: \t86 / 100\n",
      "Confirmed: \t41\n",
      "Tentative: \t37\n",
      "Rejected: \t37\n",
      "Iteration: \t87 / 100\n",
      "Confirmed: \t41\n",
      "Tentative: \t37\n",
      "Rejected: \t37\n",
      "Iteration: \t88 / 100\n",
      "Confirmed: \t42\n",
      "Tentative: \t36\n",
      "Rejected: \t37\n",
      "Iteration: \t89 / 100\n",
      "Confirmed: \t42\n",
      "Tentative: \t36\n",
      "Rejected: \t37\n",
      "Iteration: \t90 / 100\n",
      "Confirmed: \t42\n",
      "Tentative: \t36\n",
      "Rejected: \t37\n",
      "Iteration: \t91 / 100\n",
      "Confirmed: \t42\n",
      "Tentative: \t36\n",
      "Rejected: \t37\n",
      "Iteration: \t92 / 100\n",
      "Confirmed: \t42\n",
      "Tentative: \t36\n",
      "Rejected: \t37\n",
      "Iteration: \t93 / 100\n",
      "Confirmed: \t42\n",
      "Tentative: \t36\n",
      "Rejected: \t37\n",
      "Iteration: \t94 / 100\n",
      "Confirmed: \t42\n",
      "Tentative: \t36\n",
      "Rejected: \t37\n",
      "Iteration: \t95 / 100\n",
      "Confirmed: \t43\n",
      "Tentative: \t34\n",
      "Rejected: \t38\n",
      "Iteration: \t96 / 100\n",
      "Confirmed: \t43\n",
      "Tentative: \t34\n",
      "Rejected: \t38\n",
      "Iteration: \t97 / 100\n",
      "Confirmed: \t43\n",
      "Tentative: \t34\n",
      "Rejected: \t38\n",
      "Iteration: \t98 / 100\n",
      "Confirmed: \t44\n",
      "Tentative: \t33\n",
      "Rejected: \t38\n",
      "Iteration: \t99 / 100\n",
      "Confirmed: \t44\n",
      "Tentative: \t33\n",
      "Rejected: \t38\n",
      "\n",
      "\n",
      "BorutaPy finished running.\n",
      "\n",
      "Iteration: \t100 / 100\n",
      "Confirmed: \t44\n",
      "Tentative: \t30\n",
      "Rejected: \t38\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "BorutaPy(estimator=RandomForestClassifier(n_estimators=124, n_jobs=-1,\n",
       "                                          random_state=RandomState(MT19937) at 0x27C1B4A4B40),\n",
       "         n_estimators='auto',\n",
       "         random_state=RandomState(MT19937) at 0x27C1B4A4B40, verbose=2)"
      ]
     },
     "execution_count": 121,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "feat_selector.fit(X, y)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 122,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([False, False,  True, False,  True,  True, False,  True, False,\n",
       "       False,  True,  True, False, False, False,  True,  True,  True,\n",
       "        True,  True, False,  True,  True,  True, False, False, False,\n",
       "        True, False, False, False,  True,  True,  True,  True, False,\n",
       "       False,  True,  True, False, False,  True,  True,  True,  True,\n",
       "        True, False, False,  True, False, False, False, False, False,\n",
       "       False, False, False, False, False, False, False, False, False,\n",
       "       False,  True, False,  True,  True, False, False, False, False,\n",
       "       False, False, False, False, False, False, False, False, False,\n",
       "        True,  True,  True, False, False, False, False, False,  True,\n",
       "       False,  True, False, False, False, False, False, False, False,\n",
       "        True, False,  True, False,  True, False, False,  True,  True,\n",
       "        True,  True, False, False,  True,  True, False])"
      ]
     },
     "execution_count": 122,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "feat_selector.support_"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 123,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "115"
      ]
     },
     "execution_count": 123,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.size(feat_selector.support_) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 124,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 2,  2,  1,  2,  1,  1,  2,  1,  2,  2,  1,  1,  2,  2,  2,  1,  1,\n",
       "        1,  1,  1,  2,  1,  1,  1,  2,  2, 10,  1,  2,  2, 14,  1,  1,  1,\n",
       "        1,  2,  2,  1,  1,  2,  8,  1,  1,  1,  1,  1, 32,  9,  1, 32, 32,\n",
       "       32, 20, 20, 32, 32, 11, 32, 32,  7,  2, 18, 32, 32,  1,  2,  1,  1,\n",
       "        2,  3,  3, 32, 12, 20, 32, 32, 32, 16, 32, 32, 32,  1,  1,  1,  5,\n",
       "       32, 32, 32, 14,  1,  2,  1,  5, 13, 17,  2,  2, 32,  2,  1,  2,  1,\n",
       "        2,  1,  2,  2,  1,  1,  1,  1,  2,  3,  1,  1,  2])"
      ]
     },
     "execution_count": 124,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "feat_selector.ranking_"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 125,
   "metadata": {},
   "outputs": [],
   "source": [
    "X_filtered = feat_selector.transform(X)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 134,
   "metadata": {},
   "outputs": [],
   "source": [
    "#X_filtered"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 131,
   "metadata": {},
   "outputs": [],
   "source": [
    "#np.size(X_filtered)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 132,
   "metadata": {},
   "outputs": [],
   "source": [
    "test=pd.DataFrame(X_filtered)\n",
    "#test"
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
   ]
  },
  {
   "cell_type": "code",
<<<<<<< HEAD
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_final = df_final.append(res)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(104, 116)"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_final.shape"
=======
   "execution_count": 133,
   "metadata": {},
   "outputs": [],
   "source": [
    "#test.shape"
>>>>>>> c00ff1eb1125945b6a827869699b3a8975564e60
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "X = df_final.iloc[:,0:115]\n",
    "y = df_final.iloc[:,-1]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.ensemble import RandomForestClassifier\n",
    "from sklearn.model_selection import GridSearchCV"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "rf_cs = RandomForestClassifier(n_jobs=-1,random_state=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "collapsed": true
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\ProgramData\\Anaconda3\\lib\\site-packages\\sklearn\\ensemble\\forest.py:245: FutureWarning: The default value of n_estimators will change from 10 in version 0.20 to 100 in 0.22.\n",
      "  \"10 in version 0.20 to 100 in 0.22.\", FutureWarning)\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "RandomForestClassifier(bootstrap=True, class_weight=None, criterion='gini',\n",
       "                       max_depth=None, max_features='auto', max_leaf_nodes=None,\n",
       "                       min_impurity_decrease=0.0, min_impurity_split=None,\n",
       "                       min_samples_leaf=1, min_samples_split=2,\n",
       "                       min_weight_fraction_leaf=0.0, n_estimators=10,\n",
       "                       n_jobs=None, oob_score=False, random_state=None,\n",
       "                       verbose=0, warm_start=False)"
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\ProgramData\\Anaconda3\\lib\\site-packages\\sklearn\\model_selection\\_split.py:657: Warning: The least populated class in y has only 4 members, which is too few. The minimum number of members in any class cannot be less than n_splits=10.\n",
      "  % (min_groups, self.n_splits)), Warning)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'n_estimators': 10}\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "RandomForestClassifier(bootstrap=True, class_weight=None, criterion='gini',\n",
       "                       max_depth=None, max_features='auto', max_leaf_nodes=None,\n",
       "                       min_impurity_decrease=0.0, min_impurity_split=None,\n",
       "                       min_samples_leaf=1, min_samples_split=2,\n",
       "                       min_weight_fraction_leaf=0.0, n_estimators=10, n_jobs=-1,\n",
       "                       oob_score=False, random_state=0, verbose=0,\n",
       "                       warm_start=False)"
      ]
     },
     "execution_count": 35,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "param_grid = { 'n_estimators': [10,50,100,200,250,500,700],}\n",
    "    \n",
    "CV_rfc = GridSearchCV(estimator=rf_cs, param_grid=param_grid, cv= 10)\n",
    "CV_rfc.fit(X, y)\n",
    "print (CV_rfc.best_params_)\n",
    "op_RF = CV_rfc.best_estimator_\n",
    "op_RF.fit(X,y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "ename": "AttributeError",
     "evalue": "'RandomForestClassifier' object has no attribute 'test_score_'",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mAttributeError\u001b[0m                            Traceback (most recent call last)",
      "\u001b[1;32m<ipython-input-39-d24c41c48e53>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0mop_RF\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtest_score_\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[1;31mAttributeError\u001b[0m: 'RandomForestClassifier' object has no attribute 'test_score_'"
     ]
    }
   ],
   "source": []
  },
  {
   "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.8.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}